一个GDAL的读写数据例子（）应一个网友要求而发，离开GIS，RS行业快两年了，很多感触）

//

#include "stdafx.h"
#include "fangshibo.h"

#include <vector>


#include<math.h>
///gdal头文件
#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")
/


#ifdef _DEBUG
#define new DEBUG_NEW
#undef THIS_FILE
static char THIS_FILE[] = __FILE__;
#endif

/
// The one and only application object

CWinApp theApp;
using namespace std;

///
void  MaxNDVI(vector<CString> files,CString maxNDVfile)
{    
//    GDALAllRegister();
    vector<GDALDataset *> datasets;    
        //打开文件
        for(int i=0;i<files.size();i++)
        {    
            CString filename=files[i];            
            GDALDataset *tmpDataset = (GDALDataset *) GDALOpen( files[i], GA_ReadOnly);
            if(tmpDataset!=NULL)
            {                
                datasets.push_back(tmpDataset);                
            }
            else
            {