Hi Anna, many thx for your reply. You are right.
The region is set to 30meters. I was confused from the grass display that show me the pixel cell size to 15meters. Thx again. Andrea. 2015-07-27 7:43 GMT+02:00 Anna Petrášová <[email protected]>: > > On Jul 26, 2015 12:05 PM, "Andrea Peri" <[email protected]> wrote: >> >> Hi, >> >> Starting on grass 7.0. >> I try to use the >> i.fusion.hpf script to create a new panchromatic image. >> >> Starting with a raster (from landsat) 30meters and a panchromatic >> 15meters. >> With the i.fusion.hpf I produce a new raster at high resolution >> (15meters). >> >> The verify that really it has 15m , I open it with the grass display >> and measure che legth of the pixel. >> Ok, it is really 15meters. >> >> After this I try to export the new raster with the r.gdal.out. >> This is the command I run: >> >> (the LC81920302015157LGN00_B4.org_hpf@andrea is the new 15meter raster) >> >> r.out.gdal -f --verbose input=LC81920302015157LGN00_B4.org_hpf@andrea >> output=D:\temp\exported_image.tif format=GTiff type=Float64 >> createopt=INTERLEAVE=PIXEL,TFW=YES,PROFILE=BASELINE,COMPRESS=DEFLATE >> >> This is the log of execution: >> >> r.out.gdal -f --verbose input=LC81920302015157LGN00_B4.org_hpf@andrea >> output=D:\temp\exported_image.tif format=GTiff type=Float64 >> createopt=INTERLEAVE=PIXEL,TFW=YES,PROFILE=BASELINE,COMPRESS=DEFLATE >> Esportazione al tipo di dato GDAL: Float64 >> Checking GDAL data type and nodata value... >> ERROR 6: SetColorTable() only supported for Byte or UInt16 >> bands in TIFF format. >> Using GDAL data type <Float64> >> Input raster map contains cells with NULL-value (no-data). The value >> nan will be used to represent no-data values in the input map. You can >> specify a nodata value with the nodata option. >> Exporting raster data to GTiff format... >> WARNING: Too many values, color table cut to 65535 entries >> r.out.gdal completo. File <D:\temp\exported_image.tif> created. >> >> I notice the warning for too many values. >> But really I don't understand what it mean because this is a float >> raster and not an integer raster. >> >> When I go to test the new exported tif image with qgis. >> I see the pixel size is now 30meters and not 15meters. >> >> How is possibile that the r.gdal.out could reduce the resolution of the >> image ? >> >> This is really uncomprensible to me. >> >> Is this a knowed issue ? >> >> Any help is welcome. > > Is it possible your computational region is set to 30 m resolution? Try > g.region -p. r.out.gdal respects region settings. >> >> Thx. >> >> -- >> ----------------- >> Andrea Peri >> . . . . . . . . . >> qwerty àèìòù >> ----------------- >> _______________________________________________ >> grass-user mailing list >> [email protected] >> http://lists.osgeo.org/mailman/listinfo/grass-user -- ----------------- Andrea Peri . . . . . . . . . qwerty àèìòù ----------------- _______________________________________________ grass-user mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-user
