GDAL遥感影像5000标准分幅

遥感影像分幅是很常见的工作，通常也是很简单的工作，也就是影像裁剪。

但是目前我所要做的遥感影像分幅要求比较苛刻，影像分辨率为0.5米，图幅为1:5000标准图幅，要求外扩100个像素，分幅后的影像左上角像素中心坐标是图幅外扩100个像素矢量的角点坐标。

在常用的软件中裁剪功能都是影像像素的左上角为坐标，如ArcGIS，PCI；erdas可以满足该分幅要求，但是可行性不高，图幅多麻烦。

因此，我使用gdal做了一个分幅功能，简单介绍原理：图幅外扩100个像素后的矢量与影像求交集，创建交集范围大小的栅格数据，用rasterIO搞定。