作者MAGICXX (逢甲阿法)
看板Python
標題[問題] subprocess.call with gdalwrap
時間Mon May 6 11:02:25 2013
大家好 我又來麻煩大家了
我現在想要用subprocess.call來用gdal裡面的gdalwrap來進行
shapefile 切 raster的方法(在Qgis 1.8.0 的python主控台下)
我先用raster→projections→wrap(reproject)產生下面command_line
-----------------------------------------------------------------------------
gdalwarp -t_srs "+proj=tmerc +lat_0=0 +lon_0=121 +k=0.9999 +x_0=250000 +y_0=0
+ellps=GRS80 +units=m +no_defs" -dstnodata 0 -q -cutline "C:/Program Files
(x86)/Quantum GIS Wroclaw/test/area.shp"
-dstalpha -of GTiff C:/Program Files
(x86)/Quantum GIS Wroclaw/test/wrapM_0.tif "C:/Program Files
(x86)/Quantum GIS Wroclaw/test/wrap"
------------------------------------------------------------------------------
其中 gdalwrap 是他的command名稱
下面是我的code
=============================================================================
#先加入gdal,subprocess.call兩個函式
import gdal
from subprocess import call
#將變數轉換成代號，再輸入會比較清楚
crs='"+proj=tmerc +lat_0=0 +lon_0=121 +k=0.9999 +x_0=250000 +y_0=0
+ellps=GRS80 +units=m +no_defs"'
mask='"C:/Program Files(x86)/Quantum GIS Wroclaw/test/area.shp"'
input='C:/Program Files(x86)/Quantum GIS Wroclaw/test/wrapM_0.tif'
output='"C:/Program Files(x86)/Quantum GIS Wroclaw/test/wrap"'
#將上面的command_line用代號代換之後執行call
call(gdalwarp, -t_srs crs -dstnodata 0 -q -cutline msak -dstalpha -of GTiff
input output)
=============================================================================
然後會出現下面的error
=============================================================================
File "<input>", line 1
call(gdalwarp, -t_srs crs -dstnodata 0 -q -cutline msak -dstalpha -of
GTiff input output)
^
SyntaxError: invalid syntax
=============================================================================
然後如果我的code用原本的command_line的話
也是會在同樣的地方發生錯誤
請問有用過的前輩知道大概哪邊錯了嗎?
fixed issue
=============================================================================
#將call函式改成這樣 就可以用了
call('gdalwarp -t_srs ' + crs + ' -dstnodata 0 -q -cutline ' + mask + '
-dstalpha -of GTIFF ' + input + ' ' + output, shell=True)
#或者這樣也可以
call(['gdalwarp', '-t_srs ' + crs, '-dstnodata 0', '-q', '-cutline ' + mask,
'-dstalpha', '-of GTIFF', input, output])
=============================================================================
--
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.112.63.202
1F：推 os653:你參數給的方法好詭異，舉例來說，ping ptt.cc應該寫成 05/06 19:30
2F：→ os653:subprocess.call(["ping", "ptt.cc"]) 05/06 19:32
3F：→ os653:你寫成 05/06 19:41
4F：→ os653:subprocess.call('ping', ptt.cc) 的樣子 05/06 19:44
5F：→ MAGICXX:喔喔~@@ 原來是我參數給錯方法了嗎!?囧 05/06 20:50
※ 編輯: MAGICXX 來自: 140.112.63.202 (05/08 11:03)