Rough-Fuzzy Clustering and Unsupervised Feature Selection for Wavelet Based MR Image Segmentation

Image segmentation is an indispensable process in the visualization of human tissues, particularly during clinical analysis of brain magnetic resonance (MR) images. For many human experts, manual segmentation is a difficult and time consuming task, which makes an automated brain MR image segmentation method desirable. In this regard, this paper presents a new segmentation method for brain MR images, integrating judiciously the merits of rough-fuzzy computing and multiresolution image analysis technique. The proposed method assumes that the major brain tissues, namely, gray matter, white matter, and cerebrospinal fluid from the MR images are considered to have different textural properties. The dyadic wavelet analysis is used to extract the scale-space feature vector for each pixel, while the rough-fuzzy clustering is used to address the uncertainty problem of brain MR image segmentation. An unsupervised feature selection method is introduced, based on maximum relevance-maximum significance criterion, to select relevant and significant textural features for segmentation problem, while the mathematical morphology based skull stripping preprocessing step is proposed to remove the non-cerebral tissues like skull. The performance of the proposed method, along with a comparison with related approaches, is demonstrated on a set of synthetic and real brain MR images using standard validity indices.     

Automated Characterization and Parameter-Free Classification of Cell Tracks Based on Local Migration Behavior

Cell migration is the driving force behind the dynamics of many diverse biological processes. Even though microscopy experiments are routinely performed today by which populations of cells are visualized in space and time, valuable information contained in image data is often disregarded because statistical analyses are performed at the level of cell populations rather than at the single-cell level. Image-based systems biology is a modern approach that aims at quantitatively analyzing and modeling biological processes by developing novel strategies and tools for the interpretation of image data. In this study, we take first steps towards a fully automated characterization and parameter-free classification of cell track data that can be generally applied to tracked objects as obtained from image data. The requirements to achieve this aim include: (i) combination of different measures for single cell tracks, such as the confinement ratio and the asphericity of the track volume, and (ii) computation of these measures in a staggered fashion to retrieve local information from all possible combinations of track segments. We demonstrate for a population of synthetic cell tracks as well as for in vitro neutrophil tracks obtained from microscopy experiment that the information contained in the track data is fully exploited in this way and does not require any prior knowledge, which keeps the analysis unbiased and general. The identification of cells that show the same type of migration behavior within the population of all cells is achieved via agglomerative hierarchical clustering of cell tracks in the parameter space of the staggered measures. The recognition of characteristic patterns is highly desired to advance our knowledge about the dynamics of biological processes.     

Automated Classification of Fatty and Normal Liver Ultrasound Images Based on Mutual Information Feature Selection

Background

Fatty Liver Disease (FLD) is one of the most critical diseases that should be detected and cured at the earlier stage in order to decrease the mortality rate. To identify the FLD, ultrasound images have been widely used by the radiologists. However, due to poor quality of ultrasound images, they found difficulties in recognizing FLD. To resolve this problem, many researchers have developed various Computer Aided Diagnosis (CAD) systems for the classification of fatty and normal liver ultrasound images. However, the performance of existing CAD systems is not good in terms of sensitivity while classifying the FLD.

A Multi-Atlas Based Method for Automated Anatomical Rat Brain MRI Segmentation and Extraction of PET Activity

Introduction

Preclinical in vivo imaging requires precise and reproducible delineation of brain structures. Manual segmentation is time consuming and operator dependent. Automated segmentation as usually performed via single atlas registration fails to account for anatomo-physiological variability. We present, evaluate, and make available a multi-atlas approach for automatically segmenting rat brain MRI and extracting PET activies.

Methods

High-resolution 7T 2DT2 MR images of 12 Sprague-Dawley rat brains were manually segmented into 27-VOI label volumes using detailed protocols. Automated methods were developed with 7/12 atlas datasets, i.e. the MRIs and their associated label volumes. MRIs were registered to a common space, where an MRI template and a maximum probability atlas were created. Three automated methods were tested: 1/registering individual MRIs to the template, and using a single atlas (SA), 2/using the maximum probability atlas (MP), and 3/registering the MRIs from the multi-atlas dataset to an individual MRI, propagating the label volumes and fusing them in individual MRI space (propagation & fusion, PF). Evaluation was performed on the five remaining rats which additionally underwent [¹⁸F]FDG PET. Automated and manual segmentations were compared for morphometric performance (assessed by comparing volume bias and Dice overlap index) and functional performance (evaluated by comparing extracted PET measures).

Results

Only the SA method showed volume bias. Dice indices were significantly different between methods (PF>MP>SA). PET regional measures were more accurate with multi-atlas methods than with SA method.

Conclusions

Multi-atlas methods outperform SA for automated anatomical brain segmentation and PET measure’s extraction. They perform comparably to manual segmentation for FDG-PET quantification. Multi-atlas methods are suitable for rapid reproducible VOI analyses.     

An Objective Method to Optimize the MR Sequence Set for Plaque Classification in Carotid Vessel Wall Images Using Automated Image Segmentation

A typical MR imaging protocol to study the status of atherosclerosis in the carotid artery consists of the application of multiple MR sequences. Since scanner time is limited, a balance has to be reached between the duration of the applied MR protocol and the quantity and quality of the resulting images which are needed to assess the disease. In this study an objective method to optimize the MR sequence set for classification of soft plaque in vessel wall images of the carotid artery using automated image segmentation was developed. The automated method employs statistical pattern recognition techniques and was developed based on an extensive set of MR contrast weightings and corresponding manual segmentations of the vessel wall and soft plaque components, which were validated by histological sections. Evaluation of the results from nine contrast weightings showed the tradeoff between scan duration and automated image segmentation performance. For our dataset the best segmentation performance was achieved by selecting five contrast weightings. Similar performance was achieved with a set of three contrast weightings, which resulted in a reduction of scan time by more than 60%. The presented approach can help others to optimize MR imaging protocols by investigating the tradeoff between scan duration and automated image segmentation performance possibly leading to shorter scanning times and better image interpretation. This approach can potentially also be applied to other research fields focusing on different diseases and anatomical regions.     

Multimodality Multiparametric Imaging of Early Tumor Response to a Novel Antiangiogenic Therapy Based on Anticalins

Anticalins are a novel class of targeted protein therapeutics. The PEGylated Anticalin Angiocal (PRS-050-PEG40) is directed against VEGF-A. The purpose of our study was to compare the performance of diffusion weighted imaging (DWI), dynamic contrast enhanced magnetic resonance imaging (DCE)-MRI and positron emission tomography with the tracer [¹⁸F]fluorodeoxyglucose (FDG-PET) for monitoring early response to antiangiogenic therapy with PRS-050-PEG40. 31 mice were implanted subcutaneously with A673 rhabdomyosarcoma xenografts and underwent DWI, DCE-MRI and FDG-PET before and 2 days after i.p. injection of PRS-050-PEG40 (n = 13), Avastin (n = 6) or PBS (n = 12). Tumor size was measured manually with a caliper. Imaging results were correlated with histopathology. In the results, the tumor size was not significantly different in the treatment groups when compared to the control group on day 2 after therapy onset (P = 0.09). In contrast the imaging modalities DWI, DCE-MRI and FDG-PET showed significant differences between the therapeutic compared to the control group as early as 2 days after therapy onset (P<0.001). There was a strong correlation of the early changes in DWI, DCE-MRI and FDG-PET at day 2 after therapy onset and the change in tumor size at the end of therapy (r = −0.58, 0.71 and 0.67 respectively). The imaging results were confirmed by histopathology, showing early necrosis and necroptosis in the tumors. Thus multimodality multiparametric imaging was able to predict therapeutic success of PRS-050-PEG40 and Avastin as early as 2 days after onset of therapy and thus promising for monitoring early response of antiangiogenic therapy.     

基于Snake模型的医学图像分割

基于Snake模型的图像分割技术是近年来图像处理领域的研究热点之一。Snake模型承载上层先验知识并融合了图像的底层特征,针对医学图像的特殊性,能有效地应用于医学图像的分割中。本文对各种基于Snake模型的改进算法和进化模型进行了研究,并重点梳理了最新的研究成果,以利于把握基于Snake模型的医学图像分割方法的脉络和发展方向。     

Sequential Classification Based on Regression

Let us consider m general populations π₁, …,π_m. Each object belonging to these populations is represented by (p ± 1) characteristics x₁, x₂,…,x_p,y. A certain object, which is an element of one of the m general populations π₁,…,π_m has to be classified into the correct population. It will be assumed that knowledge of the value of the characteristic y would permit its correct classification, but that the observation of this characteristic is expensive, difficult or dangerous, as e.g. in medical applications. y is correlated with a set of p characteristics x₁,x₂,…,x_p, which are observed sequentially. The classification procedure is based on the division of the space of the observed value of characteristics x₁,x₂,…,x_p into nonintersecting areas determined so as to minimize the value of BAYES' risk given by equation (3).     

Translational Validation of Personalized Treatment Strategy Based on Genetic Characteristics of Glioblastoma

Glioblastoma (GBM) heterogeneity in the genomic and phenotypic properties has potentiated personalized approach against specific therapeutic targets of each GBM patient. The Cancer Genome Atlas (TCGA) Research Network has been established the comprehensive genomic abnormalities of GBM, which sub-classified GBMs into 4 different molecular subtypes. The molecular subtypes could be utilized to develop personalized treatment strategy for each subtype. We applied a classifying method, NTP (Nearest Template Prediction) method to determine molecular subtype of each GBM patient and corresponding orthotopic xenograft animal model. The models were derived from GBM cells dissociated from patient''s surgical sample. Specific drug candidates for each subtype were selected using an integrated pharmacological network database (PharmDB), which link drugs with subtype specific genes. Treatment effects of the drug candidates were determined by in vitro limiting dilution assay using patient-derived GBM cells primarily cultured from orthotopic xenograft tumors. The consistent identification of molecular subtype by the NTP method was validated using TCGA database. When subtypes were determined by the NTP method, orthotopic xenograft animal models faithfully maintained the molecular subtypes of parental tumors. Subtype specific drugs not only showed significant inhibition effects on the in vitro clonogenicity of patient-derived GBM cells but also synergistically reversed temozolomide resistance of MGMT-unmethylated patient-derived GBM cells. However, inhibitory effects on the clonogenicity were not totally subtype-specific. Personalized treatment approach based on genetic characteristics of each GBM could make better treatment outcomes of GBMs, although more sophisticated classifying techniques and subtype specific drugs need to be further elucidated.     

Online Adaptation of a c-VEP Brain-Computer Interface(BCI) Based on Error-Related Potentials and Unsupervised Learning

The goal of a Brain-Computer Interface (BCI) is to control a computer by pure brain activity. Recently, BCIs based on code-modulated visual evoked potentials (c-VEPs) have shown great potential to establish high-performance communication. In this paper we present a c-VEP BCI that uses online adaptation of the classifier to reduce calibration time and increase performance. We compare two different approaches for online adaptation of the system: an unsupervised method and a method that uses the detection of error-related potentials. Both approaches were tested in an online study, in which an average accuracy of 96% was achieved with adaptation based on error-related potentials. This accuracy corresponds to an average information transfer rate of 144 bit/min, which is the highest bitrate reported so far for a non-invasive BCI. In a free-spelling mode, the subjects were able to write with an average of 21.3 error-free letters per minute, which shows the feasibility of the BCI system in a normal-use scenario. In addition we show that a calibration of the BCI system solely based on the detection of error-related potentials is possible, without knowing the true class labels.     

Automated Protein Subcellular Localization Based on Local Invariant Features

To understand the function of the encoded proteins, we need to be able to know the subcellular location of a protein. The most common method used for determining subcellular location is fluorescence microscopy which allows subcellular localizations to be imaged in high throughput. Image feature calculation has proven invaluable in the automated analysis of cellular images. This article proposes a novel method named LDPs for feature extraction based on invariant of translation and rotation from given images, the nature which is to count the local difference features of images, and the difference features are given by calculating the D-value between the gray value of the central pixel c and the gray values of eight pixels in the neighborhood. The novel method is tested on two image sets, the first set is which fluorescently tagged protein was endogenously expressed in 10 sebcellular locations, and the second set is which protein was transfected in 11 locations. A SVM was trained and tested for each image set and classification accuracies of 96.7 and 92.3 % were obtained on the endogenous and transfected sets respectively.     

Double-Exposure Optical Sectioning Structured Illumination Microscopy Based on Hilbert Transform Reconstruction

Structured illumination microscopy (SIM) with axially optical sectioning capability has found widespread applications in three-dimensional live cell imaging in recent years, since it combines high sensitivity, short image acquisition time, and high spatial resolution. To obtain one sectioned slice, three raw images with a fixed phase-shift, normally 2π/3, are generally required. In this paper, we report a data processing algorithm based on the one-dimensional Hilbert transform, which needs only two raw images with arbitrary phase-shift for each single slice. The proposed algorithm is different from the previous two-dimensional Hilbert spiral transform algorithm in theory. The presented algorithm has the advantages of simpler data processing procedure, faster computation speed and better reconstructed image quality. The validity of the scheme is verified by imaging biological samples in our developed DMD-based LED-illumination SIM system.     

Improved Estimation of Cardiac Function Parameters Using a Combination of Independent Automated Segmentation Results in Cardiovascular Magnetic Resonance Imaging

This work aimed at combining different segmentation approaches to produce a robust and accurate segmentation result. Three to five segmentation results of the left ventricle were combined using the STAPLE algorithm and the reliability of the resulting segmentation was evaluated in comparison with the result of each individual segmentation method. This comparison was performed using a supervised approach based on a reference method. Then, we used an unsupervised statistical evaluation, the extended Regression Without Truth (eRWT) that ranks different methods according to their accuracy in estimating a specific biomarker in a population. The segmentation accuracy was evaluated by estimating six cardiac function parameters resulting from the left ventricle contour delineation using a public cardiac cine MRI database. Eight different segmentation methods, including three expert delineations and five automated methods, were considered, and sixteen combinations of the automated methods using STAPLE were investigated. The supervised and unsupervised evaluations demonstrated that in most cases, STAPLE results provided better estimates than individual automated segmentation methods. Overall, combining different automated segmentation methods improved the reliability of the segmentation result compared to that obtained using an individual method and could achieve the accuracy of an expert.     

Motor Oil Classification Based on Time-Resolved Fluorescence

A time-resolved fluorescence (TRF) technique is presented for classifying motor oils. The system is constructed with a third harmonic Nd:YAG laser, a spectrometer, and an intensified charge coupled device (ICCD) camera. Steady-state and time-resolved fluorescence (TRF) measurements are reported for several motor oils. It is found that steady-state fluorescence is insufficient to distinguish the motor oil samples. Then contour diagrams of TRF intensities (CDTRFIs) are acquired to serve as unique fingerprints to identify motor oils by using the distinct TRF of motor oils. CDTRFIs are preferable to steady-state fluorescence spectra for classifying different motor oils, making CDTRFIs a particularly choice for the development of fluorescence-based methods for the discrimination and characterization of motor oils. The two-dimensional fluorescence contour diagrams contain more information, not only the changing shapes of the LIF spectra but also the relative intensity. The results indicate that motor oils can be differentiated based on the new proposed method, which provides reliable methods for analyzing and classifying motor oils.     

Palm-Vein Classification Based on Principal Orientation Features

Personal recognition using palm–vein patterns has emerged as a promising alternative for human recognition because of its uniqueness, stability, live body identification, flexibility, and difficulty to cheat. With the expanding application of palm–vein pattern recognition, the corresponding growth of the database has resulted in a long response time. To shorten the response time of identification, this paper proposes a simple and useful classification for palm–vein identification based on principal direction features. In the registration process, the Gaussian-Radon transform is adopted to extract the orientation matrix and then compute the principal direction of a palm–vein image based on the orientation matrix. The database can be classified into six bins based on the value of the principal direction. In the identification process, the principal direction of the test sample is first extracted to ascertain the corresponding bin. One-by-one matching with the training samples is then performed in the bin. To improve recognition efficiency while maintaining better recognition accuracy, two neighborhood bins of the corresponding bin are continuously searched to identify the input palm–vein image. Evaluation experiments are conducted on three different databases, namely, PolyU, CASIA, and the database of this study. Experimental results show that the searching range of one test sample in PolyU, CASIA and our database by the proposed method for palm–vein identification can be reduced to 14.29%, 14.50%, and 14.28%, with retrieval accuracy of 96.67%, 96.00%, and 97.71%, respectively. With 10,000 training samples in the database, the execution time of the identification process by the traditional method is 18.56 s, while that by the proposed approach is 3.16 s. The experimental results confirm that the proposed approach is more efficient than the traditional method, especially for a large database.