关于GDAL读取MODIS卫星数据

        由于modis卫星数据跟我们经常遇到的geotif数据组织方式不一样，读取的时候一定要特别注意。geotif数据，一般是一个文件，包含了多个波段的数据；而modis呢，一个文件包含了多各SUBDATASETSGDAL，每个SUBDATASETS又包含多个波段数据。另外默认编译的GDAL并不包含对MODIS数据支持，需要单独下载针对HDF4，HDF5的源码，再修改下make.opt文件，这时再编译GDAL，就支持modis数据的读写了。如果嫌麻烦就下载一个别人编译好了的来用，网上有很多基于GDAL的开源项目，搜一下就可以了。
       下面给一个例子，希望对大家用modis数据研究大尺度的地理现象有所帮助。modis数据毕竟免费嘛，呵呵……

#include <stdio.h>
#include <vector>
#include<string>
using namespace std;

#include<math.h>
#include "..//include//gdal.h"
#include "..//include//gdal_priv.h"
#include "..//include//ogr_srs_api.h"
#include "..//include//cpl_string.h"
#include "..//include//cpl_conv.h"
#pragma comment (lib,"..//lib//gdal_i.lib")

int main(int argc, char *argv[])
{