CN108535403B - 数据处理的方法及装置 - Google Patents

数据处理的方法及装置 Download PDF

Info

Publication number
CN108535403B
CN108535403B CN201810346349.9A CN201810346349A CN108535403B CN 108535403 B CN108535403 B CN 108535403B CN 201810346349 A CN201810346349 A CN 201810346349A CN 108535403 B CN108535403 B CN 108535403B
Authority
CN
China
Prior art keywords
ion
parent
daughter ion
specified
specified daughter
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
CN201810346349.9A
Other languages
English (en)
Other versions
CN108535403A (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.)
Southwest University
Original Assignee
Southwest University
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 Southwest University filed Critical Southwest University
Priority to CN201810346349.9A priority Critical patent/CN108535403B/zh
Publication of CN108535403A publication Critical patent/CN108535403A/zh
Application granted granted Critical
Publication of CN108535403B publication Critical patent/CN108535403B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/96—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 using ion-exchange
    • 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/8603—Signal analysis with integration or differentiation
    • G01N2030/862—Other mathematical operations for data preprocessing

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Mathematical Physics (AREA)
  • Algebra (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Abstract

本发明提供了一种数据处理的方法及装置，涉及质谱分析领域。该方法包括：获取指定子离子和所有母离子，对指定子离子进行小波处理得到降噪后的指定子离子色谱图，基于指定子离子色谱图确定指定子离子的出峰时间，对所有母离子进行小波处理得到降噪后的所有母离子时间序列，基于所有母离子时间序列确定所有母离子的出峰时间，基于指定子离子的出峰时间和所有母离子的出峰时间获取出峰时间位于指定子离子出峰时间预设范围内的多个备选母离子，计算指定子离子和多个备选母离子的小波相关性并生成小波相关性数据，基于小波相关性数据判断小波相关性数据高于预设数据的母离子，以获取与指定子离子对应的母离子。本发明能实现查找数据中子离子对应的母离子。

Description

数据处理的方法及装置
技术领域
本发明涉及质谱分析领域，具体而言，涉及一种数据处理的方法及装置。
背景技术
质谱分析能用于检测化合物的灵敏性和诊断化合物的准确性，使得其被广泛应用于生物，医药，环境，和化学等领域。硬件的革新使得质谱中质量分析器的串联方式增多，质谱分辨率大幅提升，但同时也使得研究者对处理质谱数据的需求和要求不断增加。
四级杆串联时间飞行质谱仪是一种串联质谱。它的其中一种工作方式是，在四级杆筛选完母离子后，母离子会进入一个反应室被裂解成为子离子，所有的子离子进入到时间飞行仪中并被记录质量信息。这种质谱数据采集方式被称为独立于信息的数据采集方式(DIA)。DIA的优点在于它能够记录所有子离子的信息，但缺点在于不能直接得知母离子和子离子的一一对应关系。
Sequential Windowed Acquisition of all Theoretical fragment ions(SWATH)是DIA中的一种，主要被应用于蛋白质组学。在进行蛋白质组学的分析时，为了克服母离子和子离子的对应关系问题，研究者会事先建立好标准蛋白质的质谱库，通过将SWATH数据与标准库比对，可得出有哪些蛋白质存在于样品中。而对于没有标准品的化合物，没有数据库可以比对，从而无法鉴定样品中存在的物质。
发明内容
有鉴于此，本发明实施例的目的在于提供一种数据处理的方法及装置，以实现查找数据中子离子对应的母离子。
第一方面，本发明实施例提供了一种数据处理的方法，所述方法包括：获取指定子离子和所有母离子；对所述指定子离子进行小波处理得到降噪后的指定子离子色谱图，并基于所述指定子离子色谱图确定所述指定子离子的出峰时间；对所述所有母离子进行小波处理得到降噪后的所有母离子时间序列，并基于所述所有母离子时间序列确定所述所有母离子的出峰时间；基于所述指定子离子的出峰时间和所述所有母离子的出峰时间获取出峰时间位于所述指定子离子的出峰时间预设范围内的多个备选母离子；计算所述指定子离子和所述多个备选母离子的小波相关性并生成小波相关性数据；基于所述小波相关性数据判断所述小波相关性数据高于预设数据的母离子，以获取与所述指定子离子对应的母离子。
在本发明较佳的实施例中，上述获取指定子离子和所有母离子之前，还包括：采集荷质比为1-1000的多个母离子和荷质比为1-1000且只带一个电荷的多个子离子；将所述多个子离子和所述多个母离子以预定格式进行存储。
在本发明较佳的实施例中，上述对所述指定子离子进行小波处理得到降噪后的指定子离子色谱图，并基于所述指定子离子色谱图确定所述指定子离子的出峰时间，包括：对所述指定子离子进行小波处理得到降噪后的指定子离子色谱图，并基于所述指定子离子色谱图确定所述指定子离子的出峰个数和每个峰出峰时间；根据所述出峰个数确定所述能生成所述指定子离子的母离子个数，其中，所述出峰个数与所述能生成所述指定子离子的母离子个数相等。
在本发明较佳的实施例中，上述基于所述小波相关性数据判断所述小波相关性数据高于预设数据的母离子，以获取与所述指定子离子对应的母离子之后，还包括：根据所述指定子离子和与所述指定子离子对应的母离子，绘制所述指定子离子和与所述指定子离子对应的母离子的色谱图、小波图以及小波相关性图。
在本发明较佳的实施例中，上述方法还包括：显示所述指定子离子色谱图和所述所有母离子色谱图，以供用户直观地看到所述指定子离子色谱图、所述所有母离子色谱图以及所述指定子离子色谱图的变化情况。
第二方面，本发明实施例提供了一种数据处理的装置，所述装置包括：子母离子获取模块，用于获取指定子离子和所有母离子；第一出峰时间获取模块，用于对所述指定子离子进行小波处理得到降噪后的指定子离子色谱图，并基于所述指定子离子色谱图确定所述指定子离子的出峰时间；第二出峰时间获取模块，用于对所述所有母离子进行小波处理得到降噪后的所有母离子时间序列，并基于所述所有母离子时间序列确定所述所有母离子的出峰时间；备选母离子获取模块，用于基于所述指定子离子的出峰时间和所述所有母离子的出峰时间获取出峰时间位于所述指定子离子的出峰时间预设范围内的多个备选母离子；相关性数据生成模块，用于计算所述指定子离子和所述多个备选母离子的小波相关性并生成小波相关性数据；对应母离子获取模块，用于基于所述小波相关性数据判断所述小波相关性数据高于预设数据的母离子，以获取与所述指定子离子对应的母离子。
在本发明较佳的实施例中，上述装置还包括：子母离子采集模块，用于采集荷质比为1-1000的多个母离子和荷质比为1-1000且只带一个电荷的多个子离子；子母离子存储模块，用于将所述多个子离子和所述多个母离子以预定格式进行存储。
在本发明较佳的实施例中，上述第一出峰时间获取模块包括：出峰个数获取子模块，用于对所述指定子离子进行小波处理得到降噪后的指定子离子色谱图，并基于所述指定子离子色谱图确定所述指定子离子的出峰个数和每个峰出峰时间；母离子个数确定子模块，用于根据所述出峰个数确定所述能生成所述指定子离子的母离子个数，其中，所述出峰个数与所述能生成所述指定子离子的母离子个数相等。
在本发明较佳的实施例中，上述装置还包括：图形绘制模块，用于根据所述指定子离子和与所述指定子离子对应的母离子，绘制所述指定子离子和与所述指定子离子对应的母离子的色谱图、小波图以及小波相关性图。
在本发明较佳的实施例中，上述装置还包括：显示模块，用于显示所述指定子离子色谱图和所述所有母离子色谱图，以供用户直观地看到所述指定子离子色谱图、所述所有母离子色谱图以及所述指定子离子色谱图的变化情况。
与现有技术相比，本发明各实施例提供的数据处理的方法及装置首先获取指定子离子和所有母离子，对指定子离子进行小波处理得到降噪后的指定子离子色谱图，并基于指定子离子色谱图确定指定子离子的出峰时间，对所有母离子进行小波处理得到降噪后的所有母离子时间序列，并基于所有母离子时间序列确定所有母离子的出峰时间，然后基于指定子离子的出峰时间和所有母离子的出峰时间获取出峰时间位于指定子离子出峰时间预设范围内的多个备选母离子，计算指定子离子和多个备选母离子的小波相关性并生成小波相关性数据，最后基于小波相关性数据判断小波相关性数据高于预设数据的母离子，以获取与指定子离子对应的母离子，从而能实现查找数据中子离子对应的母离子。
为使本发明的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
图1为本发明实施例提供的电子设备的结构框图；
图2为本发明第一实施例提供的数据处理的方法的流程示意图；
图3为本发明第二实施例提供的数据处理的方法的流程示意图；
图4为本发明第三实施例提供的数据处理的装置的结构框图。
具体实施方式
下面将结合本发明实施例中附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。因此，以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施例。基于本发明的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。同时，在本发明的描述中，术语“第一”、“第二”等仅用于区分描述，而不能理解为指示或暗示相对重要性。
如图1所示，是电子设备100的方框示意图。所述电子设备100包括：数据处理的装置、存储器110、存储控制器120、处理器130、外设接口140、输入输出单元150、音频单元160、显示单元170。
所述存储器110、存储控制器120、处理器130、外设接口140、输入输出单元150、音频单元160以及显示单元170各元件相互之间直接或间接地电性连接，以实现数据的传输或交互。例如，这些元件相互之间可通过一条或多条通讯总线或信号线实现电性连接。所述数据处理的装置包括至少一个可以软件或固件(firmware)的形式存储于所述存储器中或固化在所述客户端设备的操作***(operating system，OS)中的软件功能模块。所述处理器130用于执行存储器110中存储的可执行模块，例如所述数据处理的装置包括的软件功能模块或计算机程序。
其中，存储器110可以是，但不限于，随机存取存储器(Random Access Memory，RAM)，只读存储器(Read Only Memory，ROM)，可编程只读存储器(Programmable Read-OnlyMemory，PROM)，可擦除只读存储器(Erasable Programmable Read-Only Memory，EPROM)，电可擦除只读存储器(Electric Erasable Programmable Read-Only Memory，EEPROM)等。其中，存储器110用于存储程序，所述处理器130在接收到执行指令后，执行所述程序，前述本发明实施例任一实施例揭示的流过程定义的服务器所执行的方法可以应用于处理器130中，或者由处理器130实现。
处理器130可能是一种集成电路芯片，具有信号的处理能力。上述的处理器130可以是通用处理器，包括中央处理器(Central Processing Unit，简称CPU)、网络处理器(Network Processor，简称NP)等；还可以是数字信号处理器(DSP)、专用集成电路(ASIC)、现成可编程门阵列(FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本发明实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
所述外设接口140将各种输入/输入装置耦合至处理器130以及存储器110。在一些实施例中，外设接口140，处理器130以及存储控制器120可以在单个芯片中实现。在其他一些实例中，他们可以分别由独立的芯片实现。
输入输出单元150用于提供给用户输入数据实现用户与电子设备100的交互。所述输入输出单元150可以是，但不限于，鼠标和键盘等。
音频单元160向用户提供音频接口，其可包括一个或多个麦克风、一个或者多个扬声器以及音频电路。
显示单元170在电子设备100与用户之间提供一个交互界面(例如用户操作界面)或用于显示图像数据给用户参考。在本实施例中，所述显示单元170可以是液晶显示器或触控显示器。若为触控显示器，其可为支持单点和多点触控操作的电容式触控屏或电阻式触控屏等。支持单点和多点触控操作是指触控显示器能感应到来自该触控显示器上一个或多个位置处同时产生的触控操作，并将该感应到的触控操作交由处理器130进行计算和处理。
第一实施例
请参照图2，图2是本发明第一实施例提供的一种数据处理的方法的流程图。下面将对图2所示的流程进行详细阐述，所述方法包括：
步骤S101：获取指定子离子和所有母离子。
作为本实施例的一种实施方式，在电子设备中安装有小波寻找质谱母离子数据分析软件，该软件用于对SWATH数据进行分析。可以理解的，SWATH数据包括多个子离子和多个母离子，获取指定子离子和所有母离子，即该软件读取SWATH数据中的指定子离子和所有母离子。
步骤S102：对所述指定子离子进行小波处理得到降噪后的指定子离子色谱图，并基于所述指定子离子色谱图确定所述指定子离子的出峰时间。
在本实施例中，在获取指定子离子后，对该指定子离子进行小波处理得到降噪后的指定子离子色谱图，优选的，采用mexican hat小波处理该指定子离子的信号，得到降噪后的指定子离子色谱图。
进一步的，基于该指定子离子色谱图确定该指定子离子的出峰时间，可以理解的，当该指定子离子的出峰时间为多个时，分别计算每个峰的出峰时间，以确定每个峰的出峰时间。
步骤S103：对所述所有母离子进行小波处理得到降噪后的所有母离子时间序列，并基于所述所有母离子时间序列确定所述所有母离子的出峰时间。
同理，在获取所有母离子后，对所有母离子进行小波处理得到降噪后的所有母离子时间序列，优选的，采用mexican hat小波处理该所有母离子的信号，得到降噪后的所有母离子时间序列。
进一步的，基于该所有母离子时间序列确定所有母离子的出峰时间，同样的，分别计算所有母离子中每个母离子的出峰时间，以确定所有母离子的出峰时间。
作为本发明实施例优选的一种方式，在获取子离子色谱图和母离子时间序列时，显示该指定子离子色谱图，以供用户直观地看到指定子离子色谱图，同时，在用户调节参数，优化指定子离子的峰形的选峰数量时，显示子离子色谱图的变化情况，实现可视化。
步骤S104：基于所述指定子离子的出峰时间和所述所有母离子的出峰时间获取出峰时间位于所述指定子离子的出峰时间预设范围内的多个备选母离子。
其中，在确定指定子离子的出峰时间和所有母离子的出峰时间后，对指定子离子的出峰时间和所有母离子的出峰时间进行判断，以获取出峰时间相近的指定子离子和其对应的母离子。优选的，在本实施例中，基于指定子离子的出峰时间和所有母离子的出峰时间获取出峰时间位于该指定子离子的出峰时间预设范围内的多个备选母离子。
具体的，作为一种方式，预先设置一个时间范围，当确定指定子离子的出峰时间和所有母离子的出峰时间后，计算指定子离子的出峰时间和所有母离子中每个母离子的出峰时间的差值，当差值位于预先设置的时间范围之外时，判定该指定子离子的出峰时间和所有母离子中的该母离子的出峰时间不相近；当该差值位于预先设置的时间范围之内时，判定该指定子离子的出峰时间和所有母离子中的该母离子的出峰时间相近，将差值位于预先设置的时间范围内的母离子作为备选母离子，从而获取多个备选母离子。
步骤S105：计算所述指定子离子和所述多个备选母离子的小波相关性并生成小波相关性数据。
进一步的，计算指定子离子和多个备选母离子的小波相关性，根据计算结果生成小波相关性数据，作为一种方式，采用morlet小波计算指定子离子和多个备选母离子的小波相关性。
步骤S106：基于所述小波相关性数据判断所述小波相关性数据高于预设数据的母离子，以获取与所述指定子离子对应的母离子。
其中，作为一种方式，预先设置一个阙值，将该阙值作为预设数据，然后基于小波相关性数据和预设数据对备选母离子进行筛选，当小波相关性数据高于预设数据时，则判定该小波相关性数据对应的母离子为该指定子离子对应的母离子。进一步的，将该数据报告进行存储。
本发明第一实施例提供的数据处理的方法首先获取指定子离子和所有母离子，对指定子离子进行小波处理得到降噪后的指定子离子色谱图，并基于指定子离子色谱图确定指定子离子的出峰时间，对所有母离子进行小波处理得到降噪后的所有母离子时间序列，并基于所有母离子时间序列确定所有母离子的出峰时间，然后基于指定子离子的出峰时间和所有母离子的出峰时间获取出峰时间位于指定子离子出峰时间预设范围内的多个备选母离子，计算指定子离子和多个备选母离子的小波相关性并生成小波相关性数据，最后基于小波相关性数据判断小波相关性数据高于预设数据的母离子，以获取与指定子离子对应的母离子，从而能实现查找数据中子离子对应的母离子。
第二实施例
请参照图3，图3是本发明第二实施例提供的一种数据处理的方法的流程图。下面将对图3所示的流程进行详细阐述，所述方法包括：
步骤S201：采集荷质比为1-1000的多个母离子和荷质比为1-1000且只带一个电荷的多个子离子。
作为本本实施例的一种实施方式，在质谱采集数据的时候，选取质荷比(m/z)的范围为1-1000的，多个母离子和只带一个电荷的多个子离子进行扫描，从而获取小分子的多个子离子和多个母离子，使得该软件能针对小分子的子离子和母离子的对应关系进行查找。
需要说明的是，本实施例没有对同位素峰的分析，避免了因为分析同位素峰而产生的错误。
步骤S202：将所述多个子离子和所述多个母离子以预定格式进行存储。
将采集的多个子离子和多个母离子作为SWATH数据导入该小波寻找质谱母离子数据分析软件，并将多个子离子和多个母离子的数据以预定格式进行存储。
步骤S203：获取指定子离子和所有母离子。
步骤S204：对所述指定子离子进行小波处理得到降噪后的指定子离子色谱图，并基于所述指定子离子色谱图确定所述指定子离子的出峰个数和每个峰出峰时间。
作为一种方式，在基于该指定子离子色谱图确定指定子离子的每个峰的出峰时间的同时，还基于该指定子离子色谱图确定指定子离子的出峰个数，可以理解的，该出峰个数可以为一个、可以为两个、可以为三个，也可以为多个等等，在此不做限定。
步骤S205：根据所述出峰个数确定所述能生成所述指定子离子的母离子个数，其中，所述出峰个数与所述能生成所述指定子离子的母离子个数相等。
在获取指定子离子的出峰个数后，根据该出峰个数确定能生成指定子离子的母离子个数，其中，该出峰个数与能生成指定子离子的母离子个数相等，即一个峰代表一个母离子。
步骤S206：对所述所有母离子进行小波处理得到降噪后的所有母离子时间序列，并基于所述所有母离子时间序列确定所述所有母离子的出峰时间。
步骤S207：基于所述指定子离子的出峰时间和所述所有母离子的出峰时间获取出峰时间位于所述指定子离子的出峰时间预设范围内的多个备选母离子。
步骤S208：计算所述指定子离子和所述多个备选母离子的小波相关性并生成小波相关性数据。
步骤S209：基于所述小波相关性数据判断所述小波相关性数据高于预设数据的母离子，以获取与所述指定子离子对应的母离子。
步骤S210：根据所述指定子离子和与所述指定子离子对应的母离子，绘制所述指定子离子和与所述指定子离子对应的母离子的色谱图、小波图以及小波相关性图。
其中，在获取与指定子离子对应的母离子之后，根据该指定子离子和与该指定子离子对应的母离子，绘制指定子离子和与指定子离子对应的母离子的色谱图、小波图以及小波相关性图。
本发明第二实施例提供的数据处理的方法首先采集荷质比为1-1000的多个母离子和荷质比为1-1000且只带一个电荷的多个子离子，将多个子离子和多个母离子以预定格式进行存储，然后获取指定子离子和所有母离子，对指定子离子进行小波处理得到降噪后的指定子离子色谱图，并基于指定子离子色谱图确定指定子离子的出峰个数和每个峰的出峰时间，根据出峰个数确定能成圣该指定子离子的母离子个数，再对母离子进行小波处理得到降噪后的所有母离子时间序列，并基于所有母离子时间序列确定所有母离子的出峰时间，基于指定子离子的出峰时间和所有母离子的出峰时间获取出峰时间位于指定子离子出峰时间预设范围内的多个备选母离子，最后计算指定子离子和多个备选母离子的小波相关性并生成小波相关性数据，基于小波相关性数据判断小波相关性数据高于预设数据的母离子，以获取与指定子离子对应的母离子，绘制指定子离子和与指定子离子对应的母离子的色谱图、小波图以及小波相关性图，从而现实查找数据中小分子的子离子对应的母离子。
第三实施例
请参照图4，图4是本发明第三实施例提供的一种数据处理的装置200的结构框图。下面将对图4所示的结构框图进行阐述，所示数据处理的装置200包括：子母离子采集模块210、子母离子存储模块220、子母离子获取模块230、第一出峰时间获取模块240、第二出峰时间获取模块250、备选母离子获取模块260、相关性数据生成模块270、对应母离子获取模块280、图形绘制模块290以及显示模块300，其中：
子母离子采集模块210，采集荷质比为1-1000的多个母离子和荷质比为1-1000且只带一个电荷的多个子离子。
子母离子存储模块220，用于将所述多个子离子和所述多个母离子以预定格式进行存储。
子母离子获取模块230，用于获取指定子离子和所有母离子。
第一出峰时间获取模块240，用于对所述指定子离子进行小波处理得到降噪后的指定子离子色谱图，并基于所述指定子离子色谱图确定所述指定子离子的出峰时间。作为一种方式，该第一出峰时间获取模块240包括：出峰个数获取子模块242和母离子个数确定子模块244，其中：
出峰个数获取子模块242，用于对所述指定子离子进行小波处理得到降噪后的指定子离子色谱图，并基于所述指定子离子色谱图确定所述指定子离子的出峰个数和每个峰出峰时间。
母离子个数确定子模块244，用于根据所述出峰个数确定所述能生成所述指定子离子的母离子个数，其中，所述出峰个数与所述能生成所述指定子离子的母离子个数相等。
第二出峰时间获取模块250，用于对所述所有母离子进行小波处理得到降噪后的所有母离子时间序列，并基于所述所有母离子时间序列确定所述所有母离子的出峰时间。
备选母离子获取模块260，用于基于所述指定子离子的出峰时间和所述所有母离子的出峰时间获取出峰时间位于所述指定子离子的出峰时间预设范围内的多个备选母离子。
相关性数据生成模块270，用于计算所述指定子离子和所述多个备选母离子的小波相关性并生成小波相关性数据。
对应母离子获取模块280，用于基于所述小波相关性数据判断所述小波相关性数据高于预设数据的母离子，以获取与所述指定子离子对应的母离子。
图形绘制模块290，用于根据所述指定子离子和与所述指定子离子对应的母离子，绘制所述指定子离子和与所述指定子离子对应的母离子的色谱图、小波图以及小波相关性图。
显示模块300，用于显示所述指定子离子色谱图和所述所有母离子色谱图，以供用户直观地看到所述指定子离子色谱图、所述所有母离子色谱图以及所述指定子离子色谱图和所述所有母离子色谱图的变化情况。
本实施例对数据处理的装置200的各功能模块实现各自功能的过程，请参见上述图1至图3所示实施例中描述的内容，此处不再赘述。
综上所述，本发明实施例提供的数据处理的方法及装置首先获取指定子离子和所有母离子，对指定子离子进行小波处理得到降噪后的指定子离子色谱图，并基于指定子离子色谱图确定指定子离子的出峰时间，对所有母离子进行小波处理得到降噪后的所有母离子时间序列，并基于所有母离子时间序列确定所有母离子的出峰时间，然后基于指定子离子的出峰时间和所有母离子的出峰时间获取出峰时间位于指定子离子出峰时间预设范围内的多个备选母离子，计算指定子离子和多个备选母离子的小波相关性并生成小波相关性数据，最后基于小波相关性数据判断小波相关性数据高于预设数据的母离子，以获取与指定子离子对应的母离子，从而能实现查找数据中子离子对应的母离子。
在本申请所提供的几个实施例中，应该理解到，所揭露的装置和方法，也可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的，例如，附图中的流程图和框图显示了根据本发明的多个实施例的装置、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分，所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现方式中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的***来实现，或者可以用专用硬件与计算机指令的组合来实现。
另外，在本发明各个实施例中的各功能模块可以集成在一起形成一个独立的部分，也可以是各个模块单独存在，也可以两个或两个以上模块集成形成一个独立的部分。
所述功能如果以软件功能模块的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。

Claims (10)

1.一种数据处理的方法，其特征在于，所述方法包括：
获取指定子离子和所有母离子；
对所述指定子离子进行小波处理得到降噪后的指定子离子色谱图，并基于所述指定子离子色谱图确定所述指定子离子的出峰时间；
对所述所有母离子进行小波处理得到降噪后的所有母离子时间序列，并基于所述所有母离子时间序列确定所述所有母离子的出峰时间；
基于所述指定子离子的出峰时间和所述所有母离子的出峰时间获取出峰时间位于所述指定子离子的出峰时间预设范围内的多个备选母离子；
计算所述指定子离子和所述多个备选母离子的小波相关性并生成小波相关性数据；
基于所述小波相关性数据判断所述小波相关性数据高于预设数据的母离子，以获取与所述指定子离子对应的母离子。
2.根据权利要求1所述的方法，其特征在于，所述获取指定子离子和所有母离子之前，还包括：
采集荷质比为1-1000的多个母离子和荷质比为1-1000且只带一个电荷的多个子离子；
将所述多个子离子和所述多个母离子以预定格式进行存储。
3.根据权利要求2所述的方法，其特征在于，所述对所述指定子离子进行小波处理得到降噪后的指定子离子色谱图，并基于所述指定子离子色谱图确定所述指定子离子的出峰时间，包括：
对所述指定子离子进行小波处理得到降噪后的指定子离子色谱图，并基于所述指定子离子色谱图确定所述指定子离子的出峰个数和每个峰出峰时间；
根据所述出峰个数确定能生成所述指定子离子的母离子个数，其中，所述出峰个数与所述能生成所述指定子离子的母离子个数相等。
4.根据权利要求3所述的方法，其特征在于，所述基于所述小波相关性数据判断所述小波相关性数据高于预设数据的母离子，以获取与所述指定子离子对应的母离子之后，还包括：
根据所述指定子离子和与所述指定子离子对应的母离子，绘制所述指定子离子和与所述指定子离子对应的母离子的色谱图、小波图以及小波相关性图。
5.根据权利要求4所述的方法，其特征在于，所述方法还包括：
显示所述指定子离子色谱图和所述所有母离子色谱图，以供用户直观地看到所述指定子离子色谱图、所述所有母离子色谱图以及所述指定子离子色谱图的变化情况。
6.一种数据处理的装置，其特征在于，所述装置包括：
子母离子获取模块，用于获取指定子离子和所有母离子；
第一出峰时间获取模块，用于对所述指定子离子进行小波处理得到降噪后的指定子离子色谱图，并基于所述指定子离子色谱图确定所述指定子离子的出峰时间；
第二出峰时间获取模块，用于对所述所有母离子进行小波处理得到降噪后的所有母离子时间序列，并基于所述所有母离子时间序列确定所述所有母离子的出峰时间；
备选母离子获取模块，用于基于所述指定子离子的出峰时间和所述所有母离子的出峰时间获取出峰时间位于所述指定子离子的出峰时间预设范围内的多个备选母离子；
相关性数据生成模块，用于计算所述指定子离子和所述多个备选母离子的小波相关性并生成小波相关性数据；
对应母离子获取模块，用于基于所述小波相关性数据判断所述小波相关性数据高于预设数据的母离子，以获取与所述指定子离子对应的母离子。
7.根据权利要求6所述的装置，其特征在于，所述装置还包括：
子母离子采集模块，用于采集荷质比为1-1000的多个母离子和荷质比为1-1000且只带一个电荷的多个子离子；
子母离子存储模块，用于将所述多个子离子和所述多个母离子以预定格式进行存储。
8.根据权利要求7所述的装置，其特征在于，所述第一出峰时间获取模块包括：
出峰个数获取子模块，用于对所述指定子离子进行小波处理得到降噪后的指定子离子色谱图，并基于所述指定子离子色谱图确定所述指定子离子的出峰个数和每个峰出峰时间；
母离子个数确定子模块，用于根据所述出峰个数确定能生成所述指定子离子的母离子个数，其中，所述出峰个数与所述能生成所述指定子离子的母离子个数相等。
9.根据权利要求8所述的装置，其特征在于，所述装置还包括：
图形绘制模块，用于根据所述指定子离子和与所述指定子离子对应的母离子，绘制所述指定子离子和与所述指定子离子对应的母离子的色谱图、小波图以及小波相关性图。
10.根据权利要求9所述的装置，其特征在于，所述装置还包括：
显示模块，用于显示所述指定子离子色谱图和所述所有母离子色谱图，以供用户直观地看到所述指定子离子色谱图、所述所有母离子色谱图以及所述指定子离子色谱图的变化情况。
CN201810346349.9A 2018-04-17 2018-04-17 数据处理的方法及装置 Active CN108535403B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810346349.9A CN108535403B (zh) 2018-04-17 2018-04-17 数据处理的方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810346349.9A CN108535403B (zh) 2018-04-17 2018-04-17 数据处理的方法及装置

Publications (2)

Publication Number Publication Date
CN108535403A CN108535403A (zh) 2018-09-14
CN108535403B true CN108535403B (zh) 2019-07-02

Family

ID=63481328

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810346349.9A Active CN108535403B (zh) 2018-04-17 2018-04-17 数据处理的方法及装置

Country Status (1)

Country Link
CN (1) CN108535403B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113237985B (zh) * 2021-03-24 2023-10-20 迈克医疗电子有限公司 血红蛋白中组分类型的确定方法、装置和计算机设备
CN114137103B (zh) * 2021-11-03 2023-06-02 湖北省农业科学院农业质量标准与检测技术研究所 一种液相色谱二极管阵列数据转化为指纹图谱的方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104779132A (zh) * 2009-05-06 2015-07-15 Mks仪器公司 静电离子陷阱
CN105518455A (zh) * 2013-09-09 2016-04-20 株式会社岛津制作所 波峰检测方法
CN105718723A (zh) * 2016-01-18 2016-06-29 吉林大学 一种质谱数据处理中谱峰位置检测方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104779132A (zh) * 2009-05-06 2015-07-15 Mks仪器公司 静电离子陷阱
CN105518455A (zh) * 2013-09-09 2016-04-20 株式会社岛津制作所 波峰检测方法
CN105718723A (zh) * 2016-01-18 2016-06-29 吉林大学 一种质谱数据处理中谱峰位置检测方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
Arctic sea ice drift from wavelet analysis of N SCAT and special sensor microwave imager data;Antony K. Liu 等;《JOURNAL OF GEOPHYSICAL RESEARCH》;19990515;第104卷(第C5期);第11529-11538页
Comparison of public peak detection algorithms for MALDI mass spectrometry data analysis;Chao Yang 等;《BMC Bioinformatics》;20091231;第10卷(第4期);第1-13页
Improved peak detection in mass spectrum by incorporating continuous wavelet transform-based pattern matching;Pan Du 等;《BIOINFORMATICS》;20061231;第22卷(第17期);第2059-2065页
基于质谱的定量蛋白质组学策略和方法研究进展;常乘 等;《中国科学》;20151231;第45卷(第5期);第425-438页
液质联用技术用于复杂混合物体系中小分子化合物的分析;余玲玲 等;《中国科学》;20171231;第47卷(第12期);第1379-1391页
连续小波变换-纯变量法解析气相色谱-质谱重叠信号;*** 等;《郑州轻工业学院学报(自然科学版)》;20060228;第21卷(第1期);第12-15页

Also Published As

Publication number Publication date
CN108535403A (zh) 2018-09-14

Similar Documents

Publication Publication Date Title
JP6494588B2 (ja) 滞留時間の決定または確認のための窓処理質量分析データの使用
Bonner et al. SWATH data independent acquisition mass spectrometry for metabolomics
Castillo et al. Algorithms and tools for the preprocessing of LC–MS metabolomics data
Vaudel et al. Peptide and protein quantification: a map of the minefield
Kwiecien et al. High-resolution filtering for improved small molecule identification via GC/MS
CN105917221B (zh) 串联质谱分析数据处理装置
CN109643635A (zh) 用于在扫描swath数据中识别前体及产物离子对的***及方法
Bereman Tools for monitoring system suitability in LC MS/MS centric proteomic experiments
CN108535403B (zh) 数据处理的方法及装置
EP3572805A1 (en) Analysis data analytics method and analysis data analytics device
CN108982729A (zh) 用于提取质量迹线的***和方法
Lundgren et al. Protein identification using Sorcerer 2 and SEQUEST
CN107209151B (zh) 干扰检测及所关注峰值解卷积
EP3497709B1 (en) Automated spectral library retention time correction
Hoopmann et al. Identification of peptide features in precursor spectra using Hardklör and Krönik
CN105051526B (zh) 质量分析装置
Yates III Mass spectrometry as an emerging tool for systems biology
JP2018504607A (ja) 断片化エネルギーを切り替えながらの幅広い四重極ｒｆ窓の高速スキャニング
EP3335237B1 (en) Deconvolution of mixed spectra
CN109964300B (zh) 用于实时同位素识别的***和方法
Meng et al. LipidMiner: a software for automated identification and quantification of lipids from multiple liquid chromatography-mass spectrometry data files
JP6668188B2 (ja) 質量分析方法
Willems et al. Ion-networks: a sparse data format capturing full data integrity of data independent acquisition mass spectrometry
Fröhlich et al. Data-independent acquisition: A milestone and prospect in clinical mass spectrometry-based proteomics
Charles Development of a Software Package for the Quantitative Analysis of Proteomic Mass Spectrometry Datasets Labelled with Nitrogen-15

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