Menghasilkan shapefile grid di R?

Saya memiliki shapefile berikut dari garis lintang / bujur 10x10 derajat yang saya buat di QGIS. Saya bisa membacanya ke R menggunakan paket rgdal.

Grid<-readOGR(".","GridShapeFile")

Ini memiliki atribut dan struktur berikut.

summary(Grid)

Object of class SpatialPolygonsDataFrame
Coordinates:
   min max
x -180 190
y -100  90
Is projected: FALSE 
proj4string :
[+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0]
Data attributes:
       ID             XMIN           XMAX           YMIN           YMAX    
 Min.   :  0.0   Min.   :-180   Min.   :-170   Min.   :-100   Min.   :-90  
 1st Qu.:175.5   1st Qu.: -90   1st Qu.: -80   1st Qu.: -60   1st Qu.:-50  
 Median :351.0   Median :   0   Median :  10   Median : -10   Median :  0  
 Mean   :351.0   Mean   :   0   Mean   :  10   Mean   : -10   Mean   :  0  
 3rd Qu.:526.5   3rd Qu.:  90   3rd Qu.: 100   3rd Qu.:  40   3rd Qu.: 50  
 Max.   :702.0   Max.   : 180   Max.   : 190   Max.   :  80   Max.   : 90 

# An example row of the data
Grid[50,]
class       : SpatialPolygonsDataFrame 
features    : 1 
extent      : -60, -50, 70, 80  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84+towgs84=0,0,0 
variables   : 5
names       : ID, XMIN, XMAX, YMIN, YMAX 
min values  : 49,  -60,  -50,   70,   80 
max values  : 49,  -60,  -50,   70,   80 

Saya ingin belajar cara membuat kisi / shapefile yang sama ini sepenuhnya di dalam R, daripada di QGIS, serta kisi dengan ukuran lain (1x1,5x5, dll.). Tujuan utamanya adalah untuk menggunakan over () dalam paket sp untuk overlay grid ini ke file poligon lain dan menghitung jumlah grid yang berpotongan dengan setiap poligon. Saya sudah tahu bagaimana melakukan langkah selanjutnya ini (saya pikir), saya hanya ingin tahu cara membuat layer grid di R.

Lihatlah fungsi raster dalam paket raster. Ini akan memungkinkan Anda membuat raster dengan batas tertentu, jumlah baris / kolom dan resolusi.

Di sini saya akan menggunakan karakteristik ringkasan data Anda untuk membuat raster 100x100 dalam batas yang ditentukan. Saya melewati objek sejauh untuk mendefinisikan batas x dan y. Anda juga dapat menggunakan argumen spesifik (xmn, xmx, ymn, ymx) di dalam fungsi raster.

library(raster) 
library(sp) 

r <- raster(extent(matrix( c(-180, -100, 190,  90), nrow=2)), nrow=100, ncol=100, 
            crs = "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0")            
  r[] <- 1:ncell(r)
  summary(r)
  print(r)
  plot(r)

Sangat mudah untuk memaksa objek raster ke objek sp grid menggunakan;

sp.r <- as(r, "SpatialPixelsDataFrame")
  class(sp.r)   
  spplot(sp.r, "layer")