CN107250788B - 噪声水平估计方法和测定数据处理装置 - Google Patents
噪声水平估计方法和测定数据处理装置 Download PDFInfo
- Publication number
- CN107250788B CN107250788B CN201580076316.1A CN201580076316A CN107250788B CN 107250788 B CN107250788 B CN 107250788B CN 201580076316 A CN201580076316 A CN 201580076316A CN 107250788 B CN107250788 B CN 107250788B
- Authority
- CN
- China
- Prior art keywords
- section
- frequency
- noise level
- wave data
- value
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/86—Signal analysis
- G01N30/8675—Evaluation, i.e. decoding of the signal into analytical information
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/86—Signal analysis
- G01N30/8624—Detection of slopes or peaks; baseline correction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/86—Signal analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/86—Signal analysis
- G01N30/8603—Signal analysis with integration or differentiation
- G01N30/8617—Filtering, e.g. Fourier filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/86—Signal analysis
- G01N30/8624—Detection of slopes or peaks; baseline correction
- G01N30/8631—Peaks
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/62—Detectors specially adapted therefor
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/86—Signal analysis
- G01N30/8603—Signal analysis with integration or differentiation
- G01N30/861—Differentiation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/86—Signal analysis
- G01N30/8624—Detection of slopes or peaks; baseline correction
- G01N30/8641—Baseline
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/88—Integrated analysis systems specially adapted therefor, not covered by a single one of the groups G01N30/04 - G01N30/86
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01J—ELECTRIC DISCHARGE TUBES OR DISCHARGE LAMPS
- H01J49/00—Particle spectrometers or separator tubes
- H01J49/0027—Methods for using particle spectrometers
- H01J49/0036—Step by step routines describing the handling of the data generated during a measurement
Abstract
在本发明中,对色谱等测定数据进行时间频率分析并获取表现多个频率成分的强度随时间的变化的波形数据(步骤S1),在多个规定频率各自的波形数据中,将在时间轴方向上正值连续的部分和负值连续的部分分别设为一个区段来分割为多个区段(步骤S2),针对多个规定频率成分各自的波形数据求出多个区段各自的面积并计算区段值(步骤S3),制作由从多个规定频率成分各自的波形数据的多个区段中排除具有超过规定基准值的区段值的区段后得到的区段组成的选择区段组(步骤S4),基于选择区段组中包含的区段的区段值的平均值来求出多个规定频率成分各自的噪声水平(步骤S11)。
Description
技术领域
本发明涉及对色谱或光谱等测定数据中包含的噪声成分的大小(噪声水平)进行估计的方法、这种测定数据的处理装置以及测定数据处理用的程序。
背景技术
作为对液体试样中含有的成分进行分析的装置之一,存在一种液相色谱仪。在液相色谱仪中,使液体试样随着流动相的流动导入到柱内,在柱内将试样中的各成分按时间顺序分离之后用吸光光度计等检测器进行检测并制作色谱。然后,根据色谱上的峰位置来鉴定各成分,根据峰高度或面积来决定该成分的浓度(例如专利文献1)。
在通过测定得到的色谱中,一般不仅包含峰成分还包含噪声成分。峰成分的大小随着液体试样中含有的各种成分的洗脱而发生变动。另一方面,噪声成分的大小由于各种因素而发生变动。确定所有的这种变动因素来计算噪声成分的大小是不可能的,因此以往用白噪声对噪声成分进行近似并拟合于测定数据,或者根据通过测定获取到的色谱中的被估计为不存在峰成分的部分来估计整个色谱的噪声成分,由此从色谱中去除噪声成分后求出峰成分。
专利文献1:日本特开平7-98270号公报
专利文献2:日本特开2006-163614号公报
发明内容
发明要解决的问题
如上所述,实际的色谱的噪声成分由于各种因素而发生变动,因此有时利用白噪声来进行近似是不够的,无法高精度地估计噪声成分的大小。另外,为了根据通过测定获取到的色谱的一部分来估计噪声成分,优选对在色谱中完全不包含峰成分的部分进行估计,但这是困难的,难以高精度地估计噪声成分的大小。
在此以色谱为例进行了说明,但除此以外,在如分光分析时的分光光谱或质谱分析时的质谱那样包含峰成分和噪声成分的各种测定数据中也存在同样的问题。
本发明要解决的问题在于,提供能够高精度地对色谱和光谱等测定数据中包含的噪声成分的大小(噪声水平)进行估计的噪声水平估计方法、测定数据处理装置以及测定数据处理用程序。
用于解决问题的方案
图1示出峰成分(a)的波形、噪声成分(b)的波形、这些成分之和(c)的波形(Profile)以及对它们分别进行傅立叶变换所得到的功率谱(Power Spectrum)的一例。当将图1的各功率谱进行比较时,可知当在色谱中存在峰时,在比某个频率低的低频侧谱强度变大。本发明人着眼于这一点想到了本发明。
为了解决上述问题而完成的本发明的第一方式是一种噪声水平估计方法,根据通过测定相对于规定物理量发生变化的信号的强度来获取到的包含峰成分和噪声成分的测定数据,来估计噪声成分的大小,该方法的特征在于,
a)对所述测定数据进行时间频率分析,针对多个规定频率获取表现构成所述信号的该频率成分的强度相对于所述规定物理量的变化的波形数据;
b)在所述多个规定频率各自的波形数据中,将在所述物理量的变化方向上正值连续的部分和负值连续的部分分别设为一个区段,或者将在所述物理量的变化方向上相邻的极大值与极小值之间的部分设为一个区段,来将该波形数据分割为多个区段;
c)针对所述多个规定频率各自的波形数据求出所述多个区段各自的大小;
d)针对所述多个规定频率各自的波形数据来制作选择区段组,该选择区段组由从该波形数据的多个区段中排除大小超过规定基准值的区段后得到的区段组成;以及
e)通过求出所述选择区段组中包含的区段的大小的统计值,来求出所述多个规定频率各自的噪声水平。
上述规定物理量例如是时间、波长或者质荷比,上述测定数据例如是色谱、分光光谱或者质谱。
时间频率分析是在图像处理等领域中使用的分析法,具体地说,已知连续小波变换、离散小波变换以及滤波器组等(例如专利文献2)。在使用连续小波变换的情况下,获取连续的多个频率下的波形数据并从中取出多个频率的波形数据来使用。另外,在使用滤波器组或离散小波变换的情况下,预先决定多个频率,并获取这些多个频率的波形数据。上述多个规定频率既可以由分析者每次设定,或者也可以预先设定多个标准频率。
上述的所谓“时间频率分析”这一表述并非将上述规定物理量限定为时间。即,在上述规定物理量为波长的分光光谱数据、或上述规定物理量为质荷比的质谱数据等的分析中也能够使用本发明的方法。
上述区段的大小例如能够根据区段的面积或高度来求出。另外,上述规定基准值例如能够设为各区段的大小的平均值+Nσ(N为正整数,σ为无偏标准偏差),或者设为各区段的大小的中央值+M×MAD(Median Absolute Deviation:绝对中位差)(M为正整数)。或者,还能够根据同一波形数据中包含的多个区段的大小的分布求出包含该多个区段中的规定比例(例如90%)的区段的范围,并将该范围的上限值设为基准值。并且,上述选择区段组中包含的区段的大小的统计值例如能够设为构成选择区段组的多个区段的大小的平均值或中央值。
在本发明所涉及的噪声水平估计方法中,将超过规定基准值的大小的区段作为源自峰成分的区段而从波形数据中包含的多个区段中排除。因而,能够高精度地估计测定数据的噪声水平。
色谱中的噪声水平并不限于始终固定。例如,在梯度分析中,使构成流动相的多个溶液的混合比随时间而变化,因此有时噪声水平随之发生变动。另外,当在色谱的测定中装置周边的温度发生变化时,有时噪声水平也随之变动。在这种测定数据中,如果噪声水平变大,则区段变大,有可能将源自噪声的区段错误地排除掉。
因此,在本发明所涉及的噪声水平估计方法中,优选的是,
使用所述测定数据中的与所述物理量的变化方向上的噪声成分的大小的变动有关的指标数据,在制作所述选择区段组之前将所述多个区段各自的大小进行标准化。
作为上述指标数据,例如能够使用梯度分析中的多个溶液的混合比的变化的数据、或获取测定数据时的温度的变化的数据。通过将假定对噪声水平的增减造成影响的这些数据用作指标数据,能够制作更加准确的选择区段组来高精度地估计噪声水平。
另外,在本发明所涉及的噪声水平估计方法中,能够构成为:
排除如下的区段来制作选择区段组:该区段在所述规定物理量的变化方向上处于与已排除的所述区段对应的位置,并且属于频率比已排除的所述区段所属的波形数据的频率低的波形数据。
根据图1的(a)的峰成分的功率谱可知,当着眼于存在峰成分的某个频率时,在比该频率低的频带中必然存在峰成分。因此,在规定物理量的变化方向(例如时间轴方向)上,能够将处于与作为源自峰的区段而被排除的区段对应的位置的、属于频率比已排除的区段所属的波形数据的频率低的波形数据的区段也判断为源自峰的区段。因而,通过排除这些区段,能够无遗漏地排除源自峰的区段,来以更高的精度估计噪声水平。
并且,在本发明所涉及的噪声水平估计方法中,能够构成为:
将所述多个规定频率的噪声水平相互比较,以使低频侧的噪声水平成为高频侧的噪声水平以上的值的方式校正所述噪声水平。
在测定相对于规定物理量发生变化的信号强度时,一般使用具有电容器的电气电路或具有以规定的时间常数响应的A/D变换器的检测器。已知用于蓄积规定时间电荷的电容器或以规定的时间常数响应的A/D变换器与低通滤波器同样地减少高频带的信号。因此,当经由这种检测器获取信号时,存在噪声水平随着从高频侧趋向低频侧而变大的倾向。因而,以反映该倾向的方式校正噪声水平,由此能够更高精度地求出噪声水平。
为了解决上述问题而完成的本发明的第二方式是一种测定数据处理装置,用于根据通过测定相对于规定物理量发生变化的信号的强度来获取到的包含峰成分和噪声成分的测定数据,来估计噪声成分的大小,该测定数据处理装置的特征在于,具备:
a)时间频率分析部,其对所述测定数据进行时间频率分析,针对多个规定频率获取表现构成所述信号的该频率成分的强度相对于所述规定物理量的变化的波形数据;
b)区段分割部,其在所述多个规定频率各自的波形数据中,将在所述物理量的变化方向上正值连续的部分和负值连续的部分分别设为一个区段,或者将在所述物理量的变化方向上相邻的极大值与极小值之间的部分设为一个区段,来将该波形数据分割为多个区段;
c)区段值计算部,其针对所述多个规定频率各自的波形数据求出所述多个区段各自的大小;
d)选择区段组制作部,其针对所述多个规定频率各自的波形数据来制作选择区段组,该选择区段组由从该波形数据的多个区段中排除大小超过规定基准值的区段后得到的区段组成;以及
e)噪声水平计算部,其通过求出所述选择区段组中包含的区段的大小的统计值,来求出所述多个规定频率各自的噪声水平。
为了解决上述问题而完成的本发明的第三方式是一种测定数据处理用程序,用于根据通过测定相对于规定物理量发生变化的信号的强度来获取到的包含峰成分和噪声成分的测定数据,来估计噪声成分的大小,该测定数据处理用程序的特征在于,使计算机作为上述第二方式所述的测定数据处理装置而发挥功能。
发明的效果
通过使用本发明所涉及的噪声水平估计方法、测定数据处理装置或者测定数据处理用程序,能够高精度地估计色谱或光谱等测定数据中包含的噪声成分的大小(噪声水平)。
附图说明
图1是峰成分的波形、噪声成分的波形、峰成分与噪声成分之和的波形以及它们的功率谱。
图2是本发明所涉及的测定数据处理装置的一个实施例的结构图。
图3是本实施例的噪声水平估计方法的流程图。
图4是说明对色谱的原始数据(a)进行时间频率分析而获取到的波形数据(b)的图。
图5是说明将波形数据分割为多个区段的工序的图。
图6是说明排除低频侧的区段的工序的图。
图7是对噪声水平的校正进行说明的图。
具体实施方式
下面,参照附图对本发明所涉及的噪声水平估计方法、测定数据处理装置以及测定数据处理用程序的实施例进行说明。在以下的各实施例中,对使用液相色谱仪获取到的色谱中含有的噪声成分的大小即噪声水平进行估计。
图2示出本实施例的测定数据处理装置10的结构。测定数据处理装置10的实体是通用的个人计算机,具备:中央运算装置(CPU)11、存储器12、显示部(监视器)13、输入部14、由硬盘等大容量存储装置构成的存储部16以及通信接口(通信I/F)17等。测定数据处理装置10通过通信I/F 17连接于未图示的液相色谱仪。
在存储部16中保存有在获取色谱数据时一并制作的指标数据。指标数据是例如记录有如梯度分析时的溶液混合比随时间的变化、测定室的室温随时间的变化那样对噪声水平的增减造成影响的参数随时间的变化的数据。另外,在存储部16中保存有OS(OperatingSystem:操作系统)和测定数据处理用程序18,通过启动该测定数据处理用程序18,CPU 11作为后述的时间频率分析部18a、区段分割部18b、区段值计算部18c、选择区段组制作部18d、噪声水平计算部18e以及噪声水平校正部18f而发挥功能。
以下,参照图3的流程图对使用了本实施例的测定数据处理装置10的噪声水平估计方法进行说明。
首先,时间频率分析部18a基于使用者进行了分析频率决定和分析开始的指示,对测定数据进行时间频率分析,针对由使用者指定的多个频率获取表现构成色谱的该频率成分的强度随时间的变化的波形数据(步骤S1)。图4的(b)示出对色谱的一部分(图4的(a))进行mexican-hat连续小波变换所得到的波形数据。mexican-hat连续小波变换是如下一种方法:将mexican-hat小波用作母小波(局部波)来进行频率轴方向上连续的时间频率分析。还能够使用其它小波(例如复数Morlet小波或Haar小波等)作为母小波,或者使用在频率轴方向上离散的时间频率分析即离散小波变换。或者也可以设置将色谱在时间轴方向上分割为多个的窗口,并使用按每个该窗口进行傅立叶变换的滤波器组等。但是,基于能够获取在频率轴方向和时间轴方向上连续的数据这一点,优选使用连续小波变换。
接着,区段分割部18b将在时间轴方向上正值连续的时间段的部分和负值连续的时间段的部分分别设为一个区段,从而将通过时间频率分析获取到的多个频率的波形数据中的每一个波形数据分割为多个区段(步骤S2)。图5示出将一个波形数据分割为多个区段的一例。
接着,区段值计算部18c求出最高频率的波形数据的多个区段各自的面积。然后,基于存储部16中保存的指标数据将面积值进行标准化并设为区段值(步骤S3)。即,利用指标数据中记录的噪声因素随时间的变化来计算将由该噪声因素导致的噪声水平的增减消除后的区段值。并且,选择区段组制作部18d求出属于同一波形数据的多个区段的区段值的平均值和无偏标准偏差σ,从构成波形数据的多个区段中排除区段值超过平均值+Nσ(N为正整数)的区段(图4的(b)),从而制作选择区段组(步骤S4)。
在此,说明了将各区段的面积进行标准化并求出区段值,排除超过该区段值的平均值+Nσ的值的区段的例子,但也能够使用高度来替代面积,或者将面积或高度比中央值+M×MAD(Median Absolute Deviation:绝对中位差)(M为正整数)大的区段排除。此外,对于上式中的N或M的值,能够考虑区段值的分布地使用适于各测定数据的值。
如上所述,如果在色谱中存在峰,则在时间频率分析后的波形中区段的面积和高度变大。因而,能够通过步骤S4从波形数据中排除峰成分。
在存在已被选择区段组制作部18d排除的区段的情况下(在步骤S5中为“是”),排除频率更低的波形数据的在时间轴方向上存在于相同位置的区段(步骤S6)。在不存在已被选择区段组制作部18d排除的区段的情况下(在步骤S5中为“否”),区段值计算部18c再次根据下一个最高频率的波形数据(即,在未处理的波形数据中为最高频率的波形数据)的多个区段各自的面积来求出区段值(步骤S3),与上述同样地制作选择区段组(步骤S4),在存在已被排除的区段的情况下(在步骤S5中为“是”),排除频率更低的波形数据的在时间轴方向上存在于相同位置的区段(参照步骤S6、图6)。
如上所述,在某个频率处存在峰成分的情况下,在比该频率低的频带中必然存在峰成分。因而,通过进行步骤S6,能够更加可靠地排除峰成分来制作选择区段组。
当从高频侧起依次制作选择区段组并针对所有波形数据都制作出选择区段组时(在步骤S7中为“是”),噪声水平计算部18e根据各频率的选择区段组的面积的平均值来求出噪声水平(步骤S8)。
如上所述,普通的检测器中包括的电容器或A/D变换器与低通滤波器同样地发挥功能,因此在由这种检测器获取到的测定数据中减少了高频带的信号。在此,为了反映这种倾向,噪声水平校正部18f确认低频侧的噪声水平是否为高频侧的噪声水平的值以上,在低频侧的噪声水平不为高频侧的噪声水平的值以上的情况下(在步骤S9中为“否”),对噪声水平的计算值进行校正(步骤S10),来决定各频率的噪声水平(步骤S11)。图7示出噪声水平的校正例。
此外,上述步骤S5、S6、S9以及S10是用于高精度地计算噪声水平的附加工序,并非本发明的必要工序。另外,仅在必要的情况下使用指标数据进行区段的面积的标准化即可,在一定的条件下获取到色谱的情况下(即在认为不存在使噪声水平随时间变化的特定因素的情况下),也可以不进行该工序。即,也可以根据各频率的波形数据分别独立地制作选择区段组,将构成该选择区段组的区段的区段值的平均值或中央值直接设为该频率的噪声水平。
在上述实施例中,说明了对由液相色谱仪获取到的色谱进行处理的例子,但除了能够用于对利用液相色谱仪或气相色谱仪获取到的色谱进行处理以外,还能够用于决定通过分光测定获取到的分光光谱、和通过质谱分析获取到的质谱等各种测定数据的噪声水平。
除此以外,在上述实施例中,将正值连续的时间段的部分和负值连续的时间段的部分设为一个区段来对波形数据进行了分割,但也可以将被极大值和极小值夹着的时间段的部分设为一个区段来对波形数据进行分割。
附图标记说明
10:测定数据处理装置;11:CPU;12:存储器;14:输入部;16:存储部;17:通信I/F;18:测定数据处理用程序;18a:时间频率分析部;18b:区段分割部;18c:区段值计算部;18d:选择区段组制作部;18e:噪声水平计算部;18f:噪声水平校正部。
Claims (6)
1.一种噪声水平估计方法,根据通过测定相对于规定物理量发生变化的信号的强度来获取到的包含峰成分和噪声成分的测定数据,来估计噪声成分的大小,该方法的特征在于包括如下步骤:
a)对所述测定数据进行时间频率分析,针对多个规定频率获取表现构成所述信号的频率成分的强度相对于所述规定物理量的变化的波形数据;
b)在所述多个规定频率各自的波形数据中,将在所述物理量的变化方向上正值连续的部分和负值连续的部分分别设为一个区段,或者将在所述物理量的变化方向上相邻的极大值与极小值之间的部分设为一个区段,来将该波形数据分割为多个区段;
c)针对所述多个规定频率各自的波形数据求出所述多个区段各自的大小;
d)针对所述多个规定频率各自的波形数据来制作选择区段组,该选择区段组由从该波形数据的多个区段中排除大小超过规定基准值的区段后得到的区段组成;以及
e)通过求出所述选择区段组中包含的区段的大小的统计值,来求出所述多个规定频率各自的噪声水平。
2.根据权利要求1所述的噪声水平估计方法,其特征在于,还包括如下步骤:
使用所述测定数据中的与所述物理量的变化方向上的噪声成分的大小的变动有关的指标数据,在制作所述选择区段组之前将所述多个区段各自的大小进行标准化。
3.根据权利要求1所述的噪声水平估计方法,其特征在于,还包括如下步骤:
排除如下的区段来制作选择区段组:该区段在所述规定物理量的变化方向上处于与已排除的所述区段对应的位置,并且属于频率比已排除的所述区段所属的波形数据的频率低的波形数据。
4.根据权利要求2所述的噪声水平估计方法,其特征在于,还包括如下步骤:
排除如下的区段来制作选择区段组:该区段在所述规定物理量的变化方向上处于与已排除的所述区段对应的位置,并且属于频率比已排除的所述区段所属的波形数据的频率低的波形数据。
5.根据权利要求1至4中的任一项所述的噪声水平估计方法,其特征在于,还包括如下步骤:
将所述多个规定频率的噪声水平相互比较,以使低频侧的噪声水平成为高频侧的噪声水平以上的值的方式校正所述噪声水平。
6.一种测定数据处理装置,用于根据通过测定相对于规定物理量发生变化的信号的强度来获取到的包含峰成分和噪声成分的测定数据,来估计噪声成分的大小,该测定数据处理装置的特征在于,具备:
a)时间频率分析部,其对所述测定数据进行时间频率分析,针对多个规定频率获取表现构成所述信号的频率成分的强度相对于所述规定物理量的变化的波形数据;
b)区段分割部,其在所述多个规定频率各自的波形数据中,将在所述物理量的变化方向上正值连续的部分和负值连续的部分分别设为一个区段,或者将在所述物理量的变化方向上相邻的极大值与极小值之间的部分设为一个区段,来将该波形数据分割为多个区段;
c)区段值计算部,其针对所述多个规定频率各自的波形数据求出所述多个区段各自的大小;
d)选择区段组制作部,其针对所述多个规定频率各自的波形数据来制作选择区段组,该选择区段组由从该波形数据的多个区段中排除大小超过规定基准值的区段后得到的区段组成;以及
e)噪声水平计算部,其通过求出所述选择区段组中包含的区段的大小的统计值,来求出所述多个规定频率各自的噪声水平。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2015/054068 WO2016132422A1 (ja) | 2015-02-16 | 2015-02-16 | ノイズレベル推定方法、測定データ処理装置、及び測定データ処理用プログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107250788A CN107250788A (zh) | 2017-10-13 |
CN107250788B true CN107250788B (zh) | 2019-06-07 |
Family
ID=56692019
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201580076316.1A Active CN107250788B (zh) | 2015-02-16 | 2015-02-16 | 噪声水平估计方法和测定数据处理装置 |
Country Status (5)
Country | Link |
---|---|
US (1) | US11187685B2 (zh) |
EP (1) | EP3260858A4 (zh) |
JP (1) | JP6332545B2 (zh) |
CN (1) | CN107250788B (zh) |
WO (1) | WO2016132422A1 (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB201502447D0 (en) * | 2015-02-13 | 2015-04-01 | Univ Liverpool | Method and apparatus for sample analysis |
JP7430153B2 (ja) | 2021-03-25 | 2024-02-09 | 株式会社日立製作所 | 電気泳動に係るデータを補正する方法、ピークが試料由来ピークであるかスパイクであるかを判定する方法、装置、およびプログラム |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0296781A2 (en) * | 1987-06-19 | 1988-12-28 | Applied Biosystems, Inc. | Quantitation of chromatographic information |
CN101002460A (zh) * | 2004-06-23 | 2007-07-18 | 兰达解决方案公司 | 输入数据的自动背景去除 |
US7340375B1 (en) * | 1999-02-15 | 2008-03-04 | Her Majesty The Queen In Right Of Canada, As Represented By The Minister Of Industry Through The Communications Research Centre | Method and apparatus for noise floor estimation |
CN102472732A (zh) * | 2009-07-31 | 2012-05-23 | 西门子公司 | 色谱滤波方法 |
CN102998296A (zh) * | 2012-11-28 | 2013-03-27 | 重庆绿色智能技术研究院 | 一种去除背景噪声影响的拉曼光谱预处理方法 |
Family Cites Families (43)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS62291562A (ja) * | 1986-06-11 | 1987-12-18 | Hitachi Ltd | クロマトグラフのデ−タ処理方法 |
JPH0798270A (ja) | 1993-09-29 | 1995-04-11 | Shimadzu Corp | フロースルー型紫外可視分光光度計 |
JPH08327612A (ja) * | 1995-05-30 | 1996-12-13 | Sumitomo Metal Ind Ltd | 波形信号処理方法 |
US5966684A (en) * | 1995-11-09 | 1999-10-12 | Richardson; Thomas L. | Method and apparatus for cancelling periodic electrical interference |
US5748491A (en) * | 1995-12-20 | 1998-05-05 | The Perkin-Elmer Corporation | Deconvolution method for the analysis of data resulting from analytical separation processes |
US6112161A (en) * | 1997-09-17 | 2000-08-29 | Hewlett-Packard | Method, apparatus, and article of manufacture for enhanced intergration of signals |
US6519486B1 (en) * | 1998-10-15 | 2003-02-11 | Ntc Technology Inc. | Method, apparatus and system for removing motion artifacts from measurements of bodily parameters |
US6449584B1 (en) * | 1999-11-08 | 2002-09-10 | Université de Montréal | Measurement signal processing method |
US7409298B2 (en) * | 2002-04-12 | 2008-08-05 | Northeastern University | Matched filtration with experimental noise determination for denoising, peak picking and quantitation in LC-MS |
US20040204866A1 (en) * | 2003-04-09 | 2004-10-14 | Allington Robert W. | Method and apparatus to enhance the signal to noise ratio in chromatography |
GB2410123B (en) * | 2003-09-25 | 2006-11-01 | Thermo Finnigan Llc | Method of processing and storing mass spectrometry data |
JP4726900B2 (ja) * | 2004-07-27 | 2011-07-20 | ナティビス インコーポレイテッド | 非常に低振幅な信号を収集、記憶、処理、伝送および提示するためのシステムおよび方法 |
GB2472951B (en) * | 2004-11-29 | 2011-04-27 | Thermo Finnigan Llc | Method of processing mass spectrometry data |
JP2006163614A (ja) | 2004-12-03 | 2006-06-22 | Matsushita Electric Ind Co Ltd | 画像処理装置、および画像処理方法 |
US7546208B2 (en) * | 2005-03-07 | 2009-06-09 | Mks Instruments, Inc. | Method and apparatus of signal processing for use in spectrometry using an improved apodization function |
US20060241937A1 (en) * | 2005-04-21 | 2006-10-26 | Ma Changxue C | Method and apparatus for automatically discriminating information bearing audio segments and background noise audio segments |
JP4614002B2 (ja) * | 2006-06-08 | 2011-01-19 | 株式会社島津製作所 | クロマトグラフ質量分析用データ処理装置 |
US8530828B2 (en) * | 2007-02-02 | 2013-09-10 | Dh Technologies Development Pte. Ltd. | Systems and methods for reducing noise from mass spectra |
JP2009008582A (ja) * | 2007-06-29 | 2009-01-15 | Shimadzu Corp | クロマトグラムデータ処理装置 |
US8352274B2 (en) * | 2007-09-11 | 2013-01-08 | Panasonic Corporation | Sound determination device, sound detection device, and sound determination method for determining frequency signals of a to-be-extracted sound included in a mixed sound |
US8027743B1 (en) * | 2007-10-23 | 2011-09-27 | Adobe Systems Incorporated | Adaptive noise reduction |
JP5211753B2 (ja) * | 2008-02-27 | 2013-06-12 | 株式会社島津製作所 | クロマトグラフ用データ処理装置 |
WO2010038385A1 (ja) * | 2008-09-30 | 2010-04-08 | パナソニック株式会社 | 音判定装置、音判定方法、及び、音判定プログラム |
WO2010038386A1 (ja) * | 2008-09-30 | 2010-04-08 | パナソニック株式会社 | 音判定装置、音検知装置及び音判定方法 |
US7983852B2 (en) * | 2008-10-21 | 2011-07-19 | Thermo Finnigan Llc | Methods of automated spectral peak detection and quantification without user input |
US20100283785A1 (en) * | 2009-05-11 | 2010-11-11 | Agilent Technologies, Inc. | Detecting peaks in two-dimensional signals |
EP2322922B1 (en) * | 2009-08-26 | 2015-02-25 | Thermo Fisher Scientific (Bremen) GmbH | Method of improving the resolution of compounds eluted from a chromatography device |
US8428889B2 (en) * | 2010-10-07 | 2013-04-23 | Thermo Finnigan Llc | Methods of automated spectral peak detection and quantification having learning mode |
JP2012177568A (ja) * | 2011-02-25 | 2012-09-13 | Arkray Inc | データ処理装置、データ処理方法、及びデータ処理プログラム |
US20130131998A1 (en) * | 2011-11-18 | 2013-05-23 | David A. Wright | Methods and Apparatus for Identifying Mass Spectral Isotope Patterns |
US9177559B2 (en) * | 2012-04-24 | 2015-11-03 | Tom Stephenson | Method and apparatus for analyzing animal vocalizations, extracting identification characteristics, and using databases of these characteristics for identifying the species of vocalizing animals |
JP6109927B2 (ja) * | 2012-05-04 | 2017-04-05 | カオニックス ラブス リミテッド ライアビリティ カンパニー | 源信号分離のためのシステム及び方法 |
US10497381B2 (en) * | 2012-05-04 | 2019-12-03 | Xmos Inc. | Methods and systems for improved measurement, entity and parameter estimation, and path propagation effect measurement and mitigation in source signal separation |
US10840073B2 (en) * | 2012-05-18 | 2020-11-17 | Thermo Fisher Scientific (Bremen) Gmbh | Methods and apparatus for obtaining enhanced mass spectrometric data |
US9760540B2 (en) * | 2012-11-21 | 2017-09-12 | National Central University | Methods for processing sequential data to identify possible peak points and to estimate peak to noise ratio of sequential data |
EP2747081A1 (en) * | 2012-12-18 | 2014-06-25 | Oticon A/s | An audio processing device comprising artifact reduction |
JP2014123011A (ja) * | 2012-12-21 | 2014-07-03 | Sony Corp | 雑音検出装置および方法、並びに、プログラム |
JP2014137350A (ja) * | 2013-01-18 | 2014-07-28 | Hitachi High-Technologies Corp | クロマトグラムのノイズ除去方法およびクロマトグラフィー装置 |
US20140252218A1 (en) * | 2013-03-05 | 2014-09-11 | David A. Wright | Methods and Apparatus for Decomposing Tandem Mass Spectra Generated by All-Ions Fragmentation |
CN105518455B (zh) * | 2013-09-09 | 2017-05-03 | 株式会社岛津制作所 | 波峰检测方法 |
JP6256162B2 (ja) * | 2014-04-07 | 2018-01-10 | 株式会社島津製作所 | 信号波形データ処理装置 |
CN107076712B (zh) * | 2014-09-03 | 2019-01-11 | 株式会社岛津制作所 | 色谱数据处理方法以及装置 |
JP6458812B2 (ja) * | 2015-01-26 | 2019-01-30 | 株式会社島津製作所 | ノイズレベル推定方法、測定データ処理装置、及び測定データ処理用プログラム |
-
2015
- 2015-02-16 EP EP15882527.3A patent/EP3260858A4/en not_active Withdrawn
- 2015-02-16 US US15/549,432 patent/US11187685B2/en active Active
- 2015-02-16 JP JP2017500489A patent/JP6332545B2/ja active Active
- 2015-02-16 CN CN201580076316.1A patent/CN107250788B/zh active Active
- 2015-02-16 WO PCT/JP2015/054068 patent/WO2016132422A1/ja active Application Filing
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0296781A2 (en) * | 1987-06-19 | 1988-12-28 | Applied Biosystems, Inc. | Quantitation of chromatographic information |
US7340375B1 (en) * | 1999-02-15 | 2008-03-04 | Her Majesty The Queen In Right Of Canada, As Represented By The Minister Of Industry Through The Communications Research Centre | Method and apparatus for noise floor estimation |
CN101002460A (zh) * | 2004-06-23 | 2007-07-18 | 兰达解决方案公司 | 输入数据的自动背景去除 |
CN102472732A (zh) * | 2009-07-31 | 2012-05-23 | 西门子公司 | 色谱滤波方法 |
CN102998296A (zh) * | 2012-11-28 | 2013-03-27 | 重庆绿色智能技术研究院 | 一种去除背景噪声影响的拉曼光谱预处理方法 |
Also Published As
Publication number | Publication date |
---|---|
JP6332545B2 (ja) | 2018-05-30 |
US11187685B2 (en) | 2021-11-30 |
EP3260858A1 (en) | 2017-12-27 |
CN107250788A (zh) | 2017-10-13 |
US20180003683A1 (en) | 2018-01-04 |
EP3260858A4 (en) | 2017-12-27 |
JPWO2016132422A1 (ja) | 2017-07-13 |
WO2016132422A1 (ja) | 2016-08-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10416134B2 (en) | Chromatogram data processing method and chromatogram data processing apparatus | |
US10121643B2 (en) | Chromatography/mass spectrometry data processing device | |
US10935526B2 (en) | Liquid chromatograph and method for correcting detector output value fluctuation of liquid chromatograph | |
CN109863558A (zh) | 质谱数据的评估方法和质谱法以及maldi tof质谱仪 | |
CN106018331B (zh) | 多通道光谱系统的稳定性评价方法及预处理优化方法 | |
CN107250788B (zh) | 噪声水平估计方法和测定数据处理装置 | |
CN106643835A (zh) | 光纤法珀腔解调方法及装置和光纤法珀干涉仪 | |
US9285444B2 (en) | Method for substance identification from NMR spectrum | |
CN104990895A (zh) | 一种基于局部区域的近红外光谱信号标准正态校正方法 | |
US10359404B2 (en) | Noise level estimation method, measurement data processing device and program for processing measurement data | |
CN104155262A (zh) | 一种烟草水溶性糖近红外定量模型中光谱范围的选择方法 | |
CN102472732B (zh) | 色谱滤波方法 | |
US20080073499A1 (en) | Peak finding in low-resolution mass spectrometry by use of chromatographic integration routines | |
EP2483641A1 (en) | Systems and methods for maintaining the precision of mass measurement | |
US10236167B1 (en) | Peak waveform processing device | |
US20230296572A1 (en) | Training Method | |
CN109324017A (zh) | 一种提高近红外光谱分析技术建模光谱质量的方法 | |
US11913919B2 (en) | Techniques for acquisition-stage peak width determination and baseline offset estimation | |
CN115219618B (zh) | 一种噪声处理方法、装置、液质联用仪及可读存储介质 | |
CN114609081A (zh) | 一种基于多组分样品的多传感器光谱数据处理方法 | |
CN116306960A (zh) | 一种光谱的任意阶微积分方法 | |
CN117724037A (zh) | 一种基于相关干涉的ads_b信号的测向方法 | |
Rasaruddin et al. | Signal processing strategies in FT-NIR and FTIR spectra of palm oils | |
CN111256819A (zh) | 一种光谱仪器的降噪方法 | |
CN116660207A (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 |