CN113917418B - 一种用于斜视星载sar地平面二维分辨率评估方法 - Google Patents
一种用于斜视星载sar地平面二维分辨率评估方法 Download PDFInfo
- Publication number
- CN113917418B CN113917418B CN202111526812.6A CN202111526812A CN113917418B CN 113917418 B CN113917418 B CN 113917418B CN 202111526812 A CN202111526812 A CN 202111526812A CN 113917418 B CN113917418 B CN 113917418B
- Authority
- CN
- China
- Prior art keywords
- vector
- plane
- projection
- ground
- azimuth
- 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
- 238000000034 method Methods 0.000 title claims abstract description 31
- 208000004350 Strabismus Diseases 0.000 title claims description 18
- 239000013598 vector Substances 0.000 claims abstract description 125
- 238000003384 imaging method Methods 0.000 claims abstract description 47
- 230000000007 visual effect Effects 0.000 claims abstract description 10
- 238000006243 chemical reaction Methods 0.000 claims abstract description 5
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 238000005286 illumination Methods 0.000 claims 1
- 238000013461 design Methods 0.000 description 6
- 238000011156 evaluation Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 238000011161 development Methods 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000009466 transformation 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
- 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
- G01S7/40—Means for monitoring or calibrating
-
- 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
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
本发明公开了一种用于斜视星载SAR地平面二维分辨率评估方法,包括如下步骤:步骤1、根据卫星俯仰向视角和方位向视角确定成像时的地面瞄准点的三维坐标;步骤2、根据成像起止时刻的卫星轨道信息和地面瞄准点坐标,建立斜距平面的法向矢量;步骤3、在地面瞄准点处,根据地球椭球模型建立局部切平面,该切平面即为在常规SAR成像中斜距‑地距转化的地平面;步骤4、根据照射中心时间的卫星位置和瞄准点坐标,构建斜距矢量;步骤5、根据成像几何得到斜距矢量在地平面上的投影,将投影矢量与斜距矢量求幅度之比得距离向分辨率的投影因子;步骤6、根据斜距平面法向矢量和斜距矢量得到斜距平面的方位分辨单元矢量,计算得到方位向分辨率的投影因子。
Description
技术领域
本发明涉及雷达测量领域,尤其是一种星载合成孔径雷达 (Synthetic ApertureRadar, SAR)斜视成像几何下的地平面二维分辨率评估。
背景技术
随着星载SAR技术的不断发展,星载SAR图像的分辨率与成像幅宽不断提升,并在地球环境监测、军/民目标监测、分类与识别等领域发挥了重大作用。近年来,随着星载SAR天线技术的快速发展、星载SAR应用不断提升,星载SAR的成像模式从传统的波束固定模式发展到二维波束扫描。从传统的正侧视成像模式向斜视成像模式发展,基于星载SAR方位扫描能力的发展,将逐步演变出斜视条带、多角度聚束、多角度滑聚及马赛克等成像模式,卫星的观测能力将大幅提升。
然而,星载SAR在斜视成像模式下面临着许多新的问题,其中地平面的二维分辨率是斜视模式星载SAR系统设计的一个基础问题。在斜距成像平面,二维旁瓣在正侧视及斜视条件下均是正交的。然而将目标投影到地平面后,一方面,二维旁瓣的方向将不再正交,且旁瓣夹角随着斜视角的增加而增加;另一方面,在正侧视条件下,方位向分辨单元从斜距平面投影到地平面时分辨率不会损失。而在斜视条件下,方位向分辨单元从斜距平面投影到地平面时分辨率将会下降,即斜距平面到地平面的方位向投影因子大于1。对于距离向分辨单元,地平面视线方向分辨率会随着局部入射角的变化而变化。在二维扫描情况下,由斜距平面到地距平面投影的距离向分辨率投影因子和方位向分辨率投影因子将会是空变的。不同俯仰向视角和方位向扫描角下的距离和方位投影因子在星载SAR斜视成像模式的系统设计中具有重要的作用,对于指导星载SAR系统成像带宽设计、合成孔径时间设计等具有重要的指导意义。传统基于回波数据的评估方法和基于模糊函数的评估方法所需计算量大,无法在短时间内得到二维变化的距离向和方位向分辨率投影因子。由此,为提升系统设计效率,有必要研究星载SAR斜视成像模式下的地平面二维分辨率评估方法。
发明内容
有鉴于此,本发明的主要目的在于提供一种斜视星载SAR地平面二维分辨率评估方法,能够快速且精确的计算星载SAR在斜视成像模式下的距离向和方位向投影因子。本发明根据卫星俯仰向视角和方位向视角确定地面瞄准点的三维坐标;根据波束照射时间内的卫星轨道和瞄准点坐标建立斜距平面的法向矢量;在瞄准点处,根据地球椭球模型建立局部切平面(该切平面即为在常规SAR成像中斜距-地距转化的地平面);求解距离向投影因子:根据照射中心时间的卫星位置和瞄准点坐标,构建斜距矢量;根据成像几何得到斜距矢量在地平面上的投影,将投影矢量与斜距矢量求幅度之比得距离向分辨率的投影因子;求解方位投影因子:根据斜距平面法向矢量和斜距矢量得到斜距平面的方位分辨单元矢量,进而得到该矢量在地平面的投影矢量,将方位投影矢量与斜距平面内的方位分辨单元矢量求幅度之比得方位向分辨率的投影因子;对于不同的俯仰角和方位角可计算得到全部视角内的二维分辨率在地平面内的变化规律,从而为星载SAR的系统设计提供指导。
本发明的技术方案具体为:一种用于斜视星载SAR地平面二维分辨率评估方法,包括如下步骤:
步骤1、根据卫星俯仰向视角和方位向视角确定成像时的地面瞄准点的三维坐标;
步骤2、根据成像起止时刻的卫星轨道信息和地面瞄准点坐标,建立斜距平面的法向矢量;
步骤3、在地面瞄准点处,根据地球椭球模型建立局部切平面,该切平面即为在常规SAR成像中斜距-地距转化的地平面;
步骤4、根据照射中心时间的卫星位置和瞄准点坐标,构建斜距矢量;
步骤5、根据成像几何得到斜距矢量在地平面上的投影,将投影矢量与斜距矢量求幅度之比得距离向分辨率的投影因子;
步骤6、根据斜距平面法向矢量和斜距矢量得到斜距平面的方位分辨单元矢量,进而得到该方位分辨单元矢量在地平面的投影矢量,即方位投影矢量,将方位投影矢量与斜距平面内的方位分辨单元矢量求幅度之比得方位向分辨率的投影因子。
有益效果:
对于星载SAR分辨率的评估问题,不同的成像模式、场景大小等因素都会导致分辨率评估方法的复杂性。常规的分辨率评估方法包括基于回波仿真并二维聚焦的方法、基于二维模糊函数的方法。而基于回波仿真的方法需要模拟雷达回波数据并进行成像处理,运算量大、时间复杂度高,且无法获得整个场景的地距分辨率变化。其次,基于二维模糊函数的方法虽能获取整个场景的地距分辨率,但其时间消耗大,在实际工程参数优化过程中,需要大量的时间。本方法基于几何成像关系计算二维地距分辨率,能在短时间内获得整个成像场景的二维分辨率变化,效率与准确度高。
附图说明
图1为斜视星载SAR地平面二维分辨率评估流程框图;
图2为星载SAR成像几何模型;
图3为距离分辨率投影因子计算几何;
图4为一组实例斜视星载SAR方位向分辨率投影因子;
图5为一组实例斜视星载SAR距离向分辨率投影因子。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整的描述,显然,所描述的实施例仅为本发明的一部分实施例,而不是全部的实施例,基于本发明中的实施例,本领域的普通技术人员在不付出创造性劳动的前提下所获得的所有其他实施例,都属于本发明的保护范围。
根据本发明的实施例,如图1所示,为一种用于斜视星载SAR地平面二维分辨率评估方法的流程图,包括如下步骤:
步骤1、根据卫星俯仰向视角和方位向视角确定成像时的地面瞄准点的三维坐标;
步骤2、根据成像起止时刻的卫星轨道信息和地面瞄准点坐标,建立斜距平面的法向矢量;
步骤3、在地面瞄准点处,根据地球椭球模型建立局部切平面,该切平面即为在常规SAR成像中斜距-地距转化的地平面;
步骤4、根据照射中心时间的卫星位置和瞄准点坐标,构建斜距矢量;
步骤5、根据成像几何得到斜距矢量在地平面上的投影,将投影矢量与斜距矢量求幅度之比得距离向分辨率的投影因子;
步骤6、根据斜距平面法向矢量和斜距矢量得到斜距平面的方位分辨单元矢量,进而得到该方位分辨单元矢量在地平面的投影矢量,即方位投影矢量,将方位投影矢量与斜距平面内的方位分辨单元矢量求幅度之比得方位向分辨率的投影因子。
具体的,所述步骤1、根据轨道信息和卫星的俯仰、方位扫描角得到成像时的地面瞄准点坐标。具体如下:根据卫星的轨道历程(其中N a 表示N a 个脉冲发射时刻)或轨道根数,Xs、Ys、Zs为卫星三维位置坐标,星载SAR平台视角θ、方位向扫描角φ等,采用零多普勒姿态导引方法,得到对应于成像时间的地面瞄准点,记为,、、为地面瞄准点位置三维位置坐标。
所述步骤2中,根据成像起止时刻的卫星轨道信息和地面瞄准点坐标,建立斜距平面的法向矢量;具体为:
其中‘×’表示矩阵叉乘。
所述步骤3、在地面瞄准点处,根据地球椭球模型建立局部切平面具体为:
步骤5、根据成像几何得到斜距矢量在地平面上的投影,将投影矢量与斜距矢量求幅度之比得距离向分辨率的投影因子;如图3所示为距离向分辨率由斜距平面到地平面的投影几何。设斜距矢量在地平面的投影为,则该投影可通过如下矢量计算得到:
其中,‘· ’表示矢量点乘。从而可解算常数k为:
所述步骤6、根据斜距平面法向矢量和斜距矢量得到斜距平面的方位分辨单元矢量,进而得到该方位分辨单元矢量在地平面的投影矢量,即方位投影矢量,将方位投影矢量与斜距平面内的方位分辨单元矢量求幅度之比得方位向分辨率的投影因子。具体为:
从而可得到公式(9)中的常数系数为:
从而,从斜距平面到地平面转化过程中,方位分辨单元的投影因子为:
实施例1
本处采用一个典型的星载SAR斜视成像模式验证本发明的技术方案的有效性。
如图2所示,卫星沿轨道ACB运行,其中C点为轨迹中心点,x-y-z组成场景的地距平面坐标系,x’-y’-z’组成斜平面坐标系。
图3为距离分辨率投影因子计算几何;如图3所示,在斜距平面中,方位向分辨单元矢量和距离向分辨单元矢量相互垂直,投影矢量的方向为斜平面的法向方向,将上述两个分辨单元矢量投影到地距平面坐标系,此时,两个方向的分辨单元大小将发生变化。
图4为一组实例斜视星载SAR方位向分辨率投影因子;
图5为一组实例斜视星载SAR距离向分辨率投影因子。图4和图5所示的两个方向的分辨率投影因子表示了二维分辨率从斜距平面到地距平面的变换系数。以某一个方位扫描角和视角为例,假设斜距平面的方位向和距离向分辨率分别为ρ a 、ρ b ,投影因子分别为a和b,那么对应的地距平面分辨率则为a·ρ a ,b·ρ b ,在系统迭代优化的过程中,加快了二维分辨率计算效率。
以上所述,仅为本发明的一个实施例而已,并非用于限定本发明的保护范围。
尽管上面对本发明说明性的具体实施方式进行了描述,以便于本技术领域的技术人员理解本发明,且应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
Claims (7)
1.一种用于斜视星载SAR地平面二维分辨率评估方法,其特征在于,包括如下步骤:
步骤1、根据卫星俯仰向视角和方位向视角确定成像时的地面瞄准点的三维坐标;
步骤2、根据成像起止时刻的卫星轨道信息和地面瞄准点坐标,建立斜距平面的法向矢量;
步骤3、在地面瞄准点处,根据地球椭球模型建立局部切平面,该切平面即为在常规SAR成像中斜距-地距转化的地平面;
步骤4、根据照射中心时间的卫星位置和瞄准点坐标,构建斜距矢量;
步骤5、根据成像几何得到斜距矢量在地平面上的投影,将投影矢量与斜距矢量求幅度之比得距离向分辨率的投影因子;
步骤6、根据斜距平面法向矢量和斜距矢量得到斜距平面的方位分辨单元矢量,进而得到该方位分辨单元矢量在地平面的投影矢量,即方位投影矢量,将方位投影矢量与斜距平面内的方位分辨单元矢量求幅度之比得方位向分辨率的投影因子。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111526812.6A CN113917418B (zh) | 2021-12-15 | 2021-12-15 | 一种用于斜视星载sar地平面二维分辨率评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111526812.6A CN113917418B (zh) | 2021-12-15 | 2021-12-15 | 一种用于斜视星载sar地平面二维分辨率评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113917418A CN113917418A (zh) | 2022-01-11 |
CN113917418B true CN113917418B (zh) | 2022-02-22 |
Family
ID=79249204
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111526812.6A Active CN113917418B (zh) | 2021-12-15 | 2021-12-15 | 一种用于斜视星载sar地平面二维分辨率评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113917418B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116502476B (zh) * | 2023-06-28 | 2023-09-08 | 中国科学院空天信息创新研究院 | 一种基于非线性频率扫描响应的sar系统设计方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012020439A1 (en) * | 2010-08-11 | 2012-02-16 | Selex Galileo S.P.A. | Multi-grazing isar imaging method and system providing isar side-view images with improved cross-range resolution |
CN104597446A (zh) * | 2014-10-10 | 2015-05-06 | 北京理工大学 | 一种星载合成孔径雷达地距分辨率表征与参数设计方法 |
CN104765023A (zh) * | 2015-03-20 | 2015-07-08 | 北京理工大学 | 一种星载合成孔径雷达高度分辨率计算方法 |
CN111127319A (zh) * | 2019-12-09 | 2020-05-08 | 航天东方红卫星有限公司 | 一种动中推扫成像的地面像元分辨率计算方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102009030672B3 (de) * | 2009-06-25 | 2010-08-19 | Eads Deutschland Gmbh | Verfahren zur Bestimmung der geographischen Koordinaten von Bildpunkten in SAR Bildern |
-
2021
- 2021-12-15 CN CN202111526812.6A patent/CN113917418B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012020439A1 (en) * | 2010-08-11 | 2012-02-16 | Selex Galileo S.P.A. | Multi-grazing isar imaging method and system providing isar side-view images with improved cross-range resolution |
CN104597446A (zh) * | 2014-10-10 | 2015-05-06 | 北京理工大学 | 一种星载合成孔径雷达地距分辨率表征与参数设计方法 |
CN104765023A (zh) * | 2015-03-20 | 2015-07-08 | 北京理工大学 | 一种星载合成孔径雷达高度分辨率计算方法 |
CN111127319A (zh) * | 2019-12-09 | 2020-05-08 | 航天东方红卫星有限公司 | 一种动中推扫成像的地面像元分辨率计算方法 |
Non-Patent Citations (3)
Title |
---|
Generalized approach to resolution analysis in BSAR;Tao Zeng 等;《IEEE Transactions on Aerospace and Electronic Systems》;20050430;第41卷(第2期);全文 * |
The Accurate Focusing and Resolution Analysis Method in Geosynchronous SAR;Cheng Hu 等;《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》;20111031;第49卷(第10期);全文 * |
雷达下视环扫成像分辨率研究;李静 等;《系统工程与电子技术》;20090228;第31卷(第2期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113917418A (zh) | 2022-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9927513B2 (en) | Method for determining the geographic coordinates of pixels in SAR images | |
CN107367716A (zh) | 一种高精度星载sar几何定标方法 | |
CN102654576B (zh) | 基于sar图像和dem数据的图像配准方法 | |
EP2442275A2 (en) | Method and apparatus for three-dimensional image reconstruction | |
CN102565797A (zh) | 一种针对聚束模式星载sar图像的几何校正方法 | |
CN107300699B (zh) | 基于敏捷合成孔径雷达卫星姿态机动的马赛克模式实现方法 | |
CN105677942A (zh) | 一种重复轨道星载自然场景sar复图像数据快速仿真方法 | |
CN115792907B (zh) | 星载sar斜视滑动聚束模式方位向成像参数设计方法 | |
CN110673143A (zh) | 一种子孔径大斜视sar俯冲成像的两步处理方法 | |
KR102151362B1 (ko) | 극좌표변환을 이용한 항공기기반 영상복원장치 및 이를 이용한 영상복원방법 | |
CN111102981A (zh) | 一种基于ukf的高精度卫星相对导航方法 | |
CN113917418B (zh) | 一种用于斜视星载sar地平面二维分辨率评估方法 | |
CN104807477A (zh) | 一种基于靶标控制点的卫星ccd阵列影像几何检校方法 | |
CN110109104B (zh) | 一种阵列sar等距离切片成像几何畸变校正方法 | |
CN113409400A (zh) | 一种基于自动跟踪的机载光电系统目标地理定位方法 | |
CN115980873A (zh) | 一种基于光学观测的快速精密轨道确定方法及系统 | |
CN114660606A (zh) | 低信噪比isar图像序列匹配搜索的空间目标姿态反演方法 | |
CN106353756A (zh) | 基于图像匹配的下降轨聚束合成孔径雷达定位方法 | |
CN112179314B (zh) | 一种基于三维网格投影的多角度sar高程测量方法及系统 | |
KR102028323B1 (ko) | 영상 레이더의 영상 보정 장치 및 시스템 | |
CN115932823B (zh) | 基于异源区域特征匹配的飞行器对地面目标定位方法 | |
CN107907881B (zh) | 一种长孔径星载sar地形估计方法 | |
Shao et al. | Model-data co-driven integration of detection and imaging for geosynchronous targets with wideband radar | |
US7142149B2 (en) | Mensuration for the conformal range migration algorithm | |
CN114047486A (zh) | 一种雷达导引头挂飞试验安装误差角标定方法及存储介质 |
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 |