[发明专利]一种基于能量对应的嫦娥二号伽玛谱非线性谱漂校正方法

说明书

技术领域

本发明涉及一种基于能量对应的嫦娥二号伽玛谱非线性谱漂校正方法。

背景技术

嫦娥二号伽玛谱仪为嫦娥二号卫星的主要载荷之一，其科学目标为探测月表的主要物质成分。嫦娥二号伽玛谱仪获取的谱数据在处理过程中发现存在明显的非线性谱漂现象。由于谱线在解析之前需要进行全谱数据的累加平均，谱线漂移引起的峰位和峰区偏移造成了累加后元素特征峰信息被平均湮没。

传统的伽玛谱谱漂校正方法分为有参考源和无参考源两类，有参考源的谱漂校正根据参考源放出的γ射线在仪器谱中形成一个相对稳定的参考峰，当谱线发生漂移后以此参考峰位基准完成校正；无参考源的谱漂校正主要是利用计算机软件实现漂移校正，其中校正方法包含有线性校正和非线性校正两类，线性校正算法主要是针对谱线的线性漂移，而目前的非线性漂移的算法对于大数据量的嫦娥二号伽玛谱适用性不强。

发明内容

为解决上述技术问题，本发明的目的是提供一种基于能量对应的嫦娥二号伽玛谱非线性谱漂校正方法。

本发明的目的通过以下的技术方案来实现：

一种基于能量对应的嫦娥二号伽玛谱非线性谱漂校正方法，该方法包括：

(1)读入待校正谱线y(n)，n为道址数；新建数组p(n)，用以存入校正后谱线；

(2)根据谱线y(n)的特征峰峰位和能量信息，对y(n)的能量-道址对应关系进行线性拟合，获得关系式En_y＝a₂*CH_y+b₂，其中CH_y为待校正谱的道址；En_y为道址对应的能量值；a₂为直线的斜率，表示每道谱线对应的能量间隔；b₂为直线的截距，即第零道谱线对应的能量值；

(3)计算谱线y(n)的每道对应标准谱的新的道址，第i道对应的标准谱的道址为：

chi=[(a2*i+b2)-b1a1];]]>

其中a₁为标准谱能量-道址线性关系的斜率；b₁标准谱能量-道址线性关系的截距；

(4)判断p(ch_i)是否为零，若为零，则p(ch_i)＝y(i)；否则

(5)循环执行上述步骤(3)和步骤(4)，直至全谱数据处理完毕。

与现有技术相比，本发明的一个或多个实施例可以具有如下优点：

该方法能够有效地校正谱线的非线性漂移，使得获取谱线的元素特征峰峰位和峰区处于一致的范围内，为后续嫦娥二号伽玛谱的解析以及元素月球分布有用信息的获取提供了有力的帮助。

附图说明

图1是标准伽玛谱和待校正谱对比图；

图2是标准谱和已校正后伽玛谱对比图；

图3是Fe峰校正前后峰位对比变化图；

图4是Mg/K峰校正前后峰位对比变化图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚，下面将结合实施例及附图对本发明作进一步详细的描述。

本实施例提供了一种基于能量对应的嫦娥二号伽玛谱非线性谱漂校正方法，首先设x(n)为伽玛谱标准谱线，谱线的能量-道址对应关系为：

En_x＝a₁*CH_x+b₁(1)

其中CH_x为标准谱的道址，En_y为该道址对应能量值。

设y(n)为待校正的谱漂谱线，根据y(n)谱线的特征峰信息，对y(n)进行能量-道址校正，可得y(n)谱线的能量-道址关系为：