[GRASS-dev] making r.out.gdal region sensitive

Jachym Cepicky jachym.cepicky at centrum.cz
Thu Oct 19 05:14:51 EDT 2006


hi,
On Thu, Oct 19, 2006 at 12:22:14PM +1300, Hamish wrote:
> > >> In order for the behaviour of r.out.gdal to be consistent with
> > >> other r.out.* modules,
> 
> just a check re. region bounds & grid vs. pixel conventions, does a
> full cycle result in the same exact data?
> 
> 
> r.out.gdal -r in=orig_map format=[...] out=map.ext
> 
> r.in.gdal in=map.ext out=map2
> 
> r.resample in=orig_map out=orig_map.resamp
> 
> r.mapcalc diff="orig_map.resamp - map2"
> 
> r.univar diff
> 
> ?
> Hamish

if i understand this well, you are asking for result like this:
    
d.zoom
r.out.gdal -r in=soils type=Byte out=map.ext
r.in.gdal in=map.ext out=map2 --o
r.resample in=soils out=soils.resamp --o
r.mapcalc diff="soils.resamp - map2"
r.univar diff

total null and non-null cells: 29382
total null cells: 0

Of the non-null cells:
----------------------
n: 29382
minimum: 0
maximum: 0
range: 0
mean: 0
mean of absolute values: 0
standard deviation: 0
variance: 0
variation coefficient: nan %
sum: 0

i would say, this looks good(?)

jachym

-- 
Jachym Cepicky
e-mail: jachym.cepicky at centrum.cz
URL: http://les-ejk.cz
GPG: http://les-ejk.cz/gnupg_public_key/jachym_cepicky-gpg_public_key.asc
-----------------------------------------   
OFFICE:                                     
Department of Geoinformation Technologies 
Zemedelska 3
613 00, Brno
Czech Republick
e-mail: xcepicky at node.mendelu.cz
URL:    http://mapserver.mendelu.cz
Tel.:   +420 545 134 514
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 191 bytes
Desc: Digital signature
Url : http://lists.osgeo.org/pipermail/grass-dev/attachments/20061019/1e9ffb6c/attachment.bin


More information about the grass-dev mailing list