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