[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:  closed                        
  Priority:  major        |   Milestone:                                
 Component:  Vector       |     Version:  7.0.0                         
Resolution:  fixed        |    Keywords:  v.overlay and polygons damaged
  Platform:  MSWindows 7  |         Cpu:  x86-64                        
--------------------------+-------------------------------------------------
Changes (by martinl):

  * status:  new => closed
  * resolution:  => fixed


Old description:

> 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).

New description:

 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#comment:3>
GRASS GIS <http://grass.osgeo.org>