Toward optimal X-ray flux utilization in breast CT

A realistic computer-simulation of a breast computed tomography (CT) system and subject is constructed. The model is used to investigate the optimal number of views for the scan given a fixed total X-ray fluence. The reconstruction algorithm is based on accurate solution to a constrained, TV-minimization problem, which has received much interest recently for sparse-view CT data.Comment: accepted to the 11th International Meeting on Fully Three-Dimensional Image Reconstruction in Radiology and Nuclear Medicine 201

Enhanced imaging of microcalcifications in digital breast tomosynthesis through improved image-reconstruction algorithms

PURPOSE: We develop a practical, iterative algorithm for image-reconstruction in under-sampled tomographic systems, such as digital breast tomosynthesis (DBT). METHOD: The algorithm controls image regularity by minimizing the image total $p$-variation (TpV), a function that reduces to the total variation when $p=1.0$ or the image roughness when $p=2.0$. Constraints on the image, such as image positivity and estimated projection-data tolerance, are enforced by projection onto convex sets (POCS). The fact that the tomographic system is under-sampled translates to the mathematical property that many widely varied resultant volumes may correspond to a given data tolerance. Thus the application of image regularity serves two purposes: (1) reduction of the number of resultant volumes out of those allowed by fixing the data tolerance, finding the minimum image TpV for fixed data tolerance, and (2) traditional regularization, sacrificing data fidelity for higher image regularity. The present algorithm allows for this dual role of image regularity in under-sampled tomography. RESULTS: The proposed image-reconstruction algorithm is applied to three clinical DBT data sets. The DBT cases include one with microcalcifications and two with masses. CONCLUSION: Results indicate that there may be a substantial advantage in using the present image-reconstruction algorithm for microcalcification imaging.Comment: Submitted to Medical Physic

Simulation of an 1857-like Mw 7.9 San Andreas Fault Earthquake and the Response of Tall Steel Moment Frame Buildings in Southern California – A Prototype Study

In 1857, an earthquake of magnitude 7.9 occurred on the San Andreas fault, starting at Parkfield and rupturing in a southeasterly direction for more than 360 km. Such a unilateral rupture produces significant directivity toward the San Fernando and Los Angeles basins. The strong shaking in the basins due to this earthquake would have had significant long-period content (2-8 s), and the objective of this study is to quantify the impact of such an earthquake on two 18-story steel moment frame building models, hypothetically located at 636 sites on a 3.5 km grid in southern California. End-to-end simulations include modeling the source and rupture of a fault at one end, numerically propagating the seismic waves through the earth structure, simulating the damage to engineered structures and estimating the economic impact at the other end using high-performance computing. In this prototype study, we use an inferred finite source model of the magnitude 7.9, 2002 Denali fault earthquake in Alaska, and map it onto the San Andreas fault with the rupture originating at Parkfield and propagating southward over a distance of 290 km. Using the spectral element seismic wave propagation code, SPECFEM3D, we simulate an 1857-like earthquake on the San Andreas fault and compute ground motions at the 636 analysis sites. Using the nonlinear structural analysis program, FRAME3D, we subsequently analyze 3-D structural models of an existing tall steel building designed using the 1982 Uniform Building Code (UBC), as well as one designed according to the 1997 UBC, subjected to the computed ground motion at each of these sites. We summarize the performance of these structural models on contour maps of peak interstory drift. We then perform an economic loss analysis for the two buildings at each site, using the Matlab Damage and Loss Analysis (MDLA) toolbox developed to implement the PEER loss-estimation methodology. The toolbox includes damage prediction and repair cost estimation for structural and non-structural components and allows for the computation of the mean and variance of building repair costs conditional on engineering demand parameters (i.e. inter-story drift ratios and peak floor accelerations). Here, we modify it to treat steel-frame high-rises, including aspects such as mechanical, electrical and plumbing systems, traction elevators, and the possibility of irreparable structural damage. We then generate contour plots of conditional mean losses for the San Fernando and the Los Angeles basins for the pre-Northridge and modern code-designed buildings, allowing for comparison of the economic effects of the updated code for the scenario event. In principle, by simulating multiple seismic events, consistent with the probabilistic seismic hazard for a building site, the same basic approach could be used to quantify the uncertain losses from future earthquakes

Age-dependent differences in human brain activity using a face- and location-matching task: An fMRI study

Purpose: To evaluate the differences of cortical activation patterns in young and elderly healthy subjects for object and spatial visual processing using a face- and location-matching task. Materials and Methods: We performed a face- and a location-matching task in 15 young (mean age: 28 +/- 9 years) and 19 elderly (mean age: 71 +/- 6 years) subjects. Each experiment consisted of 7 blocks alternating between activation and control condition. For face matching, the subjects had to indicate whether two displayed faces were identical or different. For location matching, the subjects had to press a button whenever two objects had an identical position. For control condition, we used a perception task with abstract images. Functional imaging was performed on a 1.5-tesla scanner using an EPI sequence. Results: In the face-matching task, the young subjects showed bilateral (right 1 left) activation in the occipito-temporal pathway (occipital gyrus, inferior and middle temporal gyrus). Predominantly right hemispheric activations were found in the fusiform gyrus, the right dorsolateral prefrontal cortex (inferior and middle frontal gyrus) and the superior parietal gyrus. In the elderly subjects, the activated areas in the right fronto-lateral cortex increased. An additional activated area could be found in the medial frontal gyrus (right > left). In the location-matching task, young subjects presented increased bilateral (right > left) activation in the superior parietal lobe and precuneus compared with face matching. The activations in the occipito-temporal pathway, in the right fronto-lateral cortex and the fusiform gyrus were similar to the activations found in the face-matching task. In the elderly subjects, we detected similar activation patterns compared to the young subjects with additional activations in the medial frontal gyrus. Conclusion: Activation patterns for object-based and spatial visual processing were established in the young and elderly healthy subjects. Differences between the elderly and young subjects could be evaluated, especially by using a face-matching task. Copyright (c) 2007 S. Karger AG, Basel

Transverse effects in multifrequency Raman generation

The theory of ultrabroadband multifrequency Raman generation is extended, for the first time, to allow for beam-propagation effects in one and two transverse dimensions. We show that a complex transverse structure develops even when diffraction is neglected. In the general case, we examine how the ultrabroadband multifrequency Raman generation process is affected by the intensity, phase quality, and width of the input beams, and by the length of the Raman medium. The evolution of power spectra, intensity profiles, and global characteristics of the multifrequency beams are investigated and explained. In the two-dimensional transverse case, bandwidths comparable to the optical carrier frequency, spanning the whole visible spectrum and beyond, are still achievable

Chaotic Orbits in Thermal-Equilibrium Beams: Existence and Dynamical Implications

Phase mixing of chaotic orbits exponentially distributes these orbits through their accessible phase space. This phenomenon, commonly called ``chaotic mixing'', stands in marked contrast to phase mixing of regular orbits which proceeds as a power law in time. It is operationally irreversible; hence, its associated e-folding time scale sets a condition on any process envisioned for emittance compensation. A key question is whether beams can support chaotic orbits, and if so, under what conditions? We numerically investigate the parameter space of three-dimensional thermal-equilibrium beams with space charge, confined by linear external focusing forces, to determine whether the associated potentials support chaotic orbits. We find that a large subset of the parameter space does support chaos and, in turn, chaotic mixing. Details and implications are enumerated.Comment: 39 pages, including 14 figure

Ultrabroad-bandwidth multifrequency Raman generation

We report on the modeling of transient stimulated rotational Raman scattering in H2 gas. We predict a multifrequency output, spanning a bandwidth greater than the pump frequency, that may be generated without any significant delay with respect to the pump pulses. The roles of dispersion and transiency are quantified

Production of Enhanced Beam Halos via Collective Modes and Colored Noise

We investigate how collective modes and colored noise conspire to produce a beam halo with much larger amplitude than could be generated by either phenomenon separately. The collective modes are lowest-order radial eigenmodes calculated self-consistently for a configuration corresponding to a direct-current, cylindrically symmetric, warm-fluid Kapchinskij-Vladimirskij equilibrium. The colored noise arises from unavoidable machine errors and influences the internal space-charge force. Its presence quickly launches statistically rare particles to ever-growing amplitudes by continually kicking them back into phase with the collective-mode oscillations. The halo amplitude is essentially the same for purely radial orbits as for orbits that are initially purely azimuthal; orbital angular momentum has no statistically significant impact. Factors that do have an impact include the amplitudes of the collective modes and the strength and autocorrelation time of the colored noise. The underlying dynamics ensues because the noise breaks the Kolmogorov-Arnol'd-Moser tori that otherwise would confine the beam. These tori are fragile; even very weak noise will eventually break them, though the time scale for their disintegration depends on the noise strength. Both collective modes and noise are therefore centrally important to the dynamics of halo formation in real beams.Comment: For full resolution pictures please go to http://www.nicadd.niu.edu/research/beams

Efficient computation of matched solutions of the Kapchinskij-Vladimirskij envelope equations for periodic focusing lattices

A new iterative method is developed to numerically calculate the periodic, matched beam envelope solution of the coupled Kapchinskij-Vladimirskij (KV) equations describing the transverse evolution of a beam in a periodic, linear focusing lattice of arbitrary complexity. Implementation of the method is straightforward. It is highly convergent and can be applied to all usual parameterizations of the matched envelope solutions. The method is applicable to all classes of linear focusing lattices without skew couplings, and also applies to all physically achievable system parameters -- including where the matched beam envelope is strongly unstable. Example applications are presented for periodic solenoidal and quadrupole focusing lattices. Convergence properties are summarized over a wide range of system parameters.Comment: 20 pages, 5 figures, Mathematica source code provide