We investigate whether it is possible to use the bidomain model and body surface potential maps (BSPMs) to compute the size and position of ischemic regions in the human heart. This leads to a severely ill posed inverse problem for a potential equation. We do not use the classical inverse problems of electrocardiography, in which the unknown sources are the epicardial potential distribution or the activation sequence. Instead we employ the bidomain theory to obtain a model that also enables identification of ischemic regions transmurally. This approach makes it possible to distinguish between subendocardial and transmural cases, only using the BSPM data. The main focus is on testing a previously published algorithm on clinical data, and the results are compared with images taken with perfusion scintigraphy. For the four patients involved in this study, the two modalities produce results that are rather similar: The relative differences between the center of mass and the size of the ischemic regions, suggested by the two modalities, are 10.8% ± 4.4% and 7.1% ± 4.6%, respectively. We also present some simulations which indicate that the methodology is robust with respect to uncertainties in important model parameters. However, in contrast to what has been observed in investigations only involving synthetic data, inequality constraints are needed to obtain sound results.

In this paper, we propose a PDE-based level set method. Traditionally, interfaces are represented by the zero level set of continuous level set functions. Instead, we let the interfaces be represented by discontinuities of piecewise constant level set functions. Each level set function can at convergence only take two values, i.e., it can only be 1 or -1; thus, our method is related to phase-field methods. Some of the properties of standard level set methods are preserved in the proposed method, while others are not. Using this new method for interface problems, we need to minimize a smooth convex functional under a quadratic constraint. The level set functions are discontinuous at convergence, but the minimization functional is smooth. We show numerical results using the method for segmentation of digital images.

In this work, we use partial differential equation techniques to remove noise from digital images. The removal is done in two steps. We first use a total-variation filter to smooth the normal vectors of the level curves of a noise image. After this, we try to find a surface to fit the smoothed normal vectors. For each of these two stages, the problem is reduced to a nonlinear partial differential equation. Finite difference schemes are used to solve these equations. A broad range of numerical examples are given in the paper.

In this paper, we introduce a new method for image smoothing based on a fourth-order PDE model. The method is tested on a broad range of real medical magnetic resonance images, both in space and time, as well as on nonmedical synthesized test images. Our algorithm demonstrates good noise suppression without destruction of important anatomical or functional detail, even at poor signal-to-noise ratio. We have also compared our method with related PDE models.