申请/专利权人：浙江迪谱诊断技术有限公司

申请日：2021-12-27

公开（公告）日：2024-04-12

公开（公告）号：CN114487072B

主分类号：G01N27/62

分类号：G01N27/62;G06F17/15

优先权：

专利状态码：有效-授权

法律状态：2024.04.12#授权;2024.03.26#著录事项变更;2022.05.31#实质审查的生效;2022.05.13#公开

摘要：本发明公开了一种飞行时间质谱峰拟合方法，包括：将整条飞行时间核酸质谱图按照一定间隔长度进行分段，并计算出各分段强度阈值和二阶导数阈值：在每个分段内找到离子强度过零点的一阶导数，确定可能的高斯峰位置；从可能的高斯峰位置中，选取符合条件的离子强度，得到确定的高斯峰位置，确定出所有的高斯峰个数l，及高斯峰的质荷比mzi和峰高Hi；根据半经验公式，计算得到高斯峰的峰宽Wi；将高斯峰的质荷比mzi和峰宽Wi，取4Wi作为高斯拟合区域，作为高斯峰最小二乘法迭代初始值；将高斯拟合区域4Wi内的所有采样点的核酸质谱数据带入高斯拟合函数fxi中进行迭代求解，得到精确的拟合高斯峰的质荷比mzfit、峰宽Wfit、峰高Hfit。

主权项：1.一种飞行时间质谱峰拟合方法，其特征在于：包括以下步骤：步骤1：获取整条飞行时间核酸质谱图，该飞行时间核酸质谱图，以采样点的质荷比为横坐标，以采样点的离子强度为纵坐标；对于第i个采样点，i∈{1,2,3,…，n}，其核酸质谱数据表示为Ii,mzi，其中，Ii表示离子强度，mzi表示质荷比；步骤2：将整条飞行时间核酸质谱图按照一定间隔长度进行分段，并计算出各分段强度阈值和二阶导数阈值，对于第j分段，其强度阈值表示为Ithresholdj和二阶导数阈值表示为I″thresholdj：对于第j分段，按照下式计算强度阈值Ithresholdj和二阶导数阈值I″thresholdj： 式中，为第j分段内离子强度的平均值，σj为第j分段内离子强度的标准差；a为信号强度阈值系数，a∈{1,2,3}；d为特定系数，取值为2，Hj和Wj为第j分段内峰高最小的特征峰的峰高和半高半宽；步骤3：在每个分段内找到离子强度过零点的一阶导数，该一阶导数对应的原始的质荷比，为可能的高斯峰位置；步骤4：从可能的高斯峰位置中，选取符合以下条件的离子强度，该离子强度对应的原始的质荷比，为确定的高斯峰位置： 式中，I′i+1j为第j分段内第i+1采样点的离子强度Ii+1j的一阶导数I′i+1j；确定出所有的高斯峰个数l，及高斯峰的质荷比mzi和峰高Hi，峰高Hi由以下公式计算得到：对于第i个采样点，当质谱基线ci为零，峰高Hi＝Ii；当质谱基线ci不为零，峰高Hi＝Ii-ci；步骤5：根据半经验公式，计算得到高斯峰的峰宽Wi＝2×0.001×mzi；步骤6：将高斯峰的质荷比mzi和峰宽Wi，取4Wi作为高斯拟合区域，作为高斯峰最小二乘法迭代初始值；步骤7：将高斯拟合区域4Wi内的所有采样点的核酸质谱数据带入高斯拟合函数fxi中进行迭代求解： 其中，Hfit、mzfit、Wfit分别为拟合高斯峰的峰高、质荷比和峰宽，α为特征峰宽系数，取值为0.60056120439323；质谱基线cfit根据下式进行求解：cfit＝A\I9其中，I为由高斯拟合区域内所有采样点的离子强度Ii组成的向量，表示为I＝I1,I2,I3,…,IN，N为高斯拟合区域内所有采样点的总数，A为向量内部所有数值均为1且长度与向量I相同，表示为A＝11,12,13,…,1N；将求解得到的拟合高斯峰的峰高Hfit和质荷比mzfit与作为迭代初始值的质荷比mzi和峰高Hi进行匹配，判断以下条件是否成立： 式中，为第k次和第k+1次迭代拟合求解得到的质荷比，k为大于1的正整数，Δ1、Δ2和Δ3为微小量常数；若不成立，则将拟合高斯峰的峰高Hfit和质荷比mzfit作为下一次迭代初始值，进行新一轮的求解直到满足式10的条件；若成立，则拟合求解结束，得到精确的拟合高斯峰的质荷比mzfit、峰宽Wfit、峰高Hfit。

全文数据：

权利要求：

百度查询： 浙江迪谱诊断技术有限公司 一种飞行时间质谱峰拟合方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD