[GRASS-dev] [GRASS GIS] #1527: vector projection over wrapping boundary is split

GRASS GIS trac at osgeo.org
Fri Dec 30 05:53:09 EST 2011

#1527: vector projection over wrapping boundary is split
 Reporter:  pertusus  |       Owner:  grass-dev@…              
     Type:  defect    |      Status:  new                      
 Priority:  normal    |   Milestone:  7.0.0                    
Component:  Vector    |     Version:  svn-trunk                
 Keywords:            |    Platform:  Linux                    
      Cpu:  x86-64    |  
 When I project a vector from a location thet overlaps with a wrapping
 boundary of a map, the vector is split in two.  For example, I have a
 longitude latitude location covering the whole world, that wraps somewhere
 in the Pacific (that is points on the eastern boundary are also on the
 western boundary).  When I project a vector constructed by doing a grid on
 a Lambert equal area location centered somewhere in the Pacific, the grid
 cells (that are, in that case, no more rectangular) that fall on the
 boundary are cut in two, with warnings like:

 WARNING: Number of centroids exceeds number of areas: 900 > 847
 WARNING: Number of incorrect boundaries: 507
 WARNING: Number of centroids outside area: 53

 I attach a tarball to reproduce the issue (there is a png showing the
 issue in the tarball too).  The script to do everything (setup the lambert
 location, and then the world map and do the png) is ./bug_wrap.sh.

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

More information about the grass-dev mailing list