<html style="direction: ltr;">
<head>
<meta http-equiv="Content-Type" content="text/html;
charset=windows-1252">
<style type="text/css">body p { margin-bottom: 0cm; margin-top: 0pt; } </style>
</head>
<body bidimailui-detected-decoding-type="latin-charset"
style="direction: ltr;" text="#000000" bgcolor="#FFFFFF">
If I understand, you have DEM tiles in one projection (location) and
you're trying to use r.proj to transform the tiles into a different
location.<br>
If that's the case, typically you would use the r.proj flag '-g' to
get the region settings in the target location from the raster
extent in it's source location. So you might do something like:<br>
<tt> g.region -p `r.proj -g loc=<source location>
map=<source mapset> input=<source dem>`</tt><br>
<br>
then run r.proj without the -g flag and you should get the
reprojected raster.<br>
<br>
Finally, once you have all the tiles reprojected, use<br>
<tt>g.list raster pattern=<....> sep=comma</tt> <br>
to get a list of all the tiles, then <tt>g.region -p rast=<list
of tiles></tt><br>
to set the computational region to cover all tiles.<br>
<br>
One additional point: consider to patch all DEM tiles together in
their native projection *first* then project the single patched DEM
once to the target projection.<br>
<br>
Regards,<br>
Micha <br>
<br>
<div class="moz-cite-prefix">On 07/06/2018 08:25 PM, Rich Shepard
wrote:<br>
</div>
<blockquote type="cite"
cite="mid:alpine.LNX.2.20.1807061016050.3631@salmo.appl-ecosys.com">
I am trying to add additional DEM maps to the project location.
There is
<br>
one DEM overlaid on a soil mapping unit map for the county in the
attached
<br>
screenshot. My reading of the g.region manual page did not let me
find the
<br>
answers I need.
<br>
<br>
The computational region is set to the soil map. Attempting to
re-project
<br>
the DEM to the immediate west of the existing one fails;
grass7.5svn tells
<br>
me, "ERROR: Input raster map is outside current region." I don't
understand
<br>
why this is since the soil map is much larger than the DEMs and I
used
<br>
'g.region vect=soils' to set the size.
<br>
<br>
Another question is how I can expand the region to add a DEM
immediately
<br>
below the one on the attached map. The eastern edge of the DEM is
beyond the
<br>
limit of the soil map.
<br>
<br>
Learning more about manipulating regions while assembling
project maps
<br>
will help me now and in the future. Pointers to more resources for
me to
<br>
read will be very helpful.
<br>
<br>
Rich
<br>
<br>
<br>
<!--'"--><br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
grass-user mailing list
<a class="moz-txt-link-abbreviated" href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a>
<a class="moz-txt-link-freetext" href="https://lists.osgeo.org/mailman/listinfo/grass-user">https://lists.osgeo.org/mailman/listinfo/grass-user</a></pre>
</blockquote>
<br>
<pre class="moz-signature" cols="72">--
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918</pre>
</body>
</html>