gdal of c++ day01

1.基础指令值

CPLErr GDALRasterBand::RasterIO (   GDALRWFlag eRWFlag, // 读入还是写入
int     nXOff, // 起始坐标(左上角x值）
int     nYOff, // 起始坐标(左上角y值）
int     nXSize, // （读入/写入）图像的宽
int     nYSize, // （读入/写入）图像的高
void * pData, // 指向（读入/写入）的指针
int     nBufXSize, // 缓冲区宽（真正显示的宽）
int     nBufYSize, // 缓冲区高（真正显示的高）
GDALDataType    eBufType, // （读入/写入）的数据的类型
int     nPixelSpace, // 像素间隔，默认0
int     nLineSpace  // 像素行间隔，默认0
)

2.具体实现
实例1

#include<iostream>
#include<algorithm>
#include<vector>
#include "gdal_priv.h"

BYTE * loadRasterData()
{		
		const char * fileName = "D:\\Programme\\Viual Studio\\C++\\learning\\point_to_line\\Project1_test\\images\\grid.tiff";
		GDALAllRegister(); // 注册驱动
		// 支持中文路径
		CPLSetConfigOption("GDAL_FILENAME_IS_UTF8", "NO");
		// 读取数据
		GDALDataset * pDataset = (GDALDataset *)GDALOpen(fileName, GA_ReadOnly);

		// 判断是否读取成功
		if (NULL == pDataset)
		{
			cout << "未能成功读取！" << endl;
			return 0;
		}

	// 图像波段信息
	GDALRasterBand * pRasterBand;
	pRasterBand = pDataset->GetRasterBand(1); // 第一个波段指针
	CPLString DataType = GDALGetDataTypeName(pRasterBand->GetRasterDataType()); // 获取数据类型名字
	int imgX = pRasterBand->GetXSize();
	int imgY = pRasterBand->GetYSize();
	cout << "x,y分别为：" << imgX << "  " << imgY << endl;
	cout << "数据类型：" << DataType << endl;
	cout << "-------------------------------------------------------" << endl;

	// 读取数据
	int size_m = sizeof(int)*imgX*imgY;
	cout << size_m << endl;
	BYTE * mat = (BYTE *)CPLMalloc(size_m); // 分配内存
	pRasterBand->RasterIO(GF_Read, 0, 0, imgX, imgY, mat, imgX, imgY, GDT_Byte, 0, 0);

	// 打印数组
	for (int i = 0; i < imgX; i++)
	{
		for (int j = 0; j < imgY; j++)
		{
			cout << (int) mat[i*imgX+j] << " ";
		}
	}

	return mat; // 返回数组的首地址
}


int main()
{
	// 加载数据
	BYTE *mat = loadRasterData();

}

实例2

#include<iostream>
#include "gdal_priv.h"

using namespace std;

// 读取栅格数据
void loadRasterData() {

	const char * fileName = "D:\\Programme\\Viual Studio\\C++\\learning\\point_to_line\\Project1_test\\images\\grid.tiff";

	GDALAllRegister(); // 注册驱动
	// 支持中文路径
	CPLSetConfigOption("GDAL_FILENAME_IS_UTF8", "NO");
	// 读取数据
	GDALDataset * pDataset = (GDALDataset *)GDALOpen(fileName, GA_ReadOnly);

	// 判断是否读取成功
	if (NULL == pDataset)
	{
		cout << "未能成功读取！" << endl;
		return;
	}


	// 输出图像格式信息
	CPLString strDrive = pDataset->GetDriver()->GetDescription();
	CPLString strInfo = pDataset->GetDriver()->GetMetadataItem(GDAL_DMD_LONGNAME);
	cout << "文件格式：" << strDrive << " " << strInfo << endl;

	// 获取波段，长宽和波段数
	int x = pDataset->GetRasterXSize();
	int y = pDataset->GetRasterYSize();
	int band = pDataset->GetRasterCount(); // 高
	cout << "影像x: " << x << endl;
	cout << "影像y: " << y << endl;
	cout << "波段数：" << band << endl;
	cout << "-------------------------------------------------------" << endl;

	// 输出投影方式
	CPLString  strProjection = pDataset->GetProjectionRef();
	if (strProjection != "") // 如果没有投影信息，返回一个空字符串，不是NULL;
	{
		cout << "投影方式：" << strProjection << endl;
		cout << "-------------------------------------------------------" << endl;

	}

	//输出几何信息
	double GeoTransform[6];
	if (pDataset->GetGeoTransform(GeoTransform) == CE_None) // CE_None 表示成功获取
	{
		cout << "几何信息为：" << endl;
		for (int i = 0; i < 6; i++)
		{
			cout << GeoTransform[i] << " ";
		}
		cout << endl;
		cout << "-------------------------------------------------------" << endl;

	}

	// 图像波段信息
	GDALRasterBand * pRasterBand;
	pRasterBand = pDataset->GetRasterBand(1); // 第一个波段指针
	int nxBlock, nyBlock;
	pRasterBand->GetBlockSize(&nxBlock, &nyBlock); // TODO
	CPLString DataType = GDALGetDataTypeName(pRasterBand->GetRasterDataType()); // 获取数据类型名字
	CPLString ColorType = GDALGetColorInterpretationName(pRasterBand->GetColorInterpretation()); // 获取颜色类型名称
	cout << "分块大小为：" << nxBlock << " / " << nyBlock << endl;
	cout << "数据类型：" << DataType << endl;
	cout << "颜色类型：" << ColorType << endl;
	cout << "-------------------------------------------------------" << endl;

	// 统计像素的最值
	int arrMaxMin[2];
	arrMaxMin[0] = pRasterBand->GetMaximum();
	arrMaxMin[1] = pRasterBand->GetMinimum();
	cout << "像素的最大值：" << arrMaxMin[0] << endl;
	cout << "像素的最小值" << arrMaxMin[1] << endl;

	// 图像金字塔信息
	if(pRasterBand->GetOverviewCount()!=0) // 获取当前波段的金字塔层数，如果没有金字塔返回0
	{ 
		cout << "金字塔数目：" << pRasterBand->GetOverviewCount() << endl;
		cout << "-------------------------------------------------------" << endl;

	}
	return;
}

int main() {

	// 读取数据
	loadRasterData();
	system("pause");
	return 0;
}

参考文章01
参考文章02