发明内容
本发明提供了一种用于颜色三刺激值计算的光谱预处理方法,其目的是为了解决传统的光谱预处理方法不适合处理无明显特征峰值,曲线变化趋势比较平缓的,用于颜色三刺激值计算的光谱响应数据的问题。
为了达到上述目的,本发明的实施例提供了一种用于颜色三刺激值计算的颜色测量平台,包括:
反射测量支架;
待测样品,所述待测样品放置在所述反射测量支架的顶部;
第一准直镜,所述第一准直镜设置在所述待测样品的正上方;
第二准直镜,所述第二准直镜倾斜设置在所述待测样品的侧上方,所述第二准直镜与所述第一准直镜的之间的夹角呈45°;
光纤光源,所述光纤光源通过入射光纤与所述第二准直镜连接;
光纤光谱仪,所述光纤光谱仪的第一端通过出射光纤与所述第一准直镜连接;
计算机,所述计算机与所述光纤光谱仪的第二端电连接。
本发明的实施例还提供了一种用于颜色三刺激值计算的光谱预处理方法,包括:
步骤1,将参考标准放置在反射测量支架上,采用光纤光谱仪对参考标准进行光谱响应数据采集,将参考标准替换为待测样品,采用光纤光谱仪对待测样品进行光谱响应数据采集,将可见光截止档板放置在反射测量支架与光纤光谱仪之间,采用光纤光谱仪对暗噪声进行光谱响应数据采集,得到参考标准光谱响应数据、待测样品光谱响应数据和暗噪声的光谱响应数据;
步骤2,在每个波长采样点分别计算参考标准光谱响应数据的标准差、待测样品光谱响应数据的标准差和暗噪声的光谱响应数据的标准差,在颜色三刺激值计算波长间隔内,分别根据参考标准光谱响应数据的标准差、待测样品光谱响应数据的标准差和暗噪声的光谱响应数据的标准差,设定异常值剔除准则和方法,并剔除参考标准光谱响应数据的标准差、待测样品光谱响应数据的标准差和暗噪声的光谱响应数据的标准差的异常值相对应参考标准光谱响应数据、待测样品光谱响应数据和暗噪声的光谱响应数据;
步骤3,分别计算剔除异常值后的参考标准光谱响应数据的平均值、剔除异常值后的待测样品光谱响应数据的平均值和剔除异常值后的暗噪声的光谱响应数据的平均值;
步骤4,计算剔除异常值后的参考标准光谱响应数据的平均值的一阶导数,在颜色三刺激值计算波长间隔内,剔除一阶导数最大值对应的两个波长点λk和λk-1的光谱响应数据,得到剔除后的参考标准光谱响应数据的平均值;
步骤5,计算剔除异常值后的待测样品光谱响应数据的平均值的一阶导数,在颜色三刺激值计算波长间隔内,剔除一阶导数最大值对应的两个波长点λk和λk-1的光谱响应数据,得到剔除后的待测样品光谱响应数据的平均值;
步骤6,计算剔除异常值后的暗噪声的光谱响应数据的平均值的一阶导数,在颜色三刺激值计算波长间隔内,剔除一阶导数最大值对应的两个波长点λk和λk-1的光谱响应数据,得到剔除后的暗噪声的光谱响应数据的平均值;
步骤7,采用光谱平滑方法对剔除后的参考标准光谱响应数据的平均值、剔除后的待测样品光谱响应数据的平均值和剔除后的暗噪声的光谱响应数据的平均值进行去噪处理,得到去噪后的参考标准光谱响应数据的平均值、去噪后的待测样品光谱响应数据的平均值和去噪后的暗噪声的光谱响应数据的平均值;
步骤8,在可见光范围内,以颜色三刺激值计算波长间隔为插值间隔,采用插值算法将去噪后的参考标准光谱响应数据的平均值、去噪后的待测样品光谱响应数据的平均值和去噪后的暗噪声的光谱响应数据的平均值进行转换,得到转换后的参考标准光谱响应数据的平均值、转换后的待测样品光谱响应数据的平均值和转换后的暗噪声的光谱响应数据的平均值;
步骤9,根据转换后的参考标准光谱响应数据的平均值、转换后的待测样品光谱响应数据的平均值和转换后的暗噪声的光谱响应数据的平均值,计算待测样品的光谱反射率;
步骤10,根据计算出的待测样品的光谱反射率计算待测样品的颜色三刺激值。
其中,所述步骤1具体包括:
采用光纤光谱仪分别对参考标准、待测样品和暗噪声进行光谱响应数据采集,采集次数不小于6次,得到参考标准光谱响应数据Vr(λ)、待测样品光谱响应数据Vs(λ)和暗噪声的光谱响应数据Vd(λ),光纤光谱仪的积分时间为参考标准光谱响应数据达到最大的积分时间,且光谱采样波长间隔小于颜色三刺激值计算波长间隔的十分之一。
其中,所述步骤2具体包括:
在每个波长采样点计算参考标准光谱响应数据的标准差sr(λ),如下所示:
其中,s
r(λ)表示参考标准光谱响应数据的标准差,n表示测量次数,i表示第i次测量,λ表示波长,V
r-i(λ)表示第i次测量的参考标准光谱响应数据,
表示参考标准光谱响应数据的平均值;
在颜色三刺激值计算波长间隔内,查找参考标准光谱响应数据的标准差sr(λ)最大的波长点λmax_r,并剔除异常波长点λmax_r的参考标准光谱响应数据;
在每个波长采样点计算待测样品光谱响应数据的标准差ss(λ),如下所示:
其中,s
s(λ)表示待测样品光谱响应数据的标准差,n表示测量次数,i表示第i次测量,λ表示波长,V
s-i(λ)表示第i次测量的待测样品光谱响应数据,
表示待测样品光谱响应数据的平均值;
在颜色三刺激值计算波长间隔内,查找待测样品光谱响应数据的标准差ss(λ)最大的波长点λmax_s,并剔除异常波长点λmax_s的待测样品光谱响应数据;
在每个波长采样点计算暗噪声的光谱响应数据的标准差sd(λ),如下所示:
其中,s
d(λ)表示暗噪声的光谱响应数据的标准差,n表示测量次数,i表示第i次测量,λ表示波长,V
d-i(λ)表示第i次测量的暗噪声的光谱响应数据,
表示暗噪声的光谱响应数据的平均值;
在颜色三刺激值计算波长间隔内,查找暗噪声的光谱响应数据的标准差sd(λ)最大的波长点λmax_d,并剔除异常波长点λmax_d的暗噪声的光谱响应数据。
其中,所述步骤3具体包括:
分别计算剔除异常值后的参考标准光谱响应数据的平均值
剔除异常值后的待测样品光谱响应数据的平均值/>
和剔除异常值后的暗噪声的光谱响应数据的平均值/>
其中,所述步骤4、所述步骤5和所述步骤6具体包括:
计算剔除异常值后的参考标准光谱响应数据的平均值
的一阶导数d
r(λ),如下所示:
其中,d
r(λ
j)表示剔除异常值后的参考标准光谱响应数据的平均值的一阶导数,
表示λ
j波长点参考标准光谱响应数据的平均值,/>
表示λ
j-1波长点参考标准光谱响应数据的平均值,λ表示波长,λ
j表示第j个波长点,λ
j-1表示第j-1个波长点;
在颜色三刺激值计算波长间隔内,将剔除异常值后的参考标准光谱响应数据的平均值
的一阶导数d
r(λ)中最大值对应的两个波长点λ
k和λ
k-1的光谱响应数据剔除,得到剔除后的参考标准光谱响应数据的平均值/>
计算剔除异常值后的待测样品光谱响应数据的平均值
的一阶导数d
s(λ),如下所示:/>
其中,d
s(λ
j)表示剔除异常值后的待测样品光谱响应数据的平均值的一阶导数,
表示λ
j波长点剔除异常值后的待测样品光谱响应数据的平均值,/>
表示λ
j-1波长点剔除异常值后的待测样品光谱响应数据的平均值,λ表示波长,λ
j表示第j个波长点,λ
j-1表示第j-1个波长点;
在颜色三刺激值计算波长间隔内,将剔除异常值后的待测样品光谱响应数据的平均值
的一阶导数d
s(λ)中最大值对应的两个波长点λ
k和λ
k-1的光谱响应数据剔除,得到剔除后的待测样品光谱响应数据的平均值/>
计算剔除异常值后的暗噪声的光谱响应数据的平均值
的一阶导数d
d(λ),如下所示:
其中,d
d(λ
j)表示剔除异常值后的暗噪声的光谱响应数据的平均值的一阶导数,
表示λ
j波长点剔除异常值后的暗噪声的光谱响应数据的平均值,/>
表示λ
j-1波长点剔除异常值后的暗噪声的光谱响应数据的平均值,λ表示波长,λ
j表示第j个波长点,λ
j-1表示第j-1个波长点;
在颜色三刺激值计算波长间隔内,将剔除异常值后的暗噪声的光谱响应数据的平均值
的一阶导数d
d(λ)中最大值对应的两个波长点λ
k和λ
k-1的光谱响应数据剔除,得到剔除后的暗噪声的光谱响应数据的平均值/>
其中,所述步骤7具体包括:
采用光谱平滑方法对剔除后的参考标准光谱响应数据的平均值
剔除后的待测样品光谱响应数据的平均值/>
和剔除后的暗噪声的光谱响应数据的平均值/>
进行去噪处理,得到去噪后的参考标准光谱响应数据的平均值/>
去噪后的待测样品光谱响应数据的平均值/>
和去噪后的暗噪声的光谱响应数据的平均值/>
其中,所述步骤8具体包括:
在可见光范围内,以颜色三刺激值计算波长间隔为插值间隔,采用插值算法将去噪后的参考标准光谱响应数据的平均值
去噪后的待测样品光谱响应数据的平均值
和去噪后的暗噪声的光谱响应数据的平均值/>
进行转换,得到相同的波长范围和相同波长间隔的转换后的参考标准光谱响应数据的平均值/>
转换后的待测样品光谱响应数据的平均值/>
和转换后的暗噪声的光谱响应数据的平均值/>
其中,所述步骤9具体包括:
根据得到的转换后的参考标准光谱响应数据的平均值
转换后的待测样品光谱响应数据的平均值/>
和转换后的暗噪声的光谱响应数据的平均值/>
计算待测样品的光谱反射率ρ
s(λ),如下所示:
其中,ρ
s(λ)表示待测样品的光谱反射率,ρ
r(λ)表示参考标准的光谱反射率,λ表示波长,
表示转换后的参考标准光谱响应数据的平均值,/>
表示转换后的待测样品光谱响应数据的平均值,/>
表示转换后的暗噪声的光谱响应数据的平均值。
其中,所述步骤10具体包括:
根据得到的待测样品的光谱反射率ρs(λ),计算待测样品的颜色三刺激值,如下所示:
其中,ρ
s(λ)表示待测样品的光谱反射率,S(λ)表示CIE规定的标准照明体相对光谱功率分布,
为CIE规定的标准色度观察者数据,Δλ表示波长间隔,k表示归一化系数,λ表示波长。
本发明的上述方案有如下的有益效果:
本发明的上述实施例所述的用于颜色三刺激值计算的颜色测量平台及光谱预处理方法,可以消除光谱异常峰值和高频随机噪声,可以剔除光谱异常峰值,适合对无明显特征峰值,曲线变化趋势比较平缓的光谱数据进行预处理,可提高光纤光谱仪测量颜色的准确性和重复性。
具体实施方式
为使本发明要解决的技术问题、技术方案和优点更加清楚,下面将结合附图及具体实施例进行详细描述。
本发明针对现有的光谱预处理方法不适合处理无明显特征峰值,曲线变化趋势比较平缓的,用于颜色三刺激值计算的光谱响应数据的问题,提供了一种用于颜色三刺激值计算的颜色测量平台及光谱预处理方法。
如图1至图6所示,本发明的实施例提供了一种用于颜色三刺激值计算的颜色测量平台,包括:反射测量支架1;待测样品2,所述待测样品2放置在所述反射测量支架1的顶部;第一准直镜3,所述第一准直镜3设置在所述待测样品2的正上方;第二准直镜4,所述第二准直镜4倾斜设置在所述待测样品2的侧上方,所述第二准直镜4与所述第一准直镜3的之间的夹角呈45°;光纤光源,所述光纤光源通过入射光纤8与所述第二准直镜4连接;光纤光谱仪6,所述光纤光谱仪6的第一端通过出射光纤9与所述第一准直镜3连接;计算机7,所述计算机7与所述光纤光谱仪6的第二端电连接。
本发明的上述实施例所述的用于颜色三刺激值计算的颜色测量平台及光谱预处理方法,所述光纤光源5采用光纤卤素灯,采用45°/0°反射测量支架1、可见光所光纤光谱仪6和装有光谱采集和分析软件的所述计算机7,搭建所述用于颜色三刺激值计算的颜色测量平台。
本发明的实施例还提供了一种用于颜色三刺激值计算的光谱预处理方法,包括:步骤1,将参考标准放置在反射测量支架1上,采用光纤光谱仪6对参考标准进行光谱响应数据采集,将参考标准替换为待测样品2,采用光纤光谱仪6对待测样品2进行光谱响应数据采集,将可见光截止档板放置在反射测量支架1与光纤光谱仪6之间,采用光纤光谱仪6对暗噪声进行光谱响应数据采集,得到参考标准光谱响应数据、待测样品2光谱响应数据和暗噪声的光谱响应数据;步骤2,在每个波长采样点分别计算参考标准光谱响应数据的标准差、待测样品2光谱响应数据的标准差和暗噪声的光谱响应数据的标准差,在颜色三刺激值计算波长间隔内,分别根据参考标准光谱响应数据的标准差、待测样品2光谱响应数据的标准差和暗噪声的光谱响应数据的标准差,设定异常值剔除准则和方法,并剔除参考标准光谱响应数据的标准差、待测样品2光谱响应数据的标准差和暗噪声的光谱响应数据的标准差的异常值相对应参考标准光谱响应数据、待测样品2光谱响应数据和暗噪声的光谱响应数据;步骤3,分别计算剔除异常值后的参考标准光谱响应数据的平均值、剔除异常值后的待测样品2光谱响应数据的平均值和剔除异常值后的暗噪声的光谱响应数据的平均值;步骤4,计算剔除异常值后的参考标准光谱响应数据的平均值的一阶导数,在颜色三刺激值计算波长间隔内,剔除一阶导数最大值对应的两个波长点λk和λk-1的光谱响应数据,得到剔除后的参考标准光谱响应数据的平均值;步骤5,计算剔除异常值后的待测样品2光谱响应数据的平均值的一阶导数,在颜色三刺激值计算波长间隔内,剔除一阶导数最大值对应的两个波长点λk和λk-1的光谱响应数据,得到剔除后的待测样品2光谱响应数据的平均值;步骤6,计算剔除异常值后的暗噪声的光谱响应数据的平均值的一阶导数,在颜色三刺激值计算波长间隔内,剔除一阶导数最大值对应的两个波长点λk和λk-1的光谱响应数据,得到剔除后的暗噪声的光谱响应数据的平均值;步骤7,采用光谱平滑方法对剔除后的参考标准光谱响应数据的平均值、剔除后的待测样品2光谱响应数据的平均值和剔除后的暗噪声的光谱响应数据的平均值进行去噪处理,得到去噪后的参考标准光谱响应数据的平均值、去噪后的待测样品2光谱响应数据的平均值和去噪后的暗噪声的光谱响应数据的平均值;步骤8,在可见光范围内,以颜色三刺激值计算波长间隔为插值间隔,采用插值算法将去噪后的参考标准光谱响应数据的平均值、去噪后的待测样品2光谱响应数据的平均值和去噪后的暗噪声的光谱响应数据的平均值进行转换,得到转换后的参考标准光谱响应数据的平均值、转换后的待测样品2光谱响应数据的平均值和转换后的暗噪声的光谱响应数据的平均值;步骤9,根据转换后的参考标准光谱响应数据的平均值、转换后的待测样品2光谱响应数据的平均值和转换后的暗噪声的光谱响应数据的平均值,计算待测样品2的光谱反射率;步骤10,根据计算出的待测样品2的光谱反射率计算待测样品2的颜色三刺激值。
其中,所述步骤1具体包括:采用光纤光谱仪6分别对参考标准、待测样品2和暗噪声进行光谱响应数据采集,采集次数不小于6次,得到参考标准光谱响应数据Vr(λ)、待测样品2光谱响应数据Vs(λ)和暗噪声的光谱响应数据Vd(λ),光纤光谱仪6的积分时间为参考标准光谱响应数据达到最大的积分时间,且光谱采样波长间隔小于颜色三刺激值计算波长间隔的十分之一。
本发明的上述实施例所述的用于颜色三刺激值计算的颜色测量平台及光谱预处理方法,确定测量参数:所述光纤光谱仪6采样波长范围为380nm-780nm,采样波长间隔为0.5nm,采样次数为6次;参考标准采用BCRA标准白色陶瓷板,所述待测样品2采用BCRA标准红色陶瓷板,在所述反射测量支架1上放置参考标准,采用所述光纤光谱仪6测量参考标准的光谱响应数据,确定参考标准光谱响应值达到最大的积分时间为155毫秒,并设置155毫秒为光纤光谱仪6积分时间。在所述反射测量支架1上放置参考标准,以155毫秒为光纤光谱仪6积分时间,测量参考标准6次,获得参考标准光谱响应数据Vr(λ),如图3所示。在所述反射测量支架1上放置所述待测样品2,以155毫秒为光纤光谱仪6积分时间,测量6次,获得所述待测样品2光谱响应数据Vs(λ),如图3所示。在所述反射测量支架1与所述光纤光谱仪6间放置可见光截止档板,以155毫秒为光纤光谱仪6积分时间,测量6次,获得暗噪声光谱响应数据Vd(λ),如图3所示。
其中,所述步骤2具体包括:在每个波长采样点计算参考标准光谱响应数据的标准差sr(λ),如下所示:
其中,s
r(λ)表示参考标准光谱响应数据的标准差,n表示测量次数,i表示第i次测量,λ表示波长,V
r-i(λ)表示第i次测量的参考标准光谱响应数据,
表示参考标准光谱响应数据的平均值;
在颜色三刺激值计算波长间隔内,查找参考标准光谱响应数据的标准差sr(λ)最大的波长点λmax_r,并剔除异常波长点λmax_r的参考标准光谱响应数据;
在每个波长采样点计算待测样品2光谱响应数据的标准差ss(λ),如下所示:
其中,s
s(λ)表示待测样品2光谱响应数据的标准差,n表示测量次数,i表示第i次测量,λ表示波长,V
s-i(λ)表示第i次测量的待测样品2光谱响应数据,
表示待测样品2光谱响应数据的平均值;
在颜色三刺激值计算波长间隔内,查找待测样品2光谱响应数据的标准差ss(λ)最大的波长点λmax_s,并剔除异常波长点λmax_s的待测样品2光谱响应数据;
在每个波长采样点计算暗噪声的光谱响应数据的标准差sd(λ),如下所示:
其中,s
d(λ)表示暗噪声的光谱响应数据的标准差,n表示测量次数,i表示第i次测量,λ表示波长,V
d-i(λ)表示第i次测量的暗噪声的光谱响应数据,
表示暗噪声的光谱响应数据的平均值;
在颜色三刺激值计算波长间隔内,查找暗噪声的光谱响应数据的标准差sd(λ)最大的波长点λmax_d,并剔除异常波长点λmax_d的暗噪声的光谱响应数据。
本发明的上述实施例所述的用于颜色三刺激值计算的颜色测量平台及光谱预处理方法,计算参考标准光谱响应数据的标准差sr(λ),如图4所示;在颜色三刺激值计算波长间隔内,将各波长点的标准差从大到小排序,并剔除前五位波长点的样品光谱响应数据。
表1600nm-610nm间参考标准光谱响应数据的标准差
例如表1所示为颜色三刺激值计算波长间隔600nm-610nm之间的从大到小排序的标准差,根据上述的异常值剔除准则和方法,判定波长点607nm,607.5nm,609nm,609.5nm,605nm的光谱响应值异常,并剔除这五个波长点的光谱响应值,其他波长点的数据采用相同方法剔除异常;计算待测样品2光谱响应数据的标准差ss(λ),在颜色三刺激值计算波长间隔内,将各波长点的标准差从大到小排序,并剔除前五位波长点的待测样品2光谱响应数据中的异常值;计算暗噪声的光谱响应数据的标准差sd(λ),在颜色三刺激值计算波长间隔内,将各波长点的标准差从大到小排序,并剔除前五位波长点的暗噪声的光谱响应数据中的异常值。
其中,所述步骤3具体包括:分别计算剔除异常值后的参考标准光谱响应数据的平均值
剔除异常值后的待测样品2光谱响应数据的平均值/>
和剔除异常值后的暗噪声的光谱响应数据的平均值/>
其中,所述步骤4、所述步骤5和所述步骤6具体包括:计算剔除异常值后的参考标准光谱响应数据的平均值
的一阶导数d
r(λ),如下所示:
其中,d
r(λ
j)表示剔除异常值后的参考标准光谱响应数据的平均值的一阶导数,
表示λ
j波长点参考标准光谱响应数据的平均值,/>
表示λ
j-1波长点参考标准光谱响应数据的平均值,λ表示波长,λ
j表示第j个波长点,λ
j-1表示第j-1个波长点;
在颜色三刺激值计算波长间隔内,将剔除异常值后的参考标准光谱响应数据的平均值
的一阶导数d
r(λ)中最大值对应的两个波长点λ
k和λ
k-1的光谱响应数据剔除,得到剔除后的参考标准光谱响应数据的平均值/>
计算剔除异常值后的待测样品2光谱响应数据的平均值
的一阶导数d
s(λ),如下所示:
其中,d
s(λ
j)表示剔除异常值后的待测样品光谱响应数据的平均值的一阶导数,
表示λ
j波长点剔除异常值后的待测样品光谱响应数据的平均值,/>
表示λ
j-1波长点剔除异常值后的待测样品光谱响应数据的平均值,λ表示波长,λ
j表示第j个波长点,λ
j-1表示第j-1个波长点;
在颜色三刺激值计算波长间隔内,将剔除异常值后的待测样品2光谱响应数据的平均值
的一阶导数d
s(λ)中最大值对应的两个波长点λ
k和λ
k-1的光谱响应数据剔除,得到剔除后的待测样品2光谱响应数据的平均值/>
计算剔除异常值后的暗噪声的光谱响应数据的平均值
的一阶导数d
d(λ),如下所示:
其中,d
d(λ
j)表示剔除异常值后的暗噪声的光谱响应数据的平均值的一阶导数,
表示λ
j波长点剔除异常值后的暗噪声的光谱响应数据的平均值,/>
表示λ
j-1波长点剔除异常值后的暗噪声的光谱响应数据的平均值,λ表示波长,λ
j表示第j个波长点,λ
j-1表示第j-1个波长点;
在颜色三刺激值计算波长间隔内,将剔除异常值后的暗噪声的光谱响应数据的平均值
的一阶导数d
d(λ)中最大值对应的两个波长点λ
k和λ
k-1的光谱响应数据剔除,得到剔除后的暗噪声的光谱响应数据的平均值/>
本发明的上述实施例所述的用于颜色三刺激值计算的颜色测量平台及光谱预处理方法,计算剔除异常值后参考标准光谱响应数据的平均值
并根据式(6)计算剔除异常值后的参考标准光谱响应数据的平均值/>
的一阶导数d
r(λ),如图5所示;在颜色三刺激值计算波长间隔内,将剔除异常值后的参考标准光谱响应数据的平均值/>
的一阶导数值的最大值对应的两个波长点λ
k和λ
k-1的光谱响应数据剔除,得到剔除后的参考标准光谱响应数据的平均值/>
表2 600nm-610nm间参考标准光谱响应数据的一阶导数
如表2所示,波长为600nm-610nm之间的从大到小排序的剔除异常值后的参考标准光谱响应数据的平均值
的一阶导数,波长点602.5nm的一阶导数最大,因此剔除一阶导数最大值对应的602.5nm和602nm的光谱响应数据,其他波长点的数据采用相同方法剔除异常;计算剔除异常值后的待测样品2光谱响应数据的平均值/>
并根据式(4)计算剔除异常值后的待测样品2光谱响应数据的平均值/>
的一阶导数d
s(λ),如图5所示;在颜色三刺激值计算波长间隔内,将剔除异常值后的待测样品2光谱响应数据的平均值/>
的一阶导数值的最大值对应的两个波长点λ
k和λ
k-1的光谱响应数据剔除,得到剔除后的待测样品2光谱响应数据的平均值/>
计算剔除异常值后的暗噪声的光谱响应数据的平均值/>
并根据式(4)计算剔除异常值后的暗噪声的光谱响应数据的平均值/>
的一阶导数d
d(λ),如图5所示;在颜色三刺激值计算波长间隔内,将剔除异常值后的暗噪声的光谱响应数据的平均值/>
的一阶导数值的最大值对应的两个波长点λ
k和λ
k-1的光谱响应数据剔除,得到剔除后的待测样品2光谱响应数据的平均值/>
其中,所述步骤7具体包括:采用光谱平滑方法对剔除后的参考标准光谱响应数据的平均值
剔除后的待测样品2光谱响应数据的平均值/>
和剔除后的暗噪声的光谱响应数据的平均值/>
进行去噪处理,得到去噪后的参考标准光谱响应数据的平均值/>
去噪后的待测样品2光谱响应数据的平均值/>
和去噪后的暗噪声的光谱响应数据的平均值/>
其中,所述步骤8具体包括:在可见光范围内,以颜色三刺激值计算波长间隔为插值间隔,采用插值算法将去噪后的参考标准光谱响应数据的平均值
去噪后的待测样品2光谱响应数据的平均值/>
和去噪后的暗噪声的光谱响应数据的平均值/>
进行转换,得到相同的波长范围和相同波长间隔的转换后的参考标准光谱响应数据的平均值
转换后的待测样品2光谱响应数据的平均值/>
和转换后的暗噪声的光谱响应数据的平均值/>
本发明的上述实施例所述的用于颜色三刺激值计算的颜色测量平台及光谱预处理方法,采用插值算法将去噪后的参考标准光谱响应数据的平均值
去噪后的待测样品2光谱响应数据的平均值/>
和去噪后的暗噪声的光谱响应数据的平均值/>
转换成波长范围为380nm~780nm,波长间隔为5nm的转换后的参考标准光谱响应数据的平均值/>
转换后的待测样品2光谱响应数据的平均值/>
和转换后的暗噪声的光谱响应数据的平均值/>
如图6所示。
其中,所述步骤9具体包括:根据得到的转换后的参考标准光谱响应数据的平均值
转换后的待测样品2光谱响应数据的平均值/>
和转换后的暗噪声的光谱响应数据的平均值/>
计算待测样品的光谱反射率ρ
s(λ),如下所示:
其中,ρ
s(λ)表示待测样品的光谱反射率,ρ
r(λ)表示参考标准的光谱反射率,λ表示波长,
表示转换后的参考标准光谱响应数据的平均值,/>
表示转换后的待测样品光谱响应数据的平均值,/>
表示转换后的暗噪声的光谱响应数据的平均值。
其中,所述步骤10具体包括:根据得到的待测样品的光谱反射率ρs(λ),计算待测样品的颜色三刺激值,如下所示:
其中,ρ
s(λ)表示待测样品的光谱反射率,S(λ)表示CIE规定的标准照明体相对光谱功率分布,
为CIE规定的标准色度观察者数据,Δλ表示波长间隔,k表示归一化系数,λ表示波长。
本发明的上述实施例所述的用于颜色三刺激值计算的颜色测量平台及光谱预处理方法,重复步骤2至步骤10,对参考标准、待测样品2和暗噪声重复测量四次,并计算颜色三刺激值X值、Y值和Z值;
表3预处理前后光谱响应数据计算X值、Y值和Z值对比
从表3可知,预处理后的光谱响应数据计算X值、Y值和Z值的标准差明显低于预处理前原始光谱响应数据计算X值、Y值和Z值的标准差,说明所述的用于颜色三刺激值计算的颜色测量平台及光谱预处理方法可以有效消除测量过程中的噪声影响,提高测量准确性和重复性。
本发明的上述实施例所述的用于颜色三刺激值计算的颜色测量平台及光谱预处理方法,根据CIE规定的颜色测量几何条件,采用所述光纤光源5、所述反射测量支架1、所述光纤光谱仪6和带光谱采集及分析软件的所述计算机7搭建所述用于颜色三刺激值计算的颜色测量平台;确定测量参数:所述光纤光谱仪6采集光谱响应数据次数不少于6次,采用所述光纤光谱仪6测量参考标准光谱响应数据,确定参考标准光谱响应数据达到最大的积分时间为155毫秒,并设置155毫秒为所述光纤光谱仪6积分时间,所述光纤光谱仪6积分时间是参考标准光谱响应数据达到最大的积分时间,光谱采样波长范围为380nm-780nm,光谱采样波长间隔为0.5nm,光谱采样波长间隔应小于颜色三刺激值计算波长间隔的十分之一。采用所述用于颜色三刺激值计算的颜色测量平台,分别测量参考标准、所述待测样品2和暗噪声的光谱响应数据。确定基于单个波长点光谱响应数据的标准差和剔除光谱响应数据异常值的准则和方法,计算每个波长点多次测量光谱响应数据的标准差,根据异常值准则和方法,查找光谱响应数据异常的波长点,并予以剔除。计算剔除异常值后的参考标准、所述待测样品2和暗噪声的光谱响应数据的平均值。计算参考标准、所述待测样品2和暗噪声的光谱响应数据的平均值的一阶导数,剔除参考标准、所述待测样品2和暗噪声的光谱响应数据的平均值的一阶导数中最大值对应的光谱响应数据的平均值,采用光谱平滑方法对剔除后的参考标准、所述待测样品2和暗噪声的光谱响应数据的平均值进行去噪处理,采用插值算法将去噪后的参考标准、所述待测样品2和暗噪声的光谱响应数据的平均值转换为相同的波长范围和相同波长间隔的转换后的参考标准、所述待测样品2和暗噪声的光谱响应数据的平均值,根据转换后的参考标准、所述待测样品2和暗噪声的光谱响应数据的平均值计算所述待测样品2的光谱反射率,根据所述待测样品2的光谱反射率计算所述待测样品2的颜色三刺激值。
本发明的上述实施例所述的用于颜色三刺激值计算的颜色测量平台及光谱预处理方法,提高了所述光纤光谱仪6测量颜色的准确性和重复性,可以剔除光谱异常峰值和高频随机噪声,适合对无明显特征峰值,曲线变化趋势比较平缓的光谱数据进行预处理。
以上所述是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明所述原理的前提下,还可以作出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。