【QGIS】pythonコンソール上でsubprocessを用いてGDALコマンドを実行したい

Question

実現したいこと

QGIS上で「マスクレイヤによる切り抜き」をする際に表示されるgdalのコードをpythonのコードに読み替えて、これをsubprocessを用いてコンソール上で実行したいのですがうまくいきません(実行はされるが何も起きない)
。
subprocess.call内の表記が悪いと思うのですが、正しい表記をご教示いただけないでしょうか。

該当のソースコード

gdalwarp -of GTiff 
　　　　　-cutline マスクレイヤのパス
　　　　　-crop_to_cutline 入力レイヤのパス 出力先

#これを以下のpythonコードに読み替え

from osgeo import gdal
import subprocess

shp = "マスクレイヤのパス"
ori_tif = "入力レイヤのパス"
cre_tif = "出力先のパス"

subprocess.call(['gdalwarp','-of GTiff','-cutline'+ shp, '-crop_to_cutline'+ ori_tif,cre_tif])

補足情報（FW/ツールのバージョンなど）

QGIS: ver 3.6.1

Answer 1

gdalwarp使ったことありませんが、一般的なコマンドの使い方として

• コマンド文字列の間は空白文字で区切る
• パスに空白を含む場合は"で囲む

ことが必要です。
さらに、何か問題が起きターミナル上で実行した結果と比較したいときに備え、コマンドをひとつの文字列としてshell=Trueで呼び出す方がよいでしょう。
これであればコマンド文字列をコピペするだけでよいです。

cmd = 'gdalwarp -of GTiff -cutline "マスクレイヤのパス" -crop_to_cutline "入力レイヤのパス" "出力先"'
subprocess.call(cmd, shell=True)

Comments

[kurage]

ありがとうございます。ご教示いただいたコードで無事実行することができました。
ただ、今回のマスクレイヤには複数の地物（切り抜き用の枠）が入っており、上記コードを実行するとすべての地物（枠）でtifをクリップしてしまいます。このうち選択した一つの地物（枠）のみを使ってtifをクリッピングをするコマンドなどはあるのでしょうか。もしご存じでしたらご教示いただけないでしょうか。

[8524ba23]

すみません、GIS自体使ったことないのでちょっと分かりません。

[kurage]

ご返信ありがとうございます。
大変助かりました。