3D Wavelet-based Fusion Techniques for Biomedical Imaging

Hoy en día las técnicas de adquisición de imágenes tridimensionales son comunes en diversas áreas, pero cabe destacar la relevancia que han adquirido en el ámbito de la imagen biomédica, dentro del cual encontramos una amplia gama de técnicas como la microscopía confocal, microscopía de dos fotones, microscopía de fluorescencia mediante lámina de luz, resonancia magnética nuclear, tomografía por emisión de positrones, tomografía de coherencia óptica, ecografía 3D y un largo etcétera. Un denominador común de todas esas aplicaciones es la constante necesidad por aumentar la resolución y la calidad de las imágenes adquiridas. En algunas de dichas técnicas de imagen tridimensional se da una interesante situación: aunque que cada volumen adquirido no contiene información suficiente para representar el objeto bajo estudio dentro de los parámetros de calidad requeridos por algunas aplicaciones finales, el esquema de adquisición permite la obtención de varios volúmenes que representan diferentes vistas de dicho objeto, de tal forma que cada una de las vistas proporciona información complementaria acerca del mismo. En este tipo de situación es posible, mediante la combinación de varias de esas vistas, obtener una mejor comprensión del objeto que a partir de cada una de ellas por separado. En el contexto de esta Tesis Doctoral se ha propuesto, desarrollado y validado una nueva metodología de proceso de imágenes basada en la transformada wavelet disc¬reta para la combinación, o fusión, de varias vistas con información complementaria de un mismo objeto. El método de fusión propuesto aprovecha la capacidad de descom¬posición en escalas y orientaciones de la transformada wavelet discreta para integrar en un solo volumen toda la información distribuida entre el conjunto de vistas adquiridas. El trabajo se centra en dos modalidades diferentes de imagen biomédica que per¬miten obtener tales adquisiciones multi-vista. La primera es una variante de la micro¬scopía de fluorescencia, la microscopía de fluorescencia mediante lámina de luz, que se utiliza para el estudio del desarrollo temprano de embriones vivos en diferentes modelos animales, como el pez cebra o el erizo de mar. La segunda modalidad es la resonancia magnética nuclear con realce tardío, que constituye una valiosa herramienta para evaluar la viabilidad del tejido miocárdico en pacientes con diversas miocardiopatías. Como parte de este trabajo, el método propuesto ha sido aplicado y validado en am¬bas modalidades de imagen. En el caso de la aplicación a microscopía de fluorescencia, los resultados de la fusión muestran un mejor contraste y nivel de detalle en comparación con cualquiera de las vistas individuales y el método no requiere de conocimiento previo acerca la función de dispersión puntual del sistema de imagen. Además, los resultados se han comparado con otros métodos existentes. Con respecto a la aplicación a imagen de resonancia magnética con realce tardío, los volúmenes fusionados resultantes pre-sentan una mejora cuantitativa en la nitidez de las estructuras relevantes y permiten una interpretación más sencilla y completa de la compleja estructura tridimensional del tejido miocárdico en pacientes con cardiopatía isquémica. Para ambas aplicaciones los resultados de esta tesis se encuentran actualmente en uso en los centros clínicos y de investigación con los que el autor ha colaborado durante este trabajo. Además se ha puesto a libre disposición de la comunidad científica la implementación del método de fusión propuesto. Por último, se ha tramitado también una solicitud de patente internacional que cubre el método de visualización desarrollado para la aplicación de Resonancia Magnética Nuclear. Abstract Nowadays three dimensional imaging techniques are common in several fields, but es-pecially in biomedical imaging, where we can find a wide range of techniques including: Laser Scanning Confocal Microscopy, Laser Scanning Two Photon Microscopy, Light Sheet Fluorescence Microscopy, Magnetic Resonance Imaging, Positron Emission To-mography, Optical Coherence Tomography, 3D Ultrasound Imaging, etc. A common denominator of all those applications being the constant need for further increasing resolution and quality of the acquired images. Interestingly, in some of the mentioned three-dimensional imaging techniques a remarkable situation arises: while a single volume does not contain enough information to represent the object being imaged within the quality parameters required by the final application, the acquisition scheme allows recording several volumes which represent different views of a given object, with each of the views providing complementary information. In this kind of situation one can get a better understanding of the object by combining several views instead of looking at each of them separately. Within such context, in this PhD Thesis we propose, develop and test new image processing methodologies based on the discrete wavelet transform for the combination, or fusion, of several views containing complementary information of a given object. The proposed fusion method exploits the scale and orientation decomposition capabil¬ities of the discrete wavelet transform to integrate in a single volume all the available information distributed among the set of acquired views. The work focuses in two different biomedical imaging modalities which provide such multi-view datasets. The first one is a particular fluorescence microscopy technique, Light-Sheet Fluorescence Microscopy, used for imaging and gaining understanding of the early development of live embryos from different animal models (like zebrafish or sea urchin). The second is Delayed Enhancement Magnetic Resonance Imaging, which is a valuable tool for assessing the viability of myocardial tissue on patients suffering from different cardiomyopathies. As part of this work, the proposed method was implemented and then validated on both imaging modalities. For the fluorescence microscopy application, the fusion results show improved contrast and detail discrimination when compared to any of the individual views and the method does not rely on prior knowledge of the system’s point spread function (PSF). Moreover, the results have shown improved performance with respect to previous PSF independent methods. With respect to its application to Delayed Enhancement Magnetic Resonance Imaging, the resulting fused volumes show a quantitative sharpness improvement and enable an easier and more complete interpretation of complex three-dimensional scar and heterogeneous tissue information in ischemic cardiomyopathy patients. In both applications, the results of this thesis are currently in use in the clinical and research centers with which the author collaborated during his work. An imple¬mentation of the fusion method has also been made freely available to the scientific community. Finally, an international patent application has been filed covering the visualization method developed for the Magnetic Resonance Imaging application.

Are Human Peripheral Nerves Sensitive To X-ray Imaging?

Diagnostic imaging techniques play an important role in assessing the exact location, cause, and extent of a nerve lesion, thus allowing clinicians to diagnose and manage more effectively a variety of pathological conditions, such as entrapment syndromes, traumatic injuries, and space-occupying lesions. Ultrasound and nuclear magnetic resonance imaging are becoming useful methods for this purpose, but they still lack spatial resolution. In this regard, recent phase contrast x-ray imaging experiments of peripheral nerve allowed the visualization of each nerve fiber surrounded by its myelin sheath as clearly as optical microscopy. In the present study, we attempted to produce high-resolution x-ray phase contrast images of a human sciatic nerve by using synchrotron radiation propagation-based imaging. The images showed high contrast and high spatial resolution, allowing clear identification of each fascicle structure and surrounding connective tissue. The outstanding result is the detection of such structures by phase contrast x-ray tomography of a thick human sciatic nerve section. This may further enable the identification of diverse pathological patterns, such as Wallerian degeneration, hypertrophic neuropathy, inflammatory infiltration, leprosy neuropathy and amyloid deposits. To the best of our knowledge, this is the first successful phase contrast x-ray imaging experiment of a human peripheral nerve sample. Our long-term goal is to develop peripheral nerve imaging methods that could supersede biopsy procedures.

Fast Transforms for Acoustic Imaging-Part I: Theory

The classical approach for acoustic imaging consists of beamforming, and produces the source distribution of interest convolved with the array point spread function. This convolution smears the image of interest, significantly reducing its effective resolution. Deconvolution methods have been proposed to enhance acoustic images and have produced significant improvements. Other proposals involve covariance fitting techniques, which avoid deconvolution altogether. However, in their traditional presentation, these enhanced reconstruction methods have very high computational costs, mostly because they have no means of efficiently transforming back and forth between a hypothetical image and the measured data. In this paper, we propose the Kronecker Array Transform ( KAT), a fast separable transform for array imaging applications. Under the assumption of a separable array, it enables the acceleration of imaging techniques by several orders of magnitude with respect to the fastest previously available methods, and enables the use of state-of-the-art regularized least-squares solvers. Using the KAT, one can reconstruct images with higher resolutions than was previously possible and use more accurate reconstruction techniques, opening new and exciting possibilities for acoustic imaging.

Simultaneous X-ray imaging of plant root growth and water uptake in thin-slab systems

Direct and simultaneous observation of root growth and plant water uptake is difficult because soils are opaque. X-ray imaging techniques such as projection radiography or Computer Tomography (CT) offer a partial alternative to such limitations. Nevertheless, there is a trade-off between resolution, large field-of-view and 3-dimensionality: With the current state of the technology, it is possible to have any two. In this study, we used X-ray transmission through thin-slab systems to monitor transient saturation fields that develop around roots as plants grow. Although restricted to 2-dimensions, this approach offers a large field-of-view together with high spatial and dynamic resolutions. To illustrate the potential of this technology, we grew peas in 1 cm thick containers filled with soil and imaged them at regular intervals. The dynamics of both the root growth and the water content field that developed around the roots could be conveniently monitored. Compared to other techniques such as X-ray CT, our system is relatively inexpensive and easy to implement. It can potentially be applied to study many agronomic problems, such as issues related to the impact of soil constraints (physical, chemical or biological) on root development.

MR Imaging Findings of Iron Overload

Hemochromatosis can be classified as (a) primary, when it originates from a genetic disturbance that promotes the increase of iron absorption, or (b) secondary, when it relates to chronic diseases or to multiple transfusions. The distribution of iron accumulation differs between these two forms; therefore, they can be distinguished by using imaging methods in the majority of cases. Magnetic resonance (MR) imaging is the most sensitive and specific imaging modality in the diagnosis of hemochromatosis. The susceptibility effect caused by the accumulation of iron leads to signal loss in the affected tissues, particularly with the T2*-weighted sequences, which makes the diagnosis of iron overload possible. By using MR imaging techniques, it is possible to estimate the hepatic iron concentration in a noninvasive way, thereby avoiding repeated biopsies. Hemochromatosis can lead to complications, such as a higher frequency of neoplasia, particularly the development of hepatocellular carcinoma. Other neoplasms, such as colorectal tumors, are also associated. Complications related to the treatment of chronic anemia include the appearance of peliosis hepatis and tumors, which can regress after the suspension of treatment with drugs. Knowledge of the disease and of the patterns of iron deposition in patients with iron overload enables not only diagnosis, but also treatment, follow-up, and the detection of possible complications by using imaging methods. (C) RSNA, 2009 . radiographics.rsna.org

Comparison of Quantitative Imaging Devices and Subjective Optic Nerve Head Assessment by General Ophthalmologists to Differentiate Normal From Glaucomatous Eyes

Purpose: To compare the ability of Subjective assessment of optic nerve head (ONH) and retinal nerve fiber layer (RNFL) by general ophthalmologists and by a glaucoma expert with objective measurements by optical coherence tomography (Stratus OCT, Carl Zeiss Meditec Inc), confocal scanning laser ophthalmoscope (HRT III; Heidelberg Engineering, Heidelberg. Germany), and scanning laser polarimetry (GDx enhanced corneal compensation; Carl Zeiss Meditec Inc, Dublin, CA) in discriminating glaucomatous and normal eyes. Methods: Sixty-one glaucomatous and 57 normal eyes or 118 subjects Were included in the study. Three independent general ophthalmologists and I glaucoma expert evaluated ONH stereo-photographs. Receiver operating characteristic curves were constructed for each imaging technique and sensitivity at fixed specificity was estimated. Comparisons or areas under these curves (aROCs) and agreement (k) were determined between stereophoto grading and best parameter from each technique. Results: Best parameter from each technique showed larger aROC (Stratus OCT RNFL 0.92; Stratus OCT ONH vertical integrated area = 0.86; Stratus OCT macular thickness = 0.82; GDx enhanced corneal compensation = 0.91, HRT3 global cup-to-disc ratio = 0.83; HRT3 glaucoma probability score numeric area score 0.83) compared with stereophotograph grading by general ophthalmologists (0.80) in separating glaucomatous and normal eyes. Glaucoma expert stereophoto grading provided equal or larger aROC (0.92) than best parameter of each computerized imaging device. Stereophoto evaluated by a glaucoma expert showed better agreement with best parameter of each quantitative imaging technique in classifying eyes either as glaucomatous or normal compared with stereophoto grading by general ophthalmologists, The combination Of Subjective assessment of the optic disc by general ophthalmologists with RNFL objective parameters improved identification of glaucoma patients in a larger proportion than the combination of these objective parameters with Subjective assessment of the optic disc by a glaucoma expert (29.5% vs. 19.7%, respectively). Conclusions: Diagnostic ability of all imaging techniques showed better performance than subjective assessment of the ONH by general ophthalmologists, but not by It glaucoma expert, Objective RNFL measurements may provide improvement in glaucoma detection when combined with subjective assessment of the optic disc by general ophthalmologists or by a glaucoma expert.

Findings of Pelvic Endometriosis at Transvaginal US, MR Imaging, and Laparoscopy

Endometriosis is a common multifocal gynecologic disease that manifests during the reproductive years, often causing chronic pelvic pain and infertility. It may occur as invasive peritoneal fibrotic nodules and adhesions or as ovarian cysts with hemorrhagic content. Although findings at physical examination may be suggestive, imaging is necessary for definitive diagnosis, patient counseling, and treatment planning. The imaging techniques that are most useful for preoperative disease mapping are transvaginal ultrasonography (US) after bowel preparation, and magnetic resonance (MR) imaging. Initial transvaginal US is a reliable technique for detecting rectosigmoid endometriotic lesions. MR imaging is indicated as a complementary examination in complex cases of endometriosis with extensive adhesions and ureteral involvement. Peritoneal endometriotic implants are typically hypoechoic on transvaginal US images and demonstrate low signal intensity on T2-weighted MR images. Endometriotic implants most commonly are found in retrocervical and rectosigmoid sites, followed by the vagina, bladder, and ureters. Cysts with low-level internal echoes and echogenic peripheral foci at transvaginal US are suggestive of endometriomas. MR imaging has high specificity for identifying endometriomas, which are characterized by high signal intensity on T1-weighted images and low signal intensity on T2-weighted images. Correlation of the radiologic imaging features of endometriotic lesions with their laparoscopic appearances may help improve individual proficiency in the radiologic diagnosis of endometriosis. Supplemental material available at http:// radiographics.rsna.org/lookup/suppl/doi: 10.1148/rg.314105193/-/DC1. (C) RSNA, 2011 radiographics.rsna.org

Advances in Imaging of Osteoarthritis and Cartilage

Osteoarthritis (OA) is the most frequent form of arthritis, with major implications for individual and public health care without effective treatment available. The field of joint imaging, and particularly magnetic resonance (MR) imaging, has evolved rapidly owing to technical advances and the application of these to the field of clinical research. Cartilage imaging certainly is at the forefront of these developments. In this review, the different aspects of OA imaging and cartilage assessment, with an emphasis on recent advances, will be presented. The current role of radiography, including advances in the technology for joint space width assessment, will be discussed. The development of various MR imaging techniques capable of facilitating assessment of cartilage morphology and the methods for evaluating the biochemical composition of cartilage will be presented. Advances in quantitative morphologic cartilage assessment and semiquantitative whole-organ assessment will be reviewed. Although MR imaging is the most important modality in imaging of OA and cartilage, others such as ultrasonography play a complementary role that will be discussed briefly.

Machine Learning Algorithms at Myocardial Perfusion Imaging - a Preliminary Study

Introduction: A major focus of data mining process - especially machine learning researches - is to automatically learn to recognize complex patterns and help to take the adequate decisions strictly based on the acquired data. Since imaging techniques like MPI – Myocardial Perfusion Imaging on Nuclear Cardiology, can implicate a huge part of the daily workflow and generate gigabytes of data, there could be advantages on Computerized Analysis of data over Human Analysis: shorter time, homogeneity and consistency, automatic recording of analysis results, relatively inexpensive, etc.Objectives: The aim of this study relates with the evaluation of the efficacy of this methodology on the evaluation of MPI Stress studies and the process of decision taking concerning the continuation – or not – of the evaluation of each patient. It has been pursued has an objective to automatically classify a patient test in one of three groups: “Positive”, “Negative” and “Indeterminate”. “Positive” would directly follow to the Rest test part of the exam, the “Negative” would be directly exempted from continuation and only the “Indeterminate” group would deserve the clinician analysis, so allowing economy of clinician’s effort, increasing workflow fluidity at the technologist’s level and probably sparing time to patients. Methods: WEKA v3.6.2 open source software was used to make a comparative analysis of three WEKA algorithms (“OneR”, “J48” and “Naïve Bayes”) - on a retrospective study using the comparison with correspondent clinical results as reference, signed by nuclear cardiologist experts - on “SPECT Heart Dataset”, available on University of California – Irvine, at the Machine Learning Repository. For evaluation purposes, criteria as “Precision”, “Incorrectly Classified Instances” and “Receiver Operating Characteristics (ROC) Areas” were considered. Results: The interpretation of the data suggests that the Naïve Bayes algorithm has the best performance among the three previously selected algorithms. Conclusions: It is believed - and apparently supported by the findings - that machine learning algorithms could significantly assist, at an intermediary level, on the analysis of scintigraphic data obtained on MPI, namely after Stress acquisition, so eventually increasing efficiency of the entire system and potentially easing both roles of Technologists and Nuclear Cardiologists. In the actual continuation of this study, it is planned to use more patient information and significantly increase the population under study, in order to allow improving system accuracy.

OPTIMAX 2014 - Radiation dose and image quality optimisation in medical imaging

Medical imaging is a powerful diagnostic tool. Consequently, the number of medical images taken has increased vastly over the past few decades. The most common medical imaging techniques use X-radiation as the primary investigative tool. The main limitation of using X-radiation is associated with the risk of developing cancers. Alongside this, technology has advanced and more centres now use CT scanners; these can incur significant radiation burdens compared with traditional X-ray imaging systems. The net effect is that the population radiation burden is rising steadily. Risk arising from X-radiation for diagnostic medical purposes needs minimising and one way to achieve this is through reducing radiation dose whilst optimising image quality. All ages are affected by risk from X-radiation however the increasing population age highlights the elderly as a new group that may require consideration. Of greatest concern are paediatric patients: firstly they are more sensitive to radiation; secondly their younger age means that the potential detriment to this group is greater. Containment of radiation exposure falls to a number of professionals within medical fields, from those who request imaging to those who produce the image. These staff are supported in their radiation protection role by engineers, physicists and technicians. It is important to realise that radiation protection is currently a major European focus of interest and minimum competence levels in radiation protection for radiographers have been defined through the integrated activities of the EU consortium called MEDRAPET. The outcomes of this project have been used by the European Federation of Radiographer Societies to describe the European Qualifications Framework levels for radiographers in radiation protection. Though variations exist between European countries radiographers and nuclear medicine technologists are normally the professional groups who are responsible for exposing screening populations and patients to X-radiation. As part of their training they learn fundamental principles of radiation protection and theoretical and practical approaches to dose minimisation. However dose minimisation is complex – it is not simply about reducing X-radiation without taking into account major contextual factors. These factors relate to the real world of clinical imaging and include the need to measure clinical image quality and lesion visibility when applying X-radiation dose reduction strategies. This requires the use of validated psychological and physics techniques to measure clinical image quality and lesion perceptibility.

Development of procedures for the design, optimization and manufacturing of customized orthopaedic and trauma implants: Geometrical/anatomical modelling from 3D medical imaging

Tese de Doutoramento (Programa Doutoral em Engenharia Biomédica)

Multimodality Imaging of Heart Valve Disease

Unidentified heart valve disease is associated with a significant morbidity and mortality. It has therefore become important to accurately identify, assess and monitor patients with this condition in order that appropriate and timely intervention can occur. Although echocardiography has emerged as the predominant imaging modality for this purpose, recent advances in cardiac magnetic resonance and cardiac computed tomography indicate that they may have an important contribution to make. The current review describes the assessment of regurgitant and stenotic heart valves by multimodality imaging (echocardiography, cardiac computed tomography and cardiac magnetic resonance) and discusses their relative strengths and weaknesses.

Stem cell transplantation in brain tumors: a new field for molecular imaging?

Neural stem cells have been proposed as a new and promising treatment modality in various pathologies of the central nervous system, including malignant brain tumors. However, the underlying mechanism by which neural stem cells target tumor areas remains elusive. Monitoring of these cells is currently done by use of various modes of molecular imaging, such as optical imaging, magnetic resonance imaging and positron emission tomography, which is a novel technology for visualizing metabolism and signal transduction to gene expression. In this new context, the microenvironment of (malignant) brain tumors and the blood-brain barrier gains increased interest. The authors of this review give a unique overview of the current molecular-imaging techniques used in different therapeutic experimental brain tumor models in relation to neural stem cells. Such methods for molecular imaging of gene-engineered neural stem/progenitor cells are currently used to trace the location and temporal level of expression of therapeutic and endogenous genes in malignant brain tumors, closing the gap between in vitro and in vivo integrative biology of disease in neural stem cell transplantation.

Electroencephalographic source imaging: a prospective study of 152 operated epileptic patients.

Electroencephalography is mandatory to determine the epilepsy syndrome. However, for the precise localization of the irritative zone in patients with focal epilepsy, costly and sometimes cumbersome imaging techniques are used. Recent small studies using electric source imaging suggest that electroencephalography itself could be used to localize the focus. However, a large prospective validation study is missing. This study presents a cohort of 152 operated patients where electric source imaging was applied as part of the pre-surgical work-up allowing a comparison with the results from other methods. Patients (n = 152) with >1 year postoperative follow-up were studied prospectively. The sensitivity and specificity of each imaging method was defined by comparing the localization of the source maximum with the resected zone and surgical outcome. Electric source imaging had a sensitivity of 84% and a specificity of 88% if the electroencephalogram was recorded with a large number of electrodes (128-256 channels) and the individual magnetic resonance image was used as head model. These values compared favourably with those of structural magnetic resonance imaging (76% sensitivity, 53% specificity), positron emission tomography (69% sensitivity, 44% specificity) and ictal/interictal single-photon emission-computed tomography (58% sensitivity, 47% specificity). The sensitivity and specificity of electric source imaging decreased to 57% and 59%, respectively, with low number of electrodes (<32 channels) and a template head model. This study demonstrated the validity and clinical utility of electric source imaging in a large prospective study. Given the low cost and high flexibility of electroencephalographic systems even with high channel counts, we conclude that electric source imaging is a highly valuable tool in pre-surgical epilepsy evaluation.

Motion compensation strategies in magnetic resonance imaging.

Image quality in magnetic resonance imaging (MRI) is considerably affected by motion. Therefore, motion is one of the most common sources of artifacts in contemporary cardiovascular MRI. Such artifacts in turn may easily lead to misinterpretations in the images and a subsequent loss in diagnostic quality. Hence, there is considerable research interest in strategies that help to overcome these limitations at minimal cost in time, spatial resolution, temporal resolution, and signal-to-noise ratio. This review summarizes and discusses the three principal sources of motion: the beating heart, the breathing lungs, and bulk patient movement. This is followed by a comprehensive overview of commonly used compensation strategies for these different types of motion. Finally, a summary and an outlook are provided.