scispace - formally typeset
Search or ask a question

Showing papers on "Attenuation published in 2011"


Journal ArticleDOI
TL;DR: This first measurement of ultrasound absorption in bone can be used to estimate the amount of heat deposition based on knowledge of the acoustic field and it is demonstrated that only a small part of the attenuation is due to absorption inBone and that the majority of the dBs are due to reflection, scattering, and mode conversion.
Abstract: (Received 27 June 2011; revised 26 October 2011; accepted for publication 22 November 2011;published 21 December 2011)Purpose: Measured values of ultrasound attenuation in bone represent a combination of differentloss mechanisms. As a wave is transmitted from a fluid into bone, reflections occur at the interface.In the bone, mode conversion occurs between longitudinal and shear modes and the mechanicalwave is scattered by its complex internal microstructure. Finally, part of the wave energy isabsorbed by the bone and converted into heat. Due to the complexity of the wave propagation andthe difficulty in performing measurements that are capable of separating the various loss mecha-nisms, there are currently no estimates of the absorption in bone. The aim of this paper is, thus, toquantify the attenuation, scattering, and thermal absorption in bone.Methods: An attenuating model of wave propagation in bone is established and used to develop athree-dimensional finite difference time domain numerical algorithm. Hydrophone and optical het-erodyne interferometer measurements of the acoustic field as well as a x-ray microtomography ofthe bone sample are used to drive the simulations and to measure the attenuation. The acousticmeasurements are performed concurrently with an infrared camera that can measure the tempera-ture elevation during insonication. A link between the temperature and the absorption via a three-dimensional thermal simulation is then used to quantify the absorption coefficients for longitudinaland shear waves in cortical bone.Results: We demonstrate that only a small part of the attenuation is due to absorption in bone andthat the majority of the attenuation is due to reflection, scattering, and mode conversion. In the ninesamples of a human used for the study, the absorption time constant for cortical bone was deter-mined to be 1.04 ls628%. This corresponds to a longitudinal absorption of 2.7 dB/cm and a shearabsorption of 5.4 dB/cm. The experimentally measured attenuation across the approximately 8 mmthick samples was 13.360.97 dB/cm.Conclusions: This first measurement of ultrasound absorption in bone can be used to estimate theamount of heat deposition based on knowledge of the acoustic field.

258 citations


Journal ArticleDOI
TL;DR: In this article, the authors measure the monochromatic transmission matrix (TM) of an optical multiple scattering medium using a spatial light modulator together with a phase-shifting interferometry measurement method.
Abstract: We experimentally measure the monochromatic transmission matrix (TM) of an optical multiple scattering medium using a spatial light modulator together with a phase-shifting interferometry measurement method. The TM contains all the information needed to shape the scattered output field at will or to detect an image through the medium. We confront theory and experiment for these applications and study the effect of noise on the reconstruction method. We also extracted from the TM information about the statistical properties of the medium and the light transport within it. In particular, we are able to isolate the contributions of the memory effect and measure its attenuation length.

249 citations


Proceedings ArticleDOI
25 Jul 2011
TL;DR: Optimization of tomographic techniques for image synthesis on displays composed of compact volumes of light-attenuating material allows optimal construction of high dynamic range displays, confirming existing heuristics and providing the first extension to multiple, disjoint layers.
Abstract: We develop tomographic techniques for image synthesis on displays composed of compact volumes of light-attenuating material. Such volumetric attenuators recreate a 4D light field or high-contrast 2D image when illuminated by a uniform backlight. Since arbitrary oblique views may be inconsistent with any single attenuator, iterative tomographic reconstruction minimizes the difference between the emitted and target light fields, subject to physical constraints on attenuation. As multi-layer generalizations of conventional parallax barriers, such displays are shown, both by theory and experiment, to exceed the performance of existing dual-layer architectures. For 3D display, spatial resolution, depth of field, and brightness are increased, compared to parallax barriers. For a plane at a fixed depth, our optimization also allows optimal construction of high dynamic range displays, confirming existing heuristics and providing the first extension to multiple, disjoint layers. We conclude by demonstrating the benefits and limitations of attenuation-based light field displays using an inexpensive fabrication method: separating multiple printed transparencies with acrylic sheets.

228 citations


Journal ArticleDOI
TL;DR: Both the speed of sound and attenuation showed dependence on the skull density that varied with the frequency, and that these effects are different according to the type of bone in the skull (cortical and trabecular).
Abstract: For medical applications of ultrasound inside the brain, it is necessary to understand the relationship between the apparent density of skull bone and its corresponding speed of sound and attenuation coefficient. Although there have been previous studies exploring this phenomenon, there is still a need to extend the measurements to cover more of the clinically relevant frequency range. The results of measurements of the longitudinal speed of sound and attenuation coefficient are presented for specimens of human calvaria. The study was performed for the frequencies of 0.27, 0.836, 1.402, 1.965 and 2.525 MHz. Specimens were obtained from fresh cadavers through a protocol with the Division of Anatomy of the University of Toronto. The protocol was approved by the Research Ethics Board of Sunnybrook Health Sciences Centre. The specimens were mounted in polycarbonate supports that were marked for stereoscopic positioning. Computer tomography (CT) scans of the skulls mounted on their supports were performed, and a three-dimensional skull surface was reconstructed. This surface was used to guide a positioning system to ensure the normal sound incidence of an acoustic signal. This signal was produced by a focused device with a diameter of 5 cm and a focal length of 10 cm. Measurements of delay in time of flight were carried out using a needle hydrophone. Measurements of effective transmitted energy were carried out using a radiation force method with a 10 µg resolution scale. Preliminary functions of speed of sound and attenuation coefficient, both of which are related to apparent density, were established using a multi-layer propagation model that takes into account speed of sound, density and thickness of the layer. An optimization process was executed from a large set of random functions and the best functions were chosen for those ones that closest reproduced the experimental observations. The final functions were obtained after a second pass of the optimization process was executed, but this time using a finite-difference time-difference solution of the Westervelt equation, which is more precise than the multi-layer model but much more time consuming for computation. For six of seven specimens, measurements were carried out on five locations on the calvaria, and for the other specimen three measurements were made. In total, measurements were carried out on 33 locations. Results indicated the presence of dispersion effects and that these effects are different according to the type of bone in the skull (cortical and trabecular). Additionally, both the speed of sound and attenuation showed dependence on the skull density that varied with the frequency. Using the optimal functions and the information of density from the CT scans, the average values (±s.d.) of the speed of sound for cortical bone were estimated to be 2384(± 130), 2471(± 90), 2504(± 120), 2327(± 90) and 2053(± 40) m s(-1) for the frequencies of 270, 836, 1402, 1965 and 2526 kHz, respectively. For trabecular bone, and in the same order of frequency values, the speeds of sound were 2140(± 130), 2300(± 100), 2219(± 200), 2133(± 130) and 1937(± 40) m s(-1), respectively. The average values of the attenuation coefficient for cortical bone were 33(± 9), 240(± 9) and 307(± 30) Np m(-1) for the frequencies of 270, 836, and 1402, respectively. For trabecular bone, and in the same order of frequency values, the average values of the attenuation coefficient were 34(± 13), 216(± 16) and 375(± 30) Np m(-1), respectively. For frequencies of 1.965 and 2.525 MHz, no measurable radiation force was detected with the setup used.

213 citations


Journal ArticleDOI
TL;DR: In x-ray computed tomography (CT), more than two types of materials can be effectively separated with the use of an energy resolved photon-counting detector and classification methodology, and this applies to the case when contrast agents that contain K-absorption edges in the energy range of interest are present in the object.
Abstract: Purpose: The objective of the study was to demonstrate that, in x-ray computed tomography (CT), more than two types of materials can be effectively separated with the use of an energy resolved photon-counting detector and classification methodology. Specifically, this applies to the case when contrast agents that contain K-absorption edges in the energy range of interest are present in the object. This separation is enabled via the use of recently developed energy resolved photon-counting detectors with multiple thresholds, which allow simultaneous measurements of the x-ray attenuation at multiple energies. Methods: To demonstrate this capability, we performed simulations and physical experiments using a six-threshold energy resolved photon-counting detector. We imaged mouse-sized cylindrical phantoms filled with several soft-tissue-like and bone-like materials and with iodine-based and gadolinium-based contrast agents. The linear attenuation coefficients were reconstructed for each material in each energy window and were visualized as scatter plots between pairs of energy windows. For comparison, a dual-kVp CT was also simulated using the same phantom materials. In this case, the linear attenuation coefficients at the lower kVp were plotted against those at the higher kVp. Results: In both the simulations and the physical experiments, the contrast agents were easily separable from other soft-tissue-like and bone-like materials, thanks to the availability of the attenuation coefficient measurements at more than two energies provided by the energy resolved photon-counting detector. In the simulations, the amount of separation was observed to be proportional to the concentration of the contrast agents; however, this was not observed in the physical experiments due to limitations of the real detector system. We used the angle between pairs of attenuation coefficient vectors in either the 5-D space (for non-contrast-agent materials using energy resolved photon-counting acquisition) or a 2-D space (for contrast agents using energy resolved photon-counting acquisition and all materials using dual-kVp acquisition) as a measure of the degree of separation. Compared to dual-kVp techniques, an energy resolved detector provided a larger separation and the ability to separate different target materials using measurements acquired in different energy window pairs with a single x-ray exposure. Conclusions: We concluded that x-ray CT with an energy resolved photon-counting detector with more than two energy windows allows the separation of more than two types of materials, e.g., soft-tissue-like, bone-like, and one or more materials with K-edges in the energy range of interest. Separating material types using energy resolved photon-counting detectors has a number of advantages over dual-kVp CT in terms of the degree of separation and the number of materials that can be separated simultaneously.

184 citations


Journal ArticleDOI
TL;DR: In this paper, a general method for reducing white-beam X-ray total scattering raw data to the differential scattering cross section and pair distribution function is described, and a rigorous and revised Lorch function is proposed for removing the effects of Fourier truncation oscillations.
Abstract: A general method is described for reducing white-beam X-ray total scattering raw data to the differential scattering cross section and pair distribution function. The method incorporates corrections for X-ray fluorescence, Bremsstrahlung radiation, polarization, attenuation, multiple scattering and sample container scattering, and invokes the Krogh-Moe and Norman method to put the data on an absolute scale. An accurate method to convert the differential scattering cross section to the pair distribution function is also described, and a rigorous and revised Lorch function is proposed for removing the effects of Fourier truncation oscillations. The method can be equally applied to synchrotron X-ray data, where the data analysis can be simpler than at a laboratory source.

172 citations


Journal ArticleDOI
TL;DR: In this article, the total mass attenuation and linear attenuation coefficients, half-value thicknesses, effective atomic numbers, effective electron densities and atomic cross-sections at photons energies of 59.5 and 661 keV have been measured and calculated.

170 citations


Journal ArticleDOI
Andre Salomon1, A. Goedicke1, B Schweizer1, Til Aach1, Volkmar Schulz1 
TL;DR: A generic iterative reconstruction approach to simultaneously estimate the local tracer concentration and the attenuation distribution using the segmented MR image as anatomical reference, which indicates a robust and reliable alternative to other MR-AC approaches targeting patient specific quantitative analysis in time-of-flight PET/MR.
Abstract: Medical investigations targeting a quantitative analysis of the position emission tomography (PET) images require the incorporation of additional knowledge about the photon attenuation distribution in the patient. Today, energy range adapted attenuation maps derived from computer tomography (CT) scans are used to effectively compensate for image quality degrading effects, such as attenuation and scatter. Replacing CT by magnetic resonance (MR) is considered as the next evolutionary step in the field of hybrid imaging systems. However, unlike CT, MR does not measure the photon attenuation and thus does not provide an easy access to this valuable information. Hence, many research groups currently investigate different technologies for MR-based attenuation correction (MR-AC). Typically, these approaches are based on techniques such as special acquisition sequences (alone or in combination with subsequent image processing), anatomical atlas registration, or pattern recognition techniques using a data base of MR and corresponding CT images. We propose a generic iterative reconstruction approach to simultaneously estimate the local tracer concentration and the attenuation distribution using the segmented MR image as anatomical reference. Instead of applying predefined attenuation values to specific anatomical regions or tissue types, the gamma attenuation at 511 keV is determined from the PET emission data. In particular, our approach uses a maximum-likelihood estimation for the activity and a gradient-ascent based algorithm for the attenuation distribution. The adverse effects of scattered and accidental gamma coincidences on the quantitative accuracy of PET, as well as artifacts caused by the inherent crosstalk between activity and attenuation estimation are efficiently reduced using enhanced decay event localization provided by time-of-flight PET, accurate correction for accidental coincidences, and a reduced number of unknown attenuation coefficients. First results achieved with measured whole body PET data and reference segmentation from CT showed an absolute mean difference of 0.005 cm in the lungs, 0.0009 cm in case of fat, and 0.0015 cm for muscles and blood. The proposed method indicates a robust and reliable alternative to other MR-AC approaches targeting patient specific quantitative analysis in time-of-flight PET/MR.

165 citations


Journal ArticleDOI
TL;DR: The most accurate measurement of the attenuation of the atmosphere in the range from low frequencies up to 2 THz is reported, which indicates the potential of THz applications in imaging, sensing, ranging and communications.
Abstract: The attenuation of electromagnetic wave propagation in the clear atmosphere from low frequencies up to 2 THz is mainly caused by water vapor. Although there have been many numerical simulations and excellent early sub-mm and far-infrared measurements of this attenuation, there has remained controversy about the background absorption in the most promising windows of transparency below 1 THz. Here, we report an accurate terahertz time-domain spectroscopy (THz-TDS) characterization of water vapor from 0.2 to 2 THz. Our results agree with the previous predicted and measured attenuations for the weak water lines, but show more attenuation for the relatively transparent windows between these lines.

164 citations


Journal ArticleDOI
TL;DR: In this paper, the influence of diffuse scattering at high frequency on the characteristics of the communication channel and its implications on the non-line-of-sight propagation path is investigated.
Abstract: Recent years have seen a tremendous increase in the demand for wireless bandwidth. To support this demand by innovative and resourceful use of technology, future communication systems will have to shift towards higher carrier frequencies. Due to the tight regulatory situation, frequencies in the atmospheric attenuation window around 300 GHz appear very attractive to facilitate an indoor, short range, ultra high speed THz communication system. In this paper, we investigate the influence of diffuse scattering at such high frequencies on the characteristics of the communication channel and its implications on the non-line-of-sight propagation path. The Kirchhoff approach is verified by an experimental study of diffuse scattering from randomly rough surfaces commonly encountered in indoor environments using a fiber-coupled terahertz time-domain spectroscopy system to perform angle- and frequency-dependent measurements. Furthermore, we integrate the Kirchhoff approach into a self-developed ray tracing algorithm to model the signal coverage of a typical office scenario.

163 citations


Journal ArticleDOI
TL;DR: In this article, a quasi-static creep test with mesoscopic-scale heterogeneities is proposed to calculate the complex and frequency-dependent P wave moduli from the modeled stress-strain relations.
Abstract: [1] The finite element method is used to solve Biot's equations of consolidation in the displacement-pressure (u − p) formulation. We compute one-dimensional (1-D) and two-dimensional (2-D) numerical quasi-static creep tests with poroelastic media exhibiting mesoscopic-scale heterogeneities to calculate the complex and frequency-dependent P wave moduli from the modeled stress-strain relations. The P wave modulus is used to calculate the frequency-dependent attenuation (i.e., inverse of quality factor) and phase velocity of the medium. Attenuation and velocity dispersion are due to fluid flow induced by pressure differences between regions of different compressibilities, e.g., regions (or patches) saturated with different fluids (i.e., so-called patchy saturation). Comparison of our numerical results with analytical solutions demonstrates the accuracy and stability of the algorithm for a wide range of frequencies (six orders of magnitude). The algorithm employs variable time stepping and an unstructured mesh which make it efficient and accurate for 2-D simulations in media with heterogeneities of arbitrary geometries (e.g., curved shapes). We further numerically calculate the quality factor and phase velocity for 1-D layered patchy saturated porous media exhibiting random distributions of patch sizes. We show that the numerical results for the random distributions can be approximated using a volume average of White's analytical solution and the proposed averaging method is, therefore, suitable for a fast and transparent prediction of both quality factor and phase velocity. Application of our results to frequency-dependent reflection coefficients of hydrocarbon reservoirs indicates that attenuation due to wave-induced flow can increase the reflection coefficient at low frequencies, as is observed at some reservoirs.

Journal ArticleDOI
TL;DR: The ground motion prediction equations (GMPEs) developed as part of the Next Generation Attenuation of Ground Motions (NGA-West) project in 2008 are becoming widely used in seismic hazard analyses as mentioned in this paper.
Abstract: The ground-motion prediction equations (GMPEs) developed as part of the Next Generation Attenuation of Ground Motions (NGA-West) project in 2008 are becoming widely used in seismic hazard analyses. However, these new models are considerably more complicated than previous GMPEs, and they require several more input parameters. When employing the NGA models, users routinely face situations in which some of the required input parameters are unknown. In this paper, we present a framework for estimating the unknown source, path, and site parameters when implementing the NGA models in engineering practice, and we derive geometrically-based equations relating the three distance measures found in the NGA models. Our intent is for the content of this paper not only to make the NGA models more accessible, but also to help with the implementation of other present or future GMPEs. [DOI: 10.1193/1.3650372]

Patent
10 Nov 2011
TL;DR: In this paper, the SLM imparts a programmable pattern of attenuation that may be used to correct for asymmetries between the first and second modes of illumination or imaging.
Abstract: Methods are disclosed for measuring target structures formed by a lithographic process on a substrate. A grating structure within the target is smaller than an illumination spot and field of view of a measurement optical system. The optical system has a first branch leading to a pupil plane imaging sensor and a second branch leading to a substrate plane imaging sensor. A spatial light modulator is arranged in an intermediate pupil plane of the second branch of the optical system. The SLM imparts a programmable pattern of attenuation that may be used to correct for asymmetries between the first and second modes of illumination or imaging. By use of specific target designs and machine-learning processes, the attenuation patterns may also be programmed to act as filter functions, enhancing sensitivity to specific parameters of interest, such as focus.

Journal ArticleDOI
TL;DR: The demonstrated performance of the spatial-convolution method shows it is a powerful tool for reducing reconstruction artifacts originating from the detector finite size and improving the quality of optoacoustic reconstructions.
Abstract: Purpose: Optoacousticimaging enables mapping the optical absorption of biological tissue using optical excitation and acoustic detection. Although most image-reconstruction algorithms are based on the assumption of a detector with an isotropic sensitivity, the geometry of the detector often leads to a response with spatially dependent magnitude and bandwidth. This effect may lead to attenuation or distortion in the recorded signal and, consequently, in the reconstructed image. Methods: Herein, an accurate numerical method for simulating the spatially dependent response of an arbitrary-shape acoustic transducer is presented. The method is based on an analytical solution obtained for a two-dimensional line detector. The calculated response is incorporated in the forward model matrix of an optoacousticimaging setup using temporal convolution, and image reconstruction is performed by inverting the matrix relation. Results: The method was numerically and experimentally demonstrated in two dimensions for both flat and focused transducers and compared to the spatial-convolution method. In forward simulations, the developed method did not suffer from the numerical errors exhibited by the spatial-convolution method. In reconstruction simulations and experiments, the use of both temporal-convolution and spatial-convolution methods lead to an enhancement in resolution compared to a reconstruction with a point detectormodel. However, because of its higher modeling accuracy, the temporal-convolution method achieved a noise figure approximated three times lower than the spatial-convolution method. Conclusions: The demonstrated performance of the spatial-convolution method shows it is a powerful tool for reducing reconstruction artifacts originating from the detector finite size and improving the quality of optoacousticreconstructions. Furthermore, the method may be used for assessing new system designs. Specifically, detectors with nonstandard shapes may be investigated.

Journal ArticleDOI
M Guillot1, Luc Gingras1, Louis Archambault1, Sam Beddar, Luc Beaulieu1 
TL;DR: The spectral method can be used to accurately correct the Cerenkov light effect in plastic scintillation detectors (PSDs) and two new procedures to determine the calibration factors of PSDs are proposed.
Abstract: Purpose: The purposes of this work were: (1) To determine if a spectral method can accurately correct the Cerenkov light effect in plastic scintillation detectors (PSDs) for situations where the Cerenkov light is dominant over the scintillation light and (2) to develop a procedural guideline for accurately determining the calibration factors of PSDs. Methods: The authors demonstrate, by using the equations of the spectral method, that the condition for accurately correcting the effect of Cerenkov light is that the ratio of the two calibration factors must be equal to the ratio of the Cerenkov light measured within the two different spectral regions used for analysis. Based on this proof, the authors propose two new procedures to determine the calibration factors of PSDs, which were designed to respect this condition. A PSD that consists of a cylindrical polystyrene scintillating fiber ( 1.6 mm 3 ) coupled to a plastic optical fiber was calibrated by using these new procedures and the two reference procedures described in the literature. To validate the extracted calibration factors, relative dose profiles and output factors for a 6 MV photon beam from a medical linac were measured with the PSD and an ionization chamber. Emphasis was placed on situations where the Cerenkov light is dominant over the scintillation light and on situations dissimilar to the calibration conditions. Results: The authors found that the accuracy of the spectral method depends on the procedure used to determine the calibration factors of the PSD and on the attenuation properties of the optical fiber used. The results from the relative dose profile measurements showed that the spectral method can correct the Cerenkov light effect with an accuracy level of 1%. The results obtained also indicate that PSDsmeasure output factors that are lower than those measured with ionization chambers for square field sizes larger than 25 × 25 cm 2 , in general agreement with previously published Monte Carlo results. Conclusions: The authors conclude that the spectral method can be used to accurately correct the Cerenkov light effect in PSDs. The authors confirmed the importance of maximizing the difference of Cerenkov light production between calibrationmeasurements. The authors also found that the attenuation of the optical fiber, which is assumed to be constant in the original formulation of the spectral method, may cause a variation of the calibration factors in some experimental setups.

Journal ArticleDOI
TL;DR: In this paper, the applicability of the x-ray attenuation method to monitor corrosion products as well as the formation and propagation of cracks in cementitious materials, reinforced mortar samples were tested under accelerated corrosion conditions.

Journal ArticleDOI
TL;DR: In this paper, the authors presented some theoretical simulations of the specular scattering coefficient of a forested area, with the aim of demonstrating the potentiality of GNSS-R in monitoring forest biomass.

Journal ArticleDOI
TL;DR: A correction algorithm for accounting for attenuation effects is proposed, applicable in both the time and frequency domains, and which part of the optoacoustic signal spectrum is practically affected by those effects in realistic imaging scenarios is investigated.
Abstract: In this paper, it is demonstrated that the effects of acoustic attenuation may play a significant role in establishing the quality of tomographic optoacoustic reconstructions. Accordingly, spatially dependent reduction of signal amplitude leads to quantification errors in the reconstructed distribution of the optical absorption coefficient while signal broadening causes loss of image resolution. Here we propose a correction algorithm for accounting for attenuation effects, which is applicable in both the time and frequency domains. It is further investigated which part of the optoacoustic signal spectrum is practically affected by those effects in realistic imaging scenarios. The validity and benefits of the suggested modelling and correction approaches are experimentally validated in phantom measurements.

Journal ArticleDOI
01 Jan 2011
TL;DR: In this article, the bulk viscosity of millipore water over a temperature range of 7 to 50?C. was measured and the results showed that the volume of the volume is approximately three times larger than its shear counterpart.
Abstract: Despite its fundamental role in the dynamics of compressible fluids, bulk viscosity has received little experimental attention and there remains a paucity of measured data. Acoustic spectroscopy provides a robust and accurate approach to measuring this parameter. Working from the Navier-Stokes model of a compressible fluid one can show that the bulk viscosity makes a significant and measurable contribution to the frequency-squared acoustic attenuation. Here we employ this methodology to determine the bulk viscosity of Millipore water over a temperature range of 7 to 50?C. The measured attenuation spectra are consistent with the theoretical predictions, while the bulk viscosity of water is found to be approximately three times larger than its shear counterpart, reinforcing its significance in acoustic propagation. Moreover, our results demonstrate that this technique can be readily and generally applied to fluids to accurately determine their temperature dependent bulk viscosities.

Proceedings ArticleDOI
01 Jan 2011
TL;DR: In this paper, a maximum-a-posteriori reconstruction algorithm for jointly estimating the attenuation and activity distributions from TOF-PET data was proposed, which does not suffer from the cross-talk that was observed in non-TOF PET.
Abstract: In positron emission tomography (PET) and single photon emission tomography (SPECT), attenuation correction is necessary for quantitative reconstruction of the tracer distribution. Previously, several attempts have been undertaken to estimate the attenuation coefficients from emission data only. These attempts had limited success, because the problem does not have a unique solution, and severe and persistent “cross-talk” between the estimated activity and attenuation distributions was observed. In this paper, we show that the availability of TOF-information eliminates the cross-talk problem by destroying symmetries in the associated Fisher information matrix. We propose a maximum-a-posteriori reconstruction algorithm for jointly estimating the attenuation and activity distributions from TOF-PET data. The performance of the algorithm is studied with 2D simulations, and further illustrated with phantom experiments and with a patient scan. The estimated attenuation image is robust to noise, and does not suffer from the cross-talk that was observed in non-TOF PET. However, some constraining is still mandatory, because the TOF-data determine the attenuation sinogram only up to a constant offset.

Journal ArticleDOI
TL;DR: A full-wave equation that describes nonlinear propagation in a heterogeneous attenuating medium is solved numerically with finite differences in the time domain (FDTD) to simulate propagation of a diagnostic ultrasound pulse through a measured representation of the human abdomen with heterogeneities in speed of sound, attenuation, density, and nonlinearity.
Abstract: A full-wave equation that describes nonlinear propagation in a heterogeneous attenuating medium is solved numerically with finite differences in the time domain (FDTD). This numerical method is used to simulate propagation of a diagnostic ultrasound pulse through a measured representation of the human abdomen with heterogeneities in speed of sound, attenuation, density, and nonlinearity. Conventional delay-andsum beamforming is used to generate point spread functions (PSF) that display the effects of these heterogeneities. For the particular imaging configuration that is modeled, these PSFs reveal that the primary source of degradation in fundamental imaging is reverberation from near-field structures. Reverberation clutter in the harmonic PSF is 26 dB higher than the fundamental PSF. An artificial medium with uniform velocity but unchanged impedance characteristics indicates that for the fundamental PSF, the primary source of degradation is phase aberration. An ultrasound image is created in silico using the same physical and algorithmic process used in an ultrasound scanner: a series of pulses are transmitted through heterogeneous scattering tissue and the received echoes are used in a delay-and-sum beamforming algorithm to generate images. These beamformed images are compared with images obtained from convolution of the PSF with a scatterer field to demonstrate that a very large portion of the PSF must be used to accurately represent the clutter observed in conventional imaging.

Journal ArticleDOI
TL;DR: In this article, a bandpass aperture type cross-dipole frequency selective surface (FSS) is designed and etched in the coatings of energy-saving glass to improve the transmission of useful signals while preserving IR attenuation as much as possible.
Abstract: Energy-saving glass is becoming very popular in building design due to their effective shielding of building interior against heat entering the building with infrared (IR) waves. This is obtained by depositing a thin layer of metallic-oxide on the glass surface using special sputtering processes. This layer attenuates IR waves and hence keeps buildings cooler in summer and warmer in winter. However, this resistive coating also attenuates useful microwave/RF signals required for mobile phone, GPS and personal communication systems etc. by as much as 30 dB. To overcome this drawback, a bandpass aperture type cross-dipole frequency selective surface (FSS) is designed and etched in the coatings of energy-saving glass to improve the transmission of useful signals while preserving IR attenuation as much as possible. With this FSS, 15-18 dB peak transmission improvement can be achieved, for waves incident with ±45° from normal for both TE and TM polarizations. Theoretical and measured results are presented.

Journal ArticleDOI
TL;DR: It is found that for a given X-ray spectrum and object size the SNR in the heavy element's basis material image peaks for a certain atomic number Z, which is the natural, exponential decrease with particularly deteriorating effects in the case where the attenuation from the object itself causes a total signal loss below the K-edge.
Abstract: The feasibility of K-edge imaging using energy-resolved, photon-counting transmission measurements in X-ray computed tomography (CT) has been demonstrated by simulations and experiments. The method is based on probing the discontinuities of the attenuation coefficient of heavy elements above and below the K-edge energy by using energy-sensitive, photon counting X-ray detectors. In this paper, we investigate the dependence of the sensitivity of K-edge imaging on the atomic number Z of the contrast material, on the object diameter D , on the spectral response of the X-ray detector and on the X-ray tube voltage. We assume a photon-counting detector equipped with six adjustable energy thresholds. Physical effects leading to a degradation of the energy resolution of the detector are taken into account using the concept of a spectral response function R(E,U) for which we assume four different models. As a validation of our analytical considerations and in order to investigate the influence of elliptically shaped phantoms, we provide CT simulations of an anthropomorphic Forbild-Abdomen phantom containing a gold-contrast agent. The dependence on the values of the energy thresholds is taken into account by optimizing the achievable signal-to-noise ratios (SNR) with respect to the threshold values. We find that for a given X-ray spectrum and object size the SNR in the heavy element's basis material image peaks for a certain atomic number Z. The dependence of the SNR in the high-Z basis-material image on the object diameter is the natural, exponential decrease with particularly deteriorating effects in the case where the attenuation from the object itself causes a total signal loss below the K-edge. The influence of the energy-response of the detector is very important. We observed that the optimal SNR values obtained with an ideal detector and with a CdTe pixel detector whose response, showing significant tailing, has been determined at a synchrotron differ by factors of about two to three. The potentially very important impact of scattered X-ray radiation and pulse pile-up occurring at high photon rates on the sensitivity of the technique is qualitatively discussed.

Journal ArticleDOI
TL;DR: In this article, the authors assess the potential of the visco-acoustic frequency domain full-waveform inversion (FWI) to reconstruct P-wave velocity (VP) and p-wave attenuation factor (Q) from surface onshore seismic data.
Abstract: SUMMARY Here we assess the potential of the visco-acoustic frequency domain full-waveform inversion (FWI) to reconstruct P-wave velocity (VP) and P-wave attenuation factor (Q) from surface onshore seismic data. First, we perform a sensitivity analysis of the FWI based upon a grid search analysis of the misfit function and several synthetic FWI examples using velocity and Q models of increasing complexity. Subsequently, we applied both the acoustic and viscoacousticFWItorealsurfacewide-apertureonshoreseismicdatafromthePolishBasin,wherea strongattenuationoftheseismicdataisobserved.Thesensitivityanalysisofthevisco-acoustic FWI suggests that the FWI can jointly reconstruct the velocity and the attenuation factor if the signatureoftheattenuationissufficientlystronginthedata.Asyntheticexamplecorresponding to a homogeneous background model with an inclusion shows a reliable reconstruction of VP and Q in the inclusion, when Q is as small as 90 and 50 in the background model and in the inclusion, respectively. A first application of acoustic FWI to real data shows that a heuristic normalization of the data with offset allows us to compensate for the effect of the attenuation in the data and reconstruct a reliable velocity model. Alternatively, we show that visco-acoustic FWI allows us to reconstruct jointly both a reliable velocity model and a Q model from the true-amplitude data. We propose a pragmatical approach based upon seismic modelling and source wavelet estimation to infer the best starting homogeneous Q model for visco-acoustic FWI. We find the source wavelet estimation quite sensitive to the quality of the velocity and attenuation models used for the estimation. For example, source-to-source wavelets are significantly more consistent when computed in the final FWI model than in the initial one. A good kinematic and amplitude match between the early-arriving phases of the real and time-domain synthetic seismograms computed in the final FWI model provides an additional evidence of the reliability of the final FWI model. We find the recovered velocity and attenuation models consistent with the expected lithology and stratigraphy in the study area. We link high-attenuation zones with the increased clay content and the presence of the mineralized fluids.

Journal ArticleDOI
TL;DR: When using a segmented attenuation map, at least five different tissue types should be considered: cortical bone, spongeous bone, soft tissue, lung, and air, and the interpatient variability of lung attenuation coefficients should be taken into account.
Abstract: Purpose: Accurate attenuation correction is important for PET quantification. Often, a segmented attenuation map is used, especially in MRI-based attenuation correction. As deriving the attenuation map from MRIimages is difficult, different errors can be present in the segmented attenuation map. The goal of this paper is to determine the effect of these errors on quantification. Methods: The authors simulated the digital XCAT phantom using the GATE Monte Carlo simulation framework and a model of the Philips Gemini TF. A whole body scan was simulated, spanning an axial field of view of 70 cm. A total of fifteen lesions were placed in the lung,liver, spine, colon, prostate, and femur. The acquired data were reconstructed with a reference attenuation map and with different attenuation maps that were modified to reflect common segmentation errors. The quantitative difference between reconstructed images was evaluated. Results: Segmentation into five tissue classes, namely cortical bone, spongeous bone, soft tissue,lung, and air yielded errors below 5%. Large errors were caused by ignoring lungtissue (up to 45%) or cortical bone (up to 17%). The interpatient variability of lung attenuation coefficients can lead to errors of 10% and more. Up to 20% tissue misclassification from bone to soft tissue yielded errors below 5%. The same applies for up to 10% misclassification from lung to air. Conclusions: When using a segmented attenuation map, at least five different tissue types should be considered: cortical bone, spongeous bone, soft tissue,lung, and air. Furthermore, the interpatient variability of lung attenuation coefficients should be taken into account. Limited misclassification from bone to soft tissue and from lung to air is acceptable, as these do not lead to relevant errors.

Journal ArticleDOI
TL;DR: In this article, the authors measured the steady state viscosity of polycrystalline materials at high homologous temperatures (T = 2 2 −48°C; T/Tm = 0.61 −0.67) and various grain sizes (3 −22 mm), the growth of which was controlled by annealing.
Abstract: attenuation QE1 were measured accurately over a broad frequency range (10 −4 ≤ f (Hz) ≤ 2.15) and at low strain amplitude (10 −5 –10 −6 ). Creep experiments were performed with the same apparatus to measure the steady state viscosity. Anelasticity and viscosity were measured at high homologous temperatures (T =2 2 –48°C; T/Tm = 0.61–0.67) and various grain sizes (3–22 mm), the growth of which was controlled by annealing. Using the measured viscosities h and the unrelaxed modulus EU determined from ultrasonic experiments, the frequency of the entire data set was normalized by the Maxwell frequency fM = EU/h, resulting in E and Q −1 master curves. The Q −1 data from previous studies on olivine‐dominated samples also collapse onto the same curve when scaled by fM,, demonstrating the universality of anelasticity for polycrystalline materials. The similitude by the Maxwell frequency scaling indicates that the dominant mechanism for the anelasticity observed in this study and in previous studies is diffusionally accommodated grain boundary sliding. A generalized formulation for this similitude is provided to extrapolate the experimental data to velocity and attenuation of seismic shear waves.

Journal ArticleDOI
TL;DR: In this article, the transmission loss and propagation loss of RF waves penetrating into fresh water are calculated for various propagation depths for air-to-water communications from air to fresh water.
Abstract: Wireless communications from air to fresh water are studied here. Our analysis relies on plane wave propagation mod-els. Specifically, the transmission loss and propagation loss of RF waves penetrating into fresh water are calculated for various propagation depths. Even though RF wireless communications are not well suited for seawater due to its high attenuation, our paper illustrates that RF communications from air to fresh water are possible. Finally, this work de-rives the optimum frequencies, which provide minimum attenuation and maximum propagation depth, for RF commu-nications from air to fresh water.

Journal ArticleDOI
TL;DR: In this article, the Fourier spectral amplitude of small-to-moderate earthquakes, recorded on the local seismic networks, is modeled by Q and κ, while apparent geometrical spreading models account for the frequency-independent decay of energy.
Abstract: SUMMARY Modelling the attenuation of shear wave energy is an important component of seismic hazard analysis. Previous studies have shown how attenuation, particularly in the uppermost layers of the crust, is regionally dependent. The impact of this is that the decay of energy radiating from an earthquake will vary from place to place. To quantify the regional attenuation in Switzerland we model the Fourier spectral amplitude of small-to-moderate earthquakes, recorded on the local seismic networks. High-frequency decay is parametrized by Q and κ, while apparent geometrical spreading models account for the frequency-independent decay of energy. We analyse ground motion encompassing the significant duration of shaking to provide models that are useful for the purpose of seismic hazard analysis. Two methods are used to estimate the whole path attenuation parameter, t*: first, a simultaneous fit of the source model and attenuation effects across the entire spectral bandwidth for earthquakes with M > 2; and secondly, a linear fit of an attenuation model to the high-frequency part of the spectrum for earthquakes with M > 3.5. The t* parameter is found to vary with hypocentral distance consistent with a weakly attenuating crust and strongly attenuating uppermost layer. 1-D tomographic inversions indicate a profile of increasing Q with depth down to the Moho. Frequency-independent decay is parametrized using a three-part model which allows for the inclusion of Moho reflection phases in the spectrum in the range of 20–140 km in the Swiss Foreland and from 70 to 140 km in the Swiss Alps.

Journal ArticleDOI
TL;DR: Computer simulations were used to compare the accuracy and the precision of the three methods that are commonly used to estimate the local ultrasonic attenuation within a region of interest (ROI) in tissue, and showed that when the ROI was homogeneous, the spectral difference method had the highest accuracy and precision.
Abstract: Accurate characterization of tissue pathologies using ultrasonic attenuation is strongly dependent on the accuracy of the algorithm that is used to obtain the attenuation coefficient estimates. In this paper, computer simulations were used to compare the accuracy and the precision of the three methods that are commonly used to estimate the local ultrasonic attenuation within a region of interest (ROI) in tissue; namely, the spectral log difference method, the spectral difference method, and the hybrid method. The effects of the inhomgeneities within the ROI on the accuracy of the three algorithms were studied, and the optimal ROI size (the number of independent echoes laterally and the number of pulse lengths axially) was quantified for each method. The three algorithms were tested for when the ROI was homogeneous, the ROI had variations in scatterer number density, and the ROI had variations in effective scatterer size. The results showed that when the ROI was homogeneous, the spectral difference method had the highest accuracy and precision followed by the spectral log difference method and the hybrid method, respectively. Also, when the scatterer number density varied, the spectral difference method completely failed, while the log difference method and hybrid method still gave good results. Lastly, when the scatterer size varied, all of the methods failed.

Journal ArticleDOI
TL;DR: A 'hybrid' imaging approach which can image both light absorption properties and acoustic transmission properties of an object in a two-dimensional slice using a computed tomography (CT) photoacoustic imager is presented.
Abstract: We present a ‘hybrid’ imaging approach which can image both light absorption properties and acoustic transmission properties of an object in a two-dimensional slice using a computed tomography (CT) photoacoustic imager. The ultrasound transmission measurement method uses a strong optical absorber of small cross-section placed in the path of the light illuminating the sample. This absorber, which we call a passive element acts as a source of ultrasound. The interaction of ultrasound with the sample can be measured in transmission, using the same ultrasound detector used for photoacoustics. Such measurements are made at various angles around the sample in a CT approach. Images of the ultrasound propagation parameters, attenuation and speed of sound, can be reconstructed by inversion of a measurement model. We validate the method on specially designed phantoms and biological specimens. The obtained images are quantitative in terms of the shape, size, location, and acoustic properties of the examined heterogeneities