发明名称 一种基于特征的实时的视频和图像抽象化方法
摘要 本发明公开了一种基于特征的实时的视频图像抽象化方法,该方法包括将实时的视频和图像的红绿蓝颜色空间转换到CIE-Lab颜色空间、采用双边网格对亮度通道进行快速双边滤波、伪量子化亮度值、采用基于特征流的各向异性的高斯差分线条抽取算法抽取线条区域并对抽取的线条区域和近似亮度值进行组合操作使线条区域的线条增强、将实时的视频和图像由CIE-Lab颜色空间重新转换回红绿蓝颜色空间五个步骤。本发明每一个步骤都可以在家用电脑的图形卡硬件中并行处理,有效地解决了现有多实时的视频和图像抽象化方法存在的抽取线条生硬,不连贯的问题,为非专业人员提供了一种简便直观的实时的视频和图像的抽象化方法。
申请公布号 CN101286228B 申请公布日期 2011.10.05
申请号 CN200810062384.4 申请日期 2008.05.15
申请人 浙江大学 发明人 赵汉理;金小刚
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 杭州天勤知识产权代理有限公司 33224 代理人 胡红娟
主权项 1.一种基于特征的实时的视频和图像抽象化方法,该方法包括以下五个步骤:(1)实时的视频和图像读取后,按照下面的公式由红绿蓝颜色空间转换到CIE-Lab颜色空间,使实时的视频和图像的亮度通道和彩色通道相分离;<maths num="0001"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><mi>L</mi><mo>=</mo><mn>0.299</mn><mo>*</mo><mi>R</mi><mo>+</mo><mn>0.587</mn><mo>*</mo><mi>G</mi><mo>+</mo><mn>0.114</mn><mo>*</mo><mi>B</mi></mtd></mtr><mtr><mtd><mi>a</mi><mo>=</mo><mn>0.713</mn><mo>*</mo><mrow><mo>(</mo><mi>R</mi><mo>-</mo><mi>L</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>b</mi><mo>=</mo><mn>0.564</mn><mo>*</mo><mrow><mo>(</mo><mi>B</mi><mo>-</mo><mi>L</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced></math>]]></maths>式中R,G,B分别表示红绿蓝颜色值,L是亮度通道值,a和b是两个彩色通道值;(2)采用双边网格对亮度通道进行快速双边滤波,得到实时的视频和图像的近似亮度值,方法如下:创建双边网格,用三维坐标来定义一个双边网格的数据结构Γ,其中前两维为图像上的二维坐标(x,y),第三维是图像的亮度,对每个双边网格(i,j,k)先初始化为:Γ(i,j,k)=(0,0),将实时的视频和图像的所有象素置入对应的网格中,对Γ按公式Γ([x/s<sub>s</sub>],[y/s<sub>s</sub>],[I(x,y)/s<sub>r</sub>])+=(I(x,y),1)进行累加,式中s<sub>s</sub>和s<sub>r</sub>分别表示空间域和亮度域的采样率,[·]为最接近的整数值操作符,(I(x,y),1)是齐次值的表示方法,第二项是齐次项,双边网格创建完成后,对双边网格的齐次值相对于三维方向分别进行高斯卷积操作,并将所得结果除去齐次项,再将该三维网格切片,重构出二维图像结果,然后利用图形硬件的双线性插值功能取得双边网格中位于坐标(x/s<sub>s</sub>,y/s<sub>s</sub>,I(x,y)/s<sub>r</sub>)的值,即实时的视频和图像的近似亮度值;(3)伪量子化实时的视频和图像的近似亮度值,方法如下:按照公式<img file="FSB00000549074100012.GIF" wi="1132" he="84" />计算每一层的亮度值,式中Q为所计算的伪量子化亮度值,Δq是离散化宽度,q<sub>nearest</sub>是最接近I(x,y)的离散值,<img file="FSB00000549074100013.GIF" wi="44" he="45" />用来控制相邻两个离散值之间变化的快慢程度;(4)采用基于特征流的各向异性的高斯差分线条抽取算法抽取线条区域并对抽取的线条区域和近似亮度值进行组合操作使线条区域的线条增强;其中,所述的线条增强的方法如下:创建一个特征流场V(x,y),初始值垂直于输入图像的梯度场<img file="FSB00000549074100021.GIF" wi="384" he="60" />对该特征流场进行滤波操作并计算出流场方向的高斯卷积值和垂直于流场方向的卷积值,再进行差分计算,最后利用各向异性的高斯差分算法抽取线条区域并对抽取的线条区域和近似亮度值进行组合操作使线条区域的线条增强;所述的滤波操作中采用如下特征流滤波器:<maths num="0002"><![CDATA[<math><mrow><msup><mi>V</mi><msup><mi>new</mi><mo>&prime;</mo></msup></msup><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><msub><mi>k</mi><mi>h</mi></msub></mfrac><munder><mi>&Sigma;</mi><mrow><mi>Y</mi><mo>&Element;</mo><msub><mi>&Omega;</mi><mi>h</mi></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow></mrow></munder><msub><mi>w</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>X</mi><mo>,</mo><mi>Y</mi><mo>)</mo></mrow><msub><mi>w</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>X</mi><mo>,</mo><mi>Y</mi><mo>)</mo></mrow><msup><mi>V</mi><mi>cur</mi></msup><mrow><mo>(</mo><mi>Y</mi><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msup><mi>V</mi><mi>new</mi></msup><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><msub><mi>k</mi><mi>v</mi></msub></mfrac><munder><mi>&Sigma;</mi><mrow><mi>Y</mi><mo>&Element;</mo><msub><mi>&Omega;</mi><mi>v</mi></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow></mrow></munder><msub><mi>w</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>X</mi><mo>,</mo><mi>Y</mi><mo>)</mo></mrow><msub><mi>w</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>X</mi><mo>,</mo><mi>Y</mi><mo>)</mo></mrow><msup><mi>V</mi><msup><mi>new</mi><mo>&prime;</mo></msup></msup><mrow><mo>(</mo><mi>Y</mi><mo>)</mo></mrow></mrow></math>]]></maths>式中,Ω<sub>h</sub>(X)和Ω<sub>v</sub>(X)分别表示象素X在水平方向及竖直方向上的邻域,k<sub>h</sub>和k<sub>v</sub>为相应的归一化项;强度权因子函数w<sub>m</sub>为:<maths num="0004"><![CDATA[<math><mrow><msub><mi>w</mi><mi>m</mi></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mi>tanh</mi><mrow><mo>(</mo><mi>&eta;</mi><mrow><mo>(</mo><mover><mi>g</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>Y</mi><mo>)</mo></mrow><mo>-</mo><mover><mi>g</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>)</mo></mrow><mo>)</mo></mrow></mrow></math>]]></maths>式中,<img file="FSB00000549074100025.GIF" wi="30" he="52" />表示梯度的强度值,参数η用来控制双曲正切函数的下降率;方向权因子函数w<sub>d</sub>定义为:w<sub>d</sub>(X,Y)=dot_product(V(X),V(Y)),w<sub>d</sub>(X,Y)表示(V(X),V(Y))的点积;V(X)表示之前计算得出的X象素的特征流向量;所述的线条区域如下:<maths num="0005"><![CDATA[<math><mrow><mi>F</mi><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mrow><mi>s</mi><mo>&Element;</mo><mi>c</mi><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow></mrow></munder><mi>G</mi><mrow><mo>(</mo><mo>|</mo><mo>|</mo><mi>s</mi><mo>-</mo><mi>X</mi><mo>|</mo><mo>|</mo><mo>,</mo><msub><mi>&sigma;</mi><mi>m</mi></msub><mo>)</mo></mrow><mi>W</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><mi>W</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mrow><mi>Y</mi><mo>&Element;</mo><mi>l</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow></munder><mrow><mo>(</mo><mi>G</mi><mrow><mo>(</mo><mo>|</mo><mo>|</mo><mi>Y</mi><mo>-</mo><mi>X</mi><mo>|</mo><mo>|</mo><mo>,</mo><msub><mi>&sigma;</mi><mi>c</mi></msub><mo>)</mo></mrow><mo>-</mo><mi>&rho;</mi><mo>&CenterDot;</mo><mi>G</mi><mrow><mo>(</mo><mo>|</mo><mo>|</mo><mi>Y</mi><mo>-</mo><mi>X</mi><mo>|</mo><mo>|</mo><mo>,</mo><msub><mi>&sigma;</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mi>I</mi><mrow><mo>(</mo><mi>Y</mi><mo>)</mo></mrow></mrow></math>]]></maths>式中c(X),l(s)分别表示沿着象素X∈I的流切线方向和垂直象素s∈c(X)流方向的邻域,因此W(s)是沿着垂直于V(s)并与s相交的线l<sub>s</sub>进行卷积操作的;阈值参数ρ影响该边缘检测器的敏感度;(5)实时的视频和图像由CIE-Lab颜色空间重新转换回红绿蓝颜色空间,按照下面的公式进行:<maths num="0007"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><mi>R</mi><mo>=</mo><mi>L</mi><mo>+</mo><mn>1.403</mn><mo>*</mo><mi>a</mi></mtd></mtr><mtr><mtd><mi>G</mi><mo>=</mo><mi>L</mi><mo>-</mo><mn>0.714</mn><mo>*</mo><mi>a</mi><mo>-</mo><mn>0.334</mn><mo>*</mo><mi>b</mi></mtd></mtr><mtr><mtd><mi>B</mi><mo>=</mo><mi>L</mi><mo>+</mo><mn>1.773</mn><mo>*</mo><mi>a</mi></mtd></mtr></mtable></mfenced></math>]]></maths>式中R,G,B分别表示红绿蓝颜色值,L是亮度通道值,a和b是两个彩色通道值。
地址 310027 浙江省杭州市西湖区浙大路38号