2016-05-23 2 views
2

Rでラスタを処理する際に問題があります。私はそれらをデータを読み取り、プロットすることができます:私はポリゴンに私のラスタを変換しようとしたときR:RasterToPolygonの動作が遅く、動作を停止します。

enter image description here

ただし、Rは

enter image description here

動作を停止し、私はパッケージを更新しようとしました、Windows 7上でOSX上にRを再インストールすると、常に同じ結果が得られます。

さらに、私は投影法をNAに設定しようとしましたが、これは役に立ちませんでした。このような問題を経験したことがありますか?あなたは問題を解決する方法をいくつか提案していますか?私は非常に任意の提案に感謝!

MYデータはここにある: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) 
+0

gdal_polygonizeを使用することをお勧めします。http://www.gdal.org/gdal_polygonize.html R rasterToPolygons関数は非常に遅く、しばしば正しく動作しません大きなデータセットについてgdal機能に関する助けが必要な場合は、教えてください。 – maRtin

答えて

0

これはスピードの問題を解決していませんが、少なくとも私のプロセスが実行され、Rは、クラッシュが...私はrasterToPolygon機能のためにすべての引数を定義する必要はありません。

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

たぶん問題が原因でした

2

私はまだ50担当者を持っていないので、私はコメントを追加することはできませんが、このweb-pageを見て。それは私がいつもメガ・スローであることがわかったrasterToPolygons関数を扱っています。著者(John Baumgartner)はプロセスをスピードアップする機能を提供しています。彼のテストによると、時間の約1/7000を要します。

+0

私はそれを試して、素晴らしい作品! ;)ありがとう – maycca

関連する問題