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