[GRASS-user] RE: 3D Streets Data [v.drape]

[Dylan Beaudette]

On Thursday 10 January 2008, Matthew Pulis wrote:
> Hi again,
>
>
>
> I am really sorry I am sending again but I have tried to read many
> documents on how to use v.drape all gave me the same conclusions I did
> here.
>
>
>
> Is there some extra information you might require that I make available in
> order to get some knowledge on this situation please?
>
>
>
> Also anyone can suggest how I can work out my problem without using v.drape
> please?
>
>
>
> Thanks again J
>
>
>
> Matthew
>


This is odd behavior indeed. The usage of v.drape is relatively simple, and 
the only problem that one might encounter is when the elevation surface does 
not extend beyond the extents of the input 2D vector.

Some questions:

1. are you sure that you are working with a vector file which contains lines 
(not boundaries) ? If this is the case v.drape will exit with an error- but 
will leave an unfinished output vector (I need to fix this!)

2. can you give some examples on how you are invoking the command?



I just gave it a try in the Spearfish dataset

# digit new line:
v.in.ascii --o format=standard out=newline << EOF
VERTI:
L  12 1
 594003.22429907 4917469.71962617
 594772.28971963 4919333.22429907
 596280.8411215 4920102.28971963
 598292.24299065 4920545.98130841
 599652.89719626 4920338.92523364
 600214.90654206 4919717.75700935
 601900.93457944 4918978.27102804
 602581.26168224 4919214.90654206
 603439.06542056 4920753.03738318
 603379.90654206 4921581.26168224
 603734.85981308 4922675.70093458
 604710.98130841 4922646.12149533
 1     1
EOF

# convert to 3D
v.drape in=newline out=newline3d rast=elevation.dem

# check:
v.info # looks ok!

# dump for testing:
VERTI:
L  12 1
 594003.22429907 4917469.71962617 1454        
 594772.28971963 4919333.22429907 1580        
 596280.8411215 4920102.28971963 1490        
 598292.24299065 4920545.98130841 1373        
 599652.89719626 4920338.92523364 1439        
 600214.90654206 4919717.75700935 1457        
 601900.93457944 4918978.27102804 1328        
 602581.26168224 4919214.90654206 1345        
 603439.06542056 4920753.03738318 1266        
 603379.90654206 4921581.26168224 1221        
 603734.85981308 4922675.70093458 1182        
 604710.98130841 4922646.12149533 1198        
 1     1 



The last column is the elevation.

I will look over the source code, and will try and get the bad behavior noted 
above fixed. Any devs have some hints on cleaning up if a module fails to 
complete?

Cheers,

Dylan





>
> From: Matthew Pulis [mailto:mpulis at gmail.com]
> Sent: Saturday, January 05, 2008 6:33 PM
> To: 'grass-user at lists.osgeo.org'
> Subject: 3D Streets Data
>
>
>
> Hi,
>
>
>
> First of all Happy New Year to all on the list J
>
>
>
> I am trying to build a 3D Street Network, thus where the streets have
> elevation so I can do routing depending on slope level.
>
>
>
> I am using the GRASS Book Dataset of North Carolina.
>
>
>
> I tried to combine : streets_wake and elev_state_500m layers together,
> using the v.drape :
>
> v.drape in=streets_wake rast=elev_state_500m method=bilinear out=streets3d
>
> But to no avail. This is the v.info on streets3d :
> |   Number of lines:        49746           Number of islands:    0
> |
>  |   Number of boundaries:   0               Number of faces:      0
>  |
>  |
>  |   Number of centroids:    0               Number of kernels:    0
>  |
>  |
>  |
>  |
>  |
>  |   Map is 3D:              1
>  |
>  |
>  |   Number of dblinks:      1
>  |
>  |
>  |
>  |
>  |
>  |         Projection: Lambert Conformal Conic
>  |
>  |
>  |               N:   258102.57210146    S:   196653.34855345
>  |
>  |
>  |               E:   676816.24830376    W:    610869.3848495
>  |
>  |
>  |               B:       49.16141063    T:      156.87280191
>
> The result of streets3d ended up exactly as streets_wake without any
> altitude data.
>
>
>
> Any idea how I can do this endeavour please? I am planning the dataset for
> routing using pGRouting afterwards.
>
>
>
> Thanks and regards
>
>
>
> Matthew
>
>
>
>   _____
>
> I am using the free version of SPAMfighter for private users.
> It has removed 19049 spam emails to date.
> Paying users do not have this message in their emails.
> Try SPAMfighter <http://www.spamfighter.com/len>  for free now!
>
>
>
>   _____
>
> I am using the free version of SPAMfighter for private users.
> It has removed 19198 spam emails to date.
> Paying users do not have this message in their emails.
> Try SPAMfighter <http://www.spamfighter.com/len>  for free now!



-- 
Dylan Beaudette
Soil Resource Laboratory
http://casoilresource.lawr.ucdavis.edu/
University of California at Davis
530.754.7341