A well-posedness framework for inpainting based on coherence transport

Image inpainting is the process of touching-up damaged or unwanted portions of a picture and is an important task in image processing. For this purpose Bornemann and März [J. Math. Imaging Vis. , 28 (2007), pp. 259– 278] introduced a very efficient method called Image Inpainting Based on Coherence Transport which fills the missing region by advecting the image information along integral curves of a coherence vector field from the boundary towards the interior of the hole. The mathematical model behind this method is a first-order functional advection PDE posed on a compact domain with all inflow boundary. We show that this problem is well-posed under certain conditions

Inpainting of Cyclic Data using First and Second Order Differences

Cyclic data arise in various image and signal processing applications such as interferometric synthetic aperture radar, electroencephalogram data analysis, and color image restoration in HSV or LCh spaces. In this paper we introduce a variational inpainting model for cyclic data which utilizes our definition of absolute cyclic second order differences. Based on analytical expressions for the proximal mappings of these differences we propose a cyclic proximal point algorithm (CPPA) for minimizing the corresponding functional. We choose appropriate cycles to implement this algorithm in an efficient way. We further introduce a simple strategy to initialize the unknown inpainting region. Numerical results both for synthetic and real-world data demonstrate the performance of our algorithm.Comment: accepted Converence Paper at EMMCVPR'1

A Second Order TV-type Approach for Inpainting and Denoising Higher Dimensional Combined Cyclic and Vector Space Data

In this paper we consider denoising and inpainting problems for higher dimensional combined cyclic and linear space valued data. These kind of data appear when dealing with nonlinear color spaces such as HSV, and they can be obtained by changing the space domain of, e.g., an optical flow field to polar coordinates. For such nonlinear data spaces, we develop algorithms for the solution of the corresponding second order total variation (TV) type problems for denoising, inpainting as well as the combination of both. We provide a convergence analysis and we apply the algorithms to concrete problems.Comment: revised submitted versio

Population of human ventricular cell models calibrated with in vivo measurements unravels ionic mechanisms of cardiac alternans

Cardiac alternansis an important risk factor in cardiac physiology, and is related to the initiation of many pathophysiological conditions. However, the mechanisms underlying the generation of alternans remain unclear. In this study, we used a population of computational human ventricle models based onthe O’Hara model [1] to explore the effect of 11 key factors experimentally reported to be related to alternans. In vivo experimental datasets coming from patients undergoing cardiac surgery were used in the calibration of our in silico population of models. The calibrated models in the population were divided into two groups (Normal and Alternans) depending on alternans occurrence. Our results showed that there were significant differences in the following 5 ionic currents between the two groups: fast sodium current, sodium calcium exchanger current, sodium potassium pump current, sarcoplasmic reticulum (SR) calcium release flux and SR calcium reuptake flux. Further analysis indicated that fast sodium current and SR calcium uptake were the two most significant currents that contributed to voltage and calcium alternans generation, respectively

Computational probabilistic quantification of pro-arrhythmic risk from scar and left-to-right heterogeneity in the human ventricles

Both scar and left-to-right ventricular (LV/RV) differences in repolarization properties have been implicated as risk factors for lethal arrhythmias. As a possible mechanism for the initiation of re-entry, a recent study has indicated that LV/RV heterogeneities in action potential duration (APD) adaptation can cause a transient increase in APD dispersion following rate acceleration, promoting unidirectional block of conduction at the LV/RV junction. In the presence of an ischemic region and ectopic stimulation, a pathological dispersion in repolarization has been suggested to increase the risk of electrical re-entry. However, the exact location and timing of the ectopic activation play a crucial role in initiation of re-entry, and certain combinations may lead to re-entry even under normal LV/RV dispersion in repolarization. This suggests that the phenomenon needs to be investigated in a quantitative way. In this study we employ a computationally efficient, phenomenological model in order to investigate the proarrhythmic properties of a range of combinations of position and timing of an ectopic activation. This allows us to probabilistically study how increasing interventricular dispersion of repolarization increases arrhythmic risk. Results indicate that a larger LV/RV dispersion in repolarization allows ectopic beats to initiate re-entry during a significantly larger time window and from a greater number of locations compared to the case of smaller LV/RV dispersion

Guidefill: GPU accelerated, artist guided geometric inpainting for 3D conversion of film

The conversion of traditional film into stereo 3D has become an important problem in the past decade. One of the main bottlenecks is a disocclusion step, which in commercial 3D conversion is usually done by teams of artists armed with a toolbox of inpainting algorithms. A current difficulty in this is that most available algorithms either are too slow for interactive use or provide no intuitive means for users to tweak the output. In this paper we present a new fast inpainting algorithm based on transporting along automatically detected splines, which the user may edit. Our algorithm is implemented on the GPU and fills the inpainting domain in successive shells that adapt their shape on the y. In order to allocate GPU resources as efficiently as possible, we propose a parallel algorithm to track the inpainting interface as it evolves, ensuring that no resources are wasted on pixels that are not currently being worked on. Theoretical analyses of the time and processor complexity of our algorithm without and with tracking (as well as numerous numerical experiments) demonstrate the merits of the latter. Our transport mechanism is similar to the one used in coherence transport [F. Bornemann and T. März, J. Math. Imaging Vision, 28 (2007), pp. 259-278; T. März, SIAM J. Imaging Sci., 4 (2011), pp. 981-1000] but improves upon it by correcting a \kinking" phenomenon whereby extrapolated isophotes may bend at the boundary of the inpainting domain. Theoretical results explaining this phenomenon and its resolution are presented. Although our method ignores texture, in many cases this is not a problem due to the thin inpainting domains in 3D conversion. Experimental results show that our method can achieve a visual quality that is competitive with the state of the art while maintaining interactive speeds and providing the user with an intuitive interface to tweak the results.The work of the first author was supported by the Cambridge Commonwealth Trust and the Cambridge Center for Analysis. The work of the third author was supported by the Leverhulme Trust project Breaking the Nonconvexity Barrier, the EPSRC grants EP/M00483X/1 and EP/N014588/1, the Cantab Capital Institute for the Mathematics of Information, the CHiPS (Horizon 2020 RISE project grant), the Global Alliance project “Statistical and Mathematical Theory of Imaging,” and the Alan Turing Institute

On the eigenvalues and eigenvectors of nonsymmetric saddle point matrices preconditioned by block triangular matrices

Block lower triangular and block upper triangular matrices are popular preconditioners for nonsymmetric saddle point matrices. In this note we show that a block lower triangular preconditioner gives the same spectrum as a block upper triangular preconditioner and that the eigenvectors of the two preconditioned systems are related