[GRASS-dev] patch for v.voronoi/v.delaunay

On 4/22/08 4:40 PM, "grass-dev-request@lists.osgeo.org"
<grass-dev-request@lists.osgeo.org> wrote:

Date: Wed, 23 Apr 2008 01:40:11 +0200
From: "Markus Neteler" <neteler@osgeo.org>
Subject: Re: [GRASS-dev] patch for v.voronoi/v.delaunay
To: "GRASS developers list" <grass-dev@lists.osgeo.org>
Message-ID:
<86782b610804221640s3c03cab0w5533af1110694349@mail.gmail.com>
Content-Type: text/plain; charset=ISO-8859-1

On Wed, Apr 23, 2008 at 1:07 AM, Markus Neteler <neteler@osgeo.org> wrote:

On Fri, Feb 1, 2008 at 2:16 PM, Benjamin Ducke wrote:

...

The output vector map will now also preserve the 3D coordinates of
the original input points. However, the centroids for the triangles
are all wrong (placed at z=0.0). I don't see a quick way of fixing
this right now.

I have fixed it using Vect_tin_get_z() for delaunay (someone needs to
do it for voronoi).

I have submitted the fixed code to SVN:
Changeset 31079 – GRASS GIS

For me it now generates a 3D Delaunay TIN.

Is there anyway to display a TIN in 3D with GRASS? In NVIZ, it just looks
like a Delaunay triangulation network unless there is something I'm not
doing.

Michael

I wonder if I could colorize the 3D TIN somehow using v.colors.
So far it does not have a DB connected.

__________________________________________
Michael Barton, Professor of Anthropology
Director of Graduate Studies
School of Human Evolution & Social Change
Center for Social Dynamics & Complexity
Arizona State University

phone: 480-965-6213
fax: 480-965-7671
www: http://www.public.asu.edu/~cmbarton

Michael Barton wrote:

Is there anyway to display a TIN in 3D with GRASS? In NVIZ, it just
looks like a Delaunay triangulation network unless there is something
I'm not doing.

Make sure your 3D triangulation output "Map is 3D".

AFAIK faces will always be based on a single color (some todo).

off on a tangent:
As I've ranted on about before, TINs are really just the point data that
define them, so in general you'd be much better off extracting the points
and using a real interpolation module like v.surf.rst to make a surface
from those points. The exception is from e.g. the poster the other day
who's point data was not uniformly distributed, density increased by
orders of magnitude as you get nearer to an object of interest (dynamic
computational mesh size)
e.g. http://gfs.sourceforge.net/examples/examples/tangaroa.html

v.surf.rst does the same sort of thing during processing for no-data
search areas, but the final product there is at the region resolution.
The problem the poster the other day had was that the overall grid at a
resolution appropriate for the coastal areas would make a massive array
when extended to include the oceanic data, and you probably have to make
a number of maps each as a different scale.

I remember in the past when talking about the GRASS 7 raster format it
was mentioned that perhaps Adaptive JPEG-like blocks could be used to
save the data instead of row-by-row. This is slightly a perversion of
that, but an interesting problem to think about.

Hamish

      ____________________________________________________________________________________
Be a better friend, newshound, and
know-it-all with Yahoo! Mobile. Try it now. http://mobile.yahoo.com/;_ylt=Ahu06i62sR8HDtDypao8Wcj9tAcJ

On Wed, Apr 23, 2008 at 6:26 AM, Michael Barton <michael.barton@asu.edu> wrote:

From: "Markus Neteler" <neteler@osgeo.org>

...

>> I have submitted the fixed code to SVN:
>> http://trac.osgeo.org/grass/changeset/31079
>>
>> For me it now generates a 3D Delaunay TIN.

Is there anyway to display a TIN in 3D with GRASS? In NVIZ, it just looks
like a Delaunay triangulation network unless there is something I'm not
doing.

Attached a v.delaunay trangulation from 3D Lidar points
as visualized in NVIZ.

Markus

(attachments)

vaihingen3d_tin.jpg

Markus:

Attached a v.delaunay trangulation from 3D Lidar points
as visualized in NVIZ.

It would be interesting to create a nviz screenshot comparing that with a
lidar -> v.in.ascii -> v.surf.rst surface placed below instead of the
plane.

Hamish

      ____________________________________________________________________________________
Be a better friend, newshound, and
know-it-all with Yahoo! Mobile. Try it now. http://mobile.yahoo.com/;_ylt=Ahu06i62sR8HDtDypao8Wcj9tAcJ

On Wed, Apr 23, 2008 at 1:12 PM, Hamish <hamish_b@yahoo.com> wrote:

Markus:

> Attached a v.delaunay trangulation from 3D Lidar points
> as visualized in NVIZ.

It would be interesting to create a nviz screenshot comparing that with a
lidar -> v.in.ascii -> v.surf.rst surface placed below instead of the
plane.

Here are the Lidar data:

North = 5403850 m; South = 5403547 m
West = 512700 m; Est = 512835 m
UTM Zone 32U, EPSG code 32632 (Germany, Stuttgart, Vaihingen)
Available data: last pulses

Website: http://www.commission3.isprs.org/wg3/
-> "Data sets" -> "downloaded"

db.in.ogr table=lidar1.dbf output=lidar1table key=id
v.in.db table=lidar1table x=N1 y=N2 z=N3 output=lidar1map3d key=id

We could use it as benchmark dataset (to have one time a data set from
Europe!).

Markus

I compiled GRASS a week ago. Is this new? Maybe I'm doing something wrong
with v.delaunay.

Michael

On 4/23/08 3:29 AM, "Markus Neteler" <neteler@osgeo.org> wrote:

On Wed, Apr 23, 2008 at 6:26 AM, Michael Barton <michael.barton@asu.edu>
wrote:

From: "Markus Neteler" <neteler@osgeo.org>

...

I have submitted the fixed code to SVN:
Changeset 31079 – GRASS GIS

For me it now generates a 3D Delaunay TIN.

Is there anyway to display a TIN in 3D with GRASS? In NVIZ, it just looks
like a Delaunay triangulation network unless there is something I'm not
doing.

Attached a v.delaunay trangulation from 3D Lidar points
as visualized in NVIZ.

Markus

__________________________________________
Michael Barton, Professor of Anthropology
Director of Graduate Studies
School of Human Evolution & Social Change
Center for Social Dynamics & Complexity
Arizona State University

phone: 480-965-6213
fax: 480-965-7671
www: http://www.public.asu.edu/~cmbarton

On Wed, Apr 23, 2008 at 4:57 PM, Michael Barton <michael.barton@asu.edu> wrote:

I compiled GRASS a week ago. Is this new? Maybe I'm doing something wrong
with v.delaunay.

I added it now 20h ago:
http://trac.osgeo.org/grass/changeset/31079

So you need to update from SVN.

Markus