QGISで作成した次の10x10度の緯度/経度グリッドのシェープファイルがあります。rgdalパッケージを使用してRに読み込むことができます。
Grid<-readOGR(".","GridShapeFile")
次の属性と構造があります。
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
QGISではなく、他のサイズのグリッド(1x1、5x5など)ではなく、完全にR内でこの同じグリッド/シェープファイルを作成する方法を学びたいと思います。最終的な目標は、spパッケージでover()を使用して、このグリッドをポリゴンの別のファイルにオーバーレイし、各ポリゴンが交差するグリッドの数を数えることです。私はすでにこの次のステップを実行する方法を知っています(私は思う)、Rでグリッドレイヤーを生成する方法を知りたいだけです。