For given data of dimension , , (usually ) this function computes the proximal map of the second order TV prior in a cyclic manner.
For the complete functional, see The second order total variation functional from the mid point model
For the cyclic proximal point algorithm we employ the same idea as for the proximal map(s) of TV just splitting with respect to three for the second order terms and in even and odd for both dimensions (leading to four inner sums) and employing hence 7 proximal maps in a cycle and apply the second order difference prox and the second order mixed difference prox, respectively.
This function can be extended by carrying a different weight
for a difference in dimension and can also be extended
to diagonal differences when is introduced, i.e. the weights are given as a(n upper triangular) matrix. The standard
is, to treat this matrix as a constant value given by
lambda. A vector for
lambda is interpreted as a diagonal matrix.
Finally the distance term can be replaced by using the optional
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 % proxTV(M,x,lambda) compute all proxima of second order finite differences % with a splitting modulo 3 for the mid point approach within the TV2 % term and a 2-2 splitting in the mixed terms in a cyclic manner % % INPUT % M : a manifold % x : data (manifold-valued) % lambda : parameter of the prox % % OPTIONAL % 'FixedMask' : binary mask the size of data items in x to indicate fixed % data items % 'SecDiffProx' : (@(x1,x2,x3,lambda) % proxAbsoluteSecondOrderDifference(M,x1,x2,x3,lambda)) % specify a prox for (mod 3) TV2 term proxes % 'SecDiffMixProx' : (@(x1,x2,x3,x4,lambda) % proxAbsoluteSecondOrderMixedDifference(M,x1,x2,x3,x4,lambda)) % specify a prox for (mod 2,2) TV2 mixed term proxes % % OUTPUT % y : result of applying the proximal maps in a cyclic manner. % --- % Manifold-valued Image Restoration Toolbox 1.2 | R. Bergmann | 2017-12-11
- The Cyclic Proximal Point Algorithm
- The proximal map of the second order difference
- The proximal map of second order mixed difference