使用C#版本的gdal库打开hdf文件

作者：朱金灿

来源：http://blog.csdn.net/clever101

最近应同事的请求帮忙研究下使用C#版的gdal库读取hdf文件，今天算是有一点成果，特地做一些记录。

首先是编译C#版的gdal库，请参考下面两篇文章：

1. GDAL源码剖析（三）之Swig编译和帮助文档生成

2. Windows下编译C#可用的GDAL附带netCDF格式支持库的一点拾遗

不过我没有遇到安全透明代码无法调用本机代码的问题，倒是遇到了OSGeo.GDAL.GdalPINVOKE”的类型初始值设定项引发异常的问题，解决办法是把gdal及其依赖的dll都拷贝到C#程序的输出目录下。

C#使用gdal读取影像数据函数如下：

Gdal.AllRegister();
Dataset ds = Gdal.Open("影像文件名", Access.GA_ReadOnly);

打开文件后，可通过ds来获取影像的信息。对于一般的影像文件比如tif等，gdal是完全可以通过这种方法来读取的，可是用这种方法来读取hdf文件是不大可能的，因为它实际数据在下一层。也就是子数据集。因此，如果想要获取hdf信息，需要操作每一个子数据集。

所以需要利用gdal自带的应用程序gdalinfo.exe来获取hdf文件的子数据集信息，比如风云三号的hdf文件子数据集信息如下：

SUBDATASET_1_NAME、SUBDATASET_2_NAME这些行的等号后面的就是子数据集名，比如要获取//Height子数据集的信息，代码如下：

       static void Main(string[] args)
        {
            OSGeo.GDAL.Gdal.AllRegister();
            OSGeo.GDAL.Gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8", "YES");
            Dataset ds = Gdal.Open("HDF5:d:\\TestData\\hdf\\FY3A_VIRRX_GBAL_L1_20130409_0315_1000M_MS.HDF://Height", Access.GA_ReadOnly);
        }

参考文献：

1. C#使用gdal读取HDF数据的正确方法