[GRASSLIST:10223] Re: Problems with r.proj - grass 6.01

Morten Hulden morten at untamo.net
Thu Feb 9 08:48:12 EST 2006


Hi Helmut,

I suspect your input map and output region are only overlapping by a few 
cells (2 colums with 4 rows each). Therefore your map appears empty, but 
I think you should see something in 8 cells in one corner.

Both the imput map and output location are in tmerc projection with 
lat_0=0. But:

Input North: 5321295.600000
Output South: 5353002.550472

shows that the northmost latitude of the input map does not reach the 
southmost latitude of the output region.

You could try to extend your output region, at least soutward, and try 
again.

rgds
Morten


Helmut Kudrnovsky wrote:
> there is no difference between the options "with region cropping" (look at 1) and "without region cropping" (look at 2)
> 
> i tried this procedure with different output-projections which are matching for germany/austria and different region-settings and different interpolation methods (cubic, bilinear).

> Input Projection Parameters: +proj=tmerc +lat_0=0 +lon_0=13.33333333333333 +k=1.000000 +x_0=450000 +y_0=0 +no_defs +a=6377397.155 +rf=299.1528128 +towgs84=577.326,90.129,463.919,5.1366,1.4742,5.2970,2.4232
> Input Unit Factor: 1
> 
> Output Projection Parameters: +proj=tmerc +lat_0=0 +lon_0=12 +k=1.000000 +x_0=4500000 +y_0=0 +no_defs +a=6377397.155 +rf=299.1528128 +towgs84=597.1,71.4,412.1,0.894,0.068,-1.563,7.58
> Output Unit Factor: 1
> Input:
> Cols:	2349 (2349)
> Rows:	2347 (2347)
> North: 5321295.600000 (5321295.600000)
> South: 5318948.400000 (5318948.400000)
> West:  4561644.000000 (4561644.000000)
> East:  4563992.800000 (4563992.800000)
> ew-res:	0.999915
> ns-res:	1.000085
> 
> Output:
> Cols:	4373 (2)
> Rows:	3210 (4)
> North: 5356212.272416 (5356212.272416)
> South: 5353002.550472 (5353002.550472)
> West:  8604787.256327 (8604787.256327)
> East:  8609160.077374 (8609160.077374)
> ew-res:	0.999959
> ns-res:	0.999913




More information about the grass-user mailing list