在分析由风暴潮引起的洪水淹没和侵蚀时,水深数据是很重要的。较高分辨率的浅海水深数据能提高在风暴潮危害评估中的建模能力。同时,水深数据在航海图更新、探索海底环境、沉积物迁移研究中也很重要,而阿拉斯加威尔士区域的水深数据非常少。本项目是通过遥感方法得到水深数据,利用实测声纳数据,使用线性回归的方法对结果进行验证,R2 = 0.7221,基本可以满足应用的精度需要。
图1:威尔士区位置及航空拍摄图
图2:2011年9月份拍摄的海岸侵蚀图像
处理流程和方法
遥感水深测量是利用可见光遥感数据反演水深,主要依据为:在可见光范围内,波长在0.4~0.58μm 之间的光信号对清澈水体的穿透深度最大,大气条件较好时,能够探测水深在30米以内的水体;(2)反映水深的辐射量与光信号在水中的衰减程度呈指数关系。
图3:太阳光水面辐射示意图
WorldView2多光谱数据包括除了一般的蓝、绿、红、近红外波段,还包括海岸蓝色波段(400-450nm)、黄色波段(585-625nm)。Lee et al.(2012), Madden( 2011), Miecznik et al.(2012), Tarantino et al.( 2012)等人已经成功在热带地区利用WorldView2测量水深。如果这种方法适合于阿拉斯加的寒带地区,将产生非常大的价值。使用处理流程如下图所示。
图4:处理流程图
(一) 底部反照率独立水深测量算法
使用底部反照率独立水深测量算法(bottom albedo-independent Bathymetry algorithm)量测水深,这个算法假设:当深度相同的时候,不管水底是被深色水草或者明亮的沙子覆盖,他们都显示为同一深度。计算公式如下:
(公式1)
其中:
D=水深
m1=常量,用于调整水深比例
n=固定的常量,保证对数是正值及是线性关系
RW=卫星观测辐射亮度值
λi=波段i
λj=波段j
m0=常量,用于补偿0m深度
(二) SPEAR Relative Water Depth Wizard
SPEAR Relative Water Depth Wizard(RWDW)是ENVI下的一个流程化工具,集成了底部反照率独立水深测量算法。在RWDW工具中分别浏览RGB532、531、432、431,最后得到RGB431分辨力最好和最少视觉噪声(如图5a所示)。RWDW根据公式1直接得到相对水深(如图5b所示)。
图5 a, b, c:相对水深制图结果(b) RGB组合对比(a) NOAA航海图 (2004)(c)
对比NOAA的航海图(图5c),标定对应区域的真实水深,如图6所示。为了降低噪声,结果图像从2米转化到16米空间分辨率, 剪出近海岸小于3m的区域。
图6:标定得到真实水深
(三) 结果
采用线性回归模型对比水深结果和声纳数据,如图8所示。使用N=1480组数据线性回归的R2=0.7221。水深结果和声纳量测数据的绝对差值如图9所示,图9显示随着深度的减小相关性增加。从反演和测量数据点(图11a-d)的横截面看,具有光谱高反射率的位置处于海槽位置(通过实测数据)。如果不存在高反射区,反演值和测量值非常相似,直到在3米这个地方开始分散(图11d)。
图7:声纳设备和实测路径
图8:测量和置信度95%之间数据的线性回归
图9:水深结果和声纳量测数据的绝对差值
图10 a, b, c & d:反演与测量值的横切面对比图,Datum: NAD 83 Projection: UTM Zone 3 North Vertical Datum: NAVD 8
结论
- 在阿拉斯加威尔士区域,WorldView2卫星影像可以用于反演小于3m的近海岸水深。
- WorldView2的蓝色波段、黄色波段可以增强光谱可分离度。
- 泥沙可能会对近海岸反射率产生影响。
《WorldView-2数据》共享方法如下:
(1)人员,限定为关注小编的用户。
(2)各类项目(包括各类科研项目)申请本数据扔享受免费政策,但需向本号捐赠一定数量的硬盘才能获取。
(3)捐赠硬盘可免留言获取数据。