<html style="direction: ltr;">
<head>
<meta http-equiv="Content-Type" content="text/html;
charset=windows-1252">
<style id="bidiui-paragraph-margins" type="text/css">body p { margin-bottom: 0cm; margin-top: 0pt; } </style>
</head>
<body bidimailui-detected-decoding-type="latin-charset"
style="direction: ltr;">
<p><br>
</p>
<div class="moz-cite-prefix">On 9/30/21 7:33 PM, Rich Shepard wrote:<br>
</div>
<blockquote type="cite"
cite="mid:f437428-858a-40a2-afd8-fdd8f55be261@appl-ecosys.com">Rather
than use raster DTM maps I decided to use depth measurement points
<br>
themselves.
<br>
<br>
The data look like this:
<br>
lon,lat,depth
<br>
7709581.42,697951.64,24
<br>
7705545.87,698624.07,16
<br>
7709732.27,699488.73,27
<br>
7712090.79,699428.24,62
<br>
7716863.99,701238.89,32
<br>
7701540.85,698790.50,31
<br>
7713676.33,699192.82,58
<br>
7705283.41,697549.00,31
<br>
7705037.10,698489.22,16
<br>
<br>
The command I'm using is:
<br>
d.vect map=lady_channel_sounding_points display=zcoor type=point
zcolor=blue size=7
<br>
<br>
and the resulting map is attached.
<br>
</blockquote>
<p><br>
</p>
<p>Here's what I did. Notice that I added the "-z=3" option to
v.in.ascii so that the new points vector is 3D. Then the "zcolor="
option to d.vect works.</p>
<p><br>
</p>
<p># Create new Location/Mapset based on the georeference tiff file<br>
</p>
<p>grass -c ~/Downloads/columbia_2010_e_dtm_35.tif --text
./columbia/PERMANENT<br>
<br>
#Points text file<br>
micha@RMS:tmp$ cat << EOF > depth_points.txt<br>
lon,lat,depth<br>
7709581.42,697951.64,24<br>
7705545.87,698624.07,16<br>
7709732.27,699488.73,27<br>
7712090.79,699428.24,62<br>
7716863.99,701238.89,32<br>
7701540.85,698790.50,31<br>
7713676.33,699192.82,58<br>
7705283.41,697549.00,31<br>
7705037.10,698489.22,16<br>
EOF<br>
<br>
# Create GRASS vector from CSV file<br>
# Note the 'z=3" parameter, to creaate a 3D vector with Z coord<br>
cat depth_points.txt | v.in.ascii input=- output=depth_points
separator=comma z=3 skip=1 columns="x DOUBLE, y DOUBLE, depth
DOUBLE" --o<br>
v.info depth_points | grep 3D<br>
| Map is 3D: Yes<br>
g.region -ap vect=depth_points<br>
<br>
# Display<br>
d.mon wx1<br>
d.vect depth_points zcolor=byr size=10 icon=basic/diamond<br>
<br>
# Result attached</p>
<p><br>
</p>
<blockquote type="cite"
cite="mid:f437428-858a-40a2-afd8-fdd8f55be261@appl-ecosys.com">
<br>
Perhaps I should make the computational region smaller to make
individual
<br>
values visible, but are the display and color values what they
should be?
<br>
The values are displayed in red, not blue.
<br>
<br>
TIA,
<br>
<br>
Rich
<br>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<pre class="moz-quote-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>
<pre class="moz-signature" cols="72">--
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918</pre>
</body>
</html>