发明名称 基于FCM聚类的非下采样Contourlet域MRI图像增强方法
摘要 本发明方法涉及在非下采样Contourlet域上,利用模糊C均值Fuzzy C-Means,FCM聚类,对磁共振成像Magnetic Resonance Imaging,即MRI图像进行增强的方法,尤其是一种基于FCM聚类的非下采样Contourlet域MRI图像增强方法。本发明方法与现有技术相比具有以下优点:1.本发明由于采用非下采样Contourlet变换,可以有效捕捉MRI图像中的边缘、轮廓信息和纹理细节信息。2.本发明由于引入FCM聚类的方法自适应地将高频方向子带系数区分为噪声、弱边缘和强边缘系数,克服了现有的基于变换域的图像增强方法对阈值选取的依赖。
申请公布号 CN102968779A 申请公布日期 2013.03.13
申请号 CN201210508617.5 申请日期 2012.12.03
申请人 北方民族大学 发明人 常霞;高岳林;黄永东;纪峰
分类号 G06T5/50(2006.01)I 主分类号 G06T5/50(2006.01)I
代理机构 宁夏专利服务中心 64100 代理人 赵明辉
主权项 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>&CenterDot;</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>&CenterDot;</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)。
地址 750021 宁夏回族自治区银川市西夏区文昌北路204号