• Ir a navegación principal
  • Ir al contenido principal
  • Ir a la barra lateral primaria
  • Ir al pie de página
Logo web Geoinnova

Territorio Geoinnova - SIG y Medio Ambiente

El Blog de SIG, Territorio y Medio Ambiente

  • Web Corporativa
  • Cursos SIG
  • Cursos Medio Ambiente
  • Consultoría
  • Blog
  • Contacto
  • Inicio Blog
  • SIG
  • Programación y Desarrollo SIG
  • Teledetección
  • Medio Ambiente
  • Corporativo
SIG

Descargar tiles de OpenStreetMap con Python para su uso offline

Ene 30, 2019 · Por Asociación Geoinnova Dejar un comentario

El teléfono móvil es uno de los dispositivos desde donde más se utilizan los mapas web por medio de conexión a una red de internet. Pero por suerte o desgracia las redes de telefonía no cubren el mundo entero. Los mapas como OpenStreetMap, tienen buena información de muchas zonas rurales y espacios naturales, senderos, pistas de tierra… Descargando los tiles desde un servidor, esas imágenes de 256*256 pixeles que forman los mapas online, y colocándolos en nuestro móvil, podemos usar esos mapas sin que tengamos cobertura. Para ello necesitaremos una aplicación que nos lea esos tiles en el móvil. La aplicación osmodroid por ejemplo.

Descarga del script de Python de Github

Para descargar los tiles usamos un script de Python de Tony Rewin ubicado en GitHub.

Debemos primero agradecer que Tony haya publicado el código y modificar ligeramente el script para hacerlo más versátil . Así queda nuestro script después de alguna pequeña modificación:

#!/usr/bin/python

from sys import argv
import os
import math
import urllib.request
import random
import os.path
import time

def deg2num(lat_deg, lon_deg, zoom):
    lat_rad = math.radians(float(lat_deg))
    n = 2.0 ** zoom
    xtile = int((float(lon_deg) + 180.0) / 360.0 * n)
    ytile = int((1.0 - math.log(float(math.tan(lat_rad) + 1 / math.cos(lat_rad))) / math.pi) / 2.0 * n)
    return xtile, ytile

def download_url(zoom, xtile, ytile):
    # Switch between otile1 - otile4
    #subdomain = random.randint(1, 4)

    url = "https://a.tile.openstreetmap.org/%d/%d/%d.png" % (zoom, xtile, ytile)
    dir_path = "tiles/%d/%d/" % (zoom, xtile)
    download_path = "tiles/%d/%d/%d.png" % (zoom, xtile, ytile)

    if not os.path.exists(dir_path):
        os.makedirs(dir_path)

    if (not os.path.isfile(download_path)):
        print ("downloading %r" % url)
        source = urllib.request.urlopen(url)
        content = source.read()
        source.close()
        destination = open(download_path,'wb')
        destination.write(content)
        destination.close()
    else: print ("skipped %r" % url)



def usage():
    print ("Usage: ")
    print ("osm_tiles_downloader <lat> <lon>")

def main(argv):
    try:
        script, latW, lonW, latE, lonE = argv
    except:
        usage()
        exit(2)

    maxzoom = 15 # redefine me if need so


    # from 0 to 6 download all
    for zoom in range(0,7,1):
        for x in range(0,2**zoom,1):
            for y in range(0,2**zoom,1):
                download_url(zoom, x, y)

    # from 7 to 15 ranges
    for zoom in range(7, int(maxzoom + 1), 1):

        xtile, ytile = deg2num(latW, lonW, zoom)
        print (xtile)
        final_xtile, final_ytile = deg2num(latE, lonE, zoom)

        print ("%d:%d-%d/%d-%d" % (zoom, xtile, final_xtile, ytile, final_ytile))
        for x in range(xtile, final_xtile + 1, 1):b
            print ("%d" % x)
            for y in range(ytile, final_ytile +1, 1):
                result = download_url(zoom, x, y)


main(argv)

Modificaciones del script para descargar las Tiles

Hemos adaptado el script de Tony Rewin a Python 3. Modificamos la manera en que importamos el módulo urllib2 que nos permite acceder a direcciones url, en la línea 6 de su código. También hay que incluir entre paréntesis el texto de la función print.

En esta página de OpenStreetMap explican detalladamente el funcionamiento de los mapas web y el almacenamiento de los tiles. La función deg2num que nos convierte los grados decimales a coordenadas de ubicación del tile está en esta página.

Otra modificación que hemos hecho al script original es la manera de determinar el área a descargar.

En el código original de GitHub, línea 59, vemos cómo hace el script para ubicar el área. El script sólo nos pide un par de coordenadas, línea 44 del código original. A esas coordenadas le resta y suma 0.1 grados decimales a la longitud y le resta y suma 0.05 a la latitud. En nuestro script indicamos el área con unas coordenadas NW y SE de un área concreta. En esta web podemos fácilmente determinar las coordenadas de un área.

En la línea 45 de nuestro código definimos las variables: latW, lonW, latE, lonE. Estas coordenadas en grados decimales son las que determinarán el área a descargar.

El script descarga todos los tiles desde el zoom 0 al 6. A partir del zoom 7 descarga solamente los del área seleccionada y hasta el zoom que le indiquemos. El zoom 0 es un tile de 256 x256 px en el que entra el mundo entero. El tile ../0/0/0.png de OSM es este:

OSM con Python TilesA partir de ahí, cada nivel de zoom dobla el número de tiles tanto en horizontal como en vertical. El número de tiles para cada zoom será:

número de tiles = 2 2xZoom

Todos los tiles de un nivel de zoom determinado se almacenan en la misma carpeta. Dentro de esta carpeta todos los tiles con la misma x se almacenan en la misma carpeta. Así, el tile ../4/3/9.png estará guardado en la carpeta 3 del zoom 4.

En esta página de maptiler se puede jugar con la ubicación y las coordenadas de los tiles y entender un poco mejor cómo funciona un mapa web.

Los mapas web usan la proyección Web Mercator. Se extendió su uso cuando Google Maps la adoptó in 2005. El sistema de coordenadas es el WGS84. La controversia con la proyección Web Mercator es debido a que asume la forma esférica y no la elipsoidal como el sistema se coordenadas WGS84.

La función download_url nos crea la estructura de carpetas donde irán ubicados los tiles según se vayan descargado.

El loop for zoom in range nos origina las url a descargar según los valores las coordenadas de nuestras variables. Le indicamos el nivel de zoom máximo que deseamos. A nivel 15 aparecen claramente visibles las pistas y caminos simbolizadas en OSM. Cada nivel cuadruplica en peso al anterior. Así que con un nivel 16 cuadruplicaríamos los datos a cargar en el móvil.

Tomando como muestra una zona como la de Bardenas Reales en Navarra el comando a ejecutar en la consola de Windows quedaría así:

python downloadOSMtiles.py 42.357339 -1.721082 41.958277 -1.276331

El script nos creará la estructura de carpetas organizadas por nivel de zoom en el directorio donde esté ubicado el script.

Y con mbutil podemos convertir los tiles que hemos generado a un único archivo .mbtile pero eso ya es otra historia.

Aprende con Geoinnova

Si quieres extender tus conocimientos de programación GIS en entornos web y diferenciarte profesionalmente, pégale un vistazo a nuestro Curso de Aplicaciones Web Mapping con GDAL, Geoserver y Leaflet, donde profesionales con gran experiencia en desarrollo te guiarán en el proceso de aprendizaje y de aplicación de las herramientas más conocidas en web mapping.

Autor: Gonzalo Echeverría

Etiquetas: web mapping

Artículos relacionados

Integración de Terre3 con OL6
Como integrar la API 3D TR3.js sobre OpenLayers 6
Ene 4, 2021
Programación y Desarrollo SIG
web mapping
Visor NATURALMENTE BENAMAHOMA
Trabajo Fin de Máster TGIS de Geoinnova: Visor web ayuntamiento de Grazalema
Nov 22, 2020
Corporativo, SIG
qgis, web mapping
Nuevo Curso de Desarrollo de aplicaciones de mapas web en 3D con tecnologías de código abierto (ThreeJS)
Nov 22, 2020
Formación
3D, threeJS, web mapping
Avatar

Asociación Geoinnova

La Asociación Geoinnova fue creada en 2009 por profesionales del medio ambiente y del territorio. Tiene como objetivos principales desarrollar su actividad en los ámbitos de la Formación en SIG y Medio Ambiente, así como en la Consultoría de proyectos ligados al Urbanismo, Ordenación del Territorio, Medio Ambiente y Sistemas de Información Geográfica.

Interacciones con los lectores

Deja una respuesta Cancelar la respuesta

Tu dirección de correo electrónico no será publicada. Los campos obligatorios están marcados con *

Barra lateral primaria

Suscríbete al boletin

    Si te apetece puedes dejarnos tu nombre para que los correos tengan un trato personalizado

    Mediante el envío de mis datos personales confirmo que he leído y acepto la política de privacidad

    RESPONSABLE: Asociación Geoinnova. FINALIDAD: envío de publicaciones, promociones e información sobre cursos y eventos. LEGITIMACIÓN: tu legítimo consentimiento. DESTINATARIOS: Active Campaign con titular Active Campaign LLC, alojada en EEUU y suscrita al EU PrivacyShield. DERECHOS: acceso, rectificación, limitación, supresión de los datos (en [email protected])y a presentar reclamación ante una autoridad de control. INFORMACIÓN ADICIONAL: Política de privacidad.

    La último

    Creación de variable

    Personalización de QGIS con Python

    Apps para identificar cimas y montañas con el teléfono móvil

    Estadísticas catastrales y tributarias

    Como trabajar con estadísticas Catastrales y Tributarias de España

    Agenda de cursos en Enero y Febrero

    R y su transformador RCaller

    Ejecutando R en FME Desktop

    La educación ambiental es fundamental para preservar el medioambiente.

    ¿Qué es la educación ambiental?

    SAGA

    Interpolación Kriging con SAGA

    Footer

    Sobre Nosotros

    Territorio Geoinnova pretende ofrecer noticias y formación sobre el medio ambiente y las tecnologías de la información geográfica con interés a dichos sectores profesionales.

    Legal

    • Nota Legal

    El blog y todo su contenido se encuentra bajo una licencia de Creative Commons Reconocimiento-NoComercial-CompartirIgual 3.0 Unported.

    Todas las categorías

    • Actualidad Ambiental (85)
    • Alimentación Saludable (17)
    • Corporativo (121)
    • Formación (81)
    • Geolibrería (40)
    • Medio Ambiente (948)
    • Patrocinado (12)
    • Programación y Desarrollo SIG (12)
    • SIG (487)
    • Tecnología (72)
    • Teledetección (19)

    © 2021 · Desarrollada por Juan María Arenas - OikosMSP