Overall Objectives
Research Program
Application Domains
Highlights of the Year
New Software and Platforms
New Results
Bilateral Contracts and Grants with Industry
Partnerships and Cooperations
XML PDF e-pub
PDF e-Pub

Section: Research Program

Discrete Biomedical Image Perception

A wide variety of tasks in medical image analysis can be formulated as discrete labeling problems. In very simple terms, a discrete optimization problem can be stated as follows: we are given a discrete set of variables 𝒱, all of which are vertices in a graph 𝒢. The edges of this graph (denoted by ℰ) encode the variables' relationships. We are also given as input a discrete set of labels ℒ. We must then assign one label from ℒ to each variable in 𝒱. However, each time we choose to assign a label, say, xp1 to a variable p1, we are forced to pay a price according to the so-called singleton potential function gp(xp), while each time we choose to assign a pair of labels, say, xp1 and xp2 to two interrelated variables p1 and p2 (two nodes that are connected by an edge in the graph 𝒢), we are also forced to pay another price, which is now determined by the so called pairwise potential function fp1p2(xp1,xp2). Both the singleton and pairwise potential functions are problem specific and are thus assumed to be provided as input.

Our goal is then to choose a labeling which will allow us to pay the smallest total price. In other words, based on what we have mentioned above, we want to choose a labeling that minimizes the sum of all the MRF potentials, or equivalently the MRF energy. This amounts to solving the following optimization problem:

The use of such a model can describe a number of challenging problems in medical image analysis. However these simplistic models can only account for simple interactions between variables, a rather constrained scenario for high-level medical imaging perception tasks. One can augment the expression power of this model through higher order interactions between variables, or a number of cliques {Ci,i∈[1,n]={{pi1,⋯,pi|Ci|}} of order |Ci| that will augment the definition of 𝒱 and will introduce hyper-vertices:

where fp1⋯pn is the price to pay for associating the labels (xpi1,⋯,pxi|Ci|) to the nodes (p1⋯pi|Ci|). Parameter inference, addressed by minimizing the problem above, is the most critical aspect in computational medicine and efficient optimization algorithms are to be evaluated both in terms of computational complexity as well as of inference performance. State of the art methods include deterministic and non-deterministic annealing, genetic algorithms, max-flow/min-cut techniques and relaxation. These methods offer certain strengths while exhibiting certain limitations, mostly related to the amount of interactions which can be tolerated among neighborhood nodes. In the area of medical imaging where domain knowledge is quite strong, one would expect that such interactions should be enforced at the largest scale possible.