CN110399904B - 一种基于全彩荧光印品的发射光谱重建方法 - Google Patents
一种基于全彩荧光印品的发射光谱重建方法 Download PDFInfo
- Publication number
- CN110399904B CN110399904B CN201910577043.9A CN201910577043A CN110399904B CN 110399904 B CN110399904 B CN 110399904B CN 201910577043 A CN201910577043 A CN 201910577043A CN 110399904 B CN110399904 B CN 110399904B
- Authority
- CN
- China
- Prior art keywords
- emission spectrum
- color
- matrix
- full
- fluorescent
- 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.)
- Active
Links
- 238000000295 emission spectrum Methods 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000007639 printing Methods 0.000 title claims abstract description 16
- 239000011159 matrix material Substances 0.000 claims abstract description 44
- 238000006243 chemical reaction Methods 0.000 claims abstract description 26
- 239000000463 material Substances 0.000 claims abstract description 11
- 238000012216 screening Methods 0.000 claims abstract description 7
- 230000005284 excitation Effects 0.000 claims abstract description 6
- 230000003595 spectral effect Effects 0.000 claims description 18
- 238000004364 calculation method Methods 0.000 claims description 15
- 238000012360 testing method Methods 0.000 claims description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 7
- 239000004814 polyurethane Substances 0.000 claims description 5
- 229920002635 polyurethane Polymers 0.000 claims description 5
- 229910052761 rare earth metal Inorganic materials 0.000 claims description 4
- 150000002910 rare earth metals Chemical class 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 3
- 238000007646 gravure printing Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 2
- 238000000513 principal component analysis Methods 0.000 claims description 2
- 238000002073 fluorescence micrograph Methods 0.000 abstract description 3
- 238000002360 preparation method Methods 0.000 abstract description 3
- 239000003086 colorant Substances 0.000 abstract 1
- 238000004020 luminiscence type Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- VHUUQVKOLVNVRT-UHFFFAOYSA-N Ammonium hydroxide Chemical compound [NH4+].[OH-] VHUUQVKOLVNVRT-UHFFFAOYSA-N 0.000 description 1
- 235000011114 ammonium hydroxide Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 239000013530 defoamer Substances 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000002189 fluorescence spectrum Methods 0.000 description 1
- 239000002655 kraft paper Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
Images
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/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
-
- 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
-
- 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/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2135—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on approximation criteria, e.g. principal component analysis
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Theoretical Computer Science (AREA)
- Analytical Chemistry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
本发明公开了一种基于全彩荧光印品的发射光谱重建方法。包括制备颜色样本，颜色样本涉及不同的网点覆盖率0%，25%，75%到100%的红、绿、蓝以及三者之间相互叠加形成的黄、品、青等336个，选择不同的承印材料和不同的加网方式处理；测量颜色样本在不同的激发光功率下的发射光谱强度；计算协方差矩阵E;计算协方差矩阵的特征值和特征向量，并将特征值从大到小进行排列；确定前m个主成分及相应的主成分样本矩阵Z；计算主成分的组成系数U，转换矩阵C；根据各自发光强度均值,主成分组成系数和转换矩阵重建发射光谱曲线。本发明的操作简单，可以准确表达全彩色荧光图像的颜色特征，评价图像色域与质量。
Description
技术领域
本发明涉及防伪印刷技术领域，更具体地，涉及一种基于全彩荧光印品的发射光谱重建方法。
背景技术
发射光谱是指用荧光油墨印刷出的样本在红外光激发下发射出来的可见光波长范围内的光强度曲线图谱，是表征无色荧光印刷图案显色的途径之一。此外，近红外光激发比起下转换发光材料，具有较低的声子能量，从而可以获得较高的发光效率，将其优越的发光性能和印刷相结合，用于防伪图案的印刷，可以加大防伪力度。同时利用不同于传统印刷的色光加色法原理，还可以用来生成保护品牌产品的独特彩色图像。为了实现整个技术流程，和光谱预测模型的建立相同，需要对荧光发光样本建立一种光谱重建方法，通过此种方法，可以预先根据想要得到的颜色重建出其发射光谱曲线。实现光谱重建，是准确表达全彩色荧光图像颜色特征的前提，也是评价图像色域与质量的标准，为指导进一步的印刷工作提供保证，所以荧光光谱重建研究受到广泛的关注。
随着假冒伪劣产品的日益渐涨，防伪产业在政府的大力支持下也得到了飞速发展，各种各样的防伪技术和手段层出不穷。隐形防伪油墨因防伪力度高，制备流程简单在众多防伪措施中脱颖而出，而上转换荧光防伪油墨与下转换相比有优越的发光性能，制备过程中因采用水性聚氨酯可以同时达到环保的目的，被作为防伪印刷领域研究热点之一。所以，对于其相应印刷工艺的研究也尤为重要，因发光原理不同，理应研究一种基于全彩荧光印品的发射光谱重建的方法。
发明内容
本发明的目的在于针对现有技术中的不足，提供了一种基于全彩荧光印品的发射光谱重建方法，该方法流程简便，成本低，可以准确表达全彩色荧光图像颜色特征。
本发明的目的通过以下技术方案实现：
一种基于全彩荧光印品的发射光谱重建方法，包括以下步骤：
S1.利用凹版印刷机和上转换荧光墨水印刷制备颜色样本；
S2.对制备完成的颜色样本进行发射光谱强度的测量；
S3.对所测量的荧光样本在351～748nm波长下的发射光谱强度均值化为E_0(λ)；
S4.对每个荧光样本计算其协方差矩阵E；
S5.计算各自协方差矩阵的特征值和特征向量，并将结果按从大到小顺序排列；
S6.对荧光样本的光谱强度值确定前m个主成分并计算对应的主成分样本矩阵Z；
S7.计算上述主成分样本矩阵的组成系数矩阵U和转换矩阵C；
S8.利用所计算出的系数矩阵、均值矩阵和转换矩阵对随机抽取的测试样本的发射光光谱值进行重建，得到重建后的发射光谱曲线。
进一步地，步骤S1中所述上转换荧光墨水中包括水性聚氨酯和稀土上转换发光材料。
本发明采用水性聚氨酯和稀土上转换发光材料作为原材料，制作的上转换荧光墨水的发光效果优于紫外发光即下转换发光材料，且更加难以复制，提高了印品的防伪性能。
进一步地，步骤S1中所述颜色样本的数量为336个。
优选地，所述颜色样本包括网点面积率分别为0％、25％、75％、100％的红、绿、蓝以及四种网点面积相互叠加后的不同明度的黄、品红、青；及其在不同的加网方式下和不同的承印材料上的荧光样本。
进一步地，步骤S2中所述对发射光谱强度的测量过程如下：利用F-4500FLSpectrophotometer荧光密度仪测量在不同的激发光功率下(1.1-1.6每间隔0.05)的光谱强度值。
进一步地，所述协方差矩阵的计算公式为：
E＝[eij]
式中的N为样本总数量336，eij即为发射光谱强度对应值，Sj为样本矩阵的标准差，其计算公式为：
进一步地，所述协方差矩阵E的特征值λ1，λ2，...，λp和特征向量在MATLAB环境下用eig函数计算，并按大小顺序排列以找出主成分个数：
进一步地，所述主成分样本矩阵的确立，计算公式为：
进一步地，所述主成分系数矩阵U通过主成分分析函数在MATLAB中进行计算，转换矩阵C的计算方法是：
U(λ)C≈E(λ)-E0(λ)
UT(λ)U(λ)C≈UT(λ)[E(λ)-E0(λ)]
C＝UT(λ)[E(λ)-E0(λ)]
进一步地，所述重建测试样本的发射光谱的计算公式是：
E(λ)≈E0(λ)+C1U1(λ)+C2U2(λ)+C3U3(λ)+C4U4(λ)+…
与现有技术相比，本发明的有益效果如下：
本发明的光谱重建方法可以预先根据想要得到的颜色重建出相对应的荧光发射光谱曲线，全程都在matlab环境下实现，所以操作简便、成本低。
本发明的计算流程简单，可以高效实现光谱重建，准确表达全彩色荧光图像的颜色特征，评价图像色域与质量。
附图说明
图1是本发明一种基于全彩荧光印品的发射光谱重建方法的流程示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对发明实施方式作进一步描述。
实施例1
如图1所示，本实施例提供了一种基于全彩荧光印品的发射光谱重建方法，包括以下具体步骤：
S1.利用凹版印刷机和上转换荧光墨水印刷提前设计好的颜色样本，颜色样本总数量为336个，包括网点面积率分别为0％、25％、75％、100％的红、绿、蓝以及四种网点面积相互叠加后的不同明度的黄、品红、青；此外，还包括其在不同的加网方式下和不同的承印材料上的荧光样本；
本实施例中，上转换荧光墨水由以下组分和重量配比组成：稀土上转换发光材料：20g、水性聚氨酯：50g、水：20g、氨水：0.6μL、水性消泡剂：0.1g、乙醇：1.5ml；
本实施例中，加网方式包括调幅加网、局部并列和点并点加网处理；承印材料包括白卡纸和牛皮纸。
S2.利用F-4500FL Spectrophotometer荧光密度仪，对制备完成的颜色样本在不同的激发光功率下(1.1-1.6每间隔0.05)的光谱强度值进行测量，对发射光谱强度的测量。
S3.对所测量的荧光样本在351～748nm波长下的发射光谱强度均值化为E0(λ)；
S4.对每个荧光样本计算其协方差矩阵E；
协方差矩阵以及光谱强度均值的计算公式为：
E＝[eij]
式中的N为样本总数量336，eij即为发射光谱强度对应值，Sj为样本矩阵的标准差，其计算公式为：
S5.计算各自协方差矩阵的特征值和特征向量，并将结果按从大到小顺序排列；特征值λ1，λ2，…，λp和特征向量在MATLAB环境下用eig函数计算，并按大小顺序排列以找出主成分个数：
S6.对336个荧光样本的光谱强度值确定前m个主成分并计算对应的主成分样本矩阵Z；计算公式为：
S7.计算上述主成分样本矩阵的组成系数矩阵U和转换矩阵C；转换矩阵C的计算方法是：
U(λ)C≈E(λ)-E0(λ)
UT(λ)U(λ)C≈UT(λ)[E(λ)-E0(λ)]
C＝UT(λ)[E(λ)-E0(λ)]
S8.利用所计算出的系数矩阵U、均值矩阵E0(λ)和转换矩阵C对随机抽取的测试样本的发射光光谱值进行重建，得到重建后的发射光谱曲线。重建测试样本的发射光谱的计算公式是：
E(λ)≈E0(λ)+C1U1(λ)+C2U2(λ)+C3U3(λ)+C4U4(λ)+…
本发明提供的光谱重建方法流程简便、通过对训练样本的发射光谱值进行统计分析，再随机选取检测样本对整个计算流程进行检验，结果证明，该方法可以有效地运用于荧光印品的光谱重建，而且因为波长间隔小，还有效提高了重建效果。
以上所述仅为本发明的较佳实施例，并不用做限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种基于全彩荧光印品的发射光谱重建方法，其特征在于，包括以下步骤：
S1.利用凹版印刷机和上转换荧光墨水印刷制备颜色样本；
S2.对制备完成的颜色样本进行发射光谱强度的测量；
S3.对所测量的荧光样本在351～748nm波长下的发射光谱强度均值化为E0(λ)；
S4.对每个荧光样本计算其协方差矩阵E；
S5.计算各自协方差矩阵的特征值和特征向量，并将结果按从大到小顺序排列；
S6.对荧光样本的光谱强度值确定前m个主成分并计算对应的主成分样本矩阵Z；
S7.计算上述主成分样本矩阵的组成系数矩阵U和转换矩阵C；
S8.利用所计算出的系数矩阵、均值矩阵和转换矩阵对随机抽取的测试样本的发射光光谱值进行重建，得到重建后的发射光谱曲线。
2.根据权利要求1所述的基于全彩荧光印品的发射光谱重建方法，其特征在于，步骤S1中所述上转换荧光墨水中包括水性聚氨酯和稀土上转换发光材料。
3.根据权利要求1所述的基于全彩荧光印品的发射光谱重建方法，其特征在于，步骤S1中所述颜色样本的数量为336个。
4.根据权利要求3所述的基于全彩荧光印品的发射光谱重建方法，其特征在于，所述颜色样本包括网点面积率分别为0％、25％、75％、100％的红、绿、蓝以及四种网点面积相互叠加后的不同明度的黄、品红、青；及其在不同的加网方式下和不同的承印材料上的荧光样本。
5.根据权利要求1所述的基于全彩荧光印品的发射光谱重建方法，其特征在于，步骤S2中所述对发射光谱强度的测量过程如下：利用F-4500FL Spectrophotometer荧光密度仪测量在1.1-1.6每间隔0.05的激发光功率下的光谱强度值。
6.根据权利要求1所述的基于全彩荧光印品的发射光谱重建方法，其特征在于，步骤S4中所述协方差矩阵以及光谱强度均值的计算公式为：
E＝[eij]
式中的N为样本总数量336，eij即为发射光谱强度对应值，Sj为样本矩阵的标准差，其计算公式为：
7.根据权利要求6所述的基于全彩荧光印品的发射光谱重建方法，其特征在于，所述协方差矩阵E的特征值λ1，λ2，…，λp和特征向量在MATLAB环境下用eig函数计算，并按大小顺序排列以找出主成分个数：
8.根据权利要求7所述的基于全彩荧光印品的发射光谱重建方法，其特征在于，所述主成分样本矩阵的确立，计算公式为：
9.根据权利要求8所述的基于全彩荧光印品的发射光谱重建方法，其特征在于，所述主成分系数矩阵U通过主成分分析函数在MATLAB中进行计算，转换矩阵C的计算方法是：
U(λ)C≈E(λ)-E0(λ)
UT(λ)U(λ)C≈UT(λ)[E(λ)-E0(λ)]
C＝UT(λ)[E(λ)-E0(λ)]。
10.根据权利要求9所述的基于全彩荧光印品的发射光谱重建方法，其特征在于，所述重建测试样本的发射光谱的计算公式是：
E(λ)≈E0(λ)+C1U1(λ)+C2U2(λ)+C3U3(λ)+C4U4(λ)+…。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910577043.9A CN110399904B (zh) | 2019-06-28 | 2019-06-28 | 一种基于全彩荧光印品的发射光谱重建方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910577043.9A CN110399904B (zh) | 2019-06-28 | 2019-06-28 | 一种基于全彩荧光印品的发射光谱重建方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110399904A CN110399904A (zh) | 2019-11-01 |
| CN110399904B true CN110399904B (zh) | 2023-06-20 |
Family
ID=68323923
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910577043.9A Active CN110399904B (zh) | 2019-06-28 | 2019-06-28 | 一种基于全彩荧光印品的发射光谱重建方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110399904B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112432934B (zh) * | 2020-11-05 | 2021-07-06 | 北京中科生仪科技有限公司 | 发射光检测方法 |
| CN112611732B (zh) * | 2020-12-21 | 2024-01-26 | 中国电力科学研究院有限公司 | 复合绝缘子用硅橡胶的中红外光谱图库建立方法及系统 |
| CN115046956B (zh) * | 2022-08-17 | 2022-12-13 | 扬州金韵乐器御工坊有限公司 | 基于光学信息的古筝面板材料检测方法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6564714B2 (en) * | 2000-12-06 | 2003-05-20 | Delaware Capital Formation, Inc. | Spectral color control method |
| CN106895916B (zh) * | 2017-01-09 | 2018-10-30 | 浙江大学 | 一种单次曝光拍摄获取多光谱图像的方法 |
| CN108020519B (zh) * | 2017-12-11 | 2020-03-10 | 齐鲁工业大学 | 一种基于颜色恒常性的虚拟多光源光谱重建方法 |
-
2019
- 2019-06-28 CN CN201910577043.9A patent/CN110399904B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN110399904A (zh) | 2019-11-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110399904B (zh) | 一种基于全彩荧光印品的发射光谱重建方法 | |
| Andres et al. | A new anti‐counterfeiting feature relying on invisible luminescent full color images printed with lanthanide‐based inks | |
| CN105835553A (zh) | 隐形防伪标识的制备方法 | |
| CA2522601A1 (en) | Print ink containing a plurality of fluorescent coloring materials and inkjet recording method | |
| JPH01152177A (ja) | インクジェットプリンタ用インク | |
| Javoršek et al. | Colour management in digital textile printing | |
| WO2016084532A1 (ja) | インク画像物生成方法 | |
| CN108501556A (zh) | 一种薄膜凹版印刷专色配色方法 | |
| CN102514408A (zh) | 使用无色荧光油墨印刷彩色防伪图像的方法 | |
| CN104316469B (zh) | 一种面向同色异谱防伪的打印印刷光谱分色方法 | |
| Aydemir et al. | Effects of ink consumption on print quality on coated cellulose-based paper surfaces | |
| WO2014152143A1 (en) | 6-color set plus achromatic(s) for subtractive color combinations | |
| CN106671594B (zh) | 一种确定印刷三原色密度的方法 | |
| US6814426B2 (en) | Color ink-jet printer with dye-based black and pigment-based color ink | |
| CN113776667A (zh) | 一种基于油墨lab值的胶印专色墨调配方法 | |
| Shang et al. | Emission spectral reconstruction based on principal component analysis applied to fluorescence full-color prints | |
| Huang et al. | Preparation and Application of Fluorescent Inkjet Ink | |
| Panák et al. | Producing two-component thermochromic pattern by means of offset printing | |
| König et al. | Assessing the lightfastness of prints by image chrominance histogram quantification | |
| JP2010059285A (ja) | オフセットインキ印刷方法 | |
| Liu | The Influence of UV Light in Color Measurement | |
| TW202024585A (zh) | 校驗色彩之方法 | |
| Možina et al. | Permanence of inkjet prints in relation to typographic and colorimetric characteristics | |
| CN103538385B (zh) | 基于胞元式Neugebauer光谱模型预测超墨量限制多墨样本光谱的方法 | |
| Andersson et al. | Color measurements on prints containing fluorescent whitening agents |
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 |