Re: [GRASS-user] Error exporting to tiff
Hi, I resolve the error. The original raster to convert was in 16bit value, and before to do the r.gdal.out is need to rescale to 8 bit. I wrongly uess the setting BYTE in r.gdal.out was to convert from 16bitto 8 bit. Instead no, is necessary a convertion using rescale before. Thx to point me in right direction. A. 2015-07-02 8:42 GMT+02:00 Andrea Peri aperi2...@gmail.com: I create it from new retrieving the settings from the original geotiff I load. Actually the region settings are this: g.region -p projection: 1 (UTM) zone: 32 datum: wgs84 ellipsoid: wgs84 north: 5058915 south: 4825185 west: 587685 east: 817515 nsres: 30 ewres: 30 rows: 7791 cols: 7661 cells: 59686851 Thx, 2015-07-02 2:24 GMT+02:00 Anna Petrášová kratocha...@gmail.com: On Wed, Jul 1, 2015 at 8:11 PM, Andrea Peri aperi2...@gmail.com wrote: Hi, I'm try-ing to export a tiff raster using this couple of commands: r.composite -c --overwrite red=refl4.hpf@PERMANENT green=refl3.hpf@PERMANENT blue=refl2.hpf@PERMANENT levels=256 output=hpf r.out.gdal -c input=hpf@PERMANENT output=/path-to/hpf.tif format=GTiff type=Byte createopt=INTERLEAVE=PIXEL,TFW=YES,PROFILE=BASELINE,COMPRESS=DEFLATE But I have an error: Checking GDAL data type and nodata value... 2..5..8..11..14..17..20..23..26..29..32..35..38..41..44..47..50..53..56..59..62..65..68..71..74..77..80..83..86..89..92..95..98..100 WARNING: Selected GDAL datatype does not cover data range. WARNING: GDAL datatype: Byte, range: 0 - 255 WARNING: Raster map hpf@PERMANENT range: 0 - 1.67114e+07 WARNING: Raster export results in data loss. Using GDAL data type Byte Input raster map contains cells with NULL-value (no-data). The value 255 will be used to represent no-data values in the input map. You can specify a nodata value with the nodata option. ERROR: Raster export aborted. The message is absolutely not clear about what is the problem. I don't know if there some log where retrieve more information. What is your region setting? Many thx for any help. Thx, -- - Andrea Peri . . . . . . . . . qwerty àèìòù - ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- - Andrea Peri . . . . . . . . . qwerty àèìòù - -- - Andrea Peri . . . . . . . . . qwerty àèìòù - ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Error exporting to tiff
I create it from new retrieving the settings from the original geotiff I load. Actually the region settings are this: g.region -p projection: 1 (UTM) zone: 32 datum: wgs84 ellipsoid: wgs84 north: 5058915 south: 4825185 west: 587685 east: 817515 nsres: 30 ewres: 30 rows: 7791 cols: 7661 cells: 59686851 Thx, 2015-07-02 2:24 GMT+02:00 Anna Petrášová kratocha...@gmail.com: On Wed, Jul 1, 2015 at 8:11 PM, Andrea Peri aperi2...@gmail.com wrote: Hi, I'm try-ing to export a tiff raster using this couple of commands: r.composite -c --overwrite red=refl4.hpf@PERMANENT green=refl3.hpf@PERMANENT blue=refl2.hpf@PERMANENT levels=256 output=hpf r.out.gdal -c input=hpf@PERMANENT output=/path-to/hpf.tif format=GTiff type=Byte createopt=INTERLEAVE=PIXEL,TFW=YES,PROFILE=BASELINE,COMPRESS=DEFLATE But I have an error: Checking GDAL data type and nodata value... 2..5..8..11..14..17..20..23..26..29..32..35..38..41..44..47..50..53..56..59..62..65..68..71..74..77..80..83..86..89..92..95..98..100 WARNING: Selected GDAL datatype does not cover data range. WARNING: GDAL datatype: Byte, range: 0 - 255 WARNING: Raster map hpf@PERMANENT range: 0 - 1.67114e+07 WARNING: Raster export results in data loss. Using GDAL data type Byte Input raster map contains cells with NULL-value (no-data). The value 255 will be used to represent no-data values in the input map. You can specify a nodata value with the nodata option. ERROR: Raster export aborted. The message is absolutely not clear about what is the problem. I don't know if there some log where retrieve more information. What is your region setting? Many thx for any help. Thx, -- - Andrea Peri . . . . . . . . . qwerty àèìòù - ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- - Andrea Peri . . . . . . . . . qwerty àèìòù - ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Exporting a raster tiff with RGB 3 bands produce a paletted tiff
Hi, I'm try-ing to export a raster RGB using this command: r.out.gdal -c input=input_image@PERMANENT output=/path-to-raster//output_image.tif format=GTiff type=Byte createopt=INTERLEAVE=PIXEL,TFW=YES,PROFILE=BASELINE,COMPRESS=DEFLATE It produce successfully a tiff , but it is a paletted tiff. Instead I like to have a true-color RGB tiff. I miss some parameter ? Thx, -- - Andrea Peri . . . . . . . . . qwerty àèìòù - ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Exporting a raster tiff with RGB 3 bands produce a paletted tiff
Hi, On Jul 2, 2015 11:58 AM, Andrea Peri aperi2...@gmail.com wrote: Hi, I'm try-ing to export a raster RGB using this command: If you have the three RGB channels, you can export via a group: http://grass.osgeo.org/grass70/manuals/r.out.gdal.html#export-r,g,b-imagery-bands-in-geotiff-format-suitable-for-esri-software : Best Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Stream extraction from pseudo-elevation map using r.watershed
Hi, I tried something that was working at least several months ago, but somehow r.watershed has changed. What I am trying to do is to extract a stream network from a pseudo elevation map. This pseudo elevation map is created by getting the distance of each stream cell from to outlet in upstream direction (using r.cost). In addition, the streams are buffered (3 cells wide) with cells of a value that is much greater than the actual stream value. This map has several features: 1) The cell values for the stream are certainly decreasing in downstream direction (=elevation is decreasing) 2) The stream is carved into a valley (important at river junctions where two streams meet). So this should actually meet the requirements for r.watershed (at least it was working some time ago) to extract the stream and get the drainage directions (e.g. for other tools like r.stream.order etc.). Now, when I try to run r.watershed (GRASS7.1) I get following warning several times: WARNING: MFD: cumulative proportion of flow distribution not 1.0 but -0.00 Although I get a drainage map as output, no output is created for the extracted stream network. I am wondering what changed and why and if this is a bug in r.watershed or if I need to change settings differently. Attached you can find the pseudo elevation map as geotiff to try it yourself! The command I am using is: r.watershed -m --o elevation=pseudo_elevation drainage=drainage_test stream=stream_test threshold=3 My system: GRASS version: 7.1.svn GRASS SVN revision: 65534 Build date: 2015-07-02 Build platform: i686-pc-linux-gnu GDAL: 1.10.0 PROJ.4: 4.8.0 GEOS: 3.4.2 SQLite: 3.7.9 Python: 2.7.3 wxPython: 2.9.4.1 Any ideas or suggestions? Best regards, Johannes pseudo_elevation.gtiff Description: Binary data ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.bioclim dataset
Tell me were can I get data to run r.bioclim module? GRASS GIS7 error: r.bioclim tmin=eurolst_clim.bio01@PERMANENT tmax=eurolst_clim.bio01@PERMANENT output=result quartals=4 ERROR: 12 maps with minimum temperatures are required regards Albert Saribekyan. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Exporting a raster tiff with RGB 3 bands produce a paletted tiff
Hi Markus, thx for hint. I tried but probably I miss something. because have an error: ERROR 6: SetColorTable() not supported for multi-sample TIFF files. ERROR 6: SetColorTable() can only be called on band 1. ERROR 6: SetColorTable() can only be called on band 1. I defined a group named gruppo and insert in it the three raster I like to have for the three bands. After I execute this command: r.out.gdal -f --overwrite input=gruppo output=D:\temp\hpf_out_rescale2.tif format=GTiff type=Byte createopt=INTERLEAVE=PIXEL,TFW=YES,PROFILE=BASELINE,COMPRESS=DEFLATE This is the results log: WARNING: Perdita di precisione: mappa raster LC81920292015157LGN00_reflect2.hpf di tipo DCELL da esportare come Byte. Può essere evitata usando Float64. WARNING: Perdita di precisione: mappa raster LC81920292015157LGN00_reflect3.hpf di tipo DCELL da esportare come Byte. Può essere evitata usando Float64. WARNING: Perdita di precisione: mappa raster LC81920292015157LGN00_reflect4.hpf di tipo DCELL da esportare come Byte. Può essere evitata usando Float64. WARNING: Si sta forzando l'esportazione del raster. Checking GDAL data type and nodata value... Using GDAL data type Byte La mappa raster di input contiene celle con valore NULL (nessun dato). Il valore 255 verrà usato per rappresentare i valori delle celle senza dati nella mappa di input. È possibile specificare un valore per le celle senza dati con l'opzione nodata. Using GDAL data type Byte La mappa raster di input contiene celle con valore NULL (nessun dato). Il valore 255 verrà usato per rappresentare i valori delle celle senza dati nella mappa di input. È possibile specificare un valore per le celle senza dati con l'opzione nodata. ERROR 6: SetColorTable() not supported for multi-sample TIFF files. Using GDAL data type Byte La mappa raster di input contiene celle con valore NULL (nessun dato). Il valore 255 verrà usato per rappresentare i valori delle celle senza dati nella mappa di input. È possibile specificare un valore per le celle senza dati con l'opzione nodata. Exporting raster data to GTiff format... ERROR 6: SetColorTable() can only be called on band 1. ERROR 6: SetColorTable() can only be called on band 1. r.out.gdal completo. File D:\temp\hpf_out_rescale2.tif created. Thx, A. 2015-07-02 12:56 GMT+02:00 Markus Neteler nete...@osgeo.org: Hi, On Jul 2, 2015 11:58 AM, Andrea Peri aperi2...@gmail.com wrote: Hi, I'm try-ing to export a raster RGB using this command: If you have the three RGB channels, you can export via a group: http://grass.osgeo.org/grass70/manuals/r.out.gdal.html#export-r,g,b-imagery-bands-in-geotiff-format-suitable-for-esri-software: Best Markus -- - Andrea Peri . . . . . . . . . qwerty àèìòù - ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.bioclim dataset
On Jul 2, 2015 5:32 PM, Albert Saribekyan albertsaribek...@rambler.ru wrote: Tell me were can I get data to run r.bioclim module? GRASS GIS7 error: r.bioclim tmin=eurolst_clim.bio01@PERMANENT tmax=eurolst_clim.bio01@PERMANENT output=result quartals=4 ERROR: 12 maps with minimum temperatures are required You can download daily data for eg Europe from http:// http://eca.knmi.nl/eca.knmi.nl/ http://eca.knmi.nl/ Likewise there are other sites for other parts of the world. Best Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Error exporting to tiff
Hi Nikos, thx for hint, I try the options Int16, UInt16, Int32 and UInt32, but always the same error. A. 2015-07-02 15:21 GMT+02:00 Nikos Alexandris n...@nikosalexandris.net: * Andrea Peri aperi2...@gmail.com [2015-07-02 09:26:31 +0200]: Hi, I resolve the error. The original raster to convert was in 16bit value, and before to do the r.gdal.out is need to rescale to 8 bit. I wrongly uess the setting BYTE in r.gdal.out was to convert from 16bitto 8 bit. Instead no, is necessary a convertion using rescale before. You should have other options (=formats) to export to so as to not sacrifice the radiometric resolution (16-bit). Did you try Int16 or UInt16? Nikos -- - Andrea Peri . . . . . . . . . qwerty àèìòù - ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user