发明名称 一种利用多尺度多结构元素数学形态学的图像增强方法
摘要 本发明涉及一种利用多尺度多结构元素数学形态学的图像增强方法,它有四大步骤:设有n个尺度、m个不同形状的结构元素参与计算,首先,利用数学形态学中的高帽变换计算每个尺度下m个不同形状结构元素对应的图像区域;然后,通过累加每个结构元素对应的图像区域得到每个尺度下由m个结构元素提取的亮区域和暗区域:随后,取所有尺度下亮区域和暗区域的最大值作为所有尺度下所有结构元素提取的亮区域和暗区域;最后,通过在原始图像上分别加上提取的亮区域并减去提取的暗区域增大原始图像上亮区域和暗区域之间的对比度,从而达到图像增强的目的。本发明广泛应用于各类基于图像的应用系统,具有广阔的市场前景与应用价值。
申请公布号 CN101807294A 申请公布日期 2010.08.18
申请号 CN201010105996.4 申请日期 2010.02.03
申请人 北京航空航天大学 发明人 白相志;周付根
分类号 G06T5/20(2006.01)I 主分类号 G06T5/20(2006.01)I
代理机构 北京慧泉知识产权代理有限公司 11232 代理人 王顺荣;唐爱华
主权项 1.一种利用多尺度多结构元素数学形态学的图像增强方法,其特征在于:该方法具体步骤如下:设有n个尺度、m个不同形状的结构元素参与计算:步骤一:计算每个尺度j(1≤j≤n)下m个不同形状结构元素B<sub>j</sub><sup>1</sup>,B<sub>j</sub><sup>2</sup>,……,B<sub>j</sub><sup>m</sup>对应的高帽变换:<img file="FSA00000014417100011.GIF" wi="390" he="52" /><img file="FSA00000014417100012.GIF" wi="390" he="53" /><maths num="0001"><![CDATA[<math><mrow><mo>.</mo><mo>.</mo><mo>.</mo></mrow></math>]]></maths><img file="FSA00000014417100014.GIF" wi="394" he="53" /><maths num="0002"><![CDATA[<math><mrow><msubsup><mi>BTH</mi><mi>j</mi><mn>1</mn></msubsup><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><mi>f</mi><mo>&CenterDot;</mo><msubsup><mi>B</mi><mi>j</mi><mn>1</mn></msubsup><mo>-</mo><mi>f</mi><mo>,</mo></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msubsup><mi>BTH</mi><mi>j</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><mi>f</mi><mo>&CenterDot;</mo><msubsup><mi>B</mi><mi>j</mi><mn>2</mn></msubsup><mo>-</mo><mi>f</mi><mo>,</mo></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><mo>.</mo><mo>.</mo><mo>.</mo></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><msubsup><mi>BTH</mi><mi>j</mi><mi>m</mi></msubsup><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><mi>f</mi><mo>&CenterDot;</mo><msubsup><mi>B</mi><mi>j</mi><mi>m</mi></msubsup><mo>-</mo><mi>f</mi><mo>.</mo></mrow></math>]]></maths>其中,WTH(x,y)=f(x,少)-f○B(x,y),BTH(x,y)=f●B(x,少)-f(x,y),<img file="FSA00000014417100019.GIF" wi="403" he="49" /><maths num="0006"><![CDATA[<math><mrow><mi>f</mi><mo>&CenterDot;</mo><mi>B</mi><mo>=</mo><mrow><mo>(</mo><mi>f</mi><mo>&CirclePlus;</mo><mi>B</mi><mo>)</mo></mrow><mi>B</mi><mo>,</mo></mrow></math>]]></maths><maths num="0007"><![CDATA[<math><mrow><mi>f</mi><mo>&CirclePlus;</mo><mi>B</mi><mo>=</mo><munder><mi>max</mi><mrow><mi>u</mi><mo>,</mo><mi>v</mi></mrow></munder><mrow><mo>(</mo><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>-</mo><mi>u</mi><mo>,</mo><mi>y</mi><mo>-</mo><mi>v</mi><mo>)</mo></mrow><mo>+</mo><mi>B</mi><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mi>v</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths><img file="FSA000000144171000112.GIF" wi="736" he="79" />f是原始图像;B是结构元素;(x,y)是图像像素坐标;(u,v)是结构元素中的像素坐标;WTH是白高帽变换;BTH是黑高帽变换;○,●,<img file="FSA000000144171000113.GIF" wi="40" he="50" />和<img file="FSA000000144171000114.GIF" wi="49" he="53" />分别是数学形态学开、闭、膨胀和腐蚀运算符;步骤二:按下式计算第j尺度下由m个结构元素提取的亮区域(MSEWTH)和暗区域(MSEBTH):<maths num="0008"><![CDATA[<math><mrow><mi>MSEWT</mi><msub><mi>H</mi><mi>j</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><msup><mi>WTH</mi><mi>i</mi></msup><mi>j</mi></msub><mo>,</mo></mrow></math>]]></maths><maths num="0009"><![CDATA[<math><mrow><mi>MSEBT</mi><msub><mi>H</mi><mi>j</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><msup><mi>BTH</mi><mi>i</mi></msup><mi>j</mi></msub><mo>;</mo></mrow></math>]]></maths>步骤三:按下式计算所有n尺度下由m个结构元素提取的亮区域(MSMSEWTH)和暗区域(MSMSEBTH):<maths num="0010"><![CDATA[<math><mrow><mi>MSMSEWTH</mi><mo>=</mo><munder><mi>max</mi><mrow><mn>0</mn><mo>&le;</mo><mi>j</mi><mo>&le;</mo><mi>n</mi></mrow></munder><mo>{</mo><msub><mi>MSEWTH</mi><mi>j</mi></msub><mo>}</mo><mo>,</mo></mrow></math>]]></maths><maths num="0011"><![CDATA[<math><mrow><mi>MSMSEBTH</mi><mo>=</mo><munder><mi>max</mi><mrow><mn>0</mn><mo>&le;</mo><mi>j</mi><mo>&le;</mo><mi>n</mi></mrow></munder><mo>{</mo><msub><mi>MSEBTH</mi><mi>j</mi></msub><mo>}</mo><mo>;</mo></mrow></math>]]></maths>步骤四:按下式通过增大原始图像上亮区域和暗区域之间的对比度达到图像增强的目的:f<sub>En</sub>=f+MSMSWTH-MSMSEBTH其中,f<sub>En</sub>是最终增强的图像,f是原始图像。
地址 100191 北京市海淀区学院路37号北航宇航学院图像中心