Forse il precedente post era troppo lungo e non ha avuto risposta o chissà
per quale altro motivo! Ci riprovo: di seguito trovate l'errore che ottengo
mentre esporto il layer in postgis

ERROR:  duplicate key value violates unique constraint
"geo_veg_stat1_pkey"
DETAIL:  Key (ogc_fid)=(217) already exists.
Command: INSERT INTO "dbtre_2015"."geo_veg_stat1"
("wkb_geometry" , "ogc_fid" , "cat", "id_tab", "ty",
"uuid_padre", "code", "definition", "id_eunis", "attend")
VALUES ('0103000020787F00000100000009000000A31BE830D31C1C416
DA1E5AEEEE6524140B962EEBA1C1C41514D3A38EDE65241FD6190EE451C1
C41F2643AE8E8E65241C9D388EE991B1C4153283A98E1E652412B3F2D3E8
71B1C4115FAB46AE2E652415400BEEED01B1C41E5733AB8E4E652415E259
5EE371C1C41412C3A28E9E65241086E6CEE9D1C1C41A0CA39E8EDE65241A
31BE830D31C1C416DA1E5AEEEE65241'::GEOMETRY, 217 , 129301,
1616704, '060106_01-01', '9fe7bdb7-ac0b-
4d09-9235-6e2af2c8d225', '060106_01-01', 'vigneti', 5178,
90)
Warning 1: You've inserted feature with an already set FID
and that's perhaps the reason for the failure. If so, this
can happen if you reuse the same feature object for
sequential insertions. Indeed, since GDAL 1.8.0, the FID of
an inserted feature is got from the server, so it is not a
good ideato reuse it afterwards... All in all, try unsetting
the FID with SetFID(-1) before calling CreateFeature()
ERROR 1: INSERT command for new feature failed.
ERROR:  duplicate key value violates unique constraint
"geo_veg_stat1_pkey"

Avete qualche consiglio da darmi?

ecco i miei dati di sistema

GRASS version:
7.0.1
GRASS SVN Revision:
00000
Build Date:
2015-08-26
Build Platform:
x86_64-unknown-linux-gnu
GDAL/OGR:
1.11.2
PROJ.4:
4.8.0
GEOS:
3.4.2
SQLite:
3.8.2
Python:
2.7.6
wxPython:
2.8.12.1
Platform: Linux-3.13.0-66-generic-x86_64-with-Ubuntu-14.04-trusty

Saluti
Romina

Il giorno 13 ottobre 2015 15:47, Romina Di Paolo <romi.dipa...@gmail.com>
ha scritto:

> Ciao,
>
> vi chiedo aiuto perchè, anche dopo aver installato la macchina "OSGEO live
> 9", persistono i miei problemi su come far comunicare grass e postgis.
>
> I miei dati risiedono dentro uno schema specifico "dbtre_2015" su due
> layer geografici con geometria poligonale:
> geo_veg2 che ha tutta la vegetazione del Piemonte
> geo_noveg che ha i restanti temi (acque,urbanizzato e viabilità)
>
> su geo_veg2 mi servono le statistiche zonali sulla base dei Piani
> Territoriali Forestali mentre geo_noveg lo uso per avere la copertura del
> suolo completa.
>
> Ho importato geo_veg2 con il comando
>
>
>
> **v.in.ogr input="PG:host=xxx.it <http://xxx.it> dbname=nina
> user=r_dipaolo password=xxx" output=geo_veg2 layer=dbtre_2015.geo_veg2
> snap=1e-08*
> ed eseguito le statistiche zonali
>
>
>
> *v.rast.stats map=geo_veg2@siti_unesco raster=ptf_rstr_recls@siti_unesco
> column_prefix=ptf
> method=number,minimum,maximum,range,first_quartile,median,third_quartile,percentile*
> mentre per importare geo_noveg ho usato la procedura riassunta qui
> <http://geo.fsv.cvut.cz/data/osgeorel/2013-04-fem-gis/landa-fem-gis-2013-04.pdf>
> che finalmente ha messo chiarezza su come linkare un layer di PG senza
> doverlo importare in grass, ho quindi lanciato prima il comando
>
> *v.external e poi v.external.out* e tutto ha funzionato usando la riga di
> comando
>
> Il mio problema sorge quando devo esportare geo_veg2 su PG perchè perdo un
> bel pò di poligoni, so che il modello dati tra grass e Postgis è diverso
> come ampiamente documentato online  ma non so proprio come fare per avere i
> miei 224.000 poligoni su PG con tanto di statistiche zonali. Tra gli errori
> più comuni che ottengo nell'export ho:
> *Boundary/centroid. Skipped*
>
> *Feature is not a polygon. Skipping*
>
> Un altro dubbio sulla consistenza dei dati dentro grass è che per geo_veg2
> PostgreSQL conta 224050 poligoni mentre l'output di v.info mi dice questo
>
> Number of points:       866             Number of centroids:  206615
> Number of lines:        10221           Number of boundaries: 631710
> Number of areas:        237132          Number of islands:    32604
>
> Come faccio ad avere la sicurezza che siano stati importati tutti i
> poligoni?
>
> *Per geo_veg2 non ho usato v.external perchè dovendo usare le statistiche
> zonali che prevedono la conversione del vettoriale in raster mi da errore
> al momento della creazione del raster.
>
> Spero di essere riuscita ad inquadrare al meglio il problema perchè devo
> trovare una soluzione al più presto!
>
> Romina
>
>
> _______________________________________________
> GRASS-Italia mailing list
> GRASS-Italia@listserv.unipr.it
> http://listserv.unipr.it/mailman/listinfo/grass-italia
>
_______________________________________________
GRASS-Italia mailing list
GRASS-Italia@listserv.unipr.it
http://listserv.unipr.it/mailman/listinfo/grass-italia

Rispondere a