Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Iterative methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Task-based detectability in CT image reconstruction by filtered backprojection and penalized likelihood estimation

Purpose: Nonstationarity is an important aspect of imaging performance in CT and cone-beam CT (CBCT), especially for systems employing iterative reconstruction. This work presents a theoretical framework for both filtered-backprojection (FBP) and penalized-likelihood (PL) reconstruction that includes explicit descriptions of nonstationary noise, spatial resolution, and task-based detectability index. Potential utility of the model was demonstrated in the optimal selection of regularization parameters in PL reconstruction. Methods: Analytical models for local modulation transfer function (MTF) and noise-power spectrum (NPS) were investigated for both FBP and PL reconstruction, including explicit dependence on the object and spatial location. For FBP, a cascaded systems analysis framework was adapted to account for nonstationarity by separately calculating fluence and system gains for each ray passing through any given voxel. For PL, the point-spread function and covariance were derived using the implicit function theorem and first-order Taylor expansion according toFessler [“Mean and variance of implicitly defined biased estimators (such as penalized maximum likelihood): Applications to tomography,” IEEE Trans. Image Process. 5(3), 493–506 (1996)]. Detectability index was calculated for a variety of simple tasks. The model for PL was used in selecting the regularization strength parameter to optimize task-based performance, with both a constant and a spatially varying regularization map. Results: Theoretical models of FBP and PL were validated in 2D simulated fan-beam data and found to yield accurate predictions of local MTF and NPS as a function of the object and the spatial location. The NPS for both FBP and PL exhibit similar anisotropic nature depending on the pathlength (and therefore, the object and spatial location within the object) traversed by each ray, with the PL NPS experiencing greater smoothing along directions with higher noise. The MTF of FBP is isotropic and independent of location to a first order approximation, whereas the MTF of PL is anisotropic in a manner complementary to the NPS. Task-based detectability demonstrates dependence on the task, object, spatial location, and smoothing parameters. A spatially varying regularization “map” designed from locally optimal regularization can improve overall detectability beyond that achievable with the commonly used constant regularization parameter. Conclusions: Analytical models for task-based FBP and PL reconstruction are predictive of nonstationary noise and resolution characteristics, providing a valuable framework for understanding and optimizing system performance in CT and CBCT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

WE-G-18A-03: Cone Artifacts Correction in Iterative Cone Beam CT Reconstruction

Purpose: For iterative reconstruction (IR) in cone-beam CT (CBCT) imaging, data truncation along the superior-inferior (SI) direction causes severe cone artifacts in the reconstructed CBCT volume images. Not only does it reduce the effective SI coverage of the reconstructed volume, it also hinders the IR algorithm convergence. This is particular a problem for regularization based IR, where smoothing type regularization operations tend to propagate the artifacts to a large area. It is our purpose to develop a practical cone artifacts correction solution. Methods: We found it is the missing data residing in the truncated cone area that leads to inconsistency between the calculated forward projections and measured projections. We overcome this problem by using FDK type reconstruction to estimate the missing data and design weighting factors to compensate the inconsistency caused by the missing data. We validate the proposed methods in our multi-GPU low-dose CBCT reconstruction system on multiple patients' datasets. Results: Compared to the FDK reconstruction with full datasets, while IR is able to reconstruct CBCT images using a subset of projection data, the severe cone artifacts degrade overall image quality. For head-neck case under a full-fan mode, 13 out of 80 slices are contaminated. It is even more severe in pelvis case under half-fan mode, where 36 out of 80 slices are affected, leading to inferior soft-tissue delineation. By applying the proposed method, the cone artifacts are effectively corrected, with a mean intensity difference decreased from ∼497 HU to ∼39HU for those contaminated slices. Conclusion: A practical and effective solution for cone artifacts correction is proposed and validated in CBCT IR algorithm. This study is supported in part by NIH (1R01CA154747-01)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Alpha image reconstruction (AIR): A new iterative CT image reconstruction approach using voxel-wise alpha blending

Purpose: Iterative image reconstruction gains more and more interest in clinical routine, as it promises to reduce image noise (and thereby patient dose), to reduce artifacts, or to improve spatial resolution. Among vendors and researchers, however, there is no consensus of how to best achieve these aims. The general approach is to incorporatea priori knowledge into iterative image reconstruction, for example, by adding additional constraints to the cost function, which penalize variations between neighboring voxels. However, this approach to regularization in general poses a resolution noise trade-off because the stronger the regularization, and thus the noise reduction, the stronger the loss of spatial resolution and thus loss of anatomical detail. The authors propose a method which tries to improve this trade-off. The proposed reconstruction algorithm is called alpha image reconstruction (AIR). One starts with generating basis images, which emphasize certain desired image properties, like high resolution or low noise. The AIR algorithm reconstructs voxel-specific weighting coefficients that are applied to combine the basis images. By combining the desired properties of each basis image, one can generate an image with lower noise and maintained high contrast resolution thus improving the resolution noise trade-off. Methods: All simulations and reconstructions are performed in native fan-beam geometry. A water phantom with resolution bar patterns and low contrast disks is simulated. A filtered backprojection (FBP) reconstruction with a Ram-Lak kernel is used as a reference reconstruction. The results of AIR are compared against the FBP results and against a penalized weighted least squares reconstruction which uses total variation as regularization. The simulations are based on the geometry of the Siemens Somatom Definition Flash scanner. To quantitatively assess image quality, the authors analyze line profiles through resolution patterns to define a contrast factor for contrast-resolution plots. Furthermore, the authors calculate the contrast-to-noise ratio with the low contrast disks and the authors compare the agreement of the reconstructions with the ground truth by calculating the normalized cross-correlation and the root-mean-square deviation. To evaluate the clinical performance of the proposed method, the authors reconstruct patient data acquired with a Somatom Definition Flash dual source CT scanner (Siemens Healthcare, Forchheim, Germany). Results: The results of the simulation study show that among the compared algorithms AIR achieves the highest resolution and the highest agreement with the ground truth. Compared to the reference FBP reconstruction AIR is able to reduce the relative pixel noise by up to 50% and at the same time achieve a higher resolution by maintaining the edge information from the basis images. These results can be confirmed with the patient data. Conclusions: To evaluate the AIR algorithm simulated and measured patient data of a state-of-the-art clinical CT system were processed. It is shown, that generating CT images through the reconstruction of weighting coefficients has the potential to improve the resolution noise trade-off and thus to improve the dose usage in clinical CT.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Investigation of statistical iterative reconstruction for dedicated breast CT

Purpose: Dedicated breast CT has great potential for improving the detection and diagnosis of breast cancer. Statistical iterative reconstruction (SIR) in dedicated breast CT is a promising alternative to traditional filtered backprojection (FBP). One of the difficulties in using SIR is the presence of free parameters in the algorithm that control the appearance of the resulting image. These parameters require tuning in order to achieve high quality reconstructions. In this study, the authors investigated the penalized maximum likelihood (PML) method with two commonly used types of roughness penalty functions: hyperbolic potential and anisotropic total variation (TV) norm. Reconstructed images were compared with images obtained using standard FBP. Optimal parameters for PML with the hyperbolic prior are reported for the task of detecting microcalcifications embedded in breast tissue.Methods: Computer simulations were used to acquire projections in a half-cone beam geometry. The modeled setup describes a realistic breast CT benchtop system, with an x-ray spectra produced by a point source and an a-Si, CsI:Tl flat-panel detector. A voxelized anthropomorphic breast phantom with 280 μm microcalcification spheres embedded in it was used to model attenuation properties of the uncompressed woman's breast in a pendant position. The reconstruction of 3D images was performed using the separable paraboloidal surrogates algorithm with ordered subsets. Task performance was assessed with the ideal observer detectability index to determine optimal PML parameters.Results: The authors' findings suggest that there is a preferred range of values of the roughness penalty weight and the edge preservation threshold in the penalized objective function with the hyperbolic potential, which resulted in low noise images with high contrast microcalcifications preserved. In terms of numerical observer detectability index, the PML method with optimal parameters yielded substantially improved performance (by a factor of greater than 10) compared to FBP. The hyperbolic prior was also observed to be superior to the TV norm. A few of the best-performing parameter pairs for the PML method also demonstrated superior performance for various radiation doses. In fact, using PML with certain parameter values results in better images, acquired using 2 mGy dose, than FBP-reconstructed images acquired using 6 mGy dose.Conclusions: A range of optimal free parameters for the PML algorithm with hyperbolic and TV norm-based potentials is presented for the microcalcification detection task, in dedicated breast CT. The reported values can be used as starting values of the free parameters, when SIR techniques are used for image reconstruction. Significant improvement in image quality can be achieved by using PML with optimal combination of parameters, as compared to FBP. Importantly, these results suggest improved detection of microcalcifications can be obtained by using PML with lower radiation dose to the patient, than using FBP with higher dose.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

TH-AB-202-10: Quantifying the Accuracy and Precision of Six Degree-Of-Freedom Motion Estimation for Use in Real-Time Tumor Motion Monitoring During Radiotherapy

Purpose: Kilovoltage intrafraction monitoring (KIM) scheme has been successfully used to simultaneously monitor 3D tumor motion during radiotherapy. Recently, an iterative closest point (ICP) algorithm was implemented in KIM to also measure rotations about three axes, enabling real-time tracking of tumor motion in six degrees-of-freedom (DoF). This study aims to evaluate the accuracy of the six DoF motion estimates of KIM by comparing it with the corresponding motion (i) measured by the Calypso; and (ii) derived from kV/MV triangulation. Methods: (i) Various motions (static and dynamic) were applied to a CIRS phantom with three embedded electromagnetic transponders (Calypso Medical) using a 5D motion platform (HexaMotion) and a rotating treatment couch while both KIM and Calypso were used to concurrently track the phantom motion in six DoF. (ii) KIM was also used to retrospectively estimate six DoF motion from continuous sets of kV projections of a prostate, implanted with three gold fiducial markers (2 patients with 80 fractions in total), acquired during the treatment. Corresponding motion was obtained from kV/MV triangulation using a closed form least squares method based on three markers’ positions. Only the frames where all three markers were present were used in the analysis. The mean differences between the corresponding motion estimates were calculated for each DoF. Results: Experimental results showed that the mean of absolute differences in six DoF phantom motion measured by Calypso and KIM were within 1.1° and 0.7 mm. kV/MV triangulation derived six DoF prostate tumor better agreed with KIM estimated motion with the mean (s.d.) difference of up to 0.2° (1.36°) and 0.2 (0.25) mm for rotation and translation, respectively. Conclusion: These results suggest that KIM can provide an accurate six DoF intrafraction tumor during radiotherapy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2012

Temporal and spectral imaging with micro-CT

Purpose: Micro-CT is widely used for small animal imaging in preclinical studies of cardiopulmonary disease, but further development is needed to improve spatial resolution, temporal resolution, and material contrast. We present a technique for visualizing the changing distribution of iodine in the cardiac cycle with dual source micro-CT. Methods: The approach entails a retrospectively gated dual energy scan with optimized filters and voltages, and a series of computational operations to reconstruct the data. Projection interpolation and five-dimensional bilateral filtration (three spatial dimensions + time + energy) are used to reduce noise and artifacts associated with retrospective gating. We reconstruct separate volumes corresponding to different cardiac phases and apply a linear transformation to decompose these volumes into components representing concentrations of water and iodine. Since the resulting material images are still compromised by noise, we improve their quality in an iterative process that minimizes the discrepancy between the original acquired projections and the projections predicted by the reconstructed volumes. The values in the voxels of each of the reconstructed volumes represent the coefficients of linear combinations of basis functions over time and energy. We have implemented the reconstruction algorithm on a graphics processing unit (GPU) with CUDA. We tested the utility of the technique in simulations and applied the technique in an in vivo scan of a C57BL/6 mouse injected with blood pool contrast agent at a dose of 0.01 ml/g body weight. Postreconstruction, at each cardiac phase in the iodine images, we segmented the left ventricle and computed its volume. Using the maximum and minimum volumes in the left ventricle, we calculated the stroke volume, the ejection fraction, and the cardiac output. Results: Our proposed method produces five-dimensional volumetric images that distinguish different materials at different points in time, and can be used to segment regions containing iodinated blood and compute measures of cardiac function. Conclusions: We believe this combined spectral and temporal imaging technique will be useful for future studies of cardiopulmonary disease in small animals.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Material Assessment for ITER's Collective Thomson Scattering first mirror

The International Thermonuclear Energy Reactor (ITER) Collective Thomson Scattering (CTS) system is a diagnostic instrument that measures plasma density and velocity through Thomson scattering of microwave radiation. Some of the key components of the CTS are quasi-optical mirrors that are used to produce astigmatic beam patterns, which have impact on the strength and spatial resolution of the diagnostic signal. The mirrors are exposed to neutron radiation, which may alter the quality of the signal received. In this work, three different materials (molybdenum (Mo), stainless steel 316 (SS-316) and tungsten (W)) are considered for the first mirror of the CTS. The objective is to access which of the material studied are best suited for this mirror, considering different neutron radiation loads simulated scenarios defined by ITER, based on the resultant stresses and temperature distributions. For it, the neutron irradiation, and subsequent heat-load on the mirrors are simulated using the Monte Carlo N-Particle (MCNP) code. Based on the MCNP heat-load results, transient thermal-structural Finite Element Analysis (FEA) of the mirror over a 400 s discharge, with and without cooling on the rear side, are conducted using in commercial FEA software ANSYS. Results show that of the tested materials Mo and W are the most suitable material for this application. Even though, this study does not yet consider the variation of the material properties with temperature, it presents a quick initial satisfactory assessment that may be considered in subsequent and more complex analysis. (authors)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Sinogram restoration in computed tomography with an edge-preserving penalty

Purpose: With the goal of producing a less computationally intensive alternative to fully iterative penalized-likelihood image reconstruction, our group has explored the use of penalized-likelihood sinogram restoration for transmission tomography. Previously, we have exclusively used a quadratic penalty in our restoration objective function. However, a quadratic penalty does not excel at preserving edges while reducing noise. Here, we derive a restoration update equation for nonquadratic penalties. Additionally, we perform a feasibility study to extend our sinogram restoration method to a helical cone-beam geometry and clinical data. Methods: A restoration update equation for nonquadratic penalties is derived using separable parabolic surrogates (SPS). A method for calculating sinogram degradation coefficients for a helical cone-beam geometry is proposed. Using simulated data, sinogram restorations are performed using both a quadratic penalty and the edge-preserving Huber penalty. After sinogram restoration, Fourier-based analytical methods are used to obtain reconstructions, and resolution-noise trade-offs are investigated. For the fan-beam geometry, a comparison is made to image-domain SPS reconstruction using the Huber penalty. The effects of varying object size and contrast are also investigated. For the helical cone-beam geometry, we investigate the effect of helical pitch (axial movement/rotation). Huber-penalty sinogram restoration is performed on 3D clinical data, and the reconstructed images are compared to those generated with no restoration. Results: We find that by applying the edge-preserving Huber penalty to our sinogram restoration methods, the reconstructed image has a better resolution-noise relationship than an image produced using a quadratic penalty in the sinogram restoration. However, we find that this relatively straightforward approach to edge preservation in the sinogram domain is affected by the physical size of imaged objects in addition to the contrast across the edge. This presents some disadvantages of this method relative to image-domain edge-preserving methods, although the computational burden of the sinogram-domain approach is much lower. For a helical cone-beam geometry, we found applying sinogram restoration in 3D was reasonable and that pitch did not make a significant difference in the general effect of sinogram restoration. The application of Huber-penalty sinogram restoration to clinical data resulted in a reconstruction with less noise while retaining resolution. Conclusions: Sinogram restoration with the Huber penalty is able to provide better resolution-noise performance than restoration with a quadratic penalty. Additionally, sinogram restoration with the Huber penalty is feasible for helical cone-beam CT and can be applied to clinical data.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

SU-F-J-36: Comparison of Ball Bearing and Iso-Cube Phantoms for KV-MV Iso-Center Coincidence Check

Purpose: To compare two different quality assurance tools for kV-MV isocenter coincidence check. Methods: Ball-Bearing device (BBD) (Elekta) provided along with a CBCT equipped LINAC and isocenter cube phantom (CP) (a commercial product of Modus Medical Devices) are utilized to check the coincidence between the MV and kV beam isocenters. The microstepping meter of the BBD allows precision adjustment to better than 0.01mm in three directions. The BBD is aligned to the MV-isocenter with the lasers and followed by taking MV-images of at two collimator angles at each cardinal gantry, then its position adjusted. This process takes in iterative steps until 0.25mm tolerance is achieved. Four planar kV-images are taken at cardinal angles and deviations of BBD position from kVisocenter are evaluated. CP is positioned on the couch with the lasers at known offset from the isocenter. CBCT is utilized to set the cube to its kV isocenter with the precision of linac’s couch (1mm). MV-images are taken at various gantry, collimator, and couch angles and evaluated with the CP manufacturer provided QA software to determine coincidence of the 6mm steel ball at cube center with MV-isocenter. Three sets of measurements, one after another, were performed on the same day for each phantom on the same linac. Results: Assuming little variation between kV and MV isocenters in a short time, BBD measurements are more reproducible than CP. With comparable conditions both methods agree within 0.5 mm in each direction, while for BBD average standard deviation was 0.07mm and for iso-cube 0.3mm. Conclusion: Since BBD is more reliable and its results are more reproducible, it should be used during the monthly QA. Since CP is a more efficient device, it should be used for daily QA. A comparison study between the two devices should be periodically performed.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

TH-C-18A-01: Is Automatic Tube Current Modulation Still Necessary with Statistical Iterative Reconstruction?

Purpose: Automatic tube current modulation (TCM) has been widely used in modern multi-detector CT to reduce noise spatial nonuniformity and streaks to improve dose efficiency. With the advent of statistical iterative reconstruction (SIR), it is expected that the importance of TCM may diminish, since SIR incorporates statistical weighting factors to reduce the negative influence of photon-starved rays. The purpose of this work is to address the following questions: Does SIR offer the same benefits as TCM? If yes, are there still any clinical benefits to using TCM? Methods: An anthropomorphic CIRS chest phantom was scanned using a state-of-the-art clinical CT system equipped with an SIR engine (Veo™, GE Healthcare). The phantom was first scanned with TCM using a routine protocol and a low-dose (LD) protocol. It was then scanned without TCM using the same protocols. For each acquisition, both FBP and Veo reconstructions were performed. All scans were repeated 50 times to generate an image ensemble from which noise spatial nonuniformity (NSN) and streak artifact levels were quantified. Monte-Carlo experiments were performed to estimate skin dose. Results: For FBP, noise streaks were reduced by 4% using TCM for both routine and LD scans. NSN values were actually slightly higher with TCM (0.25) than without TCM (0.24) for both routine and LD scans. In contrast, for Veo, noise streaks became negligible (<1%) with or without TCM for both routine and LD scans, and the NSN was reduced to 0.10 (low dose) or 0.08 (routine). The overall skin dose was 2% lower at the shoulders and more uniformly distributed across the skin without TCM. Conclusion: SIR without TCM offers superior reduction in noise nonuniformity and streaks relative to FBP with TCM. For some clinical applications in which skin dose may be a concern, SIR without TCM may be a better option. K. Li, W. Zhao, D. Gomez-Cardona: Nothing to disclose; G.-H. Chen: Research funded, General Electric Company Research funded, Siemens AG Research funded, Varian Medical Systems, Research funded, Hologic, Inc.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

3D tomodosimetry using long scintillating fibers: A feasibility study

Purpose: 3D dosimetry is recognized as an ideal for patient-specific quality assurance (QA) of highly conformal radiotherapy treatments. However, existing 3D dosimeters are not straightforward to implement in the clinic, as their read-out procedure is often tedious and their accuracy, precision, and/or sample size exhibit limitations. The purpose of this work is to develop a 3D dosimeter based on the concept of tomodosimetry inside concentric cylindrical planes using long scintillating fibers for the QA of modern radiotherapy techniques such as intensity-modulated radiation therapy (IMRT) or intensity-modulated arc therapy (IMAT).Methods: Using a model-based simulation, scintillating fibers were modeled on three concentric cylindrical planes of radii 2.5, 5.0, and 7.5 cm, inside a 10 cm radius water-equivalent cylinder phantom. The phantom was set to rotate around its central axis, made parallel to the linac gantry axis of rotation. Light acquisitions were simulated using the calculated dose from the treatment planning software and reconstructed in each cylindrical plane at a resolution of 1 mm{sup 2} using a total-variation minimization iterative reconstruction algorithm. The 3D dose was then interpolated from the reconstructed cylindrical plane doses at a resolution of 1 mm{sup 3}. Different scintillating fiber patterns were compared by varying the angle of each fiber in its cylindrical plane and introducing a light-tight cut in each fiber. The precision of the reconstructed cylindrical dose distribution was evaluated using a Poisson modeling of the acquired light signals and the accuracy of the interpolated 3D dose was evaluated using an IMRT clinical plan for a prostate case.Results: Straight scintillating fiber patterns with light-tight cuts were the most accurate in cylindrical dose reconstruction, showing less than 0.5 mm distance-to-agreement in dose gradients and a mean local dose difference of less than 0.2% in the high dose region for a 10 × 10 cm{sup 2} field. The precision attained with this fiber configuration was less than 0.9% in the high dose, low gradient region of an IMRT segment for light acquisitions of 0.1 MU over a 360 degree rotation of the cylinder phantom. 3D dose interpolation for the IMRT clinical plan yielded an overall dose difference with the reference input of less than 1%, except in high dose gradients.Conclusions: Using long scintillating fibers inside rotating, concentric cylindrical planes, the authors demonstrate that their tomodosimetry method has the potential for high resolution, precise, and accurate 3D dosimetry. Moreover, because of its water-equivalence and rotational symmetry, this design should find interesting application for both treatment QA and machine commissioning.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2012

Simplified treatment of exact resonance elastic scattering model in deterministic slowing down equation

Simplified treatment of resonance elastic scattering model considering thermal motion of heavy nuclides and the energy dependence of the resonance cross section was implemented into NJOY [1]. In order to solve deterministic slowing down equation considering the effect of up-scattering without iterative calculations, scattering kernel for heavy nuclides is pre-calculated by the formula derived by Ouisloumen and Sanchez [2], and neutron spectrum in up-scattering term is expressed by NR approximation. To check the verification of the simplified treatment, the treatment is applied to U-238 for the energy range from 4 eV to 200 eV. Calculated multi-group capture cross section of U-238 is greater than that of conventional method and the increase of the capture cross sections is remarkable as the temperature becomes high. Therefore Doppler coefficient calculated in UO{sub 2} fuel pin is calculated more negative value than that on conventional method. The impact on Doppler coefficient is equivalent to the results of exact treatment of resonance elastic scattering reported in previous studies [2-7]. The agreement supports the validation of the simplified treatment and therefore this treatment is applied for other heavy nuclide to evaluate the Doppler coefficient in MOX fuel. The result shows that the impact of considering thermal agitation in resonance scattering in Doppler coefficient comes mainly from U-238 and that of other heavy nuclides such as Pu-239, 240 etc. is not comparable in MOX fuel. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

The application of ICASPA to critical and safe mass calculations - 14586

ICASPA, an Improved Critical and Safe Parameters Algorithm, calculates the size of a fissile body for a given value of k-effective. ICASPA uses an iterative predictor/corrector algorithm in which the MONK9A code is used to derive a precise k-effective value for a given candidate body size; auxiliary equations are then used to adjust the body size. The process is repeated until MONK9A outputs the chosen target value of k-effective. EDF Energy operates eight nuclear power stations comprising one Pressurized Water Reactor (PWR) and seven Advanced Gas-cooled Reactor (AGR) stations. At these sites, the term 'fuel route' is used to encompass all of the ex-reactor arrangements for the handling and storage of fuel. The criticality safety cases for these fuel routes are based on the use of geometrically safe arrangements for the required quantities of fuel. The allowed quantities and arrangements of fuel are safely sub-critical, even under the contingency of accidental moderator ingress. The fuel route safety cases consider all credible potential fault conditions, including events such as dropped fuel. Any such events might change the geometry of the fuel elements, but it is difficult to calculate exact k-effective values for realistic models of dropped fuel. This difficulty is overcome by the use of an alternative approach in which limiting safe masses are calculated for bounding (worst-case) rearrangements of fuel rods and moderators, using geometries similar to those covered by criticality handbook data. Where required, the resulting safe masses are used to set safety case limits and conditions. For example, some fire extinguishants, such as water and ABC dry powder, are also effective moderators. Hence, in some fuel route areas, administrative criticality controls are used to either prohibit or limit the permitted quantities of these materials. This paper describes the work that has been carried out to adopt ICASPA as the approved calculation method for these safe mass calculations. The work has been carried out within the company quality assurance requirements for nuclear safety related software. It has involved the development of a configuration controlled version of ICASPA and the local verification and validation of the resulting code. The benefits from the adoption of ICASPA are also presented. These include accuracy improvements and, because the method uses the MONK9A code for k-effective calculations, the ability to use more sophisticated models for limiting case geometries. (authors)

42 ENGINEERING↗
DOE OSTI2014

Generic method for automatic bladder segmentation on cone beam CT using a patient-specific bladder shape model

Purpose: The aim of this study is to develop and validate a generic method for automatic bladder segmentation on cone beam computed tomography (CBCT), independent of gender and treatment position (prone or supine), using only pretreatment imaging data. Methods: Data of 20 patients, treated for tumors in the pelvic region with the entire bladder visible on CT and CBCT, were divided into four equally sized groups based on gender and treatment position. The full and empty bladder contour, that can be acquired with pretreatment CT imaging, were used to generate a patient-specific bladder shape model. This model was used to guide the segmentation process on CBCT. To obtain the bladder segmentation, the reference bladder contour was deformed iteratively by maximizing the cross-correlation between directional grey value gradients over the reference and CBCT bladder edge. To overcome incorrect segmentations caused by CBCT image artifacts, automatic adaptations were implemented. Moreover, locally incorrect segmentations could be adapted manually. After each adapted segmentation, the bladder shape model was expanded and new shape patterns were calculated for following segmentations. All available CBCTs were used to validate the segmentation algorithm. The bladder segmentations were validated by comparison with the manual delineations and the segmentation performance was quantified using the Dice similarity coefficient (DSC), surface distance error (SDE) and SD of contour-to-contour distances. Also, bladder volumes obtained by manual delineations and segmentations were compared using a Bland-Altman error analysis. Results: The mean DSC, mean SDE, and mean SD of contour-to-contour distances between segmentations and manual delineations were 0.87, 0.27 cm and 0.22 cm (female, prone), 0.85, 0.28 cm and 0.22 cm (female, supine), 0.89, 0.21 cm and 0.17 cm (male, supine) and 0.88, 0.23 cm and 0.17 cm (male, prone), respectively. Manual local adaptations improved the segmentation results significantly (p < 0.01) based on DSC (6.72%) and SD of contour-to-contour distances (0.08 cm) and decreased the 95% confidence intervals of the bladder volume differences. Moreover, expanding the shape model improved the segmentation results significantly (p < 0.01) based on DSC and SD of contour-to-contour distances. Conclusions: This patient-specific shape model based automatic bladder segmentation method on CBCT is accurate and generic. Our segmentation method only needs two pretreatment imaging data sets as prior knowledge, is independent of patient gender and patient treatment position and has the possibility to manually adapt the segmentation locally.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2012

Predicting target vessel location on robot-assisted coronary artery bypass graft using CT to ultrasound registration

Purpose: Although robot-assisted coronary artery bypass grafting (RA-CABG) has gained more acceptance worldwide, its success still depends on the surgeon's experience and expertise, and the conversion rate to full sternotomy is in the order of 15%-25%. One of the reasons for conversion is poor pre-operative planning, which is based solely on pre-operative computed tomography (CT) images. In this paper, the authors propose a technique to estimate the global peri-operative displacement of the heart and to predict the intra-operative target vessel location, validated via both an in vitro and a clinical study. Methods: As the peri-operative heart migration during RA-CABG has never been reported in the literatures, a simple in vitro validation study was conducted using a heart phantom. To mimic the clinical workflow, a pre-operative CT as well as peri-operative ultrasound images at three different stages in the procedure (Stage{sub 0}--following intubation; Stage{sub 1}--following lung deflation; and Stage{sub 2}--following thoracic insufflation) were acquired during the experiment. Following image acquisition, a rigid-body registration using iterative closest point algorithm with the robust estimator was employed to map the pre-operative stage to each of the peri-operative ones, to estimate the heart migration and predict the peri-operative target vessel location. Moreover, a clinical validation of this technique was conducted using offline patient data, where a Monte Carlo simulation was used to overcome the limitations arising due to the invisibility of the target vessel in the peri-operative ultrasound images. Results: For the in vitro study, the computed target registration error (TRE) at Stage{sub 0}, Stage{sub 1}, and Stage{sub 2} was 2.1, 3.3, and 2.6 mm, respectively. According to the offline clinical validation study, the maximum TRE at the left anterior descending (LAD) coronary artery was 4.1 mm at Stage{sub 0}, 5.1 mm at Stage{sub 1}, and 3.4 mm at Stage{sub 2}. Conclusions: The authors proposed a method to measure and validate peri-operative shifts of the heart during RA-CABG. In vitro and clinical validation studies were conducted and yielded a TRE in the order of 5 mm for all cases. As the desired clinical accuracy imposed by this procedure is on the order of one intercostal space (10-15 mm), our technique suits the clinical requirements. The authors therefore believe this technique has the potential to improve the pre-operative planning by updating peri-operative migration patterns of the heart and, consequently, will lead to reduced conversion to conventional open thoracic procedures.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2011

Testing and modeling of diffusion bonded prototype optical windows under ITER conditions

Glass-metal joints are a part of ITER optical diagnostics windows. These joints must be leak tight for the safety (presence of tritium in ITER) and to preserve the vacuum. They must also withstand the ITER environment: temperatures up to 220 deg.C and fast neutron fluxes of {approx}3.10{sup 9} n/cm{sup 2}.s. At the moment, little information is available about glass-metal joints suitable for ITER. Therefore, we performed mechanical and thermal tests on some prototypes of an aluminium diffusion bonded optical window. Finite element modeling with Abaqus code was used to understand the experimental results. The prototypes were helium leaking probably due to very tiny cracks in the interaction layer between the steel and the aluminium. However, they were all able to withstand a thermal cycling test up to 200 deg. C; no damage could be seen after the tests by visual inspection. The prototypes successfully passed push-out test with a 500 N load. During the destructive push-out tests the prototypes broke at a 6-12 kN load between the aluminium layer and the steel or the glass, depending on the surface quality of the glass. The microanalysis of the joints has also been performed. The finite element modeling of the push-out tests is in a reasonable agreement with the experiments. According to the model, the highest thermal stress is created in the aluminium layer. Thus, the aluminium joint seems to be the weakest part of the prototypes. If this layer is improved, it will probably make the prototype helium leak tight and as such, a good ITER window candidate. (authors)

42 ENGINEERING↗
DOE OSTI2015

A novel partial volume effects correction technique integrating deconvolution associated with denoising within an iterative PET image reconstruction

Purpose: Partial volume effect (PVE) plays an important role in both qualitative and quantitative PET image accuracy, especially for small structures. A previously proposed voxelwise PVE correction method applied on PET reconstructed images involves the use of Lucy–Richardson deconvolution incorporating wavelet-based denoising to limit the associated propagation of noise. The aim of this study is to incorporate the deconvolution, coupled with the denoising step, directly inside the iterative reconstruction process to further improve PVE correction. Methods: The list-mode ordered subset expectation maximization (OSEM) algorithm has been modified accordingly with the application of the Lucy–Richardson deconvolution algorithm to the current estimation of the image, at each reconstruction iteration. Acquisitions of the NEMA NU2-2001 IQ phantom were performed on a GE DRX PET/CT system to study the impact of incorporating the deconvolution inside the reconstruction [with and without the point spread function (PSF) model] in comparison to its application postreconstruction and to standard iterative reconstruction incorporating the PSF model. The impact of the denoising step was also evaluated. Images were semiquantitatively assessed by studying the trade-off between the intensity recovery and the noise level in the background estimated as relative standard deviation. Qualitative assessments of the developed methods were additionally performed on clinical cases. Results: Incorporating the deconvolution without denoising within the reconstruction achieved superior intensity recovery in comparison to both standard OSEM reconstruction integrating a PSF model and application of the deconvolution algorithm in a postreconstruction process. The addition of the denoising step permitted to limit the SNR degradation while preserving the intensity recovery. Conclusions: This study demonstrates the feasibility of incorporating the Lucy–Richardson deconvolution associated with a wavelet-based denoising in the reconstruction process to better correct for PVE. Future work includes further evaluations of the proposed method on clinical datasets and the use of improved PSF models.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Atomically resolved structural determination of graphene and its point defects via extrapolation assisted phase retrieval

Previously reported crystalline structures obtained by an iterative phase retrieval reconstruction of their diffraction patterns seem to be free from displaying any irregularities or defects in the lattice, which appears to be unrealistic. We demonstrate here that the structure of a nanocrystal including its atomic defects can unambiguously be recovered from its diffraction pattern alone by applying a direct phase retrieval procedure not relying on prior information of the object shape. Individual point defects in the atomic lattice are clearly apparent. Conventional phase retrieval routines assume isotropic scattering. We show that when dealing with electrons, the quantitatively correct transmission function of the sample cannot be retrieved due to anisotropic, strong forward scattering specific to electrons. We summarize the conditions for this phase retrieval method and show that the diffraction pattern can be extrapolated beyond the original record to even reveal formerly not visible Bragg peaks. Such extrapolated wave field pattern leads to enhanced spatial resolution in the reconstruction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗