一种基于遥感数据的切变锋识别方法技术

本发明专利技术提供了一种基于遥感数据的切变锋识别方法，该方法包括：首先用遥感数据反演得到海表悬浮体浓度分布的连续影像；根据A（前时刻）、B（后时刻）两个连续时刻的影像识别并计算主要水团色块的位移，确定切变锋范围；再计算B时刻海表悬浮体浓度梯度，得到最大梯度等值线分布，可能存在多条；最后根据识别的最大梯度等值线，和识别出的切变锋范围进行叠加，落入切变锋范围中的最大梯度等值线，即为切变锋锋带位置。经过实际验证，根据本发明专利技术提供的识别方法，能够有效且准确的识别切变锋。能够有效且准确的识别切变锋。能够有效且准确的识别切变锋。

【技术实现步骤摘要】
一种基于遥感数据的切变锋识别方法


[0001]本专利技术属于海洋遥感技术和图像识别
，具体是一种基于遥感数据的切变锋识别方法。

技术介绍

[0002]切变锋是一种海洋动力现象，是流向相反的水体界面或流向相同但流速相差极大的水体界面。切变锋处具有较强的物质浓度梯度，如悬浮体浓度梯度。切变锋对河口泥沙等物质的输运具有重要的阻隔作用，基于遥感观测对其开展的特征识别，不仅对河口切变锋领域的理论研究具有重要意义，对河口海岸物质输运与沉积的实践也具有重要的应用价值。
[0003]目前，切变锋的研究主要基于多船同步观测法和数值模拟方法，但船测覆盖的空间范围有限，数值模拟在河口区域仍存在一定误差。尚缺乏一种有效的、实测的、大范围的观测手段。而利用高时空分辨率的遥感资料是切变锋研究的最佳选择，但目前仍缺少基于遥感资料的切变锋识别方法。

技术实现思路

[0004]本专利技术的目的是提供一种基于遥感数据的切变锋识别方法，以弥补现有技术的不足。
[0005]本专利技术提供识别方法的原理是对遥感数据反演计算后的图像识别，首先通过遥感海表反射率数据进行反演计算，得到不同浊度分布区域的“水团色块”，由于“水团色块”的颜色和面积在短时间内具有较好的稳定性和连续性，则“水团色块”移动可以表征其随海流的移动，进而可通过“水团色块”的移动判断出海流的运动速度和方向，当两个毗邻的水团色块朝相反方向移动，则表明两个水团的流动方向相反，中间则存在切变锋带。由于切变锋是水流的剪切界面，并不能从卫星遥感影像中直接捕捉到水团移动，但由于河口切变锋处具有较高的物质浓度梯度，因此可以用悬浮体的浓度梯度来指示锋带的具体位置。
[0006]为达到上述目的，并给予上述识别原理，本专利技术采取的具体技术方案为：一种基于遥感数据的切变锋识别方法，该方法包括以下步骤：（1）利用遥感数据反演得到海表悬浮体浓度分布的连续影像；（2）根据A（前时刻）、B（后时刻）两个连续时刻的遥感影像，利用MCC（Maximum Cross Correlation）方法识别并计算主要水团色块的位移，确定切变锋范围；（3）计算B时刻海表悬浮体浓度梯度梯度，得到最大梯度等值线分布，可能存在多条；（4）根据步骤（3）识别的最大梯度等值线，和从步骤（2）识别出的切变锋范围进行叠加，落入切变锋范围中的最大梯度等值线，即为切变锋锋带位置。
[0007]上述切变锋锋带位置识别的标识可以用该梯度线的长度表示切变锋的长度，高值梯度线的运动方向表示切变锋的运动方向。
[0008]进一步的，所述步骤（2）中利用MCC方法识别并计算主要水团色块的位移，确定切变锋范围。
[0009]所述步骤（2）中水团色块位移的MCC识别方法具体为：对连续的A、B时刻遥感影像进行匹配，在A时刻影像中选定一个任意形状的区域（例如：像素的区域或者不规则形状区域），称为Source Block，记作S
sub
；在B时刻影像中选定具有相同形状的区域，称为Target Block，记作T
sub
。A时刻中每一个S
sub
都能在B时刻中选出多个T
sub
，将A时刻中的S
sub
与B时刻中所有的T
sub
逐一匹配，用相关系数CosSim来衡量S
sub
与T
sub
的相似程度，计算时使用的是图像像元的灰度值，CosSim通过下式计算：CosSim最大值表征匹配程度最好的两个区域，将此时的T
sub
作为B时刻中与A时刻中S
sub
相匹配的区域。由于水团色块的颜色和面积在短时间内具有较好的稳定性和连续性，可依据此方法识别并计算主要水团色块的位移。两个移动方向相反的水团色块，认为是位于切变锋的锋面两侧；这两个水团色块的相邻边界则为切变锋的范围。
[0010]进一步的，所述步骤（3）中海表悬浮体浓度梯度的具体计算方法为：采用Sobel算法，计算2个卷积因子，方向的卷积因子为： ，y方向的卷积因子为，将2个卷积因子分别与SSC相乘，即与SSC图像做平面卷积，分别得到x和y方向的SSC梯度，即：G
x = X
·
SSC，G
y = Y
·
SSC，考虑到相邻两点间的距离，因此还需要计算得到相邻两点的x方向的距离S
x
和y方向的距离S
y
，总的SSC梯度值为：。
[0011]本专利技术的优点和技术效果：本专利技术根据河口切变锋的定义及其物质输运的物理意义，从流速剪切和海流搬运悬浮体的角度识别切变锋。利用高时间分辨率的卫星遥感数据，基于MCC方法改进的水团色块移动判别法和物质浓度梯度计算方法，判别切变锋的位置、类型覆盖范围，通过大量数据获取其运动规律信息。
[0012]经过实际验证，根据本专利技术提供的识别方法，能够有效且准确的识别切变锋。
附图说明
[0013]图1为实施例1中黄河口悬浮体浓度在两个时间段的切变锋位置变化预测结果图，其中（a）为A时刻（2018年4月20日11:00）GOCI卫星数据反演的SSC结果图；（b）为下一时刻（B时刻，2018年4月20日12:00）GOCI卫星数据反演的SSC结果图；（c）为（b）经本专利技术计算获得的SSC梯度分布图；（d）为已有研究的可靠流场模拟结果图。
具体实施方式
[0014]以下通过具体实施例并结合附图对本专利技术进一步解释和说明。
[0015]实施例1：一种基于遥感数据的切变锋识别方法，该方法包括以下步骤：（1）利用遥感数据反演得到海表悬浮体浓度分布的连续影像；（2）根据A（前时刻）、B（后时刻）两个连续时刻的遥感影像，利用MCC（Maximum Cross Correlation）方法识别并计算主要水团色块的位移，确定切变锋范围；所述步骤（2）中水团色块位移的MCC识别方法具体为：对连续的A、B时刻遥感影像进行匹配，在A时刻影像中选定一个任意形状的区域（例如：像素的区域或者不规则形状区域），称为Source Block，记作S
sub
；在B时刻影像中选定具有相同形状的区域，称为Target Block，记作T
sub
。A时刻中每一个S
sub
都能在B时刻中选出多个T
sub
，将A时刻中的S
sub
与B时刻中所有的T
sub
逐一匹配，用相关系数CosSim来衡量S
sub
与T
sub
的相似程度，计算时使用的是图像像元的灰度值，CosSim通过下式计算：CosSim最大值表征匹配程度最好的两个区域，将此时的T
sub
作为B时刻中与A时刻中S
sub
相匹配的区域。由于水团色块的颜色和面积在短时间内具有较好的稳定性和连续性，可依据此方法识别并计算主要水团色块的位移。两个移动方向相反的水团色块，认为是位于切变锋的锋面两侧；这本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于遥感数据的切变锋识别方法，其特征在于，该方法包括以下步骤：（1）用遥感数据反演得到海表悬浮体浓度分布的连续影像；（2）根据A前时刻、B后时刻两个连续时刻的遥感影像，利用MCC方法识别并计算主要水团色块的位移，确定切变锋范围；（3）计算B时刻海表悬浮体浓度梯度，得到最大梯度等值线分布，可能存在多条；（4）根据步骤（3）识别的最大梯度等值线，和从步骤（2）识别出的切变锋范围进行叠加，落入切变锋范围中的最大梯度等值线，即为切变锋锋带位置。2.如权利要求1所述的切变锋识别方法，其特征在于，所述步骤（2）中MCC方法具体为：对连续的A、B时刻遥感影像进行匹配，在A时刻影像中选定一个任意形状的区域，称为Source Block，记作S
sub
；在B时刻影像中选定具有相同形状的区域，称为Target Block，记作T
sub
；A时刻中每一个S
sub
都能在B时刻中选出多个T
sub
，将A时刻中的S
sub
与B时刻中所有的T
sub
逐一匹配，用相关系数CosSim来衡量S
sub
与T
...

【专利技术属性】
技术研发人员：王楠，宋德海，毛政，李广雪，鲍献文，
申请(专利权)人：中国海洋大学，
类型：发明
国别省市：