NASA Astrophysics Data System (ADS)
Bai, Chao-Ying; Huang, Guo-Jiao; Li, Xiao-Ling; Zhou, Bing; Greenhalgh, Stewart
2013-11-01
To overcome the deficiency of some current grid-/cell-based ray tracing algorithms, which are only able to handle first arrivals or primary reflections (or conversions) in anisotropic media, we have extended the functionality of the multistage irregular shortest-path method to 2-D/3-D tilted transversely isotropic (TTI) media. The new approach is able to track multiple transmitted/reflected/converted arrivals composed of any kind of combinations of transmissions, reflections and mode conversions. The basic principle is that the seven parameters (five elastic parameters plus two polar angles defining the tilt of the symmetry axis) of the TTI media are sampled at primary nodes, and the group velocity values at secondary nodes are obtained by tri-linear interpolation of the primary nodes across each cell, from which the group velocities of the three wave modes (qP, qSV and qSH) are calculated. Finally, we conduct grid-/cell-based wave front expansion to trace multiple transmitted/reflected/converted arrivals from one region to the next. The results of calculations in uniform anisotropic media indicate that the numerical results agree with the analytical solutions except in directions of SV-wave triplications, at which only the lowest velocity value is selected at the singularity points by the multistage irregular shortest-path anisotropic ray tracing method. This verifies the accuracy of the methodology. Several simulation results show that the new method is able to efficiently and accurately approximate situations involving continuous velocity variations and undulating discontinuities, and that it is suitable for any combination of multiple transmitted/reflected/converted arrival tracking in TTI media of arbitrary strength and tilt. Crosshole synthetic traveltime tomographic tests have been performed, which highlight the importance of using such code when the medium is distinctly anisotropic.
ERIC Educational Resources Information Center
Gatland, Ian R.
2002-01-01
Proposes a ray tracing approach to thin lens analysis based on a vector form of Snell's law for paraxial rays as an alternative to the usual approach in introductory physics courses. The ray tracing approach accommodates skew rays and thus provides a complete analysis. (Author/KHR)
Robb, P; Pawlowski, B
1990-05-01
The results of measuring the ray trace speed and compilation speed of thirty-nine computers in fifty-seven configurations, ranging from personal computers to super computers, are described. A correlation of ray trace speed has been made with the LINPACK benchmark which allows the ray trace speed to be estimated using LINPACK performance data. The results indicate that the latest generation of workstations, using CPUs based on RISC (Reduced Instruction Set Computer) technology, are as fast or faster than mainframe computers in compute-bound situations. PMID:20563112
NASA Astrophysics Data System (ADS)
Lam, Wai Sze Tiffany
Optical components made of anisotropic materials, such as crystal polarizers and crystal waveplates, are widely used in many complex optical system, such as display systems, microlithography, biomedical imaging and many other optical systems, and induce more complex aberrations than optical components made of isotropic materials. The goal of this dissertation is to accurately simulate the performance of optical systems with anisotropic materials using polarization ray trace. This work extends the polarization ray tracing calculus to incorporate ray tracing through anisotropic materials, including uniaxial, biaxial and optically active materials. The 3D polarization ray tracing calculus is an invaluable tool for analyzing polarization properties of an optical system. The 3x3 polarization ray tracing P matrix developed for anisotropic ray trace assists tracking the 3D polarization transformations along a ray path with series of surfaces in an optical system. To better represent the anisotropic light-matter interactions, the definition of the P matrix is generalized to incorporate not only the polarization change at a refraction/reflection interface, but also the induced optical phase accumulation as light propagates through the anisotropic medium. This enables realistic modeling of crystalline polarization elements, such as crystal waveplates and crystal polarizers. The wavefront and polarization aberrations of these anisotropic components are more complex than those of isotropic optical components and can be evaluated from the resultant P matrix for each eigen-wavefront as well as for the overall image. One incident ray refracting or reflecting into an anisotropic medium produces two eigenpolarizations or eigenmodes propagating in different directions. The associated ray parameters of these modes necessary for the anisotropic ray trace are described in Chapter 2. The algorithms to calculate the P matrix from these ray parameters are described in Chapter 3 for
Ray Tracing with Virtual Objects.
ERIC Educational Resources Information Center
Leinoff, Stuart
1991-01-01
Introduces the method of ray tracing to analyze the refraction or reflection of real or virtual images from multiple optical devices. Discusses ray-tracing techniques for locating images using convex and concave lenses or mirrors. (MDH)
ERIC Educational Resources Information Center
Majewski, Mirek
1997-01-01
Ray tracing is a method that allows the creation of photo-realistic images on a computer. This article describes a shareware ray tracing program called PovRay and includes some ideas on how PovRay can be used in teaching and in 3-D geometry, physics, and other high school and university subjects. (Author/AIM)
NASA Technical Reports Server (NTRS)
Redding, David C.; Breckenridge, William G.
1990-01-01
A new, coordinate-free version of the exact ray-trace equations for optical systems consisting of conic reflecting, refracting and reference surfaces is presented. These equations are differentiated to obtain closed-form optical sensitivity dyadics. For computation, the sensitivities are evaluated in a single global coordinate frame and combined in linearized ray-trace matrix difference equations that propagate the rays and the sensitivities from element to element. One purpose of this analysis is to create optical models that can be directly integrated with models of the instrument structure and control systems for dynamic simulation.
2-D soft x-ray arrays in the EAST.
Chen, Kaiyun; Xu, Liqing; Hu, Liqun; Duan, Yanmin; Li, Xueqin; Yuan, Yi; Mao, Songtao; Sheng, Xiuli; Zhao, Jinlong
2016-06-01
A high spatial and temporal resolution soft x-ray (SXR) imaging diagnostic has been installed in EAST for the study of magnetohydrodynamics activities and core high-Z impurity transport. Up to 122 lines of sight view the poloidal plasma from three directions (two up-down symmetrical horizontal arrays and one vertical array), which renders the diagnostic able to provide detailed tomographic reconstructions under various conditions. Fourier-Bessel method based on flux coordinates was employed for 2-D SXR tomographic reconstruction. Examples of several events measured by SXR diagnostic in EAST are shown, namely the crash patterns of sawtooth, periodical burst of edge localized modes, and the transport of high-Z intrinsic impurities. PMID:27370451
2-D soft x-ray arrays in the EAST
NASA Astrophysics Data System (ADS)
Chen, Kaiyun; Xu, Liqing; Hu, Liqun; Duan, Yanmin; Li, Xueqin; Yuan, Yi; Mao, Songtao; Sheng, Xiuli; Zhao, Jinlong
2016-06-01
A high spatial and temporal resolution soft x-ray (SXR) imaging diagnostic has been installed in EAST for the study of magnetohydrodynamics activities and core high-Z impurity transport. Up to 122 lines of sight view the poloidal plasma from three directions (two up-down symmetrical horizontal arrays and one vertical array), which renders the diagnostic able to provide detailed tomographic reconstructions under various conditions. Fourier-Bessel method based on flux coordinates was employed for 2-D SXR tomographic reconstruction. Examples of several events measured by SXR diagnostic in EAST are shown, namely the crash patterns of sawtooth, periodical burst of edge localized modes, and the transport of high-Z intrinsic impurities.
Seismic ray tracing using linear traveltime interpolation
Asakawa, Eiichi; Kawanaka, Taku )
1993-01-01
A new ray-tracing method called linear traveltime interpolation (LTI) is proposed. This method computes traveltimes and raypaths in a 2D velocity structure more rapidly and accurately than other conventional methods. The LTI method is formulated for a 2D cell model, and calculations of traveltimes and raypaths are carried out only on cell boundaries. Therefore a raypath is considered to be always straight in a cell with uniform velocity. This approach is suitable to tomography analysis. The algorithm of LTI consists of two separate steps: step 1 calculates traveltimes on all cell boundaries; step 2 traces raypaths for all pairs of receivers and the shot. A traveltime at an arbitrary point on a cell boundary is assumed to be linearly interpolated between traveltimes at the adjacent discrete points at which traveltimes were calculated. Fermat's principle is used as the criterion for choosing the correct traveltimes and raypaths from several candidates routinely. The LTI method has been compared numerically with the shooting method and the finite-difference method (FDM) of the eikonal equation. The results show that the LTI method has great advantages of high speed and high accuracy in the calculation of both traveltimes and raypaths. The LTI method can be regarded as an advanced version of the conventional FDM of the eikonal equation because the formulae of FDM are independently derived from LTI. In the process of derivation, it is shown theoretically that LTI is more accurate than FDM. Moreover in the LTI method, they can avoid the numerical instability that occurs in Vidale's method where the velocity changes abruptly.
Real ray tracing in anisotropic viscoelastic media
NASA Astrophysics Data System (ADS)
Vavryčuk, Václav
2008-11-01
Ray tracing equations applicable to smoothly inhomogeneous anisotropic viscoelastic media are derived. The equations produce real rays, in contrast to previous ray-theoretical approaches, which deal with complex rays. The real rays are defined as the solutions of the Hamilton equations, with the Hamiltonian modified for viscoelastic media, and physically correspond to trajectories of high-frequency waves characterized by a real stationary phase. As a consequence, the complex eikonal equation is satisfied only approximately. The ray tracing equations are valid for weakly and moderately attenuating media. The rays are frequency-dependent and must be calculated for each frequency, separately. Solving the ray tracing equations in viscoelastic anisotropy is more time consuming than in elastic anisotropy. The main difficulty is with determining the stationary slowness vector, which is generally complex-valued and inhomogeneous and must be computed at each time step of the ray tracing procedure. In viscoelastic isotropy, the ray tracing equations considerably simplify, because the stationary slowness vector is homogeneous. The computational time for tracing rays in isotropic elastic and viscoelastic media is the same. Using numerical examples, it is shown that ray fields in weakly attenuating media (Q higher than about 30) are almost indistinguishable from those in elastic media. For moderately attenuating anisotropic media (Q between 5-20), the differences in ray fields can be visible and significant.
Reverse ray tracing for transformation optics.
Hu, Chia-Yu; Lin, Chun-Hung
2015-06-29
Ray tracing is an important technique for predicting optical system performance. In the field of transformation optics, the Hamiltonian equations of motion for ray tracing are well known. The numerical solutions to the Hamiltonian equations of motion are affected by the complexities of the inhomogeneous and anisotropic indices of the optical device. Based on our knowledge, no previous work has been conducted on ray tracing for transformation optics with extreme inhomogeneity and anisotropicity. In this study, we present the use of 3D reverse ray tracing in transformation optics. The reverse ray tracing is derived from Fermat's principle based on a sweeping method instead of finding the full solution to ordinary differential equations. The sweeping method is employed to obtain the eikonal function. The wave vectors are then obtained from the gradient of that eikonal function map in the transformed space to acquire the illuminance. Because only the rays in the points of interest have to be traced, the reverse ray tracing provides an efficient approach to investigate the illuminance of a system. This approach is useful in any form of transformation optics where the material property tensor is a symmetric positive definite matrix. The performance and analysis of three transformation optics with inhomogeneous and anisotropic indices are explored. The ray trajectories and illuminances in these demonstration cases are successfully solved by the proposed reverse ray tracing method. PMID:26191770
Validation of Ray Tracing Code Refraction Effects
NASA Technical Reports Server (NTRS)
Heath, Stephanie L.; McAninch, Gerry L.; Smith, Charles D.; Conner, David A.
2008-01-01
NASA's current predictive capabilities using the ray tracing program (RTP) are validated using helicopter noise data taken at Eglin Air Force Base in 2007. By including refractive propagation effects due to wind and temperature, the ray tracing code is able to explain large variations in the data observed during the flight test.
Symplectic ray-tracing: a new approach for nonlinear ray tracings by Hamiltonian dynamics
NASA Astrophysics Data System (ADS)
Satoh, Tetsu R.
2003-05-01
This paper describes a method of symplectic ray tracing for calculating the flows of non-linear dynamical systems. Symplectic ray tracing method traces the path of photons moving along the orbit calculated by using Hamilton's canonical equation. Using this method, we can simulate non-linear dynamical systems with various dimensions, accurate calculation, and quick implementation of scientif visualization system. This paper also demonstrates some visualization results of non-linear dynamical systems computed by using symplectic ray tracing method.
Development of 2-D-MAX-DOAS and retrievals of trace gases and aerosols optical properties
NASA Astrophysics Data System (ADS)
Ortega, Ivan
Air pollution is a major problem worldwide that adversely a_ects human health, impacts ecosystems and climate. In the atmosphere, there are hundreds of important compounds participating in complex atmospheric reactions linked to air quality and climate. Aerosols are relevant because they modify the radiation balance, a_ect clouds, and thus Earth albedo. The amount of aerosol is often characterized by the vertical integral through the entire height of the atmosphere of the logarithm fraction of incident light that is extinguished called Aerosol Optical Depth (AOD). The AOD at 550 nm (AOD550) over land is 0.19 (multi annual global mean), and that over oceans is 0.13. About 43 % of the Earth surface shows AOD550 smaller than 0.1. There is a need for measurement techniques that are optimized to measure aerosol optical properties under low AOD conditions, sample spatial scales that resemble satellite ground-pixels and atmospheric models, and help integrate remote sensing and in-situ observations to obtain optical closure on the effects of aerosols and trace gases in our changing environment. In this work, I present the recent development of the University of Colorado two dimensional (2-D) Multi-AXis Differential Optical Absorption Spectroscopy (2-D-MAX-DOAS) instrument to measure the azimuth and altitude distribution of trace gases and aerosol optical properties simultaneously with a single instrument. The instrument measures solar scattered light from any direction in the sky, including direct sun light in the hyperspectral domain. In Chapter 2, I describe the capabilities of 2-D measurements in the context of retrievals of azimuth distributions of nitrogen dioxide (NO2), formaldehyde (HCHO), and glyoxal (CHOCHO), which are precursors for tropospheric O3 and aerosols. The measurements were carried out during the Multi-Axis DOAS Comparison campaign for Aerosols and Trace gases (MAD-CAT) campaign in Mainz, Germany and show the ability to bridge spatial scales to
NASA Astrophysics Data System (ADS)
Laepple, T.; Heue, K.; Friedeburg, C. V.; Wang, P.; Knab, V.; Pundt, I.
2002-12-01
Tomographic-Differential-Optical-Absorption-Spectroscopy (Tom-DOAS) is a new application of the DOAS method designed to measure 2-3-dimensional concentration fields of different trace gases (e.g. NO2, HCHO, Ozone) in the troposphere. Numerical reconstruction techniques are used to obtain spatially resolved data from the slant column densities provided by DOAS instruments. We discuss the detection of emission plumes by AMAX (Airborne Multi AXis) DOAS Systems which measure sunlight by telescopes pointing in different directions. 2D distributions are reconstructed from slant columns by using airmass factor matrices and inversion techniques. We discuss possibilities and limitations of this technique gained with the use of simulated test fields. Therefore the effect of the parameter choice (e.g. flight track, algorithm changes) and measurement errors is investigated. Further, first results from the Partenavia aircraft measurements over Milano (Italy) during the European FORMAT campaign will be presented.
A three-dimensional ray-tracing code dedicated to x-ray laser amplification simulation
NASA Astrophysics Data System (ADS)
Temporal, M.; Jacquemot, S.; Bonnet, L.; Decoster, A.
2001-04-01
A three-dimensional (3D) ray-tracing code has been developed to simulate the x-ray intensity produced in recent experiments where a silver target was driven by two laser beams. The code is used as a postprocessor of a detailed atomic physics code, which provides emissivities and opacities for inverted transitions. The hydrodynamics of the plasma is calculated with a 1D1/2 hydrocode where transverse profiles of temperature and density follow a self-similar solution. The 3D ray-tracing code accounts for progressive target illumination and calculates the x-ray laser output by solving the eikonal equation. Once 3D paths are determined, a steady-state transport solution is used to calculate the output intensity. The ray-tracing package is discussed first, then the present 3D results are compared with 2D calculations, as well as with collected experimental data.
The method of polarized traces for the 2D Helmholtz equation
NASA Astrophysics Data System (ADS)
Zepeda-Núñez, Leonardo; Demanet, Laurent
2016-03-01
We present a solver for the 2D high-frequency Helmholtz equation in heterogeneous acoustic media, with online parallel complexity that scales optimally as O (N/L), where N is the number of volume unknowns, and L is the number of processors, as long as L grows at most like a small fractional power of N. The solver decomposes the domain into layers, and uses transmission conditions in boundary integral form to explicitly define "polarized traces", i.e., up- and down-going waves sampled at interfaces. Local direct solvers are used in each layer to precompute traces of local Green's functions in an embarrassingly parallel way (the offline part), and incomplete Green's formulas are used to propagate interface data in a sweeping fashion, as a preconditioner inside a GMRES loop (the online part). Adaptive low-rank partitioning of the integral kernels is used to speed up their application to interface data. The method uses second-order finite differences. The complexity scalings are empirical but motivated by an analysis of ranks of off-diagonal blocks of oscillatory integrals. They continue to hold in the context of standard geophysical community models such as BP and Marmousi 2, where convergence occurs in 5 to 10 GMRES iterations. While the parallelism in this paper stems from decomposing the domain, we do not explore the alternative of parallelizing the systems solves with distributed linear algebra routines.
Ray tracing package through a lens system and a spectrometer
Zurro, B.; King, P.W.; Lazarus, E.A.
1980-03-01
To study the light collection optics of the ISX-B two-dimensional (2-D) Thomson scattering system, we have implemented in the Oak Ridge National Laboratory (ORNL) Fusion Energy Division (FED) PDP-10 two computer programs, LENS and SPECT, that trace rays through a lens system and a spectrometer, respectively. The lens package follows the path of any kind of ray (meridional or skew) through a centered optical system formed by an arbitrary number of spherical surfaces. The spectrometer package performs geometrical ray tracing through a Czerney-Turner spectrometer and can be easily modified for studying any other configuration. Contained herein is a description of the procedures followed and a listing of the computer programs.
Improved backward ray tracing with stochastic sampling
NASA Astrophysics Data System (ADS)
Ryu, Seung Taek; Yoon, Kyung-Hyun
1999-03-01
This paper presents a new technique that enhances the diffuse interreflection with the concepts of backward ray tracing. In this research, we have modeled the diffuse rays with the following conditions. First, as the reflection from the diffuse surfaces occurs in all directions, it is impossible to trace all of the reflected rays. We confined the diffuse rays by sampling the spherical angle out of the reflected rays around the normal vector. Second, the traveled distance of reflected energy from the diffuse surface differs according to the object's property, and has a comparatively short reflection distance. Considering the fact that the rays created on the diffuse surfaces affect relatively small area, it is very inefficient to trace all of the sampled diffused rays. Therefore, we set a fixed distance as the critical distance and all the rays beyond this distance are ignored. The result of this research is that as the improved backward ray tracing can model the illumination effects such as the color bleeding effects, we can replace the radiosity algorithm under the limited environment.
A New 2D-Advection-Diffusion Model Simulating Trace Gas Distributions in the Lowermost Stratosphere
NASA Astrophysics Data System (ADS)
Hegglin, M. I.; Brunner, D.; Peter, T.; Wirth, V.; Fischer, H.; Hoor, P.
2004-12-01
Tracer distributions in the lowermost stratosphere are affected by both, transport (advective and non-advective) and in situ sources and sinks. They influence ozone photochemistry, radiative forcing, and heating budgets. In-situ measurements of long-lived species during eight measurement campaigns revealed relatively simple behavior of the tracers in the lowermost stratosphere when represented in an equivalent-latitude versus potential temperature framework. We here present a new 2D-advection-diffusion model that simulates the main transport pathways influencing the tracer distributions in the lowermost stratosphere. The model includes slow diabatic descent of aged stratospheric air and vertical and/or horizontal diffusion across the tropopause and within the lowermost stratosphere. The diffusion coefficients used in the model represent the combined effects of different processes with the potential of mixing tropospheric air into the lowermost stratosphere such as breaking Rossby and gravity waves, deep convection penetrating the tropopause, turbulent diffusion, radiatively driven upwelling etc. They were specified by matching model simulations to observed distributions of long-lived trace gases such as CO and N2O obtained during the project SPURT. The seasonally conducted campaigns allow us to study the seasonal dependency of the diffusion coefficients. Despite its simplicity the model yields a surprisingly good description of the small scale features of the measurements and in particular of the observed tracer gradients at the tropopause. The correlation coefficients between modeled and measured trace gas distributions were up to 0.95. Moreover, mixing across isentropes appears to be more important than mixing across surfaces of constant equivalent latitude (or PV). With the aid of the model, the distribution of the fraction of tropospheric air in the lowermost stratosphere can be determined.
What you need to know about… Tracing the X-ray Trail If you’ve just completed an x-ray, computed tomography (CT), magnetic resonance (MR) Start here! or other diagnostic imaging procedure, you probably want to know when you will ... los rayos X Si acaba de hacerse una radiografía, tomografía ¡Empezar ...
Birefringent Polarization Ray Tracing: Theory and Applications
NASA Astrophysics Data System (ADS)
McClain, Stephen Charles
Birefringent polarization ray tracing is an extension of geometric ray tracing. In addition to calculating ray paths and phases, it also analyzes the state of polarization through birefringent devices. Some systems containing birefringent elements include optical computers, radiometers, optical isolators, bar code scanners, and optical data storage systems. This dissertation derives explicit algorithms for polarization ray tracing through anisotropic media, optically active media, and anisotropic optically active media, such as quartz. The objective was to go beyond the electromagnetic relations to establish algorithms in standard ray tracing format, ready for direct inclusion into lens design software. The algorithms, derived from Maxwell's equations, constitutive relations, and boundary conditions, calculate the wavevector, ray vector, optical path length, refractive index, and polarization state of a ray. Generalized Fresnel relations govern the division of energy at each interface into two transmitted and two reflected modes. The algorithms are applied to calculate the polarization aberrations of a variety of birefringent devices. In particular, it is established that the polarization properties of quartz vary significantly (>20%) over angles of only 5 degrees. This limits the useful field of view of quartz devices. Field of view aberrations of birefringent elements can critically affect the performance of optical systems. Also, design guidelines are presented for pseudodepolarizers. These devices spatially scramble the polarization. Inserted into an instrument, a depolarizer negates the polarization sensitivity of the elements which follow it. Presented in detail is the design and analysis of a depolarizer for use in a spectrometer on NASA's Earth Observing System (EOS).
Birefringent polarization ray tracing: Theory and applications
NASA Astrophysics Data System (ADS)
McClain, Stephen Charles
1992-06-01
Birefringent polarization ray tracing is an extension of geometric ray tracing. In addition to calculating ray paths and phases, it also analyzes the state of polarization through birefringent devices. Some systems containing birefringent elements include optical computers, radiometers, optical isolators, bar code scanners, and optical data storage systems. This dissertation derives explicit algorithms for polarization ray tracing through anisotropic media, optically active media, and anisotropic optically active media, such as quartz. The objective was to go beyond the electromagnetic relations to establish algorithms in standard ray tracing format, ready for direct inclusion into lens design software. The algorithms, derived from Maxwell's equations, constitutive relations, and boundary conditions, calculate the wavevector, ray vector, optical path length, refractive index, and polarization state of a ray. Generalized Fresnel relations govern the division of energy at each interface into two transmitted and two reflected modes. The algorithms are applied to calculate the polarization aberrations of a variety of birefringent devices. In particular, it is established that the polarization properties of quartz vary significantly (greater than 20 percent) over angles of only 5 degrees. This limits the useful field of view of quartz devices. Field of view aberrations of birefringent elements can critically affect the performance of optical systems. Also, design guidelines are presented for pseudodepolarizers. These devices spatially scramble the polarization. Inserted into an instrument, a depolarizer negates the polarization sensitivity of the elements which follow it. Presented in detail is the design and analysis of a depolarizer for use in a spectrometer on NASA's Earth Observing System (EOS).
Zhou, Zhi; Liu, Xiaoxiao; Long, Brian; Peng, Hanchuan
2016-01-01
Efficient and accurate digital reconstruction of neurons from large-scale 3D microscopic images remains a challenge in neuroscience. We propose a new automatic 3D neuron reconstruction algorithm, TReMAP, which utilizes 3D Virtual Finger (a reverse-mapping technique) to detect 3D neuron structures based on tracing results on 2D projection planes. Our fully automatic tracing strategy achieves close performance with the state-of-the-art neuron tracing algorithms, with the crucial advantage of efficient computation (much less memory consumption and parallel computation) for large-scale images. PMID:26306866
AXAF FITS standard for ray trace interchange
NASA Astrophysics Data System (ADS)
Hsieh, Paul F.
1993-07-01
A standard data format for the archival and transport of x-ray events generated by ray trace models is described. Upon review and acceptance by the Advanced X-ray Astrophysics Facility (AXAF) Software Systems Working Group (SSWG), this standard shall become the official AXAF data format for ray trace events. The Flexible Image Transport System (FITS) is well suited for the purposes of the standard and was selected to be the basis of the standard. FITS is both flexible and efficient and is also widely used within the astronomical community for storage and transfer of data. In addition, software to read and write FITS format files are widely available. In selecting quantities to be included within the ray trace standard, the AXAF Mission Support team, Science Instruments team, and the other contractor teams were surveyed. From the results of this survey, the following requirements were established: (1) for the scientific needs, each photon should have associated with it: position, direction, energy, and statistical weight; the standard must also accommodate path length (relative phase), and polarization. (2) a unique photon identifier is necessary for bookkeeping purposes; (3) a log of individuals, organizations, and software packages that have modified the data must be maintained in order to create an audit trail; (4) a mechanism for extensions to the basic kernel should be provided; and (5) the ray trace standard should integrate with future AXAF data product standards.
AXAF FITS standard for ray trace interchange
NASA Technical Reports Server (NTRS)
Hsieh, Paul F.
1993-01-01
A standard data format for the archival and transport of x-ray events generated by ray trace models is described. Upon review and acceptance by the Advanced X-ray Astrophysics Facility (AXAF) Software Systems Working Group (SSWG), this standard shall become the official AXAF data format for ray trace events. The Flexible Image Transport System (FITS) is well suited for the purposes of the standard and was selected to be the basis of the standard. FITS is both flexible and efficient and is also widely used within the astronomical community for storage and transfer of data. In addition, software to read and write FITS format files are widely available. In selecting quantities to be included within the ray trace standard, the AXAF Mission Support team, Science Instruments team, and the other contractor teams were surveyed. From the results of this survey, the following requirements were established: (1) for the scientific needs, each photon should have associated with it: position, direction, energy, and statistical weight; the standard must also accommodate path length (relative phase), and polarization. (2) a unique photon identifier is necessary for bookkeeping purposes; (3) a log of individuals, organizations, and software packages that have modified the data must be maintained in order to create an audit trail; (4) a mechanism for extensions to the basic kernel should be provided; and (5) the ray trace standard should integrate with future AXAF data product standards.
Ray tracing through progressive ophthalmic lenses
NASA Astrophysics Data System (ADS)
Bourdoncle, Bernard; Chauveau, J. P.; Mercier, Jean-Louis M.
1991-01-01
Ray-tracing through Progressive Addition Lenses (PAL) has been performed. PAL is a deep non rotationally symmetric asp1ric lens used for the compensation of presbyopia. PAL and its mathematical model are presented. The special features of the ray-tracing program due to the model of the lens plus eye system are detailed. Typical results are presented showing in particular that computing conditions of contour-plots of power and astigmatism must be very strict and that coma must be taken into account for precise measurements of PAL. 1.
Ray tracing study for non-imaging daylight collectors
Wittkopf, Stephen; Oliver Grobe, Lars; Geisler-Moroder, David; Compagnon, Raphael; Kaempf, Jerome; Linhart, Friedrich; Scartezzini, Jean-Louis
2010-06-15
This paper presents a novel method to study how well non-imaging daylight collectors pipe diffuse daylight into long horizontal funnels for illuminating deep buildings. Forward ray tracing is used to derive luminous intensity distributions curves (LIDC) of such collectors centered in an arc-shaped light source representing daylight. New photometric characteristics such as 2D flux, angular spread and horizontal offset are introduced as a function of such LIDC. They are applied for quantifying and thus comparing different collector contours. (author)
A Fast Ray-Tracing Using Bounding Spheres and Frustum Rays for Dynamic Scene Rendering
NASA Astrophysics Data System (ADS)
Suzuki, Ken-Ichi; Kaeriyama, Yoshiyuki; Komatsu, Kazuhiko; Egawa, Ryusuke; Ohba, Nobuyuki; Kobayashi, Hiroaki
Ray tracing is one of the most popular techniques for generating photo-realistic images. Extensive research and development work has made interactive static scene rendering realistic. This paper deals with interactive dynamic scene rendering in which not only the eye point but also the objects in the scene change their 3D locations every frame. In order to realize interactive dynamic scene rendering, RTRPS (Ray Tracing based on Ray Plane and Bounding Sphere), which utilizes the coherency in rays, objects, and grouped-rays, is introduced. RTRPS uses bounding spheres as the spatial data structure which utilizes the coherency in objects. By using bounding spheres, RTRPS can ignore the rotation of moving objects within a sphere, and shorten the update time between frames. RTRPS utilizes the coherency in rays by merging rays into a ray-plane, assuming that the secondary rays and shadow rays are shot through an aligned grid. Since a pair of ray-planes shares an original ray, the intersection for the ray can be completed using the coherency in the ray-planes. Because of the three kinds of coherency, RTRPS can significantly reduce the number of intersection tests for ray tracing. Further acceleration techniques for ray-plane-sphere and ray-triangle intersection are also presented. A parallel projection technique converts a 3D vector inner product operation into a 2D operation and reduces the number of floating point operations. Techniques based on frustum culling and binary-tree structured ray-planes optimize the order of intersection tests between ray-planes and a sphere, resulting in 50% to 90% reduction of intersection tests. Two ray-triangle intersection techniques are also introduced, which are effective when a large number of rays are packed into a ray-plane. Our performance evaluations indicate that RTRPS gives 13 to 392 times speed up in comparison with a ray tracing algorithm without organized rays and spheres. We found out that RTRPS also provides competitive
Ray tracing on a networked processor array
NASA Astrophysics Data System (ADS)
Yang, Jungsook; Lee, Seung Eun; Chen, Chunyi; Bagherzadeh, Nader
2010-10-01
As computation costs increase to meet design requirements for computation-intensive graphics applications on today's embedded systems, the pressure to develop high-performance parallel processors on a chip will increase. Acceleration of the ray tracing computation has become a major issue as the computer graphics industry demands for rendering realistic images. Network-on-chip (NoC) techniques that interconnect multiple processing elements with routers are the solution for reducing computation time and power consumption by parallel processing on a chip. It is also essential to meet the scalability and complexity challenges for system-on-chip (SoC). In this article, we describe a parallel ray tracing application mapping on a mesh-based multicore NoC architecture. We describe an optimised ray tracing kernel and parallelisation strategies, varying the workload distribution statically and dynamically. In this work, we present results and timing performance of our parallel ray tracing application on a NoC, which are obtained through our cycle accurate multicore NoC simulator. Using a dynamic scheduling load balancing technique, we achieved a maximum speedup multiplier of 35.97 on an 8 × 8 networked processor array using a NoC as the interconnect.
Tracing Rays In Laser-Fringe Anemometers
NASA Technical Reports Server (NTRS)
Owen, Karl
1989-01-01
"OPTMAIN" is simple ray-tracing computer code developed to quantify refractive effects that result when laser-fringe anemometer used to observe flows through window. Code calculates changes for four different types of windows: flat-plate windows, simple cylindrical windows, "general" axisymmetric windows, and smooth general-surface windows. Written in FORTRAN IV.
Ray Traces Through Unsteady Jet Turbulence
NASA Technical Reports Server (NTRS)
Freund, J. B.; Fleischman, T. G.
2002-01-01
Results of an ongoing effort to quantify the role turbulence in scattering sound in jets are reported. Using a direct numerical simulation database to provide the flow data, ray paths traced through the mean flow are compared with those traced through the actual time evolving turbulent flow. Significant scattering by the turbulence is observed. The most notable effect is that upstream traveling waves that are trapped in the potential core by the mean flow, which acts as a wave guide, easily escape in the turbulent flow. A crude statistical estimate based on ray number density suggests that directivity is modified by the turbulence, but no rigorous treatment of non-uniformities in the high-frequency approximation is attempted.
Testing the ray-tracing code GYOTO
NASA Astrophysics Data System (ADS)
Grould, M.; Paumard, T.; Perrin, G.
2015-12-01
In the next few years, the near-infrared interferometer GRAVITY will observe the Galactic Center. Astrometric data will be obtained with an expected accuracy of 10 μas. In order to analyze those future data, we have developed a code named GYOTO to compute orbits and ray-trace images. We want to assess the validity and accuracy of GYOTO in a variety of contexts, in particular for stellar astrometry in the Galactic Center.
Special relativistic visualization by local ray tracing.
Müller, Thomas; Grottel, Sebastian; Weiskopf, Daniel
2010-01-01
Special relativistic visualization offers the possibility of experiencing the optical effects of traveling near the speed of light, including apparent geometric distortions as well as Doppler and searchlight effects. Early high-quality computer graphics images of relativistic scenes were created using offline, computationally expensive CPU-side 4D ray tracing. Alternate approaches such as image-based rendering and polygon-distortion methods are able to achieve interactivity, but exhibit inferior visual quality due to sampling artifacts. In this paper, we introduce a hybrid rendering technique based on polygon distortion and local ray tracing that facilitates interactive high-quality visualization of multiple objects moving at relativistic speeds in arbitrary directions. The method starts by calculating tight image-space footprints for the apparent triangles of the 3D scene objects. The final image is generated using a single image-space ray tracing step incorporating Doppler and searchlight effects. Our implementation uses GPU shader programming and hardware texture filtering to achieve high rendering speed. PMID:20975164
Nonlinear ray tracing for vessel enhanced visualization
NASA Astrophysics Data System (ADS)
Qiu, Feng; Hong, Wei
2012-02-01
3D visualization of angiography data is an important preprocessing step in diagnosis of vascular disease. This paper describes an efficient volume rendering method to emphasize feature-rich region (or focus) in the 3D angiography data. The method takes the input 3D angiography data and computes the focus with user specification or certain feature extraction algorithms. Then, a distance map is constructed based on the description of the focused region(s). While rendering the 3D angiography data, the nonlinear ray tracing method is used and the gradient of the distance volume is applied to guide ray marching. In the result image, the focused region(s) appears larger than in the normal ray-casting image, while the context (other regions of the volume) can be still preserved in the image (maybe displayed in a shrink size). This method avoids deforming the original volume to magnify focus regions, which is expensive to compute, thus improves the performance.
2D and 3D X-Ray Structural Microscopy Using Submicron-Resolution Laue Microdiffraction
Budai, John D.; Yang, Wenge; Larson, Bennett C.; Tischler, Jonathan Z.; Liu, Wenjun; Ice, Gene E.
2010-11-10
We have developed a scanning, polychromatic x-ray microscopy technique with submicron spatial resolution at the Advanced Photon Source. In this technique, white undulator radiation is focused to submicron diameter using elliptical mirrors. Laue diffraction patterns scattered from the sample are collected with an area detector and then analyzed to obtain the local crystal structure, lattice orientation, and strain tensor. These new microdiffraction capabilities have enabled both 2D and 3D structural studies of materials on mesoscopic length-scales of tenths-to-hundreds of microns. For thin samples such as deposited films, 2D structural maps are obtained by step-scanning the area of interest. For example, 2D x-ray microscopy has been applied in studies of the epitaxial growth of oxide films. For bulk samples, a 3D differential-aperture x-ray microscopy technique has been developed that yields the full diffraction information from each submicron volume element. The capabilities of 3D x-ray microscopy are demonstrated here with measurements of grain orientations and grain boundary motion in polycrystalline aluminum during 3D thermal grain growth. X-ray microscopy provides the needed, direct link between the experimentally measured 3D microstructural evolution and the results of theory and modeling of materials processes on mesoscopic length scales.
NASA Astrophysics Data System (ADS)
Clausing, Eric; Vielhauer, Claus
2015-03-01
Locksmith forensics is an important and very challenging part of classic crime scene forensics. In prior work, we propose a partial transfer to the digital domain, to effectively support forensic experts and present approaches for a full process chain consisting of five steps: Trace positioning, 2D/3D acquisition with a confocal 3D laser scanning microscope, detection by segmentation, trace type determination, and determination of the opening method. In particular the step of trace segmentation on high-resolution 3D surfaces thereby turned out to be the part most difficult to implement. The reason for that is the highly structured and complex surfaces to be analyzed. These surfaces are cluttered with a high number of toolmarks, which overlap and distort each other. In Clausing et al., we present an improved approach for a reliable segmentation of relevant trace regions but without the possibility of separating single traces out of segmented trace regions. However, in our past research, especially features based on shape and dimension turned out to be highly relevant for a fully automated analysis and interpretation. In this paper, we consequently propose an approach for this separation. To achieve this goal, we use our segmentation approach and expand it with a combination of the watershed algorithm with a graph-based analysis. Found sub-regions are compared based on their surface character and are connected or divided depending on their similarity. We evaluate our approach with a test set of about 1,300 single traces on the exemplary locking cylinder component 'key pin' and thereby are able of showing the high suitability of our approach.
Kostro, André; Geiger, Mario; Scartezzini, Jean-Louis; Schüler, Andreas
2016-07-01
Advanced optical ray tracing software, CFSpro, was developed for the study and optimization of complex fenestration systems (CFSs). Using an algorithm mixing 2D and 3D approaches, accurate computation of large numbers of rays in extruded geometries can be performed and visualized in real time. A thin film model was included to assess the spectral control provided by coatings. In this paper, the ray tracing model is described and validated. A novel glazing, engineered with this simulation tool, is presented. It combines the functions of daylight provision, glare protection, and seasonal thermal control while conserving a view to the outside at near normal incidence. PMID:27409200
Magnetospheric ray tracing studies. [Jupiter's decametric radiation
NASA Technical Reports Server (NTRS)
Six, N. F.
1982-01-01
Using a model of Jupiter's magnetized plasma environment, radiation raypaths were calculated with a three-dimension ray tracing program. It is assumed that energetic particles produce the emission in the planet's auroral zone at frequencies just above the electron gyrofrequencies. This radiation is generated in narrow sheets defined by the angle of a ray with respect to the magnetic field line. By specifying the source position: latitude, longitude, and radial distance from the planet, signatures in the spectrum of frequency versus time seen by Voyager 1 and 2 were duplicated. The frequency range and the curvature of the decametric arcs in these dynamic spectra are the result of the geometry of the radiation sheets (imposed by the plasma and by the B-field) and illumination of Voyager 1 and 2 as the rotating magnetosphere mimics a pulsar.
Application of ray tracing in radiation heat transfer
NASA Technical Reports Server (NTRS)
Baumeister, Joseph F.
1993-01-01
This collection of presentation figures displays the capabilities of ray tracing for radiation propagation calculations as compared to an analytical approach. The goal is to introduce the terminology and solution process used in ray tracing, and provide insight into radiation heat transfer principles and analysis tools. A thermal analysis working environment is introduced that solves demanding radiation heat transfer problems based on ray tracing. This information may serve as a reference for designing and building ones own analysis environment.
Powerful scriptable ray tracing package xrt
NASA Astrophysics Data System (ADS)
Klementiev, Konstantin; Chernikov, Roman
2014-09-01
We present an open source python based ray tracing tool that offers several useful features in graphical presentation, material properties, advanced calculations of synchrotron sources, implementation of diffractive and refractive elements, complex (also closed) surfaces and multiprocessing. The package has many usage examples which are supplied together with the code and visualized on its web page. We exemplify the present version by modeling (i) a curved crystal analyzer, (ii) a quarter wave plate, (iii) Bragg-Fresnel optics and (iv) multiple reflective and non-sequential optics (polycapillary). The present version implements the use of OpenCL framework that executes calculations on both CPUs and GPUs. Currently, the calculations of an undulator source on a GPU show a gain of about two orders of magnitude in computing time. The development version is successful in modelling the wavefront propagation. Two examples of diffraction on a plane mirror and a plane blazed grating are given for a beam with a finite energy band.
Probing transverse coherence of x-ray beam with 2-D phase grating interferometer
Marathe, Shashidhara; Shi, Xianbo; Wojcik, Michael J.; Kujala, Naresh G.; Divan, Ralu; Mancini, Derrick C.; Macrander, Albert T.; Assoufid, Lahsen
2014-01-01
Transverse coherence of the x-ray beam from a bending magnet source was studied along multiple directions using a 2-D π/2 phase grating by measuring interferogram visibilities at different distances behind the grating. These measurements suggest that the preferred measuring orientation of a 2-D checkerboard grating is along the diagonal directions of the square blocks, where the interferograms have higher visibility and are not sensitive to the deviation of the duty cycle of the grating period. These observations are verified by thorough wavefront propagation simulations. The accuracy of the measured coherence values was also validated by the simulation and analytical results obtained from the source parameters. In addition, capability of the technique in probing spatially resolved local transverse coherence is demonstrated. PMID:24977503
2D electron temperature diagnostic using soft x-ray imaging technique
Nishimura, K. Sanpei, A. Tanaka, H.; Ishii, G.; Kodera, R.; Ueba, R.; Himura, H.; Masamune, S.; Ohdachi, S.; Mizuguchi, N.
2014-03-15
We have developed a two-dimensional (2D) electron temperature (T{sub e}) diagnostic system for thermal structure studies in a low-aspect-ratio reversed field pinch (RFP). The system consists of a soft x-ray (SXR) camera with two pin holes for two-kinds of absorber foils, combined with a high-speed camera. Two SXR images with almost the same viewing area are formed through different absorber foils on a single micro-channel plate (MCP). A 2D T{sub e} image can then be obtained by calculating the intensity ratio for each element of the images. We have succeeded in distinguishing T{sub e} image in quasi-single helicity (QSH) from that in multi-helicity (MH) RFP states, where the former is characterized by concentrated magnetic fluctuation spectrum and the latter, by broad spectrum of edge magnetic fluctuations.
Probing transverse coherence of x-ray beam with 2-D phase grating interferometer.
Marathe, Shashidhara; Shi, Xianbo; Wojcik, Michael J; Kujala, Naresh G; Divan, Ralu; Mancini, Derrick C; Macrander, Albert T; Assoufid, Lahsen
2014-06-16
Transverse coherence of the x-ray beam from a bending magnet source was studied along multiple directions using a 2-D π/2 phase grating by measuring interferogram visibilities at different distances behind the grating. These measurements suggest that the preferred measuring orientation of a 2-D checkerboard grating is along the diagonal directions of the square blocks, where the interferograms have higher visibility and are not sensitive to the deviation of the duty cycle of the grating period. These observations are verified by thorough wavefront propagation simulations. The accuracy of the measured coherence values was also validated by the simulation and analytical results obtained from the source parameters. In addition, capability of the technique in probing spatially resolved local transverse coherence is demonstrated. PMID:24977503
2D X-ray scanner and its uses in laboratory reservoir characterization measurements
Maloney, D.; Doggett, K.
1997-08-01
X-ray techniques are used in petroleum laboratories for a variety of reservoir characterization measurements. This paper describes the configuration of a 2D X-ray scanner and many of the ways in which it simplifies and improves accuracy`s of laboratory measurements. Linear X-ray scanners are most often used to provide descriptions of fluid saturations within core plugs during flow tests. We configured our linear scanner for both horizontal and vertical movement. Samples can be scanned horizontally, vertically, or according to horizontal and vertical grids. X-ray measurements are fast, allowing measurements of two- and three-phase fluid saturations during both steady- and unsteady-state flow processes. Rock samples can be scanned while they are subjected to stress, pore pressure, and temperature conditions simulating those of a petroleum reservoir. Many types of measurements are possible by selecting appropriate X-ray power settings, dopes, filters, and collimator configurations. The scanner has been used for a variety of applications besides fluid saturation measurements. It is useful for measuring porosity distributions in rocks, concentrations of X-ray dopes within flow streams during tracer tests, gap widths in fracture flow cells, fluid interface levels in PVT cells and fluid separators, and other features and phenomena.
Ray tracing in nuclear-pumped flowing gas lasers
Mat'ev, V Yu
2003-06-30
The ray tracing in the resonators of a nuclear-pumped flowing gas lasers is considered. The refractive index profile of the medium in a direction perpendicular to the optical axis in such lasers can be considered parabolic, but the steepness of the parabola is quite nonuniform along the ray trace, and the resonator stability condition (the absolute value of the ray matrix trace for a single trip of the ray in the resonator is smaller than two) is not sufficient to confine the ray within the resonator after a large number of trips. (lasers)
Snapshot 2D tomography via coded aperture x-ray scatter imaging
MacCabe, Kenneth P.; Holmgren, Andrew D.; Tornai, Martin P.; Brady, David J.
2015-01-01
This paper describes a fan beam coded aperture x-ray scatter imaging system which acquires a tomographic image from each snapshot. This technique exploits cylindrical symmetry of the scattering cross section to avoid the scanning motion typically required by projection tomography. We use a coded aperture with a harmonic dependence to determine range, and a shift code to determine cross-range. Here we use a forward-scatter configuration to image 2D objects and use serial exposures to acquire tomographic video of motion within a plane. Our reconstruction algorithm also estimates the angular dependence of the scattered radiance, a step toward materials imaging and identification. PMID:23842254
Schmidt-Arras, Dirk; Leclercq, Olivier; Gherardini, Pier Federico; Helmer-Citterich, Manuela; Faigle, Wolfgang; Loew, Damarys; Späth, Gerald F
2011-08-24
The protozoan parasite Leishmania donovani undergoes various developmental transitions during its infectious cycle that are triggered by environmental signals encountered inside insect and vertebrate hosts. Intracellular differentiation of the pathogenic amastigote stage is induced by pH and temperature shifts that affect protein kinase activities and downstream protein phosphorylation. Identification of parasite proteins with phosphotransferase activity during intracellular infection may reveal new targets for pharmacological intervention. Here we describe an improved protocol to trace this activity in L. donovani extracts at high resolution combining in-gel kinase assay and two-dimensional gel electrophoresis. This 2D procedure allowed us to identify proteins that are associated with amastigote ATP-binding, ATPase, and phosphotransferase activities. The 2D in-gel kinase assay, in combination with recombinant phospho-protein substrates previously identified by phospho-proteomics analyses, provides a novel tool to establish specific protein kinase-substrate relationships thus improving our understanding of Leishmania signal transduction with relevance for future drug development. PMID:21443974
Some extensions of the Booker method of ray tracing
NASA Astrophysics Data System (ADS)
Budden, K. G.
1989-10-01
The Booker method of ray tracing is useful for tracing the paths of radio rays in a plane stratified magneto-plasma, and a brief summary is given. In its simplest form it is used for tracing a single ray, without taking account of the configuration of neighboring rays. Two extensions of the method are discussed. The first is the effect on the signal amplitude of the divergence or convergence of neighboring rays in a thin ray pencil. The second is the technique of complex rays and complex space, which are especially useful when electron collisions introduce attenuation of the waves. Some typical results for both extensions are presented. This paper is mainly a tutorial paper, but some new results are given, including a method for dealing with a singularity in the equations for the ray divergence effect.
Characterization of a 2D soft x-ray tomography camera with discrimination in energy bandsa)
NASA Astrophysics Data System (ADS)
Romano, A.; Pacella, D.; Mazon, D.; Murtas, F.; Malard, P.; Gabellieri, L.; Tilia, B.; Piergotti, V.; Corradi, G.
2010-10-01
A gas detector with a 2D pixel readout is proposed for a future soft x-ray (SXR) tomography with discrimination in energy bands separately per pixel. The detector has three gas electron multiplier foils for the electron amplification and it offers the advantage, compared with the single stage, to be less sensitive to neutrons and gammas. The energy resolution and the detection efficiency of the detector have been accurately studied in the laboratory with continuous SXR spectra produced by an electronic tube and line emissions produced by fluorescence (K, Fe, and Mo) in the range of 3-17 keV. The front-end electronics, working in photon counting mode with a selectable threshold for pulse discrimination, is optimized for high rates. The distribution of the pulse amplitude has been indirectly derived by means of scans of the threshold. Scans in detector gain have also been performed to assess the capability of selecting different energy ranges.
Characterization of a 2D soft x-ray tomography camera with discrimination in energy bands
Romano, A.; Pacella, D.; Gabellieri, L.; Tilia, B.; Piergotti, V.; Mazon, D.; Malard, P.
2010-10-15
A gas detector with a 2D pixel readout is proposed for a future soft x-ray (SXR) tomography with discrimination in energy bands separately per pixel. The detector has three gas electron multiplier foils for the electron amplification and it offers the advantage, compared with the single stage, to be less sensitive to neutrons and gammas. The energy resolution and the detection efficiency of the detector have been accurately studied in the laboratory with continuous SXR spectra produced by an electronic tube and line emissions produced by fluorescence (K, Fe, and Mo) in the range of 3-17 keV. The front-end electronics, working in photon counting mode with a selectable threshold for pulse discrimination, is optimized for high rates. The distribution of the pulse amplitude has been indirectly derived by means of scans of the threshold. Scans in detector gain have also been performed to assess the capability of selecting different energy ranges.
The vectorization of a ray tracing program for image generation
NASA Technical Reports Server (NTRS)
Plunkett, D. J.; Cychosz, J. M.; Bailey, M. J.
1984-01-01
Ray tracing is a widely used method for producing realistic computer generated images. Ray tracing involves firing an imaginary ray from a view point, through a point on an image plane, into a three dimensional scene. The intersections of the ray with the objects in the scene determines what is visible at the point on the image plane. This process must be repeated many times, once for each point (commonly called a pixel) in the image plane. A typical image contains more than a million pixels making this process computationally expensive. A traditional ray tracing program processes one ray at a time. In such a serial approach, as much as ninety percent of the execution time is spent computing the intersection of a ray with the surface in the scene. With the CYBER 205, many rays can be intersected with all the bodies im the scene with a single series of vector operations. Vectorization of this intersection process results in large decreases in computation time. The CADLAB's interest in ray tracing stems from the need to produce realistic images of mechanical parts. A high quality image of a part during the design process can increase the productivity of the designer by helping him visualize the results of his work. To be useful in the design process, these images must be produced in a reasonable amount of time. This discussion will explain how the ray tracing process was vectorized and gives examples of the images obtained.
Fracture network evaluation program (FraNEP): A software for analyzing 2D fracture trace-line maps
NASA Astrophysics Data System (ADS)
Zeeb, Conny; Gomez-Rivas, Enrique; Bons, Paul D.; Virgo, Simon; Blum, Philipp
2013-10-01
Fractures, such as joints, faults and veins, strongly influence the transport of fluids through rocks by either enhancing or inhibiting flow. Techniques used for the automatic detection of lineaments from satellite images and aerial photographs, LIDAR technologies and borehole televiewers significantly enhanced data acquisition. The analysis of such data is often performed manually or with different analysis software. Here we present a novel program for the analysis of 2D fracture networks called FraNEP (Fracture Network Evaluation Program). The program was developed using Visual Basic for Applications in Microsoft Excel™ and combines features from different existing software and characterization techniques. The main novelty of FraNEP is the possibility to analyse trace-line maps of fracture networks applying the (1) scanline sampling, (2) window sampling or (3) circular scanline and window method, without the need of switching programs. Additionally, binning problems are avoided by using cumulative distributions, rather than probability density functions. FraNEP is a time-efficient tool for the characterisation of fracture network parameters, such as density, intensity and mean length. Furthermore, fracture strikes can be visualized using rose diagrams and a fitting routine evaluates the distribution of fracture lengths. As an example of its application, we use FraNEP to analyse a case study of lineament data from a satellite image of the Oman Mountains.
ROBAST: ROOT-based ray-tracing library for cosmic-ray telescopes
NASA Astrophysics Data System (ADS)
Okumura, Akira
2016-03-01
ROBAST (ROOT-based simulator for ray tracing) is a non-sequential ray-tracing simulation library developed for wide use in optical simulations of gamma-ray and cosmic-ray telescopes. The library is written in C++ and fully utilizes the geometry library of the ROOT analysis framework, and can build the complex optics geometries typically used in cosmic ray experiments and ground-based gamma-ray telescopes.
Seismic wavefield propagation in 2D anisotropic media: Ray theory versus wave-equation simulation
NASA Astrophysics Data System (ADS)
Bai, Chao-ying; Hu, Guang-yi; Zhang, Yan-teng; Li, Zhong-sheng
2014-05-01
Despite the ray theory that is based on the high frequency assumption of the elastic wave-equation, the ray theory and the wave-equation simulation methods should be mutually proof of each other and hence jointly developed, but in fact parallel independent progressively. For this reason, in this paper we try an alternative way to mutually verify and test the computational accuracy and the solution correctness of both the ray theory (the multistage irregular shortest-path method) and the wave-equation simulation method (both the staggered finite difference method and the pseudo-spectral method) in anisotropic VTI and TTI media. Through the analysis and comparison of wavefield snapshot, common source gather profile and synthetic seismogram, it is able not only to verify the accuracy and correctness of each of the methods at least for kinematic features, but also to thoroughly understand the kinematic and dynamic features of the wave propagation in anisotropic media. The results show that both the staggered finite difference method and the pseudo-spectral method are able to yield the same results even for complex anisotropic media (such as a fault model); the multistage irregular shortest-path method is capable of predicting similar kinematic features as the wave-equation simulation method does, which can be used to mutually test each other for methodology accuracy and solution correctness. In addition, with the aid of the ray tracing results, it is easy to identify the multi-phases (or multiples) in the wavefield snapshot, common source point gather seismic section and synthetic seismogram predicted by the wave-equation simulation method, which is a key issue for later seismic application.
Light ray tracing through a leaf cross section
NASA Technical Reports Server (NTRS)
Kumar, R.; Silva, L. F.
1973-01-01
A light ray, incident at about 5 deg to the normal, is geometrically plotted through the drawing of the cross section of a soybean leaf using Fresnel's equations and Snell's law. The optical mediums of the leaf considered for ray tracing are: air, cell sap, chloroplast, and cell wall. The ray is also drawn through the same leaf cross section with cell wall and air as the only optical mediums. The values of the reflection and transmission found from the ray tracing tests agree closely with the experimental results obtained using a Beckman Dk-2A Spectroreflector.
Light ray tracing through a leaf cross section
NASA Technical Reports Server (NTRS)
Kumar, R.; Silva, L.
1973-01-01
A light ray, incident at about 5 deg to the normal, is geometrically plotted through the drawing of the cross section of a soybean leaf using Fresnel's equations and Snell's law. The optical mediums of the leaf considered for ray tracing are air, cell sap, chloroplast, and cell wall. The above ray is also drawn through the same leaf cross section considering cell wall and air as the only optical mediums. The values of the reflection and transmission found from ray tracing agree closely with the experimental results obtained using a Beckman DK-2A spectroreflectometer.
Registration of 2D x-ray images to 3D MRI by generating pseudo-CT data
NASA Astrophysics Data System (ADS)
van der Bom, M. J.; Pluim, J. P. W.; Gounis, M. J.; van de Kraats, E. B.; Sprinkhuizen, S. M.; Timmer, J.; Homan, R.; Bartels, L. W.
2011-02-01
Spatial and soft tissue information provided by magnetic resonance imaging can be very valuable during image-guided procedures, where usually only real-time two-dimensional (2D) x-ray images are available. Registration of 2D x-ray images to three-dimensional (3D) magnetic resonance imaging (MRI) data, acquired prior to the procedure, can provide optimal information to guide the procedure. However, registering x-ray images to MRI data is not a trivial task because of their fundamental difference in tissue contrast. This paper presents a technique that generates pseudo-computed tomography (CT) data from multi-spectral MRI acquisitions which is sufficiently similar to real CT data to enable registration of x-ray to MRI with comparable accuracy as registration of x-ray to CT. The method is based on a k-nearest-neighbors (kNN)-regression strategy which labels voxels of MRI data with CT Hounsfield Units. The regression method uses multi-spectral MRI intensities and intensity gradients as features to discriminate between various tissue types. The efficacy of using pseudo-CT data for registration of x-ray to MRI was tested on ex vivo animal data. 2D-3D registration experiments using CT and pseudo-CT data of multiple subjects were performed with a commonly used 2D-3D registration algorithm. On average, the median target registration error for registration of two x-ray images to MRI data was approximately 1 mm larger than for x-ray to CT registration. The authors have shown that pseudo-CT data generated from multi-spectral MRI facilitate registration of MRI to x-ray images. From the experiments it could be concluded that the accuracy achieved was comparable to that of registering x-ray images to CT data.
Polarization ray tracing in anisotropic optically active media
NASA Technical Reports Server (NTRS)
Mcclain, Stephen C.; Chipman, Russell A.
1992-01-01
Procedures for performing polarization ray tracing through birefringent media are presented in a form compatible with the standard methods of geometric ray tracing. The birefringent materials treated include the following: anisotropic optically active materials such as quartz, non-optically active uniaxial materials such as calcite, and isotropic optically active materials such as mercury sulfide or organic liquids. Refraction and reflection algorithms are presented which compute both ray directions and wave directions. Methods for computing polarization modes, refractive indices, optical path lengths, and Fresnel transmission and reflection coefficients are also specified.
Differential ray tracing analysis of the Schwarzschild objective
NASA Astrophysics Data System (ADS)
Prieto-Blanco, Xesús; Mouriz, Dolores; González Núñez, Héctor; Lopez Lago, Elena; de la Fuente, Raúl
2011-05-01
Differential Ray Tracing (DRT) is applied to optimize the design of a Schwarzschild objective with large aperture and for arbitrary object position. This optical system lacks of cylindrical symmetry about the non-paraxial base ray, causing astigmatism of a pencil of rays around this ray. The analysis determines the mirror radii ratio that makes the pencil anastigmatic, leading to an excellent image performance. In particular, the classical aplanatic Schwarzschild design is obtained in the limiting case where the base ray becomes paraxial. One example of a design, similar to a typical commercial objective for microscopy, is presented and the image quality is analyzed with an optical design program.
Ray-tracing code TRAVIS for ECR heating, EC current drive and ECE diagnostic
NASA Astrophysics Data System (ADS)
Marushchenko, N. B.; Turkin, Y.; Maassberg, H.
2014-01-01
A description of the recently developed ray-tracing code TRAVIS is given together with the theoretical background, results of benchmarking and examples of application. The code is written for electron cyclotron studies with emphasis on heating, current drive and ECE diagnostic. The code works with an arbitrary 3D magnetic equilibrium being applicable for both stellarators and tokamaks. The equations for ray tracing are taken in the weakly relativistic approach, i.e. with thermal effects taken into account, while the absorption, current drive and emissivity are calculated in the fully relativistic approach. For the calculation of ECCD, an adjoint technique with parallel momentum conservation is applied. The code is controlled through a specially designed graphical user interface, which allows the preparation of the input parameters and viewing the results in convenient (2D and 3D) form.
Numerical Ray-Tracing in Full Angle Spatial Compounding
NASA Astrophysics Data System (ADS)
Koch, Andreas; Koch, Ingo; Hansen, Christian; Lerch, Reinhard; Ermert, Helmut
The assumption of straight-line wave propagation is common in medical ultrasound. While sufficient for unidirectional systems, it is the main cause for degenerated FASC (Full Angle Spatial Compounding) images, where B-mode data from different viewing angles around an object, e.g. the female breast, are superimposed. To overcome this, we have implemented an eikonal equation based algorithm to perform numerical ray-tracing in inhomogeneous speed of sound distributions. Results can be used to correct ray-paths prior to FASC. Our goal was to improve FASC image quality by using numerical ray-tracing. A tissue mimicking phantom with reservoirs filled with different concentrations of saline water and correspondingly different speeds of sound was imaged with a 2.5 MHz transducer. To evaluate the isotropy of the system's spatial resolution, seven fibers were included into the phantom and the reservoirs. We compared the full width at half maximum of line scatterer images in a FASC image corrected by ray-paths from the numerical ray-tracing with an uncorrected FASC image. Results show that numerical ray-tracing improves the image contrast, eliminates double line artifacts and improves the resolution and its isotropy in FASC.
Studying the precision of ray tracing techniques with Szekeres models
NASA Astrophysics Data System (ADS)
Koksbang, S. M.; Hannestad, S.
2015-07-01
The simplest standard ray tracing scheme employing the Born and Limber approximations and neglecting lens-lens coupling is used for computing the convergence along individual rays in mock N-body data based on Szekeres swiss cheese and onion models. The results are compared with the exact convergence computed using the exact Szekeres metric combined with the Sachs formalism. A comparison is also made with an extension of the simple ray tracing scheme which includes the Doppler convergence. The exact convergence is reproduced very precisely as the sum of the gravitational and Doppler convergences along rays in Lemaitre-Tolman-Bondi swiss cheese and single void models. This is not the case when the swiss cheese models are based on nonsymmetric Szekeres models. For such models, there is a significant deviation between the exact and ray traced paths and hence also the corresponding convergences. There is also a clear deviation between the exact and ray tracing results obtained when studying both nonsymmetric and spherically symmetric Szekeres onion models.
2D/3D cryo x-ray fluorescence imaging at the bionanoprobe at the advanced photon source
NASA Astrophysics Data System (ADS)
Chen, S.; Paunesku, T.; Yuan, Y.; Deng, J.; Jin, Q.; Hong, Y. P.; Vine, D. J.; Lai, B.; Flachenecker, C.; Hornberger, B.; Brister, K.; Jacobsen, C.; Woloschak, G. E.; Vogt, S.
2016-01-01
Trace elements, particularly metals, play very important roles in biological systems. Synchrotron-based hard X-ray fluorescence microscopy offers the most suitable capabilities to quantitatively study trace metals in thick biological samples, such as whole cells and tissues. In this manuscript, we have demonstrated X-ray fluorescence imaging of frozen-hydrated whole cells using the recent developed Bionanoprobe (BNP). The BNP provides spatial resolution down to 30 nm and cryogenic capabilities. Frozen-hydrated biological cells have been directly examined on a sub-cellular level at liquid nitrogen temperatures with minimal sample preparation.
Ray-tracing-based reconstruction algorithms for digital breast tomosynthesis
NASA Astrophysics Data System (ADS)
Zhou, Weihua; Lu, Jianping; Zhou, Otto; Chen, Ying
2015-03-01
As a breast-imaging technique, digital breast tomosynthesis has great potential to improve the diagnosis of early breast cancer over mammography. Ray-tracing-based reconstruction algorithms, such as ray-tracing back projection, maximum-likelihood expectation maximization (MLEM), ordered-subset MLEM (OS-MLEM), and simultaneous algebraic reconstruction technique (SART), have been developed as reconstruction methods for different breast tomosynthesis systems. This paper provides a comparative study to investigate these algorithms by computer simulation and phantom study. Experimental results suggested that, among the four investigated reconstruction algorithms, OS-MLEM and SART performed better in interplane artifact removal with a fast speed convergence.
Modeling the effect of refraction on OCT imaging of lung tissue: a ray-tracing approach
NASA Astrophysics Data System (ADS)
Golabchi, Fatemeh N.; Golabchi, Ali; Brooks, Dana H.; Gouldstone, Andrew; DiMarzio, Charles A.
2012-03-01
Determining the structure of lung tissue is difficult in ex-vivo samples. Optical coherence tomography (OCT) can image alveoli but ignores optical effects that distort the images. For example, light refracts and changes speed at the alveolar air-tissue surface. We employ ray-tracing to model OCT imaging with directional and speed changes included, using spherical shapes in 2D. Results show apparent thickening of inter-aveolar walls and distortion of shape and depth. Our approach suggests a correction algorithm by combining the model with image analysis. Distortion correction will allow inference of tissue mechanical properties and deeper imaging.
Ray tracing a three dimensional scene using a grid
Wald, Ingo; Ize, Santiago; Parker, Steven G; Knoll, Aaron
2013-02-26
Ray tracing a three-dimensional scene using a grid. One example embodiment is a method for ray tracing a three-dimensional scene using a grid. In this example method, the three-dimensional scene is made up of objects that are spatially partitioned into a plurality of cells that make up the grid. The method includes a first act of computing a bounding frustum of a packet of rays, and a second act of traversing the grid slice by slice along a major traversal axis. Each slice traversal includes a first act of determining one or more cells in the slice that are overlapped by the frustum and a second act of testing the rays in the packet for intersection with any objects at least partially bounded by the one or more cells overlapped by the frustum.
NASA Astrophysics Data System (ADS)
Merkel, Ronny; Gruhn, Stefan; Dittmann, Jana; Vielhauer, Claus; Bräutigam, Anja
2012-03-01
Determining the age of latent fingerprint traces found at crime scenes is an unresolved research issue since decades. Solving this issue could provide criminal investigators with the specific time a fingerprint trace was left on a surface, and therefore would enable them to link potential suspects to the time a crime took place as well as to reconstruct the sequence of events or eliminate irrelevant fingerprints to ensure privacy constraints. Transferring imaging techniques from different application areas, such as 3D image acquisition, surface measurement and chemical analysis to the domain of lifting latent biometric fingerprint traces is an upcoming trend in forensics. Such non-destructive sensor devices might help to solve the challenge of determining the age of a latent fingerprint trace, since it provides the opportunity to create time series and process them using pattern recognition techniques and statistical methods on digitized 2D, 3D and chemical data, rather than classical, contact-based capturing techniques, which alter the fingerprint trace and therefore make continuous scans impossible. In prior work, we have suggested to use a feature called binary pixel, which is a novel approach in the working field of fingerprint age determination. The feature uses a Chromatic White Light (CWL) image sensor to continuously scan a fingerprint trace over time and retrieves a characteristic logarithmic aging tendency for 2D-intensity as well as 3D-topographic images from the sensor. In this paper, we propose to combine such two characteristic aging features with other 2D and 3D features from the domains of surface measurement, microscopy, photography and spectroscopy, to achieve an increase in accuracy and reliability of a potential future age determination scheme. Discussing the feasibility of such variety of sensor devices and possible aging features, we propose a general fusion approach, which might combine promising features to a joint age determination scheme
NASA Astrophysics Data System (ADS)
Powell, Leila C.; Kay, Scott T.; Babul, Arif
2009-12-01
Recent X-ray and weak-lensing observations of galaxy clusters have revealed that the hot gas does not always directly trace the dark matter within these systems. Such configurations are extremely interesting. They offer a new vista on to the complex interplay between gravity and baryonic physics, and may even be used as indicators of the clusters' dynamical state. In this paper, we undertake a study to determine what insight can be reliably gleaned from the comparison of the X-ray and the weak-lensing mass maps of galaxy clusters. We do this by investigating the two-dimensional (2D) substructure within three high-resolution cosmological simulations of galaxy clusters. Our main results focus on non-radiative gas dynamics, but we also consider the effects of radiative cooling at high redshift. For our analysis, we use a novel approach, based on unsharp-masking, to identify substructures in 2D surface mass density and X-ray surface brightness maps. At full resolution (~15h-1 kpc), this technique is capable of identifying almost all self-bound dark matter subhaloes with M > 1012h-1Msolar. We also report a correlation between the mass of a subhalo and the area of its corresponding 2D detection; such a correlation, once calibrated, could provide a useful estimator for substructure mass. Comparing our 2D mass and X-ray substructures, we find a surprising number of cases where the matching fails: around one-third of galaxy-sized substructures have no X-ray counterpart. Some interesting cases are also found at larger masses, in particular the cores of merging clusters where the situation can be complex. Finally, we degrade our mass maps to what is currently achievable with weak-lensing observations (~100h-1kpc at z = 0.2). While the completeness mass limit increases by around an order of magnitude, a mass-area correlation remains. Our paper clearly demonstrates that the next generation of lensing surveys should start to reveal a wealth of information on cluster substructure.
Ray tracing in discontinuous velocity model with implicit Interface
NASA Astrophysics Data System (ADS)
Zhang, Jianxing; Yang, Qin; Meng, Xianhai; Li, Jigang
2016-07-01
Ray tracing in the velocity model containing complex discontinuities is still facing many challenges. The main difficulty arises from the detection of the spatial relationship between the rays and the interfaces that are usually described in non-linear parametric forms. We propose a novel model representation method that can facilitate the implementation of classical shooting-ray methods. In the representation scheme, each interface is expressed as the zero contour of a signed distance field. A multi-copy strategy is adopted to describe the volumetric properties within blocks. The implicit description of the interface makes it easier to detect the ray-interface intersection. The direct calculation of the intersection point is converted into the problem of judging the signs of a ray segment's endpoints. More importantly, the normal to the interface at the intersection point can be easily acquired according to the signed distance field of the interface. The multiple storage of the velocity property in the proximity of the interface can provide accurate and unambiguous velocity information of the intersection point. Thus, the departing ray path can be determined easily and robustly. In addition, the new representation method can describe velocity models containing very complex geological structures, such as faults, salt domes, intrusions, and pinches, without any simplification. The examples on synthetic and real models validate the robustness and accuracy of the ray tracing based on the proposed model representation scheme.
Küchemann, Stefan; Mahn, Carsten; Samwer, Konrad
2014-01-15
The investigation of short time dynamics using X-ray scattering techniques is commonly limited either by the read out frequency of the detector or by a low intensity. In this paper, we present a chopper system, which can increase the temporal resolution of 2D X-ray detectors by a factor of 13. This technique only applies to amorphous or polycrystalline samples due to their circular diffraction patterns. Using the chopper, we successfully increased the temporal resolution up to 5.1 ms during synchrotron experiments. For the construction, we provide a mathematical formalism, which, in principle, allows an even higher increase of the temporal resolution.
Ray Tracing Study of Magnetospheric ULF Wave Propagation.
NASA Astrophysics Data System (ADS)
Zhang, Xinbo
1993-01-01
A semi-empirical plasma density model and Mead -Fairfield magnetic field model are incorporated into a 3-D ray tracing code to study magnetospheric ULF wave propagation from the subsolar magnetopause. The ray-tracing of Pc3 compressional waves from the magnetosheath reveals that the magnetosphere can present a major propagation barrier to the penetration of these waves to the plasmasphere. This barrier is the ion-ion cutoff between the He^+ and O ^+ gyroresonances. As a result of the frequency -dependent location of this cutoff, the magnetosphere behaves like a filter for Pc3 compressional waves, and only the low frequency components can penetrate to the inner magnetosphere. These results are in agreement with previous satellite observations. This 'filter action' strongly depends on the relative concentration of He^+ and O^+ and is, therefore, sensitive to solar and magnetic activity. The study of the propagation characteristics of Pc3 transverse Alfven waves shows that these waves cannot penetrate to low Earth altitudes for wave frequencies above about approximately 0.03 hz. The configuration of the refractive index reveals an O^+-He^+ associated cutoff located between the assumed wave source in the equatorial magnetopause and the Earth. When the O^+ concentration is removed from the plasma composition, the barrier no longer exists, and waves with much higher frequencies than 0.03 Hz can penetrate to low altitudes. The result that the 0.03 Hz or lower frequency Alfven waves can be guided to the low altitudes agrees with ground-based power spectrum observations at high latitudes. The ray tracing study of Pc 1-2 waves reproduces earlier results (Rauch and Roux, 1982) for an H ^+-He^+ two-ion-species plasma, i.e. Pc 1-2 left hand polarized Alfven mode waves originating at equatorial geostationary orbit, below He ^+ gyrofrequency, are guided to the ground. However, our ray tracing study shows that previous Pc 1-2 ray tracing results are only valid in the absence of O
Ray tracing reconstruction investigation for C-arm tomosynthesis
NASA Astrophysics Data System (ADS)
Malalla, Nuhad A. Y.; Chen, Ying
2016-04-01
C-arm tomosynthesis is a three dimensional imaging technique. Both x-ray source and the detector are mounted on a C-arm wheeled structure to provide wide variety of movement around the object. In this paper, C-arm tomosynthesis was introduced to provide three dimensional information over a limited view angle (less than 180o) to reduce radiation exposure and examination time. Reconstruction algorithms based on ray tracing method such as ray tracing back projection (BP), simultaneous algebraic reconstruction technique (SART) and maximum likelihood expectation maximization (MLEM) were developed for C-arm tomosynthesis. C-arm tomosynthesis projection images of simulated spherical object were simulated with a virtual geometric configuration with a total view angle of 40 degrees. This study demonstrated the sharpness of in-plane reconstructed structure and effectiveness of removing out-of-plane blur for each reconstruction algorithms. Results showed the ability of ray tracing based reconstruction algorithms to provide three dimensional information with limited angle C-arm tomosynthesis.
Removal of t1 noise from metabolomic 2D 1H- 13C HSQC NMR spectra by Correlated Trace Denoising
NASA Astrophysics Data System (ADS)
Poulding, Simon; Charlton, Adrian J.; Donarski, James; Wilson, Julie C.
2007-12-01
The presence of t1 noise artefacts in 2D phase-cycled Heteronuclear Single Quantum Coherence (HSQC) spectra constrains the use of this experiment despite its superior sensitivity. This paper proposes a new processing algorithm, working in the frequency-domain, for reducing t1 noise. The algorithm has been developed for use in contexts, such as metabolomic studies, where existing denoising techniques cannot always be applied. Two test cases are presented that show the algorithm to be effective in improving the SNR of peaks embedded within t1 noise by a factor of more than 2, while retaining the intensity and shape of genuine peaks.
Rigid 2D/3D registration of intraoperative digital x-ray images and preoperative CT and MR images
NASA Astrophysics Data System (ADS)
Tomazevic, Dejan; Likar, Bostjan; Pernus, Franjo
2002-05-01
This paper describes a novel approach to register 3D computed tomography (CT) or magnetic resonance (MR) images to a set of 2D X-ray images. Such a registration may be a valuable tool for intraoperative determination of the precise position and orientation of some anatomy of interest, defined in preoperative images. The registration is based solely on the information present in 2D and 3D images. It does not require fiducial markers, X-ray image segmentation, or construction of digitally reconstructed radiographs. The originality of the approach is in using normals to bone surfaces, preoperatively defined in 3D MR or CT data, and gradients of intraoperative X-ray images, which are back-projected towards the X-ray source. The registration is then concerned with finding that rigid transformation of a CT or MR volume, which provides the best match between surface normals and back projected gradients, considering their amplitudes and orientations. The method is tested on a lumbar spine phantom. Gold standard registration is obtained by fidicual markers attached to the phantom. Volumes of interest, containing single vertebrae, are registered to different pairs of X-ray images from different starting positions, chosen randomly and uniformly around the gold standard position. Target registration errors and rotation errors are in order of 0.3 mm and 0.35 degrees for the CT to X-ray registration and 1.3 mm and 1.5 degrees for MR to X-ray registration. The registration is shown to be fast and accurate.
Auroral kilometric radiation source characteristics using ray tracing techniques
NASA Astrophysics Data System (ADS)
Schreiber, R.; Santolik, O.; Parrot, M.; Lefeuvre, F.; Hanasz, J.; Brittnacher, M.; Parks, G.
2002-11-01
3-D ray tracing to the presumed auroral kilometric radiation (AKR) source region has been performed using the input data from wave distribution function (WDF) based on the AKR waveforms recorded on board the Interball 2 satellite by the French wave experiment MEMO. Both the direction of the WDF maximum and the WDF form and angular size have been taken into account. Two instances of AKR emissions were observed on 28 January 1997 at 2037 and 2107 UT. Rays traced in R-X mode out of the s/c point toward two different active regions on the auroral oval (as seen with Polar UV imager after projection of the source region along the magnetic field lines down to the ionosphere level). Source region apparent angular sizes based on WDF are compatible with sizes estimated from signal modulation produced by electric antenna system rotation.
Fast stereoscopic images with ray-traced volume rendering
Adelson, S.J.; Hansen, C.D.
1994-05-01
One of the drawbacks of standard volume rendering techniques is that is it often difficult to comprehend the three-dimensional structure of the volume from a single frame; this is especially true in cases where there is no solid surface. Generally, several frames must be generated and viewed sequentially, using motion parallax to relay depth. Another option is to generate a single spectroscopic pair, resulting in clear and unambiguous depth information in both static and moving images. Methods have been developed which take advantage of the coherence between the two halves of a stereo pair for polygon rendering and ray-tracing, generating the second half of the pair in significantly less time than that required to completely render a single image. This paper reports the results of implementing these techniques with parallel ray-traced volume rendering. In tests with different data types, the time savings is in the range of 70--80%.
Ray tracing homogenizing mirrors for synchrotron x-ray lithography
NASA Astrophysics Data System (ADS)
Homer, Michael; Rosser, Roy J.; Speer, R. J.
1991-12-01
Saddle toroid array mirrors (STAMs) are novel grazing-incidence mirrors. They have been proposed as the optical component that most efficiently matches synchrotron orbital radiation (SOR) to the needs of proximity x-ray lithography. However, STAMs have yet to be accepted by the synchrotron lithography community because of the lack of detailed data on their expected performance, due primarily to the difficulty of raytracing such mirrors using existing optical raytrace programs. A raytracing package written especially to study the design and optimization of these unusually shaped mirrors and the very encouraging results obtained with the package to date are described. The optimum STAM designs turn out to be the most effective way of homogeneously illuminating a rectangular proximity x-ray lithography mask, improving on existing scanning mirror systems by at least a factor of four. They have the added advantage of being stationary, which should lead to greater reliability--a quality of considerable value in the production environment these mirrors are intended for, namely the ultra-high vacuum of a synchrotron beamline. Based on the results of the raytracing, a prototype STAM has been constructed, and preparations are being made for an x-ray test of the device.
2D and 3D Refraction Based X-ray Imaging Suitable for Clinical and Pathological Diagnosis
Ando, Masami; Bando, Hiroko; Ueno, Ei
2007-01-19
The first observation of micro papillary (MP) breast cancer by x-ray dark-field imaging (XDFI) and the first observation of the 3D x-ray internal structure of another breast cancer, ductal carcinoma in-situ (DCIS), are reported. The specimen size for the sheet-shaped MP was 26 mm x 22 mm x 2.8 mm, and that for the rod-shaped DCIS was 3.6 mm in diameter and 4.7 mm in height. The experiment was performed at the Photon Factory, KEK: High Energy Accelerator Research Organization. We achieved a high-contrast x-ray image by adopting a thickness-controlled transmission-type angular analyzer that allows only refraction components from the object for 2D imaging. This provides a high-contrast image of cancer-cell nests, cancer cells and stroma. For x-ray 3D imaging, a new algorithm due to the refraction for x-ray CT was created. The angular information was acquired by x-ray optics diffraction-enhanced imaging (DEI). The number of data was 900 for each reconstruction. A reconstructed CT image may include ductus lactiferi, micro calcification and the breast gland. This modality has the possibility to open up a new clinical and pathological diagnosis using x-ray, offering more precise inspection and detection of early signs of breast cancer.
2D and 3D Refraction Based X-ray Imaging Suitable for Clinical and Pathological Diagnosis
NASA Astrophysics Data System (ADS)
Ando, Masami; Bando, Hiroko; Chen, Zhihua; Chikaura, Yoshinori; Choi, Chang-Hyuk; Endo, Tokiko; Esumi, Hiroyasu; Gang, Li; Hashimoto, Eiko; Hirano, Keiichi; Hyodo, Kazuyuki; Ichihara, Shu; Jheon, SangHoon; Kim, HongTae; Kim, JongKi; Kimura, Tatsuro; Lee, ChangHyun; Maksimenko, Anton; Ohbayashi, Chiho; Park, SungHwan; Shimao, Daisuke; Sugiyama, Hiroshi; Tang, Jintian; Ueno, Ei; Yamasaki, Katsuhito; Yuasa, Tetsuya
2007-01-01
The first observation of micro papillary (MP) breast cancer by x-ray dark-field imaging (XDFI) and the first observation of the 3D x-ray internal structure of another breast cancer, ductal carcinoma in-situ (DCIS), are reported. The specimen size for the sheet-shaped MP was 26 mm × 22 mm × 2.8 mm, and that for the rod-shaped DCIS was 3.6 mm in diameter and 4.7 mm in height. The experiment was performed at the Photon Factory, KEK: High Energy Accelerator Research Organization. We achieved a high-contrast x-ray image by adopting a thickness-controlled transmission-type angular analyzer that allows only refraction components from the object for 2D imaging. This provides a high-contrast image of cancer-cell nests, cancer cells and stroma. For x-ray 3D imaging, a new algorithm due to the refraction for x-ray CT was created. The angular information was acquired by x-ray optics diffraction-enhanced imaging (DEI). The number of data was 900 for each reconstruction. A reconstructed CT image may include ductus lactiferi, micro calcification and the breast gland. This modality has the possibility to open up a new clinical and pathological diagnosis using x-ray, offering more precise inspection and detection of early signs of breast cancer.
3D ultrasonic ray tracing in AutoCAD®
NASA Astrophysics Data System (ADS)
Reilly, D.; Leggat, P.; McNab, A.
2001-04-01
To assist with the design and validation of testing procedures for NDT, add-on modules have been developed for AutoCAD® 2000. One of the modules computes and displays ultrasonic 3D ray tracing. Another determines paths between two points, for instance a probe and a target or two probes. The third module displays phased array operational modes and calculates element delays for phased array operation. The modules can be applied to simple or complex solid model components.
grtrans: Polarized general relativistic radiative transfer via ray tracing
NASA Astrophysics Data System (ADS)
Dexter, Jason
2016-05-01
grtrans calculates ray tracing radiative transfer in the Kerr metric, including the full treatment of polarised radiative transfer and parallel transport along geodesics, for comparing theoretical models of black hole accretion flows and jets with observations. The code is written in Fortran 90 and parallelizes with OpenMP; the full code and several components have Python interfaces. grtrans includes Geokerr (ascl:1011.015) and requires cfitsio (ascl:1010.001) and pyfits (ascl:1207.009).
Ray tracing for point distribution in unstructured grid generation
Khamayseh, A.; Ortega, F.; Trease, H.
1995-12-31
We present a procedure by which grid points are generated on surfaces or within three-dimensional volumes to produce high quality unstructed grids for complex geometries. The virtue of this method is based on ray-tracing approach for curved polyhedra whose faces may lie on natural quadrics (planes, cylinders, cones, or spheres) or triangular faceted surfaces. We also present an efficient point location algorithm for identifying points relative to various regions with classification of inside/on/outside.
grtrans: Polarized general relativistic radiative transfer via ray tracing
NASA Astrophysics Data System (ADS)
Dexter, Jason
2016-05-01
grtrans calculates ray tracing radiative transfer in the Kerr metric, including the full treatment of polarised radiative transfer and parallel transport along geodesics, for comparing theoretical models of black hole accretion flows and jets with observations. The code is written in Fortran 90 and parallelizes with OpenMP; the full code and several components have Python interfaces. grtrans requires Geokerr (ascl:1011.015), cfitsio (ascl:1010.001), and pyfits (ascl:1207.009).
Ray tracing software application in VIP lamp design
NASA Astrophysics Data System (ADS)
Rehn, Henning
2002-08-01
In our contribution we demonstrate a wide variety of ray tracing software applications for the design of VIP short-arc discharge video projection lamps. On the basis of simulations we derive design rules for the lamp itself and for its optical environment. Light Tools software acts as a means to understand the collection efficiency of a VIP lamp with an elliptical reflector and as an instrument to prove the conclusions.
Accelerated ray tracing algorithm under urban macro cell
NASA Astrophysics Data System (ADS)
Liu, Z.-Y.; Guo, L.-X.; Guan, X.-W.
2015-10-01
In this study, an ray tracing propagation prediction model, which is based on creating a virtual source tree, is used because of their high efficiency and reliable prediction accuracy. In addition, several acceleration techniques are also adopted to improve the efficiency of ray-tracing-based prediction over large areas. However, in the process of employing the ray tracing method for coverage zone prediction, runtime is linearly proportional to the total number of prediction points, leading to large and sometimes prohibitive computation time requirements under complex geographical urban macrocell environments. In order to overcome this bottleneck, the compute unified device architecture (CUDA), which provides fine-grained data parallelism and thread parallelism, is implemented to accelerate the calculation. Taking full advantage of tens of thousands of threads in CUDA program, the decomposition of the coverage prediction problem is firstly conducted by partitioning the image tree and the visible prediction points to different sources. Then, we make every thread calculate the electromagnetic field of one propagation path and then collect these results. Comparing this parallel algorithm with the traditional sequential algorithm, it can be found that computational efficiency has been improved.
Three dimensional ray tracing of the Jovian magnetosphere in the low frequency range
NASA Technical Reports Server (NTRS)
Menietti, J. D.
1984-01-01
Ray tracing studies of Jovian low frequency emissions were studied. A comprehensive three-dimensional ray tracing computer code for examination of model Jovian decametric (DAM) emission was developed. The improvements to the computer code are outlined and described. The results of the ray tracings of Jovian emissions will be presented in summary form.
NASA Astrophysics Data System (ADS)
Dinten, Jean-Marc; Robert-Coutant, Christine; Darboux, Michel
2001-06-01
Dual Energy X-Rays Absorptiometry (DXA) is commonly used to separate soft tissues and bone contributions in radiographs. This decomposition leads to bone mineral density (BMD) measurement. Most clinical systems use pencil or fan collimated X-Rays beam with mono detectors or linear arrays. On these systems BMD is computed from bi-dimensional (2D) images obtained by scanning. Our objective is to take advantage of the newly available flat panels detectors and to propose a DXA approach without scanning, based on the use of cone beam X-Rays associated with a 2D detector. This approach yields bone densitometry systems with an equal X and Y resolution, a fast acquisition and a reduced risk of patient motion.Scatter in this case becomes an important issue. While scattering is insignificant on collimated systems, its level and geometrical structure may severely alter BMD measurement on cone beam systems. In our presentation an original DXA method taking into account scattering is proposed. This new approach leads to accurate BMD values.In order to evaluate the accuracy of our new approach, a phantom representative of the spine regions tissue composition (bone, fat , muscle) has been designed. The comparison between the expected theoretical and the reconstructed BMD values validates the accuracy of our method. Results on anthropomorphic spine and hip regions are also presented.
Twin robotic x-ray system for 2D radiographic and 3D cone-beam CT imaging
NASA Astrophysics Data System (ADS)
Fieselmann, Andreas; Steinbrener, Jan; Jerebko, Anna K.; Voigt, Johannes M.; Scholz, Rosemarie; Ritschl, Ludwig; Mertelmeier, Thomas
2016-03-01
In this work, we provide an initial characterization of a novel twin robotic X-ray system. This system is equipped with two motor-driven telescopic arms carrying X-ray tube and flat-panel detector, respectively. 2D radiographs and fluoroscopic image sequences can be obtained from different viewing angles. Projection data for 3D cone-beam CT reconstruction can be acquired during simultaneous movement of the arms along dedicated scanning trajectories. We provide an initial evaluation of the 3D image quality based on phantom scans and clinical images. Furthermore, initial evaluation of patient dose is conducted. The results show that the system delivers high image quality for a range of medical applications. In particular, high spatial resolution enables adequate visualization of bone structures. This system allows 3D X-ray scanning of patients in standing and weight-bearing position. It could enable new 2D/3D imaging workflows in musculoskeletal imaging and improve diagnosis of musculoskeletal disorders.
Ray Tracing for Complex Astrophysical High-opacity Structures
NASA Astrophysics Data System (ADS)
Steinacker, J.; Bacmann, A.; Henning, T.
2006-07-01
We present a ray-tracing technique for radiative transfer modeling of complex three-dimensional (3D) structures that include dense regions of high optical depth, such as that in dense molecular clouds, circumstellar disks, envelopes of evolved stars, and dust tori around active galactic nuclei. The corresponding continuum radiative transfer problem is described, and the numerical requirements for inverse 3D density and temperature modeling are defined. We introduce a relative intensity and transform the radiative transfer equation along the rays to solve machine precision problems and to relax strong gradients in the source term. For the optically thick regions where common ray tracers are forced to perform small trace steps, we give two criteria for making use of a simple approximative solver crossing the optically thick region quickly. Using an example of a density structure with optical depth changes of 6 orders of magnitude and sharp temperature variations, we demonstrate the accuracy of the proposed scheme using a common fifth-order Runge-Kutta ray tracer with adaptive step-size control. In our test case, the gain in computational speed is about a factor of 870. The method is applied in order to calculate the temperature distribution within a massive molecular cloud core for different boundary conditions for the radiation field.
2D grating simulation for X-ray phase-contrast and dark-field imaging with a Talbot interferometer
NASA Astrophysics Data System (ADS)
Zanette, Irene; David, Christian; Rutishauser, Simon; Weitkamp, Timm
2010-04-01
Talbot interferometry is a recently developed and an extremely powerful X-ray phase-contrast imaging technique. Besides giving access to ultra-high sensitivity differential phase contrast images, it also provides the dark field image, which is a map of the scattering power of the sample. In this paper we investigate the potentialities of an improved version of the interferometer, in which two dimensional gratings are used instead of standard line grids. This approach allows to overcome the difficulties that might be encountered in the images produced by a one dimensional interferometer. Among these limitations there are the phase wrapping and quantitative phase retrieval problems and the directionality of the differential phase and dark-field signals. The feasibility of the 2D Talbot interferometer has been studied with a numerical simulation on the performances of its optical components under different circumstances. The gratings can be obtained either by an ad hoc fabrication of the 2D structures or by a superposition of two perpendicular linear grids. Through this simulation it has been possible to find the best parameters for a practical implementation of the 2D Talbot interferometer.
Implementation of a new multiple monochromatic x-ray 2D imager at NIF
NASA Astrophysics Data System (ADS)
Kyrala, G. A.; Martinson, D.; Polk, P. J.; Gravlin, T.; Schmitt, M. J.; Johnson, R.; Murphy, T. J.; Lopez, F. E.; Oertel, J. A.; House, A.; Wood, R.; Lee, J.; Haugh, M.
2013-09-01
We will describe the installation and wavelength calibration of a multiple monochromatic imager [MMI]1 to be used on mix experiments at National Ignition Facility [NIF]2. The imager works between 8 and 13 keV, has a spatial resolution of 16 micrometers and generates many images each with an energy bandwidth of ~80 eV. The images are recorded either on image plates or on gated x-ray detectors. We will describe: how we aligned the instrument on the bench using visible light, how we checked the alignment and determined the energy range using a k-alpha x-ray source, and how we installed and aligned the instrument to the NIF target chamber.
Coronary arteries motion modeling on 2D x-ray images
NASA Astrophysics Data System (ADS)
Gao, Yang; Sundar, Hari
2012-02-01
During interventional procedures, 3D imaging modalities like CT and MRI are not commonly used due to interference with the surgery and radiation exposure concerns. Therefore, real-time information is usually limited and building models of cardiac motion are difficult. In such case, vessel motion modeling based on 2-D angiography images become indispensable. Due to issues with existing vessel segmentation algorithms and the lack of contrast in occluded vessels, manual segmentation of certain branches is usually necessary. In addition, such occluded branches are the most important vessels during coronary interventions and obtaining motion models for these can greatly help in reducing the procedure time and radiation exposure. Segmenting different cardiac phases independently does not guarantee temporal consistency and is not efficient for occluded branches required manual segmentation. In this paper, we propose a coronary motion modeling system which extracts the coronary tree for every cardiac phase, maintaining the segmentation by tracking the coronary tree during the cardiac cycle. It is able to map every frame to the specific cardiac phase, thereby inferring the shape information of the coronary arteries using the model corresponding to its phase. Our experiments show that our motion modeling system can achieve promising results with real-time performance.
Balosso, Jacques
2016-01-01
Background The advanced dose calculation algorithms implemented in treatment planning system (TPS) have remarkably improved the accuracy of dose calculation especially the modeling of electrons transport in the low density medium. The purpose of this study is to evaluate the use of 2D gamma (γ) index to quantify and evaluate the impact of the calculation of electrons transport on dose distribution for lung radiotherapy. Methods X-ray computed tomography images were used to calculate the dose for twelve radiotherapy treatment plans. The doses were originally calculated with Modified Batho (MB) 1D density correction method, and recalculated with anisotropic analytical algorithm (AAA), using the same prescribed dose. Dose parameters derived from dose volume histograms (DVH) and target coverage indices were compared. To compare dose distribution, 2D γ-index was applied, ranging from 1%/1 mm to 6%/6 mm. The results were displayed using γ-maps in 2D. Correlation between DVH metrics and γ passing rates was tested using Spearman’s rank test and Wilcoxon paired test to calculate P values. Results the plans generated with AAA predicted more heterogeneous dose distribution inside the target, with P<0.05. However, MB overestimated the dose predicting more coverage of the target by the prescribed dose. The γ analysis showed that the difference between MB and AAA could reach up to ±10%. The 2D γ-maps illustrated that AAA predicted more dose to organs at risks, as well as lower dose to the target compared to MB. Conclusions Taking into account of the electrons transport on radiotherapy plans showed a significant impact on delivered dose and dose distribution. When considering the AAA represent the true cumulative dose, a readjusting of the prescribed dose and an optimization to protect the organs at risks should be taken in consideration in order to obtain the better clinical outcome. PMID:27429908
Graphical User Interface for Interactive Seismic Ray Tracing
NASA Astrophysics Data System (ADS)
Song, Jianli; ten Brink, Uri
2005-03-01
RayGUI 2.0 is a new version of RayGUI, a graphical user interface (GUI) to the seismic travel time modeling program of Zelt and Smith [1992]. It represents a significant improvement over the previous version of RayGUI (RayGUI 1.04; Loss et al. [1998a,1998b]). RayGUI 2.0 uses an updated Java version (1.3), and can run on various operating systems (UNIX, Linux, and Mac OS X). Several new functions have been incorporated, including executing the forward and inversion codes of Zelt and Smith [1992], creating models or adding new parts of models from an ASCII file, graphically adding layers or points, graphically pinching layers, changing the velocity value of a control point, reporting point location and velocity, importing travel-time lists, generating postscript files, exporting the velocity model into an ASCII file, generating 1-D velocity profiles at specified locations, calculating root-mean-square errors between observed and calculated arrivals for selected phases, and accessing the ray trace log, as well as several other new display features.
NASA Astrophysics Data System (ADS)
Hatt, Charles R.; Speidel, Michael A.; Raval, Amish N.
2014-03-01
We present a novel 2D/ 3D registration algorithm for fusion between transesophageal echocardiography (TEE) and X-ray fluoroscopy (XRF). The TEE probe is modeled as a subset of 3D gradient and intensity point features, which facilitates efficient 3D-to-2D perspective projection. A novel cost-function, based on a combination of intensity and edge features, evaluates the registration cost value without the need for time-consuming generation of digitally reconstructed radiographs (DRRs). Validation experiments were performed with simulations and phantom data. For simulations, in silica XRF images of a TEE probe were generated in a number of different pose configurations using a previously acquired CT image. Random misregistrations were applied and our method was used to recover the TEE probe pose and compare the result to the ground truth. Phantom experiments were performed by attaching fiducial markers externally to a TEE probe, imaging the probe with an interventional cardiac angiographic x-ray system, and comparing the pose estimated from the external markers to that estimated from the TEE probe using our algorithm. Simulations found a 3D target registration error of 1.08(1.92) mm for biplane (monoplane) geometries, while the phantom experiment found a 2D target registration error of 0.69mm. For phantom experiments, we demonstrated a monoplane tracking frame-rate of 1.38 fps. The proposed feature-based registration method is computationally efficient, resulting in near real-time, accurate image based registration between TEE and XRF.
Time Resolved 2D X-Ray Densitometry of a Ventilated Partial Cavity Closure
NASA Astrophysics Data System (ADS)
Makiharju, Simo; Ceccio, Steven
2011-11-01
A time resolved x-ray densitometry system was developed to measure the spatial distribution of void fraction for nominally two-dimensional flows. The system can image a region of (15 cm)2 at a frame rate of up to 4000 fps. The source was a rotating anode type normally used for cineradiography and angiography. Supplied by a 65 kW high frequency generator with a high speed starter, it could be operated at up to 433 mA at 150 kV. The imager subsystem comprised of a high speed camera coupled with a high resolution image intensifier. The range of measured void fraction can be changed to span a desired range yielding an uncertainty on the order of 1% of the measurement range. The system is used to examine the void fraction field in the closure region of a ventilated partial cavity behind a backward facing step. The cavity has Reynolds number of O(105) based on the cavity length, and a non-dimensional gas flux of Q* = 0.0048. The bubbly flow created in the cavity wake is examined using the x-ray densitometry system, duel fiber optical probes, and high speed cinematography. The local void fraction and bubble size distributions in the cavity wake are determined, and the measurements methods are compared. The research was sponsored by ONR under grant N00014-08-1-0215, program manager Dr. L. Patrick Purtell.
2D-Omnidirectional Hard-X-Ray Scattering Sensitivity in a Single Shot.
Kagias, Matias; Wang, Zhentian; Villanueva-Perez, Pablo; Jefimovs, Konstantins; Stampanoni, Marco
2016-03-01
X-ray scattering imaging can provide complementary information to conventional absorption based radiographic imaging about the unresolved microstructures of a sample. The scattering signal can be accessed with various methods based on coherent illumination, which span from self-imaging to speckle scanning. The directional sensitivity of the existing real space imaging methods is limited to a few directions on the imaging plane and requires scanning of the optical components, or the rotation of either the sample or the imaging setup, in order to cover the full range of possible scattering directions. In this Letter the authors propose a new method that allows the simultaneous acquisition of scattering images in all possible directions in a single shot. This is achieved by a specialized phase grating and a detector with sufficient spatial resolution to record the generated interference fringe. The structural length scale sensitivity of the system can be tuned by varying its geometry for a fixed grating design. Taking into account ongoing developments in the field of compact x-ray sources that allow high brightness and sufficient spatial coherence, the applicability of omnidirectional scattering imaging in industrial and medical settings is boosted significantly. PMID:26991177
2D x-ray imaging spectroscopic diagnostics using convex bent crystal
NASA Astrophysics Data System (ADS)
Papp, Daniel; Presura, Radu; Wallace, Matt; Largent, Billy; Haque, Showera; Arias, Angel; Khanal, Vijay; Ivanov, Vladimir
2013-10-01
A new 2-dimensional time-integrated x-ray spectroscopic diagnostics technique was developed to create multi-monochromatic images of high-energy density Al plasmas. 2-dimensional is an advanced spectroscopic tool, providing a way to determine the spatial dependence of plasma temperature and density (Te and ne) in hot plasmas. The new technique uses the strong source broadening of convex cylindrically bent KAP crystal spectrometers, which contains spatial information along the dispersive direction. The perpendicular direction is imaged using a slit. The spatial resolution of the method is improved by the deconvolution of the source broadened line profiles from the lineshapes (recorded by the convex crystal spectrometer) with lineshapes of minimum instrumental broadening. The latter spectra were recorded with a concave cylindrically bent KAP crystal spectrometer, based on the Johann geometry. Spectroscopic model of the plasma x-ray emission was developed using the PrismSPECT code. The identification of suitable spectral features allows deriving Te and ne from line intensities. We applied this model to get temperature and density distribution maps for wire array z-pinch plasmas. Work supported by the DOE/NNSA under grant DE-NA0001834 and Cooperative Agreement DE-FC52-06NA27616.
2D-Omnidirectional Hard-X-Ray Scattering Sensitivity in a Single Shot
NASA Astrophysics Data System (ADS)
Kagias, Matias; Wang, Zhentian; Villanueva-Perez, Pablo; Jefimovs, Konstantins; Stampanoni, Marco
2016-03-01
X-ray scattering imaging can provide complementary information to conventional absorption based radiographic imaging about the unresolved microstructures of a sample. The scattering signal can be accessed with various methods based on coherent illumination, which span from self-imaging to speckle scanning. The directional sensitivity of the existing real space imaging methods is limited to a few directions on the imaging plane and requires scanning of the optical components, or the rotation of either the sample or the imaging setup, in order to cover the full range of possible scattering directions. In this Letter the authors propose a new method that allows the simultaneous acquisition of scattering images in all possible directions in a single shot. This is achieved by a specialized phase grating and a detector with sufficient spatial resolution to record the generated interference fringe. The structural length scale sensitivity of the system can be tuned by varying its geometry for a fixed grating design. Taking into account ongoing developments in the field of compact x-ray sources that allow high brightness and sufficient spatial coherence, the applicability of omnidirectional scattering imaging in industrial and medical settings is boosted significantly.
Dynamic ray tracing and its application in triangulated media
Rueger, A.
1993-07-01
Hale and Cohen (1991) developed software to generate two-dimensional computer models of complex geology. Their method uses a triangulation technique designed to support efficient and accurate computation of seismic wavefields for models of the earth`s interior. Subsequently, Hale (1991) used this triangulation approach to perform dynamic ray tracing and create synthetic seismograms based on the method of Gaussian beams. Here, I extend this methodology to allow an increased variety of ray-theoretical experiments. Specifically, the developed program GBmod (Gaussian Beam MODeling) can produce arbitrary multiple sequences and incorporate attenuation and density variations. In addition, I have added an option to perform Fresnel-volume ray tracing (Cerveny and Soares, 1992). Corrections for reflection and transmission losses at interfaces, and for two-and-one-half-dimensional (2.5-D) spreading are included. However, despite these enhancements, difficulties remain in attempts to compute accurate synthetic seismograms if strong lateral velocity inhomogeneities are present. Here, these problems are discussed and, to a certain extent, reduced. I provide example computations of high-frequency seismograms based on the method of Gaussian beams to exhibit the advantages and disadvantages of the proposed modeling method and illustrate new features for both surface and vertical seismic profiling (VSP) acquisition geometries.
Adaptive sample map for Monte Carlo ray tracing
NASA Astrophysics Data System (ADS)
Teng, Jun; Luo, Lixin; Chen, Zhibo
2010-07-01
Monte Carlo ray tracing algorithm is widely used by production quality renderers to generate synthesized images in films and TV programs. Noise artifact exists in synthetic images generated by Monte Carlo ray tracing methods. In this paper, a novel noise artifact detection and noise level representation method is proposed. We first apply discrete wavelet transform (DWT) on a synthetic image; the high frequency sub-bands of the DWT result encode the noise information. The sub-bands coefficients are then combined to generate a noise level description of the synthetic image, which is called noise map in the paper. This noise map is then subdivided into blocks for robust noise level metric calculation. Increasing the samples per pixel in Monte Carlo ray tracer can reduce the noise of a synthetic image to visually unnoticeable level. A noise-to-sample number mapping algorithm is thus performed on each block of the noise map, higher noise value is mapped to larger sample number, and lower noise value is mapped to smaller sample number, the result of mapping is called sample map. Each pixel in a sample map can be used by Monte Carlo ray tracer to reduce the noise level in the corresponding block of pixels in a synthetic image. However, this block based scheme produces blocky artifact as appeared in video and image compression algorithms. We use Gaussian filter to smooth the sample map, the result is adaptive sample map (ASP). ASP serves two purposes in rendering process; its statistics information can be used as noise level metric in synthetic image, and it can also be used by a Monte Carlo ray tracer to refine the synthetic image adaptively in order to reduce the noise to unnoticeable level but with less rendering time than the brute force method.
Microseismic network design assessment based on 3D ray tracing
NASA Astrophysics Data System (ADS)
Näsholm, Sven Peter; Wuestefeld, Andreas; Lubrano-Lavadera, Paul; Lang, Dominik; Kaschwich, Tina; Oye, Volker
2016-04-01
There is increasing demand on the versatility of microseismic monitoring networks. In early projects, being able to locate any triggers was considered a success. These early successes led to a better understanding of how to extract value from microseismic results. Today operators, regulators, and service providers work closely together in order to find the optimum network design to meet various requirements. In the current study we demonstrate an integrated and streamlined network capability assessment approach. It is intended for use during the microseismic network design process prior to installation. The assessments are derived from 3D ray tracing between a grid of event points and the sensors. Three aspects are discussed: 1) Magnitude of completeness or detection limit; 2) Event location accuracy; and 3) Ground-motion hazard. The network capability parameters 1) and 2) are estimated at all hypothetic event locations and are presented in the form of maps given a seismic sensor coordinate scenario. In addition, the ray tracing traveltimes permit to estimate the point-spread-functions (PSFs) at the event grid points. PSFs are useful in assessing the resolution and focusing capability of the network for stacking-based event location and imaging methods. We estimate the performance for a hypothetical network case with 11 sensors. We consider the well-documented region around the San Andreas Fault Observatory at Depth (SAFOD) located north of Parkfield, California. The ray tracing is done through a detailed velocity model which covers a 26.2 by 21.2 km wide area around the SAFOD drill site with a resolution of 200 m both for the P-and S-wave velocities. Systematic network capability assessment for different sensor site scenarios prior to installation facilitates finding a final design which meets the survey objectives.
Development of 2D soft X-ray measurement system in the large helical device.
Takemura, Y; Ohdachi, S; Watanabe, K Y; Du, X D
2014-11-01
A fast two-dimensional soft X-ray camera using silicon photo diode array is being developed in order to investigate high frequency MHD instability with high mode number. The advantage of the adopted diode is a large sensor area of 10 mm × 10 mm and small diode capacitance which enable us to measure signals with the short response time. The characteristic of the prototype is summarized as follows: Channel number is 6 × 8 = 48, detection range 1∼10 keV, the spatial resolution 128 mm at the plasma location, and frequency range DC∼100 kHz. Synthetic image of the prototype in the Large Helical Device is estimated by using perturbation model of MHD mode. PMID:25430317
2D X-ray radiography of imploding capsules at the national ignition facility.
Rygg, J R; Jones, O S; Field, J E; Barrios, M A; Benedetti, L R; Collins, G W; Eder, D C; Edwards, M J; Kline, J L; Kroll, J J; Landen, O L; Ma, T; Pak, A; Peterson, J L; Raman, K; Town, R P J; Bradley, D K
2014-05-16
First measurements of the in-flight shape of imploding inertial confinement fusion (ICF) capsules at the National Ignition Facility (NIF) were obtained by using two-dimensional x-ray radiography. The sequence of area-backlit, time-gated pinhole images is analyzed for implosion velocity, low-mode shape and density asymmetries, and the absolute offset and center-of-mass velocity of the capsule shell. The in-flight shell is often observed to be asymmetric even when the concomitant core self-emission is round. A ∼ 15 μm shell asymmetry amplitude of the Y(40) spherical harmonic mode was observed for standard NIF ICF hohlraums at a shell radius of ∼ 200 μm (capsule at ∼ 5× radial compression). This asymmetry is mitigated by a ∼ 10% increase in the hohlraum length. PMID:24877944
A GPU Simulation Tool for Training and Optimisation in 2D Digital X-Ray Imaging
Gallio, Elena; Rampado, Osvaldo; Gianaria, Elena; Bianchi, Silvio Diego; Ropolo, Roberto
2015-01-01
Conventional radiology is performed by means of digital detectors, with various types of technology and different performance in terms of efficiency and image quality. Following the arrival of a new digital detector in a radiology department, all the staff involved should adapt the procedure parameters to the properties of the detector, in order to achieve an optimal result in terms of correct diagnostic information and minimum radiation risks for the patient. The aim of this study was to develop and validate a software capable of simulating a digital X-ray imaging system, using graphics processing unit computing. All radiological image components were implemented in this application: an X-ray tube with primary beam, a virtual patient, noise, scatter radiation, a grid and a digital detector. Three different digital detectors (two digital radiography and a computed radiography systems) were implemented. In order to validate the software, we carried out a quantitative comparison of geometrical and anthropomorphic phantom simulated images with those acquired. In terms of average pixel values, the maximum differences were below 15%, while the noise values were in agreement with a maximum difference of 20%. The relative trends of contrast to noise ratio versus beam energy and intensity were well simulated. Total calculation times were below 3 seconds for clinical images with pixel size of actual dimensions less than 0.2 mm. The application proved to be efficient and realistic. Short calculation times and the accuracy of the results obtained make this software a useful tool for training operators and dose optimisation studies. PMID:26545097
A GPU Simulation Tool for Training and Optimisation in 2D Digital X-Ray Imaging.
Gallio, Elena; Rampado, Osvaldo; Gianaria, Elena; Bianchi, Silvio Diego; Ropolo, Roberto
2015-01-01
Conventional radiology is performed by means of digital detectors, with various types of technology and different performance in terms of efficiency and image quality. Following the arrival of a new digital detector in a radiology department, all the staff involved should adapt the procedure parameters to the properties of the detector, in order to achieve an optimal result in terms of correct diagnostic information and minimum radiation risks for the patient. The aim of this study was to develop and validate a software capable of simulating a digital X-ray imaging system, using graphics processing unit computing. All radiological image components were implemented in this application: an X-ray tube with primary beam, a virtual patient, noise, scatter radiation, a grid and a digital detector. Three different digital detectors (two digital radiography and a computed radiography systems) were implemented. In order to validate the software, we carried out a quantitative comparison of geometrical and anthropomorphic phantom simulated images with those acquired. In terms of average pixel values, the maximum differences were below 15%, while the noise values were in agreement with a maximum difference of 20%. The relative trends of contrast to noise ratio versus beam energy and intensity were well simulated. Total calculation times were below 3 seconds for clinical images with pixel size of actual dimensions less than 0.2 mm. The application proved to be efficient and realistic. Short calculation times and the accuracy of the results obtained make this software a useful tool for training operators and dose optimisation studies. PMID:26545097
Comment on "Improved ray tracing air mass numbers model"
NASA Astrophysics Data System (ADS)
van der Werf, Siebren Y.
2008-01-01
Air mass numbers have traditionally been obtained by techniques that use height as the integration variable. This introduces an inherent singularity at the horizon, and ad hoc solutions have been invented to cope with it. A survey of the possible options including integration by height, zenith angle, and horizontal distance or path length is presented. Ray tracing by path length is shown to avoid singularities both at the horizon and in the zenith. A fourth-order Runge-Kutta numerical integration scheme is presented, which treats refraction and air mass as path integrals. The latter may optionally be split out into separate contributions of the atmosphere's constituents.
A complete ray-trace analysis of the Mirage toy
NASA Astrophysics Data System (ADS)
Adhya, Sriya; Noé, John W.
2007-06-01
The `Mirage' (Opti-Gone International) is a well-known optics demonstration (PIRA index number 6A20.35) that uses two opposed concave mirrors to project a real image of a small object into space. We studied image formation in the Mirage by standard 2x2 matrix methods and by exact ray tracing, with particular attention to additional real images that can be observed when the mirror separation is increased beyond one focal length. We find that the three readily observed secondary images correspond to 4, 6, or 8 reflections, respectively, contrary to previous reports.
Photorealistic ray tracing to visualize automobile side mirror reflective scenes.
Lee, Hocheol; Kim, Kyuman; Lee, Gang; Lee, Sungkoo; Kim, Jingu
2014-10-20
We describe an interactive visualization procedure for determining the optimal surface of a special automobile side mirror, thereby removing the blind spot, without the need for feedback from the error-prone manufacturing process. If the horizontally progressive curvature distributions are set to the semi-mathematical expression for a free-form surface, the surface point set can then be derived through numerical integration. This is then converted to a NURBS surface while retaining the surface curvature. Then, reflective scenes from the driving environment can be virtually realized using photorealistic ray tracing, in order to evaluate how these reflected images would appear to drivers. PMID:25401606
NASA Astrophysics Data System (ADS)
van der Bom, I. M. J.; Klein, S.; Staring, M.; Homan, R.; Bartels, L. W.; Pluim, J. P. W.
2011-03-01
The advantage of 2D-3D image registration methods versus direct image-to-patient registration, is that these methods generally do not require user interaction (such as manual annotations), additional machinery or additional acquisition of 3D data. A variety of intensity-based similarity measures has been proposed and evaluated for different applications. These studies showed that the registration accuracy and capture range are influenced by the choice of similarity measure. However, the influence of the optimization method on intensity-based 2D-3D image registration has not been investigated. We have compared the registration performance of seven optimization methods in combination with three similarity measures: gradient difference, gradient correlation, and pattern intensity. Optimization methods included in this study were: regular step gradient descent, Nelder-Mead, Powell-Brent, Quasi-Newton, nonlinear conjugate gradient, simultaneous perturbation stochastic approximation, and evolution strategy. Registration experiments were performed on multiple patient data sets that were obtained during cerebral interventions. Various component combinations were evaluated on registration accuracy, capture range, and registration time. The results showed that for the same similarity measure, different registration accuracies and capture ranges were obtained when different optimization methods were used. For gradient difference, largest capture ranges were obtained with Powell-Brent and simultaneous perturbation stochastic approximation. Gradient correlation and pattern intensity had the largest capture ranges in combination with Powell-Brent, Nelder-Mead, nonlinear conjugate gradient, and Quasi-Newton. Average registration time, expressed in the number of DRRs required for convergence, was the lowest for Powell-Brent. Based on these results, we conclude that Powell-Brent is a reliable optimization method for intensity-based 2D-3D registration of x-ray images to CBCT
Ray-tracing software comparison for linear focusing solar collectors
NASA Astrophysics Data System (ADS)
Osório, Tiago; Horta, Pedro; Larcher, Marco; Pujol-Nadal, Ramón; Hertel, Julian; van Rooyen, De Wet; Heimsath, Anna; Schneider, Simon; Benitez, Daniel; Frein, Antoine; Denarie, Alice
2016-05-01
Ray-Tracing software tools have been widely used in the optical design of solar concentrating collectors. In spite of the ability of these tools to assess the geometrical and material aspects impacting the optical performance of concentrators, their use in combination with experimental measurements in the framework of collector testing procedures as not been implemented, to the date, in none of the current solar collector testing standards. In the latest revision of ISO9806 an effort was made to include linear focusing concentrating collectors but some practical and theoretical difficulties emerged. A Ray-Tracing analysis could provide important contributions to overcome these issues, complementing the experimental results obtained through thermal testing and allowing the achievement of more thorough testing outputs with lower experimental requirements. In order to evaluate different available software tools a comparison study was conducted. Taking as representative technologies for line-focus concentrators the Parabolic Trough Collector and the Linear Fresnel Reflector Collector, two exemplary cases with predefined conditions - geometry, sun model and material properties - were simulated with different software tools. This work was carried out within IEA/SHC Task 49 "Solar Heat Integration in Industrial Processes".
Ray Tracing Modeling of Gravity Wave Propagation and Dissipation
NASA Astrophysics Data System (ADS)
Vadas, Sharon; Crowley, Geoff
In this paper, we describe a ray trace model which calculates the wavevector, location and phase of a gravity wave (GW) as it propagates in the lower atmosphere and thermosphere. If used for a discreet transient source (such as a deep convective plume), we describe how this model can calculate the body forcing and the heat/cooling that are created when the GWs within a wave packet dissipate in the thermosphere from kinematic viscosity and thermal diffusivity. Although the body force calculation requires only the divergence of the momentum flux, the heat/cooling calculation requires the reconstructed GW field (e.g., density, velocity perturbations), which in turn requires the GW dissipative polarization relations. We describe these relations. We then describe the results of a recent study involving GWs identified from TIDDBIT HF Doppler sounder data taken at Wallops Island, VI, USA. Using this ray trace model, we determine if the unusual neutral wind profile measured by a rocket experiment at high altitudes (~290-370 km) could have been caused by the propagation and dissipation of several waves observed by TIDDBIT at lower altitudes.
Personalized x-ray reconstruction of the proximal femur via a non-rigid 2D-3D registration
NASA Astrophysics Data System (ADS)
Yu, Weimin; Zysset, Philippe; Zheng, Guoyan
2015-03-01
In this paper we present a new approach for a personalized X-ray reconstruction of the proximal femur via a non-rigid registration of a 3D volumetric template to 2D calibrated C-arm images. The 2D-3D registration is done with a hierarchical two-stage strategy: the global scaled rigid registration stage followed by a regularized deformable b-spline registration stage. In both stages, a set of control points with uniform spacing are placed over the domain of the 3D volumetric template and the registrations are driven by computing updated positions of these control points, which then allows to accurately register the 3D volumetric template to the reference space of the C-arm images. Comprehensive experiments on simulated images, on images of cadaveric femurs and on clinical datasets are designed and conducted to evaluate the performance of the proposed approach. Quantitative and qualitative evaluation results are given, which demonstrate the efficacy of the present approach.
NASA Astrophysics Data System (ADS)
Engelbrecht, N. E.; Burger, R. A.
2015-12-01
In this study, a novel ab initio cosmic ray (CR) modulation code that solves a set of stochastic transport equations equivalent to the Parker transport equation, and that uses output from a turbulence transport code as input for the diffusion tensor, is introduced. This code is benchmarked with a previous approach to ab initio modulation. The sensitivity of computed galactic CR proton spectra at Earth to assumptions made as to the low-wavenumber behavior of the two-dimensional (2D) turbulence power spectrum is investigated using perpendicular mean free path expressions derived from two different scattering theories. Constraints on the low-wavenumber behavior of the 2D power spectrum are inferred from the qualitative comparison of computed CR spectra with spacecraft observations at Earth. Another key difference from previous studies is that observed and inferred CR intensity spectra at 73 AU are used as boundary spectra instead of the usual local interstellar spectrum. Furthermore, the results presented here provide a tentative explanation as to the reason behind the unusually high galactic proton intensity spectra observed in 2009 during the recent unusual solar minimum.
Wolter X-Ray Microscope Computed Tomography Ray-Trace Model with Preliminary Simulation Results
Jackson, J A
2006-02-27
It is proposed to build a Wolter X-ray Microscope Computed Tomography System in order to characterize objects to sub-micrometer resolution. Wolter Optics Systems use hyperbolic, elliptical, and/or parabolic mirrors to reflect x-rays in order to focus or magnify an image. Wolter Optics have been used as telescopes and as microscopes. As microscopes they have been used for a number of purposes such as measuring emission x-rays and x-ray fluoresce of thin biological samples. Standard Computed Tomography (CT) Systems use 2D radiographic images, from a series of rotational angles, acquired by passing x-rays through an object to reconstruct a 3D image of the object. The x-ray paths in a Wolter X-ray Microscope will be considerably different than those of a standard CT system. There is little information about the 2D radiographic images that can be expected from such a system. There are questions about the quality, resolution and focusing range of an image created with such a system. It is not known whether characterization information can be obtained from these images and whether these 2D images can be reconstructed to 3D images of the object. A code has been developed to model the 2D radiographic image created by an object in a Wolter X-ray Microscope. This code simply follows the x-ray through the object and optics. There is no modeling at this point of other effects, such as scattering, reflection losses etc. Any object, of appropriate size, can be used in the model code. A series of simulations using a number of different objects was run to study the effects of the optics. The next step will be to use this model to reconstruct an object from the simulated data. Funding for the project ended before this goal could be accomplished. The following documentation includes: (1) background information on current X-ray imaging systems, (2) background on Wolter Optics, (3) description of the Wolter System being used, (4) purpose, limitations and development of the modeling
Weeratunga, S K
2008-11-06
Ares and Kull are mature code frameworks that support ALE hydrodynamics for a variety of HEDP applications at LLNL, using two widely different meshing approaches. While Ares is based on a 2-D/3-D block-structured mesh data base, Kull is designed to support unstructured, arbitrary polygonal/polyhedral meshes. In addition, both frameworks are capable of running applications on large, distributed-memory parallel machines. Currently, both these frameworks separately support assorted collections of physics packages related to HEDP, including one for the energy deposition by laser/ion-beam ray tracing. This study analyzes the options available for developing a common laser/ion-beam ray tracing package that can be easily shared between these two code frameworks and concludes with a set of recommendations for its development.
Ray-traced tropospheric total slant delays for GNSS processing
NASA Astrophysics Data System (ADS)
Hobiger, T.; Ichikawa, R.; Hatanaka, Y.; Yutsudo, T.; Iwashita, C.; Miyahara, B.; Koyama, Y.; Kondo, T.
2007-12-01
Numerical weather models have undergone an improvement of spatial and temporal resolution in the recent years, which made their use for GNSS applications feasible. Ray-tracing through such models permits the computation of total troposphere delays and ray-bending angles. At the National Institute of Information and Communications Technology (NICT), Japan the so-called KAshima RAy-tracing Tools (KARAT) have been developed which allow to obtain troposphere delay corrections in real-time. Together with fine-mesh weather models from the Japanese Meteorological Agency (JMA) huge parts of the East Asian region, including Japan, Korea, Taiwan and East China, can be covered. The Japanese GEONET with its more than 1300 GNSS receivers represent an ideal test-bed for the evaluation of the performance of KARAT. In cooperation with the Geographical Survey Institute (GSI), Japan more than 1.6 billion observations, covering measurements from July 1st until August 31st, 2006, were processed and the corresponding troposphere delays were used to modify the original RINEX files by subtraction of code- and phase delays. These modified observations were processed by a dedicated analysis run of the GEONET operation center, taking advantage of the computer cluster at GSI. First results from this study, together with an in-depth discussion about the assets and drawbacks of the reduction of troposphere total slant delays will be given in this presentation. Additionally an overview about KARAT, the treatment of observational data and the impact of future refined numerical weather models on GNSS analysis will be included in this contribution.
High performance dosimetry calculations using adapted ray-tracing
NASA Astrophysics Data System (ADS)
Perrotte, Lancelot; Saupin, Guillaume
2010-11-01
When preparing interventions on nuclear sites, it is interesting to study different scenarios, to identify the most appropriate one for the operator(s). Using virtual reality tools is a good way to simulate the potential scenarios. Thus, taking advantage of very efficient computation times can help the user studying different complex scenarios, by immediately evaluating the impact of any changes. In the field of radiation protection, people often use computation codes based on the straight line attenuation method with build-up factors. As for other approaches, geometrical computations (finding all the interactions between radiation rays and the scene objects) remain the bottleneck of the simulation. We present in this paper several optimizations used to speed up these geometrical computations, using innovative GPU ray-tracing algorithms. For instance, we manage to compute every intersectionbetween 600 000 rays and a huge 3D industrial scene in a fraction of second. Moreover, our algorithm works the same way for both static and dynamic scenes, allowing easier study of complex intervention scenarios (where everything moves: the operator(s), the shielding objects, the radiation sources).
Ray Tracing to Predict Optical Behaviour of Shock Compressed Dielectrics
NASA Astrophysics Data System (ADS)
Tear, Gareth R.; Proud, William G.
2015-06-01
In order to investigate the optical response of dielectric materials under shock compression, a characteristics model has been combined with a three dimensional optical ray tracing model. A general biaxial optical model is used along with a first order photoelastic model which couples the characteristics component to the optical component. This optical model is three dimensional and as such can be used to investigate small deviations from the perfect one dimensional shock wave which is typically assumed in plate impact experiments. A detailed description of the model will be presented, and comparison to available literature as well as recent experiments on the optical behaviour of shock compressed a-cut calcite and a-cut sapphire. The authors would like to thank Dr D E Eakins and Dr D J Chapman for fruitful discussions. The Institute of Shock Physics acknowledges the continued support of AWE and Imperial College London.
Ray tracing in FLRW flat space-times
NASA Astrophysics Data System (ADS)
Acquaviva, Giovanni; Bonetti, Luca; Cognola, Guido; Zerbini, Sergio
2013-12-01
In this work we take moves from the debate triggered by Melia et al. in [J. Cosmol. Astropart. Phys. 09 (2012) 029; Mon. Not. R. Astron. Soc. 421, 3356 (2012)] and followed by opposite comments by Lewis and Oirschot in [Mon. Not. R. Astron. Soc. Lett. 423, 26 (2012); 431, 25 (2013)]. The point in question regards the role of the Hubble horizon as a limit for observability in a cosmological setting. We propose to tackle the issue in a broader way by relating it to the causal character of the Hubble surface and to the tracing of null trajectories, focusing on both three-fluids and generalized Chaplygin gas models. The results should make clear that for quite reasonable and physically motivated models, light rays reaching a comoving observer at R(t0)=0 have never traveled a distance greater than the proper radius of the horizon until t0.
Calculation of material properties and ray tracing in transformation media.
Schurig, D; Pendry, J B; Smith, D R
2006-10-16
Complex and interesting electromagnetic behavior can be found in spaces with non-flat topology. When considering the properties of an electromagnetic medium under an arbitrary coordinate transformation an alternative interpretation presents itself. The transformed material property tensors may be interpreted as a different set of material properties in a flat, Cartesian space. We describe the calculation of these material properties for coordinate transformations that describe spaces with spherical or cylindrical holes in them. The resulting material properties can then implement invisibility cloaks in flat space. We also describe a method for performing geometric ray tracing in these materials which are both inhomogeneous and anisotropic in their electric permittivity and magnetic permeability. PMID:19529371
Polarization Ray Trace Model of the MODIS Instrument
NASA Technical Reports Server (NTRS)
Waluschka, Eugene; Xiong, Jack; Esaias, Wayne E.; Voss, Kenneth; Souaidia, Nordine; Pellicori, Samuel; Moyer, David; Guenther, Bruce; Barnes, William
2004-01-01
Sunlight reflected from the earth is, to a certain extent, polarized. Radiometers, such as the MODIS instrument on board the TERRA and AQUA spacecraft, are to a certain extent polarizers. Accurate radiometric measurements must take into account both the polarization state of the scene and the polarization sensitivity of the measuring instrument. The measured polarization characteristics of the MODIS instruments are contained in various radiometric models. Continued use of these radiometric math models, over a number of years, have shown where these models can be improved. Currently a MODIS polarization ray trace model has been created which models the thin film structure on the optical elements. This approach is described and modeled and measured instrument polarization sensitivity results presented.
A FORMALISM FOR COVARIANT POLARIZED RADIATIVE TRANSPORT BY RAY TRACING
Gammie, Charles F.; Leung, Po Kin
2012-06-20
We write down a covariant formalism for polarized radiative transfer appropriate for ray tracing through a turbulent plasma. The polarized radiation field is represented by the polarization tensor (coherency matrix) N{sup {alpha}{beta}} {identical_to} (a{sup {alpha}}{sub k} a*{sup {beta}}{sub k}), where a{sub k} is a Fourier coefficient for the vector potential. Using Maxwell's equations, the Liouville-Vlasov equation, and the WKB approximation, we show that the transport equation in vacuo is k{sup {mu}}{nabla}{sub {mu}} N{sup {alpha}{beta}} = 0. We show that this is equivalent to Broderick and Blandford's formalism based on invariant Stokes parameters and a rotation coefficient, and suggest a modification that may reduce truncation error in some situations. Finally, we write down several alternative approaches to integrating the transfer equation.
Ray trace calculation of ionospheric propagation at lower frequencies
NASA Astrophysics Data System (ADS)
Reilly, Michael H.
2006-10-01
The Raytrace/Ionospheric Conductivity and Electron Density-Bent-Gallagher model has been revised to make it applicable to ionospheric propagation at low radio frequencies (0.5-5.0 MHz), where the ionosphere and magnetic anisotropy drastically alter propagation paths and provide a severe test of propagation model algorithms. The necessary revisions are discussed, and the model is applied to the problem of ionospheric penetration from a source below the ionosphere to a receiver above the ionosphere. It is necessary to include the electron collision frequency in the Appleton-Hartree index of refraction in order to permit ionospheric penetration for radio frequencies below the maximum plasma frequency (e.g., whistler modes). The associated reformulation of the ray trace equations for a complex index of refraction is straightforward. Difficulties with numerical methods are cited for the lowest frequencies, and future improvements are indicated.
Ionospheric Plasma Drift Analysis Technique Based On Ray Tracing
NASA Astrophysics Data System (ADS)
Ari, Gizem; Toker, Cenk
2016-07-01
Ionospheric drift measurements provide important information about the variability in the ionosphere, which can be used to quantify ionospheric disturbances caused by natural phenomena such as solar, geomagnetic, gravitational and seismic activities. One of the prominent ways for drift measurement depends on instrumentation based measurements, e.g. using an ionosonde. The drift estimation of an ionosonde depends on measuring the Doppler shift on the received signal, where the main cause of Doppler shift is the change in the length of the propagation path of the signal between the transmitter and the receiver. Unfortunately, ionosondes are expensive devices and their installation and maintenance require special care. Furthermore, the ionosonde network over the world or even Europe is not dense enough to obtain a global or continental drift map. In order to overcome the difficulties related to an ionosonde, we propose a technique to perform ionospheric drift estimation based on ray tracing. First, a two dimensional TEC map is constructed by using the IONOLAB-MAP tool which spatially interpolates the VTEC estimates obtained from the EUREF CORS network. Next, a three dimensional electron density profile is generated by inputting the TEC estimates to the IRI-2015 model. Eventually, a close-to-real situation electron density profile is obtained in which ray tracing can be performed. These profiles can be constructed periodically with a period of as low as 30 seconds. By processing two consequent snapshots together and calculating the propagation paths, we estimate the drift measurements over any coordinate of concern. We test our technique by comparing the results to the drift measurements taken at the DPS ionosonde at Pruhonice, Czech Republic. This study is supported by TUBITAK 115E915 and Joint TUBITAK 114E092 and AS CR14/001 projects.
Mapping gray-scale image to 3D surface scanning data by ray tracing
NASA Astrophysics Data System (ADS)
Li, Peng; Jones, Peter R. M.
1997-03-01
The extraction and location of feature points from range imaging is an important but difficult task in machine vision based measurement systems. There exist some feature points which are not able to be detected from pure geometric characteristics, particularly in those measurement tasks related to the human body. The Loughborough Anthropometric Shadow Scanner (LASS) is a whole body surface scanner based on structured light technique. Certain applications of LASS require accurate location of anthropometric landmarks from the scanned data. This is sometimes impossible from existing raw data because some landmarks do not appear in the scanned data. Identification of these landmarks has to resort to surface texture of the scanned object. Modifications to LASS were made to allow gray-scale images to be captured before or after the object was scanned. Two-dimensional gray-scale image must be mapped to the scanned data to acquire the 3D coordinates of a landmark. The method to map 2D images to the scanned data is based on the colinearity conditions and ray-tracing method. If the camera center and image coordinates are known, the corresponding object point must lie on a ray starting from the camera center and connecting to the image coordinate. By intersecting the ray with the scanned surface of the object, the 3D coordinates of a point can be solved. Experimentation has demonstrated the feasibility of the method.
Bi-planar 2D-to-3D registration in Fourier domain for stereoscopic x-ray motion tracking
NASA Astrophysics Data System (ADS)
Zosso, Dominique; Le Callennec, Benoît; Bach Cuadra, Meritxell; Aminian, Kamiar; Jolles, Brigitte M.; Thiran, Jean-Philippe
2008-03-01
In this paper we present a new method to track bone movements in stereoscopic X-ray image series of the knee joint. The method is based on two different X-ray image sets: a rotational series of acquisitions of the still subject knee that allows the tomographic reconstruction of the three-dimensional volume (model), and a stereoscopic image series of orthogonal projections as the subject performs movements. Tracking the movements of bones throughout the stereoscopic image series means to determine, for each frame, the best pose of every moving element (bone) previously identified in the 3D reconstructed model. The quality of a pose is reflected in the similarity between its theoretical projections and the actual radiographs. We use direct Fourier reconstruction to approximate the three-dimensional volume of the knee joint. Then, to avoid the expensive computation of digitally rendered radiographs (DRR) for pose recovery, we develop a corollary to the 3-dimensional central-slice theorem and reformulate the tracking problem in the Fourier domain. Under the hypothesis of parallel X-ray beams, the heavy 2D-to-3D registration of projections in the signal domain is replaced by efficient slice-to-volume registration in the Fourier domain. Focusing on rotational movements, the translation-relevant phase information can be discarded and we only consider scalar Fourier amplitudes. The core of our motion tracking algorithm can be implemented as a classical frame-wise slice-to-volume registration task. Results on both synthetic and real images confirm the validity of our approach.
Calibration model of a dual gain flat panel detector for 2D and 3D x-ray imaging
Schmidgunst, C.; Ritter, D.; Lang, E.
2007-09-15
The continuing research and further development in flat panel detector technology have led to its integration into more and more medical x-ray systems for two-dimensional (2D) and three-dimensional (3D) imaging, such as fixed or mobile C arms. Besides the obvious advantages of flat panel detectors, like the slim design and the resulting optimum accessibility to the patient, their success is primarily a product of the image quality that can be achieved. The benefits in the physical and performance-related features as opposed to conventional image intensifier systems (e.g., distortion-free reproduction of imaging information or almost linear signal response over a large dynamic range) can be fully exploited, however, only if the raw detector images are correctly calibrated and postprocessed. Previous procedures for processing raw data contain idealizations that, in the real world, lead to artifacts or losses in image quality. Thus, for example, temperature dependencies or changes in beam geometry, as can occur with mobile C arm systems, have not been taken into account up to this time. Additionally, adverse characteristics such as image lag or aging effects have to be compensated to attain the best possible image quality. In this article a procedure is presented that takes into account the important dependencies of the individual pixel sensitivity of flat panel detectors used in 2D or 3D imaging and simultaneously minimizes the work required for an extensive recalibration. It is suitable for conventional detectors with only one gain mode as well as for the detectors specially developed for 3D imaging with dual gain read-out technology.
Image transfer through cirrus clouds. I. Ray trace analysis and wave-front reconstruction.
Landesman, B T; Kindilien, P J; Matson, C L; Caudill, T R
2000-10-20
A new technique for modeling image transfer through cirrus clouds is presented. The technique uses a ray trace to model beam propagation through a three-dimensional volume of polydisperse, hexagonal ice crystals. Beyond the cloud, the technique makes use of standard Huygens-Fresnel propagation methods. At the air-cloud interface, each wave front is resolved into a ray distribution for input to the ray trace software. Similarly, a wave front is reconstructed from the output ray distribution at the cloud-air interface. Simulation output from the ray trace program is presented and the modulation transfer function for stars imaged through cirrus clouds of varying depths is discussed. PMID:18354542
General Relativistic Ray Tracing for X-ray Reverberation and Polarimetry Studies of Black Holes
NASA Astrophysics Data System (ADS)
Hoormann, Janie; Krawczynski, Henric
2015-01-01
We present the results of General Relativistic (GR) ray tracing calculations of the X-ray emission from mass accreting stellar mass and supermassive black holes. Our study aims at exploring the X-ray reverberation and X-ray polarimetry signatures of different accretion flow geometries and different spacetime backgrounds (GR and non-GR backgrounds). We present first results derived for the well-known lamp-post model, where a point source of continuum emission illuminates an accretion disk with high energy photons which are tracked by parallel transporting the photon wave and polarization vectors. The simulation code models the reprocessing and reflection by of photons impinging on the accretion disk. We study the degeneracy of astrophysical parameters (parametrizing the geometry of the accretion disk and the location and properties of the lamppost photon source) and the parameters describing the underlying metrics. We emphasize furthermore the difference of the observational signatures for stellar mass and supermassive black holes.
RAY-RAMSES: a code for ray tracing on the fly in N-body simulations
NASA Astrophysics Data System (ADS)
Barreira, Alexandre; Llinares, Claudio; Bose, Sownak; Li, Baojiu
2016-05-01
We present a ray tracing code to compute integrated cosmological observables on the fly in AMR N-body simulations. Unlike conventional ray tracing techniques, our code takes full advantage of the time and spatial resolution attained by the N-body simulation by computing the integrals along the line of sight on a cell-by-cell basis through the AMR simulation grid. Moroever, since it runs on the fly in the N-body run, our code can produce maps of the desired observables without storing large (or any) amounts of data for post-processing. We implemented our routines in the RAMSES N-body code and tested the implementation using an example of weak lensing simulation. We analyse basic statistics of lensing convergence maps and find good agreement with semi-analytical methods. The ray tracing methodology presented here can be used in several cosmological analysis such as Sunyaev-Zel'dovich and integrated Sachs-Wolfe effect studies as well as modified gravity. Our code can also be used in cross-checks of the more conventional methods, which can be important in tests of theory systematics in preparation for upcoming large scale structure surveys.
Laser ray tracing in a parallel arbitrary Lagrangian-Eulerian adaptive mesh refinement hydrocode
NASA Astrophysics Data System (ADS)
Masters, N. D.; Kaiser, T. B.; Anderson, R. W.; Eder, D. C.; Fisher, A. C.; Koniges, A. E.
2010-08-01
ALE-AMR is a new hydrocode that we are developing as a predictive modeling tool for debris and shrapnel formation in high-energy laser experiments. In this paper we present our approach to implementing laser ray tracing in ALE-AMR. We present the basic concepts of laser ray tracing and our approach to efficiently traverse the adaptive mesh hierarchy.
Fast Ray Tracing of Lunar Digital Elevation Models
NASA Technical Reports Server (NTRS)
McClanahan, Timothy P.; Evans, L. G.; Starr, R. D.; Mitrofanov, I.
2009-01-01
Ray-tracing (RT) of Lunar Digital Elevation Models (DEM)'s is performed to virtually derive the degree of radiation incident to terrain as a function of time, orbital and ephemeris constraints [I- 4]. This process is an integral modeling process in lunar polar research and exploration due to the present paucity of terrain information at the poles and mission planning activities for the anticipated spring 2009 launch of the Lunar Reconnaissance Orbiter (LRO). As part of the Lunar Exploration Neutron Detector (LEND) and Lunar Crater Observation and Sensing Satellite (LCROSS) preparations RI methods are used to estimate the critical conditions presented by the combined effects of high latitude, terrain and the moons low obliquity [5-7]. These factors yield low incident solar illumination and subsequently extreme thermal, and radiation conditions. The presented research uses RT methods both for radiation transport modeling in space and regolith related research as well as to derive permanently shadowed regions (PSR)'s in high latitude topographic minima, e.g craters. These regions are of scientific and human exploration interest due to the near constant low temperatures in PSRs, inferred to be < 100 K. Hydrogen is thought to have accumulated in PSR's through the combined effects of periodic cometary bombardment and/or solar wind processes, and the extreme cold which minimizes hydrogen sublimation [8-9]. RT methods are also of use in surface position optimization for future illumination dependent on surface resources e.g. power and communications equipment.
NASA Astrophysics Data System (ADS)
Atzeni, Stefano; Marocchino, Alberto; Schiavi, Angelo
2016-03-01
Accurate descriptions of laser power coupling to the plasma and electron energy transport are crucial for designing shock-ignition targets and assessing their robustness (in particular with regard to laser and positioning errors). To this purpose, the 2D DUED laser fusion code has been improved with the inclusion of a 3D laser ray-tracing scheme and a model for non-local electron transport. 2D simulations with the upgraded code are presented; the dependence of the fusion yield vs target displacement is studied. Two different irradiation configurations are considered.
Ray-tracing simulations of coupled dark energy models
NASA Astrophysics Data System (ADS)
Pace, Francesco; Baldi, Marco; Moscardini, Lauro; Bacon, David; Crittenden, Robert
2015-02-01
Dark matter and dark energy are usually assumed to couple only gravitationally. An extension to this picture is to model dark energy as a scalar field coupled directly to cold dark matter. This coupling leads to new physical effects, such as a fifth force and a time-dependent dark matter particle mass. In this work we examine the impact that coupling has on weak lensing statistics by constructing realistic simulated weak lensing maps using ray-tracing techniques through N-body cosmological simulations. We construct maps for different lensing quantities, covering a range of scales from a few arcminutes to several degrees. The concordance Λ cold dark matter (ΛCDM) model is compared to different coupled dark energy models, described either by an exponential scalar field potential (standard coupled dark energy scenario) or by a SUGRA potential (bouncing model). We analyse several statistical quantities and our results, with sources at low redshifts are largely consistent with previous work on cosmic microwave background lensing by Carbone et al. The most significant differences from the ΛCDM model are due to the enhanced growth of the perturbations and to the effective friction term in non-linear dynamics. For the most extreme models, we see differences in the power spectra up to 40 per cent compared to the ΛCDM model. The different time evolution of the linear matter overdensity can account for most of the differences, but when controlling for this using a ΛCDM model having the same normalization, the overall signal is smaller due to the effect of the friction term appearing in the equation of motion for dark matter particles.
NASA Astrophysics Data System (ADS)
Kaiser, Markus; John, Matthias; Borsdorf, Anja; Mountney, Peter; Ionasec, Razvan; Nöttling, Alois; Kiefer, Philipp; Seeburger, Jörg; Neumuth, Thomas
2013-03-01
For transcatheter-based minimally invasive procedures in structural heart disease ultrasound and X-ray are the two enabling imaging modalities. A live fusion of both real-time modalities can potentially improve the workflow and the catheter navigation by combining the excellent instrument imaging of X-ray with the high-quality soft tissue imaging of ultrasound. A recently published approach to fuse X-ray fluoroscopy with trans-esophageal echo (TEE) registers the ultrasound probe to X-ray images by a 2D-3D registration method which inherently provides a registration of ultrasound images to X-ray images. In this paper, we significantly accelerate the 2D-3D registration method in this context. The main novelty is to generate the projection images (DRR) of the 3D object not via volume ray-casting but instead via a fast rendering of triangular meshes. This is possible, because in the setting for TEE/X-ray fusion the 3D geometry of the ultrasound probe is known in advance and their main components can be described by triangular meshes. We show that the new approach can achieve a speedup factor up to 65 and does not affect the registration accuracy when used in conjunction with the gradient correlation similarity measure. The improvement is independent of the underlying registration optimizer. Based on the results, a TEE/X-ray fusion could be performed with a higher frame rate and a shorter time lag towards real-time registration performance. The approach could potentially accelerate other applications of 2D-3D registrations, e.g. the registration of implant models with X-ray images.
Tabuchi, M.; Tatsumi, M.; Yamamoto, A.; Endo, T.
2013-07-01
A new correction model for ray tracing of the method of characteristics is proposed in order to reduce discretization error. As the ray tracing parameters such as azimuthal angle division, polar angle division and ray separation are considered in this study. In the method of characteristics, region average scalar fluxes can be implicitly expressed by collision probabilities, although these collision probabilities are not directly treated in the ordinary calculation scheme. From this viewpoint, difference between a coarse ray tracing condition and a detailed one can be interpreted as the difference in the estimation of collision probabilities. In other words, the discretization error for ray tracing can be recognized as a consequence of inaccurate collision probabilities caused by coarse ray tracing. This discussion suggests that accurate region average scalar flux can be obtained through an appropriate correction on collision probabilities. In this paper, a correction model on collision probabilities is theoretically derived based on the neutron balance equation, and its validity is confirmed through typical single assembly calculations. The effectiveness of the present correction method is also discussed in this paper. It is confirmed that discretization error for ray tracing can be significantly reduced by the present correction method in a multi-assembly calculation, though the correction factor is estimated in single assembly geometry. (authors)
NASA Astrophysics Data System (ADS)
Vidal, F.; de Assis, J. T.; Lopes, R. T.; Lima, I.
2014-02-01
In recent years, bone quantification led to a deeper knowledge of the 3D microarchitecture. In this study the bone architecture of rats was investigated based on 2D/3D morphometric analysis using microcomputed tomography, aiming at determining the effect of the image acquisition pixel on the quality of some 2D/3D morphometric parameters, such as porosity and trabecular density.Six pairs of bone samples were used and the scans were carried out using high microcomputed tomography system, operating at three different pixel sizes of 33.3 μm, 15.0 μm and 9.5 μm. The results showed 2D parameters values lower than those obtained in the 3D analysis, mainly for trabecular density, separation and thickness.
Three-dimensional ray tracing on Delaunay-based reconstructed surfaces.
Ortiz, Sergio; Siedlecki, Damian; Remon, Laura; Marcos, Susana
2009-07-10
A method of ray tracing for free-form optical surfaces has been developed. The ray tracing through such surfaces is based on Delaunay triangulation of the discrete data of the surface and is related to finite-element modeling. Some numerical examples of applications to analytical, noisy, and experimental free-form surfaces (in particular, a corneal topography map) are presented. Ray-tracing results (i.e., spot diagram root-mean-square error) with the new method are in agreement with those obtained using a modal fitting of the surface, for sampling densities higher than 40 x 40 elements. The method competes in flexibility, simplicity, and computing times with standard methods for surface fitting and ray tracing. PMID:19593339
Harmon, Katherine J.; Bennett, Eric E.; Gomella, Andrew A.; Wen, Han
2014-01-01
The ability to map the phase distribution and lateral coherence of an x-ray wavefront offers the potential for imaging the human body through phase contrast, without the need to deposit significant radiation energy. The classic means to achieve this goal is structured illumination, in which a periodic intensity modulation is introduced into the image, and changes in the phase distribution of the wavefront are detected as distortions of the modulation pattern. Two-dimensional periodic patterns are needed to fully characterize a transverse wavefront. Traditionally, the information in a 2D pattern is retrieved at high resolution by acquiring multiple images while shifting the pattern over a 2D matrix of positions. Here we describe a method to decode 2D periodic patterns with single-axis phase stepping, without either a loss of information or increasing the number of sampling steps. The method is created to reduce the instrumentation complexity of high-resolution 2D wavefront sensing in general. It is demonstrated with motionless electromagnetic phase stepping and a flexible processing algorithm in x-ray dark-field and phase contrast imaging. PMID:24489853
Harmon, Katherine J; Bennett, Eric E; Gomella, Andrew A; Wen, Han
2014-01-01
The ability to map the phase distribution and lateral coherence of an x-ray wavefront offers the potential for imaging the human body through phase contrast, without the need to deposit significant radiation energy. The classic means to achieve this goal is structured illumination, in which a periodic intensity modulation is introduced into the image, and changes in the phase distribution of the wavefront are detected as distortions of the modulation pattern. Two-dimensional periodic patterns are needed to fully characterize a transverse wavefront. Traditionally, the information in a 2D pattern is retrieved at high resolution by acquiring multiple images while shifting the pattern over a 2D matrix of positions. Here we describe a method to decode 2D periodic patterns with single-axis phase stepping, without either a loss of information or increasing the number of sampling steps. The method is created to reduce the instrumentation complexity of high-resolution 2D wavefront sensing in general. It is demonstrated with motionless electromagnetic phase stepping and a flexible processing algorithm in x-ray dark-field and phase contrast imaging. PMID:24489853
Paraxial ray-tracing approach for the simulation of ultrasonic inspection of welds
Gardahaut, Audrey; Jezzine, Karim; Cassereau, Didier
2014-02-18
On-site inspection of bimetallic or austenitic welds can be very difficult to interpret owing to their internal structures. Skewing and splitting of the ultrasonic beam may occur due to the anisotropic and inhomogeneous properties of the welding material. In this paper, we present a ray-based method to simulate the propagation of ultrasonic waves in such structures. The formalism is based on dynamic ray tracing system in Cartesian coordinates along a reference ray. Standard ray tracing consists in the solution of a system of linear ordinary differential equations of the first order and is used to determine the trajectory of the ray. Likewise, dynamic ray tracing (DRT) also called paraxial ray tracing consists in the solution of an additional system of linear ordinary differential equations along the ray allowing paraxial quantities to be computed. It is used to evaluate the geometrical spreading and amplitude along the ray and in its vicinity. DRT is applied on a smooth representation of the elastic properties of the weld obtained thanks to an image processing technique applied on a macrograph of the weld. Simulation results are presented and compared to finite elements and experimental results.
Paraxial ray-tracing approach for the simulation of ultrasonic inspection of welds
NASA Astrophysics Data System (ADS)
Gardahaut, Audrey; Jezzine, Karim; Cassereau, Didier
2014-02-01
On-site inspection of bimetallic or austenitic welds can be very difficult to interpret owing to their internal structures. Skewing and splitting of the ultrasonic beam may occur due to the anisotropic and inhomogeneous properties of the welding material. In this paper, we present a ray-based method to simulate the propagation of ultrasonic waves in such structures. The formalism is based on dynamic ray tracing system in Cartesian coordinates along a reference ray. Standard ray tracing consists in the solution of a system of linear ordinary differential equations of the first order and is used to determine the trajectory of the ray. Likewise, dynamic ray tracing (DRT) also called paraxial ray tracing consists in the solution of an additional system of linear ordinary differential equations along the ray allowing paraxial quantities to be computed. It is used to evaluate the geometrical spreading and amplitude along the ray and in its vicinity. DRT is applied on a smooth representation of the elastic properties of the weld obtained thanks to an image processing technique applied on a macrograph of the weld. Simulation results are presented and compared to finite elements and experimental results.
Labate, L; Köster, P; Levato, T; Gizzi, L A
2012-10-01
A novel x-ray diagnostic of laser-fusion plasmas is described, allowing 2D monochromatic images of hot, dense plasmas to be obtained in any x-ray photon energy range, over a large domain, on a single-shot basis. The device (named energy-encoded pinhole camera) is based upon the use of an array of many pinholes coupled to a large area CCD camera operating in the single-photon mode. The available x-ray spectral domain is only limited by the quantum efficiency of scientific-grade x-ray CCD cameras, thus extending from a few keV up to a few tens of keV. Spectral 2D images of the emitting plasma can be obtained at any x-ray photon energy provided that a sufficient number of photons had been collected at the desired energy. Results from recent inertial confinement fusion related experiments will be reported in order to detail the new diagnostic. PMID:23126763
Ray tracing a three-dimensional scene using a hierarchical data structure
Wald, Ingo; Boulos, Solomon; Shirley, Peter
2012-09-04
Ray tracing a three-dimensional scene made up of geometric primitives that are spatially partitioned into a hierarchical data structure. One example embodiment is a method for ray tracing a three-dimensional scene made up of geometric primitives that are spatially partitioned into a hierarchical data structure. In this example embodiment, the hierarchical data structure includes at least a parent node and a corresponding plurality of child nodes. The method includes a first act of determining that a first active ray in the packet hits the parent node and a second act of descending to each of the plurality of child nodes.
Laser Ray Tracing in a Parallel Arbitrary Lagrangian-Eulerian Adaptive Mesh Refinement Hydrocode
Masters, N D; Kaiser, T B; Anderson, R W; Eder, D C; Fisher, A C; Koniges, A E
2009-09-28
ALE-AMR is a new hydrocode that we are developing as a predictive modeling tool for debris and shrapnel formation in high-energy laser experiments. In this paper we present our approach to implementing laser ray-tracing in ALE-AMR. We present the equations of laser ray tracing, our approach to efficient traversal of the adaptive mesh hierarchy in which we propagate computational rays through a virtual composite mesh consisting of the finest resolution representation of the modeled space, and anticipate simulations that will be compared to experiments for code validation.
Comparison of a 3-D GPU-Assisted Maxwell Code and Ray Tracing for Reflectometry on ITER
NASA Astrophysics Data System (ADS)
Gady, Sarah; Kubota, Shigeyuki; Johnson, Irena
2015-11-01
Electromagnetic wave propagation and scattering in magnetized plasmas are important diagnostics for high temperature plasmas. 1-D and 2-D full-wave codes are standard tools for measurements of the electron density profile and fluctuations; however, ray tracing results have shown that beam propagation in tokamak plasmas is inherently a 3-D problem. The GPU-Assisted Maxwell Code utilizes the FDTD (Finite-Difference Time-Domain) method for solving the Maxwell equations with the cold plasma approximation in a 3-D geometry. Parallel processing with GPGPU (General-Purpose computing on Graphics Processing Units) is used to accelerate the computation. Previously, we reported on initial comparisons of the code results to 1-D numerical and analytical solutions, where the size of the computational grid was limited by the on-board memory of the GPU. In the current study, this limitation is overcome by using domain decomposition and an additional GPU. As a practical application, this code is used to study the current design of the ITER Low Field Side Reflectometer (LSFR) for the Equatorial Port Plug 11 (EPP11). A detailed examination of Gaussian beam propagation in the ITER edge plasma will be presented, as well as comparisons with ray tracing. This work was made possible by funding from the Department of Energy for the Summer Undergraduate Laboratory Internship (SULI) program. This work is supported by the US DOE Contract No.DE-AC02-09CH11466 and DE-FG02-99-ER54527.
SolTrace: A Ray-Tracing Code for Complex Solar Optical Systems
Wendelin, Tim; Dobos, Aron; Lewandowski, Allan
2013-10-01
SolTrace is an optical simulation tool designed to model optical systems used in concentrating solar power (CSP) applications. The code was first written in early 2003, but has seen significant modifications and changes since its inception, including conversion from a Pascal-based software development platform to C++. SolTrace is unique in that it can model virtually any optical system utilizingthe sun as the source. It has been made available for free and as such is in use worldwide by industry, universities, and research laboratories. The fundamental design of the code is discussed, including enhancements and improvements over the earlier version. Comparisons are made with other optical modeling tools, both non-commercial and commercial in nature. Finally, modeled results are shownfor some typical CSP systems and, in one case, compared to measured optical data.
Can We Trace "Arbitrary" Rays to Locate an Image Formed by a Thin Lens?
ERIC Educational Resources Information Center
Suppapittayaporn, Decha; Panijpan, Bhinyo; Emarat, Narumon
2010-01-01
After learning how to trace the principal rays [Fig. 1(i)] through a thin lens in order to form the image in the conventional way, students sometimes ask whether it is possible to use other rays emanating from the object to form exactly the same image--for example, the two arbitrary rays shown in Fig. 1(ii). The answer is a definite yes, and this…
Integrating 2-D position sensitive X-ray detectors with low-density alkali halide storage targets
NASA Astrophysics Data System (ADS)
Haubold, H.-G.; Hoheisel, W.; Hiller, P.
1986-05-01
For the use in scattering experiments with synchrotron radiation, integrating position sensitive X-ray detectors are discussed. These detectors store the photon number equivalent charge (PNEC) in low-density alkali halide targets. Performance tests are given for a detector which uses a Gd 2O 2S fluorescence screen for X-ray detection and the low-density KCl storage target of a television SEC vidicon tube for photon integration. Rather than directly by X-rays, this target is charged by 6 keV electrons from the image intensifier section of the vidicon. Its excellent storage capability allows measurements of extremely high-contrast, high-flux X-ray patterns with the same accuracy as achieved with any single photon detection system if the discussed readout techniques are applied.
Yuk, Sunwoo; Park, Shin-Woong; Yi, Yun
2006-01-01
A 2D array radiation sensor, consisting of an array of PIN photodiodes bump bonded to readout integrated circuit (IC), has been developed for operation with low energy X-rays. The PIN photodiode array and readout IC for this system have been fabricated. The main performance measurements are the following: a few pA-scale leakage current, 350 pF junction capacitance, 30 microm-depth depletion layer and a 250 microm intrinsic layer at zero bias. This PIN photodiode array and readout IC were fabricated using a PIN photodiode process and standard 0.35 microm CMOS technology, respectively. The readout circuit is operated from a 3.3 V single power supply. Finally, a 2D array radiation sensor has been developed using bump bonding between the PIN photodiode and the readout electronics. PMID:17946079
NASA Astrophysics Data System (ADS)
Alizadeh, Mohamad Mahdi; Schuh, Harald; Schmidt, Michael
2015-06-01
For space geodetic techniques, operating in microwave band, ionosphere is a dispersive medium; thus, signals traveling through this medium are in the first approximation, affected proportional to the inverse of the square of their frequencies. This effect allows gaining information about the parameters of the ionosphere in terms of total electron content (TEC) or the electron density (Ne). Making use of this phenomenon, space geodetic techniques have turned into a capable tool for studying the ionosphere in the last decades. Up to now, two-dimensional (2-D) models of Vertical TEC (VTEC) have been widely developed and used by different communities; however, due to the fact that these models provide information about the integral of the whole electron content along the vertical or slant raypath, these maps are not useful when information about the ionosphere at different altitude is required. This paper presents a recent study which aims at developing a global 3-D model of the electron density, using measurements from Global Navigation Satellite Systems and by applying the ray tracing technique to the upper atmosphere. The developed modeling approach represents the horizontal variations of the electron density, with two sets of spherical harmonic expansions of degree and order 15. The height dependency of the electron density is represented by a multilayered Chapman profile function for the bottomside and topside ionosphere, and an appropriate model for the plasmasphere. In addition to the geodetic applications of the developed models, within this study, the 3-D models of electron density can include geophysical parameters like maximum electron density and its corresponding height. High-resolution modeling of these parameters allows an improved geophysical interpretation, which is essential in all studies of the upper atmosphere, space weather, and for the solar-terrestrial environment.
Ray-tracing method for creeping waves on arbitrarily shaped nonuniform rational B-splines surfaces.
Chen, Xi; He, Si-Yuan; Yu, Ding-Feng; Yin, Hong-Cheng; Hu, Wei-Dong; Zhu, Guo-Qiang
2013-04-01
An accurate creeping ray-tracing algorithm is presented in this paper to determine the tracks of creeping waves (or creeping rays) on arbitrarily shaped free-form parametric surfaces [nonuniform rational B-splines (NURBS) surfaces]. The main challenge in calculating the surface diffracted fields on NURBS surfaces is due to the difficulty in determining the geodesic paths along which the creeping rays propagate. On one single parametric surface patch, the geodesic paths need to be computed by solving the geodesic equations numerically. Furthermore, realistic objects are generally modeled as the union of several connected NURBS patches. Due to the discontinuity of the parameter between the patches, it is more complicated to compute geodesic paths on several connected patches than on one single patch. Thus, a creeping ray-tracing algorithm is presented in this paper to compute the geodesic paths of creeping rays on the complex objects that are modeled as the combination of several NURBS surface patches. In the algorithm, the creeping ray tracing on each surface patch is performed by solving the geodesic equations with a Runge-Kutta method. When the creeping ray propagates from one patch to another, a transition method is developed to handle the transition of the creeping ray tracing across the border between the patches. This creeping ray-tracing algorithm can meet practical requirements because it can be applied to the objects with complex shapes. The algorithm can also extend the applicability of NURBS for electromagnetic and optical applications. The validity and usefulness of the algorithm can be verified from the numerical results. PMID:23595326
An hybrid detector GEM-ASIC for 2-D soft X-ray imaging for laser produced plasma and pulsed sources
NASA Astrophysics Data System (ADS)
Pacella, D.; Claps, G.; De Angelis, R.; Murtas, F.
2016-03-01
The following paper presents a new 2-D detector (`GEMpix') in the soft X-ray range, having a wide dynamic range thanks to its intrisic gain, working in charge integration mode to be used for diagnosing laser produced plasma (LPP) or X-ray pulsed sources. It is a gas detector based on the Gas Electron Multiplier (GEM) technology with a quad-medipix chip as read-out electronics. In our prototype, the substitution of semiconductor material with a gas triple-GEM allows several advantages with respect to the detectors commonly used in LPP, as X-ray CCDs and Micro Channel Plates or Image Plates. In these experiments the configuration Time-over-Threshold (ToT) has been used, to measure the total charge released to the gas and collected by each pixel, integrated over the X-ray burst duration. Intensity response and spatial resolution has been measured first in laboratory for calibration, as function of the voltage applied to the GEMs, in single photon regime with energies between 3.7 and 17 keV. Subsequently it has been tested at the ABC laser facility (ENEA, Frascati). In this case, we measured the X-rays produced when the ABC neodymium laser, with pulse of 50 J and 3 ns time width, hits plane targets of aluminum. 2-D images have been acquired by means of a pinhole configuration with magnification 1.5 and 50 μ m of spatial resolution. The results are encouraging regarding the capability of this imaging detector to work in experiments where soft X-ray emissivity varies over many orders of magnitude.
Lee, Yong Geun; Lee, Jeongjin; Shin, Yeong-Gil; Kang, Ho Chul
2016-01-01
Enhancing 2D angiography while maintaining a low radiation dose has become an important research topic. However, it is difficult to enhance images while preserving vessel-structure details because X-ray noise and contrast blood vessels in 2D angiography have similar intensity distributions, which can lead to ambiguous images of vessel structures. In this paper, we propose a novel and fast vessel-enhancement method for 2D angiography. We apply filtering in the principal component analysis domain for vessel regions and background regions separately, using assumptions based on energy compaction. First, we identify an approximate vessel region using a Hessian-based method. Vessel and non-vessel regions are then represented sparsely by calculating their optimal bases separately. This is achieved by identifying periodic motion in the vessel region caused by the flow of the contrast medium through the blood vessels when viewed on the time axis. Finally, we obtain noise-free images by removing noise in the new coordinate domain for the optimal bases. Our method was validated for an X-ray system, using 10 low-dose sets for training and 20 low-dose sets for testing. The results were compared with those for a high-dose dataset with respect to noise-free images. The average enhancement rate was 93.11±0.71%. The average processing time for enhancing video comprising 50-70 frames was 0.80±0.35s, which is much faster than the previously proposed technique. Our method is applicable to 2D angiography procedures such as catheterization, which requires rapid and natural vessel enhancement. PMID:26483302
Optimizing detector geometry for trace element mapping by X-ray fluorescence
Sun, Yue; Gleber, Sophie -Charlotte; Jacobsen, Chris; Kirz, Janos; Vogt, Stefan
2015-01-01
We report that trace metals play critical roles in a variety of systems, ranging from cells to photovoltaics. X-Ray Fluorescence (XRF) microscopy using X-ray excitation provides one of the highest sensitivities available for imaging the distribution of trace metals at sub-100 nm resolution. With the growing availability and increasing performance of synchrotron light source based instruments and X-ray nanofocusing optics, and with improvements in energy-dispersive XRF detectors, what are the factors that limit trace element detectability? To address this question, we describe an analytical model for the total signal incident on XRF detectors with various geometries, including the spectral response of energy dispersive detectors. This model agrees well with experimentally recorded X-ray fluorescence spectra, and involves much shorter calculation times than with Monte Carlo simulations. With such a model, one can estimate the signal when a trace element is illuminated with an X-ray beam, and when just the surrounding non-fluorescent material is illuminated. From this signal difference, a contrast parameter can be calculated and this can in turn be used to calculate the signal-to-noise ratio (S/N) for detecting a certain elemental concentration. We apply this model to the detection of trace amounts of zinc in biological materials, and to the detection of small quantities of arsenic in semiconductors. In conclusion, we conclude that increased detector collection solid angle is (nearly) always advantageous even when considering the scattered signal. However, given the choice between a smaller detector at 90° to the beam versus a larger detector at 180° (in a backscatter-like geometry), the 90° detector is better for trace element detection in thick samples, while the larger detector in 180° geometry is better suited to trace element detection in thin samples.
Optimizing detector geometry for trace element mapping by X-ray fluorescence
Sun, Yue; Gleber, Sophie -Charlotte; Jacobsen, Chris; Kirz, Janos; Vogt, Stefan
2015-01-01
We report that trace metals play critical roles in a variety of systems, ranging from cells to photovoltaics. X-Ray Fluorescence (XRF) microscopy using X-ray excitation provides one of the highest sensitivities available for imaging the distribution of trace metals at sub-100 nm resolution. With the growing availability and increasing performance of synchrotron light source based instruments and X-ray nanofocusing optics, and with improvements in energy-dispersive XRF detectors, what are the factors that limit trace element detectability? To address this question, we describe an analytical model for the total signal incident on XRF detectors with various geometries, including the spectral responsemore » of energy dispersive detectors. This model agrees well with experimentally recorded X-ray fluorescence spectra, and involves much shorter calculation times than with Monte Carlo simulations. With such a model, one can estimate the signal when a trace element is illuminated with an X-ray beam, and when just the surrounding non-fluorescent material is illuminated. From this signal difference, a contrast parameter can be calculated and this can in turn be used to calculate the signal-to-noise ratio (S/N) for detecting a certain elemental concentration. We apply this model to the detection of trace amounts of zinc in biological materials, and to the detection of small quantities of arsenic in semiconductors. In conclusion, we conclude that increased detector collection solid angle is (nearly) always advantageous even when considering the scattered signal. However, given the choice between a smaller detector at 90° to the beam versus a larger detector at 180° (in a backscatter-like geometry), the 90° detector is better for trace element detection in thick samples, while the larger detector in 180° geometry is better suited to trace element detection in thin samples.« less
Optimizing detector geometry for trace element mapping by X-ray fluorescence
Sun, Yue; Gleber, Sophie-Charlotte; Jacobsen, Chris; Kirz, Janos; Vogt, Stefan
2016-01-01
Trace metals play critical roles in a variety of systems, ranging from cells to photovoltaics. X-Ray Fluorescence (XRF) microscopy using X-ray excitation provides one of the highest sensitivities available for imaging the distribution of trace metals at sub-100 nm resolution. With the growing availability and increasing performance of synchrotron light source based instruments and X-ray nanofocusing optics, and with improvements in energy-dispersive XRF detectors, what are the factors that limit trace element detectability? To address this question, we describe an analytical model for the total signal incident on XRF detectors with various geometries, including the spectral response of energy dispersive detectors. This model agrees well with experimentally recorded X-ray fluorescence spectra, and involves much shorter calculation times than with Monte Carlo simulations. With such a model, one can estimate the signal when a trace element is illuminated with an X-ray beam, and when just the surrounding non-fluorescent material is illuminated. From this signal difference, a contrast parameter can be calculated and this can in turn be used to calculate the signal-to-noise ratio (S/N) for detecting a certain elemental concentration. We apply this model to the detection of trace amounts of zinc in biological materials, and to the detection of small quantities of arsenic in semiconductors. We conclude that increased detector collection solid angle is (nearly) always advantageous even when considering the scattered signal. However, given the choice between a smaller detector at 90° to the beam versus a larger detector at 180° (in a backscatter-like geometry), the 90° detector is better for trace element detection in thick samples, while the larger detector in 180° geometry is better suited to trace element detection in thin samples. PMID:25600825
A data distributed parallel algorithm for ray-traced volume rendering
NASA Technical Reports Server (NTRS)
Ma, Kwan-Liu; Painter, James S.; Hansen, Charles D.; Krogh, Michael F.
1993-01-01
This paper presents a divide-and-conquer ray-traced volume rendering algorithm and a parallel image compositing method, along with their implementation and performance on the Connection Machine CM-5, and networked workstations. This algorithm distributes both the data and the computations to individual processing units to achieve fast, high-quality rendering of high-resolution data. The volume data, once distributed, is left intact. The processing nodes perform local ray tracing of their subvolume concurrently. No communication between processing units is needed during this locally ray-tracing process. A subimage is generated by each processing unit and the final image is obtained by compositing subimages in the proper order, which can be determined a priori. Test results on both the CM-5 and a group of networked workstations demonstrate the practicality of our rendering algorithm and compositing method.
Long gamma-ray bursts trace the star formation history
Dado, Shlomo; Dar, Arnon
2014-04-10
We show that if the broad-line supernova explosions of Type Ic (SNeIc) produce the bulk of the observed long duration gamma-ray bursts (LGRBs), including high- and low-luminosity LGRBs and X-ray flashes, and if the LGRBs have the geometry assumed in the cannonball model of LGRBs, then their rate, measured by Swift, and their redshift distribution are consistent with the star formation rate (SFR) over the entire range of redshifts where the SFR has been measured with sufficient accuracy.
Ray Tracing Through Non-Rotationally Symmetrical Systems With A Desktop Computer
NASA Astrophysics Data System (ADS)
Mackay, R. M.; Busse lle, F. J.
1986-10-01
A general ray-trace program has been developed for use on a desktop computer which traces finite rays through any non-rotationally symmetrical system. In particular any combination of decentred, tilted and rotated surface has been considered. Surface types such as Conic sections with and without Aspherics, Toric surfaces, surfaces of S and T Cylindrical sections, and Axicons, may be ray-traced. Each surface is defined in terms of a local rectangular co-ordinate system and has a particular aperture shape attributed to it. Aperture shapes may be defined as circular, elliptical, rectangular or quadrilateral. Also the centre of any aperture shape may be displaced from its local coordinate origin to facilitate the tracing of off-axis paraboloids. Before transferring to the next surface, the local coordinates are referred back to an initial reference coordinate system. Finally a means of assessing aberrations has been included. The main task here was to get a mathematical model of a non-rotationally symmetrical finite ray-trace running on an inexpensive desk top computer. The program was written for the BBC MICRO in order to investigate devices such as scanning systems for modern Thermal Imagers etc.
NASA Astrophysics Data System (ADS)
Pouyet, E.; Cotte, M.; Fayard, B.; Salomé, M.; Meirer, F.; Mehta, A.; Uffelman, E. S.; Hull, A.; Vanmeert, F.; Kieffer, J.; Burghammer, M.; Janssens, K.; Sette, F.; Mass, J.
2015-11-01
The chemical and physical alterations of cadmium yellow (CdS) paints in Henri Matisse's The Joy of Life (1905-1906, The Barnes Foundation) have been recognized since 2006, when a survey by portable X-ray fluorescence identified this pigment in all altered regions of the monumental painting. This alteration is visible as fading, discoloration, chalking, flaking, and spalling of several regions of light to medium yellow paint. Since that time, synchrotron radiation-based techniques including elemental and spectroscopic imaging, as well as X-ray scattering have been employed to locate and identify the alteration products observed in this and related works by Henri Matisse. This information is necessary to formulate one or multiple mechanisms for degradation of Matisse's paints from this period, and thus ensure proper environmental conditions for the storage and the display of his works. This paper focuses on 2D full-field X-ray Near Edge Structure imaging, 2D micro-X-ray Diffraction, X-ray Fluorescence, and Fourier Transform Infra-red imaging of the altered paint layers to address one of the long-standing questions about cadmium yellow alteration—the roles of cadmium carbonates and cadmium sulphates found in the altered paint layers. These compounds have often been assumed to be photo-oxidation products, but could also be residual starting reagents from an indirect wet process synthesis of CdS. The data presented here allow identifying and mapping the location of cadmium carbonates, cadmium chlorides, cadmium oxalates, cadmium sulphates, and cadmium sulphides in thin sections of altered cadmium yellow paints from The Joy of Life and Matisse's Flower Piece (1906, The Barnes Foundation). Distribution of various cadmium compounds confirms that cadmium carbonates and sulphates are photo-degradation products in The Joy of Life, whereas in Flower Piece, cadmium carbonates appear to have been a [(partially) unreacted] starting reagent for the yellow paint, a role
NASA Astrophysics Data System (ADS)
Ortega, I.; Coburn, S.; Kassianov, E.; Barnard, J.; Berg, L. K.; Hostetler, C. A.; Hair, J. W.; Ferrare, R. A.; Volkamer, R. M.
2012-12-01
The two Column Aerosol Project (TCAP) investigates uncertainties in the aerosol direct effect in the northern hemisphere mid-latitudes. The DOE Atmospheric Radiation Measurement (ARM) Mobile Facility (AMF) and Mobile Aerosol Observing System (MAOS) provide an opportunity for 1) atmospheric radiation closure studies, and 2) test retrievals of aerosol optical properties in the presence and absence of clouds. This presentation discusses innovative means to access column information about aerosol optical properties in the lower atmosphere from ground based measurements of solar stray light spectra in the hyperspectral domain, i.e., measurements of the Raman Scattering Probability (RSP, the probability that an observed photon has undergone a rotational Raman scattering event), and oxygen dimer slant column densities (O4 SCD) by means of the University of Colorado 2D scanning ground Multi AXis Differential Optical Absorption Spectroscopy (2D-GMAX-DOAS) instrument that was located at the ARM/MAOS site at Cape Cod, MA. We compare retrievals of aerosol optical properties with those retrieved from the MFRSR and the Cimel Sunphotometer, for case studies in the presence/absence of clouds, and assess the need for atmospheric correction of NO2. 2D-GMAX-DOAS also facilitates a link between the ground-based ARM/MAOS dataset and DoE's G1 aircraft, NASA's King Air aircraft, and NASA's OMI satellite (i.e., NO2 vertical column). Early results that explore these linkages are presented for a case study that combines ground based MFRSR, in-situ observations aboard the G1 aircraft, as well as High Spectral Resolution LIDAR aboard the King Air aircraft.
Using Divergent Δ12CH2D2 and Δ13CH3D to Trace the Provenance and Evolution of Methane Gas
NASA Astrophysics Data System (ADS)
Young, E. D.; Freedman, P.; Mills, M.; Rumble, D.
2015-12-01
Measurements of Δ13CH3D (deviations in Δ13CH3D/12CH4 from stochastic; Ono et al. Anal. Chem. v.86, p.6487, 2014) or Δ18 (from (12CH2D2 + 13CH3D)/12CH4; Stolper et al. Science, v.344, p.1500, 2014, ) have been used to infer temperatures of formation of methane gas. However, departures from thermodynamic equilibrium isotopic bond ordering will result from any fractionating process that do not include bond rupture and reformation, including mixing, diffusion, and kinetic processing. This is because the isotopic bond ordering no longer reflects the bulk isotopic composition once fractionation occurs. A direct measure of departures from thermodynamic equilibrium isotopic bond ordering in methane comes from both Δ12CH2D2 and Δ13CH3D in the same gas. Until now, this has not been possible due to instrumental limitations. We have carried out measurements of Δ12CH2D2 and Δ13CH3D in methane gas mixtures using a unique, large-geometry double-focusing isotope ratio mass spectrometer (IRMS), the Panorama, in order to investigate the usefulness of these two mass-18 isotopologues as tracers of mixing of methane sources. This instrument has a dispersion/magnification ratio, the parameter of merit for mass resolving power, of ~ 1400 mm that exceeds that of any other gas-source IRMS by more than 3.5x and is slightly larger than that for large-geometry SIMS instruments. With this geometry we routinely operate with mass resolving power (M/ΔM, 5% and 95%) of 40,000 or greater with useful sensitivity for isotope ratio analysis. For these experiments we mixed two gases with bulk D/H differing by 100 ‰. The results follow theoretical expectations within uncertainties of 0.5 ‰ for Δ12CH2D2 and 0.1 ‰ for Δ13CH3D. Precision is sufficient to detect as little as 10% mixing in this system. This precision would also be capable of detecting subtle departures from equilibrium caused by diffusion and kinetic bond rupture (e.g. CH4 + OH).
Ray-tracing in a two-dimensional ionosphere
NASA Astrophysics Data System (ADS)
Labahn, R. W.
1985-08-01
The quasi-parabolic method is adapted to ray-track through a medium with horizontal gradients in electron density. The resulting method is applicable to any model ionosphere without the requirements for numerical derivatives. Example calculations are given for model ledges and troughs and an average worldwide ionospheric model.
Time resolved, 2-D hard X-ray imaging of relativistic electron-beam target interactions on ETA-II
Crist, C.E.; Sampayan, S.; Westenskow, G.; Caporaso, G.; Houck, T.; Weir, J.; Trimble, D.; Krogh, M.
1998-11-01
Advanced radiographic applications require a constant source size less than 1 mm. To study the time history of a relativistic electron beam as it interacts with a bremsstrahlung converter, one of the diagnostics they use is a multi-frame time-resolved hard x-ray camera. They are performing experiments on the ETA-II accelerator at Lawrence Livermore National Laboratory to investigate details of the electron beam/converter interactions. The camera they are using contains 6 time-resolved images, each image is a 5 ns frame. By starting each successive frame 10 ns after the previous frame, they create a 6-frame movie from the hard x-rays produced from the interaction of the 50-ns electron beam pulse.
NASA Astrophysics Data System (ADS)
Dinten, Jean M.; Robert-Coutant, Christine; Darboux, Michel; Gonon, Georges; Bordy, Thomas
2003-06-01
In a previous paper (SPIE Medical Imaging 2001), a dual energy method for bone densitometry using a 2D digital radiographic detector has been presented. In this paper, calcium content quantification performance of the approach is precised. The main challenge is to achieve quantification using scatter-corrected dual energy acquisitions. Therefore a scatter estimation approach, based on an expression of scatter as a functional of the primary flux, has been developed. This expression is derived from the Klein and Nishina equation and includes tabulated scatter level values. The calcium quantification performances are validated on two configurations. A first one is issued from criteria developed by the French "Groupe de Recherche et d'Information sur les Osteoporoses." It is based on the use of a phantom made of five 3mm thick PVC sheets in the form of five steps, representing five different bone mineral density values, included in a lucite container filled with water. Additional lucite plates can be put over the phantom. This phantom has been used for evaluation of quantification robustness versus patient thickness and composition variations, and for accuracy evaluation. The second configuration, composed of small calcified objects (representative of lung nodules), is used for evaluating capacities to differentiate calcified from non calcified nodules and to test calcium content quantification performance.
NASA Astrophysics Data System (ADS)
Ichikawa, Ryuichi; Hobiger, Thomas; Shoji, Yoshinori; Koyama, Yasuhiro; Kondo, Tesuro
2010-05-01
We have been developing a state-of-art tool to estimate the atmospheric path delays by ray-tracing through meso-scale analysis (MANAL data) data with 10km grid interval, which is operationally used for numerical weather prediction by Japan Meteorological Agency (JMA). The tools, which we have named 'KAshima RAytracing Tools (KARAT)', are capable of calculating total slant delays and ray-bending angles considering real atmospheric phenomena. The KARAT can estimate atmospheric slant delays by three different calculation schemes. These are (1) a piece-wise linear propagation, (2) an analytical 2-D ray-propagation model by Thayer, and (3) a 3-D Eikonal solver. By computing GPS PPP solutions for 57 GPS sites of the GEONET (GPS Earth Observation Network System) operated by Geographical Survey Institute (GSI) of Japan it could be shown that KARAT performs slightly better than results based on the Global Mapping Function (GMF) and the Vienna Mapping Function 1 (VMF1), whereas for the latter two also linear gradient models had to be applied. The grid interval of the MANAL data was updated from 10km to 5km on April 7, 2009. In addition, on October 27, 2009 the JMA started data assimilation of zenith wet delay obtained by the GEONET for meso-scale numerical weather prediction. We are now evaluating impacts of data scheme improvements and assimilation strategy change on the slant delay reduction. We will include these preliminary results in our presentation.
NASA Astrophysics Data System (ADS)
Lee, Kyoung Jin
Understanding and modeling seismic wave propagation is important in regional and exploration seismology. Ray tracing is a powerful and popular method for this purpose. Wavefront construction (WFC) method handles wavefronts instead of individual rays, thereby controlling proper ray density on the wavefront. By adaptively controlling rays over a wavefront, it efficiently models wave propagation. Algorithms for a quasi-P wave wavefront construction method and a new coordinate system used to generate wavefront construction mesh are proposed and tested for numerical properties and modeling capabilities. Traveltimes, amplitudes, and other parameters, which can be used for seismic imaging such as migrations and synthetic seismograms, are computed from the wavefront construction method. Modeling with wavefront construction code is applied to anisotropic media as well as isotropic media. Synthetic seismograms are computed using the wavefront construction method as a new way of generating synthetics. To incorporate layered velocity models, the model based interpolation (MBI) ray tracing method, which is designed to take advantage of the wavefront construction method as well as conventional ray tracing methods, is proposed and experimental codes are developed for it. Many wavefront construction codes are limited to smoothed velocity models for handling complicated problems in layered velocity models and the conventional ray tracing methods suffer from the inability to control ray density during wave propagation. By interpolating the wavefront near model boundaries, it is possible to handle the layered velocity model as well as overcome ray density control problems in conventional methods. The test results revealed this new method can be an effective modeling tool for accurate and effective computing.
Magnetospheric Whistler Mode Ray Tracing with the Inclusion of Finite Electron and Ion Temperature
NASA Astrophysics Data System (ADS)
Maxworth, A. S.; Golkowski, M.
2015-12-01
Ray tracing is an important technique for the study of whistler mode wave propagation in the Earth's magnetosphere. In numerical ray tracing the trajectory of a wave packet is calculated at each point in space by solving the Haselgrove equations, assuming a smooth, loss-less medium with no mode coupling. Previous work on ray tracing has assumed a cold plasma environment with negligible electron and ion temperatures. In this work we present magnetospheric whistler mode wave ray tracing results with the inclusion of finite ion and electron temperature. The inclusion of finite temperature effects makes the fourth order dispersion relation become sixth order. We compare our results with the work done by previous researchers for cold plasma environments, using two near earth space models (NGO and GCPM). Inclusion of finite temperature closes the otherwise open refractive index surface near the lower hybrid resonance frequency and affects the magnetospheric reflection of whistler waves. We also asses the main changes in the ray trajectory and implications for cyclotron resonance wave particle interactions including energetic particle precipitation.
Automatic creation of object hierarchies for ray tracing
NASA Technical Reports Server (NTRS)
Goldsmith, Jeffrey; Salmon, John
1987-01-01
Various methods for evaluating generated trees are proposed. The use of the hierarchical extent method of Rubin and Whitted (1980) to find the objects that will be hit by a ray is examined. This method employs tree searching; the construction of a tree of bounding volumes in order to determine the number of objects that will be hit by a ray is discussed. A tree generation algorithm, which uses a heuristic tree search strategy, is described. The effects of shuffling and sorting on the input data are investigated. The cost of inserting an object into the hierarchy during the construction of a tree algorithm is estimated. The steps involved in estimating the number of intersection calculations are presented.
A boundary integral formalism for stochastic ray tracing in billiards
Chappell, David J.; Tanner, Gregor
2014-12-15
Determining the flow of rays or non-interacting particles driven by a force or velocity field is fundamental to modelling many physical processes. These include particle flows arising in fluid mechanics and ray flows arising in the geometrical optics limit of linear wave equations. In many practical applications, the driving field is not known exactly and the dynamics are determined only up to a degree of uncertainty. This paper presents a boundary integral framework for propagating flows including uncertainties, which is shown to systematically interpolate between a deterministic and a completely random description of the trajectory propagation. A simple but efficient discretisation approach is applied to model uncertain billiard dynamics in an integrable rectangular domain.
PixFEL: developing a fine pitch, fast 2D X-ray imager for the next generation X-FELs
NASA Astrophysics Data System (ADS)
Ratti, L.; Comotti, D.; Fabris, L.; Grassi, M.; Lodola, L.; Malcovati, P.; Manghisoni, M.; Re, V.; Traversi, G.; Vacchi, C.; Bettarini, S.; Casarosa, G.; Forti, F.; Morsani, F.; Paladino, A.; Paoloni, E.; Rizzo, G.; Benkechkache, M. A.; Dalla Betta, G.-F.; Mendicino, R.; Pancheri, L.; Verzellesi, G.; Xu, H.
2015-10-01
The PixFEL project is conceived as the first stage of a long term research program aiming at the development of advanced X-ray imaging instrumentation for applications at the free electron laser (FEL) facilities. The project aims at substantially advancing the state-of-the-art in the field of 2D X-ray imaging by exploring cutting-edge solutions for sensor development, for integration processes and for readout channel architectures. The main focus is on the development of the fundamental microelectronic building blocks for detector readout and on the technologies for the assembly of a multilayer module with minimum dead area. This work serves the purpose of introducing the main features of the project, together with the simulation results leading to the first prototyping run.
NASA Astrophysics Data System (ADS)
General, S.; Pöhler, D.; Sihler, H.; Bobrowski, N.; Frieß, U.; Zielcke, J.; Horbanski, M.; Shepson, P. B.; Stirm, B. H.; Simpson, W. R.; Weber, K.; Fischer, C.; Platt, U.
2014-03-01
Many relevant processes in tropospheric chemistry take place on rather small scales (e.g. tens to hundreds of meters) but often influence areas of several square kilometer. Thus, measurements of the involved trace gases with high spatial resolution are of great scientific interest. In order to identify individual sources and sinks and ultimately to improve chemical transport models, we developed a new airborne instrument, which is based on the well established DOAS method. The Heidelberg Airborne Imaging Differential Optical Absorption Spectrometer Instrument (HAIDI) is a passive imaging DOAS spectrometer, which is capable of recording horizontal and vertical trace gas distributions with a resolution of better than 100 m. Observable species include NO2, HCHO, C2H2O2, H2O, O3, O4, SO2, IO, OClO and BrO. Here we report a technical description of the instrument including its custom build spectrographs and CCD detectors. Also first results from measurements with the new instrument are presented. These comprise spatial resolved SO2 and BrO in volcanic plumes, mapped at Mt. Etna (Sicily, Italy), NO2 emissions in the metropolitan area of Indianapolis (Indiana, USA) as well as BrO and NO2 distributions measured during arctic springtime in context of the BROMEX campaign, which was performed 2012 in Barrow (Alaska, USA).
NASA Astrophysics Data System (ADS)
General, S.; Pöhler, D.; Sihler, H.; Bobrowski, N.; Frieß, U.; Zielcke, J.; Horbanski, M.; Shepson, P. B.; Stirm, B. H.; Simpson, W. R.; Weber, K.; Fischer, C.; Platt, U.
2014-10-01
Many relevant processes in tropospheric chemistry take place on rather small scales (e.g., tens to hundreds of meters) but often influence areas of several square kilometer. Thus, measurements of the involved trace gases with high spatial resolution are of great scientific interest. In order to identify individual sources and sinks and ultimately to improve chemical transport models, we developed a new airborne instrument, which is based on the well established Differential Optical Absorption Spectroscopy (DOAS) method. The Heidelberg Airborne Imaging DOAS Instrument (HAIDI) is a passive imaging DOAS spectrometer, which is capable of recording horizontal and vertical trace gas distributions with a resolution of better than 100 m. Observable species include NO2, HCHO, C2H2O2, H2O, O3, O4, SO2, IO, OClO and BrO. Here we give a technical description of the instrument including its custom-built spectrographs and CCD detectors. Also first results from measurements with the new instrument are presented. These comprise spatial resolved SO2 and BrO in volcanic plumes, mapped at Mt. Etna (Sicily, Italy), NO2 emissions in the metropolitan area of Indianapolis (Indiana, USA) as well as BrO and NO2 distributions measured during arctic springtime in context of the BRomine, Ozone, and Mercury EXperiment (BROMEX) campaign, which was performed 2012 in Barrow (Alaska, USA).
NASA Astrophysics Data System (ADS)
Ortega, I.; Coburn, S.; Kassianov, E.; Barnard, J.; Berg, L. K.; Hostetler, C. A.; Hair, J. W.; Ferrare, R. A.; Hodges, G.; Lantz, K. O.; Volkamer, R.
2013-12-01
The two Column Aerosol Project (TCAP) investigates uncertainties in the aerosol direct effect in the northern hemisphere mid-latitudes. The DOE Atmospheric Radiation Measurement (ARM) Mobile Facility (AMF) and Mobile Aerosol Observing System (MAOS) provide an opportunity for 1) atmospheric radiation closure studies, and 2) test retrievals of aerosol optical and microphysical properties in the presence and absence of clouds. The University of Colorado 2D-MAX-DOAS instrument was deployed during the first intensive period of the TCAP field project. This presentation presents an inversion algorithm to obtain both aerosol extinction profiles and column from off axis measurements, and infer - for the first time- aerosol microphysical properties from the solar angular distribution of sky radiances. An innovative aspect of the method is that the optical and microphysical properties of aerosols are inferred without the need for an absolute radiance calibration. We compare retrievals of aerosol optical properties with those retrieved from the MFRSR and the Cimel Sunphotometer, for case studies in the presence/absence of clouds, and assess the need for atmospheric correction of NO2. 2D-GMAX-DOAS also facilitates a link between the ground-based ARM/MAOS dataset and DoE's G1 aircraft, NASA's King Air aircraft, and NASA's OMI satellite (i.e., NO2 vertical column). Early results that explore these linkages are presented for a case study that combines ground based MFRSR, in-situ observations aboard the G1 aircraft, as well as High Spectral Resolution LIDAR aboard the King Air aircraft.
Odyssey: Ray tracing and radiative transfer in Kerr spacetime
NASA Astrophysics Data System (ADS)
Pu, Hung-Yi; Yun, Kiyun; Younsi, Ziri; Yoon, Suk-Jin
2016-01-01
Odyssey is a GPU-based General Relativistic Radiative Transfer (GRRT) code for computing images and/or spectra in Kerr metric describing the spacetime around a rotating black hole. Odyssey is implemented in CUDA C/C++. For flexibility, the namespace structure in C++ is used for different tasks; the two default tasks presented in the source code are the redshift of a Keplerian disk and the image of a Keplerian rotating shell at 340GHz. Odyssey_Edu, an educational software package for visualizing the ray trajectories in the Kerr spacetime that uses Odyssey, is also available.
Modeling and Measurement of 3D Deformation of Scoliotic Spine Using 2D X-ray Images
NASA Astrophysics Data System (ADS)
Li, Hao; Leow, Wee Kheng; Huang, Chao-Hui; Howe, Tet Sen
Scoliosis causes deformations such as twisting and lateral bending of the spine. To correct scoliotic deformation, the extents of 3D spinal deformation need to be measured. This paper studies the modeling and measurement of scoliotic spine based on 3D curve model. Through modeling the spine as a 3D Cosserat rod, the 3D structure of a scoliotic spine can be recovered by obtaining the minimum potential energy registration of the rod to the scoliotic spine in the x-ray image. Test results show that it is possible to obtain accurate 3D reconstruction using only the landmarks in a single view, provided that appropriate boundary conditions and elastic properties are included as constraints.
An Energy Conservative Ray-Tracing Method With a Time Interpolation of the Force Field
Yao, Jin
2015-02-10
A new algorithm that constructs a continuous force field interpolated in time is proposed for resolving existing difficulties in numerical methods for ray-tracing. This new method has improved accuracy, but with the same degree of algebraic complexity compared to Kaisers method.
Comparing TID simulations using 3-D ray tracing and mirror reflection
NASA Astrophysics Data System (ADS)
Huang, X.; Reinisch, B. W.; Sales, G. S.; Paznukhov, V. V.; Galkin, I. A.
2016-04-01
Measuring the time variations of Doppler frequencies and angles of arrival (AoA) of ionospherically reflected HF waves has been proposed as a means of detecting the occurrence of traveling ionospheric disturbances (TIDs). Simulations are made using ray tracing through the International Reference Ionosphere (IRI) electron density model in an effort to reproduce measured signatures. The TID is represented by a wavelike perturbation of the 3-D electron density traveling horizontally in the ionosphere with an amplitude that varies sinusoidally with time. By judiciously selecting the TID parameters the ray tracing simulation reproduces the observed Doppler frequencies and AoAs. Ray tracing in a 3-D realistic ionosphere is, however, excessively time consuming considering the involved homing procedures. It is shown that a carefully selected reflecting corrugated mirror can reproduce the time variations of the AoA and Doppler frequency. The results from the ray tracing through the IRI model ionosphere and the mirror model reflections are compared to assess the applicability of the mirror-reflection model.
GPU-based ray tracing algorithm for high-speed propagation prediction in typical indoor environments
NASA Astrophysics Data System (ADS)
Guo, Lixin; Guan, Xiaowei; Liu, Zhongyu
2015-10-01
A fast 3-D ray tracing propagation prediction model based on virtual source tree is presented in this paper, whose theoretical foundations are geometrical optics(GO) and the uniform theory of diffraction(UTD). In terms of typical single room indoor scene, taking the geometrical and electromagnetic information into account, some acceleration techniques are adopted to raise the efficiency of the ray tracing algorithm. The simulation results indicate that the runtime of the ray tracing algorithm will sharply increase when the number of the objects in the single room is large enough. Therefore, GPU acceleration technology is used to solve that problem. As is known to all, GPU is good at calculation operation rather than logical judgment, so that tens of thousands of threads in CUDA programs are able to calculate at the same time, in order to achieve massively parallel acceleration. Finally, a typical single room with several objects is simulated by using the serial ray tracing algorithm and the parallel one respectively. It can be found easily from the results that compared with the serial algorithm, the GPU-based one can achieve greater efficiency.
NASA Astrophysics Data System (ADS)
Inoue, Koji; Saito, Haruo; Nagashima, Yasuyuki; Hyodo, Toshio; Nagai, Yasuyoshi; Muramatsu, Shinichi; Nagai, Shota; Masuda, Keisuke
2002-07-01
A new position-sensitive γ-ray detector to be used in a two-dimensional angular correlation of positron annihilation radiation (2D-ACAR) apparatus has been developed. It consists of 36 compact position-sensitive photomultiplier tubes (PS-PMT: HAMAMATSU R5900-00-C8), a light guide, and 2676 Bi 4Ge 3O 12 (BGO) scintillator pieces of size 2.6 mm×2.6 mm×18 mm. A high detection efficiency for 511 keV γ-ray is achieved with the length of BGO scintillators used. The detection area is about 160 mm×160 mm. The 288 anode outputs of the PS-PMTs are wired and connected to resistor chains from which 16 outputs (8 outputs each along the X and Y directions) are taken to identify the incident position of the γ-ray. The spatial resolution is about 3 mm (FWHM). The timing signal taken from the last dynodes of the PS-PMTs gives a timing resolution of 7.7 ns (FWHM) for 511 keV positron annihilation γ-rays.
Tracing Chromospheric Evaporation in Radio and Soft X-rays
NASA Technical Reports Server (NTRS)
Aschwanden, Markus J.
1997-01-01
There are three publications in refereed journals and several presentations at scientific conferences resulted from this work, over a period of 6 months during 1995/1996. In the first paper, the discovery of the chromospheric evaporation process at radio wavelengths is described. In the second paper, the radio detection is used to quantify electron densities in the upflowing heated plasma in flare loops, which is then compared with independent other density measurements from soft X-rays, or the plasma frequency of electron beams originating in the acceleration region. In the third paper, the diagnostic results of the chromospheric evaporation process are embedded into a broader picture of a standard flare scenario. Abstracts of these three papers are attached.
Leong, Andrew F. T.; Islam, M. Sirajul; Kitchen, Marcus J.; Fouras, Andreas; Wallace, Megan J.; Hooper, Stuart B.
2013-04-15
Purpose: Described herein is a new technique for measuring regional lung air volumes from two-dimensional propagation-based phase contrast x-ray (PBI) images at very high spatial and temporal resolution. Phase contrast dramatically increases lung visibility and the outlined volumetric reconstruction technique quantifies dynamic changes in respiratory function. These methods can be used for assessing pulmonary disease and injury and for optimizing mechanical ventilation techniques for preterm infants using animal models. Methods: The volumetric reconstruction combines the algorithms of temporal subtraction and single image phase retrieval (SIPR) to isolate the image of the lungs from the thoracic cage in order to measure regional lung air volumes. The SIPR algorithm was used to recover the change in projected thickness of the lungs on a pixel-by-pixel basis (pixel dimensions {approx}16.2 {mu}m). The technique has been validated using numerical simulation and compared results of measuring regional lung air volumes with and without the use of temporal subtraction for removing the thoracic cage. To test this approach, a series of PBI images of newborn rabbit pups mechanically ventilated at different frequencies was employed. Results: Regional lung air volumes measured from PBI images of newborn rabbit pups showed on average an improvement of at least 20% in 16% of pixels within the lungs in comparison to that measured without the use of temporal subtraction. The majority of pixels that showed an improvement was found to be in regions occupied by bone. Applying the volumetric technique to sequences of PBI images of newborn rabbit pups, it is shown that lung aeration at birth can be highly heterogeneous. Conclusions: This paper presents an image segmentation technique based on temporal subtraction that has successfully been used to isolate the lungs from PBI chest images, allowing the change in lung air volume to be measured over regions as small as the pixel size. Using
Irick, S.C.; Jung, C.R.
1997-07-01
There are a number of ray trace programs currently used for the design of synchrotron beamlines. While several of these programs have been written and used mostly within the programmer`s institution, many have also been available to the general public. This paper discusses three such programs. One is a commercial product oriented for the general optical designer (not specifically for synchrotron beamlines). One is designed for synchrotron beamlines and is free with restricted availability. Finally, one is designed for synchrotron beamlines and is used primarily in one institution. The wealth of information from general optical materials and components catalogs is readily available in the commercial program for general optical designs. This makes the design of an infrared beamline easier from the standpoint of component selection. However, this program is not easily configured for synchrotron beamline designs, particularly for a bending magnet source. The synchrotron ray trace programs offer a variety of sources, but generally are not as easy to use from the standpoint of the user interface. This paper shows ray traces of the same beamline Optikwerks, SHADOW, and RAY, and compares the results.
Ray tracing in a finite-element domain using nodal basis functions.
Schrader, Karl N; Subia, Samuel R; Myre, John W; Summers, Kenneth L
2014-08-20
A method is presented for tracing rays through a medium discretized as finite-element volumes. The ray-trajectory equations are cast into the local element coordinate frame, and the full finite-element interpolation is used to determine instantaneous index gradient for the ray-path integral equation. The finite-element methodology is also used to interpolate local surface deformations and the surface normal vector for computing the refraction angle when launching rays into the volume, and again when rays exit the medium. The procedure is applied to a finite-element model of an optic with a severe refractive-index gradient, and the results are compared to the closed-form gradient ray-path integral approach. PMID:25321137
Reverse Monte Carlo ray-tracing for radiative heat transfer in combustion systems
NASA Astrophysics Data System (ADS)
Sun, Xiaojing
Radiative heat transfer is a dominant heat transfer phenomenon in high temperature systems. With the rapid development of massive supercomputers, the Monte-Carlo ray tracing (MCRT) method starts to see its applications in combustion systems. This research is to find out if Monte-Carlo ray tracing can offer more accurate and efficient calculations than the discrete ordinates method (DOM). Monte-Carlo ray tracing method is a statistical method that traces the history of a bundle of rays. It is known as solving radiative heat transfer with almost no approximation. It can handle nonisotropic scattering and nongray gas mixtures with relative ease compared to conventional methods, such as DOM and spherical harmonics method, etc. There are two schemes in Monte-Carlo ray tracing method: forward and backward/reverse. Case studies and the governing equations demonstrate the advantages of reverse Monte-Carlo ray tracing (RMCRT) method. The RMCRT can be easily implemented for domain decomposition parallelism. In this dissertation, different efficiency improvements techniques for RMCRT are introduced and implemented. They are the random number generator, stratified sampling, ray-surface intersection calculation, Russian roulette, and important sampling. There are two major modules in solving the radiative heat transfer problems: the RMCRT RTE solver and the optical property models. RMCRT is first fully verified in gray, scattering, absorbing and emitting media with black/nonblack, diffuse/nondiffuse bounded surface problems. Sensitivity analysis is carried out with regard to the ray numbers, the mesh resolutions of the computational domain, optical thickness of the media and effects of variance reduction techniques (stratified sampling, Russian roulette). Results are compared with either analytical solutions or benchmark results. The efficiency (the product of error and computation time) of RMCRT has been compared to DOM and suggest great potential for RMCRT's application
NASA Astrophysics Data System (ADS)
Tichý, Vladimír; Hudec, René; Němcová, Šárka
2016-06-01
The algorithm presented is intended mainly for lobster eye optics. This type of optics (and some similar types) allows for a simplification of the classical ray-tracing procedure that requires great many rays to simulate. The method presented performs the simulation of a only few rays; therefore it is extremely effective. Moreover, to simplify the equations, a specific mathematical formalism is used. Only a few simple equations are used, therefore the program code can be simple as well. The paper also outlines how to apply the method to some other reflective optical systems.
NASA Astrophysics Data System (ADS)
Tichý, Vladimír; Hudec, René; Němcová, Šárka
2016-03-01
The algorithm presented is intended mainly for lobster eye optics. This type of optics (and some similar types) allows for a simplification of the classical ray-tracing procedure that requires great many rays to simulate. The method presented performs the simulation of a only few rays; therefore it is extremely effective. Moreover, to simplify the equations, a specific mathematical formalism is used. Only a few simple equations are used, therefore the program code can be simple as well. The paper also outlines how to apply the method to some other reflective optical systems.
Solar Proton Transport Within an ICRU Sphere Surrounded by a Complex Shield: Ray-trace Geometry
NASA Technical Reports Server (NTRS)
Slaba, Tony C.; Wilson, John W.; Badavi, Francis F.; Reddell, Brandon D.; Bahadori, Amir A.
2015-01-01
A computationally efficient 3DHZETRN code with enhanced neutron and light ion (Z is less than or equal to 2) propagation was recently developed for complex, inhomogeneous shield geometry described by combinatorial objects. Comparisons were made between 3DHZETRN results and Monte Carlo (MC) simulations at locations within the combinatorial geometry, and it was shown that 3DHZETRN agrees with the MC codes to the extent they agree with each other. In the present report, the 3DHZETRN code is extended to enable analysis in ray-trace geometry. This latest extension enables the code to be used within current engineering design practices utilizing fully detailed vehicle and habitat geometries. Through convergence testing, it is shown that fidelity in an actual shield geometry can be maintained in the discrete ray-trace description by systematically increasing the number of discrete rays used. It is also shown that this fidelity is carried into transport procedures and resulting exposure quantities without sacrificing computational efficiency.
Exact ray tracing formulas based on a nontrigonometric alternative to Snell's law.
Elagha, Hassan A
2012-12-01
In this work, Fermat's principle is applied to derive a simple exact formula for refraction (reflection) in terms of the lengths of the incident and refracted rays. This formula is a nontrigonometric alternative to Snell's law and is general for all optical surfaces. It is used to derive the paraxial optics equations in a more simple and direct way than that often used in the literature. It's also applied to derive a new single, exact ray tracing formula for the nonparaxial refraction (reflection) at a single optical surface. The obtained formulas are used to develop a simple ray tracing procedure for meridional refraction through systems of spherical surfaces without the need to use any form of Snell's law. Numerical examples are provided and discussed. PMID:23455919
NASA Astrophysics Data System (ADS)
Huang, Yong; Shi, Guo-Dong; Zhu, Ke-Yong
2016-06-01
This paper adopts the Runge-Kutta ray tracing method to obtain the ray-trajectory numerical solution in a two-dimensional gradient index medium. The emitting, absorbing and scattering processes are simulated by the Monte Carlo method. The temperature field and ray trajectory in the medium are obtained by the three methods, the Runge-Kutta ray tracing method, the ray tracing method with the cell model and the discrete curved ray tracing method with the linear refractive index cell model. Comparing the results of the three methods, it is found that the results by the Monte Carlo Runge-Kutta ray tracing method are of the highest accuracy. To improve the computational speed, the variable step-size Runge-Kutta ray tracing method is proposed, and the maximum relative error between the temperature field in the nonscattering medium by this method and the benchmark solution is less than 0.5%. The results also suggest that the Runge-Kutta ray tracing method would make the radiative transfer solution in the three-dimensional graded index media much easier.
Zheng, Guoyan
2009-01-01
The widely used procedure of evaluation of cup orientation following total hip arthroplasty using single standard anteroposterior (AP) radiograph is known inaccurate, largely due to the wide variability in individual pelvic orientation relative to X-ray plate. 2D/3D rigid image registration methods have been introduced for an accurate determination of the post-operative cup alignment with respect to an anatomical reference extracted from the CT data. Although encouraging results have been reported, their extensive usage in clinical routine is still limited. This may be explained by their requirement of a CAD model of the prosthesis, which is often difficult to be organized from the manufacturer due to the proprietary issue, and by their requirement of a pre-operative CT scan, which is not available for most retrospective studies. To address these issues, we developed and validated a statistically deformable 2D/3D registration approach for accurate determination of post-operative cup orientation. No CAD model and pre-operative CT data is required any more. Quantitative and qualitative results evaluated on cadaveric and clinical datasets are given, which indicate the validity of the approach. PMID:20426064
Trace metal content in aspirin and women's cosmetics via proton induced x-ray emission (PIXE)
Hichwa, B.P.; Pun, D.D.; Wang, D.
1981-04-01
A multielemental analysis to determine the trace metal content of generic and name-brand aspirins and name-brand lipsticks was done via proton induced x-ray (PIXE) measurements. The Hope College PIXE system is described as well as the target preparation methods. The trace metal content of twelve brands of aspirin and aspirin substitutes and fourteen brands of lipstick are reported. Detection limits for most elements are in the range of 100 parts per billion (ppb) to 10 parts per million (ppm).
Yang, Bing-jun; Chao, Keng-hsing; Tsai, Jui-che
2012-09-01
In this paper we develop a three-dimensional (3D) ray tracing tool based on the ABCD ray transfer matrices. With symmetric optical components and under paraxial approximation, two sets of 2×2 ABCD matrices, each for a two-dimensional subspace, can be used to describe the 3D ray propagation completely. Compared to commercial ray-tracing software packages, our tool requires no tedious drawing, and the results for various conditions, such as different device dimensions and incident angles, can be easily obtained by simply changing the parameter values used for the calculation. We have employed this matrix-based 3D ray tracing tool to model cat's eye retroreflectors. The cat's eye performance, including the retroreflection efficiency, acceptance angle (i.e., field of view), and beam divergence and deviation, is fully studied. The application of this 3D ray tracing technique can be further extended to other optical components. PMID:22945148
Wang, Hongkai; Stout, David B; Chatziioannou, Arion F
2013-01-01
The development of sophisticated and high throughput whole body small animal imaging technologies has created a need for improved image analysis and increased automation. The registration of a digital mouse atlas to individual images is a prerequisite for automated organ segmentation and uptake quantification. This paper presents a fully-automatic method for registering a statistical mouse atlas with individual subjects based on an anterior-posterior X-ray projection and a lateral optical photo of the mouse silhouette. The mouse atlas was trained as a statistical shape model based on 83 organ-segmented micro-CT images. For registration, a hierarchical approach is applied which first registers high contrast organs, and then estimates low contrast organs based on the registered high contrast organs. To register the high contrast organs, a 2D-registration-back-projection strategy is used that deforms the 3D atlas based on the 2D registrations of the atlas projections. For validation, this method was evaluated using 55 subjects of preclinical mouse studies. The results showed that this method can compensate for moderate variations of animal postures and organ anatomy. Two different metrics, the Dice coefficient and the average surface distance, were used to assess the registration accuracy of major organs. The Dice coefficients vary from 0.31±0.16 for the spleen to 0.88±0.03 for the whole body, and the average surface distance varies from 0.54±0.06 mm for the lungs to 0.85±0.10 mm for the skin. The method was compared with a direct 3D deformation optimization (without 2D-registration-back-projection) and a single-subject atlas registration (instead of using the statistical atlas). The comparison revealed that the 2D-registration-back-projection strategy significantly improved the registration accuracy, and the use of the statistical mouse atlas led to more plausible organ shapes than the single-subject atlas. This method was also tested with shoulder xenograft
Ionospheric differential error determination using ray tracing for a short baseline
NASA Astrophysics Data System (ADS)
Abdullah, M.; Strangeways, H. J.; Zulkifli, S. S. N.
2010-11-01
Since the United States government discontinued Selective Availability (SA) on 1 May 2000, ionospheric effects have been responsible for the largest errors in GPS systems. The standard Differential GPS (DGPS) method is incapable of completely eliminating the ionospheric error. This paper describes a new approach to determine the differential ionospheric error between geographically distributed receiver stations. The ray paths of GPS signals were simulated using a modified Jones 3D ray tracing programme that includes the effect of the geomagnetic field. A Nelder-Mead optimisation algorithm was embedded in the program to precisely determine the satellite-to-station path. A realistic ionospheric model is essential for accurate ray tracing results and for estimates of differential error that are accurate on sub-centimetre scales. Here, the ionospheric model used in the ray tracing programme was developed by fitting realistic ionosphere profiles with a number of exponential functions. Results were compared to the theoretical approach. Results show that the differential delay is about 1-5 cm at low elevation angles for a short baseline of 10 km, as reported in other literature. This delay is often neglected in DGPS application. The differential delay also shows a pattern similar to that predicted by the Klobuchar model. The method proposed here can be used to improve future GPS applications.
Weiland, C.M.; Steck, L.K.; Dawson, P.B.
1995-10-10
The authors explore the impact of three-dimensional minimum travel time ray tracing on nonlinear teleseismic inversion. This problem has particular significance when trying to image strongly contrasting low-velocity bodies, such as magma chambers, because strongly refracted/and/or diffracted rays may precede the direct P wave arrival traditionally used in straight-ray seismic tomography. They use a simplex-based ray tracer to compute the three-dimensional, minimum travel time ray paths and employ an interative technique to cope with nonlinearity. Results from synthetic data show that their algorithm results in better model reconstructions compared with traditional straight-ray inversions. The authors reexamine the teleseismic data collected at Long Valley caldera by the U.S. Geological Survey. The most prominent feature of their result is a 25-30% low-velocity zone centered at 11.5 km depth beneath the northwestern quandrant of the caldera. Beneath this at a depth of 24.5 km is a more diffuse 15% low-velocity zone. In general, the low velocities tend to deepen to the south and east. The authors interpret the shallow feature to be the residual Long Valley caldera magma chamber, while the deeper feature may represent basaltic magmas ponded in the midcrust. The deeper position of the prominent low-velocity region in comparison to earlier tomographic images is a result of using three-dimensional rays rather than straight rays in the ray tracing. The magnitude of the low-velocity anomaly is a factor of {approximately}3 times larger than earlier models from linear arrival time inversions and is consistent with models based on observations of ray bending at sites within the caldera. These results imply the presence of anywhere from 7 to 100% partial melt beneath the caldera. 40 refs., 1 fig., 1 tab.
Improvements of the Ray-Tracing Based Method Calculating Hypocentral Loci for Earthquake Location
NASA Astrophysics Data System (ADS)
Zhao, A. H.
2014-12-01
Hypocentral loci are very useful to reliable and visual earthquake location. However, they can hardly be analytically expressed when the velocity model is complex. One of methods numerically calculating them is based on a minimum traveltime tree algorithm for tracing rays: a focal locus is represented in terms of ray paths in its residual field from the minimum point (namely initial point) to low residual points (referred as reference points of the focal locus). The method has no restrictions on the complexity of the velocity model but still lacks the ability of correctly dealing with multi-segment loci. Additionally, it is rather laborious to set calculation parameters for obtaining loci with satisfying completeness and fineness. In this study, we improve the ray-tracing based numerical method to overcome its advantages. (1) Reference points of a hypocentral locus are selected from nodes of the model cells that it goes through, by means of a so-called peeling method. (2) The calculation domain of a hypocentral locus is defined as such a low residual area that its connected regions each include one segment of the locus and hence all the focal locus segments are respectively calculated with the minimum traveltime tree algorithm for tracing rays by repeatedly assigning the minimum residual reference point among those that have not been traced as an initial point. (3) Short ray paths without branching are removed to make the calculated locus finer. Numerical tests show that the improved method becomes capable of efficiently calculating complete and fine hypocentral loci of earthquakes in a complex model.
TIM, a ray-tracing program for METATOY research and its dissemination
NASA Astrophysics Data System (ADS)
Lambert, Dean; Hamilton, Alasdair C.; Constable, George; Snehanshu, Harsh; Talati, Sharvil; Courtial, Johannes
2012-03-01
TIM (The Interactive METATOY) is a ray-tracing program specifically tailored towards our research in METATOYs, which are optical components that appear to be able to create wave-optically forbidden light-ray fields. For this reason, TIM possesses features not found in other ray-tracing programs. TIM can either be used interactively or by modifying the openly available source code; in both cases, it can easily be run as an applet embedded in a web page. Here we describe the basic structure of TIM's source code and how to extend it, and we give examples of how we have used TIM in our own research. Program summaryProgram title: TIM Catalogue identifier: AEKY_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEKY_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: GNU General Public License No. of lines in distributed program, including test data, etc.: 124 478 No. of bytes in distributed program, including test data, etc.: 4 120 052 Distribution format: tar.gz Programming language: Java Computer: Any computer capable of running the Java Virtual Machine (JVM) 1.6 Operating system: Any; developed under Mac OS X Version 10.6 RAM: Typically 145 MB (interactive version running under Mac OS X Version 10.6) Classification: 14, 18 External routines: JAMA [1] (source code included) Nature of problem: Visualisation of scenes that include scene objects that create wave-optically forbidden light-ray fields. Solution method: Ray tracing. Unusual features: Specifically designed to visualise wave-optically forbidden light-ray fields; can visualise ray trajectories; can visualise geometric optic transformations; can create anaglyphs (for viewing with coloured "3D glasses") and random-dot autostereograms of the scene; integrable into web pages. Running time: Problem-dependent; typically seconds for a simple scene.
Jones, M. W. M.; Elgass, K. D.; Junker, M. D.; de Jonge, M. D.; van Riessen, G. A.
2016-01-01
Recent developments in biological X-ray microscopy have allowed structural information and elemental distribution to be simultaneously obtained by combining X-ray ptychography and X-ray fluorescence microscopy. Experimentally, these methods can be performed simultaneously; however, the optimal conditions for each measurement may not be compatible. Here, we combine two distinct measurements of ultrastructure and elemental distribution, with each measurement performed under optimised conditions. By combining optimised ptychography and fluorescence information we are able to determine molar concentrations from two-dimensional images, allowing an investigation into the interactions between the environment sensing filopodia in fibroblasts and extracellular calcium. Furthermore, the biological ptychography results we present illustrate a point of maturity where the technique can be applied to solve significant problems in structural biology. PMID:27067957
NASA Astrophysics Data System (ADS)
Jones, M. W. M.; Elgass, K. D.; Junker, M. D.; de Jonge, M. D.; van Riessen, G. A.
2016-04-01
Recent developments in biological X-ray microscopy have allowed structural information and elemental distribution to be simultaneously obtained by combining X-ray ptychography and X-ray fluorescence microscopy. Experimentally, these methods can be performed simultaneously; however, the optimal conditions for each measurement may not be compatible. Here, we combine two distinct measurements of ultrastructure and elemental distribution, with each measurement performed under optimised conditions. By combining optimised ptychography and fluorescence information we are able to determine molar concentrations from two-dimensional images, allowing an investigation into the interactions between the environment sensing filopodia in fibroblasts and extracellular calcium. Furthermore, the biological ptychography results we present illustrate a point of maturity where the technique can be applied to solve significant problems in structural biology.
Combining ray tracing and CFD in the thermal analysis of a parabolic dish tubular cavity receiver
NASA Astrophysics Data System (ADS)
Craig, Ken J.; Marsberg, Justin; Meyer, Josua P.
2016-05-01
This paper describes the numerical evaluation of a tubular receiver used in a dish Brayton cycle. In previous work considering the use of Computational Fluid Dynamics (CFD) to perform the calculation of the absorbed radiation from the parabolic dish into the cavity as well as the resulting conjugate heat transfer, it was shown that an axi-symmetric model of the dish and receiver absorbing surfaces was useful in reducing the computational cost required for a full 3-D discrete ordinates solution, but concerns remained about its accuracy. To increase the accuracy, the Monte Carlo ray tracer SolTrace is used to perform the calculation of the absorbed radiation profile to be used in the conjugate heat transfer CFD simulation. The paper describes an approach for incorporating a complex geometry like a tubular receiver generated using CFD software into SolTrace. The results illustrate the variation of CFD mesh density that translates into the number of elements in SolTrace as well as the number of rays used in the Monte Carlo approach and their effect on obtaining a resolution-independent solution. The conjugate heat transfer CFD simulation illustrates the effect of applying the SolTrace surface heat flux profile solution as a volumetric heat source to heat up the air inside the tube. Heat losses due to convection and thermal re-radiation are also determined as a function of different tube absorptivities.
Analytical calculation of spectral phase of grism pairs by the geometrical ray tracing method
NASA Astrophysics Data System (ADS)
Rahimi, L.; Askari, A. A.; Saghafifar, H.
2016-07-01
The most optimum operation of a grism pair is practically approachable when an analytical expression of its spectral phase is in hand. In this paper, we have employed the accurate geometrical ray tracing method to calculate the analytical phase shift of a grism pair, at transmission and reflection configurations. As shown by the results, for a great variety of complicated configurations, the spectral phase of a grism pair is in the same form of that of a prism pair. The only exception is when the light enters into and exits from different facets of a reflection grism. The analytical result has been used to calculate the second-order dispersions of several examples of grism pairs in various possible configurations. All results are in complete agreement with those from ray tracing method. The result of this work can be very helpful in the optimal design and application of grism pairs at various configurations.
Nonparaxial geometrical Ronchi test for spherical mirrors: an inverse ray-tracing approach.
Juarez-Salazar, Rigoberto
2016-08-01
A geometrical model based on an inverse ray-tracing approach to describe the Ronchi test for a concave spherical mirror is presented. In contrast to the conventional ray-tracing method, which refers to information unavailable in ronchigrams, the proposed model provides an explicit relation between the available information in the ronchigram and the parameters of the setup (radius of the sphere, position of the source, position and orientation of the observation, and grating planes). This allows for extracting the parameters of interest by a simple fitting procedure, as demonstrated by an application. The derived model exhibits new unexplored potential applications of the Ronchi test, establishing it as a very useful, simple, and universal tool for optical evaluation. PMID:27505380
Three-dimensional ray tracing of the Jovian magnetosphere in the low-frequency range
NASA Technical Reports Server (NTRS)
Menietti, J. D.; Green, J. L.; Gulkis, S.; Six, F.
1984-01-01
Three-dimensional ray tracing of the Jovian DAM emission has been performed utilizing the O-4 magnetic field model (Acuna and Ness, 1979) and a realistic plasma model. Minimal assumptions about the emission mechanism have been made that include radiation in the right-hand extraordinary mode, propagating nearly perpendicular to the field line at source points located just above the RX cutoff frequency along Io flux tubes. Ray tracing has been performed in the frequency range from 2-35 MHz from successive Io flux tubes separated by ten degrees of central meridian longitude for a full circumference of northern hemisphere sources. The results show unusual complexity in model arc spectra that is displayed in a constant Io phase format with many similarities to the Voyager PRA data. The results suggest much of the variation in observed DAM spectral features is a result of propagation effects rather than emission process differences.
NASA Technical Reports Server (NTRS)
Baumeister, Joseph F.
1990-01-01
Analysis of energy emitted from simple or complex cavity designs can lead to intricate solutions due to nonuniform radiosity and irradiation within a cavity. A numerical ray tracing technique was applied to simulate radiation propagating within and from various cavity designs. To obtain the energy balance relationships between isothermal and nonisothermal cavity surfaces and space, the computer code NEVADA was utilized for its statistical technique applied to numerical ray tracing. The analysis method was validated by comparing results with known theoretical and limiting solutions, and the electrical resistance network method. In general, for nonisothermal cavities the performance (apparent emissivity) is a function of cylinder length-to-diameter ratio, surface emissivity, and cylinder surface temperatures. The extent of nonisothermal conditions in a cylindrical cavity significantly affects the overall cavity performance. Results are presented over a wide range of parametric variables for use as a possible design reference.
Stochastic ray tracing for simulation of high intensity focal ultrasound therapy.
Koskela, Julius; Vahala, Erkki; de Greef, Martijn; Lafitte, Luc P; Ries, Mario
2014-09-01
An algorithm is presented for rapid simulation of high-intensity focused ultrasound (HIFU) fields. Essentially, the method combines ray tracing with Monte Carlo integration to evaluate the Rayleigh-Sommerfeld integral. A large number of computational particles, phonons, are distributed among the elements of a phase-array transducer. The phonons are emitted into random directions and are propagated along trajectories computed with the ray tracing method. As the simulation progresses, an improving stochastic estimate of the acoustic field is obtained. The method can adapt to complicated geometries, and it is well suited to parallelization. The method is verified against reference simulations and pressure measurements from an ex vivo porcine thoracic tissue sample. Results are presented for acceleration with graphics processing units (GPUs). The method is expected to serve in applications, where flexibility and rapid computation time are crucial, in particular clinical HIFU treatment planning. PMID:25190416
Hopp, T; Duric, N; Ruiter, N V
2015-03-01
Ultrasound Computer Tomography (USCT) is a promising breast imaging modality under development. Comparison to a standard method like mammography is essential for further development. Due to significant differences in image dimensionality and compression state of the breast, correlating USCT images and X-ray mammograms is challenging. In this paper we present a 2D/3D registration method to improve the spatial correspondence and allow direct comparison of the images. It is based on biomechanical modeling of the breast and simulation of the mammographic compression. We investigate the effect of including patient-specific material parameters estimated automatically from USCT images. The method was systematically evaluated using numerical phantoms and in-vivo data. The average registration accuracy using the automated registration was 11.9mm. Based on the registered images a method for analysis of the diagnostic value of the USCT images was developed and initially applied to analyze sound speed and attenuation images based on X-ray mammograms as ground truth. Combining sound speed and attenuation allows differentiating lesions from surrounding tissue. Overlaying this information on mammograms, combines quantitative and morphological information for multimodal diagnosis. PMID:25456144
Mathematic models for a ray tracing method and its applications in wireless optical communications.
Zhang, Minglun; Zhang, Yangan; Yuan, Xueguang; Zhang, Jinnan
2010-08-16
This paper presents a new ray tracing method, which contains a whole set of mathematic models, and its validity is verified by simulations. In addition, both theoretical analysis and simulation results show that the computational complexity of the method is much lower than that of previous ones. Therefore, the method can be used to rapidly calculate the impulse response of wireless optical channels for complicated systems. PMID:20721238
Detection and quantification of trace elements in rice and rice products using x-ray fluorescence
NASA Astrophysics Data System (ADS)
Foran, Kelly A.; Fleming, David E. B.
2015-12-01
We used X-ray fluorescence (XRF) to examine the presence of arsenic (As) and other trace elements (manganese, iron, nickel, copper, and zinc) in rice and rice products. A portable XRF analyzer was used to test samples, and amplitudes for the analyzed elements were identified in the resulting data. The detection limit of the system was sufficiently low to detect As in some rice and rice product samples.
Leiner, Claude; Nemitz, Wolfgang; Schweitzer, Susanne; Kuna, Ladislav; Wenzl, Franz P; Hartmann, Paul; Satzinger, Valentin; Sommer, Christian
2016-03-20
We show that with an appropriate combination of two optical simulation techniques-classical ray-tracing and the finite difference time domain method-an optical device containing multiple diffractive and refractive optical elements can be accurately simulated in an iterative simulation approach. We compare the simulation results with experimental measurements of the device to discuss the applicability and accuracy of our iterative simulation procedure. PMID:27140556
Modeling pyramidal sensors in ray-tracing software by a suitable user-defined surface
NASA Astrophysics Data System (ADS)
Antichi, Jacopo; Munari, Matteo; Magrin, Demetrio; Riccardi, Armando
2016-04-01
Following the unprecedented results in terms of performances delivered by the first light adaptive optics system at the Large Binocular Telescope, there has been a wide-spread and increasing interest on the pyramid wavefront sensor (PWFS), which is the key component, together with the adaptive secondary mirror, of the adaptive optics (AO) module. Currently, there is no straightforward way to model a PWFS in standard sequential ray-tracing software. Common modeling strategies tend to be user-specific and, in general, are unsatisfactory for general applications. To address this problem, we have developed an approach to PWFS modeling based on user-defined surface (UDS), whose properties reside in a specific code written in C language, for the ray-tracing software ZEMAX™. With our approach, the pyramid optical component is implemented as a standard surface in ZEMAX™, exploiting its dynamic link library (DLL) conversion then greatly simplifying ray tracing and analysis. We have utilized the pyramid UDS DLL surface-referred to as pyramidal acronyms may be too risky (PAM2R)-in order to design the current PWFS-based AO system for the Giant Magellan Telescope, evaluating tolerances, with particular attention to the angular sensitivities, by means of sequential ray-tracing tools only, thus verifying PAM2R reliability and robustness. This work indicates that PAM2R makes the design of PWFS as simple as that of other optical standard components. This is particularly suitable with the advent of the extremely large telescopes era for which complexity is definitely one of the main challenges.
NASA Astrophysics Data System (ADS)
Lavrinenko, Ya S.; Morozov, I. V.; Pikuz, S. A.; Skobelev, I. Yu
2015-11-01
Spherically bent crystals are widely used in focusing monochromators, spectrometers and other x-ray optical systems. In particular, they are used in focusing spectrometers with spatial resolution, applied in high energy density diagnostics and warm dense matter studies. In this case, plasma parameters are obtained via measurements of relative intensities of characteristic spectral emission lines for multiply charged ions, which are affected by an instrumental function. Here we develop and use the ray-tracing computer simulations to study reflectivity properties of spherically bent crystals in a particular experimental conditions and to provide the method to adjust and validate the measured spectral line intensities on quantitative basis.
A new kind of splines and their use for fast ray-tracing in reflective cavities
NASA Astrophysics Data System (ADS)
Pantelic, Dejan V.; Janevski, Zoran D.
1989-08-01
In this paper we are presenting a new kind of splines that are very effective in ray-tracing applications. They are designed in such a way to enable the fast and efficient computation of line-spline intersections (line representing the light ray, and spline representing the reflective cavity). These splines are piecewise parabolic polynomials, but with additional degrees of freedom. Polynomial sections of the spline can be rotated to a certain angle (each section has its own angle of rotation), enabling thus the continuity of the first derivative.
Stress optical path difference analysis of off-axis lens ray trace footprint
NASA Astrophysics Data System (ADS)
Hsu, Ming-Ying; Chan, Chia-Yen; Lin, Wei-Cheng; Wu, Kun-Huan; Chen, Chih-Wen; Chan, Shenq-Tsong; Huang, Ting-Ming
2013-06-01
The mechanical and thermal stress on lens will cause the glass refractive index different, the refractive index of light parallel and light perpendicular to the direction of stress. The refraction index changes will introduce Optical Path Difference (OPD). This study is applying Finite Element Method (FEM) and optical ray tracing; calculate off axis ray stress OPD. The optical system stress distribution result is calculated from finite element simulation, and the stress coordinate need to rotate to optical path direction. Meanwhile, weighting stress to each optical ray path and sum the ray path OPD. The Z-direction stress OPD can be fitted by Zernike polynomial, the separated to sag difference, and rigid body motion. The fitting results can be used to evaluate the stress effect on optical component.
Chorus wave-normal statistics in the Earth's radiation belts from ray tracing technique
NASA Astrophysics Data System (ADS)
Breuillard, H.; Zaliznyak, Y.; Krasnoselskikh, V.; Agapitov, O.; Artemyev, A.; Rolland, G.
2012-08-01
Discrete ELF/VLF (Extremely Low Frequency/Very Low Frequency) chorus emissions are one of the most intense electromagnetic plasma waves observed in radiation belts and in the outer terrestrial magnetosphere. These waves play a crucial role in the dynamics of radiation belts, and are responsible for the loss and the acceleration of energetic electrons. The objective of our study is to reconstruct the realistic distribution of chorus wave-normals in radiation belts for all magnetic latitudes. To achieve this aim, the data from the electric and magnetic field measurements onboard Cluster satellite are used to determine the wave-vector distribution of the chorus signal around the equator region. Then the propagation of such a wave packet is modeled using three-dimensional ray tracing technique, which employs K. Rönnmark's WHAMP to solve hot plasma dispersion relation along the wave packet trajectory. The observed chorus wave distributions close to waves source are first fitted to form the initial conditions which then propagate numerically through the inner magnetosphere in the frame of the WKB approximation. Ray tracing technique allows one to reconstruct wave packet properties (electric and magnetic fields, width of the wave packet in k-space, etc.) along the propagation path. The calculations show the spatial spreading of the signal energy due to propagation in the inhomogeneous and anisotropic magnetized plasma. Comparison of wave-normal distribution obtained from ray tracing technique with Cluster observations up to 40° latitude demonstrates the reliability of our approach and applied numerical schemes.
A data distributed, parallel algorithm for ray-traced volume rendering
Ma, Kwan-Liu; Painter, J.S.; Hansen, C.D.; Krogh, M.F.
1993-03-30
This paper presents a divide-and-conquer ray-traced volume rendering algorithm and its implementation on networked workstations and a massively parallel computer, the Connection Machine CM-5. This algorithm distributes the data and the computational load to individual processing units to achieve fast, high-quality rendering of high-resolution data, even when only a modest amount of memory is available on each machine. The volume data, once distributed, is left intact. The processing nodes perform local ray-tracing of their subvolume concurrently. No communication between processing units is needed during this locally ray-tracing process. A subimage is generated by each processing unit and the final image is obtained by compositing subimages in the proper order, which can be determined a priori. Implementations and tests on a group of networked workstations and on the Thinking Machines CM-5 demonstrate the practicality of our algorithm and expose different performance tuning issues for each platform. We use data sets from medical imaging and computational fluid dynamics simulations in the study of this algorithm.
NASA Technical Reports Server (NTRS)
Flynn, G. J.; Sutton, S. R.
1989-01-01
Trace element analyses were performed on bulk cosmic dust particles by Proton Induced X Ray Emission (PIXE) and Synchrotron X Ray Fluorescence (SXRF). When present at or near chondritic abundances the trace elements K, Ti, Cr, Mn, Cu, Zn, Ga, Ge, Se, and Br are presently detectable by SXRF in particles of 20 micron diameter. Improvements to the SXRF analysis facility at the National Synchrotron Light Source presently underway should increase the range of detectable elements and permit the analysis of smaller samples. In addition the Advanced Photon Source will be commissioned at Argonne National Laboratory in 1995. This 7 to 8 GeV positron storage ring, specifically designed for high-energy undulator and wiggler insertion devices, will be an ideal source for an x ray microprobe with one micron spatial resolution and better than 100 ppb elemental sensitivity for most elements. Thus trace element analysis of individual micron-sized grains should be possible by the time of the comet nucleus sample return mission.
Ray tracing in the human eye: measurement and modeling of optical aberrations
NASA Astrophysics Data System (ADS)
Navarro, Rafael M.; Rodriguez, P.; Gonzalez, L.; Aporta, J.; Hdez-Matamoros, J. L.
2004-10-01
The rapid development of cataract and refractive surgery requires new methods to assess the optical quality of the eye. The optimized optical design of custom treatments to improve the optical performance of individual eyes requires, at least, to have the technology to (1) measure the geometry (anatomy) of the optics of the eye; (2) measure the optical performance (refractive state, aberrations, etc); (3) Build a custom optical and anatomical model of the individual eye to treat; (4) Optimal design of custom treatments. In this communication we will present the work carried out by our group to develop methods for measuring and modeling the optical performance of the eye. In particular, we will focus, first, on the Laser Ray Tracing method that we have developed to measure the optical aberrations of the eye, as a physical in vivo implementation of the classical numerical ray tracing used by optical designers; and second, on the development of custom optical models of the eye to perform that numerical ray tracing which predicts with a high fidelity experimental measurements. The methods developed have been applied to design both custom surgery and optical aids to improve optical performance.
Ray-tracing and physical-optics analysis of the aperture efficiency in a radio telescope.
Olmi, Luca; Bolli, Pietro
2007-07-01
The performance of telescope systems working at microwave or visible-IR wavelengths is typically described in terms of different parameters according to the wavelength range. Most commercial ray-tracing packages have been specifically designed for use with visible-IR systems and thus, though very flexible and sophisticated, do not provide the appropriate parameters to fully describe microwave antennas and to compare with specifications. We demonstrate that the Strehl ratio is equal to the phase efficiency when the apodization factor is taken into account. The phase efficiency is the most critical contribution to the aperture efficiency of an antenna and the most difficult parameter to optimize during the telescope design. The equivalence between the Strehl ratio and the phase efficiency gives the designer/user of the telescope the opportunity to use the faster commercial ray-tracing software to optimize the design. We also discuss the results of several tests performed to check the validity of this relationship that we carried out using a ray-tracing software, ZEMAX, and a full Physical Optics software, GRASP9.3, applied to three different telescope designs that span a factor of approximately 10 in terms of D/lambda. The maximum measured discrepancy between phase efficiency and Strehl ratio varies between approximately 0.4% and 1.9% up to an offset angle of >40 beams, depending on the optical configuration, but it is always less than 0.5% where the Strehl ratio is >0.95. PMID:17571151
Inversion of anisotropic inner core structure from three dimensional ray tracing
NASA Astrophysics Data System (ADS)
Sun, X.; Song, X.
2005-12-01
Seismological studies have generally suggest that the Earth's inner core is anisotropic and the anisotropic structure change significantly both laterally and with depth. Previous body-wave studies of the inner core have relied on 1-D ray tracing or waveform modeling, which do not account fully the 3D anisotropic structure. Here we adopt a pseudo-bending ray tracing (PBR) method in spherical coordinates (Koketsu and Sekine, 1998) for seismic rays that traverse the inner core (PKP-DF phase). The method iteratively perturbs each discontinuity points and continuous segment of the ray through 3D (but isotropic) earth structure so that its travel time is minimum. Our implementation also includes a flexible scheme in calculating the velocity gradient needed to perturb the ray. A large volume is included in calculating the velocity gradient initially to find the global minimum, but a small volume surrounding the ray is used eventually to obtain the precise local velocity gradient that is sampled by the ray. Tests show that our implementation is very stable, reliable, and fast. We have traced the rays for over 3000 event-station pairs that we have differential PKP travel-time measurements using both the PBR method and a shooting method for a 1D model (AK135). The travel-time difference from the two methods is generally within 0.05 s with a few up to 0.07 s and the largest path difference is within 24 km; Even with a model of strong velocity gradient, the travel time difference is still less than 0.08s and the largest path difference is within 40km. Because the ray direction in the inner core does not change much (within 10 degrees even with a strong velocity gradient in the inner core), the 3D anisotropic structure of the inner core can be approximated to the first order as 3D heterogeneous (but isotropic) structure for a given ray, assuming the inner core anisotropy is axisymmetric. We are implementing the PBR method and B-spline interpolation to invert for 3D anisotropic
NASA Astrophysics Data System (ADS)
Pawlowski, Jason M.; Ding, George X.
2014-04-01
A new model-based dose calculation algorithm is presented for kilovoltage x-rays and is tested for the cases of calculating the radiation dose from kilovoltage cone-beam CT (kV-CBCT) and 2D planar projected radiographs. This algorithm calculates the radiation dose to water-like media as the sum of primary and scattered dose components. The scatter dose is calculated by convolution of a newly introduced, empirically parameterized scatter dose kernel with the primary photon fluence. Several approximations are introduced to increase the scatter dose calculation efficiency: (1) the photon energy spectrum is approximated as monoenergetic; (2) density inhomogeneities are accounted for by implementing a global distance scaling factor in the scatter kernel; (3) kernel tilting is ignored. These approximations allow for efficient calculation of the scatter dose convolution with the fast Fourier transform. Monte Carlo simulations were used to obtain the model parameters. The accuracy of using this model-based algorithm was validated by comparing with the Monte Carlo method for calculating dose distributions for real patients resulting from radiotherapy image guidance procedures including volumetric kV-CBCT scans and 2D planar projected radiographs. For all patients studied, mean dose-to-water errors for kV-CBCT are within 0.3% with a maximum standard deviation error of 4.1%. Using a medium-dependent correction method to account for the effects of photoabsorption in bone on the dose distribution, mean dose-to-medium errors for kV-CBCT are within 3.6% for bone and 2.4% for soft tissues. This algorithm offers acceptable accuracy and has the potential to extend the applicability of model-based dose calculation algorithms from megavoltage to kilovoltage photon beams.
Hapca, Simona; Baveye, Philippe C.; Wilson, Clare; Lark, Richard Murray; Otten, Wilfred
2015-01-01
There is currently a significant need to improve our understanding of the factors that control a number of critical soil processes by integrating physical, chemical and biological measurements on soils at microscopic scales to help produce 3D maps of the related properties. Because of technological limitations, most chemical and biological measurements can be carried out only on exposed soil surfaces or 2-dimensional cuts through soil samples. Methods need to be developed to produce 3D maps of soil properties based on spatial sequences of 2D maps. In this general context, the objective of the research described here was to develop a method to generate 3D maps of soil chemical properties at the microscale by combining 2D SEM-EDX data with 3D X-ray computed tomography images. A statistical approach using the regression tree method and ordinary kriging applied to the residuals was developed and applied to predict the 3D spatial distribution of carbon, silicon, iron, and oxygen at the microscale. The spatial correlation between the X-ray grayscale intensities and the chemical maps made it possible to use a regression-tree model as an initial step to predict the 3D chemical composition. For chemical elements, e.g., iron, that are sparsely distributed in a soil sample, the regression-tree model provides a good prediction, explaining as much as 90% of the variability in some of the data. However, for chemical elements that are more homogenously distributed, such as carbon, silicon, or oxygen, the additional kriging of the regression tree residuals improved significantly the prediction with an increase in the R2 value from 0.221 to 0.324 for carbon, 0.312 to 0.423 for silicon, and 0.218 to 0.374 for oxygen, respectively. The present research develops for the first time an integrated experimental and theoretical framework, which combines geostatistical methods with imaging techniques to unveil the 3-D chemical structure of soil at very fine scales. The methodology presented
Pacella, D.; Pizzicaroli, G.; Romano, A.; Gabellieri, L.; Bellazzini, R.; Brez, A.
2008-03-12
Soft-X ray 2-D imaging on ITER is not considered yet. We propose a new approach, based on a gas detector with a gas electron multiplier (GEM) as amplifying structure and with a two-dimensional readout fully integrated with the front end electronics, through an ASIC developed on purpose. The concept has been already tested by means of a prototype, with 128 pixels, carried out in Frascati in collaboration with INFN-Pisa and tested on FTU in 2001 and NSTX in 2002-2004. Thanks to the photon counting mode, it provides 2-D imaging with high time resolution (sub millisecond), high sensitivity and signal to noise ratio. Its capability of energy discrimination allows the acquisition of pictures in X-ray energy bands or to perform a spectral scan in the full energy interval. We propose the realisation of such kind a detector with a readout microchip (ASIC) equipped with 105600 hexagonal pixels arranged at 70 {mu}m pitch in a 300x352 honeycomb matrix, corresponding to an active area of 2.1x2.1 cm{sup 2}, with a pixel density of 240 pixels/ mm{sup 2}. Each pixel is connected to a charge sensitive amplifier followed by a discriminator of pulse amplitude and counter. The chip integrates more than 16.5 million transistors and it is subdivided in 64 identical clusters, to be read independently each other. An important part of the work will be also the design of the whole detector to fulfil all the constraints and requirements as plasma diagnostic in a tokamak machine. Since the detector has high and controllable intrinsic gain, it works well even at very low photon energy, ranging from 0.2 keV to 10 keV (X-VUV region). This range appears therefore particularly suitable for ITER to monitor the outer part of the plasma. In particular pedestal physics, edge modes, localization and effects of additional heating, boundary plasma control etc. The capability of this proposed detector to work in this energy range is further valuable because solid state detectors are not favorite at low
Fast robust non-sequential optical ray-tracing with implicit algebraic surfaces
NASA Astrophysics Data System (ADS)
Greynolds, Alan W.
2015-09-01
The fastest, most robust, general technique for non-sequentially ray-tracing a large class of imaging and non-imaging optical systems is by geometric modeling with algebraic (i.e. polynomial) implicit surfaces. The basic theory of these surfaces with special attention to optimizing their precise intersection with a ray (even at grazing incidence) is outlined for an admittedly limited software implementation. On a couple of "tame" examples, a 64-bit Windows 7 version is significantly faster than the fastest commercial design software (all multi-threaded). Non-sequential ray-surface interactions approaching 30M/sec are achieved on a 12-core 2.67 GHz Mac Pro desktop computer. For a more exotic example of a 6th degree Wood's horn beam dump (light trap), a 32-bit Windows single thread version traces rays nearly 4 times faster than the commercial ASAP software's implicit algebraic surface and over 13 times faster than its equivalent NURBS surface. However, implicit surfaces are foreign to most CAD systems and thus unfortunately, don't easily fit into a modern workflow.
Comparing FDTD and Ray-Tracing Models in Numerical Simulation of HgCdTe LWIR Photodetectors
NASA Astrophysics Data System (ADS)
Vallone, Marco; Goano, Michele; Bertazzi, Francesco; Ghione, Giovanni; Schirmacher, Wilhelm; Hanna, Stefan; Figgemeier, Heinrich
2016-04-01
We present a simulation study of HgCdTe-based long-wavelength infrared detectors, focusing on methodological comparisons between the finite-difference time-domain (FDTD) and ray-tracing optical models. We performed three-dimensional simulations to determine the absorbed photon density distributions and the corresponding photocurrent and quantum efficiency spectra of isolated n-on-p uniform-composition pixels, systematically comparing the results obtained with FDTD and ray tracing. Since ray tracing is a classical optics approach, unable to describe interference effects, its applicability has been found to be strongly wavelength dependent, especially when reflections from metallic layers are relevant. Interesting cavity effects around the material cutoff wavelength are described, and the cases where ray tracing can be considered a viable approximation are discussed.
Comparing FDTD and Ray-Tracing Models in Numerical Simulation of HgCdTe LWIR Photodetectors
NASA Astrophysics Data System (ADS)
Vallone, Marco; Goano, Michele; Bertazzi, Francesco; Ghione, Giovanni; Schirmacher, Wilhelm; Hanna, Stefan; Figgemeier, Heinrich
2016-09-01
We present a simulation study of HgCdTe-based long-wavelength infrared detectors, focusing on methodological comparisons between the finite-difference time-domain (FDTD) and ray-tracing optical models. We performed three-dimensional simulations to determine the absorbed photon density distributions and the corresponding photocurrent and quantum efficiency spectra of isolated n-on- p uniform-composition pixels, systematically comparing the results obtained with FDTD and ray tracing. Since ray tracing is a classical optics approach, unable to describe interference effects, its applicability has been found to be strongly wavelength dependent, especially when reflections from metallic layers are relevant. Interesting cavity effects around the material cutoff wavelength are described, and the cases where ray tracing can be considered a viable approximation are discussed.
Line shape and ray trace calculations in saturated X-ray lasers: Application to Ni-like silver
NASA Astrophysics Data System (ADS)
Benredjem, D.; Guilbaud, O.; Möller, C.; Klisnick, A.; Ros, D.; Dubau, J.; Calisti, A.; Talin, B.
2006-05-01
Longitudinal coherence length in X-ray lasers depends strongly on the shape of the amplified line. We have modelled an experiment performed at the LULI facility of Ecole Polytechnique. The experiment was devoted to the study of the temporal (longitudinal) coherence of the transient Ni-like silver 4d 4p transition X-ray laser at 13.9 nm. Accurate line shape calculations using PPP, a spectral line shape code, confirm that the Voigt profile is a good approximation for this X-ray laser line. This allows us to extensively use the Voigt shape in conditions where the amplifier, i.e. the plasma produced by the interaction of a high intensity laser with a slab target, is neither stationary nor homogeneous. Our calculations involve a ray trace code which is a post-processor to the hydrodynamic simulation EHYBRID. As the effect of saturation is important for the level populations and gains we include the interaction between the amplified beam and the medium using the Maxwell-Bloch formalism. While the FWHM of the spontaneous emission profile is ˜10 mÅ, the amplified X-ray line exhibits gain narrowing leading to the smaller width ˜3 mÅ. Comparison with experiment is discussed.
NASA Astrophysics Data System (ADS)
Jha, R. M.; Wiesbeck, W.
1995-04-01
A generalized approach to analytical surface-ray tracing in three dimensions, and a review of its application to convex conducting bodies, is presented, using the Eisenhart Coordinate System. The ray-parameters so obtained, for quadric cylinders (QUACYLs) and surfaces of revolution (QUASORs), are in a one-parameter form for UTD mutual-coupling applications. The ray analysis is also extended to the hybrid QUACYLs (e.g., aircraft wings) and hybrid QUASORs (e.g., satellite-launch vehicles), by introducing Hertz's principle of particle dynamics to EM theory. This mathematical formulation is applicable even to other important non-Eisenhart surfaces, such as the ogive. A summary of the mathematical formulations is included.
R-LODs: fast LOD-based ray tracing of massive models
Yoon, Sung-Eui; Lauterbach, Christian; Manocha, Dinesh
2006-08-25
We present a novel LOD (level-of-detail) algorithm to accelerate ray tracing of massive models. Our approach computes drastic simplifications of the model and the LODs are well integrated with the kd-tree data structure. We introduce a simple and efficient LOD metric to bound the error for primary and secondary rays. The LOD representation has small runtime overhead and our algorithm can be combined with ray coherence techniques and cache-coherent layouts to improve the performance. In practice, the use of LODs can alleviate aliasing artifacts and improve memory coherence. We implement our algorithm on both 32bit and 64bit machines and able to achieve up to 2.20 times improvement in frame rate of rendering models consisting of tens or hundreds of millions of triangles with little loss in image quality.
A computer program to trace seismic ray distribution in complex two-dimensional geological models
Yacoub, Nazieh K.; Scott, James H.
1970-01-01
A computer program has been developed to trace seismic rays and their amplitudes and energies through complex two-dimensional geological models, for which boundaries between elastic units are defined by a series of digitized X-, Y-coordinate values. Input data for the program includes problem identification, control parameters, model coordinates and elastic parameter for the elastic units. The program evaluates the partitioning of ray amplitude and energy at elastic boundaries, computes the total travel time, total travel distance and other parameters for rays arising at the earth's surface. Instructions are given for punching program control cards and data cards, and for arranging input card decks. An example of printer output for a simple problem is presented. The program is written in FORTRAN IV language. The listing of the program is shown in the Appendix, with an example output from a CDC-6600 computer.
GRay: A Massively Parallel GPU-based Code for Ray Tracing in Relativistic Spacetimes
NASA Astrophysics Data System (ADS)
Chan, Chi-kwan; Psaltis, Dimitrios; Özel, Feryal
2013-11-01
We introduce GRay, a massively parallel integrator designed to trace the trajectories of billions of photons in a curved spacetime. This graphics-processing-unit (GPU)-based integrator employs the stream processing paradigm, is implemented in CUDA C/C++, and runs on nVidia graphics cards. The peak performance of GRay using single-precision floating-point arithmetic on a single GPU exceeds 300 GFLOP (or 1 ns per photon per time step). For a realistic problem, where the peak performance cannot be reached, GRay is two orders of magnitude faster than existing central-processing-unit-based ray-tracing codes. This performance enhancement allows more effective searches of large parameter spaces when comparing theoretical predictions of images, spectra, and light curves from the vicinities of compact objects to observations. GRay can also perform on-the-fly ray tracing within general relativistic magnetohydrodynamic algorithms that simulate accretion flows around compact objects. Making use of this algorithm, we calculate the properties of the shadows of Kerr black holes and the photon rings that surround them. We also provide accurate fitting formulae of their dependencies on black hole spin and observer inclination, which can be used to interpret upcoming observations of the black holes at the center of the Milky Way, as well as M87, with the Event Horizon Telescope.
GRay: A MASSIVELY PARALLEL GPU-BASED CODE FOR RAY TRACING IN RELATIVISTIC SPACETIMES
Chan, Chi-kwan; Psaltis, Dimitrios; Özel, Feryal
2013-11-01
We introduce GRay, a massively parallel integrator designed to trace the trajectories of billions of photons in a curved spacetime. This graphics-processing-unit (GPU)-based integrator employs the stream processing paradigm, is implemented in CUDA C/C++, and runs on nVidia graphics cards. The peak performance of GRay using single-precision floating-point arithmetic on a single GPU exceeds 300 GFLOP (or 1 ns per photon per time step). For a realistic problem, where the peak performance cannot be reached, GRay is two orders of magnitude faster than existing central-processing-unit-based ray-tracing codes. This performance enhancement allows more effective searches of large parameter spaces when comparing theoretical predictions of images, spectra, and light curves from the vicinities of compact objects to observations. GRay can also perform on-the-fly ray tracing within general relativistic magnetohydrodynamic algorithms that simulate accretion flows around compact objects. Making use of this algorithm, we calculate the properties of the shadows of Kerr black holes and the photon rings that surround them. We also provide accurate fitting formulae of their dependencies on black hole spin and observer inclination, which can be used to interpret upcoming observations of the black holes at the center of the Milky Way, as well as M87, with the Event Horizon Telescope.
Ray-tracing simulation and SABER satellite observations of convective gravity waves
NASA Astrophysics Data System (ADS)
Kalisch, Silvio; Eckermann, Stephen; Ern, Manfred; Preusse, Peter; Riese, Martin; Trinh, Quang Thai; Kim, Young-Ha; Chun, Hye-Yeong
Gravity waves (GWs) are known as a coupling mechanism between different atmospheric layers. They contribute to the wave-driving of the QBO and are also responsible for driving large scale circulations like the Brewer-Dobson circulation. One major and highly variable source of GWs is convection. Deep convection in the tropics excites GWs with prominent amplitudes and horizontal phase speeds of up to 90 m/s. These GWs propagate upward and, when breaking, release the wave's momentum, thus accelerate the background flow. Direction and magnitude of the acceleration strongly depends on wind filtering between the convective GW source and the considered altitude. Both, the generation mechanism of GWs close to the top of deep convective towers and the wind filtering process during GW propagation largely influence the GW spectrum found in the tropical middle atmosphere and therefore magnitude and direction of the acceleration. We present the results of GW ray-tracing calculations from tropospheric (convective) sources up to the mesosphere. The Gravity wave Regional Or Global RAy-Tracer (GROGRAT) was used to perform the GW trajectory calculations. The convective GW source scheme from Yonsei University (South Korea) served as the lower boundary condition to quantify the GW excitation from convection. Heating rates, cloud top data, and atmospheric background data were provided by the MERRA dataset for the calculation of convective forcing from deep convection and for the atmospheric background of the ray-tracing calculations afterwards. In order to validate our ray-tracing simulation results, we compare them to satellite measurements of temperature amplitudes and momentum fluxes from the SABER instrument. Therefore, observational constrains from limb-sounding instruments have been quantified. Influences of orbit geometry, the instrument's observational filter, and the wavelength shift in the observed GW spectrum are discussed. Geographic structures in the observed global
NASA Astrophysics Data System (ADS)
Boudousq, V.; Bordy, T.; Gonon, G.; Dinten, J. M.
2005-04-01
The LEXXOS (DMS, Montpellier, France) is the first axial and total body cone beam bone densitometer using a 2D digital radiographic detector. Technical principles and performances for BMD measurements have been presented in previous papers. Bone densitometers are also used on small animals for drug development. In this paper, we show how the LEXXOS system can be adapted to small animals examinations, and its performances are evaluated. At first, in order to take advantage of the whole area of the digital flat panel X-ray detector, the geometrical configuration has been adapted. Secondly, as small animals present low BMD, a specific dual energy calibration has been defined. This adapted system has then been evaluated on two sets of mice: six reference mice and six ovariectomized mice. Each month, these two populations have been examined and the total body BMD has been measured. This evaluation has shown that the right order of BMD magnitude has been obtained and, as expected, BMD increases on the two sets until age of puberty and after this period, decreases significantly for the ovariectomized set. Moreover, the bone image obtained by dual energy processing on LEXXOS presents a radiographic image quality providing with useful complementary information on bone morphometry and architecture.
Ashley, Carlee E; Dunphy, Darren R; Jiang, Zhang; Carnes, Eric C; Yuan, Zhen; Petsev, Dimiter N; Atanassov, Plamen B; Velev, Orlin D; Sprung, Michael; Wang, Jin; Peabody, David S; Brinker, C Jeffrey
2011-04-18
The rapid assembly of icosohedral virus-like particles (VLPs) into highly ordered (domain size > 600 nm), oriented 2D superlattices directly onto a solid substrate using convective coating is demonstrated. In-situ grazing-incidence small-angle X-ray scattering (GISAXS) is used to follow the self-assembly process in real time to characterize the mechanism of superlattice formation, with the ultimate goal of tailoring film deposition conditions to optimize long-range order. From water, GISAXS data are consistent with a transport-limited assembly process where convective flow directs assembly of VLPs into a lattice oriented with respect to the water drying line. Addition of a nonvolatile solvent (glycerol) modified this assembly pathway, resulting in non-oriented superlattices with improved long-range order. Modification of electrostatic conditions (solution ionic strength, substrate charge) also alters assembly behavior; however, a comparison of in-situ assembly data between VLPs derived from the bacteriophages MS2 and Qβ show that this assembly process is not fully described by a simple Derjaguin-Landau-Verwey-Overbeek model alone. PMID:21425464
X-ray study of femtosecond structural dynamics in the 2D charge density wave compound 1T-TaS2
NASA Astrophysics Data System (ADS)
Laulhé, C.; Cario, L.; Corraze, B.; Janod, E.; Huber, T.; Lantz, G.; Boulfaat, S.; Ferrer, A.; Mariager, S. O.; Johnson, J. A.; Grübel, S.; Lübcke, A.; Ingold, G.; Beaud, P.; Johnson, S. L.; Ravy, S.
2015-03-01
1T-TaS2 is a 2D metallic compound which undergoes a series of electronically driven phase transitions toward charge density wave and Mott phases. Its intricate electron-phonon interactions and electron-electron correlations have been promising peculiar out-of-equilibrium dynamics. In this paper, we provide the first direct information on the atomic structure response to an ultra-fast infrared laser pulse in the commensurate phase of 1T-TaS2, by using femtosecond time-resolved X-ray diffraction. We show that ultra-fast excitation with near-infrared photons drives a displacive excitation of the amplitude mode of the commensurate charge density wave. About 3 ps after laser excitation, the system reaches a new, photo-induced state that is maintained for at least 10 ps. We give evidence that this long-lived state exhibits the same structural modulation as in the thermodynamically stable commensurate phase, with a large correlation length. Only the average amplitude of the modulation is found to decrease. We propose that the long-lived state is formed from the commensurate phase by reducing the modulation amplitude on few superlattice nodes. The underlying mechanism proposed is the annihilation of self-trapped polarons.
An Efficient Ray-Tracing Method for Determining Terrain Intercepts in EDL Simulations
NASA Technical Reports Server (NTRS)
Shidner, Jeremy D.
2016-01-01
The calculation of a ray's intercept from an arbitrary point in space to a prescribed surface is a common task in computer simulations. The arbitrary point often represents an object that is moving according to the simulation, while the prescribed surface is fixed in a defined frame. For detailed simulations, this surface becomes complex, taking the form of real-world objects such as mountains, craters or valleys which require more advanced methods to accurately calculate a ray's intercept location. Incorporation of these complex surfaces has commonly been implemented in graphics systems that utilize highly optimized graphics processing units to analyze such features. This paper proposes a simplified method that does not require computationally intensive graphics solutions, but rather an optimized ray-tracing method for an assumed terrain dataset. This approach was developed for the Mars Science Laboratory mission which landed on the complex terrain of Gale Crater. First, this paper begins with a discussion of the simulation used to implement the model and the applicability of finding surface intercepts with respect to atmosphere modeling, altitude determination, radar modeling, and contact forces influencing vehicle dynamics. Next, the derivation and assumptions of the intercept finding method are presented. Key assumptions are noted making the routines specific to only certain types of surface data sets that are equidistantly spaced in longitude and latitude. The derivation of the method relies on ray-tracing, requiring discussion on the formulation of the ray with respect to the terrain datasets. Further discussion includes techniques for ray initialization in order to optimize the intercept search. Then, the model implementation for various new applications in the simulation are demonstrated. Finally, a validation of the accuracy is presented along with the corresponding data sets used in the validation. A performance summary of the method will be shown using
A comprehensive ray tracing study on the impact of solar reflections from glass curtain walls.
Wong, Justin S J
2016-01-01
To facilitate the investigation of the impact of solar reflection from the façades of skyscrapers to surrounding environment, a comprehensive ray tracing model has been developed using the International Commerce Centre (ICC) in Hong Kong as an example. Taking into account the actual physical dimensions of buildings and meteorological data, the model simulates and traces the paths of solar reflections from ICC to the surrounding buildings, assessing the impact in terms of hit locations, light intensity and the hit time on each day throughout the year. Our analyses show that various design and architectural features of ICC have amplified the intensity of reflected solar rays and increased the hit rates of surrounding buildings. These factors include the high reflectivity of glass panels, their upward tilting angles, the concave profile of the 'Dragon Tail' (glass panels near the base), the particular location and orientation of ICC, as well as the immense height of ICC with its large reflective surfaces. The simulation results allow us to accurately map the date and time when the ray projections occur on each of the target buildings, rendering important information such as the number of converging (overlapping) projections, and the actual light intensity hitting each of the buildings at any given time. Comparisons with other skyscrapers such as Taipei 101 in Taiwan and 2-IFC (International Finance Centre) Hong Kong are made. Remedial actions for ICC and preventive measures are also discussed. PMID:26646546
NASA Astrophysics Data System (ADS)
Knowlton, S. F.; Hartwell, G. J.; Maurer, D. A.; Marushchenko, N. B.; Turkin, Y.; Bigelow, T.
2015-11-01
Plasmas in the Compact Toroidal Hybrid (CTH), a five field period, l = 2 torsatron (B0 = 0 . 5 T R0 = 0 . 75 m, ap ~ 0 . 2 m) will be heated by second harmonic X-mode electron cyclotron heating with power provided by a 28 GHz gyrotron capable of producing up to 200 kW. Ray-tracing calculations that will guide the selection of the launching position, antenna focal length, and beam-steering characteristics are performed with the TRAVIS code. Non-axisymmetric vacuum and current-carrying CTH equilibria for the ray tracing are modeled with the V3FIT code. The calculated absorption is highest for vertically propagating rays that traverse the region where a saddle of resonant field strength exists. However, the absorption for top-launched waves is more sensitive to variations in the magnetic equilibria than for a radial side launch where the magnetic field profile is tokamak-like. This work is supported by U.S. Department of Energy Grant No. DE-FG02-00ER54610.
Photorealistic ray tracing of free-space invisibility cloaks made of uniaxial dielectrics.
Halimeh, Jad C; Wegener, Martin
2012-12-17
The design rules of transformation optics generally lead to spatially inhomogeneous and anisotropic impedance-matched magneto-dielectric material distributions for, e.g., free-space invisibility cloaks. Recently, simplified anisotropic non-magnetic free-space cloaks made of a locally uniaxial dielectric material (calcite) have been realized experimentally. In a two-dimensional setting and for in-plane polarized light propagating in this plane, the cloaking performance can still be perfect for light rays. However, for general views in three dimensions, various imperfections are expected. In this paper, we study two different purely dielectric uniaxial cylindrical free-space cloaks. For one, the optic axis is along the radial direction, for the other one it is along the azimuthal direction. The azimuthal uniaxial cloak has not been suggested previously to the best of our knowledge. We visualize the cloaking performance of both by calculating photorealistic images rendered by ray tracing. Following and complementing our previous ray-tracing work, we use an equation of motion directly derived from Fermat's principle. The rendered images generally exhibit significant imperfections. This includes the obvious fact that cloaking does not work at all for horizontal or for ordinary linear polarization of light. Moreover, more subtle effects occur such as viewing-angle-dependent aberrations. However, we still find amazingly good cloaking performance for the purely dielectric azimuthal uniaxial cloak. PMID:23263067
Spin tracking simulations in AGS based on ray-tracing methods - bare lattice, no snakes -
Meot, F.; Ahrens, L.; Gleen, J.; Huang, H.; Luccio, A.; MacKay, W. W.; Roser, T.; Tsoupas, N.
2009-09-01
This Note reports on the first simulations of and spin dynamics in the AGS using the ray-tracing code Zgoubi. It includes lattice analysis, comparisons with MAD, DA tracking, numerical calculation of depolarizing resonance strengths and comparisons with analytical models, etc. It also includes details on the setting-up of Zgoubi input data files and on the various numerical methods of concern in and available from Zgoubi. Simulations of crossing and neighboring of spin resonances in AGS ring, bare lattice, without snake, have been performed, in order to assess the capabilities of Zgoubi in that matter, and are reported here. This yields a rather long document. The two main reasons for that are, on the one hand the desire of an extended investigation of the energy span, and on the other hand a thorough comparison of Zgoubi results with analytical models as the 'thin lens' approximation, the weak resonance approximation, and the static case. Section 2 details the working hypothesis : AGS lattice data, formulae used for deriving various resonance related quantities from the ray-tracing based 'numerical experiments', etc. Section 3 gives inventories of the intrinsic and imperfection resonances together with, in a number of cases, the strengths derived from the ray-tracing. Section 4 gives the details of the numerical simulations of resonance crossing, including behavior of various quantities (closed orbit, synchrotron motion, etc.) aimed at controlling that the conditions of particle and spin motions are correct. In a similar manner Section 5 gives the details of the numerical simulations of spin motion in the static case: fixed energy in the neighboring of the resonance. In Section 6, weak resonances are explored, Zgoubi results are compared with the Fresnel integrals model. Section 7 shows the computation of the {rvec n} vector in the AGS lattice and tuning considered. Many details on the numerical conditions as data files etc. are given in the Appendix Section
Fitting of NWM Ray-traced Slant Factors to Closed-form Tropospheric Mapping Functions
NASA Astrophysics Data System (ADS)
Urquhart, L.; Nievinski, F. G.; Santos, M. C.
2009-05-01
Ray-tracing in numerical weather models (NWM) is a promising solution for describing the elevation angle- and azimuth-dependence of tropospheric delay, especially at very low elevation angles, in an attempt to de- correlate vertical position and zenith tropospheric delay during GPS estimation. On the other hand, mapping functions expressed in closed form remain imperative, demanded by the need for (i) fast processing and (ii) convenient distribution to end-users, who employ a variety of third-party GPS processing packages. We investigate the fitting of ray-tracing results to closed-form expressions. We neglect the variation of the tropospheric delay with latitude, longitude, and height, offering a mapping function valid for a specific station site (similarly as done for VMF1-Site [Boehm et al., 1996]). We focus on the variation of the delay with time, elevation angle, and azimuth. For the time-dependence, we choose to work with slant factors instead of slant delays, because the former are more stable in time than the latter; that is a consequence of the normalization by zenith delays which removes the bulk of the variation with time. For the elevation angle-dependence we compare the continued form fraction of Yan and Ping [1995] with that of Marini [1972] (normalized to yield unity at zenith, as given by Herring [1992]). The latter is more commonly used, but the former is expected to provide a better fit at elevation angles below five degrees. Since the ray-tracing results do not necessarily assume azimuthal symmetry, we have to account for the azimuth-dependence. For that we compare the single-direction model of Davis et al. [1993] with the inclusion of secondary directions [Seko et al., 2004] and arbitrary spherical harmonics [Böhm and Schuh, 2001]). We also assess whether physically-oblivious models (i.e., not derived from analytical idealized atmospheric models), such spline or polynomials, as suggested by Rocken et al. [2001], are adequate.
Modeling and analysis of novel laser weld joint designs using optical ray tracing.
Milewski, J. O.
2002-01-01
Reflection of laser energy presents challenges in material processing that can lead to process inefficiency or process instability. Understanding the fundamentals of non-imaging optics and the reflective propagation of laser energy can allow process and weld joint designs to take advantage of these reflections to enhance process efficiency or mitigate detrimental effects. Optical ray tracing may be used within a 3D computer model to evaluate novel joint and fixture designs for laser welding that take advantage of the reflective propagation of laser energy. This modeling work extends that of previous studies by the author and provides comparison with experimental studies performed on highly reflective metals. Practical examples are discussed.
NASA Technical Reports Server (NTRS)
Gardner, C. S.; Rowlett, J. R.; Hendrickson, B. E.
1978-01-01
Errors may be introduced in satellite laser ranging data by atmospheric refractivity. Ray tracing data have indicated that horizontal refractivity gradients may introduce nearly 3-cm rms error when satellites are near 10-degree elevation. A correction formula to compensate for the horizontal gradients has been developed. Its accuracy is evaluated by comparing it to refractivity profiles. It is found that if both spherical and gradient correction formulas are employed in conjunction with meteorological measurements, a range resolution of one cm or less is feasible for satellite elevation angles above 10 degrees.
X-Ray fluorescence analysis of trace elements in fruit juice
NASA Astrophysics Data System (ADS)
Bao, Sheng-Xiang; Wang, Zhi-Hong; Liu, Jing-Song
1999-12-01
X-Ray fluorescence spectrometry is applied to the determination of trace elements in fruit juice characterized by a high content of sugar and other soluble solid substances. Samples are prepared by evaporation, carbonization and pressing into discs. The synthesis of standards is described in detail. All element concentrations are directly estimated from linear calibration curves obtained without any matrix correction. The results of the analysis are in good agreement with those given by inductively coupled plasma-atomic emission spectrometry and atomic absorption spectrometry techniques.
Infrasonic ray tracing applied to mesoscale atmospheric structures: refraction by hurricanes.
Bedard, Alfred J; Jones, R Michael
2013-11-01
A ray-tracing program is used to estimate the refraction of infrasound by the temperature structure of the atmosphere and by hurricanes represented by a Rankine-combined vortex wind plus a temperature perturbation. Refraction by the hurricane winds is significant, giving rise to regions of focusing, defocusing, and virtual sources. The refraction of infrasound by the temperature anomaly associated with a hurricane is small, probably no larger than that from uncertainties in the wind field. The results are pertinent to interpreting ocean wave generated infrasound in the vicinities of tropical cyclones. PMID:24180755