FPGA-based architecture for hyperspectral endmember extraction

Hyperspectral instruments have been incorporated in satellite missions, providing data of high spectral resolution of the Earth. This data can be used in remote sensing applications, such as, target detection, hazard prevention, and monitoring oil spills, among others. In most of these applications, one of the requirements of paramount importance is the ability to give real-time or near real-time response. Recently, onboard processing systems have emerged, in order to overcome the huge amount of data to transfer from the satellite to the ground station, and thus, avoiding delays between hyperspectral image acquisition and its interpretation. For this purpose, compact reconfigurable hardware modules, such as field programmable gate arrays (FPGAs) are widely used. This paper proposes a parallel FPGA-based architecture for endmember’s signature extraction. This method based on the Vertex Component Analysis (VCA) has several advantages, namely it is unsupervised, fully automatic, and it works without dimensionality reduction (DR) pre-processing step. The architecture has been designed for a low cost Xilinx Zynq board with a Zynq-7020 SoC FPGA based on the Artix-7 FPGA programmable logic and tested using real hyperspectral data sets collected by the NASA’s Airborne Visible Infra-Red Imaging Spectrometer (AVIRIS) over the Cuprite mining district in Nevada. Experimental results indicate that the proposed implementation can achieve real-time processing, while maintaining the methods accuracy, which indicate the potential of the proposed platform to implement high-performance, low cost embedded systems, opening new perspectives for onboard hyperspectral image processing.

Signal subspace identification in hyperspectral imagery

Terrestrial remote sensing imagery involves the acquisition of information from the Earth's surface without physical contact with the area under study. Among the remote sensing modalities, hyperspectral imaging has recently emerged as a powerful passive technology. This technology has been widely used in the fields of urban and regional planning, water resource management, environmental monitoring, food safety, counterfeit drugs detection, oil spill and other types of chemical contamination detection, biological hazards prevention, and target detection for military and security purposes [2-9]. Hyperspectral sensors sample the reflected solar radiation from the Earth surface in the portion of the spectrum extending from the visible region through the near-infrared and mid-infrared (wavelengths between 0.3 and 2.5 µm) in hundreds of narrow (of the order of 10 nm) contiguous bands [10]. This high spectral resolution can be used for object detection and for discriminating between different objects based on their spectral xharacteristics [6]. However, this huge spectral resolution yields large amounts of data to be processed. For example, the Airbone Visible/Infrared Imaging Spectrometer (AVIRIS) [11] collects a 512 (along track) X 614 (across track) X 224 (bands) X 12 (bits) data cube in 5 s, corresponding to about 140 MBs. Similar data collection ratios are achieved by other spectrometers [12]. Such huge data volumes put stringent requirements on communications, storage, and processing. The problem of signal sbspace identification of hyperspectral data represents a crucial first step in many hypersctral processing algorithms such as target detection, change detection, classification, and unmixing. The identification of this subspace enables a correct dimensionality reduction (DR) yelding gains in data storage and retrieval and in computational time and complexity. Additionally, DR may also improve algorithms performance since it reduce data dimensionality without losses in the useful signal components. The computation of statistical estimates is a relevant example of the advantages of DR, since the number of samples required to obtain accurate estimates increases drastically with the dimmensionality of the data (Hughes phnomenon) [13].

Quantificação por imagem ponderada em difusão (DWI) das lesões tumorais da mama

Introdução – A técnica de Difusão por Ressonância Magnética (RM), ao avaliar o movimento das moléculas de água nos tecidos, permite inferir sobre a arquitetura dos mesmos, em particular relativamente à celularidade, volume celular e permeabilidade das membranas. O Coeficiente de Difusão Aparente (ADC) é um parâmetro quantificável da imagem ponderada em difusão (DWI). A sua análise poderá fornecer informação clínica adicional sobre estas lesões, sobretudo em relação à sua caracterização histológica. Objetivos – Caracterizar e diferenciar tipos e alguns subtipos de lesões mamárias através da análise do ADC. Metodologia – 20 Mulheres com 23 lesões mamárias foram submetidas a RM mamária: 3 lesões benignas (3 Fibroadenomas-FA) e 20 malignas (16 Carcinomas Ductais Invasivos-CDI, 2 Carcinomas Ductais In Situ-CDIS e 2 Carcinomas Invasivos Lobulares-CLI). Num equipamento 1.5T aplicou-se uma sequência de Difusão (b=0,50,250,500,750,1000 s/mm2). Obteve-se o ADC através do ajuste exponencial da intensidade de sinal das lesões em função do valor de b, fazendo-se corresponder os valores de ADC à respetiva caracterização histológica e posterior comparação com a literatura. Resultados e Discussão – As lesões malignas apresentaram ADCs significativamente (p=0,014) inferiores [(0,94±0,22)x10-3 mm2/s] aos das benignas [(1,43±0,25)x10-3 mm2/s]. A justificação pode residir no aumento da celularidade e consequente restrição da Difusão que se observa nas lesões malignas. Os CDI apresentaram ADCs baixos [(0,88±0,17)x10-3 mm2/s], enquanto que os CDIS apresentaram ADCs mais elevados [(1,33±0,29)x10-3 mm2/s]. Estes resultados estão de acordo com o facto dos CDIS estarem limitados aos ductos mamários, mantendo-se menos alterada a estrutura do tecido adjacente e resultando numa menor restrição à difusão que nos CDI. Verificaram-se diferenças significativas entre FA e CDI (p=0,010) e entre CDI e CDIS (p=0,049). Conclusões – O ADC possibilita a diferenciação entre lesões mamárias benignas e malignas, bem como entre alguns tipos histológicos. O desenvolvimento deste conceito pode representar um avanço no papel da RM na avaliação destas neoplasias. ABSTRACT - Introduction – The Magnetic Resonance (MR) diffusion technique measures the movement of water molecules in tissues. Therefore, it provides useful information about tissue architecture, specially regarding tissue cellularity, cell volume and membrane permeability. The quantification of diffusion weighted imaging (DWI) data is done by measuring the so-called. Apparent Diffusion Coefficient (ADC). This parameter provides additional clinical information about breast lesions, and potentially allows for in-vivo histological characterization. Purpose – To characterize and differentiate breast lesions through ADC analysis. Methodology – The study comprised 20 women, with 23 breast lesions: 3 benign lesions - 3 Fibroadenomas (FA); and 20 malignant - 16 Invasive Ductal Carcinomas (CDI), 2 Ductal Carcinomas In Situ (CDIS), 2 Invasive Lobular Carcinoma (CLI). On a 1.5T equipment a diffusion-weighted sequence with 6 b-values (b=0,50,250,500,750,1000 s/mm2) was used to examine the patients. ADC was obtained by fitting an exponential to data of lesion signal intensity vs. b values. A correspondence of ADC values to histological lesion characterization was done and finally, the results were comparison with the literature. Results and Discussion – Malignant lesions showed inferior ADCs significantly (p=0.014) lower ((0.94±0.22)x10-3 mm2/s) than the benign lesions ((1.43±0.25)x10-3 mm2/s). This may be associated to increasead cellularity in malignant lesions, which result in higher tissue restriction to diffusion. CDI showed low ADC values ((0.88±0.17)x10-3 mm2/s), while the CDIS showed higher ADCs ((1.33±0.29)x10-3 mm2/s). These results agree with the fact that CDIS are limited to mammary ducts, maintaining a less altered neighboring tissue structure, which results in a lower restriction to diffusion than observed in CDI. Significant differences between FA and CDI (p=0.010) and between CDI and CDIS (p=0.049) were observed. Conclusion – The ADC parameter is able to differentiate between malignant and benign breast lesions, as well as between some histological types.

Lesion detection performance in an anthropomorphic chest phantom: comparative analysis of low-dose CT data on two hybrid imaging systems

Introduction: multimodality environment; requirement for greater understanding of the imaging technologies used, the limitations of these technologies, and how to best interpret the results; dose optimization; introduction of new techniques; current practice and best practice; incidental findings, in low-dose CT images obtained as part of the hybrid imaging process, are an increasing phenomenon with advancing CT technology; resultant ethical and medico-legal dilemmas; understanding limitations of these procedures important when reporting images and recommending follow-up; free-response observer performance study was used to evaluate lesion detection in low-dose CT images obtained during attenuation correction acquisitions for myocardial perfusion imaging, on two hybrid imaging systems.

Lesion detection performance: comparative analysis of low-dose CT data of the chest on two hybrid imaging systems

Incidental findings on low-dose CT images obtained during hybrid imaging are an increasing phenomenon as CT technology advances. Understanding the diagnostic value of incidental findings along with the technical limitations is important when reporting image results and recommending follow-up, which may result in an additional radiation dose from further diagnostic imaging and an increase in patient anxiety. This study assessed lesions incidentally detected on CT images acquired for attenuation correction on two SPECT/CT systems. Methods: An anthropomorphic chest phantom containing simulated lesions of varying size and density was imaged on an Infinia Hawkeye 4 and a Symbia T6 using the low-dose CT settings applied for attenuation correction acquisitions in myocardial perfusion imaging. Twenty-two interpreters assessed 46 images from each SPECT/CT system (15 normal images and 31 abnormal images; 41 lesions). Data were evaluated using a jackknife alternative free-response receiver-operating-characteristic analysis (JAFROC). Results: JAFROC analysis showed a significant difference (P < 0.0001) in lesion detection, with the figures of merit being 0.599 (95% confidence interval, 0.568, 0.631) and 0.810 (95% confidence interval, 0.781, 0.839) for the Infinia Hawkeye 4 and Symbia T6, respectively. Lesion detection on the Infinia Hawkeye 4 was generally limited to larger, higher-density lesions. The Symbia T6 allowed improved detection rates for midsized lesions and some lower-density lesions. However, interpreters struggled to detect small (5 mm) lesions on both image sets, irrespective of density. Conclusion: Lesion detection is more reliable on low-dose CT images from the Symbia T6 than from the Infinia Hawkeye 4. This phantom-based study gives an indication of potential lesion detection in the clinical context as shown by two commonly used SPECT/CT systems, which may assist the clinician in determining whether further diagnostic imaging is justified.

Fluorescence Microscopy Imaging Denoising with Log-Euclidean Priors and Photobleaching Compensation

Fluorescent protein microscopy imaging is nowadays one of the most important tools in biomedical research. However, the resulting images present a low signal to noise ratio and a time intensity decay due to the photobleaching effect. This phenomenon is a consequence of the decreasing on the radiation emission efficiency of the tagging protein. This occurs because the fluorophore permanently loses its ability to fluoresce, due to photochemical reactions induced by the incident light. The Poisson multiplicative noise that corrupts these images, in addition with its quality degradation due to photobleaching, make long time biological observation processes very difficult. In this paper a denoising algorithm for Poisson data, where the photobleaching effect is explicitly taken into account, is described. The algorithm is designed in a Bayesian framework where the data fidelity term models the Poisson noise generation process as well as the exponential intensity decay caused by the photobleaching. The prior term is conceived with Gibbs priors and log-Euclidean potential functions, suitable to cope with the positivity constrained nature of the parameters to be estimated. Monte Carlo tests with synthetic data are presented to characterize the performance of the algorithm. One example with real data is included to illustrate its application.

Molecular imaging agents for detection of β-amyloid plaques in Alzheimer’s disease

The formation of amyloid structures is a neuropathological feature that characterizes several neurodegenerative disorders, such as Alzheimer´s and Parkinson´s disease. Up to now, the definitive diagnosis of these diseases can only be accomplished by immunostaining of post mortem brain tissues with dyes such Thioflavin T and congo red. Aiming at early in vivo diagnosis of Alzheimer´s disease (AD), several amyloid-avid radioprobes have been developed for b-amyloid imaging by positron emission tomography (PET) and single-photon emission computed tomography (SPECT). The aim of this paper is to present a perspective of the available amyloid imaging agents, special those that have been selected for clinical trials and are at the different stages of the US Food and Drugs Administration (FDA) approval.

Classification and staging of chronic liver disease from multimodal data

Chronic liver disease (CLD) is most of the time an asymptomatic, progressive, and ultimately potentially fatal disease. In this study, an automatic hierarchical procedure to stage CLD using ultrasound images, laboratory tests, and clinical records are described. The first stage of the proposed method, called clinical based classifier (CBC), discriminates healthy from pathologic conditions. When nonhealthy conditions are detected, the method refines the results in three exclusive pathologies in a hierarchical basis: 1) chronic hepatitis; 2) compensated cirrhosis; and 3) decompensated cirrhosis. The features used as well as the classifiers (Bayes, Parzen, support vector machine, and k-nearest neighbor) are optimally selected for each stage. A large multimodal feature database was specifically built for this study containing 30 chronic hepatitis cases, 34 compensated cirrhosis cases, and 36 decompensated cirrhosis cases, all validated after histopathologic analysis by liver biopsy. The CBC classification scheme outperformed the nonhierachical one against all scheme, achieving an overall accuracy of 98.67% for the normal detector, 87.45% for the chronic hepatitis detector, and 95.71% for the cirrhosis detector.

Data mining framework for fatty liver disease classification in ultrasound: a hybrid feature extraction paradigm

PURPOSE: Fatty liver disease (FLD) is an increasing prevalent disease that can be reversed if detected early. Ultrasound is the safest and ubiquitous method for identifying FLD. Since expert sonographers are required to accurately interpret the liver ultrasound images, lack of the same will result in interobserver variability. For more objective interpretation, high accuracy, and quick second opinions, computer aided diagnostic (CAD) techniques may be exploited. The purpose of this work is to develop one such CAD technique for accurate classification of normal livers and abnormal livers affected by FLD. METHODS: In this paper, the authors present a CAD technique (called Symtosis) that uses a novel combination of significant features based on the texture, wavelet transform, and higher order spectra of the liver ultrasound images in various supervised learning-based classifiers in order to determine parameters that classify normal and FLD-affected abnormal livers. RESULTS: On evaluating the proposed technique on a database of 58 abnormal and 42 normal liver ultrasound images, the authors were able to achieve a high classification accuracy of 93.3% using the decision tree classifier. CONCLUSIONS: This high accuracy added to the completely automated classification procedure makes the authors' proposed technique highly suitable for clinical deployment and usage.

Chronic liver disease staging classification based on ultrasound, clinical and laboratorial data

In this work the identification and diagnosis of various stages of chronic liver disease is addressed. The classification results of a support vector machine, a decision tree and a k-nearest neighbor classifier are compared. Ultrasound image intensity and textural features are jointly used with clinical and laboratorial data in the staging process. The classifiers training is performed by using a population of 97 patients at six different stages of chronic liver disease and a leave-one-out cross-validation strategy. The best results are obtained using the support vector machine with a radial-basis kernel, with 73.20% of overall accuracy. The good performance of the method is a promising indicator that it can be used, in a non invasive way, to provide reliable information about the chronic liver disease staging.

Diffuse liver disease classification from ultrasound surface characterization, clinical and laboratorial data

In this work liver contour is semi-automatically segmented and quantified in order to help the identification and diagnosis of diffuse liver disease. The features extracted from the liver contour are jointly used with clinical and laboratorial data in the staging process. The classification results of a support vector machine, a Bayesian and a k-nearest neighbor classifier are compared. A population of 88 patients at five different stages of diffuse liver disease and a leave-one-out cross-validation strategy are used in the classification process. The best results are obtained using the k-nearest neighbor classifier, with an overall accuracy of 80.68%. The good performance of the proposed method shows a reliable indicator that can improve the information in the staging of diffuse liver disease.

An experimental study to compare the interface pressure and experience of healthy participants when lying still for 20 minutes in a supine position on two different imaging surfaces

Introduction: Pressure ulcers are a high cost, high volume issue for health and medical care providers, affecting patients’ recovery and psychological wellbeing. The current research of support surfaces on pressure as a risk factor in the development of pressure ulcers is not relevant to the specialised, controlled environment of the radiological setting. Method: 38 healthy participants aged 19-51 were placed supine on two different imaging surfaces. The XSENSOR pressure mapping system was used to measure the interface pressure. Data was acquired over a time of 20 minutes preceded by 6 minutes settling time to reduce measurement error. Qualitative information regarding participants’ opinion on pain and comfort was recorded using a questionnaire. Data analysis was performed using SPSS 22. Results: Data was collected from 30 participants aged 19 to 51 (mean 25.77, SD 7.72), BMI from 18.7 to 33.6 (mean 24.12, SD 3.29), for two surfaces, following eight participant exclusions due to technical faults. Total average pressure, average pressure for jeopardy areas (head, sacrum & heels) and peak pressure for jeopardy areas were calculated as interface pressure in mmHg. Qualitative data showed that a significant difference in experiences of comfort and pain was found in the jeopardy areas (P<0.05) between the two surfaces. Conclusion: A significant difference is seen in average pressure between the two surfaces. Pain and comfort data also show a significant difference between the surfaces, both findings support the proposal for further investigation into the effects of radiological surfaces as a risk factor for the formation of pressure ulcers.

Constrains on the structure of Maio Island (Cape Verde) by a three-dimensional gravity model: Imaging partially exhumed magma chambers

We propose a 3-D gravity model for the volcanic structure of the island of Maio (Cape Verde archipelago) with the objective of solving some open questions concerning the geometry and depth of the intrusive Central Igneous Complex. A gravity survey was made covering almost the entire surface of the island. The gravity data was inverted through a non-linear 3-D approach which provided a model constructed in a random growth process. The residual Bouguer gravity field shows a single positive anomaly presenting an elliptic shape with a NWSE trending long axis. This Bouguer gravity anomaly is slightly off-centred with the island but its outline is concordant with the surface exposure of the Central Igneous Complex. The gravimetric modelling shows a high-density volume whose centre of mass is about 4500 m deep. With increasing depth, and despite the restricted gravimetric resolution, the horizontal sections of the model suggest the presence of two distinct bodies, whose relative position accounts for the elongated shape of the high positive Bouguer gravity anomaly. These bodies are interpreted as magma chambers whose coeval volcanic counterparts are no longer preserved. The orientation defined by the two bodies is similar to that of other structures known in the southern group of the Cape Verde islands, thus suggesting a possible structural control constraining the location of the plutonic intrusions.

Visual function assessment in medical imaging research

Background - Medical image perception research relies on visual data to study the diagnostic relationship between observers and medical images. A consistent method to assess visual function for participants in medical imaging research has not been developed and represents a significant gap in existing research. Methods - Three visual assessment factors appropriate to observer studies were identified: visual acuity, contrast sensitivity, and stereopsis. A test was designed for each, and 30 radiography observers (mean age 31.6 years) participated in each test. Results - Mean binocular visual acuity for distance was 20/14 for all observers. The difference between observers who did and did not use corrective lenses was not statistically significant (P = .12). All subjects had a normal value for near visual acuity and stereoacuity. Contrast sensitivity was better than population norms. Conclusion - All observers had normal visual function and could participate in medical imaging visual analysis studies. Protocols of evaluation and populations norms are provided. Further studies are necessary to understand fully the relationship between visual performance on tests and diagnostic accuracy in practice.