CN107133953A - A kind of intrinsic image decomposition method learnt based on partial differential equation - Google Patents
A kind of intrinsic image decomposition method learnt based on partial differential equation Download PDFInfo
- Publication number
- CN107133953A CN107133953A CN201710138216.8A CN201710138216A CN107133953A CN 107133953 A CN107133953 A CN 107133953A CN 201710138216 A CN201710138216 A CN 201710138216A CN 107133953 A CN107133953 A CN 107133953A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- mtd
- function
- 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
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 21
- 230000011514 reflex Effects 0.000 claims abstract description 8
- 238000002939 conjugate gradient method Methods 0.000 claims abstract description 7
- 238000013519 translation Methods 0.000 claims description 13
- 238000012549 training Methods 0.000 claims description 7
- 230000004069 differentiation Effects 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 4
- 230000010365 information processing Effects 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 2
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 claims 1
- 238000005286 illumination Methods 0.000 abstract description 7
- 238000012545 processing Methods 0.000 abstract description 5
- 238000002945 steepest descent method Methods 0.000 abstract description 5
- 230000008859 change Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000004446 light reflex Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20081—Training; Learning
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
The present invention provides a kind of intrinsic image decomposition method learnt based on partial differential equation, and the present invention, independent of the prior-constrained of determination, and builds by the way of data-driven partial differential equation in the processing intrinsic Composition Estimation of image;The direction of search is determined using conjugate gradient method, relative to steepest descent method and Newton method, conjugacy is combined by this method with steepest descent method;One group of conjugate direction is constructed using the gradient at known point, and along this prescription to scanning for, obtains the minimal point of object function, to determine the optimal direction of search;This method can effectively realize that the image of different illumination conditions carries out reflex components and shade composition that eigen decomposition obtains the image.
Description
Technical field
The present invention relates to digital image processing field, more particularly, to a kind of based on the intrinsic of partial differential equation study
Picture breakdown method.
Background technology
The image that object in real world is presented in human eye depends on the intrinsic attribute of scene, the light of such as scene
According to, the shape of body surface, the material of body surface etc..Intrinsic image decomposition is a underlying issue in computer vision,
A given width input picture is, it is necessary to decomposite corresponding reflex components and shade composition.With the hair of digital image processing techniques
Exhibition, to image decompose obtaining its intrinsic composition, is played in computer vision and image processing field more and more important
Effect.Reflex components and shade that eigen decomposition obtains the image are carried out to the image with different illumination conditions under Same Scene
Composition is relatively difficult.
At present, on intrinsic image task resolution, traditional method is usually to design corresponding according to various prioris
Constrained optimization equation, acquisition reflex components and shade composition are solved by the optimization to target equation.Conventional a priori assumption
Have that the surface with the images of different illumination, input picture under Same Scene is smooth, the illumination of imaging surface color balancing, image
It is more natural etc..In addition, some current research algorithms are dependent on multiple input pictures, user's input or Depth cue
Extraneous information.Mainly for the input picture of the different illumination conditions under Same Scene, recover the anti-of input picture with this
Penetrate composition and shade composition.In place of above method comes with some shortcomings, existing method is typically based on a priori assumption, generally not
Possesses generality.
The content of the invention
The present invention provides a kind of intrinsic image decomposition method learnt based on partial differential equation, and this method, which can be realized, does not share the same light
Reflex components and shade composition that eigen decomposition obtains the image are carried out according to the image of condition.
In order to reach above-mentioned technique effect, technical scheme is as follows:
A kind of intrinsic image decomposition method learnt based on partial differential equation, is comprised the following steps:
S1:Input training data pair;
S2:Control function is initialized, because object function is non-convex, the convergence for minimizing process is initial towards depending on
The local minimum of change, step S3 is turned to when initialization procedure is not restrained, circulation is performed;Otherwise step S8 is turned to;
S3:Solve the optimal control equation constrained with PDE;
S4:Solve adjoint equation when adjoint function takes particular value;
S5:J=0,1 is calculated using following formula ..., the derivative of translation rotational invariants when 16:
Wherein, J is translation rotational invariants, and j is the number for translating rotational invariants, aj, bjIt is control function, λjAnd uj
It is positive weighting parameters,And φmIt is adjoint function, u is output image, and Ω is the rectangular area occupied by input picture, fΩ
Ω initial function, m=1,2 ..., M, M be the data sample logarithm of input, inv (u, v) represents to seek matrix (u, v)
Inverse, v is indicator function, and it is to collect the extensive information in image to be introduced into the purpose of indicator function, to instruct u differentiation.
S6:The direction of search is determined using conjugate gradient method;
S7:Golden section search is performed along the direction of search, and constantly updates system function, next circulation is carried out, directly
To j=16, it is trained;
S8:Terminate circulation, output system function;
S9:Prepare application data, data picture feature is background black, and object is single and prominent;
S10:Using obtained system function, to give data to carry out eigen decomposition application, obtain the reflex components of image
With shade composition.
Further, by solving (1) formula to control function a in the step S2j(t), t=0, Δ t, 1-
Δ t is initialized, and now fixes bj(t), j=0,1,16,
Fa (t)=d (1)
Wherein, ajAnd b (t)j(t) it is control function, a (t)={ aj} and b (t)={ b (t)j(t) } it is defined on Q
Function set, is respectively intended to control u and v differentiation, fuAnd fvIt is u and v initial function respectively.
Further, j-th of translation rotational invariants is calculated by introducing the adjoint equation of (2) formula in the step S3
Gateaux derivative, therefore local optimumWithCalculated and obtained by the algorithm based on gradient, for um
And vm, work as m=1, (2) formula solved during 2 ..., M:
Wherein, T is the time that PDE systems complete Vision information processing and output result, and Q is Ω × (0, T), and Γ is
Further, adjoint equation during adjoint function particular value is solved in the step S4:
ForAnd φm, m=1, during 2 ..., M, solution adjoint equation, whereinAnd φmAdjoint equation such as (3) formula:
Wherein, OmIt is desired output image, p, q belongs to the part of { (0,0), (0,1), (0,2), (1,1), (2,0) }
The indexed set of change.
Further, the derivative of translation rotational invariants during formula (4) calculating j=0,1 ..., 16 is utilized;Pass through public affairs
Formula (4) calculates j=0,1 ..., the derivative of the translation rotational invariants of control function when 16WithIn adjoint equation
With the help of, for ajAnd b (t)j(t), J derivative is as follows in each iteration:
Wherein, adjoint functionAnd φmIt is the answer of equation (3).
Further, golden section search is performed along the direction of search in the step S7, and constantly updates system function,
Carry out next circulation;
Golden section search is performed along the direction of search, control function a is updatedjAnd b (t)j(t), j=0,1 ..., 16, after
Continue next circulation, until j=16, be trained.
Further, in the step S8, circulation, output system function are terminatedWith
Compared with prior art, the beneficial effect of technical solution of the present invention is:
The present invention in the processing intrinsic Composition Estimation of image independent of the prior-constrained of determination, and using data-driven
Mode builds partial differential equation;The direction of search is determined using conjugate gradient method, relative to steepest descent method and Newton method, the party
Conjugacy is combined by method with steepest descent method;One group of conjugate direction is constructed using the gradient at known point, and along this prescription
To scanning for, the minimal point of object function is obtained, to determine the optimal direction of search;This method can effectively realize different illumination
The image of condition carries out the reflex components and shade composition that eigen decomposition obtains the image.
Brief description of the drawings
Fig. 1 is flow chart of the method for the present invention.
Embodiment
Accompanying drawing being given for example only property explanation, it is impossible to be interpreted as the limitation to this patent;
In order to more preferably illustrate the present embodiment, some parts of accompanying drawing have omission, zoomed in or out, and do not represent actual product
Size;
To those skilled in the art, it is to be appreciated that some known features and its explanation, which may be omitted, in accompanying drawing
's.
Technical scheme is described further with reference to the accompanying drawings and examples.
Embodiment 1
As shown in figure 1, a kind of intrinsic image decomposition method learnt based on partial differential equation, is comprised the steps of:
1), the training stage:
Step 1:Input training data pair;
The training data of input to being MIT-Intrinsic Images storehouses in image, training input data include 20 classes
Picture, the picture of the different illumination of 10 Same Scenes is included per class, totally 220 pictures, and training output data is correspondence input figure
The shade composition that piece is gathered.
Step 2:Control function is initialized, because object function is non-convex, the convergence direction for minimizing process is depended on
The local minimum of initialization, step 3 is turned to when initialization procedure is not restrained, circulation is performed;Otherwise step 8 is turned to;
By solving (1) formula to control function aj(t), t=0, Δ t, 1- Δs t is initialized, and is now fixed
bj(t), j=0,1,16,
Fa (t)=d (1)
Wherein, u is output image, and v is indicator function, and it is the extensive letter searched in image to be introduced into the purpose of indicator function
Breath, in order to correctly instruct u differentiation.Because object function is non-convex, the convergence for minimizing process is first towards depending on
The local minimum of beginningization, when control function is not restrained, turns to step 3, performs circulation, otherwise turn to step 8.
Step 3:Solve the optimal control equation constrained with PDE;
By introducing the adjoint equation of (2) formula, the gateaux derivative of j-th of translation rotational invariants is calculated, therefore part is most
The figure of meritWithIt can be calculated and obtained by the algorithm based on gradient, for umAnd vm, work as m=1,2 ..., M
When solve (2) formula:
Wherein Ω is the rectangular area occupied by input picture, and T is that PDE systems complete Vision information processing and output result
Time, and be u and v initial function respectively.For computational problem and it is related to deduction mathematically, the present invention will be around it
Fill the null value of several pixel wides.Due to chronomere can be changed, T=1, and it is fixed as being defined on one group of letter on Q
Number, is respectively intended to control u and v differentiation.
Step 4:Solve adjoint equation when adjoint function takes particular value;
ForAnd φm, m=1, during 2 ..., M, solution adjoint equation, whereinAnd φmAdjoint equation such as (3) formula:
Step 5:J=0,1 is calculated using formula (4), the derivative of translation rotational invariants when 16;
J=0,1 is calculated by formula (4), the derivative of the translation rotational invariants of control function when 16
WithWith the help of adjoint equation, for ajAnd b (t)j(t), J derivative is as follows in each iteration:
Wherein, λjAnd ujIt is positive weighting parameters, adjoint functionAnd φmIt is the answer of equation (3).
Step 6:The direction of search is determined using conjugate gradient method;
Determine the direction of search using conjugate gradient method, wherein conjugate gradient method is mutually to tie conjugacy with steepest descent method
Close.One group of conjugate direction is constructed using the gradient at known point, and along this prescription to scanning for, obtains the pole of object function
Dot, to determine the optimal direction of search.
Step 7:Golden section search is performed along the direction of search, and constantly updates system function, next circulation is carried out;
Golden section search is performed along the direction of search, system function a is updatedjAnd b (t)j(t), j=0,1 ..., 16, after
Continue next circulation, until j=16, be trained;
Step 8:Terminate circulation, output system function.
2), the application stage:
Step 9:Prepare application data, data picture feature is background black, and object is single and prominent;
Step 10:To give data carry out eigen decomposition application, obtain the reflex components and shade composition of image.
The same or analogous part of same or analogous label correspondence;
Position relationship is used for being given for example only property explanation described in accompanying drawing, it is impossible to be interpreted as the limitation to this patent;
Obviously, the above embodiment of the present invention is only intended to clearly illustrate example of the present invention, and is not pair
The restriction of embodiments of the present invention.For those of ordinary skill in the field, may be used also on the basis of the above description
To make other changes in different forms.There is no necessity and possibility to exhaust all the enbodiments.It is all this
Any modifications, equivalent substitutions and improvements made within the spirit and principle of invention etc., should be included in the claims in the present invention
Protection domain within.
Claims (7)
1. a kind of intrinsic image decomposition method learnt based on partial differential equation, it is characterised in that comprise the following steps:
S1:Input training data pair;
S2:Control function is initialized, because object function is non-convex, the convergence of process is minimized towards depending on initialization
Local minimum, step S3 is turned to when initialization procedure is not restrained, circulation is performed;Otherwise step S8 is turned to;
S3:Solve the optimal control equation constrained with PDE;
S4:Solve adjoint equation when adjoint function takes particular value;
S5:J=0,1 is calculated using following formula ..., the derivative of translation rotational invariants when 16:
Wherein, J is translation rotational invariants, and j is the number for translating rotational invariants, aj, bjIt is control function, λjAnd ujIt is positive
Weighting parameters,And φmIt is adjoint function, u is output image, and Ω is the rectangular area occupied by input picture, fΩIt is Ω
Initial function, m=1,2 ..., M, M for input data sample logarithm, inv (u, v) represent matrix (u, v) is inverted, v refers to
Show function, it is to collect the extensive information in image to be introduced into the purpose of indicator function, to instruct u differentiation.
S6:The direction of search is determined using conjugate gradient method;
S7:Golden section search is performed along the direction of search, and constantly updates system function, next circulation is carried out, until j=
16, it is trained;
S8:Terminate circulation, output system function;
S9:Prepare application data, data picture feature is background black, and object is single and prominent;
S10:Using obtained system function, to give data to carry out eigen decomposition application, obtain the reflex components and the moon of image
Shadow composition.
2. the intrinsic image decomposition method according to claim 1 learnt based on partial differential equation, it is characterised in that described
By solving (1) formula to control function a in step S2j(t), t=0, Δ t ..., 1- Δ t is initialized, and now fixes bj
(t), j=0,1 ..., 16,
Fa (t)=d (1)
<mrow>
<mi>F</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msub>
<mo>&Integral;</mo>
<mi>&Omega;</mi>
</msub>
<msub>
<mi>f</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>u</mi>
<mi>m</mi>
</msub>
<mo>,</mo>
<msub>
<mi>v</mi>
<mi>m</mi>
</msub>
<mo>)</mo>
</mrow>
<msubsup>
<mi>f</mi>
<mi>m</mi>
<mi>T</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>u</mi>
<mi>m</mi>
</msub>
<mo>,</mo>
<msub>
<mi>v</mi>
<mi>m</mi>
</msub>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>&Omega;</mi>
</mrow>
<mrow>
<mi>d</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msub>
<mo>&Integral;</mo>
<mi>&Omega;</mi>
</msub>
<msub>
<mi>f</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>u</mi>
<mi>m</mi>
</msub>
<mo>,</mo>
<msub>
<mi>v</mi>
<mi>m</mi>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>d</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>&Omega;</mi>
</mrow>
Wherein, ajAnd b (t)j(t) it is control function, a (t)={ aj} and b (t)={ b (t)j(t) } it is defined in the collection of functions on Q
Close, be respectively intended to control u and v differentiation, fuAnd fvIt is u and v initial function respectively.
3. the intrinsic image decomposition method according to claim 2 learnt based on partial differential equation, it is characterised in that described
By introducing the adjoint equation of (2) formula in step S3, the gateaux derivative of j-th of translation rotational invariants is calculated, therefore part is most
The figure of meritWithCalculated and obtained by the algorithm based on gradient, for umAnd vm, work as m=1, solved during 2 ..., M
(2) formula:
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>u</mi>
<mi>m</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<msub>
<mi>F</mi>
<mi>u</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>u</mi>
<mi>m</mi>
</msub>
<mo>,</mo>
<msub>
<mi>v</mi>
<mi>m</mi>
</msub>
<mo>,</mo>
<msubsup>
<mrow>
<mo>{</mo>
<msub>
<mi>a</mi>
<mi>j</mi>
</msub>
<mo>}</mo>
</mrow>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mn>16</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
<mo>&Element;</mo>
<mi>Q</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
<mo>&Element;</mo>
<mi>&Gamma;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mi>m</mi>
</msub>
<msub>
<mo>|</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>f</mi>
<msub>
<mi>u</mi>
<mi>m</mi>
</msub>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
<mo>&Element;</mo>
<mi>&Omega;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>u</mi>
<mi>m</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<msub>
<mi>F</mi>
<mi>v</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>u</mi>
<mi>m</mi>
</msub>
<mo>,</mo>
<msub>
<mi>v</mi>
<mi>m</mi>
</msub>
<mo>,</mo>
<msubsup>
<mrow>
<mo>{</mo>
<msub>
<mi>a</mi>
<mi>j</mi>
</msub>
<mo>}</mo>
</mrow>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mn>16</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
<mo>&Element;</mo>
<mi>Q</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>v</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
<mo>&Element;</mo>
<mi>&Gamma;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>v</mi>
<mi>m</mi>
</msub>
<msub>
<mo>|</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>f</mi>
<msub>
<mi>v</mi>
<mi>m</mi>
</msub>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
<mo>&Element;</mo>
<mi>&Omega;</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein, T is the time that PDE systems complete Vision information processing and output result, and Q is Ω × (0, T), and Γ is
4. the intrinsic image decomposition method according to claim 2 learnt based on partial differential equation, it is characterised in that described
Adjoint equation during adjoint function particular value is solved in step S4:
ForAnd φm, m=1, during 2 ..., M, solution adjoint equation, whereinAnd φmAdjoint equation such as (3) formula:
<mrow>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>F</mi>
<mi>u</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>u</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mn>16</mn>
</msubsup>
<msub>
<mi>a</mi>
<mi>j</mi>
</msub>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>inv</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>,</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>u</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>F</mi>
<mi>v</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>u</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mn>16</mn>
</msubsup>
<msub>
<mi>b</mi>
<mi>j</mi>
</msub>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>inv</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mo>,</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>u</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mover>
<mi>&sigma;</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>F</mi>
<mi>u</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>u</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mn>16</mn>
</msubsup>
<msub>
<mi>a</mi>
<mi>j</mi>
</msub>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>inv</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>,</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>u</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mover>
<mi>&sigma;</mi>
<mo>~</mo>
</mover>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>F</mi>
<mi>v</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>u</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mn>16</mn>
</msubsup>
<msub>
<mi>b</mi>
<mi>j</mi>
</msub>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>inv</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mo>,</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>u</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>u</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mrow>
<mi>p</mi>
<mo>+</mo>
<mi>q</mi>
</mrow>
</msup>
<mi>u</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>x</mi>
<mi>p</mi>
</msup>
<mo>&part;</mo>
<msup>
<mi>y</mi>
<mi>q</mi>
</msup>
</mrow>
</mfrac>
<mo>,</mo>
<msub>
<mi>v</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mrow>
<mi>p</mi>
<mo>+</mo>
<mi>q</mi>
</mrow>
</msup>
<mi>v</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>x</mi>
<mi>p</mi>
</msup>
<mo>&part;</mo>
<msup>
<mi>y</mi>
<mi>q</mi>
</msup>
</mrow>
</mfrac>
</mrow>
Wherein, OmIt is desired output image, p, q belongs to the localized variation of { (0,0), (0,1), (0,2), (1,1), (2,0) }
Indexed set.
5. the intrinsic image decomposition method according to claim 4 learnt based on partial differential equation, it is characterised in that utilize
Formula (4) calculates j=0,1 ..., the derivative of translation rotational invariants when 16;J=0,1 is calculated by formula (4) ...,
The derivative of the translation rotational invariants of control function when 16WithWith the help of adjoint equation, for ajAnd b (t)j
(t), J derivative is as follows in each iteration:
Wherein, adjoint functionAnd φmIt is the answer of equation (3).
6. the intrinsic image decomposition method according to claim 5 learnt based on partial differential equation, it is characterised in that described
Golden section search is performed along the direction of search in step S7, and constantly updates system function, next circulation is carried out;
Golden section search is performed along the direction of search, control function a is updatedjAnd b (t)j(t), j=0,1 ..., 16, under continuation
One circulation, until j=16, is trained.
7. the intrinsic image decomposition method according to claim 6 learnt based on partial differential equation, it is characterised in that described
In step S8, circulation, output system function are terminatedWith
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710138216.8A CN107133953B (en) | 2017-03-09 | 2017-03-09 | Intrinsic image decomposition method based on partial differential equation learning |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710138216.8A CN107133953B (en) | 2017-03-09 | 2017-03-09 | Intrinsic image decomposition method based on partial differential equation learning |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107133953A true CN107133953A (en) | 2017-09-05 |
CN107133953B CN107133953B (en) | 2020-04-24 |
Family
ID=59720752
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710138216.8A Active CN107133953B (en) | 2017-03-09 | 2017-03-09 | Intrinsic image decomposition method based on partial differential equation learning |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107133953B (en) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104700109A (en) * | 2015-03-24 | 2015-06-10 | 清华大学 | Method and device for decomposing hyper-spectral intrinsic images |
US20160171755A1 (en) * | 2013-08-08 | 2016-06-16 | Adobe Systems Incorporated | Automatic geometry and lighting inference for realistic image editing |
CN106355601A (en) * | 2016-08-31 | 2017-01-25 | 上海交通大学 | Intrinsic image decomposition method |
-
2017
- 2017-03-09 CN CN201710138216.8A patent/CN107133953B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160171755A1 (en) * | 2013-08-08 | 2016-06-16 | Adobe Systems Incorporated | Automatic geometry and lighting inference for realistic image editing |
CN104700109A (en) * | 2015-03-24 | 2015-06-10 | 清华大学 | Method and device for decomposing hyper-spectral intrinsic images |
CN106355601A (en) * | 2016-08-31 | 2017-01-25 | 上海交通大学 | Intrinsic image decomposition method |
Non-Patent Citations (1)
Title |
---|
胡可鑫 等: "基于先验知识的快速人脸光照迁移算法", 《计算机辅助设计与图形学学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN107133953B (en) | 2020-04-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Lee et al. | Single-image depth estimation based on fourier domain analysis | |
Gupta et al. | Texas 3D face recognition database | |
Zhang et al. | NLDN: Non-local dehazing network for dense haze removal | |
CN105069746B (en) | Video real-time face replacement method and its system based on local affine invariant and color transfer technology | |
US20200265597A1 (en) | Method for estimating high-quality depth maps based on depth prediction and enhancement subnetworks | |
CN104318569B (en) | Space salient region extraction method based on depth variation model | |
CN104408701B (en) | A kind of large scene video image joining method | |
Chai et al. | Multifocus image fusion scheme using focused region detection and multiresolution | |
CN107463954B (en) | A kind of template matching recognition methods obscuring different spectrogram picture | |
CN109376611A (en) | A kind of saliency detection method based on 3D convolutional neural networks | |
Legg et al. | Feature neighbourhood mutual information for multi-modal image registration: an application to eye fundus imaging | |
CN106886986B (en) | Image interfusion method based on adaptive group structure sparse dictionary study | |
Wei et al. | Learning flow-based feature warping for face frontalization with illumination inconsistent supervision | |
CN106952341A (en) | The underwater scene three-dimensional point cloud method for reconstructing and its system of a kind of view-based access control model | |
CN103020933B (en) | A kind of multisource image anastomosing method based on bionic visual mechanism | |
CN111882489A (en) | Super-resolution graph recovery method for simultaneously enhancing underwater images | |
CN105976395B (en) | A kind of video target tracking method based on rarefaction representation | |
CN106920214A (en) | Spatial target images super resolution ratio reconstruction method | |
CN107194912A (en) | The brain CT/MR image interfusion methods of improvement coupling dictionary learning based on rarefaction representation | |
CN108921783A (en) | A kind of satellite image super resolution ratio reconstruction method based on losses by mixture function constraint | |
CN109493309A (en) | A kind of infrared and visible images variation fusion method keeping conspicuousness information | |
Chen et al. | Blood vessel enhancement via multi-dictionary and sparse coding: Application to retinal vessel enhancing | |
CN102034247A (en) | Motion capture method for binocular vision image based on background modeling | |
Jing et al. | Cloud removal for optical remote sensing imagery using the SPA-CycleGAN network | |
CN109218706A (en) | A method of 3 D visual image is generated by single image |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |