[發明專利]基于壓縮感知的超聲成像方法和終端

在審

• 申請號: 201611091236.6
• 申請日: 2016-12-01
• 公開（公告）號: CN106725594A
• 公開（公告）日: 2017-05-31
• 發明（設計）人: 張瓊
• 申請（專利權）人: 汕頭大學醫學院
• 主分類號: A61B8/00
• 分類號: A61B8/00
• 代理公司: 北京路浩知識產權代理有限公司11002
• 代理人: 湯財寶

1.一種基于壓縮感知的超聲成像方法，其特征在于，包括以下步驟：

將超聲陣元接收的信號與預設的線性測量矩陣進行內積計算，得到實際接收的后向散射信號y；

根據分辨率要求，離散化目標區域，生成對應的超完備字典矩陣V_DCS；

結合所述超完備字典矩陣V_DCS和所述后向散射信號y，通過如下公式進行最大似然估計的求解，得到目標向量f：

y＝V_DCS·f；

將所述目標向量f轉換為圖像矩陣。

2.根據權利要求1所述的基于壓縮感知的超聲成像方法，其特征在于，所述將超聲陣元接收的信號與預設的線性測量矩陣進行內積計算，得到實際接收的后向散射信號y，包括：

第i個超聲陣元在時間t接收到信號x_i(t)：

$x_{i\left(t\right)}={\mathrm{\Σ}}_{j=1}^Q{s}_j(t-{\mathrm{\τ}}_i({\mathrm{\ρ}}_j\left)\right),i=1,2,...K---\left(1\right);$

其中，

Q表示散射點的數目，K表示超聲陣元的數目，s_j表示第j個散射點的后向散射信號，τ_i(ρ_j)表示散射點ρ_j到第i個超聲陣元的傳輸延遲，表示第i個超聲陣元的空間位置，表示散射點ρ_j的空間位置，c表示聲速；

將第i個超聲陣元接收到的信號x_i表示為：

x_i＝ψ*θ_i (3)；

其中，‖θ_i‖₀＝S；

S為信號x_i(t)在基ψ中的稀疏度；

將x_i與大小為M_i×N_t的預設的線性測量矩陣進行內積計算，得到第i個超聲陣元實際接收的后向散射信號y_i：

其中M_i＜＜N_t，M_i為第i個超聲陣元實際接收的后向散射信號y_i的個數，且M_i在其中一個陣元中的稀疏度為2S，在其他陣元的稀疏度為S；

則K個超聲陣元實際接收的后向散射信號y表示為：

y＝[y₁^T,y₂^T,…y_k^T]^T (5)。