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 值过多导致的问题..