Fwd: [GRASS-user] 3D Streets Data
Matthew Pulis
mpulis at gmail.com
Sun Jan 13 12:06:36 EST 2008
Hi Markus,
Thanks alot for your help :)
v.out.ascii on my streets3d gave me :
GRASS 6.3.0RC1 (nc_spm_06):/usr/include > v.out.ascii streets3d
format=standard | head -50
ORGANIZATION: NC OneMap
DIGIT DATE: Nov 7 6
DIGIT NAME: helena
MAP NAME: Wake County roads and streets (lines map)
MAP DATE: Tue Nov 7 20:04:56 2006
MAP SCALE: 1
OTHER INFO:
ZONE: 0
MAP THRESH: 0.000000
VERTI:
L 4 1
625304.64965949 214113.54427121 110.63827327
625301.19968723 214078.52097425 110.53516127
625299.48976068 214056.17660237 110.46902923
625299.2898505 214050.64903112 110.45167028
1 1
L 4 1
670282.21275792 227601.93222423 77.30562497
670374.57454089 227594.34858999 77.75031661
670445.44075096 227585.49027471 78.06294622
670480.4596783 227576.2012396 78.16904676
1 2
L 2 1
640859.16790248 227574.18138204 101.43534006
640807.37009659 227572.59333673 102.82421757
1 3
L 3 1
641027.66472275 226005.86943382 101.61656721
641029.28632052 226018.73890562 101.82055249
641027.87805392 226096.09665666 103.0005327
1 4
L 4 1
650813.12379857 226001.6206011 63.07461956
650817.8816922 226025.06362068 62.85028877
650825.19690683 226054.93404143 62.54284598
650837.89284223 226126.89484171 61.86218497
1 5
L 4 1
657988.42877208 239117.08507621 101.59256415
657999.49858403 239106.02572014 101.24871757
658006.91367573 239097.03756468 105.86902948
658014.55349082 239087.3753951 105.64730814
1 6
L 10 1
658852.96011106 240354.2803688 123.12467481
658856.76542719 240347.57596483 123.13784494
658860.93295374 240341.59644987 123.16069965
658865.82536929 240335.97914532 123.1962029
658871.44267384 240331.08672976 123.24659675
658878.14707781 240326.01313098 123.31148144
I followed your instructions (except the nviz one - since I do not have it
installed) and this is the v.ascii on the new streets3d_new level I
recreated using your method :
GRASS 6.3.0RC1 (nc_spm_06):/usr/include > v.out.ascii streets3d_new
format=standard | head -50
ORGANIZATION: NC OneMap
DIGIT DATE: Nov 7 6
DIGIT NAME: helena
MAP NAME: Wake County roads and streets (lines map)
MAP DATE: Tue Nov 7 20:04:56 2006
MAP SCALE: 1
OTHER INFO:
ZONE: 0
MAP THRESH: 0.000000
VERTI:
L 4 1
625304.64965949 214113.54427121 111.15730809
625301.19968723 214078.52097425 111.15730809
625299.48976068 214056.17660237 111.15730809
625299.2898505 214050.64903112 111.15730809
1 1
L 4 1
670282.21275792 227601.93222423 80.19490021
670374.57454089 227594.34858999 80.19490021
670445.44075096 227585.49027471 80.19490021
670480.4596783 227576.2012396 82.07733239
1 2
L 2 1
640859.16790248 227574.18138204 107.09736996
640807.37009659 227572.59333673 107.09736996
1 3
L 3 1
641027.66472275 226005.86943382 102.38978817
641029.28632052 226018.73890562 103.15115211
641027.87805392 226096.09665666 107.85202063
1 4
L 4 1
650813.12379857 226001.6206011 67.55276255
650817.8816922 226025.06362068 67.11642843
650825.19690683 226054.93404143 65.72167216
650837.89284223 226126.89484171 65.72167216
1 5
L 4 1
657988.42877208 239117.08507621 110.0050774
657999.49858403 239106.02572014 109.75625556
658006.91367573 239097.03756468 111.837334
658014.55349082 239087.3753951 111.66560996
1 6
L 10 1
658852.96011106 240354.2803688 115.85651881
658856.76542719 240347.57596483 115.85651881
658860.93295374 240341.59644987 115.85651881
658865.82536929 240335.97914532 115.85651881
658871.44267384 240331.08672976 115.85651881
658878.14707781 240326.01313098 115.85651881
As you can see the Z coordinate did
change. However, the problem is that when I exported this layer to a
shape file to load it from QGis, and when I exported this to a PostGIS
database there is no showings of the Altitude layer.
Technically I do not need to see the altitude however it is important that
altitude can be read since I need to make routing using the Altitude
constraint.
Any help will be extremely appreciated
Thanks for your help already :)
On 1/11/08, Markus Neteler <neteler at osgeo.org> wrote:
>
> On Jan 5, 2008 6:32 PM, Matthew Pulis <mpulis at gmail.com> wrote:
> >
> >
> >
> >
> > 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
> > |
>
> As you see - the data are 3D!
>
>
> I run
> g.region vect=streets_wake res=100 -pa
> v.drape in=streets_wake rast=elev_state_500m method=bilinear out=streets3d
> #
> r.mapcalc sealevel=0
> nviz el=sealevel vect=streets3d
>
> Looks pretty much 3D (sorry for flat terrain there...).
>
> Check also
> v.out.ascii streets3d format=standard | head -50
> ORGANIZATION: NC OneMap
> DIGIT DATE: Nov 7 6
> DIGIT NAME: helena
> MAP NAME: Wake County roads and streets (lines map)
> MAP DATE: Tue Nov 7 20:04:56 2006
> MAP SCALE: 1
> OTHER INFO:
> ZONE: 0
> MAP THRESH: 0.000000
> VERTI:
> L 4 1
> 625304.64965949 214113.54427121 111.15730809
> 625301.19968723 214078.52097425 111.15730809
> 625299.48976068 214056.17660237 111.15730809
> 625299.2898505 214050.64903112 111.15730809
> 1 1
> L 4 1
> 670282.21275792 227601.93222423 80.19490021
> 670374.57454089 227594.34858999 80.19490021
> 670445.44075096 227585.49027471 80.19490021
> 670480.4596783 227576.2012396 82.07733239
> 1 2
> L 2 1
> 640859.16790248 227574.18138204 107.09736996
> 640807.37009659 227572.59333673 107.09736996
> 1 3
> ...
>
> Markus
>
--
Matthew Pulis
URL : http://www.solutions-lab.net
MSN : pulis_matthew[@]hotmail.com
ICQ : 145951110
Skype : solutions-lab.net
--
Matthew Pulis
URL : http://www.solutions-lab.net
MSN : pulis_matthew[@]hotmail.com
ICQ : 145951110
Skype : solutions-lab.net
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/grass-user/attachments/20080113/e2a262fa/attachment.html
More information about the grass-user
mailing list