CN114002666A - 任意天线构型下星载ati-sar洋流流速提取方法及设备 - Google Patents
任意天线构型下星载ati-sar洋流流速提取方法及设备 Download PDFInfo
- Publication number
- CN114002666A CN114002666A CN202111248507.5A CN202111248507A CN114002666A CN 114002666 A CN114002666 A CN 114002666A CN 202111248507 A CN202111248507 A CN 202111248507A CN 114002666 A CN114002666 A CN 114002666A
- Authority
- CN
- China
- Prior art keywords
- channel
- velocity
- ati
- antenna configuration
- ocean current
- 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
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000003384 imaging method Methods 0.000 claims abstract description 30
- 238000012545 processing Methods 0.000 claims abstract description 23
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 20
- 238000012937 correction Methods 0.000 claims abstract description 17
- 238000004458 analytical method Methods 0.000 claims abstract description 11
- 238000012821 model calculation Methods 0.000 claims abstract description 7
- 230000006870 function Effects 0.000 claims description 22
- 238000010586 diagram Methods 0.000 claims description 16
- 238000000605 extraction Methods 0.000 claims description 6
- 230000035485 pulse pressure Effects 0.000 claims description 6
- 230000003044 adaptive effect Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000009826 distribution Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 238000012935 Averaging Methods 0.000 claims description 2
- 238000004590 computer program Methods 0.000 description 7
- 208000004350 Strabismus Diseases 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000004088 simulation Methods 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 230000001427 coherent effect Effects 0.000 description 2
- 238000002592 echocardiography Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000010587 phase diagram Methods 0.000 description 2
- 230000001902 propagating effect Effects 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 108010003272 Hyaluronate lyase Proteins 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011160 research Methods 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
- 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/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
- G01S13/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
- G01S13/589—Velocity or trajectory determination systems; Sense-of-movement determination systems measuring the velocity vector
-
- 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/9021—SAR image post-processing techniques
- G01S13/9023—SAR image post-processing techniques combined with interferometric 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/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/417—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section involving the use of neural networks
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)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种任意天线构型下星载ATI‑SAR洋流流速提取方法及设备，方法包括：根据参考区域通道A、通道B回波信号计算通道间幅度相位误差，对待处理目标区域进行补偿，实现通道校正；通过星载ATI‑SAR空间几何关系分析得到任意天线构型下的距离历程，对通过通道校正后的通道A、通道B回波信号采用BP成像算法得到两个通道的成像结果；将通道A、通道B的成像结果进行共轭相乘，获得场景回波的干涉相位；通过星载ATI‑SAR空间几何关系分析得到任意天线构型下的等效模型计算方法，根据该等效模型得到洋流混合速度；根据洋流混合速度提取得到洋流的径向速度。本发明可以实现任意天线构型下星载ATI‑SAR回波的洋流流速反演，且无需进行通道间的配准，提高处理效率。
Description
技术领域
本发明涉及合成孔径雷达对海遥感，具体涉及一种任意天线构型下的星载ATI-SAR洋流流速提取方法及设备。
背景技术
海洋流场是海洋动力环境最为重要的参数之一，利用干涉合成孔径雷达(Synthetic Aperture Radar，SAR)获取高分辨率、大测绘面积、高观测频率、高测量精度的海洋流场具有巨大的潜力和优势。利用干涉相位图可以直接计算得到雷达视线向速度分量后通过各种方法分离出流场速度，因此获得各种复杂构型下的干涉相位也是必不可少的部分。
同时对于两通道的斜视回波，采用常规成像算法(如wk算法、RD算法)时，目标成像在方位零频位置，导致两个通道方位向存在偏差以及距离向存在偏差，在提取干涉相位之前需要进行通道间的配准(频率相位补偿)。
发明内容
发明目的：为了解决上述背景技术提到的技术问题，本发明提出了任意天线构型下星载顺轨干涉SAR(Along-Track Interferometric-SAR，ATI-SAR)洋流流速提取方法及设备，可以实现任意天线构型下星载ATI-SAR回波的洋流流速反演，同时无需进行通道间的配准，提高处理效率。
技术方案：一种任意天线构型下星载ATI-SAR洋流流速提取方法，包括以下步骤：
根据参考区域通道A、通道B回波信号计算通道间幅度相位误差，利用通道间幅度相位误差对待处理目标区域进行补偿，实现通道校正；
通过星载ATI-SAR空间几何关系分析得到任意天线构型下的距离历程，基于距离历程对通过通道校正后的通道A、通道B回波信号采用BP成像算法得到两个通道的成像结果；
将通道A、通道B的成像结果进行共轭相乘，获得场景回波的干涉相位；
通过星载ATI-SAR空间几何关系分析得到任意天线构型下的等效模型计算方法，根据该等效模型得到洋流混合速度；
将洋流混合速度中的Bragg相速度与大尺度波轨道速度去除，得到洋流的径向速度。
进一步地，根据参考区域通道A、通道B回波信号计算通道间幅度相位误差包括：对参考区域通道A、通道B的回波信号进行二维频域自适应校正处理，获得通道A和通道B的通道间幅相误差。
进一步地，二维频域自适应校正处理包括：首先进行方位向迭代处理，然后进行距离向迭代处理，求解如下目标函数的最优解：
Z1(f,Ω)、Z2(f,Ω)分别表示通道A、通道B在二维频域的回波数据，f表示方位频域，Ω表示距离频域，h1,2(Ω)表示距离频域的补偿函数，D1,2(f)表示方位频域的补偿函数。
进一步地，通过星载ATI-SAR空间几何关系分析得到任意天线构型下的等效模型计算方法，根据该等效模型得到洋流混合速度包括：
基于天线的空间斜平面几何关系，构建新型的天线构型图，其中空间斜平面几何关系为：卫星沿y轴方向运动，速度为v，两个天线沿坐标平面原点呈点对称分布，距离为2d，且两天线之间的连线与x轴之间的夹角为θ，t＝0时刻运动目标和坐标中心的起始距离为R，且运动目标和原点的连线与x轴之间的夹角为γ，运动目标沿雷达视线向的速度为Vr；构建的新型的天线构型图为：新的x轴方向为目标与两个天线对称分布的原点之间的斜距方向，新的y轴方向为成像平面上与新的x轴垂直的方向；
等效模型表示为：等效基线长度为2d·sin(θ-γ)，等效卫星速度为v·cos(γ)，等效目标径向速度为
基于该等效模型的时间延迟表示为：
相应的干涉相位为：
根据干涉相位，得到洋流混合速度：
λ为雷达波长。
进一步地，Bragg相速度的近似计算公式为：
cp表示Bragg波相速度，g为重力加速度；G(θω)表示风向的扩展函数，θω为雷达视线向和风向的夹角，l为扩展因子。
本发明还提供一种计算机设备，包括：
一个或多个处理器；
存储器；以及
一个或多个程序，其中所述一个或多个程序被存储在所述存储器中，并且被配置为由所述一个或多个处理器执行，所述程序被处理器执行时实现本发明的任意天线构型下的星载ATI-SAR洋流流速提取方法。
有益效果：本发明可以实现任意天线构型下星载ATI-SAR回波的洋流流速反演，同时也为复杂基线条件下(天线斜视及混合基线存在的情况)星载ATI-SAR洋流测速性能分析提供解决方法与途径。同时对于两通道的斜视回波，采用不同的成像算法，同一个目标最终的成像位置会有差异。采用wk算法时，目标成像在方位零频位置，导致两个通道方位向存在偏差以及距离向存在偏差，在提取干涉相位之前需要进行通道间的配准(频率相位补偿)。采用BP算法时，如果选择地面网格作为参考点，则在任何的通道夹角、雷达斜视角下，同一个目标在两个通道中成像后的方位位置和距离位置完全一致，无需进行配准，大大提高了实验效率。
附图说明
图1是本发明实施例的任意天线构型下星载ATI-SAR洋流流速提取方法总体流程图；
图2是本发明实施例的任意天线构型的几何关系图；
图3是本发明实施例的任意天线构型的等效几何关系图；
图4是本发明实施例的通道A的距离脉压后的回波；
图5是本发明实施例的通道B的距离脉压后的回波；
图6是本发明实施例的通道A的BP成像结果；
图7是本发明实施例的通道B的BP成像结果。
具体实施方式
下面结合附图对本发明的技术方案作进一步说明。
本发明提出了任意天线构型下星载ATI-SAR洋流流速提取方法，在实施例中，利用仿真的星载SAR回波数据进行验证和分析。仿真实验所用计算机CPU为[email protected]八核，内存为16GB。具体的SAR***仿真参数如表1所示。
表1
参照图1，任意天线构型下星载ATI-SAR洋流流速提取方法包括以下步骤：
步骤1，通道校正：以参考区域两个通道回波信号获得通道间幅相误差，对目标区域进行补偿。
首先，选取一片海岛作为参考区域，由于该区域海拔高度近似为0，且无径向速度，因此没有沿航迹和垂直航迹引入的混合干涉相位。然后，将通道A、通道B的回波信号进行二维频域自适应校正处理，获得通道A和通道B的通道间幅相误差。最后，将估计的通道幅度相位误差补偿到待处理目标区域，实现通道校正。通道A、通道B有时也称为通道一、通道二或第一个通道、第二个通道。
具体而言，步骤1可以通过以下方式实现：
假设Zi(f,Ω)(i＝1,2)表示两个通道在二维频域的回波数据，其中f表示方位频域，Ω表示距离频域，则Zi(f,Ω)(i＝1,2)可以表示为：
其中，a(f)表示两个通道间的共同幅相部分，h1(Ω)和h2(Ω)分别表示两个通道的传递函数，D1(f)和D2(f)表示两个通道间的双程天线方向图，f·ΔT表示由通道间隔引起的在方位频域的相位差，ΔT表示通道间的时延。
雷达回波数据变换到距离-多普勒域，这时***方位向相位误差主要包括两个部分，第一个部分是由通道间的间隔引起的线性相位部分，第二个部分是非线性相位部分，这是由不同通道间接收天线方向图以及天线后端通道传递函数的差异引起的。第一部分理论上可以表示为：
其中fm表示第m个多普勒单元对应的多普勒频率，d表示两个通道之间的间隔，V表示雷达平台的运动速度。
根据式(3)知ΔT可以表示为：
将通道A、通道B的回波信号进行二维频域自适应校正处理，二维频域自适应校正算法实质上就是求解下面的目标函数的最优解：
该目标函数可以理解为：通过寻找距离频域的补偿函数h1,2(Ω)和方位频域的补偿函数D1,2(f)，使得两个通道间的差异最小化。同时可以看出二维自适应校正算法认为通道误差在方位频域和距离频域是相互独立的。关于二维自适应校正算法详见参考文献“Gierull,C.H.Digital Channel Balancing of Along-Track Interferometric SARData.2003,42.”。
该目标函数的解可以表示为：
其中，上标符号*表示共轭操作，上述解可以进一步表示成迭代的结构：
其中，
表示第二个通道的数据进行n次迭代后的数据。
通过式(8)和(9)可以看出，二维自适应校正算法首先是针对参考区域的两通道数据进行方位向迭代处理，然后进行距离向迭代处理，如此反复。一般情况下经过3-4次迭代后就可达到收敛。通道校正的实质就是以通道一作为参考，通过迭代处理后，将通道二的数据校正到与通道一的数据完全相同，即补偿量为
之后针对目标区域的二维频域双通道数据点乘补偿量，实现通道校正。
步骤2，将BP成像算法引入星载ATI-SAR信号处理中：
通过星载ATI-SAR空间几何关系分析得到一种任意天线构型下的距离历程，基于距离历程对通过通道校正后的通道A、通道B回波信号采用BP成像算法得到两个通道的成像结果。
在本实施例中，步骤2采用如下优选方案实现：
(2-1)距离向压缩：将雷达回波脉冲信号通过匹配滤波器进行脉压处理；
(2-2)距离向插值：计算成像区域每一像素点与每一脉冲时刻雷达天线的斜距，通常可以使用频域补零后进行逆傅里叶变换的方法进行距离向插值计算；
(2-3)划分网格：将成像区域划分成网格，获得所有网格点目标，最终形成的SAR图像中每个像素代表一个网格；
(2-4)后向投影：计算雷达在每个方位时刻(发射脉冲的时刻)与每个网格点的距离Ruv，并计算出双程时延
以及所有网格点相对于最近距离参考点的延迟时间△t。Ruv中的u、v分别表示为成像区域距离向和方位向上的对应点，c表示为光速。
(2-5)相干叠加：利用每一网格点的延迟时间△t，通过插值计算所对应的回波值，与上一方位向该网格点的回波值进行相位补偿后叠加，补偿量为
λ表示为雷达波长。
BP算法参考了“时延-叠加”的思想，在雷达应用中，其对雷达接收天线接收到的回波信号进行距离向匹配，获取回波数据中包含的相幅信息，再通过IFFT进行逆傅里叶变换，获取收发天线组合的时延，最后累积信号相干相加得到目标函数。详见参考文献“基于BP的InSAR成像算法及多基线相位解缠算法研究[D].电子科技大学,2014.”。
步骤3，提取干涉相位：将通道A、通道B的成像结果进行共轭相乘，获得场景回波的干涉相位；
步骤4，提取任意天线构型下的洋流混合速度：通过星载ATI-SAR空间几何关系分析得到一种任意天线构型下的等效模型计算方法，根据该等效模型得到洋流混合速度。
在本实施例中，步骤4采用如下方法实现：
根据项目情况考虑一般的空间斜平面几何关系，可以得到新型的天线构型图。图2为任意天线构型的几何关系图，考虑斜平面几何关系，卫星沿y轴方向运动，速度为v。两个天线沿图中的原点呈点对称分布，距离为2d，且两天线之间的连线与x轴之间的夹角为θ。运动目标和坐标中心的起始距离为R(t＝0时刻)，且运动目标和原点的连线与x轴之间的夹角为γ。运动目标沿雷达视线向的速度为Vr。
由于采用了斜视处理，所以可以构建新型天线的等效平面。图3为任意天线构型的等效几何关系图。新的x轴方向即为目标和两个天线对称分布的原点之间的斜距方向，新的y轴方向为成像平面上与新的x轴垂直的方向(即图3中的Xc和Yc轴)。则等效基线长度为2d·sin(θ-γ)，等效卫星速度为v·cos(γ)，等效目标径向速度为
在上述的等效模型下，时间延迟可以表示为：
相应的干涉相位为：
上式中，Vr为雷达视线向的目标速度，
为雷达视线向投影到海面的速度，这个速度包含了海浪运动等速度的多普勒速度。
通过步骤3与式(11)，可知提取出来的洋流混合速度为：
步骤5，任意天线构型下的洋流径向速度分离：将洋流混合速度中的Bragg波相速度与大尺度波轨道速度去除后可以得到洋流的径向速度。
洋流混合速度可以表示为：
Vr＝Vc+Vwind+Vo+Vb (13)
上式中，Vc为海表面流场本身的速度，Vwind是海表面风场引起的流速，Vo是大尺度波轨道速度，Vb是Bragg相速度。其中，海表面风场引起的表面流速和海表面流场本身的速度之和被认为是所需要的海表面流速，即Vs＝Vc+Vwind。
为了获得海表面流场信息，必须去除大尺度波轨道速度和Bragg相速度。其中大尺度波轨道速度是一个快变过程，具有周期性，而且它的均值为0，因此可以通过空间平均来去除，但是洋流和Bragg相速度是慢变的过程，且处于同一量级，不易去除。
Bragg波相速度包含了两个方向相反的Bragg波，一个是接近雷达方向传播的Bragg波，另一个是远离雷达方向传播的Bragg波。顺轨干涉相位图获得的多普勒速度中包含的Bragg波相速度本质上就是这两个不同方向传播的Bragg波相速度的加权矢量和。Bragg波相速度的计算公式为：
其中，g为重力加速度，λ为雷达波长，θ为入射角，即前述斜平面几何关系中两天线之间的连线与x轴之间的夹角。两个不同方向传播的Bragg波所占的比重可以用风向的扩展函数表示为：
其中，θω为雷达视线向和风向的夹角，l为扩展因子，与雷达频率相关，通常取2-5。因此，多普勒速度中的Bragg相速度可以表示为：
如图4所示为基于BP算法的通道A的距离脉压后的回波，图5为基于BP算法的通道B的距离脉压后的回波，图6为通道A的BP成像结果，图7为通道B的BP成像结果，因此双通道目标所在处的干涉相位差值为-0.0736弧度。基于任意天线构型下的等效模型计算方法所得到的理论干涉相位差为0弧度。通过表2的五种复杂构型下的回波仿真，将理论值与实际值进行对比可以验证任意天线构型下的等效模型计算方法的有效性。
表2
根据本发明的另一实施例，提供一种计算机设备，包括：
一个或多个处理器；
存储器；以及
一个或多个程序，其中所述一个或多个程序被存储在所述存储器中，并且被配置为由所述一个或多个处理器执行，所述程序被处理器执行时实现如方法实施例中所述的任意天线构型下星载ATI-SAR洋流流速提取方法。
本领域内的技术人员应明白，本发明的实施例可提供为方法、***、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
最后应当说明的是：以上实施例仅用于说明本发明的技术方案而非对其保护范围的限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：本领域技术人员阅读本发明后依然可对发明的具体实施方式进行种种变更、修改或者等同替换，但这些变更、修改或者等同替换，均在发明待批的权利要求保护范围之内。
Claims (10)
1.一种任意天线构型下星载ATI-SAR洋流流速提取方法，其特征在于，包括以下步骤：
根据参考区域通道A、通道B回波信号计算通道间幅度相位误差，利用通道间幅度相位误差对待处理目标区域进行补偿，实现通道校正；
通过星载ATI-SAR空间几何关系分析得到任意天线构型下的距离历程，基于距离历程对通过通道校正后的通道A、通道B回波信号采用BP成像算法得到两个通道的成像结果；
将通道A、通道B的成像结果进行共轭相乘，获得场景回波的干涉相位；
通过星载ATI-SAR空间几何关系分析得到任意天线构型下的等效模型计算方法，根据该等效模型得到洋流混合速度；
将洋流混合速度中的Bragg相速度与大尺度波轨道速度去除，得到洋流的径向速度。
2.根据权利要求1所述的任意天线构型下星载ATI-SAR洋流流速提取方法，其特征在于，根据参考区域通道A、通道B回波信号计算通道间幅度相位误差包括：对参考区域通道A、通道B的回波信号进行二维频域自适应校正处理，获得通道A和通道B的通道间幅相误差。
3.根据权利要求2所述的任意天线构型下星载ATI-SAR洋流流速提取方法，其特征在于，二维频域自适应校正处理包括：首先进行方位向迭代处理，然后进行距离向迭代处理，求解如下目标函数的最优解：
Z1(f,Ω)、Z2(f,Ω)分别表示通道A、通道B在二维频域的回波数据，f表示方位频域，Ω表示距离频域，h1,2(Ω)表示距离频域的补偿函数，D1,2(f)表示方位频域的补偿函数。
4.根据权利要求3所述的任意天线构型下星载ATI-SAR洋流流速提取方法，其特征在于，所述目标函数的解的迭代形式表示为：
上标符号*表示共轭操作，
表示通道B的数据进行n次迭代后的数据。
5.根据权利要求1所述的任意天线构型下星载ATI-SAR洋流流速提取方法，其特征在于，参考区域为海拔高度近似为0且无径向速度的海洋区域。
6.根据权利要求1所述的任意天线构型下星载ATI-SAR洋流流速提取方法，其特征在于，通过星载ATI-SAR空间几何关系分析得到任意天线构型下的距离历程包括：
将雷达回波脉冲信号通过匹配滤波器进行脉压处理；
计算成像区域每一像素点与每一脉冲时刻雷达天线的斜距；
将成像区域划分成网格，获得所有网格点目标，最终形成的SAR图像中每个像素代表一个网格；
计算雷达在每个方位时刻与每个网格点的距离Ruv，并计算出双程时延tuv以及所有网格点相对于最近距离参考点的延迟时间△t；
利用每一网格点的延迟时间△t，通过插值计算所对应的回波值，与上一方位向该网格点的回波值进行相位补偿后叠加。
7.根据权利要求1所述的任意天线构型下星载ATI-SAR洋流流速提取方法，其特征在于，通过星载ATI-SAR空间几何关系分析得到任意天线构型下的等效模型计算方法，根据该等效模型得到洋流混合速度包括：
基于天线的空间斜平面几何关系，构建新型的天线构型图，其中空间斜平面几何关系为：卫星沿y轴方向运动，速度为v，两个天线沿坐标平面原点呈点对称分布，距离为2d，且两天线之间的连线与x轴之间的夹角为θ，t＝0时刻运动目标和坐标中心的起始距离为R，且运动目标和原点的连线与x轴之间的夹角为γ，运动目标沿雷达视线向的速度为Vr；构建的新型的天线构型图为：新的x轴方向为目标与两个天线对称分布的原点之间的斜距方向，新的y轴方向为成像平面上与新的x轴垂直的方向；
等效模型表示为：等效基线长度为2d·sin(θ-γ)，等效卫星速度为v·cos(γ)，等效目标径向速度为
基于该等效模型的时间延迟表示为：
相应的干涉相位为：
根据干涉相位，得到洋流混合速度：
λ为雷达波长。
8.根据权利要求1所述的任意天线构型下星载ATI-SAR洋流流速提取方法，其特征在于，大尺度波轨道速度通过空间平均来去除。
9.根据权利要求7所述的任意天线构型下星载ATI-SAR洋流流速提取方法，其特征在于，Bragg相速度的近似计算公式为：
cp表示Bragg波相速度，g为重力加速度；G(θω)表示风向的扩展函数，θω为雷达视线向和风向的夹角，l为扩展因子。
10.一种计算机设备，其特征在于，包括：
一个或多个处理器；
存储器；以及
一个或多个程序，其中所述一个或多个程序被存储在所述存储器中，并且被配置为由所述一个或多个处理器执行，所述程序被处理器执行时实现如权利要求1-9中任一项所述的任意天线构型下星载ATI-SAR洋流流速提取方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111248507.5A CN114002666B (zh) | 2021-10-26 | 2021-10-26 | 任意天线构型下星载ati-sar洋流流速提取方法及设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111248507.5A CN114002666B (zh) | 2021-10-26 | 2021-10-26 | 任意天线构型下星载ati-sar洋流流速提取方法及设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114002666A true CN114002666A (zh) | 2022-02-01 |
| CN114002666B CN114002666B (zh) | 2024-04-12 |
Family
ID=79924137
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111248507.5A Active CN114002666B (zh) | 2021-10-26 | 2021-10-26 | 任意天线构型下星载ati-sar洋流流速提取方法及设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114002666B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116577772A (zh) * | 2023-07-14 | 2023-08-11 | 无锡航征科技有限公司 | 流速测量方法、流速流量仪、计算机设备及存储介质 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2004019056A1 (en) * | 2002-08-23 | 2004-03-04 | Dstl | Interferometric synthetic aperture radar for imaging the ocean surface |
| CN103235304A (zh) * | 2013-03-26 | 2013-08-07 | 中国科学院电子学研究所 | 一种基于误差等效rd模型修正的sar几何校正方法 |
| CN104730518A (zh) * | 2015-03-30 | 2015-06-24 | 北京空间飞行器总体设计部 | 一种基于高斯拟合的雷达多普勒谱估计海面流场的方法 |
| CN107505614A (zh) * | 2017-07-20 | 2017-12-22 | 上海无线电设备研究所 | 一种基于ati‑sar图像幅相信息融合的目标检测方法 |
| CN108152817A (zh) * | 2017-11-20 | 2018-06-12 | 上海无线电设备研究所 | 流速测量中顺轨干涉sar的多普勒中心误差补偿方法 |
| CN110007302A (zh) * | 2019-03-27 | 2019-07-12 | 北京空间飞行器总体设计部 | 一种星载双天线斜视顺轨干涉sar洋流速度测量方法 |
| CN110823191A (zh) * | 2019-10-08 | 2020-02-21 | 北京空间飞行器总体设计部 | 混合基线双天线斜视干涉sar洋流测量性能确定方法及*** |
-
2021
- 2021-10-26 CN CN202111248507.5A patent/CN114002666B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2004019056A1 (en) * | 2002-08-23 | 2004-03-04 | Dstl | Interferometric synthetic aperture radar for imaging the ocean surface |
| CN103235304A (zh) * | 2013-03-26 | 2013-08-07 | 中国科学院电子学研究所 | 一种基于误差等效rd模型修正的sar几何校正方法 |
| CN104730518A (zh) * | 2015-03-30 | 2015-06-24 | 北京空间飞行器总体设计部 | 一种基于高斯拟合的雷达多普勒谱估计海面流场的方法 |
| CN107505614A (zh) * | 2017-07-20 | 2017-12-22 | 上海无线电设备研究所 | 一种基于ati‑sar图像幅相信息融合的目标检测方法 |
| CN108152817A (zh) * | 2017-11-20 | 2018-06-12 | 上海无线电设备研究所 | 流速测量中顺轨干涉sar的多普勒中心误差补偿方法 |
| CN110007302A (zh) * | 2019-03-27 | 2019-07-12 | 北京空间飞行器总体设计部 | 一种星载双天线斜视顺轨干涉sar洋流速度测量方法 |
| CN110823191A (zh) * | 2019-10-08 | 2020-02-21 | 北京空间飞行器总体设计部 | 混合基线双天线斜视干涉sar洋流测量性能确定方法及*** |
Non-Patent Citations (1)
| Title |
|---|
| 闫贺;朱岱寅;: "星载ATI-SAR***洋流测速要点分析与仿真", ***仿真学报, no. 10, 8 October 2018 (2018-10-08) * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116577772A (zh) * | 2023-07-14 | 2023-08-11 | 无锡航征科技有限公司 | 流速测量方法、流速流量仪、计算机设备及存储介质 |
| CN116577772B (zh) * | 2023-07-14 | 2023-09-12 | 无锡航征科技有限公司 | 流速测量方法、流速流量仪、计算机设备及存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114002666B (zh) | 2024-04-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104833972B (zh) | 一种双基地调频连续波合成孔径雷达频率变标成像方法 | |
| CN111443339A (zh) | 一种双基sar空变校正成像方法、装置、设备及存储介质 | |
| CN103543453B (zh) | 一种地球同步轨道合成孔径雷达干涉的高程反演方法 | |
| CN104698457B (zh) | 一种迭代曲面预测InSAR成像及高度估计方法 | |
| CN108051809A (zh) | 基于Radon变换的运动目标成像方法、装置及电子设备 | |
| CN103454638B (zh) | 一种圆迹合成孔径雷达三维层析成像方法 | |
| KR102151362B1 (ko) | 극좌표변환을 이용한 항공기기반 영상복원장치 및 이를 이용한 영상복원방법 | |
| CN110673143A (zh) | 一种子孔径大斜视sar俯冲成像的两步处理方法 | |
| US20200150265A1 (en) | Close-range microwave imaging method and system | |
| CN112415515B (zh) | 一种机载圆迹sar对不同高度目标分离的方法 | |
| CN103616688A (zh) | 一种提升三维干涉逆合成孔径雷达图像质量的方法 | |
| CN113589285A (zh) | 一种飞行器sar实时成像方法 | |
| CN102004250A (zh) | 基于频域展开的星机联合双基地合成孔径雷达成像方法 | |
| Yang et al. | A novel three-step image formation scheme for unified focusing on spaceborne SAR data | |
| CN114545411A (zh) | 一种基于工程实现的极坐标格式多模高分辨sar成像方法 | |
| CN108132466B (zh) | 一种机载阵列天线下视三维成像方法和*** | |
| CN111487614B (zh) | 基于子孔径的曲线航迹弹载sar波前重建成像方法及*** | |
| CN103308914B (zh) | 一站固定双基地干涉sar处理方法 | |
| CN113156435A (zh) | 一种基于嵌入式gpu的弹载sar前侧视时域成像方法 | |
| CN112558070B (zh) | 圆周扫描地基sar的频域成像方法及装置 | |
| CN114002666B (zh) | 任意天线构型下星载ati-sar洋流流速提取方法及设备 | |
| CN103344958A (zh) | 基于星历数据的星载sar高阶多普勒参数估算方法 | |
| CN103792534B (zh) | 一种基于先验相位结构知识的sar两维自聚焦方法 | |
| CN113671497B (zh) | 基于圆柱对称模型的单通道sar目标三维坐标提取方法 | |
| CN116559905A (zh) | 一种双基sar海面舰船运动目标无畸变三维图像重构方法 |
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 | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20230109 Address after: No. 29, Qinhuai District, Qinhuai District, Nanjing, Jiangsu Applicant after: Nanjing University of Aeronautics and Astronautics Applicant after: NANJING HENGDIAN ELECTRONICS CO.,LTD. Address before: No. 29, Qinhuai District, Qinhuai District, Nanjing, Jiangsu Applicant before: Nanjing University of Aeronautics and Astronautics |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |