vector projection over wrapping boundary is split
|Reported by:||pertusus||Owned by:|
|Keywords:||LL, dateline, v.proj||Cc:|
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.