GDAL讀取tiff文件/C++源碼

// gdal_geotiff.cpp : 定義控制台應用程序的入口點。
//

#include "stdafx.h"
#include "gdal_priv.h"
#include <iostream>
#include <iomanip>
#pragma comment(lib, "gdal_i.lib")
#define BYTE short      //方便數據類型的修改
using namespace std;

int main(int argc, char* argv[])
{
    int num_iamge_size=0;

    BYTE *pafScanblock1;  //開辟緩存區
    char *file_path_name="C:/Users/afa/Desktop/srtm_59_04/srtm_59_04.tif";
    GDALDataset *poDataset;   //GDAL數據集
    GDALAllRegister();  //注冊所有的驅動
    poDataset = (GDALDataset *) GDALOpen(file_path_name, GA_ReadOnly );
    if( poDataset == NULL )
    {
        cout<<"fail in open files!!!"<<endl;
        return 0;
    }

    //獲取圖像波段
    GDALRasterBand *poBand1;
    poBand1=poDataset->GetRasterBand(1);

    //獲取圖像的尺寸
    int nImgSizeX=poDataset->GetRasterXSize();
    int nImgSizeY=poDataset->GetRasterYSize();

    //獲取坐標變換系數
    double trans[6];
    CPLErr aaa=poDataset->GetGeoTransform(trans);

    //讀取圖像高程數據
    double Xgeo,Ygeo;
    pafScanblock1 = (BYTE *) CPLMalloc(sizeof(BYTE)*(nImgSizeX)*(nImgSizeY));
    poBand1->RasterIO( GF_Read, 0, 0,nImgSizeX,nImgSizeY,pafScanblock1,nImgSizeX,nImgSizeY,GDALDataType(poBand1->GetRasterDataType()),0, 0 );
    for(int i=0;i<(nImgSizeX-5990);i++)
    {
        for(int j=0;j<(nImgSizeY-5990);j++)
        {
            BYTE elevation=*pafScanblock1;
            Xgeo=trans[0]+i*trans[1]+j*trans[2];
            Ygeo=trans[3]+i*trans[4]+j*trans[5];
            num_iamge_size++;
            pafScanblock1++;
            cout<<setprecision(15)<<Xgeo<<"  "<<Ygeo<<"  "<<elevation<<endl;
        }
    }
    cout<<"總計："<<num_iamge_size<<endl;

    delete poDataset;
    return 0;
}