R: RasterToPolygon arbeitet langsam, dann aufhören zu arbeiten

Question

Ich habe ein Problem, Raster in R zu verarbeiten, auch wenn ich es vorher getan habe. Ich kann die Daten lesen und zeichnen sie:

aber als ich versuchte, mein Raster zu Polygon zu konvertieren, stoppt der R

Ich habe versucht, Arbeitspakete zu aktualisieren , R, unter OSX, unter Windows 7, immer mit dem gleichen Ergebnis neu installieren.

Außerdem habe ich versucht, meine Projektion auf NA zu setzen, aber das half auch nicht. Hast du jemals ein Problem wie dieses erlebt? Haben Sie Vorschläge, wie Sie das Problem lösen können? Ich schätze alle Vorschläge sehr!

Meine Daten sind hier: https://ulozto.sk/x4pHuyra/data-to-stack-zip

Und mein Code:

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

Dies löst nicht das Geschwindigkeitsproblem, aber zumindest meine Prozesse laufen und R nicht abstürzt ... ich alle Argumente für rasterToPolygon Funktion definieren müssen:

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

Vielleicht ist das Problem aufgrund wurde zu vielen NA-Werten im ursprünglichen Raster.

Answer 2

I 50 rep noch nicht haben, so kann ich einen Kommentar nicht hinzufügen, haben aber einen Blick auf diese web-page . Es handelt sich um die rasterToPolygons Funktion, die ich immer als mega-langsam empfunden habe, und der Autor (John Baumgartner) bietet eine Funktion von ihm, um den Prozess zu beschleunigen. Es dauert ungefähr 1/7000 der Zeit, nach seinen Tests.