CN103675760B - 一种星载地球同步轨道合成孔径雷达姿态导引方法 - Google Patents
一种星载地球同步轨道合成孔径雷达姿态导引方法 Download PDFInfo
- Publication number
- CN103675760B CN103675760B CN201310651459.3A CN201310651459A CN103675760B CN 103675760 B CN103675760 B CN 103675760B CN 201310651459 A CN201310651459 A CN 201310651459A CN 103675760 B CN103675760 B CN 103675760B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- omega
- angle
- msup
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 15
- 238000005096 rolling process Methods 0.000 claims description 9
- 230000005484 gravity Effects 0.000 claims description 3
- 238000005457 optimization Methods 0.000 abstract 2
- 239000011159 matrix material Substances 0.000 description 5
- 230000008859 change Effects 0.000 description 3
- 230000000007 visual effect Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 208000004350 Strabismus Diseases 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
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/9058—Bistatic or multistatic SAR
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种星载地球同步轨道合成孔径雷达姿态导引方法。使用本发明能够利用较小的姿态导引角实现最优地距分辨率。本发明首先根据已知的轨道参数、下视角等条件获得卫星速度方向与地距分辨椭圆面积的表达式,然后利用地距分辨椭圆面积最小时地距分辨率最优,解算出最优地距分辨率的最优方位角,最后通过俯仰横滚导引或横滚俯仰导引使波束指向最优方位角所在方向,完成最优分辨率姿态导引,姿态导引角度小。
Description
技术领域
本发明涉及合成孔径雷达技术领域,具体涉及一种星载地球同步轨道合成孔径雷达姿态导引方法。
背景技术
合成孔径雷达(SAR)是一种全天时、全天候的高分辨率的微波遥感成像雷达,可安装在飞机、卫星、宇宙飞船等飞行平台上。在环境监测、海洋观测、资源勘探、农作物估产、测绘和军事等方面的应用上具有独特的优势,可发挥其他遥感手段难以发挥的作用。对于SAR如此广泛的应用,它的应用质量好坏一般取决于SAR的图像分辨率。
星载地球同步轨道(GEO)合成孔径雷达是指位于地球同步轨道、以卫星为载体的合成孔径雷达。星载地球同步轨道合成孔径雷达(GEOSAR)具有覆盖区域广、重访时间短等优点。在一个GEOSAR的轨道周期中,地球中心转动(ECR)坐标系速度的方向变化范围大。当GEOSAR位于赤道附近时,ECR速度方向接近距离方向。此时,模糊函数投影在地面的二维旁瓣将不再是正交的。这导致方位向地距分辨率远大于方位向斜距分辨率,地距分辨率严重恶化。
合成孔径雷达可以通过提高发射带宽提高距离向分辨率,可以通过延长合成孔径时间提高方位向分辨率。但由于距离向、方位向两维不再正交,甚至可能重合,因而仅仅通过提高发射带宽与延长合成孔径时间难以有效的提高分辨率。
合成孔径雷达还可以通过姿态导引改变波束指向,进而改善地距分辨率。姿态导引是指为达到某种理想特性而进行的卫星姿态的调整。在地球低轨道(LEO)合成孔径雷达中,姿态导引的目的通常是降低多普勒频率。目前LEOSAR常用的姿态导引方法有一维偏航导引与全零多普勒导引,前者可降低多普勒频率,后者可使多普勒频率为零。在GEOSAR中,偏航导引能够有效降低地距分辨率的恶化,但所需的导引角大,对平台姿态控制***的要求高。全零多普勒导引不仅需要较大的导引角,而且不能降低地距分辨率的恶化。
发明内容
有鉴于此,本发明提供了一种星载地球同步轨道合成孔径雷达姿态导引方法,能够利用较小的姿态导引角实现最优地距分辨率。
本发明的星载地球同步轨道合成孔径雷达姿态导引方法,包括以下步骤:
步骤1,确定地距分辨椭圆面积:
其中,S是-3dB地距分辨率椭圆的面积,ρa为斜距方位向分辨率;ρr为斜距距离向分辨率;α为雷达的俯仰角;β为雷达的下视角;为雷达的方位角;
其中,
其中,ωs为轨道角速度,ωe为地球自转角速度,Rs为卫星距离地心的高度,为Rs对时间的导数,i为轨道倾角,u为纬度幅角,μ为地球引力常数,a为轨道半长轴,e为轨道偏心率,ω为近地点幅角;
步骤2,根据地距分辨率椭圆面积确定最优方位角:
求解得出最优方位角为
其中,正号表示右侧视,负号表示左侧视;
步骤3,根据最优方位角确定姿态导引角,进行姿态导引:
若采用先横滚后俯仰的引导方式,则横滚导引角θR,俯仰导引角θP分别为
若采用先俯仰后横滚的引导方式,则横滚导引角θR,俯仰导引角θP分别为
其中,
有益效果:
本发明利用地距分辨椭圆面积最小时地距分辨率最优,解算出最优地距分辨率的最优方位角,然后通过俯仰横滚导引或横滚俯仰导引使波束指向最优方位角所在方向,完成最优分辨率姿态导引,姿态导引角度小。
附图说明
图1为本发明流程图。
图2为本发明坐标角度示意图。
图3为采用先横滚后俯仰引导方式的姿态导引横滚角变化。
图4为采用先横滚后俯仰引导方式的姿态导引俯仰角变化。
图5为姿态导引后地距模糊函数。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
本发明提供了一种星载地球同步轨道合成孔径雷达姿态导引方法,首先根据已知的轨道参数、下视角等条件获得卫星速度方向与地距分辨椭圆面积的表达式。然后根据地距分辨椭圆面积的表达式,求解实现最小分辨椭圆面积所需的最优方位角。最后,通过俯仰横滚导引或横滚俯仰导引使波束指向上述最优方位角所在方向,如图1所示。具体步骤如下。
步骤1,根据轨道参数与下视角确定卫星速度方向与地距分辨椭圆面积。
步骤1.1,确定卫星速度方向。
根据轨道半长轴a、轨道倾角i、轨道偏心率e、近地点幅角ω、纬度幅角u确定卫星在地球固定坐标系中的速度方向。其中,地球固定坐标系原点位于地球质心,x轴在地球赤道面内指向格林尼治子午线,z轴沿地球自转方向指向北极,y轴与x轴、z轴构成右手直角坐标系。
地球固定坐标系下卫星速度的方向可使用偏航角φ与俯仰角α表示,偏航角φ与俯仰角α计算公式如下:
其中,ωe为地球自转角速度,ωs为轨道角速度,Rs为卫星距离地心的高度,为Rs对时间的导数,可以表示为
其中,μ是地球引力常数,数值为398600.5km3/s2;sign(·)为符号函数,可以表示为
步骤1.2,确定地距分辨椭圆面积
地距分辨率可以使用-3dB分辨率椭圆的面积表示,其公式为
其中,S是-3dB分辨率椭圆的面积,ρa为斜距方位向分辨率,ρr为斜距距离向分辨率。斜视角γ为波束指向与ECR速度方向的夹角,方位角为波束指向在水平面的投影与ECR速度在水平面的投影间的夹角,β为下视角,方位角和下视角β根据实际的波束指向确定。
上述角度如图2所示,其中V为ECR速度,R为波束指向,O点为地心。
斜视角γ与方位角之间存在如下关系:
因而,地距分辨率面积可以表示为
步骤2,根据地距分辨率椭圆面积确定最优方位角。
若地距分辨率面积最小,需要满足方程
求解方程(7),可获得对应于最优地距分辨率的最优方位角:
其中,是最优方位角,正号表示右侧视,负号表示左侧视。
步骤3,根据最优方位角确定姿态导引角。
在初始时刻波束指向星下点,姿态导引后波束指向最优方位角所在位置。姿态导引策略为横滚俯仰二维导引。姿态导引前的波束指向可以表示
u'=[0,0,-1]T(9)
姿态导引后的波束指向可以表示为
姿态导引过程可以表示为
u'=Au(11)
其中,A为姿态导引对应的旋转矩阵。
若采用先横滚后俯仰的策略,旋转矩阵可以表示为
其中,θR为横滚导引角,θP为俯仰导引角。将旋转矩阵代入方程(11),求解方程(11)可得到姿态导引角度为
若采用先俯仰后横滚的策略,旋转矩阵可以表示为
旋转矩阵代入方程(11),求解方程(11)可得到姿态导引角度为
自此,就实现了一种星载合成孔径雷达姿态导引方法。
卫星运行在椭圆轨道上,轨道的各参数具体如下所示:地球半径为6371.004km;轨道半长轴为42100km;轨道倾角为50°;轨道偏心率为0.1;轨道近地点幅角为90°;下视角为7°;斜距分辨率为1m;地球自转角速度为7.292115×10-5rads。姿态导引采用先横滚后俯仰的方式。姿态导引所需的横滚角如图3所示,俯仰角如图4所示。姿态引导后的地距分辨率模糊函数如图5所示。由图3、图4可知,地球同步轨道合成孔径雷达最优分辨率姿态导引方法所需的姿态导引角度较小。由图5可知,在进行地球同步轨道合成孔径雷达最优分辨率姿态导引后,地距分辨两位旁瓣垂直,地距分辨率最优。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种星载地球同步轨道合成孔径雷达姿态导引方法,其特征在于,包括以下步骤:
步骤1,确定地距分辨椭圆面积:
其中,S是-3dB地距分辨率椭圆的面积,ρa为斜距方位向分辨率;ρr为斜距距离向分辨率;α为雷达的俯仰角;β为雷达的下视角;为雷达的方位角;
其中,sign(·)为符号函数,可以表示为
其中,ωs为轨道角速度,ωe为地球自转角速度,Rs为卫星距离地心的高度,为Rs对时间的导数,i为轨道倾角,u为纬度幅角,μ为地球引力常数,a为轨道半长轴,e为轨道偏心率,ω为近地点幅角;
步骤2,根据地距分辨率椭圆面积确定最优方位角:
求解得出最优方位角为
其中,正号表示右侧视,负号表示左侧视;
步骤3,根据最优方位角确定姿态导引角,进行姿态导引:
若采用先横滚后俯仰的引导方式,则横滚导引角θR,俯仰导引角θP分别为
若采用先俯仰后横滚的引导方式,则横滚导引角θR,俯仰导引角θP分别为
其中,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310651459.3A CN103675760B (zh) | 2013-12-03 | 2013-12-03 | 一种星载地球同步轨道合成孔径雷达姿态导引方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310651459.3A CN103675760B (zh) | 2013-12-03 | 2013-12-03 | 一种星载地球同步轨道合成孔径雷达姿态导引方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103675760A CN103675760A (zh) | 2014-03-26 |
CN103675760B true CN103675760B (zh) | 2015-12-02 |
Family
ID=50313903
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310651459.3A Expired - Fee Related CN103675760B (zh) | 2013-12-03 | 2013-12-03 | 一种星载地球同步轨道合成孔径雷达姿态导引方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103675760B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109001762A (zh) * | 2018-06-08 | 2018-12-14 | 上海微小卫星工程中心 | 一种同步带目标光学观测卫星的地气光抑制方法及*** |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104597446B (zh) * | 2014-10-10 | 2017-02-15 | 北京理工大学 | 一种星载合成孔径雷达地距分辨率表征与参数设计方法 |
CN104375511A (zh) * | 2014-10-28 | 2015-02-25 | 上海卫星工程研究所 | 波束协同控制的地球同步轨道sar卫星偏航导引方法 |
CN104765023B (zh) * | 2015-03-20 | 2017-03-29 | 北京理工大学 | 一种星载合成孔径雷达高度分辨率计算方法 |
CN105403888A (zh) * | 2015-12-30 | 2016-03-16 | 中国科学院电子学研究所 | 采用波束指向控制的地球同步轨道圆迹sar成像方法 |
CN108958272A (zh) * | 2018-06-15 | 2018-12-07 | 上海卫星工程研究所 | 静止卫星成像偏航导引方法 |
CN111796277B (zh) * | 2020-02-28 | 2022-05-20 | 北京理工大学 | 一种基于无人机平台的穿墙雷达快速成像方法 |
CN112255606A (zh) * | 2020-09-29 | 2021-01-22 | 北京空间飞行器总体设计部 | 基于单反射面天线的Geo-SAR卫星正侧视成像姿态角计算方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE3114600C2 (zh) * | 1981-04-10 | 1990-06-21 | Telefunken Systemtechnik Gmbh, 7900 Ulm, De | |
EP1965223A1 (en) * | 2007-03-02 | 2008-09-03 | Saab Ab | Subsurface Imaging radar |
CN101464511B (zh) * | 2009-01-16 | 2011-07-27 | 北京航空航天大学 | 一种星载合成孔径雷达的工作波位判定方法 |
CN102608576A (zh) * | 2012-03-20 | 2012-07-25 | 北京理工大学 | 一种用于大前斜合成孔径雷达回波图像的几何校正方法 |
CN103364782A (zh) * | 2013-07-04 | 2013-10-23 | 北京理工大学 | 一种地球同步轨道合成孔径雷达速度空变性的补偿方法 |
-
2013
- 2013-12-03 CN CN201310651459.3A patent/CN103675760B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE3114600C2 (zh) * | 1981-04-10 | 1990-06-21 | Telefunken Systemtechnik Gmbh, 7900 Ulm, De | |
EP1965223A1 (en) * | 2007-03-02 | 2008-09-03 | Saab Ab | Subsurface Imaging radar |
CN101256235A (zh) * | 2007-03-02 | 2008-09-03 | Saab公司 | 表面下成像雷达 |
CN101464511B (zh) * | 2009-01-16 | 2011-07-27 | 北京航空航天大学 | 一种星载合成孔径雷达的工作波位判定方法 |
CN102608576A (zh) * | 2012-03-20 | 2012-07-25 | 北京理工大学 | 一种用于大前斜合成孔径雷达回波图像的几何校正方法 |
CN103364782A (zh) * | 2013-07-04 | 2013-10-23 | 北京理工大学 | 一种地球同步轨道合成孔径雷达速度空变性的补偿方法 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109001762A (zh) * | 2018-06-08 | 2018-12-14 | 上海微小卫星工程中心 | 一种同步带目标光学观测卫星的地气光抑制方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN103675760A (zh) | 2014-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103675760B (zh) | 一种星载地球同步轨道合成孔径雷达姿态导引方法 | |
CN113439220B (zh) | 一种基于卫星编队的遥感***及星座*** | |
CN102508243B (zh) | 一种倾斜地球同步轨道合成孔径雷达的波位设计方法 | |
US9927513B2 (en) | Method for determining the geographic coordinates of pixels in SAR images | |
CN107607947B (zh) | 基于卡尔曼滤波的星载雷达成像参数在线估计方法 | |
CN111965640B (zh) | 一种大斜视中高轨sar二维波束扫描方法 | |
CN102169173B (zh) | 一种倾斜地球同步轨道合成孔径雷达模糊度分析方法 | |
CN103364782B (zh) | 一种地球同步轨道合成孔径雷达速度空变性的补偿方法 | |
CN110488292B (zh) | 一种基于卫星编队的遥感*** | |
CN102621994B (zh) | 一种地球同步轨道合成孔径雷达覆盖中国全境的控制方法 | |
CN101270993A (zh) | 一种远程高精度自主组合导航定位方法 | |
CN101226058B (zh) | 一种卫星侧摆机动图像实时传输的实现方法 | |
CN103076607B (zh) | 一种基于sar卫星姿态控制实现滑动聚束模式的方法 | |
CN102866393B (zh) | 一种基于pos与dem数据的sar多普勒参数估计方法 | |
CN110456350B (zh) | 一种星载sar星座*** | |
CN110823191B (zh) | 混合基线双天线斜视干涉sar洋流测量性能确定方法及*** | |
US6677884B2 (en) | Satellite configuration for interferometric and/or tomographic remote sensing by means of synthetic aperture radar (SAR) | |
CN112255606A (zh) | 基于单反射面天线的Geo-SAR卫星正侧视成像姿态角计算方法 | |
CN104597446B (zh) | 一种星载合成孔径雷达地距分辨率表征与参数设计方法 | |
CN111856423A (zh) | 星载sar回波仿真处理方法、装置和设备 | |
CN103760562B (zh) | 一种机载圆迹合成孔径雷达航线的获取方法 | |
CN107064935A (zh) | 一种星载sar***及其构建方法 | |
CN103245948B (zh) | 双区成像合成孔径雷达图像匹配导航的方法 | |
CN112130147B (zh) | 一种基于海陆目标位置信息的成像波位确定方法 | |
Ze et al. | A new satellite attitude steering approach for zero Doppler centroid |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
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: 20151202 Termination date: 20161203 |
|
CF01 | Termination of patent right due to non-payment of annual fee |