[发明专利]一种PROPELLER技术中参考位置图像的非迭代生成方法

权利要求书

1.一种PROPELLER技术中参考位置图像的非迭代生成方法，其特征在于包括以下步骤：

步骤01：采集数据，采集N个数据带A，所述N个数据带A为旋转的同中心数据带，所述N个数据带A在K空间的中心区域重叠，N为正整数；

步骤02：处理数据，将步骤01采集到的N个数据带A先进行相位矫正，再进行网格化处理，处理后得到N个数据带B；

步骤03：转换数据，将N个数据带B在K空间中心重叠区域的二维数据转换为一维数据，所述一维数据为一列；每个数据带B含有D个数据，N个数据带B转换为N个D*1的列向量，转换后得到N个数据带M，D为正整数；

步骤04：构建矩阵，构建矫正矩阵C，其大小为D*N，即D行N列；

将步骤03中所述的N个D*1列向量形式的数据带M分别放置于矫正矩阵C的对应列，将第1个数据带M放置在矫正矩阵C的第1列，第i个数据带M放置在矫正矩阵C的第i列，i为正整数；

步骤05：奇异值分解，对矫正矩阵C进行奇异值分解，分解过程采用以下方程：

C＝UΣV^H (Ⅰ)

在式(Ⅰ)中，C为D*N大小的矫正矩阵，U为D*D大小的酉矩阵，V为N*N大小的酉矩阵，V^H为V的共轭转置，Σ为D*N大小的对角阵，对角阵中的每个元素均为矩阵C的奇异值；

步骤06：合成数据，将所述对角阵Σ中的奇异值元素按照降序排列，选择对应于最大奇异值的右奇异向量R中的对应列R_max，所述R_max为N*1大小的列向量；R_max为权重因子，使用R_max对C中的N列数据进行合成，采用以下方程：

D_cop＝C*R_max (Ⅱ)

在式(Ⅱ)中，C为D*N大小的矫正矩阵，R_max为N*1大小列向量形式的权重因子，D_cop为D*1大小形式的合成数据；

步骤07：图像重建，将步骤06合成的一维数据D_cop转换为二维数据，将所述一维数据D_cop转换后的二维数据执行傅立叶变换，得到的图像为基于奇异值分解生成的自适应参考位置图像；

步骤02中所述的相位矫正采用金字塔窗函数滤波的方法；步骤02中所述的网格化处理为将非笛卡尔坐标系下的数据转换到笛卡尔坐标系下的数据，所述网格化处理采用Jackson网格化处理算法；

步骤03中所述二维数据转换为一维数据采用列优先原则或者行优先原则，所述N个数据带B均采用相同的转换原则。