采用自适应处理的地下目标二维高分辨成像方法及系统

本发明属于雷达高分辨成像，尤其涉及一种采用自适应处理的地下目标二维高分辨成像方法。

背景技术：

1、地下成像技术在地球物理探测应用中具有重要作用。尽管合成孔径雷达(synthetic aperture radar,sar)在遥感成像应用中得到广泛认可，但它对射频干扰非常敏感，而且地下成像过程中土壤和空气之间的介电常数不同，会使信号的传播路径发生弯曲。此外，地面多径反射会带来非平稳干扰导致成像性能下降。rosenberg等人验证了非平稳干扰可以通过快时间空时自适应处理(fast-time space time adaptive processing,stap)来抑制，而不会造成明显的成像性能损失。

2、在sar成像技术中，高距离分辨率可采用宽带波形来实现，而方位向分辨率则受限于天线的有效孔径，这意味着要获得高方位角分辨率，就必须使用大孔径天线，但在地下成像中难以实现，这是因为地下sar成像中应采用低于1ghz的相对较低载频，否则由于土壤的低通电磁传播效应会造成严重的信号衰减。而且，不同于机载和星载sar，地下sar成像时土壤和空气之间不同的介电常数会使电磁波在其界面处的传播路径发生弯曲。更为严重的是，复杂地下介质中的多径传播效应会降低地下成像性能。

3、频控阵可看作是一种综合空间复用和频率复用的复合技术，与传统相控阵和mimo不同，频控阵在其阵元间采用小频偏来产生与距离有关的波束方向图。cetinepe等人研究表明，在多径环境下，如果适当设置频控阵的带宽，可以在期望的范围内形成准平台干涉效应。jingwei xu等人研究表明，前视频控阵sar可以缓解成像中的多普勒与距离关系。因此，频控阵的这些特性使频控阵sar在地下高分辨成像中具有应用优势。

4、通过上述分析，现有技术存在的问题及缺陷为：复杂地下介质中的多径传播效应会降低地下成像性能。

技术实现思路

1、针对现有技术存在的问题，本发明将快时间stap和频控阵相结合，用于sar地下高分辨成像，主要采用快时间stap抑制成像中的不良干扰，并采用频控阵距离依赖性波束方向图减少多径干扰环境下的回波信号波动效应，从而提高sar地下成像方位分辨率。

2、本发明是这样实现的，一种采用自适应处理的地下目标二维高分辨成像方法，包括以下步骤：

3、s1：构建频控阵sar地下成像信号模型；

4、s2：快时间自适应波束形成；

5、s3：快时间stap处理：对频控阵sar回波信号做快时间stap处理，对频控阵sar的每个脉冲回波数据进行快时间维度的自适应波束形成处理；

6、s4：基于目标散射点定位的目标反射系数估计：

7、s5：地表与地下折射效应估计与补偿；

8、s6：地下高分辨成像折射效应补偿，对地下目标高分辨成像时的折射效应进行补偿；

9、s7：地下目标二维高分辨成像，到此完成折射效应补偿后的快时间-慢时间二维图像，即最后获得的地下目标二维高分辨成像结果。

10、进一步，s1中具体包括：第m个阵元发射载频为：

11、fm＝f0+mδf,m＝0,…,m-1；

12、其中，f0是频控阵sar参考载频，δf是阵元间频偏。假设频控阵sar的基带信号为u(t),0≤t≤tp，其中tp为脉冲宽度，则m阵元频控阵sar的发射信号可表征为：

13、

14、则第n个接收阵元接收到的频控阵sar回波信号可表示为：

15、

16、其中，qn(t,τ)和vn(t,τ)分别表示理想的目标回波信号、不期望的多径干扰和不可避免的接收机噪声；t和τ分别为快时间和慢时间，将全部阵元接收信号统一写成矩阵形式：

17、x(t,τ)＝[x1(t,τ),x2(t,τ),.....,xn(t,τ)]t；

18、t为转置算子。

19、进一步，s2具体包括：为了抑制多径干扰，在快时间维采用自适应波束形成：

20、y(t,τ)＝wh(τ)x(t,τ)；

21、其中，h是共轭转置算子，方位慢时间τ依赖的接收阵列加权向量由下式计算：

22、

23、其中ri+n(τ)是干扰与噪声协方差矩阵，是理想的目标信号向量，相应地，可求解得：

24、

25、其中-1为矩阵求逆运算，多径干扰可以逐脉冲地采用快时间stap来抑制。

26、进一步，s4具体包括：假设目标由k个散射点构成，第k个散射点到第n个接收阵元在方位时间τ时的距离为rkn(τ)时，则有：

27、

28、地下高分辨成像等效为目标散射点在p个方位位置的反射系数，假设高斯噪声，则采用最大似然估计算法：

29、

30、相应地，目标散射系数可估计为：

31、

32、式中，为发射信号在一个脉冲内的能量。

33、进一步，s5具体包括：根据snell定律，ra、rb、φa和φb可由以下方程组求解：

34、

35、ra sinφa+rb sinφb＝xb-xa；

36、ra cosφa＝ya；

37、rb cosφb＝-yb；

38、则折射效应引起的时延可计算为：

39、

40、本发明的另一目的在于提供一种应用所述采用自适应处理的地下目标二维高分辨成像方法的采用自适应处理的地下目标二维高分辨成像系统，包括：

41、信号模型构建模块，用于构建频控阵sar地下成像信号模型；

42、波束形成模块，用于进行快时间自适应波束形成；

43、快时间stap处理模块，用于对频控阵sar回波信号做快时间stap处理，对频控阵sar的每个脉冲回波数据进行快时间维度的自适应波束形成处理；

44、反射系数估计模块，用于基于目标散射点定位进行目标反射系数估计：

45、折射效应估计模块，用于进行地表与地下折射效应估计与补偿；

46、折射效应补偿模块，对地下目标高分辨成像时的折射效应进行补偿；

47、高分辨成像模块，用于完成折射效应补偿后的快时间-慢时间二维图像，即最后获得的地下目标二维高分辨成像结果。

48、本发明的另一目的在于提供一种计算机设备，计算机设备包括存储器和处理器，存储器存储有计算机程序，计算机程序被处理器执行时，使得处理器执行所述的采用自适应处理的地下目标二维高分辨成像方法的步骤。

49、本发明的另一目的在于提供一种计算机可读存储介质，存储有计算机程序，计算机程序被处理器执行时，使得处理器执行所述的采用自适应处理的地下目标二维高分辨成像方法的步骤。

50、本发明的另一目的在于提供一种信息数据处理终端，信息数据处理终端用于实现所述的采用自适应处理的地下目标二维高分辨成像系统。

51、结合上述的技术方案和解决的技术问题，本发明所要保护的技术方案所具备的优点及积极效果为：

52、第一，本发明将新体制频控阵引入地下目标雷达二维高分辨成像，并结合频控阵发射波束的距离依赖性，利用快时间stap有效抑制地下目标成像时的多径传播干扰；同时，利用频控阵波束的距离依赖性，可提升成像分辨率，并可有效补偿地下成像中的地表与地下折射效应带来的成像误差问题。

53、第二，本发明方法相对现有技术在抑制多径传播干扰、提升成像分辨率和补偿地表与地下折射效应方面具有技术优势。

54、第三，本发明提供的新的地下目标二维高分辨成像方法采用了自适应处理的概念，具有以下的优点：

55、1)优化成像信号模型：通过构建频控阵sar地下成像信号模型，可以精确描述地下目标的回波信号特性，从而提高成像结果的准确性。

56、2)提高信号处理效率：利用快时间自适应波束形成和快时间stap处理，可以提高频控阵sar回波信号的处理效率，从而缩短成像时间。

57、3)提高目标定位精度：基于目标散射点定位的目标反射系数估计，可以提高目标的定位精度，从而提高成像结果的精度。

58、4)提高成像质量：通过地表与地下折射效应估计与补偿以及地下高分辨成像折射效应补偿，可以消除成像过程中的折射效应，从而提高成像质量。

59、5)实现高分辨率成像：最后，完成折射效应补偿后的快时间-慢时间二维图像，即最后获得的地下目标二维高分辨成像结果，实现了高分辨率的地下目标成像。

60、综上，这种采用自适应处理的地下目标二维高分辨成像方法，由于其在成像信号模型优化、信号处理效率提高、目标定位精度提高、成像质量提高以及高分辨率成像等方面的显著优点，取得了显著的技术进步。