基于多光谱LED照明的物体光谱反射率重建方法与流程

本发明涉及物体光谱反射率重建领域，特别涉及一种基于多光谱led照明的物体光谱反射率重建方法。

背景技术：

颜色信息的传递和真实再现与人们的日常生活密切相关。数字图像在人们日常生活中的应用越来越广泛，人们在享受数字图像设备带来巨大便利的同时也不同程度上受到彩色图像与原物不一致的困扰。基于光谱的颜色复制技术能有效避免同色异谱现象带来的各种困扰一些对颜色重现要求较高的应用领域需要获取目标表面上各点的光谱反射率。分光光度计能直接测量物体的光谱反射率，但其对测量面积、平整度都有限制，且效率较低，光谱扫描仪效率较高但不能进行非接触式测量且测量面积有限，显然都不能满足国画、油画等艺术品基于光谱的颜色复制要求。多光谱相机能高效获取到物体表面的光谱反射率，但目前的多光谱相机均在不同程度上存在结构复杂、价格昂贵，效率低等问题限制了其在实际中应用。数码相机作为一种主要的图像获取工具，在人们的日常生活中得到了广泛应用，但普通的商用数码相机只有三通道，只能采集特定照明条件下的颜色信息，也不能满足高精度颜色复制的需要。因此，如何高效地获取物体表面的光谱反射率成为颜色科学研究领域的热点。

技术实现要素：

本发明要解决的技术问题是提供一种在给定照明条件和观测条件的情况下，由ccd(电荷耦合器件)数字相机的响应值重建物体表面光谱反射率的方法。

为解决上述技术问题，本发明是通过以下技术方案实现的：一种基于多光谱led照明的物体光谱反射率重建方法，包括以下步骤：

步骤一：根据cdd成像系统的光谱响应曲线参数和光线波长对应光散射能力选择光源，优先选择偏红色光源。

步骤二：从现有颜色数据库中提取p个颜色基函数，通过该p个颜色基函数的线性叠加重建被测物体反射率，设bl(λ)表示物体反射率的正交基函数，σl表示权重，λ为光波长，重建的被测物体反射率可如下得到：

步骤三：选取n种led发光芯片使它们的光谱线性叠加后覆盖波长为300纳米至830纳米的范围，按照光谱密度函数中心波长从小到大依次将n种led发光芯片排列布置在铝基板上形成led阵列，再分别对每种led发光芯片的光谱密度函数取出n个离散波长点，n种led发光芯片的光谱密度函数记为si(k),i＝1,2,…,n；k＝1,2,…,n；

步骤四：记n种led发光芯片的点亮权值为wi,i＝1,2,…n，0≤wi≤1；wi组成权值矩阵w＝[w1,w2,…,wn]^t，则第k个混合光源的光谱可表示为：

式中s(k)＝[s1(k),s2(k),…,sn(k)]，s(k)为led阵列在第k个离散波长点光谱矩阵，将所有led光谱密度函数si(k)在n个离散波长点的功率构造成n行n列的离散光谱密度函数矩阵s，

根据脊回归模型构造无约束最优化问题可得：

记性能函数为：

令i为单位矩阵，获得唯一解：

w＝(s^tds+i)^-1ys^ty

计算获得所需要目标光谱的照明光源；

将权值矩阵w发送至驱动电路板控制led芯片的亮度，亮度的范围为0-100％，从而产生测量系统的照明光源，根据系统需要选取p种led照明光源配方；

步骤五：固定相机的曝光量，分别在l种照明光源条件下，通过已知光谱反射率的标准色卡的过光谱图像优化获得相机响应缩放系数αl，其中l＝1,2,…,p；

步骤六：设相机传感器灵敏度为s(λ)，则l种照明光源的相对光谱功率分布为：

式中bl(λ)表示物体反射率的正交基函数；

步骤七：tl时刻相机的输出为ol(tl)，则σl＝αlol(tl)，αl为第l种照明光源对应的相机响应缩放系数，计算p个颜色基函数权值σl，并代入式(1)重建光谱反射率。

本发明的有益效果是在给定照明条件和观测条件的情况下，由ccd(电荷耦合器件)数字相机的响应值重建物体表面光谱反射率，从而高效地获取物体表面的光谱反射率。

附图说明

图1为本实施例的流程图；

图2为ccd感光单元对不同波长的光线的响应能力的示意图。

具体实施方式

为了使本发明的技术方案更加清楚明白，以下结合附图及实施例，对本发明作进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。

参考图1至图2，本发明提供一种基于多光谱led照明的物体光谱反射率重建方法，根据ccd的感光单元对于不同波长的光线的响应能力是不同的和瑞利定律“线变小于光的波长的微粒，对入射光散射所遵循的规律是：散射光和入射光波长相同，散射光的强度和散射方向有关，并和波长的四次方成反比”，所以需针对不同的cdd成像系统的光谱响应曲线的具体参数以及光线波长对应光散射能力进行光源的折衷选择，优先选择偏红色光源。可以按照不同照射角度、不同颜色组合，从而更能突显出物体信息；节省安装空间、可解决对角照射阴影问题、光线均匀扩散可解决一致性的问题。

本发明包括以下步骤：

s1：根据cdd成像系统的光谱响应曲线参数和光线波长对应光散射能力选择光源，优先选择偏红色光源。

s2：从现有颜色数据库中提取p个颜色基函数，通过该p个颜色基函数的线性叠加重建被测物体反射率，设bl(λ)表示物体反射率的正交基函数，σl表示权重，λ为光波长，重建的被测物体反射率可如下得到：

s3：选取n种led发光芯片使它们的光谱线性叠加后覆盖波长为300纳米至830纳米的范围，按照光谱密度函数中心波长从小到大依次将n种led发光芯片排列布置在铝基板上形成led阵列，再分别对每种led发光芯片的光谱密度函数取出n个离散波长点，n种led发光芯片的光谱密度函数记为si(k),i＝1,2,…,n；k＝1,2,…,n；

s4：记n种led发光芯片的点亮权值为wi,i＝1,2,…n，0≤wi≤1；wi组成权值矩阵w＝[w1,w2,…,wn]^t，则第k个混合光源的光谱可表示为：

式中s(k)＝[s1(k),s2(k),…,sn(k)]，s(k)为led阵列在第k个离散波长点光谱矩阵，将所有led光谱密度函数si(k)在n个离散波长点的功率构造成n行n列的离散光谱密度函数矩阵s，

根据脊回归模型构造无约束最优化问题可得：

记性能函数为：

令i为单位矩阵，获得唯一解：

w＝(s^tds+i)^-1ys^ty

计算获得所需要目标光谱的照明光源；

将权值矩阵w发送至驱动电路板控制led芯片的亮度，亮度的范围为0-100％，从而产生测量系统的照明光源，根据系统需要选取p种led照明光源配方；

s5：固定相机的曝光量，分别在l种照明光源条件下，通过已知光谱反射率的标准色卡的过光谱图像优化获得相机响应缩放系数αl，其中l＝1,2,…,p；

s6：设相机传感器灵敏度为s(λ)，则l种照明光源的相对光谱功率分布为：

式中bl(λ)表示物体反射率的正交基函数；

s7：tl时刻的相机输出为ol(tl)，则σl＝αlol(tl)，αl为第l种照明光源对应的相机响应缩放系数，计算p个颜色基函数权值σl，并代入式(1)重建光谱反射率。

以上实施例仅仅是对本发明的解释，其并不是对本发明的限制，对于本领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围，因此本发明专利的保护范围应以权利要求为准。