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

Question

我在 R 中处理栅格时遇到问题，即使我以前已经这样做过。我可以读取数据并绘制它们:

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

我尝试过在 OSX、Windows 7 上更新软件包、重新安装 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 的时间。