基于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注册