基于光谱色差权重函数的多光谱降维方法研究

引用本文

曹前, 麻祥才, 白春燕, 苏娜, 崔庆斌. 基于光谱色差权重函数的多光谱降维方法研究[J]. 光谱学与光谱分析, 2023,43(9): 2679-2686.
CAO Qian, MA Xiang-cai, BAI Chun-yan, SU Na, CUI Qing-bin. Research on Multispectral Dimension Reduction Method Based on Weight Function Composed of Spectral Color Difference[J]. Spectroscopy and Spectral Analysis, 2023,43(9): 2679-2686.
Doi:10.3964/j.issn.1000-0593(2023)09-2679-08 复制到剪切板

Permissions

《光谱学与光谱分析》期刊社所有

基于光谱色差权重函数的多光谱降维方法研究

曹前, 麻祥才, 白春燕, 苏娜, 崔庆斌

上海出版印刷高等专科学校印刷包装工程系, 上海 200093

作者简介: 曹前, 1982年生, 上海出版印刷高等专科学校讲师 e-mail: caoqian8889@163.com

收稿日期: 2022-03-01 修回日期: 2022-11-11 接受日期: 2022-11-11

基金: 上海市自然科学基金项目(21ZR1422100)资助

摘要

光谱反射率维度高, 与光照和观察条件无关, 能够真实、客观地描述物体的颜色信息, 由物体本身特性决定, 因此被称为物体的“指纹”。但是, 光谱反射率数据量超过传统三色系统十倍以上, 这些巨大的光谱数据在存储、数据处理及传递等方面造成巨大的负担, 花费太多的计算时间。如果高维光谱可以通过数学变换方法映射到低维空间, 并确保低维空间能够更好地表示原始光谱所覆盖的信息, 可以有效地压缩多光谱数据, 提高基于光谱的颜色复制的处理效率。针对主成分分析法平等地对待可见光范围的所有波长, 重建光谱仅仅是对原始光谱的数学逼近, 由于波长对颜色的重要性不同, 经常会导致光谱重建误差较小而色度误差较大的问题, 提出了一种基于光谱色差权重函数的多光谱降维算法。使用主成分分析法将孟塞尔颜色系统Munsell维度降低到1维, 再恢复重建到31维, Munsell的原始光谱和重建光谱的平均光谱色差作为权重函数。以NCS为训练样本, 分别以NCS、 Munsell和3张多光谱图像为测试样本, 分析和比较本文推荐的权重主成分分析法与主成分分析法以及另外4种权重主成分的性能。以D65/2°和A/2°照明观察条件下的CIELAB色差和均方根误差(RMSE)分别评价测试样本的原始光谱和重建光谱之间的色度重建精度和光谱重建精度。实验统计结果表明: 相对于主成分分析法, 无论测试样本是多光谱数据还是多光谱图像, 推荐的方法在牺牲一定光谱重建精度的情况下, 在D65/2°和在A/2°两种照明观察条件下的色度重建精度得到显著的提高, 而色度重建精度提高对于目前广受关注的基于光谱的颜色复制研究具有非常重要意义。实验统计结果也表明本文推荐的方法的色度重建精度优于目前已经存在的另外4种权重主成分分析法。

关键词: 多光谱降维; 多光谱压缩; 主成分分析法; 权重函数; 光谱色差

中图分类号:O432.3 文献标志码:A

Research on Multispectral Dimension Reduction Method Based on Weight Function Composed of Spectral Color Difference

CAO Qian, MA Xiang-cai, BAI Chun-yan, SU Na, CUI Qing-bin

Department of Printing and Packaging Engineering, Shanghai Publishing and Printing College, Shanghai 200093, China

Abstract

The reflectance spectrum has a high dimension and has nothing to do with illumination and observation. It can truly and objectively describe the color information of an object. The characteristics of the object itself determine it, so it is called the "fingerprint" of the object. However, the amount of reflectance spectra data is more than ten times that of the traditional three-color system, and these huge spectral data cause huge burdens in terms of storage, data processing, and data transfer, and they spend too much computing time. Suppose the high-dimensional spectrum can be mapped to the low-dimensional space through mathematical transformation methods, and ensure that the low-dimensional space data can better represent the information covered by the original spectrum. In that case, the multi-spectral data can be effectively compressed and the processing efficiency of spectrum-based color reproduction can be improved. PCA treats all wavelengths in the visible range equally, and the reconstructed spectrum is only a mathematical approximation to the original spectrum, which often leads to the problem of small spectral reconstruction error and significant colorimetric reconstruction error. A multispectral dimension reduction algorithm based on the weight function of spectral color difference is proposed in this paper. The dimensionality of Munsell was reduced to one dimension by PCA and then restored to 31 dimensions, and the average spectral color difference between Munsell's original spectrum and its reconstructed spectrum was used as a weight function. Taking NCS as the training sample and NCS, Munsell and 3 multispectral images as the test samples respectively, the performance of the proposed method of this paper and the classical PCA and the other four weighted PCA are analyzed and compared. CIELAB color difference under the conditions of multiple Lighting and viewing (D65/2° and A/2°) and root mean square error (RMSE) evaluate the colorimetric and spectral reconstruction accuracy between the original spectra and the reconstructed spectra of the test sample respectively. The experimental results show that: compared with PCA, the proposed method has greatly improved colorimetric reconstruction accuracy at the expense of a small amount of spectral reconstruction accuracy. The improvement of colorimetric reconstruction accuracy is very important for spectral color reproduction. The results also show that the colorimetric reconstruction accuracy of the proposed method is better than that of the other four existing weighted PCA.

Keyword: Multispectral dimensionality reduction; Multispectral compression; Principal component analysis; Weight function; Spectral color difference

文章图片

引言

光谱反射率维度高, 与光照和观察条件无关, 能够客观、真实地描述物体的颜色信息, 是物体本身特性决定的, 因此被称为物体的“ 指纹” , 在光谱颜色复制^[1]、纺织^[2]等领域受到广泛关注。但是, 用于光谱颜色复制领域的光谱数据通常选取10 nm间隔的400~700 nm波段的共计31维数据组成。其数据量超过传统三色颜色系统十倍以上, 这些巨大的光谱数据在存储、数据处理和数据传递等方面造成巨大的负担, 花费太多的计算时间。如果高维光谱可以通过数学变换方法映射到低维空间, 并确保低维空间数据能够更好地表示原始光谱所覆盖的信息, 可以有效地压缩多光谱数据, 提高多光谱色彩再现的处理效率。因此, 光谱数据压缩或者光谱数据维度降低成为光谱颜色复制的一项关键技术。已经有大量数据维度降低算法应用到高维光谱降维, 如主成分分析(principal component analysis, PCA)、压缩感知算法^[3]、 box-cox变换^[4], 神经网络^[5]等。物体表面的反射率光谱通常是由可见光范围内波长的光滑函数组成, 这些光谱强烈相关并可以由少量的正交基函数组成, 主成分分析法正是基于这个事实来实现多光谱数据降维。因此, 基于PCA的多光谱数据压缩算法引起了广泛关注, 国内外开展了大量相关研究。但是, 主成分分析法平等地对待可见光范围的所有波长, 重建光谱仅仅是对原始光谱的数学逼近, 由于根据波长对颜色的重要性不同, 经常会导致光谱误差较小而颜色色差较大。为了克服这个缺点, 主成分分析法应该根据波长对颜色的重要性进行修改, 更确切地说, 如果光谱数据的某些部分对颜色值不太重要, 就可以增加相应的权重来减少这部分光谱对颜色影响。已经有一些研究者提出权重的主成分分析法来进行光谱数据降维, 具体来说, 就是在形成相关矩阵和计算特征值之前将光谱数据乘以适当的系数或者权重函数。为了尽可能减少光谱维度降低过程中颜色信息的损失, 不同的学者从理论推导和大量实验基础上, 提出了基于颜色匹配函数、视锥细胞响应函数和光谱变化引起的色差等方面提出了不同的权重函数。 Laamanen等^[6]于2008年提出基于颜色匹配函数的权重函数的权重主成分多光谱降维算法。王莹等^[7]于2011年提出一个基于代表人眼平均视觉特性的颜色匹配函数的权重函数。 Tian等^[8]于2013年建议基于颜色匹配函数的不同组合的多光谱降维算法的权重函数; Wu等^[9]于2015年提议基于人眼视觉特性的颜色匹配函数之和的平方根作为多光谱降维算法的权重函数; Cao等^[10]于2016年也提出颜色匹配函数的平方根之和作为多光谱降维算法的权重函数; 刘士伟等^[11]于2017年构建基于视锥细胞响应函数不同组合的权重函数。梁金星等^[12]提出基于视觉特征函数的权重函数; Ma等^[13]于2020年构建了另外一个基于色差变化的权重函数用于保留更多的颜色信息。

上述学者提出的权重函数在确保光谱精度可以接受的条件下, 显著地提高了光谱降维和压缩过程的色度精度, 保留了更多的颜色信息, 对于多光谱色彩复制有非常积极意义。可是, 上面提到的权重主成分分析法, 当改变测试样本或是使用其他光谱图像数据时, 其色度误差仍然不可忽视, 对于在光谱颜色复制领域中的实际应用还有很大的距离。在多光谱图像通过降低维度进行压缩过程中寻找一个保留更多颜色信息的权重函数对于光谱颜色再现具有重要现实意义。

本文提出了一种基于光谱色差的权重函数的多光谱降维算法, 权重函数由光谱数据集原始光谱和该光谱数据集重建光谱之间的平均光谱色差构成。本研究以NCS为训练样本, 分别以NCS、 Munsell和3张多光谱图像测试和对比本文推荐的方法和主成分分析法以及另外4种权重主成分在多光谱降维过程中的性能; 利用公式将测试样本压缩到低维光谱空间, 再利用公式将低维光谱数据恢复还原到高维光谱。以CIELAB色差和均方根误差(RMSE)来评价测试样本的原始光谱和恢复光谱之间的色度误差和光谱误差。

1 实验部分

1.1 权重主成分分析

一组光谱数据矩阵, 设采样维度为n, 样本数为p, 用矩阵S表示这个光谱数据集。定义对角矩阵W, 其对角线元素由权重函数w(λ )组成。光谱矩阵S乘以矩阵W, 权重光谱矩阵S_W由式(1)表示

$S_{W} = WS (1)$

权重光谱矩阵S_W的相关系数矩阵Σ _W如式(2)所示

$Σ_{W} = S_{W} {S^{T}}_{W} / p = (WS) {(WS)}^{T} / p = WΣ W^{T} (2)$

如果相关系数矩阵Σ _W的特征值按照从大到小的顺序, 前m特征值对应的特征向量组成矩阵 ${V^{T}}_{W}$ , 权重光谱在低维空间的投影值可以用式(3)表示。

$C_{W} = {V^{T}}_{W} (S_{W} - {\overset{̅}{S}}_{W}) (3)$

式(3)中, ${\overset{̅}{S}}_{W}$ 表示权重光谱矩阵S_W在维度方向的平均值。重建的光谱由式(4)计算得到

$\dot{S} = inv (W) ({\dot{V}}_{W} C_{W} + {\overset{̅}{S}}_{W}) (4)$

式(4)中, inv(W)表示方阵W的逆矩阵。

1.2 光谱色差公式

根据文献[14, 15], 本文定义光谱色差如下: 任意一个光谱反射率在可见光范围内各个波段的色差, 其计算公式如式(5)所示。式中, L(λ )、 a(λ )和b(λ )的计算公式如式(6)、式(7)和式(8)所示。

$\begin{array}{l} Δ E (λ) = \\ \sqrt[]{(L_{1} (λ) - L_{2} {(λ))}^{2} + (a_{1} (λ) - a_{2} {(λ))}^{2} + (b_{1} (λ) - b_{2} {(λ))}^{2}} \end{array}$ (5)

$\begin{array}{l} L (λ) = 116 kS (λ) \overset{̅}{y} (λ) f (\frac{Y}{Y_{n}}) (6) \\ a (λ) = 500 kS (λ) [\overset{̅}{x} (λ) f (\frac{X}{X_{n}}) - \overset{̅}{z} (λ) f (\frac{Z}{Z_{n}})] (7) \\ b (λ) = 200 kS (λ) [\overset{̅}{z} (λ) f (\frac{Z}{Z_{n}}) - \overset{̅}{y} (λ) f (\frac{Y}{Y_{n}})] (8) \end{array}$

其中,

$f (\frac{u}{u_{n}}) = \{\begin{array}{l} 7.787 (\frac{u}{u_{n}}), & \frac{u}{u_{n}} \leq 0.008856 \\ \frac{1}{3 u} {(\frac{u}{u_{n}})}^{\frac{1}{3}}, & \frac{u}{u_{n}} > 0.008856 \end{array} (9)$

其中, u依次替换为物体的三刺激值X、 Y、 Z, u_n依次替换为参考光源的三刺激值X_n、 Y_n和Z_n, S(λ )表示光源的相对光谱功率分布, k表示使得参考光源的Y_n等于100时的系数。

1.3 权重函数的构造

1.3.1 权重函数构成思路和方法

定义光谱数据集R的采样维度是n, 样本个数为p, 由光谱反射率r₁, r₂, …, r_p组成, 如式(10), 采用主成分分析对其降低维度到m (m< n)维, 再将低维数据重建到n维光谱数据R', 如式(11)所示。

$\begin{array}{l} R = [r_{1} r_{2} \dots r_{p}] (10) \\ R' = [r'_{1} r'_{2} \dots r'_{p}] (11) \end{array}$

光谱反射率r₁的L₁(λ ), a₁(λ )和b₁(λ )是分别使用式(6)、式(7)和式(8)计算得到, 光谱反射率r'₁的L'₁(λ ), a'₁(λ )和b'₁(λ )是分别使用式(6)、式(7)和式(8)计算得到。光谱反射率r₁和r'₁的光谱色差Δ E₁(λ )就可以使用式(5)计算。类似的, 可以计算光谱反射率r₂和r'₂的光谱色差Δ E₂(λ ), …, 光谱反射率r_p和r'_p的光谱色差Δ E_p(λ )。最后, 使用Δ E₁(λ ), Δ E₂(λ ), …, Δ E_p(λ )的平均值作为本文推荐的权重函数。

1.3.2 权重函数的计算

通过大量的测试, 本文选取有1 600个光谱组成的Munsell数据集作为R, 使用主成分分析法降低维度到1维, 再恢复到31维, 得到光谱数据集R', 全色差公式计算使用的参考光源为等能白光E, 根据权重函数构成思路和方法, 经过归一化后, 权重函数值如表1所示。为了便于表达, 本文提出的权重函数命名为SCDWF, 对应的权重主成分命名为SCDPCA。

表1 由Munsell光谱数据集的原始光谱和恢复光谱的平均光谱色差构成的权重函数值(SCDWF) Table 1 The weight function value (SCDWF) composed of the average spectral color difference between the original spectra and the recovered spectra of Munsell

1.4 其他研究者提出的权重主成分多光谱降维算法

为了光谱压缩过程中保留更多的颜色信息, Laamanen等^[6]于2008年提出一个基于颜色匹配函数的权重函数, 命名为WF1, 如式(12)所示, 其对应的权重主成分命名为W1PCA。

$WF 1 = \overset{̅}{x} (λ) + \overset{̅}{y} (λ) + \overset{̅}{z} (λ) + 1 (12)$

式(12)中, $\overset{̅}{x}$ (λ )、 $\overset{̅}{y}$ (λ )和 $\overset{̅}{z}$ (λ )表示颜色匹配函数。

Tian等^[8]于2013年也提议基于人眼视觉特性的颜色匹配函数的不同组合的权重函数, 命名为WF2, 如式(13)所示, 其对应的权重主成分命名为W2PCA。

$WF 2 = \overset{̅}{x} (λ) + \overset{̅}{y} (λ) + \overset{̅}{z} (λ) (13)$

刘士伟等^[11]于2017年构建基于视锥细胞响应函数不同组合的权重函数, 命名为WF3, 如式(14)所示, 其对应的权重主成分命名为W3PCA。

$WF 3 = \sqrt[]{L (λ) + M (λ) + S (λ)} (14)$

式(14)中, L(λ ), M(λ )和S(λ )分别表示感红、感绿和感蓝三种视锥细胞响应函数。

Ma等^[13]于2020年构建了另外一个基于微小波长变化引起的色差的权重函数, 命名为WF4, 其对应的权重主成分命名为W4PCA。本文推荐的权重函数和其他4种权重函数曲线如图1所示。

	Figure Option View Download New Window
	图1 本文推荐的权重函数SCDWF和其他4种权重函数曲线Fig.1 Curves of the recommended weight function (SCDWF) and other four weight functions

1.5 实验数据

本文使用两个著名的显色系统孟塞尔系统Munsell^[16]和自然颜色系统NCS的色卡的光谱数据作为实验数据, 它们分别是有1 600和1 950色卡构成, 这些光谱数据都是用分光光度计测量得到, 它们的光谱反射率和对应的色度坐标如图2和图3所示。来自哥伦比亚大学^[17]的3张多光谱图像也作为本实验的测试样本, 图4显示这3张多光谱图像, 每张多光谱图像的分辨率都是512像素× 512像素, 共计262 144个像素构成。

	Figure Option View Download New Window
	图2 照明观察D65/ 2° 条件下(a)NCS的光谱反射率及其在CIELAB空间(b)中对应的a* b* 图Fig.2 Spectral reflectance of (a) NCS along with their corresponding a* b* plots in CIELAB space (b) under D65 and the CIE 1931 2° Standard Observer

	Figure Option View Download New Window
	图3 照明观察D65/2° 条件下(a)Munsell的光谱反射率及其在CIELAB空间(b)中对应的a* b* 图Fig.3 Spectral reflectance of (a) Munsell along with their corresponding a* b* plots in CIELAB space (b) under D65 and the CIE 1931 2° Standard Observer

	Figure Option View Download New Window
	图4 来自哥伦比亚大学的3张多光谱图像 (a): Image 1; (b): Image 2; (c): Image 3Fig.4 Three multispectral images from Columbia University (a): Image 1; (b): Image 2; (c): Image 3

首先, 使用本文推荐的方法和已经存在的方法分别通过训练样本NCS得到6维特征向量; 其次, 通过式(3)将31维的测试样本NCS、 MUNSELL和3张多光谱图像压缩到6维低维空间, 再通过式(4)将测试样本的6维低维空间数据还原到31维。然后, 使用光谱精度和色度精度评价测试样本的原始光谱数据和压缩后再恢复的光谱数据之间的误差。最后, 比较本文推荐的方法和已经存在的方法性能优劣。

1.6 评价指标

原始光谱和压缩后恢复光谱之间存在误差, 该误差主要通过色度和光谱两个指标来评价。原始光谱和压缩后恢复光谱的光谱精度的色度精度采用1976CIELAB色差评估, 色差小表示色度精度高。原始光谱和压缩后恢复光谱的光谱精度通过均方根误差(root mean squared error, RMSE)^{[18, 19]}来评价, 如式(15)所示。

$RMSE = \sqrt[]{\frac{\sum (S (λ_{i}) - \dot{S} (λ_{i} {))}^{2}}{n}} (15)$

式(15)中, S(λ _i)和 $\dot{S}$ (λ _i)分别表示原始光谱和压缩后恢复光谱在波长λ _i处的光谱反射率, n表示光谱样本数。

2 结果与讨论

对比了我们推荐的方法和一些已经存在的方法的色度精度和光谱精度。以NCS光谱数据集为训练样本, 使用本文推荐的方法和一些已经存在的方法得到6维主成分, 分别以NCS、 Munsell和3种多光谱图像为测试样本, 利用式(3)将测试样本的维度降低到6维低维空间, 再利用式(4)还原恢复到31维光谱数据, 使用式(13)和式(14)评价测试样本的原始光谱数据和恢复光谱之间的色度精度和光谱精度。

2.1 色度精度比较

表2统计了以NCS为训练样本, 测试样本NCS的原始光谱和重建光谱数据在D65/2° 和A/2° 照明观察条件下的平均色差和最大色差。图5表示以NCS为训练样本, 测试样本NCS的原始光谱和重建光谱数据之间色差分布的箱形图(boxplot)。从表2和图5可以看出, 本文推荐的权重主成分在D65/2° 和A/2° 照明观察条件下的色度精度不仅优于主成分分析法, 而且也优于其他4种权重主成分分析法。

表2 以NCS为测试样本时推荐的方法和已经存在一些方法在D65/2° 和A/2° 照明观察条件下色度精度对比 Table 2 Comparison of colorimetric accuracy under the illuminats of D65/2 ° and A/2° between the proposed method and some existing methods when the test sample is NCS

	Figure Option View Download New Window
	图5 测试样本NCS的原始光谱和重建光谱之间的色差分布的箱形图(boxplot) (a): D65/2° ; (b): A/2° Fig.5 Boxplots of the color difference distribution between the original spectra and the recovered spectra of NCS (a): D65/2° ; (b): A/2°

表3统计了以NCS为训练样本, 测试样本Munsell的原始光谱和重建光谱数据在在D65/2° 和A/2° 照明观察条件下的平均色差和最大色差。图6表示以NCS为训练样本, 测试样本Munsell的原始光谱和重建光谱数据之间色差分布的箱形图(boxplot)。从表3和图6可以看出, 以NCS为训练样本, 以Munsell测试样本时的色度精度的统计结果和色差分布情况再一次验证以NCS为训练样本, 以NCS为测试样本时在色度精度方面的结论。

表3 以Munsell为测试样本时推荐的方法和已经存在一些方法在D65/2° 和A/2° 照明观察条件下色度精度对比 Table 3 Comparison of colorimetric accuracy under the illuminats of D65/2 ° and A/2° between the proposed method and some existing methods when the test sample is Munsell

	Figure Option View Download New Window
	图6 测试样本Munsell的原始光谱和重建光谱之间的色差分布的箱形图(boxplot) (a): D65/2° ; (b): A/2° Fig.6 Boxplots of the color difference distribution between the original spectra and the recovered spectra of Munsell (a): D65/2° ; (b): A/2°

2.2 光谱精度比较

表4统计了测试样本NCS和Munsell的原始光谱和恢复光谱之间的平均均方根误差(mean RMSE)和最大均方根误差(max RMSE)。图7表示测试样本NCS和Munsell的原始光谱和恢复光谱之间的均方根误差(RMSE)分布的箱形图(boxplot); 从表4和图7可以看出, 测试样本为NCS或Munsell时, 主成分分析法的光谱精度(RMSE)优于本文推荐的方法和其他4种权重主成分。结合色度精度结论, 权重主成分的色度精度提高是以牺牲一定的光谱精度为代价。但是, 色度精度提高对于光谱颜色复制有重要的现实意义。

表4 测试样本为NCS和Munsell时推荐的方法和已经存在一些方法的光谱精度对比 Table 4 Comparison of spectral accuracy between the proposed method and some existing methods when the test samples are NCS and Munsell

	Figure Option View Download New Window
	图7 测试样本的原始光谱和恢复光谱之间的均方根误差(RMSE)的箱形图(boxplot) (a): NCS; (b): MunsellFig.7 Boxplot of the root mean square error (RMSE) between the original spectra of the test sample and the recovered spectra of the testing samples (a): NCS; (b): Munsell

2.3 测试样本为3张多光谱图像时各种降维方法的色度精度和光谱精度对比

表5统计了以NCS为训练样本, 3张测试样本的原始光谱和重建光谱数据在D65/2° 照明观察条件下的平均色差和最大色差。表6统计了以NCS为训练样本, 3张测试样本的原始光谱和重建光谱数据在A/2° 照明观察条件下的平均色差和最大色差。表7统计了3张测试样本的原始光谱和重建光谱数据之间的平均均方根误差(mean RMSE)和最大均方根误差(max RMSE)。

表5 测试样本为3张多光谱图像时推荐的方法和已经存在一些方法在D65/2° 照明观察条件下色度精度对比 Table 5 Comparison of colorimetric accuracy under the illuminats of D65/2° between the proposed method and some existing methods when the test samples are three multispectral images

表6 测试样本为3张多光谱图像时推荐的方法和已经存在一些方法在A/2° 照明观察条件下色度精度对比 Table 6 Comparison of colorimetric accuracy under the illuminats of A/2° between the proposed method and some existing methods when the test samples are three multispectral images

表7 测试样本为3张多光谱图像时推荐的方法和已经存在一些方法的光谱精度对比 Table 7 Comparison of spectral accuracy between the proposed method and some existing methods when the test samples are three multispectral images

从表5、表6和表7可以看出, 当测试样本为3张多光谱图像时色度精度和光谱精度的统计结果与测试样本为NCS或者Munsell时的结论大致一致, 更进一步验证了本文推荐的权重函数的色度精度不仅优于主成分分析法也优于其他4种权重主成分分析法。

3 结论

提出了基于光谱色差权重函数的多光谱降维方法, 权重函数由光谱数据集和该光谱数据集经由主成分重建光谱的平均光谱色差构成。本研究以NCS为训练样本, 分别以NCS、 Munsell和3张多光谱图像验证本文推荐的方法和主成分分析法以及另外4种权重主成分在多光谱降维过程中的性能优劣; 结果表明: 对比于主成分分析法, 测试样本无论是多光谱数据还是多光谱图像, 本文推荐的方法在牺牲了一定光谱精度的情况下, 在D65/2° 和在A/2° 两种照明观察条件下的色度精度得到显著的提高。实验统计结果也显示本文推荐的方法的色度精度优于目前已经存在的另外4种权重主成分。本文提出的多光谱降维方法保留了更多的颜色信息, 适合于高维光谱的降维、压缩和存储等, 为高维光谱数据到低维连接空间数据的映射方法提供参考, 对于基于光谱的颜色复制研究有很重要的意义。

参考文献

文献列表

[1]	Habib T, Green P, Nussbaum P. In Proc. IS&T Int'. Symp. on Electronic Imaging: Color Imaging XXV: Displaying, In Processing, Hardcopy, and Applications, 2020: 121. [本文引用:1]
[2]	Mohtasham J, Nateri A S, Khalili H. Color. Technol. , 2012, 128(3): 199. [本文引用:1]
[3]	Shen D, Zhang L, Liang D, et al. Laser Phys. , 2017, 27(7): 075201. [本文引用:1]
[4]	Rayat A, Amirshahi S, Agahian F. ColorRes. Appl. , 2014, 39(2): 136. [本文引用:1]
[5]	Hajipour A, Shams Nateri A. Color Res. Appl. , 2017, 42(2): 182. [本文引用:1]
[6]	Laamanen H, Jetsu T, Jaaskelainen T, et al. JOSA A, 2008, 25(6): 1383. [本文引用:2]
[7]	WANG Ying, WANG Zhong-min, WANG Yi-feng, et al(王莹, 王忠民, 王义峰, 等). Optics and Precision Engineering(光学精密工程), 2011, 19(5): 1171. [本文引用:1]
[8]	Tian J, Tang Y. Opt. Lett. , 2013, 38(15): 2818. [本文引用:2]
[9]	Wu G, Liu Z, Fang E, et al. Optik, 2015, 126(11-12): 1249. [本文引用:1]
[10]	Cao Q, Wan X, Li J, et al. Opt. Rev. , 2016, 23(5): 753. [本文引用:1]
[11]	LIU Shi-wei, LIU Zhen, TIAN Quan-hui, et al(刘士伟, 刘真, 田全慧, 等). Spectroscopy and Spectral Analysis(光谱学与光谱分析), 2017, 37(5): 1601. [本文引用:2]
[12]	LIANG Jin-xing, WAN Xiao-xia, LU Wei-pent(梁金星, 万晓霞, 卢玮朋). Spectroscopy and Spectral Analysis(光谱学与光谱分析), 2017, 37(1): 177. [本文引用:1]
[13]	Ma X, Wang, X, Zhang, D. Optik, 2020, 203: 163959. [本文引用:2]
[14]	Viggiano J A S. Advanced Printing of Conference Summaries: SPSE's 43rd Annual Conference, 1990: 222. [本文引用:1]
[15]	Viggiano J A S. 9th Congress of the International Colour Association. International Society for Optics and Photonics, 2002, 4421: 701. [本文引用:1]
[16]	University of Eastern Finland , Spectral Color Research Group: https: //www. uef. fi/spectral/spectral-databas. [本文引用:1]
[17]	Yasuma F, Mitsunaga T, Iso D, et al. IEEE Transactions on Image Processing, 2010, 19(9): 2241. [本文引用:1]
[18]	Imai F H, Berns R S, Tzeng D Y. J. Imaging Sci. Techn. , 2000, 44(4): 280. [本文引用:1]
[19]	Imai F H, Rosen M R, Berns R S. First European Conference on Color in Graphics, Imaging, and Vision. Conference Proceedings, 2002, 2002(1): 492. [本文引用:1]

2020

0.0

... , 在光谱颜色复制^[1]、纺织^[2]等领域受到广泛关注 ...

2012

0.0

... , 在光谱颜色复制^[1]、纺织^[2]等领域受到广泛关注 ...

2017

0.0

... 已经有大量数据维度降低算法应用到高维光谱降维, 如主成分分析(principal component analysis, PCA)、压缩感知算法^[3]、 box-cox变换^[4], 神经网络^[5]等 ...

2014

0.0

... 已经有大量数据维度降低算法应用到高维光谱降维, 如主成分分析(principal component analysis, PCA)、压缩感知算法^[3]、 box-cox变换^[4], 神经网络^[5]等 ...

2017

0.0

... 已经有大量数据维度降低算法应用到高维光谱降维, 如主成分分析(principal component analysis, PCA)、压缩感知算法^[3]、 box-cox变换^[4], 神经网络^[5]等 ...

2008

0.0

... Laamanen等^[6]于2008年提出基于颜色匹配函数的权重函数的权重主成分多光谱降维算法 ...

... 4 其他研究者提出的权重主成分多光谱降维算法为了光谱压缩过程中保留更多的颜色信息, Laamanen等^[6]于2008年提出一个基于颜色匹配函数的权重函数, 命名为WF1, 如式(12)所示, 其对应的权重主成分命名为W1PCA ...

2011

0.0

... 王莹等^[7]于2011年提出一个基于代表人眼平均视觉特性的颜色匹配函数的权重函数 ...

2013

0.0

... Tian等^[8]于2013年建议基于颜色匹配函数的不同组合的多光谱降维算法的权重函数 ...

... Tian等^[8]于2013年也提议基于人眼视觉特性的颜色匹配函数的不同组合的权重函数, 命名为WF2, 如式(13)所示, 其对应的权重主成分命名为W2PCA ...

2015

0.0

... Wu等^[9]于2015年提议基于人眼视觉特性的颜色匹配函数之和的平方根作为多光谱降维算法的权重函数 ...

2016

0.0

... Cao等^[10]于2016年也提出颜色匹配函数的平方根之和作为多光谱降维算法的权重函数 ...

2017

0.0

... 刘士伟等^[11]于2017年构建基于视锥细胞响应函数不同组合的权重函数 ...

... 刘士伟等^[11]于2017年构建基于视锥细胞响应函数不同组合的权重函数, 命名为WF3, 如式(14)所示, 其对应的权重主成分命名为W3PCA ...

2017

0.0

... 梁金星等^[12]提出基于视觉特征函数的权重函数 ...

2020

0.0

... Ma等^[13]于2020年构建了另外一个基于色差变化的权重函数用于保留更多的颜色信息 ...

... Ma等^[13]于2020年构建了另外一个基于微小波长变化引起的色差的权重函数, 命名为WF4, 其对应的权重主成分命名为W4PCA ...

1990

0.0

... 2 光谱色差公式根据文献[14,15], 本文定义光谱色差如下: 任意一个光谱反射率在可见光范围内各个波段的色差, 其计算公式如式(5)所示 ...

2002

0.0

... 2 光谱色差公式根据文献[14,15], 本文定义光谱色差如下: 任意一个光谱反射率在可见光范围内各个波段的色差, 其计算公式如式(5)所示 ...

0.0

... 5 实验数据本文使用两个著名的显色系统孟塞尔系统Munsell^[16]和自然颜色系统NCS的色卡的光谱数据作为实验数据, 它们分别是有1 600和1 950色卡构成, 这些光谱数据都是用分光光度计测量得到, 它们的光谱反射率和对应的色度坐标如图2和图3所示 ...

2010

0.0

... 来自哥伦比亚大学^[17]的3张多光谱图像也作为本实验的测试样本, 图4显示这3张多光谱图像, 每张多光谱图像的分辨率都是512像素#cod#x000D7 ...

2000

0.0

... 原始光谱和压缩后恢复光谱的光谱精度通过均方根误差(root mean squared error, RMSE)^[18,19]来评价, 如式(15)所示 ...

2002

0.0

... 原始光谱和压缩后恢复光谱的光谱精度通过均方根误差(root mean squared error, RMSE)^[18,19]来评价, 如式(15)所示 ...