Correccion Geometrica PDF

Correccin geomtrica
z Proceso por el cual la geometra de la

imagen se hace planimtrica.
z Generar funciones que transformen las

coordenadas de imagen (lnea y columna)
a coordenadas de mapa (UTM u otras).
z La referencia puede ser un mapa o una

imagen de otra fecha.
Tipos de correccin geomtrica

z Orbitales:
z Erroressistemticos.
z Conocidas a priori.
z Automatizables.
z Estadsticas (por PUNTOS DE CONTROL)

z Modelizacin inductiva.
z Errores aleatorios.
z Necesitan cartografa de apoyo / GPS.
Chuvieco, 2002
1
Correccin por puntos de control:
fases
z Localizacin de los puntos de control
z Interpolacin espacial: clculo de las
funciones que transformarn las
coordenadas de la imagen en las del
sistema de referencia.
z Interpolacin radiomtrica: transferencia
de los valores digitales a su nueva
posicin.
Localizacin de puntos de control

z Muestra de las distorsiones geomtricas que presenta la
imagen. La calidad de la correccin depende de una
correcta seleccin de los PC.
z Localizacin:
z Identificable en la imagen y el documento de referencia (mapa o
imagen o sobre el terreno con GPS).
z Lugares estables.
z Distribucin:
z Que cubran toda la zona.
z A distintas alturas si las hay.
z Nmero (cuantos ms, mejor):
z Depende del tamao de la imagen y el grado de distorsin
geomtrica.
z El mnimo necesario viene dado por el grado de la ecuacin de
ajuste (3 para grado 1, 6 para grado 2, 10 para grado 3, etc).
etc).
2
Ejemplo de puntos de control
Chuvieco, 2002
Clculo de las funciones de

transformacin (Interpolacin espacial)
z Funciones matemticas que ponen en
relacin las coordenadas de la imagen
con las del sistema cartogrfico de
referencia.
z Se calculan por regresin mltiple a partir

de las coordenadas de los puntos de
control.
3
z Ejemplo ms simple: transformacin de primer
grado:
X = a + bC + dL
Y = e + fC + gL
A, b, d, e, f y g son los coeficientes de regresin, C y L son las
las variables
independientes (columnas y lneas de la imagen), X e Y son las
variables dependientes)
z Slo produce rotacin o traslacin de la imagen.

Corrige slo deformaciones lineales.

z Un poco ms complejo: transformacin de

segundo grado:
X = a + bC + cL + dC2 + eL2 + fCL
Y = g + hC + iL + jC2 + kL2 + lCL
z Los coeficientes se derivan tambin por
regresin mltiple, pero ahora hay ms
variables independientes (C, C2, L, L2 y CL)
z Corrige deformaciones ms complejas.
4
Valoracin del error
z El error cometido por la correccin se evala
comparando las coordenadas reales de los PC
con las que ha calculado la ecuacin de
transformacin, a travs del
Error medio cuadrtico (RMS) o residual para

cada punto y el error medio para toda la
imagen (media de los errores para cada PC):
Error medio cuadrtico o residual

para cada PC
RMS = (ccorr - corig )2 + ( l corr - l orig )2
Corregido
Lcorr
RMS
Lorig
Original
Lneas
Corig Ccorr Chuvieco, 2002

Columnas
5
Residuales
Col. Col. est. Resid Fila Fila Est. Resid Error
992,00 990,96 1,04 102,00 102,60 -0,60 1,20
936,00 936,65 -0,65 195,00 194,05 0,95 1,15
923,00 923,34 -0,34 431,00 431,33 -0,33 0,47
810,00 810,02 -0,02 498,00 497,90 0,10 0,10
651,00 650,79 0,21 453,00 453,16 -0,16 0,27
653,00 653,33 -0,33 264,00 263,31 0,69 0,76
677,00 676,56 0,44 117,00 117,85 -0,85 0,95
464,00 465,04 -1,04 230,00 229,50 0,50 1,16
253,00 252,31 0,69 512,00 512,30 -0,30 0,75
Col. = 15.177,79 + 0,032452 X - 0,006443 Y

Lin = 147.365,5 - 0,006590 X - 0,032202 Y
Chuvieco, 2002
Residuales
z Se elimina el punto con mayor RMS, y se
calculan de nuevo las ecuaciones de ajuste
y los RMS para cada PC y el RMS medio
para toda la imagen. Si el error ya es
aceptable, se procede con la correccin y si
no lo es, se repite el proceso hasta que lo
sea.
6
Interpolacin radiomtrica
z Se ha creado una imagen vaca, con la
geometra correcta. Ahora hay que
transferir a cada nueva posicin el valor
digital que le corresponda: interpolacin
radiomtrica o resampling.
z Se puede hacer por distintos mtodos:
z Vecino ms prximo
z Interpolacin bilineal
z Interpolacin cbica
Criterios de interpolacin
IMAGEN ORIGINAL IMAGEN CORREGIDA

Chuvieco, 2002
VECINO MAS PRXIMO INTERPOLACIN BILINEAL CONVOLUCIN CBICA
7
Vecino ms prximo
+ No modifica los valores.

+ Es el ms rpido
- Aparecen elementos lineales en escalera
- Pueden perderse valores
Chuvieco, 2002
Interpolacin bilineal
Convolucin cbica
A cada celdilla de la imagen corregida se le asigna como
ND la media de los ND de las cuatro (interpolacin bilineal)
o diecisis (convolucin cbica) celdillas de la imagen
original que estn ms cerca. Es una media ponderada en
funcin de la distancia del pixel original al corregido.
Estos mtodos:
- Interpolan valores.
- Son ms lentos.
+ Se suavizan las transiciones lineales.
+ Mejor aspecto visual.
8
Efecto del tipo
de
interpolacin
a) original;
b) vecino ms prximo;
c) Interpolacin bilineal;
d) convolucin cbica
Chuvieco, 2002
En Arc View, la correccin geomtrica de

imgenes se hace con la extensin
ImageWarp
9
EJERCICIO: Imagewarp.
Vamos a corregir la geometra de una imagen de satlite procedente del sensor

Landsat TM para que se ajuste a nuestra base cartogrfica (UTM huso 30, ED-50). El
sistema de adquisicin de una imagen de satlite implica la existencia de una serie de
errores geomtricos no sistemticos que solamente pueden corregirse de una forma
emprica mediante una transformacin rubber sheet o polinmica, a travs de puntos
de control. Esta operacin la realizan todos los programas de tratamiento digital de
imgenes y SIGs raster. En este caso, Arc View la realiza a travs de una extensin:
ImagenWarp, muy similar a Shape Warp que acabamos de ver.
Para hacer este ejercicio podis elegir entre la imagen Landsat TM de la hoja de
Almera (hoja 1045, archivo tmalme.bip) y la de la hoja de Ubrique (hoja 1050, archivo
tmubri.bip).
1. Abre Arc View y crea una vista nueva. Aade la imagen sin corregir de Ubrique o
Almera (por ejemplo, tmalme.bip). Mira sus coordenadas: son filas y columnas, no
est georreferenciada y su geometra no se ajusta a nuestra base cartogrfica.
2. Aade como fondo la imagen SPOT de la misma hoja (sp1045.bil para almera, por
ejemplo). La geometra de esta imagen se corresponde con la proyeccin UTM del
huso 30, datum ED-50. Como ves, esta imagen s tiene coordenadas, y por supuesto
no se superpone con la anterior.
3. Cierra Arc View.
4. Instala la extensin imagewarp, descomprimiendo el archivo imagewarp.zip y

ejecutando el archivo setup.exe. As instalars la extensin ImageWarp en el directorio
ext32.
5. Abre Arc View y activa la extensin ImageWarp (ProjectExtensions).
6. Para iniciar la sesin, selecciona la opcin ImageWarp Session del men

ImageWarp.
7. Seleccionaremos la imagen que vamos a corregir y la que vamos a utilizar como
referencia (la referencia no tienen porqu ser una imagen, puede ser un shape o
cualquier otro tipo de tema que Arc View pueda mostrar)
8. A la pregunta de si queremos indicar la proyeccin para la vista TO (de referencia),

responderemos no (porque nuestra referencia es una imagen tambin y ya sabemos
que Arc View no proyecta imgenes al dibujarlas en una vista).
9. A continuacin nos pregunta el nombre de la tabla que almacenar las coordenadas

de los puntos de control en la imagen a corregir y en el sistema de referencia.
Podemos indicar si la tabla ya existe (porque hayamos dejado la sesin a medias
anteriormente, por ejemplo).
10. Ahora aparecer la pantalla principal, que tendr este aspecto, con las
herramientas en la ventana central y tres vistas, que corresponden al tema a corregir
(FROM) y al tema de referencia (TO, en la que se puede hacer zoom y TO**ROAM,
que sirve como referencia):
11. Ahora haremos zoom en la vista TO (podemos hacerlo en la TO**ROAM y el
efecto tendr lugar en la vista TO) y en la FROM para localizar algn punto en comn
(punto de control), y despus lo seleccionaremos primero en TO y luego en FROM,
con la herramienta .
12. Repetiremos esta operacin tantas veces como deseemos, sabiendo que cuantos
ms puntos de control encontremos mejor ser la transformacin y que el mnimo
matemtico para realizar una transformacin de primer grado es de 3 puntos, 6 para
una de segundo grado y 10 para una de tercer grado.
13. Si nos equivocamos, podemos borrar un punto de control seleccionndolo con la

herramienta y utilizando despus el botn .
14. Una vez que hayamos encontrado 3 puntos de control, podemos usar el botn
para que el sistema, mediante una ecuacin polinomial provisional, site el punto
de control en la vista FROM una vez que nosotros lo hayamos localizado en la vista
TO. Entonces nosotros slo tendremos que desplazarlo para que se site en el lugar
exacto. Para ello, lo seleccionaremos y luego apuntaremos a su posicin correcta con
la herramienta .
15. Cuando tengamos un nmero de puntos de control suficiente, podemos calcular el

error medio (RMS), en unidades del sistema de referencia (en nuestro caso metros,
que son las unidades en que estn expresadas la coordenadas en la referencia que es
la imagen SPOT) que resultara si hiciramos la transformacin con esos puntos de
control con una ecuacin de orden 1 a 12. Esto lo haremos con el botn .
Normalmente, al corregir una imagen respecto a otra, en zonas donde no hay mucho
relieve, se consiguen buenos resultados con una ecuacin de orden 1. Para grados
superiores hay que asegurarse de que el nmero de puntos de control es elevado y su
distribucin muy regular por toda la imagen.
16. Tambin podemos examinar el error (en unidades de referencia, metros en este
caso) que se producira en la posicin de cada punto de control si realizramos la
transformacin en este momento. Para ello se usa el botn , que muestra la tabla
de los puntos de control, con el error calculado para cada uno de ellos en la ltima
transformacin que hayamos seleccionado.
17. Si algn punto tiene un error demasiado grande, podemos desactivarlo para que
no sea usado en la transformacin y ver si as mejora el RMS general, con el botn
.
18 Cuando finalmente estemos satisfechos con el error medio cuadrtico de la
transformacin, pulsaremos el botn para crear la nueva imagen, con la

geometra transformada. Se aplicar la ltima ecuacin que se haya seleccionado.
Hay que indicar el tamao del pixel en la imagen corregida (debera ser el mismo que
en la imagen original, en este caso 30 metros) y tambin el nombre y ubicacin para la
nueva imagen. Podemos llamarla cgalme.bip. En este formato, solamente corrige las
bandas de la imagen que se estn usando en la vista FROM.
19. Cierra el proyecto, y abre una nueva sesin de Arc View. Aade a la vista la
imagen SPOT (1045.bil), la imagen corregida (cgalme) y algn mapa de la zona (un
10.000 50.000). Has conseguido un buen ajuste?

Correccion Geometrica PDF

Cargado por

Copyright:

Formatos disponibles

Correccion Geometrica PDF

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Correccion Geometrica PDF

Cargado por

Copyright:

Formatos disponibles

Correccin geomtrica

z Proceso por el cual la geometra de la

z Generar funciones que transformen las

z La referencia puede ser un mapa o una

Tipos de correccin geomtrica

z Estadsticas (por PUNTOS DE CONTROL)

Localizacin de puntos de control

Clculo de las funciones de

z Se calculan por regresin mltiple a partir

z Slo produce rotacin o traslacin de la imagen.

Clculo de las funciones de

z Un poco ms complejo: transformacin de

Error medio cuadrtico (RMS) o residual para

Error medio cuadrtico o residual

RMS = (ccorr - corig )2 + ( l corr - l orig )2

Corig Ccorr Chuvieco, 2002

Col. = 15.177,79 + 0,032452 X - 0,006443 Y

IMAGEN ORIGINAL IMAGEN CORREGIDA

VECINO MAS PRXIMO INTERPOLACIN BILINEAL CONVOLUCIN CBICA

+ No modifica los valores.

En Arc View, la correccin geomtrica de

Vamos a corregir la geometra de una imagen de satlite procedente del sensor

3. Cierra Arc View.

4. Instala la extensin imagewarp, descomprimiendo el archivo imagewarp.zip y

5. Abre Arc View y activa la extensin ImageWarp (ProjectExtensions).

6. Para iniciar la sesin, selecciona la opcin ImageWarp Session del men

8. A la pregunta de si queremos indicar la proyeccin para la vista TO (de referencia),

9. A continuacin nos pregunta el nombre de la tabla que almacenar las coordenadas

13. Si nos equivocamos, podemos borrar un punto de control seleccionndolo con la

15. Cuando tengamos un nmero de puntos de control suficiente, podemos calcular el

18 Cuando finalmente estemos satisfechos con el error medio cuadrtico de la

transformacin, pulsaremos el botn para crear la nueva imagen, con la

También podría gustarte