将shapefile转换为光栅.tif文件R
我添加了一个新变量来对形状文件的单元格进行编码,我想将该文件另存为将shapefile转换为光栅.tif文件R,r,geospatial,raster,spatial,sf,R,Geospatial,Raster,Spatial,Sf,我添加了一个新变量来对形状文件的单元格进行编码,我想将该文件另存为.tif,并获得新变量的所有信息。我将它保存为一个新的.shp文件,然后使用shp2master函数,但我得到了一个错误 beijing10<-sf::read_sf("Landuse/Beijing2010.shp") beijing10<-mutate(beijing10,land=ifelse(CODE%in%c(41,42),"Cultivated Land","
.tif
,并获得新变量的所有信息。我将它保存为一个新的.shp
文件,然后使用shp2master
函数,但我得到了一个错误
beijing10<-sf::read_sf("Landuse/Beijing2010.shp")
beijing10<-mutate(beijing10,land=ifelse(CODE%in%c(41,42),"Cultivated Land","Unused Land"))
st_write(beijing10,"beijing2010new.shp")
shp2raster("beijing2010new.shp",column = "land")
是否有其他选项可以获取
.tif
文件。有人能帮助我理解转换是如何工作的,我应该如何添加正确的代码,这样我就不需要使用ArcMap来获取文件了。谢谢你当问一个R问题时,请提供一个最小的、可复制的、自包含的示例,如下所示
library(raster)
p <- shapefile(system.file("external/lux.shp", package="raster"))
r <- raster(p, res=0.01)
p$land <- ifelse(p$ID_2 > 6, "Cultivated Land", "Unused Land")
库(光栅)
P
library(raster)
p <- shapefile(system.file("external/lux.shp", package="raster"))
r <- raster(p, res=0.01)
p$land <- ifelse(p$ID_2 > 6, "Cultivated Land", "Unused Land")
p$land <- as.factor(p$land)
r <- rasterize(p, r, "land")
r
#class : RasterLayer
#dimensions : 73, 78, 5694 (nrow, ncol, ncell)
#resolution : 0.01, 0.01 (x, y)
#extent : 5.74414, 6.52414, 49.45162, 50.18162 (xmin, xmax, ymin, ymax)
#crs : +proj=longlat +datum=WGS84 +no_defs
#source : memory
#names : layer
#values : 1, 2 (min, max)