WO2015029254A1 - クロマトグラフ用データ処理装置及び方法 - Google Patents

クロマトグラフ用データ処理装置及び方法 Download PDF

Info

Publication number
WO2015029254A1
WO2015029254A1 PCT/JP2013/073555 JP2013073555W WO2015029254A1 WO 2015029254 A1 WO2015029254 A1 WO 2015029254A1 JP 2013073555 W JP2013073555 W JP 2013073555W WO 2015029254 A1 WO2015029254 A1 WO 2015029254A1
Authority
WO
WIPO (PCT)
Prior art keywords
wavelength
standard sample
intensity
chromatogram
designated
Prior art date
Application number
PCT/JP2013/073555
Other languages
English (en)
French (fr)
Inventor
悦輔 鎌田
年伸 柳沢
康敬 水戸
賢一 三嶋
Original Assignee
株式会社島津製作所
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 株式会社島津製作所 filed Critical 株式会社島津製作所
Priority to CN201380079317.2A priority Critical patent/CN105518454B/zh
Priority to PCT/JP2013/073555 priority patent/WO2015029254A1/ja
Priority to JP2015533921A priority patent/JP6065981B2/ja
Priority to US14/915,794 priority patent/US10151734B2/en
Priority to CN201810180594.7A priority patent/CN108508122B/zh
Publication of WO2015029254A1 publication Critical patent/WO2015029254A1/ja

Links

Images

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/8624—Detection of slopes or peaks; baseline correction
    • G01N30/8631—Peaks
    • G01N30/8634—Peak quality criteria
    • 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
    • G01N30/74—Optical detectors

Definitions

  • the present invention relates to a data processing apparatus and a data processing method for a chromatographic apparatus such as a liquid chromatograph and a gas chromatograph.
  • chromatogram data data representing a chromatogram
  • signal intensity output voltage, etc.
  • chromatograph data processing device a peak appearing in such a chromatogram is detected, a substance corresponding to the peak is identified from the peak position (retention time) with reference to a preset identification table, Further, the concentration and amount of the substance are calculated from the peak height and area.
  • the detector of the chromatographic apparatus has a problem that the reliability of the detection result differs depending on the level of the signal.
  • the component concentration in the sample is too low, the quantification accuracy deteriorates due to the effects of noise in the detection signal, and vice versa. If the component concentration is too high, the non-linearity of the detection signal becomes remarkable, and the quantitative accuracy is deteriorated. Therefore, analysis in a conventional chromatographic apparatus is performed by appropriately diluting the sample so that the component concentration in the sample falls within a predetermined range (dynamic range).
  • the sample concentration and detector sensitivity should be adjusted so that all target components are within the dynamic range. You only have to set it. However, if the concentration difference of the target component is large, the signal of the maximum concentration component (principal component) will be distorted or saturated if the minimum concentration component (impurity) is detected correctly, no matter what setting is made. Therefore, if the maximum concentration component (principal component) is to be detected correctly, the minimum concentration component (impurity) is buried in noise, and neither of them can be analyzed correctly.
  • Non-Patent Document 1 “Acetylcysteine Purity Test (6) Related Substances” (pp. 311-312), a liquid using an ultraviolet absorptiometer with a measurement wavelength of 220 nm It is described that when the test is performed with a chromatograph, the area of peaks other than acetylcysteine is 0.3% or less per peak, and 0.6% or less in total, compared with the area of the peak of acetylcysteine.
  • An object of the present invention is to provide a chromatographic data processing apparatus and data capable of performing data processing in one analysis and using one detector without being affected by noise and non-linearity of a detection signal. It is to provide a processing method.
  • the chromatographic data processing apparatus which has been made to solve the above-mentioned problems, creates a chromatogram based on a temporal change of a spectrum acquired within a wavelength range including a target wavelength corresponding to a target component.
  • a standard sample sensitivity coefficient calculating means for calculating a standard sample sensitivity coefficient divided by a second intensity that is an intensity at a specified wavelength; a-3) corrected standard sample chromatogram intensity calculating means for calculating the peak intensity of the corrected standard sample chromatogram obtained by multiplying the chromatogram at the specified wavelength by the standard sample sensitivity coefficient for each standard sample; a-4) index value calculation means for obtaining an index value indicating the correlation between the concentration of the standard sample and the peak intensity of the corrected standard sample chromatogram; a-5) One or both of the designated retention time and the designated wavelength are changed, and the set of the designated retention time and the designation is within a predetermined range indicating that the index value is highly correlated.
  • a specific designated holding time / specific designated wavelength setting means for selecting a wavelength and setting a specific designated holding time and a specific designated wavelength; and b-1) measurement sample data storage means for storing spectrum data for a measurement sample whose concentration of the target component is unknown; b-2) Calculate the measurement sample sensitivity coefficient obtained by dividing the first intensity, which is the intensity at the target wavelength, of the spectrum of the measurement sample at the specific designated retention time by the second intensity, which is the intensity at the specific designated wavelength.
  • a measurement sample sensitivity coefficient calculating means b-3)
  • a corrected measurement sample chromatogram creation means for creating a corrected measurement sample chromatogram obtained by multiplying the chromatogram at the specific designated wavelength by the measurement sample sensitivity coefficient.
  • the chromatographic data processing apparatus prepares a chromatogram using a plurality of standard samples (concentrations known) having different concentrations of the target component described in a-1) to a-5).
  • a configuration for specifying parameters to be used, and a configuration for creating a chromatogram of a measurement sample whose concentration of a target component is unknown using the parameters described in b-1) to b-3) Have
  • the configurations of a-1) to a-5) are collectively referred to as “parameter setting means”
  • the configurations of b-1) to b-3) are collectively referred to as “measurement sample chromatogram creation means”.
  • the measurement sample chromatogram creation means will be described first. Because the spectrum of a component has a shape unique to the component, the spectrum of the component has the same shape and is different only in intensity if no distortion or saturation occurs in the retention time within the peak of the chromatogram by that component. Have similarities. When the chromatogram is created at each wavelength belonging to one spectrum peak due to the similarity of the spectrum shape, the intensity of the chromatogram is proportional to the intensity of the spectrum at the wavelength. Therefore, in the measurement sample chromatogram creation means, for the spectrum of the measurement sample in a predetermined retention time (this is called the designated retention time), the intensity of the spectrum at the target wavelength is different from that in the same spectrum peak.
  • An operation of multiplying the chromatogram intensity at the designated wavelength by the sensitivity coefficient divided by the intensity of the spectrum in this (referred to as the designated wavelength), that is, the intensity ratio of the spectrum at the target wavelength and the designated wavelength, is performed.
  • the intensity ratio of such a spectrum matches the intensity ratio of the chromatogram at the target wavelength and the chromatogram at the designated wavelength due to the above similarity unless the spectrum is distorted or saturated. Therefore, even if the spectrum is distorted or saturated at the retention time near the peak top of the chromatogram, the target wavelength can be determined by specifying the specified wavelength and the specified retention time at which such saturation does not occur.
  • a chromatogram corrected sample chromatogram after correction
  • the parameter setting means performs the setting.
  • the parameter setting means calculates the sensitivity coefficient (standard sample sensitivity coefficient) for each standard sample after setting the specified retention time and the specified wavelength (temporarily) for a plurality of standard samples with different concentrations of the target component. ) Calculate the peak intensity of the corrected standard sample chromatogram obtained by multiplying the chromatogram at the specified wavelength by the standard sample sensitivity coefficient.
  • an index value (described later) indicating the correlation between the concentration of the standard sample and the peak intensity of the corrected standard sample chromatogram is obtained while changing either one or both of the designated holding time and the designated wavelength.
  • a set of designated holding time and designated wavelength when this index value is within a predetermined range is selected and specified as the specified designated holding time and specified designated wavelength.
  • obtaining the relationship between the concentration of the standard sample and the peak intensity of the chromatogram (normal, not the standard sample chromatogram after correction) in the parameter setting means corresponds to creating a calibration curve.
  • the specific designated retention time and the specific designated wavelength are specified based on the data of the standard sample measured in such a wide concentration range.
  • the index value includes (1) absolute value of correlation coefficient, (2) average value of absolute value of deviation in peak intensity of corrected standard sample chromatogram, and (3) maximum value of absolute value of deviation, etc. These can be used, and two or more of these can be used in combination. Hereinafter, these index values will be described.
  • the correlation coefficient C is the concentration d i in the standard sample of n types (n is a natural number of 2 or more) having different concentrations of the target component and the peak intensity I of the corrected standard sample chromatogram.
  • n a natural number of 2 or more
  • I the peak intensity of the corrected standard sample chromatogram.
  • d av and I av are average values of n pieces of data ⁇ d i ⁇ and ⁇ I i ⁇ , respectively.
  • the correlation coefficient C can take a value from ⁇ 1 to +1.
  • the correlation coefficient C means that the correlation between two variables increases as the correlation coefficient C approaches +1 or ⁇ 1 (that is, the absolute value approaches 1).
  • a correlation coefficient of -1 means that when one variable increases, the other variable decreases, so when the concentration and peak intensity are variables as in the present invention, the correlation The number must approach +1. Therefore, in the present invention, the specified holding time and the specified wavelength when the correlation coefficient C is closer to +1 than the predetermined value, that is, the predetermined value or more are set as the specific specified holding time and the specific specified wavelength.
  • the chromatographic data processing apparatus has been described as having both the parameter setting means and the measurement sample chromatogram creation means, but the chromatographic data processing apparatus has only the parameter setting means (ie, A calibration curve may only be created), and a chromatogram of the measurement sample may be separately created based on the parameters set by the parameter setting means.
  • a chromatographic data processing method is a chromatographic data processing method for creating a chromatogram based on a time change of a spectrum acquired within a wavelength range including a target wavelength corresponding to a target component.
  • a first intensity which is an intensity at the target wavelength, of the spectrum of the standard sample at a designated retention time common to each standard sample, Calculating a standard sample sensitivity coefficient divided by a second intensity which is an intensity at a specified wavelength within the wavelength range different from the wavelength;
  • the first intensity, which is the intensity at the target wavelength, of the spectrum of the measurement sample in the specific designated retention time is the second intensity, which is the intensity at the specific designated wavelength.
  • analysis data processing using a chromatograph can be performed in one analysis and using one detector without being affected by noise and non-linearity of the detection signal. . And since the designated holding time and the designated wavelength, which are parameters necessary for the data processing, can be automatically specified, even a beginner can easily handle them.
  • FIG. 1 is a schematic configuration diagram of an analysis system including an embodiment of a chromatographic data processing apparatus according to the present invention.
  • the graph which shows the example of the calibration curve (a) when not correcting chromatogram, and the calibration curve (b) calculated
  • the schematic block diagram which shows the other example of the data processor for chromatographs which concerns on this invention.
  • Embodiments of the chromatograph data processing apparatus according to the present invention will be described with reference to FIGS.
  • the case of a liquid chromatograph (LC) will be described as an example, but the same applies to a gas chromatograph.
  • the chromatographic data processing apparatus of the present embodiment constitutes a part of the analysis system 10 shown in FIG.
  • the analysis system 10 processes LC 11 for separating components contained in a liquid sample in time, a detector 12 for analyzing each separated component within a predetermined wavelength band, and data output from the detector 12.
  • a chromatograph data processing device device of this embodiment; hereinafter referred to as “data processing device”.
  • the data processing device 13 includes a general computer (hardware) having a CPU (central processing unit), a storage device (memory, hard disk, solid state drive, etc.), a display, an input device (keyboard, mouse, etc.), etc. Consists of dedicated data processing software installed on the computer.
  • the data processing device 13 has functions as a parameter setting unit 131 and a measurement sample chromatogram creation unit 132.
  • the parameter setting unit 131 is used to set the specific designated retention time and the specific designated wavelength using measurement results of a plurality of standard samples whose concentrations of the target component are known and different from each other.
  • Parameter search condition input unit 1310 for a user to input various conditions to be described later, and a standard sample data storage unit for storing spectrum measurement data for a plurality of holding times for a plurality of types of standard samples having different concentrations 1311.
  • the parameter setting unit 131 includes a standard sample sensitivity coefficient calculation unit 1312, a corrected standard sample chromatogram intensity calculation unit 1313, an index value calculation unit 1314, a specific designated retention time / specific designated wavelength setting unit 1315, and a calibration curve creation unit. 1316. These units will be described when the operation of the data processing device 13 is described.
  • the measurement sample chromatogram creation unit 132 performs a process using the measurement result of the measurement sample in which the concentration of the target component is unknown, the specific designated retention time and the specific designated wavelength, and finally corrects the measurement sample.
  • a measurement sample data storage unit 1321, a measurement sample sensitivity coefficient calculation unit 1322, and a corrected measurement sample chromatogram creation unit 1323 for creating a chromatogram are provided.
  • the measurement sample data storage unit 1321 stores spectrum measurement data for a number of holding times for each measurement sample. The functions of the measurement sample sensitivity coefficient calculation unit 1322 and the corrected measurement sample chromatogram creation unit 1323 will be described when the operation of the data processing device 13 is described.
  • the display of the data processing device 13 displays two windows, a quantitative browser window 20 shown in FIG. 2 and a data analysis window 40 shown in FIG.
  • a quantitative browser window 20 shown in FIG. 2 displays two windows
  • the user clicks a predetermined icon in the quantitative browser window 20 or performs a predetermined operation such as selecting a predetermined item from “Window (W)” in the pull-down menu the parameters shown in FIG.
  • a search condition input window 30 is displayed.
  • the quantitative browser window 20, the parameter search condition input window 30, and the data analysis window 40 will be described in this order.
  • the quantification browser window 20 is a window used when setting parameters using a standard sample, and has display areas such as a quantification result view 21, a chromatogram view 22, and a calibration curve / spectrum view 23.
  • the quantification result view 21 shows the area and height obtained from the name of the data file and the chromatogram for each standard sample for a plurality of standard samples whose concentration of the target component is known (the numerical values of these areas and heights will be described later).
  • the information can be changed by data processing by the data processing device 13), and information such as concentration is displayed.
  • the chromatogram view 22 displays a chromatogram for one type of standard sample selected by the user using a mouse device or the like in the quantitative result view 21.
  • the chromatogram view 22 can also display detailed information of the sample recorded in the data file regarding the one type of standard sample by switching the tabs displayed therein.
  • the calibration curve / spectrum view 23 displays either a calibration curve or a spectrum at one standard sample and a retention time by switching tabs.
  • the calibration curve is a graph showing the relationship between the concentration of the standard sample and the intensity (area or height) of the chromatogram, in order to obtain the concentration of the measurement sample from the intensity of the chromatogram of the measurement sample whose concentration is unknown. Used for. This calibration curve can be changed by data processing by the data processing device 13 described later.
  • the parameter search condition input window 30 is a window for executing automatic parameter search after a user inputs parameter setting conditions using a standard sample.
  • correction wavelength setting method correction wavelength (manual), correction wavelength (automatic) intensity, correction wavelength (automatic) moving direction, sensitivity correction spectrum extraction intensity, and background correction are shown.
  • An input field is displayed.
  • a button for the user to select whether the above-mentioned specific designated wavelength is set automatically or manually is displayed.
  • the correction wavelength (manual) input column allows the user to input a numerical value of a specific designated wavelength when “manual” is selected in the correction wavelength setting method input column.
  • a column for selecting a calculation method for obtaining an index value indicating the correlation between the concentration of the standard sample and the intensity of the standard sample chromatogram after correction is displayed.
  • one or a plurality of index values can be selected from three types of index values: a correlation coefficient of a calibration curve, an average value of deviation, and a maximum value of deviation.
  • the parameter search condition input window 30 displays a search start button for starting an automatic parameter search after the input of these conditions is completed.
  • the data analysis window 40 is a window used when processing the data of the measurement sample, and has display areas such as a contour line view 41, a spectrum view 42, and a chromatogram view 43.
  • the contour line view 41 shows the signal intensity of the detector 12 as a contour line on a graph in which the horizontal axis represents the holding time and the vertical axis represents the wavelength of the light detected by the detector 12.
  • the spectrum view 42 displays the spectrum of the measurement sample at a specific holding time
  • the chromatogram view 43 displays the corrected chromatogram of the measurement sample.
  • one retention time T s in which saturation does not occur in the pre-correction chromatogram 511 is designated (FIG. 5 (a)).
  • This holding time T s is referred to as “designated holding time”.
  • the designated holding time T s may be designated directly, but the saturation of the pre-correction chromatogram 511 does not occur when the intensity is equal to or lower than a predetermined value.
  • the designated holding time T s corresponding to the intensity may be designated.
  • the value input in the sensitivity correction spectrum extraction intensity input field of the parameter search condition input window 30 described above corresponds to the intensity Is1 specified here.
  • the spectrum (referred to as “second spectrum”) 522 in the designated retention time T s does not show saturation or distortion (FIG. 5B).
  • a wavelength ⁇ s at which no saturation or distortion is observed in the first spectrum 521 is specified (FIG. 5B).
  • This wavelength ⁇ s is referred to as a "specified wavelength”.
  • a designated wavelength lambda s instead of specifying the specified wavelength lambda s directly, by selecting the intensity I s2 below a predetermined value saturation spectrum does not occur, may be designated a designated wavelength lambda s which corresponds to the intensity.
  • the value input in the correction wavelength (automatic) intensity input field of the parameter search condition input window 30 described above corresponds to the intensity Is2 specified here.
  • the first intensity I a that is the intensity of the second spectrum 522 at ⁇ t where the target component has the maximum spectrum (in general, if the target component is known, this wavelength is also known)
  • a second intensity I b that is the intensity of the second spectrum 522 at the specified wavelength ⁇ s is obtained (FIG. 5C).
  • a sensitivity coefficient K I a / I b which is a ratio between the first intensity I a and the second intensity I b is obtained.
  • the ratio of the first intensity I a and the second intensity in the second spectrum 522 is equal to the ratio of the intensity at the target wavelength ⁇ t of the first spectrum 521 and the intensity at the designated wavelength ⁇ s . Therefore, by creating a chromatogram 531 at the designated wavelength ⁇ s where no saturation / distortion occurs in the first spectrum 521, the saturation / distortion at the target wavelength ⁇ t is obtained by multiplying the chromatogram 531 by the sensitivity coefficient K. No corrected chromatogram 532 is obtained.
  • the user selects “automatic” in the correction wavelength setting method input field, and the sensitivity correction spectrum described above in the sensitivity correction spectrum extraction intensity input field and correction wavelength (automatic) intensity input field, respectively.
  • these values are input as initial values, and finally the optimum values are obtained by the data processing device 13 of this embodiment.
  • the user corrects the sensitivity correction spectrum extraction intensity Is1 and the correction wavelength (Automatic) Enter the numerical value of the search range and search step (interval) for changing the intensity Is2 .
  • the user inputs whether or not to perform background correction in the background correction input field. Since background correction is generally performed during chromatogram data processing, detailed description thereof is omitted here.
  • the user further selects a calculation method for obtaining an index value indicating the correlation between the standard sample concentration and the intensity of the corrected standard sample chromatogram in the lower column of the parameter search condition input window 30.
  • the wavelength corresponding to the compensation wavelength (auto) intensity I s2 is present on both the long wavelength side and the shorter wavelength side than the target wavelength lambda t.
  • the user inputs which wavelength to select in the correction wavelength (automatic) movement direction input column.
  • one specified holding time T s may be input instead of the sensitivity correction spectrum extraction intensity I s1
  • one specified wavelength ⁇ s may be input instead of the corrected wavelength (automatic) intensity I s2. It may be.
  • step S ⁇ b> 3 the standard sample sensitivity coefficient calculation unit 1312 acquires spectrum data for each standard sample for a plurality of types of standard samples from the standard sample data storage unit 1311.
  • the standard sample sensitivity coefficient calculation unit 1312 calculates the sensitivity coefficient K described in (2) for each standard sample using one set of designated holding time T s and designated wavelength ⁇ s (step S4).
  • the corrected standard sample chromatogram intensity calculation unit 1313 creates a corrected chromatogram 532 as described in (2) based on the sensitivity coefficient K obtained in step S4 for each standard sample.
  • the peak intensity (peak area intensity or peak top intensity) of the corrected chromatogram 532 is calculated (step S5).
  • the index value calculation unit 1314 calculates index values for the set of designated holding times T s and designated wavelengths ⁇ s (step S6).
  • the index value is calculated from the correlation coefficient of the calibration curve, the average value of the deviation, and the median value of the deviation, which is selected by the user in the parameter search condition input window 30.
  • the correlation coefficient of the calibration curve is calculated using the above-described calculation formula of the correlation coefficient C, and the average value and median value of the deviation are obtained by calculating the average value and median value by a known method. calculate.
  • the index value calculation unit 1314 corresponds to the search range and search step entered by the user in the parameter search condition input window 30. It is checked whether or not the calculation of the index value is completed for all combinations of the designated holding time T s and the designated wavelength ⁇ s to be performed (step S7). If it is completed, the process proceeds to step S9. On the other hand, if it has not been completed, the designated holding time T s and the designated wavelength ⁇ s are changed to those for which calculation of the index value has not been completed (step S8), and the operations of steps S4 to S6 are repeated.
  • the specific designated retention time / specific designated wavelength setting unit 1315 indicates that the index value has the highest correlation among the combinations of the plural designated retention times T s and the designated wavelengths ⁇ s obtained so far.
  • the designated holding time T s and the designated wavelength ⁇ s indicating are selected as the specific designated holding time and the specific designated wavelength.
  • the intensity value corresponding to the specific designated holding time and the specific designated wavelength, and the calculation result of the correlation coefficient are displayed in the corresponding column in the search condition input window 30.
  • calibration curve creating unit 1316 determines the particular specified retention time and data row obtained in certain specified wavelengths ⁇ (d i, I i) ⁇ to prepare a calibration curve based on.
  • the data point deviates from a linear straight line due to saturation of the chromatogram peak (FIG. 7 (a)), whereas the calibration curve obtained here shows the data The points agree well with the straight line of the linear function (FIG. 7 (b)).
  • the calibration curve creation unit 1316 displays the corrected calibration curve thus obtained on the calibration curve / spectrum view 23 of the quantitative browser window 20. Thereby, the operation of the parameter setting unit 131 ends.
  • (3-2) Operation of the measurement sample chromatogram creation unit 132 When the user performs a predetermined operation for specifying one of the measurement samples whose concentration of the target component is unknown in the data analysis window 40, the measurement sample chromatogram creation is performed.
  • the unit 132 starts operation. First, the measurement sample sensitivity coefficient calculation unit 1322 acquires spectrum data of the measurement sample from the measurement sample data storage unit 1321 (step S21). Next, the measurement sample sensitivity coefficient calculation unit 1322 calculates the sensitivity coefficient K described in (2), where T s is the specific designated retention time selected by the parameter setting unit 131 and ⁇ s is the specific designated wavelength ( Step S22).
  • the corrected measurement sample chromatogram creation unit 1323 creates the corrected chromatogram 532 for the measurement sample as described in (2) based on the sensitivity coefficient K obtained in step S23 (step 2). S23). Thus, the operation of the measurement sample chromatogram creation unit 132 is finished.
  • the calibration curve creation unit 1316 creates a calibration curve based on the data obtained at the specific designated retention time and the specific designated wavelength.
  • the calibration curve creation unit 1316 is not essential.
  • it when it is intended to create a calibration curve that is not affected by chromatogram saturation or distortion, it has a parameter setting unit 131 as in the analysis system 10A shown in FIG. A configuration without the chromatogram creation unit 132 may also be adopted.

Landscapes

  • Physics & Mathematics (AREA)
  • Biochemistry (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Engineering & Computer Science (AREA)
  • Quality & Reliability (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Treatment Of Liquids With Adsorbents In General (AREA)

Abstract

　クロマトグラフ用データ処理装置は、目的成分の濃度が異なる複数の標準試料について、標準試料毎にスペクトルデータを記憶する標準試料データ記憶部１３１１と、各標準試料に共通の指定保持時間における該標準試料のスペクトルの、目的波長における強度である第１強度を、該目的波長とは異なる所定波長範囲内の指定波長における強度である第２強度で除した標準試料感度係数を算出する標準試料感度係数算出部１３１２と、指定波長におけるクロマトグラムに標準試料感度係数を乗算した補正後標準試料クロマトグラムのピーク強度を算出する補正後標準試料クロマトグラム強度算出部１３１３と、標準試料の濃度と補正後標準試料クロマトグラムのピーク強度の相関性を示す指標値を求める指標値算出部１３１４と、指定保持時間及び指定波長のいずれか一方又は両方を変化させ、指標値が、相関性が高いことを示す所定の範囲内にある１組の該指定保持時間及び該指定波長を選択して特定指定保持時間及び特定指定波長を設定する特定指定保持時間・特定指定波長設定部１３１５と、目的成分の濃度が未知である測定試料についてスペクトルのデータを記憶する測定試料データ記憶部１３２１と、特定指定保持時間における測定試料感度係数を算出する測定試料感度係数算出部１３２２と、特定指定波長におけるクロマトグラムに測定試料感度係数を乗算した補正後測定試料クロマトグラムを作成する補正後測定試料クロマトグラム作成部１３２３を備える。

Description

クロマトグラフ用データ処理装置及び方法
　本発明は、液体クロマトグラフ、ガスクロマトグラフなどのクロマトグラフ装置用のデータ処理装置及びデータ処理方法に関する。
　クロマトグラフ装置では、試料を分析することによって横軸に時間、縦軸に信号強度（出力電圧など）をとったクロマトグラムを表すデータ（以下、クロマトグラムデータという）を取得する。クロマトグラフ用データ処理装置では、このようなクロマトグラムに出現しているピークを検出し、予め設定された同定テーブルを参照してピーク位置（保持時間）からそのピークに対応する物質を同定し、さらにピークの高さや面積からその物質の濃度や量を算出する。
　こうしたデータ処理装置では、一般に、A/D変換器を含む信号処理回路のハードウエア上での制約から処理可能な信号の大きさに限界があり、上限以上又は下限以下の大きさの信号が入力されても正確な演算処理を実行することができない。
　また、こうした信号処理上での限界とは別に、クロマトグラフ装置の検出器には、その信号のレベルによって検出結果の信頼度が相違するという問題がある。例えば液体クロマトグラフの検出器として使用される紫外可視分光光度計、フォトダイオードアレイ検出器などでは、試料中の成分濃度が低すぎると検出信号におけるノイズの影響により、定量精度が悪化し、その反対に成分濃度が高すぎると検出信号の非直線性が顕著になって、やはり定量精度が悪化してしまう。そのため、従来のクロマトグラフ装置における分析は、試料中の成分濃度が所定の範囲（ダイナミックレンジ）内となるように、試料を適宜希釈して行っていた。
第十六改正日本薬局方,[online],2011年3月24日,厚生労働省,[2013年6月21日検索],インターネット<URL: http://jpdb.nihs.go.jp/jp16/>
　試料中に複数の分析目的成分が存在する場合において、それら複数の分析目的成分の濃度が大きく異ならない場合には、全目的成分がダイナミックレンジに収まるように試料の濃度や検出器の感度などを設定すればよい。しかし、目的成分の濃度差が大きい場合、それらについてどのような設定を行っても、最小濃度の成分（不純物）を正しく検出しようとすると最大濃度の成分（主成分）の信号が歪み又は飽和してしまい、最大濃度の成分（主成分）を正しく検出しようとすると最小濃度の成分（不純物）がノイズに埋もれてしまい、いずれも正しい分析を行うことができない。
　例えば、第十六改正日本薬局方（非特許文献１）の「アセチルシステイン　純度試験（6）類縁物質」（pp. 311-312）では、測定波長を220nmとした紫外吸光光度計を用いた液体クロマトグラフで試験を行ったとき、アセチルシステインのピークの面積と比較して、アセチルシステイン以外のピークの面積がピーク1つあたり0.3%以下、合計で0.6%以下となることが記載されている。このような大きな濃度差がある場合、濃度の低いアセチルシステイン以外の成分のピークの面積を正確に求めるために試料全体の濃度を高く（希釈率を小さく）すると、アセチルシステインの信号の歪み又は飽和が生じてしまう。
　そのため、従来は、(i)濃度が異なる複数の試料を用意してそれぞれ分析を行い、得られた結果に希釈率の補正を行うことで各目的成分の濃度（又はそれらの濃度比）を求めたり、あるいは、(ii)光路長の異なるセルを２台の検出器に設置し、光路長補正を行うことで、１回の分析で各目的成分の濃度（又はそれらの濃度比）を求めていた。しかし、(i)の方法は測定に要する時間が長くなるため非効率であり、(ii)の方法は検出器を増設する必要がありコストを要するという欠点がある。
　本発明の課題は、1回の分析で且つ1台の検出器を用いて、検出信号のノイズ及び非直線性の影響を受けることなくデータ処理を行うことができるクロマトグラフ用データ処理装置及びデータ処理方法を提供することである。
　上記課題を解決するために成された本発明に係るクロマトグラフ用データ処理装置は、目的成分に対応した目的波長を含む波長範囲内で取得したスペクトルの時間変化に基づいてクロマトグラムを作成するためのクロマトグラフ用データ処理装置であって、
　a-1) 前記目的成分の濃度が異なる複数の標準試料について、標準試料毎に前記スペクトルのデータを記憶する標準試料データ記憶手段と、
　a-2) 前記標準試料毎に、各標準試料に共通の指定保持時間における該標準試料のスペクトルの、前記目的波長における強度である第１強度を、該目的波長とは異なる前記波長範囲内の指定波長における強度である第２強度で除した標準試料感度係数を算出する標準試料感度係数算出手段と、
　a-3) 前記標準試料毎に、前記指定波長におけるクロマトグラムに前記標準試料感度係数を乗算した補正後標準試料クロマトグラムのピーク強度を算出する補正後標準試料クロマトグラム強度算出手段と、
　a-4) 前記標準試料の濃度と前記補正後標準試料クロマトグラムのピーク強度の相関性を示す指標値を求める指標値算出手段と、
　a-5) 前記指定保持時間及び前記指定波長のいずれか一方又は両方を変化させ、前記指標値が、相関性が高いことを示す所定の範囲内にある１組の該指定保持時間及び該指定波長を選択して特定指定保持時間及び特定指定波長を設定する特定指定保持時間・特定指定波長設定手段と、
　b-1) 前記目的成分の濃度が未知である測定試料についてスペクトルのデータを記憶する測定試料データ記憶手段と、
　b-2) 前記特定指定保持時間における前記測定試料のスペクトルの、前記目的波長における強度である第１強度を、前記特定指定波長における強度である第２強度で除した測定試料感度係数を算出する測定試料感度係数算出手段と、
　b-3) 前記特定指定波長におけるクロマトグラムに前記測定試料感度係数を乗算した補正後測定試料クロマトグラムを作成する補正後測定試料クロマトグラム作成手段と
を備えることを特徴とする。
　本発明に係るクロマトグラフ用データ処理装置は、a-1)～a-5)に記載の、目的成分の濃度が異なる複数の（濃度が既知である）標準試料を用いて、クロマトグラムを作成するのに用いるパラメータを特定するための構成と、b-1)～b-3)に記載の、前記パラメータを用いて目的成分の濃度が未知である測定試料のクロマトグラムを作成するための構成を有する。以下では、a-1)～a-5)の構成を合わせて「パラメータ設定手段」、b-1)～b-3)の構成を合わせて「測定試料クロマトグラム作成手段」と呼ぶ。
　説明の都合上、まず、測定試料クロマトグラム作成手段について説明する。
　ある成分のスペクトルは成分固有の形状を有するため、その成分によるクロマトグラムのピーク内の保持時間では、その成分のスペクトルは、歪みや飽和が生じなければ、形状は同じであって強度のみが相違するという相似性を有する。このスペクトル形状の相似性により、１つのスペクトルピークに属する各波長においてクロマトグラムを作成すると、クロマトグラムの強度は、当該波長におけるスペクトルの強度に比例することとなる。そこで、測定試料クロマトグラム作成手段では、所定の保持時間（これを指定保持時間と呼ぶ。）における測定試料のスペクトルにつき、目的波長における当該スペクトルの強度を、同じスペクトルピーク内のそれとは別の波長（これを指定波長と呼ぶ。）における当該スペクトルの強度で除した感度係数、すなわち目的波長と指定波長におけるスペクトルの強度比を、該指定波長におけるクロマトグラム強度に乗算する操作を行う。このようなスペクトルの強度比は、当該スペクトルに歪みや飽和が生じていない限り、上記の相似性により、目的波長におけるクロマトグラムと指定波長におけるクロマトグラムの強度比と一致する。そのため、たとえ、クロマトグラムのピークトップ付近の保持時間においてスペクトルに歪み又は飽和が生じた場合であっても、そのような飽和が生じていない指定波長及び指定保持時間を特定することにより、目的波長におけるクロマトグラムと同じ強度で、歪みや飽和を補正したクロマトグラム（補正後測定試料クロマトグラム）を得ることができる。
　従って、歪みや飽和のないクロマトグラムを作成するためには、適切な指定保持時間及び指定波長（これらを「特定指定保持時間」及び「特定指定波長」と呼ぶ。）を設定することが必要となる。その設定を行うのがパラメータ設定手段である。以下、パラメータ設定手段について説明する。
　パラメータ設定手段では、目的成分の濃度が異なる複数の標準試料について、指定保持時間及び指定波長を（仮に）定めたうえで標準試料毎に感度係数（標準試料感度係数）を算出し、（仮の）指定波長におけるクロマトグラムに標準試料感度係数を乗算した補正後標準試料クロマトグラムのピーク強度を算出する。そして、指定保持時間及び指定波長のいずれか一方又は両方を変化させながら、標準試料の濃度と補正後標準試料クロマトグラムのピーク強度の相関性を示す指標値（後述）を求める。そして、この指標値が所定の範囲内にあるときの１組の指定保持時間及び指定波長を選択して特定指定保持時間及び特定指定波長として特定する。ここで、指標値が所定の範囲内にある指定保持時間及び指定波長の組は多数存在し得る。それらの中で、指標値が最適なものを選択することが望ましいが、最適なものを選択するための演算を省略して任意の１組を選択するようにしてもよい。
　上記のように、パラメータ設定手段において標準試料の濃度と（補正後標準試料クロマトグラムではなく、通常の）クロマトグラムのピーク強度の関係を求めることは、検量線を作成することに相当する。検量線を作成する際には、未知試料の濃度の特定に用いることができるように、当該未知試料で想定される濃度よりも広い濃度範囲で標準試料の測定を行うのが一般的である。本発明においても、このような広い濃度範囲で測定した標準試料のデータに基づいて特定指定保持時間及び特定指定波長を特定する。
　前記指標値には、(1)相関係数の絶対値、(2)補正後標準試料クロマトグラムのピーク強度における偏差の絶対値の平均値、(3)前記偏差の絶対値の最大値などを用いることができ、これらのうちの2つ以上を組み合わせて用いることもできる。以下、これらの指標値について説明する。
(1)相関係数
　本発明では、相関係数Cは、目的成分の濃度が異なるn種（nは2以上の自然数）の標準試料における濃度diと補正後標準試料クロマトグラムのピーク強度Iiのデータ列{(di, Ii)}(i=1, 2, …n)について、
Figure JPOXMLDOC01-appb-M000002
 で与えられる。ここで、dav, Iavはそれぞれ、n個のデータ{di}, {Ii}の平均値である。相関係数Cは、-1から+1までの値を取り得る。一般的に、相関係数Cは+1又は-1に（即ち、絶対値が1に）近づくほど2つの変数の相関性が高いことを意味する。しかし、相関係数が-1であることは、一方の変数が増加すると他方の変数が減少することを意味するため、本発明のように濃度とピーク強度が変数である場合には、相関係数は+1に近づかなければならない。そのため、本発明では、相関係数Cが所定値よりも+1に近い、すなわち当該所定値以上のときの指定保持時間及び指定波長を、前記特定指定保持時間及び前記特定指定波長として設定する。
(2)偏差の絶対値の平均値
　上述のデータ列{(di, Ii)}(i=1, 2, …n)について、最小二乗法等を用いて濃度dとピーク強度Iの関係を１次関数I=kd（kは係数）で求め、各濃度diにおいて、その１次関数の値kdiとピーク強度Iの実測値Iiとの差である偏差を求める。そして、この偏差の絶対値の平均値が所定値以下であるときの指定保持時間及び指定波長を、前記特定指定保持時間及び前記特定指定波長としてとして設定する。
(3)偏差の絶対値の最大値
　上記(2)と同様に求められる、各濃度diにおける偏差の絶対値につき、その最大値を前記指標値として用いる。この最大値が所定値以下であるときの指定保持時間及び指定波長を、前記特定指定保持時間及び前記特定指定波長として設定する。
　ここまで、本発明に係るクロマトグラフ用データ処理装置はパラメータ設定手段と測定試料クロマトグラム作成手段の双方を有するものとして説明したが、クロマトグラフ用データ処理装置はパラメータ設定手段のみを有し（すなわち検量線の作成のみを行い）、測定試料のクロマトグラムの作成は当該パラメータ設定手段で設定されたパラメータに基づいて別途行うようにしてもよい。
　本発明に係るクロマトグラフ用データ処理方法は、目的成分に対応した目的波長を含む波長範囲内で取得したスペクトルの時間変化に基づいてクロマトグラムを作成するためのクロマトグラフ用データ処理方法であって、
　前記目的成分の濃度が異なる複数の標準試料について、該標準試料毎に、各標準試料に共通の指定保持時間における該標準試料のスペクトルの、前記目的波長における強度である第１強度を、該目的波長とは異なる前記波長範囲内の指定波長における強度である第２強度で除した標準試料感度係数を算出し、
　前記標準試料毎に、前記指定波長におけるクロマトグラムに前記標準試料感度係数を乗算した補正後標準試料クロマトグラムのピーク強度を算出し、
　前記標準試料の濃度と前記補正後標準試料クロマトグラムのピーク強度の相関性を示す指標値を求め、
　前記指定保持時間及び前記指定波長のいずれか一方又は両方を変化させながら、前記指標値が、相関性が高いことを示す所定の範囲内にある１組の該指定保持時間及び該指定波長を選択することにより特定指定保持時間及び特定指定波長を設定し、
　前記目的成分の濃度が未知である測定試料について、前記特定指定保持時間における該測定試料のスペクトルの、前記目的波長における強度である第１強度を、前記特定指定波長における強度である第２強度で除した測定試料感度係数を算出し、
　前記特定指定波長におけるクロマトグラムに前記測定試料感度係数を乗算した補正後測定試料クロマトグラムを作成する
ことを特徴とする。
　本発明によれば、検出信号のノイズ及び非直線性の影響を受けることなく、1回の分析で且つ1台の検出器を用いて、クロマトグラフを用いた分析のデータ処理を行うことができる。そして、このデータ処理の際に必要となるパラメータである指定保持時間及び指定波長を自動的に特定することができるため、初心者でも容易に取り扱うことができる。
本発明に係るクロマトグラフ用データ処理装置の一実施例を含む分析システムの概略構成図。 パラメータ設定に用いる定量ブラウザウインドウの一例を示す図。 パラメータ設定に用いるパラメータ探索条件入力ウインドウの一例を示す図。 測定試料のクロマトグラムの作成に用いるデータ解析ウインドウの一例を示す図。 本実施例のクロマトグラフ用データ処理装置装置で用いる、クロマトグラフ用データにおける検出信号のノイズの除去方法を説明するための図。 本実施例のクロマトグラフ用データ処理装置におけるパラメータ設定部の動作を説明するためのフローチャート。 クロマトグラムの補正を行わない場合の検量線(a)、及び本実施例の方法によりクロマトグラムの補正を行ったうえで求めた検量線(b)の例を示すグラフ。 本実施例のクロマトグラフ用データ処理装置における測定試料クロマトグラム作成部の動作を説明するためのフローチャート。 本発明に係るクロマトグラフ用データ処理装置の他の例を示す概略構成図。
　本発明に係るクロマトグラフ用データ処理装置の実施例を、図１～図９を用いて説明する。なお、以下では液体クロマトグラフ（ＬＣ）の場合を例に説明するが、ガスクロマトグラフの場合でも同様である。
(1) 本実施例のクロマトグラフ用データ処理装置の構成
　本実施例のクロマトグラフ用データ処理装置は、図１に示した分析システム１０の一部を構成する。分析システム１０は、液体試料中の含有成分を時間的に分離するＬＣ１１と、分離された各成分を所定の波長帯内で分析する検出器１２と、検出器１２から出力されるデータを処理するクロマトグラフ用データ処理装置（本実施例の装置。以下、「データ処理装置」とする）１３とを備える。
　データ処理装置１３は、ＣＰＵ（中央演算装置）、記憶装置（メモリ、ハードディスク、ソリッドステートドライブ等）、ディスプレイ、入力装置（キーボード、マウス等）等を有する一般的なコンピュータ（ハードウエア）と、このコンピュータにインストールされた専用のデータ処理ソフトウエアから構成される。このデータ処理装置１３は、パラメータ設定部１３１と、測定試料クロマトグラム作成部１３２としての機能を有する。
　パラメータ設定部１３１は、目的成分の濃度が既知であって互いに相違する複数の標準試料の測定結果を用いて、前記特定指定保持時間及び前記特定指定波長を設定するためのものであり、当該設定のためにユーザが後述の各種条件を入力するためのパラメータ探索条件入力部１３１０、及び、濃度が異なる複数種の標準試料毎に多数の保持時間におけるスペクトルの測定データを記憶する標準試料データ記憶部１３１１を有する。また、パラメータ設定部１３１は、標準試料感度係数算出部１３１２、補正後標準試料クロマトグラム強度算出部１３１３、指標値算出部１３１４、特定指定保持時間・特定指定波長設定部１３１５、及び検量線作成部１３１６を有するが、これら各部に関しては、データ処理装置１３の動作を述べる際に説明する。
　測定試料クロマトグラム作成部１３２は、目的成分の濃度が未知である測定試料の測定結果、並びに前記特定指定保持時間及び前記特定指定波長を用いた処理を行って、最終的に測定試料の補正後クロマトグラムを作成するための、測定試料データ記憶部１３２１、測定試料感度係数算出部１３２２、及び、補正後測定試料クロマトグラム作成部１３２３を有する。測定試料データ記憶部１３２１は、測定試料毎に、多数の保持時間におけるスペクトルの測定データを記憶するものである。測定試料感度係数算出部１３２２及び補正後測定試料クロマトグラム作成部１３２３の機能は、データ処理装置１３の動作を述べる際に説明する。
　データ処理装置１３のディスプレイには、図２に示す定量ブラウザウインドウ２０と、図４に示すデータ解析ウインドウ４０、の２つのウインドウが表示される。また、定量ブラウザウインドウ２０において、ユーザが所定のアイコンをクリックするか、プルダウンメニューの「ウインドウ(W)」中から所定の項目を選択する等の所定の操作を行うことにより、図３に示すパラメータ探索条件入力ウインドウ３０が表示される。以下、定量ブラウザウインドウ２０、パラメータ探索条件入力ウインドウ３０、データ解析ウインドウ４０の順に説明する。
　定量ブラウザウインドウ２０は、標準試料を用いたパラメータ設定の際に用いるウインドウであり、定量結果ビュー２１と、クロマトグラムビュー２２と、検量線／スペクトルビュー２３等の表示領域を有する。定量結果ビュー２１は、目的成分の濃度が既知である複数の標準試料について、標準試料毎に、データファイルの名称、クロマトグラムから求められる面積及び高さ（これら面積及び高さの数値は、後述のデータ処理装置１３によるデータ処理によって、変更され得る）、並びに濃度等の情報が表示される。クロマトグラムビュー２２には、定量結果ビュー２１においてユーザがマウスデバイスなどを用いて選択した１種類の標準試料について、クロマトグラムを表示するものである。ここでクロマトグラムビュー２２は、その中に表示されるタブを切り替えることにより、当該１種類の標準試料に関してデータファイルに収録されている試料の詳細情報を表示することもできる。検量線／スペクトルビュー２３は、タブを切り替えることにより、検量線、又は１つの標準試料及び保持時間におけるスペクトルのいずれか一方を表示するものである。ここで検量線は、標準試料の濃度とクロマトグラムの強度（面積又は高さ）の関係を示すグラフであり、濃度が未知である測定試料のクロマトグラムの強度から該測定試料の濃度を求めるために用いられる。この検量線は、後述のデータ処理装置１３によるデータ処理によって、変更され得る。
　パラメータ探索条件入力ウインドウ３０は、標準試料を用いたパラメータ設定の条件をユーザが入力したうえで、パラメータの自動探索を実行させるためのウインドウである。
　パラメータ探索条件入力ウインドウ３０中の上段には、補正波長設定方法、補正波長（手動）、補正波長（自動）強度、補正波長（自動）移動方向、感度補正スペクトル抽出強度、及びバックグラウンド補正の各入力欄が表示される。補正波長設定方法入力欄には、上述の特定指定波長を自動で設定するか、手動で設定するかをユーザが選択するボタンが表示される。このうち「自動」を選択した場合に、データ処理装置１３は本発明の実施例として機能する。補正波長（手動）入力欄は、補正波長設定方法入力欄において「手動」を選択した場合に、ユーザが特定指定波長の数値を入力することが可能になるものであるが、本発明の実施例とは無関係であるため、詳細な説明は省略する。補正波長（自動）強度入力欄、補正波長（自動）移動方向入力欄、感度補正スペクトル抽出強度入力欄、及びバックグラウンド補正入力欄の詳細は、データ処理装置１３の動作を説明する際に詳述する。
　パラメータ探索条件入力ウインドウ３０中の下段には、標準試料の濃度と、前述の補正後標準試料クロマトグラムの強度の相関性を示す指標値を求めるための計算方法を選択する欄が表示される。本実施例では、検量線の相関係数、偏差の平均値、及び偏差の最大値という3種の指標値から1つ又は複数を選択可能としている。その他、パラメータ探索条件入力ウインドウ３０には、これらの条件の入力が完了した後に、パラメータの自動探索を開始させるための探索開始ボタン等が表示される。
　データ解析ウインドウ４０は、測定試料のデータを処理する際に用いるウインドウであり、等高線ビュー４１、スペクトルビュー４２、クロマトグラムビュー４３等の表示領域を有する。等高線ビュー４１は、横軸に保持時間、縦軸に検出器１２で検出した光の波長を取ったグラフ上において検出器１２の信号強度を等高線で示すものである。スペクトルビュー４２は特定の保持時間における測定試料のスペクトルを表示し、クロマトグラムビュー４３は測定試料の補正後クロマトグラムを表示するものである。
(2) クロマトグラムの飽和・歪みの補正方法
　本実施例のクロマトグラフ用データ処理装置の動作を説明する前に、当該装置で用いる、クロマトグラムにおいて生じる飽和・歪みを補正する方法を、図５を用いて説明する。この補正方法は、標準試料、測定試料のいずれのデータにおいても用いられる。
　一般に、試料における目的成分の濃度が高すぎると、図５(a)に示すように、一定の範囲内にある保持時間T0において、補正前クロマトグラム５１１に飽和や歪みが見られる。このような補正前クロマトグラム５１１の飽和は、図５(b)に示すように、当該保持時間T0におけるスペクトル（「第１スペクトル」とする）５２１に飽和が生じていることに起因する。ここで、１つの目的成分のスペクトルは、成分固有の形状を有することから、飽和・歪みが生じなければ、形状は同じであって強度のみが相違するという相似性を有する。この相似性を利用して、以下のように補正前クロマトグラム５１１の飽和・歪みを補正する。
　まず、補正前クロマトグラム５１１において飽和が生じていない1つの保持時間Tsを指定する（図５(a)）。この保持時間Tsを「指定保持時間」と呼ぶ。その際、指定保持時間Tsを直接指定してもよいが、補正前クロマトグラム５１１の飽和が生じないのはその強度が所定の値以下のときであるため、補正前クロマトグラム５１１の強度Is1を選択することによって、その強度に対応する指定保持時間Tsを指定してもよい。上述のパラメータ探索条件入力ウインドウ３０の感度補正スペクトル抽出強度入力欄で入力する値が、ここで指定する強度Is1に対応する。この指定保持時間Tsにおけるスペクトル（「第２スペクトル」とする）５２２は、飽和や歪みが見られない（図５(b)）。
　続いて、第１スペクトル５２１において飽和や歪みが見られない波長λsを指定する（図５(b)）。この波長λsを「指定波長」と呼ぶ。ここにおいても、指定波長λsを直接指定する代わりに、スペクトルの飽和が生じない所定値以下の強度Is2を選択することによって、その強度に対応する指定波長λsを指定してもよい。上述のパラメータ探索条件入力ウインドウ３０の補正波長（自動）強度入力欄で入力する値が、ここで指定する強度Is2に対応する。
　そして、目的成分においてスペクトルが最大となる目的波長（目的成分が既知のものであれば、一般にこの波長も既知である）λtにおける第２スペクトル５２２の強度である第１強度Iaと、前記指定波長λsにおける第２スペクトル５２２の強度である第２強度Ibを求める（図５(c)）。そして、第１強度Iaと第２強度Ibの比である感度係数K=Ia/Ibを求める。
　上述のスペクトルの相似性により、第２スペクトル５２２における第１強度Iaと第２強度の比は、第１スペクトル５２１の目的波長λtにおける強度と指定波長λsにおける強度の比に等しい。したがって、第１スペクトル５２１において飽和・歪みが生じていない指定波長λsにおけるクロマトグラム５３１を作成したうえで、そのクロマトグラム５３１に感度係数Kを乗じることにより、目的波長λtにおいて飽和・歪みがない、補正後クロマトグラム５３２が得られる。
(3) 本実施例のクロマトグラフ用データ処理装置の動作
　次に、データ処理装置１３の動作を、図６及び図８のフローチャート、並びに図７の検量線の図を用いて説明する。
(3-1) パラメータ設定部１３１の動作
　まず、ユーザが定量ブラウザウインドウ２０において上記所定の操作を行うと、パラメータ探索条件入力部１３１０は、ディスプレイにパラメータ探索条件入力ウインドウ３０を表示し（ステップＳ１）、ユーザが各種条件の入力を行った後に探索開始ボタンを押下するまで待機する（ステップＳ２）。
　パラメータ探索条件入力ウインドウ３０では、ユーザは、補正波長設定方法入力欄において「自動」を選択すると共に、感度補正スペクトル抽出強度入力欄及び補正波長（自動）強度入力欄にそれぞれ、上述の感度補正スペクトル抽出強度Is1及び補正波長（自動）強度Is2の設定値を入力する。但し、これらの値は初期値として入力されるものであって、最終的には、本実施例のデータ処理装置１３により最適値が求められる。また、ユーザは、感度補正スペクトル抽出強度入力欄及び補正波長（自動）強度入力欄に、データ処理装置１３がこの最適値を求める計算を行う際に、感度補正スペクトル抽出強度Is1及び補正波長（自動）強度Is2を変化させる探索範囲及び探索ステップ（間隔）の数値を入力する。それと共に、ユーザは、バックグラウンド補正入力欄に、バックグラウンド補正を行うか否かを入力する。なお、バックグラウンド補正は、クロマトグラムのデータ処理の際に一般に行われているものであるため、ここではその詳細な説明を省略する。ユーザはさらに、パラメータ探索条件入力ウインドウ３０中の下段の欄において、標準試料の濃度と、前述の補正後標準試料クロマトグラムの強度の相関性を示す指標値を求めるための計算方法を選択する。
　なお、補正波長（自動）強度Is2に対応する波長が、目的波長λtよりも長波長側及び短波長側の双方に存在する場合がある。その場合には、ユーザは補正波長（自動）移動方向入力欄に、いずれの波長を選択するかを入力する。また、感度補正スペクトル抽出強度Is1の代わりに1つの指定保持時間Tsを入力するようにしてもよいし、補正波長（自動）強度Is2の代わりに1つの指定波長λsを入力するようにしてもよい。
　これらの入力が完了した後に、ユーザが探索開始ボタンを押下することにより、ステップＳ３に移る。ステップＳ３では、標準試料感度係数算出部１３１２が標準試料データ記憶部１３１１から、複数種の標準試料について、標準試料毎にスペクトルのデータを取得する。次に、標準試料感度係数算出部１３１２は、１組の指定保持時間Ts及び指定波長λsを用いて、標準試料毎に、(2)で説明した感度係数Kの計算を行う（ステップＳ４）。続いて、補正後標準試料クロマトグラム強度算出部１３１３は、標準試料毎に、ステップＳ４で得られた感度係数Kに基づいて、(2)で説明したように補正後クロマトグラム５３２を作成し、当該補正後クロマトグラム５３２のピーク強度（ピークの面積強度、又はピークトップの強度）を算出する（ステップＳ５）。
　ここまでの動作により、１組の指定保持時間Ts及び指定波長λsに対して、n種の標準試料の濃度diと標準試料の補正後クロマトグラム５３２のピーク強度Iiのデータ列{(di, Ii)}(i=1, 2, …n)が得られる。そこで、指標値算出部１３１４は、当該１組の指定保持時間Ts及び指定波長λsにおける指標値を算出する（ステップＳ６）。本実施例では、指標値は、検量線の相関係数、偏差の平均値、及び偏差の中央値のうち、ユーザがパラメータ探索条件入力ウインドウ３０で選択したものを算出する。これらのうち検量線の相関係数は上述の相関係数Cの算出式を用いて算出し、偏差の平均値や中央値は周知の方法によって偏差を求めたうえでそれら平均値や中央値を算出する。相関係数Cは大きいほど相関性が高く、偏差の平均値や中央値は小さいほど相関性が高いことを意味する。
　こうして１組の指定保持時間Ts及び指定波長λsに対して指標値が求められた後、指標値算出部１３１４は、パラメータ探索条件入力ウインドウ３０にユーザが入力した探索範囲及び探索ステップに対応する指定保持時間Ts及び指定波長λsの全て組み合わせについて、指標値の算出が完了したか否かを確認する（ステップＳ７）。それが完了したならばステップＳ９に移る。一方、それが完了していなければ、指定保持時間Ts及び指定波長λsを指標値の算出が完了していないものに変更し（ステップＳ８）、ステップＳ４～Ｓ６の操作を繰り返す。
　ステップＳ９では、特定指定保持時間・特定指定波長設定部１３１５は、これまでに得られた複数組の指定保持時間Ts及び指定波長λsの組み合わせにおいて、指標値が、相関性が最も高いことを示す指定保持時間Ts及び指定波長λsを、特定指定保持時間及び特定指定波長として選択する。なお、指標値が所定の基準を超えた指定保持時間Ts及び指定波長λsの組み合わせから、任意の１組を特定指定保持時間及び特定指定波長として選択するようにしてもよい。また、特定指定保持時間及び特定指定波長に対応する強度値、並びに相関係数の計算結果を探索条件入力ウインドウ３０中の該当欄に表示する。
　続いて、検量線作成部１３１６は、決定された特定指定保持時間及び特定指定波長において求めたデータ列{(di, Ii)}に基づいて検量線を作成する。クロマトグラムの補正前には、クロマトグラムのピークの飽和に起因して、データ点が一次関数の直線からずれる（図７(a)）のに対して、ここで得られた検量線は、データ点が一次関数の直線によく一致する（図７(b)）。検量線作成部１３１６は、定量ブラウザウインドウ２０の検量線／スペクトルビュー２３に、こうして得られた補正後の検量線を表示する。
　これにより、パラメータ設定部１３１の動作を終了する。
(3-2) 測定試料クロマトグラム作成部１３２の動作
　ユーザがデータ解析ウインドウ４０において、目的成分の濃度が未知である測定試料の１つを特定する所定の操作を行うと、測定試料クロマトグラム作成部１３２は動作を開始する。まず、測定試料感度係数算出部１３２２は、測定試料データ記憶部１３２１から、当該測定試料のスペクトルのデータを取得する（ステップＳ２１）。次に、測定試料感度係数算出部１３２２は、パラメータ設定部１３１において選択された特定指定保持時間をTs、特定指定波長をλsとして、(2)で説明した感度係数Kの計算を行う（ステップＳ２２）。続いて、補正後測定試料クロマトグラム作成部１３２３は、ステップＳ２３で得られた感度係数Kに基づいて、当該測定試料における補正後クロマトグラム５３２を、(2)で説明したように作成する（ステップＳ２３）。これにより測定試料クロマトグラム作成部１３２の動作を終了する。
　本発明は上記実施例のものには限られない。
　例えば、上記実施例では、パラメータ設定部１３１において、特定指定保持時間及び特定指定波長において得られたデータに基づいて検量線作成部１３１６が検量線を作成するが、単に特定指定保持時間及び特定指定波長を選択するという目的においては、検量線作成部１３１６は必須ではない。その反対に、クロマトグラムの飽和や歪みの影響のない検量線を作成することを目的とするとする場合には、図９に示す分析システム１０Ａのように、パラメータ設定部１３１を有し、測定試料クロマトグラム作成部１３２を有しない構成を取ることもできる。
１０、１０Ａ…分析システム
１１…ＬＣ
１２…検出器
１３…データ処理装置
１３１…パラメータ設定部
１３１０…パラメータ探索条件入力部
１３１１…標準試料データ記憶部
１３１２…標準試料感度係数算出部
１３１３…補正後標準試料クロマトグラム強度算出部
１３１４…指標値算出部
１３１５…特定指定保持時間・特定指定波長設定部
１３１６…検量線作成部
１３２…測定試料クロマトグラム作成部
１３２１…測定試料データ記憶部
１３２２…測定試料感度係数算出部
１３２３…補正後測定試料クロマトグラム作成部
２０…定量ブラウザウインドウ
２１…定量結果ビュー
２２…クロマトグラムビュー
２３…検量線／スペクトルビュー
３０…パラメータ探索条件入力ウインドウ
４０…データ解析ウインドウ
４１…等高線ビュー
４２…スペクトルビュー
４３…クロマトグラムビュー
５１１…補正前クロマトグラム
５２１…第１スペクトル
５２２…第２スペクトル
５３１…指定波長λsにおけるクロマトグラム
５３２…補正後クロマトグラム

Claims (6)

  1. 　目的成分に対応した目的波長を含む波長範囲内で取得したスペクトルの時間変化に基づいてクロマトグラムを作成するためのクロマトグラフ用データ処理装置であって、
    　a-1) 前記目的成分の濃度が異なる複数の標準試料について、標準試料毎に前記スペクトルのデータを記憶する標準試料データ記憶手段と、
    　a-2) 前記標準試料毎に、各標準試料に共通の指定保持時間における該標準試料のスペクトルの、前記目的波長における強度である第１強度を、該目的波長とは異なる前記波長範囲内の指定波長における強度である第２強度で除した標準試料感度係数を算出する標準試料感度係数算出手段と、
    　a-3) 前記標準試料毎に、前記指定波長におけるクロマトグラムに前記標準試料感度係数を乗算した補正後標準試料クロマトグラムのピーク強度を算出する補正後標準試料クロマトグラム強度算出手段と、
    　a-4) 前記標準試料の濃度と前記補正後標準試料クロマトグラムのピーク強度の相関性を示す指標値を求める指標値算出手段と、
    　a-5) 前記指定保持時間及び前記指定波長のいずれか一方又は両方を変化させ、前記指標値が、相関性が高いことを示す所定の範囲内にある１組の該指定保持時間及び該指定波長を選択して特定指定保持時間及び特定指定波長を設定する特定指定保持時間・特定指定波長設定手段と、
    　b-1) 前記目的成分の濃度が未知である測定試料についてスペクトルのデータを記憶する測定試料データ記憶手段と、
    　b-2) 前記特定指定保持時間における前記測定試料のスペクトルの、前記目的波長における強度である第１強度を、前記特定指定波長における強度である第２強度で除した測定試料感度係数を算出する測定試料感度係数算出手段と、
    　b-3) 前記特定指定波長におけるクロマトグラムに前記測定試料感度係数を乗算した補正後測定試料クロマトグラムを作成する補正後測定試料クロマトグラム作成手段と
    を備えることを特徴とするクロマトグラフ用データ処理装置。
  2. 　目的成分に対応した目的波長を含む波長範囲内で取得したスペクトルの時間変化に基づいてクロマトグラムを作成する際に用いるパラメータを特定するためのクロマトグラフ用データ処理装置であって、
    　1) 前記目的成分の濃度が異なる複数の標準試料について、標準試料毎に前記スペクトルのデータを記憶する標準試料データ記憶手段と、
    　2) 前記標準試料毎に、各標準試料に共通の指定保持時間における該標準試料のスペクトルの、前記目的波長における強度である第１強度を、該目的波長とは異なる前記波長範囲内の指定波長における強度である第２強度で除した標準試料感度係数を算出する標準試料感度係数算出手段と、
    　3) 前記標準試料毎に、前記指定波長におけるクロマトグラムに前記標準試料感度係数を乗算した補正後標準試料クロマトグラムのピーク強度を算出する補正後標準試料クロマトグラム強度算出手段と、
    　4) 前記標準試料の濃度と前記補正後標準試料クロマトグラムのピーク強度の相関性を示す指標値を求める指標値算出手段と、
    　5) 前記指定保持時間及び前記指定波長のいずれか一方又は両方を変化させ、前記指標値が、相関性が高いことを示す所定の範囲内にある１組の該指定保持時間及び該指定波長を選択して特定指定保持時間及び特定指定波長を設定する特定指定保持時間・特定指定波長設定手段と
    を備えることを特徴とするクロマトグラフ用データ処理装置。
  3. 　前記特定指定保持時間及び前記特定指定波長における前記標準試料の濃度と前記補正後標準試料クロマトグラムのピーク強度に基づいて、検量線を作成する検量線作成手段を備えることを特徴とする請求項１又は２に記載のクロマトグラフ用データ処理装置。
  4. 　前記指標値が、目的成分の濃度が異なるn種（nは2以上の自然数）の標準試料における濃度diと補正後標準試料クロマトグラムのピーク強度Iiのデータ列{(di, Ii)}(i=1, 2, …n)について、
    Figure JPOXMLDOC01-appb-M000001
     で与えられる相関係数Cであることを特徴とする請求項１～３のいずれかに記載のクロマトグラフ用データ処理装置。
  5. 　前記指標値が、目的成分が異なるn種（nは2以上の自然数）の標準試料における濃度diと補正後標準試料クロマトグラムのピーク強度Iiのデータ列{(di, Ii)}(i=1, 2, …n)に基づいて求められる濃度dとピーク強度Iの１次関数I=kd（kは係数）について、標準試料毎の該１次関数の値とkdiとピーク強度Iの実測値Iiとの差である偏差の絶対値の平均値又は最大値であることを特徴とする請求項１～４のいずれかに記載のクロマトグラフ用データ処理装置。
  6. 　目的成分に対応した目的波長を含む波長範囲内で取得したスペクトルの時間変化に基づいてクロマトグラムを作成するためのクロマトグラフ用データ処理方法であって、
    　前記目的成分の濃度が異なる複数の標準試料について、該標準試料毎に、各標準試料に共通の指定保持時間における該標準試料のスペクトルの、前記目的波長における強度である第１強度を、該目的波長とは異なる前記波長範囲内の指定波長における強度である第２強度で除した標準試料感度係数を算出し、
    　前記標準試料毎に、前記指定波長におけるクロマトグラムに前記標準試料感度係数を乗算した補正後標準試料クロマトグラムのピーク強度を算出し、
    　前記標準試料の濃度と前記補正後標準試料クロマトグラムのピーク強度の相関性を示す指標値を求め、
    　前記指定保持時間及び前記指定波長のいずれか一方又は両方を変化させながら、前記指標値が、相関性が高いことを示す所定の範囲内にある１組の該指定保持時間及び該指定波長を選択することにより特定指定保持時間及び特定指定波長を設定し、
    　前記目的成分の濃度が未知である測定試料について、前記特定指定保持時間における該測定試料のスペクトルの、前記目的波長における強度である第１強度を、前記特定指定波長における強度である第２強度で除した測定試料感度係数を算出し、
    　前記特定指定波長におけるクロマトグラムに前記測定試料感度係数を乗算した補正後測定試料クロマトグラムを作成する
    ことを特徴とするクロマトグラフ用データ処理方法。
PCT/JP2013/073555 2013-09-02 2013-09-02 クロマトグラフ用データ処理装置及び方法 WO2015029254A1 (ja)

Priority Applications (5)

Application Number Priority Date Filing Date Title
CN201380079317.2A CN105518454B (zh) 2013-09-02 2013-09-02 色谱仪用数据处理装置以及方法
PCT/JP2013/073555 WO2015029254A1 (ja) 2013-09-02 2013-09-02 クロマトグラフ用データ処理装置及び方法
JP2015533921A JP6065981B2 (ja) 2013-09-02 2013-09-02 クロマトグラフ用データ処理装置及び方法
US14/915,794 US10151734B2 (en) 2013-09-02 2013-09-02 Data processing system and method for chromatograph
CN201810180594.7A CN108508122B (zh) 2013-09-02 2013-09-02 色谱仪用数据处理装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2013/073555 WO2015029254A1 (ja) 2013-09-02 2013-09-02 クロマトグラフ用データ処理装置及び方法

Publications (1)

Publication Number Publication Date
WO2015029254A1 true WO2015029254A1 (ja) 2015-03-05

Family

ID=52585860

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2013/073555 WO2015029254A1 (ja) 2013-09-02 2013-09-02 クロマトグラフ用データ処理装置及び方法

Country Status (4)

Country Link
US (1) US10151734B2 (ja)
JP (1) JP6065981B2 (ja)
CN (2) CN105518454B (ja)
WO (1) WO2015029254A1 (ja)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108780073A (zh) * 2016-03-01 2018-11-09 株式会社岛津制作所 色谱仪装置
CN114428127A (zh) * 2020-10-29 2022-05-03 中国石油化工股份有限公司 一种鉴别石油产品的方法
KR102546894B1 (ko) * 2022-07-15 2023-06-23 주식회사 위드텍 크로마토그래피 기반의 단일물질 표준용액 측정 시스템

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6269803B2 (ja) * 2016-12-26 2018-01-31 株式会社島津製作所 クロマトグラフ用データ処理装置及び方法
JP7012998B2 (ja) * 2017-09-21 2022-01-31 株式会社日立ハイテクサイエンス クロマトグラフのデータ処理装置
EP3805748A4 (en) * 2018-05-30 2021-06-23 Shimadzu Corporation SPECTRAL DATA PROCESSING DEVICE AND ANALYSIS DEVICE
JP6973323B2 (ja) * 2018-08-08 2021-11-24 株式会社島津製作所 電子線マイクロアナライザー、データ処理方法及びデータ処理プログラム
CN109541100B (zh) * 2018-12-13 2021-06-29 安徽皖仪科技股份有限公司 多通道波长的信号漂移处理方法、装置及多通道检测器
US20220155270A1 (en) * 2019-03-22 2022-05-19 Shimadzu Corporation Sample Measurement Device and Measurement Sample Identification Method
CN110068646A (zh) * 2019-05-21 2019-07-30 苏州芷宁信息科技有限公司 有机混合物的等级评定方法和装置
EP3742161B1 (en) 2019-05-24 2023-07-19 Sartorius Stedim Biotech GmbH Chromatography method, method of determining the concentration of at least one compound in a chromatography method and method of obtaining at least one chromatography method parameter
EP3742159A1 (en) * 2019-05-24 2020-11-25 Sartorius Stedim Biotech GmbH Chromatography method, method of determining the influence of the interdependency of at least two parameters in a chromatography method and method of obtaining at least one chromatography method parameter
CN110632228B (zh) * 2019-09-24 2022-01-11 宁夏医科大学 一种自动实现uplc-hrms中保留时间漂移校正与化合物注册的方法
WO2021113977A1 (en) * 2019-12-13 2021-06-17 Etrace Medical Diagnostics Inc. Anomaly detection and diagnosis in chromatography applications
CN113567602A (zh) * 2020-04-28 2021-10-29 株式会社岛津制作所 检测方法以及检测设备
JP2023158945A (ja) * 2022-04-19 2023-10-31 アークレイ株式会社 制御装置、分離分析装置、分離分析方法、及び分離分析プログラム
CN115792067B (zh) * 2023-02-07 2023-04-14 河北对外经贸职业学院 一种基于工业杀菌剂毒性检测的计算机数据分析方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH04303768A (ja) * 1991-03-30 1992-10-27 Shimadzu Corp クロマトグラフ用データ処理装置
JPH10300741A (ja) * 1997-04-22 1998-11-13 Sekisui Chem Co Ltd 糖化ヘモグロビンの測定方法
JPH11326304A (ja) * 1998-05-12 1999-11-26 Shimadzu Corp クロマトグラフ用データ処理装置
JP2001056325A (ja) * 1999-08-20 2001-02-27 Shimadzu Corp ガスクロマトグラフ装置
JP2005077094A (ja) * 2003-08-29 2005-03-24 Shimadzu Corp クロマトグラフ質量分析装置
JP2010054406A (ja) * 2008-08-29 2010-03-11 Shimadzu Corp クロマトグラフ質量分析データ解析方法及び装置
JP2011185743A (ja) * 2010-03-09 2011-09-22 Shimadzu Corp クロマトグラムピーク純度判定装置

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1016737B (zh) * 1985-04-01 1992-05-20 株式会社岛津制作所 液体色谱用多通道检测数据处理装置
US4674323A (en) * 1985-07-31 1987-06-23 The Dow Chemical Company Self-diagnostic gel permeation/size exclusion chromatograph
US4941101A (en) * 1987-09-01 1990-07-10 Hewlett-Packard Company Method for analyzing chromatograms
JPH1172489A (ja) * 1997-08-28 1999-03-16 Hitachi Ltd 液体クロマトグラフ
JP5279485B2 (ja) * 2006-03-24 2013-09-04 アークレイ株式会社 グリコヘモグロビン濃度測定方法および濃度測定装置
JP5590145B2 (ja) * 2010-11-30 2014-09-17 株式会社島津製作所 質量分析データ処理装置
JP5930066B2 (ja) * 2012-12-07 2016-06-08 株式会社島津製作所 クロマトグラフ用データ処理装置及びデータ処理方法
GB201304631D0 (en) * 2013-03-14 2013-05-01 Malvern Instr Ltd Monomer detection in protein separation
CN105308451B (zh) * 2013-08-05 2017-12-08 株式会社岛津制作所 色谱仪用数据处理装置及色谱仪用数据处理方法
US9606088B2 (en) * 2014-03-17 2017-03-28 Prism Analytical Technologies, Inc. Process and system for rapid sample analysis

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH04303768A (ja) * 1991-03-30 1992-10-27 Shimadzu Corp クロマトグラフ用データ処理装置
JPH10300741A (ja) * 1997-04-22 1998-11-13 Sekisui Chem Co Ltd 糖化ヘモグロビンの測定方法
JPH11326304A (ja) * 1998-05-12 1999-11-26 Shimadzu Corp クロマトグラフ用データ処理装置
JP2001056325A (ja) * 1999-08-20 2001-02-27 Shimadzu Corp ガスクロマトグラフ装置
JP2005077094A (ja) * 2003-08-29 2005-03-24 Shimadzu Corp クロマトグラフ質量分析装置
JP2010054406A (ja) * 2008-08-29 2010-03-11 Shimadzu Corp クロマトグラフ質量分析データ解析方法及び装置
JP2011185743A (ja) * 2010-03-09 2011-09-22 Shimadzu Corp クロマトグラムピーク純度判定装置

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108780073A (zh) * 2016-03-01 2018-11-09 株式会社岛津制作所 色谱仪装置
CN114428127A (zh) * 2020-10-29 2022-05-03 中国石油化工股份有限公司 一种鉴别石油产品的方法
KR102546894B1 (ko) * 2022-07-15 2023-06-23 주식회사 위드텍 크로마토그래피 기반의 단일물질 표준용액 측정 시스템

Also Published As

Publication number Publication date
CN105518454B (zh) 2018-04-10
US20160216240A1 (en) 2016-07-28
CN108508122A (zh) 2018-09-07
US10151734B2 (en) 2018-12-11
JPWO2015029254A1 (ja) 2017-03-02
JP6065981B2 (ja) 2017-01-25
CN108508122B (zh) 2020-07-14
CN105518454A (zh) 2016-04-20

Similar Documents

Publication Publication Date Title
JP6065981B2 (ja) クロマトグラフ用データ処理装置及び方法
JP6620894B2 (ja) クロマトグラフ質量分析用データ解析装置
JP6132067B2 (ja) クロマトグラフ質量分析装置用データ処理装置及びプログラム
JP6573028B2 (ja) データ処理装置
JP5930066B2 (ja) クロマトグラフ用データ処理装置及びデータ処理方法
US20240011956A1 (en) System and Method for Background Removal in Spectrometry System
JP5954497B2 (ja) クロマトグラフ用データ処理装置及びデータ処理方法
JP6048373B2 (ja) クロマトグラフ用データ処理装置及びデータ処理方法
JP6269803B2 (ja) クロマトグラフ用データ処理装置及び方法
JP6226823B2 (ja) クロマトグラフ質量分析装置及びその制御方法
JP6583159B2 (ja) クロマトグラフによる濃度測定方法及びクロマトグラフ
JP2016017941A5 (ja)
JP2017223670A (ja) 液体クロマトグラフおよび液体クロマトグラフの検出器出力値変動補正手段

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 13892511

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2015533921

Country of ref document: JP

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 14915794

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 13892511

Country of ref document: EP

Kind code of ref document: A1