# #海洋以及湖泊离水辐亮度反演研究

#### 1.3波段运算，得出反演结果

##### 1.3.1 反演大气瑞利散射辐亮度。

IDL波段运算表达式：
Fs*(0.008569*(λ^(-4))(1+0.0113(λ^(-2))+0.00013*λ^(-4))(P/1013.25))((0.75*(1+(-cos (((B2)))cos((B3))-sin((B2))*sin((B3))*cos((B4)-(B5)))^2))+((0.5((sin((B2)-asin(sin ((B2)/1.34)))/sin ((B2)+asin (sin ((B2)/1.34))))^2+(tan ((B2)-asin (sin ((B2)/1.34)))/tan ((B2)+asin (sin ((B2)/1.34))))^2))+(0.5*((sin ((B3)-asin (sin ((B3)/1.34)))/sin ((B3)+asin (sin ((B3)/1.34))))^2+(tan ((B3)-asin (sin ((B3)/1.34)))/tan ((B3)+asin (sin ((B3)/1.34))))^2)))(0.75(1+(cos ((B2)) cos ((B3))-sin ((B2))*sin ((B3))*cos((B4)-(B5)))^2)))(1/(cos ((B2))+0.15*(93.885-(B2))^(-1.253) *(P/1013.25))) /(4*3.1415926)

Figure 1瑞利散射辐亮度

##### 1.3.2. 反演气溶胶辐射亮度值。
Gordon等人认为气溶胶的水平分布通常比较均匀，假设气溶胶类型、气溶胶光学厚度在小范围内(（50～100 km）基本不变 。气溶胶散射辐射亮度值为:


Fs*0.9*0.896*(((0.983*(1-0.82*0.82))/(1+0.82*0.82-0.82*0.82*(cos ((B2)) cos ((B3))-sin((B2))*sin((B3))*cos((B4)-(B5))))^1.5+((1-0.983)(1-(-0.55)^2))/(1+(-0.55)^2-2*(0.55)(cos ((B2)) *cos ((B3))-sin ((B2))*sin ((B3))*cos((B4)-(B5))))^1.5)+(0.5((sin ((B2)-asin (sin ((B2)/1.34)))/sin ((B2)+asin (sin ((B2)/1.34))))^2+(tan ((B2)-asin (sin ((B2)/1.34)))/tan ((B2)+asin (sin ((B2)/1.34))))^2)+0.5*((sin ((B3)-asin (sin ((B3)/1.34)))/sin ((B3)+asin (sin ((B3)/1.34))))^2+(tan ((B3)-asin (sin ((B3)/1.34)))/tan ((B3)+asin (sin ((B3)/1.34))))^2))((0.983(1-0.82*0.82))/(1+0.82*0.82-0.82*0.82*(-cos ((B2)) cos ((B3))-sin ((B2))*sin((B3))*cos((B4)-(B5))))^1.5+((1-0.983)(1-(-0.55)^2))/(1+(-0.55)^2-2*(0.55)(-cos ((B2)) *cos ((B3))-sin ((B2))*sin ((B3))*cos((B4)-(B5))))^1.5)) (1/(cos ((B2))+0.15*(93.885-(B2))^(-1.253) *(P/1013.25))) /(4*3.1415926)

Figure 2气溶胶辐射亮度

##### 1.3.3. 反演臭氧的光学厚度 和大气漫透射率t(λ,θ)
臭氧的光学厚度 求法如下：



1 620-670 645 163.00 0.07 250m
2 841 858 100.58 0 250m
8 405~420 412 180.50 0 1000m
9 438~448 443 194.50 0.003 1000m
10 483~498 488 187.50 0.019 1000m
11 526~536 531 196.00 0.064 1000m
12 546~556 551 186.50 0.085 1000m
13 662~672 667 154.00 0.049 1000m
14 673~683 678 150.50 0.04 1000m
15 743~753 748 128.05 0.0092 1000m
16 862~877 865 99.05 0 1000m
DU=346 =0.346，吸收系数如上表，臭氧光学厚度 =0.02422。

IDL波段运算表达式：

(exp(-((0.008569*λ^(-4)(1+0.0113λ^(-2)+0.00013*λ^(-4))*(P/1013.25))/2+0.896+0.02422)/cos (B2)))

Figure 3大气漫透射率

##### 1.3.4. 反演离水辐亮度
水色遥感大气校正方程：



IDL波段运算表达式：
(B1-(Fs*(0.008569*(λ^(-4))(1+0.0113(λ^(-2))+0.00013*λ^(-4))(P/1013.25))((0.75*(1+(-cos (((B2)))cos((B3))-sin((B2))*sin((B3))*cos((B4)-(B5)))^2))+((0.5((sin((B2)-asin(sin ((B2)/1.34)))/sin ((B2)+asin (sin ((B2)/1.34))))^2+(tan ((B2)-asin (sin ((B2)/1.34)))/tan ((B2)+asin (sin ((B2)/1.34))))^2))+(0.5*((sin ((B3)-asin (sin ((B3)/1.34)))/sin ((B3)+asin (sin ((B3)/1.34))))^2+(tan ((B3)-asin (sin ((B3)/1.34)))/tan ((B3)+asin (sin ((B3)/1.34))))^2)))(0.75(1+(cos ((B2)) cos ((B3))-sin ((B2))*sin ((B3))*cos((B4)-(B5)))^2)))(1/(cos ((B2))+0.15*(93.885-(B2))^(-1.253) (P/1013.25))) /(4*3.1415926))-( Fs*0.9*0.896(((0.983*(1-0.82*0.82))/(1+0.82*0.82-0.82*0.82*(cos ((B2)) cos ((B3))-sin((B2))*sin((B3))*cos((B4)-(B5))))^1.5+((1-0.983)(1-(-0.55)^2))/(1+(-0.55)^2-2*(0.55)(cos ((B2)) *cos ((B3))-sin ((B2))*sin ((B3))*cos((B4)-(B5))))^1.5)+(0.5((sin ((B2)-asin (sin ((B2)/1.34)))/sin ((B2)+asin (sin ((B2)/1.34))))^2+(tan ((B2)-asin (sin ((B2)/1.34)))/tan ((B2)+asin (sin ((B2)/1.34))))^2)+0.5*((sin ((B3)-asin (sin ((B3)/1.34)))/sin ((B3)+asin (sin ((B3)/1.34))))^2+(tan ((B3)-asin (sin ((B3)/1.34)))/tan ((B3)+asin (sin ((B3)/1.34))))^2))((0.983(1-0.82*0.82))/(1+0.82*0.82-0.82*0.82*(-cos ((B2)) cos ((B3))-sin ((B2))*sin((B3))*cos((B4)-(B5))))^1.5+((1-0.983)(1-(-0.55)^2))/(1+(-0.55)^2-2*(0.55)(-cos ((B2)) *cos ((B3))-sin ((B2))*sin ((B3))*cos((B4)-(B5))))^1.5)) (1/(cos ((B2))+0.15*(93.885-(B2))^(-1.253) (P/1013.25))) /(4*3.1415926)))/ (exp(-((0.008569λ^(-4)(1+0.0113λ^(-2)+0.00013*λ^(-4))*(P/1013.25))/2+0.896+0.02422)/cos (B2)))

在实际的反演过程中遇到了很多难题，总的来说，有数据获取上的难题，反演过程上的难题，还有实际操作过程中的难题


#### 2.1 数据获取中存在的问题

##### 2.1.1 臭氧数据

ppt里提到的臭氧的一些资料现在已经无法使用，不能获得准确的臭氧单位。

NASA以前OMI数据已经不再提供，被新的臭氧传感器OMPS所代替。但是历史数据仍然能被检索到。

Figure 4臭氧数据下载

##### 2.1.2 气溶胶数据产品。
暂时用MODIS气溶胶04级的气溶胶光学厚度产品。但是对于MODIS04级产品不熟悉，没有了解到04级各波段所代表的产品是什么。



#### 2.2 波段运算中存在的问题

最开始几次反演离水辐亮度得到的图都是无效值NaN，和老师探讨多次，终于找到一些解决的办法。


(1)MODIS中的角度都是以特定格式存储的，在波段运算时要乘比例系数scale=0.01，然后转换成弧度制。一开始对波段运算了解的不多，计算用的是角度制，得到的概率因子都是无效值NaN。

Figure 5 MODIS太阳高度角的存储格式

(2)波段运算时，系数尽量以小数形式表示。

(3)反演离水辐亮度波段运算公式较长，不能一次性将所有的公式全部输入。尽量分步走，将公式拆分成几个部分分块处理。这样即可以得到各部分的结果图，又便于检查波段运算过程中的错误。

3离水辐射亮度反演结果分析

5
Figure 6油膜影响区域

Figure 7实际反演的离水辐亮度图

Figure 8 离水辐亮度图中疑似油膜放大图

Figure 10疑似油膜b点 离水辐亮度172.77

Figure 11渤海内海水 离水辐亮度131.75

Figure 12渤海口海水 离水辐亮度95.31

Figure 13太平洋海水 离水辐亮度55.21

１通过图像目视判读，

２通过离水辐亮度值。

### 工作展望：

1.进一步完善细节.气溶胶部分的反演过程较为粗糙，只是简单的利用MODIS气溶胶产品，没有考虑水汽对大气的影响，下一步要展开更深层次的气溶胶散射机制研究。

2.考虑利用交叉验证的方法来验证该反演方法的准确程度。

#### ENVI/IDL实现HJ卫星气溶胶反演

2013-07-19 09:55:15

#### 猴子都能看懂的PBR（才怪）

2018-02-03 15:36:08

#### 游戏中基于物理的渲染简简明科普

2016-02-02 10:39:43

#### 物理渲染(PBR)-基于物理的光照模型

2015-03-28 16:59:24

#### 定量遥感和程辐射

2011-09-02 19:55:52

#### 物理渲染-基于物理的光照模型

2015-10-11 08:43:03

#### 基于物理着色

2016-09-01 15:46:22

#### PBR——基于物理的实时渲染

2017-02-18 12:41:00

#### 输入日期，计算太阳表观辐亮度

2011年03月09日 2KB 下载

#### MODIS热红外数据辐射亮度值的计算

2015-10-18 15:24:45