海岸带遥感09

1、第九章 合成孔径雷达（SAR）卫星遥感在海岸带环境监测中的应用 一、概述 二、SAR遥感研究现状与发展前景 三、SAR遥感技术在海洋以及海岸带的应用 一、概述（一）、合成孔径雷达 SAR一般由发射机、调制器、天线、接收机、数据记录和处理等系统组成。发射脉冲通过天线射向地面，接收地物后向散射信号（雷达回波），经处理生成SAR图像。 （二）、雷达波散射特性 1、后向散射系数 雷达波的极化后向散射系数 是单位面积上目标的pq极化 (p,q=v,h为垂直、水平极化)雷达散射截面，为无量纲系数或dB 。 0pq2、雷达图像强度 图像像元显示的强度 与对应地面单元的雷达后向散射系数 成直接的比例关系： 是

2、一个与距离有关的系统因子，它与发射功率、地面单元的大小及天线增益函数等有关。 ),(jiPji,0),()(),(0jijkjiP)( jk（二）、雷达波散射特性 3、雷达图像斑点 随机面元的雷达波散射是生成图像斑点的原因，表现为地表目标的不同单元所对应像元的图像色调有随机变化。 二、SAR遥感研究现状与发展前景 国际上主要星载SAR的情况美国欧空局俄罗斯日本加拿大型号SeasatSIR-ASIR-BSIR-CERS-1,2AlmazJERS-1Radarsat发射时间（年.月）1978.061981.111984.11994.041991.071995.041991.031992.02199

3、5.11任务期限3个月3a8a8a2a2a2a5a轨道倾角/()108385757或极轨98.5739898.6飞行高度/km794252250325800300570798-821卫星发射质量/kg100024001850014002750工作波段LLLL,C,XCSLC极化方式HHHHHH全极化VVHHHH视角/()204715601560239060351060距离分辨率/m25405817581730153018(4视)9100方位分辨率/m25(4视)40(6视)25(4视)25(4视)30(4视)153018(4视)9100测绘带宽/km1005010602010010075455

4、10天线尺寸/mm10.72.169.42.1610.72.16124.1101151.5122.2151.5波束操纵方式机械电扫机械电扫信号带宽/MHz196121215.551511.6,17.3或30.0发射功率/W100010001000L:3600,C:2235,X:1000480019000011005000发射脉宽/s33.430.430.4L:30.4,C:3337.123542重复频率/Hz146316401464182414641824120019001640172030001550169012701390目的及应用海洋陆地资源多参数观测资源 资源资源三、SAR遥感技术在海

5、洋以及海岸带的应用 SAR遥感图像可显示重要的海洋特征，例如:（1）大面积海域的海面波浪场,如大洋风浪与涌浪,浅海区波浪折射、绕射等;（2）海面风场、海流流场、锋面等海面运动;（3）显示大陆边缘、海山与海流切变区的海洋内波;（4）强海流区域和河口附近海域涡流;（5）内海湾潮流展示,从近岸到近海区域的水系分布;（6）浅海水下地形和浅海水深调制,海面起伏;（7）运动船舶产生的海面尾迹;（8）各类型的海冰及其变化;（9）海面溢油及海面油污染;（10）海岸带地物变迁;（11）风暴潮、热带气旋、降雨等; 三、SAR遥感技术在海洋以及海岸带的应用 由此能提供定量的海洋信息有: （1）海面风场（风速风向）、

6、海面流场定量重构；（2）浅海下地形重构与水深反演； （3）海面舰船识别；（4）海洋油污染，船漏油、海难、井喷、排污等突发性事件监测；（5）海洋水产、海洋运输作业与渔业安全；（6）海冰监测、海面冰山监测；（7）灾害性海况预报：风暴潮、巨浪、台风、海啸、暴雨等；（8）海岸带制图；海岸带资源调查、海岸线变迁及土地利用、滩涂与海湾水域面 积、河口动态及浅水地形监测等。特别是对常年云雾覆盖、植被茂密的热带海岸，微波SAR具有独特的遥感能力。 （一）、SAR图像预处理 SAR图像数据的预处理，主要是几何校正与斑噪滤除。1、几何校正 具体工作流程为： 输入试验区地形图，生成DEM； 输入成像雷达工作参数，生

7、成试验区的模拟雷达影像； 以模拟影像为参考坐标系，选取对应控制点与实际影像进行配准； 对配准后的影像实施正射投影校正。 校正精度与所用地形图的比例尺有关，如精度要求高，应采用较大比例尺的地形图。 （一）、SAR图像预处理 2、斑点噪声滤除 改善和滤除斑点噪声常用的方法大体分两类：一是成像前的多视平滑预处理；二是成像后的滤波技术。成像前的多视平滑处理会损害空间分辨率，现大都采用后者。 近年来，有不少斑点噪声分布概率模型以及滤波处理算法的研究。其中有均值滤波、中值滤波、Frost自适应滤波、改良K-均值自适应滤波等。 （二）、SAR图像海洋特征要素探测 SAR-C观测到的海浪黄海海面海浪 波段极化

8、图像CV VSIR-C020km（二）、SAR图像海洋特征要素探测 1、 海浪的遥感探测 （1）在SAR图像上选择两块子图像（如512512像元大小）的区域进行海浪谱反演；（2)在计算图像谱前，须对图像进行辐射校正，将SAR图像中给出的各点的灰度值转化为 ，不同来源的SAR图像辐射校正方法有所不同，建议按照数据提供机构推荐的辐射校正方法进行。（3）从SAR图像上提取波浪谱信息，先选择合适的子图像（如512512像元）， 进行二维快速傅里叶变换（FFT），得到SAR图像谱。 02、内波探测 SAR对内波的成像机制是内波对海表面的微尺度波动产生调制作用，使得海面粗糙度发生变化，从而在SAR影像上以

9、明暗相间的条纹显示出来。SAR内波成像主要包括三个物理过程：1、内波在传播过程中引起海表面流场发生辐聚和辐散的变化；2、变化的表层流场通过调制表面微尺度波，改变了海面粗糙度；3、海面微尺度波使雷达波发生Bragg散射，从而在SAR图像上形成明暗相间的内波条纹 。星载SAR的内波成像与SAR的极化方式、工作波段等参数密切相关，同时也受到内波的外在属性（例如：波长、波向）和内在特性（例如跃层深度、振幅）等参数的影响。 （二）、SAR图像海洋特征要素探测 黄海内波ERS-2 SAR影像（1997.7.23） 南海内波可见光遥感影像（SPOT-3,1996.7.24） 南海海面内波的L波段HH极化SI

10、R-C图像(1994.4.18) 南海海面内波的 波段极化图像LH HSIR-C(1)飞行方向照射方向北010kmAA SAR对内波成像的示意图 传播方向大陆架SAR图像海表面内波辐聚辐散21hh 2h1h21hh 21hh （1）对SAR图像进行几何校正，确定内波的位置，同时对SAR图像进行滤波和图像增强处理，改善内波图像质量；（2）从SAR图像测量得到两个由半日潮或全日潮产生的独立内波群的间距A值和单个内波中最亮点与最暗点的间距（或亮带中心点与暗带中心点的间距）D值；（3）将从SAR图像上测量的内波群的间距A带入式： ，计算内波波速，其中 为半日潮或全日潮周期。（4）将内波波速 、水深 和

11、g代入式： ，计算内波深度；（5）将内波波速 、水深 、D值和g代入式： ，计算内波振幅。 TACg/TgCh2/1222 12)4(ghCghghghggCh2/1222 2420432. 134gghCghgDgCh（二）、SAR图像海洋特征要素探测 （二）、SAR图像海洋特征要素探测 3、船舶目标的检测（1）静止船舶目标 静止船舶目标在SAR图像上可能是孤立点或小区域，在海面回波比较强的时候，区分它们和噪声斑点是较困难的，采用如下步骤可较好地区分它们。A 原始图像进行斑噪滤波处理；B 用恒虚警率（CFAR）阈值方法作分割，判断是否有适当大小的候选目标存在；C 用双阈值KSW法再次分割；D

12、 经过分割后，相邻的目标点被聚集成候选目标；E 对候选目标进行特征检测。可选择三个特征：目标的大小、强峰点占整个目标 的比例和所有船舶的对比度。F 对船舶目标进行参数提取，如几何参数包括传播的长度、宽度、朝向。 （2）运动船舶目标 运动船舶目标往往通过其尾迹来间接检测。根据尾迹的特征，可以估计船只的运动方向、速度甚至船只的类型、结构和吨位。尾迹成像条件取决于船体结构与运动速度、SAR系统参数（如船只到SAR平台速度及其到星下点的距离），以及海面风速风向、海流和波浪状况。 （二）、SAR图像海洋特征要素探测 由海洋卫星Seasat SAR观测的海面尾迹识别运动船只中国L波段机载SAR在山东烟台海

13、域观测的运动船迹烟台海域图像L-SA R飞行方向照射方向北01 km4、海冰信息提取 由Radarsat SAR 观测的渤海辽东湾海冰 2 0 004 0 06 0 08 0 01 0 0 01 2 0 01 3 6 102 0 04 0 06 0 08 0 01 0 0 0 1 2 0 0 1 4 0 0 1 6 0 0 1 8 0 0 2 0 0 0 2 2 0 0 2 4 0 0 2 6 0 0 2 7 7 0-1 8 .0 0-2 5 .5 0-3 3 .0 0区域 的图像AHH（二）、SAR图像海洋特征要素探测 5、海洋污染监测A SAR影像的预处理B 直接分析C 关联分析D 模型研究E 基于研究结果的结论 （二）、SAR图像海洋特征要素探测 6、浅海水深反演 定量的浅海水深计算与反演一般是在下垫海底与粗糙海面的物理数学模型基础上，根据辐射定标，计算归一化后向散射截面 ；通过 和雷达参数(波长、极化和入射角)等，求作用量谱 ；输入 、海面风速风向和初始流速流向等数据，利用特征线方法和迭代计算海表层流