TEMA 9

LA RECTIFICACIÓN DE FOTOGRAMAS
9.1. Introducción
9.2. Fundamento de la rectificación
9.3. Desplazamiento debido al relieve
9.4. Métodos de rectificación
9.5. De la rectificación a la ortofotogrametría
Bibliografía

9.1. Introducción

Por rectificación se conoce el proceso fotogramétrico mediante el cual se puede transformar una
imagen en una proyección del terreno, por medio del enderezamiento del eje del fotograma. La
rectificación corrige los desplazamientos existentes en la imagen fotográfica original producidos por la
inclinación del eje de la toma.

Su desarrollo se ha debido principalmente a las aplicaciones en las fotografías aéreas como un

método muy rápido para actualizar cartografía. En la actualidad puede aplicarse a un gran número de
imágenes (radar, multiespectrales, digitales, etc.). Es una técnica que permite obtener pseudomapas de
zonas en las que la fotografía clásica no puede hacerlo.

N N´

Figura 1. Enderezamiento de un fotogramas.

FOTOGRAMETRÍA III

La imagen rectificada debe tener las características geométricas de una proyección ortogonal del
objeto fotografiado sobre un cierto plano y a una determinada escala. Se presupone que el objeto es
plano y vertical.

Para poder realizarse una rectificación se han de cumplir dos condiciones:

ƒ Ángulo de inclinación de la toma tiene que ser muy pequeño (3º).

ƒ El terreno tiene que ser prácticamente horizontal y llano.

9.2. Fundamento de la rectificación

Para el caso de una fotografía aérea vertical, el valor del ángulo i (Fig. 1) suele ser del orden de 1g a
3g , y excepcionalmente puede tomar valores de 4g. Las fórmulas que definen la reproducción de la
perspectiva central son las ecuaciones de colinealidad:

m11 (X − X o ) + m12 (Y − Yo ) + m13 (Z − Z o )

x ′ = −c
m 31 (X − X o ) + m 32 (Y − Yo ) + m 33 (Z − Z o )

m 21 (X − X o ) + m 22 (Y − Yo ) + m 23 (Z − Z o )
y′ = − c
m 31 (X − X o ) + m 32 (Y − Yo ) + m 33 (Z − Z o )

Estas ecuaciones contienen seis incógnitas y relacionan los parámetros de orientación externa, lastres
traslaciones ω, ϕ, κ , función de los parámetros a11, a12,.., a13, las coordenadas del centro de proyección
X0, Y0, Z0 y las coordenadas terreno de los puntos de apoyo X, Y, Z.

En el caso de que el terreno sea horizontal, condición que implica que la coordenada Z sea constante,
el problema puede tener una solución única, dado que el conocimiento de las direcciones que nos
permite conocer el sistema de ecuaciones anteriores, tendría distintas soluciones al proyectarse la
perspectiva central sobre un plano y dentro de los infinitos planos solución habrá uno que nos de la
relación de escala buscada.

La determinación de los ocho coeficientes que permiten obtener la solución de las relaciones entre
dos planos, la podemos expresar por las ecuaciones de la transformación bidimensional proyectiva,
siendo las expresiones:

a 1 x ′ + b1 y ′ + c1
X=
a 3 x ′ + b 3 y′ + 1

a 2 x ′ + b 2 y′ + c 2
Y=
a 3 x ′ + b 3 y′ + 1

donde x´, y´ son las coordenadas imagen, X e Y las coordenadas del plano y los coeficientes a, b y c son
proporcionales a las ecuaciones generales de la perspectiva central.

Si se conocen las coordenadas de un punto sobre el plano y las correspondientes sobre la fotografía
pueden establecerse dos ecuaciones, en las cuales, los coeficientes a, b y c son las incógnitas; con
cuatro puntos se formarían un total de ocho ecuaciones que permitirían obtener una solución única del
sistema de ecuaciones que se formara.

178
 TEMA 9: RECTIFICACIÓN DE FOTOGRAMAS

9.3. Desplazamiento debido al relieve

Se ha establecido como requisito indispensable para poder realizar correctamente una rectificación el
disponer de un terreno horizontal y llano, por supuesto siempre dentro de unos ciertos límites de
tolerancia, que vienen impuestos por los errores admisibles en la proyección, debido a los
desplazamientos del relieve.

El desplazamiento del relieve se ilustra en la figura 2, donde se representa una toma fotográfica
vertical tomada a un altura media de vuelo de valor H. La distancia principal de la cámara es c y O es el
centro de proyección de la toma. La imagen del punto P sobre el terreno de altura h queda representado
sobre el plano de la imagen por p. Un punto imaginario P´ es localizado en la vertical del punto P,
sobre el plano de referencia (terreno), y sus correspondiente posición sobre la imagen viene dada por
p´.

En la figura se observa que el desplazamiento del punto P producido por el relieve h viene dado por
en la imagen por el segmento rr´= (∆r) , siendo un desplazamiento radial desde el punto principal.
Lógicamente cuanto más alejado estemos del centro de la fotografía mayor será la magnitud de dicho
desplazamiento.

O
∆r

c
p´ p

r´
r
H

O´ P´
h ∆R
R

Figura 2. Desplazamiento en el plano de la imagen debido al relieve.

De la figura se deduce las expresiones siguientes:

R H−h
= ⇒ Rc = r (H - h)
r c

R H
= ⇒ Rc = r ′H
r′ c

Igualando las dos expresiones anteriores tendremos:

r (H - h) = r ′H ⇒ H(r - r ′) = rh

179
FOTOGRAMETRÍA III

Sustituyendo en la expresión ∆r = r-r´ tendremos:

rh
∆r =
H
siendo ∆r el desplazamiento producido en la imagen por el relieve del terreno, h la altura del objeto, H
la altura de vuelo y r la distancia radial sobre la fotografía desde el punto principal al punto desplazado.

Esta ecuación expresa el desplazamiento del relieve para fotogramas verticales. Observando la
ecuación se deduce que el desplazamiento debido al relieve aumenta con el incremento de la distancia
radial y también incrementa con el aumento de la altura de los objetos.

9.4 Métodos de rectificación

La rectificación puede ser analítica o gráfica, teniendo en ambos casos la obtención de un documento
cartográfico, enderezando un solo fotograma. La rectificación puede realizarse punto a punto o de
forma global.

Básicamente hay cuatro métodos que permiten corregir el desplazamiento introducido sobre una
fotografía por la inclinación del eje de la toma, estos son:

ƒ Técnicas gráficas
ƒ Transformaciones matemáticas
ƒ Métodos mecánicos de retroproyección
ƒ Transformaciones ópticas.

La rectificación basada en transformaciones matemáticas se basa en la utilización de la

transformación bidimensional proyectiva.

a 1 x ′ + b1 y ′ + c1 a 2 x ′ + b 2 y′ + c 2
X= Y=
a 3 x ′ + b 3 y′ + 1 a 3 x ′ + b 3 y′ + 1

Las transformaciones ópticas son las que han estado siendo más utilizada en el pasado. Los otros
métodos básicamente se basan en la rectificación punto a punto, mientras que por procedimientos
ópticos la rectificación se efectuaba por medio de un rectificador y ésta era global. Con este
instrumento todos los puntos eran transformados instantáneamente

Figura 3. Rectificación gráfica.

180
 TEMA 9: RECTIFICACIÓN DE FOTOGRAMAS

Figura 4. Rectificación óptica.

Figura 5. Izquierda, rectificador Zeiss SEG-5. Derecha, efecto de los ajustes en el rectificador.

181
FOTOGRAMETRÍA III

9.5. De la rectificación a la ortofotogrametría

Se ha visto en anteriores apartados que la rectificación de fotogramas se puede realizar siempre y

cuando se partan de imágenes adquiridas con inclinaciones muy pequeñas y considerando el terreno
practicante horizontal y llano o asumiendo que los desniveles existentes no superan los errores
admisibles de la escala de proyección.

Admitidas unas tolerancias para las desviaciones producidas por el relieve, si estas son sobrepasadas,
la rectificación daría valores errónea Para que fuera válida habría que rectificar uno a uno todos los
rayos del haz perspectivo.

Ante la imposibilidad evidente de llevar a cabo esta tarea, se puede hacer una rectificación de un
conjunto de puntos (línea o superficie). Esta rectificación llamada rectificación diferencial es la idea
general de la ortoproyección.

La rectificación de imágenes satélite (Spot, Landsat, etc.) es posible en mayor número de casos que
las imágenes aéreas debido a que la altura de vuelo H es muy grande, resultando que los valores de los
desplazamientos introducidos por el relieve (∆r) son muy pequeños, haciendo posible que puedan en
mayor medida entrar en tolerancia.

En la actualidad la mayoría de las rectificaciones que se realizan esta relacionadas con la

rectificación de elementos arquitectónicos, como por ejemplo el alzado de una fachada donde no
existen grandes variaciones de distancias al punto de toma (profundidades), de los planos que forman el
mismo.

Bibliografía

ƒ American Society of Photogrammetry. Manual of Photogrammetry, ASPRS, Cuarta Edición

(1980). ISBN 0-937294-01-2.
ƒ Clavo, D., Fotogrametría II. (EPM - EUITT).
ƒ Curso sobre técnicas fotogramétricas y sus aplicaciones, CEDEX (MOPTMA), Madrid (1994).
ƒ López, S., Fotogrametría, Egraf S.A., Madrid (1980). ISBN 84-300-2559-6.
ƒ Krauss, K., Photogrammetry, Dümmler/Bonn (1993), 3-427-78684-6.
ƒ Wolf, P.R. y Dewitt, B.A., Elements of Photogrammetry with Applications in GIS, McGraw-Hill,
EEUU, Tercera Edición (2000). ISBN 0-07-292454-3.

182
 TEMA 10
ORTOFOTOGRAMETRÍA DIGITAL
10.1. Introducción
10.2. Rectificación de imágenes digitales
10.3. Métodos de rectificación de imágenes digitales
10.3.1. Rectificación polinómica
10.3.2. Rectificación proyectiva
10.3.3. Rectificación diferencial
10.4. Esquema de la generación de ortoimágenes
10.5. Técnicas de remuestreo
10.6. Calidad de la imagen ortofotográfica
10.7. Correcciones radiométricas
10.8. Confección de mosaicos
10.9. Estereo-ortofotogrametría digital
10.10. Propiedades de las ortofotos digitales
Bibliografía

10.1. Introducción

La representación fotográfica del terreno debido a su mayor riqueza informativa resulta ventajosa en
numerosas aplicaciones respecto a los mapas gráficos. La ortofotografía es un documento fotográfico
que permite aprovechar el contenido informativo de la fotografía aérea y realizar las mismas medidas
que en un plano gráfico.

Los principios en los que se basa el desarrollo de la ortofotografía, fueron enunciados por Lacmann
en el año 1931 y más tarde por Ferber (1933), pero hasta la construcción del Orthophotoscope, en 1955
por Bean, se puede decir que no empezó el desarrollo y aplicaciones de esta técnica cartográfica.

En el caso ideal teórico de ser el eje de la toma rigurosamente vertical y el terreno llano y horizontal,
la fotografía aérea coincidirá con la ortofotografía. Sí el terreno es llano y horizontal y el eje de la toma
forma un pequeño ángulo con la vertical, se corrigen las deformaciones de la imagen original a través
de la rectificación fotográfica convencional. En el caso general de disponer de una imagen
correspondiente a una zona de terreno con relieve y tomada con el eje inclinado, para obtener una
ortofotografía es preciso aplicar la técnica de la rectificación diferencial.
FOTOGRAMETRÍA III

Esta técnica corrige las deformaciones de la imagen, causadas, tanto por los desniveles del terreno
como por la falta de verticalidad del eje óptico de la cámara en el momento de la toma, transformando
geométricamente zonas elementales de la fotografía aérea en ortofotografía.

Podemos definir ortoproyección como el método fotogramétrico de obtener fotomapas precisos,

corregidos de las deformaciones producidas en la imagen fotográfica por el relieve y la inclinación de
la fotografía, mediante la rectificación diferencial de la fotografía original.

10.2. Rectificación de imágenes digitales

La rectificación digital presenta importantes ventajas respecto a los procedimientos óptico-

mecánicos. La calidad de la imagen, al no estar en estos procedimientos supeditada a complejos
dispositivos de enfoque que afecten a la nitidez fotográfica en el plano de reproducción, no constituyen
ninguna restricción al sistema, tampoco imponen condicionantes importantes al sistema.

Se mantienen las restricciones en cuanto a desplazamiento debido al relieve, tanto según las
características del terreno como de los objetos en el caso de fotogrametría no topográfica.

La rectificación diferencial de la fotografía consiste en enderezar y poner en escala, pequeñas

unidades geométricas del terreno, que perfectamente ensambladas dan una imagen fotográfica del
terreno.

Este nuevo documento fotográfico se denomina ortofotografía, y es tan preciso como un mapa.

Para que la transformación de la proyección central de la fotografía a la ortogonal del terreno, sea
teóricamente exacta, sería necesario que las unidades geométricas a rectificar fuesen un punto. Esto en
la práctica es imposible, utilizándose unidades geométricas de reducido tamaño (pequeñas zonas de la
imagen, que en el caso digital será el píxel).

Podemos diferenciar dos tipos de ortofotografías:

ƒ Ortofotografía óptica: A partir de la diapositiva de la fotografía original se transmite la imagen

ópticamente, en pequeñas unidades geométricas, al negativo final. El producto final es una
imagen fotográfica.

ƒ Ortofotografía digital: Se parte de imágenes generadas en forma digital o imágenes fotográficas

que previamente se han escaneado. Estas imágenes están dividas en elementos rectangulares
denominados pixeles. La corrección de los errores debidos al relieve y a la falta de verticalidad
del eje de torna, hasta obtener una nueva imagen digital se realizan analíticamente. El resultado es
una nueva imagen digital corregida, que puede ser reproducida fotográficamente.

En este tema nos centraremos en la ortofotografía digital por ser el método actualmente más
extendido.

10.3. Métodos de rectificación de imágenes digitales

Podemos obtener imágenes digitales del terreno por medio de sensores aerotransportados, cámara
digitales o bien transformando la imagen fotográfica en digital por medio de un escáner. Los sensores
no forman la imagen según una proyección cónica, y además, para el caso de plataformas espaciales
(satélites) dada la altura de vuelo los errores debidos al relieve son mínimos.

184
 TEMA 10: ORTOFOTOGRAMETRÍA DIGITAL

Los métodos más frecuentes de rectificación de imágenes digitales son:

ƒ Rectificación polinómica.
ƒ Rectificación proyectiva.
ƒ Rectificación diferencial.

10.3.1. Rectificación polinómica.

La transformación entre la imagen original y la rectificada se realiza mediante el empleo de

polinomios.

x = x ′ T Ay′
y = x ′ T By′
donde:

) x, y son las coordenadas de la imagen original.

) x´, y´ son las coordenadas de la imagen rectificada.

x ′ T = (1, x ′, x ′ 2 , x ′ 3 ,.........) ; y′ T = (1, y′, y′ 2 , y′ 3 ,.........)

) A, B son las matrices de los coeficientes de los polinomios

a 00 a 01 a 02 K b 00 b 01 b 02 K
   
 a 10 a 11 a 12 K  b10 b11 b12 K
A=  B= 
a 20 a 21 a 22 K b 20 b 21 b 22 K
K K K K K K K K
 

Este método corrige las distorsiones de la imagen a base a una red de puntos de control. El orden del
polinomio utilizado depende del número de puntos de control disponible, siendo mayor la precisión al
aumentar el número de éstos. Este método es completamente independiente de la geometría del sensor,
por lo que puede ser empleado tanto para fotografías aéreas como para imágenes satélite. Debido a la
transformación polinómica es posible el cambio de escala y la rotación de la imagen de la forma que
mejor se ajuste a los puntos de referencia.

Aunque el método polinómico es muy sencillo de utilizar es necesario tener en cuenta que puede
originar problemas y errores en la imagen transformada, debido a que no corrige adecuadamente los
desplazamiento debido al relieve y además no considera las características específicos del sistema de
captura de las imágenes. Es importante tener en cuenta, que el incremento del número de orden del
polinomio reduce considerablemente los residuos del ajuste pero debido a ondulaciones en el terreno
puede ocasionar que aparezcan grandes errores entre los puntos de control. No obstante, este método
presenta la gran ventaja de que corrige en un solo proceso todas las distorsiones de la imagen (debido a
la geometría del sensor, curvatura de la Tierra, etc.).

Este es un método muy empleado en la rectificación de imágenes satélite debido a que la geometría
y distorsiones de las mismas a menudo son complejas de modelizar. Además en este tipo de imágenes,
el desplazamiento debido al relieve es muy pequeño debido a las escasas diferencias de cota en relación
con la que se toma la imagen y afecta poco a los resultados finales.

185
FOTOGRAMETRÍA III

10.3.2. Transformación proyectiva

La transformación proyectiva plantea la relación existente entre dos planos (Fig 1). Esta relación es
definida por 8 parámetros que pueden calcularse a partir de cuatro puntos de control y sus
correspondientes coordenadas imagen. No es necesario el cálculo de los elementos de orientación
interna y externa puesto que están incluidos implícitamente en los coeficientes de la transformación.
Este método es de uso frecuente para rectificación de fotografías aéreas de zonas llanas o de imágenes
de fachadas de edificios.

a 1 x ′ + b1 y ′ + c1
X=
a 3 x ′ + b 3 y′ + 1

a 2 x ′ + b 2 y′ + c 2
Y=
a 3 x ′ + b 3 y′ + 1

2´

1´
3´

4´

2
1

Figura 1. Transformación proyectiva.

La expresión de la transformación proyectiva puede ser modificada para la toma de una imagen
desde un satélite. Los sensores instalados en los satélites, como por ejemplo el SPOT, realizan la
captura mediante un barrido de líneas, por tanto, su geometría es de tipo perspectiva lineal, no
existiendo un punto de vista sino una línea de vista. Cada línea de barrido tiene un centro de perspectiva
diferente que se puede asumir que están en una línea recta (órbita). Si se considera el eje y en la
dirección de vuelo y x representa el pixel en la línea de barrido, la expresión de la transformación
perspectiva se modifica de la forma:

a 1 x ′ + b1 y ′ + c1
X=
a 3 x ′ + b 3 y′ + 1

Y = a 2 x ′ + b 2 y′ + c 2

La utilización de este método en el tratamiento de imágenes satélites es compleja ya que debido a la

curvatura terrestre, el definir la superficie terrestre como un plano lleva a producir errores en la imagen

186
 TEMA 10: ORTOFOTOGRAMETRÍA DIGITAL

rectificada, incluso para terrenos llanos. No obstante, este método es valido para la rectificación de
imágenes de sensores aerotransportados.

10.3.3. Rectificación diferencial

Como sabemos la imagen digital se compone de una matriz de elementos rectangulares de pequeño
tamaño, denominados pixel. De cada uno de estos pixel conocemos su tamaño, su posición en la
imagen y el valor en la escala de grises asociado.

En imágenes digitales está claro que el elemento diferencial a rectificar es el pixel, puesto que es la
unidad mínima de que se compone.

El proceso se realiza en modo diferido, obteniendo en primer lugar el modelo digital del terreno
(datos altimétricos), procediéndose a continuación a realizar la rectificación diferencial.

Para la realización de un proceso de rectificación diferencial es necesario disponer de un modelo

digital de elevaciones almacenado de la misma forma que la imagen, es decir, con estructura raster de
una matriz de elevaciones. Este modelo debe estar referido a un sistema de proyección de coordenadas
que será en el que se obtenga la imagen rectificada.

Para entender mejor el proceso imaginemos una imagen digital y su transformación al terreno, (Fig.
2), suponiendo que cada pixel es plano se transforma en un trapecio. El conjunto de estas piezas en el
terreno perfectamente ensambladas constituiría la ortofotografía. El elemento diferencial rectificado ha
sido el pixel.

y
x
O
c

Transformación directa
Z

Y
P(X,Y,Z)

Figura 2. Transformación directa (fotografía mapa)

El caso contrario (Fig. 3), una imagen digital que sea proyección ortogonal del terreno, por ejemplo
un mapa escaneado. Si transformamos los píxeles del mapa a la proyección de la fotografía,
obtendríamos una serie de trapecios que ensamblados componen la imagen fotográfica.

187
FOTOGRAMETRÍA III

x
O
c

Transformación inversa Z

Y P(X,Y,Z)

Figura 3. Transformación inversa (mapa fotografía)

La relación entre ambas proyecciones, la ortogonal del terreno y la cónica de la fotografía viene
dada por las ecuaciones de colinealidad. El objetivo de la rectificación de las imágenes digitales es
asignar a cada uno de los elementos de la matriz del modelo digital el valor de gris correspondiente.
Para la determinación del valor de gris se utilizan estas ecuaciones. Así las coordenadas terreno X,Y,Z
de los puntos que integran el modelo digital son transformados en coordenadas imagen x, y mediante la
utilización de las ecuaciones de colinearidad.

m11 (X − X o ) + m12 (Y − Yo ) + m13 (Z − Z o )

x = −c
m 31 (X − X o ) + m 32 (Y − Yo ) + m 33 (Z − Z o )

m 21 (X − X o ) + m 22 (Y − Yo ) + m 23 (Z − Z o )
y = −c
m 31 (X − X o ) + m 32 (Y − Yo ) + m 33 (Z − Z o )

Estas ecuaciones expresan las coordenadas (x, y) en la imagen en función de las coordenadas (X, Y,
Z) del terreno y los parámetros de orientación exterior (X0, Y0, Z0, ω, ϕ, κ) . Si realizamos la orientación
interna y exterior de la fotografía y además conocemos la Z del punto terreno, podemos realizar las
anteriores transformaciones:

ƒ Proyección de la fotografía (x, y) a mapa (X, Y -terreno-).

ƒ Proyección de mapa (X, Y -terreno-) a fotografía (x, y).

En la ortoproyección digital se utiliza la transformación inversa. A cada píxel de la imagen en la

proyección del mapa, le vamos a asignar el valor de la escala de grises en la proyección cartográfica.

La posición de la imagen digital en el terreno la conocemos y por tanto las coordenadas de cada
píxel, ya que previamente habremos definido una distribución de ortos con un determinado formato
para cubrir la zona. A continuación habremos elegido una fotografía que cubra la orto a realizar.

Para realizar las transformación es necesario disponer de los siguientes parámetros:

ƒ Orientación interna de la cámara: coordenadas del punto principal de la cámara, distancia

principal, marcas fiduciales y distorsión de la lente.
ƒ Orientación externa: posición de la cámara en el espacio; coordenadas de la posición del centro de
proyección (X0, Y0, Z0) y los tres giros de la cámara (ω, ϕ, κ).

188
 TEMA 10: ORTOFOTOGRAMETRÍA DIGITAL

ƒ El tamaño del pixel de la imagen digital en unidades cámara (mm).

ƒ El tamaño de la malla del modelo digital en unidades terreno (m).
ƒ Las coordenadas de referencia del modelo digital en una determinada proyección cartográfica
(normalmente las coordenadas del punto superior izquierda).

10.4. Esquema de la generación de ortoimágenes

Se parte de la premisa de disponer de los datos vistos en el apartado anterior (imagen digital,
certificado de calibración, modelo digital, etc.).

En la figura 4 se muestran los pasos a seguir en la generación de una ortofotografía digital. Se parte
de una imagen “vacía”, formada por un conjunto de píxeles sin asignación de valor de tonalidad (Nd),
caracterizados por estar referidos al sistema de referencia terreno y por poseer un cierto tamaño físico,
que constituye la resolución de la ortoimagen. Aplicando en sentido inverso (transformación inversa)
las ecuaciones de colinealidad, se rellenará dicho valor de tonalidad.

xp

yp

5. Determinación del valor de la escala de grises

y
4. Obtención de la posición del punto sobre la imagen

3. Condición de colinealidad

2. Cordenada Z interpolada del modelo

1. Posición píxel terreno Y

6. Asignación del valor al pixel inicial correspondiente en la ortoimagen

Figura 4. Generación de una ortofotografía digital.

Supongamos que se parte de la posición de cualquier punto sobre el terreno (pixel) de coordenadas
X, Y conocidas; por interpolación del modelo digital del terreno podremos obtener la coordenada Z de
dicho punto. Con estas coordenadas, los parámetros de orientación externa y las ecuaciones de
colinealidad podrá obtenerse analíticamente la posición del punto sobre la imagen (x, y), como si
realmente se reconstruyera el rayo luminoso que lo impresionó. A partir de la orientación interna
efectuada se identifica dicha posición en la imagen y su correspondiente tonalidad, bien en modo de

189
FOTOGRAMETRÍA III

escala de grises, o bien en modo color (RGB); este valor se asignará a la posición (X, Y) de partida de la
ortoimagen.

Está claro que en el proceso de rectificación que se lleva a cabo, las posiciones de los píxeles
proyectados no coinciden con la posición de los píxeles originales, siendo necesario combinar de
alguna forma la información espectral (Nd) de los píxeles vecinos para obtener el valor final deseado de
la tonalidad, empleándose habitualmente las técnicas de remuestreo digital, asegurando así que todos
los pixeles de la ortofotografía resultante, tengan su nivel de gris.

10.5. Técnicas de remuestreo

Como se ha visto anteriormente hay que determinar el nivel de gris que le corresponde a cada píxel,
debido a que la posición proyectada no tiene porqué coincidir con la posición original del mismo ni con
su valor entero. Al realizar la transformación espacial se tiene que los valores de las coordenadas de la
imagen no coincidirán con los de la imagen original de forma que habrá puntos que no quedarán
definido por sus niveles digitales (Fig. 5). Por tanto, es necesario interpolar el valor de gris
conociéndose esta operación con el nombre de remuestreo (resampling).

Imagen original Imagen resultado

Figura 5. Remuestreo de imágenes.

Por medio de este método (remuestreo) hacemos varios muestreos sobre la imagen inicial. Cada vez
que se efectúan cambios geométricos en la imagen digital, surgen varios errores de forma inevitable.

Los tres métodos de remuestreo más frecuentemente usados son:

ƒ Vecino mas próximo

ƒ Interpolación bilineal
ƒ Convolución bicúbica

Método del vecino o píxel más cercano

Este método se llama así por situar en cada cuadrado de la malla de la ortofoto el valor de gris del
píxel más cercano en la imagen original. Esta es la solución mas rápida y su principal inconveniente
está en la distorsión que introduce en rasgos lineales de la imagen como carreteras, caminos, ....etc.

Siendo P el punto más cercano, entonces las coordenadas son:

xP = Int ( x + 0.5)
yP = Int ( y + 0.5)

siendo Inte la expresión reducida de Integer, significa un número entero.

190
 TEMA 10: ORTOFOTOGRAMETRÍA DIGITAL

Vecino más proximo

Imagen original Imagen resultado

Figura 6. Remuestreo por medio del vecino más próximo.

Método deInterpolación bilineal

Este método supone promediar los niveles de gris de los cuatro píxeles más próximos de la imagen
original. Este promedio se pondera según la distancia del píxel original al corregido; teniendo mayor
influencia los píxeles más cercanos de la imagen original. Reduce el efecto de distorsión en los rasgos
lineales, pero tiende a difuminar un tanto los contrastes espaciales de la imagen original.

Método de convolución bicúbica

Considera los niveles de gris de los 16 píxeles más próximos. El efecto visual es más correcto que
los anteriores pero supone un volumen de cálculo mucho mayor.

De los tres métodos de remuestreo mencionados, el método del vecino más próximo es el más
simple. Se distingue por su rápida ejecución en cuanto a cálculo se refiere, pero geométricamente
hablando es el menos adecuado, con una precisión de ±0.5 píxel. Los otros dos métodos restantes tienen
una mayor precisión geométrica pero se requiere mayor trabajo en el cálculo. El método de
convolución cúbica es el que absorbe mayor tiempo. En casos generales, el más utilizado es el método
de interpolación bilineal. La elección del método más idóneo dependerá de la finalidad del proceso y de
los recursos disponibles. El más preciso de todos es el método de la convolución que tiene un 0,3 % de
error de interpolación. La interpolación bilineal tiene un valor de error del 3,7 % y el procedimiento del
vecino más próximo del 15,7 % .

10.6. Calidad de la imagen ortofotográfica

Una vez obtenida la imagen ortorrectificada es fundamental la supervisión de la calidad de la

misma, los efectos fundamentales que es necesario localizar son los siguientes:

ƒ Comprobación de que la imagen esté completa. Si el área que cubre la imagen no está cubierta
adecuadamente por un MDE se pueden producir errores en la ortoimagen final.

ƒ Presencia de arrastres en la imagen. Las causas más frecuentes de aparición de arrastres (la imagen
aparece movida) es la presencia de errores en el modelo digital de elevaciones (puntos anómalos)
o bien la presencia de valores de desplazamiento debido al relieve muy elevados al encontrarnos
en el borde de la imagen. Esta última causa es muy frecuente cuando se intenta la rectificación de
una cumbre montañosa, en la que la cara orientada hacia la cámara aparece adecuadamente
rectificada pero en la cara contraria existen arrastres debido a que una pequeña porción de la
imagen original debe “estirarse" para cubrir el área correspondiente en la ortoimagen.

ƒ Aparición de dobles imágenes. Este problema era especialmente importante en los primeros
sistemas fotogramétricos digitales que trabajaban a nivel de modelo en el que errores debidos a la

191
FOTOGRAMETRÍA III

orientación de las imágenes (los parámetros de la orientación externa de la imagen derecha del
primer modelo presenta ligeras diferencias con respecto a los calculados con la misma imagen
cuando hace de imagen izquierda del segundo modelo) o bien falta de continuidad en el modelo
digital de elevaciones.

ƒ Aparición de porciones de ortoimagen no cubiertas. Este hecho tiene el mismo origen que la
aparición de dobles imágenes pero con el resultado contrario, no existe continuidad en la
ortoimagen final por un problema de orientación o del MDE empleado en la ortorrectificación.

ƒ Errores planimétricos. Se comprueban mediante la medida directa de puntos sobre el terreno que
sean perfectamente distinguibles en la imagen, comparando posteriormente sus coordenadas
planimétricas. Los principales problemas que afectan a la calidad de las ortoimágenes generadas
mediante rectificación diferencial son los debidos a los errores o falta de precisión del modelo
digital de elevaciones, dejando a un lado los problemas de ajuste radiométrico que se compensan
mediante correcciones.

10.7. Correcciones radiométricas

Una de las criticas que se hacía a los productos derivados de la ortofotografía analógica era su falta
de homogeneidad (diferencias de contrastes y de tonalidades entre las ortofotos que formaban
mosaicos). En los sistemas digitales, el procesamiento radiométrico mejora la calidad de las imágenes
por procedimientos simples y controlables. La corrección radiométrica de las imágenes tiene dos
objetivos:

ƒ Eliminar las discontinuidades de color entre ortofotos contiguas.

ƒ Crear una representación homogénea de color y tonalidad en un área geográfica grande para la
producción de mosaicos.

En las fotografías aéreas convencionales existen cinco parámetros físicos que contribuyen a estas
variaciones de color y tonalidad:

ƒ La posición del punto de vista.

ƒ La dirección de la iluminación de la escena.
ƒ La distancia focal.
ƒ Las condiciones atmosféricas.
ƒ La reflectancia de los distintos tipos de superficies topográficas.

La presencia de estos parámetros hacen que las correcciones radiométricas sean físicamente
imposibles.

Durante la conversión de la imagen de formato analógico a digital, los datos de entrada del escáner
tienen una gran influencia sobre la calidad de esta y los sistemas mas sofisticados ofrecen un amplio
espectro de corrección selectiva de colores, así como tablas de referencia de color (LUT), por medio de
las cuales se evita la pérdida de información.

El algoritmo de corrección radiométrica se aplica en técnicas de procesamiento de imágenes de bajo

nivel (que no tienen relación con problemas de interpretación), tales como ajuste o realce de contrastes,
transformaciones con las tablas de referencia de color, manipulación de histogramas etc.

La calibración se basa en la hipótesis de Helmhontz sobre la constancia del color (en el espectro
visible e infrarrojo) y permite eliminar las diferencias entre las condiciones de iluminación de las
imágenes. En radiometría (al igual que en el problema de la localización de píxeles) se necesitan
también referencias, uno de los mejores parámetros de calibración es el histograma. Por ejemplo, para
eliminar las diferencia entre los valores de gris de ortofotos contiguas, se parte M hecho que los

192
 TEMA 10: ORTOFOTOGRAMETRÍA DIGITAL

histogramas, en estas zonas de solape de imágenes adyacentes han de ser iguales, lo que implica que se
ha de llevar a cabo un ajuste radiométrico total de ambas imágenes. Las incógnitas serán las tablas de
corrección de los valores de gris de cada imagen

El proceso de corrección radiométrica de datos emplea una gran cantidad de tiempo de ordenador
pero mejora considerablemente la calidad de la imagen; solo después de esta calibración se obtiene la
ortofoto con la cobertura homogénea necesaria para la formación de mosaicos.

10.8. Confección de mosaicos

Cada ortofoto se produce con un cierto grado de recubrimiento con respecto a sus vecinas. Durante
el paso anterior se han eliminado las diferencias en el balance de color pero pueden aparecer otras
perturbaciones como son:

ƒ Diferencias en la orientación de las sombras (en edificios, bosques, etc..).

ƒ Pequeñas discontinuidades geométricas como las que se presentan en puentes y edificios,
relacionadas con efectos de paralaje. Recordemos que los modelos digitales de terreno definen
la superficie topográfica y no tienen en cuenta las altitudes de puentes y edificios.
ƒ Diferencias fenomenológicas en fotografías tomadas en las diferentes estaciones del año.
ƒ Diferencias de reflectancia en zonas de agua (embalses, ríos, lagos). arenosas o de nieve.

En esta fase del proceso de imágenes solo se tienen en cuenta las zonas de recubrimiento entre las
dos ortofotos consideradas, El área común se centra con respecto a unos límites reales comunes y a
continuación comienza un proceso interactivo mediante el cuál se define una trayectoria de corte en la
zona común de tal modo que solo se va a guardar en la base de datos la parte correspondiente a cada
una de las dos ortofotos. Esto tiene las ventajas de:

ƒ Minimizar los efectos visuales de la paralaje.

ƒ Minimizar los efectos visuales de las sombras.
ƒ Realizar una representación homogénea en los casos de diferencias fenomenológicas.

10.9. Estereo-ortofotogrametría digital

La ausencia de la tercera dimensión en las ortofotos puede suplirse con ayuda de las estero-
ortofotos, necesitamos otra imagen que haga de par estereoscópico.

El principio básico puede observarse en la figura 7. Suponemos que tenemos un modelo

tridimensional del terreno que nos sirven para generar la ortofoto y su par estereoscópico. Para la
ortofoto, que es una proyección ortogonal de la superficie terrestre, los puntos de la malla se
transforman mediante la transformación tridimensional de semejanza en puntos en la imagen creando
una malla deformada de cuadriláteros en el sistema x, y de la imagen.

El mismo modelo tridimensional del terreno se somete a una proyección paralela inclinada y los
vértices de la malla deformada en la imagen se determinan de la forma siguiente:

) Se genera una malla cuadrada en el plano XY con un tamaño de cuadrado igual que el del modelo
digital del terreno.

) Generamos un modelo digital con una dirección inclinada con respecto al eje Z pero paralela al
plano XZ.

) Realizamos la proyección de estas coordenadas XYZ en la imagen original.

193
FOTOGRAMETRÍA III

O O

Dirección de

Dirección de
Proyección

Proyección
MDE
MDE

Figura 7. Relaciones entre imagen original, modelo digital, ortofoto y modelo estéreo acoplado. (Kraus, K. Photogrammetry).

Finalmente, se les asigna el valor digital a todos los vértices de la malla deformada en la imagen con
una transformación bilineal.

Una condición necesaria para estos cálculos es el conocimiento de la orientación externa de la foto.
Al ser una proyección con un eje inclinado pero paralelo al plano XZ tendremos las mismas imágenes a
lo largo de la dirección Y, es decir, no tendremos paralajes. En el eje X existirán paralajes relacionadas
con las alturas Z del terreno, figura 8.

Dirección de proyección Dirección de proyección

para la ortofoto para el stereomate

Plano de proyección

Figura 8. Paralajes X función de las alturas terreno. (Kraus, K. Photogrammetry).

Estas dos condiciones, la falta de paralaje en Y y la existencia de paralajes en X son las condiciones
necesarias para una observación estereoscópica, por lo tanto, tenemos un par estereoscópico de una
ortofoto.

194
 TEMA 10: ORTOFOTOGRAMETRÍA DIGITAL

Si las estereo-ortofotos están producidas desde una pasada fotogramétrica con un recubrimiento de
un 60% aproximadamente, entonces la ortofoto se realiza saltando un fotograma y el par estereoscópico
se suele formar con el 50% de la derecha y de la izquierda de una ortofoto; con la foto de la derecha y
con la foto de la izquierda respectivamente, que no han sido tomadas para la ortofoto.

La ventaja de producir una ortofoto y su par con la foto vecina, a parte de poder observar el modelo
estereoscópicamente, es que los detalles tales como tejados, árboles, etc., que no están incluidos en el
modelo digital del terreno, se pueden observar estereoscópicamente.

Las estereo-ortofotos se pueden medir en instrumentos analógicos (óptico-mecánicos), analíticos o

digitales para determinar desniveles entre puntos.

10.10. Propiedades de las ortofotos digitales

ƒ Las ortofotos digitales están geocodificadas y ofrecen al menos la misma precisión que los mapas
de línea. Cada píxel tiene una posición claramente determinada en el sistema de referencia
empleado (UTM, Lambert, etc.)
ƒ Están calibradas respecto a los colores y por medio de los módulos de ajuste radiométricos y de
formación de mosaicos se pueden fundir dos o mas ortofotos y producir una imagen de color
homogénea de mejor calidad que la fotografía aérea original.
ƒ Las ortofotos digitales son multitemporales. Los vuelos fotográficos con frecuencia se
interrumpen por causas atmosféricas durante largos periodos, al cabo de los cuales pueden haber
cambiado las condiciones de iluminación y los aspectos de los cultivos. Con esta tecnología
pueden obtener documentos homogéneos sin las discontinuidades de color resultantes en los
métodos clásicos. Gracias a esta cualidad, las ortofotos digitales pueden utilizarse para poner al
día de una forma rápida los mapas existentes.
ƒ Las ortofotos digitales pueden adoptar escalas variables, así como es posible fundir dos ortofotos
digitales, también es posible.
ƒ Las ortofotos digitales pueden adoptar escalas variables; a partir de una escala de consistencia se
pueden formar ortofotos a distintas escalas (menores)
ƒ Las ortofotos digitales son multiespectrales. Al ser imágenes de tipo raster se pueden aplicar
técnicas de clasificación a las fotos ya sean en color o en infrarrojo color. Las propiedades de
multitemporal, multiescala y multiespectral son exclusivas de las ortofotos digitales, y no existen
en los documentos analógicos.
ƒ Debido a sus propiedades multiespectrales y multiescalares pueden compararse con datos
provenientes de teledetección con objeto de completar el proceso de interpretación con otros
canales multiespectrales.
ƒ Permiten la sobreimpresión de mapas de línea o temáticos para la formación de documentos más
claros y precisos.
ƒ Es obvio decir que se pueden utilizar para hacer cálculos de superficies, distancias,
desplazamientos...etc.
ƒ Se pueden adaptar a la superficie de un DTM para formar perspectivas, a menudo utilizadas en
estudios de impactos ambientales.
ƒ Debido a su carácter multidisciplinario, constituyen una fuente de datos muy útil para usuarios de
GIS; por esta razón las ortofotos constituyen un complemento esencial de cualquier tipo de
aplicación GIS al ofrecer información temática específica en formato vectorial y raster.

Bibliografía

ƒ Chueca P., M., Topografía, Dossat, Madrid (1982), Tomo II, 84-237-0591-9.
ƒ Curso sobre técnicas fotogramétricas y sus aplicaciones, CEDEX (MOPTMA), Madrid (1994).
ƒ Delgado, J., Apuntes de Fotogrametría Digital, Universidad Politécnica de Jaén, 1997.
ƒ Gómez, A., Apuntes de Fotogrametría, Universidad Politécnica de Madrid, ETSIA.

195
FOTOGRAMETRÍA III

ƒ Gómez, A., Flujo de trabajo en Fotogrametría Digital, Mapping, nº 53, Abril de 1999.
ƒ Heras, A.E., Apuntes de Fotogrametría, Universidad Politécnica de Madrid, EUITT.
ƒ Krauss, K., Photogrammetry, Dümmler/Bonn (1993), 3-427-78684-6.
ƒ Moffitt, F.H. y Mikhail, E.M., Photogrammetry, Harper & Row Publishers., New York (1980), 3ª
edición.
ƒ Sanchez Espeso, J.M., Análisis del proceso de rectficación de una imagen aérea de eje vertical
para obtener una imagen digital, Mapping, nº 64, septiembre de 2000.
ƒ Wolf, P.R.; Dewitt, B. A, Elements of Photogrammetry with Applications in GIS, McGraw-Hill,
EEUU, Tercera Edición (2000). 0-07-292454-3.
ƒ Zhizhuo, W., Principles of Photogrammetry, (1990), 7-81030-00-8.

196