[QGIS-it-user] Problema con r.viewshed
Carlo
naca84 a alice.it
Lun 8 Gen 2018 03:52:13 PST
Buongiorno a tutti
Scusate se riprendo un vecchio post, spero di aver cercato male, ma non
si era giunti ad una conclusione.
Il problema era ed è come fare l'intervisibilità da uno o più punti.
Visto che è una cosa che mi interessa (con ARCGIS andava tutto alla
perfezione), mi piacerebbe farlo con QGIS e GRASS o altre librerie.
Su GRASS la funzione che ho provato è r.viewshed
Al momento utilizzo QGIS 2.14.18, con GRASS 7 su Win8
I miei dati di partenza sono:
DEM (chiamato ritaglio, dimensione dica 9x9Km) di coordinate ED50_32N e
risoluzione 20m in formato .adf
Al momento dell'impostazione dei parametri inserisco il mio dem e le
coordinate del punto (prese in automatico sulla carta) ed il resto tutto
di default.
Questo è il risultato che ottengo:
The following layers were not correctly generated.
* Intervisibility
You can check the log messages to find more information about the
execution of the algorithm
Sottostante metto il log
Dove sbaglio?
Spero di avere una risposta
Grazie per la disponibilità
Buona giornata
Carlo
Questo è il risultato del log:
Avvio dell'algoritmo r.viewshed - Computes the viewshed of a point on an
elevation raster map....
g.proj -c proj4="+proj=utm +zone=32 +datum=WGS84 +units=m +no_defs"
r.external input="C:/Users/Carlo/Desktop/Intervisibilità/Ritaglio.adf"
band=1 output=tmp1514634485085 --overwrite -o
g.region n=4784300.41447 s=4775860.41447 e=878428.146484 w=869488.146484
res=20.0
r.viewshed input="tmp1514634485085" coordinates=874273.9051,4779833.1504
observer_elevation="1.75" target_elevation="0.0" max_distance="-1"
refraction_coeff="0.14286" memory="500"
output=output748429110e1f4714a5552789a2099a4d --overwrite
g.region raster=output748429110e1f4714a5552789a2099a4d
r.out.gdal --overwrite -c createopt="TFW=YES,COMPRESS=LZW"
input=output748429110e1f4714a5552789a2099a4d
output="C:\Users\Carlo\AppData\Local\Temp\processingbd0e790536f444528effa8279a2507f8\868aad2d005d4b458a06db9238ea81b9\output.tif"
C:\OSGeo4W64\bin>set HOME=C:\Users\Carlo
C:\OSGeo4W64\bin>set
GISRC=C:\Users\Carlo\.qgis2\\processing\processing.gisrc7
C:\OSGeo4W64\bin>set WINGISBASE=C:\OSGeo4W64\apps\grass\grass-7.2.1
C:\OSGeo4W64\bin>set GISBASE=C:\OSGeo4W64\apps\grass\grass-7.2.1
C:\OSGeo4W64\bin>set
GRASS_PROJSHARE=C:\OSGeo4W64\apps\grass\grass-7.2.1\share\proj
C:\OSGeo4W64\bin>set GRASS_MESSAGE_FORMAT=plain
C:\OSGeo4W64\bin>if "" == "" set
PATH=C:\OSGeo4W64\apps\grass\grass-7.2.1\bin;C:\OSGeo4W64\apps\grass\grass-7.2.1\lib;C:\OSGEO4~1\apps\Python27\lib\site-packages\Shapely-1.2.18-py2.7-win-amd64.egg\shapely\DLLs;C:\OSGEO4~1\apps\Python27\DLLs;C:\OSGEO4~1\apps\Python27\lib\site-packages\numpy\core;C:\OSGEO4~1\apps\qgis\bin;C:\OSGEO4~1\apps\grass\grass-7.2.2\lib;C:\OSGEO4~1\apps\grass\grass-7.2.2\bin;C:\OSGEO4~1\apps\Python27\Scripts;C:\OSGEO4~1\bin;C:\Windows\system32;C:\Windows;C:\Windows\system32\WBem
C:\OSGeo4W64\bin>if not "" == "" set
PATH=C:\OSGeo4W64\apps\grass\grass-7.2.1\bin;C:\OSGeo4W64\apps\grass\grass-7.2.1\lib;;C:\OSGeo4W64\apps\grass\grass-7.2.1\bin;C:\OSGeo4W64\apps\grass\grass-7.2.1\lib;C:\OSGEO4~1\apps\Python27\lib\site-packages\Shapely-1.2.18-py2.7-win-amd64.egg\shapely\DLLs;C:\OSGEO4~1\apps\Python27\DLLs;C:\OSGEO4~1\apps\Python27\lib\site-packages\numpy\core;C:\OSGEO4~1\apps\qgis\bin;C:\OSGEO4~1\apps\grass\grass-7.2.2\lib;C:\OSGEO4~1\apps\grass\grass-7.2.2\bin;C:\OSGEO4~1\apps\Python27\Scripts;C:\OSGEO4~1\bin;C:\Windows\system32;C:\Windows;C:\Windows\system32\WBem
C:\OSGeo4W64\bin>set GRASS_VERSION=7.2.1
C:\OSGeo4W64\bin>if not "" == "" goto langset
C:\OSGeo4W64\bin>FOR /F "usebackq delims==" %i IN
(`"C:\OSGeo4W64\apps\grass\grass-7.2.1\etc\winlocale"`) DO @set LANG=%i
C:\OSGeo4W64\bin>set
PATHEXT=.COM;.EXE;.BAT;.CMD;.VBS;.VBE;.JS;.JSE;.WSF;.WSH;.MSC;.PY
C:\OSGeo4W64\bin>set
PYTHONPATH=;C:\OSGeo4W64\apps\grass\grass-7.2.1\etc\python;C:\OSGeo4W64\apps\grass\grass-7.2.1\etc\wxpython\n
C:\OSGeo4W64\bin>g.gisenv.exe set="MAPSET=PERMANENT"
C:\OSGeo4W64\bin>g.gisenv.exe set="LOCATION=temp_location"
C:\OSGeo4W64\bin>g.gisenv.exe set="LOCATION_NAME=temp_location"
C:\OSGeo4W64\bin>g.gisenv.exe
set="GISDBASE=C:\Users\Carlo\AppData\Local\Temp\processingbd0e790536f444528effa8279a2507f8\grassdata"
C:\OSGeo4W64\bin>g.gisenv.exe set="GRASS_GUI=text"
C:\OSGeo4W64\bin>g.proj -c proj4="+proj=utm +zone=32 +datum=WGS84
+units=m +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:\OSGeo4W64\bin>r.external
input="C:/Users/Carlo/Desktop/Intervisibilità /Ritaglio.adf" band=1
output=tmp1514634485085 --overwrite -o
ERROR 4: `C:/Users/Carlo/Desktop/Intervisibilit+á/Ritaglio.adf' does not
exist in the file system,
and is not recognized as a supported dataset name.
C:\OSGeo4W64\bin>g.region n=4784300.41447 s=4775860.41447
e=878428.146484 w=869488.146484 res=20.0
C:\OSGeo4W64\bin>r.viewshed input="tmp1514634485085"
coordinates=874273.9051,4779833.1504 observer_elevation="1.75"
target_elevation="0.0" max_distance="-1" refraction_coeff="0.14286"
memory="500" output=output748429110e1f4714a5552789a2099a4d --overwrite
Computing events...
ERRORE: Raster map [tmp1514634485085] not found
C:\OSGeo4W64\bin>g.region raster=output748429110e1f4714a5552789a2099a4d
ERRORE: Mappa raster non trovata
C:\OSGeo4W64\bin>r.out.gdal --overwrite -c
createopt="TFW=YES,COMPRESS=LZW"
input=output748429110e1f4714a5552789a2099a4d
output="C:\Users\Carlo\AppData\Local\Temp\processingbd0e790536f444528effa8279a2507f8\868aad2d005d4b458a06db9238ea81b9\output.tif"
ERRORE: Mappa raster o gruppo non trovato
C:\OSGeo4W64\bin>exit
Converting outputs
Caricamento layer
The following layers were not correctly generated.
* Intervisibility
You can check the log messages to find more information about the
execution of the algorithm
-------------- parte successiva --------------
Un allegato HTML è stato rimosso...
URL: <http://lists.osgeo.org/pipermail/qgis-it-user/attachments/20180108/4c3d8005/attachment-0001.html>
Maggiori informazioni sulla lista
QGIS-it-user