6 votos

El cálculo de Landsat 7 NDVI de las imágenes utilizando QGIS

He Landsat 7 imágenes desde el 2011 hasta el 2016 para mis estudios agrícolas. Primero he corregido el slc-off problema de la banda 3 y la banda 4. A continuación, he calculado el NDVI utilizando tanto la trama de la calculadora y el índice de vegetación (pendiente) basado en el procesamiento de la caja de herramientas. Me di cuenta de que el rango de datos del NDVI Imágenes son diferentes. La trama resultados de la calculadora de entre -0.242 - 0.568 y el índice de vegetación bajo la caja de herramientas de procesamiento de resultados entre -0.225993 - 0.553768. Este puede ser debido al redondeo o hay algo que me puede estar faltando?

Yo también uso la capa de trama herramienta de estadísticas para ver si estas dos imágenes NDVI son muy diferentes. Los resultados son:

A partir de índices de vegetación algoritmo: Válido células: 42384039 No-celdas de datos: 15942162 Valor mínimo: -0.991935491562 Valor máximo: 0.982142865658 Suma: 2857838.41327 Valor medio: 0.0674272315876 Desviación estándar: 0.184433938177

De la trama de la calculadora (que yo elija de la lista desplegable de la barra de menú de la Trama) Válido células: 42696562 No-celdas de datos: 15629639 Valor mínimo: -1.0 Valor máximo: 1.0 Suma: 2862525.41327 Valor medio: 0.067043463904 Desviación estándar: 0.20274331204

This is my QGIS screen:

6voto

Boushley Puntos 1360

Comprobación y comparación de una imagen Landsat 7 (área diferente a la tuya, pero con el mismo sensor de datos) revela que los valores son, de hecho, tanto "correcta".

Usted (y yo) hemos caído en la trampa de la presentación de Trama definición de la eficacia de la Trama estilo. Este se encuentra por right clicking the layer -> properties -> Style

Raster style screen

En la pantalla puede definir cómo cargar la min/max valores de datos. Como norma, estos parecen ser definido por el número Acumulado de corte de 2.0 98.0%. Experimente con los ajustes y te darás cuenta de los valores en el Min Max cajas debajo de la Gradiente de Color Ajustes de cambio.

Al tener los datos de prueba a través de ambos métodos descritos puedo obtener los siguientes resultados cuando se consultan las capas (por ejemplo punto):

Data values at sample point

La parte superior a la parte inferior representan estos valores:

  • 0.4603175 - el valor calculado de La Raster -> Raster calc herramienta
  • 0.4603175 - el valor devuelto por la SAGA -> Vegetation Index (Slope based) NDVI
  • 34 - la Banda 3 (Rojo) de Landsat 7 ETM+
  • 92 - la Banda 4 (NIR) de Landsat 7 ETM+

Manual rápido de control de cálculo revela que (B4-B3)/(B4+B3) = 0.4603175 (redondeado).

Por último, podemos comparar los datos máximos valores de calculados los rásteres en la capa de Metadatos. Aquí podemos ver una discrepancia entre las capas, aunque. La primera imagen muestra la Trama calc valores, la segunda muestra el índice de Vegetación de valores.

Raster calc metadata

NDVI metadata

EDITAR:

Resulta que la discrepancia en los valores de los tallos de los problemas con el NoData y los bordes de la imagen. En mi configuración de prueba de los últimos 3 filas de píxeles, de hecho, fueron diferentes entre los dos rásteres. El recorte a la zona de interés eliminan la distorsión en los valores y en ambos se muestra idéntico Min/Max en los metadatos.

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X