Page 51 - 2023年第54卷第8期
P. 51
(1)基于 PRF - ANN(Parallelrandom forestandartificialneuralnetworkalgorithm)河流表面信息提取
方法 [26] 及 ARWE河宽自动提取算法,提取流域高分辨率河流表面信 息、河 流中 心线 及平 滩、最枯
河宽。
( 2)基于 Birkett算法 [38] 提取经河流表面信息去噪的研究区域 Jason - 1?2卫星的平滩及枯水水位。
( 3)确定 Jason卫星轨道与步骤(1)提取的河流表面信息的交线,交线的端点坐标为 A、B且认为
A、B处水位相等。读取该处水位 h(图 3)。
( 4)按水位递减顺序,从左岸至右岸依次连接端点得图 3中最低水位以上断面。
图 3 多源遥感观测的河流断面示意
4.2.2 河宽自动提取 基于 PRF - ANN及 ARWE法提取河宽的流程为:(1)基于 PRF - ANN法提取平
滩、最枯流量下的河流表面信息。(2)基于 ARWE算法,提取对应河流表面信息的河流中心线及水面
边界并校正。(3)ARWE基于校正后的河流中心线及水面边界,自动提取河宽。
4.2.3 平滩河宽及枯水河宽提取结果 基于 PRF - ANN法,通过预实验确定的流域丰、枯水期能较好
区分水体及其他地物的指数,依据这些指数提取流域内水体表面信息(包含河流、湖泊、水库等水体
3 +
的初步产品)。如丰水期(6—9月)选用了三价铁离子指数(Fe )、归一化水体指数(NDWI)、增强水
体指数( EVI)、色彩渲染指数 550(CRI550)4种指数、2种 Landsat - 8TIRS热红外数据和 4种灰度共生
矩阵(GLCM)纹理指数(绿光、近红外波段均值,红光、近红外波段熵)提取平滩流量下的水体表面信
3 +
息。枯水期使用形状指数(IF)、氧化铁指数(IO)、显色指数(CI)、Fe 及蓝光、绿光、红光波段的 4
种纹理指数提取最枯流量下的水体表面信息。
水体表面信息提取完成后,需基于由流域 DNET河流水系生产的缓冲区(半径 300~1000m)去除
湖泊、水库、水塘等非河流水体和其他地物,再通过形态学算法、适当人工编辑后,得到精确的河流
表面信息 [20] 。基于精确的河流表面信息,ARWE算法结合 DYRIM的子流域划分结果共提取了 1360个
平滩、枯水河宽数据对(图 4(a))。提取的平滩河宽通过 80组野外实测数据及 15组经高分辨率 DEM
断面算得的平滩河宽检验,枯水期最枯水位河宽则通过水文站及 GoogleEarth历史影像检验。平滩河
2
宽提取的 R、RMSE、平均偏差 MBE分别为 0.99、1.49、0.06,最小河宽对应 DNET的 2级河流;枯
2
水河宽提取的 R、RMSE、MBE分别为 0.91、1.61、0.08,最小河宽对应 DNET的 3级河流。提取所得
流域枯、平滩河宽之比为 0.31,实测值为 0.25,整体结果较准,提取结果的误差分析示意见图 4(b)。
平滩、枯水河宽数据对将用以提取相应水位。
4.2.4 平滩及枯水水位的自动提取 提取 Jason系列卫星 103号轨道与皇甫川流域河流交线,据交线
与河流边界交点确定水位提取点,假定每时期、每断面左右交点水位相等,概化河流断面。受限于
Jason卫星的测量幅宽和过境时间等,在 2012年 8月 6日、2017年 4月 17日仅从皇甫川流域自动提取
到 3个断面,如图 4(a),断面关键参数见表 2。
— 9 3 5 —