QGISで点群データからTifを一括作成したいのでpythonコードを組んでみる

本記事は以下の記事を参考にしています。
・三次元点群データをQGISでひり出してみた
処理の大まかな流れとしてはtxt→las→dtm→tifの順で作成します。

事前準備

QGIS 3.28.6
LAStoolsのセットアップ
FUSIONのセットアップ
※ LAStoolsとFUSIONの本体DL後のインストール先に関して
　こちらの環境だと何故かCドライブ直下に置かないと動きませんでした
もし動かないのであれば参考にされてください。
(理由は分かりません、だれか教えてください)

点群データ(txt)のフォーマット

1	-14000	69712.16	59.38	2
2	-14000	69714.98	60.93	1
3	-14000	69715.75	61.02	1
4	-14000	69728.65	53.12	1
5	-14000	69744.88	50.03	1
6	-14000	69706.35	54.81	2
7	-14000	69747.85	48	2
8	-14000	69748.55	49.81	1
9	-14000	69719.77	54.95	2
10	-14000	69741.7	51.9	1
11	-14000	69703.45	62.23	1
12	-14000	69720.27	49.39	5

上記はNo,X,Y,Z,flgのフォーマットとなります。

QGISのpythonコンソールからプログラムの実行

pythonコンソールを出す（右上側にある赤で囲まれたボタンを押下します）
エディタを表示(黄色枠) 表示されない場合は左下の赤枠のボタンを押下

3.以下のコードをエディタに張り付けます
target_pathを作業対象フォルダの絶対パスで書き換えてください
target_pathで指定した作業対象フォルダの中に不要なtxtは格納しないでください

createTif.py

import processing
import os

"""
in: target_path  txtが入ったフォルダを指定
　　　　　　　     ※target_pathで指定したフォルダには作業対象の
　　　　　　　　　　  txt以外入れないでください
out: target_pathにtifとlasとdtmが作成される
"""
target_path = r"C:\Users\user\hogehoge"

for filename in os.listdir(target_path):
    if filename.endswith('.txt'):
        input_txt_path = os.path.join(target_path, filename)
        out_name = os.path.splitext(filename)[0] + ".las"
        out_las_path = os.path.join(target_path, out_name)
        # lasデータの作成
        try:
            processing.run("LAStools:Txt2Las", {'VERBOSE':False,
                                                'CPU64':True,
                                                'GUI':False,
                                                'INPUT_GENERIC':input_txt_path,
                                                'PARSE':'sxyz',
                                                # ※1
                                                'SKIP':0,
                                                'SCALE_FACTOR_XY':0.01,
                                                'SCALE_FACTOR_Z':0.01,
                                                'PROJECTION':1,
                                                # ※3
                                                'EPSG_CODE':6670,
                                                'UTM':0,
                                                'SP':0,
                                                'OUTPUT_LASLAZ':out_las_path,
                                                'ADDITIONAL_OPTIONS':''
                                                })
        except Exception:
        # ※2
            print("it's a decoding error, so ignore it.")
        
        # dtmデータの作成
        out_name = os.path.splitext(filename)[0] + ".dtm"
        out_dtm_path = os.path.join(target_path, out_name)
        processing.run("fusion:gridsurfacecreate", {'INPUT':out_las_path,
                                                    'CELLSIZE':0.1,
                                                    'XYUNITS':0,
                                                    'ZUNITS':0,
                                                    'VERSION64':True,
                                                    'OUTPUT_DTM':out_dtm_path,
                                                    'SPIKE':'',
                                                    'MEDIAN':'',
                                                    'SMOOTH':'',
                                                    'MINIMUM':'',
                                                    'CLASS':'',
                                                    'SLOPE':'',
                                                    'ADVANCED_MODIFIERS':''}
                                                    )
        # tifデータの作成
        out_name = os.path.splitext(filename)[0] + ".tif"
        out_tif_path = os.path.join(target_path, out_name)
        processing.run("fusion:dtm2ascii", {'INPUT':out_dtm_path,
                                            'CSV':False,
                                            'RASTER':False,
                                            'MULTIPLIER':None,
                                            'OUTPUT':out_tif_path
                                            })
print("finish")

※1 txtのフォーマットが変わった場合は'PARSE'のオプション変更が必要
　　現在のオプションの"sxyz"のsはskip対象列の事、txtが先頭からxyzの順で
　　並ぶのであればPARSEのオプションはxyzのみで良い
※2 LAStoolsのTxt2Lasでlasを作成する際に何故かUnicodeDecodeErrorになります。
　　Errorは出てもlasはしっかり吐き出してくれているので問題視はしてないです。
　　UnicodeDecodeErrorを指定してcatchしようとしても何故か
　　catch出来ないので無理やりExceptionにしてます。
　　つよつよエンジニアの方対策教えてください。
※3 座標系は'EPSG_CODE':6670,ここをいじくりまわしてください

4.赤枠の実行ボタン押下

めんどくさかったので終了時のダイアログは用意してないです。
コンソールにfinishと出たらokです。

例：作成したTIF