CN102135496A - Infrared spectrum quantitative analysis method and infrared spectrum quantitative analysis device based on multi-scale regression - Google Patents

Infrared spectrum quantitative analysis method and infrared spectrum quantitative analysis device based on multi-scale regression Download PDF

Info

Publication number
CN102135496A
CN102135496A CN2010106013539A CN201010601353A CN102135496A CN 102135496 A CN102135496 A CN 102135496A CN 2010106013539 A CN2010106013539 A CN 2010106013539A CN 201010601353 A CN201010601353 A CN 201010601353A CN 102135496 A CN102135496 A CN 102135496A
Authority
CN
China
Prior art keywords
infrared spectrum
quantitative analysis
infrared
scale
model
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.)
Pending
Application number
CN2010106013539A
Other languages
Chinese (zh)
Inventor
郝勇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
East China Jiaotong University
Original Assignee
East China Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by East China Jiaotong University filed Critical East China Jiaotong University
Priority to CN2010106013539A priority Critical patent/CN102135496A/en
Publication of CN102135496A publication Critical patent/CN102135496A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Spectrometry And Color Measurement (AREA)

Abstract

本发明涉及基于多尺度回归的红外光谱定量分析方法和装置,包括数据信号线相连的光谱仪、预处理器、小波分解与重构处理器、偏最小二乘回归模型集成器;所述的红外光谱包括中红外和近红外光谱,即波长范围为:780nm-50000nm。通过小波分解和重构变换,实现了多模型的构建,克服了单模型方法对光谱信号信息提取的困难;通过对不同子模型单独确定因子数,实现了有效信息的充分提取,提高了红外光谱分析模型的预测精度和稳健性。

The invention relates to an infrared spectrum quantitative analysis method and device based on multi-scale regression, including a spectrometer connected to a data signal line, a preprocessor, a wavelet decomposition and reconstruction processor, and a partial least squares regression model integrator; the infrared spectrum Including mid-infrared and near-infrared spectra, that is, the wavelength range is: 780nm-50000nm. Through wavelet decomposition and reconstruction transformation, the construction of multiple models is realized, which overcomes the difficulty of extracting spectral signal information by the single model method; by separately determining the number of factors for different sub-models, the effective information is fully extracted and the infrared spectrum is improved. Analyze the predictive accuracy and robustness of the model.

Description

基于多尺度回归的红外光谱定量分析方法和装置Method and device for quantitative analysis of infrared spectroscopy based on multi-scale regression

技术领域technical field

本发明涉及了一种红外光谱定量分析方法和装置,特别是基于多尺度回归的红外光谱定量分析方法和装置。The invention relates to an infrared spectrum quantitative analysis method and device, in particular to an infrared spectrum quantitative analysis method and device based on multi-scale regression.

背景技术Background technique

由于光谱仪得到的光谱信号中既含有有用信息,同时也叠加着其它各种随机误差(背景和噪声)。因此,在利用偏最小二乘回归方法进行定量分析时,很难得到预测精度较高的模型。针对提高偏最小二乘回归模型的预测精度这一问题,大量的研究工作已被开展,目前主要包括光谱预处理方法和变量筛选方法研究,这些方法已经成功地用于光谱背景校正、噪声消除、无信息变量的消除。然而在利用上述方法进行信息提取时,常会面临有用信息提取不充分或提取的信息中包含有噪声的问题。因此,需要发明一种不需要对光谱信息进行提取,而是通过将光谱包含的信息进行分类,然后对分类后的光谱信息分别进行建模预测,最后对各个不同类的光谱信息的模型结果进行集成,这样既可以避免光谱信息提取的不准确,又可以提高光谱模型对最终结果的预测准确率。Because the spectral signal obtained by the spectrometer contains both useful information and other random errors (background and noise). Therefore, it is difficult to obtain a model with high prediction accuracy when using the partial least squares regression method for quantitative analysis. To improve the prediction accuracy of the partial least squares regression model, a lot of research work has been carried out, mainly including spectral preprocessing methods and variable screening methods, these methods have been successfully used in spectral background correction, noise removal, Elimination of non-informative variables. However, when using the above methods for information extraction, we often face the problem of insufficient extraction of useful information or noise in the extracted information. Therefore, it is necessary to invent a method that does not need to extract the spectral information, but classifies the information contained in the spectrum, and then models and predicts the classified spectral information respectively, and finally performs model results for different types of spectral information. Integrating, this can not only avoid inaccurate extraction of spectral information, but also improve the prediction accuracy of the final result of the spectral model.

小波变换已经成功地应用于光谱信号的预处理中,包括数据压缩、平滑滤噪、基线校正、重叠信号解析以及分析图象处理等领域。与其它预处理方法不同,小波变换具有“时频优势”。通过小波变换处理,一条光谱信号可以被分成不同频率的几个子信号,对这些子信号进行重构后,重构光谱所包含的信息几乎与原始光谱完全相同。最后,通过对不同频率下的重构子光谱信号分别建模预测,并对预测结果进行集成。Wavelet transform has been successfully applied in the preprocessing of spectral signals, including data compression, smoothing and noise filtering, baseline correction, overlapping signal analysis and image processing. Unlike other preprocessing methods, wavelet transform has a "time-frequency advantage". Through wavelet transform processing, a spectral signal can be divided into several sub-signals of different frequencies. After these sub-signals are reconstructed, the information contained in the reconstructed spectrum is almost identical to the original spectrum. Finally, the reconstructed subspectral signals at different frequencies are modeled and predicted separately, and the predicted results are integrated.

在本专利中,发明了一种基于小波分解与重构结合偏最小二乘回归方法的红外光谱定量分析方法和装置。克服了单一模型方法对光谱信息提取的难点,在分析复杂的光谱数据时,利用本发明中的方法和装置可以直接利用各尺度下重构光谱建立多个偏最小二乘子模型,这些子模型因子数可以依据子光谱包含的信息量更加灵活有效地进行选择,达到充分提取各个尺度下的有用信息的目的。In this patent, an infrared spectrum quantitative analysis method and device based on wavelet decomposition and reconstruction combined with partial least squares regression method are invented. It overcomes the difficulty of extracting spectral information by a single model method. When analyzing complex spectral data, the method and device of the present invention can directly use the reconstructed spectra at each scale to establish multiple partial least squares sub-models. These sub-models The number of factors can be selected more flexibly and effectively according to the amount of information contained in the sub-spectrum, so as to achieve the purpose of fully extracting useful information at each scale.

发明内容Contents of the invention

在应用红外光谱分析技术对待测组分进行快速分析过程中,为了解决光谱信息提取的不准确的问题。本发明提供了一种基于多尺度回归的红外光谱定量分析方法和装置。该方法和装置充分利用的小波变换的多分辨分析的特点,实现了红外光谱信息的充分利用和合理分配。In order to solve the problem of inaccurate extraction of spectral information in the process of rapid analysis of components to be measured by infrared spectral analysis technology. The invention provides an infrared spectrum quantitative analysis method and device based on multi-scale regression. The method and device make full use of the characteristics of multi-resolution analysis of wavelet transform, and realize the full use and reasonable distribution of infrared spectrum information.

实现上述技术方案的装置包括:经数据信号线相连的光谱仪、预处理器、小波分解与重构处理器、偏最小二乘回归模型集成器。The device for realizing the technical solution includes: a spectrometer connected via a data signal line, a preprocessor, a wavelet decomposition and reconstruction processor, and a partial least square regression model integrator.

所述的红外光谱包括中红外和近红外光谱,即波长范围为:780 nm - 50000 nm。The infrared spectrum includes mid-infrared and near-infrared spectra, that is, the wavelength range is: 780 nm-50000 nm.

所述的预处理器采取中心化和矢量归一化对光谱仪采集的原始信号进行处理。The preprocessor adopts centralization and vector normalization to process the original signal collected by the spectrometer.

所述的小波分解与重构处理器,具体包含的处理步骤如下:首先,在小波分解过程中,需要对分解尺度和小波基两个参数进行设定;其次对分解后各尺度下的近似或细节光谱成分进行重构,这样就把原始的光谱信号数据阵变换为与原始光谱信号数据阵维数相同的几个子光谱信号数据阵。The wavelet decomposition and reconstruction processor specifically includes the following processing steps: first, in the wavelet decomposition process, two parameters, the decomposition scale and the wavelet base, need to be set; secondly, the approximate or The detailed spectral components are reconstructed, so that the original spectral signal data array is transformed into several sub-spectral signal data arrays with the same dimensions as the original spectral signal data array.

所述的偏最小二乘回归模型集成器,其具体操作如下:分别对经小波分解与重构处理器处理后得到的子光谱信号数据阵建立相应的偏最小二乘回归子模型(在每一个子模型中,主成分数可以选择不同的数值),利用已建立的子模型分别对测试集样品指标进行预测,分别得到相应的预测结果;最后对每一个子模型的预测结果进行加权,得到最终的样品指标的预测值。通过比较不同分解尺度和小波基下,模型的预测均方根误差(Root Mean Square Error of Prediction, RMSEP)值,确定合适的分解尺度、小波基和模型参数并保存,用于后续新样品红外光谱的预测分析。The specific operation of the partial least squares regression model integrator is as follows: respectively establish corresponding partial least squares regression submodels for the subspectral signal data arrays obtained after processing by the wavelet decomposition and reconstruction processor (in each In the sub-model, the principal component scores can choose different values), use the established sub-models to predict the test set sample indicators respectively, and obtain the corresponding prediction results; finally, weight the prediction results of each sub-model to obtain the final The predicted value of the sample index of . By comparing the root mean square error (Root Mean Square Error of Prediction, RMSEP) value of the model under different decomposition scales and wavelet bases, determine the appropriate decomposition scale, wavelet base and model parameters and save them for subsequent infrared spectra of new samples predictive analysis.

由于本发明采用以上的技术方案,得到以下效果:Because the present invention adopts above technical scheme, obtain following effect:

通过小波分解和重构变换,实现了多模型的构建,克服了单模型方法对光谱信号信息提取的困难;通过对不同子模型单独确定因子数,实现了有效信息的充分提取,提高了红外光谱分析模型的预测精度和稳健性。Through wavelet decomposition and reconstruction transformation, the construction of multiple models is realized, which overcomes the difficulty of extracting spectral signal information by the single model method; by separately determining the number of factors for different sub-models, the effective information is fully extracted and the infrared spectrum is improved. Analyze the predictive accuracy and robustness of the model.

附图说明Description of drawings

图1 多尺度回归的红外光谱定量分析方法和装置示意图;Figure 1 Schematic diagram of the infrared spectroscopy quantitative analysis method and device for multi-scale regression;

图2 近红外光谱图;Figure 2 near-infrared spectrum;

图3 多尺度偏最小二乘回归方法操作示意图;Figure 3 Schematic diagram of the operation of the multi-scale partial least squares regression method;

图4 分解尺度对多尺度偏最小二乘回归方法的影响Figure 4. Effect of decomposition scale on multi-scale partial least squares regression method

图5 小波分解与重构处理器操作示意图;Fig. 5 Schematic diagram of wavelet decomposition and reconstruction processor operation;

图6 红外光谱模型的预测值和参考值的相关图。Fig. 6. The correlation diagram of the predicted value and the reference value of the infrared spectrum model.

具体实施方式Detailed ways

具体实施方式结合下面实施实例进行说明。以梨的近红外光谱为例,对梨内部的糖度指标进行多尺度回归模型的构建。The specific embodiment will be described in conjunction with the following implementation examples. Taking the near-infrared spectrum of pears as an example, a multi-scale regression model was constructed for the sugar content index inside pears.

图1为多尺度回归的红外光谱定量分析方法和装置示意图,图2为近红外光谱图,光谱范围为750 -1800 nm,每条光谱包括1051个数据点。将所有样品按照2:1的比例划分为校正集和测试集。Figure 1 is a schematic diagram of the multi-scale regression infrared spectrum quantitative analysis method and device, and Figure 2 is a near-infrared spectrum diagram with a spectral range of 750-1800 nm, and each spectrum includes 1051 data points. All samples were divided into calibration set and test set according to the ratio of 2:1.

将所有样品进行小波分解,分解时我们这里选择db4小波基,分解尺度从1到20,图3所示为多尺度偏最小二乘回归方法操作示意图,从图中可以看出光谱经小波分解重构后,得到一系列与原始光谱维数相同的子光谱信号矩阵,所有的子光谱图谱差异较大,每一个子光谱矩阵包含的有用信息也存在差异。因此建模时,原始光谱只能选择一个主成分数,而多尺度偏最小二乘回归方法可以对子光谱矩阵选取不同的主成分数,使得信息提取更加灵活充分。All samples are decomposed by wavelet. When decomposing, we choose db4 wavelet base here, and the decomposition scale is from 1 to 20. Figure 3 shows the operation diagram of the multi-scale partial least squares regression method. It can be seen from the figure that the spectrum is reconstructed by wavelet decomposition. After construction, a series of sub-spectral signal matrices with the same dimensions as the original spectrum are obtained. All sub-spectral spectra are quite different, and the useful information contained in each sub-spectral matrix is also different. Therefore, when modeling, only one principal component number can be selected for the original spectrum, while the multi-scale partial least squares regression method can select different principal component numbers for the subspectral matrix, making information extraction more flexible and sufficient.

图4所示为分解尺度对多尺度偏最小二乘回归方法的影响,图中虚线表示采用传统偏最小二乘回归方法的预测均方根误差,实线表示测试集的预测均方根误差随小波分解尺度的变化曲线,从图中可知,当采用5尺度小波分解重构时,测试集的预测均方根误差最小。从而可以确定小波分解与重构处理器的分解尺度和小波基。Figure 4 shows the influence of the decomposition scale on the multi-scale partial least squares regression method. The dotted line in the figure represents the predicted root mean square error of the traditional partial least squares regression method, and the solid line represents the predicted root mean square error of the test set with The change curve of the wavelet decomposition scale, as can be seen from the figure, when the 5-scale wavelet decomposition is used for reconstruction, the prediction root mean square error of the test set is the smallest. Therefore, the decomposition scale and wavelet basis of the wavelet decomposition and reconstruction processor can be determined.

图5所示为小波分解与重构处理器操作示意图。对图中各个子光谱信号矩阵分别建立校正模型,相应的主成分数分别为7、4、5、3、2、4,而不经小波分解与重构处理器进行处理时,模型的主成分数为4,可见采用多尺度偏最小二乘回归分析方法更利于光谱信息的提取。Figure 5 is a schematic diagram of the operation of the wavelet decomposition and reconstruction processor. Establish correction models for each sub-spectral signal matrix in the figure, and the corresponding principal component numbers are 7, 4, 5, 3, 2, 4 respectively. When not processed by wavelet decomposition and reconstruction processor, the principal components of the model The number is 4, it can be seen that the multi-scale partial least squares regression analysis method is more conducive to the extraction of spectral information.

图6所示为红外光谱多尺度偏最小二乘回归模型预测值和参考值的相关图。从图中可知,得到了较好的预测结果。Figure 6 shows the correlation diagram between the predicted value and the reference value of the infrared spectrum multi-scale partial least squares regression model. It can be seen from the figure that better prediction results are obtained.

本发明的效果是:通过采用多尺度偏最小二乘回归方法对红外光谱进行定量分析,有助于提取不同频率下光谱信号的有用信息,得到更加准确、稳定的预测结果。因此,该发明的方法和装置有望成为一种十分有应用前景的红外光谱分析方法。The effect of the present invention is: by using the multi-scale partial least squares regression method to quantitatively analyze the infrared spectrum, it is helpful to extract useful information of spectral signals at different frequencies, and obtain more accurate and stable prediction results. Therefore, the method and device of the invention are expected to become a very promising infrared spectroscopic analysis method.

Claims (1)

1.基于多尺度回归的红外光谱定量分析方法和装置,其特征在于:包括数据信号线相连的光谱仪、预处理器、小波分解与重构处理器、偏最小二乘回归模型集成器;1. The infrared spectrum quantitative analysis method and device based on multi-scale regression, characterized in that: comprising a spectrometer connected to a data signal line, a preprocessor, a wavelet decomposition and reconstruction processor, and a partial least squares regression model integrator; 所述的红外光谱包括中红外和近红外光谱,即波长范围为:780 nm - 50000 nm;The infrared spectrum includes mid-infrared and near-infrared spectra, that is, the wavelength range is: 780 nm-50000 nm; 所述的预处理器采取中心化和矢量归一化对光谱仪采集的原始信号进行处理。The preprocessor adopts centralization and vector normalization to process the original signal collected by the spectrometer.
CN2010106013539A 2010-12-23 2010-12-23 Infrared spectrum quantitative analysis method and infrared spectrum quantitative analysis device based on multi-scale regression Pending CN102135496A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010106013539A CN102135496A (en) 2010-12-23 2010-12-23 Infrared spectrum quantitative analysis method and infrared spectrum quantitative analysis device based on multi-scale regression

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010106013539A CN102135496A (en) 2010-12-23 2010-12-23 Infrared spectrum quantitative analysis method and infrared spectrum quantitative analysis device based on multi-scale regression

Publications (1)

Publication Number Publication Date
CN102135496A true CN102135496A (en) 2011-07-27

Family

ID=44295343

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010106013539A Pending CN102135496A (en) 2010-12-23 2010-12-23 Infrared spectrum quantitative analysis method and infrared spectrum quantitative analysis device based on multi-scale regression

Country Status (1)

Country Link
CN (1) CN102135496A (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102435556A (en) * 2011-09-20 2012-05-02 湖南大学 Accurate spectrum quantitative analysis method used for complex heterogeneous mixture system
CN103884670A (en) * 2014-03-13 2014-06-25 西安交通大学 Smoke component quantitative analysis method based on near infrared spectrum
CN104574341A (en) * 2013-10-11 2015-04-29 中国林业科学研究院资源信息研究所 Method and device for determining sugar degree of fruit
CN105842190A (en) * 2016-03-17 2016-08-10 浙江中烟工业有限责任公司 Near-infrared model transfer method based on spectral regression
CN109492707A (en) * 2018-11-28 2019-03-19 武汉轻工大学 Construction method, device, equipment and the storage medium of spectrum analysis model
CN114112978A (en) * 2021-12-13 2022-03-01 大连理工大学 A Wavelet Function Calibration Modeling Method for In-situ Infrared Spectroscopy to Detect Solution Concentration in Crystallization Process

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101806729A (en) * 2010-03-31 2010-08-18 中国人民解放军总后勤部油料研究所 In-use lubricating oil quality rapid testing method
CN101915744A (en) * 2010-07-05 2010-12-15 北京航空航天大学 Near-infrared spectroscopy non-destructive testing method and device for substance composition content

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101806729A (en) * 2010-03-31 2010-08-18 中国人民解放军总后勤部油料研究所 In-use lubricating oil quality rapid testing method
CN101915744A (en) * 2010-07-05 2010-12-15 北京航空航天大学 Near-infrared spectroscopy non-destructive testing method and device for substance composition content

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《化工时刊》 20080202 殷俊等 小波变换_偏最小二乘法用于苯胺_苯酚和苯甲酸的同时测定 第43-44、52页 1 第22卷, 第2期 *
褚小立等: "小波变换结合多维偏最小二乘方法用于近红外光谱定量分析", 《分析化学研究简报》, vol. 34, 30 September 2006 (2006-09-30), pages 175 - 178 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102435556A (en) * 2011-09-20 2012-05-02 湖南大学 Accurate spectrum quantitative analysis method used for complex heterogeneous mixture system
CN102435556B (en) * 2011-09-20 2013-09-04 湖南大学 Accurate spectrum quantitative analysis method used for complex heterogeneous mixture system
CN104574341A (en) * 2013-10-11 2015-04-29 中国林业科学研究院资源信息研究所 Method and device for determining sugar degree of fruit
CN104574341B (en) * 2013-10-11 2017-09-05 中国林业科学研究院资源信息研究所 A kind of determination method and apparatus of sugar degree
CN103884670A (en) * 2014-03-13 2014-06-25 西安交通大学 Smoke component quantitative analysis method based on near infrared spectrum
CN103884670B (en) * 2014-03-13 2016-01-20 西安交通大学 Based on the smoke components quantitative analysis method of near infrared spectrum
CN105842190A (en) * 2016-03-17 2016-08-10 浙江中烟工业有限责任公司 Near-infrared model transfer method based on spectral regression
CN105842190B (en) * 2016-03-17 2018-12-11 浙江中烟工业有限责任公司 A kind of method for transferring near infrared model returned based on spectrum
CN109492707A (en) * 2018-11-28 2019-03-19 武汉轻工大学 Construction method, device, equipment and the storage medium of spectrum analysis model
CN109492707B (en) * 2018-11-28 2020-10-23 武汉轻工大学 Method, device and equipment for constructing spectral analysis model and storage medium
CN114112978A (en) * 2021-12-13 2022-03-01 大连理工大学 A Wavelet Function Calibration Modeling Method for In-situ Infrared Spectroscopy to Detect Solution Concentration in Crystallization Process

Similar Documents

Publication Publication Date Title
CN101915744B (en) Near infrared spectrum nondestructive testing method and device for material component content
CN103234922A (en) Rapid soil organic matter detection method based on large sample soil visible-near infrared spectrum classification
CN101520412A (en) Near infrared spectrum analyzing method based on isolated component analysis and genetic neural network
CN113008805B (en) Radix angelicae decoction piece quality prediction method based on hyperspectral imaging depth analysis
CN104062258B (en) Method for rapid determination of soluble solids in compound ass-hide glue pulp by near infrared spectroscopy
CN104062257A (en) Method for determining total flavone content of solution based on near infrared spectroscopy
CN107247033B (en) The method of identifying the maturity of Huanghua pear based on the fast decay elimination algorithm and PLSDA
CN103854305A (en) Module transfer method based on multiscale modeling
CN102135496A (en) Infrared spectrum quantitative analysis method and infrared spectrum quantitative analysis device based on multi-scale regression
CN107478580A (en) Heavy metal content in soil evaluation method and device based on high-spectrum remote-sensing
CN101738383A (en) Mid-infrared spectrum-based method for quickly detecting contents of heavy metal elements in soil
CN109060771B (en) A Consensus Model Construction Method Based on Different Spectral Feature Sets
CN105158200A (en) Modeling method capable of improving accuracy of qualitative near-infrared spectroscopic analysis
CN102072767A (en) Wavelength similarity consensus regression-based infrared spectrum quantitative analysis method and device
CN106596456A (en) Changeable moving window based selection method of near infrared spectral region of solution
CN105486655A (en) Rapid detection method for organic matters in soil based on infrared spectroscopic intelligent identification model
WO2023123329A1 (en) Method and system for extracting net signal in near-infrared spectrum
CN102937575A (en) Watermelon sugar degree rapid modeling method based on secondary spectrum recombination
CN102128805A (en) Method and device for near infrared spectrum wavelength selection and quick quantitative analysis of fruit
CN106153561A (en) The many metal ion inspections of uv-vis spectra based on wavelength screening
CN117312968A (en) Method for predicting organic matter content of saline-alkali farmland soil
CN104897709A (en) Agricultural product element quantitative detection model building method based on X-ray fluorescence analysis
CN107632010A (en) A kind of quantitative approach of combination LIBS to steel samples
CN109001181B (en) A rapid identification method of edible oil species based on Raman spectroscopy canonical correlation analysis
CN114611582A (en) A method and system for analyzing substance concentration based on near-infrared spectroscopy

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20110727