La segmentación de una gran trama

Question

En realidad estoy tratando de convertir un mapa de bits en un SpatialPolygons. Tengo una gran trama (12 millones de píxeles) con un montón de ceros y unos valores negativos (que representan a la depresión en un DEM). Me gustaría clúster de todos estos valores negativos en los polígonos, para obtener un mapa de todas las depresiones en mi mapa. Soy un R principiante y no puedo encontrar una función adecuada para hacer el trabajo.

¿Alguien tiene una idea ?

Gracias por tu ayuda.

Answer 1

La herramienta (parte de la trama & rgeos bibliotecas) que usted necesita para tener en cuenta es:

output_shp <- rasterToPolygons(x, fun=function(x){x<0}, n=4, na.rm=TRUE, digits=12, dissolve=TRUE)

Donde 'x' es su DEM-raster.
Esto se traduce en un 'SpatialPolygonsDataFrame', que es lo que buscábamos.

Sin embargo, antes de ejecutar el rasterToPolygons-herramienta, sería aconsejable para reclasificar su DEM en -1s y 0s:

reclassified_raster <- reclassify(x,c(-Inf,0,-1, 0,Inf,1))

Esto le dará una trama con todos los valores negativos como '-1' y todos los valores positivos, '0', para su uso en la rasterToPolygons-herramienta. Si este paso no se toma, entonces cada una de las alturas en el modelo de elevación se convierta en su propio polígono.