Transformación al sistema ETRS89 de ortofotos con rejilla NTv2 para la Comunidad Valenciana con Quantum GIS 2.2.0

Transformación al sistema ETRS89 de ortofotos con rejilla NTv2 para la Comunidad Valenciana con Quantum GIS 2.2.0 En este manual vamos a explicar una serie de pasos para realizar la transformación de ED50 a ETRS89 de ortofotos con software libre, en este caso, mediante QGIS Valmiera (v. 2.2.0) con sus librerías GDAL. Para ello necesitamos: •

Descargar la Rejilla para cambio de Datum entre ED50 y ETRS89 (en formato NTV2) desde la página web del IGN: http://www.ign.es/ign/layoutIn/herramientas.do

•

Guardar el archivo sped2et.gsb o PENR2009.gsb (península) y BALR2009.gsb (Baleares) para la versión 2009 de la transformación en un directorio destinado, como por ejemplo: • Versión 2.2.0: C:\Archivos de programa\QGIS Valmiera\share\proj • En el caso de la C. Valenciana, la transformación realizada por el fichero .gsb es válida ya que en su cálculo se ha empleado una selección de datos consistentes entre si de las redes y la densificación que coexisten en la zona

• Es recomendable el uso de la versión de rejilla más consistente con los datos origen o del proyecto. • Los shapes, coberturas o capas deben tener los sistemas de referencia y proyección de los datos origen (ED50) y destino (ETRS89) definidos de forma correcta.

1.

Antes de empezar recomendamos leer el manual Transformación NTv2 para la Comunidad Valenciana con Quantum GIS (v 2.2.0 Valmiera), donde se explica cómo se integran la rejillas dentro de la última versión de Quantum GIS. [http://goo.gl/llh22n] Primero cambiamos las opciones del Proyecto, para activar la reproyección al vuelo y las transformaciones de Datum predeterminadas.

•Activar reproyección al vuelo por defecto •El SRC del proyecto debe ser: ETRS89 UTM / zone 30N (o 31N) •Preguntar SRC cuando se carga una capa •“Preguntar la transformación de DATUM cuando no haya una predefinida”.

Esto permitirá hacer la transformación cuando detecte que la capa cargada no tiene el mismo sistema de referencia que el proyecto

3. Añadimos la ortofoto en ED50, que en este caso es en formato tif, y seleccionamos el sistema de referencia ED50/UTM zone 30N, que es el sistema en el que tenemos la ortofoto.

4. En primer lugar vamos a explicar como hacer la transformación al vuelo. Una vez seleccionado el sistema de referencia, aparecerá una ventana que nos preguntará que transformación de DATUM queremos para la capa.

5. Seleccionamos la que más nos convenga, (PENR2009.gsb), para la Comunidad Valenciana, y QGIS realizará la transformación. De este modo ya está hecha la transformación al vuelo. Ahora sólo hace falta guardar la capa con su nuevo sistema de referencia.

6. Para guardar la ortofoto en el nuevo sistema de referencia, usamos las librerías GDAL, que es una biblioteca de software para la lectura y escritura de formatos geoespaciales, y la librería PROJ4, que es una librería de proyecciones cartográficas. Vamos a Ráster/Proyecciones/Combar (Reproyectar). Una vez abierta la herramienta, pulsamos el botón editar.

7. Hay que editar el código, para que realice la transformación: gdalwarp -s_srs "+init=epsg:23030 +nadgrids=C:\RUTA\REJILLA\PENR2009.gsb +wktext" t_srs "+init=epsg:25830 +nadgrids=null +wktext" -of GTiff C:/RUTA/ORTOED50/ED50Senija_2.tif C:/RUTA/ORTOETRS89/ETRS89Senija_3.tif

Por orden se indica: gdalwarp: utilidad GDAL para la transformación de archivos raster (ortofotos, mdt, etc) a nuevos sistemas de referencia. -s_srs: sistema referencia origen(ED50), debe introducirse según su código EPSG +nadgrids: ruta donde se encuentra el fichero rejilla .gsb -t_srs: sistema de referencia destino (ETRS89) , debe introducirse según su código EPSG -of Gtiff: indica el formato de salida. Las dos últimas rutas indican la ruta del fichero de entrada y el de salida respectivamente.

Para más información sobre las librerías PROJ4 y GDAL: http://trac.osgeo.org/proj/ http://www.gdal.org/

8. Una vez realizada la transformación, QGIS crea la nueva ortofoto con un SRC personalizado equivalente al ETRS89. Usando la herramienta GDAL Translate: Se le asigna el sistema destino (SRC) de forma definitiva.

El tiempo empleado es de aproximadamente 5 minutos para una una ortofoto de 40 Mb en un ordenador de características mínimas: • Intel Core 2 4300 @ 1.86GHz – 2Gb RAM • Windows XP

• Usar la herramienta GDAL Translate: 1.

2.

3. 1. Seleccionamos la imagen creada. 2. Elegimos el SRC que queremos: 25830 (ETRS89 30 N) 3. Elegimos el fichero de salida

Una vez realicemos este último paso, el proceso de transformación se habrá terminado. Para comprobar que la transformación se ha realizado correctamente, bastará con superponer capas vectoriales en ETRS89 o utilizar puntos de control en ETRS89.