viernes, 19 de septiembre de 2008

Diferencias entre ETRS89, WGS84 y ED50 S&P

En la actualidad existen infinidad de servidores WMS desde los que se pueden descargar mapas de todo tipo en diferentes proyecciones y datums. En España los mas comunes son el ETRS89 (nuevo sistema de referencia geodésico oficial de España), el ED50 S&P (este es el antiguo sistema oficial y en el que están la mayor parte de los mapas de España) y el WGS84 (en el que funcionan los satélites GPS). En la página del ministerio de fomento pone que aunque "WGS84 y ETRS89 son equivalentes para la gran mayoría de aplicaciones topográficas o cartográficas", "ETRS89 es el sistema de referencia geocéntrico oficial en Europa de precisiones mucho más elevadas que la última solución WGS84" pero ¿como se cuantifica eso?.

Para calcular las diferencias entre ETRS89 y WGS84 he generado una red regular de coordenadas geográficas (lat/lon) en el datum ETRS89 y con el programa Global Mapper lo he transformado a WGS84. Después he calculado la distancia angular entre los puntos de ambas redes y se obtiene el siguiente resultado (pinchar en la imagen para verla mas grande):

La escala está en micro segundos de arco (μ arcsec).

Como se ve en la imagen, en los polos y el ecuador no hay ninguna diferencia entre WGS84 y ETRS89, y donde la diferencia es mas grande es a 45º N y S. Si consideramos que el radio medio de la tierra es de 6371km, un segundo de arco en la superficie equivale a 30.88m y un μ arcsec seria 0.03mm. La diferencia máxima que hay en el mapa es de 3.2 μ arcsec que es 0.1mm, osea, completamente despreciable. Por lo tanto se puede decir que a efectos prácticos ETRS89 y WGS84 son completamente equivalentes en coordenadas geográficas. La transformación entre coordenadas la he hecho con Global Mapper, y tampoco puedo estar completamente seguro de su exactitud.

Si en vez de usar coordenadas geográficas uso la proyección UTM en la zona 30, los resultados salen algo diferentes pero muy parecidos:

Diferencias entre ETRS89 y WGS84:

Dentro de la zona UTM30, las diferencias son menores a 1mm. Fuera de la zona 30 (zonas 29 y 31), las diferencias no llegan a los 10cm. El hecho de que no salga lo mismo usando coordenadas geográficas o coordenadas en UTM creo que es por usar una zona UTM (en este caso la 30) en zonas diferentes (29 y 31). Como se puede ver en la imagen las diferencias van aumentando según nos alejamos de la zona 30.

El sistema ED50 S&P es muy diferente a los dos anteriores. He repetido la transformación de antes pero ahora entre ETRS89 y ED50 S&P y en la proyección UTM zona 30.


Las diferencias ahora están entre 220m y 240m según la zona de España.

Conclusión: A la hora de calibrar mapas, es muy importante saber si está en ETRS89 o ED50 S&P, pero da lo mismo si está en WGS84 o ETRS89.

sábado, 12 de julio de 2008

True Marble en ECW. Imagen del planeta de 15Gpix!!

La empresa Unearthed Outdoors ha producido un una serie de imágenes del planeta con colores realistas y a una resolución de 15m/pix a partir de fotografías de landsat. Estas imágenes son bastante caras, pero en su página web tienen también imágenes gratuitas a una resolución de 250m/pix.


Las imágenes están en formato GeoTIFF y ocupan una barbaridad. Todo el planeta esta dividido en 32 cuadrantes y al ser imágenes sin comprimir, el fichero de cada cuadrante ocupa 1GB!!! (tienen 21600x21600pix en la resolución de 250m/pix), por lo que el mapa completo son 32GB.. Para ahorrar espacio y tener algo mas manejable, he descargado todos los cuadrantes y he generado un único fichero en formato ecw de 500MB (con compresión 1:10). La imagen resultante tiene un tamaño de 172800x86400 pix, o lo que es lo mismo 15Gpix. Este único fichero se puede abrir con el visor de imágenes .ecw ER Viewer y es fácilmente manejable por programas de GPS o GIS como Oziexplorer, CompeGPS o Global Mapper. También he generado la imagen a otras resoluciones (500m, 1km, 2km, 4km, 8km, 16km y 32km/pix). He puesto todo en rapidshare para que sea fácil su descarga. Estos son los links:

http://rapidshare.com/users/B8S1R/1

El password para descomprimir los ficheros es: oxbile

En el mismo sitio también tienen mapas de usos del suelo utilizando dos esquemas de colores:

NLCD 1992



y GLC2000.



Para estos mapas también he generado los ficheros ECW correspondientes, y se pueden descargar desde los siguientes links:

Ficheros de los colores GLC:
http://rapidshare.com/users/B8S1R/2

Ficheros de los colores NLCD:
http://rapidshare.com/users/B8S1R/3

El password para descomprimir todos los ficheros es: oxbile

Las siguientes imágenes son algunos detalles de zonas concretas (pinchar para verlas a tamaño completo):

Fiordos en Noruega.


Islas Canarias.


Gran Cañón del Colorado.


Monte Everest.

miércoles, 9 de julio de 2008

Mapas Topográficos MTN50 del IGN

El Instituto Geográfico Nacional (IGN) ha puesto en su página web los mapas topográficos 1:50000 para su descarga gratuita. Según he podido leer en alguna noticia, el IGN ha liberado la información geográfica digital para uso no comercial, y lo único que hay que hacer es mencionar el origen y propiedad de los datos.

Los mapas tienen un tamaño A1 y están en formato pdf. También hay mapas en A3 pero la calidad es muy mala. En algunos de ellos hay varias versiones con fechas desde 1875 hasta el 2005. Esta serie temporal que hay en algunas hojas es interesante para ver la evolución y crecimiento de algunas ciudades, por ejemplo de Madrid. Las siguientes imágenes son un detalle de la hoja 559 en los años 1875, 1916, 1944, 1969 y 2002 (pinchar en cada imagen para verla mas grande).


En total hay 2900 ficheros pdf de las 1130 hojas que tiene la serie MTN50 y ocupan 8.3 GB

Aunque en realidad el 99% de los mapas no los voy a necesitar nunca, los he descargado todos por dos motivos:
1.- Las paginas web van apareciendo y desapareciendo con el tiempo. Puede que cuando quiera consultar alguno de los mapas, el servicio ya no esté disponible.
2.- Me gusta coleccionar mapas :-)

Dado que es bastante pesado bajarse todos los mapas, los voy a ir poniendo en rapidshare por si alguien quiere descargarselos en bloque. Los he organizado en ficheros de 200 hojas y comprimidos con winrar. En cada "pack" hay también un fichero de recuperación por si alguno de los otros ficheros desaparece o se corrompe. Los mapas los he pasado a jpg para que sea facil editarlos y calibrarlos con ozi u otro programa. Si alguien necesita solo unos pocos mapas, es mejor que se los baje directamente de la pagina del IGN.

La siguiente imagen es un ejemplo de una hoja completa. hacer click sobre la imagen para verla mas grande.



Hojas de la 1 a la 200 en:
http://rapidshare.com/users/B8S1R/4
Hojas de la 201 a la 400 en:
http://rapidshare.com/users/B8S1R/5
Hojas de la 401 a la 600 en:
http://rapidshare.com/users/B8S1R/6
Hojas de la 601 a la 800 en:
http://rapidshare.com/users/B8S1R/7
Hojas de la 801 a la 1000 en:
http://rapidshare.com/users/B8S1R/8
Hojas de la 1001 a la 1130 en:
http://rapidshare.com/users/B8S1R/9
pass: oxb