CN106153193B - 一种利用多光谱相机双光源响应值获取光谱反射比的方法 - Google Patents
一种利用多光谱相机双光源响应值获取光谱反射比的方法 Download PDFInfo
- Publication number
- CN106153193B CN106153193B CN201610590506.1A CN201610590506A CN106153193B CN 106153193 B CN106153193 B CN 106153193B CN 201610590506 A CN201610590506 A CN 201610590506A CN 106153193 B CN106153193 B CN 106153193B
- Authority
- CN
- China
- Prior art keywords
- light source
- under
- response
- mrow
- spectral reflectance
- 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.)
- Expired - Fee Related
Links
- 230000003595 spectral effect Effects 0.000 title claims abstract description 43
- 238000000034 method Methods 0.000 title claims abstract description 19
- 238000002310 reflectometry Methods 0.000 claims abstract description 8
- 239000011159 matrix material Substances 0.000 claims description 24
- 230000007704 transition Effects 0.000 claims description 8
- 238000005286 illumination Methods 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 6
- 239000000284 extract Substances 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 3
- 239000003086 colorant Substances 0.000 claims description 2
- 241000532370 Atla Species 0.000 claims 1
- 230000004298 light response Effects 0.000 claims 1
- 238000001228 spectrum Methods 0.000 abstract description 10
- 238000003384 imaging method Methods 0.000 description 6
- 230000004075 alteration Effects 0.000 description 2
- 239000004973 liquid crystal related substance Substances 0.000 description 2
- 238000002834 transmittance Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000001579 optical reflectometry Methods 0.000 description 1
- 238000000985 reflectance spectrum Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/2823—Imaging spectrometer
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/2823—Imaging spectrometer
- G01J2003/2826—Multispectral imaging, e.g. filter imaging
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
本发明公开了一种利用多光谱相机双光源响应值获取光谱反射比的方法,获得物体光谱反射比可以实现高精度的颜色测量,且能够忠实复现物体在任意光源下的颜色。首先利用训练样本建立多光谱相机每个通道在两种光源下响应值之间的函数关系,然后在一种光源下拍摄目标样本,并利用之前建立的两种光源下的通道响应值函数关系以及该光源下的所有通道响应值,计算出另一种光源下的全部通道响应值,从而利用两种光源下的多通道响应值获取物体表面的光谱反射比,显著增加了可利用的目标样本的响应值信息,提高了光谱反射比的获取精度,同时不需要在另一种光源下实际拍摄目标样本,简化了操作,提高了效率。
Description
技术领域
本发明涉及成像式获取物体表面光谱及色度信息的方法,尤其是利用具有多个可见光波段通道的多光谱相机获取物体的反射光谱信息以进一步获取物体色度信息的方法。
背景技术
由于同色异谱现象的存在,物体的色度参数难以可靠表征物体表面的本源信息,而物体表面的光谱反射比是设备无关和捕获光源无关的物理量,表征了物体表面的本源信息,因此准确获取物体光谱反射比可以忠实复现物体在任意光源下的色貌。
成像式获取物体光谱反射比可以克服分光光度计逐点式测量和接触式测量的缺陷,能够快速获取成像范围内所有物体的光谱反射比。普通商用相机可以用于成像式获取物体光谱反射比,但由于只有R、G、B三个通道而无法高精度获取光谱反射比。而多光谱相机用于成像式获取物体光谱反射比,由于增加了成像的通道数,所以可以大幅度提高光谱反射比的获取精度。
当采用更多的通道响应值获取物体光谱反射比时,获取精度会进一步提高。利用目标样本在两种光源下的通道响应值可以显著增加目标样本通道响应值的数量,但在两种光源下实际拍摄会增加操作复杂度和多光谱图像采集时间,如在另一种光源下需要重新调整每一个通道的曝光时间,对每个通道对焦,多通道图像之间的配准等。
发明内容
为了利用更多的通道响应值获取物体表面的光谱反射比,并减少多光谱相机的操作复杂度和多通道图像采集时间,本发明提出了一种利用多光谱相机双光源响应值获取光谱反射比的方法。
本发明的目的是通过以下技术方案来实现的:一种利用多光谱相机双光源响应值获取光谱反射比的方法,包括以下步骤:
(1)利用训练样本建立两种光源下多光谱相机通道响应值之间的函数关系。具体包括以下子步骤:
(1.1)在第一种光源a下和第二种光源b下分别用多光谱相机拍摄训练样本。采用包含各种颜色块的标准色卡作为训练样本,将标准色卡以45/0照明与测量几何条件放置于光源a下,用多光谱相机的所有通道拍摄标准色卡,提取光源a下标准色卡的全部通道响应值Pa为n行q列的矩阵,其中n为标准色卡上的色块数,q为多光谱相机通道的数量,每一行表示一个色块对应的q个通道的响应值,每一列表示一个通道输出的n个色块的响应值。再将标准色卡以同样的照明与测量几何条件放置于光源b下进行拍摄,提取标准色卡在光源b下全部通道的响应值Pb也为n行q列的矩阵。
(1.2)建立每一个通道在a和b两种光源下响应值之间的函数关系。假设pa1和pb1分别为矩阵Pa的第1列和Pb的第1列,分别对应标准色卡在a和b光源下第1个通道的n个色块的响应值。通过二次多项式拟合光源b下第1个通道的响应值pb1与光源a下第1个通道响应值pa1之间的关系,即
pb1=k0+k1*pa1+k2*pa1 2
利用pa1和pb1,通过最小二乘法求解上述多项式的系数k=[k0,k1,k2],即
k=(CTC)-1CTpb1
C=[o,pa1,pa1οpa1]
其中,为含有n个元素都为1的列向量,ο为矩阵或向量的Hadamard积。
用上述同样的方法,类似地可以求出其余q-1个通道在b光源下的响应值与在a光源下的响应值之间的二次多项式函数。
(2)计算光谱反射比与多通道响应值之间的转换矩阵。用分光光度计测量标准色卡的所有色块的光谱反射比R为n行m列的矩阵,其中n为标准色卡上的色块数,m为光谱反射比的维数。根据两种光源下通道响应值之间的二次多项式函数,利用标准色卡在光源a下的响应值Pa计算出其在光源b下的通道响应值则标准色卡在两种光源下的所有通道的响应值为然后,由最小二乘法计算从多通道响应值到光谱反射比的转换矩阵M,即
(3)在光源a下拍摄目标样本,并计算目标样本的光谱反射比。在光源a下,用多光谱相机拍摄目标样本,得到目标样本在光源a下的通道响应值Pat,然后根据步骤(1)中光源b下通道响应值与光源a下通道响应值之间的二次多项式函数,计算出目标样本在光源b下的通道响应值则目标样本在两种光源下的通道响应值为最后,根据步骤(2)中多通道响应值到光谱反射比的转换矩阵M,计算出目标样本的光谱反射比
本发明的有益效果是:本发明预先建立多光谱相机在两种光源下通道响应值之间的函数关系,然后仅在一种光源下对目标样本成像,并根据两种光源下通道响应值之间的函数关系,利用目标样本在一种光源下的通道响应值计算出在第二种光源下的通道响应值,从而显著增加了目标样本的通道响应值信息,提高了光谱反射比的获取精度,同时不需要在第二种光源下进行实际拍摄,大大降低了操作的复杂度和多通道图像的采集时间。
附图说明
图1是利用多光谱相机双光源响应值获取光谱反射比的流程图。
图2是实例中多光谱相机结构示意图。
图3是实例中采用双光源响应值前后获得的光谱反射比,(a)、(b)分别为两个色块的示意图。
具体实施方式
以一台滤色片轮式的多光谱相机为例,阐述利用多光谱相机双光源响应值获取物体表面光谱反射比的方法,多光谱相机的结构示意图如图2,由镜头、滤色片轮、滤色片轮控制单元和单色CCD传感器构成。该多光谱相机利用8个干涉滤色片形成8个通道。8个滤色片的峰值透过率波长为从420nm到700nm,FWHM(峰值半高宽,Full Width of Half Maximum)为20nm,因此8个通道对应的峰值透过率波长为420nm,460nm,500nm,540nm,580nm,620nm,660nm,700nm。需要说明的是,本发明不仅限于滤色片轮式的多光谱相机,对于其它形式的多光谱相机,如基于LCTF(液晶可调滤色片,Liquid Crystal Tuneable Filter)的多光谱相机也适用;同时,本发明也不仅限于通道数量为8的多光谱相机,少于或多于该滤色片数量的多光谱相机也适用。
如图1所示,本发明利用多光谱相机双光源响应值获取光谱反射比的方法,具体包括以下步骤:
(1)利用训练样本建立两种光源下多光谱相机通道响应值之间的函数关系。具体包括以下子步骤:
(1.1)采用GretagMacbeth SpectraLightⅢ灯箱中的D65光源和A光源分别作为两个照明光源,选择GretagMacbeth ColorChecker DC色卡(DC色卡)作为训练样本,将DC色卡以45/0照明与测量几何条件放置于D65照明光源下,用多光谱相机拍摄DC色卡,提取在D65光源下DC色卡的所有色块的通道响应值PD65为240行8列的矩阵,其中240为DC色卡上的色块数,8为多光谱相机的通道数量,每一行代表一个色块对应于8个通道的响应值,每一列则代表一个通道下240个色块的响应值。再将DC色卡以同样的照明与测量几何条件放置于A光源下进行拍摄,提取A光源下DC色卡的所有色块的通道响应值PA也为240行8列的矩阵。
(1.2)建立每一个通道在D65和A两种光源下响应值之间的函数关系。假设pD651和pA1分别为矩阵PD65的第1列和PA的第1列,分别对应于DC色卡在D65和A光源下第1个通道(420nm)对240个色块的响应值。通过二次多项式拟合A光源下第1个通道响应值pA1与D65光源下第1个通道响应值pD651之间的关系,即
pA1=k0+k1*pD651+k2*pD651 2
利用pD651和pA1,通过最小二乘法求解上述多项式的系数k=[k0,k1,k2],即
k=(CTC)-1CTpA1
C=[o,pD651,pD651οpD651]
其中,为含有240个元素都为1的列向量,ο为矩阵或向量的Hadamard积。
用上述同样的方法,类似地可以求出其余7个通道(460nm,500nm,540nm,580nm,620nm,660nm,700nm)在A光源下响应值与在D65光源下响应值之间的二次多项式函数。
(2)计算光谱反射比与多通道响应值之间的转换矩阵。用分光光度计测量DC色卡的所有色块的光谱反射比R为240行31列的矩阵,其中240为标准色卡上的色块数,31为光谱反射比的维数,即400nm到700nm,10nm为间隔。根据两种光源下的二次多项式函数,利用DC色卡在光源D65下的响应值计算出其在光源A下的通道响应值则DC色卡在两种光源下的通道的响应值为然后,由最小二乘法计算从多通道响应值到光谱反射比的转换矩阵
(3)在光源D65下拍摄目标样本,并计算目标样本的光谱反射比。采用X-Rite公司的ColorChecker Color Rendition Chart色卡(CRC色卡)作为目标样本。在光源D65下,用多光谱相机拍摄CRC色卡,提取CRC色卡在D65光源下的通道响应值其中24为CRC色卡上的色块数,8为多光谱相机通道数。根据步骤(1)中A光源下通道响应值与D65光源下通道响应值之间的二次多项式函数,计算出CRC色卡在光源A下的通道响应值则CRC色卡在两种光源下的通道响应值为最后,根据步骤(2)中多通道响应值到光谱反射比的转换矩阵M,计算出目标样本的光谱反射比 即
实验结果:未采用CRC色卡双光源通道响应值即仅用D65光源下的通道响应值时,所获得的CRC色卡24个色块光谱反射比的平均光谱差(RMSE,均方根误差)为0.0135,平均色差(ΔE00,CIEDE2000色差)为1.23。采用双光源通道响应值即包含D65光源下的实测响应值以及A光源下计算得到的响应值后,所获得CRC色卡24个色块光谱反射比的平均光谱差(RMSE)为0.0108,平均色差(ΔE00)为0.80。由此可见,采用双光源下通道响应值所获取的光谱反射比,其光谱差和色差均明显减小。附图3展示了CRC色卡中第一行左起第1个和第3个色块采用双光源响应值前后所获取的光谱反射比与其实际反射光谱曲线的比较,可见其采用双光源响应值所获取的光谱反射比更接近于实际值。
Claims (1)
1.一种利用多光谱相机双光源响应值获取光谱反射比的方法,其特征在于,包括以下步骤:
(1)利用训练样本建立两种光源下多光谱相机通道响应值之间的函数关系;具体包括以下子步骤:
(1.1)在第一种光源a下和第二种光源b下分别用多光谱相机拍摄训练样本;采用包含各种颜色块的标准色卡作为训练样本,将标准色卡以45/0照明与测量几何条件放置于光源a下,用多光谱相机的所有通道拍摄标准色卡,提取光源a下标准色卡的全部通道响应值Pa为n行q列的矩阵,其中n为标准色卡上的色块数,q为多光谱相机通道的数量,每一行表示一个色块对应的q个通道的响应值,每一列表示一个通道输出的n个色块的响应值;再将标准色卡以同样的照明与测量几何条件放置于光源b下进行拍摄,提取标准色卡在光源b下全部通道的响应值Pb也为n行q列的矩阵;
(1.2)建立每一个通道在a和b两种光源下响应值之间的函数关系;假设pa1和pb1分别为矩阵Pa的第1列和Pb的第1列,分别表示标准色卡在a和b光源下第1个通道的n个色块的响应值;通过二次多项式拟合光源b下第1个通道的响应值pb1与光源a下第1个通道响应值pa1之间的关系,即
pb1=k0+k1*pa1+k2*pa1 2
利用pa1和pb1,通过最小二乘法求解上述多项式的系数k=[k0,k1,k2],即
k=(CTC)-1CTpb1
其中,为含有n个元素都为1的列向量,为矩阵或向量的Hadamard积;
用上述同样的方法,类似地可以求出其余q-1个通道在b光源下的响应值与在a光源下的响应值之间的二次多项式函数;
(2)计算光谱反射比与多通道响应值之间的转换矩阵;用分光光度计测量标准色卡中所有色块的光谱反射比R为n行m列的矩阵,其中n为标准色卡上的色块数,m为光谱反射比的维数;根据两种光源下通道响应值之间的二次多项式函数,利用标准色卡在光源a下的响应值Pa计算出其在光源b下的通道响应值则标准色卡在两种光源下的所有通道的响应值为然后,由最小二乘法计算从多通道响应值到光谱反射比的转换矩阵M,即
<mrow>
<mi>M</mi>
<mo>=</mo>
<msup>
<mi>R</mi>
<mi>T</mi>
</msup>
<msub>
<mover>
<mi>P</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mover>
<mi>P</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mi>T</mi>
</msup>
<msub>
<mover>
<mi>P</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
(3)在光源a下拍摄目标样本,并计算目标样本的光谱反射比;在光源a下,用多光谱相机拍摄目标样本,得到目标样本在光源a下的通道响应值Pat,然后根据步骤(1)中光源b下通道响应值与光源a下通道响应值之间的二次多项式函数,计算出目标样本在光源b下的通道响应值则目标样本在两种光源下的通道响应值为最后,根据步骤(2)中多通道响应值到光谱反射比的转换矩阵M,计算出目标样本的光谱反射比
<mrow>
<msub>
<mover>
<mi>R</mi>
<mo>^</mo>
</mover>
<mi>t</mi>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>M</mi>
<msubsup>
<mover>
<mi>P</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>a</mi>
<mi>b</mi>
<mi>t</mi>
</mrow>
<mi>T</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>.</mo>
</mrow>
1
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610590506.1A CN106153193B (zh) | 2016-07-22 | 2016-07-22 | 一种利用多光谱相机双光源响应值获取光谱反射比的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610590506.1A CN106153193B (zh) | 2016-07-22 | 2016-07-22 | 一种利用多光谱相机双光源响应值获取光谱反射比的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106153193A CN106153193A (zh) | 2016-11-23 |
CN106153193B true CN106153193B (zh) | 2017-12-15 |
Family
ID=58060096
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610590506.1A Expired - Fee Related CN106153193B (zh) | 2016-07-22 | 2016-07-22 | 一种利用多光谱相机双光源响应值获取光谱反射比的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106153193B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106895916B (zh) * | 2017-01-09 | 2018-10-30 | 浙江大学 | 一种单次曝光拍摄获取多光谱图像的方法 |
CN106896068B (zh) * | 2017-04-12 | 2023-07-18 | 齐鲁工业大学 | 解决网购颜色失真的商品颜色获取与再现方法及其系统 |
CN110726536B (zh) * | 2019-09-25 | 2021-08-06 | 宁波永新光学股份有限公司 | 一种彩色数码反射显微镜颜色校正方法 |
CN110660112B (zh) * | 2019-09-29 | 2021-09-24 | 浙江大学 | 一种基于专用色卡和多光谱成像的绘画光谱重构方法 |
CN111750995B (zh) * | 2020-06-04 | 2022-10-14 | 武汉纺织大学 | 一种面向开放测量环境应用的光谱测量方法 |
-
2016
- 2016-07-22 CN CN201610590506.1A patent/CN106153193B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN106153193A (zh) | 2016-11-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106153193B (zh) | 一种利用多光谱相机双光源响应值获取光谱反射比的方法 | |
CN106153192B (zh) | 一种利用多光谱相机虚拟响应值获取光谱反射比的方法 | |
US20180270462A1 (en) | Imaging processing device and imaging processing method | |
EP3054273B1 (en) | Colorimetry system for display testing | |
CN103201602B (zh) | 具有至少两个独立数字照相机的数字多光谱照相机系统 | |
WO2017193738A1 (zh) | 图像传感器、成像方法和成像装置 | |
CN104390703A (zh) | 用于确定光谱仪的校准参数的方法 | |
CN110017904B (zh) | 一种基于ccd相机的多光谱辐射测温方法 | |
CN1971927A (zh) | 物理信息获取方法、物理信息获取装置和半导体器件 | |
TW200539690A (en) | Image pickup device | |
CN101874401A (zh) | 用于测量物体活力的单芯片图像传感器 | |
KR100323503B1 (ko) | 촬상장치의 분광응답도 특성 측정방법 및 촬상데이타의 구성방법 | |
Shrestha et al. | One-shot multispectral color imaging with a stereo camera | |
US9188841B2 (en) | Imaging device | |
CN101294867B (zh) | Ccd成像设备相对光谱响应特性的标定方法及系统 | |
CN106841055B (zh) | 一种重构艺术绘画光谱图像的训练样本选择方法 | |
US9432642B2 (en) | Imaging system and method having extended depth of field | |
CN109738068B (zh) | 一种校正多光谱相机响应值非线性的方法 | |
EP2672719A2 (en) | Color calibration of a image capture device in a way that is adaptive to the scene to be captured | |
Tominaga | CIC@ 20: multispectral imaging | |
CN105628190A (zh) | 一种基于滤光单元的光辐射测量方法及其装置 | |
JP4344257B2 (ja) | 撮影装置及び色分解光学系 | |
CN105744267B (zh) | 基于齐次多项式数码相机参数可变的获取三刺激值方法 | |
US7046288B1 (en) | Image recording apparatus employing a single CCD chip to record two digital optical images | |
KR102459343B1 (ko) | 치과용 컬러 카메라를 교정하기 위한 컬러 기준 부재 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171215 |