El proyecto inició en el mes de enero de 2015 con una capacitación rápida en la que se repasaron las características básicas del equipo GNSS que posee la ECG y los tipos de medición que pueden realizarse con el mismo, además, se revisó la licencia del software de procesamiento, se actualizó el firmware de los receptores y el programa Survey Controller y la configuración de los métodos de observación definidos en el software del controlador.

La ECG ya contaba con un equipo Trimble® R4 GNSS (ver figura 1), con el software de campo Survey Controller, instalado en el controlador TSC2 y una licencia del programa de post-procesamiento Trimble® Business Center (TBC, versión 2.5).

Figura 1:
equipo Trimble® R4 GNSS y controlador TSC2.
Fuente: Foto de equipo Trimble® R4 GNSS y controlador TSC2 de la Escuela Ciencias Geográficas, 2015.

Según la hoja de datos de Trimble®, el equipo R4 rastrea simultáneamente señales de satélite GPS L1, L2, L2C y QZSS, además, permite actualizar a GLONASS, Galileo, y BeiDou (COMPASS), facilitando mediciones GNSS con precisión horizontal 3 mm +0,1 ppm RMS y vertical de 3,5 mm + 0,4 ppm RMS en estáticos de alta precisión, con una precisión horizontal de 8 mm + 1 ppm RMS y vertical de 15 mm + 1 ppm RMS en redes RTK.

Este equipo había sido utilizado únicamente en prácticas de algunos cursos del Programa de Maestría en Sistemas de Información Geográfica y Teledetección, por lo que se consideró que estaba subutilizado.

El proyecto “Aplicación de la herramienta GPS en la enseñanza de Cartografía” realizó las primeras pruebas de medición en el Campus Omar Dengo durante el mes de marzo de 2015. El ejercicio fue de gran aprendizaje para identificar problemas en el trabajo de campo y en el procesamiento de los datos, entre los que se encontraban que los sitios de medición seleccionados resultaron no ser los ideales, dificultades en la descarga de los datos por inexperiencia en el manejo de los formatos, y principalmente, en el manejo del software de post-procesamiento, debido a desconocimiento de nociones básicas de geodesia.

Ante las dificultades enfrentadas, se decidió solicitar capacitación y asesoramiento al Centro Nacional de Procesamiento de Datos Geoespaciales (CNPDG) de la Escuela de Topografía Catastro y Geodesia (ETCG) de la UNA.

Los datos generados durante la medición realizada en el mes de marzo fueron de gran utilidad para esta capacitación, permitiendo que los participantes del proyecto se concentraran en aprender sobre datos de alta precisión y equipos GNSS, definición de redes para la medición, la planificación del trabajo de campo, las mediciones y post-proceso de los datos brutos.

Se redefinieron las redes de medición en el Campus Omar Dengo y en el distrito primero del cantón Central de Heredia, considerando la cantidad y características del equipo disponible, las distancias y desplazamientos entre los sitios de medición, la cantidad de personas que participarían del trabajo de campo, el uso que se le daría a las coordenadas en el resto del proyecto, entre otros aspectos. Se generaron una serie de plantillas a utilizar en el trabajo de campo como fichas para la planificación y control del trabajo de campo y las mediciones.

En el 2015 se llevaron a cabo dos nuevos levantamientos de datos, utilizando mediciones de tipo estático, el de la red del Campus Omar Dengo que se realizó en los meses de julio y agosto, y el de la red del distrito Heredia efectuada en agosto. Se levantaron un total de 14 puntos de control distribuidos de la siguiente forma:

• Red Heredia: se establecieron siete puntos de medición en el distrito primero del cantón Central de Heredia. Se definió como base un punto ubicado en El Fortín y los otros se ubicaron en el Colegio Humanista de la UNA, la Escuela Cleto González Víquez, el Polideportivo de Barrio Fátima, el Palacio de los Deportes, el Polideportivo de San Francisco y el parque infantil de la Urbanización Miraflores (ver mapa 1). El trabajo de campo se realizó los días 24 y 31 de julio de 2015.

• Red Campus Omar Dengo: también se definieron 7 sitios de medición, con una base establecida en la Plaza de la Diversidad. Los restantes seis puntos de ubicaron en el parqueo de la Facultad de Ciencias Sociales, el Centro de Estudios Generales, el puente de la Escuela de Ciencias Geográficas, el parqueo de la Escuela de Ciencias Biológicas, el parqueo de la Escuela de Música y en la Estación Meteorológica de la Escuela de Física. (ver mapa 2). Este levantamiento se realizó los días 7 y 21 de agosto de 2015.

Para ambas redes se realizaron mediciones bajo el método de observación estático, estableciendo para la base tiempos de medición de aproximadamente 8 horas y de 30 minutos para cada uno de los otros seis puntos. Los levantamientos se realizaron en dos días distintos, pues las redes fueron configuradas de forma tal que en un día se midiera desde la base hacia los puntos periféricos, y que en el segundo día las mediciones se realizan entre puntos.

Los datos brutos obtenidos fueron descargados en el controlador tras cada medición, y posteriormente en una computadora, en la cual fueron procesados.

Mapa 1:
Distribución de los puntos de control de la Red Heredia

Mapa 2:
Distribución de los puntos de control de la Red Campus Omar Dengo.

El post-procesamiento de los datos se realizó utilizando el software Trimble® Business Center (TBC), versión 2.5, que funciona mediante una llave de hardware HASP® conectada a un puerto USB del equipo.

Se trabajaron de forma separada para los datos de la Red Heredia y la Red Campus Omar Dengo, es decir, el procedimiento que a continuación se detalla se ejecutó en dos ocasiones.

Se inició con la creación de un nuevo proyecto en el programa TBC y la definición del sistema de coordenadas, coincidente con el sistema en el que se encontraban los datos adquiridos en el levantamiento. En el caso de este proyecto los datos fueron adquiridos en el sistema CRTM05, cuyos parámetros estaban incluidos en el software del controlador, no así en el TBC.

De seguido se realizó la importación de los datos brutos hacia el programa, haciendo una conexión a la carpeta que contenía los datos y seleccionando los archivos de interés (figura 2) (formato *.T01, *.T02, ambos nativos de Trimble®).

Figura 2:
Resumen de la importación de los datos brutos del 24 de julio de 2015, correspondientes a la Red Heredia.
Fuente: Elaboración propia a partir del resumen de importación de datos del programa TBC.

Durante el proceso de importación se prestó atención a aspectos reseñados en las cejillas Punto, Antena y Receptor, pues enlistan los archivos importados y sus características respectivas, ofreciendo la posibilidad de detectar errores o faltantes de información, y de pensar aquí en formas de solventarlo, mediante nuevos levantamientos o utilizando datos de estaciones de apoyo.

El producto de las mediciones es una red estática que “describe la geometría y el orden en el que se organizan y procesan las líneas base GPS capturadas mediante las técnicas estática y FastStatic […] El resultado es un conjunto de líneas base que se deriva utilizando coordenadas iniciales precisas” (Trimble®, 2005: 103).

Figura 3:
Puntos y líneas base de la Red Heredia desplegados en la vista del TBC.
Fuente: Elaboración propia a partir de la visualización de los puntos y las líneas base en el programa TBC.

Lo que se busca con el post-procesamiento es corregir los datos brutos generados en el trabajo de campo, a partir de otros datos que sean aún más precisos y confiables. Estos datos provienen de las estaciones de medición continua más próximas al área de estudio, para el caso de Costa Rica se cuenta con información de las estaciones del Instituto Geográfico Nacional (IGN) y de la red SIRGAS.

A continuación en la tabla 1 se detallan las estaciones de medición continúa utilizadas en este proyecto:

Tabla 1:
Estaciones de medición continua utilizadas en el proyecto.

Fuente: Elaboración propia con datos de estaciones SIRGAS de operación continua, http://www.sirgas.org/es/sirgas-con-network/stations/station-list/

Las coordenadas de cada estación para las fechas de levantamiento fueron descargadas utilizando el sistema SpiderWeb del Registro Nacional y mediante un servicio FTP que ofrece el CeNAT para las estaciones SIRGAS de Costa Rica.

Los datos de las estaciones de medición continua deben corresponder a las fechas de las mediciones realizadas en campo, por lo que fue necesario conocer el día y semana GPS de los levantamientos. Esta consulta se realizó en el calendario GNSS (Jahic, 2016). Por ejemplo, (figura 4) una de las mediciones en el campo se llevó a cabo el 24 de julio de 2015, fecha que corresponde al día GPS 205 y a la semana GPS 1854.

Figura 4:
Visualización de consulta al Calendario GNSS.
Fuente: Calendario GNSS http://www.gnsscalendar.com/

Tabla 2:
Fechas de las mediciones y sus correspondientes días según calendario GNSS, junto con el nombre de cada archivo descargado de las estaciones de medición continúa.

Fuente: Elaboración propia a partir de datos descargados.

A continuación en la figura 5 se muestra un resumen de la importación de los datos de estaciones de medición continúa del 24 y 31 de julio de 2015, correspondientes a la Red Heredia y en la figura 6 se presentan los puntos y líneas base de la Red Heredia y estaciones de medición continua desplegados en la vista del TBC.

Figura 5:
Resumen de la importación de los datos de estaciones de medición continúa del 24 y 31 de julio de 2015, correspondientes a la Red Heredia.
Fuente: Elaboración propia a partir del resumen de importación de datos del programa TBC.

Figura 6:
Puntos y líneas base de la Red Heredia y estaciones de medición continua desplegados en la vista del TBC.
Fuente: Elaboración propia a partir de la visualización de los puntos y las líneas base en el programa TBC.

Adicionalmente en la tabla 3, y con el mismo objetivo de incrementar la precisión de las coordenadas de los sitios o puntos medidos, se incorporaron las efemérides al post-procesamiento. Las efemérides describen la órbita de cada satélite por un espacio de tiempo, calculada por el satélite en base a una serie de parámetros y ecuaciones (García, 2008). Las mismas pueden ser de navegación (transmitidas) o precisas y se definen como:

• Efemérides de navegación: posiciones teóricas de los satélites, las cuales se insertan dentro del mensaje de navegación, siendo posible que sean transmitidas en tiempo real.

• Efemérides precisas: describen las posiciones reales de los satélites, obtenidas a partir del seguimiento continuo de las estaciones del segmento control. Tales efemérides son obtenidas a partir de un post-procesado y están disponibles unos días más tarde. Son más exactas que las de navegación, pues se corresponden con posiciones reales y no teóricas de los satélites (Comisión Interdepartamental de Estadística y Cartografía de Andalucía, 2016: 45)

Para cada día de medición se identificó y descargó el archivo de efeméride precisa, también denominado órbita final.

Tabla 3:
Fechas de las mediciones y sus correspondientes días según calendario GNSS, junto con el nombre de cada archivo descargado de las órbitas precisas.

Fuente: Elaboración propia a partir de datos descargados.

Las órbitas finales (IGS, 2017), se descargaron en formato sp3 (automáticamente legible por el TBC) y se almacenaron en una carpeta desde la que se importaron.

En resumen, en el post-procesamiento se utilizaron los datos GNSS resultado de las mediciones, los puntos con las coordenadas de las estaciones de medición continua y las efemérides precisas. Como se mencionó unos párrafos atrás, la finalidad de realizar este post-procesamiento es aumentar la precisión de las coordenadas de los sitios medidos, lo que se logró mediante la aplicación de dos herramientas del TBC: el procesamiento de líneas base y el ajuste de red.

De acuerdo con la Guía del usuario Trimble® Business Center (2005), el procesamiento de líneas base busca superposiciones en los tiempos de medición para calcular la media de las coordenadas, determinando las coordenadas de más alta calidad para cada punto del proyecto.

Tras revisar, corregir y aceptar los resultados del procesamiento de líneas base, y antes de continuar con el ajuste final de la red, se editaron las coordenadas de las de las estaciones de medición continua utilizadas (AACR, CRCP, ETCG y RIDC) a partir de las soluciones semanales finales de las estaciones pertenecientes a la red SIRGAS que pueden ser consultadas y descargadas en formato *.CRD (SIRGAS, 2017).

Las coordenadas semanales del marco de referencia SIRGAS (red continental + redes nacionales de referencia) son obtenidas de la combinación de las soluciones individuales proporcionadas por los Centros de Procesamiento SIRGAS (SIRGAS, 2017), entre los que se encuentra el CNPDG - UNA, y “detallan coordenadas de precisión milimétrica (asociadas a una época específica de referencia) y sus cambios a través del tiempo (velocidades de las estaciones) para cada estación, lo cual garantiza su orientación permanente dentro del mismo sistema de coordenadas al que se refieren los satélites GNSS” (IGAC, 2016).

Al tratarse de soluciones semanales es necesario conocer a que semana GPS corresponde el día en que se realizó la medición, y como se mencionó unas páginas atrás, esto se determinó consultando el calendario GNSS. (ver tabla 4)

Tabla 4:
Fechas de las mediciones y sus correspondientes días y semana según calendario GNSS.

Fuente: Elaboración propia a partir del Calendario GNSS, http://www.gnsscalendar.com/

Las coordenadas que ofrece esta solución semanal están dadas en coordenadas cartesianas que deben pasar una serie de transformaciones para llegar a las coordenadas elipsoidicas, para WGS84 que se puedan ingresar como coordenadas locales. Este proceso se realizó con la ayuda de las hojas de cálculo para transformaciones del país elaboradas por el CNPDG,

Los datos descargados de las soluciones semanales de SIRGAS están en coordenadas cartesianas que deben transformarse a coordenadas elipsoídicas para WGS84, utilizando la hoja de cálculo denominada “xyz_flh_wgs84” (CNPDG, 2017).

A continuación en tabla 5 y 6, se muestran un ejemplo de coordenadas obtenidas a partir de la solución semanal SIRGAS para cuatro estaciones de medición continua, correspondientes a la semana GPS 1854.

Tabla 5:
Coordenadas cartesianas de la solución semanal SIRGAS para la semana 1854.

Fuente: Coordenadas semanales de las estaciones SIRGAS-CON, http://www.sirgas.org/es/sirgas-con-network/coordinates/weekly-positions/

Tabla 6:
Coordenadas elipsoídicas para WGS84 para la semana 1854.

Fuente: Elaboración propia a partir de coordenadas semanales de las estaciones SIRGAS-CON, http://www.sirgas.org/es/sirgas-con-network/coordinates/weekly-positions/

Estas coordenadas se incluyeron antes del ajuste de red con el fin de aumentar el nivel de precisión de las estaciones de medición continua y con esto las coordenadas finales que se obtendrán para los datos brutos tras el post-procesamiento. Se seleccionaron una a una cada estación y se digitaron en el cuadro de Propiedades correspondiente. Se establecieron como coordenadas locales y se indicó que poseen calidad de control.

Para finalizar el proceso de cálculo de posiciones más precisas se debe hacer el ajuste de la red, proceso que da como resultado una coordenada adicional para cada punto, misma que pasa a ser el valor final del punto. Según la Guía del usuario Trimble® Business Center (2005), la herramienta realiza un ajuste de mínimos cuadrados que asegura que no hay errores sistemáticos o severos en los datos, sino más bien errores pequeños, aleatorios y distribuidos adecuadamente.

Terminado el proceso de ajuste de red, se observó en la vista del TBC una serie de elipses de error que muestran gráficamente el error posicional (magnitud y dirección del error) de cada punto ajustado en la red. (ver figura 7)

Figura 7:
Puntos, líneas base y elipses de error de la Red Heredia y estaciones de medición continua, con órbitas finales, tras el ajuste de red.
Fuente: Elaboración propia a partir de la visualización del resultado del ajuste de red en el programa TBC.

El ajuste de red ofrece también un informe en el que se muestra información general del proyecto (nombre, fecha, sistema de coordenadas), y se indican elementos estadísticos del ajuste, los valores de las coordenadas ajustadas (de cuadrícula, geodésicas y cartesianas), los errores estimados para cada punto, dimensiones de las elipses de error y una descripción del ajuste realizado a cada línea base.

Se decidió que las coordenadas obtenidas tras el post-procesamiento eran satisfactorias para el proyecto por lo que los puntos se exportaron a formato SHP, utilizando una herramienta que ofrece el TBC, lo que permitió contar con 14 puntos de control de alta precisión para el área de estudio.

Sin embargo, a esta altura del proyecto, resulto evidente que la cantidad de puntos adquiridos en la Red Heredia y Campus Omar Dengo no serían suficientes para la georreferencia y ortorectificación de las fotografías aéreas, por lo que se decidió realizar una tercera campaña de medición de una segunda red para el distrito Heredia.

El trabajo de campo se realizó los días 24 y 25 de febrero del 2016, utilizando mediciones de tipo estático, donde se levantaron un total de 13 puntos de control distribuidos de la siguiente forma: Parque de Urbanización La Lilliana, Hospital San Vicente de Paul (nuevo), Cementerio de Heredia, Parque Infantil de la Urbanización La Esperanza (3º etapa), Parque Infantil del Residencial Don Eloy, Parque Infantil de Urbanización La Amelia (La Puebla), puente sobre el río Pirro, Escuela de Enseñanza Especial (San Rafael de Heredia), Bajo Los Molinos, Escuela de Fátima, Universidad Estatal a Distancia (UNED, Mercedes Norte), plaza de Cubujuquí y Antiguo Hospital San Vicente de Paul.

Como se observa en el Mapa 3, los sitios seleccionados para esta red de apoyo se ubican en la periferia del distrito Heredia y con ellos se buscó aumentar la cantidad y mejorar la distribución de los puntos para la elaboración de un Modelo de Elevación Digital y para el trabajo con las fotografías aéreas.

Los datos crudos obtenidos tras este levantamiento estático se post-procesaron siguiendo los mismos pasos que se reseñaron para las otras dos redes, es decir, se utilizaron datos de estaciones de medición continua (ver tabla 7) y efemérides precisas (tabla 8) para ajustar las coordenadas finales de alta precisión.

Tabla 7:
Fechas de las mediciones y sus correspondientes días según calendario GNSS, junto con el nombre de cada archivo descargado de las estaciones de medición continúa.

Fuente: Elaboración propia a partir de datos descargados.

Mapa 3:
Distribución de los puntos de control de apoyo para la Red Heredia.

Tabla 8:
Fechas de las mediciones y sus correspondientes días según calendario GNSS, junto con el nombre de cada archivo descargado de las órbitas precisas.

Fuente: Elaboración propia a partir de datos descargados.

Para esta etapa del proyecto se solicitó al IGN el material fotográfico aéreo crudo escala 1:1000 y 1:5000, utilizado como base para la elaboración del mosaico de ortofotos disponibles en el visor del Sistema Nacional de Información Territorial (SNIT). El IGN facilitó un total de 13 fotografías aéreas escala 1:6000, acompañadas de los reportes de calibración de las cámaras utilizadas en estos vuelos fotogramétricos y la información de posicionamiento cinemático de las líneas de vuelo PRCR-CR-017 y PRCR-CR-026.

De las 13 fotografías aéreas proporcionadas por el IGN se identificaron y seleccionaron cuatro imágenes que cubrían la totalidad del área de estudio, a saber: PRCR-CR-017_166, PRCR-CR-017_168, PRCR-CR-026_055 y PRCR-CR-026_057.

Estas fotografías se georreferenciaron utilizando el software ArcGIS 10.3, específicamente, la herramienta Georreferencia, mediante la cual se logró establecer una relación de correspondencia entre las imágenes aéreas y los puntos de control generados en las anteriores etapas de este proyecto.

Para cada fotografía aérea se ubicaron entre 10 y 12 puntos de control, por lo que fue posible utilizar un polinomio de segundo orden para la asignación de las coordenadas y en promedio el RMS con el que se georreferenciaron las imágenes fue inferior a 1 metro.

Posteriormente, se delimitaron como polígono en formato shp, áreas de corte en cada fotografía que permitieran cubrir todo el ámbito de estudio y mantener la sobre posición entre ellas, con el fin de generar un mosaico que integre la cuatro imágenes. El corte y generación del mosaico se elaboraron en el programa ERDAS 2013, utilizando las herramientas Subset y Mosaic Pro del menú Raster.

La ortorectificación de las fotografías aéreas se realizó en el software ERDAS 2013, mediante el módulo IMAGINE Photogrammetry, en el cual se generaron bloques de trabajo con pares de fotografía aéreas organizados de la siguiente forma:

La creación del bloque fotogramétrico inicia con la selección del modelo geométrico que en este caso corresponde a una cámara tradicional (Frame Camera) bajo la categoría de Cámaras (Camera).

Se definió también el sistema de referencia horizontal y vertical de las imágenes aéreas, después de repasar las características de los datos con los que se cuenta para el proceso de ortorectificación, a saber:

• Datos procesados en el TBC: en el informe del ajuste de red se ofrecen varios juegos de coordenadas, las primeras denominadas de cuadrícula se indican las coordenadas mN, mE para el sistema horizontal CTRM05 (variación SIRGAS) y en vertical la elevación (en metros) que corresponde a la altura ortométrica (H). En un segundo cuadro se encuentran las coordenadas geodésicas ajustadas, en las que se visualizan coordenadas de latitud, longitud para el sistema horizontal y en vertical la altura (en metros) que se designa altura elipsoidica (h) porque depende del elipsoide de referencia, en este caso el EGM08.

• Datos de posicionamiento cinemático de las líneas de vuelo: indican las coordenadas Norte y Este para el sistema de referencia horizontal CRTM05 y la altura (en m) elipsoidica, pero asociada al elipsoide EGM96.

Es posible notar que ambos juegos de datos coinciden en el sistema de referencia horizontal, no así en el vertical. Con la intención de trabajar con una única altura elipsoidica se calculó la ondulación del geoide (N), utilizando la herramienta calculadora EGM96 (National Geospatial – Intelligence Agency, 2015).

En esa herramienta se introdujeron los valores de coordenadas geodésicas (latitud y longitud en grados, minutos y segundos), lo que devuelve un valor correspondiente a la ondulación del geoide (N) EGM96.

Para obtener la nueva altura ortométrica (H) asociada al elipsoide EGM96 es necesario restarle a la altura elipsoidica (h) la ondulación del geoide, es decir, H = h-n.3 Lo anterior permitirá que se pueda utilizar CRTM05 como sistema de referencia vertical y el EGM96 como sistema horizontal.

Tras analizar y resolver la discordancia entre datos, se procedió a definir el sistema de proyección: CRTM05, que al no estar incluido dentro de las opciones predeterminadas del programa ERDAS obligó a ingresar los siguientes parámetros:

• Tipo de proyección: Transversal Mercator

• Esferoide: WGS84

• Datum: Costa Rica 2005 (To WGS 84 1)

• Factor de escala: 0.999900

• Longitud del meridiano central: 84:00:00.0000 W

• Latitud del origen de la proyección: 0:00:00.0000 N

• Falso Este: 0.000000 metros

• Falso Norte: 500000.000000 metros

Inmediatamente, se definió el sistema vertical de los datos, que como se explicó anteriormente se homogenizó a EGM96, sin embargo, ERDAS 2013 no cuenta con ese sistema dentro de sus opciones, por lo que se seleccionó GRS80 que es un sistema de referencia horizontal muy similar.

De seguido se completó el cuadro de información específica de las tomas fotográficas, manteniéndose los valores por defecto para sistema de rotación (Omega, Phy, Kappa), las unidades angulares (grados) y la dirección del ángulo de toma (Z). Se colocó la altura promedio de vuelo y se ingresó a la edición de la cámara.

De los reportes de calibración de cada cámara se extrajo la informacion del nombre y descripción, el valor de la distancia focal y del puntos principal Xo,Yo , de las marcas fiduciales y la distorsión de los lentes; datos que se colocaron en la edición de la cámara.

Finalizada la creación del bloque, se visualizó la ventana principal del módulo de fotogrametría de ERDAS, donde se agregaron las imágenes a ortorectificar. Se seleccionaron tres imágenes de cada línea de vuelo, específicamente, la PCRC-CR-026-055, PCRC-CR-026-056 y PCRC-CR-026-055. En un segundo bloque se ortorectificaron las fotografías aéreas PCRC-CR-017-166, PCRC-CR-017-167 y PCRC-CR-017-168.

Se inició con la definición de la orientación interna, posicionando cada una de las marcas fiduciales sobre las imágenes y manteniendo el orden indicado en las fotografías. Luego se continuó con la orientación externa del bloque de imágenes, proceso realizado con la herramienta Classic Point Measurement, ventana de trabajo en la que las fotografías aéreas se trabajan en pares, es decir, se seleccionaron primeros dos fotografías aéreas continuas y luego el otro par.

Se identificaron en las imágenes los sitios en los que se había realizado al medición GNSS estática, pues para estos sitios se contaba con puntos de calidad control generados específicamente para la ortorectificación. Se colocaron puntos en el primer par, otros en el segundo par y para cada uno se anotó el valor de las coordenada X, Y y Z ajustada y se les definió como puntos de tipo FULL (certeza de la coordenada X,Y,Z) y de Control (por la calidad).

Adicionalmente, se permitió que el programa agregará automáticamente una serie de puntos en común (TIE points) en las tres imágenes, de los que se eliminaron aquellos localizados en techos o sombras, indicando que correspondían a tipo HORIZONTAL (confianza solo en las coordenadas X, Y) y manteniendo la clasificación en TIE.

Dentro del mismo módulo de definición de la orientación externa, se accedió a la Aerotriangulación, logrando la convergencia de los puntos y en la que se alcanzó un RMS de 21.0108 pixeles. Así, se completó la orientación externa de las imágenes y se observó cómo las fotografías aéreas se posicionaron de acuerdo a los puntos de control colocados, traslapándose.

Figura 8:
Ventana principal del módulo IMAGINE Photogrammetry de ERDAS, después de completar la orientación interna y externa.
Fuente: Visualización del proceso de ortorectificación en el módulo IMAGINE Photogrammetry de ERDAS.

Se avanzó hacia el proceso de generación del Modelo de Elevación Digital (MDE) que complementó la información de elevación necesaria para la ortorectificación de las fotografías aéreas. El MDE se generó con una resolución espacial de 2m, utilizando como base los 27 puntos de control generados en el levantamiento GNSS.

Finalmente, se llegó a la generación de las imágenes ortorectificadas, revisando que las tres imágenes de cada bloque se incluyeran y manteniendo como sistema de proyección de salida el CRTM05. El MDE generado para el primer bloque fue el mismo utilizado para la ortorectificación del segundo bloque.

Se obtuvieron así dos mosaicos de fotografías aéreas ortorectificadas que se integraron en un único mosaico utilizando la herramienta Mosaic Pro del menú Raster.