[FOSSGIS-Talk] r.watherhed - SingleFlow
Lisa Römer
lisa.roemer at fischer-teamplan.de
Di Apr 13 09:45:30 CEST 2021
Hallo Liste,
ich habe ein Problem bei einer Fließwegberechnung mit der QGIS-Version 3.10.12 bzw. den Grass-Tools. Nachdem ich mit r.fill.dir die Senken gefüllt habe führe ich das Tool r.watherhed für die Fließwegeberechnung aus.
Standartmäßig ist hier Multiflow eingestellt. Damit funktioniert es auch.
Ich möchte aber gerne die SingleFlow/D8 Methode verwenden. Wenn ich das Häkchen dafür setze gibt es eine Fehlermeldung die ich nicht verstehe:
WARNUNG: G_spawn: unable to execute command
WARNUNG: Unterprozess fehlgeschlagen mit dem Fehlercode -1073741571
WARNUNG: Kategorie-Informationen für [accumulation300e15fd0c4948778d4120fd077831ae] in [PERMANENT] fehlen oder sind ungültig
......
.....
FEHLER: Rasterkarte <accumulation300e15fd0c4948778d4120fd077831ae> konnte nicht gefunden werden.
C:\Users\admin\Documents>r.out.gdal -t -m input="accumulation300e15fd0c4948778d4120fd077831ae" output="C:\temp\processing_fa1eef3b85314ede8bd2d19baccebc9b\c1def506b67e42ae9a1a9b4c5ecc052f\accumulation.tif" format="GTiff" createopt="TFW=YES,COMPRESS=LZW" --overwrite
....
....
Das gesamte Protokoll habe ich unten beigefügt.
Vielleicht kann mir da jemand einen Tipp geben.
Vielen Dank und viele Grüße,
Lisa
QGIS-Version: 3.10.12-A Coruña
QGIS-Codeversion: 75c848ffb1
Qt-Version: 5.11.2
GDAL-Version: 3.1.4
GEOS-Version: 3.8.1-CAPI-1.13.3
PROJ-Version: Rel. 6.3.2, May 1st, 2020
Verarbeite Algorithmus...
Algorithmus r.watershed startet...
Eingabeparameter:
{ '-4' : False, '-a' : False, '-b' : False, '-m' : False, '-s' : True, 'GRASS_RASTER_FORMAT_META' : '', 'GRASS_RASTER_FORMAT_OPT' : '', 'GRASS_REGION_CELLSIZE_PARAMETER' : 0, 'GRASS_REGION_PARAMETER' : None, 'accumulation' : 'TEMPORARY_OUTPUT', 'basin' : 'TEMPORARY_OUTPUT', 'blocking' : None, 'convergence' : 5, 'depression' : None, 'disturbed_land' : None, 'drainage' : 'TEMPORARY_OUTPUT', 'elevation' : 'C:/TMP/skibbe/Hydrologie/DepressionlessDEM2.tif', 'flow' : None, 'half_basin' : 'TEMPORARY_OUTPUT', 'length_slope' : 'TEMPORARY_OUTPUT', 'max_slope_length' : 1e+06, 'memory' : 300, 'slope_steepness' : 'TEMPORARY_OUTPUT', 'spi' : 'TEMPORARY_OUTPUT', 'stream' : 'TEMPORARY_OUTPUT', 'tci' : 'TEMPORARY_OUTPUT', 'threshold' : 10 }
g.proj -c proj4="+proj=utm +zone=32 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs +type=crs"
r.external input="C:\TMP\skibbe\Hydrologie\DepressionlessDEM2.tif" band=1 output="rast_607548a8584693" --overwrite -o
g.region n=5714000.0 s=5708000.0 e=382000.0 w=374000.0 res=1.0
r.watershed elevation=rast_607548a8584693 threshold=10 max_slope_length=1e+06 convergence=5 memory=300 -s accumulation=accumulation300e15fd0c4948778d4120fd077831ae drainage=drainage300e15fd0c4948778d4120fd077831ae basin=basin300e15fd0c4948778d4120fd077831ae stream=stream300e15fd0c4948778d4120fd077831ae half_basin=half_basin300e15fd0c4948778d4120fd077831ae length_slope=length_slope300e15fd0c4948778d4120fd077831ae slope_steepness=slope_steepness300e15fd0c4948778d4120fd077831ae tci=tci300e15fd0c4948778d4120fd077831ae spi=spi300e15fd0c4948778d4120fd077831ae --overwrite
g.region raster=accumulation300e15fd0c4948778d4120fd077831ae
r.out.gdal -t -m input="accumulation300e15fd0c4948778d4120fd077831ae" output="C:\temp\processing_fa1eef3b85314ede8bd2d19baccebc9b\c1def506b67e42ae9a1a9b4c5ecc052f\accumulation.tif" format="GTiff" createopt="TFW=YES,COMPRESS=LZW" --overwrite
g.region raster=drainage300e15fd0c4948778d4120fd077831ae
r.out.gdal -t -m input="drainage300e15fd0c4948778d4120fd077831ae" output="C:\temp\processing_fa1eef3b85314ede8bd2d19baccebc9b\020f69e886ce4560880d52cefe1f261b\drainage.tif" format="GTiff" createopt="TFW=YES,COMPRESS=LZW" --overwrite
g.region raster=basin300e15fd0c4948778d4120fd077831ae
r.out.gdal -t -m input="basin300e15fd0c4948778d4120fd077831ae" output="C:\temp\processing_fa1eef3b85314ede8bd2d19baccebc9b\0867bfdbd7d747aabe3ff96e33272d23\basin.tif" format="GTiff" createopt="TFW=YES,COMPRESS=LZW" --overwrite
g.region raster=stream300e15fd0c4948778d4120fd077831ae
r.out.gdal -t -m input="stream300e15fd0c4948778d4120fd077831ae" output="C:\temp\processing_fa1eef3b85314ede8bd2d19baccebc9b\4f820e95ee6441888a962be2231a4e9b\stream.tif" format="GTiff" createopt="TFW=YES,COMPRESS=LZW" --overwrite
g.region raster=half_basin300e15fd0c4948778d4120fd077831ae
r.out.gdal -t -m input="half_basin300e15fd0c4948778d4120fd077831ae" output="C:\temp\processing_fa1eef3b85314ede8bd2d19baccebc9b\69e271483af442a896f84c4284e42518\half_basin.tif" format="GTiff" createopt="TFW=YES,COMPRESS=LZW" --overwrite
g.region raster=length_slope300e15fd0c4948778d4120fd077831ae
r.out.gdal -t -m input="length_slope300e15fd0c4948778d4120fd077831ae" output="C:\temp\processing_fa1eef3b85314ede8bd2d19baccebc9b\6c728ab1d0e24d4aa2ec8fe48dc31022\length_slope.tif" format="GTiff" createopt="TFW=YES,COMPRESS=LZW" --overwrite
g.region raster=slope_steepness300e15fd0c4948778d4120fd077831ae
r.out.gdal -t -m input="slope_steepness300e15fd0c4948778d4120fd077831ae" output="C:\temp\processing_fa1eef3b85314ede8bd2d19baccebc9b\1540087275ef49c9a802adca684d7813\slope_steepness.tif" format="GTiff" createopt="TFW=YES,COMPRESS=LZW" --overwrite
g.region raster=tci300e15fd0c4948778d4120fd077831ae
r.out.gdal -t -m input="tci300e15fd0c4948778d4120fd077831ae" output="C:\temp\processing_fa1eef3b85314ede8bd2d19baccebc9b\89283763b9bb48178e0592f3656a2780\tci.tif" format="GTiff" createopt="TFW=YES,COMPRESS=LZW" --overwrite
g.region raster=spi300e15fd0c4948778d4120fd077831ae
r.out.gdal -t -m input="spi300e15fd0c4948778d4120fd077831ae" output="C:\temp\processing_fa1eef3b85314ede8bd2d19baccebc9b\7d57306837bc4e5fbb10fa28a3cd7809\spi.tif" format="GTiff" createopt="TFW=YES,COMPRESS=LZW" --overwrite
Starting GRASS GIS...
WARNUNG: Sperren gleichzeitiger Zugriffe auf ein Mapset ist unter Windows nicht möglich.
Cleaning up temporary files...
Executing <C:\temp\processing_fa1eef3b85314ede8bd2d19baccebc9b\grassdata\grass_batch_job.cmd> ...
C:\Users\admin\Documents>chcp 1252 1>NUL
C:\Users\admin\Documents>g.proj -c proj4="+proj=utm +zone=32 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs +type=crs"
Die Standard Region wurde auf die neue Projektion aktualisiert. Wenn Sie aber mehrere Mapsets haben, sollten Sie `g.region -d` in jedem ausführen, um die Einstellungen von der Standardregion zu übernehmen.
Projektionsinformationen aktualisiert
C:\Users\admin\Documents>r.external input="C:\TMP\skibbe\Hydrologie\DepressionlessDEM2.tif" band=1 output="rast_607548a8584693" --overwrite -o
Übersteuere Projektionsprüfung
Reading band 1 of 1...
Link to raster map <rast_607548a8584693> created.
C:\Users\admin\Documents>g.region n=5714000.0 s=5708000.0 e=382000.0 w=374000.0 res=1.0
C:\Users\admin\Documents>r.watershed elevation=rast_607548a8584693 threshold=10 max_slope_length=1e+06 convergence=5 memory=300 -s accumulation=accumulation300e15fd0c4948778d4120fd077831ae drainage=drainage300e15fd0c4948778d4120fd077831ae basin=basin300e15fd0c4948778d4120fd077831ae stream=stream300e15fd0c4948778d4120fd077831ae half_basin=half_basin300e15fd0c4948778d4120fd077831ae length_slope=length_slope300e15fd0c4948778d4120fd077831ae slope_steepness=slope_steepness300e15fd0c4948778d4120fd077831ae tci=tci300e15fd0c4948778d4120fd077831ae spi=spi300e15fd0c4948778d4120fd077831ae --overwrite
WARNUNG: G_spawn: unable to execute command
WARNUNG: Unterprozess fehlgeschlagen mit dem Fehlercode -1073741571
WARNUNG: Kategorie-Informationen für [accumulation300e15fd0c4948778d4120fd077831ae] in [PERMANENT] fehlen oder sind ungültig.
WARNUNG: Kategorie-Informationen für [tci300e15fd0c4948778d4120fd077831ae] in [PERMANENT] fehlen oder sind ungültig.
WARNUNG: Kategorie-Informationen für [spi300e15fd0c4948778d4120fd077831ae] in [PERMANENT] fehlen oder sind ungültig.
WARNUNG: Kategorie-Informationen für [drainage300e15fd0c4948778d4120fd077831ae] in [PERMANENT] fehlen oder sind ungültig.
WARNUNG: Kategorie-Informationen für [basin300e15fd0c4948778d4120fd077831ae] in [PERMANENT] fehlen oder sind ungültig.
WARNUNG: Kategorie-Informationen für [stream300e15fd0c4948778d4120fd077831ae] in [PERMANENT] fehlen oder sind ungültig.
WARNUNG: Kategorie-Informationen für [half_basin300e15fd0c4948778d4120fd077831ae] in [PERMANENT] fehlen oder sind ungültig.
WARNUNG: Kategorie-Informationen für [length_slope300e15fd0c4948778d4120fd077831ae] in [PERMANENT] fehlen oder sind ungültig.
WARNUNG: Kategorie-Informationen für [slope_steepness300e15fd0c4948778d4120fd077831ae] in [PERMANENT] fehlen oder sind ungültig.
C:\Users\admin\Documents>g.region raster=accumulation300e15fd0c4948778d4120fd077831ae
FEHLER: Rasterkarte <accumulation300e15fd0c4948778d4120fd077831ae> konnte nicht gefunden werden.
C:\Users\admin\Documents>r.out.gdal -t -m input="accumulation300e15fd0c4948778d4120fd077831ae" output="C:\temp\processing_fa1eef3b85314ede8bd2d19baccebc9b\c1def506b67e42ae9a1a9b4c5ecc052f\accumulation.tif" format="GTiff" createopt="TFW=YES,COMPRESS=LZW" --overwrite
FEHLER: Rasterkarte oder Gruppe <accumulation300e15fd0c4948778d4120fd077831ae> nicht gefunden.
C:\Users\admin\Documents>g.region raster=drainage300e15fd0c4948778d4120fd077831ae
FEHLER: Rasterkarte <drainage300e15fd0c4948778d4120fd077831ae> konnte nicht gefunden werden.
C:\Users\admin\Documents>r.out.gdal -t -m input="drainage300e15fd0c4948778d4120fd077831ae" output="C:\temp\processing_fa1eef3b85314ede8bd2d19baccebc9b\020f69e886ce4560880d52cefe1f261b\drainage.tif" format="GTiff" createopt="TFW=YES,COMPRESS=LZW" --overwrite
FEHLER: Rasterkarte oder Gruppe <drainage300e15fd0c4948778d4120fd077831ae> nicht gefunden.
C:\Users\admin\Documents>g.region raster=basin300e15fd0c4948778d4120fd077831ae
FEHLER: Rasterkarte <basin300e15fd0c4948778d4120fd077831ae> konnte nicht gefunden werden.
C:\Users\admin\Documents>r.out.gdal -t -m input="basin300e15fd0c4948778d4120fd077831ae" output="C:\temp\processing_fa1eef3b85314ede8bd2d19baccebc9b\0867bfdbd7d747aabe3ff96e33272d23\basin.tif" format="GTiff" createopt="TFW=YES,COMPRESS=LZW" --overwrite
FEHLER: Rasterkarte oder Gruppe <basin300e15fd0c4948778d4120fd077831ae> nicht gefunden.
C:\Users\admin\Documents>g.region raster=stream300e15fd0c4948778d4120fd077831ae
FEHLER: Rasterkarte <stream300e15fd0c4948778d4120fd077831ae> konnte nicht gefunden werden.
C:\Users\admin\Documents>r.out.gdal -t -m input="stream300e15fd0c4948778d4120fd077831ae" output="C:\temp\processing_fa1eef3b85314ede8bd2d19baccebc9b\4f820e95ee6441888a962be2231a4e9b\stream.tif" format="GTiff" createopt="TFW=YES,COMPRESS=LZW" --overwrite
FEHLER: Rasterkarte oder Gruppe <stream300e15fd0c4948778d4120fd077831ae> nicht gefunden.
C:\Users\admin\Documents>g.region raster=half_basin300e15fd0c4948778d4120fd077831ae
FEHLER: Rasterkarte <half_basin300e15fd0c4948778d4120fd077831ae> konnte nicht gefunden werden.
C:\Users\admin\Documents>r.out.gdal -t -m input="half_basin300e15fd0c4948778d4120fd077831ae" output="C:\temp\processing_fa1eef3b85314ede8bd2d19baccebc9b\69e271483af442a896f84c4284e42518\half_basin.tif" format="GTiff" createopt="TFW=YES,COMPRESS=LZW" --overwrite
FEHLER: Rasterkarte oder Gruppe <half_basin300e15fd0c4948778d4120fd077831ae> nicht gefunden.
C:\Users\admin\Documents>g.region raster=length_slope300e15fd0c4948778d4120fd077831ae
FEHLER: Rasterkarte <length_slope300e15fd0c4948778d4120fd077831ae> konnte nicht gefunden werden.
C:\Users\admin\Documents>r.out.gdal -t -m input="length_slope300e15fd0c4948778d4120fd077831ae" output="C:\temp\processing_fa1eef3b85314ede8bd2d19baccebc9b\6c728ab1d0e24d4aa2ec8fe48dc31022\length_slope.tif" format="GTiff" createopt="TFW=YES,COMPRESS=LZW" --overwrite
FEHLER: Rasterkarte oder Gruppe <length_slope300e15fd0c4948778d4120fd077831ae> nicht gefunden.
C:\Users\admin\Documents>g.region raster=slope_steepness300e15fd0c4948778d4120fd077831ae
FEHLER: Rasterkarte <slope_steepness300e15fd0c4948778d4120fd077831ae> konnte nicht gefunden werden.
C:\Users\admin\Documents>r.out.gdal -t -m input="slope_steepness300e15fd0c4948778d4120fd077831ae" output="C:\temp\processing_fa1eef3b85314ede8bd2d19baccebc9b\1540087275ef49c9a802adca684d7813\slope_steepness.tif" format="GTiff" createopt="TFW=YES,COMPRESS=LZW" --overwrite
FEHLER: Rasterkarte oder Gruppe <slope_steepness300e15fd0c4948778d4120fd077831ae> nicht gefunden.
C:\Users\admin\Documents>g.region raster=tci300e15fd0c4948778d4120fd077831ae
FEHLER: Rasterkarte <tci300e15fd0c4948778d4120fd077831ae> konnte nicht gefunden werden.
C:\Users\admin\Documents>r.out.gdal -t -m input="tci300e15fd0c4948778d4120fd077831ae" output="C:\temp\processing_fa1eef3b85314ede8bd2d19baccebc9b\89283763b9bb48178e0592f3656a2780\tci.tif" format="GTiff" createopt="TFW=YES,COMPRESS=LZW" --overwrite
FEHLER: Rasterkarte oder Gruppe <tci300e15fd0c4948778d4120fd077831ae> nicht gefunden.
C:\Users\admin\Documents>g.region raster=spi300e15fd0c4948778d4120fd077831ae
FEHLER: Rasterkarte <spi300e15fd0c4948778d4120fd077831ae> konnte nicht gefunden werden.
C:\Users\admin\Documents>r.out.gdal -t -m input="spi300e15fd0c4948778d4120fd077831ae" output="C:\temp\processing_fa1eef3b85314ede8bd2d19baccebc9b\7d57306837bc4e5fbb10fa28a3cd7809\spi.tif" format="GTiff" createopt="TFW=YES,COMPRESS=LZW" --overwrite
FEHLER: Rasterkarte oder Gruppe <spi300e15fd0c4948778d4120fd077831ae> nicht gefunden.
C:\Users\admin\Documents>exit
Execution of <C:\temp\processing_fa1eef3b85314ede8bd2d19baccebc9b\grassdata\grass_batch_job.cmd> finished.
Cleaning up temporary files...
Drücken Sie eine beliebige Taste . . .
Ausführung nach 64.54 Sekunden abgeschlossen
Ergebnisse:
{'accumulation': <QgsProcessingOutputLayerDefinition {'sink':TEMPORARY_OUTPUT, 'createOptions': {'fileEncoding': 'System'}}>,
'basin': <QgsProcessingOutputLayerDefinition {'sink':TEMPORARY_OUTPUT, 'createOptions': {'fileEncoding': 'System'}}>,
'drainage': <QgsProcessingOutputLayerDefinition {'sink':TEMPORARY_OUTPUT, 'createOptions': {'fileEncoding': 'System'}}>,
'half_basin': <QgsProcessingOutputLayerDefinition {'sink':TEMPORARY_OUTPUT, 'createOptions': {'fileEncoding': 'System'}}>,
'length_slope': <QgsProcessingOutputLayerDefinition {'sink':TEMPORARY_OUTPUT, 'createOptions': {'fileEncoding': 'System'}}>,
'slope_steepness': <QgsProcessingOutputLayerDefinition {'sink':TEMPORARY_OUTPUT, 'createOptions': {'fileEncoding': 'System'}}>,
'spi': <QgsProcessingOutputLayerDefinition {'sink':TEMPORARY_OUTPUT, 'createOptions': {'fileEncoding': 'System'}}>,
'stream': <QgsProcessingOutputLayerDefinition {'sink':TEMPORARY_OUTPUT, 'createOptions': {'fileEncoding': 'System'}}>,
'tci': <QgsProcessingOutputLayerDefinition {'sink':TEMPORARY_OUTPUT, 'createOptions': {'fileEncoding': 'System'}}>}
Lade Ergebnis Layer
Die folgenden Layer wurden nicht erzeugt.<ul><li>C:/temp/processing_fa1eef3b85314ede8bd2d19baccebc9b/020f69e886ce4560880d52cefe1f261b/drainage.tif</li><li>C:/temp/processing_fa1eef3b85314ede8bd2d19baccebc9b/0867bfdbd7d747aabe3ff96e33272d23/basin.tif</li><li>C:/temp/processing_fa1eef3b85314ede8bd2d19baccebc9b/c1def506b67e42ae9a1a9b4c5ecc052f/accumulation.tif</li></ul>Im 'Protokoll-Fenster' im QGIS-Hauptfenster sind mehr Informationen zur Ausführung des Algorithmus zu finden.
Mehr Informationen über die Mailingliste FOSSGIS-Talk-Liste