CN106202522B - A kind of multiplexing method and system of flow field integral curve - Google Patents
A kind of multiplexing method and system of flow field integral curve Download PDFInfo
- Publication number
- CN106202522B CN106202522B CN201610584306.5A CN201610584306A CN106202522B CN 106202522 B CN106202522 B CN 106202522B CN 201610584306 A CN201610584306 A CN 201610584306A CN 106202522 B CN106202522 B CN 106202522B
- Authority
- CN
- China
- Prior art keywords
- integral curve
- space
- flow field
- curve
- dimensional
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/90—Details of database functions independent of the retrieved data types
- G06F16/903—Querying
- G06F16/9038—Presentation of query results
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/90—Details of database functions independent of the retrieved data types
- G06F16/901—Indexing; Data structures therefor; Storage structures
- G06F16/9027—Trees
Landscapes
- Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computational Linguistics (AREA)
- Software Systems (AREA)
- Processing Or Creating Images (AREA)
Abstract
The invention discloses the multiplexing method and system of a kind of flow field integral curve, belongs to and be related to flow field data analysis technique field.This method comprises: obtaining the spatial mesh structure of flow field data；Initial seed point is placed on all mesh points of spatial mesh structure, calculates the integral curve of each seed point, and the integral curve is stored in storage equipment；The integral curve is trace or streamline；The flow field data are analyzed or are visualized with the integral curve in application, the initial seed point seed point at corresponding grid point locations required for calling directly in the storage device.This method and system are by precalculating enough integral curves and preserving, in the analysis of flow field data or visualization application later, regional space directly specified according to flow field analysis or visualization application, required integral curve can be rapidly obtained, treatment effeciency is greatly improved, the analysis and visualization processing of extensive flow field data are particularly suitable for.
Description
Technical field
The present invention relates to flow field data analysis technique fields, and in particular to a kind of multiplexing method of flow field integral curve and is
System.
Background technique
Space occupied by the fluid motion of flow field (Flow Field) i.e., flow field data are a kind of important science data,
It is typically defined in two dimension or three-dimensional space, and contains at least one velocity field.The most common flow field data include ocean
The wind field data etc. in ocean current data, atmospheric simulation in simulation.Many existing flow field visualization methods are all based on integral
Curve, for example, integral curve it is direct rendering, source converge inquiry, FTLE (finite time lyapunov index field) calculating, with
And some other visualization and analysis method specially designed.In Unsteady Flow (i.e. velocity field changes over time), most often
The integral curve seen be trace (in flow field certain position place a massless particle (referred to as initial seed point), allow its
It is moved in velocity field, being formed by track is exactly a trace)；In steady flow field (i.e. velocity field does not change over time), then
Correspond to streamline (all tangent curve is known as streamline with velocity vector on every bit in flow field).It is practical to calculate, trace and
Streamline is acquired by numerical integration, so they are referred to as integral curve.
Common flow field visualization method places seed point generally according to needs with position in different times, and calculates big
The trace of amount.For biggish flow field data, this calculating process devotes a tremendous amount of time cost, and resulting due to calculating
Trace data amount is too big, it will usually abandon the trace as intermediate result, more cause the waste of computing resource in this way.In order to
Overcome the deficiencies in the prior art, the present invention provide the multiplexing method and system of a kind of flow field integral curve.
Summary of the invention
In view of the deficiencies in the prior art, the purpose of the present invention is to provide a kind of multiplexing sides of flow field integral curve
Method and system, by this method and system, integral curve is obtained when can effectively improve the analysis of flow field data or visualization application
Take efficiency.
To achieve the above object, the technical solution adopted by the present invention are as follows: a kind of multiplexing method of flow field integral curve, including
Following steps:
(1) spatial mesh structure of flow field data is obtained；
(2) initial seed point is placed on all mesh points of the spatial mesh structure, calculates the product of each seed point
Component curve, and the integral curve is stored in storage equipment；The integral curve is trace or streamline；
(3) needed for the flow field data being analyzed or are visualized in application, calling directly in the storage device
The integral curve for the initial seed point at correspondence grid point locations wanted.
Further, a kind of multiplexing method of flow field integral curve as described above, in step (1), the space lattice knot
Structure is regular grid, when the spatial mesh structure of the flow field data of acquisition is irregular grid, on the basis of irregular grid
On, construct the rule space network of flow field data；The regular grid is two-dimensional regular grid or three-dimensional regular grid.
Further, a kind of multiplexing method of flow field integral curve as described above, in step (2), by the integral curve
It is stored in storage equipment, comprising:
2.1) area of space of the spatial mesh structure is divided, is divided into multiple subspaces；
2.2) index for establishing every sub-spaces, the corresponding integral of initial seed point included in every sub-spaces is bent
Line is associated with the index of its subspace to be stored in storage equipment.
Further, a kind of multiplexing method of flow field integral curve as described above, when the spatial mesh structure is two dimension
When regular grid, the integral curve is stored in storage equipment, comprising:
1) the four of two-dimensional space grid are established as the root node of quaternary tree using the whole region of the two-dimensional regular grid
Tree construction is pitched, the space of the two-dimensional regular grid is divided into multiple two-dimensional sub-spaces, each leaf section of quad-tree structure
The corresponding two-dimensional sub-spaces of point；
2) index for establishing each leaf node of quad-tree structure, corresponding two dimension of each leaf node is empty
Between included the corresponding integral curve of initial seed point and its node index be associated with to be stored in and stores in equipment；
When the spatial mesh structure is three-dimensional regular grid, the integral curve is stored in storage equipment, comprising:
1. establishing the eight of three-dimensional regular grid as the root node of Octree using the whole region of the three-dimensional regular grid
Tree construction is pitched, the space of the three-dimensional regular grid is divided into multiple three n-dimensional subspace ns, each leaf of octree structure
Node corresponds to three n-dimensional subspace ns；
2. the index of each leaf node in octree structure is established, by of three-dimensional corresponding to each leaf node
The corresponding integral curve of initial seed point included in space is associated with its node index and is stored in storage equipment.
Further, a kind of multiplexing method of flow field integral curve as described above, in step (2), by initial seed point pair
It is stored in storage equipment after the integral curve compression answered.
Further, a kind of multiplexing method of flow field integral curve as described above, by the space of the spatial mesh structure
Region division is multiple subspaces, comprising:
A., the compression ratio threshold value of integral curve is set；
B. the spatial mesh structure is divided for the first time, obtains several sub-spaces；
C. every sub-spaces are again divided into M sub-spaces；
D. for the subspace that each is divided again, initial seed point included in the subspace before division is corresponding
Integral curve compressed, compressed byte number is denoted as N1, by initial seed point included in the subspace M after division
Corresponding integral curve is compressed respectively, and the sum of M compression postbyte number is denoted as N2, judges whether the ratio of N1 and N2 is less than
Or it is equal to compression ratio threshold value, it completes if so, the space of the two-dimensional regular grid or three-dimensional regular grid divides, divides again
Preceding division result is final division result, if it is not, then return step c.
Further, a kind of multiplexing method of flow field integral curve as described above, in step (3), in the storage equipment
In call directly required for corresponding mesh point corresponding position initial seed point integral curve, comprising:
3.1) it determines and is analyzed or visualized using specified flow field data space region, and in the space region
The space density of initial seed point is placed in domain；
3.2) subspace intersected with the area of space is searched in all subspaces of spatial mesh structure；
3.3) search what association in storage equipment saved according to the index of the subspace intersected with the area of space
Integral curve；
3.4) corresponding integral curve is extracted in the integral curve found according to the space density.
A kind of multiplex system of flow field integral curve, comprising:
Spatial mesh structure obtains module, for obtaining the spatial mesh structure of flow field data；
Integral curve computing module, for placing initial seed point on all mesh points of the spatial mesh structure,
Calculate the integral curve of each seed point；The integral curve is trace or streamline；
Integral curve preserving module, for the integral curve to be stored in storage equipment；
Integral curve calling module, for the flow field data to be analyzed or visualized in application, in the storage
The integral curve of initial seed point at corresponding grid point locations required for being called directly in equipment.
Further, a kind of multiplex system of flow field integral curve as described above, the integral curve preserving module include:
Space lattice division unit is divided into multiple for dividing the area of space of the spatial mesh structure
Subspace；
Curve associated storage unit will be initial included in every sub-spaces for establishing the index of every sub-spaces
The corresponding integral curve of seed point is associated with the index of its subspace to be stored in storage equipment.
Further, a kind of multiplex system of flow field integral curve as described above, when the spatial mesh structure is two dimension
When regular grid, the space lattice division unit, comprising:
First grid dividing subelement is built for the root node using the whole region of two-dimensional regular grid as quaternary tree
The space of the two-dimensional regular grid is divided into multiple two-dimensional sub-spaces, four forks by the quad-tree structure of vertical two-dimensional space grid
The corresponding two-dimensional sub-spaces of each leaf node of tree construction；
The curve associated storage unit includes:
First associated storage subelement, the index of each leaf node for establishing quad-tree structure, by each
The corresponding integral curve of initial seed point that the corresponding two-dimensional sub-spaces of leaf node are included is associated with preservation with its node index
In storage equipment；
When the spatial mesh structure is three-dimensional regular grid, the space lattice division unit, comprising:
Second grid dividing subelement, for the root section using the whole region of the three-dimensional regular grid as Octree
Point establishes the octree structure of three-dimensional regular grid, and the space of the three-dimensional regular grid is divided into multiple three n-dimensional subspace ns,
Corresponding three n-dimensional subspace ns of each leaf node of octree structure；
The curve associated storage unit includes:
Second associated storage subelement, for establishing the index of each leaf node in octree structure, by each
The corresponding integral curve of initial seed point included in three n-dimensional subspace ns corresponding to leaf node is associated with its node index
It is stored in storage equipment.
Further, a kind of multiplex system of flow field integral curve as described above, the integral curve preserving module are also wrapped
It includes:
Integral curve compression unit, for storage equipment will to be stored in after the corresponding integral curve compression of initial seed point
In.
Further, a kind of multiplex system of flow field integral curve as described above, the space lattice division unit include:
First division subelement obtains several sub-spaces for being divided the spatial mesh structure for the first time；
Subelement is divided again, and every sub-spaces for dividing first division subelement are again divided into M
Subspace；
Division result determines subelement, the subspace for dividing again for each, judges the son before division is empty
Between included in the corresponding integral curve of initial seed point carry out compressed byte number N1 with will divide after the subspace M in
The corresponding integral curve of the initial seed point for being included carries out the ratio of the sum of compressed M compressed byte numbers N2 respectively
Whether it is less than or equal to compression ratio threshold value, is completed if so, the space of the two-dimensional regular grid or three-dimensional regular grid divides,
Division result before dividing again is final division result, divides subelement again if it is not, then entering.
Further, a kind of multiplex system of flow field integral curve as described above, the integral curve calling module include:
Area determination unit to be analyzed is analyzed or is visualized using specified flow field data space area for determining
Domain, and in the area of space place initial seed point space density；
Subspace determination unit, for searching the son intersected with the area of space in the subspace of spatial mesh structure
Space；
Integral curve searching unit, for searching storage according to the index of the subspace intersected with the area of space
The integral curve of preservation is associated in equipment；
Integral curve extraction unit, it is corresponding for being extracted in the integral curve found according to the space density
Integral curve.
The beneficial effects of the present invention are: the multiplexing method and system of flow field integral curve provided by the present invention pass through
It precalculates enough integral curves and preserves, in the analysis of flow field data or visualization application later, direct root
According to the specified regional space of flow field analysis or visualization application, it can be rapidly obtained in corresponding subspace and correspond to seed point
Integral curve, the greatly less acquisition time of integral curve, improves treatment effeciency, is particularly suitable for extensive flow field number
According to analysis and visualization processing.
Detailed description of the invention
Fig. 1 is a kind of flow chart of the multiplexing method of the flow field integral curve provided in specific embodiment；
Fig. 2 is the flow chart that the spatial mesh structure of flow field data in specific embodiment divides；
Fig. 3 is a kind of structural block diagram of the multiplex system of the flow field integral curve provided in specific embodiment；
Fig. 4 is the structural block diagram of integral curve preserving module in specific embodiment；
Fig. 5 is the structural block diagram of space grid dividing unit in specific embodiment；
Fig. 6 is the placement schematic diagram of three-dimensional regular network and initial seed point in embodiment.
Specific embodiment
The present invention is described in further detail with specific embodiment with reference to the accompanying drawings of the specification.
The multiplexing method and system of flow field integral curve provided by the present invention, basic ideas are to precalculate subsequent needs
Enough integral curves of the flow field data of analysis simultaneously store, in analysis or visualization later, directly from being deposited
Data in extract required integral curve.
Fig. 1 shows a kind of flow chart of the multiplexing method of the flow field integral curve provided in present embodiment, You Tuzhong
As can be seen that this method mainly may include the following three steps:
Step S100: the spatial mesh structure of flow field data is obtained；
For the subsequent flow field data for being analyzed or being visualized application, its spatial mesh structure is obtained first.
In general, space defined in the data of flow field is typically expressed as two dimension or three-dimensional regular grid, such as rectilinear grid
(rectilinear grid), rectangular mesh (erctilinear grid) etc., it is also possible to be curvilinear grid
(curvilinear grid) or unstrctured grid (unstructured grid).
In present embodiment, the spatial mesh structure is regular grid, when the space lattice knot of the flow field data of acquisition
When structure is irregular grid, on the basis of irregular grid, the rule space network of flow field data is constructed；The rule
Grid is two-dimensional regular grid or three-dimensional regular grid.
If flow field data inherently regular network, can be directly used its regular grid structure,
If the network of flow field data itself is irregular grid, as long as the space where its data itself is three-dimensional space,
The regular grid of a set of flow field data can be constructed, specific building mode uses the space lattice building mode of existing flow field data
.
Step S200: placing initial seed point on all mesh points of the spatial mesh structure of flow field data, calculates every
The integral curve of a seed point and preservation；
Initial seed point is placed on all mesh points of the spatial mesh structure of flow field data, calculates each seed point
Integral curve, and the integral curve is stored in storage equipment；The integral curve is trace or streamline.
The integral curve is the trace or streamline of initial seed point, Unsteady Flow (i.e. velocity field changes over time)
In, the most common integral curve is trace, in steady flow field (i.e. velocity field does not change over time), then corresponds to streamline.Institute
Initial seed point is stated, i.e., the particle of a massless is placed in certain position in flow field, allows it to move in velocity field, is formed by
Track is exactly a trace or streamline.
In present embodiment, the integral curve for the initial seed point densely placed in space defined in the data of flow field,
Track should be grown as far as possible, until seed point is walked out defined space and stopped.The calculation method of integral curve uses existing trace
Or streamline calculation or directly generated using existing visualization tool, in practical applications, when flow field data compared with
Hour, it can be calculated in single computer, it when the data is larger, can be in parallel computation under cluster environment.
In present embodiment, the integral curve of all initial seed points is stored in the mode of storage equipment are as follows:
1) area of space of the spatial mesh structure is divided, is divided into multiple subspaces；
2) index for establishing every sub-spaces, by the corresponding integral curve of initial seed point included in every sub-spaces
It is associated with and is stored in storage equipment with the index of its subspace.
Through the storage mode, in subsequent lookup integral curve, it is only necessary to can stored according to the index of subspace
Corresponding section component curve is found in equipment.
For different spatial mesh structures, present embodiments provide for two different integral curve preserving types.When
When the spatial mesh structure is two-dimensional regular grid, the integral curve is stored in the concrete mode of storage equipment are as follows:
1) the four of two-dimensional space grid are established as the root node of quaternary tree using the whole region of the two-dimensional regular grid
Tree construction is pitched, the space of the two-dimensional regular grid is divided into multiple two-dimensional sub-spaces, each leaf section of quad-tree structure
The corresponding two-dimensional sub-spaces of point；
2) index for establishing each leaf node of quad-tree structure, corresponding two dimension of each leaf node is empty
Between included the corresponding integral curve of initial seed point and its node index be associated with to be stored in and stores in equipment.
When the spatial mesh structure is three-dimensional regular grid, the integral curve is stored in the specific of storage equipment
Mode are as follows:
1. establishing the eight of three-dimensional regular grid as the root node of Octree using the whole region of the three-dimensional regular grid
Tree construction is pitched, the space of the three-dimensional regular grid is divided into multiple three n-dimensional subspace ns, each leaf of octree structure
Node corresponds to three n-dimensional subspace ns；
2. the index of each leaf node in octree structure is established, by of three-dimensional corresponding to each leaf node
The corresponding integral curve of initial seed point included in space is associated with its node index and is stored in storage equipment.
Since the integral curve of calculative initial seed point is usually primary flow field data hundreds of or thousands of times, it is
Guarantee all flow field integral curves can store, and in present embodiment, set integral curve is stored in storage
When standby, compression processing carried out to the corresponding integral curve of initial seed point first, then by the rope of compressed result and subspace
Draw (or node index) associated storage.
In practical applications, the compress mode of specific integral curve can be selected according to the actual situation.For example, can
It is compressed together with and by the corresponding integral curve of initial seed point included in every sub-spaces, it can also be in subspace
The corresponding integral curve of each initial seed point for being included is compressed one by one, is pressed again after compressed result is merged
Contracting.Specific compression algorithm also can according to need selection, and the compression that parameter curve is fitted integral curve can be used
Method, as Bezier is fitted to；It is also possible to the compression algorithm for any floating number array, such as fpzip algorithm, certainly
It is also possible to a variety of compression algorithms to be applied in combination.Exist to reduce the corresponding initial seed point of integral curve as caused by compression
Location error in grid generally requires the setting compressed error upper limit when carrying out curve compression.
Division for the spatial mesh structure of flow field data, if the number of subspace is fewer, institute in every sub-spaces
The initial seed point for including is more, and corresponding integral curve is also more, compression effectiveness better (compressed word when being compressed
Joint number is less), but in the integral curve that subsequent calls need, it can be slower when decompression, in order to balance these two aspects,
A kind of division mode of spatial mesh structure is additionally provided in present embodiment, which is as shown in Fig. 2, process is as follows:
A., the compression ratio threshold value of integral curve is set；The compression ratio threshold value is empirical value.
B. the spatial mesh structure is divided for the first time, obtains several sub-spaces；
C. every sub-spaces are again divided into M sub-spaces；
D. for the subspace that each is divided again, initial seed point included in the subspace before division is corresponding
Integral curve compressed, compressed byte number is denoted as N1, by initial seed point included in the subspace M after division
Corresponding integral curve is compressed respectively, and the sum of M compression postbyte number is denoted as N2, judges whether the ratio of N1 and N2 is less than
Or it is equal to compression ratio threshold value, it completes if so, the space of the two-dimensional regular grid or three-dimensional regular grid divides, divides again
Preceding division result is final division result, if it is not, then return step c.
Which passes through the corresponding product of all initial seed points included in the space before judging sub-spaces division
The compressed byte number of component curve compresses respectively with by the corresponding integral curve in all subspaces corresponding after the Subspace partition
Relationship between the ratio of the sum of several byte numbers afterwards and the compression ratio threshold value of setting determines whether to need in antithetical phrase sky
Between continue to divide.
For example, when being divided for above-mentioned three regular grids using octree structure, it first will three-dimensional grid really
Space is divided into 8 the first subspaces, is denoted as space S 1-S8 respectively, judges whether to need to 8 later through the above way
One subspace is divided again, such as space S 1, is again divided into S11-S18 first, to included in S1
The corresponding integral curve of initial seed point is compressed, and compressed byte number is denoted as T1, is compressed respectively to S11-S18,
8 compressed byte numbers are added to obtain T2, if T1/T2≤compression ratio threshold value T, decision space S1 does not need dividing
At S11-S18.
Step S300: calling directly required integral curve in storage equipment, complete flow field data analysis or can
It is applied depending on changing.
Required for the flow field data being analyzed or are visualized in application, calling directly in the storage device
The integral curve of initial seed point at corresponding grid point locations.
In present embodiment, corresponding mesh point corresponding position is first required for calling directly in the storage device
The concrete mode of the integral curve of beginning seed point are as follows:
1) it determines and is analyzed or visualized using specified flow field data space region, and in the area of space
The middle space density for placing initial seed point；Space density, which refers to, places an initial seed point every several mesh points；
2) subspace intersected with the area of space is searched in all subspaces of spatial mesh structure；Two spaces
Region intersection just refers to that two spaces have the part of intersection, is exactly the part that plane domain has intersection for two-dimensional grid section,
For three-dimensional grid section, just refer to that cuboid has the part of intersection；
3) product that preservation is associated in storage equipment is searched according to the index of the subspace intersected with the area of space
Component curve；
4) corresponding integral curve is extracted in the integral curve found according to the space density.
Certainly, if integral curve is stored in storage equipment after compressing, the index first according to subspace is needed
Corresponding compressed data is found in storage equipment, extracts corresponding integral curve again after then decompressing.
Wherein, convection current field data visualization application include but is not limited to integral curve it is direct rendering, source converge inquire,
FTLE (finite time lyapunov index field) calculating and some other visualization and analysis method specially designed.
Using method provided by the present invention, the acquisition efficiency that can greatly improve flow field data integral curve can be in resource
Extensive flow field data are handled in limited situation.
It is corresponding with method shown in Fig. 1, a kind of multiplexing system of flow field integral curve is additionally provided in present embodiment
System, as shown in figure 3, the system includes that spatial mesh structure obtains module 100, integral curve computing module 200, integral curve guarantor
Storing module 300 and integral curve calling module 400.Wherein:
Spatial mesh structure obtains module 100, for obtaining the spatial mesh structure of flow field data；
Integral curve computing module 200, for placing initial seed on all mesh points of the spatial mesh structure
Point calculates the integral curve of each seed point；The integral curve is trace or streamline；
Integral curve preserving module 300, for the integral curve to be stored in storage equipment；
Integral curve calling module 400, for the flow field data to be analyzed or visualized in application, depositing described
The integral curve of initial seed point at corresponding grid point locations required for being called directly in storage equipment.
In present embodiment, the integral curve preserving module 300 may include space lattice division unit 310 and curve
Associated storage unit 320.
Space lattice division unit 310 is divided into more for dividing the area of space of the spatial mesh structure
Sub-spaces；
Curve associated storage unit 320 will be first included in every sub-spaces for establishing the index of every sub-spaces
The corresponding integral curve of beginning seed point is associated with the index of its subspace to be stored in storage equipment.
When the spatial mesh structure of flow field data is two-dimensional regular grid, the space lattice division unit 310 can be with
Including the first grid dividing subelement 311, the curve associated storage unit 320 may include the first associated storage subelement
321.When the spatial mesh structure is three-dimensional regular grid, the space lattice division unit 310 may include the second net
Lattice divide subelement 312, and the curve associated storage unit 320 may include the second associated storage subelement 322.Such as Fig. 4 institute
Show.
First grid dividing subelement 311, for the root node using the whole region of two-dimensional regular grid as quaternary tree,
The space of the two-dimensional regular grid is divided into multiple two-dimensional sub-spaces by the quad-tree structure for establishing two-dimensional space grid, and four
Pitch the corresponding two-dimensional sub-spaces of each leaf node of tree construction；
First associated storage subelement 321, the index of each leaf node for establishing quad-tree structure, will be each
The corresponding integral curve of initial seed point that the corresponding two-dimensional sub-spaces of a leaf node are included is associated with guarantor with its node index
It is stored in storage equipment.
Second grid dividing subelement 312, for the root using the whole region of the three-dimensional regular grid as Octree
Node establishes the octree structure of three-dimensional regular grid, and it is empty that the space of the three-dimensional regular grid is divided into multiple three-dimensional sons
Between, corresponding three n-dimensional subspace ns of each leaf node of octree structure；
Second associated storage subelement 322 will be each for establishing the index of each leaf node in octree structure
The corresponding integral curve of initial seed point included in three n-dimensional subspace ns corresponding to a leaf node and its node index close
UNPROFOR exists in storage equipment.
The integral curve preserving module can also include integral curve compression unit 330.
Integral curve compression unit 330 is set for will be stored in storage after the corresponding integral curve compression of initial seed point
In standby.
In present embodiment, the space lattice division unit 310 includes first division subelement 3101, divides son again
Unit 3102 and division result determine subelement 3103, as shown in Figure 5.
It is first to divide subelement 3101, for being divided the spatial mesh structure for the first time, it is empty to obtain several height
Between；
Subelement 3102 is divided again, and every sub-spaces for dividing first division subelement are again divided into
M sub-spaces；
Division result determines subelement 3103, the subspace for dividing again for each, and judging will be before division
M after the corresponding integral curve of initial seed point included in subspace carries out compressed byte number N1 and will divide is empty
Between included in the corresponding integral curve of initial seed point carry out the sum of compressed M compressed byte number N2's respectively
Whether ratio is less than or equal to compression ratio threshold value, if so, the space of the two-dimensional regular grid or three-dimensional regular grid divides
It completes, the division result before dividing again is final division result, divides subelement 3102 again if it is not, then entering.
The integral curve calling module 400 may include area determination unit 410 to be analyzed, subspace determination unit
420, integral curve searching unit 430 and integral curve extraction unit 440.
Area determination unit 410 to be analyzed is analyzed or is visualized using specified flow field data sky for determining
Between region, and in the area of space place initial seed point space density；
Subspace determination unit 420 intersects for searching in the subspace of spatial mesh structure with the area of space
Subspace；
Integral curve searching unit 430, for being searched according to the index of the subspace intersected with the area of space
The integral curve of preservation is associated in storage equipment；
Integral curve extraction unit 440, for being extracted in the integral curve found according to the space density pair
The integral curve answered.
For a better understanding of the present invention, invention is further explained combined with specific embodiments below.
Embodiment
Flow field data in the present embodiment are hurricane Isabel data, the spatial mesh structure of hurricane Isabel data
It is defined on 500 × 500 × 100 three-dimensional space grid, as shown in Figure 6.
Firstly, initial seed point is placed on all mesh points of three-dimensional space grid, such as the dot that grid is pointed out in Fig. 6
It is shown, allow these seed points to move in the velocity field of hurricane Isabel data, and the integral for calculating each seed point is bent
Line.
Later, using entire three-dimensional spatial area as the root node of Octree, by three-dimensional space grid according to Octree knot
Structure is divided into multiple three n-dimensional subspace ns, and the index of each three n-dimensional subspace n is arranged, each leaf node pair of octree structure
Answer three n-dimensional subspace ns (can simply regard a cuboid as)
The corresponding integral curve of initial seed point included in every sub-spaces is compressed, compressed data with
The index of the corresponding leaf node in subspace is associated with.The integral curve of about 160GB is compressed to 8GB in the present embodiment, and will be each
The location error (compressed error) of seed point controls within 1 lattice point on integral curve.
User can star repeatedly for the analysis of the flow field data or visualization application, according to octree structure later
In corresponding leaf node index, from compressed data decompression extract needed for integral curve.For example, answering
With needing to calculate FTLE, need to place an initial seed points every 4 lattice points in three directions, it is total up to 125 ×
125 × 25 seed points, at this time, it is only necessary to which the seed point of corresponding position is extracted in decompression from stored data
Integral curve, which is brought, carries out subsequent FTLE of calculating, substantially increases treatment effeciency.
Obviously, various changes and modifications can be made to the invention without departing from essence of the invention by those skilled in the art
Mind and range.In this way, if these modifications and changes of the present invention belongs to the range of the claims in the present invention and its equivalent technology
Within, then the present invention is also intended to include these modifications and variations.
Claims (9)
1. a kind of multiplexing method of flow field integral curve, comprising the following steps:
(1) spatial mesh structure of flow field data is obtained；
(2) initial seed point is placed on all mesh points of the spatial mesh structure, the integral for calculating each seed point is bent
Line, and will be stored in storage equipment after the corresponding integral curve compression of initial seed point；The integral curve is trace or stream
Line；
Wherein, the area of space of the spatial mesh structure is divided, is divided into multiple subspaces, comprising: a. setting product
The compression ratio threshold value of component curve；
B. the spatial mesh structure is divided for the first time, obtains several sub-spaces；
C. every sub-spaces are again divided into M sub-spaces；
D. for the subspace that each is divided again, by the corresponding product of initial seed point included in the subspace before division
Component curve is compressed, and compressed byte number is denoted as N1, and initial seed point included in the subspace M after division is corresponding
Integral curve compressed respectively, the sum of M compression postbyte number is denoted as N2, judges whether the ratio of N1 and N2 is less than or waits
In compression ratio threshold value, if so, the Spacial domain decomposition is completed, the division result before dividing again is final division result,
If it is not, then return step c；
(3) required for the flow field data being analyzed or are visualized in application, calling directly in the storage device
The integral curve of initial seed point at corresponding grid point locations.
2. a kind of multiplexing method of flow field integral curve according to claim 1, it is characterised in that: described in step (1)
Spatial mesh structure is regular grid, when the spatial mesh structure of the flow field data of acquisition is irregular grid, irregular
On the basis of grid, the rule space network of flow field data is constructed；The regular grid is two-dimensional regular grid or three-dimensional
Regular grid.
3. a kind of multiplexing method of flow field integral curve according to claim 2, it is characterised in that: in step (2), by institute
It states integral curve and is stored in storage equipment, comprising:
The index for establishing every sub-spaces, by the corresponding integral curve of initial seed point included in every sub-spaces and its son
The index association in space is stored in storage equipment.
4. a kind of multiplexing method of flow field integral curve according to claim 3, it is characterised in that: when the space lattice
When structure is two-dimensional regular grid, the integral curve is stored in storage equipment, comprising:
1) quaternary tree of two-dimensional space grid is established as the root node of quaternary tree using the whole region of the two-dimensional regular grid
The space of the two-dimensional regular grid is divided into multiple two-dimensional sub-spaces, each leaf node pair of quad-tree structure by structure
Answer a two-dimensional sub-spaces；
2) index for establishing each leaf node of quad-tree structure, by the corresponding two-dimensional sub-spaces institute of each leaf node
The corresponding integral curve of the initial seed point for including is associated with its node index and is stored in storage equipment；
When the spatial mesh structure is three-dimensional regular grid, the integral curve is stored in storage equipment, comprising:
1. establishing the Octree of three-dimensional regular grid as the root node of Octree using the whole region of the three-dimensional regular grid
The space of the three-dimensional regular grid is divided into multiple three n-dimensional subspace ns, each leaf node of octree structure by structure
Corresponding three n-dimensional subspace ns；
2. the index of each leaf node in octree structure is established, by three n-dimensional subspace ns corresponding to each leaf node
Included in the corresponding integral curve of initial seed point and its node index be associated with to be stored in and stores in equipment.
5. a kind of multiplexing method of flow field integral curve according to one of claim 3 to 4, it is characterised in that: step (3)
In, the integral of the initial seed point of corresponding mesh point corresponding position required for calling directly in the storage device is bent
Line, comprising:
3.1) it determines and is analyzed or visualized using specified flow field data space region, and in the area of space
Place the space density of initial seed point；
3.2) subspace intersected with the area of space is searched in all subspaces of spatial mesh structure；
3.3) integral that preservation is associated in storage equipment is searched according to the index of the subspace intersected with the area of space
Curve；
3.4) corresponding integral curve is extracted in the integral curve found according to the space density.
6. a kind of multiplex system of flow field integral curve, comprising:
Spatial mesh structure obtains module, and for obtaining the spatial mesh structure of flow field data, the spatial mesh structure is rule
Then grid, the regular grid are two-dimensional regular grid or three-dimensional regular grid；
Integral curve computing module is calculated for placing initial seed point on all mesh points of the spatial mesh structure
The integral curve of each seed point；The integral curve is trace or streamline；
It is empty to be divided into multiple sons for dividing the area of space of the spatial mesh structure for space lattice division unit
Between, comprising:
First division subelement obtains several sub-spaces for being divided the spatial mesh structure for the first time；
Subelement is divided again, and it is empty that every sub-spaces for dividing first division subelement are again divided into M son
Between；
Division result determines subelement, the subspace for dividing again for each, and judging will be in the subspace before division
The corresponding integral curve of the initial seed point for being included carries out compressed byte number N1 and will be wrapped in the subspace M after dividing
The corresponding integral curve of the initial seed point contained carry out respectively the sum of compressed M compressed byte number N2 ratio whether
Less than or equal to compression ratio threshold value, if so, the two-dimensional regular grid or three-dimensional regular grid dividing are completed, before dividing again
Division result be final division result, if it is not, then enter divide subelement again；
Integral curve preserving module, including integral curve compression unit, for compressing the corresponding integral curve of initial seed point
After be stored in storage equipment in；
Integral curve calling module, for the flow field data being analyzed or being visualized in application, in the storage equipment
In call directly required for initial seed point at corresponding grid point locations integral curve.
7. a kind of multiplex system of flow field integral curve according to claim 6, it is characterised in that: the integral curve is protected
Storing module includes:
Curve associated storage unit, for establishing the index of every sub-spaces, by initial seed included in every sub-spaces
The corresponding integral curve of point is associated with the index of its subspace to be stored in storage equipment.
8. a kind of multiplex system of flow field integral curve according to claim 7, it is characterised in that: when the space lattice
When structure is two-dimensional regular grid, the space lattice division unit, comprising:
First grid dividing subelement establishes two for the root node using the whole region of two-dimensional regular grid as quaternary tree
The space of the two-dimensional regular grid is divided into multiple two-dimensional sub-spaces, quaternary tree knot by the quad-tree structure of dimension space grid
The corresponding two-dimensional sub-spaces of each leaf node of structure；
The curve associated storage unit includes:
First associated storage subelement, the index of each leaf node for establishing quad-tree structure, by each leaf
The corresponding integral curve of initial seed point that the corresponding two-dimensional sub-spaces of node are included is associated with to be stored in and deposit with its node index
It stores up in equipment；
When the spatial mesh structure is three-dimensional regular grid, the space lattice division unit, comprising:
Second grid dividing subelement, for building using the whole region of the three-dimensional regular grid as the root node of Octree
The space of the three-dimensional regular grid is divided into multiple three n-dimensional subspace ns, eight forks by the octree structure of vertical three-dimensional regular grid
Corresponding three n-dimensional subspace ns of each leaf node of tree construction；
The curve associated storage unit includes:
Second associated storage subelement, for establishing the index of each leaf node in octree structure, by each leaf
The corresponding integral curve of initial seed point included in three n-dimensional subspace ns corresponding to node is associated with preservation with its node index
In storage equipment.
9. the multiplex system of any one flow field integral curve according to claim 6-8, it is characterised in that: the integral
Curve calling module includes:
Area determination unit to be analyzed is analyzed or is visualized for determining using specified flow field data space region,
And the space density of initial seed point is placed in the area of space；
Subspace determination unit, it is empty for searching the son intersected with the area of space in the subspace of spatial mesh structure
Between；
Integral curve searching unit, for searching storage equipment according to the index of the subspace intersected with the area of space
The integral curve that middle association saves；
Integral curve extraction unit, for extracting corresponding integral in the integral curve found according to the space density
Curve.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610584306.5A CN106202522B (en) | 2016-07-22 | 2016-07-22 | A kind of multiplexing method and system of flow field integral curve |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610584306.5A CN106202522B (en) | 2016-07-22 | 2016-07-22 | A kind of multiplexing method and system of flow field integral curve |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106202522A CN106202522A (en) | 2016-12-07 |
| CN106202522B true CN106202522B (en) | 2019-11-08 |
Family
ID=57492295
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610584306.5A Active CN106202522B (en) | 2016-07-22 | 2016-07-22 | A kind of multiplexing method and system of flow field integral curve |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106202522B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109117433B (en) * | 2017-06-23 | 2022-05-24 | 菜鸟智能物流控股有限公司 | Index tree object creation and index method and related device thereof |
| CN110941894B (en) * | 2019-11-05 | 2024-02-06 | 天津大学 | In-situ visualization method for large-scale fluid simulation based on supercomputer |
| CN112948643B (en) * | 2021-05-13 | 2021-08-06 | 中国空气动力研究与发展中心计算空气动力研究所 | Structured grid streamline integration method based on thread parallelism |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102521854A (en) * | 2011-12-29 | 2012-06-27 | 北京理工大学 | Parallel flow line placing method applicable to two-dimensional flow field |
| CN102646288A (en) * | 2012-02-23 | 2012-08-22 | 北京理工大学 | Three-dimension flow line placing method with basically uniform intervals |
| CN104658027A (en) * | 2015-02-11 | 2015-05-27 | 中国海洋大学 | Three-dimensional streamline dynamic visualization algorithm facing irregular ocean flow field data |
| CN105653881A (en) * | 2016-03-04 | 2016-06-08 | 天津大学 | Flow field visualization method based on multi-density levels |
-
2016
- 2016-07-22 CN CN201610584306.5A patent/CN106202522B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102521854A (en) * | 2011-12-29 | 2012-06-27 | 北京理工大学 | Parallel flow line placing method applicable to two-dimensional flow field |
| CN102646288A (en) * | 2012-02-23 | 2012-08-22 | 北京理工大学 | Three-dimension flow line placing method with basically uniform intervals |
| CN104658027A (en) * | 2015-02-11 | 2015-05-27 | 中国海洋大学 | Three-dimensional streamline dynamic visualization algorithm facing irregular ocean flow field data |
| CN105653881A (en) * | 2016-03-04 | 2016-06-08 | 天津大学 | Flow field visualization method based on multi-density levels |
Non-Patent Citations (4)
| Title |
|---|
| SECONDO:An Extensible DBMS Architecture and Prototype;Ralf Hartmut Guting et al;《Collaborative Design》;20041231;1-30 * |
| 基于MATLAB的振动压路机计算机仿真;何建和等;《计算机仿真》;20040630;第21卷(第6期);第3.3节 * |
| 复杂向量场数据可视化技术研究与实现;周璐;《中国博士学位论文全文数据库 信息科技辑》;20040315;第2004年卷(第1期);第1.3、2.2-2.3、2.5、7.3节 * |
| 轨迹数据可视分析研究;王祖超等;《计算机辅助设计与图形学学报》;20150131;第27卷(第1期);9-25 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106202522A (en) | 2016-12-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105446979B (en) | Data digging method and node | |
| CN108205703B (en) | Multi-input multi-output matrix average value pooling vectorization implementation method | |
| CN106202522B (en) | A kind of multiplexing method and system of flow field integral curve | |
| CN107967516A (en) | A kind of acceleration of neutral net based on trace norm constraint and compression method | |
| CN105739951B (en) | A kind of L1 minimization problem fast solution methods based on GPU | |
| CN111612143A (en) | Compression method and system of deep convolutional neural network | |
| CN110874636B (en) | Neural network model compression method and device and computer equipment | |
| CN109299298A (en) | Construction method, device, application method and the system of image fusion model | |
| CN108304926B (en) | Pooling computing device and method suitable for neural network | |
| CN111523713A (en) | Method and device for predicting residual oil saturation distribution in oil field | |
| CN111008631B (en) | Image association method and device, storage medium and electronic device | |
| CN110796251A (en) | Image compression optimization method based on convolutional neural network | |
| CN109117742A (en) | Gestures detection model treatment method, apparatus, equipment and storage medium | |
| CN107644063A (en) | Time series analysis method and system based on data parallel | |
| CN112052940B (en) | Social network feature dynamic extraction method based on vector compression and reconstruction | |
| CN104951442B (en) | A kind of method and apparatus of definitive result vector | |
| CN113256793A (en) | Three-dimensional data processing method and system | |
| CN116167425B (en) | Neural network acceleration method, device, equipment and medium | |
| CN107506572A (en) | The method and apparatus for obtaining the height of target point | |
| CN107908696A (en) | A kind of parallel efficiently multidimensional space data clustering algorithm GRIDEN based on grid and density | |
| CN106990913A (en) | A kind of distributed approach of extensive streaming collective data | |
| CN108052755B (en) | Vector space based on completely random forest calculates intensity prediction method and system | |
| CN115983343A (en) | YOLOv4 convolutional neural network lightweight method based on FPGA | |
| CN113095279B (en) | Intelligent visual recognition method, device and system for flower quantity of fruit tree and storage medium | |
| CN112149826B (en) | Profile graph-based optimization method in deep neural network inference calculation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |