<html style="direction: ltr;">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
<style id="bidiui-paragraph-margins" type="text/css">body p { margin-bottom: 0cm; margin-top: 0pt; } </style>
</head>
<body bidimailui-charset-is-forced="true" style="direction: ltr;">
<p><br>
</p>
<div class="moz-cite-prefix">On 01/05/2023 19:52, Rich Shepard
wrote:<br>
</div>
<blockquote type="cite"
cite="mid:b155fd40-939-9531-a0cc-3421814e90e8@appl-ecosys.com">I've
started a new thread on this issue: it's more specific.
<br>
<br>
The script that imports the LiDAR topographic map works:
<br>
<br>
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
<br>
</blockquote>
<p><br>
</p>
<p>First comment: I think you should not have GRASS map names
starting with a digit. <br>
</p>
<p>Can you change the above to <br>
</p>
<p>output=topo_2010</p>
<p>??</p>
<p><br>
</p>
<blockquote type="cite"
cite="mid:b155fd40-939-9531-a0cc-3421814e90e8@appl-ecosys.com">Importing
raster map <2010_topo>...
<br>
100%
<br>
GRASS Dixie_Mountain_2010_LiDAR/PERMANENT:~ > g.list
type=raster
<br>
2010_topo
<br>
<br>
The imported 2010_topo exists in the designated location and
mapset.
<br>
<br>
But, the script to re-project this map fails. I think it has
something to do
<br>
with how I'm mis-using Micha's suggestion of assigning the
re-projected
<br>
output to an environmental variable, NSR.
<br>
<br>
grass /data1/grassdata/Oregon/PERMANENT
<br>
<br>
g.mapset -c NSR_2010_DEM
<br>
<br>
GRASS Oregon/NSR_2010_DEM:~ > NSR=`r.proj -g
loc=Dixie_Mountain_2010_LiDAR map=PERMANENT in=2010_topo
method=lanczos mem=2000 --o`
<br>
Input map <2010_topo@PERMANENT> in location
<Dixie_Mountain_2010_LiDAR>:
<br>
Selected PROJ pipeline:
<br>
+proj=pipeline +step +inv +proj=utm +zone=10 +ellps=GRS80 +step
+proj=lcc
<br>
+lat_0=43.6666666666667 +lon_0=-120.5 +lat_1=46
+lat_2=44.3333333333333
<br>
+x_0=2500000 +y_0=0 +ellps=GRS80
<br>
************************
<br>
GRASS Oregon/NSR_2010_DEM:~ > g.list type=rast
<br>
GRASS Oregon/NSR_2010_DEM:~ >
<br>
<br>
</blockquote>
<p><br>
</p>
<p>The first run of r.proj (with the -g flag) only gives you the
target region settings. You now need to set the region with those
settings, and then run r.proj again (without -g) do actually do
the work.</p>
<p><font face="monospace"><br>
</font></p>
<p><font face="monospace">g.region -ap $NSR</font></p>
<p><font face="monospace">r.proj loc=Dixie_Mountain_2010_LiDAR
map=PERMANENT in=2010_topo method=lanczos mem=2000
output=topo_2010 --o</font></p>
<p><font face="monospace"><br>
</font></p>
<blockquote type="cite"
cite="mid:b155fd40-939-9531-a0cc-3421814e90e8@appl-ecosys.com">There's
no output and I'm not seeing why that is.
<br>
<br>
TIA,
<br>
<br>
Rich
<br>
_______________________________________________
<br>
grass-user mailing list
<br>
<a class="moz-txt-link-abbreviated" href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a>
<br>
<a class="moz-txt-link-freetext" href="https://lists.osgeo.org/mailman/listinfo/grass-user">https://lists.osgeo.org/mailman/listinfo/grass-user</a>
<br>
</blockquote>
<pre class="moz-signature" cols="72">--
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918</pre>
</body>
</html>