Creación de un mapa-mosaico en formato ráster ECW a partir de documentos PDF

Se explica en este artículo el procedimiento para obtener a partir archivos PDF y del mosaico de imágenes obtenidas de ellos, un único archivo de mapa georreferenciado en formato ráster ECW.

Para ello, he tomado como ejemplo el Mapa Geológico de Catalunya 1:50.000 del IGC (Institut Geològic de Catalunya) publicado por comarcas y en formato PDF de baja resolución.

Los documentos PDF del Mapa geològic comarcal de Catalunya a escala 1:50.000 se pueden descargar de la web del IGC:

http://www.igc.cat

Una vez descargados, comprobaremos que se encuentran "bloqueados" para su edición. La primera acción es la de desbloquear dichos PDF para poder extraer de ellos las imágenes (leyendas, mapas en sí, etc.). En la web se encuentra un amplio abanico de sencillas utilidades que realizan dicha acción (por ejemplo la de la siguiente imagen, que si os la guardais y cambiais la extensión a zip ya la tendreis a mano para utilizarla...;-p)



Una vez desbloqueados, la mayoría de visores/editores de documentos PDF permiten extraer las imágenes embebidas en ellos. En estos PDFs en concreto, todo el documento es una única imagen JPG.

Una vez extraidas las imágenes JPG, nos interesa recortarlas para obtener tan sólo la zona correspondiente al mapa geológico en sí. Para ello utiizo la herramienta gratuita "jpegcrop" que se descarga de aquí:

http://jpegclub.org



Se selecciona el botón resaltado de la imagen y con el ratón se desplazan los límites del marco de la imagen hasta obtener la escena requerida. Guardamos la escena correspondiente al mapa como un nuevo archivo JPG.

Una vez tenemos todos los mosaicos del "Mapa geològic comarcal de Catalunya" como imágenes JPG, es el momento de "georreferenciarlos". Para ello utilizo el programa GlobalMapper, pues resulta muy sencilla dicha labor.

Existen 2 métodos en GlobalMapper para georreferenciar mapas:

1) Cargando primero un mapa georreferenciado (ráster o vectorial) y a continuación el otro mapa sin georreferenciar. Entonces en la ventana que aparece "Image Rectifier", se trata de escoger en el marco central con la imagen a georreferenciar, un punto que podamos localizar sin problemas en el marco izquierdo con la imagen georreferenciada, de manera que una vez seleccionado en una y otra imagen se le asignan las coordenadas correspondientes. Le damos al botón [Add Point to List] para tener al menos una lista de 4 puntos georreferenciados ("Ground Control Points") que abarquen todo el mapa, preferiblemente.

¡OJO! Debemos tener bién configurado el dátum y el tipo de coordenadas que estamos utilizando para georreferenciar el mapa, para lo cual utilizaremos el botón [Select Projection...] si tuviera que cambiarse...

Una vez tengamos la lista con los "Ground Control Points", le damos al botón [Apply] y [OK], y ya tenemos el mapa georreferenciado. Basta con acceder al "Control Panel" del GlobalMapper para cerrar el mapa que hemos utilizado para obtener las coordenadas, de manera que nos quede tan solo nuestro nuevo mapa georreferenciado. Exportamos entonces a formato ráster (recomiendo el ECW), y ¡ya está!. Para crear un único mapa con todos los mapas geológicos comarcales, los abrimos en GlobalMapper y exportamos el mosaico como nuevo mapa ECW.

2) El segundo método es el de cargar en GlobalMapper tan solo la imagen a georreferenciar, y crear la lista con los "Ground Control Points" introduciendo las coordenadas manualmente. Éste es el método a utilizar si tenemos las coordenadas disponibles en un mapa clásico "en papel" o en formato electrónico sin georreferenciar, como es el caso de los documentos PDF mencionados.

Abrimos así el archivo JPG en GlobalMapper. Un cuadro de diálogo nos indica que la imagen no está georreferenciada, y nos da la opción de si queremos georreferenciarla. Presionamos el botón [YES] y se abre la ventana para realizar el procedimiento.

Lo primero de todo es configurar el dátum y tipo de coordenadas que vamos a utilizar. En este caso vamos a utilizar coordenadas UTM de la zona 31N con dátum ED50, que son los disponibles en el PDF, seleccionándolos en el botón [Select projection...]

Abro el documento PDF con el visor correspondiente. Como los mapas incorporan la malla de la proyección UTM, resulta muy sencillo localizar 4 puntos donde se intersecten las líneas de la malla y asignarles manualmente las coordenadas UTM correspondientes. Le damos al botón [Add Point to List] para ir creando la lista con los al menos 4 "Ground Control Points".



El resto del procedimiento es igual al descrito en el método anterior... ;-p

¡Una cosa que se me olvidaba de explicar al exportar los mapas al formato ECW!

Una vez georreferenciados los mapas JPG según la trama UTM, éstos quedan entonces alineados según dicha trama y en consecuencia ligeramente rotados respecto a su forma rectangular inicial. Eso significa que al exportar de nuevo el mapa como ECW, el GlobalMapper "encuadra" el nuevo mapa para que sea rectangular, añadiendo "fondo" (negro en este caso) al mapa en sí (indicado en la imagen con las flechas amarillas), de forma que si intentamos hacer un mosaico con ellos, los mapas ECW no nos quedaran perfectamente ensamblados, en el sentido de que nos aparecerán zonas con "fondo".



Para evitarlo, el exportar a ECW (también a otros formatos ráster) existe una opción que deberemos habilitar:

[x]Make Background (Void)Pixels Transparent

Con ello conseguimos que el "fondo" añadido a los mapas sea "transparente", de manera que al ensamblar los mapas para hacer el mosaico, dicho "fondo" desaparece.



Seleccionar dicha opción tiene un inconveniente: el tamaño de cada uno de los mapas ECW resultantes puede llegar a ser el doble del que se obtiene si no se habilita dicha opción, así que recomiendo seleccionarla tan sólo si vamos a realizar un posterior ensamblado con ellos...