[QGIS-it-user] bug interpolazione processing?

QGIS QGIS-it-user
,
1

Salve, ho scaricato un file grezzo in formato LAZ del lidar messo a
disposizione dal geocatalogo della provincia di Trento [0].

L'ho convertito nel formato XYZ con PDAL da riga di comando, perchè ho
seguito quanto viene spiegato in questo blog [1], e mi sembra ben fatto,
anche se la versione di Qgis è vecchiotta e differisce un pò dalla mia
attuale.

Tra l'altro credo l'autore frequenti questa lista.... :slight_smile:

Il mio problema riguarda l'interpolazione, in particolare, dopo aver
caricato come layer testuale delimitato, e creatolo impostando le
coordinate nei due menu a discesa per scegliere i campi dedicati e
sistema di riferimento, avviando lo strumento di processing di
interpolazione, in particolare interpolazione TIN, carico il file XYZ
come vettore normalmente, però nel menu di "attributo interpolazione"
non succede nulla, menu a discesa bloccato, di qui il mio dubbio.

Provando a generare i raster, poi provando per esempio da processing,
Analisi Geomorfologica --> Rilievo, genera un raster nero senza rilievi.

Le informazioni del mio PC, Debian testing/experimental 64bit, server
grafico KDE, e versione Qgis:

Versione di QGIS
3.4.4-Madeira
Revisione codice QGIS
3331959
Compilato con Qt
5.11.3
Esecuzione con Qt
5.11.3
Compilato con GDAL/OGR
2.4.0
Esecuzione con GDAL/OGR
2.4.0
Compilato con GEOS
3.7.1-CAPI-1.11.1
Esecuzione con GEOS
3.7.0-CAPI-1.11.0 673b9939
Versione client PostgreSQL
11.1 (Debian 11.1-2)
Versione SpatiaLite
4.4.0-RC1
Versione QWT
6.1.3
Versione QScintilla2
2.10.4
Compiled against PROJ
520
Running against PROJ
5.2.0

Quindi sono a chiedere se vi risulta come bug o è, per forza di cosa
visto son niubbo, un mio limite, grazie.

[0]
https://siat.provincia.tn.it/geonetwork/srv/ita/catalog.search;jsessionid=60382A743F67AB616EC968D0258FB787#/metadata/p_tn:3203e1b3-5d73-47d8-b04e-9ab27c8538cf

[1]
https://massimilianomoraca.it/blog/usare-i-file-lidar-con-qgis-articolo-esercitazione?fbclid=IwAR2M12WJA0X4hJwNF7VxZvi5Kr89JJfkk69sT7dcha09a5CyBQt029l8mXY
--
Simone Girardelli
_|_|_|_|_|_|_|_|_|_
|_|_|_|_|_|_|_|_|_|_|

2

Non posso provare, érché il sito con i dati lidar non é accessibile, almeno
dall'estero, in questo momento.

--
Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html

3

Il 28/01/19 19:39, Luca Bellani ha scritto:

Non posso provare, érché il sito con i dati lidar non é accessibile, almeno
dall'estero, in questo momento.

Strano, prova qui direttamente nell'archivio dei singoli file?

https://siat.provincia.tn.it/IDT/raster/.public/lidar_2006_pat_laz/

--
Simone Girardelli
_|_|_|_|_|_|_|_|_|_
|_|_|_|_|_|_|_|_|_|_|

4

girarsi AT posteo DOT eu wrote

Il mio problema riguarda l'interpolazione, in particolare, dopo aver
caricato come layer testuale delimitato, e creatolo impostando le
coordinate nei due menu a discesa per scegliere i campi dedicati e
sistema di riferimento, avviando lo strumento di processing di
interpolazione, in particolare interpolazione TIN, carico il file XYZ
come vettore normalmente, però nel menu di "attributo interpolazione"
non succede nulla, menu a discesa bloccato, di qui il mio dubbio.

Ho appena fatto una prova utilizzando il primo file in elenco
https://siat.provincia.tn.it/IDT/raster/.public/lidar_2006_pat_laz/

funziona tutto senza problemi.
PS per trasformare il file .laz ho utilizzato QGIS 2.18.28 che ha ancora gli
strumenti lidar in processing
<http://osgeo-org.1560.x6.nabble.com/file/t378656/las.png&gt;

-----
https://pigrecoinfinito.wordpress.com/
--
Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html

5

Il 28/01/19 21:15, Totò ha scritto:

Ho appena fatto una prova utilizzando il primo file in elenco
https://siat.provincia.tn.it/IDT/raster/.public/lidar_2006_pat_laz/

funziona tutto senza problemi.
PS per trasformare il file .laz ho utilizzato QGIS 2.18.28 che ha ancora gli
strumenti lidar in processing
<http://osgeo-org.1560.x6.nabble.com/file/t378656/las.png&gt;

Col 2.18.28, non so qual'è in ordine di aggiornamento, mi pare
funzionava il menu di attributo interpolazione, mi ricordo di aver
aggiornato Debian e la stabile al 2.18.XX, non mi ricordo il numero, e
di lì mi si presentò il problema della disattivazione di quel menu, a
naso penso sia un bug, almeno su linux.

--
Simone Girardelli
_|_|_|_|_|_|_|_|_|_
|_|_|_|_|_|_|_|_|_|_|

6

poi far tutto da dentro pdal senza passare da altri backend via qgis

o sia… pdal ti permette di rasterizzare alla risoluzione e algoritmo di interpolazioen che desideri (writers.gdal)… ha anche dei filtri morfologici implementati direttametne o usabili via altri backend. PMF e SPMF son gia’ abbastanza potenti trovando gli opportuni parametri.

se poi il raster e’ nero… sciuro che la palette non sia scalata?

···

Luigi Pirelli

7

Il 28/01/19 21:44, Luigi Pirelli ha scritto:

se poi il raster e' nero... sciuro che la palette non sia scalata?

Luigi Pirelli

Sono sicuro solo della mia incompetenza, per cui non ti so rispondere di
preciso, avevo messo la spunta sul generatore automatico di colori o
quel che l'è nell'interfaccia del rilievo, ma veniva tutto nero.

Poi, con la 2.18, penultimo aggiornamento ricordo di esser riuscito a
creare il rilievo con il pacchetto dati messo a disposizione nel blog da
me linkato;
poi all'ultimo aggiornamento della 2.18, ricordo che si era disattivato
quel menu dell'interpolazione, e non esser riuscito a creare più niente.

Ora sono alla 3.4.4, installazione pulita, per cui non so, e non ti so
rispondere nel merito in maniera precisa, so solo che quel menu è ancora
disattivato/bloccato e non mi permette di scegliere ELEV.

--
Simone Girardelli
_|_|_|_|_|_|_|_|_|_
|_|_|_|_|_|_|_|_|_|_|

8

condividi un pezzo di raster e velocizzi le verifiche… dalle prove che hai fatto io capisco poco.

···

Luigi Pirelli

ps. che dice la statistica del raster?

9

Buongiorno,
secondo me è un problema di dato o di caricamento del dato in QGIS:
CIOè:
1. forse è stato convertito male da laz a xyz;
2. il file in xyz non ha intestazione e il separatore è uno spazio; quindi
l'importazione in QGIS, tramite aggiungi layer testo delimitato, deve essere
ben settata: EPSG 32632

allego un file .xyz che è la trasformazione del file laz; usa questo file
[0] e facci sapere:
quando importi in QGIS occhio all'intestazione e al separatore, inoltre stai
attendo se usi separatore decimale il punto o la virgola.

[0] https://mega.nz/#!Qc5gRQSY!HDAxfSIB9esGEIMaNIL5butwCkYhsBAQ6e5T7OD_dx8

-----
https://pigrecoinfinito.wordpress.com/
--
Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html

10

Totó, anche per la 3 ci sono, come plug.
Peró non riesco mai ad usarli, sará la mia incapacitá!

--
Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html

11

Ho eseguito in QGIS 3.4.4 la procedura spiegata da Totó, e tutto ok.
Solamente si é ripresentato l'errore nella visualizzazione del hillshade al
volo.
Devo ancora aprire il ticket, ma la foto lo dimostra
<http://osgeo-org.1560.x6.nabble.com/file/t381967/1.png&gt;

--
Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html

12

Luca Bellani wrote

Devo ancora aprire il ticket, ma la foto lo dimostra
&lt;http://osgeo-org.1560.x6.nabble.com/file/t381967/1.png&gt;

ciao Luca, confermo bug per la tematizzazione ombreggiatura.

per favore apri prima possibile un ticket e spiega che è un problema
relativo solo alla tematizzazione in quanto l'algoritmo gdal ombreggiatura
funziona bene.

saluti

-----
https://pigrecoinfinito.wordpress.com/
--
Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html

13

Aperto il ticket

https://issues.qgis.org/issues/21121

ciao

--
Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html

14

Il 29/01/19 08:34, Totò ha scritto:

Buongiorno,
secondo me è un problema di dato o di caricamento del dato in QGIS:
CIOè:
1. forse è stato convertito male da laz a xyz;
2. il file in xyz non ha intestazione e il separatore è uno spazio; quindi
l'importazione in QGIS, tramite aggiungi layer testo delimitato, deve essere
ben settata: EPSG 32632

allego un file .xyz che è la trasformazione del file laz; usa questo file
[0] e facci sapere:
quando importi in QGIS occhio all'intestazione e al separatore, inoltre stai
attendo se usi separatore decimale il punto o la virgola.

[0] https://mega.nz/#!Qc5gRQSY!HDAxfSIB9esGEIMaNIL5butwCkYhsBAQ6e5T7OD_dx8

Intanto grazie Totò, come scrivi sopra, facevo esattamente così, perchè
la pagina del blog che ho linkato nel primo commento alla discussione
aperta da me, così chiedeva per l'esercizio.

Ho provato il tuo file, e continuava a darmi problemi uguali a ieri.

Così ho deciso di disinstallare, pulire cercando cartelle, file ed
eventuali pacchetti orfani, e poi reinstallato la versione sperimentale,
e adesso funziona, probabilmente avevo qualche file o configurazione
vecchia che disturbava.

Non ho capito perchè dopo l'interpolazione lo strumento rilievo da
processing --> analisi geomorfologica, non mi permette di creare un
raster senza colorazione, ovvero con la classica scala di grigi, mi da
errore di mancato inserimento dei colori, se invece metto la spunta su "
Genera automaticamente classi di rilievo", allora funziona senza
problemi, oppure sono io che traviso lo strumento?

Se provo strumenti "ombreggiatura" o "esposizione (aspect)", funzionano,
mi generano un raster con la montagna.

Intanto grazie per la pazienza e l'aiuto, evidentemente non avevo pulito
bene la volta precedente, domani o dopodomani provo uno dei file da me
convertiti, e vi so dire.

grazie ancora.

--
Simone Girardelli
_|_|_|_|_|_|_|_|_|_
|_|_|_|_|_|_|_|_|_|_|

15

Il 29/01/19 08:34, Totò ha scritto:

Buongiorno,
secondo me è un problema di dato o di caricamento del dato in QGIS:
CIOè:
1. forse è stato convertito male da laz a xyz;
2. il file in xyz non ha intestazione e il separatore è uno spazio; quindi
l'importazione in QGIS, tramite aggiungi layer testo delimitato, deve essere
ben settata: EPSG 32632

allego un file .xyz che è la trasformazione del file laz; usa questo file
[0] e facci sapere:
quando importi in QGIS occhio all'intestazione e al separatore, inoltre stai
attendo se usi separatore decimale il punto o la virgola.

[0] https://mega.nz/#!Qc5gRQSY!HDAxfSIB9esGEIMaNIL5butwCkYhsBAQ6e5T7OD_dx8

Solo a segnalare che il caricamento del layer testuale del file da me
convertito, è avvenuta tranquillamente, mentre per il resto, per l'uso
degli strumenti interpolazione e rilievo, ho deciso di migliorare il mio
hardware, perchè ho visto che per file di grosse dimensioni, almeno per
me (1,4 GB), il mio pc ci mette una vita (elaboro con scheda grafica
integrata nel processore, AMD APU quad-core 3,7 Ghz, 8 GB di ram), circa
25 minuti solo per caricare il layer testuale.

Quindi vi saprò dire più avanti che risultati ottengo, intanto grazie
ancora a tutti.

--
Simone Girardelli
_|_|_|_|_|_|_|_|_|_
|_|_|_|_|_|_|_|_|_|_|