2019-11-08 使用shp文件批量裁剪GeoTif遥感图像

需求与目的

在批量将MODIS数据hdf转为GeoTif格式之后（具体操作可参考我的另一篇博客MODIS数据hdf转tif），因为是全球数据，而我只需要某一块区域的遥感影像，故而需要进行批量裁剪。（如下图所示）
数据按照月份文件夹存储，每个月份文件夹下有当月的tif文件。如下图所示：

运行环境

Windows 7、10（其它系统版本未尝试）；
ArcGIS10.5自带Python 2.7.12，因为要调用ArcGIS中的ArcPy模块，所以只能使用ArcGIS自带的Python 2.7.12版本；
ArcGIS需要在扩展模块中勾选 Spatial Analyst模块。如下图所示。

Python代码实现

因为是调用arcpy的代码，所以最重要的就一行(arcpy.gp.ExtractByMask_SA(file2_path, inMaskShp, outfile_path))。
另外，在创建shp文件时，请保证shp文件的地理坐标系与tif的地理坐标系保持一致。
完整代码如下：

结果

用ArcGIS打开看一下：

OK，搞定！