Optimizing the number of equivalent iterations of 3D OSEM in SPECT reconstruction of I-131 focal activities

International Nuclear Information System (INIS)

Koral, Kenneth F.; Kritzmaan, James N.; Rogers, Virginia E.; Ackermann, Robert J.; A Fessler, Jeffrey

2007-01-01

To externally estimate the radiation dose to a tumor during therapy with I-131 radiopharmaceuticals, and its distribution, one must accurately estimate activity, and its distribution, by means of SPECT imaging. Our objective is to characterize the quantification of the total activity in focal targets and in their uniform background, and of the activity distribution within the targets, after 3D Ordered Subsets Expectation Maximization (OSEM) reconstruction with attenuation and scatter correction and no post smoothing, in the good-counting-statistics case. A cylindrical phantom containing seven spheres simulating tumors, ranging in volume from 209 to 4.2 cm 3 , and filled with an I-131 water solution containing background, was imaged. A Siemens Symbia SPECT/CT scanner was used to acquire 128x128 projection images, employing 60 angles over 360 o . With dynamic SPECT, 10 sequential acquisitions of 15 min duration each were obtained and each was reconstructed with particular values of the number of subsets and the number of iterations. Let the product of the number of subsets times the number of iterations equal the equivalent number of iterations, EI. The counts-to-activity conversion factor was derived from the average ratio of total count divided by true activity for the largest sphere at the largest value of EI. Then, for the activity of each sphere at each of the values of EI, we evaluated (1) the fractional variance (variance in estimate over true activity), (2) the fractional bias (average estimate bias over true activity) and (3) the fractional error (the root mean square error (RMSE) in the estimate divided by the true activity). The fractional bias and fractional variance were smaller for the larger spheres compared to the smaller (the fractional bias decreased faster with an increase in the fractional variance for them as well). The RMSE was dominated by the bias. The fractional error decreased as EI increased for all sphere sizes. The minimum average value

Studies oriented to optimize the image quality of the small animal PET: Clear PET, modifying some of the parameters of the reconstruction algorithm IMF-OSEM 3D on the data acquisition simulated with GAMOS; Estudios para la optimizaciOn de la calidad de imagen en el escaner ClearPET, modifi cando parametros del algoritmo IMF-OSEM 3D sobre adquisiciones simuladas con GAMOS

Energy Technology Data Exchange (ETDEWEB)

Canadas, M.; Mendoza, J.; Embid, M.

2007-09-27

This report presents studies oriented to optimize the image quality of the small animal PET: Clear- PET. Certain figures of merit (FOM) were used to assess a quantitative value of the contrast and delectability of lesions. The optimization was carried out modifying some of the parameters in the reconstruction software of the scanner, imaging a mini-Derenzo phantom and a cylinder phantom with background activity and two hot spheres. Specifically, it was evaluated the incidence of the inter-update Metz filter (IMF) inside the iterative reconstruction algorithm 3D OSEM. The data acquisition was simulated using the GAMOS framework (Monte Carlo simulation). Integrating GAMOS output with the reconstruction software of the scanner was an additional novelty of this work, to achieve this, data sets were written with the list-mode format (LMF) of ClearPET. In order to verify the optimum values obtained, we foresee to make real acquisitions in the ClearPET of CIEMAT. (Author) 17 refs.

Accelerating image reconstruction in dual-head PET system by GPU and symmetry properties.

Directory of Open Access Journals (Sweden)

Cheng-Ying Chou

Full Text Available Positron emission tomography (PET is an important imaging modality in both clinical usage and research studies. We have developed a compact high-sensitivity PET system that consisted of two large-area panel PET detector heads, which produce more than 224 million lines of response and thus request dramatic computational demands. In this work, we employed a state-of-the-art graphics processing unit (GPU, NVIDIA Tesla C2070, to yield an efficient reconstruction process. Our approaches ingeniously integrate the distinguished features of the symmetry properties of the imaging system and GPU architectures, including block/warp/thread assignments and effective memory usage, to accelerate the computations for ordered subset expectation maximization (OSEM image reconstruction. The OSEM reconstruction algorithms were implemented employing both CPU-based and GPU-based codes, and their computational performance was quantitatively analyzed and compared. The results showed that the GPU-accelerated scheme can drastically reduce the reconstruction time and thus can largely expand the applicability of the dual-head PET system.

High-resolution 3D-GRE imaging of the abdomen using controlled aliasing acceleration technique - a feasibility study

Energy Technology Data Exchange (ETDEWEB)

AlObaidy, Mamdoh; Ramalho, Miguel; Busireddy, Kiran K.R.; Liu, Baodong; Burke, Lauren M.; Altun, Ersan; Semelka, Richard C. [University of North Carolina at Chapel Hill, Department of Radiology, Chapel Hill, NC (United States); Dale, Brian M. [Siemens Medical Solutions, MR Research and Development, Morrisville, NC (United States)

2015-12-15

To assess the feasibility of high-resolution 3D-gradient-recalled echo (GRE) fat-suppressed T1-weighted images using controlled aliasing acceleration technique (CAIPIRINHA-VIBE), and compare image quality and lesion detection to standard-resolution 3D-GRE images using conventional acceleration technique (GRAPPA-VIBE). Eighty-four patients (41 males, 43 females; age range: 14-90 years, 58.8 ± 15.6 years) underwent abdominal MRI at 1.5 T with CAIPIRINHA-VIBE [spatial resolution, 0.76 ± 0.04 mm] and GRAPPA-VIBE [spatial resolution, 1.17 ± 0.14 mm]. Two readers independently reviewed image quality, presence of artefacts, lesion conspicuity, and lesion detection. Kappa statistic was used to assess interobserver agreement. Wilcoxon signed-rank test was used for image qualitative pairwise comparisons. Logistic regression with post-hoc testing was used to evaluate statistical significance of lesions evaluation. Interobserver agreement ranged between 0.45-0.93. Pre-contrast CAIPIRINHA-VIBE showed significantly (p < 0.001) sharper images and lesion conspicuity with decreased residual aliasing, but more noise enhancement and inferior image quality. Post-contrast CAIPIRINHA-VIBE showed significantly (p < 0.001) sharper images and higher lesion conspicuity, with less respiratory motion and residual aliasing artefacts. Inferior fat-suppression was noticeable on CAIPIRINHA-VIBE sequences (p < 0.001). High in-plane resolution abdominal 3D-GRE fat-suppressed T1-weighted imaging using controlled-aliasing acceleration technique is feasible and yields sharper images compared to standard-resolution images using standard acceleration, with higher post-contrast image quality and trend for improved hepatic lesions detection. (orig.)

High-resolution 3D-GRE imaging of the abdomen using controlled aliasing acceleration technique - a feasibility study

International Nuclear Information System (INIS)

AlObaidy, Mamdoh; Ramalho, Miguel; Busireddy, Kiran K.R.; Liu, Baodong; Burke, Lauren M.; Altun, Ersan; Semelka, Richard C.; Dale, Brian M.

2015-01-01

To assess the feasibility of high-resolution 3D-gradient-recalled echo (GRE) fat-suppressed T1-weighted images using controlled aliasing acceleration technique (CAIPIRINHA-VIBE), and compare image quality and lesion detection to standard-resolution 3D-GRE images using conventional acceleration technique (GRAPPA-VIBE). Eighty-four patients (41 males, 43 females; age range: 14-90 years, 58.8 ± 15.6 years) underwent abdominal MRI at 1.5 T with CAIPIRINHA-VIBE [spatial resolution, 0.76 ± 0.04 mm] and GRAPPA-VIBE [spatial resolution, 1.17 ± 0.14 mm]. Two readers independently reviewed image quality, presence of artefacts, lesion conspicuity, and lesion detection. Kappa statistic was used to assess interobserver agreement. Wilcoxon signed-rank test was used for image qualitative pairwise comparisons. Logistic regression with post-hoc testing was used to evaluate statistical significance of lesions evaluation. Interobserver agreement ranged between 0.45-0.93. Pre-contrast CAIPIRINHA-VIBE showed significantly (p < 0.001) sharper images and lesion conspicuity with decreased residual aliasing, but more noise enhancement and inferior image quality. Post-contrast CAIPIRINHA-VIBE showed significantly (p < 0.001) sharper images and higher lesion conspicuity, with less respiratory motion and residual aliasing artefacts. Inferior fat-suppression was noticeable on CAIPIRINHA-VIBE sequences (p < 0.001). High in-plane resolution abdominal 3D-GRE fat-suppressed T1-weighted imaging using controlled-aliasing acceleration technique is feasible and yields sharper images compared to standard-resolution images using standard acceleration, with higher post-contrast image quality and trend for improved hepatic lesions detection. (orig.)

Studies oriented to optimize the image quality of the small animal PET: Clear PET, modifying some of the parameters of the reconstruction algorithm IMF-OSEM 3D on the data acquisition simulated with GAMOS

International Nuclear Information System (INIS)

Canadas, M.; Mendoza, J.; Embid, M.

2007-01-01

This report presents studies oriented to optimize the image quality of the small animal PET: Clear- PET. Certain figures of merit (FOM) were used to assess a quantitative value of the contrast and delectability of lesions. The optimization was carried out modifying some of the parameters in the reconstruction software of the scanner, imaging a mini-Derenzo phantom and a cylinder phantom with background activity and two hot spheres. Specifically, it was evaluated the incidence of the inter-update Metz filter (IMF) inside the iterative reconstruction algorithm 3D OSEM. The data acquisition was simulated using the GAMOS framework (Monte Carlo simulation). Integrating GAMOS output with the reconstruction software of the scanner was an additional novelty of this work, to achieve this, data sets were written with the list-mode format (LMF) of ClearPET. In order to verify the optimum values obtained, we foresee to make real acquisitions in the ClearPET of CIEMAT. (Author) 17 refs

Characterization and simulation of noise in PET images reconstructed with OSEM: Development of a method for the generation of synthetic images.

Science.gov (United States)

Castro, P; Huerga, C; Chamorro, P; Garayoa, J; Roch, M; Pérez, L

2018-04-17

The goals of the study are to characterize imaging properties in 2D PET images reconstructed with the iterative algorithm ordered-subset expectation maximization (OSEM) and to propose a new method for the generation of synthetic images. The noise is analyzed in terms of its magnitude, spatial correlation, and spectral distribution through standard deviation, autocorrelation function, and noise power spectrum (NPS), respectively. Their variations with position and activity level are also analyzed. This noise analysis is based on phantom images acquired from 18 F uniform distributions. Experimental recovery coefficients of hot spheres in different backgrounds are employed to study the spatial resolution of the system through point spread function (PSF). The NPS and PSF functions provide the baseline for the proposed simulation method: convolution with PSF as kernel and noise addition from NPS. The noise spectral analysis shows that the main contribution is of random nature. It is also proven that attenuation correction does not alter noise texture but it modifies its magnitude. Finally, synthetic images of 2 phantoms, one of them an anatomical brain, are quantitatively compared with experimental images showing a good agreement in terms of pixel values and pixel correlations. Thus, the contrast to noise ratio for the biggest sphere in the NEMA IEC phantom is 10.7 for the synthetic image and 8.8 for the experimental image. The properties of the analyzed OSEM-PET images can be described by NPS and PSF functions. Synthetic images, even anatomical ones, are successfully generated by the proposed method based on the NPS and PSF. Copyright © 2018 Sociedad Española de Medicina Nuclear e Imagen Molecular. Publicado por Elsevier España, S.L.U. All rights reserved.

LOR-OSEM: statistical PET reconstruction from raw line-of-response histograms

International Nuclear Information System (INIS)

Kadrmas, Dan J

2004-01-01

Iterative statistical reconstruction methods are becoming the standard in positron emission tomography (PET). Conventional maximum-likelihood expectation-maximization (MLEM) and ordered-subsets (OSEM) algorithms act on data which have been pre-processed into corrected, evenly-spaced histograms; however, such pre-processing corrupts the Poisson statistics. Recent advances have incorporated attenuation, scatter and randoms compensation into the iterative reconstruction. The objective of this work was to incorporate the remaining pre-processing steps, including arc correction, to reconstruct directly from raw unevenly-spaced line-of-response (LOR) histograms. This exactly preserves Poisson statistics and full spatial information in a manner closely related to listmode ML, making full use of the ML statistical model. The LOR-OSEM algorithm was implemented using a rotation-based projector which maps directly to the unevenly-spaced LOR grid. Simulation and phantom experiments were performed to characterize resolution, contrast and noise properties for 2D PET. LOR-OSEM provided a beneficial noise-resolution tradeoff, outperforming AW-OSEM by about the same margin that AW-OSEM outperformed pre-corrected OSEM. The relationship between LOR-ML and listmode ML algorithms was explored, and implementation differences are discussed. LOR-OSEM is a viable alternative to AW-OSEM for histogram-based reconstruction with improved spatial resolution and noise properties

Clinical evaluation of 2D versus 3D whole-body PET image quality using a dedicated BGO PET scanner

International Nuclear Information System (INIS)

Visvikis, D.; Griffiths, D.; Costa, D.C.; Bomanji, J.; Ell, P.J.

2005-01-01

Three-dimensional positron emission tomography (3D PET) results in higher system sensitivity, with an associated increase in the detection of scatter and random coincidences. The objective of this work was to compare, from a clinical perspective, 3D and two-dimensional (2D) acquisitions in terms of whole-body (WB) PET image quality with a dedicated BGO PET system. 2D and 3D WB emission acquisitions were carried out in 70 patients. Variable acquisition parameters in terms of time of emission acquisition per axial field of view (aFOV) and slice overlap between sequential aFOVs were used during the 3D acquisitions. 3D and 2D images were reconstructed using FORE+WLS and OSEM respectively. Scatter correction was performed by convolution subtraction and a model-based scatter correction in 2D and 3D respectively. All WB images were attenuation corrected using segmented transmission scans. Images were blindly assessed by three observers for the presence of artefacts, confidence in lesion detection and overall image quality using a scoring system. Statistically significant differences between 2D and 3D image quality were only obtained for 3D emission acquisitions of 3 min. No statistically significant differences were observed for image artefacts or lesion detectability scores. Image quality correlated significantly with patient weight for both modes of operation. Finally, no differences were seen in image artefact scores for the different axial slice overlaps considered, suggesting the use of five slice overlaps in 3D WB acquisitions. 3D WB imaging using a dedicated BGO-based PET scanner offers similar image quality to that obtained in 2D considering similar overall times of acquisitions. (orig.)

Image based cardiac acceleration map using statistical shape and 3D+t myocardial tracking models; in-vitro study on heart phantom

Science.gov (United States)

Pashaei, Ali; Piella, Gemma; Planes, Xavier; Duchateau, Nicolas; de Caralt, Teresa M.; Sitges, Marta; Frangi, Alejandro F.

2013-03-01

It has been demonstrated that the acceleration signal has potential to monitor heart function and adaptively optimize Cardiac Resynchronization Therapy (CRT) systems. In this paper, we propose a non-invasive method for computing myocardial acceleration from 3D echocardiographic sequences. Displacement of the myocardium was estimated using a two-step approach: (1) 3D automatic segmentation of the myocardium at end-diastole using 3D Active Shape Models (ASM); (2) propagation of this segmentation along the sequence using non-rigid 3D+t image registration (temporal di eomorphic free-form-deformation, TDFFD). Acceleration was obtained locally at each point of the myocardium from local displacement. The framework has been tested on images from a realistic physical heart phantom (DHP-01, Shelley Medical Imaging Technologies, London, ON, CA) in which the displacement of some control regions was known. Good correlation has been demonstrated between the estimated displacement function from the algorithms and the phantom setup. Due to the limited temporal resolution, the acceleration signals are sparse and highly noisy. The study suggests a non-invasive technique to measure the cardiac acceleration that may be used to improve the monitoring of cardiac mechanics and optimization of CRT.

GPU accelerated generation of digitally reconstructed radiographs for 2-D/3-D image registration.

Science.gov (United States)

Dorgham, Osama M; Laycock, Stephen D; Fisher, Mark H

2012-09-01

Recent advances in programming languages for graphics processing units (GPUs) provide developers with a convenient way of implementing applications which can be executed on the CPU and GPU interchangeably. GPUs are becoming relatively cheap, powerful, and widely available hardware components, which can be used to perform intensive calculations. The last decade of hardware performance developments shows that GPU-based computation is progressing significantly faster than CPU-based computation, particularly if one considers the execution of highly parallelisable algorithms. Future predictions illustrate that this trend is likely to continue. In this paper, we introduce a way of accelerating 2-D/3-D image registration by developing a hybrid system which executes on the CPU and utilizes the GPU for parallelizing the generation of digitally reconstructed radiographs (DRRs). Based on the advancements of the GPU over the CPU, it is timely to exploit the benefits of many-core GPU technology by developing algorithms for DRR generation. Although some previous work has investigated the rendering of DRRs using the GPU, this paper investigates approximations which reduce the computational overhead while still maintaining a quality consistent with that needed for 2-D/3-D registration with sufficient accuracy to be clinically acceptable in certain applications of radiation oncology. Furthermore, by comparing implementations of 2-D/3-D registration on the CPU and GPU, we investigate current performance and propose an optimal framework for PC implementations addressing the rigid registration problem. Using this framework, we are able to render DRR images from a 256×256×133 CT volume in ~24 ms using an NVidia GeForce 8800 GTX and in ~2 ms using NVidia GeForce GTX 580. In addition to applications requiring fast automatic patient setup, these levels of performance suggest image-guided radiation therapy at video frame rates is technically feasible using relatively low cost PC

GPU-accelerated denoising of 3D magnetic resonance images

Energy Technology Data Exchange (ETDEWEB)

Howison, Mark; Wes Bethel, E.

2014-05-29

The raw computational power of GPU accelerators enables fast denoising of 3D MR images using bilateral filtering, anisotropic diffusion, and non-local means. In practice, applying these filtering operations requires setting multiple parameters. This study was designed to provide better guidance to practitioners for choosing the most appropriate parameters by answering two questions: what parameters yield the best denoising results in practice? And what tuning is necessary to achieve optimal performance on a modern GPU? To answer the first question, we use two different metrics, mean squared error (MSE) and mean structural similarity (MSSIM), to compare denoising quality against a reference image. Surprisingly, the best improvement in structural similarity with the bilateral filter is achieved with a small stencil size that lies within the range of real-time execution on an NVIDIA Tesla M2050 GPU. Moreover, inappropriate choices for parameters, especially scaling parameters, can yield very poor denoising performance. To answer the second question, we perform an autotuning study to empirically determine optimal memory tiling on the GPU. The variation in these results suggests that such tuning is an essential step in achieving real-time performance. These results have important implications for the real-time application of denoising to MR images in clinical settings that require fast turn-around times.

Application of the OS-EM method to the 123I-IMP ARG method. Comparison between FBP and OS-EM methods

International Nuclear Information System (INIS)

Sasaki, Kazuko; Satou, Ayuko; Oomura, Tomomi; Ono, Madoka; Hachiya, Takenori

2004-01-01

We investigated application of the ordered subsets-expectation maximization (OS-EM) method to the 123 I-IMP autoradiography (ARG) method to measure regional cerebral blood flow (rCBF). First, scan time and subsets were fixed at 20 min and 16, respectively, and the influence of iteration on the cross calibration factor (CCF) and quantitative rCBF values obtained by the ARG method was investigated when the iteration number was set at 2, 4, 8, 16, 32, and 90. Next, with the number of iterations set at 4, we compared the scanning times of OS-EM and filtered back projection (FBP). We determined that the CCF values remained at the same level irrespective of iteration number. Quantitative rCBF values had no association with iteration number, either. Using the quantitative rCBF values obtained by 20-min. scanning with FBP as a standard, the time period for collecting SPECT data was 10 min, without sacrificing image quality or quantification. Quantitative rCBF obtained by OS-EM was estimated to be higher than that by FBP. (author)

Comparing an accelerated 3D fast spin-echo sequence (CS-SPACE) for knee 3-T magnetic resonance imaging with traditional 3D fast spin-echo (SPACE) and routine 2D sequences

Energy Technology Data Exchange (ETDEWEB)

Altahawi, Faysal F.; Blount, Kevin J.; Omar, Imran M. [Northwestern University Feinberg School of Medicine, Department of Radiology, Chicago, IL (United States); Morley, Nicholas P. [Marshfield Clinic, Department of Radiology, Marshfield, WI (United States); Raithel, Esther [Siemens Healthcare GmbH, Erlangen (Germany)

2017-01-15

To compare a faster, new, high-resolution accelerated 3D-fast-spin-echo (3D-FSE) acquisition sequence (CS-SPACE) to traditional 2D and high-resolution 3D sequences for knee 3-T magnetic resonance imaging (MRI). Twenty patients received knee MRIs that included routine 2D (T1, PD ± FS, T2-FS; 0.5 x 0.5 x 3 mm{sup 3}; ∝10 min), traditional 3D FSE (SPACE-PD-FS; 0.5 x 0.5 x 0.5 mm{sup 3}; ∝7.5 min), and accelerated 3D-FSE prototype (CS-SPACE-PD-FS; 0.5 x 0.5 x 0.5 mm{sup 3}; ∝5 min) acquisitions on a 3-T MRI system (Siemens MAGNETOM Skyra). Three musculoskeletal radiologists (MSKRs) prospectively and independently reviewed the studies with graded surveys comparing image and diagnostic quality. Tissue-specific signal-to-noise ratios (SNR) and contrast-to-noise ratios (CNR) were also compared. MSKR-perceived diagnostic quality of cartilage was significantly higher for CS-SPACE than for SPACE and 2D sequences (p < 0.001). Assessment of diagnostic quality of menisci and synovial fluid was higher for CS-SPACE than for SPACE (p < 0.001). CS-SPACE was not significantly different from SPACE but had lower assessments than 2D sequences for evaluation of bones, ligaments, muscles, and fat (p ≤ 0.004). 3D sequences had higher spatial resolution, but lower overall assessed contrast (p < 0.001). Overall image quality from CS-SPACE was assessed as higher than SPACE (p = 0.007), but lower than 2D sequences (p < 0.001). Compared to SPACE, CS-SPACE had higher fluid SNR and CNR against all other tissues (all p < 0.001). The CS-SPACE prototype allows for faster isotropic acquisitions of knee MRIs over currently used protocols. High fluid-to-cartilage CNR and higher spatial resolution over routine 2D sequences may present a valuable role for CS-SPACE in the evaluation of cartilage and menisci. (orig.)

Influence of the partial volume correction method on (18)F-fluorodeoxyglucose brain kinetic modelling from dynamic PET images reconstructed with resolution model based OSEM.

Science.gov (United States)

Bowen, Spencer L; Byars, Larry G; Michel, Christian J; Chonde, Daniel B; Catana, Ciprian

2013-10-21

Kinetic parameters estimated from dynamic (18)F-fluorodeoxyglucose ((18)F-FDG) PET acquisitions have been used frequently to assess brain function in humans. Neglecting partial volume correction (PVC) for a dynamic series has been shown to produce significant bias in model estimates. Accurate PVC requires a space-variant model describing the reconstructed image spatial point spread function (PSF) that accounts for resolution limitations, including non-uniformities across the field of view due to the parallax effect. For ordered subsets expectation maximization (OSEM), image resolution convergence is local and influenced significantly by the number of iterations, the count density, and background-to-target ratio. As both count density and background-to-target values for a brain structure can change during a dynamic scan, the local image resolution may also concurrently vary. When PVC is applied post-reconstruction the kinetic parameter estimates may be biased when neglecting the frame-dependent resolution. We explored the influence of the PVC method and implementation on kinetic parameters estimated by fitting (18)F-FDG dynamic data acquired on a dedicated brain PET scanner and reconstructed with and without PSF modelling in the OSEM algorithm. The performance of several PVC algorithms was quantified with a phantom experiment, an anthropomorphic Monte Carlo simulation, and a patient scan. Using the last frame reconstructed image only for regional spread function (RSF) generation, as opposed to computing RSFs for each frame independently, and applying perturbation geometric transfer matrix PVC with PSF based OSEM produced the lowest magnitude bias kinetic parameter estimates in most instances, although at the cost of increased noise compared to the PVC methods utilizing conventional OSEM. Use of the last frame RSFs for PVC with no PSF modelling in the OSEM algorithm produced the lowest bias in cerebral metabolic rate of glucose estimates, although by less than 5% in

Geometrically undistorted MRI in the presence of field inhomogeneities using compressed sensing accelerated broadband 3D phase encoded turbo spin-echo imaging

International Nuclear Information System (INIS)

Van Gorp, Jetse S; Bakker, Chris J G; Bouwman, Job G; Zijlstra, Frank; Seevinck, Peter R; Smink, Jouke

2015-01-01

In this study, we explore the potential of compressed sensing (CS) accelerated broadband 3D phase-encoded turbo spin-echo (3D-PE-TSE) for the purpose of geometrically undistorted imaging in the presence of field inhomogeneities. To achieve this goal 3D-PE-SE and 3D-PE-TSE sequences with broadband rf pulses and dedicated undersampling patterns were implemented on a clinical scanner. Additionally, a 3D multi-spectral spin-echo (ms3D-SE) sequence was implemented for reference purposes. First, we demonstrated the influence of susceptibility induced off-resonance effects on the spatial encoding of broadband 3D-SE, ms3D-SE, 3D-PE-SE and 3D-PE-TSE using a grid phantom containing a titanium implant (Δχ = 182 ppm) with x-ray CT as a gold standard. These experiments showed that the spatial encoding of 3D-PE-(T)SE was unaffected by susceptibility induced off-resonance effects, which caused geometrical distortions and/or signal hyper-intensities in broadband 3D-SE and, to a lesser extent, in ms3D-SE frequency encoded methods. Additionally, an SNR analysis was performed and the temporally resolved signal of 3D-PE-(T)SE sequences was exploited to retrospectively decrease the acquisition bandwidth and obtain field offset maps. The feasibility of CS acceleration was studied retrospectively and prospectively for the 3D-PE-SE sequence using an existing CS algorithm adapted for the reconstruction of 3D data with undersampling in all three phase encoded dimensions. CS was combined with turbo-acceleration by variable density undersampling and spherical stepwise T 2 weighting by randomly sorting consecutive echoes in predefined spherical k-space layers. The CS-TSE combination resulted in an overall acceleration factor of 60, decreasing the original 3D-PE-SE scan time from 7 h to 7 min. Finally, CS accelerated 3D-PE-TSE in vivo images of a titanium screw were obtained within 10 min using a micro-coil demonstrating the feasibility of geometrically undistorted MRI near severe

Whole-brain background-suppressed pCASL MRI with 1D-accelerated 3D RARE Stack-Of-Spirals readout.

Directory of Open Access Journals (Sweden)

Marta Vidorreta

Full Text Available Arterial Spin Labeled (ASL perfusion MRI enables non-invasive, quantitative measurements of tissue perfusion, and has a broad range of applications including brain functional imaging. However, ASL suffers from low signal-to-noise ratio (SNR, limiting image resolution. Acquisitions using 3D readouts are optimal for background-suppression of static signals, but can be SAR intensive and typically suffer from through-plane blurring. In this study, we investigated the use of accelerated 3D readouts to obtain whole-brain, high-SNR ASL perfusion maps and reduce SAR deposition. Parallel imaging was implemented along the partition-encoding direction in a pseudo-continuous ASL sequence with background-suppression and 3D RARE Stack-Of-Spirals readout, and its performance was evaluated in three small cohorts. First, both non-accelerated and two-fold accelerated single-shot versions of the sequence were evaluated in healthy volunteers during a motor-photic task, and the performance was compared in terms of temporal SNR, GM-WM contrast, and statistical significance of the detected activation. Secondly, single-shot 1D-accelerated imaging was compared to a two-shot accelerated version to assess benefits of SNR and spatial resolution for applications in which temporal resolution is not paramount. Third, the efficacy of this approach in clinical populations was assessed by applying the single-shot 1D-accelerated version to a larger cohort of elderly volunteers. Accelerated data demonstrated the ability to detect functional activation at the subject level, including cerebellar activity, without loss in the perfusion signal temporal stability and the statistical power of the activations. The use of acceleration also resulted in increased GM-WM contrast, likely due to reduced through-plane partial volume effects, that were further attenuated with the use of two-shot readouts. In a clinical cohort, image quality remained excellent, and expected effects of age and sex

3D IBFV : Hardware-Accelerated 3D Flow Visualization

NARCIS (Netherlands)

Telea, Alexandru; Wijk, Jarke J. van

2003-01-01

We present a hardware-accelerated method for visualizing 3D flow fields. The method is based on insertion, advection, and decay of dye. To this aim, we extend the texture-based IBFV technique for 2D flow visualization in two main directions. First, we decompose the 3D flow visualization problem in a

3D IBFV : hardware-accelerated 3D flow visualization

NARCIS (Netherlands)

Telea, A.C.; Wijk, van J.J.

2003-01-01

We present a hardware-accelerated method for visualizing 3D flow fields. The method is based on insertion, advection, and decay of dye. To this aim, we extend the texture-based IBFV technique presented by van Wijk (2001) for 2D flow visualization in two main directions. First, we decompose the 3D

A new reconstruction strategy for image improvement in pinhole SPECT

International Nuclear Information System (INIS)

Zeniya, Tsutomu; Watabe, Hiroshi; Kim, Kyeong Min; Teramoto, Noboru; Hayashi, Takuya; Iida, Hidehiro; Aoi, Toshiyuki; Sohlberg, Antti; Kudo, Hiroyuki

2004-01-01

Pinhole single-photon emission computed tomography (SPECT) is able to provide information on the biodistribution of several radioligands in small laboratory animals, but has limitations associated with non-uniform spatial resolution or axial blurring. We have hypothesised that this blurring is due to incompleteness of the projection data acquired by a single circular pinhole orbit, and have evaluated a new strategy for accurate image reconstruction with better spatial resolution uniformity. A pinhole SPECT system using two circular orbits and a dedicated three-dimensional ordered subsets expectation maximisation (3D-OSEM) reconstruction method were developed. In this system, not the camera but the object rotates, and the two orbits are at 90 and 45 relative to the object's axis. This system satisfies Tuy's condition, and is thus able to provide complete data for 3D pinhole SPECT reconstruction within the whole field of view (FOV). To evaluate this system, a series of experiments was carried out using a multiple-disk phantom filled with 99m Tc solution. The feasibility of the proposed method for small animal imaging was tested with a mouse bone study using 99m Tc-hydroxymethylene diphosphonate. Feldkamp's filtered back-projection (FBP) method and the 3D-OSEM method were applied to these data sets, and the visual and statistical properties were examined. Axial blurring, which was still visible at the edge of the FOV even after applying the conventional 3D-OSEM instead of FBP for single-orbit data, was not visible after application of 3D-OSEM using two-orbit data. 3D-OSEM using two-orbit data dramatically reduced the resolution non-uniformity and statistical noise, and also demonstrated considerably better image quality in the mouse scan. This system may be of use in quantitative assessment of bio-physiological functions in small animals. (orig.)

Image-based point spread function implementation in a fully 3D OSEM reconstruction algorithm for PET.

Science.gov (United States)

Rapisarda, E; Bettinardi, V; Thielemans, K; Gilardi, M C

2010-07-21

The interest in positron emission tomography (PET) and particularly in hybrid integrated PET/CT systems has significantly increased in the last few years due to the improved quality of the obtained images. Nevertheless, one of the most important limits of the PET imaging technique is still its poor spatial resolution due to several physical factors originating both at the emission (e.g. positron range, photon non-collinearity) and at detection levels (e.g. scatter inside the scintillating crystals, finite dimensions of the crystals and depth of interaction). To improve the spatial resolution of the images, a possible way consists of measuring the point spread function (PSF) of the system and then accounting for it inside the reconstruction algorithm. In this work, the system response of the GE Discovery STE operating in 3D mode has been characterized by acquiring (22)Na point sources in different positions of the scanner field of view. An image-based model of the PSF was then obtained by fitting asymmetric two-dimensional Gaussians on the (22)Na images reconstructed with small pixel sizes. The PSF was then incorporated, at the image level, in a three-dimensional ordered subset maximum likelihood expectation maximization (OS-MLEM) reconstruction algorithm. A qualitative and quantitative validation of the algorithm accounting for the PSF has been performed on phantom and clinical data, showing improved spatial resolution, higher contrast and lower noise compared with the corresponding images obtained using the standard OS-MLEM algorithm.

3D-LSI technology for image sensor

International Nuclear Information System (INIS)

Motoyoshi, Makoto; Koyanagi, Mitsumasa

2009-01-01

Recently, the development of three-dimensional large-scale integration (3D-LSI) technologies has accelerated and has advanced from the research level or the limited production level to the investigation level, which might lead to mass production. By separating 3D-LSI technology into elementary technologies such as (1) through silicon via (TSV) formation, (2) bump formation, (3) wafer thinning, (4) chip/wafer alignment, and (5) chip/wafer stacking and reconstructing the entire process and structure, many methods to realize 3D-LSI devices can be developed. However, by considering a specific application, the supply chain of base wafers, and the purpose of 3D integration, a few suitable combinations can be identified. In this paper, we focus on the application of 3D-LSI technologies to image sensors. We describe the process and structure of the chip size package (CSP), developed on the basis of current and advanced 3D-LSI technologies, to be used in CMOS image sensors. Using the current LSI technologies, CSPs for 1.3 M, 2 M, and 5 M pixel CMOS image sensors were successfully fabricated without any performance degradation. 3D-LSI devices can be potentially employed in high-performance focal-plane-array image sensors. We propose a high-speed image sensor with an optical fill factor of 100% to be developed using next-generation 3D-LSI technology and fabricated using micro(μ)-bumps and micro(μ)-TSVs.

Influence of OSEM and segmented attenuation correction in the calculation of standardised uptake values for [18F]FDG PET

International Nuclear Information System (INIS)

Visvikis, D.; Costa, D.C.; Bomanji, J.; Gacinovic, S.; Ell, P.J.; Cheze-LeRest, C.

2001-01-01

Standardised Uptake Values (SUVs) are widely used in positron emission tomography (PET) as a semi-quantitative index of fluorine-18 labelled fluorodeoxyglucose uptake. The objective of this study was to investigate any bias introduced in the calculation of SUVs as a result of employing ordered subsets-expectation maximisation (OSEM) image reconstruction and segmented attenuation correction (SAC). Variable emission and transmission time durations were investigated. Both a phantom and a clinical evaluation of the bias were carried out. The software implemented in the GE Advance PET scanner was used. Phantom studies simulating tumour imaging conditions were performed. Since a variable count rate may influence the results obtained using OSEM, similar acquisitions were performed at total count rates of 34 kcps and 12 kcps. Clinical data consisted of 100 patient studies. Emission datasets of 5 and 15 min duration were combined with 15-, 3-, 2- and 1-min transmission datasets for the reconstruction of both phantom and patient studies. Two SUVs were estimated using the average (SUV avg ) and the maximum (SUV max ) count density from regions of interest placed well inside structures of interest. The percentage bias of these SUVs compared with the values obtained using a reference image was calculated. The reference image was considered to be the one produced by filtered backprojection (FBP) image reconstruction with measured attenuation correction using the 15-min emission and transmission datasets for each phantom and patient study. A bias of 5%-20% was found for the SUV avg and SUV max in the case of FBP with SAC using variable transmission times. In the case of OSEM with SAC, the bias increased to 10%-30%. An overall increase of 5%-10% was observed with the use of SUV max . The 5-min emission dataset led to an increase in the bias of 25%-100%, with the larger increase recorded for the SUV max . The results suggest that OSEM and SAC with 3 and 2 min transmission may be

Reduction of artefacts due to missing projections using OSEM

International Nuclear Information System (INIS)

Hutton, B.F.; Kyme, A.; Choong, K.

2002-01-01

Full text: It is well recognised that missing or corrupted projections can result in artefacts. This occasionally occurs due to errors in data transfer from acquisition memory to disk. A possible approach for reducing these artefacts was investigated, Using ordered subsets expectation maximization (OSEM) the iterative reconstruction proceeds by progressively including additional projections until a single iteration is complete. Clinically useful results can be obtained using a small subset size in a single iteration. Stopping prior to the complete iteration so as to avoid inclusion of missing or corrupted data should provide a 'partial' reconstruction with minimal artefacts. To test this hypothesis projections were selectively removed from a complete data set (2, 4, 8, 12 adjacent projections) and reconstructions were performed using both filtered back projection (FBP) and OSEM. To maintain a constant number of sub-iterations in OSEM an equal number of duplicate projections were substituted for the missing projections. Both 180 and 360 degrees reconstructions with missing data were compared with reconstruction for the complete data using sum of absolute differences. Results indicate that missing data causes artefacts for both FBP and OSEM however the severity of artefacts is significantly reduced using OSEM. The effect of missing data is generally greater for 180 degrees acquisition. OSEM is recommended for minimising reconstruction artefacts due to missing projections. Copyright (2002) The Australian and New Zealand Society of Nuclear Medicine Inc

Study of CT-based positron range correction in high resolution 3D PET imaging

Energy Technology Data Exchange (ETDEWEB)

Cal-Gonzalez, J., E-mail: jacobo@nuclear.fis.ucm.es [Grupo de Fisica Nuclear, Dpto. Fisica Atomica, Molecular y Nuclear, Universidad Complutense de Madrid (Spain); Herraiz, J.L. [Grupo de Fisica Nuclear, Dpto. Fisica Atomica, Molecular y Nuclear, Universidad Complutense de Madrid (Spain); Espana, S. [Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Boston, MA (United States); Vicente, E. [Grupo de Fisica Nuclear, Dpto. Fisica Atomica, Molecular y Nuclear, Universidad Complutense de Madrid (Spain); Instituto de Estructura de la Materia, Consejo Superior de Investigaciones Cientificas (CSIC), Madrid (Spain); Herranz, E. [Grupo de Fisica Nuclear, Dpto. Fisica Atomica, Molecular y Nuclear, Universidad Complutense de Madrid (Spain); Desco, M. [Unidad de Medicina y Cirugia Experimental, Hospital General Universitario Gregorio Maranon, Madrid (Spain); Vaquero, J.J. [Dpto. de Bioingenieria e Ingenieria Espacial, Universidad Carlos III, Madrid (Spain); Udias, J.M. [Grupo de Fisica Nuclear, Dpto. Fisica Atomica, Molecular y Nuclear, Universidad Complutense de Madrid (Spain)

2011-08-21

Positron range limits the spatial resolution of PET images and has a different effect for different isotopes and positron propagation materials. Therefore it is important to consider it during image reconstruction, in order to obtain optimal image quality. Positron range distributions for most common isotopes used in PET in different materials were computed using the Monte Carlo simulations with PeneloPET. The range profiles were introduced into the 3D OSEM image reconstruction software FIRST and employed to blur the image either in the forward projection or in the forward and backward projection. The blurring introduced takes into account the different materials in which the positron propagates. Information on these materials may be obtained, for instance, from a segmentation of a CT image. The results of introducing positron blurring in both forward and backward projection operations was compared to using it only during forward projection. Further, the effect of different shapes of positron range profile in the quality of the reconstructed images with positron range correction was studied. For high positron energy isotopes, the reconstructed images show significant improvement in spatial resolution when positron range is taken into account during reconstruction, compared to reconstructions without positron range modeling.

Study of CT-based positron range correction in high resolution 3D PET imaging

International Nuclear Information System (INIS)

Cal-Gonzalez, J.; Herraiz, J.L.; Espana, S.; Vicente, E.; Herranz, E.; Desco, M.; Vaquero, J.J.; Udias, J.M.

2011-01-01

Positron range limits the spatial resolution of PET images and has a different effect for different isotopes and positron propagation materials. Therefore it is important to consider it during image reconstruction, in order to obtain optimal image quality. Positron range distributions for most common isotopes used in PET in different materials were computed using the Monte Carlo simulations with PeneloPET. The range profiles were introduced into the 3D OSEM image reconstruction software FIRST and employed to blur the image either in the forward projection or in the forward and backward projection. The blurring introduced takes into account the different materials in which the positron propagates. Information on these materials may be obtained, for instance, from a segmentation of a CT image. The results of introducing positron blurring in both forward and backward projection operations was compared to using it only during forward projection. Further, the effect of different shapes of positron range profile in the quality of the reconstructed images with positron range correction was studied. For high positron energy isotopes, the reconstructed images show significant improvement in spatial resolution when positron range is taken into account during reconstruction, compared to reconstructions without positron range modeling.

A 3D HIDAC-PET camera with sub-millimeter resolution for imaging small animals

International Nuclear Information System (INIS)

Jeavons, A.P.; Chandler, R.A.; Dettmar, C.A.R.

1999-01-01

A HIDAC-PET camera consisting essentially of 5 million 0.5 mm gas avalanching detectors has been constructed for small-animal imaging. The particular HIDAC advantage--a high 3D spatial resolution--has been improved to 0.95 mm fwhm and to 0.7 mm fwhm when reconstructing with 3D-OSEM methods incorporating resolution recovery. A depth-of-interaction resolution of 2.5 mm is implicit, due to the laminar construction. Scatter-corrected sensitivity, at 8.9 cps/kBq (i.e. 0.9%) from a central point source, or 7.2 cps/kBq (543 cps/kBq/cm 3 ) from a distributed (40 mm diameter, 60 mm long) source is now much higher than previous, and other, work. A field-of-view of 100 mm (adjustable to 200 mm) diameter by 210 mm axially permits whole-body imaging of small animals, containing typically 4MBqs of activity, at 40 kcps of which 16% are random coincidences, with a typical scatter fraction of 44%. Throughout the field-of-view there are no positional distortions and relative quantitation is uniform to ± 3.5%, but some variation of spatial resolution is found. The performance demonstrates that HIDAC technology is quite appropriate for small-animal PET cameras

3D dynamic pituitary MR imaging with CAIPIRINHA: Initial experience and comparison with 2D dynamic MR imaging

Energy Technology Data Exchange (ETDEWEB)

Fushimi, Yasutaka, E-mail: yfushimi@kuhp.kyoto-u.ac.jp [Department of Diagnostic Imaging and Nuclear Medicine, Kyoto University Graduate School of Medicine, Kyoto 606-8507 (Japan); Okada, Tomohisa; Kanagaki, Mitsunori; Yamamoto, Akira; Kanda, Yumiko; Sakamoto, Ryo [Department of Diagnostic Imaging and Nuclear Medicine, Kyoto University Graduate School of Medicine, Kyoto 606-8507 (Japan); Hojo, Masato; Takahashi, Jun C.; Miyamoto, Susumu [Department of Neurosurgery, Kyoto University Graduate School of Medicine, Kyoto 606-8507 (Japan); Togashi, Kaori [Department of Diagnostic Imaging and Nuclear Medicine, Kyoto University Graduate School of Medicine, Kyoto 606-8507 (Japan)

2014-10-15

Objectives: To evaluate the validity of 3D dynamic pituitary MR imaging with controlled aliasing in parallel imaging results in higher acceleration (CAIPIRINHA), with special emphasis on demarcation of pituitary posterior lobe and stalk. Methods: Participants comprised 32 patients who underwent dynamic pituitary MR imaging due to pituitary or parasellar lesions. 3D dynamic MR with CAIPIRINHA was performed at 3 T with 20-s-interval, precontrast, 1st to 5th dynamic images. Normalized values and enhanced ratios (dynamic postcontrast image values divided by precontrast ones) were compared between 3D and 2D dynamic MR imaging for patients with visual identification of posterior lobe and stalk. Results: In 3D, stalk was identified in 29 patients and unidentified in 3, and posterior lobe was identified in 28 and unidentified in 4. In 2D, stalk was identified in 26 patients and unidentified in 6 patients, and posterior lobe was identified in 15 and unidentified in 17. Normalized values of pituitary posterior lobe and stalk were higher in 3D than 2D (P < 0.001). No significant difference in enhancement ratio was seen between 3D and 2D. Conclusions: 3D dynamic pituitary MR provided better identification and higher normalized values of pituitary posterior lobe and stalk than 2D.

3D tomographic imaging with the γ-eye planar scintigraphic gamma camera

Science.gov (United States)

Tunnicliffe, H.; Georgiou, M.; Loudos, G. K.; Simcox, A.; Tsoumpas, C.

2017-11-01

γ-eye is a desktop planar scintigraphic gamma camera (100 mm × 50 mm field of view) designed by BET Solutions as an affordable tool for dynamic, whole body, small-animal imaging. This investigation tests the viability of using γ-eye for the collection of tomographic data for 3D SPECT reconstruction. Two software packages, QSPECT and STIR (software for tomographic image reconstruction), have been compared. Reconstructions have been performed using QSPECT’s implementation of the OSEM algorithm and STIR’s OSMAPOSL (Ordered Subset Maximum A Posteriori One Step Late) and OSSPS (Ordered Subsets Separable Paraboloidal Surrogate) algorithms. Reconstructed images of phantom and mouse data have been assessed in terms of spatial resolution, sensitivity to varying activity levels and uniformity. The effect of varying the number of iterations, the voxel size (1.25 mm default voxel size reduced to 0.625 mm and 0.3125 mm), the point spread function correction and the weight of prior terms were explored. While QSPECT demonstrated faster reconstructions, STIR outperformed it in terms of resolution (as low as 1 mm versus 3 mm), particularly when smaller voxel sizes were used, and in terms of uniformity, particularly when prior terms were used. Little difference in terms of sensitivity was seen throughout.

Accelerated time-of-flight (TOF) PET image reconstruction using TOF bin subsetization and TOF weighting matrix pre-computation

International Nuclear Information System (INIS)

Mehranian, Abolfazl; Kotasidis, Fotis; Zaidi, Habib

2016-01-01

FDG-PET study also revealed that for the same noise level, a higher contrast recovery can be obtained by increasing the number of TOF subsets. It can be concluded that the proposed TOF weighting matrix pre-computation and subsetization approaches enable to further accelerate and improve the convergence properties of OSEM and MLEM algorithms, thus opening new avenues for accelerated TOF PET image reconstruction. (paper)

3D accelerator magnet calculations using MAGNUS-3D

International Nuclear Information System (INIS)

Pissanetzky, S.; Miao, Y.

1989-01-01

The steady trend towards increased magnetic and geometric complexity in the design of accelerator magnets has caused a need for reliable 3D computer models and a better understanding of the behavior of magnetic system in three dimensions. The capabilities of the MAGNUS-3D family of programs are ideally suited to solve this class of problems and provide insight into 3D effects. MAGNUS-3D can solve any problem of magnetostatics involving permanent magnets, nonlinear ferromagnetic materials and electric conductors. MAGNUS-3D uses the finite element method and the two-scalar-potentials formulation of Maxwell's equations to obtain the solution, which can then be used interactively to obtain tables of field components at specific points or lines, plots of field lines, function graphs representing a field component plotted against a coordinate along any line in space (such as the beam line), and views of the conductors, the mesh and the magnetic bodies. The magnetic quantities that can be calculated include the force or torque on conductors or magnetic parts, the energy, the flux through a specified surface, line integrals of any field component along any line in space, and the average field or potential harmonic coefficients. We describe the programs with emphasis placed on their use for accelerator magnet design, and present an advanced example of actual calculations. (orig.)

GPU acceleration towards real-time image reconstruction in 3D tomographic diffractive microscopy

Science.gov (United States)

Bailleul, J.; Simon, B.; Debailleul, M.; Liu, H.; Haeberlé, O.

2012-06-01

Phase microscopy techniques regained interest in allowing for the observation of unprepared specimens with excellent temporal resolution. Tomographic diffractive microscopy is an extension of holographic microscopy which permits 3D observations with a finer resolution than incoherent light microscopes. Specimens are imaged by a series of 2D holograms: their accumulation progressively fills the range of frequencies of the specimen in Fourier space. A 3D inverse FFT eventually provides a spatial image of the specimen. Consequently, acquisition then reconstruction are mandatory to produce an image that could prelude real-time control of the observed specimen. The MIPS Laboratory has built a tomographic diffractive microscope with an unsurpassed 130nm resolution but a low imaging speed - no less than one minute. Afterwards, a high-end PC reconstructs the 3D image in 20 seconds. We now expect an interactive system providing preview images during the acquisition for monitoring purposes. We first present a prototype implementing this solution on CPU: acquisition and reconstruction are tied in a producer-consumer scheme, sharing common data into CPU memory. Then we present a prototype dispatching some reconstruction tasks to GPU in order to take advantage of SIMDparallelization for FFT and higher bandwidth for filtering operations. The CPU scheme takes 6 seconds for a 3D image update while the GPU scheme can go down to 2 or > 1 seconds depending on the GPU class. This opens opportunities for 4D imaging of living organisms or crystallization processes. We also consider the relevance of GPU for 3D image interaction in our specific conditions.

Quantification of dopaminergic neurotransmission SPECT studies with 123I-labelled radioligands. A comparison between different imaging systems and data acquisition protocols using Monte Carlo simulation

International Nuclear Information System (INIS)

Crespo, Cristina; Aguiar, Pablo; Gallego, Judith; Cot, Albert; Falcon, Carles; Ros, Domenec; Bullich, Santiago; Pareto, Deborah; Sempau, Josep; Lomena, Francisco; Calvino, Francisco; Pavia, Javier

2008-01-01

123 I-labelled radioligands are commonly used for single-photon emission computed tomography (SPECT) imaging of the dopaminergic system to study the dopamine transporter binding. The aim of this work was to compare the quantitative capabilities of two different SPECT systems through Monte Carlo (MC) simulation. The SimSET MC code was employed to generate simulated projections of a numerical phantom for two gamma cameras equipped with a parallel and a fan-beam collimator, respectively. A fully 3D iterative reconstruction algorithm was used to compensate for attenuation, the spatially variant point spread function (PSF) and scatter. A post-reconstruction partial volume effect (PVE) compensation was also developed. For both systems, the correction for all degradations and PVE compensation resulted in recovery factors of the theoretical specific uptake ratio (SUR) close to 100%. For a SUR value of 4, the recovered SUR for the parallel imaging system was 33% for a reconstruction without corrections (OSEM), 45% for a reconstruction with attenuation correction (OSEM-A), 56% for a 3D reconstruction with attenuation and PSF corrections (OSEM-AP), 68% for OSEM-AP with scatter correction (OSEM-APS) and 97% for OSEM-APS plus PVE compensation (OSEM-APSV). For the fan-beam imaging system, the recovered SUR was 41% without corrections, 55% for OSEM-A, 65% for OSEM-AP, 75% for OSEM-APS and 102% for OSEM-APSV. Our findings indicate that the correction for degradations increases the quantification accuracy, with PVE compensation playing a major role in the SUR quantification. The proposed methodology allows us to reach similar SUR values for different SPECT systems, thereby allowing a reliable standardisation in multicentric studies. (orig.)

Quantification of dopaminergic neurotransmission SPECT studies with {sup 123}I-labelled radioligands. A comparison between different imaging systems and data acquisition protocols using Monte Carlo simulation

Energy Technology Data Exchange (ETDEWEB)

Crespo, Cristina; Aguiar, Pablo [Universitat de Barcelona - IDIBAPS, Unitat de Biofisica i Bioenginyeria, Departament de Ciencies Fisiologiques I, Facultat de Medicina, Barcelona (Spain); Gallego, Judith [Universitat Politecnica de Catalunya, Institut de Tecniques Energetiques, Barcelona (Spain); Institut de Bioenginyeria de Catalunya, Barcelona (Spain); Cot, Albert [Universitat de Barcelona - IDIBAPS, Unitat de Biofisica i Bioenginyeria, Departament de Ciencies Fisiologiques I, Facultat de Medicina, Barcelona (Spain); Universitat Politecnica de Catalunya, Seccio d' Enginyeria Nuclear, Departament de Fisica i Enginyeria Nuclear, Barcelona (Spain); Falcon, Carles; Ros, Domenec [Universitat de Barcelona - IDIBAPS, Unitat de Biofisica i Bioenginyeria, Departament de Ciencies Fisiologiques I, Facultat de Medicina, Barcelona (Spain); CIBER en Bioingenieria, Biomateriales y Nanomedicina (CIBER-BBN), Barcelona (Spain); Bullich, Santiago [Hospital del Mar, Center for Imaging in Psychiatry, CRC-MAR, Barcelona (Spain); Pareto, Deborah [CIBER en Bioingenieria, Biomateriales y Nanomedicina (CIBER-BBN), Barcelona (Spain); PRBB, Institut d' Alta Tecnologia, Barcelona (Spain); Sempau, Josep [Universitat Politecnica de Catalunya, Institut de Tecniques Energetiques, Barcelona (Spain); CIBER en Bioingenieria, Biomateriales y Nanomedicina (CIBER-BBN), Barcelona (Spain); Lomena, Francisco [IDIBAPS, Servei de Medicina Nuclear, Hospital Clinic, Barcelona (Spain); Calvino, Francisco [Universitat Politecnica de Catalunya, Institut de Tecniques Energetiques, Barcelona (Spain); Universitat Politecnica de Catalunya, Seccio d' Enginyeria Nuclear, Departament de Fisica i Enginyeria Nuclear, Barcelona (Spain); Pavia, Javier [CIBER en Bioingenieria, Biomateriales y Nanomedicina (CIBER-BBN), Barcelona (Spain); IDIBAPS, Servei de Medicina Nuclear, Hospital Clinic, Barcelona (Spain)

2008-07-15

{sup 123}I-labelled radioligands are commonly used for single-photon emission computed tomography (SPECT) imaging of the dopaminergic system to study the dopamine transporter binding. The aim of this work was to compare the quantitative capabilities of two different SPECT systems through Monte Carlo (MC) simulation. The SimSET MC code was employed to generate simulated projections of a numerical phantom for two gamma cameras equipped with a parallel and a fan-beam collimator, respectively. A fully 3D iterative reconstruction algorithm was used to compensate for attenuation, the spatially variant point spread function (PSF) and scatter. A post-reconstruction partial volume effect (PVE) compensation was also developed. For both systems, the correction for all degradations and PVE compensation resulted in recovery factors of the theoretical specific uptake ratio (SUR) close to 100%. For a SUR value of 4, the recovered SUR for the parallel imaging system was 33% for a reconstruction without corrections (OSEM), 45% for a reconstruction with attenuation correction (OSEM-A), 56% for a 3D reconstruction with attenuation and PSF corrections (OSEM-AP), 68% for OSEM-AP with scatter correction (OSEM-APS) and 97% for OSEM-APS plus PVE compensation (OSEM-APSV). For the fan-beam imaging system, the recovered SUR was 41% without corrections, 55% for OSEM-A, 65% for OSEM-AP, 75% for OSEM-APS and 102% for OSEM-APSV. Our findings indicate that the correction for degradations increases the quantification accuracy, with PVE compensation playing a major role in the SUR quantification. The proposed methodology allows us to reach similar SUR values for different SPECT systems, thereby allowing a reliable standardisation in multicentric studies. (orig.)

TU-F-12A-05: Sensitivity of Textural Features to 3D Vs. 4D FDG-PET/CT Imaging in NSCLC Patients

Energy Technology Data Exchange (ETDEWEB)

Yang, F; Nyflot, M; Bowen, S; Kinahan, P; Sandison, G [University of Washington Medical Center, Seattle, WA (United States)

2014-06-15

Purpose: Neighborhood Gray-level difference matrices (NGLDM) based texture parameters extracted from conventional (3D) 18F-FDG PET scans in patients with NSCLC have been previously shown to associate with response to chemoradiation and poorer patient outcome. However, the change in these parameters when utilizing respiratory-correlated (4D) FDG-PET scans has not yet been characterized for NSCLC. The Objectives: of this study was to assess the extent to which NGLDM-based texture parameters on 4D PET images vary with reference to values derived from 3D scans in NSCLC. Methods: Eight patients with newly diagnosed NSCLC treated with concomitant chemoradiotherapy were included in this study. 4D PET scans were reconstructed with OSEM-IR in 5 respiratory phase-binned images and corresponding CT data of each phase were employed for attenuation correction. NGLDM-based texture features, consisting of coarseness, contrast, busyness, complexity and strength, were evaluated for gross tumor volumes defined on 3D/4D PET scans by radiation oncologists. Variation of the obtained texture parameters over the respiratory cycle were examined with respect to values extracted from 3D scans. Results: Differences between texture parameters derived from 4D scans at different respiratory phases and those extracted from 3D scans ranged from −30% to 13% for coarseness, −12% to 40% for contrast, −5% to 50% for busyness, −7% to 38% for complexity, and −43% to 20% for strength. Furthermore, no evident correlations were observed between respiratory phase and 4D scan texture parameters. Conclusion: Results of the current study showed that NGLDM-based texture parameters varied considerably based on choice of 3D PET and 4D PET reconstruction of NSCLC patient images, indicating that standardized image acquisition and analysis protocols need to be established for clinical studies, especially multicenter clinical trials, intending to validate prognostic values of texture features for NSCLC.

Fully 3-D list-mode positron emission tomography image reconstruction on a multi-GPU cluster

Energy Technology Data Exchange (ETDEWEB)

Cui, Jingyu [Stanford Univ., CA (United States). Dept. of Electrical Engineering; Prevrhal, Sven; Shao, Lingxiong [Philips Healthcare, San Jose, CA (United States); Pratx, Guillem [Stanford Univ., CA (United States). Dept. of Radiation Oncology; Levin, Craig S. [Stanford Univ., CA (United States). Dept. of Radiology, Electrical Engineering, and Physics; Stanford Univ., CA (United States). Molecular Imaging Program at Stanford (MIPS); Stanford Univ., CA (United States). School of Medicine

2011-07-01

List-mode processing is an efficient way of dealing with the sparse nature of PET data sets, and is the processing method of choice for time-of-flight (ToF) PET. We present a novel method of computing line projection operations required for list-mode ordered subsets expectation maximization (OSEM) for fully 3-D PET image reconstruction on a graphics processing unit (GPU) using the compute unified device architecture (CUDA) framework. Our method overcomes challenges such as compute thread divergence, and exploits GPU capabilities such as shared memory and atomic operations. When applied to line projection operations for list-mode time-of-flight PET, this new GPU-CUDA reformulation is 188X faster than a single-threaded reference CPU implementation. When embedded in a multi-process environment on a GPU-equipped small cluster, a speedup of 4X was observed over the same configuration but without GPU support. Image quality is preserved with root mean squared (RMS) deviation of 0.05% between CPU and GPU-generated images, which has negligible effect in typical clinical applications. (orig.)

Evaluation of bias and variance in low-count OSEM list mode reconstruction

International Nuclear Information System (INIS)

Jian, Y; Carson, R E; Planeta, B

2015-01-01

Statistical algorithms have been widely used in PET image reconstruction. The maximum likelihood expectation maximization reconstruction has been shown to produce bias in applications where images are reconstructed from a relatively small number of counts. In this study, image bias and variability in low-count OSEM reconstruction are investigated on images reconstructed with MOLAR (motion-compensation OSEM list-mode algorithm for resolution-recovery reconstruction) platform. A human brain ([ 11 C]AFM) and a NEMA phantom are used in the simulation and real experiments respectively, for the HRRT and Biograph mCT. Image reconstructions were repeated with different combinations of subsets and iterations. Regions of interest were defined on low-activity and high-activity regions to evaluate the bias and noise at matched effective iteration numbers (iterations × subsets). Minimal negative biases and no positive biases were found at moderate count levels and less than 5% negative bias was found using extremely low levels of counts (0.2 M NEC). At any given count level, other factors, such as subset numbers and frame-based scatter correction may introduce small biases (1–5%) in the reconstructed images. The observed bias was substantially lower than that reported in the literature, perhaps due to the use of point spread function and/or other implementation methods in MOLAR. (paper)

Fast MR Imaging of the Paediatric Abdomen with CAIPIRINHA-Accelerated T1w 3D FLASH and with High-Resolution T2w HASTE: A Study on Image Quality

Directory of Open Access Journals (Sweden)

Mengxia Li

2015-01-01

Full Text Available The aim of this study was to explore the applicability of fast MR techniques to routine paediatric abdominopelvic MRI at 1.5 Tesla. “Controlled Aliasing in Parallel Imaging Results in Higher Acceleration-” (CAIPIRINHA- accelerated contrast-enhanced-T1w 3D FLASH imaging was compared to standard T1w 2D FLASH imaging with breath-holding in 40 paediatric patients and to respiratory-triggered T1w TSE imaging in 10 sedated young children. In 20 nonsedated patients, we compared T2w TIRM to fat-saturated T2w HASTE imaging. Two observers performed an independent and blinded assessment of overall image quality. Acquisition time was reduced by the factor of 15 with CAIPIRINHA-accelerated T1w FLASH and by 7 with T2w HASTE. With CAIPIRINHA and with HASTE, there were significantly less motion artefacts in nonsedated patients. In sedated patients, respiratory-triggered T1w imaging in general showed better image quality. However, satisfactory image quality was achieved with CAIPIRINHA in two sedated patients where respiratory triggering failed. In summary, fast scanning with CAIPIRINHA and HASTE presents a reliable high quality alternative to standard sequences in paediatric abdominal MRI. Paediatric patients, in particular, benefit greatly from fast image acquisition with less breath-hold cycles or shorter sedation.

GPU-accelerated Kernel Regression Reconstruction for Freehand 3D Ultrasound Imaging.

Science.gov (United States)

Wen, Tiexiang; Li, Ling; Zhu, Qingsong; Qin, Wenjian; Gu, Jia; Yang, Feng; Xie, Yaoqin

2017-07-01

Volume reconstruction method plays an important role in improving reconstructed volumetric image quality for freehand three-dimensional (3D) ultrasound imaging. By utilizing the capability of programmable graphics processing unit (GPU), we can achieve a real-time incremental volume reconstruction at a speed of 25-50 frames per second (fps). After incremental reconstruction and visualization, hole-filling is performed on GPU to fill remaining empty voxels. However, traditional pixel nearest neighbor-based hole-filling fails to reconstruct volume with high image quality. On the contrary, the kernel regression provides an accurate volume reconstruction method for 3D ultrasound imaging but with the cost of heavy computational complexity. In this paper, a GPU-based fast kernel regression method is proposed for high-quality volume after the incremental reconstruction of freehand ultrasound. The experimental results show that improved image quality for speckle reduction and details preservation can be obtained with the parameter setting of kernel window size of [Formula: see text] and kernel bandwidth of 1.0. The computational performance of the proposed GPU-based method can be over 200 times faster than that on central processing unit (CPU), and the volume with size of 50 million voxels in our experiment can be reconstructed within 10 seconds.

Integrated optical 3D digital imaging based on DSP scheme

Science.gov (United States)

Wang, Xiaodong; Peng, Xiang; Gao, Bruce Z.

2008-03-01

We present a scheme of integrated optical 3-D digital imaging (IO3DI) based on digital signal processor (DSP), which can acquire range images independently without PC support. This scheme is based on a parallel hardware structure with aid of DSP and field programmable gate array (FPGA) to realize 3-D imaging. In this integrated scheme of 3-D imaging, the phase measurement profilometry is adopted. To realize the pipeline processing of the fringe projection, image acquisition and fringe pattern analysis, we present a multi-threads application program that is developed under the environment of DSP/BIOS RTOS (real-time operating system). Since RTOS provides a preemptive kernel and powerful configuration tool, with which we are able to achieve a real-time scheduling and synchronization. To accelerate automatic fringe analysis and phase unwrapping, we make use of the technique of software optimization. The proposed scheme can reach a performance of 39.5 f/s (frames per second), so it may well fit into real-time fringe-pattern analysis and can implement fast 3-D imaging. Experiment results are also presented to show the validity of proposed scheme.

MAGNUS-3D: Accelerator magnet calculations in 3-dimensions

Science.gov (United States)

Pissanetzky, S.

1988-12-01

MAGNUS-3D is a professional finite element code for nonlinear magnetic engineering. MAGNUS-3D can solve numerically any general problem of linear or nonlinear magnetostatics in three dimensions. The problem is formulated in a domain with Dirichlet, Neumann or periodic boundary conditions, that can contain any combination of conductors of any shape in space, nonlinear magnetic materials with magnetic properties specified by magnetization tables, and nonlinear permanent magnets with any given demagnetization curve. MAGNUS-3D uses the two-scalar-potentials formulation of Magnetostatics and the finite element method, has an automatic 3D mesh generator, and advanced post-processing features that include graphics on a variety of supported devices, tabulation, and calculation of design quantities required in Magnetic Engineering. MAGNUS-3D is a general purpose 3D code, but it has been extensively used for accelerator work and many special features required for accelerator engineering have been incorporated into the code. One of such features is the calculation of field harmonic coefficients averaged in the direction of the beam, so important for the design of magnet ends. Another feature is its ability to calculate line integrals of any field component along the direction of the beam, or plot the field as a function of the z coordinate. MAGNUS-3D has found applications to the design of accelerator magnets and spectrometers, steering magnets, wigglers and undulators for free electron lasers, microtrons and magnets for synchrotron light sources, as well as magnets for NMR and medical applications, recording heads and various magnetic devices. There are three more programs closely associated with MAGNUS-3D. MAGNUS-GKS is the graphical postprocessor for the package; it supports a numer of output devices, including color vector or bit map devices. WIRE is an independent program that can calculate the field produced by any configuration of electric conductors in space, at any

3D Metallic Lattices for Accelerator Applications

CERN Document Server

Shapiro, Michael A; Sirigiri, Jagadishwar R; Temkin, Richard J

2005-01-01

We present the results of research on 3D metallic lattices operating at microwave frequencies for application in (1) accelerator structures with higher order mode suppression, (2) Smith-Purcell radiation beam diagnostics, and (3) polaritonic materials for laser acceleration. Electromagnetic waves in a 3D simple cubic lattice formed by metal wires are calculated using HFSS. The bulk modes in the lattice are determined using single cell calculations with different phase advances in all three directions. The Brillouin diagram for the bulk modes is presented and indicates the absence of band gaps in simple lattices except the band below the cutoff. Lattices with thin wires as well as with thick wires have been analyzed. The Brillouin diagram also indicates the presence of low frequency 3D plasmon mode as well as the two degenerate photon modes analogous to those in a 2D lattice. Surface modes for a semi-infinite cubic lattice are modeled as a stack of cells with different phase advances in the two directions alon...

Advanced 3-D Ultrasound Imaging

DEFF Research Database (Denmark)

Rasmussen, Morten Fischer

The main purpose of the PhD project was to develop methods that increase the 3-D ultrasound imaging quality available for the medical personnel in the clinic. Acquiring a 3-D volume gives the medical doctor the freedom to investigate the measured anatomy in any slice desirable after the scan has...... been completed. This allows for precise measurements of organs dimensions and makes the scan more operator independent. Real-time 3-D ultrasound imaging is still not as widespread in use in the clinics as 2-D imaging. A limiting factor has traditionally been the low image quality achievable using...... a channel limited 2-D transducer array and the conventional 3-D beamforming technique, Parallel Beamforming. The first part of the scientific contributions demonstrate that 3-D synthetic aperture imaging achieves a better image quality than the Parallel Beamforming technique. Data were obtained using both...

3D-2D Deformable Image Registration Using Feature-Based Nonuniform Meshes.

Science.gov (United States)

Zhong, Zichun; Guo, Xiaohu; Cai, Yiqi; Yang, Yin; Wang, Jing; Jia, Xun; Mao, Weihua

2016-01-01

By using prior information of planning CT images and feature-based nonuniform meshes, this paper demonstrates that volumetric images can be efficiently registered with a very small portion of 2D projection images of a Cone-Beam Computed Tomography (CBCT) scan. After a density field is computed based on the extracted feature edges from planning CT images, nonuniform tetrahedral meshes will be automatically generated to better characterize the image features according to the density field; that is, finer meshes are generated for features. The displacement vector fields (DVFs) are specified at the mesh vertices to drive the deformation of original CT images. Digitally reconstructed radiographs (DRRs) of the deformed anatomy are generated and compared with corresponding 2D projections. DVFs are optimized to minimize the objective function including differences between DRRs and projections and the regularity. To further accelerate the above 3D-2D registration, a procedure to obtain good initial deformations by deforming the volume surface to match 2D body boundary on projections has been developed. This complete method is evaluated quantitatively by using several digital phantoms and data from head and neck cancer patients. The feature-based nonuniform meshing method leads to better results than either uniform orthogonal grid or uniform tetrahedral meshes.

3D-2D Deformable Image Registration Using Feature-Based Nonuniform Meshes

Directory of Open Access Journals (Sweden)

Zichun Zhong

2016-01-01

Full Text Available By using prior information of planning CT images and feature-based nonuniform meshes, this paper demonstrates that volumetric images can be efficiently registered with a very small portion of 2D projection images of a Cone-Beam Computed Tomography (CBCT scan. After a density field is computed based on the extracted feature edges from planning CT images, nonuniform tetrahedral meshes will be automatically generated to better characterize the image features according to the density field; that is, finer meshes are generated for features. The displacement vector fields (DVFs are specified at the mesh vertices to drive the deformation of original CT images. Digitally reconstructed radiographs (DRRs of the deformed anatomy are generated and compared with corresponding 2D projections. DVFs are optimized to minimize the objective function including differences between DRRs and projections and the regularity. To further accelerate the above 3D-2D registration, a procedure to obtain good initial deformations by deforming the volume surface to match 2D body boundary on projections has been developed. This complete method is evaluated quantitatively by using several digital phantoms and data from head and neck cancer patients. The feature-based nonuniform meshing method leads to better results than either uniform orthogonal grid or uniform tetrahedral meshes.

3D composite image, 3D MRI, 3D SPECT, hydrocephalus

International Nuclear Information System (INIS)

Mito, T.; Shibata, I.; Sugo, N.; Takano, M.; Takahashi, H.

2002-01-01

The three-dimensional (3D)SPECT imaging technique we have studied and published for the past several years is an analytical tool that permits visual expression of the cerebral circulation profile in various cerebral diseases. The greatest drawback of SPECT is that the limitation on precision of spacial resolution makes intracranial localization impossible. In 3D SPECT imaging, intracranial volume and morphology may vary with the threshold established. To solve this problem, we have produced complimentarily combined SPECT and helical-CT 3D images by means of general-purpose visualization software for intracranial localization. In hydrocephalus, however, the key subject to be studied is the profile of cerebral circulation around the ventricles of the brain. This suggests that, for displaying the cerebral ventricles in three dimensions, CT is a difficult technique whereas MRI is more useful. For this reason, we attempted to establish the profile of cerebral circulation around the cerebral ventricles by the production of combined 3D images of SPECT and MRI. In patients who had shunt surgery for hydrocephalus, a difference between pre- and postoperative cerebral circulation profiles was assessed by a voxel distribution curve, 3D SPECT images, and combined 3D SPECT and MRI images. As the shunt system in this study, an Orbis-Sigma valve of the automatic cerebrospinal fluid volume adjustment type was used in place of the variable pressure type Medos valve currently in use, because this device requires frequent changes in pressure and a change in pressure may be detected after MRI procedure. The SPECT apparatus used was PRISM3000 of the three-detector type, and 123I-IMP was used as the radionuclide in a dose of 222 MBq. MRI data were collected with an MAGNEXa+2 with a magnetic flux density of 0.5 tesla under the following conditions: field echo; TR 50 msec; TE, 10 msec; flip, 30ueK; 1 NEX; FOV, 23 cm; 1-mm slices; and gapless. 3D images are produced on the workstation TITAN

GPU-Based Block-Wise Nonlocal Means Denoising for 3D Ultrasound Images

Directory of Open Access Journals (Sweden)

Liu Li

2013-01-01

Full Text Available Speckle suppression plays an important role in improving ultrasound (US image quality. While lots of algorithms have been proposed for 2D US image denoising with remarkable filtering quality, there is relatively less work done on 3D ultrasound speckle suppression, where the whole volume data rather than just one frame needs to be considered. Then, the most crucial problem with 3D US denoising is that the computational complexity increases tremendously. The nonlocal means (NLM provides an effective method for speckle suppression in US images. In this paper, a programmable graphic-processor-unit- (GPU- based fast NLM filter is proposed for 3D ultrasound speckle reduction. A Gamma distribution noise model, which is able to reliably capture image statistics for Log-compressed ultrasound images, was used for the 3D block-wise NLM filter on basis of Bayesian framework. The most significant aspect of our method was the adopting of powerful data-parallel computing capability of GPU to improve the overall efficiency. Experimental results demonstrate that the proposed method can enormously accelerate the algorithm.

Enhanced FIB-SEM systems for large-volume 3D imaging

Science.gov (United States)

Xu, C Shan; Hayworth, Kenneth J; Lu, Zhiyuan; Grob, Patricia; Hassan, Ahmed M; García-Cerdán, José G; Niyogi, Krishna K; Nogales, Eva; Weinberg, Richard J; Hess, Harald F

2017-01-01

Focused Ion Beam Scanning Electron Microscopy (FIB-SEM) can automatically generate 3D images with superior z-axis resolution, yielding data that needs minimal image registration and related post-processing. Obstacles blocking wider adoption of FIB-SEM include slow imaging speed and lack of long-term system stability, which caps the maximum possible acquisition volume. Here, we present techniques that accelerate image acquisition while greatly improving FIB-SEM reliability, allowing the system to operate for months and generating continuously imaged volumes > 106 µm3. These volumes are large enough for connectomics, where the excellent z resolution can help in tracing of small neuronal processes and accelerate the tedious and time-consuming human proofreading effort. Even higher resolution can be achieved on smaller volumes. We present example data sets from mammalian neural tissue, Drosophila brain, and Chlamydomonas reinhardtii to illustrate the power of this novel high-resolution technique to address questions in both connectomics and cell biology. DOI: http://dx.doi.org/10.7554/eLife.25916.001 PMID:28500755

Contributions in compression of 3D medical images and 2D images; Contributions en compression d'images medicales 3D et d'images naturelles 2D

Energy Technology Data Exchange (ETDEWEB)

Gaudeau, Y

2006-12-15

The huge amounts of volumetric data generated by current medical imaging techniques in the context of an increasing demand for long term archiving solutions, as well as the rapid development of distant radiology make the use of compression inevitable. Indeed, if the medical community has sided until now with compression without losses, most of applications suffer from compression ratios which are too low with this kind of compression. In this context, compression with acceptable losses could be the most appropriate answer. So, we propose a new loss coding scheme based on 3D (3 dimensional) Wavelet Transform and Dead Zone Lattice Vector Quantization 3D (DZLVQ) for medical images. Our algorithm has been evaluated on several computerized tomography (CT) and magnetic resonance image volumes. The main contribution of this work is the design of a multidimensional dead zone which enables to take into account correlations between neighbouring elementary volumes. At high compression ratios, we show that it can out-perform visually and numerically the best existing methods. These promising results are confirmed on head CT by two medical patricians. The second contribution of this document assesses the effect with-loss image compression on CAD (Computer-Aided Decision) detection performance of solid lung nodules. This work on 120 significant lungs images shows that detection did not suffer until 48:1 compression and still was robust at 96:1. The last contribution consists in the complexity reduction of our compression scheme. The first allocation dedicated to 2D DZLVQ uses an exponential of the rate-distortion (R-D) functions. The second allocation for 2D and 3D medical images is based on block statistical model to estimate the R-D curves. These R-D models are based on the joint distribution of wavelet vectors using a multidimensional mixture of generalized Gaussian (MMGG) densities. (author)

MRCP and 3D LAVA imaging of extrahepatic cholangiocarcinoma at 3 T MRI

Energy Technology Data Exchange (ETDEWEB)

Li, N.; Liu, C.; Bi, W.; Lin, X.; Jiao, H. [Shandong Medical Imaging Research Institute, Shandong University, Jinan (China); Zhao, P., E-mail: Gavinsdu@163.com [Department of Radiology, Shandong Provincial Hospital, Shandong University, Jinan (China)

2012-06-15

Extrahepatic cholangiocarcinoma (CCA) is a primary bile duct malignant tumour with poor prognosis. Familiarity with their varied imaging characteristics can be helpful in developing a correct diagnosis and in optimal treatment planning, and thus contribute to a better prognosis. The purpose of this article is to illustrate the typical appearances of extrahepatic CCA on magnetic resonance cholangiopancreatography (MRCP) and three-dimensional (3D) LAVA (liver acquisition with volume acceleration) sequences at 3 T magnetic resonance imaging (MRI), and to discuss the superiority of the two techniques in the diagnosis of CCA.

Accelerated median root prior reconstruction for pinhole single-photon emission tomography (SPET)

Energy Technology Data Exchange (ETDEWEB)

Sohlberg, Antti [Department of Clinical Physiology and Nuclear Medicine, Kuopio University Hospital, PO Box 1777 FIN-70211, Kuopio (Finland); Ruotsalainen, Ulla [Institute of Signal Processing, DMI, Tampere University of Technology, PO Box 553 FIN-33101, Tampere (Finland); Watabe, Hiroshi [National Cardiovascular Center Research Institute, 5-7-1 Fujisihro-dai, Suita City, Osaka 565-8565 (Japan); Iida, Hidehiro [National Cardiovascular Center Research Institute, 5-7-1 Fujisihro-dai, Suita City, Osaka 565-8565 (Japan); Kuikka, Jyrki T [Department of Clinical Physiology and Nuclear Medicine, Kuopio University Hospital, PO Box 1777 FIN-70211, Kuopio (Finland)

2003-07-07

Pinhole collimation can be used to improve spatial resolution in SPET. However, the resolution improvement is achieved at the cost of reduced sensitivity, which leads to projection images with poor statistics. Images reconstructed from these projections using the maximum likelihood expectation maximization (ML-EM) algorithms, which have been used to reduce the artefacts generated by the filtered backprojection (FBP) based reconstruction, suffer from noise/bias trade-off: noise contaminates the images at high iteration numbers, whereas early abortion of the algorithm produces images that are excessively smooth and biased towards the initial estimate of the algorithm. To limit the noise accumulation we propose the use of the pinhole median root prior (PH-MRP) reconstruction algorithm. MRP is a Bayesian reconstruction method that has already been used in PET imaging and shown to possess good noise reduction and edge preservation properties. In this study the PH-MRP algorithm was accelerated with the ordered subsets (OS) procedure and compared to the FBP, OS-EM and conventional Bayesian reconstruction methods in terms of noise reduction, quantitative accuracy, edge preservation and visual quality. The results showed that the accelerated PH-MRP algorithm was very robust. It provided visually pleasing images with lower noise level than the FBP or OS-EM and with smaller bias and sharper edges than the conventional Bayesian methods.

Optimization of PET image quality by means of 3D data acquisition and iterative image reconstruction

International Nuclear Information System (INIS)

Doll, J.; Zaers, J.; Trojan, H.; Bellemann, M.E.; Adam, L.E.; Haberkorn, U.; Brix, G.

1998-01-01

The experiments were performed at the latest-generation whole-body PET system ECAT EXACT HR + . For 2D data acquisition, a collimator of thin tungsten septa was positioned in the field-of-view. Prior to image reconstruction, the measured 3D data were sorted into 2D sinograms by using the Fourier rebinning (FORE) algorithm developed by M. Defrise. The standard filtered backprojection (FBP) method and an optimized ML/EM algorithm with overrelaxation for accelerated convergence were employed for image reconstruction. The spatial resolution of both methods as well as the convergence and noise properties of the ML/EM algorithm were studied in phantom measurements. Furthermore, patient data were acquired in the 2D mode as well as in the 3D mode and reconstructed with both techniques. At the same spatial resolution, the ML/EM-reconstructed images showed fewer and less prominent artefacts than the FBP-reconstructed images. The resulting improved detail conspicuously was achieved for the data acquired in the 2D mode as well as in the 3D mode. The best image quality was obtained by iterative 2D reconstruction of 3D data sets which were previously rebinned into 2D sinograms with help of the FORE algorithm. The phantom measurements revealed that 50 iteration steps with the otpimized ML/EM algorithm were sufficient to keep the relative quantitation error below 5%. (orig./MG) [de

3D Hyperpolarized C-13 EPI with Calibrationless Parallel Imaging

DEFF Research Database (Denmark)

Gordon, Jeremy W.; Hansen, Rie Beck; Shin, Peter J.

2018-01-01

With the translation of metabolic MRI with hyperpolarized 13C agents into the clinic, imaging approaches will require large volumetric FOVs to support clinical applications. Parallel imaging techniques will be crucial to increasing volumetric scan coverage while minimizing RF requirements and tem...... strategies to accelerate and undersample hyperpolarized 13C data using 3D blipped EPI acquisitions and multichannel receive coils, and demonstrated its application in a human study of [1-13C]pyruvate metabolism....

Three-dimensional cardiac cine imaging using the kat ARC acceleration: Initial experience in clinical adult patients at 3T.

Science.gov (United States)

Okuda, Shigeo; Yamada, Yoshitake; Tanimoto, Akihiro; Fujita, Jun; Sano, Motoaki; Fukuda, Keiichi; Kuribayashi, Sachio; Jinzaki, Masahiro; Nozaki, Atsushi; Lai, Peng

2015-09-01

Three-dimensional cardiac cine imaging has demonstrated promising clinical 1.5-Tesla results; however, its application to 3T scanners has been limited because of the higher sensitivity to off-resonance artifacts. The aim of this study was to apply 3D cardiac cine imaging during a single breath hold in clinical patients on a 3T scanner using the kat ARC (k- and adaptive-t auto-calibrating reconstruction for Cartesian sampling) technique and to evaluate the interchangeability between 2D and 3D cine imaging for cardiac functional analysis and detection of abnormalities in regional wall motion. Following institutional review board approval, we obtained 2D cine images with an acceleration factor of two during multiple breath holds and 3D cine images with a net scan acceleration factor of 7.7 during a single breath hold in 20 patients using a 3T unit. Two readers independently evaluated the wall motion of the left ventricle (LV) using a 5-point scale, and the consistency in the detection of regional wall motion abnormality between 2D and 3D cine was analyzed by Cohen's kappa test. The LV volume was calculated at end-diastole and end-systole (LVEDV, LVESV); the ejection fraction (LVEF) and myocardial weight (LVmass) were also calculated. The relationship between functional parameters calculated for 2D and 3D cine images was analyzed using Pearson's correlation analysis. The bias and 95% limit of agreement (LA) were calculated using Bland-Altman plots. In addition, a qualitative evaluation of image quality was performed with regard to the myocardium-blood contrast, noise level and boundary definition. Despite slight degradation in image quality for 3D cine, excellent agreement was obtained for the detection of wall motion abnormalities between 2D and 3D cine images (κ=0.84 and 0.94 for each reader). Excellent correlations between the two imaging methods were shown for the evaluation of functional parameters (r>0.97). Slight differences in LVEDV, LVESV, LVEF and LVmass

3D vector flow imaging

DEFF Research Database (Denmark)

Pihl, Michael Johannes

The main purpose of this PhD project is to develop an ultrasonic method for 3D vector flow imaging. The motivation is to advance the field of velocity estimation in ultrasound, which plays an important role in the clinic. The velocity of blood has components in all three spatial dimensions, yet...... are (vx, vy, vz) = (-0.03, 95, 1.0) ± (9, 6, 1) cm/s compared with the expected (0, 96, 0) cm/s. Afterwards, 3D vector flow images from a cross-sectional plane of the vessel are presented. The out of plane velocities exhibit the expected 2D circular-symmetric parabolic shape. The experimental results...... verify that the 3D TO method estimates the complete 3D velocity vectors, and that the method is suitable for 3D vector flow imaging....

Contributions in compression of 3D medical images and 2D images; Contributions en compression d'images medicales 3D et d'images naturelles 2D

Energy Technology Data Exchange (ETDEWEB)

Gaudeau, Y

2006-12-15

The huge amounts of volumetric data generated by current medical imaging techniques in the context of an increasing demand for long term archiving solutions, as well as the rapid development of distant radiology make the use of compression inevitable. Indeed, if the medical community has sided until now with compression without losses, most of applications suffer from compression ratios which are too low with this kind of compression. In this context, compression with acceptable losses could be the most appropriate answer. So, we propose a new loss coding scheme based on 3D (3 dimensional) Wavelet Transform and Dead Zone Lattice Vector Quantization 3D (DZLVQ) for medical images. Our algorithm has been evaluated on several computerized tomography (CT) and magnetic resonance image volumes. The main contribution of this work is the design of a multidimensional dead zone which enables to take into account correlations between neighbouring elementary volumes. At high compression ratios, we show that it can out-perform visually and numerically the best existing methods. These promising results are confirmed on head CT by two medical patricians. The second contribution of this document assesses the effect with-loss image compression on CAD (Computer-Aided Decision) detection performance of solid lung nodules. This work on 120 significant lungs images shows that detection did not suffer until 48:1 compression and still was robust at 96:1. The last contribution consists in the complexity reduction of our compression scheme. The first allocation dedicated to 2D DZLVQ uses an exponential of the rate-distortion (R-D) functions. The second allocation for 2D and 3D medical images is based on block statistical model to estimate the R-D curves. These R-D models are based on the joint distribution of wavelet vectors using a multidimensional mixture of generalized Gaussian (MMGG) densities. (author)

Highly-accelerated self-gated free-breathing 3D cardiac cine MRI: validation in assessment of left ventricular function.

Science.gov (United States)

Liu, Jing; Feng, Li; Shen, Hsin-Wei; Zhu, Chengcheng; Wang, Yan; Mukai, Kanae; Brooks, Gabriel C; Ordovas, Karen; Saloner, David

2017-08-01

This work presents a highly-accelerated, self-gated, free-breathing 3D cardiac cine MRI method for cardiac function assessment. A golden-ratio profile based variable-density, pseudo-random, Cartesian undersampling scheme was implemented for continuous 3D data acquisition. Respiratory self-gating was achieved by deriving motion signal from the acquired MRI data. A multi-coil compressed sensing technique was employed to reconstruct 4D images (3D+time). 3D cardiac cine imaging with self-gating was compared to bellows gating and the clinical standard breath-held 2D cine imaging for evaluation of self-gating accuracy, image quality, and cardiac function in eight volunteers. Reproducibility of 3D imaging was assessed. Self-gated 3D imaging provided an image quality score of 3.4 ± 0.7 vs 4.0 ± 0 with the 2D method (p = 0.06). It determined left ventricular end-systolic volume as 42.4 ± 11.5 mL, end-diastolic volume as 111.1 ± 24.7 mL, and ejection fraction as 62.0 ± 3.1%, which were comparable to the 2D method, with bias ± 1.96 × SD of -0.8 ± 7.5 mL (p = 0.90), 2.6 ± 3.3 mL (p = 0.84) and 1.4 ± 6.4% (p = 0.45), respectively. The proposed 3D cardiac cine imaging method enables reliable respiratory self-gating performance with good reproducibility, and provides comparable image quality and functional measurements to 2D imaging, suggesting that self-gated, free-breathing 3D cardiac cine MRI framework is promising for improved patient comfort and cardiac MRI scan efficiency.

Clinical performance of a free-breathing spatiotemporally accelerated 3-D time-resolved contrast-enhanced pediatric abdominal MR angiography.

Science.gov (United States)

Zhang, Tao; Yousaf, Ufra; Hsiao, Albert; Cheng, Joseph Y; Alley, Marcus T; Lustig, Michael; Pauly, John M; Vasanawala, Shreyas S

2015-10-01

Pediatric contrast-enhanced MR angiography is often limited by respiration, other patient motion and compromised spatiotemporal resolution. To determine the reliability of a free-breathing spatiotemporally accelerated 3-D time-resolved contrast-enhanced MR angiography method for depicting abdominal arterial anatomy in young children. With IRB approval and informed consent, we retrospectively identified 27 consecutive children (16 males and 11 females; mean age: 3.8 years, range: 14 days to 8.4 years) referred for contrast-enhanced MR angiography at our institution, who had undergone free-breathing spatiotemporally accelerated time-resolved contrast-enhanced MR angiography studies. A radio-frequency-spoiled gradient echo sequence with Cartesian variable density k-space sampling and radial view ordering, intrinsic motion navigation and intermittent fat suppression was developed. Images were reconstructed with soft-gated parallel imaging locally low-rank method to achieve both motion correction and high spatiotemporal resolution. Quality of delineation of 13 abdominal arteries in the reconstructed images was assessed independently by two radiologists on a five-point scale. Ninety-five percent confidence intervals of the proportion of diagnostically adequate cases were calculated. Interobserver agreements were also analyzed. Eleven out of 13 arteries achieved acceptable image quality (mean score range: 3.9-5.0) for both readers. Fair to substantial interobserver agreement was reached on nine arteries. Free-breathing spatiotemporally accelerated 3-D time-resolved contrast-enhanced MR angiography frequently yields diagnostic image quality for most abdominal arteries in young children.

Clinical performance of a free-breathing spatiotemporally accelerated 3-D time-resolved contrast-enhanced pediatric abdominal MR angiography

Energy Technology Data Exchange (ETDEWEB)

Zhang, Tao; Cheng, Joseph Y. [Stanford University, Department of Radiology, Stanford, CA (United States); Stanford University, Department of Electrical Engineering, Stanford, CA (United States); Yousaf, Ufra; Alley, Marcus T.; Vasanawala, Shreyas S. [Stanford University, Department of Radiology, Stanford, CA (United States); Hsiao, Albert [University of California, San Diego, Department of Radiology, San Diego, CA (United States); Lustig, Michael [Stanford University, Department of Electrical Engineering, Stanford, CA (United States); University of California, Berkeley, Department of Electrical Engineering and Computer Sciences, Berkeley, CA (United States); Pauly, John M. [Stanford University, Department of Electrical Engineering, Stanford, CA (United States)

2015-10-15

Pediatric contrast-enhanced MR angiography is often limited by respiration, other patient motion and compromised spatiotemporal resolution. To determine the reliability of a free-breathing spatiotemporally accelerated 3-D time-resolved contrast-enhanced MR angiography method for depicting abdominal arterial anatomy in young children. With IRB approval and informed consent, we retrospectively identified 27 consecutive children (16 males and 11 females; mean age: 3.8 years, range: 14 days to 8.4 years) referred for contrast-enhanced MR angiography at our institution, who had undergone free-breathing spatiotemporally accelerated time-resolved contrast-enhanced MR angiography studies. A radio-frequency-spoiled gradient echo sequence with Cartesian variable density k-space sampling and radial view ordering, intrinsic motion navigation and intermittent fat suppression was developed. Images were reconstructed with soft-gated parallel imaging locally low-rank method to achieve both motion correction and high spatiotemporal resolution. Quality of delineation of 13 abdominal arteries in the reconstructed images was assessed independently by two radiologists on a five-point scale. Ninety-five percent confidence intervals of the proportion of diagnostically adequate cases were calculated. Interobserver agreements were also analyzed. Eleven out of 13 arteries achieved acceptable image quality (mean score range: 3.9-5.0) for both readers. Fair to substantial interobserver agreement was reached on nine arteries. Free-breathing spatiotemporally accelerated 3-D time-resolved contrast-enhanced MR angiography frequently yields diagnostic image quality for most abdominal arteries in young children. (orig.)

Clinical performance of a free-breathing spatiotemporally accelerated 3-D time-resolved contrast-enhanced pediatric abdominal MR angiography

International Nuclear Information System (INIS)

Zhang, Tao; Cheng, Joseph Y.; Yousaf, Ufra; Alley, Marcus T.; Vasanawala, Shreyas S.; Hsiao, Albert; Lustig, Michael; Pauly, John M.

2015-01-01

Pediatric contrast-enhanced MR angiography is often limited by respiration, other patient motion and compromised spatiotemporal resolution. To determine the reliability of a free-breathing spatiotemporally accelerated 3-D time-resolved contrast-enhanced MR angiography method for depicting abdominal arterial anatomy in young children. With IRB approval and informed consent, we retrospectively identified 27 consecutive children (16 males and 11 females; mean age: 3.8 years, range: 14 days to 8.4 years) referred for contrast-enhanced MR angiography at our institution, who had undergone free-breathing spatiotemporally accelerated time-resolved contrast-enhanced MR angiography studies. A radio-frequency-spoiled gradient echo sequence with Cartesian variable density k-space sampling and radial view ordering, intrinsic motion navigation and intermittent fat suppression was developed. Images were reconstructed with soft-gated parallel imaging locally low-rank method to achieve both motion correction and high spatiotemporal resolution. Quality of delineation of 13 abdominal arteries in the reconstructed images was assessed independently by two radiologists on a five-point scale. Ninety-five percent confidence intervals of the proportion of diagnostically adequate cases were calculated. Interobserver agreements were also analyzed. Eleven out of 13 arteries achieved acceptable image quality (mean score range: 3.9-5.0) for both readers. Fair to substantial interobserver agreement was reached on nine arteries. Free-breathing spatiotemporally accelerated 3-D time-resolved contrast-enhanced MR angiography frequently yields diagnostic image quality for most abdominal arteries in young children. (orig.)

CAIPIRINHA accelerated SPACE enables 10-min isotropic 3D TSE MRI of the ankle for optimized visualization of curved and oblique ligaments and tendons.

Science.gov (United States)

Kalia, Vivek; Fritz, Benjamin; Johnson, Rory; Gilson, Wesley D; Raithel, Esther; Fritz, Jan

2017-09-01

To test the hypothesis that a fourfold CAIPIRINHA accelerated, 10-min, high-resolution, isotropic 3D TSE MRI prototype protocol of the ankle derives equal or better quality than a 20-min 2D TSE standard protocol. Following internal review board approval and informed consent, 3-Tesla MRI of the ankle was obtained in 24 asymptomatic subjects including 10-min 3D CAIPIRINHA SPACE TSE prototype and 20-min 2D TSE standard protocols. Outcome variables included image quality and visibility of anatomical structures using 5-point Likert scales. Non-parametric statistical testing was used. P values ≤0.001 were considered significant. Edge sharpness, contrast resolution, uniformity, noise, fat suppression and magic angle effects were without statistical difference on 2D and 3D TSE images (p > 0.035). Fluid was mildly brighter on intermediate-weighted 2D images (p acceleration enables high-spatial resolution oblique and curved planar MRI of the ankle and visualization of ligaments, tendons and joints equally well or better than a more time-consuming anisotropic 2D TSE MRI. • High-resolution 3D TSE MRI improves visualization of ankle structures. • Limitations of current 3D TSE MRI include long scan times. • 3D CAIPIRINHA SPACE allows now a fourfold-accelerated data acquisition. • 3D CAIPIRINHA SPACE enables high-spatial-resolution ankle MRI within 10 min. • 10-min 3D CAIPIRINHA SPACE produces equal-or-better quality than 20-min 2D TSE.

Absolute quantification of myocardial blood flow with 13N-ammonia and 3-dimensional PET.

Science.gov (United States)

Schepis, Tiziano; Gaemperli, Oliver; Treyer, Valerie; Valenta, Ines; Burger, Cyrill; Koepfli, Pascal; Namdar, Mehdi; Adachi, Itaru; Alkadhi, Hatem; Kaufmann, Philipp A

2007-11-01

The aim of this study was to compare 2-dimensional (2D) and 3-dimensional (3D) dynamic PET for the absolute quantification of myocardial blood flow (MBF) with (13)N-ammonia ((13)N-NH(3)). 2D and 3D MBF measurements were collected from 21 patients undergoing cardiac evaluation at rest (n = 14) and during standard adenosine stress (n = 7). A lutetium yttrium oxyorthosilicate-based PET/CT system with retractable septa, enabling the sequential acquisition of 2D and 3D images within the same patient and study, was used. All 2D studies were performed by injecting 700-900 MBq of (13)N-NH(3). For 14 patients, 3D studies were performed with the same injected (13)N-NH(3) dose as that used in 2D studies. For the remaining 7 patients, 3D images were acquired with a lower dose of (13)N-NH(3), that is, 500 MBq. 2D images reconstructed by use of filtered backprojection (FBP) provided the reference standard for MBF measurements. 3D images were reconstructed by use of Fourier rebinning (FORE) with FBP (FORE-FBP), FORE with ordered-subsets expectation maximization (FORE-OSEM), and a reprojection algorithm (RP). Global MBF measurements derived from 3D PET with FORE-FBP (r = 0.97), FORE-OSEM (r = 0.97), and RP (r = 0.97) were well correlated with those derived from 2D FBP (all Ps measurements between 3D FORE-FBP and 2D FBP and between 3D FORE-OSEM and 2D FBP were 0.01 +/- 0.14 and 0.01 +/- 0.15 mL/min/g, respectively. The mean +/- SD difference in global MBF measurements between 3D RP and 2D FBP was 0.00 +/- 0.16 mL/min/g. The best correlation between 2D PET and 3D PET performed with the lower injected activity was found for the 3D FORE-FBP reconstruction algorithm (r = 0.95, P measurements of MBF with 3D PET and (13)N-NH(3) were in excellent agreement with those obtained with the 2D technique, even when a lower activity was injected.

Accelerating image reconstruction in three-dimensional optoacoustic tomography on graphics processing units.

Science.gov (United States)

Wang, Kun; Huang, Chao; Kao, Yu-Jiun; Chou, Cheng-Ying; Oraevsky, Alexander A; Anastasio, Mark A

2013-02-01

Optoacoustic tomography (OAT) is inherently a three-dimensional (3D) inverse problem. However, most studies of OAT image reconstruction still employ two-dimensional imaging models. One important reason is because 3D image reconstruction is computationally burdensome. The aim of this work is to accelerate existing image reconstruction algorithms for 3D OAT by use of parallel programming techniques. Parallelization strategies are proposed to accelerate a filtered backprojection (FBP) algorithm and two different pairs of projection/backprojection operations that correspond to two different numerical imaging models. The algorithms are designed to fully exploit the parallel computing power of graphics processing units (GPUs). In order to evaluate the parallelization strategies for the projection/backprojection pairs, an iterative image reconstruction algorithm is implemented. Computer simulation and experimental studies are conducted to investigate the computational efficiency and numerical accuracy of the developed algorithms. The GPU implementations improve the computational efficiency by factors of 1000, 125, and 250 for the FBP algorithm and the two pairs of projection/backprojection operators, respectively. Accurate images are reconstructed by use of the FBP and iterative image reconstruction algorithms from both computer-simulated and experimental data. Parallelization strategies for 3D OAT image reconstruction are proposed for the first time. These GPU-based implementations significantly reduce the computational time for 3D image reconstruction, complementing our earlier work on 3D OAT iterative image reconstruction.

An implementation of the NiftyRec medical imaging library for PIXE-tomography reconstruction

Science.gov (United States)

Michelet, C.; Barberet, P.; Desbarats, P.; Giovannelli, J.-F.; Schou, C.; Chebil, I.; Delville, M.-H.; Gordillo, N.; Beasley, D. G.; Devès, G.; Moretto, P.; Seznec, H.

2017-08-01

A new development of the TomoRebuild software package is presented, including ;thick sample; correction for non linear X-ray production (NLXP) and X-ray absorption (XA). As in the previous versions, C++ programming with standard libraries was used for easier portability. Data reduction requires different steps which may be run either from a command line instruction or via a user friendly interface, developed as a portable Java plugin in ImageJ. All experimental and reconstruction parameters can be easily modified, either directly in the ASCII parameter files or via the ImageJ interface. A detailed user guide in English is provided. Sinograms and final reconstructed images are generated in usual binary formats that can be read by most public domain graphic softwares. New MLEM and OSEM methods are proposed, using optimized methods from the NiftyRec medical imaging library. An overview of the different medical imaging methods that have been used for ion beam microtomography applications is presented. In TomoRebuild, PIXET data reduction is performed for each chemical element independently and separately from STIMT, except for two steps where the fusion of STIMT and PIXET data is required: the calculation of the correction matrix and the normalization of PIXET data to obtain mass fraction distributions. Correction matrices for NLXP and XA are calculated using procedures extracted from the DISRA code, taking into account a large X-ray detection solid angle. For this, the 3D STIMT mass density distribution is used, considering a homogeneous global composition. A first example of PIXET experiment using two detectors is presented. Reconstruction results are compared and found in good agreement between different codes: FBP, NiftyRec MLEM and OSEM of the TomoRebuild software package, the original DISRA, its accelerated version provided in JPIXET and the accelerated MLEM version of JPIXET, with or without correction.

Measurable realistic image-based 3D mapping

Science.gov (United States)

Liu, W.; Wang, J.; Wang, J. J.; Ding, W.; Almagbile, A.

2011-12-01

Maps with 3D visual models are becoming a remarkable feature of 3D map services. High-resolution image data is obtained for the construction of 3D visualized models.The3D map not only provides the capabilities of 3D measurements and knowledge mining, but also provides the virtual experienceof places of interest, such as demonstrated in the Google Earth. Applications of 3D maps are expanding into the areas of architecture, property management, and urban environment monitoring. However, the reconstruction of high quality 3D models is time consuming, and requires robust hardware and powerful software to handle the enormous amount of data. This is especially for automatic implementation of 3D models and the representation of complicated surfacesthat still need improvements with in the visualisation techniques. The shortcoming of 3D model-based maps is the limitation of detailed coverage since a user can only view and measure objects that are already modelled in the virtual environment. This paper proposes and demonstrates a 3D map concept that is realistic and image-based, that enables geometric measurements and geo-location services. Additionally, image-based 3D maps provide more detailed information of the real world than 3D model-based maps. The image-based 3D maps use geo-referenced stereo images or panoramic images. The geometric relationships between objects in the images can be resolved from the geometric model of stereo images. The panoramic function makes 3D maps more interactive with users but also creates an interesting immersive circumstance. Actually, unmeasurable image-based 3D maps already exist, such as Google street view, but only provide virtual experiences in terms of photos. The topographic and terrain attributes, such as shapes and heights though are omitted. This paper also discusses the potential for using a low cost land Mobile Mapping System (MMS) to implement realistic image 3D mapping, and evaluates the positioning accuracy that a measureable

Applications of 3-D Maxwell solvers to accelerator design

International Nuclear Information System (INIS)

Chou, W.

1990-01-01

This paper gives a brief discussion on various applications of 3-D Maxwell solvers to accelerator design. The work is based on our experience gained during the design of the storage ring of the 7-GeV Advanced Photon Source (APS). It shows that 3-D codes are not replaceable in many cases, and that a lot of work remains to be done in order to establish a solid base for 3-D simulations

MO-G-17A-07: Improved Image Quality in Brain F-18 FDG PET Using Penalized-Likelihood Image Reconstruction Via a Generalized Preconditioned Alternating Projection Algorithm: The First Patient Results

International Nuclear Information System (INIS)

Schmidtlein, CR; Beattie, B; Humm, J; Li, S; Wu, Z; Xu, Y; Zhang, J; Shen, L; Vogelsang, L; Feiglin, D; Krol, A

2014-01-01

Purpose: To investigate the performance of a new penalized-likelihood PET image reconstruction algorithm using the 1 1 -norm total-variation (TV) sum of the 1st through 4th-order gradients as the penalty. Simulated and brain patient data sets were analyzed. Methods: This work represents an extension of the preconditioned alternating projection algorithm (PAPA) for emission-computed tomography. In this new generalized algorithm (GPAPA), the penalty term is expanded to allow multiple components, in this case the sum of the 1st to 4th order gradients, to reduce artificial piece-wise constant regions (“staircase” artifacts typical for TV) seen in PAPA images penalized with only the 1st order gradient. Simulated data were used to test for “staircase” artifacts and to optimize the penalty hyper-parameter in the root-mean-squared error (RMSE) sense. Patient FDG brain scans were acquired on a GE D690 PET/CT (370 MBq at 1-hour post-injection for 10 minutes) in time-of-flight mode and in all cases were reconstructed using resolution recovery projectors. GPAPA images were compared PAPA and RMSE-optimally filtered OSEM (fully converged) in simulations and to clinical OSEM reconstructions (3 iterations, 32 subsets) with 2.6 mm XYGaussian and standard 3-point axial smoothing post-filters. Results: The results from the simulated data show a significant reduction in the 'staircase' artifact for GPAPA compared to PAPA and lower RMSE (up to 35%) compared to optimally filtered OSEM. A simple power-law relationship between the RMSE-optimal hyper-parameters and the noise equivalent counts (NEC) per voxel is revealed. Qualitatively, the patient images appear much sharper and with less noise than standard clinical images. The convergence rate is similar to OSEM. Conclusions: GPAPA reconstructions using the 1 1 -norm total-variation sum of the 1st through 4th-order gradients as the penalty show great promise for the improvement of image quality over that currently

MO-G-17A-07: Improved Image Quality in Brain F-18 FDG PET Using Penalized-Likelihood Image Reconstruction Via a Generalized Preconditioned Alternating Projection Algorithm: The First Patient Results

Energy Technology Data Exchange (ETDEWEB)

Schmidtlein, CR; Beattie, B; Humm, J [Memorial Sloan Kettering Cancer Center, New York, NY (United States); Li, S; Wu, Z; Xu, Y [Sun Yat-sen University, Guangzhou, Guangdong (China); Zhang, J; Shen, L [Syracuse University, Syracuse, NY (United States); Vogelsang, L [VirtualScopics, Rochester, NY (United States); Feiglin, D; Krol, A [SUNY Upstate Medical University, Syracuse, NY (United States)

2014-06-15

Purpose: To investigate the performance of a new penalized-likelihood PET image reconstruction algorithm using the 1{sub 1}-norm total-variation (TV) sum of the 1st through 4th-order gradients as the penalty. Simulated and brain patient data sets were analyzed. Methods: This work represents an extension of the preconditioned alternating projection algorithm (PAPA) for emission-computed tomography. In this new generalized algorithm (GPAPA), the penalty term is expanded to allow multiple components, in this case the sum of the 1st to 4th order gradients, to reduce artificial piece-wise constant regions (“staircase” artifacts typical for TV) seen in PAPA images penalized with only the 1st order gradient. Simulated data were used to test for “staircase” artifacts and to optimize the penalty hyper-parameter in the root-mean-squared error (RMSE) sense. Patient FDG brain scans were acquired on a GE D690 PET/CT (370 MBq at 1-hour post-injection for 10 minutes) in time-of-flight mode and in all cases were reconstructed using resolution recovery projectors. GPAPA images were compared PAPA and RMSE-optimally filtered OSEM (fully converged) in simulations and to clinical OSEM reconstructions (3 iterations, 32 subsets) with 2.6 mm XYGaussian and standard 3-point axial smoothing post-filters. Results: The results from the simulated data show a significant reduction in the 'staircase' artifact for GPAPA compared to PAPA and lower RMSE (up to 35%) compared to optimally filtered OSEM. A simple power-law relationship between the RMSE-optimal hyper-parameters and the noise equivalent counts (NEC) per voxel is revealed. Qualitatively, the patient images appear much sharper and with less noise than standard clinical images. The convergence rate is similar to OSEM. Conclusions: GPAPA reconstructions using the 1{sub 1}-norm total-variation sum of the 1st through 4th-order gradients as the penalty show great promise for the improvement of image quality over that

Single-breath-hold 3-D CINE imaging of the left ventricle using Cartesian sampling.

Science.gov (United States)

Wetzl, Jens; Schmidt, Michaela; Pontana, François; Longère, Benjamin; Lugauer, Felix; Maier, Andreas; Hornegger, Joachim; Forman, Christoph

2018-02-01

Our objectives were to evaluate a single-breath-hold approach for Cartesian 3-D CINE imaging of the left ventricle with a nearly isotropic resolution of [Formula: see text] and a breath-hold duration of [Formula: see text]19 s against a standard stack of 2-D CINE slices acquired in multiple breath-holds. Validation is performed with data sets from ten healthy volunteers. A Cartesian sampling pattern based on the spiral phyllotaxis and a compressed sensing reconstruction method are proposed to allow 3-D CINE imaging with high acceleration factors. The fully integrated reconstruction uses multiple graphics processing units to speed up the reconstruction. The 2-D CINE and 3-D CINE are compared based on ventricular function parameters, contrast-to-noise ratio and edge sharpness measurements. Visual comparisons of corresponding short-axis slices of 2-D and 3-D CINE show an excellent match, while 3-D CINE also allows reformatting to other orientations. Ventricular function parameters do not significantly differ from values based on 2-D CINE imaging. Reconstruction times are below 4 min. We demonstrate single-breath-hold 3-D CINE imaging in volunteers and three example patient cases, which features fast reconstruction and allows reformatting to arbitrary orientations.

Minimally invasive vascular imaging using 3D-CTA and 3D-MRA. Update

International Nuclear Information System (INIS)

Hayashi, Hiromitsu; Kawamata, Hiroshi; Takagi, Ryo; Amano, Yasuo; Wakabayashi, Hiroyuki; Ichikawa, Kazuo; Kumazaki, Tatsuo

1998-01-01

Conventional angiography is considered the standard of reference for diagnostic imaging of vascular diseases with respect to its temporal and spatial resolution. This procedure, however is invasive and repeated studies are difficult, and arterial complications are occasionally associated in catheter-based conventional angiography. Recent advances in diagnostic imaging have facilitated three-dimensional CT angiography (3D-CTA) using the volumetric acquisition capabilities inherent in spiral CT and three-dimensional MR angiography (3D-MRA) using the 3D gradient-echo sequence with a bolus injection of Gd-DTPA. These techniques can provide vascular images exceedingly similar to conventional angiograms within a short acquisition time. 3D-CTA and 3D-MRA are considered to be promising, minimally invasive methods for obtaining images of the vasculature, and alternatives to catheter angiography. This study reviews the current status of 3D-CTA and 3D-MRA, with emphasis on the clinical usefulness of three-dimensional diagnostic imaging for the evaluation of diverse vascular pathologies. (author)

3D images and expert system

International Nuclear Information System (INIS)

Hasegawa, Jun-ichi

1998-01-01

This paper presents an expert system called 3D-IMPRESS for supporting applications of three dimensional (3D) image processing. This system can automatically construct a 3D image processing procedure based on a pictorial example of the goal given by a user. In the paper, to evaluate the performance of the system, it was applied to construction of procedures for extracting specific component figures from practical chest X-ray CT images. (author)

Tomographic spectral imaging: microanalysis in 3D

International Nuclear Information System (INIS)

Kotula, P.G.; Keenan, M.R.; Michael, J.R.

2003-01-01

Full text: Spectral imaging, where a series of complete x-ray spectra are typically collected from a 2D area, holds great promise for comprehensive near-surface microanalysis. There are however numerous microanalysis problems where 3D chemical information is needed as well. In the SEM, some sort of sectioning (either mechanical or with a focused ion beam (FIB) tool) followed by x-ray mapping has, in the past, been utilized in an attempt to perform 3D microanalysis. Reliance on simple mapping has the potential to miss important chemical features as well as misidentify others. In this paper we will describe the acquisition of serial-section tomographic spectral images (TSI) with a dual-beam FIB/SEM equipped with an EDS system. We will also describe the application of a modified version of our multivariate statistical analysis algorithms to TSIs. Serial sectioning was performed with a FEI DB-235 FIB/SEM. Firstly, the specimen normal was tilted to the optic axis of the FIB column and a trench was milled into the surface of the specimen. A second trench was then milled perpendicular to the first to provide visibility of the entire analysis surface to the x-ray detector. In addition, several fiducial markers were milled into the surface to allow for alignment from slice to slice. The electron column is at an angle of 52 deg to the ion column so the electron beam can 'see' the analysis surface milled by the FIB with no additional specimen tilting or rotation. Likewise the x-ray detector is at a radial angle of 45 deg to the plane of the electron and ion columns (about the electron column) and a take-off-angle of 35 deg with respect to an untilted specimen so it can 'see' the analysis surface as well with no additional sample tilting or rotation. Spectral images were acquired from regions 40 μm wide and 20μm deep for each slice. Approximately 1μm/slice was milled and 10-12 total slices were cut. Spectral images were acquired with a Thermo NORAN Vantage (Digital imaging

MO-DE-210-06: Development of a Supercompounded 3D Volumetric Ultrasound Image Guidance System for Prone Accelerated Partial Breast Irradiation (APBI)

Energy Technology Data Exchange (ETDEWEB)

Chiu, T; Hrycushko, B; Zhao, B; Jiang, S; Gu, X [UT Southwestern Medical Center, Dallas, TX (United States)

2015-06-15

Purpose: For early-stage breast cancer, accelerated partial breast irradiation (APBI) is a cost-effective breast-conserving treatment. Irradiation in a prone position can mitigate respiratory induced breast movement and achieve maximal sparing of heart and lung tissues. However, accurate dose delivery is challenging due to breast deformation and lumpectomy cavity shrinkage. We propose a 3D volumetric ultrasound (US) image guidance system for accurate prone APBI Methods: The designed system, set beneath the prone breast board, consists of a water container, an US scanner, and a two-layer breast immobilization cup. The outer layer of the breast cup forms the inner wall of water container while the inner layer is attached to patient breast directly to immobilization. The US transducer scans is attached to the outer-layer of breast cup at the dent of water container. Rotational US scans in a transverse plane are achieved by simultaneously rotating water container and transducer, and multiple transverse scanning forms a 3D scan. A supercompounding-technique-based volumetric US reconstruction algorithm is developed for 3D image reconstruction. The performance of the designed system is evaluated with two custom-made gelatin phantoms containing several cylindrical inserts filled in with water (11% reflection coefficient between materials). One phantom is designed for positioning evaluation while the other is for scaling assessment. Results: In the positioning evaluation phantom, the central distances between the inserts are 15, 20, 30 and 40 mm. The distances on reconstructed images differ by −0.19, −0.65, −0.11 and −1.67 mm, respectively. In the scaling evaluation phantom, inserts are 12.7, 19.05, 25.40 and 31.75 mm in diameter. Measured inserts’ sizes on images differed by 0.23, 0.19, −0.1 and 0.22 mm, respectively. Conclusion: The phantom evaluation results show that the developed 3D volumetric US system can accurately localize target position and determine

CAIPIRINHA accelerated SPACE enables 10-min isotropic 3D TSE MRI of the ankle for optimized visualization of curved and oblique ligaments and tendons

Energy Technology Data Exchange (ETDEWEB)

Kalia, Vivek [University of Vermont Medical Center, Department of Radiology, Burlington, VT (United States); Johns Hopkins University School of Medicine, Russell H. Morgan Department of Radiology and Radiological Science, Section of Musculoskeletal Radiology, Baltimore, MD (United States); Fritz, Benjamin [University Medical Center Freiburg, Department of Radiology, Freiburg im Breisgau (Germany); Johnson, Rory [Siemens Healthcare USA, Inc, Cary, NC (United States); Gilson, Wesley D. [Siemens Healthcare USA, Inc, Baltimore, MD (United States); Raithel, Esther [Siemens Healthcare GmbH, Erlangen (Germany); Fritz, Jan [Johns Hopkins University School of Medicine, Russell H. Morgan Department of Radiology and Radiological Science, Section of Musculoskeletal Radiology, Baltimore, MD (United States)

2017-09-15

To test the hypothesis that a fourfold CAIPIRINHA accelerated, 10-min, high-resolution, isotropic 3D TSE MRI prototype protocol of the ankle derives equal or better quality than a 20-min 2D TSE standard protocol. Following internal review board approval and informed consent, 3-Tesla MRI of the ankle was obtained in 24 asymptomatic subjects including 10-min 3D CAIPIRINHA SPACE TSE prototype and 20-min 2D TSE standard protocols. Outcome variables included image quality and visibility of anatomical structures using 5-point Likert scales. Non-parametric statistical testing was used. P values ≤0.001 were considered significant. Edge sharpness, contrast resolution, uniformity, noise, fat suppression and magic angle effects were without statistical difference on 2D and 3D TSE images (p > 0.035). Fluid was mildly brighter on intermediate-weighted 2D images (p < 0.001), whereas 3D images had substantially less partial volume, chemical shift and no pulsatile-flow artifacts (p < 0.001). Oblique and curved planar 3D images resulted in mildly-to-substantially improved visualization of joints, spring, bifurcate, syndesmotic, collateral and sinus tarsi ligaments, and tendons (p < 0.001, respectively). 3D TSE MRI with CAIPIRINHA acceleration enables high-spatial resolution oblique and curved planar MRI of the ankle and visualization of ligaments, tendons and joints equally well or better than a more time-consuming anisotropic 2D TSE MRI. (orig.)

CAIPIRINHA accelerated SPACE enables 10-min isotropic 3D TSE MRI of the ankle for optimized visualization of curved and oblique ligaments and tendons

International Nuclear Information System (INIS)

Kalia, Vivek; Fritz, Benjamin; Johnson, Rory; Gilson, Wesley D.; Raithel, Esther; Fritz, Jan

2017-01-01

To test the hypothesis that a fourfold CAIPIRINHA accelerated, 10-min, high-resolution, isotropic 3D TSE MRI prototype protocol of the ankle derives equal or better quality than a 20-min 2D TSE standard protocol. Following internal review board approval and informed consent, 3-Tesla MRI of the ankle was obtained in 24 asymptomatic subjects including 10-min 3D CAIPIRINHA SPACE TSE prototype and 20-min 2D TSE standard protocols. Outcome variables included image quality and visibility of anatomical structures using 5-point Likert scales. Non-parametric statistical testing was used. P values ≤0.001 were considered significant. Edge sharpness, contrast resolution, uniformity, noise, fat suppression and magic angle effects were without statistical difference on 2D and 3D TSE images (p > 0.035). Fluid was mildly brighter on intermediate-weighted 2D images (p < 0.001), whereas 3D images had substantially less partial volume, chemical shift and no pulsatile-flow artifacts (p < 0.001). Oblique and curved planar 3D images resulted in mildly-to-substantially improved visualization of joints, spring, bifurcate, syndesmotic, collateral and sinus tarsi ligaments, and tendons (p < 0.001, respectively). 3D TSE MRI with CAIPIRINHA acceleration enables high-spatial resolution oblique and curved planar MRI of the ankle and visualization of ligaments, tendons and joints equally well or better than a more time-consuming anisotropic 2D TSE MRI. (orig.)

Fast GPU-based Monte Carlo code for SPECT/CT reconstructions generates improved 177Lu images.

Science.gov (United States)

Rydén, T; Heydorn Lagerlöf, J; Hemmingsson, J; Marin, I; Svensson, J; Båth, M; Gjertsson, P; Bernhardt, P

2018-01-04

Full Monte Carlo (MC)-based SPECT reconstructions have a strong potential for correcting for image degrading factors, but the reconstruction times are long. The objective of this study was to develop a highly parallel Monte Carlo code for fast, ordered subset expectation maximum (OSEM) reconstructions of SPECT/CT images. The MC code was written in the Compute Unified Device Architecture language for a computer with four graphics processing units (GPUs) (GeForce GTX Titan X, Nvidia, USA). This enabled simulations of parallel photon emissions from the voxels matrix (128 3 or 256 3 ). Each computed tomography (CT) number was converted to attenuation coefficients for photo absorption, coherent scattering, and incoherent scattering. For photon scattering, the deflection angle was determined by the differential scattering cross sections. An angular response function was developed and used to model the accepted angles for photon interaction with the crystal, and a detector scattering kernel was used for modeling the photon scattering in the detector. Predefined energy and spatial resolution kernels for the crystal were used. The MC code was implemented in the OSEM reconstruction of clinical and phantom 177 Lu SPECT/CT images. The Jaszczak image quality phantom was used to evaluate the performance of the MC reconstruction in comparison with attenuated corrected (AC) OSEM reconstructions and attenuated corrected OSEM reconstructions with resolution recovery corrections (RRC). The performance of the MC code was 3200 million photons/s. The required number of photons emitted per voxel to obtain a sufficiently low noise level in the simulated image was 200 for a 128 3 voxel matrix. With this number of emitted photons/voxel, the MC-based OSEM reconstruction with ten subsets was performed within 20 s/iteration. The images converged after around six iterations. Therefore, the reconstruction time was around 3 min. The activity recovery for the spheres in the Jaszczak phantom was

Nonlaser-based 3D surface imaging

Energy Technology Data Exchange (ETDEWEB)

Lu, Shin-yee; Johnson, R.K.; Sherwood, R.J. [Lawrence Livermore National Lab., CA (United States)

1994-11-15

3D surface imaging refers to methods that generate a 3D surface representation of objects of a scene under viewing. Laser-based 3D surface imaging systems are commonly used in manufacturing, robotics and biomedical research. Although laser-based systems provide satisfactory solutions for most applications, there are situations where non laser-based approaches are preferred. The issues that make alternative methods sometimes more attractive are: (1) real-time data capturing, (2) eye-safety, (3) portability, and (4) work distance. The focus of this presentation is on generating a 3D surface from multiple 2D projected images using CCD cameras, without a laser light source. Two methods are presented: stereo vision and depth-from-focus. Their applications are described.

3-D image reconstruction in radiology

International Nuclear Information System (INIS)

Grangeat, P.

1999-01-01

In this course, we present highlights on fully 3-D image reconstruction algorithms used in 3-D X-ray Computed Tomography (3-D-CT) and 3-D Rotational Radiography (3-D-RR). We first consider the case of spiral CT with a one-row detector. Starting from the 2-D fan-beam inversion formula for a circular trajectory, we introduce spiral CT 3-D image reconstruction algorithm using axial interpolation for each transverse slice. In order to improve the X-ray detection efficiency and to speed the acquisition process, the future is to use multi-row detectors associated with small angle cone-beam geometry. The generalization of the 2-D fan-beam image reconstruction algorithm to cone beam defined direct inversion formula referred as Feldkamp's algorithm for a circular trajectory and Wang's algorithm for a spiral trajectory. However, large area detectors does exist such as Radiological Image Intensifiers or in a near future solid state detectors. To get a larger zoom effect, it defines a cone-beam geometry associated with a large aperture angle. For this case, we introduce indirect image reconstruction algorithm by plane re-binning in the Radon domain. We will present some results from a prototype MORPHOMETER device using the RADON reconstruction software. Lastly, we consider the special case of 3-D Rotational Digital Subtraction Angiography with a restricted number of views. We introduce constraint optimization algorithm using quadratic, entropic or half-quadratic constraints. Generalized ART (Algebraic Reconstruction Technique) iterative reconstruction algorithm can be derived from the Bregman algorithm. We present reconstructed vascular trees from a prototype MORPHOMETER device. (author)

Image-based RSA: Roentgen stereophotogrammetric analysis based on 2D-3D image registration.

Science.gov (United States)

de Bruin, P W; Kaptein, B L; Stoel, B C; Reiber, J H C; Rozing, P M; Valstar, E R

2008-01-01

Image-based Roentgen stereophotogrammetric analysis (IBRSA) integrates 2D-3D image registration and conventional RSA. Instead of radiopaque RSA bone markers, IBRSA uses 3D CT data, from which digitally reconstructed radiographs (DRRs) are generated. Using 2D-3D image registration, the 3D pose of the CT is iteratively adjusted such that the generated DRRs resemble the 2D RSA images as closely as possible, according to an image matching metric. Effectively, by registering all 2D follow-up moments to the same 3D CT, the CT volume functions as common ground. In two experiments, using RSA and using a micromanipulator as gold standard, IBRSA has been validated on cadaveric and sawbone scapula radiographs, and good matching results have been achieved. The accuracy was: |mu |RSA but higher than in vivo standard RSA. Because IBRSA does not require radiopaque markers, it adds functionality to the RSA method by opening new directions and possibilities for research, such as dynamic analyses using fluoroscopy on subjects without markers and computer navigation applications.

Complex adaptation-based LDR image rendering for 3D image reconstruction

Science.gov (United States)

Lee, Sung-Hak; Kwon, Hyuk-Ju; Sohng, Kyu-Ik

2014-07-01

A low-dynamic tone-compression technique is developed for realistic image rendering that can make three-dimensional (3D) images similar to realistic scenes by overcoming brightness dimming in the 3D display mode. The 3D surround provides varying conditions for image quality, illuminant adaptation, contrast, gamma, color, sharpness, and so on. In general, gain/offset adjustment, gamma compensation, and histogram equalization have performed well in contrast compression; however, as a result of signal saturation and clipping effects, image details are removed and information is lost on bright and dark areas. Thus, an enhanced image mapping technique is proposed based on space-varying image compression. The performance of contrast compression is enhanced with complex adaptation in a 3D viewing surround combining global and local adaptation. Evaluating local image rendering in view of tone and color expression, noise reduction, and edge compensation confirms that the proposed 3D image-mapping model can compensate for the loss of image quality in the 3D mode.

Multiplexed phase-space imaging for 3D fluorescence microscopy.

Science.gov (United States)

Liu, Hsiou-Yuan; Zhong, Jingshan; Waller, Laura

2017-06-26

Optical phase-space functions describe spatial and angular information simultaneously; examples of optical phase-space functions include light fields in ray optics and Wigner functions in wave optics. Measurement of phase-space enables digital refocusing, aberration removal and 3D reconstruction. High-resolution capture of 4D phase-space datasets is, however, challenging. Previous scanning approaches are slow, light inefficient and do not achieve diffraction-limited resolution. Here, we propose a multiplexed method that solves these problems. We use a spatial light modulator (SLM) in the pupil plane of a microscope in order to sequentially pattern multiplexed coded apertures while capturing images in real space. Then, we reconstruct the 3D fluorescence distribution of our sample by solving an inverse problem via regularized least squares with a proximal accelerated gradient descent solver. We experimentally reconstruct a 101 Megavoxel 3D volume (1010×510×500µm with NA 0.4), demonstrating improved acquisition time, light throughput and resolution compared to scanning aperture methods. Our flexible patterning scheme further allows sparsity in the sample to be exploited for reduced data capture.

2D-3D image registration in diagnostic and interventional X-Ray imaging

NARCIS (Netherlands)

Bom, I.M.J. van der

2010-01-01

Clinical procedures that are conventionally guided by 2D x-ray imaging, may benefit from the additional spatial information provided by 3D image data. For instance, guidance of minimally invasive procedures with CT or MRI data provides 3D spatial information and visualization of structures that are

Contributions in compression of 3D medical images and 2D images

International Nuclear Information System (INIS)

Gaudeau, Y.

2006-12-01

The huge amounts of volumetric data generated by current medical imaging techniques in the context of an increasing demand for long term archiving solutions, as well as the rapid development of distant radiology make the use of compression inevitable. Indeed, if the medical community has sided until now with compression without losses, most of applications suffer from compression ratios which are too low with this kind of compression. In this context, compression with acceptable losses could be the most appropriate answer. So, we propose a new loss coding scheme based on 3D (3 dimensional) Wavelet Transform and Dead Zone Lattice Vector Quantization 3D (DZLVQ) for medical images. Our algorithm has been evaluated on several computerized tomography (CT) and magnetic resonance image volumes. The main contribution of this work is the design of a multidimensional dead zone which enables to take into account correlations between neighbouring elementary volumes. At high compression ratios, we show that it can out-perform visually and numerically the best existing methods. These promising results are confirmed on head CT by two medical patricians. The second contribution of this document assesses the effect with-loss image compression on CAD (Computer-Aided Decision) detection performance of solid lung nodules. This work on 120 significant lungs images shows that detection did not suffer until 48:1 compression and still was robust at 96:1. The last contribution consists in the complexity reduction of our compression scheme. The first allocation dedicated to 2D DZLVQ uses an exponential of the rate-distortion (R-D) functions. The second allocation for 2D and 3D medical images is based on block statistical model to estimate the R-D curves. These R-D models are based on the joint distribution of wavelet vectors using a multidimensional mixture of generalized Gaussian (MMGG) densities. (author)

Optimization of image reconstruction method for SPECT studies performed using [⁹⁹mTc-EDDA/HYNIC] octreotate in patients with neuroendocrine tumors.

Science.gov (United States)

Sowa-Staszczak, Anna; Lenda-Tracz, Wioletta; Tomaszuk, Monika; Głowa, Bogusław; Hubalewska-Dydejczyk, Alicja

2013-01-01

Somatostatin receptor scintigraphy (SRS) is a useful tool in the assessment of GEP-NET (gastroenteropancreatic neuroendocrine tumor) patients. The choice of appropriate settings of image reconstruction parameters is crucial in interpretation of these images. The aim of the study was to investigate how the GEP NET lesion signal to noise ratio (TCS/TCB) depends on different reconstruction settings for Flash 3D software (Siemens). SRS results of 76 randomly selected patients with confirmed GEP-NET were analyzed. For SPECT studies the data were acquired using standard clinical settings 3-4 h after the injection of 740 MBq 99mTc-[EDDA/HYNIC] octreotate. To obtain final images the OSEM 3D Flash reconstruction with different settings and FBP reconstruction were used. First, the TCS/TCB ratio in voxels was analyzed for different combinations of the number of subsets and the number of iterations of the OSEM 3D Flash reconstruction. Secondly, the same ratio was analyzed for different parameters of the Gaussian filter (with FWHM = 2-4 times greater from the pixel size). Also the influence of scatter correction on the TCS/TCB ratio was investigated. With increasing number of subsets and iterations, the increase of TCS/TCB ratio was observed. With increasing settings of Gauss [FWHM coefficient] filter, the decrease of TCS/TCB ratio was reported. The use of scatter correction slightly decreases the values of this ratio. OSEM algorithm provides a meaningfully better reconstruction of the SRS SPECT study as compared to the FBP technique. A high number of subsets improves image quality (images are smoother). Increasing number of iterations gives a better contrast and the shapes of lesions and organs are sharper. The choice of reconstruction parameters is a compromise between image qualitative appearance and its quantitative accuracy and should not be modified when comparing multiple studies of the same patient.

Advanced 3-D Ultrasound Imaging: 3-D Synthetic Aperture Imaging using Fully Addressed and Row-Column Addressed 2-D Transducer Arrays

DEFF Research Database (Denmark)

Bouzari, Hamed

the important diagnostic information in a noninvasive manner. Diagnostic and therapeutic decisions often require accurate estimates of e.g., organ, cyst, or tumor volumes. 3-D ultrasound imaging can provide these measurements without relying on the geometrical assumptions and operator-dependent skills involved...... is one of the factors for the widespread use of ultrasound imaging. The high price tag on the high quality 3-D scanners is limiting their market share. Row-column addressing of 2-D transducer arrays is a low cost alternative to fully addressed 2-D arrays, for 3-D ultrasound imaging. Using row....... Based on a set of acoustical measurements the center frequency, bandwidth, surface pressure, sensitivity, and acoustical cross-talks were evaluated and discussed. The imaging quality assessments were carried out based on Field II simulations as well as phantom measurements. Moreover, an analysis...

Extracting 3D layout from a single image using global image structures.

Science.gov (United States)

Lou, Zhongyu; Gevers, Theo; Hu, Ninghang

2015-10-01

Extracting the pixel-level 3D layout from a single image is important for different applications, such as object localization, image, and video categorization. Traditionally, the 3D layout is derived by solving a pixel-level classification problem. However, the image-level 3D structure can be very beneficial for extracting pixel-level 3D layout since it implies the way how pixels in the image are organized. In this paper, we propose an approach that first predicts the global image structure, and then we use the global structure for fine-grained pixel-level 3D layout extraction. In particular, image features are extracted based on multiple layout templates. We then learn a discriminative model for classifying the global layout at the image-level. Using latent variables, we implicitly model the sublevel semantics of the image, which enrich the expressiveness of our model. After the image-level structure is obtained, it is used as the prior knowledge to infer pixel-wise 3D layout. Experiments show that the results of our model outperform the state-of-the-art methods by 11.7% for 3D structure classification. Moreover, we show that employing the 3D structure prior information yields accurate 3D scene layout segmentation.

Evaluation of an accelerated 3D SPACE sequence with compressed sensing and free-stop scan mode for imaging of the knee.

Science.gov (United States)

Henninger, B; Raithel, E; Kranewitter, C; Steurer, M; Jaschke, W; Kremser, C

2018-05-01

To prospectively evaluate a prototypical 3D turbo-spin-echo proton-density-weighted sequence with compressed sensing and free-stop scan mode for preventing motion artefacts (3D-PD-CS-SPACE free-stop) for knee imaging in a clinical setting. 80 patients underwent 3T magnetic resonance imaging (MRI) of the knee with our 2D routine protocol and with 3D-PD-CS-SPACE free-stop. In case of a scan-stop caused by motion (images are calculated nevertheless) the sequence was repeated without free-stop mode. All scans were evaluated by 2 radiologists concerning image quality of the 3D-PD-CS-SPACE (with and without free-stop). Important knee structures were further assessed in a lesion based analysis and compared to our reference 2D-PD-fs sequences. Image quality of the 3D-PD-CS-SPACE free-stop was found optimal in 47/80, slightly compromised in 21/80, moderately in 10/80 and severely in 2/80. In 29/80, the free-stop scan mode stopped the 3D-PD-CS-SPACE due to subject motion with a slight increase of image quality at longer effective acquisition times. Compared to the 3D-PD-CS-SPACE with free-stop, the image quality of the acquired 3D-PD-CS-SPACE without free-stop was found equal in 6/29, slightly improved in 13/29, improved with equal contours in 8/29, and improved with sharper contours in 2/29. The lesion based analysis showed a high agreement between the results from the 3D-PD-CS-SPACE free-stop and our 2D-PD-fs routine protocol (overall agreement 96.25%-100%, Cohen's Kappa 0.883-1, p SPACE free-stop is a reliable alternative for standard 2D-PD-fs protocols with acceptable acquisition times. Copyright © 2018 Elsevier B.V. All rights reserved.

Preliminary examples of 3D vector flow imaging

DEFF Research Database (Denmark)

Pihl, Michael Johannes; Stuart, Matthias Bo; Tomov, Borislav Gueorguiev

2013-01-01

This paper presents 3D vector flow images obtained using the 3D Transverse Oscillation (TO) method. The method employs a 2D transducer and estimates the three velocity components simultaneously, which is important for visualizing complex flow patterns. Data are acquired using the experimental ult...... as opposed to magnetic resonance imaging (MRI). The results demonstrate that the 3D TO method is capable of performing 3D vector flow imaging.......This paper presents 3D vector flow images obtained using the 3D Transverse Oscillation (TO) method. The method employs a 2D transducer and estimates the three velocity components simultaneously, which is important for visualizing complex flow patterns. Data are acquired using the experimental...... ultrasound scanner SARUS on a flow rig system with steady flow. The vessel of the flow-rig is centered at a depth of 30 mm, and the flow has an expected 2D circular-symmetric parabolic prole with a peak velocity of 1 m/s. Ten frames of 3D vector flow images are acquired in a cross-sectional plane orthogonal...

Fast implementations of reconstruction-based scatter compensation in fully 3D SPECT image reconstruction

International Nuclear Information System (INIS)

Kadrmas, Dan J.; Karimi, Seemeen S.; Frey, Eric C.; Tsui, Benjamin M.W.

1998-01-01

Accurate scatter compensation in SPECT can be performed by modelling the scatter response function during the reconstruction process. This method is called reconstruction-based scatter compensation (RBSC). It has been shown that RBSC has a number of advantages over other methods of compensating for scatter, but using RBSC for fully 3D compensation has resulted in prohibitively long reconstruction times. In this work we propose two new methods that can be used in conjunction with existing methods to achieve marked reductions in RBSC reconstruction times. The first method, coarse-grid scatter modelling, significantly accelerates the scatter model by exploiting the fact that scatter is dominated by low-frequency information. The second method, intermittent RBSC, further accelerates the reconstruction process by limiting the number of iterations during which scatter is modelled. The fast implementations were evaluated using a Monte Carlo simulated experiment of the 3D MCAT phantom with 99m Tc tracer, and also using experimentally acquired data with 201 Tl tracer. Results indicated that these fast methods can reconstruct, with fully 3D compensation, images very similar to those obtained using standard RBSC methods, and in reconstruction times that are an order of magnitude shorter. Using these methods, fully 3D iterative reconstruction with RBSC can be performed well within the realm of clinically realistic times (under 10 minutes for 64x64x24 image reconstruction). (author)

GPU-Based 3D Cone-Beam CT Image Reconstruction for Large Data Volume

Directory of Open Access Journals (Sweden)

Xing Zhao

2009-01-01

Full Text Available Currently, 3D cone-beam CT image reconstruction speed is still a severe limitation for clinical application. The computational power of modern graphics processing units (GPUs has been harnessed to provide impressive acceleration of 3D volume image reconstruction. For extra large data volume exceeding the physical graphic memory of GPU, a straightforward compromise is to divide data volume into blocks. Different from the conventional Octree partition method, a new partition scheme is proposed in this paper. This method divides both projection data and reconstructed image volume into subsets according to geometric symmetries in circular cone-beam projection layout, and a fast reconstruction for large data volume can be implemented by packing the subsets of projection data into the RGBA channels of GPU, performing the reconstruction chunk by chunk and combining the individual results in the end. The method is evaluated by reconstructing 3D images from computer-simulation data and real micro-CT data. Our results indicate that the GPU implementation can maintain original precision and speed up the reconstruction process by 110–120 times for circular cone-beam scan, as compared to traditional CPU implementation.

3D FaceCam: a fast and accurate 3D facial imaging device for biometrics applications

Science.gov (United States)

Geng, Jason; Zhuang, Ping; May, Patrick; Yi, Steven; Tunnell, David

2004-08-01

Human faces are fundamentally three-dimensional (3D) objects, and each face has its unique 3D geometric profile. The 3D geometric features of a human face can be used, together with its 2D texture, for rapid and accurate face recognition purposes. Due to the lack of low-cost and robust 3D sensors and effective 3D facial recognition (FR) algorithms, almost all existing FR systems use 2D face images. Genex has developed 3D solutions that overcome the inherent problems in 2D while also addressing limitations in other 3D alternatives. One important aspect of our solution is a unique 3D camera (the 3D FaceCam) that combines multiple imaging sensors within a single compact device to provide instantaneous, ear-to-ear coverage of a human face. This 3D camera uses three high-resolution CCD sensors and a color encoded pattern projection system. The RGB color information from each pixel is used to compute the range data and generate an accurate 3D surface map. The imaging system uses no moving parts and combines multiple 3D views to provide detailed and complete 3D coverage of the entire face. Images are captured within a fraction of a second and full-frame 3D data is produced within a few seconds. This described method provides much better data coverage and accuracy in feature areas with sharp features or details (such as the nose and eyes). Using this 3D data, we have been able to demonstrate that a 3D approach can significantly improve the performance of facial recognition. We have conducted tests in which we have varied the lighting conditions and angle of image acquisition in the "field." These tests have shown that the matching results are significantly improved when enrolling a 3D image rather than a single 2D image. With its 3D solutions, Genex is working toward unlocking the promise of powerful 3D FR and transferring FR from a lab technology into a real-world biometric solution.

Nonthermal Particle Acceleration in 3D Relativistic Magnetic Reconnection in Pair Plasma

Energy Technology Data Exchange (ETDEWEB)

Werner, Gregory R.; Uzdensky, Dmitri A., E-mail: Greg.Werner@colorado.edu [Center for Integrated Plasma Studies, Physics Department, University of Colorado, 390 UCB, Boulder, CO 80309 (United States)

2017-07-10

As a fundamental process converting magnetic to plasma energy in high-energy astrophysical plasmas, relativistic magnetic reconnection is a leading explanation for the acceleration of particles to the ultrarelativistic energies that are necessary to power nonthermal emission (especially X-rays and gamma-rays) in pulsar magnetospheres and pulsar wind nebulae, coronae and jets of accreting black holes, and gamma-ray bursts. An important objective of plasma astrophysics is therefore the characterization of nonthermal particle acceleration (NTPA) effected by reconnection. Reconnection-powered NTPA has been demonstrated over a wide range of physical conditions using large 2D kinetic simulations. However, its robustness in realistic 3D reconnection—in particular, whether the 3D relativistic drift-kink instability (RDKI) disrupts NTPA—has not been systematically investigated, although pioneering 3D simulations have observed NTPA in isolated cases. Here, we present the first comprehensive study of NTPA in 3D relativistic reconnection in collisionless electron–positron plasmas, characterizing NTPA as the strength of 3D effects is varied systematically via the length in the third dimension and the strength of the guide magnetic field. We find that, while the RDKI prominently perturbs 3D reconnecting current sheets, it does not suppress particle acceleration, even for zero guide field; fully 3D reconnection robustly and efficiently produces nonthermal power-law particle spectra closely resembling those obtained in 2D. This finding provides strong support for reconnection as the key mechanism powering high-energy flares in various astrophysical systems. We also show that strong guide fields significantly inhibit NTPA, slowing reconnection and limiting the energy available for plasma energization, yielding steeper and shorter power-law spectra.

Accelerated high-resolution 3D magnetic resonance spectroscopic imaging in the brain At 7 T

International Nuclear Information System (INIS)

Hangel, G.

2015-01-01

With the announcement of the first series of magnetic resonance (MR) scanners with a field strength of 7 Tesla (T) intended for clinical practice, the development of high-performance sequences for higher field strengths has gained importance. Magnetic resonance spectroscopic imaging (MRSI) in the brain currently offers the unique ability to spatially resolve the distribution of multiple metabolites simultaneously. Its big diagnostic potential could be applied to many clinical protocols, for example the assessment of tumour treatment or progress of Multiple Sclerosis. Moving to ultra-high fields like 7 T has the main benefits of increased signal-to-noise ratio (SNR) and improved spectral quality, but brings its own challenges due to stronger field inhomogeneities. Necessary for a robust, flexible and useful MRSI sequence in the brain are high resolutions, shortened measurement times, the possibility for 3D-MRSI and the suppression of spectral contamination by trans-cranial lipids. This thesis addresses these limitations and proposes Hadamard spectroscopic imaging (HSI) as solution for multi-slice MRSI, the application of generalized autocalibrating partially parallel acquisition (GRAPPA) and spiral trajectories for measurement acceleration, non-selective inversion recovery (IR) lipid-suppression as well as combinations of these methods. Further, the optimisation of water suppression for 7 T systems and the acquisition of ultra-high resolution (UHR)-MRSI are discussed. In order to demonstrate the clinical feasibility of these approaches, MRSI measurement results of a glioma patient are presented. The discussion of the obtained results in the context of the state-of-art in 7 T MRSI in the brain, possible future applications as well as potential further improvements of the MRSI sequences conclude this thesis. (author) [de

D3D augmented reality imaging system: proof of concept in mammography.

Science.gov (United States)

Douglas, David B; Petricoin, Emanuel F; Liotta, Lance; Wilson, Eugene

2016-01-01

The purpose of this article is to present images from simulated breast microcalcifications and assess the pattern of the microcalcifications with a technical development called "depth 3-dimensional (D3D) augmented reality". A computer, head display unit, joystick, D3D augmented reality software, and an in-house script of simulated data of breast microcalcifications in a ductal distribution were used. No patient data was used and no statistical analysis was performed. The D3D augmented reality system demonstrated stereoscopic depth perception by presenting a unique image to each eye, focal point convergence, head position tracking, 3D cursor, and joystick fly-through. The D3D augmented reality imaging system offers image viewing with depth perception and focal point convergence. The D3D augmented reality system should be tested to determine its utility in clinical practice.

Efficient methodologies for system matrix modelling in iterative image reconstruction for rotating high-resolution PET

Energy Technology Data Exchange (ETDEWEB)

Ortuno, J E; Kontaxakis, G; Rubio, J L; Santos, A [Departamento de Ingenieria Electronica (DIE), Universidad Politecnica de Madrid, Ciudad Universitaria s/n, 28040 Madrid (Spain); Guerra, P [Networking Research Center on Bioengineering, Biomaterials and Nanomedicine (CIBER-BBN), Madrid (Spain)], E-mail: juanen@die.upm.es

2010-04-07

A fully 3D iterative image reconstruction algorithm has been developed for high-resolution PET cameras composed of pixelated scintillator crystal arrays and rotating planar detectors, based on the ordered subsets approach. The associated system matrix is precalculated with Monte Carlo methods that incorporate physical effects not included in analytical models, such as positron range effects and interaction of the incident gammas with the scintillator material. Custom Monte Carlo methodologies have been developed and optimized for modelling of system matrices for fast iterative image reconstruction adapted to specific scanner geometries, without redundant calculations. According to the methodology proposed here, only one-eighth of the voxels within two central transaxial slices need to be modelled in detail. The rest of the system matrix elements can be obtained with the aid of axial symmetries and redundancies, as well as in-plane symmetries within transaxial slices. Sparse matrix techniques for the non-zero system matrix elements are employed, allowing for fast execution of the image reconstruction process. This 3D image reconstruction scheme has been compared in terms of image quality to a 2D fast implementation of the OSEM algorithm combined with Fourier rebinning approaches. This work confirms the superiority of fully 3D OSEM in terms of spatial resolution, contrast recovery and noise reduction as compared to conventional 2D approaches based on rebinning schemes. At the same time it demonstrates that fully 3D methodologies can be efficiently applied to the image reconstruction problem for high-resolution rotational PET cameras by applying accurate pre-calculated system models and taking advantage of the system's symmetries.

Optimization, evaluation, and comparison of standard algorithms for image reconstruction with the VIP-PET.

Science.gov (United States)

Mikhaylova, E; Kolstein, M; De Lorenzo, G; Chmeissani, M

2014-07-01

A novel positron emission tomography (PET) scanner design based on a room-temperature pixelated CdTe solid-state detector is being developed within the framework of the Voxel Imaging PET (VIP) Pathfinder project [1]. The simulation results show a great potential of the VIP to produce high-resolution images even in extremely challenging conditions such as the screening of a human head [2]. With unprecedented high channel density (450 channels/cm 3 ) image reconstruction is a challenge. Therefore optimization is needed to find the best algorithm in order to exploit correctly the promising detector potential. The following reconstruction algorithms are evaluated: 2-D Filtered Backprojection (FBP), Ordered Subset Expectation Maximization (OSEM), List-Mode OSEM (LM-OSEM), and the Origin Ensemble (OE) algorithm. The evaluation is based on the comparison of a true image phantom with a set of reconstructed images obtained by each algorithm. This is achieved by calculation of image quality merit parameters such as the bias, the variance and the mean square error (MSE). A systematic optimization of each algorithm is performed by varying the reconstruction parameters, such as the cutoff frequency of the noise filters and the number of iterations. The region of interest (ROI) analysis of the reconstructed phantom is also performed for each algorithm and the results are compared. Additionally, the performance of the image reconstruction methods is compared by calculating the modulation transfer function (MTF). The reconstruction time is also taken into account to choose the optimal algorithm. The analysis is based on GAMOS [3] simulation including the expected CdTe and electronic specifics.

Quantitative evaluation of temporal partial coherence using 3D Fourier transforms of through-focus TEM images

International Nuclear Information System (INIS)

Kimoto, Koji; Sawada, Hidetaka; Sasaki, Takeo; Sato, Yuta; Nagai, Takuro; Ohwada, Megumi; Suenaga, Kazu; Ishizuka, Kazuo

2013-01-01

We evaluate the temporal partial coherence of transmission electron microscopy (TEM) using the three-dimensional (3D) Fourier transform (FT) of through-focus images. Young's fringe method often indicates the unexpected high-frequency information due to non-linear imaging terms. We have already used the 3D FT of axial (non-tilted) through-focus images to reduce the effect of non-linear terms on the linear imaging term, and demonstrated the improvement of monochromated lower-voltage TEM performance [Kimoto et al., Ultramicroscopy 121 (2012) 31–39]. Here we apply the 3D FT method with intentionally tilted incidence to normalize various factors associated with a TEM specimen and an imaging device. The temporal partial coherence of two microscopes operated at 30, 60 and 80 kV is evaluated. Our method is applicable to such cases where the non-linear terms become more significant in lower acceleration voltage or aberration-corrected high spatial resolution TEM. - Highlights: • We assess the temporal partial coherence of TEM using a 3-dimensional (3D) Fourier transform (FT) of through-focus images. • We apply the 3D FT method with intentionally tilted incidence to normalize various factors associated with a TEM specimen and an imaging device. • The spatial frequency at which information transfer decreases to 1/e 2 (13.5%) is determined for two lower-voltage TEM systems

Comparison of modern 3D and 2D MR imaging sequences of the wrist at 3 Tesla

International Nuclear Information System (INIS)

Rehnitz, C.; Klaan, B.; Amarteifio, E.; Kauczor, H.U.; Weber, M.A.; Stillfried, F. von; Burkholder, I.

2016-01-01

To compare the image quality of modern 3 D and 2 D sequences for dedicated wrist imaging at 3 Tesla (T) MRI. At 3 T MRI, 18 patients (mean age: 36.2 years) with wrist pain and 16 healthy volunteers (mean age: 26.4 years) were examined using 2 D proton density-weighted fat-saturated (PDfs), isotropic 3 D TrueFISP, 3 D MEDIC, and 3 D PDfs SPACE sequences. Image quality was rated on a five-point scale (0 - 4) including overall image quality (OIQ), visibility of important structures (cartilage, ligaments, TFCC) and degree of artifacts. Signal-to-noise ratios (SNR) and contrast-to-noise ratios (CNR) of cartilage/bone/muscle/fluid as well as the mean overall SNR/CNR were calculated using region-of-interest analysis. ANOVA, paired t-, and Wilcoxon-signed-rank tests were applied. The image quality of all tested sequences was superior to 3 D PDfs SPACE (p < 0.01). 3 D TrueFISP had the highest combined cartilage score (mean: 3.4) and performed better in cartilage comparisons against 3 D PDfs SPACE in both groups and 2 D PDfs in volunteers (p < 0.05). 3 D MEDIC performed better in 7 of 8 comparisons (p < 0.05) regarding ligaments and TFCC. 2 D PDfs provided constantly high scores. The mean overall SNR/CNR for 2 D PDfs, 3 D PDfs SPACE, 3 D TrueFISP, and 3 D MEDIC were 68/65, 32/27, 45/47, and 57/45, respectively. 2 D PDfs performed best in most SNR/CNR comparisons (p < 0.05) and 3 D MEDIC performed best within the 3 D sequences (p < 0.05). Except 3 D PDfs SPACE, all tested 3 D and 2 D sequences provided high image quality. 3 D TrueFISP was best for cartilage imaging, 3 D MEDIC for ligaments and TFCC and 2 D PDfs for general wrist imaging.

Dual matrix ordered subsets reconstruction for accelerated 3D scatter compensation in single-photon emission tomography

International Nuclear Information System (INIS)

Kamphuis, C.; Beekman, F.J.; Van Rijk, P.P.; Viergever, M.A.

1998-01-01

Three-dimensional (3D) iterative maximum likelihood expectation maximization (ML-EM) algorithms for single-photon emission tomography (SPET) are capable of correcting image-degrading effects of non-uniform attenuation, distance-dependent camera response and patient shape-dependent scatter. However, the resulting improvements in quantitation, resolution and signal-to-noise ratio (SNR) are obtained at the cost of a huge computational burden. This paper presents a new acceleration method for ML-EM: dual matrix ordered subsets (DM-OS). DM-OS combines two acceleration methods: (a) different matrices for projection and back-projection and (b) ordered subsets of projections. DM-OS was compared with ML-EM on simulated data and on physical thorax phantom data, for both 180 and 360 orbits. Contrast, normalized standard deviation and mean squared error were calculated for the digital phantom experiment. DM-OS resulted in similar image quality to ML-EM, even for speed-up factors of 200 compared to ML-EM in the case of 120 projections. The thorax phantom data could be reconstructed 50 times faster (60 projections) using DM-OS with preservation of image quality. ML-EM and DM-OS with scatter compensation showed significant improvement of SNR compared to ML-EM without scatter compensation. Furthermore, inclusion of complex image formation models in the computer code is simplified in the case of DM-OS. It is thus shown that DM-OS is a fast and relatively simple algorithm for 3D iterative scatter compensation, with similar results to conventional ML-EM, for both 180 and 360 acquired data. (orig.)

On Alternative Approaches to 3D Image Perception: Monoscopic 3D Techniques

Science.gov (United States)

Blundell, Barry G.

2015-06-01

In the eighteenth century, techniques that enabled a strong sense of 3D perception to be experienced without recourse to binocular disparities (arising from the spatial separation of the eyes) underpinned the first significant commercial sales of 3D viewing devices and associated content. However following the advent of stereoscopic techniques in the nineteenth century, 3D image depiction has become inextricably linked to binocular parallax and outside the vision science and arts communities relatively little attention has been directed towards earlier approaches. Here we introduce relevant concepts and terminology and consider a number of techniques and optical devices that enable 3D perception to be experienced on the basis of planar images rendered from a single vantage point. Subsequently we allude to possible mechanisms for non-binocular parallax based 3D perception. Particular attention is given to reviewing areas likely to be thought-provoking to those involved in 3D display development, spatial visualization, HCI, and other related areas of interdisciplinary research.

D3D augmented reality imaging system: proof of concept in mammography

Directory of Open Access Journals (Sweden)

Douglas DB

2016-08-01

Full Text Available David B Douglas,1 Emanuel F Petricoin,2 Lance Liotta,2 Eugene Wilson3 1Department of Radiology, Stanford University, Palo Alto, CA, 2Center for Applied Proteomics and Molecular Medicine, George Mason University, Manassas, VA, 3Department of Radiology, Fort Benning, Columbus, GA, USA Purpose: The purpose of this article is to present images from simulated breast microcalcifications and assess the pattern of the microcalcifications with a technical development called “depth 3-dimensional (D3D augmented reality”. Materials and methods: A computer, head display unit, joystick, D3D augmented reality software, and an in-house script of simulated data of breast microcalcifications in a ductal distribution were used. No patient data was used and no statistical analysis was performed. Results: The D3D augmented reality system demonstrated stereoscopic depth perception by presenting a unique image to each eye, focal point convergence, head position tracking, 3D cursor, and joystick fly-through. Conclusion: The D3D augmented reality imaging system offers image viewing with depth perception and focal point convergence. The D3D augmented reality system should be tested to determine its utility in clinical practice. Keywords: augmented reality, 3D medical imaging, radiology, depth perception

3D Reconstruction of NMR Images

Directory of Open Access Journals (Sweden)

Peter Izak

2007-01-01

Full Text Available This paper introduces experiment of 3D reconstruction NMR images scanned from magnetic resonance device. There are described methods which can be used for 3D reconstruction magnetic resonance images in biomedical application. The main idea is based on marching cubes algorithm. For this task was chosen sophistication method by program Vision Assistant, which is a part of program LabVIEW.

Dependence of image quality on acquisition time for the PET/CT Biograph mCT

Energy Technology Data Exchange (ETDEWEB)

Molina-Duran, Flavia; Glatting, Gerhard [Heidelberg Univ., Mannheim (Germany). Medical Radiation Physics/Radiation Protection; Dinter, Dietmar; Schoenberg, Stefan O. [Heidelberg Univ., Mannheim (Germany). Inst. of Clinical Radiology and Nuclear Medicine; Schoenahl, Frederic [Siemens Healthcare Molecular Imaging, Erlangen (Germany)

2014-03-01

The impact of acquisition time on reconstructed PET image quality is analyzed for different acquisition times (1, 2, 3 and 4 min). Image quality was tested according to the National Electrical Manufacturers Association (NEMA) NU 2-2007, the evaluation for the signal to noise ratio (SNR) and the reconstructed activity ratio (RAR) for three algorithms, i.e. OSEM, TrueX and TOF applying different effective iteration numbers. The present work shows that the image quality of 3 and 4 min acquisition time for spherical lesions of 10 mm diameter are not significantly different between TrueX, TOF and OSEM. The 2 min acquisition time should be used carefully for the TrueX and OSEM algorithms in small lesions, because the levels of background noise are high compared to 3 or 4 min measurements. Also, the reconstructed activity ratio is underestimated to be approximately half of the expected value. For large lesions the three algorithms perform similarly for all acquisition durations, however, OSEM has the advantage of a more accurately reconstructed activity ratio compared to TrueX and TOF, which are more strongly influenced by noise. (orig.)

3D widefield light microscope image reconstruction without dyes

Science.gov (United States)

Larkin, S.; Larson, J.; Holmes, C.; Vaicik, M.; Turturro, M.; Jurkevich, A.; Sinha, S.; Ezashi, T.; Papavasiliou, G.; Brey, E.; Holmes, T.

2015-03-01

3D image reconstruction using light microscope modalities without exogenous contrast agents is proposed and investigated as an approach to produce 3D images of biological samples for live imaging applications. Multimodality and multispectral imaging, used in concert with this 3D optical sectioning approach is also proposed as a way to further produce contrast that could be specific to components in the sample. The methods avoid usage of contrast agents. Contrast agents, such as fluorescent or absorbing dyes, can be toxic to cells or alter cell behavior. Current modes of producing 3D image sets from a light microscope, such as 3D deconvolution algorithms and confocal microscopy generally require contrast agents. Zernike phase contrast (ZPC), transmitted light brightfield (TLB), darkfield microscopy and others can produce contrast without dyes. Some of these modalities have not previously benefitted from 3D image reconstruction algorithms, however. The 3D image reconstruction algorithm is based on an underlying physical model of scattering potential, expressed as the sample's 3D absorption and phase quantities. The algorithm is based upon optimizing an objective function - the I-divergence - while solving for the 3D absorption and phase quantities. Unlike typical deconvolution algorithms, each microscope modality, such as ZPC or TLB, produces two output image sets instead of one. Contrast in the displayed image and 3D renderings is further enabled by treating the multispectral/multimodal data as a feature set in a mathematical formulation that uses the principal component method of statistics.

Fully 3D PET image reconstruction using a fourier preconditioned conjugate-gradient algorithm

International Nuclear Information System (INIS)

Fessler, J.A.; Ficaro, E.P.

1996-01-01

Since the data sizes in fully 3D PET imaging are very large, iterative image reconstruction algorithms must converge in very few iterations to be useful. One can improve the convergence rate of the conjugate-gradient (CG) algorithm by incorporating preconditioning operators that approximate the inverse of the Hessian of the objective function. If the 3D cylindrical PET geometry were not truncated at the ends, then the Hessian of the penalized least-squares objective function would be approximately shift-invariant, i.e. G'G would be nearly block-circulant, where G is the system matrix. We propose a Fourier preconditioner based on this shift-invariant approximation to the Hessian. Results show that this preconditioner significantly accelerates the convergence of the CG algorithm with only a small increase in computation

Software for 3D diagnostic image reconstruction and analysis

International Nuclear Information System (INIS)

Taton, G.; Rokita, E.; Sierzega, M.; Klek, S.; Kulig, J.; Urbanik, A.

2005-01-01

Recent advances in computer technologies have opened new frontiers in medical diagnostics. Interesting possibilities are the use of three-dimensional (3D) imaging and the combination of images from different modalities. Software prepared in our laboratories devoted to 3D image reconstruction and analysis from computed tomography and ultrasonography is presented. In developing our software it was assumed that it should be applicable in standard medical practice, i.e. it should work effectively with a PC. An additional feature is the possibility of combining 3D images from different modalities. The reconstruction and data processing can be conducted using a standard PC, so low investment costs result in the introduction of advanced and useful diagnostic possibilities. The program was tested on a PC using DICOM data from computed tomography and TIFF files obtained from a 3D ultrasound system. The results of the anthropomorphic phantom and patient data were taken into consideration. A new approach was used to achieve spatial correlation of two independently obtained 3D images. The method relies on the use of four pairs of markers within the regions under consideration. The user selects the markers manually and the computer calculates the transformations necessary for coupling the images. The main software feature is the possibility of 3D image reconstruction from a series of two-dimensional (2D) images. The reconstructed 3D image can be: (1) viewed with the most popular methods of 3D image viewing, (2) filtered and processed to improve image quality, (3) analyzed quantitatively (geometrical measurements), and (4) coupled with another, independently acquired 3D image. The reconstructed and processed 3D image can be stored at every stage of image processing. The overall software performance was good considering the relatively low costs of the hardware used and the huge data sets processed. The program can be freely used and tested (source code and program available at

2D sparse array transducer optimization for 3D ultrasound imaging

International Nuclear Information System (INIS)

Choi, Jae Hoon; Park, Kwan Kyu

2014-01-01

A 3D ultrasound image is desired in many medical examinations. However, the implementation of a 2D array, which is needed for a 3D image, is challenging with respect to fabrication, interconnection and cabling. A 2D sparse array, which needs fewer elements than a dense array, is a realistic way to achieve 3D images. Because the number of ways the elements can be placed in an array is extremely large, a method for optimizing the array configuration is needed. Previous research placed the target point far from the transducer array, making it impossible to optimize the array in the operating range. In our study, we focused on optimizing a 2D sparse array transducer for 3D imaging by using a simulated annealing method. We compared the far-field optimization method with the near-field optimization method by analyzing a point-spread function (PSF). The resolution of the optimized sparse array is comparable to that of the dense array.

2D-Driven 3D Object Detection in RGB-D Images

KAUST Repository

Lahoud, Jean

2017-12-25

In this paper, we present a technique that places 3D bounding boxes around objects in an RGB-D scene. Our approach makes best use of the 2D information to quickly reduce the search space in 3D, benefiting from state-of-the-art 2D object detection techniques. We then use the 3D information to orient, place, and score bounding boxes around objects. We independently estimate the orientation for every object, using previous techniques that utilize normal information. Object locations and sizes in 3D are learned using a multilayer perceptron (MLP). In the final step, we refine our detections based on object class relations within a scene. When compared to state-of-the-art detection methods that operate almost entirely in the sparse 3D domain, extensive experiments on the well-known SUN RGB-D dataset [29] show that our proposed method is much faster (4.1s per image) in detecting 3D objects in RGB-D images and performs better (3 mAP higher) than the state-of-the-art method that is 4.7 times slower and comparably to the method that is two orders of magnitude slower. This work hints at the idea that 2D-driven object detection in 3D should be further explored, especially in cases where the 3D input is sparse.

Suitability of tile-based rendering for low-power 3d graphics accelerators

NARCIS (Netherlands)

Antochi, I.

2007-01-01

In this dissertation, we address low-power high performance 3D graphics accelerator architectures. The purpose of these accelerators is to relieve the burden of graphical computations from the main processor and also to achieve a better energy efficiency than can be achieved by executing these

Validity of computational hemodynamics in human arteries based on 3D time-of-flight MR angiography and 2D electrocardiogram gated phase contrast images

Science.gov (United States)

Yu, Huidan (Whitney); Chen, Xi; Chen, Rou; Wang, Zhiqiang; Lin, Chen; Kralik, Stephen; Zhao, Ye

2015-11-01

In this work, we demonstrate the validity of 4-D patient-specific computational hemodynamics (PSCH) based on 3-D time-of-flight (TOF) MR angiography (MRA) and 2-D electrocardiogram (ECG) gated phase contrast (PC) images. The mesoscale lattice Boltzmann method (LBM) is employed to segment morphological arterial geometry from TOF MRA, to extract velocity profiles from ECG PC images, and to simulate fluid dynamics on a unified GPU accelerated computational platform. Two healthy volunteers are recruited to participate in the study. For each volunteer, a 3-D high resolution TOF MRA image and 10 2-D ECG gated PC images are acquired to provide the morphological geometry and the time-varying flow velocity profiles for necessary inputs of the PSCH. Validation results will be presented through comparisons of LBM vs. 4D Flow Software for flow rates and LBM simulation vs. MRA measurement for blood flow velocity maps. Indiana University Health (IUH) Values Fund.

Denoising imaging polarimetry by adapted BM3D method.

Science.gov (United States)

Tibbs, Alexander B; Daly, Ilse M; Roberts, Nicholas W; Bull, David R

2018-04-01

In addition to the visual information contained in intensity and color, imaging polarimetry allows visual information to be extracted from the polarization of light. However, a major challenge of imaging polarimetry is image degradation due to noise. This paper investigates the mitigation of noise through denoising algorithms and compares existing denoising algorithms with a new method, based on BM3D (Block Matching 3D). This algorithm, Polarization-BM3D (PBM3D), gives visual quality superior to the state of the art across all images and noise standard deviations tested. We show that denoising polarization images using PBM3D allows the degree of polarization to be more accurately calculated by comparing it with spectral polarimetry measurements.

3D object-oriented image analysis in 3D geophysical modelling

DEFF Research Database (Denmark)

Fadel, I.; van der Meijde, M.; Kerle, N.

2015-01-01

Non-uniqueness of satellite gravity interpretation has traditionally been reduced by using a priori information from seismic tomography models. This reduction in the non-uniqueness has been based on velocity-density conversion formulas or user interpretation of the 3D subsurface structures (objects......) based on the seismic tomography models and then forward modelling these objects. However, this form of object-based approach has been done without a standardized methodology on how to extract the subsurface structures from the 3D models. In this research, a 3D object-oriented image analysis (3D OOA......) approach was implemented to extract the 3D subsurface structures from geophysical data. The approach was applied on a 3D shear wave seismic tomography model of the central part of the East African Rift System. Subsequently, the extracted 3D objects from the tomography model were reconstructed in the 3D...

3D Simulations for a Micron-Scale, Dielectric-Based Acceleration Experiment

International Nuclear Information System (INIS)

Yoder, R. B.; Travish, G.; Xu Jin; Rosenzweig, J. B.

2009-01-01

An experimental program to demonstrate a dielectric, slab-symmetric accelerator structure has been underway for the past two years. These resonant devices are driven by a side-coupled 800-nm laser and can be configured to maintain the field profile necessary for synchronous acceleration and focusing of relativistic or nonrelativistic particles. We present 3D simulations of various versions of the structure geometry, including a metal-walled structure relevant to ongoing cold tests on resonant properties, and an all-dielectric structure to be constructed for a proof-of-principle acceleration experiment.

Rainbow Particle Imaging Velocimetry for Dense 3D Fluid Velocity Imaging

KAUST Repository

Xiong, Jinhui

2017-04-11

Despite significant recent progress, dense, time-resolved imaging of complex, non-stationary 3D flow velocities remains an elusive goal. In this work we tackle this problem by extending an established 2D method, Particle Imaging Velocimetry, to three dimensions by encoding depth into color. The encoding is achieved by illuminating the flow volume with a continuum of light planes (a “rainbow”), such that each depth corresponds to a specific wavelength of light. A diffractive component in the camera optics ensures that all planes are in focus simultaneously. For reconstruction, we derive an image formation model for recovering stationary 3D particle positions. 3D velocity estimation is achieved with a variant of 3D optical flow that accounts for both physical constraints as well as the rainbow image formation model. We evaluate our method with both simulations and an experimental prototype setup.

3D reconstruction based on light field images

Science.gov (United States)

Zhu, Dong; Wu, Chunhong; Liu, Yunluo; Fu, Dongmei

2018-04-01

This paper proposed a method of reconstructing three-dimensional (3D) scene from two light field images capture by Lytro illium. The work was carried out by first extracting the sub-aperture images from light field images and using the scale-invariant feature transform (SIFT) for feature registration on the selected sub-aperture images. Structure from motion (SFM) algorithm is further used on the registration completed sub-aperture images to reconstruct the three-dimensional scene. 3D sparse point cloud was obtained in the end. The method shows that the 3D reconstruction can be implemented by only two light field camera captures, rather than at least a dozen times captures by traditional cameras. This can effectively solve the time-consuming, laborious issues for 3D reconstruction based on traditional digital cameras, to achieve a more rapid, convenient and accurate reconstruction.

Image based 3D city modeling : Comparative study

Directory of Open Access Journals (Sweden)

S. P. Singh

2014-06-01

Full Text Available 3D city model is a digital representation of the Earth’s surface and it’s related objects such as building, tree, vegetation, and some manmade feature belonging to urban area. The demand of 3D city modeling is increasing rapidly for various engineering and non-engineering applications. Generally four main image based approaches were used for virtual 3D city models generation. In first approach, researchers were used Sketch based modeling, second method is Procedural grammar based modeling, third approach is Close range photogrammetry based modeling and fourth approach is mainly based on Computer Vision techniques. SketchUp, CityEngine, Photomodeler and Agisoft Photoscan are the main softwares to represent these approaches respectively. These softwares have different approaches & methods suitable for image based 3D city modeling. Literature study shows that till date, there is no complete such type of comparative study available to create complete 3D city model by using images. This paper gives a comparative assessment of these four image based 3D modeling approaches. This comparative study is mainly based on data acquisition methods, data processing techniques and output 3D model products. For this research work, study area is the campus of civil engineering department, Indian Institute of Technology, Roorkee (India. This 3D campus acts as a prototype for city. This study also explains various governing parameters, factors and work experiences. This research work also gives a brief introduction, strengths and weakness of these four image based techniques. Some personal comment is also given as what can do or what can’t do from these softwares. At the last, this study shows; it concluded that, each and every software has some advantages and limitations. Choice of software depends on user requirements of 3D project. For normal visualization project, SketchUp software is a good option. For 3D documentation record, Photomodeler gives good

Whole-heart coronary MRA with 3D affine motion correction using 3D image-based navigation.

Science.gov (United States)

Henningsson, Markus; Prieto, Claudia; Chiribiri, Amedeo; Vaillant, Ghislain; Razavi, Reza; Botnar, René M

2014-01-01

Robust motion correction is necessary to minimize respiratory motion artefacts in coronary MR angiography (CMRA). The state-of-the-art method uses a 1D feet-head translational motion correction approach, and data acquisition is limited to a small window in the respiratory cycle, which prolongs the scan by a factor of 2-3. The purpose of this work was to implement 3D affine motion correction for Cartesian whole-heart CMRA using a 3D navigator (3D-NAV) to allow for data acquisition throughout the whole respiratory cycle. 3D affine transformations for different respiratory states (bins) were estimated by using 3D-NAV image acquisitions which were acquired during the startup profiles of a steady-state free precession sequence. The calculated 3D affine transformations were applied to the corresponding high-resolution Cartesian image acquisition which had been similarly binned, to correct for respiratory motion between bins. Quantitative and qualitative comparisons showed no statistical difference between images acquired with the proposed method and the reference method using a diaphragmatic navigator with a narrow gating window. We demonstrate that 3D-NAV and 3D affine correction can be used to acquire Cartesian whole-heart 3D coronary artery images with 100% scan efficiency with similar image quality as with the state-of-the-art gated and corrected method with approximately 50% scan efficiency. Copyright © 2013 Wiley Periodicals, Inc.

3D EIT image reconstruction with GREIT.

Science.gov (United States)

Grychtol, Bartłomiej; Müller, Beat; Adler, Andy

2016-06-01

Most applications of thoracic EIT use a single plane of electrodes on the chest from which a transverse image 'slice' is calculated. However, interpretation of EIT images is made difficult by the large region above and below the electrode plane to which EIT is sensitive. Volumetric EIT images using two (or more) electrode planes should help compensate, but are little used currently. The Graz consensus reconstruction algorithm for EIT (GREIT) has become popular in lung EIT. One shortcoming of the original formulation of GREIT is its restriction to reconstruction onto a 2D planar image. We present an extension of the GREIT algorithm to 3D and develop open-source tools to evaluate its performance as a function of the choice of stimulation and measurement pattern. Results show 3D GREIT using two electrode layers has significantly more uniform sensitivity profiles through the chest region. Overall, the advantages of 3D EIT are compelling.

Aerial 3D display by use of a 3D-shaped screen with aerial imaging by retro-reflection (AIRR)

Science.gov (United States)

Kurokawa, Nao; Ito, Shusei; Yamamoto, Hirotsugu

2017-06-01

The purpose of this paper is to realize an aerial 3D display. We design optical system that employs a projector below a retro-reflector and a 3D-shaped screen. A floating 3D image is formed with aerial imaging by retro-reflection (AIRR). Our proposed system is composed of a 3D-shaped screen, a projector, a quarter-wave retarder, a retro-reflector, and a reflective polarizer. Because AIRR forms aerial images that are plane-symmetric of the light sources regarding the reflective polarizer, the shape of the 3D screen is inverted from a desired aerial 3D image. In order to expand viewing angle, the 3D-shaped screen is surrounded by a retro-reflector. In order to separate the aerial image from reflected lights on the retro- reflector surface, the retro-reflector is tilted by 30 degrees. A projector is located below the retro-reflector at the same height of the 3D-shaped screen. The optical axis of the projector is orthogonal to the 3D-shaped screen. Scattered light on the 3D-shaped screen forms the aerial 3D image. In order to demonstrate the proposed optical design, a corner-cube-shaped screen is used for the 3D-shaped screen. Thus, the aerial 3D image is a cube that is floating above the reflective polarizer. For example, an aerial green cube is formed by projecting a calculated image on the 3D-shaped screen. The green cube image is digitally inverted in depth by our developed software. Thus, we have succeeded in forming aerial 3D image with our designed optical system.

Methodological development of topographic correction in 2D/3D ToF-SIMS images using AFM images

Science.gov (United States)

Jung, Seokwon; Lee, Nodo; Choi, Myungshin; Lee, Jungmin; Cho, Eunkyunng; Joo, Minho

2018-02-01

Time-of-flight secondary-ion mass spectrometry (ToF-SIMS) is an emerging technique that provides chemical information directly from the surface of electronic materials, e.g. OLED and solar cell. It is very versatile and highly sensitive mass spectrometric technique that provides surface molecular information with their lateral distribution as a two-dimensional (2D) molecular image. Extending the usefulness of ToF-SIMS, a 3D molecular image can be generated by acquiring multiple 2D images in a stack. These imaging techniques by ToF-SIMS provide an insight into understanding the complex structures of unknown composition in electronic material. However, one drawback in ToF-SIMS is not able to represent topographical information in 2D and 3D mapping images. To overcome this technical limitation, topographic information by ex-situ technique such as atomic force microscopy (AFM) has been combined with chemical information from SIMS that provides both chemical and physical information in one image. The key to combine two different images obtained from ToF-SIMS and AFM techniques is to develop the image processing algorithm, which performs resize and alignment by comparing the specific pixel information of each image. In this work, we present methodological development of the semiautomatic alignment and the 3D structure interpolation system for the combination of 2D/3D images obtained by ToF-SIMS and AFM measurements, which allows providing useful analytical information in a single representation.

2-D and 3-D computations of curved accelerator magnets

International Nuclear Information System (INIS)

Turner, L.R.

1991-01-01

In order to save computer memory, a long accelerator magnet may be computed by treating the long central region and the end regions separately. The dipole magnets for the injector synchrotron of the Advanced Photon Source (APS), now under construction at Argonne National Laboratory (ANL), employ magnet iron consisting of parallel laminations, stacked with a uniform radius of curvature of 33.379 m. Laplace's equation for the magnetic scalar potential has a different form for a straight magnet (x-y coordinates), a magnet with surfaces curved about a common center (r-θ coordinates), and a magnet with parallel laminations like the APS injector dipole. Yet pseudo 2-D computations for the three geometries give basically identical results, even for a much more strongly curved magnet. Hence 2-D (x-y) computations of the central region and 3-D computations of the end regions can be combined to determine the overall magnetic behavior of the magnets. 1 ref., 6 figs

Accelerated three-dimensional cine phase contrast imaging using randomly undersampled echo planar imaging with compressed sensing reconstruction.

Science.gov (United States)

Basha, Tamer A; Akçakaya, Mehmet; Goddu, Beth; Berg, Sophie; Nezafat, Reza

2015-01-01

The aim of this study was to implement and evaluate an accelerated three-dimensional (3D) cine phase contrast MRI sequence by combining a randomly sampled 3D k-space acquisition sequence with an echo planar imaging (EPI) readout. An accelerated 3D cine phase contrast MRI sequence was implemented by combining EPI readout with randomly undersampled 3D k-space data suitable for compressed sensing (CS) reconstruction. The undersampled data were then reconstructed using low-dimensional structural self-learning and thresholding (LOST). 3D phase contrast MRI was acquired in 11 healthy adults using an overall acceleration of 7 (EPI factor of 3 and CS rate of 3). For comparison, a single two-dimensional (2D) cine phase contrast scan was also performed with sensitivity encoding (SENSE) rate 2 and approximately at the level of the pulmonary artery bifurcation. The stroke volume and mean velocity in both the ascending and descending aorta were measured and compared between two sequences using Bland-Altman plots. An average scan time of 3 min and 30 s, corresponding to an acceleration rate of 7, was achieved for 3D cine phase contrast scan with one direction flow encoding, voxel size of 2 × 2 × 3 mm(3) , foot-head coverage of 6 cm and temporal resolution of 30 ms. The mean velocity and stroke volume in both the ascending and descending aorta were statistically equivalent between the proposed 3D sequence and the standard 2D cine phase contrast sequence. The combination of EPI with a randomly undersampled 3D k-space sampling sequence using LOST reconstruction allows a seven-fold reduction in scan time of 3D cine phase contrast MRI without compromising blood flow quantification. Copyright © 2014 John Wiley & Sons, Ltd.

Influence of high magnetic field strengths and parallel acquisition strategies on image quality in cardiac 2D CINE magnetic resonance imaging: comparison of 1.5 T vs. 3.0 T

International Nuclear Information System (INIS)

Gutberlet, Matthias; Schwinge, Kerstin; Freyhardt, Patrick; Spors, Birgit; Grothoff, Matthias; Denecke, Timm; Luedemann, Lutz; Felix, Roland; Noeske, Ralph; Niendorf, Thoralf

2005-01-01

The aim of this paper is to examine signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR) and image quality of cardiac CINE imaging at 1.5 T and 3.0 T. Twenty volunteers underwent cardiac magnetic resonance imaging (MRI) examinations using a 1.5-T and a 3.0-T scanner. Three different sets of breath-held, electrocardiogram-gated (ECG) CINE imaging techniques were employed, including: (1) unaccelerated SSFP (steady state free precession), (2) accelerated SSFP imaging and (3) gradient-echo-based myocardial tagging. Two-dimensional CINE SSFP at 3.0 T revealed an SNR improvement of 103% and a CNR increase of 19% as compared to the results obtained at 1.5 T. The SNR reduction in accelerated 2D CINE SSFP imaging was larger at 1.5 T (37%) compared to 3.0 T (26%). The mean SNR and CNR increase at 3.0 T obtained for the tagging sequence was 88% and 187%, respectively. At 3.0 T, the duration of the saturation bands persisted throughout the entire cardiac cycle. For comparison, the saturation bands were significantly diminished at 1.5 T during end-diastole. For 2D CINE SSFP imaging, no significant difference in the left ventricular volumetry and in the overall image quality was obtained. For myocardial tagging, image quality was significantly improved at 3.0 T. The SNR reduction in accelerated SSFP imaging was overcompensated by the increase in the baseline SNR at 3.0 T and did not result in any image quality degradation. For cardiac tagging techniques, 3.0 T was highly beneficial, which holds the promise to improve its diagnostic value. (orig.)

Filters in 2D and 3D Cardiac SPECT Image Processing

Directory of Open Access Journals (Sweden)

Maria Lyra

2014-01-01

Full Text Available Nuclear cardiac imaging is a noninvasive, sensitive method providing information on cardiac structure and physiology. Single photon emission tomography (SPECT evaluates myocardial perfusion, viability, and function and is widely used in clinical routine. The quality of the tomographic image is a key for accurate diagnosis. Image filtering, a mathematical processing, compensates for loss of detail in an image while reducing image noise, and it can improve the image resolution and limit the degradation of the image. SPECT images are then reconstructed, either by filter back projection (FBP analytical technique or iteratively, by algebraic methods. The aim of this study is to review filters in cardiac 2D, 3D, and 4D SPECT applications and how these affect the image quality mirroring the diagnostic accuracy of SPECT images. Several filters, including the Hanning, Butterworth, and Parzen filters, were evaluated in combination with the two reconstruction methods as well as with a specified MatLab program. Results showed that for both 3D and 4D cardiac SPECT the Butterworth filter, for different critical frequencies and orders, produced the best results. Between the two reconstruction methods, the iterative one might be more appropriate for cardiac SPECT, since it improves lesion detectability due to the significant improvement of image contrast.

3D surface reconstruction using optical flow for medical imaging

International Nuclear Information System (INIS)

Weng, Nan; Yang, Yee-Hong; Pierson, R.

1996-01-01

The recovery of a 3D model from a sequence of 2D images is very useful in medical image analysis. Image sequences obtained from the relative motion between the object and the camera or the scanner contain more 3D information than a single image. Methods to visualize the computed tomograms can be divided into two approaches: the surface rendering approach and the volume rendering approach. A new surface rendering method using optical flow is proposed. Optical flow is the apparent motion in the image plane produced by the projection of the real 3D motion onto 2D image. In this paper, the object remains stationary while the scanner undergoes translational motion. The 3D motion of an object can be recovered from the optical flow field using additional constraints. By extracting the surface information from 3D motion, it is possible to get an accurate 3D model of the object. Both synthetic and real image sequences have been used to illustrate the feasibility of the proposed method. The experimental results suggest that the proposed method is suitable for the reconstruction of 3D models from ultrasound medical images as well as other computed tomograms

An L1-norm phase constraint for half-Fourier compressed sensing in 3D MR imaging.

Science.gov (United States)

Li, Guobin; Hennig, Jürgen; Raithel, Esther; Büchert, Martin; Paul, Dominik; Korvink, Jan G; Zaitsev, Maxim

2015-10-01

In most half-Fourier imaging methods, explicit phase replacement is used. In combination with parallel imaging, or compressed sensing, half-Fourier reconstruction is usually performed in a separate step. The purpose of this paper is to report that integration of half-Fourier reconstruction into iterative reconstruction minimizes reconstruction errors. The L1-norm phase constraint for half-Fourier imaging proposed in this work is compared with the L2-norm variant of the same algorithm, with several typical half-Fourier reconstruction methods. Half-Fourier imaging with the proposed phase constraint can be seamlessly combined with parallel imaging and compressed sensing to achieve high acceleration factors. In simulations and in in-vivo experiments half-Fourier imaging with the proposed L1-norm phase constraint enables superior performance both reconstruction of image details and with regard to robustness against phase estimation errors. The performance and feasibility of half-Fourier imaging with the proposed L1-norm phase constraint is reported. Its seamless combination with parallel imaging and compressed sensing enables use of greater acceleration in 3D MR imaging.

Handbook of 3D machine vision optical metrology and imaging

CERN Document Server

Zhang, Song

2013-01-01

With the ongoing release of 3D movies and the emergence of 3D TVs, 3D imaging technologies have penetrated our daily lives. Yet choosing from the numerous 3D vision methods available can be frustrating for scientists and engineers, especially without a comprehensive resource to consult. Filling this gap, Handbook of 3D Machine Vision: Optical Metrology and Imaging gives an extensive, in-depth look at the most popular 3D imaging techniques. It focuses on noninvasive, noncontact optical methods (optical metrology and imaging). The handbook begins with the well-studied method of stereo vision and

High spatial resolution CT image reconstruction using parallel computing

International Nuclear Information System (INIS)

Yin Yin; Liu Li; Sun Gongxing

2003-01-01

Using the PC cluster system with 16 dual CPU nodes, we accelerate the FBP and OR-OSEM reconstruction of high spatial resolution image (2048 x 2048). Based on the number of projections, we rewrite the reconstruction algorithms into parallel format and dispatch the tasks to each CPU. By parallel computing, the speedup factor is roughly equal to the number of CPUs, which can be up to about 25 times when 25 CPUs used. This technique is very suitable for real-time high spatial resolution CT image reconstruction. (authors)

Wide area 2D/3D imaging development, analysis and applications

CERN Document Server

Langmann, Benjamin

2014-01-01

Imaging technology is an important research area and it is widely utilized in a growing number of disciplines ranging from gaming, robotics and automation to medicine. In the last decade 3D imaging became popular mainly driven by the introduction of novel 3D cameras and measuring devices. These cameras are usually limited to indoor scenes with relatively low distances. Benjamin Langmann introduces medium and long-range 2D/3D cameras to overcome these limitations. He reports measurement results for these devices and studies their characteristic behavior. In order to facilitate the application o

Accelerating 3D Elastic Wave Equations on Knights Landing based Intel Xeon Phi processors

Science.gov (United States)

Sourouri, Mohammed; Birger Raknes, Espen

2017-04-01

In advanced imaging methods like reverse-time migration (RTM) and full waveform inversion (FWI) the elastic wave equation (EWE) is numerically solved many times to create the seismic image or the elastic parameter model update. Thus, it is essential to optimize the solution time for solving the EWE as this will have a major impact on the total computational cost in running RTM or FWI. From a computational point of view applications implementing EWEs are associated with two major challenges. The first challenge is the amount of memory-bound computations involved, while the second challenge is the execution of such computations over very large datasets. So far, multi-core processors have not been able to tackle these two challenges, which eventually led to the adoption of accelerators such as Graphics Processing Units (GPUs). Compared to conventional CPUs, GPUs are densely populated with many floating-point units and fast memory, a type of architecture that has proven to map well to many scientific computations. Despite its architectural advantages, full-scale adoption of accelerators has yet to materialize. First, accelerators require a significant programming effort imposed by programming models such as CUDA or OpenCL. Second, accelerators come with a limited amount of memory, which also require explicit data transfers between the CPU and the accelerator over the slow PCI bus. The second generation of the Xeon Phi processor based on the Knights Landing (KNL) architecture, promises the computational capabilities of an accelerator but require the same programming effort as traditional multi-core processors. The high computational performance is realized through many integrated cores (number of cores and tiles and memory varies with the model) organized in tiles that are connected via a 2D mesh based interconnect. In contrary to accelerators, KNL is a self-hosted system, meaning explicit data transfers over the PCI bus are no longer required. However, like most

Preparing diagnostic 3D images for image registration with planning CT images

International Nuclear Information System (INIS)

Tracton, Gregg S.; Miller, Elizabeth P.; Rosenman, Julian; Chang, Sha X.; Sailer, Scott; Boxwala, Azaz; Chaney, Edward L.

1997-01-01

Purpose: Pre-radiotherapy (pre-RT) tomographic images acquired for diagnostic purposes often contain important tumor and/or normal tissue information which is poorly defined or absent in planning CT images. Our two years of clinical experience has shown that computer-assisted 3D registration of pre-RT images with planning CT images often plays an indispensable role in accurate treatment volume definition. Often the only available format of the diagnostic images is film from which the original 3D digital data must be reconstructed. In addition, any digital data, whether reconstructed or not, must be put into a form suitable for incorporation into the treatment planning system. The purpose of this investigation was to identify all problems that must be overcome before this data is suitable for clinical use. Materials and Methods: In the past two years we have 3D-reconstructed 300 diagnostic images from film and digital sources. As a problem was discovered we built a software tool to correct it. In time we collected a large set of such tools and found that they must be applied in a specific order to achieve the correct reconstruction. Finally, a toolkit (ediScan) was built that made all these tools available in the proper manner via a pleasant yet efficient mouse-based user interface. Results: Problems we discovered included different magnifications, shifted display centers, non-parallel image planes, image planes not perpendicular to the long axis of the table-top (shearing), irregularly spaced scans, non contiguous scan volumes, multiple slices per film, different orientations for slice axes (e.g. left-right reversal), slices printed at window settings corresponding to tissues of interest for diagnostic purposes, and printing artifacts. We have learned that the specific steps to correct these problems, in order of application, are: Also, we found that fast feedback and large image capacity (at least 2000 x 2000 12-bit pixels) are essential for practical application

A new tool fixation for external 3D head tracking using the Polaris Vicra system with the HRRT PET scanner

DEFF Research Database (Denmark)

Olesen, Oline Vinter; Andersen, Flemming; Holm, Søren

-water studies for up to 75 min (3-8 injections) were registered by the Polaris system in 4 volunteers. The tracking tool was fixed. Scans were divided into subframes based on the registered movements and reconstructed using the 3D-OSEM PSF method. The reconstructed subframes were repositioned to a reference......Objectives: The Polaris Vicra system (Northern Digital Inc.) is used for external 3D head registration with the Siemens HRRT brain PET. Our new tool fixation using a standard bandaid with a velcro-strap implies an improved frame repositioning. Methods: Head movements during serial PET 15O...... position and pairwise similarity of subframes was evaluated before and after the repositioning. Results: Registered movements during scans were less than 4.3mm with. Images were compared before/after motion correction. Conclusions: Our new velcro band-aid fixation is suitable for clinical use: easy to use...

Segmentation of the lumen and media-adventitia boundaries of the common carotid artery from 3D ultrasound images

Science.gov (United States)

Ukwatta, E.; Awad, J.; Ward, A. D.; Samarabandu, J.; Krasinski, A.; Parraga, G.; Fenster, A.

2011-03-01

Three-dimensional ultrasound (3D US) vessel wall volume (VWV) measurements provide high measurement sensitivity and reproducibility for the monitoring and assessment of carotid atherosclerosis. In this paper, we describe a semiautomated approach based on the level set method to delineate the media-adventitia and lumen boundaries of the common carotid artery from 3D US images to support the computation of VWV. Due to the presence of plaque and US image artifacts, the carotid arteries are challenging to segment using image information alone. Our segmentation framework combines several image cues with domain knowledge and limited user interaction. Our method was evaluated with respect to manually outlined boundaries on 430 2D US images extracted from 3D US images of 30 patients who have carotid stenosis of 60% or more. The VWV given by our method differed from that given by manual segmentation by 6.7% +/- 5.0%. For the media-adventitia and lumen segmentations, respectively, our method yielded Dice coefficients of 95.2% +/- 1.6%, 94.3% +/- 2.6%, mean absolute distances of 0.3 +/- 0.1 mm, 0.2 +/- 0.1 mm, maximum absolute distances of 0.8 +/- 0.4 mm, 0.6 +/- 0.3 mm, and volume differences of 4.2% +/- 3.1%, 3.4% +/- 2.6%. The realization of a semi-automated segmentation method will accelerate the translation of 3D carotid US to clinical care for the rapid, non-invasive, and economical monitoring of atherosclerotic disease progression and regression during therapy.

Objective and subjective comparison of standard 2-D and fully 3-D reconstructed data on a PET/CT system.

Science.gov (United States)

Strobel, Klaus; Rüdy, Matthias; Treyer, Valerie; Veit-Haibach, Patrick; Burger, Cyrill; Hany, Thomas F

2007-07-01

The relative advantage of fully 3-D versus 2-D mode for whole-body imaging is currently the focus of considerable expert debate. The nature of 3-D PET acquisition for FDG PET/CT theoretically allows a shorter scan time and improved efficiency of FDG use than in the standard 2-D acquisition. We therefore objectively and subjectively compared standard 2-D and fully 3-D reconstructed data for FDG PET/CT on a research PET/CT system. In a total of 36 patients (mean 58.9 years, range 17.3-78.9 years; 21 male, 15 female) referred for known or suspected malignancy, FDG PET/CT was performed using a research PET/CT system with advanced detector technology with improved sensitivity and spatial resolution. After 45 min uptake, a low-dose CT (40 mAs) from head to thigh was performed followed by 2-D PET (emission 3 min per field) and 3-D PET (emission 1.5 min per field) with both seven slices overlap to cover the identical anatomical region. Acquisition time was therefore 50% less (seven fields; 21 min vs. 10.5 min). PET data was acquired in a randomized fashion, so in 50% of the cases 2-D data was acquired first. CT data was used for attenuation correction. 2-D (OSEM) and 3-D PET images were iteratively reconstructed. Subjective analysis of 2-D and 3-D images was performed by two readers in a blinded, randomized fashion evaluating the following criteria: sharpness of organs (liver, chest wall/lung), overall image quality and detectability and dignity of each identified lesion. Objective analysis of PET data was investigated measuring maximum standard uptake value with lean body mass (SUV(max,LBM)) of identified lesions. On average, per patient, the SUV(max) was 7.86 (SD 7.79) for 2-D and 6.96 (SD 5.19) for 3-D. On a lesion basis, the average SUV(max) was 7.65 (SD 7.79) for 2-D and 6.75 (SD 5.89) for 3-D. The absolute difference on a paired t-test of SUV 3-D-2-D based on each measured lesion was significant with an average of -0.956 (P=0.002) and an average of -0.884 on a

GPU-based, parallel-line, omni-directional integration of measured acceleration field to obtain the 3D pressure distribution

Science.gov (United States)

Wang, Jin; Zhang, Cao; Katz, Joseph

2016-11-01

A PIV based method to reconstruct the volumetric pressure field by direct integration of the 3D material acceleration directions has been developed. Extending the 2D virtual-boundary omni-directional method (Omni2D, Liu & Katz, 2013), the new 3D parallel-line omni-directional method (Omni3D) integrates the material acceleration along parallel lines aligned in multiple directions. Their angles are set by a spherical virtual grid. The integration is parallelized on a Tesla K40c GPU, which reduced the computing time from three hours to one minute for a single realization. To validate its performance, this method is utilized to calculate the 3D pressure fields in isotropic turbulence and channel flow using the JHU DNS Databases (http://turbulence.pha.jhu.edu). Both integration of the DNS acceleration as well as acceleration from synthetic 3D particles are tested. Results are compared to other method, e.g. solution to the Pressure Poisson Equation (e.g. PPE, Ghaemi et al., 2012) with Bernoulli based Dirichlet boundary conditions, and the Omni2D method. The error in Omni3D prediction is uniformly low, and its sensitivity to acceleration errors is local. It agrees with the PPE/Bernoulli prediction away from the Dirichlet boundary. The Omni3D method is also applied to experimental data obtained using tomographic PIV, and results are correlated with deformation of a compliant wall. ONR.

Effects of intra-operative fluoroscopic 3D-imaging on peri-operative imaging strategy in calcaneal fracture surgery.

Science.gov (United States)

Beerekamp, M S H; Backes, M; Schep, N W L; Ubbink, D T; Luitse, J S; Schepers, T; Goslings, J C

2017-12-01

Previous studies demonstrated that intra-operative fluoroscopic 3D-imaging (3D-imaging) in calcaneal fracture surgery is promising to prevent revision surgery and save costs. However, these studies limited their focus to corrections performed after 3D-imaging, thereby neglecting corrections after intra-operative fluoroscopic 2D-imaging (2D-imaging). The aim of this study was to assess the effects of additional 3D-imaging on intra-operative corrections, peri-operative imaging used, and patient-relevant outcomes compared to 2D-imaging alone. In this before-after study, data of adult patients who underwent open reduction and internal fixation (ORIF) of a calcaneal fracture between 2000 and 2014 in our level-I Trauma center were collected. 3D-imaging (BV Pulsera with 3D-RX, Philips Healthcare, Best, The Netherlands) was available as of 2007 at the surgeons' discretion. Patient and fracture characteristics, peri-operative imaging, intra-operative corrections and patient-relevant outcomes were collected from the hospital databases. Patients in whom additional 3D-imaging was applied were compared to those undergoing 2D-imaging alone. A total of 231 patients were included of whom 107 (46%) were operated with the use of 3D-imaging. No significant differences were found in baseline characteristics. The median duration of surgery was significantly longer when using 3D-imaging (2:08 vs. 1:54 h; p = 0.002). Corrections after additional 3D-imaging were performed in 53% of the patients. However, significantly fewer corrections were made after 2D-imaging when 3D-imaging was available (Risk difference (RD) -15%; 95% Confidence interval (CI) -29 to -2). Peri-operative imaging, besides intra-operative 3D-imaging, and patient-relevant outcomes were similar between groups. Intra-operative 3D-imaging provides additional information resulting in additional corrections. Moreover, 3D-imaging probably changed the surgeons' attitude to rely more on 3D-imaging, hence a 15%-decrease of

2D vs. 3D imaging in laparoscopic surgery-results of a prospective randomized trial.

Science.gov (United States)

Buia, Alexander; Stockhausen, Florian; Filmann, Natalie; Hanisch, Ernst

2017-12-01

3D imaging is an upcoming technology in laparoscopic surgery, and recent studies have shown that the modern 3D technique is superior in an experimental setting. However, the first randomized controlled clinical trial in this context dates back to 1998 and showed no significant difference between 2D and 3D visualization using the first 3D generation technique, which is now more than 15 years old. Positive results measured in an experimental setting considering 3D imaging on surgical performance led us to initiate a randomized controlled pragmatic clinical trial to validate our findings in daily clinical routine. Standard laparoscopic operations (cholecystectomy, appendectomy) were preoperatively randomized to a 2D or 3D imaging system. We used a surgical comfort scale (Likert scale) and the Raw NASA Workload TLX for the subjective assessment of 2D and 3D imaging; the duration of surgery was also measured. The results of 3D imaging were statistically significant better than 2D imaging concerning the parameters "own felt safety" and "task efficiency"; the difficulty level of the procedures in the 2D and 3D groups did not differ. Overall, the Raw NASA Workload TLX showed no significance between the groups. 3D imaging could be a possible advantage in laparoscopic surgery. The results of our clinical trial show increased personal felt safety and efficiency of the surgeon using a 3D imaging system. Overall of the procedures, the findings assessed using Likert scales in terms of own felt safety and task efficiency were statistically significant for 3D imaging. The individually perceived workload assessed with the Raw NASA TLX shows no difference. Although these findings are subjective impressions of the performing surgeons without a clear benefit for 3D technology in clinical outcome, we think that these results show the capability that 3D laparoscopy can have a positive impact while performing laparoscopic procedures.

Effect of Post-Reconstruction Gaussian Filtering on Image Quality and Myocardial Blood Flow Measurement with N-13 Ammonia PET

International Nuclear Information System (INIS)

Kim, Hyeon Sik; Cho, Sang-Geon; Kim, Ju Han; Kwon, Seong Young; Lee, Byeong-il; Bom, Hee-Seung

2014-01-01

In order to evaluate the effect of post-reconstruction Gaussian filtering on image quality and myocardial blood flow (MBF) measurement by dynamic N-13 ammonia positron emission tomography (PET), we compared various reconstruction and filtering methods with image characteristics. Dynamic PET images of three patients with coronary artery disease (male-female ratio of 2:1; age: 57, 53, and 76 years) were reconstructed, using filtered back projection (FBP) and ordered subset expectation maximization (OSEM) methods. OSEM reconstruction consisted of OSEM-2I, OSEM-4I, and OSEM-6I with 2, 4, and 6 iterations, respectively. The images, reconstructed and filtered by Gaussian filters of 5, 10, and 15 mm, were obtained, as well as non-filtered images. Visual analysis of image quality (IQ) was performed using a 3-grade scoring system by 2 independent readers, blinded to the reconstruction and filtering methods of stress images. Then, signal-to-noise ratio (SNR) was calculated by noise and contrast recovery (CR). Stress and rest MBF and coronary flow reserve (CFR) were obtained for each method. IQ scores, stress and rest MBF, and CFR were compared between the methods, using Chi-square and Kruskal-Wallis tests. In the visual analysis, IQ was significantly higher by 10 mm Gaussian filtering, compared to other sizes of filter (P<0.001 for both readers). However, no significant difference of IQ was found between FBP and various numbers of iteration in OSEM (P=0.923 and 0.855 for readers 1 and 2, respectively). SNR was significantly higher in 10 mm Gaussian filter. There was a significant difference in stress and rest MBF between several vascular territories. However CFR was not significantly different according to various filtering methods. Post-reconstruction Gaussian filtering with a filter size of 10 mm significantly enhances the IQ of N-13 ammonia PET-CT, without changing the results of CFR calculation

Advancement in PET quantification using 3D-OP-OSEM point spread function reconstruction with the HRRT

Energy Technology Data Exchange (ETDEWEB)

Varrone, Andrea; Sjoeholm, Nils; Gulyas, Balazs; Halldin, Christer; Farde, Lars [Karolinska Hospital, Karolinska Institutet, Department of Clinical Neuroscience, Psychiatry Section and Stockholm Brain Institute, Stockholm (Sweden); Eriksson, Lars [Karolinska Hospital, Karolinska Institutet, Department of Clinical Neuroscience, Psychiatry Section and Stockholm Brain Institute, Stockholm (Sweden); Siemens Molecular Imaging, Knoxville, TN (United States); University of Stockholm, Department of Physics, Stockholm (Sweden)

2009-10-15

Image reconstruction including the modelling of the point spread function (PSF) is an approach improving the resolution of the PET images. This study assessed the quantitative improvements provided by the implementation of the PSF modelling in the reconstruction of the PET data using the High Resolution Research Tomograph (HRRT). Measurements were performed on the NEMA-IEC/2001 (Image Quality) phantom for image quality and on an anthropomorphic brain phantom (STEPBRAIN). PSF reconstruction was also applied to PET measurements in two cynomolgus monkeys examined with [{sup 18}F]FE-PE2I (dopamine transporter) and with [{sup 11}C]MNPA (D{sub 2} receptor), and in one human subject examined with [{sup 11}C]raclopride (D{sub 2} receptor). PSF reconstruction increased the recovery coefficient (RC) in the NEMA phantom by 11-40% and the grey to white matter ratio in the STEPBRAIN phantom by 17%. PSF reconstruction increased binding potential (BP{sub ND}) in the striatum and midbrain by 14 and 18% in the [{sup 18}F]FE-PE2I study, and striatal BP{sub ND} by 6 and 10% in the [{sup 11}C]MNPA and [{sup 11}C]raclopride studies. PSF reconstruction improved quantification by increasing the RC and thus reducing the partial volume effect. This method provides improved conditions for PET quantification in clinical studies with the HRRT system, particularly when targeting receptor populations in small brain structures. (orig.)

3D and 4D magnetic susceptibility tomography based on complex MR images

Science.gov (United States)

Chen, Zikuan; Calhoun, Vince D

2014-11-11

Magnetic susceptibility is the physical property for T2*-weighted magnetic resonance imaging (T2*MRI). The invention relates to methods for reconstructing an internal distribution (3D map) of magnetic susceptibility values, .chi. (x,y,z), of an object, from 3D T2*MRI phase images, by using Computed Inverse Magnetic Resonance Imaging (CIMRI) tomography. The CIMRI technique solves the inverse problem of the 3D convolution by executing a 3D Total Variation (TV) regularized iterative convolution scheme, using a split Bregman iteration algorithm. The reconstruction of .chi. (x,y,z) can be designed for low-pass, band-pass, and high-pass features by using a convolution kernel that is modified from the standard dipole kernel. Multiple reconstructions can be implemented in parallel, and averaging the reconstructions can suppress noise. 4D dynamic magnetic susceptibility tomography can be implemented by reconstructing a 3D susceptibility volume from a 3D phase volume by performing 3D CIMRI magnetic susceptibility tomography at each snapshot time.

PET/CT detectability and classification of simulated pulmonary lesions using an SUV correction scheme

Science.gov (United States)

Morrow, Andrew N.; Matthews, Kenneth L., II; Bujenovic, Steven

2008-03-01

Positron emission tomography (PET) and computed tomography (CT) together are a powerful diagnostic tool, but imperfect image quality allows false positive and false negative diagnoses to be made by any observer despite experience and training. This work investigates PET acquisition mode, reconstruction method and a standard uptake value (SUV) correction scheme on the classification of lesions as benign or malignant in PET/CT images, in an anthropomorphic phantom. The scheme accounts for partial volume effect (PVE) and PET resolution. The observer draws a region of interest (ROI) around the lesion using the CT dataset. A simulated homogenous PET lesion of the same shape as the drawn ROI is blurred with the point spread function (PSF) of the PET scanner to estimate the PVE, providing a scaling factor to produce a corrected SUV. Computer simulations showed that the accuracy of the corrected PET values depends on variations in the CT-drawn boundary and the position of the lesion with respect to the PET image matrix, especially for smaller lesions. Correction accuracy was affected slightly by mismatch of the simulation PSF and the actual scanner PSF. The receiver operating characteristic (ROC) study resulted in several observations. Using observer drawn ROIs, scaled tumor-background ratios (TBRs) more accurately represented actual TBRs than unscaled TBRs. For the PET images, 3D OSEM outperformed 2D OSEM, 3D OSEM outperformed 3D FBP, and 2D OSEM outperformed 2D FBP. The correction scheme significantly increased sensitivity and slightly increased accuracy for all acquisition and reconstruction modes at the cost of a small decrease in specificity.

3D Printing Openable Imaging Phantom Design

International Nuclear Information System (INIS)

Kim, Myoung Keun; Won, Jun Hyeok; Lee, Seung Wook

2017-01-01

The purpose of this study is to design an openable phantom that can replace the internal measurement bar used for contrast comparison in order to increase the efficiency of manufacturing imaging phantom used in the medical industry and to improve convenience using 3D printer. Phantom concept design, 3D printing, and Image reconstruction were defined as the scope of the thesis. Also, we study metal artifact reduction with openable phantom. We have designed a Openable phantom using 3D printing, and have investigated metal artifact reduction after inserting a metallic material inside the phantom. The openable phantom can be adjusted at any time to suit the user's experiment and can be easily replaced and useful.

R and D of Nb(3)Sn accelerator magnets at Fermilab

International Nuclear Information System (INIS)

Zlobin, A.V.; Ambrosio, G.; Andreev, N.; Barzi, E; Bordini, B.; Bossert, R.; Carcagno, R.; Chichili, D.R.; DiMarco, J.; Elementi, L.; Feher, S.; Kashikhin, V.S.; Kashikhin, V.V.; Kephart, R.; Lamm, M.; Limon, P.J.; Novitski, I.; Orris, D.; Pischalnikov, Yu.; Schlabach, P.; Stanek, R.

2004-01-01

Fermilab is developing and investigating different high-field magnet designs for present and future accelerators. The magnet RandD program was focused on the 10-12 T accelerator magnets based on Nb 3 Sn superconductor and explored both basic magnet technologies for brittle superconductors--wind-and-react and react-and-wind. Magnet design studies in support of LHC upgrades and VLHC are being performed. A series of 1-m long single-bore models of cos-theta Nb 3 Sn dipoles based on wind-and-react technique was fabricated and tested. Three 1-m long flat racetracks and the common coil dipole model, based on a single-layer coil and wide reacted Nb 3 Sn cable, have also been fabricated and tested. Extensive theoretical studies of magnetic instabilities in Nb 3 Sn strands, cable and magnet were performed which led to successful 10 T dipole model. This paper presents the details of the Fermilab's high field accelerator magnet program, reports its status and major results, and formulates the program next steps

R and D of Nb(3)Sn accelerator magnets at Fermilab

Energy Technology Data Exchange (ETDEWEB)

Zlobin, A.V.; Ambrosio, G.; Andreev, N.; Barzi, E; Bordini, B.; Bossert, R.; Carcagno, R.; Chichili, D.R.; DiMarco, J.; Elementi, L.; Feher, S.; Kashikhin, V.S.; Kashikhin, V.V.; Kephart, R.; Lamm, M.; Limon, P.J.; Novitski, I.; Orris, D.; Pischalnikov, Yu.; Schlabach, P.; Stanek, R.; /Fermilab

2004-11-01

Fermilab is developing and investigating different high-field magnet designs for present and future accelerators. The magnet R&D program was focused on the 10-12 T accelerator magnets based on Nb{sub 3}Sn superconductor and explored both basic magnet technologies for brittle superconductors--wind-and-react and react-and-wind. Magnet design studies in support of LHC upgrades and VLHC are being performed. A series of 1-m long single-bore models of cos-theta Nb{sub 3}Sn dipoles based on wind-and-react technique was fabricated and tested. Three 1-m long flat racetracks and the common coil dipole model, based on a single-layer coil and wide reacted Nb{sub 3}Sn cable, have also been fabricated and tested. Extensive theoretical studies of magnetic instabilities in Nb{sub 3}Sn strands, cable and magnet were performed which led to successful 10 T dipole model. This paper presents the details of the Fermilab's high field accelerator magnet program, reports its status and major results, and formulates the program next steps.

Optical 3D watermark based digital image watermarking for telemedicine

Science.gov (United States)

Li, Xiao Wei; Kim, Seok Tae

2013-12-01

Region of interest (ROI) of a medical image is an area including important diagnostic information and must be stored without any distortion. This algorithm for application of watermarking technique for non-ROI of the medical image preserving ROI. The paper presents a 3D watermark based medical image watermarking scheme. In this paper, a 3D watermark object is first decomposed into 2D elemental image array (EIA) by a lenslet array, and then the 2D elemental image array data is embedded into the host image. The watermark extraction process is an inverse process of embedding. The extracted EIA through the computational integral imaging reconstruction (CIIR) technique, the 3D watermark can be reconstructed. Because the EIA is composed of a number of elemental images possesses their own perspectives of a 3D watermark object. Even though the embedded watermark data badly damaged, the 3D virtual watermark can be successfully reconstructed. Furthermore, using CAT with various rule number parameters, it is possible to get many channels for embedding. So our method can recover the weak point having only one transform plane in traditional watermarking methods. The effectiveness of the proposed watermarking scheme is demonstrated with the aid of experimental results.

From medical imaging data to 3D printed anatomical models.

Directory of Open Access Journals (Sweden)

Thore M Bücking

Full Text Available Anatomical models are important training and teaching tools in the clinical environment and are routinely used in medical imaging research. Advances in segmentation algorithms and increased availability of three-dimensional (3D printers have made it possible to create cost-efficient patient-specific models without expert knowledge. We introduce a general workflow that can be used to convert volumetric medical imaging data (as generated by Computer Tomography (CT to 3D printed physical models. This process is broken up into three steps: image segmentation, mesh refinement and 3D printing. To lower the barrier to entry and provide the best options when aiming to 3D print an anatomical model from medical images, we provide an overview of relevant free and open-source image segmentation tools as well as 3D printing technologies. We demonstrate the utility of this streamlined workflow by creating models of ribs, liver, and lung using a Fused Deposition Modelling 3D printer.

Development of 3-D Medical Image VIsualization System

African Journals Online (AJOL)

User

uses standard 2-D medical imaging inputs and generates medical images of human body parts ... light wave from points on the 3-D object(s) in ... tools, and communication bandwidth cannot .... locations along the track that correspond with.

3D ultrasound imaging for prosthesis fabrication and diagnostic imaging

Energy Technology Data Exchange (ETDEWEB)

Morimoto, A.K.; Bow, W.J.; Strong, D.S. [and others

1995-06-01

The fabrication of a prosthetic socket for a below-the-knee amputee requires knowledge of the underlying bone structure in order to provide pressure relief for sensitive areas and support for load bearing areas. The goal is to enable the residual limb to bear pressure with greater ease and utility. Conventional methods of prosthesis fabrication are based on limited knowledge about the patient`s underlying bone structure. A 3D ultrasound imaging system was developed at Sandia National Laboratories. The imaging system provides information about the location of the bones in the residual limb along with the shape of the skin surface. Computer assisted design (CAD) software can use this data to design prosthetic sockets for amputees. Ultrasound was selected as the imaging modality. A computer model was developed to analyze the effect of the various scanning parameters and to assist in the design of the overall system. The 3D ultrasound imaging system combines off-the-shelf technology for image capturing, custom hardware, and control and image processing software to generate two types of image data -- volumetric and planar. Both volumetric and planar images reveal definition of skin and bone geometry with planar images providing details on muscle fascial planes, muscle/fat interfaces, and blood vessel definition. The 3D ultrasound imaging system was tested on 9 unilateral below-the- knee amputees. Image data was acquired from both the sound limb and the residual limb. The imaging system was operated in both volumetric and planar formats. An x-ray CT (Computed Tomography) scan was performed on each amputee for comparison. Results of the test indicate beneficial use of ultrasound to generate databases for fabrication of prostheses at a lower cost and with better initial fit as compared to manually fabricated prostheses.

Computer assisted determination of acetabular cup orientation using 2D-3D image registration

International Nuclear Information System (INIS)

Zheng, Guoyan; Zhang, Xuan

2010-01-01

2D-3D image-based registration methods have been developed to measure acetabular cup orientation after total hip arthroplasty (THA). These methods require registration of both the prosthesis and the CT images to 2D radiographs and compute implant position with respect to a reference. The application of these methods is limited in clinical practice due to two limitations: (1) the requirement of a computer-aided design (CAD) model of the prosthesis, which may be unavailable due to the proprietary concerns of the manufacturer, and (2) the requirement of either multiple radiographs or radiograph-specific calibration, usually unavailable for retrospective studies. In this paper, we propose a new method to address these limitations. A new formulation for determination of post-operative cup orientation, which couples a radiographic measurement with 2D-3D image matching, was developed. In our formulation, the radiographic measurement can be obtained with known methods so that the challenge lies in the 2D-3D image matching. To solve this problem, a hybrid 2D-3D registration scheme combining a landmark-to-ray 2D-3D alignment with a robust intensity-based 2D-3D registration was used. The hybrid 2D-3D registration scheme allows computing both the post-operative cup orientation with respect to an anatomical reference and the pelvic tilt and rotation with respect to the X-ray imaging table/plate. The method was validated using 2D adult cadaver hips. Using the hybrid 2D-3D registration scheme, our method showed a mean accuracy of 1.0 ± 0.7 (range from 0.1 to 2.0 ) for inclination and 1.7 ± 1.2 (range from 0.0 to 3.9 ) for anteversion, taking the measurements from post-operative CT images as ground truths. Our new solution formulation and the hybrid 2D-3D registration scheme facilitate estimation of post-operative cup orientation and measurement of pelvic tilt and rotation. (orig.)

3D Seismic Imaging using Marchenko Methods

Science.gov (United States)

Lomas, A.; Curtis, A.

2017-12-01

Marchenko methods are novel, data driven techniques that allow seismic wavefields from sources and receivers on the Earth's surface to be redatumed to construct wavefields with sources in the subsurface - including complex multiply-reflected waves, and without the need for a complex reference model. In turn, this allows subsurface images to be constructed at any such subsurface redatuming points (image or virtual receiver points). Such images are then free of artefacts from multiply-scattered waves that usually contaminate migrated seismic images. Marchenko algorithms require as input the same information as standard migration methods: the full reflection response from sources and receivers at the Earth's surface, and an estimate of the first arriving wave between the chosen image point and the surface. The latter can be calculated using a smooth velocity model estimated using standard methods. The algorithm iteratively calculates a signal that focuses at the image point to create a virtual source at that point, and this can be used to retrieve the signal between the virtual source and the surface. A feature of these methods is that the retrieved signals are naturally decomposed into up- and down-going components. That is, we obtain both the signal that initially propagated upwards from the virtual source and arrived at the surface, separated from the signal that initially propagated downwards. Figure (a) shows a 3D subsurface model with a variable density but a constant velocity (3000m/s). Along the surface of this model (z=0) in both the x and y directions are co-located sources and receivers at 20-meter intervals. The redatumed signal in figure (b) has been calculated using Marchenko methods from a virtual source (1200m, 500m and 400m) to the surface. For comparison the true solution is given in figure (c), and shows a good match when compared to figure (b). While these 2D redatuming and imaging methods are still in their infancy having first been developed in

Efficient Sample Delay Calculation for 2-D and 3-D Ultrasound Imaging.

Science.gov (United States)

Ibrahim, Aya; Hager, Pascal A; Bartolini, Andrea; Angiolini, Federico; Arditi, Marcel; Thiran, Jean-Philippe; Benini, Luca; De Micheli, Giovanni

2017-08-01

Ultrasound imaging is a reference medical diagnostic technique, thanks to its blend of versatility, effectiveness, and moderate cost. The core computation of all ultrasound imaging methods is based on simple formulae, except for those required to calculate acoustic propagation delays with high precision and throughput. Unfortunately, advanced three-dimensional (3-D) systems require the calculation or storage of billions of such delay values per frame, which is a challenge. In 2-D systems, this requirement can be four orders of magnitude lower, but efficient computation is still crucial in view of low-power implementations that can be battery-operated, enabling usage in numerous additional scenarios. In this paper, we explore two smart designs of the delay generation function. To quantify their hardware cost, we implement them on FPGA and study their footprint and performance. We evaluate how these architectures scale to different ultrasound applications, from a low-power 2-D system to a next-generation 3-D machine. When using numerical approximations, we demonstrate the ability to generate delay values with sufficient throughput to support 10 000-channel 3-D imaging at up to 30 fps while using 63% of a Virtex 7 FPGA, requiring 24 MB of external memory accessed at about 32 GB/s bandwidth. Alternatively, with similar FPGA occupation, we show an exact calculation method that reaches 24 fps on 1225-channel 3-D imaging and does not require external memory at all. Both designs can be scaled to use a negligible amount of resources for 2-D imaging in low-power applications and for ultrafast 2-D imaging at hundreds of frames per second.

3D imaging and characterisation of strengthening particles in inconel 718 using FIB tomography

Energy Technology Data Exchange (ETDEWEB)

Kruk, Adam; Gruszczynski, Adam; Czyrska-Filemonowicz, Aleksandra [AGH University of Science and Technology, Faculty of Metals Engineering and Industrial Computer Science, Al. A. Mickiewicza 30, 30-059 Krakow (Poland)

2011-07-01

The Inconel 718 is a commercial nickel-base superalloy, widely used for critical pieces in turbine engines. Its microstructure consists of the {gamma} matrix and strengthening coherent nanoparticles {gamma}' and {gamma}''. In the present work FIB tomography technique was used for imaging and characterisation of strengthening particles. FIB tomography is based on a serial sectioning procedure using a FIB/SEM dual beam workstation. Repeated removal of layers as thin as several nm for some hundred times allows to investigate at total a volume of some {mu}m3 with a voxel size as 2.5 nm x 2.5 nm x 2.5 nm. 3D mapping of nanoparticles with high Z-resolution by serial FIB slicing (in a distance of about 2.5 nm) and SEM imaging was performed. Ga ion beam at 30 kV was used to perform a precise in-situ milling. The SEM images at accelerating voltage 1.5 kV were taken with using ESB detector. The real 3D-data of precipitates obtained by FIB tomography, open a new possibility for microstructure analysis of materials for industrial applications.

Hands-on guide for 3D image creation for geological purposes

Science.gov (United States)

Frehner, Marcel; Tisato, Nicola

2013-04-01

Geological structures in outcrops or hand specimens are inherently three dimensional (3D), and therefore better understandable if viewed in 3D. While 3D models can easily be created, manipulated, and looked at from all sides on the computer screen (e.g., using photogrammetry or laser scanning data), 3D visualizations for publications or conference posters are much more challenging as they have to live in a 2D-world (i.e., on a sheet of paper). Perspective 2D visualizations of 3D models do not fully transmit the "feeling and depth of the third dimension" to the audience; but this feeling is desirable for a better examination and understanding in 3D of the structure under consideration. One of the very few possibilities to generate real 3D images, which work on a 2D display, is by using so-called stereoscopic images. Stereoscopic images are two images of the same object recorded from two slightly offset viewpoints. Special glasses and techniques have to be used to make sure that one image is seen only by one eye, and the other image is seen by the other eye, which together lead to the "3D effect". Geoscientists are often familiar with such 3D images. For example, geomorphologists traditionally view stereographic orthophotos by employing a mirror-steroscope. Nowadays, petroleum-geoscientists examine high-resolution 3D seismic data sets in special 3D visualization rooms. One of the methods for generating and viewing a stereoscopic image, which does not require a high-tech viewing device, is to create a so-called anaglyph. The principle is to overlay two images saturated in red and cyan, respectively. The two images are then viewed through red-cyan-stereoscopic glasses. This method is simple and cost-effective, but has some drawbacks in preserving colors accurately. A similar method is used in 3D movies, where polarized light or shuttering techniques are used to separate the left from the right image, which allows preserving the original colors. The advantage of red

3D printing of gas jet nozzles for laser-plasma accelerators

Energy Technology Data Exchange (ETDEWEB)

Döpp, A.; Guillaume, E.; Thaury, C.; Gautier, J.; Ta Phuoc, K.; Malka, V. [LOA, ENSTA ParisTech, CNRS, École Polytechnique, Université Paris-Saclay, 828 Boulevard des Maréchaux, 91762 Palaiseau Cedex (France)

2016-07-15

Recent results on laser wakefield acceleration in tailored plasma channels have underlined the importance of controlling the density profile of the gas target. In particular, it was reported that the appropriate density tailoring can result in improved injection, acceleration, and collimation of laser-accelerated electron beams. To achieve such profiles, innovative target designs are required. For this purpose, we have reviewed the usage of additive layer manufacturing, commonly known as 3D printing, in order to produce gas jet nozzles. Notably we have compared the performance of two industry standard techniques, namely, selective laser sintering (SLS) and stereolithography (SLA). Furthermore we have used the common fused deposition modeling to reproduce basic gas jet designs and used SLA and SLS for more sophisticated nozzle designs. The nozzles are characterized interferometrically and used for electron acceleration experiments with the SALLE JAUNE terawatt laser at Laboratoire d’Optique Appliquée.

Highly accelerated cardiovascular MR imaging using many channel technology: concepts and clinical applications

International Nuclear Information System (INIS)

Niendorf, Thoralf; Sodickson, Daniel K.

2008-01-01

Cardiovascular magnetic resonance imaging (CVMRI) is of proven clinical value in the non-invasive imaging of cardiovascular diseases. CVMRI requires rapid image acquisition, but acquisition speed is fundamentally limited in conventional MRI. Parallel imaging provides a means for increasing acquisition speed and efficiency. However, signal-to-noise (SNR) limitations and the limited number of receiver channels available on most MR systems have in the past imposed practical constraints, which dictated the use of moderate accelerations in CVMRI. High levels of acceleration, which were unattainable previously, have become possible with many-receiver MR systems and many-element, cardiac-optimized RF-coil arrays. The resulting imaging speed improvements can be exploited in a number of ways, ranging from enhancement of spatial and temporal resolution to efficient whole heart coverage to streamlining of CVMRI work flow. In this review, examples of these strategies are provided, following an outline of the fundamentals of the highly accelerated imaging approaches employed in CVMRI. Topics discussed include basic principles of parallel imaging; key requirements for MR systems and RF-coil design; practical considerations of SNR management, supported by multi-dimensional accelerations, 3D noise averaging and high field imaging; highly accelerated clinical state-of-the art cardiovascular imaging applications spanning the range from SNR-rich to SNR-limited; and current trends and future directions. (orig.)

Motion robust high resolution 3D free-breathing pulmonary MRI using dynamic 3D image self-navigator.

Science.gov (United States)

Jiang, Wenwen; Ong, Frank; Johnson, Kevin M; Nagle, Scott K; Hope, Thomas A; Lustig, Michael; Larson, Peder E Z

2018-06-01

To achieve motion robust high resolution 3D free-breathing pulmonary MRI utilizing a novel dynamic 3D image navigator derived directly from imaging data. Five-minute free-breathing scans were acquired with a 3D ultrashort echo time (UTE) sequence with 1.25 mm isotropic resolution. From this data, dynamic 3D self-navigating images were reconstructed under locally low rank (LLR) constraints and used for motion compensation with one of two methods: a soft-gating technique to penalize the respiratory motion induced data inconsistency, and a respiratory motion-resolved technique to provide images of all respiratory motion states. Respiratory motion estimation derived from the proposed dynamic 3D self-navigator of 7.5 mm isotropic reconstruction resolution and a temporal resolution of 300 ms was successful for estimating complex respiratory motion patterns. This estimation improved image quality compared to respiratory belt and DC-based navigators. Respiratory motion compensation with soft-gating and respiratory motion-resolved techniques provided good image quality from highly undersampled data in volunteers and clinical patients. An optimized 3D UTE sequence combined with the proposed reconstruction methods can provide high-resolution motion robust pulmonary MRI. Feasibility was shown in patients who had irregular breathing patterns in which our approach could depict clinically relevant pulmonary pathologies. Magn Reson Med 79:2954-2967, 2018. © 2017 International Society for Magnetic Resonance in Medicine. © 2017 International Society for Magnetic Resonance in Medicine.

Accelerated gradient methods for total-variation-based CT image reconstruction

Energy Technology Data Exchange (ETDEWEB)

Joergensen, Jakob H.; Hansen, Per Christian [Technical Univ. of Denmark, Lyngby (Denmark). Dept. of Informatics and Mathematical Modeling; Jensen, Tobias L.; Jensen, Soeren H. [Aalborg Univ. (Denmark). Dept. of Electronic Systems; Sidky, Emil Y.; Pan, Xiaochuan [Chicago Univ., Chicago, IL (United States). Dept. of Radiology

2011-07-01

Total-variation (TV)-based CT image reconstruction has shown experimentally to be capable of producing accurate reconstructions from sparse-view data. In particular TV-based reconstruction is well suited for images with piecewise nearly constant regions. Computationally, however, TV-based reconstruction is demanding, especially for 3D imaging, and the reconstruction from clinical data sets is far from being close to real-time. This is undesirable from a clinical perspective, and thus there is an incentive to accelerate the solution of the underlying optimization problem. The TV reconstruction can in principle be found by any optimization method, but in practice the large scale of the systems arising in CT image reconstruction preclude the use of memory-intensive methods such as Newton's method. The simple gradient method has much lower memory requirements, but exhibits prohibitively slow convergence. In the present work we address the question of how to reduce the number of gradient method iterations needed to achieve a high-accuracy TV reconstruction. We consider the use of two accelerated gradient-based methods, GPBB and UPN, to solve the 3D-TV minimization problem in CT image reconstruction. The former incorporates several heuristics from the optimization literature such as Barzilai-Borwein (BB) step size selection and nonmonotone line search. The latter uses a cleverly chosen sequence of auxiliary points to achieve a better convergence rate. The methods are memory efficient and equipped with a stopping criterion to ensure that the TV reconstruction has indeed been found. An implementation of the methods (in C with interface to Matlab) is available for download from http://www2.imm.dtu.dk/~pch/TVReg/. We compare the proposed methods with the standard gradient method, applied to a 3D test problem with synthetic few-view data. We find experimentally that for realistic parameters the proposed methods significantly outperform the standard gradient method. (orig.)

MR imaging in epilepsy with use of 3D MP-RAGE

International Nuclear Information System (INIS)

Tanaka, Akio; Ohno, Sigeru; Sei, Tetsuro; Kanazawa, Susumu; Yasui, Koutaro; Kuroda, Masahiro; Hiraki, Yoshio; Oka, Eiji

1996-01-01

The patients were 40 males and 33 females; their ages ranged from 1 month to 39 years (mean: 15.7 years). The patients underwent MR imaging, including spin-echo T 1 -weighted, turbo spin-echo proton density/T 2 -weighted, and 3D magnetization-prepared rapid gradient-echo (3D MP-RAGE) images. These examinations disclosed 39 focal abnormalities. On visual evaluation, the boundary of abnormal gray matter in the neuronal migration disorder (NMD) cases was most clealy shown on 3D MP-RAGE images as compared to the other images. This is considered to be due to the higher spatial resolution and the better contrast of the 3D MP-RAGE images than those of the other techniques. The relative contrast difference between abnormal gray matter and the adjacent white matter was also assessed. The results revealed that the contrast differences on the 3D MP-RAGE images were larger than those on the other images; this was statistically significant. Although the sensitivity of 3D MP-RAGE for NMD was not specifically evaluated in this study, the possibility of this disorder, in cases suspected on other images, could be ruled out. Thus, it appears that the specificity with respect to NMD was at least increased with us of 3D MP-RAGE. 3D MP-RAGE also enabled us to build three-dimensional surface models that were helpful in understanding the three-dimensional anatomy. Furthermore. 3D MP-RAGE was considered to be the best technique for evaluating hippocampus atrophy in patients with MTS. On the other hand, the sensitivity in the signal change of the hippocampus was higher on T 2 -weighted images. In addition, demonstration of cortical tubers of tuberous sclerosis in neurocutaneous syndrome was superior on T 2 -weighted images than on 3D MP-RAGE images. (K.H.)

Accuracy of volume measurement using 3D ultrasound and development of CT-3D US image fusion algorithm for prostate cancer radiotherapy

International Nuclear Information System (INIS)

Baek, Jihye; Huh, Jangyoung; Hyun An, So; Oh, Yoonjin; Kim, Myungsoo; Kim, DongYoung; Chung, Kwangzoo; Cho, Sungho; Lee, Rena

2013-01-01

Purpose: To evaluate the accuracy of measuring volumes using three-dimensional ultrasound (3D US), and to verify the feasibility of the replacement of CT-MR fusion images with CT-3D US in radiotherapy treatment planning. Methods: Phantoms, consisting of water, contrast agent, and agarose, were manufactured. The volume was measured using 3D US, CT, and MR devices. A CT-3D US and MR-3D US image fusion software was developed using the Insight Toolkit library in order to acquire three-dimensional fusion images. The quality of the image fusion was evaluated using metric value and fusion images. Results: Volume measurement, using 3D US, shows a 2.8 ± 1.5% error, 4.4 ± 3.0% error for CT, and 3.1 ± 2.0% error for MR. The results imply that volume measurement using the 3D US devices has a similar accuracy level to that of CT and MR. Three-dimensional image fusion of CT-3D US and MR-3D US was successfully performed using phantom images. Moreover, MR-3D US image fusion was performed using human bladder images. Conclusions: 3D US could be used in the volume measurement of human bladders and prostates. CT-3D US image fusion could be used in monitoring the target position in each fraction of external beam radiation therapy. Moreover, the feasibility of replacing the CT-MR image fusion to the CT-3D US in radiotherapy treatment planning was verified.

Influence of image reconstruction methods on statistical parametric mapping of brain PET images

International Nuclear Information System (INIS)

Yin Dayi; Chen Yingmao; Yao Shulin; Shao Mingzhe; Yin Ling; Tian Jiahe; Cui Hongyan

2007-01-01

Objective: Statistic parametric mapping (SPM) was widely recognized as an useful tool in brain function study. The aim of this study was to investigate if imaging reconstruction algorithm of PET images could influence SPM of brain. Methods: PET imaging of whole brain was performed in six normal volunteers. Each volunteer had two scans with true and false acupuncturing. The PET scans were reconstructed using ordered subsets expectation maximization (OSEM) and filtered back projection (FBP) with 3 varied parameters respectively. The images were realigned, normalized and smoothed using SPM program. The difference between true and false acupuncture scans was tested using a matched pair t test at every voxel. Results: (1) SPM corrected multiple comparison (P corrected uncorrected <0.001): SPM derived from the images with different reconstruction method were different. The largest difference, in number and position of the activated voxels, was noticed between FBP and OSEM re- construction algorithm. Conclusions: The method of PET image reconstruction could influence the results of SPM uncorrected multiple comparison. Attention should be paid when the conclusion was drawn using SPM uncorrected multiple comparison. (authors)

3D Reconstruction of NMR Images by LabVIEW

Directory of Open Access Journals (Sweden)

Peter IZAK

2007-01-01

Full Text Available This paper introduces the experiment of 3D reconstruction NMR images via virtual instrumentation - LabVIEW. The main idea is based on marching cubes algorithm and image processing implemented by module of Vision assistant. The two dimensional images shot by the magnetic resonance device provide information about the surface properties of human body. There is implemented algorithm which can be used for 3D reconstruction of magnetic resonance images in biomedical application.

PHOTOGRAMMETRIC 3D BUILDING RECONSTRUCTION FROM THERMAL IMAGES

Directory of Open Access Journals (Sweden)

E. Maset

2017-08-01

Full Text Available This paper addresses the problem of 3D building reconstruction from thermal infrared (TIR images. We show that a commercial Computer Vision software can be used to automatically orient sequences of TIR images taken from an Unmanned Aerial Vehicle (UAV and to generate 3D point clouds, without requiring any GNSS/INS data about position and attitude of the images nor camera calibration parameters. Moreover, we propose a procedure based on Iterative Closest Point (ICP algorithm to create a model that combines high resolution and geometric accuracy of RGB images with the thermal information deriving from TIR images. The process can be carried out entirely by the aforesaid software in a simple and efficient way.

Iterative reconstruction or filtered backprojection for semi-quantitative assessment of dopamine D2 receptor SPECT studies?

International Nuclear Information System (INIS)

Koch, Walter; Suessmair, Christine; Tatsch, Klaus; Poepperl, Gabriele

2011-01-01

In routine clinical practice striatal dopamine D 2 receptor binding is generally assessed using data reconstructed by filtered backprojection (FBP). The aim of this study was to investigate the use of an iterative reconstruction algorithm (ordered subset expectation maximization, OSEM) and to assess whether it may provide comparable or even better results than those obtained by standard FBP. In 56 patients with parkinsonian syndromes, single photon emission computed tomography (SPECT) scans were acquired 2 h after i.v. application of 185 MBq [ 123 I]iodobenzamide (IBZM) using a triple-head gamma camera (Siemens MS 3). The scans were reconstructed both by FBP and OSEM (3 iterations, 8 subsets) and filtered using a Butterworth filter. After attenuation correction the studies were automatically fitted to a mean template with a corresponding 3-D volume of interest (VOI) map covering striatum (S), caudate (C), putamen (P) and several reference VOIs using BRASS software. Visual assessment of the fitted studies suggests a better separation between C and P in studies reconstructed by OSEM than FBP. Unspecific background activity appears more homogeneous after iterative reconstruction. The correlation shows a good accordance of dopamine receptor binding using FBP and OSEM (intra-class correlation coefficients S: 0.87; C: 0.88; P: 0.84). Receiver-operating characteristic (ROC) analyses show comparable diagnostic power of OSEM and FBP in the differentiation between idiopathic parkinsonian syndrome (IPS) and non-IPS. Iterative reconstruction of IBZM SPECT studies for assessment of the D 2 receptors is feasible in routine clinical practice. Close correlations between FBP and OSEM data suggest that iteratively reconstructed IBZM studies allow reliable quantification of dopamine receptor binding even though a gain in diagnostic power could not be demonstrated. (orig.)

3D Interpolation Method for CT Images of the Lung

Directory of Open Access Journals (Sweden)

Noriaki Asada

2003-06-01

Full Text Available A 3-D image can be reconstructed from numerous CT images of the lung. The procedure reconstructs a solid from multiple cross section images, which are collected during pulsation of the heart. Thus the motion of the heart is a special factor that must be taken into consideration during reconstruction. The lung exhibits a repeating transformation synchronized to the beating of the heart as an elastic body. There are discontinuities among neighboring CT images due to the beating of the heart, if no special techniques are used in taking CT images. The 3-D heart image is reconstructed from numerous CT images in which both the heart and the lung are taken. Although the outline shape of the reconstructed 3-D heart is quite unnatural, the envelope of the 3-D unnatural heart is fit to the shape of the standard heart. The envelopes of the lung in the CT images are calculated after the section images of the best fitting standard heart are located at the same positions of the CT images. Thus the CT images are geometrically transformed to the optimal CT images fitting best to the standard heart. Since correct transformation of images is required, an Area oriented interpolation method proposed by us is used for interpolation of transformed images. An attempt to reconstruct a 3-D lung image by a series of such operations without discontinuity is shown. Additionally, the same geometrical transformation method to the original projection images is proposed as a more advanced method.

Accelerated two-dimensional cine DENSE cardiovascular magnetic resonance using compressed sensing and parallel imaging.

Science.gov (United States)

Chen, Xiao; Yang, Yang; Cai, Xiaoying; Auger, Daniel A; Meyer, Craig H; Salerno, Michael; Epstein, Frederick H

2016-06-14

Cine Displacement Encoding with Stimulated Echoes (DENSE) provides accurate quantitative imaging of cardiac mechanics with rapid displacement and strain analysis; however, image acquisition times are relatively long. Compressed sensing (CS) with parallel imaging (PI) can generally provide high-quality images recovered from data sampled below the Nyquist rate. The purposes of the present study were to develop CS-PI-accelerated acquisition and reconstruction methods for cine DENSE, to assess their accuracy for cardiac imaging using retrospective undersampling, and to demonstrate their feasibility for prospectively-accelerated 2D cine DENSE imaging in a single breathhold. An accelerated cine DENSE sequence with variable-density spiral k-space sampling and golden angle rotations through time was implemented. A CS method, Block LOw-rank Sparsity with Motion-guidance (BLOSM), was combined with sensitivity encoding (SENSE) for the reconstruction of under-sampled multi-coil spiral data. Seven healthy volunteers and 7 patients underwent 2D cine DENSE imaging with fully-sampled acquisitions (14-26 heartbeats in duration) and with prospectively rate-2 and rate-4 accelerated acquisitions (14 and 8 heartbeats in duration). Retrospectively- and prospectively-accelerated data were reconstructed using BLOSM-SENSE and SENSE. Image quality of retrospectively-undersampled data was quantified using the relative root mean square error (rRMSE). Myocardial displacement and circumferential strain were computed for functional assessment, and linear correlation and Bland-Altman analyses were used to compare accelerated acquisitions to fully-sampled reference datasets. For retrospectively-undersampled data, BLOSM-SENSE provided similar or lower rRMSE at rate-2 and lower rRMSE at rate-4 acceleration compared to SENSE (p cine DENSE provided good image quality and expected values of displacement and strain. BLOSM-SENSE-accelerated spiral cine DENSE imaging with 2D displacement encoding can be

3D IMAGING USING COHERENT SYNCHROTRON RADIATION

Directory of Open Access Journals (Sweden)

Peter Cloetens

2011-05-01

Full Text Available Three dimensional imaging is becoming a standard tool for medical, scientific and industrial applications. The use of modem synchrotron radiation sources for monochromatic beam micro-tomography provides several new features. Along with enhanced signal-to-noise ratio and improved spatial resolution, these include the possibility of quantitative measurements, the easy incorporation of special sample environment devices for in-situ experiments, and a simple implementation of phase imaging. These 3D approaches overcome some of the limitations of 2D measurements. They require new tools for image analysis.

AUTOMATED CELL SEGMENTATION WITH 3D FLUORESCENCE MICROSCOPY IMAGES.

Science.gov (United States)

Kong, Jun; Wang, Fusheng; Teodoro, George; Liang, Yanhui; Zhu, Yangyang; Tucker-Burden, Carol; Brat, Daniel J

2015-04-01

A large number of cell-oriented cancer investigations require an effective and reliable cell segmentation method on three dimensional (3D) fluorescence microscopic images for quantitative analysis of cell biological properties. In this paper, we present a fully automated cell segmentation method that can detect cells from 3D fluorescence microscopic images. Enlightened by fluorescence imaging techniques, we regulated the image gradient field by gradient vector flow (GVF) with interpolated and smoothed data volume, and grouped voxels based on gradient modes identified by tracking GVF field. Adaptive thresholding was then applied to voxels associated with the same gradient mode where voxel intensities were enhanced by a multiscale cell filter. We applied the method to a large volume of 3D fluorescence imaging data of human brain tumor cells with (1) small cell false detection and missing rates for individual cells; and (2) trivial over and under segmentation incidences for clustered cells. Additionally, the concordance of cell morphometry structure between automated and manual segmentation was encouraging. These results suggest a promising 3D cell segmentation method applicable to cancer studies.

Comparison of post-contrast 3D-T1-MPRAGE, 3D-T1-SPACE and 3D-T2-FLAIR MR images in evaluation of meningeal abnormalities at 3-T MRI.

Science.gov (United States)

Jeevanandham, Balaji; Kalyanpur, Tejas; Gupta, Prashant; Cherian, Mathew

2017-06-01

This study was to assess the usefulness of newer three-dimensional (3D)-T 1 sampling perfection with application optimized contrast using different flip-angle evolutions (SPACE) and 3D-T 2 fluid-attenuated inversion recovery (FLAIR) sequences in evaluation of meningeal abnormalities. 78 patients who presented with high suspicion of meningeal abnormalities were evaluated using post-contrast 3D-T 2 -FLAIR, 3D-T 1 magnetization-prepared rapid gradient-echo (MPRAGE) and 3D-T 1 -SPACE sequences. The images were evaluated independently by two radiologists for cortical gyral, sulcal space, basal cisterns and dural enhancement. The diagnoses were confirmed by further investigations including histopathology. Post-contrast 3D-T 1 -SPACE and 3D-T 2 -FLAIR images yielded significantly more information than MPRAGE images (p evaluation of meningeal abnormalities and when used in combination have the maximum sensitivity for leptomeningeal abnormalities. The negative-predictive value is nearly 100%, where no leptomeningeal abnormality was detected on these sequences. Advances in knowledge: Post-contrast 3D-T 1 -SPACE and 3D-T 2 -FLAIR images are more useful than 3D-T 1 -MPRAGE images in evaluation of meningeal abnormalities.

Potential Cost Savings with 3D Printing Combined With 3D Imaging and CPLM for Fleet Maintenance and Revitalization

Science.gov (United States)

2014-05-01

1 Potential Cost Savings with 3D Printing Combined With 3D Imaging and CPLM for Fleet Maintenance and Revitalization David N. Ford...2014 4. TITLE AND SUBTITLE Potential Cost Savings with 3D Printing Combined With 3D Imaging and CPLM for Fleet Maintenance and Revitalization 5a...Manufacturing ( 3D printing ) 2 Research Context Problem: Learning curve savings forecasted in SHIPMAIN maintenance initiative have not materialized

3D quantitative phase imaging of neural networks using WDT

Science.gov (United States)

Kim, Taewoo; Liu, S. C.; Iyer, Raj; Gillette, Martha U.; Popescu, Gabriel

2015-03-01

White-light diffraction tomography (WDT) is a recently developed 3D imaging technique based on a quantitative phase imaging system called spatial light interference microscopy (SLIM). The technique has achieved a sub-micron resolution in all three directions with high sensitivity granted by the low-coherence of a white-light source. Demonstrations of the technique on single cell imaging have been presented previously; however, imaging on any larger sample, including a cluster of cells, has not been demonstrated using the technique. Neurons in an animal body form a highly complex and spatially organized 3D structure, which can be characterized by neuronal networks or circuits. Currently, the most common method of studying the 3D structure of neuron networks is by using a confocal fluorescence microscope, which requires fluorescence tagging with either transient membrane dyes or after fixation of the cells. Therefore, studies on neurons are often limited to samples that are chemically treated and/or dead. WDT presents a solution for imaging live neuron networks with a high spatial and temporal resolution, because it is a 3D imaging method that is label-free and non-invasive. Using this method, a mouse or rat hippocampal neuron culture and a mouse dorsal root ganglion (DRG) neuron culture have been imaged in order to see the extension of processes between the cells in 3D. Furthermore, the tomogram is compared with a confocal fluorescence image in order to investigate the 3D structure at synapses.

Progressive attenuation fields: Fast 2D-3D image registration without precomputation

International Nuclear Information System (INIS)

Rohlfing, Torsten; Russakoff, Daniel B.; Denzler, Joachim; Mori, Kensaku; Maurer, Calvin R. Jr.

2005-01-01

Computation of digitally reconstructed radiograph (DRR) images is the rate-limiting step in most current intensity-based algorithms for the registration of three-dimensional (3D) images to two-dimensional (2D) projection images. This paper introduces and evaluates the progressive attenuation field (PAF), which is a new method to speed up DRR computation. A PAF is closely related to an attenuation field (AF). A major difference is that a PAF is constructed on the fly as the registration proceeds; it does not require any precomputation time, nor does it make any prior assumptions of the patient pose or limit the permissible range of patient motion. A PAF effectively acts as a cache memory for projection values once they are computed, rather than as a lookup table for precomputed projections like standard AFs. We use a cylindrical attenuation field parametrization, which is better suited for many medical applications of 2D-3D registration than the usual two-plane parametrization. The computed attenuation values are stored in a hash table for time-efficient storage and access. Using clinical gold-standard spine image data sets from five patients, we demonstrate consistent speedups of intensity-based 2D-3D image registration using PAF DRRs by a factor of 10 over conventional ray casting DRRs with no decrease of registration accuracy or robustness

Extracting 3D Parametric Curves from 2D Images of Helical Objects.

Science.gov (United States)

Willcocks, Chris G; Jackson, Philip T G; Nelson, Carl J; Obara, Boguslaw

2017-09-01

Helical objects occur in medicine, biology, cosmetics, nanotechnology, and engineering. Extracting a 3D parametric curve from a 2D image of a helical object has many practical applications, in particular being able to extract metrics such as tortuosity, frequency, and pitch. We present a method that is able to straighten the image object and derive a robust 3D helical curve from peaks in the object boundary. The algorithm has a small number of stable parameters that require little tuning, and the curve is validated against both synthetic and real-world data. The results show that the extracted 3D curve comes within close Hausdorff distance to the ground truth, and has near identical tortuosity for helical objects with a circular profile. Parameter insensitivity and robustness against high levels of image noise are demonstrated thoroughly and quantitatively.

Convolution of large 3D images on GPU and its decomposition

Science.gov (United States)

Karas, Pavel; Svoboda, David

2011-12-01

In this article, we propose a method for computing convolution of large 3D images. The convolution is performed in a frequency domain using a convolution theorem. The algorithm is accelerated on a graphic card by means of the CUDA parallel computing model. Convolution is decomposed in a frequency domain using the decimation in frequency algorithm. We pay attention to keeping our approach efficient in terms of both time and memory consumption and also in terms of memory transfers between CPU and GPU which have a significant inuence on overall computational time. We also study the implementation on multiple GPUs and compare the results between the multi-GPU and multi-CPU implementations.

Extracting 3D parametric curves from 2D images of helical objects.

OpenAIRE

Willcocks, Chris; Jackson, Philip T.G.; Nelson, Carl J.; Obara, Boguslaw

2016-01-01

Helical objects occur in medicine, biology, cosmetics, nanotechnology, and engineering. Extracting a 3D parametric curve from a 2D image of a helical object has many practical applications, in particular being able to extract metrics such as tortuosity, frequency, and pitch. We present a method that is able to straighten the image object and derive a robust 3D helical curve from peaks in the object boundary. The algorithm has a small number of stable parameters that require little tuning, and...

Simultaneous Multislice Echo Planar Imaging With Blipped Controlled Aliasing in Parallel Imaging Results in Higher Acceleration: A Promising Technique for Accelerated Diffusion Tensor Imaging of Skeletal Muscle.

Science.gov (United States)

Filli, Lukas; Piccirelli, Marco; Kenkel, David; Guggenberger, Roman; Andreisek, Gustav; Beck, Thomas; Runge, Val M; Boss, Andreas

2015-07-01

The aim of this study was to investigate the feasibility of accelerated diffusion tensor imaging (DTI) of skeletal muscle using echo planar imaging (EPI) applying simultaneous multislice excitation with a blipped controlled aliasing in parallel imaging results in higher acceleration unaliasing technique. After federal ethics board approval, the lower leg muscles of 8 healthy volunteers (mean [SD] age, 29.4 [2.9] years) were examined in a clinical 3-T magnetic resonance scanner using a 15-channel knee coil. The EPI was performed at a b value of 500 s/mm2 without slice acceleration (conventional DTI) as well as with 2-fold and 3-fold acceleration. Fractional anisotropy (FA) and mean diffusivity (MD) were measured in all 3 acquisitions. Fiber tracking performance was compared between the acquisitions regarding the number of tracks, average track length, and anatomical precision using multivariate analysis of variance and Mann-Whitney U tests. Acquisition time was 7:24 minutes for conventional DTI, 3:53 minutes for 2-fold acceleration, and 2:38 minutes for 3-fold acceleration. Overall FA and MD values ranged from 0.220 to 0.378 and 1.595 to 1.829 mm2/s, respectively. Two-fold acceleration yielded similar FA and MD values (P ≥ 0.901) and similar fiber tracking performance compared with conventional DTI. Three-fold acceleration resulted in comparable MD (P = 0.199) but higher FA values (P = 0.006) and significantly impaired fiber tracking in the soleus and tibialis anterior muscles (number of tracks, P DTI of skeletal muscle with similar image quality and quantification accuracy of diffusion parameters. This may increase the clinical applicability of muscle anisotropy measurements.

Developing 3D Imaging Programmes-Workflow and Quality Control

OpenAIRE

Hess, M.; Robson, S.; Serpico, M.; Amati, G.; Pridden, I.; Nelson, T.

2016-01-01

This article reports on a successful project for 3D imaging research, digital applications, and use of new technologies in the museum. The article will focus on the development and implementation of a viable workflow for the production of high-quality 3D models of museum objects, based on the 3D laser scanning and photogrammetry of selected ancient Egyptian artefacts. The development of a robust protocol for the complete process chain for imaging cultural heritage artefacts, from the acquisit...

3D ultrasound imaging : Fast and cost-effective morphometry of musculoskeletal tissue

NARCIS (Netherlands)

Weide, Guido; Van Der Zwaard, Stephan; Huijing, Peter A.; Jaspers, Richard T.; Harlaar, Jaap

2017-01-01

The developmental goal of 3D ultrasound imaging (3DUS) is to engineer a modality to perform 3D morphological ultrasound analysis of human muscles. 3DUS images are constructed from calibrated freehand 2D B-mode ultrasound images, which are positioned into a voxel array. Ultrasound (US) imaging allows

Flatbed-type 3D display systems using integral imaging method

Science.gov (United States)

Hirayama, Yuzo; Nagatani, Hiroyuki; Saishu, Tatsuo; Fukushima, Rieko; Taira, Kazuki

2006-10-01

We have developed prototypes of flatbed-type autostereoscopic display systems using one-dimensional integral imaging method. The integral imaging system reproduces light beams similar of those produced by a real object. Our display architecture is suitable for flatbed configurations because it has a large margin for viewing distance and angle and has continuous motion parallax. We have applied our technology to 15.4-inch displays. We realized horizontal resolution of 480 with 12 parallaxes due to adoption of mosaic pixel arrangement of the display panel. It allows viewers to see high quality autostereoscopic images. Viewing the display from angle allows the viewer to experience 3-D images that stand out several centimeters from the surface of the display. Mixed reality of virtual 3-D objects and real objects are also realized on a flatbed display. In seeking reproduction of natural 3-D images on the flatbed display, we developed proprietary software. The fast playback of the CG movie contents and real-time interaction are realized with the aid of a graphics card. Realization of the safety 3-D images to the human beings is very important. Therefore, we have measured the effects on the visual function and evaluated the biological effects. For example, the accommodation and convergence were measured at the same time. The various biological effects are also measured before and after the task of watching 3-D images. We have found that our displays show better results than those to a conventional stereoscopic display. The new technology opens up new areas of application for 3-D displays, including arcade games, e-learning, simulations of buildings and landscapes, and even 3-D menus in restaurants.

3D Imaging with Structured Illumination for Advanced Security Applications

Energy Technology Data Exchange (ETDEWEB)

Birch, Gabriel Carisle [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Dagel, Amber Lynn [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Kast, Brian A. [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Smith, Collin S. [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)

2015-09-01

Three-dimensional (3D) information in a physical security system is a highly useful dis- criminator. The two-dimensional data from an imaging systems fails to provide target dis- tance and three-dimensional motion vector, which can be used to reduce nuisance alarm rates and increase system effectiveness. However, 3D imaging devices designed primarily for use in physical security systems are uncommon. This report discusses an architecture favorable to physical security systems; an inexpensive snapshot 3D imaging system utilizing a simple illumination system. The method of acquiring 3D data, tests to understand illumination de- sign, and software modifications possible to maximize information gathering capability are discussed.

EISCAT Aperture Synthesis Imaging (EASI _3D) for the EISCAT_3D Project

Science.gov (United States)

La Hoz, Cesar; Belyey, Vasyl

2012-07-01

Aperture Synthesis Imaging Radar (ASIR) is one of the technologies adopted by the EISCAT_3D project to endow it with imaging capabilities in 3-dimensions that includes sub-beam resolution. Complemented by pulse compression, it will provide 3-dimensional images of certain types of incoherent scatter radar targets resolved to about 100 metres at 100 km range, depending on the signal-to-noise ratio. This ability will open new research opportunities to map small structures associated with non-homogeneous, unstable processes such as aurora, summer and winter polar radar echoes (PMSE and PMWE), Natural Enhanced Ion Acoustic Lines (NEIALs), structures excited by HF ionospheric heating, meteors, space debris, and others. The underlying physico-mathematical principles of the technique are the same as the technique employed in radioastronomy to image stellar objects; both require sophisticated inversion techniques to obtain reliable images.

Rainbow particle imaging velocimetry for dense 3D fluid velocity imaging

KAUST Repository

Xiong, Jinhui

2017-07-21

Despite significant recent progress, dense, time-resolved imaging of complex, non-stationary 3D flow velocities remains an elusive goal. In this work we tackle this problem by extending an established 2D method, Particle Imaging Velocimetry, to three dimensions by encoding depth into color. The encoding is achieved by illuminating the flow volume with a continuum of light planes (a

A fast, accurate, and automatic 2D-3D image registration for image-guided cranial radiosurgery

International Nuclear Information System (INIS)

Fu Dongshan; Kuduvalli, Gopinath

2008-01-01

The authors developed a fast and accurate two-dimensional (2D)-three-dimensional (3D) image registration method to perform precise initial patient setup and frequent detection and correction for patient movement during image-guided cranial radiosurgery treatment. In this method, an approximate geometric relationship is first established to decompose a 3D rigid transformation in the 3D patient coordinate into in-plane transformations and out-of-plane rotations in two orthogonal 2D projections. Digitally reconstructed radiographs are generated offline from a preoperative computed tomography volume prior to treatment and used as the reference for patient position. A multiphase framework is designed to register the digitally reconstructed radiographs with the x-ray images periodically acquired during patient setup and treatment. The registration in each projection is performed independently; the results in the two projections are then combined and converted to a 3D rigid transformation by 2D-3D geometric backprojection. The in-plane transformation and the out-of-plane rotation are estimated using different search methods, including multiresolution matching, steepest descent minimization, and one-dimensional search. Two similarity measures, optimized pattern intensity and sum of squared difference, are applied at different registration phases to optimize accuracy and computation speed. Various experiments on an anthropomorphic head-and-neck phantom showed that, using fiducial registration as a gold standard, the registration errors were 0.33±0.16 mm (s.d.) in overall translation and 0.29 deg. ±0.11 deg. (s.d.) in overall rotation. The total targeting errors were 0.34±0.16 mm (s.d.), 0.40±0.2 mm (s.d.), and 0.51±0.26 mm (s.d.) for the targets at the distances of 2, 6, and 10 cm from the rotation center, respectively. The computation time was less than 3 s on a computer with an Intel Pentium 3.0 GHz dual processor

Silhouette-based approach of 3D image reconstruction for automated image acquisition using robotic arm

Science.gov (United States)

Azhar, N.; Saad, W. H. M.; Manap, N. A.; Saad, N. M.; Syafeeza, A. R.

2017-06-01

This study presents the approach of 3D image reconstruction using an autonomous robotic arm for the image acquisition process. A low cost of the automated imaging platform is created using a pair of G15 servo motor connected in series to an Arduino UNO as a main microcontroller. Two sets of sequential images were obtained using different projection angle of the camera. The silhouette-based approach is used in this study for 3D reconstruction from the sequential images captured from several different angles of the object. Other than that, an analysis based on the effect of different number of sequential images on the accuracy of 3D model reconstruction was also carried out with a fixed projection angle of the camera. The effecting elements in the 3D reconstruction are discussed and the overall result of the analysis is concluded according to the prototype of imaging platform.

Signal alteration of the cochlear perilymph on 3 different sequences after intratympanic Gd-DTPA administration at 3 tesla. Comparison of 3D-FLAIR, 3D-T1-weighted imaging, and 3D-CISS

International Nuclear Information System (INIS)

Yamazaki, Masahiro; Naganawa, Shinji; Kawai, Hisashi; Nihashi, Takashi; Nakashima, Tsutomu

2010-01-01

Three-dimensional fluid-attenuated inversion recovery (3D-FLAIR) imaging after intratympanic gadolinium injection is useful for pathophysiologic and morphologic analysis of the inner ear. However, statistical analysis of differences in inner ear signal intensity among 3D-FLAIR and other sequences has not been reported. We evaluated the signal intensity of cochlear fluid on each of 3D-FLAIR, 3D-T 1 -weighted imaging (T 1 WI), and 3D-constructive interference in the steady state (CISS) to clarify the differences in contrast effect among these 3 sequences using intratympanic gadolinium injection. Twenty-one patients underwent 3D-FLAIR, 3D-T 1 WI, and 3D-CISS imaging at 3 tesla 24 hours after intratympanic injection of gadolinium. We determined regions of interest of the cochleae (C) and medulla oblongata (M) on each image, evaluated the signal intensity ratio between C and M (CM ratio), and determined the ratio of cochlear signal intensity of the injected side to that of the non-injected side (contrast value). The CM ratio of the injected side (3.00±1.31, range, 0.53 to 4.88, on 3D-FLAIR; 0.83±0.30, range, 0.36 to 1.58 on 3D-T 1 WI) was significantly higher than that of the non-injected side (0.52±0.14, range, 0.30 to 0.76 on 3D-FLAIR; 0.49±0.11, range, 0.30 to 0.71 on 3D-T 1 WI) on 3D-FLAIR and 3D-T 1 WI (P 1 WI (1.73±0.60 range, 0.98 to 3.09) (P<0.001). The 3D-FLAIR sequence is the most sensitive for observing alteration in inner ear fluid signal after intratympanic gadolinium injection. Our results warrant use of 3D-FLAIR as a sensitive imaging technique to clarify the pathological and morphological mechanisms of disorders of the inner ear. (author)

Effect of Post-Reconstruction Gaussian Filtering on Image Quality and Myocardial Blood Flow Measurement with N-13 Ammonia PET

Directory of Open Access Journals (Sweden)

Hyeon Sik Kim

2014-10-01

Full Text Available Objective(s: In order to evaluate the effect of post-reconstruction Gaussian filtering on image quality and myocardial blood flow (MBF measurement by dynamic N-13 ammonia positron emission tomography (PET, we compared various reconstruction and filtering methods with image characteristics. Methods: Dynamic PET images of three patients with coronary artery disease (male-female ratio of 2:1; age: 57, 53, and 76 years were reconstructed, using filtered back projection (FBP and ordered subset expectation maximization (OSEM methods. OSEM reconstruction consisted of OSEM_2I, OSEM_4I, and OSEM_6I with 2, 4, and 6 iterations, respectively. The images, reconstructed and filtered by Gaussian filters of 5, 10, and 15 mm, were obtained, as well as non-filtered images. Visual analysis of image quality (IQ was performed using a 3-grade scoring system by 2 independent readers, blinded to the reconstruction and filtering methods of stress images. Then, signal-to-noise ratio (SNR was calculated by noise and contrast recovery (CR. Stress and rest MBF and coronary flow reserve (CFR were obtained for each method. IQ scores, stress and rest MBF, and CFR were compared between the methods, using Chi-square and Kruskal-Wallis tests. Results: In the visual analysis, IQ was significantly higher by 10 mm Gaussian filtering, compared to other sizes of filter (PP=0.923 and 0.855 for readers 1 and 2, respectively. SNR was significantly higher in 10 mm Gaussian filter. There was a significant difference in stress and rest MBF between several vascular territories. However CFR was not significantly different according to various filtering methods. Conclusion: Post-reconstruction Gaussian filtering with a filter size of 10 mm significantly enhances the IQ of N-13 ammonia PET-CT, without changing the results of CFR calculation. .

Multithreaded real-time 3D image processing software architecture and implementation

Science.gov (United States)

Ramachandra, Vikas; Atanassov, Kalin; Aleksic, Milivoje; Goma, Sergio R.

2011-03-01

Recently, 3D displays and videos have generated a lot of interest in the consumer electronics industry. To make 3D capture and playback popular and practical, a user friendly playback interface is desirable. Towards this end, we built a real time software 3D video player. The 3D video player displays user captured 3D videos, provides for various 3D specific image processing functions and ensures a pleasant viewing experience. Moreover, the player enables user interactivity by providing digital zoom and pan functionalities. This real time 3D player was implemented on the GPU using CUDA and OpenGL. The player provides user interactive 3D video playback. Stereo images are first read by the player from a fast drive and rectified. Further processing of the images determines the optimal convergence point in the 3D scene to reduce eye strain. The rationale for this convergence point selection takes into account scene depth and display geometry. The first step in this processing chain is identifying keypoints by detecting vertical edges within the left image. Regions surrounding reliable keypoints are then located on the right image through the use of block matching. The difference in the positions between the corresponding regions in the left and right images are then used to calculate disparity. The extrema of the disparity histogram gives the scene disparity range. The left and right images are shifted based upon the calculated range, in order to place the desired region of the 3D scene at convergence. All the above computations are performed on one CPU thread which calls CUDA functions. Image upsampling and shifting is performed in response to user zoom and pan. The player also consists of a CPU display thread, which uses OpenGL rendering (quad buffers). This also gathers user input for digital zoom and pan and sends them to the processing thread.

3D fingerprint imaging system based on full-field fringe projection profilometry

Science.gov (United States)

Huang, Shujun; Zhang, Zonghua; Zhao, Yan; Dai, Jie; Chen, Chao; Xu, Yongjia; Zhang, E.; Xie, Lili

2014-01-01

As an unique, unchangeable and easily acquired biometrics, fingerprint has been widely studied in academics and applied in many fields over the years. The traditional fingerprint recognition methods are based on the obtained 2D feature of fingerprint. However, fingerprint is a 3D biological characteristic. The mapping from 3D to 2D loses 1D information and causes nonlinear distortion of the captured fingerprint. Therefore, it is becoming more and more important to obtain 3D fingerprint information for recognition. In this paper, a novel 3D fingerprint imaging system is presented based on fringe projection technique to obtain 3D features and the corresponding color texture information. A series of color sinusoidal fringe patterns with optimum three-fringe numbers are projected onto a finger surface. From another viewpoint, the fringe patterns are deformed by the finger surface and captured by a CCD camera. 3D shape data of the finger can be obtained from the captured fringe pattern images. This paper studies the prototype of the 3D fingerprint imaging system, including principle of 3D fingerprint acquisition, hardware design of the 3D imaging system, 3D calibration of the system, and software development. Some experiments are carried out by acquiring several 3D fingerprint data. The experimental results demonstrate the feasibility of the proposed 3D fingerprint imaging system.

Numerical modeling of time domain 3-D problems in accelerator physics

International Nuclear Information System (INIS)

Harfoush, F.A.; Jurgens, T.G.

1990-06-01

Time domain analysis is relevant in particle accelerators to study the electromagnetic field interaction of a moving source particle on a lagging test particle as the particles pass an accelerating cavity or some other structure. These fields are called wake fields. The travelling beam inside a beam pipe may undergo more complicated interactions with its environment due to the presence of other irregularities like wires, thin slots, joints and other types of obstacles. Analytical solutions of such problems is impossible and one has to resort to a numerical method. In this paper we present results of our first attempt to model these problems in 3-D using our finite difference time domain (FDTD) code. 10 refs., 9 figs

Rapid 3D Reconstruction for Image Sequence Acquired from UAV Camera

Directory of Open Access Journals (Sweden)

Yufu Qu

2018-01-01

Full Text Available In order to reconstruct three-dimensional (3D structures from an image sequence captured by unmanned aerial vehicles’ camera (UAVs and improve the processing speed, we propose a rapid 3D reconstruction method that is based on an image queue, considering the continuity and relevance of UAV camera images. The proposed approach first compresses the feature points of each image into three principal component points by using the principal component analysis method. In order to select the key images suitable for 3D reconstruction, the principal component points are used to estimate the interrelationships between images. Second, these key images are inserted into a fixed-length image queue. The positions and orientations of the images are calculated, and the 3D coordinates of the feature points are estimated using weighted bundle adjustment. With this structural information, the depth maps of these images can be calculated. Next, we update the image queue by deleting some of the old images and inserting some new images into the queue, and a structural calculation of all the images can be performed by repeating the previous steps. Finally, a dense 3D point cloud can be obtained using the depth–map fusion method. The experimental results indicate that when the texture of the images is complex and the number of images exceeds 100, the proposed method can improve the calculation speed by more than a factor of four with almost no loss of precision. Furthermore, as the number of images increases, the improvement in the calculation speed will become more noticeable.

Rapid 3D Reconstruction for Image Sequence Acquired from UAV Camera.

Science.gov (United States)

Qu, Yufu; Huang, Jianyu; Zhang, Xuan

2018-01-14

In order to reconstruct three-dimensional (3D) structures from an image sequence captured by unmanned aerial vehicles' camera (UAVs) and improve the processing speed, we propose a rapid 3D reconstruction method that is based on an image queue, considering the continuity and relevance of UAV camera images. The proposed approach first compresses the feature points of each image into three principal component points by using the principal component analysis method. In order to select the key images suitable for 3D reconstruction, the principal component points are used to estimate the interrelationships between images. Second, these key images are inserted into a fixed-length image queue. The positions and orientations of the images are calculated, and the 3D coordinates of the feature points are estimated using weighted bundle adjustment. With this structural information, the depth maps of these images can be calculated. Next, we update the image queue by deleting some of the old images and inserting some new images into the queue, and a structural calculation of all the images can be performed by repeating the previous steps. Finally, a dense 3D point cloud can be obtained using the depth-map fusion method. The experimental results indicate that when the texture of the images is complex and the number of images exceeds 100, the proposed method can improve the calculation speed by more than a factor of four with almost no loss of precision. Furthermore, as the number of images increases, the improvement in the calculation speed will become more noticeable.

3D T2-weighted imaging to shorten multiparametric prostate MRI protocols.

Science.gov (United States)

Polanec, Stephan H; Lazar, Mathias; Wengert, Georg J; Bickel, Hubert; Spick, Claudio; Susani, Martin; Shariat, Shahrokh; Clauser, Paola; Baltzer, Pascal A T

2018-04-01

To determine whether 3D acquisitions provide equivalent image quality, lesion delineation quality and PI-RADS v2 performance compared to 2D acquisitions in T2-weighted imaging of the prostate at 3 T. This IRB-approved, prospective study included 150 consecutive patients (mean age 63.7 years, 35-84 years; mean PSA 7.2 ng/ml, 0.4-31.1 ng/ml). Two uroradiologists (R1, R2) independently rated image quality and lesion delineation quality using a five-point ordinal scale and assigned a PI-RADS score for 2D and 3D T2-weighted image data sets. Data were compared using visual grading characteristics (VGC) and receiver operating characteristics (ROC)/area under the curve (AUC) analysis. Image quality was similarly good to excellent for 2D T2w (mean score R1, 4.3 ± 0.81; R2, 4.7 ± 0.83) and 3D T2w (mean score R1, 4.3 ± 0.82; R2, 4.7 ± 0.69), p = 0.269. Lesion delineation was rated good to excellent for 2D (mean score R1, 4.16 ± 0.81; R2, 4.19 ± 0.92) and 3D T2w (R1, 4.19 ± 0.94; R2, 4.27 ± 0.94) without significant differences (p = 0.785). ROC analysis showed an equivalent performance for 2D (AUC 0.580-0.623) and 3D (AUC 0.576-0.629) T2w (p > 0.05, respectively). Three-dimensional acquisitions demonstrated equivalent image and lesion delineation quality, and PI-RADS v2 performance, compared to 2D in T2-weighted imaging of the prostate. Three-dimensional T2-weighted imaging could be used to considerably shorten prostate MRI protocols in clinical practice. • 3D shows equivalent image quality and lesion delineation compared to 2D T2w. • 3D T2w and 2D T2w image acquisition demonstrated comparable diagnostic performance. • Using a single 3D T2w acquisition may shorten the protocol by 40%. • Combined with short DCE, multiparametric protocols of 10 min are feasible.

Remote laboratory for phase-aided 3D microscopic imaging and metrology

Science.gov (United States)

Wang, Meng; Yin, Yongkai; Liu, Zeyi; He, Wenqi; Li, Boqun; Peng, Xiang

2014-05-01

In this paper, the establishment of a remote laboratory for phase-aided 3D microscopic imaging and metrology is presented. Proposed remote laboratory consists of three major components, including the network-based infrastructure for remote control and data management, the identity verification scheme for user authentication and management, and the local experimental system for phase-aided 3D microscopic imaging and metrology. The virtual network computer (VNC) is introduced to remotely control the 3D microscopic imaging system. Data storage and management are handled through the open source project eSciDoc. Considering the security of remote laboratory, the fingerprint is used for authentication with an optical joint transform correlation (JTC) system. The phase-aided fringe projection 3D microscope (FP-3DM), which can be remotely controlled, is employed to achieve the 3D imaging and metrology of micro objects.

Research on GPU-accelerated algorithm in 3D finite difference neutron diffusion calculation method

International Nuclear Information System (INIS)

Xu Qi; Yu Ganglin; Wang Kan; Sun Jialong

2014-01-01

In this paper, the adaptability of the neutron diffusion numerical algorithm on GPUs was studied, and a GPU-accelerated multi-group 3D neutron diffusion code based on finite difference method was developed. The IAEA 3D PWR benchmark problem was calculated in the numerical test. The results demonstrate both high efficiency and adequate accuracy of the GPU implementation for neutron diffusion equation. (authors)

Subnuclear foci quantification using high-throughput 3D image cytometry

Science.gov (United States)

Wadduwage, Dushan N.; Parrish, Marcus; Choi, Heejin; Engelward, Bevin P.; Matsudaira, Paul; So, Peter T. C.

2015-07-01

Ionising radiation causes various types of DNA damages including double strand breaks (DSBs). DSBs are often recognized by DNA repair protein ATM which forms gamma-H2AX foci at the site of the DSBs that can be visualized using immunohistochemistry. However most of such experiments are of low throughput in terms of imaging and image analysis techniques. Most of the studies still use manual counting or classification. Hence they are limited to counting a low number of foci per cell (5 foci per nucleus) as the quantification process is extremely labour intensive. Therefore we have developed a high throughput instrumentation and computational pipeline specialized for gamma-H2AX foci quantification. A population of cells with highly clustered foci inside nuclei were imaged, in 3D with submicron resolution, using an in-house developed high throughput image cytometer. Imaging speeds as high as 800 cells/second in 3D were achieved by using HiLo wide-field depth resolved imaging and a remote z-scanning technique. Then the number of foci per cell nucleus were quantified using a 3D extended maxima transform based algorithm. Our results suggests that while most of the other 2D imaging and manual quantification studies can count only up to about 5 foci per nucleus our method is capable of counting more than 100. Moreover we show that 3D analysis is significantly superior compared to the 2D techniques.

3D-SIFT-Flow for atlas-based CT liver image segmentation.

Science.gov (United States)

Xu, Yan; Xu, Chenchao; Kuang, Xiao; Wang, Hongkai; Chang, Eric I-Chao; Huang, Weimin; Fan, Yubo

2016-05-01

In this paper, the authors proposed a new 3D registration algorithm, 3D-scale invariant feature transform (SIFT)-Flow, for multiatlas-based liver segmentation in computed tomography (CT) images. In the registration work, the authors developed a new registration method that takes advantage of dense correspondence using the informative and robust SIFT feature. The authors computed the dense SIFT features for the source image and the target image and designed an objective function to obtain the correspondence between these two images. Labeling of the source image was then mapped to the target image according to the former correspondence, resulting in accurate segmentation. In the fusion work, the 2D-based nonparametric label transfer method was extended to 3D for fusing the registered 3D atlases. Compared with existing registration algorithms, 3D-SIFT-Flow has its particular advantage in matching anatomical structures (such as the liver) that observe large variation/deformation. The authors observed consistent improvement over widely adopted state-of-the-art registration methods such as ELASTIX, ANTS, and multiatlas fusion methods such as joint label fusion. Experimental results of liver segmentation on the MICCAI 2007 Grand Challenge are encouraging, e.g., Dice overlap ratio 96.27% ± 0.96% by our method compared with the previous state-of-the-art result of 94.90% ± 2.86%. Experimental results show that 3D-SIFT-Flow is robust for segmenting the liver from CT images, which has large tissue deformation and blurry boundary, and 3D label transfer is effective and efficient for improving the registration accuracy.

3D visualization of medical images for personalized learning of human anatomy

NARCIS (Netherlands)

Laurence Alpay; Jelle Scheurleer; Harmen Bijwaard

2015-01-01

to be held in Lisbon/Portugal on October 15-17, 2015 Medical imaging nowadays often yields high definition 3D images (from CT, PET, MRI, etc.). Usually these images need to be evaluated on 2D monitors. In the transition from 3D to 2D the image becomes more difficult to interpret as a whole. To aid

2D VS 3D imaging of brain tumours with 18F-Fluoromisonidazole (FMISO) and positron emission tomography (PET)

International Nuclear Information System (INIS)

Pathmaraj, K.; Scott, A.M.; Egan, G.F.; Hannah, A.; Tauro, A.; Tochon-Danguy, A.; Sachinidis, J.; Berlangieri, S.U.; Fabinyi, G.; McKay, W.J.; Cher, L.

1998-01-01

Full text: 18 F-FMISO accumulates in hypoxic cells and can be used in the PET imaging of brain tumours containing viable but hypoxic cells. The limited activity (typically 130 MBq) of injected 18 F-FMISO yield poor statistics, requiring prolonged imaging in the conventional 2D mode of PET scanning. 3D (septa retracted) imaging allows for more counts to be collected over a shorter time period making it a more practical alternative. This study investigates the contrast resolution that can be obtained from 3D PET scans compared to the corresponding 2D scan. A patient recently diagnosed with brain tumour was injected with 18 -FMISO 2 hours prior to scanning and imaged supine on a 951/31R PET scanner with the head secured firmly in a head holder. The imaging protocol consisted of a 3 min emission rectilinear scan to position the brain in the FOV, a 10 min post-emission transmission scan, a 20 min 2D emission scan and a 5X10 min frames 3D emission scan. Both the 2D and 3D scans were reconstructed with filtered backprojection algorithm. The first 10 min frame of the 3D acquisition was reconstructed. The total true counts were 3 million and 6.06 million in the 2D image and 3D images respectively. The random events were 0.24 million and 0.96 million in the 2D and 3D images respectively. The Noise Equivalent Counts (NEC) were 2.2 million and 2.02 million for the 2D and 3D images respectively indicating that the 2D and 3D scans (in spite of the nominal true events being vastly different in the 2 scans) had similar Signal to Noise Ratio (SNR). Circular ROI's were defined in the tumour and the contralateral cortex in comparable transaxial slices of the 2D and 3D images. Contrast resolution of the tumour to the background was calculated as 1.4 and 1.38 in the 2D and 3D images respectively. Thus comparable contrast resolution is obtained in the brain with both 3D and 2D images, making 3D imaging a viable alternative to 2D imaging and greatly reducing imaging time. Optimum time

Comparison of 3D cube FLAIR with 2D FLAIR for multiple sclerosis imaging at 3 tesla

Energy Technology Data Exchange (ETDEWEB)

Patzig, M.; Brueckmann, H.; Fesl, G. [Muenchen Univ. (Germany). Dept. of Neuroradiology; Burke, M. [GE Healthcare, Solingen (Germany)

2014-05-15

Purpose: Three-dimensional (3 D) MRI sequences allow improved spatial resolution with good signal and contrast properties as well as multiplanar reconstruction. We sought to compare Cube, a 3 D FLAIR sequence, to a standard 2 D FLAIR sequence in multiple sclerosis (MS) imaging. Materials and Methods: Examinations were performed in the clinical routine on a 3.0 Tesla scanner. 12 patients with definite MS were included. Lesions with MS-typical properties on the images of Cube FLAIR and 2 D FLAIR sequences were counted and allocated to different brain regions. Signal-to-noise ratios (SNR) and contrast-to-noise ratios (CNR) were calculated. Results: With 384 the overall number of lesions found with Cube FLAIR was significantly higher than with 2 D FLAIR (N = 221). The difference was mostly accounted for by supratentorial lesions (N = 372 vs. N = 216) while the infratentorial lesion counts were low in both sequences. SNRs and CNRs were significantly higher in CUBE FLAIR with the exception of the CNR of lesion to gray matter, which was not significantly different. Conclusion: Cube FLAIR showed a higher sensitivity for MS lesions compared to a 2 D FLAIR sequence. 3 D FLAIR might replace 2 D FLAIR sequences in MS imaging in the future. (orig.)

An efficient implementation of 3D high-resolution imaging for large-scale seismic data with GPU/CPU heterogeneous parallel computing

Science.gov (United States)

Xu, Jincheng; Liu, Wei; Wang, Jin; Liu, Linong; Zhang, Jianfeng

2018-02-01

De-absorption pre-stack time migration (QPSTM) compensates for the absorption and dispersion of seismic waves by introducing an effective Q parameter, thereby making it an effective tool for 3D, high-resolution imaging of seismic data. Although the optimal aperture obtained via stationary-phase migration reduces the computational cost of 3D QPSTM and yields 3D stationary-phase QPSTM, the associated computational efficiency is still the main problem in the processing of 3D, high-resolution images for real large-scale seismic data. In the current paper, we proposed a division method for large-scale, 3D seismic data to optimize the performance of stationary-phase QPSTM on clusters of graphics processing units (GPU). Then, we designed an imaging point parallel strategy to achieve an optimal parallel computing performance. Afterward, we adopted an asynchronous double buffering scheme for multi-stream to perform the GPU/CPU parallel computing. Moreover, several key optimization strategies of computation and storage based on the compute unified device architecture (CUDA) were adopted to accelerate the 3D stationary-phase QPSTM algorithm. Compared with the initial GPU code, the implementation of the key optimization steps, including thread optimization, shared memory optimization, register optimization and special function units (SFU), greatly improved the efficiency. A numerical example employing real large-scale, 3D seismic data showed that our scheme is nearly 80 times faster than the CPU-QPSTM algorithm. Our GPU/CPU heterogeneous parallel computing framework significant reduces the computational cost and facilitates 3D high-resolution imaging for large-scale seismic data.

Phase aided 3D imaging and modeling: dedicated systems and case studies

Science.gov (United States)

Yin, Yongkai; He, Dong; Liu, Zeyi; Liu, Xiaoli; Peng, Xiang

2014-05-01

Dedicated prototype systems for 3D imaging and modeling (3DIM) are presented. The 3D imaging systems are based on the principle of phase-aided active stereo, which have been developed in our laboratory over the past few years. The reported 3D imaging prototypes range from single 3D sensor to a kind of optical measurement network composed of multiple node 3D-sensors. To enable these 3D imaging systems, we briefly discuss the corresponding calibration techniques for both single sensor and multi-sensor optical measurement network, allowing good performance of the 3DIM prototype systems in terms of measurement accuracy and repeatability. Furthermore, two case studies including the generation of high quality color model of movable cultural heritage and photo booth from body scanning are presented to demonstrate our approach.

2D imaging and 3D sensing data acquisition and mutual registration for painting conservation

Science.gov (United States)

Fontana, Raffaella; Gambino, Maria Chiara; Greco, Marinella; Marras, Luciano; Pampaloni, Enrico M.; Pelagotti, Anna; Pezzati, Luca; Poggi, Pasquale

2005-01-01

We describe the application of 2D and 3D data acquisition and mutual registration to the conservation of paintings. RGB color image acquisition, IR and UV fluorescence imaging, together with the more recent hyperspectral imaging (32 bands) are among the most useful techniques in this field. They generally are meant to provide information on the painting materials, on the employed techniques and on the object state of conservation. However, only when the various images are perfectly registered on each other and on the 3D model, no ambiguity is possible and safe conclusions may be drawn. We present the integration of 2D and 3D measurements carried out on two different paintings: "Madonna of the Yarnwinder" by Leonardo da Vinci, and "Portrait of Lionello d'Este", by Pisanello, both painted in the XV century.

An automated 3D reconstruction method of UAV images

Science.gov (United States)

Liu, Jun; Wang, He; Liu, Xiaoyang; Li, Feng; Sun, Guangtong; Song, Ping

2015-10-01

In this paper a novel fully automated 3D reconstruction approach based on low-altitude unmanned aerial vehicle system (UAVs) images will be presented, which does not require previous camera calibration or any other external prior knowledge. Dense 3D point clouds are generated by integrating orderly feature extraction, image matching, structure from motion (SfM) and multi-view stereo (MVS) algorithms, overcoming many of the cost, time limitations of rigorous photogrammetry techniques. An image topology analysis strategy is introduced to speed up large scene reconstruction by taking advantage of the flight-control data acquired by UAV. Image topology map can significantly reduce the running time of feature matching by limiting the combination of images. A high-resolution digital surface model of the study area is produced base on UAV point clouds by constructing the triangular irregular network. Experimental results show that the proposed approach is robust and feasible for automatic 3D reconstruction of low-altitude UAV images, and has great potential for the acquisition of spatial information at large scales mapping, especially suitable for rapid response and precise modelling in disaster emergency.

AUTOMATIC 3D MAPPING USING MULTIPLE UNCALIBRATED CLOSE RANGE IMAGES

Directory of Open Access Journals (Sweden)

M. Rafiei

2013-09-01

Full Text Available Automatic three-dimensions modeling of the real world is an important research topic in the geomatics and computer vision fields for many years. By development of commercial digital cameras and modern image processing techniques, close range photogrammetry is vastly utilized in many fields such as structure measurements, topographic surveying, architectural and archeological surveying, etc. A non-contact photogrammetry provides methods to determine 3D locations of objects from two-dimensional (2D images. Problem of estimating the locations of 3D points from multiple images, often involves simultaneously estimating both 3D geometry (structure and camera pose (motion, it is commonly known as structure from motion (SfM. In this research a step by step approach to generate the 3D point cloud of a scene is considered. After taking images with a camera, we should detect corresponding points in each two views. Here an efficient SIFT method is used for image matching for large baselines. After that, we must retrieve the camera motion and 3D position of the matched feature points up to a projective transformation (projective reconstruction. Lacking additional information on the camera or the scene makes the parallel lines to be unparalleled. The results of SfM computation are much more useful if a metric reconstruction is obtained. Therefor multiple views Euclidean reconstruction applied and discussed. To refine and achieve the precise 3D points we use more general and useful approach, namely bundle adjustment. At the end two real cases have been considered to reconstruct (an excavation and a tower.

Performance evaluation of 3-D enhancement filters for detection of lung cancer from 3-D chest X-ray CT images

International Nuclear Information System (INIS)

Shimizu, Akinobu; Hagai, Makoto; Toriwaki, Jun-ichiro; Hasegawa, Jun-ichi.

1995-01-01

This paper evaluates the performance of several three dimensional enhancement filters used in procedures for detecting lung cancer shadows from three dimensional (3D) chest X-ray CT images. Two dimensional enhancement filters such as Min-DD filter, Contrast filter and N-Quoit filter have been proposed for enhancing cancer shadows in conventional 2D X-ray images. In this paper, we extend each of these 2D filters to a 3D filter and evaluate its performance experimentally by using CT images with artificial and true lung cancer shadows. As a result, we find that these 3D filters are effective for determining the position of a lung cancer shadow in a 3D chest CT image, as compared with the simple procedure such as smoothing filter, and that the performance of these filters become lower in the hilar area due to the influence of the vessel shadows. (author)

Reconstruction, Processing and Display of 3D-Images

International Nuclear Information System (INIS)

Lenz, R.

1986-01-01

In the last few years a number of methods have been developed which can produce true 3D images, volumes of density values. We review two of these techniques (confocal microscopy and X-ray tomography) which were used in the reconstruction of some of our images. The other images came from transmission electron microscopes, gammacameras and magnetic resonance scanners. A new algorithm is suggested which uses projection onto convex sets to improve the depth resolution in the microscopy case. Since we use a TV-monitor as display device we have to project 3D volumes to 2D images. We use the following type of projections: reprojections, range images, colorcoded depth and shaded surface displays. Shaded surface displays use the surface gradient to compute the gray value in the projection. We describe how this gradient can be computed from the range image and from the original density volume. Normally we compute a whole series of projections where the volume is rotated some degrees between two projections. In a separate display session we can display these images in stereo and motion. We describe how noise reduction filters, gray value transformations, geometric manipulations, gradient filters, texture filters and binary techniques can be used to remove uninteresting points from the volume. Finally, a filter design strategy is developed which is based on the optimal basis function approach by Hummel. We show that for a large class of patterns, in images of arbitrary dimensions, the optimal basis functions are rotation-invariant operators as introduced by Danielsson in the 2D case. We also describe how the orientation of a pattern can be computed from its feature vector. (With 107 refs.) (author)

Image-Based 3D Face Modeling System

Directory of Open Access Journals (Sweden)

Vladimir Vezhnevets

2005-08-01

Full Text Available This paper describes an automatic system for 3D face modeling using frontal and profile images taken by an ordinary digital camera. The system consists of four subsystems including frontal feature detection, profile feature detection, shape deformation, and texture generation modules. The frontal and profile feature detection modules automatically extract the facial parts such as the eye, nose, mouth, and ear. The shape deformation module utilizes the detected features to deform the generic head mesh model such that the deformed model coincides with the detected features. A texture is created by combining the facial textures augmented from the input images and the synthesized texture and mapped onto the deformed generic head model. This paper provides a practical system for 3D face modeling, which is highly automated by aggregating, customizing, and optimizing a bunch of individual computer vision algorithms. The experimental results show a highly automated process of modeling, which is sufficiently robust to various imaging conditions. The whole model creation including all the optional manual corrections takes only 2∼3 minutes.

3D-SIFT-Flow for atlas-based CT liver image segmentation

Energy Technology Data Exchange (ETDEWEB)

Xu, Yan, E-mail: xuyan04@gmail.com [State Key Laboratory of Software Development Environment and Key Laboratory of Biomechanics and Mechanobiology of Ministry of Education, Beihang University, Beijing 100191, China and Research Institute of Beihang University in Shenzhen and Microsoft Research, Beijing 100080 (China); Xu, Chenchao, E-mail: chenchaoxu33@gmail.com; Kuang, Xiao, E-mail: kuangxiao.ace@gmail.com [School of Biological Science and Medical Engineering, Beihang University, Beijing 100191 (China); Wang, Hongkai, E-mail: wang.hongkai@gmail.com [Department of Biomedical Engineering, Dalian University of Technology, Dalian 116024 (China); Chang, Eric I-Chao, E-mail: eric.chang@microsoft.com [Microsoft Research, Beijing 100080 (China); Huang, Weimin, E-mail: wmhuang@i2r.a-star.edu.sg [Institute for Infocomm Research (I2R), Singapore 138632 (Singapore); Fan, Yubo, E-mail: yubofan@buaa.edu.cn [Key Laboratory of Biomechanics and Mechanobiology of Ministry of Education, Beihang University, Beijing 100191 (China)

2016-05-15

Purpose: In this paper, the authors proposed a new 3D registration algorithm, 3D-scale invariant feature transform (SIFT)-Flow, for multiatlas-based liver segmentation in computed tomography (CT) images. Methods: In the registration work, the authors developed a new registration method that takes advantage of dense correspondence using the informative and robust SIFT feature. The authors computed the dense SIFT features for the source image and the target image and designed an objective function to obtain the correspondence between these two images. Labeling of the source image was then mapped to the target image according to the former correspondence, resulting in accurate segmentation. In the fusion work, the 2D-based nonparametric label transfer method was extended to 3D for fusing the registered 3D atlases. Results: Compared with existing registration algorithms, 3D-SIFT-Flow has its particular advantage in matching anatomical structures (such as the liver) that observe large variation/deformation. The authors observed consistent improvement over widely adopted state-of-the-art registration methods such as ELASTIX, ANTS, and multiatlas fusion methods such as joint label fusion. Experimental results of liver segmentation on the MICCAI 2007 Grand Challenge are encouraging, e.g., Dice overlap ratio 96.27% ± 0.96% by our method compared with the previous state-of-the-art result of 94.90% ± 2.86%. Conclusions: Experimental results show that 3D-SIFT-Flow is robust for segmenting the liver from CT images, which has large tissue deformation and blurry boundary, and 3D label transfer is effective and efficient for improving the registration accuracy.

3D-SIFT-Flow for atlas-based CT liver image segmentation

International Nuclear Information System (INIS)

Xu, Yan; Xu, Chenchao; Kuang, Xiao; Wang, Hongkai; Chang, Eric I-Chao; Huang, Weimin; Fan, Yubo

2016-01-01

Purpose: In this paper, the authors proposed a new 3D registration algorithm, 3D-scale invariant feature transform (SIFT)-Flow, for multiatlas-based liver segmentation in computed tomography (CT) images. Methods: In the registration work, the authors developed a new registration method that takes advantage of dense correspondence using the informative and robust SIFT feature. The authors computed the dense SIFT features for the source image and the target image and designed an objective function to obtain the correspondence between these two images. Labeling of the source image was then mapped to the target image according to the former correspondence, resulting in accurate segmentation. In the fusion work, the 2D-based nonparametric label transfer method was extended to 3D for fusing the registered 3D atlases. Results: Compared with existing registration algorithms, 3D-SIFT-Flow has its particular advantage in matching anatomical structures (such as the liver) that observe large variation/deformation. The authors observed consistent improvement over widely adopted state-of-the-art registration methods such as ELASTIX, ANTS, and multiatlas fusion methods such as joint label fusion. Experimental results of liver segmentation on the MICCAI 2007 Grand Challenge are encouraging, e.g., Dice overlap ratio 96.27% ± 0.96% by our method compared with the previous state-of-the-art result of 94.90% ± 2.86%. Conclusions: Experimental results show that 3D-SIFT-Flow is robust for segmenting the liver from CT images, which has large tissue deformation and blurry boundary, and 3D label transfer is effective and efficient for improving the registration accuracy.

Simultaneous 3D–2D image registration and C-arm calibration: Application to endovascular image-guided interventions

Energy Technology Data Exchange (ETDEWEB)

Mitrović, Uroš [Faculty of Electrical Engineering, University of Ljubljana, Tržaška 25, Ljubljana 1000, Slovenia and Cosylab, Control System Laboratory, Teslova ulica 30, Ljubljana 1000 (Slovenia); Pernuš, Franjo [Faculty of Electrical Engineering, University of Ljubljana, Tržaška 25, Ljubljana 1000 (Slovenia); Likar, Boštjan; Špiclin, Žiga, E-mail: ziga.spiclin@fe.uni-lj.si [Faculty of Electrical Engineering, University of Ljubljana, Tržaška 25, Ljubljana 1000, Slovenia and Sensum, Computer Vision Systems, Tehnološki Park 21, Ljubljana 1000 (Slovenia)

2015-11-15

Purpose: Three-dimensional to two-dimensional (3D–2D) image registration is a key to fusion and simultaneous visualization of valuable information contained in 3D pre-interventional and 2D intra-interventional images with the final goal of image guidance of a procedure. In this paper, the authors focus on 3D–2D image registration within the context of intracranial endovascular image-guided interventions (EIGIs), where the 3D and 2D images are generally acquired with the same C-arm system. The accuracy and robustness of any 3D–2D registration method, to be used in a clinical setting, is influenced by (1) the method itself, (2) uncertainty of initial pose of the 3D image from which registration starts, (3) uncertainty of C-arm’s geometry and pose, and (4) the number of 2D intra-interventional images used for registration, which is generally one and at most two. The study of these influences requires rigorous and objective validation of any 3D–2D registration method against a highly accurate reference or “gold standard” registration, performed on clinical image datasets acquired in the context of the intervention. Methods: The registration process is split into two sequential, i.e., initial and final, registration stages. The initial stage is either machine-based or template matching. The latter aims to reduce possibly large in-plane translation errors by matching a projection of the 3D vessel model and 2D image. In the final registration stage, four state-of-the-art intrinsic image-based 3D–2D registration methods, which involve simultaneous refinement of rigid-body and C-arm parameters, are evaluated. For objective validation, the authors acquired an image database of 15 patients undergoing cerebral EIGI, for which accurate gold standard registrations were established by fiducial marker coregistration. Results: Based on target registration error, the obtained success rates of 3D to a single 2D image registration after initial machine-based and

3D visualisation of the middle ear and adjacent structures using reconstructed multi-slice CT datasets, correlating 3D images and virtual endoscopy to the 2D cross-sectional images

International Nuclear Information System (INIS)

Rodt, T.; Ratiu, P.; Kacher, D.F.; Anderson, M.; Jolesz, F.A.; Kikinis, R.; Becker, H.; Bartling, S.

2002-01-01

The 3D imaging of the middle ear facilitates better understanding of the patient's anatomy. Cross-sectional slices, however, often allow a more accurate evaluation of anatomical structures, as some detail may be lost through post-processing. In order to demonstrate the advantages of combining both approaches, we performed computed tomography (CT) imaging in two normal and 15 different pathological cases, and the 3D models were correlated to the cross-sectional CT slices. Reconstructed CT datasets were acquired by multi-slice CT. Post-processing was performed using the in-house software ''3D Slicer'', applying thresholding and manual segmentation. 3D models of the individual anatomical structures were generated and displayed in different colours. The display of relevant anatomical and pathological structures was evaluated in the greyscale 2D slices, 3D images, and the 2D slices showing the segmented 2D anatomy in different colours for each structure. Correlating 2D slices to the 3D models and virtual endoscopy helps to combine the advantages of each method. As generating 3D models can be extremely time-consuming, this approach can be a clinically applicable way of gaining a 3D understanding of the patient's anatomy by using models as a reference. Furthermore, it can help radiologists and otolaryngologists evaluating the 2D slices by adding the correct 3D information that would otherwise have to be mentally integrated. The method can be applied to radiological diagnosis, surgical planning, and especially, to teaching. (orig.)

A multi-frequency electrical impedance tomography system for real-time 2D and 3D imaging

Science.gov (United States)

Yang, Yunjie; Jia, Jiabin

2017-08-01

This paper presents the design and evaluation of a configurable, fast multi-frequency Electrical Impedance Tomography (mfEIT) system for real-time 2D and 3D imaging, particularly for biomedical imaging. The system integrates 32 electrode interfaces and the current frequency ranges from 10 kHz to 1 MHz. The system incorporates the following novel features. First, a fully adjustable multi-frequency current source with current monitoring function is designed. Second, a flexible switching scheme is developed for arbitrary sensing configuration and a semi-parallel data acquisition architecture is implemented for high-frame-rate data acquisition. Furthermore, multi-frequency digital quadrature demodulation is accomplished in a high-capacity Field Programmable Gate Array. At last, a 3D imaging software, visual tomography, is developed for real-time 2D and 3D image reconstruction, data analysis, and visualization. The mfEIT system is systematically tested and evaluated from the aspects of signal to noise ratio (SNR), frame rate, and 2D and 3D multi-frequency phantom imaging. The highest SNR is 82.82 dB on a 16-electrode sensor. The frame rate is up to 546 fps at serial mode and 1014 fps at semi-parallel mode. The evaluation results indicate that the presented mfEIT system is a powerful tool for real-time 2D and 3D imaging.

Real-time beam profile imaging system for actinotherapy accelerator

International Nuclear Information System (INIS)

Lin Yong; Wang Jingjin; Song Zheng; Zheng Putang; Wang Jianguo

2003-01-01

This paper describes a real-time beam profile imaging system for actinotheraphy accelerator. With the flash X-ray imager and the technique of digital image processing, a real-time 3-dimension dosage image is created from the intensity profile of the accelerator beam in real time. This system helps to obtain all the physical characters of the beam in any section plane, such as FWHM, penumbra, peak value, symmetry and homogeneity. This system has been used to acquire a 3-dimension dosage distribution of dynamic wedge modulator and the transient process of beam dosage. The system configure and the tested beam profile images are also presented

Estimating 3D tilt from local image cues in natural scenes

OpenAIRE

Burge, Johannes; McCann, Brian C.; Geisler, Wilson S.

2016-01-01

Estimating three-dimensional (3D) surface orientation (slant and tilt) is an important first step toward estimating 3D shape. Here, we examine how three local image cues from the same location (disparity gradient, luminance gradient, and dominant texture orientation) should be combined to estimate 3D tilt in natural scenes. We collected a database of natural stereoscopic images with precisely co-registered range images that provide the ground-truth distance at each pixel location. We then ana...

3D Point Cloud Reconstruction from Single Plenoptic Image

Directory of Open Access Journals (Sweden)

F. Murgia

2016-06-01

Full Text Available Novel plenoptic cameras sample the light field crossing the main camera lens. The information available in a plenoptic image must be processed, in order to create the depth map of the scene from a single camera shot. In this paper a novel algorithm, for the reconstruction of 3D point cloud of the scene from a single plenoptic image, taken with a consumer plenoptic camera, is proposed. Experimental analysis is conducted on several test images, and results are compared with state of the art methodologies. The results are very promising, as the quality of the 3D point cloud from plenoptic image, is comparable with the quality obtained with current non-plenoptic methodologies, that necessitate more than one image.

3D/3D registration of coronary CTA and biplane XA reconstructions for improved image guidance

Energy Technology Data Exchange (ETDEWEB)

Dibildox, Gerardo, E-mail: g.dibildox@erasmusmc.nl; Baka, Nora; Walsum, Theo van [Biomedical Imaging Group Rotterdam, Departments of Radiology and Medical Informatics, Erasmus Medical Center, 3015 GE Rotterdam (Netherlands); Punt, Mark; Aben, Jean-Paul [Pie Medical Imaging, 6227 AJ Maastricht (Netherlands); Schultz, Carl [Department of Cardiology, Erasmus Medical Center, 3015 GE Rotterdam (Netherlands); Niessen, Wiro [Quantitative Imaging Group, Faculty of Applied Sciences, Delft University of Technology, 2628 CJ Delft, The Netherlands and Biomedical Imaging Group Rotterdam, Departments of Radiology and Medical Informatics, Erasmus Medical Center, 3015 GE Rotterdam (Netherlands)

2014-09-15

Purpose: The authors aim to improve image guidance during percutaneous coronary interventions of chronic total occlusions (CTO) by providing information obtained from computed tomography angiography (CTA) to the cardiac interventionist. To this end, the authors investigate a method to register a 3D CTA model to biplane reconstructions. Methods: The authors developed a method for registering preoperative coronary CTA with intraoperative biplane x-ray angiography (XA) images via 3D models of the coronary arteries. The models are extracted from the CTA and biplane XA images, and are temporally aligned based on CTA reconstruction phase and XA ECG signals. Rigid spatial alignment is achieved with a robust probabilistic point set registration approach using Gaussian mixture models (GMMs). This approach is extended by including orientation in the Gaussian mixtures and by weighting bifurcation points. The method is evaluated on retrospectively acquired coronary CTA datasets of 23 CTO patients for which biplane XA images are available. Results: The Gaussian mixture model approach achieved a median registration accuracy of 1.7 mm. The extended GMM approach including orientation was not significantly different (P > 0.1) but did improve robustness with regards to the initialization of the 3D models. Conclusions: The authors demonstrated that the GMM approach can effectively be applied to register CTA to biplane XA images for the purpose of improving image guidance in percutaneous coronary interventions.

Imaging system for creating 3D block-face cryo-images of whole mice

Science.gov (United States)

Roy, Debashish; Breen, Michael; Salvado, Olivier; Heinzel, Meredith; McKinley, Eliot; Wilson, David

2006-03-01

We developed a cryomicrotome/imaging system that provides high resolution, high sensitivity block-face images of whole mice or excised organs, and applied it to a variety of biological applications. With this cryo-imaging system, we sectioned cryo-preserved tissues at 2-40 μm thickness and acquired high resolution brightfield and fluorescence images with microscopic in-plane resolution (as good as 1.2 μm). Brightfield images of normal and pathological anatomy show exquisite detail, especially in the abdominal cavity. Multi-planar reformatting and 3D renderings allow one to interrogate 3D structures. In this report, we present brightfield images of mouse anatomy, as well as 3D renderings of organs. For BPK mice model of polycystic kidney disease, we compared brightfield cryo-images and kidney volumes to MRI. The color images provided greater contrast and resolution of cysts as compared to in vivo MRI. We note that color cryo-images are closer to what a researcher sees in dissection, making it easier for them to interpret image data. The combination of field of view, depth of field, ultra high resolution and color/fluorescence contrast enables cryo-image volumes to provide details that cannot be found through in vivo imaging or other ex vivo optical imaging approaches. We believe that this novel imaging system will have applications that include identification of mouse phenotypes, characterization of diseases like blood vessel disease, kidney disease, and cancer, assessment of drug and gene therapy delivery and efficacy and validation of other imaging modalities.

Clinical evaluation of PET image reconstruction using a spatial resolution model

DEFF Research Database (Denmark)

Andersen, Flemming Littrup; Klausen, Thomas Levin; Loft, Annika

2013-01-01

PURPOSE: PET image resolution is variable across the measured field-of-view and described by the point spread function (PSF). When accounting for the PSF during PET image reconstruction image resolution is improved and partial volume effects are reduced. Here, we evaluate the effect of PSF......-based reconstruction on lesion quantification in routine clinical whole-body (WB) PET/CT imaging. MATERIALS AND METHODS: 41 oncology patients were referred for a WB-PET/CT examination (Biograph 40 TruePoint). Emission data were acquired at 2.5min/bed at 1hpi of 400 MBq [18F]-FDG. Attenuation-corrected PET images were...... reconstructed on 336×336-matrices using: (R1) standard AW-OSEM (4 iter, 8 subsets, 4mm Gaussian) and (R2) AW-OSEM with PSF (3 iter, 21 subsets, 2mm). Blinded and randomised reading of R1- and R2-PET images was performed. Individual lesions were located and counted independently on both sets of images...

OpenCL-Based FPGA Accelerator for 3D FDTD with Periodic and Absorbing Boundary Conditions

Directory of Open Access Journals (Sweden)

Hasitha Muthumala Waidyasooriya

2017-01-01

Full Text Available Finite difference time domain (FDTD method is a very poplar way of numerically solving partial differential equations. FDTD has a low operational intensity so that the performances in CPUs and GPUs are often restricted by the memory bandwidth. Recently, deeply pipelined FPGA accelerators have shown a lot of success by exploiting streaming data flows in FDTD computation. In spite of this success, many FPGA accelerators are not suitable for real-world applications that contain complex boundary conditions. Boundary conditions break the regularity of the data flow, so that the performances are significantly reduced. This paper proposes an FPGA accelerator that computes commonly used absorbing and periodic boundary conditions in many 3D FDTD applications. Accelerator is designed using a “C-like” programming language called OpenCL (open computing language. As a result, the proposed accelerator can be customized easily by changing the software code. According to the experimental results, we achieved over 3.3 times and 1.5 times higher processing speed compared to the CPUs and GPUs, respectively. Moreover, the proposed accelerator is more than 14 times faster compared to the recently proposed FPGA accelerators that are capable of handling complex boundary conditions.

Estimating 3D Object Parameters from 2D Grey-Level Images

NARCIS (Netherlands)

Houkes, Z.

2000-01-01

This thesis describes a general framework for parameter estimation, which is suitable for computer vision applications. The approach described combines 3D modelling, animation and estimation tools to determine parameters of objects in a scene from 2D grey-level images. The animation tool predicts

SU-D-BRA-03: Analysis of Systematic Errors with 2D/3D Image Registration for Target Localization and Treatment Delivery in Stereotactic Radiosurgery

International Nuclear Information System (INIS)

Xu, H; Chetty, I; Wen, N

2016-01-01

Purpose: Determine systematic deviations between 2D/3D and 3D/3D image registrations with six degrees of freedom (6DOF) for various imaging modalities and registration algorithms on the Varian Edge Linac. Methods: The 6DOF systematic errors were assessed by comparing automated 2D/3D (kV/MV vs. CT) with 3D/3D (CBCT vs. CT) image registrations from different imaging pairs, CT slice thicknesses, couch angles, similarity measures, etc., using a Rando head and a pelvic phantom. The 2D/3D image registration accuracy was evaluated at different treatment sites (intra-cranial and extra-cranial) by statistically analyzing 2D/3D pre-treatment verification against 3D/3D localization of 192 Stereotactic Radiosurgery/Stereotactic Body Radiation Therapy treatment fractions for 88 patients. Results: The systematic errors of 2D/3D image registration using kV-kV, MV-kV and MV-MV image pairs using 0.8 mm slice thickness CT images were within 0.3 mm and 0.3° for translations and rotations with a 95% confidence interval (CI). No significant difference between 2D/3D and 3D/3D image registrations (P>0.05) was observed for target localization at various CT slice thicknesses ranging from 0.8 to 3 mm. Couch angles (30, 45, 60 degree) did not impact the accuracy of 2D/3D image registration. Using pattern intensity with content image filtering was recommended for 2D/3D image registration to achieve the best accuracy. For the patient study, translational error was within 2 mm and rotational error was within 0.6 degrees in terms of 95% CI for 2D/3D image registration. For intra-cranial sites, means and std. deviations of translational errors were −0.2±0.7, 0.04±0.5, 0.1±0.4 mm for LNG, LAT, VRT directions, respectively. For extra-cranial sites, means and std. deviations of translational errors were - 0.04±1, 0.2±1, 0.1±1 mm for LNG, LAT, VRT directions, respectively. 2D/3D image registration uncertainties for intra-cranial and extra-cranial sites were comparable. Conclusion: The Varian

SU-D-BRA-03: Analysis of Systematic Errors with 2D/3D Image Registration for Target Localization and Treatment Delivery in Stereotactic Radiosurgery

Energy Technology Data Exchange (ETDEWEB)

Xu, H [Wayne State University, Detroit, MI (United States); Chetty, I; Wen, N [Henry Ford Health System, Detroit, MI (United States)

2016-06-15

Purpose: Determine systematic deviations between 2D/3D and 3D/3D image registrations with six degrees of freedom (6DOF) for various imaging modalities and registration algorithms on the Varian Edge Linac. Methods: The 6DOF systematic errors were assessed by comparing automated 2D/3D (kV/MV vs. CT) with 3D/3D (CBCT vs. CT) image registrations from different imaging pairs, CT slice thicknesses, couch angles, similarity measures, etc., using a Rando head and a pelvic phantom. The 2D/3D image registration accuracy was evaluated at different treatment sites (intra-cranial and extra-cranial) by statistically analyzing 2D/3D pre-treatment verification against 3D/3D localization of 192 Stereotactic Radiosurgery/Stereotactic Body Radiation Therapy treatment fractions for 88 patients. Results: The systematic errors of 2D/3D image registration using kV-kV, MV-kV and MV-MV image pairs using 0.8 mm slice thickness CT images were within 0.3 mm and 0.3° for translations and rotations with a 95% confidence interval (CI). No significant difference between 2D/3D and 3D/3D image registrations (P>0.05) was observed for target localization at various CT slice thicknesses ranging from 0.8 to 3 mm. Couch angles (30, 45, 60 degree) did not impact the accuracy of 2D/3D image registration. Using pattern intensity with content image filtering was recommended for 2D/3D image registration to achieve the best accuracy. For the patient study, translational error was within 2 mm and rotational error was within 0.6 degrees in terms of 95% CI for 2D/3D image registration. For intra-cranial sites, means and std. deviations of translational errors were −0.2±0.7, 0.04±0.5, 0.1±0.4 mm for LNG, LAT, VRT directions, respectively. For extra-cranial sites, means and std. deviations of translational errors were - 0.04±1, 0.2±1, 0.1±1 mm for LNG, LAT, VRT directions, respectively. 2D/3D image registration uncertainties for intra-cranial and extra-cranial sites were comparable. Conclusion: The Varian

Simultaneous multislice echo planar imaging with blipped controlled aliasing in parallel imaging results in higher acceleration: a promising technique for accelerated diffusion tensor imaging of skeletal muscle

OpenAIRE

Filli, Lukas; Piccirelli, Marco; Kenkel, David; Guggenberger, Roman; Andreisek, Gustav; Beck, Thomas; Runge, Val M; Boss, Andreas

2015-01-01

PURPOSE The aim of this study was to investigate the feasibility of accelerated diffusion tensor imaging (DTI) of skeletal muscle using echo planar imaging (EPI) applying simultaneous multislice excitation with a blipped controlled aliasing in parallel imaging results in higher acceleration unaliasing technique. MATERIALS AND METHODS After federal ethics board approval, the lower leg muscles of 8 healthy volunteers (mean [SD] age, 29.4 [2.9] years) were examined in a clinical 3-T magnetic ...

Demo: Distributed Real-Time Generative 3D Hand Tracking using Edge GPGPU Acceleration

DEFF Research Database (Denmark)

Qammaz, Ammar; Kosta, Sokol; Kyriazis, Nikolaos

2018-01-01

computations locally. The network connection takes the place of a GPGPU accelerator and sharing resources with a larger workstation becomes the acceleration mechanism. The unique properties of a generative optimizer are examined and constitute a challenging use-case, since the requirement for real......This work demonstrates a real-time 3D hand tracking application that runs via computation offloading. The proposed framework enables the application to run on low-end mobile devices such as laptops and tablets, despite the fact that they lack the sufficient hardware to perform the required...

Automated curved planar reformation of 3D spine images

International Nuclear Information System (INIS)

Vrtovec, Tomaz; Likar, Bostjan; Pernus, Franjo

2005-01-01

Traditional techniques for visualizing anatomical structures are based on planar cross-sections from volume images, such as images obtained by computed tomography (CT) or magnetic resonance imaging (MRI). However, planar cross-sections taken in the coordinate system of the 3D image often do not provide sufficient or qualitative enough diagnostic information, because planar cross-sections cannot follow curved anatomical structures (e.g. arteries, colon, spine, etc). Therefore, not all of the important details can be shown simultaneously in any planar cross-section. To overcome this problem, reformatted images in the coordinate system of the inspected structure must be created. This operation is usually referred to as curved planar reformation (CPR). In this paper we propose an automated method for CPR of 3D spine images, which is based on the image transformation from the standard image-based to a novel spine-based coordinate system. The axes of the proposed spine-based coordinate system are determined on the curve that represents the vertebral column, and the rotation of the vertebrae around the spine curve, both of which are described by polynomial models. The optimal polynomial parameters are obtained in an image analysis based optimization framework. The proposed method was qualitatively and quantitatively evaluated on five CT spine images. The method performed well on both normal and pathological cases and was consistent with manually obtained ground truth data. The proposed spine-based CPR benefits from reduced structural complexity in favour of improved feature perception of the spine. The reformatted images are diagnostically valuable and enable easier navigation, manipulation and orientation in 3D space. Moreover, reformatted images may prove useful for segmentation and other image analysis tasks

Augmented Reality Imaging System: 3D Viewing of a Breast Cancer.

Science.gov (United States)

Douglas, David B; Boone, John M; Petricoin, Emanuel; Liotta, Lance; Wilson, Eugene

2016-01-01

To display images of breast cancer from a dedicated breast CT using Depth 3-Dimensional (D3D) augmented reality. A case of breast cancer imaged using contrast-enhanced breast CT (Computed Tomography) was viewed with the augmented reality imaging, which uses a head display unit (HDU) and joystick control interface. The augmented reality system demonstrated 3D viewing of the breast mass with head position tracking, stereoscopic depth perception, focal point convergence and the use of a 3D cursor and joy-stick enabled fly through with visualization of the spiculations extending from the breast cancer. The augmented reality system provided 3D visualization of the breast cancer with depth perception and visualization of the mass's spiculations. The augmented reality system should be further researched to determine the utility in clinical practice.

Subsurface Profile Mapping using 3-D Compressive Wave Imaging

Directory of Open Access Journals (Sweden)

Hazreek Z A M

2017-01-01

Full Text Available Geotechnical site investigation related to subsurface profile mapping was commonly performed to provide valuable data for design and construction stage based on conventional drilling techniques. From past experience, drilling techniques particularly using borehole method suffer from limitations related to expensive, time consuming and limited data coverage. Hence, this study performs subsurface profile mapping using 3-D compressive wave imaging in order to minimize those conventional method constraints. Field measurement and data analysis of compressive wave (p-wave, vp was performed using seismic refraction survey (ABEM Terraloc MK 8, 7 kg of sledgehammer and 24 units of vertical geophone and OPTIM (SeisOpt@Picker & SeisOpt@2D software respectively. Then, 3-D compressive wave distribution of subsurface studied was obtained using analysis of SURFER software. Based on 3-D compressive wave image analyzed, it was found that subsurface profile studied consist of three main layers representing top soil (vp = 376 – 600 m/s, weathered material (vp = 900 – 2600 m/s and bedrock (vp > 3000 m/s. Thickness of each layer was varied from 0 – 2 m (first layer, 2 – 20 m (second layer and 20 m and over (third layer. Moreover, groundwater (vp = 1400 – 1600 m/s starts to be detected at 2.0 m depth from ground surface. This study has demonstrated that geotechnical site investigation data related to subsurface profiling was applicable to be obtained using 3-D compressive wave imaging. Furthermore, 3-D compressive wave imaging was performed based on non destructive principle in ground exploration thus consider economic, less time, large data coverage and sustainable to our environment.

3D-2D registration in endovascular image-guided surgery: evaluation of state-of-the-art methods on cerebral angiograms.

Science.gov (United States)

Mitrović, Uroš; Likar, Boštjan; Pernuš, Franjo; Špiclin, Žiga

2018-02-01

Image guidance for minimally invasive surgery is based on spatial co-registration and fusion of 3D pre-interventional images and treatment plans with the 2D live intra-interventional images. The spatial co-registration or 3D-2D registration is the key enabling technology; however, the performance of state-of-the-art automated methods is rather unclear as they have not been assessed under the same test conditions. Herein we perform a quantitative and comparative evaluation of ten state-of-the-art methods for 3D-2D registration on a public dataset of clinical angiograms. Image database consisted of 3D and 2D angiograms of 25 patients undergoing treatment for cerebral aneurysms or arteriovenous malformations. On each of the datasets, highly accurate "gold-standard" registrations of 3D and 2D images were established based on patient-attached fiducial markers. The database was used to rigorously evaluate ten state-of-the-art 3D-2D registration methods, namely two intensity-, two gradient-, three feature-based and three hybrid methods, both for registration of 3D pre-interventional image to monoplane or biplane 2D images. Intensity-based methods were most accurate in all tests (0.3 mm). One of the hybrid methods was most robust with 98.75% of successful registrations (SR) and capture range of 18 mm for registrations of 3D to biplane 2D angiograms. In general, registration accuracy was similar whether registration of 3D image was performed onto mono- or biplanar 2D images; however, the SR was substantially lower in case of 3D to monoplane 2D registration. Two feature-based and two hybrid methods had clinically feasible execution times in the order of a second. Performance of methods seems to fall below expectations in terms of robustness in case of registration of 3D to monoplane 2D images, while translation into clinical image guidance systems seems readily feasible for methods that perform registration of the 3D pre-interventional image onto biplanar intra

Coupling 2D/3D registration method and statistical model to perform 3D reconstruction from partial x-rays images data.

Science.gov (United States)

Cresson, T; Chav, R; Branchaud, D; Humbert, L; Godbout, B; Aubert, B; Skalli, W; De Guise, J A

2009-01-01

3D reconstructions of the spine from a frontal and sagittal radiographs is extremely challenging. The overlying features of soft tissues and air cavities interfere with image processing. It is also difficult to obtain information that is accurate enough to reconstruct complete 3D models. To overcome these problems, the proposed method efficiently combines the partial information contained in two images from a patient with a statistical 3D spine model generated from a database of scoliotic patients. The algorithm operates through two simultaneous iterating processes. The first one generates a personalized vertebra model using a 2D/3D registration process with bone boundaries extracted from radiographs, while the other one infers the position and the shape of other vertebrae from the current estimation of the registration process using a statistical 3D model. Experimental evaluations have shown good performances of the proposed approach in terms of accuracy and robustness when compared to CT-scan.

The impact of reconstruction method on the quantification of DaTSCAN images

Energy Technology Data Exchange (ETDEWEB)

Dickson, John C.; Erlandsson, Kjell; Hutton, Brian F. [UCLH NHS Foundation Trust and University College London, Institute of Nuclear Medicine, London (United Kingdom); Tossici-Bolt, Livia [Southampton University Hospitals NHS Trust, Department of Medical Physics, Southampton (United Kingdom); Sera, Terez [University of Szeged, Department of Nuclear Medicine and Euromedic Szeged, Szeged (Hungary); Varrone, Andrea [Psychiatry Section and Stockholm Brain Institute, Karolinska Institute, Department of Clinical Neuroscience, Stockholm (Sweden); Tatsch, Klaus [EANM/European Network of Excellence for Brain Imaging, Vienna (Austria)

2010-01-15

Reconstruction of DaTSCAN brain studies using OS-EM iterative reconstruction offers better image quality and more accurate quantification than filtered back-projection. However, reconstruction must proceed for a sufficient number of iterations to achieve stable and accurate data. This study assessed the impact of the number of iterations on the image quantification, comparing the results of the iterative reconstruction with filtered back-projection data. A striatal phantom filled with {sup 123}I using striatal to background ratios between 2:1 and 10:1 was imaged on five different gamma camera systems. Data from each system were reconstructed using OS-EM (which included depth-independent resolution recovery) with various combinations of iterations and subsets to achieve up to 200 EM-equivalent iterations and with filtered back-projection. Using volume of interest analysis, the relationships between image reconstruction strategy and quantification of striatal uptake were assessed. For phantom filling ratios of 5:1 or less, significant convergence of measured ratios occurred close to 100 EM-equivalent iterations, whereas for higher filling ratios, measured uptake ratios did not display a convergence pattern. Assessment of the count concentrations used to derive the measured uptake ratio showed that nonconvergence of low background count concentrations caused peaking in higher measured uptake ratios. Compared to filtered back-projection, OS-EM displayed larger uptake ratios because of the resolution recovery applied in the iterative algorithm. The number of EM-equivalent iterations used in OS-EM reconstruction influences the quantification of DaTSCAN studies because of incomplete convergence and possible bias in areas of low activity due to the nonnegativity constraint in OS-EM reconstruction. Nevertheless, OS-EM using 100 EM-equivalent iterations provides the best linear discriminatory measure to quantify the uptake in DaTSCAN studies. (orig.)

3D Membrane Imaging and Porosity Visualization

KAUST Repository

Sundaramoorthi, Ganesh; Hadwiger, Markus; Ben Romdhane, Mohamed; Behzad, Ali Reza; Madhavan, Poornima; Nunes, Suzana Pereira

2016-01-01

Ultrafiltration asymmetric porous membranes were imaged by two microscopy methods, which allow 3D reconstruction: Focused Ion Beam and Serial Block Face Scanning Electron Microscopy. A new algorithm was proposed to evaluate porosity and average pore

Dual optimization based prostate zonal segmentation in 3D MR images.

Science.gov (United States)

Qiu, Wu; Yuan, Jing; Ukwatta, Eranga; Sun, Yue; Rajchl, Martin; Fenster, Aaron

2014-05-01

Efficient and accurate segmentation of the prostate and two of its clinically meaningful sub-regions: the central gland (CG) and peripheral zone (PZ), from 3D MR images, is of great interest in image-guided prostate interventions and diagnosis of prostate cancer. In this work, a novel multi-region segmentation approach is proposed to simultaneously segment the prostate and its two major sub-regions from only a single 3D T2-weighted (T2w) MR image, which makes use of the prior spatial region consistency and incorporates a customized prostate appearance model into the segmentation task. The formulated challenging combinatorial optimization problem is solved by means of convex relaxation, for which a novel spatially continuous max-flow model is introduced as the dual optimization formulation to the studied convex relaxed optimization problem with region consistency constraints. The proposed continuous max-flow model derives an efficient duality-based algorithm that enjoys numerical advantages and can be easily implemented on GPUs. The proposed approach was validated using 18 3D prostate T2w MR images with a body-coil and 25 images with an endo-rectal coil. Experimental results demonstrate that the proposed method is capable of efficiently and accurately extracting both the prostate zones: CG and PZ, and the whole prostate gland from the input 3D prostate MR images, with a mean Dice similarity coefficient (DSC) of 89.3±3.2% for the whole gland (WG), 82.2±3.0% for the CG, and 69.1±6.9% for the PZ in 3D body-coil MR images; 89.2±3.3% for the WG, 83.0±2.4% for the CG, and 70.0±6.5% for the PZ in 3D endo-rectal coil MR images. In addition, the experiments of intra- and inter-observer variability introduced by user initialization indicate a good reproducibility of the proposed approach in terms of volume difference (VD) and coefficient-of-variation (CV) of DSC. Copyright © 2014 Elsevier B.V. All rights reserved.

3-D computations and measurements of accelerator magnets for the APS

International Nuclear Information System (INIS)

Turner, L.R.; Kim, S.H.; Kim, K.

1993-01-01

The Advanced Photon Source (APS), now under construction at Argonne National Laboratory (ANL), requires dipole, quadrupole, sextupole, and corrector magnets for each of its circular accelerator systems. Three-dimensional (3-D) field computations are needed to eliminate unwanted multipole fields from the ends of long quadrupole and dipole magnets and to guarantee that the flux levels in the poles of short magnets will not cause saturation. Measurements of the magnets show good agreement with the computations

Naked-eye 3D imaging employing a modified MIMO micro-ring conjugate mirrors

Science.gov (United States)

Youplao, P.; Pornsuwancharoen, N.; Amiri, I. S.; Thieu, V. N.; Yupapin, P.

2018-03-01

In this work, the use of a micro-conjugate mirror that can produce the 3D image incident probe and display is proposed. By using the proposed system together with the concept of naked-eye 3D imaging, a pixel and a large volume pixel of a 3D image can be created and displayed as naked-eye perception, which is valuable for the large volume naked-eye 3D imaging applications. In operation, a naked-eye 3D image that has a large pixel volume will be constructed by using the MIMO micro-ring conjugate mirror system. Thereafter, these 3D images, formed by the first micro-ring conjugate mirror system, can be transmitted through an optical link to a short distance away and reconstructed via the recovery conjugate mirror at the other end of the transmission. The image transmission is performed by the Fourier integral in MATLAB and compares to the Opti-wave program results. The Fourier convolution is also included for the large volume image transmission. The simulation is used for the manipulation, where the array of a micro-conjugate mirror system is designed and simulated for the MIMO system. The naked-eye 3D imaging is confirmed by the concept of the conjugate mirror in both the input and output images, in terms of the four-wave mixing (FWM), which is discussed and interpreted.

SU-E-J-100: Reconstruction of Prompt Gamma Ray Three Dimensional SPECT Image From Boron Neutron Capture Therapy(BNCT)

Energy Technology Data Exchange (ETDEWEB)

Yoon, D; Jung, J; Suh, T [The Catholic University of Korea, College of medicine, Department of biomedical engineering (Korea, Republic of)

2014-06-01

Purpose: Purpose of paper is to confirm the feasibility of acquisition of three dimensional single photon emission computed tomography (SPECT) image from boron neutron capture therapy (BNCT) using Monte Carlo simulation. Methods: In case of simulation, the pixelated SPECT detector, collimator and phantom were simulated using Monte Carlo n particle extended (MCNPX) simulation tool. A thermal neutron source (<1 eV) was used to react with the boron uptake region (BUR) in the phantom. Each geometry had a spherical pattern, and three different BURs (A, B and C region, density: 2.08 g/cm3) were located in the middle of the brain phantom. The data from 128 projections for each sorting process were used to achieve image reconstruction. The ordered subset expectation maximization (OSEM) reconstruction algorithm was used to obtain a tomographic image with eight subsets and five iterations. The receiver operating characteristic (ROC) curve analysis was used to evaluate the geometric accuracy of reconstructed image. Results: The OSEM image was compared with the original phantom pattern image. The area under the curve (AUC) was calculated as the gross area under each ROC curve. The three calculated AUC values were 0.738 (A region), 0.623 (B region), and 0.817 (C region). The differences between length of centers of two boron regions and distance of maximum count points were 0.3 cm, 1.6 cm and 1.4 cm. Conclusion: The possibility of extracting a 3D BNCT SPECT image was confirmed using the Monte Carlo simulation and OSEM algorithm. The prospects for obtaining an actual BNCT SPECT image were estimated from the quality of the simulated image and the simulation conditions. When multiple tumor region should be treated using the BNCT, a reasonable model to determine how many useful images can be obtained from the SPECT could be provided to the BNCT facilities. This research was supported by the Leading Foreign Research Institute Recruitment Program through the National Research