[GRASS-user] Region doesn't match EPSG

[Rich Shepard]

On Sun, 30 Apr 2023, Micha Silver wrote:

> What you do is run the r.proj module with the '-g' flag. That gives you
> the region settings in the current mapset that would match the region of
> the input raster from the other mapset. Then you use those setting as
> input to g.region. So:
>
> new_region=`r.proj -g location=<original_location> mapset=<original_mapset> 
> input=<original_raster>`
>
> g.region -ap $new_region

Micha,

(Reformatted to fit the scrsen)

It's still not working:
GRASS NorthsideRock/PERMANENT:~ > NSR=`r.proj -g
loc=Dixie_Mountain_2010_LiDAR map=Northside_rock_2010_topo in=2010_DEM
out=2010_topo method=lanczos_f mem=2000 --o`

Input map <2010_DEM at Northside_rock_2010_topo> in location
<Dixie_Mountain_2010_LiDAR>:

Selected PROJ pipeline:
+proj=pipeline +step +inv +proj=utm +zone=10 +ellps=GRS80 +step +proj=lcc
+lat_0=43.6666666666667 +lon_0=-120.5 +lat_1=46 +lat_2=44.3333333333333
+x_0=2500000 +y_0=0 +ellps=GRS80
************************
GRASS NorthsideRock/PERMANENT:~ > g.region -ap NSR
ERROR: Region <NSR> not found
GRASS NorthsideRock/PERMANENT:~ > g.list type=raster
GRASS NorthsideRock/PERMANENT:~ >

Regards,

Rich