CN113049104B - 三维物体色度值计算方法及装置、三维物体色度值计算系统 - Google Patents
三维物体色度值计算方法及装置、三维物体色度值计算系统 Download PDFInfo
- Publication number
- CN113049104B CN113049104B CN202110242021.4A CN202110242021A CN113049104B CN 113049104 B CN113049104 B CN 113049104B CN 202110242021 A CN202110242021 A CN 202110242021A CN 113049104 B CN113049104 B CN 113049104B
- Authority
- CN
- China
- Prior art keywords
- dimensional object
- value
- model
- value calculation
- fitted
- 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.)
- Active
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 166
- 230000003595 spectral effect Effects 0.000 claims abstract description 41
- 238000005286 illumination Methods 0.000 claims abstract description 23
- 238000002310 reflectometry Methods 0.000 claims abstract description 15
- 238000000034 method Methods 0.000 claims description 37
- 230000006870 function Effects 0.000 claims description 17
- 238000005259 measurement Methods 0.000 claims description 15
- 239000003086 colorant Substances 0.000 claims description 14
- 238000004590 computer program Methods 0.000 claims description 14
- 238000012216 screening Methods 0.000 claims description 14
- 238000012545 processing Methods 0.000 claims description 12
- 238000001514 detection method Methods 0.000 claims description 6
- 238000010146 3D printing Methods 0.000 claims description 4
- 238000005516 engineering process Methods 0.000 claims description 4
- 238000000691 measurement method Methods 0.000 claims description 3
- 239000007787 solid Substances 0.000 claims description 3
- 238000007639 printing Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 7
- 230000004456 color vision Effects 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 4
- 230000004075 alteration Effects 0.000 description 3
- 230000000007 visual effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000002798 spectrophotometry method Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
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/46—Measurement of colour; Colour measuring devices, e.g. colorimeters
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
本申请涉及三维物体色度值计算方法及装置、三维物体色度值计算系统，其中，三维物体色度值计算方法包括：获取多个三维物体模型以及同一色系的多个平面模型；获取置于标准观察箱中的三维物体模型与目标平面模型的光谱反射率以及标准观察箱中的照明光源的相对光谱能量分布值；根据计算得到的三维物体模型的第一计算色度值、目标平面模型的第二计算色度值构建得到线性拟合函数，线性拟合函数为拟合后的色度值计算公式；根据拟合后的色度值计算公式对待评价三维物体进行色度值计算。本申请提供的三维物体色度值计算方法及装置、系统，能够提高三维物体的色度值的准确度。
Description
技术领域
本申请涉及三维物体颜色评价技术领域，尤其涉及三维物体色度值计算方法及装置、三维物体色度值计算系统。
背景技术
平面物体颜色测量包括光源颜色测量与物体颜色测量两大类，其中物体颜色测量分为接触式测量和非接触式测量。现有平面物体颜色测量方法有目视法、光电积分法和分光光度法三种。专利CN107421468A提出了无需标记点的彩色三维扫描系统，将编码图案投影到扫描物体上，用两个工业相机采集投影编码图案进行几何三维模型重建，通过一个彩色相机获取三维物体的彩色照片；这种方法需要用到标定算法，将三维数据与颜色数据配对起来，过程复杂，且由于采样点的色彩是由扫描系统提供的环境光照射下呈现的颜色，获得的物体表面色彩与自然光下视觉效果不会完全一样。目前在工业界，对三维物体的颜色测量，较多沿用对二维平面物体颜色测量的方法。然而，三维物体由于具有不同于二维平面物体的形状，在散射光或定向光源照明下，其外观颜色受光源照明角度、物体形状、半透明度和纹理等因素的影响，即便是具有同一色度值的颜色也会出现观察者颜色感觉不一致的现象。因此，现有的三维物体的色度值计算方法准确度低。
发明内容
本申请提供三维物体色度值计算方法及装置、三维物体色度值计算系统，以提高三维物体的色度值计算准确度。
第一方面，本申请提供一种三维物体色度值计算方法，所述方法包括：
获取多个三维物体模型；基于所述三维物体模型的颜色获取同一色系的多个平面模型；获取置于标准观察箱中的三维物体模型与目标平面模型的光谱反射率以及标准观察箱中的照明光源的相对光谱能量分布值；其中，所述目标平面模型为所述同一色系的多个平面模型中与所述三维物体模型的颜色接近度最高的一个平面模型；根据所述光谱反射率及所述光谱能量分布值分别计算得到所述三维物体模型的第一计算色度值、所述目标平面模型的第二计算色度值；基于多个所述三维物体模型的第一计算色度值与多个所述目标平面模型的第二计算色度值构建得到线性拟合函数，所述线性拟合函数为拟合后的色度值计算公式；根据拟合后的色度值计算公式对待评价三维物体进行色度值计算。
结合第一方面，在一种可行的实施方式中，所述第一计算色度值和第二计算色度值包括明度值、饱和度值及色调角值。
结合第一方面，在一种可行的实施方式中，所述方法满足下列特征(1)至(7)中的至少一种：
(1)所述三维物体模型、所述平面模型均采用三维打印技术得到；
(2)所述三维物体模型的最小维度大于或等于4cm；
(3)所述三维物体模型为规则形状且单色的立体模型；
(4)同一色系的所述三维物体模型的数量为m，m为大于或等于4的整数；
(5)所述多个三维物体模型的颜色选自国际照明委员会推荐的颜色中心的颜色进行打印，至少包括5种不同色系，且同一色系的多个所述三维物体模型的色度值不同；
(6)与每个所述三维物体模型对应的同一色系的平面模型的数量为n，n为大于或等于10的整数；
(7)同一色系的三维物体模型、平面模型的维度一致。
结合第一方面，在一种可行的实施方式中，所述多个三维物体模型的颜色分别选自灰色、红色、黄色、绿色和蓝色进行打印，其中，灰色的
色度值为(62.0，0.0，0.0)，红色的
色度值为(44.0，37.0，23.0)，黄色的
色度值为(87.0，-7.0，47.0)，绿色的
色度值为(56.0，-32.0，0.0)，蓝色的
色度值为(36.0，5.0，-31.0)。
结合第一方面，在一种可行的实施方式中，所述平面模型的厚度小于或等于1mm，且所述三维物体模型在所述平面模型上的正投影面积等于所述平面模型的面积。
结合第一方面，在一种可行的实施方式中，在所述获取置于标准观察箱中的三维物体模型与目标平面模型的光谱反射率以及标准观察箱中的照明光源的相对光谱能量分布值之前，所述方法还包括：
对多个平面模型进行筛选，得到有效的平面模型；基于每个所述三维物体模型的颜色从所述同一色系的多个有效的平面模型中筛选得到目标平面模型，其中，所述目标平面模型与同一色系的所述三维物体模型的颜色接近度最高。
结合第一方面，在一种可行的实施方式中，所述对多个平面模型进行筛选，得到有效的平面模型具体包括：
测量得到每个所述三维物体模型的第一测量色度值；测量得到每个所述平面模型的第二测量色度值；基于所述第一测量色度值和第二测量色度值获取与所述三维物体模型的颜色属于同一色系的多个有效的平面模型。
结合第一方面，在一种可行的实施方式中，所述测量得到每个所述三维物体模型的第一测量色度值，包括：测量每个所述三维物体模型至少5个不同位置处的色度值，并计算所述色度值的算术平均值，得到所述三维物体模型的第一测量色度值。
结合第一方面，在一种可行的实施方式中，所述至少5个不同位置位于所述三维物体模型同一平面或同一弧面上。
结合第一方面，在一种可行的实施方式中，所述测量得到每个所述平面模型的第二测量色度值，包括：测量每个所述平面模型至少5个不同位置处的色度值，并计算所述色度值的算术平均值，得到所述平面模型的第二测量色度值。
结合第一方面，在一种可行的实施方式中，所述三维物体模型为立方体，所述立方体包括上表面，所述拟合后的上表面色度值计算公式包括明度值计算公式、饱和度值计算公式及色调角值计算公式；其中，x为待评价三维物体的初始色度值，y为待评价三维物体拟合后的色度值；
所述明度值计算公式如式(A-1)为：yL上＝1.0797xL上-5.6981 (A-1)；
所述饱和度值计算公式如式(B-1)为：yc上＝1.0077xC上+0.666 (B-1)；
所述色调角值计算公式如式(C-1)为：yh上＝0.9853xh上+1.587 (C-1)。
结合第一方面，在一种可行的实施方式中，所述三维物体模型为立方体，所述立方体包括前表面，所述拟合后的前表面色度值计算公式包括明度值计算公式、饱和度值计算公式及色调角值计算公式；其中，x为待评价三维物体的初始色度值，y为待评价三维物体拟合后的色度值；
所述明度值计算公式如式(A-2)为：yL前＝1.0304xL前-10.3 (A-2)；
所述饱和度值计算公式如式(B-2)为：yc前＝1.0082xC前-1.7921 (B-2)；
所述色调角值计算公式如式(C-2)为：yh前＝0.9582xh前+3.8288 (C-2)。
结合第一方面，在一种可行的实施方式中，根据拟合后的色度值计算公式对待评价三维物体进行色度值计算，包括：
获取待评价三维物体，测量得到所述待评价三维物体上表面和/或前表面的测量色度值，所述上表面和/或前表面的测量色度值包括明度值、饱和度值及色调角值；将所述上表面和/或前表面的明度值代入拟合后的上表面和/或前表面的明度值计算公式，得到拟合后的明度值；将所述上表面和/或前表面的饱和度值代入拟合后的上表面和/或前表面的饱和度值计算公式，得到拟合后的饱和度值；将所述上表面和/或前表面的色调角值代入拟合后的上表面和/或前表面的色调角值计算公式，得到拟合后的色调角值；根据拟合后的明度值、饱和度值及色调角值得到所述待评价三维物体拟合后的色度值。
第二方面，本申请实施例提供一种三维物体色度值计算装置，所述装置包括：
第一获取单元，用于获取多个三维物体模型；
第二获取单元，用于基于所述三维物体模型的颜色获取同一色系的多个平面模型；
第三获取单元，用于获取置于标准观察箱中的三维物体模型与目标平面模型的光谱反射率以及标准观察箱中的照明光源的相对光谱能量分布值；
第一计算单元，用于根据所述光谱反射率及所述光谱能量分布值分别计算得到所述三维物体模型的第一计算色度值、所述目标平面模型的第二计算色度值；其中，所述目标平面模型为所述同一色系的多个平面模型中与所述三维物体模型的颜色接近度最高的一个平面模型；
构建单元，用于基于多个所述三维物体模型的第一计算色度值与多个所述目标平面模型的第二计算色度值构建得到线性拟合函数，所述线性拟合函数为拟合后的色度值计算公式；
第二计算单元，用于根据拟合后的色度值计算公式对待评价三维物体进行色度值计算。
结合第二方面，在一种可行的实施方式中，所述第一计算色度值和第二计算色度值包括明度值、饱和度值及色调角值。
结合第二方面，在一种可行的实施方式中，所述三维物体模型为立方体，所述立方体包括上表面，所述拟合后的上表面色度值计算公式包括明度值计算公式、饱和度值计算公式及色调角值计算公式；其中，x为待评价三维物体的初始色度值，y为待评价三维物体拟合后的色度值；
所述明度值计算公式如式(A-1)为：yL上＝1.0797xL上-5.6981 (A-1)；
所述饱和度值计算公式如式(B-1)为：yc上＝1.0077xC上+0.666 (B-1)；
所述色调角值计算公式如式(C-1)为：yh上＝0.9853xh上+1.587 (C-1)。
结合第二方面，在一种可行的实施方式中，所述三维物体模型为立方体，所述立方体包括前表面，所述拟合后的前表面色度值计算公式包括明度值计算公式、饱和度值计算公式及色调角值计算公式；其中，x为待评价三维物体的初始色度值，y为待评价三维物体拟合后的色度值；
所述明度值计算公式如式(A-2)为：yL前＝1.0304xL前-10.3 (A-2)；
所述饱和度值计算公式如式(B-2)为：yc前＝1.0082xC前-1.7921 (B-2)；
所述色调角值计算公式如式(C-2)为：yh前＝0.9582xh前+3.8288 (C-2)。
结合第二方面，在一种可行的实施方式中，所述第二计算单元，包括：
获取子单元，用于获取待评价三维物体，测量得到所述待评价三维物体上表面和/或前表面的测量色度值，所述上表面和/或前表面的测量色度值包括明度值、饱和度值及色调角值；
第一处理子单元，用于将所述上表面和/或前表面的明度值代入拟合后的上表面和/或前表面的明度值计算公式，得到拟合后的明度值；
第二处理子单元，用于将所述上表面和/或前表面的饱和度值代入拟合后的上表面和/或前表面的饱和度值计算公式，得到拟合后的饱和度值；
第三处理子单元，用于将所述上表面和/或前表面的色调角值代入拟合后的上表面和/或前表面的色调角值计算公式，得到拟合后的色调角值；及
输出子单元，用于根据拟合后的明度值、饱和度值及色调角值得到所述待评价三维物体拟合后的色度值。
第三方面，本申请提供一种计算机非易失性存储介质，所述存储介质包括存储的程序，在所述程序运行时控制所述存储介质所在设备执行上述第一方面所述的三维物体色度值计算方法中拟合后的色度值计算公式。
第四方面，本申请提供一种计算机设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现上述第一方面所述的三维物体色度值计算方法中拟合后的色度值计算公式。
第五方面，本申请提供一种三维物体色度值计算系统，包括检测部件、存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述第一方面所述的三维物体色度值计算方法中拟合后的色度值计算公式。
有益效果：
在本方案中，通过利用与三维物体模型颜色最接近的目标平面模型，将两者置于同一标准观察箱中，得到三维物体模型与目标平面模型的计算色度值，通过多个三维物体模型的计算色度值以及多个目标平面模型的计算色度值来构建拟合函数，从而得到拟合后的色度值计算公式，使得通过拟合后的色度值计算公式能够提高得到的三维物体的色度值的准确度。
附图说明
下面结合附图和实施例对本申请进一步说明。
图1为本申请具体实施例中三维物体色度值计算方法流程示意图；
图2-1为本申请具体实施例中三维物体模型的上表面与目标平面模型的颜色L* 10的散点分布和数学拟合关系图；
图2-2为本申请具体实施例中三维物体模型的上表面与目标平面模型的颜色C* 10,ab的散点分布和数学拟合关系图；
图2-3为本申请具体实施例中三维物体模型的上表面与目标平面模型的颜色h* 10,ab的散点分布和数学拟合关系图；
图3-1为本申请实施例中三维物体模型的前表面与目标平面模型的颜色L* 10的散点分布和数学拟合关系图；
图3-2为本申请实施例中三维物体模型的前表面与目标平面模型的颜色C* 10,ab的散点分布和数学拟合关系图；
图3-3为本申请实施例中三维物体模型的前表面与目标平面模型的颜色h* 10,ab的散点分布和数学拟合关系图；
图4为本申请另一具体实施例中三维物体色度值计算装置的结构框图示意图；
图5为本申请另一具体实施例中计算机设备的结构框图示意图；
图6为本申请另一具体实施例中三维物体色度值计算系统的结构框图示意图。
具体实施方式
为了更好的理解本发明的技术方案，下面结合附图对本发明实施例进行详细描述。
应当明确，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。
在本发明实施例中使用的术语是仅仅出于描述特定实施例的目的，而非旨在限制本发明。在本发明实施例和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义。
应当理解，本文中使用的术语“和/或”仅仅是一种描述关联对象的关联关系，表示可以存在三种关系，例如，A和/或B，可以表示：单独存在A，同时存在A和B，单独存在B这三种情况。另外，本文中字符“/”，一般表示前后关联对象是一种“或”的关系。
图1为本申请具体实施例中三维物体色度值计算方法的流程示意图，本申请三维物体色度值计算方法包括以下步骤S10至步骤S60：
步骤S10，获取多个三维物体模型；
步骤S20，基于所述三维物体模型的颜色获取同一色系的多个平面模型；
步骤S30，获取置于标准观察箱中的三维物体模型与目标平面模型的光谱反射率以及标准观察箱中的照明光源的相对光谱能量分布值；其中，所述目标平面模型为所述同一色系的多个平面模型中与所述三维物体模型的颜色接近度最高的一个平面模型；
步骤S40，根据所述光谱反射率及所述光谱能量分布值分别计算得到所述三维物体模型的第一计算色度值、所述目标平面模型的第二计算色度值；
步骤S50，基于多个所述三维物体模型的第一计算色度值与多个所述目标平面模型的第二计算色度值构建得到线性拟合函数，所述线性拟合函数为拟合后的色度值计算公式；
步骤S60，根据拟合后的色度值计算公式对待评价三维物体进行色度值计算。
在本方案中，通过利用与三维物体模型颜色最接近的目标平面模型，将两者置于同一标准观察箱中，得到三维物体模型与目标平面模型的计算色度值，通过多个三维物体模型的计算色度值以及多个目标平面模型的计算色度值来构建拟合函数，从而得到拟合后的色度值计算公式，使得通过拟合后的色度值计算公式能够提高得到的三维物体的色度值的准确度。
以下结合实施例以及本申请提供的计算方法详细阐述本方案：
步骤S10，获取多个三维物体模型。
在本实施例中，三维物体模型采用三维打印技术打印成型。所述三维物体模型的最小维度大于或等于4cm，所述三维物体模型为规则形状且单色的立体模型，如球体、立方体、锥体等。示例性地，三维物体模型可以是长4cm×宽4cm×高4cm的立方体。
多个三维物体模型的颜色不同，具体地，可以打印至少5种不同色系的三维物体模型。同一色系的三维物体模型的颜色不同且同一色系的三维物体模型的数量为m，m为大于或等于4的整数。示例性地，可以打印4个红色的三维物体模型，4个蓝色的三维物体模型，4个黄色的三维物体模型，4个绿色的三维物体模型，4个灰色的三维物体模型。
所述三维物体模型的颜色选自国际照明委员会推荐的颜色中心的颜色进行打印。在本实施例中，选取的每种颜色中心的颜色打印得到4个立方体的三维物体模型且该4个立方体的三维物体模型的颜色不同，立方体的尺寸为长4cm×宽4cm×高4cm，选取的5种颜色中心的颜色分别为灰色、红色、黄色、绿色和蓝色。
在本实施例中，所述第一计算色度值和第二计算色度值均包括明度值、饱和度值及色调角值。具体地，灰色的三维物体模型在打印时设置的
色度值为(62.0，0.0，0.0)，红色的三维物体模型在打印时设置的
色度值为(44.0，37.0，23.0)，黄色的三维物体模型在打印时设置的
色度值为(87.0，-7.0，47.0)，绿色的三维物体模型在打印时设置的
色度值为(56.0，-32.0，0.0)，蓝色的三维物体模型在打印时设置的
色度值为(36.0，5.0，-31.0)。
进一步地，为了保证三维物体模型的颜色均匀性以及样本有效性，选用的三维物体模型的表面应具有均匀的颜色。所述方法还包括：
对多个拟选用的三维物体模型进行筛查，得到有效的三维物体模型。
在本实施例中，可以通过计算单个三维物体模型不同位置处的离均色差来判断三维物体模型表面颜色是否均匀，从而确定该三维物体模型是否为有效的三维物体模型。
具体地，所述离均色差的计算步骤，包括：计算单个三维物体模型任意一个位置的色度值与所述三维物体模型不同位置处的色度值的平均值之间的色差。
确定三维物体模型的离均色差小于或等于第一预设值。在本实施例中，采用CIELAB色差计算公式进行离均色差计算，第一预设值为1.0。当三维物体模型的离均色差小于或等于1.0时，则该三维物体模型的表面颜色满足均匀性要求，可以作为有效的三维物体模型。
步骤S20，基于所述三维物体模型的颜色获取同一色系的多个平面模型。
在本实施例中，平面模型采用三维打印技术打印成型，平面模型为形状规则、颜色均匀的平面图形，如圆形、矩形、三角形等模型。基于三维物体模型打印中选取的颜色中心的颜色，分别以灰色
色度值(50.0，0.0，0.0)、红色
色度值(44.0，37.0，23.0)，黄色
色度值(87.0，-7.0，47.0)，绿色
色度值(56.0，-32.0，0.0)，蓝色
色度值(36.0，5.0，-31.0)为输入值打印多个平面模型，其中，平面模型的厚度小于或等于1mm，且所述三维物体模型在所述平面模型上的正投影面积等于所述平面模型的面积。示例性地，三维物体模型长4cm×宽4cm×高4cm，平面模型为长4cm×宽4cm×高1mm的长方体；或者，三维物体模型长4cm×宽5cm×高4cm，与三维物体模型长4cm×宽5cm所对应的平面模型为长4cm×宽5cm。
并且，与每个所述三维物体模型对应的同一色系的平面模型的数量为n，n为大于或等于10的整数。在本实施例中，与每个所述三维物体模型对应的同一色系的平面模型的数量为16，即1个灰色的三维物体模型配置16个灰色的平面模型。
需要说明的是，步骤S10以及步骤S20可以同步进行，也可以异步进行，在此不做限定。打印三维物体模型的设备和打印平面模型的设备相同，所述三维物体模型在所述平面模型上的正投影面积等于所述平面模型的面积。
在步骤S30之前，所述方法还包括：
对多个平面模型进行筛查，得到有效的平面模型；
基于每个所述三维物体模型的颜色从所述同一色系的多个有效的平面模型中筛选得到目标平面模型，其中，所述目标平面模型与同一色系的所述三维物体模型的颜色接近度最高。
在本实施例中，所述对多个平面模型进行筛查，得到有效的平面模型的步骤包括：
测量得到每个所述三维物体模型的第一测量色度值；
测量得到每个所述平面模型的第二测量色度值；
基于所述第一测量色度值和第二测量色度值获取与所述三维物体模型的颜色属于同一色系的多个有效的平面模型。
为了提高实验数据的准确性，三维物体模型的第一测量色度值的测量条件、测量方法与平面模型的第二测量色度值的测量条件、测量方法相同。
具体地，使用X-Rite eXact分光光度计在D65、10°视场的测量条件下，测量得到每个所述三维物体模型的第一测量色度值，具体测量方法为：
测量每个所述三维物体模型至少5个不同位置处的色度值，并计算所述色度值的算术平均值，得到所述三维物体模型的第一测量色度值。
其中，所述至少5个不同位置位于所述三维物体模型同一平面或同一弧面上。在本实施例中，打印得到并进行有效性筛选后的三维物体模型的第一测量色度值如表1所示：
表1.三维物体模型的第一测量色度值分布
其中，L* 10表示三维物体模型的明度值；a* 10，b* 10表示三维物体模型的颜色的色品参数。
具体地，使用X-Rite eXact分光光度计在D65、10°视场的测量条件下，所述测量得到每个所述平面模型的第二测量色度值，包括：
测量每个所述平面模型至少5个不同位置处的色度值，并计算所述色度值的算术平均值，得到所述平面模型的第二测量色度值。
以4个灰色三维物体模型的上表面为例，对应的平面模型的L*10，a*10，b*10如表2所示:
表2与灰色三维物体模型的上表面对应的平面模型的第二测量色度值
基于所述第一测量色度值和第二测量色度值对平面模型的有效性筛查，可以得到与三维物体模型颜色接近的有效的平面模型。
在一个可选的实施方式中，在确定所述平面模型的第二测量色度值之前，通过计算单个平面模型不同位置处的离均色差来判断平面模型表面的颜色是否均匀，具体的平面模型离均色差的计算方法与三维物体模型离均色差的计算方法相同；确定平面模型的离均色差小于或等于第一预设值。在本实施例中，采用CIELAB色差计算公式进行离均色差计算，第一预设值为1.0。当平面模型的离均色差小于或等于1.0时，则该平面模型的表面颜色满足均匀性要求，可以作为满足有效的平面模型的前提条件。
在平面模型的表面颜色满足均匀性的前提下，将计算的平面模型的色度值的算术平均值，确定为所述平面模型的第二测量色度值。
在其他实施方式中，还可以组织奇数名观察者进行观察，人数要求是越多越好，观察者人数最少5人，要求所有观察者中要有50％以上的观察者认为平面模型与对应的三维物体模型的颜色视觉上接近，则初步认定该平面模型与该三维物体模型为同一色系且颜色接近的平面模型。示例性地，如果7名观察者中，至少有4名观察者认为平面模型的颜色与三维物体模型的颜色一致，此时认为该平面模型的颜色与三维物体模型的颜色视觉上接近，为有效的平面模型。观察者可通过标准观察箱进行观察，具体的标准观察箱的观察条件如下文所述，在此不再赘述。
在另一可选的实施方式中，还可以在组织奇数名观察者对三维物体模型的颜色和平面模型的颜色进行观察并同时参考第一测量色度值和第二测量色度值以得到有效的平面模型。
进一步地，基于每个所述三维物体模型的颜色从所述同一色系的多个有效的平面模型中筛选得到目标平面模型。
具体地，可以获取至少10名色觉正常的男性和女性观察者从所述多个平面模型中筛选得到目标平面模型。其中女性的人数多于男性，且观察者的年龄分布在20-32岁。在本实施例中，组织了31名，其中男9名，女22名，年龄分布在20～26岁的观察者，逐一比对在标准观察箱中的多个平面模型以及对应的三维物体模型在颜色上的差异，筛选得到目标平面模型。要求每个三维物体模型和对应的平面模型在观察箱中观察次数不少于30次。当目标平面模型存在多个时，选取被挑选概率最高的目标平面模型作为目标平面模型。
观察者距离标准观察箱的观察距离为50cm-60cm，实验的观察环境为暗室，实验在LEDView标准观察箱中进行，实验光源的色温为6496K，照度800lux，显色指数93.5。
每组待评价的平面模型和三维物体模型随机呈现给观察者，每个观察者在观察实验中单次持续时间为20min左右；本实施例中所述的色觉正常是基于李春慧，李育宏编著，辽宁科学技术出版社出版的《新编色觉检查图》(1994年第2版)进行观察者色觉检测，检测结果满足要求，即判断为观察者色觉正常。观察者还要求具有颜色科学的专业背景。
步骤S30，获取置于标准观察箱中的三维物体模型与目标平面模型的光谱反射率以及标准观察箱中的照明光源的相对光谱能量分布值；其中，所述目标平面模型为所述同一色系的多个平面模型中与所述三维物体模型的颜色接近度最高的一个平面模型。
在具体实施例中，可以将所述三维物体模型与所述目标平面模型置于同一标准观察箱内，测量得到所述三维物体模型的第一光谱反射率ρ1(λ)及所述目标平面模型的第二光谱反射率ρ2(λ)，其中，所述标准观察箱的照明光源的相对光谱能量分布值为S(λ)；
步骤S40，根据所述光谱反射率及所述光谱能量分布值分别计算得到所述三维物体模型的第一计算色度值、所述目标平面模型的第二计算色度值。
在本实施例中，使用光谱辐射度PR655测量标准观察箱中照明光源的相对光谱能量分布S(λ)、三维物体模型的第一光谱反射率ρ1(λ)和目标平面模型的第二光谱反射率ρ2(λ)，代入色度值计算公式，式(1)可得到三维物体模型和平面模型的三刺激值XYZ，进而可由式(2)计算得到明度值(L*)、色品参数(a*，b*)：
L*＝116(Y/Yn)1/3-16
a*＝500[(X/Xn)1/3-(Y/Yn)1/3]
b*＝200[(Y/Yn)1/3-(Z/Zn)1/3])； 式(2)；
其中，式(2)中Xn、Yn、Zn分别为标准观察箱中照明光源的三刺激值，X、Y、Z分别为三维物体模型或平面模型的三刺激值；
通过式(1)计算照明光源的三刺激值Xn、Yn、Zn时，式(1)中ρ(λ)＝1。
进一步根据式(3)及式(4)计算得到色彩饱和度值C* ab、色调角值h* ab。
对于20个三维物体模型，标准观察箱中光源照明下的三维物体模型上表面及前表面的第一计算色度值如表3所示。标准观察箱中光源照明下的目标平面模型的第二计算色度值如表4所示：
表3标准观察箱中光源照明下三维物体模型上表面及前表面的第一计算色度值
表4与三维物体模型上表面和前表面颜色视觉上最接近的平面模型的第二计算色度值
步骤S50，基于多个所述三维物体模型的第一计算色度值与多个所述目标平面模型的第二计算色度值构建得到线性拟合函数，所述线性拟合函数为拟合后的色度值计算公式。
在本实施例中，分别以上述多个三维物体模型的上表面和前表面的第一计算色度值以及对应的目标平面模型的第二计算色度值为基础，绘制散点图，并采用线性拟合的方法，拟合得到拟合后的色度值计算公式。所述拟合后的上表面色度值计算公式包括明度值计算公式、饱和度值计算公式及色调角值计算公式。
具体地，步骤S50，包括：
根据所述三维物体模型上表面和/或前表面的明度值以及目标平面模型的明度值拟合得到拟合后的上表面和/或前表面明度值计算公式；
根据所述三维物体模型上表面和/或前表面的饱和度值以及目标平面模型的饱和度值拟合得到拟合后的上表面和/或前表面饱和度值计算公式；
根据所述三维物体模型上表面和/或前表面的色调角值以及目标平面模型的色调角值拟合得到拟合后的上表面和/或前表面色调角值计算公式。
图2-1为本申请具体实施例中三维物体模型的上表面与目标平面模型的颜色L* 10的散点分布和数学拟合关系图，如图2-1所示，三维物体模型的上表面拟合后的明度值的计算公式为：
yL上＝1.0797xL上-5.6981 (A-1)，其中，xL上表示三维物体模型上表面的初始明度值，yL上表示拟合后的三维物体模型上表面的明度值。
图2-2为本申请具体实施例中三维物体模型的上表面与目标平面模型的颜色C* 10,ab的散点分布和数学拟合关系图，如图2-2所示，三维物体模型的上表面拟合后的饱和度值的计算公式为：
yC上＝1.0077xC上+0.666 (B-1)，其中，xC上表示三维物体模型上表面的初始饱和度值，yC上表示拟合后的三维物体模型上表面的饱和度值。
图2-3为本申请具体实施例中三维物体模型的上表面与目标平面模型的颜色h* 10,ab的散点分布和数学拟合关系图，如图2-3所示，三维物体模型的上表面拟合后的色调角值的计算公式为：
yh上＝0.9853xh上+1.587 (C-1)，其中，xh上表示三维物体模型上表面的初始色调角值，yh上表示拟合后的三维物体模型的上表面的色调角值。
图3-1为本申请实施例中三维物体模型的前表面与目标平面模型的颜色L* 10的散点分布和数学拟合关系图，如图3-1所示，三维物体模型的前表面拟合后的明度值的计算公式为：
yL前＝1.0304xL前-10.3 (A-2)，其中，xL前表示三维物体模型前表面的初始明度值，yL前表示拟合后的三维物体模型前表面的明度值。
图3-2为本申请实施例中三维物体模型的前表面与目标平面模型的颜色C* 10,ab的散点分布和数学拟合关系图，如图3-2所示，三维物体模型的前表面拟合后的饱和度值的计算公式为：
yC前＝1.0082xC前-1.7921 (B-2)，其中，xC前表示三维物体模型前表面的初始饱和度值，yC前表示拟合后的三维物体模型前表面的饱和度值。
图3-3为本申请实施例中三维物体模型的前表面与目标平面模型的颜色h* 10,ab的散点分布和数学拟合关系图，如图3-3所示，三维物体模型的前表面拟合后的色调角值的计算公式为：
yh前＝0.9582xh前+3.8288 (C-2)，其中，xh前表示三维物体模型上表面的初始色调角值，yh上表示拟合后的三维物体模型的前表面的色调角值。
步骤S60，根据拟合后的色度值计算公式对待评价三维物体进行色度值计算，具体包括：
获取待评价三维物体，测量得到所述待评价三维物体上表面和/或前表面的测量色度值，所述上表面和/或前表面的测量色度值包括明度值、饱和度值及色调角值；
将所述上表面和/或前表面的明度值代入拟合后的上表面和/或前表面的明度值计算公式，得到拟合后的明度值；
将所述上表面和/或前表面的饱和度值代入拟合后的上表面和/或前表面的饱和度值计算公式，得到拟合后的饱和度值；
将所述上表面和/或前表面的色调角值代入拟合后的上表面和/或前表面的色调角值计算公式，得到拟合后的色调角值；
根据拟合后的明度值、饱和度值及色调角值得到所述待评价三维物体拟合后的色度值。
在具体实施例中，可以按照上述测量三维物体模型的色度值的测量条件，测量待评价三维物体的测量色度值，待评价三维物体的测量色度值包括明度值、饱和度值及色调角值，并分别代入对应的计算公式中，从而得到待评价三维物体拟合后的色度值。本申请中待评价三维物体的测量色度值包括明度值、饱和度值及色调角值可以是直接测量或间接测量，可以是通过测量明度值、色品参数进而根据饱和度值及色调角值的计算公式得出的饱和度值和色调角值。
本申请中，使用不同形状的三维物体模型、和/或在不同的测量条件下，得到的三维物体的色度值的计算公式中系数存在差异，但这也在本申请的保护范围内。
实施例2
本申请提供一种三维物体色度值计算装置，其特征在于，如图4所示，所述装置包括：
第一获取单元，用于获取多个三维物体模型；
第二获取单元，用于基于所述三维物体模型的颜色获取同一色系的多个平面模型；
第三获取单元，用于获取置于标准观察箱中的三维物体模型与目标平面模型的光谱反射率以及标准观察箱中的照明光源的相对光谱能量分布值；
第一计算单元，用于根据所述光谱反射率及所述光谱能量分布值分别计算得到所述三维物体模型的第一计算色度值、所述目标平面模型的第二计算色度值；其中，所述目标平面模型为所述同一色系的多个平面模型中与所述三维物体模型的颜色接近度最高的一个平面模型；
构建单元，用于基于多个所述三维物体模型的第一计算色度值与多个所述目标平面模型的第二计算色度值构建得到线性拟合函数，所述线性拟合函数为拟合后的色度值计算公式；
第二计算单元，用于根据拟合后的色度值计算公式对待评价三维物体进行色度值计算。
在本实施例中，所述第一计算色度值和第二计算色度值包括明度值、饱和度值及色调角值。
具体地，所述三维物体模型为立方体，所述立方体包括上表面，所述拟合后的上表面色度值计算公式包括明度值计算公式、饱和度值计算公式及色调角值计算公式；其中，x为待评价三维物体的初始色度值，y为待评价三维物体拟合后的色度值；
所述明度值计算公式如式(A-1)为：yL上＝1.0797xL上-5.6981 (A-1)；
所述饱和度值计算公式如式(B-1)为：yc上＝1.0077xC上+0.666 (B-1)；
所述色调角值计算公式如式(C-1)为：yh上＝0.9853xh上+1.587 (C-1)。
进一步地，所述三维物体模型为立方体，所述立方体包括前表面，所述拟合后的前表面色度值计算公式包括明度值计算公式、饱和度值计算公式及色调角值计算公式；其中，x为待评价三维物体的初始色度值，y为待评价三维物体拟合后的色度值；
所述明度值计算公式如式(A-2)为：yL前＝1.0304xL前-10.3 (A-2)；
所述饱和度值计算公式如式(B-2)为：yc前＝1.0082xC前-1.7921 (B-2)；
所述色调角值计算公式如式(C-2)为：yh前＝0.9582xh前+3.8288 (C-2)。
可选地，所述第二计算单元，包括：
获取子单元，用于获取待评价三维物体，测量得到所述待评价三维物体上表面和/或前表面的测量色度值，所述上表面和/或前表面的测量色度值包括明度值、饱和度值及色调角值；
第一处理子单元，用于将所述上表面和/或前表面的明度值代入拟合后的上表面和/或前表面的明度值计算公式，得到拟合后的明度值；
第二处理子单元，用于将所述上表面和/或前表面的饱和度值代入拟合后的上表面和/或前表面的饱和度值计算公式，得到拟合后的饱和度值；
第三处理子单元，用于将所述上表面和/或前表面的色调角值代入拟合后的上表面和/或前表面的色调角值计算公式，得到拟合后的色调角值；及
输出子单元，用于根据拟合后的明度值、饱和度值及色调角值得到所述待评价三维物体拟合后的色度值。
需要说明的是本申请中三维物体色度值计算装置可以是一个集成的装置，也可以是多个分开的装置，多个分开的装置执行不同的功能。
实施例3
本申请实施例还提供一种计算机非易失性存储介质，所述存储介质包括存储的程序，在所述程序运行时控制所述存储介质所在设备执行实施例1所述的三维物体色度值计算方法中拟合后的色度值计算公式。在本实施例中，所述拟合后的色度值计算公式如实施例1中所示，在此不再赘述。
实施例4
本申请实施例还提供一种计算机设备，如图5所示，计算机设备200包括存储器202、处理器201以及存储在所述存储器202中并可在所述处理器上运行的计算机程序203，所述处理器执行所述计算机程序203时实现实施例1中的三维物体色度值计算方法中拟合后的色度值计算公式。在本实施例中，所述拟合后的色度值计算公式如实施例1中所示，在此不再赘述。
实施例5
本申请实施例还提供一种三维物体色度值计算系统，如图6所示，三维物体色度值计算系统300包括检测部件303、存储器302、处理器301以及存储在所述存储器302中并可在所述处理器上运行的计算机程序，所述检测部件303用于测量待评价三维物体的色度值(如L*2，C* 2，h* 2)，例如可以是X-Rite eXact分光光度计，所述处理器301执行所述计算机程序时实现实施例1所述的三维物体色度值计算方法中拟合后的色度值计算公式。即所述计算机程序包括拟合后的色度值计算公式，从而能够计算得到拟合后的色度值。
在本实施例中，所述拟合后的色度值计算公式如实施例1中所示，在此不再赘述。
需要说明的是，计算机设备可以是桌上型计算机、笔记本、掌上电脑及云端服务器等计算设备。计算机设备可包括，但不仅限于，处理器、存储器。本领域技术人员可以理解，图5仅仅是计算机设备的示例，并不构成对计算机设备的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件，例如计算机设备还可以包括输入输出设备、网络接入设备、总线等。
所称处理器可以是中央处理单元(Central Processing Unit，CPU)，还可以是其他通用处理器、数字信号处理器(Digital Signal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现场可编程门阵列(Field-Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
存储器可以是计算机设备的内部存储单元，例如计算机设备的硬盘或内存。存储器也可以是计算机设备的外部存储设备，例如计算机设备上配备的插接式硬盘，智能存储卡(Smart Media Card,SMC)，安全数字(Secure Digital,SD)卡，闪存卡(Flash Card)等。进一步地，存储器还可以既包括计算机设备的内部存储单元也包括外部存储设备。存储器用于存储计算机程序以及计算机设备所需的其他程序和数据。存储器还可以用于暂时地存储已经输出或者将要输出的数据。
以上仅为本申请的较佳实施例而已，并不用以限制本申请，凡在本申请的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本申请保护的范围之内。
Claims (18)
1.一种三维物体色度值计算方法，其特征在于，所述方法包括：
获取多个三维物体模型，所述三维物体模型的最小维度大于或等于4cm，所述三维物体模型为规则形状且单色的立体模型；
基于所述三维物体模型的颜色获取同一色系的多个平面模型，所述平面模型为形状规则且颜色均匀的平面模型；
获取置于标准观察箱中的三维物体模型与目标平面模型的光谱反射率以及标准观察箱中的照明光源的相对光谱能量分布值；其中，所述目标平面模型为所述同一色系的多个平面模型中与所述三维物体模型的颜色接近度最高的一个平面模型；
根据所述光谱反射率及所述光谱能量分布值分别计算得到所述三维物体模型的三刺激值以及所述目标平面模型的三刺激值，并将所述三维物体模型的三刺激值代入色度值计算公式计算得到所述三维物体模型的第一计算色度值，将所述平面模型的三刺激值代入色度值计算公式计算得到所述目标平面模型的第二计算色度值；其中，所述第一计算色度值和第二计算色度值均包括明度值、饱和度值及色调角值；
基于多个所述三维物体模型的第一计算色度值与多个所述目标平面模型的第二计算色度值构建得到线性拟合函数，所述线性拟合函数为拟合后的色度值计算公式，所述拟合后的色度值计算公式包括明度值计算公式、饱和度值计算公式及色调角值计算公式；
根据拟合后的色度值计算公式对待评价三维物体进行色度值计算。
2.根据权利要求1所述的方法，其特征在于，所述方法满足下列特征(1)至(4)中的至少一种：
(1)所述三维物体模型、所述平面模型均采用三维打印技术得到；
(2)同一色系的所述三维物体模型的数量为m，m为大于或等于4的整数；
(3)所述多个三维物体模型的颜色选自国际照明委员会推荐的颜色中心的颜色进行打印，至少包括5种不同色系，且同一色系的多个所述三维物体模型的色度值不同；
(4)与每个所述三维物体模型对应的同一色系的平面模型的数量为n，n为大于或等于10的整数。
3.根据权利要求2所述的方法，其特征在于，所述多个三维物体模型的颜色分别选自灰色、红色、黄色、绿色和蓝色进行打印，其中，灰色的
色度值为(62.0，0.0，0.0)，红色的
色度值为(44.0，37.0，23.0)，黄色的
色度值为(87.0，-7.0，47.0)，绿色的
色度值为(56.0，-32.0，0.0)，蓝色的
色度值为(36.0，5.0，-31.0)。
4.根据权利要求2或3所述的方法，其特征在于，所述平面模型的厚度小于或等于1mm，且所述三维物体模型在所述平面模型上的正投影面积等于所述平面模型的面积。
5.根据权利要求1所述的方法，其特征在于，在所述获取置于标准观察箱中的三维物体模型与目标平面模型的光谱反射率以及标准观察箱中的照明光源的相对光谱能量分布值之前，所述方法还包括：
对多个平面模型进行筛选，得到筛选后的平面模型；
基于每个所述三维物体模型的颜色从所述同一色系的多个筛选后的平面模型中筛选得到目标平面模型，其中，所述目标平面模型与同一色系的所述三维物体模型的颜色接近度最高。
6.根据权利要求5所述的方法，其特征在于，所述对多个平面模型进行筛选，得到筛选后的平面模型具体包括：
测量得到每个所述三维物体模型的第一测量色度值；
测量得到每个所述平面模型的第二测量色度值；
基于所述第一测量色度值和第二测量色度值获取与所述三维物体模型的颜色属于同一色系的多个筛选后的平面模型，其中，所述三维物体模型的第一测量色度值的测量条件、测量方法与所述平面模型的第二测量色度值的测量条件、测量方法相同。
7.根据权利要求6所述的方法，其特征在于，所述测量得到每个所述三维物体模型的第一测量色度值，包括：
测量每个所述三维物体模型至少5个不同位置处的色度值，并计算所述色度值的算术平均值，得到所述三维物体模型的第一测量色度值，其中，所述至少5个不同位置位于所述三维物体模型同一平面或同一弧面上。
8.根据权利要求6所述的方法，其特征在于，所述测量得到每个所述平面模型的第二测量色度值，包括：
测量每个所述平面模型至少5个不同位置处的色度值，并计算所述色度值的算术平均值，得到所述平面模型的第二测量色度值。
9.根据权利要求1所述的方法，其特征在于，所述三维物体模型为立方体，所述立方体包括上表面，所述拟合后的上表面色度值计算公式包括明度值计算公式、饱和度值计算公式及色调角值计算公式；其中，x为待评价三维物体的初始色度值，y为待评价三维物体拟合后的色度值；
所述明度值计算公式如式(A-1)为：yL上＝1.0797xL上-5.6981 (A-1)；
所述饱和度值计算公式如式(B-1)为：yc上＝1.0077xC上+0.666 (B-1)；
所述色调角值计算公式如式(C-1)为：yh上＝0.9853xh上+1.587 (C-1)。
10.根据权利要求1所述的方法，其特征在于，所述三维物体模型为立方体，所述立方体包括前表面，所述拟合后的前表面色度值计算公式包括明度值计算公式、饱和度值计算公式及色调角值计算公式；其中，x为待评价三维物体的初始色度值，y为待评价三维物体拟合后的色度值；
所述明度值计算公式如式(A-2)为：yL前＝1.0304xL前-10.3 (A-2)；
所述饱和度值计算公式如式(B-2)为：yc前＝1.0082xC前-1.7921 (B-2)；
所述色调角值计算公式如式(C-2)为：yh前＝0.9582xh前+3.8288 (C-2)。
11.根据权利要求9或10所述的方法，其特征在于，根据拟合后的色度值计算公式对待评价三维物体进行色度值计算，包括：
获取待评价三维物体，测量得到所述待评价三维物体上表面和/或前表面的测量色度值，所述上表面和/或前表面的测量色度值包括明度值、饱和度值及色调角值；
将所述上表面和/或前表面的明度值代入拟合后的上表面和/或前表面的明度值计算公式，得到拟合后的明度值；
将所述上表面和/或前表面的饱和度值代入拟合后的上表面和/或前表面的饱和度值计算公式，得到拟合后的饱和度值；
将所述上表面和/或前表面的色调角值代入拟合后的上表面和/或前表面的色调角值计算公式，得到拟合后的色调角值；
根据拟合后的明度值、饱和度值及色调角值得到所述待评价三维物体拟合后的色度值。
12.一种三维物体色度值计算装置，其特征在于，所述装置包括：
第一获取单元，用于获取多个三维物体模型，所述三维物体模型的最小维度大于或等于4cm，所述三维物体模型为规则形状且单色的立体模型；
第二获取单元，用于基于所述三维物体模型的颜色获取同一色系的多个平面模型，所述平面模型为形状规则且颜色均匀的平面模型；
第三获取单元，用于获取置于标准观察箱中的三维物体模型与目标平面模型的光谱反射率以及标准观察箱中的照明光源的相对光谱能量分布值；其中，所述目标平面模型为所述同一色系的多个平面模型中与所述三维物体模型的颜色接近度最高的一个平面模型；
第一计算单元，用于根据所述光谱反射率及所述光谱能量分布值分别计算得到所述三维物体模型的三刺激值以及所述目标平面模型的三刺激值，并将所述三维物体模型的三刺激值代入色度值计算公式计算得到所述三维物体模型的第一计算色度值，将所述目标平面模型的三刺激值代入色度值计算公式计算得到所述目标平面模型的第二计算色度值；其中，所述第一计算色度值和第二计算色度值均包括明度值、饱和度值及色调角值；
构建单元，用于基于多个所述三维物体模型的第一计算色度值与多个所述目标平面模型的第二计算色度值构建得到线性拟合函数，所述线性拟合函数为拟合后的色度值计算公式，所述拟合后的色度值计算公式包括明度值计算公式、饱和度值计算公式及色调角值计算公式；
第二计算单元，用于根据拟合后的色度值计算公式对待评价三维物体进行色度值计算。
13.根据权利要求12所述的装置，其特征在于，所述三维物体模型为立方体，所述立方体包括上表面，所述拟合后的上表面色度值计算公式包括明度值计算公式、饱和度值计算公式及色调角值计算公式；其中，x为待评价三维物体的初始色度值，y为待评价三维物体拟合后的色度值；
所述明度值计算公式如式(A-1)为：yL上＝1.0797xL上-5.6981 (A-1)；
所述饱和度值计算公式如式(B-1)为：yc上＝1.0077xC上+0.666 (B-1)；
所述色调角值计算公式如式(C-1)为：yh上＝0.9853xh上+1.587 (C-1)。
14.根据权利要求12所述的装置，其特征在于，所述三维物体模型为立方体，所述立方体包括前表面，所述拟合后的前表面色度值计算公式包括明度值计算公式、饱和度值计算公式及色调角值计算公式；其中，x为待评价三维物体的初始色度值，y为待评价三维物体拟合后的色度值；
所述明度值计算公式如式(A-2)为：yL前＝1.0304xL前-10.3 (A-2)；
所述饱和度值计算公式如式(B-2)为：yc前＝1.0082xC前-1.7921 (B-2)；
所述色调角值计算公式如式(C-2)为：yh前＝0.9582xh前+3.8288 (C-2)。
15.根据权利要求13或14所述的装置，其特征在于，所述第二计算单元，包括：
获取子单元，用于获取待评价三维物体，测量得到所述待评价三维物体上表面和/或前表面的测量色度值，所述上表面和/或前表面的测量色度值包括明度值、饱和度值及色调角值；
第一处理子单元，用于将所述上表面和/或前表面的明度值代入拟合后的上表面和/或前表面的明度值计算公式，得到拟合后的明度值；
第二处理子单元，用于将所述上表面和/或前表面的饱和度值代入拟合后的上表面和/或前表面的饱和度值计算公式，得到拟合后的饱和度值；
第三处理子单元，用于将所述上表面和/或前表面的色调角值代入拟合后的上表面和/或前表面的色调角值计算公式，得到拟合后的色调角值；及
输出子单元，用于根据拟合后的明度值、饱和度值及色调角值得到所述待评价三维物体拟合后的色度值。
16.一种计算机非易失性存储介质，其特征在于，所述存储介质包括存储的程序，在所述程序运行时控制所述存储介质所在设备执行权利要求1至11任意一项所述的三维物体色度值计算方法中拟合后的色度值计算公式。
17.一种计算机设备，其特征在于，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至11任意一项所述的三维物体色度值计算方法中拟合后的色度值计算公式。
18.一种三维物体色度值计算系统，包括检测部件、存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至11任意一项所述的三维物体色度值计算方法中拟合后的色度值计算公式。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110242021.4A CN113049104B (zh) | 2021-03-04 | 2021-03-04 | 三维物体色度值计算方法及装置、三维物体色度值计算系统 |
| PCT/CN2021/111939 WO2022183683A1 (zh) | 2021-03-04 | 2021-08-11 | 三维物体色度值计算方法及装置、三维物体色度值计算系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110242021.4A CN113049104B (zh) | 2021-03-04 | 2021-03-04 | 三维物体色度值计算方法及装置、三维物体色度值计算系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113049104A CN113049104A (zh) | 2021-06-29 |
| CN113049104B true CN113049104B (zh) | 2022-07-12 |
Family
ID=76510009
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110242021.4A Active CN113049104B (zh) | 2021-03-04 | 2021-03-04 | 三维物体色度值计算方法及装置、三维物体色度值计算系统 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN113049104B (zh) |
| WO (1) | WO2022183683A1 (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113049104B (zh) * | 2021-03-04 | 2022-07-12 | 珠海赛纳三维科技有限公司 | 三维物体色度值计算方法及装置、三维物体色度值计算系统 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5495429A (en) * | 1993-02-12 | 1996-02-27 | West Virginia University | Method and apparatus for measuring the color of three dimensional objects |
| CN105674912A (zh) * | 2016-01-26 | 2016-06-15 | 中国科学院上海光学精密机械研究所 | 结合达曼光栅的多光刀彩色三维测量装置和方法 |
| CN105825020A (zh) * | 2016-03-23 | 2016-08-03 | 天津师范大学 | 三维可感知色域计算方法 |
| JP2019032187A (ja) * | 2017-08-04 | 2019-02-28 | ４Ｄセンサー株式会社 | カラー物体の３次元形状とカラー情報とを同時に取得可能な計測方法、計測装置、計測プログラムを記録した、コンピュータ読み取り可能な記録媒体 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5186457B2 (ja) * | 2009-09-14 | 2013-04-17 | 富士フイルム株式会社 | 測色値算出方法、プロファイル生成方法並びに色変換方法、その装置及びプログラム |
| CN113049104B (zh) * | 2021-03-04 | 2022-07-12 | 珠海赛纳三维科技有限公司 | 三维物体色度值计算方法及装置、三维物体色度值计算系统 |
-
2021
- 2021-03-04 CN CN202110242021.4A patent/CN113049104B/zh active Active
- 2021-08-11 WO PCT/CN2021/111939 patent/WO2022183683A1/zh active Application Filing
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5495429A (en) * | 1993-02-12 | 1996-02-27 | West Virginia University | Method and apparatus for measuring the color of three dimensional objects |
| CN105674912A (zh) * | 2016-01-26 | 2016-06-15 | 中国科学院上海光学精密机械研究所 | 结合达曼光栅的多光刀彩色三维测量装置和方法 |
| CN105825020A (zh) * | 2016-03-23 | 2016-08-03 | 天津师范大学 | 三维可感知色域计算方法 |
| JP2019032187A (ja) * | 2017-08-04 | 2019-02-28 | ４Ｄセンサー株式会社 | カラー物体の３次元形状とカラー情報とを同時に取得可能な計測方法、計測装置、計測プログラムを記録した、コンピュータ読み取り可能な記録媒体 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113049104A (zh) | 2021-06-29 |
| WO2022183683A1 (zh) | 2022-09-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Mokrzycki et al. | Colour difference∆ EA survey | |
| JP6039109B2 (ja) | 着色検査装置および着色検査方法 | |
| JP6039008B2 (ja) | 着色評価装置及び着色評価方法 | |
| JP2006215756A (ja) | 画像処理装置および画像処理方法ならびにそのプログラム | |
| JP3809838B2 (ja) | 照明源固有ｈｓｖ色座標を用いた画像ハイライト補正方法、画像ハイライト補正プログラム、および画像取得システム | |
| Menk et al. | Visualisation techniques for using spatial augmented reality in the design process of a car | |
| Nam et al. | Multispectral photometric stereo for acquiring high-fidelity surface normals | |
| CN113049104B (zh) | 三维物体色度值计算方法及装置、三维物体色度值计算系统 | |
| JPH08297054A (ja) | 色感測定装置 | |
| US20230260236A1 (en) | Displaying a virtual object in a real-life scene | |
| Macdonald | Realistic visualisation of cultural heritage objects | |
| MacDonald | Colour and directionality in surface reflectance | |
| Haneishi et al. | Goniospectral imaging of three-dimensional objects | |
| JP2014187558A (ja) | 画像色分布検査装置および画像色分布検査方法 | |
| JP2022177166A (ja) | 検査方法、プログラム、及び、検査システム | |
| CN116057576A (zh) | 在异质测量环境中将至少两种材料的外观可视化 | |
| JP7235038B2 (ja) | 光沢値算出装置、光沢値測定装置、光沢色の色調定量化装置および光沢値算出方法 | |
| WO2022242608A1 (zh) | 物体喜好记忆色的获取方法及喜好记忆色标准色卡 | |
| WO2022170747A1 (zh) | 三维物体色差计算方法及装置、三维物体色差计算系统 | |
| JP2016194449A (ja) | 着色検査装置および着色検査方法 | |
| McCann et al. | Pixel and spatial mechanisms of color constancy | |
| Elkhuizen et al. | Reproducing oil paint gloss in print for the purpose of creating reproductions of Old Masters | |
| US20150124273A1 (en) | Tone definitional equation generation method, tone type identification method, a method for calculating munsell value corresponding to arbitrarily defined tone, image formation method, and interface screen display device | |
| Wang et al. | Retinex-based color correction for displaying high dynamic range images | |
| Galuza et al. | Modeling of an object chromaticity with a given emission spectrum |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |