[发明专利]一种磁共振灌注成像后处理方法及系统

摘要

本发明提供了一种磁共振灌注成像后处理方法及系统，方法包括：对原始脑灌注图像滤波，去除脑部边缘，得到滤波后脑灌注图像；接收用户的区域选择指令，获取在滤波后脑灌注图像的中动脉区的选定区域，输出评估动脉输入函数的加权动脉输入曲线；根据改进的Gamma函数分别对加权动脉输入曲线和浓度时间曲线进行拟合，得到优化的动脉输入曲线和浓度时间曲线；根据改进的基于非参数模型的奇异矩阵分解方法求解灌注的相对定量参数。本发明通过对AIF和求解矩阵进行加权优化，降低对噪声的敏感度，采用简化且有效的Gamma函数进行拟合，将非线性问题转化为线性求解，加快了后处理速度。

主权项

1.一种磁共振灌注成像后处理方法，其特征在于，所述方法包括以下步骤：A、对原始脑灌注图像滤波，去除脑部边缘，得到滤波后脑灌注图像；B、接收用户的区域选择指令，获取在滤波后脑灌注图像的中动脉区的选定区域，输出评估动脉输入函数的加权动脉输入曲线；C、根据改进的Gamma函数分别对加权动脉输入曲线和浓度时间曲线进行拟合，得到优化的动脉输入曲线和浓度时间曲线；D、根据改进的基于非参数模型的奇异矩阵分解方法求解灌注的相对定量参数；所述步骤C具体包括：C1、根据浓度时间曲线的最大值、上升期、下降期估算改进的Gamma函数拟合的初始值；C2、在[T0,Tmax+n]曲线范围内结合Gamma函数拟合的初始值得到改进的Gamma函数的优化参数；其中T0是初始到达时间，Tmax是最大峰值对应的时间，n是峰值后参与拟合的点数；C3、根据该改进的Gamma函数对加权动脉输入曲线和浓度时间曲线拟合，得到优化的动脉输入曲线和浓度时间曲线；Gamma函数表达式如下：G(t)＝A*(t‑D)B*e(‑(t‑D)/C)其中A，D，B，C是需要估算的优化参数；t是时间变量；Gamma函数的初始值[A0,B0,C0,D0]估计如下：A0取浓度时间曲线Ct的最大峰值；D0取浓度时间曲线Ct的第一个拐点位置即对比剂开始进入组织的时间；C0取浓度时间曲线Ct在下降时刻与D0相对应的时间；B0＝(Tmax‑T0)/C0；Tmax是A0对应的时间。