CN113850894B - Global vortex track dynamic visualization method based on leaf programs - Google Patents

Global vortex track dynamic visualization method based on leaf programs Download PDF

Info

Publication number
CN113850894B
CN113850894B CN202111427269.4A CN202111427269A CN113850894B CN 113850894 B CN113850894 B CN 113850894B CN 202111427269 A CN202111427269 A CN 202111427269A CN 113850894 B CN113850894 B CN 113850894B
Authority
CN
China
Prior art keywords
vortex
track
motion
movement
calculating
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202111427269.4A
Other languages
Chinese (zh)
Other versions
CN113850894A (en
Inventor
马纯永
郑杰
梁达
陈戈
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Ocean University of China
Original Assignee
Ocean University of China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Ocean University of China filed Critical Ocean University of China
Priority to CN202111427269.4A priority Critical patent/CN113850894B/en
Publication of CN113850894A publication Critical patent/CN113850894A/en
Application granted granted Critical
Publication of CN113850894B publication Critical patent/CN113850894B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T13/00Animation
    • G06T13/802D [Two Dimensional] animation, e.g. using sprites

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Processing Or Creating Images (AREA)

Abstract

A global vortex track dynamic visualization method based on a leaf xlet comprises the following steps: generating vortex motion data; drawing vortex dynamic tracks of different levels: calculating the service life of the vortex, grading the service life of the vortex according to the transparency of the color, and storing the service life of the vortex in an array; generating a vortex track through multithread calculation, drawing the vortex track through a VectorField visualization class based on a leaf let, and gradually rendering to form a data layer, namely completing visualization; and then carrying out visual analysis on the obtained vortex motion track dynamics. The invention provides a vortex track dynamic visualization algorithm on the basis of a flow field visualization effect. Compared with the traditional method for analyzing the vortex track motion, the method can more intuitively see the overall vortex motion and the detail parts of some vortex motions, and lays a foundation for the system analysis and research of the interrelation between the vortex motion and various environmental factors based on the global angle.

Description

Global vortex track dynamic visualization method based on leaf programs
Technical Field
The invention belongs to the technical field of marine observation, and particularly relates to a global vortex track dynamic visualization method based on a leaf let.
Background
Mesoscale vortices are widespread in the oceans of the world, and their movement plays a very important role in the transportation of global substances and energy, with profound effects on environmental and climatic changes in the ocean. In recent years, with the development of ocean observation technologies such as satellite altimeters, Argo buoys and the like, the identification and tracking method of the mesoscale vortex is continuously improved, and researchers at home and abroad increasingly study the vortex movement law. The current research methods of the vortex movement law mainly fall into two categories. The first type is an observation statistical analysis based on vortex trajectories: morrow et al systematically studied the vortex trajectories in the southeast part of the indian ocean, the southeast part of the atlantic ocean, and the northeast part of the pacific ocean and found that warm vortices tend to deflect in the equatorial direction, while cold vortices tend to deflect in the polar direction; chelton et al found that vortices move primarily westward and eastward in the south pole streaming region, but the meridional directions of warm and cold vortices move in opposite directions. The second type is to apply artificial intelligence method to do space-time analysis to the vortex track. Du et al applied a density-based trajectory spatiotemporal clustering method to perform spatiotemporal clustering analysis on the south-sea vortex trajectory to obtain the spatial distribution and temporal characteristics of the typical vortex movement pattern.
The visualization method is various, and the common three-dimensional vector field visualization methods are streamline tracking, line integral convolution and volume rendering; the basic principle is to extend the flow line from a two-dimensional plane to a three-dimensional plane. Common two-dimensional vector field visualization methods include an arrow mark method, a texture method and a vector line method. The method of the arrow mark method is simple and intuitive, but when the method is applied to a large-area with disordered directions, more vectors are disordered, and less vectors are free of details, so that the method is not suitable for sampling dense and changing abrupt vector fields; the texture method is that when the colors of the points in the texture are arranged according to a certain rule, the texture has a shape and can reflect the structure of the whole vector field, but the blurring and the confusion are easy to generate; the vector line method uses a curve trend in space to describe a vector field. In the vector field, the instantaneous velocity of all points on the streamline are tangent to the line. The trace is a vector line formed by the release of a particle, without mass, in a vector field along its trajectory in the flow field. The flow lines are suitable for studying the stable field and the traces are suitable for studying the unstable field, in which both are coincident. The instantaneous characteristics of the vortex motion are expressed in a streamline form in space, the continuous characteristics of the vortex motion are expressed in a track form in time, the density, the service life and the smoothness of the vortex motion can be controlled, and various characteristics of vectors can be displayed intuitively and efficiently.
Disclosure of Invention
The invention aims to provide a global vortex track dynamic visualization method based on a leaf, so as to make up for the defects of the prior art and realize the global vortex track dynamic visualization.
The research method of vortex motion mainly comprises the steps of comprehensively analyzing the motion characteristics of vortex according to the motion of all vortices based on a vortex tracking data set, mainly based on regional analysis, or carrying out systematic analysis on the motion process of individual vortex in certain regions. When the method is applied to a global layer surface, only the overall movement trend can be observed, and the method is not easy to be perceived in a micro area with obvious change. Therefore, the invention also discloses a flow field visualization method, which comprises the steps of preprocessing the vortex track data, generating vortex track vector field data which is consistent with the flow field visualization data format, associating the motion characteristics of ocean vortices, processing the vortex track vector field data according to the visualization principle, and realizing the visualization of the global vortex track dynamics.
In the traditional vortex motion analysis, most researches independently analyze the overall motion trend of a vortex in a certain area and the specific motion of a single vortex in the area according to the position information of a target vortex, but the researches on the global vortex motion trend and the motion detail change are lacked. The method calculates the average movement direction and speed of the vortex on each resolution grid, thereby obtaining refined vortex movement characteristics; secondly, initializing the position of the vortex on a two-dimensional plane, and defining the distribution density and the life cycle of the vortex; then, according to the movement direction and speed of the vortex in the grid, adopting a Longge Kutta function to calculate the next position to which the vortex moves; and (4) performing iterative calculation, and finally displaying the motion track of the vortex in a streamline form. The invention aims to comprehensively observe the movement of vortex based on global-regional level, and represents the movement direction of the vortex in the form of streamline.
The invention provides a vortex track dynamic visualization method on the basis of flow field visualization, wherein the motion direction of a vortex is calculated by utilizing a fourth-order Runge Kutta integral function according to the average motion of all vortices; the vortex motion is characterized by being instantaneous in the form of streamline in space and continuous in the form of track in time, and the density, the service life and the smoothness of the vortex motion can be controlled, so that various characteristics of vectors can be displayed intuitively and efficiently.
Based on the principle, in order to achieve the purpose, the invention adopts the following specific technical scheme:
a global vortex track dynamic visualization method based on a leaf let comprises the following steps:
(1) generating vortex motion data: performing gridding processing by extracting relevant information of vortex (extracting the service life of the vortex and the longitude and latitude of the vortex in each day), calculating the direction and the speed of the vortex motion, and generating vector field data and speed field data of the vortex motion;
(2) calculate the path of movement (trajectory) of the vortex: firstly, initializing vortex parameters, calculating the position of vortex motion by utilizing a fourth-order Runge Kutta function, and generating a vortex track through iterative calculation;
(3) drawing vortex dynamic tracks of different levels: calculating the service life of the vortex, grading the service life of the vortex according to the transparency of the color, and storing the service life of the vortex in an array;
(4) and (3) dynamically displaying vortex tracks: generating a vortex track by multithread calculation (3), drawing the vortex track through a VectorField visualization class based on the leaf, and gradually rendering to form a data layer, namely completing visualization;
(5) and carrying out visual analysis on the obtained vortex motion track dynamics.
Further, in the step (1), the vortex tracing data set of at least 20 years around the world is preprocessed, specifically:
(a) extracting relevant information of the vortex; extracting the service life and the longitude and latitude of each vortex in the vortex tracking data set respectively;
(b) gridding vortex information; setting the grid resolution of a two-dimensional plane to be 1 degree multiplied by 1 degree, regarding the position change of a vortex from the nth day to the (n + 1) th day as a straight line a, calculating grids passed by the straight line a, and determining the intersection point of the straight line a and each passed grid;
(c) calculating the movement direction of the vortex; calculating the movement direction of the vortex in each grid according to two intersection points of the straight line and the grid, wherein the movement direction is 0-360 degrees, the north direction is 0 degrees, and the movement direction is the movement direction of the corresponding grid through which the straight line a passes;
(d) calculating the movement speed of the vortex; the distance of the straight line a represents the movement displacement of the vortex every day, so that the distance of the straight line a represents the movement speed of the vortex from the nth day to the (n + 1) th day, and the speed is the speed of the grid passed by the straight line a;
(e) generating vortex motion direction vector field data; averaging the motion directions on each grid calculated in c), so as to obtain an average motion direction, and calculating a sine and cosine angle according to the average motion direction, wherein the sine angle is defined as a u component, and the cosine angle is defined as a v component, namely the vector field data of the vortex motion direction;
(f) generating vortex motion velocity field data; averaging the motion velocity on each grid calculated in d), vortex motion velocity field data can be obtained.
Further, the step (2) is specifically as follows:
(a) initializing vortex parameters; initializing the position of the vortex on a two-dimensional plane, and defining parameters such as the distribution density, the life cycle, the motion track color, the streamline transparency and the like of the vortex.
(b) Calculating the vortex motion position; according to the movement speed and direction of the grid where the vortex is located, if the movement speed of the vortex is not available in the current grid, the vortex stops moving, and a new vortex is randomly generated again at the position; otherwise, determining the coordinate of the next position of the vortex through a fourth-order Runge Kutta function, and determining the grid where the new position of the vortex is located according to the obtained coordinate;
(c) generating a vortex track; determining the number of times of vortex movement according to the life cycle of the vortex, iteratively calculating the position of the vortex movement until the vortex stops moving, and connecting the position of the vortex movement calculated each time into a track, namely a vortex track.
Further, the step (3) is specifically as follows:
(a) calculating the vortex life; determining the service life of the vortex according to the frequency of the vortex movement, wherein one time of the vortex movement is added to the service life of the position where the vortex moves;
(b) grading track colors; and grading the service life of the vortex according to the transparency of the color, wherein the service life of the vortex is stored in an array, the color corresponding to each item of the array is the color of the vortex in the service life, and the longer the service life is, the lighter the color of the vortex track is.
Further, the step (4) is specifically as follows:
(a) generating a track multithread; calculating a reasonable thread number k, calculating the number n/k of vortexes distributed by each thread according to the number n of vortexes defined by initialization, constructing multiple threads, and transmitting vortex motion vector data to the multiple threads, so as to realize the simultaneous generation of vortex tracks;
(b) drawing a track; judging whether all the vortex tracks are generated in real time, and if not, not drawing; otherwise, returning the vortex track data to the visual class, and drawing the vortex track according to the vortex position calculated on the vortex track and the service life of the vortex track;
(c) rendering a track; acquiring corresponding color values according to the vortex service lives of the vortex track at the current position and the next position, and performing color rendering on the vortex track in the interval;
(d) interface visualization is realized; and extracting a global map, and matching the longitude and latitude of the global map with the longitude and latitude of the vortex track layer, so that the rendered vortex track layer is loaded and drawn on the global map.
Further, the step (5) is specifically as follows:
(a) the vortex is divided into warm vortex and cold vortex; respectively processing the vortex track data of the cold vortex and the warm vortex according to the steps (1) to (4), thereby realizing the visualization of the tracks of the warm vortex and the cold vortex;
(b) dynamically visualizing the movement track of the regional vortex; and (3) selecting a region to be amplified, processing the range coordinate of the region to be amplified, reinitializing the vortex parameters in the range coordinate, and processing according to the steps (2) to (3), so that the dynamic visualization of the region vortex motion track is realized.
The invention has the advantages and beneficial effects that:
(1) the vortex motion is combined with a visualization principle, and dynamic display of a vortex motion track is achieved.
(2) The traditional analysis of the vortex track movement is based on the average movement direction of the vortex, and only the movement direction of the vortex track can be observed macroscopically, but the details of the vortex track movement can be seen by adopting a visual method, and the vortex can be found to be in a convoluted movement in some places.
(3) The visualization method is used for analyzing by combining the movement direction and the movement speed of the vortex track, the difference of the movement of the vortex at different places can be seen on a visualization interface, and the vortex is analyzed for a density field by comparing with the vortex speed, so that the visualization method is more intuitive.
(4) The movement condition of the vortex can be analyzed macroscopically according to the movement direction of the vortex, and the method has extremely important reference value for the research of the global vortex movement.
(5) By utilizing the invention, the global overall movement trend can be observed, the region can be enlarged, and the tiny change of the region can be observed; the movement trend of the global vortex is analyzed more intuitively, and for some detail changes, the movement of the vortex can be seen through an enlarged area.
The invention provides a vortex track dynamic visualization algorithm on the basis of a flow field visualization effect. Compared with the traditional method for analyzing the vortex track motion, the method can more intuitively see the overall vortex motion and the detail parts of some vortex motions, and lays a foundation for the system analysis and research of the interrelation between the vortex motion and various environmental factors based on the global angle.
Drawings
FIG. 1 is a flow chart of a global vortex trajectory dynamic visualization method based on a leaf xlet.
Fig. 2 is a graph showing a dynamic visualization of a track of global warm vortexes (AEs).
FIG. 3 is a dynamic visualization of the trajectory of global cold vortexes (CEs).
FIG. 4 is a dynamic partial magnified view of the vortex trajectory.
Detailed Description
The invention will be further explained and illustrated by means of specific embodiments and with reference to the drawings.
Example 1:
and generating motion direction vector field and speed field data by using the 20-year vortex motion data. Randomly generating initialization vortexes in a global range, enabling the initialization vortexes distributed densely to move along with a vortex track motion direction vector field, controlling the motion speed of the initialization vortexes through the motion speed of the vortex track, and expanding the initialization vortexes into a streamline through integration; and finally, changing the color depth of the streamline according to the number of grids through which the initial vortex motion passes, and realizing dynamic visualization of the global vortex track through the streamlines of different detail levels.
Referring to fig. 1, a global vortex trajectory dynamic visualization method based on a leaf let includes the specific steps:
1. vortex motion data is generated. Preprocessing a global 20-year vortex tracking data set, specifically:
(1) and meshing vortex information. Setting the grid resolution of the two-dimensional plane to 1 ° × 1 °, 180 × 360 grids around the world, vortex from day n ((n))
Figure 48474DEST_PATH_IMAGE001
) By day n + 1: (
Figure 152565DEST_PATH_IMAGE002
) The position change of the straight line a is regarded as a straight line a, the grids passed by the straight line a are calculated, and the intersection point of the straight line a and each passed grid is determined
Figure 869985DEST_PATH_IMAGE003
(2) The direction of vortex motion is calculated. From the coordinates of two points of intersection of a straight line with the grid
Figure 953610DEST_PATH_IMAGE004
Calculating the movement direction of the vortex in each grid
Figure DEST_PATH_IMAGE005
The moving direction is between 0 degrees and 360 degrees, wherein the north direction is 0 degrees, and the moving direction is the moving direction of the corresponding grid through which the straight line a passes. In which the direction of movement
Figure 564720DEST_PATH_IMAGE006
The calculation formula is as follows:
Figure 72668DEST_PATH_IMAGE007
(3) the speed of the vortex motion is calculated. The distance of the straight line a represents the movement displacement of the vortex every day, therefore, the distance of the straight line a represents the movement speed of the vortex from the nth day to the (n + 1) th day, and the speed is the speed of the grid passed by the straight line a, wherein the speed is calculated as follows:
Figure 277384DEST_PATH_IMAGE008
(4) vortex motion direction vector field data is generated. And (3) averaging the motion directions on each grid calculated in the step (2) to obtain an average motion direction, and calculating a sine and cosine angle according to the average motion direction, wherein the sine angle is defined as a U component, and the cosine angle is defined as a V component, namely the vortex motion direction vector field data.
Figure 741864DEST_PATH_IMAGE009
Figure 597693DEST_PATH_IMAGE010
Wherein
Figure 528740DEST_PATH_IMAGE011
Is the average direction of motion of each vortex in the (x, y) grid, and n is the number of vortex trajectories through the grid.
(5) Vortex motion velocity field data is generated. And (4) averaging the motion speed on each grid calculated in the step (3), and obtaining vortex motion speed field data.
(6) Storing the generated data, as shown in table 1:
TABLE 1 vortex motion velocity field data sheet
Figure 971485DEST_PATH_IMAGE013
Where Lon is the global longitude distribution and the parameter values are obtained with a resolution of 1 deg. over the world (between 0 deg. -360 deg.); lat is the distribution of the latitude around the world, and the parameter values are obtained with a resolution of 1 ° around the world (-90 ° -90 °). Two-dimensional data generated with a resolution of 1 ° × 1 ° on a two-dimensional plane of global extent (0 ° -360 °, -90 ° -90 °). The U component is a movement component of the vortex in the dimension direction, and the parameter value is a component of the movement direction calculated in a grid of 1 degree multiplied by 1 degree in the latitude direction according to the vortex; the V component is a movement component of the vortex in the longitudinal direction, and the parameter value is a component of the movement direction calculated in a grid of 1 degree multiplied by 1 degree in the longitudinal direction according to the vortex; speed represents the vortex motion Speed and the parameter values are determined from the calculated motion Speed of the vortex in a 1 x 1 grid.
2. The position of the swirling motion is calculated.
(1) And initializing vortex parameters and initializing the position of vortex on a two-dimensional plane. That is, an initial array is generated, the initial longitude and latitude of each vortex are stored in the array, the life cycle of each vortex is defined to be 50, the color of the motion track is white, the width of the vortex track is 1.2, the transparency is 0.99, and the vortex track disappears when the life of each vortex track is 20.
(2) The vortex motion position is calculated. According to the movement speed and the direction of the grid where the vortex is located, if the movement speed of the current grid is 0, the movement is stopped, and a new vortex is randomly generated again at the position; otherwise, determining the coordinate of the next position of the vortex through a fourth-order Runge Kutta function, and determining the grid where the new position of the vortex is located according to the obtained coordinate. The fourth-order Runge Kutta function is that the next position of the vortex is calculated by using the weighted average of four slopes, a step length streamline is constructed according to the movement time, and the real vortex movement track can be accurately fitted. The formula is as follows:
Figure 911759DEST_PATH_IMAGE014
(1)
Figure 622095DEST_PATH_IMAGE015
(2)
Figure 143950DEST_PATH_IMAGE016
(3)
Figure 854417DEST_PATH_IMAGE017
(4)
Figure 519753DEST_PATH_IMAGE018
(5)
wherein
Figure 100908DEST_PATH_IMAGE019
Is four slopes and h is the distance the vortex moves at this location.
(3) A vortex track is generated. And determining the number of times of vortex movement according to the life cycle of the vortex, iteratively calculating the position of the vortex movement, and stopping the vortex movement if the movement speed of the position of the vortex is 0 or the life of the vortex starting point is 50. And connecting the positions of each vortex motion calculation into a track after stopping the motion, namely the vortex track.
(4) A new vortex is regenerated. And (3) when the life of each vortex track defined in the step (2) is 20, the vortex track disappears, and when the vortex track at the position disappears, a new vortex is regenerated at the position, and the movement of the vortex iterates the above (1) - (4).
3. Drawing vortex dynamic tracks of different levels, specifically:
(1) vortex life was calculated. The service life of the vortex is determined according to the number of times of vortex movement, and the service life of the position where the vortex moves is increased by one time of vortex movement.
(2) The track is color graded. And grading the service life of the vortex according to the transparency of the color, and storing the service life of the vortex in an array, wherein the color corresponding to each item of the array is the color of the vortex in the service life. The longer the life, the lighter the vortex track color, and the vortex track disappears until the life of each vortex track is 20.
4. And generating vortex motion tracks of different detail levels.
(1) Calculating the vortex life; determining the service life of the vortex according to the frequency of the vortex movement, wherein one time of the vortex movement is added to the service life of the position where the vortex moves;
(2) grading track colors; grading the service life of the vortex according to the transparency of the color, and storing the service life of the vortex in an array, wherein the color corresponding to each item of the array is the color of the vortex in the service life, and the longer the service life is, the lighter the color of the vortex track is;
(3) therefore, it is possible to generate a swirling motion trajectory with different motion speeds, motion directions, and color changes. The white streamline represents the movement track of the vortex, and as the vortex moves, the service life of the vortex is prolonged, and the track color gradually becomes lighter until disappearing.
5. Vortex track dynamic display
(1) And generating multiple threads of the track. Firstly, detecting the hardware configuration of a machine, calculating a reasonable thread number k, calculating the number n/k of vortexes distributed by each thread according to the number n of vortexes defined by initialization, constructing multiple threads, and transmitting vortex motion vector data to the multiple threads, thereby realizing the simultaneous generation of vortex tracks;
(2) and (6) drawing a track. Judging whether all the vortex tracks are generated in real time, and if not, not drawing; otherwise, returning the vortex track data to the visual class, and drawing the vortex track according to the vortex position calculated on the vortex track and the service life of the vortex track;
(3) and rendering the track. Dividing the colors of the vortex track into 50 types, acquiring corresponding color values according to the vortex service life of the current position, and performing color rendering on the vortex track;
(4) and realizing interface visualization. And extracting a global map, and matching the longitude and latitude of the global map with the longitude and latitude of the vortex track layer, so that the rendered vortex track layer is loaded and drawn on the global map.
6. And (4) dynamically and visually analyzing the vortex motion track.
(1) The vortex trajectory is visually analyzed. The warm vortex trajectory visualization is shown in fig. 2, and the cold vortex trajectory visualization is shown in fig. 3. According to the visual graph of the vortex tracks of the cold vortex and the warm vortex, the main motion characteristics of the warm vortex and the cold vortex are similar, and the warm vortex and the cold vortex move in the western direction in a southern ocean area; the motion is east-oriented near the equator, and the motion has equatorial character; moves along the west direction of land near 60 degrees north latitude. The major movement trend of global vortex can be better analyzed.
And (4) dynamically and visually analyzing the regional vortex motion track. The effect of amplification of the swirling orbital motion in the region near australia is shown in figure 4. The vortex in this area can be seen globally with east and west motion, but its specific details are not visible; by enlarging the area, it can be seen that the vortex moves in south along the bank in australia, and the vortex moves in various directions near the island below australia. Through regional amplification, details which are not seen in the visualization of the global vortex track can be found, and the movement characteristics and detail changes of the vortex can be better analyzed.

Claims (5)

1. A global vortex track dynamic visualization method based on a leaf xlet is characterized by comprising the following steps:
(1) generating vortex motion data: performing gridding processing by extracting relevant information of vortex, calculating the direction and speed of vortex motion, and generating vector field data of the direction of the vortex motion and speed field data of the vortex motion; the step (1) is specifically as follows:
(a) extracting relevant information of the vortex; extracting the service life and the longitude and latitude of each vortex in the vortex tracking data set respectively;
(b) gridding vortex information; setting the grid resolution of a two-dimensional plane to be 1 degree multiplied by 1 degree, regarding the position change of a vortex from the nth day to the (n + 1) th day as a straight line a, calculating grids passed by the straight line a, and determining the intersection point of the straight line a and each passed grid;
(c) calculating the movement direction of the vortex; calculating the movement direction of the vortex in each grid according to two intersection points of the straight line and the grid, wherein the movement direction is 0-360 degrees, the north direction is 0 degrees, and the movement direction is the movement direction of the corresponding grid through which the straight line a passes;
(d) calculating the movement speed of the vortex; the distance of the straight line a represents the movement displacement of the vortex every day, so that the distance of the straight line a represents the movement speed of the vortex from the nth day to the (n + 1) th day, and the speed is the speed of the grid passed by the straight line a;
(e) generating vortex motion direction vector field data; averaging the motion directions on each grid calculated in the step (c) to obtain an average motion direction, and calculating a sine and cosine angle according to the average motion direction, wherein the sine angle is defined as a u component, and the cosine angle is defined as a v component, namely vortex motion direction vector field data;
(f) generating vortex motion velocity field data; averaging the motion velocities on each grid calculated in (d) to obtain vortex motion velocity field data;
(2) calculating the motion trajectory of the vortex): firstly, initializing vortex parameters, calculating the position of vortex motion by utilizing a fourth-order Runge Kutta function, and generating a vortex track through iterative calculation;
(3) drawing vortex dynamic tracks of different levels: calculating the service life of the vortex, grading the service life of the vortex according to the transparency of the color, and storing the service life of the vortex in an array;
(4) and (3) dynamically displaying vortex tracks: generating a vortex track by multithread calculation (3), drawing the vortex track through a VectorField visualization class based on the leaf, and gradually rendering to form a data layer, namely completing visualization;
(5) and then carrying out visual analysis on the obtained vortex motion track dynamics.
2. The global vortex trajectory dynamic visualization method based on a leaf xlet according to claim 1, wherein the step (2) is specifically:
(a) initializing vortex parameters; initializing the position of a vortex on a two-dimensional plane, and defining the distribution density, the life cycle, the motion track color and the streamline transparency parameters of the vortex;
(b) calculating the vortex motion position; according to the movement speed and direction of the grid where the vortex is located, if the movement speed of the vortex does not exist in the current grid, the movement is stopped, and a new vortex is randomly generated again at the position; otherwise, determining the coordinate of the next position of the vortex through a fourth-order Runge Kutta function, and determining the grid where the new position of the vortex is located according to the obtained coordinate;
(c) generating a vortex track; determining the number of times of vortex movement according to the life cycle of the vortex, iteratively calculating the position of the vortex movement until the vortex stops moving, and connecting the position of the vortex movement calculated each time into a track, namely a vortex track.
3. The global vortex trajectory dynamic visualization method based on a leaf xlet according to claim 1, wherein the step (3) is specifically:
(a) calculating the vortex life; determining the service life of the vortex according to the frequency of the vortex movement, wherein one time of the vortex movement is added to the service life of the position where the vortex moves;
(b) grading track colors; and grading the service life of the vortex according to the transparency of the color, wherein the service life of the vortex is stored in an array, the color corresponding to each item of the array is the color of the vortex in the service life, and the longer the service life is, the lighter the color of the vortex track is.
4. The global vortex trajectory dynamic visualization method based on a leaf xlet according to claim 1, wherein the step (4) is specifically:
(a) generating a track multithread; calculating a reasonable thread number k, calculating the number n/k of vortexes distributed by each thread according to the number n of vortexes defined by initialization, constructing multiple threads, and transmitting vortex motion vector data to the multiple threads, so as to realize the simultaneous generation of vortex tracks;
(b) drawing a track; judging whether all the vortex tracks are generated in real time, and if not, not drawing; otherwise, returning the vortex track data to the visual class, and drawing the vortex track according to the vortex position calculated on the vortex track and the service life of the vortex track;
(c) rendering a track; acquiring corresponding color values according to the vortex service lives of the vortex track at the current position and the next position, and performing color rendering on the vortex track in the interval;
(d) interface visualization is realized; and extracting a global map, and matching the longitude and latitude of the global map with the longitude and latitude of the vortex track layer, so that the rendered vortex track layer is loaded and drawn on the global map.
5. The global vortex trajectory dynamic visualization method based on a leaf xlet according to claim 1, wherein the step (5) is specifically:
(a) the vortex is divided into warm vortex and cold vortex; respectively processing the vortex track data of the cold vortex and the warm vortex according to the steps (1) to (4), thereby realizing the visualization of the tracks of the warm vortex and the cold vortex;
(b) dynamically visualizing the movement track of the regional vortex; and (3) selecting a region to be amplified, processing the range coordinate of the region to be amplified, reinitializing the vortex parameters in the range coordinate, and processing according to the steps (2) to (3), so that the dynamic visualization of the region vortex motion track is realized.
CN202111427269.4A 2021-11-29 2021-11-29 Global vortex track dynamic visualization method based on leaf programs Active CN113850894B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111427269.4A CN113850894B (en) 2021-11-29 2021-11-29 Global vortex track dynamic visualization method based on leaf programs

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111427269.4A CN113850894B (en) 2021-11-29 2021-11-29 Global vortex track dynamic visualization method based on leaf programs

Publications (2)

Publication Number Publication Date
CN113850894A CN113850894A (en) 2021-12-28
CN113850894B true CN113850894B (en) 2022-03-08

Family

ID=78982162

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111427269.4A Active CN113850894B (en) 2021-11-29 2021-11-29 Global vortex track dynamic visualization method based on leaf programs

Country Status (1)

Country Link
CN (1) CN113850894B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115964546B (en) * 2023-01-06 2023-05-30 中国海洋大学 Vortex migration channel extraction and visualization method based on edge binding
CN115826089B (en) * 2023-02-15 2023-05-12 中国海洋大学 East-west vortex boundary extraction method based on vortex tracking data

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20120137827A (en) * 2011-06-13 2012-12-24 한국과학기술원 Sub-particle scale turbulence simulation method and system for sph fluid, recording medium for the same
CN104867186A (en) * 2015-04-29 2015-08-26 中国海洋大学 GPU-based interactive ocean three-dimensional flow field dynamic visual algorithm
CN110781446A (en) * 2019-09-20 2020-02-11 中国海洋大学 Method for rapidly calculating average vorticity deviation of ocean mesoscale vortex Lagrange
CN110779500A (en) * 2019-11-14 2020-02-11 中国人民解放军国防科技大学 Mesoscale vortex detection method for incremental deployment sensor

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105844036B (en) * 2016-03-30 2018-12-21 中国科学院地理科学与资源研究所 A kind of vortex flow field dynamic and visual method based on random particles

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20120137827A (en) * 2011-06-13 2012-12-24 한국과학기술원 Sub-particle scale turbulence simulation method and system for sph fluid, recording medium for the same
CN104867186A (en) * 2015-04-29 2015-08-26 中国海洋大学 GPU-based interactive ocean three-dimensional flow field dynamic visual algorithm
CN110781446A (en) * 2019-09-20 2020-02-11 中国海洋大学 Method for rapidly calculating average vorticity deviation of ocean mesoscale vortex Lagrange
CN110779500A (en) * 2019-11-14 2020-02-11 中国人民解放军国防科技大学 Mesoscale vortex detection method for incremental deployment sensor

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
A daily global mesoscale ocean eddy dataset from satellite altimetry;James H. Faghmous 等;《SCIENTIFIC DATA》;20150609;1-17 *
基于传输函数的中尺度涡旋时空连续可视化;田丰林 等;《海洋学报》;20200930;第42卷(第9期);119-133 *

Also Published As

Publication number Publication date
CN113850894A (en) 2021-12-28

Similar Documents

Publication Publication Date Title
CN113850894B (en) Global vortex track dynamic visualization method based on leaf programs
Potter et al. From quantification to visualization: A taxonomy of uncertainty visualization approaches
Wang et al. Integration of Google Maps/Earth with microscale meteorology models and data visualization
CN107168516B (en) Global climate vector field data method for visualizing based on VR and gesture interaction technology
CN111639426B (en) Wind characteristic space-time process visualization method
Jenny et al. Interactive local terrain deformation inspired by hand-painted panoramas
Burger et al. Importance-driven particle techniques for flow visualization
Liu et al. A framework for interactive visual analysis of heterogeneous marine data in an integrated problem solving environment
Peng et al. Higher Dimensional Vector Field Visualization: A Survey.
Li et al. Meteovis: Visualizing meteorological events in virtual reality
Bachthaler et al. Space-time Visualization of Dynamics in Lagrangian Coherent Structures of Time-dependent 2D Vector Fields.
Zhang et al. An interactive 4D spatio-temporal visualization system for hydrometeorological data in natural disasters
Zahan et al. Contour line stylization to visualize multivariate information
Rocha et al. Multivariate Visualization of Oceanography Data Using Decals.
Hood et al. Inferring airflow across martian dunes from ripple patterns and dynamics
CN106097426B (en) A kind of weather radar data three-dimensional visualization space leaping acceleration light projection method
Neuroth et al. Scalable visualization of discrete velocity decompositions using spatially organized histograms
Cuntz et al. GPU-based Dynamic Flow Visualization for Climate Research Applications.
Li et al. GPU accelerated marine data visualization method
Zhang et al. Fast simulation method for ocean wave base on ocean wave spectrum and improved gerstner model with GPU
CN115964546B (en) Vortex migration channel extraction and visualization method based on edge binding
CN101937576A (en) Dynamic texture waterfall modeling method combined with multiple physical attributes
Ahmed et al. Contourdiff: Revealing differential trends in spatiotemporal data
Singh et al. Front and skeleton features based methods for tracking salinity propagation in the ocean
Tong et al. Crystal Glyph: Visualization of Directional Distributions Based on the Cube Map.

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant