A solution? (including resolution) - Re: [GRASS-user] Re: Help: Converting a raster map between locations (from wgs84 lat / long to UTM)

Glynn Clements glynn at gclements.plus.com
Mon Nov 3 14:55:03 EST 2008


Corrado wrote:

> this is what I have done:

> 9) set the resolution to the resolution calculated previously:
> 
> g.region nsres=927.90359439 ewres=508.14568663
> 
> 10) project the raster file:
> 
> r.proj input=dummy01 location=test12 mapset=PERMANENT 
> output=test_fromtest12_02 method=cubic
> 
> 11) set the resolution I want:
> 
> g.region nsres=2000 ewres=2000
> 
> 12) resample:
> r.resample input=test_fromtest12_02 
> output=test_fromtest12_02_resample --overwrite
> 
> What do you think?

You would probably get better results if you replaced the above four
steps with just:

	g.region nsres=2000 ewres=2000
	r.proj ...

Resampling is invariably lossy, so fewer resampling steps are
preferable to more steps.

-- 
Glynn Clements <glynn at gclements.plus.com>


More information about the grass-user mailing list