Modelling of errors due to speed of sound variations in photoacoustic tomography using a Bayesian framework

Inverse problem of estimating initial pressure in photoacoustic tomography is ill-posed and thus sensitive to errors in modelling and measurements. In practical experiments, accurate knowledge of the speed of sound of the imaged target is commonly not available, and therefore an approximate speed of sound is used in the computational model. This can result in errors in the solution of the inverse problem that can appear as artefacts in the reconstructed images. In this paper, the inverse problem of photoacoustic tomography is approached in a Bayesian framework. Errors due to uncertainties in the speed of sound are modelled using Bayesian approximation error modelling. Estimation of the initial pressure distribution together with information on the reliability of these estimates are considered. The approach was studied using numerical simulations. The results show that uncertainties in the speed of sound can cause significant errors in the solution of the inverse problem. However, modelling of these uncertainties improves the accuracy of the solution

Bayesian approach to image reconstruction in photoacoustic tomography

Photoacoustic tomography is a hybrid imaging method that has a variety of biomedical applications. In photoacoustic tomography, the image reconstruction problem (inverse problem) is to resolve the initial pressure distribution from detected ultrasound waves generated within an object due to an illumination of a short light pulse. In this work, this problem is approached in Bayesian framework. Image reconstruction is investigated with numerical simulations in different detector geometries, including limited view setup, and utilizing different prior information. Furthermore, assessing the reliability of the estimates is investigated. The simulations show that the Bayesian approach can produce accurate estimates of the initial pressure distribution and uncertainty information even in a limited view setup if proper prior information is utilized

Modelling of errors and uncertainties in photoacoustic tomography using a Bayesian framework

Photoacoustic tomography is studied in the framework of Bayesian inverse problems. Modelling of errors and uncertainties using Bayesian approximation error modelling is investigated. The approach is tested with simulation

Computationally Efficient Forward Operator for Photoacoustic Tomography Based on Coordinate Transformations

IEEE Photoacoustic tomography (PAT) is an imaging modality that utilizes the photoacoustic effect. In PAT, a photoacoustic image is computed from measured data by modeling ultrasound propagation in the imaged domain and solving an inverse problem utilizing a discrete forward operator. However, in realistic measurement geometries with several ultrasound transducers and relatively large imaging volume, an explicit formation and use of the forward operator can be computationally prohibitively expensive. In this work, we propose a transformation based approach for efficient modeling of photoacoustic signals and reconstruction of photoacoustic images. In the approach, the forward operator is constructed for a reference ultrasound transducer and expanded into a general measurement geometry using transformations that map the formulated forward operator in local coordinates to the global coordinates of the measurement geometry. The inverse problem is solved using a Bayesian framework. The approach is evaluated with numerical simulations and experimental data. The results show that the proposed approach produces accurate three-dimensional photoacoustic images with a significantly reduced computational cost both in memory requirements and in time. In the studied cases, depending on the computational factors such as discretization, over 30-fold reduction in memory consumption and was achieved without a reduction in image quality compared to a conventional approach

Fluorescence measurements show stronger cold inhibition of photosynthetic light reactions in Scots pine compared to Norway spruce as well as during spring compared to autumn

We studied the photosynthetic activity of Scots pine (Pinus sylvestris L.) and Norway spruce (Picea abies [L.] Karst) in relation to air temperature changes from March 2013 to February 2014. We measured the chlorophyll fluorescence of approximately 50 trees of each species growing in southern Finland. Fluorescence was measured 13 times per week. We began by measuring shoots present in late winter (i.e., March 2013) before including new shoots once they started to elongate in spring. By July, when the spring shoots had achieved similar fluorescence levels to the older ones, we proceeded to measure the new shoots only.We analysed the data by fitting a sigmoidal model containing four parameters to link sliding averages of temperature and fluorescence. A parameter defining the temperature range over which predicted fluorescence increased most rapidly was the most informative with in describing temperature dependence of fluorescence.The model generated similar fluorescence patterns for both species, but differences were observed for critical temperature and needle age. Down regulation of the light reaction was stronger in spring than in autumn. Pine showed more conservative control of the photosynthetic light reactions, which were activated later in spring and more readily attenuated in autumn. Under the assumption of a close correlation of fluorescence and photosynthesis, spruce should therefore benefit more than pine from the increased photosynthetic potential during warmer springs, but be more likely to suffer frost damage with a sudden cooling following a warm period. The winter of 20132014 was unusually mild and similar to future conditions predicted by global warming models. During the mild winter, the activity of photosynthetic light reactions of both conifers, especially spruce, remained high. Because light levels during winter are too low for photosynthesis, this activity may translate to a net carbon loss due to respiration

Substorm occurrence during quiet solar wind driving

Peer reviewe

Characterization of ultrasound fields using a potential optical flow based synthetic schlieren tomography

Synthetic schlieren tomography is an optical imaging method for characterization of ultrasound fields based on observing bending of light due to acousto-optic effect. In this work, potential optical flow based pressure estimation method is introduced

A study of omega bands and Ps6 pulsations on the ground, at low altitude and at geostationary orbit

We investigate the electrodynamic coupling between auroral omega bands and the inner magnetosphere. The goal of this study is to determine the features to which omega bands map in the magnetosphere. To establish the auroral-magnetosphere connection, we appeal to the case study analysis of the data rich event of September 26, 1989. At 6 magnetic local time (MLT), two trains of Ps6 pulsations (ground magnetic signatures of omega bands) were observed to drift over the Canadian Auroral Network For the OPEN Program Unified Study (CANOPUS) chain. At the same time periodic ionospheric flow patterns moved through the collocated Bistatic Auroral Radar System (BARS) field of view. Similar coincident magnetic variations were observed by GOES 6, GOES 7 and SCATHA, all of which had magnetic foot points near the CANOPUS/BARS stations. SCATHA, which was located at 6 MLT, 0.5 RE earthward of GOES 7 observed the 10 min period pulsations, whereas GOES 7 did not. In addition, DMSP F6 and F8 were over-flying the region and observed characteristic precipitation and flow signatures. From this fortunate constellation of ground and space observations, we conclude that auroral omega bands are the electrodynamic signature of a corrugated current sheet (or some similar spatially localized magnetic structure) in the near-Earth geostationary magnetosphere

On the characterization of magnetic reconnection in global MHD simulations

The conventional definition of reconnection rate as the electric field parallel to an x-line is problematic in global MHD simulations for several reasons: the x-line itself may be hard to find in a non-trivial geometry such as at the magnetopause, and the lack of realistic resistivity modelling leaves us without reliable non-convective electric field. In this article we describe reconnection characterization methods that avoid those problems and are practical to apply in global MHD simulations. We propose that the reconnection separator line can be identified as the region where magnetic field lines of different topological properties meet, rather than by local considerations. The global convection associated with reconnection is then quantified by calculating the transfer of mass, energy or magnetic field across the boundary of closed and open field line regions. The extent of the diffusion region is determined from the destruction of electromagnetic energy, given by the divergence of the Poynting vector. Integrals of this energy conversion provide a way to estimate the total reconnection efficiency

The magnetotail reconnection region in a global MHD simulation

International audienceThis work investigates the nature and the role of magnetic reconnection in a global magnetohydrodynamic simulation of the magnetosphere. We use the Gumics-4 simulation to study reconnection that occurs in the near-Earth region of the current sheet in the magnetotail. We locate the current sheet surface and the magnetic x-line that appears when reconnection starts. We illustrate the difference between quiet and active states of the reconnection region: variations in such quantities as the current sheet thickness, plasma flow velocities, and Poynting vector divergence are strong. A characteristic feature is strong asymmetry caused by non-perpendicular inflows. We determine the reconnection efficiency by the net rate of Poynting flux into the reconnection region. The reconnection efficiency in the simulation is directly proportional to the energy flux into the magnetosphere through the magnetopause: about half of all energy flowing through the magnetosphere is converted from an electromagnetic into a mechanical form in the reconnection region. Thus, the tail reconnection that is central to the magnetospheric circulation is directly driven; the tail does not exhibit a cycle of storage and rapid release of magnetic energy. We find similar behaviour of the tail in both synthetic and real event runs