ORTORECTIFICACIÓN DE IMÁGENES SATELITALES GEOEYE POR RCP Y

GCP MEDIANTE PCI GEOMATICS Y ERDAS 2020

PAOLA CRISTINA MORALES MARTÍNEZ

PABLO ANDRÉS LOMBANA ZORRO

UNIVERSIDAD MILITAR NUEVA GRANADA

FACULTAD DE INGENIERÍA
ESPECIALIZACIÓN EN GEOMÁTICA
BOGOTÁ D.C.
2020
ORTORECTIFICACIÓN DE IMÁGENES SATELITALES GEOEYE POR RCP Y
GCP MEDIANTE PCI GEOMATICS Y ERDAS 2020

Autores:
PAOLA CRISTINA MORALES MARTÍNEZ 3101502
PABLO ANDRÉS LOMBANA ZORRO 3101509

Docente:
Oscar Forigua
Ingeniero Catastral y Geodesta

UNIVERSIDAD MILITAR NUEVA GRANADA

FACULTAD DE INGENIERÍA
ESPECIALIZACIÓN EN GEOMÁTICA
BOGOTÁ D.C.
2020
TABLA DE CONTENIDO

Introducción ................................................................................................................................4

Marco Teórico.............................................................................................................................5

Marco Teórico Conceptual ..........................................................................................................5

1.1. Imagen Satelital ............................................................................................................5

1.2. Irradiancia y Radiancia..................................................................................................6

1.3. Nivel Digital .................................................................................................................6

1.4. Ground Sample Distance (GSD)....................................................................................6

1.5. Resolución Espacial, Radiométrica, Espectral ...............................................................6

1.6. Orientación Interna, Externa y Absoluta ........................................................................7

Objetivos.....................................................................................................................................8

4.1. General .........................................................................................................................8

4.2. Específicos ....................................................................................................................8

Metodología ................................................................................................................................9

5.1. Área de estudio .............................................................................................................9

5.2. Materiales .....................................................................................................................9

5.2.1. Insumos .................................................................................................................9

5.2.2. Software .............................................................................................................. 10

5.3. Métodos ...................................................................................................................... 10

Resultados................................................................................................................................. 19

6.1. Generación de Ortofoto. .............................................................................................. 19

6.2. Verificación con puntos de control posicional ............................................................. 21

6.3. Calidad de los Datos ................................................................................................... 21

Conclusiones ............................................................................................................................. 24

2
Bibliografía ............................................................................................................................... 25

Anexos ...................................................................................................................................... 27

Índice de figuras

Figura 1. Área de estudio, Departamento de Cundinamarca, Colombia ........................................9

Figura 2. Ejemplo de ingreso de información para la creación de un nuevo proyecto. ................ 11
Figura 3. Imagen satelital original. ............................................................................................ 12
Figura 4. Barra, Menú de ortorectificación en PCI GEOMATICS. ............................................ 12
Figura 5. Generación de la Orthoimage con software PCI GEOMATICS. ................................. 14
Figura 6. Matriz de valores de pixeles en el canal 2 de la imagen DEM 1 (DEM_1_2). ............. 15
Figura 7. Configuración modelo geométrico .............................................................................. 16
Figura 8. Cargue de archivo *.TIL............................................................................................. 17
Figura 9. Ventana de activación para ajuste de puntos de control terrestre de manera manual. ... 17
Figura 10. Ventana Ortho Resample para creación de la Ortho Imagen...................................... 18
Figura 11. Imagen Ortorectificada mediante PCI GEOMATICS 2018 ....................................... 19
Figura 12. Imagen Ortorectificada mediante ERDAS IMAGINE 2020 ...................................... 20

Índice de tablas

Tabla 1. Insumos recopilados .................................................................................................... 10

Tabla 2. Reporte de residuales Erdas Imagine............................................................................ 18
Tabla 3. Comparación entre la imagen original y la ortofoto realizada en ERDAS IMAGINE 2020.
................................................................................................................................................. 21
Tabla 4. Análisis estadístico de resultados del proceso de ortorectificación. .............................. 21
Tabla 5. Análisis de calidad de datos en la ortoimagen procedente de PCI GEOMATICS y ERDAS
IMAGINE 2020. ....................................................................................................................... 22

3
 Introducción

Los diferentes cambios que presenta la superficie terrestre con el pasar del tiempo, promueven a
que cada vez más se busquen metodologías científicas, en donde, se permita clarificar y estudiar
plenamente la dinámica de los objetos espaciales y su posterior desarrollo en las diferentes
regiones de Colombia y el mundo. La toma de imágenes de la superficie mediante satélites, drones
o vuelos tripulados permite detectar alteraciones en el territorio bien sea de tipo natural o antrópico,
por ende, estos insumos requieren que se presenten ajustes con un alto grado de precisión en el
cual se evite o se minimice el desplazamiento de las orto-imágenes y se cumplan los estándares de
calidad exigidos por las autoridades competentes (Arevalo et al., 2004)

Por otra parte, la producción de la información temática es indispensable puesto que permite
obtener elementos espaciales representados correctamente georeferenciados. Lo anterior toma
como base fundamental que dichos datos cartográficos sean extraidos mediante la orto-imagen, la
cual, se constituye de una fotografía que representa la superficie del terreno y en donde las
deformaciones causadas por el terreno y los problemas de inclinación de la cámara han sido
correctamente ajustados (Jauregui, 2006)

En consideración, con la importancia de la generación de productos cartográficos altamente

confiables, se procede a trabajar con una imagen GeoEye con fecha de julio 18 de 2016, que
recubre los municipios de San Cayetano, Carmen de Carupa, Ubate, Sutatausa y Tausa en el
departamento de Cundinamarca, con la cual posteriormente se va a desarrollar una metodología
de generación de la orto-imagen utilizando coeficientes racionales polinomiales y puntos de
control terrestre, en el software PCI Geomatics y Erdas Imagine 2020.

4
 Marco Teórico

La fotogrametría análoga es una disciplina que permite generar modelos en 3 dimensiones

tomando como base imágenes que se encuentran en 2 dimensiones, con el objetivo de poder
obtener características de los objetos que se representan en el espacio (Jauregui, 2006). El paso de
la fotogrametría análoga a digital, no solo implica un avance en el proceso de interpretación de
imágenes si no también, conlleva la llegada de nuevas técnicas en extracción de información
geométrica, mejora gráfica de las imágenes, correspondencia automática en la identificación de
puntos de control entre otros, que esencialmente requieren métodos más sofisticados y
matemáticos que permitan el correcto procesamiento de la imagen (Gómez Lahoz, 2008), llevando
a todo un procesamiento de la imagen mediante técnicas científicas y con software especializado.

El problema que se presenta cuando se pretende trabajar con fotografías áreas sin algún tipo de
rectificación, es que estás no encajen entre ellas, este problema puede deberse a varios factores: a)
La óptica de los elementos de captura de imagen; b) la orografía provoca distorsión relacionado
con la pendiente del terreno; c) errores provocados por la posición o movimiento del satélite entre
otros, lo que obliga a realizar el proceso de ortorectificación para obtener precisión en el elementos
capturados de la imagen (Benitez, 2002). “El ortofotomosaico u ortoimagen es una presentación
fotográfica en proyección ortogonal de una zona de la superficie terrestre, en la que se ha removido
las distorsiones geométricas causadas por el sistema de captura y el relieve, en donde, todos los
elementos que están al nivel del terreno presentan la misma escala, con la misma validez de la
cartografía básica vectorial” (IGAC, 2016).

Marco Teórico Conceptual

Esta sección se enfoca en una revisión de los conceptos principales sobre las imágenes de satélite
y los aspectos propios sobre la información que se obtiene mediante un sensor.

1.1. Imagen Satelital

Una imagen de satélite es una representación gráfica de los datos reflejados en la superficie de la

5
tierra, que es posible capturar mediante un sensor óptico artificial montado en un satélite, en donde
los datos son enviados a una estación en tierra que los almacena y los procesa para generar el
producto final para el usuario (INEGI, 2015).

1.2. Irradiancia y Radiancia

La irradiancia es la magnitud que describe la potencia por unidad de superficie que relaciona todo
tipo de energía electromagnética y la radiancia, corresponde a la cantidad de energía
electromagnética que es emitida por una unidad de área particular (Oliver & Guevara, 2017).

1.3. Nivel Digital

Según los autores (Chavez & Bautista, 2005), el nivel digital corresponde al valor de tipo numérico
que no es visual y que puede ser fácilmente asimilado a una intensidad, en donde, un convertidor
digital registra esta información obtenida de sensores en las imágenes.

1.4. Ground Sample Distance (GSD)

Correspondiente con estudios realizados, el GSD es una dimensión lineal que representa el tamaño
de un pixel respecto a la huella que se toma en el suelo (Comer et al., 1998). Esta información
representa la mínima unidad espacial en la cual los objetos son claramente separables uno de otro
dentro de la imagen satelital.

1.5. Resolución Espacial, Radiométrica, Espectral

Según Labrador et al., 2012, los sensores que se instalan en los satélites para la teledetección
poseen una serie de características de las imágenes y vienen definidas por los tipos de resolución:

 Espacial: Es una medida de la distancia angular mínima que puede capturar un sensor de
la superficie terrestre.
 Radiométrica: Cantidad mínima de energía requerida que se almacena en un pixel que se

6
 denomina como nivel digital.
 Espectral: Número y anchura de las bandas espectrales que tiene un sensor y de las cuales
se pueda discriminar.

1.6. Orientación Interna, Externa y Absoluta

El autor Rosado, 2014 define a orientación interna, externa y absoluta de la siguiente manera: 1)
transformación de un sistema de coordenadas pixel al sistema de coordenadas cartesiano con un
origen previamente definido; 2) determina la posición exacta del sensor al instante en que se realizó
la captura de la imagen; 3) consiste en trasladar, nivelar, y escalar un modelo a su verdadero
posición en el espacio.

7
 Objetivos

4.1. General

Generar la ortorectificación de una imagen satelital Geoeye, a partir de dos programas especialistas
en procesos fotogramétricos, realizar el análisis y comparación de los RMS obtenidos para cada
caso.
4.2. Específicos

 Realizar la ortorectificación de una imagen satelital Geoeye por medio del software PCI
GEOMATICS, teniendo como base la captura de coordenadas de 6 puntos de control
tomados en campo.
 Realizar la ortorectificación de una imagen satelital Geoeye por medio del software
ERDAS IMAGE, teniendo como base la captura de coordenadas de 6 puntos de control
tomados en campo.
 Calcular la precisión de rectificación de una imagen satelital Geoeye a través del análisis
de los RMS y los intervalos de confianza obtenidos a partir de la captura de coordenadas
de algunos puntos de chequeo.
 Analizar y comparar los resultados obtenidos a partir de lo dispuesto en la Resolución 1392
de 2016 del IGAC, sobre los parámetros de conformidad y precisión para las ortoimagenes.

8
 Metodología

5.1. Área de estudio

El departamento de Cundinamarca está ubicado en la zona central de Colombia perteneciente a la

región Andina, rodeado por los departamentos de Boyacá, Tolima, Caldas y el Meta. La imagen
(delimitada en color morado) tiene un recubrimiento en los municipios de Ubaté, San Cayetano,
Tausa, Sutatausa y Carme de Carupa, en la Figura 1 se representa la ubicación geográfica del
departamento y el conjunto de municipios que lo integran.

Figura 1. Área de estudio, Departamento de Cundinamarca, Colombia

5.2. Materiales

5.2.1. Insumos

Las imágenes y documentos fueron suministrados por el docente y se describen a continuación:

9
 Tabla 1. Insumos recopilados
Año o
Insumo Unidad Formato Fuente Descripción
Periodo
Imagen Satélite 2016 Pixel Raster Digital Imagen de una zona de Cundinamarca de
Globle resolución espacial de 2 metros, resolución
(Geoeye) radiométrica de 16 bit sin signo, imagen
multiespectral de 4 bandas.
DEM 1 2017 Pixel Raster Intermap Imagen que contiene las elevaciones para el
recubrimiento de la zona, cuenta con 3 canales y
una resolución espacial de 9 metros, resolución
radiométrica de 32 bits
DEM 2 2017 Pixel Raster Intermap Imagen que contiene las elevaciones para el
recubrimiento de la zona, cuenta con 1 canal y una
resolución espacial de 30 metros, resolución
radiométrica de 32 bits
Trabajo Contiene los 6 puntos de control terrestre tomados
GCP 2016 - Doc
Campo para la ortorectificación de la imagen.
Cartografía Planchas
209IA2, 209IA4, 2017 - Vector IGAC Planchas a escala 1:10.000 del IGAC
209IB1y 209IB3

5.2.2. Software

Para el desarrollo de la investigación se utilizó el software PCI Geomatics 2018, mediante el

módulo de OrthoEngine y Focus. Por otra parte, también se utilizó el software Erdas Imagine 2020
con licencia estudiantil.

5.3. Métodos

El desarrollo del proceso de ortorectificación de la imagen, procedió en primera instancia a la

verificación de las propiedades de la imagen, en donde, estuvo dado por la consulta de las
características de los pixeles, resolución espacial (GSD=2 metros), resolución radiométrica (16
bit) y en general una vista previa de la información que contiene la imagen. Posteriormente, se
desarrolla el proceso de ortorectificación el cual en esencia corresponde al mismo, pero que
dependiendo del programa, varia en ciertos aspectos.

10
5.3.1. Proceso de Ortorectificación de la imagen en PCI GEOMATICS

 Creación de un nuevo proyecto de Ortorectificación:

Para la creación de un nuevo proyecto de ortorectificación, en primera instancia

utilizaremos el software GEOMATICA BANFF (2020). Desde la barra principal que se
presenta al abrir el programa, seleccionamos el ícono de “ORTHOENGINE”, luego en el
botón “FILE” seleccionamos nuevo proyecto que para el presente informe lo nombraremos
“Ortho_PCI”.

El modelo matemático empleado para la ortorrectificación de la imagen será el RFM

(RATIONAL FUNCTIONS MODEL) que extrae los coeficientes racionales directamente
de la imagen satelital (RPC. RATIONAL POLYNOMIAL COEFFICIENTS).

La proyección empleada para el proyecto de ortorectificación será UTM con datum

WGS84 en zona 18 que corresponde a Bogotá, de igual manera es la misma proyección
empleada para la captura en campo de los puntos de control. Es necesario recordar que el
tamaño de los pixeles de salida es la misma que contiene la imagen original,
correspondiente a 2 x 2 pixeles.

Figura 2. Ejemplo de ingreso de información para la creación de un nuevo proyecto.

Luego procedemos a cargar la imagen satelital a partir del parámetro “DATA SET”,
agregamos la imagen con formato de archivo .imd, el cual es el formato perteneciente al

11
 software PCI GEOMATICS (16FEB11154004-M2AS-05406576110_01_P001.IMD)

Figura 3. Imagen satelital original.

 Toma de Puntos de Control de forma manual:

En el campo desplegable de “PROCESSING STEP” seleccionamos la opción de “GCP/TP

COLLECTION” que nos permite cargar los puntos de control tomados en campo, a
continuación, se muestra el procedimiento:

Figura 4. Barra, Menú de ortorectificación en PCI GEOMATICS.

Al seleccionar el botón para ingresar puntos manuales, nos despliega una ventana dentro
de la cual ingresamos las coordenadas de cada uno de los puntos (x,y,z) y paralelamente
podemos perfeccionar la ubicación del punto, según el reporte de los puntos de control
registrado en campo.

Luego procedemos al cálculo de los residuales por cada uno de los 6 puntos de control y el
RMS obtenido al momento de la ubicación de los mismos. Para lo anterior, damos check
en la casilla de “COMPUTE MODEL”, el cual nos calcula los residuales individuales. Al
obtener estos datos se pueden observar los puntos que presentan residuales altos y de esta

12
 forma permite corregir su ubicación con el fin de minimizar el error.

Por último podemos analizar los RMS resultantes al momento de comparar los valores de
0 a 2 del campo de RPC, permitiendo ajustar de mejor manera el error medio cuadrático
final.

 Selección de la imagen DEM (Modelo Digital de Elevación):

Dentro de los insumos encontramos dos tipos de imágenes DEM, cada una cuenta con
características distintas. En primera instancia, la imagen DEM_1 cuenta con 3 canales que
describen características distintas, en cuanto a tamaño de pixel, datos de Nivel Digital y
completitud de información, así mismo, la imagen DEM_2 cuenta con un solo canal. A
continuación, describiremos las diferentes características para cada una de las imágenes:

- Canal 1, imagen DEM_1: El primer canal cuenta con información de pixeles en

números enteros positivos y con un tamaño de pixel de 9 metros.

- Canal 2, imagen DEM_1_2: El segundo canal cuenta con valores decimales en su

matriz de niveles digitales, lo cual permite tener mayor exactitud a la hora de realizar
la rectificación de la imagen satelital. De igual manera, el tamaño de sus pixeles es de
9x9 metros y por último cuenta con gran cantidad de información almacenada como
ND.

- Canal 3, imagen DEM_1_3: El tercer canal contiene en su matriz valores enteros,

además, no presenta gran cantidad de información geográfica, en su mayoría valores
nulos, lo que generaría que grandes campos de información digital de la GeoEye no
pueda ser corregida.

 Creación de la Orthoimage:

En el desplegable “PROCESSING STEP” seleccionamos “ORTHO GENERATION”, el

cual nos permite generar la ortorectificación de la imagen satelital a partir de los puntos de
control ingresados anteriormente junto con la imagen DEM seleccionada.

13
 Figura 5. Generación de la Orthoimage con software PCI GEOMATICS.

En la ventana que se despliega ingresamos la imagen original a procesar y luego

procedemos a verificar los datos en cada uno de los campos. El proceso de
georeferenciación de imágenes de satélite requiere introducir algunos parámetros
necesarios en la corrección. Independientemente del modelo de corrección utilizado existen
dos parámetros muy importantes: el tamaño de píxel en la imagen corregida y el método
de interpolación espacial utilizado durante el remuestreo (habitualmente el vecino más
próximo, la interpolación bilineal o la interpolación bicúbica).

En la sección de Ortho Image procedemos a seleccionar la carpeta donde se va a almacenar

la imagen ortorectificada y procedemos a renombrarla. De igual forma, verificamos que el
formato de salida de la imagen sea archivo .tiff y que el NoData sea 0.

En el campo de opciones de la orto generación, cargamos la imagen DEM que se escogió

con anterioridad (DEM_1_2) la cual presenta las mejores características permitiendo
obtener un intervalo de confianza conforme y un RMS aceptable.

14
 Figura 6. Matriz de valores de pixeles en el canal 2 de la imagen DEM 1 (DEM_1_2).

En el campo “opciones de procesamiento” observamos los diferentes métodos de

interpolación para el remuestreo. En este caso escogemos bilineal.

Luego de verificar todos los campos necesarios, procedemos a generar la imagen

ortorectificada.

 Toma de Puntos de chequeo:

Para la evaluación de precisión de la imagen ortorectificada, en el presente trabajo fue

necesario la captura de coordenadas de 20 puntos de chequeo que han sido tomados en
campo previamente.

Para lo anterior, empleamos como insumos los datos de las cartas del IGAC (Instituto
Geografico Agustin Codazzi) pertenecientes a shape files que contienen información sobre
vías, cuencas, ríos y demás objetos pertenecientes a la zona del departamento de
Cundinamarca y que hacen parte de la zona captada en la imagen satelital.

Luego de la captura de las coordenadas planas (x,y,z) de los puntos de chequeo

directamente de las capas del IGAC, procedemos a capturar las coordenadas planas de los

15
 mismo puntos pero en la imagen ortorectificada que ha sido cargada al software ARGIS
previamente. De esta forma podemos verificar las diferencias entre coordenadas para los
mismos puntos y a partir de esto, verificar la precisión horizontal obtenida a partir de los
RMS (Error Medio Cuadratico) para el 90% y el 95%.

5.3.2. Proceso de ortorectificación de la imagen en ERDAS IMAGINE

 Creación del proyecto para la ortorectificación:

Primero es necesario ejecutar Erdas Imagine y luego para crear el proyecto se dirige a la
pestaña File, luego seleccionando la opción new y photogrammetric Project.
Posteriormente, se debe nombrar el archivo con extensión *.blk con el nombre de
“Ortho_Erdas.blk”.

Se deben configurar los parámetros del modelo para la ortorectificación, sobre el cual se
selecciona un modelo que utilice los coeficientes racionales polinomiales y donde el
modelo geométrico se asocie al Geoeye RCP. Luego, se establece la proyección de salida,
que como ya se mencionó anteriormente corresponde a UTM con datum WGS84 en zona
18.

Figura 7. Configuración modelo geométrico

Una vez realizado el proceso anterior, ERDAS crea un proyecto asociado con dos visores

16
 en dos dimensiones, el cual en la parte izquierda en la tabla de contenidos asocia los grupos
de datos que se mantienen durante el proyecto (Images, Orthos, Terrains, Points).

Consecuentemente, procedemos a cargar la imagen de satélite sobre la opción de images y

se busca el espacio de trabajo en donde se encuentra ubicado la imagen, luego se carga al
archivo de cabecera con extensión *.TIL.

Figura 8. Cargue de archivo *.TIL

 Adición de los puntos de control de terreno:

Para cargar los puntos de control manualmente, se debe ir la barra de herramientas y

seleccionar la opción Point Measurement y luego Classic Point Measurement

Figura 9. Ventana de activación para ajuste de puntos de control terrestre de manera

manual.

Una vez se despliega la pantalla se procede a cargar los puntos de control terrestre
provenientes del archivo llamado “GCPs_Informe.doc”, esto se realiza de manera manual y
se refina verificando el punto en la imagen y la fotografía del punto en campo. Los puntos
se cargan en su totalidad y deben ser de tipo Full y Control en las opciones que se habilitan
para cada punto de control, por último, se procede a realizar la verificación de los residuales.

Para la verificación del RSM se selecciona la opción Rational Function Refinement en la

cual se deben configurar el número de iteraciones a un total de 10 y el order polinomial
activarlo en 1, para que luego nos genere el reporte de los residuales tal y como se observa
en la siguiente tabla.

17
 Control point image residuals

The image ID = 1
Point ID Vx Vy
1 -0.667 -1.035
2 0.203 2.058
3 -0.887 0.126
4 0.826 -0.406
5 -0.357 -0.010
6 0.881 -0.732
meanx=-0.000, meany=-0.000
rmsex=0.690, rmsey=1.002

Tabla 2. Reporte de residuales Erdas Imagine

Como se relacionó en el proceso de ortorectificación de PCI Geomatics, para la selección

del DEM, se deja definido el archivo correspondiente a DEM-1 canal 2. Por último, se
procede a realizar la ortofoto.

Se dirige a la pestaña Ortho Resample, allí se debe configurar la ruta del archivo de salida
con el nombre de Ortho_Erdas.img, adicionalmente, se configura el DTM Source con el
DEM seleccionado en el paso anterior. Ya para finalizar se configura el tamaño de celda de
2 en X y en Y.

Figura 10. Ventana Ortho Resample para creación de la Ortho Imagen.

18
Resultados

6.1. Generación de Ortofoto.

En el presente capitulo se realizará un análisis de los resultados arrojados en cada uno de los
programas utilizados para la ortorectificación de la imagen satelital, a su vez se hará la
comparación de los RMS arrojados en el momento del cálculo de la precisión a partir de las
coordenadas de los puntos de chequeo y por último, se analizarán los valores obtenidos en los
intervalos de confianza de 90 y 95%.

A continuación, analizamos los resultados obtenidos de las imágenes ortorectificadas para cada
software

Figura 11. Imagen Ortorectificada mediante PCI GEOMATICS 2018

19
 Figura 12. Imagen Ortorectificada mediante ERDAS IMAGINE 2020

Las Figura 12 y 11 muestran el resultado final de la ortorectificación de la imagen Geoeye en

donde se puede observar como el relieve sobresale para la zona más montañosa, además, se
elimina la zona de background que tenía la imagen original. Las imágenes se presenten en
composición de color natural RGB con un realce mediante el método de histograma ecualizado.

El siguiente cuadro muestra la diferencia que se presenta en algunas zonas de la imagen

delimitadas por el área de interés. La figura de la izquierda se encuentra en falso color para
efectos de la visualización y la de la derecha en color natural ortorectificada respecto de la
imagen original:
Imagen Original Imagen Ortorectificada Análisis de resultado
Con el AOI previamente
definido se puede observar que
tomando la misma área la
imagen de la izquierda se
encuentra bastante desplazada
respecto de la ortofoto, por otra
parte, se alcanzan a notar las
variaciones que se presentan en
la pendiente del terreo ya que
algunas zonas se ven más
planas en la imagen original,
pero en la ortofoto si se
delimitan bien las zonas
montañosas, lo que podría

20
 deducir que ya se estableció
correctamente la corrección
geométrica.

Tabla 3. Comparación entre la imagen original y la ortofoto realizada en ERDAS IMAGINE

2020.

6.2. Verificación con puntos de control posicional

Los puntos de verificación posicional o de chequeo (20 en total, ver anexo 1) arrojaron los
siguientes resultados:

PARÁMETROS Erdas Imagine PCI Geomatics

RESULTADOS Valor Unidad Valor Unidad
NUMERO DE PUNTOS (n) 20 puntos 20 Puntos
MEDIA 2,348136444 metros 1,986113995 metros
ERROR MEDIO CUADRÁTICO 2,944227094 metros 2,488302433 metros
DESV. STANDARD 1,822297112 1,537941274
ERROR STANDARD 0,407478022 0,343894123
LIMITE ERROR GRUESO 7,815027778 metros 6,599937817 metros
Z 1,28 1,28
INT. DE CONFIANZA PARA UN 90% 2,869708311 metros 2,426298473 metros
INT. DE CONFIANZA PARA UN 95% 5,075861903 metros 4,250676083 metros
Tabla 4. Análisis estadístico de resultados del proceso de ortorectificación.

La Tabla 4 contiene el resumen de los datos estadísticos obtenidos. A partir de esto, podemos
observar que la imagen ortorectificada con mejor precisión es la obtenida con el programa PCI
GEOMATICS, donde su RMS es de 2.48 metros de acuerdo a la precisión del satelite cercano
a 2 metros. Esto sumado a el resultado que arroja el intervalo de confianza al 95 % (4,25), indica
que los objetos pueden estar desplazados aproximadamente 4 metros de su posición real, por
otra parte, es importante aclarar que ambas imágenes pueden ser utilizadas para interpretación
bien sea de cobertura, cartografía temática o cartografía convencional.

6.3. Calidad de los Datos

Según el IGAC en su documento ESPECIFICACIONES TÉCNICAS CARTOGRAFÍA

BÁSICA del año 2016, “Los elementos y subelementos de calidad y las medidas
correspondientes del ortofotomosaico/ortoimagen están definidos detalladamente en la ficha

21
 del perfil de especificación técnica para cada escala conforme a la norma técnica colombiana
NTC 5043”, este documento define los parámetros que deben ser utilizados para la validación
de la ortofoto generada. En la siguiente tabla se muestra los resultados del proceso de medición
de la calidad de los datos con base a la metodología del IGAC:

Software
Utilizado PCI Geomatics y Erdas Imagine 2020
Nombre Archivo Ortho PCI.tif y Ortho_ERDAS.img
Fecha 21/02/2020
Valor Nivel de Producto
Parámetro Resultado Observación
exigido conformidad conforme
La imagen recubre una
zona en específico, pero
Área
Comisión <= 3% N/A Conforme actualmente no se tiene
omitida
delimitado un área de
proyecto
Totalidad
La imagen recubre una
zona en específico, pero
Área
Omisión <= 20% N/A Conforme actualmente no se tiene
Omisión
delimitado un área de
proyecto
Este parametro no se
Consistencia Consistencia mide ya que es para el
Empalme <= 5% N/A N/A
Lógica Topológica caso en que se tengan
más imágenes
Resolución Resolución espacial de
Boleano VERDADERO Conforme
Espacial 2 metros
Resolución
Consistencia Boleano VERDADERO Conforme 4 bandas
Espectral
Conceptual Resolución
Boleano VERDADERO Conforme 16 bits
Radiométrica
De Formato Boleano VERDADERO Conforme
Exactitud 4.25 m y
Real 0.6 m Conforme
Exactitud de Absoluta 5.08 m
Posición Exactitud
Boleano VERDADERO Conforme
Temática
Tabla 5. Análisis de calidad de datos en la ortoimagen procedente de PCI GEOMATICS y
ERDAS IMAGINE 2020.
Los resultados de calidad de datos muestran que los parámetros de exactitud en posición y
consistencia conceptual arrojan valores de conformidad como es el caso de la exactitud absoluta
que da un valor de 4.25 m y 5.08 m. Lo anterior permite establecer un tamaño de GSD de 50
centímetros por cuanto los resultados del RMS son menores a 3 m y su exactitud horizontal
menor a 5.88 para las dos imágenes ortorectificadas.

De acuerdo a la tabla 25 del documento de “ESPECIFICACIONES TÉCNICAS PARA

CARTOGRAFÍA BÁSICA” del IGAC, donde condensa los lineamientos para la exactitud

22
absoluta, podemos decir que ya que el error medio cuadrático arrojado en los dos ortoimagenes
(2.48 y 2.94) es inferior al establecido para su resolución, los resultados cumplen con las
medidas de exactitud posicional.

Los valores calculados en el intervalo de confianza de 95% (4.25 m y 5.08m) nos permiten
ingerir que el tamaño del GSD podría aplicar para el de 50 centímetros, a su vez y de acuerdo
a la tabla de resolución espacial para la GSD, podría generar cartografía para escalas de
1:10.000 y 1:25.000.

23
Conclusiones

 El proceso de ortorectificación de las imágenes de satélite mediante coeficientes racionales

polinomiales y puntos de control terrestre, es una buena forma de corregir la imagen para
obtener correctamente la orientación absoluta, ya que, con puntos correctamente tomados
de campo y un buen proceso de rectificación en oficina es posible generar productos
cartográficos altamente confiables.
 El uso de parámetros estadísticos como el RMS permite validar que también se está
ejecutando el proceso de validación de los puntos de control, para que en el proceso de la
ortorectificación de la imagen, se puedan corregir dichos errores que puedan afectar el
producto final y que también se puedan cumplir los requerimientos de la autoridad estatal
competente.
 Una correcta especificación del DEM, en la cual se tienen valores flotantes para los niveles
digitales, una resolución radiométrica adecuada y una proyección acorde, permiten que el
desarrollo de la metodología y en general el producto cartográfico pueda cumplir con los
estándares requeridos por la autoridad estatal, ya que esto es una base para cualquier
proyecto fotogramétrico que se pretenda desarrollar.
 Como resultado de todo el proceso de ortorectificación, disponemos de imágenes que
facilitarían la elaboración de cartografía fidedigna. Lo anterior, evitando trabajos
cartográficos poco veraces que presenten información desligada a la realidad del terreno,
eliminando deformaciones del relieve y perfeccionando la posición de los objetos.
 La ortocorrección de una imagen Geoeye realizada en el presente informe presenta una
buena calidad de los datos, con la cual se pueden generar productos cartográficos a escalas
de 1:10.000 desde donde la forma y posición de los objetos no presentarán deformaciones.
 Es recomendable emplear DEM que presenten tamaños de pixeles acordes y que la
información almacenada en sus celdas contenga valores decimales que permitan minimizar
el error.

24
 Bibliografía

Arevalo, V., Gonzales - Jimenez, J., & Ambrosio, G. (2004). (PDF) Corrección Geométrica de

Imágenes de Satélite de Distinta Resolución (in spanish).

https://www.researchgate.net/publication/266478543_Correccion_Geometrica_de_Image

nes_de_Satelite_de_Distinta_Resolucion_in_spanish

Benitez, J. C. (2002). ORTORECTIFICACIÓN DE FOTOGRAFÍAS ÁREAS UTILIZANDO EL

SOFTWARE ERDAS IMAGINE. http://www.cartografia.cl/download/manuales/erdas.pdf

Chavez, O. E. R., & Bautista, H. A. A. (2005). MANUAL PARA EL MANEJO Y

PROCESAMIENTO DE IMÁGENES SATELITALES OBTENIDAS DEL SENSOR

REMOTO MODIS DE LA NASA, APLICADO EN ESTUDIOS DE INGENIERIA CIVIL.

230.

Comer, R., Kinn, G., Light, D., & Mondello, C. (1998). Talking Digital.

http://www.asprs.org/a/publications/pers/98journal/december/12-98-digital.pdf

Gómez Lahoz, J. (2008). INTRODUCCIÓN A LA FOTOGRAMETRÍA DIGITAL: EL MÉTODO

GENERAL DE LA FOTOGRAMETRÍA DIGITAL.

http://ocw.usal.es/eduCommons/ensenanzas-tecnicas/procesamiento-avanzado-de-

imagenes-digitales/contenidos/Tema4.pdf

IGAC. (2016). Especificaciones técnicas cartografía básica.

https://www.igac.gov.co/sites/igac.gov.co/files/especificacionesv1.pdf#overlay-

context=es/contenido/areas-estrategicas/especificaciones-tecnicas-para-la-generacion-de-

cartografia-basica

INEGI. (2015). Imágenes de Satélite Elementos de la Percepción Remota.

http://www3.inegi.org.mx/contenidos/temas/mapas/imagenes/imgsatelite/metadatos/elem

25
 _per_rem.pdf

Jauregui, L. (2006). Capítulo 1 Introducción a la Fotogrametría. En Introducción a la

Fotogrametría.

http://webdelprofesor.ula.ve/ingenieria/iluis/publicaciones/Fotogrametr%EDa/CAPITUL

O1.pdf

Labrador, G. M., Évora, J. A., & Arbelo, M. (2012). Satelites de Teledetección para la gestión del

territorio. http://www.gmrcanarias.com/wp-

content/uploads/2016/01/20_catalogo_satelites_es.pdf

Moreno, R. M., & Vayá, E. V. (2000). Técnicas econométricas para el tratamiento de datos

espaciales: La econometría espacial. Edicions Universitat Barcelona.

Oliver, A., & Guevara, A. B. (2017). Irradiancia y Radiancia. http://leias.fa.unam.mx/wp-

content/uploads/2018/07/180515_Practica15_LES.pdf

Rosado, E. Q. (2014). Introducción a la Fotogrametría y Cartografía aplicadas a la Ingeniería

Civil. 139.

26
 Anexos

Anexo 1 Informe control de Precisión

MÉTODO No. E-Imagen N-Imagen E-Vector N-Vector E N E² N² Xi Xi² (Xi-Xm)²

1 617919,05 579375,49 617920,14 579375,89 -1,091237 -0,399244 1,1907982 0,1593958 1,1619785 1,350194 0,6791994

2 622888,7 579878,74 622890,83 579877,19 -2,129007 1,548756 4,5326708 2,3986451 2,6327392 6,931316 0,4181242

3 622782,51 579742,46 622781,3 579741,31 1,210451 1,145104 1,4651916 1,3112632 1,6662697 2,7764548 0,1023004

4 622883,81 580898,99 622885,22 580897 -1,406201 1,990706 1,9774013 3,9629104 2,4372755 5,9403116 0,2035467

5 622840,84 581062,41 622841,25 581060,82 -0,410718 1,591115 0,1686893 2,5316469 1,64327 2,7003362 0,117542

6 623048,68 580804,01 623052,41 580802,47 -3,724759 1,538625 13,87383 2,3673669 4,0300368 16,241196 4,1776204

7 620738,92 582054,26 620738,93 582048,51 -0,008909 5,75407 7,937E-05 33,109322 5,7540769 33,109401 14,197544

8 619071,69 581983,37 619071,29 581983,03 0,39997 0,347787 0,159976 0,1209558 0,53003 0,2809318 2,1201806

9 621843,51 582045,84 621843,88 582046,02 -0,37654 -0,183127 0,1417824 0,0335355 0,4187098 0,1753179 2,456756

PCI 10 622420,46 582088,26 622423,41 582087,83 -2,94554 0,432 8,6762059 0,186624 2,9770505 8,8628299 0,9819552
GEOMATICS
11 624479,69 581631,77 624479,44 581631,07 0,247261 0,69896 0,061138 0,4885451 0,7414062 0,5496831 1,5492976

12 623612,23 582216,15 623612,92 582212,31 -0,69858 3,846409 0,488014 14,794862 3,9093319 15,282876 3,6987673

13 619632,51 582417,02 619631,96 582417,52 0,553965 -0,499476 0,3068772 0,2494763 0,7458911 0,5563535 1,5381529

14 620136,79 581699,02 620137,21 581697,89 -0,415775 1,13626 0,1728689 1,2910868 1,2099403 1,4639556 0,6024455

15 621148,55 578553,14 621149,34 578553,15 -0,792972 -0,011938 0,6288046 0,0001425 0,7930619 0,6289471 1,4233734

16 621441,65 578936,92 621442,6 578937,55 -0,955792 -0,63118 0,9135383 0,3983882 1,1453936 1,3119265 0,7068108

17 617734,82 582004,65 617734,95 582004,68 -0,13563 -0,025573 0,0183955 0,000654 0,1380198 0,0190495 3,415452

18 619003,98 580441,28 619006,16 580438,79 -2,177654 2,486434 4,7421769 6,182354 3,3052278 10,924531 1,7400613

19 618271,4 579535,86 618272,1 579535,96 -0,703967 -0,103076 0,4955695 0,0106247 0,7114733 0,5061942 1,624709

20 618476,99 579924,75 618478,56 579921,32 -1,567502 3,429885 2,4570625 11,764111 3,7710971 14,221174 3,1861648

1 617917,53 579375,64 617920,14 579375,89 -2,606237 -0,254244 6,7924713 0,06464 2,6186087 6,8571113 0,0731552

2 622892,16 579878,94 622890,83 579877,19 1,332993 1,750756 1,7768703 3,0651466 2,2004583 4,8420169 0,0218088

3 622782,57 579742,49 622781,3 579741,31 1,265451 1,176104 1,6013662 1,3832206 1,7275957 2,9845869 0,3850708

4 622883,72 580899,34 622885,22 580897 -1,500201 2,347706 2,250603 5,5117235 2,7860952 7,7623265 0,1918079

5 622839,97 581063,5 622841,25 581060,82 -1,276718 2,678115 1,6300089 7,1723 2,9668685 8,8023088 0,3828294

6 623047,69 580803,16 623052,41 580802,47 -4,717759 0,689625 22,25725 0,4755826 4,767896 22,732833 5,8552365

7 620738,38 582055,41 620738,93 582048,51 -0,548909 6,89907 0,3013011 47,597167 6,9208719 47,898468 20,90991

ERDAS 8 619070,34 581984,26 619071,29 581983,03 -0,94603 1,233787 0,8949728 1,5222304 1,5547357 2,4172031 0,6294847
IMAGINE
9 621843,86 582045,93 621843,88 582046,02 -0,02054 -0,089127 0,0004219 0,0079436 0,0914632 0,0083655 5,0925742

10 622419,95 582087,85 622423,41 582087,83 -3,45654 0,016 11,947669 0,000256 3,456577 11,947925 1,2286405

11 624480,19 581631,94 624479,44 581631,07 0,751261 0,87296 0,5643931 0,7620592 1,1517171 1,3264523 1,4314193

12 623609,84 582215,81 623612,92 582212,31 -3,08858 3,501409 9,5393264 12,259865 4,6689604 21,799191 5,3862239

13 619631,86 582417,52 619631,96 582417,52 -0,099035 0,004524 0,0098079 2,047E-05 0,0991383 0,0098284 5,0579928

14 620136 581699,79 620137,21 581697,89 -1,211775 1,90326 1,4683987 3,6223986 2,2562795 5,0907973 0,0084377

15 621149,34 578552,98 621149,34 578553,15 0,006028 -0,173938 3,634E-05 0,0302544 0,1740424 0,0302908 4,7266848

16 621441,71 578936,49 621442,6 578937,55 -0,887792 -1,05918 0,7881746 1,1218623 1,3820408 1,9100369 0,9333407

27
17 617735,32 582004,85 617734,95 582004,68 0,37337 0,168427 0,1394052 0,0283677 0,4096008 0,1677728 3,7579205

18 619003,78 580441,61 619006,16 580438,79 -2,374654 2,820434 5,6389816 7,9548479 3,6869811 13,59383 1,792505

19 618271,67 579535,88 618272,1 579535,96 -0,429967 -0,082076 0,1848716 0,0067365 0,4377306 0,1916081 3,6496504

20 618476,52 579924,28 618478,56 579921,32 -2,046502 2,967885 4,1881704 8,8083414 3,6050675 12,996512 1,5798757

28