CN106202522A - The multiplexing method of a kind of flow field integral curve and system - Google Patents
The multiplexing method of a kind of flow field integral curve and system Download PDFInfo
- Publication number
- CN106202522A CN106202522A CN201610584306.5A CN201610584306A CN106202522A CN 106202522 A CN106202522 A CN 106202522A CN 201610584306 A CN201610584306 A CN 201610584306A CN 106202522 A CN106202522 A CN 106202522A
- Authority
- CN
- China
- Prior art keywords
- integral curve
- flow field
- curve
- space
- 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.)
- Granted
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 multiplexing method and the system of a kind of flow field integral curve, belong to and relate to data analysis technique field, flow field.The method includes: obtain the spatial mesh structure of flow field data；All mesh points of spatial mesh structure are placed initial seed point, calculates the integral curve of each seed points, and described integral curve is saved in storage device；Described integral curve is trace or streamline；When described flow field data are analyzed or visualize application, described storage device directly invokes the integral curve of initial seed point seed points at required corresponding grid point locations.The method and system are by precalculating abundant integral curve and preserving, in flow field data analysis later or visualization application, the regional space directly specified according to flow field analysis or visualization application, required integral curve can be rapidly obtained, drastically increase treatment effeciency, be particularly suited for analysis and the visualization processing of extensive flow field data.
Description
Technical field
The present invention relates to data analysis technique field, flow field, be specifically related to the multiplexing method of a kind of flow field integral curve and be
System.
Background technology
Space occupied by the motion of flow field (Flow Field) i.e. fluid, flow field data are the science data that a class is important,
It is typically defined in two dimension or three dimensions, and including at least a velocity field.Modal flow field data include ocean
Ocean current data in simulation, the wind field data in atmospheric simulation etc..Many existing flow field visualization methods are all based on integration
Directly the rendering of curve, such as integral curve, source converge inquiry, FTLE field (finite time lyapunov index field) calculate, with
And some other custom-designed visualization and the method for analysis.In Unsteady Flow (i.e. velocity field changes over), the most often
The integral curve seen be trace (in flow field certain position place a particle without quality (referred to as initial seed point), allow its
Moving in velocity field, the track formed is exactly a trace)；In steady flow field (i.e. velocity field does not changes over time), then
Correspond to streamline (curve all tangent with velocity on every bit in flow field is referred to as streamline).In Practical Calculation, trace and
Streamline is all to be tried to achieve by numerical integration, so they are called integral curve.
Common flow field visualization method is generally according to needing to place seed points in different time and position, and calculates big
The trace of amount.For bigger flow field data, this calculating process devotes a tremendous amount of time cost, and owing to calculating gained
Trace data amount is the biggest, it will usually abandoned by the trace as intermediate object program, the most more causes the waste calculating resource.In order to
The deficiencies in the prior art, the present invention are overcome to provide multiplexing method and the system of a kind of flow field integral curve.
Summary of the invention
For defect present in prior art, it is an object of the invention to provide the multiplexing side of a kind of flow field integral curve
Method and system, by the method and system, it is possible to when being effectively improved flow field data analysis or visualization application, integral curve obtains
Take efficiency.
For achieving the above object, the technical solution used in the present invention is: the multiplexing method of a kind of flow field integral curve, including
Following steps:
(1) spatial mesh structure of flow field data is obtained；
(2) on all mesh points of described spatial mesh structure, place initial seed point, calculate the long-pending of each seed points
Component curve, and described integral curve is saved in storage device；Described integral curve is trace or streamline；
(3), when described flow field data being analyzed or visualize application, directly invoke required in described storage device
The integral curve of the initial seed point at corresponding grid point locations wanted.
Further, the multiplexing method of a kind of flow field as above integral curve, in step (1), described space lattice is tied
Structure is regular grid, when the spatial mesh structure of the flow field data obtained is irregular grid, on the basis of irregular grid
On, build the rule space network of flow field data；Described regular grid is two-dimensional regular grid or three-dimensional regular grid.
Further, the multiplexing method of a kind of flow field as above integral curve, in step (2), by described integral curve
It is saved in storage device, including:
2.1) area of space of described spatial mesh structure is divided, be divided into many sub spaces；
2.2) index of every sub spaces is set up, integration corresponding for the initial seed point included in every sub spaces is bent
Line associates with the index of its subspace and is stored in storage device.
Further, the multiplexing method of a kind of flow field as above integral curve, when described spatial mesh structure is two dimension
During regular grid, described integral curve is saved in storage device, including:
1) using the whole region of described two-dimensional regular grid as the root node of quaternary tree, the four of two-dimensional space grid are set up
Fork tree construction, is divided into multiple two-dimensional sub-spaces by the space of described two-dimensional regular grid, each leaf joint of quad-tree structure
The corresponding two-dimensional sub-spaces of point；
2) index of each leaf node of quad-tree structure is set up, by empty for two dimension corresponding for each leaf node
Between integral curve corresponding to the initial seed point that comprised associate with its node index and be stored in storage device；
When described spatial mesh structure is three-dimensional regular grid, described integral curve is saved in storage device, including:
1. using the whole region of described three-dimensional regular grid as the root node of Octree, the eight of three-dimensional regular grid are set up
Fork tree construction, is divided into multiple three n-dimensional subspace ns by the space of described three-dimensional regular grid, each leaf of octree structure
Corresponding three n-dimensional subspace ns of node；
2. the index of each leaf node in octree structure is set up, by the three-dimensional son corresponding to each leaf node
Integral curve corresponding to initial seed point included in space associates with its node index and is saved in storage device.
Further, the multiplexing method of a kind of flow field as above integral curve, in step (2), by initial seed point pair
It is saved in storage device after the integral curve compression answered.
Further, the multiplexing method of a kind of flow field as above integral curve, by the space of described spatial mesh structure
Region is divided into many sub spaces, including:
A., the compression ratio threshold value of integral curve is set；
B. described spatial mesh structure is divided for the first time, obtain several subspaces；
C. every sub spaces is again divided into M sub spaces；
D. the subspace again divided for each, the initial seed point included in subspace before dividing is corresponding
Integral curve be compressed, the byte number after compression is designated as N1, will divide after M subspace included in initial seed point
Corresponding integral curve is compressed respectively, and M compression postbyte number sum is designated as N2, it is judged that whether the ratio of N1 with N2 is less than
Or equal to compression ratio threshold value, the space of the most described two-dimensional regular grid or three-dimensional regular grid has divided, and again divides
Front division result is final division result, if it is not, then return step c.
Further, the multiplexing method of a kind of flow field as above integral curve, in step (3), in described storage device
In directly invoke the integral curve of initial seed point of required corresponding mesh point corresponding position, including:
3.1) the data space region, flow field being analyzed or visualize specified by application is determined, and in described space region
Territory is placed the spatial density of initial seed point；
3.2) in all subspaces of spatial mesh structure, search the subspace intersected with described area of space；
3.3) association preservation in the index search storage device of the subspace intersected according to described and described area of space
Integral curve；
3.4) in the integral curve found, the integral curve of correspondence is extracted according to described spatial density.
A kind of multiplex system of flow field integral curve, including:
Spatial mesh structure acquisition 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 described spatial mesh structure,
Calculate the integral curve of each seed points；Described integral curve is trace or streamline；
Integral curve preserves module, for being saved in storage device by described integral curve；
Integral curve calling module, in time being analyzed described flow field data or visualize application, in described storage
Equipment directly invokes the integral curve of initial seed point at required corresponding grid point locations.
Further, the multiplex system of a kind of flow field as above integral curve, described integral curve preserves module and includes:
Space lattice division unit, for being divided by the area of space of described spatial mesh structure, is divided into multiple
Subspace；
Curve association memory element, for setting up the index of every sub spaces, initial by included in every sub spaces
Integral curve corresponding to seed points associates with the index of its subspace and is stored in storage device.
Further, the multiplex system of a kind of flow field as above integral curve, when described spatial mesh structure is two dimension
During regular grid, described space lattice division unit, including:
First stress and strain model subelement, for the whole region using two-dimensional regular grid as the root node of quaternary tree, builds
The quad-tree structure of vertical two-dimensional space grid, is divided into multiple two-dimensional sub-spaces, four forks by the space of described two-dimensional regular grid
The corresponding two-dimensional sub-spaces of each leaf node of tree construction；
Described curve association memory element includes:
First association storing sub-units, for setting up the index of each leaf node of quad-tree structure, by each
Integral curve corresponding to initial seed point that the two-dimensional sub-spaces that leaf node is corresponding is comprised associates preservation with its node index
In storage device；
When described spatial mesh structure is three-dimensional regular grid, described space lattice division unit, including:
Second stress and strain model subelement, saves for the root using the whole region of described three-dimensional regular grid as Octree
Point, sets up the octree structure of three-dimensional regular grid, and the space of described 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；
Described curve association memory element includes:
Second association storing sub-units, for setting up the index of each leaf node in octree structure, by each
Integral curve corresponding to the initial seed point included in three n-dimensional subspace ns corresponding to leaf node associates with its node index
It is saved in storage device.
Further, the multiplex system of a kind of flow field as above integral curve, described integral curve preserves module and also wraps
Include:
Integral curve compression unit, for being saved in storage device after integral curve compression corresponding for initial seed point
In.
Further, the multiplex system of a kind of flow field as above integral curve, described space lattice division unit includes:
Divide subelement for the first time, for being divided for the first time by described spatial mesh structure, obtain several subspaces；
Again divide subelement, for dividing every sub spaces of obtaining be again divided into M by dividing subelement for the first time
Subspace；
Division result determines subelement, for the subspace again divided for each, it is judged that the son before dividing is empty
Integral curve corresponding to initial seed point included between be compressed after byte number N1 with will divide after M subspace in
Integral curve corresponding to the initial seed point that comprised be compressed respectively after M compression after the ratio of byte number sum N2
Whether less than or equal to compression ratio threshold value, the space of the most described two-dimensional regular grid or three-dimensional regular grid has divided,
Division result before again dividing is final division result, again divides subelement if it is not, then enter.
Further, the multiplex system of a kind of flow field as above integral curve, described integral curve calling module includes:
Area determination unit to be analyzed, for determining the data space district, flow field being analyzed or visualize specified by application
Territory, and in described area of space, place the spatial density of initial seed point；
Subspace determines unit, for searching the son intersected with described area of space in the subspace of spatial mesh structure
Space；
Integral curve searches unit, the index search storage of the subspace for intersecting according to described and described area of space
Equipment associates the integral curve of preservation；
Integral curve extraction unit, for extracting correspondence according to described spatial density in the integral curve found
Integral curve.
The beneficial effects of the present invention is: the multiplexing method of flow field provided by the present invention integral curve and system, pass through
Precalculate abundant integral curve and preserve, in flow field data analysis later or visualization application, direct root
The regional space specified according to flow field analysis or visualization application, corresponding seed points in corresponding subspace can be rapidly obtained
Integral curve, the most less acquisition time of integral curve, improves treatment effeciency, is particularly suited for extensive flow field number
According to analysis and visualization processing.
Accompanying drawing explanation
Fig. 1 is the flow chart of the multiplexing method of a kind of flow field integral curve of offer in detailed description of the invention；
Fig. 2 is the flow chart of the spatial mesh structure division of flow field data in detailed description of the invention；
Fig. 3 is the structured flowchart of the multiplex system of a kind of flow field integral curve of offer in detailed description of the invention；
Fig. 4 is the structured flowchart that in detailed description of the invention, integral curve preserves module；
Fig. 5 is the structured flowchart of space stress and strain model unit in detailed description of the invention；
Fig. 6 is the placement schematic diagram of three-dimensional regular network and initial seed point in embodiment.
Detailed description of the invention
Below in conjunction with Figure of description, the present invention is described in further detail with detailed description of the invention.
The multiplexing method of flow field provided by the present invention integral curve and system, basic ideas are to precalculate follow-up needs
The abundant integral curve of flow field data analyzed also stores, in analysis later or visualization, directly from being deposited
Extracting data needed for integral curve.
Fig. 1 shows the flow chart of the multiplexing method of a kind of flow field integral curve provided in present embodiment, by figure
It can be seen that the method mainly can include three below step:
Step S100: obtain the spatial mesh structure of flow field data；
Follow-up needs are analyzed or are visualized the flow field data of application, first obtain its spatial mesh structure.
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, described spatial mesh structure is regular grid, when the space lattice knot of the flow field data obtained
When structure is irregular grid, on the basis of irregular grid, build the rule space network of flow field data；Described rule
Grid is two-dimensional regular grid or three-dimensional regular grid.
If the most regular network of flow field data, can directly use its regular grid structure,
If the network of flow field data itself is irregular grid, as long as the space at the place of its data own is three dimensions, just
Can build the regular grid of a set of flow field data, concrete building mode uses the space lattice building mode of existing flow field data
?.
Step S200: place initial seed point on all mesh points of the spatial mesh structure of flow field data, calculates every
The integral curve of individual seed points also preserves；
All mesh points of the spatial mesh structure of flow field data are placed initial seed point, calculates each seed points
Integral curve, and described integral curve is saved in storage device；Described integral curve is trace or streamline.
Described integral curve is trace or the streamline of initial seed point, Unsteady Flow (i.e. velocity field changes over)
In, modal integral curve is trace, in steady flow field (i.e. velocity field does not changes over time), then corresponds to streamline.Institute
Stating initial seed point, i.e. in flow field, a particle without quality is placed in certain position, allows it move in velocity field, is formed
Track is exactly a trace or streamline.
In present embodiment, the integral curve of the initial seed point placed thick and fast in space defined in the data of flow field,
Track should be the longest, until seed points is walked out defined space and stopped.The computational methods of integral curve use existing trace
Or the calculation of streamline or use existing visualization tool to directly generate, in actual applications, when flow field data relatively
Hour, can calculate in single computer, when the data is larger, can be in parallel computation under cluster environment.
In present embodiment, the mode that the integral curve of all initial seed point is saved in storage device is:
1) area of space of described spatial mesh structure is divided, be divided into many sub spaces；
2) index of every sub spaces is set up, by integral curve corresponding for the initial seed point included in every sub spaces
Associate with the index of its subspace and be stored in storage device.
By this storage mode, when follow-up lookup integral curve, it is only necessary to the index according to subspace can be in storage
Equipment finds the joint component curve of correspondence.
For different spatial mesh structures, present embodiments provide for two kinds of different integral curve preserving types.When
When described spatial mesh structure is two-dimensional regular grid, the concrete mode that described integral curve is saved in storage device is:
1) using the whole region of described two-dimensional regular grid as the root node of quaternary tree, the four of two-dimensional space grid are set up
Fork tree construction, is divided into multiple two-dimensional sub-spaces by the space of described two-dimensional regular grid, each leaf joint of quad-tree structure
The corresponding two-dimensional sub-spaces of point；
2) index of each leaf node of quad-tree structure is set up, by empty for two dimension corresponding for each leaf node
Between integral curve corresponding to the initial seed point that comprised associate with its node index and be stored in storage device.
When described spatial mesh structure is three-dimensional regular grid, described integral curve is saved in the concrete of storage device
Mode is:
1. using the whole region of described three-dimensional regular grid as the root node of Octree, the eight of three-dimensional regular grid are set up
Fork tree construction, is divided into multiple three n-dimensional subspace ns by the space of described three-dimensional regular grid, each leaf of octree structure
Corresponding three n-dimensional subspace ns of node；
2. the index of each leaf node in octree structure is set up, by the three-dimensional son corresponding to each leaf node
Integral curve corresponding to initial seed point included in space associates with its node index and is saved in storage device.
Owing to the integral curve of calculative initial seed point is typically hundreds of or thousands of times of primary flow field data, for
Guarantee all of flow field integral curve can store, and in present embodiment, sets integral curve is saved in storage
Time standby, first it is compressed the integral curve that initial seed point is corresponding processing, then the rope by the result after compression with subspace
Draw (or node index) association storage.
In actual applications, the compress mode of concrete integral curve can select according to practical situation.Such as, may be used
With and integral curve corresponding for the initial seed point included in every sub spaces is compressed together, it is also possible in subspace
Integral curve corresponding to each initial seed point of being comprised is compressed one by one, presses after the result after compression being merged again
Contracting.Concrete compression algorithm can also select as required, it is possible to use the compression that integral curve is fitted by parameter curve
Method, as Bezier fits to；Can also be the compression algorithm for any floating number array, such as fpzip algorithm, certainly
Can also be that multiple compression algorithm is applied in combination.In order to reduce owing to the initial seed point compressing the integral curve caused corresponding exists
Site error in grid, when march wire compression, it is generally required to arrange the compressed error upper limit.
For the division of the spatial mesh structure of flow field data, if the number of subspace is the fewest, institute in every sub spaces
The initial seed point comprised is the most, and corresponding integral curve is the most, the compression effectiveness the best (word after compression when being compressed
Joint number is less), but when the integral curve that subsequent calls needs, can be slower when of decompression, in order to balance these two aspects,
Additionally providing the dividing mode of a kind of spatial mesh structure in present embodiment, which is as in figure 2 it is shown, flow process is as follows:
A., the compression ratio threshold value of integral curve is set；Described compression ratio threshold value is empirical value.
B. described spatial mesh structure is divided for the first time, obtain several subspaces；
C. every sub spaces is again divided into M sub spaces；
D. the subspace again divided for each, the initial seed point included in subspace before dividing is corresponding
Integral curve be compressed, the byte number after compression is designated as N1, will divide after M subspace included in initial seed point
Corresponding integral curve is compressed respectively, and M compression postbyte number sum is designated as N2, it is judged that whether the ratio of N1 with N2 is less than
Or equal to compression ratio threshold value, the space of the most described two-dimensional regular grid or three-dimensional regular grid has divided, and again divides
Front division result is final division result, if it is not, then return step c.
Which is by corresponding long-pending of all initial seed point included in the space before judging a sub spaces and dividing
Byte number after component curve compression compresses respectively with by integral curve corresponding for all subspaces corresponding after this Subspace partition
After several byte number sums ratio and the compression ratio threshold value of setting between relation, determine whether to need at antithetical phrase empty
Between proceed divide.
Such as, when using octree structure to divide for above-mentioned three regular grids, first by three-dimensional grid really
Space is divided into 8 the first subspaces, is designated as space S 1-S8 respectively, judges whether the most by the way to need 8
One subspace divides again, such as, for space S 1, be first again divided into S11-S18, included in S1
Integral curve corresponding to initial seed point is compressed, and the byte number after compression is designated as T1, is compressed S11-S18 respectively,
Byte number after compressing 8 is added and obtains T2, if T1/T2≤compression ratio threshold value T, then decision space S1 need not dividing
Become S11-S18.
Step S300: directly invoke required integral curve in storage device, complete flow field data analysis or can
Depending on changing application.
When described flow field data are analyzed or visualize application, directly invoke required in described storage device
The integral curve of the initial seed point at corresponding grid point locations.
In present embodiment, directly invoke in described storage device at the beginning of required corresponding mesh point corresponding position
The concrete mode of the integral curve of beginning seed points is:
1) the data space region, flow field being analyzed or visualize specified by application is determined, and at described area of space
The spatial density of middle placement initial seed point；Spatial density refers to place an initial seed point every several mesh points；
2) in all subspaces of spatial mesh structure, search the subspace intersected with described area of space；Two spaces
Region is intersected and is referred to two spaces and have crossing part, interval for two-dimensional grid, it is simply that plane domain has crossing part,
Interval for three-dimensional grid, it is simply that to refer to that cuboid has crossing part；
3) the index search storage device of the subspace intersected according to described and described area of space associates the long-pending of preservation
Component curve；
4) in the integral curve found, the integral curve of correspondence is extracted according to described spatial density.
Certainly, if integral curve be compression after be saved in storage device, need first according to the index of subspace
In storage device, find the compression data of correspondence, after then decompressing, extract the integral curve of correspondence again.
Wherein, the visualization application of stream field data include but not limited to directly the rendering of integral curve, source converge inquiry,
FTLE field (finite time lyapunov index field) calculate and some other custom-designed visualization with analyze method.
Use method provided by the present invention, it is possible to the acquisition efficiency being greatly improved flow field data integral curve can be in resource
Extensive flow field data are processed in the case of limited.
Corresponding with the method shown in Fig. 1, present embodiment additionally provides the multiplexing system of a kind of flow field integral curve
System, as it is shown on figure 3, this system includes that spatial mesh structure acquisition module 100, integral curve computing module 200, integral curve are protected
Storing module 300 and integral curve calling module 400.Wherein:
Spatial mesh structure acquisition 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 described spatial mesh structure
Point, calculates the integral curve of each seed points；Described integral curve is trace or streamline；
Integral curve preserves module 300, for being saved in storage device by described integral curve；
Integral curve calling module 400, in time being analyzed described flow field data or visualize application, deposits described
Storage equipment directly invokes the integral curve of initial seed point at required corresponding grid point locations.
In present embodiment, described integral curve preserves module 300 can include space lattice division unit 310 and curve
Association memory element 320.
Space lattice division unit 310, for being divided by the area of space of described spatial mesh structure, is divided into many
Sub spaces；
Curve association memory element 320, for setting up the index of every sub spaces, at the beginning of included in every sub spaces
Integral curve corresponding to beginning seed points associates with the index of its subspace and is stored in storage device.
When the spatial mesh structure of flow field data is two-dimensional regular grid, described space lattice division unit 310 is permissible
Including the first stress and strain model subelement 311, described curve association memory element 320 can include the first association storing sub-units
321.When described spatial mesh structure is three-dimensional regular grid, described space lattice division unit 310 can include the second net
Lattice divide subelement 312, and described curve association memory element 320 can include the second association storing sub-units 322.Such as Fig. 4 institute
Show.
First stress and strain model subelement 311, for the root node using the whole region of two-dimensional regular grid as quaternary tree,
Set up the quad-tree structure of two-dimensional space grid, the space of described two-dimensional regular grid is divided into multiple two-dimensional sub-spaces, four
The corresponding two-dimensional sub-spaces of each leaf node of fork tree construction；
First association storing sub-units 321, for setting up the index of each leaf node of quad-tree structure, by each
Integral curve corresponding to initial seed point that the two-dimensional sub-spaces that individual leaf node is corresponding is comprised closes UNPROFOR with its node index
It is stored in storage device.
Second stress and strain model subelement 312, for the root using the whole region of described three-dimensional regular grid as Octree
Node, sets up the octree structure of three-dimensional regular grid, the space of described three-dimensional regular grid is divided into multiple three-dimensional son empty
Between, corresponding three n-dimensional subspace ns of each leaf node of octree structure；
Second association storing sub-units 322, for setting up the index of each leaf node in octree structure, by each
Integral curve corresponding to the initial seed point included in three n-dimensional subspace ns corresponding to individual leaf node closes with its node index
UNPROFOR exists in storage device.
Described integral curve preserves module can also include integral curve compression unit 330.
Integral curve compression unit 330, for setting being saved in storage after integral curve compression corresponding for initial seed point
In Bei.
In present embodiment, described space lattice division unit 310 includes dividing for the first time subelement 3101, again divides son
Unit 3102 and division result determine subelement 3103, as shown in Figure 5.
First division subelement 3101, for being divided for the first time by described spatial mesh structure, obtains some height empty
Between；
Again divide subelement 3102, for dividing every sub spaces of obtaining be again divided into dividing subelement for the first time
M sub spaces；
Division result determines subelement 3103, for the subspace again divided for each, it is judged that before dividing
Integral curve corresponding to initial seed point included in subspace be compressed after byte number N1 with will divide after M empty
Integral curve corresponding to initial seed point included between be compressed respectively after M compression after byte number sum N2
Whether ratio is less than or equal to compression ratio threshold value, and the space of the most described two-dimensional regular grid or three-dimensional regular grid divides
Completing, the division result before again dividing is final division result, if it is not, then enter division subelement 3102 again.
Described integral curve calling module 400 can include that area determination unit 410 to be analyzed, subspace determine unit
420, integral curve searches unit 430 and integral curve extraction unit 440.
Area determination unit 410 to be analyzed, for determining the flow field data sky being analyzed or visualize specified by application
Between region, and in described area of space, place the spatial density of initial seed point；
Subspace determines unit 420, intersects with described area of space for searching in the subspace of spatial mesh structure
Subspace；
Integral curve searches unit 430, the index search of the subspace for intersecting according to described and described area of space
Storage device associates the integral curve of preservation；
Integral curve extraction unit 440, right for extracting in the integral curve found according to described spatial density
The integral curve answered.
For a better understanding of the present invention, below in conjunction with specific embodiment, the present invention is further described.
Embodiment
Flow field data in the present embodiment are hurricane Isabel data, the spatial mesh structure of hurricane Isabel data
It is defined on the three-dimensional space grid of 500 × 500 × 100, as shown in Figure 6.
First, all mesh points of three-dimensional space grid are placed initial seed point, the round dot pointed out such as grid in Fig. 6
Shown in, allow these seed points move in the velocity field of hurricane Isabel data, and the integration calculating each seed points is bent
Line.
Afterwards, using whole three-dimensional spatial area as the root node of Octree, three-dimensional space grid is tied according to Octree
Structure is divided into multiple three n-dimensional subspace ns, and arranges the index of each three n-dimensional subspace ns, each leaf node pair of octree structure
Answer three n-dimensional subspace ns (can simply regard a cuboid as)
Integral curve corresponding for initial seed point included in every sub spaces is compressed, the data after compression with
The index association of the leaf node that subspace is corresponding.The integral curve of about 160GB is compressed to 8GB by the present embodiment, and will be each
On integral curve, the site error (compressed error) of seed points controls within 1 lattice point.
Afterwards user can start repeatedly for these flow field data analysis or visualization application, according to octree structure
The index of the leaf node of middle correspondence, decompresses the integral curve needed for extracting from the data after compression.Such as, having should
Calculating FTLE field with needs, it needs to place initial seed point every 4 lattice points in three directions, it is total up to 125 ×
125 × 25 seed points, now, it is only necessary to from the most stored data, the seed points of corresponding position is extracted in decompression
Integral curve brings the calculating carrying out follow-up FTLE field, substantially increases treatment effeciency.
Obviously, those skilled in the art can carry out various change and the modification essence without deviating from the present invention to the present invention
God and scope.So, if these amendments of the present invention and modification belong to the scope of the claims in the present invention and equivalent technology thereof
Within, then the present invention is also intended to comprise these change and modification.
Claims (13)
1. a multiplexing method for flow field integral curve, comprises the following steps:
(1) spatial mesh structure of flow field data is obtained；
(2) placing initial seed point on all mesh points of described spatial mesh structure, the integration calculating each seed points is bent
Line, and described integral curve is saved in storage device；Described integral curve is trace or streamline；
(3), when described flow field data being analyzed or visualize application, directly invoke required in described storage device
The integral curve of the initial seed point at corresponding grid point locations.
The multiplexing method of a kind of flow field the most according to claim 1 integral curve, it is characterised in that: in step (1), described
Spatial mesh structure is regular grid, when the spatial mesh structure of the flow field data obtained is irregular grid, irregular
On the basis of grid, build the rule space network of flow field data；Described regular grid is two-dimensional regular grid or three-dimensional
Regular grid.
The multiplexing method of a kind of flow field the most according to claim 2 integral curve, it is characterised in that: in step (2), by institute
State integral curve and be saved in storage device, including:
2.1) area of space of described spatial mesh structure is divided, be divided into many sub spaces；
2.2) set up the index of every sub spaces, by integral curve corresponding for the initial seed point included in every sub spaces with
The index association of its subspace is stored in storage device.
The multiplexing method of a kind of flow field the most according to claim 3 integral curve, it is characterised in that: when described space lattice
When structure is two-dimensional regular grid, described integral curve is saved in storage device, including:
1) using the whole region of described two-dimensional regular grid as the root node of quaternary tree, the quaternary tree of two-dimensional space grid is set up
Structure, is divided into multiple two-dimensional sub-spaces by the space of described two-dimensional regular grid, each leaf node pair of quad-tree structure
Answer a two-dimensional sub-spaces；
2) index of each leaf node of quad-tree structure is set up, by two-dimensional sub-spaces institute corresponding for each leaf node
Integral curve corresponding to the initial seed point that comprises associates with its node index and is stored in storage device；
When described spatial mesh structure is three-dimensional regular grid, described integral curve is saved in storage device, including:
1. using the whole region of described three-dimensional regular grid as the root node of Octree, the Octree of three-dimensional regular grid is set up
Structure, is divided into multiple three n-dimensional subspace ns by the space of described three-dimensional regular grid, each leaf node of octree structure
Corresponding three n-dimensional subspace ns；
2. the index of each leaf node in octree structure is set up, by three n-dimensional subspace ns corresponding to each leaf node
Included in integral curve corresponding to initial seed point associate with its node index and be saved in storage device.
The multiplexing method of a kind of flow field the most according to claim 3 integral curve, it is characterised in that: in step (2), will just
It is saved in storage device after the integral curve compression that beginning seed points is corresponding.
The multiplexing method of a kind of flow field the most according to claim 5 integral curve, it is characterised in that: by described space lattice
The Spacial domain decomposition of structure is many sub spaces, including:
A., the compression ratio threshold value of integral curve is set；
B. described spatial mesh structure is divided for the first time, obtain several subspaces；
C. every sub spaces is again divided into M sub spaces；
D. the subspace each divided again, corresponding the amassing of the initial seed point included in subspace before dividing
Component curve is compressed, and the byte number after compression is designated as N1, and the initial seed point included in M subspace after dividing is corresponding
Integral curve be compressed respectively, M compression postbyte number sum is designated as N2, it is judged that the ratio of N1 Yu N2 whether be less than or etc.
In compression ratio threshold value, the space of the most described two-dimensional regular grid or three-dimensional regular grid has divided, before again dividing
Division result is final division result, if it is not, then return step c.
7. according to the multiplexing method of a kind of flow field integral curve one of claim 3 to 6 Suo Shu, it is characterised in that: step (3)
In, the integration of the initial seed point directly invoking required corresponding mesh point corresponding position in described storage device is bent
Line, including:
3.1) the data space region, flow field being analyzed or visualize specified by application is determined, and in described area of space
Place the spatial density of initial seed point；
3.2) in all subspaces of spatial mesh structure, search the subspace intersected with described area of space；
3.3) the index search storage device of the subspace intersected according to described and described area of space associates the integration of preservation
Curve；
3.4) in the integral curve found, the integral curve of correspondence is extracted according to described spatial density.
8. a multiplex system for flow field integral curve, including:
Spatial mesh structure acquisition 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 described spatial mesh structure, calculates
The integral curve of each seed points；Described integral curve is trace or streamline；
Integral curve preserves module, for being saved in storage device by described integral curve；
Integral curve calling module, in time being analyzed described flow field data or visualize application, in described storage device
In directly invoke the integral curve of initial seed point at required corresponding grid point locations.
The multiplex system of a kind of flow field the most according to claim 8 integral curve, it is characterised in that: described integral curve is protected
Storing module includes:
Space lattice division unit, for being divided by the area of space of described spatial mesh structure, is divided into many height empty
Between；
Curve association memory element, for setting up the index of every sub spaces, by the initial seed included in every sub spaces
The integral curve that point is corresponding associates with the index of its subspace and is stored in storage device.
The multiplex system of a kind of flow field the most according to claim 9 integral curve, it is characterised in that: when described space networks
When lattice structure is two-dimensional regular grid, described space lattice division unit, including:
First stress and strain model subelement, for the whole region using two-dimensional regular grid as the root node of quaternary tree, sets up two
The quad-tree structure of dimension space grid, is divided into multiple two-dimensional sub-spaces by the space of described two-dimensional regular grid, and quaternary tree is tied
The corresponding two-dimensional sub-spaces of each leaf node of structure；
Described curve association memory element includes:
First association storing sub-units, for setting up the index of each leaf node of quad-tree structure, by each leaf
Integral curve corresponding to initial seed point that the two-dimensional sub-spaces that node is corresponding is comprised associates to be stored in its node index to be deposited
In storage equipment；
When described spatial mesh structure is three-dimensional regular grid, described space lattice division unit, including:
Second stress and strain model subelement, for the whole region using described three-dimensional regular grid as the root node of Octree, builds
The octree structure of vertical three-dimensional regular grid, is divided into multiple three n-dimensional subspace ns, eight forks by the space of described three-dimensional regular grid
Corresponding three n-dimensional subspace ns of each leaf node of tree construction；
Described curve association memory element includes:
Second association storing sub-units, for setting up the index of each leaf node in octree structure, by each leaf
Integral curve corresponding to the initial seed point included in three n-dimensional subspace ns corresponding to node associates preservation with its node index
In storage device.
The multiplex system of 11. a kind of flow field according to claim 9 integral curves, it is characterised in that: described integral curve
Preserve module also to include:
Integral curve compression unit, for being saved in storage device after integral curve compression corresponding for initial seed point.
The multiplex system of 12. a kind of flow field according to claim 11 integral curves, it is characterised in that: described space lattice
Division unit includes:
Divide subelement for the first time, for being divided for the first time by described spatial mesh structure, obtain several subspaces；
Again divide subelement, for dividing every sub spaces of obtaining be again divided into M son sky by dividing subelement for the first time
Between；
Division result determines subelement, for the subspace again divided for each, it is judged that in the subspace before dividing
Integral curve corresponding to the initial seed point that comprised be compressed after byte number N1 with will divide after M subspace in wrapped
Integral curve corresponding to the initial seed point that contains be compressed respectively after M compression after byte number sum N2 ratio whether
Less than or equal to compression ratio threshold value, the space of the most described two-dimensional regular grid or three-dimensional regular grid has divided, again
Division result before division is final division result, again divides subelement if it is not, then enter.
13. according to the multiplex system of a kind of flow field integral curve one of claim 9 to 12 Suo Shu, it is characterised in that: described
Integral curve calling module includes:
Area determination unit to be analyzed, for determining the data space region, flow field being analyzed or visualize specified by application,
And in described area of space, place the spatial density of initial seed point；
Subspace determines unit, empty for searching the son intersected with described area of space in the subspace of spatial mesh structure
Between；
Integral curve searches unit, the index search storage device of the subspace for intersecting according to described and described area of space
The integral curve that middle association preserves；
Integral curve extraction unit, for extracting the integration of correspondence in the integral curve found according to described spatial 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 true CN106202522A (en) | 2016-12-07 |
| CN106202522B 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) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109117433A (en) * | 2017-06-23 | 2019-01-01 | 菜鸟智能物流控股有限公司 | Index tree object creation method and index method and related device thereof |
| CN110941894A (en) * | 2019-11-05 | 2020-03-31 | 天津大学 | In-situ visualization method for large-scale fluid simulation based on supercomputer |
| CN112948643A (en) * | 2021-05-13 | 2021-06-11 | 中国空气动力研究与发展中心计算空气动力研究所 | 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 |
|---|
| RALF HARTMUT GUTING ET AL: "SECONDO:An Extensible DBMS Architecture and Prototype", 《COLLABORATIVE DESIGN》 * |
| 何建和等: "基于MATLAB的振动压路机计算机仿真", 《计算机仿真》 * |
| 周璐: "复杂向量场数据可视化技术研究与实现", 《中国博士学位论文全文数据库 信息科技辑》 * |
| 王祖超等: "轨迹数据可视分析研究", 《计算机辅助设计与图形学学报》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109117433A (en) * | 2017-06-23 | 2019-01-01 | 菜鸟智能物流控股有限公司 | Index tree object creation method and index method and related device thereof |
| CN110941894A (en) * | 2019-11-05 | 2020-03-31 | 天津大学 | In-situ visualization method for large-scale fluid simulation based on supercomputer |
| CN110941894B (en) * | 2019-11-05 | 2024-02-06 | 天津大学 | In-situ visualization method for large-scale fluid simulation based on supercomputer |
| CN112948643A (en) * | 2021-05-13 | 2021-06-11 | 中国空气动力研究与发展中心计算空气动力研究所 | Structured grid streamline integration method based on thread parallelism |
| CN112948643B (en) * | 2021-05-13 | 2021-08-06 | 中国空气动力研究与发展中心计算空气动力研究所 | Structured grid streamline integration method based on thread parallelism |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106202522B (en) | 2019-11-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105824974B (en) | The method and system of Data Analysis Services | |
| CN107967516A (en) | A kind of acceleration of neutral net based on trace norm constraint and compression method | |
| CN106202522A (en) | The multiplexing method of a kind of flow field integral curve and system | |
| CN109766949A (en) | Convolutional neural networks light weight method, device and electronic equipment | |
| CN107644063A (en) | Time series analysis method and system based on data parallel | |
| CN111553471A (en) | Data analysis processing method and device | |
| RU2019108264A (en) | SAVING AND EXTRACTING SPATIAL DATA TO / FROM DATABASE | |
| CN111160516B (en) | Convolutional layer sparsification method and device for deep neural network | |
| CN114820279B (en) | Distributed deep learning method and device based on multiple GPUs and electronic equipment | |
| CN108320019B (en) | Convolution calculation method and device for deep convolution neural network | |
| KR20150114950A (en) | Increasing frame rate of an image stream | |
| CN111860771A (en) | Convolutional neural network computing method applied to edge computing | |
| CN105404611A (en) | Matrix model based multi-calculation-engine automatic selection method | |
| CN103049561B (en) | A kind of data compression method, storage engines and storage system | |
| CN108682047B (en) | Three-dimensional building model file conversion method and system | |
| CN112149826B (en) | Profile graph-based optimization method in deep neural network inference calculation | |
| CN104618715A (en) | Method and device for obtaining minimal rate-distortion cost | |
| CN114036319A (en) | Power knowledge extraction method, system, device and storage medium | |
| CN111723905A (en) | Automatic convolution parameter configuration method based on linear programming | |
| CN114065923A (en) | Compression method, system and accelerating device of convolutional neural network | |
| CN110888860A (en) | Model feature processing method and device, computer equipment and storage medium | |
| CN116089095B (en) | Deployment method for ReRAM neural network computing engine network | |
| CN110162021A (en) | A kind of cybernetics system method for testing performance | |
| CN117829242B (en) | Model processing method and related equipment | |
| CN102902876B (en) | A kind of calculated crosswise method and apparatus of raster data |
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 |