CN116087953B - 一种星载sar多目标成像参数设计方法 - Google Patents
一种星载sar多目标成像参数设计方法 Download PDFInfo
- Publication number
- CN116087953B CN116087953B CN202310265593.3A CN202310265593A CN116087953B CN 116087953 B CN116087953 B CN 116087953B CN 202310265593 A CN202310265593 A CN 202310265593A CN 116087953 B CN116087953 B CN 116087953B
- Authority
- CN
- China
- Prior art keywords
- imaging
- target
- satellite
- target area
- sar
- 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.)
- Active
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 280
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000013461 design Methods 0.000 title claims abstract description 32
- 208000004350 Strabismus Diseases 0.000 claims abstract description 77
- 238000012216 screening Methods 0.000 claims abstract description 14
- 239000013598 vector Substances 0.000 claims description 104
- 238000004364 calculation method Methods 0.000 claims description 38
- 238000005070 sampling Methods 0.000 claims description 12
- 238000012937 correction Methods 0.000 claims description 5
- 239000011159 matrix material Substances 0.000 description 21
- 238000010586 diagram Methods 0.000 description 6
- 239000011295 pitch Substances 0.000 description 4
- 238000013508 migration Methods 0.000 description 3
- 230000005012 migration Effects 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000012821 model calculation Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 102100029469 WD repeat and HMG-box DNA-binding protein 1 Human genes 0.000 description 1
- 101710097421 WD repeat and HMG-box DNA-binding protein 1 Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 238000002620 method output Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
- G01S13/9041—Squint mode
Abstract
本发明公开了一种星载SAR多目标成像参数设计方法。所述方法以高分辨率和高时效性成像为设计目标,以滑动聚束模式为基本成像模式,满足高分辨率成像需求,利用斜视成像联合高优先级目标先成像,相同优先级目标先方位向由近及远再距离向由近及远的目标筛选原则,提高观测时效性。方法根据卫星精确轨道预报数据和目标区域信息,结合星载SAR斜视滑动聚束模式成像原理,快速确定可成像多目标位置,准确获取各目标区域成像需要的关键时间和斜视角参数,使波束指向在指定的成像时刻精确指向目标区域,满足高分辨率SAR卫星多目标成像工程实现和高时效性观测应用需求。
Description
技术领域
本发明属于星载合成孔径雷达成像领域,具体涉及一种星载SAR多目标成像参数设计方法。
背景技术
星载滑动聚束模式是近年来国内外星载SAR获取高分辨率宽覆盖图像的主要模式,我国多颗星载SAR卫星采用滑动聚束模式获得了大批高质量的高分辨率遥感图像,用于军民领域对地观测。随着两维大扫描相控阵天线技术和卫星大姿态角机动技术发展水平的提高,大斜视星载SAR成像成为可能,从而使星载SAR成像模式设计更加多样,应用更加灵活。
近年来随着对地高效观测应用需求的提出,尤其是在陆地军事目标和海洋舰船监视方面,多目标成像成为星载SAR设计的一个重要方向。陆地军事目标或海洋舰船监视一般要求SAR图像处理具有目标检测和识别功能,要实现目标识别功能要求SAR图像具有较高的分辨率,一般来说要优于0.5m。到目前为止,国内外高分辨率星载SAR都是采用凝视聚束或滑动聚束模式实现,和凝视聚束模式相比,滑动聚束模式由于图像信噪比均匀得到更广泛应用。和正侧视成像相比,前后斜视成像能够充分利用卫星姿态机动能力或者相控阵天线波束电扫能力,在卫星过顶时刻之前或之后对目标成像,增加了可成像目标的数量,提高SAR卫星观测时效性。综上所述,基于高分辨率大斜视成像模式的星载SAR多目标成像既能提高观测效能,又能满足目标检测和识别的应用需求,已成为未来星载SAR发展和应用的主流方向之一。
到目前为止,可能是由于国内外具备大斜视成像能力的SAR卫星较少,采用斜视滑动聚束模式实现高分辨率星载SAR多目标成像的设计方法在公开文献中未见有报道。
发明内容
为解决上述技术问题,本发明提出一种星载SAR多目标成像参数设计方法,方法同时兼顾高分辨率成像和高时效性观测两方面性能。方法包含两部分功能,一部分是单个目标区域成像参数计算,另一部分是可成像目标区域筛选;方法输出设计参数包括可成像目标数量,可成像目标对应的经纬度,所有可成像目标区域对应的成像开始和结束时刻,成像开始和结束斜视角。
进一步地,所述的高时效性含义为成像期间可成像目标数量尽可能多,观测时效性是多目标成像设计要考虑的关键指标。
进一步地,所述的输出设计参数包括可成像目标数量和可成像目标对应的经纬度,用于为地面应用系统提供图像目标信息;所述的输出设计参数包含所有目标区域对应的成像开始和结束时刻,所有目标区域对应的成像开始和结束斜视角,一方面用于SAR卫星生成精确的波束指向控制策略,另一方面用于计算其它SAR成像参数如脉冲重复周期,发射脉冲宽度,回波采样时间,回波帧长等。
为了实现本发明的目的,本发明采用以下技术方案:
设定多目标成像含义,多目标成像含义为在一段连续时间内在不同时刻对不同目标区域进行成像。
设定单个目标区域成像模式,对任意单个目标区域采用斜视滑动聚束模式成像,满足星载SAR高分辨率需求,包含前斜视成像和后斜视成像两种情况,充分利用卫星姿态机动能力或者相控阵天线波束电扫能力,增加多目标观测数量,提高观测时效性。
设定单个目标区域大小按照距离向成像幅宽和方位向成像长度两维成像范围描述,目标处于成像区域中心。不同目标区域之间在地理位置和成像时刻上都没有重叠,目标区域大小由任务规划系统根据目标几何特征确定。
设定目标区域开始成像时刻和结束成像时刻是相对于卫星轨道预报数据起始时刻定义。
设定可成像目标区域筛选原则,所述的可成像目标区域筛选原则为按照优先级高的目标先成像,相同优先级目标先方位向由近及远再距离向由近及远的原则,从预置观测目标区域筛选可成像目标区域。
进一步地,所述的成像目标区域筛选原则是多目标成像设计的关键内容,优先级高的目标先成像是从实际应用需求考虑,首先确保重点关注目标能够成像;相同优先级目标先方位向由近及远再距离向由近及远的目标筛选原则是从提高成像时效性方面考虑而提出,实现过程为沿方位向分布的目标成像时间有冲突时舍弃时间靠后的目标;沿方位向位置分布相同但沿距离向位置分布不同的目标优先对距离近的目标成像,因为,在方位向成像长度相同的情况下,距离近的目标需要的成像时间更短。一段连续时间内单个目标区域成像时间越短,可成像目标区域越多,成像时效性越高。
本发明的星载SAR多目标成像参数设计方法按照顺序包含以下计算步骤:
步骤1,设置输入条件,所述输入条件包括一段连续成像期间卫星精确轨道预报数据,预置观测目标区域信息,SAR系统参数,输入条件由任务规划系统给出。
步骤2,根据空间矢量距离关系,计算成像期间所有预置观测目标对应的最近斜距;
步骤3,根据星载SAR斜视滑动聚束模式成像原理,计算各目标区域成像时间长度;
步骤4,成像期间任一时刻作为成像开始时刻,根据各目标区域对应的成像时间长度,计算对应的成像结束时刻;利用星地斜距矢量和卫星速度矢量空间标量积公式,计算任一目标区域对应的成像开始和结束斜视角;
步骤5,按照设定的可成像目标区域筛选原则,以开始斜视角和结束斜视角绝对量都不大于SAR系统最大可实现斜视角为判断依据,判断当前预置观测目标区域是否为可成像目标区域;
步骤6,输出筛选出的可成像目标区域的数量,所有可成像目标区域对应的目标经纬度,成像开始时刻和成像结束时刻,成像开始斜视角和成像结束斜视角。
进一步地,所述的步骤1中,设置输入条件包含成像期间卫星精确轨道预报数据,要求数据时间长度要包含所有预置观测目标对应的过顶时刻轨道预报数据;要求数据包含WGS84坐标系下的卫星位置矢量和速度矢量;要求数据做过偏航角和俯仰角两维姿态角修正,以保证非斜视成像时卫星速度方向和波束视线方向垂直。
进一步地,所述的步骤1,设置输入条件包含预置观测目标区域信息,所述的预置观测目标区域由任务规划系统给出,一般数量多于可成像目标区域。预置观测目标区域信息包含区域中心点的经纬度,测绘带中心沿方位向起始点、中心点和终止点在WGS84下的位置矢量,方位向成像长度,固定波束指向对应的波束地面速度;要求目标位于场景中心,即目标区域中心点经纬度为目标点经纬度。
进一步地,所述的步骤1中,设置输入条件包含SAR系统参数,所述的SAR系统参数包含方位向分辨率,用表示,方位向分辨率展宽系数,用/>表示,方位向有效天线长度,用/>表示,方位向波束宽度,用/>表示,SAR系统最大可实现斜视角,用/>表示。
进一步地,所述的步骤2中,根据空间矢量距离关系,计算成像期间所有预置观测目标对应的最近斜距,是为了计算单点目标对应的波束地面长度,采用逐点遍历求最小值的方法确定成像期间所有预置观测目标对应的最近斜距。
进一步地,所述的步骤3中,根据星载SAR斜视滑动聚束模式成像原理,计算各目标区域成像时间长度,不同目标区域成像时间长度是影响多目标成像数量的关键影响因素,一般来说,单个目标区域成像时间越短,可成像目标数量越多,SAR卫星成像时效性越高。
星载SAR滑动聚束模式成像时间长度几何含义为地面波束按照一定的速度沿卫星飞行方向滑过方位向成像长度加上一个方位向波束宽度对应的地面长度。第个目标区域对应的方位向成像长度用/>表示,一个方位向波束宽度对应的地面长度用/>表示,滑动聚束下的波束地面速度用/>表示,成像时间长度用/>表示,计算公式为:
进一步地,所述的第个目标区域对应的一个波束宽度在方位向地面上的长度的几何含义为斜视成像下方位向波束宽度在地面上的投影长度,也是为了满足方位分辨率要求的点目标对应的方位向波束有效照射长度。斜视成像下斜距的计算要考虑斜视角引起的距离徙动影响,斜视角越大距离徙动越大,/>的计算在此用SAR系统最大可实现斜视角计算斜距,以保证方位分辨率设计值优于指标要求,/>计算公式为;
进一步地,所述的第个目标区域对应的滑动聚束下的波束地面速度/>用于计算对应目标区域的成像时间长度。根据星载SAR成像原理,和固定波束指向下的波束地面速度相比,滑动聚束模式成像通过降低波束地面速度,增加单个点目标合成孔径时间获得方位向高分辨率,滑动聚束模式下方位分辨率比固定波束指向下的分辨率提高几倍,则对应的波束地面速度降低几倍。/>
进一步地,所述的SAR方位向分辨率展宽系数是星载SAR斜视滑动聚束模式点目标仿真的主要内容,也是SAR方位向参数设计的关键输入参数,决定SAR图像产品方位向分辨率是否能够达到设计要求,同时/>的取值影响目标区域成像时间长度计算,并最终影响多目标成像数量。方位向分辨率展宽系数和多种影响因素有关,主要考虑方位向方向图加权,成像处理加权,多普勒参数估计误差,成像算法近似误差。
进一步地,根据星载SAR成像原理,单点目标合成孔径时间定义为从方位向方向图右边缘进入到左边缘离开的这段时间为波束有效照射时间,由此定义所述的步骤4中目标区域成像开始时刻为方位向方向图右边缘进入目标区域起始点时刻,目标区域成像结束时刻为方位向方向图左边缘离开目标区域终止点时刻。
进一步地,为了有利于工程上实现SAR成像波束指向控制,所述的步骤4中目标区域成像开始斜视角定义为WGS84坐标系下成像开始时刻天线波束中心视线矢量与卫星速度矢量的夹角,目标区域成像结束斜视角定义为WGS84坐标系下成像结束时刻天线波束中心视线矢量与卫星速度矢量的夹角。
因为天线波束中心视线矢量无法作为输入条件给出,不能通过该参数直接计算出成像开始和结束斜视角,根据星载SAR滑动聚束模式星地几何关系,成像开始斜视角可通过成像开始时刻卫星到目标区域起始点的距离矢量与卫星速度矢量的夹角加上半个方位向波束宽度计算得到;成像结束斜视角通过成像结束时刻卫星到目标区域终止点的距离矢量与卫星速度矢量的夹角减去半个方位向波束宽度计算得到。
进一步地,所述的步骤5 ,以开始斜视角和结束斜视角绝对量都不大于SAR系统最大可实现斜视角为判断依据,从预置观测目标区域筛选可成像目标区域。在开始成像时刻为时刻和对应的成像结束时刻为/>时刻,当前第/>个目标区域为可成像目标区域的条件为:
有益效果:
本发明旨在提出一种高分辨率星载SAR多目标成像关键参数的设计方法,方法从实际应用需求出发,以高分辨率和高时效性成像为设计目标。为满足高分辨率成像需求,本发明提出以滑动聚束模式为基本成像模式。为提高观测时效性,本发明从两方面考虑设计,一方面是充分利用卫星姿态机动能力或者相控阵天线波束电扫能力,提出采用前后斜视滑动聚束模式成像,增加多目标观测数量,提高观测时效性;另一方面采用高优先级目标先成像,相同优先级目标先方位向由近及远再距离向由近及远的目标筛选原则,提高成像时效性。
本发明根据卫星精确轨道预报数据和目标区域位置信息,结合星载SAR斜视滑动聚束模式成像原理,可快速确定可成像多目标位置,准确获取各目标成像需要的关键时间和斜视角参数,满足高分辨率SAR卫星多目标成像工程实现和高时效性观测应用需求。
到目前为止,采用斜视滑动聚束模式实现高分辨率星载SAR多目标成像的设计方法在公开文献中未见有报道,因此,本发明在理论研究和工程实现方面都具有创新意义。本发明已在我国在轨运行的高分12号遥感卫星上得到有效验证,分辨率优于0.5m,基于斜视滑动聚束模式的高分辨率多目标成像模式由此获得推广应用,因此,本发明具有重要的工程实用价值。
附图说明
图1为星载SAR前后斜视滑动聚束模式几何关系示意图;
图2为星载SAR多目标成像星地几何关系示意图;
图3为星载SAR多目标成像第 个目标区域前斜视滑动聚束模式几何关系详细示意图;
图4为WGS84坐标系下星载SAR星地斜距矢量关系示意图;
图5为本发明的星载SAR多目标成像参数设计方法流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图,对本发明进行进一步详细说明,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
本发明的第一个方面是明确本发明所采用的成像模式为斜视滑动聚束模式,用于高分辨率星载SAR成像,包含前斜视和后斜视两种情况,如图1所示,左图为前斜视星载SAR滑动聚束模式成像示意图,成像开始时刻天线波束中心偏向卫星前进方向,右图为后斜视星载SAR滑动聚束模式成像示意图,成像开始时刻天线波束中心偏向卫星前进反方向;卫星轨迹为卫星沿轨道运行的位置;星下点轨迹为卫星在不同位置时刻与地心的连线在地球表面上的交点;波束足印中心轨迹为成像期间波束中心在地面上移动的轨迹;测绘带中心为成像场景中心点。
本发明的第二个方面是明确本发明所包含的星载SAR多目标成像含义,以避免因概念理解不同造成设计内容过于宽泛。本发明的星载SAR多目标成像含义为在一段连续时间内在不同时刻对不同目标区域进行成像,如图2所示,轨道数据起始时刻到结束时刻期间,可对个目标进行成像。
本发明的第三个方面是明确从重点目标先成像和提高成像时效性的实际应用需求考虑,提出一种可成像目标区域筛选原则,即按照优先级高的目标先成像,相同优先级按照目标先方位向由近及远再距离向由近及远的原则,从预置观测目标区域筛选可成像目标区域。所述的目标区域筛选原则实现过程为,沿方位向分布的目标成像时间有冲突时舍弃时间靠后的目标,沿方位向位置分布相同但沿距离向位置分布不同的目标优先对距离近的目标成像;以开始斜视角和结束斜视角绝对量都不大于SAR系统最大可实现斜视角为判断依据,确定预置观测目标区域是否为可成像目标区域。
本发明的第四个方面是明确星载SAR多目标成像参数设计方法的关键参数仅包括可成像多目标数量,每个可成像目标区域对应的目标经纬度,成像开始时刻,成像结束时刻,成像开始斜视角,成像结束斜视角,不涉及星载SAR成像工程实现所必须的其他参数。
如图5所示,本发明的第五个方面是约定三项输入参数作为已知条件,第一项输入参数为一段连续成像时间内的高精度卫星轨道预报数据,时间包括所有预置观测目标过顶时刻,轨道预报数据含按照0.001s为时间间隔的WGS84坐标系下的卫星位置矢量矩阵和速度矢量矩阵; 第二项输入参数为所有预置观测目标区域信息,包含区域中心点的经纬度,测绘带中心沿方位向起始点、中心点和终止点在WGS84下的位置矢量,方位向成像长度,固定波束指向对应的波束地面速度;第三项输入参数为SAR系统参数,包含方位向分辨率,方位向分辨率展宽系数,方位向有效天线长度,方位向波束宽度,SAR系统最大可实现斜视角。
本发明的第六个方面是利用所述的已知条件,根据目标区域分布和如图2、图3和图4所示的星载SAR斜视滑动聚束模式几何关系、成像原理以及星地几何关系推导相关参数计算公式。
假设高精度卫星轨道预报数据共有个采用时刻,/>为大于1的整数;假设所述的预置观测目标区域数量为/>个,/>为大于或等于1的整数,任一目标区域用/>()目标区域表示。假设/>,/>,/>分别为高精度卫星轨道预报数据中的三个时刻。
如图2和3所示,本发明的星载SAR多目标成像星地几何关系,所需变量为:
为第/>个目标卫星过顶时刻对应的最近斜距,/>为第/>个目标与轨道预报数据中第/>个采样时刻卫星位置之间的距离,/>为第/>个目标的位置矢量,为卫星在第/>时刻的位置矢量,/>为第/>个目标区域对应的方位向成像长度,为第/>个目标区域对应的一个波束宽度在方位向地面长度,/>为第/>个目标区域的成像时间长度,/>为第/>个目标区域在/>时刻的斜距矢量,/>为第/>个目标区域在/>时刻的斜距矢量,/>和/>分别为第/>个目标区域起始和终止位置矢量,/>和/>分别为为/>时刻和/>时刻卫星位置矢量,/>为/>时刻卫星速度矢量,/>为/>时刻卫星速度矢量,/>为第/>个目标区域在/>时刻对应的成像开始斜视角,/>为第/>个目标区域在/>时刻对应的成像结束斜视角。
如图2、图3和图4所示,多目标成像相关计算公式推导如下:
(1)目标最近斜距:
第个目标卫星过顶时刻对应的最近斜距用/>表示。第/>个目标与轨道预报数据中第/>个采样时刻卫星位置之间的距离用/>表示。根据空间两点间距离矢量关系计算卫星任一位置时刻与目标之间的斜距矢量,对斜距矢量求模得到卫星与目标间的斜距,/>计算公式为:
(2)一个波束宽度在方位向地面上的长度:
第个目标区域对应的一个波束宽度在方位向地面长度用/>表示,含义为方位向一个波束宽度在地面上的投影长度。斜视成像下斜距的计算要考虑斜视角引起的距离徙动影响, 本发明的/>的计算用SAR系统最大可实现斜视角计算斜距,以保证方位分辨率设计值优于指标要求,/>计算公式为;
(3)滑动聚束模式成像对应的波束地面速度:
第个目标区域对应的滑动聚束模式下的波束地面速度用/>表示,根据星载SAR成像原理,和固定波束指向下的波束地面速度相比,滑动聚束模式成像通过降低波束地面速度,增加单个点目标合成孔径时间获得方位向高分辨率。第/>个目标区域对应的固定波束指向下波束地面速度用/>表示,考虑方位向分辨率展宽系数,用/>表示,/>计算公式为:
(4)目标区域成像时间长度:
第个目标区域的成像时间长度用/>表示,星载SAR滑动聚束模式成像时间长度含义为地面波束按照一定的速度沿卫星飞行方向滑过方位向成像长度加上一个方位向波束宽度对应的地面长度所用的时间长度,第/>个目标区域对应的方位向成像长度用表示,/>计算公式为:/>
(5)目标区域成像结束时刻:
(6)目标区域成像开始和结束时刻星地斜距矢量:
(7)目标区域成像开始斜视角:
根据星载SAR滑动聚束模式星地几何关系,成像开始斜视角可通过成像开始时刻卫星到目标区域起始点的距离矢量与卫星速度矢量的夹角加上半个方位向波束宽度计算得到。任一卫星位置采样时刻作为第/>个目标区域成像开始时刻,第/>个目标区域在开始时刻的成像开始斜视角用/>表示。利用空间矢量标量积公式计算两矢量间的夹角,计算公式为:
(8)目标区域成像结束斜视角:
根据星载SAR滑动聚束模式星地几何关系,成像结束斜视角通过成像结束时刻卫星到目标区域终止点的距离矢量与卫星速度矢量的夹角减去半个方位向波束宽度计算得到。任一卫星位置采样时刻作为第/>个目标区域成像结束时刻,第/>个目标区域在结束时刻对应的成像结束斜视角用/>表示,计算公式为:
(9)可成像目标区域判断依据:
一般SAR卫星都具备方位向波束向前或向后扫描的能力,用正负号来区分,星载SAR成像需要的开始和结束斜视角绝对量不能超过SAR系统最大设计斜视角。第个目标区域在/>时刻作为成像开始时刻,对应的成像结束时刻为/>时刻,第/>个目标区域的开始和结束斜视角分别为/>和/>,目标区域可成像需要满足的条件为:
为了实现本发明的目的,本发明的第七个方面是根据所述的第一到六方面的技术方案,已知条件,相关计算公式,确定基于斜视滑动聚束模式的高分辨率星载SAR多目标成像参数设计方法实现流程,下面结合图5对本发明的星载SAR设计方法相关计算步骤进行详细说明。
如图5所示,本发明的一种星载SAR多目标成像参数设计方法从输入已知条件开始到设计参数输出,共包含14步计算或判断。
第1步为输入已知条件,所述的已知条件中包含星载SAR多目标成像参数设计方法需要预知的基本相关参数,所述的基本相关参数分为三项,第一项是卫星轨道预报数据,第二项是预置观测目标区域信息,第三项是SAR系统参数,第一项和第二项由任务规划系统给出,第三项由SAR系统给出。所述的基本相关参数包括:
预置观测目标区域优先级按照从高到低对应正整数从1开始排序。
以上任务规划输入参数由任务规划系统给出。
根据目标区域优先级排序,按照优先级高的目标先成像,相同优先级目标按照先方位向由近及远再距离向由近及远的原则,筛选可成像目标区域;
第8步,以时刻作为第/>个预置观测目标区域成像开始时刻,按照公式(6)计算成像结束时刻/>;在此需要说明的是,对于不同的预置观测目标区域,/>和/>取值范围不同,需要重新设置,优先级低的预置观测目标区域在设置/>和/>时不能取优先级高的可成像目标区域对应的成像开始和结束时刻范围内的时刻点,以保证不同可成像目标区域在成像时间上没有冲突。
第12步, 如果第个预置观测目标区域为可成像目标区域,可成像目标区域数量加1,保存当前目标经纬度,成像开始和结束时刻,成像开始和结束斜视角;当前目标经纬度从任务规划输入参数的预置观测目标区域中心点经纬度/>和/>矩阵中查找得到。
输出全部参数包括可成像多目标数量,每个可成像目标区域对应的目标经纬度,成像开始时刻,成像结束时刻,成像开始斜视角,成像结束斜视角。
本发明未详细阐述的部分属于本领域公知技术。本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (9)
1.一种星载SAR多目标成像参数设计方法,其特征在于,包括:
设定单个目标区域成像模式,所述单个目标区域成像模式采用斜视滑动聚束模式;
设定多目标成像含义,所述多目标成像含义为在一段连续时间内在不同时刻对不同目标区域进行成像;
设定多目标成像设计关键参数,所述多目标成像设计关键参数包括可成像目标数量、每个可成像目标区域对应的目标经纬度、成像开始时刻、成像结束时刻、成像开始斜视角和成像结束斜视角;
设定可成像目标区域筛选原则,所述的可成像目标区域筛选原则为优先级高的目标先成像,相同优先级目标按照先方位向由近及远再距离向由近及远的原则,从预置观测目标区域筛选可成像目标区域;所述预置观测目标区域的信息包含区域中心点的经纬度,测绘带中心沿方位向起始点、中心点和终止点在 WGS84 下的位置矢量,方位向成像长度,固定波束指向对应的波束地面速度;目标位于场景中心,即目标区域中心点经纬度为目标点经纬度;
具体包括如下步骤:
步骤1,设置输入条件,所述输入条件包括一段连续成像期间卫星精确轨道预报数据、预置观测目标区域信息和SAR系统参数;所述卫星精确轨道预报数据的数据时间长度包含所有预置观测目标对应的过顶时刻轨道预报数据,且包含 WGS84 坐标系下的卫星位置矢量和速度矢量,并做过偏航角和俯仰角两维姿态角修正,以保证非斜视成像时卫星速度方向和波束视线方向垂直;
步骤2,根据空间矢量距离关系,计算成像期间所有预置观测目标对应的最近斜距;
步骤3,根据星载SAR斜视滑动聚束模式成像原理,计算各目标区域成像时间长度;
步骤4,利用星地斜距矢量和卫星速度矢量空间标量积公式,计算成像期间任一时刻、任一目标区域对应的成像开始斜视角和成像结束斜视角;
步骤5,按照设定的可成像目标区域筛选原则,以开始斜视角和结束斜视角绝对量都不大于SAR系统最大可实现斜视角为判断依据,判定当前预置观测目标区域是否为可成像目标区域;
步骤6,输出多目标成像设计关键参数。
2.如权利要求1所述一种星载SAR多目标成像参数设计方法,其特征在于,所述的单个目标区域成像模式包含前斜视成像和后斜视成像两种情况。
3.如权利要求1所述一种星载SAR多目标成像参数设计方法,其特征在于,单个目标区域的大小按照距离向成像幅宽和方位向成像长度描述,目标处于成像区域中心;不同目标区域之间在地理位置和成像时刻上均没有重叠,目标区域的大小由地面任务规划系统根据目标几何特征确定。
4.如权利要求1所述一种星载SAR多目标成像参数设计方法,其特征在于,所述的可成像目标区域筛选原则实现过程为:沿方位向分布的目标区域成像时刻有冲突时舍弃时间靠后的目标,沿方位向位置分布相同但沿距离向位置分布不同的目标区域优先对距离近的目标区域成像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310265593.3A CN116087953B (zh) | 2023-03-20 | 2023-03-20 | 一种星载sar多目标成像参数设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310265593.3A CN116087953B (zh) | 2023-03-20 | 2023-03-20 | 一种星载sar多目标成像参数设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116087953A CN116087953A (zh) | 2023-05-09 |
CN116087953B true CN116087953B (zh) | 2023-06-09 |
Family
ID=86204717
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310265593.3A Active CN116087953B (zh) | 2023-03-20 | 2023-03-20 | 一种星载sar多目标成像参数设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116087953B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102007031020B3 (de) * | 2007-07-04 | 2008-12-24 | Deutsches Zentrum für Luft- und Raumfahrt e.V. | Verfahren zur Verarbeitung von TOPS(Terrain Observation by Progressive Scan)-SAR(Synthetic Aperture Radar)-Rohdaten und Verwendung des Verfahrens |
CN115792907A (zh) * | 2023-02-07 | 2023-03-14 | 中国科学院空天信息创新研究院 | 星载sar斜视滑动聚束模式方位向成像参数设计方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102176016B (zh) * | 2011-01-25 | 2012-12-19 | 北京航空航天大学 | 一种大斜视滑动聚束sar成像处理方法 |
CN103576152B (zh) * | 2013-11-22 | 2016-04-06 | 中国科学院电子学研究所 | 一种滑动聚束合成孔径雷达及其实现方法和装置 |
CN107505615A (zh) * | 2017-08-22 | 2017-12-22 | 西安电子科技大学 | 星载sar滑动聚束成像模式星体姿态设计方法 |
CN110058232B (zh) * | 2019-04-19 | 2021-04-13 | 北京空间飞行器总体设计部 | 一种星载sar大斜视滑动聚束模式回波信号方位预处理方法及系统 |
CN114137519A (zh) * | 2021-09-28 | 2022-03-04 | 北京空间飞行器总体设计部 | 一种高分辨率sar成像参数计算方法 |
-
2023
- 2023-03-20 CN CN202310265593.3A patent/CN116087953B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102007031020B3 (de) * | 2007-07-04 | 2008-12-24 | Deutsches Zentrum für Luft- und Raumfahrt e.V. | Verfahren zur Verarbeitung von TOPS(Terrain Observation by Progressive Scan)-SAR(Synthetic Aperture Radar)-Rohdaten und Verwendung des Verfahrens |
CN115792907A (zh) * | 2023-02-07 | 2023-03-14 | 中国科学院空天信息创新研究院 | 星载sar斜视滑动聚束模式方位向成像参数设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN116087953A (zh) | 2023-05-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA1212166A (en) | Range/azimuth/elevation angle ship imaging for ordnance control | |
US9110170B1 (en) | Terrain aided navigation using multi-channel monopulse radar imaging | |
US7532150B1 (en) | Restoration of signal to noise and spatial aperture in squint angles range migration algorithm for SAR | |
CN115792907B (zh) | 星载sar斜视滑动聚束模式方位向成像参数设计方法 | |
CN109507665A (zh) | 一种基于星载ais实时信息引导的星上自主成像方法 | |
CN107300699B (zh) | 基于敏捷合成孔径雷达卫星姿态机动的马赛克模式实现方法 | |
JPS62231188A (ja) | 船舶の合成開口レ−ダ結像方法および装置 | |
US8212714B1 (en) | Using doppler radar images to estimate aircraft navigational heading error | |
CN107390181B (zh) | 一种基于多波束扫描的雷达高分辨成像方法 | |
CN107064880A (zh) | 分布式多基雷达收发波束同步和波束控制的高精度方法 | |
Kraus et al. | TerraSAR-X staring spotlight mode optimization and global performance predictions | |
CN114545411B (zh) | 一种基于工程实现的极坐标格式多模高分辨sar成像方法 | |
KR101958547B1 (ko) | 단일 경로와 단일 안테나를 이용한 영상 레이다의 고도 측정 시스템 및 방법 | |
CN114137519A (zh) | 一种高分辨率sar成像参数计算方法 | |
US6987479B1 (en) | Conformal range migration algorithm (CRMA) “KARMA” | |
CN116087953B (zh) | 一种星载sar多目标成像参数设计方法 | |
Bezvesilniy et al. | Synthetic aperture radar systems for small aircrafts: Data processing approaches | |
KR102185307B1 (ko) | Sar 영상의 객체 응답 초해상도화 방법 및 객체 응답 초해상도화 장치 | |
CN109738890B (zh) | 一种基于弹载双基sar距离多普勒图像生成地距图的方法 | |
CN116148891A (zh) | 基于顺轨干涉的星地双基海面二维流场测量方法 | |
Wang et al. | Study on motion compensation method for W-band UAV MISAR real-time imaging | |
CN114185047B (zh) | 一种基于最优极坐标变换的双基sar动目标重聚焦方法 | |
Fan et al. | On the processing of Gaofen-3 spaceborne dual-channel sliding spotlight SAR data | |
Nie et al. | A quadtree beam-segmenting based wide-swath SAR polar format algorithm | |
CN113156435A (zh) | 一种基于嵌入式gpu的弹载sar前侧视时域成像方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |