Estoy tratando de generar mosaicos con gdal2tiles.py y los mosaicos resultantes se almacenan como 18/132639/164219. Si hago lo mismo con MapTiler (la variante GUI de gdal2tiles), el resultado del mosaico es 18/132639/97923. Quiero los resultados de MapTiler (97923 ), pero solo con gdal2tiles, ya que quiero hacer esto en un proceso por lotes y no con una GUI.

Ya probé muchas cosas, deformando, traduciendo, lo que sea, y se está poniendo muy molesto, especialmente porque dice que MapTiler es solo una interfaz gráfica para gdal2tiles y hace básicamente lo mismo.

gdal2tiles.py -z 17-19 -s EPSG:900913 -p raster overlay-level1.png gdal2tiles.py -z 17-20 -s EPSG:3785 overlay-level1.png gdal2tiles.py -p geodetic -s EPSG:4326 -z 18-21 -w none overlay-level1.png 

La mayoría de ellos dan como resultado un ERROR 5 Access window out of range in RasterIO() o ERROR 5: Illegal values for buffer size.

Este funciona, pero da como resultado valores de Y incorrectos:

gdal2tiles.py -s EPSG:4326 -z 16-19 -w none overlay-level1.png 

¿Alguna idea de qué pasa ahí?

Comentarios

  • Intenté producir mosaicos para publicar mi mapa en línea. Quiere obtener la misma estructura que produjo la última versión 0.5.2 de maptiler pero quiere hacerlo usando el script gdal. He probado – último maptiler para comparar cómo debería verse – versión antigua de maptiler – gdal2tiles.py 25611 2013-02-07 10: 24: 56Z bishop $ (del paquete qgis valmiera) – gdal2tiles.py 19288 2010-04- 02 18: 36: 17Z rouault $ + su modificación (Andre Joost) http://gis.stackexchange.com/questions/63024/gdal2tiles-maptiles-from-bsb-kap-are-switched

Respuesta

La solución es mucho más simple de lo esperado. Simplemente cambie la configuración de TMS a true en el folleto.

L.tileLayer("http://{s}.somedomain.com/{foo}/{z}/{x}/{y}.png", {tms: true}); 

Comentarios

  • Tenga en cuenta que gdal2 gdal2tiles.py (GDAL 2.1.3, released 2017/20/01) utiliza el " correcto " o la convención de Google / OSM para x vs y, por lo que los mosaicos generados pueden no necesitar tales correcciones.
  • @AhmedFasih, todavía no funciona. probado desde GDAL 2.3.1

Respuesta

gdal2tiles usa la especificación TMS para nombrar mosaicos (de abajo hacia arriba) , mientras que Maptiler le permite usar el sistema de nombres de mosaicos de Google / Openstreetmap (de arriba a abajo).

Vea mi respuesta a este tema sobre cómo modificar gdal2tiles.py para los nombres de mosaicos de Google / OSM:

GDAL2Tiles: MapTiles de BSB / KAP se cambian

Comentarios

  • ¿Hice los cambios, pero los números Y siguen siendo los mismos? Estoy ' trabajando con CRC EPSG: 4326, que probablemente también sea un problema, ya que el folleto requiere EPSG: 3857.
  • Debería haber reemplazado ty y y por ty2 e y2 en las líneas existentes citadas también. Mi fuente y mosaicos estaban todos en EPSG: 3857.
  • Oh, ok. Bueno, ya que no ' t quiero cambiar gdal2tiles.py en todas las máquinas involucradas, la variante tms: true es la mejor solución. ¡Gracias!

Responder

Una contribución tardía, pero quizás ayude …

Como alternativa a cambiar el código en gdal2tiles.py o agregar la bandera especial al folleto, se pueden cambiar los nombres de archivo de los mosaicos. Esta fue una solución fácil para crear mosaicos para usar con ArcGIS online. Aquí hay una secuencia de comandos bash aproximada que alternará la convención z / x / y entre OSM / Google-Maps (gdal2tiles predeterminado) y TMS (usado en ArcGIS Online)

 # Each time it runs it will toggle the convention between OSM/Google-Maps and TMS tiledir=$1 for x in $(find $tiledir -name "*.png") do #echo $x # isolate where tile folders begin tiles="${x#$tiledir/}" #echo $tiles # separate each z,x,y value zdir=$(echo $tiles | cut -d"/" -f1-1) xdir=$(echo $tiles | cut -d"/" -f2-2) ypng=$(echo $tiles | cut -d"/" -f3-3) # remove extension from y.png files y="${ypng%.*}" # calculate new y using formula referenced here: # http://osgeo-org.1560.x6.nabble.com/gdal2tiles-tiles-in-wrong-hemisphere-and-or-Openlayers-problem-td3742809.html (( newy=(2**zdir - 1)-y )) #echo $zdir,$xdir,$y $newy is newy # rename the files mv -v $x $tiledir/$zdir/$xdir/$newy.png done 

Respuesta

Puede descargar el gdal2tiles.py para Google / OSM mosaicos por vínculo

http://gisfile.com/user/admin/gdal2tiles.py

Comentarios

  • Bienvenido a GIS @SE, gracias por su respuesta, sin embargo, debe explicar en su respuesta cuál es la diferencia en su enlace de descarga recomendado en comparación con la pregunta original, para que los futuros usuarios también lo entiendan. Saludos cordiales

Deja una respuesta

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