发明名称 骨表微血管血流灌注分离检测与成像方法
摘要 骨表微血管血流灌注分离检测与成像方法,先编程满足骨表多模式漏Lamb波条件下发射与接收波束,再构建骨表多模式漏Lamb波微泡母小波,然后采用骨表多模式漏Lamb波微泡小波方法对入射、反射、多模式漏Lamb波混叠射频回波信号进行模式分离检测,最后实现对骨表多模式漏Lamb波微泡小波变换微血管造影成像与血流灌注参量成像,本发明可有效提高分辨率、造影组织比,抑制了骨表多模式漏Lamb波对骨表组织微血管造影成像的混叠干扰,这对于恶性骨肿瘤早期诊断与评估具有重要临床意义。
申请公布号 CN103381096B 申请公布日期 2015.04.15
申请号 CN201310139512.1 申请日期 2013.04.19
申请人 西安交通大学 发明人 万明;王弟亚;翟宇;钟徽
分类号 A61B8/00(2006.01)I 主分类号 A61B8/00(2006.01)I
代理机构 西安智大知识产权代理事务所 61215 代理人 弋才富
主权项 骨表微血管血流灌注分离检测与成像方法,其特征在于,包括以下步骤:步骤一、在可编程控制超声成像系统上,依据肌骨系统具体目标和临床要求,编程控制窄/宽波束发射接收模式、同相位/脉冲逆转相位设置模式;步骤二、在主控计算机平台上,基于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>&omega;</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>&omega;</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>&omega;</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>&omega;</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>&omega;</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>&omega;</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>&omega;</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>&omega;</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>&omega;</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>&omega;</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>&rho;R</mi><mover><mi>R</mi><mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mrow></mover><mo>+</mo><mfrac><mn>3</mn><mn>2</mn></mfrac><mi>&rho;</mi><msup><mover><mi>R</mi><mo>&CenterDot;</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>&sigma;</mi></mrow><msub><mi>R</mi><mn>0</mn></msub></mfrac><mo>+</mo><mfrac><mrow><mn>2</mn><mi>&chi;</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>&gamma;</mi></mrow></msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mfrac><mrow><mn>3</mn><mi>&gamma;</mi></mrow><mi>c</mi></mfrac><mover><mi>R</mi><mo>&CenterDot;</mo></mover><mo>)</mo></mrow><mo>-</mo><mfrac><mrow><mn>4</mn><mi>&mu;</mi><mover><mi>R</mi><mo>&CenterDot;</mo></mover></mrow><mi>R</mi></mfrac><mo>-</mo><mfrac><mrow><mn>2</mn><mi>&sigma;</mi></mrow><mi>R</mi></mfrac><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mfrac><mover><mi>R</mi><mo>&CenterDot;</mo></mover><mi>c</mi></mfrac><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>-</mo><mfrac><mrow><mn>2</mn><mi>&chi;</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>&CenterDot;</mo></mover><mo>)</mo></mrow><mo>-</mo><mn>12</mn><msub><mi>&mu;</mi><mi>sh</mi></msub><mi>&epsiv;</mi><mfrac><mover><mi>R</mi><mo>&CenterDot;</mo></mover><mrow><mi>R</mi><mrow><mo>(</mo><mi>R</mi><mo>-</mo><mi>&epsiv;</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>&rho;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>&CenterDot;</mo><mo>&CenterDot;</mo></mrow></mover><mo>+</mo><mi>R</mi><msup><mover><mi>R</mi><mo>&CenterDot;</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的造影图像。
地址 710049 陕西省西安市咸宁路28号