遥感影像云检测-传统方法-: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

相关文章

SpringBoot日志配置SLF4J和Logback的方法实现

《SpringBoot日志配置SLF4J和Logback的方法实现》日志记录是不可或缺的一部分,本文主要介绍了SpringBoot日志配置SLF4J和Logback的方法实现,文中通过示例代码介绍的非... 目录一、前言二、案例一:初识日志三、案例二:使用Lombok输出日志四、案例三:配置Logback一

Python实现无痛修改第三方库源码的方法详解

《Python实现无痛修改第三方库源码的方法详解》很多时候,我们下载的第三方库是不会有需求不满足的情况,但也有极少的情况,第三方库没有兼顾到需求,本文将介绍几个修改源码的操作,大家可以根据需求进行选择... 目录需求不符合模拟示例 1. 修改源文件2. 继承修改3. 猴子补丁4. 追踪局部变量需求不符合很

mysql出现ERROR 2003 (HY000): Can‘t connect to MySQL server on ‘localhost‘ (10061)的解决方法

《mysql出现ERROR2003(HY000):Can‘tconnecttoMySQLserveron‘localhost‘(10061)的解决方法》本文主要介绍了mysql出现... 目录前言:第一步:第二步:第三步:总结:前言:当你想通过命令窗口想打开mysql时候发现提http://www.cpp

Mysql删除几亿条数据表中的部分数据的方法实现

《Mysql删除几亿条数据表中的部分数据的方法实现》在MySQL中删除一个大表中的数据时,需要特别注意操作的性能和对系统的影响,本文主要介绍了Mysql删除几亿条数据表中的部分数据的方法实现,具有一定... 目录1、需求2、方案1. 使用 DELETE 语句分批删除2. 使用 INPLACE ALTER T

MySQL INSERT语句实现当记录不存在时插入的几种方法

《MySQLINSERT语句实现当记录不存在时插入的几种方法》MySQL的INSERT语句是用于向数据库表中插入新记录的关键命令,下面:本文主要介绍MySQLINSERT语句实现当记录不存在时... 目录使用 INSERT IGNORE使用 ON DUPLICATE KEY UPDATE使用 REPLACE

CentOS 7部署主域名服务器 DNS的方法

《CentOS7部署主域名服务器DNS的方法》文章详细介绍了在CentOS7上部署主域名服务器DNS的步骤,包括安装BIND服务、配置DNS服务、添加域名区域、创建区域文件、配置反向解析、检查配置... 目录1. 安装 BIND 服务和工具2.  配置 BIND 服务3 . 添加你的域名区域配置4.创建区域

mss32.dll文件丢失怎么办? 电脑提示mss32.dll丢失的多种修复方法

《mss32.dll文件丢失怎么办?电脑提示mss32.dll丢失的多种修复方法》最近,很多电脑用户可能遇到了mss32.dll文件丢失的问题,导致一些应用程序无法正常启动,那么,如何修复这个问题呢... 在电脑常年累月的使用过程中,偶尔会遇到一些问题令人头疼。像是某个程序尝试运行时,系统突然弹出一个错误提

电脑提示找不到openal32.dll文件怎么办? openal32.dll丢失完美修复方法

《电脑提示找不到openal32.dll文件怎么办?openal32.dll丢失完美修复方法》openal32.dll是一种重要的系统文件,当它丢失时,会给我们的电脑带来很大的困扰,很多人都曾经遇到... 在使用电脑过程中,我们常常会遇到一些.dll文件丢失的问题,而openal32.dll的丢失是其中比较

python中字符串拼接的几种方法及优缺点对比详解

《python中字符串拼接的几种方法及优缺点对比详解》在Python中,字符串拼接是常见的操作,Python提供了多种方法来拼接字符串,每种方法有其优缺点和适用场景,以下是几种常见的字符串拼接方法,需... 目录1. 使用 + 运算符示例:优缺点:2. 使用&nbsjsp;join() 方法示例:优缺点:3

Mysql中深分页的五种常用方法整理

《Mysql中深分页的五种常用方法整理》在数据量非常大的情况下,深分页查询则变得很常见,这篇文章为大家整理了5个常用的方法,文中的示例代码讲解详细,大家可以根据自己的需求进行选择... 目录方案一:延迟关联 (Deferred Join)方案二:有序唯一键分页 (Cursor-based Paginatio