Fast Microwave Tomography Algorithm for Breast Cancer Imaging

Microwave tomography has shown promise for breast cancer imaging. The microwaves are harmless to body tissues, which makes microwave tomography a safe adjuvant screening to mammography. Although many clinical studies have shown the effectiveness of regular screening for the detection of breast cancer, the anatomy of the breast and its critical tissues challenge the identification and diagnosis of tumors in this region. Detection of tumors in the breast is more challenging in heterogeneously dense and extremely dense breasts, and microwave tomography has the potential to be effective in such cases. The sensitivity of microwaves to various breast tissues and the comfort and safety of the screening method have made microwave tomography an attractive imaging technique. Despite the need for an alternative screening technique, microwave tomography has not yet been introduced as a screening modality in regular health care, and is still subject to research. The main obstacles are imperfect hardware systems and inefficient imaging algorithms. The immense computational costs for the image reconstruction algorithm present a crucial challenge. 2D imaging algorithms are proposed to reduce the amount of hardware resources required and the imaging time. Although 2D microwave tomography algorithms are computationally less expensive, few imaging groups have been successful in integrating the acquired 3D data into the 2D tomography algorithms for clinical applications. The microwave tomography algorithms include two main computation problems: the forward problem and the inverse problem. The first part of this thesis focuses on a new fast forward solver, the 2D discrete dipole approximation (DDA), which is formulated and modeled. The effect of frequency, sampling number, target size, and contrast on the accuracy of the solver are studied. Additionally, the 2D DDA time efficiency and computation time as a single forward solver are investigated.\ua0 The second part of this thesis focuses on the inverse problem. This portion of the algorithm is based on a log-magnitude and phase transformation optimization problem and is formulated as the Gauss-Newton iterative algorithm. The synthetic data from a finite-element-based solver (COMSOL Multiphysics) and the experimental data acquired from the breast imaging system at Chalmers University of Technology are used to evaluate the DDA-based image reconstruction algorithm. The investigations of modeling and computational complexity show that the 2D DDA is a fast and accurate forward solver that can be embedded in tomography algorithms to produce images in seconds. The successful development and implementation in this thesis of 2D tomographic breast imaging with acceptable accuracy and high computational cost efficiency has provided significant savings in time and in-use memory and is a dramatic improvement over previous implementations

Sparse-Representation-Based Direct Minimum L

A sparse-representation-based direct minimum Lp-norm algorithm is proposed for a two-dimensional MRI phase unwrapping. First, the algorithm converts the weighted-Lp-norm-minimization-based phase unwrapping problem into a linear system problem whose system (coefficient) matrix is a large, symmetric one. Then, the coefficient-matrix is represented in the sparse structure. Finally, standard direct solvers are employed to solve this linear system. Several wrapped phase datasets, including simulated and MR data, were used to evaluate this algorithm’s performance. The results demonstrated that the proposed algorithm for unwrapping MRI phase data is reliable and robust

Microwave Tomography Using Stochastic Optimization And High Performance Computing

This thesis discusses the application of parallel computing in microwave tomography for detection and imaging of dielectric objects. The main focus is on microwave tomography with the use of a parallelized Finite Difference Time Domain (FDTD) forward solver in conjunction with non-linear stochastic optimization based inverse solvers. Because such solvers require very heavy computation, their investigation has been limited in favour of deterministic inverse solvers that make use of assumptions and approximations of the imaging target. Without the use of linearization assumptions, a non-linear stochastic microwave tomography system is able to resolve targets of arbitrary permittivity contrast profiles while avoiding convergence to local minima of the microwave tomography optimization space. This work is focused on ameliorating this computational load with the use of heavy parallelization. The presented microwave tomography system is capable of modelling complex, heterogeneous, and dispersive media using the Debye model. A detailed explanation of the dispersive FDTD is presented herein. The system uses scattered field data due to multiple excitation angles, frequencies, and observation angles in order to improve target resolution, reduce the ill-posedness of the microwave tomography inverse problem, and improve the accuracy of the complex permittivity profile of the imaging target. The FDTD forward solver is parallelized with the use of the Common Unified Device Architecture (CUDA) programming model developed by NVIDIA corporation. In the forward solver, the time stepping of the fields are computed on a Graphics Processing Unit (GPU). In addition the inverse solver makes use of the Message Passing Interface (MPI) system to distribute computation across multiple work stations. The FDTD method was chosen due to its ease of parallelization using GPU computing, in addition to its ability to simulate wideband excitation signals during a single forward simulation. We investigated the use of distributed Particle Swarm Optimization (PSO) and Differential Evolution (DE) methods in the inverse solver for this microwave tomography system. In these optimization algorithms, candidate solutions are farmed out to separate workstations to be evaluated. As fitness evaluations are returned asynchronously, the optimization algorithm updates the population of candidate solutions and gives new candidate solutions to be evaluated to open workstations. In this manner, we used a total of eight graphics processing units during optimization with minimal downtime. Presented in this thesis is a microwave tomography algorithm that does not rely on linearization assumptions, capable of imaging a target in a reasonable amount of time for clinical applications. The proposed algorithm was tested using numerical phantoms that with material parameters similar to what one would find in normal or malignant human tissue

Engineering Parallel Transmit/Receive Radio-Frequency Coil Arrays for Human Brain MRI at 7 Tesla

Magnetic resonance imaging is widely used in medical diagnosis to obtain anatomical details of the human body in a non-invasive way. Clinical MR scanners typically operate at a static magnetic field strength (B0) of 1.5T or 3T. However, going to higher field is of great interest since the signal-to-noise ratio is proportional to B0. Therefore, higher image resolution and better contrast between the human tissues could be achieved. Nevertheless, new challenges arise when increasing B0. The wavelength associated with the radio-frequency field B1+ has smaller dimensions - approx. 12 cm for human brain tissues - than the human brain itself (20 cm in length), the organ of interest in this thesis. The main consequence is that the transmit field distribution pattern (B1+) is altered and the final MR images present bright and dark signal spots. These effects prevent the ultra-high field MR scanners (>= 7T) to be used for routine clinical diagnosis. Parallel-transmit is one approach to address these new challenges. Instead of using an RF coil connected to a single power input as it is commonly done at lower magnetic fields, multiple RF coils are used with independent power inputs. The subsequent distinct RF signals can be manipulated separately, which provides an additional degree of freedom to generate homogeneous B1+-field distributions over large or specific regions in the human body. A transmit/receive RF coil array optimized for whole-brain MR imaging was developed and is described in this thesis. Dipoles antennas were used since they could provide a large longitudinal (vertical axis-head to neck) coverage and high transmit field efficiency. Results demonstrated a complete coverage of the human brain, and particularly high homogeneity over the cerebellum. However, since the receive sensitivity over large field-of-views is related to the number of channels available to detect the NMR signal, the next work was to add a 32-channel receive loop coil array to the transmit coil array. The complete coverage of the human brain was assessed with a substantial increase in signal-to-noise compared to the transmit/receive dipole coil array alone. Moreover, acquisition time was shortened since higher acceleration factors could be used. To optimize the individual RF fields and generate an homogeneous B1+-field, a method was developed making use of the particle-swarm algorithm. A user-friendly graphical interface was implemented. Good homogeneity could be achieved over the whole-brain after optimization with the coil array built in this study. Moreover, the optimization was shown to be robust across multiple subjects. The last project was focused on the single transmit system. Local volume coils (single transmit) present pronounced transmit field inhomogeneities in specific regions of the human brain such as the temporal lobes. A widely used approach to address locally these challenges is to add dielectric pads inside the volume coils to enhance the local transmit field efficiency. It was shown in this thesis that constructing dedicated surface coils is a valuable alternative to the dielectric pads in terms of transmit field efficiency and MR spectroscopy results. Two RF coil setups were developed for the temporal and frontal lobes of the human brain, respectively. This thesis provides extensive insight on MR engineering of RF coils at ultra-high field and the potential of parallel-transmit to address the future needs in clinical applications

DEVELOPMENT OF PHASE RETRIEVAL TECHNIQUES IN DISPLACEMENT AND DISPLACEMENT DERIVATIVES MEASUREMENT

Ph.DDOCTOR OF PHILOSOPH

MATLAB

A well-known statement says that the PID controller is the "bread and butter" of the control engineer. This is indeed true, from a scientific standpoint. However, nowadays, in the era of computer science, when the paper and pencil have been replaced by the keyboard and the display of computers, one may equally say that MATLAB is the "bread" in the above statement. MATLAB has became a de facto tool for the modern system engineer. This book is written for both engineering students, as well as for practicing engineers. The wide range of applications in which MATLAB is the working framework, shows that it is a powerful, comprehensive and easy-to-use environment for performing technical computations. The book includes various excellent applications in which MATLAB is employed: from pure algebraic computations to data acquisition in real-life experiments, from control strategies to image processing algorithms, from graphical user interface design for educational purposes to Simulink embedded systems

Steady-state anatomical and quantitative magnetic resonance imaging of the heart using RF-frequencymodulated techniques

Cardiovascular disease (CVD) is the leading cause of death in the United States and Europe and generates healthcare costs of hundreds of billions of dollars annually. Conventional methods of diagnosing CVD are often invasive and carry risks for the patient. For example, the gold standard for diagnosing coronary artery disease, a major class of CVD, is x-ray coronary angiography, which has the disadvantages of being invasive, being expensive, using ionizing radiation, and having a ris k of complications. Conversely, coronary MR angiography (MRA) does not use ionizing radiation, can effectively visualize tissues without the need for exogenous contrast agents, and benefits from an adaptable temporal resolution. However, the acquisition time of cardiac MRI is far longer than the temporal scales of cardiac and respiratory motion, necessitating some method of compensating for this motion. The free-running framework is a novel development in our lab, benefitting from advances over the past three decades, that attempts to address disadvantages of previous cardiac MRI approaches: it provides fully self-gated 5D cardiac MRI with a simplified workflow, improved ease-of-use, reduced operator dependence, and automatic patient-specific motion detection. Free-running imaging increases the amount of information available to the clinician and is flexible enough to be translated to different app lications within cardiac MRI. Moreover, the self-gating of the free-running framework decoupled the acquisition from the motion compensation and thereby opened up cardiac MRI to the wider class of steady-state-based techniques utilizing balanced steady-state free precession (bSSFP) sequences, which have the benefits of practical simplicity and high signal-to-noise ratio. The focus of this thesis was therefore on the application of steady- state techniques to cardiac MRI. The first part addressed the long acquisition time of the current free-running framework and focused on anatomical coronary imaging. The published protocol of the free- running framework used an interrupted bSSFP acquisition where CHESS fat saturation modules were inserted to provide blood-fat contrast, as they suppress the signal of fat tissue surrounding the coronary arteries, and were followed by ramp-up pulses to reduce artefacts arising from the return to steady-state. This interrupted acquisition, however, suffered from an interrupted steady-state, reduced time efficiency, and higher specific absorption rate (SAR). Using novel lipid-insensitive binomial off-resonant RF excitation (LIBRE) pulses developed in our lab, the first project showed that LIBRE pulses incorporated into an uninterrupted free-running bSSFP sequence could be successfully used for 5D cardiac MRI at 1.5T. The free-running LIBRE approach reduced the acquisition time and SAR relative to the previous interrupted approach while maintaining image quality and vessel conspicuity. Furthermore, this had been the first successful use of a fat-suppressing RF excitation pulse in an uninterrupted bSSFP sequence for cardiac imaging, demonstrating that uninterrupted bSSFP can be used for cardiac MRI and addressing the problem of clinical sequence availability. Inspired by the feasibility of uninterrupted bSSFP for cardiac MRI, the second part investigated the potential of PLANET, a novel 3D multiparametric mapping technique, for free-running 5D myocardial mapping. PLANET utilizes a phase-cycled bSSFP acquisition and a direct ellipse-fitting algorithm to calculate T1 and T2 relaxation times, which suggested that it could be readily integrated into the free-running framework without interrupting the steady-state. After initially calibrating the acquisition, the possibility of accelerating the static PLANET acquisition was explored prior to applying it to the moving heart. It was shown that PLANET accuracy and precision could be maintained with two-fold acceleration with a 3D Cartesian spiral trajectory, suggesting that PLANET for myocardial mapping with the free-running 5D radial acquisition is feasible. Further work should investigate optimizing the reconstruction scheme, improving the coil sensitivity estimate, and examining the use of the radial trajectory with a view to implementing free-running 5D myocardial T1 and T2 mapping. This thesis presents two approaches utilizing RF-frequency-modulated steady-state techniques for cardiac MRI. The first approach involved the novel application of an uninterrupted bSSFP acquisition with off-resonant RF excitation for anatomical coronary imaging. The second approach investigated the use of phase-cycled bSSFP for free-running 5D myocardial T1 and T2 mapping. Both methods addressed the challenge of clinical availability of sequences in cardiac MRI, by showing that a common and simple sequence like bSSFP can be used for acquisition while the steps of motion compensation and reconstruction can be handled offline, and thus have the potential to improve adoption of cardiac MRI. -- Les maladies cardiovasculaires (MCV) représentent la principale cause de décès aux États-Unis et en Europe et génèrent des coûts de santé de plusieurs centaines de milliards de dollars par an. Les méthodes conventionnelles de diagnostic des MCV sont souvent invasives et comportent des risques pour le patient. Par exemple, la méthode de référence pour le diagnostic de la maladie coronarienne, une catégorie majeure de MCV, est la coronarographie par rayons X qui a comme inconvénients son caractère invasif, son coût, l’utilisation de rayonnements ionisants et le risque de complications. A l’inverse, l'angiographie coronarienne par résonance magnétique (ARM) n'utilise pas de rayonnements ionisants, permet de visualiser efficacement les tissus sans avoir recours à des agents de contraste exogènes et bénéficie d'une résolution temporelle ajustable. Cependant, le temps d'acquisition en IRM cardiaque est bien plus long que les échelles temporelles des mouvements cardiaques et respiratoires en jeu, ce qui rend la compensation de ces mouvements indispensable. Le cadre dit de « free -running » est un nouveau développement de notre laboratoire qui bénéficie des progrès réalisés au cours des trois dernières décennies et tente de remédier aux inconvénients des approches précédentes pour l'IRM cardiaque : il fournit une IRM cardiaque en cinq dimensions (5D) complètement « self-gated » , c’est-à-dire capable de détecter les mouvements cardiaques et respiratoires, forte d’une implémentation simplifiée, d’une plus grande facilité d'utilisation, d’une dépendance réduite vis-à-vis de l'opérateur et d’une détection automatique des mouvements spécifiques du patient. L'imagerie « free- running » augmente la quantité d'informations à disposition du clinicien et est suffisamment flexible pour être appliquée à différents domaines de l'IRM cardiaque. De plus, le « self-gating » du cadre « free-running » a découplé l'acquisition de la compensation de mouvement et a ainsi ouvert l'IRM cardiaque à la classe plus large des techniques basées sur l'état stationnaire utilisant des séquences de précession libre équilibrée en état stationnaire (bSSFP), qui se distinguent par leur simplicité d’utilisation et leur rapport signal sur bruit élevé. Le thème de cette thèse est donc l'application des techniques basées sur l'état stationnaire à l'IRM cardiaque. La première partie porte sur le long temps d'acquisition de l'actuel cadre « free-running» et se concentre sur l'imagerie anatomique coronaire. Le protocole publié utilise une acquisition bSSFP interrompue où des modules de saturation de graisse (CHESS) sont insérés de façon à fournir un contraste sang-graisse puisqu’ils suppriment le signal du tissu graisseux entourant les artères coronaires, et sont suivis par des impulsions en rampe pour réduire les artefacts résultant du retour à l'état stable. Cette acquisition interrompue souffre cependant d'un état d'équilibre interrompu, d'une efficacité temporelle réduite et d'un débit d'absorption spécifique (DAS) plus élevé. En utilisant les nouvelles impulsions d'excitation radiofréquence (RF) binomiales hors -résonance insensibles aux lipides (LIBRE) développées dans notre laboratoi re, ce premier projet montre que les impulsions LIBRE incorporées dans une séquence bSSFP ininterrompue et « free-running » peuvent être utilisées avec succès pour l'IRM cardiaque 5D à 1,5 T. L'approche « free-running LIBRE » permet de réduire le temps d'acquisition et le DAS par rapport à l'approche interrompue précédente, tout en maintenant la perceptibilité des artères coronariennes. En outre, il s'agit de la première utilisation réussie d'une impulsion d'excitation RF supprimant la graisse dans une séquence bSSFP ininterrompue pour l'imagerie cardiaque, ce qui démontre le potentiel d’utilisation de la séquence bSSFP ininterrompue pour l'IRM cardiaque et résout le problème de la disponibilité de la séquence en clinique. Inspirée par la faisabilité d’utilisation de la séquence bSSFP ininterrompue pour l'IRM cardiaque, la deuxième partie étudie le potentiel de PLANET, une nouvelle technique de cartographie 3D multiparamétrique, pour la cartographie 5D du myocarde via l’imagerie « free-running ». PLANET utilise une acquisition bSSFP à cycle de phase et un algorithme d'ajustement d'ellipse direct pour calculer les temps de relaxation T1 et T2, ce qui suggère que cette méthode pourrait être facilement intégrée au cadre « free - running » sans interruption de l’état d'équilibre. Après calibration de l'acquisition, nous explorons la possibilité d'accélérer l'acquisition statique de PLANET pour l'appliquer au cœur. Nous démontrons que l'exactitude et la précision de PLANET peuvent être maintenues pour une accélération double avec une trajectoire 3D cartésienne en spirale, ce qui suggère que PLANET est réalisable pour la cartographie du myocarde avec une acquisition radiale 5D « free-running ». D'autres travaux devraient porter sur l'optimisation du schéma de reconstruction, l'amélioration de l'estimation de la sensibilité de l’antenne et l'examen de l'utilisation de la trajectoire radiale en vue de la mise en œuvre de la cartographie 5D « free-running » T1 et T2 du myocarde. Cette thèse présente deux approches utilisant des techniques de modulation de fréquence radio en état stationnaire pour l'IRM cardiaque. La première approche implique l'application nouvelle d'une acquisition bSSFP ininterrompue avec une excitation RF hors résonance pour l'imagerie anatomique coronaire. La seconde approche porte sur l'utilisation d’une séquence bSSFP à cycle de phase pour la cartographie 5D T1 et T2 du myocarde. Ces deux méthodes permettent de répondre au défi posé par la disponibilité des séquences en IRM cardiaque en montrant qu'une séquence commune et simple comme la bSSFP peut être utilisée pour l'acquisition, tandis que les étapes de compensation du mouvement et de reconstruction peuvent être traitées hors ligne. Ainsi, ces méthodes ont le potentiel de favoriser l'adoption de l'IRM cardiaque

Quantification and segmentation of breast cancer diagnosis: efficient hardware accelerator approach

The mammography image eccentric area is the breast density percentage measurement. The technical challenge of quantification in radiology leads to misinterpretation in screening. Data feedback from society, institutional, and industry shows that quantification and segmentation frameworks have rapidly become the primary methodologies for structuring and interpreting mammogram digital images. Segmentation clustering algorithms have setbacks on overlapping clusters, proportion, and multidimensional scaling to map and leverage the data. In combination, mammogram quantification creates a long-standing focus area. The algorithm proposed must reduce complexity and target data points distributed in iterative, and boost cluster centroid merged into a single updating process to evade the large storage requirement. The mammogram database's initial test segment is critical for evaluating performance and determining the Area Under the Curve (AUC) to alias with medical policy. In addition, a new image clustering algorithm anticipates the need for largescale serial and parallel processing. There is no solution on the market, and it is necessary to implement communication protocols between devices. Exploiting and targeting utilization hardware tasks will further extend the prospect of improvement in the cluster. Benchmarking their resources and performance is required. Finally, the medical imperatives cluster was objectively validated using qualitative and quantitative inspection. The proposed method should overcome the technical challenges that radiologists face