[发明专利]液质色谱峰噪声估计方法、装置、存储介质、系统

说明书

本发明涉及种液质色谱峰噪声估计方法、装置、存储介质、系统。该方法包括；获取待处理的色谱峰谱图；利用分割算法对谱图进行分割，得到若干分割图像区域块；根据若干分割图像区域块中的图像数据获取噪声的期望估计值；根据期望估计值计算噪声的估计值。通过对色谱图进行处理，以对输入色谱峰的噪声进行估计，无需输入任何形如期望、方差等的先验参数，解决应用人员使用门槛高、应用场景窄的问题。

技术领域

本发明涉及检测技术领域，特别涉及一种液质色谱峰噪声估计方法、装置、存储介质、系统。

背景技术

目前在检验领域中液质联用是一项重要的检验手段，而液质色谱图中往往噪声信号较多，对于噪声估计存在一定困难。如果对液质色谱图的噪声估计不恰当，会导致信噪比计算和寻峰计算出现偏差，往往会影响实验结果，甚至会影响后续峰面积的计算，干扰待测样本的浓度准确性。

本申请旨在建立一种系统性解决液质色谱峰噪声估计方法及实施系统。

发明内容

为了实现根据本发明的上述目的和其他优点，本发明的第一目的是提供一种液质色谱峰噪声估计方法，包括以下步骤：

获取待处理的色谱峰谱图；

利用分割算法对所述谱图进行分割，得到若干分割图像区域块；

根据所述若干分割图像区域块中的图像数据获取噪声的期望估计值；

根据所述期望估计值计算噪声的估计值。

优选地，所述分割算法还包括步骤：

对获取的所述谱图进行分割，若所述谱图的数据长度为N，将数据分为若干段，每一段的数据长度均为L，若最后一段的数据点不足L，则舍去最后一段，基于此原则可以得到谱图被分为K段，其中K表示为：

符号表示向下取整；

所述谱图被分段后，第一段数据记作S(1,t)，依此类推，第n段数据记作S(n,t)，最后一段数据记作S(K,t)，其中t表示所述段数据的离散时间顺序，取值为t＝1,…,L。

优选地，还包括步骤：

分别对K段数据做离散傅里叶变换，公式如下：

F(n,w)＝FFT[S(n,t)]for n＝1,…,K (2)

其中频率w的取值为

依次提取F(1,w),…,F(K,w)中的第一个点组成新的集合，该集合表示为：A(0)＝{F²(1,0),…,F²(K,0)}；F²(1,w),…,F²(K,w)的所有第i+1个点组成的集合为这样的集合有L个，分别为A(0),…,A(L-1)；

分别计算出每个集合中的最小值，即：

F²(w)＝min{F²(1,w),…,F²(K,w)} (3)

其中

此时序列F(0),…,即被认为是噪声的频域估计，而F(0)为噪声的期望估计值。

优选地，还包括步骤：

根据所述F(0)计算噪声的估计值公式如下：

优选地，还包括平衡参数η，以对所述噪声的估计值进行修正；公式如下所示：