LAS FVG sovrapposizione tile bordi sfumati

[Mario Crevato]

Buongiorno a tutti, ho scaricato tutto una serie di LAS da Eagle FVG.

Con sagagis ho creato i DTM in formato ma effettuando il "mosaic" (merge) e facendo il rendering (hill shade, positive openess o qualsiasi altra elaborazione) i confini tra i vari tile vengono una schifezza, parecchi pixel sovrapposti ma sfocati. Ho elaborazioni batch via sagagis,  massive che uso da anni con i las sloveni (ARSO) e vengono perfetti (le elaborazioni grafiche non mostrano quasi le giunture tra i  tile ARSO sloveni che hanno hanno una sovrapposizione di un paio di pixel al max ma abbastanza affini).

Qualcuno di voi ci ha "sbattuto" su questo problema ?

grazie per l'attenzione

Mario

[Alberto Beinat]

Buongiorno Crevato,

ritengo che sia una questione di Saga, ma in ogni caso potrebbe cortesemente indicare alcuni dei LAZ in cui riscontra questo problema?

Le "tile" LAZ derivano dalla sovrapposizione e dal semplice "taglio" delle strisciate originali che altrimenti sarebbero troppo estese per essere gestite dagli utenti.

Come responsabile del collaudo della nuova cartografia, né io né i colleghi di UniUD, UniTS ed e-Laser abbiamo fin qui avuto evidenze di lacune a seguito del taglio dei LAS. Può comunque accadere che in quasi 8 TByte di dati Lidar possa annidarsi qualche problema. Se risulterà che il dato è difettoso in tal caso siamo in tempo per fare emendare i prodotti.

In relazione alle necessità degli utenti, segnalo inoltre che oltre ai LAZ sono già o saranno resi a breve disponibili su Eagle tutti i prodotti cartografici derivati dal Lidar: il DTM, due DSM (last e first pulse), un DBM - tutti a 0.5 metri pixel - come pure le mappe "hillshade", sia con quote ellissoidiche che sul livello mare (ortomteriche o quasi ...)

Buona giornata a tutte e a tutti

Alberto Beinat

--
Hai ricevuto questo messaggio perché sei iscritto al gruppo "QGIS_Utenti_FVG" di Google Gruppi.
Per annullare l'iscrizione a questo gruppo e non ricevere più le sue email, invia un'email a qgis_utenti_f...@googlegroups.com.
Per visualizzare questa discussione sul Web, visita https://groups.google.com/d/msgid/qgis_utenti_fvg/68ad95a9-2a02-4329-b444-9e842cb61ab1n%40googlegroups.com.

-- 
dott. ing. Alberto Beinat, Dott. Ric. (Ph.D.)
Ricercatore area 08/A4 - GEOMATICA
Dipartimento Politecnico di Ingegneria e Architettura (DPIA)
Università degli Studi di Udine
Via delle Scienze 206 - 33100 Udine UD
Tel. 0432 558708 - Fax. 0432 558052
Email: alberto...@uniud.it

[Mario Crevato]

Grazie della risposta, indagherò a livello parametri SagaGIS. I problema
si verifica su tutti i tile. Uso i stessi come peri i LAS sloveni ARSO
(per i quali non ci sono problemi ma lì la sovrapposizione è di 1-2
unità (metri), per questi del FVG noto che c'è una sovrapposizione e una
rotazione di parecchi metri. Ripeto, invece che fare da batch la
trasformazione ZLAS=>LAS=>POINTS=>GRID 0,25m la farò "a mano" in maniera
interattiva da SagaGIS per indagare dove sta l'inghippo.

Vi terrò informati sui risultati nel caso vi fossero altri utenti ad
avere il medesimo rpoblema.

grazie ancora

Mario Crevato

Il 20/09/2021 11:30, Alberto Beinat ha scritto:
> one a

[Mario Crevato]

Alberto grazie,

allora, ecco i step che eseguo in saga (riporto il comando richiamato usando PHP come motore batch):

1) converto il las in punti saga_cmd.exe io_shapes_las 1 -POINTS $THISDIR\\TEMP.subset.sg-pts -FILES $DIRIN\\$f
2) estraggo solo le classi 2,6,9,17 saga_cmd.exe pointcloud_tools 6 -INPUT $THISDIR\\TEMP.subset.sg-pts -ATTRIB classification -RESULT $THISDIR\\TEMP.subset1.sg-pts -MODE 1 -METHOD 2 -RETAB $THISDIR\\subset-class-table.txt -TOPERATOR 1
3) converto i punti in griglia interpolando i punti mancanti saga_cmd.exe grid_spline 4 -SHAPES $THISDIR\\TEMP.subset1.sg-pts -FIELD Z -METHOD 1 -TARGET_USER_SIZE $CELLSIZE -TARGET_OUT_GRID $DIRSGRD\\$F.sgrd
E qui (al punto 3) nasce il problema. L'area dei dati utili LAS è leggermente ruotata ripetto al proprio sistema di riferimento e quindi anche alla griglia di Saga e i "no-data" fuori bordo (nei triangoli risultanti dalla "rotazione") vengono semplicemente considerati da saga come dati mancanti e quindi interpolati secondo gli ultimi dati utili nell'area.

Devo consultarmi con l'autore di SagaGIS (Volker Wichmann) su come risolvere il problema e dal quale attendo una risposta

per ora grazie ancora

Mario Crevato
LIMES Arheologija prostora  Archeologia per il territorio

Il 20/09/2021 11:30, Alberto Beinat ha scritto:

Per visualizzare questa discussione sul Web, visita https://groups.google.com/d/msgid/qgis_utenti_fvg/6306f5c7-af3a-867e-580b-7cdb4156052b%40uniud.it.