本节书摘来自异步社区出版社《GDAL源码剖析与开发指南》一书中的第1章，第1.9节，作者：李民录 更多章节内容可以访问云栖社区“异步社区”公众号查看。

1.9 简单的调用

GDAL源码剖析与开发指南
关于GDAL的使用，关键是要熟悉GDAL的组织结构、类以及类的函数等。最常用的就是动态库的GDAL，当然我们也可以使用静态库，这里只是简单地介绍使用动态GDAL库来做开发。

下面我们简单地使用C++、C#、Python和Java 4种语言分别进行示例，来说明如何使用GDAL进行开发。

1.9.1 C ++使用GDAL
首先我们打开Visual Studio2008，新建一个Win32的控制台工程（名字叫GDALCppTest），然后在工程的属性对话框中，找到“配置属性”→“C/C++”→“常规”，在右侧的“附加包含目录”中将GDAL的include文件夹路径填写到输入框，如图1-14所示。

我们继续在属性对话框中找到“配置属性”→“链接器”→“常规”，在右侧的“附加库目录”中将GDAL的lib文件夹路径填写完整，如图1-15所示。

最后在“配置属性”→“链接器”→“输入”，在右侧的“附加依赖项”中将gdal_i.lib填写完整。然后点击“确定”按钮即可。至此，使用GDAL的环境全部搭建完成，剩下的就是在代码中使用GDAL了。

我们将以下代码复制到刚才新建工程中的cpp文件中，保存后编译。正常情况下系统会提示生成成功，然后我们将GDAL编译出来的gdal110.dll及其依赖的dll文件复制到编译出来的exe所在目录，最后运行，会在控制台上输出图像的一些信息。

#include "gdal\_priv.h"
int main()
{//注册文件格式GDALAllRegister();const char* pszFile = "C:\\Test.img";//使用只读方式打开图像GDALDataset *poDataset = (GDALDataset*) GDALOpen( pszFile, GA\_ReadOnly );if( poDataset == NULL ){printf( "File: %s不能打开！\n",pszFile);return 0;}//输出图像的格式信息printf( "Driver:%s/%s\n",poDataset->GetDriver()->GetDescription(),poDataset->GetDriver()->GetMetadataItem( GDAL\_DMD\_LONGNAME) );//输出图像的大小和波段个数printf( "Size is%dx%dx%d\n",poDataset->GetRasterXSize(),poDataset->GetRasterYSize(),poDataset->GetRasterCount());//输出图像的投影信息if( poDataset->GetProjectionRef() != NULL )printf( "Projectionis `%s'\n", poDataset->GetProjectionRef() );//输出图像的坐标和分辨率信息double adfGeoTransform[6];if( poDataset->GetGeoTransform( adfGeoTransform) == CE\_None ){printf( "Origin =(%.6f,%.6f)\n",adfGeoTransform[0], adfGeoTransform[3]);printf( "PixelSize = (%.6f,%.6f)\n",adfGeoTransform[1], adfGeoTransform[5]);}//读取第一个波段GDALRasterBand *poBand = poDataset->GetRasterBand( 1 );//获取该波段的最大值最小值，如果获取失败，则进行统计int            bGotMin, bGotMax;double         adfMinMax[2];adfMinMax[0] = poBand->GetMinimum( &bGotMin);adfMinMax[1] = poBand->GetMaximum( &bGotMax);if( ! (bGotMin&& bGotMax) )GDALComputeRasterMinMax((GDALRasterBandH)poBand, TRUE, adfMinMax);printf( "Min=%.3fd,Max=%.3f\n", adfMinMax[0], adfMinMax[1] );int nXSize = poBand->GetXSize();float *pafScanline = new float[nXSize];//读取图像的第一行数据poBand->RasterIO(GF\_Read, 0, 0, nXSize,1, pafScanline, nXSize,1, GDT\_Float32, 0, 0 );delete []pafScanline;//关闭文件GDALClose((GDALDatasetH)poDataset);
}

执行的结果如图1-16所示。

1.9.2 C#使用GDAL
我们首先打开Visual Studio 2008，新建一个Visual C#的控制台应用程序（名字叫GDALCSTest），然后在工程的“引用”点击右键，打开“添加引用”对话框，找到浏览标签页，将GDAL以_cshrp.dll结尾的四个dll添加到工程中，如图1-17所示。

接下来我们将下面的代码复制到刚才新建工程中的Program.cs文件中，保存后编译。正常情况下系统会提示生成成功，然后我们将GDAL编译出来的gdal110.dll及其依赖的dll文件复制到编译出来的exe所在目录，最后运行，会在控制台上输出图像的一些信息。

需要值得注意的是，在64位的系统下，如果GDAL编译的时候是32位的，那么选择编译C#工程的时候，我们要使用X86编译才能正常运行，否则会报错。

using System;
using OSGeo.GDAL;
using OSGeo.OSR;
namespace GDALCSTest
{class Program{static void Main(string[] args){Console.WriteLine("");try{//注册文件格式Gdal.AllRegister();string strFile = "F:\\Work\\GDAL源码剖析\\GDALTest\\Test.tif";//使用只读方式打开图像Dataset ds = Gdal.Open(strFile, Access.GA\_ReadOnly);if (ds == null){Console.WriteLine("不能打开：" + strFile);System.Environment.Exit(-1);}Driver drv = ds.GetDriver();if (drv == null){Console.WriteLine("Can't get driver.");System.Environment.Exit(-1);}//输出图像的格式信息Console.WriteLine("Dirver :" + drv.GetDescription());Console.WriteLine("Using driver " + drv.LongName);//输出图像的大小和波段个数Console.WriteLine("  RasterCount: " + ds.RasterCount);Console.WriteLine("  RasterSize (" + ds.RasterXSize + "," + ds.RasterYSize + ")");//输出图像的投影信息Console.WriteLine("  Projection: " + ds.GetProjectionRef());//输出图像的坐标和分辨率信息double[] adfGeoTransform = new double[6];ds.GetGeoTransform(adfGeoTransform);Console.WriteLine("Origin =(" + adfGeoTransform[0] + "," + adfGeoTransform[3] + ")");Console.WriteLine("PixelSize = (" + adfGeoTransform[1] + "," + adfGeoTransform[5] + ")");Console.WriteLine("");//读取第一个波段Band band = ds.GetRasterBand(1);double val;int hasval;band.GetMinimum(out val, out hasval);if (hasval != 0)Console.WriteLine("   Minimum: " + val.ToString());band.GetMaximum(out val, out hasval);if (hasval != 0)Console.WriteLine("   Maximum: " + val.ToString());int nXSize = ds.RasterXSize;float[] pafScanLine;pafScanLine = new float[nXSize];band.ReadRaster(0,0, nXSize, 1, pafScanLine, nXSize, 1, 0, 0);}catch (Exception e){Console.WriteLine("Application error: " + e.Message);}}}
}

1.9.3 Python使用GDAL
在使用Python版本GDAL之前，我们首先要确保GDAL的Python版本安装成功，否则无法按照下面的步骤进行安装。

按照1.7.3小节的python编译步骤进行编译，之后我们将GDAL所有的dll文件复制到python版本的GDAL安装目录中，如图1-18所示。

接下来我们使用记事本等工具，将以下代码复制进去，保存为GDALPythonTest.py。

from osgeo import gdal
from osgeo.gdalconst import *
gdal.AllRegister()
#使用只读方式打开图像
ds = gdal.Open("F:\Work\GDAL源码剖析\GDALTest\Test.tif",GA\_ReadOnly)
#输出图像的格式信息
print ds.GetDriver().ShortName
#输出图像的大小和波段个数
print ds.RasterXSize
print ds.RasterYSize
print ds.RasterCount
#输出图像的投影信息
print ds.GetProjectionRef()
#输出图像的坐标和分辨率信息
adfGeoTransform = ds.GetGeoTransform()
print adfGeoTransform[0]
print adfGeoTransform[3]
print adfGeoTransform[1]
print adfGeoTransform[5]
band = ds.GetRasterBand(1)
#读取图像的第一行数据
data = band.ReadRaster( 0, 0, ds.RasterXSize, 1,
ds.RasterXSize, 1, band.DataType )

我们在开始菜单中的Python目录中打开“IDLE（Python GUI）”，在弹出的程序中，使用文件菜单中的打开命令打开刚才保存的GDALPythonTest.py文件，如图1-19所示，然后按“F5”运行，即可输出影像的一些信息，如图1-20所示。

**
1.9.4 Java使用GDAL**
Java使用GDAL，这里使用的IDE是Eclipse。我们使用Eclipse新建一个工程，名字叫GDALJavaTest，然后在这个工程中新建一个名为GDALTestApp的类，打开GDALTestApp.java，将里面的代码修改为以下代码内容，然后保存，如图1-21所示。

import org.gdal.gdal.Band;
import org.gdal.gdal.Dataset;
import org.gdal.gdal.Driver;
import org.gdal.gdal.gdal;
import org.gdal.gdalconst.gdalconstConstants;
public class GDALTestApp {/*** @param args*/public static void main(String[] args) {gdal.AllRegister();String fileName\_tif = "F:\\Work\\GDAL源码剖析\\GDALTest\\Test.tif";  //使用只读方式打开图像Dataset hDataset = gdal.Open(fileName\_tif, gdalconstConstants. GA\_ReadOnly);if (hDataset == null){System.err.println("GDALOpen failed - " + gdal.GetLastErrorNo());System.err.println(gdal.GetLastErrorMsg());System.exit(1);}//输出图像的格式信息Driver hDriver = hDataset.GetDriver();System.out.println("Driver: " + hDriver.getShortName() + "/" + hDriver.getLongName());//输出图像的大小和波段个数int iXSize = hDataset.getRasterXSize();int iYSize = hDataset.getRasterYSize();int iBandCount = hDataset.getRasterCount();System.out.println("Size is " + iXSize + ", " + iYSize + "," + iBandCount);Band band = hDataset.GetRasterBand(1);int itype = band.GetRasterDataType();System.out.println(itype + "\n");//读取一行数据int buf[] = new int[iXSize]; band.ReadRaster(0, 0, iXSize, 1, buf); hDataset.delete();// 可选gdal.GDALDestroyDriverManager();}
}

我们将之前编译Java版本GDAL时生成的五个文件（dll和jar文件）复制到工程所在文件夹中，然后点击“工程”右键，选择“Build Path->configure build path”，弹出“属性”对话框，然后切换到“libraries”标签，点击“Add Jars…”，选择gdal.jar和名字中有ini的dll文件，添加到lib中，如图1-22所示，然后确定，编译运行即可。

**
1.9.5 调试GDAL源代码**
很多时候我们需要查看GDAL的源代码以了解它的内部实现，为了更好地理解源代码，可以通过调试查看源代码中变量的内容。

调试GDAL的源代码，需要GDAL的debug版本以及编译GDAL的时候生成的pdb等调试文件。此外我们也可以把GDAL的源代码加入到自己的工程中，但是这样太费时费力。

这里介绍一个很简单的方法：首先编译一下GDAL的debug版本（具体步骤可以参考1.6节的内容），将编译生成的文件（主要有gdal110.dll，gdal_i.exp，gdal_i.lib，gdal.lib，gdal110.pdb，gdal110.ilk，gdal110.exp等）复制到自己工程的生成目录中，然后调试自己的程序，在执行到GDALOpen函数（或者其他GDAL的函数）时按“F11”键，就可以进入到GDAL的源代码中进行调试。

对于使用C#语言来说，也可以跨语言来对GDAL进行调试。首先我们要确保GDAL库编译时使用debug版本，同时编译的C#版本也是debug版本。在执行调试时我们需要将GDAL的C#版本中八个dll文件及其对应的pdf等文件全部复制到测试程序输出目录，然后在测试目录的属性中勾选“启用非托管代码调试”，最后启动，按“F11”键即可进入GDAL源码。