Search

In 2016, we published a method for processing dual-energy computed tomography (DECT) data called system-independent rho-e/Z-e ( $\rho _{\mathrm {e}}/Z_{\mathrm {e}}$ ) or SIRZ. Using data from multiple DECT systems and spectra, SIRZ estimated the electron density $\rho _{\mathrm {e}}$ and effective atomic number $Z_{\mathrm {e}}$ (based on published X-ray cross sections), for a set of known specimens. However, the decomposition process required complex spectral modeling of the DECT system, which made SIRZ difficult to implement and automate. This paper describes the subsequent work on “SIRZ-2” to simplify the spectral modeling, automate the process, and improve its range and versatility. The SIRZ-2 basis functions are more accurate for the X-ray energy range ( $Z_{\mathrm {e}}$ of 6–20) were scanned with multiple spectral pairs (up to 200 keV) on seven different DECT systems, including a commercial airport luggage scanner. For these varied tests, the average SIRZ-2 relative errors for $\rho _{\mathrm {e}}$ estimates were four times lower (0.7% compared to 3.0% for SIRZ), while average $Z_{\mathrm {e}}$ relative errors were comparably low for both methods ( $\rho _{\mathrm {e}}$ , $Z_{\mathrm {e}}$ ) that can be directly compared across DECT systems and over time.

In this paper, based on an invited talk at SORMA (May 2016), we present an overview of x-ray sources, detectors and system configurations for non-intrusive inspection (NII) of cargo containers. Our emphasis is on dual-energy x-ray NII for detecting high-atomic-number ( $\text {Z}\geq 72$ ) materials such as tungsten shielding and special nuclear materials (SNM). Standard single-energy (MeV and above) x-rays needed to penetrate and image cargo provide little SNM contrast, whereas dual-energy x-ray NII is demonstrated as a way to improve the selectivity of materials with Z $\text {Z}\geq 72$ . For two possible dual-energy x-ray source technologies—polyenergetic dual-energy (PDE) and quasi-monoenergetic x-ray sources (QMXS)—we investigate their trade-offs and future prospects using experimental and simulated results. The reduced scatter and larger separation of low- and high-energy photons provided by QMXS offers improved high-Z material contrast, but practical considerations such as flux and pulse rate need to be solved before making a deployable system. Straight-ray simulations show factor of four increases in contrast for QMXS over PDE scans of tin (Z=50) and iron (Z=26) relative to a uranium plate (Z=92) behind 20 cm of iron simulated cargo.

We present a new decomposition approach for dual-energy computed tomography (DECT) called SIRZ that provides precise and accurate material description, independent of the scanner, over diagnostic energy ranges (30 to 200 keV). System independence is achieved by explicitly including a scanner-specific spectral description in the decomposition method, and a new X-ray-relevant feature space. The feature space consists of electron density, ${\rho _{\rm e}}$ , and a new effective atomic number, ${{\rm Z}_{\rm e}}$ , which is based on published X-ray cross sections. Reference materials are used in conjunction with the system spectral response so that additional beam-hardening correction is not necessary. The technique is tested against other methods on DECT data of known specimens scanned by diverse spectra and systems. Uncertainties in accuracy and precision are less than 3% and 2% respectively for the ( ${\rho _{\rm e}}$ , ${{\rm Z}_{\rm e}}$ ) results compared to prior methods that are inaccurate and imprecise (over 9%).

Recent advances in cadmium telluride (CdTe) energy-discriminating pixelated detectors have enabled the possibility of Multi-Spectral X-ray Computed Tomography (MSXCT) to incorporate spectroscopic information into CT. MultiX ME 100 V2 is a CdTe-based spectroscopic x-ray detector array capable of recording energies from 20 to 160 keV in 1.1 keV energy bin increments. Hardware and software have been designed to perform radiographic and computed tomography tasks with this spectroscopic detector. Energy calibration is examined using the end-point energy of a bremsstrahlung spectrum and radioisotope spectral lines. When measuring the spectrum from Am-241 across 500 detector elements, the standard deviation of the peak-location and FWHM measurements are ± 0.4 and ± 0.6 keV, respectively. As these values are within the energy bin size (1.1 keV), detector elements are consistent with each other. The count rate is characterized, using a nonparalyzable model with a dead time of 64 ± 5 ns. This is consistent with the manufacturer’s quoted per detector-element linear-deviation at 2 Mpps (million photons per sec) of 8.9 % (typical) and 12 % (max). When comparing measured and simulated spectra, a low-energy tail is visible in the measured data due to the spectral response of the detector. If no valid photon detections are expected in the low-energy tail, then a background subtraction may be applied to allow for a possible first-order correction. If photons are expected in the low-energy tail, a detailed model must be implemented. A radiograph of an aluminum step wedge with a maximum height of 20 mm shows an underestimation of attenuation by about 10 % at 60 keV. This error is due to partial energy deposition from higher energy (>60 keV) photons into a lower-energy (∼60 keV) bin, reducing the apparent attenuation. A radiograph of a polytetrafluoroethylene (PTFE) cylinder taken using a bremsstrahlung spectrum from an x-ray voltage of 100 kV filtered by 1.3 mm Cu is reconstructed using Abel inversion. As no counts are expected in the low energy tail, a first order background correction is applied to the spectrum. The measured linear attenuation coefficient (LAC) is within 10% of the expected value in the 60 to 100 keV range. Below 60 keV, low counts in the corrected spectrum and partial energy deposition from incident photons of energy greater than 60 keV into energy bins below 60 keV impact the LAC measurements. This report ends with a demonstration of the tomographic capability of the system. The quantitative understanding of the detector developed in this report will enable further study in evaluating the system for characterization of an object’s chemical make-up for industrial and security purposes.

In a recent journal article [IEEE Trans. Nuc. Sci., 63(1), 341-350, 2016], we introduced a novel method that decomposes dual-energy X-ray CT (DECT) data into electron density (ρe) and a new effective-atomic-number called Ze in pursuit of system-independent characterization of materials. The Ze of a material, unlike the traditional Zeff, is defined relative to the actual X-ray absorption properties of the constituent atoms in the material, which are based on published X-ray cross sections. Our DECT method, called SIRZ (System-Independent ρe, Ze), uses a set of well-known reference materials and an understanding of the system spectral response to produce accurate and precise estimates of the X-ray-relevant basis variables (ρe, Ze) regardless of scanner or spectra in diagnostic energy ranges (30 to 200 keV). Potentially, SIRZ can account for and correct spectral changes in a scanner over time and, because the system spectral response is included in the technique, additional beam-hardening correction is not needed. Results show accuracy (

The National Ignition Facility (NIF) at the Lawrence Livermore National Laboratory is the world's most energetic laser, providing a scientific research center to study inertial confinement fusion and matter at extreme energy densities and pressures. A target shot involves over 30 specialized diagnostics measuring critical x-ray, optical and nuclear phenomena to quantify ignition results for comparison with computational models. The Shot Analysis and Visualization System (SAVI) acquires and analyzes target diagnostic data for display within a time-budget of 30 min. Laser and target diagnostic data are automatically loaded into the NIF archive database through clustered software data collection agents. The SAVI Analysis Engine distributes signal and image processing tasks to a Linux cluster where computation is performed. Intermediate results are archived at each step of the analysis pipeline. Data is archived with metadata and pedigree. Experiment results are visualized through a web-based user interface in interactive dashboards tailored to single or multiple shot perspectives. The SAVI system integrates open-source software, commercial workflow tools, relational database and messaging technologies into a service-oriented and distributed software architecture that is highly parallel, scalable, and flexible. The architecture and functionality of the SAVI system will be presented along with examples.

Some diagnostics at the National Ignition Facility (NIF), including the Gamma Reaction History (GRH) diagnostic, require multiple channels of data to achieve the required dynamic range. These channels need to be stitched together into a single time series, and they may have non-uniform and redundant time samples. We chose to apply the popular cubic smoothing spline technique to our stitching problem because we needed a general non-parametric method. We adapted one of the algorithms in the literature, by Hutchinson and deHoog, to our needs. The modified algorithm and the resulting code perform a cubic smoothing spline fit to multiple data channels with redundant time samples and missing data points. The data channels can have different, time-varying, zero-mean white noise characteristics. The method we employ automatically determines an optimal smoothing level by minimizing the Generalized Cross Validation (GCV) score. In order to automatically validate the smoothing level selection, the Weighted Sum-Squared Residual (WSSR) and zero-mean tests are performed on the residuals. Further, confidence intervals, both analytical and Monte Carlo, are also calculated. In this paper, we describe the derivation of our cubic smoothing spline algorithm. We outline the algorithm and test it with simulated and experimental data.

X-ray computerized tomography (CAT or CT) is a sophisticated imaging technique that provides cross-sectional views of materials, components and assemblies for industrial nondestructive evaluation (NDE). We have studied the feasibility of using CT as an inspection tool for reinforced concrete and the use of multi-energy, linear, attenuation techniques to deduce variations in density (ρ) and/or atomic number (Z) that could be caused by varying the types of concrete mixes and/or compaction in the concrete itself. To perform this study, we designed and built a prototype medium-/high-energy (200- to 2000 keV) CT scanner — ZCAT — to image small concrete samples (± 30 cm in diameter and ± 75 cm in height) with a spatial resolution of about 2 mm. We used ZCAT to quantitatively inspect a 20 cm concrete cube with 1.27 cm diameter reinforcing bars (rebars) and to measure p and/or Z variations in a 20 cm diameter concrete cylinder. We describe the ZCAT scanner design, some of its physical limitations and the data-acquisition parameters used in our study. Our results and those of others [1,2] show that CT can be used to inspect reinforced concrete and to distinguish material p and/or Z variations within concrete.

An efficient method for accurately calculating the center-of-rotation, or projection center, for parallel computed tomography projection data, or sinograms, is described. This method uses all the data in the sinogram to estimate the center by a least-squares technique and requires no previous calibration scan. The method also finds the object's center-of-mass without reconstructing its image. Since the method uses the measured data, it is sensitive to noise in the measurements, but that sensitivity is relatively small compared to other techniques. Examples of its use on simulated and actual data are included. For fan-beam data over 360 degrees , two related methods are described to find the center in the presence or absence of a midline offset. >

x-Ray computed axial tomography (CT) provides cross-sectional views of materials, components, and assemblies for industrial non-destructive evaluation. We have applied CT imaging to quantitatively measure the 3-D distribution ogf x-ray attenuation at reasonably high resolutions. In our industrial x-ray CT-studies, we have centered on two technical approaches: a first-generation translate/rotate CT system that consist of well-collimated (∼ 0.55 mm) photon source detector, and a third-generation scanner that uses a fluoroscopy detector.

Ultrasound imaging is widely used in medicine because of its benign characteristics and real-time capabilities. Physics theory suggests that the application of tomographic techniques may allow ultrasound imaging to reach its full potential as a diagnostic tool allowing it to compete with other tomographic modalities such as x-ray computer tomography, and MRI. This paper describes the construction and use of a prototype tomographic scanner and reports on the feasibility of implementing tomographic theory in practice and the potential of ultrasound (US) tomography in diagnostic imaging. Data were collected with the prototype by scanning two types of phantoms and a cadaveric breast. A specialized suite of algorithms was developed and utilized to construct images of reflectivity and sound speed from the phantom data. The basic results can be summarized as follows. (i) A fast, clinically relevant US tomography scanner can be built using existing technology. (ii) The spatial resolution, deduced from images of reflectivity, is 0.4 mm. The demonstrated 10 cm depth-of-field is superior to that of conventional ultrasound and the image contrast is improved through the reduction of speckle noise and overall lowering of the noise floor. (iii) Images of acoustic properties such as sound speed suggest that it is possible to measure variations in the sound speed of 5 m/s. An apparent correlation with x-ray attenuation suggests that the sound speed can be used to discriminate between various types of soft tissue. (iv) Ultrasound tomography has the potential to improve diagnostic imaging in relation to breast cancer detection.

A high-speed Radon transform computer (RTC) has been designed and simulated that operates on the principles of the parallel projection pipeline engine (PPPE) processor suggested by J.L.C. Sanz, et al. (1988). The computer consists of an array of identical processing elements with fast communication links to each of its neighbors and uses forward- and back-projection formulas tests linear interpolation. Simulation and experimental timing tests have been performed on both a systolic array of 64 transputer processors and on arrays of commercial DSP (digital signal processor) chips. A custom VLSI implementation is currently under design which is expected to perform forward and inverse transforms of 512 by 512 images at or near video rates. >

A class of hybrid domain-iterative algorithms for 2-D image reconstruction is described. The consistent iterative reconstruction-reprojection (CIRR) algorithm was designed and implemented along with other classical domain-iterative methods. The algorithms are tested on the hollow projections problem, where simple object models are assumed known. Results show that the CIRR algorithm is robust and accurate, and converges quickly. >

Our Computerized Ultrasound Risk Evaluation (CURE) system has been developed to the engineering prototype stage and generated unique data sets of both transmission and reflection ultrasound (US). This paper will help define the clinical underpinnings of the developmental process and interpret the imaging results from a similar perspective. The CURE project was designed to incorporate numerous diagnostic parameters to improve upon two major areas of early breast cancer detection. CURE may provide improved tissue characterization of breast masses and reliable detection of abnormal microcalcifications found in some breast cancers and ductal carcinoma in situ (DCIS). Current breast US is limited to mass evaluation, whereas mammography also detects and guides biopsy of malignant calcifications. Screening with CURE remains a distant goal, but improved follow-up of mammographic abnormalities may represent a feasible breakthrough.

A novel handheld time-domain array GPR antipersonnel mine detection system prototype has been developed. Using an offset paraboloidal reflector antenna to collimate rays form an ultra-wideband feed, the transmitted microwave impulse is concentrated forward, in front of the antenna structure. The resulting wave is a non-uniform plane wave over the portion of ground be investigated, and is incident at 45 degrees to normal. As such, much of the ground reflect wave is directed further forward, away from the operator, the reflector, and the receiving antennas, thereby reducing clutter. However, the wave transmitted into the ground, which interacts with the target, tends to have significant backscatter returning toward the receiving antennas. These receiving antennas are configured in a 2 by 2 array to provide spatial focusing in both along and cross-track directions. This is accomplished by measuring and comparing the backscattered signal at each receiver in the narrow time window between the times when the ground reflected wave passes the receiver and before this wave re-reflects from the reflector components. 2D FDTD simulation of this parabolic reflector transmitter indicates that it generates a beam with a non-uniform planar wavefront, which scatters form rough ground primarily in the forward direction. The wave transmitted into the ground is also planar, propagating at the angle of refraction, and scattering fairly isotropically from a small penetrable target. This system has been built and tested at LLNL, using a very narrow pulse shape. LLNL's Micro-Impulse Radar (MIR) and custom-built wideband antenna elements operate in the 1.5 to 5 GHz range. One particular advantage of using the MIR module is its low cost: an important feature for mine detectors used in developing countries. Preliminary measured data indicates that the surface clutter is indeed reduced relative to the target signal, and that small non-metallic anti-personnel mines can be reliably detected at burial depths as shallow as 1 inch in both dry.

Current ground penetrating radars (GPR) have been tested for land mine detection, but they have generally been costly and have poor performance. Comprehensive modeling and experimentation must be done to predict the electromagnetic (EM) signatures of mines to access the effect of clutter on the EM signature of the mine, and to understand the merit and limitations of using radar for various mine detection scenarios. This modeling can provide a basis for advanced radar design and detection techniques leading to superior performance. Lawrence Livermore National Laboratory (LLNL) has developed a radar technology that when combined with comprehensive modeling and detection methodologies could be the basis of an advanced mine detection system. Micropower Impulse Radar (MIR) technology exhibits a combination of properties, including wideband operation, extremely low power consumption, extremely small size and low cost, array configurability, and noise encoded pulse generation. LLNL is in the process of developing an 'optimal' processing algorithm to use with the MIR sensor. In this paper, we use classical numerical models to obtain the signature of mine-like targets and examine the effect of surface roughness on the reconstructed signals. These results are then qualitatively compared to experimental data.

Corrosion of rebar in concrete bridges causes subsurface cracks and is a major cause of structural degradation that necessitates repair or replacement. Early detection of corrosion effects can limit the location and extent of necessary repairs, while providing long-term information about the infrastructure status. Most current detection methods, however, are destructive of the road surface and require closing or restricting traffic while the tests are performed. A ground-penetrating radar imaging system has been designed and developed that will perform the nondestructive evaluation of road-bed cracking at traffic speeds; i.e., without the need to restrict traffic flow. The first-generation system (called the HERMES bridge inspector), consists of an offset-linear array of 64 impulse radar transceivers and associated electronics housed in a trailer. Computers in the trailer and in the towing vehicle control the data acquisition, processing, and display. Cross-road resolution is three centimeters at up to 30 cm in depth, while down-road resolution depends on speed; 3 cm below 20 mph up to 8 cm at 50 mph. A two-meter- wide path is inspected on each pass over the roadway. This paper, describes the design of this system, shows preliminary results, and lays out its deployment schedule.

The image reconstruction problem, also known as the inverse Radon transform, for x-ray computed tomography (CT) is found in numerous applications in medicine and industry. The most common algorithm used in these cases is filtered backprojection (FBP), which, while a simple procedure, is time-consuming for large images on any type of computational engine. Specially designed, dedicated parallel processors are commonly used in medical CT scanners, whose results are then passed to a graphics workstation for rendering and analysis. However, a fast direct FBP algorithm can be implemented on modern texture-mapping hardware in current high-end workstation platforms. This is done by casting the FBP algorithm as an image warping operation with summing. Texture- mapping hardware, such as that on the silicon Graphics Reality Engine, shows around 600 times speedup of backprojection over a CPU-based implementation (a 100 Mhz R4400 in our case). This technique has the further advantages of flexibility and rapid programming. In addition, the same hardware can be used for both image reconstruction and for volumetric rendering. Our technique can also be used to accelerate iterative reconstruction algorithms. The hardware architecture also allows more complex operations than straight-ray backprojection if they are required, including fan-beam, cone-beam, and curved ray paths, with little or no speed penalties.© (1994) COPYRIGHT SPIE--The International Society for Optical Engineering. Downloading of the abstract is permitted for personal use only.

Computed Tomography (CT) using penetrating radiation (x- or gamma-rays) can be used in a number of aircraft applications. This technique results in 3D volumetric attenuation data that is related to density and effective atomic number. CT is a transmission scanning method that must allow complete access to both sides of the object under inspection; the radiation source and detection systems must surround the object. This normally precludes the inspection of some large or planar (large aspect ratio) parts of the aircraft. However, we are pursuing recent limited-data techniques using object model information to obtain useful data from the partial information acquired. As illustrative examples, we describe how CT was instrumental in the analysis of particular aircraft components. These include fuselage panels, single crystal turbine blades, and aluminum-lithium composites.© (1993) COPYRIGHT SPIE--The International Society for Optical Engineering. Downloading of the abstract is permitted for personal use only.

A method and system for examining tissue are provided in which the tissue is maintained in a position so that it may be insonified with a plurality of pulsed spherical or cylindrical acoustic waves. The insonifying acoustic waves are scattered by the tissue so that scattered acoustic radiation including a mix of reflected and transmitted acoustic waves is received. A representation of a portion of the tissue is then derived from the received scattered acoustic radiation.

The amount of data generated by computed tomography (CT) scanners is enormous, making the image reconstruction operation slow, especially for 3-D and limited-data scans requiring iterative algorithms. The inverse Radon transform, commonly used for CT image reconstructions from projections, and the forward Radon transform are computationally burdensome for single-processor computer architectures. Fortunately, the forward Radon transform and the back projection operation (involved in the inverse Radon transform) are easily calculated using a parallel pipelined processor array. Using this array the processing time for the Radon transform and the back projection can be reduced dramatically. This paper describes a unified, pipelined architecture for an integrated circuit that computes both the forward Radon transform and the back projection operation at a 10 MHz data rate in a pipelined processor array. The trade-offs between computational complexity and reconstruction error of different interpolation schemes are presented along with an evaluation of the architecture's noise characteristics due to finite word lengths. The fully pipelined architecture is designed to reconstruct 1024 pixel by 1024 pixel images using up to 1024 projections over 180 degrees. The chip contains three pipelined data-paths, each five stages long, and uses a single multiplier.© (1990) COPYRIGHT SPIE--The International Society for Optical Engineering. Downloading of the abstract is permitted for personal use only.

X-ray computed axial tomography (CAT or CT) is an advanced imaging technique which provides three-dimensional nondestructive inspection and characterization of materials, components and assemblies. The Lawrence Livermore National Laboratory (LLNL) and the Mason Hanger-Silas Mason Co., Pantex Plant are cooperating to examine the use of CT Technology to inspect high-explosive (HE) pressings. The goals of this joint project are to study the HE pressing process using computed tomography scanners, to provide a foundation for the further development of industrial CT scanners, and to work with the private sector on the design and implementation of an automated CT inspection scanner to be used in the manufacturing process of high explosives at the Pantex Plant. The early detection of defective HE pressings will reduce the cost and improve the safety of the high-explosives production cycle. At present, all pressings are bulk density tested. In addition, samples are randomly selected from each pressing lot for destructive core sampling to obtain an estimate of internal density gradients. CT is a viable technology for a high-explosives inspection system. CT measures density variations, detects foreign inclusions and voids with a probability of detection as good as or better than present inspection techniques. In addition, CT canmore » be fully automated to allow for a 100% sample inspection. LLNL has the advantage of several different in-house built CT scanners to bring to bear on the inspection of high explosives. A preliminary study was performed to evaluate the LLNL CT scanners and prove the characterization principles associated with CT on small PBX9502 high-explosive samples. This paper summarizes the results obtained, and introduces the capabilities of LLNL's Nondestructive Evaluation Section CT scanners. 36 refs., 50 figs.« less

At Lawrence Livermore National Laboratory, there are experimental applications requiring digital signal acquisition as well as data reduction and analysis. A prototype Signal Acquisition/Processing Station (SAPS) has been constructed and is currently undergoing tests. The system employs an LSI-11/23 computer with Data Translation analog-to-digital hardware. SAPS is housed in a roll-around cart which has been designed to withstand most subtle EMI/RFI environments. A user-friendly menu allows a user to access powerful data acquisition packages with a minimum of training. The software architecture of SAPS involves two operating systems, each being transparent to the user. Since this is a general purpose workstation with several units being utilized, an emphasis on low cost, reliability, and maintenance was stressed during conception and design. The system is targeted fur mid-range frequency data acquisition; between a data logger and a transient digitizer.

The Systolic Processor with a Reconfigurable Interconnection Network of Transputers (SPRINT) [1] is a sixty-four-element multiprocessor developed at Lawrence Livermore National Laboratory to evaluate systolic algorithms and architectures experimentally. The processors are interconnected in a reconfigurable network which can emulate networks such as the two-dimensional mesh, the triangular mesh, the tree, and the shuffle-exchange network. New systolic algorithms and architectures are described which perform the Radon transform [8] and inverse Radon transform with efficiency arbitrarily close to 100%. High efficiency is possible with any connected network topology, even with low communication bandwidth. The results of the algorithms executed on the SPRINT compare closely with theory.