主权项 |
一种磁共振参数重建方法,包括以下步骤:按随机降采矩阵采集得到弥散加权成像的K空间数据;建立所述弥散加权成像的数学模型,所述数学模型含有弥散张量系数;定义目标函数,采用最优化问题求解的方法拟合所述弥散张量系数,并加入稀疏约束,所述目标函数加入稀疏约束定义如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>Φ</mi><mrow><mo>(</mo><mi>D</mi><mo>)</mo></mrow><mo>=</mo><mo>{</mo><munder><mi>Σ</mi><mi>j</mi></munder><msubsup><mrow><mo>|</mo><mo>|</mo><msub><mi>d</mi><mi>j</mi></msub><mo>-</mo><msub><mi>P</mi><mi>j</mi></msub><msub><mi>Ff</mi><mi>j</mi></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><mi>λ</mi><msub><mrow><mo>|</mo><mo>|</mo><mi>C</mi><mo>|</mo><mo>|</mo></mrow><mn>1,2</mn></msub><mo>}</mo><mo>,</mo><mi>C</mi><mo>=</mo><mi>Ψ</mi><mo>[</mo><msub><mi>f</mi><mn>1</mn></msub><mo>,</mo><msub><mi>f</mi><mn>2</mn></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>f</mi><mi>J</mi></msub><mo>]</mo></mrow>]]></math><img file="FDA0000594835470000011.GIF" wi="930" he="137" /></maths>其中,d<sub>j</sub>为降采的K空间数据,P为降采矩阵,F是傅立叶变换,C为稀疏系数矩阵,所述稀疏系数矩阵为图像像素×弥散方向个数大小的矩阵;稀疏约束||C||<sub>1,2</sub>为求取稀疏系数矩阵C的L1‑L2模,Ψ为稀疏变换,所述稀疏变换为对所述弥散加权图像的有限差分;根据降采数据,求解所述弥散张量系数,使得所述目标函数最小。 |