T Espe 040400
T Espe 040400
T Espe 040400
CONSTRUCCIÓN
SANGOLQUÍ
2018
i
AUTORÍA DE RESPONSABILIDAD
iii
DEDICATORIA
Esta tesis se la dedico a mis abuelitos María Natividad, Manuel e Isabel, quienes en vida me
dieron los consejos y la fuerza necesaria para continuar ante toda adversidad.
A mi padre René y madre Susana, por su apoyo incondicional y sus incontables muestras de
afecto a lo largo del camino de culminar una carrera; y a mis hermanos Santiago y Jenny por ser
AGRADECIMIENTO
director de esta investigación Ing. Eduardo Kirby, por sus largas jornadas de enseñanza y ser un
A mis grandes amigos Oscar y Andrés, que me acompañaron a sobrellevar el estrés con su
A mis compañeras Sisa, Damariz y Erika por el apoyo brindado en las salidas de campo, en la
A mis diferentes profesores en el transcurso de la carrera que supieron guiarme con sus
por mantenerme firme, constante y darme la fuerza todos los días por seguir dando lo mejor de mí.
vi
ÍNDICE DE CONTENIDO
Dedicatoria ................................................................................................................................. iv
Agradecimiento ........................................................................................................................... v
Índice de contenido..................................................................................................................... vi
Abstract.................................................................................................................................... xxix
CAPITULO I ............................................................................................................................... 1
1. GENERALIDADES ................................................................................................................ 1
1.1. Introducción.......................................................................................................................... 1
1.6. Objetivos............................................................................................................................... 7
CAPÍTULO II.............................................................................................................................. 9
2.6.1. ASPRS “Positional Accuracy Standards for Digital Geospatial Data” ........................... 67
x
3. METODOLOGÍA.................................................................................................................. 82
3.2.5. Recopilación de monografías y bloque fotogramétrico con puntos Von Gruber. ........... 91
xi
3.3.3. Procesamiento.................................................................................................................. 96
3.4.1. Generación de ortofotos sin matriz de pesos y 13 GCP en LPS ................................... 100
3.4.1.7. Generación de los MDS sin matriz de pesos y 13 GCP ............................................. 107
3.4.2. Generación de ortofotos con matriz de pesos y 13 GCP en LPS .................................. 113
3.4.2.2. Generación de los MDS con matriz de pesos y 13GCP ............................................. 115
3.4.3. Generación de ortofotos sin matriz de pesos y 19 GCP en LPS ................................... 115
3.4.3.2. Configuración de la aerotriangulación sin matriz de pesos en LPS y 19GCP ........... 117
3.4.3.3. Generación de los MDS sin matriz de pesos y 19 GCP ............................................. 117
3.4.4. Generación de ortofotos con matriz de pesos y 19 GCP en LPS .................................. 118
3.4.4.1. Configuración de la aerotriangulación con matriz de pesos en LPS y 19GCP .......... 118
3.4.4.2. Generación de los MDS con matriz de pesos y 19 GCP ............................................ 119
3.5.1. Generación de ortofotos con ajuste en bloque y 13 GCP en Orima .............................. 120
3.5.1.7. Generación de los MDS con ajuste en bloque y 13 GCP ........................................... 133
3.5.2. Generación de ortofotos con ajuste de líneas de vuelo y 13 GCP en Orima ................. 134
3.5.2.1. Configuración de la aerotriangulación con ajuste de líneas de vuelo en Orima. ....... 134
3.5.2.2. Generación de los MDS con ajuste de líneas de vuelo y con 13 GCP ....................... 136
3.5.2.3. Generación de Ortofotos con ajuste de líneas de vuelo y con 13 GCP ...................... 136
3.5.2.4. Generación de Mosaicos con ajuste de líneas de vuelo y con 13 GCP ...................... 137
3.5.3. Generación de ortofotos con ajuste en bloque y 19 GCP en Orima .............................. 137
3.5.3.2. Configuración de la aerotriangulación con ajuste de bloque en Orima y 19GCP ...... 138
3.5.3.3. Generación de los MDS con ajuste en bloque y 19 GCP ........................................... 138
3.5.4. Generación de ortofotos con ajuste en líneas de vuelo y 19 GCP en Orima ................. 139
3.5.4.2. Generación de los MDS con ajuste en líneas de vuelo y 19 GCP .............................. 140
3.5.4.3. Generación de Ortofotos con ajuste en líneas de vuelo y 19 GCP ............................. 140
xiv
3.5.4.4. Generación de Mosaicos con ajuste en líneas de vuelo y 19 GCP ............................. 140
3.6.5. Generación del mosaico de MDS del ajuste en Z/I Imagine ......................................... 148
4. Resultados............................................................................................................................ 162
4.2. Resultados de las distribución de los GCP y puntos de exactitud posicional .................. 164
4.3.1.1. Resultado de la aerotriangulación sin matriz de pesos y 13GCP en LPS .................. 169
4.3.1.2. Resultado de la aerotriangulación con matriz de pesos y 13 GCP en LPS ................ 170
4.3.1.3. Resultado de la aerotriangulación sin matriz de pesos y 19GCP en LPS .................. 171
4.3.1.4. Resultado de la aerotriangulación con matriz de pesos y 19GCP en LPS ................. 172
4.4.1. Resultado del mosaico MDS sin matriz de pesos y 13 GCP en LPS ............................ 179
xvi
4.4.2. Resultado del mosaico MDS con matriz de pesos y 13 GCP en LPS ........................... 179
4.4.3. Resultado del mosaico MDS sin matriz de pesos y 19 GCP en LPS ............................ 180
4.4.4. Resultado del mosaico MDS con matriz de pesos y 19 GCP en LPS ........................... 181
4.4.5. Resultado del mosaico MDS con ajuste en bloque y 13 GCP en Orima....................... 182
4.4.6. Resultado del mosaico MDS con ajuste en líneas de vuelo y 13 GCP en Orima.......... 183
4.4.7. Resultado del mosaico MDS con ajuste en bloque y 19 GCP en Orima....................... 184
4.4.8. Resultado del mosaico MDS con ajuste en líneas de vuelo y 19GCP en Orima........... 185
4.4.9. Resultado del mosaico MDS con ajuste en bloque y 13 GCP en Z/I Imagine .............. 186
4.5.1. Resultado del ortomosaico sin matriz de pesos y 13GCP en LPS ................................ 187
4.5.2. Resultado del ortomosaico con matriz de pesos y 13 GCP en LPS .............................. 188
4.5.3. Resultado del ortomosaico sin matriz de pesos y 19 GCP en LPS ............................... 189
4.5.4. Resultado del ortomosaico con matriz de pesos y 19 GCP en LPS .............................. 189
4.5.5. Resultado del ortomosaico con ajuste en bloque y 13 GCP en Orima .......................... 190
4.5.6. Resultado del ortomosaico con ajuste de líneas de vuelo y 13 GCP en Orima ............. 191
4.5.7. Resultado del ortomosaico con ajuste en bloque y 19 GCP en Orima .......................... 191
4.5.8. Resultado del ortomosaico con ajuste en líneas de vuelo y 19 GCP en Orima ............. 192
4.5.9. Resultado del ortomosaico con ajuste en bloque y 13 GCP en Z/I Imagine ................. 193
4.6.4.1. Resultados del error planimétrico del ortomosaico E_X_i_VII_X ............................ 200
4.6.4.2. Resultados del error planimétrico del ortomosaico E_x_ii_vii_x .............................. 201
4.6.4.3. Resultados del error planimétrico del ortomosaico E_Y_i_VII_X ............................ 202
4.6.4.4. Resultados del error planimétrico del ortomosaico E_y_ii_VII_X ............................ 203
4.6.4.5. Resultados del error planimétrico del ortomosaico O_X_iii_VII_X ......................... 203
4.6.4.7. Resultados del error planimétrico del ortomosaico O_Y_iii_VII_X ......................... 205
4.6.4.9. Resultados del error planimétrico del ortomosaico en Z/I Imagine ........................... 206
ÍNDICE DE TABLAS
Tabla 11. Estándares de precisión vertical para datos de elevación digital .......................... 75
Tabla 12. Precisión horizontal/ejemplos de calidad para datos planimétricos digitales ....... 76
Tabla 14. Relación del GSD y escala de la foto para imágenes digitales............................... 86
Tabla 15. Relación entre escala de la foto y escala del mapa derivado. ................................ 87
Tabla 20. Disposición de las plantillas con los puntos de von Gruber. ................................ 124
xix
Tabla 21. Identificadores de los GCP en el ajuste Z/I Imagine. ........................................... 145
Tabla 25. Elenco de coordenadas de las monografías de Barragán (2015) ......................... 163
Tabla 35. Nivel de confianza de los ortomosaicos con la normativa IGM. .......................... 209
Tabla 36. Exactitud de los ortomosaicos con la normativa ASPRS. ..................................... 210
xx
ÍNDICE DE FIGURAS
Figura 11. Distorición del pixel por los parámetros de orientación. ........................................ 17
Figura 38. Diferencias entre DTM (en español MDT) y DSM (en español MDS).................. 50
Figura 99. Ajuste preliminar con 13GCP y puntos Von Gruber. ........................................... 130
Figura 100. Generación de los tie points para todo el bloque. ............................................... 131
Figura 103. Generación de los tie points por línea de vuelo. ................................................. 135
Figura 118. Puntos de exactitud posicional para productos con 13GCP................................ 155
Figura 119. Puntos de exactitud posicional para productos con 19GCP................................ 155
Figura 120. Puntos de exactitud posicional para productos del ajuste en Z/I Imagine .......... 156
Figura.121 Extracción de valores del MDS respecto de los puntos de control ...................... 157
Figura 125. Bloque fotogramétrico con los puntos de Von Gruber ....................................... 164
Figura 130. Reporte A. sin matriz de pesos y 13GCP en LPS ............................................... 169
xxvi
Figura 132. Reporte A. con matriz de pesos y 13GCP en LPS .............................................. 171
Figura 133. Reporte A. sin matriz de pesos y 19GCP en LPS ............................................... 171
Figura 134. Reporte A. con matriz de pesos y 19GCP en LPS .............................................. 172
Figura 141. Mosaico de MDS sin matriz de pesos y 13GCP en LPS .................................... 179
Figura 142. Mosaico MDS con matriz de pesos y 13GCP en LPS ........................................ 180
Figura 143. Mosaico MDS sin matriz de pesos y 19 GCP en LPS ........................................ 181
Figura 144. Mosaico MDS con matriz de pesos y 19 GCP en LPS ....................................... 182
Figura 145. Mosaico MDS con ajuste en bloque y 13 GCP en Orima................................... 183
Figura 146. Mosaico MDS con ajuste en líneas de vuelo y 13 GCP en Orima. .................... 184
Figura 147. Mosaico MDS con ajuste en bloque y 19 GCP en Orima................................... 185
Figura 148. Mosaico MDS con ajuste en líneas de vuelo y 19 GCP en Orima. .................... 186
Figura 149. Mosaico MDS con 13GCP y ajuste bloque en Z/I Imagine................................ 187
Figura 150. Ortomosaico sin matriz de pesos y 13 GCP en LPS ........................................... 188
Figura 152. Ortomosaico sin matriz de pesos y 19 GCP en LPS ........................................... 189
xxvii
Figura 154. Ortomosaico con ajuste en bloque y 13 GCP en Orima ..................................... 190
Figura 155. Ortomosaico con ajuste de líneas de vuelo y 13 GCP en Orima. ....................... 191
Figura 156. Ortomosaico con ajuste en bloque y19 GCP en Orima ...................................... 192
Figura 157. Ortomosaico con ajuste en líneas de vuelo y 19 GCP en Orima ........................ 193
Figura 158. Ortomosaico con 13 GCP y ajuste en bloque con Z/I Imagine........................... 193
Figura 159. Orden ascendende de la precisión de los mosaicos MDS. .................................. 200
Figura 169. Orden ascendende del RMSEr de todos los ortomosaicos.................................. 211
xxviii
RESUMEN
La cartografía digital generada por métodos fotogramétricos, ofrece una información espacial
precisa y a detalle de una extensión considerablemente amplia del terreno; dicho procedimiento se
denomina proceso fotogramétrico y a su vez cuenta con la aerotriangulación, la misma que utiliza
los puntos de control de tierra o puntos de apoyo fotogramétrico (GCP) para densificar y así obtener
puntos de amarre con coordenadas del terreno; consecuentemente, del ajuste se generan productos
como: modelos digitales de superficie (MDS) y ortofotos, siendo estos últimos evaluados con la
exactitud posicional; por tanto, se analizó los métodos de aerotriangulación en los programas
fotogramétricos Erdas LPS CORE, Orima y Z/I Imagine, aplicando el concepto de la ecuación de
colinealidad y varias funcionalidades, con distribuciones y cantidades diferentes de GCP (13 y
19GCP), utilizando la normativa de exactitud posicional propuesta por el Instituto Geográfico
Militar (IGM) y la nueva normativa ASPRS; obteniendo así, que el mejor resultado está vinculado
a la cantidad de 19GCP y las teorías de Ackermann (1990) e IGN de Perú (2011) en su distribuc ió n;
al tener un producto de 64cm precisión al 90% de confianza con la normativa IGM y de 73cm al
95% de confianza con la normativa ASPRS; por tanto, dicha ortofoto se puede utilizar en la
generación de cartografía 2D a escala 1:5000, hasta escala 1:3400, cuyo umbral corresponde a la
normativa ISO (0.2mm); a su vez, la ortofoto se puede utilizar con fines catastrales, análisis SIG y
exploratorios a trabajos de mayor precisión.
Palabras Clave:
• AEROTRIANGULACIÓN
• CANTIDAD Y DISTRIBUCIÓN PUNTOS DE CONTROL
• EXACTITUD POSICIONAL
• NORMATIVA ASPRS
• ORIMA
xxix
ABSTRACT
Keywords:
• AEROTRIANGULATION
• QUANTITY AND DISTRIBUTION CONTROL POINTS
• POSITIONAL ACCURACY
• NORMATIVE ASPRS
• ORIMA
1
CAPITULO I
1. GENERALIDADES
planteamiento del problema, justificación, descripción del área de estudio, objetivo general,
1.1. Introducción
confiable que ayude en la toma de decisiones; además, dicha información debe contar con
lineamientos técnicos y metodológicos acordes a cada uso (CLIRSEN, 2011). La cartografía digita l
generada por procesos fotogramétricos, ofrece una información espacial precisa y a detalle de una
extensión considerablemente amplia del terreno (Buill, Nuñez, & Rodriguez, 2003), reduciendo así
costos, espacio de almacenamiento y su manipulación resulta mucho más amigable que las
fotogramétrico y a su vez cuenta con un paso denominado aerotriangulación (Hernández L., 2007);
por tanto, el avance de la tecnología hace que la aerotriangulación actual permita generar una
cantidad significativa de puntos de apoyo en los pares estereoscópicos (Falkner & Morgan, 2002),
(Hernández, 2006); sin embargo, los puntos de control de tierra no dejan de ser una parte primordia l
dentro del proceso fotogramétrico, pues son el apoyo para mejorar los resultados de
georreferenciación.
A saber, las ortofotos son un producto del proceso de fotogramétrico que incluye a la
fotogramétrico (Soria, 2014). La evaluación de la calidad las ortofotos según Tsarovski (2015), es
necesaria para dar una representación precisa de la superficie de la tierra tanto como si se midiera
en una representación lineal de un mapa; en dicho sentido, la Universidad de las Fuerzas Armadas
- ESPE en el año 2014 adquiere imágenes digitales, con la finalidad de validar las precisiones
obtenidas en trabajos con cámara digital y la necesidad de investigar los procesos fotogramétr icos
a su época; varias investigaciones se han llevado al cabo en otros países, de las cuales se presenta
a continuación:
3
especificaciones técnicas para escalas pequeñas 1:25000 hasta cartografía urbana escala 1:1000
Con respecto a la calidad de los productos Tsarovski (2015), evaluó la precisión y la calidad de
las ortoimágenes (ortofotos) con la norma NSSDA, generadas a partir de un modelo digital del
terreno plano de 240Km2 con el programa ERDAS Imagine, para generar cartografía 1:2500 con
control para cartografía básica del Instituto Geográfico Agustín Codazzi de Colombia (IGAC),
utilizando la resolución 1392 del IGAG y la “Norma para la exactitud posicional de los datos
Teledetección (ASPRS) para la calidad de las ortofotos y MDS, con la que recomienda una
“Optimizing the distribution of tie points for the bundle adjusment of HRSC image mosaics” siendo
HRSC “High Resolution Stereo Camara”; radica en una distribución homogénea de los puntos de
confiabilidad de las ortofotos y modelos digitales del terreno para lugares donde es imposible la
4
toma de puntos de control tierra, pues esto se implementa en la exploración extra planetaria en
Marte.
su compensación (Hernández L., 2007); así pues, repercute directamente sobre la calidad de los
productos como: modelo digital de superficies (MDS) y la ortofoto; otros aspectos que influyen en
las características de las ortofotos son: desplazamientos del terreno, zonas débiles de correlación,
surcos en edificios y objetos poligonales de altura, no datos (gaps) en MDS, entre otros.
de las ortofotos.
Por tanto, se plantea analizar el método o los métodos adecuados de aerotriangulación para
mejorar la calidad del producto (Ortofoto) en función de la cantidad de GCP y su distribución con
diferentes programas fotogramétricos, como es el caso de ORIMA, LPS Core y Z/I Imagine.
5
generación de la información, entre ellas la aerotriangulación, que con diferentes métodos trata de
dar información confiable para la generación de MDS y posteriormente la ortofoto, dando una
representación fiel del terreno, con medidas y distancias en una imagen digital. El desempeño de
los métodos ha sido estudiado por diferentes autores, pero solo uno ha comparado varios a la vez,
(Liba, Järve, & Rand, 2013) y considerando tanto la cantidad como la distribución de los puntos
programas fotogramétricos como: ORIMA, LPS CORE y Z/I Imagine; así también, en la
Fuerzas Armadas – ESPE, la empresa privada, pública, los profesionales que utilizan o que generan
cartografía; además, de los estudiantes en general, que se beneficiarán de este trabajo con los
La zona de estudio se encuentra al sur de la provincia de Pichincha, dentro de los cantones Quito
Las coordenadas de la zona de estudio se presentan en la Tabla 1, cuyas imágenes del proyecto
han sido tomadas con cámara digital, teniendo un área aproximadamente 102Km2 .
Tabla 1.
Coordenadas en sistema WGS 84 UTM, zona 17S de la zona de estudio
Número Este [m] Sur [m]
W1 776974 9,961680
X1 776980 9,968458
Y1 791972 9,961680
Z1 791978 9,968455
Fuente: Google Earth, Pro 6.0 Fecha: 04/10/2017
7
1.6. Objetivos
geométrica de las ortofotos con ERDAS LPS CORE y ORIMA, en un área comprendida entre los
cantones de Quito y Rumiñahui tomadas con cámara digital, entre el periodo agosto 2012 y enero
2013.
control con una distribución pertinente y uniforme, a través de los métodos de aerotriangulac ió n
• Comparar la calidad geométrica entre los productos generados mediante el uso de tecnología
1.7. Metas
CAPÍTULO II
2. MARCO TEÓRICO
En este capítulo se presenta los conceptos teóricos utilizados en el transcurso del proyecto,
“La Fotogrametría y detección remota es el arte, la ciencia y la tecnología para obtener informac ió n
representación”
de una fotografía o imagen aérea entre dos puntos de vista diferentes (Sánchez, 2007) (Ver Figura
2).
10
Así mismo, los puntos A, B y C en una superficie ∑, se proyectan en un plano llamado S con un
rayo perspectivo Г, formando semirrectas SA, SB, SC; al introducir un segundo planos S2 , los
Al obtener las coordenadas de los objetos A, B, C la fotogrametría siguie los algunos pasos
detallados en la figura 3.
un computador. Las fotografías tradicionales pueden ser escaneadas o a su vez capturadas por una
cámara digital, siendo así otras imágenes digitales (Leica Geosystems, 2005).
Según Schenk (2002), una imagen digital puede describirse como una función continua f(x,y)
(Ver Ecuación 1); donde, las coordenadas x, y son variables espaciales y la función valor
denominado informalmente brillo; en definitiva, una imagen digital se puede expresar por:
en la figura 4.
Según Vozmediano (2006), a cada celda se denominada píxel y está caracterizada por la
intensidad de tono de gris (resolución radiométrica); además, si la imagen está definida por la
combinación de tres bandas del espectro visible: Rojo, Verde y Azul (RGB), se le conoce como
color verdadero.
Según Zapata (2003), existen diferentes tipos de resoluciones en las cámara digitales, entre ellas:
continuación:
- Resolución espacial
Según Zapata (2003), “la resolución espacial indica la mínima medida que se puede distinguir
o apreciar sobre una imagen (ver Figura 5), cuando la resolución espacial de la imagen es mayor,
almacenamiento”
- Resolución radiométrica
de gris que son recogidos por cada banda. En imágenes digitales lo más habitual es que contengan
256 niveles por píxel, cuyo valor de 0 corresponde al color negro y el máximo valor de 255 al
Según Sánchez (2007), las cámaras fotogramétricas tienen una serie de lentes perfectamente
centrados, cuya finalidad es la de medir objetos; por tanto, es necesario corregir algunas
es reemplazado por un sensor con células fotoeléctricas muy pequeñas que registran las imágenes
cuyas siglas son CCD (Change Couple Device); por tanto, según Quirós (2014), la disposición de
los CCD hace posible que exista dos clases de cámaras fotogramétricas digitales: cámaras
En una cámara matricial las imágenes son de geometría estable y conocida; además, se puede
eliminar efectos como la refracción atmosférica, la niebla y es posible trabajar con un programa
Según Quirós (2014), la cámara lineal dispone de tres líneas paralelas con sensores trasversales
a la dirección del vuelo, cuyas inclinaciones son: delantera, nadiral y posterior (ver Figura 7); esta
característica ayuda a que todos los puntos de una imagen se encuentre en tres imágenes.
2.1.6. Paralaje
Según Zapata (2003), el paralaje es “El desplazamiento relativo de objetos inmóviles, desde el
punto de vista de un observador en movimiento que visualiza dichos puntos, en vistas diferentes ”;
superficie de toma y los puntos de vista por las imágenes consecutivas, ver Figura 8.
estereoscópico consecutivo calculado como: diferencia entre el “eje x” de las dos imágenes
consecutivas en un sistema de coordenadas definido por la línea de vuelo y el eje de ordenadas “y”
La condición de colinealidad establece: “Un punto sobre el terreno (A), un punto en la imagen
(a) y el centro de proyección (O) se encuentran en la misma recta” (ver Figura 9) (García & Ortíz,
2018).
16
Con respecto, a la coplanaridad según García y Ortiz (2018), “En esta condición se cumple que
los centros de proyección de dos imágenes contiguas (O1 y O2), dos puntos de imágenes
homólogas (a1 y a2 ) y el punto homólogo de estos en el terreno, A (X’, Y’, Z’), forman parte del
Según Schenk (2002), una imagen digital está afectada por el ruido de la señal, el cambio de
otros; produciendo diferentes niveles de gris, que repercute en la orientación automática de las
Según Schenk (2002), el tamaño del pixel es afectado por la diferencia de escala, producida de
las distintas alturas de vuelo (a); así también, la diferencia de los ángulos de rotación entre los ejes
𝜔𝜔, 𝜑𝜑, 𝜅𝜅 (eje x, eje y, eje z) (b) afecta al tamaño del pixel; dichas afectaciones se representan en la
Figura 11.
(a) (b)
Según Schenk (2002), la distorsión geométrica por inclinación, se produce debido al declive del
puede llegar a ser tan grande que la superficie del objeto se reduce a una línea (Ver figura 12).
Según Schenk (2002) y Sánchez (2007), indican que la distorsión debido al relieve se debe a la
topografía del terreno, como zonas montañosas; a pesar de que la toma sea verticalmente al terreno
Figura 13).
- Refracción atmosférica
que aumenta la altitud, en consecuencia los rayos de luz no viajan de forma rectilínea, sino curvados
según la ley de Snell; por tanto, las coordenadas de imagen deben corregirse para obtener su
refracción atmosférica.
- Curvatura terrestre
Según Lerma (2002), la corrección por curvatura terrestre es un problema de coordenadas del
cartografía se suelen usar sistemas de proyección con referencias altimétricas al nivel medio del
proceso que incluye la toma de imágenes, toma del control terrestre, ajuste fotogramétrico y la
planteamientos metodológicos, que estarán plasmados en un documento llamado contrato (Sato &
2.2.1. Contrato
Según la RAE, 2018: “Un contrato es un convenio, oral o escrito, entre partes que se obligan
sobre materia o cosa determinada, y a cuyo cumplimiento pueden ser obligatorio”. La estructura
depende de la entidad emisora y los requerimientos del cliente; algunos de sus elementos se detallan
en la tabla 2.
22
Tabla 2.
Elementos de un contrato fotogramétrico
No ELEMENTO DESCRIPCIÓN
1. Nombre de quien Es el nombre de una empresa o de una persona natural que representa
realiza el proyecto legalmente al proyecto.
2. Objetivos del Es la actividad a desempeñar, así como el lugar y el fin.
proyecto
3. Marco de referencia Estará incluido el área a cubrir; la escala, tanto de la imagen, como del
levantamiento; el sistema de referencia a utilizar; sistema cartográfico de
representación; equidistancia entre curvas de nivel; elementos a representar;
formatos; precisiones altimétricas y planimétricos (IGM, 2005).
4. Planteamiento Indica todos los detalles metodológicos, técnicos, instrumentales, etc. que serán
Metodológico usados para el trabajo.
4.1. Obtención de la Incluye el diseño plan de vuelo, almacenamiento de las imágenes, un pre vuelo,
fotografía aérea vuelo, post vuelo, procesamiento, evaluación y obtención de productos.
4.2. Control geodésico Indica el control de referencia, planificación de rastreo, control horizontal -
vertical y el procesamiento de la información.
4.3. Clasificación de Todos los detalles discernibles en la fotografía de acuerdo a las normas y
especificaciones vigentes en conformidad de la escala a restituir. Ejemplo:
campo Elementos hidrográficos, avenidas, calles parques, infraestructura, uso de suelo
en áreas no construidas, etc. (IGM, 2005)
4.4. Fotogrametría Incluye la realización de las diferentes orientaciones; restitución; revisión
fotogramétrica; edición y simbolización de cartografía.
4.5. Productos Señala los productos directos del proyecto ajustado por el método planificado,
fotogramétricos así como los productos indirectos que son: el MDE y la ortofoto (IGM, 2008).
5. Valor de los trabajos Incluye el valor del proyecto así como las formas de pago, primeros valores a
liquidar del proyecto y el costo de cada proceso.
6. Plazos Representa el tiempo que durará el proyecto, tomando en cuenta la planificación
de toma de imágenes, control terrestre, procesamiento de la información y
generación de resultados (IGM, 2008).
7. Garantía La entidad que asume la responsabilidad del proyecto y compromiso de seguir
la metodología planteada.
8. Productos y Detalla los productos, formatos, entregables e informes de cada proceso
documentos técnicos realizado, entre otros
9. Anexos Puede ser el cronograma de actividades, mapas base, entre otros.
Una vez realizado el acuerdo entre el beneficiario y el productor del contrato en un proyecto
información como: el uso de las imágenes, la escala del producto, las precisiones y a su vez se
considera la altura de vuelo, base en el aire, separación entre líneas de vuelo, entre otras (Reuter,
Además, un plan de vuelo tiene, entre otras cosas, un bosquejo de lo que será la toma de las
imágenes, la cantidad de imágenes, escala de la foto, altura del terreno, límites de toma, cota
máxima, cota mínima, distancia entre fotos, número de líneas de vuelo, reportes metodológicos de
la toma, sistema de referencia para el apoyo GPS, el método a utilizar, recubrimiento longitudina l,
Tabla 3.
Recubrimiento longitudinal y transversal
Tipo de terreno Traslapo longitudinal Traslapo transversal
El control geodésico o control terrestre, tiene por objeto “Brindar el apoyo terrestre necesario
para realizar la aerotriangulación y ajuste de los bloques de fotografía aérea; en base a éste se
realizará posteriormente la generación del modelo digital del terreno, insumo necesario para la
elaboración de las ortofotos” (MAGAP, 2008); además, el control terrestre se puede determinar
usando procedimientos de topografía como: una red básica obtenida de triangulacio nes,
emplear tecnología como “Global Navigation Satellite System (GNSS)”, en la obtención de las
coordenadas de los centros de las imágenes y los puntos de control. Las consideraciones de los
“Sistema global de navegación por satélite (GNSS) se refiere a un conjunto de satélites que
proporcionan señales desde el espacio, que transmite datos de posicionamiento y tiempo” (ver
Según el Instituto Geográfico Militar (2013), una estación de monitoreo continuo es un conjunto
de infraestructura física y técnica, que se encuentra fija y recoge los datos de varios GNSS a la vez
del punto en el que está ubicado; y en conjunto con otros equipos distribuidos a nivel nacional,
constituye la Red GNSS de monitoreo continuo del Ecuador (REGME) (Ver Figura 19).
Según González (2015), existen diferentes criterios para clasificar los métodos de
posicionamiento, entre ellos: según el sistema de referencia, movimiento del receptor, momento de
obtención de las coordenadas, entre otras. La clasificación más usada se presenta en la Figura 20.
Posicionamiento
Absolut
Relativo
o
Post- Tiempo
DGPS SPS
proceso Real
Movimien
Estático RTK
to
Estático Stop
Cinemático
Rápido and Go
- Posicionamiento absoluto
“Se calcula la posición de un punto utilizando las medidas de pseudodistancias por código (C/A,
L2C o P) con un solo receptor. La precisión del método está bajo los 10 metros” (González, 2015).
- Posicionamiento relativo
Los autores Berrocoso et al. (2009) y la Universidad de Sonora (2008), indican que en el
posicionamiento relativo se mantiene una base fija rastreando, al mismo tiempo uno o varios
receptores toman las coordenadas de los puntos, emitidas por los satélites en un lapso de tiempo
27
determinado; las coordenadas de la base son pre establecidas y las del levantamiento requiere un
- Método estático
Según Farjas (2006), el método estático también es conocido como método estático relativo
estándar, que proporciona precisiones de 5mm + 1ppm, en la que dos o más receptores se
estacionan durante un periodo de tiempo (mínimo media hora) y generalmente se utiliza para medir
distancias mayores a 20 km. Este método se aplica en: redes geodésicas de grandes distancias,
Es una variante del método anterior que reduce los tiempos de observación, generalmente de 5
a 10 minutos por estación, manteniendo las precisiones anteriores, pero la distancia de aplicació n
es de menos de 20 kilómetros (Farjas, 2006); este método es ideal para redes topográficas locales,
Ground Control Point, por sus siglas en inglés GCP, se denomina a los puntos con coordenadas
del terreno y son usados con varios fines, entre ellos, apoyar al ajuste fotogramétrico en la fase de
orientación absoluta (Pérez J. , 2001); además, dichos puntos pueden ser usados en la verificac ió n
dentro del proceso de exactitud posicional por lo que deben aparecer claramente en las imágenes
28
(Ver Figura 21.), dicho proceso es conocido como foto identificación y la obtención usualmente es
con métodos GNSS (Zapata, 2003); además, la cantidad y distribución de dichos puntos puede
variar; tal y como se indicará en el apartado 2.5 “Distribución y cantidad de puntos de apoyo
fotogramétrico (GCP)”.
para la reconstrucción de la geometría de los objetos, cuya precisión geométrica depende en gran
medida de la similitud del haz de imágenes que ingresan durante el momento de la toma. Para un
mismo procedimiento, pueden existir varios modelos matemáticos que se diferencian en qué tan
cerca describen los procesos físicos. En este capítulo se describen estos procedimientos y los
“Aerotriangulación”.
Una vez definidos los objetivos del proyecto, la toma de imágenes y tomados los GCPs de apoyo
terrestre, normalmente se procede a ajustar del proyecto; por tanto, se requiere de datos de
29
calibración de la cámara para realizar la orientación interior y posteriormente las siguie ntes
orientaciones (Schenk, 2002); al mismo tiempo, se debe considerar el tipo de cámara a usar según
su características (Ver apartado 2.1.5), a continuación una breve descripción del contenido de un
Según Crouse et al. (2018), el certificado de calibración es un documento con las constantes de
calibración, con las cuales se evalúa la calidad de la toma de una cámara sobre un objetivo, un
ejemplo es este certificado se ilustra en la Figura 22, cuyos elementos son: número de filas y
columnas con unidad de pixel (1), tamaño del pixel en µm de la cámara (2), distancia focal (3),
coordenadas del punto principal de autocolimación (4), calibración geométrica de los ases de luz
Lerma (2002), define que la orientación interna es “La transformación entre el sistema de
que implica una traslación del sistema de coordenadas”(Citado en Angulo, 2014, p. 27), los
elementos de un sistema de coordenadas imagen – cámara son los siguientes: centro fiducial (FC),
punto principal (PP), punto de simetría (PS), distancia focal (c) y vector de la imagen (p), (ver
Figura 24.).
31
interna de la cámara en el momento de captura de los datos, teniendo en cuenta las variables de la
imagen desde la perspectiva de una cámara aérea” (Leica Geosystems, 2005). En las imágenes
tomadas con cámaras digitales, no se encuentran las marcas fiduciales y los parámetros de
orientación son determinados con el certificado de calibración (ver Figura 25) (Quirós, 2014).
Según Coelho y Nunes (2007), la orientación exterior es “La obtención de la posición y altitud
referencia”. Una imagen se encuentra orientada exteriormente cuando se conocen los seis
parámetros de orientación: coordenadas espacio-objeto para el centro de perspectiva (𝑋𝑋�, 𝑌𝑌�, 𝑍𝑍̅ ) y
ángulos de rotación del sensor (φ,ω,κ). (Citado en Angulo, 2014, p. 28) (Ver Figura 26).
Dicho de otra manera, la orientación exterior implica la orientación del sensor con los
parámetros de donde se tomó la imagen (𝑋𝑋�, 𝑌𝑌�, 𝑍𝑍̅ ) y aquellos que dan la inclinación de la cámara
(φ,ω,κ); según Cheli (2011), “esta georreferenciación del espacio-imagen es conocida como
La orientación relativa orienta cada objeto de una imagen con su homólogo posterior,
reconstruyendo la posición exacta de un par estereoscópico (Coelho & Nunes, 2007). Para la
formación del modelo, es necesario medir una serie de puntos homólogos entre imágenes (Puntos
de paralaje o puntos von Gruber) y seguidamente calcular los cinco parámetros de orientación (φ,
ω, κ, Tx, Ty, Tz) (Schenk, 2002). Al terminar de identificar dichos puntos se procede a reconstruir
los haces mediante la condición de coplanaridad, (ver Figura 27) (Quirós, 2014).
Son puntos de paso entre modelos estereoscópicos, mencionados en 1968 por Von Gruber en su
medir al menos seis puntos bien distribuidos por el modelo (Springer, 2009) (ver Figura 28).
34
La distribución de los puntos sugiere colocar dos en la parte superior, dos en la parte central y
dos en la parte inferior del modelo; a partir del segundo modelo, la imagen tendrá nueve puntos de
paso.
con respecto al terreno, determinado por siete parámetros, resultantes de las siete incógnitas que se
resuelven con ayuda de los puntos de control, medidos en el terreno e identificables en los
fotogramas y según Coelho y Nunes (2007), los parámetros se determinan de la siguiente manera:
“X0 , Y0 y Z0 representan la posición del centro de perspectiva en los ejes X, Y y Z del sistema
de coordenadas del espacio-objeto. Los ángulos (𝜙𝜙, 𝜔𝜔, 𝜅𝜅), a su vez, son rotaciones del sistema
(phi) representa la rotación del eje y en relación a Y; estos ángulos deben ser pequeños, no
35
último, κ (kappa) representa la rotación del eje z con respecto a Z, dichos ángulos son medidos
2.3. Aerotriangulación
También es conocida como triangulación aérea, tiene sus orígenes con el invento de la marca
en 1901, que fue la base de la fotogrametría analógica; no obstante, los costos de producción y el
uso de instrumentos ópticos-mecánicos, hacían del proceso de triangulación aérea algo engorroso
algunos de sus procesos fue la aerotriangulación, pues mejoraba sustancialmente la medición óptica
36
y el tiempo en los cálculos de los modelos para densificar puntos; esta etapa es conocida como
El paso de fotogramas a imágenes digitales hizo que los restituidores analíticos sean
remplazados por las estaciones fotogramétricas digitales (EFD), pues su desempeño era
La aerotriangulación es una buena demostración del gran potencial de una EFD, que ha
automatizado varios de sus procesos, haciendo que la misma sea aceptada y demostrada; sin
(Schenk, 2002). Con tales antecedentes se presenta a continuación varias definiciones y modelos
para la extensión del control horizontal y/o vertical, mediante las mediciones de ángulos y
Con el paso de los años, el término aerotriangulación o triangulación aérea se hizo más común
puntos; la cantidad de triángulos hace que sea un método de densificación de puntos (Lerma,
37
2002); así mismo, estos puntos se conocen como puntos de unión o de paso, que al estar medidos
en cada foto, se asocian con los puntos de control y sus coordenadas pueden ser determinadas con
fotografías aéreas o imágenes digitales, con el mínimo de trabajo de campo, aprovechando sus
propiedades geométricas; estas mediciones se realizan con equipos fotogramétricos digitales que
aumentar el número de puntos de control y de apoyo, para su uso en diversas aplicacio nes
criterios como: instrumentación, apoyo GPS, tipo de bloque, cantidad de fotos/imáge nes,
ecuaciones, resolución de ecuaciones, entre otros (Cheli, 2011); varias clasificaciones se presentan
a continuación:
38
Por su traducción al español de “American Society of Photogrammetry (1980)”, indica que los
A. Triangulación Radial.- Es un método gráfico para establecer coordenadas del terreno (ver
Figura 30), que utiliza la fotografía aérea como la unidad; se basa en las diferencias de
elevación del terreno, y el cambio de escala que producen los desplazamientos radiales de
Photogrammetry, 1980).
operador a realizar su trabajo; a su vez, se subdivide según las siguientes consideracio nes :
tanto, se divide en: modelos independientes y haces de rayos; en el modelo independiente considera
por separado el ajuste planimétrico (X, Y) y altimétrico (Z), cuya combinación se lo denomina
formado por dos a tres imágenes (ver Figura 31) y las coordenadas de los puntos de paso se
combinan para compensar el bloque. Matemáticamente pueden ser sistemas libres o ligados;
es decir, puntos de apoyo matemáticamente perfectos y exentos de error o con cierto nivel
B. Método de haces de rayos.- Considera el rayo (óptico) espacial como una unidad de trabajo,
cuyo sistema de ecuaciones utiliza la condición de colinealidad, (Lerma, 1999), ver Figura
32.
instrumental apropiado; tal es el caso del uso de otras tecnologías como los sistemas GNSS;
por tanto su clasificación se debe a con control geodésico y sin control geodésico.
C. Respecto al modo del operador.- considera una clasificación planimétrica y una espacial;
haces resuelve con las ecuaciones de colinealidad en todo un bloque fotogramétrico (Sánchez,
todo el bloque completo, con una transformación de 7 parámetros y el cálculo por métodos
iterativos.
plantean para todo un bloque fotogramétrico y la unidad base es el fotograma; es decir, las
que ver con su pasado (fotogrametría analítica); algunos de ellos se resumen en la tabla 4.
43
Tabla 4.
Clasificaciones de la aerotriangulación por diferentes autores
Autor Clasificación Descripción
AMERICAN SOCIETY 1. Triangulación Radial Método gráfico utilizando en la fotografía aérea como
OF base.
PHOTOGRAMMETRY, 2. Triangulación A su vez se puede dividir en:
1980 analítica. 1. Procedimientos de triangulación
2. Tipo de ecuaciones de condición
3. Método de resolver las ecuaciones de
condición
3. Estereotriangulación Considera la instrumentación utilizada.
lerma, 1999 1. Modelos Se considera la unidad de trabajo el modelo
independientes estereoscópico.
1.1. Aerotriangulación Utiliza puntos de paso y de apoyo, considerando un
planimétrica sistema libre o ligado, para su compensación en
coordenadas X,Y.
1.2. Aerotriangulación Utiliza puntos de paso y de apoyo, considerando un
altimétrica sistema libre o ligado, en su compensación solo en la
coordenadas Z.
1.3. Aerotriangulación Combina la aerotriangulación planimétrica y altimétrica
tridimensional para su compensación. Método resuelto con matrices y
métodos iterativos.
2. Aerotriangulación Considera el rayo (óptico) espacial como una unidad y
por el método de formula el sistema de ecuaciones utilizando la condición
los haces de rayos de colinealidad.
Përez j., 2001 1. Respecto a la Considera la disposición del fotograma: modelos
unidad de trabajo aislados, pasada simple, bloque de pasadas y bloque de
secciones.
2. Respecto al control Considera el apoyo de otras tecnologías: con control
terrestre geodésico y sin control geodésico.
3. Respecto al modo Es la clasificación más utilizada.
del operador.
3.1. Planimetría o A su vez se puede clasificar en: gráfica, numérica y
Radial mecánica
Continúa
44
Según Sánchez (2007), el ajuste por modelos independientes utiliza como unidad de medición
el modelo estereoscópico (ver Figura 34), cuyo sistema de coordenadas es único para todos los
modelos y suele empezar con el primer modelo; al principio se realiza el paso de puntos de un
modelo a otro, aprovechando el recubrimiento; además, se requieren los centros de proyección que
ayudan la disminución de la deformación por los giros. Las coordenadas se pasan consecutivame nte
El ajuste por modelos independientes se puede realizar de tres maneras: ajuste planímetr ico,
ajuste altimétrico y ajuste tridimensional (ver Figura 35). En el ajuste planímetro y altimétrico se
usan puntos de paso, puntos de apoyo y las coordenadas de los centro de proyección de las
En el ajuste tridimensional, cada modelo puede trasladarse, rotarse y cambiar de tamaño; pero
la estructura interna del modelo, permanecerá invariable (Sánchez, 2007); además, según Pérez J.
(2001) un ajuste secuencial de un bloque debe conocerse las coordenadas de los cuatro puntos de
las esquinas y los dos centros de proyección del modelo; las coordenadas requeridas se obtienen
46
𝑋𝑋 𝑥𝑥 𝑖𝑖 𝑇𝑇𝑥𝑥
[ ] 𝑦𝑦 𝑇𝑇
� 𝑌𝑌� = 𝜆𝜆𝑗𝑗 𝑅𝑅 𝑗𝑗 � 𝑖𝑖 � + � 𝑦𝑦 �
𝑍𝑍 𝑧𝑧𝑖𝑖 𝑇𝑇𝑧𝑧 [ 2]
Donde:
coordenadas.
Estas ecuaciones permiten definir la función F (ecuación 3), que integra los siete parámetros de
la transformación (λ, Ω, Φ, Κ, TX, TY, TZ), con las observaciones (xi, yi, zj).
𝐹𝐹𝑥𝑥 𝑥𝑥 𝑖𝑖 𝑇𝑇𝑥𝑥 𝑋𝑋
𝐹𝐹 [ ] 𝑦𝑦 𝑇𝑇
� 𝑦𝑦 � = 𝜆𝜆𝑗𝑗 𝑅𝑅 𝑗𝑗 � 𝑖𝑖 � + � 𝑦𝑦 � − �𝑌𝑌 �
𝐹𝐹𝑧𝑧 𝑧𝑧𝑖𝑖 𝑇𝑇𝑧𝑧 𝑍𝑍 [ 3]
Estas ecuaciones nos son lineales y por tanto, se tienen que someter a un proceso de
linealización.
47
Según Coelho y Nunes (2007), el método de ajuste de haces es el único que puede ser
considerado totalmente analítico, debido a que el operador interviene muy poco en el proceso;
además, es un método que utilizan la mayoría de los restituidores analíticos y soluciones digita les
actuales; así pues, se define como un método que permite la reconstrucción óptima de la geometría
de los haces de rayos luminosos formados por el sistema espacio-imagen en el momento de la toma
de las mismas.
El proceso utiliza el Método de los Mínimos Cuadrados (MMCC) para minimizar los residuos
del ajuste, de forma que cada radio ajustado pase lo más cerca posible de su homólogo, en el
espacio-imagen y del centro de perspectiva de la imagen (ver Figura 36) (Coelho & Nunes, 2007).
Esta condición de colinealidad puede ser expresada por las ecuaciones 4 y 5, por cada eje.
48
Donde:
• (ξ, η), son las coordenadas del punto en el sistema Espacio-Imagen [mm].
• (rnm), es el elemento de la matriz de rotación entre los sistemas XYZ y X'Y'Z ' (este
Así mismo, Coelho y Nunes (2007), indican que las funciones de colinealidad a pesar de
describir una recta, no son lineales, pues combinan parámetros lineales y angulares. El proceso de
forma iterativa, hasta que las diferencias estén dentro del límite deseado; es decir, de tolerancias
fijadas a priori, tanto para los valores lineales, como para los angulares.
2.3.6. Convergencia
o valor. “En las matemáticas la convergencia es una propiedad de las sucesiones y series; la una
49
tiende a un límite, mientras que en la otra, la suma de sus términos, se aproxima cada vez a una
determinada cantidad” (RAE, 2018). El valor de convergencia se usa como un umbral para
2.3.7. Restitución
Según Sánchez (2007), Coelho y Nunes (2007), la restitución fotogramétrica no es más que la
o artificiales presentes en el terreno, a fin de generar cartografía; todo esto con la ayuda del
principio de la marca flotante, realizadas en un restituidor digital, que tiene como característica
Los resultados de un ajuste fotogramétrico según Quirós (2014), son: cartografía 2D y 3D;
producto antes que los otros, pues su obtención se basa en productos anteriores. A continuación se
Según Durón (2007), un modelo digital de superficie (MDS) es una matriz de datos de elevación,
que representa los elementos del terreno y del suelo son obtenidos por métodos de interpolac ió n;
además, representa objetos como vegetación, infraestructura (edificios, puentes, viaductos, entre
Así mismo, según Brito, Prado y Augusto (1999), indican que el MDS puede definirse como:
“Cualquier representación numérica de una determinada superficie física sobre el del terreno”
Cuando un MDS expresa altitud, se llama el modelo numérico de elevaciones (MNE) o el modelo
digital de elevaciones (MDE o DEM, Digital Elevation Model) (Citado en Coelho y Nunes, 2007,
p. 157).
Figura 38. Diferencias entre DTM (en español MDT) y DSM (en español MDS)
Fuente: (GIS & Remote Sensing Service, 2018)
51
Según Felicísimo (1999), un MDE es un grupo de datos, cuyo valor numérico representa la
superficie topográfica, mediante la repartición de puntos de altitud en el terreno, cuya unidad mínima
está dada por la variable Z (Citando en Barragán, 2015, p. 39) (ver Figura 39).
Según Ormeño (2004), Niño (2009) y Arranz (2010), indican que un modelo digital del terreno
agua, expresado en coordenadas x,y,z de los puntos que lo definen; además, considera la
información de la restitución fotogramétrica como ríos, líneas de cresta, líneas de perfil del terreno,
entre otros elementos; en el caso del terreno desierto un MDS es muy cercano (sino el mismo) que
Barragán (2015), indica que una imagen corregida en fotogrametría tiene un ajuste geométrico,
a continuación:
• Ortofoto digital
Según Schenk (2002) y Lerma (2002), indican que una Ortofoto Digital u Ortoimagen es una
representación ortogonal del terreno y resulta del uso de un MDT generado automáticamente; al
mismo tiempo, permite ajustes radiométricos de las imágenes que intervienen en el proceso; en
otras palabras, es una representación digital de una superficie de terreno, que permite realizar
• Ortomosaico
necesario juntarlas, sobre todo cuando tienen cierto grado de recubrimiento entre ellas (ver Figura
41); esto se realiza con el fin de tener una sola información, a este proceso se lo denomina mosaico;
no obstante, requiere de un ajuste radiométrico para homogenizar sus colores y evitar que dicha
Según Pozo (2002) y Gómez (2011), indican que el avance de la fotogrametría ha llevado a la
que realizan cálculos de una manera rápida y sencilla; además, permiten la detección de errores por
métodos estadísticos.
Coelho y Nunes (2007), señalan algunos programas como: Erdas Imagine y PCI OrthoEngine,
de corta distancia, aerotriangulación; sin embargo, existen otras opciones como: LH-Systems, ISM,
KLT, DAT/EM, VinGeo, DVP, TopoL, VirtuoZo, Autometric/Boeing, E-foto, Z/I Imaging, y
Orima.
La influencia de estos programas hace que se priorice la investigación de éstos para dar mayor
confiabilidad al usuario sobre el producto a obtener (Topoequipos S.A, 2015). A continuación una
2.4.1. Erdas-LPS
Erdas Imagine, es un programa que es utilizado en diferentes ramas del conocimiento como: la
fotogrametría y los sensores remotos, aplicando tecnología LiDAR, análisis vectorial, entre otras
herramientas (Hexagon, 2018); además, cuenta con Leica Photogrammetry System (LPS), que es
un paquete de herramientas fotogramétrico del mismo Erdas Image, diseñado para medir GCP en
múltiples imágenes, realizar productos como: ortofoto, MDS, aerotriangulación de múltip les
imágenes (ver Figura 42), entre otras funcionalidades (Erdas, 2011). La versión 2014 al módulo
Según Erdas (2010), la versión 11.0.2 de LPS contiene un paquete de aerotriangulación que se
Donde:
• (A), matriz de las derivadas parciales con respecto a los parámetros desconocidos;
Además, existen varios modelos matemáticos al considerar la matriz de pesos (P) en un ajuste
por bloque resuelto por la ecuación 6; por ejemplo: modelo simple de Bauer (3 parámetros), modelo
simple Jacobsen´s (4 parámetros), modelo ortogonal Ebner´s (12 parámetros), etc. (Ver tabla 5).
Tabla 5.
Parámetros de adicionales de ajuste.
Nombre del modelo Descripción
No additional Sin uso de parámetros, esta es la opción predeterminada.
parameters
Bauer’s simple model Un modelo de tres parámetros del Dr. Bauer, en donde, dos parámetros son para la
deformación afín (causada por la no ortogonalidad) y parámetro para la distorsión
simétrica de la lente.
Jacobsen’s simple Un modelo de cuatro parámetros del Dr. Jacobsen, en el que se compensan la distorsión
model afín, distorsiones de primero y sengundo; además de la distorsión de la lente.
Ebner’s orthogonal Es un modelo de doce parámetros del Dr. Ebner, que compensa la mayoría de errores
model sistemáticos, pero también requiere más observaciones por lo que requiere un mayor
número de GCP y al menos nueve puntos de amarre por imagen.
Brown’s physical model Es un modelo de catorce parámetros del Dr. Brown, y puede compensar la mayoría de
distorsiones lineales y no lineales de la película y la lente, pero requiere de más puntos
de control y de enlace.
Lens distortion model Es un modelo de dos parámetros diseñado para compensar errores de auto calibración
de la lente, generalmente para cámaras de aficionados.
Fuente: (Erdas, 2011).
Así mismo, un indicador de la calidad del ajuste es el valor del error cuadrático medio total
(RMSE), que está basado en los residuales de las coordenadas de imagen y los residuales de las
2.4.2. Orima
Según Hinsken (2010), autor de Orima, indica que dicho programa es capaz de gestionar algunos
colinealidad, pues tiene en cuenta la integración con GPS/IMU (ver apartado 3.2.4) a los sensores
58
(ejemplo: cámara ADS40), esto con el fin de obtener la mayor precisión posible (Hinsken & Miller,
2002); el modelo antes descrito presenta el trabajo de Müller (1991);donde, las coordenadas de
imagen se expresan como una función del punto en tierra (Pi) y los parámetros de orientación de
Por tanto, las ecuaciones en sus ejes están definidos como la ecuación 7 en el eje x y la ecuación
𝑥𝑥 𝑖𝑖𝑖𝑖 = 𝐹𝐹𝑖𝑖𝑖𝑖 (𝑋𝑋𝑖𝑖 , 𝑌𝑌𝑖𝑖 , 𝑍𝑍𝑖𝑖 , 𝑋𝑋𝑘𝑘 , 𝑌𝑌𝑘𝑘 , 𝑍𝑍𝑘𝑘 , 𝜔𝜔𝑘𝑘 , 𝜑𝜑𝑘𝑘 , 𝜅𝜅𝑘𝑘 , 𝑋𝑋𝑘𝑘 +1 , 𝑌𝑌𝑘𝑘+1 , 𝑍𝑍𝑘𝑘+1 , 𝑋𝑋𝑘𝑘 +1 , 𝑌𝑌𝑘𝑘+1 , 𝑍𝑍𝑘𝑘+1 , 𝜔𝜔𝑘𝑘 +1 , 𝜑𝜑𝑘𝑘+1 , 𝜅𝜅𝑘𝑘+1 ) [ 7]
𝑦𝑦𝑖𝑖𝑖𝑖 = 𝐺𝐺𝑖𝑖𝑖𝑖 (𝑋𝑋𝑖𝑖 , 𝑌𝑌𝑖𝑖 , 𝑍𝑍𝑖𝑖 , 𝑋𝑋𝑘𝑘 , 𝑌𝑌𝑘𝑘 , 𝑍𝑍𝑘𝑘 , 𝜔𝜔𝑘𝑘 , 𝜑𝜑𝑘𝑘 , 𝜅𝜅𝑘𝑘 , 𝑋𝑋𝑘𝑘 +1 , 𝑌𝑌𝑘𝑘+1 , 𝑍𝑍𝑘𝑘+1 , 𝑋𝑋𝑘𝑘 +1 , 𝑌𝑌𝑘𝑘+1 , 𝑍𝑍𝑘𝑘+1 , 𝜔𝜔𝑘𝑘 +1 , 𝜑𝜑𝑘𝑘+1 , 𝜅𝜅𝑘𝑘 +1 ) [ 8]
Donde:
Además, el modelo matemático se acopla muy bien en áreas sin mediciones suficientes de
puntos de enlace (por ejemplo lagos o bosques), al determinar los parámetros de orientación con
El programa Orima cuenta con diferentes licencias de trabajo detalladas en la Figura 45 y que
son utilizadas para cámara ADS40, imágenes con marcas fiduciales e imágenes con apoyo
IMU/GPS incorporado.
Z/I Imaging Corporation es una empresa que ofrece desarrollo, soporte y soluciones de
imágenes terrestres, incluidas cámaras aéreas, estaciones de trabajo estéreo, analítico y escáneres
fotogramétricos. La compañía fue fundada en 1999 y tiene sede en Huntsville, Alabama (Global
Market, 2018).
60
triangulation”, que parte de los estudios de la empresa Carl Zeiss de la división fotogramétrica de
la empresa Z/I Imaging GmbH; cuyos procesos se basan en la medición de puntos de manera
La distribución y cantidad de puntos de control varían según las especificaciones del proyecto
y de la tecnología de soporte (apoyo aéreo cinemático) (Muñoz, 2004); así pues, varios autores
como: Otero y Sevilla (1989), Pérez J. (2001) y Cartesia (2010), indican que un plano se define
con tres puntos y concuerdan en que un modelo se encuentra bien orientado teóricamente, cuando
utilizan dos puntos con coordenadas X,Y,Z y para escalar el modelo se requiere de un punto con
coordenada Z; sin embargo, una mejor orientación requiere de dos puntos más, que aumentará la
redundancia de las ecuaciones; así como, sus parámetros (Obregón, 2007), siendo su distribuc ió n
Punto X, Y, Z
Punto Z
modelos; así pues, los puntos de control según Hussain y Bethel (2004), se clasifican en horizonta l
y vertical, en donde los puntos de control horizontal, normalmente, se ubican a lo largo del
perímetro del bloque, con un espaciamiento promedio de 3 a 5 modelos o más; por otro lado, la
distribución del control vertical se realiza a lo largo de la línea de vuelo, debido a la deformació n
Figura 47.
IMU incorporado en el aire (para medir la posición y altitud de cada imagen) en la etapa de
McGlone, 2001), logrando una distribución como la Figura 48; donde, se muestra puntos de control
No obstante, existen varias formas de distribuir los puntos de control, así como su cantidad, los
Tabla 6.
Distribución y cantidad de puntos de apoyo fotogramétrico.
Descripción Ilustración
Continúa
63
(Ackermann, 1990)
A B
(Ackermann, 1994)
Continúa
64
tener por lo menos un punto en cada tercera imagen y (Leica Geosystems, 2005)
principalmente en zonas de múltiple traslapo.
(▲), puntos de control.
(IGN, 2011)
Continúa
65
(PIX4D, 2018)
aproximada en una imagen, para su posterior obtención en campo (ver Figura 49); se apoya con
Los criterios para elegir dichos puntos, según Pérez J. (2001) y Heipke (1997), son:
• Si es posible, puntos ubicados en líneas rectas, esquinas definidas, cruces de caminos, bien
definidos.
66
• Puntos visibles desde varias direcciones bien distribuido a través de las imágenes y de fácil
los diferentes metadatos de un punto característico de la zona, que está presente en las imágenes
foto identificables. Algunos de los datos que debe contener son: coordenadas geográficas,
coordenadas planas, altura, ubicación geopolítica, orden, descripción geográfica del punto, tipo de
monumentación, vías de acceso, fecha de determinación, proyecto, autores, gráficos, croquis, entre
Según Ariza, Emerson, Chicaiza y Buenaño (2017), la exactitud posicional es parte de la calidad
Geomática; así mismo, un control de calidad del proceso fotogramétrico, disminuye los posibles
errores que se acumulan a medida que se avanza en las diferentes fases (Quirós, 2014).
Algunos de los estándares más conocidos al evaluar cartografía son: EMAS, NMAS, NSSDA,
denominado ASPRS “Positional Accuracy Standards for Digital Geospatial Data” (2014); por
ASPRS, es una organización responsable del estándar “ASPRS Positional Accuracy Standards for
Digital Geospatial Data”, entre sus objetivos tiene: reemplazar los estándares de precisión ASPRS
existentes para mapas a gran escala (1990) y preparar informes verticales de precisión para datos
Lidar (2004) (ASPRS, 2014); así mismo, la ASPRS incluye estándares de precisión posicional de
ortoimágenes digitales, datos planimétricos digitales y datos de elevación digitales, basadas en los
valores del error cuadrático medio (root mean square error - RMSE).
68
Vale destacar, que la norma incluye tablas con umbrales de RMSE para generar productos con
las tendencias actuales de precisiones basadas en el GSD; además, de ecuaciones que determina n
mismo, la cantidad de puntos de muestreo en la exactitud posicional según el área a cubrir. Algunos
Según el IGM (2006), las imágenes deben ser exentas de manchas, partículas de polvo, cabellos,
rayas, marcas fiduciales visibles, entre otros; además, las características radiométricas del
histograma de color de cada banda, deben seguir una distribución normal y la saturación de los
orientaciones; por tanto, existen diversas tolerancias al realizar la orientación interior, relativa,
absoluta y aerotriangulación.
Según IGM (2006), las imágenes digitales obtenidas con cámara analógica y escaneadas, el error
estándar en la orientación interior es de 0.6 de pixel; también, Quiroz (2014), indica que si el
fotograma ha sido obtenido con cámara analógica el error medio máximo no debe ser mayor de un
69
pixel y cuando se utiliza cámara digital la orientación interior no precisa de tolerancias, pues su
orientación es automática.
relativa y orientación absoluta; según el IGM (2006), el error estándar para el ajuste relativo es
El Instituto de Información Territorial del Estado de Jalisco (2008), explica que los puntos
aerotriangulados no debe excederse más allá de la cubierta de los puntos de control y un ajuste por
Además, la norma Positional Accuracy Standards for Digital Geospatial Data de la ASPRS
(2014) indica que la precisión de la aerotriangulación para datos planimétricos digita les
(ortoimágenes o mapas planimétricos digitales), utilizando el RMS del mapa y el tamaño del GSD
Tabla 7.
Ejemplos de calidad para datos planimétricos digitales de alta precisión.
ASPRS Positional Accuracy Standards for Digital Geospatial Data 2014 Equivalente a
Clase de precisión RMSEr Precisión horizontal en el GSD aproximado la escala del
horizontal RMSEx y (cm) nivel de confianza del del origen de las mapa en
RMSEy (cm) 95% (cm) imágenes (cm) NMAS
15.0 21.2 36.7 7.5 to 15.0 1:380
25.0 35.4 61.2 12.5 to 25.0 1:634
30.0 42.4 73.4 15.0 to 30.0 1:760
45.0 63.6 110.1 22.5 to 45.0 1:1,141
60.0 84.9 146.9 30.0 to 60.0 1:1,521
Continúa
70
se reemplaza los valores en las ecuaciones 9 y 10, tanto para aerotriangulación planimétrica o
1 [ 9]
𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑋𝑋 (𝐴𝐴𝐴𝐴 ) 𝑜𝑜 𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑦𝑦 (𝐴𝐴𝐴𝐴 ) = ∗ [𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑋𝑋 (𝑀𝑀𝑀𝑀𝑀𝑀𝑀𝑀) 𝑜𝑜 𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑦𝑦 (𝑀𝑀𝑀𝑀𝑀𝑀𝑀𝑀)]
2
Donde:
dirección Este.
dirección Norte.
altura.
• RMSE(Mapa), representa el error cuadrático medio máximo del producto final a obtener.
71
Según IGM (2006), los puntos de apoyo fotogramétrico (GCP) deben ser obtenidos con
receptores GPS geodésicos de una frecuencia L1 o de doble frecuencia L1/L2 y que son de código
Las precisiones de estos puntos en planimetría será del error medio cuadrático de las
coordenadas, cuyo valor máximo será de 0.1mm por el módulo de la escala del mapa a obtener
(IGM, 2013); sin embargo, la norma ASPRS (2014), precisa el valor de los puntos de control
1 [ 11]
𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑋𝑋 𝑜𝑜 𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑦𝑦 = ∗ [𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑋𝑋 (𝑀𝑀𝑀𝑀𝑀𝑀𝑀𝑀) 𝑜𝑜 𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑦𝑦 (𝑀𝑀𝑀𝑀𝑀𝑀𝑀𝑀)]
4
1 [ 12]
𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑧𝑧 = [𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑋𝑋 (𝑀𝑀𝑀𝑀𝑀𝑀𝑀𝑀) 𝑜𝑜 𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑦𝑦 (𝑀𝑀𝑀𝑀𝑀𝑀𝑀𝑀)]
2
Donde:
• (RMSEx), representa el error cuadrático medio máximo de los GCP en la dirección Este.
• (RMSEy), representa el error cuadrático medio máximo de los GCP en la dirección Norte.
• (RMSE(Mapa)), representa el error cuadrático medio máximo del producto final a obtener.
72
Según el IGM (2006), las ortofotos digitales que se generan con un MDE de cobertura del 5%
Tabla 8.
Especificaciones técnicas de ortofotos digitales.
Escala de Tamaño del Escala de Tamaño del Precisión Precisión
la pixel en la toma ortofotos pixel en el horizontal (m) vertical (m)
fotografía (m) producto (m)
1:5000 0.07 1:1000 0.10 0.30 0.25
1:10000 0.14 1:2000 0.20 0.60 0.50
1:20000 0.28 1:2500 0.25 0.80 0.60
1:30000 0.42 1:5000 0.50 1.50 1.25
1:60000 0.84 1:10000 1.00 3.00 2.50
Fuente: (IGM, 2006).
Según varios autores como: Mikhail (1976), Goodchild y Gopa (1989), Li (1991), Leung y Yang
(1998), entre otros autores, detallados por Ariza y Atkinson (2006), el tamaño de la muestra para
datos cartográficos y posicionales, asumen una distribución normal; además, dicho autor señala
muestras recomendados por diferentes normas, entre ellas NMAS, EMAS, NSSDA, ASPRS
indican un n=20; sin embargo, la ASPRS “Positional Accuracy Standards for Digital Geospatial
Data” indica: “el número de puntos varía según el área a cubrir” lo que se detalla en la tabla 9.
73
Tabla 9.
Número recomendado de puntos de control basado en la cantidad de área.
Area del Proyecto Prueba de precisión horizontal de una Ortoimagen y planimetría.
(km2 ) Número total de puntos de control 2D / 3D estáticos (puntos
claramente definidos)
≤500 20
501-750 25
751-1000 30
1001-1250 35
1251-1500 40
1501-1750 45
1751-2000 50
2001-2250 55
2251-2500 60
Fuente: traducido de (ASPRS, 2014)
en una hoja cuadrada de 15 segundos de latitud y longitud; además, el total de hojas se conoce
como tamaño de la población (IGM, 2017). La precisión de estos puntos deben ser obtenidos por
medios más precisos de hasta un cuarto del producto a obtener (ASPRS, 2014).
según IGM (2017), para generar cartografía con fines catastrales y a escala 1:1000, se debe
determinar el residuo o diferencia, de los puntos obtenidos en la Ortoimagen respecto del control
efectuado en campo con métodos más precisos, tal y como indica la norma NSSDA, al aplicar el
𝑛𝑛
1 2
𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅 = � ��𝑥𝑥 𝑖𝑖(𝑚𝑚𝑚𝑚𝑚𝑚𝑚𝑚) − 𝑥𝑥 𝑖𝑖(𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 ) �
𝑛𝑛
𝑖𝑖 =1 [ 13]
Donde:
conjunto de datos.
A continuación, según IGM (2017), se debe calcular el RMSEx y el RMSEy considerando dos
casos; el primero caso implica RMSEx=RMSEy; así mismo, tiene en cuenta el valor de 2.1460 al
90% de confianza, que al dividir para la raíz de dos resulta la ecuación 14 (FGDC, 1998).
Por tanto, el segundo caso implica RMSEx ≠ RMSEy; y al 90% de confianza el mismo valor de
2.1460; por tanto, se multiplica por 0.5, la suma del RMSEx y RMSEy (la suma de una
aproximación del error circular, considerando si el RMSEmin /RMSEmax , se encuentra entre 0.6 y 1
Si es el caso que los valores no se encuentran entre dicho límite, se aplica el caso uno; donde, a
Sin embargo, la norma “Positional Accuracy Standards for Digital Geospatial Data” de la
posicional en unidades de centímetros (ver tabla 10), que se evalúan con el RMSE en el plano
Tabla 10.
Estándares de precisión horizontal para datos geoespaciales
Componente Exactitud Absoluta
Horizontal de RMSEx y RMSEr (cm) Precisión horizontal al 95%
precisión RMSEy (cm) de nivel de confianza (cm)
Además, existe un criterio para evaluar la componente de altura tal y como indica la tabla 11.
Tabla 11.
Estándares de precisión vertical para datos de elevación digital
Clase de precisión vertical Exactitud Absoluta
RMSEz Sin Vegetación NVA at 95% Nivel de confianza (cm)
(cm)
X-cm ≤X ≤1.96*X
Fuente: Traducido de (ASPRS, 2014).
76
A continuación se presentan ejemplos de RMSE con la norma ASPRS con confiabilidad al 95%,
tamaño del pixel aproximado y la escala del mapa a obtener (ver tabla 12).
Tabla 12.
Precisión horizontal/ejemplos de calidad para datos planimétricos digitales
ASPRS 2014 Equivalencia del Equivalenci
mapa en: a de la
escala del
Clase de RMSEr Precisión Tamaño del GSD ASPRS ASPRS mapa en
precisión (cm) horizontal en el aproximado en 1990 Clase 1990 norma
horizontal nivel de imágenes de origen 1 Clase 2 NMAS
RMSEx y confianza del (cm)
RMSEy (cm) 95% (cm)
• Población.- “se refiere a todo el grupo o universo, los cuales tienen una variable común de
estudio” (Murray & Larry, 2009); además, se puede distinguir entre población infinita (un
• Margen de error (e).- Según Lara, Cerpa, Rodríguez y Núñez (2006), el margen de error es
• Nivel de confianza.- Según Triola (2004), “es un valor expresado en porcentaje que expresa
la certeza de que el dato se encuentre dentro del margen de error; las opciones más comunes
• Error.- en fotogrametría se define como “diferencia entre el valor medido y real; de los
cuales se considera valor real a las coordenadas de los puntos de control fotogramétr ico ”
(Quirós, 2014).
Donde:
• Error media cuadrático.- (RMSE).- La raíz cuadrada del promedio del conjunto de
diferencias cuadradas entre los valores de coordenadas del conjunto de datos y los valores
(ASPRS, 2014).
• Intervalo de confianza para la media.- Según Bianco y Martínez (2004), se define como “un
rango de valores, calculado en una muestra, que señala el valor del parámetro estimado con
𝑡𝑡 𝛼𝛼�2 . 𝑆𝑆 [ 17]
𝑥𝑥̅ ±
√ 𝑛𝑛
Donde:
79
- 𝑡𝑡 𝛼𝛼�2, representa el valor del área bajo la curva de la distribución t-student, a cierto nivel
de confianza.
• Sigma 0.- en la aerotriangulación indica la precisión general del ajuste; además, refleja la
errores en los GCP y errores sistemáticos; al mismo tiempo dicho valor es utilizado como
factor de escala en otros criterios estadísticos como desviación estándar, elipses de error, y
otras pruebas estadísticas; así mismo, el valor está influenciado por: calidad de las
posibilidad de compensar los efectos sistemáticos y observaciones GPS (Hinsken & Miller,
2002).
que permanece dentro de la elipse con una probabilidad del 95%; dicho criterio ayuda a
evaluar la precisión de los puntos y los centros de proyección después del ajuste de
• Error medio.- se define como el error posicional en un conjunto de valores para una
dimensión (x, y o z); obtenida al sumar todos los errores en una sola dimensión y luego
dividir por el número total de errores para esas dimensiones (ASPRS, 2014), tal y como
𝑛𝑛
1
𝑥𝑥̅ = � 𝑥𝑥 𝑖𝑖
(𝑛𝑛)
𝑖𝑖 =1 [ 18]
Donde:
coordenada de control.
• Desviación estándar.- es una medida de dispersión, que indica cuanto se puede alejar un
valor respecto a la media; representa el rango que puede tener un determinado valor
𝑛𝑛
1
𝑆𝑆𝑥𝑥 = � �(𝑥𝑥 𝑖𝑖 − 𝑥𝑥̅ )2
(𝑛𝑛 − 1)
𝑖𝑖=1 [ 19]
81
Donde:
CAPITULO III
3. METODOLOGÍA
el marco teórico, excluyendo el contrato y la toma de imágenes, debido a contar con las imágenes
3.1. Generalidades
fotogramétrico (GCP) en base a modificaciones de algunos autores, todo esto en programas como:
LPS, ORIMA y ZI/Imagine, a fin de analizar su influencia sobre la exactitud posicional de las
ortofotos, con estándares nacionales y la propuesta por ASPRS; por tanto, el desarrollo del
- Fase 3.- Generación de ortofotos con distribución homogénea de los GCP y métodos de
aerotriangulación.
En primer lugar, este estudio parte de una investigación previa del año 2015 titulada
Intergraph Erdas Imagine” cuyo autor es Barragán G., en cuya investigación se validan varias
cámara; en segundo lugar, la obtención de los GCP y puntos de exactitud posicional se realiza con
tecnología GNSS y considerando una distribución homogénea dentro de la zona de estudio; así
programa. Finalmente, el análisis sobre los productos fotogramétricos se realizará en base a los
estándares antes expuestos; así pues, cada paso se describe de forma general en la Figura 51.
Continúa
84
3.2.1. Imágenes
adquirieron 24 imágenes tomadas con cámara fotogramétrica matricial “Vexcel UltraCam XP”,
con las bandas R, G, B, IR, cuya fecha de toma y línea de vuelo a la que pertenecen, se detallan en
la tabla 13.
Tabla 13.
Fecha y línea de vuelo de la imágenes utilizadas en el proyecto
Imágenes Entregadas del Proyecto Mejía Rumiñahui 30 cm
Línea Número de foto: Fecha de toma
Desde Hasta Total
18 537 542 6 9 de Agosto 2012
19 829 834 6 19 de Septiembre 2012
20 930 935 6 19 de Septiembre 2012
21 1240 1245 6 19 de Enero 2013
24
Fuente: (Barragán, 2015).
En el entorno de trabajo de LPS, ORIMA y Z/I Imagen se decidió trabajar solo con las bandas
de exportación fue de TIFF con sus layers piramidales; además, los traslapos de las imágenes según
imágenes digitales no tienen una escala definida; por tanto, se realiza una recopilación de
información del tamaño de GSD en función de la escala y altura de vuelo con normativas
nacionales, IGN de Perú y el Instituto Geográfico Agustín Codazzi de Colombia, para la toma de
Tabla 14.
Relación del GSD y escala de la foto para imágenes digitales.
Altura de vuel o GSD [cm] Escala de la Autor
aproximada foto
5000 m 30 1: 50 000 (IGM, 2010)
-------- 28 1: 20 000 (IGM, 2006)
-------- 42 1: 30 000
-------- 30 1: 15 000 (IGN, 2011) de Perú
-------- 15 1: 10 000 (IGAC, 2007)
Fuente: (Barragán, 2015)
Además, según la ASPRS (2014), existe una relación entre escala de la foto, altura de vuelo,
GSD y la escala del mapa; en donde, para obtener un GSD de 30 cm la altura de vuelo es de
Tabla 15.
Relación entre escala de la foto y escala del mapa derivado.
Escalas comunes de fotografía (con cámara de formato de película de 9" y lente de
6")
Escala de la foto 1:3,600 1:7,200 1:14,400 1:28,800
Alturas de vuelo 550 m 1,100 2,200 m 4,400 m
m
Distancia aproximada de 7.5 cm 0.15 m 0.3 m 0.6 m
muestreo de tierra (GSD) de
escaneo
Escalas compatibles de mapas / ortoimágenes e intervalos de curva
GSD 7.5 cm 15 cm 30 cm 60 cm
Interval o de curva 30 cm 60 cm 1.2 m 2.4 m
Escala del mapa 1:600 1:1,200 1:2,400 1:4,800
Fuente: Modificado de (ASPRS, 2014)
Así mismo, existe una ecuación para calcular la altura de vuelo propuesta por Comer, Kinn,
Light y Mondello (1998), en la que relaciona el producto entre el tamaño del pixel en el terreno
(GSD), la distancia focal (f) de la cámara, sobre el tamaño del pixel en la cámara (CCD),
de vuelo; por tanto, en este estudio se siguen las recomendaciones de Barragán (2015) y la
a la cámara que se empleó en la toma de las imágenes, cuyas características se han visto en el
apartado 2.2.4 “Certificado de Calibración” y los valores más relevantes son: tamaño del pixel
utilizados por el IGM; sin embargo, el sentido de la rotación no se encuentra claramente; por tanto,
según Barragán (2015), el valor del punto principal de autocolimación pertenece a una rotación de
3.2.4. IMU/GPS
Se obtuvieron varios archivos, entregados por el IGM a la ESPE, llamados archivos APLANIX
(ver Figura 55), que son los metadatos de las imágenes y que contienen información como: versión
de uso, path de las imágenes descargadas, parámetros geodésicos, parámetros de la toma de vuelo,
sistema de referencia, coordenadas de los centros de las imágenes, ángulos de rotación, altura
depuración y así obtener las coordenadas de los centros con sus rotaciones; además, se debe
corregir el tipo de altura; que en base a una entrevista con el director de aerotriangulación del IGM,
Ing. Lenin Jaramillo manifiesta que “la altura expuesta en el archivo Aplanix es de carácter
ortométrico corregida con el modelo EGM96”; con respecto a la información de las coordenadas,
ángulos y línea de vuelo, son tal y como indica Barragán (2015) en su trabajo.
91
Considerando todos los parámetros anteriores, se procedió a seleccionar las coordenadas de los
centros, rotaciones de la cámara (IMU/GPS) y la altura ortométrica, siendo esta última sujeto a
depuración, debido a que utiliza el modelo EGM96 y según Tierra (2009), el modelo geoidal
mencionado puede llegar a tener errores en el Ecuador de hasta 4m al reducir de altura elipsoidal a
altura ortométrica; así pues, se obtuvo la ondulación geoidal con el programa F477 geoid, propuesto
por “The National Geospatila Intelligence (NGA)/NASA”, respecto del elipsoide de referencia
WGS84 (NGA/NASA, 2018), a fin de obtener la altura elipsoidal, proceso que es detallado en el
Anexo1; así pues, el archivo IMU/GPS con las coordenadas y los ángulos de rotación, se muestran
en la tabla 24 del apartado 4.1.1 del capítulo de resultados; posteriormente, al archivo IMU/GPS
bloque fotogramétrico con la orientación relativa generada por Barragán (2015), con el fin de
generar una ortofoto preliminar con parámetros por defecto (Molina, Redondo, Zamorano,
92
Hermosilla, & Caturla, 2010), que servirán como insumos para la navegación y posterior captura
El resultado de las monografías (ver tabla 25) y el bloque fotogramétrico con los puntos Von
Gruber se presenta en el apartado 4.1.2, cuya ortofoto preliminar se muestra en la Figura 57.
GCP; así como, el uso que tendrán los mismos. Los puntos fueron planificados en base las teorías
3.3.1. Planificación
detalladas en el apartado 2.2.3 “Control geodésico o control terrestre” y las teorías expresadas en
Las distribuciones de los GCP fueron elegidas en base a las teorías de los siguientes autores:
según Ackermann (1990), al ubicar los puntos en las esquinas del proyecto de una manera múltip le
y según el IGN de Perú (2011), con puntos distribuidos homogéneamente cada tercer modelo
perpendicular a las líneas de vuelo y en zonas de triple traslapo; (ver Figura 58), a la que se
Segundo, considerando lo expuesto por Ackermann (1987), se distribuye los puntos al contorno
del proyecto; además de puntos en el interior en zonas de triple traslapo, expuesto por el
Departamento de Transporte de Colorado (2015) y Milkhail et al. (2001), dando un total de 13 GCP
A pesar de seguir las consideraciones del apartado 2.5.1 “Ubicación de los puntos de control
inexistencia del objeto planificado, lugares con accesos privados y temporalidad muy larga entre
las líneas de vuelo 20 y 21; por tanto, las muestras de campo de los GCP obtenidos se presentan en
el apartado 4.2.1.
La navegación a los puntos se realizó mediante un equipo GPS navegador marca Magellan
ArcGis, el cual contiene la ortofoto preliminar, además de la cartografía básica de la zona de acceso
libre obtenida del Geoportal del IGM; este proceso es detallado por Kirby (2009), en su trabajo
“Metodología para navegación en Arcgis a través del uso del GPS”; por tanto, en la Figura 60 se
Así pues, se realizaron ocho salidas de campo, para obtener un total de 48 puntos (ver tabla 16),
rastreados con un equipo GNSS marca Trimble R3 de una frecuencia y las mediciones al ápice de
la muesca de la antena.
Tabla 16.
Fechas y cantidad de puntos obtenidos en campo.
No Fechas Cantidad de puntos obtenidos
1 20 de abril de 2016 3
2 22 de abril de 2016 5
3 30 de julio de 2016 4
4 31 de julio de 2016 2
5 13 de marzo de 2017 1
6 14 de marzo de 2017 11
7 15 de marzo de 2017 8
8 25 de septiembre de 2017 7
Total de puntos 41
• Método de posicionamiento estático rápido (Fast Static), con una base (Estación REGME –
EPEC) y un móvil (Trimble R3). La máxima distancia a la base fue de 7950 m y la distancia
mínima de 40 m.
• Puntos foto identificables (ver Figura 61 y Figura 62) en varias zonas como: cruces de vías,
3.3.3. Procesamiento
Los puntos fueron procesados con el software comercial Trimble Business Center, con el uso
de candado de procesamiento, considerando la constelación GPS, sin efemérides finales siguie ndo
97
las recomendaciones de Portilla, Tierra y Godoy (2016). La base fue fijada con las coordenadas de
la base EPEC, la cual garantiza la consistencia de las coordenadas a partir del 17 de abril de 2016,
ajustadas al ITRF08, época de referencia 2016.4; todo esto debido al evento sísmico que sufrió el
Ecuador el 15 de abril de 2016; sin embargo, la vigencia de las coordenadas variará en función de
Además, los GCP deben cumplir con ciertas características como: un PDOP menor a 4, precisión
mayor a la cartografía a generar (ver apartado 2.6.4); por tanto, al utilizar el criterio del IGM, se
tiene 0.1 por el módulo de la escala (0.1mm * 5000=500mm), dando un valor de 50cm al evaluar
𝑃𝑃𝑃𝑃
𝑃𝑃𝐺𝐺𝐺𝐺𝐺𝐺 =
4
[ 21]
Donde:
De la misma forma, la norma Positional Accuracy Standards for Digital Geospatial Data de la
ASPRS de la tabla 12 el valor del RMSE del mapa es de 63.6cm (tamaño del GSD de 30 dentro
98
del intervalo de la tabla 12) y al aplicar la ecuación 11 se tiene, de la misma norma resulta 16cm
1 45 [𝑐𝑐𝑐𝑐] [11]
𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑋𝑋 𝑜𝑜 𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑦𝑦 = ∗ [𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑋𝑋−𝑌𝑌 (𝑀𝑀𝑀𝑀𝑀𝑀𝑀𝑀)] = = 11.25 [𝑐𝑐𝑐𝑐]
4 4
Donde:
• [𝑅𝑅𝑀𝑀𝑀𝑀𝑀𝑀𝑋𝑋−𝑌𝑌 (𝑀𝑀𝑀𝑀𝑀𝑀𝑀𝑀)], representa el valor del RMSE del mapa, de la tabla 12 con un GSD de
Entonces, se tiene que el procesamiento de los datos se descarta 5 puntos, debido a un PDOP
Tabla 17.
Informe de procesamiento de puntos con inconvenientes.
ID Norte [m] Este [m] Altura Tiempo P. Horizontal P. Vertical PDOP
Elip. [minut.] [m] [m]
[m]
X100 9963877.71 784411.76 2524.69 25 0.59 1.55 15.45
I045 9962449.41 781145.81 2537.25 25 0.01 0.01 11.41
G026 9962828.05 782527.12 2522.94 20 0.02 0.03 8.02
K051 9966866.13 777719.97 2801.54 20 0.01 0.02 4.32
K050 9962416.33 785065.98 2544.39 20 0.02 0.03 4.16
Nota: P (Precisión), Elip (Elipsoidal).
99
En total se obtuvieron 36 puntos con las condiciones descritas de GNSS y foto identificab les
aceptables, a las que se les añadió los puntos de control obtenidos por Barragán (2015) (ver tabla
25); sin embargo, se descartaron los puntos GCPG y GCPI al no contar con información del PDOP.
posicional fueron de 45 puntos, cuyos detalles de precisiones están bajo los límites de la norma
más estricta (Positional Accuracy Standards for Digital Geospatial Data=16cm) se indican en la
tabla 26 y su localización en la Figura 128 respectivamente del apartado 4.2.2 de los resultados;
además, la documentación de cada punto está representada en una monografía, cuyas características
mismo apartado.
matemáticos y funcionalidades son detalladas en el apartado 2.4.1 “Erdas-LPS”; por tanto, en LPS
se realiza la creación del proyecto, importación de los centros, las diferentes orientaciones, la
aerotriangulación por diferentes métodos, generación de MDS, las ortofotos y los diferentes
fotogramétrico (GCP).
100
En el programa fotogramétrico ERDAS Imagine 2011 se debe elegir el módulo LPS para la
creación del proyecto fotogramétrico y escoger en “Model Setup”, el sensor de cámara digita l;
posteriormente, ingresar el sistema de coordenadas WGS84 proyección UTM zona 17 Sur y una
altura de vuelo de 5000m del apartado 3.2.2 “Escala y altura de vuelo” y ejemplificado en la Figura
63. La explicación detallada de este proceso, así como la orientación interior se manifiesta a detalle
en el Anexo 3.
Posteriormente, se debe insertar el IMU/GPS del apartado 4.1.1 en formato txt (ver Figura 65)
Culminados los pasos anteriores, se procede a direccionar todas imágenes del archivo
IMU/GPS, que contienen las tres bandas (RGB) y cuentan con sus layers piramidales, teniendo un
entorno de trabajo como la Figura 66, donde el color verde representa los pasos terminados.
Como se indicó en el apartado 2.2.5, la orientación interna de la cámara digital viene dada por
el certificado de calibración; por tanto, se incluyen los valores del apartado 3.2.3 “Información de
la cámara” pertenecientes a un tamaño de pixel de 6µm y se aplica para todas las imágenes activas
en el menú “Interior Orientation” de la ventana “Digital Camera Frame Editor”; así pues, en la
ventana general de LPS en la columna “Int” se vuelve de color verde, indicando que el proceso ha
Como se indicó en el apartado 2.2.6.1, la orientación relativa se relaciona cada objeto de una
imagen con su homólogo en la imagen posterior, usando la distribución de la teoría de los puntos
de Von Gruber visto en el apartado 2.2.6.2. Este proceso en LPS Core es de forma monoscópica
(vista 2D) y se deben tener en cuenta las distorsiones de las imágenes tales como: distorsión por
parámetros de orientación, por el relieve, por la inclinación de la imagen, entre otras; detalladas en
Por tanto, se procede a importar los puntos de la orientación relativa de Barragán (2015), que
fueron pinchados manualmente por el autor y con ayuda de la herramienta Point Measurement,
siguiendo la conceptualización del apartado 2.2.6.2 “Puntos de Von Gruber (ver Figura 68),
Como se indicó en el apartado 2.2.6.3 “Orientación absoluta”, los puntos de control escalarán
el modelo con respecto al terreno; por tanto, se procede a pinchar los 13 puntos de control (GCP)
Tabla 18.
Identificadores de los GCP en la Distribución 13.
NO 13 GCP
1 A001 D019 I044
2 A004 D020 I046
3 A006 G030 I048
4 A008 H036
5 A010 H040
105
La distribución de los puntos de control pinchados vistos desde la ventana general de LPS se
indica la Figura 69; además, vale destacar que el valor de altura de los GCP en el los programas
LPS y Orima son de carácter elipsoidal con el fin analizar la influencia directa de la
aerotriangulación sobre los productos; debido a que la transformación de altura elipsoidal a altura
ortométrica puede tener errores de hasta 4m dependiendo del modelo geoidal usado (EGM96 o
Una vez pinchados los GCP e importados los puntos de Von Gruber, se procede a la
densificación de puntos (Tie Points), con la ayuda de la herramienta Automatic Tie Point
estándar del punto principal (X=Y=0.33, por defecto) y las opciones avanzadas de eliminac ió n
robusta de los errores groseros “Advanced Robust Checking”, sin parámetros adicionales por
Al correr el programa (Run) se genera una ventana con el valor del RMS final del ajuste presente
en el capítulo de resultados 4.3.1.1 (ver Figura 130); el cual se refiere al “Método de haces de luz
en ajuste de bloques”, sin matriz de pesos en su compensación, tal y como se definió en el apartado
2.3.5; además, se ejemplifica el reporte del ajuste de la aerotriangulación en la Figura 131. Una
107
vez aceptada la aerotriangulación aparece en color verde la columna “Ext” en la ventana de LPS
(Figura 72).
superficie del terreno con sus elementos y como tal, en este proyecto representa el valor de altura
de carácter elipsoidal.
Al crear un MDS en el entorno de Erdas – LPS se utiliza la opción “Start DTM Extraction –
Classic ATE” donde se configura “Advanced Properties” al activar los puntos de control y los
Una vez aceptados los cambios de Advanced Properties, en la ventana DTM Extraction, se
deben ingresar los valores de Cell Size que pueden tener varias consideraciones, entre ellas: criterio
de la cartografía 1:5000, con un intervalo de curva de 5 metros o el tamaño del pixel como en el
caso de Barragán (2015), donde se introduce un valor de 0.33 m; sin embargo, estudios como el de
Mora (2013), indican que al colocar valores más pequeños de celda, crean grandes variaciones en
las curvas de nivel; por tanto, se procede a insertar un valor de 5m (ver Figura 74), del criterio de
la cartografía a obtener.
Este procedimiento se realiza solo una vez y el programa genera un MDS por cada par de
imágenes de la misma línea de vuelo; por tanto, se generan cinco MDS por cada línea de vuelo,
dando un total de 20 archivos, este proceso puede tardar de entre 3 a 5 horas, dependiendo de la
Una vez generado los MDS se procede a rectificar cada imagen, con la herramienta “Ortho
Resampling”, considerando los estudios de Serra, Moré y Pons (2010); Fallas (2011); donde,
indican que el método de resampleo bilineal es el más recomendable al trabajar con fotografías
aéreas digitales y el resto de valores son por defecto (ver Figura 75).
En la pestaña “General” de la ventana anterior, se debe seleccionar cada MDS creado, del
apartado 3.4.1.7, correspondiente a la imagen del traslapo; además, insertar el tamaño del pixel en
cada una, cuyo valor corresponde a la teoría del apartado 2.1.4 “Resolución imagen-sensor”; donde
110
el mínimo objeto a distinguir en cartografía 1:5000 es de 1m; por tanto, se necesita de 3 a 4 píxeles
para visualizar dicho objeto; por tanto, se inserta un valor de 30cm como tamaño de celda (ver
Figura 76).
Dicho proceso se realiza para las 24 imágenes y al culminar con todas, la columna Ortho aparece
Como se indicó en el apartado 2.3.8.4, un mosaico es la unión entre imágenes con cierto
Erdas, tanto en el MDS como en la Ortofoto, la elaboración del mosaico es detallado en el Anexo
con la herramienta “Display Add Images” (ver Figura 78), se procede a realizar el mosaico con
parámetros por defecto, pues es de interés conocer sus valores sin modificaciones digitales.
apartado 2.1.4, a fin de dar continuidad a los elementos de las ortofotos; por tanto, se utiliza la
herramienta “Color Corrections” y el método “Automatic Color Balancing” (ver Figura 79).
112
Además, se utiliza el histograma de unión por áreas de traslapo (ver Figura 80).
En las líneas de unión del mosaico se utilizó la herramienta Seamline Generation Options, que
indica la vista más perpendicular a la cámara (Most Nadir); posteriormente, se procede a generar
el mosaico, ilustrado en la Figura 81 (para mayor detalle de este proceso se puede consultar el
Anexo 5). El resultado del mosaico del MDS se presenta en el apartado 4.4.1 y el resultado del
necesario seguir los mismos pasos desde el apartado 3.4.1.1 “Creación del bloque en LPS” hasta
Con los puntos de control pinchados y los puntos de Von Gruber importados, se procede a
generar los Tie points con los parámetros descritos en el apartado 3.4.1.6, pero la compensación se
realizará considerando la matriz de “Ebner’s” (ver Figura 82), pues es la ideal de los modelos
Al ejecutar el programa (Run) se genera un ventana con el valor del RMS final de la
aerotriangulación presentado en el apartado 4.3.1.2 (ver Figura 133) y al aceptar los resultados, se
actualizan las coordenadas en los puntos tie y aparece de color verde la columna “Ext”, que indica
Así mismo, al realizar los MDS con 13 GCP y considerando la matriz de pesos (modelo de
Ebner’s) se siguieron los mismos pasos del apartado 3.4.1.7, que serán utilizados en la elaboración
Al realizar las 24 ortofotos de las 24 imágenes, con 13 GCP y considerando la matriz de pesos
(modelo de Ebner’s), se siguieron los mismos pasos del apartado 3.4.1.8, cuyas ortofotos serán
En la realización del mosaico del MDS con 13 GCP y considerando el modelo de Ebner’s se
siguieron los mismos pasos del apartado 3.4.1.9; además, sus resultados se muestran en el apartado
4.4.2 “Resultado del mosaico MDS con matriz de pesos y 13 GCP en LPS”; mientras que, el
resultado del ortomosaico se presenta en el apartado 4.5.2 “Resultado del ortomosaico con matriz
los puntos de control, es necesario seguir los mismos pasos del apartado 3.4.1.1 “Creación del
116
Tabla 19.
Identificadores de los GCP en la Distribución 19.
NO 19 GCP
1 A001 B012 D020 I044
2 A004 C014 G025 I046
3 A005 C015 G030 I047
4 A010 C017 G035 I048
5 A011 D019 H042
La distribución de los 19 GCP a lo largo del bloque en el entorno de LPS se ilustra la Figura 84.
En la realización de la aerotriangulación sin matriz de pesos, se debe realizar los mismos pasos
del apartado 3.4.1.6 “Configuración de la aerotriangulación sin matriz de pesos en LPS” hasta
ajustar el proyecto, cuyos resultados se indican en del apartado 4.3.1.3; así mismo, la columna
En la realización de los 20 MDS con 19 GCP y sin considerar la matriz de pesos, se siguiero n
los mismos pasos del apartado 3.4.1.7; además sus resultados se utilizarán en la elaboración del
Al realizar las 24 ortofotos de las 24 imágenes, con 19 GCP y sin considerar la matriz de pesos,
se siguieron los mismos pasos del apartado 3.4.1.8; además, sus resultados se utilizarán en la
En la realización del mosaico del MDS con 19 GCP y sin considerar la matriz de pesos se
siguieron los mismos pasos del apartado 3.4.1.9, cuyos resultados se encuentran el apartado 4.4.3
“Resultado del mosaico MDS sin matriz de pesos y 19 GCP en LPS”; mientras que el resultado del
ortomosaico se presenta en el apartado 4.5.3 “Resultado del ortomosaico sin matriz de pesos y 19
GCP en LPS”.
necesario seguir los mismos pasos del apartado 3.4.1.1 “Creación del bloque”, hasta el apartado
3.4.1.4 “Configuración de la orientación relativa”; así mismo, considerar el mismo pinchado de los
19 puntos de control del apartado 3.4.3.1 “Configuración de la orientación absoluta con 19 GCP”,
Al generar los tie points se siguieron los mismos pasos de apartado 3.4.2.1 “Configuración de
la aerotriangulación con matriz de pesos”, con una distribución de 19 GCP, cuyos resultados se
encuentran en el apartado 4.3.1.4. En la ventana general de LPS la casilla “Ext”, muestra de color
Ebner’s), se siguieron los mismos pasos del apartado 3.4.1.7; además, sus resultados se utilizará n
Al realizar las 24 ortofotos de las 24 imágenes del proyecto, con 19 GCP y considerando la
matriz de pesos (Modelo Ebner’s), se siguieron los mismos pasos del apartado 3.4.1.8; además, las
En la realización del mosaico del MDS con 19 GCP y sin considerar la matriz de pesos (modelo
de Ebner’s) se siguieron los mismos pasos del apartado 3.4.1.9, cuyos resultados se encuentran el
apartado 4.4.4 “Resultado del mosaico MDS con matriz de pesos y 19 GCP en LPS”; mientras que,
120
el resultado del ortomosaico se presenta en el apartado 4.5.4 “Resultado del ortomosaico con matriz
depende de LPS-Core en cuanto a la creación del proyecto y la orientación interna; por tanto, se
realizan los mismos pasos desde el apartado 3.4.1.1 “Creación del bloque” hasta el apartado 3.4.1.3
elipsoidal; así mismo, en la aerotriangulación con Orima es la denomina ajuste en bloque y ajuste
en líneas de vuelo.
Al trabajar en el entorno de Orima es necesario conocer todas las funcionalidades del programa,
La versión 11.00 de Orima se abre desde Erdas-LPS; donde, se debe escoger el tipo de licencia
para trabajar con cámara digital y apoyo GPS denominado Orima DP-TE/GPS (ver Figura 87);
En la ventana siguiente de Orima se encuentra una opción que permite insertar el nombre del
usuario, unidades y otros parámetros; al aceptar dicha preferencia, se abren dos ventanas de trabajo,
(ver Figura 88), que se usa en los procesos de cálculo, inserción de parámetros, mediciones, entre
Orima permite la visualización de los puntos que están fuera de un rango o umbral previame nte
insertados; por tanto, se insertan primero los parámetros antes de las mediciones, con el fin de tener
Los primero parámetros a ser insertados son: “Accuracy Requirements” en la que se introduce
el valor de 2 en Relative Orientation, que se refiere al tamaño del pixel de la cámara de 6µm entre
tres (Sigma0 = 1/3 tamaño del pixel) (Hinsken, 2010), el valor de Absolute Orientation se calcula
considerando la calidad del producto a obtener (Cartografía 1:5000 con norma ISO 0.2mm) y el
micro procesos (orientación interna, relativa y absoluta); es decir, dividir 1m entre todos los
procesos que intervienen, dando un valor aproximado de 0.10m al el resto de valores por defecto
(ver Figura 89), para mayores detalles consultar Anexo 7 (Tema: Pasos introductorios).
(secuencialmente) y lectura de modelos (al finalizar las mediciones) (ver Figura 90).
Los parámetros de cálculo en la orientación absoluta tiene cierta similitud con la relativa, pero
difiere en cuanto a los puntos de control siguiendo las siguientes configuraciones: puntos de enlace
proyectados independientemente (Not Measured) y proyectar todos los puntos de control (All) (ver
Figura 91).
Orima utiliza la medición monoscópica con plantillas en el pinchado de puntos Tie o de amarre;
donde, se conduce una cruz a la posición escrita en la plantilla; ahora, el operador puede escoger
que objeto o detalle es adecuado para pinchar; a todo este proceso se le conoce como orientació n
relativa. Las coordenadas de la plantilla son coordenadas locales de la imagen cuya distribución y
cantidad dependen del operador; se puede utilizar el criterio de los puntos de Von Gruber o el
Este criterio de plantillas permite una repetición uniforme del mismo proceso en todos los
modelos estereoscópicos; así pues, en este proyecto se generaron dos plantillas llamadas: Left1_tes
y Right1_tes, con tres puntos por plantilla y una codificación en múltiplos de 10; aquellos IDs de
número impar, pertenecen a la plantilla izquierda y con número par pertenecen a la plantilla
derecha. Las coordenadas (X, Y) de la plantilla se detallan en la tabla 20, al mismo tiempo la
Estas plantillas son importadas directamente al directorio de Orima en Project Files (Scheme) en
formato *sch.
Tabla 20.
Disposición de las plantillas con los puntos de von Gruber.
Archivo Colum. 1 Descripción Colum. 2 Descripción Colum. 3 Descripción Colum. 4 Desc.
Left1_tes 10 ID de los 0 Coordenada 100 Coordenada -100.0 Coordenad
30 puntos 0 del eje x en 0 del eje y en -100.0 a del eje Z
50 0 mm -100 mm -100.0 en mm
Right1_Tes 20 80 0.0 -100.0
60 80 -100.0 -100.0
40 80 100.0 -100.0
125
Una vez creada la plantilla e importada dentro de Project Files, es necesario codificar los puntos
de Von Gruber para su medición en Orima, el cual como ejemplo de un punto codificado es:
i_18_0537_10 (Sufijo i, 18 línea de vuelo, 0537 Imagen, 10 ID del punto de von Gruber), (ver
Figura 92); los detalles de este proceso se detallan en el Anexo 7 (Tema: Orientación Relativa).
ventana Orima Orientation Management para transferir los puntos y de la ventana Leica
Este proceso se realiza en todas las imágenes hasta que el par estereoscópico cumpla con los
parámetros ingresados en el apartado 3.5.1.2. El color verde representa los puntos dentro de los
umbrales del proyecto y el color rosado el punto con el máximo de los valores (Ver Figura 94).
En el entorno de Orima se requiere importar las coordenadas (Este, Norte y Altura) y los IDs de
los puntos de control de la tabla 18, en formato txt, con la herramienta Import Control Point;
posteriormente, se debe cambiar el tipo de vista de Split a Mono Left, que indica una vista por
En el pinchado de los puntos de control, el programa Orima localiza el GCP automáticame nte
de una manera aproximada, con la herramienta Measure Multiple y el operador debe pinchar uno
a uno, los puntos de control, en cada imagen con la ayuda visual de las monografías (Anexo 2) (ver
Cabe mencionar que el programa también cuenta con el pinchado automático de los puntos, con
pruebas que se realizaron en su pinchado generaban mucho error (ver Figura 97).
128
realizar un primer ajuste, que disminuirá la propagación de errores en el pinchado de los puntos
Tie (von Gruber) o de los puntos de control (GCP); además, este proceso da una apreciación
preliminar del estado de los puntos de control, que son los puntos que mayormente aportan a la
aerotriangulación.
Adjustment Program for Aerial Triangulation), cuyos parámetros se insertan en la ventana Input
Parameters.
(Recomendación del Dr. L. Hinsken creador de Orima para observaciones medidas manualmente),
129
Sigma a priori de 2 (1/3 del tamaño del pixel de 6µm), corregir las coordenadas por refracción
en metros y el resto de valores por defecto (ver Figura 98); así mismo, el detalle del contenido de
Al realizar el primer ajuste, el programa crea un reporte indicando cantidades de: puntos de
control, puntos tie, imágenes, observaciones, cámaras, archivos IMU/GPS, entre otros; los valores
más relevantes son de Im. Blunders que indica el número de puntos tanto de control como Tie que
tienen algún error, los valores de Add. Blunders que indica los posibles errores de: cámara, centros
de proyección, parámetros GPS, entre otros; y el valor de Sigma 0 (ver apartado 2.7.1) que es un
valor aproximado al insertado en Sigma a priori = 2 (ver Figura 99 recuadros de color rojo)
130
En este proyecto se corrigieron los valores de Im. Blunders hasta llevar al menor valor posible
siguiendo las consideraciones del Dr. L Hinsken; la corrección de dichos puntos se detallan en el
A continuación, se debe realizar la densificación de puntos tie, en la que se utiliza la herramie nta
APM (Automatic Point Measurement) de Orima (ejemplificado en los cuadros de color rojo), en la
cual se considera la cantidad de puntos por imagen (81), debido a varias pruebas realizadas con
distribuciones de 25, 49, 81 y 121, siendo el más conveniente el de 81 (A); sin embargo, el usuario
puede crear una distribución de puntos tie propia (ver Anexo 6 Tema: APM). El siguie nte
parámetro considerado fue la estrategia del terreno, que en este caso para todo el proyecto es de
tipo: Alto Urbano (High Urban) (B); es decir, que para cada imagen buscará 81 puntos tie según el
131
algoritmo para terreno urbano, evitando las sombras (Avoid Shadow)(C) (ver Figura 100). A esta
La generación de los puntos tie puede tardar varios minutos u horas dependiendo de la capacidad
del computador; al finalizar, se debe ejecutar el programa CAP-A, cambiando el valor de las
iteraciones por 50 en Input Parameters, que es la recomendación del autor de Orima con medicio nes
la cantidad de puntos y el Sigma0, visto en el apartado 3.5.1.5; así mismo, se debe tener en cuenta
La edición es importante debido a la existencia de errores como: puntos fuera del área de triple
traslapo (ver Figura 101), puntos no pinchados en todas las imágenes posibles, mala correlación
132
entre imágenes, puntos en zonas erróneas, objetos que no se encuentran en todas las imágenes,
entre otras.
Estos errores se pueden visualizar utilizando las elipses de error (ver apartado 2.7.1) y son de
de mayor cantidad de datos en su compensación; por tanto, al terminar una edición de puntos se
debe seguir ejecutando CAP-A, hasta que el valor de Im. Blunder disminuya o el operador
el apartado 4.3.2.1. Culminados estos pasos, la orientación relativa de LPS se coloca en verde; por
el entorno de Erdas – LPS y como tal se siguen los mismos pasos del apartado 3.4.1.7, cuyos MDS
en bloque, se trabaja en el entorno de Erdas – LPS y como tal se siguen los mismos pasos del
siguieron los mismos pasos del apartado 3.4.1.9, cuyos resultados se encuentran el apartado 4.4.5
“Resultado del mosaico MDS con ajuste en bloque y 13 GCP en Orima”; mientras que el resultado
134
del ortomosaico se presenta en el apartado 4.5.5 “Resultado del ortomosaico con ajuste en bloque
y 13 GCP en Orima”
En la realización de un nuevo ajuste en el entorno de Orima, realizaron los mismos pasos desde
de Orima, con un valor de 81 por imagen de 81 (A), la estrategia del terreno es por líneas de vuelo,
siendo las líneas 18 y 19 de tipo: Montaña (Mountain) (B); es decir, que para cada imagen de la
línea 18 y 19 (12 imágenes) (C) buscará 81 puntos tie siguiendo el algoritmo para terreno
montañoso, y las líneas de vuelo 20 y 21 con estrategia de terreno urbano (High Urban); además,
se mantiene la configuración de evitar las sombras (Avoid Shadow) (D) (ver Figura 103). A este
Cabe mencionar que este proceso se realiza primero para la línea 18 y 19, y al finalizar se
continua con las líneas 20 y 21; posteriormente, se debe cambiar el valor de las iteraciones por 50
y ejecutar el programa CAP-A, con ello se genera el reporte típico con los Im Blunders y el Sigma
0 (ver Figura 104); por tanto, es necesario modificar los puntos mencionados y continuar
Después de las ediciones, visualizar el estado de los puntos, proceso detallado en Anexo 7
apartado 4.3.2.2. Con dicho ajuste el proceso de la orientación relativa en LPS se culmina y la
3.5.2.2. Generación de los MDS con ajuste de líneas de vuelo y con 13 GCP
La realización de los 20 MDS con 13 GCP y considerando un ajuste de líneas de vuelo se siguen
los mismos pasos del apartado 3.4.1.7, cuyos MDS se utilizaron en la elaboración del mosaico del
apartado 3.5.2.4.
en líneas de vuelo, se siguieron los mismos pasos del apartado 3.4.1.8, cuyas ortofotos se utilizaro n
La realización del mosaico del MDS con 13 GCP y considerando un ajuste en líneas de vuelo
se siguieron los mismos pasos del apartado 3.4.1.9, cuyos resultados se encuentran el apartado
4.4.6 “Resultado del mosaico MDS con ajuste en líneas de vuelo y 13 GCP en Orima”; mientras
que el resultado del ortomosaico se presenta en el apartado 4.5.6 “Resultado del ortomosaico con
En la realización de un nuevo ajuste en el entorno de Orima, se siguen los mismos pasos del
apartado 3.5.1.1 “Pasos introductorios de Orima” hasta el apartado 3.5.1.3 “Orientación relativa en
Orima” y el mismo tipo de aerotriangulación del apartado 3.5.1.6, pero cambiando la orientació n
largo del bloque, como en la Figura 106, siguiendo las mismas recomendaciones del apartado
3.5.1.4.
138
Al realizar el nuevo ajuste con los 19 GCP, se realizaron mismos pasos desde el apartado 3.5.1.5
pero con la distribución de 19 GCP y considerando la edición de puntos, hasta ajustar el proyecto,
cuyos resultados se presentan en el apartado 4.3.2.3; y como tal se puede continuar con la
En la realización de los 20 MDS con 19 GCP y un ajuste en bloque de Orima, se siguieron los
mismos pasos del apartado 3.4.1.7, cuyos MDS se utilizaron en la elaboración del mosaico del
apartado 3.5.3.5.
139
Orima, se siguen los mismos pasos del apartado 3.4.1.8, cuyas ortofotos se utilizaron en el mosaico
La realización del mosaico de MDS con 19 GCP y un ajuste en bloque, se siguieron los mismos
pasos del apartado 3.4.1.9, cuyos resultados se encuentran el apartado 4.4.7 “Resultado del mosaico
MDS con ajuste en bloque y 19 GCP en Orima”; mientras que el resultado del ortomosaico se
presenta en el apartado 4.5.7 “Resultado del ortomosaico con ajuste en bloque y 19 GCP en Orima”.
En este nuevo ajuste se realizaron los mismos pasos del apartado 3.5.1.1 “Pasos introducto r ios
de Orima” hasta el apartado 3.5.1.3 “Orientación relativa en Orima” y el mismo marcado de los 19
puntos de control del apartado 3.5.3.1 “Configuración de la orientación absoluta con 19 GCP en
En la generación de los tie points se siguen los mismos pasos del apartado 3.5.2.1, para una
con ajuste de líneas de vuelo en Orima.”, hasta la edición de los puntos; cuyos resultados se
140
siguieron los mismos pasos del apartado 3.4.1.7, cuyos MDS se utilizaron en la elaboración del
La realización de las 24 ortofotos de las 24 imágenes, con 19 GCP y un ajuste en líneas de vuelo,
se siguieron los mismos pasos del apartado 3.4.1.8, cuyas ortofotos se utilizaron en el mosaico del
apartado 3.5.4.4.
En la realización del mosaico de MDS con 19 GCP y un ajuste en líneas de vuelo en Orima, se
siguieron los mismos pasos del apartado 3.4.1.9, cuyos resultados se encuentran el apartado 4.4.8
“Resultado del mosaico MDS con ajuste en líneas de vuelo y 19GCP en Orima”; mientras que el
resultado del ortomosaico se presenta en el apartado 4.5.8 “Resultado del ortomosaico con ajuste
Por último se realiza una aerotriangulación en el programa fotogramétrico Z/I Imagine, el cual
realiza un ajuste de aerotriangulación por el modelo de haces de luz, cuyas ecuaciones utilizan el
(GCP); sin embargo, al realizar varias pruebas en el ajuste, se decidió no utilizar ciertos puntos,
pues con ellos se incrementaba el RMSE y no pasaba el ajuste; por tanto, se determinaron 13 puntos
de control ideales para el pinchado en 3D, cuya distribución se ilustra en la Figura 107, al considerar
Automatic Triangulation versión 13.00 (2013), en el cual se configura el nombre del proyecto, el
142
tipo (Aerial Photography), sistema de referencia (UTM WGS84 Zona 17 Sur), unidades de trabajo
Antes de insertar los datos de precisiones en Z/I Imagine, cabe señalar lo siguiente: el Ecuador
presenta diversidad climática y morfológica; algo muy beneficioso desde el punto de vista de la
biodiversidad; sin embargo, para fines fotogramétricos la presencia de nubosidad, brumas, clima
muy variable en la toma, topografía muy irregular en zonas cortas, entre otros; ocasiona errores en
la adquisición de las imágenes, llevando incluso a una nueva toma, es así que las precisiones se
deben adaptar a las realidades del medio, a fin de tener amplitud al cumplir las características
técnicas del proyecto, es por este motivo que se trabaja con los valores máximos, de las precisiones
Considerando lo anterior, al generar cartografía 1:5000 con la norma IGM (0.3mm), la precisión
del producto es de 1.5m (producto entre el módulo de la escala y el factor en metros); además, la
precisión altimétrica es ¼ del intervalo de curva (5m/4=1.25m), dicho valor se inserta para todos
los procesos, pues es el valor máximo que se debería cumplir el producto; sin embargo, los cálculos
para generar cartografía a escala 1:5000 considerando un ajuste ideal; es decir, sin nubosidad,
buenos centros de las imágenes, sin atemporalidad entre líneas de vuelo, topografía homogénea,
entre otros factores; se realiza de la siguiente manera: 5000 * 0.3mm = 1500mm; obteniendo 1.5m
aerotriangulación; en altimetría se tiene que es un cuarto del intervalo de curva obteniendo un valor
de 1.25m.
143
Además, los valores de Sigma deberían ser de ½ del tamaño del pixel de la cámara (Intergraph,
2013) (6um/2=3um); sin embargo, se utiliza los valores máximos de 10um en todos las
apartado 3.2.3 y los datos de IMU/GPS, obtenidos del apartado 3.2.4, incluyendo en el archivo txt
la ubicación de la carpeta de las imágenes (ver Figura 109), este proceso es detallado en el Anexo
8.
144
importación de los valores del certificado de calibración de la cámara; tal como, se definió en el
apartado 2.2.5 “Orientación interna”; así pues, se puede continuar con las siguientes orientacio nes.
orientación absoluta; que en el caso del programa Z/I Imagine, la orientación relativa no se realiza
con todos los puntos de von Gruber, debido a que las imágenes se encuentran orientadas de manera
aproximada con ayudan de los centros de las mismas; por esta razón, se aplica el concepto de los
puntos de von Gruber solo en zonas donde la correlación entre imágenes es débil y no se han
La orientación absoluta, busca orientar el modelo con respecto al terreno, con ayuda de los
centros IMU/GPS, la correlación entre imágenes y los puntos de control fotogramétrico; por tanto,
se importa los 13 puntos de control de la tabla 21, para realizar la orientación absoluta y a partir de
Tabla 21.
Identificadores de los GCP en el ajuste Z/I Imagine.
NO 13 GCP
1 A006 E023 H040
2 A009 G025 I046
3 A010 G029 I048
4 C014 G030
5 D020 G034
Los puntos de control (GCP), se importan con la opción Control Points del menú ISAT, en
formato txt y siguen el orden de: ID, Este, Norte y Altura ortométrica, corregida con el modelo
geoidal del IGM, pues es un modelo para el Ecuador; posteriormente, se debe pinchar cada GCP
en vista estereoscópica (ver Figura 110), hasta tener los 13 GCP marcados.
146
Con los puntos de control marcados, se procede a realizar la densificación de puntos de paso,
con una grilla de 3x3 en cada imagen, con ayuda de la herramienta Thinning And Weak Area
como: mala correlación, atemporalidad entre líneas de vuelo (línea 20 y 21), desplazamientos en
147
zonas limítrofes del proyecto, entre otros; su valor alto es representado de color rojo en el programa
Z/I Imagine (ver Figura 112); además, los puntos en color verde y con símbolo de asterisco son
puntos de amarre bien correlacionados y los puntos rojos con triángulo son puntos de control.
el apartado 2.3.5 “Método de haces de luz en ajuste de bloques” en el programa Z/I Imagine, cuyo
los puntos densificados pasan a tener la simbología de triángulos de color rojo, ilustrados en la
Figura 113.
148
pues en este programa no se realiza un MDS por cada par estereoscópico como anteriormente se
ha visto. Las principales configuraciones se refieren a la extención del MDS, la altura media del
terreno (2500m), el formato de salida (.dtm) y el tamaño del pixel, que el programa genera
automáticamente de 9m; sin embargo, se puede considerar el intervalo de curva para generar
cartografía 1:5000 (5m). Una vez configurados los parámetros se procesa el mosaico de MDS y su
Al presentar el mosaico de MDS de los resultados es necesario cortar dichas zonas que están
fuera de los puntos de control, como los recuadros de la figura anterior; cabe aclarar, que este
resultado presenta alturas ortométricas correguidas con el modelo geoidal del IGM. El resultado
final se presenta en 4.4.9 “Resultado del mosaico MDS con ajuste en bloque y 13 GCP en Z/I
La rectificación de las imágenes se realiza en el programa Geomedia el cual necesita del mosaico
MDS del apartado 3.6.5 y de las imágenes con layers piramidales con sus tres bandas (RGB); cabe
aclarar que no se realizó el MDS en Geomedia; debido a no contar con la licencia instalada.
En Geomedia el método de interpolación utilizado es del vecino más cercano, con un tamaño
de pixel de 30cm (ver apartado 3.2.2), y con el uso del MDS anterior, tal como indica la Figura
150
115. Este proceso es de manera automática en las 24 imágenes del proyecto, cuyos resultados se
todas las imágenes; dicho ortomosaico, debe ser sometido a un proceso de edición manual, para
eliminar errores como: homogenización de colores, vista más vertical, edificios con
desplazamientos, entre otros. En la Figura 116 se indica con líneas azules los límites del MDS y de
Al finalizar la edición se inicia procesamiento del ortomosaico, el cual debe ser cortadas y
eliminadas las zonas que no tienen triple traslapo (rectángulos), pues están fuera de la zona de
El resultado del ortomosaico se presenta en el apartado 4.5.9 “Resultado del ortomosaico con
A continuación se detallan los procedimientos del control de calidad de los mosaicos de MDS
y la exactitud posicional de los ortomosaicos generados con diferentes puntos de control y métodos
de aerotriangulación.
programas fotogramétricos distintos (LPS, Orima, ZI/Imagine); por tanto, llevan una nomencla tura
individual, a fin de identificar mucho más rápido en una Geodatabase creada en Arcgis.
La primera estructura de la información con los resultados de los mosaicos de MDS, los cuales
Tabla 22.
Nomenclatura de los mosaicos MDS en la geodatabase
No Programa GCP Aerotriangulación Proceso Nomenclatura
Erdas Orima 19 13 Sin Con Ajuste Líneas M DS final
pesos pesos en de M osaico
Erdas bloque vuelo
Nomenclatura E O Y X I ii iii iv VI_X
1 √ √ √ √ E_X_i_VI_X
2 √ √ √ √ E_X_ii_VI_X
3 √ √ √ √ E_Y_i_VI_X
4 √ √ √ √ E_Y_ii_VI_X
Continúa
153
5 √ √ √ √ O_X_iii_VI_X
6 √ √ √ √ O_X_iv_VI_X
7 √ √ √ √ O_Y_iii_VI_X
8 √ √ √ √ O_Y_iv_VI_X
Nota: El mosaico MDS de Z/I Imagine se denomina Z_x_iii_vi_x
métodos, puntos de control y programas fotogramétricos, que a su vez llevan una nomencla tura
parecida a la tabla 22; sin embargo, cambian en su parte final como se observa en la tabla 23
“Base_datos_Ortomosaicos”.
Tabla 23.
Nomenclatura de los ortomosaicos en la geodatabase.
No Programa GCP Aerotriangulación Proceso Nomenclatura
Erdas Orima 19 13 Sin Con Ajuste Líneas Orto final
pesos pesos en de vuelo M osaico
Erdas bloque
Nomenclatura E O Y X i ii iii iv VII_X
1 √ √ √ √ E_X_i_VII_X
2 √ √ √ √ E_X_ii_VII_X
3 √ √ √ √ E_Y_i_VII_ X
4 √ √ √ √ E_Y_ii_VII_X
5 √ √ √ √ O_X_iii_VII_X
6 √ √ √ √ O_X_iv_VII_X
7 √ √ √ √ O_Y_iii_VII_X
8 √ √ √ √ O_Y_iv_VII_ X
Nota: El ortomosaico de Z/I Imagine se denomina Z_x_iii_vii_ x
estructuración de la información”.
154
Según la revisión bibliográfica del apartado 2.6.6, el tamaño de la muestra en los puntos de
for Digital Geospatial Data de la ASPRS para un área menor a 500km2 ; criterio que es compartido
por diferentes normas en la cantidad mínima de muestreo; así mismo, coincide con lo establecido
en el Ecuador por el IGM; así pues, al evaluar la exactitud posicional en este proyecto, se estableció
Los puntos utilizados en la exactitud posicional, son diferentes a los GCP del ajuste
fotogramétrico y no deben participar en ningún proceso del ajuste; además, están distribuidos en
distintos lugares de los GCP; y tienen las siguientes características: promedio en PDOP de 2.44
con una desviación estándar de 0.47; además, el promedio de las precisiones horizontales de 0.6cm
con un máximo de 3.3cm y un promedio de 0.7cm en las precisiones verticales con un máximo de
A continuación, se eligen los puntos de la base de datos de la tabla 26, los mismos que fueron
organizados en tres grupos para evaluar la exactitud posicional; el primer grupo, de 23 puntos, se
El segundo grupo de puntos, consta de 23 puntos para evaluar productos generados con 19GCP,
El tercer grupo de puntos consta de 23 GCP, para productos generados en Z/I Imagine ilustrados
en la Figura 120.
156
Figura 120. Puntos de exactitud posicional para productos del ajuste en Z/I Imagine
Además, los identificadores se presentan en el apartado 4.6.2 de los resultados tabla 28.
Al realizar la calidad de los mosaicos de MDS, se debe importar cada MDS, uno a la vez, de la
base de datos “Base_datos_MDS” e importar los puntos de control de la tabla 28 (con sus
Dependiendo del producto a evaluar, por ejemplo: 13GCP o 19GCP, se utilizarán los 23 puntos
de calidad para cada caso y se extrae el valor del pixel (valor de altura elipsoidal en el caso de LPS
u Orima); así también, se extrae el valor de altura ortométrica en el caso del mosaico MDS producto
del ajuste en Z/I Imagine, pues los puntos de control tenían el modelo geoidal propuesto por el
IGM.
157
Según Kirby (2007) y Mancerto et al (2015), la extracción del valor del pixel se realiza con la
herramienta “extract values to points”; a fin de determinar las diferencias entre la altura del raster
Una vez obtenidos los valores de la diferencia de altura entre el valor del raster y el del control
(altura elipsoidal en productos LPS u Orima y altura ortométrica en productos del ajuste Z/I
Imagine), se procede a calcular los valores estadísticos como: mínimo, máximo, media y desviación
estándar, RMSEz vistos en las “Definiciones de estadísticos” del apartado 2.7.1; así pues, se realiza
el proceso de evaluación para cada mosaico de MDS detallado a continuación, cuyos resultados se
La validación de los mosaicos de MDS, representa una dificultad extra, pues en este proyecto se
utilizó la altura elipsoidal y los estándares de confiabilidad trabajan con MDE; es decir, eliminan
puntos, como LIDAR para mejorar el producto (Soria, 2014); sin embargo, los puntos de exactitud
están tomados a nivel de piso, eliminando así ciertas equivocaciones que pudieran ocurrir; además, se
puede calcular el límite de aceptación con el criterio expuesto por el IGM, el cual indica el intervalo de
para escala 1:5000, considerando la norma “Positional Accuracy Standards for Digital Geospatial
Data” de la ASPRS y la propuesta por el IGM para fiscalizar proyectos obtenidos por métodos
fotogramétricos; en los cuales, se utilizan las definiciones de RMSEr, así como la utilización de las
ecuaciones 14 y 15 de los casos vistos en el apartado 2.6.7, con cierto nivel de confianza. Dicho
de puntos, mismos ID que fueron utilizados en el MDS) por cada ortomosaico; posteriormente, con
una escala adecuada (1:150, para este estudio) y fija se procede a pinchar cada detalle de acuerdo
a la monografía, hasta obtener todos los puntos pinchados (ver Figura 122).
159
Una vez marcados los puntos en la imagen, se procede a calcular las coordenadas, Este (x) con
la herramienta “Calculate Geometry”, sistema de referencia WGS84 UTM zona 17 Sur y luego, la
Posteriormente, se realiza la diferencia entre las coordenadas Este y Norte calculadas del
ortomosaico, con las obtenidas en los puntos de exactitud posicional (ver tabla 26), cuyos ID se
indican en la tabla 28; generando así, un Δx y un Δy, a los cuales se les aplica la raíz cuadrada de
160
sus componentes al cuadrado para representar el error planimétrico (Pintos, 2016); tal y como se
Dichos valores se pueden representar en forma de círculos a lo largo del proyecto y para pada
resultado, indicados en el apartado 4.6.4 “Resultados del error planimétrico de los ortomosaicos”,
mínimo, error medio (ecuación 18), desviación estándar (ecuación 19), RMSEx, RMSEy (ecuación
13 para X e Y) y RMSEr, indicados en la tabla 32 del apartado 4.6.5. Este es un proceso secuencial,
RMSEmin/RMSEmax, para su aplicación en el tipo de ecuación, según los casos del apartado 2.6.7;
donde, los valores menores a 0.6 indican la aplicabilidad de la ecuación 14, (ver tabla 33).
la normativa IGM (ver apartado 2.6.7), considera dos posibles casos; donde, RMSEx ≠ RMSEy, aplicar
la ecuación 15 bajo las condiciones de RMSEmin/RMSEmax entre 0.6 y 1 (FGDC, 1998), con un 90%
de confiabilidad; si dicho valor no se encuentra en dicho rango se aplica el RMSEr con la ecuación 14,
Los resultados de dichos cálculos se presentan en la tabla 34 del apartado 4.6.5; así como, la
escala máxima que podrían ser utilizados dichos ortomosaicos; además, el límite de la exactitud
posicional se calcula con la normativa ISO (0.2mm), al realizar el producto del módulo de la escala
En el caso de la normativa ASPRS, el límite de la exactitud posicional viene dado por la tabla 11
(ver apartado 2.6.7) cuyo valor máximo al 95% de confiabilidad es de 1.101m, para un GSD de 30cm
aproximadamente y un RMSEr de 63.6cm; eso quiere decir, que el cálculo se realiza de la siguiente
𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅
𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸𝐸 = 2.448 ∗ = 1.731 ∗ 𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅𝑅
√2
Los resultados de dichos cálculos se presentan en tablas del apartado 4.6.5 “Resultados estadísticos
de los ortomosaicos”.
162
CAPITULO IV
4. RESULTADOS
En la tabla 24 se presenta los resultados de la depuración de los datos del archivo APLANIX,
Tabla 24.
Resultados del Aplanix a archivo IMU/GPS
ID ID Este [m] Norte [m] Altura Omega [°] Phi [°] Kappa [°]
Imagen Elipsoidal
[m]
1 18_537 779600.63 9967843.6 8169.182 0.0416 -0.04853 90.37356
2 18_538 779601.626 9966757.46 8173.792 -0.01069 -0.03016 90.36288
3 18_539 779606.87 9965671.8 8176.02 0.00764 0.02138 90.77048
4 18_540 779608.866 9964587.24 8182.397 -0.00806 0.01311 90.72516
5 18_541 779611.025 9963501.22 8188.567 -0.01906 -0.02381 90.49753
Continúa
163
fotogramétrico tomados por Barragán (2015) con equipo GNSS de doble frecuencia (L1/L2).
Tabla 25.
Elenco de coordenadas de las monografías de Barragán (2015)
ID NORTE [m] ESTE [m] Precisión Precisión Altura elipsoidal
Horizontal [m] Vertical [m] [m]
GCPC 9967656.246 781210.625 0.012 0.014 2565.285
GCPG 9963447.427 787229.463 0.005 0.009 2560.211
GCPI 9967190.456 787605.817 0.007 0.010 2554.092
GCPX 9968275.126 777028.821 0.013 0.028 3007.168
GCPY 9962161.490 791186.294 0.014 0.025 2749.505
GCPE 9965033.915 784443.902 0.001 0.002 2520.095
Continúa
164
“Bloque_ruminahui_gruber.blk” con la orientación relativa (72 puntos de Von Gruber) (ver Figura
125).
A continuación se presentan las distribuciones de GCP, los resultados del procesamiento y sus
monografías.
165
En las campañas de campo se logró obtener 19 puntos de control modificados, de los planteados
A continuación se presenta un resumen de los 45 puntos de control (ver tabla 26) en donde la
columna “uso” representa la participación del punto: 1300 en “Distribución 13”, 1900
“Distribución 19”, 1319 puntos compartidos en ambas distribuciones y “9999” puntos solamente
en control de calidad; además, los puntos obtenidos de Barragán (2015), tienen una codificac ió n
Tabla 26.
Detalles de los puntos a utilizar en el proyecto
No Mon. Norte [m] Este [m] Altura Elip P. P. Uso PDOP RMS
[m] Horizontal Vertical [m]
[m] [m]
1 A001 9967656.246 781210.625 2565.285 0.012 0.014 1319 3.901 0.001
2 A004 9968275.126 777028.821 3007.168 0.013 0.028 1319 2.738 0.001
3 A005 9962161.490 791186.294 2749.505 0.014 0.025 1900 2.262 0.002
4 A006 9965033.915 784443.902 2520.095 0.001 0.002 1300 1.831 0.001
5 A007 9963512.875 780761.586 2538.293 0.005 0.006 9999 1.792 0.001
6 A008 9963076.165 784612.941 2551.316 0.003 0.005 1300 2.645 0.001
7 A009 9966469.118 784048.411 2500.598 0.002 0.003 9999 1.795 0.001
8 A010 9961837.692 777372.520 2839.422 0.005 0.011 1319 1.85 0.001
9 A011 9968378.800 791497.809 2566.699 0.007 0.009 1900 2.55 0.001
10 B012 9965493.415 784449.746 2536.315 0.003 0.004 1900 2.31 0.035
11 B013 9964861.857 784593.351 2533.752 0.002 0.004 9999 2.943 0.019
12 C014 9967772.804 777584.835 2855.623 0.005 0.006 1900 2.351 0.017
13 C015 9962820.938 777846.513 2739.132 0.008 0.012 1900 3.65 0.05
14 C016 9965824.004 779092.843 2635.156 0.003 0.007 9999 3.055 0.032
15 C017 9966273.223 782002.918 2495.664 0.003 0.007 1900 2.234 0.028
16 D018 9967035.860 789791.892 2597.935 0.004 0.005 9999 1.661 0.008
17 D019 9968418.687 791501.879 2566.403 0.003 0.006 1319 4.006 0.022
18 D020 9962277.748 791243.853 2745.658 0.005 0.007 1319 2.191 0.025
Continúa
167
La distribución espacial de todos los puntos se presenta en la Figura 128; en el cual el color
verde representa los puntos de exactitud posicional y los triángulos rojos los puntos de control
(GCP).
168
El detalle de las características de cada punto se representa en una monografía, cuyo ejemplo se
las monografías de los puntos mencionados en la tabla 17, exceptuando el punto X100, debido a
un PDOP excesivo.
dos ajustes con la distribución de 19 puntos de control, considerando la matriz de pesos y sin
I044 e I048; además, tiene un RMSE de 0.4479m en el eje X de las imágenes, un RMSE de 0.4079m
Así mismo, en el reporte del ajuste se incluyen varios datos, entre ellos el total de puntos
considerados en el ajuste (313 tie points), de los cuales 148 de ellos han sido pinchados en más de
4 imágenes. En la Figura 131 se presenta un ejemplo del reporte de LPS, cuyo contenido completo
del reporte; así como, todos los reportes de aerotriangulación de los diferentes programas se
presentan en el Anexo 4.
A010, I044 e I048; además, tiene un RMSE de 0.3787m en el eje X de las imágenes, un RMSE de
0.3334m en el eje Y de las imágenes y un RMSE total de 0.2732m (ver Figura 132).
171
También, se genera un reporte en el que incluye el total de los puntos considerados en el ajuste
(307 tie points); de los cuales 143 de ellos se han pinchados en más de 4 imágenes; el resto del
A004, A005 A010, A001, C014, D019, e I044; además, tiene un RMSE de 0.3146m en el eje X
de las imágenes, un RMSE de 0.2961m en el eje Y de las imágenes y un RMSE total de 0.2627m
considerados en el ajuste (310 tie points); de los cuales 143 de ellos, se encuentran pinchados en
A004, A005 A010, A001, C014, D019, I044 e I048; además, tiene un RMSE de 0.2853m en el eje
X de las imágenes, un RMSE de 0.2546m en el eje Y de las imágenes y un RMSE total de 0.2595m
De la misma forma, se puede generar un reporte de aerotriangulación, que incluye el total de los
puntos utilizado en el ajuste (310 tie points), de los cuales 144 de ellos han sido pinchados en más
dos ajustes con la distribución de 13 GCP, considerando un ajuste en bloque y un ajuste en línea
173
de vuelo; así mismo, dos ajustes con la distribución de 19 GCP, considerando un ajuste en bloque
se consideró los puntos de control: A010, I044 e I048, debido a residuales muy altos (RMS de
0.37m en cada componente). El RMSE de las componentes del ajuste resultante son: RMSE X
Al correr el programa CAP-A del ajuste de Orima genera un reporte, con la cantidad de puntos
con errores que sobrepasan los límites y no han sido eliminados (Im. Blunders = 2); además,
(1.6µm) con el número de iteración de convergencia (ver Figura 135), entre otros valores.
modelo, codificación de los puntos de amarre, coordenadas del centro de las imágenes, ángulos de
rotación, entre otros datos (ver Figura 136), detallados en el Anexo 6 (Tema: Apertura de Orima);
A su vez, la cantidad total de puntos utilizados en el ajuste fue de 781, de los cuales 410 de ellos
programa CAP-A de ajuste de Orima genera un reporte, con la cantidad de puntos con errores que
175
sobrepasan los límites y no han sido eliminados (Im. Blunders = 5); además, presenta la cantidad
de otros errores no corregidos (Add. Blunders = 8) y el valor de Sigma 0 (1.6 um) con el número
de iteración de convergencia (ver Figura 137), entre otros valores. Los valores de RMSE por
Además, el número total de puntos utilizados en el ajuste es de 890 tie points; donde, 391 de
ellos han sido pinchados en más de 4 imágenes. El contenido completo del reporte de
de puntos manual, al ejecutar el programa CAP-A de ajuste que genera una ventana con la cantidad
de puntos con errores que sobrepasan los límites y no han sido eliminados (Im. Blunders = 20);
Sigma0 (1.7 µm) con el número de iteración de convergencia (ver Figura 138), entre otros valores.
El RMS por componentes del ajuste son: RMSE X (0.053m), RMSE Y (0.071m) y RMSE Z
(0.018m).
Así mismo, en el reporte se presenta la cantidad de puntos utilizados en el ajuste (779 tie points),
de los cuales 417 han sido pinchados en más de 4 imágenes. El contenido completo del reporte se
presenta en el Anexo 4.
el programa CAP-A de ajuste se genera un reporte, con la cantidad de puntos con errores que
sobrepasan los límites y no han sido eliminados (Im. Blunders = 25); además, presenta la cantidad
de otros errores no corregidos (Add. Blunders = 6) y el valor de Sigma 0 (1.6µm) con el número
177
de iteración de convergencia (ver Figura139), entre otros valores. El valor de RMSE por
componentes son los siguientes: RMSE X (0.056m), RMSE Y (0.07m) y RMSE Z (0.023m).
de los cuales 400 de ellos han sido pinchados en más de 4 imágenes. El contenido completo del
Imagine con los 13 puntos de control; sin embargo, el programa eliminó automáticamente los
puntos: A010, D020, G030 y H040; además, se presenta un RMS por componentes cuyos valores
son: RMS X (0.082m), RMS Y (0.228m) y RMS Z (0.063m) (ver Figura 140) y un sigma de 2um.
178
proyecto no sobrepasa dichos valores (1.25m); además, los valores máximos de los residuales son:
0.082m, 0.228m y 0.063m en cada componente; también, en el reporte del ajuste en Z/I Imagine
se indica que se utilizaron 284 puntos; donde, 183 fueron pinchados en más de cuatro imágenes.
4.4.1. Resultado del mosaico MDS sin matriz de pesos y 13 GCP en LPS
Se generó un mosaico con los 20 MDS del ajuste con 13 GCP y sin considerar la matriz de pesos
además, se representa de color café las elevaciones grandes y de color azul las zonas más bajas (ver
Figura 141); además se presenta una zona representativa de color verde (vegetación) de ejemplo.
4.4.2. Resultado del mosaico MDS con matriz de pesos y 13 GCP en LPS
matriz de pesos (modelo Ebners) en LPS-Core, cuyo valor máximo de altura elipsoidal es de
180
3030.32m y mínimo de 2440.24m; además, se representa de color café las elevaciones grandes y
amarillo las zonas más bajas (ver Figura 142); en el recuadro se presenta una zona de ejemplo
(quebrada).
4.4.3. Resultado del mosaico MDS sin matriz de pesos y 19 GCP en LPS
Se generó un mosaico con los 20 MDS del ajuste con 19 GCP y sin considerar la matriz de pesos
además, se representa de color café las elevaciones grandes y de color azul las zonas más bajas (ver
Figura 143); en el recuadro se presenta una zona de ejemplo (estructura - edificio San Luis de
Sangolquí).
181
4.4.4. Resultado del mosaico MDS con matriz de pesos y 19 GCP en LPS
Se generó un mosaico con los 20 MDS del ajuste con 19 GCP y considerando la matriz de pesos
(modelo Ebners) en LPS-Core, cuyo valor máximo de altura elipsoidal es de 3031.04m y mínimo
de 2438.96m; además, se representa de color verde las elevaciones grandes y de color azul las
zonas más bajas (ver Figura 144); además, se presenta un ejemplo de edificios (color verde oscuro).
182
4.4.5. Resultado del mosaico MDS con ajuste en bloque y 13 GCP en Orima
Se generó un mosaico con los 20 MDS del ajuste en bloque de Orima y con 13 GCP, cuyo valor
color café las elevaciones grandes y de color azul las zonas más bajas (ver Figura 145); además, se
4.4.6. Resultado del mosaico MDS con ajuste en líneas de vuelo y 13 GCP en Orima
Se generó un mosaico con los 20 MDS del ajuste en líneas de vuelo en Orima y con 13 GCP,
representa de color café las elevaciones grandes y de color azul las zonas más bajas (ver Figura
146); así también, las zonas de color verde claro representan quebradas y la zona cafe una zona de
alta elevación.
184
Figura 146. Mosaico MDS con ajuste en líneas de vuelo y 13 GCP en Orima.
4.4.7. Resultado del mosaico MDS con ajuste en bloque y 19 GCP en Orima
Se generó un mosaico con los 20 MDS del ajuste en bloque con Orima y con 19 GCP, cuyo
color café las elevaciones grandes y de color azul las zonas uniformes (ver Figura 147.); de la
misma forma, el color rojo oscuro representa alta vegetación y el anaranjado zonas de cultivo.
185
4.4.8. Resultado del mosaico MDS con ajuste en líneas de vuelo y 19GCP en Orima
Se generó un mosaico con los 20 MDS del ajuste en líneas de vuelo en Orima y con 19 GCP,
representa de color café las elevaciones grandes y de color azul las zonas más bajas (ver Figura
Figura 148. Mosaico MDS con ajuste en líneas de vuelo y 19 GCP en Orima.
4.4.9. Resultado del mosaico MDS con ajuste en bloque y 13 GCP en Z/I Imagine
Se generó un mosaico del ajuste en bloque y con 13 GCP del programa ZI-Imagine, cuyo valor
color café las elevaciones grandes y de color amarillo las zonas más bajas (ver Figura 149) (zonas
de ladera).
187
Figura 149. Mosaico MDS con 13GCP y ajuste bloque en Z/I Imagine.
del ajuste de Erdas-LPS, 96 ortofotos generadas del ajuste de ORIMA y 24 ortofotos generadas
Figura 152; donde, en el recuadro verde se observa una ligera curvatura de la edificación color
blanco oeste; sin embargo las edificaciones del recuadro rojo y del recuadro amarillo no presentan
Se generaron 24 ortofotos con el ajuste de 19 GCP y considerando una matriz de pesos (modelo
recuadro verde de la Figura 154) presenta curvaturas en sus extremos; entre otros ejemplos.
4.5.6. Resultado del ortomosaico con ajuste de líneas de vuelo y 13 GCP en Orima
Se generaron 24 ortofotos con 13 GCP y un ajuste en líneas de vuelo con el entorno de Orima,
estructuras.
Del recuadro verde anterior, se visualiza la curvatura de un edificio en sus extremos (bloque de
así también, en el recuadro de color rojo se presenta un edificio (centro comercial San Luis –
4.5.8. Resultado del ortomosaico con ajuste en líneas de vuelo y 19 GCP en Orima
Se generaron 24 ortofotos con 19 GCP y un ajuste en líneas de vuelo con el entorno de Orima,
la Figura 157.
193
4.5.9. Resultado del ortomosaico con ajuste en bloque y 13 GCP en Z/I Imagine
la Figura 158.
Figura 158. Ortomosaico con 13 GCP y ajuste en bloque con Z/I Imagine.
194
test de calidad aplicando la normativa IGM y la propuesta por norma “Positional Accuracy
El modelo lógico de las bases de datos es de tipo raster tanto en los resultados de los mosaicos
de MDS y los resultados de los ortomosaicos, cuyos nombres se detallan en la tabla 27.
Tabla 27.
Modelo lógico de las geodatabases
GDB Raster Dataset GDB Raster Dataset
Base_datos_MDS E_X_i_VI_X Base_datos_Ortomosaico E_X_i_VII_X
E_X_ii_VI_X s E_X_ii_VII_X
E_Y_i_VI_X E_Y_i_VII_X
E_Y_ii_VI_ X E_Y_ii_VII_X
O_X_iii_VI_X O_X_iii_VII_X
O_X_iv_VI_X O_X_iv_VII_X
O_Y_iii_VI_ X O_Y_iii_VII_X
O_Y_iv_VI_X O_Y_iv_VII_X
Z_x_iii_vi_x Z_x_iii_vii_x
195
23 puntos de control en la calidad de los MDS al extraer su valor de altura; además, se considera
que los puntos no deben intervenir en ningún proceso del ajuste fotogramétrico; por tanto, se
Tabla 28.
Puntos utilizados para el control de calidad
Puntos de calidad compartidos por todos los mosaicos No puntos Total de puntos
A007 D021 G032 H039 4 15
B013 E022 G033 H041 4
C016 G028 H037 I043 4
D018 G031 H038 3
Puntos de calidad solo para productos generados con 13 GCP 15 23
A005 A011 C015 G029 4
A009 C014 E023 G034 4
Puntos de calidad solo para productos generados con 19 GCP 15 23
A006 A009 G029 H036 4
A008 E023 G034 H040 4
Puntos de calidad solo para productos generados en 15 23
ZI/Imagine
A005 C015 G035 I047 4
B012 C017 H042 A011 3
el programa fotogramétrico (LSP, Orima y Z/I Imagine) y el tipo de ajuste, cuyos parámetros
estadísticos utilizados fueron: máximo, mínimo, error medio (ecuación 18), media aritmética
196
absoluta, desviación estándar (ecuación 19); además, cabe mencionar que el análisis del mosaico
de MDS del ajuste en Z/I Imaging, se eliminó un punto, debido a encontrarse fuera de la zona de
estudio (punto A011); por tanto, cada uno de los valores calculados se muestran en la tabla 29.
Tabla 29.
Parámetros estadísticos de los mosaicos de MDS
MDS Evaluado con Elementos Máximo Mínimo Error Media Desviación desviación
altura elipsoidal (n) [m] [m] Medio arimetica estándar estándar2
(x� )[m] absoluta (S) [m] (S) [m]
(x� )1 [m]
E_X_i_VI_X 23 1.563 -2.362 -0.245 0.871 1.106 0.701
E_X_ii_VI_X 23 1.403 -2.091 -0.145 0.809 1.001 0.582
E_Y_i_VI_X 23 1.3 -3.638 -0.515 0.900 1.129 0.841
E_Y_ii_VI_X 23 0.99 -2.581 -0.500 0.790 0.929 0.680
O_X_iv_VI_X 23 1.416 -2.294 0.072 0.751 1.017 0.671
O_X_iii_VI_X 23 1.109 -2.035 0.051 0.637 0.841 0.535
O_Y_iv_VI_X 23 1.048 -2.163 0.155 0.632 0.760 0.430
O_Y_iii_VI_X 23 1.406 -2.281 0.225 0.678 0.835 0.518
z_x_iii_vi_x3 22 3.895 -5.192 -0.202 1.682 2.249 1.463
Promedio 1.57 -2.737 -0.123 0.861 0.270 0.322
1. Representa el valor de la media aritmética del valor absoluto de las diferencias.
2. Representa el valor de la desviación estándar con el valor absoluto de las diferencias entre el valor del MDS
producido y el Control
3. Representa el MDS con correcciones de altura elipsoidal a ortométrica con el modelo IGM
A continuación se calcula el nivel de confianza del MDS, el cual resulta de despejar 𝑡𝑡𝛼𝛼�2 de la
ecuación 16 (ver apartado 2.7.1); además, como los datos son menores a 30 se utiliza la tabla t –
𝜏𝜏 . √ 𝑛𝑛 [ 23]
𝑡𝑡 𝛼𝛼�2 =
𝜎𝜎
Donde:
- 𝜏𝜏, representa el valor del error permito o tolerado en el proceso y resulta del interva lo
Al obtener el valor 𝑡𝑡 𝛼𝛼�2 de la ecuación anterior, se busca el valor interpolado de t con n-1, grados
de libertad (n=22 y n=21) en la tabla t-Student (t0 ); donde, t% indica el valor en porcentaje de t0 ,
Tabla 30.
Nivel de confianza para los MDS
MDS Evaluado con t t0 1 t% %
altura elipsoidal
E_X_i_VI_X 5.42 2.82 0.995 99.5
E_X_ii_VI_X 5.99 2.82 0.995 99.5
E_Y_i_VI_X 5.31 2.82 0.995 99.5
E_Y_ii_VI_X 6.45 2.82 0.995 99.5
O_X_iv_VI_X 5.89 2.82 0.995 99.5
O_X_iii_VI_X 7.13 2.82 0.995 99.5
O_Y_iv_VI_X 7.89 2.82 0.995 99.5
O_Y_iii_VI_X 7.18 2.82 0.995 99.5
z_x_iii_vi_x2 2.61 2.52 0.99 99.0
[1] Representa el valor del área bajo la curva de la tabla t-
Student
Continúa
198
Así también, se puede encontrar el intervalo de confianza para cada uno de los MDS, al utilizar
la ecuación 17 (ver apartado 2.7.1) de la siguiente manera, cuyos resultados se presenta en la tabla
31:
𝑡𝑡 𝛼𝛼�2 . 𝑆𝑆 [17]
𝑥𝑥̅ ±
√ 𝑛𝑛
Donde:
- 𝑡𝑡 𝛼𝛼�2, representa el valor del área bajo la curva de la distribución t-student, a cierto nivel
Tabla 31.
Intervalo de confianza de los MDS
MDS Evaluado con Error Intervalo de
altura elipsoidal �)
Medio (𝒙𝒙 confianza al
[m] 95% [m]
E_X_i_VI_X -0.245 ± 0.40
E_X_ii_VI_X -0.145 ± 0.36
E_Y_i_VI_X -0.515 ± 0.40
E_Y_ii_VI_X -0.500 ± 0.33
O_X_iv_VI_X 0.072 ± 0.36
O_X_iii_VI_X 0.051 ± 0.30
Continúa
199
Al analizar el valor del RMSEz con la ecuación 13 (ver apartado 2.6.7) y multiplicar por el valor
de “Z” (tabla z=1.645), con un nivel de confianza del 90%, se tienen los valores de la tabla 32:
Tabla 32.
Precisión de los MDS al 90% de confianza
MDS Evaluado con RMSEz Precisión al 90%
altura elipsoidal [m] de confianza [m]
E_X_i_VI_X 1.109 1.82
E_X_ii_VI_X 0.989 1.62
E_Y_i_VI_X 1.219 2.00
E_Y_ii_VI_X 1.037 1.70
O_X_iv_VI_X 0.997 1.64
O_X_iii_VI_X 0.824 1.35
O_Y_iv_VI_X 0.759 1.24
O_Y_iii_VI_X 0.847 1.39
z_x_iii_vi_x 2.207 3.62
Promedio 1.109 1.821
[1] Resultados de multiplicar el promedio obtenido con
el valor de 1.645 para el 90% de confianza.
El límite de un mosaico de MDS a escala 1:5000, es de 1.25m, pues implica el ¼ del interva lo
de curva; por tanto, de los valores de la tabla 32, se tiene que a un nivel de confianza del 90% solo
un producto cumple con la precisión indicada (O_Y_iii_VI_X), el cual representa al resultado con
19GCP, ajuste en líneas de vuelo en Orima; al mismo tiempo, el valor más alto es el realizado en
Z/I Imagine; debiéndose a utilizar el modelo geoidal del IGM, pues para el resto de los resultados
no se presentan dichos extremos del valor; no obstante, existe otro producto que se acerca al límite
200
del ¼ del intervalo de curva, el cual representa al mosaico de MDS con 13GCP, ajuste en bloque
3
2.5
2
1.5
1
0.5
0
O_Y_iv_VI O_X_iii_VI O_Y_iii_VI E_X_ii_VI O_X_iv_VI E_Y_ii_VI E_X_i_VI_ E_Y_i_VI_ z_x_iii_vi_
_X _X _X _X _X _X X X x
RMSEr [m] 1.24 1.35 1.39 1.62 1.64 1.7 1.82 2 3.62
uno en círculos; donde, el radio de las mismas está dado por la ecuación 22.
cual representa el resultado del apartado 4.5.1 “Resultado del ortomosaico sin matriz de pesos y
13GCP en LPS”, cuyo error planimétrico de los puntos de exactitud posicional van de un valor
mínimo de 0.07m y un máximo de 1.06m; además, la mayor parte de los errores se encuentra entre
las coordenadas Este [77600m; 77800m], pertenecientes a una zona de alta pendiente, el resto de
9969000
0.74
9968000 0.49
1.06 0.27
0.10
0.63 0.24
9967000
0.49 0.31 0.62
NORTE [m]
9966000 0.67
0.39 0.46
9965000 0.07
0.85 0.25
0.25
9964000 0.07
0.46 0.30
9963000 0.36
0.51
9962000 0.77
9961000
776000 778000 780000 782000 784000 786000 788000 790000 792000
ESTE [m]
cual representa el resultado del apartado 4.5.2 “Resultado del ortomosaico con matriz de pesos y
13 GCP en LPS”, cuyo error planimétrico de los puntos de exactitud posicional, van de un valor
mínimo es de 0.1m y un máximo de 1.9m; además, la mayor parte de los errores se encuentra entre
las coordenadas Este [77600m; 78000m], pertenecientes a una zona de alta pendiente, el resto de
representa el resultado del apartado 4.5.3 “Resultado del ortomosaico sin matriz de pesos y 19 GCP
círculos), van de un valor mínimo es de 0.04m y un máximo de 0.88m; además, los valores del
0.88 0.67
0.750.61 0.07
9965000 0.80 0.69
0.23
0.34 0.06
0.59 0.22
0.36
9963000 0.04
0.63
9961000
776000 778000 780000 782000 784000 786000 788000 790000 792000
ESTE [m]
representa el resultado del apartado 4.5.4 “Resultado del ortomosaico con matriz de pesos y 19
GCP en LPS”; así mismo, los errores planimétricos de los puntos de exactitud posicional van de
un valor mínimo es de 0.04m y un máximo de 0.91m; además, los errores se distribuyen en zonas
de alta pendiente (0.91m) y zonas extremas del proyecto (0.59m), tal y como se indican los círculos
en la Figura 163.
0.72 0.19
9965000 0.38
0.04
0.91 0.07
0.36
9964000 0.26
0.18 0.49
9963000 0.10
0.74
9962000
9961000
776000 778000 780000 782000 784000 786000 788000 790000 792000
ESTE [m]
que representa el resultado del apartado 4.5.5 “Resultado del ortomosaico con ajuste en bloque y
204
13 GCP en Orima”, cuyo error planimétrico máximo es de 1.6m y un mínimo de 0.1m; además,
los mayores errores se encuentran en las esquinas del proyecto y en las zonas de alta pendiente, en
0.39 0.50
9966000 0.44
0.64 0.30
9965000 0.53
0.84 0.07
0.13
9964000 0.20
0.14 0.62
9963000 0.89 0.57
9962000 0.72
9961000
776000 778000 780000 782000 784000 786000 788000 790000 792000
ESTE [m]
que representa el resultado del apartado 4.5.6 “Resultado del ortomosaico con ajuste de líneas de
0.02m, el resto de valores son menores a 0.8m representados en círculos de la Figura 165.
205
que representa el resultado del apartado 4.5.7 “Resultado del ortomosaico con ajuste en bloque y
19 GCP en Orima”, cuyo error planimétrico máximo es de 0.9m, el resto de errores son menores a
0.7m representados en círculos; vale destacar, los mayores desplazamientos se ubican entre la línea
de vuelo 20 y 21 del proyecto, cuyas coordenadas Este son [78600m; 792000m] (ver Figura 166).
0.90 0.65
0.660.22 0.25
9965000 0.63 0.38
0.17
0.21 0.18
0.66 0.13
0.13
9963000 0.23
0.68
9961000
776000 778000 780000 782000 784000 786000 788000 790000 792000
ESTE [m]
que representa el resultado del apartado 4.5.8 “Resultado del ortomosaico con ajuste en líneas de
vuelo y 19 GCP en Orima”, cuyo error planimétrico máximo es de 0.8m y se ubica en el extremo
0.32 0.14
9965000 0.68
0.21
0.65 0.11
0.12
9964000 0.38
0.28 0.55
9963000 0.10
0.84
9962000
9961000
776000 778000 780000 782000 784000 786000 788000 790000 792000
ESTE [m]
representa el resultado del apartado 4.5.9; donde, el valor máximo es de 1.1m ubicado en el extremo
Este del proyecto; además, en la zona de alta pendiente los errores planimétricos son menores a
0.64 0.44
0.13
9965000 0.29
0.28 0.18 0.75
9964000 0.56
0.24
9963000 0.87
0.91
9962000 0.36
9961000
776000 778000 780000 782000 784000 786000 788000 790000 792000
ESTE [m]
A continuación se detallan los siguientes estadísticos: máximo, mínimo, error medio, desviación
estándar, RMSE (RMSEx y RMSEy) y RMSEr, indicados en la tabla 33 de los nueve ortomosaicos.
Tabla 33.
Estadísticos generales de los ortomosaicos.
No Ortomosaico Coordenada Estadísticos generales
No Máx. Mín. Error Desviación RMSE RMSEr
elemen. [m] [m] Medio Estándar [m] [m]
[m] [m]
1 E_X_i_VI_X ΔEste [E] 23 0.377 -0.615 -0.155 0.272 0.308 0.515
ΔNorte [N] 23 0.714 -1.032 -0.125 0.403 0.413
2 E_X_ii_VI_X ΔEste [E] 23 0.377 -0.642 -0.026 0.264 0.26 0.721
ΔNorte [N] 23 0.754 -1.865 -0.318 0.607 0.673
3 E_Y_i_VI_X ΔEste [E] 23 0.655 -0.761 -0.138 0.406 0.421 0.498
ΔNorte [N] 23 0.476 -0.516 -0.105 0.250 0.266
4 E_Y_ii_VI_ X ΔEste [E] 23 0.496 -0.734 -0.132 0.346 0.363 0.477
Continúa
208
aplicación en el tipo de ecuación, según los casos del apartado 2.6.7; donde, el valor 0.6 de la tabla
Tabla 34.
Indicador del RMSEmin / RMSEmax
Productos RMSEx RMSEy RMSEmin /
[m] [m] RMSEmax
Continúa
209
El límite de aceptación del producto (ortomosaico) a escala 1:5000 con la normativa ISO es de
1m; donde, el único producto fuera de dicho valor es el realizado con 13GCP y considerando una
matriz de pesos (Ebner’s) en LPS al 90% de confianza; en contra posición a este, el mejor resultado
se obtiene con 19 GCP y ajuste en líneas de vuelo del programa Orima al 90% de confianza;
además, en la columna “Escala máxima de ampliación”, se indica a que escala podría llegar si se
Tabla 35.
Nivel de confianza de los ortomosaicos con la normativa IGM.
No Productos RMSEx RMSEy RMSEr IGM al Escala
[m] [m] [m] 90% máxima de
[m] ampliación
1 E_X_i_VI_X 0.308 0.413 0.515 0.774 1:4000
2 E_x_ii_vi_x 0.26 0.673 0.721 1.095 1:5600
3 E_Y_i_VI_X 0.421 0.266 0.498 0.737 1:3800
4 E_y_ii_VI_X 0.363 0.309 0.477 0.721 1:3800
5 O_X_iii_VI_X 0.273 0.578 0.639 0.970 1:5000
6 O_X_iv_VI_X 0.287 0.386 0.481 0.722 1:3800
7 O_Y_iii_VI_ X 0.309 0.341 0.460 0.697 1:3600
8 O_Y_iv_VI_X 0.283 0.314 0.423 0.641 1:3400
9 Z_X_iii_VI_X 0.414 0.371 0.556 0.842 1:4400
Promedio 0.324 0.406 0.324 0.783 1:4000
En la norma “Positional Accuracy Standards for Digital Geospatial Data” de la ASPRS, existen
dos resultados fuera del rango esperado (1.101m), los resultados 2 y 5 de la tabla 35, de los producto
210
con 13 GCP en LPS y Orima que consideran un ajuste con matriz de pesos y ajuste en bloque
respectivamente; así también, el resultado con mejor exactitud posicional representa el ajuste
realizado con 19GCP en ajuste de líneas de vuelo en Orima, tal y como se muestra en la tabla 36.
Tabla 36.
Exactitud de los ortomosaicos con la normativa ASPRS.
No Productos RMSEr [m] ASPRS al 95%
[m]
1 E_X_i_VI_X 0.515 0.892
2 e_x_ii_vi_x 0.721 1.249
3 E_Y_i_VI_X 0.498 0.862
4 E_y_ii_VI_X 0.477 0.825
5 O_X_iii_VI_X 0.639 1.107
6 O_X_iv_VI_X 0.481 0.833
7 O_Y_iii_VI_ X 0.460 0.797
8 O_Y_iv_VI_X 0.423 0.732
9 Z_X_iii_VI_X 0.556 0.962
Promedio 0.519 0.899
Vale destacar que según el orden del RMSEr, en forma ascendente se tiene tal y como se ilustra
en la en la Figura 169.
211
RMSEr [m]
0.800
0.700
0.600
0.500
RMSEZ [M]
0.400
0.300
0.200
0.100
0.000
O_Y_iv_VI O_Y_iii_VI E_y_ii_VI O_X_iv_VI E_Y_i_VI_ E_X_i_VI_ Z_X_iii_VI O_X_iii_VI e_x_ii_vi_
_X _X _X _X X X _X _X x
RMSEr [m] 0.423 0.460 0.477 0.481 0.498 0.515 0.556 0.639 0.721
CAPITULO IV
5. CONCLUSIONES Y RECOMENDACIONES
5.1. Conclusiones
En efecto, se cumplieron todos los objetivos y metas propuestas en este proyecto; es decir, se
en diferentes programas fotogramétricos (ver apartados 3.4, 3.5, 3.6 y anexos 3 al 8), con su
y distribución de GCP, en Erdas LPS Core, Orima y Z/I Imagine; se puede decir, que cada programa
(ajuste con matriz de pesos, ajuste en bloque, ajuste en líneas de vuelo), pero utilizando el concepto
los mejores resultados de exactitud posicional, están vinculados a las teorías de Ackermann (1990)
213
e IGN de Perú (2011), con ligeras variaciones en la ubicación de los puntos, con una cantidad de
19 GCP; pues, consideran puntos dobles en las esquinas del proyecto, puntos perpendiculares a las
líneas de vuelo y en zonas de triple traslapo distribuidos homogéneamente; además, se puede decir
que la distribución de 13 GCP alcanza buenas precisiones (ver tabla 35 y tabla 36) siguiendo la
metodología planteada.
Al mismo tiempo al comparar la exactitud posicional de las ortofotos con el umbral de la norma
ISO (1m), solo el resultado con 13GCP en LPS Core (ver apartado 4.5.2) se encuentra fuera del
rango permitido (1.10m al 90% de confianza); a su vez, al utilizar el límite de la norma propuesta
por la ASPRS (1.1m), las ortofotos generadas con 13GCP en LPS Core y Orima, sobrepasan el
límite con 1.25m y 1.11m respectivamente, al 95% de confianza; sin embargo, si se considera la
Además, cabe mencionar que la mejor exactitud posicional de los productos obtenidos, es de
64cm al 90% de confianza con la normativa IGM y de 73cm al 95% de confianza con la normativa
ASPRS; por tanto, se puede utilizar dicho ortomosaico en la generación de cartografía 2D para
escala 1:5000, hasta escala 1:3400 (ver tabla 35), pudiendo utilizarlo con fines catastrales, anális is
En cuanto a los modelos digitales del terreno, solo un producto, que considera 19GCP en ajuste
de Orima, es admitido dentro de la tolerancia de 1.25m a escala 1:5000 con el estadístico del
RMSEz, si se considera el límite a escala 1:8000 (2m), los productos generados en Z/I Imagine son
los únicos no aceptados; cabe mencionar que no se debe utilizar la media o la media absoluta pues
214
no es un buen estimador estadístico al observar que algunos productos entrarían en los límites
1:5000.
en un área menor a 500km2 , son suficientes para tener una confiabilidad del 99% (ver tabla 30).
Z/I Imagine); se tiene que, si se considera las precisiones en cada proceso (ver apartados 3.5.1.2 y
En importante agregar que la revisión de la calidad radiométrica es de una manera visual y todos
los ortomosaicos no presentan cambios extremos en sus niveles digitales; sin embargo, aún
persisten los desplazamientos en: edificaciones, zonas de alta pendiente y extremos de las ortofotos
5.2. Recomendaciones
tecnologías (LIDAR, Vehículos Aéreos no Tripulados), con el fin de obtener un MDT que mejore
de estandarizar los procesos y localizar de una manera más efectiva los posibles errores que puedan
ocurrir; así también, se debe determinar los umbrales de cada proceso y documentarlos; de modo
que se tenga un control más exhaustivo. Vale destacar que sería de gran importancia la aplicació n
de flujogramas, manuales, documentos con las directrices técnicas del área fotogramétr ica,
En el caso de aumentar el número de GCP, estos deberían realizarse sobre zonas de alta
pendiente, zonas de fiscalización, áreas urbanas, laderas, extremos del proyecto y control vertical
minutos; así también, evitar objetos cercanos que puedan afectar la señal del equipo por ejemplo:
la orientación exterior; así como, otros modelos matemáticos para compensar un proyecto.
la relación entre RMSEx y RMSEy son iguales con pruebas de hipótesis (Prueba para diferenc ias
6. BIBLIOGRAFÍA
Photogrammetria.
Ackermann, F. (1990). Kinematic GPS Camera Positioning for aerial triangulation. (págs. 1155 -
Agromundo. (27 de Mayo de 2018). Rg3V.LVN. Obtenido de 2.Modelos digitales del terreno
mde
Society of Photogrammetry.
Angulo. (2014). Generación de cartografía básica a detalle mediante una metodología de toma con
Ariza, F., & Atkinson, A. (2006). Metodologías de control posicional: Visión General y Análisis
Ariza, F., Emerson, X., Chicaiza, E., & Buenaño, X. (2017). Métodos de evaluación de la calidad
Arranz, J. (2010). Modelos digitales del terreno. E.T.S. de ingenieros en topografía, geodesia y
ASPRS. (2014). ASPRS Positional Accuracy Standards for Digital Geospatial Data. En
Photogrmmetric Engineering & Remote Sensing (págs. A1 - A26). American Society for
Berrocoso, M., Sánchez, A., Pérez, A., Jigena, B., Fernández, J., & Páez, R. (s.f.). Sistemas GNSS:
http://lagc.uca.es/web_lagc/docs/curso_rap/Presentacion_II.pdf
Bianco, M., & Martínez, J. (24 de Junio de 2004). Inferencia estadística. Obtenido de Intervalos
de confianza :
http://www.dm.uba.ar/materias/probabilidades_estadistica_C/2005/1/PyEC014.pdf
218
Bostelmann, J., Breitkopf, U., & Heipke, C. (2017). Optimizing the distribution of tie points for
Heipke, The international Archives of the photogrammetry, Remote sensing and Spatial
Buill, Nuñez, & Rodriguez. (2003). Fotogrametría analítica. Barcelona: Centre de Publicacio nes
http://www.cartesia.org:
http://www.cartesia.org/data/apuntes/fotogrametria/Introduccion_a_la_Fotogrametria.pdf
CLIRSEN. (2011). Generación de Geoinformación para la gestión del territorio a nivel nacional
Cliver, E. (1997). A history of the use of photogrammetry at the Historic American Buildings
Coelho, & Nunes. (2007). Fotogrametría Digital. Río de Janeiro: Universidade do Estado do Rio
de Janeiro.
219
Comer, R., Kinn, G., Light, D., & Mondello, C. (1998). Talking digital. Photogrammetr ic
Crouse R. & Associates. (17 de Abril de 2018). Quality aerial photography and Remote Sensing
http://www.richardcrouse.com/pdf/CalibProtocol_DMCIIe230-23522.pdf
Cruz, E. (2011). El uso del GPS en restitución fotogramétrica y de las poligonales de referencias
De la Torre, A. (18 de Agosto de 2018). Escuela Técnica Superior de Ingeniería Civil, Alternativa
Fotogrametría Digita l:
http://coello.ujaen.es/asignaturas/fotodigital/descargas/fd_tema3.pdf
Delta 95. (18 de Junio de 2018). Topografía cartografía medición de fincas. Obtenido de
http://www.delta95.com/web/imatges%20index_web/mediostecnicos/restitdigit.html
220
Durón, J. (2007). Modelos digitales de elevación para el continuo de elevación mexicano. Obtenido
de INEGI:
http://mapserver.inegi.org.mx/eventos/cng2007/cng2007/extcemjavierduron.pdf
https://www.egnos-portal.eu/discover-egnos/about-egnos/what-gnss
básica del instituto geográfico Agustín Codazzi. Bogotá, Colombia.: Universidad distrita l
Falkner, E., & Morgan, D. (2002). Aerial Mapping. Methods and Applications (Second Edition).
Lewis Publisher.
Farjas, M. (2006). Aplicaciones topográficas del GPS. Obtenido de OpenCourseWare Univers idad
fotogrametria/topografia- ii/Teoria_GPS_Tema_12.pdf
221
FGDC. (1998). Geospatial Positioning Accuracy Standards. Part 3: National Standard for Spatial
García, P., & Ortíz, E. (26 de Marzo de 2018). Influencia del GSD en la orientación de un bloque
fotogramétrico. Obtenido de
https://addi.ehu.es/bitstream/10810/9005/1/Influencia%20del%20GSD%20en%20la%20o
rientaci%C3%B3n%20de%20un%20bloque%20fotogram%C3%A9trico.pdf
GIS & Remote Sensing Service. (27 de Mayo de 2018). Vector and Raster Mapping. Obtenido de
Services: https://www.overallvision.com.au/our-services/gis-remote-sensing-services/
Gis Ibérica. (17 de Junio de 2018). Paralaje / Estereoscópica. Obtenido de Equipos de precisión:
http://www.gisiberica.com/estereoscopos/PARALAJE%20ESTEREOSC%D3PICA.htm
Global Market. (02 de Junio de 2018). Technology Hardware, Storage and Peripherals Company
https://www.bloomberg.com/research/stocks/private/snapshot.asp?privcapid=3235439
http://ocw.usal.es/ensenanzas-tecnicas/procesamiento-avanzado-de-imagenes-
digitales/contenidos/Tema4.pdf
222
González, P. (2015). Levantamiento mediante GPS de una red de puntos establezidos para
Gualcamayo: http://www.gpacsrl.com/index.php?pr=Gualcamayo
GPS.gov. (14 de Agosto de 2018). National Coordination Office for Space-Based Positioning,
https://www.gps.gov/spanish.php
Heipke, C. (1997). Automation of interior, relative, and absolute orientation. ISPRS Journal of
Mancha.
Hexagon. (30 de Mayo de 2018). Erdas Imagine: World-Class Remote Sensing Software. Obtenido
Hinojosa, A., Rodríguez, V., Munguía, L., & Meillón, O. (Mayo de 27 de 2011). Scielo. Obtenido
33222011000100003&lng=es&tlng=es.
Hinsken, L., & Miller, S. (2002). Triangulation fo the LH systems’ ADS40 using ORIMA GPS/IMU.
Konstanz: LH Systems.
Hussain, M., & Bethel, J. (2004). Project and Mission Planning in Manual of Photogrammetry,
Fitth edition. Maryland, USA: American Society for Photogrammetry and Remote Sensing.
Idainature. (17 de Junio de 2018). Satélites Sentinel 2 para teledetección en agricultura. Obtenido
sentinel-2-para-teledeteccion-en-agricultura/
ortofotos digitales a escala 1:5000 del área rural del cantón Ambato. Quito: Divis ió n
Cartográfica.
224
ortofotos digitales de la cuenca baja del río Guayas. Quito: Gestión mercadotécnia y
comunicación.
IGM. (2010). Altura de vuelo en vuelos fotogramétricos. Quito, Ecuador: Instituto Geográfico
Militar.
IGM. (2012). Archivos APLANIX, fecha de toma. Quio: Instituto Geográfico Militar.
IGM. (18 de Marzo de 2013). Geoportal IGM. Obtenido de Evolución de la geodésia en el Ecuador:
http://www.geoportaligm.gob.ec/wordpress/?page_id=248
IGM. (2013). Metodología utilizada para la generación de cartografía básica del Ecuador
IGM. (2017). Protocolo de fiscalización para proyectos de generación de cartografía base con
Militar.
225
IGM. (12 de Junio de 2017). RED GNSS de monitoreo continuo del Ecuador. Obtenido de Ficha
http://www.geoportaligm.gob.ec/wordpress/?wpfb_dl=49
IGM. (21 de Febrero de 2018). Ecuador_XP. Camara. Quito, Pichincha, Ecuador: Instituto
Geográfico Militar.
http://www.geoportaligm.gob.ec/visor_regme/
IGN. (2011). Especificaciones técnicas para la producción de cartografía básica escala 1:1000.
Instituto de Información Territorial del Estado de Jalisco. (2008). Norma técnica para la
ISPRS. (13 de Marzo de 2018). International Society for Photogrammetry and Remote Sensing.
Kimera. (25 de Julio de 2018). Biblioteca virtual. Obtenido de Aspectos Financie ros :
http://www.kimera.com/RLCF/RECURSOS/BIBLIOTECA%20CAFETERA/N%20-
%20ASPECTOS%20FINANCIEROS/desviacion%20estandar.pdf
Kirby, E. (2007). Metodología para el cálculo de las anomalías de Bouger completa. Sangolquí:
Kirby, E. (2009). Metodología para navegación en Arcgis a través del uso del GPS. Caso de
Kirby, E. (2012). Cálculo del número de puntos GPS para la corrección geométrica y generación
Lara, A., Cerpa, G., Rodríguez, M., & Núñez, H. (2006). Física para bachillerato cinemática. En
de Valencia.
Liba, N., Järve, I., & Rand, M. (2013). Quality of orthophoto mosaics made with different methods.
Lisintuña, J., & Molina, G. (2017). Análisis y aplicación de ingeniería inversa en objetos con
Mancero, H., Toctaguano, D., Tacuri, C., Kirby, E., & Tierra, A. (2015). Evaluación de modelos
Milkhail, E., Bethel, J., & McGlone, J. (2001). Introduction to modern photogrammetry. John
Molina, S., Redondo, F., Zamorano, M., Hermosilla, A., & Caturla, C. (2010). Producción de
Sevilla.
228
Muñoz, P. (2004). Apoyo aéreo cinemático y aerotriangulación digital frente a los sistemas de
Madrid.
Murray, S., & Larry, J. (2009). Estadística Schaum. México: The McGraw - Hill Companies.
NGA/NASA. (16 de Junio de 2018). EGM96 Earth Gravitational Model. Obtenido de Geoid
info.nga.mil/GandG/wgs84/gravitymod/egm96/egm96.html
Argentina:
http://webcache.googleusercontent.com/search?q=cache:2h2NAxWZzcwJ:ftp://ftp.unsj.e
du.ar/agrimensura/Fotogrametria/Unidad6/Tema_6_Concepto_de_ortofoto_y_ortofotopla
no.doc+&cd=3&hl=es-419&ct=clnk&gl=ec
229
Otero, J., & Sevilla, M. (1989). Modelo matemático para el ajuste simultanea mínimos cuadrados
http://agrimensores.org.uy/descargas/Evaluacion_de_Calidad_Posicional_Producto.pdf
PIX4D. (05 de Junio de 2018). Getting GCPs on the field or through other Sources. Obtenido de
Starting-a-Project-4-Getting-GCPs- in-the-Field-or-Through-Other-Sources-optional-but-
recommended-Ram%F3n
230
Portilla, O., Tierra, A., & Godoy, L. (2016). Variables que influyen en el rastreo óptimo GNSS
Pozo, M. (2002). Desarrollo de equipos de fotogrametría: Explotación del dato digital. Santander,
http://dle.rae.es/?id=XfDUOlU
http://fcf.unse.edu.ar/archivos/lpr/pdf/p14.PDF
Sáez, A. (03 de Junio de 2012). Apuntes de Estadística para Ingenieros. Obtenido de Departamento
http://www4.ujaen.es/~ajsaez/recursos/EstadisticaIngenieros.pdf
Sato, & Silva. (2004). Brazilian quality control systems for digital photogrammetry mapping
Soria, B. (2014). Metodología para ejecución y control de calidad para ortofotos, mediante
combinación de sensor digital con lidar en zonas rurales. En B. M. Soria, Metodología para
ejecución y control de calidad para ortofotos, mediante combinación de sensor digital con
Spiller, R. (1999). Z/I Imaging: A New System Provider for Photogrammetry and GIS. Oberkochen,
3-540-92725-9/espanol/foto.pdf
Tierra, A. (2009). Evaluación del EGM08 y EGM96 en el Ecuador a partir de datos de GPS y
http://www.topoequipos.com/topoequipos2.0/fotogrametraa-digital/fotogrametria-digital
232
Tsarovski, S. (2015). Accuracy Evaluation and Quality Control of Digital Orthomap-Sheets. Sofia,
Bulgaria: From the Wisdom of the Ages to the Challenges of the Modern World.
15.pdf
Vexcel Imaging. (2010). Calibration Report, Geometric Calibration, UltraCam Xp. Graz, Austria:
Microsoft Company.
Vizuete. (2012). Análisis costo beneficio de la utilización del método alternativo para la
y proyectos.