主权项 |
1.一种基于FCM聚类的非下采样Contourlet域MRI图像增强方法,其特征在于,包括如下步骤:(1)输入大小为M×N的MRI图像I(m,n),其中1≤m≤M,1≤n≤N,其中M和N均为大于1的自然数,并对其进行L层非下采样Contourlet变换,得到各个尺度上的高频方向子带系数D<sub>l,i</sub>(m,n)和低频子带系数,其中0≤l≤L-1,1≤i≤k<sub>l</sub>,k<sub>l</sub>表示在尺度2<sup>-l</sup>上的高频方向子带的数目,D<sub>l,i</sub>(m,n)表示在尺度2<sup>-l</sup>上的第i个高频方向子带系数,L为3~5,L为自然数;(2)通过对各个尺度上的高频方向子带系数D<sub>l,i</sub>(m,n)分别进行聚类,将其区分为噪声、弱边缘和强边缘系数;(3)利用增益规则公式,对获得了类别属性的高频方向子带系数D<sub>l,i</sub>(m,n)进行修正,即:<img file="FDA00002515660800011.GIF" wi="1464" he="255" />式中,通过设定参数a<sub>0</sub>=0,对噪声系数置零,通过设定参数a<sub>2</sub>=1.25,对强边缘系数进行适度增强,D′<sub>l,i</sub>(m,n)为修正后的高频方向子带系数,A<sub>l</sub>为在尺度2<sup>-l</sup>上的高频方向子带系数的最大模值,a<sub>1</sub>[·]为非线性增益函数,有:<maths num="0001"><![CDATA[<math><mrow><msub><mi>a</mi><mn>1</mn></msub><mo>[</mo><msub><mi>D</mi><mrow><mi>l</mi><mo>,</mo><mi>i</mi></mrow></msub><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>/</mo><msub><mi>A</mi><mi>l</mi></msub><mo>]</mo><mo>=</mo><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><msup><mi>e</mi><mrow><mrow><mo>(</mo><mo>-</mo><mi>t</mi><mo>·</mo><msub><mi>D</mi><mrow><mi>l</mi><mo>,</mo><mi>i</mi></mrow></msub><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>/</mo><msub><mi>A</mi><mi>l</mi></msub></mrow></msup></mrow></mfrac><mo>-</mo><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><msup><mi>e</mi><mrow><mrow><mo>(</mo><mi>t</mi><mo>·</mo><msub><mi>D</mi><mrow><mi>l</mi><mo>,</mo><mi>i</mi></mrow></msub><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>/</mo><msub><mi>A</mi><mi>l</mi></msub></mrow></msup></mrow></mfrac><mo>;</mo></mrow></math>]]></maths>式中,参数t用于控制MRI图像增强的程度,t值越大,MRI图像增强程度越大,为避免MRI图像过增强,设置t=7,e为自然常数,e=2.71828;(4)对步骤(1)中得到的低频子带系数和步骤(3)中得到的修正过的高频方向子带系数作非下采样Contourlet逆变换,即可得到增强的MRI图像I′(m,n)。 |