一种磁共振成像匀场方法

摘要

本发明提出一种磁共振成像匀场方法，根据自适应Ostu阈值将成像区域分为信噪比较高的可靠区域(R_reliable区域)和信噪比较低的非可靠区域（R_fg区域），通过插值计算根据可靠区域的像素点对非可靠区域的像素点进行恢复，计算出包含可靠区域和非可靠区域内所有像素点的缠绕相位图，排除低信噪比区域噪声对匀场计算过程的影响；搜索相位变化最缓慢的点为起始点进行相位解缠，提高了相位解缠和匀场计算的精度，从而提高了匀场的可靠性。

主权项

一种磁共振成像匀场方法，包括以下步骤：选定扫描序列对被成像组织扫描，获得磁共振信号，进行图像重建；根据重建的图像计算Harmonic系数；使用Harmonic系数控制匀场线圈对磁场进行补偿；其特征在于，所述Harmonic系数的计算方法包括以下步骤：a)根据阈值将图像分为可靠区域、非可靠区域和背景区域；b)根据可靠区域内的像素点对非可靠区域内的像素点进行插值计算，得到包括可靠区域和非可靠区域内所有像素点的缠绕相位图；c)以解缠种子点作为相位解缠的起始点，对所述缠绕相位图作相位解缠，根据解缠后的相位图得到以解缠种子点为参考点的Harmonic系数；所述解缠种子点通过以下方式获取：计算缠绕相位图中每个像素点与其邻域内其他各像素点的相位变化，以所述相位变化最小的像素点作为解缠种子点；d)将以解缠种子点为参考点的Harmonic系数转换为以磁体中心为参考点的Harmonic系数，所述将以解缠种子点为参考点的Harmonic系数转化为以磁体中心为参考点的Harmonic系数的可以通过以下公式获取：$\begin{array}{c}{A}_{10}={A^{\prime }}_{10}-{2z}_s{A^{\prime }}_{20}-\frac{\sqrt{15}}{2}{x}_s{A^{\prime }}_{21}-\frac{\sqrt{15}}{2}{y}_s{B^{\prime }}_{21}\\ A_{11}={A^{\prime }}_{11}+{A^{\prime }}_{20}{x}_s-\sqrt{5}{z}_s{A^{\prime }}_{21}-\sqrt{5}{x}_s{A^{\prime }}_{22}-\sqrt{5}{y}_s{B^{\prime }}_{22}\\ B_{11}={B^{\prime }}_{11}+{A^{\prime }}_{20}{y}_s-\sqrt{5}{z}_s{B^{\prime }}_{21}-\sqrt{5}{x}_s{A^{\prime }}_{22}-\sqrt{5}{y}_s{B^{\prime }}_{22}\\ A_{20}={A^{\prime }}_{20}\\ A_{21}={A^{\prime }}_{21}\\ A_{22}={A^{\prime }}_{22}\\ B_{21}={B^{\prime }}_{21}\\ B_{22}={B^{\prime }}_{22}\end{array}$式中：(x_s,y_s,z_s)为所述解缠种子点，以解缠种子点为参考点的Harmonic系数为A'₁₀、A'₁₁、B'₁₁、A'₂₀、A'₂₁、A'₂₂、B'₂₁、B'₂₂；以磁场中心为参考点的Harmonic系数为A₁₀、A₁₁、B₁₁、A₂₀、A₂₁、A₂₂、B₂₁、B₂₂。