Born-ratio type data normalization improves quantitation in photoacoustic tomography

In this report, we present a Born-ratio type of data normalization for reconstruction of initial acoustic pressure distribution in photoacoustic tomography (PAT). The normalized Born-ratio type of data is obtained as a ratio of photoacoustic pressure obtained with tissue sample in a coupling medium to the one obtained using purely coupling medium. It is shown that this type of data normalization improves the quantitation (intrinsic contrast) of the reconstructed images in comparison to the traditional techniques (unnormalized) that are currently available in PAT. Studies are carried out using various tissue samples. The robustness of the proposed method is studied at various noise levels added to the collected data. The improvement in quantitation can enable accurate estimation of pathophysiological parameter (optical absorption coefficient, a) of tissue sample under investigation leading to better sensitivity in PAT.

Intense Acoustic Burst Ultrasound Modulated Optical Tomography for Elasticity Mapping of Soft Biological Tissue Mimicking Phantom: A Laser Speckle Contrast Analysis Study

This report addresses the assessment of variation in elastic property of soft biological tissues non-invasively using laser speckle contrast measurement. The experimental as well as the numerical (Monte-Carlo simulation) studies are carried out. In this an intense acoustic burst of ultrasound (an acoustic pulse with high power within standard safety limits), instead of continuous wave, is employed to induce large modulation of the tissue materials in the ultrasound insonified region of interest (ROI) and it results to enhance the strength of the ultrasound modulated optical signal in ultrasound modulated optical tomography (UMOT) system. The intensity fluctuation of speckle patterns formed by interference of light scattered (while traversing through tissue medium) is characterized by the motion of scattering sites. The displacement of scattering particles is inversely related to the elastic property of the tissue. We study the feasibility of laser speckle contrast analysis (LSCA) technique to reconstruct a map of the elastic property of a soft tissue-mimicking phantom. We employ source synchronized parallel speckle detection scheme to (experimentally) measure the speckle contrast from the light traversing through ultrasound (US) insonified tissue-mimicking phantom. The measured relative image contrast (the ratio of the difference of the maximum and the minimum values to the maximum value) for intense acoustic burst is 86.44 % in comparison to 67.28 % for continuous wave excitation of ultrasound. We also present 1-D and 2-D image of speckle contrast which is the representative of elastic property distribution.

Fienup Algorithm With Sparsity Constraints: Application to Frequency-Domain Optical-Coherence Tomography

We address the problem of reconstructing a sparse signal from its DFT magnitude. We refer to this problem as the sparse phase retrieval (SPR) problem, which finds applications in tomography, digital holography, electron microscopy, etc. We develop a Fienup-type iterative algorithm, referred to as the Max-K algorithm, to enforce sparsity and successively refine the estimate of phase. We show that the Max-K algorithm possesses Cauchy convergence properties under certain conditions, that is, the MSE of reconstruction does not increase with iterations. We also formulate the problem of SPR as a feasibility problem, where the goal is to find a signal that is sparse in a known basis and whose Fourier transform magnitude is consistent with the measurement. Subsequently, we interpret the Max-K algorithm as alternating projections onto the object-domain and measurement-domain constraint sets and generalize it to a parameterized relaxation, known as the relaxed averaged alternating reflections (RAAR) algorithm. On the application front, we work with measurements acquired using a frequency-domain optical-coherence tomography (FDOCT) experimental setup. Experimental results on measured data show that the proposed algorithms exhibit good reconstruction performance compared with the direct inversion technique, homomorphic technique, and the classical Fienup algorithm without sparsity constraint; specifically, the autocorrelation artifacts and background noise are suppressed to a significant extent. We also demonstrate that the RAAR algorithm offers a broader framework for FDOCT reconstruction, of which the direct inversion technique and the proposed Max-K algorithm become special instances corresponding to specific values of the relaxation parameter.

Isotropic finite volume discretization

Finite volume methods traditionally employ dimension by dimension extension of the one-dimensional reconstruction and averaging procedures to achieve spatial discretization of the governing partial differential equations on a structured Cartesian mesh in multiple dimensions. This simple approach based on tensor product stencils introduces an undesirable grid orientation dependence in the computed solution. The resulting anisotropic errors lead to a disparity in the calculations that is most prominent between directions parallel and diagonal to the grid lines. In this work we develop isotropic finite volume discretization schemes which minimize such grid orientation effects in multidimensional calculations by eliminating the directional bias in the lowest order term in the truncation error. Explicit isotropic expressions that relate the cell face averaged line and surface integrals of a function and its derivatives to the given cell area and volume averages are derived in two and three dimensions, respectively. It is found that a family of isotropic approximations with a free parameter can be derived by combining isotropic schemes based on next-nearest and next-next-nearest neighbors in three dimensions. Use of these isotropic expressions alone in a standard finite volume framework, however, is found to be insufficient in enforcing rotational invariance when the flux vector is nonlinear and/or spatially non-uniform. The rotationally invariant terms which lead to a loss of isotropy in such cases are explicitly identified and recast in a differential form. Various forms of flux correction terms which allow for a full recovery of rotational invariance in the lowest order truncation error terms, while preserving the formal order of accuracy and discrete conservation of the original finite volume method, are developed. Numerical tests in two and three dimensions attest the superior directional attributes of the proposed isotropic finite volume method. Prominent anisotropic errors, such as spurious asymmetric distortions on a circular reaction-diffusion wave that feature in the conventional finite volume implementation are effectively suppressed through isotropic finite volume discretization. Furthermore, for a given spatial resolution, a striking improvement in the prediction of kinetic energy decay rate corresponding to a general two-dimensional incompressible flow field is observed with the use of an isotropic finite volume method instead of the conventional discretization. (C) 2014 Elsevier Inc. All rights reserved.

Super-Resolution Reconstruction in Frequency-Domain Optical-Coherence Tomography Using the Finite-Rate-of-Innovation Principle

The standard approach to signal reconstruction in frequency-domain optical-coherence tomography (FDOCT) is to apply the inverse Fourier transform to the measurements. This technique offers limited resolution (due to Heisenberg's uncertainty principle). We propose a new super-resolution reconstruction method based on a parametric representation. We consider multilayer specimens, wherein each layer has a constant refractive index and show that the backscattered signal from such a specimen fits accurately in to the framework of finite-rate-of-innovation (FRI) signal model and is represented by a finite number of free parameters. We deploy the high-resolution Prony method and show that high-quality, super-resolved reconstruction is possible with fewer measurements (about one-fourth of the number required for the standard Fourier technique). To further improve robustness to noise in practical scenarios, we take advantage of an iterated singular-value decomposition algorithm (Cadzow denoiser). We present results of Monte Carlo analyses, and assess statistical efficiency of the reconstruction techniques by comparing their performance against the Cramer-Rao bound. Reconstruction results on experimental data obtained from technical as well as biological specimens show a distinct improvement in resolution and signal-to-reconstruction noise offered by the proposed method in comparison with the standard approach.

Effect of hydrogen on the yielding behavior and shear transformation zone volume in metallic glass ribbons

The effect of hydrogen (H) charging on the shear yield strength (tau(max)) and shear transformation zone volume (Omega) of Ni-Nb-Zr metallic glass ribbons, with varying Zr content, were studied through the first pop-in loads during nanoindentation. Weight gain measurements after H charging and desorption studies were utilized to identify how the total H absorbed during charging is partitioned into mobile and immobile (or trapped) parts. These, in turn, indicate the significant role of H mobility in the amorphous structure on the yielding behavior. In high-Zr alloys, tau(max) increases significantly whereas Omega decreases. In low-Zr alloys, a slight decrease in tau(max) and increase in Omega were noted. These experimental observations are rationalized in terms of the mobility of the absorbed H in the amorphous structure and the possible role of it in the shear transformation zone dynamics during deformation of the metallic glass. (C) 2014 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Estimating the volume of glaciers in the Himalayan-Karakoram region using different methods

Ice volume estimates are crucial for assessing water reserves stored in glaciers. Due to its large glacier coverage, such estimates are of particular interest for the Himalayan-Karakoram (HK) region. In this study, different existing methodologies are used to estimate the ice reserves: three area-volume relations, one slope-dependent volume estimation method, and two ice-thickness distribution models are applied to a recent, detailed, and complete glacier inventory of the HK region, spanning over the period 2000-2010 and revealing an ice coverage of 40 775 km(2). An uncertainty and sensitivity assessment is performed to investigate the influence of the observed glacier area and important model parameters on the resulting total ice volume. Results of the two ice-thickness distribution models are validated with local ice-thickness measurements at six glaciers. The resulting ice volumes for the entire HK region range from 2955 to 4737 km(3), depending on the approach. This range is lower than most previous estimates. Results from the ice thickness distribution models and the slope-dependent thickness estimations agree well with measured local ice thicknesses. However, total volume estimates from area-related relations are larger than those from other approaches. The study provides evidence on the significant effect of the selected method on results and underlines the importance of a careful and critical evaluation.

A Gold Sensors Array for Imaging The Real Tissue Phantom in Electrical Impedance Tomography

Surface electrodes in Electrical Impedance Tomography (EIT) phantoms usually reduce the SNR of the boundary potential data due to their design and development errors. A novel gold sensors array with high geometric precision is developed for EIT phantoms to improve the resistivity image quality. Gold thin films are deposited on a flexible FR4 sheet using electro-deposition process to make a sixteen electrode array with electrodes of identical geometry. A real tissue gold electrode phantom is developed with chicken tissue paste and the fat cylinders as the inhomogeneity. Boundary data are collected using a USB based high speed data acquisition system in a LabVIEW platform for different inhomogeneity positions. Resistivity images are reconstructed using EIDORS and compared with identical stainless steel electrode systems. Image contrast parameters are calculated from the resistivity matrix and the reconstructed images are evaluated for both the phantoms. Image contrast and image resolution of resistivity images are improved with gold electrode array.

A new approach to depict bone surfaces in finger imaging using photoacoustic tomography

Imaging the vasculature close around the finger joints is of interest in the field of rheumatology. Locally increased vasculature in the synovial membrane of these joints can be a marker for rheumatoid arthritis. In previous work we showed that part of the photoacoustically induced ultrasound from the epidermis reflects on the bone surface within the finger. These reflected signals could be wrongly interpreted as new photoacoustic sources. In this work we show that a conventional ultrasound reconstruction algorithm, that considers the skin as a collection of ultrasound transmitters and the PA tomography probe as the detector array, can be used to delineate bone surfaces of a finger. This can in the future assist in the localization of the joint gaps. This can provide us with a landmark to localize the region of the inflamed synovial membrane. We test the approach on finger mimicking phantoms.

Photoacoustic tomography of the human finger: towards the assessment of inflammatory joint diseases

Inflammatory arthritis is often manifested in finger joints. The growth of new or withdrawal of old blood vessels can be a sensitive marker for these diseases. Photoacoustic (PA) imaging has great potential in this respect since it allows the sensitive and highly resolved visualization of blood. We systematically investigated PA imaging of finger vasculature in healthy volunteers using a newly developed PA tomographic system. We present the PA results which show excellent detail of the vasculature. Vessels with diameters ranging between 100 mu m and 1.5 mm are visible along with details of the skin, including the epidermis and the subpapillary plexus. The focus of all the studies is at the proximal and distal interphalangeal joints, and in the context of ultimately visualizing the inflamed synovial membrane in patients. This work is important in laying the foundation for detailed research into PA imaging of the phalangeal vasculature in patients suffering from rheumatoid arthritis.

Volume conservation issues in incompressible smoothed particle hydrodynamics

A divergence-free velocity field is usually sought in numerical simulations of incompressible fluids. We show that the particle methods that compute a divergence-free velocity field to achieve incompressibility suffer from a volume conservation issue when a finite time-step position update scheme is used. Further, we propose a deformation gradient based approach to arrive at a velocity field that reduces the volume conservation issues in free surface flows and maintains density uniformity in internal flows while retaining the simplicity of first order time updates. (C) 2015 Elsevier Inc. All rights reserved.

Biocompatibility property of 100% strontium-substituted SiO2-Al2O3-P2O5-CaO-CaF2 glass ceramics over 26 weeks implantation in rabbit model: Histology and micro-Computed Tomography analysis

One of the desired properties for any new biomaterial composition is its long-term stability in a suitable animal model and such property cannot be appropriately assessed by performing short-term implantation studies. While hydroxyapatite (HA) or bioglass coated metallic biomaterials are being investigated for in vivo biocompatibility properties, such study is not extensively being pursued for bulk glass ceramics. In view of their inherent brittle nature, the implant stability as well as impact of long-term release of metallic ions on bone regeneration have been a major concern. In this perspective, the present article reports the results of the in vivo implantation experiments carried out using 100% strontium (Sr)-substituted glass ceramics with the nominal composition of 4.5 SiO2-3Al(2)O(3)-1.5P(2)O(5)-3SrO-2SrF(2) for 26 weeks in cylindrical bone defects in rabbit model. The combination of histological and micro-computed tomography analysis provided a qualitative and quantitative understanding of the bone regeneration around the glass ceramic implants in comparison to the highly bioactive HA bioglass implants (control). The sequential polychrome labeling of bone during in vivo osseointegration using three fluorochromes followed by fluorescence microscopy observation confirmed homogeneous bone formation around the test implants. The results of the present study unequivocally confirm the long-term implant stability as well as osteoconductive property of 100% Sr-substituted glass ceramics, which is comparable to that of a known bioactive implant, that is, HA-based bioglass. (c) 2014 Wiley Periodicals, Inc. J Biomed Mater Res Part B: Appl Biomater, 103B: 1168-1179, 2015.

Role of the Molar Volume on Estimated Diffusion Coefficients

The role of the molar volume on the estimated diffusion parameters has been speculated for decades. The Matano-Boltzmann method was the first to be developed for the estimation of the variation of the interdiffusion coefficients with composition. However, this could be used only when the molar volume varies ideally or remains constant. Although there are no such systems, this method is still being used to consider the ideal variation. More efficient methods were developed by Sauer-Freise, Den Broeder, and Wagner to tackle this problem. However, there is a lack of research indicating the most efficient method. We have shown that Wagner's method is the most suitable one when the molar volume deviates from the ideal value. Similarly, there are two methods for the estimation of the ratio of intrinsic diffusion coefficients at the Kirkendall marker plane proposed by Heumann and van Loo. The Heumann method, like the Matano-Boltzmann method, is suitable to use only when the molar volume varies more or less ideally or remains constant. In most of the real systems, where molar volume deviates from the ideality, it is safe to use the van Loo method. We have shown that the Heumann method introduces large errors even for a very small deviation of the molar volume from the ideal value. On the other hand, the van Loo method is relatively less sensitive to it. Overall, the estimation of the intrinsic diffusion coefficient is more sensitive than the interdiffusion coefficient.

Optimum Parameter Selection in Sparse Reconstruction of Frequency-Domain Optical-Coherence Tomography Signals

For a multilayered specimen, the back-scattered signal in frequency-domain optical-coherence tomography (FDOCT) is expressible as a sum of cosines, each corresponding to a change of refractive index in the specimen. Each of the cosines represent a peak in the reconstructed tomogram. We consider a truncated cosine series representation of the signal, with the constraint that the coefficients in the basis expansion be sparse. An l(2) (sum of squared errors) data error is considered with an l(1) (summation of absolute values) constraint on the coefficients. The optimization problem is solved using Weiszfeld's iteratively reweighted least squares (IRLS) algorithm. On real FDOCT data, improved results are obtained over the standard reconstruction technique with lower levels of background measurement noise and artifacts due to a strong l(1) penalty. The previous sparse tomogram reconstruction techniques in the literature proposed collecting sparse samples, necessitating a change in the data capturing process conventionally used in FDOCT. The IRLS-based method proposed in this paper does not suffer from this drawback.

Morphological and functional characterisation of the burrow systems of six earthworm species (Lumbricidae)

Earthworm burrow systems are generally described based on postulated behaviours associated with the three ecological types. In this study, we used X-ray tomography to obtain 3D information on the burrowing behaviour of six very common anecic (Aporrectodea nocturna and Lumbricus terrestris) and endogeic (Aporrectodea rosea, Allolobophora chlorotica, Aporrectodea caliginosa, Aporrectodea icterica) earthworm species, introduced into repacked soil cores for 6 weeks. A simple water infiltration test, the Beerkan method, was also used to assess some functional properties of these burrow systems. Endogeic worms make larger burrow systems, which are more highly branched, less continuous and of smaller diameter, than those of anecic worms. Among the anecic species, L. terrestris burrow systems are shorter (9.2 vs 21.2 m) with a higher number (14.5 vs 23.5) of less branched burrows (12.2 vs 20.2 branches m(-1)), which are also wider (7.78 vs 5.16 mm) than those of A. nocturna. In comparison, the burrow systems made by endogeic species appeared similar to each other. However, A. rosea burrows were short and narrow, whereas A. icterica had a longer burrow system (15.7 m), more intense bioturbation intensity (refilled macropores or soil lateral compaction around them) and thus a greater number of burrows. Regarding water infiltration, anecic burrow systems were far more efficient due to open burrows linking the top and bottom of the cores. For endogeic species, we observed a linear relationship between burrow length and the water infiltration rate (R (2) = 0.49, p < 0.01). Overall, the three main characteristics significantly influencing water infiltration were burrow length, burrow number and bioturbation volume. This last characteristic highlighted the effect of burrow refilling by casts.