¿Polygonize área de vegetación de Landsat7 RGB utilizando QGis / Bilko?

Question

Estoy tratando de extraer las áreas con vegetación, como un polígono en forma de un landsat7-imagen, a 30m de resolución, RGB, utilizando QGis 1.8 o Bilko o (esperemos que en el futuro cercano) de ArcGIS. El área se encuentra en el Norte de Perú, temporada seca, por lo que, básicamente quiero extraer todos los matices de verde.

Nunca he trabajado realmente con datos ráster y no sé qué hacer. He encontrado que "hay índices de vegetación", que puede ser calculado, y han encontrado Bilko como un software que se podía hacer esto con. Sin embargo, si lo tengo necesito más bandas de RGB a ello, ¿es esto cierto? No estoy seguro.

En este sitio web también he encontrado hilos sobre cómo calcular el área de vegetación en QGis - uso de GDAL o herramientas de SEXTANTE: calcular el área de las vegetaciones en un archivo raster

Sin embargo, me gustaría no sólo calcular pero extraer los datos como un polígono.

¿Cuál sería la mejor manera de hacer esto? Nunca he trabajado con Bilko, GDAL o SEXTANTE, por supuesto voy a entrar en eso y tratar de aprender, pero me gustaría saber lo que va a trabajar antes de que me pase horas y horas sólo para descubrir que no puedo hacerlo...

P. S.: yo probablemente han de ARCGis en un tiempo, así que si es mucho más fácil con ARCGis podía esperar para que...

Answer 1

También hay otra posibilidad para resolver el problema. El software "de Monteverdi", que se basa en la OTB de la Biblioteca, podría manejar eso. Es una fuente abierta de la detección remota de software, desarrollado por el CNES (agencia espacial francés) para procesar sus Pléyades-Imágenes.

Sugiero la Versión 1.18.

http://sourceforge.net/projects/orfeo-toolbox/files/Monteverdi/Monteverdi-1.18/

Hay un filtrado del módulo incluido, el cual es llamado "Componente Conectado Segmentación". Básicamente hace lo siguiente (más detalles se pueden encontrar aquí: http://wiki.orfeo-toolbox.org/index.php/Connected_component_segmentation_module):

1. Enmascaramiento de la imagen de acuerdo a sus reglas (en este caso, debería ocultar todo con un NDVI por debajo de 0.2 por ejemplo. El NDVI es una función integrada)
2. La segmentación de los píxeles restantes (creación de polígonos)
3. Si es necesario descartar pequeñas áreas (conjunto de un mínimo tamaño)
4. Exportar el resultado como un *.shp

He utilizado este exactamente para ese propósito. No en los datos de Landsat, pero las imágenes aéreas. Una NIR canal es obligatorio, aunque.

Answer 2

Si no tienes experiencia en cualquiera de las herramientas mencionaste que le sugiero use R, tiene una muy interfaz de usuario para las herramientas gdal. Podría hacer lo que quiera de esa manera:

library(raster)
b.red = 1 # according to your input image
b.nir = 2 # according to your input image
threshold = 0.2 # for vegetation detection
inDs <- 'your_raster.tif'
inRst.red <- raster(inDs, b.red)
inRst.nir <- raster(inDs, b.nir)
NDVI <- (nir-red)/(nir+red)
NDVI[NDVI>=threshold] <- 1
NDVI[NDVI<threshold] <- 0
setNA(NDVI) <- 0
p <- rasterToPolygons(NDVI, dissolve=TRUE)

Answer 3

Muchas gracias por las respuestas, pregunté la pregunta hace un tiempo y todo resultó diferente. He intentado hacer esto para hacer una pasantía en Perú sin tener acceso a una buena computadora - no funcionaba. Así lo hice en Alemania con otros programas disponibles.

Finalmente usé una imagen ASTER (15m de resolución, verde / rojo / banda de IR) y calcula un NDVI, que luego clasificado, usando los programas Erdas y ArcMap.