CN105445730B - 一种基于角度分集的海洋流场反演星载sar***及其方法 - Google Patents
一种基于角度分集的海洋流场反演星载sar***及其方法 Download PDFInfo
- Publication number
- CN105445730B CN105445730B CN201510852862.1A CN201510852862A CN105445730B CN 105445730 B CN105445730 B CN 105445730B CN 201510852862 A CN201510852862 A CN 201510852862A CN 105445730 B CN105445730 B CN 105445730B
- Authority
- CN
- China
- Prior art keywords
- mrow
- mfrac
- radar
- msubsup
- ocean
- 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 32
- 208000004350 Strabismus Diseases 0.000 claims abstract description 11
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 230000008859 change Effects 0.000 abstract description 5
- 238000005259 measurement Methods 0.000 description 7
- 238000006243 chemical reaction Methods 0.000 description 4
- 230000002452 interceptive effect Effects 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 229930002875 chlorophyll Natural products 0.000 description 2
- 235000019804 chlorophyll Nutrition 0.000 description 2
- ATNHDLDRLWWWCB-AENOIHSZSA-M chlorophyll a Chemical compound C1([C@@H](C(=O)OC)C(=O)C2=C3C)=C2N2C3=CC(C(CC)=C3C)=[N+]4C3=CC3=C(C=C)C(C)=C5N3[Mg-2]42[N+]2=C1[C@@H](CCC(=O)OC\C=C(/C)CCC[C@H](C)CCC[C@H](C)CCCC(C)C)[C@H](C)C2=C5 ATNHDLDRLWWWCB-AENOIHSZSA-M 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000008020 evaporation Effects 0.000 description 1
- 238000001704 evaporation Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
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/9027—Pattern recognition for feature extraction
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
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)
- Computer Vision & Pattern Recognition (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于角度分集的海洋流场反演星载SAR***及其方法，该星载SAR***包括两个沿着卫星飞行方向放置的雷达天线，该两个雷达天线分别产生两个雷达波束，两个波束具有不同的斜视角且能够在不同的时间照射海面上的同一块区域。利用该星载SAR***以多普勒中心频率作为观测变量对海洋流场进行反演；得到海流的完整速度矢量，即反演得到的海流信息不仅包括海流的径向速度信息还包括方位向速度信息。相比于传统的海面流场反演手段，如现场观测、岸基雷达观测，本发明提出的基于角度分集的海洋流场星载SAR***能够反演全球范围的海洋流场，全面了解海流的时空变化；且能够全天时、全天候工作，不受气象条件限制。
Description
技术领域
本发明属于海洋遥感技术领域，具体涉及一种基于角度分集的海洋流场反演星载SAR***。
背景技术
海流对海洋上空的气候和天气的形成及变化以及海洋中多种物理、化学、生物以及地质过程有重要的影响和制约作用。比如，由低纬度流向高纬度的海流会造成水汽向上输送，使得空气湿度增大而产生降水；而由高纬度流向低纬度的海流会产生逆温，水汽不易向上输送，蒸发较弱而不易成雨。此外，海流在军事上也有重要的意义。对于海军作战，海流是作战考虑的重要因素之一，合理地利用海流会对作战起到降低战争消耗，增加胜利的筹码，反之，不考虑海流的影响会对战争起到不利影响，甚至导致失败。海流信息的获取在民用方面也有重要的作用，比如，贸易船只航线的选择要参考海流的状况等。因此，全球范围的海洋流场信息的反演无论是对于研究海洋本身的某些特征，还是与海洋密切相关的气候、天气，甚至军事应用均起着至关重要的作用。然而，目前海流是卫星遥感海洋动力参数中最不成熟的一个参数，还没有一个专门的卫星遥感测流业务化传感器。
海面流场的观测方法主要包括现场观测、岸基雷达观测、卫星遥感观测等方式。首先，测量海流的比较直接的一种方式就是“现场观测”。这种方法的突出优点是其具有较高的测量精度，然而，这种方法的空间覆盖范围有限，难以全面了解海流的时空变化，并且每次观测的成本较高，这些缺点使得“现场观测”不能满足实际应用的需求。其次，岸基雷达测量海流，主要包含X波段雷达和高频地波雷达。X波段雷达测量近岸3公里以内的流场，但它只有在海浪大于1米，且波浪较为均匀的地方才能获得相应的流场。高频地波雷达是专门用来测量沿岸海流的雷达设备，然而高频地波雷达测量的空间范围有限，不能进行全球观测。其三是卫星遥感测流技术，卫星遥感技术主要包括雷达高度计获取地转流、多时相遥感图像反演海面流场、合成孔径雷达(SAR)技术获取海面流场等。卫星高度计通常仅适用于测量大范围的地转流，且存在重复观测周期较长的问题；高度计接收到的雷达回波信号容易受到来自陆地回波信号的“污染”以及大气误差的影响，离海岸线50公里以内的高度计数据通常不能实用。多时相光学遥感图像反演海面流场也是卫星遥感海面流场的一种方法，该方法是利用遥感图像中的示踪物(如海表温度、叶绿素浓度等)在海流作用下的变化来测量海流。遗憾的是，海表温度数据以及叶绿素数据的获得受到“云层覆盖”的制约，并且这种方法对于“等温区域”或“等叶绿素区域”通常是失效的。星载SAR流场反演主要有三种方法。第一种方法即多时相SAR图像反演海面流场，其方法同多时相光学图像反演算法。第二种方法称为沿轨干涉法(ATI),它利用两个沿着轨道放置的天线，通过计算两天线的干涉相位即可获取径向海流信息。第三种方法利用单天线的多普勒中心偏移进行流场反演，这种方法的海流径向速度测量精度依赖于由卫星轨道和姿态测量数据计算得到的多普勒中心频率。对于上述第一种方法，它能克服光学图像受云雨等天气的影响，但同样存在与光学图像类似的困难，难以获取粗糙度较均匀海面的流场，业务化应用受到很大的限制。对于第二种和第三种方法，它们都只能获取径向流场，难以获取海面二维矢量流场。此外，ATI方法还容易受到“场景相干时间”等因素的制约，而“场景相干时间”直接与Bragg散射波浪有关系；当“场景相干时间”小于“干涉SAR图像对的时间间隔”时，ATI方法测流的精度将会明显变差。
发明内容
为了克服现有海流反演方法的缺点：比如，现有的星载SAR测流方法无法反演海流的完整二维速度矢量；“现场观测”方法空间覆盖范围有限，难以全面了解海流的时空变化；X波段雷达只有在海浪大于1米，且波浪较为均匀的地方才能获得相应的流场。同时，为了充分发挥星载SAR在海洋流场反演中的优势，比如全球范围，全天时、全天候工作的能力，不受气象条件限制等，本发明的目的是提供一种基于角度分集的海洋流场反演星载SAR***，该星载SAR***利用“多普勒中心频率”作为海流信息的“载体”，能够反演海流的二维速度矢量，即包括海流的径向速度同时也包括海流的方位向速度。
为实现上述目的，本发明采用的技术方案为：
一种基于角度分集的海洋流场反演星载SAR***，包括两个沿着卫星飞行方向放置的雷达天线，该两个雷达天线分别产生两个雷达波束，两个波束具有不同的斜视角且能够在不同的时间照射海面上的同一块区域。
两个雷达天线分别将接收到的回波信号依次传输到低噪声放大器、下变频器或解调器、射频模拟波束形成模块、模拟信号/数字信号转换模块；在射频模拟波束形成模块，两个天线所用的阵元加权系数不同，从而产生两个具有不同斜视角的波束。
两个雷达天线交替发射脉冲信号并交替接收回波信号，从而使得两天线接收的回波信号避免相互干扰。
一种利用上述的星载SAR***的反演海洋流场的方法，利用基于角度分集的海洋流场反演星载SAR***的两个雷达天线分别产生两个雷达波束，两个雷达波束具有不同的斜视角且能够在不同的时间照射海面上的同一块区域，从而实现角度分集；并以多普勒中心频率作为观测变量对海洋流场进行反演；反演流场时，同一区域的海流将会在两个雷达天线接收信号中表现为两个不同的多普勒中心频率，从而根据两个方程反演出海流的完整速度矢量，即反演得到的海流信息不仅包括海流的径向速度信息还包括方位向速度信息。
本发明的有益效果是：
1)本发明提出的基于角度分集的星载SAR***为一全新的星载SAR体制，在海洋环境探测领域将具有非常好的应用前景。
2)相比于传统的海面流场反演手段，如现场观测、岸基雷达观测，本发明提出的基于角度分集的海洋流场星载SAR***能够反演全球范围的海洋流场，全面了解海流的时空变化；且能够全天时、全天候工作，不受气象条件限制。
3)SAR-ATI方法以及单天线SAR多普勒中心偏移方法只能得到海流的径向速度而无法得到完整的二维速度矢量。本发明提出的基于角度分集的星载SAR***很好地解决了这一难点问题，能够同时得到海流的径向速度信息与方位向速度信息。
4)对于已有的SAR-ATI方法，通常要求“场景相干时间”大于“干涉SAR图像对的间隔时间”；当此条件不满足时，ATI方法的流场反演精度将会明显下降。比较起来，当利用本发明提出的基于角度分集的星载SAR***反演海面流场时，不受上述条件的限制。
5)相比于已有流场反演手段，本发明提出的SAR***的流场反演步骤简单易行、运算量较小：仅需要SAR原始数据的一阶相位信息。
附图说明
图1为基于角度分集的海洋流场反演星载SAR***
图2为角度分集的实现方案流程；
图3为海洋流场反演流程；
图4为角度分集星载SAR反演流场仿真结果。
具体实施方式
下面结合附图和实例进一步说明本发明的具体实施方式，但本发明的内容并不限于此。
如图1所示为一种基于角度分集的海洋流场反演星载SAR***，包括两个沿着卫星飞行方向放置的雷达天线，该两个雷达天线分别产生两个雷达波束，两个波束具有相同的下视角(图中的α)，但具有不同的斜视角，两个雷达波束能够在不同的时间照射海面上的同一块区域。两个雷达天线分别将接收到的回波信号依次传输到低噪声放大器、下变频器或解调器、射频模拟波束形成模块、模拟信号/数字信号转换模块。
下面根据图2详细说明实现角度分集的方案。如图2所示，对于第一个雷达天线(或第二个雷达天线)，每个天线阵元接收到的回波信号经过低噪声放大、下变频(或解调)后，接着进行射频模拟波束形成，最后进行模拟信号到数字信号的转换。射频模拟波束形成的作用是为了形成具有一定波束中心指向的天线方向图，其实现方式就是对各个阵元解调后的信号进行加权，然后求和。为了使第一个雷达天线与第二个雷达天线产生两个具有不同角度指向的波束，应使得各自的加权系数不同。如图2所示，两个具有不同中心指向的波束就构成了一个角度-分集SAR***，此时，该SAR***具有两个角度自由度，其中第一波束的斜视角为第二波束的斜视角为两个波束斜视的方向相反。
两个SAR天线交替发射脉冲信号并交替接收回波信号，从而使得两天线接收的回波信号避免相互干扰。
在观测变量的选取方面，本发明采用SAR多普勒中心频率这个变量。
下面根据图3说明利用角度-分集星载SAR反演海洋流场的原理与步骤。
a)首先将第一个雷达天线接收的SAR回波数据在慢时间(η)维平移如下时间量(见图1)：
其中，H为卫星平台高度，V为卫星平台的速度，α为两个波束的下视角，与分别为第一个SAR波束与第二个SAR波束的斜视角。
b)根据平均互相关系数(ACCC)法分别估计两天线接收信号的基带多普勒中心频率：
其中，f'Dc1与f'Dc2分别为第一个雷达天线与第二个雷达天线的基带多普勒中心频率；PRF为雷达脉冲重复频率；η表示慢时间(即方位时间)；s1(η)与s2(η)分别为第一个雷达天线与第二个雷达天线接收海面场景回波方位维原始SAR信号；∠{·}表示取信号的相位。
c)利用多视互相关算法(MLCC)分别估计第一个雷达天线与第二个雷达天线接收信号的多普勒模糊数Mamb1与Mamb2。
d)计算两天线接收SAR信号的绝对多普勒中心频率：
fDc1＝f'Dc1+Mamb1·PRF (4)
fDc2＝f'Dc2+Mamb2·PRF (5)
其中，fDc1与fDc2分别为第一个雷达天线与第二个雷达天线接收回波SAR信号的绝对多普勒中心频率。
e)从fDc1与fDc2中去除由***斜视所引起的多普勒频率，分别得到与
其中，λ为雷达波长。
f)海流的二维速度矢量，可以由下式计算得到：
其中，与分别为海流的方位向速度与距离向速度，H2×2为一2×2的矩阵，其表达式为：
其中，[]-1表示矩阵(方阵)的逆矩阵。
利用本发明提出的基于角度分集星载SAR***反演得到的海洋流场的仿真结果如图4所示。
综上，本发明有如下特点：
1)基于角度分集的海洋流场反演星载SAR***，用于反演海洋流场的完整二维速度矢量。采用角度分集星载SAR***并以多普勒中心频率作为观测变量对海洋流场进行反演，反演得到的海流信息不仅包括海流的径向速度信息还包括方位向速度信息。
2)利用两个雷达天线分别产生两个雷达波束；这两个波束“同时”照射海面，但它们的斜视角不同，从而实现SAR***的“角度分集”。
3)角度分集SAR***的两个波束的下视角相同，从而使得两个SAR波束能够在不同的时间照射海面上的同一块区域。
4)两个雷达天线交替发射脉冲信号并交替接收回波信号，从而使得两天线接收的回波信号避免相互干扰。
5)为了反演海洋流场，角度分集SAR***利用多普勒中心频率(即雷达波束中心经过某海面散射单元时对应的瞬时多普勒频率)作为观测变量，即作为海流信息的载体。
6)利用角度分集SAR反演流场时，同一区域的海流将会在两SAR天线接收信号中表现为两个“不同”的多普勒中心频率，从而可以根据两个方程反演出海流的完整速度矢量，即反演得到的海流信息不仅包括海流的径向速度信息还包括方位向速度信息。
以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
Claims (1)
1.一种利用星载SAR***的反演海洋流场的方法，其特征在于：利用基于角度分集的海洋流场反演星载SAR***的两个雷达天线分别产生两个雷达波束，两个波束具有不同的斜视角且能够在不同的时间照射海面上的同一块区域，从而实现角度分集；并以多普勒中心频率作为观测变量对海洋流场进行反演；反演流场时，同一区域的海流将会在两个雷达天线接收信号中表现为两个不同的多普勒中心频率，从而根据两个方程反演出海流的完整速度矢量；具体步骤如下：
(1)首先将第一个雷达天线接收的SAR回波数据在慢时间(η)维平移如下时间量：
其中，H为卫星平台高度，V为卫星平台的速度，α为两个波束的下视角，与分别为两个波束的斜视角；
(2)根据平均互相关系数法分别估计两天线接收信号的基带多普勒中心频率：
<mrow>
<msubsup>
<mi>f</mi>
<mrow>
<mi>D</mi>
<mi>c</mi>
<mn>1</mn>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<mi>P</mi>
<mi>R</mi>
<mi>F</mi>
</mrow>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</mfrac>
<mo>&angle;</mo>
<mo>{</mo>
<munder>
<mo>&Sigma;</mo>
<mi>&eta;</mi>
</munder>
<msubsup>
<mi>s</mi>
<mn>1</mn>
<mo>*</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&eta;</mi>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>P</mi>
<mi>R</mi>
<mi>F</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>}</mo>
<mo>,</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mi>P</mi>
<mi>R</mi>
<mi>F</mi>
</mrow>
<mn>2</mn>
</mfrac>
<mo>&le;</mo>
<msubsup>
<mi>f</mi>
<mrow>
<mi>D</mi>
<mi>c</mi>
<mn>1</mn>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>&le;</mo>
<mfrac>
<mrow>
<mi>P</mi>
<mi>R</mi>
<mi>F</mi>
</mrow>
<mn>2</mn>
</mfrac>
</mrow>
<mrow>
<msubsup>
<mi>f</mi>
<mrow>
<mi>D</mi>
<mi>c</mi>
<mn>2</mn>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<mfrac>
<mrow>
<mi>P</mi>
<mi>R</mi>
<mi>F</mi>
</mrow>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</mfrac>
<mo>&angle;</mo>
<mo>{</mo>
<munder>
<mo>&Sigma;</mo>
<mi>&eta;</mi>
</munder>
<msubsup>
<mi>s</mi>
<mn>2</mn>
<mo>*</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&eta;</mi>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>P</mi>
<mi>R</mi>
<mi>F</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>}</mo>
<mo>,</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mi>P</mi>
<mi>R</mi>
<mi>F</mi>
</mrow>
<mn>2</mn>
</mfrac>
<mo>&le;</mo>
<msubsup>
<mi>f</mi>
<mrow>
<mi>D</mi>
<mi>c</mi>
<mn>2</mn>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>&le;</mo>
<mfrac>
<mrow>
<mi>P</mi>
<mi>R</mi>
<mi>F</mi>
</mrow>
<mn>2</mn>
</mfrac>
</mrow>
其中，f′Dc1与f′Dc2分别为第一个雷达天线与第二个雷达天线的基带多普勒中心频率；PRF为雷达脉冲重复频率；η表示慢时间，即方位时间；s1(η)与s2(η)分别为第一个雷达天线与第二个雷达天线接收海面场景回波方位维原始SAR信号；∠{·}表示取信号的相位；
(3)利用多视互相关算法分别估计第一个雷达天线与第二个雷达天线接收信号的多普勒模糊数Mamb1与Mamb2；
(4)计算两天线接收SAR信号的绝对多普勒中心频率：
fDc1＝f′Dc1+Mamb1·PRF
fDc2＝f′Dc2+Mamb2·PRF
其中，fDc1与fDc2分别为第一个雷达天线与第二个雷达天线接收回波SAR信号的绝对多普勒中心频率；
(5)从fDc1与fDc2中去除由***斜视所引起的多普勒频率，分别得到与
其中，λ为雷达波长；
(6)海流的二维速度矢量，可以由下式计算得到：
其中，与分别为海流的方位向速度与距离向速度，H2×2为一2×2的矩阵，其表达式为：
其中，[]-1表示矩阵的逆矩阵。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510852862.1A CN105445730B (zh) | 2015-11-27 | 2015-11-27 | 一种基于角度分集的海洋流场反演星载sar***及其方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510852862.1A CN105445730B (zh) | 2015-11-27 | 2015-11-27 | 一种基于角度分集的海洋流场反演星载sar***及其方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105445730A CN105445730A (zh) | 2016-03-30 |
| CN105445730B true CN105445730B (zh) | 2017-09-15 |
Family
ID=55556148
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510852862.1A Active CN105445730B (zh) | 2015-11-27 | 2015-11-27 | 一种基于角度分集的海洋流场反演星载sar***及其方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105445730B (zh) |
Families Citing this family (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106301511B (zh) * | 2016-08-12 | 2019-10-29 | 上海卫星工程研究所 | 星载数传天线对地面接力传输***及其控制方法 |
| CN109283495B (zh) * | 2018-09-06 | 2022-05-06 | 南京信息工程大学 | 一种基于交叉谱分析的x波段雷达海流反演方法 |
| CN110109116B (zh) * | 2019-05-29 | 2021-03-30 | 国家卫星海洋应用中心 | 一种海洋表面流场值的确定方法、确定装置、设备和介质 |
| CN110221299B (zh) * | 2019-07-04 | 2023-02-10 | 南京信息工程大学 | 一种星载双通道双波束InSAR测流*** |
| CN110618403B (zh) * | 2019-08-26 | 2021-11-16 | 西安空间无线电技术研究所 | 一种基于双波束雷达的着陆飞行器参数测量方法 |
| CN110554377B (zh) * | 2019-09-05 | 2021-04-09 | 中国科学院电子学研究所 | 基于多普勒中心偏移的单通道sar二维流场反演方法及*** |
| CN110823191B (zh) * | 2019-10-08 | 2021-12-07 | 北京空间飞行器总体设计部 | 混合基线双天线斜视干涉sar洋流测量性能确定方法及*** |
| CN111487621B (zh) * | 2020-05-08 | 2023-06-13 | 宁波大学 | 一种基于雷达图像的海表流场反演方法及电子设备 |
| CN113064129B (zh) * | 2021-03-03 | 2023-06-16 | 湖北中南鹏力海洋探测***工程有限公司 | 一种高频地波雷达海流合成方法 |
| CN113608216B (zh) * | 2021-06-25 | 2024-04-09 | 航天恒星科技有限公司 | 一种星载多波段共口径sar及目标联合在轨检测***及方法 |
| CN113466854B (zh) * | 2021-06-29 | 2022-09-30 | 哈尔滨工业大学 | 基于海洋动力模型的高频地波雷达反演矢量流速方法 |
| CN114140697B (zh) * | 2021-09-02 | 2022-11-08 | 广东海启星海洋科技有限公司 | 海表流场遥感探测方法及装置 |
| CN114705884B (zh) * | 2022-06-06 | 2022-09-23 | 自然资源部第一海洋研究所 | 海洋风风向和海洋洋流流向监测方法 |
| CN115204058B (zh) * | 2022-08-02 | 2023-04-21 | 南方海洋科学与工程广东省实验室(珠海) | 基于bp神经网络的地波雷达流场计算方法及装置 |
| CN115291213B (zh) * | 2022-09-28 | 2022-12-20 | 中国电子科技集团公司第十四研究所 | 一种基于双波束sar的动目标方位速度估计方法 |
| CN116087949B (zh) * | 2022-12-13 | 2023-11-10 | 宁波大学 | 一种针对海岛海洋环境的sar海表流场反演方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102662164A (zh) * | 2012-03-20 | 2012-09-12 | 哈尔滨工程大学 | 基于x波段雷达图像和粒子群优化的海表面流信息提取方法 |
| CN103105603A (zh) * | 2013-01-25 | 2013-05-15 | 武汉大学 | 一种x波段测波雷达海流反演预处理方法 |
-
2015
- 2015-11-27 CN CN201510852862.1A patent/CN105445730B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102662164A (zh) * | 2012-03-20 | 2012-09-12 | 哈尔滨工程大学 | 基于x波段雷达图像和粒子群优化的海表面流信息提取方法 |
| CN103105603A (zh) * | 2013-01-25 | 2013-05-15 | 武汉大学 | 一种x波段测波雷达海流反演预处理方法 |
Non-Patent Citations (4)
| Title |
|---|
| "Mapping the Agulhas Current from space: An assessment of ASAR surface current velocities";M. J. Rouault et al.;《JOURNAL OF GEOPHYSICAL RESEARCH》;20101231;第115卷;第1-14页 * |
| "Retrieval of Sea Surface Range Velocities From Envisat ASAR Doppler Centroid Measurements";Morten Wergeland Hansen et al.;《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》;20111031;第49卷(第10期);第3582-3592页 * |
| "海面流场卫星遥感方法综述";何宜军 等;《广西科学》;20150630;第22卷(第3期);第3.1节，第4节，图4 * |
| "高海况海洋遥感信息提取技术研究进展";张彪 等;《海洋技术学报》;20150630;第34卷(第3期);第16-20页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105445730A (zh) | 2016-03-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105445730B (zh) | 一种基于角度分集的海洋流场反演星载sar***及其方法 | |
| CN110609287B (zh) | 一种双频雷达散射计及同时测量海面风场和流场的方法 | |
| CN110221299B (zh) | 一种星载双通道双波束InSAR测流*** | |
| CN101551450B (zh) | 一种星载极化sar法拉第旋转效应校正平台的构建方法 | |
| CN109507668B (zh) | 一种基于导航卫星信号的双基成像方法 | |
| CN101697011B (zh) | 双站合成孔径雷达海浪方向谱的仿真方法 | |
| JP2015180870A (ja) | マルチスタティックｆｍｃｗレーダーを用いるマイナス擬似レンジ処理 | |
| CN115336420B (zh) | 底视合成孔径与近底视双刈幅干涉成像的一体化探测装置 | |
| Shearman | Radio science and oceanography | |
| Li et al. | Ocean surface current extraction scheme with high-frequency distributed hybrid sky-surface wave radar system | |
| CN102565772B (zh) | 基于sar子孔径序列图像的海洋动态信息提取方法 | |
| CN109061640B (zh) | 一种用于顺轨干涉sar海流反演的方位模糊抑制方法 | |
| CN108008392B (zh) | 一种基于船载高频地波雷达的海洋表面风场测量方法 | |
| CN107607919B (zh) | 一种浮标基/船基高频地波雷达接收天线方向图测量方法 | |
| Romeiser et al. | Direct surface current field imaging from space by along-track InSAR and conventional SAR | |
| CN112162282A (zh) | 一种基于合成孔径雷达海表流速反演方法 | |
| CN107728144A (zh) | 一种基于前视双基模式的干涉sar成像技术 | |
| Farquharson et al. | An update on the capella space radar constellation | |
| CN110673128B (zh) | 一种基于间断上下调频波的x波段岸基雷达测流方法 | |
| Rosen et al. | A dual-frequency spaceborne SAR mission concept | |
| KR102289433B1 (ko) | 지향각 sar 기반 영상복원장치 및 이를 이용한 영상복원방법 | |
| Barrick et al. | Ocean surface features observed by HF coastal ground-wave radars: A progress review | |
| Kolpuke et al. | Airborne Multi-channel UWB FMCW Radar for Snow Depth Measurements | |
| Peters et al. | Measurements of ocean surface currents from a moving ship using VHF radar | |
| Maresca Jr et al. | High-frequency skywave radar measurements of hurricane Anita |
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 | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20200805 Address after: Room 2091, Heying building, No. 99, Tuanjie Road, yanchuangyuan, Jiangbei new district, Nanjing, Jiangsu Province Patentee after: Nanjing Xinda Environmental Technology Co.,Ltd. Address before: The Olympic Avenue in Jianye District of Nanjing city of Jiangsu Province, No. 69 210019 Patentee before: Nanjing University of Information Science and Technology |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20220630 Address after: 804, block B, No. 9, Shangdi Third Street, Haidian District, Beijing 100085 Patentee after: Zhongkexing Tuwei Tianxin (Beijing) Technology Co.,Ltd. Address before: Room 2091, hatch Eagle building, No. 99, Tuanjie Road, yanchuangyuan, Jiangbei new district, Nanjing, Jiangsu 211800 Patentee before: Nanjing Xinda Environmental Technology Co.,Ltd. |
|
| TR01 | Transfer of patent right | ||
| CP03 | Change of name, title or address |
Address after: Room 1801, 18th floor, building 1, yard 1, No. 81, Beiqing Road, Haidian District, Beijing 100094 Patentee after: Zhongkexing Tuwei Tianxin Technology Co.,Ltd. Address before: 804, block B, No. 9, Shangdi Third Street, Haidian District, Beijing 100085 Patentee before: Zhongkexing Tuwei Tianxin (Beijing) Technology Co.,Ltd. |
|
| CP03 | Change of name, title or address |