[GRASS-dev] [GRASS GIS] #818: v.proj: split_q.c:69: RTreeClassify: Assertion `!p->taken[i]' failed.

[GRASS GIS]

#818: v.proj: split_q.c:69: RTreeClassify: Assertion `!p->taken[i]' failed.
----------------------+-----------------------------------------------------
 Reporter:  pertusus  |       Owner:  grass-dev at lists.osgeo.org
     Type:  defect    |      Status:  new                      
 Priority:  major     |   Milestone:  6.4.0                    
Component:  Vector    |     Version:  6.4.0 RCs                
 Keywords:            |    Platform:  Linux                    
      Cpu:  x86-32    |  
----------------------+-----------------------------------------------------
 On a location with a longitude latitude projection I do a map that
 corresponds with the grid, like:
 {{{
 rows_nr=`g.region -g |grep '^rows=' |sed 's/^rows=//'`
 r.mapcalc "grid_map=col() + ((row()-1) * $rows_nr)"
 r.to.vect input=grid_map out=grid_map_v feature=area
 }}}

 {{{
 g.region -p
 projection: 3 (Latitude-Longitude)
 zone:       0
 datum:      ** unknown (default: WGS84) **
 ellipsoid:  wgs84
 north:      90N
 south:      90S
 west:       180W
 east:       180E
 nsres:      0:30
 ewres:      0:30
 rows:       360
 cols:       720
 cells:      259200
 }}}

 In another location, where the region is:
 {{{
 projection: 99 (Lambert Azimuthal Equal Area)
 zone:       0
 datum:      ** unknown (default: WGS84) **
 ellipsoid:  sphere
 north:      4149500
 south:      -5044500
 west:       -4368500
 east:       4367500
 nsres:      1000
 ewres:      1000
 rows:       9194
 cols:       8736
 cells:      80318784
 }}}

 I try to project the vector and I get an assertion:
 {{{
 v.proj input="grid_map_v" location="cru_maps"
 Building topology for vector map <grid_map_v>...
 Registering primitives...
 v.proj: split_q.c:69: RTreeClassify: Assertion `!p->taken[i]' failed.
 Aborted
 }}}

 In the svn version the file and the assertion is at
 lib/vector/rtree/split.c and the bug may be fixed. I can try to reproduce
 on the svn version if you give me some directions on how to test in
 source.

-- 
Ticket URL: <http://trac.osgeo.org/grass/ticket/818>
GRASS GIS <http://grass.osgeo.org>