Fast Transforms for Acoustic Imaging-Part I: Theory

The classical approach for acoustic imaging consists of beamforming, and produces the source distribution of interest convolved with the array point spread function. This convolution smears the image of interest, significantly reducing its effective resolution. Deconvolution methods have been proposed to enhance acoustic images and have produced significant improvements. Other proposals involve covariance fitting techniques, which avoid deconvolution altogether. However, in their traditional presentation, these enhanced reconstruction methods have very high computational costs, mostly because they have no means of efficiently transforming back and forth between a hypothetical image and the measured data. In this paper, we propose the Kronecker Array Transform ( KAT), a fast separable transform for array imaging applications. Under the assumption of a separable array, it enables the acceleration of imaging techniques by several orders of magnitude with respect to the fastest previously available methods, and enables the use of state-of-the-art regularized least-squares solvers. Using the KAT, one can reconstruct images with higher resolutions than was previously possible and use more accurate reconstruction techniques, opening new and exciting possibilities for acoustic imaging.

Fast Transforms for Acoustic Imaging-Part II: Applications

In Part I [""Fast Transforms for Acoustic Imaging-Part I: Theory,"" IEEE TRANSACTIONS ON IMAGE PROCESSING], we introduced the Kronecker array transform (KAT), a fast transform for imaging with separable arrays. Given a source distribution, the KAT produces the spectral matrix which would be measured by a separable sensor array. In Part II, we establish connections between the KAT, beamforming and 2-D convolutions, and show how these results can be used to accelerate classical and state of the art array imaging algorithms. We also propose using the KAT to accelerate general purpose regularized least-squares solvers. Using this approach, we avoid ill-conditioned deconvolution steps and obtain more accurate reconstructions than previously possible, while maintaining low computational costs. We also show how the KAT performs when imaging near-field source distributions, and illustrate the trade-off between accuracy and computational complexity. Finally, we show that separable designs can deliver accuracy competitive with multi-arm logarithmic spiral geometries, while having the computational advantages of the KAT.

Effect of irradiation distance on image contrast in epi-optoacoustic imaging of human volunteers

In combined clinical optoacoustic (OA) and ultrasound (US) imaging, epi-mode irradiation and detection integrated into one single probe offers flexible imaging of the human body. The imaging depth in epi-illumination is, however, strongly affected by clutter. As shown in previous phantom experiments, the location of irradiation plays an important role in clutter generation. We investigated the influence of the irradiation geometry on the local image contrast of clinical images, by varying the separation distance between the irradiated area and the acoustic imaging plane of a linear ultrasound transducer in an automated scanning setup. The results for different volunteers show that the image contrast can be enhanced on average by 25% and locally by more than a factor of two, when the irradiated area is slightly separated from the probe. Our findings have an important impact on the design of future optoacoustic probes for clinical application.

Polymer optical fibre sensors for endoscopic optoacoustic imaging

Opto-acoustic imaging (OAI) shows particular promise for in-vivo biomedical diagnostics. Its applications include cardiovascular, gastrointestinal and urogenital systems imaging. Opto-acoustic endoscopy (OAE) allows the imaging of body parts through cavities permitting entry. The critical parameter is the physical size of the device, allowing compatibility with current technology, while governing flexibility of the distal end of the endoscope based on the needs of the sensor. Polymer optical fibre (POF) presents a novel approach for endoscopic applications and has been positively discussed and compared in existing publications. A great advantage can be obtained for endoscopy due to a small size and array potential to provide discrete imaging speed improvements. Optical fibre exhibits numerous advantages over conventional piezo-electric transducers, such as immunity from electromagnetic interference and a higher resolution at small sizes. Furthermore, micro structured polymer optical fibres offer over 12 times the sensitivity of silica fibre. We present a polymer fibre Bragg grating ultrasound detector with a core diameter of 125 microns. We discuss the ultrasonic signals received and draw conclusions on the opportunities and challenges of applying this technology in biomedical applications.

Microstructured polymer optical fibre sensors for opto-acoustic endoscopy

Opto-acoustic imaging is a growing field of research in recent years, providing functional imaging of physiological biomarkers, such as the oxygenation of haemoglobin. Piezo electric transducers are the industry standard detector for ultrasonics, but their limited bandwidth, susceptibility to electromagnetic interference and their inversely proportional sensitivity to size all affect the detector performance. Sensors based on polymer optical fibres (POF) are immune to electromagnetic interference, have lower acoustic impedance and a reduced Young's Modulus compared to silica fibres. Furthermore, POF enables the possibility of a wideband sensor and a size appropriate to endoscopy. Micro-structured POF (mPOF) used in an interferometric detector has been shown to be an order of magnitude more sensitive than silica fibre at 1 MHz and 3 times more sensitive at 10 MHz. We present the first opto-acoustic measurements obtained using a 4.7mm PMMA mPOF Bragg grating with a fibre diameter of 130 μm and present the lateral directivity pattern of a PMMA mPOF FBG ultrasound sensor over a frequency range of 1-50 MHz. We discuss the impact of the pattern with respect to the targeted application and draw conclusions on how to mitigate the problems encountered.

Interferometric technique for measuring terahertz antenna phase patterns

A quasi-optical interferometric technique capable of measuring antenna phase patterns without the need for a heterodyne receiver is presented. It is particularly suited to the characterization of terahertz antennas feeding power detectors or mixers employing quasi-optical local oscillator injection. Examples of recorded antenna phase patterns at frequencies of 1.4 and 2.5 THz using homodyne detectors are presented. To our knowledge, these are the highest frequency antenna phase patterns ever recovered. Knowledge of both the amplitude and phase patterns in the far field enable a Gauss-Hermite or Gauss-Laguerre beam-mode analysis to be carried out for the antenna, of importance in performance optimization calculations, such as antenna gain and beam efficiency parameters at the design and prototype stage of antenna development. A full description of the beam would also be required if the antenna is to be used to feed a quasi-optical system in the near-field to far-field transition region. This situation could often arise when the device is fitted directly at the back of telescopes in flying observatories. A further benefit of the proposed technique is simplicity for characterizing systems in situ, an advantage of considerable importance as in many situations, the components may not be removable for further characterization once assembled. The proposed methodology is generic and should be useful across the wider sensing community, e.g., in single detector acoustic imaging or in adaptive imaging array applications. Furthermore, it is applicable across other frequencies of the EM spectrum, provided adequate spatial and temporal phase stability of the source can be maintained throughout the measurement process. Phase information retrieval is also of importance to emergent research areas, such as band-gap structure characterization, meta-materials research, electromagnetic cloaking, slow light, super-lens design as well as near-field and virtual imaging applications.

Electromagnetic, rock magnetic, and geochemical properties of surficial sediments in Eckernförde Bay

Submarine groundwater discharge in coastal settings can massively modify the hydraulic and geochemical conditions of the seafloor. Resulting local anomalies in the morphology and physical properties of surface sediments are usually explored with seismo-acoustic imaging techniques. Controlled source electromagnetic imaging offers an innovative dual approach to seep characterization by its ability to detect pore-water electrical conductivity, hence salinity, as well as sediment magnetic susceptibility, hence preservation or diagenetic alteration of iron oxides. The newly developed electromagnetic (EM) profiler Neridis II successfully realized this concept for a first time with a high-resolution survey of freshwater seeps in Eckernförde Bay (SW Baltic Sea). We demonstrate that EM profiling, complemented and validated by acoustic as well as sample-based rock magnetic and geochemical methods, can create a crisp and revealing fingerprint image of freshwater seepage and related reductive alteration of near-surface sediments. Our findings imply that (1) freshwater penetrates the pore space of Holocene mud sediments by both diffuse and focused advection, (2) pockmarks are marked by focused freshwater seepage, underlying sand highs, reduced mud thickness, higher porosity, fining of grain size, and anoxic conditions, (3) depletion of Fe oxides, especially magnetite, is more pervasive within pockmarks due to higher concentrations of organic and sulfidic reaction partners, and (4) freshwater advection reduces sediment magnetic susceptibility by a combination of pore-water injection (dilution) and magnetite reduction (depletion). The conductivity vs. susceptibility biplot resolves subtle lateral litho- and hydrofacies variations.

Investigação de aquífero fraturado para entendimento de fluxo e transporte de contaminantes clorados: estudo de caso em Valinhos, SP

Casos de contaminação de aquíferos fraturados são bastante complexos, tendo em vista a heterogeneidade das redes de fraturas, e no geral, sua investigação demanda a utilização de técnicas pouco usuais, como por exemplo o imageamento acústico e a perfilagem de velocidade de fluxo de água. Na área de estudo, localizada em Valinhos/SP, o uso inadequado de solventes organoclorados no passado ocasionou a contaminação do aquífero raso em duas áreas, e o aparecimento de concentrações no aquifero profundo levaram a condução do atual trabalho, que teve como principal objetivo a elaboração de um modelo conceitual de fluxo de água e transporte de contaminantes no aquífero cristalino. Previamente à investigação do aquífero fraturado, foi realizada uma análise de trabalhos existentes, incluindo a interpretação de lineamentos, levantamentos geológicos além de perfilagens geofísicas de superfície. Em cada área investigada, foi realizada a perfuração de um poço profundo e aplicadas as técnicas de perfilagens de raios gama, cáliper, flowmeter, imageamento acústico, além da filmagem do poço e realização de ensaios hidráulicos nos dois pontos perfurados. Para caracterização química do aquífero fraturado, foram realizadas coletas de água subterrânea em intervalos selecionados com a utilização de obturadores pneumáticos. As cargas hidráulicas medidas durante a amostragem também auxiliaram no entendimento da direção do fluxo de água. O aquífero cristalino é formado por rochas gnáissicas e se encontra bastante fraturado e intemperizado, principalmente na porção superficial da rocha (até aproximadamente 65,0 m) onde as maiores velocidades de fluxo de água também foram observadas. A rocha sã possui uma menor densidade de fraturas e predominância de minerais mais claros. As fraturas de baixo a médio angulo de mergulho (Grupo 1) são as mais frequentes em ambas as perfurações e possuem direção principal N-S a NE-SW. São observadas, no geral, exercendo grande influência sobre o fluxo de água, principalmente na porção alterada do gnaisse. Fraturas com ângulo elevado de mergulho, classificadas como Grupo 2 (paralelas à foliação) e Grupo 3 (direção NW à W), são também observadas ao longo de toda a perfuração estabelecendo a conexão hidráulica entre as fraturas do Grupo 1. Em menor proporção, são ainda verificadas fraturas com ângulos de mergulho >40 ° pertencente aos Grupos 4 (NE-SW), 5 (E-W), 6 (NW-SE) e 7 (E-W). O fluxo de água subterrânea se mostrou descendente na porção superior da rocha alterada e ascendente na porção mais profunda, possivelmente direcionando a água subterrânea para a região de transição da rocha mais alterada para a rocha sã (entre 61 a 65 m de profundidade). Apesar do fluxo ascendente em profundidade, o bombeamento de poços tubulares existentes no entorno ao longo dos anos, favoreceu a migração dos contaminantes para porções mais profundas. Os contaminantes observados no poço tubular P6 possuem maior semelhança com os contaminantes observados na Área 2, e ambos estão localizados entre lineamentos NW-SE, indicando uma possível influência dos lineamentos no controle sobre o fluxo de água. No entanto, para entendimento do transporte dos contaminantes em área, é necessário um adensamento da rede de monitoramento, levando em consideração a heterogeneidade do meio e as incertezas relacionadas à extrapolação dos dados para áreas não investigadas.

Swimming against the flow: A mechanism of zooplankton aggregation

Zooplankton reside in a constantly flowing environment. However, information about their response to ambient flow has remained elusive, because of the difficulties of following the individual motions of these minute, nearly transparent animals in the ocean. Using a three-dimensional acoustic imaging system, we tracked >375,000 zooplankters at two coastal sites in the Red Sea. Resolution of their motion from that of the water showed that the animals effectively maintained their depth by swimming against upwelling and downwelling currents moving at rates of up to tens of body lengths per second, causing their accumulation at frontal zones. This mechanism explains how oceanic fronts become major feeding grounds for predators and targets for fishermen.

Coding Strategies and Implementations of Compressive Sensing

This dissertation studies the coding strategies of computational imaging to overcome the limitation of conventional sensing techniques. The information capacity of conventional sensing is limited by the physical properties of optics, such as aperture size, detector pixels, quantum efficiency, and sampling rate. These parameters determine the spatial, depth, spectral, temporal, and polarization sensitivity of each imager. To increase sensitivity in any dimension can significantly compromise the others.

This research implements various coding strategies subject to optical multidimensional imaging and acoustic sensing in order to extend their sensing abilities. The proposed coding strategies combine hardware modification and signal processing to exploiting bandwidth and sensitivity from conventional sensors. We discuss the hardware architecture, compression strategies, sensing process modeling, and reconstruction algorithm of each sensing system.

Optical multidimensional imaging measures three or more dimensional information of the optical signal. Traditional multidimensional imagers acquire extra dimensional information at the cost of degrading temporal or spatial resolution. Compressive multidimensional imaging multiplexes the transverse spatial, spectral, temporal, and polarization information on a two-dimensional (2D) detector. The corresponding spectral, temporal and polarization coding strategies adapt optics, electronic devices, and designed modulation techniques for multiplex measurement. This computational imaging technique provides multispectral, temporal super-resolution, and polarization imaging abilities with minimal loss in spatial resolution and noise level while maintaining or gaining higher temporal resolution. The experimental results prove that the appropriate coding strategies may improve hundreds times more sensing capacity.

Human auditory system has the astonishing ability in localizing, tracking, and filtering the selected sound sources or information from a noisy environment. Using engineering efforts to accomplish the same task usually requires multiple detectors, advanced computational algorithms, or artificial intelligence systems. Compressive acoustic sensing incorporates acoustic metamaterials in compressive sensing theory to emulate the abilities of sound localization and selective attention. This research investigates and optimizes the sensing capacity and the spatial sensitivity of the acoustic sensor. The well-modeled acoustic sensor allows localizing multiple speakers in both stationary and dynamic auditory scene; and distinguishing mixed conversations from independent sources with high audio recognition rate.

Acoustic hypersensitivity in adult rats after neonatal ventral hippocampus lesions

Rats with a bilateral neonatal ventral hippocampus lesion (NVHL) are used as models of neurobiological aspects of schizophrenia. In view of their decreased number of GABAergic interneurons, we hypothesized that they would show increased reactivity to acoustic stimuli. We systematically characterized the acoustic reactivity of NVHL rats and sham operated controls. They were behaviourally observed during a loud white noise. A first cohort of 7 months` old rats was studied. Then the observations were reproduced in a second cohort of the same age after characterizing the reactivity of the same rats to dopaminergic drugs. A third cohort of rats was studied at 2, 3, 4, 5 and 6 months. In subsets of lesioned and control rats, inferior colliculus auditory evoked potentials were recorded. A significant proportion of rats (50-62%) showed aberrant audiogenic responses with explosive wild running resembling the initial phase of audiogenic seizures. This was not correlated with their well-known enhanced reactivity to dopaminergic drugs. The proportion of rats showing this strong reaction increased with rats` age. After the cessation of the noise, NVHL rats showed a long freezing period that did neither depend on the size of the lesion nor on the rats` age. The initial negative deflection of the auditory evoked potential was enhanced in the inferior colliculus of only NVHL rats that displayed wild running. Complementary anatomical investigations using X-ray scans in the living animal, and alizarin red staining of brain slices, revealed a thin layer of calcium deposit close to the medial geniculate nuclei in post-NVHL rats, raising the possibility that this may contribute to the hyper-reactivity to sounds seen in these animals. The findings of this study provide complementary information with potential relevance for the hyper-reactivity noted in patients with schizophrenia, and therefore a tool to investigate the underlying biology of this endophenotype. (C) 2009 Elsevier B.V. All rights reserved.

Magnetic resonance imaging of the vocal tract: techniques and applications

Magnetic resonance (MR) imaging has been used to analyse and evaluate the vocal tract shape through different techniques and with promising results in several fields. Our purpose is to demonstrate the relevance of MR and image processing for the vocal tract study. The extraction of contours of the air cavities allowed the set - up of a number of 3D reconstruction image stacks by means of the combination of orthogonally oriented sets of slices for e ach articulatory gesture, as a new approach to solve the expected spatial under sampling of the imaging process. In result these models give improved information for the visualization of morphologic and anatomical aspects and are useful for partial measure ments of the vocal tract shape in different situations. Potential use can be found in Medical and therapeutic applications as well as in acoustic articulatory speech modelling.

Fracture imaging within a granitic rock aquifer using multiple-offset single-hole and cross-hole GPR reflection data

The sparsely spaced highly permeable fractures of the granitic rock aquifer at Stang-er-Brune (Brittany, France) form a well-connected fracture network of high permeability but unknown geometry. Previous work based on optical and acoustic logging together with single-hole and cross-hole flowmeter data acquired in 3 neighbouring boreholes (70-100 m deep) has identified the most important permeable fractures crossing the boreholes and their hydraulic connections. To constrain possible flow paths by estimating the geometries of known and previously unknown fractures, we have acquired, processed and interpreted multifold, single- and cross-hole GPR data using 100 and 250 MHz antennas. The GPR data processing scheme consisting of timezero corrections, scaling, bandpass filtering and F-X deconvolution, eigenvector filtering, muting, pre-stack Kirchhoff depth migration and stacking was used to differentiate fluid-filled fracture reflections from source generated noise. The final stacked and pre-stack depth-migrated GPR sections provide high-resolution images of individual fractures (dipping 30-90°) in the surroundings (2-20 m for the 100 MHz antennas; 2-12 m for the 250 MHz antennas) of each borehole in a 2D plane projection that are of superior quality to those obtained from single-offset sections. Most fractures previously identified from hydraulic testing can be correlated to reflections in the single-hole data. Several previously unknown major near vertical fractures have also been identified away from the boreholes.

Validity of the acoustic approximation in full-waveform seismic crosshole tomography

Acoustic waveform inversions are an increasingly popular tool for extracting subsurface information from seismic data. They are computationally much more efficient than elastic inversions. Naturally, an inherent disadvantage is that any elastic effects present in the recorded data are ignored in acoustic inversions. We investigate the extent to which elastic effects influence seismic crosshole data. Our numerical modeling studies reveal that in the presence of high contrast interfaces, at which P-to-S conversions occur, elastic effects can dominate the seismic sections, even for experiments involving pressure sources and pressure receivers. Comparisons of waveform inversion results using a purely acoustic algorithm on synthetic data that is either acoustic or elastic, show that subsurface models comprising small low-to-medium contrast (?30%) structures can be successfully resolved in the acoustic approximation. However, in the presence of extended high-contrast anomalous bodies, P-to-S-conversions may substantially degrade the quality of the tomographic images. In particular, extended low-velocity zones are difficult to image. Likewise, relatively small low-velocity features are unresolved, even when advanced a priori information is included. One option for mitigating elastic effects is data windowing, which suppresses later arriving seismic arrivals, such as shear waves. Our tests of this approach found it to be inappropriate because elastic effects are also included in earlier arriving wavetrains. Furthermore, data windowing removes later arriving P-wave phases that may provide critical constraints on the tomograms. Finally, we investigated the extent to which acoustic inversions of elastic data are useful for time-lapse analyses of high contrast engineered structures, for which accurate reconstruction of the subsurface structure is not as critical as imaging differential changes between sequential experiments. Based on a realistic scenario for monitoring a radioactive waste repository, we demonstrated that acoustic inversions of elastic data yield substantial distortions of the tomograms and also unreliable information on trends in the velocity changes.