主权项 |
1、一种融合上下文信息的数字血管造影图像增强方法,其步骤为:(1)用哈勃滤波器将灰度反转的数字血管造影图像分解至K个方向L个尺度,K的取值范围为4~16,L根据最大血管宽度W<sub>max</sub>与最小血管宽度W<sub>min</sub>确定;(2)将哈勃滤波器的输出作为增强模型的初始输入,并对每个方向和尺度的哈勃滤波响应实施基于共圆约束的增强;其过程为:(2.1)将中心位于(x,y)的像素在方向α、中心频率ω分量上受到来自环境的增强作用F(x,y,α,ω)表示为:<maths num="0001"><![CDATA[<math><mrow><mi>F</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>α</mi><mo>,</mo><mi>ω</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>Σ</mi><mrow><mrow><mo>(</mo><msup><mi>x</mi><mo>′</mo></msup><mo>,</mo><msup><mi>y</mi><mo>′</mo></msup><mo>)</mo></mrow><mo>∈</mo><msub><mi>A</mi><mi>e</mi></msub></mrow></munder><munder><mi>Σ</mi><mi>β</mi></munder><msub><mi>W</mi><mi>c</mi></msub><mrow><mo>(</mo><msup><mi>x</mi><mo>′</mo></msup><mo>,</mo><msup><mi>y</mi><mo>′</mo></msup><mo>,</mo><mi>β</mi><mo>;</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>α</mi><mo>)</mo></mrow><msub><mi>W</mi><mi>d</mi></msub><mrow><mo>(</mo><msup><mi>x</mi><mo>′</mo></msup><mo>,</mo><msup><mi>y</mi><mo>′</mo></msup><mo>;</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mi>R</mi><mrow><mo>(</mo><msup><mi>x</mi><mo>′</mo></msup><mo>,</mo><msup><mi>y</mi><mo>′</mo></msup><mo>,</mo><mi>β</mi><mo>,</mo><mi>ω</mi><mo>)</mo></mrow></mrow></math>]]></maths>式中A<sub>e</sub>表示满足增强特性的区域,R(x′,y′,β,ω)表示位于(x′,y′)的像素在中心频率为ω、方向为β上的响应,式中的<maths num="0002"><![CDATA[<math><mrow><msub><mi>W</mi><mi>c</mi></msub><mrow><mo>(</mo><msup><mi>x</mi><mo>′</mo></msup><mo>,</mo><msup><mi>y</mi><mo>′</mo></msup><mo>,</mo><mi>β</mi><mo>;</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>α</mi><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><msup><mi>Q</mi><mn>2</mn></msup><msubsup><mi>σ</mi><mi>c</mi><mn>2</mn></msubsup></mfrac><mo>)</mo></mrow></mrow></math>]]></maths>为曲率加权函数,Q为曲率,且<maths num="0003"><![CDATA[<math><mrow><mi>Q</mi><mo>=</mo><mfrac><mn>1</mn><mi>r</mi></mfrac><mo>=</mo><mfrac><mn>2</mn><mi>d</mi></mfrac><mi>sin</mi><mrow><mo>(</mo><mi>φ</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mfrac><mn>2</mn><mi>d</mi></mfrac><mi>sin</mi><mo>|</mo><mfrac><mrow><mi>β</mi><mo>-</mo><mi>α</mi></mrow><mn>2</mn></mfrac><mo>|</mo></mtd><mtd><mn>0</mn><mo>≤</mo><mn>2</mn><mi>γ</mi><mo>-</mo><mi>α</mi><mo><</mo><mi>π</mi></mtd></mtr><mtr><mtd><mfrac><mn>2</mn><mi>d</mi></mfrac><mi>cos</mi><mo>|</mo><mfrac><mrow><mi>β</mi><mo>-</mo><mi>α</mi></mrow><mn>2</mn></mfrac><mo>|</mo></mtd><mtd><mn>2</mn><mi>γ</mi><mo>-</mo><mi>α</mi><mo><</mo><mn>0</mn><mi>or</mi><mn>2</mn><mi>γ</mi><mo>-</mo><mi>α</mi><mo>≥</mo><mi>π</mi></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>式中<maths num="0004"><![CDATA[<math><mrow><mi>d</mi><mo>=</mo><msqrt><msup><mrow><mo>(</mo><msup><mi>x</mi><mo>′</mo></msup><mo>-</mo><mi>x</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msup><mi>y</mi><mo>′</mo></msup><mo>-</mo><mi>y</mi><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><msub><mi>W</mi><mi>d</mi></msub><mrow><mo>(</mo><msup><mi>x</mi><mo>′</mo></msup><mo>,</mo><msup><mi>y</mi><mo>′</mo></msup><mo>;</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mo>[</mo><mo>-</mo><mfrac><mrow><msup><mrow><mo>(</mo><msup><mi>x</mi><mo>′</mo></msup><mo>-</mo><mi>x</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msup><mi>y</mi><mo>′</mo></msup><mo>-</mo><mi>y</mi><mo>)</mo></mrow><mn>2</mn></msup></mrow><msubsup><mi>σ</mi><mi>d</mi><mn>2</mn></msubsup></mfrac><mo>]</mo></mrow></math>]]></maths>为距离加权函数;(2.2)通过S形变换进行迭代处理:<img file="S2007101686634C00016.GIF" wi="896" he="132" />式中s表示增强作用的上限饱和系数,取值范围为20~80,<img file="S2007101686634C00017.GIF" wi="246" he="58" />的初始值为图像经过哈勃滤波器的输出值;(2.3)选每个像素所有方向和尺度上的响应值之中最大的一个作为模型的输出,即:<maths num="0006"><![CDATA[<math><mrow><mi>R</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>max</mi><mrow><mi>ω</mi><mo>,</mo><mi>α</mi></mrow></munder><mo>[</mo><mi>R</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>α</mi><mo>,</mo><mi>ω</mi><mo>)</mo></mrow><mo>]</mo></mrow></math>]]></maths>(3)运用各向异性扩散去噪方法对步骤(2)的结果进行平滑去噪处理,得到增强后的图像。 |