CN113567604A - Detection and analysis method of chromatographic spectrogram and electronic equipment - Google Patents
Detection and analysis method of chromatographic spectrogram and electronic equipment Download PDFInfo
- Publication number
- CN113567604A CN113567604A CN202110830534.7A CN202110830534A CN113567604A CN 113567604 A CN113567604 A CN 113567604A CN 202110830534 A CN202110830534 A CN 202110830534A CN 113567604 A CN113567604 A CN 113567604A
- Authority
- CN
- China
- Prior art keywords
- point
- peak
- chromatographic
- detection
- curve
- 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.)
- Granted
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 102
- 238000004458 analytical method Methods 0.000 title claims abstract description 30
- 238000000034 method Methods 0.000 claims abstract description 55
- 238000012937 correction Methods 0.000 claims abstract description 20
- 238000012545 processing Methods 0.000 claims abstract description 11
- 238000003909 pattern recognition Methods 0.000 claims abstract description 8
- 238000004422 calculation algorithm Methods 0.000 claims description 17
- 238000011156 evaluation Methods 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 11
- 238000011208 chromatographic data Methods 0.000 claims description 6
- 238000007619 statistical method Methods 0.000 claims description 4
- 238000001228 spectrum Methods 0.000 claims description 2
- 238000003556 assay Methods 0.000 claims 3
- 238000010205 computational analysis Methods 0.000 claims 1
- 230000003595 spectral effect Effects 0.000 abstract 1
- 230000014759 maintenance of location Effects 0.000 description 10
- 230000000875 corresponding effect Effects 0.000 description 5
- 230000002159 abnormal effect Effects 0.000 description 4
- 238000004587 chromatography analysis Methods 0.000 description 4
- 230000007547 defect Effects 0.000 description 4
- 239000012071 phase Substances 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 239000007791 liquid phase Substances 0.000 description 3
- AYFVYJQAPQTCCC-GBXIJSLDSA-N L-threonine Chemical compound C[C@@H](O)[C@H](N)C(O)=O AYFVYJQAPQTCCC-GBXIJSLDSA-N 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 230000005856 abnormality Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
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/8627—Slopes
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Treatment Of Liquids With Adsorbents In General (AREA)
Abstract
The application relates to a detection and analysis method of a chromatographic spectrogram and electronic equipment, wherein the method comprises the steps of adopting a time window to sequentially slide from the starting point of a chromatographic curve to detect spectral peaks aiming at obtained original spectrogram data until all the chromatographic peaks in the curve are detected; wherein, in the detection processing process of each chromatographic peak, the method comprises the following steps: comparing the slope of each point on the curve with a threshold value, and determining a reference point of a chromatographic peak based on the comparison result; performing pattern recognition detection on the chromatographic curve near the position by taking the Gaussian wave as a matching wave, and taking the peak top in a recognition detection result as a peak top correction point; correcting the temporary peak point reference point in the reference points according to the peak point correction point to obtain a peak point reference point; and determining the chromatographic peak characterized by the peak top point reference point, the peak starting point reference point and the peak end point reference point in the reference points as the detected chromatographic peak. The spectrogram detection and analysis method can realize better comprehensive performance.
Description
Technical Field
The application belongs to the technical field of chromatographic analysis, and particularly relates to a detection and analysis method of a chromatographic spectrogram and electronic equipment.
Background
In the related art, a computer is used for qualitative and quantitative analysis of a chromatographic spectrum in chromatographic analysis, wherein the important basis of the analysis is to identify chromatographic peaks, and the peak types of the chromatographic peaks are complex and diverse, including overlapping peaks, front shoulder peaks, rear shoulder peaks, trailing peaks, negative peaks and the like.
At present, the detection and identification of chromatographic peaks generally adopt a time window method and a derivative method. The time window method obtains the chromatographic peak position by finding the extreme value in the retention time range of the components; the derivative method obtains a series of extreme values by taking the derivative of the original chromatographic data, and then finds the chromatographic peak position based on the retention time range of the components.
In the implementation process of the time window method and the derivative method, amplitude threshold values are required to be combined to judge whether the chromatographic peaks are real chromatographic peaks, so that whether the threshold values are set reasonably is critical, if the threshold values are set too large, the real chromatographic peaks can be missed, the resolution of chromatographic components is reduced, and if the threshold values are set too small, small peaks formed by some noises can be regarded as chromatographic peaks, and misjudgment can be caused. And the time window method and the derivative both need the retention time range of the components to search the chromatographic peak position, so that if the peak shifts out of the range, the problem of incorrect found chromatographic peak or component misjudgment can be caused. In addition, the method of using the first derivative and the second derivative to be combined finds that the recognition rate of chromatographic peaks to a single peak is higher, but the recognition rate is greatly reduced when multiple peaks are overlapped, particularly when a shoulder peak appears.
In the related art, there is a method for identifying a chromatographic peak based on pattern matching, which is independent of retention time of chromatographic components when identifying the peak and is insensitive to noise and changes in width and amplitude of the chromatographic peak, but since the algorithm calculates similarity between a matching waveform and a chromatographic curve by using a statistical method, i.e., a correlation coefficient, a threshold value of correlation is important, which may result in poor universality of the pattern matching method.
The above is only for the purpose of assisting understanding of the technical aspects of the present invention, and does not represent an admission that the above is prior art.
Disclosure of Invention
In order to overcome the problems in the related art at least to a certain extent, the application provides a method for detecting and analyzing a chromatographic spectrogram and an electronic device, a chromatographic peak is detected based on pattern matching and the slope of a chromatographic curve, and analysis is performed based on the detected chromatographic peak.
In order to achieve the purpose, the following technical scheme is adopted in the application:
in a first aspect,
the application provides a detection and analysis method of a chromatographic spectrogram, which comprises the following steps:
acquiring original spectrogram data to be processed;
for the original spectrogram data, sequentially sliding from the starting point of a chromatogram curve by adopting a time window with a preset size to perform detection processing of chromatographic peaks until the original spectrogram data is processed, detecting all chromatographic peaks in the curve, and generating a detection result based on the detected chromatographic peaks; wherein, in the detection processing process of each chromatographic peak, the method comprises the following steps:
comparing the slopes of the points on the chromatographic curve with a threshold value, determining reference points of chromatographic peaks based on the comparison result, wherein the reference points comprise a peak start point reference point, a temporary peak top point reference point and a peak end point reference point, and the threshold value is determined based on the calculation and analysis of the no-load output signal of the chromatographic instrument generating the original spectrogram data;
based on the position of the chromatographic peak represented by the reference point on the chromatographic curve, performing pattern recognition detection on the chromatographic curve near the position by taking a Gaussian wave as a matching wave, and taking the peak top in a recognition detection result as a peak top correction point;
correcting the temporary peak point reference point according to the peak point correction point, and taking the corrected point as a peak point reference point;
and determining the chromatographic peak characterized by the peak starting point reference point, the peak end point reference point and the peak top point reference point as the detected chromatographic peak.
Optionally, the comparing the slope of each point on the chromatographic curve with a threshold value, and determining the reference point of the chromatographic peak based on the comparison result includes:
calculating the slope of a point on the chromatographic curve during the sliding of the time window in real time,
comparing the calculated slope with a predetermined first threshold, comparing ordinate values of two consecutive points when the slopes of the two points are greater than the first threshold, determining a point having a smaller ordinate value as the peak start point reference point, and
analyzing and judging the positive and negative changes of the slope of the point after the peak starting point reference point, when the slope of one point is negative and the slope of the point before the point is positive, comparing the vertical coordinate values of the two points, determining the point with the larger vertical coordinate value as the temporary peak top point reference point, and
and comparing the slope of the point behind the temporary peak top point reference point with a predetermined second threshold, when the slopes of two continuous points are smaller than the second threshold, comparing the ordinate values of the two points, and determining the point with the smaller ordinate value as the peak terminal point reference point.
Optionally, the process of computing and analyzing an idle output signal of a chromatography instrument generating the raw spectrogram data includes:
carrying out statistical analysis on the slope change of an output baseline of the chromatographic instrument in no-load, calculating and determining the variance of the slope change, and further determining the standard deviation of the slope change;
three times the standard deviation of the slope change is taken as the first threshold, and minus three times the standard deviation of the slope change is taken as the second threshold.
Optionally, the correcting the temporary peak top reference point according to the peak top correction point, and taking the corrected point as the peak top reference point specifically includes:
and comparing the longitudinal coordinate values of the peak top point correction point and the temporary peak top point reference point, and determining a point with a larger longitudinal coordinate value as a peak top point reference point.
Optionally, the performing pattern recognition detection on the chromatographic curve near the position by using a gaussian wave as a matching wave, and using a peak top in a recognition detection result as a peak top correction point includes:
taking the chromatographic curve near the position as a curve to be detected, sliding the waveform of the Gaussian wave on the curve to be detected from a left end point to a right end point, simultaneously calculating the correlation coefficient of the chromatographic curve and the Gaussian wave, and obtaining a correlation coefficient array of chromatographic data of the curve to be detected relative to the Gaussian wave based on the calculation result;
and comparing and analyzing each correlation coefficient in the correlation coefficient group with a preset value, determining the position of a Gaussian peak position based on the correlation coefficient of which the coefficient value is greater than the preset value, and determining a point at the position on the curve to be detected as the peak top point correction point.
Optionally, the generating a detection result based on the detected chromatographic peak comprises:
and integrating the detected chromatographic peaks, and calculating to determine the areas and the heights of the chromatographic peaks.
Optionally, the method further comprises the step of,
detecting and processing the original spectrogram data by adopting a reference chromatogram spectrogram detection algorithm to obtain a reference detection result;
and comparing and analyzing the detection result with the reference detection result, generating a detection evaluation report, and displaying and outputting the detection evaluation report.
Optionally, the comparing and analyzing the detection result and the reference detection result includes:
matching the chromatographic peaks detected in the detection result and the reference detection result, determining the matched chromatographic peaks, and generating a detection evaluation report based on the proportion of the matched chromatographic peaks in the detected chromatographic peaks and the difference of the matched chromatographic peaks.
In a second aspect of the present invention,
the application provides an electronic device, including:
a memory having an executable program stored thereon;
a processor for executing the executable program in the memory to implement the steps of the method described above.
This application adopts above technical scheme, possesses following beneficial effect at least:
according to the technical scheme, the liquid phase chromatogram and the gas phase chromatogram are actually and specifically detected based on the mode matching and by combining the slope of the chromatogram curve, so that the detection reliability is integrally improved. In the method, the threshold value in the slope detection mode is automatically calculated and determined based on the self signal of the instrument, and the peak characteristic point detected by mode matching is only used as a reference point to correct the peak characteristic point obtained based on slope detection, so that the defects of the two modes in the prior art are overcome.
Additional advantages, objects, and features of the invention will be set forth in part in the description which follows and in part will become apparent to those having ordinary skill in the art upon examination of the following or may be learned from practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
Drawings
The accompanying drawings are included to provide a further understanding of the technology or prior art of the present application and are incorporated in and constitute a part of this specification. The drawings expressing the embodiments of the present application are used for explaining the technical solutions of the present application, and should not be construed as limiting the technical solutions of the present application.
FIG. 1 is a schematic flowchart illustrating a method for detecting and analyzing a chromatogram according to an embodiment of the present application;
FIG. 2 is a schematic illustration of a process for detection of a chromatographic peak in the embodiment of FIG. 1;
FIG. 3 is a schematic illustration of a process for performing peak matching between different test results according to an embodiment of the present disclosure;
FIG. 4 is a schematic illustration of a flow chart for displaying and outputting a detection and evaluation report according to an embodiment of the present application;
fig. 5 is a schematic structural diagram of an electronic device according to an embodiment of the present application.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application more apparent, the technical solutions of the present application will be described in detail below. It is to be understood that the embodiments described are only a few embodiments of the present application and not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the examples given herein without making any creative effort, shall fall within the protection scope of the present application.
As described in the background art, in the related art of detecting and analyzing a chromatogram, methods such as a time window method, a derivative method, and pattern matching have the defects of good multimodal overlap identification or poor universality.
In view of the above, the present application provides a method for detecting and analyzing a chromatogram, which detects a chromatogram peak based on pattern matching and combining a slope of a chromatogram curve, and the method helps to make up for the deficiencies in the prior art and realize the detection and analysis of the chromatogram with better comprehensive performance.
As shown in fig. 1 and fig. 2, in an embodiment, the method for detecting and analyzing a chromatogram provided by the present application includes:
step S110, acquiring original spectrogram data to be processed;
step S120, for original spectrogram data, adopting a preset size time window (the size of the window is set based on a standard retention time parameter of an instrument, for example, the window is 10% of the standard retention time parameter) to sequentially slide from the starting point of a chromatographic curve to perform chromatographic peak detection processing, detecting all chromatographic peaks in the curve until the original spectrogram data is processed, and generating a detection result based on the detected chromatographic peaks;
it is easy to understand that, for a certain complete spectrogram data, there are usually a plurality of chromatographic peaks, in other words, the detection process of step S120 is performed by sliding the time window to sequentially detect the individual chromatographic peaks. Specifically, as shown in fig. 2, in step S120, in the process of detecting each chromatographic peak, the method includes:
step S121, comparing the slope of each point on the chromatographic curve with a threshold value, and determining a reference point of a chromatographic peak based on the comparison result, wherein the reference point comprises a peak starting point reference point, a temporary peak top point reference point and a peak end point reference point, and the threshold value is determined based on the calculation and analysis of a no-load output signal of a chromatographic instrument for generating original spectrogram data;
in step S121, unlike the prior art, in the process of performing detection based on the slope, the threshold is not manually set, but determined based on calculation and analysis of the no-load output signal of the chromatography apparatus that generates the original spectrogram data.
Step S122, based on the position of the chromatographic peak represented by the reference point on the chromatographic curve, performing pattern recognition detection on the chromatographic curve near the position by taking a Gaussian wave as a matching wave, and taking the peak top in the recognition detection result as a peak top correction point;
step S123, correcting the temporary peak point reference point according to the peak point correction point, and taking the corrected point as the peak point reference point;
specifically, in this step, by comparing the vertical coordinate values of the peak-top correction point and the temporary peak-top reference point, a point having a larger vertical coordinate value is determined as the peak-top reference point.
And step S124, determining the chromatographic peak represented by the peak starting point reference point, the peak end point reference point and the peak top point reference point as the detected chromatographic peak.
In step S120, the time window is slid, and steps S121 to S124 are repeated in each chromatographic peak detection process until all raw spectrogram data are processed.
According to the technical scheme, the liquid phase chromatogram and the gas phase chromatogram are actually and specifically detected based on the mode matching and by combining the slope of the chromatogram curve, so that the detection reliability is integrally improved. In the method, the threshold value in the slope detection mode is automatically calculated and determined based on the self signal of the instrument, and the peak characteristic point detected by mode matching is only used as a reference point to correct the peak characteristic point obtained based on slope detection, so that the defects of the two modes in the prior art are overcome.
To facilitate understanding of the technical solutions of the present application, the technical solutions of the present application will be described below with reference to another embodiment.
In this embodiment, similarly, step S210 is performed first, and raw spectrogram data to be processed is acquired.
Then, step S220 is performed, for the original spectrogram data, a preset size time window is adopted to sequentially slide from the starting point of the chromatographic curve to perform chromatographic peak detection processing, until the original spectrogram data is processed, all chromatographic peaks in the curve are detected, and a detection result is generated based on the detected chromatographic peaks;
for example, generating the detection result based on the detected chromatographic peak includes integrating the detected chromatographic peak and computationally determining the area and height of the chromatographic peak.
Similarly, in this embodiment, in the step S220, during the detection process of each chromatographic peak, the method includes:
step S221, comparing the slope of each point on the chromatogram curve with a threshold value, and determining a reference point of a chromatogram peak based on the comparison result, wherein the reference point comprises a peak starting point reference point, a temporary peak top point reference point and a peak end point reference point, and the threshold value is determined based on the calculation and analysis of the no-load output signal of the chromatogram instrument for generating the original spectrogram data;
specifically, comparing the slope of each point on the chromatographic curve with a threshold, and determining the reference point of the chromatographic peak based on the comparison result includes:
calculating the slope G of a point i on a chromatographic curve in real time during the sliding of a time windowi，
The calculated slope GiWith a predetermined first threshold value TthreComparing, and when the slope of two consecutive points is greater than the first threshold value TthreI.e. Gi>Tthre,Gi-1>TthreComparing ordinate values of the two points, determining a point having a smaller ordinate value as a peak start point reference point, an
The change of the positive or negative slope of the point after the peak start point reference point is analytically determined, when the slope of one point is negative and the slope of the point before the point is positive, i.e. Gi-1>0,Gi<0, comparing ordinate values of the two points, determining a point having a larger ordinate value as a temporary peak top point reference point, and
slope of point after the temporary peak-top reference pointWith a predetermined second threshold value T ″threComparing, when the slope of two consecutive points is less than the second threshold T ″threWhen is at Gi<T`thre,Gi-1<T`threAnd comparing the ordinate values of the two points, and determining the point with the smaller ordinate value as the peak terminal point reference point.
First threshold value T in step S221threAnd a second threshold value TthreThe method is determined based on calculation and analysis of no-load output signals of a chromatographic instrument for generating original spectrogram data, and the specific determination process comprises the following steps:
carrying out statistical analysis on the slope change of an output baseline of the chromatographic instrument in no-load, calculating and determining the variance of the slope change, and further determining the standard deviation of the slope change;
in the field of analytical instrumentation, it is generally accepted that the slope change of random noise and baseline drift follows a normal distribution and has a zero mean, and therefore, here, it is only necessary to find its variance based on the following expression (1),
in the expression (1), e2Representing variance, e standard deviation, fiRepresenting the difference between the samples and the mean, the number of samples n is typically greater than 100.
Then, based on the statistical characteristics of normal distribution (3e can ensure that the slope of 97.3% of the baseline falls within a zero slope interval), the standard deviation of the slope change of three times is taken as a first threshold, and the standard deviation of the slope change of minus three times is taken as a second threshold, namely the first threshold Tthre3e, the second threshold value T ″thre＝-3e。
After the step S221, performing a step S222, performing pattern recognition detection on the chromatographic curve near the position by using a gaussian wave as a matching wave based on the position of the chromatographic peak represented by the reference point on the chromatographic curve, and using a peak top in a recognition detection result as a peak top correction point;
specifically, in step S222, similarly to the prior art, first, a chromatographic curve near the position is used as a curve to be detected, a gaussian wave waveform is slid on the curve to be detected from a left end point to a right end point while calculating correlation coefficients of the two, a correlation coefficient group of chromatographic data of the curve to be detected with respect to the gaussian wave is obtained based on the calculation result, and a formula of a calculation process is expressed as,
in the expression (2), R represents a correlation coefficient, xi,yiChromatographic peak data and matched wave data respectively representing the currently calculated correlation,
the mean values of the two sets of data are shown separately.
And then, comparing and analyzing each correlation coefficient in the correlation array with a preset value, determining the position of a Gaussian peak position based on the correlation coefficient of which the coefficient value is greater than the preset value, and determining a point at the position on the curve to be detected as a peak top point correction point.
For example, the predetermined value is 0.8, and when the correlation coefficient is greater than 0.8, which indicates that the two are strongly correlated, the process of determining the position of the gaussian peak is the same as the prior art, and the detailed description thereof is omitted here.
After step S222, step S223 is performed to determine a point at which the vertical coordinate value is large as the peak top reference point by comparing the vertical coordinate values of the peak top correction point and the temporary peak top reference point.
In this embodiment, step S224 is finally performed to determine the chromatographic peak characterized by the peak start point reference point, the peak end point reference point, and the peak top point reference point as the detected chromatographic peak.
In this embodiment, in step S220, the time window is slid, and in each chromatographic peak detection process, steps S221 to S224 are repeated until all raw spectrogram data are processed.
According to the technical scheme, the liquid phase chromatogram spectrogram and the gas phase chromatogram spectrogram are actually and specifically detected based on pattern matching and in combination with the slope of a chromatogram curve, and the advantages of two modes (such as insensitivity of pattern recognition to noise, chromatographic peak width and amplitude change, good anti-interference performance, fault tolerance and robustness) are combined, so that the detection reliability is integrally improved. In the analysis method, spectrogram data is identified based on a time window sliding mode, detection analysis processing can be performed while spectrogram data is output, and therefore detection analysis results can be output more quickly.
In addition, the threshold value in the slope detection mode in the method is automatically calculated and determined based on the self signal of the instrument, and the peak characteristic point detected by the mode matching is only used as a reference point to correct the peak characteristic point obtained based on the slope detection, so that the defects of the two modes in the prior art are also overcome.
In addition, in order to facilitate a user to quickly understand and evaluate the performance of the detection and analysis method of the present application, in a specific application scenario, on the basis of the above embodiments, the technical solution of the present application further includes:
detecting and processing original spectrogram data by adopting a reference chromatographic spectrogram detection algorithm to obtain a reference detection result; comparing and analyzing the detection result with the reference detection result, generating a detection evaluation report, and displaying and outputting the detection evaluation report; the reference chromatogram detection algorithm refers to other detection and analysis methods which are different from the detection and analysis method and have the same function and purpose as the detection and analysis method.
The comparing and analyzing the detection result with the reference detection result comprises: matching the detected chromatographic peaks in the detection result and the reference detection result, determining the matched chromatographic peaks (the process is shown in fig. 3), and generating a detection evaluation report based on the proportion of the matched chromatographic peaks in the detected chromatographic peaks and the difference of the matched chromatographic peaks (the process is shown in fig. 4).
In other words, in the evaluation and analysis process of the present application, instead of comparing the results in a time series in a one-to-one manner, the results of the detected peaks are first arranged in a descending manner according to the peak heights, where the peak results include the start point, the end point, the peak height, the area, the retention time (the time corresponding to the peak top), the start point of the baseline, and the end point of the baseline, and then the results after the sorting are compared and matched (one peak result corresponds to one piece of data, and the whole piece of data will change along with the change of the peak height sequence).
As shown in fig. 3, in this implementation, the comparison matching is first considered in combination with the dual conditions of peak height and retention time. Firstly, comparing the peak heights, if the difference value of the peak heights of the comparison algorithm is in a certain range, comparing the retention time, if the difference value is also in a certain range, indicating that the corresponding peak energy of the comparison algorithm is correspondingly matched, and storing the comparison result for subsequently calculating the index values of the peak heights.
If any two conditions are not met, the current compared data are not matched, one index is fixed, and the other index is moved backwards until the data which can be matched are found or the other index is moved to the last index. It should be noted that the calculation results of different algorithms for the same peak of the same chromatogram data are not always the same, and due to the difference of the algorithms, the obtained results have a certain difference, but the phase difference value is not very large, so the difference value is within a certain range. Of course, similar peak heights may exist in the same chromatographic data. Therefore, when the results are compared, the application also can increase the comparison of retention time according to the corresponding situation to ensure the accuracy of the results.
After the matching process is completed, storing the matched peak result information in a one-to-one correspondence manner, then calculating a corresponding index according to each peak result information, entering a visual output stage of a detection evaluation report shown in fig. 4, specifically detecting an evaluation report in the stage, and including:
A. looking at approximate matching results
The number of matched and unmatched peaks is plotted for visual display, for example, a bar chart is used for display, and assuming that the number of peaks detected by method1 (the detection and analysis method of the present application) is n, the number of peaks detected by method2 (the reference chromatogram detection algorithm) is m, and the number of matched peaks is s, wherein s < min (n, m). A histogram contains three sections, the lower method1 being the number of unmatched peaks n-s, the middle being the number of two algorithmically matched peaks s, and the upper section being the number of unmatched peaks m-s of method 2. If the proportion of the middle part is larger, the number of the detected peaks of the two algorithms is larger, and the detection performance is closer.
B. In most cases, chromatographic data is noisy due to the conditions of the instrument and experiment, and different algorithms are not sensitive to noise. In matching, the small peak often appears not to match, but this situation may mislead the result obtained in a to some extent. Thus here is added a matching sum of the peak areas
Occupying the total area of all peaks
In percentage (b)
Bar diagram. If the two percentage results are very close and the values are both large, the two groups of result large peaks are considered to be matched, and the difference of the detection performances of the two algorithms is verified from another index.
C. Checking whether the matching result is abnormal or not
Specifically, for example, std (root mean square error) of the start-end point difference value of the matched upper peak is calculated. The std of the peak area difference, if std is floated within a certain range, the peak matched by the two comparison algorithms is considered to have no more prominent abnormality at the start point and the end point. Otherwise, the abnormal point is considered to be present, and then the abnormal data can be used for checking the abnormal condition of the detection algorithm.
Furthermore, it is easy to understand that, for different original data, the algorithm is contrasted and analyzed based on corresponding results, one data can correspond to one evaluation report, and a comprehensive evaluation is formed for the overall data result according to the reports to comprehensively evaluate the performance of the algorithm.
Fig. 5 is a schematic structural diagram of an electronic device according to an embodiment of the present application, and as shown in fig. 5, the electronic device 400 includes:
a memory 401 having an executable program stored thereon;
a processor 402 for executing the executable program in the memory 401 to implement the steps of the above method.
With respect to the electronic device 400 in the above embodiment, the specific manner of executing the program in the memory 401 by the processor 402 thereof has been described in detail in the embodiment related to the method, and will not be elaborated herein.
The above description is only a preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are included in the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (9)
1. A method for detecting and analyzing a chromatographic spectrum is characterized by comprising the following steps:
acquiring original spectrogram data to be processed;
for the original spectrogram data, sequentially sliding from the starting point of a chromatogram curve by adopting a time window with a preset size to perform detection processing of chromatographic peaks until the original spectrogram data is processed, detecting all chromatographic peaks in the curve, and generating a detection result based on the detected chromatographic peaks; wherein, in the detection processing process of each chromatographic peak, the method comprises the following steps:
comparing the slopes of the points on the chromatographic curve with a threshold value, and determining reference points of chromatographic peaks based on the comparison result, wherein the reference points comprise a peak starting point reference point, a temporary peak top point reference point and a peak end point reference point, and the threshold value is determined based on the calculation and analysis of the no-load output signal of the chromatographic instrument generating the original spectrogram data;
based on the position of the chromatographic peak represented by the reference point on the chromatographic curve, performing pattern recognition detection on the chromatographic curve near the position by taking a Gaussian wave as a matching wave, and taking the peak top in a recognition detection result as a peak top correction point;
correcting the temporary peak point reference point according to the peak point correction point, and taking the corrected point as a peak point reference point;
and determining the chromatographic peak characterized by the peak starting point reference point, the peak end point reference point and the peak top point reference point as the detected chromatographic peak.
2. The detection analysis method of claim 1, wherein comparing the slope of each point on the chromatographic curve with a threshold value and determining the reference point of the chromatographic peak based on the comparison comprises:
calculating the slope of a point on the chromatographic curve during the sliding of the time window in real time,
comparing the calculated slope with a predetermined first threshold, comparing ordinate values of two consecutive points when the slopes of the two points are greater than the first threshold, determining a point having a smaller ordinate value as the peak start point reference point, and
analyzing and judging the positive and negative changes of the slope of the point after the peak starting point reference point, when the slope of one point is negative and the slope of the point before the point is positive, comparing the vertical coordinate values of the two points, determining the point with the larger vertical coordinate value as the temporary peak top point reference point, and
and comparing the slope of the point behind the temporary peak top point reference point with a predetermined second threshold, when the slopes of two continuous points are smaller than the second threshold, comparing the ordinate values of the two points, and determining the point with the smaller ordinate value as the peak terminal point reference point.
3. The detection analysis method of claim 2, wherein the computational analysis of the no-load output signal of the chromatographic instrument generating the raw spectrogram data comprises:
carrying out statistical analysis on the slope change of an output baseline of the chromatographic instrument in no-load, calculating and determining the variance of the slope change, and further determining the standard deviation of the slope change;
three times the standard deviation of the slope change is taken as the first threshold, and minus three times the standard deviation of the slope change is taken as the second threshold.
4. The detection analysis method according to claim 1, wherein the temporary peak-top reference point is corrected according to the peak-top correction point, and the corrected point is used as the peak-top reference point, and specifically:
and comparing the longitudinal coordinate values of the peak top point correction point and the temporary peak top point reference point, and determining a point with a larger longitudinal coordinate value as a peak top point reference point.
5. The detection analysis method according to claim 1, wherein performing pattern recognition detection on the chromatogram curve near the position by using a gaussian wave as a matching wave, and using a peak top in a recognition detection result as a peak top correction point comprises:
taking the chromatographic curve near the position as a curve to be detected, sliding the waveform of the Gaussian wave on the curve to be detected from a left end point to a right end point, simultaneously calculating the correlation coefficient of the chromatographic curve and the Gaussian wave, and obtaining a correlation coefficient array of chromatographic data of the curve to be detected relative to the Gaussian wave based on the calculation result;
and comparing and analyzing each correlation coefficient in the correlation coefficient group with a preset value, determining the position of a Gaussian peak position based on the correlation coefficient of which the coefficient value is greater than the preset value, and determining a point at the position on the curve to be detected as the peak top point correction point.
6. The detection analysis method according to claim 1, wherein the generating a detection result based on the detected chromatographic peak comprises:
and integrating the detected chromatographic peaks, and calculating to determine the areas and the heights of the chromatographic peaks.
7. The detection analysis method according to claim 6, further comprising,
detecting and processing the original spectrogram data by adopting a reference chromatogram spectrogram detection algorithm to obtain a reference detection result;
and comparing and analyzing the detection result with the reference detection result, generating a detection evaluation report, and displaying and outputting the detection evaluation report.
8. The assay method of claim 7, wherein the comparing the assay result with the reference assay result comprises:
matching the chromatographic peaks detected in the detection result and the reference detection result, determining the matched chromatographic peaks, and generating a detection evaluation report based on the proportion of the matched chromatographic peaks in the detected chromatographic peaks and the difference of the matched chromatographic peaks.
9. An electronic device, comprising:
a memory having an executable program stored thereon;
a processor for executing the executable program in the memory to implement the steps of the method of any one of claims 1-9.
Priority Applications (7)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110830534.7A CN113567604B (en) | 2021-07-22 | 2021-07-22 | Detection and analysis method of chromatographic spectrogram and electronic equipment |
| US17/640,767 US20230152286A1 (en) | 2021-02-25 | 2022-02-18 | Chromatographic analysis system, chromatogram detection and analysis method and electronic device |
| EP22707004.2A EP4080514A4 (en) | 2021-02-25 | 2022-02-18 | Chromatographic analysis system, method for detecting and analyzing chromatogram, and electronic device |
| JP2022516419A JP7369279B2 (en) | 2021-02-25 | 2022-02-18 | Chromatography analysis system, chromatogram detection and analysis method, and electronic equipment |
| PCT/CN2022/076808 WO2022179444A1 (en) | 2021-02-25 | 2022-02-18 | Chromatographic analysis system, method for detecting and analyzing chromatogram, and electronic device |
| BE20225122A BE1029100B1 (en) | 2021-02-25 | 2022-02-23 | CHROMATOGRAPHIC ANALYSIS SYSTEM, DETECTION AND ANALYSIS METHOD FOR CHROMATOGRAMS AND ELECTRONIC DEVICE |
| CH00179/22A CH718313A2 (en) | 2021-02-25 | 2022-02-23 | Chromatographic analysis system, detection and analysis method for chromatograms. |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110830534.7A CN113567604B (en) | 2021-07-22 | 2021-07-22 | Detection and analysis method of chromatographic spectrogram and electronic equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113567604A true CN113567604A (en) | 2021-10-29 |
| CN113567604B CN113567604B (en) | 2022-09-30 |
Family
ID=78166393
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110830534.7A Active CN113567604B (en) | 2021-02-25 | 2021-07-22 | Detection and analysis method of chromatographic spectrogram and electronic equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113567604B (en) |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114062280A (en) * | 2021-11-03 | 2022-02-18 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | Transformer tap switch fault monitoring method, system and device |
| CN114137103A (en) * | 2021-11-03 | 2022-03-04 | 湖北省农业科学院农业质量标准与检测技术研究所 | Method for converting liquid chromatography diode array data into fingerprint |
| CN114157683A (en) * | 2022-02-10 | 2022-03-08 | 华谱科仪(北京)科技有限公司 | Online chromatograph data storage method and device and storage medium |
| CN114166988A (en) * | 2022-02-09 | 2022-03-11 | 华谱科仪(北京)科技有限公司 | Chromatograph cluster fault detection method and device |
| CN114218244A (en) * | 2022-02-23 | 2022-03-22 | 华谱科仪(北京)科技有限公司 | Online chromatograph database updating method, data identification method and device |
| CN114218993A (en) * | 2022-02-22 | 2022-03-22 | 华谱科仪(北京)科技有限公司 | Online chromatograph control method and device, storage medium and server |
| CN114280212A (en) * | 2022-01-10 | 2022-04-05 | 华谱科仪(北京)科技有限公司 | Chromatographic detection correction method, storage medium and electronic equipment |
| CN114371232A (en) * | 2021-12-22 | 2022-04-19 | 天津国科医工科技发展有限公司 | Chromatographic filtering method, device, medium and system based on Kalman filtering algorithm |
| WO2022179444A1 (en) * | 2021-02-25 | 2022-09-01 | 华谱科仪(大连)科技有限公司 | Chromatographic analysis system, method for detecting and analyzing chromatogram, and electronic device |
| CN115078617A (en) * | 2022-06-15 | 2022-09-20 | 思凡(上海)石油设备有限公司 | Method and system for identifying chromatographic peak retention time, terminal device and storage medium |
| CN115201395A (en) * | 2022-09-16 | 2022-10-18 | 华谱科仪(北京)科技有限公司 | Chromatograph supervision method and device based on big data technology |
| CN116387208A (en) * | 2023-06-02 | 2023-07-04 | 合肥喆塔科技有限公司 | Chamber matching analysis method, system, equipment and medium based on threshold control |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3714813A (en) * | 1971-01-12 | 1973-02-06 | Dow Chemical Co | Apparatus for stabilizing component elution times in gas chromatographs |
| US4633169A (en) * | 1983-07-08 | 1986-12-30 | Kontron Holding Ag | Blood particle analyzer |
| EP0395481A2 (en) * | 1989-04-25 | 1990-10-31 | Spectra-Physics, Inc. | Method and apparatus for estimation of parameters describing chromatographic peaks |
| US20030110000A1 (en) * | 2001-11-30 | 2003-06-12 | Quimby Bruce D. | Method of matching retention times among multiple chromatographic system |
| JP2004271422A (en) * | 2003-03-11 | 2004-09-30 | Shimadzu Corp | Data processor for chromatograph |
| CN1712955A (en) * | 2004-06-25 | 2005-12-28 | 中国科学院大连化学物理研究所 | Precisive measurement for parameter of chromatography spike and area of overlapped peak |
| US20140129169A1 (en) * | 2012-11-05 | 2014-05-08 | Shimadzu Corporation | Chromatogram data processing method and device |
| CN103808492A (en) * | 2014-02-13 | 2014-05-21 | 上海温光自动化技术有限公司 | Curve display method of OTDR collected data |
| US20160252484A1 (en) * | 2013-10-09 | 2016-09-01 | Spectrosense Ltd. | System and method for modified gas chromatographic data analysis |
| CN109993155A (en) * | 2019-04-23 | 2019-07-09 | 北京理工大学 | For the characteristic peak extracting method of low signal-to-noise ratio uv raman spectroscopy |
| JP2020091700A (en) * | 2018-12-06 | 2020-06-11 | 株式会社日立製作所 | Data extracting method and data analyzing system |
| CN111595992A (en) * | 2020-06-30 | 2020-08-28 | 浙江三青环保科技有限公司 | Rapid peak searching method for online gas chromatographic peak |
| CN112863615A (en) * | 2021-02-25 | 2021-05-28 | 华谱科仪(大连)科技有限公司 | Chromatographic analysis system based on double-database architecture |
-
2021
- 2021-07-22 CN CN202110830534.7A patent/CN113567604B/en active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3714813A (en) * | 1971-01-12 | 1973-02-06 | Dow Chemical Co | Apparatus for stabilizing component elution times in gas chromatographs |
| US4633169A (en) * | 1983-07-08 | 1986-12-30 | Kontron Holding Ag | Blood particle analyzer |
| EP0395481A2 (en) * | 1989-04-25 | 1990-10-31 | Spectra-Physics, Inc. | Method and apparatus for estimation of parameters describing chromatographic peaks |
| US20030110000A1 (en) * | 2001-11-30 | 2003-06-12 | Quimby Bruce D. | Method of matching retention times among multiple chromatographic system |
| JP2004271422A (en) * | 2003-03-11 | 2004-09-30 | Shimadzu Corp | Data processor for chromatograph |
| CN1712955A (en) * | 2004-06-25 | 2005-12-28 | 中国科学院大连化学物理研究所 | Precisive measurement for parameter of chromatography spike and area of overlapped peak |
| US20140129169A1 (en) * | 2012-11-05 | 2014-05-08 | Shimadzu Corporation | Chromatogram data processing method and device |
| US20160252484A1 (en) * | 2013-10-09 | 2016-09-01 | Spectrosense Ltd. | System and method for modified gas chromatographic data analysis |
| CN103808492A (en) * | 2014-02-13 | 2014-05-21 | 上海温光自动化技术有限公司 | Curve display method of OTDR collected data |
| JP2020091700A (en) * | 2018-12-06 | 2020-06-11 | 株式会社日立製作所 | Data extracting method and data analyzing system |
| CN109993155A (en) * | 2019-04-23 | 2019-07-09 | 北京理工大学 | For the characteristic peak extracting method of low signal-to-noise ratio uv raman spectroscopy |
| CN111595992A (en) * | 2020-06-30 | 2020-08-28 | 浙江三青环保科技有限公司 | Rapid peak searching method for online gas chromatographic peak |
| CN112863615A (en) * | 2021-02-25 | 2021-05-28 | 华谱科仪(大连)科技有限公司 | Chromatographic analysis system based on double-database architecture |
Non-Patent Citations (3)
| Title |
|---|
| PAUL G. STEVENSON 等: "Automated methods for the location of the boundaries of chromatographic peaks", 《JOURNAL OF CHROMATOGRAPHY A》 * |
| 朱凯汉 等: "关于谱峰起始值及求积的研讨", 《浙江师范大学学报(自然科学版)》 * |
| 陈淑珍 等: "一种色谱谱峰识别算法的实现", 《计算机应用与软件》 * |
Cited By (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2022179444A1 (en) * | 2021-02-25 | 2022-09-01 | 华谱科仪(大连)科技有限公司 | Chromatographic analysis system, method for detecting and analyzing chromatogram, and electronic device |
| CN114062280B (en) * | 2021-11-03 | 2023-10-27 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | Transformer tap switch fault monitoring method, system and device |
| CN114137103A (en) * | 2021-11-03 | 2022-03-04 | 湖北省农业科学院农业质量标准与检测技术研究所 | Method for converting liquid chromatography diode array data into fingerprint |
| CN114062280A (en) * | 2021-11-03 | 2022-02-18 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | Transformer tap switch fault monitoring method, system and device |
| CN114371232A (en) * | 2021-12-22 | 2022-04-19 | 天津国科医工科技发展有限公司 | Chromatographic filtering method, device, medium and system based on Kalman filtering algorithm |
| CN114371232B (en) * | 2021-12-22 | 2024-03-22 | 天津国科医工科技发展有限公司 | Chromatographic filtering method, device, medium and system based on Kalman filtering algorithm |
| CN114280212A (en) * | 2022-01-10 | 2022-04-05 | 华谱科仪(北京)科技有限公司 | Chromatographic detection correction method, storage medium and electronic equipment |
| CN114280212B (en) * | 2022-01-10 | 2022-06-03 | 华谱科仪(北京)科技有限公司 | Chromatographic detection correction method, storage medium and electronic equipment |
| CN114166988A (en) * | 2022-02-09 | 2022-03-11 | 华谱科仪(北京)科技有限公司 | Chromatograph cluster fault detection method and device |
| CN114157683A (en) * | 2022-02-10 | 2022-03-08 | 华谱科仪(北京)科技有限公司 | Online chromatograph data storage method and device and storage medium |
| CN114157683B (en) * | 2022-02-10 | 2022-05-13 | 华谱科仪(北京)科技有限公司 | Online chromatograph data storage method and device and storage medium |
| CN114218993A (en) * | 2022-02-22 | 2022-03-22 | 华谱科仪(北京)科技有限公司 | Online chromatograph control method and device, storage medium and server |
| CN114218244A (en) * | 2022-02-23 | 2022-03-22 | 华谱科仪(北京)科技有限公司 | Online chromatograph database updating method, data identification method and device |
| CN115078617A (en) * | 2022-06-15 | 2022-09-20 | 思凡(上海)石油设备有限公司 | Method and system for identifying chromatographic peak retention time, terminal device and storage medium |
| CN115078617B (en) * | 2022-06-15 | 2024-06-21 | 思凡(上海)石油设备有限公司 | Identification method, identification system, identification terminal device and identification storage medium for chromatographic peak retention time |
| CN115201395A (en) * | 2022-09-16 | 2022-10-18 | 华谱科仪(北京)科技有限公司 | Chromatograph supervision method and device based on big data technology |
| CN116387208A (en) * | 2023-06-02 | 2023-07-04 | 合肥喆塔科技有限公司 | Chamber matching analysis method, system, equipment and medium based on threshold control |
| CN116387208B (en) * | 2023-06-02 | 2023-08-18 | 合肥喆塔科技有限公司 | Chamber matching analysis method, system, equipment and medium based on threshold control |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113567604B (en) | 2022-09-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113567604B (en) | Detection and analysis method of chromatographic spectrogram and electronic equipment | |
| CN113567603B (en) | Detection and analysis method of chromatographic spectrogram and electronic equipment | |
| CN109239360B (en) | Reaction curve abnormity detection method and device | |
| JP6610678B2 (en) | Peak detection method and data processing apparatus | |
| US10198630B2 (en) | Peak detection method | |
| CN102445544B (en) | Method and system for increasing judgment accuracy of monoisotopic peaks | |
| CN109541401B (en) | Cable detection method and device and electronic equipment | |
| CN104305991B (en) | The method of detection noise and equipment from electrocardiosignal | |
| CN115982602B (en) | Photovoltaic transformer electrical fault detection method | |
| US20170356889A1 (en) | Three-dimensional spectral data processing device and processing method | |
| WO2016004687A1 (en) | Method for distinguishing initial time point of ultra-high-frequency partial discharge signal | |
| CN113076986B (en) | Photovoltaic fault arc characteristic selection method combining filtering type and packaging type evaluation strategies | |
| JP2017027076A (en) | Method and apparatus for detecting correctness of pitch period | |
| CN109425894A (en) | A kind of seismic anomaly road detection method and device | |
| CN111191671A (en) | Electrical appliance waveform detection method and system, electronic equipment and storage medium | |
| JPS6375573A (en) | Method of discriminating waveform | |
| CN102542284B (en) | Method for identifying spectrum | |
| CN110673125B (en) | Sound source positioning method, device, equipment and storage medium based on millimeter wave radar | |
| US20050061057A1 (en) | Interference detection systems and methods | |
| CN116660176A (en) | Fourier spectrum automatic baseline correction method, device and storage medium | |
| US20110130989A1 (en) | System and method for identifying a peripheral component interconnect express signal | |
| CN113723207B (en) | Acoustic emission signal mutation detection method based on histogram distance | |
| EP3561718A1 (en) | Method for identifying raman spectrogram, and electronic device | |
| CN115047400A (en) | Method and system for checking accuracy of three-phase electric energy meter, terminal equipment and medium | |
| CN110376436B (en) | Multi-scale noise power spectral line spectrum detection method |
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 | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20240105 Address after: 101407 101, Zone 4, 1st Floor, Building 1, No. 6, Yanqi South 4th Street, Yanqi Economic Development Zone, Huairou District, Beijing Patentee after: China spectrum Tech (Beijing) Technology Co.,Ltd. Address before: 116000 Room 102, 1f, front building, No. 909c, Huangpu Road, high tech Industrial Park, Dalian, Liaoning Patentee before: Huapu Keyi (Dalian) Technology Co.,Ltd. |
|
| TR01 | Transfer of patent right |