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

GRASS GRASS user
,
1

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.

Thanks for your help.

--
Camille Bezzina
Geophom

2

On Thu, Apr 12, 2018 at 3:04 PM, Camille Bezzina <camille.bezzina@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@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

3

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.

Thanks.

(attachments)

v.import_verbose.txt (6.13 KB)

···

Le 12/04/2018 à 16:39, Markus Metz a écrit :

On Thu, Apr 12, 2018 at 3:04 PM, Camille Bezzina <camille.bezzina@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@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

4

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@geophom.fr <mailto:camille.bezzina@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@lists.osgeo.org <mailto:grass-user@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@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

5

On Thu, Apr 12, 2018 at 6:20 PM, Moritz Lennert <mlennert@club.worldonline.be> wrote:

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.

I would rather suggest snapping with a smaller value: snap=1e-4.

The main problem is close to the end, when building topology for the final output bati_D90_decoupe:
WARNING: Nombre de contours incorrects: 8772

This is the reason for disappearing polygons. This is a bug that has been fixed a while ago, please update your GRASS version to the latest 7.2 or 7.4 release.

Markus M

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@geophom.fr mailto:[camille.bezzina@geophom.fr](mailto:camille.bezzina@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@lists.osgeo.org mailto:[grass-user@lists.osgeo.org](mailto:grass-user@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@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

6

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

(attachments)

v.import_verbose.txt (4.48 KB)
difference.pdf (256 KB)

···

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@geophom.fr mailto:camille.bezzina@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@lists.osgeo.org mailto:grass-user@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@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

7

On Wed, Apr 18, 2018 at 5:27 PM, Camille Bezzina <camille.bezzina@geophom.fr> wrote:

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.

Can you provide off-list (a link to) the shapefile causing problems? Then we can test.

Markus M

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@geophom.fr mailto:[camille.bezzina@geophom.fr](mailto:camille.bezzina@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@lists.osgeo.org mailto:[grass-user@lists.osgeo.org](mailto:grass-user@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@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

8

You can download the shapefile causing problems here : https://we.tl/qByZnpMNvU

···

Le 18/04/2018 à 17:30, Markus Metz a écrit :

On Wed, Apr 18, 2018 at 5:27 PM, Camille Bezzina <camille.bezzina@geophom.fr> wrote:

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.

Can you provide off-list (a link to) the shapefile causing problems? Then we can test.

Markus M

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@geophom.fr mailto:[camille.bezzina@geophom.fr](mailto:camille.bezzina@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@lists.osgeo.org mailto:[grass-user@lists.osgeo.org](mailto:grass-user@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@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

–

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

9

On Wed, Apr 18, 2018 at 5:39 PM, Camille Bezzina <camille.bezzina@geophom.fr> wrote:

You can download the shapefile causing problems here : https://we.tl/qByZnpMNvU

Thanks a lot!

The input is 3D, and pairs of consecutive vertices are exactly on top of each other. The data can be imported successfully with v.in.ogr -2, forcing 2D output.

In this case I don’t understand (yet) why forcing 2D is necessary.

Markus M

Le 18/04/2018 à 17:30, Markus Metz a écrit :

On Wed, Apr 18, 2018 at 5:27 PM, Camille Bezzina <camille.bezzina@geophom.fr> wrote:

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.

Can you provide off-list (a link to) the shapefile causing problems? Then we can test.

Markus M

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@geophom.fr mailto:[camille.bezzina@geophom.fr](mailto:camille.bezzina@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@lists.osgeo.org mailto:[grass-user@lists.osgeo.org](mailto:grass-user@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@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

–
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

10

Thanks a lot Markus for your help !

It works. :slight_smile:

Last question about that, the “verbose” says me that I have 57851input polygons but 58049 areas in output. Why this difference ?

How can I check if all polygons are imported ?

···

Le 18/04/2018 à 18:55, Markus Metz a écrit :

On Wed, Apr 18, 2018 at 5:39 PM, Camille Bezzina <camille.bezzina@geophom.fr> wrote:

You can download the shapefile causing problems here : https://we.tl/qByZnpMNvU

Thanks a lot!

The input is 3D, and pairs of consecutive vertices are exactly on top of each other. The data can be imported successfully with v.in.ogr -2, forcing 2D output.

In this case I don’t understand (yet) why forcing 2D is necessary.

Markus M

Le 18/04/2018 à 17:30, Markus Metz a écrit :

On Wed, Apr 18, 2018 at 5:27 PM, Camille Bezzina <camille.bezzina@geophom.fr> wrote:

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.

Can you provide off-list (a link to) the shapefile causing problems? Then we can test.

Markus M

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@geophom.fr mailto:[camille.bezzina@geophom.fr](mailto:camille.bezzina@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@lists.osgeo.org mailto:[grass-user@lists.osgeo.org](mailto:grass-user@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@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

–
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

–

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

11

On Fri, Apr 20, 2018 at 10:47 AM, Camille Bezzina <camille.bezzina@geophom.fr> wrote:

Thanks a lot Markus for your help !

It works. :slight_smile:

Last question about that, the “verbose” says me that I have 57851input polygons but 58049 areas in output. Why this difference ?

The number of areas is larger than the number of polygons because holes inside polygons are also represented as areas.

How can I check if all polygons are imported ?

What counts is the number of centroids which is ideally identical to the number of input polygons, as long as input polygons don’t overlap. In case polygons overlap, the number of centroids can be different from the number of input polygons.

Markus M

Le 18/04/2018 à 18:55, Markus Metz a écrit :

On Wed, Apr 18, 2018 at 5:39 PM, Camille Bezzina <camille.bezzina@geophom.fr> wrote:

You can download the shapefile causing problems here : https://we.tl/qByZnpMNvU

Thanks a lot!

The input is 3D, and pairs of consecutive vertices are exactly on top of each other. The data can be imported successfully with v.in.ogr -2, forcing 2D output.

In this case I don’t understand (yet) why forcing 2D is necessary.

Markus M

Le 18/04/2018 à 17:30, Markus Metz a écrit :

On Wed, Apr 18, 2018 at 5:27 PM, Camille Bezzina <camille.bezzina@geophom.fr> wrote:

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.

Can you provide off-list (a link to) the shapefile causing problems? Then we can test.

Markus M

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@geophom.fr mailto:[camille.bezzina@geophom.fr](mailto:camille.bezzina@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@lists.osgeo.org mailto:[grass-user@lists.osgeo.org](mailto:grass-user@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@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

–
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

–
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