张晓晓+方玉萍+李蒙+吕宪魁+王远方舟+董明+张洋+朱嘉玮+杨卫平+罗艳琳

摘要： 研究了基于选择最优滤光片多光谱成像系统优化问题，采用了优化的滤光片组合，组建了一个六通道的多光谱成像系统。利用提出的基向量和多项式组合的方式重构光谱数据，成功重建出博物馆、美术馆等展览馆常用的两种典型光源模拟器D65和A光源下油画的光谱反射率，并以光谱均方误差和色差综合误差评价成像质量。对两种光源下重建结果分析可得，最优滤光片的选取与光照环境密切相关。根据光照环境不同采用选取出的最优滤光片可得到最好光谱成像图，并且能在光谱和色度上都达到较高精度。

关键词： 多光谱成像； 最优滤光片； 基向量法； 多项式模型； 数码相机

中图分类号：O 433 文献标志码： A doi： 10.3969/j.issn.1005-5630.2016.05.004

文章编号： 1005-5630（2016）05-0393-009

引 言

多光谱成像技术早在20世纪70年代就应用于航空航天遥感领域，并且不断得到发展[1]。由于其克服了传统三原色数字成像设备固有的同色异谱问题，实现了物体表面光谱信息的真实记录。该技术逐步在博物馆、美术馆和计算机图形学得到应用[2]，为文化遗产真实的数字化典藏、艺术品的高保真复制提供基础[3]。当前在颜色科学领域，光谱成像设备主要分为由窄带滤光片组成的窄带光谱成像系统[4]和由宽带滤光片配合三通道或单色CCD数码相机组成的宽带成像系统[5]。宽带成像系统通过光谱重建算法计算得到物体表面的光谱反射率，并且只需要对具有高分辨率的数码相机进行简单的改造，即可获取具有高分辨率和高空间分辨率的光谱图像，较好地克服了窄带光谱成像系统的低空间分辨率的缺点，光谱图像获取周期短，所以得到了广泛应用[6]。

美国孟塞尔颜色科学实验室Berns等对Sinar54型号数码相机进行改造，配合从肖特（scott）滤光片中所选取出的两个宽带滤光片schottBG39和schottGG475组成6通道宽带光谱相机，重建出D50光源下油画的光谱反射率。所选的滤光片只是为了得到与改造后的数码相机有近似的光谱灵敏度，并未形成一个使用滤光片的选择方法[7]。Day在Berns等的指导下也进行了相应的滤光片选择实验，结合伪逆法与相机数学模型成功地选择出在光谱上达到最佳精度的滤光片柯达雷登55（Kodak Wratten 55）和unfiltered （NF）。在色度上达到较高精度的是柯达雷登60（Kodak Wratten 60）和unfiltered （NF），但只是对不同型号滤光片组合的性能进行了简单分析，且选取出的滤光片在短波段重建效果较差，因为实验中采用的光源主要分布在长波段[8]。由此可见，不同光源下对于滤光片的选择也是决定光谱成像系统能否精确获取数据的重要因素之一。

国内武汉大学参考美国孟塞尔颜色科学实验室的研究成果，建立了一套完整的宽带光谱成像系统，并在文物保护领域得到了有效应用[9-10]。而对于宽带光谱图像滤光片的选择研究工作，目前还尚未深入。本文针对展览馆或美术馆进行艺术品高保真复制时用到的模拟日光D65和偏红A光源这两种典型的光照环境，基于普通的商用数码相机与普通的吸收宽带滤光片搭建出多光谱成像平台，对不同光源下滤光片的选择进行研究。

3 实验样本的获取和数据处理

3.1 实验环境

搭建的宽带多光谱成像系统平台如图1所示。实验是在Image Quality Labs提供的反射测试光源下进行。左右两边照明光束的光轴方向与样品法线成45°，相机拍摄方向的法线和样品法线之间的夹角尽量重合，误差不超过2°，以确保待测样品在不同空间位置所受到的光辐射量是均匀的。

实验中采用国际照明委员会（CIE）针对颜色的测量和计算推荐的标准光源A和标准光源D65，这两种常见的光源的相对光谱功率分布用PR-715辐射分光光度计测得，如图2所示。

3.2 实验仪器及样本的选取

3.2.1实验仪器

为了检验光谱重建精度，需要预先测量训练样本和测试样本的光谱反射率。本实验使用X-RiteSP64分光光度计来测量实验样本的光谱反射率，测量范围为400～700 nm，采样间隔为10 nm，排除镜面反射。采用佳能EOS400D配合不同颜色的宽带滤光片来采集实验样本的颜色信息。通过反复实验找到最佳拍摄方法。采用光学变焦镜头，焦距为38 mm，光圈值为1/2.8，ISO设置为400，曝光时间为1/160 s，聚焦设置为手动模式，且拍摄过程中相机保持静止不动。拍摄完后，通过在MATLAB平台上编程提取实验样本的RGB值。

3.2.2滤光片及实验样本

滤光片的选择决定了多光谱成像系统获取数据的精度，实验采用高透光率的卓美P系列方形滤光片，每个滤光片厚度为1.76 mm，尺寸为83 mm×95 mm。滤光片共11种颜色，分别为全粉、全蓝、全紫、全橙、全红、全黄、全茶、全绿、灰ND2、灰ND4、灰ND8。11种颜色滤光片的透过率如图3所示。

选用爱色丽公司的标准色卡Digital Color Checker SG（简称SG）中间的96个色块（外围的44色是三种中性色的重复）作为训练样本集，如图4所示其彩色色块大都代表自然界的真实颜色。待测样本为D65光源和A光源下的油画，针对油画具有代表性的不同颜色参考区域进行研究。

4 颜色再现的分析与评价

最优的滤光片组合应该具有色差ΔE00最小，光谱均方根误差（RSME）最小以及颜色主观因素影响的评分值最大。评价流程如图5所示。

4.1 D65和A光源下的颜色再现研究

在标准光源D65和标准光源A下拍摄的测试样本如图6所示，其中选取了8种（在图中用1、2、3、4、5、6、7、8标出）具有代表性的不同颜色的部分区域进行研究[20]。

对55组滤光片组合再现的油画颜色，利用主观评价的方法，通过视觉判断将误差太大的组合筛选出去，在光源D65和A光源下分别选取出6组重构图像效果较好的数字油画，如图7、图8所示。

将两种光源下选取的6种组合进行进一步的分析研究，分别采用6组滤光片的组合对其重建，结果如表1、表2所示。ΔE00是国际照明委员会于2001年正式推荐的CIE2000色差公式的色差符号。

基于两种光源下分别选出的6种组合，需要选出一组最优滤光片。其分析结果如表3所示。

4.2 D65和A光源下的颜色再现分析

在D65光源下，从表1重构的数据可看出，每个色块在蓝绿滤光片组合下的误差相对较小，由于所研究的油画色调偏暖，对其中冷色调的5号和8号色块进行分析，在蓝绿组合的滤光片下其误差也比较小，这进一步说明，蓝绿滤光片的组合不仅可以较好地重建与其成对立色的颜色，同时也可以照顾到一些与其成补色的颜色。所选取的最优滤光片同样适用于冷色调的油画。同时对表3所选取的色块进行综合误差分析可知，全蓝和全绿滤光片组合的SymbolDA@ E00和光谱均方根误差（RMSE）的综合误差优于其他组合。蓝和绿滤光片的光谱透射率如图9所示。

对表2及表3分析可知，在A光源下8个色块在全红和全蓝这两个滤光片组合下，SymbolDA@ E00和RMSE都很小，包括冷色调的5号和8号在红蓝组合的滤光片下误差也比较小。所以所选取的最优滤光片同样适用于冷色调的油画。其综合误差分析结果也是红蓝滤光片组合优于其它组合。滤光片红和蓝的光谱透射率如图10所示。

图9 滤光片蓝和绿的光谱透射率

Fig.9 The spectral transmittance of blue and green filter

深入分析可知，测试样本再现效果的好坏与光源密切相关。A光源偏暖，光谱功率分布主要在长波波段，滤光片1的长波波段透过率较高，并由D65光源下滤光片的选取分析可知，蓝色滤光片也可以较好地重建与其成补色的颜色，所以在红蓝滤光片下重现效果最好。D65光源色温相比A光源偏冷，所以在蓝绿滤光片组合下重现效果最好。

4.3 D65和A光源下的颜色再现主观评价

主观评价要求所有的评测者均无色盲和色弱现象且视力正常，平均年龄在20岁左右，男性25人，女性25人，其中有颜色技术背景的15人，没有颜色技术背景的35人。将显示器的背景设置为灰色，色温设置为6 500 K，色空间设置为sRGB色空间，源场景和显示器均使用可调光照强度的标准D65光源和A光源照明，调节光照强度使人眼感受到的源场景亮度与显示图像亮度持平，然后所有的评测者对图7（a）～（f）6幅测试图颜色再现的效果进行主观评价。同样，将显示器的色温设置为2 856 K（A光源），重复上述步骤再对图8（a）～（f）6幅测试图颜色再现的效果进行主观评价。0～0.25表示颜色再现质量非常差，0.25～0.5表示颜色再现质量一般，0.5～0.75表示颜色再现质量好，0.75～1.00表示颜色再现质量非常好。表4为主观评价的实验结果。

由表4可知，所有场景再现的效果都在可接受的范围之内，在D65和A光源下再现效果最好的分别是蓝绿滤光片和红蓝滤光片下的再现图像，主观评价与客观评价结果一致。实验表明，所选取的滤光片组合为55种组合中最优。

5 结 论

本文是基于普通数码相机结合优化的滤光片搭建了可用于实际颜色复制应用的多光谱图像重建系统。通过实验分析知，滤光片的选取与光源密切相关，在D65光源下选取出的最优滤光片组合为蓝、绿，重构光谱的色差最小可以达到1.53，RMSE最小可达到0.045 1。而Berns等在D50光源下所选取出的两个宽带滤光片schottBG39和schottGG475的光谱特性曲线与本文所研究的蓝色和绿色滤光片光谱特性曲线基本一致。标准照明体D50和D65都是日光模拟器，所以在模拟日光环境下，选用蓝、绿滤光片搭建的多光谱成像系统基本会达到最佳精度。在A光源下选取出的最优滤光片组合为红、蓝，重构光谱的色差最小可以达到0.66，RMSE最小可达到0.046 0，都具有较高的重构精度。并对再现颜色的质量进行了主观评价，与客观评价的结果相一致。初步实验结果表明，根据光照环境不同采用选出的最优滤光片可得到最好的光谱成像图，能满足人们在日常生活中对于油画等艺术品高保真图像再现的基本需求，可以在实际中得到广泛应用。

参考文献：

[1] CURRAN P J.Principle of remote sencing[M].Harlow，UK：Addison Wesley Longman，1985.

[2] ZHAO Y，LI P S，WANG C S，et al.A novel fiber-optic sensor used for small internal curved surface measurement[J].Sensors and Actuators A，2000，86（3）：211-215.

[3] CREAGH D，BRADLEY D.Physical techniques in the study of art，archaeology and cultural heritage[M].Amsterdam：Elsevier，2007.

[4] HARDEBERG J Y，SCHMITT F，BRETTEL H.Multispectral color image capture using a liquid crystal tunable filter[J].Optical Engineering，2002，41（10）：2532-2548.

[5] SUGIURA H，KUNO T，WATANABE N，et al.Development of a multispectral camera system[C]//Proceedings of SPIE3965，Sensors and Camera Systems for Scientific，Industrial，and Digital Photography Applications.San Jose，CA：SPIE，2000，3965：331-339.

[6] BERNS R S，BYRNS S，CASADIO F，et al.Rejuvenating the color palette of Georges Seurats A Sunday on La Grande Jatte—1884：a simulation[J].Color Research & Application，2006，31（4）：278-293.

[7] BERNS R S，TAPLIN L A.Practical spectral imaging using a color-filter array digital camera[R].Rochester，New York，United States：Rochester Institute of Technology，Munsell Color Science Laboratory，2006.

[8] DAY D C.Filter selection for spectral estimation using a trichromatic camera[D].New York：Rochester Institute of Technology，2001.

[9] 刘振.面向复制的多光谱图像获取技术研究[D].武汉：武汉大学，2012.

[10] 刘振，万晓霞，黄新国，等.基于宽带多通道的光谱反射率重建方法研究[J].光谱学与光谱分析，2013，33（4）：1076-1081.

[11] MALONEY L T.Evaluation of linear models of surface spectral reflectance with small numbers of parameters[J].Journal of the Optical Society of America A，1986，3（10）：1673-1683.

[12] FAIRMAN H S，BRILL M H.The principal components of reflectances[J].Color Research & Application，2004，29（2）：104-110.

[13] TZENG D Y，BERNS R S.A review of principal component analysis and its applications to color technology[J].Color Research & Application，2005，30（2）：84-98.

[14] AGAHIAN F，AMIRSHAHI S A，AMIRSHAHI S H.Reconstruction of reflectance spectra using weighted principal component analysis[J].Color Research & Application，2008，33（5）：360-371.

[15] HONG G W，LUO M R，RHODES P A.A study of digital camera colorimetric characterization based on polynomial modeling[J].Color Research &Application，2001，26（1）：76-84.

[16] 王莹，王忠民，王义峰，等.面向色彩再现的多光谱图像非线性降维方法[J].光学 精密工程，2011，19（5）：1171-1178.

[17] 宓超波，黄小乔，石俊生.普通办公照明环境下数码照相机的颜色特性化[J].光学学报，2013，33（9）：0933001.

[18] PARK J I，LEE M H，GROSSBERG M D，et al.Multispectral imaging using multiplexed illumination[C]//Proceedings ofthe IEEE 11th International Conference on Computer Vision.Rio de Janeiro，Brazil：IEEE，2007：1-8.

[19] IMAI F H，BERNS R S，TZENG D Y.Acomparative analysis of spectral reflectance estimated in various spaces using a trichromatic camera system[J].Journal of Imaging Science and Technology，2000，44（4）：280-287.

[20] 秦峰，杨卫平，杨葭，等.基于亮度信息匹配的国画艺术品图像重建研究[J].光学学报，2014，34（10）：1033001.