[GRASS-user] problem with g.in.ogr

Camille Bezzina camille.bezzina at geophom.fr
Wed Apr 18 08:27:12 PDT 2018


Hi,

I tried to import my vector with the latest version 7.4 and testing 
snap=0.01or other values but there is always the same problem.

I also tried with another shapefile (smaller) but it doen't change 
anything.

You can find the result in the attached file (verbose + .pdf). As you 
can see, when I import a .shp, a lot of polygons desappear.

Camille


Le 12/04/2018 à 18:20, Moritz Lennert a écrit :
> On 12/04/18 17:03, Camille Bezzina wrote:
>> Hi Markus,
>>
>> In the file attached you can find the full output of v.import with 
>> verbose.
>>
>> I launch this one with v.in.ogr and this is the answer.
>
> At the end you can see:
>
> "Some input polygons are overlapping each other.
> If overlapping is not desired, the data need to be cleaned.
> The input could be cleaned by snapping vertices to each other.
> Estimated range of snapping threshold: [1e-08, 1]
> Try to import again, snapping with at least 1e-08: 'snap=1e-08'"
>
> Try importing with snap=0.01 or something like that (depending on the 
> nature and precision of your data.
>
> Moritz
>
>>
>> Thanks.
>>
>>
>>
>> Le 12/04/2018 à 16:39, Markus Metz a écrit :
>>>
>>>
>>> On Thu, Apr 12, 2018 at 3:04 PM, Camille Bezzina 
>>> <camille.bezzina at geophom.fr <mailto:camille.bezzina at geophom.fr>> wrote:
>>> >
>>> > Hello all,
>>> >
>>> > I have a problem with v.in.ogr.
>>> > I would like to import a shape file (with 57851 polygons) with 
>>> v.in.ogr.
>>> >
>>> > v.import --verbose 
>>> input=/media/hdd1/phom/eolien/AIP_UNESCO_RONCHAMP/data/qgis/MNE/bati_D90_decoupe.shp 
>>> layer=bati_D90_decoupe output=bati_D90_decoupe -o
>>> >
>>> > Visibly the order is going well, my attribute table contains 57851 
>>> attributes but when I check my vector with display, some polygons 
>>> are missing.
>>>
>>> Can you post the full output of v.import with --verbose?
>>>
>>> Considering that you override the projection check with -o, you can 
>>> also use v.in.ogr directly and then post the full output of v.in.ogr 
>>> -o --verbose.
>>>
>>> Markus M
>>>
>>> >
>>> > Thanks for your help.
>>> >
>>> > --
>>> > Camille Bezzina
>>> > Geophom
>>> > _______________________________________________
>>> > grass-user mailing list
>>> > grass-user at lists.osgeo.org <mailto:grass-user at lists.osgeo.org>
>>> > https://lists.osgeo.org/mailman/listinfo/grass-user
>>>
>>
>> -- 
>> Camille BEZZINA
>> Chargé d'études géomatiques et photomontages
>> Geophom
>> 57 rue du Chemin Neuf 44521 OUDON
>> Standard: +33(0)2 85 52 02 59
>> Ligne directe: +33(0)9 72 56 81 71
>> www.geophom.fr
>>
>>
>> _______________________________________________
>> grass-user mailing list
>> grass-user at lists.osgeo.org
>> https://lists.osgeo.org/mailman/listinfo/grass-user
>>
>
>

-- 
Camille BEZZINA
Chargé d'études géomatiques et photomontages
Geophom
57 rue du Chemin Neuf 44521 OUDON
Standard: +33(0)2 85 52 02 59
Ligne directe: +33(0)9 72 56 81 71
www.geophom.fr
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20180418/c7fad1fd/attachment-0001.html>
-------------- next part --------------
(Wed Apr 18 17:23:37 2018)                                                      
v.import input=/media/hdd1/phom/eolien/AIP_UNESCO_RONCHAMP/data/qgis/MNE/bati_D90_decoupe_test_reduit1.shp layer=bati_D90_decoupe_test_reduit1 output=bati_D90_decoupe_test_reduit1 -o --overwrite
WARNING: Datum <Reseau_Geodesique_Francais_1993> non reconnu par GRASS et aucuns paramètres trouvés
WARNING: Datum <Reseau_Geodesique_Francais_1993> non reconnu par GRASS et aucuns paramètres trouvés
Over-riding projection check
Check if OGR layer <bati_D90_decoupe_test_reduit1> contains polygons...
WARNING: La couche vectorielle <bati_D90_decoupe_test_reduit1> existe déjà et va être écrasé
Importing 50 features (OGR layer <bati_D90_decoupe_test_reduit1>)...
-----------------------------------------------------
Enregistrement des primitives ...
592 primitives registered
1403 vertices registered
Nombre de nœuds: 559
Nombre de primitives: 592
Nombre de points: 0
Nombre de lignes :0
Nombre de contours: 592
Nombre de centroïdes: 0
Nombre de surfaces: -
Nombre d'îles: -
-----------------------------------------------------
Cleaning polygons
-----------------------------------------------------
Snapping boundaries (threshold = 1.000e-13)...
Reading features...
Accrochage des sommets passe 1 : sélection des points
Accrochage des sommets passe 2 : assignation des ancrages de sommets
Accrochage des sommets passe 3 : accrochage sur les points assignés
-----------------------------------------------------
Breaking polygons...
Breaking polygons (pass 1: select break points)...
Breaking polygons (pass 2: break at selected points)...
-----------------------------------------------------
Removing duplicates...
-----------------------------------------------------
Breaking boundaries...
-----------------------------------------------------
Removing duplicates...
-----------------------------------------------------
Cleaning boundaries at nodes...
-----------------------------------------------------
Breaking boundaries...
-----------------------------------------------------
Removing duplicates...
-----------------------------------------------------
Cleaning boundaries at nodes...
-----------------------------------------------------
Merging boundaries...
-----------------------------------------------------
Removing dangles...
-----------------------------------------------------
Construction des surfaces ...
51 areas built
38 isles built
Nombre de nœuds: 590
Nombre de primitives: 674
Nombre de points: 0
Nombre de lignes :0
Nombre de contours: 674
Nombre de centroïdes: 0
Nombre de surfaces: 51
Nombre d'îles: 38
WARNING: Nombre de contours incorrects: 22
-----------------------------------------------------
Removing bridges...
-----------------------------------------------------
Enregistrement des primitives ...
101 primitives registered
911 vertices registered
Construction des surfaces ...
30 areas built
30 isles built
Attachement des îles ...
Nombre de nœuds: 102
Nombre de primitives: 101
Nombre de points: 0
Nombre de lignes :0
Nombre de contours: 101
Nombre de centroïdes: 0
Nombre de surfaces: 30
Nombre d'îles: 30
WARNING: Nombre de contours incorrects: 71
-----------------------------------------------------
Finding centroids for OGR layer <bati_D90_decoupe_test_reduit1>...
-----------------------------------------------------
Writing centroids...
-----------------------------------------------------
50 input polygons
Total area: 142222 (30 areas)
Area without category: 26.48 (1 areas)
-----------------------------------------------------
Copie des objets ...
Construction de la topologie pour la carte vectorielle <bati_D90_decoupe_test_reduit1 at AIP_RONCHAMP_MNE>...
Enregistrement des primitives ...
130 primitives registered
940 vertices registered
Construction des surfaces ...
51 areas built
38 isles built
Attachement des îles ...
Attachement des centroïdes ...
Nombre de nœuds: 66
Nombre de primitives: 130
Nombre de points: 0
Nombre de lignes :0
Nombre de contours: 101
Nombre de centroïdes: 29
Nombre de surfaces: 51
Nombre d'îles: 38
WARNING: Nombre de contours incorrects: 22
-----------------------------------------------------
21 input polygons got lost during import.
The snapping threshold 1e-13 might be too large.
Estimated range of snapping threshold: [1e-08, 1]
Manual cleaning may be needed.
Input </media/hdd1/phom/eolien/AIP_UNESCO_RONCHAMP/data/qgis/MNE/bati_D90_decoupe_test_reduit1.shp> successfully imported without reprojection
(Wed Apr 18 17:23:39 2018) La commande s'est terminée (1 sec)                   

-------------- next part --------------
A non-text attachment was scrubbed...
Name: difference.pdf
Type: application/pdf
Size: 262578 bytes
Desc: not available
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20180418/c7fad1fd/attachment-0001.pdf>


More information about the grass-user mailing list