#include "gdal_priv.h"

#include<iostream>
#include <iostream>
#include <list>

#include "gdal_priv.h"   
#include "cpl_conv.h" //for CPLMalloc()   

using namespace std;
void ReadTiffFile(double lon, double lat)
{
	//注冊文件格式   
	GDALAllRegister();

	const char* pszFile = "./fromglc10v01_0_100.tif"; // 讀取tiff文件
	//const char* pszFile = "./01C_20210101-20220101.tif";//"C:\\Test.img";  //1.jpg
	//使用只讀方式打開圖像 
	GDALDataset* poDataset = (GDALDataset*)GDALOpen(pszFile, GA_ReadOnly);
	if (poDataset == NULL)
	{
		printf("File: %s不能打開！\n", pszFile);
		return;
	}

	//輸出圖像的格式信息   
	std::cout << "Driver:" << poDataset->GetDriver()->GetDescription() << "/"
		<< poDataset->GetDriver()->GetMetadataItem(GDAL_DMD_LONGNAME) << std::endl;

	//輸出圖像的大小和波段個數(shù)   
	std::cout << "Size is" << poDataset->GetRasterXSize() << "*"
		<< poDataset->GetRasterYSize() << "*" << poDataset->GetRasterCount() << std::endl;

	//輸出圖像的投影信息   
	if (poDataset->GetProjectionRef() != NULL)
	{
		std::cout << "Projectionis" << poDataset->GetProjectionRef() << std::endl;
	}

	GDALRasterBand* poBand;
	int            nBlockXSize, nBlockYSize;
	int            bGotMin, bGotMax;
	double         adfMinMax[2];

	//讀取第一個波段   
	poBand = poDataset->GetRasterBand(1);

	int   nXSize = poBand->GetXSize();
	int nYSize = poBand->GetYSize();

	//輸出圖像的坐標和分辨率信息   
	double adfGeoTransform[6];
	if (poDataset->GetGeoTransform(adfGeoTransform) == CE_None)
	{
		printf("Origin =(%.6f,%.6f)\n",
			adfGeoTransform[0], adfGeoTransform[3]);

		printf("PixelSize = (%.6f,%.6f)\n",
			adfGeoTransform[1], adfGeoTransform[5]);
	}


	CPLFree(poDataset);
	//關閉文件   
	GDALClose((GDALDatasetH)poDataset);
}

int main()
{
	double lon = 101.8615835729;
	double lat = 1.7389243989;
	ReadTiffFile(lon, lat);

	return 0;
}