基于GDAL的栅格图像空间插值预处理——C语言版

基于GDAL的栅格图像预处理

前言

        栅格数据和矢量数据构成空间数据的主要来源，怎样以开源方式读取并处理这些空间数据？目前有多种开源支持包，这里只介绍GDAL包。GDAL包的优点是支持库简洁、支持栅格和矢量、与多种开发平台结合。OpenGis方式读取空间数据，有利于自己编写程序进行图像预处理和智能识别等等，比如：遥感影像的降噪、锐化；红外图像的林火识别；工厂监控视频识别等等。本文中利用GDAL包读取高程栅格DEM，并添加气象自动站点的数据，进行空间插值研究。

 

一、程序主要程序功能实现过程

第一步：读取栅格数据，包含坡向和DEM

第二步：读入站点信息数据

第三步：按照行列号读取栅格单元到内存，不考虑高程为0的单元

第三步：每一个坡向情况中，参考固定数目的气象站点。首先确定搜索范围，获取定量数目的监测站点。

第四步:在同一个计算窗口内，通过给各个因子赋权重，依据海拔高程回归关系与加权回归分析得到

温度递减率。

第五步：计算单个站点的温度预测值，然后计算所有站点的距离权重因子，根据因子大小确定综合

影响后的温度值

第四步：开辟新内存存储处理后的栅格数据，然后新建一个tiff格式的文件，把内存

数据导出到该文件中

 

二、代码示例

int _tmain(int argc, _TCHAR* argv[])  
{
	/*
	DEM、坡向栅格数据的数据框大小
	*/
	int XsizeDEM;
	int YsizeDEM;
	int XsizeAspect;
	int YsizeAspect;
	//double geoTransform[6];
	//double Xp,Yp;


	/*
	DEM、坡向栅格单元对象的VALUE值
	*/
	short int *pmemDEM;
	float *pmemAspect;
	float *pmemNew;

	//GDAL注册