主权项 |
骨表微血管血流灌注分离检测与成像方法,其特征在于,包括以下步骤:步骤一、在可编程控制超声成像系统上,依据肌骨系统具体目标和临床要求,编程控制窄/宽波束发射接收模式、同相位/脉冲逆转相位设置模式;步骤二、在主控计算机平台上,基于Rayleigh‑Lamb频散方程即公式(1)、(2)与Morgen模型修正的Herring–Trilling微泡振动方程即公式(3)、(4),构建骨表多模式漏Lamb波微泡母小波;<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><mi>tan</mi><mrow><mo>(</mo><mfrac><mi>d</mi><mn>2</mn></mfrac><msqrt><msup><mrow><mo>(</mo><mi>ω</mi><mo>/</mo><msub><mi>C</mi><mi>l</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><msubsup><mi>k</mi><mi>z</mi><mn>2</mn></msubsup></msqrt><mo>)</mo></mrow><mo>/</mo><mi>tan</mi><mrow><mo>(</mo><mfrac><mi>d</mi><mn>2</mn></mfrac><msqrt><msup><mrow><mo>(</mo><mi>ω</mi><mo>/</mo><msub><mi>C</mi><mi>s</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><msubsup><mi>k</mi><mi>z</mi><mn>2</mn></msubsup></msqrt><mo>)</mo></mrow><mo>=</mo></mtd></mtr><mtr><mtd><mo>-</mo><mn>4</mn><msubsup><mi>k</mi><mi>z</mi><mn>2</mn></msubsup><msqrt><mrow><mo>(</mo><msup><mrow><mo>(</mo><mi>ω</mi><mo>/</mo><msub><mi>C</mi><mi>l</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><msubsup><mi>k</mi><mi>z</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mrow><mo>(</mo><msup><mrow><mo>(</mo><mi>ω</mi><mo>/</mo><msub><mi>C</mi><mi>s</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><msubsup><mi>k</mi><mi>z</mi><mn>2</mn></msubsup><mo>)</mo></mrow></msqrt><mo>/</mo><msup><mrow><mo>(</mo><msubsup><mi>k</mi><mi>z</mi><mn>2</mn></msubsup><mo>-</mo><msup><mrow><mo>(</mo><mi>ω</mi><mo>/</mo><msub><mi>C</mi><mi>s</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msubsup><mi>k</mi><mi>z</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mn>2</mn></msup></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000632220300000011.GIF" wi="1317" he="279" /></maths> (1,S<sub>P</sub>模式)<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><mi>tan</mi><mrow><mo>(</mo><mfrac><mi>d</mi><mn>2</mn></mfrac><msqrt><msup><mrow><mo>(</mo><mi>ω</mi><mo>/</mo><msub><mi>C</mi><mi>l</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><msubsup><mi>k</mi><mi>z</mi><mn>2</mn></msubsup></msqrt><mo>)</mo></mrow><mo>/</mo><mi>tan</mi><mrow><mo>(</mo><mfrac><mi>d</mi><mn>2</mn></mfrac><msqrt><msup><mrow><mo>(</mo><mi>ω</mi><mo>/</mo><msub><mi>C</mi><mi>s</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><msubsup><mi>k</mi><mi>z</mi><mn>2</mn></msubsup></msqrt><mo>)</mo></mrow><mo>=</mo></mtd></mtr><mtr><mtd><mo>-</mo><msup><mrow><mo>(</mo><msubsup><mi>k</mi><mi>z</mi><mn>2</mn></msubsup><mo>-</mo><msup><mrow><mo>(</mo><mi>ω</mi><mo>/</mo><msub><mi>C</mi><mi>s</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msubsup><mi>k</mi><mi>z</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>/</mo><mn>4</mn><msubsup><mi>k</mi><mi>z</mi><mn>2</mn></msubsup><msqrt><mrow><mo>(</mo><msup><mrow><mo>(</mo><mi>ω</mi><mo>/</mo><msub><mi>C</mi><mi>l</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><msubsup><mi>k</mi><mi>z</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mrow><mo>(</mo><msup><mrow><mo>(</mo><mi>ω</mi><mo>/</mo><msub><mi>C</mi><mi>s</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><msubsup><mi>k</mi><mi>z</mi><mn>2</mn></msubsup><mo>)</mo></mrow></msqrt></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000632220300000012.GIF" wi="1295" he="277" /></maths> (2,A<sub>P</sub>模式)其中,ω为角频率;C<sub>l</sub>、C<sub>s</sub>分别为介质中纵波声速和横波声速;k<sub>z</sub>为沿波导轴(即z轴)的波数,k<sub>z</sub>=2π/λ;d为板厚;<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mfenced open='' close=''><mtable><mtr><mtd><mi>ρR</mi><mover><mi>R</mi><mrow><mo>·</mo><mo>·</mo></mrow></mover><mo>+</mo><mfrac><mn>3</mn><mn>2</mn></mfrac><mi>ρ</mi><msup><mover><mi>R</mi><mo>·</mo></mover><mn>2</mn></msup><mo>=</mo><mrow><mo>(</mo><msub><mi>P</mi><mn>0</mn></msub><mo>+</mo><mfrac><mrow><mn>2</mn><mi>σ</mi></mrow><msub><mi>R</mi><mn>0</mn></msub></mfrac><mo>+</mo><mfrac><mrow><mn>2</mn><mi>χ</mi></mrow><msub><mi>R</mi><mn>0</mn></msub></mfrac><mo>)</mo></mrow><msup><mrow><mo>(</mo><mfrac><msub><mi>R</mi><mn>0</mn></msub><mi>R</mi></mfrac><mo>)</mo></mrow><mrow><mn>3</mn><mi>γ</mi></mrow></msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mfrac><mrow><mn>3</mn><mi>γ</mi></mrow><mi>c</mi></mfrac><mover><mi>R</mi><mo>·</mo></mover><mo>)</mo></mrow><mo>-</mo><mfrac><mrow><mn>4</mn><mi>μ</mi><mover><mi>R</mi><mo>·</mo></mover></mrow><mi>R</mi></mfrac><mo>-</mo><mfrac><mrow><mn>2</mn><mi>σ</mi></mrow><mi>R</mi></mfrac><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mfrac><mover><mi>R</mi><mo>·</mo></mover><mi>c</mi></mfrac><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>-</mo><mfrac><mrow><mn>2</mn><mi>χ</mi></mrow><mi>R</mi></mfrac><msup><mrow><mo>(</mo><mfrac><msub><mi>R</mi><mn>0</mn></msub><mi>R</mi></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mfrac><mn>3</mn><mi>c</mi></mfrac><mover><mi>R</mi><mo>·</mo></mover><mo>)</mo></mrow><mo>-</mo><mn>12</mn><msub><mi>μ</mi><mi>sh</mi></msub><mi>ϵ</mi><mfrac><mover><mi>R</mi><mo>·</mo></mover><mrow><mi>R</mi><mrow><mo>(</mo><mi>R</mi><mo>-</mo><mi>ϵ</mi><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mrow><mo>(</mo><msub><mi>P</mi><mn>0</mn></msub><mo>+</mo><msub><mi>P</mi><mi>driv</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000632220300000013.GIF" wi="1597" he="348" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>P</mi><mo>=</mo><msup><mi>ρr</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><msup><mi>R</mi><mn>2</mn></msup><mover><mi>R</mi><mrow><mo>·</mo><mo>·</mo></mrow></mover><mo>+</mo><mi>R</mi><msup><mover><mi>R</mi><mo>·</mo></mover><mn>2</mn></msup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000632220300000014.GIF" wi="1596" he="112" /></maths>式中,ρ为液体密度;R为瞬态微泡半径;<img file="FDA0000632220300000015.GIF" wi="52" he="76" />为加速度;<img file="FDA0000632220300000016.GIF" wi="56" he="75" />为速度;P<sub>0</sub>为液体静态压;<img file="FDA0000632220300000021.GIF" wi="90" he="134" />为界面张力;R<sub>0</sub>为微泡的静态半径;χ为薄包膜的弹性模量;γ为气体的多归系数;c为液体中的声速;<img file="FDA0000632220300000022.GIF" wi="130" he="144" />为液体的黏滞作用;<img file="FDA0000632220300000023.GIF" wi="328" he="153" />近似为膜的黏滞阻尼,ε、μ<sub>sh</sub>分别为薄包膜的厚度及材料黏滞系数;P<sub>driv</sub>(t)为驱动声压,r为液体与微泡的距离;步骤三、对于采用团式灌注模式进行造影剂灌注后的肌骨系统目标,依次对各模式混叠作用造影微泡后的射频回波信号进行骨表多模式漏Lamb波微泡小波分析、信号重建,进行模式分离与检测;步骤四、对分离与检测的信号进行微血管造影成像处理,且依据所定义的灌注参量,进行骨表组织微血管血流灌注参量成像,获取高造影组织比CTR、造影骨导比CBR、信噪比SNR的造影图像。 |