SAR基础知识课件(PPT 63页)

1、SAR基础知识技术支持邮箱：ENVI-技术支持热线：400-819-2881-5官方技术博客：1学习交流PPT第1页，共63页。1、SAR基础2、InSAR基础3、InSAR技术与地表形变主要内容2学习交流PPT第2页，共63页。1、SAR基础3学习交流PPT第3页，共63页。Radar概念Radar（Radio Detection And Ranging），一个Radar系统主要包括三个功能：发射微波信号到场景接收从场景中传回的部分后向散射能量观测返回的强度（检测）和延时（测距）信号RAR真实孔径雷达成像分辨率 雷达天线长度4学习交流PPT第4页，共63页。SAR（Synthetic Ape

2、rture Radar）：合成孔径雷达用一个小天线作为单个辐射单元，将此单元沿一直线不断移动，在不同位置上接收同一地物的回波信号并进行相关解调压缩处理的侧视雷达。可以获取高分辨率的地球表面图像，是目前广泛使用的雷达系统。SAR合成孔径雷达5学习交流PPT第5页，共63页。几个波的概念6学习交流PPT第6页，共63页。SAR参数波长雷达遥感使用的微波部分的电磁频谱，频率从0.3GHz至300 GHz，波长从1米到1毫米。微波中常用的波段有：P-band = 65 cm AIRSAR L-band = 23 cm JERS-1 SAR, ALOS PALSAR S-band = 10 cm Alm

3、az-1 C-band = 5 cm ERS-1/2 SAR, RADARSAT-1/2, ENVISAT ASAR, RISAT-1 X-band = 3 cm TerraSAR-X-1 , COSMO-SkyMed K-band = 1.2 cm 军事领域 波长越长穿透能力就越强，如波长大于2cm的雷达系统不会受到云的影响。7学习交流PPT第7页，共63页。SAR参数波长雷达频率的应用：冰雪识别，小型特征，使用X-band地质制图，大型特征，使用L-band叶面渗透，最好使用低频率，如P-band一般情况，C-band是折中波段波段名称的字母是早起军事应用遗留下来的8学习交流PPT第8页，

4、共63页。SAR参数极化极化：电磁波振动的矢量方向，即电场的方向9学习交流PPT第9页，共63页。SAR参数极化特点当雷达作用于地球表面时，其极化方式可能改变，产生随机极化反射信号，其中包含水平和垂直两种分量极化方式是否改变取决于目标的物理和电特性雷达可以接收反射信号的水平和垂直极化分量四种模式：HH、VV、HV、VH对于同一区域，不同极化方式获取的图像不同10学习交流PPT第10页，共63页。SAR参数极化美国航天飞机SAR获取的XVV、CHV、LHV的图像不同波长、 不同极化的散射机理差异11学习交流PPT第11页，共63页。SAR参数入射角(视角)入射角是雷达波束与垂直表面直线之间的夹角

5、（）。微波与表面的相互作用复杂，不同的角度区域会产生不同的回波。小入射角通常返回较强的信号，随入射角增加，返回信号逐渐减弱根据雷达距离地表高度的情况，入射角会随着近距离到远距离的改变而改变，从而影响成像几何。12学习交流PPT第12页，共63页。SAR参数入射角反射率随着入射角的增大而减小不同地物随入射角变化的反射率变化情况入射角与雷达后向散射的关系不同的应用选择相应装置的依据13学习交流PPT第13页，共63页。距离向(Range)：像平面内垂直于飞行方向斜距(Slant range)：雷达到目标的距离方向，雷达探测斜距方向的回波信号地距(Ground range)：将斜距投影到地球表面，是

6、地面物体间的真实距离方位向(Azimuth)：平行于飞行方向SAR的观测几何14学习交流PPT第14页，共63页。SAR图像几何分辨率距离分辨率Range侧视方向上的分辨率称为距离分辨率包括斜距分辨率和地距分辨率方位分辨率Azimuth沿航线方向上的分辨率，也称沿迹分辨率15学习交流PPT第15页，共63页。SAR图像几何特征雷达图像的构像几何学属于斜距投影类型斜距投影是以天线为中心，以斜距为半径的同心圆在像面空间上的投影方位向的比例尺是个常量距离向的比例尺由地面目标的位置由该目标到雷达天线的距离决定斜距影像地距影像SR=GR*sin16学习交流PPT第16页，共63页。SAR图像几何特征在雷

7、达成像中，地物目标的位置在方位向是按飞行平台的时序记录成像的在距离向上是按照地物目标反射信息的先后记录成像的在高程上即使微小变化都可造成相当大范围的扭曲这些诱导因子包括透视收缩、叠掩、阴影17学习交流PPT第17页，共63页。SAR图像几何特征：透视收缩透视收缩山上面向雷达的一面在图像上被压缩，这一部分往往表现为较高的亮度18学习交流PPT第18页，共63页。SAR图像几何特征需要提供高分辨率的DEM数据，透视收缩可通过几何校正和辐射定标进行校正用DEM校正不用DEM校正19学习交流PPT第19页，共63页。SAR图像几何特征：叠掩叠掩当面向雷达的山坡很陡时，出现山顶比山底更接近雷达，因此在图

8、像的距离方向，山顶与山底的相对位置出现颠倒20学习交流PPT第20页，共63页。SAR图像几何特征：叠掩叠掩21学习交流PPT第21页，共63页。SAR图像几何特征：阴影阴影地面上雷达信号照射不到的部分22学习交流PPT第22页，共63页。SAR图像几何特征在地形起伏的区域一般迎面坡是前向收缩坡度较大时，顶底叠置背面坡坡度较大时出现阴影阴影顶底叠置23学习交流PPT第23页，共63页。SAR散射机制雷达图像表示的是地面雷达后向散射的估算值主要可分为5种散射：表面和体散射双回波(Double Bounce)组合散射穿透散射介电属性散射。24学习交流PPT第24页，共63页。SAR散射机制表面和体

9、散射粗糙的表面能得到更高的后向散射，平整表面在雷达图像上经常表现暗区域。 双回波25学习交流PPT第25页，共63页。SAR散射机制组合散射一般发生在低频SAR系统（如L、P波段）包括表面、体散射、双回波等 森林的组合散射（上-林冠层，中-树干层，下-地面层）26学习交流PPT第26页，共63页。SAR散射机制穿透散射根据极化方式和波长情况，微波可以透入植被、裸土（干雪或沙地）一般情况，波长越长，穿透能力越强。交叉极化(VH/HV)相比同极化（HH/VV）的渗透能力弱。 27学习交流PPT第27页，共63页。SAR散射机制SAR穿透地表，可以清晰的看到沙漠下的地下河道28学习交流PPT第28页

10、，共63页。SAR图像的理解图像亮度代表后向散射强度像元内表面越粗糙，后向散射越强。光滑表面镜面反射，后向散射很弱与散射体的复介电常数有关含水量越大，后向散射越强29学习交流PPT第29页，共63页。水面上的溢油30学习交流PPT第30页，共63页。SAR图像的斑点噪声SAR是相干系统，斑点噪声是其固有特性均匀的区域，图像表现出明显的亮度随机变化，与分辨率、极化、入射角没有直接关系，属于乘机噪声多视和滤波可以抑制斑点噪声31学习交流PPT第31页，共63页。SAR数据类型单视复数（Single Look Complex-SLC）SAR数据由实部和虚部构成相位分布在-之间InSAR处理必须采用这

11、种数据振幅（Amplitude）、强度（Intensity）数据对单通道SAR系统来说，相位没有提供任何信息，振幅（或强度、能量）是唯一有用的信息In-Phase：同相分量Quadrature：正交相位分量32学习交流PPT第32页，共63页。SAR数据类型强度数据对经过聚焦处理的SAR数据进行多视处理，多视得到的强度图像是距离向和/或方位向像元分辨率的平均值一个L视数的图像，本质上是L的指数分布，33学习交流PPT第33页，共63页。星载SAR获取模式条带模式Stripmap当运行Stripmap SAR时，雷达天线可以灵活的调整，改变入射角以获取不同的成像宽幅。最新的SAR系统都具有这种成

12、像模式，包括RADARSAT-1/2， ENVISAT ASAR， ALOS PALSAR，TerraSAR-X-1， COSMOSkyMed和RISAT-1。34学习交流PPT第34页，共63页。星载SAR获取模式扫描模式ScanSAR 扫描模式是共享多个独立sub-swaths的操作时间，最后获取一个完整的图像覆盖区域。35学习交流PPT第35页，共63页。星载SAR获取模式聚束模式Spotlight 当执行聚束模式采集数据时，传感器控制天线不停向成像区域发射微波束。与条带模式主要区别为：在使用相同物理天线时，聚束模式提供更好的方位分辨率；在可能成像的以一个区域内，聚束模式在单通道上的提供

13、更多的视角；聚束模式可以更有效的获取多个小区域。36学习交流PPT第36页，共63页。系统发射时间波段极化图幅宽度（KM）分辨率重复周期轨道精度(cm)接受模式国家ERS-2(2011年退役)1995 CVV10025 m3530Stripmap欧洲RADASAT1（2013年4月出现故障）1995CVV10-50010-30-10024100StripmapScanSAR加拿大ENVISAT-ASAR（2012失去联系）2002CVV100-400203530StripmapScanSAR欧洲ALOS（2011已停止运行）2006LFull40-3507-14-10046100Stripma

14、p日本TerraSAR-XTandem-X20072010XFull5-10-30-1001-3-161110SpotlightStripmapScanSAR德国Cosmo-skymed2007X,LFull10-30-2001-3-151-1610SpotlightStripmapScanSAR意大利RADASAT22007CFull10-5003-1001-2410SpotlightStripmapScanSAR加拿大主要星载SAR系统37学习交流PPT第37页，共63页。主要星载SAR系统系统发射时间波段极化图幅宽度（KM）分辨率重复周期轨道精度(cm)接受模式国家ALOS-22014L

15、Full25/35/60/70/3501/3/6/10 /100 14100-日本哨兵”-1A(Sentinel-1A)2014cFull20/80/100/250/4005/20/4012-欧空局其他：RiSAT1（印度C波段）、Kompsat5（韩国X波段）38学习交流PPT第38页，共63页。SAR特性与光学遥感相比，SAR具有如下特性：全天候，不受云雾雪的影响，雨的影响有限全天时，主动遥感系统对地表有一定的穿透能力，与土壤含水量有关，依赖于波长对植被有一定的穿透能力，依赖于波长和入射角高分辨率，分辨率与距离无关独特的辐射和几何特性干涉测量能力多极化观测能力39学习交流PPT第39页，共

16、63页。雷达的应用雷达数据合成孔径雷达SARSAR图像基本应用雷达可以进行全天候观测，并且可以透过云层覆盖。自然应急救灾、农业估产、森林资源调查、军事应用、干旱监测InSAR（干涉SAR）合成孔径雷达干涉测量(InSAR)及差分InSAR(D-InSAR)技术是近十几年来发展非常迅速的微波遥感技术。由于它具有全天候、全天时、覆盖面广和高精度获取地表形变信息的能力。能达到厘米级精度，采用短基线干涉像对序列，能达到毫米级精度。InSAR两个方面应用：DEM提取和地表形变监测40学习交流PPT第40页，共63页。2、InSAR基础41学习交流PPT第41页，共63页。InSARSynthetic A

17、perture Radar InterferometryInSAR技术是利用雷达系统获取同一地区两幅SAR影像所提供的相位信息进行干涉处理，来获取地表的三维信息。InSAR（合成孔径）雷达干涉测量42学习交流PPT第42页，共63页。InSAR技术基本原理计算每次观测的相位，通过相位干涉求解相位差，进而反算地形及其地表形变信息。in SAR image #1:1 = R in SAR image #2:2 =( R + R) 43学习交流PPT第43页，共63页。InSAR测量获取DEM时实际上假设地表没有变化干涉相位组成44学习交流PPT第44页，共63页。InSAR技术获取地表形变信息（D

18、-InSAR）雷达两次不同位置获取同一监测区域的相位，差分干涉得到形变信息地表形变diff =Rmov+ N + A45学习交流PPT第45页，共63页。大范围DEM数据获取30/90米SRTM DEM，覆盖全球80%陆地Tandem-X 全球覆盖的5米高分辨率DEM数据获取高精度地表形变监测InSAR技术应用最为广泛且最能发挥优势的领域缓慢微小形变与突发形变大覆盖范围短周期、连续监测高精度、低成本InSAR技术应用46学习交流PPT第46页，共63页。1960s，InSAR概念提出；1974， Graham最早提出应用InSAR测量地形；1986，Goldstein等首先获得机载InSAR的

19、结果；1989，Gabriel等获得土壤表面形变结果；1993年，Massonet首先得到Landers地震的星载差分干涉结果；1995年，ERS2发射并与ERS1（1991）组合运行，提供了大量高质量InSAR数据，全球范围内掀起了InSAR研究热潮；1994、1996年SIR C都进行了重复轨道干涉试验；2000年，美国SRTM项目的成功；2001年，Envisat卫星发射；2006年，ALOS卫星发射；2007年，TerraSARX，COSMO SkyMed，RADASAT2，星载SAR进入米级分辨率时代。2010年，TanDEM-X，创新型雷达干涉仪InSAR技术发展历程47学习交流P

20、PT第47页，共63页。以前：专家系统开源软件非图形化界面.现在/未来工具软件简化操作.商业化软件的发展长弓，1-2年才能掌握滑膛枪，1-2月才能掌握自动来复枪，几小时内掌握48学习交流PPT第48页，共63页。3、InSAR技术与地表形变49学习交流PPT第49页，共63页。常规D-InSAR技术单时相差分干涉处理获取地表形变信息主要用于突发形地表形变，如地震形变场分析监测等干涉叠加技术基于时序分析方法，可以获得形变速率，达到毫米精度主要用于缓慢地表形变，如地面沉降主要包括PS、SBAS角反射器InSAR人工目标网络，形变体上不存在或存在较少相干目标的情况需要地表安装角反射器山体滑坡监测等，

21、空间监测范围较小，形变体变形幅度较小。地表形变测量的主要InSAR技术50学习交流PPT第50页，共63页。重复轨道InSAR测量DEM时实际上假设地表没有变化。实际上在发生地震、火山活动或者地壳运动的情况下，地表会有或大或小的形变。在InSAR技术的基础上，如果重复进行干涉成像或结合已有的精细DEM数据来消除干涉图中地形因素的影响，可以检测出地表的微小形变，这是D-InSAR的技术基础三种技术方法：双过差分、三过差分、四过差分。从可靠性上讲，双过差分干涉最可靠，而且目前全球大部分地区都有免费的SRTM3的DEM，可以满足很多应用需求。D-InSAR51学习交流PPT第51页，共63页。地形信

22、息来自已有的高精度DEM假设两幅SAR图像获取的时间段中存在SAR观测斜距方向的形变r ，干涉位相1可表示为D-InSAR的方法：双过（2pass）差分0是由DEM按照干涉基线和入射角模拟的位相差分干涉结果只与波长有关，与基线无关实际差分干涉处理时，干涉位相仍然经过去平，而模拟位相只需模拟地形位相即可。52学习交流PPT第52页，共63页。需要三幅SAR图像，其中两幅应该是在事件发生前，其中一幅作为公共主图像，其干涉结果得到本底DEM信息；公共主图像与事件后的图像进行干涉，再与得到DEM的干涉图进行差分，得到形变。其原理与二过差分干涉基本相同，但是处理流程有一些不同。D-InSAR的方法：三过

23、差分干涉53学习交流PPT第53页，共63页。四过差分干涉的数据要求较高，时空基线对适合的情况比较少。4幅SAR图像中3幅是事件发生前的，1幅是事件发生后的，但是无法找到公共主图像，只能分二组，一组生产DEM，一组包含事件形变信息。差分干涉流程类似于二过差分干涉，不同的是DEM是由另一组干涉生成。D-InSAR的方法：四过差分干涉SLC 1 dataSLC 2 dataSLC 3 dataSLC 4 data干涉图干涉图差分干涉相干变化54学习交流PPT第54页，共63页。技术特点算法简单，参数较少成本低，最少2景SAR图像即可可监测较大的地表形变，在LOS方向上能探测的SAR图像相邻分辨单元

24、之间的最大形变值不大于/2(如C波段为2.8cm，L波段为12cm)可连续、大范围监测可以监测突发性质的地表形变滑坡、地震、冰川移动、火山活动、矿区塌陷等制约因素相干条件时间基线和空间基线限制了可进行InSAR处理的数据量，以及地表覆盖物如浓密植被容易造成失相干轨道误差依赖于精密轨道数据大气影响非均一大气延迟对大区域处理的影响D-InSAR应用特点55学习交流PPT第55页，共63页。干涉叠加技术（Interferometric Stacking）D-InSAR技术适用于大范围地表形变监测单次测量对大的型变量监测有效，厘米级的（测量瞬间变化）时相上的失相干减少干涉相位精度大气影响解决方案干涉叠

25、加技术挖掘时间序列SAR图像可以识别区域（像素）范围内，一定时间内地面位移表现在信号相关和一致性，获取地表毫米级的形变信息。主要包括两个应用方向：PS (Persistent Scatterers，永久散射体)点目标特征，较合适监测SBAS (Small Baselines，短基线)分布式特征t56学习交流PPT第56页，共63页。PS技术将InSAR技术扩展到多时相的数据。假设研究区内有N+1副SAR图像，可得到N副相对于主图像的干涉图，再与DEM进行差分后生成N副差分干涉图，第K副图像差分相位为：PS-InSAR基本原理TBPS起源于2001年57学习交流PPT第57页，共63页。相干目标密集的地区城市、乡镇等植被覆盖地区适用性降低可以监测缓慢性质的地表形变城市地表沉降监测、线性交通工程（地铁、高铁等）、人工建筑（港口、机场）等数据量的要求大于20-30以上景等，基线要求尽可能小获取时间周期尽量连续，如一月一景测量精度形变速率1mmPS-InSAR应用特点58学习交流PPT第58页，共63页。连接由长基线造成的相