v.generalize: hermite, boyle corrupt data

[msieczka]

1. import the attached shapefile

2. run v.generalize method=hermite on it

The resulting vector map renders fine in GRASS, but makes QGIS crash and say:

qgsgrass.cpp: 374: (error_routine) error_routine (fatal = 1): BUG (Vect_build_sidx_from_topo): node does not exist

After re-building topology the map renders fine in QGIS too. Topology dump of the 2 maps differ like this:

Before topology re-build:

---------- TOPOLOGY DUMP ----------
N,S,E,W,T,B: 36.999861, 35.889244, -6.178626, -7.133487, 0.000000, 0.000000
Nodes (2 nodes, alive + dead ):
node = 1, n_lines = 2, xy = -7.000139, 36.999861
  line =  -1, type = 4, angle = -1.781112
  line =   1, type = 4, angle = -0.417989
node = 2, n_lines = 1, xy = -6.662233, 36.459010
  line =   2, type = 8, angle = -9.000000
Lines (2 lines, alive + dead ):
line = 1, type = 4, offset = 14 n1 = 1, n2 = 1, left/area = -1, right = 1
N,S,E,W,T,B: 36.999861, 35.889244, -6.178626, -7.133487, 0.000000, 0.000000
line = 2, type = 8, offset = 316 n1 = 2, n2 = 2, left/area = 1, right = 0
N,S,E,W,T,B: 36.459010, 36.459010, -6.662233, -6.662233, 0.000000, 0.000000
Areas (1 areas, alive + dead ):
area = 1, n_lines = 1, n_isles = 0 centroid = 2
N,S,E,W,T,B: 36.999861, 35.889244, -6.178626, -7.133487, 0.000000, 0.000000
  line =   1
Islands (1 islands, alive + dead ):
isle = 1, n_lines = 1 area = 0
N,S,E,W,T,B: 36.999861, 35.889244, -6.178626, -7.133487, 0.000000, 0.000000
  line =  -1

After topology re-build:

---------- TOPOLOGY DUMP ----------
N,S,E,W,T,B: 36.999861, 35.889244, -6.178626, -7.133487, 0.000000, 0.000000
Nodes (7 nodes, alive + dead ):
node = 1, n_lines = 2, xy = -7.000139, 36.999861
  line =  -1, type = 4, angle = -1.781112
  line =   1, type = 4, angle = -0.417989
node = 7, n_lines = 1, xy = -6.662233, 36.459010
  line =   7, type = 8, angle = -9.000000
Lines (7 lines, alive + dead ):
line = 1, type = 4, offset = 14 n1 = 1, n2 = 1, left/area = -1, right = 1
N,S,E,W,T,B: 36.999861, 35.889244, -6.178626, -7.133487, 0.000000, 0.000000
line = 7, type = 8, offset = 316 n1 = 7, n2 = 7, left/area = 1, right = 0
N,S,E,W,T,B: 36.459010, 36.459010, -6.662233, -6.662233, 0.000000, 0.000000
Areas (1 areas, alive + dead ):
area = 1, n_lines = 1, n_isles = 0 centroid = 7
N,S,E,W,T,B: 36.999861, 35.889244, -6.178626, -7.133487, 0.000000, 0.000000
  line =   1
Islands (1 islands, alive + dead ):
isle = 1, n_lines = 1 area = 0
N,S,E,W,T,B: 36.999861, 35.889244, -6.178626, -7.133487, 0.000000, 0.000000
  line =  -1

[msieczka]

Same happens for:

v.generalize input=w007n36e output=w007n36e_boyle method=boyle look_ahead=2 iterations=100

[mmetz]

Fixed in 6.4.2. Close ticket?

Markus M

[neteler]

Closing since fixed.