分幅输出影像瓦片

　　最近在做瓦片下载拼接的功能，但是如果将中国范围内的影像按照最高精度下载，拼接后的影像可能就太大了。假设0.5m分辨率，10km*10km的范围就是2w*2w像素。这个还是可以输出的，如果范围再大，50km*50km呢？就10w*10w像素了。用GDAL输出的Tiff格式上限是4G，Img无上限。但是太大总是很难搞。还是要分幅。将框选矩形范围内的的影像输出，同时保证输出图像的大小不至于太大，不超过2w*2w像素。

首先确定每一个分幅的大小假设为size=256*78≈20000。

基本算法流程：

int m = (int)((row2 - row1) / size);
int n = (int)((col2 - col1) / size);
//2000像素以拼接
for (int i = 0; i < m; i++)
{
　　for (int j = 0; j < n; j++)
　　{
　　　　计算左上角点行列号和坐标
　　　　写入size*size大小的影像
　　}
}
//右侧剩余
for (int i = 0; i <= row2 - row1 - m * size; i++)
{
　　for (int j = 0; j < n; j++)
　　{
　　　　计算左上角点行列号和坐标
　　　　写入（row2 - row1 - m * size）*size大小的影像
　　}
}

//底部剩余
for (int i = 0; i < m; i++)
{
　　for (int j = 0; j <= col2 - col1 - n * size; j++)
　　{
　　　　计算左上角点行列号和坐标
　　　　写入size*（col2 - col1 - n * size）大小的影像
　　}
}
//右下角剩余部分，特殊情况m=n=0，输出一个图幅
for (int ii = 0; ii <= row2 - row1 - m * size; ii++)
{
　　for (int jj = 0; jj <= col2 -col1 - n * size; jj++)
　　{
　　　　　计算左上角点行列号和坐标

　　　　　　　写入（row2 - row1 - m * size）*（col2 - col1 - n * size）大小的影像

　　} 
}