CN110542658A - A classification method of tobacco non-smoke substances based on hyperspectral imaging technology - Google Patents

A classification method of tobacco non-smoke substances based on hyperspectral imaging technology Download PDF

Info

Publication number
CN110542658A
CN110542658A CN201910850623.0A CN201910850623A CN110542658A CN 110542658 A CN110542658 A CN 110542658A CN 201910850623 A CN201910850623 A CN 201910850623A CN 110542658 A CN110542658 A CN 110542658A
Authority
CN
China
Prior art keywords
tobacco
hyperspectral
sample
spectral
sundries
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
CN201910850623.0A
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.)
China Tobacco Fujian Industrial Co Ltd
Hefei Institutes of Physical Science of CAS
Zhengzhou Tobacco Research Institute of CNTC
Original Assignee
China Tobacco Fujian Industrial Co Ltd
Hefei Institutes of Physical Science of CAS
Zhengzhou Tobacco Research Institute of CNTC
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 China Tobacco Fujian Industrial Co Ltd, Hefei Institutes of Physical Science of CAS, Zhengzhou Tobacco Research Institute of CNTC filed Critical China Tobacco Fujian Industrial Co Ltd
Priority to CN201910850623.0A priority Critical patent/CN110542658A/en
Publication of CN110542658A publication Critical patent/CN110542658A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands

Landscapes

  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (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)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

一种基于高光谱成像技术的烟草非烟物质分类方法,其特征在于:是利用短波成像高光谱技术先对烟叶和杂物进行分类,建立包括不同物质的光谱库,然后对待测样品采集图像数据,利用光谱库中的参考光谱匹配所测样品,并对其进行判断,进而完成烟叶和杂物的有效分类识别。本发明结合了光谱技术和二维图像成像技术。其与现有技术相比,具有如下显著的进步:1、本发明利用短波成像高光谱技术对烟叶和杂物进行分类,建立包含有不同物质的光谱库。2、本发明实验过程不使用有毒有害化学品,简便、快捷、对样品无破坏性、对环境无污染。3、本发明具有操作简便、快速、准确、成本低、效率高的优点。

A method for classifying tobacco non-smoke substances based on hyperspectral imaging technology, which is characterized in that: using short-wave imaging hyperspectral technology to first classify tobacco leaves and sundries, establish a spectral library including different substances, and then collect image data from samples to be tested , use the reference spectrum in the spectral library to match the measured sample, and judge it, and then complete the effective classification and identification of tobacco leaves and sundries. The invention combines spectral technology and two-dimensional image imaging technology. Compared with the prior art, it has the following significant improvements: 1. The present invention uses the short-wave imaging hyperspectral technology to classify tobacco leaves and sundries, and establishes a spectral library containing different substances. 2. The experimental process of the present invention does not use toxic and harmful chemicals, which is simple, fast, non-destructive to the sample, and non-polluting to the environment. 3. The present invention has the advantages of simple operation, rapidity, accuracy, low cost and high efficiency.

Description

一种基于高光谱成像技术的烟草非烟物质分类方法A classification method of tobacco non-smoke substances based on hyperspectral imaging technology

技术领域technical field

本发明属于高光谱农业生产应用领域,尤其涉及一种基于高光谱成像技术的烟草非烟物质分类方法。The invention belongs to the field of hyperspectral agricultural production and application, in particular to a method for classifying tobacco non-smoke substances based on hyperspectral imaging technology.

背景技术Background technique

高光谱成像(Hyperspectral Image)是集探测器技术、精密光学机械、微弱信号检测、计算机技术、信息处理技术于一体的综合性技术。高光谱成像同时探测目标的二维几何空间与一维光谱信息、获取高光谱分辨率的连续、窄波段的图像数据,是一种将成像技术与光谱技术相结合的多维信息获取技术。因此,利用高光谱影像技术可以对物体进行精细反映。在烟叶行业中,烟叶原料的质量和纯净度直接关系到卷烟产品的质量。烟叶若在生产、收购、运输、加工过程中混入非烟叶杂物(一类杂物:金属、羽毛、塑料等;二类杂物:纸、石头、麻绳、玻璃等;三类杂物:非烟叶叶子、种子、竹签等),将对卷烟工业加工的卷烟成品造成极大的质量隐患。Hyperspectral imaging is a comprehensive technology that integrates detector technology, precision optical machinery, weak signal detection, computer technology, and information processing technology. Hyperspectral imaging simultaneously detects the two-dimensional geometric space and one-dimensional spectral information of the target, and obtains continuous and narrow-band image data with high spectral resolution. It is a multi-dimensional information acquisition technology that combines imaging technology with spectral technology. Therefore, the use of hyperspectral imaging technology can accurately reflect the object. In the tobacco industry, the quality and purity of tobacco raw materials are directly related to the quality of cigarette products. If tobacco leaves are mixed with non-tobacco sundries during production, acquisition, transportation and processing (Class I sundries: metal, feather, plastic, etc.; Class II sundries: paper, stone, hemp rope, glass, etc.; Class III sundries: Non-tobacco leaves, seeds, bamboo sticks, etc.) will cause great hidden dangers to the quality of cigarette products processed by the cigarette industry.

目前卷烟生产加工线上采用的除杂方式主要有:风力除杂、光电除杂、磁力除杂和人工挑选除杂,其中风力除杂、光电除杂、磁力除杂均是针对杂物的特定属性差异进行有针对性的识别并剔除,如比重的差异、颜色的差异和磁性的差异,因此一种除杂方式只能识别和剔除某一类的杂物。人工除杂可以识别绝大部分非烟物质杂物,但人眼容易疲劳且工作效率较低。At present, the methods of impurity removal used in cigarette production and processing lines mainly include: wind impurity removal, photoelectric impurity removal, magnetic impurity removal and manual selection of impurity removal. Attribute differences are identified and eliminated in a targeted manner, such as differences in specific gravity, color differences and magnetic differences. Therefore, a method of removing impurities can only identify and remove a certain type of debris. Manual removal of impurities can identify most of the non-smokers, but the human eyes are easily fatigued and work efficiency is low.

烟叶独特的化学组分与物理检测特性,导致其具有独特的诊断性特征吸收谱带,这些特征谱带中具有较稳定的波长位置和独特波形。从烟叶化学角度分析,烟叶内主要的化学成分为碳水、含氮化合物、有机酸和矿物质。其中,碳水的吸收特征表现在1000-2500nm,氮化合物的吸收特征表现在1500nm-1750nm,作为烟叶的重要化学成分,二者在可见光波段(400-1000nm)均无体现。即烟叶光谱在碳水与氮化合物附近会形成较强的吸收峰,这些吸收峰可以作为判别烟叶的重要依据。利用光谱信息技术,非烟叶杂物与烟叶在相同光学环境下形成不同特定光谱,识别出不同物质,完成烟叶与杂物的分类。The unique chemical composition and physical detection characteristics of tobacco leaves lead to unique diagnostic characteristic absorption bands, which have relatively stable wavelength positions and unique waveforms. From the chemical point of view of tobacco leaves, the main chemical components in tobacco leaves are carbohydrates, nitrogen compounds, organic acids and minerals. Among them, the absorption characteristics of carbon water are at 1000-2500nm, and the absorption characteristics of nitrogen compounds are at 1500nm-1750nm. As important chemical components of tobacco leaves, the two are not reflected in the visible light band (400-1000nm). That is, the spectrum of tobacco leaves will form strong absorption peaks near carbon, water and nitrogen compounds, and these absorption peaks can be used as an important basis for distinguishing tobacco leaves. Using spectral information technology, non-tobacco leaf debris and tobacco leaves form different specific spectra in the same optical environment, identify different substances, and complete the classification of tobacco leaves and debris.

中国专利(200910059486.5)公开了一种基于近红外光谱仪的烟丝参配比例的快速测定方法。与本专利提出方法相比,二者主要有两点不同:其一,二者的适用范围不同。中国专利(200910059486.5)主要用于烟丝参配比例测定,本文方法主要用于烟叶中杂物的识别;其二,由于所使用实验仪器的不同(中国专利(200910059486.5)使用近红外光谱仪,本方法使用高光谱成像仪)本方法操作更为简便快捷,无需进行样品粉碎、筛选等操作,可进行在线标注。Chinese Patent (200910059486.5) discloses a method for rapid determination of the proportion of tobacco cut tobacco based on a near-infrared spectrometer. Compared with the method proposed in this patent, there are two main differences between the two: First, the scope of application of the two is different. The Chinese patent (200910059486.5) is mainly used for the determination of the proportion of shredded tobacco, and the method in this paper is mainly used for the identification of impurities in tobacco leaves; Hyperspectral imager) This method is more convenient and quick to operate, without the need for sample crushing, screening and other operations, and can be marked online.

中国专利(201410491816.9)主要公开了一种基于光谱成像技术的烟丝组分识别方法。相较于本专利所提方法,有以下几点区别。1、本文所提方法使用的仪器为高光谱成像仪,分辨率为12nm,辨识精度更高。2、该专利所使用光谱相机由带通滤光片和带有定焦镜头的CCD相机构成,其带通滤光片所定制波段仅针对烟丝进行成像和图像处理分析,无法对其他物质进行分析和辨识。3、该专利所提方法仅用于测定烟丝组成成分(叶丝、梗丝、再造烟叶丝、膨胀叶丝),本文所提方法用于烟叶与非烟物质的识别,二者在功能上不一致。Chinese patent (201410491816.9) mainly discloses a method for identifying components of cut tobacco based on spectral imaging technology. Compared with the method proposed in this patent, there are the following differences. 1. The instrument used in the method proposed in this paper is a hyperspectral imager with a resolution of 12 nm and a higher identification accuracy. 2. The spectroscopic camera used in this patent is composed of a band-pass filter and a CCD camera with a fixed-focus lens. The customized wavelength band of the band-pass filter is only for imaging and image processing analysis of cut tobacco, and cannot be analyzed for other substances. and identification. 3. The method proposed in this patent is only used to determine the components of cut tobacco (cut leaf, cut stem, reconstituted tobacco leaf, and expanded cut leaf), and the method proposed in this paper is used for the identification of tobacco leaves and non-tobacco substances, and the two are inconsistent in function. .

综上,利用光谱技术测定烟丝含量、高光谱定量测定烟草中化学物质含量在行业内有一定的研究和应用。而利用光谱信息技术,根据非烟叶杂物与烟叶在相同光学环境下形成不同特定光谱,识别出不同物质,完成烟叶与非烟物质的分类尚属空白。In summary, the use of spectroscopic techniques to determine the content of cut tobacco and the quantitative determination of chemical substances in tobacco by hyperspectral have certain research and applications in the industry. However, using spectral information technology, according to the formation of different specific spectra of non-tobacco leaf debris and tobacco leaves in the same optical environment, different substances are identified, and the classification of tobacco leaves and non-tobacco substances is still blank.

发明内容SUMMARY OF THE INVENTION

本发明的目的正是基于上述现有技术状况而提供的一种基于高光谱成像技术的烟草非烟物质分类方法,能够快速且准确区分烟叶与杂物,杜绝非烟杂物进入片烟成品当中。The purpose of the present invention is to provide a method for classifying tobacco non-smoke substances based on hyperspectral imaging technology based on the above state of the art, which can quickly and accurately distinguish tobacco leaves and sundries, and prevent non-smoke sundries from entering the finished tobacco products. .

本发明的目的是通过以下技术方案来实现的:The purpose of this invention is to realize through the following technical solutions:

一种基于高光谱成像的烟叶杂物分类识别的方法,是利用短波高光谱成像技术先对烟叶和杂物进行分类,建立包括不同物质的光谱库,然后对待测样品采集高光谱图像数据,利用光谱库中的参考光谱匹配待测样品,并对其进行判断,进而完成烟叶和杂物的有效分类识别。A method for classifying and identifying tobacco leaf debris based on hyperspectral imaging is to use short-wave hyperspectral imaging technology to first classify tobacco leaves and debris, establish a spectral library including different substances, and then collect hyperspectral image data of the sample to be tested, and use The reference spectrum in the spectral library matches the sample to be tested and judges it, thereby completing the effective classification and identification of tobacco leaves and sundries.

具体步骤如下:Specific steps are as follows:

1)样本采集:获取包括纯净烟叶在内,一类杂物、二类杂物、三类杂物的样本;1) Sample collection: Obtain samples of Class I, Class II and Class III sundries including pure tobacco leaves;

样本表面保持干燥清洁,无其他附着物。The surface of the sample is kept dry and clean with no other attachments.

所述一类杂物包括金属、羽毛、塑料;二类杂物包括纸、石头、麻绳、玻璃;三类杂物包括非烟叶叶子、种子、竹签。The first-class sundries include metals, feathers, and plastics; the second-class sundries include paper, stone, hemp rope, and glass; and the third-class sundries include non-tobacco leaves, seeds, and bamboo sticks.

2)样品制备与高光谱成像并进行黑白帧校正;2) Sample preparation and hyperspectral imaging with black and white frame correction;

使用卤钨灯作为照明光源,对获取的烟叶和杂物样本进行高光谱图像采集,获取样本的高光谱图像。为减小噪声影响,对高光谱图像进行黑白帧校正。黑白帧校正公式如下:Using a tungsten halogen lamp as the illumination light source, hyperspectral image acquisition was performed on the obtained tobacco leaf and debris samples to obtain hyperspectral images of the samples. To reduce the effect of noise, black and white frame correction is performed on the hyperspectral image. The black and white frame correction formula is as follows:

式中:R-校正后高光谱图像;I-原始高光谱图像;B-关闭相机镜头采集的全黑图像;W-扫描白色校正板得到的全白图像。In the formula: R-corrected hyperspectral image; I-original hyperspectral image; B-completely black image collected by turning off the camera lens; W-completely white image obtained by scanning the white correction plate.

3)高光谱图像预处理及特征图像的获取;3) Hyperspectral image preprocessing and acquisition of characteristic images;

为提高数据信噪比,对高光谱图像数据进行预处理,预处理方法包括但不仅限于:In order to improve the signal-to-noise ratio of the data, the hyperspectral image data is preprocessed, and the preprocessing methods include but are not limited to:

通过使用Savitzky-Golay卷积平滑滤波算法,去除基线漂移和倾斜,去除噪声,提高光谱曲线平滑度;By using the Savitzky-Golay convolution smoothing filtering algorithm, the baseline drift and tilt are removed, noise is removed, and the smoothness of the spectral curve is improved;

再通过多元散射校正(MSC)消减物体表面散射效应,增强相同物质间光谱吸收信息。Then, the scattering effect on the surface of the object is reduced by multivariate scattering correction (MSC), and the spectral absorption information between the same substances is enhanced.

4)提取样本光谱信息,建立光谱库文件;4) Extract the spectral information of the sample and establish a spectral library file;

具体过程包括但不限于:The specific process includes but is not limited to:

所述烟叶和杂物样本包括:纯净烟叶样本,一类杂物样本,二类杂物样本,以及三类杂物样本;The tobacco leaf and foreign matter samples include: pure tobacco leaf samples, first-class foreign matter samples, second-class foreign matter samples, and third-class foreign matter samples;

分别选取烟叶和一、二、三类杂物感兴趣区(ROI),对所选样本感兴趣区域处光谱特征进行提取,并取得平均光谱;The tobacco leaves and the regions of interest (ROI) of the first, second and third types of debris were selected respectively, and the spectral features of the selected sample regions of interest were extracted, and the average spectrum was obtained;

建立光谱信息库文件,将样本ROI区域所取得的平均光谱导入库文件中进行保存。A spectral information library file is established, and the average spectrum obtained in the sample ROI area is imported into the library file for saving.

5)烟叶杂物分类的实现:5) The realization of tobacco leaf debris classification:

采用主成分分析法(PCA)对采集样品高光谱图像进行降维处理,然后利用光谱角匹配法(SAM)对目标光谱进行匹配,根据光谱特征向量判断该烟叶样品中是否混有杂物,并对不同样本进行标记。Principal component analysis (PCA) was used to reduce the dimension of the hyperspectral image of the collected sample, and then the target spectrum was matched by the spectral angle matching method (SAM). Label different samples.

具体过程如下:将混合有一、二、三类杂物的烟叶用短波高光谱成像仪进行图像采集,将所采集数据进行预处理后与光谱库内收录的光谱信息进行比对,通过主成分分析法(PCA)和光谱角匹配(SAM)算法对高光谱图像进行降维和烟叶、杂物的特征识别,其中包括:The specific process is as follows: the tobacco leaves mixed with one, two, and three types of impurities are imaged with a short-wave hyperspectral imager, and the collected data is preprocessed and compared with the spectral information recorded in the spectral library. Dimensionality reduction and feature recognition of tobacco leaves and sundries on hyperspectral images by using spectral angle matching (PCA) and spectral angle matching (SAM) algorithms, including:

首先,扫描并获得该样品的短波高光谱成像信息;First, scan and obtain short-wave hyperspectral imaging information of the sample;

其次,采用主成分分析法对采集样本进行降维处理;Secondly, adopt the principal component analysis method to reduce the dimension of the collected samples;

最后,通过光谱角匹配算法计算,根据光谱特征向量判断该烟叶样品中是否混合有杂物,并对不同样本进行标记。Finally, through the calculation of the spectral angle matching algorithm, it is judged whether the tobacco leaf sample is mixed with foreign matter according to the spectral characteristic vector, and different samples are marked.

步骤3)中的预处理方法还可采用均值中心化(mean centering)、标准化(autoscaling)、归一化(normalization)、标准正态变量变换(SNV)、导数、平滑去噪算法、小波变换中的一种The preprocessing method in step 3) can also adopt mean centering (mean centering), standardization (autoscaling), normalization (normalization), standard normal variable transformation (SNV), derivative, smoothing denoising algorithm, wavelet transform. a kind of

本发明提供的一种基于短波成像高光谱技术的烟叶杂物分类方法,结合了光谱技术和二维图像成像技术。其与现有技术相比,具有如下显著的进步:The invention provides a tobacco leaf debris classification method based on short-wave imaging hyperspectral technology, which combines spectral technology and two-dimensional image imaging technology. Compared with the existing technology, it has the following significant progress:

1.本发明利用短波成像高光谱技术对烟叶和杂物进行分类,建立包含有不同物质的光谱库。1. The present invention uses short-wave imaging hyperspectral technology to classify tobacco leaves and sundries, and establishes a spectral library containing different substances.

2.本发明实验过程不使用有毒有害化学品,简便、快捷、对样品无破坏性、对环境无污染。2. The experimental process of the present invention does not use toxic and harmful chemicals, which is simple, fast, non-destructive to the sample, and non-polluting to the environment.

3.本发明具有操作简便、快速、准确、成本低、效率高的优点。3. The present invention has the advantages of simple operation, rapidity, accuracy, low cost and high efficiency.

附图说明Description of drawings

图1为本发明所提供的一种基于短波高光谱成像技术的烟叶杂物识别分类方法的流程图;1 is a flowchart of a method for identifying and classifying sundries in tobacco leaves based on short-wave hyperspectral imaging technology provided by the present invention;

图2为烟叶混合杂物的光谱图。Figure 2 is the spectrum of the mixed impurities of tobacco leaves.

具体实施方式Detailed ways

下面结合附图对本发明的具体实施方式作进一步详细说明。The specific embodiments of the present invention will be further described in detail below with reference to the accompanying drawings.

实验仪器及参数Experimental equipment and parameters

四川双利合谱科技有限公司的GaiaSorter-N25E生产的SWIR型高光谱成像仪,波长范围为1000-2500nm,光谱分辨率为12nm,图像分辨率为384*288个像素,光谱仪帧数为400,曝光时间为20ms;实验平台使用GaiaSorter-Dual型暗箱系统,其包括卤钨灯光源,电控平台等。为避免外界杂散光源影响,高光谱图像采集过程在暗箱系统中进行。The SWIR hyperspectral imager produced by GaiaSorter-N25E of Sichuan Shuangli Hepu Technology Co., Ltd. has a wavelength range of 1000-2500nm, a spectral resolution of 12nm, an image resolution of 384*288 pixels, and a spectrometer frame number of 400. The exposure time is 20ms; the experimental platform uses a GaiaSorter-Dual type dark box system, which includes a halogen tungsten light source, an electronic control platform, and the like. In order to avoid the influence of external stray light sources, the hyperspectral image acquisition process is carried out in a dark box system.

步骤1:样本的采集Step 1: Collection of samples

为使所建立的光谱数据库更具有广泛的适用性,本实施例选取了2018年某中烟烤机出口片烟、某中烟提供的羽毛、布条、塑料绳等一、二类杂物。In order to make the established spectral database more widely applicable, this example selects the first and second types of sundries, such as cigarettes from a certain China Tobacco flue-curing machine in 2018, and feathers, cloth strips, and plastic ropes provided by a certain China Tobacco.

步骤2:样品制备与高光谱成像并进行黑白帧校正Step 2: Sample Preparation and Hyperspectral Imaging with Black and White Frame Correction

选择经除尘处理后表面无附着物的纯净烟叶、布条、羽毛、塑料绳平铺至黑色背板上,并在其所在位置进行标记。采用双利合谱公司生产的SWIR型高光谱成像仪,首先关闭镜头盖,采集反射率为零的全黑标定图像。再打开镜头盖,扫描标准白板,采集发射率为99.9%的全白标定图像,然后对烟叶和杂物样品进行高光谱成像。应用SpecView软件对图像进行采集,并根据所采集的全黑、全白标定图像对采集到的样本高光谱图像进行黑白帧校正,以减小外界杂散光所带来的噪声,最后将其作为原始光谱数据格式进行存储。Select pure tobacco leaves, cloth strips, feathers, and plastic ropes with no attachments on the surface after dust removal, lay them flat on the black backboard, and mark their locations. Using the SWIR hyperspectral imager produced by Shuangli Hepu Company, firstly close the lens cover and collect the full black calibration image with zero reflectivity. Then open the lens cover, scan a standard white board, acquire an all-white calibration image with an emissivity of 99.9%, and then perform hyperspectral imaging on tobacco leaves and debris samples. The SpecView software is used to collect the image, and black and white frame correction is performed on the collected hyperspectral image of the sample according to the collected all-black and all-white calibration images to reduce the noise caused by external stray light, and finally it is used as the original image. Spectral data format for storage.

黑白帧校正公式如下:The black and white frame correction formula is as follows:

式中:R-校正后高光谱图像;I-原始高光谱图像;B-关闭相机镜头采集的全黑标定图像;W-扫描白色校正板得到的全白标定图像。In the formula: R-corrected hyperspectral image; I-original hyperspectral image; B-all-black calibration image collected by turning off the camera lens; W-all-white calibration image obtained by scanning the white calibration plate.

步骤3:高光谱图像预处理及特征图像的获取Step 3: Hyperspectral image preprocessing and acquisition of characteristic images

为去除基线漂移和倾斜,去除噪声的影响,提高光谱曲线平滑度,使用Savitzky-Golay卷积平滑滤波算法。In order to remove baseline drift and tilt, remove the influence of noise, and improve the smoothness of the spectral curve, the Savitzky-Golay convolution smoothing filtering algorithm is used.

Savitzky-Golay滤波可以提高光谱曲线的平滑性,降低噪音的干扰。其卷积平滑的关键在于对于矩阵算子的求解。引入均方误差(MSE)选取合适窗宽n,滤波核中心点数m,其MSE值越小,噪声越小。设滤波窗口的宽度为n=2m+1,各测量点为x=(-m,-m+1…0…m-1,m),采用k-1次多项式对窗口内的数据点进行拟合:Savitzky-Golay filtering can improve the smoothness of the spectral curve and reduce the interference of noise. The key to its smooth convolution lies in the solution of the matrix operator. The mean square error (MSE) is introduced to select the appropriate window width n and the number of filter kernel center points m. The smaller the MSE value, the smaller the noise. Let the width of the filter window be n=2m+1, and each measurement point is x=(-m,-m+1...0...m-1,m), and use k-1 polynomial to fit the data points in the window. combine:

y=a0+a1x+a2x2+…+ak-1xk-1 y=a 0 +a 1 x+a 2 x 2 +…+a k-1 x k-1

存在n个上述方程,构成k元线性方程组。为使方程组有解,令n>k,通过最小二乘法求解拟合参数A,由此可以得到:There are n above equations, which form a system of k-element linear equations. In order to make the system of equations have a solution, let n>k, and solve the fitting parameter A by the least square method, which can be obtained:

用矩阵表示为:It is represented by a matrix as:

Y(2m+1)·1=X(2m+1)·k·Ak·1+E(2m+1)·1 Y (2m+1)·1 =X (2m+1)·k ·A k·1 +E (2m+1)·1

A的最小二乘解为:The least squares solution of A for:

Y的模型预测值或者滤波值为:Model predicted or filtered value of Y for:

利用多元散射校正(MSC)可以对采集样品光谱特征进行修正,提高光谱信噪比。The spectral characteristics of the collected samples can be corrected by using multiple scattering correction (MSC) to improve the spectral signal-to-noise ratio.

计算需校正光谱的平均光谱:Calculate the average spectrum of the spectrum to be corrected:

一元线性回归:Unary linear regression:

多元散射校正:Multivariate Scatter Correction:

式中,A-定标光谱数据矩阵,n-定标样品个数,p-采集光谱时的波长点数;-平均光谱矢量。In the formula, A is the calibration spectrum data matrix, n is the number of calibration samples, and p is the number of wavelength points when the spectrum is collected; - Average spectral vector.

步骤4:提取样本光谱信息,建立光谱库文件Step 4: Extract the spectral information of the sample and create a spectral library file

对包括纯净烟叶样本在内的一、二、三类杂物进行高光谱图像采集,在所采集的不同样本间选取感兴趣区域,提取样本间平均光谱信息,建立光谱库文件。The hyperspectral images were collected for the first, second and third types of debris including pure tobacco leaf samples, and the region of interest was selected between the different samples collected, the average spectral information between the samples was extracted, and the spectral library file was established.

步骤5:烟叶杂物分类的实现Step 5: Realization of Tobacco Leaf Debris Classification

采用主成分分析法(PCA)对采集样品高光谱图像进行降维处理,然后利用光谱角匹配法(SAM)算法对样本光谱进行匹配,从而实现对烟叶杂物分类的目的。Principal component analysis (PCA) was used to reduce the dimension of the collected sample hyperspectral images, and then the spectral angle matching (SAM) algorithm was used to match the sample spectra, so as to achieve the purpose of classifying tobacco leaf debris.

主成分分析法(PCA)选取特征向量实现多指标向量的降维。其降维模型如下:Principal Component Analysis (PCA) selects eigenvectors to achieve dimensionality reduction of multi-index vectors. Its dimensionality reduction model is as follows:

y1=l11x1+l12x2+…+l1nxn y 1 =l 11 x 1 +l 12 x 2 +...+l 1n x n

y2=l21x1+l22x2+…+l2nxn y 2 =l 21 x 1 +l 22 x 2 +...+l 2n x n

ym=lm1x1+lm2x2+…+lmnxn y m =l m1 x 1 +l m2 x 2 +…+l mn x n

式中x为n维多指标向量,y为处理后所得m维主成分向量。In the formula, x is the n-dimensional multi-index vector, and y is the m-dimensional principal component vector obtained after processing.

实施例中主成分分析特征值及方差贡献率见下表1。The principal component analysis eigenvalues and variance contribution rates in the examples are shown in Table 1 below.

表1主成分分析特征值及方差贡献率Table 1 Principal component analysis eigenvalues and variance contribution rate

光谱角匹配法利用光谱库中参考光谱匹配未知样品,利用设定狭义光谱角阈值大小对其进行判断。其公式为:The spectral angle matching method uses the reference spectrum in the spectral library to match the unknown sample, and judges it by setting the narrow-sense spectral angle threshold. Its formula is:

式中T为标准主成分得分矢量,R为参考主成分得分矢量为(T、R均非零向量)。where T is the standard principal component score vector, and R is the reference principal component score vector (both T and R are non-zero vectors).

进一步的,将线推扫成像的高光谱图像光谱信息与监督分类所建立的光谱库进行比对,得到识别结果。Further, the spectral information of the hyperspectral image of the line push-broom imaging is compared with the spectral library established by the supervised classification to obtain the identification result.

本实施例提供的基于高光谱成像的分类识别方法,至少包括如下有益效果:The classification and identification method based on hyperspectral imaging provided by this embodiment at least includes the following beneficial effects:

1.利用Savitzky-Golay滤波,提高光谱曲线的平滑性,降低了噪音干扰;利用多元散射校正,提升光谱信号信噪比;二者为之后的图像处理工作提供基础1. Use Savitzky-Golay filtering to improve the smoothness of the spectral curve and reduce noise interference; use multiple scattering correction to improve the signal-to-noise ratio of the spectral signal; both provide the basis for subsequent image processing work

2.通过高光谱图像采集,得到不同物质的光谱信息。建立光谱库文件,保存样本光谱信息。2. Obtain spectral information of different substances through hyperspectral image acquisition. Create a spectral library file to save sample spectral information.

3.利用主成分分析的方法进行降维,降低图像计算量,消除非重要特征所带来的影响,大大减少计算时间。通过光谱角匹配算法,利用光谱库中参考光谱匹配未知样品,利用设定狭义光谱角阈值大小对其进行判断,完成烟叶和杂物有效的分类识别,从而得到较为准确的分类效果。3. The method of principal component analysis is used to reduce the dimension, reduce the amount of image calculation, eliminate the influence of non-important features, and greatly reduce the calculation time. Through the spectral angle matching algorithm, the reference spectrum in the spectral library is used to match the unknown sample, and the narrow-sense spectral angle threshold is used to judge it, so as to complete the effective classification and identification of tobacco leaves and sundries, so as to obtain a more accurate classification effect.

4.相较于传统人工除杂来讲,提升了工作效率与检测速度。同时,利用高光谱成像技术可以获得更为精细的成分信息,提高分辨的准确率。4. Compared with the traditional manual impurity removal, the work efficiency and detection speed are improved. At the same time, the use of hyperspectral imaging technology can obtain finer composition information and improve the accuracy of resolution.

上面结合实施例对本发明作了详细说明,而并非是对本发明的实施方式的限定,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明原理和宗旨的前提下做出各种变化、修改、替换与变形。本发明的范围由权利要求及其等同物限定。The present invention has been described in detail above in conjunction with the embodiments, rather than limiting the embodiments of the present invention. Within the scope of knowledge possessed by those of ordinary skill in the art, it can also be made without departing from the principles and purposes of the present invention. Various changes, modifications, substitutions and deformations. The scope of the invention is defined by the claims and their equivalents.

Claims (10)

1. A tobacco non-smoke substance classification method based on a hyperspectral imaging technology is characterized by comprising the following steps: tobacco leaves and impurities are classified by utilizing a short-wave hyperspectral imaging technology, a spectrum library comprising different substances is established, hyperspectral image data are collected for a sample to be detected, a reference spectrum in the spectrum library is used for matching the sample to be detected and judging the sample to be detected, and then effective classification and identification of the tobacco leaves and the impurities are completed.
2. The hyperspectral imaging technology-based tobacco non-smoke substance classification method according to claim 1, characterized in that: the method comprises the following specific steps:
1) Collecting samples: obtaining samples of first-class sundries, second-class sundries and third-class sundries including pure tobacco leaves;
2) Sample preparation and hyperspectral imaging are carried out, and black and white frame correction is carried out;
3) preprocessing a hyperspectral image and acquiring a characteristic image;
4) Extracting sample spectrum information and establishing a spectrum library file;
5) the tobacco leaf sundries are classified:
performing dimensionality reduction on a hyperspectral image of an acquired sample by adopting a Principal Component Analysis (PCA), matching a target spectrum by utilizing a spectral angle matching method (SAM), judging whether impurities are mixed in the tobacco sample or not according to a spectral feature vector, and marking different samples.
3. The tobacco non-smoke substance classification method based on the hyperspectral imaging technology according to claim 2, characterized in that: and 2) using a halogen tungsten lamp as an illumination light source, and performing hyperspectral image acquisition on the obtained tobacco leaf and sundry samples to obtain hyperspectral images of the samples. In order to reduce the noise influence, the hyperspectral image is corrected by a black and white frame correction formula as follows:
In the formula: r-a corrected hyperspectral image; i-an original hyperspectral image; b, closing the all-black image collected by the camera lens; w-scanning the white correction plate to obtain a full white image.
4. The tobacco non-smoke substance classification method based on the hyperspectral imaging technology according to claim 2, characterized in that: in step 3), in order to improve the signal-to-noise ratio of the data, the hyperspectral image data is preprocessed, and the preprocessing method includes but is not limited to:
By using a Savitzky-Golay convolution smoothing filtering algorithm, baseline drift and inclination are removed, noise is removed, and smoothness of a spectral curve is improved;
And then reducing the scattering effect on the surface of the object by Multivariate Scattering Correction (MSC) to enhance the spectral absorption information among the same substances.
5. the tobacco non-smoke substance classification method based on the hyperspectral imaging technology according to claim 2, characterized in that: in step 4), the specific processes include but are not limited to:
The tobacco leaf and sundries sample comprises: a pure tobacco leaf sample, a first class sundry sample, a second class sundry sample and a third class sundry sample;
Respectively selecting tobacco leaves and regions of interest (ROI) of the first, second and third sundries, extracting spectral characteristics at the selected region of interest of the sample, and obtaining an average spectrum;
and establishing a spectrum information library file, and importing the average spectrum obtained from the ROI area of the sample into the library file for storage.
6. The tobacco non-smoke substance classification method based on the hyperspectral imaging technology according to claim 2, characterized in that: the specific process in step 5) is as follows: the tobacco leaves mixed with the first, second and third impurities are subjected to image acquisition by a short wave hyperspectral imager, the acquired data are preprocessed and then compared with spectral information recorded in a spectral library, and the hyperspectral images are subjected to dimensionality reduction and characteristic identification of the tobacco leaves and the impurities by a Principal Component Analysis (PCA) and a Spectral Angle Matching (SAM) algorithm, wherein the method comprises the following steps:
Firstly, scanning and obtaining short-wave hyperspectral imaging information of the sample;
secondly, performing dimensionality reduction on the collected sample by adopting a principal component analysis method;
And finally, calculating by a spectral angle matching algorithm, judging whether impurities are mixed in the tobacco leaf sample according to the spectral feature vector, and marking different samples.
7. The hyperspectral imaging technology-based tobacco non-smoke substance classification method according to claim 4, characterized in that: the preprocessing method in step 3 may further adopt one of mean centering (mean centering), normalization (normalization), standard normal variable transformation (SNV), derivative, smooth denoising algorithm, and wavelet transformation.
8. The tobacco non-smoke substance classification method based on hyperspectral imaging technology according to claim 1 or 2, characterized in that: the sample surface remained dry and clean with no other attachments.
9. the tobacco non-smoke substance classification method based on hyperspectral imaging technology according to claim 1 or 2, characterized in that: the hyperspectral imager is used for hyperspectral imaging, the wavelength range is 1000-2500nm, the spectral resolution is 12nm, the image resolution is 384 × 288 pixels, and the spectrometer frame number is 400.
10. A tobacco non-smoke substance classification method based on hyperspectral imaging technology according to claim 2 or 5, characterized in that: the sundries comprise metal, feather and plastic; the second class of sundries comprises paper, stones, hemp ropes and glass; the three kinds of impurities comprise non-tobacco leaves, seeds and bamboo sticks.
CN201910850623.0A 2019-09-10 2019-09-10 A classification method of tobacco non-smoke substances based on hyperspectral imaging technology Pending CN110542658A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910850623.0A CN110542658A (en) 2019-09-10 2019-09-10 A classification method of tobacco non-smoke substances based on hyperspectral imaging technology

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910850623.0A CN110542658A (en) 2019-09-10 2019-09-10 A classification method of tobacco non-smoke substances based on hyperspectral imaging technology

Publications (1)

Publication Number Publication Date
CN110542658A true CN110542658A (en) 2019-12-06

Family

ID=68713013

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910850623.0A Pending CN110542658A (en) 2019-09-10 2019-09-10 A classification method of tobacco non-smoke substances based on hyperspectral imaging technology

Country Status (1)

Country Link
CN (1) CN110542658A (en)

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112189877A (en) * 2020-10-13 2021-01-08 燕山大学 On-line detection method for tobacco shred impurities in tobacco production line
CN112613413A (en) * 2020-12-25 2021-04-06 平安国际智慧城市科技股份有限公司 Perishable garbage classification quality determination method and device and computer readable storage medium
CN112697725A (en) * 2020-12-09 2021-04-23 云南省烟草农业科学研究院 Method and system for predicting total sugar content of tobacco leaves in field based on hyperspectrum of unmanned aerial vehicle
CN112697726A (en) * 2020-12-09 2021-04-23 云南省烟草农业科学研究院 Hyperspectral field tobacco nicotine content prediction method and system based on unmanned aerial vehicle
CN113210264A (en) * 2021-05-19 2021-08-06 江苏鑫源烟草薄片有限公司 Method and device for removing tobacco impurities
CN113237836A (en) * 2021-05-10 2021-08-10 福建省烟草公司南平市公司 Flue-cured tobacco leaf moisture content estimation method based on hyperspectral image
CN113390799A (en) * 2021-06-24 2021-09-14 中国烟草总公司郑州烟草研究院 Method for identifying and detecting stems in tobacco leaves
CN113655072A (en) * 2021-08-20 2021-11-16 福建中烟工业有限责任公司 Method, apparatus and computer readable medium for detecting contaminants on a sample surface
CN113681156A (en) * 2021-06-24 2021-11-23 中国烟草总公司郑州烟草研究院 Cigarette raw material slitting device and method based on ultraviolet laser
CN114155385A (en) * 2021-11-30 2022-03-08 江苏鑫源烟草薄片有限公司 A RGB and Hyperspectral Tobacco Recognition Method Based on Deep Learning
CN114486774A (en) * 2021-12-31 2022-05-13 中科谱光(郑州)应用科学技术研究院有限公司 Artificial intelligence algorithm matching method based on hyperspectral big data
CN114659994A (en) * 2022-05-23 2022-06-24 武汉大学 Tobacco foreign matter extraction and detection system
CN115128033A (en) * 2022-07-04 2022-09-30 厦门烟草工业有限责任公司 Tobacco leaf detection method, device and system and storage medium
CN115165772A (en) * 2022-07-25 2022-10-11 江苏大学 High-spectrum-based detection device and detection method for impurity content of grains
CN115165795A (en) * 2022-06-28 2022-10-11 中国烟草总公司郑州烟草研究院 A method and device for predicting chemical composition of tobacco based on spectral conversion
CN115356297A (en) * 2022-08-26 2022-11-18 中国烟草总公司郑州烟草研究院 A method and device for identifying slightly green flue-cured tobacco leaves
CN115356275A (en) * 2022-08-26 2022-11-18 中国烟草总公司郑州烟草研究院 Tobacco leaf grade evaluation method
CN117611828A (en) * 2024-01-19 2024-02-27 云南烟叶复烤有限责任公司 Non-smoke sundry detection method based on hyperspectral image segmentation technology
CN118837302A (en) * 2024-06-28 2024-10-25 湖南烟叶复烤有限公司 Specific multispectral information acquisition system for detecting tobacco flakes and sundries and manufacturing method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105842173A (en) * 2016-06-06 2016-08-10 南京大学 Method for identifying hyperspectral material
CN106501208A (en) * 2016-09-20 2017-03-15 广西中烟工业有限责任公司 A kind of tobacco style similitude sorting technique based near infrared light spectrum signature
CN107843593A (en) * 2017-10-13 2018-03-27 上海工程技术大学 A kind of textile material recognition methods and system based on high light spectrum image-forming technology
CN207379913U (en) * 2017-10-19 2018-05-18 新疆农业大学 A cotton impurity detection device based on reflection and transmission hyperspectral imaging

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105842173A (en) * 2016-06-06 2016-08-10 南京大学 Method for identifying hyperspectral material
CN106501208A (en) * 2016-09-20 2017-03-15 广西中烟工业有限责任公司 A kind of tobacco style similitude sorting technique based near infrared light spectrum signature
CN107843593A (en) * 2017-10-13 2018-03-27 上海工程技术大学 A kind of textile material recognition methods and system based on high light spectrum image-forming technology
CN207379913U (en) * 2017-10-19 2018-05-18 新疆农业大学 A cotton impurity detection device based on reflection and transmission hyperspectral imaging

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
赵冬娥等: "高光谱成像的垃圾分类识别研究", 《光谱学与光谱分析》 *

Cited By (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112189877A (en) * 2020-10-13 2021-01-08 燕山大学 On-line detection method for tobacco shred impurities in tobacco production line
CN112697725A (en) * 2020-12-09 2021-04-23 云南省烟草农业科学研究院 Method and system for predicting total sugar content of tobacco leaves in field based on hyperspectrum of unmanned aerial vehicle
CN112697726A (en) * 2020-12-09 2021-04-23 云南省烟草农业科学研究院 Hyperspectral field tobacco nicotine content prediction method and system based on unmanned aerial vehicle
CN112613413A (en) * 2020-12-25 2021-04-06 平安国际智慧城市科技股份有限公司 Perishable garbage classification quality determination method and device and computer readable storage medium
CN113237836A (en) * 2021-05-10 2021-08-10 福建省烟草公司南平市公司 Flue-cured tobacco leaf moisture content estimation method based on hyperspectral image
CN113210264A (en) * 2021-05-19 2021-08-06 江苏鑫源烟草薄片有限公司 Method and device for removing tobacco impurities
CN113210264B (en) * 2021-05-19 2023-09-05 江苏鑫源烟草薄片有限公司 Tobacco sundry removing method and device
CN113390799A (en) * 2021-06-24 2021-09-14 中国烟草总公司郑州烟草研究院 Method for identifying and detecting stems in tobacco leaves
CN113681156A (en) * 2021-06-24 2021-11-23 中国烟草总公司郑州烟草研究院 Cigarette raw material slitting device and method based on ultraviolet laser
CN113681156B (en) * 2021-06-24 2023-08-08 中国烟草总公司郑州烟草研究院 Device and method for cutting raw materials for cigarettes based on ultraviolet laser
CN113655072B (en) * 2021-08-20 2023-12-26 福建中烟工业有限责任公司 Method, apparatus and computer readable medium for detecting contaminants on a sample surface
CN113655072A (en) * 2021-08-20 2021-11-16 福建中烟工业有限责任公司 Method, apparatus and computer readable medium for detecting contaminants on a sample surface
CN114155385A (en) * 2021-11-30 2022-03-08 江苏鑫源烟草薄片有限公司 A RGB and Hyperspectral Tobacco Recognition Method Based on Deep Learning
CN114155385B (en) * 2021-11-30 2022-09-27 江苏鑫源烟草薄片有限公司 A RGB and Hyperspectral Tobacco Recognition Method Based on Deep Learning
CN114486774A (en) * 2021-12-31 2022-05-13 中科谱光(郑州)应用科学技术研究院有限公司 Artificial intelligence algorithm matching method based on hyperspectral big data
CN114659994B (en) * 2022-05-23 2022-08-05 武汉大学 Tobacco foreign matter extraction and detection system
CN114659994A (en) * 2022-05-23 2022-06-24 武汉大学 Tobacco foreign matter extraction and detection system
CN115165795A (en) * 2022-06-28 2022-10-11 中国烟草总公司郑州烟草研究院 A method and device for predicting chemical composition of tobacco based on spectral conversion
CN115128033A (en) * 2022-07-04 2022-09-30 厦门烟草工业有限责任公司 Tobacco leaf detection method, device and system and storage medium
CN115165772A (en) * 2022-07-25 2022-10-11 江苏大学 High-spectrum-based detection device and detection method for impurity content of grains
CN115356297A (en) * 2022-08-26 2022-11-18 中国烟草总公司郑州烟草研究院 A method and device for identifying slightly green flue-cured tobacco leaves
CN115356275A (en) * 2022-08-26 2022-11-18 中国烟草总公司郑州烟草研究院 Tobacco leaf grade evaluation method
CN117611828A (en) * 2024-01-19 2024-02-27 云南烟叶复烤有限责任公司 Non-smoke sundry detection method based on hyperspectral image segmentation technology
CN117611828B (en) * 2024-01-19 2024-05-24 云南烟叶复烤有限责任公司 Non-smoke sundry detection method based on hyperspectral image segmentation technology
CN118837302A (en) * 2024-06-28 2024-10-25 湖南烟叶复烤有限公司 Specific multispectral information acquisition system for detecting tobacco flakes and sundries and manufacturing method

Similar Documents

Publication Publication Date Title
CN110542658A (en) A classification method of tobacco non-smoke substances based on hyperspectral imaging technology
Zhao et al. Application of hyperspectral imaging and chemometrics for variety classification of maize seeds
Singh et al. Nondestructive identification of barley seeds variety using near‐infrared hyperspectral imaging coupled with convolutional neural network
Yu et al. Identification of crack features in fresh jujube using Vis/NIR hyperspectral imaging combined with image processing
Lu et al. Hyperspectral imaging for safety inspection of food and agricultural products
Yu et al. Application of visible and near-infrared hyperspectral imaging for detection of defective features in loquat
CN110763698B (en) Hyperspectral citrus leaf disease identification method based on characteristic wavelength
CN109540831B (en) Identification method of wolfberry varieties based on hyperspectral imaging technology
CN110243806B (en) Mixture component identification method based on similarity under Raman spectrum
Femenias et al. Hyperspectral imaging for the classification of individual cereal kernels according to fungal and mycotoxins contamination: A review
CN108802002B (en) Silkworm egg Raman spectrum model construction method for rapidly identifying and removing diapause without damage
CN113390799B (en) Method for identifying and detecting stems in tobacco leaves
CN107527326A (en) A kind of wheat scab diagnostic method based on high light spectrum image-forming
CN106290238A (en) A kind of apple variety method for quick identification based on high light spectrum image-forming
CN103954570A (en) Food insect attack degree distinguishing method based on spectral imaging technology
CN110376202B (en) Tea tree anthracnose lesion identification method based on imaging hyperspectral technology
CN109001218A (en) Apple surface defect quick nondestructive recognition methods based on high light spectrum image-forming technology
CN111443043A (en) A method for quality detection of walnut kernels based on hyperspectral images
Guo et al. Detection of foreign materials on surface of ginned cotton by hyper-spectral imaging
Chen et al. Quality detection and variety classification of pecan seeds using hyperspectral imaging technology combined with machine learning
CN104914052B (en) A kind of light spectrum image-forming detection method for the peanut that goes mouldy
Cheng et al. Identification of wine grape varieties based on near-infrared hyperspectral imaging
CN117054358A (en) Fresh tobacco leaf pesticide residue type detection method and system based on hyperspectral technology
CN113791040B (en) A method and system for detecting heavy metals in soil
CN107944426B (en) Wheat leaf powdery mildew spot marking method based on combination of texture filtering and two-dimensional spectral feature space discrimination

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
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20191206

WD01 Invention patent application deemed withdrawn after publication