我有一个在国家公园的qgis中制作的shapefile,该文件具有土地覆盖类型。我剪裁了泰国的一个大shapefile和一个只是公园的较小shapefile。该文件(DPKY.lc5)现在为10.5 Mb。当我运行代码时,它要花很多时间,而且还没有成功。为什么会这样?
library(raster)
DPKY.lc5 <- shapefile("dpky.lc5.shp")
DPKY.lc5 <- spTransform(DPKY.lc5, CRS('+init=EPSG:4326'))
DPKY.lc5$VALUE<-as.numeric(DPKY.lc5$VALUE)
rr <- raster(DPKY.lc5, res=0.01)
rr1 <- rasterize(DPKY.lc5, rr, field="VALUE")
rr1
我希望它能正常工作,因为它只有10.5Mb,但需要花费很长时间。我需要这个栅格才能工作,因此我可以使用另一个数据框中的数据点来显示公园中大象栖息地的使用频率。我成功使用了其他栅格文件来获取高程和坡度纵横比,但这次却不喜欢它。