本发明涉及核磁共振成像技术领域,尤其涉及一种用于核磁共振成像质量控制的数据处理方法。
背景技术:
核磁共振成像(magneticresonanceimaging)是目前最先进的医学成像技术之一,其主要原理是通过不同的射频脉冲信号对生物组织进行激励,利用线圈梯度场对组织进行空间定位,并利用接收线圈检测组织的弛豫时间和质子密度信息,从而形成组织图像。它能够从人体分子内部反映出人体器官失常和早期病变。临床医师通过mri图像对患者的疾病进行诊断,除了和医师本人的临床经验有关外,还与图像中所包含的信息量密切相关。因此对核磁共振成像(qualitycontrol,qc)的质量控制是确保每一个磁共振检查者的生命安全以及疾病得到治疗的根本保证。
目前最常见的质量控制方法是通过磁共振水模体在核磁共振下生成的影像,工作人员通过测试图像得到的数据进行计算分析,进而得到核磁共振成像设备的各种性能指标。而磁共振成像设备得到的图像,一般具有较大的信息量,那么如何从这些数据中得到所需要的信息与参数,是一个长期困扰众多专家学者的难题。
近几年发展起来的支持向量机(supportvectormachines,svm)为解决该类问题提供了一种新的思路。支持向量机建立在统计学习理论的vc维理论和结构风险最小原理基础上,作为支持向量机的扩展,最小二乘支持向量机有效地解决了svm处理大样本数据时速度慢等问题。
核磁共振成像设备所探测的核磁共振信号与每个像素内被测对象某元素的质子密度ρ、表观扩散(弥散)系数、豫驰时间t1及豫驰时间t2有关,这些信息全部存在于核磁共振成像设备所输出的图像中,一般可以表示为:
在实际操作中,上述这些信息均会呈现在mri图像的灰度信息中,却并不能在磁共振设备中给出以上全部的数据,从而给磁共振设备的质量检查造成困难。
技术实现要素:
本发明的目的在于提供一种用于核磁共振成像质量控制的数据处理方法,能够根据磁共振水模体的mri图像,对磁共振成像设备的质量进行检测,保证每一个磁共振检查者的生命安全。
为实现上述目的,本发明采用如下技术方案:
一种用于核磁共振成像质量控制的数据处理方法,依次包括以下步骤:
step1:选择参数区域及层厚,并选定各个参数,包括驰豫时间t1、驰豫时间t2、质子密度ρ及表观扩散系数adc;
step2:获取磁共振水模体的mri图像,并通过图像边缘检测算法提取出磁共振水模体中小球的灰度值;
step3:根据提取出的小球的灰度值分别进行以下计算:
(1)将step2中提取出的小球的灰度值进行小波阈值滤波,并计算小球的平均灰度值与方差,该部分主要目的是获取小球的空域信息;
(2)将step2中提取出的小球的灰度值进行傅里叶变换,得到小球图像的频谱图,根据该频谱图的频率分布情况,截取频谱图中心位置,并求取该部分的平均模长,以便获取小球的频域低频信息;
(3)对小球成像的区域首先求取质心位置,根据质心位置将小球的成像区域分成3x3共9个区域,其中,该9个区域中包含质心位置所在的区域及其灰度值矩阵内数据,将质心所在灰度值矩阵内的数据分别与剩余的八个区域所对应的灰度值矩阵内的数据求相关系数,最后求8个相关系数的平均值;
对于数字图像而言每一个区域对应一个灰度值矩阵,小球所成的像的封闭区域外的灰度值默认为零,封闭区域内的灰度为原mri图像对应坐标的灰度值;
step4:将step3中得到的空域信息、频域低频信息及相关系数,作为最小二乘支持向量机的输入向量,通过由粒子群算法优化的最小二乘支持向量机进行计算,最终输出记录相应数据;
step5:根据以下公式求出偏离率:
其中,e为标准值,表示豫驰时间t1、豫驰时间t2和质子密度ρ中的某一个参数,e′为step4中的计算结果,e′与标准值e表示同一参数;若p小于设定值,则该磁共振成像设备合格,若p大于等于设定值,则该磁共振成像设备不合格,并给出不合格的参数。
优选地,所述step4中,最小二乘支持向量机的惩罚因子和核函数的计算过程如下:
假设输入数据为实数集
则拟合的一般目标是使得残差平方和最小,即求下式的最小值及其参数:
下式为本发明使用的惩罚函数:
这里利用拉格朗日数乘法求解,最终可以得到如下形式:
其中,αk和b均为拉格朗日参数,k(x,xi)表示核函数。
优选地,所述step4中,在粒子群算法中,假设种群中每一个粒子所在搜索空间为n维(n=1,2,…,n),每一个粒子的速度向量和位置向量分别为vi=(vi1,vi2,vi3…vin)和xi=(xi1,xi2,xi3…xin),其中,i表示种群中的第i个粒子,i=1,2,…i,则对粒子i的第d(d=1,2,3......n)维进行更新时的公式为:
xid(t+1)=vid(t)+xid(t)(2)
其中,w为惯性系数,c1和c1分别表示粒子的自身学习因子和社会学习因子,r1和r2是介于0~1之间相互独立的两个随机数,p_bestid(t)和g_bestd(t)分别是粒子i的最好位置和群体最好位置的第d维分量,t表示更新次数;
优选地,所述惯性系数w采用自适应变化的惯性权重系数:
w=(tmax-t)2(wmax-wmin)+wmin
其中,tmax是粒子群算法的最大更新次数,t为当前更新次数,wmax和wmin根据经验选取,wmax=0.9,wmin=0.35。
本发明利用mri图像的空域信息、频域低频信息及相关系数,建立初始最小二乘法支持向量机的模型,然后通过由粒子群算法优化的最小二乘支持向量机进行计算,并记录输出的数据,最后将其结果与对应的核磁共振参数做对比,进而判断核磁共振成像设备是否合格,具有计算量小、计算速度快、计算精度高的优点,能够准确地判断核磁共振成像设备是否合格,从而保证每一个磁共振检查者的生命安全。
附图说明
图1为本发明的流程图;
图2为本发明所述adc层的mri图像;
图3为通过图像边缘检测算法提取后的结果;
图4为本发明所述小球图像的频谱图;
图5为本发明所述小球的成像区域分割示意图。
具体实施方式
以下结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明的部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的其他所有实施例,都属于本发明的保护范围。
如图1至图5所示,本发明所述的一种用于核磁共振成像质量控制的数据处理方法,依次包括以下步骤:
step1:选择参数区域及层厚,并选定各个参数,包括驰豫时间t1、驰豫时间t2、质子密度ρ及表观扩散系数adc;
step2:获取磁共振水模体的mri图像,并通过图像边缘检测算法提取出磁共振水模体中小球的灰度值;
利用图像边缘检测算法提取出磁共振水模体中小球的灰度值的过程为现有技术,不再赘述。
step3:根据提取出的小球的灰度值分别进行以下计算:
(1)将step2中提取出的小球的灰度值进行小波阈值滤波,并计算小球的平均灰度值与方差,该部分主要目的是获取小球的空域信息;
(2)将step2中提取出的小球的灰度值进行傅里叶变换,得到小球图像的频谱图,根据该频谱图的频率分布情况,截取频谱图中心位置,并求取该部分的平均模长,以便获取小球的频域低频信息;
如图4所示,圆框内为截取出的频谱图的中心位置。
(3)对小球成像的区域首先求取质心位置,根据质心位置将小球的成像区域分成3x3共9个区域,其中,该9个区域中包含质心位置所在的区域及其灰度值矩阵内数据,将质心所在灰度值矩阵内的数据分别与剩余的八个区域所对应的灰度值矩阵内的数据求相关系数,最后求8个相关系数的平均值;
对于数字图像而言每一个区域对应一个灰度值矩阵,小球所成的像的封闭区域外的灰度值默认为零,封闭区域内的灰度为原mri图像对应坐标的灰度值;如图5所示,圆点所在位置为小球的质心。
step4:将step3中得到的空域信息、频域低频信息及相关系数,作为最小二乘支持向量机的输入向量,通过由粒子群算法优化的最小二乘支持向量机进行计算,最终输出并记录相应数据;
最小二乘支持向量机的主要思想是利用某种非线性变换,将输入空间数据变换至高维特征空间的问题,本发明主要利用最小二乘支持向量机在数据拟合的作用,最小二乘支持向量机的惩罚因子和核函数的计算过程如下:
假设输入数据为实数集
则拟合的一般目标是使得残差平方和最小,即求下式的最小值及其参数:
下式为本发明使用的惩罚函数:
这里利用拉格朗日数乘法求解,最终可以得到如下形式:
其中,αk和b均为拉格朗日参数,k(x,xi)表示核函数,表1为几种核函数的形式。
表1
在最小二乘支持向量机中核函数的选取和惩罚因子的选取是非常重要的,而一般的做法是通过网格搜索,交叉验证的方式来达到搜寻最优参数的目的,这种做法具有极大的盲目性,因此会带来较大的计算量。而粒子群算法在众多优化算法中具有结构简单,全局能力强等诸多优点,因此本发明采用了一种自适应调整惯性权重的粒子群算法来搜索最小二乘支持向量机的最优参数。
粒子群优化算法是受鸟类觅食这类社会行为的启发而提出的,主要通过解空间中各个解向量之间合作-竞争搜索机制来搜索全局最优解。在计算时,假设种群中每一个粒子所在搜索空间为n维(n=1,2,…,n),每一个粒子的速度向量和位置向量分别为vi=(vi1,vi2,vi3…vin)和xi=(xi1,xi2,xi3…xin),其中,i表示种群中的第i个粒子,i=1,2,…i,则对粒子i的第d(d=1,2,3......n)进行更新时的公式为:
xid(t+1)=vid(t)+xid(t)(2)
其中,w为惯性系数,c1和c1分别表示粒子的自身学习因子和社会学习因子,r1和r2是介于0~1之间相互独立的两个随机数,
惯性系数w较大时,全局搜索能力强而局部搜索能力弱,较小时反之,故此在粒子群算法中将其设置为一个常量是不可取的。本发明中,惯性系数w采用自适应变化的惯性权重系数:
w=(tmax-t)2(wmax-wmin)+wmin
其中,tmax是粒子群算法的最大更新次数,t为当前更新次数,wmax和wmin根据经验选取,wmax=0.9,wmin=0.35。
其中,最小二乘支持向量机、粒子群优化算法,以及由粒子群算法优化的最小二乘支持向量机均为现有技术,以上均为简单概述,在此不再赘述。
step5:根据以下公式求出偏离率:
其中,e为标准值,表示豫驰时间t1、豫驰时间t2和质子密度ρ中的某一个参数,e′为step4中的计算结果,e′与标准值e表示同一参数;若p小于设定值,则该磁共振成像设备合格,若p大于等于设定值,则该磁共振成像设备不合格,并给出不合格的参数。
本发明具有计算量小、计算速度快、计算精度高的优点,能够准确地判断核磁共振成像设备是否合格,从而保证每一个磁共振检查者的生命安全。