[GRASS-dev] Building a vector topology (and failing)

[Michael Barton]

I've been having topology problems too and don't quite understand. I'm getting a coor file too large error and v. build is not fixing it the way it used to. This is with current svn GRASS 7

Michael


> Message: 1
> Date: Mon, 12 Mar 2012 09:05:54 -0600
> From: Seth Price <seth at pricepages.org>
> Subject: [GRASS-dev] Building a vector topology (and failing)
> To: grass-dev <grass-dev at lists.osgeo.org>
> Message-ID: <FB3BC20D-073E-46A0-B8B0-696AD54E22B8 at pricepages.org>
> Content-Type: text/plain; charset="us-ascii"
> 
> Hey all, I've been trying to get my topology to build for a while now in my module, and I'm not getting anywhere. I currently have it to the point where I'm simply trying to draw a triangle with corners (0,0), (0,100), & (100,0). It calls Vect_append_point() three times, then Vect_write_line(,GV_BOUNDARY,,) to create the polygon. When building the topology, it fails to create an area with "WARNING: Number of incorrect boundaries: 1".
> 
> I'm used to raster programming, but I would think that this would work. You can see in the debug that the error shouldn't happen. There aren't duplicate angles. I'm using a recent SVN of GRASS 7.0. Anyone have any idea what I'm doing wrong? Below is the relevant debug output.
> Thanks,
> Seth
> 
> [...]
> Registering primitives...
> D3/3: Vect_read_next_line()
> D3/3: V1_read_next_line_nat()
> D3/3: Vect__Read_line_nat: offset = 18
> D3/3:     type = 4, do_cats = 1 dead = 0
> D3/3:     n_cats = 1
> D3/3:     n_points = 3
> D3/3:     off = 107
> D3/3: Register line: offset = 18
> D3/3: dig_spidx_add_line(): line = 1
> D3/3: Register node: type = 4,  589971.509972,5300028.490028
> D3/3: dig_find_node()
> D3/3: node = 0
> D3/3: dig_add_node(): n_nodes = 0, alloc_nodes = 0
> D3/3: dig_spidx_add_node(): node = 1, x,y,z = 589971.509972, 5300028.490028, 0.000000
> D3/3: new node = 1, n_nodes = 1, alloc_nodes = 1000
> D3/3: Add new node: 1
> D3/3: dig_node_add_line(): node = 1 line = 1
> D3/3: dig_node_alloc_line(): add = 1
> D3/3:     angle = -1.570796
> D3/3: dig_node_add_line(): line 1 added position 0 n_lines: 1 angle -1.570796
> D3/3: Register node 591396.011396,5300028.490028
> D3/3: dig_find_node()
> D3/3: node = 0
> D3/3: dig_add_node(): n_nodes = 1, alloc_nodes = 1000
> D3/3: dig_spidx_add_node(): node = 2, x,y,z = 591396.011396, 5300028.490028, 0.000000
> D3/3: new node = 2, n_nodes = 2, alloc_nodes = 1000
> D3/3: Add new node: 2
> D3/3: dig_node_add_line(): node = 2 line = -1
> D3/3: dig_node_alloc_line(): add = 1
> D3/3:     angle = -2.356194
> D3/3: dig_node_add_line(): line -1 added position 0 n_lines: 1 angle -2.356194
> D3/3: dig_cidx_add_cat(): field = 4 cat = 1 line = 1 type = 4
> D3/3: Vect_read_next_line()
> D3/3: V1_read_next_line_nat()
> D3/3: Vect__Read_line_nat: offset = 107
> 1 primitives registered
> 3 vertices registered
> Building areas...
> 100%
> D3/3: Build area for line = 1, side = 1
> D3/3: Vect_build_line_area() line = 1, side = 1
> D3/3: dig_line_get_area(): line = 1, side = 1 (left), area = 0
> D3/3: dig_build_area_with_line(): first_line = 1, side = 1
> D3/3: dig_node_line_angle: node = 1 line = 1
> D3/3: dig__angle_next_line: line = 1, side = 2, type = 4
> D3/3:  node = 1
> D3/3:   n_lines = 1
> D3/3:   i = 0 line = 1 angle = -1.570796
> D3/3:   current position = 0
> D3/3:   next = 0 line = 1 angle = -1.570796
> D3/3:   this one
> D3/3: next_line = 1
> D3/3: dig_node_angle_check: line = 1, type = 4
> D3/3: dig_node_line_angle: node = 1 line = 1
> D3/3: dig__angle_next_line: line = 1, side = 2, type = 4
> D3/3:  node = 1
> D3/3:   n_lines = 1
> D3/3:   i = 0 line = 1 angle = -1.570796
> D3/3:   current position = 0
> D3/3:   next = 0 line = 1 angle = -1.570796
> D3/3:   this one
> D3/3: dig_node_line_angle: node = 1 line = 1
> D3/3:   The line to the right has the same angle: node = 1, line = 1
> D3/3: Cannot build area, a neighbour of the line 1 has the same angle at the node
> D3/3:   n_lines = 0
> D3/3: Build area for line = 1, side = 2
> D3/3: Vect_build_line_area() line = 1, side = 2
> D3/3: dig_line_get_area(): line = 1, side = 2 (right), area = 0
> D3/3: dig_build_area_with_line(): first_line = 1, side = 2
> D3/3: dig_node_line_angle: node = 1 line = 1
> D3/3: dig__angle_next_line: line = -1, side = 2, type = 4
> D3/3:  node = 2
> D3/3:   n_lines = 1
> D3/3:   i = 0 line = -1 angle = -2.356194
> D3/3:   current position = 0
> D3/3:   next = 0 line = -1 angle = -2.356194
> D3/3:   this one
> D3/3: next_line = -1
> D3/3: dig_node_angle_check: line = -1, type = 4
> D3/3: dig_node_line_angle: node = 2 line = -1
> D3/3: dig__angle_next_line: line = -1, side = 2, type = 4
> D3/3:  node = 2
> D3/3:   n_lines = 1
> D3/3:   i = 0 line = -1 angle = -2.356194
> D3/3:   current position = 0
> D3/3:   next = 0 line = -1 angle = -2.356194
> D3/3:   this one
> D3/3: dig_node_line_angle: node = 2 line = -1
> D3/3:   The line to the right has the same angle: node = 2, line = -1
> D3/3: Cannot build area, a neighbour of the line -1 has the same angle at the node
> D3/3:   n_lines = 0
> 0 areas built
> 0 isles built
> [...]
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: http://lists.osgeo.org/pipermail/grass-dev/attachments/20120312/e2fe1efa/attachment-0001.html
> 
> ------------------------------
> 
> _______________________________________________
> grass-dev mailing list
> grass-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-dev
> 
> End of grass-dev Digest, Vol 73, Issue 14
> *****************************************

_____________________
C. Michael Barton
Visiting Scientist, Integrated Science Program
National Center for Atmospheric Research &
University Consortium for Atmospheric Research
303-497-2889 (voice)

Director, Center for Social Dynamics & Complexity 
Professor of Anthropology, School of Human Evolution & Social Change
Arizona State University
www: http://www.public.asu.edu/~cmbarton, http://csdc.asu.edu