一种机载合成孔径雷达多普勒中心偏移计算方法技术

本发明专利技术公开了一种机载合成孔径雷达多普勒中心偏移计算方法，首先读取SAR探测回波数据，进行方位向FFT，获取多普勒功率谱；然后对多普勒功率谱进行小波分析滤波，确定基于回波数据的多普勒中心频率；接着读取机载SAR惯性导航数据和雷达性能参数，计算基于惯导数据的多普勒中心频率；最后根据基于回波数据的多普勒中心频率和基于惯导数据的多普勒中心频率，确定多普勒中心偏移。本发明专利技术计算方法简单有效，提高了多普勒中心偏移的计算精度，为未来的海洋环境应用提供了有利条件。

A method for calculating Doppler center offset of airborne SAR

The invention discloses a method for calculating Doppler center shift of airborne synthetic aperture radar. Firstly, SAR detection echo data are read, azimuth FFT is performed to obtain Doppler power spectrum, then the Doppler power spectrum is filtered by wavelet analysis to determine Doppler center frequency based on echo data, and then the airborne SAR inertia is read. The Doppler center frequency based on INS data is calculated from the data of sexual navigation and radar performance parameters. Finally, the Doppler center offset is determined according to the Doppler center frequency based on echo data and the Doppler center frequency based on INS data. The calculation method of the invention is simple and effective, improves the calculation accuracy of Doppler center migration, and provides favorable conditions for future marine environmental applications.

【技术实现步骤摘要】
一种机载合成孔径雷达多普勒中心偏移计算方法
本专利技术涉及微波遥感技术，特别是涉及一种机载合成孔径雷达多普勒中心偏移计算方法。
技术介绍
在合成孔径雷达(SyntheticApertureRadar,SAR)探测数据中，多普勒中心频率代表探测回波的相位信息，直接影响图像信噪比、方位模糊信号比以及成像后的目标位置，主要应用于雷达成像领域。近几年，随着搭载平台测量设备测量速度和精度的提高，人们在处理SAR探测回波数据时发现，平台运动状态数据计算的多普勒中心频率和海洋探测回波数据计算的多普勒中心频率始终存在特定差异，并开始认识到这种差异可能来自于海洋表层运动的影响。基于这一假设，学者们将SAR海洋探测回波数据和平台运动状态数据计算多普勒中心频率存在差异定义为多普勒中心偏移(DopplerCentroidShift，DCS)，进而开展海面风场、海表流场等海洋环境参数反演研究。目前，在多普勒中心频率计算研究方面，针对平台运动状态数据，Raney等和Bezvesilniy等[BezvesilniyOO,VavrivDM.Syntheticapertureradarsystemsforsmallaircrafts:Dataprocessingapproaches[M].INTECHOpenAccessPublisher,2012,p465-498.]分别提出了星载、机载SAR多普勒中心频率计算方法。针对探测回波数据，研究人员按照探测回波数据处理的不同阶段，提出了多种多普勒中心频率计算方法[HerlandEA.SeasatSARprocessingattheNorwegianDefenceResearchEstablishment(designandoperationofadigitalprocess)[J].ESAAppl.ofRemoteSensingDataontheContinentalShelfp247-253(SEEN82-1455305-43),1981.]。但是针对多普勒中心偏移计算方法的研究尚比较少，只有Hansen等利用EnvisatASAR探测数据计算了多普勒中心偏移，并对处理误差进行分析[HansenMW,CollardF,DagestadK,etal.RetrievalofseasurfacerangevelocitiesfromEnvisatASARDopplercentroidmeasurements[J].IEEETransactionsonGeoscienceandRemoteSensing,2011,49(10):3582-3592.]。然而，该方法主要针对星载SAR，不适用于机载SAR多普勒中心偏移的计算，且存在着计算精度低的不足。
技术实现思路
本专利技术所解决的技术问题在于提供一种机载合成孔径雷达多普勒中心偏移计算方法。实现本专利技术目的的技术解决方案为：一种机载合成孔径雷达多普勒中心偏移计算方法，包括以下步骤：步骤1、读取SAR探测回波数据，进行方位向FFT，获取多普勒功率谱；步骤2、对多普勒功率谱进行小波分析滤波，确定基于回波数据的多普勒中心频率；步骤3、读取机载SAR惯性导航数据和雷达性能参数，计算基于惯导数据的多普勒中心频率；步骤4、根据基于回波的多普勒中心频率和基于惯导数据的多普勒中心频率，确定多普勒中心偏移。本专利技术与现有技术相比，其显著优点为：1)本专利技术利用基于回波的多普勒中心频率和基于惯导的多普勒中心频率计算多普勒中心偏移，方法简单有效；2)本专利技术在多普勒谱分析中加入小波分析去除噪声的影响，提高了计算精度；3)本专利技术解决了机载SAR的多普勒中心偏移计算问题，为未来的海洋环境应用提供了有利条件。附图说明图1是本专利技术机载SAR探测几何关系图。图2是本专利技术多普勒中心频率的计算示例图，其中图(a)是一个脉冲序列的全极化回波数据排列，图(b)是读取的VV极化探测回波，图(c)是FFT后的探测回波，图是(d)多普勒功率谱。图3是本专利技术探测数据20120726-011的姿态信息图，其中图(a)是载机飞行探测的俯仰角、偏航角和横滚角，图(b)是载机飞行速度，图(c)是载机飞行高度。图4是本专利技术机载合成孔径雷达多普勒中心偏移计算方法的流程图。具体实施方式下面结合附图和具体实施例，进一步描述本专利技术方案。机载合成孔径雷达多普勒中心偏移计算方法，包括以下步骤：步骤1、读取SAR探测回波数据，进行方位向FFT，获取多普勒功率谱。如果目标场景均匀(如海洋、草地)，后向散射系数ξ(t,r)满足高斯分布，方位向回波信号可表示为：式中，t表示接收时间，r表示探测斜距，ha(t,r)是方位向系统响应函数。将上式变换到频域，方位向功率谱(又称多普勒功率谱)可表示为：式中，f表示频率，r表示探测斜距，ξ0表示后向散射系数，Ha(f,r)表示方位向系统响应函数。步骤2、对多普勒功率谱进行小波分析滤波，确定基于回波数据的多普勒中心频率。小波基函数通常定义为：式中，ψa,b(x)是小波母函数ψ(x)经过平移和缩放的结果，a为尺度因子，b为缩放因子，x表示自变量。小波函数选取是小波分析应用研究中重要的科学问题。SAR数据处理中，我们发现Mexican-Hat母小波函数具有任意阶的正则性，同时具有较好的时域和频域局部化特性及信号能量集中特性，有利于去除噪声对多普勒功率谱的影响。Mexican-Hat小波函数为高斯函数的二阶导数，可表示为：式中，ψ(x)表示小波母函数，x表示自变量。步骤3、读取机载SAR惯性导航数据和雷达性能参数，计算基于惯导数据的多普勒中心频率。机载SAR姿态变化时构成的探测几何关系如图1所示。载机沿着x轴以速度V在H高度平行于平面(xoy)开展探测飞行。雷达天线指向用矢量表示，为天线方向图与平面(xoy)的切线，α为俯仰角，β为偏航角，(xR,yR)为天线波束中心到探测目标的指向点。此时，天线波束中心到地面指向点的坐标为：则多普勒中心频率计算公式为：式中，λ表示雷达波长，V表示搭载平台飞行速度，r表示探测斜距，xR表示天线波束指向在x轴的投影，其中，H表示载机飞行高度，α表示俯仰角，β表示偏航角。步骤4、将基于回波的多普勒中心频率和基于惯导数据的多普勒中心频率作差，确定多普勒中心偏移：式中，fdcs表示多普勒中心偏移，表示基于回波数据计算的多普勒中心频率，表示基于惯导数据计算的多普勒中心频率。为了验证本专利的实用性，利用机载SAR探测数据开展计算效果检验。机载SAR的飞行探测时间为2012年7月26日14时到21时，载机飞行高度为3000m，探测分辨率为1.5m、3m，雷达中心入射角为30°、40°。惯性导航参数及性能指标、机载SAR工作参数及指标分别如表1,2所示。表1HJG-5惯性导航系统输出参数及性能指标表2机载SAR工作参数及指标机载SAR的探测回波为全极化数据，本实施例采用VV极化回波数据计算多普勒中心频率。本专利以探测数据20120726-011为例，介绍机载SAR多普勒中心偏移的计算方法。首先，利用探测回波数据计算多普勒中心频率。读取机载SAR探测回波数据，一个脉冲序列的全极化回波数据排列如图2(a)所示，图中幅度显著变化部分为回波信号，极化方式依次为HH、HV、VH和VV，其余本文档来自技高网...

【技术保护点】
1.一种机载合成孔径雷达多普勒中心偏移计算方法，其特征在于，包括以下步骤：步骤1、读取SAR探测回波数据，进行方位向FFT，获取多普勒功率谱；步骤2、对多普勒功率谱进行小波分析滤波，确定基于回波数据的多普勒中心频率；步骤3、读取机载SAR惯性导航数据和雷达性能参数，计算基于惯导数据的多普勒中心频率；步骤4、根据基于回波数据的多普勒中心频率和基于惯导数据的多普勒中心频率，确定多普勒中心偏移。

【技术特征摘要】
1.一种机载合成孔径雷达多普勒中心偏移计算方法，其特征在于，包括以下步骤：步骤1、读取SAR探测回波数据，进行方位向FFT，获取多普勒功率谱；步骤2、对多普勒功率谱进行小波分析滤波，确定基于回波数据的多普勒中心频率；步骤3、读取机载SAR惯性导航数据和雷达性能参数，计算基于惯导数据的多普勒中心频率；步骤4、根据基于回波数据的多普勒中心频率和基于惯导数据的多普勒中心频率，确定多普勒中心偏移。2.根据权利要求1所述的机载合成孔径雷达多普勒中心偏移计算方法，其特征在于，步骤1确定多普勒功率谱的公式为：式中，Sa(f,r)表示多普勒功率谱，f表示频率，r表示探测斜距，ξ0表示后向散射系数，Ha(f,r)表示方位向系统响应函数。3.根据权利要求1所述机载合成孔...

【专利技术属性】
技术研发人员：赵现斌，严卫，艾未华，陆文，王蕊，马烁，胡申森，陈冠宇，王宇训，谭仲辉，
申请(专利权)人：中国人民解放军国防科技大学，
类型：发明
国别省市：湖南,43