GDAL处理地理图像坐标计算
本文讲解如何使用GDAL处理地理图像时,通过使用行列号计算和转换成tiff图像的地理坐标:
tif中坐标计算的方法如下,其中Col表示该坐标点处图像的列号,ROW表示该坐标点处图像的行号。
比如图像左上角Col为0,ROW为0,图像右下角Col为图像宽度,ROW为图像高度。
Xgeo = GT(0) + Col*GT(1) + Row*GT(2)
Ygeo = GT(3) + Col*GT(4) + Row*GT(5)
#include <gdal_priv.h>
#include <ogr_spatialref.h>
#include <ogr_geometry.h>
#include <fstream>
#include <iostream>
#include <iomanip>
#define BYTE float //方便数据类型的修改
int main()
{
//tif文件读取
GDALAllRegister(); //注册所有的驱动
CPLSetConfigOption("GDAL_FILENAME_IS_UTF8", "NO");//以防中文名不能正常读取
int num_iamge_size = 0;
char *file_path_name = "../高程.tif";
GDALDataset *poDataset = (GDALDataset *)GDALOpen(file_path_name, GA_ReadOnly);//GDAL数据集
if (poDataset == NULL)
{
s