Detección de cambios de valor de píxel

La detección de cambios mediante datos continuos se basa en el cálculo de diferencias entre imágenes ópticas o datos modelados, como la temperatura o la velocidad del viento, píxel por píxel.

El propósito de comparar datos modelados continuos es normalmente el de identificar las áreas que han cambiado en magnitud o en una dirección determinada, normalmente a lo largo de un periodo de tiempo. ArcGIS AllSource proporciona varios métodos para comparar varios datasets ráster continuos.

Detección de cambio de valor de píxel entre dos rásteres de índice de vegetación

La función ráster Calcular cambio calcula la diferencia entre dos capas ráster sobre la marcha. Solo se procesan los píxeles dentro de la extensión visible y, si se aleja más allá de la resolución de origen de los datos de entrada, los píxeles se remuestrean a un tamaño mayor para acelerar el procesamiento. Esta es una forma de calcular el cambio para una rápida visualización y evaluación del cambio. Puede combinar esta función con otras funciones ráster para un flujo de trabajo más complejo.

La herramienta de geoprocesamiento Calcular ráster de cambio calcula la diferencia entre dos datasets ráster y genera un nuevo dataset ráster que contiene la información del cambio. Puede combinar esta herramienta con otras herramientas para un flujo de trabajo más complejo.

El Asistente de detección de cambios combina herramientas y funciones para guiarle a través del proceso de realización de detección de cambios categóricos, detección de cambios de valor de píxel, detección de cambios espectrales o detección de cambios de serie temporal. Puede comparar dos datasets ráster o dos divisiones en zonas en un dataset ráster multidimensional.

En la siguiente sección se proporcionan detalles sobre cada panel del Asistente para detección de cambio al realizar una detección de cambio continua.

Asistente de detección de cambios

El Asistente de detección de cambios se inicia desde el botón desplegable de detección de cambios de la pestaña Imágenes, en el grupo Análisis. El botón no está disponible si no está trabajando en una escena de mapa 2D o si no tiene la extensión de Image Analyst .

Panel Configurar

El primer panel del Asistente para detección de cambio es el panel Configurar, donde puede establecer el valor de Método de detección de cambio. Para comparar datasets continuos, establezca el valor de Método de detección de cambios en Cambio del valor del píxel.

Si se selecciona un ráster continuo en el panel Contenido al iniciar el asistente, la opción Método de detección de cambios se establece en Cambio del valor de píxel de forma predeterminada.

ParámetroDescripción

De ráster

El primer ráster que se utilizará en el cómputo. Para evaluar el cambio de la hora 1 (anterior) a la hora 2 (posterior), introduzca el ráster de la hora 1.

A ráster

El segundo ráster que se utilizará en el cómputo. Para evaluar el cambio de la hora 1 (anterior) a la hora 2 (posterior), introduzca el ráster de la hora 2.

Para comparar dos divisiones en zonas en un dataset multidimensional, introduzca el ráster multidimensional como el parámetro De ráster. El panel Configurar se actualiza con los parámetros de la tabla siguiente.

ParámetroDescripción

Variable

La variable que se va a analizar.

Dimensión

El campo de dimensión que se va a utilizar para comparar las divisiones en zonas.

Desde división en zonas

La primera división en zonas que se va a utilizar en el cálculo. Para evaluar el cambio de la hora 1 (anterior) a la hora 2 (posterior), introduzca la división en zonas de la hora 1.

A división en zonas

La segunda división en zonas que se va a utilizar en el cálculo. Para evaluar el cambio de la hora 1 (anterior) a la hora 2 (posterior), introduzca la división en zonas de la hora 2.

Extensión de procesamiento

La extensión de procesamiento para la salida.

Factor de omisión X

El factor de omisión X utilizado para calcular estadísticas y el histograma para el ráster de cambio.

Factor de omisión Y

El factor de omisión Y utilizado para calcular estadísticas y el histograma para el ráster de cambio.

Panel Diferencia de banda

El panel Diferencia de banda le permite especificar cómo comparar los datasets ráster continuos de entrada. Si los datos son imágenes multibanda, puede comparar una banda única entre las dos imágenes o calcular primero un índice de banda antes de comparar.

ParámetroDescripción

Tipos de diferencia

El método de diferencia a utilizar en el cálculo.

  • Absoluta: diferencia matemática, o resta, entre los valores de píxel del parámetro De ráster y los valores de píxel del parámetro A ráster.
    Salida = (A ráster) - (De ráster)
    Esta es la opción predeterminada.
  • Relativa: diferencia en los valores de píxel teniendo en cuenta las cantidades de los valores que se comparan.
    Salida = (A ráster - De ráster) / máx(A ráster, De ráster)

Método de diferencia de banda

Elija cómo desea comparar los rásteres continuos.

  • Diferencia de banda única: calcula la diferencia entre los valores de banda para una sola banda en cada ráster.
  • Diferencia de índice de banda: calcula un índice de banda en cada ráster y, a continuación, calcula la diferencia entre los valores de índice.

De banda de ráster y A banda de ráster

Seleccione la banda en el parámetro De ráster y la banda del parámetro A ráster para comparar en el cálculo. Este parámetro es obligatorio si Método de diferencia de banda es Diferencia de banda única.

Normalmente, se comparan valores de píxel de la misma banda en ambas imágenes. Por ejemplo, puede comparar los valores de banda Infrarroja cercana entre dos fechas. En ocasiones, los índices de banda difieren entre las dos imágenes, por ejemplo, si se comparan imágenes de Landsat 5 con imágenes de Landsat 8, por lo que es importante especificar la banda de cada imagen.

Índice de banda

Seleccione el índice de banda que desea calcular para cada ráster antes de calcular la diferencia entre ellos. Este parámetro es obligatorio si Método de diferencia de banda es Diferencia de banda única.

La lista de opciones es la misma lista de índices disponibles en la Galería de índices.

Debe especificar las bandas necesarias para generar el índice para cada imagen.

Tipo de tamaño de celda

Elija el tamaño de celda que desee usar en el ráster de salida. Si todos los tamaños de celda de entrada son iguales, todas las opciones generarán los mismos resultados.

  • Primero de: utiliza el primer tamaño de celda de los rásters de entrada.
  • Mínimo de: utiliza el tamaño de celda más pequeño de todos los rásters de entrada.
  • Máximo de: utiliza el tamaño de celda más grande de todos los rásters de entrada. Esta es la opción predeterminada.
  • Valor medio de: utiliza el tamaño de celda medio de todos los rásters de entrada.
  • Último de: utiliza el último tamaño de celda de los rásteres de entrada.

Tipo de extensión

Elija qué extensión se debería usar en el ráster de salida.

  • Primero de: utiliza la extensión del primer ráster de entrada para determinar la extensión del procesamiento.

  • Intersección de: utiliza la extensión de los píxeles superpuestos para determinar la extensión del procesamiento. Esta es la opción predeterminada.
  • Unión de: utiliza la extensión de todos los rásters para determinar la extensión del procesamiento.
  • Último de: utiliza la extensión del último ráster de entrada para determinar la extensión del procesamiento.

Haga clic en Vista previa para agregar una capa de vista previa al mapa, utilizando las opciones que especificó. Puede modificar las opciones y hacer clic en Actualizar vista previa para actualizar la capa de vista previa con las nuevas opciones.

Nota:

Las capas de vista previa se generan utilizando funciones ráster. Cuando se aleja más allá de la resolución de origen, las funciones ráster procesan los datos utilizando un tamaño de píxel remuestreado. Para asegurarse de que la vista previa tenga el aspecto que desea ver en el resultado final, acérquese a la resolución de origen de los datos. Haga clic con el botón derecho en la capa de vista previa y haga clic en Zoom a la resolución de origen.

Al hacer clic en Siguiente en el panel Diferencia de banda, el asistente calcula el índice de banda (si procede), calcula la diferencia entre los rásteres y calcula el histograma de los valores de diferencia. Este cálculo puede tardar un tiempo en función del tamaño de los datasets.

Panel Clasificar diferencia

El panel Clasificar diferencia le permite explorar las diferencias en los valores de píxel al vuelo en el mapa, y usar la información para clasificar los resultados. Por ejemplo, puede clasificar los cambios de NDVI en las clases Pérdida o Ganancia, o puede clasificar las diferencias en los valores de banda NIR para extraer un área quemada y clasificarla como clase Área quemada.

La sección Explorar diferencias ocupa la mitad superior del panel Clasificar diferencia.

La sección Explorar diferencias y el histograma interactivo
El histograma de valores de diferencia se muestra junto con controladores interactivos para visualizar el cambio. La capa Preview_Mask se agrega al mapa de forma predeterminada al abrir este panel. Al arrastrar los controladores de histograma mínimo y máximo a diferentes valores en el histograma, los valores de píxel entre el mínimo y el máximo estarán visibles en la capa Preview_Mask de su mapa. El resto de valores de píxel se enmascaran y no son visibles. Puede resultar más fácil ver los valores si el resto de capas del mapa están desactivadas. También puede modificar la simbología de la capa Preview_Mask para ver los resultados con mayor facilidad.

Para omitir la clasificación y conservar las diferencias de valores de píxel originales, desactive el parámetro Clasificar la diferencia en los valores y haga clic en Siguiente para pasar al siguiente panel.

Utilizar el histograma interactivo para generar clases

Para usar el histograma interactivo para generar clases, siga estos pasos:

  1. Arrastre los controladores mínimo y máximo en el histograma hasta que haya identificado la información de cambio que sea importante para su análisis (por ejemplo, un fuerte cambio negativo).
  2. Haga clic en el botón Agregar nueva claseInsertar flecha abajo.

    Los valores mínimo y máximo seleccionados en el histograma se agregan como valores de clase mínimo y máximo a la tabla Clasificar salida en la parte inferior del panel.

  3. En la tabla Clasificar salida, agregue el valor de clase (un entero) al campo Output y el nombre de la clase al campo Class Name. Estos campos son obligatorios.
  4. Especifique el color que se utilizará para simbolizar las clases de salida.

Generar clases manualmente

Para introducir de forma manual la información de la clase, sin utilizar el histograma interactivo, siga estos pasos:

  1. Haga clic en el botón Generar Esquema de colores.
  2. En el panel Generar, defina el parámetro Valor máximo con el número de clases que desee generar.

    Puede especificar un nombre base para las clases asociadas a cada valor único. También puede especificar el esquema de color y la transparencia de la simbología de las clases.

  3. Cuando haya terminado, haga clic en Aceptar.
  4. En la tabla Clasificar salida, escriba de forma manual los valores de los campos de diferencia Minimum y Maximum para cada clase.

Haga clic en Vista previa para agregar una capa de vista previa al mapa, utilizando las opciones que especificó. Puede modificar las opciones y hacer clic en Actualizar vista previa para actualizar la capa de vista previa con las nuevas opciones.

Nota:

Las capas de vista previa se generan utilizando funciones ráster. Cuando se aleja más allá de la resolución de origen, las funciones ráster procesan los datos utilizando un tamaño de píxel remuestreado. Para asegurarse de que la vista previa tenga el aspecto que desea ver en el resultado final, acérquese a la resolución de origen de los datos. Haga clic con el botón derecho en la capa de vista previa y haga clic en Zoom a la resolución de origen.

Panel Posprocesamiento

El panel Posprocesamiento permite suavizar los resultados y generar una salida.

ParámetroDescripción

Vecindad de suavizado

El tamaño de la vecindad focal que se va a utilizar, en filas y columnas de píxeles, para suavizar los resultados.

  • Ninguno: no se realiza el suavizado. Esta es la opción predeterminada.
  • 3 × 3: los resultados de suavizado se basan en una vecindad de 3 filas de píxeles por 3 columnas de píxeles.
  • 5 × 5: los resultados de suavizado se basan en una vecindad de 5 filas de píxeles por 5 columnas de píxeles.
  • 7 × 7: los resultados de suavizado se basan en una vecindad de 7 filas de píxeles por 7 columnas de píxeles.

Método de relleno de estadísticas

El método estadístico que se utiliza para recalcular los valores de píxel para un resultado más suave.

  • Mínimo: calcula el valor mínimo de los píxeles dentro de la vecindad.
  • Máximo: calcula el valor máximo de los píxeles dentro de la vecindad.
  • Media: calcula el valor medio de los píxeles dentro de la vecindad. Esta es la opción predeterminada.
  • Desviación estándar: calcula el valor de la desviación estándar de los píxeles dentro de la vecindad.
  • Valor medio: calcula el valor medio de los píxeles dentro de la vecindad.
  • Mayoría: calcula el valor de mayoría o el valor que se produce con mayor frecuencia de los píxeles dentro de la vecindad.
  • Minoría: calcula el valor de minoría o el valor que se produce con menor frecuencia de los píxeles dentro de la vecindad.

Guardar resultado como

Especifique el tipo de salida que se genera.

  • Dataset ráster: los resultados se guardan como un dataset ráster.
  • Plantilla de función ráster: los resultados se guardan como una plantilla de función ráster para procesamiento adicional.
  • Clase de entidad: los resultados se guardan como una clase de entidad poligonal.

Dataset de salida

El nombre del dataset de salida. Este parámetro solo está disponible si Guardar resultado como está establecido en Dataset ráster o Clase de entidad. Si la salida es un dataset ráster, especifique la extensión de archivo para generar un ráster en un formato compatible.

Para opciones de suavizado adicionales, puede utilizar la función ráster Estadísticas en los resultados. Elija Ninguno como Vecindad de suavizado en el asistente y utilice la función ráster Estadísticas en el dataset ráster de salida.

Haga clic en Vista previa para agregar una capa de vista previa al mapa, utilizando las opciones que especificó. Puede modificar las opciones y hacer clic en Actualizar vista previa para actualizar la capa de vista previa con las nuevas opciones.

Nota:

Las capas de vista previa se generan utilizando funciones ráster. Cuando se aleja más allá de la resolución de origen, las funciones ráster procesan los datos utilizando un tamaño de píxel remuestreado. Para asegurarse de que la vista previa tenga el aspecto que desea ver en el resultado final, acérquese a la resolución de origen de los datos. Haga clic con el botón derecho en la capa de vista previa y haga clic en Zoom a la resolución de origen.

Ejemplo

El siguiente ejemplo extrae la pérdida en los valores del Índice de vegetación ajustada del suelo (SAVI) entre dos imágenes de Landsat 8.

  1. Agregue las dos imágenes de Landsat 8 a su mapa.
  2. Con la capa de imagen anterior seleccionada en el panel Contenido, inicie el Asistente de detección de cambios en la pestaña Imágenes del grupo Análisis.

    En el panel Configurar, el parámetro De ráster se define como la capa ráster anterior.

  3. Establezca el parámetro A ráster como la capa ráster posterior.
  4. Haga clic en Siguiente.
  5. En el panel Diferencia de banda, configure el análisis para calcular los cambios en los valores SAVI.
    1. Cambie el parámetro Tipo de diferencia a Absoluta.
    2. Fije el parámetro Método de diferencia de bandas a Diferencia de índice de banda.
    3. En la sección Parámetros de índice de banda - De ráster, defina Índice de banda infrarroja cercana como Banda 5 y defina Índice de banda roja como Banda 4. Mantenga el valor predeterminado del factor de corrección de brillo del suelo en 0,5.
    4. En la sección Parámetros del índice de banda - A ráster, defina Índice de banda cercana a Infrarrojo como Banda 5 y defina Índice de banda roja como Banda 4. Mantenga el valor predeterminado del factor de corrección de brillo del suelo en 0,5.
    5. Cambie el parámetro Tipo de tamaño de celda a Máx. de y el parámetro Tipo de extensión a Intersección de.
  6. Haga clic en Siguiente.
  7. En el panel Clasificar diferencia, arrastre el controlador máximo del histograma a 0 y deje el controlador mínimo del histograma en el valor negativo más bajo. Haga clic en Agregar nueva clase.
  8. En la tabla Clasificar salida, defina el valor del campo Output en 1, el campo Class Name en Pérdidas y utilice el selector de color para seleccionar un color rojo.
  9. Haga clic en Siguiente.
  10. En el panel Posprocesamiento, establezca el parámetro Vecindad suavizada en 3 x 3 y el parámetro Método de relleno de estadísticas en Mayoría.
  11. En el parámetro Dataset de salida, escriba SAVI_Loss.tif. Haga clic en Ejecutar.
  12. Cuando se agrega el dataset al mapa, haga clic en Finalizar para cerrar el Asistente de detección de cambios.