基于GDAL的遥感影像显示（C#版）

        说明：本文章转载自：http://blog.csdn.net/rsyaoxin/article/details/9220735

        接触GDAL有四五年多时间了，平时都是在C++下使用，最近需要在C#下调用GDAL，所以就开始学习了下，相比C++调用，C#下使用GDAL做图像处理的效率有点低，但是其简单易学，适合菜鸟上手，现把自己刚学到的心得跟大伙分享下，以遥感影像的显示为例。

1、  程序环境搭建

       GDAL是一个非常强悍的遥感数据格式解析库，支持多种遥感数据格式的读写，而且还有一些算法实现。然而，它只是一个开源库，并不是一个单独的软件，而是C++写的源代码，需要编译成动态链接库后才能为我们程序调用，如果需要在C#/.Net环境下调用，还需特别编译C#版的DLL。关于编译的方法网上很多教程，我就不介绍了，具体编译步骤可以参照民录大哥的博客：http://blog.csdn.net/liminlu0314/article/details/6937194

      如果很难编译成功，需要编译后的版本，可以留下邮箱，或者给我发邮件RSyaoxin@163.com。

      编译完成后，我们会得到9个DLL文件（以GDAL1.10版为例）：gdal110.dll、gdal_csharp.dll、gdal_wrap.dll、gdalconst_csharp.dll、gdalconst_wrap.dll、ogr_csharp.dll、ogr_wrap.dll、osr_csharp.dll、osr_wrap.dll。

注意：需要强调的一点是，如果编译GDAL时添加了额外的依赖库，需要将其动态链接库一并拷贝过来。比如，我编译GDAL库时就添加了HDF4、HDF5、JPEG2000、NetCDF、Proj4等文件格式支持，那么调用的时候就需要将hd425m.dll等拷贝到同上面几个DLL一起。

        我们首先建立一个Windows窗体应用程序，将上面得到的9个DLL文件拷贝到项目所在的文件夹，然后在解决方案面板下找到引用项，右键添加引用，选择gdal_csharp.dll将它添加进来，这样就可以在C#调用GDAL函数了。

         再次，初始化GDAL环境。具体做法是：

        using  OSGeo.GDAL;

        将上面这一句添加到namespace这一行前面。

        然后在构造函数或者Form_Load函数中添加下面两行：

        OSGeo.GDAL.Gdal.AllRegister();

        OSGeo.GDAL.Gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8","YES");

        说明：第一行是注册所有的格式驱动，第二行是支持中文路径和名称，由于GDAL默认不支持中文路径，所以在编译的时候会修改源码让它支持中文路径，C++代码就可以直接支持了，但是会发现C#版还是会不支持，所以需要加上第二句。

2、  关键代码

         C#下图像的显示方法有多种，最简单的方法就是构建位图。我们可以把GDAL下的Dataset转换为Bitmap供C#调用。在C#里面调用GDAL读取栅格数据的主要函数是ReadRaster，它相当于C++下的RasterIO函数，其有多重形式：

public CPLErr ReadRaster(int xOff, int yOff, int xSize, int ySize, byte[] buffer, 
    int buf_xSize, int buf_ySize, int pixelSpace, int lineSpace)
public CPLErr WriteRaster(int xOff, int yOff, int xSize, int ySize, byte[] buffer, 
    int buf_xSize, int buf_ySize, int pixelSpace, int lineSpace)
public CPLErr ReadRaster(int xOff, int yOff, int xSize, int ySize, short[] buffer, 
    int buf_xSize, int buf_ySize, int pixelSpace, int lineSpace)
public CPLErr WriteRaster(int xOff, int yOff, int xSize, int ySize, short[] buffer, 
    int buf_xSize, int buf_ySize, int pixelSpace, int lineSpace)
public CPLErr ReadRaster(int xOff, int yOff, int xSize, int ySize, int[] buffer, 
    int buf_xSize, int buf_ySize, int pixelSpace, int lineSpace)
public CPLErr WriteRaster(int xOff, int yOff, int xSize, int ySize, int[] buffer, 
    int buf_xSize, int buf_ySize, int pixelSpace, int lineSpace)
public CPLErr ReadRaster(int xOff, int yOff, int xSize, int ySize, float[] buffer, 
    int buf_xSize, int buf_ySize, int pixelSpace, int lineSpace)
public CPLErr WriteRaster(int xOff, int yOff, int xSize, int ySize, float[] buffer, 
    int buf_xSize, int buf_ySize, int pixelSpace, int lineSpace)
public CPLErr ReadRaster(int xOff, int yOff, int xSize, int ySize, double[] buffer, 
    int buf_xSize, int buf_ySize, int pixelSpace, int lineSpace)
public CPLErr WriteRaster(int xOff, int yOff, int xSize, int ySize, double[] buffer, 
    int buf_xSize, int buf_ySize, int pixelSpace, int lineSpace)
public CPLErr ReadRaster(int xOff, int yOff, int xSize, int ySize, IntPtr buffer, i
    nt buf_xSize, int buf_ySize, DataType buf_type, int pixelSpace, int lineSpace)
public CPLErr WriteRaster(int xOff, int yOff, int xSize, int ySize, IntPtr buffer, 
    int buf_xSize, int buf_ySize, DataType buf_type, int pixelSpace, int lineSpace)

       关键代码如下：

3、  主函数调用

  得到Bitmap，我们就可以在程序中调用它了。我们可以在窗体上加一个PictureBox控件来显示图像，其name设为pictureBox1。主要调用代码如下：

4、  结果展示

    这样我们就可以很容易地显示遥感影像了。