CN102998655B - 一种合成孔径雷达扫描工作模式最优波位选择方法 - Google Patents

一种合成孔径雷达扫描工作模式最优波位选择方法 Download PDF

Info

Publication number
CN102998655B
CN102998655B CN201210480438.5A CN201210480438A CN102998655B CN 102998655 B CN102998655 B CN 102998655B CN 201210480438 A CN201210480438 A CN 201210480438A CN 102998655 B CN102998655 B CN 102998655B
Authority
CN
China
Prior art keywords
ripple
prime
child
wave position
ripple position
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201210480438.5A
Other languages
English (en)
Other versions
CN102998655A (zh
Inventor
丁泽刚
龙腾
曾涛
尹伟
汪精华
杨文付
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201210480438.5A priority Critical patent/CN102998655B/zh
Publication of CN102998655A publication Critical patent/CN102998655A/zh
Application granted granted Critical
Publication of CN102998655B publication Critical patent/CN102998655B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明针对扫描工作模式下,波位的选取仅能通过手动选择、波位选择结果非参数最优结果的问题,提出一种合成孔径雷达扫描工作模式最优波位选择方法。步骤一、根据距离向测绘区域、等效后向散射系数最大值、模糊比最大值参数的限制,以及可选用的i组条带宽度,分别进行条带工作模式波位设计,得出i组不同条带宽度的波位,并将其并集作为扫描工作模式的待选波位集合;步骤二、初始波位选择;步骤三、波位树的构建;步骤四、波位树的遍历;步骤五、选取最优波位组合。

Description

一种合成孔径雷达扫描工作模式最优波位选择方法
技术领域
本发明涉及一种合成孔径雷达扫描工作模式最优波位选择方法,属于合成孔径雷达技术领域。
背景技术
合成孔径雷达(SAR)是一种全天时、全天候的高分辨率的微波遥感成像雷达,可安装在飞机、卫星、宇宙飞船等飞行平台上。在环境监测、海洋观测、资源勘探、农作物估产、测绘和军事等方面的应用上具有独特的优势,可发挥其他遥感手段难以发挥的作用。
扫描工作模式是SAR的基础工作模式之一,主要用于距离向宽测绘带成像。当SAR工作在扫描模式时,天线在一个波束指向上发射脉冲并接收相应回波,形成当前子条带的回波数据块,然后跳转到另一波束指向,继续发射并接收回波,获取新子条带的回波数据块。不断重复此过程便可获得相互平行的子条带的回波数据块,然后经过成像处理和拼接,就可以获得较宽的测绘带结果。
在进行扫描工作模式的波位设计时,首先需要进行的条带模式的波位设计,得出几组不同测绘带宽度的波位设计结果,然后在从中选取能够完全覆盖测绘区域的波位。目前,上述选择过程是利用斑马图手动进行的,选择过程中不能直观的比较不同选择方式的性能参数优劣,选择的结果依赖于操作人员的个人判断,因此波位选择的最终结果具有很大的随意性,并不一定是性能参数最优的结果。
发明内容
本发明的目的是为了弥补现有技术的不足,针对扫描工作模式下,波位的选取仅能通过手动选择、波位选择结果非参数最优结果的问题,提出一种合成孔径雷达扫描工作模式最优波位选择方法。
本发明方法是通过下述技术方案实现的:
一种合成孔径雷达扫描工作模式最优波位选择方法,包括如下步骤:
步骤一、根据距离向测绘区域、等效后向散射系数(NESZ)最大值、模糊比最大值参数的限制,以及可选用的i组条带宽度{w1 w2...wi},分别进行条带工作模式波位设计,得出i组不同条带宽度的波位 x 11 x 12 . . . x 1 n 1 , · · · , x i 1 x i 2 . . . x in i , 并将其并集 x 11 x 12 · · · x 1 n 1 · · · x i 1 x i 2 · · · x in i 作为扫描工作模式的待选波位集合;
步骤二、初始波位选择:在所述的待选波位集合 x 11 x 12 · · · x 1 n 1 · · · x i 1 x i 2 · · · x in i 中选取波位覆盖区域包含距离向测绘区域近端的波位集合 y 1 y 2 . . . y n root , 作为扫描工作模式波位选择的待选起始波位;其余波位的集合 z 1 z 2 . . . z n child 作为待选后续波位;
步骤三、波位树的构建:从起始波位集合 y 1 y 2 . . . y n root 中依次选取元素yk,k=1,2…nroot作为根节点,分别构建nroot个波位树,其中子节点与父节点波位含有重合的覆盖区域;
步骤四、波位树的遍历:对于已构建的nroot个波位树依次进行遍历,在波位树的各个支路中找出能够满足测绘带覆盖要求的波位组合,即该支路的根节点波位的近端下视角小于要求覆盖区域的近端下视角、叶节点波位的远端下视角大于要求覆盖区域的远端下视角,并分别计算对应波位的NESZ、方位向模糊度、距离向模糊度;
步骤五、选取最优波位组合:从步骤四的遍历结果中,选取参数最优的波位组合,作为扫描工作模式的波位设计结果。
步骤三中构建波位树采用以下方法:首先将待选后续波位 z 1 z 2 . . . z n child 分为互补的两个集合 y k 1 y k 2 . . . y kn k z 1 ′ z 2 ′ . . . z ′ n child , 其中 y k 1 y k 2 . . . y kn k 作为根节点yk的子节点集合, z 1 ′ z 2 ′ . . . z ′ n child 作为新的待选后续波位集合,将 z 1 z 2 . . . z n child 中波位覆盖区域与根节点yk波位覆盖区域有重叠的元素集合作为 y k 1 y k 2 . . . y kn k , 剩余元素集合作为 z 1 ′ z 2 ′ . . . z ′ n child , 然后依照上述原则从新的待选后续波位集合 z 1 ′ z 2 ′ . . . z ′ n child 中分别选取子节点 y k 1 y k 2 . . . y kn k 的子节点,并产生新的待选后续波位集合,重复上述过程,直至待选后续波位集合为空集。
本发明的有益效果:
本方法摆脱了手动选择扫描工作模式的限制,能够利用程序对覆盖测绘区域的所有波位组合进行遍历,从中选取参数最优的波位组合。
附图说明
图1条带模式波位设计结果;
图2以X1为根节点建立的波位树;
图3以X2为根节点建立的波位树;
图4为满足覆盖要求的波位示意图。
具体实施方式
下面结合附图和实施例对本发明方法进一步详细说明。
本实施例中,卫星运行在椭圆轨道上,轨道的各参数具体如下所示:
地球半径:6371.004km
轨道高度:800km
入射角范围:20°~40°
下视角:30°
条带模式分辨率:25m
条带宽度:135km/105km
地球自转角速度:7.292115×10-5rad/s
采用本发明所述的一种合成孔径雷达扫描工作模式最优波位选择方法选择该卫星扫描工作模式的最优波位并仿真,其具体步骤为:
步骤一,根据距离向测绘区域、NESZ、模糊度等参数的限制,以及可选用的i组条带宽度{w1 w2...wi},分别进行条带工作模式波位设计,得出i组不同条带宽度的波位设计结果 x 11 x 12 . . . x 1 n 1 , · · · , x i 1 x i 2 . . . x in i , 并将其并集 x 11 x 12 · · · x 1 n 1 · · · x i 1 x i 2 · · · x in i 作为扫描工作模式的待选波位集合。根据轨道高度800km、入射角范围20°~40°可以确定距离向测绘带区域,NESZ小于-19dB,模糊度小于-20dB,根据上述条件分别进行条带宽度135km与105km的条带工作模式波位设计,得到的波位如图1所示。条带宽度135km的波位集合为{x2 x3 x5 x8},条带宽度105km的波位集合为{x1 x4 x6 x7},扫描工作模式的待选波位集合为{x1 x2 x3 x4 x5 x6 x7 x8}。
步骤二,初始波位选择。在待选波位集合 x 11 x 12 · · · x 1 n 1 · · · x i 1 x i 2 · · · x in i 中选取波位覆盖区域包含距离向测绘区域近端的波位集合 y 1 y 2 . . . y n root , 作为扫描工作模式波位选择的待选起始波位;其余波位的集合 z 1 z 2 . . . z n child 作为待选后续波位。根据覆盖特性,可以选定初始波位集合为{x1 x2},待选后续波位集合为{x3 x4 x5 x6 x7 x8}。
步骤三,波位树的构建。从起始波位集合 y 1 y 2 . . . y n root 中依次选取元素yk(k=1,2…nroot)作为根节点,分别构建nroot个波位树。具体的做法是,根据波位覆盖区域是否与根节点yk波位覆盖区域有重叠,将后续波位集合 z 1 z 2 . . . z n child 分为互补的两个集合,有覆盖区域重叠的波位集合 y k 1 y k 2 . . . y kn k 作为根节点yk的子节点,无覆盖区域重叠的波位集合 z 1 ′ z 2 ′ . . . z ′ n child 作为新的待选后续波位集合。然后依照上述原则从新的待选后续波位集合 z 1 ′ z 2 ′ . . . z ′ n child 中分别选取子节点 y k 1 y k 2 . . . y kn k 的子节点,并产生新的待选后续波位集合,重复上述过程,直至待选后续波位集合为空集。
以起始波位为x1为例,在待选后续波位集合{x3 x4 x5 x6 x7 x8}中,与x1波位有覆盖区域重叠的波位集合为{x3 x4},以上述两个波位作为根节点x1的子节点;与x1波位无覆盖区域重叠的波位集合为{x5 x6 x7 x8},以上述四个波位作为新的待选后续波位集合。对于波位集合{x3 x4},待选后续波位集合{x5 x6 x7 x8}中与x3波位有覆盖区域重合的为波位{x5 x6},与x4波位有重合的为{x5 x6},故而x3与x4的子节点均为{x5 x6},新的待选后续波位集合为{x7 x8}。对于波位集合{x5 x6},待选后续波位集合{x7 x8}中与x5波位有覆盖区域重合的为波位{x7 x8},与x6波位有重合的为波位{x7},故而x5子节点集合为{x7 x8}、x6的子节点集合为{x7},新的待选后续波位集合为空集。以x1为根节点的波位树建立完毕,建立的波位树如图2所示。同理可以建立以x2为起始波位的波位树,如图3所示。
步骤四,波位树的遍历。对于已构建的nroot个波位树依次进行遍历,找出能够满足测绘带覆盖要求的波位组合,并分别计算其相应的参数,如NESZ、方位向模糊度、距离向模糊度等。
分别遍历以x1为根节点和以x2为根节点建立的波位树,得到满足测绘带覆盖要求的十二组波位组合,分别为:1、{x1,x3,x5,x7},2、{x1,x3,x5,x8},3、{x1,x3,x6,x7},4、{x1,x4,x5,x7},5、{x1,x4,x5,x8},6、{x1,x4,x6,x7},7、{x2,x3,x5,x7},8、{x2,x3,x5,x8},9、{x2,x3,x6,x7},10、{x2,x4,x5,x7},11、{x2,x4,x5,x8},12、{x2,x4,x6,x7}。
计算以上十二组波位组合的等效后向散射系数,方位向模糊比和距离向模糊比等参数。表1为各项性能指标的平均值。
表1 各项参数平均值列表
波位号 重叠宽度(km) AASR(dB) RASR(dB) NESZ中心(dB)
1 16587.5 -22.2275 -30.44 -28.8425
2 16712.5 -22.4125 -31.15 -29.4
3 20687.5 -21.2125 -32.48 -30.0225
4 4712.5 -22.6975 -30.9075 -29.235
5 4587.5 -22.5125 -30.1975 -28.6775
6 8687.5 -21.4975 -32.2375 -29.8575
7 16587.5 -23.6375 -27.5125 -28.34
8 18400 -22.375 -26.8775 -28.1925
9 20687.5 -22.6225 -29.5525 -29.52
10 6116.667 -23.4433 -26.3167 -28.4933
11 6283.333 -23.69 -27.2633 -29.2367
12 11583.33 -22.09 -29.0367 -30.0667
步骤五,选取最优波位组合。从步骤四的遍历结果中,选取性能指标最优的波位组合,作为扫描工作模式的波位设计结果。
参考表1各项参数结果,本实例依照重叠率最小的标准选取第四组波位作为最佳组合。各项工程对波位要求不同,具体标准可按实际情况确定。
以上所述的具体描述,对发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (1)

1.一种合成孔径雷达扫描工作模式最优波位选择方法,其特征在于,包括如下步骤:
步骤一、根据距离向测绘区域、等效后向散射系数NESZ最大值、模糊比最大值参数的限制,以及可选用的i组条带宽度 w 1 w 2 . . . w i ,分别进行条带工作模式波位设计,得出i组不同条带宽度的波位 x 11 x 12 . . . x 1 n 1 , . . . , x i 1 x i 2 . . . x in i ,并将其并集 x 11 x 12 . . . x 1 n 1 . . . x i 1 x i 2 . . . x in i 作为扫描工作模式的待选波位集合;
步骤二、初始波位选择:在所述的待选波位集合 x 11 x 12 . . . x 1 n 1 . . . x i 1 x i 2 . . . x in i 中选取波位覆盖区域包含距离向测绘区域近端的波位集合 y 1 y 2 . . . y n root ,作为扫描工作模式波位选择的待选起始波位;其余波位的集合 z 1 z 2 . . . z n child 作为待选后续波位;
步骤三、波位树的构建:首先将待选后续波位 z 1 z 2 . . . z n child 分为互补的两个集合 y k 1 y k 2 . . . y kn k z 1 ′ z 2 ′ . . . z ′ n child ,其中 y k 1 y k 2 . . . y kn k 作为根节点yk的子节点集合, z 1 ′ z 2 ′ . . . z ′ n child 作为新的待选后续波位集合,将 z 1 z 2 . . . z n child 中波位覆盖区域与根节点yk波位覆盖区域有重叠的元素集合作为 y k 1 y k 2 . . . y kn k ,剩余元素集合作为 z 1 ′ z 2 ′ . . . z ′ n child ,然后依照上述原则从新的待选后续波位集合 z 1 ′ z 2 ′ . . . z ′ n child 中分别选取子节点 y k 1 y k 2 . . . y k n k 的子节点,并产生新的待选后续波位集合,重复上述过程,直至待选后续波位集合为空集;
步骤四、波位树的遍历:对于已构建的nroot个波位树依次进行遍历,在波位树的各个支路中找出能够满足测绘带覆盖要求的波位组合,即该支路的根节点波位的近端下视角小于要求覆盖区域的近端下视角、叶节点波位的远端下视角大于要求覆盖区域的远端下视角,并分别计算对应波位的等效后向散射系数NESZ、方位向模糊度、距离向模糊度;
步骤五、选取最优波位组合:从步骤四的遍历结果中,选取参数最优的波位组合,作为扫描工作模式的波位设计结果。
CN201210480438.5A 2012-11-22 2012-11-22 一种合成孔径雷达扫描工作模式最优波位选择方法 Expired - Fee Related CN102998655B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210480438.5A CN102998655B (zh) 2012-11-22 2012-11-22 一种合成孔径雷达扫描工作模式最优波位选择方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210480438.5A CN102998655B (zh) 2012-11-22 2012-11-22 一种合成孔径雷达扫描工作模式最优波位选择方法

Publications (2)

Publication Number Publication Date
CN102998655A CN102998655A (zh) 2013-03-27
CN102998655B true CN102998655B (zh) 2014-07-30

Family

ID=47927464

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210480438.5A Expired - Fee Related CN102998655B (zh) 2012-11-22 2012-11-22 一种合成孔径雷达扫描工作模式最优波位选择方法

Country Status (1)

Country Link
CN (1) CN102998655B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105487051B (zh) * 2015-11-10 2019-01-25 中国空空导弹研究院 雷达角度搜索的扫描方法
CN107144835B (zh) * 2017-04-28 2020-01-21 安徽四创电子股份有限公司 一种低空目标监视方法
CN112198482B (zh) * 2020-09-18 2023-09-19 北京理工大学 星载超高分辨率sar徙动精确控制的连续变脉冲间隔生成方法
CN112684448B (zh) * 2020-12-22 2022-06-03 上海卫星工程研究所 面向多参数组合的星载sar波位计算及生成方法及系统
CN112946651B (zh) * 2021-04-23 2023-10-27 成都汇蓉国科微系统技术有限公司 一种基于分布式sar的空中协同感知系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5548290A (en) * 1994-03-23 1996-08-20 National Space Development Agency Of Japan Method of designing transmission power in synthetic aperture radar
DE102004061486B4 (de) * 2004-12-21 2007-04-05 Eads Deutschland Gmbh Verfahren zur Signaldatenverarbeitung eines fluggerätegetragenen Radars mit synthetischer Apertur und eine Vorrichtung dafür
CN101526614B (zh) * 2009-04-03 2011-09-07 北京理工大学 基于子孔径和等效散射体的sar回波快速模拟方法
CN102520403A (zh) * 2011-11-30 2012-06-27 北京理工大学 一种改进的基于频域频谱重建的频率步进sar成像方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5548290A (en) * 1994-03-23 1996-08-20 National Space Development Agency Of Japan Method of designing transmission power in synthetic aperture radar
DE102004061486B4 (de) * 2004-12-21 2007-04-05 Eads Deutschland Gmbh Verfahren zur Signaldatenverarbeitung eines fluggerätegetragenen Radars mit synthetischer Apertur und eine Vorrichtung dafür
CN101526614B (zh) * 2009-04-03 2011-09-07 北京理工大学 基于子孔径和等效散射体的sar回波快速模拟方法
CN102520403A (zh) * 2011-11-30 2012-06-27 北京理工大学 一种改进的基于频域频谱重建的频率步进sar成像方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
丁泽刚等.基于方位向去斜的星载滑动聚束两步成像处理算法研究.《兵工学报》.2012,第33卷(第1期),95-102.
基于方位向去斜的星载滑动聚束两步成像处理算法研究;丁泽刚等;《兵工学报》;20120131;第33卷(第1期);95-102 *
数字馈相条件下相控阵天线最优扫描波位分析;李秋生等;《制导与引信》;20090930;第30卷(第3期);6-10,35 *
李秋生等.数字馈相条件下相控阵天线最优扫描波位分析.《制导与引信》.2009,第30卷(第3期),6-10,35.

Also Published As

Publication number Publication date
CN102998655A (zh) 2013-03-27

Similar Documents

Publication Publication Date Title
CN102998655B (zh) 一种合成孔径雷达扫描工作模式最优波位选择方法
CN110132238B (zh) 用于地形影像数字高程模型的无人机测绘方法
Fujii et al. An overview of developments and applications of oceanographic radar networks in Asia and Oceania countries
CN102313885B (zh) 多维度微波成像系统及方法
Felfili et al. Diversity, floristic and structural patterns of cerrado vegetation in Central Brazil
CN107368094A (zh) 一种无人机植保作业航线规划方法及装置
CN106504362A (zh) 基于无人机的输变电系统巡检方法
CN106503433A (zh) 基于无人机作业的喷洒区域确定方法及装置
CN106657882A (zh) 基于无人机的输变电系统实时监控方法
CN103559493A (zh) 一种线状地物的提取方法
CN104036108A (zh) 一种海面风场和海浪联合反演的方法
CN104714225B (zh) 一种基于广义似然比的动态规划检测前跟踪方法
CN104361338A (zh) 一种基于ENVISAT ASAR、Landsat TM与DEM数据的泥炭沼泽信息提取方法
CN114003981B (zh) 一种基于时空一体数字地球的电磁频谱可视化分析方法
CN106959700A (zh) 一种基于上限置信区间算法的无人机群协同巡逻追踪轨迹规划方法
CN104331588A (zh) 一种海面风场和海浪联合反演的装置
CN103176168A (zh) 一种机载非正侧视阵雷达近程杂波对消方法
CN107102329A (zh) Ka波段TOPS InSAR系统及其工作方法
CN103440489A (zh) 基于像素级sar影像时间序列相似性分析的水体提取方法
CN106872973A (zh) 利用cinrad/sa天气雷达监测和跟踪林火回波的方法
CN109541586A (zh) 一种探地雷达快速探测道路工程混凝土垫层厚度分析方法
Jin et al. [Retracted] Mine Geological Environment Monitoring and Risk Assessment in Arid and Semiarid Areas
CN107871323A (zh) 一种基于全固态雷达的海冰运动信息获取方法
CN103576153A (zh) 一种方位向多波束合成孔径雷达及其实现方法和装置
CN103886582A (zh) 一种利用特征点Voronoi图优选的星载合成孔径干涉雷达影像配准方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C53 Correction of patent for invention or patent application
CB03 Change of inventor or designer information

Inventor after: Ding Zegang

Inventor after: Long Teng

Inventor after: Zeng Tao

Inventor after: Yin Wei

Inventor after: Wang Jinghua

Inventor after: Yang Wenfu

Inventor before: Ding Zegang

Inventor before: Long Teng

Inventor before: Zeng Tao

Inventor before: Yin Wei

Inventor before: Yang Wenfu

COR Change of bibliographic data

Free format text: CORRECT: INVENTOR; FROM: DING ZEGANG LONG TENG CENG TAO YIN WEI YANG WENFU TO: DING ZEGANG LONG TENG CENG TAO YIN WEI WANG JINGHUA YANG WENFU

C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20140730

Termination date: 20141122

EXPY Termination of patent right or utility model