对于很多情况下，批量栅格数据的掩膜或者裁剪都是很不方便的，这种问题我也遇到过。

目前，app_store中只有批量的矢量裁剪栅格，对于栅格掩膜栅格，网上也没有很好的代码

所以如何利用envi+idl实现批量掩膜或者裁剪呢？

1 栅格掩膜介绍

简单来说，就是用一个栅格裁剪另一个栅格，分别叫做被裁剪栅格和裁剪栅格，那么最终的结果大小和裁剪栅格一样。

2 举例

图中为同一区域各种水体指数图，那么在我研究特定区域时，需要对其进行裁剪，原先的做法就是一幅一幅裁剪，但是这个效率太低下了。我当时裁剪了好久，所以，将其批量进行处理就很方便了。

3 裁剪代码envi+idl 二次开发

;+

; :Author: 15845

;-date:2019/3/7

;********************栅格文件批量裁剪*******************************************

pro batch_raster

compile_opt idl2

envi,/restore_base_save_files

envi_batch_init

;***********裁剪文件的打开**********

cd,"d:\"

file=dialog_pickfile(title="打开裁剪文件")

envi_open_file,file,r_fid=fid

envi_file_query,fid,ns=ns,nl=nl,nb=nb,dims=dims,data_type=type,interleave=interleave

map_info=envi_get_map_info(fid=fid)

base_raster=envi.openraster(file)

;获取的裁剪文件的左上角的地理坐标

base_map_coordinate=base_raster.spatialref.tie_point_map

x=base_map_coordinate[0]

y=base_map_coordinate[1]

; ****批量裁剪文件的打开*********

fn_subset=dialog_pickfile(title="打开栅格数据文件",/directory) ;被裁剪的文件

fn_name=file_search(fn_subset,'*.{dat,tif}',count=count_subset)

fn_basename=file_basename(fn_name) ;获取被裁减文件的名称

for i=0,count_subset-1 do begin

temp=envi.openraster(fn_name[i])

temp.spatialref.convertmaptofile,x,y,x1,y1 ;x1,y1为裁剪文件在被裁剪文件中的左上角位置

x_base=fix(x1)

y_base=fix(y1)

subset=envisubsetraster(temp,sub_rect=[x_base,y_base,x_base+ns-1,y_base+nl-1])

data="subset_"+fn_basename[i]

subset.export,data,"envi"

endfor

end说明:同时打开代码，运行处理。

### 裁剪结果：

4 总结这个栅格数据批量掩膜使用的是envi+idl 新的开发模式，即面向对象思想的编写。代码量减少很多。

因知乎排版问题，同步我的Github