Deconvolution-based deblurring of reconstructed images in photoacoustic/thermoacoustic tomography

Photoacoustic/thermoacoustic tomography is an emerging hybrid imaging modality combining optical/microwave imaging with ultrasound imaging. Here, a k-wave MATLAB toolbox was used to simulate various configurations of excitation pulse shape, width, transducer types, and target object sizes to see their effect on the photoacoustic/thermoacoustic signals. A numerical blood vessel phantom was also used to demonstrate the effect of various excitation pulse waveforms and pulse widths on the reconstructed images. Reconstructed images were blurred due to the broadening of the pressure waves by the excitation pulse width as well as by the limited transducer bandwidth. The blurring increases with increase in pulse width. A deconvolution approach is presented here with Tikhonov regularization to correct the photoacoustic/thermoacoustic signals, which resulted in improved reconstructed images by reducing the blurring effect. It is observed that the reconstructed images remain unaffected by change in pulse widths or pulse shapes, as well as by the limited bandwidth of the ultrasound detectors after the use of the deconvolution technique. (C) 2013 Optical Society of America

Monte Carlo simulation of light transport in tissue for optimizing light delivery in photoacoustic imaging of the sentinel lymph node

Noninvasive or minimally invasive identification of sentinel lymph node (SLN) is essential to reduce the surgical effects of SLN biopsy. Photoacoustic (PA) imaging of SLN in animal models has shown its promise for clinical use in the future. Here, we present a Monte Carlo simulation for light transport in the SLN for various light delivery configurations with a clinical ultrasound probe. Our simulation assumes a realistic tissue layer model and also can handle the transmission/reflectance at SLN-tissue boundary due to the mismatch of refractive index. Various light incidence angles show that for deeply situated SLNs the maximum absorption of light in the SLN is for normal incidence. We also show that if a part of the diffused reflected photons is reflected back into the skin using a reflector, the absorption of light in the SLN can be increased significantly to enhance the PA signal. (C) 2013 Society of Photo-Optical Instrumentation Engineers (SPIE)

Improving tangential resolution with a modified delay-and-sum reconstruction algorithm in photoacoustic and thermoacoustic tomography

Spatial resolution in photoacoustic and thermoacoustic tomography is ultrasound transducer (detector) bandwidth limited. For a circular scanning geometry the axial (radial) resolution is not affected by the detector aperture, but the tangential (lateral) resolution is highly dependent on the aperture size, and it is also spatially varying (depending on the location relative to the scanning center). Several approaches have been reported to counter this problem by physically attaching a negative acoustic lens in front of the nonfocused transducer or by using virtual point detectors. Here, we have implemented a modified delay-and-sum reconstruction method, which takes into account the large aperture of the detector, leading to more than fivefold improvement in the tangential resolution in photoacoustic (and thermoacoustic) tomography. Three different types of numerical phantoms were used to validate our reconstruction method. It is also shown that we were able to preserve the shape of the reconstructed objects with the modified algorithm. (C) 2014 Optical Society of America

Basis pursuit deconvolution for improving model-based reconstructed images in photoacoustic tomography

The model-based image reconstruction approaches in photoacoustic tomography have a distinct advantage compared to traditional analytical methods for cases where limited data is available. These methods typically deploy Tikhonov based regularization scheme to reconstruct the initial pressure from the boundary acoustic data. The model-resolution for these cases represents the blur induced by the regularization scheme. A method that utilizes this blurring model and performs the basis pursuit deconvolution to improve the quantitative accuracy of the reconstructed photoacoustic image is proposed and shown to be superior compared to other traditional methods via three numerical experiments. Moreover, this deconvolution including the building of an approximate blur matrix is achieved via the Lanczos bidagonalization (least-squares QR) making this approach attractive in real-time. (C) 2014 Optical Society of America

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.

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.

A Method for Delineation of Bone Surfaces in Photoacoustic Computed Tomography of the Finger

Photoacoustic (PA) imaging of interphalangeal peripheral joints is of interest in the context of using the synovial membrane as a surrogate marker of rheumatoid arthritis. Previous work has shown that ultrasound (US) produced by absorption of light at the epidermis reflects on the bone surfaces within the finger. When the reflected signals are backprojected in the region of interest, artifacts are produced, confounding interpretation of the images. In this work, we present an approach where the PA signals known to originate from the epidermis are treated as virtual US transmitters, and a separate reconstruction is performed as in US reflection imaging. This allows us to identify the bone surfaces. Furthermore, the identification of the joint space is important as this provides a landmark to localize a region-of-interest in seeking the inflamed synovial membrane. The ability to delineate bone surfaces allows us to identify not only the artifacts but also the interphalangeal joint space without recourse to new US hardware or a new measurement. We test the approach on phantoms and on a healthy human finger.

Studies of ambient and chamber aerosol composition using the aerosol mass spectrometer

This thesis presents composition measurements for atmospherically relevant inorganic and organic aerosol from laboratory and ambient measurements using the Aerodyne aerosol mass spectrometer. Studies include the oxidation of dodecane in the Caltech environmental chambers, and several aircraft- and ground-based field studies, which include the quantification of wildfire emissions off the coast of California, and Los Angeles urban emissions.

The oxidation of dodecane by OH under low NO conditions and the formation of secondary organic aerosol (SOA) was explored using a gas-phase chemical model, gas-phase CIMS measurements, and high molecular weight ion traces from particle- phase HR-TOF-AMS mass spectra. The combination of these measurements support the hypothesis that particle-phase chemistry leading to peroxyhemiacetal formation is important. Positive matrix factorization (PMF) was applied to the AMS mass spectra which revealed three factors representing a combination of gas-particle partitioning, chemical conversion in the aerosol, and wall deposition.

Airborne measurements of biomass burning emissions from a chaparral fire on the central Californian coast were carried out in November 2009. Physical and chemical changes were reported for smoke ages 0 – 4 h old. CO₂ normalized ammonium, nitrate, and sulfate increased, whereas the normalized OA decreased sharply in the first 1.5 - 2 h, and then slowly increased for the remaining 2 h (net decrease in normalized OA). Comparison to wildfire samples from the Yucatan revealed that factors such as relative humidity, incident UV radiation, age of smoke, and concentration of emissions are important for wildfire evolution.

Ground-based aerosol composition is reported for Pasadena, CA during the summer of 2009. The OA component, which dominated the submicron aerosol mass, was deconvolved into hydrocarbon-like organic aerosol (HOA), semi-volatile oxidized organic aerosol (SVOOA), and low-volatility oxidized organic aerosol (LVOOA). The HOA/OA was only 0.08–0.23, indicating that most of Pasadena OA in the summer months is dominated by oxidized OA resulting from transported emissions that have undergone photochemistry and/or moisture-influenced processing, as apposed to only primary organic aerosol emissions. Airborne measurements and model predictions of aerosol composition are reported for the 2010 CalNex field campaign.