해부학적 유도 PET 재구성: 매끄럽지 않은 사전 함수부터 딥러닝 접근까지

학위논문 (박사) -- 서울대학교 대학원 : 의과대학 의과학과, 2021. 2. 이재성.Advances in simultaneous positron emission tomography/magnetic resonance imaging (PET/MRI) technology have led to an active investigation of the anatomy-guided regularized PET image reconstruction algorithm based on MR images. Among the various priors proposed for anatomy-guided regularized PET image reconstruction, Bowsher’s method based on second-order smoothing priors sometimes suffers from over-smoothing of detailed structures. Therefore, in this study, we propose a Bowsher prior based on the l1 norm and an iteratively reweighting scheme to overcome the limitation of the original Bowsher method. In addition, we have derived a closed solution for iterative image reconstruction based on this non-smooth prior. A comparison study between the original l2 and proposed l1 Bowsher priors were conducted using computer simulation and real human data. In the simulation and real data application, small lesions with abnormal PET uptake were better detected by the proposed l1 Bowsher prior methods than the original Bowsher prior. The original l2 Bowsher leads to a decreased PET intensity in small lesions when there is no clear separation between the lesions and surrounding tissue in the anatomical prior. However, the proposed l1 Bowsher prior methods showed better contrast between the tumors and surrounding tissues owing to the intrinsic edge-preserving property of the prior which is attributed to the sparseness induced by l1 norm, especially in the iterative reweighting scheme. Besides, the proposed methods demonstrated lower bias and less hyper-parameter dependency on PET intensity estimation in the regions with matched anatomical boundaries in PET and MRI. Moreover, based on the formulation of l1 Bowsher prior, the unrolled network containing the conventional maximum-likelihood expectation-maximization (ML-EM) module was also proposed. The convolutional layers successfully learned the distribution of anatomically-guided PET images and the EM module corrected the intermediate outputs by comparing them with sinograms. The proposed unrolled network showed better performance than ordinary U-Net, where the regional uptake is less biased and deviated. Therefore, these methods will help improve the PET image quality based on the anatomical side information.양전자방출단층촬영 / 자기공명영상 (PET/MRI) 동시 획득 기술의 발전으로 MR 영상을 기반으로 한 해부학적 사전 함수로 정규화 된 PET 영상 재구성 알고리즘에 대한 심도있는 평가가 이루어졌다. 해부학 기반으로 정규화 된 PET 이미지 재구성을 위해 제안 된 다양한 사전 중 2차 평활화 사전함수에 기반한 Bowsher의 방법은 때때로 세부 구조의 과도한 평활화로 어려움을 겪는다. 따라서 본 연구에서는 원래 Bowsher 방법의 한계를 극복하기 위해 l1 norm에 기반한 Bowsher 사전 함수와 반복적인 재가중치 기법을 제안한다. 또한, 우리는 이 매끄럽지 않은 사전 함수를 이용한 반복적 이미지 재구성에 대해 닫힌 해를 도출했다. 원래 l2와 제안 된 l1 Bowsher 사전 함수 간의 비교 연구는 컴퓨터 시뮬레이션과 실제 데이터를 사용하여 수행되었다. 시뮬레이션 및 실제 데이터에서 비정상적인 PET 흡수를 가진 작은 병변은 원래 Bowsher 이전보다 제안 된 l1 Bowsher 사전 방법으로 더 잘 감지되었다. 원래의 l2 Bowsher는 해부학적 영상에서 병변과 주변 조직 사이에 명확한 분리가 없을 때 작은 병변에서의 PET 강도를 감소시킨다. 그러나 제안 된 l1 Bowsher 사전 방법은 특히 반복적 재가중치 기법에서 l1 노름에 의해 유도된 희소성에 기인한 특성으로 인해 종양과 주변 조직 사이에 더 나은 대비를 보여주었다. 또한 제안된 방법은 PET과 MRI의 해부학적 경계가 일치하는 영역에서 PET 강도 추정에 대한 편향이 더 낮고 하이퍼 파라미터 종속성이 적음을 보여주었다. 또한, l1Bowsher 사전 함수의 닫힌 해를 기반으로 기존의 ML-EM (maximum-likelihood expectation-maximization) 모듈을 포함하는 펼쳐진 네트워크도 제안되었다. 컨볼루션 레이어는 해부학적으로 유도 재구성된 PET 이미지의 분포를 성공적으로 학습했으며, EM 모듈은 중간 출력들을 사이노그램과 비교하여 결과 이미지가 잘 들어맞게 수정했다. 제안된 펼쳐진 네트워크는 지역의 흡수선량이 덜 편향되고 편차가 적어, 일반 U-Net보다 더 나은 성능을 보여주었다. 따라서 이러한 방법들은 해부학적 정보를 기반으로 PET 이미지 품질을 향상시키는 데 유용할 것이다.Chapter 1. Introduction 1 1.1. Backgrounds 1 1.1.1. Positron Emission Tomography 1 1.1.2. Maximum a Posterior Reconstruction 1 1.1.3. Anatomical Prior 2 1.1.4. Proposed l_1 Bowsher Prior 3 1.1.5. Deep Learning for MR-less Application 4 1.2. Purpose of the Research 4 Chapter 2. Anatomically-guided PET Reconstruction Using Bowsher Prior 6 2.1. Backgrounds 6 2.1.1. PET Data Model 6 2.1.2. Original Bowsher Prior 7 2.2. Methods and Materials 8 2.2.1. Proposed l_1 Bowsher Prior 8 2.2.2. Iterative Reweighting 13 2.2.3. Computer Simulations 15 2.2.4. Human Data 16 2.2.5. Image Analysis 17 2.3. Results 19 2.3.1. Simulation with Brain Phantom 19 2.3.2.Human Data 20 2.4. Discussions 25 Chapter 3. Deep Learning Approach for Anatomically-guided PET Reconstruction 31 3.1. Backgrounds 31 3.2. Methods and Materials 33 3.2.1. Douglas-Rachford Splitting 33 3.2.2. Network Architecture 34 3.2.3. Dataset and Training Details 35 3.2.4. Image Analysis 36 3.3. Results 37 3.4. Discussions 38 Chapter 4. Conclusions 40 Bibliography 41 Abstract in Korean (국문 초록) 52Docto

Advanced industrial X-ray computed tomography for defect detection and characterisation of composite structures

X-ray Computer Tomography (CT) is well suited to the inspection of Fibre-Reinforced-Plastic (FRP) composite materials. However, a range of limitations currently restrict its uptake. The aim of the present research was to develop advanced inspection procedures that overcome these limitations and increase the scope of composite structures that can be inspected by industrial cone beam CT. Region of Interest (ROI) CT inspection of FRP laminated panels was investigated and two data completion methods developed to overcome reconstruction errors caused by truncated projection data. These allow accurate, highly magnified regions to be reconstructed on objects that extend beyond the Field-of-View (FOV) of the detector. The first method extended the truncated projection data using a cosine signal tailing off to zero attenuation. This method removed the strong 'glowing' artefacts but an inherent error existed across the reconstructed ROI. This did not affect the defect detectability of the inspection but was viewed as problematic for applications requiring accurate density measurements. The second method used prior knowledge of the test object so that a model could be created to estimate the missing data. This technique removed errors associated with ROI reconstruction thus significantly improving the accuracy. Techniques for extending the FOV were developed and applied to the inspection of FRP wind turbine blades; over 1.5X larger than the conventional scanning FOV. Two data completion methods were developed requiring an asymmetrically positioned detector. The first was based on the cosine tailing technique and the second used fan beam ray redundancy properties to estimate the missing data. Both produced accurate reconstructions for the 'offset' projection data, demonstrating that it was possible to approximately double the FOV. The cosine tailing method was found to be the more reliable. A dual energy image CT technique was developed to extend the optimum dynamic range and improve defect detectability for multi-density objects. This was applied to FRP composite/Titanium lap joints showing improved detectability of both volumetric and planar defects within the low density FRP. The dual energy procedure was validated using statistical performance measures on a specially fabricated multi-density phantom. The results showed a significant improvement in the detail SNR when compared to conventional CT scans.EThOS - Electronic Theses Online ServiceTWI LtdThe Engineering and Physical Sciences Research Board (EPSRC)GBUnited Kingdo

Reduction of Limited Angle Artifacts in Medical Tomography via Image Reconstruction

Artifacts are unwanted effects in tomographic images that do not reflect the nature of the object. Their widespread occurrence makes their reduction and if possible removal an important subject in the development of tomographic image reconstruction algorithms. Limited angle artifacts are caused by the limited angular measurements, constraining the available tomographic information. This thesis focuses on reducing these artifacts via image reconstruction in two cases of incomplete measurements from: (1) the gaps left after the removal of high density objects such as dental fillings, screws and implants in computed tomography (CT) and (2) partial ring scanner configurations in positron emission tomography (PET). In order to include knowledge about the measurement and noise, prior terms were used within the reconstruction methods. Careful consideration was given to the trade-off between image blurring and noise reduction upon reconstruction of low-dose measurements.Development of reconstruction methods is an incremental process starting with testing on simple phantoms towards more clinically relevant ones by modeling the respective physical processes involved. In this work, phantoms were constructed to ensure that the proposed reconstruction methods addressed to the limited angle problem. The reconstructed images were assessed qualitatively and quantitatively in terms of noise reduction, edge sharpness and contrast recovery.Maximum a posteriori (MAP) estimation with median root prior (MRP) was selected for the reconstruction of limited angle measurements. MAP with MRP successfully reduced the artifacts caused by limited angle data in various datasets, tested with the reconstruction of both list-mode and projection data. In all cases, its performance was found to be superior to conventional reconstruction methods such as total-variation (TV) prior, maximum likelihood expectation maximization (MLEM) and filtered backprojection (FBP). MAP with MRP was also more robust with respect to parameter selection than MAP with TV prior.This thesis demonstrates the wide-range applicability of MAP with MRP in medical tomography, especially in low-dose imaging. Furthermore, we emphasize the importance of developing and testing reconstruction methods with application-specific phantoms, together with the properties and limitations of the measurements in mind

Multiplexed fluorescence diffuse optical tomography

Fluorescence tomography (FT) is an emerging non-invasive in vivo molecular imaging modality that aims at quantification and three-dimensional (3D) localization of fluorescent tagged inclusions, such as cancer lesions and drug molecules, buried deep in human and animal subjects. Depth-resolved 3D reconstruction of fluorescent inclusions distributed over the volume of optically turbid biological tissue using the diffuse fluorescent photons detected on the skin poses a highly ill-conditioned problem, as depth information must be extracted from boundary data. Due to this ill-posed nature of FT reconstructions, noise and errors in the data can severely impair the accuracy of the 3D reconstructions. Consequently, improvements in the signal-to-noise ratio (SNR) of the data significantly enhance the quality of the FT reconstructions. Furthermore, enhancing the SNR of the FT data can greatly contribute to the speed of FT scans. The pivotal factor in the SNR of the FT data is the power of the radiation illuminating the subject and exciting the administered fluorescent agents. In existing single-point illumination FT systems, the illumination power level is limited by the skin maximum radiation exposure levels. In this research, a multiplexed architecture governed by the Hadamard transform was conceptualized, developed, and experimentally implemented for orders-of-magnitude enhancement of the SNR and the robustness of FT reconstructions. The multiplexed FT system allows for Hadamard-coded multi-point illumination of the subject while maintaining the maximal information content of the FT data. The significant improvements offered by the multiplexed FT system were validated by numerical and experimental studies carried out using a custom-built multiplexed FT system developed exclusively in this work. The studies indicate that Hadamard multiplexing offers significantly enhanced robustness in reconstructing deep fluorescent inclusions from low-SNR FT data.Ph.D

Improved Modeling and Image Generation for Fluorescence Molecular Tomography (FMT) and Positron Emission Tomography (PET)

In this thesis, we aim to improve quantitative medical imaging with advanced image generation algorithms. We focus on two specific imaging modalities: fluorescence molecular tomography (FMT) and positron emission tomography (PET). For FMT, we present a novel photon propagation model for its forward model, and in addition, we propose and investigate a reconstruction algorithm for its inverse problem. In the first part, we develop a novel Neumann-series-based radiative transfer equation (RTE) that incorporates reflection boundary conditions in the model. In addition, we propose a novel reconstruction technique for diffuse optical imaging that incorporates this Neumann-series-based RTE as forward model. The proposed model is assessed using a simulated 3D diffuse optical imaging setup, and the results demonstrate the importance of considering photon reflection at boundaries when performing photon propagation modeling. In the second part, we propose a statistical reconstruction algorithm for FMT. The algorithm is based on sparsity-initialized maximum-likelihood expectation maximization (MLEM), taking into account the Poisson nature of data in FMT and the sparse nature of images. The proposed method is compared with a pure sparse reconstruction method as well as a uniform-initialized MLEM reconstruction method. Results indicate the proposed method is more robust to noise and shows improved qualitative and quantitative performance. For PET, we present an MRI-guided partial volume correction algorithm for brain imaging, aiming to recover qualitative and quantitative loss due to the limited resolution of PET system, while keeping image noise at a low level. The proposed method is based on an iterative deconvolution model with regularization using parallel level sets. A non-smooth optimization algorithm is developed so that the proposed method can be feasibly applied for 3D images and avoid additional blurring caused by conventional smooth optimization process. We evaluate the proposed method using both simulation data and in vivo human data collected from the Baltimore Longitudinal Study of Aging (BLSA). Our proposed method is shown to generate images with reduced noise and improved structure details, as well as increased number of statistically significant voxels in study of aging. Results demonstrate our method has promise to provide superior performance in clinical imaging scenarios

Numerical Approaches for Solving the Combined Reconstruction and Registration of Digital Breast Tomosynthesis

Heavy demands on the development of medical imaging modalities for breast cancer detection have been witnessed in the last three decades in an attempt to reduce the mortality associated with the disease. Recently, Digital Breast Tomosynthesis (DBT) shows its promising in the early diagnosis when lesions are small. In particular, it offers potential benefits over X-ray mammography - the current modality of choice for breast screening - of increased sensitivity and specificity for comparable X-ray dose, speed, and cost. An important feature of DBT is that it provides a pseudo-3D image of the breast. This is of particular relevance for heterogeneous dense breasts of young women, which can inhibit detection of cancer using conventional mammography. In the same way that it is difficult to see a bird from the edge of the forest, detecting cancer in a conventional 2D mammogram is a challenging task. Three-dimensional DBT, however, enables us to step through the forest, i.e., the breast, reducing the confounding effect of superimposed tissue and so (potentially) increasing the sensitivity and specificity of cancer detection. The workflow in which DBT would be used clinically, involves two key tasks: reconstruction, to generate a 3D image of the breast, and registration, to enable images from different visits to be compared as is routinely performed by radiologists working with conventional mammograms. Conventional approaches proposed in the literature separate these steps, solving each task independently. This can be effective if reconstructing using a complete set of data. However, for ill-posed limited-angle problems such as DBT, estimating the deformation is difficult because of the significant artefacts associated with DBT reconstructions, leading to severe inaccuracies in the registration. The aim of my work is to find and evaluate methods capable of allying these two tasks, which will enhance the performance of each process as a result. Consequently, I prove that the processes of reconstruction and registration of DBT are not independent but reciprocal. This thesis proposes innovative numerical approaches combining reconstruction of a pair of temporal DBT acquisitions with their registration iteratively and simultaneously. To evaluate the performance of my methods I use synthetic images, breast MRI, and DBT simulations with in-vivo breast compressions. I show that, compared to the conventional sequential method, jointly estimating image intensities and transformation parameters gives superior results with respect to both reconstruction fidelity and registration accuracy

Signal sampling and processing in magnetic resonance applications

In this thesis, signal sampling and processing techniques are developed for magnetic resonance applications, to improve the estimation of magnetic resonance parameters and to reduce experimental acquisition times. Two processing techniques are developed for Nuclear Magnetic Resonance (NMR) relaxation and diffusion experiments. $L_1$ regularization is recommended for extracting parameter distributions which are a priori known to be composed of sparse features. $L_1$ regularization is shown to be stable at signal-to-noise ratios < 20 and capable of resolving relaxation time constants and diffusion coefficients which differ by as little as 10%, such as in relaxation and diffusion studies of hexane/dodecane in porous media. Modified Total Generalized Variation (MTGV) regularization is recommended for extracting parameter distributions for which there is no prior knowledge of whether they are composed of sparse or smooth features. MTGV regularization is shown to perform better than conventional processing techniques and $L_1$ regularization over a range of simulated distributions. A method for optimising sampling patterns for relaxation and diffusion experiments, based on the Cramér-Rao Lower Bound theory, is presented. The method is validated against pulsed field gradient NMR diffusion data of two experimental systems. In the first experimental system, the sampling pattern is optimised for the most accurate estimation of the lognormal distribution parameters of an emulsion droplet size distribution of toluene in water. In the second experimental system, the sampling pattern is optimised for the most accurate estimation of the bi-exponential model parameters of a binary mixture of methane/ethane adsorbed in a zeolite. The proposed method predicts an uncertainty in estimating the model parameters which is < 10% different from the uncertainty estimated from the experimental data sampled using the same sampling pattern. Signal sampling and processing techniques are subsequently combined to reduce experimental acquisition times, which opens opportunities for studying unsteady systems over a long acquisition time and investigating fast-changing phenomena. A 32-fold decrease in the experimental acquisition time is achieved in extracting 3D spatially resolved spin spin relaxation maps. This is expected to be useful in investigating porous media systems. Three-component velocity maps on a 2D image, acquired every 4 ms, are used to capture, for the first time, the hydrodynamics of a bubble burst event. The experimental data are used to validate the predictions of numerical works.Gates Cambridge Scholarshi

Abstracts on Radio Direction Finding (1899 - 1995)

The files on this record represent the various databases that originally composed the CD-ROM issue of "Abstracts on Radio Direction Finding" database, which is now part of the Dudley Knox Library's Abstracts and Selected Full Text Documents on Radio Direction Finding (1899 - 1995) Collection. (See Calhoun record https://calhoun.nps.edu/handle/10945/57364 for further information on this collection and the bibliography). Due to issues of technological obsolescence preventing current and future audiences from accessing the bibliography, DKL exported and converted into the three files on this record the various databases contained in the CD-ROM. The contents of these files are: 1) RDFA_CompleteBibliography_xls.zip [RDFA_CompleteBibliography.xls: Metadata for the complete bibliography, in Excel 97-2003 Workbook format; RDFA_Glossary.xls: Glossary of terms, in Excel 97-2003 Workbookformat; RDFA_Biographies.xls: Biographies of leading figures, in Excel 97-2003 Workbook format]; 2) RDFA_CompleteBibliography_csv.zip [RDFA_CompleteBibliography.TXT: Metadata for the complete bibliography, in CSV format; RDFA_Glossary.TXT: Glossary of terms, in CSV format; RDFA_Biographies.TXT: Biographies of leading figures, in CSV format]; 3) RDFA_CompleteBibliography.pdf: A human readable display of the bibliographic data, as a means of double-checking any possible deviations due to conversion