CN106908655B - 光电信号峰值检测方法与系统 - Google Patents

光电信号峰值检测方法与系统 Download PDF

Info

Publication number
CN106908655B
CN106908655B CN201710128877.2A CN201710128877A CN106908655B CN 106908655 B CN106908655 B CN 106908655B CN 201710128877 A CN201710128877 A CN 201710128877A CN 106908655 B CN106908655 B CN 106908655B
Authority
CN
China
Prior art keywords
segment data
photosignal
peak
maximum point
crest
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201710128877.2A
Other languages
English (en)
Other versions
CN106908655A (zh
Inventor
郭成祥
吴本涛
陈绪兵
师伟展
陈俊璘
易雁飞
李鑫鑫
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangdong Shunde Industrial Design Institute
Original Assignee
Guangdong Shunde Industrial Design Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Guangdong Shunde Industrial Design Institute filed Critical Guangdong Shunde Industrial Design Institute
Priority to CN201710128877.2A priority Critical patent/CN106908655B/zh
Publication of CN106908655A publication Critical patent/CN106908655A/zh
Application granted granted Critical
Publication of CN106908655B publication Critical patent/CN106908655B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
    • G01R19/25—Arrangements for measuring currents or voltages or for indicating presence or sign thereof using digital measurement techniques
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64—Fluorescence; Phosphorescence
    • G01N21/6428—Measuring fluorescence of fluorescent products of reactions or of fluorochrome labelled reactive substances, e.g. measuring quenching effects, using measuring "optrodes"

Landscapes

  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Immunology (AREA)
  • General Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Optics & Photonics (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Biological Materials (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明提供一种光电信号峰值检测方法与系统，对原始光电信号分段读取，并对分段读取的每段数据前端和后端插入预设量平均值参数，选取母小波函数对处理后的分段数据进行一阶求导和二阶求导，提取一阶导数极大值点与第二分段数据的二阶导数极大值点，确定各第二分段数据中波峰位置，根据各第二分段数据中波峰位置，计算各第二分段数据中波峰宽度以及波峰高度，获取原始光电信号的波峰峰值。整个过程中，对原始光电信号的分段前端和后端分别补入平均值参数，提高对原始光电信号的处理精度、采用母小波函数对补入之后的分段数据进行一阶导数和二阶导数处理，准确确定波峰位置，另外还无需对原始光电数据进行复杂的滤波处理，显著提高原始光电信号峰值检测的效率。

Description

光电信号峰值检测方法与系统
技术领域
本发明涉及信号检测技术领域，特别是涉及光电信号峰值检测方法与系统。
背景技术
目前在信号检测技术领域主要分为模拟信号检测和数字信号检测两类，某些时候针对相同设备(仪器)输出的信号可以采用模拟信号检测，也可以采用数字信号检测，例如针对数字PCR(Polymerase Chain Reaction，聚合酶链式反应)仪检测的信号检测中，模拟信号检测和数字信号检测均可采用。
采用不同的信号检测方式其可以带来不同的效果，一般来说，针对数字类型的设备(例如PCR仪)采用数字信号检测方式会更好，但是在使用数字信号检测方法时，就会产生大量采集的光电信号数据，因此在处理原始信号数据时，需要一种高精确度并且能够迅速提取出光电信号峰值的方法。
发明内容
基于此，有必要针对目前尚无一种高精确度且高效的光电信号峰值检测方法的问题，提供一种高精确度且高效的光电信号峰值检测方法与系统。
一种光电信号峰值检测方法，包括步骤：
分段读取原始光电信号，获取原始光电信号的第一分段数据集合；
获取各第一分段数据的平均值参数，并对各第一分段数据的最前端和最后端分别补入预设量对应的第一分段数据的平均值参数，获得第二分段数据集合；
选取母小波函数对各第二分段数据进行一阶求导和二阶求导；
提取各第二分段数据的一阶导数极大值点和第二分段数据的二阶导数极大值点；
根据各第二分段数据的一阶导数极大值点和第二分段数据的二阶导数极大值点，确定各第二分段数据中波峰位置；
根据各第二分段数据中波峰位置，计算各第二分段数据中波峰宽度以及波峰高度，获取原始光电信号的波峰峰值。
一种光电信号峰值检测系统，包括：
分段读取模块，用于分段读取原始光电信号，获取原始光电信号的第一分段数据集合；
补入模块，用于获取各第一分段数据的平均值参数，并对各第一分段数据的最前端和最后端分别补入预设量对应的第一分段数据的平均值参数，获得第二分段数据集合；
求导模块，用于选取母小波函数对各第二分段数据进行一阶求导和二阶求导；
极大值提取模块，用于提取各第二分段数据的一阶导数极大值点和第二分段数据的二阶导数极大值点；
波峰确定模块，用于根据各第二分段数据的一阶导数极大值点和第二分段数据的二阶导数极大值点，确定各第二分段数据中波峰位置；
波峰值确定模块，用于根据各第二分段数据中波峰位置，计算各第二分段数据中波峰宽度以及波峰高度，获取原始光电信号的波峰峰值。
本发明光电信号峰值检测方法与系统，对原始光电信号分段读取，并对分段读取的每段数据前端和后端插入预设量平均值参数，选取母小波函数对处理后的分段数据进行一阶求导和二阶求导，提取一阶导数极大值点与第二分段数据的二阶导数极大值点，确定各第二分段数据中波峰位置，根据各第二分段数据中波峰位置，计算各第二分段数据中波峰宽度以及波峰高度，获取原始光电信号的波峰峰值。整个过程中，对原始光电信号的分段前端和后端分别补入平均值参数，提高对原始光电信号的处理精度、采用母小波函数对补入之后的分段数据进行一阶导数和二阶导数处理，准确确定波峰位置，另外还无需对原始光电数据进行复杂的滤波处理，显著提高原始光电信号峰值检测的效率。
附图说明
图1为本发明光电信号峰值检测方法第一个实施例的流程示意图；
图2为微滴荧光采集系统示意图；
图3为原始光电信号数据波形示意图；
图4为原始光电信号数据局部(第二分段数据)放大波形示意图；
图5为第二分段数据及其一阶导数波形示意图；
图6为第二分段数据及其一阶导数与二阶导数波形示意图；
图7为选取的一组相邻极大值点示意图；
图8为由一组相邻的一阶导数极大值和二阶导数极大值求解出来的峰值点(局部放大)示意图；
图9为本发明光电信号峰值检测系统第一个实施例的结构示意图。
具体实施方式
如图1所示，一种光电信号峰值检测方法，包括步骤：
S100：分段读取原始光电信号，获取原始光电信号的第一分段数据集合。
原始光电信号是从仪器设备上直接采集到的光电信号。在实际应用中，仪器设备中的光电信号系统将光信号转换为电信号，此时电信号属于模拟信号，再经过模数转换变为数字信号，以便进行下一步的处理。为便于解释说明原始光电信号的来源，下面将以微滴式数字PCR仪输出的信号为例详细解释。微滴式数字PCR系统是在一般的PCR扩增前对样品进行微滴化处理，即将含有核酸分子的反应体系分成上万个纳升级的微滴，其中每个微滴不含待检测核酸靶分子，或者含有一个至数个待检测核酸靶分子。经过PCR扩增后，逐个对每个微滴进行检测，有荧光信号的微滴判断为1，没有荧光信号的微滴判断为0，然后根据泊松分布原理及阳性微滴的个数与比例即可得出靶分子的起始拷贝数或者浓度。具体如图2所示，微滴荧光采集系统中滴液生成器生成直径几乎相同的N个微滴，发光系统激发微滴荧光，1号光信号采集系统和2号光信号采集系统采集微滴荧光信号→通过光电转换系统，将光信号转换为电信号→通过模数转换信号，将模拟信号转换为数字信号输出，即得到原始光电信号。
一般来说原始光电信号是一个数据量庞大的数据，在对其进行处理过程中需要对进行分段读取，获得一段一段的原始光电信号，再针对每一段原始光电信号单独进行后续相同的处理，最终将每一段原始光电信号处理结果汇聚，得到原始光电信号最终的处理结果。非必要的，针对每一段原始光电信号可以采用并行处理的方式进行处理，以便有效提高数据处理效率。另外，在第一分段数据集合中第一分段数据的数量可以为单个或为多个，所有的第一分段数据汇聚即为原始光电信号，具体如图3与图4所示，图4所示为图3中原始光电信号的一段(第一分段数据)。
S200：获取各第一分段数据的平均值参数，并对各第一分段数据的最前端和最后端分别补入预设量对应的第一分段数据的平均值参数，获得第二分段数据集合。
针对每一个第一分段数据，分别获取其对应的平均值参数，在每个第一分段数据的最前端和最后端分别补入一定量对应的第一分段数据的平均值参数，获得第二分段数据集合。在这里，补入量的大小可以基于历史经验数据或者试验数据进行设定，最终得到的第二分段数据的数量与第一分段数据的数量相等。例如步骤S100对原始光电数据分段处理之后得到1号第一分段数据、2号第一分段数据和3号第一分段数据，步骤S200即获取1号第一分段数据的平均值参数O、2号第一分段数据的平均值参数P以及3号第一分段数据的平均值参数Q，并在1号第一分段数据的前端和后端补入一定量平均值参数O、在2号第一分段数据的前端和后端补入一定量平均值参数P、在3号第一分段数据的前端和后端补入一定量平均值参数Q，得到1号第二分段数据、2号第二分段数据以及3号第二分段数据。在每段信号都加入一定数量的原始光电信号的平均值，在后续使用小波变换求导的时候并不会出现旁瓣效应，因此不会丢失原始光电信号的最前端和最后端部分数据的处理，从而提高了对原始光电信号的处理精度。第二分段数据具体如图4所示。
S300：选取母小波函数对各第二分段数据进行一阶求导和二阶求导。
选用母小波函数对每个第二分段数据进行一阶求导和二阶求导，具体来说，母小波函数可以选择多贝西(Daubeichies)小波，高斯(Gaussian)小波及其导数等满足母小波条件的函数。具体来说，可以采用以下方式选取最优母小波函数。
由小波变换的性质可得公式(1)、公式(2)，公式(1)中：Wf(a,b)为小波系数、f(x)为原始信号，a为缩放尺度，b为平移尺度，公式(2)中：为所选择的小波类型，其中*表示该函数的共轭函数，为卷积运算符。
由公式(3)可得，当选择合适的函数便可对函数f(x)求其n阶导数，因此当选择的g(x)函数的n阶导数满足构成母小波函数的条件时，就可以求出原始信号的n阶导数。具体来说，在实际应用中例如微滴式PCR仪器检测系统的原始光电信号数据的波形类似为高斯波形，并且高斯函数的导数都可以满足构成母小波的条件，所以在对微滴式PCR仪光信号进行峰值检测时可以选择高斯函数的一阶导数和二阶导数作为求解原始数据一阶导数和二阶导数的母小波函数。
在公式(1)中缩放尺度a与平移尺度b的参数值可以基于历史经验数据获得。更一步的，在一些情况下还可以获得缩放尺度a与平移尺度b的最优解。例如针对本发明光电信号峰值检测方法应用于微滴式PCR仪光信号检测时，由于微滴式PCR的微滴直径都是固定的参数，因此每个正常波峰的宽度是大致相等的，所以在选择小波的尺度参数时，该值是能够通过实验结果来确定最优的尺度参数a。当确定尺度参数之后，就可以直接对原始信号求取一阶导数和二阶导数。针对原始光电信号进行一阶求导和二阶求导的结果具体如图5以及图6所示。
S400：提取各第二分段数据的一阶导数极大值点和第二分段数据的二阶导数极大值点。
如图5以及图6所示针对步骤S300一阶求导和二阶求导的结果，提取各第二分段数据的一阶导数极大值点和第二分段数据的二阶导数极大值点。非必要的，还可以基于历史经验数据对一阶导数极大值设置阈值，并针对二阶导数极大值设置阈值，剔除不在阈值范围内一阶导数极大值点和二阶导数极大值点。
S500：根据各第二分段数据的一阶导数极大值点和第二分段数据的二阶导数极大值点，确定各第二分段数据中波峰位置。
基于步骤S400得到原始一阶导数和二阶导数的极大值点，即可确定第二分段数据中波峰位置。具体来说，确定过程主要包括：
步骤一：针对单个第二分段数据，选取一组分别与一阶导数极大值点A相邻的二阶导数极大值点B与二阶导数极大值点C；
步骤二：当一阶导数极大值点A对应的横坐标参数大于二阶导数极大值点B对应的横坐标参数、且小于二阶导数极大值点C对应的横坐标参数时，确定第二分段数据中波峰的横坐标参数为二阶导数极大值点B的横坐标参数与二阶导数极大值点C的横坐标参数之和的二分之一
在具体实施例中，具体如图7所示，提取一组相邻的一阶导数极大值点为(x1，y1)，和二阶导数极大值点分别为(x2，y2)、(x2’，y2’)，当这三个点满足条件：x2<x1<x2’,那么就可以判断在(x2，x2’)中间存才一个原始信号的波峰，即x2、x2’的平均值(x2+x2’)/2就是该波峰的位置Xp，若不满足直接剔除三个极值点。
S600：根据各第二分段数据中波峰位置，计算各第二分段数据中波峰宽度以及波峰高度，获取原始光电信号的波峰峰值。
原始光电信号的波峰峰值还包括波峰宽度和波峰高度，在步骤S500中已经求得波峰位置，基于波峰位置以及第二分段数据中原始光电信号的整个波形数据，即可计算出波峰宽度以及波峰高度，具体如图8所示。获得当前第二分段数据的波峰峰值，将各个第二分段数据的波峰峰值汇聚，最终获得原始光电信号的波峰峰值。非必要的，步骤300至步骤S600中的计算过程可以通过FPGA(Field-Programmable Gate Array，现场可编程门阵列)实现，另外还可以通过其他软件编程语言(譬如C语言，C#，C++，Java，Matlab，R语言等)编写的程序来实现。
本发明光电信号峰值检测方法，对原始光电信号分段读取，并对分段读取的每段数据前端和后端插入预设量平均值参数，选取母小波函数对处理后的分段数据进行一阶求导和二阶求导，提取一阶导数极大值点与第二分段数据的二阶导数极大值点，确定各第二分段数据中波峰位置，根据各第二分段数据中波峰位置，计算各第二分段数据中波峰宽度以及波峰高度，获取原始光电信号的波峰峰值。整个过程中，对原始光电信号的分段前端和后端分别补入平均值参数，提高对原始光电信号的处理精度、采用母小波函数对补入之后的分段数据进行一阶导数和二阶导数处理，准确确定波峰位置，另外还无需对原始光电数据进行复杂的滤波处理，显著提高原始光电信号峰值检测的效率。
在其中一个实施例中，计算各第二分段数据中波峰宽度包括步骤：
步骤一：针对单个第二分段数据，选取一组对应的波峰位置与一阶导数极大值点。
步骤二：根据选取的对应的波峰位置以及对应波峰的一阶导数极大值点，计算单个第二分段数据中波峰宽度。
当求出各个波峰位置之后，可以根据一组对应的波峰位置Xp和一阶导数极大值x1，求出该波峰的宽度Xd＝2(Xp-x1)。
在其中一个实施例中，计算各第二分段数据中波峰高度包括步骤：
步骤一：针对单个第二分段数据，在第二分段数据中波峰位置两侧各提取预设数量的原始数据点，获得第一原始数据点集合。
步骤二：剔除第一原始数据点集合中最大值和最小值，获得第二原始数据点集合。
步骤三：计算第二原始数据点中各原始数据点纵坐标的平均值参数，获得单个第二分段数据中波峰高度。
在波峰位置Xp处左右各取m(m可以根据实验情况来设置合适的值)个原始数据点，然后剔除这2m+1个点中的最大值和最小值，最后求出的2m-1(包含Xp)个点纵坐标的平均值作为该波峰高度，设波峰高度为Xh。
在其中一个实施例中，根据各第二分段数据中波峰位置，计算各第二分段数据中波峰宽度以及波峰高度，获取原始光电信号的波峰峰值的步骤包括：
步骤一：计算各第二分段数据中波峰宽度以及波峰高度。
步骤二：根据切比雪夫不等式设定波峰宽度阈值与波峰高度阈值。
步骤三：分别剔除各第二分段数据中波峰宽度不在波峰宽度阈值范围内的波峰以及波峰高度不在波峰高度阈值范围内的波峰峰值，获得原始光电信号的波峰峰值。
峰值评估。经过以上处理可以得到一系列的波峰峰值(Xp，Xd，Xh)，Xp为波峰位置，Xd为波峰宽度，Xh为波峰高度，然后可以根据需要设定波峰宽度和波峰高度阈值(根据切比雪夫不等式设定阈值)，并剔除波峰宽度Xd和波峰高度Xh为异常值(不在阈值范围之内)的波峰，从而保留符合要求的波峰峰值。在具体应用于针对微滴式PCR仪光电信号波峰检测时，此处的峰值评估是为了剔除不合格微滴产生的峰值，从而提高了微滴式PCR仪检测系统中微滴的假阳性和假阴性的识别率。
如图9所示，一种光电信号峰值检测系统，包括：
分段读取模块100，用于分段读取原始光电信号，获取原始光电信号的第一分段数据集合。
补入模块200，用于获取各第一分段数据的平均值参数，并对各第一分段数据的最前端和最后端分别补入预设量对应的第一分段数据的平均值参数，获得第二分段数据集合。
求导模块300，用于选取母小波函数对各第二分段数据进行一阶求导和二阶求导。
极大值提取模块400，用于提取各第二分段数据的一阶导数极大值点和第二分段数据的二阶导数极大值点。
波峰确定模块500，用于根据各第二分段数据的一阶导数极大值点和第二分段数据的二阶导数极大值点，确定各第二分段数据中波峰位置。
波峰值确定模块600，用于根据各第二分段数据中波峰位置，计算各第二分段数据中波峰宽度以及波峰高度，获取原始光电信号的波峰峰值。
本发明光电信号峰值检测系统，分段读取模块100对原始光电信号分段读取，补入模块200对分段读取的每段数据前端和后端插入预设量平均值参数，求导模块300选取母小波函数对处理后的分段数据进行一阶求导和二阶求导，极大值提取模块400提取一阶导数极大值点与第二分段数据的二阶导数极大值点，波峰确定模块500确定各第二分段数据中波峰位置，波峰值确定模块600根据各第二分段数据中波峰位置，计算各第二分段数据中波峰宽度以及波峰高度，获取原始光电信号的波峰峰值。整个过程中，对原始光电信号的分段前端和后端分别补入平均值参数，提高对原始光电信号的处理精度、采用母小波函数对补入之后的分段数据进行一阶导数和二阶导数处理，准确确定波峰位置，另外还无需对原始光电数据进行复杂的滤波处理，显著提高原始光电信号峰值检测的效率
在其中一个实施例中，波峰确定模块500包括：
选取单元，用于针对单个第二分段数据，选取一组分别与一阶导数极大值点A相邻的二阶导数极大值点B与二阶导数极大值点C。
确定单元，用于当一阶导数极大值点A对应的横坐标参数大于二阶导数极大值点B对应的横坐标参数、且小于二阶导数极大值点C对应的横坐标参数时，确定第二分段数据中波峰的横坐标参数为二阶导数极大值点B的横坐标参数与二阶导数极大值点C的横坐标参数之和的二分之一。
在其中一个实施例中，波峰值确定模块600还用于：
针对单个第二分段数据，选取一组对应的波峰位置与一阶导数极大值点；
根据选取的对应的波峰位置以及对应波峰的一阶导数极大值点，计算单个第二分段数据中波峰宽度。
在其中一个实施例中，波峰值确定模块600还用于：
针对单个第二分段数据，在第二分段数据中波峰位置两侧各提取预设数量的原始数据点，获得第一原始数据点集合；
剔除第一原始数据点集合中最大值和最小值，获得第二原始数据点集合；
计算第二原始数据点中各原始数据点纵坐标的平均值参数，获得单个第二分段数据中波峰高度。
在其中一个实施例中，波峰值确定模块600包括：
参数计算单元，用于计算各第二分段数据中波峰宽度以及波峰高度；
阈值计算单元，用于根据切比雪夫不等式设定波峰宽度阈值与波峰高度阈值；
筛选单元，用于分别剔除各第二分段数据中波峰宽度不在波峰宽度阈值范围内的波峰以及波峰高度不在波峰高度阈值范围内的波峰峰值，获得原始光电信号的波峰峰值。
以上实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。

Claims (10)

1.一种光电信号峰值检测方法，其特征在于，包括步骤：
分段读取原始光电信号，获取所述原始光电信号的第一分段数据集合；
获取各第一分段数据的平均值参数，并对各所述第一分段数据的最前端和最后端分别补入预设个数对应的第一分段数据的平均值参数，获得第二分段数据集合；
选取母小波函数对各第二分段数据进行一阶求导和二阶求导；
提取各所述第二分段数据的一阶导数极大值点和所述第二分段数据的二阶导数极大值点；
根据各所述第二分段数据的一阶导数极大值点和所述第二分段数据的二阶导数极大值点，确定各所述第二分段数据中波峰位置；
根据各所述第二分段数据中波峰位置，计算各所述第二分段数据中波峰宽度以及波峰高度，获取所述原始光电信号的波峰峰值。
2.根据权利要求1所述的光电信号峰值检测方法，其特征在于，所述根据各所述第二分段数据的一阶导数极大值点和所述第二分段数据的二阶导数极大值点，确定各所述第二分段数据中波峰位置的步骤包括：
针对单个所述第二分段数据，选取一组分别与一阶导数极大值点A相邻的二阶导数极大值点B与二阶导数极大值点C；
当所述一阶导数极大值点A对应的横坐标参数大于所述二阶导数极大值点B对应的横坐标参数、且小于所述二阶导数极大值点C对应的横坐标参数时，确定所述第二分段数据中波峰的横坐标参数为所述二阶导数极大值点B的横坐标参数与所述二阶导数极大值点C的横坐标参数之和的二分之一；
根据所述第二分段数据中波峰的横坐标参数，确定所述第二分段数据中波峰位置。
3.根据权利要求1所述的光电信号峰值检测方法，其特征在于，计算各所述第二分段数据中波峰宽度包括步骤：
针对单个所述第二分段数据，选取一组对应的波峰位置与一阶导数极大值点；
根据选取的对应的波峰位置以及对应波峰的一阶导数极大值点，计算单个所述第二分段数据中波峰宽度。
4.根据权利要求1所述的光电信号峰值检测方法，其特征在于，计算各所述第二分段数据中波峰高度包括步骤：
针对单个所述第二分段数据，在所述第二分段数据中波峰位置两侧各提取预设数量的原始数据点，获得第一原始数据点集合；
剔除所述第一原始数据点集合中最大值和最小值，获得第二原始数据点集合；
计算所述第二原始数据点集合中原始数据点纵坐标的平均值参数，获得单个所述第二分段数据中波峰高度。
5.根据权利要求1所述的光电信号峰值检测方法，其特征在于，所述根据各所述第二分段数据中波峰位置，计算各所述第二分段数据中波峰宽度以及波峰高度，获取所述原始光电信号的波峰峰值的步骤包括：
计算各所述第二分段数据中波峰宽度以及波峰高度；
根据切比雪夫不等式设定波峰宽度阈值与波峰高度阈值；
分别剔除各所述第二分段数据中波峰宽度为异常值即不在波峰宽度阈值范围内的波峰以及波峰高度为异常值即不在波峰高度阈值范围内的波峰峰值，获得所述原始光电信号的波峰峰值。
6.一种光电信号峰值检测系统，其特征在于，包括：
分段读取模块，用于分段读取原始光电信号，获取所述原始光电信号的第一分段数据集合；
补入模块，用于获取各第一分段数据的平均值参数，并对各所述第一分段数据的最前端和最后端分别补入预设个数对应的第一分段数据的平均值参数，获得第二分段数据集合；
求导模块，用于选取母小波函数对各第二分段数据进行一阶求导和二阶求导；
极大值提取模块，用于提取各所述第二分段数据的一阶导数极大值点和所述第二分段数据的二阶导数极大值点；
波峰确定模块，用于根据各所述第二分段数据的一阶导数极大值点和所述第二分段数据的二阶导数极大值点，确定各所述第二分段数据中波峰位置；
波峰值确定模块，用于根据各所述第二分段数据中波峰位置，计算各所述第二分段数据中波峰宽度以及波峰高度，获取所述原始光电信号的波峰峰值。
7.根据权利要求6所述的光电信号峰值检测系统，其特征在于，所述波峰确定模块包括：
选取单元，用于针对单个所述第二分段数据，选取一组分别与一阶导数极大值点A相邻的二阶导数极大值点B与二阶导数极大值点C；
确定单元，用于当所述一阶导数极大值点A对应的横坐标参数大于所述二阶导数极大值点B对应的横坐标参数、且小于所述二阶导数极大值点C对应的横坐标参数时，确定所述第二分段数据中波峰的横坐标参数为所述二阶导数极大值点B的横坐标参数与所述二阶导数极大值点C的横坐标参数之和的二分之一。
8.根据权利要求6所述的光电信号峰值检测系统，其特征在于，所述波峰值确定模块还用于：
针对单个所述第二分段数据，选取一组对应的波峰位置与一阶导数极大值点；
根据选取的对应的波峰位置以及对应波峰的一阶导数极大值点，计算单个所述第二分段数据中波峰宽度。
9.根据权利要求6所述的光电信号峰值检测系统，其特征在于，所述波峰值确定模块还用于：
针对单个所述第二分段数据，在所述第二分段数据中波峰位置两侧各提取预设数量的原始数据点，获得第一原始数据点集合；
剔除所述第一原始数据点集合中最大值和最小值，获得第二原始数据点集合；
计算所述第二原始数据点集合中原始数据点纵坐标的平均值参数，获得单个所述第二分段数据中波峰高度。
10.根据权利要求6所述的光电信号峰值检测系统，其特征在于，所述波峰值确定模块包括：
参数计算单元，用于计算各所述第二分段数据中波峰宽度以及波峰高度；
阈值计算单元，用于根据切比雪夫不等式设定波峰宽度阈值与波峰高度阈值；
筛选单元，用于分别剔除各所述第二分段数据中波峰宽度为异常值即不在波峰宽度阈值范围内的波峰以及波峰高度为异常值即不在波峰高度阈值范围内的波峰峰值，获得所述原始光电信号的波峰峰值。
CN201710128877.2A 2017-03-06 2017-03-06 光电信号峰值检测方法与系统 Active CN106908655B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710128877.2A CN106908655B (zh) 2017-03-06 2017-03-06 光电信号峰值检测方法与系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710128877.2A CN106908655B (zh) 2017-03-06 2017-03-06 光电信号峰值检测方法与系统

Publications (2)

Publication Number Publication Date
CN106908655A CN106908655A (zh) 2017-06-30
CN106908655B true CN106908655B (zh) 2019-06-14

Family

ID=59186792

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710128877.2A Active CN106908655B (zh) 2017-03-06 2017-03-06 光电信号峰值检测方法与系统

Country Status (1)

Country Link
CN (1) CN106908655B (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107312850A (zh) * 2017-07-19 2017-11-03 华东医药（杭州）基因科技有限公司 一种pcr无效扩增的检测方法
CN108444954B (zh) * 2017-12-08 2020-09-11 广东顺德工业设计研究院(广东顺德创新设计研究院) 光谱信号峰值检测方法、装置以及系统
CN108318735B (zh) * 2018-01-19 2020-04-07 河海大学常州校区 高压交流或脉冲电源单次供电能量分离方法
CN109583029B (zh) * 2018-11-01 2022-02-18 郑州云海信息技术有限公司 一种消除受Ledge影响信号边沿斜率值的方法及装置
CN111272060B (zh) * 2018-12-05 2021-08-31 中国石油大学(华东) 基于单对电极电容成像检测技术的非导电材料开口缺陷宽度方向尺寸量化方法
CN111089856B (zh) * 2019-12-26 2021-05-14 厦门大学 一种拉曼光谱弱信号提取的后处理方法
CN111358454B (zh) * 2020-03-17 2022-07-29 乐普（北京）医疗器械股份有限公司 一种基于差异波占比检测信号质量的方法和装置
CN111623958B (zh) * 2020-05-18 2021-11-12 长春欧意光电技术有限公司 干涉信号中的小波峰峰值提取方法
CN112461805A (zh) * 2020-11-16 2021-03-09 三诺生物传感股份有限公司 一种用于荧光强度基底计算的方法
CN113640445B (zh) * 2021-08-11 2024-06-11 贵州中烟工业有限责任公司 基于图像处理的特征峰识别方法及计算设备、存储介质
CN115684737B (zh) * 2022-10-26 2023-11-17 常州同惠电子股份有限公司 一种用于计算波形毛刺的算法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1117832B1 (en) * 1998-10-08 2009-04-08 Hybaid Limited Detection of nucleic acid polymorphism
CN102187215A (zh) * 2008-09-09 2011-09-14 生物辐射实验室股份有限公司 基于回归的多级pcr分析系统
CN104950215A (zh) * 2015-05-08 2015-09-30 广西大学 一种微机保护方法
CN105067650A (zh) * 2015-08-03 2015-11-18 东南大学 一种利用小波计算导数检测谱特征峰的方法
CN105286815A (zh) * 2015-11-02 2016-02-03 重庆大学 一种基于波形时域特征的脉搏波信号特征点检测方法
CN105518455A (zh) * 2013-09-09 2016-04-20 株式会社岛津制作所 波峰检测方法
CN105637360A (zh) * 2013-10-17 2016-06-01 株式会社岛津制作所 波形中的峰值端点检测方法及检测装置
CN105954565A (zh) * 2016-05-26 2016-09-21 中国电力科学研究院 一种高压直流线路电晕电流信号的提取方法
CN106137184A (zh) * 2015-04-08 2016-11-23 四川锦江电子科技有限公司 基于小波变换的心电信号qrs波检测方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7228237B2 (en) * 2002-02-07 2007-06-05 Applera Corporation Automatic threshold setting and baseline determination for real-time PCR
US7124949B2 (en) * 2003-06-17 2006-10-24 Symbol Technologies, Inc. Reducing speckle noise in electro-optical readers
CA2606604A1 (en) * 2005-05-13 2006-11-23 Bio-Rad Laboratories, Inc. Baselining amplification data
WO2014094039A1 (en) * 2012-12-19 2014-06-26 Rmit University A background correction method for a spectrum of a target sample

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1117832B1 (en) * 1998-10-08 2009-04-08 Hybaid Limited Detection of nucleic acid polymorphism
CN102187215A (zh) * 2008-09-09 2011-09-14 生物辐射实验室股份有限公司 基于回归的多级pcr分析系统
CN105518455A (zh) * 2013-09-09 2016-04-20 株式会社岛津制作所 波峰检测方法
CN105637360A (zh) * 2013-10-17 2016-06-01 株式会社岛津制作所 波形中的峰值端点检测方法及检测装置
CN106137184A (zh) * 2015-04-08 2016-11-23 四川锦江电子科技有限公司 基于小波变换的心电信号qrs波检测方法
CN104950215A (zh) * 2015-05-08 2015-09-30 广西大学 一种微机保护方法
CN105067650A (zh) * 2015-08-03 2015-11-18 东南大学 一种利用小波计算导数检测谱特征峰的方法
CN105286815A (zh) * 2015-11-02 2016-02-03 重庆大学 一种基于波形时域特征的脉搏波信号特征点检测方法
CN105954565A (zh) * 2016-05-26 2016-09-21 中国电力科学研究院 一种高压直流线路电晕电流信号的提取方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Feature extraction using first and second derivative extrema (FSDE) for real-time and hardware-efficient spike sorting;Sivylla E. Paraskevopoulou et al.;《Journal of Neuroscience Methods》;20130131;第215卷(第1期);第29-37页
基于小波变换的心电图QRS波群峰值检测;张文琼等;《计算机仿真》;20051231;第22卷(第12期);第147-150页
基于改进峰值搜索法的旋转机械瞬时频率估计;胡爱军等;《振动与冲击》;20131231;第32卷(第7期);第113-117页

Also Published As

Publication number Publication date
CN106908655A (zh) 2017-06-30

Similar Documents

Publication Publication Date Title
CN106908655B (zh) 光电信号峰值检测方法与系统
Xie et al. A deep-learning-based real-time detector for grape leaf diseases using improved convolutional neural networks
CN112308158B (zh) 一种基于部分特征对齐的多源领域自适应模型及方法
CN110533084A (zh) 一种基于自注意力机制的多尺度目标检测方法
CN105488770B (zh) 一种面向对象的机载激光雷达点云滤波方法
CN107076712B (zh) 色谱数据处理方法以及装置
CN107220277A (zh) 基于手绘草图的图像检索算法
CN110414577A (zh) 一种基于深度学习的激光雷达点云多目标地物识别方法
CN110032954A (zh) 一种钢筋智能识别与计数方法及系统
CN105223110B (zh) 一种显微镜定位跟踪成像方法、装置及尿沉渣分析系统
CN109544554B (zh) 一种植物图像分割及叶片骨架提取方法及系统
CN104504389A (zh) 一种基于卷积神经网络的卫星云量计算方法
JP2008523354A (ja) フロー・アナライザからデータを管理するための装置と方法
CN105190553B (zh) 用于检测重合样本事件的方法
CN104502329B (zh) 一种libs成分分析中元素最优工作曲线的建立方法
CN116522153B (zh) 锂电池容量预测方法、装置、计算机设备和存储介质
CN109448795A (zh) 一种circRNA的识别方法及装置
CN101226133A (zh) 一种血细胞脉冲信号的分类识别方法
CN110490836A (zh) dPCR微阵列图像信息处理方法
CN107480466A (zh) 基因组数据存储方法及电子设备
CN105096293A (zh) 用于处理尿液沉渣图像的待处理区块的方法和装置
CN107862326A (zh) 一种基于全卷积神经网络的成熟苹果识别方法
US9773090B2 (en) Method computer program and system to analyze mass spectra
CN108734221A (zh) 基于子类数目自动确定的改进密度峰值聚类方法及系统
CN109902800A (zh) 基于拟反馈神经网络的多级骨干网络检测通用物体的方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant