[GRASS-user] Re-projection fails

Rich Shepard rshepard at appl-ecosys.com
Mon May 1 09:52:21 PDT 2023


I've started a new thread on this issue: it's more specific.

The script that imports the LiDAR topographic map works:

GRASS Dixie_Mountain_2010_LiDAR/PERMANENT:~ > r.in.gdal in=$HOME/projects/oregon/northside-rock/data/LDQ-45122F8/2010-USACE-Columbia-River/Bare_Earth/be45122f8/hdr.adf out=2010_topo mem=2000 --o
Importing raster map <2010_topo>...
  100%
GRASS Dixie_Mountain_2010_LiDAR/PERMANENT:~ > g.list type=raster
2010_topo

The imported 2010_topo exists in the designated location and mapset.

But, the script to re-project this map fails. I think it has something to do
with how I'm mis-using Micha's suggestion of assigning the re-projected
output to an environmental variable, NSR.

grass /data1/grassdata/Oregon/PERMANENT

g.mapset -c NSR_2010_DEM

GRASS Oregon/NSR_2010_DEM:~ > NSR=`r.proj -g loc=Dixie_Mountain_2010_LiDAR map=PERMANENT in=2010_topo method=lanczos mem=2000 --o`
Input map <2010_topo at PERMANENT> 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 Oregon/NSR_2010_DEM:~ > g.list type=rast
GRASS Oregon/NSR_2010_DEM:~ >

There's no output and I'm not seeing why that is.

TIA,

Rich


More information about the grass-user mailing list