Standardization of near infrared spectroscopies via sample spectral correlation equalization.

A novel method for near-infrared (NIR) spectroscopy spectra standardization is presented. NIR spectroscopies have been widely used in analytical chemistry, and many methods have been developed for NIR spectra standardization. To establish a robust standardization transformation, most existing methods require spectral data sets from both primal and secondary instruments for 1-1 correspondence validation. However, this limits the usage of standardization methods. This paper investigates an interesting issue, "Can spectra data in sets be arbitrarily order?" and further develops a completely different approach from existing methods in view of statistical signal processing. The key idea is to first compensate for the distortion along the wavelength and intensity of the spectra, and then transfer the second order statistic (2OS) from the primal spectra to the secondary spectra via data sphering and an inverse sphering transform so that the 2OS can be estimated regardless of the sample statistic order. To further demonstrate how the developed method can extend the usage of the NIR spectra standardization, several application-driven experiments on classification and regression are conducted for demonstration, and a comparison to the piecewise direct standardization (PDS) is also studied.

Breast Tumor Detection and Classification Using Intravoxel Incoherent Motion Hyperspectral Imaging Techniques.

Breast cancer is a main cause of disease and death for women globally. Because of the limitations of traditional mammography and ultrasonography, magnetic resonance imaging (MRI) has gradually become an important radiological method for breast cancer assessment over the past decades. MRI is free of the problems related to radiation exposure and provides excellent image resolution and contrast. However, a disadvantage is the injection of contrast agent, which is toxic for some patients (such as patients with chronic renal disease or pregnant and lactating women). Recent findings of gadolinium deposits in the brain are also a concern. To address these issues, this paper develops an intravoxel incoherent motion- (IVIM-) MRI-based histogram analysis approach, which takes advantage of several hyperspectral techniques, such as the band expansion process (BEP), to expand a multispectral image to hyperspectral images and create an automatic target generation process (ATGP). After automatically finding suspected targets, further detection was attained by using kernel constrained energy minimization (KCEM). A decision tree and histogram analysis were applied to classify breast tissue via quantitative analysis for detected lesions, which were used to distinguish between three categories of breast tissue: malignant tumors (i.e., central and peripheral zone), cysts, and normal breast tissues. The experimental results demonstrated that the proposed IVIM-MRI-based histogram analysis approach can effectively differentiate between these three breast tissue types.

Anomaly Detection Outperforms Logistic Regression in Predicting Outcomes in Trauma Patients.

OBJECTIVE: Recent advancements in trauma resuscitation have shown a great benefit of early identification and control of hemorrhage, which is the most common cause of death in injured patients. We introduce a new analytical approach, anomaly detection (AD), as an alternative method to the traditional logistic regression (LR) method in predicting which injured patients receive transfusions, intensive care, and other interventions. METHODS: We abstracted routinely collected prehospital vital sign data from patient records (adult patients who survived more than 15 minutes after being directly admitted to a level 1 trauma center). The vital signs of the study cohort were analyzed using both LR and AD methods. Predictions on blood transfusions generated by these approaches were compared with hospital records using the respective areas under the receiver operating characteristic curves (AUROC). RESULTS: Of the patients seen at our trauma center between January 1, 2009, and December 31, 2010, 5,464 were included. AD significantly outperformed LR, identifying which patients would receive transfusions of uncrossmatched blood, transfusion of blood between the time of admission and 6 hours later, the need for intensive care, and in-hospital mortality (mean AUROC = 0.764 and 0.720, respectively). AD and LR provided similar predictions for the patients who would receive massive transfusion. Under the stratified 10 fold times 10 cross-validation test, AD also had significantly lower AUROC variance across subgroups than LR, suggesting AD is a more stable predictions model. CONCLUSIONS: AD provides enhanced predictions for clinically relevant outcomes in the trauma patient cohort studied and may assist providers in caring for acutely injured patients in the prehospital arena.

Classification of Tree Species in Overstorey Canopy of Subtropical Forest Using QuickBird Images.

This paper proposes a supervised classification scheme to identify 40 tree species (2 coniferous, 38 broadleaf) belonging to 22 families and 36 genera in high spatial resolution QuickBird multispectral images (HMS). Overall kappa coefficient (OKC) and species conditional kappa coefficients (SCKC) were used to evaluate classification performance in training samples and estimate accuracy and uncertainty in test samples. Baseline classification performance using HMS images and vegetation index (VI) images were evaluated with an OKC value of 0.58 and 0.48 respectively, but performance improved significantly (up to 0.99) when used in combination with an HMS spectral-spatial texture image (SpecTex). One of the 40 species had very high conditional kappa coefficient performance (SCKC ≥ 0.95) using 4-band HMS and 5-band VIs images, but, only five species had lower performance (0.68 ≤ SCKC ≤ 0.94) using the SpecTex images. When SpecTex images were combined with a Visible Atmospherically Resistant Index (VARI), there was a significant improvement in performance in the training samples. The same level of improvement could not be replicated in the test samples indicating that a high degree of uncertainty exists in species classification accuracy which may be due to individual tree crown density, leaf greenness (inter-canopy gaps), and noise in the background environment (intra-canopy gaps). These factors increase uncertainty in the spectral texture features and therefore represent potential problems when using pixel-based classification techniques for multi-species classification.

[Orthogonal Vector Projection Algorithm for Spectral Unmixing].

Spectrum unmixing is an important part of hyperspectral technologies, which is essential for material quantity analysis in hyperspectral imagery. Most linear unmixing algorithms require computations of matrix multiplication and matrix inversion or matrix determination. These are difficult for programming, especially hard for realization on hardware. At the same time, the computation costs of the algorithms increase significantly as the number of endmembers grows. Here, based on the traditional algorithm Orthogonal Subspace Projection, a new method called. Orthogonal Vector Projection is prompted using orthogonal principle. It simplifies this process by avoiding matrix multiplication and inversion. It firstly computes the final orthogonal vector via Gram-Schmidt process for each endmember spectrum. And then, these orthogonal vectors are used as projection vector for the pixel signature. The unconstrained abundance can be obtained directly by projecting the signature to the projection vectors, and computing the ratio of projected vector length and orthogonal vector length. Compared to the Orthogonal Subspace Projection and Least Squares Error algorithms, this method does not need matrix inversion, which is much computation costing and hard to implement on hardware. It just completes the orthogonalization process by repeated vector operations, easy for application on both parallel computation and hardware. The reasonability of the algorithm is proved by its relationship with Orthogonal Sub-space Projection and Least Squares Error algorithms. And its computational complexity is also compared with the other two algorithms', which is the lowest one. At last, the experimental results on synthetic image and real image are also provided, giving another evidence for effectiveness of the method.

Adrenal gland abnormality detection using random forest classification.

Adrenal abnormalities are commonly identified on computed tomography (CT) and are seen in at least 5 % of CT examinations of the thorax and abdomen. Previous studies have suggested that evaluation of Hounsfield units within a region of interest or a histogram analysis of a region of interest can be used to determine the likelihood that an adrenal gland is abnormal. However, the selection of a region of interest can be arbitrary and operator dependent. We hypothesize that segmenting the entire adrenal gland automatically without any human intervention and then performing a histogram analysis can accurately detect adrenal abnormality. We use the random forest classification framework to automatically perform a pixel-wise classification of an entire CT volume (abdomen and pelvis) into three classes namely right adrenal, left adrenal, and background. Once we obtain this classification, we perform histogram analysis to detect adrenal abnormality. The combination of these methods resulted in a sensitivity and specificity of 80 and 90 %, respectively, when analyzing 20 adrenal glands seen on volumetric CT datasets for abnormality.

Random N-finder (N-FINDR) endmember extraction algorithms for hyperspectral imagery.

N-finder algorithm (N-FINDR) has been widely used in endmember extraction. When it comes to implementation several issues need to be addressed. One is determination of endmembers, p required for N-FINDR to generate. Another is its computational complexity resulting from an exhaustive search. A third one is its requirement of dimensionality reduction. A fourth and probably the most critical issue is its use of random initial endmembers which results in inconsistent final endmember selection and results are not reproducible. This paper re-invents the wheel by re-designing the N-FINDR in such a way that all the above-mentioned issues can be resolved while making the last issue an advantage. The idea is to implement the N-FINDR as a random algorithm, called random N-FINDR (RN-FINDR) so that a single run using one set of random initial endmembers is considered as one realization. If there is an endmember present in the data, it should appear in any realization regardless of what random set of initial endmembers is used. In this case, the N-FINDR is terminated when the intersection of all realizations produced by two consecutive runs of RN-FINDR remains the same in which case the p is then automatically determined by the intersection set without appealing for any criterion. In order to substantiate the proposed RN-FINDR custom-designed synthetic image experiments with complete knowledge are conducted for validation and real image experiments are also performed to demonstrate its utility in applications.

Quantitative analysis in clinical applications of brain MRI using independent component analysis coupled with support vector machine.

PURPOSE: To effectively perform quantification of brain normal tissues and pathologies simultaneously, independent component analysis (ICA) coupled with support vector machine (SVM) is investigated and evaluated for effective volumetric measurements of normal and lesion tissues using multispectral MR images. MATERIALS AND METHODS: Synthetic and real MR data of normal brain and white matter lesion (WML) data were used to evaluate the accuracy and reproducibility of gray matter (GM), white matter (WM), and WML volume measurements by using the proposed ICA+SVM method to analyze three sets of MR images, T1-weighted, T2-weighted, and proton density/fluid-attenuated inversion recovery images. RESULTS: The Tanimoto indexes of GM/WM classification in the normal synthetic data calculated by the ICA+SVM method were 0.82/0.89 for data with 0% noise level. As for clinical MR data experiments, the ICA+SVM method clearly extracted the normal tissues and white matter hyperintensity lesions from the MR images, with low intra- and inter-operator coefficient of variations. CONCLUSION: The experiments conducted provide evidence that the ICA+SVM method has shown promise and potential in applications to classification of normal and pathological tissues in brain MRI.

An automatic computer-aided detection system for meniscal tears on magnetic resonance images.

Knee-related injuries including meniscal tears are common in both young athletes and the aging population, and require accurate diagnosis and surgical intervention when appropriate. With proper techniques and radiologists' experienced skills, confidence in detection of meniscal tears can be quite high. This paper develops a novel computer-aided detection (CAD) diagnostic system for automatic detection of meniscal tears in the knee. Evaluation of this CAD system using an archived database of images from 40 individuals with suspected knee injuries indicates that the sensitivity and specificity of the proposed CAD system are 83.87% and 75.19%, respectively, compared to the mean sensitivity and specificity of 77.41% and 81.39%, respectively, obtained by experienced radiologists in routine diagnosis without using the CAD. The experimental results suggest that the developed CAD system has great potential and promise in automatic detection of both simple and complex meniscal tears of the knee.

Band expansion-based over-complete independent component analysis for multispectral processing of magnetic resonance images.

Independent component analysis (ICA) has found great promise in magnetic resonance (MR) image analysis. Unfortunately, two key issues have been overlooked and not investigated. One is the lack of MR images to be used to unmix signal sources of interest. Another is the use of random initial projection vectors by ICA, which causes inconsistent results. In order to address the first issue, this paper introduces a band-expansion process (BEP) to generate an additional new set of images from the original MR images via nonlinear functions. These newly generated images are then combined with the original MR images to provide sufficient MR images for ICA analysis. In order to resolve the second issue, a prioritized ICA (PICA) is designed to rank the ICA-generated independent components (ICs) so that MR brain tissue substances can be unmixed and separated by different ICs in a prioritized order. Finally, BEP and PICA are combined to further develop a new ICA-based approach, referred to as PICA-BEP to perform MR image analysis.

3D localization of clustered microcalcifications using cranio-caudal and medio-lateral oblique views.

This paper presents a 3D localization method to register clustered microcalcifications on mammograms from cranio-caudal (CC) and medio-lateral oblique (MLO) views. The method consists of three major components: registration of clustered microcalcifications in CC and MLO views, 3D localization of clustered microcalcifications and 3D visualization of clustered microcalcifications. The registration is performed based on three features, gradient, energy and local entropy codes that are independent of spatial locations of microcalcifications in two different views and are prioritized by discriminability in a binary decision tree. The 3D localization is determined by a sequence of coordinate corrections of calcified pixels using the breast nipple as a controlling point. Finally, the 3D visualization implements a virtual reality modeling language viewer (VRMLV) to view the exact location of the lesion as a guide for needle biopsy. In order to validate our proposed 3D localization system, a set of breast lesions, which appear both in mammograms and in MR Images is used for experiments where the depth of clustered microcalcifications can be verified by the MR images.

New texture shape feature coding-based computer aided diagnostic methods for classification of masses on mammograms.

This work presents new texture shape feature coding (TSFC)-based computer aided diagnostic (CAD) classification methods for classification of masses on mammograms. It introduces a new concept of 1.5-order 3-neighbor 3x3 connectivity to extract texture shape features that can describe multiples of 22.5 degrees . In order to effectively utilize these shape features, two new methods of implementing TFSC are further proposed to convert these features to texture feature numbers (TFNs), TFNs in quaternary (TFNq) which expresses a TFN in quaternary expansion and TFNs in product (TFNx ) that represents a TFN in terms of a product. Both TFNq and TFNx can then produce texture shape histograms in the same way that a gray-level histogram is generated for an image. Such a texture shape histogram is further used to generate various shape features of masses on mammograms for classification. In order to demonstrate the promise of our TSFC-based CAD methods, the MiniMammographic database provided by the Mammographic Image Analysis Society (MIAS) is used for experiments.

Detection of spectral signatures in multispectral MR images for classification.

This paper presents a new spectral signature detection approach to magnetic resonance (MR) image classification. It is called constrained energy minimization (CEM) method, which is derived from the minimum variance distortionless response in passive sensor array processing. It considers a bank of spectral channels as an array of sensors where each spectral channel represents a sensor and object spectral signature in multispectral MR images are viewed as signals impinging upon the array. The strength of the CEM lies on its ability in detection of spectral signatures of interest without knowing image background. The detected spectral signatures are then used for classification. The CEM makes use of a finite impulse response (FIR) filter to linearly constrain a desired object while minimizing interfering effects caused by other unknown signal sources. Unlike most spatial-based classification techniques, the proposed CEM takes advantage of spectral characteristics to achieve object detection and classification. A series of experiments is conducted and compared with the commonly used c-means method for performance evaluation. The results show that the CEM method is a promising and effective spectral technique for MR image classification.

Classification of clustered microcalcifications using a Shape Cognitron neural network.

A new shape recognition-based neural network built with universal feature planes, called Shape Cognitron (S-Cognitron) is introduced to classify clustered microcalcifications. The architecture of S-Cognitron consists of two modules and an extra layer, called 3D figure layer lies in between. The first module contains a shape orientation layer, built with 20 cell planes of low level universal shape features to convert first-order shape orientations into numeric values, and a complex layer, to extract second-order shape features. The 3D figure layer is a feature extract-display layer that extracts the shape curvatures of an input pattern and displays them as a 3D figure. It is then followed by a second module made up of a feature formation layer and a probabilistic neural network-based classification layer. The system is evaluated by using Nijmegen mammogram database and experimental results show that sensitivity and specificity can reach 86.1 and 74.1%, respectively.