Suite

Définir la taille des cellules à partir de spsample en mètres

Définir la taille des cellules à partir de spsample en mètres


J'utilise un shapefile pour créer une grille régulière en utilisantéchantillon, desppaquet.

Deux problèmes ici: (1) Je voudrais définir la taille de mes cellules de grille en mètres, mais j'ai vérifié la documentation despet il n'est pas clair quelle est l'unité de mesure par défaut ; (2) Ce code fonctionne pour créerhexagonalgrille mais je ne peux pas faire en sorte que cela fonctionne pour créer unrégulierla grille. Des idées?

Voici un exemple reproductible :

# charger les bibliothèques bibliothèque(sp) bibliothèque(rgdal)

Télécharger un fichier de formes

URL <- "ftp://dnrftp.dnr.ne.gov/pub/data/state/Legislative2010UTM.zip" fil <- basename(URL) if (!file.exists(fil)) download.file(URL, fil ) # décompressez et récupérez la liste des fichiers fils <- unzip(fil) # trouvez le fichier de formes qu'il contient shp <- grep("shp$", fils, value=TRUE) # récupérez le premier calque à partir de celui-ci lay <- ogrListLayers(shp )[1] # lu dans le shapefile leg <- readOGR(shp, lay) plot(leg)

Définir la taille des bacs hexagonaux

HexPts <-spsample(leg,type="hexagonal", offset=c(0,0), cellsize=100000) # Créer une grille hexagonale sous forme de polygones spatiaux HexPols <- HexPoints2SpatialPolygons(HexPts) plot(HexPols, add=T)


Utilisationéchantillonavectaille de la cellule:

pts = spsample(leg, cellsize=c(10000,10000), type="regular") plot(leg) points(pts)

vérifier l'espacement des coordonnées :

> coordonnées(pts)[1:10,] x1 x2 17 244356 4437161 18 254356 4437161 19 264356 4437161

vous voyez la coordonnée x passer par 10 000 (et si vous regardez plus tard, vous verrez la coordonnée y aussi)


Voir la vidéo: Gérer le format de cellule sous Excel