[发明专利]多角度遥感水深决策融合反演方法

摘要

多角度遥感水深决策融合反演方法包括：第一步对多光谱遥感影像进行预处理，得到海表反射率；第二步现场实测水深值获取及处理；第三步单角度水深反演和水深段标识；第四步多角度水深反演融合；第五步水深反演精度验证；水深反演精度验证完成后，将最终水深值作为遥感图像实际水深值输出数据。与现有的反演方法相比，本发明利用多角度遥感对不同海况、不同底质和不同水质情形下多方面的成像能力，可为水深遥感反演带来丰富的数据信息，提高反演精度，尤其适用于复杂情况下浅水区域的海洋水深测量。

主权项

多角度遥感水深决策融合反演方法，其特征在于，包括以下步骤：第一步：对多光谱遥感影像进行预处理，得到海表反射率；所述预处理包括辐射亮度转换、大气校正和太阳耀斑去除；第二步：现场实测水深值获取及处理；获取实验区的水深数据和对应的经纬度坐标，通过潮汐表确认测量时刻的潮高值，将水深数据校正获得理论深度基准面的水深，再根据多光谱遥感影像的获取时刻，对理论深度基准面的水深数据进行瞬时水深的潮汐校正以获得瞬时水深；第三步：单角度水深反演和水深段标识；根据水深控制点处水深与对应影像像元反射率值间的关系，采用多波段模型进行统计回归出，输出该景影像水深反演的参数作为多角度水深反演融合的一项输入，并对多波段模型进行参数定标，多波段模型公式如下，$Z=A_0+\underset{i=1}{\overset{n}{\Σ}}{A}_i{X}_i---\left(1\right)$X_i＝Ln(ρ_i‑ρ_si)   (2)其中，Z为水深，n为参与反演的波段个数，A₀和A_i为待定系数，ρ_i是第i波段反射率数据，ρ_si是该波段深水处的反射率；将水深控制点分成多个水深段作为输入，输出各水深段的平均相对误差和可信度，作为多角度水深反演融合的另一项输入，即融合参数，${\mathrm{\δ}}_k=\frac{1}{n}\underset{i=1}{\overset{n}{\Σ}}\frac{|z_i-{z_i}^{\′}|}{z_i}---\left(3\right)$$L_k=0.5+{\left\{\underset{i=1}{\overset{n}{\Σ}}{w}_i{\left(L_k^i-0.5\right)}^{\mathrm{\α}}\right\}}^{\frac{1}{\mathrm{\α}}}---\left(4\right)$其中，k表示水深段，n为水深检查点个数，式3中，δ_k是第k个水深段的平均相对误差，z_i是第i个水深检查点的实测值，z_i'为其反演值，式4中，L_k表示第k个水深段的可信度，w_i表示每个信源的相对重要程度，且满足其公式为代表该i景影像控制点的Kappa系数，n表示参与反演融合影像的个数，α为奇数，利用融合参数和整景遥感影像，计算得到单角度水深反演结果，并将其潮汐校正到理论深度基准面，之后得到水深段标识影像；第四步：多角度水深反演融合；将单角度水深反演结果、水深段标识影像和融合参数作为多角度水深反演融合的输入，逐像元开展融合，具体包括：(1)判断像元对应的水深段标识，对在同一水深段且该水深段内像元个数最多的进行下一步运算，否则，按照第(3)步进行运算；(2)判断像元的水深值，若有2个或2个以上的像元值相同，赋该值为最终的水深值，若不同，进行下一步运算；(3)判断像元在5×5邻域内是否为奇异值，选非奇异值点赋值为最终水深值，若非奇异点个数有2个或2个以上，进行下一步运算，若均为奇异点，按照第(5)步进行；(4)比较像元在当前水深段的平均相对误差，以最小误差对应的像元值为最终的水深值；(5)计算像元5×5邻域的中值，并以该邻域中值为新像元值，判断其水深段归属，若相同，按照第(4)步进行，否则进行下一步运算；(6)比较像元值所在水深段的可信度大小，以可信度大的像元值为最终水深值；第五步：水深反演精度验证；所述精度验证是利用检查点开展融合前单角度反演结果和融合后多角度反演结果的比较，水深反演精度验证完成后，将最终水深值作为遥感图像实际水深值输出数据。