发明名称 一种利用四面体插值法实现测色色差计颜色校正的方法
摘要 一种利用四面体插值法实现测色色差计颜色校正的方法,首先通过传统的插值算法测得所需数量的基色点,然后,利用三维Delaunay三角剖分算法,对基色集进行四面体剖分,利用剖分得到的四面体集合,建立查找表;生成查找表后,假定输入点为P,得到A后,可以利用点p与A中每一个四面体的体积关系,进一步判断输入点p位于A中的哪个四面体内,如果B中有多于两个的四面体,说明p点是某个四面体的顶点,也就是说属于插值点集中的点,那么无需插值,可直接输出目标RGB值;假设成像系统的映射呈线性关系,则输出的标准颜色值为pn=(Rn,Gn,Bn)T可以表示为:pn-p=Λnσ与普通基于三色分析的测色方法相比较,本发明具有简单、快捷、准确的特点。
申请公布号 CN102393965B 申请公布日期 2014.01.29
申请号 CN201110156725.6 申请日期 2011.06.01
申请人 博雅智汇国际有限公司 发明人 林金龙;欧阳泳春;王果
分类号 G06T7/40(2006.01)I;G01J3/46(2006.01)I 主分类号 G06T7/40(2006.01)I
代理机构 广州市红荔专利代理有限公司 44214 代理人 吴世民
主权项 一种利用四面体插值法实现测色色差计颜色校正的方法,其特征在于:首先通过传统的插值算法测得所需数量的基色点,然后,利用三维Delaunay三角剖分算法,对基色点集进行四面体剖分,利用剖分得到的四面体集合,建立由每个四面体的外接圆圆心、外接圆半径、四面体顶点原RGB值和目标RGB值构成的查找表;生成查找表后,假定输入点为p,工作流程如下:利用输入点p与每个四面体的外接球球心O的距离d和四面体外接球的半径r的关系,判断输入点p是否落入四面体的外接球内,p落入其外接球的所有四面体构成集合A;具体判断方法如下: <mrow> <mi>p</mi> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mo>&Element;</mo> <mi>T</mi> </mtd> <mtd> <mi>d</mi> <mo>&le;</mo> <mi>r</mi> </mtd> </mtr> <mtr> <mtd> <mo>&NotElement;</mo> <mi>T</mi> </mtd> <mtd> <mi>d</mi> <mo>></mo> <mi>r</mi> </mtd> </mtr> </mtable> </mfenced> </mrow>得到A后,利用点p与A中每一个四面体的体积关系,进一步判断输入点p位于A中的哪个四面体内,只有当点p落入顶点为p0,p1,p2,p3的四面体T内时,点p与p0,p1,p2构成的四面体T1,与p0,p1,p3构成的四面体T2,与p0,p2,p3构成的四面体T3,与p1,p2,p3构成的四面体T4,才满足以下关系:V(T1)+V(T2)+V(T3)+V(T4)=V(T)其中,四面体体积的计算通过下述得到,当四面体T的四个顶点为p0(x0,y0,z0),p1(x1,y1,z1),p2(x2,y2,z2),p3(x3,y3,z3)时,该四面体的体积V(T)的计算公式为: <mrow> <mi>V</mi> <mrow> <mo>(</mo> <mi>T</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mn>6</mn> </mfrac> <mo>|</mo> <mo>|</mo> <mfenced open='' close=''> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>x</mi> <mn>0</mn> </msub> </mtd> <mtd> <msub> <mi>x</mi> <mn>1</mn> </msub> </mtd> <mtd> <msub> <mi>x</mi> <mn>2</mn> </msub> </mtd> <mtd> <msub> <mi>x</mi> <mn>3</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>y</mi> <mn>0</mn> </msub> </mtd> <mtd> <msub> <mi>y</mi> <mn>1</mn> </msub> </mtd> <mtd> <msub> <mi>y</mi> <mn>2</mn> </msub> </mtd> <mtd> <msub> <mi>y</mi> <mn>3</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>z</mi> <mn>0</mn> </msub> </mtd> <mtd> <msub> <mi>z</mi> <mn>1</mn> </msub> </mtd> <mtd> <msub> <mi>z</mi> <mn>2</mn> </msub> </mtd> <mtd> <msub> <mi>z</mi> <mn>3</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>|</mo> <mo>|</mo> </mrow>所以判断点p是否属于某个四面体T时,用如下公式判断: <mrow> <mi>p</mi> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mo>&Element;</mo> <mi>T</mi> </mtd> <mtd> <mi>V</mi> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mi>V</mi> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mi>V</mi> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mi>V</mi> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mi>V</mi> <mrow> <mo>(</mo> <mi>T</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mo>&NotElement;</mo> <mi>T</mi> </mtd> <mtd> <mi>V</mi> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mi>V</mi> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mi>V</mi> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mi>V</mi> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> <mo>&NotEqual;</mo> <mi>V</mi> <mrow> <mo>(</mo> <mi>T</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>记录点p所在的四面体T于集合B中,如果B为空,那么p点超出能够校正的范围,不在该算法的考虑范围之内;如果B中有两个四面体,则选择任意一个,利用线性回归法插值计算得到目标值;如果B中有多于两个的四面体,说明p点是某个四面体的顶点,也就是说属于插值点集中的点,那么无需插值,直接输出目标值;通过以下方式得到线性回归方程:假设成像系统的映射呈线性关系,则输出的标准颜色值为pn=(Rn,Gn,Bn)T表示为:pn‑p0=Λnσ其中,p0为零点坐标,Λ为包含了光照因子和传感器相应因子的3×n维矩阵,σ为物体表面反射率基函数的权重系数;同样在非标准光照条件下得到的颜色值pa=(Ra,Ga,Ba)T满足:pa‑p0=Λaσ联立上面两式即可得到 <mrow> <msub> <mi>p</mi> <mi>n</mi> </msub> <mo>=</mo> <msub> <mi>&Lambda;</mi> <mi>n</mi> </msub> <msubsup> <mi>&Lambda;</mi> <mi>a</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <msub> <mi>p</mi> <mi>a</mi> </msub> <mo>+</mo> <mrow> <mo>(</mo> <mi>I</mi> <mo>-</mo> <msub> <mi>&Lambda;</mi> <mi>n</mi> </msub> <msubsup> <mi>&Lambda;</mi> <mi>a</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>)</mo> </mrow> <msub> <mi>p</mi> <mn>0</mn> </msub> </mrow>令 <mrow> <mi>&Lambda;</mi> <mo>=</mo> <msub> <mi>&Lambda;</mi> <mi>n</mi> </msub> <msubsup> <mi>&Lambda;</mi> <mi>a</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>=</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>a</mi> <mn>11</mn> </msub> </mtd> <mtd> <msub> <mi>a</mi> <mn>12</mn> </msub> </mtd> <mtd> <msub> <mi>a</mi> <mn>13</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>a</mi> <mn>21</mn> </msub> </mtd> <mtd> <msub> <mi>a</mi> <mn>22</mn> </msub> </mtd> <mtd> <msub> <mi>a</mi> <mn>23</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>a</mi> <mn>31</mn> </msub> </mtd> <mtd> <msub> <mi>a</mi> <mn>32</mn> </msub> </mtd> <mtd> <msub> <mi>a</mi> <mn>33</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> <mi>B</mi> <mo>=</mo> <mrow> <mo>(</mo> <mi>I</mi> <mo>-</mo> <msub> <mi>&Lambda;</mi> <mi>n</mi> </msub> <msubsup> <mi>&Lambda;</mi> <mi>a</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>)</mo> </mrow> <msub> <mi>p</mi> <mn>0</mn> </msub> <mo>=</mo> <mfenced open='(' close=')'> <mtable> <mtr> <mtd> <msub> <mi>b</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>b</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>b</mi> <mn>3</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> </mrow>则上式变为pn=Λpa+B。
地址 中国香港九龙美孚荔欣苑荔林阁35楼12室