功能:实现单个矢量数据(shp文件)对多幅遥感影像的裁剪
pro batch_subset_raster_task
compile_opt idl2
;compole_opt是对idl编译规则的修改,关键字有defint32,strictarr,idl2等
;defint32把默认的idl整型数据16位改为32位
;strictarr强迫数组元素用中括号,不用小括号
;idl2则为上述二者的并。
e=envi(/headless)
;启动envi
start = systime(1)
image_dir = DIALOG_PICKFILE(PATH='D:\IDL_workspace', /DIRECTORY, $
TITLE="选择包含裁剪数据的文件夹")
image_files=file_search(image_dir,'*.tif',count=numfiles) ;根据相应的文件格式修改过滤条件
vectorFile = Dialog_pickfile(PATH='D:\IDL_workspace',title='选择裁剪矢量文件') ;选择矢量数据
print,image_files
for i=0,numfiles-1 do begin
image_file=image_files[i]
;print,image_file
if strlen(image_file) eq 0 then return
raster = e.OpenRaster(image_file)
; Select input vector data
vector = e.OpenVector(vectorFile)