Suite

R : convertir l'objet d'étendue en fichier de formes

R : convertir l'objet d'étendue en fichier de formes


Je veux convertir R 'library:raster"Le degré' objet à un fichier de formes. J'ai le CRS disponible aussi.

e= c(4304916, 4305325, 365216, 365439)#xmin, xmax, ymin, ymax proj4string(e) = CRS("+proj=utm +zone=18 +ellps=GRS80 +datum=NAD83 +units=m +no_defs ")

Dans l'ancienne version des données raster, il y avait une fonction spécifique "" 'extentTopolygon' qui est maintenant dépréciée. Mon approche actuelle consiste à créer un objet raster, à modifier sa projection, à affecter l'étendue actuelle et à le convertir en fichier de formes, ce qui n'est pas très efficace.


Vous pouvez effectuer les opérations suivantes :

library(raster) e <- extend( c(4304916, 4305325, 365216, 365439) ) p <- as(e, 'SpatialPolygons') crs(p) <- "+proj=utm +zone=18 +ellps=GRS80 +datum=NAD83 +units=m +no_defs" shapefile(p, 'file.shp')