[发明专利]一种优化的磁共振灌注成像图像中血容量量化方法

权利要求书

1.一种优化的磁共振灌注成像图像中血容量量化方法，其特征在于步骤如下：

步骤1：根据C(x,y,z,t)=-1TEln(S(x,y,z,t)S0(x,y,z)),]]>将图像中三维像素的时间-灰度信号曲线转换为时间-造影剂浓度曲线，其中：(x，y，z)为三维像素在图像矩阵中的索引值，C(x，y，z，t)为索引值(x，y，z)处三维像素的时间-造影剂浓度曲线，TE为磁共振灌注成像图像的回波时间，S(x，y，z，t)为索引值(x，y，z)处三维像素的时间-灰度信号曲线，S₀(x，y，z)为造影剂到达之前索引值(x，y，z)处三维像素的基准灰度值，ln(^·)为对数运算符；所述的基准灰度值为该像素时间-灰度信号曲线中造影剂达到之前的时间点的均值；

步骤2：根据造影剂渗漏效应矫正算法，矫正三维像素的时间-造影剂浓度曲线中因造影剂渗漏效应带来的附加顺磁效应，得到矫正后的时间-造影剂浓度曲线C_corr(t)；

步骤3：在图像中脑主中央动脉邻域选取一个像素，将该像素的时间-造影剂浓度曲线作为动脉输入函数f(t)；所述选取像素的时间-造影剂浓度曲线具有典型双峰形状；

步骤4：采用标准奇异值分解算法，以步骤3得到的f(t)为输入函数，对步骤2得到的C_corr(t)做去卷积运算，得到残余函数R(t)；

步骤5：采用gamma函数以非线性最小二乘算法拟合步骤3中得到的动脉输入函数f(t)，得到f(t)的首次通过f₁(t)；

步骤6：将步骤5中得到的f₁(t)与步骤4中得到的R(t)做卷积运算以重建步骤2中的C_corr(t)的首次通过分量C₁(t)；

步骤7：根据rCBV∝∫t0tC1(τ)dτ]]>，利用步骤6中得到的C₁(t)来量化该像素点处的局部血容量rCBV，其中rCBV为局部血容量，t₀为造影剂到达时刻，τ为积分常数。

2.根据权利要求1所述的优化的磁共振灌注成像图像中血容量量化方法，其特征在于：步骤2中的造影剂泄漏效应校正算法采用文献“Boxerman，J.L.，K.M.Schmainda，and R.M.Weisskoff，Relative cerebral blood volume mapscorrected for contrast agent extravasation significantly correlate withglioma tumor grade，whereas uncorrected maps do not.AJNR Am J Neuroradiol，2006.27(4)：p.859-67.”中公布的算法。