<div dir="ltr"><p class="gmail-MsoPlainText"><span lang="EN-US">Dear
all,<span></span></span></p>
<p class="gmail-MsoPlainText"><span lang="EN-US">I am
importing a fairly large polygon dataset (10 086 419 land cover polygons
covering the area of Finland) from PostGIS into GRASS 7.2.1 (Linux) using
v.import:</span></p>
<p class="gmail-MsoPlainText"><span lang="EN-US">v.import
--o input=PG:dbname=gisdata layer=zofie_cimburova.clip_finland output=clip_finland3
(see outputs of v.import, and <a href="http://v.info">v.info</a> at the end of the mail).<span></span></span></p>
<p class="gmail-MsoPlainText">The
import finishes successfully, but when visualizing the layer (i.e. Zoom to
selected maps), I get a following error:</p>
<p class="gmail-MsoPlainText">Failed
to run command 'd.vect map=clip_finland3@u_zofie.cimburova fcolor=200:200:200
color=black width=0 size=5 type=point,line,area,face icon=basic/x'. Details: d.vect:
spindex_rw.c:1362: rtree_get_node: Assertion `t->nb[level][which].n.level ==
level' failed.</p>
<p class="gmail-MsoPlainText">Do you
by any chance know what the problem might be? I thought it is caused by
overlaps in the dataset, but when I import a small subset of the original
dataset, the error does not occur. Is the problem in the amount of data? If so,
how do I deal with it?<br></p>
<p class="gmail-MsoPlainText">Thanks a
lot,</p><p class="gmail-MsoPlainText">Zofie</p>
<p class="gmail-MsoPlainText"><span lang="EN-US"> </span></p>
<p class="gmail-MsoPlainText"><span lang="EN-US"><a href="http://v.info">v.info</a> map=clip_finland3:<span></span></span></p>
<p class="gmail-MsoPlainText"><span lang="EN-US">(…)<span></span></span></p>
<p class="gmail-MsoPlainText"><span lang="EN-US">|----------------------------------------------------------------------------|<span></span></span></p>
<p class="gmail-MsoPlainText"><span lang="EN-US">| Map
format: native
|<span></span></span></p>
<p class="gmail-MsoPlainText"><span lang="EN-US">|----------------------------------------------------------------------------|<span></span></span></p>
<p class="gmail-MsoPlainText"><span lang="EN-US">| Type of map: vector (level: 2) |<span></span></span></p>
<p class="gmail-MsoPlainText"><span lang="EN-US">|
|<span></span></span></p>
<p class="gmail-MsoPlainText"><span lang="EN-US">| Number of points: 0 Number of centroids: 12653604
|<span></span></span></p>
<p class="gmail-MsoPlainText"><span lang="EN-US">| Number of lines: 0 Number of boundaries:
19187847 |<span></span></span></p>
<p class="gmail-MsoPlainText"><span lang="EN-US">| Number of areas: 12728699 Number of islands: 8610980
|<span></span></span></p>
<p class="gmail-MsoPlainText"><span lang="EN-US">|
|<span></span></span></p>
<p class="gmail-MsoPlainText"><span lang="EN-US">| Map is 3D: No
|<span></span></span></p>
<p class="gmail-MsoPlainText"><span lang="EN-US">| Number of dblinks: 1
|<span></span></span></p>
<p class="gmail-MsoPlainText"><span lang="EN-US">|
|<span></span></span></p>
<p class="gmail-MsoPlainText"><span lang="EN-US">| Projection: UTM (zone 33)
|<span></span></span></p>
<p class="gmail-MsoPlainText"><span lang="EN-US">| |<span></span></span></p>
<p class="gmail-MsoPlainText"><span lang="EN-US">| N: 7828443.03 S:
6602796.5606247 |<span></span></span></p>
<p class="gmail-MsoPlainText"><span lang="EN-US">| E: 1335791.79465162 W:
721051.43 |<span></span></span></p>
<p class="gmail-MsoPlainText"><span lang="EN-US">| |<span></span></span></p>
<p class="gmail-MsoPlainText"><span lang="EN-US">| Digitization threshold: 0
|<span></span></span></p>
<p class="gmail-MsoPlainText"><span lang="EN-US">| Comment:
|<span></span></span></p>
<p class="gmail-MsoPlainText"><span lang="EN-US">|
|<span></span></span></p>
<p class="gmail-MsoPlainText"><span lang="EN-US">+----------------------------------------------------------------------------+<span></span></span></p>
<p class="gmail-MsoPlainText"><span lang="EN-US"> </span></p>
<p class="gmail-MsoPlainText"><span lang="EN-US">v.import:<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US"> </span></p>
<p class="MsoNormal"><span lang="EN-US">Check if
OGR layer <zofie_cimburova.clip_finland> contains polygons...<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Importing
10086419 features (OGR layer <zofie_cimburova.clip_finland>)...<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">-----------------------------------------------------<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Registering
primitives...<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">19371733
primitives registered<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">519601629
vertices registered<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Number of
nodes: 12123276<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Number of
primitives: 19371733<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Number of
points: 0<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Number of
lines: 0<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Number of
boundaries: 19371733<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Number of
centroids: 0<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Number of
areas: -<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Number of
isles: -<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">-----------------------------------------------------<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Cleaning
polygons<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">-----------------------------------------------------<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Snapping
boundaries (threshold = 1.000e-13)...<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Reading
features...<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Snap
vertices Pass 1: select points<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Snap
vertices Pass 2: assign anchor vertices<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Snap
vertices Pass 3: snap to assigned points<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">-----------------------------------------------------<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Breaking
polygons...<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Breaking
polygons (pass 1: select break points)...<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Breaking
polygons (pass 2: break at selected points)...<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">-----------------------------------------------------<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Removing
duplicates...<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">-----------------------------------------------------<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Breaking
boundaries...<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">-----------------------------------------------------<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Removing
duplicates...<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">-----------------------------------------------------<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Cleaning
boundaries at nodes...<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">-----------------------------------------------------<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">(...)<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">-----------------------------------------------------<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Merging boundaries...<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">-----------------------------------------------------<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Removing
dangles...<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">-----------------------------------------------------<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Building
areas...<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">WARNING:
Area of size = 0.0 ignored<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">(...)<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">12728699
areas built<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">8610980
isles built<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Attaching
islands...<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Number of
nodes: 15022766<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Number of
primitives: 19187847<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Number of
points: 0<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Number of
lines: 0<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Number of
boundaries: 19187847<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Number of
centroids: 0<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Number of
areas: 12728699<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Number of
isles: 8610980<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">WARNING:
Number of incorrect boundaries: 52223<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">WARNING:
Vect_get_point_in_poly_isl(): collapsed area<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">WARNING:
Unable to calculate area centroid<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">(...)<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">-----------------------------------------------------<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Finding
centroids for OGR layer <zofie_cimburova.clip_finland>...<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">-----------------------------------------------------<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Writing
centroids...<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">WARNING:
2258413 areas represent more (overlapping) features, because polygons overlap
in input layer(s). Such areas are linked to more than 1 row in attribute table.
The number of features for those areas is stored as category in layer 2<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">-----------------------------------------------------<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">11478002
input polygons<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Total area:
3.93306E+11 (12728699 areas)<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Overlapping
area: 4.43546E+09 (2258413 areas)<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Area
without category: 0.1304 (6448 areas)<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">-----------------------------------------------------<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Copying
features...<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Building
topology for vector map <clip_finland3@u_zofie.cimburova>...<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Registering
primitives...<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">31841451
primitives registered<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">292562260
vertices registered<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Building
areas...<span></span></span></p>
<p class="MsoNormal">WARNING: Area of size = 0.0 ignored<span></span></p>
<p class="MsoNormal"><span lang="EN-US">(...)<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">12728699
areas built<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">8610980
isles built<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Attaching
islands...<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Attaching
centroids...<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Number of
nodes: 15022766<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Number of
primitives: 31841451<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Number of
points: 0<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Number of
lines: 0<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Number of
boundaries: 19187847<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Number of
centroids: 12653604<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Number of
areas: 12728699<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Number of
isles: 8610980<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">WARNING:
Number of incorrect boundaries: 52223<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">-----------------------------------------------------<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Some input
polygons are overlapping each other.<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">If
overlapping is not desired, the data need to be cleaned.<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">The input
could be cleaned by snapping vertices to each other.<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Estimated
range of snapping threshold: [1e-08, 1]<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Try to
import again, snapping with at least 1e-08: 'snap=1e-08'<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">Input
<PG:dbname=gisdata> successfully imported without reprojection<span></span></span></p>
<p class="MsoNormal"><span lang="EN-US">(Tue Jun 13
00:50:36 2017) Command finished (509 min 22 sec) <span></span></span></p></div>