反演离水辐亮度流程
离水辐亮度分成三部分进行求解,一是减去大气瑞利散射,然后减去气溶胶的影响,最后除以漫透射率得到离水辐亮度。流程图如下:
目前,我所做的工作主要有如下几块:
1.1遥感影像以及各类参数数据的搜集
前期的数据准备需要大量的光学遥感影像,所以选用了高时间分辨率的MODIS影像作为实验数据。我通过NASA网站,获取到合适时间尺度内(152天到175天)的MODIS光学影像数据。
反演过程中还需要用到气溶胶光学厚度AOD以及臭氧光学厚度,对于气溶胶光学厚度,我们可以用Modis的04级气溶胶产品代替。而臭氧光学厚度需要我们获取实验区域内的臭氧单位。
1.2数据预处理
数据预处理最主要的工作就是对影像几何校正,然后把各类数据重采样为同一分辨率,方便接下来的波段运算工作
1.3波段运算,得出反演结果
波段运算是将离水辐亮度的反演公式从数学表达式转换成idl可读表达式的过程,这也是反演过程中最主要的工作量之一。
波段运算我将其拆成三小部分完成。分别是反演大气瑞利散射辐亮度,气溶胶辐亮度,臭氧光学厚度,以及通过大气漫透过率反演离水辐亮度。
考虑到太阳耀斑相对于瑞丽散射以及气溶胶对离水辐亮度的影响较小,暂时没有进行太阳耀斑的提取工作。
1.3.1 反演大气瑞利散射辐亮度。
水色遥感中,瑞利散射占传感器接收的总辐射的较大份额,特别时在蓝光波段占80%-90%。太阳辐射经过大气层后,有一部分经大气分子后向散射到达卫星传感器,还有一部分到达海面,经海面漫反射后到达传感器,卫星接收到的大气瑞利散射为:
是太阳以三种不同的通道散射到卫星的概率因子, 是单次散射反照率(非水汽和氧气吸收带的波段近似为1),瑞利光学厚度 ,Fs是经太阳-地球距离校正后的辐照度。
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