遥感影像大气校正一:6S模型使用

6S介绍        

        1986年,法国里尔科技大学大气光学实验室Tanré等人为了简化大气辐射传输方程,开发了太阳光谱波段卫星信号模拟程序5S(SIMULATION OF THE SATELLITE SIGNAL IN THE SOLAR SPECTRUM),用来模拟地气系统中太阳辐射的传输过程并计算卫星入瞳处辐射亮度。

        1997年,Eric Vemote对5S进行了改进,发展到6S(SECOND SIMULATION OF THE SATELLITE SIGNAL IN THE SOLAR SPECTRUM),6S吸收了最新的散射计算方法,使太阳光谱波段的散射计算精度比5S有所提高。

        本文以GF1-PMS第一波段为例介绍6s模型使用(使用版本6Sv2.1)。

6S输入参数介绍

第一个参数:卫星的几何条件(Geometrical conditions):

        0-7中的数字

        0.use definition(用户自定义);1.meteosat observation ; 2.goes east observation ; 3.goes west ;4.avhrr ( PM noaa ) ; 5.avhrr ( AM noaa ) ;  6.hrv( spot ) ; 7.tm( landsat )

        输入0的话下一个参数要指定观测几何条件,其他1-7为内置的已经配置好的观测条件。

        这里选择0!

第二个参数:输入时间和位置(根据第一个参数确定输入内容):

几何条件

类型

输入对应时间,角度

0

输入:qs(太阳天顶角), fs(太阳方位角), qv(卫星天顶角), fv(卫星方位角), Month(月), Day(日)

1

meteosat observation

(气象卫星观测)

输入: Month, Day, Hour, Nb. cols, Nb. rows

month,day,decimal hour(universal time-hh.ddd) 

n. of column,n. of line.(full scale 5000*2500) 

2

goes east observation(GOES(东)卫星

输入: Month, Day, Hour, Nb. cols, Nb. rows

month,day,decimal hour (universal time-hh.ddd) 
n. of column,n. of line.(full scale 17000*12000)c

3

goes east observation(GOES(西)卫星))

输入: Month, Day, Hour, Nb. cols, Nb. rows

month,day,decimal hour (universal time-hh.ddd) 

n. of column,n. of line.(full scale 17000*12000)

4

avhrr ( PM noaa )

下午NOAA卫星

输入: Month, Day, Hour, Nb. cols, Nb. rows

month,day,decimal hour (universal time-hh.ddd) 
n. of column(1-2048),xlonan,hna

give long.(xlonan) and overpass hour (hna) at the ascendant node at equator 

5

avhrr ( AM noaa )

(上午NOAA卫星

输入:Month, Day, Hour, Nb. cols, Long A.N., Hour A.N.

month,day,decimal hour (universal time-hh.ddd) 
n. of column(1-2048),xlonan,hna 
give long.(xlonan) and overpass hour (hna) at the ascendant node at equator 

6

hrv   ( spot )

输入: Month, Day, Hour, Long, Lat

7

tm    ( landsat )

输入: Month, Day, Hour, Long, Lat

  这里输入:60 150 2 0 11 27

第三个参数:输入大气的模型( Atmospheric model ):

模型代号

含义

0

无气体的吸收(no gaseous absorption)

1

热带的(Tropical)

2

中纬度夏季(midlatitude summer)

3

中纬度冬季(midlatitude winter )

4

亚北极夏季(subarctic summer)

5

亚北极(subarctic winter)

6

北美标准(us standard 62 )

7

用户定义大气廓线(34层无线电探空数据)

包括:高度(km 气压( mb )  温度( k )  水汽密度( g/m3)  臭氧密度(g/m3)

输入: 共34行,每行为: altitude  pressure   temperature  h2o  o3 

单位:   altitude ( in km )  pressure(in mb )  temperature ( in k ) h2o density(in g/m3)  o3 density    (in g/m3) 

for example, altitudes are  from  0 to 25km step of 1km 
                     from 25 to 50km step of 5km 
                     and two values at 70km and 100km 
                     so you have 34*5 values to input.

8

用户自定义 输入水汽和臭氧总含量
输入:    uw  (in  g/cm2 )   uo3 (in  cm-atm) 

water vapor and ozone contents

9读入无线电探空数据文件

这里输入 3

第四个参数:输入气溶胶模型(Aerosols model):

模型代号

含义

0

无气溶胶(no aerosols )

1

大陆模型(continental model)

2

海洋模型(maritime model)

3

城市模型(urban model)

4

用户模型(using basic components)

下一行输入以下四种粒子所占体积百分比(0-1)

c(1) : 灰尘    

c(2) :水溶型

c(3) :海洋型      

c(4) :烟灰

如:0.25 0.25 0.25 0.25

5

沙漠模型(Desertic Model)

6

生物圈模型(biomass Model)

7

同温层模型(stratospheric model)

8

用户模型(size distribution)

多峰对数正态分布

9

用户模型(size distribution)

改进的gamma分布

10

用户模型(size distribution)

Junge幂指数律分布

11

用户模型(size distribution)

按太阳光度计测量结果定义气溶胶模型

需要输入参数有:粒子半径(µm) 粒径分布(d V / d (logr),cm3/cm2/micron)和复折射指数的实部和虚部谱

12

用户模型(size distribution)

利用事先计算的结果给出文件名

这里输入 1

第五个参数:气溶胶含量参数(concentration)

        在集合{-1,0,>0的正值}中选一个值:

        如果输入-1代表没有气溶胶;

        如果输入>0的正数,代表气溶胶能见度(公里,必须大于5);

        如果输入0需要在下一行输入550纳米气溶胶光学厚度。

这里输入:0

                   0.145

第六个参数:目标高度参数(altitude of target)

        xps >=0:目标在海平面高度

        xps < 0:绝对值代表目标高度(公里)

这里输入 0

第七个参数:传感器高度参数(sensor altitude)

        xpp= -1000:卫星观测

        xpp=     0:地面观测

        -100< xpp <0:飞机观测,绝对值代表飞机相对于目标的高度(公里)

        对于飞机观测,必须输入飞机和地面之间的水汽,臭氧含量和550纳米气溶胶光学厚度,如无数据则输入负值,水汽和臭氧根据62年美国标准大气内差,气溶胶则根据2公里指数廓线计算,如:

        -3.3(飞机观测)

        -1.5 -0.35(水汽和臭氧含量)

        0.25(气溶胶光学厚度)

这里我们输入 -1000

第八个参数:光谱参数(spectral conditions)

可以从以下数字中选择:

模型代号

含义

-2

用户自定义

-1

用户自定义

0

用户自定义

1

用户自定义

-2:用户输入光谱范围的下限和上限(微米),滤光片函数为1,输出文件中给出单色结果。

-1:单色计算,用户给出单色波长(微米)

0:用户输入光谱范围的下限和上限(微米),滤光片函数为1

1:用户输入光谱范围的下限和上限(微米)并输入滤光片函数(间隔为0.0025微米。

注意:因为国产卫星的光谱参数没有被内置在代码里,所以这里需要我们自定义,一般是指定光谱范围然后下一行输入光谱响应值,间隔为0.0025微米,同时一定要注意,第二行输入的响应值个数(b-a)微米/4+1个,多一个或者少一个都不行!!!!!!

这里输入(3行):

        1
        0.450000 0.520000
     0.020823 0.061016 0.085883 0.203136 0.302934 0.315089 0.439565 0.531091 0.579737 0.675640 0.725676 0.747863 0.798039 0.835303 0.855022 0.883200 0.911129 0.936508 0.963268 0.986688 0.997939 0.997788 0.980309 0.916002 0.781131 0.593427 0.400828 0.250990 0.151027

        如果选择2-70,代表下列卫星通道:

2

  vis bandof meteosat     ( 0.350-1.110 ) 

3

  vis bandof goes east    ( 0.490-0.900 )

4

  vis bandof goes west    ( 0.490-0.900 )

5

  1st bandof avhrr(noaa6) ( 0.550-0.750 )

6

  2nd     "              ( 0.690-1.120 )

7

  1st bandof avhrr(noaa7) ( 0.500-0.800 )

8

  2nd     "              ( 0.640-1.170 )

9

  1st bandof avhrr(noaa8) ( 0.540-1.010 )

10

  2nd     "              ( 0.680-1.120 )

11

  1st band of avhrr(noaa9)( 0.530-0.810 )

12

  2nd     "              ( 0.680-1.170 )

13

  1st band of avhrr(noaa10) ( 0.530-0.780 )

14

  2nd     "              ( 0.600-1.190 )

15

  1st band of avhrr(noaa11)( 0.540-0.820 )

16

  2nd     "              ( 0.600-1.120 )

17

  1st band of hrv1(spot1) ( 0.470-0.650 )

18

  2nd     "              ( 0.600-0.720 )

19

  3rd     "              ( 0.730-0.930 )

20

  pan     "              ( 0.470-0.790 )

21

  1st band of hrv2(spot1) ( 0.470-0.650 )

22

  2nd     "              ( 0.590-0.730 )

23

  3rd     "              ( 0.740-0.940 )

24

  pan     "              ( 0.470-0.790 )

25

  1st band of tm(landsat5)( 0.430-0.560 )

26

  2nd     "              ( 0.500-0.650 )

27

  3rd     "              ( 0.580-0.740 )

28

  4th     "              ( 0.730-0.950 )

29

  5th     "              ( 1.5025-1.890 )

30

  7th     "              ( 1.950-2.410 )

31

  1st band of mss(landsat5)(0.475-0.640 )

32

  2nd     "              ( 0.580-0.750 )

33

  3rd     "              ( 0.655-0.855 )

34

  4th     "              ( 0.785-1.100 )

35

  1st band of MAS(ER2)    ( 0.5025-0.5875)

36

  2nd     "              ( 0.6075-0.7000)

37

  3rd     "              ( 0.8300-0.9125)

38

  4th     "              ( 0.9000-0.9975)

39

  5th     "              ( 1.8200-1.9575)

40

  6th     "              ( 2.0950-2.1925)

41

  7th     "              ( 3.5800-3.8700)

42

  MODIS  band 1           ( 0.6100-0.6850)

43

  MODIS  band 2           ( 0.8200-0.9025)

44

  MODIS  band 3           ( 0.4500-0.4825)

45

  MODIS  band 4           ( 0.5400-0.5700)

46

  MODIS  band 5           ( 1.2150-1.2700)

47

  MODIS  band 6           ( 1.6000-1.6650)

48

  MODIS  band 7           ( 2.0575-2.1825)

49

  1st band of avhrr(noaa12) ( 0.500-1.000 )

50

  2nd     "              ( 0.650-1.120 )

51

  1st band of avhrr(noaa14) ( 0.500-1.110 )

52

  2nd     "              ( 0.680-1.100 )

53

  POLDER  band1           ( 0.4125-0.4775)

54

  POLDER  band2 (non polar( 0.4100-0.5225)

55

  POLDER  band3 (non polar( 0.5325-0.5950)

56

  POLDER  band4   P1      ( 0.6300-0.7025)

57

  POLDER  band5 (non polar( 0.7450-0.7800)

58

  POLDER  band6 (non polar( 0.7000-0.8300)

59

  POLDER  band7   P1      ( 0.8100-0.9200)

60

  POLDER  band8 (non polar( 0.8650-0.9400)

61

ETM+1

62

ETM+2

63

ETM+3

64

ETM+4

65

ETM+5

66

ETM+7

67

ETM+8

这里其实也可以修改6s源码,直接将国产卫星光谱值写入到代码中,可参考基于6S模型的国产卫星数据大气校正_hj卫星6s大气校正-CSDN博客

第九个参数:输入地表是否均匀(HOMOGENEOUS CASE):

        0代表均匀,1代表不均匀

        当输入0:下一行输入是否直接影响(0代表无方向效应1代表有方向效应):

                当输入0:则无方向效应,则输入然后输入地面反射类型GROUND REFLECTANCE:   

        -1

用户定义

        0

 用户定义

        1

 VEGETA

        2

CLEARW

        3

SAND

        4

 LAKEW

                当输入-1:输入ro(暂不明确具体含义)

                当输入0:输入ro [组]:by step of 0.0025 micron      (暂不明确具体含义) 

        当输入1:则代表有方向效应,则输入BRDF模型代号:(暂不明确具体含义)

模型代号

含义

0

in-situ measurements

1

Hapke's model

2

Verstraete's model

3

Roujean's model

4

Walthall's model

5

Minnaert's Model

6

Ocean model

7

Iaquinta & Pinty's model

8

Rahman's model

9

Kuusk's model

这里输入:

0
0
1

第十个参数:激活大气订正方式:

        如果输入<-1的数,代表不激活大气校正

        如果输入>0,反演地面反射率,反演出的地面反射率满足大气层顶的辐射亮度=rapp(w/m2/str/mic)

        如果输入-1到0,反演地面反射率,反演出的地面反射率满足表观反射率;

这里输入:

1
-0.1(这个值对最后的系数没啥影响,在范围内即可)

第十一个参数:地面未极化:

        这个参数未找到具体含义,按照默认参数4输入。

输入参数汇总:

0
60 150 2 0 11 27
3
1
0
0.145
0
-1000
1
0.450000 0.520000
0.020823 0.061016 0.085883 0.203136 0.302934 0.315089 0.439565 0.531091 0.579737 0.675640 0.725676 0.747863 0.798039 0.835303 0.855022 0.883200 0.911129 0.936508 0.963268 0.986688 0.997939 0.997788 0.980309 0.916002 0.781131 0.593427 0.400828 0.250990 0.151027
0
0
1
1
-0.5
4

执行并输出

命令行调用:

输出:

        最后就是计算地表反射率的系数了。

参考文献:

6S大气校正模型说明 - CodeAntenna

太阳天顶角和方位角计算_太阳天顶角计算公式-CSDN博客

基于6S模型的国产卫星数据大气校正_hj卫星6s大气校正-CSDN博客

6S大气传输模型修改源码添加、自定义卫星光谱响应(以HJ-1BCCD为例) - 百度文库 (baidu.com)

定量遥感——大气校正.ppt-原创力文档 (book118.com)

大气辐射校正-6s-实习 - 豆丁网 (docin.com)

高分系列 - 知乎 (zhihu.com)

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/610729.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Hive之set参数大全-3

D 是否启用本地任务调试模式 hive.debug.localtask 是 Apache Hive 中的一个配置参数&#xff0c;用于控制是否启用本地任务调试模式。在调试模式下&#xff0c;Hive 将尝试在本地模式下运行一些任务&#xff0c;以便更容易调试和分析问题。 具体来说&#xff0c;当 hive.de…

spring cloud之集成sentinel

写在前面 源码 。 本文一起看下spring cloud的sentinel组件的使用。 1&#xff1a;准备 1.1&#xff1a;理论 对于一个系统来说&#xff0c;最重要的就是高可用&#xff0c;那么如何实现高可用呢&#xff1f;你可能会说&#xff0c;集群部署不就可以了&#xff0c;但事实并…

【C#】用CefSharp.ChromiumWebBrowser做winform开发过程中遇到的一些坑

文章目录 使用键盘F5刷新、F12打开控制台浏览器默认下载弹窗禁用GPU配置 使用键盘F5刷新、F12打开控制台 public class CEFKeyBoardHander : IKeyboardHandler {public bool OnKeyEvent(IWebBrowser browserControl, IBrowser browser, KeyType type, int windowsKeyCode, int…

【3D动画】Animcraft 基础运用

骨骼包地址&#xff1a;G:\Animcraft相关资料\motion_2024 入库&#xff1a;FBX文件 Mixamo 下载一个模型。格式&#xff1a;FBX 导入到animcraft软件中。

论文封面表格制作

原文参考&#xff1a;【【论文排版】论文封面完美对齐 强迫症重度患者的经验分享】https://www.bilibili.com/video/BV18f4y1p7hc?vd_source046491acdcff4b39fed20406b36a93e2 视频里up主介绍很详细。我自己也记录一下。 介绍一下如何完成论文封面信息的填写。 创建一个3列…

MATLAB对数据隔位抽取和插值的几种方法

对于串行的数据&#xff0c;有时我们需要转成多路并行的数据进行处理&#xff0c;抽取&#xff1b;或者是需要对数据进行隔点抽取&#xff0c;或对数据进行插值处理。此处以4倍抽取或插值为例&#xff0c;MATLAB代码实现。 文章目录 抽取方法一&#xff1a;downsample函数方法…

vue:如果.vue文件内容过多,我们可以这样拆分

一、标签内容太多&#xff0c;我们可以用组件的方式拆分 二、那如果JS的内容太多呢&#xff1f;因为耦合性太高&#xff0c;拆成组件后父子组件传值不方便&#xff0c;我们可以这样&#xff1a; 子组件: export default {data() {return {};},methods: {a(){alert(1)}} };父组…

算法练习:找出数组大于 n/2 的元素

题目&#xff1a; 给定一个大小为 n 的数组nums &#xff0c;返回其中的多数元素。多数元素是指在数组中出现次数 大于 n/2 的元素。你可以假设数组是非空的&#xff0c;并且给定的数组总是存在多数元素。 示例&#xff1a; 输入&#xff1a;nums [3,2,3] 输出&#xff1a;…

[足式机器人]Part2 Dr. CAN学习笔记 - Ch02动态系统建模与分析

本文仅供学习使用 本文参考&#xff1a; B站&#xff1a;DR_CAN Dr. CAN学习笔记 - Ch02动态系统建模与分析 1. 课程介绍2. 电路系统建模、基尔霍夫定律3. 流体系统建模4. 拉普拉斯变换&#xff08;Laplace&#xff09;传递函数、微分方程4.1 Laplace Transform 拉式变换4.2 收…

Exception sending a multicast message:Network is unreachable故障

出现这个故障就是没有连接到网络&#xff0c;如果虚拟机没有连接到本机&#xff0c;那么就会出现这个情况&#xff0c;当虚拟机连接到本机就会自动消失&#xff0c;同时如果是用电脑直接安装Ubuntu运行也会出现这个情况&#xff0c;应该是要连接到一个路由器里面这个情况才会消…

Python 借助装饰器,实现父对象使用子对象的方法!

文章目录 1. 引言2. 实践2.1 装饰器说明2.2 定义一个父类2.3 定义子类 3. 整体来看4. 参考 1. 引言 这个需求的背景是这样的&#xff0c;我们在写复杂代码时候&#xff0c;可能会用到 Python 的类继承的方法&#xff0c;即子对象继承父对象的一些属性方法。 在这个过程中&…

游戏后端如何实现服务器之间的负载均衡?

网络游戏已成为人们休闲娱乐的重要方式之一。而在游戏开发中&#xff0c;如何实现服务器之间的负载均衡是一个非常关键的问题。负载均衡不仅可以提高服务器的处理能力&#xff0c;还能保证游戏的稳定性和流畅性。本文将探讨游戏后端如何实现服务器之间的负载均衡。 一、负载均…

Eigen 中的传播计算

Eigen: Reductions, visitors and broadcasting

太实用了!关于ControlNet,这篇你一定要看

一篇文章教会你&#xff0c;从入门到使用。 这里是行者AI&#xff0c;我们专注于人工智能在游戏领域的研究和应用&#xff0c;凭借自研算法&#xff0c;推出游戏AI、智能内容审核、数据平台等产品服务。 controlNet控制生成 anypaint 看到这篇文章的你&#xff0c;一定也是AI绘…

SpringBoot整合人大金仓数据库KingBase

1 去KingBase官网下载驱动jar包 2 将解压得到的所有jar包放置在libs目录下&#xff08;没有就新建一个目录&#xff09; 3 在pom文件添加相关依赖 <!--添加KingBase所需要的依赖--> <dependency><groupId>com.kingbase</groupId><artifactId>kin…

DNS解析原理和k8s DNS 实践

1. 问题背景 1.1 域名解析异常 近期开发的一个功能&#xff0c;需要在k8s集群容器环境中调用公司内部api&#xff0c;api提供了内网域名&#xff0c;解析内网域名异常导致请求超时&#xff0c;因此梳理了下DNS的知识点。 可以先看到下面&#x1f447;这段配置&#xff0c;修…

向量数据库:Milvus

特性 Milvus由Go(63.4%),Python(17.0%),C(16.6%),Shell(1.3%)等语言开发开发&#xff0c;支持python&#xff0c;go&#xff0c;java接口(C,Rust,c#等语言还在开发中)&#xff0c;支持单机、集群部署&#xff0c;支持CPU、GPU运算。Milvus 中的所有搜索和查询操作都在内存中执行…

《每天十分钟》-红宝书第4版-基本引用类型

引用值&#xff08;或者对象&#xff09;是某个特定引用类型的实例。在 ECMAScript 中&#xff0c;引用类型是把数据和功能组织到一起的结构&#xff0c;经常被人错误地称作“类”。虽然从技术上JavaScript 是一门面向对象语言&#xff0c;但ECMAScript 缺少传统的面向对象编程…

Disruptor相关问题及答案(2024)

1、Disruptor 是什么&#xff0c;它解决了什么问题&#xff1f; Disruptor是一个高性能的内存中数据交换框架&#xff0c;由LMAX(伦敦多资产交易所)开发&#xff0c;目的是用于支持交易系统中极高的并发处理。它可以说是一个用于线程间消息传递的工具&#xff0c;但与传统的队…

第32次CCF计算机软件能力认证-第一题

仓库规划 问题描述 西西艾弗岛上共有 n个仓库&#xff0c;依次编号为 1…n。每个仓库均有一个m维向量的位置编码&#xff0c;用来表示仓库间的物流运转关系。 具体来说&#xff0c;每个仓库 i均可能有一个上级仓库 j&#xff0c;满足&#xff1a;仓库 j 位置编码的每一维均大…