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 _______________________________________________ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user