CN102565772B - 基于sar子孔径序列图像的海洋动态信息提取方法 - Google Patents
基于sar子孔径序列图像的海洋动态信息提取方法 Download PDFInfo
- Publication number
- CN102565772B CN102565772B CN 201010594962 CN201010594962A CN102565772B CN 102565772 B CN102565772 B CN 102565772B CN 201010594962 CN201010594962 CN 201010594962 CN 201010594962 A CN201010594962 A CN 201010594962A CN 102565772 B CN102565772 B CN 102565772B
- Authority
- CN
- China
- Prior art keywords
- sub
- image
- grid
- aperture
- sar
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于SAR子孔径序列图像的海洋动态信息提取方法,涉及SAR技术,包括:(A)对SAR回波数据进行子孔径成像处理,得到子孔径序列图像;(B)对每个子孔径图像进行特征提取和网格划分后,对每两个相邻子孔径图像进行一次海面动态信息提取,得到每个网格的特征运动速度Vx(S,T),Vy(S,T);(C)再对多次提取的结果进行平均,减小误差和噪声,得海洋动态信息。本发明方法,简单、易行,适用范围广,可准确地提取海面的动态信息。
Description
技术领域
本发明涉及合成孔径雷达(SAR)技术领域,是一种基于SAR子孔径序列图像的海洋动态信息提取方法。
背景技术
合成孔径雷达(SAR)是20世纪50年代末发展起来的一种主动微波遥感手段,与其他遥感手段相比,有探测范围大、分辨率高、全天时、全天候、分辨率不随斜距变化等优势,在很多领域都得到了大量成功的应用。尤其在对海洋这样探测范围大、云雾干扰大的区域进行遥感探测时,SAR更是一个不可或缺的重要遥感手段。世界上第一个星载SAR——SeaSAT就是应用于海洋遥感的,后来陆续发射的多颗SAR卫星,包括ERS-1/2,RADARSAT-1/2,EnviSAT等都有大量成功的海洋遥感应用。我国正在规划的“海洋三号”卫星预计也将装备SAR载荷。
海洋场景最大的特点就在于它的复杂的动态特性,它既包含周期小至秒级的小尺度波浪振动,也包括周期长达几十分钟甚至数小时的内波、洋流等等。海洋遥感的一个重要目的就是提取海洋的动态信息。目前的SAR海洋遥感信息处理通常都是基于单幅SAR图像进行人工或机器检测和信息提取,形象的说就是用“照相机”而不是用“摄像机”来获取动态信息。基于单幅图像来提取目标动态信息需要对目标特性有很好的先验知识而且提取的动态信息还是比较有限,如测区的海面风速、水下地形环境、海水分层结构等。另一方面,经典的SAR成像理论是要求目标在合成孔径时间内不变的,而海洋很多动态信息的时变周期往往比合成孔径时间还短,目前海洋SAR数据的信号处理方法往往直接套用固定目标的SAR处理方法,这样很容易丢失一些快变信息。SAR的合成孔径时间取决于方位向天线长度、平台速度和电磁波长,随着SAR分辨率指标的不断提高,合成孔径时间也不断提高,例如目前机载SAR米级分辨率的L波段SAR的合成孔径时间已经长达30秒甚至更长,如果直接采用固定目标的SAR信号处理方法就会丢失时变周期小于30秒的动态信息,这已经不能满足数百米波长的波浪成像的要求了(深水海域500波长波浪振动周期约为18秒)。这也是机载L波段SAR图像虽然额定分辨率可以达到米量级甚至更高,但却往往无法看到百米以内波长波浪的主要原因。因此直接套用陆地目标的SAR信息处理方法来处理海洋动态目标,会丢失很多海面动态信息。
发明内容
本发明的目的是公开一种基于SAR子孔径序列图像的海洋动态信息提取方法,不用进行复杂的目标特性建模,而准确地提取海面的动态信息。
为达到上述目的,本发明的技术解决方案是:
一种基于SAR子孔径序列图像的海洋动态信息提取方法,其包括步骤:
(A)对SAR回波数据进行子孔径成像处理,得到子孔径序列图像;
(B)对每个子孔径图像进行斑点噪声抑制、特征增强和网格划分后,对每两个相邻子孔径图像进行一次海面动态信息提取,得到每个网格的特征运动速度Vx(S,T),Vy(S,T);假设有N个子孔径图像,则有N-1个运动速度提取结果。
(C)再对N-1个运动速度提取的结果进行平均,减小误差和噪声,得海洋动态信息。
所述的海洋动态信息提取方法,其所述步骤(A),其处理流程如下:
1)SAR原始信号经前期处理后,转换到方位向频域;
2)对方位向频域进行子孔径分割,得多数个子孔径,为:子孔径1、子孔径2……、子孔径N;
3)对步骤2)进行后续SAR成像处理,得多数个子孔径序列图像,为:子孔径图像1、子孔径图像2……、子孔径图像N。
所述的海洋动态信息提取方法,其所述每个子孔径图像对应的频宽为:
其中,v为雷达飞行速度,D为天线方位向长度,N为子孔径个数。
所述的海洋动态信息提取方法,其所述相邻子孔径图像之间的间隔时间为:
其中,R0为目标在合成孔径时间内的最近斜距,λ为波长。
所述的海洋动态信息提取方法,其所述步骤(B)中,对每个子孔径图像首先采用多视的方法对SAR图像斑点噪声进行抑制,然后根据目标特性进行特征增强,抑制海面背景中的无用信息,其中,特征增强所采用的方法需根据所关注的目标特性进行选择:(1)对内波、海浪的线状特征,采用小波域非线性增强、曲波域非线性增强的方法进行特征增强;(2)对强点目标舰船,采用灰度非线性增强方法进行强点增强,增强图像中的有用特征,抑制无关信息和噪声;然后对图像进行网格划分。
所述的海洋动态信息提取方法,其所述网格的大小需与关注的目标类型相适合,对于内波的大尺度特征网格为1km2左右,对于流场、波浪的中尺度特征网格为0.01km2左右,对于舰船小尺度特征网格为100m2左右。
所述的海洋动态信息提取方法,其所述步骤(B)中,对每两个相邻子孔径图像进行一次海面动态信息提取,包括:
(1)对每个网格到相邻子孔径图像中对应位置网格进行特征匹配,由于子孔径图像间隔时间短,图像间特征相关度很高,但位移很小,因此采用频域相位相关法的方法进行亚像元的高精度平移估计;
具体计算方法如下:
首先将两个网格图像进行二维傅里叶变换后进行共轭相乘,并求解缠相位:
其中F1(u,v)、F2(u,v)分别为第一和第二网格图像的二维傅里叶变换,u、V分别为二维频率,∠·表示求相位运算,unwrap(·)表示相位解缠运算;
则两个图像在x和y方向的平移量采用如下公式计算:
n遍历1到Ny得到Ny组位移量a、b,对这N组位移量进行平均,减少相位噪声的干扰;
(2)各网格的位移量除以子孔径图像时间间隔,得到每个网格中特征运动的速度:
其中,Vx、Vy分别为x、y方向速度,δx、δy为图像x、y方向像素分辨率,S,T分别为x、y方向的网格序号。
本发明方法,简单、易行,适用范围广,可准确地提取海面的动态信息。
附图说明
图1为合成孔径雷达(SAR)成像原理示意图;
图2为本发明的基于SAR子孔径序列图像的海洋动态信息提取方法采用子孔径成像形成SAR图像序列的原理示意图;
图3为本发明方法中的子孔径图像生成的信号处理流程示意图;
图4为本发明方法中的动态信息提取流程示意图;
图5为本发明的基于SAR子孔径序列图像的海洋动态信息提取方法流程图。
具体实施方式
要获取场景中的动态信息,最佳的方法就是获得同一场景在不同时刻的图像序列,通过序列图像的对比和处理可以明确图像中什么东西是动的什么东西是静止的,运动速度等信息也可以更准确地提取,而不用进行复杂的目标特性建模。目前常规的SAR图像判读和处理基本上都是基于单幅图像的,形象地说SAR目前都是“照像机”而不是“摄像机”,这就大大限制了SAR在海洋动态信息提取方面的能力和应用。
但是经过信号处理是可以将SAR从“照相机”变成“摄像机”的。SAR的原理可以理解为通过平台飞行和雷达按均匀时间间隔采样,雷达的天线中心在空间中就形成了一个虚拟的大阵列,每个采样点就构成了虚拟阵列中的一个“阵元”,SAR成像的过程就是对虚拟大阵列的各个“阵元”进行相位加权相加,形成不同的虚拟“窄方向图”,从而实现对方位向的高分辨成像,其原理示意图如下所示:
图1中各个实线三角形代表天线在每个采样位置的天线方向图,而虚线代表虚拟阵列形成的不同方向的窄天线方向图。
如果将虚拟的大阵列在飞行方向分为多个小阵列,每个小阵列各自形成虚拟方向图实现对方位向成像。由于各个小阵列对同一个目标成像,这样各个小阵列所得到的图像就存在时间差。这就是采用子孔径成像形成SAR图像序列的原理,图2是两个子孔径成像示意图。
如图2所示,地面上的每个目标都可以被子孔径1或者子孔径2观测到,但是两者的观测时间有一定差异,时间差为雷达从子孔径1的中心到子孔径2的中心所经历的时间。通过比较各个子孔径图像间的差异就可以提取出海面的动态信息。
图5为本发明的基于SAR子孔径序列图像的海洋动态信息提取方法流程图,包括:
1、子孔径图像的生成方法:
由于SAR信号的方位向多普勒瞬时频率与天线方位向视角有近似的对应关系,因此通过将SAR信号方位向多普勒频谱进行分割可以很方便地实现子孔径序列成像。合成孔径雷达与目标间距离可以表示为:
其中R0为目标在合成孔径时间内的最近斜距,v为雷达飞行速度,t为方位向时间,X0为目标方位向位置。一般情况R0>>vt-X0则瞬时多普勒频率可以表示为:
其中λ为波长。
从(2)可以看出目标的瞬时多普勒频率与雷达和目标的相对方位位置(vt-X0)成正比。因此通过对SAR方位向多普勒频谱进行分割可以很方便实现子孔径的分割。由于RD、CS、ωK等主流SAR成像算法都有将信号变换到方位向频域的步骤,因此在频域进行子孔径分割的步骤不会给成像带来太多额外的运算量。
如果将孔径分为N份,则进行进行方位向频域分割时每个子孔径图像对应的频宽为:
其中D为天线方位向长度,N为子孔径个数。这样对应的相邻子孔径图像之间的间隔时间为:
子孔径图像生成的信号处理流程如图3所示。
2、海洋动态信息提取方法
首先对所有子孔径图像进行多视处理,抑制斑点噪声。然后针对应用所关注目标的特性,对各个子孔径图像进行特征增强。例如关注内波、海浪等线状特征的可以用小波域非线性增强、曲波域非线性增强等方法,关注舰船等强目标的可以采用采用灰度非线性增强方法进行强点增强,增强图像中的有用特征,抑制无关信息和噪声。然后对第一幅图像中含有目标特征区域进行网格划分。网格的大小需要视具体应用所关注的目标类型而定,例如内波等大尺度目标网格可以为1km2右,流场等中尺度目标网格可以为0.01km2左右,舰船等小尺度目标网格可以为100m2左右。然后对每个网格到相邻子孔径图像中对应位置附近进行特征匹配,由于子孔径图像间隔时间短,图像间特征相关度很高,但位移很小,因此采用频域相位相关法的方法进行亚像元的高精度平移估计。
具体计算方法如下:
首先将两个网格图像进行二维傅里叶变换后进行共轭相乘,并求解缠相位:
其中F1(u,v)、F2(u,v)分别为第一和第二网格图像的二维傅里叶变换,u、v分别为二维频率,∠·表示求相位运算,unwrap(·)表示相位解缠运算。
则两个图像在x和y方向的平移量可以采用如下公式计算
n遍历1到Ny可以得到Ny组位移量a、b,对这N组位移量进行平均,减少相位噪声的干扰。
各网格的位移量除以子孔径图像时间间隔就可以得到每个网格中特征运动的速度:
其中Vx、Vy分别为x、y方向速度,δx、δy为图像x、y方向像素分辨率,S,T分别为x、y方向的网格序号。
每两个相邻子孔径图像做一次特征匹配后,共得到N-1次匹配结果。对N-1次匹配结果进行平均,可以使得运算结果的误差和噪声更小。
动态信息提取的流程图如图4所示。
实施步骤:
1.根据应用要求确定子孔径图像间的时间间隔,然后根据公式(4)计算要分割的子孔径数目。
2.按照图3所示的流程图进行子孔径成像处理,得到同一场景的子孔径序列图像。
3.按照应用要求对图像进行滤波和特征提取等预处理和网格划分。
4.按照图4所示流程图对每两个相邻子孔径图像进行一次海面动态信息提取,得到每个网格的特征运动速度Vx(S,T),Vy(S,T)。再对N-1次结果进行平均。
Claims (7)
1.一种基于SAR子孔径序列图像的海洋动态信息提取方法,其特征在于:包括步骤:
(A)对SAR回波数据进行子孔径成像处理,得到子孔径序列图像;
(B)对每个子孔径图像进行斑点噪声抑制、特征增强和网格划分后,对每两个相邻子孔径图像进行一次海面动态信息提取,得到每个网格的特征运动速度Vx(S,T),Vy(S,T);假设有N个子孔径图像,则有N-1个运动速度提取结果;
(C)再对N-1个运动速度提取的结果进行平均,减小误差和噪声,得海洋动态信息。
2.如权利要求1所述的海洋动态信息提取方法,其特征在于:所述步骤(A),其处理流程如下:
1)SAR原始信号经前期处理后,转换到方位向频域;
2)对方位向频域进行子孔径分割,得多数个子孔径,为:子孔径1、子孔径2......、子孔径N;
3)对步骤2)进行后续SAR成像处理,得多数个子孔径序列图像,为:子孔径图像1、子孔径图像2......、子孔径图像N。
3.如权利要求1或2所述的海洋动态信息提取方法,其特征在于:所述每个子孔径图像对应的频宽为:
其中,v为雷达飞行速度,D为天线方位向长度,N为子孔径个数。
4.如权利要求3所述的海洋动态信息提取方法,其特征在于:所述相邻子孔径图像之间的间隔时间为:
其中,R0为目标在合成孔径时间内的最近斜距,λ为波长。
5.如权利要求1所述的海洋动态信息提取方法,其特征在于:所述步骤(B)中,对每个子孔径图像首先采用多视的方法对SAR图像斑点噪声进行抑制,然后根据目标特性进行特征增强,抑制海面背景中的无用信息,其中,特征增强所采用的方法需根据所关注的目标特性进行选择:(1)对内波、海浪的线状特征,采用小波域非线性增强、曲波域非线性增强的方法进行特征增强;(2)对强点目标舰船,采用灰度非线性增强方法进行强点增强,增强图像中的有用特征,抑制无关信息和噪声;然后对图像进行网格划分。
6.如权利要求5所述的海洋动态信息提取方法,其特征在于:所述网格的大小需与关注的目标类型相适合,对于内波的大尺度特征网格为1km2左右,对于流场、波浪的中尺度特征网格为0.01km2左右,对于舰船小尺度特征网格为100m2左右。
7.如权利要求4所述的海洋动态信息提取方法,其特征在于:所述步骤(B)中,对每两个相邻子孔径图像进行一次海面动态信息提取,包括:
(1)对每个网格到相邻子孔径图像中对应位置网格进行特征匹配,由于子孔径图像间隔时间短,图像间特征相关度很高,但位移很小,因此采用频域相位相关法的方法进行亚像元的高精度平移估计;
具体计算方法如下:
首先将两个网格图像进行二维傅里叶变换后进行共轭相乘,并求解缠相位:
其中F1(u,v)、F2(u,v)分别为第一和第二网格图像的二维傅里叶变换,u、v分别为二维频率,∠·表示求相位运算,unwrap(·)表示相位解缠运算;
则两个图像在x和y方向的平移量采用如下公式计算:
n遍历1到Ny得到Ny组位移量a、b,对这N组位移量进行平均,减少相位噪声的干扰;
(2)各网格的位移量除以子孔径图像时间间隔,得到每个网格中特征运动的速度:
其中,Vx、Vy分别为x、y方向速度,δx、δy为图像x、y方向像素分辨率,S,T分别为x、y方向的网格序号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010594962 CN102565772B (zh) | 2010-12-20 | 2010-12-20 | 基于sar子孔径序列图像的海洋动态信息提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010594962 CN102565772B (zh) | 2010-12-20 | 2010-12-20 | 基于sar子孔径序列图像的海洋动态信息提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102565772A CN102565772A (zh) | 2012-07-11 |
CN102565772B true CN102565772B (zh) | 2013-07-24 |
Family
ID=46411661
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010594962 Expired - Fee Related CN102565772B (zh) | 2010-12-20 | 2010-12-20 | 基于sar子孔径序列图像的海洋动态信息提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102565772B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103675812B (zh) * | 2013-08-21 | 2016-02-03 | 中国科学院电子学研究所 | 一种基于圆迹sar数据的角度维散射信息提取及表征方法 |
CN104330796B (zh) * | 2014-10-20 | 2016-08-24 | 北京理工大学 | 一种基于子图像相干合成的地基合成孔径雷达快速成像方法 |
CN105093188B (zh) * | 2015-08-24 | 2017-05-24 | 中国科学院电子学研究所 | 基于子孔径图像的sar海洋图像中海浪杂波抑制方法 |
CN110147716B (zh) * | 2019-04-02 | 2021-02-19 | 北京理工雷科电子信息技术有限公司 | 一种基于频域与空域相结合的sar图像内波区域检测方法 |
CN110554377B (zh) * | 2019-09-05 | 2021-04-09 | 中国科学院电子学研究所 | 基于多普勒中心偏移的单通道sar二维流场反演方法及*** |
CN113514827B (zh) * | 2021-03-03 | 2023-09-05 | 南昌大学 | 无人机集群模式下合成孔径雷达成像处理方法及应用 |
US20230213648A1 (en) * | 2022-01-06 | 2023-07-06 | The Boeing Company | Radar processing method and system |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101140324A (zh) * | 2007-10-11 | 2008-03-12 | 上海交通大学 | 海域合成孔径雷达图像点目标提取的方法 |
CN101515036A (zh) * | 2008-02-20 | 2009-08-26 | 中国科学院电子学研究所 | 消除海面合成孔径雷达成像中时变模糊效应的方法 |
CN101893711A (zh) * | 2009-05-20 | 2010-11-24 | 中国科学院电子学研究所 | 基于智能体群体演化的海洋分布目标检测方法 |
US20110169686A1 (en) * | 2008-05-29 | 2011-07-14 | Gian Luca Eusebi Borzelli | Target detection in a sar-imaged sea area |
-
2010
- 2010-12-20 CN CN 201010594962 patent/CN102565772B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101140324A (zh) * | 2007-10-11 | 2008-03-12 | 上海交通大学 | 海域合成孔径雷达图像点目标提取的方法 |
CN101515036A (zh) * | 2008-02-20 | 2009-08-26 | 中国科学院电子学研究所 | 消除海面合成孔径雷达成像中时变模糊效应的方法 |
US20110169686A1 (en) * | 2008-05-29 | 2011-07-14 | Gian Luca Eusebi Borzelli | Target detection in a sar-imaged sea area |
CN101893711A (zh) * | 2009-05-20 | 2010-11-24 | 中国科学院电子学研究所 | 基于智能体群体演化的海洋分布目标检测方法 |
Non-Patent Citations (2)
Title |
---|
基于SAR的海洋内波参数的自动化检测算法与实现;樊慧; 王捷; 朱四华; 王子玲; 何高健;《海军航空工程学院学报》;20100531;第25卷(第5期);567-569 * |
基于上下层MAS的SAR图像海洋分布目标检测;胡希驰; 朱敏慧; 张永光; 王小青;《***工程与电子技术》;20091231;第31卷(第12期);2830-2834 * |
Also Published As
Publication number | Publication date |
---|---|
CN102565772A (zh) | 2012-07-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102565772B (zh) | 基于sar子孔径序列图像的海洋动态信息提取方法 | |
Chen et al. | Radon-fractional ambiguity function-based detection method of low-observable maneuvering target | |
Romeiser et al. | First analysis of TerraSAR-X along-track InSAR-derived current fields | |
CN102147469B (zh) | 一种双基地前视合成孔径雷达成像方法 | |
CN104950306B (zh) | 一种海杂波背景下前视海面目标角超分辨成像方法 | |
CN105445730B (zh) | 一种基于角度分集的海洋流场反演星载sar***及其方法 | |
CN103412310B (zh) | 双基地前视合成孔径雷达地面动目标检测方法与成像方法 | |
CN104833972B (zh) | 一种双基地调频连续波合成孔径雷达频率变标成像方法 | |
CN103207387B (zh) | 一种机载相控阵pd雷达杂波的快速模拟方法 | |
CN110568434B (zh) | 一种多通道匀加速sar动目标二维速度估计方法 | |
CN107576961B (zh) | 一种互质降采样间歇合成孔径雷达稀疏成像方法 | |
CN102749621B (zh) | 一种双基地合成孔径雷达频域成像方法 | |
CN103293521B (zh) | 一种利用x波段雷达探测近海海域水深的方法 | |
CN102914773B (zh) | 一种多航过圆周sar三维成像方法 | |
CN103091674A (zh) | 基于hrrp序列的空间目标高分辨成像方法 | |
KR102151362B1 (ko) | 극좌표변환을 이용한 항공기기반 영상복원장치 및 이를 이용한 영상복원방법 | |
CN104698457A (zh) | 一种迭代曲面预测InSAR成像及高度估计方法 | |
CN109061640B (zh) | 一种用于顺轨干涉sar海流反演的方位模糊抑制方法 | |
CN103091682B (zh) | 基于时频分析InISAR多动目标成像和运动轨迹重建法 | |
CN110554377B (zh) | 基于多普勒中心偏移的单通道sar二维流场反演方法及*** | |
CN108008387A (zh) | 一种机载阵列天线下视三维成像方法 | |
CN103630899A (zh) | 地面运动目标高分辨雷达压缩感知成像的方法 | |
CN109143236B (zh) | 适用于复杂飞行轨迹的双基聚束sar大场景成像方法 | |
CN106872979B (zh) | 基于海面起伏运动目标sar图像再聚焦的海浪参数获取方法 | |
CN107422320A (zh) | 一种消除降雨对x波段雷达观测海浪的影响的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130724 Termination date: 20181220 |
|
CF01 | Termination of patent right due to non-payment of annual fee |