CN104898100A - 基于弱布设角反射器的机载sar辐射外定标处理方法 - Google Patents
基于弱布设角反射器的机载sar辐射外定标处理方法 Download PDFInfo
- Publication number
- CN104898100A CN104898100A CN201510219899.0A CN201510219899A CN104898100A CN 104898100 A CN104898100 A CN 104898100A CN 201510219899 A CN201510219899 A CN 201510219899A CN 104898100 A CN104898100 A CN 104898100A
- Authority
- CN
- China
- Prior art keywords
- sar
- radar
- scaler
- target
- calibration
- 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
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
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辐射外定标处理方法,包括:获得定标场原始SAR图像数据;在点目标雷达方程的基础上,重建辐射外定标方程;进行定标器图像能量的稳定估计和定标器RCS的几何光学解算;进行雷达系统传递函数的拟合与校正,得到初步辐射定标SAR图像数据;进行定标器图像能量的稳定估计和定标器RCS的几何光学解算;对系统常数项进行估算,生成定标数据,得到最终辐射定标SAR图像数据。在点目标雷达方程的基础上,依据雷达系统参数的自身特性,重建包含了变化系统传递函数和不变系统常数项的SAR辐射外定标方程,简化了辐射外定标的应用过程,避免了机载SAR飞行很难获得类似星载SAR辐射定标常用的大范围同质稳定热带雨林区域的困境。
Description
技术领域
本发明涉及机载合成孔径雷达SAR数据处理技术领域,尤其是一种基于弱布设角反射器的机载SAR辐射外定标处理方法。
背景技术
SAR辐射定标的目的就是在获得的SAR图像中,建立每个图像像元强度值与目标真实后向散射系数的正确对应关系,使得SAR图像强度与被观测目标的物理参数准确关联,有利于SAR数据的定量应用。按照定标过程中是否需要借助SAR系统外部的定标参照目标,SAR辐射定标一般可以分为内定标和外定标,辐射外定标优于内定标之处在于它是直接测量端到端的系统性能,能够测量出那些很难进行测量的系统参数,如天线方向图、波束中心增益和角度以及信号传播效应等。外定标中采用的定标目标分为点目标和分布目标两类。相对于摆放背景区域,点目标要求具有相当大的并且事先已知的雷达散射截面,常见的有角反射器,雷达转发器等,目前角反射器一般采用最大RCS值作为已知值,因此需要根据载机平台飞行路线和雷达入射角来调整每个角反射器方位角和俯仰角,然而,由于距离向上雷达入射角的变化以及机载SAR平台飞行的不稳定性,实际想达到精确调整的目的是很困难的;同时,由于SAR系统聚焦程度以及图像分辨率的不同,点目标在SAR图像上并非表现为单个强像素点,而往往表现为多个强像素点的集合,因此带来了点目标像素峰值与实际点目标图像能量之间的差异性影响定标精度。而分布目标一般应具有后向散射稳定的、时不变的、各向同性的性质,如热带雨林,大片的均匀草地等,相对于星载SAR,由于机载SAR飞行范围有限,一般难以获得热带雨林这类已知散射特性稳定的大范围分布目标,而采用其它稳定性未知的草地、植被等分布目标会降低SAR辐射外定标的精度。
发明内容
本发明的目的在于提供一种简化辐射外定标的应用过程,避免了机载SAR飞行很难获得类似星载SAR辐射定标常用的大范围同质稳定热带雨林区域的困境,解决了常规方法为了满足定标器最大RCS值带来的严格布设的难题,从而大大增强了对SAR数据进行辐射外定标处理的适用性和鲁棒性的基于弱布设角反射器的机载SAR辐射外定标处理方法。
为实现上述目的,本发明采用了以下技术方案:一种基于弱布设角反射器的机载SAR辐射外定标处理方法,该方法包括下列顺序的步骤:
(1)获得定标场原始SAR图像数据;
(2)在点目标雷达方程的基础上,重建辐射外定标方程;
(3)进行定标器图像能量的稳定估计和定标器RCS的几何光学解算;
(4)进行雷达系统传递函数的拟合与校正,得到初步辐射定标SAR图像数据;
(5)进行定标器图像能量的稳定估计和定标器RCS的几何光学解算;
(6)对系统常数项进行估算,生成定标数据,得到最终辐射定标SAR图像数据。所述定标场原始SAR图像数据是通过机载SAR雷达平台对地面定标场场景成像的结果,定标场选择平坦开阔的场景,其大小大于等于SAR成像幅宽;然后在定标场中沿距离向即垂直于载机飞行方向上均匀布设多个正方形三面角反射器,避开高塔桥梁房屋和电子力线的具有强散射能力的人工目标;接着,通过SAR雷达载机平台匀速直线飞行,并保持其飞行航线与定标角反射器的指向垂直,SAR雷达波束侧视扫过地面定标场成像,最终获得含角反射器目标的定标场原始SAR图像数据。
所述点目标雷达方程的基本公式如下:
式中:
Ps——目标的雷达回波功率;
Pt——雷达的发射功率;
Gr——接收增益;
——双程天线方向图增益;
——天线俯仰角,度;
λ——雷达波长,m;
σ0——目标的雷达后向散射系数(RCS);
c——电磁波真空速度,m/s;
τp——脉冲宽度,s;
R——目标的斜距距离,m;
La——雷达的天线长度,m;
θ——目标的雷达入射角,度;
对于SAR辐射外定标,其目的是建立SAR图像中像元能量值与绝对的雷达后向散射系数之间的一一对应关系,因此,重新定义两个参数,即雷达系统传递函数 和雷达系统常数项Ks:
其中,斜距R、入射角θ是与特定目标有关的已知变量,是雷达天线方向图,除R、θ和外的所有其它参数即为雷达系统常数项:得到辐射外定标方程,即目标在SAR图像中的能量值与目标的RCS之间的关系式:
其中,目标的雷达回波功率Ps与目标在SAR图像上的能量响应相对应,因此,辐射外定标问题简化为和Ks的估计问题,通过获取实验场的SAR图像数据和布设的已知散射特性定标器来估计SAR系统传递函数和系统常数Ks,获取成像目标真实雷达后向散射系数σ0值,由此实现SAR数据的辐射外定标。
所述进行定标器图像能量的稳定估计是指:
定标器图像能量即回波功率的估计要在待估定标器的周围选择大小不同的两个处理窗口,即目标窗和背景窗:
a)目标窗即目标区,该区域要求包含待估定标器目标所有主旁瓣像素能量,目标应置于一个大的均匀背景区域;
b)背景窗要求包含目标窗,且与目标窗之间的阴影区域为背景杂波区,背景杂波应为同质均匀背景地物;
设目标区的像素个数为Npu,背景杂波区的像素个数为Nu,由于对于数字图像数据,能量即功率用求和来计算,这样目标区的能量Wpu为
式中:
apu——目标区内的像素能量值;
同样,背景杂波区的能量Wu为
式中:
au——背景杂波区内的像素能量值;
忽略系统噪声影响,那么待估定标器目标的能量Wp通过下式计算得到:
Wp=Wpu-Wu(Npu/Nu) (6)。
所述进行定标器RCS的几何光学解算是指:
正方形三面角反射器的RCS最大值为:
其中,b为三面角反射器边长,λ是雷达波长;
当雷达波长小于角反射器尺寸时,采用几何光学模型根据观测几何关系估计RCS值,定义两个观测角度,即方位角η,0<η<π/2,和俯仰角φ,0<φ<π/2,其中方位角η定义为雷达视线方向在xoy平面的投影与ox之间的夹角,俯仰角φ定义为雷达视线方向与xoy平面之间的夹角,这两个角度由各个定标器与雷达观测方向之间的实际几何关系决定,雷达入射角θ可由下式计算:
其中,H为载机平台飞行高度,R0为SAR图像的初始斜距,Y为目标在图像中沿距离向的偏移像素值,ΔY为图像距离向分辨率,arcos是反余弦函数;
因此,弱布设定标器表现在任意几何关系下,通过几何光学的方法可计算出定标器即正方形三面角反射器的RCS值σsquare(φ,η):
其中,c1,c2,c3为方位角和俯仰角的函数,分别定义如下:
由知可知,当方位角η=π/4,即45度;俯仰角即35.26度时,达到RCS最大值。
所述进行雷达系统传递函数的拟合与校正是指:
假设在距离向上布设了N个定标器,即正方形三面角反射器,那么第i个,
1≤i≤N,定标器的系统传递函数为:
σsquare,i(φ,η)表示第i个定标器的RCS值,Wpi表示第i个定标器的目标能量估值,η为方位角,φ为俯仰角;
由于是随距离向参数变化的函数,利用多个定标器估值的曲线拟合来获得,并通过下式校正原始SAR数据的系统传递函数,得到初步辐射定标SAR图像数据,其图像能量为Ps':
在得到初步辐射定标SAR图像数据基础上,再进行定标器图像能量的稳定估计和定标器RCS的几何光学解算,同样获得多个定标器图像能量的估计Wp'i以及RCS的计算值,由于初步辐射定标SAR图像数据在经过雷达系统传递函数的拟合与校正之后已消除了雷达系统传递函数的影响,因此,此时得到的是各个定标器的系 统常数Ksi:
系统常数Ks采用多个Ksi值估计的平均,最后,校正初步辐射校正SAR数据中的系统常数Ks得到最终SAR辐射定标数据,即地物目标RCS值σ0:
σ0=Ps'/Ks (14)。
由上述技术方案可知,本发明的优点如下:第一,在已知点目标雷达方程的基础上,针对辐射外定标要获得雷达后向散射系数的需求,依据雷达系统参数的自身特性,重建包含了变化系统传递函数和不变系统常数项的SAR辐射外定标方程,简化了辐射外定标的应用过程,避免了机载SAR飞行很难获得类似星载SAR辐射定标常用的大范围同质稳定热带雨林区域的困境;第二,通过SAR图像中定标器目标窗以及背景窗之间的能量计算实现定标器图像能量的稳定、准确的估计,并且最终结果不受所选窗口大小、图像分辨率、系统聚集程度的影响;第三,采用正方形三面角反射器作为已知散射特性的定标器,根据布设定标器与雷达视线方向所呈的方位角和俯仰角,利用定标器的几何光学解而非理论最大值作为其实际RCS值,不但避免了使用常规理论RCS最大值引起的实际估计误差,也弱化了机载SAR情况下对定标器的布设要求;第四,在辐射外定标方程重建的基础上,利用连续两次定标器图像能量估计和RCS几何光学解算步骤,通过距离向多个定标器样本采样,一方面采用曲线拟合方法估算系统传递函数的变化曲线,一方面采用样本均值估计辐射外定标的系统常数项,从而最终达到辐射外定标获取图像目标的真实RCS值的目的。
附图说明
图1为本发明的方法流程图;
图2为本发明定标器图像能量估计示意图;
图3为本发明正方形三面角反射器的雷达几何关系图;
图4为本发明实施例一原始SAR图像;
图5为本发明实施例一系统传递函数的曲线拟合估计;
图6为本发明实施例一系统传递函数校正后的SAR图像;
图7为本发明实施例一用于辐射精度验证的定标器图像;
图8为本发明实施例一的典型地物RCS值的图表;
图9为本发明实施例二原始SAR图像;
图10为本发明实施例二系统传递函数的曲线拟合估计;
图11为本发明实施例二系统传递函数校正后的SAR图像;
图12为本发明实施例二用于辐射精度验证的定标器图像;
图13为本发明实施例二的典型地物RCS值的图表。
具体实施方式
一种基于弱布设角反射器的机载SAR辐射外定标处理方法,包括:获得定标场原始SAR图像数据;在点目标雷达方程的基础上,重建辐射外定标方程;进行定标器图像能量的稳定估计和定标器RCS的几何光学解算;进行雷达系统传递函数的拟合与校正,得到初步辐射定标SAR图像数据;进行定标器图像能量的稳定估计和定标器RCS的几何光学解算;对系统常数项进行估算,生成定标数据,得到最终辐射定标SAR图像数据,如图1所示。
所述定标场原始SAR图像数据是通过机载SAR雷达平台对地面定标场场景成像的结果,定标场选择平坦开阔的场景,其大小大于等于SAR成像幅宽;然后在定标场中沿距离向即垂直于载机飞行方向上均匀布设多个正方形三面角反射器,避开高塔桥梁房屋和电子力线的具有强散射能力的人工目标;接着,通过SAR雷达载机平台匀速直线飞行,并保持其飞行航线与定标角反射器的指向垂直,SAR雷达波束侧视扫过地面定标场成像,最终获得含角反射器目标的定标场原始SAR图像数据。
所述点目标雷达方程的基本公式如下:
式中:
Ps——目标的雷达回波功率;
Pt——雷达的发射功率;
Gr——接收增益;
——双程天线方向图增益;
——天线俯仰角,度;
λ——雷达波长,m;
σ0——目标的雷达后向散射系数(RCS);
c——电磁波真空速度,m/s;
τp——脉冲宽度,s;
R——目标的斜距距离,m;
La——雷达的天线长度,m;
θ——目标的雷达入射角,度;
对于SAR辐射外定标,其目的是建立SAR图像中像元能量值与绝对的雷达后向散射系数之间的一一对应关系,因此,重新定义两个参数,即雷达系统传递函数 和雷达系统常数项Ks:
其中,斜距R、入射角θ是与特定目标有关的已知变量,是雷达天线方向图,除R、θ和外的所有其它参数即为雷达系统常数项:得到辐射外定标方程,即目标在SAR图像中的能量值与目标的RCS之间的关系式:
其中,目标的雷达回波功率Ps与目标在SAR图像上的能量响应相对应,因此,辐射外定标问题简化为和Ks的估计问题,通过获取实验场的SAR图像数据和布设的已知散射特性定标器来估计SAR系统传递函数和系统常数Ks,获取成像目标真实雷达后向散射系数σ0值,由此实现SAR数据的辐射外定标。
所述进行定标器图像能量的稳定估计是指:
定标器图像能量即回波功率的估计要在待估定标器的周围选择大小不同的两个处理窗口,即目标窗和背景窗:
a)目标窗即目标区,该区域要求包含待估定标器目标所有主旁瓣像素能量,目标应置于一个大的均匀背景区域;
b)背景窗要求包含目标窗,且与目标窗之间的阴影区域为背景杂波区,背景杂波应为同质均匀背景地物;
设目标区的像素个数为Npu,背景杂波区的像素个数为Nu,由于对于数字图像数据,能量即功率用求和来计算,这样目标区的能量Wpu为
式中:
apu——目标区内的像素能量值;
同样,背景杂波区的能量Wu为
式中:
au——背景杂波区内的像素能量值;
忽略系统噪声影响,那么待估定标器目标的能量Wp通过下式计算得到:
Wp=Wpu-Wu(Npu/Nu) (6)。
所述进行定标器RCS的几何光学解算是指:
正方形三面角反射器的RCS最大值为:
其中,b为三面角反射器边长,λ是雷达波长;
当雷达波长小于角反射器尺寸时,采用几何光学模型根据观测几何关系估计RCS值,定义两个观测角度,即方位角η,0<η<π/2,和俯仰角φ,0<φ<π/2,其中方位角η定义为雷达视线方向在xoy平面的投影与ox之间的夹角,俯仰角φ定义为雷达视线方向与xoy平面之间的夹角,这两个角度由各个定标器与雷达观测方向之间的实际几何关系决定,雷达入射角θ可由下式计算:
其中,H为载机平台飞行高度,R0为SAR图像的初始斜距,Y为目标在图像中沿 距离向的偏移像素值,ΔY为图像距离向分辨率,arcos是反余弦函数;
因此,弱布设定标器表现在任意几何关系下,通过几何光学的方法可计算出定标器即正方形三面角反射器的RCS值σsquare(φ,η):
其中,c1,c2,c3为方位角和俯仰角的函数,分别定义如下:
由知可知,当方位角η=π 4,即45度;俯仰角即35.26度时,达到RCS最大值。
所述进行雷达系统传递函数的拟合与校正是指:
假设在距离向上布设了N个定标器,即正方形三面角反射器,那么第i个,
1≤i≤N,定标器的系统传递函数为:
σsquare,i(φ,η)表示第i个定标器的RCS值,Wpi表示第i个定标器的目标能量估值,η为方位角,φ为俯仰角;
由于是随距离向参数变化的函数,利用多个定标器估值的曲线拟合来获得,并通过下式校正原始SAR数据的系统传递函数,得到初步辐射定标SAR图像数据,其图像能量为Ps':
在得到初步辐射定标SAR图像数据基础上,再进行定标器图像能量的稳定估计和定标器RCS的几何光学解算,同样获得多个定标器图像能量的估计Wp'i以及RCS的计算值,由于初步辐射定标SAR图像数据在经过雷达系统传递函数的拟合与校 正之后已消除了雷达系统传递函数的影响,因此,此时得到的是各个定标器的系统常数Ksi:
系统常数Ks采用多个Ksi值估计的平均,最后,校正初步辐射校正SAR数据中的系统常数Ks得到最终SAR辐射定标数据,即地物目标RCS值σ0:
σ0=Ps'/Ks (14)。
由于本发明采用的是无源角反射器作为SAR辐射定标器,因此需要事先将一定数量经过精确标定的标准正方形三面角反射器沿距离向,即垂直飞行方向均匀排列于定标场中,角反射器的一般性要求为:
各个正方形三面角反射器采用相同弱条件布设方式,无需针对每个角反射器调整,具体布设方式为:参照图3将反射器的底面正方形的对角线指向与已知的SAR飞行航线方向垂直,且底面正方形平面与地面保持平行,即45度方位角,俯仰角为各个角反射器位置的实际雷达入射角。然后,通过SAR载机平台飞行将雷达照射波束扫过定标场和成像处理得到的SAR图像中所有标准反射器目标,最后实现SAR图像辐射外定标处理。
下面结合图1至图13对本发明作进一步说明。
实施例一
C波段陆上有人机SAR数据
采用国内中国电科38所有人机SAR的陆上试验飞行在某地定标场获得的C波段机载SAR图像,图4为该定标场的原始HH极化通道的SAR图像,图像大小为3480像素*2270像素,载机飞行高度为2050m,图像初始斜距为2950m,图像距离向分辨率为0.75m,图像方位向分辨率为1.0m。由于系统传递函数的存 在,原始SAR图像表现出明显的明暗变化,即图像上部的亮度偏亮,图像底部的亮度则偏暗。
首先建立起SAR图像能量与系统传递函数、系统常数以及目标RCS之间的辐射外定标方程,然后在此方程基础上,参照图2在图4定标场中沿机场跑道方向即距离向上布设的大量正方形三面角反射器的周围选取合适窗口进行各个定标器的
图像能量估计,表1给出各个角反射器的图像目标能量Wpi的估值。
同时根据实际SAR成像几何关系计算相应各个定标器的σtriangle,i(φ,η)值;已知角反射器的边长b为0.35m,方位角η为45度,俯仰角φ=90-θ,由式(8)和式(9)可得到各个位置角反射器的RCS值,见表1:
表1实施例一的系统传递函数计算表
表2实施例一的系统常数计算表
然后在距离向上通过多项式曲线拟合出系统传递函数当获得系统传递函数的拟合曲线后,可将原始SAR图像进行系统传递函数的校正。
图5是利用距离向上布设的定标器估计拟合得到的系统传递函数曲线。图6是原始SAR图像经系统传递函数校正后的SAR图像,可以看出图4中底部的阴暗图像在图6中呈现出了原有的地物细节,增强了后期地物解译的能力。
接着,在图6中的初步校正SAR图像上进行定标器图像能量的稳定估计和定标器RCS的几何光学解算,并估计系统常数,见表2,最终的系统常数Ks采用表2中5个Ksi估值的平均,即12.4174dB。然后,对系统传递函数校正后的SAR图像进一步进行系统定标常数Ks的校正,得到完成辐射外定标后的定标场SAR图像。图8给出了完成辐射外定标后定标场中一些典型地物(如跑道、草地、植被、城区、农田)的RCS值与《各类地物雷达散射统计手册》进行的对比,吻合各种典型地物的实测RCS值。《各类地物雷达散射统计手册》由美国F.T.Ulaby和M.C.Dobson编著,全面描述了各种传感器参数下各类地物的“地杂波”测量数据,在已有大量实测数据的基础上,对数据进行了严格反复的鉴别、剔除和校准并界定精确度后编入本手册。因此,迄今为止这本手册被认为对已校准过的后向散射系数σ0数据是最全面的。
最后,利用图7中停机坪处摆放的0.4三角形三面体角反射器和0.35m正方形三面体角反射器进行了辐射精度验证,其辐射定标精度的均值为-0.5349dB,均方差为1.1597dB,满足1.5dB的辐射外定标要求。
实施例二
C波段海上无人机SAR数据
采用国内中国电科38所无人机SAR的海上实验飞行在某沿海定标场获得的C波段机载SAR图像。海上飞行实验的思路是无人机先飞过临海的地面定标场,然后进入内海进行成像飞行,对地面定标场进行外辐射定标处理获得的定标参数,再用于内海SAR数据的外辐射定标,最后得到海洋SAR辐射外定标数据可用于海洋风速、风向以及波向、波长、波高等信息的定量反演。图9为该定标场的原始HH极化通道的SAR图像,图像大小为2900像素*1470像素,载机飞行高度为2941m,图像初始斜距为3627m,图像距离向分辨率为0.7495m,图像方位向分辨率为0.6m。由于系统传递函数的存在,该原始SAR图像同样表现出明显的明暗变化,即图像上部下部的亮度偏暗亮,图像中部的亮度偏亮。
与实施例一一样,首先参照图2在定标场图像中沿距离向上布设的大量正方形三面角反射器的周围选取合适窗口进行各个定标器的图像能量估计,同时根据实际SAR成像几何关系计算相应各个定标器的RCS值,计算结果见表3。
表3实施例二的系统传递函数计算表
表4实施例二的系统常数计算表
然后在距离向上通过多项式曲线拟合出系统传递函数当获得系统传递函数的拟合曲线后,再将原始SAR图像进行系统传递函数的校正。图10是利用距离向上布设的定标器估计拟合得到的系统传递函数曲线。图11是原始SAR图像经系统传递函数校正后的SAR图像,可以看出图9中上部下部的阴暗图像在图11中呈现出了原有的地物细节。
接着,在图11中的初步校正SAR图像上进行定标器图像能量的稳定估计和定标器RCS的几何光学解算,并通估计系统常数,见表4,最终的系统常数Ks采用表4中5个Ksi估值的平均,即6.9090dB。然后,对系统传递函数校正后的SAR图像进一步进行系统定标常数Ks的校正,得到完成辐射外定标后的定标场SAR图像。
在完成辐射外定标处理后,在成像定标场景选取了跑道和城区典型地物进行RCS分析,其结果与《各类地物雷达散射统计手册》实测结果基本吻合,如图13所示。
最后,利用图12中其它位置摆放的多个正方形三面体角反射器进行了辐射精度验证,得到的各定标精度分别为-0.9407dB,-2.2676dB,0.1546dB,-1.2984dB,0.0899dB,-0.4644dB,-0.6292dB,-0.6219dB,-0.5536dB,其定标精度标准差为0.9dB。在保证海洋成像中系统稳定性的条件下,满足利用经典物理参数模型反演海洋风速、风向、波浪信息的1.0dB辐射定标精度。
综上所述,本发明在已知点目标雷达方程的基础上,针对辐射外定标要获得雷达后向散射系数的需求,依据雷达系统参数的自身特性,重建包含了变化系统传递函数和不变系统常数项的SAR辐射外定标方程,简化了辐射外定标的应用过程,避免了机载SAR飞行很难获得类似星载SAR辐射定标常用的大范围同质稳定热带雨林区域的困境,同时利用角反射器的几何光学模型解,解决了常规方法为了满足定标器最大RCS值带来的严格布设的难题,从而大大增强了对SAR数据进行辐射外定标处理的适用性和鲁棒性。
Claims (7)
1.一种基于弱布设角反射器的机载SAR辐射外定标处理方法,该方法包括下列顺序的步骤:
(1)获得定标场原始SAR图像数据;
(2)在点目标雷达方程的基础上,重建辐射外定标方程;
(3)进行定标器图像能量的稳定估计和定标器RCS的几何光学解算;
(4)进行雷达系统传递函数的拟合与校正,得到初步辐射定标SAR图像数据;
(5)进行定标器图像能量的稳定估计和定标器RCS的几何光学解算;
(6)对系统常数项进行估算,生成定标数据,得到最终辐射定标SAR图像数据。
2.根据权利要求1所述的基于弱布设角反射器的机载SAR辐射外定标处理方法,其特征在于:所述定标场原始SAR图像数据是通过机载SAR雷达平台对地面定标场场景成像的结果,定标场选择平坦开阔的场景,其大小大于等于SAR成像幅宽;然后在定标场中沿距离向即垂直于载机飞行方向上均匀布设多个正方形三面角反射器,避开高塔桥梁房屋和电子力线的具有强散射能力的人工目标;接着,通过SAR雷达载机平台匀速直线飞行,并保持其飞行航线与定标角反射器的指向垂直,SAR雷达波束侧视扫过地面定标场成像,最终获得含角反射器目标的定标场原始SAR图像数据。
3.根据权利要求1所述的基于弱布设角反射器的机载SAR辐射外定标处理方法,其特征在于:所述点目标雷达方程的基本公式如下:
式中:
Ps——目标的雷达回波功率;
Pt——雷达的发射功率;
Gr——接收增益;
——双程天线方向图增益;
——天线俯仰角,度;
λ——雷达波长,m;
σ0——目标的雷达后向散射系数(RCS);
c——电磁波真空速度,m/s;
τp——脉冲宽度,s;
R——目标的斜距距离,m;
La——雷达的天线长度,m;
θ——目标的雷达入射角,度;
对于SAR辐射外定标,其目的是建立SAR图像中像元能量值与绝对的雷达后向散射系数之间的一一对应关系,因此,重新定义两个参数,即雷达系统传递函数和雷达系统常数项Ks:
其中,斜距R、入射角θ是与特定目标有关的已知变量,是雷达天线方向图,除R、θ和外的所有其它参数即为雷达系统常数项:得到辐射外定标方程,即目标在SAR图像中的能量值与目标的RCS之间的关系式:
其中,目标的雷达回波功率Ps与目标在SAR图像上的能量响应相对应,因此,辐射外定标问题简化为和Ks的估计问题,通过获取实验场的SAR图像数据和布设的已知散射特性定标器来估计SAR系统传递函数和系统常数Ks,获取成像目标真实雷达后向散射系数σ0值,由此实现SAR数据的辐射外定标。
4.根据权利要求1所述的基于弱布设角反射器的机载SAR辐射外定标处理方法,其特征在于:所述进行定标器图像能量的稳定估计是指:
定标器图像能量即回波功率的估计要在待估定标器的周围选择大小不同的两个处理窗口,即目标窗和背景窗:
a)目标窗即目标区,该区域要求包含待估定标器目标所有主旁瓣像素能量,目标应置于一个大的均匀背景区域;
b)背景窗要求包含目标窗,且与目标窗之间的阴影区域为背景杂波区,背景杂波应为同质均匀背景地物;
设目标区的像素个数为Npu,背景杂波区的像素个数为Nu,由于对于数字图像数据,能量即功率用求和来计算,这样目标区的能量Wpu为
式中:
apu——目标区内的像素能量值;
同样,背景杂波区的能量Wu为
式中:
au——背景杂波区内的像素能量值;
忽略系统噪声影响,那么待估定标器目标的能量Wp通过下式计算得到:
Wp=Wpu-Wu(Npu/Nu) (6)。
5.根据权利要求1所述的基于弱布设角反射器的机载SAR辐射外定标处理方法,其特征在于:所述进行定标器RCS的几何光学解算是指:
正方形三面角反射器的RCS最大值为:
其中,b为三面角反射器边长,λ是雷达波长;
当雷达波长小于角反射器尺寸时,采用几何光学模型根据观测几何关系估计RCS值,定义两个观测角度,即方位角η,0<η<π/2,和俯仰角φ,0<φ<π/2,其中方位角η定义为雷达视线方向在xoy平面的投影与ox之间的夹角,俯仰角φ定义为雷达视线方向与xoy平面之间的夹角,这两个角度由各个定标器与雷达观测方向之间的实际几何关系决定,雷达入射角θ可由下式计算:
其中,H为载机平台飞行高度,R0为SAR图像的初始斜距,Y为目标在图像中沿距离向的偏移像素值,ΔY为图像距离向分辨率,arcos是反余弦函数;
因此,弱布设定标器表现在任意几何关系下,通过几何光学的方法可计算出定标器即正方形三面角反射器的RCS值σsquare(φ,η):
其中,c1,c2,c3为方位角和俯仰角的函数,分别定义如下:
由知可知,当方位角η=π/4,即45度;俯仰角即35.26度时,达到RCS最大值。
6.根据权利要求1所述的基于弱布设角反射器的机载SAR辐射外定标处理方法,其特征在于:所述进行雷达系统传递函数的拟合与校正是指:
假设在距离向上布设了N个定标器,即正方形三面角反射器,那么第i个,
1≤i≤N,定标器的系统传递函数为:
σsquare,i(φ,η)表示第i个定标器的RCS值,Wpi表示第i个定标器的目标能量估值,η为方位角,φ为俯仰角;
由于是随距离向参数变化的函数,利用多个定标器估值的曲线拟合来获得,并通过下式校正原始SAR数据的系统传递函数,得到初步辐射定标SAR图像数据,其图像能量为P′s:
7.根据权利要求1所述的基于弱布设角反射器的机载SAR辐射外定标处理方法,其特征在于:在得到初步辐射定标SAR图像数据基础上,再进行定标器图像能量的稳定估计和定标器RCS的几何光学解算,同样获得多个定标器图像能量的估计W′pi以及RCS的计算值,由于初步辐射定标SAR图像数据在经过雷达系统传递函数的拟合与校正之后已消除了雷达系统传递函数的影响,因此,此时得到的是各个定标器的系统常数Ksi:
系统常数Ks采用多个Ksi值估计的平均,最后,校正初步辐射校正SAR数据中的系统常数Ks得到最终SAR辐射定标数据,即地物目标RCS值σ0:
σ0=P′s/Ks (14)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510219899.0A CN104898100B (zh) | 2015-04-30 | 2015-04-30 | 基于弱布设角反射器的机载sar辐射外定标处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510219899.0A CN104898100B (zh) | 2015-04-30 | 2015-04-30 | 基于弱布设角反射器的机载sar辐射外定标处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104898100A true CN104898100A (zh) | 2015-09-09 |
CN104898100B CN104898100B (zh) | 2017-10-03 |
Family
ID=54030864
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510219899.0A Active CN104898100B (zh) | 2015-04-30 | 2015-04-30 | 基于弱布设角反射器的机载sar辐射外定标处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104898100B (zh) |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105301570A (zh) * | 2015-10-20 | 2016-02-03 | 中国科学院电子学研究所 | 一种机载顺轨干涉sar系统的外场定标方法 |
CN105747694A (zh) * | 2016-02-17 | 2016-07-13 | 韦醒妃 | 一种具有识别功能的文物展示柜 |
CN106950551A (zh) * | 2017-03-31 | 2017-07-14 | 中国科学院电子学研究所 | 基于辐射稳定特性挖掘的sar辐射定标方法 |
CN108983172A (zh) * | 2018-08-17 | 2018-12-11 | 中国电子科技集团公司第三十八研究所 | 一种随机进化约束的sar辐射定标方法 |
CN109541560A (zh) * | 2018-11-14 | 2019-03-29 | 辽宁工程技术大学 | 一种基于通讯信号塔的辐射定标方法 |
CN109752696A (zh) * | 2017-11-06 | 2019-05-14 | 中国人民解放军61646部队 | 一种高分辨率合成孔径雷达卫星图像中角反射器rcs校正方法 |
CN110058204A (zh) * | 2019-04-25 | 2019-07-26 | 西安空间无线电技术研究所 | 一种基于方向图匹配的星载天线波束中心定标方法 |
CN110161855A (zh) * | 2019-05-21 | 2019-08-23 | 中国电子科技集团公司第三十八研究所 | 一种基于鲁棒伺服增益调度无人机控制器的设计方法 |
CN110865346A (zh) * | 2018-08-28 | 2020-03-06 | 中国人民解放军61646部队 | 一种基于直接定位算法的星载sar时间参数标定方法 |
CN111596272A (zh) * | 2020-06-01 | 2020-08-28 | 中国科学院空天信息创新研究院 | 基于角反射器的机载多极化sar系统外定标方法 |
CN111812646A (zh) * | 2020-07-01 | 2020-10-23 | 自然资源部第二海洋研究所 | 一种利用合成孔径雷达图像反演海面风速的方法及系统 |
CN112105937A (zh) * | 2019-09-12 | 2020-12-18 | 深圳市大疆创新科技有限公司 | 一种天线方向图测试方法、设备及存储介质 |
CN112346027A (zh) * | 2020-10-29 | 2021-02-09 | 中国科学院空天信息创新研究院 | 用于确定合成孔径雷达图像的散射特性的方法和系统 |
CN112444783A (zh) * | 2019-09-05 | 2021-03-05 | 中国科学院光电研究院 | 基于自然均匀目标的合成孔径雷达在轨快速定标方法 |
CN112859020A (zh) * | 2021-01-18 | 2021-05-28 | 中国科学院空天信息创新研究院 | 一种太阳能电池板对sar定标器rcs影响分析方法及装置 |
WO2021103715A1 (zh) * | 2019-11-28 | 2021-06-03 | 中国科学院合肥物质科学研究院 | 一种机载高光谱成像激光雷达系统的辐射标定方法 |
CN113406638A (zh) * | 2021-06-28 | 2021-09-17 | 南京工程学院 | 一种基于多旋翼无人机sar成像系统的参数一体化设定方法 |
CN113640758A (zh) * | 2021-08-23 | 2021-11-12 | 中国科学院空天信息创新研究院 | 一种城市复杂环境下的sar图像定标器放置方法与系统 |
CN113702926A (zh) * | 2021-07-16 | 2021-11-26 | 中国科学院空天信息创新研究院 | 一种有源定标器雷达截面积标校方法及装置 |
CN115079114A (zh) * | 2022-08-22 | 2022-09-20 | 中国科学院空天信息创新研究院 | 基于海洋与角反射器后向散射的机载sar辐射定标方法 |
CN112444783B (zh) * | 2019-09-05 | 2024-05-24 | 中国科学院光电研究院 | 基于自然均匀目标的合成孔径雷达在轨快速定标方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3742881B2 (ja) * | 2003-10-24 | 2006-02-08 | 独立行政法人情報通信研究機構 | 偏波合成開口レーダ較正方法及び装置 |
CN101082670A (zh) * | 2006-06-02 | 2007-12-05 | 中国科学院电子学研究所 | 一种宽带合成孔径雷达的有源外定标器及其定标方法 |
EP2065731A1 (en) * | 2007-11-27 | 2009-06-03 | Politecnico di Milano | Method for radiometric calibration of SAR sensors |
CN102393513A (zh) * | 2011-08-31 | 2012-03-28 | 中国电子科技集团公司第三十八研究所 | 一种基于自然分布场景和稀少定标器的极化定标技术 |
CN103530469A (zh) * | 2013-10-23 | 2014-01-22 | 武汉大学 | 一种角反射器雷达截面积的计算方法 |
-
2015
- 2015-04-30 CN CN201510219899.0A patent/CN104898100B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3742881B2 (ja) * | 2003-10-24 | 2006-02-08 | 独立行政法人情報通信研究機構 | 偏波合成開口レーダ較正方法及び装置 |
CN101082670A (zh) * | 2006-06-02 | 2007-12-05 | 中国科学院电子学研究所 | 一种宽带合成孔径雷达的有源外定标器及其定标方法 |
EP2065731A1 (en) * | 2007-11-27 | 2009-06-03 | Politecnico di Milano | Method for radiometric calibration of SAR sensors |
CN102393513A (zh) * | 2011-08-31 | 2012-03-28 | 中国电子科技集团公司第三十八研究所 | 一种基于自然分布场景和稀少定标器的极化定标技术 |
CN103530469A (zh) * | 2013-10-23 | 2014-01-22 | 武汉大学 | 一种角反射器雷达截面积的计算方法 |
Non-Patent Citations (3)
Title |
---|
A.LAURENCE GRAY ET AL: "Synthetic Aperture Radar Calibration Using Reference Reflectors", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 * |
潘德炉 等: "《海洋微波遥感与应用》", 30 September 2013 * |
袁礼海 等: ""利用点目标进行SAR辐射定标的方法研究"", 《无线电工程》 * |
Cited By (34)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105301570A (zh) * | 2015-10-20 | 2016-02-03 | 中国科学院电子学研究所 | 一种机载顺轨干涉sar系统的外场定标方法 |
CN105747694B (zh) * | 2016-02-17 | 2018-12-11 | 盐城吉大智能终端产业研究院有限公司 | 一种具有识别功能的文物展示柜 |
CN105747694A (zh) * | 2016-02-17 | 2016-07-13 | 韦醒妃 | 一种具有识别功能的文物展示柜 |
CN106950551B (zh) * | 2017-03-31 | 2019-12-31 | 中国科学院电子学研究所 | 基于辐射稳定特性挖掘的sar辐射定标方法 |
CN106950551A (zh) * | 2017-03-31 | 2017-07-14 | 中国科学院电子学研究所 | 基于辐射稳定特性挖掘的sar辐射定标方法 |
CN109752696A (zh) * | 2017-11-06 | 2019-05-14 | 中国人民解放军61646部队 | 一种高分辨率合成孔径雷达卫星图像中角反射器rcs校正方法 |
CN109752696B (zh) * | 2017-11-06 | 2020-03-10 | 中国人民解放军61646部队 | 一种高分辨率合成孔径雷达卫星图像中角反射器rcs校正方法 |
CN108983172A (zh) * | 2018-08-17 | 2018-12-11 | 中国电子科技集团公司第三十八研究所 | 一种随机进化约束的sar辐射定标方法 |
CN110865346A (zh) * | 2018-08-28 | 2020-03-06 | 中国人民解放军61646部队 | 一种基于直接定位算法的星载sar时间参数标定方法 |
CN110865346B (zh) * | 2018-08-28 | 2021-07-20 | 中国人民解放军61646部队 | 一种基于直接定位算法的星载sar时间参数标定方法 |
CN109541560A (zh) * | 2018-11-14 | 2019-03-29 | 辽宁工程技术大学 | 一种基于通讯信号塔的辐射定标方法 |
CN109541560B (zh) * | 2018-11-14 | 2023-05-23 | 辽宁工程技术大学 | 一种基于通讯信号塔的辐射定标方法 |
CN110058204A (zh) * | 2019-04-25 | 2019-07-26 | 西安空间无线电技术研究所 | 一种基于方向图匹配的星载天线波束中心定标方法 |
CN110161855A (zh) * | 2019-05-21 | 2019-08-23 | 中国电子科技集团公司第三十八研究所 | 一种基于鲁棒伺服增益调度无人机控制器的设计方法 |
CN112444783A (zh) * | 2019-09-05 | 2021-03-05 | 中国科学院光电研究院 | 基于自然均匀目标的合成孔径雷达在轨快速定标方法 |
CN112444783B (zh) * | 2019-09-05 | 2024-05-24 | 中国科学院光电研究院 | 基于自然均匀目标的合成孔径雷达在轨快速定标方法 |
CN112105937A (zh) * | 2019-09-12 | 2020-12-18 | 深圳市大疆创新科技有限公司 | 一种天线方向图测试方法、设备及存储介质 |
WO2021046797A1 (zh) * | 2019-09-12 | 2021-03-18 | 深圳市大疆创新科技有限公司 | 一种天线方向图测试方法、设备及存储介质 |
US11867846B2 (en) * | 2019-11-28 | 2024-01-09 | Hefei Institute of Physical Science, Chinese Academy of Sciences | Method for radiation calibration of airborne hyperspectral imaging LiDAR system |
WO2021103715A1 (zh) * | 2019-11-28 | 2021-06-03 | 中国科学院合肥物质科学研究院 | 一种机载高光谱成像激光雷达系统的辐射标定方法 |
US20220146651A1 (en) * | 2019-11-28 | 2022-05-12 | Hefei Institute of Physical Science, Chinese Academy of Sciences | Method for Radiation Calibration of Airborne Hyperspectral Imaging LiDAR System |
CN111596272B (zh) * | 2020-06-01 | 2022-11-29 | 中国科学院空天信息创新研究院 | 基于角反射器的机载多极化sar系统外定标方法 |
CN111596272A (zh) * | 2020-06-01 | 2020-08-28 | 中国科学院空天信息创新研究院 | 基于角反射器的机载多极化sar系统外定标方法 |
CN111812646A (zh) * | 2020-07-01 | 2020-10-23 | 自然资源部第二海洋研究所 | 一种利用合成孔径雷达图像反演海面风速的方法及系统 |
CN112346027A (zh) * | 2020-10-29 | 2021-02-09 | 中国科学院空天信息创新研究院 | 用于确定合成孔径雷达图像的散射特性的方法和系统 |
CN112346027B (zh) * | 2020-10-29 | 2022-02-08 | 中国科学院空天信息创新研究院 | 用于确定合成孔径雷达图像的散射特性的方法和系统 |
CN112859020A (zh) * | 2021-01-18 | 2021-05-28 | 中国科学院空天信息创新研究院 | 一种太阳能电池板对sar定标器rcs影响分析方法及装置 |
CN112859020B (zh) * | 2021-01-18 | 2023-05-02 | 中国科学院空天信息创新研究院 | 一种太阳能电池板对sar定标器rcs影响分析方法及装置 |
CN113406638A (zh) * | 2021-06-28 | 2021-09-17 | 南京工程学院 | 一种基于多旋翼无人机sar成像系统的参数一体化设定方法 |
CN113702926A (zh) * | 2021-07-16 | 2021-11-26 | 中国科学院空天信息创新研究院 | 一种有源定标器雷达截面积标校方法及装置 |
CN113702926B (zh) * | 2021-07-16 | 2023-09-08 | 中国科学院空天信息创新研究院 | 一种有源定标器雷达截面积标校方法及装置 |
CN113640758A (zh) * | 2021-08-23 | 2021-11-12 | 中国科学院空天信息创新研究院 | 一种城市复杂环境下的sar图像定标器放置方法与系统 |
CN115079114B (zh) * | 2022-08-22 | 2022-11-04 | 中国科学院空天信息创新研究院 | 基于海洋与角反射器后向散射的机载sar辐射定标方法 |
CN115079114A (zh) * | 2022-08-22 | 2022-09-20 | 中国科学院空天信息创新研究院 | 基于海洋与角反射器后向散射的机载sar辐射定标方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104898100B (zh) | 2017-10-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104898100A (zh) | 基于弱布设角反射器的机载sar辐射外定标处理方法 | |
CN103091666B (zh) | 非理想定标器条件下的机载p波段极化sar定标方法 | |
CN109752696B (zh) | 一种高分辨率合成孔径雷达卫星图像中角反射器rcs校正方法 | |
CN108318871B (zh) | 一种星载微波辐射计冷空定标修正方法 | |
CN106405547A (zh) | 一种多极化sar简缩极化模式定标方法 | |
CN108562882B (zh) | 一种星载sar影像几何交叉定标方法和系统 | |
Mazouni et al. | 76.5 GHz millimeter-wave radar for foreign objects debris detection on airport runways | |
CN104502911A (zh) | 一种穿墙成像雷达的墙壁参数估计方法 | |
CN104865274A (zh) | 基于单极化雷达数据的土壤含水量反演方法及系统 | |
CN104101297A (zh) | 一种基于光电观测的空间物体尺度获取方法 | |
CN103645466B (zh) | 一种基于平台姿态时变性补偿的极化定标方法 | |
CN112798013A (zh) | 一种对光学载荷在轨绝对辐射定标结果进行验证的方法 | |
CN108375770B (zh) | 一种基于旋转对称性的全极化sar数据的极化定标方法 | |
CN110764153B (zh) | 一种星载微波成像仪热镜背瓣在轨误差校正系统及方法 | |
CN108983172B (zh) | 一种随机进化约束的sar辐射定标方法 | |
CN108983235B (zh) | 一种基于刈幅的交叉极化sar风场反演方法 | |
Li et al. | Land contamination analysis of SMOS brightness temperature error near coastal areas | |
CN115079114B (zh) | 基于海洋与角反射器后向散射的机载sar辐射定标方法 | |
Liu et al. | Retrieval of sea surface wind vectors from simulated satellite microwave polarimetric measurements | |
CN115113156B (zh) | 用于双极化相控阵气象雷达的标校方法和系统 | |
Marty et al. | From passive to active radar reflectors and beyond | |
Zhang et al. | A Spaceborne SAR Calibration Simulator Based on Gaofen-3 Data | |
Mazouni et al. | 77 GHz offset reflectarray for FOD detection on airport runways | |
CN112859020B (zh) | 一种太阳能电池板对sar定标器rcs影响分析方法及装置 | |
Xiao et al. | A semi-empirical topographic correction model for multi-source satellite images |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |