R:RasterToPolygon 工作缓慢,然后停止工作
R: RasterToPolygon works slow, then stop working
我在 R 中处理栅格时遇到问题,即使我以前做过。我可以读取数据并绘制它们:
但是,当我尝试将光栅转换为多边形时,R 停止工作
我尝试更新软件包,在 OSX 上重新安装 R,在 Windows 7 上,总是得到相同的结果。
此外,我尝试将我的投影设置为 NA,但这都无济于事。你遇到过这样的问题吗?你对如何解决这个问题有什么建议吗?我非常感谢任何建议!
我的数据在这里:https://ulozto.sk/x4pHuyra/data-to-stack-zip
还有我的代码:
library(sp)
library(raster)
library(rgeos)
library(spatstat)
library(rgdal)
library(maptools)
require(spdep)
# read final GFW raster - in JTSK
gfw_13<-raster("H:/.../gfw_ext_forest_03.tif")
projection(gfw_13)<-NA
# read shp
manag<-readOGR(dsn = "H:/...",
layer = "cleared_management_by_NP")
projection(manag)<-NA
# convert raster GFW to shp # !!!!! here it stops to work
pol.gfw<-rasterToPolygons(gfw_13, dissolve = T)
我还没有 50 个代表,所以我不能添加评论,但请看看这个 web-page。它处理 rasterToPolygons
函数,我一直发现它非常慢,作者 (John Baumgartner) 提供了他的一个函数来加速这个过程。根据他的测试,这大约需要 1/7000 的时间。
这不能解决速度问题,但至少我的进程 运行 和 R 不会崩溃...我需要为 rasterToPolygon
函数定义所有参数:
pol.gfw <- rasterToPolygons(gfw_13, fun=NULL, n=4, na.rm=TRUE, digits=12, dissolve=TRUE)
可能是由于原始栅格中的 NA 值过多导致的问题..
我在 R 中处理栅格时遇到问题,即使我以前做过。我可以读取数据并绘制它们:
但是,当我尝试将光栅转换为多边形时,R 停止工作
我尝试更新软件包,在 OSX 上重新安装 R,在 Windows 7 上,总是得到相同的结果。
此外,我尝试将我的投影设置为 NA,但这都无济于事。你遇到过这样的问题吗?你对如何解决这个问题有什么建议吗?我非常感谢任何建议!
我的数据在这里:https://ulozto.sk/x4pHuyra/data-to-stack-zip
还有我的代码:
library(sp)
library(raster)
library(rgeos)
library(spatstat)
library(rgdal)
library(maptools)
require(spdep)
# read final GFW raster - in JTSK
gfw_13<-raster("H:/.../gfw_ext_forest_03.tif")
projection(gfw_13)<-NA
# read shp
manag<-readOGR(dsn = "H:/...",
layer = "cleared_management_by_NP")
projection(manag)<-NA
# convert raster GFW to shp # !!!!! here it stops to work
pol.gfw<-rasterToPolygons(gfw_13, dissolve = T)
我还没有 50 个代表,所以我不能添加评论,但请看看这个 web-page。它处理 rasterToPolygons
函数,我一直发现它非常慢,作者 (John Baumgartner) 提供了他的一个函数来加速这个过程。根据他的测试,这大约需要 1/7000 的时间。
这不能解决速度问题,但至少我的进程 运行 和 R 不会崩溃...我需要为 rasterToPolygon
函数定义所有参数:
pol.gfw <- rasterToPolygons(gfw_13, fun=NULL, n=4, na.rm=TRUE, digits=12, dissolve=TRUE)
可能是由于原始栅格中的 NA 值过多导致的问题..