R 또는 Python을 사용하여 큰 래스터를 폴리 라인으로 변환하는 가장 빠른 방법은 무엇입니까?

전역 수역 (1 비트 값 0 및 1)이있는 큰 래스터 파일 (129600 x 64800 픽셀)을 가지고 바다와 내륙 수 해안선을 추출하려고합니다.

래스터에서 폴리 라인으로 변환하기 위해 ArcGIS 및 QGIS를 사용해 보았지만 시간이 오래 걸립니다.

누구 든지이 작업에 더 나은 / 빠른 방법 (Python 또는 R) 또는 더 나은 도구를 알고 있습니까?

최신 정보

• R : rasterToContour는 빠르고 정확할 수 있지만 광산과 같은 매우 큰 데이터 세트 (8,398,080,000 픽셀)가있는 경우 매우 많은 양의 RAM (16GB 이상)이 필요하거나 R이 하드 드라이브에서 더 많은 처리를 수행하도록 강요합니다. 또한 나이가 걸릴 것입니다.
• Python / GDAL : gdal_poligonize는 폴리 라인 대신 다각형을 만듭니다.

업데이트 2

• R rasterToContour : rasterToContour는 원하는 결과를 제공하지 않습니다. ArcGIS (래스터-다각형, 피처-라인)와 비교하면 아래 예제와 같이 정확한 픽셀 외곽선을 추출하지 않습니다.

rasterToContour 결과

ArcGIS 결과

업데이트 3

파이썬 / GDAL : 테스트 데이터 세트에서 ArcGIS에 대해 명령 줄에서 gdal_polygonize를 실행했으며 결과가 매우 명확했습니다.

• 페달 : 49 초
• ArcGIS : 1.84 초

나는 R과 함께 일하고 과거 rasterToPolygons에는 raster패키지에서 사용 했지만 이제는 gdal_polygonizeRJohn Baumgartner가 선호 합니다. 그것은 기초 gdal_polygonize.py이며 훨씬 빠릅니다. John Baumgartner는이 코드를 게시하고 블로그 에서 사용하는 예제를 제공했습니다 .

파이썬에 익숙하다면 gdal_polygonize.py물론 직접 사용할 수 있습니다 .

후손을 위해, 나는 이러한 유형의 작업을 신속하게 수행하기 위해 stars::패키지로 성공했습니다 R.

library(raster)
library(stars)
library(sf)
library(magrittr)

f <- system.file("external/test.grd", package="raster")
r <- raster(f)
r[r[] < 750] <- 0
r[r[] >= 750] <- 1

x <- st_as_stars(r) %>% 
  st_as_sf() %>% # this is the raster to polygons part
  st_cast("MULTILINESTRING") # cast the polygons to polylines

plot(x)

plot(r)
plot(x, add = TRUE)

시도 rasterToContour로부터 래스터 패키지로 제공된다.

f <- system.file("external/test.grd", package="raster")
r <- raster(f)
r[r[] < 750] <- 0
r[r[] >= 750] <- 1

x <- rasterToContour(r)
class(x)
> [1] "SpatialLinesDataFrame"
> attr(,"package")
> [1] "sp"

plot(r)
plot(x, add=TRUE)

그런 다음 아래 코드를 사용하여 파일을 로컬 폴더 (예 : 'ESRI Shapefile'(. shp))에 쉽게 쓸 수 있습니다. 한 번 봐 가지고 ogrDrivers에서 rgdal는 시스템과 호환되는 드라이버를 발견 할 수 있습니다.

library(rgdal)
writeOGR(x, dsn = getwd(), layer = "coastlines", driver = "ESRI Shapefile")

내가 GDAL을 좋아하지만 다각형 화 도구는 응용 프로그램에 너무 느 렸습니다.

빠른 대안이다 gdal_trace_outline에서 하건데 GDAL 스크립트 도 더 많은 옵션에 대한 내성을 가지고, 도넛 등

마찬가지로 gdal_polygonize이 또한 당신이 나중에 변환해야 것 다각형을 생성합니다 ogr2ogr -nlt MULTILINESTRING.

단점은 Linux 또는 Mac OsX 시스템을 사용하지 않는 한 직접 컴파일해야한다는 것입니다.