Si ho provato anche con v.out.ogr ma me ne esporta ancora meno!

Il giorno 23 ottobre 2015 15:39, Roberto Marzocchi <
roberto.marzoc...@gmail.com> ha scritto:

> Effettivamente non avevo risposto sulla lista.. re-inoltro perchè magari
> qualcun altro intuisce meglio di me... Non mi convince
> l'opzione PRIMARY_KEY=NO perchè sul manuale dice "*disable adding primary
> key on FID column"* quindi mi sembra che voglia dire che non ti crea una
> nuova colonna FID ma una qualche chiave credo che te la crei ugualmente...
> e forse è quello che gli da fastidio...
> la colonna ogc_fid c'è già nel tuo vettoriale GRASS?
>
> hai già a usare v.out.ogr ?
> R
>
>
>
> Il giorno 23 ottobre 2015 15:11, Romina Di Paolo <romi.dipa...@gmail.com>
> ha scritto:
>
>>
>> Ciao forse la tua risposta non è arrivata in lista...
>>
>> In ogni caso per rispondere alle tue domande la tabella viene creata dal
>> comando v.out.postgis ex-novo e la sequenza probabilmente è implicita nel
>> comando ho usato anche le opzioni ogr
>>  del tipo PRIMARY_KEY=NO
>>  ecco il comando completo
>>
>> v.out.postgis --overwrite input=geo_veg type=area
>> output=PG:host=XXXXXXX dbname=XXXX" output_layer=dbtre_2015.stat_pg
>> options=PRIMARY_KEY=NO
>>
>> in ogni caso non vengono esportate circa 16000 per lo stesso motivo già
>> segnalato prima
>> DETAIL:  Key (ogc_fid)=(xxx) already exists.
>>
>> Comunque grazie!
>>
>>
>> Il giorno 23 ottobre 2015 14:52, Roberto Marzocchi <
>> roberto.marzoc...@gmail.com> ha scritto:
>>
>>> Il problema è la tabella postGIS e in particolare il fatto che trova un
>>> elemento con ogc_fid = 217 che già esiste ed è la tua chiave primaria (
>>> geo_veg_stat1_pkey):
>>>
>>> 1) la tabella è già piena e va prima pulita con un truncate
>>>
>>> 2) stai usando una campo serial in postgis che potrebbe essere collegato
>>> a una sequenza () che è inizializzata male
>>>
>>> R
>>>
>>>
>>> Il giorno 23 ottobre 2015 11:06, Romina Di Paolo <romi.dipa...@gmail.com
>>> > ha scritto:
>>>
>>>> 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
>>>>
>>>
>>>
>>
>
_______________________________________________
GRASS-Italia mailing list
GRASS-Italia@listserv.unipr.it
http://listserv.unipr.it/mailman/listinfo/grass-italia

Rispondere a