J'ai le fichier de formes suivant d'une grille de latitude / longitude de 10 x 10 degrés que j'ai créée dans QGIS. Je peux le lire dans R en utilisant le paquetage rgdal.
Grid<-readOGR(".","GridShapeFile")
Il a les attributs et la structure suivants.
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
Je voudrais apprendre à créer cette même grille / fichier de formes entièrement dans R, plutôt que dans QGIS, ainsi que des grilles d'autres tailles (1x1,5x5, etc.). Le but ultime est d'utiliser over () dans le package sp pour superposer cette grille sur un autre fichier de polygones et compter le nombre de grilles coupées par chaque polygone. Je sais déjà comment faire cette prochaine étape (je pense), je veux juste savoir comment générer la couche de grille en R.