I used r.proj for a raster map. Part of the output includes these:
Input:
Cols: 46742 (original: 46742)
Rows: 18655 (original: 18655)
North: 5058076.050000 (original: 5058076.050000)
South: 5039421.050000 (original: 5039421.050000)
West: 537042.640000 (original: 537042.640000)
East: 583784.640000 (original: 583784.640000)
EW-res: 1.000000
NS-res: 1.000000
Output:
Cols: 47260 (original: 168219)
Rows: 20128 (original: 117314)
North: 40028.627671 (original: 121849.848538)
South: 19900.573338 (original: 4535.531862)
West: 341691.218417 (original: 220732.502481)
East: 388951.107430 (original: 388951.107430)
EW-res: 0.999998
NS-res: 1.000003
I assume that the number of cells
and the corners of the region are those for the source and target
projections, but to what do the 'original' values refer? Is this
information on a manual page?
TIA,
Rich
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