1.一种基于超声成像系统点扩散函数测量和压缩感知的超声成像方法,其包括以下步骤:获得所述超声成像系统的点扩散函数;
根据所述点扩散函数建立关系矩阵M;以及
建立反映超声回波射频信号s与超声图像的像素点上的散射子散射强度分布I之间关系的方程组:s=MI;
其特征在于,
设置一个阈值,将所述关系矩阵M中低于所述阈值的元素都置为0。
2.根据权利要求1所述的方法,其特征在于,通过实验测量的方式获得所述超声成像系统在某种超声发射模式下在与超声图像上所有像素点对应的空间位置处的所述点扩散函数。
3.根据权利要求1或2所述的方法,其特征在于,所述实验测量的方式包括以下步骤:(i)在与超声图像上某个像素点对应的空间位置上放置尺寸足够小的散射子,同时在其他位置上没有任何散射子的干扰;
(ii)所述超声成像系统发射某种模式的超声信号,经所述散射子反射后,形成超声回波信号,所述超声回波信号被所述超声成像系统接收,其中被接收的所述超声回波信号即所述系统在所述某种超声发射模式下在所述空间位置的点扩散函数;
(iii)对所述超声图像上的所有像素点执行步骤(i)和(ii),得到所述超声成像系统在某种超声发射模式下在与所述超声图像上所有像素点对应的空间位置处的所述点扩散函数。
4.根据权利要求2所述的方法,其特征在于,所述通过实验测量的方式获得所述超声成像系统在某种超声发射模式下在与超声图像上所有像素点对应的空间位置处的所述点扩散函数包括但不限于:在水槽中利用尼龙线测量的方式。
5.根据权利要求1所述的方法,其特征在于所述方法还包括以下步骤:通过成熟的压缩感知算法对所述方程组s=MI进行求解,得到向量I,将其从向量变换为对应图像像素点数量的矩阵,再经过调整动态范围和数字扫描变换步骤,获得希望得到的超声图像。
6.根据权利要求5所述的方法,其特征在于,所述成熟的压缩感知算法包括:匹配追踪法、Bregman算法、operator/variable splitting、FPC算法、L1-magic算法、牛顿下降法。
7.根据权利要求5或6所述的方法,其特征在于,当I不是稀疏的时,对I进行稀疏变换Ψ,令θ=ΨI,其中θ是I在稀疏变换域内的系数,θ是稀疏的,根据压缩感知的求解公式来求解I,其中β表示允许存在多少噪声成分。
8.根据权利要求7所述的方法,其特征在于,所述稀疏变换Ψ包括:离散余弦变换(DCT)、各种小波变换。
9.根据权利要求1或2所述的方法,其特征在于,所述超声成像系统采用超声平面波发射模式、超声凸面波发射模式、超声凹面波发射模式、或任意不规则波形信号激励的发射模式。