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