买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

申请/专利权人：郑州轻工业大学

摘要：本发明公开了一种基于离子峰谱线的FAIMS图谱特征提取方法，依次包括以下步骤：A：对原始谱线进行噪声去除并得到去噪后谱线，规定原始谱线的左基线范围和右基线范围；B：计算左基线波动和右基线波动；C：确定有效离子峰；D：判断离子峰的类型；E：对不同类型的离子峰分别提取离子峰特征，离子峰特征包括离子峰值特征、离子峰位置特征、半峰宽特征、离子峰峰值与半峰宽的比值特征和峰面积特征。本发明能够准确找到每个有效离子峰的形状和区域，并提取关于离子峰的5种特征。

主权项：1.一种基于离子峰谱线的FAIMS图谱特征提取方法，其特征在于，依次包括以下步骤：A：对原始谱线进行噪声去除并得到去噪后谱线，规定原始谱线的左基线范围和右基线范围；B：计算左基线波动和右基线波动；C：确定有效离子峰；D：判断离子峰的类型；E：对不同类型的离子峰分别提取离子峰特征，离子峰特征包括离子峰值特征、离子峰位置特征、半峰宽特征、离子峰峰值与半峰宽的比值特征和峰面积特征；其中：所述的步骤A中，首先对谱线上的尖脉冲噪声和锯齿波噪声进行去除，得到去噪后谱线，然后对原始谱线的左基线范围和右基线范围进行规定，当补偿电压处于最低或最高时，没有检测离子响应的两部分谱线分别作为左基线和右基线；所述的步骤B中，将原始谱线的左基线BLo的值减去去噪后谱线的左基线BLd的值，得到左基线波动FL；将原始谱线的右基线BRo的值减去去噪后谱线的右基线BRd的值，得到右基线波动FR；所述的步骤C中，在去噪后谱线中找出所有的驻点，并在驻点中确定峰点和谷点，通过设定最小离子峰峰值TP判断待判定的峰点是否为有效离子峰；判定去噪后谱线中的某一个待定点是否为驻点时，将去噪后谱线中的待定点减去待定点左侧相邻点的差值记为d1，待定点右侧相邻点减去待定点的差值记为d2；满足d1×d2＜0的待定点为驻点；驻点中满足d1＞0且d2＜0的点为峰点；驻点中满足d1＜0且d2＞0的点为谷点；最小离子峰峰值TP为：Tp＝kp×maxmedianFL,medianFR式1；如果待判定的峰点的值小于Tp，则认为该峰点是噪声；如果待判定的峰点的值大于Tp，则进入下一步的判断；计算待判定的峰点左右两个谷点的最小跨度值是否大于最小跨度阈值S，若大于最小跨度阈值S，则认为待判定的峰点为有效离子峰；若小于等于最小跨度阈值S，则认为是尖脉冲噪声，待判定的峰点不是有效离子峰；其中，kp为乘性系数，在一般正态分布假设下kp大于等于3；maxg,g表示取两个输入的最大者；mediang表示计算输入序列的中位数，FL为左基线波动，FR为右基线波动；两个谷点的跨度是指两个谷点位置的差的绝对值；所述的步骤D中，通过判断待判定有效离子峰峰顶的两个最近的谷点，与离子峰类型阈值TC的大小，确定离子峰的类型；离子峰类型阈值TC为：Tc＝meanBLd+kc×stdBLd；其中，kc为乘性系数，meang为计算输入序列的均值，stdg为计算输入序列的标准差，BLd为去噪后的左基线；若某个有效离子峰峰顶的两个最近的谷点都小于TC，则这个峰顶对应的离子峰为单峰；如果其中一个谷点小于TC，而另一个谷点大于TC，则对应的离子峰为连峰；如果其中两个谷点都大于TC，则对应的离子峰为夹峰；所述的步骤E中，根据离子峰的类型，分别提取离子峰特征：E1：对单峰提取5种离子峰特征；E11：将单峰所对应的峰值提取为离子峰值特征PV，峰值对应的补偿电压的位置提取为离子峰位置特征PI；E12：将单峰所对应的峰值的一半作为半峰值；在峰顶左侧上升的响应序列和右侧下降序列中分别确定一个与半峰值最接近的响应点，这两个响应点位置的差值的绝对值作为半峰宽特征HPW，峰值除以半峰宽作为离子峰峰值与半峰宽的比值特征VWR；E13：分别从该单峰的峰点左右最近的两个谷点出发，在朝向峰点的方向上连续的4个点中，如果3个相邻点的离子响应差值的绝对值都大于落差阈值Tf，则这4个点中的第1个点被认为是峰脚，两个峰脚内所有响应点的绝对值相加作为峰面积特征PA；E2：对连峰提取5种离子峰特征；E21：连峰包括左连峰和右连峰2种类型；对任何一种连峰，直接将该连峰所对应的峰值提取为离子峰值特征PV，峰值对应的补偿电压的位置提取为离子峰位置特征PI；E22：对于左连峰，取该左连峰右侧下降序列中与半峰值最接近的响应点位置，与该左连峰的离子峰位置差值的绝对值的2倍作为半峰宽特征HPW，峰值除以半峰宽被作为离子峰峰值与半峰宽的比值特征VWR；对于右连峰，取该右连峰左侧上升序列中与半峰值最接近的响应点位置，与该右连峰的离子峰位置差值的绝对值的2倍作为半峰宽特征HPW，峰值除以半峰宽被作为离子峰峰值与半峰宽的比值特征VWR；E23：对于左连峰，从该左连峰的峰点右侧最近的一个谷点出发，在朝向峰点的方向上连续的4个点中，如果3个相邻点的离子响应差值的绝对值都大于阈值Tf，则这4个点中的第1个点被认为是峰脚；将这个峰脚到左连峰峰点内所有响应点的绝对值相加的和的2倍作为峰面积特征PA；对于右连峰，从该右连峰的峰点左侧最近的一个谷点出发，在朝向峰点的方向上连续的4个点中，如果3个相邻点的离子响应差值的绝对值都大于阈值Tf，则这4个点中的第1个点被认为是峰脚；将这个峰脚到左连峰峰点内所有响应点的绝对值相加的和的2倍作为峰面积特征PA；E3：对夹峰提取5种离子峰特征；E31：将夹峰所对应的峰值提取为离子峰值特征PV，峰值对应的补偿电压的位置提取为离子峰位置特征PI；E32：分别在夹峰左侧上升序列和右侧下降序列中找到最大上升差值点和最大下降差值点；其中最大上升差值点为最大斜率点，最大下降差值点为最小斜率点；设某个点的差值为该点的响应值减去该点左侧相邻点的响应值；基于最大斜率点和最小斜率点之间的响应点，采用基于非线性最小二乘的高斯拟合计算出一个拟合曲线；根据夹峰拟合曲线直接提取半峰宽特征HPW，以及离子峰峰值与半峰宽的比值特征VWR特征；该夹峰拟合曲线左侧上升序列中与去噪后谱线的左基线BLd的均值最接近的点被作为左侧虚峰脚，该夹峰拟合曲线右侧下降序列中与去噪后谱线的右基线BRd的均值最接近的点被作为右侧虚峰脚，左虚峰脚与右虚峰脚之间所有响应点的绝对值相加的和作为峰面积特征PA。

全文数据：

权利要求：

百度查询： 郑州轻工业大学 一种基于离子峰谱线的FAIMS图谱特征提取方法