Implementar un algoritmo utilizando funciones en un dataset de mosaico

Complejidad: Intermedio Requisitos de datos: ArcGIS.com

En este flujo de trabajo se explica cómo combinar funciones en un dataset de mosaico para crear un NDVI. Hay disponible una función NDVI que aplicará una escala y un desplazamiento para permitir su visualización; sin embargo, si desea obtener valores NDVI brutos, esto es, comprendidos entre -1,0 y 1,0, tendrá que implementar los pasos siguientes.

Aplicará la función al dataset de mosaico, no a los rásteres individuales contenidos en el dataset de mosaico.

El algoritmo para el cálculo del NDVI es el siguiente:

NDVI = ((IR - R)/(IR+R))

dónde

IR = banda infrarroja

R = banda roja

Se obtendrá como resultado la siguiente disposición en la ficha Función.

El flujo de trabajo se divide en función de cada entrada y de cada acción. A continuación se indica el orden que utilizaría para calcular el algoritmo, pero al encadenar las funciones, empezará por el resultado e irá hacia las entradas, como puede ver en el diagrama anterior.

1. Defina las entradas de banda: IR y R
2. Sustraiga una de la otra: IR - R
3. Súmelas: IR + R
4. Divídalas: (IR - R) / (IR + R)

Preparar el dataset de mosaico

Va a aplicar las funciones a un dataset de mosaico. Por consiguiente, debe crear un dataset de mosaico y agregar todos sus datos ráster. Se supone que los datos tienen asignado el mismo número de banda para las bandas infrarroja y roja. En este ejemplo, suponga que la banda infrarroja es la banda 1 y la banda roja es la banda 2.

Pasos:

1. Cree un dataset de mosaico.
2. Agregue sus datos ráster al dataset de mosaico.

   Asegúrese de generar las vistas generales.

Encadene las funciones aritméticas

Pasos:

1. Agregue su dataset de mosaico a la pantalla de ArcMap.

   Puede agregar el dataset de mosaico ahora o al final. Al hacerlo, verá el resultado en la capa Imagen. Si acaba de crear el dataset de mosaico, puede que ya esté abierto en ArcMap.

2. Abra la ventana Catálogo y desplácese hasta el dataset de mosaico.
3. Haga clic con el botón derecho en el dataset de mosaico y, a continuación, haga clic en Propiedades.
4. Haga clic en la ficha Funciones.

   Verá la siguiente cadena de funciones, formada únicamente por la función Mosaico. La función Mosaico aparece de forma predeterminada porque estamos trabajando con un dataset de mosaico y todos los rásteres se están incorporando al mosaico. El algoritmo se aplicará tras la operación de creación del mosaico.

5. Haga clic con el botón derecho en la función Mosaico y, a continuación, haga clic en Insertar > Función aritmética.

   La inserción tiene lugar encima de la función en la que haya hecho clic.

   

   Se abre el cuadro de diálogo Propiedades de funciones de ráster, que le permite establecer las opciones para esta función.

6. Haga clic en la flecha de la lista desplegable Operación y, a continuación, en Dividir.

   Las entradas se quedarán tal como están. La entrada se agregará a la cadena de funciones al cerrar este cuadro de diálogo.

7. Haga clic en la ficha General.
8. Escriba Dividir en el cuadro de texto Nombre de la función.
9. Haga clic en la flecha de la lista desplegable Tipo de píxel de salida y, a continuación, haga clic en 32 bit flotante.

   Este algoritmo da como resultado valores entre -1 y 1. Deberá especificar el tipo de píxel; de lo contrario, puede acabar obteniendo valores de píxel iguales a las entradas, por ejemplo de 8 bits.

10. Haga clic en Aceptar.

    La función Dividir se agrega a la parte superior de la cadena de funciones y se agrega también otra función Mosaico, puesto que esta función aritmética requiere dos entradas.

11. Haga clic con el botón derecho en la función Mosaico situada en la parte superior de la cadena de funciones y, a continuación, haga clic en Insertar > Función aritmética.
12. Haga clic en la flecha de la lista desplegable Operación y, a continuación, en Resta.

    Las entradas se quedarán tal como están. La entrada se agregará a la cadena de funciones al cerrar este cuadro de diálogo.

13. Haga clic en la ficha General.
14. Escriba Resta en el cuadro de texto Nombre de la función.
15. Haga clic en la flecha de la lista desplegable Tipo de píxel de salida y, a continuación, haga clic en 32 bit flotante.
16. Haga clic en Aceptar.

    La función Resta se agrega después de la función Dividir, y se agregan debajo dos entradas denominadas Función Mosaico.

17. Haga clic con el botón derecho en la función Mosaico situada en la parte inferior de la cadena de funciones y, a continuación, haga clic en Insertar > Función aritmética.
18. Haga clic en la flecha de la lista desplegable Operación y, a continuación, en Suma.

    Las entradas se quedarán tal como están. La entrada se agregará a la cadena de funciones al cerrar este cuadro de diálogo.

19. Haga clic en la ficha General.
20. Escriba Suma en el cuadro de texto Nombre de la función.
21. Haga clic en la flecha de la lista desplegable Tipo de píxel de salida y, a continuación, haga clic en 32 bit flotante.
22. Haga clic en Aceptar.

    La función Suma se agrega en la parte inferior en el mismo nivel que la función Resta, y debajo se agregan dos entradas más denominadas Función Mosaico.

Definir las entradas de las funciones aritméticas

Ya ha agregado las funciones que se utilizarán para realizar los cálculos aritméticos; a continuación, debe identificar las entradas. Esas entradas son las bandas IR y R (bandas 1 y 2 respectivamente, tal y como se definió anteriormente). Para ello, utilice la función Extraer bandas.

Pasos:

1. Haga clic con el botón derecho en la Función Mosaico ubicada debajo de la función Resta y, a continuación, haga clic en > Función Extraer banda.

2. Escriba 1 en el cuadro de texto ID de banda.
3. Haga clic en la ficha General.
4. Escriba Extraer IR en el cuadro de texto Nombre de la función.
5. Haga clic en Aceptar.

   Se agrega a la cadena de funciones la función Extraer IR, a la que a su vez agrega una función Mosaico como entrada.

6. Haga clic con el botón derecho en la otra Función Mosaico ubicada debajo de la función Resta y, a continuación, haga clic en Insertar > Función Extraer banda.
7. Escriba 2 en el cuadro de texto ID de banda.
8. Haga clic en la ficha General.
9. Escriba Extraer Red en el cuadro de texto Nombre de la función.
10. Haga clic en Aceptar.

    Se agrega a la cadena de funciones la función Extraer Red, a la que a su vez agrega una función Mosaico como entrada.
11. Repita estos pasos para las entradas situadas debajo de la función Suma.

    Obtendrá finalmente una cadena de funciones parecida a la mostrada al principio de este tema.

12. Haga clic en Aceptar para cerrar el cuadro de diálogo Propiedades de dataset de mosaico.

Cuando haya terminado de agregar estas funciones, puede obtener una vista previa de la imagen en mosaico. Puede que necesite calcular las estadísticas del dataset de mosaico, puesto que se ha modificado para contener una salida de 32 bits. Puede utilizar para ello la herramienta Calcular Estadísticas.

Puede examinar los valores de píxel mediante las herramientas Identificar o Inspección de Píxeles.