基于决策融合的浅海水深多时相遥感影像反演方法技术

本发明专利技术公开了涉及基于决策融合的浅海水深多时相遥感影像反演方法，包括以下步骤：第一步：对多光谱遥感影像进行预处理，得到海表反射率；第二步：现场实测水深值获取及处理；第三步：单时相水深反演及云掩膜；第四步：多时相水深反演融合；第五步：水深反演精度验证；水深反演精度验证完成后，将最终水深值作为遥感图像实际水深值输出数据。与现有的遥感影像水深反演方法相比，本发明专利技术可以采用多云和高海况情况下获取的遥感影像，挖掘其中的水深数据，扩大了遥感图像的使用范围，实现多时相遥感数据的充分利用，提高了反演精度，尤其适用于浅水区域的海洋水深测量。

【技术实现步骤摘要】

本专利技术涉及一种海洋水深遥感测量方法，属于遥感
，尤其涉及基于决策融合的浅海水深多时相遥感影像反演方法。
技术介绍
海洋水深数据测量是保障船舶航行、开展港口码头和海洋工程建设、制定海岸和海岛相关规划的必要基础数据。与水深现场测量手段相比，遥感技术具有覆盖广、周期短、费用低、空间分辨率高的优势。自20世纪70年代以来，国内外开展了各种被动遥感水深反演模型的研究，常用的可见光水深反演模型主要包括分析模型、半分析半经验模型和统计模型。利用不同模型，近年来在河流、湖泊、水库、海岛和海岸带周边等水深测量领域进行了反演应用。水深可见光遥感反演是获取浅海复杂地形水深的有利解决办法，尤其可以反演获取船只无法靠近和难以进入区域的水深资料，是水深测量的一种重要手段。但由于模型难以兼顾物理机制和参数化，因此已有的可见光水深遥感反演模型精度再提高的空间有限。水深多时相遥感反演可以克服单一时相影像成像时环境条件的限制，有利于水深信息的提取，当前已有将多时相应用于水深遥感反演的研究工作，但多是利用多时相遥感影像开展空间信息的插补，没有在决策融合层面进行开发应用。而决策融合能够充分利用已有遥感影像资源和信息，为提高光学遥感水深反演精度提供了新途径。中国专利(授权公告号CN102176001B)公开了“一种基于透水波段比值因子的水深反演方法”，采用遥感图像获得的数据建立神经网络水深反演模型，通过人工神经网络的自适应学习性能和非线性映射能力，降低海洋水深测量的标准差，增强其适用性。但该方法对海洋水深的测量较为复杂，采用的仍然是单期影像作为探测数据源，因而在复杂情况下对浅海区域水深的探测效果有限。
技术实现思路
本专利技术提供了基于决策融合的浅海水深多时相遥感影像反演方法，用于解决现有技术中水深反演模型精度较低的问题，与现有的反演方法相比，基于决策融合的浅海水深多时相遥感影像反演方法可以采用多云和高海况情况下获取的遥感影像，挖掘其中的水深数据，提高了反演精度，尤其适用于浅水区域的海洋水深测量。基于决策融合的浅海水深多时相遥感影像反演方法，包括以下步骤：第一步：对多光谱遥感影像进行预处理，得到海表反射率；所述预处理包括对辐射亮度转换、大气校正和太阳耀斑去除；第二步：现场实测水深值获取及处理；获取实验区的水深数据和对应的经纬度坐标，通过潮汐表确认潮高值，将水深数据校正到理论深度基准面上，并根据采用的多时相影像获取时刻，分别对水深值进行对应影像的潮汐校正；第三步：单时相水深反演及云掩膜；对每景处理后的多光谱影像，进行分类处理得到云掩膜影像，将云及云影区赋值为1，非云区为0，作为多时相反演融合输入的一部分，根据水深控制点处水深与对应影像像元反射率值间的关系，统计回归出多波段模型在该景影像水深反演的参数，多波段模型公式如下，Z=A0+Σi=1nAiXi...(1)]]>其中，Z为水深，n为参与反演的波段个数，A0和Ai为待定系数，Xi＝Ln(ρi-ρsi)，ρi是第i波段反射率数据，ρsi是该波段深水处的反射率；将水深控制点分成数个的水深段，并计算出各水深段的平均相对误差和可信度，作为多时相水深反演融合的另一项输入，即融合参数，δk=1nΣi=1n|zi-zi′|zi...(2)]]>Lk0.5+{Σi=1nwj(Lki-0.5)α本文档来自技高网...

【技术保护点】
基于决策融合的浅海水深多时相遥感影像反演方法，其特征在于，包括以下步骤：第一步：对多光谱遥感影像进行预处理，得到海表反射率；所述预处理包括对辐射亮度转换、大气校正和太阳耀斑去除；第二步：现场实测水深值获取及处理；获取实验区的水深数据和对应的经纬度坐标，通过潮汐表确认潮高值，将水深数据校正到理论深度基准面上，并根据采用的多时相影像获取时刻，分别对水深值进行对应影像的潮汐校正；第三步：单时相水深反演及云掩膜；对每景处理后的多光谱影像，进行分类处理得到云掩膜影像，将云及云影区赋值为1，非云区为0，作为多时相反演融合输入的一部分，根据水深控制点处水深与对应影像像元反射率值间的关系，统计回归出多波段模型在该景影像水深反演的参数，多波段模型公式如下，Z=A0+Σi=1nAiXi---(1)]]>其中，Z为水深，n为参与反演的波段个数，A0和Ai为待定系数，Xi＝Ln(ρi‑ρsi)，ρi是第i波段反射率数据，ρsi是该波段深水处的反射率；将水深控制点分成数个的水深段，并计算出各水深段的平均相对误差和可信度，作为多时相水深反演融合的另一项输入，即融合参数，δk=1nΣi=1n|zi-zi′|zi---(2)]]>Lk=0.5+{Σi=1nwj(Lki-0.5)α}1α---(3)]]>其中，k均表示水深段，式2中，zi是第i个水深检查点的实测值，zi'为其反演值，n为水深检查点个数，式3中，Lk表示第k水深段的可信度，wj表示每个信源的相对重要程度，且满足其公式为代表该i景影像控制点的Kappa系数，n表示参与反演融合影像的个数，α为奇数；利用得到的参数和整景遥感影像，计算得单时相水深反演结果，并将其潮汐校正到理论深度基准面，之后得到水深段标识影像；第四步：多时相水深反演融合；将单时相水深反演结果、水深段标识影像、云掩膜和融合参数作为多时相水深反演融合的输入，逐像元开展融合，具体包括以下几个步骤；(1)判断是否为云及云影区，非云区的像元参与下一步融合运算，当非云区像元个数为1时，赋该值为最终的水深值，仅当参与下一步运算的像元数为0时，所有像元都继续下一步；(2)判断像元对应的水深段标识，对在同一水深段且该水深段内像元个数最多的进行下一步运算，否则，按照第(4)步进行运算；(3)判断像元的水深值，若有2个或2个以上的像元值相同，赋该值为最终的水深值，若不同，进行下一步运算；(4)判断像元在5×5邻域内是否为奇异值，选非奇异值点赋值为最终水深值，若非奇异点个数有2个或2个以上，进行下一步运算，若均为奇异点，按照第(6)步进行；(5)比较像元在当前水深段的平均相对误差，以最小误差对应的像元值为最终的水深值；(6)计算像元5×5邻域的中值，并以该邻域中值为新像元值，判断其水深段归属，若相同，按照第(5)步进行，否则进行下一步运算；(7)比较像元值所在水深段的可信度大小，以可信度大的像元值为最终水深值；第五步：水深反演精度验证；所述精度验证是利用检查点开展融合前单时相反演结果和融合后多时相反演结果的比较。...

【技术特征摘要】
1.基于决策融合的浅海水深多时相遥感影像反演方法，其特征在于，包括
以下步骤：
第一步：对多光谱遥感影像进行预处理，得到海表反射率；
所述预处理包括对辐射亮度转换、大气校正和太阳耀斑去除；
第二步：现场实测水深值获取及处理；
获取实验区的水深数据和对应的经纬度坐标，通过潮汐表确认潮高值，将
水深数据校正到理论深度基准面上，并根据采用的多时相影像获取时刻，分别
对水深值进行对应影像的潮汐校正；
第三步：单时相水深反演及云掩膜；
对每景处理后的多光谱影像，进行分类处理得到云掩膜影像，将云及云影
区赋值为1，非云区为0，作为多时相反演融合输入的一部分，根据水深控制点
处水...

【专利技术属性】
技术研发人员：张靖宇，马毅，张震，
申请(专利权)人：国家海洋局第一海洋研究所，
类型：发明
国别省市：山东;37