CN103247058B - A kind of quick the Computation of Optical Flow based on error Distributed-tier grid - Google Patents
A kind of quick the Computation of Optical Flow based on error Distributed-tier grid Download PDFInfo
- Publication number
- CN103247058B CN103247058B CN201310174158.6A CN201310174158A CN103247058B CN 103247058 B CN103247058 B CN 103247058B CN 201310174158 A CN201310174158 A CN 201310174158A CN 103247058 B CN103247058 B CN 103247058B
- Authority
- CN
- China
- Prior art keywords
- alpha
- error
- layer
- formula
- image
- 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.)
- Expired - Fee Related
Links
- 230000003287 optical effect Effects 0.000 title claims abstract description 17
- 238000000034 method Methods 0.000 claims abstract description 27
- 238000012937 correction Methods 0.000 claims abstract description 9
- 238000005457 optimization Methods 0.000 claims abstract description 6
- 238000004800 variational method Methods 0.000 claims abstract description 5
- 230000005540 biological transmission Effects 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 20
- 230000008569 process Effects 0.000 claims description 7
- 238000013507 mapping Methods 0.000 claims description 4
- 238000012804 iterative process Methods 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 2
- 238000006073 displacement reaction Methods 0.000 claims description 2
- 238000011156 evaluation Methods 0.000 claims description 2
- 230000001105 regulatory effect Effects 0.000 claims 1
- 230000000007 visual effect Effects 0.000 abstract description 2
- 238000012360 testing method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 2
- 230000004075 alteration Effects 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 230000033001 locomotion Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000009012 visual motion Effects 0.000 description 1
Landscapes
- Image Processing (AREA)
Abstract
The invention belongs to field of machine vision, disclose a kind of quick the Computation of Optical Flow based on error Distributed-tier grid, comprising: step one, input picture, build linear equation Ax=f; Step 2, sets up multi-layer image pyramid; Step 3, carries out front optimization, eliminates high fdrequency component; Step 4, carries out residual error transmission, eliminates low frequency component; Step 5, repeats step 3, four, until residual error is passed to most coarse layer; Step 7, carries out error passback from most coarse layer; Step 8, carries out refined net error correction; Step 9, carries out rear Optimized Iterative, improves stability of solution; Step 10, repeats step 7 ~ nine, until residual error is passed to most sub-layers.Method of the present invention is a kind of effective ways of accelerate equation Optimization Solution, and high frequency error convergence is very fast, and can significantly improve the computing velocity in visual light flow field, compared with the variational method, the speed of convergence of the method can bring up to more than 3.5 times.
Description
Technical field
The invention belongs to field of machine vision, relate to a kind of quick the Computation of Optical Flow based on error Distributed-tier grid.
Background technology
In visual motion analysis theory, when having a relative motion when between camera and scene objects, viewed brightness of image pattern further is referred to as light stream (optical flow).Because the light stream vector field between sequence video image have expressed the change of image, can matching relationship quantitatively between descriptor frame between corresponding image points, and provide image planes side-play amount, therefore for observer provides motion about target and structural information.
The variational method is the main stream approach solving light stream.Its guiding theory is level and smooth in gradient and the meaning of data fidelity builds energy functional, and emphasis obtains fine and close light stream vector field by the Euler-Lagrange equation (partial differential equation) solving optimized energy function.The advantage of variation light stream is, the versatility of model is good, and the degree of accuracy solved is high, and optical flow field is fine and close.Under variation framework, the method that Horn and Schunck proposes can obtain different optical flow computation models according to different assumed condition, its energy functional model comprises level and smooth item and data item, can mathematics expression be:
In formula, I
0, I
1represent adjacent two two field pictures of camera acquisition respectively, x=(x, y)
trepresent certain picture point on image, with u (x)=(u (x), v (x))
trepresent the light stream vector of this position, wherein u (x) and v (x) is respectively the horizontal and vertical component of this vector.Symbol ▽ is gradient operator, and Ω represents imaging plane.Integration item is above level and smooth item, and integration item is below data item, and α is the constant coefficient of both adjustments weight.
By above formula gray scale conservation item I
1(x+u (x))-I
0x () carries out Taylor expansion, obtain:
I
1(x+u(x))-I
0(x)=I
xu+I
yv+I
t(2)
In formula, x, y represent level and the vertical component of image I respectively, and t represents the time.Lower right corner subscript represents corresponding partial derivative, namely
Bring formula (2) into formula (1), above-mentioned energy functional is converted to:
Utilize variational method, obtaining its corresponding Euler-Lagrange equation is:
In formula, u and v is respectively the horizontal and vertical component of light stream vector, and symbol Δ represents Laplace operator.After this system of equations discretize, if with traditional numerical method, as Gauss – Seidel, the methods such as SOR need iteration just can try to achieve comparatively ideal result thousands of times, and therefore the real-time of algorithm is poor.
In conventional numeric iterative algorithm, the convergence process of analytical error can find, high frequency error component is decayed rapidly in an iterative process, and the error component being in low frequency is decayed slowly.Therefore, the same numerical problem of rapid solving can be carried out with the multi-layer net that a series of resolution is different.This is because, refined net utilize less iterations can eliminate high frequency error; The low frequency aberration left over changes high frequency error into relative to coarse grid, will map on thicker grid, can continue on coarse grid, utilize less iterations to eliminate this part low frequency aberration, reach the object of accelerating convergence.Usual setting multi-layer net, transmits layer by layer by error, realizes the Distributed fusion of error.When specific implementation, formula (3) need be converted to system of linear equations.Formula (3) is out of shape a little, following form can be obtained:
β(x)=f (4)
Wherein:
Visible, not only comprise constant coefficient in function β, also comprise nonlinear Laplace operator
the calculating of optical flow field is caused can not directly to apply multi-layer net Algorithm for Solving thus.
Summary of the invention
According to the problems referred to above existed in optical flow computation method, the present invention proposes a kind of quick the Computation of Optical Flow based on error Distributed-tier grid, the speed of convergence of light stream iterative computation can be accelerated, improve the real-time of computer vision system.
Multi-layer net is a kind of effective ways of accelerate equation Optimization Solution.Wherein, equationof structure Ax=f is a primary step.Due in optical flow equation (4), containing nonlinear Laplace operator in function β, therefore function β is converted into not containing the coefficient matrices A of operator, to solve under multi-layer net algorithm frame by the present invention.
Based on a quick the Computation of Optical Flow for error Distributed-tier grid, it is characterized in that comprising the steps:
Step one, input picture, build linear equation Ax=f, method is as follows:
1) calculate spatio-temporal gradient tensor matrix, solving equations discretize will be treated;
Suppose in a certain two field picture plane, (i, j) represents the position of the i-th row jth row pixel, u
i,jand v
i,jrepresent the displacement field in this horizontal x-axis occurred between consecutive frame and longitudinal y-axis respectively.The space partial derivative f of computed image
x, f
yand time partial derivative f
t, the spatio-temporal gradient tensor J matrix of each pixel in image can be obtained:
If with [J]
i,jrepresent the value of J matrix at pixel (i, j) place, by h representation space step-length, then Euler-Lagrange equation (formula (3)) can discretely be:
2) difference equation (6) is changed into linear constant coefficient system of equations A
1x=f
1;
Abbreviation equation (6), obtains the DIFFERENCE EQUATIONS in the D of domain:
Always total N+1 is capable to suppose image, and remove zero row and N-th row, the DIFFERENCE EQUATIONS of each pixel of first row has following form:
In formula, g
ijrepresent the borderline element u of domain D
i, i.e. g
ij=u
ij.
If introduce optical flow components u
ijand v
ijon column vector:
And order:
For the J of the J matrix element at pixel (m, n) place
mn, matrix [J in formula (9)
mn]
ibe equivalent to unit matrix and the J on N-1 rank
mnthe product of i-th this column vector of row.Therefore, the DIFFERENCE EQUATIONS (8) of each pixel of first row can turn to:
In formula, I is the unit matrix on N-1 rank, b
1for the column vector of definition in formula (6), G is N-1 rank matrix below:
In like manner, the difference equation of each pixel of secondary series is:
In formula:
All transformed by all for image row, then system of equations can be written as A
1x=f
1form, in formula:
3) difference equation (7) is changed into linear constant coefficient system of equations A
2x=f
2;
4) by A
1x=f
1with A
2x=f
2common formation system of equations:
Ax=f (10)
Step 2, sets up multi-layer image pyramid;
The pyramidal schematic diagram of multi-layer image as shown in Figure 2.A uppermost tomographic image is former image in different resolution, each tomographic image below represent successively reduce resolution images falls.Because pixel itself is uniform discrete, therefore can multi-layer image be regarded as multi-layer net.
Usually, the size (long and wide) of image is between 300-800.According to the scale factor η (0.5< η <0.95) of setting, former image in different resolution is successively reduced, and size is rounded, obtain multi-layer image pyramid.Usually, number of plies N=4 ~ 5 layer.The initial value of every layer of light stream is set to zero.
Step 3, carries out front optimization, eliminates high fdrequency component;
Current layer i (i=1,2 ..., N-1), take null matrix as light stream initial value to refined net Equation Iterative m time, to eliminate high fdrequency component.If x
ifor initial value, obtain approximate evaluation value
be designated as:
In formula, Relax
mrepresent Gauss-Seidel iterative process, m is the iterations of setting.
A
idetermined by formula (10); As i=1, r is determined by formula (10), i.e. r
1=f.Otherwise r equals the residual error that last layer is passed to this layer.The corresponding residual error of this layer is updated to:
Step 4, carries out residual error transmission, eliminates low frequency component;
Front optimizing process is intended to eliminate high frequency error component, and inherited error is mainly low frequency component.Therefore, residual error is passed on coarse grid, to eliminate low frequency component.Be restricted to the surplus r on coarse grid
i+1for:
In formula,
for refined net is to the mapping operator on coarse grid.
Step 5, as i=1 ~ N-1, repeats step 3, four, until residual error is passed to most coarse layer;
Step 6, solves the system of equations on coarse grid:
A
Ne
N=r
N
In formula, e
nfor the light stream error of Exact Solution in most coarse layer.
Obtained by matrix operation:
e
N=(A
N)
-1r
N
Step 7, carries out error passback from most coarse layer (j=N);
Current layer j, passes back to last layer comparatively on refined net j-1, that is: by the amount of error correction of trying to achieve
In formula,
for coarse grid is to the mapping operator on refined net.
Step 8, carries out refined net error correction;
The initial value of calculating is added the error passed back by coarse grid, obtain the solution after refined net correction
that is:
Step 9, carries out rear Optimized Iterative, improves stability of solution;
After refined net corrects, with the solution after correcting for initial value, Optimized Iterative n time after performing, obtains the corrected value after renewal
its expression formula is:
Step 10, works as j=N, N-1 ..., when 2, repeat step 7 ~ nine, until residual error is passed to most sub-layers.
Compared with prior art, the present invention has following beneficial effect:
Light stream solving equation launches by row by the present invention, is turned to the form of system of linear equations, makes the complicated calculations process of light stream can accelerate to realize by linear multi-layer trellis algorithm.Secondly, have employed the multi-layer net of different resolution, the error component of different frequency is distributed on different grids, make error become high frequency error relative to this layer of grid.Very fast due to high frequency error convergence, thus the method can significantly improve the computing velocity in visual light flow field, and compared with the variational method, the speed of convergence of the method can bring up to more than 3.5 times.
Accompanying drawing explanation
Fig. 1 is the process flow diagram of method involved in the present invention;
Fig. 2 is multi-layer image pyramid schematic diagram;
Fig. 3 is first group of test pattern RubberWhale of application example of the present invention;
Fig. 4 is first group of test pattern of application example of the present invention.
Embodiment
Below in conjunction with the drawings and specific embodiments, the present invention will be further described.
The microcomputer one implementing to rely on P3800M dominant frequency, 512M internal memory and above configuration of the present invention, as hardware device platform.Running environment is Windows XP operating system and Matlab software platform.
The detailed process of present embodiment is, first according to optical flow estimation to be solved, and structure system of linear equations Ax=f.Then set up 4 ~ 5 tomographic image pyramids, be optimized iteration respectively at every one deck, and iteration residual error is successively mapped on the thicker grid of lower one deck.In most coarse layer, solve the system of equations on most coarse grid, obtain accurate estimation of error.After this, return error in opposite direction, thinner grid successively carries out error correction.Fig. 1 is the process flow diagram of method involved in the present invention, specifically comprises the following steps:
Step one, input picture, builds linear equation Ax=f;
Step 2, sets up multi-layer image pyramid;
Step 3, carries out front optimization, eliminates high fdrequency component;
Step 4, carries out residual error transmission, eliminates low frequency component;
Step 5, as i=1 ~ N-1, repeats step 3, four, until residual error is passed to most coarse layer;
Step 6, solves the system of equations on coarse grid: A
ne
n=r
n
Step 7, carries out error passback from most coarse layer (j=N);
Step 8, carries out refined net error correction;
Step 9, carries out rear Optimized Iterative, improves stability of solution;
Step 10, works as j=N, N-1 ..., when 2, repeat step 7 ~ nine, until residual error is passed to most sub-layers.
Provide an application example of the present invention below.
Two groups of test patterns are selected to verify the fast light flow field algorithm based on error Distributed-tier grid that the present invention proposes.One group is a pair test pattern RubberWhale in the Middlebury java standard library extensively adopted in the world, as shown in Figure 3; Another group is scientific paper " the An improved algorithm for TV-L that the people such as Wedel deliver at them
1optical flow " middle a pair test pattern adopted, as shown in Figure 4.
During experiment, the fast light flow field algorithm based on error Distributed-tier grid adopting traditional variational algorithm and the present invention to propose respectively calculate often organize image between optical flow field.By the working time of software statistics two kinds of methods, experimental result is as shown in table 1.In table 1, experiment one adopts image shown in Fig. 3, and experiment two adopts image shown in Fig. 4.As shown in Table 1, the operation time that the fast light flow field algorithm based on error Distributed-tier grid compares traditional variational algorithm shortens greatly, and the speed of convergence of two experiments all brings up to more than 3.5 times, significantly improves counting yield.
The contrast of the traditional variational algorithm of table 1 and the method for the invention working time
Claims (1)
1., based on a quick the Computation of Optical Flow for error Distributed-tier grid, it is characterized in that comprising the steps:
Step one, input picture, build linear equation Ax=f, method is as follows:
1) calculate spatio-temporal gradient tensor matrix, solving equations discretize will be treated;
Suppose in a certain two field picture plane, (i, j) represents the position of the i-th row jth row pixel, u
i,jand v
i,jrepresent the displacement field in this horizontal x-axis occurred between consecutive frame and longitudinal y-axis respectively; The space partial derivative f of computed image
x, f
yand time partial derivative f
t, the spatio-temporal gradient tensor J matrix of each pixel in image can be obtained:
If with [J]
i,jrepresent the value of J matrix at pixel (i, j) place, by h representation space step-length, the Euler-Lagrange equation variational method obtained by energy functional number is discrete is:
In formula, α is for regulating the coefficient of both level and smooth item and data item weight;
2) difference equation (2) is changed into linear constant coefficient system of equations A
1x=f
1;
Abbreviation equation (2), obtains the DIFFERENCE EQUATIONS in the D of domain:
Always total N+1 is capable to suppose image, and remove zero row and N-th row, the DIFFERENCE EQUATIONS of each pixel of first row has following form:
In formula, g
ijrepresent the borderline element u of domain D
i, i.e. g
ij=u
ij;
Introduce optical flow components u
ijand v
ijon column vector:
And order:
For the J of the J matrix element at pixel (m, n) place
mn, matrix [J in formula (5)
mn]
ibe equivalent to unit matrix and the J on N-1 rank
mnthe product of i-th this column vector of row; Therefore, the DIFFERENCE EQUATIONS (4) of each pixel of first row can turn to:
In formula, I is N-1 rank unit matrix, b
1for the column vector of definition in formula (2), G is N-1 rank matrix below:
The difference equation of each pixel of secondary series is:
In formula:
All transformed by all for image row, then system of equations can be written as A
1x=f
1form, in formula:
3) difference equation (3) is changed into linear constant coefficient system of equations A
2x=f
2;
4) by A
1x=f
1with A
2x=f
2common formation system of equations:
Ax=f (6)
Step 2, set up multi-layer image pyramid, method is as follows:
The uppermost tomographic image of multi-layer image pyramid is former image in different resolution, each tomographic image below represent successively reduce resolution images falls; Because pixel itself is uniform discrete, therefore can multi-layer image be regarded as multi-layer net;
Usually, the length of image and wide between 300 ~ 800; According to the scale factor η of setting, 0.5< η <0.95, successively reduces former image in different resolution, and size is rounded, and obtains multi-layer image pyramid; Usually, number of plies N=4 ~ 5 layer, the initial value of every layer of light stream is set to zero;
Step 3, carries out front optimization, and eliminate high fdrequency component, method is as follows:
At current layer i, i=1,2 ..., N-1 take null matrix as light stream initial value to refined net Equation Iterative m time, to eliminate high fdrequency component; If x
ifor initial value, obtain approximate evaluation value
be designated as:
In formula, Relax
mrepresent Gauss-Seidel iterative process, m is the iterations of setting;
A
idetermined by formula (6); As i=1, r is determined by formula (6), i.e. r
1=f; Otherwise r equals the residual error that last layer is passed to this layer; The corresponding residual error of this layer is updated to:
Step 4, carries out residual error transmission, and eliminate low frequency component, method is as follows:
Front optimizing process is intended to eliminate high frequency error component, and inherited error is mainly low frequency component, is therefore passed on coarse grid by residual error, to eliminate low frequency component; Be restricted to the surplus r on coarse grid
i+1for:
In formula,
for refined net is to the mapping operator on coarse grid;
Step 5, as i=1 ~ N-1, repeats step 3, four, until residual error is passed to most coarse layer;
Step 6, solves the system of equations on coarse grid:
A
Ne
N=r
N
In formula, e
nfor the light stream error of Exact Solution in most coarse layer;
Obtained by matrix operation:
e
N=(A
N)
-1r
N
Step 7, from most coarse layer, carry out error passback, method is as follows:
Current layer is j, the amount of error correction of trying to achieve is passed back to last layer comparatively on refined net j-1, that is:
In formula,
for coarse grid is to the mapping operator on refined net;
Step 8, carry out refined net error correction, method is as follows:
The initial value of calculating is added the error passed back by coarse grid, obtain the solution after refined net correction
that is:
Step 9, carries out rear Optimized Iterative, and improve stability of solution, method is as follows:
After refined net corrects, with the solution after correcting for initial value, Optimized Iterative n time after performing, obtains the corrected value after renewal
its expression formula is:
Step 10, works as j=N, N-1 ..., when 2, repeat step 7 ~ nine, until residual error is passed to most sub-layers.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310174158.6A CN103247058B (en) | 2013-05-13 | 2013-05-13 | A kind of quick the Computation of Optical Flow based on error Distributed-tier grid |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310174158.6A CN103247058B (en) | 2013-05-13 | 2013-05-13 | A kind of quick the Computation of Optical Flow based on error Distributed-tier grid |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103247058A CN103247058A (en) | 2013-08-14 |
CN103247058B true CN103247058B (en) | 2015-08-19 |
Family
ID=48926563
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310174158.6A Expired - Fee Related CN103247058B (en) | 2013-05-13 | 2013-05-13 | A kind of quick the Computation of Optical Flow based on error Distributed-tier grid |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103247058B (en) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103618904B (en) * | 2013-11-20 | 2017-02-22 | 华为技术有限公司 | Motion estimation method and device based on pixels |
CN103700117A (en) * | 2013-11-21 | 2014-04-02 | 北京工业大学 | Robust optical flow field estimating method based on TV-L1 variation model |
CN103761749B (en) * | 2014-02-12 | 2016-08-17 | 北京工业大学 | A kind of optical flow field method for quick estimating based on non-linear multi-grid method |
CN104915968B (en) * | 2015-05-19 | 2017-07-21 | 北京航空航天大学 | A kind of optical flow velocity method of estimation based on navier stokes equations |
US10257501B2 (en) * | 2016-04-06 | 2019-04-09 | Facebook, Inc. | Efficient canvas view generation from intermediate views |
CN106385541A (en) * | 2016-09-30 | 2017-02-08 | 虹软(杭州)科技有限公司 | Method for realizing zooming through wide-angle photographing component and long-focus photographing component |
CN108399629B (en) * | 2018-02-05 | 2020-06-02 | 西南交通大学 | Image array optical flow estimation method of artificial compound eye camera |
CN109061554B (en) * | 2018-06-26 | 2022-07-15 | 哈尔滨工程大学 | Target arrival angle estimation method based on dynamic update of spatial discrete grid |
CN110322477B (en) * | 2019-06-10 | 2022-01-04 | 广州视源电子科技股份有限公司 | Feature point observation window setting method, tracking method, device, equipment and medium |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101770568A (en) * | 2008-12-31 | 2010-07-07 | 南京理工大学 | Target automatically recognizing and tracking method based on affine invariant point and optical flow calculation |
CN102509071A (en) * | 2011-10-14 | 2012-06-20 | 江南大学 | Light stream calculation system and method |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8553943B2 (en) * | 2011-06-14 | 2013-10-08 | Qualcomm Incorporated | Content-adaptive systems, methods and apparatus for determining optical flow |
-
2013
- 2013-05-13 CN CN201310174158.6A patent/CN103247058B/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101770568A (en) * | 2008-12-31 | 2010-07-07 | 南京理工大学 | Target automatically recognizing and tracking method based on affine invariant point and optical flow calculation |
CN102509071A (en) * | 2011-10-14 | 2012-06-20 | 江南大学 | Light stream calculation system and method |
Non-Patent Citations (5)
Title |
---|
Pierre-Marc Jodoin,Max Mignotte.Optical-flow based on an edge-avoidance procedure.《Computer Vision and Image Understanding》.2009,第113卷(第4期), * |
Variational Optical Flow Computation in Real Time;Andrés Bruhn, Joachim Weickert, Christian Feddern, Timo Kohlberg;《IEEE TRANSACTIONS ON IMAGE PROCESSING》;20050531;第14卷(第5期);608-615 * |
基于多重网格法的变分光流计算;王骁,王卫卫,聂秀志;《现代电子技术》;20080229(第2期);145-148,151 * |
昌猛,吴新根,罗立民,舒华忠,傅瑶.一种快速的光流场计算方法.《东南大学学报》.1999,第29卷(第5期), * |
谢剑斌,王晖,刘通,李沛秦.基于多重网格的光流场快速计算方法.《仪器仪表学报》.2008,第29卷(第4期), * |
Also Published As
Publication number | Publication date |
---|---|
CN103247058A (en) | 2013-08-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103247058B (en) | A kind of quick the Computation of Optical Flow based on error Distributed-tier grid | |
CN111523546B (en) | Image semantic segmentation method, system and computer storage medium | |
CN107229920B (en) | Behavior identification method based on integration depth typical time warping and related correction | |
CN106101535A (en) | A kind of based on local and the video stabilizing method of mass motion disparity compensation | |
CN107016649A (en) | A kind of vision data complementing method estimated based on local low-rank tensor | |
CN111881773A (en) | Event camera human body posture estimation method and system based on position deviation | |
CN104091350B (en) | A kind of object tracking methods of utilization motion blur information | |
CN104134196A (en) | Split Bregman weight iteration image blind restoration method based on non-convex higher-order total variation model | |
CN106408653A (en) | Real-time robust cluster adjustment method for large-scale three-dimensional reconstruction | |
CN110139046A (en) | A kind of video frame synthetic method based on tensor | |
CN109887029A (en) | A kind of monocular vision mileage measurement method based on color of image feature | |
CN113570658A (en) | Monocular video depth estimation method based on depth convolutional network | |
CN107862733A (en) | Large scale scene real-time three-dimensional method for reconstructing and system based on sight more new algorithm | |
Héas et al. | Bayesian estimation of turbulent motion | |
CN115797808A (en) | Unmanned aerial vehicle inspection defect image identification method, system, device and medium | |
CN115100294A (en) | Event camera calibration method, device and equipment based on linear features | |
CN103839280B (en) | A kind of human body attitude tracking of view-based access control model information | |
CN106934822A (en) | A kind of automobile workpiece non-rigid three-dimensional point cloud method for registering based on linear hybrid deformation | |
CN106874881A (en) | A kind of anti-joint sparse of part of multi-template space time correlation represents method for tracking target | |
CN104091364A (en) | Single-image super-resolution reconstruction method | |
CN105608719B (en) | A kind of rapid CT image rebuilding method based on two benches projection adjustment | |
CN104424625A (en) | GPU (graphics processing unit) acceleration CBCT image reconstruction method and device | |
CN117274349A (en) | Transparent object reconstruction method and system based on RGB-D camera consistency depth prediction | |
CN111553954A (en) | Direct method monocular SLAM-based online luminosity calibration method | |
CN116797640A (en) | Depth and 3D key point estimation method for intelligent companion line inspection device |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150819 |