Bayesian mapping of brain regions using compound Markov random field priors

Human brain mapping, i.e. the detection of functional regions and their connections, has experienced enormous progress through the use of functional magnetic resonance imaging (fMRI). The massive spatio-temporal data sets generated by this imaging technique impose challenging problems for statistical analysis. Many approaches focus on adequate modeling of the temporal component. Spatial aspects are often considered only in a separate postprocessing step, if at all, or modeling is based on Gaussian random fields. A weakness of Gaussian spatial smoothing is possible underestimation of activation peaks or blurring of sharp transitions between activated and non-activated regions. In this paper we suggest Bayesian spatio-temporal models, where spatial adaptivity is improved through inhomogeneous or compound Markov random field priors. Inference is based on an approximate MCMC technique. Performance of our approach is investigated through a simulation study, including a comparison to models based on Gaussian as well as more robust spatial priors in terms of pixelwise and global MSEs. Finally we demonstrate its use by an application to fMRI data from a visual stimulation experiment for assessing activation in visual cortical areas

Spatial priors for tomographic reconstructions from limited data

Tomografie is het reconstrueren van het inwendige van een object a.d.h.v externe metingen, b.v. beelden verkregen met X-stralen of microgolven. Deze thesis bekijkt de specifieke aspecten van microgolftomografie en magnetische resonantie beeldvorming (Magnetic Resonance Imaging – MRI); beide technieken zijn onschadelijk voor de mens. Terwijl het gebruik van MRI wijdverspreid is voor veel klinische toepassingen, is microgolftomografie nog niet in klinisch gebruik ondanks zijn potentiële voordelen. Door de lage kost en draagbaarheid van de toestellen is het een waardevolle aanvulling aan het assortiment

Bayesian Inference for Inverse Problems

Inverse problems arise everywhere we have indirect measurement. Regularization and Bayesian inference methods are two main approaches to handle inverse problems. Bayesian inference approach is more general and has much more tools for developing efficient methods for difficult problems. In this chapter, first, an overview of the Bayesian parameter estimation is presented, then we see the extension for inverse problems. The main difficulty is the great dimension of unknown quantity and the appropriate choice of the prior law. The second main difficulty is the computational aspects. Different approximate Bayesian computations and in particular the variational Bayesian approximation (VBA) methods are explained in details

Probabilistic partial volume modelling of biomedical tomographic image data

EThOS - Electronic Theses Online ServiceGBUnited Kingdo

Inverse Problems in data-driven multi-scale Systems Medicine: application to cancer physiology

Systems Medicine is an interdisciplinary framework involving reciprocal feedback between clinical investigation and mathematical modeling/analysis. Its aim is to improve the understanding of complex diseases by integrating knowledge and data across multiple levels of biological organization. This Thesis focuses on three inverse problems, arising from three kinds of data and related to cancer physiology, at different scales: tissues, cells, molecules. The general assumption of this piece of research is that cancer is associated toa path ological glucose consumption and, in fact, its functional behavior can be assessed by nuclear medicine experiments using [18F]-fluorodeoxyglucose (FDG) as a radioactive tracer mimicking the glucose properties. At tissue-scale, this Thesis considers the Positron Emission Tomography (PET) imaging technique, and deals with two distinct issues within compartmental analysis. First, this Thesis presents a compartmental approach, referred to as reference tissue model, for the estimation of FDG kinetics inside cancer tissues when the arterial blood input of the system is unknown. Then, this Thesis proposes an efficient and reliable method for recovering the compartmental kinetic parameters for each PET image pixel in the context of parametric imaging, exploiting information on the tissue physiology. Standard models in compartmental analysis assume that phosphorylation and dephosphorylation of FDG occur in the same intracellular cytosolic volume. Advances in cell biochemistry have shown that the appropriate location of dephosphorylation is the endoplasmic reticulum (ER). Therefore, at cell-scale, this Thesis formalizes a biochemically-driven compartmental model accounting for the specific role played by the ER, and applies it to the analysis of in vitro experiments on FDG uptake by cancer cell cultures obtained with a LigandTracer (LT) device. Finally, at molecule-scale, this Thesis provides a preliminary mathematical investigation of a chemical reaction network (CRN), represented by a huge Molecular Interaction Map (MIM), describing the biochemical interactions occurring between signaling proteins in specific pathways within a cancer cell. The main issue addressed in this case is the network parameterization problem, i.e. how to determine the reaction rate coefficients from protein concentration data

Fusion and Analysis of Multidimensional Medical Image Data

Analýza medicínských obrazů je předmětem základního výzkumu již řadu let. Za tu dobu bylo v této oblasti publikováno mnoho výzkumných prací zabývajících se dílčími částmi jako je rekonstrukce obrazů, restaurace, segmentace, klasifikace, registrace (lícování) a fúze. Kromě obecného úvodu, pojednává tato disertační práce o dvou medicínsky orientovaných tématech, jež byla formulována ve spolupráci s Philips Netherland BV, divizí Philips Healthcare. První téma je zaměřeno na oblast zpracování obrazů subtrakční angiografie dolních končetin člověka získaných pomocí výpočetní X-Ray tomografie (CT). Subtrakční angiografie je obvykle využívaná při podezření na periferní cévní onemocnění (PAOD) nebo při akutním poškození dolních končetin jako jsou fraktury apod. Současné komerční metody nejsou dostatečně spolehlivé už v předzpracování, jako je například odstranění pacientského stolu, pokrývky, dlahy, apod. Spolehlivost a přesnost identifikace cév v subtrahovaných datech vedoucích v blízkosti kostí je v důsledku Partial Volume artefaktu rovněž nízká. Automatické odstranění kalcifikací nebo detekce malých cév doplňujících nezbytnou informaci o náhradním zásobení dolních končetin krví v případě přerušení hlavních zásobujících cév v současné době rovněž nesplňují kritéria pro plně automatické zpracování. Proto hlavním cílem týkající se tohoto tématu bylo vyvinout automatický systém, který by mohl současné nedostatky v CTSA vyšetření odstranit. Druhé téma je orientováno na identifikaci patologických změn na páteři člověka v CT obrazech se zaměřením na osteolytické a osteoblastické léze u jednotlivých obratlů. Tyto změny obvykle nastávají v důsledků postižení metastazujícím procesem rakovinového onemocnění. Pro detekci patologických změn je pak potřeba identifikace a segmentace jednotlivých obratlů. Přesnost analýzy jednotlivých lézí však závisí rovněž na správné identifikaci těla a zadních segmentů u jednotlivých obratlů a na segmentaci trabekulárního centra obratlů, tj. odstranění kortikální kosti. Během léčby mohou být pacienti skenováni vícekrát, obvykle s několika-mesíčním odstupem. Hodnocení případného vývoje již detekovaných patologických změn pak logicky vychází ze správné detekce patologií v jednotlivých obratlech korespondujících si v jednotlivých akvizicích. Jelikož jsou příslušné obratle v jednotlivých akvizicích obvykle na různé pozici, jejich fúze, vedoucí k analýze časového vývoje detekovaných patologií, je komplikovaná. Požadovaným výsledkem v tomto tématu je vytvoření komplexního systému pro detekci patologických změn v páteři, především osteoblastických a osteolytických lézí. Takový systém tedy musí umožnovat jak segmentaci jednotlivých obratlů, jejich automatické rozdělení na hlavní části a odstranění kortikální kosti, tak také detekci patologických změn a jejich hodnocení. Ačkoliv je tato disertační práce v obou výše zmíněných tématech primárně zaměřena na experimentální část zpracování medicínských obrazů, zabývá se všemi nezbytnými kroky, jako je předzpracování, registrace, dodatečné zpracování a hodnocení výsledků, vedoucími k možné aplikovatelnosti obou systému v klinické praxi. Jelikož oba systémy byly řešeny v rámci týmové spolupráce jako celek, u obou témat jsou pro některé konkrétní kroky uvedeny odkazy na doktorskou práci Miloše Malínského.Analysis of medical images has been subject of basic research for many years. Many research papers have been published in the field related to image analysis and focused on partial aspects such as reconstruction, restoration, segmentation and classification, registration (spatial alignment) and fusion. Besides the introduction of related general concepts used in medical image processing, this thesis deals with two specific medical problems formulated in cooperation with Philips Netherland BV, Philips Healthcare division. The first topic is focused on subtraction angiography in patients’ lower legs utilizing image data from X-Ray computed tomography (CT). CT subtraction angiography (CTSA) is typically used for indication of the Peripheral Artery Occlusive Disease (PAOD) and for examination of acute injuries of lower legs such as acute fractures, etc. Current methods in clinical praxis are not sufficient regarding the pre-processing such as masking of patient desk, cover, splint, etc. The subtraction of blood vessels adjacent to neighboring bones in lower legs is of low accuracy due to the Partial Volume artifact. Masking of calcifications and detection of tiny blood vessels complementing necessary information about the alternative blood supply in lower legs in case of obstruction in main arteries is also not reliable for fully automated process presently. Therefore, the main aim regarding this topic was to develop an automated framework that could overcome current shortcomings in CTSA examination. The second topic is oriented on the identification and evaluation of pathologic changes in human spine, focusing on osteolytic and osteoblastic lesions in individual vertebrae in CT images. Such changes occur typically as a consequence of metastasizing process of cancerous disease. For the detection of pathologic changes, an identification and segmentation of individual vertebrae is necessary. Moreover, the analysis of individual lesions in vertebrae depends also on correct identification of vertebral body and posterior segments of each vertebra, and on segmentation of their trabecular centers. Patients are typically examined more than once during their therapy. Then, the evaluation of possible tumorous progression is based on accurate detection of pathologies in individual vertebrae in the base-line and corresponding follow-up images. Since the corresponding vertebrae are in mutually different positions in the follow-up images, their fusion leading to the analysis of the lesion progression is complicated. The main aim regarding this topic is to develop a complex framework for detection of pathologic lesions on spine, with the main focus on osteoblastic and osteolystic lesions. Such system has to provide not only reliable segmentation of individual vertebrae and detection of their main regions but also the masking of their cortical bone, detection of their pathologic changes and their evaluation. Although this dissertation thesis is primarily oriented at the experimental part of medical image processing considering both the above mentioned topics, it deals with all necessary processing steps, i.e. preprocessing, image registration, post-processing and evaluation of results, leading to the future use of both frameworks in clinical practice. Since both frameworks were developed in a team, there are some chapters referring to the dissertation thesis of Milos Malinsky.

Advanced Image Acquisition, Processing Techniques and Applications

"Advanced Image Acquisition, Processing Techniques and Applications" is the first book of a series that provides image processing principles and practical software implementation on a broad range of applications. The book integrates material from leading researchers on Applied Digital Image Acquisition and Processing. An important feature of the book is its emphasis on software tools and scientific computing in order to enhance results and arrive at problem solution

Enhancement of noisy planar nuclear medicine images using mean field annealing

Abstract Nuclear Medicine (NM) images inherently suffer from large amounts of noise and blur. The purpose of this research is to reduce the noise and blur while maintaining image integrity for improved diagnosis. The proposal is to further improve image quality after the standard pre- and post-processing undertaken by a gamma camera system. Mean Field Annealing (MFA), the image processing technique used in this research is a well known image processing approach. The MFA algorithm uses two techniques to achieve image restoration. Gradient descent is used as the minimisation technique, while a deterministic approximation to Simulated Annealing (SA) is used for optimisation. The algorithm anisotropically diffuses an image, iteratively smoothing regions that are considered non-edges and still preserving edge integrity until a global minimum is obtained. A known advantage of MFA is that it is able to minimise to this global minimum, skipping over local minima while still providing comparable results to SA with significantly less computational effort. Image blur is measured using either a point or line source. Both allow for the derivation of a Point Spread Function (PSF) that is used to de-blur the image. The noise variance can be measured using a flood source. The noise is due to the random fluctuations in the environment as well as other contributors. Noisy blurred NM images can be difficult to diagnose particularly at regions with steep intensity gradients and for this reason MFA is considered suitable for image restoration. From the literature it is evident that MFA can be applied successfully to digital phantom images providing improved performance over Wiener filters. In this paper MFA is shown to yield image enhancement of planar NM images by implementing a sharpening filter as a post MFA processing technique