GDAL使用方法VC+C#分析和总结.docxVIP

GDAL 栅格图像操作 GDAL 是一个操作各种栅格和矢量（由 ogr 这个库实现）地理数据格式的开源库。包括读取、写入、转换、处理各种栅格和矢量数据格式（有些特定的格式对一些操作如写入等不支持）。 即使不是进行地理遥感方面的应用研究，GDAL 也是一个非常有用的库，因为它可以支持大量我们常 见的图像数据，比如 jpg，gif 之类的。完整的格式清单可以到此链接查看 /formats_list.html。而且已经有包括 GoogleEarth 在内的很多软件都是在使用GDAL 作为后台库。 本文就以 VC 为开发平台介绍 GDAL 对栅格数据的操作方法。 Include 目录是开发中需要的头文件，lib 中是所需要的 lib 文件，在 VC8 中应当将其存放目录添加到目录列表中，选择菜单的“工具-选项-项目和解决方案-VC++目录”，分别在“包含文件”和“库文件”中将此两个目录添加进去。 在项目的属性页中，选择“配置属性-链接器-输入”，在“附加依赖项”中添加 gdal_i-vc8.lib 和gdal_id-vc8.lib 两个使用 GDAL 中需要的静态库文件，或者在程序中添加以下两行代码也可以。#pragma comment(lib, gdal_i-vc8.lib) #pragma comment(lib, gdal_id-vc8.lib) Bin 目录下的动态链接库文件应当放置于程序能够访问的位置，比如 windows\system32 中。此外，在程序中需要引入的头文件是 gdal_priv.h。 现在开始用 C++来对图像文件进行操作。 在打开文件之前需要首先注册所需要的驱动程序，一般来说我们可以默认注册所有支持的格式驱动， 所使用的函数是 GDALAllRegister()。然后就是打开文件操作。 这里要说一个数据集的概念，也就是所谓的 Dataset。在 GDAL 中可以说数据的核心就是 Dataset， 简单来说可以将 Dataset 就理解为图像文件，比如说一个 jpeg 格式的文件就是一个数据集，当然其他一些文件格式可能在一个数据集中包含多于一个文件，比如可能除了图像数据文件外还可能会有一些附加信息文件等。 在数据集下最重要组成部分就是所谓的波段 band，波段可多可少，比如一个 RGB 真彩色的图像就有 3 个波段，分别代表红色绿色和蓝色波段，如果是灰度图，那可能就只有一个波段，而很多遥感图像可 能就会多于 3 个波段。 除了波段外，数据集中还含有图像相关的坐标系投影信息，元数据信息等数据。 文件的打开使用的是 GDALOpen ( const char * pszFilename, GDALAccess eAccess )， pszFilename 是文件路径，eAccess 是访问权限，可以是 GA_ReadOnly 只读，也可以是 GA_Update 来对文件进行修改。比如我们以只读模式打开一个 tif 文件： GDALDataset *poDataset; //数据集对象指针GDALAllRegister();//注册驱动 poDataset = (GDALDataset *) GDALOpen( c:\\terra335h_EV_250_Aggr500_RefSB_b0.tif, GA_ReadOnly ); if( poDataset != NULL /*检查是否正常打开文件*/) { //do something } delete poDataset; //释放资源 在确认 poDataset 不是 NULL 的情况下就可以对图像数据集进行操作了。首先来看一下有关这个图像的基本信息，比如长宽和波段数。 coutRasterXSize:poDataset-GetRasterXSize()endl;//x 方向长度coutRasterYSize:poDataset-GetRasterYSize()endl;//y 方向长度coutRasterCount:poDataset-GetRasterCount()endl;//波段数量 我的这个文件输出结果如下： RasterXSize:5684 RasterYSize:3655 RasterCount:1 这是个 5684×3655 的灰度图。 有关的投影信息可以由 GetProjectionRef 函数获得。 coutpoDataset-GetProjectionRef()endl; 输出一个字符串，是 WKT 格式的坐标投影信息，形如： GEOGCS[WGS 84,DATUM[WGS_1984,SPHEROID[WGS 84,6378137,298.2572235630016,AUT HORITY[EPSG,7030]],AUTHORITY[EP