[GRASS-dev] [GRASS GIS] #2613: Polygons disappear after v.overlay AND operation

[GRASS GIS]

#2613: Polygons disappear after v.overlay AND operation
--------------------------------------------+-------------------------------
 Reporter:  dido                            |       Owner:  grass-dev@…              
     Type:  defect                          |      Status:  new                      
 Priority:  major                           |   Milestone:                           
Component:  Vector                          |     Version:  7.0.0                    
 Keywords:  v.overlay and polygons damaged  |    Platform:  MSWindows 7              
      Cpu:  x86-64                          |  
--------------------------------------------+-------------------------------
 A shapefile polygon set A has been cut using a grid polygon B using a
 v.overlay AND operation. The result is expected to contain all the A
 polygons that are enclosed within the grid polygon B.

 The output is almost as expected but several A polygons within the B
 rectangle are missing. GRASS version is 7.0.0.RC1. Here is the io console
 log:

 (Fri Mar 06 11:07:12 2015)
 v.overlay ainput=K3445_NEW_FISH at BGM_POLY_UTM
 binput=GRID_K34058 at BGM_POLY_UTM operator=and output=K34058_NEW
 Copying vector features from <K3445_NEW_FISH at BGM_POLY_UTM>...
 Copying vector features from <GRID_K34058 at BGM_POLY_UTM>...
 Snapping boundaries with 1e-008 ...
 Breaking lines...
 Removing duplicates...
 Cleaning boundaries at nodes...
 Breaking lines...
 Removing duplicates...
 Cleaning boundaries at nodes...
 Breaking lines...
 Removing duplicates...
 Cleaning boundaries at nodes...
 WARNING: Number of incorrect boundaries: 949
 Merging lines...
 Attaching islands...
 Building areas...
 25262 areas built
 7375 isles built
 Attaching islands...
 Number of nodes: 41280
 Number of primitives: 60697
 Number of points: 0
 Number of lines: 0
 Number of boundaries: 60697
 Number of centroids: 0
 Number of areas: 25262
 Number of isles: 7375
 Querying vector map <K3445_NEW_FISH at BGM_POLY_UTM>...
 Querying vector map <GRID_K34058 at BGM_POLY_UTM>...
 Writing centroids...
 WARNING: Unable to delete file
 'E:\mapping\GRASS7DB/BGM_POLY_UTM/BGM_POLY_UTM/.tmp/unknown/vector/tmp_7316/coor'
 Building topology for vector map <K34058_NEW at BGM_POLY_UTM>...
 Registering primitives...
 42630 primitives registered
 642388 vertices registered
 Building areas...
 12823 areas built
 2888 isles built
 Attaching islands...
 Attaching centroids...
 Number of nodes: 19873
 Number of primitives: 42630
 Number of points: 0
 Number of lines: 0
 Number of boundaries: 29808
 Number of centroids: 12822
 Number of areas: 12823
 Number of isles: 2888
 v.overlay complete.
 (Fri Mar 06 11:07:46 2015) Command finished (34 sec)

 Files can be provided upon request (~50 Mb).

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