[QGIS-it-user] Numero 15 Analisi di Visibilità - r.viewshed (Dario Saccoccia)

Dario Saccoccia saccoccia.dario a gmail.com
Dom 12 Apr 2020 12:46:42 PDT 

Ciao,

per correttezza e completezza riporto risposta giuntami da altra sede:

CiaoDario 
<https://www.facebook.com/profile.php?id=1536881856&hc_location=ufi>, 
un'analisi viewshed si basa sulla determinazione della linea di vista 
(il percorso ottico in "linea retta" fra due punti) per ogni cella del 
dem e la sua "intersezione" con il dem stesso. Se da un punto con 
coordinate x1,y1,z1 è "visibile" un punto con coordinate x2,y2,z2, 
allora sarà vero anche l'inverso. In questo senso, un'analisi viewshed è 
un'anlisi di intervisibilità.

Tuttavia solitamente i tool per l'analisi viewshed, come r.viewshed di 
GRASS, hanno dei parametri da impostare, per effettuare un calcolo in 
base alle proprie necessità, che possono rendere l'analisi asimmetrica.

Per esempio i parametri observer_elevation=value (Viewing elevation 
above the ground - Default: 1.75) e target_elevation=value (Offset for 
target elevation above the ground Default: 0.0).
Con i valori di dafault, la quota del punto di osservazione è aumentata 
di 1.75 metri (cioé l'altezza idealizzata di un uomo) mentre non lo è la 
quota degli altri punti del dem.
Probabilmente, invertire i valori di tali parametri dovrebbe portarti ad 
ottenere il risultato che hai indicato nel punto 2.

Per ulteriori dettagli potresti consultare la descrizione di 
r.viewshedhttps://grass.osgeo.org/grass78/manuals/r.viewshed.html 
<https://grass.osgeo.org/grass78/manuals/r.viewshed.html?fbclid=IwAR0RIwJsm4Fef7xicHsbuLPXT2q7Jjp53LH13QgYh-qSISxs_hNd7z8XR5c>e, 
per esempio, le seguenti risposte su 
gis.stackexchangehttps://gis.stackexchange.com/.../ztv-viewshed-analysis... 
<https://gis.stackexchange.com/questions/239886/ztv-viewshed-analysis-plugin-in-qgis?fbclid=IwAR24FM-mxG0vOiPtirazTI75zu5fpJzvs24ecLiE43T0bWvaORGUlLWvyJ4>ehttps://gis.stackexchange.com/.../creating-a-map... 
<https://l.facebook.com/l.php?u=https%3A%2F%2Fgis.stackexchange.com%2Fquestions%2F239584%2Fcreating-a-map-highlighting-areas-of-potential-visibility-of-a-high-structure%3Ffbclid%3DIwAR2nUsWx8Zt3zjLtQt85XZCs_g6JVXX4f1F3TRWmwLVzGRm86a59YzMOB7w&h=AT1PBQbBvtUmSClDloOU7O8P1Kfb5Av1NrkOE5jYOMIZqFleSmdg9sB92zLVImtYYmtTkkevGK5TmnTrcL7ehEnFysvdXUmCh2XJQkfP9JjacgMAYIKrmG6tPE8eD-EwEQ>


Ringraziando Andrea Giudicandrea che credo sia anche su questa lista, 
nel merito aggiungo quindi che:

1) il plug :" Visibility Analisys '' risolve questo quesito e anche 
altri, ovviamente comunque si Luca, si sarebbe potuto intervenire come 
suggerivi tu, ma mi chiedevo se si potesse     caricare direttamente un 
layer intero, e con Vibility Analisys, si può fare.

2) La risposta alla prima domanda diviene quasi ovvia: analizzando la 
visibilità di un punto è vero anche l'inverso, prestando un pò di 
attenzione ai parametri.

Grazie a tutti davvero.

Il 12/04/2020 21:00, qgis-it-user-request a lists.osgeo.org ha scritto:
> Invia le richieste di iscrizione alla lista QGIS-it-user all'indirizzo
> 	qgis-it-user a lists.osgeo.org
>
> Per iscriverti o cancellarti attraverso il web, visita
> 	https://lists.osgeo.org/mailman/listinfo/qgis-it-user
> oppure, via email, manda un messaggio con oggetto `help' all'indirizzo
> 	qgis-it-user-request a lists.osgeo.org
>
> Puoi contattare la persona che gestisce la lista all'indirizzo
> 	qgis-it-user-owner a lists.osgeo.org
>
> Se rispondi a questo messaggio, per favore edita la linea dell'oggetto
> in modo che sia più utile di un semplice "Re: Contenuti del digest
> della lista QGIS-it-user..."
>
>
> Argomenti del Giorno:
>
>     1. Re: Dubbio trasformazioni EPSG con GDAL 3.0.14 e	PROJ 6.3.1
>        (Luca Bellani)
>     2. Analisi di Visibilità - r.viewshed (Dario Saccoccia)
>     3. Re: Dubbio trasformazioni EPSG con GDAL 3.0.14 e	PROJ 6.3.1
>        (Luca Bellani)
>     4. Re: Dubbio trasformazioni EPSG con GDAL 3.0.14 e	PROJ 6.3.1
>        (Luca Bellani)
>     5. Re: Analisi di Visibilità - r.viewshed (Luca Bellani)
>     6. Re: Dubbio trasformazioni EPSG con GDAL 3.0.14 e	PROJ 6.3.1
>        (Andrea Giudiceandrea)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Sat, 11 Apr 2020 13:37:34 -0700 (MST)
> From: Luca Bellani <lucabellani.arqueologia a gmail.com>
> To: qgis-it-user a lists.osgeo.org
> Subject: Re: [QGIS-it-user] Dubbio trasformazioni EPSG con GDAL 3.0.14
> 	e	PROJ 6.3.1
> Message-ID: <1586637454589-0.post a n6.nabble.com>
> Content-Type: text/plain; charset=UTF-8
>
> Classe magistrale..!! no avevo idea (nemmeno mi ero accorto) di tutti questi
> dettagli, e dei colori dei rettangoli..
> Ora tutto sará piú chiaro! grazie mille!!
>
>
>
> --
> Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html
>
>
> ------------------------------
>
> Message: 2
> Date: Sat, 11 Apr 2020 22:40:45 +0200
> From: Dario Saccoccia <saccoccia.dario a gmail.com>
> To: qgis-it-user a lists.osgeo.org
> Subject: [QGIS-it-user] Analisi di Visibilità - r.viewshed
> Message-ID: <9b3e8484-92dd-15fb-017d-67dc4d767df3 a gmail.com>
> Content-Type: text/plain; charset=iso-8859-15; format=flowed
>
> Buonasera a tutti,
>
> scusate l'ignoranza e spero possiate essermi di aiuto.
>
> Mi è capita di utilizzare l'algoritmo r.viewshed senza problemi, ma mi
> ponevo due domande:
>
> 1) E' possibile fare la stessa analisi, ma utilizzando una linea? O una
> seria di punti su un tracciato? Senza voler ricorrere ai processi in
> serie o ai modelli, impostando di volta in volta un punto e poi usando
> il calcolatore raster per 'pulire' i risultati, capisco che la logica
> vorrebbe che su punti molto vicini l'analisi non cambi molto, ma intendo
> dire impostare ad esempio una sequenza di punti su un tracciato.
>
> 2) E fare un analisi di visibilità inversa? Cioè sapere da quali aree è
> visibile un punto?
>
> Grazie mille
> Buona serata e buona pasqua
>
>
>
> ------------------------------
>
> Message: 3
> Date: Sat, 11 Apr 2020 13:50:50 -0700 (MST)
> From: Luca Bellani <lucabellani.arqueologia a gmail.com>
> To: qgis-it-user a lists.osgeo.org
> Subject: Re: [QGIS-it-user] Dubbio trasformazioni EPSG con GDAL 3.0.14
> 	e	PROJ 6.3.1
> Message-ID: <1586638250865-0.post a n6.nabble.com>
> Content-Type: text/plain; charset=UTF-8
>
> Ho un dubbio:
>
> "Vedrai, poi, che se la tua mappa è in 4326 e importi il tuo layer in 4269,
> benché QGIS ti chieda quale trasformazione vuoi applicare tra 4269 e 4326,
> il crs della mappa non resterà 4326 ma verrà cambiato in 4269. La tua scelta
> verrà comunque memorizzata."
>
> Non dovrebbe essere al contrario? ovvero il crs del progetto, rimanere
> sempre lo stesso, e y layer adeguarsi?
>
> Grazie mille Andrea.
>
>
>
> --
> Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html
>
>
> ------------------------------
>
> Message: 4
> Date: Sat, 11 Apr 2020 13:57:59 -0700 (MST)
> From: Luca Bellani <lucabellani.arqueologia a gmail.com>
> To: qgis-it-user a lists.osgeo.org
> Subject: Re: [QGIS-it-user] Dubbio trasformazioni EPSG con GDAL 3.0.14
> 	e	PROJ 6.3.1
> Message-ID: <1586638679366-0.post a n6.nabble.com>
> Content-Type: text/plain; charset=UTF-8
>
> Con le tue spiegazioni credo aver decifrato il problema di vedere selezionata
> Africa.. essendo il progetto nuov, vuoto, immagino venga preso come
> referimento la coordinata 0,0 in quella zona. Perché se prima di importare
> il layer in questione, carico una mappa di fondo e mi sposto, tutto va ok!
> Grazie ancora
>
>
>
> --
> Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html
>
>
> ------------------------------
>
> Message: 5
> Date: Sat, 11 Apr 2020 14:04:52 -0700 (MST)
> From: Luca Bellani <lucabellani.arqueologia a gmail.com>
> To: qgis-it-user a lists.osgeo.org
> Subject: Re: [QGIS-it-user] Analisi di Visibilità - r.viewshed
> Message-ID: <1586639092296-0.post a n6.nabble.com>
> Content-Type: text/plain; charset=UTF-8
>
> Ciao, ho poca dimestichezza su questo tema, ma interessandomi molto,
> condivido cose che aveo letto n passato...
>
> 1) credo di si, nel quadro di dialogo del tool non c'é un bottoncino verde,
> che serve per rendere ricorsivo l'algoritmo su tutte le geometria del layer?
>
>
> 2) sono sicuro di si...ma non ho sotto mano informazioni precise, e non
> ricordo se nei plug di QGIS, ve ne era uno che permettesse l'analisi
> contrario..
>
> Hai provato con il plug :" Visibility Analisys " di QGIS? Mi pare permetta
> punti intermedi
> Spero possnao darti altre info pié precise, e seguo con interesse.
> Luca
>
>
>
> --
> Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html
>
>
> ------------------------------
>
> Message: 6
> Date: Sat, 11 Apr 2020 15:57:58 -0700 (MST)
> From: Andrea Giudiceandrea <andreaerdna a libero.it>
> To: qgis-it-user a lists.osgeo.org
> Subject: Re: [QGIS-it-user] Dubbio trasformazioni EPSG con GDAL 3.0.14
> 	e	PROJ 6.3.1
> Message-ID: <1586645878516-0.post a n6.nabble.com>
> Content-Type: text/plain; charset=UTF-8
>
> Luca Bellani wrote
>> Ho un dubbio:
>>
>> "Vedrai, poi, che se la tua mappa è in 4326 e importi il tuo layer in
>> 4269,
>> benché QGIS ti chieda quale trasformazione vuoi applicare tra 4269 e 4326,
>> il crs della mappa non resterà 4326 ma verrà cambiato in 4269. La tua
>> scelta
>> verrà comunque memorizzata."
> Avevo dimenticato di specificare che il comportamento descritto si riferisce
> ad un nuovo progetto in cui si importa un layer con crs diverso da quello
> del progetto, quando è impostata l'opzione predefinita "When a new project
> is created - Use CRS from first layer added" (in Settings->Options...->CRS)
>
> Diversamente, invece, accade modificando l'impostazione in "When a new
> project is created - Use a default CRS".
>
> A presto.
>
> Andrea
>
>
>
> --
> Sent from: http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html
>
>
> ------------------------------
>
> Subject: Chiusura del digest
>
> _______________________________________________
> QGIS-it-user mailing list
> QGIS-it-user a lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/qgis-it-user
>
>
> ------------------------------
>
> Fine di Digest di QGIS-it-user, Volume 54, Numero 15
> ****************************************************
-------------- parte successiva --------------
Un allegato HTML è stato rimosso...
URL: <http://lists.osgeo.org/pipermail/qgis-it-user/attachments/20200412/64c378bd/attachment-0001.html>

Maggiori informazioni sulla lista QGIS-it-user