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.
* priority: normal => minor
* platform: Linux => All
* cpu: x86-32 => All
Comment:
Assertion happens only if the whole globe gets reprojected from lonlat to
laea, but not if only the current region extends of laea get reprojected
to laea.
Although it looks like the assertion prevented you from an illegal
reprojection, the assertion should not happen. The changes to the vector
spatial index in grass7 can not be backported because the spatial index in
grass7 has been nearly completely rewritten and is incompatible with the
grass6.x spatial index. The code in question has been replaced by another
algorithm written anew from scratch. In short, this is unlikely to get
fixed in grass6.x.
Replying to [comment:4 neteler]:
> Close as "wontfix" (for GRASS 6)?
Closed as "wontfix" because 1) it is quite difficult to reproject the
whole globe from latlon to laea without causing all sorts of numerical
instabilities since laea is not meant to represent the whole globe, 2) the
spatial index of GRASS 7 will not be backported.