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