[GRASS-user] Reprojection regions

[Rich Shepard]

On Mon, 27 Sep 2021, Rich Shepard wrote:

> I used r.proj for a raster map. Part of the output includes these:

I have another question about the reprojection of these topographic maps.

The input location has this projection and units:
name: NAD_1983_HARN_StatePlane_Washington_South_FIPS_4602_Feet
ellps: grs80
proj: lcc
lat_0: 45.3333333333333
lon_0: -120.5
lat_1: 45.8333333333333
lat_2: 47.3333333333333
x_0: 500000
y_0: 0
no_defs: defined
unit: U.S. Surveyor's Foot
units: U.S. Surveyor's Feet
meters: 0.3048006096012192

The output location has this projection and units:
name: NAD83(HARN) / Washington South
datum: nad83harn
ellps: grs80
proj: lcc
lat_0: 45.3333333333333
lon_0: -120.5
lat_1: 47.3333333333333
lat_2: 45.8333333333333
x_0: 500000
y_0: 0
no_defs: defined
unit: meter
units: meters
meters: 1

So why does r.proj show this pipeline (note the in and out units at the
end):
Selected PROJ pipeline:
+proj=pipeline +step +inv +proj=lcc +lat_0=45.3333333333333 +lon_0=-120.5
+lat_1=47.3333333333333 +lat_2=45.8333333333333 +x_0=500000 +y_0=0
+ellps=GRS80 +step +proj=lcc +lat_0=45.3333333333333 +lon_0=-120.5
+lat_1=45.8333333333333 +lat_2=47.3333333333333 +x_0=500000 +y_0=0
+ellps=GRS80 +step +proj=unitconvert +xy_in=m +xy_out=us-ft

Shouldn't the xy_in be us-ft and the xy_out be m?

Rich