一种用于黑土有机质信息提取的自适应高斯低通滤波算法制造技术

本发明专利技术属于信息提取技术领域，具体涉及一种用于黑土有机质信息提取的自适应高斯低通滤波算法；本发明专利技术的目的是针对现有技术的缺陷，提供一种用于黑土有机质信息提取的自适应高斯低通滤波算法；包括以下步骤：步骤一、黑土地高光谱预处理；步骤二、频率域转化；步骤三、自适应高斯滤波；步骤四、分级提取；步骤五、有机质提取结果；将各个级别含量的有机质进行综合制图，形成有机质含量的航空高光谱提取结果；将总体精度和Kappa系数最高的参数设置σ1＝0.5，ε1＝1.1应用到黑土地有机质的运算中，分级之后，根据地面分析数据，给出每一级别的有机质含量，得出黑土有机质含量的空间分布情况。

An Adaptive Gauss Low Pass Filter for Extracting Organic Matter Information from Black Soil

【技术实现步骤摘要】
一种用于黑土有机质信息提取的自适应高斯低通滤波算法
本专利技术属于信息提取
，具体涉及一种用于黑土有机质信息提取的自适应高斯低通滤波算法。
技术介绍
我国东北黑土富含有机质，引入航空高光谱数据能够为科学评估黑土地质量提供新的技术手段。在基于高光谱数据的黑土有机质含量识别中，将光谱数据理解为反映黑土成分含量目标电磁辐射特性的能量分布图，这一处理仅仅描述了地物能量分布的光谱特性，是一种黑土单像元分析的方法，能够提取的黑土信息有限。为掌握遥感光学系统信息传递规律，从遥感光学系统的频域信道矩阵推导出发，从频域的角度分析对成像系统频域信息传递参数进行数值计算，研究成果对遥感光学系统的设计起到了参考作用。常规方法是，将黑土地像元光谱的能量表达与频域频谱能量的表征进行分开处理，这种处理不能将原始光谱数据重新分布到另一个空间，无法反映空域图像中能量的变化规律，使得黑土的色调、纹理、走向和边界等空间特征未能参与到土质评估中，其评估结果常常不是一个综合的黑土特征集合，限制了黑土有机质信息提取的精度。因此，迫切需要在分析黑土频谱特征的基础上，设计一种自适应高斯低通滤波算法，通过逐步排除不同能量等级的数据，达到黑土地分类分级的目的。联合黑土地高光谱像元的空间上下文信息和光谱信息对高光谱像元进行分类，突破影像的空间信息突破传统逐像素光谱分类的性能瓶颈，为黑土有机质大范围精确制图提供支撑，来解决现有的方法的弊端。
技术实现思路
本专利技术的目的是针对现有技术的缺陷，提供一种用于黑土有机质信息提取的自适应高斯低通滤波算法。本专利技术的技术方案是：一种用于黑土有机质信息提取的自适应高斯低通滤波算法，包括以下步骤：步骤一、黑土地高光谱预处理；辐射校正使用基于MODTRAN4+辐射传输模型的FLAASH算法对高光谱数据进行预处理。计算时需要元数据包括观测视场角、太阳角度、平均海拔高度、大气模型、气溶胶类型和能见度范围；几何校正利用机载POS510系统实现，机上惯导系统和定位系统记录了每一个像元的位置参数X，Y，Z和姿态参数Roll，Pitch，Heading，通过每一帧图像的GPS时间，将坐标赋值给该象元；步骤二、频率域转化；将原始黑土地高光谱数据进行频域转换，生成频率域数据；步骤三、自适应高斯滤波；结合黑土地理空间中已知有机质数据的插值数据，设计能够自适应的高斯滤波器，使得在保留黑土图像局部特征的基础上，选择不同的σ，实现原始图像的分类；这一过程中，不断与地面化验数据进行比对，迭代运算，直至满足分级要求；步骤四、分级提取；设计一种自适应高斯低通滤波算法，通过逐步排除不同能量等级的数据，实现黑土地有机质含量分类分级的目的；步骤五、有机质提取结果；将各个级别含量的有机质进行综合制图，形成有机质含量的航空高光谱提取结果；将总体精度和Kappa系数最高的参数设置σ1＝0.5，ε1＝1.1应用到黑土地有机质的运算中，分级之后，根据地面分析数据，给出每一级别的有机质含量，得出黑土有机质含量的空间分布情况。所述步骤一包括以下步骤：步骤1.1辐射校正计算在高光谱传感器处接收的像元光谱辐射亮度：式中，L为传感器接收的总辐射亮度；ρ为像元表面反射率；ρe为像元周围平均表面反射率；S为大气球面反照率；Lα为大气后向散射辐射率；A、B为大气和几何条件的系数；方法的机理在于使用基于MODTRAN4+的大气辐射传输模型；步骤1.2几何校正原始图像中记录了每一帧图像的GPS时间，经过辐射校正后生成了一个时间标志记录文件*.att文件，该文件中记录了每一帧图像的GPS时间，该GPS时间都包含在POS系统结果文件的时间段内，通过时间比对和坐标投影变换就可以获得每一帧图像的姿态数据和位置数据，对图像进行几何校正；所述步骤二包括以下步骤：步骤2.1建立基本的滤波模型在频率域中，基本的滤波模型为：G(u,v)＝H(u,v)·F(u,v)式中，F(u,v)为被滤波的傅里叶变换图像；H(u,v)为滤波器变换函数；通过衰减后高频信息后生成平滑后的图像G(u,v)。步骤2.2高斯低通滤波器的建立能够将空间域转化为频率域的高斯低通滤波器为：式中，D(u,v)为距变换原点的距离，d为高斯曲线扩展的程度，即截止频率。设计一种自适应高斯滤波器，根据被平滑图像的黑土含量特征，自动选择不同的d，使得在处理后的结果图像中得到相应含量的有机质数据。所述步骤三包括以下步骤：步骤3.1高斯平滑的表达假设黑土影像的高斯平滑用下式表示：I0(x,y)＝Id(x,y)+ed(x,y)式中，对于(x,y)处的像素点，I0(x,y)为原始高光谱数据的灰度值，Id(x,y)为截止频率d下的低通灰度值，εd(x,y)为截止频率d下的残差值。步骤3.2自适宜高斯滤波的算法实现过程输入：黑土高光谱影像x＝(x1,x2,…,xi)∈Rd×i，滤波有机质含量级别C∈(1,2,3,4,5)，黑土影像分类集合ξ＝{1,2,…,n}，训练集Tm＝{(x1,c1),…,(xm,cm)}∈(Rd×ξ)m，其中，c∈ξ且m代表训练样本的数量；a)：输入训练集，构建高斯低通滤波器，对高光谱黑土影像进行逐像素分类获取初始分类结果c；b)：将地面化验数据插值为面域信息，结果表示为类概率图ξn，n＝1,2,…,N；c)：对高光谱影像x进行初始截止频率d1计算，获得初始方差σ1和残差值ε1；d)：ifd1＝dbestthene)：所提取信息符合某一级别有机质含量区间，则记为c1f)：endg)：elseh)：d1’＝d1++；i)：对高光谱影像x进行初始截止频率d1’计算，获得新方差σ1’和残差值ε1’；j)：直至，所提取信息符合某一级别有机质含量区间，则记为c1；k)：对高光谱数据x进行逆变换，去除空域中相应像素，生成x’；l)：endm)：对高光谱影像x’进行初始截止频率d2计算，获得初始方差σ2和残差值ε2；n)：重复d～m步骤，直至生成c1、c2、c3、c4和c5输出：有机质含量分类结果c＝(c1,c2,…,ci)，分类结果中每一个像素i被赋予一个标记ci∈ξ所述步骤四包括以下步骤：步骤4.1计算截止频率d设计一种计算截止频率d的方法，达到既能平滑黑土高光谱数据实现分类分级，又能保持与化验值最近似的目的。设计一种基于能量函数的计算方法：式中，c为常数项，根据采样点的有机质化验数据确定；σ为方差；ε为残差值。步骤4.2平滑分类与细节保持原则得出，在c为已知常数项下，最优截止频率d需要方差σ尽可能大，将高光谱数据变换为更加平滑的数据。而残差ε必须尽可能小，即经过高斯滤波后，原始像素(x,y)处的反射率值变化幅度越小越好。这样，建立了一种在因素有机质含量的综合高斯滤波方法，达到了平滑分类与保持细节特征的平衡。本专利技术的有益效果是：一是自适应高斯低通滤波算法不仅使得平滑区域像素分类精度提升，同时使得分级结果中的田块有机质含量边缘检测更加精确。地块平滑滤波能够有效的联合影像的空间与光谱信息，在有效去除类似传统分类方法造成的椒盐噪声基础上，使得分类结果中黑土地的轮廓与真实的黑土地地块轮廓基本一致；二是相比与传统高光谱分类算法，本章提出的频域识别与提取算法计算效率较高，在明确有机质含量等级后，能够快速计算出每一等级的结果数据；频域方法从新的角度对高本文档来自技高网...

【技术保护点】
1.一种用于黑土有机质信息提取的自适应高斯低通滤波算法，其特征在于：包括以下步骤：步骤一、黑土地高光谱预处理；辐射校正使用基于MODTRAN4+辐射传输模型的FLAASH算法对高光谱数据进行预处理。计算时需要元数据包括观测视场角、太阳角度、平均海拔高度、大气模型、气溶胶类型和能见度范围；几何校正利用机载POS 510系统实现，机上惯导系统和定位系统记录了每一个像元的位置参数X，Y，Z和姿态参数Roll，Pitch，Heading，通过每一帧图像的GPS时间，将坐标赋值给该象元；步骤二、频率域转化；将原始黑土地高光谱数据进行频域转换，生成频率域数据；步骤三、自适应高斯滤波；结合黑土地理空间中已知有机质数据的插值数据，设计能够自适应的高斯滤波器，使得在保留黑土图像局部特征的基础上，选择不同的σ，实现原始图像的分类；这一过程中，不断与地面化验数据进行比对，迭代运算，直至满足分级要求；步骤四、分级提取；设计一种自适应高斯低通滤波算法，通过逐步排除不同能量等级的数据，实现黑土地有机质含量分类分级的目的；步骤五、有机质提取结果；将各个级别含量的有机质进行综合制图，形成有机质含量的航空高光谱提取结果；将总体精度和Kappa系数最高的参数设置σ1＝0.5，ε1＝1.1应用到黑土地有机质的运算中，分级之后，根据地面分析数据，给出每一级别的有机质含量，得出黑土有机质含量的空间分布情况。...

【技术特征摘要】
1.一种用于黑土有机质信息提取的自适应高斯低通滤波算法，其特征在于：包括以下步骤：步骤一、黑土地高光谱预处理；辐射校正使用基于MODTRAN4+辐射传输模型的FLAASH算法对高光谱数据进行预处理。计算时需要元数据包括观测视场角、太阳角度、平均海拔高度、大气模型、气溶胶类型和能见度范围；几何校正利用机载POS510系统实现，机上惯导系统和定位系统记录了每一个像元的位置参数X，Y，Z和姿态参数Roll，Pitch，Heading，通过每一帧图像的GPS时间，将坐标赋值给该象元；步骤二、频率域转化；将原始黑土地高光谱数据进行频域转换，生成频率域数据；步骤三、自适应高斯滤波；结合黑土地理空间中已知有机质数据的插值数据，设计能够自适应的高斯滤波器，使得在保留黑土图像局部特征的基础上，选择不同的σ，实现原始图像的分类；这一过程中，不断与地面化验数据进行比对，迭代运算，直至满足分级要求；步骤四、分级提取；设计一种自适应高斯低通滤波算法，通过逐步排除不同能量等级的数据，实现黑土地有机质含量分类分级的目的；步骤五、有机质提取结果；将各个级别含量的有机质进行综合制图，形成有机质含量的航空高光谱提取结果；将总体精度和Kappa系数最高的参数设置σ1＝0.5，ε1＝1.1应用到黑土地有机质的运算中，分级之后，根据地面分析数据，给出每一级别的有机质含量，得出黑土有机质含量的空间分布情况。2.根据权利要求1所述的一种用于黑土有机质信息提取的自适应高斯低通滤波算法，其特征在于：所述步骤一包括以下步骤：步骤1.1辐射校正计算在高光谱传感器处接收的像元光谱辐射亮度：式中，L为传感器接收的总辐射亮度；ρ为像元表面反射率；ρe为像元周围平均表面反射率；S为大气球面反照率；Lα为大气后向散射辐射率；A、B为大气和几何条件的系数；方法的机理在于使用基于MODTRAN4+的大气辐射传输模型；步骤1.2几何校正原始图像中记录了每一帧图像的GPS时间，经过辐射校正后生成了一个时间标志记录文件*.att文件，该文件中记录了每一帧图像的GPS时间，该GPS时间都包含在POS系统结果文件的时间段内，通过时间比对和坐标投影变换就可以获得每一帧图像的姿态数据和位置数据，对图像进行几何校正。3.根据权利要求1所述的一种用于黑土有机质信息提取的自适应高斯低通滤波算法，其特征在于：所述步骤二包括以下步骤：步骤2.1建立基本的滤波模型在频率域中，基本的滤波模型为：G(u,v)＝H(u,v)·F(u,v)式中，F(u,v)为被滤波的傅里叶变换图像；H(u,v)为滤波器变换函数；通过衰减后高频信息后生成平滑后的图像G(u,v)。步骤2.2高斯低通滤波器的建立能够将空间域转化为频率域的高斯低通滤波器为：式中，D(u,v)为距变换原...

【专利技术属性】
技术研发人员：张东辉，赵英俊，李璇琼，秦凯，杨越超，张玉燕，喻翔，
申请(专利权)人：核工业北京地质研究院，
类型：发明
国别省市：北京,11