Using GPUs for beamforming acceleration on SAFT imaging

SAFT techniques are based on the sequential activation, in emission and reception, of the array elements and the post-processing of all the received signals to compose the image. Thus, the image generation can be divided into two stages: (1) the excitation and acquisition stage, where the signals received by each element or group of elements are stored; and (2) the beamforming stage, where the signals are combined together to obtain the image pixels. The use of Graphics Processing Units (GPUs), which are programmable devices with a high level of parallelism, can accelerate the computations of the beamforming process, that usually includes different functions such as dynamic focusing, band-pass filtering, spatial filtering or envelope detection. This work shows that using GPU technology can accelerate, in more than one order of magnitude with respect to CPU implementations, the beamforming and post-processing algorithms in SAFT imaging. ©2009 IEEE.

CRS seismic data imaging system: a case study for basib reevaluation

Este trabalho apresenta resultados práticos de uma atenção sistemática dada ao processamento e à interpretação sísmica de algumas linhas terrestres do conjunto de dados do gráben do Tacutu (Brasil), sobre os quais foram aplicadas etapas fundamentais do sistema WIT de imageamento do empilhamento CRS (Superfície de Reflexão Comum) vinculado a dados. Como resultado, esperamos estabelecer um fluxograma para a reavaliação sísmica de bacias sedimentares. Fundamentado nos atributos de frente de onda resultantes do empilhamento CRS, um macro-modelo suave de velocidades foi obtido através de inversão tomográfica. Usando este macro-modelo, foi realizado uma migração à profundidade pré- e pós-empilhamento. Além disso, outras técnicas baseadas no empilhamento CRS foram realizadas em paralelo como correção estática residual e migração de abertura-limitada baseada na zona de Fresnel projetada. Uma interpretação geológica sobre as seções empilhadas e migradas foi esboçada. A partir dos detalhes visuais dos painéis é possível interpretar desconformidades, afinamentos, um anticlinal principal falhado com conjuntos de horstes e grábens. Também, uma parte da linha selecionada precisa de processamento mais detalhado para evidenciar melhor qualquer estrutura presente na subsuperfície.

Inversão de dados sísmicos de reflexão a partir da curva do tempo de trânsito

A presente Dissertação de Mestrado tem como objetivo o estudo do problema de inversão sísmica baseada em refletores planos para arranjo fonte-comum (FC) e ponto-médiocomum (PMC). O modelo direto é descrito por camadas homogêneas, isotrópicas com interfaces plano-horizontais. O problema é relacionado ao empilhamento NMO baseado na otimização da função semblance, para seções PMC corrigidas de sobretempo normal (NMO). O estudo foi baseado em dois princípios. O primeiro princípio adotado foi de combinar dois grupos de métodos de inversão: um Método Global e um Método Local. O segundo princípio adotado foi o de cascata, segundo a teoria Wichert-Herglotz-Bateman, que estabelece que para conhecer uma camada inferior tem-se que conhecer primeiro a camada superior (dissecação). A aplicação do estudo é voltada à simulação sísmica de Bacia Sedimentar do Solimões e de Bacia Marinha para se obter uma distribuição local 1D de velocidades e espessuras para a subsuperfície em horizontes alvo. Sendo assim, limitamos a inversão entre 4 e 11 refletores, uma vez que na prática a indústria limita uma interpretação realizada apenas em número equivalente de 3 a 4 refletores principais. Ressalta-se que este modelo é aplicável como condição inicial ao imageamento de seções sísmicas em regiões geologicamente complexas com variação horizontal suave de velocidades. Os dados sintéticos foram gerados a partir dos modelos relacionados a informações geológicas, o que corresponde a uma forte informação a priori no modelo de inversão. Para a construção dos modelos relacionados aos projetos da Rede Risco Exploratório (FINEP) e de formação de recursos humanos da ANP em andamento, analisamos os seguintes assuntos relevantes: (1) Geologia de bacias sedimentares terrestre dos Solimões e ma rinha (estratigráfica, estrutural, tectônica e petrolífera); (2) Física da resolução vertical e horizontal; e (3) Discretização temporal-espacial no cubo de multi-cobertura. O processo de inversão é dependente do efeito da discretização tempo-espacial do campo de ondas, dos parâmetros físicos do levantamento sísmico, e da posterior reamostragem no cubo de cobertura múltipla. O modelo direto empregado corresponde ao caso do operador do empilhamento NMO (1D), considerando uma topografia de observação plana. O critério básico tomado como referência para a inversão e o ajuste de curvas é a norma 2 (quadrática). A inversão usando o presente modelo simples é computacionalmente atrativa por ser rápida, e conveniente por permitir que vários outros recursos possam ser incluídos com interpretação física lógica; por exemplo, a Zona de Fresnel Projetada (ZFP), cálculo direto da divergência esférica, inversão Dix, inversão linear por reparametrização, informações a priori, regularização. A ZFP mostra ser um conceito út il para estabelecer a abertura da janela espacial da inversão na seção tempo-distância, e representa a influência dos dados na resolução horizontal. A estimativa da ZFP indica uma abertura mínima com base num modelo adotado, e atualizável. A divergência esférica é uma função suave, e tem base física para ser usada na definição da matriz ponderação dos dados em métodos de inversão tomográfica. A necessidade de robustez na inversão pode ser analisada em seções sísmicas (FC, PMC) submetida a filtragens (freqüências de cantos: 5;15;75;85; banda-passante trapezoidal), onde se pode identificar, comparar e interpretar as informações contidas. A partir das seções, concluímos que os dados são contaminados com pontos isolados, o que propõe métodos na classe dos considerados robustos, tendo-se como referência a norma 2 (quadrados- mínimos) de ajuste de curvas. Os algoritmos foram desenvolvidos na linguagem de programação FORTRAN 90/95, usando o programa MATLAB para apresentação de resultados, e o sistema CWP/SU para modelagem sísmica sintética, marcação de eventos e apresentação de resultados.

Buccal Bone Plate Remodeling After Immediate Implant Placement With and Without Synthetic Bone Grafting and Flapless Surgery: Radiographic Study in Dogs

Recent studies in animals have shown pronounced resorption of the buccal bone plate after immediate implantation. The use of flapless surgical procedures prior to the installation of immediate implants, as well as the use of synthetic bone graft in the gaps represent viable alternatives to minimize buccal bone resorption and to favor osseointegration. The aim of this study was to evaluate the healing of the buccal bone plate following immediate implantation using the flapless approach, and to compare this process with sites in which a synthetic bone graft was or was not inserted into the gap between the implant and the buccal bone plate. Lower bicuspids from 8 dogs were bilaterally extracted without the use of flaps, and 4 implants were installed in the alveoli in each side of the mandible and were positioned 2.0 mm from the buccal bone plate (gap). Four groups were devised: 2.0-mm subcrestal implants (3.3 x 8 mm) using bone grafts (SCTG), 2.0-mm subcrestal implants without bone grafts (SCCG), equicrestal implants (3.3 x 10 mm) with bone grafts (EGG), and equicrestal implants without bone grafts (ECCG). One week following the surgical procedures, metallic prostheses were installed, and within 12 weeks the dogs were sacrificed. The blocks containing the individual implants were turned sideways, and radiographic imaging was obtained to analyze the remodeling of the buccal bone plate. In the analysis of the resulting distance between the implant shoulder and the bone crest, statistically significant differences were found in the SCTG when compared to the ECTG (P = .02) and ECCG (P = .03). For mean value comparison of the resulting linear distance between the implant surface and the buccal plate, no statistically significant difference was found among all groups (P > .05). The same result was observed in the parameter for presence or absence of tissue formation between the implant surface and buccal plate. Equicrestally placed implants, in this methodology, presented little or no loss of the buccal bone. The subcrestally positioned implants presented loss of buccal bone, even though synthetic bone graft was used. The buccal bone, however, was always coronal to the implant shoulder.

A comparative thermoluminescence and electron spin resonance study of synthetic carbonated A-type hydroxyapatite

Intensity of the 150 degrees C thermoluminescence peak of beta-irradiated carbonated synthetic A-type hydroxyapatite is approximately 12 times higher than that of the noncarbonated material. Deconvolution of the glow curve showed that this peak is a result of a trap distribution. An attempt was made to relate this thermoluminescence peak enhanced by carbonation with the ESR signal of the CO2- radical in natural or synthetic hydroxyapatite. (C) 2011 Elsevier Ltd. All rights reserved.

Nanoassisted Laser Desorption-Ionization-MS Imaging of Tumors

The ability of nanoassisted laser desorption-ionization mass spectrometry (NALDI-MS) imaging to provide selective chemical monitoring with proper spatial distribution of lipid profiles from tumor tissues after plate imprinting has been tested. NALDI-MS imaging identified and mapped several potential lipid biomarkers in a murine model of melanoma tumor (inoculation of B16/F10 cells). It also confirmed that the in vivo treatment of tumor bearing mice with synthetic supplement containing phosphoethanolamine (PHO-S) promoted an accentuated decrease in relative abundance of the tumor biomarkers. NALDI-MS imaging is a matrix-free LDI protocol based on the selective imprinting of lipids in the NALDI plate followed by the removal of the tissue. It therefore provides good quality and selective chemical images with preservation of spatial distribution and less interference from tissue material. The test case described herein illustrates the potential of chemically selective NALDI-MS imaging for biomarker discovery.

On the development of novel cocaine-analogues for in vivo imaging of the dopamine transporter status

The present thesis is concerned with the development of novel cocaine-derived dopamine transporter ligands for the non-invasive exploration of the striatal and extra-striatal dopamine transporter (DAT) in living systems. The presynaptic dopamine transporter acquires an important function within the mediation of dopaminergic signal transduction. Its availability can serve as a measure for the overall integrity of the dopaminergic system. The DAT is upregulated in early Parkinson’s disease (PD), resulting in an increased availability of DAT-binding sites in the striatal DAT domains. Thereby, DAT imaging has become an important routine diagnostic tool for the early diagnosis of PD in patients, as well as for the differentiation of PD from symptomatically similar medical conditions. Furthermore, the dopaminergic system is involved in a variety of psychiatric diseases. In this regard, DAT-selective imaging agents may provide detailed insights into the scientific understanding of the biochemical background of both, the progress as well as the origins of the symptoms. DAT-imaging may also contribute to the determination of the dopaminergic therapeutic response for a given medication and thereby contribute to more convenient conditions for the patient. From an imaging point of view, the former demands a high availability of the radioactive probe to facilitate broad application of the modality, whereas the latter profits from short-lived probes, suitable for multi-injection studies. Therefore, labelling with longer-lived 18F-fluoride and in particular the generator nuclide 68Ga is worthwhile for clinical routine imaging. In contrast, the introduction of a 11C-label is a prerequisite for detailed scientific studies of neuronal interactions. The development of suitable DAT-ligands for medical imaging has often been complicated by the mixed binding profile of many compounds that that interact with the DAT. Other drawbacks have included high non-specific binding, extensive metabolism and slow accumulation in the DAT-rich brain areas. However, some recent examples have partially overcome the mentioned complications. Based on the structural speciality of these leads, novel ligand structures were designed and successfully synthesised in the present work. A structure activity relationship (SAR) study was conducted wherein the new structural modifications were examined for their influence on DAT-affinity and selectivity. Two of the compounds showed improvements in in vitro affinity for the DAT as well as selectivity versus the serotonin transporter (SERT) and norepinephrine transporter (NET). The main effort was focussed on the high-affinity candidate PR04.MZ, which was subsequently labelled with 18F and 11C in high yield. An initial pharmacological characterisation of PR04.MZ in rodents revealed highly specific binding to the target brain structures. As a result of low non-specific binding, the DAT-rich striatal area was clearly visualised by autoradiography and µPET. Furthermore, the radioactivity uptake into the DAT-rich brain regions was rapid and indicated fast binding equilibrium. No radioactive metabolite was found in the rat brain. [18F]PR04.MZ and [11C]PR04.MZ were compared in the primate brain and the plasma metabolism was studied. It was found that the ligands specifically visualise the DAT in high and low density in the primate brain. The activity uptake was rapid and quantitative evaluation by Logan graphical analysis and simplified reference tissue model was possible after a scanning time of 30 min. These results further reflect the good characteristics of PR04.MZ as a selective ligand of the neuronal DAT. To pursue 68Ga-labelling of the DAT, initial synthetic studies were performed as part of the present thesis. Thereby, a concept for the convenient preparation of novel bifunctional chelators (BFCs) was developed. Furthermore, the suitability of novel 1,4,7-triazacyclononane based N3S3-type BFCs for biomolecule-chelator conjugates of sufficient lipophilicity for the penetration of the blood-brain-barrier was elucidated.

Towards the 3D attenuation imaging of active volcanoes: methods and tests on real and simulated data

The purpose of my PhD thesis has been to face the issue of retrieving a three dimensional attenuation model in volcanic areas. To this purpose, I first elaborated a robust strategy for the analysis of seismic data. This was done by performing several synthetic tests to assess the applicability of spectral ratio method to our purposes. The results of the tests allowed us to conclude that: 1) spectral ratio method gives reliable differential attenuation (dt*) measurements in smooth velocity models; 2) short signal time window has to be chosen to perform spectral analysis; 3) the frequency range over which to compute spectral ratios greatly affects dt* measurements. Furthermore, a refined approach for the application of spectral ratio method has been developed and tested. Through this procedure, the effects caused by heterogeneities of propagation medium on the seismic signals may be removed. The tested data analysis technique was applied to the real active seismic SERAPIS database. It provided a dataset of dt* measurements which was used to obtain a three dimensional attenuation model of the shallowest part of Campi Flegrei caldera. Then, a linearized, iterative, damped attenuation tomography technique has been tested and applied to the selected dataset. The tomography, with a resolution of 0.5 km in the horizontal directions and 0.25 km in the vertical direction, allowed to image important features in the off-shore part of Campi Flegrei caldera. High QP bodies are immersed in a high attenuation body (Qp=30). The latter is well correlated with low Vp and high Vp/Vs values and it is interpreted as a saturated marine and volcanic sediments layer. High Qp anomalies, instead, are interpreted as the effects either of cooled lava bodies or of a CO2 reservoir. A pseudo-circular high Qp anomaly was detected and interpreted as the buried rim of NYT caldera.

Vessel orientation-dependent sensitivity of optoacoustic imaging using a linear array transducer

For clinical optoacoustic imaging, linear probes are preferably used because they allow versatile imaging of the human body with real-time display and free-hand probe guidance. The two-dimensional (2-D) optoacoustic image obtained with this type of probe is generally interpreted as a 2-D cross-section of the tissue just as is common in echo ultrasound. We demonstrate in three-dimensional simulations, phantom experiments, and in vivo mouse experiments that for vascular imaging this interpretation is often inaccurate. The cylindrical blood vessels emit anisotropic acoustic transients, which can be sensitively detected only if the direction of acoustic radiation coincides with the probe aperture. Our results reveal for this reason that the signal amplitude of different blood vessels may differ even if the vessels have the same diameter and initial pressure distribution but different orientation relative to the imaging plane. This has important implications for the image interpretation, for the probe guidance technique, and especially in cases when a quantitative reconstruction of the optical tissue properties is required.

Morphological ripening of fluid inclusions and coupled zone-refining in quartz crystals revealed by cathodoluminescence imaging: Implications for CL-petrography, fluid inclusion analysis and trace-element geothermometry

Cathodoluminescence (CL) studies have previously shown that some secondary fluid inclusions in luminescent quartz are surrounded by dark, non-luminescent patches, resulting from fracture-sealing by late, trace-element-poor quartz. This finding has led to the tacit generalization that all dark CL patches indicate influx of low temperature, late-stage fluids. In this study we have examined natural and synthetic hydrothermal quartz crystals using CL imaging supplemented by in-situ elemental analysis. The results lead us to propose that all natural, liquid-water-bearing inclusions in quartz, whether trapped on former crystal growth surfaces (i.e., of primary origin) or in healed fractures (i.e., of pseudosecondary or secondary origin), are surrounded by three-dimensional, non-luminescent patches. Cross-cutting relations show that the patches form after entrapment of the fluid inclusions and therefore they are not diagnostic of the timing of fluid entrapment. Instead, the dark patches reveal the mechanism by which fluid inclusions spontaneously approach morphological equilibrium and purify their host quartz over geological time. Fluid inclusions that contain solvent water perpetually dissolve and reprecipitate their walls, gradually adopting low-energy euhedral and equant shapes. Defects in the host quartz constitute solubility gradients that drive physical migration of the inclusions over distances of tens of μm (commonly) up to several mm (rarely). Inclusions thus sequester from their walls any trace elements (e.g., Li, Al, Na, Ti) present in excess of equilibrium concentrations, thereby chemically purifying their host crystals in a process analogous to industrial zone refining. Non-luminescent patches of quartz are left in their wake. Fluid inclusions that contain no liquid water as solvent (e.g., inclusions of low-density H2O vapor or other non-aqueous volatiles) do not undergo this process and therefore do not migrate, do not modify their shapes with time, and are not associated with dark-CL zone-refined patches. This new understanding has implications for the interpretation of solids within fluid inclusions (e.g., Ti- and Al-minerals) and for the elemental analysis of hydrothermal and metamorphic quartz and its fluid inclusions by microbeam methods such as LA-ICPMS and SIMS. As Ti is a common trace element in quartz, its sequestration by fluid inclusions and its depletion in zone-refined patches impacts on applications of the Ti-in-quartz geothermometer.

Application of Independent Components Analysis with the JADE algorithm and NIR hyperspectral imaging for revealing food adulteration

In recent years, Independent Components Analysis (ICA) has proven itself to be a powerful signal-processing technique for solving the Blind-Source Separation (BSS) problems in different scientific domains. In the present work, an application of ICA for processing NIR hyperspectral images to detect traces of peanut in wheat flour is presented. Processing was performed without a priori knowledge of the chemical composition of the two food materials. The aim was to extract the source signals of the different chemical components from the initial data set and to use them in order to determine the distribution of peanut traces in the hyperspectral images. To determine the optimal number of independent component to be extracted, the Random ICA by blocks method was used. This method is based on the repeated calculation of several models using an increasing number of independent components after randomly segmenting the matrix data into two blocks and then calculating the correlations between the signals extracted from the two blocks. The extracted ICA signals were interpreted and their ability to classify peanut and wheat flour was studied. Finally, all the extracted ICs were used to construct a single synthetic signal that could be used directly with the hyperspectral images to enhance the contrast between the peanut and the wheat flours in a real multi-use industrial environment. Furthermore, feature extraction methods (connected components labelling algorithm followed by flood fill method to extract object contours) were applied in order to target the spatial location of the presence of peanut traces. A good visualization of the distributions of peanut traces was thus obtained

Modeling diffraction and imaging of laser beams by the mode-expansion method

We present a new method of modeling imaging of laser beams in the presence of diffraction. Our method is based on the concept of first orthogonally expanding the resultant diffraction field (that would have otherwise been obtained by the laborious application of the Huygens diffraction principle) and then representing it by an effective multimodal laser beam with different beam parameters. We show not only that the process of obtaining the new beam parameters is straightforward but also that it permits a different interpretation of the diffraction-caused focal shift in laser beams. All of the criteria that we have used to determine the minimum number of higher-order modes needed to accurately represent the diffraction field show that the mode-expansion method is numerically efficient. Finally, the characteristics of the mode-expansion method are such that it allows modeling of a vast array of diffraction problems, regardless of the characteristics of the incident laser beam, the diffracting element, or the observation plane. (C) 2005 Optical Society of America.

Synthetic Approaches to Flexible Fluorescent Conjugated Polymers

Conjugated polymers (CPs) are intrinsically fluorescent materials that have been used for various biological applications including imaging, sensing, and delivery of biologically active substances. The synthetic control over flexibility and biodegradability of these materials aids the understanding of the structure-function relationships among the photophysical properties, the self-assembly behaviors of the corresponding conjugated polymer nanoparticles (CPNs), and the cellular behaviors of CPNs, such as toxicity, cellular uptake mechanisms, and sub-cellular localization patterns. Synthetic approaches towards two classes of flexible CPs with well-preserved fluorescent properties are described. The synthesis of flexible poly(p-phenylenebutadiynylene)s (PPBs) uses competing Sonogashira and Glaser coupling reactions and the differences in monomer reactivity to incorporate a small amount (~10%) of flexible, non-conjugated linkers into the backbone. The reaction conditions provide limited control over the proportion of flexible monomer incorporation. Improved synthetic control was achieved in a series of flexible poly(p-phenyleneethynylene)s (PPEs) using modified Sonogashira conditions. In addition to controlling the degree of flexibility, the linker provides disruption of backbone conjugation that offers control of the length of conjugated segments within the polymer chain. Therefore, such control also results in the modulation of the photophysical properties of the materials. CPNs fabricated from flexible PPBs are non-toxic to cells, and exhibit subcellular localization patterns clearly different from those observed with non-flexible PPE CPNs. The subcellular localization patterns of the flexible PPEs have not yet been determined, due to the toxicity of the materials, most likely related to the side-chain structure used in this series. The study of the effect of CP flexibility on self-assembly reorganization upon polyanion complexation is presented. Owing to its high rigidity and hydrophobicity, the PPB backbone undergoes reorganization more readily than PPE. The effects are enhanced in the presence of the flexible linker, which enables more efficient π-π stacking of the aromatic backbone segments. Flexibility has minimal effects on the self-assembly of PPEs. Understanding the role of flexibility on the biophysical behaviors of CPNs is key to the successful development of novel efficient fluorescent therapeutic delivery vehicles.

Tissue Equivalent Phantom Design for Optimization of a Coherent Scatter Imaging System

Scatter in medical imaging is typically cast off as image-related noise that detracts from meaningful diagnosis. It is therefore typically rejected or removed from medical images. However, it has been found that every material, including cancerous tissue, has a unique X-ray coherent scatter signature that can be used to identify the material or tissue. Such scatter-based tissue-identification provides the advantage of locating and identifying particular materials over conventional anatomical imaging through X-ray radiography. A coded aperture X-ray coherent scatter spectral imaging system has been developed in our group to classify different tissue types based on their unique scatter signatures. Previous experiments using our prototype have demonstrated that the depth-resolved coherent scatter spectral imaging system (CACSSI) can discriminate healthy and cancerous tissue present in the path of a non-destructive x-ray beam. A key to the successful optimization of CACSSI as a clinical imaging method is to obtain anatomically accurate phantoms of the human body. This thesis describes the development and fabrication of 3D printed anatomical scatter phantoms of the breast and lung.

The purpose of this work is to accurately model different breast geometries using a tissue equivalent phantom, and to classify these tissues in a coherent x-ray scatter imaging system. Tissue-equivalent anatomical phantoms were designed to assess the capability of the CACSSI system to classify different types of breast tissue (adipose, fibroglandular, malignant). These phantoms were 3D printed based on DICOM data obtained from CT scans of prone breasts. The phantoms were tested through comparison of measured scatter signatures with those of adipose and fibroglandular tissue from literature. Tumors in the phantom were modeled using a variety of biological tissue including actual surgically excised benign and malignant tissue specimens. Lung based phantoms have also been printed for future testing. Our imaging system has been able to define the location and composition of the various materials in the phantom. These phantoms were used to characterize the CACSSI system in terms of beam width and imaging technique. The result of this work showed accurate modeling and characterization of the phantoms through comparison of the tissue-equivalent form factors to those from literature. The physical construction of the phantoms, based on actual patient anatomy, was validated using mammography and computed tomography to visually compare the clinical images to those of actual patient anatomy.

A fast parallel hyperspectral coded aperture algorithm for compressive sensing using OpenCL

In this paper, we develop a fast implementation of an hyperspectral coded aperture (HYCA) algorithm on different platforms using OpenCL, an open standard for parallel programing on heterogeneous systems, which includes a wide variety of devices, from dense multicore systems from major manufactures such as Intel or ARM to new accelerators such as graphics processing units (GPUs), field programmable gate arrays (FPGAs), the Intel Xeon Phi and other custom devices. Our proposed implementation of HYCA significantly reduces its computational cost. Our experiments have been conducted using simulated data and reveal considerable acceleration factors. This kind of implementations with the same descriptive language on different architectures are very important in order to really calibrate the possibility of using heterogeneous platforms for efficient hyperspectral imaging processing in real remote sensing missions.