CN102928829A - 空间目标参数反演方法 - Google Patents
空间目标参数反演方法 Download PDFInfo
- Publication number
- CN102928829A CN102928829A CN201210510686XA CN201210510686A CN102928829A CN 102928829 A CN102928829 A CN 102928829A CN 201210510686X A CN201210510686X A CN 201210510686XA CN 201210510686 A CN201210510686 A CN 201210510686A CN 102928829 A CN102928829 A CN 102928829A
- Authority
- CN
- China
- Prior art keywords
- target
- rcs
- angle
- extraterrestrial target
- peak value
- 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.)
- Granted
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开一种空间目标参数反演方法,该方法包含以下步骤:1、根据空间目标及探测器的空间轨道位置计算获取观测角大小的变化曲线;2、通过空间微波雷达测量或仿真获取空间目标飞行过程中的动态雷达散射截面曲线,并根据镜面反射机制寻找雷达散射截面峰值;3、利用雷达散射截面峰值对应的观测角确定空间目标姿态;4、根据雷达散射截面从峰值的下降特性反演空间目标的横向尺寸。本发明直接利用空间目标近距离观测的窄带雷达散射截面RCS变化特性数据,实现其观测面内的姿态角及横向尺寸等参数的反演;利用微波雷达观测数据,反演获取目标的姿态参数,可为空间会合等提供参考判断依据;可为空间目标普查中建立特征数据库提供有效手段。
Description
技术领域
本发明涉及一种目标电磁散射特性数据处理技术,具体涉及一种空间目标参数反演方法。
背景技术
目前,根据微波散射特性的目标参数反演主要采用成像方式开展,如利用高分辨距离像反演目标径向尺寸参数,利用二维或三维散射图像反演目标姿态参数等。但基于成像的反演方式要求提供目标一定角度范围内的宽带RCS特性,现有空间微波雷达的观测条件难以满足要求。
而在探测器与空间目标交会、对接等过程,探测器与目标距离相对较近,且观测角变化缓慢,目标RCS变化也将较慢,利于研究其小角度范围内的变化特性。同时,卫星等空间目标的形状相对规则,如本体通常为方形或圆柱形,太阳能电池帆板通常为板状,这些外形结构的雷达散射截面(RCS,Radar Cross Section)变化特性具有相似性,可建立其随尺寸、姿态等参数的关系模型,进而进行空间目标参数的反演。目前,未见有利用空间目标小角度RCS变化特性进行参数反演研究的公开技术。
发明内容
本发明提供一种空间目标参数反演方法,根据空间目标外形特征,建立了其雷达散射截面峰值附近小角度范围内的变化特性模型,给出利用动态雷达散射截面峰值位置反演目标姿态、利用雷达散射截面下降特性反演目标尺寸的方法。
为实现上述目的,本发明提供一种空间目标参数反演方法,其特点是,该方法包含以下步骤:
步骤1、根据空间目标及探测器的空间轨道位置计算获取观测角大小的变化曲线;
步骤2、通过空间微波雷达测量或仿真获取空间目标飞行过程中的动态雷达散射截面曲线,并根据镜面反射机制寻找雷达散射截面峰值;
利用微波探测器测量或仿真数据,获取空间目标随轨道位置变化的动态雷达散射截面曲线,散射机制对应的小角度雷达散射截面的变化趋势为:
(2)
步骤3、利用雷达散射截面峰值对应的观测角确定空间目标姿态;
步骤4、根据雷达散射截面从峰值的下降特性反演空间目标的横向尺寸;
先选取雷达散射截面观测曲线的一个峰值;
再查找雷达散射截面曲线下降10dB时的观测偏离角度;
目标尺寸为:
上述的空间目标为平板型目标或圆柱形目标。
上述空间目标为平板形目标,在其镜面反射雷达散射截面峰值附近,其小角度雷达散射截面变化趋势为:
(3)
上述空间目标为圆柱形目标,在其镜面反射RCS峰值附近,其小角度雷达散射截面变化趋势为:
(6)。
本发明基于小角度雷达散射截面变化特性的空间目标参数反演方法和现有技术的目标参数反演相比,其优点在于,本发明直接利用空间目标近距离观测的窄带雷达散射截面RCS变化特性数据,实现其观测面内的姿态角及横向尺寸等参数的反演;
本发明方法利用微波雷达观测数据,反演获取目标的姿态参数,可为空间会合等提供参考判断依据;
本发明中的反演方法,可为空间目标普查中建立特征数据库提供有效手段;
本发明方法反演的空间目标特征参数,还可为空间目标的识别提供依据。
附图说明
图1为本发明空间目标参数反演方法的流程图;
图2为本发明空间目标参数反演方法的观测角示意图;
图3为本发明空间目标参数反演方法的空间目标的观测角变化曲线图;
图4为本发明空间目标参数反演方法的空间目标小角度内的雷达散射截面变化趋势图;
图5为本发明空间目标参数反演方法的空间目标的雷达散射截面变化曲线图。
具体实施方式
以下结合附图,进一步说明本发明的具体实施例。
本发明公开一种基于小角度雷达散射截面(RCS,Radar Cross Section)变化特性的空间目标参数反演方法的实施例。本实施例中以某卫星目标为例,该卫星目标本体尺寸为1.2米,并假定该卫星目标轨道为圆形,探测器轨道为扁椭球形。
如图1所示,该空间目标参数反演方法包含以下步骤:
观测角的大小一般可由空间目标及探测器的空间轨道位置计算。如图2所示,这里观测角的定义选取空间目标的轨道局部坐标系,该坐标系X轴取空间目标的飞行方向,Y轴取空间目标指向地球中心的方向,观测角则取探测器至目标连线与X轴夹角,根据飞行轨道位置的不同,观测角将在一定范围内变化,这里记为。
根据两者轨道位置,利用公式(1)可计算获取观测角的变化曲线,结果如图3所示,其中横坐标为空间目标沿轨道飞行一周的角度位置。
同时,将卫星目标姿态角定义为卫星目标的轴向与图2中X轴夹角,这里记为;根据一般空间目标的形状特征,目标轴向一般与其大面垂直,根据镜面反射机制,目标轴向的后向RCS将存在局部峰值。
步骤2、空间目标动态雷达散射截面(RCS)获取及小角度特性建模,通过微波探测器,例如空间微波雷达,测量或仿真数据获取目标飞行过程中随轨道位置变化的动态雷达散射截面RCS曲线,并根据镜面反射机制寻找雷达散射截面峰值。
同时,由于空间中空间目标一般存在较多的镜面反射型强散射源,例如本体面、太阳帆板面等,这类散射机制对应的小角度RCS的一般变化趋势如图4所示,其中纵坐标为RCS下降量,横坐标为观测角偏离大小,该关系式通常可写为:
空间目标通常分为平板型目标和圆柱形目标。
1)对平板形目标,在其镜面反射RCS峰值附近,其小角度RCS变化趋势为:
(3)
2)对圆柱形目标,在其镜面反射RCS峰值附近,其小角度RCS变化趋势为:
如图5所示,本实施例中采用仿真数据获取目标飞行过程中随轨道位置变化的动态雷达散射截面RCS曲线,其中横坐标为空间目标的空间位置,纵坐标为雷达散射截面。假定观测电磁波频率为14GHz,这里采用仿真方式获取目标沿轨道飞行一周的RCS变化曲线,结果其中存在RCS峰值。
步骤3、相对姿态反演,利用雷达散射截面峰值对应的观测角确定空间目标姿态。
如图5所示,本实施例中,读取RCS峰值,约为48.9dBm2,其目标位置分别处于为13.9°及166.2°处;对应的观测角大小为-0.021°及-0.023°,即目标姿态角近似为-0.02°,而实际目标在该平面内姿态角为0°,两者十分接近。
步骤5、空间目标尺寸反演,根据雷达散射截面从峰值的下降特性反演空间目标的横向尺寸。
由于空间目标尺寸与镜面反射贡献的大小相关联,可根据RCS下降一定量时观测角变化直接反演目标尺寸。一般可先获取RCS观测曲线的峰值,再查找RCS曲线下降10dB时的观测偏离角度,目标尺寸则可近似估算为:
(7)
本实施例中,选取图5中166.2°处RCS峰进行估算,其RCS下降10dB的目标位置角约为164.8°,其对应的观测角变化为0.395°,根据公式(7),可计算目标尺寸约为1.14m,与实际尺寸1.2m相比,相对偏差约为5%,说明了本方法的有效性。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
Claims (7)
1.一种空间目标参数反演方法,其特征在于,该方法包含以下步骤:
步骤1、根据空间目标及探测器的空间轨道位置计算获取观测角大小的变化曲线;
步骤2、通过空间微波雷达测量或仿真获取空间目标飞行过程中的动态雷达散射截面曲线,并根据镜面反射机制寻找雷达散射截面峰值;
步骤3、利用雷达散射截面峰值对应的观测角确定空间目标姿态;
步骤4、根据雷达散射截面从峰值的下降特性反演空间目标的横向尺寸。
3. 如权利要求2所述的空间目标参数反演方法,其特征在于,所述的空间目标为平板型目标或圆柱形目标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210510686.XA CN102928829B (zh) | 2012-12-04 | 2012-12-04 | 空间目标参数反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210510686.XA CN102928829B (zh) | 2012-12-04 | 2012-12-04 | 空间目标参数反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102928829A true CN102928829A (zh) | 2013-02-13 |
CN102928829B CN102928829B (zh) | 2014-08-27 |
Family
ID=47643672
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210510686.XA Expired - Fee Related CN102928829B (zh) | 2012-12-04 | 2012-12-04 | 空间目标参数反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102928829B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103644916A (zh) * | 2013-11-28 | 2014-03-19 | 江西洪都航空工业集团有限责任公司 | 一种隐形飞机的航路规划方法 |
CN105353376A (zh) * | 2015-11-24 | 2016-02-24 | 北京环境特性研究所 | 一种基于isar图像的空间目标姿态反演方法 |
CN105928448A (zh) * | 2016-04-18 | 2016-09-07 | 北京理工大学 | 一种基于瑞利区解析散射建模的昆虫尺寸测量方法 |
CN107192990A (zh) * | 2017-06-20 | 2017-09-22 | 中国计量科学研究院 | 外推法测雷达散射截面积 |
CN111208485A (zh) * | 2020-01-16 | 2020-05-29 | 中国人民解放军63921部队 | 一种基于分段多项式拟合的空间目标rcs特性识别方法 |
CN113063991A (zh) * | 2021-03-15 | 2021-07-02 | 北京环境特性研究所 | 一种紧缩场极化偏转角的检测方法及装置 |
CN118501872A (zh) * | 2024-07-18 | 2024-08-16 | 中国人民解放军国防科技大学 | 基于太赫兹雷达的空间目标异动反演方法、装置及设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3378200B2 (ja) * | 1998-09-03 | 2003-02-17 | 三菱電機株式会社 | Rcs計測装置及びrcs計測方法 |
JP3668928B2 (ja) * | 2000-03-21 | 2005-07-06 | 三菱電機株式会社 | レーダ装置 |
CN102062857A (zh) * | 2010-12-10 | 2011-05-18 | 中国舰船研究设计中心 | 矩形平板近场雷达散射截面“平顶”效应预估方法 |
CN102401893A (zh) * | 2010-09-10 | 2012-04-04 | 上海无线电设备研究所 | 一种近场目标对称极化rcs测试方法 |
-
2012
- 2012-12-04 CN CN201210510686.XA patent/CN102928829B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3378200B2 (ja) * | 1998-09-03 | 2003-02-17 | 三菱電機株式会社 | Rcs計測装置及びrcs計測方法 |
JP3668928B2 (ja) * | 2000-03-21 | 2005-07-06 | 三菱電機株式会社 | レーダ装置 |
CN102401893A (zh) * | 2010-09-10 | 2012-04-04 | 上海无线电设备研究所 | 一种近场目标对称极化rcs测试方法 |
CN102062857A (zh) * | 2010-12-10 | 2011-05-18 | 中国舰船研究设计中心 | 矩形平板近场雷达散射截面“平顶”效应预估方法 |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103644916B (zh) * | 2013-11-28 | 2016-01-06 | 江西洪都航空工业集团有限责任公司 | 一种隐形飞机的航路规划方法 |
CN103644916A (zh) * | 2013-11-28 | 2014-03-19 | 江西洪都航空工业集团有限责任公司 | 一种隐形飞机的航路规划方法 |
CN105353376A (zh) * | 2015-11-24 | 2016-02-24 | 北京环境特性研究所 | 一种基于isar图像的空间目标姿态反演方法 |
CN105928448B (zh) * | 2016-04-18 | 2018-06-22 | 北京理工大学 | 一种基于瑞利区解析散射建模的昆虫尺寸测量方法 |
CN105928448A (zh) * | 2016-04-18 | 2016-09-07 | 北京理工大学 | 一种基于瑞利区解析散射建模的昆虫尺寸测量方法 |
CN107192990B (zh) * | 2017-06-20 | 2019-09-17 | 中国计量科学研究院 | 外推法测雷达散射截面积 |
CN107192990A (zh) * | 2017-06-20 | 2017-09-22 | 中国计量科学研究院 | 外推法测雷达散射截面积 |
CN111208485A (zh) * | 2020-01-16 | 2020-05-29 | 中国人民解放军63921部队 | 一种基于分段多项式拟合的空间目标rcs特性识别方法 |
CN111208485B (zh) * | 2020-01-16 | 2023-04-07 | 中国人民解放军63921部队 | 一种基于分段多项式拟合的空间目标rcs特性识别方法 |
CN113063991A (zh) * | 2021-03-15 | 2021-07-02 | 北京环境特性研究所 | 一种紧缩场极化偏转角的检测方法及装置 |
CN113063991B (zh) * | 2021-03-15 | 2022-10-21 | 北京环境特性研究所 | 一种紧缩场极化偏转角的检测方法及装置 |
CN118501872A (zh) * | 2024-07-18 | 2024-08-16 | 中国人民解放军国防科技大学 | 基于太赫兹雷达的空间目标异动反演方法、装置及设备 |
CN118501872B (zh) * | 2024-07-18 | 2024-09-20 | 中国人民解放军国防科技大学 | 基于太赫兹雷达的空间目标异动反演方法、装置及设备 |
Also Published As
Publication number | Publication date |
---|---|
CN102928829B (zh) | 2014-08-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102928829B (zh) | 空间目标参数反演方法 | |
Goddi et al. | Infall and outflow within 400 au from a high-mass protostar-3d velocity fields from methanol and water masers in aflg 5142 | |
US9495596B2 (en) | Systems and method for identifying locations of infrastructure assets using aerial imagery | |
CN101539627B (zh) | 一种电离层立体探测星载sar成像处理平台的构建方法 | |
Tahani et al. | 3D magnetic-field morphology of the Perseus molecular cloud | |
CN103744081B (zh) | 一种机载圆迹合成孔径雷达高精度三维成像补偿方法 | |
EP3789797B1 (en) | Method and apparatus for evaluating data, device, and computer-readable storage medium | |
Kusk et al. | Synthetic SAR image generation using sensor, terrain and target models | |
CN103630900A (zh) | 3-d sar波数域快速成像的方法 | |
CN110456304A (zh) | 机载测向定位方法 | |
Long et al. | Object detection research of SAR image using improved faster region-based convolutional neural network | |
CN103235304A (zh) | 一种基于误差等效rd模型修正的sar几何校正方法 | |
Liu et al. | Observational Scan-induced Artificial Cosmic Microwave Background Anisotropy | |
Wang et al. | High-precision pose measurement method based on binocular vision in dark lighting environments | |
CN107300692B (zh) | 星载sar目标定位精度指标地面评估方法 | |
Chen et al. | A new pose estimation method for non-cooperative spacecraft based on point cloud | |
Rabinovitch et al. | Full-scale supersonic parachute shape reconstruction using three-dimensional stereo imagery | |
CN103631990A (zh) | Sar照射区域的仿真场景模型建立方法和系统 | |
Zhuang et al. | Image processing and deep normalized CNN for the location measurement and reading distance prediction of RFID multi-tags | |
Chao-ying et al. | Attitude estimation method of space targets by 3D reconstruction of principal axis from ISAR image | |
LeGrand et al. | Relative multiple space object tracking using intensity filters | |
Gong et al. | Deformation measuring methods based on inertial sensors for airborne distributed POS | |
Huang et al. | Dynamic Object-aware LiDAR Odometry Aided by Joint Weightings Estimation in Urban Areas | |
Sun et al. | Simulated Gravity Field Estimation for the Main Belt Comet 133P/Elst-Pizarro Based on a Satellite-to-satellite Tracking Mode | |
Zhuang et al. | A novel simulation approach of aircraft dynamic RCS |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140827 Termination date: 20211204 |