CN110503003A - 基于led阵列与卷积神经网络的茶叶品种鉴定装置与方法 - Google Patents
基于led阵列与卷积神经网络的茶叶品种鉴定装置与方法 Download PDFInfo
- Publication number
- CN110503003A CN110503003A CN201910687893.4A CN201910687893A CN110503003A CN 110503003 A CN110503003 A CN 110503003A CN 201910687893 A CN201910687893 A CN 201910687893A CN 110503003 A CN110503003 A CN 110503003A
- Authority
- CN
- China
- Prior art keywords
- led
- neural networks
- convolutional neural
- data
- tealeaves
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
- G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
- G01N21/64—Fluorescence; Phosphorescence
- G01N21/6486—Measuring fluorescence of biological material, e.g. DNA, RNA, cells
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/13—Satellite images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/60—Type of objects
- G06V20/68—Food, e.g. fruit or vegetables
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Data Mining & Analysis (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Biomedical Technology (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Software Systems (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Mathematical Physics (AREA)
- Computational Linguistics (AREA)
- Biophysics (AREA)
- Computing Systems (AREA)
- Evolutionary Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Computational Biology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Immunology (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Chemical & Material Sciences (AREA)
- Astronomy & Astrophysics (AREA)
- Remote Sensing (AREA)
- Multimedia (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明涉及一种基于LED阵列与卷积神经网络的茶叶品种鉴定装置与方法。本发明中的LED阵列采用多只LED,其中主LED中心波长在370~380nm之间,该波长的LED能够激发出茶叶中茶多酚、黄酮、花青素、叶绿素的荧光,最大程度地提供有用信息,辅助LED中心波长在360~500nm之间,采用轮流点亮的方式,受控于计算控制单元。迷你光谱仪用于采集待鉴定茶叶的光谱数据,并将所述的光谱数据发送至计算控制单元。计算控制单元将光谱数据转化为二维数据后,送入已训练好的卷积神经网络,得到待鉴定茶叶的品种。该装置轻便可手持,且能够快速、准确地对茶叶品种进行鉴定。
Description
技术领域
本发明属于光学传感技术领域以及分析检测仪器技术领域,涉及一种基于多波长LED阵列与卷积神经网络的手持式茶叶品种鉴定方法与鉴定装置。
背景技术
茶叶是世界上最受欢迎的非酒精饮料,在世界各地均有茶叶的种植与销售。根据发酵程度的不同,茶叶分为绿茶、红茶、青茶、白茶、黄茶与黑茶。根据茶叶种类的不同,又有西湖龙井茶、洞庭碧螺春、黄山毛峰、庐山云雾茶、六安瓜片、君山银针、信阳毛尖、武夷岩茶、安溪铁观音、祁门红茶等分类。根据采摘时间的不同,又有清明茶、谷雨茶、夏茶、秋茶的分类。品种、时令、产地等都会对茶叶的品质产生影响,并最终影响茶叶的销售价格。消费者在面对市场上的茶叶时,容易受到产品包装、商家鼓吹的影响,将次茶误认为好茶,将“浙江龙井”误认为“狮峰龙井”。因此,快速、轻便、准确的茶叶品种鉴定装置具有很大的社会价值与经济价值。
申请号201811050174.3与201810444937.6的发明专利分别公布了基于tensorflow与全自动支持向量机的茶叶品种识别系统与方法,其核心是基于机器视觉与图像处理。这种方法在采样时需要准确的照明角度与照相位置,不然所得结果会有很大的偏差。
申请号201510258210.5的发明专利公布了一种基于激光诱导击穿光谱技术的茶叶分类鉴别方法。但激光诱导击穿光谱需要使用大功率激光器,有一定的安全防护要求,且多台光谱仪联用,系统体积庞大,不适用于线下消费场景。
申请号201810276203.9、201810335781.8、201710342868.3、201810468730.2公布的茶叶分类方法,基于红外光谱与不同的数据处理方法。其缺点是红外光谱采集所使用光谱仪价格在几十万,十分昂贵。
激光诱导荧光技术或者LED诱导荧光技术能够激发茶叶中叶绿素、茶多酚、黄铜等物质的荧光,因而可以通过这些物质在荧光谱中的相对强度对茶叶进行分类。申请号为201310463824.8和201510391628.3的发明专利分别公布了基于主成分分析以及二维相关光谱的茶叶品质检测方法。这两种方案的系统由于使用了光纤收光至微型光谱仪的结构,因而测量时由于光纤弯折状态的变化会导致光谱曲线的变化,进而影响测量精度。同时,这两种方法的正确率有进一步提升的空间。
发明内容
本发明的目的在于提出一种基于多波长LED阵列与卷积神经网络的快速、轻便、准确的手持式茶叶品种鉴定装置与方法。
本发明所采用的技术方案:
基于LED阵列与卷积神经网络的茶叶品种鉴定装置包括LED座、LED阵列、滤波片、迷你光谱仪和计算控制单元。
所述LED座为LED提供支撑,使得LED照射方向对准待鉴定的茶叶,LED座中心有圆孔,用于安装滤波片与迷你光谱仪。
所述LED阵列采用多只LED,其中主LED中心波长在370~380nm之间,该波长的LED能够激发出茶叶中茶多酚、黄酮、花青素、叶绿素的荧光,最大程度地提供有用信息,辅助LED中心波长在360~500nm之间,采用轮流点亮的方式,受控于计算控制单元。
所述迷你光谱仪用于采集待鉴定茶叶的光谱数据,并将所述的光谱数据发送至计算控制单元。
所述的计算控制单元将光谱数据转化为二维数据后,送入已训练好的卷积神经网络,得到待鉴定茶叶的品种。
进一步,还包括用于显示分类结果的显示屏。
进一步,所述的滤波片为长通滤波片,截止波长为500nm,因为茶多酚、黄酮、花青素的荧光在500~650nm之间,叶绿素的荧光在650~800nm之间,因而500nm以后的光需要尽可能收集,而之前的光需要加以抑制。
基于LED阵列与卷积神经网络的茶叶品种鉴定方法,上述装置,包括以下步骤:
步骤1、使用茶叶品种鉴定装置将茶叶覆盖,并启动装置。
步骤2、茶叶种类鉴定装置中的N只LED依次点亮熄灭。在单只LED点亮后,迷你光谱仪记录光谱数据,共T个,然后关闭LED。在N只LED都记录完毕后,记录一组无LED点亮情况下的背景光谱数据。所得数据维数为(N+1)*T。
步骤3、对(N+1)*T中前N行数据分别减去背景光谱数据,并进行平滑处理。对光谱数据进行归一化。所得数据维数为N*T。
步骤4、对k种茶叶进行测量,每种茶叶测量M次,每次所得为N*T二维矩阵。对所有数据使用卷积神经网络进行训练,建立判别模型。所述卷积神经网络由输入层,卷积层,函数激活层,全连接层构成。针对由多条一维光谱数据组合而成的二维光谱数据,其波长位置信息十分确定,因而其卷积神经网络结构不包括池化层。
步骤5、在鉴定时,重复步骤1~3。在计算控制单元中调用训练好的判别模型进行计算,得到分类结果,并将结果显示在屏幕上。
进一步,步骤1中的茶叶为单层或多层覆盖,但不超过3mm厚。
进一步,步骤2中迷你光谱仪记录数据范围为500~850nm。
进一步,步骤3中归一化方式为每一行数据除以本行650nm至700nm范围内的最大值。
进一步,步骤4中卷积神经网络的卷积层的卷积核采用N*N矩阵。因而训练得到的卷积核能够得到N只LED激发荧光之间的协同关系,而非小于N*N卷积核得到的离散的、局部的激发-发射荧光关系。
本发明的有益效果。
(1)LED的选择,尤其是主LED波长的确定,能够最大程度激发茶叶中特征物质,如茶多酚的荧光信息。
(2)迷你光谱仪直接对空间光进行光谱分析,不需要结合光纤等收光结构。迷你光谱仪的使用,使得系统进一步集成,可以做到手持式小型化设备,且省去了光纤的使用,避免了基于微型光谱仪的茶叶分类设备中由于光纤弯曲产生的光谱差异。
(3)针对茶叶的荧光特性改进了卷积神经网络结构,如使用多只LED诱导产生的荧光数据组合成为二维荧光数据作为输入,卷积神经网络结构不设置池化层等,使得准确度显著提高。
(4)荧光具有发散性,因而相较基于可见光/近红外光谱结合算法的茶叶分类技术,基于荧光的茶叶分类计数降低了对样本摆放状态的要求,更适用于实际生活中由非专业的消费者使用。
附图说明
图1为本发明的结构示意图;
图2为实施例中LED座的俯视图;
图3为实施例中LED的光谱。
具体实施方式
为了使专利局的审查员尤其是公众能够更加清楚地理解本发明的技术实质和有益效果,申请人将在下面以实施例的方式作详细说明,但是对实施例的描述均不是对本发明方案的限制,任何依据本发明构思所作出的仅仅为形式上的而非实质性的等效变换都应视为本发明的技术方案范畴。
如图1和图2所示,一种基于多波长LED阵列与卷积神经网络的手持式茶叶品种鉴定装置,能够实现对茶叶样品的种类进行快速准确的鉴定。该装置包括LED座1、LED阵列2、滤波片3、迷你光谱仪4、计算控制单元5、显示屏6。在本实施例中,LED阵列由7只LED组成(对应于图2中的标号11~17),其中心波长分别是371nm、381nm、394nm、404nm、412nm、423nm以及431nm。其中,371nmLED符合主LED波长要求,能够激发茶叶中茶多酚、黄酮、花青素、叶绿素的荧光。辅助LED的波长较为均匀地分布在380~430nm之间,且其光谱末端不会超过500nm,不会对特征物质的荧光带来影响;7只LED的光谱如图3所示。本实施例中的迷你光谱仪采用滨松Hamamatsu Mini-spectrometer C12880MA。该迷你光谱仪由狭缝、反射凹面闪耀光栅、高灵敏度CMOS线性图像传感器组成,其大小有别于传统的微型光谱仪,只有手指头大小。计算控制单元可以采用单片机、ARM、FPGA等嵌入式芯片以及稳流电路,主要实现LED的闪烁控制、迷你光谱仪的采样控制以及数据计算。
为实现茶叶分类,需要先对众多样品数据进行训练,得到数据库与判别模型。本次分类使用380元每斤的龙井、480元每斤的龙井、580元每斤的龙井、越南绿茶、乌牛早绿茶、黄山毛峰绿茶、慢城白茶、祁门红茶、立顿红茶作为样本。训练过程包含以下几个步骤:
步骤1、使用茶叶品种鉴定装置将茶叶覆盖,即茶叶放底下,装置罩在上面,并启动装置。
步骤2、茶叶种类鉴定装置所属7只LED依次点亮熄灭。在单只LED点亮后,迷你光谱仪记录光谱数据,共837个,然后关闭LED。在7只LED都记录完毕后,记录一组无LED点亮情况下的背景光谱数据,所得数据维数为8*837。
步骤3、对8*837中前7行数据分别减去第8行背景光谱数据,并使用s-g二阶滤波进行平滑处理。对光谱数据除以685nm处的最大值,进行归一化,所得数据维数为7*837。
步骤4、对9种茶叶进行测量,每种茶叶测量200次,每次所得为7*837二维矩阵。对所有数据使用卷积神经网络进行训练,建立判别模型。卷积神经网络由输入层,卷积层,函数激活层,全连接层构成。
步骤5、完成数据库与判别模型后,对实际样本进行检测,检测过程重复步骤1~3。在计算控制单元中调用训练好的判别模型进行计算,得到分类结果。结果会显示在屏幕上。
若使用单只LED激发荧光数据结合主成分分析算法与回归模型、7只LED激发荧光数据结合主成分分析算法与回归模型来训练,所得到的准确率只有80%。在采用本实施例的数据处理方法后,准确率可以提升到99.7%。
消费者在采用本实施例的设备后,可以很容易的区分出不同等级的龙井,又或者是不同等级的红茶,从而对所购茶叶的大体价位会有一个准确的认识,不会由于自身知识的缺乏而被商家蒙蔽。
Claims (8)
1.基于LED阵列与卷积神经网络的茶叶品种鉴定装置,包括LED座、LED阵列、滤波片、迷你光谱仪和计算控制单元,其特征在于:
所述LED座为LED提供支撑,使得LED照射方向对准待鉴定的茶叶,LED座中心有圆孔,用于安装滤波片与迷你光谱仪;
所述LED阵列采用多只LED,其中主LED中心波长在370~380nm之间,辅助LED中心波长在360~500nm之间,采用轮流点亮的方式,受控于计算控制单元;
所述迷你光谱仪用于采集待鉴定茶叶的光谱数据,并将所述的光谱数据发送至计算控制单元;
所述的计算控制单元将光谱数据转化为二维数据后,送入已训练好的卷积神经网络,得到待鉴定茶叶的品种。
2.根据权利要求1所述的基于LED阵列与卷积神经网络的茶叶品种鉴定装置,其特征在于:还包括用于显示分类结果的显示屏。
3.根据权利要求1或2所述的基于LED阵列与卷积神经网络的茶叶品种鉴定装置,其特征在于:所述的滤波片为长通滤波片,截止波长为500nm。
4.基于LED阵列与卷积神经网络的茶叶品种鉴定方法,使用权利要求1所述的装置,其特征在于包括以下步骤:
步骤1、使用茶叶品种鉴定装置将茶叶覆盖,并启动装置;
步骤2、茶叶种类鉴定装置中的N只LED依次点亮熄灭;在单只LED点亮后,迷你光谱仪记录光谱数据,共T个;
关闭LED,在N只LED都记录完毕后,记录一组无LED点亮情况下的背景光谱数据,所得数据维数为(N+1)*T;
步骤3、对(N+1)*T中前N行数据分别减去背景光谱数据,并进行平滑处理;对光谱数据进行归一化,所得数据维数为N*T;
步骤4、对k种茶叶进行测量,每种茶叶测量M次,每次所得为N*T二维矩阵;
对所有数据使用卷积神经网络进行训练,建立判别模型;
所述卷积神经网络由输入层,卷积层,函数激活层,全连接层构成;
步骤5、在鉴定时,重复步骤1~3,在计算控制单元中调用训练好的判别模型进行计算,得到分类结果,并将结果显示在屏幕上。
5.根据权利要求4所述的基于LED阵列与卷积神经网络的茶叶品种鉴定方法,其特征在于:步骤1中的茶叶为单层或多层覆盖,但不超过3mm厚。
6.根据权利要求4所述的基于LED阵列与卷积神经网络的茶叶品种鉴定方法,其特征在于:步骤2中迷你光谱仪记录数据范围为500~850nm。
7.根据权利要求4所述的基于LED阵列与卷积神经网络的茶叶品种鉴定方法,其特征在于:步骤3中归一化方式为每一行数据除以本行650nm至700nm范围内的最大值。
8.根据权利要求4所述的基于LED阵列与卷积神经网络的茶叶品种鉴定方法,其特征在于:步骤4中卷积神经网络的卷积层的卷积核采用N*N矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910687893.4A CN110503003A (zh) | 2019-07-29 | 2019-07-29 | 基于led阵列与卷积神经网络的茶叶品种鉴定装置与方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910687893.4A CN110503003A (zh) | 2019-07-29 | 2019-07-29 | 基于led阵列与卷积神经网络的茶叶品种鉴定装置与方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110503003A true CN110503003A (zh) | 2019-11-26 |
Family
ID=68587688
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910687893.4A Pending CN110503003A (zh) | 2019-07-29 | 2019-07-29 | 基于led阵列与卷积神经网络的茶叶品种鉴定装置与方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110503003A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111222543A (zh) * | 2019-12-19 | 2020-06-02 | 津海威视技术(天津)有限公司 | 物质识别方法和设备及计算机可读存储介质 |
CN112945920A (zh) * | 2021-02-01 | 2021-06-11 | 长江大学 | 一种绿茶等级鉴别方法、装置及计算机可读存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103175815A (zh) * | 2013-03-06 | 2013-06-26 | 浙江大学 | 多波长led诱导荧光的茶叶品质无损检测方法及装置 |
CN104977258A (zh) * | 2015-07-07 | 2015-10-14 | 江苏鼎云信息科技有限公司 | 基于二维相关光谱的茶叶/化妆品等品质检测方法 |
CN105973858A (zh) * | 2016-06-13 | 2016-09-28 | 宜春学院 | 一种中药品质自动检测系统 |
-
2019
- 2019-07-29 CN CN201910687893.4A patent/CN110503003A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103175815A (zh) * | 2013-03-06 | 2013-06-26 | 浙江大学 | 多波长led诱导荧光的茶叶品质无损检测方法及装置 |
CN104977258A (zh) * | 2015-07-07 | 2015-10-14 | 江苏鼎云信息科技有限公司 | 基于二维相关光谱的茶叶/化妆品等品质检测方法 |
CN105973858A (zh) * | 2016-06-13 | 2016-09-28 | 宜春学院 | 一种中药品质自动检测系统 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111222543A (zh) * | 2019-12-19 | 2020-06-02 | 津海威视技术(天津)有限公司 | 物质识别方法和设备及计算机可读存储介质 |
CN111222543B (zh) * | 2019-12-19 | 2024-03-22 | 津海威视技术(天津)有限公司 | 物质识别方法和设备及计算机可读存储介质 |
CN112945920A (zh) * | 2021-02-01 | 2021-06-11 | 长江大学 | 一种绿茶等级鉴别方法、装置及计算机可读存储介质 |
CN112945920B (zh) * | 2021-02-01 | 2023-11-28 | 长江大学 | 一种绿茶等级鉴别方法、装置及计算机可读存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11282323B2 (en) | Devices, systems, and methods for optical validation | |
KR101297702B1 (ko) | 적분된 투과 및 반사 스펙트럼 응답을 이용한 개선된 위조 화폐 검출기 | |
CA2559102C (en) | Improved fake currency detector using visual and reflective spectral response | |
US6621916B1 (en) | Method and apparatus for determining document authenticity | |
CN101147175B (zh) | 关于钞票验证的改进 | |
CN110503003A (zh) | 基于led阵列与卷积神经网络的茶叶品种鉴定装置与方法 | |
CN104977258B (zh) | 基于二维相关光谱的茶叶/化妆品等品质检测方法 | |
US20100149531A1 (en) | Apparatus and method for object authentication using taggant material | |
CZ299020B6 (cs) | Zpusob, zarízení a bezpecnostní systém, vše k overování pravosti znackování | |
CN104794800B (zh) | 一种基于光滤波的人民币光变油墨识别方法及装置 | |
CN101692052A (zh) | 基于超光谱图像技术的名优茶真伪鉴别方法及装置 | |
KR20060105427A (ko) | 아이템의 진위 판정 방법 및 장치 | |
CN206132659U (zh) | 用于检测柑橘表皮缺陷的在线图像采集系统 | |
CN103903326B (zh) | 多光谱伪钞检测系统及其检测方法 | |
CN103175815A (zh) | 多波长led诱导荧光的茶叶品质无损检测方法及装置 | |
US10664713B2 (en) | Method and system for acquiring and analysing fingerprints with fraud detection | |
EP3903097A1 (en) | Proximity triggered textile substrate classification apparatus and procedure | |
JPS5933220B2 (ja) | 螢光物質を識別する方法 | |
US6658138B1 (en) | Produce texture data collecting apparatus and method | |
CN108027278A (zh) | 照明检测装置及其方法 | |
JP2000346797A (ja) | 茶葉の品質測定方法、品質測定装置及び格付処理装置 | |
RU2460140C1 (ru) | Способ идентификации защитной метки, содержащей люминофор | |
US8203457B1 (en) | Portable currency reader and process for detection of currency value | |
WO2000041143A1 (en) | Coin discriminating device and method | |
CN205541079U (zh) | 基于led荧光光谱和反射光谱的验钞设备 |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20191126 |
|
RJ01 | Rejection of invention patent application after publication |