QGIS: ESTRARRE CURVE DI LIVELLO

Pubblicato il da Totò

In questo breve articolo vedremo due approcci differenti per estrarre le curve di livello da un DTM (Modello digitale del terreno), in un primo caso estrarremo semplici curve di livello in 2D (cioè la quota sarà solo un attributo e non una geometria) e nel secondo caso la Z (quota) sarà anche nella geometria (curve di livello 3D).

importiamo un DTM in QGIS:

Immagine 5
DTM

dal menù Raster => Estrazione => Curve di livello …

Immagine 6
finestra di dialogo Curve di livello (2D)

con questa procedura si estraggono curve di livello (in 2D) con la possibilità di creare un attributo elevazione:

Immagine 7
Curve di livello e tabella attributi

Utilizziamo ora il processing di QGIS:

cerchiamo l’algoritmo di GRASS GIS r.contour.step

Immagine 9
Processing – QGIS

doppio clic e si aprirà una finestra di dialogo:

Immagine 10
finestra di dialogo – r.contour.step

otterremo:

Immagine 11
curve livello 3D

 

Immagine 12
Vista 3D – SAGA GIS

 

EDIT:
in QGIS 3.0 le curve di livello possono essere salvate in 3D

CURVE
screenshot 3.0 vs 2.18.17

video:

Donate:

 

 

 

 

Buon lavoro!!!

 

EDIT: se il blog post vi è piaciuto cliccate su ‘Mi piace’, grazie!!!
if you liked the blog post click on ‘Like’, thank you !!!

86 pensieri su “QGIS: ESTRARRE CURVE DI LIVELLO”

  1. Salve Totó stavo seguendo il tuo post per estrarre delle curve di livello 3D da DEM per poi esportarle in DXF, purtroppo non ho capito e non riesco a mantenere il valore z delle curve nell’esportazione, il file in DXF mi risulta quindi piatto, sapresti darmi alcuni consigli su come concludere con successo il procedimento?

    Grazie

    Piace a 1 persona

    Rispondi

  2. Dopo che hai estratto le curve di livello con ‘r.contour.step’ otterrai uno shapefile 3D; per esportarlo in DXF puoi seguire due strade:
    1. da Progetto => Esportazione DXF …, questa procedura (nativa in QGIS) appiattisce sempre tutti gli shapefile;
    2. tasto destro sul layer ‘curve di livello’ => ‘salva con nome’ e selezioni come formato DXF; questa procedura usa il metodo ‘OGR’ quindi esporta le geometrie senza modificarle, quindi otterrai un DXF 3D.
    Spero di essere stato chiaro.
    Saluti!!!

    "Mi piace"

    Rispondi

      1. Ciao,
        se hai seguito il mio post e ottieni sempre la Z nulla allora significa solo una cosa:
        nei tuo dataset NON esiste la Z, purtroppo questa è una possibile soluzione.

        saluti

        "Mi piace"

      15. Ho aperto il tuo file curve 3d.dxf con nanocad,
        cliccato su una curva di livello e nelle proprietà leggo z = 700 in due casi e z = 800 in altro caso.

        "Mi piace"

  3. Salve, ho seguito questo tutorial per estrarre delle curve di livello da un DTM5. Le curve sono veramente molto dettagliate. Se le esporto in dxf (mi servono per lavorare con il Cad), risultano ingestibili. Esiste un modo per approssimare queste curve in qgis, rendendole meno dettagliate?
    Andrea

    "Mi piace"

    Rispondi

    1. Ciao, per alleggerire lo shapefile puoi semplificare la geometria cioè diminuire sensibilmente il numero di vertici (ma perderai in precisione):
      menù Vettore => Strumenti di geometria => Semplifica geometria
      varia il parametro di semplificazione per esempio mettendo 1, noterai che il file diventerà molto più leggero.

      Saluti

      "Mi piace"

      Rispondi

  4. Ciao, innanzi tutto complimenti per il sito e per i tutorial, veramente esaustivi. Premetto che non sono un’esperto di qgis, ultimamente mi ci sto dedicando con la speranza di acquisire buona manualità sopratutto per il lavoro. Nella mia ricerca su internet per un problema che non riuscivo a risolvere mi ha portato su questa pagina, ho seguito il tutorial sull’estrazione delle curve di livello dai DTM molto simile a quello che cercavo, nel mio caso mi interesserebbe sapere come si possono estrarre i punti quotati da un DTM con passo 1m, credo che la procedura sia la stessa ma non riesco proprio a cavarne piede. Grazie in anticipo per l’aiuto che potrai darmi. Pierluigi

    "Mi piace"

    Rispondi

    1. Ciao Pierluigi,
      per risponderti in modo esaustivo devo capire bene cosa chiedi; intanto parlare di curve di livello o punti quotati non vi è molta differenza in quanto le curve sono un insieme di punti alla stessa quota;
      quindi basta configurare r.contour.step con passo 1 m e successivamente usare l’algoritmo ‘estrai vertici’; cosi facendo otterrai un layer di soli punti quotati.

      Se invece partissi da un layer di punti e vuoi associare la quota puoi usare il plugin ‘Point sampling tool’, naturalmente devi avere un DTM!!!

      spero di esserti stato utile.

      "Mi piace"

      Rispondi

  6. Ciao, ma se volessi estrarre curve di livello da una quota 0 da me stabilita con incremento progressivo visto che si tratta di una strada in leggera salita come devo fare?

    "Mi piace"

    Rispondi

    1. Puoi definire la tua quota di inizio contorno solo con il r.contour.step sostituendo lo ‘0’ con il tuo valore di partenza nel campo minimum contour level;
      oppure ritagli il DTM in modo da contenere solo la strada.

      "Mi piace"

      Rispondi

      1. ho risolto cosi, ho inserito il valore minimo nel mio caso -35 nel contour level, una volta inserite le etichette sono andato su proprietà -> campi -> elevation -> calcolatore campi ho selezionato il valore “elevation” (doppio Click) e inserito nel calcolare “ID elevetion” + 35 ( nel mio caso interessava che la quota – 35 fosse la quota 0) in automatico il calcolatore ha calcolato tutte le quote progressive con +35 cambiando la quota di riferimento

        "Mi piace"

  7. Buongiorno sto provando ad estrarre le curve in formato 3d attraverso il comando r.contour.step, ma durante il processo si blocca e da il seguente errore:

    The following layers were not correctly generated.
    Contours
    You can check the log messages to find more information about the execution of the algorithm

    "Mi piace"

    Rispondi

  8. Buongiorno,
    probabilmente sono sull’argomento sbagliato ma avrei bisogno di un aiuto.
    Ho un DTM di maglia 2m, ma vorrei modificarlo/ricrearlo con una maglia più grande (10m) poichè il livello di precisione è piuttosto basso e, soprattutto, l’area molto estesa quindi il programma ci mette un’eternità a girare. Ho provato (come scritto sulla guida di QGIS) a creare curve di livello e punti quotati con maglia 10 e successivamente a fare un’interpolazione triangolare (TIN) attraverso raster –> interpolazione, ma il programma mi si blocca ogni volta che provo a fare questa operazione. C’è mica un altro modo per creare o modificare la maglia del DTM?

    Ringrazio in anticipo per l’attenzione.
    Simona.

    "Mi piace"

    Rispondi

    1. Ciao,
      potresti provare questa procedura:
      Proprietà => Salva con nome; tra le varie opzioni trovi ‘Risoluzione’; scrivi 10 sia in orizzontale che verticale;
      avrai il tuo DTM in poco tempo. 😉

      "Mi piace"

      Rispondi

  9. ciao ti faccio una domanda: ho estratto un TIN dalle curve di livello. poi ho una base senza elevazione contente strade in shp. vorrei inserire l’elevazione del TIN dentro il dbf dello sdh strade. come posso fare?

    "Mi piace"

    Rispondi

    1. Per trasferire la terza dimensione ad uno shapefile (2D) in termini di geometria, quindi avere per ogni vertice la x,y e z potresti usare GRASS GIS (che è già installato in QGIS) tramite il processing e l’algoritmo ‘v.drape’; ti serve un DTM che contiene le quote; successivamente puoi popolare (dopo aver creato tre campi nuovi,x,y,z ) il dbf usando il calcolatore di campi e le seguenti formule:
      per i lcampo x: $x
      per il campo y: $y
      per la z: z($geometry)

      ciao

      "Mi piace"

      Rispondi

  10. Buon giorno, ho provato ad estrarre le curve di livello3D seguendo il tutorial ma ecco cosa mi dice il log message:

    Avvio dell’algoritmo r.contour.step – Create vector contours from raster at specified steps…
    g.proj -c proj4=”+proj=longlat +datum=WGS84 +no_defs”
    r.external input=”C:\Users\feded\Documents\Lavoro\Progetti\Smalzi Franco\DTM\dtm000171.tif” band=1 output=tmp14930275790411 –overwrite -o
    g.region n=5142010.0 s=5139990.0 e=640010.0 w=637990.0 res=1.0
    r.contour input=”tmp14930275790411″ minlevel=”0″ maxlevel=”10000″ step=”5″ cut=”0″ output=outputbf929147a93941269529e8b90d4140c1 –overwrite
    v.out.ogr -s -e input=outputbf929147a93941269529e8b90d4140c1 type=auto output=”C:/Users/feded/Documents/Lavoro/Progetti/Smalzi Franco/DTM” format=ESRI_Shapefile output_layer=3d_curve_5m –overwrite

    C:\PROGRA~1\QGIS2~1.18\bin>set HOME=C:\Users\feded

    C:\PROGRA~1\QGIS2~1.18\bin>set GISRC=C:\Users\feded\.qgis2\\processing\processing.gisrc7

    C:\PROGRA~1\QGIS2~1.18\bin>set WINGISBASE=C:\PROGRA~1\QGIS2~1.18\apps\grass\grass-7.2.0

    C:\PROGRA~1\QGIS2~1.18\bin>set GISBASE=C:\PROGRA~1\QGIS2~1.18\apps\grass\grass-7.2.0

    C:\PROGRA~1\QGIS2~1.18\bin>set GRASS_PROJSHARE=C:\PROGRA~1\QGIS2~1.18\apps\grass\grass-7.2.0\share\proj

    C:\PROGRA~1\QGIS2~1.18\bin>set GRASS_MESSAGE_FORMAT=plain

    C:\PROGRA~1\QGIS2~1.18\bin>if “” == “” set PATH=C:\PROGRA~1\QGIS2~1.18\apps\grass\grass-7.2.0\bin;C:\PROGRA~1\QGIS2~1.18\apps\grass\grass-7.2.0\lib;C:\PROGRA~1\QGIS2~1.18\apps\Python27\lib\site-packages\Shapely-1.2.18-py2.7-win-amd64.egg\shapely\DLLs;C:\PROGRA~1\QGIS2~1.18\apps\Python27\DLLs;C:\PROGRA~1\QGIS2~1.18\apps\qgis\bin;C:\PROGRA~1\QGIS2~1.18\apps\Python27\Scripts;C:\PROGRA~1\QGIS2~1.18\bin;C:\WINDOWS\system32;C:\WINDOWS;C:\WINDOWS\WBem

    C:\PROGRA~1\QGIS2~1.18\bin>if not “” == “” set PATH=C:\PROGRA~1\QGIS2~1.18\apps\grass\grass-7.2.0\bin;C:\PROGRA~1\QGIS2~1.18\apps\grass\grass-7.2.0\lib;;C:\PROGRA~1\QGIS2~1.18\apps\grass\grass-7.2.0\bin;C:\PROGRA~1\QGIS2~1.18\apps\grass\grass-7.2.0\lib;C:\PROGRA~1\QGIS2~1.18\apps\Python27\lib\site-packages\Shapely-1.2.18-py2.7-win-amd64.egg\shapely\DLLs;C:\PROGRA~1\QGIS2~1.18\apps\Python27\DLLs;C:\PROGRA~1\QGIS2~1.18\apps\qgis\bin;C:\PROGRA~1\QGIS2~1.18\apps\Python27\Scripts;C:\PROGRA~1\QGIS2~1.18\bin;C:\WINDOWS\system32;C:\WINDOWS;C:\WINDOWS\WBem

    C:\PROGRA~1\QGIS2~1.18\bin>set GRASS_VERSION=7.2.0

    C:\PROGRA~1\QGIS2~1.18\bin>if not “” == “” goto langset

    C:\PROGRA~1\QGIS2~1.18\bin>FOR /F “usebackq delims==” %i IN (`”C:\PROGRA~1\QGIS2~1.18\apps\grass\grass-7.2.0\etc\winlocale”`) DO @set LANG=%i

    C:\PROGRA~1\QGIS2~1.18\bin>set PATHEXT=.COM;.EXE;.BAT;.CMD;.VBS;.VBE;.JS;.JSE;.WSF;.WSH;.MSC;.PY

    C:\PROGRA~1\QGIS2~1.18\bin>set PYTHONPATH=;C:\PROGRA~1\QGIS2~1.18\apps\grass\grass-7.2.0\etc\python;C:\PROGRA~1\QGIS2~1.18\apps\grass\grass-7.2.0\etc\wxpython\n

    C:\PROGRA~1\QGIS2~1.18\bin>g.gisenv.exe set=”MAPSET=PERMANENT”

    C:\PROGRA~1\QGIS2~1.18\bin>g.gisenv.exe set=”LOCATION=temp_location”

    C:\PROGRA~1\QGIS2~1.18\bin>g.gisenv.exe set=”LOCATION_NAME=temp_location”

    C:\PROGRA~1\QGIS2~1.18\bin>g.gisenv.exe set=”GISDBASE=C:\Users\feded\AppData\Local\Temp\processinga92e796ac86f433daff9fc06c35e4c4c\grassdata”

    C:\PROGRA~1\QGIS2~1.18\bin>g.gisenv.exe set=”GRASS_GUI=text”

    C:\PROGRA~1\QGIS2~1.18\bin>g.proj -c proj4=”+proj=longlat +datum=WGS84 +no_defs”
    Default region was updated to the new projection, but if you have multiple mapsets `g.region -d` should be run in each to update the region from the default
    Informazioni sulla proiezione aggiornate

    C:\PROGRA~1\QGIS2~1.18\bin>r.external input=”C:\Users\feded\Documents\Lavoro\Progetti\Smalzi Franco\DTM\dtm000171.tif” band=1 output=tmp14930275790411 –overwrite -o
    ATTENZIONE: La verifica della proiezione verrà ignorata
    ERRORE: Valore di latitudine Nord non valido

    C:\PROGRA~1\QGIS2~1.18\bin>g.region n=5142010.0 s=5139990.0 e=640010.0 w=637990.0 res=1.0
    ERRORE: Input non valido

    C:\PROGRA~1\QGIS2~1.18\bin>r.contour input=”tmp14930275790411″ minlevel=”0″ maxlevel=”10000″ step=”5″ cut=”0″ output=outputbf929147a93941269529e8b90d4140c1 –overwrite
    ERRORE: Raster non trovato

    C:\PROGRA~1\QGIS2~1.18\bin>v.out.ogr -s -e input=outputbf929147a93941269529e8b90d4140c1 type=auto output=”C:/Users/feded/Documents/Lavoro/Progetti/Smalzi Franco/DTM” format=ESRI_Shapefile output_layer=3d_curve_5m –overwrite
    ERRORE: Vettoriale non trovato

    C:\PROGRA~1\QGIS2~1.18\bin>exit
    Converting outputs
    Caricamento layer

    The following layers were not correctly generated.
    Contours
    You can check the log messages to find more information about the execution of the algorithm

    Che cosa devo fare?!?
    Grazie!

    "Mi piace"

    Rispondi

  12. Salve, forse l’argomento non è attinente con questa sezione ma lo espongo ugualmente, lavorando su un vettoriale CTR al 10000 e con i DTM oramai ho acquisito dimistichezza, però ho notato una cosa, quando vado a inserire il plugin di Google Maps da Web>openLayersPlugin>Google Maps>Google Satellite la sovrapposizione non avviene correttamente, il satellite e il CTR al 10000 come pure il DTM sono sfasati e come effettuo lo zoom la cosa si amplifica. Mi sapete dire il perché? Grazie

    "Mi piace"

    Rispondi

  13. Il Plugin OpenLayer ormai non è più sviluppato e quindi lo shift è noto da tempo; prova ad usare il plugin QuickMapsSercices; poi da setting -> more services aggiungi le altre mappe (google, bing ecc..)
    ciao

    "Mi piace"

    Rispondi

  15. Ciao, ho un dsm con pixel di 3cm circa come posso evitare che le contours create da questa base intersechino alcuni file vettoriali lineari (es: strade) o poligonali (es:aree di rispetto) da me creati ?
    Se parto da punti quotati e creo un dtm con lo strumento interpolazione ho la possibilità di distinguere aree di restrizione, ecc. caricando i vari layer e specificandoli, ma se invece ho già il DTM/DSM ? Grazie mille!

    "Mi piace"

    Rispondi

    1. Ciao, la prima cosa che mi viene in mente è di estrarre dal DTM/DSM i punti e successivamente fare l’interpolazione.
      Per estrarre i punti dal raster potresti usare la ‘conversione di formato’ e salvare l’output in ASCII Gridded XYZ.
      Non ho mai avuto questa esigenza quindi non ho mai fatto prove!!!

      "Mi piace"

      Rispondi

  16. Ciao, e complimenti per questo sito.
    Pongo questo problema nella speranza che qualcuno mi possa dare una mano:
    Tramite il sito del Geoportale nazionale ho convertito il Lidar regione Lazio WMS (solo una parte) in file rester (TIF) al fine successivo di ricavare dal raster le curve di livello.
    Nelle proprietà dei TIF però non ritrovo più l’informazione corretta della quota (che invece interrogando il geoportale appare correttamente) ma trovo invece informazioni su 4 bande di colore. Estraendo le curve dai tali tif le info sulle quote sono per forza di cose errate.
    C’è un modo per recuperare la info quota giusta?
    Grazie

    Andrea

    "Mi piace"

    Rispondi

    1. A differenza di un WMS, che restituisce i dati territoriali come mappe statiche (rappresentate direttamente come immagini dal server), il WCS rende disponibili dati territoriali insieme con i loro attributi che possono essere interpretati, analizzati, elaborati, ecc.

      "Mi piace"

      Rispondi

  17. Ciao Totò , Grazie mille per il tutorial, anch’io provando ad estrarre le curve 3d da dem sono incappato nello stesso errore di FEDE. Ho provato a cambiare epsg ma senza successo. cosa può essere secondo te? Grazie mille

    "Mi piace"

    Rispondi

  18. Ciao Totò, innanzitutto grazie mille per il video che hai pubblicato! Ho un problema riguardo la lettura delle quote delle curve di livello: dopo aver esportato le isoipse estratte dal DTM in formato .dxf per poi utilizzarle su Autocad, ho notato che le quote ricavate sono diverse dallo shapefile che QGIS mi ha generato! Ad es., mentre su una curva di livello formato .shp ho una quota pari a 205 m, sul file .dxf aperto con Autocad la stessa curva di livello ha quota 57 m. Come mai?

    Piace a 1 persona

    Rispondi

  19. Salve,
    premetto che sono alle prime armi con QGIS, io ho seguito il tuo tutorial (veramente chiarissimo complimenti) per estrarre le Curve di livello 3D dal DTM Lombardia ma non genera le curve di livello, dopo aver eseguito il comando mi compare questo:
    Algoritmo in esecuzione…
    Algorithm ‘Curve di livello’ starting…
    Input parameters:
    { ‘BAND’ : 1, ‘CREATE_3D’ : True, ‘FIELD_NAME’ : ‘ELEV’, ‘IGNORE_NODATA’ : False, ‘INPUT’ : ‘contextualWMSLegend=0&crs=EPSG:4326&dpiMode=7&featureCount=10&format=image/jpeg&layers=DTM_5X5&styles&url=http://www.cartografia.servizirl.it/arcgis/services/wms/DTM5_RL_wms/MapServer/WMSServer’, ‘INTERVAL’ : 100, ‘NODATA’ : None, ‘OFFSET’ : 0, ‘OPTIONS’ : ”, ‘OUTPUT’ : ‘C:/Users/Utente/AppData/Local/Temp/processing_3cfc91b8b450413183183692e32ecffe/0a35093a63f640339188b337a0b4777f/OUTPUT.shp’ }

    GDAL command:
    gdal_contour -b 1 -a ELEV -i 100.0 -3d -f “ESRI Shapefile” contextualWMSLegend=0&crs=EPSG:4326&dpiMode=7&featureCount=10&format=image/jpeg&layers=DTM_5X5&styles&url=http://www.cartografia.servizirl.it/arcgis/services/wms/DTM5_RL_wms/MapServer/WMSServer C:/Users/Utente/AppData/Local/Temp/processing_3cfc91b8b450413183183692e32ecffe/0a35093a63f640339188b337a0b4777f/OUTPUT.shp
    GDAL command output:

    FAILURE: Missing destination filename.

    Usage: gdal_contour [-b ] [-a ] [-3d] [-inodata]

    [-snodata n] [-f ] [-i ]

    [[-dsco NAME=VALUE] …] [[-lco NAME=VALUE] …]

    [-off ] [-fl …]

    [-nln ] [-q]

    Traceback (most recent call last):
    File “C:/OSGEO4~1/apps/qgis/./python/plugins\processing\algs\gdal\GdalAlgorithm.py”, line 133, in processAlgorithm
    GdalUtils.runGdal(commands, feedback)
    File “C:/OSGEO4~1/apps/qgis/./python/plugins\processing\algs\gdal\GdalUtils.py”, line 110, in runGdal
    for line in proc.stdout:
    File “C:\OSGEO4~1\apps\Python37\lib\codecs.py”, line 322, in decode
    (result, consumed) = self._buffer_decode(data, self.errors, final)
    UnicodeDecodeError: ‘utf-8’ codec can’t decode byte 0x8a in position 10: invalid start byte

    Execution failed after 0.11 seconds

    Caricamento layer
    Algorithm ‘Curve di livello’ finished

    Da cosa potrebbe dipendere?
    Grazie mille davvero

    Piace a 1 persona

    Rispondi

  20. Si sto usando un WMS. Il portarle della Regione Lombardia mi fornisce solo quell’indirizzo URL per il DTM. Avrei la necessità di ottenere le curve di livello da quello ma non so come fare?
    Grazie

    Piace a 1 persona

    Rispondi

  21. premetto che sei un grande, sono un neofita e vorrei estrarre le curve di livello per poter fare carte della pendenza, ombreggiatura ecc. mi chiedo dove posso scaricare i DEM dell’area di mio interesse , zone limitate? ho provato a scaricare dal sito .usgs.gov/gdex/ per esempio la regione Molise ma purtroppo mi scarica tutta l’Italia.
    grazie

    Piace a 1 persona

    Rispondi

      1. Scusa se ti disturbo ancora, quando estraggo le curve di livello da un dem risultano un po, anzi molto, spigolose c’è un sistema per renderle più armoniche? inoltre se salvo le curve di livelli , etichettate con la quota, e le salvo come KLM quando le importo in google earth la quota sparisce., come faccio a far comparire anche la quota?
        Grazie Michele

        "Mi piace"

      2. Per rendere le linee più ‘lisce’ cerca in processing ‘Lisciatura’ e poi gioca con i parametri.
        Per il KML e le etichette non saprei consigliarti.
        saluti

        "Mi piace"

  22. Salve,
    Sto cercando di creare le curve di livello dal DTM della Regione Toscana (epsg: 3003) però continuo ad riscontrare questo errore:

    QGIS version: 3.6.2-Noosa
    Qt version: 5.11.2
    GDAL version: 2.4.1
    GEOS version: 3.7.1-CAPI-1.11.1 27a5e771
    PROJ version: 520
    Algoritmo in esecuzione…
    In avvio l’algoritmo ‘Curve di livello’ …
    Input parameters:
    { ‘BAND’ : 1, ‘CREATE_3D’ : False, ‘FIELD_NAME’ : ‘ELEV’, ‘IGNORE_NODATA’ : False, ‘INPUT’ : ‘/Users/alessandromasini/Desktop/DTM da regione.tif’, ‘INTERVAL’ : 10, ‘NODATA’ : None, ‘OFFSET’ : 0, ‘OPTIONS’ : ”, ‘OUTPUT’ : ‘TEMPORARY_OUTPUT’ }

    GDAL command:
    gdal_contour -b 1 -a ELEV -i 10.0 -f “ESRI Shapefile” “/Users/alessandromasini/Desktop/DTM da regione.tif” /var/folders/l6/9l0qvr094d9bnh3kkfxmjp_80000gn/T/processing_dbe35101022b4534a3fe5df35ca4a2f6/a43953cac6654ad3be9c8f9d9e75cdc4/OUTPUT.shp
    GDAL command output:
    /bin/sh: gdal_contour: command not found

    Execution completed in 0.17 seconds
    Results:
    {‘OUTPUT’: ‘/var/folders/l6/9l0qvr094d9bnh3kkfxmjp_80000gn/T/processing_dbe35101022b4534a3fe5df35ca4a2f6/a43953cac6654ad3be9c8f9d9e75cdc4/OUTPUT.shp’}

    Caricamento layer
    I seguenti layer non sono stati generati correttamente/var/folders/l6/9l0qvr094d9bnh3kkfxmjp_80000gn/T/processing_dbe35101022b4534a3fe5df35ca4a2f6/a43953cac6654ad3be9c8f9d9e75cdc4/OUTPUT.shpPuoi verificare il “Pannello dei messaggi” nella finestra principale di QGIS per trovare maggiori informazioni circa l’esecuzione dell’algoritmo

    Cosa potrebbe essere??
    Grazie per l’aiuto.

    Piace a 1 persona

    Rispondi

    1. Ciao una cosa che noto è che il nome del file ha degli spazi (DTM da regione.tif) prova a toglierli; poi l’errore ti consiglia di verificare nel “Pannello dei messaggi”.
      Prova a vedere cosa c’è scritto

      "Mi piace"

      Rispondi

      1. Ho provato a togliere gli spazi dal nome ma continua ha restituire il medesimo errore.
        Nel pannello messaggi trovo questo:

        2019-05-03T18:34:34 INFO gdal_contour -b 1 -a ELEV -i 10.0 -3d -f “ESRI Shapefile” /Users/alessandromasini/Desktop/DTM_da_regione.tif /var/folders/l6/9l0qvr094d9bnh3kkfxmjp_80000gn/T/processing_dbe35101022b4534a3fe5df35ca4a2f6/a54114dff2f74801a91cac1922b9d168/OUTPUT.shp
        2019-05-03T18:34:34 INFO GDAL execution console output
        /bin/sh: gdal_contour: command not found

        Sinceramente non capisco perchè..

        "Mi piace"

      2. L’ultima riga dell’errore c’è scritto che non trova il comando gdal_contour, forse non l’installazione di QGIS non è andata a buon fine.
        usando grass e r.contour.step ti funziona??

        "Mi piace"

  23. Si con il comando grass e r.contour.step funziona…
    il problema sembra nel Gdal perché tutti i comandi raster mi ridanno lo stesso avviso… forse perché sono su Mac???

    Piace a 1 persona

    Rispondi

  24. buongiorno, intanto complimenti per l’ottima e chiarissima spiegazione sia sui video che nelle risposte ai commenti.
    Sono una studentessa di architettura e per un progetto di paesaggio ho necessità di estrapolare curve da un file DTM. il file in questione è molto esteso, quindi ho tentato di estrarre una parte del DTM usando il comando Raster-estrazione-ritaglia raster con maschera. Ho inserito l’informazione di quale fosse il DTM sorgente e quale il layer da ritagliare, ma non riesco ad ottenere risultati poichè mi viene fuori un messaggio di errore.
    Aggiungo che mi sono assicurata di salvare il layer con il poligono di ritaglio in una cartella nota su cui mi è possibile scrivere, ma l’errore persiste. Forse sbaglio la procedura.

    qui il messaggio di errore nei Messaggi di log – processing

    2019-09-30T13:57:12 INFO gdalwarp -of GTiff -cutline “C:\\Users\\aless\\Google Drive\\00_ALESSIA\\ACCADEMIA\\02_SEMESTRE CORRENTE\\MSC1_SA_ATELIER COLLOMB II\\CAD Aufbau\\00_materiale iniziale\\planimetria da Qgis\\layer qgis lavoro\\test.shp” -cl test -crop_to_cutline “C:/Users/aless/Google Drive/00_ALESSIA/ACCADEMIA/02_SEMESTRE CORRENTE/MSC1_SA_ATELIER COLLOMB II/CAD da cattedra/00_Layer gis/Terreno raster qgis/nuova batimetria estuario_da alessandro/mntue0869_01m_file batimetria” C:/Users/aless/AppData/Local/Temp/processing_75c132ac091148dab3cc53f081d053d8/523ac1663c404a1ca16c18215564df27/OUTPUT.tif
    2019-09-30T13:57:12 INFO GDAL execution console output
    ERROR 1: Did not get any cutline features.

    Legato sempre a questo, vorrei capire se è normale che per estrarre curve ogni 20 cm da un DTM molto dettagliato (pesa più di 1 giga) che si estende per circa 50 Km può volerci più di un’ora (l’ho mandato un’ora fa e sta ancora calcolando, ma senza alcun avviso di errore)

    Grazie mille in anticipo, buona giornata
    Alessia

    Piace a 1 persona

    Rispondi

    1. Ciao Alessia,
      i problemi potrebbero essere tanti e il messaggio di log non mi aiuta tanto.
      Scrivimi una email dettaliata e se puoi allega i dati zippati (magari mettendo il dtm in un cloud e poimi alleghi link).
      Vai nella sezione ‘dovemi trovo’ e scrivi mail.
      ciao

      "Mi piace"

      Rispondi

      1. Alessia segui il mio consiglio invia tutto a Totò perchè è bravissimo e disponibilissimo. A me ha risolto tanti problemi.
        Saluti Michele

        Piace a 1 persona

      2. Salve,
        Ringrazio molto per la disponibilità e mi scuso per la tardività della mia risposta. Non so bene per quale ragione ma salvando semplicemente il file con un nuovo nome il problema si è risolto da solo.
        Mi scuso ancora per il disturbo e per non aver potuto rispondere prima, a causa di urgenti impegni.
        Grazie molte per la disponibilità,
        Saluti
        Alessia

        Piace a 1 persona

  26. Salve, sto cercando di estrarre dal DTM della Sicilia le curve di livello. Sto utilizzando il collegamento WCS (che inserisco tra i lager direttamente dal browser). Quando eseguo i passaggi, però, mi da come errore il seguente:

    Caricamento layer
    I seguenti layer non sono stati generati correttamente/private/var/folders/yg/dxz830ks29x52czssk2cmc400000gn/T/processing_gUwBqN/aca71c4b552e4ae798f96c5063035010/OUTPUT.gpkgPuoi verificare il “Pannello dei messaggi” nella finestra principale di QGIS per trovare maggiori informazioni circa l’esecuzione dell’algoritmo

    che nel pannello log è il seguente:

    2020-05-05T08:56:43 INFO GDAL execution console output
    /bin/sh: /private/var/folders/yg/dxz830ks29x52czssk2cmc400000gn/T/processing_gUwBqN/aca71c4b552e4ae798f96c5063035010/OUTPUT.gpkg: No such file or directory

    FAILURE: Missing destination filename.

    Usage: gdal_contour [-b ] [-a ] [-amin ] [-amax ]

    [-3d] [-inodata] [-snodata n] [-f ] [-i ]

    [[-dsco NAME=VALUE] …] [[-lco NAME=VALUE] …]

    [-off ] [-fl …] [-e ]

    [-nln ] [-q] [-p]

    Non saprei come risolvere.

    "Mi piace"

    Rispondi

    1. Ciao,
      non ha spiegato se hai salvato il DTM in locale oppure hai tentato di eseguire l’algoritmo direttamente usando il wcs.
      Per favore spiegati meglio, grazie.

      "Mi piace"

      Rispondi

  27. Buongiorno, quale è la procedura più efficiente e logica per estrarre curve di livello 3D di una data area da Qgis e importarle in Revit ? Io ho cercato di estrarre curve da Qgis in dxf (secondo il filmato da Lei proposto) ma cercando di caricarle in autocad quest’ultimo restituisce un messaggio di errore. Grazie per eventuali consigli e complimenti per la sua inarrivabile competenza.

    Piace a 1 persona

    Rispondi

    1. Ciao Luigi, non voglio polemizzare, ma i dati di input e il settaggio lo decidi tu, sei tu a dargli l’area da cui estrapolare le isolinee, sei tu a decidere il passo delle isolinee.
      Morale della favola: usa aree più piccole oppure imposta un passo alto.

      "Mi piace"

      Rispondi

  29. Scusami ma sono veramente inesperto. Ho caricato in Qgis un DTM dal geoportale nazionale. Ho sovrapposto ad esso la mappa di google. Ora vorrei ritagliare un pezzo del DTM riferito solo all’area di mio interesse ed estrarre le curve. Come si fa in Qgis a ritagliare il DTM o a fare in modo che il comando di estrazione cirve agisca solo su una zona ? Grazie

    Piace a 1 persona

    Rispondi

    1. Un modo per ritagliare un DTM:
      crea uno shapefile poligonale e traccia un rettangolo sull’area che ti interessa, poi, da strumenti di processing cerca Ritaglia il raster con maschera, usa il poligono appena creato come maschera

      "Mi piace"

      Rispondi

  30. Esiste qualche tutorial per capire il flusso di lavoro per la creazione dello Shape poligonale ? Ho fatto qualche prova ma non riesco ad applicare il ritaglio. Sono riuscito, invece, a caricare le curve nel sistema di riferimento proiettato. A tal proposito ti chiedo come fare per esportare oltre che il dxf delle curve anche l’ortofoto di google che in Qgis ho posto in trasparenza sotto il DTM. Saluti

    "Mi piace"

    Rispondi

Lascia un commento

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

Sono un ingegnere civile, nel 2008 seguii un corso di "Alta formazione per progettisti di Sistemi Informativi Territoriali" (SIT) dove iniziai a conoscere SO linux, SW GIS ( Autodesk, ESRI e FOSS) e mi 'innamorai' della 'Geomatica'.