[发明专利]磁共振弥散加权成像自适应修正方法

权利要求书

1.磁共振弥散加权成像自适应修正方法，其特征在于：包括以下步骤：

步骤1，重复采集弥散加权图像N次，N为自然数，N≥3；

步骤2，基于原始图像或者压缩后的图像逐点构造相关性矩阵；步骤2包括以下步骤：

步骤2.1，对于第n次采集的图像中任意一个像素点(x,y)，取周围相邻的K个点，构成一个邻域向量Xn；

步骤2.2，对于N次重复采集的图像，每个像素点对应了N个邻域向量，按公式(1)计算其中第n个向量Xn和第m个向量Xm之间的相关性；

公式(1)中，x_i为向量Xn中的第i个元素，y_i为向量Xm中的第i个元素，为向量Xn的均值，为向量Xm的均值；

步骤2.3，任一像素点(x,y)对应一个N*N的相关性矩阵R(x,y)；

其中，r_1,1…r_1,N为按公式(1)计算得到的两两向量之间的相关性系数；

步骤3，主成分分析；获得每个相关性矩阵的最大特征值对应的特征向量；步骤3包括以下步骤：

步骤3.1，计算矩阵R(x,y)的特征值，找出最大的特征值；

步骤3.2，计算矩阵R(x,y)的最大特征值对应的特征向量γ；

步骤4，根据特征向量计算权重；步骤4中通过公式(2)计算权重：

公式(2)中，γ_n是特征向量γ的第n个元素，γ_min是特征向量γ的最小元素，γ_max是特征向量γ的最大元素,a和p为参数控制因子；其中，a＝0.2，p＝1；

步骤5，根据步骤4获得的权重对步骤1中采集到的原始图像进行加权合成，获得修正后的弥散加权图像。

2.如权利要求1所述的磁共振弥散加权成像自适应修正方法，其特征在于：步骤2之前利用插值算法对所有采集的原始图像进行压缩。

3.如权利要求1所述的磁共振弥散加权成像自适应修正方法，其特征在于：步骤3之前先对相关性矩阵做平滑滤波处理。

4.如权利要求3所述的磁共振弥散加权成像自适应修正方法，其特征在于：平滑滤波处理包括以下步骤；

步骤a，从每个像素点(x,y)对应的相关性矩阵R(x,y)中，取出第i个相关性系数，构成一个和图像矩阵大小相同的矩阵Ri；

步骤b，对矩阵Ri进行二维低通滤波；

步骤c，将滤波后的结果取代R(x,y)中对应的元素；

步骤d，重复a-c，直到R(x,y)中所有元素均处理完毕。

5.如权利要求1所述的磁共振弥散加权成像自适应修正方法，其特征在于：步骤5中通过公式(3)对原始图像进行加权合成；

公式(3)中，M_n为第n次采集得到弥散加权原始图像，w_n为权重。

6.如权利要求1所述的磁共振弥散加权成像自适应修正方法，其特征在于：步骤1中以相同的扫描参数重复采集弥散加权图像N次。