R：RasterToPolygon工作缓慢，然后停止工作

Question

我在R中处理栅格有问题，即使我之前已经完成了。我可以读取数据并绘制它们：

但是，当我尝试将光栅转换为多边形时，R停止工作

我试图在Windows 7上的OSX上更新软件包，重新安装R，总是得到相同的结果。

此外，我试图将我的预测设定为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)

Answer 1

我还没有50个代表，所以我无法添加评论，但请查看此web-page。它涉及rasterToPolygons函数，我总是发现它很慢，而作者（John Baumgartner）提供了他的功能来加速这个过程。根据他的测试，它大约需要1/7000的时间。

Answer 2

这不能解决速度问题，但至少我的进程运行且R没有崩溃......我需要为rasterToPolygon函数定义所有参数：

pol.gfw <- rasterToPolygons(gfw_13, fun=NULL, n=4, na.rm=TRUE, digits=12, dissolve=TRUE) 

可能问题是由于原始栅格中存在大量的NA值。