[GRASS-user] Re: tin to raster dem: SOLVED

[ivan marchesini]

Dear Thomas, dear all...
my application is, at the moment, really simple...
I have the stratigraphy of tree points very near each other...
I was able to note the same stratigraphic sequence in the tree points,
but at different deeps....
My needs were to rebuild the bedding of the bottom and of the top of
some geologic layers....
in this way I can estimate the deep of these layer in different points
from the tree already sampled...
this is interesting for geothecnical reasons... because the thickness of
the layers is important for the calculation regarding the basement of
tree building
I have used the modified version of v.surf.nnbathy to find the plane
passing through tree 3d point (I had a vector 3d point map and I have
converted it into a simple raster map)
in this way I have obtained a raster map representing a 3d triangular
plane (see fig1.jpg)
using r.slope.aspect I have been able to obtain deep and direction of
the plane
then using r.plane I have been able to create a plane passing through
the tree points and extending outside of them (see fig2.jpg)
I hope this help to understand!!!
best regards
Ivan



Il giorno dom, 18/06/2006 alle 06.33 -0400, Thomas Adams ha scritto:
> Ivan,
> 
> Do you mind posting an example of your results? Also, out of curiosity, 
> how do you plan on making use of this capability?
> 
> Regards,
> Tom
> 
> 
> ivan marchesini wrote:
> > Dear all, dear Maciek,
> > many thanks to all for the answers..
> > the solution proposed by Maciek works really well...
> > it seems that the module v.surf.nnbathy modified changing in line 138
> > "-P alg=nn" with "-P alg=l" is able to create a plane passing for tree
> > rasterized points!!!!
> > many thanks!!!
> > Ivan
> >
> >
> >
> > Il giorno sab, 17/06/2006 alle 13.00 +0200, Maciek Sieczka ha scritto:
> >   
> >> On Fri, 16 Jun 2006 10:50:23 +0200
> >> ivan marchesini <marchesini at unipg.it> wrote:
> >>
> >>     
> >>> Dear all...
> >>> I have been able to create a tin starting from the v.delaunay module
> >>> and then manual editing the v.delaunay output exported with v.in.ascii
> >>> (format=standard) 
> >>> now I can see, using nviz, the some filled triangles in the 3d
> >>> space...!!!!
> >>> now the problem:
> >>> I would like to transfrom the tin surface in a raster surface...
> >>> is it possible in some way!!!!
> >>> any idea????
> >>>       
> >> This *might* work:
> >>
> >> Get the apices of your triangles into points vector, rasterize it
> >> (v.to.rast) and feed to my r.surf.nnbathy
> >> (http://grass.gdf-hannover.de/wiki/GRASS_AddOns).
> >>
> >> You will only have to modify the script slightly. In line 138 put "-P
> >> alg=l" instead of "-P alg=nn". This will make nnbathy triangulate
> >> only, without any smoothing.
> >>
> >> Let me know if it worked.
> >>
> >> Maciek
> >>
> >> P.S.
> >>
> >> Please don't post to grasslist at baylor.edu. Post to 
> >> grassuser at grass.itc.it.
> >>
> >> ------------------------------------------------------------------------
> >> CIEP?E KRAJE - CIEP?E MORZA. Szukasz atrakcyjnego wypoczynku w przyst?pnej cenie, zapoznaj si? z nasz? ofert?.
> >> ZAPRASZAMY
> >>
> >> www.skarpatravel.pl
> >>
> >>
> >>     
> 
> 
-- 
Ivan Marchesini
Department of Civil and Environmental Engineering
University of Perugia
Via G. Duranti 93/a 
06125
Perugia (Italy)
e-mail: marchesini at unipg.it
        ivan.marchesini at gmail.com
tel: +39(0)755853760
fax: +39(0)755853756
jabber: geoivan73 at jabber.org

-------------- next part --------------
A non-text attachment was scrubbed...
Name: fig1.jpg
Type: image/jpeg
Size: 23976 bytes
Desc: not available
Url : http://lists.osgeo.org/pipermail/grass-user/attachments/20060618/3d1ea0b1/fig1.jpg
-------------- next part --------------
A non-text attachment was scrubbed...
Name: fig2.jpg
Type: image/jpeg
Size: 18691 bytes
Desc: not available
Url : http://lists.osgeo.org/pipermail/grass-user/attachments/20060618/3d1ea0b1/fig2.jpg