<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html style="direction: ltr;">
<head>
<meta content="text/html;charset=UTF-8" http-equiv="Content-Type">
<title></title>
</head>
<body style="direction: ltr;" bgcolor="#ffffff" text="#000000">
georgew wrote:
<blockquote cite="mid:20803902.post@talk.nabble.com" type="cite">
<pre wrap="">Micha, I followed your instructions with the following results:
----
Import data from 2 shapefiles, Contour (layer 1)and Track (layer 2): (no
problems reported)
v.in.ogr dsn=/home/george/GRASSDATA/ output=testvectmap min_area=0.0001
snap=-1
----
</pre>
</blockquote>
How many shapefiles did you have in the GRASSDATA directory?<br>
What does testvectmap contain? You should make two separate GRASS
vectors, one for the tracks and one for the contours. You can use the
'layer' parameter each time to specify each shapefile, like:<br>
v.in.ogr dsn=/home/george/GRASSDATA/ layer=contours.shp output=contours<br>
<br>
<blockquote cite="mid:20803902.post@talk.nabble.com" type="cite">
<pre wrap="">Display vector map: (no problems, map displayes correctly)
d.vect map=testvectmap@testvector color=0:0:0 lcolor=0:0:0
fcolor=170:170:170 display=shape type=point,line,boundary,area
icon=basic/circle size=5 layer=-1 lsize=8 xref=left yref=center llayer=1
----
</pre>
</blockquote>
At this point I'd check to make sure the elevation data came thru:<br>
v.info -c # shows the column headers, to make sure you have
ELEVATION<br>
v.db.select # shows the actual ELEVATION values <br>
<blockquote cite="mid:20803902.post@talk.nabble.com" type="cite">
<pre wrap="">Convert vector to raster and thin: (no problems reported)
v.to.rast input=testvectmap@testvector output=testrastmap use=attr
type=point,line,area layer=1 column=ELEVATION value=1 rows=100000
r.thin input=testrastmap@testvector output=thinrastmap iterations=200
----
Display raster map (PROBLEM:only displays large, brightly coloured
rectangles, no contours
d.rast map=testrastmap@testvector -o
----
</pre>
</blockquote>
How did you set up the Location/Mapset ? You might have a problem with
the resolution.<br>
What do you get from:<br>
r.info testrastmap<br>
g.region -p<br>
<br>
<br>
<blockquote cite="mid:20803902.post@talk.nabble.com" type="cite">
<pre wrap="">At this point I stopped as I got the impression something was already wrong
but could not work out what.
Also, should I not rasterize the track data as well?
Any help would be appreciated. If you require the source shapefiles I have
no problems sending them to you.
</pre>
</blockquote>
Let's keep stabbing at it from your end. That's always the best way to
understand the methodology. <br>
<blockquote cite="mid:20803902.post@talk.nabble.com" type="cite">
<pre wrap="">Micha Silver wrote:
</pre>
<blockquote type="cite">
<pre wrap="">
Yes.
I think you should be looking at converting the shapefile contour lines
to a raster....
v.in.org to import the contour lines and tracks into grass
v.to.rast and r.thin to convert the contour lines to a raster map.
r.surf.contour to create an elevation surface from the rasterized
contours
r.profile to get the elevations along the tracks.
Post back with specific questions if you hit a snag.
Regards,
Micha
___________________________________________
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="http://lists.osgeo.org/mailman/listinfo/grass-user">http://lists.osgeo.org/mailman/listinfo/grass-user</a>
</pre>
</blockquote>
<pre wrap=""><!---->
</pre>
</blockquote>
<br>
</body>
</html>