遥感影像云检测-传统方法-:Haze Optimized Transformation(HOT)方法

2024-01-30 19:18

本文主要是介绍遥感影像云检测-传统方法-:Haze Optimized Transformation(HOT)方法,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

An image transform to characterize and compensate for spatial variationsin thin cloud contamination of Landsat images

  • 1.BackGround
    • 1.1.存在问题
    • 1.2.现有方法
    • 1.3.本文贡献
  • 2.HOT算法
    • 2.1.基本思路
    • 2.2.问题
  • 3.对云和气溶胶类型的敏感性
  • 4.示例结果
  • 5.HOT在可见光波段应用的调整
    • 5.1.CL定义
    • 5.2.HOT影像生成
    • 5.3.直方图生成
    • 5.4.辐射平差估计
    • 5.5.图像调整
  • 说实话,没完全看懂

1.BackGround

1.1.存在问题

  • 即使在陆地卫星场景地理覆盖范围有限的情况下(例如185-185公里),感兴趣的时间窗口也很少有大气清晰的场景可用。
  • 对于朦胧的场景,通常缺乏辅助数据来进行绝对大气校正。
  • 空间变化的半透明云层和气溶胶层(以下简称“霾”)的污染是一个常见问题,影响到现有陆地卫星档案中的大部分场景。

光学深度(optical depth)指的是物质层不透明性的量度。即介质吸收辐射的能力的量度,等于入射辐射强度与出射辐射强度之比。

1.2.现有方法

- Dark Target Subtraction Methods
从直方图下界估计特定于波段的灰度偏移。然后在场景范围内应用这些偏移。可以通过将场景划分为子区域并估计每个子区域的直方图下限来扩展到场景内烟雾变化的建模
参考:An improved dark-object subtraction technique for atmospheric scattering correction of multispectral data
问题:精确的直方图下限估计需要大量像素,这种扩展只能用于模拟粗略尺度的烟雾结构(即,在几十公里的尺度上)。
改进: 结合DTSM和大气的物理特性.
参考:A dynamic regression algorithm for incorporating atmospheric models into image correction procedures.
**- 通过辐射测量变换分离出Haze **
在将Tasselled Cap(TC)变换扩展到TM传感器的过程中,烟雾主要显示在第四个TC分量中
问题:原始TC变换是专门为分离土壤和植被的表面辐射贡献(即“亮度”和“绿色”成分)而设计的,因此没有针对烟雾检测进行优化。
改进:开发了一种基于该TC变换的简化双波段版本(TM波段1和3)的气溶胶烟雾去除方法。

1.3.本文贡献

  • 提出一种稳健的算法,提取霾和云在图像上的空间分布
  • 提出一种利用烟雾变换图像对可见波段图像进行辐射调整的方法。对TM场景进行了处理和定量评估。

2.HOT算法

2.1.基本思路

找到一种变量可以量化大气对卫星图像辐射测量的影响,但同时对表面反射效应不敏感。必须满足:

  • 选择一个光谱空间(即光谱波段的混合),在该空间内,在晴朗大气条件下,不同土地覆盖类型的光谱响应高度相关。这将在光谱空间中形成定义良好的表面响应向量。此后,我们将此矢量称为**“清晰线”(CL)。**
  • 在选定光谱空间中,雾在不同波段的影响应当是不同的,因此当大气污染增加时,这个光谱空间的表示应当远离CL值
  • 合适的变换是其系数定义为与CL正交的方向,并且响应幅度与偏离该线成比例的变换。因为光谱空间中对地表覆盖和大气变化的光谱响应不是正交的,所以变换响应将仅与大气诱导的辐射成比例。

由两个可见波段组成的简单光谱空间满足前两个标准,例如,陆地卫星TM(以下简称TM1和TM3)的波段1和3:
以一组具有代表性的加拿大景观地表覆盖类型生成了TM1和TM3的晴空辐射估计值

TM1–TM3光谱空间示意图,在晴朗的天空条件下,编码为A–K的常见地表覆盖类型的辐射显示出高度相关性,并定义了“清晰线”(CL)。
数值序列1-18表示0.5μm处的云深度在0-6.7之间被均等分为18分,说明了光学深度增加的雾度对像素“迁移”离开CL的影响。HOT通过其在光谱空间中与CL的垂直距离量化像素位置处的大气污染水平。

可以看到,不同的地表覆盖对应的TM1TM3在图上的点呈现出强烈的相关性,也就存在定义明确的CL,当存在大气污染(即云雾霾),就会偏离这条线
HOT量化像素从该线的垂直位移的变换如下:
HOT=B1sinθ-B3cosθ

2.2.问题

大多数场景中大气条件的绝对信息将不可用。
解决:角度必须根据从视觉上被认为是最清晰的场景区域中选择的像素进行估计,并且在场景之间会有所不同。然后使用HOT将场景的其余部分调整为与这些最清晰区域相同的大气条件,
对于每个表面类,“云雾”所在点产生的轨迹随着云雾增加(光学深度增加)会从CL发生迁移,并且:

  • 恒定云线的斜率类似CL的斜率,光学深度类似于真实 CL 的斜率。如,图 1 中晴空 CL 的斜率是 1.03,而19 级云线的斜率为 0.98。 这意味着雾度的一阶辐射调整可以是使用适当缩放.
  • 对于上述选定波段,那些绝对反射率增加的表面,烟雾轨迹的斜率单调减小,说明所剔除的HOT对表面反射率存在一定的敏感度,或许需要二阶项对表面进行拟合
    如何解决表面反射率变化带来的影响?
    以HOT和TC转换的第四分量对应的相应水平来看(后一个波段是为了说明对光谱可见部分以外的波段应用热效应的可能性。),以下是模型的预测结果:
    TM1和近红外波段TM4的HOT()和TC(小写字母)值与DN水平的关系图。
  • 对于相同的雾度等级,不同表面覆盖等级的HOT是相似的,因此对于一阶而言,这种HOT转换与覆盖等级无关。
  • TC值表现出显著的变化,表明TC变换对表面反射率更为敏感。
  • 并且由于云层光学深度增加和地表反射率增加,TC响应趋于增加。这意味着在TC变换图像中,表面引起的增加可能会被误认为是雾。
  • 在变换(如HOT或TC)的实际应用中,必须指定一个阈值,低于该阈值时,不尝试进行雾度补偿,假设较低的变换后的值是雾度和表面成分的未分辨混合物。为了排除表面效应,TC的HAZE等效阈值必须远高于HOT。(因为上一条)
    从TM场景P21/R26的清晰区域中选择的像素的HOT和TC变换值与TM1 DN级别的散点图3。模型预测(见图2)表明,与TC变换相比,HOT对地表反射率明显不敏感,这一预测得到了真实图像的证实。
    在这里插入图片描述在这里插入图片描述

3.对云和气溶胶类型的敏感性

讨论HOT对不同类型烟雾的依赖性问题。为此,选择了一个表面类别(针叶林),并为其生成了各种云和气溶胶条件下的烟雾轨迹
在这里插入图片描述

比较MODTRAN模型预测的不同纯云和纯气溶胶类型的烟雾“轨迹”。对于五种仅含气溶胶的情况,数字绘制在数据点位置。对于相同的地表覆盖类型,这些轨迹几乎重合,这表明HOT以类似的方式对不同类型的大气污染作出响应是稳健的。另一方面,它不能用于独立估计光学深度。地表覆盖类型为针叶林。
VIS表示能见度。为了进行比较,还提出了一种没有云和气溶胶的理想大气(1型)。

对于四种类型的光学深度为0.55AM的大气颗粒(冰晶、水滴、城市气溶胶和农村气溶胶)(d0.55=1.27),由于散射和吸收特性的差异,每个TM波段中相应的辐射发生变化。TM可见光波段的视在辐射率随大气条件的变化而增加,依次为冰云、城市气溶胶、水云和农村气溶胶,和HOT的响应类似。在所有TM波段,冰云的反射比水云少,因此,在TM图像中,冰云比具有相同光学深度的水云显得更暗。在相同云层光学厚度下,冰云的HOT也低于水云。
此外,模型结果还表明,在云和气溶胶的低光学深度范围内,TM1和TM3 DNs增加,但TM4随着HOT响应的增加而保持不变。这意味着,至少对于针叶林上方光学非常薄的云层或气溶胶层,不需要对TM4数据进行辐射补偿。由于大气气溶胶的吸收和散射特性可能因气溶胶成分的不同而有很大差异,因此我们预计,表现出相同光学深度的不同气溶胶不会在给定TM波段引起相同的光谱辐射变化。

4.示例结果

为了测试HOT的有效性,我们将其应用于七个场景中,这些场景覆盖了典型的加拿大景观,涵盖了从晴朗到半透明薄雾到不透明云的各种大气条件。
场景及其CL描述符(即斜率和相关系数)述

七个TM测试场景的热参数摘要,包括CL斜率(SCL)、TM1和TM3变换系数(分别为sinHand cosH)和相关系数rCL

在这里插入图片描述

JC:急流尾迹和尾迹诱导云,Ci:卷云,St:薄层云或雾,Cu:积云,Mi:雾或薄卷云。A–M的表面类型与表1中的表面类型相对应。

在这里插入图片描述

a)薄薄的薄雾层覆盖了同一图像左侧的大片区域。图像的右侧部分看起来最清晰,该区域被用作晴空区域,用于定义HOT的“CL”。根据粗略的检查,虽然总体上看是清晰的,但HOT图像和更详细的视觉分析显示,场景左侧存在薄薄的薄雾层。
b)TC变换图像中存在明亮的特征,这些特征显然与雾度无关,而是与高可见波段反射地面特征相关,如建筑群、道路和裸露地面。另一方面,HOT对这些目标不太敏感。
c)覆盖了主要由成熟林分和采伐后再生区域组成的混交林山区。山峰通常是贫瘠的岩石,部分被冰雪覆盖。一层薄薄的层云覆盖着一个宽阔的山谷,占据了场景的中部。HOT详细地捕捉到了这一特征,
d)层云或卷云层覆盖了山谷景观,包括森林和更新的低植被斑块。
e)针叶林/落叶林和农田上方有两个轨迹和积云团
f)在图像右侧,一个广泛、薄、透明的卷云系统与明亮的白色积云团重叠。
g)有发育良好的厚积云团和薄冰云顶。薄薄的卷云层与这些对流云团重叠。
通过所提出的转换,可以详细地检测不同的大气条件。

HOT是基于这样一个假设,即TM1和TM3辐射将与场景最清晰部分内的像素高度相关,并且这种关系适用于所有地表覆盖类型。 当某些表面类型在不同程度上违反了这一假设,就可能触发虚假的、非大气的HOT。比如:

  • 积雪和雪上的阴影。这些区域很容易与多云像素混淆。 幸运的是,这个问题可以克服,因为雪和云可以通过其他方式分离,例如,通过基于TM波段5数据的分析生成雪遮罩
  • 水体。某些但并非所有水面上的HOT相对较高。 与雪一样,可以通过创建水体遮罩(例如,通过对近红外TM4图像进行阈值化)轻松克服此问题。
  • 裸露的土壤。 裸露土壤区域的HOT响应似乎相对较低。在非常低的雾霾条件下,与邻近作物相比,休耕地在热图像中是暗的。

5.HOT在可见光波段应用的调整

一种使用热图像对陆地卫星图像的可见波段进行辐射调整的方法。简言之,本程序包括以下步骤。

  • CL定义
  • HOT影像生成
  • 直方图生成
  • 辐射平差估计
  • 图像调整

5.1.CL定义

对场景进行目视检查,以确定其最清晰的区域。该区域中的示例子图像及其组成像素的TM1和TM3 DNs被回归以定义CL,从而定义HOT系数。

5.2.HOT影像生成

变换应用于每个像素以生成热图像(图7c的左面板)。主导场景中心部分的主要雾度特征显示的热值范围为30–40,而更高的值与周围的孤立云相关。

5.3.直方图生成

对于每个可见波段,将为根据其HOT级别分组的像素生成DN直方图。这一分组以从最清晰到最模糊的HOT级别的窄增量进行。
在这里插入图片描述在这里插入图片描述

根据热级别分组的像素的TM3 DN级别直方图。观察到直方图下限值随着热度的增加而增加,这与烟雾对图像辐射测量的加性效应是一致的。
由于在可见光带中,霾增加了观测到的辐射,预计随着HOT辐射水平的增加,这些直方图将系统地迁移到更高的DN水平。

5.4.辐射平差估计

利用一种“暗目标”减法将图像标准化为最清晰区域的辐射水平。以TM1为例。对于这个波段,清晰像素(即热=30)的直方图下限约为20 DNs。考虑一个模糊的像素,观察到的HOT水平为40。它是具有下限27的直方图的成员。这意味着在辐射测量调整阶段,像素的波段1 DN级别应减少7。

5.5.图像调整

说实话,没完全看懂

这篇关于遥感影像云检测-传统方法-:Haze Optimized Transformation(HOT)方法的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/661278

相关文章

Python实现图片分割的多种方法总结

《Python实现图片分割的多种方法总结》图片分割是图像处理中的一个重要任务,它的目标是将图像划分为多个区域或者对象,本文为大家整理了一些常用的分割方法,大家可以根据需求自行选择... 目录1. 基于传统图像处理的分割方法(1) 使用固定阈值分割图片(2) 自适应阈值分割(3) 使用图像边缘检测分割(4)

Java中Switch Case多个条件处理方法举例

《Java中SwitchCase多个条件处理方法举例》Java中switch语句用于根据变量值执行不同代码块,适用于多个条件的处理,:本文主要介绍Java中SwitchCase多个条件处理的相... 目录前言基本语法处理多个条件示例1:合并相同代码的多个case示例2:通过字符串合并多个case进阶用法使用

Python中__init__方法使用的深度解析

《Python中__init__方法使用的深度解析》在Python的面向对象编程(OOP)体系中,__init__方法如同建造房屋时的奠基仪式——它定义了对象诞生时的初始状态,下面我们就来深入了解下_... 目录一、__init__的基因图谱二、初始化过程的魔法时刻继承链中的初始化顺序self参数的奥秘默认

html5的响应式布局的方法示例详解

《html5的响应式布局的方法示例详解》:本文主要介绍了HTML5中使用媒体查询和Flexbox进行响应式布局的方法,简要介绍了CSSGrid布局的基础知识和如何实现自动换行的网格布局,详细内容请阅读本文,希望能对你有所帮助... 一 使用媒体查询响应式布局        使用的参数@media这是常用的

Spring 基于XML配置 bean管理 Bean-IOC的方法

《Spring基于XML配置bean管理Bean-IOC的方法》:本文主要介绍Spring基于XML配置bean管理Bean-IOC的方法,本文给大家介绍的非常详细,对大家的学习或工作具有一... 目录一. spring学习的核心内容二. 基于 XML 配置 bean1. 通过类型来获取 bean2. 通过

基于Python实现读取嵌套压缩包下文件的方法

《基于Python实现读取嵌套压缩包下文件的方法》工作中遇到的问题,需要用Python实现嵌套压缩包下文件读取,本文给大家介绍了详细的解决方法,并有相关的代码示例供大家参考,需要的朋友可以参考下... 目录思路完整代码代码优化思路打开外层zip压缩包并遍历文件:使用with zipfile.ZipFil

Python处理函数调用超时的四种方法

《Python处理函数调用超时的四种方法》在实际开发过程中,我们可能会遇到一些场景,需要对函数的执行时间进行限制,例如,当一个函数执行时间过长时,可能会导致程序卡顿、资源占用过高,因此,在某些情况下,... 目录前言func-timeout1. 安装 func-timeout2. 基本用法自定义进程subp

Python列表去重的4种核心方法与实战指南详解

《Python列表去重的4种核心方法与实战指南详解》在Python开发中,处理列表数据时经常需要去除重复元素,本文将详细介绍4种最实用的列表去重方法,有需要的小伙伴可以根据自己的需要进行选择... 目录方法1:集合(set)去重法(最快速)方法2:顺序遍历法(保持顺序)方法3:副本删除法(原地修改)方法4:

Python中判断对象是否为空的方法

《Python中判断对象是否为空的方法》在Python开发中,判断对象是否为“空”是高频操作,但看似简单的需求却暗藏玄机,从None到空容器,从零值到自定义对象的“假值”状态,不同场景下的“空”需要精... 目录一、python中的“空”值体系二、精准判定方法对比三、常见误区解析四、进阶处理技巧五、性能优化

C++中初始化二维数组的几种常见方法

《C++中初始化二维数组的几种常见方法》本文详细介绍了在C++中初始化二维数组的不同方式,包括静态初始化、循环、全部为零、部分初始化、std::array和std::vector,以及std::vec... 目录1. 静态初始化2. 使用循环初始化3. 全部初始化为零4. 部分初始化5. 使用 std::a