Search NASA⌕ Search

SEARCH · Search NASA

Results for “CONES”

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

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

At least 775 records · Page 43

HELIOSPHERIC PROPAGATION OF CORONAL MASS EJECTIONS: COMPARISON OF NUMERICAL WSA-ENLIL+CONE MODEL AND ANALYTICAL DRAG-BASED MODEL

Real-time forecasting of the arrival of coronal mass ejections (CMEs) at Earth, based on remote solar observations, is one of the central issues of space-weather research. In this paper, we compare arrival-time predictions calculated applying the numerical ''WSA-ENLIL+Cone model'' and the analytical ''drag-based model'' (DBM). Both models use coronagraphic observations of CMEs as input data, thus providing an early space-weather forecast two to four days before the arrival of the disturbance at the Earth, depending on the CME speed. It is shown that both methods give very similar results if the drag parameter Γ = 0.1 is used in DBM in combination with a background solar-wind speed of w = 400 km s{sup –1}. For this combination, the mean value of the difference between arrival times calculated by ENLIL and DBM is Δ-bar =0.09±9.0 hr with an average of the absolute-value differences of |Δ|-bar =7.1 hr. Comparing the observed arrivals (O) with the calculated ones (C) for ENLIL gives O – C = –0.3 ± 16.9 hr and, analogously, O – C = +1.1 ± 19.1 hr for DBM. Applying Γ = 0.2 with w = 450 km s{sup –1} in DBM, one finds O – C = –1.7 ± 18.3 hr, with an average of the absolute-value differences of 14.8 hr, which is similar to that for ENLIL, 14.1 hr. Finally, we demonstrate that the prediction accuracy significantly degrades with increasing solar activity.

79 ASTRONOMY AND ASTROPHYSICS↗

What's inside the cone? Numerically reconstructing the metric from observations

We investigate the possibility of using Gaussian process regression to smooth data on the current past null-cone for use as the input to a relativistic integration scheme. The algorithm we present is designed to reconstruct the metric of spacetime within the class of spherically symmetric dust universes, with or without a cosmological constant. Assuming that gravity is well described by General Relativity, we demonstrate how the algorithm can be employed to test the Copernican principle based on currently available observations. It is shown that currently available data is not sufficient for a conclusive result. The intrinsic noise present in realistic data presents a challenge for our smoothing algorithm and we discuss some of its limitations as well as possible extensions to it. We conclude by demonstrating how a direct determination of the cosmological constant is possible using redshift drift data.

79 ASTRONOMY AND ASTROPHYSICS↗

Synthetic light-cone catalogues of modern redshift and weak lensing surveys with abacussummit

ABSTRACT The joint analysis of different cosmological probes, such as galaxy clustering and weak lensing, can potentially yield invaluable insights into the nature of the primordial Universe, dark energy, and dark matter. However, the development of high-fidelity theoretical models is a necessary stepping stone. Here, we present public high-resolution weak lensing maps on the light-cone, generated using the N-body simulation suite abacussummit, and accompanying weak lensing mock catalogues, tuned to the Early Data Release small-scale clustering measurements of the Dark Energy Spectroscopic Instrument. Available in this release are maps of the cosmic shear, deflection angle, and convergence fields at source redshifts ranging from z = 0.15 to 2.45 as well as cosmic microwave background convergence maps for each of the 25 base-resolution simulations ($L_{\rm box} = 2000\, h^{-1}\, {\rm Mpc}$ and Npart = 69123) as well as for the two huge simulations ($L_{\rm box} = 7500\, h^{-1}\, {\rm Mpc}$ and Npart = 86403) at the fiducial abacussummit cosmology. The pixel resolution of each map is 0.21 arcmin, corresponding to a healpix Nside of 16 384. The sky coverage of the base simulations is an octant until z ≈ 0.8 (decreasing to about 1800 deg2 at z ≈ 2.4), whereas the huge simulations offer full-sky coverage until z ≈ 2.2. Mock lensing source catalogues are sampled matching the ensemble properties of the Kilo-Degree Survey, Dark Energy Survey, and Hyper Suprime-Cam data sets. The mock catalogues are validated against theoretical predictions for various clustering and lensing statistics, such as correlation multipoles, galaxy–shear, and shear–shear, showing excellent agreement. All products can be downloaded via a Globus endpoint (see Data Availability section).

79 ASTRONOMY AND ASTROPHYSICS↗

Full forward model of galaxy clustering statistics with AbacusSummit light cones

ABSTRACT Novel summary statistics beyond the standard 2-point correlation function (2PCF) are necessary to capture the full astrophysical and cosmological information from the small-scale (r < 30h−1Mpc) galaxy clustering. However, the analysis of beyond-2PCF statistics on small scales is challenging because we lack the appropriate treatment of observational systematics for arbitrary summary statistics of the galaxy field. In this paper, we develop a full forward modelling pipeline for a wide range of summary statistics using the large high-fidelity AbacusSummit light cones that account for many systematic effects as well as remain flexible and computationally efficient to enable posterior sampling. We apply our forward model approach to a fully realistic mock galaxy catalog and demonstrate that we can recover unbiased constraints on the underlying galaxy–halo connection model using two separate summary statistics: the standard 2PCF and the novel k-th nearest neighbour (kNN) statistics, which are sensitive to correlation functions of all orders. We will demonstrate its strong constraining power on extended galaxy–halo connection models and cosmology in follow up papers. We expect this to become a powerful approach when applying to upcoming surveys such as DESI where we can leverage a multitude of summary statistics across a wide redshift range to maximally extract information from the non-linear scales.

79 ASTRONOMY AND ASTROPHYSICS↗

Taking apart the enhanced backscattering cone: Interference fringes from reciprocal paths in multiple light scattering

We report the decomposition of the enhanced backscattering cone into its constitutive interference fringes. These fringes are due to the constructive interference between reciprocal paths of any multiply scattered wave after ensemble averaging. An optical setup combining a two-point continuous-wave illumination and matching detection allows the observation of the fringes and, therefore, the quantitative characterization of the Green's function for light propagation between the two points in a multiple-scattering media.

74 ATOMIC AND MOLECULAR PHYSICS↗

Far-infrared absorption of undoped and Br-doped carbon nanofiber powder in stacked-cup cone configuration

In this study, we performed room-temperature far-infrared (40–650cm –1 ) transmission measurements on undoped and bromine-doped powder samples of carbon nanofibers in stacked-cup cone geometry. The transmission spectra show enhanced transmittance after bromine doping and all spectra were fit to a Drude-Lorentz (DL) model. A decreased metallic conductivity along with a redshift of the lowest semiconducting gap was found with doped samples. A significant decrease in the scattering rate upon heavy doping has been qualitatively explained as partial ordering of intercalated dopant ions. Absorption spectra were derived from the transmission spectra under the assumption of nondispersive reflectance and subsequently compared with DL-model-derived spectra. The free-carrier density of the n -type powder and the electronic mean free path were estimated and compared with reported values for single-walled nanotubes and pyrolytic graphite.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Type-II Dirac cones and electron-phonon interaction in monolayer biphenylene from first-principles calculations

In this work, we report a first-principles investigation of electronic structure, topological bands, and electron-phonon interactions in metallic biphenylene sheets. Biphenylene is a recently synthesized sp 2 -bonded carbon allotrope. We find coupling of electrons at the Fermi surface to very high frequency carbon-derived phonons, analogous to superconducting MgB 2 . This leads to low-temperature weak coupling superconductivity due to an unusual combination of exceptionally large logarithmically averaged phonon frequency ω log =1369 K and moderate electron-phonon coupling. The electronic structure shows a two-band Fermi surface dominated by C p z orbitals and a pair of type-II tilted Dirac cones along the Γ-Y line at the Brillouin zone boundary. Berry curvature and edge-state calculations show that monolayer biphenylene is a two-dimensional $\mathbb{Z}$ 2 topological material. Thus, monolayer biphenylene is predicted to be a topological superconductor based on C p orbitals and high-frequency phonons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Velocity reconstruction in the era of DESI and Rubin/LSST. II. Realistic samples on the light cone

Reconstructing the galaxy peculiar velocity field from the distribution of large-scale structure plays an important role in cosmology. On one hand, it gives us an insight into structure formation and gravity; on the other, it allows us to selectively extract the kinetic Sunyaev-Zel’dovich (kSZ) effect from cosmic microwave background maps. In this work, we employ high-accuracy synthetic galaxy catalogs on the light cone to investigate how well we can recover the velocity field when utilizing the three-dimensional spatial distribution of the galaxies in a modern large-scale structure experiment such as the Dark Energy Spectroscopic Instrument (DESI) and the Rubin Observatory Legacy Survey of Space and Time. In particular, we adopt the standard technique used in baryon acoustic oscillation analysis for reconstructing the Zel’dovich displacements of galaxies through the continuity equation, which yields a first-order approximation to their large-scale velocities. We investigate variations in the number density, bias, mask, area, redshift noise, and survey depth, as well as modifications to the settings of the standard reconstruction algorithm. Since our main goal is to provide guidance for planned kSZ analysis between DESI and the Atacama Cosmology Telescope, we apply velocity reconstruction to a faithful representation of DESI spectroscopic and photometric targets. We report the cross-correlation coefficient between the reconstructed and the true velocities along the line of sight. For the DESI Y1 spectroscopic survey, we expect the correlation coefficient to be r ≈ 0.64, while for a photometric survey with δ z /(1+z) = 0.02, as is approximately the case for the Legacy Survey used in the target selection of DESI galaxies, r shrinks by half to r ≈ 0.31. Here, we hope the results in this paper can be used to inform future kSZ stacking studies and other velocity reconstruction analyses planned with the next generation of cosmology experiments.

79 ASTRONOMY AND ASTROPHYSICS↗

Transport of an intense proton beam from a cone-structured target through plastic foam with unique proton source modeling

Laser-accelerated proton beams are applicable to several research areas within high-energy density science, including warm dense matter generation, proton radiography, and inertial confinement fusion, which all involve transport of the beam through matter. Here, we report on experimental measurements of intense proton beam transport through plastic foam blocks. The intense proton beam was accelerated by the 10 ps, 700 $\textit{J}$ OMEGA EP laser irradiating a curved foil target, and focused by an attached hollow cone. The protons then entered the foam block of density 0.38 g/cm 3 and thickness 0.55 or 1.00 mm. At the rear of the foam block, a Cu layer revealed the cross section of the intense beam via proton- and hot electron-induced Cu-K α emission. Images of x-ray emission show a bright spot on the rear Cu film indicative of a forward-directed beam without major breakup. 2D fluid-PIC simulations of the transport were conducted using a unique multi-injection source model incorporating energy-dependent beam divergence. Along with postprocessed calculations of the Cu – $K_α$ emission profile, simulations showed that protons retain their ballistic transport through the foam and are able to heat the foam up to several keV in temperature. The total experimental emission profile for the 1.0 mm foam agrees qualitatively with the simulated profile, suggesting that the protons indeed retain their beamlike qualities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Cone beam neutron interferometry: From modeling to applications

Phase-grating moiré interferometers (PGMIs) have emerged as promising candidates for the next generation of neutron interferometry, enabling the use of a polychromatic beam and manifesting interference patterns that can be directly imaged by existing neutron cameras. However, the modeling of the various PGMI configurations is limited to cumbersome numerical calculations and backward propagation models which often do not enable one to explore the setup parameters. Here we generalize the Fresnel scaling theorem to introduce a k -space model for PGMI setups illuminated by a cone beam, thus enabling an intuitive forward propagation model for a wide range of parameters and experimental setups. The interference manifested by a PGMI is shown to be a special case of the Talbot effect, and the optimal fringe visibility is shown to occur at the moiré location of the Talbot distances. We derive analytical expressions for the contrast and the propagating intensity profiles in various conditions and provide the first analysis of the PGMI dark-field imaging signal when considering sample characterization. The model's predictions are compared to experimental measurements and good agreement is found between them. Last, we propose and experimentally verify a method to recover contrast at typically inaccessible PGMI autocorrelation lengths. The presented work provides a toolbox for analyzing and understanding existing PGMI setups and their future applications, for example extensions to two-dimensional PGMIs and characterization of samples with nontrivial structures. Published by the American Physical Society 2024

Sarenac, D. (ORCID:0000000185753367)↗

Anatomical background and generalized detectability in tomosynthesis and cone-beam CT

Purpose: Anatomical background presents a major impediment to detectability in 2D radiography as well as 3D tomosynthesis and cone-beam CT (CBCT). This article incorporates theoretical and experimental analysis of anatomical background ''noise'' in cascaded systems analysis of 2D and 3D imaging performance to yield ''generalized'' metrics of noise-equivalent quanta (NEQ) and detectability index as a function of the orbital extent of the (circular arc) source-detector orbit. Methods: A physical phantom was designed based on principles of fractal self-similarity to exhibit power-law spectral density ({kappa}/f{sup {beta}}) comparable to various anatomical sites (e.g., breast and lung). Background power spectra [S{sub B}(f)] were computed as a function of source-detector orbital extent, including tomosynthesis ({approx}10 deg. - 180 deg.) and CBCT (180 deg. +fan to 360 deg.) under two acquisition schemes: (1) Constant angular separation between projections (variable dose) and (2) constant total number of projections (constant dose). The resulting S{sub B} was incorporated in the generalized NEQ, and detectability index was computed from 3D cascaded systems analysis for a variety of imaging tasks. Results: The phantom yielded power-law spectra within the expected spatial frequency range, quantifying the dependence of clutter magnitude ({kappa}) and correlation ({beta}) with increasing tomosynthesis angle. Incorporation of S{sub B} in the 3D NEQ provided a useful framework for analyzing the tradeoffs among anatomical, quantum, and electronic noise with dose and orbital extent. Distinct implications are posed for breast and chest tomosynthesis imaging system design--applications varying significantly in {kappa} and {beta}, and imaging task and, therefore, in optimal selection of orbital extent, number of projections, and dose. For example, low-frequency tasks (e.g., soft-tissue masses or nodules) tend to benefit from larger orbital extent and more fully 3D tomographic imaging, whereas high-frequency tasks (e.g., microcalcifications) require careful, application-specific selection of orbital extent and number of projections to minimize negative effects of quantum and electronic noise. Conclusions: The complex tradeoffs among anatomical background, quantum noise, and electronic noise in projection imaging, tomosynthesis, and CBCT can be described by generalized cascaded systems analysis, providing a useful framework for system design and optimization.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

GPU-based fast cone beam CT reconstruction from undersampled and noisy projection data via total variation

Purpose: Cone-beam CT (CBCT) plays an important role in image guided radiation therapy (IGRT). However, the large radiation dose from serial CBCT scans in most IGRT procedures raises a clinical concern, especially for pediatric patients who are essentially excluded from receiving IGRT for this reason. The goal of this work is to develop a fast GPU-based algorithm to reconstruct CBCT from undersampled and noisy projection data so as to lower the imaging dose. Methods: The CBCT is reconstructed by minimizing an energy functional consisting of a data fidelity term and a total variation regularization term. The authors developed a GPU-friendly version of the forward-backward splitting algorithm to solve this model. A multigrid technique is also employed. Results: It is found that 20-40 x-ray projections are sufficient to reconstruct images with satisfactory quality for IGRT. The reconstruction time ranges from 77 to 130 s on an NVIDIA Tesla C1060 (NVIDIA, Santa Clara, CA) GPU card, depending on the number of projections used, which is estimated about 100 times faster than similar iterative reconstruction approaches. Moreover, phantom studies indicate that the algorithm enables the CBCT to be reconstructed under a scanning protocol with as low as 0.1 mA s/projection. Comparing with currently widely used full-fan head and neck scanning protocol of {approx}360 projections with 0.4 mA s/projection, it is estimated that an overall 36-72 times dose reduction has been achieved in our fast CBCT reconstruction algorithm. Conclusions: This work indicates that the developed GPU-based CBCT reconstruction algorithm is capable of lowering imaging dose considerably. The high computation efficiency in this algorithm makes the iterative CBCT reconstruction approach applicable in real clinical environments.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Balancing dose and image registration accuracy for cone beam tomosynthesis (CBTS) for breast patient setup

Purpose: To balance dose reduction and image registration accuracy in breast setup imaging. In particular, the authors demonstrate the relationship between scan angle and dose delivery for cone beam tomosynthesis (CBTS) when employed for setup verification of breast cancer patients with surgical clips. Methods: The dose measurements were performed in a female torso phantom for varying scan angles of CBTS. Setup accuracy was measured using three registration methods: Clip centroid localization accuracy and the accuracy of two semiautomatic registration algorithms. The dose to the organs outside of the ipsilateral breast and registration accuracy information were compared to determine the optimal scan angle for CBTS for breast patient setup verification. Isocenter positions at the center of the patient and at the breast-chest wall interface were considered. Results: Image registration accuracy was within 1 mm for the CBTS scan angles {theta} above 20 deg. for some scenarios and as large as 80 deg. for the worst case, depending on the imaged breast and registration algorithm. Registration accuracy was highest based on clip centroid localization. For left and right breast imaging with the isocenter at the chest wall, the dose to the contralateral side of the patient was very low (<0.5 cGy) for all scan angles considered. For central isocenter location, the optimal scan angles were 30 deg. - 50 deg. for the left breast imaging and 40 deg. - 50 deg. for the right breast imaging, with the difference due to the geometric asymmetry of the current clinical imaging system. Conclusions: The optimal scan angles for CBTS imaging were found to be between 10 deg. and 50 deg., depending on the isocenter location and ipsilateral breast. Use of the isocenter at the breast-chest wall locations always resulted in greater accuracy of image registration (<1 mm) at smaller angles (10 deg. - 20 deg.) and at lower doses (<0.1 cGy) to the contralateral organs. For chest wall isocenters, doses delivered to organs outside of the target breast were much smaller than the scattered and leakage doses of the treatment beams. The complete volumetric information of all clips in the region of interest, combined with the small dose to the contralateral organs and the small scan angle, could result in an advantage for small angle CBTS with off center isocenters over simple orthogonal pairs.

60 APPLIED LIFE SCIENCES↗

Shading correction for on-board cone-beam CT in radiation therapy using planning MDCT images

Purpose: Applications of cone-beam CT (CBCT) to image-guided radiation therapy (IGRT) are hampered by shading artifacts in the reconstructed images. These artifacts are mainly due to scatter contamination in the projections but also can result from uncorrected beam hardening effects as well as nonlinearities in responses of the amorphous silicon flat panel detectors. While currently, CBCT is mainly used to provide patient geometry information for treatment setup, more demanding applications requiring high-quality CBCT images are under investigation. To tackle these challenges, many CBCT correction algorithms have been proposed; yet, a standard approach still remains unclear. In this work, we propose a shading correction method for CBCT that addresses artifacts from low-frequency projection errors. The method is consistent with the current workflow of radiation therapy. Methods: With much smaller inherent scatter signals and more accurate detectors, diagnostic multidetector CT (MDCT) provides high quality CT images that are routinely used for radiation treatment planning. Using the MDCT image as ''free'' prior information, we first estimate the primary projections in the CBCT scan via forward projection of the spatially registered MDCT data. Since most of the CBCT shading artifacts stem from low-frequency errors in the projections such as scatter, these errors can be accurately estimated by low-pass filtering the difference between the estimated and raw CBCT projections. The error estimates are then subtracted from the raw CBCT projections. Our method is distinct from other published correction methods that use the MDCT image as a prior because it is projection-based and uses limited patient anatomical information from the MDCT image. The merit of CBCT-based treatment monitoring is therefore retained. Results: The proposed method is evaluated using two phantom studies on tabletop systems. On the Catphan(c)600 phantom, our approach reduces the reconstruction error from 348 Hounsfield unit (HU) without correction to 4 HU around the object center after correction, and from 375 HU to 17 HU in the high-contrast regions. In the selected regions of interest (ROIs), the average image contrast is increased by a factor of 3.3. When noise suppression is implemented, the proposed correction substantially improves the contrast-to-noise ratio (CNR) and therefore the visibility of low-contrast objects, as seen in a more challenging pelvis phantom study. Besides a significant improvement in image uniformity, a low-contrast object of {approx}25 HU, which is otherwise buried in the shading artifacts, can be clearly identified after the proposed correction due to a CNR increase of 3.1. Compared to a kernel-based scatter correction method coupled with an analytical beam hardening correction, our approach also shows an overall improved performance with some residual artifacts. Conclusions: By providing effective shading correction, our approach has the potential to improve the accuracy of more advanced CBCT-based clinical applications for IGRT, such as tumor delineation and dose calculation.

60 APPLIED LIFE SCIENCES↗

A novel off-axis scanning method for an enlarged ellipse cone-beam computed tomography field of view

Purpose: Current on-board imaging systems commonly used by modern linear accelerators (LINACs) have a limited field of view (FOV) for a cone-beam CT (CBCT) scan, which is typically less than 50 cm. Consequently, truncation artifacts often occur for large patients. The goal of this work is to investigate a novel method to increase the FOV for current on-board CBCT systems. Methods: When a large patient is scanned with CBCT, any region outside the FOV is only partially sampled within a short range of projection angles, and at other angles no x-ray beams may pass through that region. To increase the sampling rate for the region outside the FOV, we have designed a new source trajectory by shifting the center of rotation during a CBCT scan. This resulted in a reduced sampling rate at the central area and increased sampling rate at the edges. The tradeoff led to a more balanced sampling for an enlarged FOV. An iterative algorithm was also developed to reconstruct the CT image under the new sampling scheme using a compressed sensing technique. Results: The method was validated by numerical simulations mimicking a Varian Trilogy CBCT system, and it was found that artifact-free images could be obtained with the FOV as large as 80 cm. Conclusions: The new CT scanning trajectory can be easily realized under current clinical setup with little modification of the control system, and this can be useful for treating obese patients.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Reconstruction of a cone-beam CT image via forward iterative projection matching

Purpose: To demonstrate the feasibility of reconstructing a cone-beam CT (CBCT) image by deformably altering a prior fan-beam CT (FBCT) image such that it matches the anatomy portrayed in the CBCT projection data set. Methods: A prior FBCT image of the patient is assumed to be available as a source image. A CBCT projection data set is obtained and used as a target image set. A parametrized deformation model is applied to the source FBCT image, digitally reconstructed radiographs (DRRs) that emulate the CBCT projection image geometry are calculated and compared to the target CBCT projection data, and the deformation model parameters are adjusted iteratively until the DRRs optimally match the CBCT projection data set. The resulting deformed FBCT image is hypothesized to be an accurate representation of the patient's anatomy imaged by the CBCT system. The process is demonstrated via numerical simulation. A known deformation is applied to a prior FBCT image and used to create a synthetic set of CBCT target projections. The iterative projection matching process is then applied to reconstruct the deformation represented in the synthetic target projections; the reconstructed deformation is then compared to the known deformation. The sensitivity of the process to the number of projections and the DRR/CBCT projection mismatch is explored by systematically adding noise to and perturbing the contrast of the target projections relative to the iterated source DRRs and by reducing the number of projections. Results: When there is no noise or contrast mismatch in the CBCT projection images, a set of 64 projections allows the known deformed CT image to be reconstructed to within a nRMS error of 1% and the known deformation to within a nRMS error of 7%. A CT image nRMS error of less than 4% is maintained at noise levels up to 3% of the mean projection intensity, at which the deformation error is 13%. At 1% noise level, the number of projections can be reduced to 8 while maintaining CT image and deformation errors of less than 4% and 13%, respectively. The method is sensitive to contrast mismatch between the simulated projections and the target projections when the soft-tissue contrast in the projections is low. Conclusions: By using prior knowledge available in a FBCT image, the authors show that a CBCT image can be iteratively reconstructed from a comparatively small number of projection images, thus saving acquisition time and reducing imaging dose. This will enable more frequent daily imaging during radiation therapy. Because the process preserves the CT numbers of the FBCT image, the resulting 3D image intensities will be more accurate than a CBCT image reconstructed via conventional backprojection methods. Reconstruction errors are insensitive to noise at levels beyond what would typically be found in CBCT projection data, but are sensitive to contrast mismatch errors between the CBCT projection data and the DRRs.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Reconstruction of brachytherapy seed positions and orientations from cone-beam CT x-ray projections via a novel iterative forward projection matching method

Purpose: To generalize and experimentally validate a novel algorithm for reconstructing the 3D pose (position and orientation) of implanted brachytherapy seeds from a set of a few measured 2D cone-beam CT (CBCT) x-ray projections. Methods: The iterative forward projection matching (IFPM) algorithm was generalized to reconstruct the 3D pose, as well as the centroid, of brachytherapy seeds from three to ten measured 2D projections. The gIFPM algorithm finds the set of seed poses that minimizes the sum-of-squared-difference of the pixel-by-pixel intensities between computed and measured autosegmented radiographic projections of the implant. Numerical simulations of clinically realistic brachytherapy seed configurations were performed to demonstrate the proof of principle. An in-house machined brachytherapy phantom, which supports precise specification of seed position and orientation at known values for simulated implant geometries, was used to experimentally validate this algorithm. The phantom was scanned on an ACUITY CBCT digital simulator over a full 660 sinogram projections. Three to ten x-ray images were selected from the full set of CBCT sinogram projections and postprocessed to create binary seed-only images. Results: In the numerical simulations, seed reconstruction position and orientation errors were approximately 0.6 mm and 5 deg., respectively. The physical phantom measurements demonstrated an absolute positional accuracy of (0.78{+-}0.57) mm or less. The {theta} and {phi} angle errors were found to be (5.7{+-}4.9) deg. and (6.0{+-}4.1) deg., respectively, or less when using three projections; with six projections, results were slightly better. The mean registration error was better than 1 mm/6 deg. compared to the measured seed projections. Each test trial converged in 10-20 iterations with computation time of 12-18 min/iteration on a 1 GHz processor. Conclusions: This work describes a novel, accurate, and completely automatic method for reconstructing seed orientations, as well as centroids, from a small number of radiographic projections, in support of intraoperative planning and adaptive replanning. Unlike standard back-projection methods, gIFPM avoids the need to match corresponding seed images on the projections. This algorithm also successfully reconstructs overlapping clustered and highly migrated seeds in the implant. The accuracy of better than 1 mm and 6 deg. demonstrates that gIFPM has the potential to support 2D Task Group 43 calculations in clinical practice.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Feasibility of using cone-beam CT to verify and reposition the optically guided target localization of linear accelerator based stereotactic radiosurgery

Purpose: The optically guided target localization had been developed for linear accelerator based stereotactic radiosurgery (SRS). Unlike the traditional laser localization, the optical guided target localization utilizes a digital system to position patient. Although the system has been proven accurate and robust, it takes away the capability of physicist to directly double check the target position prior to irradiation. Any error from system calibration, data transformation, or head ring position maintenance will not be caught. The purpose of this work is to investigate the possibility of using cone-beam CT (CBCT) to double check the optically guided SRS target localization and reposition the patient. Methods: A SRS quality assurance (QA) phantom was used in the study. The phantom mounted with SRS head frame was scanned by computer tomography (CT) and planned according to the SRS radiation treatment planning process. A target isocenter is defined and transferred to the optically guided target localization system. The phantom was then transported to the linear accelerator room and localized at the initial position agreed by the optically guided target localization system and the CBCT system. Tests were conducted by moving/rotating the phantom to a set of preset offsets and taking CBCT images. Shifts detected by CBCT were compared with the preset offsets. Agreements between them were studied to see how well the CBCT was in discovering the optically guided target localization error. Results: Experiment results demonstrated good agreement between the CBCT detected phantom shift and the preset offset, when the offset is above 1 mm shift or 0.2 degree rotation. Offset less than 1 mm shift or 0.2 degree rotation was not detectable by CBCT. Conclusions: The study concludes that the CBCT is able to discover the optically guided target localization error due to the system calibration or had ring migration. It is a valuable second check tool for SRS target localization quality assurance. The accuracy of CBCT in estimating patient positioning deviation satisfies the SRS procedures with generous tumor size and margin that can tolerate 1 mm or 0.2 degree accuracy. This avoids sending patient home without treatment. CBCT can be neither used as a primary SRS target localization nor can it be used to reposition the patient that cannot tolerate 1 mm shift or 0.2 degree rotation.

60 APPLIED LIFE SCIENCES↗