[GRASS-user] Importing the GSHHS data set using v.in.ogr (grass_trunk)

Hi!

I am trying to import the GSHHS data set (full resolution, level 1, ESRI
Shapefile) using v.in.ogr in G7.

While it starts nicely, and progresses, slowly though in my machince (after
all, the data set is really large), it get's stuck at "Finding centroids
for...". Does it really require *too* much time?

Anybody imported this data set recently?

Best, Nikos

--%<--- --%<--- --%<--- --%<--- --%<--- --%<--- --%<--- --%<--- --%<--- --%<-

v.in.ogr dsn=/geo/geodata/GSHHS_shp/f/GSHHS_f_L1.shp output=gshhs_fl1Counting
polygons for 180446 features (OGR layer <GSHHS_f_L1>)...
100%
Importing 180446 features (OGR layer <GSHHS_f_L1>)...
100%
-----------------------------------------------------
Registering primitives...
191290 primitives registered
9471840 vertices registered
Number of nodes: 191290
Number of primitives: 191290
Number of points: 0
Number of lines: 0
Number of boundaries: 191290
Number of centroids: 0
Number of areas: -
Number of isles: -
-----------------------------------------------------
Cleaning polygons
-----------------------------------------------------
Breaking polygons...
Breaking polygons (pass 1: select break points)...
Warning: QMainWindow::saveState(): 'objectName' not set for QToolBar 0x36a97d0
'Vector field controller'
100%
Breaking polygons (pass 2: break at selected points)...
100%
-----------------------------------------------------
Removing duplicates...
100%
-----------------------------------------------------
Breaking boundaries...
100%
-----------------------------------------------------
Removing duplicates...
100%
-----------------------------------------------------
Cleaning boundaries at nodes...
100%
-----------------------------------------------------
Merging boundaries...
100%
-----------------------------------------------------
Removing dangles...
100%
-----------------------------------------------------
Building areas...
100%
180444 areas built
180443 isles built
Number of nodes: 192267
Number of primitives: 192269
Number of points: 0
Number of lines: 0
Number of boundaries: 192269
Number of centroids: 0
Number of areas: 180444
Number of isles: 180443
-----------------------------------------------------
Removing bridges...
100%
-----------------------------------------------------
Registering primitives...
180445 primitives registered
9458031 vertices registered
Building areas...
100%
180444 areas built
180443 isles built
Attaching islands...
100%
Number of nodes: 180444
Number of primitives: 180445
Number of points: 0
Number of lines: 0
Number of boundaries: 180445
Number of centroids: 0
Number of areas: 180444
Number of isles: 180443
-----------------------------------------------------
Finding centroids for OGR layer <GSHHS_f_L1>...
   0%
--->%-- --->%-- --->%-- --->%-- --->%-- --->%-- --->%-- --->%-- --->%-- --->%

The rest of the log pasted below...

On Wednesday 05 of June 2013 18:35:52 Nikos Alexandris wrote:

Hi!

I am trying to import the GSHHS data set (full resolution, level 1, ESRI
Shapefile) using v.in.ogr in G7.

While it starts nicely, and progresses, slowly though in my machince (after
all, the data set is really large), it get's stuck at "Finding centroids
for...". Does it really require *too* much time?

--%<--- --%<--- --%<--- --%<--- --%<--- --%<--- --%<--- --%<--- --%<---
v.in.ogr dsn=/geo/geodata/GSHHS_shp/f/GSHHS_f_L1.shp
output=gshhs_fl1Counting polygons for 180446 features (OGR layer
<GSHHS_f_L1>)...
100%
Importing 180446 features (OGR layer <GSHHS_f_L1>)...
100%
-----------------------------------------------------
Registering primitives...
191290 primitives registered
9471840 vertices registered
Number of nodes: 191290
Number of primitives: 191290
Number of points: 0
Number of lines: 0
Number of boundaries: 191290
Number of centroids: 0
Number of areas: -
Number of isles: -
-----------------------------------------------------
Cleaning polygons
-----------------------------------------------------
Breaking polygons...
Breaking polygons (pass 1: select break points)...
Warning: QMainWindow::saveState(): 'objectName' not set for QToolBar
0x36a97d0 'Vector field controller'
100%
Breaking polygons (pass 2: break at selected points)...
100%
-----------------------------------------------------
Removing duplicates...
100%
-----------------------------------------------------
Breaking boundaries...
100%
-----------------------------------------------------
Removing duplicates...
100%
-----------------------------------------------------
Cleaning boundaries at nodes...
100%
-----------------------------------------------------
Merging boundaries...
100%
-----------------------------------------------------
Removing dangles...
100%
-----------------------------------------------------
Building areas...
100%
180444 areas built
180443 isles built
Number of nodes: 192267
Number of primitives: 192269
Number of points: 0
Number of lines: 0
Number of boundaries: 192269
Number of centroids: 0
Number of areas: 180444
Number of isles: 180443
-----------------------------------------------------
Removing bridges...
100%
-----------------------------------------------------
Registering primitives...
180445 primitives registered
9458031 vertices registered
Building areas...
100%
180444 areas built
180443 isles built
Attaching islands...
100%
Number of nodes: 180444
Number of primitives: 180445
Number of points: 0
Number of lines: 0
Number of boundaries: 180445
Number of centroids: 0
Number of areas: 180444
Number of isles: 180443
-----------------------------------------------------
Finding centroids for OGR layer <GSHHS_f_L1>...
   0%
--->%-- --->%-- --->%-- --->%-- --->%-- --->%-- --->%-- --->%-- --->%--

-----------------------------------------------------
Writing centroids...
100%
-----------------------------------------------------
180446 input polygons
Total area: 1.49054E+14 (180444 areas)
-----------------------------------------------------
Copying features...
100%
Building topology for vector map <gshhs_fl1@PERMANENT>...
Registering primitives...
360889 primitives registered
9638475 vertices registered
Building areas...
100%
180444 areas built
180443 isles built
Attaching islands...
100%
Attaching centroids...
100%
Number of nodes: 180444
Number of primitives: 360889
Number of points: 0
Number of lines: 0
Number of boundaries: 180445
Number of centroids: 180444
Number of areas: 180444
Number of isles: 180443
-----------------------------------------------------
WARNING: Errors were encountered during the import
Try to import again, snapping with at least 1e-13: 'snap=1e-13'

O-K, I am trying v.in.gshhs (in G64) with the binary data.

Nikos

Nikos Alexandris wrote:

> I am trying to import the GSHHS data set (full resolution, level 1, ESRI
> Shapefile) using v.in.ogr in G7.
>
> While it starts nicely, and progresses, slowly though in my machince
> (after
> all, the data set is really large), it get's stuck at "Finding centroids
> for...". Does it really require *too* much time?

..

O-K, I am trying v.in.gshhs (in G64) with the binary data.

WoW! That was lightening fasssst!!

Thanks for that add-on, Nikos