一种基于高斯滤波的高维脑核磁图像多阈值分割方法技术

本发明专利技术提供了一种基于高斯滤波的高维脑核磁图像多阈值分割方法，属于智慧医疗技术领域。其技术方案为：首先，获取脑核磁图像数据F；其次，通过标准差为脑核磁图像的粗糙性度量的高斯滤波进行自适应平滑滤波预处理；再次，对预处理后的脑核磁图像进行灰度直方图统计，并根据灰度直方图的峰值设定脑部组织的灰质、白质、脑脊液和背景的三个初始分割阈值(t1,t2,t3)；最后，将脑核磁图像中三个目标组织和背景的四类间方差σ2(t1,t2,t3)作为混合蛙跳算法的适应度函数寻找三个最优分割阈值并输出分割后的二值化图像。本发明专利技术的有益效果为：降低了噪声对脑核磁图像分割的影响，提高了对脑核磁图像中三个目标组织的分割精度，对脑核磁图像智能辅助分割和诊断具有较强的应用价值。用价值。用价值。

【技术实现步骤摘要】
一种基于高斯滤波的高维脑核磁图像多阈值分割方法


[0001]本专利技术涉及智慧医疗
，尤其涉及一种基于高斯滤波的高维脑核磁图像多阈值分割方法。

技术介绍

[0002]核磁共振成像技术已经成为脑疾病临床诊断的重要辅助手段。脑核磁与电子计算断层扫描相比，它具有无放射线损害，无骨性伪影，能多方面、多参数成像，有高度的软组织分辨能力，不需要使用对比剂即可显示血管结构等独特的优点。然而在实际应用中，由于射频场的不均匀性等因素，导致脑核磁图像的灰度均匀性变差；其次，脑核磁图像成像过程中由于受仪器设备等物理原因影响，使得图像中经常含有噪声。这给医生的临床诊断工作增加了不确定性因素和工作量。同时，这也影响了脑核磁图像的分割精度。如何有效地提高脑核磁图像的分割精度和帮助医生有效地分析患者脑部组织的临床症状，亟需一种新的方法能有效地减少脑核磁图像中的噪声和提高分割精度，有效进行脑核磁图像临床症状的筛查和分析。
[0003]如何解决上述技术问题为本专利技术面临的课题。

技术实现思路

[0004]本专利技术的目的在于提供一种基于高斯滤波的高维脑核磁图像多阈值分割方法。
[0005]本专利技术是通过如下措施实现的：一种基于高斯滤波的高维脑核磁图像多阈值分割方法，包括以下步骤：
[0006]步骤1、读取脑核磁图像F，其中图像的大小为256
×
256，单张的脑核磁图像中共有256
×
256＝65536个像素点，可表示为论域T＝{F(0,0),F(0,1),...,F(255,255)}，其中F(i,j)(i,j∈[0，255])表示第(i,j)个像素点的灰度级且F(i,j)∈[0,255]；
[0007]步骤2、计算脑核磁图像的粗糙性度量ρ，并将ρ作为高斯滤波的标准差σ对脑核磁图像进行自适应平滑滤波预处理；
[0008]步骤3、对预处理后的脑核磁图像进行灰度直方图统计，根据灰度直方图的峰值设定脑部组织的灰质、白质、脑脊液和背景的三个初始分割阈值(t1,t2,t3)；
[0009]步骤4、计算脑核磁图像中三个目标组织和背景的四类间方差σ2(t1,t2,t3)，并将其作为混合蛙跳算方法的适应度函数；
[0010]步骤5、通过混合蛙跳算法寻找脑核磁图像中脑部组织的灰质、白质、脑脊液和背景的三个最优分割阈值
[0011]步骤6、对脑核磁图像进行多阈值分割，输出分割后的二值化图像。
[0012]作为本专利技术提供的一种基于高斯滤波的高维脑核磁图像多阈值分割方法进一步优化方案，所述步骤2的具体步骤如下：
[0013]步骤2.1、计算脑核磁图像F在不同灰度级上分布信息下近似H(l)，其计算公式如下：
[0014][0015]其中，F(i,j)为像素(i,j)的灰度级，δ(
·
)近似为零点处高度为1的冲激函数，L为图像的灰度级数目；
[0016]步骤2.2、计算脑核磁图像中两个像素点F1：F(i1,j1)和F2：F(i2,j2)之间的空间距离d(F1,F2)，其计算公式如下：
[0017][0018]步骤2.3、计算脑核磁图像中像素F(i,j)的3＝3规模邻域，邻域中的各邻接像素与F(i,j)之间的空间距离d3×3(i,j)，其计算公式如下：
[0019][0020]步骤2.4、计算脑核磁图像中3
×
3规模邻域相对于像素F(i,j)具有的同质程度S(i,j)，其计算公式如下：
[0021][0022]其中，邻域同质程度函数S(i,j)是一个Cauchy型分布函数，r参数表示不可区分色差的阈值；
[0023]步骤2.5、根据像素邻域通知性度量，计算脑核磁图像F在不同灰度级上分布信息的上近似其计算公式如下：
[0024][0025]步骤2.6、根据脑核磁图像F在不同灰度级上分布信息的下近似H(l)和上近似信息，在不同灰度级上的粗糙量化信息的计算公式如下：
[0026][0027]其中，为灰度值l上的分布粗糙度，W(H(l))为权重函数，设参数0＜b＜c，权重函数的计算公式如下：
[0028][0029]其中，k1＝c2/2b2，k2＝c2/2(c
‑
b)2，参数b，c的取值与统计直方图分布均值的比例分别为0.05和0.5；
[0030]步骤2.7、将脑核磁图像F在不同灰度级上的粗糙量化信息ρ(l)作为高斯滤波的标准差σ对不同灰度级的像素点进行自适应平滑滤波处理，得到滤波后的图像F'。
[0031]作为本专利技术提供的一种基于高斯滤波的高维脑核磁图像多阈值分割方法进一步优化方案，所述步骤5的具体步骤如下：
[0032]步骤5.1、根据脑核磁图像的灰度级数目随机生成青蛙群体，对其进行8位二进制编码，随机选取P＝M
×
N只作为初始青蛙群体，第k只青蛙个体记为x
k
＝(x
k1
,x
k2
,...,x
ks
)，其中s＝3为分割阈值的个数，M为族群数，N为族群内青蛙个数；
[0033]步骤5.2、将脑核磁图像的四类间方差作为混合蛙跳算法的适应度函数f(x
k
)，并计算每只青蛙个体的适应值，其计算公式如下：
[0034][0035]其中，ω
p
(p＝0,1,2,3)表示脑核磁图像中脑组织的灰质、白质、脑脊液和背景四个部分的概率分布，μ
p
(p＝0,1,2,3)表示脑核磁图像中脑组织的灰质、白质、脑脊液和背景四个部分的灰度均值，表示脑核磁图像总体的灰度均值；
[0036]步骤5.3、按照所有青蛙个体的适应度f(x
k
)降序排序，记全局最优青蛙个体为X
g
，并逐一循环分配给M个族群；
[0037]步骤5.4、对每个青蛙种群进行局部搜索，种群中最优适应值和最差适应值的青蛙个体分别记为X
mb
和X
mw
，其中m∈[0,M]表示第m个青蛙种群，然后，对青蛙种群最差适应值的青蛙个体进行更新，更新规则如下：
[0038]D＝rand()
×
(X
mb
‑
X
mw
)
ꢀꢀꢀ
(9)
[0039]X'
mw
＝X
mw
+D,||D||≤D
max
ꢀꢀꢀ
(10)
[0040]其中，rand()为0到1之间的随机数，D为青蛙个体的移动步长，D
max
表示青蛙个体的最大移动步长；
[0041]步骤5.5、判断更新生成的新青蛙个体X'
mw
与原青蛙个体X
mw
之间的大小关系；如果
更新后X'
mw
的适应值优于X
mw
，那么更新成功，否则，用X
g
代替X本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于高斯滤波的高维脑核磁图像多阈值分割方法，其特征在于，包括以下步骤：步骤1、读取脑核磁图像F，其中图像的大小为256
×
256，单张的脑核磁图像中共有256
×
256＝65536个像素点，可表示为论域T＝{F(0,0),F(0,1),...,F(255,255)}，其中F(i,j)(i,j∈[0，255])表示第(i,j)个像素点的灰度级且F(i,j)∈[0,255]；步骤2、计算脑核磁图像的粗糙性度量ρ，并将ρ作为高斯滤波的标准差σ对脑核磁图像进行自适应平滑滤波预处理；步骤3、对预处理后的脑核磁图像进行灰度直方图统计，根据灰度直方图的峰值设定脑部组织的灰质、白质、脑脊液和背景的三个初始分割阈值(t1,t2,t3)；步骤4、计算脑核磁图像中三个目标组织和背景的四类间方差σ2(t1,t2,t3)，并将其作为混合蛙跳算方法的适应度函数；步骤5、通过混合蛙跳算法寻找脑核磁图像中脑部组织的灰质、白质、脑脊液和背景的三个最优分割阈值步骤6、对脑核磁图像进行多阈值分割，输出分割后的二值化图像。2.根据权利要求1所述的基于高斯滤波的高维脑核磁图像多阈值分割方法，其特征在于，所述步骤2的具体步骤如下：步骤2.1、计算脑核磁图像F在不同灰度级上分布信息下近似H(l)，其计算公式如下：其中，F(i,j)为像素(i,j)的灰度级，δ(
·
)近似为零点处高度为1的冲激函数，L为图像的灰度级数目；步骤2.2、计算脑核磁图像中两个像素点F1：F(i1,j1)和F2：F(i2,j2)之间的空间距离d(F1,F2)，其计算公式如下：步骤2.3、计算脑核磁图像中像素F(i,j)的3
×
3规模邻域，邻域中的各邻接像素与F(i,j)之间的空间距离d3×3(i,j)，其计算公式如下：步骤2.4、计算脑核磁图像中3
×
3规模邻域相对于像素F(i,j)具有的同质程度S(i,j)，其计算公式如下：其中，邻域同质程度函数S(i,j)是一个Cauchy型分布函数，r参数表示不可区分色差的阈值；
步骤2.5、根据像素邻域通知性度量，计算脑核磁图像F在不同灰度级上分布信息的上近似其计算公式如下：步骤2.6、根据脑核磁图像F在不同灰度级上分布信息的下近似H(l)和上近似信息，在不同灰度级上的粗糙量化信息的计算公式如下：其中，为灰度值l上的分布粗糙度，W(H(l))为权重函数，设参数0＜b＜c，权重函数的计算公式如下：其中，k1＝c2/2b2，k2＝c2/2(c
‑
b)2，参数b，c的取值与统计直方图分布均值的比例分别为0.05和0.5；步骤2.7、将脑核磁图像F在不同灰度级上的粗糙量化信息ρ(l)作为高斯滤波的标准差σ对不同灰度级的像素点进行自适应平滑滤波处理，得到滤波后的图像F'。3.根据权利要求1所述的基于高斯滤波的高维脑核磁图像多阈值分割方法，其特征在于，所述步骤5的具体步骤如下：步骤5.1、...

【专利技术属性】
技术研发人员：丁卫平，李铭，沈子宁，孙颖，鞠恒荣，黄嘉爽，曹金鑫，程纯，沈鑫杰，王海鹏，孙楚迪，
申请(专利权)人：南通大学，
类型：发明
国别省市：