[gdal-dev] gdal2tiles.py crashes

[Klokan Petr Přidal]

Hi Sven,

It seems you are using old version of GDAL2Tiles, from GDAL 1.5, there
is a new version available in GDAL 1.6, which supports Google Maps
tiles.

It is not clear how your data look-like, I guess the Erdas image you
are trying to tile is using 16-bit values in the raster.
(Seeing gdalinfo on your input would help me to be sure.)

The solution in such case is to preprocess your input raster by calling:

gdal_translate -ot Byte -of vrt input.img temp.vrt

and run gdal2tiles on the temp.vrt.

Best regards

Klokan Petr Pridal

On Sun, Dec 7, 2008 at 3:42 PM, Sven Geggus <lists at fuchsschwanzdomain.de> wrote:
> I try to generate Google compatible tiles from one big rasterimage (EPSG:4326).
>
> This does however not work:
>
> Input (out-4326.img):
> ================================================================================
>  Driver: HFA / Erdas Imagine Images (.img)
>  Size: 87231 x 57321 x 3
>  Projection: GEOGCS["WGS_1984",DATUM["WGS_1984",SPHEROID["WGS
>  84",6378137,298.2572235630016],TOWGS84[0,0,0,0,0,0,0]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]]
>  NSEW:  (50.068567784495805, 48.740330987425509, 13.196098581605034,
>  11.174790090035813)
> Projection detected as EPSG:4326
> Generating of KML is possible
> Output (../lvg-by-slippymap/):
> ================================================================================
>  Format of tiles: PNG / Portable Network Graphics
>  Size of a tile: 256 x 256 pixels
>  Count of tiles: 101974
>  Zoom levels of the pyramid: 9
>  Pixel resolution by zoomlevels: [0.011864015633013961,
>  0.0059320078165069805, 0.0029660039082534903,
>  0.0014830019541267451, 0.00074150097706337256,
>  0.00037075048853168628, 0.00018537524426584314,
>  9.268762213292157e-05, 4.6343811066460785e-05,
>  2.3171905533230393e-05]
> --------------------------------------------------------------------------------
> Zoom 9 - pixel 0.00002317190553323039 131072
> --------------------------------------------------------------------------------
> 1/101974 ../lvg-by-slippymap/9/0/0.png [0, 0] [0, 57065] [256, 256]
> [256, 256]
> Traceback (most recent call last):
>  File "/usr/bin/gdal2tiles.py", line 875, in ?
>    writetile( filename, data, dxsize, dysize, bands)
>  File "/usr/bin/gdal2tiles.py", line 75, in writetile
>    alpha.WriteArray( alphaarray, 0, tilesize-dysize )
>  File "/usr/lib/python2.4/site-packages/osgeo/gdal.py", line 782, in WriteArray
>    return gdalnumeric.BandWriteArray( self, array, xoff, yoff )
>  File "/usr/lib/python2.4/site-packages/osgeo/gdal_array.py", line 157, in BandWriteArray
>    raise ValueError, "array does not have corresponding GDAL data type"
> ValueError: array does not have corresponding GDAL data type
>
> Regards
>
> Sven
>
> --
> "In my opinion MS is a lot better at making money than it is at making good
> operating systems" (Linus Torvalds, August 1997)
>
> /me is giggls at ircnet, http://sven.gegg.us/ on the Web
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev
>