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

#### 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 bandmath的地表温度反演

• esricd
• 2012年05月04日 10:51
• 5181

## ETM蚀变信息提取3-- Flaash大气纠正

• kupe87826
• 2012年10月30日 21:39
• 5039

• cmoster
• 2016年06月17日 11:14
• 3056

• lucky51222
• 2014年03月30日 21:38
• 3612

## Security solution for TO C TO B

Security solution for TO C TO B
• cnbird2008
• 2016年12月14日 21:11
• 663

## Check Soduku Solution

Check Soduku Solution: http://bo-yang.github.io/2014/06/14/soduku
• bonny95
• 2014年09月13日 12:21
• 458

## PND Solution(Design House) (China)

PND Solution(Design House) (China)
• king_sundi
• 2011年04月13日 13:55
• 2074

## 海洋内波简介

• xuanzi_eli
• 2016年08月27日 15:19
• 1592

## HDOJ 2719 The Seven Percent Solution

The Seven Percent Solution Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (J...
• A_ice_
• 2016年07月18日 17:24
• 109

## 富春江水污染事件

• erkeguoji8
• 2014年05月19日 11:30
• 318

举报原因： 您举报文章：#海洋以及湖泊离水辐亮度反演研究 色情 政治 抄袭 广告 招聘 骂人 其他 (最多只允许输入30个字)