[gdal-dev] Gdaladdo spoils rasterlite with multible raster tables

[Jukka Rahkonen]

Hi,

I made some further tests with GDAL 1.9.0, gdal_translate and 
gdaladdo and I start to be pretty sure that gdaladdo is poisoning 
the rasterlite databases which contain several raster tables. 
I can put the original sample images somewhere for download if requested.

QGis shows all the rasters OK after the initial gdal_translate phase
but not any more after running gdaladdo. By looking at the gdalinfo 
listins it is not a wonder.

Building rasterlite DB from 3 originals
=======================================
gdal_translate -of rasterlite -a_srs epsg:2393 18m.tif
rasterlite:multitest.sqlite,table=raster_1
Input file size is 4000, 4000

gdal_translate -of rasterlite -a_srs epsg:2393 18n.tif
rasterlite:multitest.sqlite,table=raster_2
Input file size is 4000, 4000

gdal_translate -of rasterlite -a_srs epsg:2393 18mb.tif
rasterlite:multitest.sqlite,table=raster_3
Input file size is 2000, 2000

Check gdalinfo; everything OK 
=============================
gdalinfo rasterlite:multitest.sqlite,table=raster_1
Driver: Rasterlite/Rasterlite
Files: multitest.sqlite
Size is 4000, 4000
Origin = (3180000.000000000000000,6640000.000000000000000)
Pixel Size = (20.000000000000000,-20.000000000000000)

gdalinfo rasterlite:multitest.sqlite,table=raster_2
Driver: Rasterlite/Rasterlite
Files: multitest.sqlite
Size is 4000, 4000
Origin = (3260000.000000000000000,6640000.000000000000000)
Pixel Size = (20.000000000000000,-20.000000000000000)

gdalinfo rasterlite:multitest.sqlite,table=raster_3
Driver: Rasterlite/Rasterlite
Files: multitest.sqlite
Size is 2000, 2000
Origin = (3180000.000000000000000,6640000.000000000000000)
Pixel Size = (40.000000000000000,-40.000000000000000)

Add overviews
=============

gdaladdo rasterlite:multitest.sqlite,table=raster_1 2 4 8
0...10...20...30...40...50...60...70...80...90...100 - done.
0...10...20...30...40...50...60...70...80...90...100 - done.
0...10...20...30...40...50...60...70...80...90...100 - done.

gdaladdo rasterlite:multitest.sqlite,table=raster_2 2 4 8
0...10...20...30...40...50...60...70...80...90...100 - done.
0...10...20...30...40...50...60...70...80...90...100 - done.
0...10...20...30...40...50...60...70...80...90...100 - done.

gdaladdo rasterlite:multitest.sqlite,table=raster_3 2 4 8
0...10...20...30...40...50...60...70...80...90...100 - done.
0...10...20...30...40...50...60...70...80...90...100 - done.

Check gdalinfo; all but first table have wrong pixel size
=========================================================

gdalinfo rasterlite:multitest.sqlite,table=raster_1
Driver: Rasterlite/Rasterlite
Files: multitest.sqlite
Size is 4000, 4000
Origin = (3180000.000000000000000,6640000.000000000000000)
Pixel Size = (20.000000000000000,-20.000000000000000)

gdalinfo rasterlite:multitest.sqlite,table=raster_2
Driver: Rasterlite/Rasterlite
Files: multitest.sqlite
Size is 4000, 4000
Origin = (3260000.000000000000000,6640000.000000000000000)
Pixel Size = (40.000000000000000,-40.000000000000000)

gdalinfo rasterlite:multitest.sqlite,table=raster_3
Driver: Rasterlite/Rasterlite
Files: multitest.sqlite
Size is 2000, 2000
Origin = (3180000.000000000000000,6640000.000000000000000)
Pixel Size = (80.000000000000000,-80.000000000000000)

-Jukka Rahkonen-