DEMの平滑化処理

[sasayama]

質問失礼します。
　
現在、点群から滑らかな等高線を生成できないかをQGIS version3.22.2で試行錯誤しています。
　
具体的には「QGISで滑らかな等高線を作成する」を参照しながら、277012点の点群からTIN内挿(不規則三角網)してDEMを生成し、r.mapcalc.simple(プロセシング > プロセシングツールボックス > GRASS > ラスタ > r.mapcalc.simple)でDEMを平滑化しようとしました。
　
しかし、TIN内挿(不規則三角網)で生成したDEMをr.mapcalc.simpleで処理すると「次のレイヤは正しく生成されませんでした」と表示され、何も生成されません。
　
TIN内挿とr.mapcalc.simpleのパラメータは添付写真を見て頂けると幸いです。
　
解決方法をご教示いただけると幸甚です。よろしくお願いします。

[mooya]

sasayamaさま

こんばんは

　不具合は環境によって異なりますので、OSとログをお示しください

　Windowdのmsi版でgrassが動かないという投稿が複数あり、私の環境でもエラーで動きません。

今の所対処法は、3.16のexe版を使うくらいです。

2022年1月15日土曜日 14:41:06 UTC+9 sasayama:

[sasayama]

Mooyaさま

　
ご回答いただき、ありがとうございます。ご返信、遅れて申し訳ございません。
　
Mooyaさまのご指摘通り、Windowsのmsi版を使っているので、上手くいかないのかもしれません。exe版を使ってみます。
　
OSはWindows 10 Proです。
　
ログは以下の通りになっています。
　
--------------------------

QGIS version: 3.20.3-Odense
QGIS code revision: 495fbaec
Qt version: 5.15.2
Pythonバージョン: 3.9.5
GDAL version: 3.3.1
GEOS version: 3.9.1-CAPI-1.14.2
PROJ version: Rel. 8.1.1, September 1st, 2021
PDAL version: 2.3.0 (git-version: eac774)
プロセシングアルゴリズム...
アルゴリズム 'r.mapcalc.simple' を開始しています...
入力パラメータ:
{ 'GRASS_RASTER_FORMAT_META' : '', 'GRASS_RASTER_FORMAT_OPT' : '', 'GRASS_REGION_CELLSIZE_PARAMETER' : 0, 'GRASS_REGION_PARAMETER' : None, 'a' : 'D:/3D/20210806_hotaka_e6_2019.0.0.1457-sokuryouki/20210806_hotaka_e6_contour/20210806_hotaka_e6_raster_raster_tin_from-subsampled-10cm_expect-trench/20210806_hotaka_e6_raster__raster_tin_from-subsampled-10cm_expect-trench.tif', 'b' : None, 'c' : None, 'd' : None, 'e' : None, 'expression' : '(A[-1,-1]+A[0,-1]+A[1,-1]+A[-1,0]+A[0,0]+A[1,0]+A[-1,1]+A[0,1]+A[1,1])/9', 'f' : None, 'output' : 'TEMPORARY_OUTPUT' }

g.proj -c wkt="C:/Users/owner/AppData/Local/Temp/processing_FlXHDk/b888bccd0a4140b9bcba6521ee5912d7/crs.prj"
r.in.gdal input="D:\3D\20210806_hotaka_e6_2019.0.0.1457-sokuryouki\20210806_hotaka_e6_contour\20210806_hotaka_e6_raster_raster_tin_from-subsampled-10cm_expect-trench\20210806_hotaka_e6_raster__raster_tin_from-subsampled-10cm_expect-trench.tif" band=1 output="rast_61e648dd166e82" --overwrite -o
g.region n=4023228.1375 s=4023101.5 e=755184.574 w=755086.82 res=0.09995296523524905
r.mapcalc.simple a=rast_61e648dd166e82 expression="(A[-1,-1]+A[0,-1]+A[1,-1]+A[-1,0]+A[0,0]+A[1,0]+A[-1,1]+A[0,1]+A[1,1])/9" output=outputc3d0c2be46ac47bbb83d68e7df14d312 --overwrite
g.region raster=outputc3d0c2be46ac47bbb83d68e7df14d312
r.out.gdal -t -m input="outputc3d0c2be46ac47bbb83d68e7df14d312" output="C:\Users\owner\AppData\Local\Temp\processing_FlXHDk\15532f93c7374aa1b1697c3a94177a4a\output.tif" format="GTiff" createopt="TFW=YES,COMPRESS=LZW" --overwrite
Starting GRASS GIS...
警告: マップセットの平行ロックはウィンドウズではサポートされていません。
Traceback (most recent call last):
File "C:\PROGRA~1/QGIS32~1.3/apps/qgis/./python/plugins\processing\algs\grass7\Grass7Algorithm.py", line 433, in processAlgorithm
Grass7Utils.executeGrass(self.commands, feedback, self.outputCommands)
File "C:\PROGRA~1/QGIS32~1.3/apps/qgis/./python/plugins\processing\algs\grass7\Grass7Utils.py", line 397, in executeGrass
for line in iter(proc.stdout.readline, ''):
UnicodeDecodeError: 'cp932' codec can't decode byte 0x84 in position 35: illegal multibyte sequence
Execution failed after 1.04 秒（seconds）

出力レイヤの読み込み
次のレイヤは正しく生成されませんでした。
• C:/Users/owner/AppData/Local/Temp/processing_FlXHDk/15532f93c7374aa1b1697c3a94177a4a/output.tif
QGISメインウィンドウの"ログメッセージパネル"をチェックすると、アルゴリズムの実行に関する詳細情報が表示されます。

--------------------------
　
お手数をお掛けしても申し訳ございませんが、ご指導のほどよろしくお願いします。
　
sasayama

2022年1月15日土曜日 22:23:44 UTC+9 mooya:

[sasayama]

mooyaさま
　
3.16のexe版で試してみたところ、うまくいきました。
　
ご指導いただき、ありがとうございます。
　
sasayama

2022年1月18日火曜日 14:16:02 UTC+9 sasayama: