CN106154268A - 一种基于子带分割的低波段isar电离层校正方法 - Google Patents
一种基于子带分割的低波段isar电离层校正方法 Download PDFInfo
- Publication number
- CN106154268A CN106154268A CN201610804971.0A CN201610804971A CN106154268A CN 106154268 A CN106154268 A CN 106154268A CN 201610804971 A CN201610804971 A CN 201610804971A CN 106154268 A CN106154268 A CN 106154268A
- Authority
- CN
- China
- Prior art keywords
- distance
- signal
- ionospheric
- echo
- orientation
- 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.)
- Withdrawn
Links
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
-
- 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/9064—Inverse SAR [ISAR]
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
本发明公开了一种基于子带分割的低波段ISAR电离层校正方法,对距离脉压和距离走动校正后的数据,利用距离向信号带宽和方位合成孔径时间这二个维度,将回波数据划分成多个子数据块。利用电离层干扰和信号斜距之间的关系,利用估计拟合计算得到的残留斜距误差推导出对应的电离层参数,从而建立电离层干扰校正函数。将回波数据与校正函数相乘,恢复出屋电离层干扰下的回波信号。本发明在有效减弱电离层对雷达回波信号的影响的同时进一步提高ISAR的成像质量,加强了地基低波段ISAR雷达的空间目标监测的有效性。
Description
技术领域
本发明涉及ISAR成像信号补偿方法领域,具体是一种基于子带分割的低波段ISAR电离层误差校正方法。
背景技术
几乎所有的SAR/ISAR成像算法都假设雷达或者目标在真空的理想路径中运动,雷达或者目标的运动位置和运动状态在观测时间内都精确已知。然而在雷达信号发射、传播、接收、处理的各个过程里,存在方方面面的误差干扰影响最终获得的雷达图像的效果。
雷达发射的电磁波在空间传播中,受到地球大气、地球电离层的干扰,使得成像质量发生图像偏移、几何畸变和分辨率降低等不利影响。早在20世纪中期,美国NASA、欧空局等就已发布了包含电离层对雷达性能影响的研究项目。Quegan等(1986)利用合成天线增益方向图和相位屏模拟研究了电离层对SAR的影响。而后,电离层对200MHz-2GHz频率星载SAR、电离层色散对距离和方位向分辨率、电离层吸收对信号能量、电离层闪烁对雷达性能的影响,以及Faraday旋转对极化SAR和干涉SAR成像的影响也迅速开展。
自二十世纪八九十年代开始,低频段、大带宽成为雷达***发展的一个重要方向,低频段信号(VHF/UHF/L/P)的穿透特性在军事和民用领域具有重要应用前景。但低频段信号受电离层的影响十分明显,回波信号幅度被削弱、相位严重恶化,从而难以形成高质量图像。国内外学者对恢复受电离层影响的电磁波成像已进行了大量研究,但对P波段的实测数据进行的校正研究目前较少。Jun Liu等(2003)针对甚高频段(Ultra High Frequency,UHF)信号分析了电离层对SAR雷达信号分辨率和成像位置的影响。Arlen D Schmidt(2000)和唐辉(2003)分别研究了电离层对超高频段(Very High Frequency,VHF)波段和L波段的地基雷达对空间目标ISAR成像时的电离层影响,并提出了校正方法。翟文帅(2008)等给出了一种通过测量P波段回波时宽展宽量的校正方法,而这种方法需要微秒级的测量精度,在存在杂波时难以实际应用。
发明内容
本发明的目的是提供一种基于子带分割的低波段ISAR电离层校正方法,以解决现有技术在存在电离层干扰时不能正常成像的问题,提出的校正方法能够有效去除信号中的电离层误差,还原出无干扰数据,提高雷达图像的成像效果。
为了达到上述目的,本发明所采用的技术方案为:
(1)、二维子带数据分割,具体过程如下:
1a)、按信号带宽划分:
将距离脉压和距离走动校正后的数据转换到距离频域-方位时域,将距离向信号带宽划分成M个距离子带,每个子带的中心频率(fm)互不相同;
1b)、按方位合成孔径时间划分:
按方位合成孔径时间将回波数据继续划分成N个方位子孔径,每个方位子孔径对应的数据获取时间(n)互不相同;之内可以看作电离层参数TEC近似不变;
(2)、电离层干扰误差计算,具体过程如下:
电离层对回波信号产生的干扰误差体现在回波的相位对应的斜距中,得到斜距就可以反推出电离层干扰的大小。对每一个二维子带划分后的块数据,进行回波脉压处理,得到脉压峰值位置,由峰值位置的平均值与距离分辨单元大小的乘积得到该块数据对应的残留斜距误差
(3)、电离层误差校正,信号恢复,具体步骤如下:
3a)、误差向量拟合:
对于同一个方位子孔径时间内的M个距离子带中估计得到的M个残留斜距误差记成误差向量:
对误差向量在中心载频(fc)处进行多项式拟合,结果可以写成:
其中,f=fc+fr,fc是信号的中心载频,fr是距离频率且fr∈[-B/2,B/2],B是信号带宽,p0、p1和p2分别是各次多项式的系数;
3b)、电离层参数计算:
在研究某个方位慢时间对应的电离层影响时,斜距误差和电离层参数TEC之间的关系可近似推导写为:
其中,K=40.28m3/s2为电离层常数,tm为慢时间,c为真空中的光速,上式中距离频率fr对应项的系数与拟合多项式的系数应一一对应相等,利用fr的一次项系数可得每个方位子孔径对应的电离层参数:
3c)、无干扰信号恢复:
利用每个方位子孔径对应的电离层参数,得到电离层参数TEC在整个方位观测时间内的变化,建立推导公式中的各项相位的共轭补偿函数,将回波数据与共轭补偿函数相乘,消除回波信号中由电离层带来的误差。从而,恢复出无电离层干扰下的回波信号,再通过常规成像算法即得到最终的聚焦良好的目标图像。
与已有技术相比,本发明的有益效果体现在:
1)本发明针对地基ISAR对空间目标成像时,电离层对雷达回波信号会产生影响导致成像模糊的问题,提出一种能够有效去除信号中电离层带来的误差的处理方法;
2)本发明相较于传统ISAR成像质量发生图像偏移、几何畸变和分辨率降低,通过电离层影响补偿,成像结果的图像熵降低而图像对比度提高,成像效果有明显提升;
3)本发明提出对地基ISAR成像中电离层影响校正,有利于发展地基低波段ISAR雷达进行空间目标监测,同时对低波段SAR成像雷达平台的建设具有重要的理论借鉴意义。
附图说明
图1是国际空间站的光学图像。
图2是存在电离层干扰情况下国际空间站利用常规距离多普勒算法成像的结果。
图3是对电离层校正后的数据利用距离多普勒算法成像的结果。
具体实施方式
基于子带分割的低波段ISAR电离层校正方法,包括以下步骤:
步骤一:子带分割。将距离脉压和距离走动校正后的数据转换到距离频域-方位时域,将距离向信号带宽划分成距离子带;再按方位合成孔径时间将回波数据继续划分成多个方位子孔径;
步骤二:电离层干扰误差计算。在距离时域方位时域中,标出每块数据的每一次方位慢时间回波脉压峰值所处的位置,然后由该块中峰值位置的平均值乘以距离分辨单元大小得到该块数据对应的残留斜距误差;
步骤三:误差校正,信号恢复。对比拟合多项式和推导多项式可得每个方位子孔径对应的电离层参数,建立电离层干扰的校正函数,将回波数据与校正函数相乘,就可以消除主要由电离层带来的干扰误差,恢复无电离层影响下的回波信号。
至此,基于子带分割的低波段ISAR电离层校正处理基本完成。
以下通过实测的地基ISAR实验数据进一步说明本发明的有效性。
实测地基ISAR实验:
1.实验设置:
实验中,P波段地基ISAR的参数为:信号载频560MHz,带宽40MHz,PRF1000MHz,初始斜距500km。以国际空间站目标(the International Space Station,ISS)为观测目标,用回波数据中估计得到的TEC对回波数据进行补偿后成像。
2.结果分析:
图1是本次实验观测目标——国际空间站的光学图像。
图2是在包络对齐和方位相位自聚焦技术下,利用常规距离多普勒算法成像的结果。从图2的成像结果中可见,虽然相位自聚焦技术使方位向完成聚焦,但破坏了原始的回波相位,无法从最终的成像结果中辨别出国际空间站目标的姿态。
图3是进行二维子带划分估计并对电离层影响进行补偿之后,利用距离多普勒算法成像的结果。
对比图2和图3可见,补偿掉电离层影响之后,成像结果的质量大大提高了,可以清晰辨别出国际空间站的姿态和太阳能电池板。
Claims (1)
1.一种基于子带分割的低波段ISAR电离层校正方法,其特征在于:包括以下步骤:
(1)、二维子带数据分割,具体过程如下:
1a)、按信号带宽划分:
将距离脉压和距离走动校正后的数据转换到距离频域-方位时域,将距离向信号带宽划分成M个距离子带,每个子带的中心频率(fm)互不相同;
1b)、按方位合成孔径时间划分:
按方位合成孔径时间将回波数据继续划分成N个方位子孔径,每个方位子孔径对应的数据获取时间(n)互不相同;之内可以看作电离层参数TEC近似不变;
(2)、电离层干扰误差计算,具体过程如下:
电离层对回波信号产生的干扰误差体现在回波的相位对应的斜距中,得到斜距就可以反推出电离层干扰的大小。对每一个二维子带划分后的块数据,进行回波脉压处理,得到脉压峰值位置,由峰值位置的平均值与距离分辨单元大小的乘积得到该块数据对应的残留斜距误差
(3)、电离层误差校正,信号恢复,具体步骤如下:
3a)、误差向量拟合:
对于同一个方位子孔径时间内的M个距离子带中估计得到的M个残留斜距误差记成误差向量:
对误差向量在中心载频(fc)处进行多项式拟合,结果可以写成:
其中,f=fc+fr,fc是信号的中心载频,fr是距离频率且fr∈[-B/2,B/2],B是信号带宽,p0、p1和p2分别是各次多项式的系数;
3b)、电离层参数计算:
在研究某个方位慢时间对应的电离层影响时,斜距误差和电离层参数TEC之间的关系可近似推导写为:
其中,K=40.28m3/s2为电离层常数,tm为慢时间,c为真空中的光速,上式中距离频率fr对应项的系数与拟合多项式的系数应一一对应相等,利用fr的一次项系数可得每个方位子孔径对应的电离层参数:
3c)、无干扰信号恢复:
利用每个方位子孔径对应的电离层参数,得到电离层参数TEC在整个方位观测时间内的变化,建立电离层干扰的校正函数,将回波数据与校正函数相乘,消除回波信号中由电离层带来的误差。从而,恢复出无电离层干扰下的回波信号,再通过常规成像算法即得到最终的聚焦良好的目标图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610804971.0A CN106154268A (zh) | 2016-09-06 | 2016-09-06 | 一种基于子带分割的低波段isar电离层校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610804971.0A CN106154268A (zh) | 2016-09-06 | 2016-09-06 | 一种基于子带分割的低波段isar电离层校正方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106154268A true CN106154268A (zh) | 2016-11-23 |
Family
ID=57340695
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610804971.0A Withdrawn CN106154268A (zh) | 2016-09-06 | 2016-09-06 | 一种基于子带分割的低波段isar电离层校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106154268A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108008367A (zh) * | 2017-12-01 | 2018-05-08 | 中国人民解放军国防科技大学 | 星载单航过InSAR***电离层误差校正方法 |
CN108931761A (zh) * | 2018-06-05 | 2018-12-04 | 中国科学院电子学研究所 | 基于星载合成孔径雷达的电离层定标方法及*** |
CN108983230A (zh) * | 2018-05-28 | 2018-12-11 | 长安大学 | 一种基于sar方位向偏移的电离层层析构建方法 |
CN109143188A (zh) * | 2018-10-10 | 2019-01-04 | 中国科学院遥感与数字地球研究所 | Tops哨兵-1数据电离层校正方法 |
CN111679258A (zh) * | 2020-05-18 | 2020-09-18 | 中国科学院空天信息创新研究院 | 一种雷达回波数据处理方法及装置、设备 |
CN117289277A (zh) * | 2023-11-27 | 2023-12-26 | 中山大学 | 一种基于子带分割合成的多频雷达三维成像方法及*** |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000310674A (ja) * | 1999-04-09 | 2000-11-07 | Vsis Inc | 2衛星を用いた単一周波数のgps受信機に対するイオン層修正法 |
CN101487883A (zh) * | 2009-03-09 | 2009-07-22 | 北京航空航天大学 | 一种适用于多模卫星导航***的电离层格网校正方法 |
CN103760534A (zh) * | 2014-01-19 | 2014-04-30 | 中国人民解放军国防科学技术大学 | 一种星载sar数据的电离层色散效应校正方法 |
-
2016
- 2016-09-06 CN CN201610804971.0A patent/CN106154268A/zh not_active Withdrawn
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000310674A (ja) * | 1999-04-09 | 2000-11-07 | Vsis Inc | 2衛星を用いた単一周波数のgps受信機に対するイオン層修正法 |
CN101487883A (zh) * | 2009-03-09 | 2009-07-22 | 北京航空航天大学 | 一种适用于多模卫星导航***的电离层格网校正方法 |
CN103760534A (zh) * | 2014-01-19 | 2014-04-30 | 中国人民解放军国防科学技术大学 | 一种星载sar数据的电离层色散效应校正方法 |
Non-Patent Citations (2)
Title |
---|
周芳: "P波段ISAR信号电离层污染相位校正方法", 《***工程与电子技术》 * |
周芳: "高分辨SAR/ISAR成像信号补偿新技术研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108008367A (zh) * | 2017-12-01 | 2018-05-08 | 中国人民解放军国防科技大学 | 星载单航过InSAR***电离层误差校正方法 |
CN108008367B (zh) * | 2017-12-01 | 2019-09-13 | 中国人民解放军国防科技大学 | 星载单航过InSAR***电离层误差校正方法 |
CN108983230A (zh) * | 2018-05-28 | 2018-12-11 | 长安大学 | 一种基于sar方位向偏移的电离层层析构建方法 |
CN108931761A (zh) * | 2018-06-05 | 2018-12-04 | 中国科学院电子学研究所 | 基于星载合成孔径雷达的电离层定标方法及*** |
CN109143188A (zh) * | 2018-10-10 | 2019-01-04 | 中国科学院遥感与数字地球研究所 | Tops哨兵-1数据电离层校正方法 |
CN109143188B (zh) * | 2018-10-10 | 2020-07-31 | 中国科学院遥感与数字地球研究所 | Tops哨兵-1数据电离层校正方法 |
CN111679258A (zh) * | 2020-05-18 | 2020-09-18 | 中国科学院空天信息创新研究院 | 一种雷达回波数据处理方法及装置、设备 |
CN111679258B (zh) * | 2020-05-18 | 2023-09-15 | 中国科学院空天信息创新研究院 | 一种雷达回波数据处理方法及装置、设备 |
CN117289277A (zh) * | 2023-11-27 | 2023-12-26 | 中山大学 | 一种基于子带分割合成的多频雷达三维成像方法及*** |
CN117289277B (zh) * | 2023-11-27 | 2024-01-30 | 中山大学 | 一种基于子带分割合成的多频雷达三维成像方法及*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106154268A (zh) | 一种基于子带分割的低波段isar电离层校正方法 | |
EP3144702B1 (en) | Method and device for synthethic aperture radar imaging based on non-linear frequency modulation signal | |
CN105785327B (zh) | 频率分集阵列合成孔径雷达高分辨宽测绘带成像方法 | |
CN105759263B (zh) | 一种高分辨率大场景下的星载斜视sar雷达成像方法 | |
CN102176016B (zh) | 一种大斜视滑动聚束sar成像处理方法 | |
CN104977582A (zh) | 一种实现扫描雷达方位超分辨成像的解卷积方法 | |
CN103454635B (zh) | 基于高超声速飞行器平飞段的前斜视sar成像方法 | |
CN104950305A (zh) | 一种基于稀疏约束的实波束扫描雷达角超分辨成像方法 | |
CN104133215B (zh) | 基于距离徙动微调和子带分割的同步轨道雷达成像方法 | |
CN104316924A (zh) | 机载超高分辨率sar反投影图像的自聚焦运动补偿方法 | |
CN102183762A (zh) | 一种压缩感知合成孔径雷达数据获取与成像方法 | |
US20140009324A1 (en) | Autofocus-based compensation (abc) system and method for a hovering ground moving target indication (gmti) sensor | |
CN103969628B (zh) | 一种基于压缩感知的合成孔径雷达pfa成像算法 | |
CN101581780A (zh) | 一种用于侧视层析合成孔径雷达的三维聚焦成像方法 | |
CN104678386B (zh) | 一种利用gnss海面反射信号相关功率探测目标的方法 | |
CN103091674A (zh) | 基于hrrp序列的空间目标高分辨成像方法 | |
KR102142674B1 (ko) | 위성레이더 오프셋트래킹 지표변위 모델 기반 위성레이더 간섭영상의 불구속화 방법 및 그 장치 | |
CN107918124A (zh) | 带有方位空变校正的机载大斜视高分辨sar成像方法 | |
CN105372657A (zh) | 基于回波数据的视频合成孔径雷达运动补偿成像方法 | |
CN102707283A (zh) | 一种固定站双基地合成孔径雷达成像方法 | |
CN103792535A (zh) | 一种利用sar卫星测量电离层tec值的方法 | |
CN106950565A (zh) | 星载sar成像抖动补偿方法、成像方法 | |
CN105676190A (zh) | 一种校正合成孔径雷达回波数据的方法和装置 | |
CN107607945A (zh) | 一种基于空间嵌入映射的扫描雷达前视成像方法 | |
CN103576148B (zh) | 模拟星载sar距离模糊噪声图像的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20161123 |
|
WW01 | Invention patent application withdrawn after publication |