基于参数估计的压缩感知成像方法与流程

文档序号:12457027阅读:来源:国知局

技术特征:

1.基于参数估计的压缩感知成像方法,其特征在于:基于参数估计的压缩感知成像方法具体过程为:

步骤一、将每个距离单元的回波信号表示为线性调频信号模型;

步骤二、根据线性调频信号模型构造压缩感知稀疏基矩阵,根据压缩感知稀疏基矩阵构造图像Θ(γ0,α);

步骤三、依据图像对比度的定义选取检索范围,结合压缩感知的成像方法进行参数估计,得到每个距离单元的回波信号的初始参考频率和每个距离单元的回波信号频率的倾斜度两个参数的估计值;

步骤四、根据每个距离单元的回波信号的初始参考频率和每个距离单元的回波信号频率的倾斜度两个参数的估计值,结合压缩感知的成像方法进行成像,得到回波缺失情况下目标的高分辨图像。

2.根据权利要求1所述基于参数估计的压缩感知成像方法,其特征在于:所述步骤一中将每个距离单元的回波信号表示为线性调频信号模型;具体过程为:

线性调频信号模型的建立:

<mrow> <msub> <mi>s</mi> <mi>m</mi> </msub> <mrow> <mo>(</mo> <mi>T</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>&mu;</mi> <mi>m</mi> </msub> </munderover> <msub> <mi>&sigma;</mi> <mi>i</mi> </msub> <mi>exp</mi> <mo>&lsqb;</mo> <mi>j</mi> <mn>2</mn> <mi>&pi;</mi> <mrow> <mo>(</mo> <msub> <mi>f</mi> <msub> <mi>d</mi> <mi>i</mi> </msub> </msub> <mi>T</mi> <mo>+</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msub> <mi>&gamma;</mi> <mi>m</mi> </msub> <msup> <mi>T</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> <mo>&rsqb;</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>

其中,μm为散射点的个数,σi为第i个散射点的反射系数,i取值为正整数;为第i个散射点的多普勒频率,j是虚数单位,j2=-1;T为时间;γm为信号线性调频参数;sm(T)为线性调频信号模型;

γm的具体表示形式如下

<mrow> <msub> <mi>&gamma;</mi> <mi>m</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mn>2</mn> <msubsup> <mi>y&omega;</mi> <mn>0</mn> <mn>2</mn> </msubsup> </mrow> <mi>&lambda;</mi> </mfrac> <mo>=</mo> <mfrac> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mn>0</mn> </msub> <mo>+</mo> <mi>N</mi> <mi>&Delta;</mi> <mi>r</mi> <mo>)</mo> </mrow> <msubsup> <mi>&omega;</mi> <mn>0</mn> <mn>2</mn> </msubsup> </mrow> <mi>&lambda;</mi> </mfrac> <mo>=</mo> <msub> <mi>&gamma;</mi> <mn>0</mn> </msub> <mo>+</mo> <mi>N</mi> <mi>&alpha;</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>

其中,△r为两个连续距离单元间的距离,即为距离向分辨率;y为目标位置;为角速度;λ为波长;y0为目标初始位置;N为距离单元个数,取值为正整数;γ0为每个距离单元的回波信号的初始参考频率,α为每个距离单元的回波信号频率的倾斜度。

3.根据权利要求2所述基于参数估计的压缩感知成像方法,其特征在于:所述步骤二中根据线性调频信号模型构造压缩感知的稀疏基矩阵,根据压缩感知的稀疏基矩阵构造图像Θ(γ0,α);具体过程为:

依据线性调频信号模型,设K为多普勒频率的离散度,构造压缩感知的稀疏基矩阵如下:

其中,K为多普勒频率的离散度,fk为离散多普勒频率,fk=k·fr/K,k=1,2,...,K,K取值为正整数;

则每个距离单元的回波表示为:

Sm=Φ(γmm (4)

式中,Sm=[sm(1),sm(2),...,sm(N)]T,θm为包含了所有散射点的反射系数;为Φ(γm)中每一列;图像Θ(γ0,α)为所有距离单元θm的集合,即Θ(γ0,α)={θ12,...,θM}。

4.根据权利要求3所述基于参数估计的压缩感知成像方法,其特征在于:所述步骤三中依据图像对比度的定义选取检索范围,结合压缩感知的成像方法进行参数估计,得到每个距离单元的回波信号的初始参考频率和每个距离单元的回波信号频率的倾斜度两个参数的估计值;

依据图像对比度的定义选取检索范围,结合压缩感知的成像方法进行参数估计,得到参数;具体过程为:

图像对比度的定义,具体如下

<mrow> <msub> <mi>C</mi> <mrow> <mi>&Theta;</mi> <mrow> <mo>(</mo> <msub> <mi>&gamma;</mi> <mn>0</mn> </msub> <mo>,</mo> <mi>&alpha;</mi> <mo>)</mo> </mrow> </mrow> </msub> <mo>=</mo> <mo>&lt;</mo> <mi>&Theta;</mi> <msup> <mrow> <mo>(</mo> <msub> <mi>&gamma;</mi> <mn>0</mn> </msub> <mo>,</mo> <mi>&alpha;</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>&gt;</mo> <mo>/</mo> <mo>&lt;</mo> <mi>&Theta;</mi> <mrow> <mo>(</mo> <msub> <mi>&gamma;</mi> <mn>0</mn> </msub> <mo>,</mo> <mi>&alpha;</mi> <mo>)</mo> </mrow> <msup> <mo>&gt;</mo> <mn>2</mn> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>

其中,<·>表示平均算子;

故在参数检索范围内,当图像对比度达到最大值时表明所选择的参数最合适,即

<mrow> <mo>{</mo> <msub> <mi>&gamma;</mi> <mn>0</mn> </msub> <mo>,</mo> <mi>&alpha;</mi> <mo>}</mo> <mo>=</mo> <mi>arg</mi> <mi> </mi> <msub> <mi>maxC</mi> <mrow> <mi>&Theta;</mi> <mrow> <mo>(</mo> <msub> <mi>&gamma;</mi> <mn>0</mn> </msub> <mo>,</mo> <mi>&alpha;</mi> <mo>)</mo> </mrow> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> <mo>.</mo> </mrow>

5.根据权利要求4所述基于参数估计的压缩感知成像方法,其特征在于:所述α=0.65,γ0=-5.9。

当前第2页1 2 3 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1