发明名称 一种基于颜色差分的彩色图像边缘检测方法
摘要 本发明公开了一种基于颜色差分的彩色图像边缘检测方法,属于数字图像处理领域,该方法对平滑滤波操作具有更强的鲁棒性,可以排除大量的噪声边缘,只保留反映物体基本结构的真实边缘。平滑滤波去噪。构成邻域点对。对于平滑滤波后的图像I<sup>C</sup>(x,y),将每个像素(x,y)相邻的八个像素分别按照水平、竖直、45°和135°总共四个方向分成四个邻域点对。计算颜色相离结果。计算颜色差分图,得到的颜色差分图CDM(x,y)即为最终的边缘检测结果。能够有效抑制噪声干扰,排除大量人们不感兴趣的噪声边缘,仅保留反映物体基本结构变化的真实边缘;得到的边缘对于平滑滤波具有较好的鲁棒性,在滤波窗口尺寸较大时仍然能够具有较强的边缘响应,不易造成重要边缘丢失和断裂。
申请公布号 CN105957067A 申请公布日期 2016.09.21
申请号 CN201610258532.4 申请日期 2016.04.23
申请人 北京工业大学 发明人 卓力;胡笑尘;张菁;李晓光
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 北京思海天达知识产权代理有限公司 11203 代理人 沈波
主权项 一种基于颜色差分的彩色图像边缘检测方法,其特征在于:该方法具体包括以下步骤:,步骤1,平滑滤波去噪;对输入彩色图像进行平滑滤波,抑制噪声干扰;该过程由下式表示:<maths num="0001"><math><![CDATA[<mrow><msup><mi>I</mi><mi>C</mi></msup><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><munder><mrow><mi>S</mi><mi>m</mi><mi>o</mi><mi>o</mi><mi>t</mi><mi>h</mi></mrow><mrow><mi>&Omega;</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></munder><mo>&lsqb;</mo><msup><mi>H</mi><mi>C</mi></msup><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000972787850000011.GIF" wi="1342" he="94" /></maths>式中,H<sup>C</sup>(x,y)为输入彩色图像,C代表图像的颜色通道,C={R,G,B},R、G和B分别表示红绿蓝三种颜色;<img file="FDA0000972787850000012.GIF" wi="165" he="87" />代表以Ω(n)为窗口尺寸的平滑滤波,所述平滑滤波的种类包括均值滤波、高斯滤波、中值滤波;I<sup>C</sup>(x,y)为平滑滤波后的图像;步骤2,构成邻域点对;对于平滑滤波后的图像I<sup>C</sup>(x,y),将每个像素(x,y)相邻的八个像素分别按照水平、竖直、45°和135°总共四个方向分成四个邻域点对;这四个邻域点对分别为(x‑1,y)与(x+1,y)、(x,y‑1)与(x,y+1)、(x‑1,y+1)与(x+1,y‑1)、(x‑1,y‑1)与(x+1,y+1);得到像素(x,y)与其相邻八个域像素的位置对应关系;步骤3,计算颜色相离结果;分别计算每个像素(x,y)的四个邻域点对的颜色相离结果,这四个邻域点对的颜色分别进行相离运算,具体表示如下:I<sup>C</sup>(x‑1,y)ΘI<sup>C</sup>(x+1,y)                (2)I<sup>C</sup>(x,y‑1)ΘI<sup>C</sup>(x,y+1)                 (3)I<sup>C</sup>(x‑1,y+1)ΘI<sup>C</sup>(x+1,y‑1)             (4)I<sup>C</sup>(x‑1,y‑1)ΘI<sup>C</sup>(x+1,y+1)              (5)其中,Θ表示颜色相离运算,颜色相离运算是一种双目运算,对于坐标分别为(x<sub>1</sub>,y<sub>1</sub>)和(x<sub>2</sub>,y<sub>2</sub>)的两个像素点,其颜色相离运算定义如下:I<sup>C</sup>(x<sub>1</sub>,y<sub>1</sub>)ΘI<sup>C</sup>(x<sub>2</sub>,y<sub>2</sub>)=DR∪DG∪DB             (6)式中,DR、DG、DB分别为<maths num="0002"><math><![CDATA[<mrow><mi>D</mi><mi>R</mi><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mn>1</mn><mo>,</mo></mrow></mtd><mtd><mrow><msup><mi>I</mi><mi>R</mi></msup><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><msub><mi>y</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><msup><mi>I</mi><mi>R</mi></msup><mrow><mo>(</mo><msub><mi>x</mi><mn>2</mn></msub><mo>,</mo><msub><mi>y</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>&GreaterEqual;</mo><mi>T</mi><mi>H</mi></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><msup><mi>I</mi><mi>R</mi></msup><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><msub><mi>y</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><msup><mi>I</mi><mi>R</mi></msup><mrow><mo>(</mo><msub><mi>x</mi><mn>2</mn></msub><mo>,</mo><msub><mi>y</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>&lt;</mo><mi>T</mi><mi>H</mi></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000972787850000013.GIF" wi="1326" he="131" /></maths><maths num="0003"><math><![CDATA[<mrow><mi>D</mi><mi>G</mi><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mn>1</mn><mo>,</mo></mrow></mtd><mtd><mrow><msup><mi>I</mi><mi>G</mi></msup><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><msub><mi>y</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><msup><mi>I</mi><mi>G</mi></msup><mrow><mo>(</mo><msub><mi>x</mi><mn>2</mn></msub><mo>,</mo><msub><mi>y</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>&GreaterEqual;</mo><mi>T</mi><mi>H</mi></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><msup><mi>I</mi><mi>G</mi></msup><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><msub><mi>y</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><msup><mi>I</mi><mi>G</mi></msup><mrow><mo>(</mo><msub><mi>x</mi><mn>2</mn></msub><mo>,</mo><msub><mi>y</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>&lt;</mo><mi>T</mi><mi>H</mi></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000972787850000014.GIF" wi="1334" he="134" /></maths><maths num="0004"><math><![CDATA[<mrow><mi>D</mi><mi>B</mi><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mn>1</mn><mo>,</mo></mrow></mtd><mtd><mrow><msup><mi>I</mi><mi>B</mi></msup><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><msub><mi>y</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><msup><mi>I</mi><mi>B</mi></msup><mrow><mo>(</mo><msub><mi>x</mi><mn>2</mn></msub><mo>,</mo><msub><mi>y</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>&GreaterEqual;</mo><mi>T</mi><mi>H</mi></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><msup><mi>I</mi><mi>B</mi></msup><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><msub><mi>y</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><msup><mi>I</mi><mi>B</mi></msup><mrow><mo>(</mo><msub><mi>x</mi><mn>2</mn></msub><mo>,</mo><msub><mi>y</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>&lt;</mo><mi>T</mi><mi>H</mi></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000972787850000015.GIF" wi="1327" he="127" /></maths>其中,I<sup>R</sup>(x,y)、I<sup>G</sup>(x,y)和I<sup>B</sup>(x,y)分别代表R、G、B三个颜色通道的像素值,TH是阈值;步骤4,计算颜色差分图;根据每个像素(x,y)的四个邻域点对的颜色相离结果,计算颜色差分图CDM(x,y),计算方法如下:<maths num="0005"><math><![CDATA[<mrow><mi>C</mi><mi>D</mi><mi>M</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mn>255</mn><mo>,</mo></mrow></mtd><mtd><mrow><mi>C</mi><mi>D</mi><mo>=</mo><mn>1</mn></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><mi>C</mi><mi>D</mi><mo>=</mo><mn>0</mn></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000972787850000021.GIF" wi="1294" he="134" /></maths>式中,CD的计算公式为<maths num="0006"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mi>C</mi><mi>D</mi><mo>=</mo><msup><mi>I</mi><mi>C</mi></msup><mrow><mo>(</mo><mi>x</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>y</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><msup><mi>&Theta;I</mi><mi>C</mi></msup><mrow><mo>(</mo><mi>x</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>y</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>&cup;</mo><msup><mi>I</mi><mi>C</mi></msup><mrow><mo>(</mo><mi>x</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>y</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><msup><mi>&Theta;I</mi><mi>C</mi></msup><mrow><mo>(</mo><mi>x</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>y</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>&cup;</mo><msup><mi>I</mi><mi>C</mi></msup><mrow><mo>(</mo><mi>x</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>y</mi><mo>)</mo></mrow><msup><mi>&Theta;I</mi><mi>C</mi></msup><mrow><mo>(</mo><mi>x</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>&cup;</mo><msup><mi>I</mi><mi>C</mi></msup><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><msup><mi>&Theta;I</mi><mi>C</mi></msup><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000972787850000022.GIF" wi="1621" he="135" /></maths>式中,颜色相离运算“Θ”的优先级高于或运算“∪”;计算得到的颜色差分图CDM(x,y)即为最终的边缘检测结果。
地址 100124 北京市朝阳区平乐园100号