阅读说明：本技术 一种基于卫星资料的海洋锋识别方法 (Ocean front identification method based on satellite data ) 是由 郭海龙 王宁 杜辉 史剑 张雪艳 于 2021-07-16 设计创作，主要内容包括：本发明公开了一种基于卫星资料的海洋锋识别方法,本发明方法在梯度阈值法部分,采用了基于概率密度函数的分区段阈值,尽可能地降低了主观阈值选取对识别结果的影响,提高了阈值选取的合理性；在锋生动力函数法部分,将锋生因子分解为不同动力学框架下的平流作用项、流场散度项、变形场项和外强迫项,提出动力学判定标准,并对像素点进行再次筛选,对阈值进行动力学订正,提高锋面识别的准确性。本发明结合了海洋锋要素特征及动力学特征,兼具阈值法的高计算效率特征和动力学法的科学可靠性特征,能够基于卫星资料实时进行海洋锋面识别；输入数据为卫星观测的海表温度场、盐度场、海表风场及动力学高度数据,具备易收集、高分辨率和全公开等特征。(The invention discloses a sea front identification method based on satellite data, which adopts a subsection threshold based on a probability density function in a gradient threshold method part, reduces the influence of subjective threshold selection on an identification result as much as possible, and improves the rationality of threshold selection; in the frontal generation power function part, the frontal generation factor is decomposed into a advection action item, a flow field divergence item, a deformation field item and an external forcing item under different dynamic frames, a dynamic judgment standard is put forward, pixel points are screened again, a threshold value is subjected to dynamic correction, and the accuracy of frontal surface identification is improved. The ocean front feature identification method combines the ocean front feature and the dynamics feature, has the high calculation efficiency feature of a threshold value method and the scientific reliability feature of the dynamics method, and can identify the ocean front in real time based on satellite data; the input data are sea surface temperature field, salinity field, sea surface wind field and dynamic height data observed by a satellite, and the sea surface wind field and dynamic height data have the characteristics of easy collection, high resolution, full disclosure and the like.)

1. A sea front identification method based on satellite data is characterized by comprising the following steps:

(1) the first part adopts a gradient threshold value method module, for a specific sea area, temperature and salinity data observed by a satellite are input, gradient distribution of the elements is obtained based on a gradient operator, a threshold value range is given based on an accumulated probability density function, when the element gradient is smaller than a maximum gradient threshold value of the non-frontal area, the pixel points are judged as the non-frontal area, when the element gradient is larger than a minimum gradient threshold value of the frontal area, the pixel points are judged as the frontal area, and when the element gradient is within the threshold value range, the pixel points are classified as possible frontal area pixel points; classifying the pixel point into a category corresponding to the maximum probability value by using a Bayes principle, judging based on a judgment standard 2, classifying the pixel point into a frontal region pixel point if a condition is met, and classifying the pixel point into a non-frontal region pixel point if the condition is not met; and outputting the identification result to the frontal generating power function module;

(2) and the second part adopts an edge generation power function module based on an edge generation tendency equation and comprises three types, namely an edge generation tendency equation, a half-turn edge generation tendency equation and a full-current edge generation tendency equation under the quasi-rotation frame, which are used for identifying ocean edges with different scales and forming mechanisms, and taking the planet scale front as an example, the edge generation tendency equation under the quasi-rotation frame is adopted, and the concrete form is shown as the following formula:

wherein the differential operatorFor gradient operators, F characterizes temperature, salinity, or densityAnd the like, the parameters of the system,horizontal component of the current field, F, for satellite altimeter data inversion_outCharacterizing the frontal contribution of the external forcing term to the factor; the method comprises the steps of decomposing an edge function to obtain contribution items of edge factors, namely a advection action item of a rotating flow field, a divergence action item of the flow field, a deformation action item and an external forcing item, judging a basic mechanism of the edge generation by utilizing contribution weights of the factors, identifying edge pixel points based on the edge factors, comparing element gradient values corresponding to the non-pixel points with a non-edge maximum gradient threshold value if the non-edge pixel points are in an edge factor high-value area, modifying the threshold value if the non-edge pixel points are less than the threshold value, re-judging the first part, and classifying the first part as the edge pixel points if the non-edge pixel points are not in the edge factor high-value area.

2. The method according to claim 1, wherein the method has both high computational efficiency characteristics of a threshold method and scientific reliability characteristics of a dynamics method, and the satellite observation data is one or more of temperature, salinity, dynamics altitude, sea surface wind field and dynamics altitude.

Technical Field

The invention relates to a sea front identification method, in particular to a sea front identification method based on satellite data.

Background

Ocean fronts are a phenomenon commonly found in the ocean and refer to narrow transition zones where ocean elements change dramatically. Ocean fronts are of various types and vary significantly in their origin of formation and mechanism of formation. For example, the temperature front refers to a transition zone with a large gradient where the temperature of seawater changes strongly in the horizontal direction, and is one of the most active areas of interaction of seawater and air due to the sharp change of the temperature on both sides of the front, often accompanied by strong momentum, heat and water vapor exchange; due to the frontal action, the ocean frontal area is often accompanied by strong secondary circulation, and is an area with abundant ocean primary productivity; in addition, the ocean front has great influence on underwater acoustic signal propagation, and is one of the small-scale ocean phenomena which are most concerned in military application. Therefore, accurately identifying the position of the ocean frontal area is a necessary premise for researching the characteristics of the ocean frontal area, such as strength, dynamic environment and the like.

The existing ocean front identification technology focuses on ocean temperature fronts and is mainly identified based on ocean front temperature gradient thresholds, and the method is simple and effective, has the advantages that the position of the temperature front in the ocean can be quickly identified, and has the defect that the selection of the temperature gradient thresholds is seriously influenced by subjective factors; based on this method, histogram thresholding was proposed by Cayula and Cornillon (1992, 1995) et al, which first obtained the location of the ocean temperature front based on satellite images; the specific scheme is shown in figure 1. Subsequently, Diehl develops the method, improves the application range and the recognition capability of the method, and effectively reduces the influence of subjective factors in the selection of the gradient threshold, but the method is limited by the temperature gradient in a statistical region, and does not consider the dynamic factors in the marine frontal process, such as the characteristics of a frontal flow field and the like.

Hideyuki Nakano et al (2018) recognize the ocean front in the northwest Pacific area from the angle of a kinetic factor based on satellite altimeter data, define that when the average kinetic height is a specific value, the position of the isoline is the position of the ocean front, and the theoretical basis is that the rotating velocity in the ocean front is high, namely the characteristics of the flow field, such as shear or convergence, are accompanied by the frontal process, and the specific technical scheme is shown in FIG. 2.

The core technology is the second step, the average speed is calculated by adopting an accumulative integral method, and the formula is as follows:

the average speed is the average speed, the average dynamics height value is the average dynamics height value, the distribution of the average speed about the dynamics height value is obtained based on the formula, the average dynamics height value corresponding to the local maximum value of the average speed is obtained by combining a local maximum value method, the ocean front corresponding to the average dynamics height value is determined based on the prior ocean front information, and the position of the average dynamics height contour line is the position of the ocean front.

The method has the advantages that the axis of the identified ocean front is continuous and smooth, the method is mainly suitable for the permanently or semi-permanently existing ocean front, the calculated amount is large, certain ocean front prior information needs to be provided, the position error of the ocean front in the identification result is large, the ocean front with a short life cycle is difficult to identify, and the identified ocean front has false information.

In summary, with the improvement of the resolution of the satellite data, the invention provides a sea front identification method combining element gradients and dynamic factors, which comprehensively considers the sea front characteristics of different life cycles to realize the sea front target identification function under different resolution requirements.

Disclosure of Invention

Aiming at the defects in the prior art, the invention aims to provide a method for identifying the ocean front based on satellite data, which realizes the real-time identification of the ocean front based on the satellite data; the accuracy of threshold selection in the gradient threshold method is improved; the accuracy of the ocean front identification result is improved.

In order to achieve the purpose, the invention is realized by the following technical scheme: a sea front identification method based on satellite data comprises the following steps:

1. the first part adopts a gradient threshold value method module, for a specific sea area, temperature and salinity data observed by a satellite are input, gradient distribution of the element is obtained based on a gradient operator, a threshold value range is given based on an accumulated probability density function, when the element gradient is smaller than a maximum gradient threshold value of the non-frontal area, the pixel point is judged as the non-frontal area, when the element gradient is larger than a minimum gradient threshold value of the frontal area, the pixel point is judged as the frontal area, and when the element gradient is within the threshold value range, the pixel point is classified as a possible frontal area pixel point; classifying the pixel point into a category corresponding to the maximum probability value by using a Bayes principle, judging based on a judgment standard 2, classifying the pixel point into a frontal region pixel point if a condition is met, and classifying the pixel point into a non-frontal region pixel point if the condition is not met; and outputting the identification result to the frontal generating power function module.

2. The second part adopts an edge generation power function module based on an edge generation tendency equation and is composed of three types, namely an edge generation tendency equation, a half-turn edge generation tendency equation and a full-current edge generation tendency equation under the quasi-rotation frame, which are used for identifying ocean edges with different scales and forming mechanisms, and taking the planet scale front as an example, the edge generation tendency equation under the quasi-rotation frame is adopted, and the concrete form is shown as the following formula:

wherein the differential operatorIs a gradient operator, F represents parameters such as temperature, salinity or density and the like,horizontal component of the current field, F, for satellite altimeter data inversion_outThe sharp contribution of the external forcing term to the factor is characterized. By decomposing the frontal function, the contribution term of the frontal factor can be obtained as the advection action term of the rotating flow field and the divergence of the flow fieldAnd judging a basic mechanism for generating the frontal surface by utilizing contribution weights of the factors, identifying frontal region pixel points based on the frontal region factors, comparing element gradient values corresponding to the non-pixel points with a non-frontal region maximum gradient threshold if the non-frontal region pixel points are in a frontal region with a large value of the frontal region factors, modifying the threshold if the non-frontal region pixel points are less than the threshold, and judging the first part again, otherwise, classifying the non-frontal region pixel points into the frontal region pixel points.

The invention has the following beneficial effects:

(1) the identification method combines the ocean front element characteristics and the dynamic characteristics, has high calculation efficiency, and can identify the ocean front in real time based on satellite data; the input data of the method are sea surface temperature field, salinity field, sea surface wind field and dynamic height data observed by a satellite, and the method has the characteristics of easy collection, high resolution, full disclosure and the like.

(2) In the method, in the part of a gradient threshold value method, aiming at the defect of the traditional single threshold value selection, the segmented threshold value based on the probability density function is adopted, so that the influence of subjective threshold value selection on the identification result is eliminated as much as possible, and the reasonability of the threshold value selection is improved.

(3) In the method, the frontal factor is decomposed into a advection action item, a flow field divergence item, a deformation field item and an external forcing item under different dynamic frames in the frontal dynamic function part, a judgment standard is put forward, pixel points are screened again, a threshold value is subjected to dynamic correction, and the accuracy of frontal surface identification is improved.

Drawings

The invention is described in detail below with reference to the drawings and the detailed description;

FIG. 1 is a flow chart of a technical scheme for Cayula and Cornillon ocean front identification in the background art of the present invention;

FIG. 2 is a flow chart of a Hideyuki Nakano ocean front identification technical scheme in the background art of the present invention;

fig. 3 is a flow chart of the ocean front identification of the present invention.

Detailed Description

In order to make the technical means, the creation characteristics, the achievement purposes and the effects of the invention easy to understand, the invention is further described with the specific embodiments.

Referring to fig. 3, the following technical solutions are adopted in the present embodiment: a sea front identification method based on satellite data comprises the following specific steps:

1. the first part adopts a gradient threshold value method module, for a specific sea area, temperature and salinity data observed by a satellite are input, gradient distribution of the element is obtained based on a gradient operator, a threshold value range is given based on an accumulated probability density function, when the element gradient is smaller than a maximum gradient threshold value of the non-frontal area, the pixel point is judged as the non-frontal area, when the element gradient is larger than a minimum gradient threshold value of the frontal area, the pixel point is judged as the frontal area, and when the element gradient is within the threshold value range, the pixel point is classified as a possible frontal area pixel point; classifying the pixel point into a category corresponding to the maximum probability value by using a Bayes principle, judging based on a judgment standard 2, classifying the pixel point into a frontal region pixel point if a condition is met, and classifying the pixel point into a non-frontal region pixel point if the condition is not met; and outputting the identification result to the frontal generating power function module.

2. The second part adopts an edge generation power function module based on an edge generation tendency equation and is composed of three types, namely an edge generation tendency equation, a half-turn edge generation tendency equation and a full-current edge generation tendency equation under the quasi-rotation frame, which are used for identifying ocean edges with different scales and forming mechanisms, and taking the planet scale front as an example, the edge generation tendency equation under the quasi-rotation frame is adopted, and the concrete form is shown as the following formula:

wherein the differential operatorIs a gradient operator, F represents parameters such as temperature, salinity or density and the like,for satellite altimeter dataInverted horizontal component of the current field, F_outThe sharp contribution of the external forcing term to the factor is characterized. The method comprises the steps of decomposing an edge function to obtain contribution items of edge factors, namely a advection action item of a rotating flow field, a divergence action item of the flow field, a deformation action item and an external forcing item, judging a basic mechanism of the edge generation by utilizing contribution weights of the factors, identifying edge pixel points based on the edge factors, comparing element gradient values corresponding to the non-pixel points with a non-edge maximum gradient threshold value if the non-edge pixel points are in an edge factor high-value area, modifying the threshold value if the non-edge pixel points are less than the threshold value, re-judging the first part, and classifying the first part as the edge pixel points if the non-edge pixel points are not in the edge factor high-value area.

Wherein the satellite observation data comprise temperature, salinity, dynamic altitude, sea surface wind field, dynamic altitude and the like.

This embodiment is based on the satellite data, realizes the real-time identification of ocean frontal: with the development of marine observation satellite data, the method can provide the thermodynamic parameters of the sea surface, such as temperature and salinity in real time; and corresponding dynamic elements such as flow field, sea surface wind field and the like can be provided at the same time, and based on the method, thermodynamic parameters and dynamic characteristics are combined to realize real-time identification of different types of ocean fronts, so that support is provided for scientific research and economic production.

The specific implementation mode improves the accuracy of threshold selection in the gradient threshold method: the method adopts the Bayesian principle to select the sea-front element gradient threshold, eliminates the influence of subjective factors in the threshold selection step by step, and provides data input for a subsequent dynamics detection module.

This embodiment improves the degree of accuracy of ocean front recognition result: the method adopts frontal dynamic factor detection as a post-processing part of the ocean front gradient identification module, and improves the accuracy of the ocean front identification result based on the frontal dynamic theory.

The foregoing shows and describes the general principles and broad features of the present invention and advantages thereof. It will be understood by those skilled in the art that the present invention is not limited to the embodiments described above, which are described in the specification and illustrated only to illustrate the principle of the present invention, but that various changes and modifications may be made therein without departing from the spirit and scope of the present invention, which fall within the scope of the invention as claimed. The scope of the invention is defined by the appended claims and equivalents thereof.