ETM ＋数据的云掩膜产品生成算法.pdfVIP

第 2期 矿 山 测 量 NO．2 2015年 4月 MINE SURVEYING Apr．2015 doi：10．3969／j．issn．100l一358X．2015．02．09 LandsatTM／ETM+数据的云掩膜产品生成算法 郭 秋 ，刘贺春 (1．晋 中职业技术学院，山西 晋 中 030600；2．河南理工大学，河南 焦作 454150； 3．山西华晋岩土工程勘察有限公司，山西 太原 030021) 摘要 ：文 中介 绍 了美 国地质勘探 调查局 (UnitedStatesGeologicalSurvey：USGS)针 对 LandsatTM／ ETM+数据所开发的一种 自动化厚云检 测算法 (AutomaticCloudCoverAssessment：ACCA)。该算法 以TM／ETM+传感器的热波段数据为基础 ，并通过两个以热波段数据为处理核心的过滤步骤来识别 出最终的云像元。ACCA算法的主要特点是 能够针对厚云高亮低温的特性并结合 多个综合性特征 阈值来有效检测 出图像 中的云团，此外对于云边界廓线的定位精度也较 高。 关键词：厚云检 测；LandsatTM／ETM+；ACCA算法 中图分类号 ：P23 文献标识码 ：B 文章编号：1001—358X(2015)02—0023—03 遥感影像上的云类型大致可以分为薄云和厚云 ～ 5需转换为大气顶层 (TopofAtmosphere：TOA)反 两种。薄云虽然使 图像变得模糊 ，但仍包含一定的 射率数据 ，波段6(热红外)需转换为星上温度数据。 有效信 息。厚云在可见光波段具有很高 的反射特 对于光谱波段 ，其图像量级水平 Q；与 TOA辐 性 ，同时在近红外波段吸收作用 明显，因而其主要的 射率 ；的关系可表达为 ： 物理特性是高反射与低温度 。厚云检测是利用厚云 Qi=Gi×Li+Q0i (1) 在反射光谱 、温度表现、形状等方面 的特征将 图像 中 其 中：Gi和 Q。分别为波段 i上传感器的响应度 的云像元进行有效标记。当前 的厚云检测方法按处 以及零辐射偏差。 理方式可大致分为三类 ，即基于热波段和多光谱数 波段 2～6可按下式将 TOA辐射率转换成 TOA 据的阈值提取技术…、多图像分析方法 I3以及 以 反射率 ： 分类技术为代表 的统计分析方法。多图像分析方法 Pi=仃 i ／(E0，cos0) (2) 的检测精度 尚佳 ，但需借助一幅或 多幅辅助参考数 上式 中：E。i为波段 i上的大气层外太 阳辐照度 ， 据且应经过一定 的相对辐射校正处理 ，如时间序 0是太阳天顶角 ，d 是太 阳一地球间的天文距离。 列分析法 ；基于分类技术的云检测方法主要是针 热红外波段 (波段 6)的星上温度 T可根据其 对那些不含热波段数据的高空间分辨率 图像发展起 TOA辐射率 L 得到 ： 来的，如基于最大似然估计。和模糊逻辑学 的方 T=K／ln((K ／L：)+1) (3) 法等 ，但这类方法的整体检测精度较低 ；由于厚云的 其中：K，、K2为定标常量。 高亮低温特性 ，基于热波段和多光谱数据 的阈值提 1．2 ACCA算法的处理步骤 取方法一般具有很高的检测精度 ，这类方法以ACCA ACCA算法综合使用 了5个光谱波段 的数据来 算法 为代表 。本文即在介绍 A