Adam Wang

All Publications

• Markerless Tumor Tracking using Fast-kV Switching Dual Energy Fluoroscopy on a Benchtop System. Medical physics Haytmyradov, M., Mostafavi, H., Wang, A., Zhu, L., Surucu, M., Patel, R., Ganguly, A., Richmond, M., Cassetta, R., Harkenrider, M. M., Roeske, J. C. 2019

  Abstract

  PURPOSE: To evaluate markerless tumor tracking (MTT) using fast-kV switching dual energy (DE) fluoroscopy on a bench top system.METHODS: Fast-kV-switching DE fluoroscopy was implemented on a bench top which includes a turntable stand, flat panel detector, and x-ray tube. The customized generator firmware enables consecutive x-ray pulses that alternate between programmed high and low energies (e.g. 60 and 120 kVp) with a maximum frame rate of 15 Hz. In-house software was implemented to perform weighted DE subtraction of consecutive images to create an image sequence that removes bone and enhances soft tissues. The weighting factor was optimized based on gantry angle. To characterize this system, a phantom was used that simulates the chest anatomy and tumor motion in the lung. Five clinically relevant tumor sizes (5-25 mm diameter) were considered. The targets were programmed to move in the inferior-superior direction of the phantom, perpendicular to the x-ray beam, using a cos4 waveform to mimic respiratory motion. Target inserts were then tracked with MTT software using a template matching method. The optimal computed tomography (CT) slice thickness for template generation was also evaluated. Tracking success rate and accuracy were calculated in regions of the phantom where the target overlapped ribs vs. spine, to compare the performance of single energy (SE) and DE imaging methods.RESULTS: For the 5 mm target, a CT slice thickness of 0.75 mm resulted in the lowest tracking error. For the larger targets (≥ 10 mm) a CT slice thickness ≤ 2 mm resulted in comparable tracking errors for SE and DE images. Overall DE imaging improved MTT accuracy, relative to SE imaging, for all tumor targets in a rotational acquisition. Compared to SE, DE imaging increased tracking success rate of small target inserts (5 and 10 mm). For fast motion tracking, success rates improved from 23% to 64% and 74% to 90% for 5 and 10 mm targets inserts overlapping ribs, respectively. For slow moving targets success rates improved from 19% to 59% and 59% to 91% in 5 and 10 mm targets overlapping the ribs, respectively. Similar results were observed when the targets overlapped the spine. For larger targets (≥ 15 mm) tracking success rates were comparable using SE and DE imaging.CONCLUSION: This work presents the first results of MTT using fast-kV switching DE fluoroscopy. Using DE imaging has improved the tracking accuracy of MTT, especially for small targets. The results of this study will guide the future implementation of fast-kV switching DE imaging using the on-board imager of a linear accelerator. This article is protected by copyright. All rights reserved.

  View details for PubMedID 31059124

• A novel method for fast image simulation of flat panel detectors. Physics in medicine and biology Shi, M., Myronakis, M. E., Hu, Y., Jacobson, M. W., Lehmann, M., Fueglistaller, R., Huber, P., Baturin, P., Wang, A. S., Ferguson, D., Harris, T., Morf, D., Berbeco, R. I. 2019

  Abstract

  We have developed a novel method for fast image simulation of flat panel detectors, based on the photon energy deposition efficiency and the optical spread function (OSF). The proposed method, FastEPID, determines the photon detection using photon energy deposition and replaces particle transport within the detector with precalculated OSFs. The FastEPID results are validated against experimental measurement and conventional Monte Carlo simulation in terms of modulation transfer function (MTF), signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR), contrast, and relative difference of pixel value, obtained with a slanted slit image, Las Vegas phantom, and anthropomorphic pelvis phantom. Excellent agreement is observed between simulation and measurement in all cases. Without degrading image quality, the FastEPID method is capable of reducing simulation time up to a factor of 150. Multiple applications, such as imager design optimization for planar and volumetric imaging, are expected to benefit from the implementation of the FastEPID method.

  View details for DOI 10.1088/1361-6560/ab12aa

  View details for PubMedID 30901759

• Fast shading correction for cone-beam CT via partitioned tissue classification. Physics in medicine and biology Shi, L., Wang, A. S., Wei, J., Zhu, L. 2019

  Abstract

  The quantitative use of cone beam computed tomography (CBCT) in radiation therapy is limited by severe shading artifacts, even with system embedded correction. We recently proposed effective shading correction methods, using planning CT (pCT) as prior information to estimate low-frequency errors in either the projection domain or the image domain. In this work, we further improve the clinical practicality of our previous methods by removing the requirement of prior pCT images. Clinical CBCT images are typically composed of a limited number of tissue types. By utilizing the low-frequency characteristic of shading distribution, we first generate a "shading-free" template image by enforcing uniformity on CBCT voxels of the same tissue type via a technique named partitioned tissue classification. Only a small subset of voxels on the template image is used to generate sparse samples of shading errors. Local filtration, a Fourier transform based algorithm, is employed to efficiently process the sparse errors to compute a full-field distribution of shading errors for CBCT correction. We evaluate the method performance on an anthropomorphic pelvis phantom and 6 pelvis patients. The proposed method improves the image quality of CBCT on both phantom and patients to a level matching that of pCT. On phantom, the signal non-uniformity (SNU) is reduced from 12.11 to 3.11% and 8.40 to 2.21% on fat and muscle, respectively. The maximum CT number error is reduced from 70 to 10 HU and 73 to 11 HU on fat and muscle, respectively. On patients, the average SNU is reduced from 9.22% to 1.06% and 11.41% to 1.67% on fat and muscle, respectively. The maximum CT number error is reduced from 95 to 9 HU and 88 to 8 HU on fat and muscle, respectively. The typical processing time for one CBCT dataset is about 45 seconds on a standard PC. .

  View details for PubMedID 30721886

• Characterizing a novel scintillating glass for application to megavoltage cone-beam computed tomography. Medical physics Hu, Y., Shedlock, D., Wang, A., Rottmann, J., Baturin, P., Myronakis, M., Huber, P., Fueglistaller, R., Shi, M., Morf, D., Star-Lack, J., Berbeco, R. I. 2018

  Abstract

  PURPOSE: The purpose of the present study was to evaluate the performance of a prototype electric portal imaging device (EPID) with a high detective quantum efficiency (DQE) scintillator, LKH-5. Specifically, image quality in context of both planar and megavoltage (MV) cone-beam computed tomography (CBCT) is analyzed.METHODS: Planar image quality in terms of modulation transfer function (MTF), noise power spectrum (NPS), and detective quantum efficiency (DQE) are measured and compared to an existing EPID (AS-1200) using the 6 MV beamline for a Varian TrueBeam linac. Imager performance is contextualized for three-dimensional (3D), MV-CBCT performance by measuring imager lag and analyzing the expected degradation of the DQE as a function of dose. Finally, comparisons between reconstructed images of the Catphan phantom in terms of qualitative quality and signal-difference-to-noise ratio (SDNR) are made for 6 MV images using both conventional and LKH-5 EPIDs as well as for the kilovoltage (kV) on-board imager (OBI).RESULTS: Analysis of the NPS reveals linearity at all measured doses using the prototype LKH-5 detector. While the first zero of the MTF is much lower for the LKH-5 detector than the conventional EPID (0.6 cycles/mm vs. 1.6 cycles/mm), the normalized NPS (NNPS) multiplied by total quanta (qNNPS) of the LKH-5 detector is roughly a factor of 7-8 times lower, yielding a DQE(0) of approximately 8%. First, second, and third frame lag was measured at approximately 23%, 5%, and 1% respectively, although no noticeable image artifacts were apparent in reconstructed volumes. Analysis of low dose performance reveals that DQE(0) remains at 80% of its maximum value at a dose as low as 7.5x10-6 MU. For a 400 projection technique, this represents a total scan dose of 0.0030 MU, suggesting that if imaging doses are increased to a value typical of kV-CBCT scans (~2.7 cGy), the LKH-5 detector will retain quantum noise limited performance. Finally, comparing Catphan scans, the prototype detector exhibits much lower image noise than the conventional EPID, resulting in improved small object representation. Further, SDNR of H2 O and polystyrene cylinders improved from -1.95 and 2.94 to -15 and 18.7, respectively.CONCLUSIONS: Imaging performance of the prototype LKH-5 detector was measured and analyzed for both planar and 3D contexts. Improving noise transfer of the detector results in concurrent improvement of DQE(0). For 3D imaging, temporal characteristics were adequate for artifact-free performance and at relevant doses, the detector retained quantum-noise limited performance. Although quantitative MTF measurements suggest poorer resolution, small object representation of the prototype imager is qualitatively improved over the conventional detector due to the measured reduction in noise. This article is protected by copyright. All rights reserved.

  View details for PubMedID 30586163

• A novel phantom for characterization of dual energy imaging using an on-board imaging system. Physics in medicine and biology Haytmyradov, M., Patel, R., Mostafavi, H., Surucu, M., Wang, A. S., Harkenrider, M. M., Roeske, J. C. 2018

  Abstract

  Dual-energy (DE) imaging using an on-board imager (OBI) is being considered for real-time tumor tracking purposes. We describe here a custom phantom designed to optimize DE imaging parameters using the OBI of a commercial linear accelerator. The phantom was constructed of lung-, tissue- and bone-equivalent material slabs. Five simulated tumors located at two different depths were encased in the lung-equivalent materials. Two slabs with bone-equivalent material inserts were constructed to simulate ribs, which overlap the simulated tumors. DE bone suppression was performed using a weighted logarithmic subtraction based on an iterative method that minimized the contrast between simulated bone- and lung-equivalent materials. The phantom was subsequently used to evaluate different combinations of high-low energy pairs based on the signal-difference-to-noise ratio (SDNR) metric. The results show a strong correlation between tumor visibility and selected energy pairs, where higher energy separation leads to larger SDNR values. To evaluate the effect of image post-processing methods on tumor visibility, an anti-correlated noise reduction (ACNR) and adaptive kernel scatter correction methods were applied to subsequent DE images. Application of the ACNR technique approximately doubled the SDNR values, hence increasing tumor visibility, while scatter correction had little effect on SDNR values. This phantom allows for quick image acquisition and optimization of imaging parameters and weighting factors. Optimized DE imaging increases soft tissue visibility and may enhance automated lung tumor tracking allowing for real-time adaptive radiotherapy.

  View details for DOI 10.1088/1361-6560/aaf9dd

  View details for PubMedID 30566913

• Feasibility of closed-MLC tracking using high sensitivity and multi-layer electronic portal imagers PHYSICS IN MEDICINE AND BIOLOGY Hu, Y., Jacobson, M. W., Shi, M., Myronakis, M., Wang, A., Baturin, P., Huber, P., Fueglistaller, R., Morf, D., Star-Lack, J., Berbeco, R. 2018; 63 (23): 235030

  Abstract

  In radiation therapy, improvements in treatment conformality are often limited by movement of target tissue. To better treat the target, tumor tracking strategies involving beam's-eye-view (BEV) have been explored. However, localization surrogates like implanted fiducial markers may sometimes leave the field-of-view (FOV), as defined by the linear accelerator (LINAC) multi-leaf collimator (MLC). Radiation leakage through the MLC has been measured previously at approximately 1%-2%. High sensitivity prototype detectors imagers may improve the ability to visualize objects outside of the MLC FOV during treatment. The present study presents a proof-of-concept for tracking fiducial markers outside the MLC FOV by employing high sensitivity detectors using a high-efficiency, prototype scintillating glass called LKH-5 and also investigates the impact of multi-layer imager (MLI) architecture. It was found that by improving the detector efficiency, using either of these methods results in a reduction of dose required for fiducial marker visibility. Further, image correction by a rectangular median filter will improve fiducial marker representation in the MLC blocked images. Quantified by measuring the peak-to-sidelobe ratio (PSR) of the normalized cross correlation (NCC) between a template of the fiducial marker with the blocked MLC acquisition, visibility has been found at a threshold of roughly 5 for all configurations with a 3  ×  3 cm2 ROI. For typical gadolinium oxysulfide (GOS) detectors in single and simulated 4-layer configurations, the minimum dose required for visualization was 20 and 10 MU, respectively. For LKH-5 detectors in single and simulated 4-layer configurations, this minimum dose was reduced to 4 and 2 MU, respectively. With a 6 MV flattening filter free (FFF) beam dose rate of 1400 MU min-1, the maximum detector frame rate while maintaining fiducial visibility is approximately 12 fps for a 4-layer LKH-5 configuration.

  View details for DOI 10.1088/1361-6560/aaef60

  View details for Web of Science ID 000452388800004

  View details for PubMedID 30520416

• A fast, linear Boltzmann transport equation solver for computed tomography dose calculation (Acuros CTD). Medical physics Wang, A., Maslowski, A., Wareing, T., Star-Lack, J., Schmidt, T. G. 2018

  Abstract

  PURPOSE: To improve dose reporting of CT scans, patient-specific organ doses are highly desired. However, estimating the dose distribution in a fast and accurate manner remains challenging, despite advances in Monte Carlo methods. In this work, we present an alternative method that deterministically solves the linear Boltzmann transport equation (LBTE), which governs the behavior of x-ray photon transport through an object.METHODS: Our deterministic solver for CT dose (Acuros CTD) is based on the same approach used to estimate scatter in projection images of a CT scan (Acuros CTS). A deterministic method is used to compute photon fluence within the object, which is then converted to deposited energy by multiplying by known, material-specific conversion factors. To benchmark Acuros CTD, we used the AAPM Task Group 195 test for CT dose, which models an axial, fan-beam scan (10 mm thick beam) and calculates energy deposited in each organ of an anthropomorphic phantom. We also validated our own Monte Carlo implementation of Geant4 to use as a reference to compare Acuros against for other common geometries like an axial, cone-beam scan (160 mm thick beam) and a helical scan (40 mm thick beam with table motion for a pitch of 1).RESULTS: For the fan-beam scan, Acuros CTD accurately estimated organ dose, with a maximum error of 2.7% and RMSE of 1.4% when excluding organs with < 0.1% of the total energy deposited. The cone-beam and helical scans yielded similar levels of accuracy compared to Geant4. Increasing the number of source positions beyond 18 or decreasing the voxel size below 5*5*5 mm3 provided marginal improvement to the accuracy for the cone-beam scan but came at the expense of increased run time. Across the different scan geometries, run time of Acuros CTD ranged from 8 to 23 seconds.CONCLUSIONS: In this digital phantom study, a deterministic LBTE solver was capable of fast and accurate organ dose estimates. This article is protected by copyright. All rights reserved.

  View details for PubMedID 30471131

• A modified McKinnon-Bates (MKB) algorithm for improved 4D cone-beam computed tomography (CBCT) of the lung MEDICAL PHYSICS Star-Lack, J., Sun, M., Oelhafen, M., Berkus, T., Pavkovich, J., Brehm, M., Arheit, M., Paysan, P., Wang, A., Munro, P., Seghers, D., Carvalho, L., Verbakel, W. R. 2018; 45 (8): 3783–99

  Abstract

  Four-dimensional (4D) cone-beam computed tomography (CBCT) of the lung is an effective tool for motion management in radiotherapy but presents a challenge because of slow gantry rotation times. Sorting the individual projections by breathing phase and using an established technique such as Feldkamp-Davis-Kress (FDK) to generate corresponding phase-correlated (PC) three-dimensional (3D) images results in reconstructions (FDK-PC) that often contain severe streaking artifacts due to the sparse angular sampling distributions. These can be reduced by further slowing down the gantry at the expense of incurring unwanted increases in scan times and dose. A computationally efficient alternative is the McKinnon-Bates (MKB) reconstruction algorithm that has shown promise in reducing view aliasing-induced streaking but can produce ghosting artifacts that reduce contrast and impede the determination of motion trajectories. The purpose of this work was to identify and correct shortcomings in the MKB algorithm.In the general MKB approach, a time-averaged 3D prior image is first reconstructed. The prior is then forward-projected at the same angles as the original projection data creating time-averaged reprojections. These reprojections are subsequently subtracted from the original (unblurred) projections to create motion-encoded difference projections. The difference projections are reconstructed into PC difference images that are added to the well-sampled 3D prior to create the higher quality 4D image. The cause of the ghosting in the traditional 4D MKB images was studied and traced to motion-induced streaking in the prior that, when reprojected, has the undesirable effect of re-encoding for motion in what should be a purely time-averaged reprojection. A new method, designated as the modified McKinnon-Bates (mMKB) algorithm, was developed based on destreaking the prior. This was coupled with a postprocessing 4D bilateral filter for noise suppression and edge preservation (mMKBbf ). The algorithms were tested with the 4D XCAT phantom using four simulated scan times (57, 60, 120, 180 s) and with two in vivo thorax studies (acquisition time of 60 and 90 s). Contrast-to-noise ratios (CNRs) of the target lesions and overall visual quality of the images were assessed.Prior destreaking (mMKB algorithm) reduced ghosting artifacts and increased CNRs for all cases, with the biggest impacts seen in the end inhale (EI) and end exhale (EE) phases of the respiratory cycle. For the XCAT phantom, mMKB lesion CNR was 44% higher than the MKB lesion CNR and was 81% higher than the FDK-PC lesion CNR (EI and EE phases). The bilateral filter provided a further average CNR improvement of 87% with the highest increases associated with longer scan times. Across all phases and scan times, the maximum mMKBbf -to-FDK-PC CNR improvement was over 300%. In vivo results agreed with XCAT results. Significantly less ghosting was observed throughout the mMKB images including near the lesions-of-interest and the diaphragm allowing for, in one case, visualization of a small tumor with nearly 30 mm of motion. The maximum FDK-PC-to-MKBbf CNR improvement for Patient 1's lesion was 261% and for Patient 2's lesion was 318%.The 4D mMKB algorithm yields good quality coronal and sagittal images in the thorax that may provide sufficient information for patient verification.

  View details for PubMedID 29869784

• Physics considerations in MV-CBCT multi-layer imager design PHYSICS IN MEDICINE AND BIOLOGY Hu, Y., Fueglistaller, R., Myronakis, M., Rottmann, J., Wang, A., Shedlock, D., Morf, D., Baturin, P., Huber, P., Star-Lack, J., Berbeco, R. 2018; 63 (12): 125016

  Abstract

  Megavoltage (MV) cone-beam computed tomography (CBCT) using an electronic portal imaging (EPID) offers advantageous features, including 3D mapping, treatment beam registration, high-z artifact suppression, and direct radiation dose calculation. Adoption has been slowed by image quality limitations and concerns about imaging dose. Developments in imager design, including pixelated scintillators, structured phosphors, inexpensive scintillation materials, and multi-layer imager (MLI) architecture have been explored to improve EPID image quality and reduce imaging dose. The present study employs a hybrid Monte Carlo and linear systems model to determine the effect of detector design elements, such as multi-layer architecture and scintillation materials. We follow metrics of image quality including modulation transfer function (MTF) and noise power spectrum (NPS) from projection images to 3D reconstructions to in-plane slices and apply a task based figure-of-merit, the ideal observer signal-to-noise ratio (d') to determine the effect of detector design on object detectability. Generally, detectability was limited by detector noise performance. Deploying an MLI imager with a single scintillation material for all layers yields improvement in noise performance and d' linear with the number of layers. In general, improving x-ray absorption using thicker scintillators results in improved DQE(0). However, if light yield is low, performance will be affected by electronic noise at relatively high doses, resulting in rapid image quality degradation. Maximizing image quality in a heterogenous MLI detector (i.e. multiple different scintillation materials) is most affected by limiting total noise. However, while a second-order effect, maximizing total spatial resolution of the MLI detector is a balance between the intensity contribution of each layer against its individual MTF. So, while a thinner scintillator may yield a maximal individual-layer MTF, its quantum efficiency will be relatively low in comparison to a thicker scintillator and thus, intensity contribution may be insufficient to noticeably improve the total detector MTF.

  View details for DOI 10.1088/1361-6560/aac8c6

  View details for Web of Science ID 000435938800001

  View details for PubMedID 29846180

  View details for PubMedCentralID PMC6042214

• Multi-layer imager design for mega-voltage spectral imaging PHYSICS IN MEDICINE AND BIOLOGY Myronakis, M., Hu, Y., Fueglistaller, R., Wang, A., Baturin, P., Huber, P., Morf, D., Star-Lack, J., Berbeco, R. 2018; 63 (10): 105002

  Abstract

  The architecture of multi-layer imagers (MLIs) can be exploited to provide megavoltage spectral imaging (MVSPI) for specific imaging tasks. In the current work, we investigated bone suppression and gold fiducial contrast enhancement as two clinical tasks which could be improved with spectral imaging. A method based on analytical calculations that enables rapid investigation of MLI component materials and thicknesses was developed and validated against Monte Carlo computations. The figure of merit for task-specific imaging performance was the contrast-to-noise ratio (CNR) of the gold fiducial when the CNR of bone was equal to zero after a weighted subtraction of the signals obtained from each MLI layer. Results demonstrated a sharp increase in the CNR of gold when the build-up component or scintillation materials and thicknesses were modified. The potential for low-cost, prompt implementation of specific modifications (e.g. composition of the build-up component) could accelerate clinical translation of MVSPI.

  View details for DOI 10.1088/1361-6560/aabe21

  View details for Web of Science ID 000431949400002

  View details for PubMedID 29652670

  View details for PubMedCentralID PMC5991631

• Acuros CTS: A fast, linear Boltzmann transport equation solver for computed tomography scatter - Part I: Core algorithms and validation MEDICAL PHYSICS Maslowski, A., Wang, A., Sun, M., Wareing, T., Davis, I., Star-Lack, J. 2018; 45 (5): 1899–1913

  Abstract

  To describe Acuros® CTS, a new software tool for rapidly and accurately estimating scatter in x-ray projection images by deterministically solving the linear Boltzmann transport equation (LBTE).The LBTE describes the behavior of particles as they interact with an object across spatial, energy, and directional (propagation) domains. Acuros CTS deterministically solves the LBTE by modeling photon transport associated with an x-ray projection in three main steps: (a) Ray tracing photons from the x-ray source into the object where they experience their first scattering event and form scattering sources. (b) Propagating photons from their first scattering sources across the object in all directions to form second scattering sources, then repeating this process until all high-order scattering sources are computed using the source iteration method. (c) Ray-tracing photons from scattering sources within the object to the detector, accounting for the detector's energy and anti-scatter grid responses. To make this process computationally tractable, a combination of analytical and discrete methods is applied. The three domains are discretized using the Linear Discontinuous Finite Elements, Multigroup, and Discrete Ordinates methods, respectively, which confer the ability to maintain the accuracy of a continuous solution. Furthermore, through the implementation in CUDA, we sought to exploit the parallel computing capabilities of graphics processing units (GPUs) to achieve the speeds required for clinical utilization. Acuros CTS was validated against Geant4 Monte Carlo simulations using two digital phantoms: (a) a water phantom containing lung, air, and bone inserts (WLAB phantom) and (b) a pelvis phantom derived from a clinical CT dataset. For these studies, we modeled the TrueBeam® (Varian Medical Systems, Palo Alto, CA) kV imaging system with a source energy of 125 kVp. The imager comprised a 600 μm-thick Cesium Iodide (CsI) scintillator and a 10:1 one-dimensional anti-scatter grid. For the WLAB studies, the full-fan geometry without a bowtie filter was used (with and without the anti-scatter grid). For the pelvis phantom studies, a half-fan geometry with bowtie was used (with the anti-scatter grid). Scattered and primary photon fluences and energies deposited in the detector were recorded.The Acuros CTS and Monte Carlo results demonstrated excellent agreement. For the WLAB studies, the average percent difference between the Monte Carlo- and Acuros-generated scattered photon fluences at the face of the detector was -0.7%. After including the detector response, the average percent differences between the Monte Carlo- and Acuros-generated scatter fractions (SF) were -0.1% without the grid and 0.6% with the grid. For the digital pelvis simulation, the Monte Carlo- and Acuros-generated SFs agreed to within 0.1% on average, despite the scatter-to-primary ratios (SPRs) being as high as 5.5. The Acuros CTS computation time for each scatter image was ~1 s using a single GPU.Acuros CTS enables a fast and accurate calculation of scatter images by deterministically solving the LBTE thus offering a computationally attractive alternative to Monte Carlo methods. Part II describes the application of Acuros CTS to scatter correction of CBCT scans on the TrueBeam system.

  View details for PubMedID 29509970

  View details for PubMedCentralID PMC5948176

• Acuros CTS: A fast, linear Boltzmann transport equation solver for computed tomography scatter - Part II: System modeling, scatter correction, and optimization MEDICAL PHYSICS Wang, A., Maslowski, A., Messmer, P., Lehmann, M., Strzelecki, A., Yu, E., Paysan, P., Brehm, M., Munro, P., Star-Lack, J., Seghers, D. 2018; 45 (5): 1914–25

  Abstract

  To correct for scatter in kV cone-beam CT (CBCT) projection data on a clinical system using a new tool, Acuros® CTS, that estimates scatter images rapidly and accurately by deterministically solving the linear Boltzmann transport equation.Phantom and patient CBCT scans were acquired on TrueBeam® radiotherapy machines. A first-pass reconstruction was used to create water and bone density maps of the imaged object, which was updated to include a more accurate representation of the patient couch. The imaging system model accounted for the TrueBeam x-ray source (polychromatic spectrum, beam filtration, bowtie filter, and collimation hardware) and x-ray detection system (antiscatter grid, flat-panel imager). Acuros CTS then used the system and object models to estimate the scatter component of each projection image, which was subtracted from the measured projections. The corrected projections were then reconstructed to produce the final result. We examined the tradeoff between run time and accuracy using a Pareto optimization of key parameters, including the voxel size of the down-sampled object model, the number of pixels in the down-sampled detector, and the number of scatter images (angular down-sampling). All computations and reconstructions were performed on a research workstation containing two graphics processing units (GPUs). In addition, we established a method for selecting a subset of projections for which scatter images were calculated. The projections were selected to minimize interpolation errors in the remaining projections. Image quality improvement was assessed by measuring the accuracy of the reconstructed phantom and patient images.The Pareto optimization yielded a set of parameters with an average run time of 26 seconds for scatter correction while maintaining high accuracy of scatter estimation. This was achieved in part by means of optimizing the projection angles that were processed, thus favoring the use of more angles in the lateral (i.e., horizontal) direction and fewer angles in the AP direction. In a 40 cm solid water phantom reconstruction, nonuniformities were decreased from 217 HU without scatter correction to 51 HU with conventional (kernel-based) scatter correction to 17 HU with Acuros CTS-based scatter correction. In clinical pelvis scans, nonuniformities in the bladder were reduced from 85 HU with conventional scatter correction to 14 HU with Acuros CTS.Acuros CTS is a promising new tool for fast and accurate scatter correction for CBCT imaging. By carefully modeling the imaging chain and optimizing several parameters, we achieved high correction accuracies with computation times compatible with the clinical workflow. The improvement in image quality enables better soft-tissue visualization and potentially enables applications such as adaptive radiotherapy.

  View details for PubMedID 29509973

• Leveraging multi-layer imager detector design to improve low-dose performance for megavoltage cone-beam computed tomography PHYSICS IN MEDICINE AND BIOLOGY Hu, Y., Rottmann, J., Fueglistaller, R., Myronakis, M., Wang, A., Huber, P., Shedlock, D., Morf, D., Baturin, P., Star-Lack, J., Berbeco, R. 2018; 63 (3): 035022

  Abstract

  While megavoltage cone-beam computed tomography (CBCT) using an electronic portal imaging device (EPID) provides many advantages over kilovoltage (kV) CBCT, clinical adoption is limited by its high doses. Multi-layer imager (MLI) EPIDs increase DQE(0) while maintaining high resolution. However, even well-designed, high-performance MLIs suffer from increased electronic noise from each readout, degrading low-dose image quality. To improve low-dose performance, shift-and-bin addition (ShiBA) imaging is proposed, leveraging the unique architecture of the MLI. ShiBA combines hardware readout-binning and super-resolution concepts, reducing electronic noise while maintaining native image sampling. The imaging performance of full-resolution (FR); standard, aligned binned (BIN); and ShiBA images in terms of noise power spectrum (NPS), electronic NPS, modulation transfer function (MTF), and the ideal observer signal-to-noise ratio (SNR)-the detectability index (d')-are compared. The FR 4-layer readout of the prototype MLI exhibits an electronic NPS magnitude 6-times higher than a state-of-the-art single layer (SLI) EPID. Although the MLI is built on the same readout platform as the SLI, with each layer exhibiting equivalent electronic noise, the multi-stage readout of the MLI results in electronic noise 50% higher than simple summation. Electronic noise is mitigated in both BIN and ShiBA imaging, reducing its total by ~12 times. ShiBA further reduces the NPS, effectively upsampling the image, resulting in a multiplication by a sinc2 function. Normalized NPS show that neither ShiBA nor BIN otherwise affects image noise. The LSF shows that ShiBA removes the pixilation artifact of BIN images and mitigates the effect of detector shift, but does not quantifiably improve the MTF. ShiBA provides a pre-sampled representation of the images, mitigating phase dependence. Hardware binning strategies lower the quantum noise floor, with 2  ×  2 implementation reducing the dose at which DQE(0) degrades by 10% from 0.01 MU to 0.004 MU, representing 20% improvement in d'.

  View details for DOI 10.1088/1361-6560/aaa160

  View details for Web of Science ID 000423829300005

  View details for PubMedID 29235440

  View details for PubMedCentralID PMC5824638

• Investigation of combined kV/MV CBCT imaging with a high-DQE MV detector. Medical physics Lindsay, C., Bazalova-Carter, M., Wang, A., Shedlock, D., Wu, M., Newson, M., Xing, L., Ansbacher, W., Fahrig, R., Star-Lack, J. 2018

  Abstract

  Combined kV-MV cone-beam tomography (CBCT) imaging has been proposed for two potentially important image-guided radiotherapy applications: (a) scan time reduction (STR) and (b) metal artifact reduction (MAR). However, the feasibility of these techniques has been in question due to the low detective quantum efficiencies (DQEs) of commercially available electronic portal imagers (EPIDs). The goal of the work was to test whether a prototype high DQE MV detector can be used to generate acceptable quality pretreatment CBCT images at acceptable dose levels.6MV and 100 kVp projection data were acquired on a Truebeam system (Varian, Palo Alto, CA). The MV data were acquired using a prototype EPID containing two scintillators (a) a standard copper-gadolinium oxysulfide (Cu-GOS) screen having a zero-frequency DQE (DQE(0)) value of 1.4%, and (b) a prototype-focused cadmium tungstate (CWO) pixelated "strip" with a DQE(0) = 22%. The kV data were acquired using the standard onboard imager (DQE(0) = 70%). The angular spacing of the MV projections was 0.81° and the source output was 0.03 MU/projection while the kV projections were acquired with an angular spacing of 0.4° at 0.3 mAs/projection. Image quality was evaluated using (a) an 18-cm diameter electron density phantom (CIRS, Norfolk, VA) with nine contrast inserts and (b) the resolution section of the 20-cm diameter Catphan phantom (The Phantom Laboratory, Greenwich, NY). For the MAR studies, two opposing CIRS phantom inserts were replaced by steel rods. The reconstruction methods were based on combining MV and kV data into one sinogram. The MAR reconstruction utilized mostly kV raw data with only those rays corrupted by metal requiring replacement with MV data (total absorbed dose = 0.7 cGy). For the STR study, projections from partially overlapping 105°kV and MV acquisitions were combined to create a complete dataset that could have been acquired in 18 sec (absorbed dose = 2.5 cGy). MV-only (4.3 cGy) and kV-only (0.3 cGy) images were also reconstructed.The average signal-to-noise ratio (SNR) of the inserts in the MV-only CWO and GOS CIRS phantom images were 0.62× and 0.12× the SNR of the inserts in kV-only image, respectively. The limiting spatial resolutions in the MV-only GOS, MV-only CWO, and kV-only Catphan images were 3, 6, and 8 lp/cm, respectively. In the combined kV/CWO STR reconstruction, all contrast inserts were visible while only two were detectable in the kV/Cu-GOS image due to high levels of noise (average SNRs of kV/CWO and kV/GOS inserts were 0.97× and 0.18× the SNR of the kV-only inserts, respectively). In the kV-MV MAR reconstructions, streaking artifacts were substantially reduced with all inserts becoming clearly visible in the kV/CWO image while only two were visible in the kV/Cu-GOS image (average SNRs of the kV/CWO and kV/Cu-GOS CIRS with metal inserts were 0.94× and 0.35× the SNRs of the kV-only CIRS without metal inserts).We have demonstrated that a high-DQE MV detector can be applied to generating high-quality combined kV-MV images for SRT and MAR. Clinically acceptable doses were utilized.

  View details for DOI 10.1002/mp.13291

  View details for PubMedID 30428131

• Spectral imaging using clinical megavoltage beams and a novel multi-layer imager PHYSICS IN MEDICINE AND BIOLOGY Myronakis, M., Fueglistaller, R., Rottmann, J., Hu, Y., Wang, A., Baturin, P., Huber, P., Morf, D., Star-Lack, J., Berbeco, R. 2017; 62 (23): 9127–39

  Abstract

  We assess the feasibility of clinical megavoltage (MV) spectral imaging for material and bone separation with a novel multi-layer imager (MLI) prototype. The MLI provides higher detective quantum efficiency and lower noise than conventional electronic portal imagers. Simulated experiments were performed using a validated Monte Carlo model of the MLI to estimate energy absorption and energy separation between the MLI components. Material separation was evaluated experimentally using solid water and aluminum (Al), copper (Cu) and gold (Au) for 2.5 MV, 6 MV and 6 MV flattening filter free (FFF) clinical photon beams. An anthropomorphic phantom with implanted gold fiducials was utilized to further demonstrate bone/gold separation. Weighted subtraction imaging was employed for material and bone separation. The weighting factor (w) was iteratively estimated, with the optimal w value determined by minimization of the relative signal difference ([Formula: see text]) and signal-difference-to-noise ratio (SDNR) between material (or bone) and the background. Energy separation between layers of the MLI was mainly the result of beam hardening between components with an average energy separation between 34 and 47 keV depending on the x-ray beam energy. The minimum average energy of the detected spectrum in the phosphor layer was 123 keV in the top layer of the MLI with the 2.5 MV beam. The w values that minimized [Formula: see text] and SDNR for Al, Cu and Au were 0.89, 0.76 and 0.64 for 2.5 MV; for 6 MV FFF, w was 0.98, 0.93 and 0.77 respectively. Bone suppression in the anthropomorphic phantom resulted in improved visibility of the gold fiducials with the 2.5 MV beam. Optimization of the MLI design is required to achieve optimal separation at clinical MV beam energies.

  View details for DOI 10.1088/1361-6560/aa94f9

  View details for Web of Science ID 000415302000003

  View details for PubMedID 29053107

  View details for PubMedCentralID PMC5724525

• A novel method for quantification of beam's-eye-view tumor tracking performance MEDICAL PHYSICS Hu, Y., Myronakis, M., Rottmann, J., Wang, A., Morf, D., Shedlock, D., Baturin, P., Star-Lack, J., Berbeco, R. 2017; 44 (11): 5650–59

  Abstract

  In-treatment imaging using an electronic portal imaging device (EPID) can be used to confirm patient and tumor positioning. Real-time tumor tracking performance using current digital megavolt (MV) imagers is hindered by poor image quality. Novel EPID designs may help to improve quantum noise response, while also preserving the high spatial resolution of the current clinical detector. Recently investigated EPID design improvements include but are not limited to multi-layer imager (MLI) architecture, thick crystalline and amorphous scintillators, and phosphor pixilation and focusing. The goal of the present study was to provide a method of quantitating improvement in tracking performance as well as to reveal the physical underpinnings of detector design that impact tracking quality. The study employs a generalizable ideal observer methodology for the quantification of tumor tracking performance. The analysis is applied to study both the effect of increasing scintillator thickness on a standard, single-layer imager (SLI) design as well as the effect of MLI architecture on tracking performance.The present study uses the ideal observer signal-to-noise ratio (d') as a surrogate for tracking performance. We employ functions which model clinically relevant tasks and generalized frequency-domain imaging metrics to connect image quality with tumor tracking. A detection task for relevant Cartesian shapes (i.e., spheres and cylinders) was used to quantitate trackability of cases employing fiducial markers. Automated lung tumor tracking algorithms often leverage the differences in benign and malignant lung tissue textures. These types of algorithms (e.g., soft-tissue localization - STiL) were simulated by designing a discrimination task, which quantifies the differentiation of tissue textures, measured experimentally and fit as a power-law in trend (with exponent β) using a cohort of MV images of patient lungs. The modeled MTF and NPS were used to investigate the effect of scintillator thickness and MLI architecture on tumor tracking performance.Quantification of MV images of lung tissue as an inverse power-law with respect to frequency yields exponent values of β = 3.11 and 3.29 for benign and malignant tissues, respectively. Tracking performance with and without fiducials was found to be generally limited by quantum noise, a factor dominated by quantum detective efficiency (QDE). For generic SLI construction, increasing the scintillator thickness (gadolinium oxysulfide - GOS) from a standard 290 μm to 1720 μm reduces noise to about 10%. However, 81% of this reduction is appreciated between 290 and 1000 μm. In comparing MLI and SLI detectors of equivalent individual GOS layer thickness, the improvement in noise is equal to the number of layers in the detector (i.e., 4) with almost no difference in MTF. Further, improvement in tracking performance was slightly less than the square-root of the reduction in noise, approximately 84-90%. In comparing an MLI detector with an SLI with a GOS scintillator of equivalent total thickness, improvement in object detectability is approximately 34-39%.We have presented a novel method for quantification of tumor tracking quality and have applied this model to evaluate the performance of SLI and MLI EPID designs. We showed that improved tracking quality is primarily limited by improvements in NPS. When compared to very thick scintillator SLI, employing MLI architecture exhibits the same gains in QDE, but by mitigating the effect of optical Swank noise, results in more dramatic improvements in tracking performance.

  View details for PubMedID 28887836

  View details for PubMedCentralID PMC5689096

• A novel multilayer MV imager computational model for component optimization MEDICAL PHYSICS Myronakis, M., Star-Lack, J., Baturin, P., Rottmann, J., Morf, D., Wang, A., Hu, Y., Shedlock, D., Berbeco, R. I. 2017; 44 (8): 4213–22

  Abstract

  A novel Megavoltage (MV) multilayer imager (MLI) design featuring higher detective quantum efficiency and lower noise than current conventional MV imagers in clinical use has been recently reported. Optimization of the MLI design for multiple applications including tumor tracking, MV-CBCT and portal dosimetry requires a computational model that will provide insight into the physics processes that affect the overall and individual components' performance. The purpose of the current work was to develop and validate a comprehensive computational model that can be used for MLI optimization.The MLI model was built using the Geant4 Application for Tomographic Emission (GATE) application. The model includes x-ray and charged-particle interactions as well as the optical transfer within the phosphor. A first prototype MLI device featuring a stack of four detection layers was used for model validation. Each layer of the prototype contains a copper buildup plate, a phosphor screen and photodiode array. The model was validated against measured data of Modulation Transfer Function (MTF), Noise-Power Spectrum (NPS), and Detective Quantum Efficiency (DQE). MTF was computed using a slanted slit with 2.3° angle and 0.1 mm width. NPS was obtained using the autocorrelation function technique. DQE was calculated from MTF and NPS data. The comparison metrics between simulated and measured data were the Pearson's correlation coefficient (r) and the normalized root-mean-square error (NRMSE).Good agreement between measured and simulated MTF and NPS values was observed. Pearson's correlation coefficient for the combined signal from all layers of the MLI was equal to 0.9991 for MTF and 0.9992 for NPS; NRMSE was 0.0121 for MTF and 0.0194 for NPS. Similarly, the DQE correlation coefficient for the combined signal was 0.9888 and the NRMSE was 0.0686.A comprehensive model of the novel MLI design was developed using the GATE toolkit and validated against measured MTF, NPS, and DQE data acquired with a prototype device featuring four layers. This model will be used for further optimization of the imager components and configuration for clinical radiotherapy applications.

  View details for PubMedID 28555935

  View details for PubMedCentralID PMC5553708

• Accuracy of patient-specific organ dose estimates obtained using an automated image segmentation algorithm JOURNAL OF MEDICAL IMAGING Schmidt, T., Wang, A. S., Coradi, T., Haas, B., Star-Lack, J. 2016; 3 (4): 043502

  Abstract

  The overall goal of this work is to develop a rapid, accurate, and automated software tool to estimate patient-specific organ doses from computed tomography (CT) scans using simulations to generate dose maps combined with automated segmentation algorithms. This work quantified the accuracy of organ dose estimates obtained by an automated segmentation algorithm. We hypothesized that the autosegmentation algorithm is sufficiently accurate to provide organ dose estimates, since small errors delineating organ boundaries will have minimal effect when computing mean organ dose. A leave-one-out validation study of the automated algorithm was performed with 20 head-neck CT scans expertly segmented into nine regions. Mean organ doses of the automatically and expertly segmented regions were computed from Monte Carlo-generated dose maps and compared. The automated segmentation algorithm estimated the mean organ dose to be within 10% of the expert segmentation for regions other than the spinal canal, with the median error for each organ region below 2%. In the spinal canal region, the median error was [Formula: see text], with a maximum absolute error of 28% for the single-atlas approach and 11% for the multiatlas approach. The results demonstrate that the automated segmentation algorithm can provide accurate organ dose estimates despite some segmentation errors.

  View details for PubMedID 27921070

  View details for PubMedCentralID PMC5126039

• Non-local total-variation (NLTV) minimization combined with reweighted L1-norm for compressed sensing CT reconstruction PHYSICS IN MEDICINE AND BIOLOGY Kim, H., Chen, J., Wang, A., Chuang, C., Held, M., Pouliot, J. 2016; 61 (18): 6878–91

  Abstract

  The compressed sensing (CS) technique has been employed to reconstruct CT/CBCT images from fewer projections as it is designed to recover a sparse signal from highly under-sampled measurements. Since the CT image itself cannot be sparse, a variety of transforms were developed to make the image sufficiently sparse. The total-variation (TV) transform with local image gradient in L1-norm was adopted in most cases. This approach, however, which utilizes very local information and penalizes the weight at a constant rate regardless of different degrees of spatial gradient, may not produce qualified reconstructed images from noise-contaminated CT projection data. This work presents a new non-local operator of total-variation (NLTV) to overcome the deficits stated above by utilizing a more global search and non-uniform weight penalization in reconstruction. To further improve the reconstructed results, a reweighted L1-norm that approximates the ideal sparse signal recovery of the L0-norm is incorporated into the NLTV reconstruction with additional iterates. This study tested the proposed reconstruction method (reweighted NLTV) from under-sampled projections of 4 objects and 5 experiments (1 digital phantom with low and high noise scenarios, 1 pelvic CT, and 2 CBCT images). We assessed its performance against the conventional TV, NLTV and reweighted TV transforms in the tissue contrast, reconstruction accuracy, and imaging resolution by comparing contrast-noise-ratio (CNR), normalized root-mean square error (nRMSE), and profiles of the reconstructed images. Relative to the conventional NLTV, combining the reweighted L1-norm with NLTV further enhanced the CNRs by 2-4 times and improved reconstruction accuracy. Overall, except for the digital phantom with low noise simulation, our proposed algorithm produced the reconstructed image with the lowest nRMSEs and the highest CNRs for each experiment.

  View details for DOI 10.1088/0031-9155/61/18/6878

  View details for Web of Science ID 000384317800002

  View details for PubMedID 27589006

• Striped Ratio Grids for Scatter Estimation Hsieh, S. S., Wang, A. S., Star-Lack, J., Kontos, D., Flohr, T. G., Lo, J. Y. SPIE-INT SOC OPTICAL ENGINEERING. 2016

  View details for DOI 10.1117/12.2216896

  View details for Web of Science ID 000378352900018

• Accuracy of patient specific organ-dose estimates obtained using an automated image segmentation algorithm Schmidt, T., Wang, A., Coradi, T., Haas, B., Star-Lack, J., Kontos, D., Flohr, T. G., Lo, J. Y. SPIE-INT SOC OPTICAL ENGINEERING. 2016

  View details for DOI 10.1117/12.2217374

  View details for Web of Science ID 000378352900035

• A piecewise-focused high DQE detector for MV imaging MEDICAL PHYSICS Star-Lack, J., Shedlock, D., Swahn, D., Humber, D., Wang, A., Hirsh, H., Zentai, G., Sawkey, D., Kruger, I., Sun, M., Abel, E., Virshup, G., Shin, M., Fahrig, R. 2015; 42 (9): 5084-5099

  Abstract

  Electronic portal imagers (EPIDs) with high detective quantum efficiencies (DQEs) are sought to facilitate the use of the megavoltage (MV) radiotherapy treatment beam for image guidance. Potential advantages include high quality (treatment) beam's eye view imaging, and improved cone-beam computed tomography (CBCT) generating images with more accurate electron density maps with immunity to metal artifacts. One approach to increasing detector sensitivity is to couple a thick pixelated scintillator array to an active matrix flat panel imager (AMFPI) incorporating amorphous silicon thin film electronics. Cadmium tungstate (CWO) has many desirable scintillation properties including good light output, a high index of refraction, high optical transparency, and reasonable cost. However, due to the 0 1 0 cleave plane inherent in its crystalline structure, the difficulty of cutting and polishing CWO has, in part, limited its study relative to other scintillators such as cesium iodide and bismuth germanate (BGO). The goal of this work was to build and test a focused large-area pixelated "strip" CWO detector.A 361  ×  52 mm scintillator assembly that contained a total of 28 072 pixels was constructed. The assembly comprised seven subarrays, each 15 mm thick. Six of the subarrays were fabricated from CWO with a pixel pitch of 0.784 mm, while one array was constructed from BGO for comparison. Focusing was achieved by coupling the arrays to the Varian AS1000 AMFPI through a piecewise linear arc-shaped fiber optic plate. Simulation and experimental studies of modulation transfer function (MTF) and DQE were undertaken using a 6 MV beam, and comparisons were made between the performance of the pixelated strip assembly and the most common EPID configuration comprising a 1 mm-thick copper build-up plate attached to a 133 mg/cm(2) gadolinium oxysulfide scintillator screen (Cu-GOS). Projection radiographs and CBCT images of phantoms were acquired. The work also introduces the use of a lightweight edge phantom to generate MTF measurements at MV energies and shows its functional equivalence to the more cumbersome slit-based method.Measured and simulated DQE(0)'s of the pixelated CWO detector were 22% and 26%, respectively. The average measured and simulated ratios of CWO DQE(f) to Cu-GOS DQE(f) across the frequency range of 0.0-0.62 mm(-1) were 23 and 29, respectively. 2D and 3D imaging studies confirmed the large dose efficiency improvement and that focus was maintained across the field of view. In the CWO CBCT images, the measured spatial resolution was 7 lp/cm. The contrast-to-noise ratio was dramatically improved reflecting a 22 × sensitivity increase relative to Cu-GOS. The CWO scintillator material showed significantly higher stability and light yield than the BGO material.An efficient piecewise-focused pixelated strip scintillator for MV imaging is described that offers more than a 20-fold dose efficiency improvement over Cu-GOS.

  View details for DOI 10.1118/1.4927786

  View details for Web of Science ID 000360645000011

  View details for PubMedID 26328960

  View details for PubMedCentralID PMC4529442

• Accelerated statistical reconstruction for C-arm cone-beam CT using Nesterov's method MEDICAL PHYSICS Wang, A. S., Stayman, J., Otake, Y., Vogt, S., Kleinszig, G., Siewerdsen, J. H. 2015; 42 (5): 2699–2708

  Abstract

  To accelerate model-based iterative reconstruction (IR) methods for C-arm cone-beam CT (CBCT), thereby combining the benefits of improved image quality and/or reduced radiation dose with reconstruction times on the order of minutes rather than hours.The ordered-subsets, separable quadratic surrogates (OS-SQS) algorithm for solving the penalized-likelihood (PL) objective was modified to include Nesterov's method, which utilizes "momentum" from image updates of previous iterations to better inform the current iteration and provide significantly faster convergence. Reconstruction performance of an anthropomorphic head phantom was assessed on a benchtop CBCT system, followed by CBCT on a mobile C-arm, which provided typical levels of incomplete data, including lateral truncation. Additionally, a cadaveric torso that presented realistic soft-tissue and bony anatomy was imaged on the C-arm, and different projectors were assessed for reconstruction speed.Nesterov's method provided equivalent image quality to OS-SQS while reducing the reconstruction time by an order of magnitude (10.0 ×) by reducing the number of iterations required for convergence. The faster projectors were shown to produce similar levels of convergence as more accurate projectors and reduced the reconstruction time by another 5.3 ×. Despite the slower convergence of IR with truncated C-arm CBCT, comparison of PL reconstruction methods implemented on graphics processing units showed that reconstruction time was reduced from 106 min for the conventional OS-SQS method to as little as 2.0 min with Nesterov's method for a volumetric reconstruction of the head. In body imaging, reconstruction of the larger cadaveric torso was reduced from 159 min down to 3.3 min with Nesterov's method.The acceleration achieved through Nesterov's method combined with ordered subsets reduced IR times down to a few minutes. This improved compatibility with clinical workflow better enables broader adoption of IR in CBCT-guided procedures, with corresponding benefits in overcoming conventional limits of image quality at lower dose.

  View details for PubMedID 25979068

• Automatic Localization of Target Vertebrae in Spine Surgery SPINE Lo, S. L., Otake, Y., Puvanesarajah, V., Wang, A. S., Uneri, A., De Silva, T., Vogt, S., Kleinszig, G., Elder, B. D., Goodwin, C., Kosztowski, T. A., Liauw, J. A., Groves, M., Bydon, A., Sciubba, D. M., Witham, T. F., Wolinsky, J., Aygun, N., Gokaslan, Z. L., Siewerdsen, J. H. 2015; 40 (8): E476–E483

  Abstract

  A 3-dimensional-2-dimensional (3D-2D) image registration algorithm, "LevelCheck," was used to automatically label vertebrae in intraoperative mobile radiographs obtained during spine surgery. Accuracy, computation time, and potential failure modes were evaluated in a retrospective study of 20 patients.To measure the performance of the LevelCheck algorithm using clinical images acquired during spine surgery.In spine surgery, the potential for wrong level surgery is significant due to the difficulty of localizing target vertebrae based solely on visual impression, palpation, and fluoroscopy. To remedy this difficulty and reduce the risk of wrong-level surgery, our team introduced a program (dubbed LevelCheck) to automatically localize target vertebrae in mobile radiographs using robust 3D-2D image registration to preoperative computed tomographic (CT) scan.Twenty consecutive patients undergoing thoracolumbar spine surgery, for whom both a preoperative CT scan and an intraoperative mobile radiograph were available, were retrospectively analyzed. A board-certified neuroradiologist determined the "true" vertebra levels in each radiograph. Registration of the preoperative CT scan to the intraoperative radiograph was calculated via LevelCheck, and projection distance errors were analyzed. Five hundred random initializations were performed for each patient, and algorithm settings (viz, the number of robust multistarts, ranging 50-200) were varied to evaluate the trade-off between registration error and computation time. Failure mode analysis was performed by individually analyzing unsuccessful registrations (>5 mm distance error) observed with 50 multistarts.At 200 robust multistarts (computation time of ∼26 s), the registration accuracy was 100% across all 10,000 trials. As the number of multistarts (and computation time) decreased, the registration remained fairly robust, down to 99.3% registration accuracy at 50 multistarts (computation time ∼7 s).The LevelCheck algorithm correctly identified target vertebrae in intraoperative mobile radiographs of the thoracolumbar spine, demonstrating acceptable computation time, compatibility with routinely obtained preoperative CT scans, and warranting investigation in prospective studies.N/A.

  View details for PubMedID 25646750

  View details for PubMedCentralID PMC4433144

• 3D-2D registration in mobile radiographs: algorithm development and preliminary clinical evaluation PHYSICS IN MEDICINE AND BIOLOGY Otake, Y., Wang, A. S., Uneri, A., Kleinszig, G., Vogt, S., Aygun, N., Lo, S. L., Wolinsky, J., Gokaslan, Z. L., Siewerdsen, J. H. 2015; 60 (5): 2075–90

  Abstract

  An image-based 3D-2D registration method is presented using radiographs acquired in the uncalibrated, unconstrained geometry of mobile radiography. The approach extends a previous method for six degree-of-freedom (DOF) registration in C-arm fluoroscopy (namely 'LevelCheck') to solve the 9-DOF estimate of geometry in which the position of the source and detector are unconstrained. The method was implemented using a gradient correlation similarity metric and stochastic derivative-free optimization on a GPU. Development and evaluation were conducted in three steps. First, simulation studies were performed that involved a CT scan of an anthropomorphic body phantom and 1000 randomly generated digitally reconstructed radiographs in posterior-anterior and lateral views. A median projection distance error (PDE) of 0.007 mm was achieved with 9-DOF registration compared to 0.767 mm for 6-DOF. Second, cadaver studies were conducted using mobile radiographs acquired in three anatomical regions (thorax, abdomen and pelvis) and three levels of source-detector distance (~800, ~1000 and ~1200 mm). The 9-DOF method achieved a median PDE of 0.49 mm (compared to 2.53 mm for the 6-DOF method) and demonstrated robustness in the unconstrained imaging geometry. Finally, a retrospective clinical study was conducted with intraoperative radiographs of the spine exhibiting real anatomical deformation and image content mismatch (e.g. interventional devices in the radiograph that were not in the CT), demonstrating a PDE = 1.1 mm for the 9-DOF approach. Average computation time was 48.5 s, involving 687 701 function evaluations on average, compared to 18.2 s for the 6-DOF method. Despite the greater computational load, the 9-DOF method may offer a valuable tool for target localization (e.g. decision support in level counting) as well as safety and quality assurance checks at the conclusion of a procedure (e.g. overlay of planning data on the radiograph for verification of the surgical product) in a manner consistent with natural surgical workflow.

  View details for PubMedID 25674851

  View details for PubMedCentralID PMC4640192

• Asymmetric Scatter Kernels for Software-Based Scatter Correction of Gridless Mammography Wang, A., Shapiro, E., Yoon, S., Ganguly, A., Proano, C., Colbeth, R., Lehto, E., Star-Lack, J., Hoeschen, C., Kontos, D. SPIE-INT SOC OPTICAL ENGINEERING. 2015

  View details for DOI 10.1117/12.2081501

  View details for Web of Science ID 000355581700050

• Known-Component 3D-2D Registration for Image Guidance and Quality Assurance in Spine Surgery Pedicle Screw Placement Uneri, A., Stayman, J. W., De Silva, T., Wang, A. S., Kleinszig, G., Vogt, S., Khanna, A. J., Wolinsky, J., Gokaslan, Z. L., Siewerdsen, J. H., Yaniv, Z. R., Webster, R. J. SPIE-INT SOC OPTICAL ENGINEERING. 2015

  View details for DOI 10.1117/12.2082210

  View details for Web of Science ID 000354365300049

• Evaluation of low-dose limits in 3D-2D rigid registration for surgical guidance Uneri, A., Wang, A. S., Otake, Y., Kleinszig, G., Vogt, S., Khanna, A. J., Gallia, G. L., Gokaslan, Z. L., Siewerdsen, J. H. IOP PUBLISHING LTD. 2014: 5329–45

  Abstract

  An algorithm for intensity-based 3D-2D registration of CT and C-arm fluoroscopy is evaluated for use in surgical guidance, specifically considering the low-dose limits of the fluoroscopic x-ray projections. The registration method is based on a framework using the covariance matrix adaptation evolution strategy (CMA-ES) to identify the 3D patient pose that maximizes the gradient information similarity metric. Registration performance was evaluated in an anthropomorphic head phantom emulating intracranial neurosurgery, using target registration error (TRE) to characterize accuracy and robustness in terms of 95% confidence upper bound in comparison to that of an infrared surgical tracking system. Three clinical scenarios were considered: (1) single-view image+guidance, wherein a single x-ray projection is used for visualization and 3D-2D guidance; (2) dual-view image+guidance, wherein one projection is acquired for visualization, combined with a second (lower-dose) projection acquired at a different C-arm angle for 3D-2D guidance; and (3) dual-view guidance, wherein both projections are acquired at low dose for the purpose of 3D-2D guidance alone (not visualization). In each case, registration accuracy was evaluated as a function of the entrance surface dose associated with the projection view(s). Results indicate that images acquired at a dose as low as 4 μGy (approximately one-tenth the dose of a typical fluoroscopic frame) were sufficient to provide TRE comparable or superior to that of conventional surgical tracking, allowing 3D-2D guidance at a level of dose that is at most 10% greater than conventional fluoroscopy (scenario #2) and potentially reducing the dose to approximately 20% of the level in a conventional fluoroscopically guided procedure (scenario #3).

  View details for DOI 10.1088/0031-9155/59/18/5329

  View details for Web of Science ID 000341381900011

  View details for PubMedID 25146673

• dPIRPLE: a joint estimation framework for deformable registration and penalized-likelihood CT image reconstruction using prior images PHYSICS IN MEDICINE AND BIOLOGY Dang, H., Wang, A. S., Sussman, M. S., Siewerdsen, J. H., Stayman, J. W. 2014; 59 (17): 4799–4826

  Abstract

  Sequential imaging studies are conducted in many clinical scenarios. Prior images from previous studies contain a great deal of patient-specific anatomical information and can be used in conjunction with subsequent imaging acquisitions to maintain image quality while enabling radiation dose reduction (e.g., through sparse angular sampling, reduction in fluence, etc). However, patient motion between images in such sequences results in misregistration between the prior image and current anatomy. Existing prior-image-based approaches often include only a simple rigid registration step that can be insufficient for capturing complex anatomical motion, introducing detrimental effects in subsequent image reconstruction. In this work, we propose a joint framework that estimates the 3D deformation between an unregistered prior image and the current anatomy (based on a subsequent data acquisition) and reconstructs the current anatomical image using a model-based reconstruction approach that includes regularization based on the deformed prior image. This framework is referred to as deformable prior image registration, penalized-likelihood estimation (dPIRPLE). Central to this framework is the inclusion of a 3D B-spline-based free-form-deformation model into the joint registration-reconstruction objective function. The proposed framework is solved using a maximization strategy whereby alternating updates to the registration parameters and image estimates are applied allowing for improvements in both the registration and reconstruction throughout the optimization process. Cadaver experiments were conducted on a cone-beam CT testbench emulating a lung nodule surveillance scenario. Superior reconstruction accuracy and image quality were demonstrated using the dPIRPLE algorithm as compared to more traditional reconstruction methods including filtered backprojection, penalized-likelihood estimation (PLE), prior image penalized-likelihood estimation (PIPLE) without registration, and prior image penalized-likelihood estimation with rigid registration of a prior image (PIRPLE) over a wide range of sampling sparsity and exposure levels.

  View details for DOI 10.1088/0031-9155/59/17/4799

  View details for Web of Science ID 000341328200004

  View details for PubMedID 25097144

  View details for PubMedCentralID PMC4142353

• Deformable image registration with local rigidity constraints for cone-beam CT-guided spine surgery PHYSICS IN MEDICINE AND BIOLOGY Reaungamornrat, S., Wang, A. S., Uneri, A., Otake, Y., Khanna, A. J., Siewerdsen, J. H. 2014; 59 (14): 3761–87

  Abstract

  Image-guided spine surgery (IGSS) is associated with reduced co-morbidity and improved surgical outcome. However, precise localization of target anatomy and adjacent nerves and vessels relative to planning information (e.g., device trajectories) can be challenged by anatomical deformation. Rigid registration alone fails to account for deformation associated with changes in spine curvature, and conventional deformable registration fails to account for rigidity of the vertebrae, causing unrealistic distortions in the registered image that can confound high-precision surgery. We developed and evaluated a deformable registration method capable of preserving rigidity of bones while resolving the deformation of surrounding soft tissue. The method aligns preoperative CT to intraoperative cone-beam CT (CBCT) using free-form deformation (FFD) with constraints on rigid body motion imposed according to a simple intensity threshold of bone intensities. The constraints enforced three properties of a rigid transformation-namely, constraints on affinity (AC), orthogonality (OC), and properness (PC). The method also incorporated an injectivity constraint (IC) to preserve topology. Physical experiments involving phantoms, an ovine spine, and a human cadaver as well as digital simulations were performed to evaluate the sensitivity to registration parameters, preservation of rigid body morphology, and overall registration accuracy of constrained FFD in comparison to conventional unconstrained FFD (uFFD) and Demons registration. FFD with orthogonality and injectivity constraints (denoted FFD+OC+IC) demonstrated improved performance compared to uFFD and Demons. Affinity and properness constraints offered little or no additional improvement. The FFD+OC+IC method preserved rigid body morphology at near-ideal values of zero dilatation (D = 0.05, compared to 0.39 and 0.56 for uFFD and Demons, respectively) and shear (S = 0.08, compared to 0.36 and 0.44 for uFFD and Demons, respectively). Target registration error (TRE) was similarly improved for FFD+OC+IC (0.7 mm), compared to 1.4 and 1.8 mm for uFFD and Demons. Results were validated in human cadaver studies using CT and CBCT images, with FFD+OC+IC providing excellent preservation of rigid morphology and equivalent or improved TRE. The approach therefore overcomes distortions intrinsic to uFFD and could better facilitate high-precision IGSS.

  View details for DOI 10.1088/0031-9155/59/14/3761

  View details for Web of Science ID 000338771300008

  View details for PubMedID 24937093

  View details for PubMedCentralID PMC4118832

• Low-dose preview for patient-specific, task-specific technique selection in cone-beam CT MEDICAL PHYSICS Wang, A. S., Stayman, J., Otake, Y., Vogt, S., Kleinszig, G., Khanna, A., Gallia, G. L., Siewerdsen, J. H. 2014; 41 (7): 071915

  Abstract

  A method is presented for generating simulated low-dose cone-beam CT (CBCT) preview images from which patient- and task-specific minimum-dose protocols can be confidently selected prospectively in clinical scenarios involving repeat scans.In clinical scenarios involving a series of CBCT images, the low-dose preview (LDP) method operates upon the first scan to create a projection dataset that accurately simulates the effects of dose reduction in subsequent scans by injecting noise of proper magnitude and correlation, including both quantum and electronic readout noise as important components of image noise in flat-panel detector CBCT. Experiments were conducted to validate the LDP method in both a head phantom and a cadaveric torso by performing CBCT acquisitions spanning a wide dose range (head: 0.8-13.2 mGy, body: 0.8-12.4 mGy) with a prototype mobile C-arm system. After injecting correlated noise to simulate dose reduction, the projections were reconstructed using both conventional filtered backprojection (FBP) and an iterative, model-based image reconstruction method (MBIR). The LDP images were then compared to real CBCT images in terms of noise magnitude, noise-power spectrum (NPS), spatial resolution, contrast, and artifacts.For both FBP and MBIR, the LDP images exhibited accurate levels of spatial resolution and contrast that were unaffected by the correlated noise injection, as expected. Furthermore, the LDP image noise magnitude and NPS were in strong agreement with real CBCT images acquired at the corresponding, reduced dose level across the entire dose range considered. The noise magnitude agreed within 7% for both the head phantom and cadaveric torso, and the NPS showed a similar level of agreement up to the Nyquist frequency. Therefore, the LDP images were highly representative of real image quality across a broad range of dose and reconstruction methods. On the other hand, naïve injection ofuncorrelated noise resulted in strong underestimation of the true noise, which would lead to overly optimistic predictions of dose reduction.Correlated noise injection is essential to accurate simulation of CBCT image quality at reduced dose. With the proposed LDP method, the user can prospectively select patient-specific, minimum-dose protocols (viz., acquisition technique and reconstruction method) suitable to a particular imaging task and to the user's own observer preferences for CBCT scans following the first acquisition. The method could provide dose reduction in common clinical scenarios involving multiple CBCT scans, such as image-guided surgery and radiotherapy.

  View details for PubMedID 24989393

• Efficacy of fixed filtration for rapid kVp-switching dual energy x-ray systems. Medical physics Yao, Y., Wang, A. S., Pelc, N. J. 2014; 41 (3): 031914-?

  Abstract

  Dose efficiency of dual kVp imaging can be improved if the two beams are filtered to remove photons in the common part of their spectra, thereby increasing spectral separation. While there are a number of advantages to rapid kVp-switching for dual energy, it may not be feasible to have two different filters for the two spectra. Therefore, the authors are interested in whether a fixed added filter can improve the dose efficiency of kVp-switching dual energy x-ray systems.The authors hypothesized that a K-edge filter would provide the energy selectivity needed to remove overlap of the spectra and hence increase the precision of material separation at constant dose. Preliminary simulations were done using calcium and water basis materials and 80 and 140 kVp x-ray spectra. Precision of the decomposition was evaluated based on the propagation of the Poisson noise through the decomposition function. Considering availability and cost, the authors chose a commercial Gd2O2S screen as the filter for their experimental validation. Experiments were conducted on a table-top system using a phantom with various thicknesses of acrylic and copper and 70 and 125 kVp x-ray spectra. The authors kept the phantom exposure roughly constant with and without filtration by adjusting the tube current. The filtered and unfiltered raw data of both low and high energy were decomposed into basis material and the variance of the decomposition for each thickness pair was calculated. To evaluate the filtration performance, the authors measured the ratio of material decomposition variance with and without filtration.Simulation results show that the ideal filter material depends on the object composition and thickness, and ranges across the lanthanide series, with higher atomic number filters being preferred for more attenuating objects. Variance reduction increases with filter thickness, and substantial reductions (40%) can be achieved with a 2× loss in intensity. The authors' experimental results validate the simulations, yet were overall slightly worse than expectation. For large objects, conventional (non-K-edge) beam hardening filters perform well.This study demonstrates the potential of fixed K-edge filtration to improve the dose efficiency and material decomposition precision for rapid kVp-switching dual energy systems.

  View details for DOI 10.1118/1.4866381

  View details for PubMedID 24593732

• Soft-tissue imaging with C-arm cone-beam CT using statistical reconstruction PHYSICS IN MEDICINE AND BIOLOGY Wang, A. S., Stayman, J., Otake, Y., Kleinszig, G., Vogt, S., Gallia, G. L., Khanna, A., Siewerdsen, J. H. 2014; 59 (4): 1005–26

  Abstract

  The potential for statistical image reconstruction methods such as penalized-likelihood (PL) to improve C-arm cone-beam CT (CBCT) soft-tissue visualization for intraoperative imaging over conventional filtered backprojection (FBP) is assessed in this work by making a fair comparison in relation to soft-tissue performance. A prototype mobile C-arm was used to scan anthropomorphic head and abdomen phantoms as well as a cadaveric torso at doses substantially lower than typical values in diagnostic CT, and the effects of dose reduction via tube current reduction and sparse sampling were also compared. Matched spatial resolution between PL and FBP was determined by the edge spread function of low-contrast (∼ 40-80 HU) spheres in the phantoms, which were representative of soft-tissue imaging tasks. PL using the non-quadratic Huber penalty was found to substantially reduce noise relative to FBP, especially at lower spatial resolution where PL provides a contrast-to-noise ratio increase up to 1.4-2.2 × over FBP at 50% dose reduction across all objects. Comparison of sampling strategies indicates that soft-tissue imaging benefits from fully sampled acquisitions at dose above ∼ 1.7 mGy and benefits from 50% sparsity at dose below ∼ 1.0 mGy. Therefore, an appropriate sampling strategy along with the improved low-contrast visualization offered by statistical reconstruction demonstrates the potential for extending intraoperative C-arm CBCT to applications in soft-tissue interventions in neurosurgery as well as thoracic and abdominal surgeries by overcoming conventional tradeoffs in noise, spatial resolution, and dose.

  View details for PubMedID 24504126

  View details for PubMedCentralID PMC4046706

• Dual-energy cone-beam CT with a flat-panel detector: Effect of reconstruction algorithm on material classification MEDICAL PHYSICS Zbijewski, W., Gang, G. J., Xu, J., Wang, A. S., Stayman, J. W., Taguchi, K., Carrino, J. A., Siewerdsen, J. H. 2014; 41 (2): 021908

  Abstract

  Cone-beam CT (CBCT) with a flat-panel detector (FPD) is finding application in areas such as breast and musculoskeletal imaging, where dual-energy (DE) capabilities offer potential benefit. The authors investigate the accuracy of material classification in DE CBCT using filtered backprojection (FBP) and penalized likelihood (PL) reconstruction and optimize contrast-enhanced DE CBCT of the joints as a function of dose, material concentration, and detail size.Phantoms consisting of a 15 cm diameter water cylinder with solid calcium inserts (50-200 mg/ml, 3-28.4 mm diameter) and solid iodine inserts (2-10 mg/ml, 3-28.4 mm diameter), as well as a cadaveric knee with intra-articular injection of iodine were imaged on a CBCT bench with a Varian 4343 FPD. The low energy (LE) beam was 70 kVp (+0.2 mm Cu), and the high energy (HE) beam was 120 kVp (+0.2 mm Cu, +0.5 mm Ag). Total dose (LE+HE) was varied from 3.1 to 15.6 mGy with equal dose allocation. Image-based DE classification involved a nearest distance classifier in the space of LE versus HE attenuation values. Recognizing the differences in noise between LE and HE beams, the LE and HE data were differentially filtered (in FBP) or regularized (in PL). Both a quadratic (PLQ) and a total-variation penalty (PLTV) were investigated for PL. The performance of DE CBCT material discrimination was quantified in terms of voxelwise specificity, sensitivity, and accuracy.Noise in the HE image was primarily responsible for classification errors within the contrast inserts, whereas noise in the LE image mainly influenced classification in the surrounding water. For inserts of diameter 28.4 mm, DE CBCT reconstructions were optimized to maximize the total combined accuracy across the range of calcium and iodine concentrations, yielding values of ∼ 88% for FBP and PLQ, and ∼ 95% for PLTV at 3.1 mGy total dose, increasing to ∼ 95% for FBP and PLQ, and ∼ 98% for PLTV at 15.6 mGy total dose. For a fixed iodine concentration of 5 mg/ml and reconstructions maximizing overall accuracy across the range of insert diameters, the minimum diameter classified with accuracy >80% was ∼ 15 mm for FBP and PLQ and ∼ 10 mm for PLTV, improving to ∼ 7 mm for FBP and PLQ and ∼ 3 mm for PLTV at 15.6 mGy. The results indicate similar performance for FBP and PLQ and showed improved classification accuracy with edge-preserving PLTV. A slight preference for increased smoothing of the HE data was found. DE CBCT discrimination of iodine and bone in the knee was demonstrated with FBP and PLTV at 6.2 mGy total dose.For iodine concentrations >5 mg/ml and detail size ∼ 20 mm, material classification accuracy of >90% was achieved in DE CBCT with both FBP and PL at total doses <10 mGy. Optimal performance was attained by selection of reconstruction parameters based on the differences in noise between HE and LE data, typically favoring stronger smoothing of the HE data, and by using penalties matched to the imaging task (e.g., edge-preserving PLTV in areas of uniform enhancement).

  View details for DOI 10.1118/1.4863598

  View details for Web of Science ID 000331213300041

  View details for PubMedID 24506629

  View details for PubMedCentralID PMC3977791

• 3D-2D registration for surgical guidance: effect of projection view angles on registration accuracy PHYSICS IN MEDICINE AND BIOLOGY Uneri, A., Otake, Y., Wang, A. S., Kleinszig, G., Vogt, S., Khanna, A. J., Siewerdsen, J. H. 2014; 59 (2): 271–87

  Abstract

  An algorithm for intensity-based 3D-2D registration of CT and x-ray projections is evaluated, specifically using single- or dual-projection views to provide 3D localization. The registration framework employs the gradient information similarity metric and covariance matrix adaptation evolution strategy to solve for the patient pose in six degrees of freedom. Registration performance was evaluated in an anthropomorphic phantom and cadaver, using C-arm projection views acquired at angular separation, Δθ, ranging from ∼0°-180° at variable C-arm magnification. Registration accuracy was assessed in terms of 2D projection distance error and 3D target registration error (TRE) and compared to that of an electromagnetic (EM) tracker. The results indicate that angular separation as small as Δθ ∼10°-20° achieved TRE <2 mm with 95% confidence, comparable or superior to that of the EM tracker. The method allows direct registration of preoperative CT and planning data to intraoperative fluoroscopy, providing 3D localization free from conventional limitations associated with external fiducial markers, stereotactic frames, trackers and manual registration.

  View details for DOI 10.1088/0031-9155/59/2/271

  View details for Web of Science ID 000332842000003

  View details for PubMedID 24351769

  View details for PubMedCentralID PMC4927006

• Dual-Projection 3D-2D Registration for Surgical Guidance: Preclinical Evaluation of Performance and Minimum Angular Separation Uneri, A., Otake, Y., Wang, A. S., Kleinszig, G., Vogt, S., Gallia, G. L., Rigamonti, D., Wolinsky, J., Gokaslan, Z. L., Khanna, A. J., Siewerdsen, J. H., Yaniv, Z. R., Holmes, D. R. SPIE-INT SOC OPTICAL ENGINEERING. 2014

  View details for DOI 10.1117/12.2043561

  View details for Web of Science ID 000348029400084

• Patient-Specific Minimum-Dose Imaging Protocols for Statistical Image Reconstruction in C-arm Cone-Beam CT Using Correlated Noise Injection Wang, A. S., Stayman, J. W., Otake, Y., Khanna, A. J., Gallia, G. L., Siewerdsen, J. H., Whiting, B. R., Hoeschen, C., Kontos, D. SPIE-INT SOC OPTICAL ENGINEERING. 2014

  View details for DOI 10.1117/12.2043083

  View details for Web of Science ID 000338775800058

• Deformable Registration for Image-Guided Spine Surgery: Preserving Rigid Body Vertebral Morphology in Free-Form Transformations Reaungamornrat, S., Wang, A. S., Uneri, A., Otake, Y., Zhao, Z., Khanna, A. J., Siewerdsen, J. H., Yaniv, Z. R., Holmes, D. R. SPIE-INT SOC OPTICAL ENGINEERING. 2014

  View details for DOI 10.1117/12.2043474

  View details for Web of Science ID 000348029400027

• Robust 3D-2D image registration: application to spine interventions and vertebral labeling in the presence of anatomical deformation PHYSICS IN MEDICINE AND BIOLOGY Otake, Y., Wang, A. S., Stayman, J., Uneri, A., Kleinszig, G., Vogt, S., Khanna, A., Gokaslan, Z. L., Siewerdsen, J. H. 2013; 58 (23): 8535–53

  Abstract

  We present a framework for robustly estimating registration between a 3D volume image and a 2D projection image and evaluate its precision and robustness in spine interventions for vertebral localization in the presence of anatomical deformation. The framework employs a normalized gradient information similarity metric and multi-start covariance matrix adaptation evolution strategy optimization with local-restarts, which provided improved robustness against deformation and content mismatch. The parallelized implementation allowed orders-of-magnitude acceleration in computation time and improved the robustness of registration via multi-start global optimization. Experiments involved a cadaver specimen and two CT datasets (supine and prone) and 36 C-arm fluoroscopy images acquired with the specimen in four positions (supine, prone, supine with lordosis, prone with kyphosis), three regions (thoracic, abdominal, and lumbar), and three levels of geometric magnification (1.7, 2.0, 2.4). Registration accuracy was evaluated in terms of projection distance error (PDE) between the estimated and true target points in the projection image, including 14 400 random trials (200 trials on the 72 registration scenarios) with initialization error up to ±200 mm and ±10°. The resulting median PDE was better than 0.1 mm in all cases, depending somewhat on the resolution of input CT and fluoroscopy images. The cadaver experiments illustrated the tradeoff between robustness and computation time, yielding a success rate of 99.993% in vertebral labeling (with 'success' defined as PDE <5 mm) using 1,718 664 ± 96 582 function evaluations computed in 54.0 ± 3.5 s on a mid-range GPU (nVidia, GeForce GTX690). Parameters yielding a faster search (e.g., fewer multi-starts) reduced robustness under conditions of large deformation and poor initialization (99.535% success for the same data registered in 13.1 s), but given good initialization (e.g., ±5 mm, assuming a robust initial run) the same registration could be solved with 99.993% success in 6.3 s. The ability to register CT to fluoroscopy in a manner robust to patient deformation could be valuable in applications such as radiation therapy, interventional radiology, and an assistant to target localization (e.g., vertebral labeling) in image-guided spine surgery.

  View details for PubMedID 24246386

  View details for PubMedCentralID PMC4915373

• Noise Reduction in Material Decomposition for Low-Dose Dual-Energy Cone-Beam CT Zbijewski, W., Gang, G., Wang, A. S., Stayman, J. W., Taguchi, K., Carrino, J. A., Siewerdsen, J. H., Nishikawa, R. M., Whiting, B. R., Hoeschen, C. SPIE-INT SOC OPTICAL ENGINEERING. 2013

  View details for DOI 10.1117/12.2008431

  View details for Web of Science ID 000322002700041

• Intraoperative Imaging for Patient Safety and QA: Detection of Intracranial Hemorrhage Using C-Arm Cone-Beam CT Schafer, S., Wang, A., Otake, Y., Stayman, J., Zbijewski, W., Kleinszig, G., Xia, X., Gallia, G. L., Siewerdsen, J. H., Holmes, D. R., Yaniv, Z. R. SPIE-INT SOC OPTICAL ENGINEERING. 2013

  View details for DOI 10.1117/12.2008043

  View details for Web of Science ID 000321905800068

• Soft-Tissue Imaging in Low-Dose, C-Arm Conbe-Beam CT Using Statistical Image Reconstruction Wang, A. S., Schafer, S., Stadyman, J., Otake, Y., Sussman, M. S., Khanna, A., Gallia, G. L., Siewerdsen, J. H., Nishikawa, R. M., Whiting, B. R., Hoeschen, C. SPIE-INT SOC OPTICAL ENGINEERING. 2013

  View details for DOI 10.1117/12.2008421

  View details for Web of Science ID 000322002700047

• Efficacy of Fixed Filtration for Rapid kVp-Switching Dual Energy X-ray Systems: Experimental Verification Conference on Medical Imaging - Physics of Medical Imaging Yao, Y., Wang, A. S., Pelc, N. J. SPIE-INT SOC OPTICAL ENGINEERING. 2012

  View details for DOI 10.1117/12.913222

  View details for Web of Science ID 000304768000047

• A comparison of dual kV energy integrating and energy discriminating photon counting detectors for dual energy x-ray imaging Conference on Medical Imaging - Physics of Medical Imaging Wang, A. S., Pelc, N. J. SPIE-INT SOC OPTICAL ENGINEERING. 2012

  View details for DOI 10.1117/12.912030

  View details for Web of Science ID 000304768000029

• Image-based Synthetic CT: simulating arbitrary low dose single and dual energy protocols from dual energy images Conference on Medical Imaging - Physics of Medical Imaging Wang, A. S., Feng, C., Pelc, N. J. SPIE-INT SOC OPTICAL ENGINEERING. 2012

  View details for DOI 10.1117/12.912163

  View details for Web of Science ID 000304768000048

• Synthetic CT: Simulating low dose single and dual energy protocols from a dual energy scan MEDICAL PHYSICS Wang, A. S., Pelc, N. J. 2011; 38 (10): 5551-5562

  Abstract

  The choice of CT protocol can greatly impact patient dose and image quality. Since acquiring multiple scans at different techniques on a given patient is undesirable, the ability to predict image quality changes starting from a high quality exam can be quite useful. While existing methods allow one to generate simulated images of lower exposure (mAs) from an acquired CT exam, the authors present and validate a new method called synthetic CT that can generate realistic images of a patient at arbitrary low dose protocols (kVp, mAs, and filtration) for both single and dual energy scans.The synthetic CT algorithm is derived by carefully ensuring that the expected signal and noise are accurate for the simulated protocol. The method relies on the observation that the material decomposition from a dual energy CT scan allows the transmission of an arbitrary spectrum to be predicted. It requires an initial dual energy scan of the patient to either synthesize raw projections of a single energy scan or synthesize the material decompositions of a dual energy scan. The initial dual energy scan contributes inherent noise to the synthesized projections that must be accounted for before adding more noise to simulate low dose protocols. Therefore, synthetic CT is subject to the constraint that the synthesized data have noise greater than the inherent noise. The authors experimentally validated the synthetic CT algorithm across a range of protocols using a dual energy scan of an acrylic phantom with solutions of different iodine concentrations. An initial 80/140 kVp dual energy scan of the phantom provided the material decomposition necessary to synthesize images at 100 kVp and at 120 kVp, across a range of mAs values. They compared these synthesized single energy scans of the phantom to actual scans at the same protocols. Furthermore, material decompositions of a 100/120 kVp dual energy scan are synthesized by adding correlated noise to the initial material decompositions. The aforementioned noise constraint also allows us to compute feasible mAs values that can be synthesized for each kVp.The single energy synthesized and actual reconstructed images exhibit identical signal and noise properties at 100 kVp and at 120 kVp, and across a range of mAs values. For example, the noise in both the synthesized and actual images at 100 kVp increases by 2 when the mAs is halved. The synthesized and actual material decompositions of a dual energy protocol show excellent agreement when the decomposition images are linearly weighted to form monoenergetic images at energies from 40 to 100 keV. For simulated single energy protocols with kVp between 80 and 140, the highest feasible mAs exceeds that of either initial scan.This work describes and validates the synthetic CT theory and algorithm by comparing its results to actual scans. Synthetic CT is a powerful new tool that allows users to realistically see how protocol selection affects CT images and enables radiologists to retrospectively identify the lowest dose protocol achievable that provides diagnostic quality images on real patients.

  View details for DOI 10.1118/1.3633895

  View details for Web of Science ID 000295617400030

  View details for PubMedID 21992373

• Pulse pileup statistics for energy discriminating photon counting x-ray detectors MEDICAL PHYSICS Wang, A. S., Harrison, D., Lobastov, V., Tkaczyk, J. E. 2011; 38 (7): 4265-4275

  Abstract

  Purpose: Energy discriminating photon counting x-ray detectors can be subject to a wide range of flux rates if applied in clinical settings. Even when the incident rate is a small fraction of the detector's maximum periodic rate No, pulse pileup leads to count rate losses and spectral distortion. Although the deterministic effects can be corrected, the detrimental effect of pileup on image noise is not well understood and may limit the performance of photon counting systems. Therefore, the authors devise a method to determine the detector count statistics and imaging performance.The detector count statistics are derived analytically for an idealized pileup model with delta pulses of a nonparalyzable detector. These statistics are then used to compute the performance (e.g., contrast-to-noise ratio) for both single material and material decomposition contrast detection tasks via the Cramdr-Rao lower bound (CRLB) as a function of the detector input count rate. With more realistic unipolar and bipolar pulse pileup models of a nonparalyzable detector, the imaging task performance is determined by Monte Carlo simulations and also approximated by a multinomial method based solely on the mean detected output spectrum. Photon counting performance at different count rates is compared with ideal energy integration, which is unaffected by count rate.The authors found that an ideal photon counting detector with perfect energy resolution outperforms energy integration for our contrast detection tasks, but when the input count rate exceeds 20% N0, many of these benefits disappear. The benefit with iodine contrast falls rapidly with increased count rate while water contrast is not as sensitive to count rates. The performance with a delta pulse model is overoptimistic when compared to the more realistic bipolar pulse model. The multinomial approximation predicts imaging performance very close to the prediction from Monte Carlo simulations. The monoenergetic image with maximum contrast-to-noise ratio from dual energy imaging with ideal photon counting is only slightly better than with dual kVp energy integration, and with a bipolar pulse model, energy integration outperforms photon counting for this particular metric because of the count rate losses. However, the material resolving capability of photon counting can be superior to energy integration with dual kVp even in the presence of pileup because of the energy information available to photon counting.A computationally efficient multinomial approximation of the count statistics that is based on the mean output spectrum can accurately predict imaging performance. This enables photon counting system designers to directly relate the effect of pileup to its impact on imaging statistics and how to best take advantage of the benefits of energy discriminating photon counting detectors, such as material separation with spectral imaging.

  View details for DOI 10.1118/1.3592932

  View details for Web of Science ID 000292521100043

  View details for PubMedID 21859028

• Sufficient Statistics as a Generalization of Binning in Spectral X-ray Imaging IEEE TRANSACTIONS ON MEDICAL IMAGING Wang, A. S., Pelc, N. J. 2011; 30 (1): 84-93

  Abstract

  It is well known that the energy dependence of X-ray attenuation can be used to characterize materials. Yet, even with energy discriminating photon counting X-ray detectors, it is still unclear how to best form energy dependent measurements for spectral imaging. Common ideas include binning photon counts based on their energies and detectors with both photon counting and energy integrating electronics. These approaches can be generalized to energy weighted measurements, which we prove can form a sufficient statistic for spectral X-ray imaging if the weights used, which we term μ-weights, are basis attenuation functions that can also be used for material decomposition. To study the performance of these different methods, we evaluate the Cramér-Rao lower bound (CRLB) of material estimates in the presence of quantum noise. We found that the choice of binning and weighting schemes can greatly affect the performance of material decomposition. Even with optimized thresholds, binning condenses information but incurs penalties to decomposition precision and is not robust to changes in the source spectrum or object size, although this can be mitigated by adding more bins or removing photons of certain energies from the spectrum. On the other hand, because μ-weighted measurements form a sufficient statistic for spectral imaging, the CRLB of the material decomposition estimates is identical to the quantum noise limited performance of a system with complete energy information of all photons. Finally, we show that μ-weights lead to increased conspicuity over other methods in a simulated calcium contrast experiment.

  View details for DOI 10.1109/TMI.2010.2061862

  View details for Web of Science ID 000285844900008

  View details for PubMedID 20682470

• Contrast-to-Noise of a Non-Ideal, Multi-bin, Photon Counting X-ray Detector Tkaczyk, J., Lobastov, V., Harrison, D. D., Wang, A. S., Pelc, N. J., Samei, E., Nishikawa, R. M. SPIE-INT SOC OPTICAL ENGINEERING. 2011

  View details for DOI 10.1117/12.878290

  View details for Web of Science ID 000294178500122

• Synthetic CT: simulating arbitrary low dose single and dual energy protocols Conference on Medical Imaging 2011 - Physics of Medical Imaging Wang, A. S., Pelc, N. J. SPIE-INT SOC OPTICAL ENGINEERING. 2011

  View details for DOI 10.1117/12.878771

  View details for Web of Science ID 000294178500058

• Impact of Photon Counting Detector Spectral Response on Dual Energy Techniques Conference on Medical Imaging - Physics of Medical Imaging Wang, A. S., Pelc, N. J. SPIE-INT SOC OPTICAL ENGINEERING. 2010

  View details for DOI 10.1117/12.840052

  View details for Web of Science ID 000285047200122

• Understanding and controlling the effect of lossy raw data compression on CT images MEDICAL PHYSICS Wang, A. S., Pelc, N. J. 2009; 36 (8): 3643-3653

  Abstract

  The requirements for raw data transmission through a CT scanner slip ring, through the computation system, and for storage of raw CT data can be quite challenging as scanners continue to increase in speed and to collect more data per rotation. Although lossy compression greatly mitigates this problem, users must be cautious about how errors introduced manifest themselves in the reconstructed images. This paper describes two simple yet effective methods for controlling the effect of errors in raw data compression and describe the impact of each stage on the image errors. A CT system simulator (CATSIM, GE Global Research Center, Niskayuna, NY) was used to generate raw CT datasets that simulate different regions of human anatomy. The raw data are digitized by a 20-bit ADC and companded by a log compander. Lossy compression is performed by quantization and is followed by JPEG-LS (lossless), which takes advantage of the correlations between neighboring measurements in the sinogram. Error feedback, a previously proposed method that controls the spatial distribution of reconstructed image errors, and projection filtering, a newly proposed method that takes advantage of the filtered backprojection reconstruction process, are applied independently (and combined) to study their intended impact on the control and behavior of the additional noise due to the compression methods used. The log compander and the projection filtering method considerably reduce image error levels, while error feedback pushes image errors toward the periphery of the field of view. The results for the images are a compression ratio (CR) of 3 that keeps peak compression errors under 1 HU and a CR of 9 that increases image noise by only 1 HU in common CT applications. Lossy compression can substantially reduce raw CT data size at low computational cost. The proposed methods have the flexibility to operate at a wide range of compression ratios and produce predictable, object-independent, and often imperceptible image artifacts.

  View details for DOI 10.1118/1.3158738

  View details for Web of Science ID 000268440600029

  View details for PubMedID 19746798

• Lossy raw data compression in computed tomography with noise shaping to control image effects Medical Imaging 2008 Conference Xie, Y., Wang, A. S., Pelc, N. J. SPIE-INT SOC OPTICAL ENGINEERING. 2008

  View details for DOI 10.1117/12.769954

  View details for Web of Science ID 000256660300105

• Effect of the frequency content and spatial location of raw data errors on CT images MEDICAL IMAGING 2008: PHYSICS OF MEDICAL IMAGING, PTS 1-3 Wang, A. S., Xie, Y., Pelc, N. J. 2008; 6913

  View details for DOI 10.1117/12.770552

  View details for Web of Science ID 000256660300107

• Detection of flicker caused by high-frequency interharmonics Kim, T., Wang, A., Powers, E. J., Grady, W., Arapostathis, A., IEEE IEEE. 2007: 336-+

  View details for Web of Science ID 000251296800066