《基于单边拟合策略合成孔径雷达图像变化检测》由会员分享,可在线阅读,更多相关《基于单边拟合策略合成孔径雷达图像变化检测(33页珍藏版)》请在金锄头文库上搜索。
1、第一章 绪论1第一章 绪论1.1 SAR 成像基本原理及图像特点雷达主要是利用电磁波的散射进行目标探测和测距,具有全天候和强透射性等优点,使其在各个领域都发挥着重要作用。传统雷达主要解决目标有无的问题,获得回波脉冲及强弱,并不产生图像。后来人们利用多普勒锐化技术,实现高分辨率的雷达成像。SAR 是在雷达基础上发展起来的一种探测技术,它是一种利用雷达与目标的相对运动把尺寸较小的真实天线孔径用数据处理的方法合成较大的等效天线孔径的雷达。为了获得二维图像,在距离上采用线性调频信号技术,可以从其回波中将较近的目标分离出来,从而达到高分辨率;在方向上,通过平台的运动构建一个虚拟的大孔径天线,可以获得对另
2、一维的高分辨率1。合成孔径技术的基本思想是对在平台前进方向的不同位置上所接收的包含相位信息进行记录和处理,从而得到比实际天线更长的假设天线进行观测的结果2。SAR 获取图像的几何原理如图 1.1 所示。图 1.1 SAR 成像原理示意图成像处理一般包括三个步骤:1、距离向脉冲压缩,对一个合成孔径内的每第一章 绪论2个回波脉冲均进行压缩;2、距离徙动矫正,把属于同一点目标的脉冲校正到等距离线上、3、方位向脉冲压缩。难点在于方位向,需要利用回波信号的相位提取多普勒频率。平台轨道和姿态的不稳定及测量误差、雷达系统的不稳定等,均会引起信号的幅度和相位误差,从而影响匹配滤波的处理精度,很难实现理论的分辨
3、率(D/2) 。方位压缩的滤波匹配需要估计出两个参数:多普勒质心和多普勒频率。常用的算法有:距离-多普勒(Range-Doppler,R-D)算法。Chirp-Scaling(C-S)算法SAR 主要采用侧视雷达,以有源主动方式工作,其特殊的成像方式,形成了图像的独有特征,包括:辐射特征、噪声特征、目标集合特征、统计特征、纹理特征等。雷达图像的分辨率有距离向和方位向两种。雷达图像的分辨率和像元大小是两个概念,雷达图像中每个像元大小通常低于分辨率。SAR 图像的灰度统计特征:雷达图像上色调的变化,主要取决于目标物的后向散射截面。每一个接受到的回波被转换成电信号,并以某一特定的灰度色调记录在光学胶
4、片上或转换成一个具有特定值、由于表示亮度的数字化像元。坡度的变化,复介电常数(含水量) 、表面的粗糙度 是影响雷达图像色调的三个主要因素。其中表面的粗糙度在决定雷达图像的灰度(也就是回波强度)上起着决定性的作用。SAR 图像与传统的光学图像不同,它是通过天线发射有规律的无线电信号,并接收从地面返回的回波信号然后采取用成像处理的方法将回波信号转化为实际图片。雷达成像是依靠对回波信号提取目标的后向散射系数,反映所探测区域的微波特性,光学成像依靠的是普通的光学反射1。相对于与光学图像,SAR 图像具有许多重要特性:1 全天候、全天时,受风雨天气影响非常小,不受太阳照射影响,白天黑夜均可工作2 对电介
5、质(水分量、生物量、冰雪等) 、地表粗糙度(海洋风速等敏感)3 精确测量距离第一章 绪论34 对人工目标以及目标结构较敏感5 强透射性,能穿透遮盖物,比如植被、树叶6 分辨率不受雷达波长、载体飞行高度、雷达作用距离无关1.2 SAR 技术发展历程雷达卫星是载有合成孔径雷达(SAR)的对地观测遥感卫星的统称。尽管迄今为止,已在一些发射的卫星上携有 SAR,如 Seasat SAR, Almaz SAR, JERS-1 SAR, ERS-1/2 SAR, 与它们搭载在同一遥感平台上还装载着其他传感器。而 1995年 11 月发射的加拿大雷达卫星(Radarsat)则是一个兼顾商用及科学试验用途的雷
6、达系统,其主要探测目标为海冰, 同时还考虑到陆地成像,以便应用于农业、地质等领域。该系统有 5 种波束工作模式,即:(1)标准波束模式,入射角 2049 ,成像宽度 100 公里,距离及方位分辨率为 25 米 x28 米;(2)宽辐射波束,入射角 2040 ,成像宽度及空间分辨率分别为 150 公里和 28 米 x35 米;(3)高分辨率波束,三种参数依此为 3748,45 公里及 10 米 x10 米;(4)扫描雷达波束,该模式具有对全球快速成像能力,成像宽度大(300 公里或 500 公里) ,分辨率较低(50 米 x50 米或 100 米 x100 米) ,入射角为2049;(5)试验波
7、束,该模式最大特点为入射角大,且变化幅度小 4959,成像宽度及分辨率分别为 75 公里及 28 米 x 30 米。合成孔径的概念产生于 20 世纪 50 年代初期。1951 年 6 月,美国 Goodyear航空公司的 Carl Wiley 首先提出利用频率分析的方法来改善雷达角分辨率,并称其为多普勒波束锐化。1952 年 C.W.Shervin 提出了全聚焦陈列概念以及运动补偿概念。这些思想促进了 X 波段相干雷达的研制。1957 年密歇根大学雷达与光学实验室研制的 SAR 系统获取了第一张全聚焦的 SAR 图像。50 年代末,美国成功研制第一批可供军事侦察用的机载高分辨力合成孔径雷达。1
8、967 年 Greeberg 提出在卫星上安装 SAR 的设想。从此,人们开始了一航天飞机、卫星等为载体的空载 SAR 研究。60 年代中期,随着遥感技术的发展,合成孔径雷达技术从军用推广到民用方面,并发挥了重要作用。1978 年,NASA 发第一章 绪论4射了全球第一颗装载了合成孔径雷达的人造卫星 SEASAT,此卫星标志着 SAR 进入太空安装 SAR 对地观测的新时代,星载 SAR 由实验阶段转向实际应用的阶段。1981 年,美国“哥伦比亚”号航天飞机搭载 SAR-A 顺利升空以及 1984 年进行了“挑战者”航天飞机搭载 SAR-B 的实验。1987 年前苏联发射了第一部长期运行的空间
9、合成孔径雷达,名为 ALMAZ-I。1988 年 12 月美国用“阿特兰蒂斯”号航天飞机将世界上第一颗高分辨率雷达成像卫星“长曲棍球”(LACROSSE)SAR 卫星送入太空,其分辨率已达 1m,在海湾战争中发挥了重要作用。欧空局于 1991 年 7 月和 1995 年 4 月分别发射成功 ERS-1 和 ERS-2 卫星。日本于 1992 年发射了 JERS-1 空间合成空间雷达,工作于 L 波段。1994 年 9 月,美国航天局、德国空间局和意大利空间局共同进行了航天飞机成像雷达飞行任务。1995 年,加拿大发射了第一颗资源勘测卫星 Padarsat-1,该系统为商业和科学研究提供了全球的
10、冰情、海洋和地球资源等数据。SAR 成像技术随着时代不断前进发展,且功能实现更加完善,包括频率方面由单频段到多频段、成像空间由单极化到多极化,视角角度方面升级到多视角,性能的大大升级扩大了对全地球的覆盖率,提供的全日时的合成孔径雷达探测数据信息准确度更加准确,SAR 逐渐向小型化、SAR 星座组网及干涉测量等方向发展。目前合成空间雷达分辨率已达到 0.1m 数量级,纵观国外 SAR 的发展历程,可以看出随着科学技术的不断进步,SAR 的水平和功能也在不断提高。可以相信,科学家们将不断地挖掘 SAR 的技术潜力,为人类的需求服务。我国 SAR 研究开始于 20 世纪 70 年代后期,起步较晚。经
11、过几代人的不懈努力,许多技术已经走在世界前列23。1.3 SAR 图像变化检测介绍SAR 图像变化检测是指对不同时段同一地区的遥感图像,经过一定的分析和运算得出地面的变化区域、特征以及过程,涉及变化的类型、分布状况与变化量。SAR 图像变化检测的主要任务是检测局部文理和辐射值的变化。目前 SAR 图像变化检测思想主要分为两大类,一种是非监督的变化检测思路。常用非监督的变化检测方法常用的例如图像比值法4、图像差值法。此类方第一章 绪论5法优点在于简单、直接。如图 1.2 所示差,异图分类法的关键在于差异图的构造和阈值的选取。图 1.2 非监督的变化检测方法另一种主流思想是基于分类的变化检测方法思
12、路。如图 1.3 所示首先对不同时段的图像分别进行分类处理,然后根据分类结果比较,得出类别的不同并获得变化的区域。后分类比较法的关键在于多时相图像的准确分割和分割结果与差异图变化信息的提取与反馈。时相 1的图像时相 2的图像图像预处理图像预处理自动分类处理变化检测结果图分类图分类图比 较处 理自动分类处理图 1.3 基于分类的变化检测方法由于星载 SAR 遥感术后局具有全天候、全天时的独特优势,同事具有固定的重访周期,所以它比光学遥感更适合于变化检测。针对 SAR 图像的特点,目前比较通用的方法主要有差值法、相关系数法、直方图比较法、图像熵方法来对SAR 图像来进行变化检测2。归结起来,近些年
13、逐步形成的 SAR 图像变化检测方法可以分为:(1)基于简时相 1 的图 像时相 2 的图 像图像 预处 理图像 预处 理自动分 类处理变化检 测结果 图分类 图分类 图比较 处理自动分 类处理时相 1图像时相 2的图像图 像预处理图 像预处理比 较处 理变化检测结果图差 异图 像差异图像分析第一章 绪论6单代数运算的变化检测,经典的变化检测方法包括图像差值法、图像比值法、对数比值法;(2)基于图像变换的变换检测,经典的变化检测方法包括主成分分析、变化向量分析法、相关分析法图像变换法;(3)基于图像分类的变化检测方法。1.4 SAR 图像变化检测发展现状及应用国外在应用 SAR 遥感进行变化检
14、测方面起步较早,进行了一系列研究。Cyril 采用模糊隐马尔可夫链方法进行 SAR 变化检测5,通过计算两幅配准数据的对数比值,并引入模糊隐马尔可夫链,最大限度地提取相关性进行变化检测。Ramin6提出来通过模拟部分极化电磁场引入向量和观察量而建立变化空间,从而对相关或变化特征分解,更深入理解图像变化的本质。Andreas7在小波变换和位置上衣脊状特征重建原始图像,在曲波空间进行变化检测。Hachicha 和Chaabane8应用 Dezert-Smarandache 理论进行了变化检测的探索。Eric 采用 ERS-1 数据对 SAR 变化检测进行多种方法的初步研究和对比9。Onana 针对
15、 ERS-SAR多时相图像的特点,采用了数据融合的方法解决城镇地区的变化检测10。Patricia和 Kenneth 分别对 SAR 相关系数的变化检测方法进行了深入研究11,12。Francesca将小波变换的多尺度方法引入到 SAR 图像变化检测中13。Jordi 将聚合交叉熵计算方法引入到 SAR 图像变化检测中,取得了初步的结果13,14。在国内,由于 SAR 数据源不足的问题,开展 SAR 变化检测的研究相对数少,近几年来的主要成果有:基于二阶灰度统计特性的 SAR 图像变化检测15,该方法着重研究了统计特征中的二阶灰度特征如二维直方图、自相关、能量、协方差、惯性矩、熵、绝对值等,先
16、对多幅 SAR 图像分别计算灰度特征,然后用差值法对不同时相图像内相应的特征进行变化检测,最后用数学期望最大算法获得阈值。基于分形维数的 SAR 图像变化检测16,引入盒子维计算法,该方法最大的优点就是对 SAR 图像固有的斑点早说不敏感。也有学者引入马尔科夫随机场理论,研究了一种顾及空间领域关系的多时相 SAR 图像变化检测方法,该方法能有效提高漏检率,实现变化特征不明显区域的有效提取17。对于现在 SAR 图像变化检测方法大多根据 SAR 幅度的不足,有学者针对 SAR 图像的极化特征及其分布模型构建极化似然比检测模型18,实验证明该方法能有效地区分地物变化情况。针对传统分类比较法存在的分类误差累积问题,根据不同时相图像的不变信息的相关性,有学者提出了一种基于两时相图像联合分类的 SAR 图像变化检测方法19,第一章 绪论7该方法以灰度值作为输入信息,通过相似度计算得到两时相图像对应位置像素的灰度相似度,然后求解全局相似度阈值,并用 K 均值对两幅时相图像进行联合分