[GRASS-user] QGIS + GRASS GIS console

GRASS GRASS user
, , ,
1

Hi,
QGIS + GRASS GIS 7.8.2, in Windows 10, 64 OS, not working, when computed a DTM from QGIS “TIN interpolation”, followed bi “r.fill.dir”:

QGIS version: 3.12.0-București

QGIS code revision: cd141490ec

Qt version: 5.11.2

GDAL version: 3.0.4

GEOS version: 3.8.1-CAPI-1.13.3

PROJ version: Rel. 6.3.1, February 10th, 2020

Processing algorithm…

Algorithm ‘r.fill.dir’ starting…

Input parameters:

{ ‘-f’ : False, ‘GRASS_RASTER_FORMAT_META’ : ‘’, ‘GRASS_RASTER_FORMAT_OPT’ : ‘’, ‘GRASS_REGION_CELLSIZE_PARAMETER’ : 0, ‘GRASS_REGION_PARAMETER’ : None, ‘areas’ : ‘TEMPORARY_OUTPUT’, ‘direction’ : ‘TEMPORARY_OUTPUT’, ‘format’ : 0, ‘input’ : ‘H:/FICHEIROS_DADOS/CARTOGRAFIA/MNT/MNT_1415/MNT_10K_2011/EPSG_3763/MDE_QGIS/mde2geral.tif’, ‘output’ : ‘TEMPORARY_OUTPUT’ }

g.proj -c proj4=“+proj=tmerc +lat_0=39.6682583333333 +lon_0=-8.13310833333333 +k=1 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs”

r.in.gdal input=“H:\FICHEIROS_DADOS\CARTOGRAFIA\MNT\MNT_1415\MNT_10K_2011\EPSG_3763\MDE_QGIS\mde2geral.tif” band=1 output=“rast_5e78c95a53bf92” --overwrite -o

g.region n=-55000.37998979999 s=-80000.43 e=-31998.2197489 w=-63998.17 res=4.9992111

r.fill.dir input=rast_5e78c95a53bf92 format=“grass” output=output9c34ebfb2d204df1afd1d27147a3522b direction=direction9c34ebfb2d204df1afd1d27147a3522b areas=areas9c34ebfb2d204df1afd1d27147a3522b --overwrite

g.region raster=output9c34ebfb2d204df1afd1d27147a3522b

r.out.gdal -t -m input=“output9c34ebfb2d204df1afd1d27147a3522b” output=“C:\Users\Particular\AppData\Local\Temp\processing_Dhabro\0081688cca90464486f44241eca07462\output.tif” format=“GTiff” createopt=“TFW=YES,COMPRESS=LZW” --overwrite

g.region raster=direction9c34ebfb2d204df1afd1d27147a3522b

r.out.gdal -t -m input=“direction9c34ebfb2d204df1afd1d27147a3522b” output=“C:\Users\Particular\AppData\Local\Temp\processing_Dhabro\54bb8b8a97534d53a90cfaa02ea9a703\direction.tif” format=“GTiff” createopt=“TFW=YES,COMPRESS=LZW” --overwrite

g.region raster=areas9c34ebfb2d204df1afd1d27147a3522b

r.out.gdal -t -m input=“areas9c34ebfb2d204df1afd1d27147a3522b” output=“C:\Users\Particular\AppData\Local\Temp\processing_Dhabro\734c1979cfb44de098ebea2a6ab5aad3\areas.tif” format=“GTiff” createopt=“TFW=YES,COMPRESS=LZW” --overwrite

Starting GRASS GIS…

WARNING: Concurrent mapset locking is not supported on Windows

Cleaning up temporary files…

Executing <C:\Users\Particular\AppData\Local\Temp\processing_Dhabro\grassdata\grass_batch_job.cmd> …

C:\Users\Particular\Documents>chcp 1252 1>NUL

C:\Users\Particular\Documents>g.proj -c proj4=“+proj=tmerc +lat_0=39.6682583333333 +lon_0=-8.13310833333333 +k=1 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +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

Projection information updated

C:\Users\Particular\Documents>r.in.gdal input=“H:\FICHEIROS_DADOS\CARTOGRAFIA\MNT\MNT_1415\MNT_10K_2011\EPSG_3763\MDE_QGIS\mde2geral.tif” band=1 output=“rast_5e78c95a53bf92” --overwrite -o

Over-riding projection check

Importing raster map <rast_5e78c95a53bf92>…

0…3…6…9…12…15…18…21…24…27…30…33…36…39…42…45…48…51…54…57…60…63…66…69…72…75…78…81…84…87…90…93…96…99…100

C:\Users\Particular\Documents>g.region n=-55000.37998979999 s=-80000.43 e=-31998.2197489 w=-63998.17 res=4.9992111

C:\Users\Particular\Documents>r.fill.dir input=rast_5e78c95a53bf92 format=“grass” output=output9c34ebfb2d204df1afd1d27147a3522b direction=direction9c34ebfb2d204df1afd1d27147a3522b areas=areas9c34ebfb2d204df1afd1d27147a3522b --overwrite

C:\Users\Particular\Documents>g.region raster=output9c34ebfb2d204df1afd1d27147a3522b

ERROR: Raster map not found

C:\Users\Particular\Documents>r.out.gdal -t -m input=“output9c34ebfb2d204df1afd1d27147a3522b” output=“C:\Users\Particular\AppData\Local\Temp\processing_Dhabro\0081688cca90464486f44241eca07462\output.tif” format=“GTiff” createopt=“TFW=YES,COMPRESS=LZW” --overwrite

ERROR: Raster map or group not found

C:\Users\Particular\Documents>g.region raster=direction9c34ebfb2d204df1afd1d27147a3522b

ERROR: Raster map not found

C:\Users\Particular\Documents>r.out.gdal -t -m input=“direction9c34ebfb2d204df1afd1d27147a3522b” output=“C:\Users\Particular\AppData\Local\Temp\processing_Dhabro\54bb8b8a97534d53a90cfaa02ea9a703\direction.tif” format=“GTiff” createopt=“TFW=YES,COMPRESS=LZW” --overwrite

ERROR: Raster map or group not found

C:\Users\Particular\Documents>g.region raster=areas9c34ebfb2d204df1afd1d27147a3522b

ERROR: Raster map not found

C:\Users\Particular\Documents>r.out.gdal -t -m input=“areas9c34ebfb2d204df1afd1d27147a3522b” output=“C:\Users\Particular\AppData\Local\Temp\processing_Dhabro\734c1979cfb44de098ebea2a6ab5aad3\areas.tif” format=“GTiff” createopt=“TFW=YES,COMPRESS=LZW” --overwrite

ERROR: Raster map or group not found

C:\Users\Particular\Documents>exit

Execution of <C:\Users\Particular\AppData\Local\Temp\processing_Dhabro\grassdata\grass_batch_job.cmd> finished.

Cleaning up temporary files…

Press any key to continue . . .

Execution completed in 62.54 seconds

Results:

{‘areas’: <QgsProcessingOutputLayerDefinition {‘sink’:TEMPORARY_OUTPUT, ‘createOptions’: {‘fileEncoding’: ‘System’}}>,
‘direction’: <QgsProcessingOutputLayerDefinition {‘sink’:TEMPORARY_OUTPUT, ‘createOptions’: {‘fileEncoding’: ‘System’}}>,
‘output’: <QgsProcessingOutputLayerDefinition {‘sink’:TEMPORARY_OUTPUT, ‘createOptions’: {‘fileEncoding’: ‘System’}}>}

Loading resulting layers

The following layers were not correctly generated.

  • C:/Users/Particular/AppData/Local/Temp/processing_Dhabro/0081688cca90464486f44241eca07462/output.tif
  • C:/Users/Particular/AppData/Local/Temp/processing_Dhabro/54bb8b8a97534d53a90cfaa02ea9a703/direction.tif
  • C:/Users/Particular/AppData/Local/Temp/processing_Dhabro/734c1979cfb44de098ebea2a6ab5aad3/areas.tif
You can check the ‘Log Messages Panel’ in QGIS main window to find more information about the execution of the algorithm.

What should be?

Thanks in advance

···

Cumprimentos,
Valter Albino - Geógrafo Físico, M.Sc.

Modelação H&H / Riscos ambientais / OT&U

www.valteralbino.wixsite.com/hydrodynamics

2

In the first time (previous email), tried to compute “r.fill.dir” from “Processing Toolbox”. I also tried from GRASS GIS console: the output is:
"
C:/OSGEO4~1/apps/grass/grass78/bin/r.fill.dir.exe input=mde2geral@TRABALHOS format=grass output=mde2geralfill direction=mde2geralflowdir

Reading input elevation raster map…

Module crashed or killed

"

Valter Albino <valteralbino@gmail.com> escreveu no dia segunda, 23/03/2020 à(s) 14:51:

···

Cumprimentos,
Valter Albino - Geógrafo Físico, M.Sc.

Modelação H&H / Riscos ambientais / OT&U

www.valteralbino.wixsite.com/hydrodynamics

3

Looks like the “TIN interpolation” raster can’t be read by the “r.fill.dir”.

Just tested in another DTM and it worked.

I used the “TIN interpolation” from QGIS because i can choose the from “structured line” and/or “break line”.

What are the alternatives?

Valter Albino <valteralbino@gmail.com> escreveu no dia segunda, 23/03/2020 à(s) 15:16:

···

Cumprimentos,
Valter Albino - Geógrafo Físico, M.Sc.

Modelação H&H / Riscos ambientais / OT&U

www.valteralbino.wixsite.com/hydrodynamics