15 pensieri su “Importare shapefile in un db.sqlite utilizzando la shell”

1. Come al solito “stupefacente”!! ma …domanda posta da un mio collega su FB…per le vestizioni come dovremmo fare?

   "Mi piace""Mi piace"

   Rispondi

   1. Ciao,
      le vestizioni dipendono dal software GIS che usi: es: se usi QGIS è possibile creare un file .qml che contiene la vestizione del layer, ma è un file che legge solo QGIS; esistono altri formati di intescambio tipo SLD ma non sempre sono compatibili tra i vari software GIS in quanto NON esiste uno standard sulle vestizioni.
      Detto ciò è possibile salvare in un database (spatilite o postgis) lo stile dei vari layer in una tabella particolare ‘layer_styles’ in formato xml; per far ciò basti tematizzare un layer (proveniente da un database) e poi, dalle proprietà del layer, salva stile nel database.
      Rendere automatico questo processo, come descritto per gli shapefile, non ho mai avuto la necessità di farlo.
      ciao

      "Mi piace""Mi piace"

      Rispondi

2. Bellissimo. Ho però un problema: non trovo MSYS. Ricordo che quando minstallai le prime versioni di QGis in effetti mi figurava questa misteriosa icona di MSys, ma all’epoca ero alle primissime armi e QGis mi avanzava e mi bastava. Adesso ho installato la versione 2.18.9 e mi sono accorto che Msys non c’è più. Come posso fare? E’ nascosto da qualche parte o devo installarlo come cosa aggiuntiva? Grazie.

   "Mi piace""Mi piace"

   Rispondi

   1. Strano che non trovi l’applicativo, forse non compare più l’icona: prova a cercarlo in QGIS/apps/msys

      io ormai uso la versione OSGeo4W e non so più cosa accade nella standalone 😉

      "Mi piace""Mi piace"

      Rispondi

3. Per chi utilizzasse ubuntu o derivate può provare con questo script (da utilizzare nello stesso modo spiegato dall’incomparabile Totò, che ringrazio per il suo prezioso servizio di divulgazione):

   #!/bin/bash
   for f in `ls *.shp`
   do
   echo “Append shapefile ” $f ” …..”
   ogr2ogr -update -append output.sqlite $f -f “SQLite” -dsco SPATIALITE=YES -a_srs “EPSG:3003”
   done

   "Mi piace"Piace a 1 persona

   Rispondi

   1. Grazie,
      ma cosa significa “….”?

      "Mi piace""Mi piace"

      Rispondi

      1. Viene visualizzata la scritta “Append shapefile nomedelfile …..”

         "Mi piace"Piace a 1 persona

4. ciao
   io ho il problema contrario. devo convertire da sqlite a shp (ho una serie di data base spaziali sqlite e devo convertirli)
   mi potete aiutare?

   "Mi piace""Mi piace"

   Rispondi

   1. ciao,
      ho creato un altro script che esporta tutte le tabelle di un db sqlite in shapefile.
      segui tutti i passaggi dell articolo utilizzando questo script:

      
      

      	#!/bin/bash
      	# attivo il debug così vedo cosa avviene nella shell
      	set -x
      	# per ogni file sqlite contenuto contenuto nella cartella dello script, esegue
      	# l'export in shp nella cartella output
      	for i in *.sqlite;
      	do ogr2ogr -f "ESRI Shapefile" output ./"$i";
      	done

      view raw
      
      export_in_shp.sh
      
      hosted with ❤ by GitHub

      ciao

      "Mi piace""Mi piace"

      Rispondi

5. Ciao, molto interessante!
   Una variante possibile sempre tramite la shell, per la FUSIONE di tutti gli shp in un unico layer?
   (Gli shp ovviamente con la medesima struttura tabellare).
   Es: ho centinaia di shp di particelle catastali da fondere in un unico layer sqlite.

   "Mi piace"Piace a 1 persona

   Rispondi

   1. Ciao,
      prova con questo:
      https://github.com/pigreco/CLI-riga-di-comando-/blob/master/merge_shp.sh

      "Mi piace""Mi piace"

      Rispondi

6. Ho trovato un metodo alternativo (problemi di shell) usando solo sql spatialite:

   –FUSIONE MASSIVA DI PIU SHAPEFILES IN UN LAYER SPATIALITE

   –1) crea la tabella geometrica contenitore:
   create table fusion
   (pt_type text,
   pt_code text,
   pt_ccat integer,
   pt_fabb text);
   select addgeometrycolumn(‘fusion’,’geom’,3044,’polygon’,2);

   –2) elenca tutti gli shp che vuoi importare e mergiare:
   create virtual table “001_vl_parcel_poly”
   using virtualshape(‘c:\test\001_vl_parcel_poly’, ‘utf-8′, 3044);
   insert into fusion (pt_type, pt_code, pt_ccat, pt_fabb, geom)
   select pt_type, pt_code, pt_ccat, pt_fabb, geometry
   from “001_vl_parcel_poly”;
   delete from virts_geometry_columns where virt_name=’001_vl_parcel_poly’;
   drop table “001_vl_parcel_poly”;

   create virtual table “002_vl_parcel_poly”
   using virtualshape(‘c:\test\002_vl_parcel_poly’, ‘utf-8′, 3044);
   insert into fusion (pt_type, pt_code, pt_ccat, pt_fabb, geom)
   select pt_type, pt_code, pt_ccat, pt_fabb, geometry
   from “002_vl_parcel_poly”;
   delete from virts_geometry_columns where virt_name=’002_vl_parcel_poly’;
   drop table “002_vl_parcel_poly”;

   create virtual table “003_vl_parcel_poly”
   using virtualshape(‘c:\test\003_vl_parcel_poly’, ‘utf-8′, 3044);
   insert into fusion (pt_type, pt_code, pt_ccat, pt_fabb, geom)
   select pt_type, pt_code, pt_ccat, pt_fabb, geometry
   from “003_vl_parcel_poly”;
   delete from virts_geometry_columns where virt_name=’003_vl_parcel_poly’;
   drop table “003_vl_parcel_poly”;

   –3) crea indice spaziale per velocizzare la visualizzazione:
   select createspatialindex(‘fusion’,’geom’);
   vacuum

   "Mi piace"Piace a 1 persona

   Rispondi

   1. Grazie.

      "Mi piace""Mi piace"

      Rispondi

7. A me viene fuori sempre la scritta “ogr2ogr: comando non trovato”
   Come si aggiunge questo comando?

   "Mi piace"Piace a 1 persona

   Rispondi

   1. Ciao, in che ambiente lanci lo script suggerito??
      Normalmente ogr2ogr è già presente nella shell di MSYS

      Poi, scrivi maggiori dettagli altrimenti non posso aiutarti:(sistema operativo usato, shell usata, versione di QGIS ecc…)

      "Mi piace""Mi piace"

      Rispondi

Questo sito utilizza Akismet per ridurre lo spam. Scopri come vengono elaborati i dati derivati dai commenti.