[发明专利]一种基于视网膜视觉机制的彩色图像颜色恒常方法

权利要求书

1.一种基于视网膜视觉机制的彩色图像颜色恒常方法，包括以下步骤：

步骤1：设定近似感受野模板大小及确定相应的模型参数；

设定3×3像素的近似感受野模板，其中感受野中心区为近似感受野模板的中心像素点，其半径近似为0.5；感受野外周抑制区由近似感受野模板的周围8个像素点构成，即由周围8个抑制亚区构成，感受野外周抑制区的半径近似为设定高斯核函数和的参数

步骤2：对原始彩色色偏图像的每一个像素点分别提取红色R、绿色G、蓝色B和黄色Y颜色分量I_R、I_G、I_B和I_Y，并计算各分量的n阶导数n＝0、1或2；

步骤3：利用公式

AR1=AR2=AR3=(IR(n))2+(IG(n))2+(IB(n))2+(IY(n))2/IR(n)]]>

AG1=AG2=AG3=(IR(n))2+(IG(n))2+(IB(n))2+(IY(n))2/IG(n)]]>

AB1=AB2=AB3=(IR(n))2+(IG(n))2+(IB(n))2+(IY(n))2/IB(n)]]>

AY1=AY2=AY3=(IR(n))2+(IG(n))2+(IB(n))2+(IY(n))2/IY(n)]]>

计算出近似感受野模板的各个敏感度参数A_R1、A_R2、A_R3、A_G1、A_G2、A_G3、A_B1、A_B2、A_B3、A_Y1、A_Y2和A_Y3；其中：A_R1表示近似感受野模板中心区的红光敏感度，A_R2表示近似感受野模板中某一抑制亚区对中心区的红光抑制性作用的敏感度，A_R3表示近似感受野模板中另一抑制亚区对某一抑制亚区的红光抑制性作用的敏感度，A_G1表示近似感受野模板中心区的绿光敏感度，A_G2表示近似感受野模板中某一抑制亚区对中心区的绿光抑制性作用的敏感度，A_G3表示近似感受野模板中另一抑制亚区对某一抑制亚区的绿光抑制性作用的敏感度，A_B1表示近似感受野模板中心区的蓝光敏感度，A_B2表示近似感受野模板中某一抑制亚区对中心区的蓝光抑制性作用的敏感度，A_B3表示近似感受野模板中另一抑制亚区对某一抑制亚区的蓝光抑制性作用的敏感度，A_Y1表示近似感受野模板中心区的黄光敏感度，A_Y2表示近似感受野模板中某一抑制亚区对中心区的黄光抑制性作用的敏感度，A_Y3表示近似感受野模板中另一抑制亚区对某一抑制亚区的黄光抑制性作用的敏感度；

依据步骤1所确定的近似感受野模板，将原始彩色色偏图像的每一个像素点作为一个近似感受野模板的中心区，与该像素点相邻的周围8个像素点作为近似感受野外周区的8个抑制亚区，按照从左至右、从上到下的顺序依次对每一个近似感受野模板进行下述步骤4至步骤9的操作：

步骤4：根据公式

TRc=IR(x0,y0)AR1G1]]>

TGc=IG(x0,y0)AG1G1]]>

TBc=IB(x0,y0)AB1G1]]>

TYc=IY(x0,Y0)AY1G1]]>

计算出每个近似感受野模板中心区(x₀，y₀)像素点对红、绿、蓝、黄四色光兴奋反应大小T_Rc、T_Gc、T_Bc、T_Yc；其中为中心区(x₀，y₀)像素点上的红光分量，为中心区(x₀，y₀)像素点上的绿光分量，为中心区(x₀，y₀)像素点上的蓝光分量，为中心区(x₀，y₀)像素点上的黄光分量；

步骤5：根据公式

II_{R[(x，y)-(x+m，y+n)]}＝I_R(x+m，y+n)A_R3G₃

II_{G[(x，y)-(x+m，y+n)]}＝I_G(x+m，y+n)A_G3G₃

II_{B[(x，y)-(x+m，y+n)]}＝I_B(x+m，y+n)A_B3G₃

II_{Y[(x，y)-(x+m，y+n)]}＝I_Y(x+m，y+n)A_Y3G₃

计算出每个近似感受野模板外周抑制区内某一抑制亚区(x，y)受到周围另一抑制亚区(x+m，y+n)的红、绿、蓝、黄四色光的抑制性作用强度II_{R[(x，y)-(x+m，y+n)]}、II_{G[(x，y)-(x+m，y+n)]}、II_{B[(x，y)-(x+m，y+n)]}、II_{Y[(x，y)-(x+m，y+n)]}；其中I_R(x+m，y+n)为另一抑制亚区(x+m，y+n)的红光分量，I_G(x+m，y+n)为另一抑制亚区(x+m，y+n)的绿光分量，I_B(x+m，y+n)为另一抑制亚区(x+m，y+n)的蓝光分量，I_Y(x+m，y+n)为另一抑制亚区(x+m，y+n)的黄光分量；m表示另一抑制亚区(x+m，y+n)在x方向上距某一抑制亚区(x，y)的距离，n表示另一抑制亚区(x+m，y+n)在y方向上距某一抑制亚区(x，y)的距离；

步骤6：根据公式

IIR(x,y)=Σ(x+m,y+n)∈N/(x,y)IR(x+m,y+n)AR3G3]]>

IIG(x,y)=Σ(x+m,y+n)∈N/(x,y)IG(x+m,y+n)AG3G3]]>

IIB(x,y)=Σ(x+m,y+n)∈N/(x,y)IB(x,+m,y+n)AB3G3]]>

IIY(x,y)=Σ(x+m,y+n)∈N/(x,y)IY(x+m,y+n)AY3G3]]>

计算出每个近似感受野模板外周抑制区内某一抑制亚区(x，y)受到周围其它所有抑制亚区的红、绿、蓝、黄四色光的抑制性作用强度的线性总和II_R(x，y)、II_G(x，y)、II_R(x，y)、II_Y(x，y)；其中N/(x，y)表示近似感受野模板外周区中除去(x，y)抑制亚区的其他所有抑制亚区；

步骤7：根据公式

TRN=Σ(x,y)∈Nmax{AR2G2(IR(x,y)-IIR(x,y)),0}]]>

TGN=Σ(x,y)∈Nmax{AG2G2(IG(x,y)-IIG(x,y)),0}]]>

TBN=Σ(x,y)∈Nmax{AB2G2(IB(x,y)-IIB(x,y)),0}]]>

TYN=Σ(x,y)∈Nmax{AY2G2(IY(x,y)-IIY(x,y)),0}]]>

计算出每个近似感受野模板外周抑制区内各抑制亚区对中心区总的红、绿、蓝、黄四色光的抑制性作用强度T_RN、T_GN、T_BN、T_YN；其中N表示近似感受野模板外周区中所有抑制亚区，max表示二者之中取较大者；

步骤8：根据公式

T_R＝max(T_Rc-T_GN，0)

T_G＝max(T_Gc-T_RN，0)

T_B＝max(T_Bc-T_YN，0)

T_Y＝max(T_Yc-T_GN，0)

计算出每一个近似感受野模板的红-绿互拮抗输出T_R和T_G，以及蓝-黄互拮抗的输出T_B和T_Y，取T_R、T_G和T_B三个输出作为每个近似感受野模板中心区(x₀，y₀)像素点的新的红色R、绿色G和蓝色B的分量；

步骤9：将步骤8所得的所有中心区像素点的新的红色R、绿色G和蓝色B的分量按像素点一一合成，就得到一幅校正了色偏，即颜色恒常的彩色图像。