Having created a new location using EPSG:2856 I tried to reproject a mapset that defines the region I want to use for this project: r.proj loc=<old_location_name> map=bathymetry in=local_depths
The result was not as I expected: Found 3 possible transformations ************************ Operation 1: Description: Inverse of SPCS83 Washington South zone (meters) + Inverse of NAD83 to NAD83(HARN) (28) + UTM zone 10N Area of use: United States (USA) - Oregon and Washington. Accuracy within area of use: 0.05 m PROJ string: +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 +inv +proj=hgridshift +grids=us_noaa_WO.tif +step +proj=utm +zone=10 +ellps=GRS80 ************************ This is the most appropriate as the other two are used in the whole world and states other than Washington and Oregon. But, I don't know how to use this information now that GRASS gave it to me. The second problem is: ERROR: Input raster map is outside current region
g.region -p
projection: 99 (NAD83(HARN) / Washington South) zone: 0 datum: nad83harn ellipsoid: grs80 north: 1 south: 0 west: 0 east: 1 nsres: 1 ewres: 1 rows: 1 cols: 1 cells: 1 I need help to get out of the weeds in which I floundered. Apparently creating the location with the EPSG code of 2856 was insufficient to define all CRS parameters and the region. TIA, Rich _______________________________________________ grass-user mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/grass-user
