航空γ能谱数据NASVD降噪有效奇异值数目确定方法技术

本发明专利技术涉及一种航空γ能谱数据NASVD降噪有效奇异值数目确定方法，构建输入谱线数据矩阵后，利用NASVD算法得到主成分矩阵和幅度矩阵，接着利用“滑窗法”对获取的谱线主成分进行分类，对分类后的谱线主成分进行谱线重构，得到系列重构谱线。然后对系列重构谱线中相邻谱线之间的相关性进行分析计算，得到重构谱线的系列相关系数。以相关系数的绝对值作为判断依据，确定有效奇异值的数目。本发明专利技术通过一种简单、高效的方法能够快速确定航空γ能谱数据NASVD降噪中有效奇异值的数目，提高了降噪质量、缩短了处理时间周期，避免人工靠经验或推荐值确定有效奇异值数目产生过降噪或产生结果系统性误差，避免返工和重复计算，降低操作者的工作强度。者的工作强度。者的工作强度。

【技术实现步骤摘要】
航空
γ
能谱数据NASVD降噪有效奇异值数目确定方法


[0001]本专利技术涉及一种航空γ能谱降噪方法，具体地说是一种航空γ能谱数据NASVD降噪有效奇异值数目确定方法。

技术介绍

[0002]航空γ能谱测量是根据天然或人工放射性核素的γ射线能量差异，利用航空γ能谱测量系统确定地表岩石、土壤和大气中的钾、铀、钍和其它放射性核素含量及其分布的一种方法，是开展战略性矿产勘查、环境辐射监测与评价、核设施监测、核应急航空监测的一种重要技术手段，在众多应用领域中发挥着重要作用。近年来随着资源与环境对航空测量的急速增长，尤其对航空γ能谱测量技术提出了新的更高的要求。
[0003]γ能谱降噪是航空γ能谱测量中的关键技术之一，是提高放射性核素含量估算精度和减小谱定性定量分析误差的关键步骤，多年里一直是国内外研究者的研究热点。目前，航空γ能谱采用的最主要、最有效的谱线降噪方法是国际原子能机构IAEA
‑
1363报告中推荐的NASVD（Noise Adjusted Singular Value Decomposition，噪声调整奇异值分解）。NASVD是一种基于多元统计分析的多道能谱降噪方法，即通过主成分分析技术从原始γ能谱数据集中提取出主成分谱，通过选择低阶主成分谱重构谱线来实现降噪。自1997年NASVD方法提出以来，国内外学者对其降噪效果以及与传统方法比对进行了深入研究，大量研究表明该方法不仅能够显著地减少γ能谱中的噪声，而且在克服传统降噪方法不足方面具有显著的优势，因而在国内外大量的航测项目中得到广泛应用。
[0004]NASVD方法降噪的关键是如何选取最佳的低阶主成分个数（即有效奇异值数目）来重构谱线。有效奇异值如果选择个数过少，则有可能将能谱信号分量当作噪声去掉，造成重要能谱异常的丢失；选择个数过多，可能导致能谱信号中噪声信息成分过多，造成重要微弱异常淹没在噪声中而不能被有效识别。因此，有效奇异值数目选取至关重要，选取不当将可能造成整个测区或区域航测结果分析和解释的偏差，由此可能造成不可挽回的损失，尤其是拥有海量γ能谱数据的重特大航测项目。
[0005]目前，有效奇异值数目确定方法主要是通过观察主成分谱线形状进行人工判断。在“国际原子能机构IAEA
‑
1363报告”中，关于有效奇异值个数选取问题仅提出不能仅依靠特征值的大小进行判断，同时还要考虑参考空间信息，且也仅给出了推荐值。在“EJ/T 1032
‑
2005航空伽马能谱测量规范”中，关于自行设计的滤波器，要求滤波类型和参数的选择要经过反复试验进行确定。查阅的国内外相关文献中，涉及内容大多数是降噪效果比对及改进方法，而关于如何选取有效奇异值个数内容，只简单描述为观察确定，而未给出确定的理论依据。在国际专业处理软件Praga4中，开发了NASVD方法处理模块，但有效奇异值数目还是需要依靠观察主成分谱线形状由人工确定。
[0006]而仅靠人工判断方法在项目实践中却存在很大不足，主要表现在：1、易发生过降噪现象。降噪效果取决于操作者的工作经验，人工判断缺乏理论依据，难以避免过降噪现象的发生，往往造成重要的航测成果可能丢失或难以识别而被忽略。
[0007]2、易发生结果系统性偏差，造成项目降噪处理工作全部返工。项目实践中，NASVD方法通常采用一条/多条测线、单/多架次或单/多/全部测区数据集进行降噪处理，涉及成百万甚至几千万条谱线数据。由于测区环境的复杂性，导致噪声信号大量存在，原始信号的主成分谱线形状往往难以区分，致使有效奇异值数目难以确定，仅依靠人工判断达到最佳降噪效果，对操作人员专业素质要求高且工作效率十分低下，易产生结果系统性偏差且难以溯源，甚至会造成项目降噪处理工作全部返工。
[0008]3、重复计算工作量大，操作人员工作强度高，降噪质量难以保证。
[0009]4、自动化程度低，处理时间周期长。难以满足堪称海量航空γ能谱数据快速处理和分析的要求。目前，运用新技术对历史航测资料的二次开发涉及近500万平方公里的海量能谱数据，工作量巨大且不能满足数据快速处理和分析的要求。
[0010]所以，如果能够通过一种方法准确确定航空γ能谱数据NASVD降噪有效奇异值的数目，将有效避免由于人为因素造成的过降噪现象的发生，提高NASVD方法降噪的自动化水平和工作效率，进一步保障放射性核素含量估计精度和减小航空γ能谱定性定量分析的误差。

技术实现思路

[0011]本专利技术的目的就是提供一种航空γ能谱数据NASVD降噪有效奇异值数目确定方法，以解决现在航空γ能谱数据NASVD降噪中人工判断有效奇异值数目容易发生过降噪现象、产生结果系统性偏差、重复计算工作量大以及降噪质量难以保证的问题。
[0012]本专利技术是这样实现的：一种航空γ能谱数据NASVD降噪有效奇异值数目确定方法，包括以下步骤。
[0013]a.读入原始航空γ能谱数据，构建输入谱线数据矩阵A。
[0014]b.利用NASVD算法对输入谱线数据矩阵A进行计算，得到主成分矩阵S和幅度矩阵B，主成分矩阵S的各个列向量依次为各主成分，且各主成分按照对应特征值的大小降序排列。
[0015]c.对主成分按照前1个、前2个、前3个
…
前n个进行分类。
[0016]d.分别对前1个、前2个、前3个
…
前n个主成分进行重构，分别得到重构谱矩阵A1，A
2 ，A3，
…
A
n
。
[0017]e.依次计算相邻两个重构谱矩阵的相关系数|r
i
|。
[0018]f.判断|r
i
|的大小，当|r
i
|第一次等于1时，i值即被确定为有效奇异值的数目。
[0019]进一步地，在步骤b中，主成分矩阵S和幅度矩阵B的求解步骤如下。
[0020]1）对输入谱线数据矩阵A各行元素求和，公式如下：R=(r
11
，r
21
，
…ꢀ
r
m1
)
T
式中，R表示输入谱线数据矩阵A各行元素之和组成的列向量，r
i1
表示输入谱线数据矩阵A中第i行元素之和，其中i=1，2，
…
，m。
[0021]2）对输入谱线数据矩阵A各列元素求和，公式如下：C=(c
11
，c
12
，
…
，c
1n
)式中，C表示输入谱线数据矩阵A中各列元素之和组成的行向量，c
1j
表示输入谱线数据矩阵A中第j列元素之和，其中j=1，2，
…
，n。
[0022]3）归一化单位谱，公式如下：式中，N表示对C进行归一化处理后的行向量。
[0023]4）求取单位方差矩阵，公式如下：式中，A
new
表示m
×
n单位方差矩阵。本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种航空γ能谱数据NASVD降噪有效奇异值数目确定方法，其特征在于，包括以下步骤：a.读入原始航空γ能谱数据，构建输入谱线数据矩阵A；b.利用NASVD算法对输入谱线数据矩阵A进行计算，得到主成分矩阵S和幅度矩阵B，主成分矩阵S的各个列向量依次为各主成分，且各主成分按照对应特征值的大小降序排列；c.对主成分按照前1个、前2个、前3个
…
前n个进行分类；d.分别对前1个、前2个、前3个
…
前n个主成分进行重构，分别得到重构谱矩阵A1，A
2 ，A3，
…
A
n
；e.依次计算相邻两个重构谱矩阵的相关系数|r
i
|；f.判断|r
i
|的大小，当|r
i
|第一次等于1时，i值即被确定为有效奇异值的数目。2.根据权利要求1所述的航空γ能谱数据NASVD降噪有效奇异值数目确定方法，其特征在于，在步骤b中，主成分矩阵S和幅度矩阵B的求解步骤如下：1）对输入谱线数据矩阵A各行元素求和，公式如下：R=(r
11
，r
21
，
…ꢀ
r
m1
)
T
式中，R表示输入谱线数据矩阵A各行元素之和组成的列向量，r
i1
表示输入谱线数据矩阵A中第i行元素之和，其中i=1，2，
…
，m；2）对输入谱线数据矩阵A各列元素求和，公式如下：C=(c
11
，c
12
，
…
，c
1n
)式中，C表示输入谱线数据矩阵A中各列元素之和组成的行向量，c
1j
表示输入谱线数据矩阵A中第j列元素之和，其中j=1，2，
…
，n；3）归一化单位谱，公式如下：式中，N表示对C进行归一化处理后的行向量；4）求取单位方差矩阵，公式如下：式中，A
new
表示m
×
n单位方差矩阵；5）对A
new
进行奇异值分解，公式如下：A
new
=UΣV
T
式中，...

【专利技术属性】
技术研发人员：张胜，王强，王永军，王瑞军，唐首普，王虎，王会波，张春雷，孟凡兴，
申请(专利权)人：核工业航测遥感中心，
类型：发明
国别省市：