Tutorial: Crear un producto de modelo digital de superficie a partir de imágenes de satélite en ArcGIS Reality for ArcGIS Pro

Disponible con una licencia Advanced.

En ArcGIS Pro, puede corregir las imágenes de satélite fotogramétricamente para eliminar las distorsiones geométricas causadas por la plataforma y el desplazamiento del terreno. Después de eliminar estas distorsiones, puede generar productos ArcGIS Reality for ArcGIS Pro, lo que incluye el modelo digital de superficie (DSM) y la malla DSM. En este tutorial, generará un DSM de alta resolución.

En primer lugar, configurará un espacio de trabajo de representación cartográfica de realidad para administrar la colección de imágenes de satélite. A continuación, realizará un ajuste de bloques, seguido de un ajuste afinado utilizando puntos de control del terreno. Para terminar, generará un DSM.

ArcGIS Pro puede procesar imágenes de satélite de muchas plataformas de sensores, siempre y cuando la orientación de la imagen se describa con un modelo de coeficientes polinómicos racionales (RPC) o un modelo de sensor riguroso. Este modelo viene típicamente incorporado en el archivo de imagen o incluido como archivo de metadatos separado.

Crear un espacio de trabajo de representación cartográfica de realidad

Un espacio de trabajo de representación cartográfica de realidad es un subproyecto de ArcGIS Pro que se dedica a los flujos de trabajo de representación cartográfica de realidad. Es un contenedor en una carpeta de proyecto de ArcGIS Pro que almacena los recursos y archivos derivados que pertenecen a una única colección de imágenes en una tarea de representación cartográfica de realidad.

Maxar Technologies ha recopilado y proporcionado las imágenes empaquetadas para este tutorial. Incluye un par de imágenes multiespectrales y pancromáticas, una tabla de puntos de control del terreno y un DEM.

1. Descargue el dataset del tutorial y guárdelo en C:\sampledata\RM_satellite_tutorial.
2. Descomprima los datos.
3. En ArcGIS Pro, cree un proyecto con la plantilla Mapa e inicie sesión en su cuenta de ArcGIS Online si es necesario.
4. En la pestaña Imágenes, en el grupo Representación cartográfica de realidad, haga clic en el menú desplegable Nuevo espacio de trabajo y elija Nuevo espacio de trabajo.
5. En la página Configuración de espacio de trabajo, escriba un nombre para el espacio de trabajo.
6. Asegúrese de que la opción Tipo de espacio de trabajo esté establecida en Representación cartográfica de realidad.
7. En la lista desplegable Tipo de datos del sensor, elija Satélite.
8. Opcionalmente, en la lista desplegable Mapa base, elija un mapa base como telón de fondo para la colección de imágenes.
9. También puede activar la casilla de verificación Permitir restablecimiento de ajuste, que ofrece la posibilidad de revertir el espacio de trabajo a un estado anterior.
10. Haga clic en Siguiente.

    

11. En la ventana Colección de imágenes, en el menú desplegable Tipo de sensor, elija GeoEye-1.
12. En Carpeta que contiene imágenes, haga clic en el botón Examinar, busque la carpeta de datos del tutorial en su equipo y elija la carpeta de imágenes (GeoEye_NAD83_UTM10N).
13. En Referencia espacial, haga clic en el botón Examinar .
14. En la ventana Referencia espacial, en XY actual, establezca la referencia espacial en NAD 83 UTM Zone 10N (WKID 26910).
15. En Z actual, expanda Sistema de coordenadas verticales y en Relacionado con gravedad, expanda Mundo y elija EGM96 Geoid.
16. Haga clic en Aceptar para cerrar la ventana Referencia espacial y en Siguiente.

    

17. En la ventana Cargador de datos, en Fuente de elevación, elija DEM. En DEM, vaya al DEM proporcionado con el dataset del tutorial.

    Nota:

    Este DEM se utilizará para permitir el proceso de ajuste de bloques.

18. La mayoría de datos de elevación utilizan alturas ortométricas, así que tendrá que aplicar una corrección geoide. En Corrección de geoide, asegúrese de que EGM96 esté seleccionado.
19. En Plantilla de procesamiento, elija Pancromático.
20. Expanda Opciones avanzadas.
21. Asegúrese de que la opción Estimar estadísticas está activada.
22. Expanda Preprocesamiento y active la casilla situada junto a Calcular estadísticas.

    Esto hará que se muestren las opciones de Calcular estadísticas.

23. Defina los valores de Número de columnas para omitir y Número de filas para omitir como 1.
24. Acepte el resto de valores predeterminados y haga clic en Finalizar.

    

    Una vez creado el espacio de trabajo, se muestran las imágenes y las huellas de imagen en el mapa. También se agrega una categoría Reality Mapping al panel Contenido, donde se almacenarán los datos de imágenes de origen y los productos de representación cartográfica de realidad derivada.

    La visualización inicial de las imágenes en el espacio de trabajo confirma que se proporcionaron todas las imágenes y los metadatos necesarios para iniciar el espacio de trabajo. Las imágenes no se alinearon ni se ajustaron, por lo que el mosaico podría no parecer geométricamente correcto.

    Se agregará una pestaña Representación cartográfica de realidad nueva al menú principal ArcGIS Pro. Al hacer clic en esta pestaña aparece una serie de herramientas y flujos de trabajo dedicada a la representación cartográfica de realidad. En la categoría Producto, ninguno de los botones está disponible porque todavía no se han ajustado las imágenes.

    

    Imagen de satélite © 2020 Maxar Technologies

Eliminar NoData

A continuación, eliminará las áreas NoData que circundan las imágenes de satélite válidas.

1. En el menú principal, haga clic en la pestaña Análisis y en Herramientas en el grupo Geoprocesamiento.

   Se abre el panel Geoprocesamiento.

2. En la ventana de búsqueda del panel Geoprocesamiento, escriba crear huella.

   Aparece el cuadro de diálogo Crear huellas .

3. En Dataset de mosaico, elija la colección de imágenes en la lista desplegable.
4. Acepte el resto de valores predeterminados y haga clic en Ejecutar.

   Al final, las áreas NoData se eliminan y la huella de imagen se actualiza como corresponde.

Ajuste de bloques

Después de crear un espacio de trabajo de representación cartográfica de realidad, el paso siguiente consiste en realizar un ajuste de bloques con las herramientas de los grupos Ajustar y Perfeccionar. El ajuste de bloques calculará en primer lugar los puntos de enlace, que son puntos comunes en áreas de superposición de imágenes. Los puntos de enlace se utilizarán para calcular la orientación de cada imagen, conocida como orientación exterior en fotogrametría.

1. En la pestaña Representación cartográfica de realidad, en el grupo Ajustar, haga clic en Ajustar .
2. En la ventana Ajustar, en Tipo de transformación, elija RPC. La transformación de coeficientes polinómicos racionales (RPC) se aplicará al ajuste, que se utiliza con las imágenes de satélite que contienen información de RPC en los metadatos.
3. En la ventana Ajustar, en Umbral de punto de equivocación, defina un valor de 3. Los puntos de enlace que presentan un error residual superior a este valor no se usarán en el cálculo del ajuste. La unidad de medida son píxeles.
4. Active la casilla situada junto a Volver a proyectar puntos de enlace.

   Esto garantizará el cálculo de las coordenadas del mapa de puntos.

5. En Coincidencia de puntos de enlace del panel Adjutar, asegúrese de que la configuración de parámetros coincida con el ejemplo mostrado.

   

6. Haga clic en Ejecutar para realizar un ajuste de bloques.
7. Una vez completado el ajuste, active la capa Puntos de enlace en el panel Contenido para ver la distribución de los puntos de enlace generados en el mapa. La distribución de los puntos de enlace puede diferir de la mostrada a continuación.

   

8. En el archivo de registro puede ver residuales de los puntos de enlace o informes de precisión. En la pestaña Representación cartográfica de realidad, en el grupo Revisar, haga clic en Registros para acceder al archivo de registros. Los residuales de puntos de enlace se muestran en la fila denominada RMSE_Tie_Image(x,y). Las unidades para el punto de enlace RSME están en píxeles.

   El residual puede ser distinto de lo que se muestra abajo.

   

Tras el ajuste inicial se activan dos botones en la categoría Producto. Estos botones resaltan los productos de imagen que se pueden generan. Antes de generar productos, se utilizarán los puntos de control del terreno para mejorar la precisión absoluta de las imágenes.

Agregar puntos de control del terreno

Los puntos de control del terreno (GCP) son puntos con coordenadas terrestres x,y,z conocidas. A menudo se obtienen a partir de reconocimientos terrestres o datos existentes y se utilizan para garantizar que las imágenes se georreferencien con precisión en el sistema de coordenadas terrestres. El ajuste de bloques se puede aplicar sin los GCP y aun así garantiza una precisión relativa, pero agregar los GCP aumenta la precisión absoluta de las imágenes ajustadas. Si no dispone de GCP de un estudio del terreno, pero dispone de una capa ráster georreferenciada (dataset ráster, dataset de mosaico o servicio de imágenes), puede agregarla como referencia para calcular los GCP.

Los GCP seleccionados en una capa de referencia y almacenados en un archivo de texto también pueden importarse y utilizarse para mejorar la precisión del ajuste. En este tutorial, utilizará este método de importación para agregar GCP al proyecto.

Importar GCP

Para importar GCP, siga estos pasos:

1. En la pestaña Representación cartográfica de realidad, en el grupo Perfeccionar, haga clic en Administrar GCP para abrir el Administrador de GCP.
2. En la ventana Administrador de GCP, haga clic en el botón Importar GCP .
3. En la ventana Importar GCP, busque y seleccione el archivo GCP (Vancouver_NAD83-UTM10N.csv). Haga clic en Aceptar.
4. En Definir referencia espacial GCP, haga clic en el botón Examinar . En XY actual, expanda Capas y elija NAD 1983 UTM Zone 10N. Acepte todos los demás valores predeterminados y haga clic en Aceptar para aceptar los cambios y cerrar la ventana Referencia espacial.

   El sistema de coordenadas verticales (VCS) o Z actual, no se ha definido porque el modelo de elevación digital utilizado para extraer los valores de altura correspondiente a los GCP no tienen un VCS definido. Si el DEM que ha utilizado tiene un VCS definido, en Z actual (VCS) se habrán establecido las coordenadas correspondientes.

5. En Transformaciones geográficas, acepte los valores predeterminados. No se requieren transformaciones porque la referencia espacial horizontal de los GCP y el espacio de trabajo es la misma.
6. Asegúrese de que las asignaciones de campo son correctas.
7. Haga clic en el botón de exploración situado debajo de Ubicación de fotos de GCP, vaya a la carpeta que contiene las imágenes de las ubicaciones de GCP y selecciónela. Haga clic en Aceptar.
8. En Extensión de archivo fotográfico, elija PNG en la lista desplegable.

   

9. Haga clic en Aceptar para importar los GCP.

   Una vez importados los GCP, la tabla del Administrador de GCP se rellena.

   

   En el Administrador de GCP, el primer GCP de la lista se selecciona de manera predeterminada y en la sección Vista previa se muestra una vista previa de la imagen en la que se encuentra el GCP.

10. Para agregar puntos de enlace para el GCP seleccionado, en la sección Vista previa haga clic en el botón Ver foto de GCP para mostrar el chip de imagen de GCP. Use la rueda del ratón para acercar el chip de imagen para ver la ubicación GCP indicada mediante una flecha roja.
11. Haga clic en el botón Agregar punto de enlace para agregar un punto de enlace en el visor de imágenes de cada imagen. El algoritmo de coincidencia de imágenes calcula automáticamente cuando sea posible, los puntos de enlace de otras imágenes, pero la precisión de cada punto de enlace debe verificarse. Si el punto de enlace no se identifica automáticamente, agregue manualmente el punto de enlace seleccionando la ubicación adecuada en la imagen.
12. Repita los pasos 10 y 11 para seleccionar y agregar puntos de enlace para los GCP restantes.
13. Una vez que se haya medido cada GCP con puntos de enlace, seleccione el punto PP_GCP02 y haga clic con el botón derecho para que se convierta en un punto de verificación. Este punto no se utilizará en el proceso de ajuste debido a la incertidumbre en su ubicación exacta.
14. Después de agregar GCP y puntos de verificación, hay que ejecutar de nuevo el ajuste para incorporar estos puntos. Haga clic en Ajustar.

    

Revisar los resultados del ajuste

Los resultados de calidad del ajuste se pueden visualizar en el Administrador de GCP analizando los residuales de cada GCP. Los residuales representan la diferencia entre la posición medida y calculada de un punto. Se miden en las unidades del sistema de referenciación espacial del proyecto. Después de completar el ajuste con GCP, se agregan tres campos nuevos (dX, dY y dZ) a la tabla Administrador de GCP y se muestran los residuales para cada GCP. La calidad del ajuste entre el bloque ajustado y el sistema de coordenadas del mapa se puede evaluar utilizando estos valores. El error cuadrático medio (RMSE) de los residuales se puede visualizar expandiendo la sección Vista general residual del Administrador de GCP.

También se proporcionan estadísticas de precisión adicionales en el informe de ajuste. Para generar el informe, en la pestaña Representación cartográfica de realidad, en el grupo Revisar, haga clic en Informe de ajuste.

Generar un modelo digital de superficie (DSM)

Una vez completado el ajuste de bloques, los productos de imágenes 2D se pueden generar utilizando las herramientas del grupo Producto de la pestaña Representación cartográfica de realidad. Se pueden generar varios productos simultáneamente utilizando el asistente Múltiples productos o individualmente seleccionando la herramienta de producto correspondiente en el grupo Producto. Los tipos de productos que pueden generarse a partir de los datos del sensor incluyen un DSM y una malla DSM. En este tutorial solo se generará un DSM.

El DSM es un primer producto de superficie que incluye las elevaciones de árboles, edificios y otras entidades que están encima del suelo.

Para generar un DSM a partir de imágenes de satélite mediante el asistente de DSM, siga los pasos que se describen a continuación.

1. En la pestaña Representación cartográfica de realidad, en el grupo Producto, haga clic en el botón DSM para abrir el asistente Productos de representación cartográfica de realidad.
2. Haga clic en Configuración avanzada compartida para abrir el cuadro de diálogo donde puede definir los parámetros que afectarán a los productos de Representación cartográfica de realidad que se generarán.

   Los valores de Calidad y Escenario se establecen automáticamente. No se deben cambiar para garantizar un rendimiento y una calidad óptimos del producto. Si desea generar un producto de resolución reducida, el valor de Calidad se puede reducir. Consulte Configuración avanzada compartida para obtener más información sobre el impacto de varios ajustes de calidad en la generación de productos.

3. Confirme que Calidad esté definida como Ultra.
4. Acepte los valores predeterminados de Tamaño de píxel para generar el DSM con la resolución de la imagen de origen.
5. En Límite de producto, acepte la configuración predeterminada para que se procese el área del proyecto completa.

   Se recomienda introducir un límite de producto por los siguientes motivos:

   • Definir la extensión de salida adecuada: cuando no se define un límite de producto, la aplicación define automáticamente una extensión en función de varios parámetros de dataset que podrían no coincidir con la extensión del proyecto.
   • Reducir el tiempo de procesamiento: si la extensión requerida del producto es menor que la extensión de la colección de imágenes, definir un límite de producto reduce la duración del procesamiento y recorta automáticamente la salida a la extensión del límite.

6. Para Masas de agua, seleccione SatTut_WaterBody.shp en la ubicación de datos del tutorial.
7. En Carpeta de procesamiento, haga clic en el botón Examinar y navegue hasta una ubicación en el disco que tenga un almacenamiento disponible mínimo de 10 veces el tamaño total de las imágenes que se están procesando.

   En este caso, se recomiendan 40 GB o más de almacenamiento disponible.

   Nota:

   La carpeta de procesamiento almacena los archivos temporales generados durante el procesamiento de Reality. Se recomienda que la carpeta de procesamiento se encuentre en una unidad rápida con gran espacio de almacenamiento disponible.

8. Acepte el resto de valores predeterminados y haga clic en Aceptar.

   El cuadro de diálogo Configuración avanzada de productos se cierra y le devuelve al panel Configuración de generación de productos del asistente Productos de representación cartográfica de realidad.

9. Haga clic en Siguiente para ir al panel Configuración de DSM. Asegúrese de que la configuración coincide con lo siguiente:
   1. Tipo de salida: Mosaico
   2. Formato: Formato de ráster de nube
   3. Compresión: Ninguna
   4. Remuestreo: Bilineal
10. Haga clic en Finalizar para iniciar el proceso de generación de productos.

    Una vez completado el procesamiento, el producto DSM se agrega al panel Contenido en la categoría Productos de datos y la vista de mapa 2D. También se agrega al panel Catálogo del contenedor Reality Mapping y a la carpeta de DEMs.

    

Nota:

El sistema de coordenadas verticales (VCS) predeterminado del DSM de salida es WGS84 (altura elipsoidal). Si necesita convertir la salida en una proyección o sistema de coordenadas diferente, utilice la herramienta Proyecto.

En este tutorial, creó un espacio de trabajo de representación cartográfica de realidad para imágenes de satélite y utilizó herramientas en la pestaña Representación cartográfica de realidad para aplicar un ajuste fotogramétrico con puntos de control del terreno. Luego ha utilizado el asistente de DSM de representación cartográfica de realidad para generar un DSM de alta resolución.

Maxar Technologies ha adquirido y proporcionado las imágenes de satélite utilizadas en este tutorial.