GDAL仿射变换
一、仿射变换(Affine GeoTransform)有六个参数，描述的是栅格行列号和地理坐标之间的关系。
在GDAL中可以通过函数GDALDataset::GetGeoTransform()函数获得，存放在一个double类型的数组里面。
每个元素下标都代表着不同的含义：
- 0：图像左上角的X坐标；
- 1：图像东西方向分辨率；
- 2：旋转角度，如果图像北方朝上，该值为0；
- 3：图像左上角的Y坐标；
- 4：旋转角度，如果图像北方朝上，该值为0；
- 5：图像南北方向分辨率；
二、基于遥感影像具有的特点，这六个参数可以描述遥感影像起点坐标、分辨率、旋转角度。
以C++实现代码为例：
GDALAllRegister();
GDALDataset *pSrcDS = (GDALDataset*)GDALOpen(“源数据的路径”, GA_ReadOnly);
double GeoTransform[6] = {0};//初始化
pSrcDS->GetGeoTransform(GeoTransform);
//这样就得到了读取数据的仿射变换参数，并且已经放入了GeoTransform数组中。
GeoTransform[0]和GeoTransform[3]代表图像左上角坐标（地理坐标）
GeoTransform[1]和GeoTransform[5];代表影像分辨率（东西方向为正，南北方向为负）
GeoTransform[2]和GeoTransform[4];一般值都是0
三、得到这六个参数之后就可以进行图像行列号和地理坐标之间的变换，变换关系如下：
Xgeo = GeoTransform[0] + Xpixel*GeoTransform[1] + Ypixel*GeoTransform[2]
Ygeo = GeoTransform[3] + Xpixel*GeoTransform[4] + Ypixel*GeoTransform[5]
这样就能得到图像任一点的地理坐标了。
爱被打了一巴掌
拼搏现实的明天。
柔情只为你懂
左手的ㄟ右手
还没有评论，来说两句吧...