Ciao,
per correttezza e completezza riporto risposta giuntami da altra sede:
Ciao Dario, 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.viewshed https://grass.osgeo.org/grass78/manuals/r.viewshed.html e, per esempio, le seguenti risposte su gis.stackexchange https://gis.stackexchange.com/…/ztv-viewshed-analysis… e https://gis.stackexchange.com/…/creating-a-map…
Ringraziando Andrea Giudicandrea che credo sia anche su questa lista, nel merito aggiungo quindi che:
-
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.
-
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@lists.osgeo.org ha scritto:
Invia le richieste di iscrizione alla lista QGIS-it-user all'indirizzo [qgis-it-user@lists.osgeo.org](mailto:qgis-it-user@lists.osgeo.org) Per iscriverti o cancellarti attraverso il web, visita [https://lists.osgeo.org/mailman/listinfo/qgis-it-user](https://lists.osgeo.org/mailman/listinfo/qgis-it-user) oppure, via email, manda un messaggio con oggetto `help' all'indirizzo [qgis-it-user-request@lists.osgeo.org](mailto:qgis-it-user-request@lists.osgeo.org) Puoi contattare la persona che gestisce la lista all'indirizzo [qgis-it-user-owner@lists.osgeo.org](mailto:qgis-it-user-owner@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@gmail.com>](mailto:lucabellani.arqueologia@gmail.com) To: [qgis-it-user@lists.osgeo.org](mailto:qgis-it-user@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@n6.nabble.com>](mailto:1586637454589-0.post@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](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@gmail.com>](mailto:saccoccia.dario@gmail.com) To: [qgis-it-user@lists.osgeo.org](mailto:qgis-it-user@lists.osgeo.org) Subject: [QGIS-it-user] Analisi di Visibilità - r.viewshed Message-ID: [<9b3e8484-92dd-15fb-017d-67dc4d767df3@gmail.com>](mailto:9b3e8484-92dd-15fb-017d-67dc4d767df3@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@gmail.com>](mailto:lucabellani.arqueologia@gmail.com) To: [qgis-it-user@lists.osgeo.org](mailto:qgis-it-user@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@n6.nabble.com>](mailto:1586638250865-0.post@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](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@gmail.com>](mailto:lucabellani.arqueologia@gmail.com) To: [qgis-it-user@lists.osgeo.org](mailto:qgis-it-user@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@n6.nabble.com>](mailto:1586638679366-0.post@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](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@gmail.com>](mailto:lucabellani.arqueologia@gmail.com) To: [qgis-it-user@lists.osgeo.org](mailto:qgis-it-user@lists.osgeo.org) Subject: Re: [QGIS-it-user] Analisi di Visibilità - r.viewshed Message-ID: [<1586639092296-0.post@n6.nabble.com>](mailto:1586639092296-0.post@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](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@libero.it>](mailto:andreaerdna@libero.it) To: [qgis-it-user@lists.osgeo.org](mailto:qgis-it-user@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@n6.nabble.com>](mailto:1586645878516-0.post@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](http://osgeo-org.1560.x6.nabble.com/QGIS-Italian-User-f5250612.html) ------------------------------ Subject: Chiusura del digest _______________________________________________ QGIS-it-user mailing list [QGIS-it-user@lists.osgeo.org](mailto:QGIS-it-user@lists.osgeo.org) [https://lists.osgeo.org/mailman/listinfo/qgis-it-user](https://lists.osgeo.org/mailman/listinfo/qgis-it-user) ------------------------------ Fine di Digest di QGIS-it-user, Volume 54, Numero 15 ****************************************************