快了后处理速度。
[0096] 基于上述方法实施例，本发明还提供了一种磁共振灌注成像后处理系统。如图4所 示，所述磁共振灌注成像后处理系统包括：
[0097] 滤波模块100,用于对原始脑灌注图像滤波，去除脑部边缘，得到滤波后脑灌注图 像；
[0098] 区域选定及输出模块200,用于接收用户的区域选择指令，获取在滤波后脑灌注图 像的中动脉区的选定区域，输出评估动脉输入函数的加权动脉输入曲线；
[0099] 曲线拟合模块300,用于根据改进的Gamma函数分别对加权动脉输入曲线和浓度时 间曲线进行拟合，得到优化的动脉输入曲线和浓度时间曲线；
[0100] 参数获取模块400,用于根据改进的基于非参数模型的奇异矩阵分解系统求解灌 注的相对定量参数。
[0101] 进一步的，在所述磁共振灌注成像后处理系统中，所述区域选定及输出模块200具 体包括：
[0102] 区域选定单元，用于接收用户的区域选择指令，获取在滤波后脑灌注图像的中动 脉区的选定区域；
[0103] 加权输出单元，用于对中动脉区的选定区域的浓度时间信号进行加权平滑，得到 用于参数计算的动脉输入函数，即输出加权动脉输入曲线。
[0104] 进一步的，在所述磁共振灌注成像后处理系统中，所述曲线拟合模块300具体包 括：
[0105] 估算单元，用于根据浓度时间曲线的最大值、上升期、下降期估算改进的Gamma函 数拟合的初始值；
[010W Gamma函数参数获取单元，用于化To,Tmax+n]范围内，结合估算单元输出的初始值 获取用于曲线拟合的Gamma函数参数。；其中To是初始到达时间，Tmax是最大峰值对应的时 间，n是峰值后参与拟合的点数。
[0107] 曲线拟合单元，用于对加权动脉输入曲线和浓度时间曲线拟合，得到优化的动脉 输入曲线和浓度时间曲线。
[0108] 进一步的，在所述磁共振灌注成像后处理系统中，所述参数获取模块400具体包 括：
[0109] 分解单元，用于根据改进的基于非参数模型的奇异矩阵分解方法中的矩阵A进行 加权组合，并对矩阵A进行SVD分解得到对角矩阵W，正交矩阵V，S角矩阵U;
[0110] 相对定量参数获取单元，用于根据阔值策略得到处理后的对角矩阵S，并根据对角 矩阵S确定灌注的相对定量参数;其中S = l/W。
[0111] 综上所述，本发明提供了一种磁共振灌注成像后处理方法及系统，方法包括:对原 始脑灌注图像滤波，去除脑部边缘，得到滤波后脑灌注图像;接收用户的区域选择指令，获 取在滤波后脑灌注图像的中动脉区的选定区域，输出评估动脉输入函数的加权动脉输入曲 线;根据改进的Gamma函数分别对加权动脉输入曲线和浓度时间曲线进行拟合，得到优化的 动脉输入曲线和浓度时间曲线根据改进的基于非参数模型的奇异矩阵分解方法求解灌注 的相对定量参数。本发明通过对AIF和求解矩阵进行加权优化，降低对噪声的敏感度，采用 简化且有效的Gamma函数进行拟合，将非线性问题转化为线性求解，加快了后处理速度。
[0112] 可W理解的是，对本领域普通技术人员来说，可W根据本发明的技术方案及本发 明构思加 W等同替换或改变，而所有运些改变或替换都应属于本发明所附的权利要求的保 护范围。
【主权项】
1. 一种磁共振灌注成像后处理方法，其特征在于，所述方法包括以下步骤： A、 对原始脑灌注图像滤波，去除脑部边缘，得到滤波后脑灌注图像； B、 接收用户的区域选择指令，获取在滤波后脑灌注图像的中动脉区的选定区域，输出 评估动脉输入函数的加权动脉输入曲线； C、 根据改进的Gamma函数分别对加权动脉输入曲线和浓度时间曲线进行拟合，得到优 化的动脉输入曲线和浓度时间曲线； D、 根据改进的基于非参数模型的奇异矩阵分解方法求解灌注的相对定量参数。2. 根据权利要求1所述磁共振灌注成像后处理方法，其特征在于，所述步骤B具体包括： B1、接收用户的区域选择指令，获取在滤波后脑灌注图像的中动脉区的选定区域； B2、对中动脉区的选定区域的浓度时间信号进行加权平滑，得到用于参数计算的动脉 输入函数，并根据动脉输入函数输出加权动脉输入曲线。3. 根据权利要求1所述磁共振灌注成像后处理方法，其特征在于，所述步骤C具体包括： C1、根据浓度时间曲线的最大值、上升期、下降期估算改进的Gamma函数拟合的初始值； C2、在[To，Tmax+n ]曲线范围内结合Gamma函数拟合的初始值得到改进的Gamma函数的优 化参数;其中To是初始到达时间，Tmax是最大峰值对应的时间，η是峰值后参与拟合的点数； C3、根据该改进的Gamma函数对加权动脉输入曲线和浓度时间曲线拟合，得到优化的动 脉输入曲线和浓度时间曲线。4. 根据权利要求1所述磁共振灌注成像后处理方法，其特征在于，所述步骤D具体包括： D1、根据改进的基于非参数模型的奇异矩阵分解方法中的矩阵A进行加权组合，并对矩 阵A进行SVD分解得到对角矩阵W，正交矩阵V，三角矩阵U; D2、根据阈值策略得到处理后的对角矩阵S，并根据对角矩阵S确定灌注的相对定量参 数;其中S = l/W。5. -种磁共振灌注成像后处理系统，其特征在于，包括： 滤波模块，用于对原始脑灌注图像滤波，去除脑部边缘，得到滤波后脑灌注图像； 区域选定及输出模块，用于接收用户的区域选择指令，获取在滤波后脑灌注图像的中 动脉区的选定区域，输出评估动脉输入函数的加权动脉输入曲线； 曲线拟合模块，用于根据改进的Gamma函数分别对加权动脉输入曲线和浓度时间曲线 进行拟合，得到优化的动脉输入曲线和浓度时间曲线； 参数获取模块，用于根据改进的基于非参数模型的奇异矩阵分解系统求解灌注的相对 定量参数。6. 根据权利要求5所述磁共振灌注成像后处理系统，其特征在于，所述区域选定及输出 模块具体包括： 区域选定单元，用于接收用户的区域选择指令，获取在滤波后脑灌注图像的中动脉区 的选定区域； 加权输出单元，用于对中动脉区的选定区域的浓度时间信号进行加权平滑，得到用于 参数计算的动脉输入函数;其中，所述动脉输入函数为加权动脉输入曲线。7. 根据权利要求5所述磁共振灌注成像后处理系统，其特征在于，所述曲线拟合模块具 体包括： 估算单元，用于根据浓度时间曲线的最大值、上升期、下降期估算改进的Gamma函数拟 合的初始值； Gamma函数参数获取单元，用于在[To, Wn]范围内，结合估算单元输出的初始值获取 用于曲线拟合的Gamma函数参数;其中To是初始到达时间，Tmax是最大峰值对应的时间，η是 峰值后参与拟合的点数。 曲线拟合单元，用于根据该改进的Gamma函数对加权动脉输入曲线和浓度时间曲线拟 合，得到优化的动脉输入曲线和浓度时间曲线。8.根据权利要求5所述磁共振灌注成像后处理系统，其特征在于，所述参数获取模块具 体包括： 分解单元，用于根据改进的基于非参数模型的奇异矩阵分解方法中的矩阵A进行加权 组合，并对矩阵A进行SVD分解得到对角矩阵W，正交矩阵V，三角矩阵U; 相对定量参数获取单元，用于根据阈值策略得到处理后的对角矩阵S，并根据对角矩阵 S确定灌注的相对定量参数;其中S = l/W。
【专利摘要】本发明提供了一种磁共振灌注成像后处理方法及系统，方法包括：对原始脑灌注图像滤波，去除脑部边缘，得到滤波后脑灌注图像；接收用户的区域选择指令，获取在滤波后脑灌注图像的中动脉区的选定区域，输出评估动脉输入函数的加权动脉输入曲线；根据改进的Gamma函数分别对加权动脉输入曲线和浓度时间曲线进行拟合，得到优化的动脉输入曲线和浓度时间曲线；根据改进的基于非参数模型的奇异矩阵分解方法求解灌注的相对定量参数。本发明通过对AIF和求解矩阵进行加权优化，降低对噪声的敏感度，采用简化且有效的Gamma函数进行拟合，将非线性问题转化为线性求解，加快了后处理速度。
【IPC分类】G06T7/00
【公开号】CN105701815
【申请号】CN201610019251
【发明人】刘素娟, 胡曙光
【申请人】深圳安科高技术股份有限公司
【公开日】2016年6月22日
【申请日】2016年1月12日