CN110223328A - A kind of improvement particle image velocimetry robustness optical flow approach based on physics - Google Patents
A kind of improvement particle image velocimetry robustness optical flow approach based on physics Download PDFInfo
- Publication number
- CN110223328A CN110223328A CN201910397141.4A CN201910397141A CN110223328A CN 110223328 A CN110223328 A CN 110223328A CN 201910397141 A CN201910397141 A CN 201910397141A CN 110223328 A CN110223328 A CN 110223328A
- Authority
- CN
- China
- Prior art keywords
- optical flow
- building
- physics
- robustness
- energy functional
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
- G06T7/269—Analysis of motion using gradient-based methods
Landscapes
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
Abstract
The improvement particle image velocimetry robustness optical flow approach based on physics that the invention discloses a kind of, includes the following steps: the building of S1, steady energy functional;S2, the building for checking data item;The selection of S3, smooth item;The building of S4, penalty;S5, the processing of model energy functional minimization;S6, optical flow computation;The present invention provides a kind of improvement particle image velocimetry robustness optical flow approach based on physics, and present invention encompasses the buildings of steady energy functional;Check the building of data item;The selection of smooth item;The building of penalty;The processing of model energy functional minimization;Optical flow computation loose particles image speed measurement robustness is high;In the present invention further in the building process of the steady energy functional, the optical flow constraint equation based on physics is constructed first;Then change for optical flow constraint equation partial differential;Last to generate steady energy functional according to result of variations, energy functional versatility is high.
Description
Technical field
The invention belongs to Particle Image Velocimetry field, specially a kind of improvement particle image velocimetry based on physics
Robustness optical flow approach.
Background technique
Particle image velocimetry is that a kind of use repeatedly images to record the position of particle in flow field, and analyzes the image taken the photograph,
Method to measure flowing velocity.The new development tradition FLOW VISUALIZATION for being substantially Flow visualisation technology is Experimental Flowing Object power
The important component learned, its main task is that certain properties of flowing are subject to visual representation, to obtain to flowing
All-round developing understanding is obtained, because having formed a flourishing long time project in experimental fluid mechanics;
For application aspect, fluid velocity field measurement is for understanding complicated stream in hydrodynamics and aerodynamic studies
Body is of great significance.Fluid motion is a kind of typical non-rigid motion, and calculating need to be based on image processing techniques.It can pass through
Partial fluid motion vector size, direction and distribution situation are obtained to dynamic image sequence analysis, and then can be obtained such as sticky
And the physical characteristics such as vortex field distribution.Under normal conditions, due to moving object is transparent or not easily pass through optical device observation, need by
Visible particles are placed in measured object, and to obtain fluid motion feature indirectly by estimation particle motion vector, this calls particle picture
It tests the speed.
However traditional particle image velocimetry optical flow approach, robustness is poor, optical flow computation low efficiency and accuracy is undesirable,
For this purpose, we are proposed a kind of improvement particle image velocimetry robustness optical flow approach based on physics.
Summary of the invention
It is an object of the invention to: in order to solve traditional particle image velocimetry optical flow approach, robustness is poor, optical flow computation
Low efficiency and the undesirable technical problem of accuracy, provide a kind of improvement particle image velocimetry robustness light stream based on physics
Method.
The technical solution adopted by the invention is as follows:
A kind of improvement particle image velocimetry robustness optical flow approach based on physics, includes the following steps:
The building of S1, steady energy functional;
S2, the building for checking data item;
The selection of S3, smooth item;
The building of S4, penalty;
S5, the processing of model energy functional minimization;
S6, optical flow computation.
Wherein, in the S1 steady energy functional building, the specific steps of which are as follows:
S101, optical flow constraint equation of the building based on physics;
S102, change for optical flow constraint equation partial differential;
S103, steady energy functional is generated according to result of variations.
Wherein, the building of data item is checked in the S2, the specific steps of which are as follows:
S201, luminance picture data are obtained;
S202, gradient image data are obtained;
S203, color difference image data are obtained;
S204, the complex light stream data model for constructing brightness, gradient and color difference.
Wherein, in the S3 smooth item selection, the specific steps of which are as follows:
S301, global smoothing model building;
The verifying and amendment of S302, global smoothing model.
Wherein, in the S4 penalty building, the specific steps of which are as follows:
S401, secondary penalty is introduced;
The full variation transformation of S402, secondary penalty;
S403, robust variation penalty is obtained.
Wherein, model energy functional minimization is handled in the S5, the specific steps of which are as follows:
S501, basic minimization algorithm equation is introduced;
S502, robust variation penalty is imported;
S503, overrelaxation iterative calculation obtain optical flow field.
Wherein, optical flow computation in the S6, the specific steps of which are as follows:
S601, optical flow field is layered using pyramid decimation factor;
S602, up-sampling interpolation calculation;
S603, image transformation and inter-layer data processing are approached using linearisation.
Wherein, the pyramid decimation factor value range is 0.9~0.95.
Wherein, the up-sampling interpolation calculation is calculated using bicubic interpolation.
In conclusion by adopting the above-described technical solution, the beneficial effects of the present invention are:
1, the present invention provides a kind of improvement particle image velocimetry robustness optical flow approach based on physics, and the present invention covers
The building of steady energy functional;Check the building of data item;The selection of smooth item;The building of penalty;Model energy is general
The processing of letter minimization;Optical flow computation loose particles image speed measurement robustness is high.
2, in the present invention further in the building process of the steady energy functional, light of the building based on physics first
Flow constraint equation;Then change for optical flow constraint equation partial differential;It is last to generate steady energy functional, energy according to result of variations
It is high to measure functional versatility.
3, secondary penalty is firstly introduced into further in the building process of the penalty in the present invention;Then
The full variation of secondary penalty converts;Finally obtain robust variation penalty, penalty high treating effect.
4, in the present invention further during the optical flow computation, first with pyramid decimation factor to optical flow field point
Layer;Then up-sampling interpolation calculation;Image transformation finally is approached using linearisation and inter-layer data is handled;The wherein pyramid
Decimation factor value range is 0.9~0.95;The up-sampling interpolation calculation is calculated using bicubic interpolation;Computational efficiency it is high and
Accuracy is high.
Detailed description of the invention
Fig. 1 is process simplified schematic diagram of the invention;
Fig. 2 is the process simplified schematic diagram of the building of steady energy functional in the present invention;
Fig. 3 is the process simplified schematic diagram that the building of data item is checked in the present invention;
Fig. 4 is the process simplified schematic diagram of the selection of smooth item in the present invention;
Fig. 5 is the process simplified schematic diagram of the building of penalty in the present invention;
Fig. 6 is the process simplified schematic diagram of model energy functional minimization processing in the present invention;
Fig. 7 is the process simplified schematic diagram of optical flow computation in the present invention.
Specific embodiment
In order to make the objectives, technical solutions, and advantages of the present invention clearer, with reference to the accompanying drawings and embodiments, right
The present invention is further elaborated.It should be appreciated that described herein, specific examples are only used to explain the present invention, not
For limiting the present invention.
In the description of the present invention, it should be noted that term " center ", "upper", "lower", "left", "right", "vertical",
The orientation or positional relationship of the instructions such as "horizontal", "inner", "outside" be based on the orientation or positional relationship shown in the drawings, merely to
Convenient for description the present invention and simplify description, rather than the device or element of indication or suggestion meaning must have a particular orientation,
It is constructed and operated in a specific orientation, therefore is not considered as limiting the invention;Term " first ", " second ", " third "
It is used for descriptive purposes only and cannot be understood as indicating or suggesting relative importance;In addition, unless otherwise specific regulation and limit
Fixed, term " installation ", " connected ", " connection " shall be understood in a broad sense, for example, it may be being fixedly connected, be also possible to detachably connect
It connects, or is integrally connected;It can be mechanical connection, be also possible to be electrically connected;It can be directly connected, intermediate matchmaker can also be passed through
Jie is indirectly connected, and can be the connection inside two elements.It for the ordinary skill in the art, can be with concrete condition
Understand the concrete meaning of above-mentioned term in the present invention.
Embodiment one, referring to Fig.1~3, a kind of improvement particle image velocimetry robustness optical flow approach based on physics, packet
Include following steps:
The building of S1, steady energy functional;
S2, the building for checking data item;
The selection of S3, smooth item;
The building of S4, penalty;
S5, the processing of model energy functional minimization;
S6, optical flow computation.
Wherein, in the S1 steady energy functional building, the specific steps of which are as follows:
S101, optical flow constraint equation of the building based on physics;
S102, change for optical flow constraint equation partial differential;
S103, steady energy functional is generated according to result of variations.
Wherein, the building of data item is checked in the S2, the specific steps of which are as follows:
S201, luminance picture data are obtained;
S202, gradient image data are obtained;
S203, color difference image data are obtained;
S204, the complex light stream data model for constructing brightness, gradient and color difference.
Embodiment two, referring to Fig.1~4, a kind of improvement particle image velocimetry robustness optical flow approach based on physics, packet
Include following steps:
The building of S1, steady energy functional;
S2, the building for checking data item;
The selection of S3, smooth item;
The building of S4, penalty;
S5, the processing of model energy functional minimization;
S6, optical flow computation.
Wherein, in the S1 steady energy functional building, the specific steps of which are as follows:
S101, optical flow constraint equation of the building based on physics;
S102, change for optical flow constraint equation partial differential;
S103, steady energy functional is generated according to result of variations.
Wherein, the building of data item is checked in the S2, the specific steps of which are as follows:
S201, luminance picture data are obtained;
S202, gradient image data are obtained;
S203, color difference image data are obtained;
S204, the complex light stream data model for constructing brightness, gradient and color difference.
Wherein, in the S3 smooth item selection, the specific steps of which are as follows:
S301, global smoothing model building;
The verifying and amendment of S302, global smoothing model.
Embodiment three, referring to Fig.1~5, a kind of improvement particle image velocimetry robustness optical flow approach based on physics, packet
Include following steps:
The building of S1, steady energy functional;
S2, the building for checking data item;
The selection of S3, smooth item;
The building of S4, penalty;
S5, the processing of model energy functional minimization;
S6, optical flow computation.
Wherein, in the S1 steady energy functional building, the specific steps of which are as follows:
S101, optical flow constraint equation of the building based on physics;
S102, change for optical flow constraint equation partial differential;
S103, steady energy functional is generated according to result of variations.
Wherein, the building of data item is checked in the S2, the specific steps of which are as follows:
S201, luminance picture data are obtained;
S202, gradient image data are obtained;
S203, color difference image data are obtained;
S204, the complex light stream data model for constructing brightness, gradient and color difference.
Wherein, in the S3 smooth item selection, the specific steps of which are as follows:
S301, global smoothing model building;
The verifying and amendment of S302, global smoothing model.
Wherein, in the S4 penalty building, the specific steps of which are as follows:
S401, secondary penalty is introduced;
The full variation transformation of S402, secondary penalty;
S403, robust variation penalty is obtained.
Example IV, referring to Fig.1~6, a kind of improvement particle image velocimetry robustness optical flow approach based on physics, packet
Include following steps:
The building of S1, steady energy functional;
S2, the building for checking data item;
The selection of S3, smooth item;
The building of S4, penalty;
S5, the processing of model energy functional minimization;
S6, optical flow computation.
Wherein, in the S1 steady energy functional building, the specific steps of which are as follows:
S101, optical flow constraint equation of the building based on physics;
S102, change for optical flow constraint equation partial differential;
S103, steady energy functional is generated according to result of variations.
Wherein, the building of data item is checked in the S2, the specific steps of which are as follows:
S201, luminance picture data are obtained;
S202, gradient image data are obtained;
S203, color difference image data are obtained;
S204, the complex light stream data model for constructing brightness, gradient and color difference.
Wherein, in the S3 smooth item selection, the specific steps of which are as follows:
S301, global smoothing model building;
The verifying and amendment of S302, global smoothing model.
Wherein, in the S4 penalty building, the specific steps of which are as follows:
S401, secondary penalty is introduced;
The full variation transformation of S402, secondary penalty;
S403, robust variation penalty is obtained.
Wherein, model energy functional minimization is handled in the S5, the specific steps of which are as follows:
S501, basic minimization algorithm equation is introduced;
S502, robust variation penalty is imported;
S503, overrelaxation iterative calculation obtain optical flow field.
Embodiment five, referring to Fig.1~7, a kind of improvement particle image velocimetry robustness optical flow approach based on physics, packet
Include following steps: the building of S1, steady energy functional;
S2, the building for checking data item;
The selection of S3, smooth item;
The building of S4, penalty;
S5, the processing of model energy functional minimization;
S6, optical flow computation.
Wherein, in the S1 steady energy functional building, the specific steps of which are as follows:
S101, optical flow constraint equation of the building based on physics;
S102, change for optical flow constraint equation partial differential;
S103, steady energy functional is generated according to result of variations.
Wherein, the building of data item is checked in the S2, the specific steps of which are as follows:
S201, luminance picture data are obtained;
S202, gradient image data are obtained;
S203, color difference image data are obtained;
S204, the complex light stream data model for constructing brightness, gradient and color difference.
Wherein, in the S3 smooth item selection, the specific steps of which are as follows:
S301, global smoothing model building;
The verifying and amendment of S302, global smoothing model.
Wherein, in the S4 penalty building, the specific steps of which are as follows:
S401, secondary penalty is introduced;
The full variation transformation of S402, secondary penalty;
S403, robust variation penalty is obtained.
Wherein, model energy functional minimization is handled in the S5, the specific steps of which are as follows:
S501, basic minimization algorithm equation is introduced;
S502, robust variation penalty is imported;
S503, overrelaxation iterative calculation obtain optical flow field.
Wherein, optical flow computation in the S6, the specific steps of which are as follows:
S601, optical flow field is layered using pyramid decimation factor;
S602, up-sampling interpolation calculation;
S603, image transformation and inter-layer data processing are approached using linearisation.
Wherein, the pyramid decimation factor value is 0.9;Wherein, the up-sampling interpolation calculation is inserted using bicubic
Value calculates.
Working principle: the present invention provides a kind of improvement particle image velocimetry robustness optical flow approach based on physics, this
Invention covers the building of steady energy functional;Check the building of data item;The selection of smooth item;The building of penalty;Mould
The processing of type energy functional minimization;Optical flow computation loose particles image speed measurement robustness is high;Further described steady in the present invention
In the building process of strong energy functional, the optical flow constraint equation based on physics is constructed first;Then it is directed to optical flow constraint equation
Partial differential variation;Last to generate steady energy functional according to result of variations, energy functional versatility is high;Further exist in the present invention
In the building process of the penalty, it is firstly introduced into secondary penalty;The then full variation transformation of secondary penalty;Most
Robust variation penalty, penalty high treating effect are obtained afterwards;In the present invention further during the optical flow computation,
Optical flow field is layered first with pyramid decimation factor;Then up-sampling interpolation calculation;Finally image is approached using linearisation
Transformation and inter-layer data processing;Wherein the pyramid decimation factor value range is 0.9~0.95;The up-sampling interpolation meter
It calculates and is calculated using bicubic interpolation;Computational efficiency is high and accuracy is high.
The foregoing is merely illustrative of the preferred embodiments of the present invention, is not intended to limit the invention, all in essence of the invention
Made any modifications, equivalent replacements, and improvements etc., should all be included in the protection scope of the present invention within mind and principle.
Claims (9)
1. a kind of improvement particle image velocimetry robustness optical flow approach based on physics, characterized by the following steps:
The building of S1, steady energy functional;
S2, the building for checking data item;
The selection of S3, smooth item;
The building of S4, penalty;
S5, the processing of model energy functional minimization;
S6, optical flow computation.
2. a kind of improvement particle image velocimetry robustness optical flow approach based on physics as described in claim 1, feature
It is: the building of steady energy functional in the S1, the specific steps of which are as follows:
S101, optical flow constraint equation of the building based on physics;
S102, change for optical flow constraint equation partial differential;
S103, steady energy functional is generated according to result of variations.
3. a kind of improvement particle image velocimetry robustness optical flow approach based on physics as claimed in claim 2, feature
It is: checks the building of data item in the S2, the specific steps of which are as follows:
S201, luminance picture data are obtained;
S202, gradient image data are obtained;
S203, color difference image data are obtained;
S204, the complex light stream data model for constructing brightness, gradient and color difference.
4. a kind of improvement particle image velocimetry robustness optical flow approach based on physics as claimed in claim 3, feature
It is: the selection of smooth item in the S3, the specific steps of which are as follows:
S301, global smoothing model building;
The verifying and amendment of S302, global smoothing model.
5. a kind of improvement particle image velocimetry robustness optical flow approach based on physics as described in claim 1, feature
It is: the building of penalty in the S4, the specific steps of which are as follows:
S401, secondary penalty is introduced;
The full variation transformation of S402, secondary penalty;
S403, robust variation penalty is obtained.
6. a kind of improvement particle image velocimetry robustness optical flow approach based on physics a method as claimed in any one of claims 1 to 5,
It is characterized by: model energy functional minimization is handled in the S5, the specific steps of which are as follows:
S501, basic minimization algorithm equation is introduced;
S502, robust variation penalty is imported;
S503, overrelaxation iterative calculation obtain optical flow field.
7. a kind of improvement particle image velocimetry robustness optical flow approach based on physics as claimed in claim 6, feature
It is: optical flow computation in the S6, the specific steps of which are as follows:
S601, optical flow field is layered using pyramid decimation factor;
S602, up-sampling interpolation calculation;
S603, image transformation and inter-layer data processing are approached using linearisation.
8. a kind of improvement particle image velocimetry robustness optical flow approach based on physics as claimed in claim 7, feature
Be: the pyramid decimation factor value range is 0.9~0.95.
9. a kind of improvement particle image velocimetry robustness optical flow approach based on physics as claimed in claim 7, feature
Be: the up-sampling interpolation calculation is calculated using bicubic interpolation.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910397141.4A CN110223328A (en) | 2019-05-14 | 2019-05-14 | A kind of improvement particle image velocimetry robustness optical flow approach based on physics |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910397141.4A CN110223328A (en) | 2019-05-14 | 2019-05-14 | A kind of improvement particle image velocimetry robustness optical flow approach based on physics |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110223328A true CN110223328A (en) | 2019-09-10 |
Family
ID=67821019
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910397141.4A Pending CN110223328A (en) | 2019-05-14 | 2019-05-14 | A kind of improvement particle image velocimetry robustness optical flow approach based on physics |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110223328A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116070550A (en) * | 2023-03-07 | 2023-05-05 | 浙江大学 | Improved reconstruction flow field pressure field method based on time-resolved PIV |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101629966A (en) * | 2009-08-18 | 2010-01-20 | 清华大学深圳研究生院 | Particle image velocimetry (PIV) processing method |
CN107818572A (en) * | 2017-10-16 | 2018-03-20 | 大连大学 | Improvement particle image velocimetry robustness optical flow approach based on physics |
CN107958464A (en) * | 2017-12-05 | 2018-04-24 | 大连大学 | The quick Robust Estimate Method of optical flow field under noise circumstance |
EP3367106A1 (en) * | 2017-02-28 | 2018-08-29 | King Abdullah University Of Science And Technology | Rainbow particle imaging velocimetry for dense 3d fluid velocity imaging |
CN109615650A (en) * | 2018-11-22 | 2019-04-12 | 浙江工商大学 | It is a kind of based on the calculus of variations and to block complementary light field stream estimation method |
-
2019
- 2019-05-14 CN CN201910397141.4A patent/CN110223328A/en active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101629966A (en) * | 2009-08-18 | 2010-01-20 | 清华大学深圳研究生院 | Particle image velocimetry (PIV) processing method |
EP3367106A1 (en) * | 2017-02-28 | 2018-08-29 | King Abdullah University Of Science And Technology | Rainbow particle imaging velocimetry for dense 3d fluid velocity imaging |
CN107818572A (en) * | 2017-10-16 | 2018-03-20 | 大连大学 | Improvement particle image velocimetry robustness optical flow approach based on physics |
CN107958464A (en) * | 2017-12-05 | 2018-04-24 | 大连大学 | The quick Robust Estimate Method of optical flow field under noise circumstance |
CN109615650A (en) * | 2018-11-22 | 2019-04-12 | 浙江工商大学 | It is a kind of based on the calculus of variations and to block complementary light field stream estimation method |
Non-Patent Citations (1)
Title |
---|
杜军平,梁美玉,訾玲玲: "《跨尺度运动图像的插值、增强和重建》", 30 April 2019 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116070550A (en) * | 2023-03-07 | 2023-05-05 | 浙江大学 | Improved reconstruction flow field pressure field method based on time-resolved PIV |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Bellani et al. | Shape effects on turbulent modulation by large nearly neutrally buoyant particles | |
Shan et al. | Effects of the orifice to pipe diameter ratio on orifice flows | |
CN106596037B (en) | Model in wind tunnel streams the video measuring method of Intensity Projection field | |
Fujita et al. | Unseeded and seeded PIV measurements of river flows videotaped from a helicopter | |
JP5165121B2 (en) | How to measure flow status | |
Örlü et al. | Instantaneous wall-shear-stress measurements: advances and application to near-wall extreme events | |
CN111693999B (en) | Multi-sensor fusion wind speed and direction estimation method based on radar wind measurement combination strategy | |
CN113689473B (en) | River channel surface flow velocity calculation method and device based on variational principle | |
CN110223328A (en) | A kind of improvement particle image velocimetry robustness optical flow approach based on physics | |
Gargiulo et al. | Examination of flow sensitivities in turbulence model validation experiments | |
Shi et al. | 3D surface pressure measurement with single light-field camera and pressure-sensitive paint | |
Zhong et al. | Skin-friction topology of wing–body junction flows | |
CN110207927A (en) | A kind of wind tunnel wall interference correction correction method of airfoil wind tunnel test | |
CN109489745B (en) | Flow measuring method based on data iteration | |
Jones et al. | Ionospheric effects of magneto-acoustic-gravity waves: dispersion relation | |
Du et al. | Study on the reconstruction method of stereo vision in glass flume | |
KR101721164B1 (en) | Method and Apparatus for Determining Leak Location of Waterfront Structures | |
Johnson | Computation of flow in a vortex-shedding flowmeter | |
Chang et al. | Portable Particle Image Velocimetry Measurement Using a Laser-Based Technique | |
Huang et al. | Velocity Measurements of Hydrogen Jets Using the Optical Flow Method | |
Hu et al. | Adaptive measurement for automated field reconstruction and calibration of magnetic systems | |
Zappa et al. | Uncertainty budget in PSV technique measurements | |
JP2014153093A (en) | Calibration system of electromagnetic current meter, computation system of calibration coefficient, and method thereof | |
Sous et al. | A Tsai’s model based S-PIV method for velocity measurements in a turbulent Ekman layer | |
Liu et al. | Inverse problems in global flow diagnostics |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190910 |
|
RJ01 | Rejection of invention patent application after publication |