光谱图像处理方法技术

本发明专利技术为光谱图像处理方法，将人眼三种视锥光谱响应构成一个加权矩阵，与原始光谱图像相乘后进行主成分分析，取前三个特征向量作为新的六维线性模型的前三个向量，接着以残余光谱进行主成分分析，取其前三特征向量作为新的六维线性模型的后三个向量，并用这六个基向量将原始高维光谱图降维处理得到六维多通道图像，再将该六维多通道图像代替原始高维光谱图像用于颜色复制过程；处理光谱图像时巧妙地将人眼视觉特性与光谱线性降维模型相结合，并兼顾光谱和色度两方面的匹配精度，能够还原和复制出高精度、高仿度、最佳色彩匹配和还原效果的光谱图像。

【技术实现步骤摘要】
光谱图像处理方法
本专利技术涉及光谱图像压缩和复制领域，尤其涉及一种采用将人眼视锥响应与光谱图像的压缩相结合的处理方式完成复制的方法。
技术介绍
将光谱图像用于颜色复制的处理方法主要分为两大类，一类是直接将高维的原始图像用于颜色复制，在这个复制过程中需要经过复杂的计算才能将全部光谱信息用于颜色复制，由于计算复杂度高，会需要占用系统内存进行运算，而且整个运算的耗时长。更为重要的是高维光谱图像中很大一部分光谱信息的对于颜色复制来说是烦冗而多余的信息。在直接将高维光谱图像用于颜色复制时，耗费了大量的计算过程在不需要的多余信息上，存在对时间、系统和资源的浪费。另外一类方法是将高维光谱图像进行降维压缩处理，去掉对颜色复制多余的光谱信息后，再用于颜色复制。该方法保存的降维光谱图像存储空间小，运算效率高。但是这一类方法存在的缺陷是采用数学统计方式对光谱图像进行降维压缩，重构之后的光谱图像是对原始光谱图像的数学逼近，而非完全精准复制，在视觉上存在失真，并非色彩上的最佳匹配和还原。在对原始光谱图像进行降维压缩处理时，一般会采用主成分分析、独立成分分析等方法，原始信息数据被统计结果取代，在进行还原的时候就无法完全重建原始信息，而造成仅能重建光谱图像的数学近似结果。
技术实现思路
本专利技术针对以上情况提出了一种将高维光谱图像首先结合人眼视觉特性进行线性降维压缩再用于颜色复制的光谱图像处理方法，其在用于颜色复制时可最大限度保留原始光谱图像色度信息和光谱信息。本专利技术所涉及的用于颜色复制的光谱图像处理方法包括以下步骤：通过人眼视锥响应构建加权矩阵Σ；将该加权矩阵Σ与原始光谱图像R相乘，得到加权后光谱图像G，再对加权光谱图像G进行主成分分析，得到第一组三个特征向量A1/A2/A3；取第一组的三个特征向量Ag＝(A1,A2,A3)作为六维线性模型的前三个基向量，求与这个三基向量对应的系数矩阵Wg，并用这三个基向量重构光谱，得到重构光谱图像通过原始光谱图像R减去三个基向量Ag＝(A1,A2,A3)重构的光谱图像得到残余光谱Rlost；对残余光谱Rlost进行主成分分析，得到第二组三个特征向量B1/B2/B3；取第二组的三个特征向量Blost＝(B1,B2,B3)，作为六维线性模型的后三个基向量，求与特征向量组B1/B2/B3对应的系数矩阵Wlost；利用这六个基向量将原始高维光谱R进行降维处理得到六维多通道图像；六维多通道图像的保存；最后将六维多通道图像代替原始光谱图像用于颜色复制过程。所述加权矩阵Σ，满足通过人眼三种视锥细胞长波L(λ)，中波M(λ)，短波S(λ)对重构光谱图像反射率和原始光谱图像反射率R的视锥响应总的偏差最小情况下，由人眼长波L(λ)、中波M(λ)和短波S(λ)对角化后再进行矩阵相加构成的对角矩阵Σ＝L+M+S。所述重构光谱是通过所得，Wg是由加权后光谱图像G通过最小二乘法求得。所述六维多通道图像由对应系数矩阵Wg和Wlost构成。系数矩阵Wg对应于六维多通道图像的前三个通道数据，系数矩阵Wlost对应于六维多通道图像的后三个通道数据。所述六维多通道图像的系数矩阵Wlost由残余光谱Rlost通过最小二乘法求得。所述六维多通道图像的保存是对系数矩阵Wg和Wlost保存。所述六维多通道图像在用于颜色复制时其对原始高维光谱图像的重构光谱为前三个基向量Ag＝(A1,A2,A3)重构的光谱图像加上后三个基向量Blost＝(B1,B2,B3)重构的残余光谱图像一种用于颜色复制的光谱图像处理方法，其特征在于，该处理方法包括以下步骤：通过人眼视锥响应构建加权矩阵Σ，设ε1，ε2，ε3分别为人眼长波L，人眼中波M，人眼短波S三种视锥细胞对重构光谱图像反射率和原始光谱图像反射率R的视锥响应偏差，而ε为总的视锥响应偏差。则有：ε＝ε1+ε2+ε3(4)上述公式中||||1表示矩阵的1范数。将式(1)、(2)、(3)代入到(4)，则有如果重构光谱与原始光谱在总的视锥响应上具有最小的偏差ε，则可最大程度实现变光源下的颜色匹配。这里，L、M、S是视锥响应灵敏度矢量的对角阵形式，因此定义加权矩阵ΣΣ＝L+M+S(6)矩阵Σ与原始光谱图像R相乘，得到加权后光谱图像G＝ΣR，对相乘得到的结果进行主成分分析，得到第一组三个特征向量Ag＝(A1,A2,A3)；G＝ΣR(7)取第一组的三个特征向量作为六维线性模型的前三个基向量，并用第一组的三个特征向量Ag＝(A1,A2,A3)重构光谱，得到重构光谱图像Wg是对应的系数矩阵，由加权后光谱图像G通过最小二乘法求得。通过原始光谱图像R减去三个基向量重构的光谱图像得到残余光谱Rlost；残余光谱Rlost进行主成分分析，得到第二组三个特征向量B1/B2/B3；第二组的三个特征向量Blost＝(B1,B2,B3)，作为六维线性模型的后三个基向量，得到重构残余光谱将重构光谱与重构残余光谱相加得到原始高维光谱图像R的重构光谱图像本专利技术将人眼三种视锥光谱响应构成一个加权矩阵，与原始光谱图像相乘后进行主成分分析，取前三个特征向量作为新的六维线性模型的前三个向量，接着以残余光谱进行主成分分析，取其前三特征向量作为新的六维线性模型的后三个向量，并用这六个基向量将原始高维光谱图降维处理得到六维多通道图像，再将该六维多通道图像代替原始高维光谱图像用于颜色复制过程；处理光谱图像时巧妙地将人眼视觉特性与光谱线性降维模型相结合，所得到的六维多通道图像兼顾光谱和色度两方面的匹配精度，能够还原和复制出高精度、高仿度、最佳色彩匹配和还原效果的光谱图像。【附图说明】图1是本专利技术用于颜色复制的光谱图像处理方法流程图；【具体实施方式】下面将结合本专利技术附图和具体实施方式对本专利技术用于颜色复制的光谱图像处理方法进行进一步的详细说明。本专利技术所涉及的用于颜色复制的光谱图像处理方法以400nm-700nm的光谱范围并按10nm的均匀采样得到的31维光谱图像为例，包括以下步骤：通过人眼视锥响应构建加权矩阵Σ：人眼的三个视锥响应灵敏度L、M、S本身是31x1维的矢量，在用其构建加权矩阵Σ时，先将L、M、S矢量对角化，变成一个31x31维的对角矩阵(除了主对角线的元素不为零，其余元素全为零)。再将对角化后的L、M、S矩阵相加，得到加权矩阵Σ，则矩阵Σ也是一个31x31维的对角阵，只是该对角阵主角线上的元素是L、M、S对角阵主对角线上相应元素相加求和的结果。将加权矩阵Σ与光谱图像R相乘，得到加权后光谱图像G：如果某图像的像素数是m，则其光谱图像R是一个31xm维的矩阵，两者相乘得到加权后光谱图像G是一个31xm维的矩阵；再对加权光谱图像G进行主成分分析，得到第一组三个特征向量A1/A2/A3：在对光谱图像G进行主成分分析时，可得到31个特征向量，每个特征向量均是一个31x1维的矢量。而在本方法中只需在31个特征向量中取其前三个特征向量作为重构光谱图像的第一组特征向量A1/A2/A3，这组特征向量互相正交。以第一组的三个特征向量Ag＝(A1,A2,A3)作为六维线性模型的前三个基向量，求与基向量组A1/A2/A3对应的系数矩阵Wg：由光谱图像R通过最小二乘法求得，这样求得的系数矩阵Wg是一个3xm维的矩阵，m是图像像素数。用第一组的三个特征向量Ag＝本文档来自技高网...

【技术保护点】
一种用于颜色复制的光谱图像处理方法，其特征在于，该处理方法包括以下步骤： 通过人眼视锥响应构建加权矩阵Σ； 将该加权矩阵Σ与原始光谱图像R相乘，得到加权后光谱图像G，再对加权光谱图像G进行主成分分析，得到第一组三个特征向量A1/A2/A3； 取第一组的三个特征向量Ag＝(A1,A2,A3)作为六维线性模型的前三个基向量，求与这个三基向量对应的系数矩阵并用这三个基向量重构光谱，得到重构光谱图像通过原始光谱图像R减去三个基向量Ag＝(A1,A2,A3)重构的光谱图像得到残余光谱Rlost； 对残余光谱Rlost进行主成分分析，得到第二组三个特征向量B1/B2/B3； 取第二组的三个特征向量Blost＝(B1,B2,B3)，作为六维线性模型的后三个基向量，求与特征向量组B1/B2/B3对应的系数矩阵Wlost； 利用这六个基向量将原始高维光谱R进行降维处理得到六维多通道图像； 六维多通道图像的保存； 最后将六维多通道图像代替原始光谱图像用于颜色复制过程。 FDA0000424620710000011.jpg,FDA0000424620710000012.jpg,FDA0000424620710000013.jpg...

【技术特征摘要】
1.一种用于颜色复制的光谱图像处理方法，其特征在于，该处理方法包括以下步骤：通过人眼视锥响应构建加权矩阵Σ；将该加权矩阵Σ与原始光谱图像R相乘，得到加权后光谱图像G，再对加权后光谱图像G进行主成分分析，得到第一组三个特征向量A1/A2/A3；取第一组的三个特征向量Ag＝(A1,A2,A3)作为六维线性模型的前三个基向量，求与这三个基向量对应的系数矩阵Wg，并用这三个基向量重构光谱，得到重构光谱图像通过原始光谱图像R减去三个基向量Ag＝(A1,A2,A3)重构的光谱图像得到残余光谱Rlost；对残余光谱Rlost进行主成分分析，得到第二组三个特征向量B1/B2/B3；取第二组的三个特征向量Blost＝(B1,B2,B3)，作为六维线性模型的后三个基向量，求与特征向量组B1/B2/B3对应的系数矩阵Wlost；利用这六个基向量将原始高维光谱R进行降维处理得到六维多通道图像；六维多通道图像的保存；最后将六维多通道图像代替原始光谱图像用于颜色复制过程；所述加权矩阵Σ，满足通过人眼三种视锥细胞长波L(λ)，中波M(λ)，短波S(λ)对重构光谱图像反射率和原始光谱图像反射率R的视锥响应总的偏差最小情况下，由人眼长波L(λ)、中波M(λ)和短波S(λ)对角化后再进行矩阵相加构成加权矩阵Σ＝L+M+S。2.根据权利要求1所述用于颜色复制的光谱图像处理方法，其特征在于，所述重构光谱是通过所得，Wg是由加权后光谱图像G通过最小二乘法求得。3.根据权利要求1所述用于颜色复制的光谱图像处理方法，其特征在于，所述六维多通道图像由对应系数矩阵Wg和Wlost构成，系数矩阵Wg对应于六维多通道图像的前三个通道数据，系数矩阵Wlost对应于六维多通道图像的后三个通道数据。4.根据权利要求3所述用于颜色复制的光谱图像处理方法，其特征在于，所述六维多通道图像的系数矩阵Wlost由残余光谱Rlost通过最小...

【专利技术属性】
技术研发人员：何颂华，
申请(专利权)人：深圳职业技术学院，
类型：发明
国别省市：广东;44