[GRASS-dev] Exporting maps as multilayer GeoTIFF - layer limit?

GRASS GIS GRASS GIS developer
,
1

Hi,

I struggle to export a time series as multilayer GeoTIFF: it stops at layer 145.
I have generated a time series with 365 maps and added all into a
group. Exporting that in order to export as multilayer GeoTIFF, it
fails:

GRASS 7.0.svn (nc_spm_08):~ > r.out.gdal ...
...
100%
100%
100%
Exporting raster data to GTiff format...
ERROR 6: SetColorTable() not supported for multi-sample TIFF files.
100%
...
ERROR: Unable to open input file
       </home/neteler/grassdata/nc_spm_08/user1/cell_misc/testmap.145/f_format>

Where is the trick? I have added a small script to generate a
synthetic time series for debugging.

If it matters:
GRASS 7.0.svn (nc_spm_08):~ > ulimit
unlimited

It happens in a different location as well with my original dataset.

Hints welcome,
Markus

(attachments)

grass7_multilayer_test.sh (443 Bytes)

2

Markus Neteler wrote:

I struggle to export a time series as multilayer GeoTIFF: it stops at layer
145.

Markus, your script works here! The exported file:

gdalinfo multilayer.tif -nogcp -noct -nomd | tail -5

Band 361 Block=1x1 Type=Float64, ColorInterp=Undefined
Band 362 Block=1x1 Type=Float64, ColorInterp=Undefined
Band 363 Block=1x1 Type=Float64, ColorInterp=Undefined
Band 364 Block=1x1 Type=Float64, ColorInterp=Undefined
Band 365 Block=1x1 Type=Float64, ColorInterp=Undefined

Nikos

3

Nikos,

the file is also generated here but the layers are of value 0.0
starting with layer 145...
I used gdalinfo -mm geotiff to find out. Also using latest SVN version
here, recompiled after distclean.

Best
Markus

On 11/25/13, Nikos Alexandris <nik@nikosalexandris.net> wrote:

Markus Neteler wrote:

I struggle to export a time series as multilayer GeoTIFF: it stops at
layer
145.

Markus, your script works here! The exported file:

gdalinfo multilayer.tif -nogcp -noct -nomd | tail -5

Band 361 Block=1x1 Type=Float64, ColorInterp=Undefined
Band 362 Block=1x1 Type=Float64, ColorInterp=Undefined
Band 363 Block=1x1 Type=Float64, ColorInterp=Undefined
Band 364 Block=1x1 Type=Float64, ColorInterp=Undefined
Band 365 Block=1x1 Type=Float64, ColorInterp=Undefined

Nikos

4

Markus Neteler wrote:

the file is also generated here but the layers are of value 0.0
starting with layer 145...
I used gdalinfo -mm geotiff to find out. Also using latest SVN version
here, recompiled after distclean.

Here it works:

gdalinfo multilayer.tif -nogcp -noct -nomd -mm | head -35 | tail -4

Band 1 Block=1x1 Type=Float64, ColorInterp=Gray
    Computed Min/Max=1.001,1.001
Band 2 Block=1x1 Type=Float64, ColorInterp=Undefined
    Computed Min/Max=1.002,1.002
..

gdalinfo multilayer.tif -nogcp -noct -nomd -mm | tail -5

    Computed Min/Max=1.363,1.363
Band 364 Block=1x1 Type=Float64, ColorInterp=Undefined
    Computed Min/Max=1.364,1.364
Band 365 Block=1x1 Type=Float64, ColorInterp=Undefined
    Computed Min/Max=1.365,1.365

Nikos

5

Markus Neteler wrote:

Hi,

I struggle to export a time series as multilayer GeoTIFF: it stops at layer 145.
I have generated a time series with 365 maps and added all into a
group. Exporting that in order to export as multilayer GeoTIFF, it
fails:

GRASS 7.0.svn (nc_spm_08):~ > r.out.gdal ...
...
100%
100%
100%
Exporting raster data to GTiff format...
ERROR 6: SetColorTable() not supported for multi-sample TIFF files.
100%
...
ERROR: Unable to open input file
       </home/neteler/grassdata/nc_spm_08/user1/cell_misc/testmap.145/f_format>

You were bitten by the ulimit of your system (ulimit -n). Actually it
was a bug in r.out.gdal which did not close raster maps after export.
Fixed in r58304.

Markus M

Where is the trick? I have added a small script to generate a
synthetic time series for debugging.

If it matters:
GRASS 7.0.svn (nc_spm_08):~ > ulimit
unlimited

It happens in a different location as well with my original dataset.

Hints welcome,
Markus

_______________________________________________
grass-dev mailing list
grass-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-dev