Hair and bare skin discrimination for laser-assisted hair removal systems.

Laser-assisted hair removal devices aim to remove body hair permanently. In most cases, these devices irradiate the whole area of the skin with a homogenous power density. Thus, a significant portion of the skin, where hair is not present, is burnt unnecessarily causing health risks. Therefore, methods that can distinguish hair regions automatically would be very helpful avoiding these unnecessary applications of laser. This study proposes a new system of algorithms to detect hair regions with the help of a digital camera. Unlike previous limited number of studies, our methods are very fast allowing for real-time application. Proposed methods are based on certain features derived from histograms of hair and skin regions. We compare our algorithm with competing methods in terms of localization performance and computation time and show that a much faster real-time accurate localization of hair regions is possible with the proposed method. Our results show that the algorithm we have developed is extremely fast (around 45 milliseconds) allowing for real-time application with high accuracy hair localization ( 96.48 %).

Fully automated blink detection for uveal melanoma radiotherapy.

Uveal melanoma is a widespread neoplasm in the eye that cause vision defects. The disease may result in loss of eye, even loss of life with a risk of metastasis to other tissues. There have been proposed various treatment methods to overcome the disease including enucleation and radiotherapy. Radiotherapy treatment is certainly preferable since it avoids the loss of eye. The treatment is applied with radioactive beam directed to the tumor, and thus the eye must be kept steady usually with local anesthesia. However, anesthesia has many health risks; therefore it is preferable to perform this operation without it. This paper presents an approach where a camera will be recording the eye and the radiotherapy will immediately be halted when the eye is closed to avoid unnecessary beams directed to healthy regions of the eye. System is based on the classification of the acquired image frames as "eye open" and "eye closed". The speed of the algorithm is aimed as close as possible to real time so that it can respond instantly to any change in present eye state. A total of 480 video frames, 120 frames for each open and closed eyes in bright and dark lighting conditions separately, were used for experimentation. Test results indicate that developed features are successful in terms of being very fast providing no classification error.

Cross-device automated prostate cancer localization with multiparametric MRI.

Prostate cancer localization using supervised classification techniques has aroused considerable interest in medical imaging community in recent years. However, it is crucial to have an accurate training data set for supervised classification techniques. Since different devices with, e.g., different protocols and/or field strengths cause different intensity profiles, each device/protocol must have an accompanying training data set, which is very costly to obtain. It is highly desirable to adapt the existing classifier(s) trained for one device/protocol to help classify data coming from another device/protocol. In this paper, we propose a novel method that has the ability to design classifiers obtained from one imaging protocol and/or MRI device to be used on a data set from another protocol and/or imaging device. As an example problem, we consider prostate cancer localization with multiparametric MRI. We show that simple normalization techniques such as z-score are not sufficient for cross-device automated cancer localization. On the other hand, the method we have originally developed based on relative intensity allows us to successfully use a classifier obtained from one device to be applied on a test patient imaged with another device. Proposed method also allows us to employ T2-weighted MR images directly instead of an additional step to normalize T2-weighted images usually performed in an ad hoc manner when T2 maps are not available. To demonstrate the effectiveness of the proposed method, we use a multiparametric MRI data set acquired from 18 biopsy-confirmed cancer patients with two separate scanners: 1) 1.5-T (Excite HD) GE and 2) 1.5-T (Achieva) Philips Healthcare scanners. A comprehensive visual, quantitative, and statistical analysis of the results show that methods we have developed allow us to: 1) perform cross-device automated classification and 2) use T2-weighted images without an ad hoc subject-specific normalization.

Prostate cancer localization using multiparametric MRI based on semi-supervised techniques with automated seed initialization.

In this paper, we propose a novel and efficient semisupervised technique for automated prostate cancer localization using multiparametric magnetic resonance imaging (MRI). This method can be used in guiding biopsy, surgery, and therapy. We systematically present a new segmentation technique by developing a multiparametric graph based random walker (RW) algorithm with automated seed initialization to perform prostate cancer segmentation using multiparametric MRI. RW algorithm has proved to be accurate and fast in segmentation applications; however it requires a set of (user provided) seed points in order to perform segmentation. In this study, we first developed a novel RW method, which can be used with multiparametric MR images and then devised alternative methods that can determine seed points in an automated manner using discriminative classifiers such as support vector machines (SVM). Proposed RW method with automated seed initialization is able to produce improved segmentation results by assigning more weights to the images with more discriminative power.We applied the proposed method to a multiparametric dataset obtained from biopsy confirmed prostate cancer patients. Proposed method produces a sensitivity/ specificity rate of 0.76 and 0.86, respectively. Both visual, quantitative as well as statistical results are presented to show the significant performance improvements. Fisher sign test is used to demonstrate the statistical significance of our results by achieving p-values less than 0.05. This method outperforms available RW and SVM based methods by achieving a high specificity rate while not reducing sensitivity.

A new ROC analysis method considering the correlation between neighboring pixels.

In this paper, we introduce a novel receiver operating characteristic (ROC) analysis method that considers spatial correlation between pixels to evaluate classification algorithms. ROC analysis is one of the most important tools in the evaluation of medical images and computer aided diagnosis (CAD) systems. It provides a comprehensive description of the detection accuracy of the test image. To evaluate the localization performance, operating points of ROC curves are obtained based on the classification results of individual pixels. To this date, the confidence level or intensity value of each pixel is assumed to be independent within the image. However, this assumption is not satisfied in real problems. In this paper, a new ROC analysis algorithm that considers the correlation between neighboring pixels is proposed. Our results show that the new ROC curves provide a more accurate evaluation of the test image.

Cross-device automated prostate cancer localization with multiparametric MRI.

Automated cancer localization with supervised techniques plays an important role in guiding biopsy, surgery and treatment. It is crucial to have an accurate training dataset for supervised techniques. Since different devices with e.g. different protocols and/or field strengths cause different intensity profiles, each device/protocol must have an accompanying training dataset which is very costly to obtain. In this paper, we propose a novel method that has the ability to design classifiers obtained from one imaging protocol and/or MRI device to be used on a dataset from another protocol and/or imaging device. As an example problem we consider prostate cancer localization with multiparametric MRI. We show that simple normalization techniques such as z-score are not sufficient to allow for cross-device automated cancer localization. On the other hand, the methods we have originally developed based on relative intensity allows us to successfully use a classifier obtained from one device to be applied on a test patient imaged with another device.

Automated prostate cancer localization without the need for peripheral zone extraction using multiparametric MRI.

PURPOSE: Multiparametric magnetic resonance imaging (MRI) has been shown to have higher localization accuracy than transrectal ultrasound (TRUS) for prostate cancer. Therefore, automated cancer segmentation using multiparametric MRI is receiving a growing interest, since MRI can provide both morphological and functional images for tissue of interest. However, all automated methods to this date are applicable to a single zone of the prostate, and the peripheral zone (PZ) of the prostate needs to be extracted manually, which is a tedious and time-consuming job. In this paper, our goal is to remove the need of PZ extraction by incorporating the spatial and geometric information of prostate tumors with multiparametric MRI derived from T2-weighted MRI, diffusion-weighted imaging (DWI) and dynamic contrast enhanced MRI (DCE-MRI). METHODS: In order to remove the need of PZ extraction, the authors propose a new method to incorporate the spatial information of the cancer. This is done by introducing a new feature called location map. This new feature is constructed by applying a nonlinear transformation to the spatial position coordinates of each pixel, so that the location map implicitly represents the geometric position of each pixel with respect to the prostate region. Then, this new feature is combined with multiparametric MR images to perform tumor localization. The proposed algorithm is applied to multiparametric prostate MRI data obtained from 20 patients with biopsy-confirmed prostate cancer. RESULTS: The proposed method which does not need the masks of PZ was found to have prostate cancer detection specificity of 0.84, sensitivity of 0.80 and dice coefficient value of 0.42. CONCLUSIONS: The authors have found that fusing the spatial information allows us to obtain tumor outline without the need of PZ extraction with a considerable success (better or similar performance to methods that require manual PZ extraction). Our experimental results quantitatively demonstrate the effectiveness of the proposed method, depicting that the proposed method has a slightly better or similar localization performance compared to methods which require the masks of PZ.

Propagation of Blood Function Errors to the Estimates of Kinetic Parameters with Dynamic PET.

Dynamic PET, in contrast to static PET, can identify temporal variations in the radiotracer concentration. Mathematical modeling of the tissue of interest in dynamic PET can be simplified using compartment models as a linear system where the time activity curve of a specific tissue is the convolution of the tracer concentration in the plasma and the impulse response of the tissue containing kinetic parameters. Since the arterial sampling of blood to acquire the value of tracer concentration is invasive, blind methods to estimate both blood input function and kinetic parameters have recently drawn attention. Several methods have been developed, but the effect of accuracy of the estimated blood function on the estimation of the kinetic parameters is not studied. In this paper, we present a method to compute the error in the kinetic parameter estimates caused by the error in the blood input function. Computer simulations show that analytical expressions we derive are sufficiently close to results obtained from numerical methods. Our findings are important to observe the effect of the blood function on kinetic parameter estimation, but also useful to evaluate various blind methods and observe the dependence of kinetic parameter estimates to certain parts of the blood function.

Prostate cancer localization with multispectral MRI using cost-sensitive support vector machines and conditional random fields.

Prostate cancer is a leading cause of cancer death for men in the United States. Fortunately, the survival rate for early diagnosed patients is relatively high. Therefore, in vivo imaging plays an important role for the detection and treatment of the disease. Accurate prostate cancer localization with noninvasive imaging can be used to guide biopsy, radiotherapy, and surgery as well as to monitor disease progression. Magnetic resonance imaging (MRI) performed with an endorectal coil provides higher prostate cancer localization accuracy, when compared to transrectal ultrasound (TRUS). However, in general, a single type of MRI is not sufficient for reliable tumor localization. As an alternative, multispectral MRI, i.e., the use of multiple MRI-derived datasets, has emerged as a promising noninvasive imaging technique for the localization of prostate cancer; however almost all studies are with human readers. There is a significant inter and intraobserver variability for human readers, and it is substantially difficult for humans to analyze the large dataset of multispectral MRI. To solve these problems, this study presents an automated localization method using cost-sensitive support vector machines (SVMs) and shows that this method results in improved localization accuracy than classical SVM. Additionally, we develop a new segmentation method by combining conditional random fields (CRF) with a cost-sensitive framework and show that our method further improves cost-sensitive SVM results by incorporating spatial information. We test SVM, cost-sensitive SVM, and the proposed cost-sensitive CRF on multispectral MRI datasets acquired from 21 biopsy-confirmed cancer patients. Our results show that multispectral MRI helps to increase the accuracy of prostate cancer localization when compared to single MR images; and that using advanced methods such as cost-sensitive SVM as well as the proposed cost-sensitive CRF can boost the performance significantly when compared to SVM.

Regularized estimation of retinal vascular oxygen tension from phosphorescence images.

The level of retinal oxygenation is potentially an important cue to the onset or presence of some common retinal diseases. An improved method for assessing oxygen tension in retinal blood vessels from phosphorescence lifetime imaging data is reported in this paper. The optimum estimate for phosphorescence lifetime and oxygen tension is obtained by regularizing the least-squares (LS) method. The estimation method is implemented with an iterative algorithm to minimize a regularized LS cost function. The effectiveness of the proposed method is demonstrated by applying it to simulated data as well as image data acquired from rat retinas. The method is shown to yield estimates that are robust to noise and whose variance is lower than that obtained with the classical LS method.

Prostate cancer segmentation with simultaneous estimation of Markov random field parameters and class.

Prostate cancer is one of the leading causes of death from cancer among men in the United States. Currently, high-resolution magnetic resonance imaging (MRI) has been shown to have higher accuracy than trans-rectal ultrasound (TRUS) when used to ascertain the presence of prostate cancer. As MRI can provide both morphological and functional images for a tissue of interest, some researchers are exploring the uses of multispectral MRI to guide prostate biopsies and radiation therapy. However, success with prostate cancer localization based on current imaging methods has been limited due to overlap in feature space of benign and malignant tissues using any one MRI method and the interobserver variability. In this paper, we present a new unsupervised segmentation method for prostate cancer detection, using fuzzy Markov random fields (fuzzy MRFs) for the segmentation of multispectral MR prostate images. Typically, both hard and fuzzy MRF models have two groups of parameters to be estimated: the MRF parameters and class parameters for each pixel in the image. To date, these two parameters have been treated separately, and estimated in an alternating fashion. In this paper, we develop a new method to estimate the parameters defining the Markovian distribution of the measured data, while performing the data clustering simultaneously. We perform computer simulations on synthetic test images and multispectral MR prostate datasets to demonstrate the efficacy and efficiency of the proposed method and also provide a comparison with some of the commonly used methods.

Parametric surface-source modeling and estimation with electroencephalography.

Electroencephalography (EEG) is an important tool for studying the brain functions and is becoming popular in clinical practice. In this paper, we develop four parametric EEG models to estimate current sources that are spatially distributed on a surface. Our models approximate the source shape and extent explicitly and can be applied to localize extended sources which are often encountered, e.g., in epilepsy diagnosis. We assume a realistic head model and solve the EEG forward problem using the boundary element method. We present the source models with increasing degrees of freedom, provide the forward solutions, and derive the maximum-likelihood estimates as well as Cramér-Rao bounds of the unknown source parameters. In order to evaluate the applicability of the proposed models, we first compare their estimation performances with the dipole model's using several known source distributions. We then discuss the conditions under which we can distinguish between the proposed extended sources and the focal dipole using the generalized likelihood ratio test. We also apply our models to the electric measurements obtained from a phantom body in which an extended electric source is imbedded. We observe that the proposed model can capture the source extent information satisfactorily and the localization accuracy is better than the dipole model.

Estimating parametric line-source models with electroencephalography.

We develop three parametric models for electroencephalography (EEG) to estimate current sources that are spatially distributed on a line. We assume a realistic head model and solve the EEG forward problem using the boundary element method (BEM). We present the models with increasing degrees of freedom, provide the forward solutions, and derive the maximum-likelihood estimates as well as Cramér-Rao bounds of the unknown source parameters. A series of experiments are conducted to evaluate the applicability of the proposed models. We use numerical examples to demonstrate the usefulness of our line-source models in estimating extended sources. We also apply our models to the real EEG data of N20 response that is known to have an extended source. We observe that the line-source models explain the N20 measurements better than the dipole model.

Surface-source modeling and estimation using biomagnetic measurements.

We propose a number of electric source models that are spatially distributed on an unknown surface for biomagnetism. These can be useful to model, e.g., patches of electrical activity on the cortex. We use a realistic head (or another organ) model and discuss the special case of a spherical head model with radial sensors resulting in more efficient computations of the estimates for magnetoencephalography. We derive forward solutions, maximum likelihood (ML) estimates, and Cramér-Rao bound (CRB) expressions for the unknown source parameters. A model selection method is applied to decide on the most appropriate model. We also present numerical examples to compare the performances and computational costs of the different models and illustrate when it is possible to distinguish between surface and focal sources or line sources. Finally, we apply our methods to real biomagnetic data of phantom human torso and demonstrate the applicability of them.

Line-source modeling and estimation with magnetoencephalography.

We propose a number of source models that are spatially distributed on a line for magnetoencephalography (MEG) using both a spherical head with radial sensors for more efficient computation and a realistic head model for more accurate results. We develop these models with increasing degrees of freedom, derive forward solutions, maximum-likelihood (ML) estimates, and Cramér-Rao bound (CRB) expressions for the unknown source parameters. A model selection method is applied to select the most appropriate model. We also present numerical examples to compare the performances and computational costs of the different models, to determine the regions where better estimates are possible and when it is possible to distinguish between line and focal sources. We demonstrate the usefulness of the proposed line-source models over the previously available focal source model in certain distributed source cases. Finally, we apply our methods to real MEG data, the N2O response after electric stimulation of the median nerve known to be an extended source.

Distinguishing between moving and stationary sources using EEG/MEG measurements with an application to epilepsy.

Performances of electroencephalography (EEG) and magnetoencephalography (MEG) source estimation methods depend on the validity of the assumed model. In many cases, the model structure is related to physical information. We discuss a number of statistical selection methods to distinguish between two possible models using least-squares estimation and assuming a spherical head model. The first model has a single moving source whereas the second has two stationary sources; these may result in similar EEG/MEG measurements. The need to decide between such models occurs for example in Jacksonian seizures (e.g., epilepsy) or in intralobular activities, where a model with either two stationary dipole sources or a single moving dipole source may be possible. We also show that all of the selection methods discussed choose the correct model with probability one when the number of trials goes to infinity. Finally we present numerical examples and compare the performances of the methods by varying parameters such as the signal-to-noise ratio, source depth, and separation of sources, and also apply the methods to real MEG data for epilepsy.

Optimization of orders in multichannel fractional Fourier-domain filtering circuits and its application to the synthesis of mutual-intensity distributions.

Owing to the nonlinear nature of the problem, the transform orders in fractional Fourier-domain filtering configurations have usually not been optimized but chosen uniformly. We discuss the optimization of these orders for multi-channel-filtering configurations by first finding the optimal filter coefficients for a larger number of uniformly chosen orders, and then maintaining the most important ones. The method is illustrated with the problem of synthesizing desired mutual-intensity distributions. The method we propose allows those fractional Fourier domains, which add little benefit to the filtering process but increase the overall cost, to be pruned, so that comparable performance can be attained with less cost, or higher performance can be obtained with the same cost. The method we propose is more likely to be useful when confronted with low-cost rather than high-performance applications, because larger improvements are obtained when the use of a smaller number of filters is desired.