发明名称 一种提高二维阵列实时超声体成像质量的方法
摘要 本发明公开了一种提高二维阵列实时超声体成像质量的方法,包括以下步骤:设置面阵阵列的发射和接收阵列布局,对面阵阵列执行参数初始化,将超声波的扫描区域的极坐标形式转换为直角坐标形式,利用远场连续波对面阵阵列的发射超声声场进行近似,以获取面阵阵列的空间脉冲响应,对接收到的超声波信号进行CF处理,根据处理后的超声波信号获取面阵阵列的双向脉冲回波响应,对双向脉冲回波响应<img file="DDA00001960887700011.GIF" wi="126" he="47" />进行归一化对数压缩变换。本发明能够降低面阵的阵元数量和通道数,降低三维超声系统的前端电路复杂程度,并提高三维图像质量,达到实时体成像。
申请公布号 CN103565469A 申请公布日期 2014.02.12
申请号 CN201210270396.2 申请日期 2012.08.01
申请人 华中科技大学 发明人 尉迟明;丁明跃;贾艳平;许梦玲;李旭
分类号 A61B8/00(2006.01)I 主分类号 A61B8/00(2006.01)I
代理机构 华中科技大学专利中心 42201 代理人 朱仁玲
主权项 1.一种提高二维阵列实时超声体成像质量的方法,其特征在于,包括以下步骤:(1)设置面阵阵列的发射和接收阵列布局;(2)对面阵阵列执行参数初始化;(3)将超声波的扫描区域的极坐标形式<img file="FDA00001960887400011.GIF" wi="168" he="53" />转换为直角坐标形式,其中r<sub>0</sub>为扫描深度,θ为扫描仰角,<img file="FDA00001960887400012.GIF" wi="31" he="34" />为扫描方位角,0&lt;θ≤π/2,<img file="FDA00001960887400013.GIF" wi="220" he="44" />直角坐标的横、纵坐标u和v分别为:<img file="FDA00001960887400014.GIF" wi="293" he="46" /><img file="FDA00001960887400015.GIF" wi="258" he="46" />(4)利用远场连续波对面阵阵列的发射超声声场进行近似,以获取面阵阵列的空间脉冲响应,具体而言,包括以下子步骤:(4-1)根据扫描区域的聚集中心对面阵阵列中的各阵元进行聚焦延时,具体而言,聚焦中心为<img file="FDA00001960887400016.GIF" wi="222" he="53" />θ<sub>0</sub>为聚焦中心的仰角,<img file="FDA00001960887400017.GIF" wi="43" he="38" />为聚焦中心的方位角。各个阵元到聚焦中心的聚焦延时Δt<sub>nm</sub>为:<maths num="0001"><![CDATA[<math><mrow><mi>&Delta;</mi><msub><mi>t</mi><mi>nm</mi></msub><mo>=</mo><mfrac><msub><mi>r</mi><mn>0</mn></msub><mi>c</mi></mfrac><mo>-</mo><mfrac><msqrt><msup><mrow><mo>(</mo><msub><mi>r</mi><mn>0</mn></msub><msub><mi>u</mi><mn>0</mn></msub><mo>-</mo><mi>nd</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>r</mi><mn>0</mn></msub><msub><mi>v</mi><mn>0</mn></msub><mo>-</mo><mi>md</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>cos</mi><msub><mi>&theta;</mi><mn>0</mn></msub><msub><mi>r</mi><mn>0</mn></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mi>c</mi></mfrac></mrow></math>]]></maths>其中n=0..N-1,m=0...M-1,(n,m)用于表示面阵阵列中的各个阵元,N和M分别为面阵阵列中一行和一列的阵元数量,nd,md分别代表阵元(n,m)的坐标位置,<img file="FDA00001960887400019.GIF" wi="343" he="50" /><img file="FDA000019608874000110.GIF" wi="326" he="50" />(r<sub>0</sub>u<sub>0</sub>,r<sub>0</sub>v<sub>0</sub>)代表聚焦中心对直角坐标的投影;(4-2)获取面阵阵列中各阵元到达空间区域中各点的时间:<maths num="0002"><![CDATA[<math><mrow><msub><mi>t</mi><mi>nm</mi></msub><mo>=</mo><mfrac><msqrt><msup><mrow><mo>(</mo><mi>ru</mi><mo>-</mo><mi>nd</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>rv</mi><mo>-</mo><mi>md</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>cos</mi><mi>&theta;r</mi><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mi>c</mi></mfrac></mrow></math>]]></maths>(4-3)获取面阵阵列的发射空间脉冲响应<img file="FDA000019608874000112.GIF" wi="160" he="50" /><img file="FDA000019608874000113.GIF" wi="889" he="114" />其中<img file="FDA00001960887400021.GIF" wi="169" he="105" />代表发射超声波的波数,N<sub>1</sub>,M<sub>1</sub>为发射阵列中一行和一列的阵元数量,N<sub>1</sub>×M<sub>1</sub>为发射阵列中总的阵元数量;(4-4)获取面阵阵列的接收空间脉冲响应<img file="FDA00001960887400022.GIF" wi="260" he="50" /><img file="FDA00001960887400023.GIF" wi="938" he="114" />其中N<sub>2</sub>,M<sub>2</sub>为接收阵列中一行和一列的阵元数量,N<sub>2</sub>×M<sub>2</sub>为接收阵列中总的阵元数量;(5)对接收到的超声波信号进行CF处理,具体包括以下子步骤:(5-1)计算接收阵列中各个阵元所接收到的超声波信号所需的CF值:<maths num="0003"><![CDATA[<math><mrow><mi>CF</mi><mo>=</mo><mfrac><mrow><mo>|</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><msub><mi>N</mi><mn>2</mn></msub><mo>-</mo><mn>1</mn></mrow></munderover><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>0</mn></mrow><mrow><msub><mi>M</mi><mn>2</mn></msub><mo>-</mo><mn>1</mn></mrow></munderover><mi>S</mi><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mi>m</mi><mo>)</mo></mrow><mo>|</mo></mrow><mrow><msub><mi>N</mi><mn>2</mn></msub><mo>&times;</mo><msub><mi>M</mi><mn>2</mn></msub><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><msub><mi>N</mi><mn>2</mn></msub></munderover><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>0</mn></mrow><msub><mi>M</mi><mn>2</mn></msub></munderover><msup><mrow><mo>|</mo><mi>S</mi><mrow><mo>(</mo><mi>n</mi><mo>,</mo><mi>m</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup></mrow></mfrac></mrow></math>]]></maths>其中<img file="FDA00001960887400025.GIF" wi="422" he="128" />表示各个阵元接收到的超声波信号;(5-2)根据CF值对接收到的超声波信号进行CF处理,具体采用以下公式:<img file="FDA00001960887400026.GIF" wi="654" he="50" />(6)根据处理后的超声波信号获取面阵阵列的双向脉冲回波响应<img file="FDA00001960887400027.GIF" wi="150" he="47" /><img file="FDA00001960887400028.GIF" wi="809" he="50" />(7)对双向脉冲回波响应<img file="FDA00001960887400029.GIF" wi="126" he="47" />进行归一化对数压缩变换:<img file="FDA000019608874000210.GIF" wi="711" he="145" />
地址 430074 湖北省武汉市洪山区珞瑜路1037号