[GRASS-user] Projection Units and Values Negatively Affect r.topidx Results

[Rich Shepard]

On Thu, 18 Feb 2010, Micha Silver wrote:

> If I followed your post correctly, you should not be converting anything 
> during the inport. You import into a Location defined with projection params 
> exactly like the original ArcInfo file.

Micha,

   I realized that I need to use r.proj (and v.proj) to make the conversions.

> Then create a new, empty Location with units meters.

location=oregonM, mapset=PERMANENT PROJ_UNITS:
oregonM/PERMANENT/PROJ_UNITS 
unit: meter
units: meters
meters: 1

> Switch to the new location, import the region vector polygon (from
> v.in.region in the original Location), and set the region to that vector.

Source location:

GRASS 6.4.0svn (Oregon):/usr4/grassbase > v.in.region out=coverage type=area
Building topology for vector map <coverage>...
Registering primitives...
2 primitives registered
6 vertices registered
Building areas...
  100%
1 areas built
1 isles built

New destination location:

GRASS 6.4.0svn (oregonM):/usr4/grassbase > v.proj in=coverage loc=Oregon
map=PERMANENT out=coverage --o

Number of areas: 1
Number of isles: 1

GRASS 6.4.0svn (oregonM):/usr4/grassbase > g.region vect=coverage -p
projection: 99 (Lambert Conformal Conic)
zone:       0
datum:      nad83
ellipsoid:  grs80
north:      562245.46557333
south:      374923.21748089
west:       43789.75941523
east:       290468.06913619
nsres:      32.82324305
ewres:      32.82479171
rows:       5707
cols:       7515
cells:      42888105

> Now you should be ready to do r.proj to re-project the original (feet)
> raster to the new region with units=meters.

   That's the theory, but it just is not working. No matter how many times I
re-do the process the resolutions remain in International Feet, not meters.
ARRGGGHHH-H-H!

   The immediate question is: why is this not working? It's probably user
error, but I am not seeing what I'm doing incorrectly.

Toda,

Rich