Estudio Piloto para la Generaci�n de Cartograf�a por Sensores Remotos
"Venezuelan Cartographic Sensor Technology Study"
Lic. Ramiro Salcedo G., MSc.
Centro de Procesamiento Digital de Im�genes/FII
Email: ramiros@fii.org
RESUMEN
Se describe un estudio piloto, el cual se realiz� con financiamiento de la Agencia de Comercio y Desarrollo de los Estados Unidos (TDA), para CVG-Tecm�n.
El prop�sito del estudio fu�, en primer lugar, proveer un an�lisis de las diversas tecnolog�as de sensores remotos con el fin de evaluar las fortalezas y debilidades de cada uno de estos para la generaci�n de cartograf�a b�sica digital al sur de Venezuela; y en segundo lugar producir mapas topogr�ficos a escala 1:50.000 con curvas de nivel de 40 metros de intervalo.
El �rea de estudio seleccionada, tiene una superficie aproximada de 500 km2 (carta 7929-I) y se localiza en los alrededores de la poblaci�n de Santa Elena de Uair�n, al sureste del Estado Bol�var, en la frontera con Brasil (ver Fig. 1).
Fig 1: Localizaci�n Area Estudio
Introducci�n: La producci�n de cartograf�a topogr�fica b�sica en medios tropicales enfrenta importantes limitaciones naturales y log�sticas. La cobertura de nubes se cuenta entre las limitaciones naturales determinantes por cuanto impide la toma de datos, ya sea desde plataformas aerotransportadas u orbitales. En este sentido, se plantea de manera perentoria, la utilizaci�n de sensores y t�cnicas distintas a las convencionales para la generaci�n de mapas en aquellas regiones donde, como es el caso del sur de Venezuela, no se ha podido completar el cubrimiento cartogr�fico por medios convencionales.
Objetivo: Utilizar diversos sensores remotos �pticos y de radar; tanto aerotransportados como orbitales, con el objeto de evaluar sus posibilidadess de aplicaci�n con fines de generaci�n de cartograf�a topogr�fica a escala 1:50.000 con curvas de nivel cada 40 metros. El �rea de estudio escogida es Santa Elena de Uair�n en el Estado Bol�var.
Para cumplir dicho cometido, se constituy� un equipo de trabajo formado por las siguientes instituciones y empresas:
Miembros del equipo de trabajo
Instituci�n/Empresa |
Rol |
Raytheon Systems Company, Inc. | Gerencia del proyecto, procesamiento de las fotograf�as a�reas y evaluaci�n de los datos. |
Instituto de Ingenier�a | Recolectar los PCT y generaci�n de los mapas |
Vexcel Corporation | Procesamiento de las im�genes SAR |
Organizaciones de Apoyo Adicional
Empresas de apoyo |
Rol |
Ashtech | Apoyo al levantamiento GPS |
OMNI Ltd. | Adquisici�n de Fotograf�as a�reas digitales |
The ABCs of XYZ | Procesamiento Diferencial de datos GPS |
SIGIS | Digitalizaci�n de las fotos a�reas |
Aerograf�as de Venezuela | Adquisici�n de misi�n aerofotogr�fica. |
El equipo de trabajo realiz� las tareas relativas a la selecci�n de los sensores, recolecci�n y procesamiento de datos, generaci�n de los mapas y an�lisis final de los resultados.
ETAPAS DEL ESTUDIO
1.1 Selecci�n de los sensores
Basados en un an�lisis comparativo de los sistemas disponibles, se decidi� la adquisici�n de un grupo de sensores diferentes. La lista de los mismos se limit� a aquellos comercialmente disponibles que ya estuviesen adquiridos y a sensores que podr�an ser programados para adquisiciones nuevas. La lista inicial incluy� tanto sistemas aerotransportados como orbitales disponibles. Se analizaron y compararon las especificaciones t�cnicas del funcionamiento de estos sensores para determinar cuales podr�an cumplir con los requerimientos cartogr�ficos de Venezuela; ya fuese individualmente o en una combinaci�n de dos o m�s de estos.
1.2 Adquisici�n de datos
La recolecci�n de datos consisti� en la obtenci�n de las im�genes y los puntos de control terrestre con GPS en el �rea de estudio de Santa Elena.
Se contact� a los proveedores de datos que cumplieron con los criterios de selecci�n, a fin de iniciar la adquisici�n. Se investig� la disponibilidad y utilidad de datos existentes en los archivos. Se programaron nuevas adquisiciones de aquellos sensores de los que no hubiese datos adquiridos. Se ordenaron las im�genes que tuviesen coberturas adecuadas. La cobertura de nubes predominante en la regi�n de Santa Elena jug� un papel importante en la adquisici�n de nuevos datos.
Se requiri� la adquisici�n de Puntos de Control Terrestre (PCT) sobre el �rea piloto, para apoyar la aerotriangulaci�n y la evaluaci�n de la precisi�n de cada tipo de sensor. Se realiz� el trabajo de campo necesario para colectar dichos PCT. Esta informaci�n consisti� de coordenadas X,Y,Z grabadas por un posicionador GPS sobre localidades espec�ficas del �rea de estudio. Los datos de los PCT fueron usados para geo-referenciar en forma precisa las im�genes con respecto a puntos conocidos en el terreno, as� como proveer un est�ndar comparativo para medir la precisi�n vertical y horizontal de los mapas.
1.3 Procesamiento de datos
Se planific� el procesamiento de los datos con el objetivo de elaborar mapas de la zona de estudio, a partir de cada uno de los sensores disponibles. Al irse recibiendo las im�genes, los correspondientes grupos de trabajo se dedicaron a generar Modelos de Elevaci�n Digital (DEM) e im�genes orto-rectificadas por cada tipo de sensor. En cada caso, el procesamiento de los datos se hizo de acuerdo a los requerimientos de precisi�n de la imagen. Los DEM y las orto-im�genes sirvieron posteriormente para la extracci�n de elementos planim�tricos y altim�tricos.
1.4 Generaci�n de mapas
La extracci�n de elementos planim�tricos de las orto-im�genes permiti� generar las capas vectoriales relativas a cobertura vegetal, �reas urbanas y tipos de vialidad. Los DEM a su vez, fueron procesados para producir las curvas de nivel con un intervalo de 40 metros. Todos los detalles relativos a leyendas, t�tulos, toponimia y marginalia, fueron seleccionados de acuerdo a est�ndares cartogr�ficos oficiales. Todos los DEM y ortoim�genes fueron enviados al Instituto de Ingenier�a (CPDI) para ser convertidos en mapas, de acuerdo a los requerimientos de Cartograf�a Nacional.
Sensores y datos utilizados:
Se evaluaron todos los sensores comerciales posibles, con el fin de determinar los candidatos m�s factibles a ser utilizados en el marco de este proyecto. As� tenemos: �stos y sus caracter�sticas:
Sensores orbitales y aerotransportados
Sensor |
Descripci�n |
LANDSAT |
LANDSAT es operado por el gobierno de los Estados Unidos y distribu�do por Spaceimaging. Las escenas adquiridas por Landsat son multiespectrales y monosc�picas, y por lo tanto no pueden ser usadas para extracci�n de datos de elevaci�n. |
SPOT |
SPOT (Satellite Pour l'Observation de la Terre) es operado por SPOT Image en Toulouse, Francia. SPOT puede adquirir im�genes �pticas en modos monosc�pico y estereosc�pico. Las im�genes est�reo pueden ser usadas para producir Modelos de Elevaci�n Digital (DEM) e im�genes orto-rectificadas. |
IRS-1C |
El Instituto Indio de Investigaci�n Espacial y la Agencia Nacional de Sensores Remotos operan el sat�lite Indio (IRS). Las im�genes son distribu�das por Spaceimaging. El IRS-1C,D colecta im�genes �pticas y monosc�picas solamente. |
KVR-1000 |
KVR-1000 es parte del sistema satelital Ruso de c�maras �pticas SPIN-2 que opera con pel�cula. Provee im�genes monosc�picas y est�reo. Despu�s de su exposici�n a bordo del sat�lite, el cartucho de pel�cula es enviado a tierra por paraca�das, donde es recuperado. |
Fotograf�as A�reas |
Las fotograf�as a�reas se adquieren desde aviones a diversas alturas y con c�maras m�tricas especiales, usando un lente de 6" y formato de pel�cula de 9"x 9". Los negativos fotogr�ficos deben ser convertidos a formato digital usando scanners de alta resoluci�n, para luego realizar el procesamiento automatizado. La fotograf�a puede ser adquirida en un amplio rango de escalas y resoluciones, seleccionando las alturas de vuelo apropiadas. |
Sensores orbitales y aerotransportados (continuaci�n)
Sensor |
Descripci�n |
Fotograf�a A�rea Digital |
La fotograf�a a�rea Digital es adquirida desde aviones a distintas alturas usando sistemas electro-�pticos para capturar la imagen digitalmente, y as� evitar el proceso de escaneo. Al igual que con la fotograf�a a�rea convencional, �stas pueden ser adquiridas en un amplio rango de escalas y resoluciones, seleccionando las alturas de vuelo apropiadas. |
ERS-1, 2 |
ERS-1, 2 (European Earth Observation Radar Satellite) Es un par de sat�lites SAR en banda C operados por European Space Agency. El sistema puede ser usado para producir DEM v�a pares tandem de interferometr�a. |
JERS-1 |
Japanese Earth Resources Satellite (JERS) es un sat�lite SAR en banda L y �ptico multiesprectral que ha sido usado extensivamente para adquirir im�genes del Amazonas Brasile�o y de los bosques tropicales alrededor del mundo. Se pueden generar DEM y ortoim�genes a partir de interferometr�a de �rbitas repetidas. |
RADARSAT |
RADARSAT es un sat�lite SAR en banda C lanzado en 1995 y operado por RADARSAT International (RSI). El sistema puede adquirir im�genes est�reo para producci�n de DEM y ortoim�genes en distintas resoluciones espaciales y �ngulos de incidencia. |
IFSAR |
Interferometric Synthetic Aperture Radar (IFSAR) es un sistema aerotransportado desarrollado bajo contrato con US DARPA, actualmente operado por Intermap Corp. Usa un sistema SAR dual en banda X para producir DEM de alta resoluci�n mediante interferometr�a SAR. |
Sensores Hiper-espectrales |
Los sensores hiper-espectrales (Satelitales y aerotransportados) tales como AVIRIS (250 bandas), proveen datos de reflectancia sobre un amplio rango de bandas espectrales cont�guas y de ancho uniforme. Teniendo canales de un ancho tan peque�o como algunos nanometros, los espectros resultantes pueden ser usados para caracterizar vegetaci�n y geolog�a con m�s alta precisi�n que utilizando los sensores multi-espectrales existentes. |
Sensores futuros |
El gobierno de los EEUU ha licenciado recientemente, cinco empresas comerciales para lanzar sat�lites de observaci�n de la tierra, como el IKONOS-1 de Space Imaging Inc. Pancrom�tico y con resoluci�n espacial de 1 metro. |
Producci�n cartogr�fica
Dos mapas topogr�ficos con sus elementos altim�tricos y planim�tricos (curvas de nivel, carreteras, r�os, centros poblados, y cobertura natural), fueron producidos a partir de los datos de los sat�lites SPOT y Radarsat.
Mapa a partir de SPOT: Una escena Spot-Pancrom�tica orto-rectificada fue usada para hacer la interpretaci�n de aquellos elementos planim�tricos requeridos en un mapa topogr�fico standard a escala 1:50.000, tales como red hidrogr�fica, cuerpos de agua, cobertura natural, vialidad y areas urbanas. (ver figura N� 2).
Fig. 2: Imagen SPOT-P con planimetr�a
Este procedimiento involucr� el uso del software PCI para visualizar la imagen y hacer la captura en pantalla de dichos elementos. Se utiliz� la leyenda establecida por Cartograf�a Nacional-MARNR como el modelo a seguir.
Durante el proceso de interpretaci�n o captura de los elementos planim�tricos antes se�alados, se asign� a cada uno de �stos un c�digo, tal que cuando los datos fueran transferidos hacia el ambiente de edici�n cartogr�fica digital, �stos ya tuviesen su identificaci�n respectiva y as� hacer el proceso de edici�n m�s eficiente. Se demostr� que una de las ventajas de interpretar en pantalla directamente sobre la imagen, est� en el hecho de que permite mejorar los contrastes de esta donde y cuando se requiere, as� como hacer los acercamientos o ampliaciones necesarias para obtener el mejor detalle posible.
Adicionalmente, utilizando un Modelo de Elevaci�n Digital (DEM) generado con dos im�genes estereosc�picas del sat�lite SPOT-P, se produjeron las curvas de nivel con intervalos de curvas cada 40 metros (ver figura N� 3) . Es evidente la presencia de "huecos" los cuales se deben a la presencia de nubes en las im�genes.
Fig. 3: DEM generado a partir de im�genes est�reo SPOT-P
La integraci�n de las curvas de nivel (altimetr�a) y de la planimetr�a (r�os, v�as, ciudades, vegetaci�n) permiti�, a trav�s de un proceso de edici�n cartogr�fica digital, hacer la producc��n del mapa final (ver figura 4).
Fig. 4: Mapa con planimetr�a a partir de im�genes SPOT-P
Mapa a partir de Radarsat: Se utiliz� un par de im�genes adquiridas en modo standard-3 (ver figura 5).
Fig. 5: Imagen Radarsat del �rea de estudio
Los datos de este sat�lite son sustancialmente distintos a aquellos de SPOT. Por tratarse de un sensor de microondas, las im�genes no se ven afectadas por las condiciones clim�ticas o de iluminaci�n. En este sentido, las im�genes de radar captan datos tanto en horas diurnas (orbita descendente), como nocturnas (�rbita ascendente) y tienen la capacidad de penetrar la cobertura de nubes. Esta caracter�stica las hace particularmente �tiles en medios como los del �rea de estudio, donde la presencia cuasi-permanente de �stas impide la toma de datos por sensores �pticos durante la mayor parte del a�o.
La adquisici�n de dos im�genes estereosc�picas del sat�lite Radarsat permiti� generar un DEM (ver figura 6),
Fig. 6: DEM de Radarsat
del cual se produjo a su vez, el archivo de curvas de nivel con intervalo de 40 metros, usando el m�dulo TIN (triangular Irregular Network) de Arc/Info (ver figura 7). Aqu� se observa la ventaja derivada de la utilizaci�n de datos de radar, los cuales no se ven afectados por la interferencia nubosa.
Fig. 7: curvas de nivel a partir del DEM Radarsat
La planimetr�a para este mapa no se extrajo de la imagen de radar por cuanto este sensor no permite una �ptima identificaci�n de este tipo de elementos. En consecuencia, se tom� la producida a partir de SPOT-P y se integr� a las curvas de nivel derivadas del DEM de Radarsat. El mapa as� producido es el resultado de la utilizaci�n de lo mejor de las dos plataformas (ver figura 8).
Fig. 8: Mapa resultante de Radarsat y SPOT
1.5 An�lisis de resultados
Se realiz� un an�lisis para evaluar las distintas fuentes de datos, en t�rminos de calidad de las im�genes y su interpretabilidad; errores verticales y horizontales de los DEM y ortoim�genes.
Resumen del procesamiento de los datos
Sensor |
Resultados |
Acci�n |
Fotograf�a a�rea |
Inadecuados procedimientos de adquisici�n de la misi�n aerofotogr�fica y del escaneo impidieron el procesamiento de toda el �rea piloto. Se gener� el DEM y una ortoimagen en un �rea limitada. | No se gener� mapa. |
SPOT |
Se produjo el DEM, las curvas de nivel y la ortoimagen. La cobertura de nubes caus� "huecos" en el DEM. Se extrajo toda la planimetr�a. | Generar y analizar mapa usando las curvas y la planimetr�a. |
RADARSAT |
Se produjo el DEM, las curvas de nivel y la orto-imagen. No fue posible extraer las caracter�sticas planim�tricas. | Generar y analizar el mapa usando la planimetr�a de SPOT y las curvas de nivel de RADARSAT. |
ERS-1,2 |
Baja correlaci�n sobre �reas de bosque caus� amplios "huecos" en el DEM. No se produjeron ortoim�genes ni curvas. | No se gener� mapa. |
JERS-1 |
L�nea de base corta caus� una sensibilidad inmanejable. No se produjeron curvas ni ortoim�genes. | No se gener� mapa. |
Evaluaci�n de errores Horizontales y Verticales
Como se recomend� en la secci�n de procesamiento de datos, se realiz� una evaluaci�n de errores en los datos de SPOT, Radarsat y fotos a�reas.
Horizontales:
Estos errores son la diferencia en localizaci�n de un elemento observado en una imagen (puntos de chequeo) y su localizaci�n real en X y Y (latitude y longitud) a partir de un origen controlado. En este estudio, los PCT procesados, fueron usados como la fuente de control y los elementos en las otras im�genes fueron usados como puntos de chequeo.
Verticales:
Los errores verticales son la diferencia en la elevaci�n (Z) de una localizaci�n en un DEM y la elevaci�n observada por una fuente de control. En este estudio, Los PCT procesados fueron usados como control, y los DEM derivados de Radarsat, SPOT y fotos a�reas, usados como puntos de chequeo. Dos m�todos fueron usados para localizar el punto de chequeo en el DEM. El primer y m�s com�n m�todo fue muestrear la elevaci�n en el DEM en la localizaci�n X,Y del PCT. El segundo m�todo tom� en cuenta el l�mite observado en el an�lisis de errores horizontales. En este m�todo, la elevaci�n fu� muestreada en la localizaci�n observada del elemento, en lugar de la localizaci�n absoluta de latitud, longitud. Como se oberva en la tablas inferiores, el �ltimo m�todo di� lecturas de elevaci�n m�s precisas.
Res�men de evaluaci�n de errores horizontales
RADARSAT |
SPOT |
Fotos a�reasl |
||||
Base (Metros) |
Desplazam. (Metros) |
Base (Metros) |
Desplazam (Metros) |
Base (Metros) |
Desplazam |
|
D E promedio |
34.0 |
0.0 |
(62.2) |
0.0 |
(4.3) |
- |
D N promedio |
22.8 |
0.0 |
(25.3) |
0.0 |
4.2 |
- |
CE promedio |
145.4 |
143.3 |
72.7 |
30.9 |
12.5 |
- |
CE Minimo | 49.7 |
81.5 |
18.9 |
5.4 |
6.7 |
- |
CE Maximo | 236.5 |
236.2 |
96.2 |
69.9 |
24.9 |
- |
CE Desv. Est. |
81.9 |
72.7 |
19.9 |
14.3 |
6.2 |
- |
CE RMSE | 161.8 |
156.5 |
75.2 |
33.8 |
13.8 |
- |
Res�men de evaluaci�n de errores en elevaci�n
RADARSAT |
SPOT |
Fotos a�reas |
||||
En PCT (Metros) |
Elemento (Metros) |
En PCT (Metros) |
Elemento (Metros) |
En PCT (Metros) |
Elemento (Metros) |
|
Error promedio | 48.2 |
28.6 |
33.2 |
34.9 |
2.6 |
2.7 |
Error m�nimo | 1.4 |
0.4 |
3.6 |
4.6 |
0.3 |
0.3 |
Error m�ximo | 136.7 |
51.4 |
88.8 |
104.8 |
7.8 |
7.5 |
Desv. Est�ndar | 34.4 |
21.7 |
23.2 |
26.5 |
2.2 |
2.2 |
RMSE | 58.9 |
34.6 |
40.2 |
43.4 |
3.3 |
3.4 |
CONCLUSIONES: