基于FUI分区和Ransac的浅海遥感水深反演方法技术

本发明专利技术涉及一种基于FUI分区和Ransac的浅海遥感水深反演方法，本发明专利技术首先使用多光谱遥感数据计算得到FUI水色指数，并基于FUI水色指数的分布，结合人工目视解译，对浅海地区进行分区；然后使用Ransac算法对每个分区的实地水深点数据集（实地水深点的水深值及其对应的R、G、B波段上的反射率）进行筛选，保留正确的数据，去除存在较大误差的数据；最后，使用保存下来的实地水深点数据集进行分区反演。本发明专利技术将分区思想与抗差估计思想有机地结合起来，在一定程度上提升了反演精度。本发明专利技术能够广泛地应用于浅海多光谱遥感水深反演领域。

【技术实现步骤摘要】
基于FUI分区和Ransac的浅海遥感水深反演方法
本专利技术涉及海洋遥感领域，具体涉及一种基于FUI分区和Ransac的浅海遥感水深反演方法。
技术介绍
水深是浅海地区的重要地理要素之一，对于海岸带开发与保护、海上航行、海洋科学研究等具有重要意义。遥感技术具有大尺度、低成本、重复观测等优势，其为浅海水深探测提供了一种行之有效的新手段。由于技术方法简单、获取数据相对容易等优势，利用多光谱遥感数据来反演水深最为普遍。但是，使用多光谱遥感数据反演浅海水深存在着两个难点。其一，由于浅海地区环境因素复杂，使用单一的半理论半经验模型无法很好地拟合整个区域。因此，众多学者提出了许多分区反演的方法，例如根据底质分区、根据反演残差分区等等。这些方法虽取得了较好的反演结果，但实施起来较为复杂，需要大量的先验知识及工作。其二，多光谱遥感数据和源于海图的实地水深点数据不可避免地存在着各种测量误差，甚至是粗差，这个问题目前还没有引起其他学者的过多关注和研究。上述难点严重影响了浅海多光谱遥感水深反演的精度、效率。
技术实现思路
本专利技术要解决的技术问题是针对以上不足，提供一种基于FUI分区和Ransac的浅海遥感水深反演方法。为解决以上技术问题，本专利技术采用以下技术方案：步骤1、获取目标水域的多光谱遥感数据，并对多光谱遥感数据进行预处理；依据多光谱遥感数据的成像时间，对源于海图的实地水深点进行潮汐校正，得到实地水深点潮汐校正后的水深值；步骤2、使用预处理后的多光谱遥感数据的R、G和B波段，利用CIE-XYZ色度系统中三原色RGB与三刺激值XYZ之间的转换关系，计算得到多光谱遥感数据的R、G和B波段对应的三刺激值X、Y和Z：将得到的刺激值X和Y归一化，得到（x，y）：在二维CIE色度空间里，将得到的（x，y）进行转换，得到以等能白光点为原点的坐标（x’,y’）：（x’,y’）向量与二维CIE色度空间中的x’轴负方向之间的顺时针夹角即为色度角α，然后根据FUI水色指数查找表得到FUI水色指数；步骤3、根据步骤2得到的FUI水色指数分布,结合人工目视解译，对不同水色的水域进行人工分区，得到若干个水域分区；步骤4、在步骤3得到的每个水域分区里，将步骤1的实地水深点潮汐校正后的水深值及其对应的R、G和B波段的反射率导出；每个水域分区的实地水深点潮汐校正后的水深值和对应的R、G和B波段上的反射率组成每个水域分区的实地水深点数据集；步骤5、在步骤3得到的每个水域分区里，使用蓝绿双波段对数线性模型、红蓝绿三波段对数线性模型和Stumpf对数转换比值模型三种反演模型中的任意一种，利用对应的Ransac算法对步骤4得到的实地水深点数据集进行筛选，保留数据误差小于误差阈值的实地水深点数据及其在R、G和B波段上的反射率，去除数据误差大于误差阈值的数据，得到筛选后的实地水深点数据集；步骤6、在步骤3得到的每个水域分区里，将上一步得到的实地水深点数据集分成拟合数据和验证数据，然后利用拟合数据对步骤5用到的反演模型进行拟合，得到水深反演模型；步骤7、将步骤1得到的预处理后的多光谱遥感数据的R、G和B波段输入上一步得到的水深反演模型，得到水深图。进一步的，所述步骤1中，预处理包括辐射校正和大气校正，根据需要选择是否包括几何校正和太阳耀斑校正；所述潮汐校正的方法为获取多光谱遥感数据成像的具体时间，查询此时刻浅海水域的潮高，然后对实地水深点进行潮汐校正。进一步的，所述步骤3中的分区方法为：人工目视解译，根据水域FUI水色指数的分布进行人工分区。进一步的，所述步骤5中，筛选时使用对应模型的Ransac算法，Ransac算法的主要参数设置方法如下：（1）采样次数。此参数需要根据具体问题确定，蓝绿双波段对数线性模型、红蓝绿三波段对数线性模型、Stumpf对数转换比值模型分别至少采样3、4和2次。（2）使用最小二乘法计算反演模型的参数。（3）内外点阈值。此参数的确定无法具体量化，只能根据具体问题的需求来确定。在较浅的分区中，选用0.5m作为阈值；在较深的分区中，选用1.0m作为阈值。（4）每组参数的支持度。每组参数对应模型的内点个数。（5）循环次数。此参数的确定过程较为复杂，使其数值很大，保证迭代的充分即可。进一步的，所述步骤6中，拟合反演模型使用的是最小二乘法，反演精度的评价指标包括平均绝对误差、平均相对误差、验证数据的真实水深值与拟合水深值之间的相关系数和箱线图工具。进一步的，所述步骤6中，拟合数据和验证数据的比例为2：1或3：1。本专利技术的有益效果为：在传感器的入瞳辐射中，包含着来自水体底部、水柱、水面、大气的贡献，但其中包含水深信息的只有底部反射；然而，底部反射不仅与水体深度有关，还与水体底质、水体光学性质有关。FUI水体颜色是太阳光照与水中物质相互作用的结果，在20m以浅的水域，FUI水体颜色还受底部反射的影响。研究表明，FUI水色指数涵盖了非常广泛的自然水体光学特征，与水体透明度、水体营养状态等有着密切的相关关系。因此，利用FUI水色指数对浅海水域进行事先分区，能够综合各种水质水况信息以及部分底质信息，更好地、更全面地进行分区。此外，基于FUI水色指数的分区方法较其他分区方法来讲，具有操作简单、不需要过多其他资料等优点。Ransac算法具有良好的抗粗差性能，能够将不利于水深反演的外点提前剔除出去，提升实地水深点数据集的准确性，进而提升水深反演的精度。本专利技术将两种策略结合起来运用于水深反演，充分考虑了水体的空间异质性、多光谱遥感数据和实地水深点数据的不确定性，显著提升了反演精度。下面结合附图和实施例对本专利技术进行详细说明。附图说明图1为本专利技术实施例的方法流程示意图；图2为本专利技术实施例的多光谱遥感影像图，同时展示了实地水深点的分布；图3为本专利技术实施例的FUI水色指数分布图；图4为本专利技术实施例中使用蓝绿双波段对数线性模型对应的Ransac算法筛选后实地水深点的分布图；图5为本专利技术实施例中使用红绿蓝三波段对数线性模型对应的Ransac算法筛选后实地水深点的分布图；图6为本专利技术实施例中使用Stumpf对数转换比值模型对应的Ransac算法筛选后实地水深点的分布图；图7为本专利技术实施例中使用蓝绿双波段对数线性模型对应的Ransac算法筛选后拟合数据、验证数据的分布图；图8为本专利技术实施例中使用红绿蓝三波段对数线性模型对应的Ransac算法筛选后拟合数据、验证数据的分布图；图9为本专利技术实施例中使用Stumpf对数转换比值模型对应的Ransac算法筛选后拟合数据、验证数据的分布图；图10为本专利技术实施例中使用蓝绿双波段对数线性模型进行反演得到的MAE、MRE折线图；图11为本专利技术实施例中使用红绿蓝三波段对数线性模型进行反演得到的MAE、MRE折线图；图12为本专利技术实施例中使用Stu本文档来自技高网...

【技术保护点】
1.基于FUI分区和Ransac的浅海遥感水深反演方法，其特征在于，包括以下步骤：/n步骤1、获取目标水域的多光谱遥感数据，并对多光谱遥感数据进行预处理；依据多光谱遥感数据的成像时间，对源于海图的实地水深点进行潮汐校正，得到实地水深点潮汐校正后的水深值；/n步骤2、使用预处理后的多光谱遥感数据的R、G和B波段，利用CIE-XYZ色度系统中三原色RGB与三刺激值XYZ之间的转换关系，计算得到多光谱遥感数据的R、G和B波段对应的三刺激值X、Y和Z,并将X和Y归一化得到（x，y），在二维CIE色度空间里，将（x，y）进行转换，得到以等能白光点为原点的坐标（x’,y’），原点到（x’,y’）对应的向量与二维CIE色度空间中的x’轴负方向之间的顺时针夹角即为色度角α，然后根据FUI水色指数查找表得到FUI水色指数；/n步骤3、根据步骤2得到的FUI水色指数分布,结合人工目视解译，对不同水色的水域进行人工分区，得到若干个水域分区；/n步骤4、在步骤3得到的每个水域分区里，将步骤1的实地水深点潮汐校正后的水深值及其对应的R、G和B波段的反射率导出；每个水域分区的实地水深点潮汐校正后的水深值和对应的R、G和B波段上的反射率组成每个水域分区的实地水深点数据集；/n步骤5、在步骤3得到的每个水域分区里，使用蓝绿双波段对数线性模型、红蓝绿三波段对数线性模型和Stumpf对数转换比值模型三种反演模型中的任意一种，利用对应的Ransac算法对步骤4得到的实地水深点数据集进行筛选，保留数据误差小于误差阈值的实地水深点数据及其在R、G和B波段上的反射率，去除数据误差大于误差阈值的数据，得到筛选后的实地水深点数据集；/n步骤6、在步骤3得到的每个水域分区里，将上一步得到的实地水深点数据集分成拟合数据和验证数据，然后利用拟合数据对步骤5用到的反演模型进行拟合，得到水深反演模型；/n步骤7、将步骤1得到的预处理后的多光谱遥感数据的R、G和B波段输入上一步得到的水深反演模型，得到水深图。/n...

【技术特征摘要】
1.基于FUI分区和Ransac的浅海遥感水深反演方法，其特征在于，包括以下步骤：
步骤1、获取目标水域的多光谱遥感数据，并对多光谱遥感数据进行预处理；依据多光谱遥感数据的成像时间，对源于海图的实地水深点进行潮汐校正，得到实地水深点潮汐校正后的水深值；
步骤2、使用预处理后的多光谱遥感数据的R、G和B波段，利用CIE-XYZ色度系统中三原色RGB与三刺激值XYZ之间的转换关系，计算得到多光谱遥感数据的R、G和B波段对应的三刺激值X、Y和Z,并将X和Y归一化得到（x，y），在二维CIE色度空间里，将（x，y）进行转换，得到以等能白光点为原点的坐标（x’,y’），原点到（x’,y’）对应的向量与二维CIE色度空间中的x’轴负方向之间的顺时针夹角即为色度角α，然后根据FUI水色指数查找表得到FUI水色指数；
步骤3、根据步骤2得到的FUI水色指数分布,结合人工目视解译，对不同水色的水域进行人工分区，得到若干个水域分区；
步骤4、在步骤3得到的每个水域分区里，将步骤1的实地水深点潮汐校正后的水深值及其对应的R、G和B波段的反射率导出；每个水域分区的实地水深点潮汐校正后的水深值和对应的R、G和B波段上的反射率组成每个水域分区的实地水深点数据集；
步骤5、在步骤3得到的每个水域分区里，使用蓝绿双波段对数线性模型、红蓝绿三波段对数线性模型和Stumpf对数转换比值模型三种反演模型中的任意一种，利用对应的Ransac算法对步骤4得到的实地水深点数据集进行筛选，保留数据误差小于误差阈值的实地水深点数据及其在R、G和B波段上的反射率，去除数据误差大于误差阈值的数据，得到筛选后的实地水深点数据集；
步骤6、在步骤3得到的每个水域分区里，将上一步得到的实地水深点数据集分成拟合数据和验证数据，然后利用拟合数...

【专利技术属性】
技术研发人员：宋妍，徐皓冉，李贝贝，
申请(专利权)人：中国地质大学武汉，
类型：发明
国别省市：湖北;42