Search

Objective: We propose a MATLAB-based tool to convert electrocardiography (ECG) waveforms from paper-based ECG records into digitized ECG signals that is vendor-agnostic. The tool is packaged as an open source standalone graphical user interface (GUI) based application. Methods and procedures: To reach this objective we: (1) preprocess the ECG records, which includes skew correction, background grid removal and linear filtering; (2) segment ECG signals using Connected Components Analysis (CCA); (3) implement Optical Character Recognition (OCR) for removal of overlapping ECG lead characters and for interfacing of patients’ demographic information with their research records or their electronic medical record (EMR). The ECG digitization results are validated through a reader study where clinically salient features, such as intervals of QRST complex, between the paper ECG records and the digitized ECG records are compared. Results: Comparison of clinically important features between the paper-based ECG records and the digitized ECG signals, reveals intra- and inter-observer correlations of 0.86–0.99 and 0.79–0.94, respectively. The kappa statistic was found to average at 0.86 and 0.72 for intra- and inter-observer correlations, respectively. Conclusion: The clinically salient features of the ECG waveforms such as the intervals of QRST complex, are preserved during the digitization procedure. Clinical and Healthcare Impact: This open-source digitization tool can be used as a research resource to digitize paper ECG records thereby enabling development of new prediction algorithms to risk stratify individuals with cardiovascular disease, and/or allow for development of ECG-based cardiovascular diagnoses relying upon automated digital algorithms.

This paper reviews the recent literature on technologies and methodologies for quantitative human gait analysis in the context of neurodegnerative diseases. The use of technological instruments can be of great support in both clinical diagnosis and severity assessment of these pathologies. In this paper, sensors, features and processing methodologies have been reviewed in order to provide a highly consistent work that explores the issues related to gait analysis. First, the phases of the human gait cycle are briefly explained, along with some non-normal gait patterns (gait abnormalities) typical of some neurodegenerative diseases. The work continues with a survey on the publicly available datasets principally used for comparing results. Then the paper reports the most common processing techniques for both feature selection and extraction and for classification and clustering. Finally, a conclusive discussion on current open problems and future directions is outlined.

Rehabilitation level evaluation is an important part of the automatic rehabilitation training system. As a general rule, this process is manually performed by rehabilitation doctors using chart-based ordinal scales which can be both subjective and inefficient. In this paper, a novel approach based on ensemble learning is proposed which automatically evaluates stroke patients' rehabilitation level using multi-channel sEMG signals to this problem. The correlation between rehabilitation levels and rehabilitation training actions is investigated and actions suitable for rehabilitation assessment are selected. Then, features are extracted from the selected actions. Finally, the features are used to train the stacking classification model. Experiments using sEMG data collected from 24 stroke patients have been carried out to examine the validity and feasibility of the proposed method. The experiment results show that the algorithm proposed in this paper can improve the classification accuracy of 6 Brunnstrom stages to 94.36%, which can promote the application of home-based rehabilitation training in practice.

In clinical conversational applications, extracted entities tend to capture the main subject of a patient's complaint, namely symptoms or diseases. However, they mostly fail to recognize the characterizations of a complaint such as the time, the onset, and the severity. For example, if the input is "I have a headache and it is extreme", state-of-the-art models only recognize the main symptom entity - headache, but ignore the severity factor of "extreme", that characterizes headache. In this paper, we design a two-stage approach to detect the characterizations of entities like symptoms presented by general users in contexts where they would describe their symptoms to a clinician. We use Word2Vec and BERT to encode clinical text given by the patients. We transform the output and re-frame the task as multi-label classification problem. Finally, we combine the processed encodings with the Linear Discriminant Analysis (LDA) algorithm to classify the characterizations of the main entity. Experimental results demonstrate that our method achieves 40-50% improvement on the accuracy over the state-of-the-art models., Comment: This paper has been communicated to 42nd IEEE Annual Conference of Engineering in Medicine and Biology Society, Montreal, Canada. It is 5 pages long. It has 2 figures

Despite all the recent developments of using the surface electromyography (sEMG) as a control signal, reliable classifications still remain an arduous task due to overlapping classes and classification ripples. In this paper, we present a straightforward approach to avoid classification ripple based on smoothing the arg max value of an Extreme Learning Machine (ELM) classifier. We compare the baseline accuracy of the classifier with an arg max filtered by a traditional Exponential Smoothing Filter (ESF) and our adaptation of Antonyan Vardan Transform (AVT). The classifiers were evaluated using sEMG data acquired through 12 channels from four subjects performing 17 different movements of forearm and fingers with three repetitions each. In the best scenario, our methods reached results higher than 96% and 82% of overall and weighted accuracy, respectively. Those results match or outperform similar papers of the literature using a simpler model, which may help the application of the techniques on embedded platforms and make the practical use of such devices more feasible.

Clinical assessment of abnormal neuromechanics is typically performed by manipulation of the affected limbs; a process with low inter- and intra-rater reliability. This paper aims at formalizing a framework that closes the loop between a clinician's expertise and computational algorithms, to enhance the clinician's diagnostic capabilities during physical manipulation. The framework's premise is that the dynamics that can be measured by manipulation of a limb are distinct between movement disorders. An a priori database contains measurements encoded in a space called the information map. Based on this map, a computational algorithm identifies which probing motions are more likely to yield distinguishing information about a patient's movement disorder. The clinician executes this movement and the resulting dynamics, combined with clinician input, is used by the algorithm to estimate which of the movement disorders in the database are most probable. This is recursively repeated until a diagnosis can be confidently made. The main contributions of this paper are the formalization of the framework and the addition of the information map to select informative movements. The establishment of the framework provides a foundation for a standardized assessment of movement disorders and future work will aim at testing the framework's efficacy.

This paper presents a voice control interface prototype for assistive robots aiming to help people living with upper limb disabilities to perform daily activities autonomously. Assistive robotic devices can be used to help people with upper-body disabilities gain more autonomy in their daily life. However, it is very difficult or even impossible for certain users to control the robot with conventional control systems (e.g. joystick, sip-and-puff). This paper presents the design and preliminary evaluation of a voice command system prototype for the control of assistive robotic arms' movements. This work aims at making the control of assistive robots more intuitive and fluid, and to perform various tasks in less time and with a lesser effort. The prototype of the voice command interface developed is first presented, followed by two experiments with five able-bodied subjects in order to assess the system's performance and guide future development.

This paper presents the development of a new active assistive eating device, which aims to stabilize the movements of people living with movement disorders, such as spasticity and ataxia. Many people living with upper-body incapacities are unable to eat on their own, due to movement disorders (ex. tremors, spastic motions, lack of muscular tone), resulting from various ailments like Cerebral palsy, Parkinson's disease, Dystonia, Multiple sclerosis, strokes, and Muscular dystrophy). Our past work focused on the development of a purely mechanical device, which involved damping of the system via passive mechanical dampers. This paper extends said work by using active stabilization of user movements. The active assistance enables the design of intelligent algorithms that can assist human movements more efficiently. This active version has the benefits of being easily adjustable; the level of damping can be adjusted in real-time, depending on the user movement; different control modes are offered, and the guiding of user movements is also allowed. Firstly, the mechanical design of the device is presented, followed by the damping arrangement, the electronic design, the control algorithms and finally, the preliminary experiments are mentioned.

Alzheimer's disease (AD) accounts for 60%-70% of all dementia cases, and clinical diagnosis at its early stage is extremely difficult. As several new drugs aiming to modify disease progression or alleviate symptoms are being developed, to assess their efficacy, novel robust biomarkers of brain function are urgently required. This paper aims to explore a routine to gain such biomarkers using the quantitative analysis of electroencephalography (QEEG). This paper proposes a supervised classification framework that uses EEG signals to classify healthy controls (HC) and AD participants. The framework consists of data augmentation, feature extraction, K-nearest neighbor (KNN) classification, quantitative evaluation, and topographic visualization. Considering the human brain either as a stationary or a dynamical system, both the frequency-based and time-frequency-based features were tested in 40 participants. The results show that: 1) the proposed method can achieve up to a 99% classification accuracy on short (4s) eyes open EEG epochs, with the KNN algorithm that has best performance when compared with alternative machine learning approaches 2) the features extracted using the wavelet transform produced better classification performance in comparison to the features based on FFT and 3) in the spatial domain, the temporal and parietal areas offer the best distinction between healthy controls and AD. The proposed framework can effectively classify HC and AD participants with high accuracy, meanwhile offering identification and the localization of significant QEEG features. These important findings and the proposed classification framework could be used for the development of a biomarker for the diagnosis and monitoring of disease progression in AD. 27 5 826 835

Statistical analysis using large-scale brain magnetic resonance (MR) image databases has examined that brain tissues have age-related morphological changes. The age of a subject can be estimated from the brain MR image by evaluating morphological changes with healthy aging. This paper proposes an age estimation method using local features of T1-weighted MR images. The brain local features are defined by volumes of brain tissues parcellated into 1,024 local regions defined by the automated anatomical labeling atlas. This paper also proposes the effective local feature selection method to improve the accuracy of age estimation. We evaluate the accuracy of the proposed method using 1,099 T1-weighted images from a Japanese MR image database. We also analyze effectiveness of each local region for age estimation and discuss its medical implication.

Support Vector Machine (SVM) is a common classifier used for efficient classification with high accuracy. SVM shows high accuracy for classifying melanoma (skin cancer) clinical images within computer-aided diagnosis systems used by skin cancer specialists to detect melanoma early and save lives. We aim to develop a medical low-cost handheld device that runs a real-time embedded SVM- based diagnosis system for use in primary care for early detection of melanoma. In this paper, an optimized SVM classifier is implemented onto a recent FPGA platform using the latest design methodology to be embedded into the proposed device for realizing online efficient melanoma detection on a single system on chip/device. The hardware implementation results demonstrate a high classification accuracy of 97.9% and a significant acceleration factor of 26 from equivalent software implementation on an embedded processor, with 34% of resources utilization and 2 watts for power consumption. Consequently, the implemented system meets crucial embedded systems constraints of high performance and low cost, resources utilization and power consumption, while achieving high classification accuracy., Conference paper, 5 pages, 4 figures, 1 tables

This paper considers the auditory attention detection (AAD) paradigm, where the goal is to determine which of two simultaneous speakers a person is attending to. The paradigm relies on recordings of the listener's brain activity, e.g., from electroencephalography (EEG). To perform AAD, decoded EEG signals are typically correlated with the temporal envelopes of the speech signals of the separate speakers. In this paper, we study how the inclusion of various degrees of auditory modelling in this speech envelope extraction process affects the AAD performance, where the best performance is found for an auditory-inspired linear filter bank followed by power law compression. These two modelling stages are computationally cheap, which is important for implementation in wearable devices, such as future neuro-steered auditory prostheses. We also introduce a more natural way to combine recordings (over trials and subjects) to train the decoder, which reduces the dependence of the algorithm on regularization parameters. Finally, we investigate the simultaneous design of the EEG decoder and the audio subband envelope recombination weights vector using either a norm-constrained least squares or a canonical correlation analysis, but conclude that this increases computational complexity without improving AAD performance. ispartof: IEEE Transactions on Neural Systems and Rehabilitation Engineering vol:25 issue:2 pages:1-11 ispartof: location:United States status: published

A novel sensor placement method for better removal of motion artifacts (MA) from photoplethysmography (PPG) signal using Fast Independent Component Analysis (ICA) is proposed in this paper. The method enhances the determination of pulse transit time (PTT) of PPG signals. The design makes use of double reflectance mode based PPG probes, which are placed complementary to each other and on the two sides of a single finger. Furthermore, a novel indicator denoted as Separating Factor is proposed as well. It helps evaluating the performance of ICA with different sensors configuration. This paper then addresses the effectiveness of FastICA in motion artifacts reduction by using the novel method and normal sensors placement method to capture PPG. Results indicate that better independent source separation can be achieved and morphology of PPG signal is perfectly restored when using the method proposed in this paper.

Impulse radio ultra-wide band (IR-UWB) radar has recently emerged as a promising candidate for non-contact monitoring of respiration and heart rate. Different studies have reported various radar based algorithms for estimation of these physiological parameters. The radar can be placed under a subject's mattress as he lays stationary on his back or it can be attached to the ceiling directly above the subject's bed. However, advertent or inadvertent movement on part of the subject and different postures can affect the radar returned signal and also the accuracy of the estimated parameters from it. The detection and analysis of these postural changes can not only lead to improvement in estimation algorithms but also towards prevention of bed sores and ulcers in patients who require periodic posture changes. In this paper, we present an algorithm that detects and quantifies different types of motion events using an under-the-mattress IR-UWB radar. The algorithm also indicates a change in posture after a macro-movement event. Based on the findings of this paper, we anticipate that IR-UWB radar can be used for extracting posture related information in non-clinical enviroments for patients who are bed-ridden.

This paper investigates how to train a recurrent neural network (RNN) using the Levenberg–Marquardt (LM) algorithm as well as how to implement optimal control of a grid-connected converter (GCC) using an RNN. To successfully and efficiently train an RNN using the LM algorithm, a new forward accumulation through time (FATT) algorithm is proposed to calculate the Jacobian matrix required by the LM algorithm. This paper explores how to incorporate FATT into the LM algorithm. The results show that the combination of the LM and FATT algorithms trains RNNs better than the conventional backpropagation through time algorithm. This paper presents an analytical study on the optimal control of GCCs, including theoretically ideal optimal and suboptimal controllers. To overcome the inapplicability of the optimal GCC controller under practical conditions, a new RNN controller with an improved input structure is proposed to approximate the ideal optimal controller. The performance of an ideal optimal controller and a well-trained RNN controller was compared in close to real-life power converter switching environments, demonstrating that the proposed RNN controller can achieve close to ideal optimal control performance even under low sampling rate conditions. The excellent performance of the proposed RNN controller under challenging and distorted system conditions further indicates the feasibility of using an RNN to approximate optimal control in practical applications.

Estimation of neural models based on observed spike timing faces challenges as the amount of recorded units increases, especially when identifying detailed model features. Given that neural regions are generally sparsely connected, input selection is a critical step in model estimation but oftentimes computationally and theoretically challenging. In this paper, we detail an efficient methodology for estimating a sparse, nonlinear dynamical multiple-input, single-output model (MISO) applicable to large-scale (n > 50) single-unit recorded activity. The main contribution of this paper is the complete implementation of a principled group-lasso and local coordinate descent (LCD) algorithm into a generalized Volterra model (GVM) framework to achieve efficient sparse model estimation. Input selection is achieved with group-lasso by simultaneously selecting groups of parameters that are associated with each input. LCD yields efficient computation as the amount of inputs and parameters increase. We investigate and validate the performance of this estimation procedure with the application to a 64 input simulated model.

This paper addresses the problem of automatic detection of tumoral frames in endoscopic capsule videos by using features directly extracted from the color space. We show that tumor can be appropriately discriminated from normal tissue by using only color information histogram measures from the Lab color space and that light saturated regions are usually classified as tumoral regions when color based discriminative procedures are used. These regions are correctly classified if lightening is discarded becoming the tissue classifier based only on the color differences a and b of the Lab color space. While current state of the art systems for small bowel tumor detection usually rely on the processing of the whole frame regarding features extraction this paper proposes the use of fully automatic segmentation in order to select regions likely to contain tumoral tissue. Classification is performed by using Support Vector Machine (SVM) and Multilayer Perceptron (MLP) by using features from color channels a and b of the Lab color space. The proposed algorithm outperforms in more than 5% a series of other algorithms based on features obtained from the higher frequency components selected from Wavelets and Curvelets transforms while saving important computational resources. In a matter of fact the proposed algorithm is more than 25 times faster than algorithms requiring wavelet/curvelet and co-occurrence computations.

This paper presents QRS complex detection algorithm based on dual slope technique, which is suitable for wearable electrocardiogram (ECG) applications. For cardiac patients of different arrhythmias, ECG signals are needed to be monitored over an extensive period of time. Thus, the wearable heart monitoring system needs computationally efficient QRS detection technique with good accuracy. In this paper, a method of QRS detection based on two slopes on both sides of an R peak is presented which is computationally efficient. Based on the slopes, first, a variable measuring steepness is developed, then by introducing an adjustable R-R interval based window and adaptive thresholding techniques, depending on the number of peaks detected in such window, R peaks are detected. The algorithm was evaluated against MIT/BIH arrhythmia database and achieved 99.16% detection rate with sensitivity of 0.9935 and positive predictivity of 0.9981. The method was compared with two widely used R peaks detection algorithms.

The paper presents an analytical model to estimate blood glucose level from measurements made non-invasively and in real time by an antenna strapped to a patient's wrist. Some promising success has been shown by the RIT ETA Lab research group that an antenna's resonant frequency can track, in real time, changes in glucose concentration. Based on an in-vitro study of blood samples of diabetic patients, the paper presents a modified Cole-Cole model that incorporates a factor to represent the change in glucose level. A calibration technique using the input impedance technique is discussed and the results show a good estimation as compared to the glucose meter readings. An alternate calibration methodology has been developed that is based on the shift in the antenna resonant frequency using an equivalent circuit model containing a shunt capacitor to represent the shift in resonant frequency with changing glucose levels. Work under progress is the optimization of the technique with a larger sample of patients.

This paper presents a segmentation methodology of abdominal axial CT images. The aim of the study is to determine the location of mesenteric area from the axial images so the organs enclosed within can be localized precisely for diagnostic purposes. The challenge confronted here is that there is no a certain deterministic shape of abdominal organs. The methodology implemented here utilizes a curvelets stage followed by morphological image processing to achieve a contour emphasized segmentation from the gestalts of surrounding organs. This paper gives a detailed analysis of approach taken with the problems faced and a brief comparison wrt to other wavelet approaches.

With the increasing utilization of X-ray Computed Tomography (CT) in medical diagnosis, obtaining higher quality image with lower exposure to radiation has become a highly challenging task in image processing. In this paper, a novel sparse fusion algorithm is proposed to address the problem of lower Signal to Noise Ratio (SNR) in low dose CT images. Initial fused image is obtained by combining low dose and medium dose images in sparse domain, utilizing the Dual Tree Complex Wavelet Transform (DTCWT) dictionary which is trained by high dose image. And then, the strongly focused image is obtained by determining the pixels of source images which have high similarity with the pixels of the initial fused image. Final denoised image is obtained by fusing strongly focused image and decomposed sparse vectors of source images, thereby preserving the edges and other critical information needed for diagnosis. This paper demonstrates the effectiveness of the proposed algorithm both quantitatively and qualitatively.

This paper presents a computationally efficient QRS detection algorithm for wearable electrocardiogram (ECG) applications based on dual-slope analysis. In general, ECG signals of arrhythmias are pseudo-periodic and contaminated with noises like the patient's contraction muscles, respiration, 60 Hz interference and other types which impede correct QRS detection. To resolve this problem, in this paper, a technique is presented which is based on two slopes on both sides of a peak in ECG signal. Based on these slopes, a variable measuring steepness is developed and R peaks are detected. The algorithm was evaluated against MIT/BIH arrhythmia database and achieved 99.38% detection rate. This method was compared with one of the recently developed dual-slope based QRS detection methods. The results showed that the proposed method has 12.48 times faster runtime than the old method.

Perimetry is essential for identifying visual field defects due to disorders of the eye and brain. However, young children are often unable to reliably perform the preferred method of visual field assessment known as automated static perimetry (ASP). This paper introduces a novel method of ASP specifically developed for children called Saccadic Vector Optokinetic Perimetry (SVOP). SVOP uses eye tracking to detect the natural saccadic eye response of gaze orientation towards visual field stimuli if they are seen. In this paper, the direction and magnitude of a sample of subject gaze responses to visual field stimuli is used to construct a software decision algorithm for use in SVOP. SVOP was clinically evaluated in a group of 24 subjects, comprising children and adults, with and without visual field defects, by comparison with an equivalent test on the Humphrey Field Analyser (HFA). SVOP provides promising visual field test results when compared with the reference HFA test, and has proven extremely useful in detecting visual field defects in children unable to perform traditional ASP.

Performance of robotics based hand-held surgical devices in real-time is mainly dependent on accurate filtering of physiological tremor. The presence of phase delay in sensors (hardware) and filtering (software) processes affects the cancellation accuracy. This paper focuses on developing an estimation algorithm to improve the estimation accuracy in the presence of phase delay for real-time implementations. Moving window based online training approach for least squares-support vector machines (LSSVM) is employed in this paper for tremor estimation. A study is conducted with tremor data recorded from the subjects to analyze the suitability of proposed approach for both single-step and multi-step prediction.

Nowadays the interpretation of the images acquired by wireless endoscopy system is a tedious job for doctors. A viable solution is to construct a map, which is the 2D representation of gastrointestinal (GI) tract to reduce the redundancy of images and improve the understandability of them. The work reported in this paper addresses the problem of the 2D representation of GI tract based on a new wireless Micro-Ball endoscopy system with multiple image sensors. This paper firstly models the problem of constructing the map, and then discusses mainly on the issues of perspective distortion correction, image preprocessing and image registration, which lie in the whole problem. The perspective distortion correction algorithm is realized based on attitude angles, while the image registration is based on phase correlation method (PCM) and scale invariant feature transform (SIFT) combined with particular image preprocessing methods. Based on R channels of images, the algorithm can deal with 26.3% to 100% of image registration when the ratio of overlap varies from 25% to 80%. The performance and effectiveness of the algorithms are verified by experiments.

Wireless body area networks (WBANs) will take on more diverse forms in terms of their sensor combinations and communication protocols as their presence is extended to a greater number of monitoring scenarios. This paper presents an application layer protocol that solves issues caused by sensor nodes that must compete for high speed, real-time communication with the receiver. Such applications emphasize the delivery of large amounts of raw data from different sensor nodes in a time-synchronized manner, rather than channels that experience intermittent operation. An example of tracking pulse wave velocity (PWV) is introduced in this paper, where high-precision PWVs are estimated with the help of timeline recovery and feature extraction processes in MATLAB.

In this paper, we propose a Meta-Cognitive Neuro-Fuzzy Inference System (McFIS) for accurate detection of human actions from video sequences. In this paper, we employ optical flow based features as they can represent information from local pixel level to global object level between two consecutive image planes. The functional relationship between these optical flow based features and action classes is approximated using McFIS classifier. The sequential learning algorithm is developed based on the principles of self-regulation observed in human meta-cognition. McFIS decides on what-to-learn, when-to-learn and how-to-learn based on the knowledge stored in the classifier and the information contained in the new training sample. The sequential learning algorithm of McFIS is controlled and monitored by the meta-cognitive components which uses class-specific and knowledge based criteria along with self-regulatory thresholds to decide on one of the following strategies: a) sample deletion b) sample learning and c) sample reserve. Performance of proposed McFIS based human action recognition system is evaluated using benchmark Weizmann and KTH video sequences. The simulation results are compared with well known support vector machine classifier and also with state-of-the-art action recognition results reported in the literature. The results clearly indicates McFIS action recognition system achieves better performances with minimal computational effort.

In this paper, particle swarm optimization (PSO) algorithm is used to estimate the permittivities of the tissue layers at microwave frequency band. According to the literature, microwave radiometry (MWR) is potentially a promising cancer detection technique. In addition, breast cancer is an appropriate candidate of MWR due to the breast's exclusive physiology. Several algorithms have been evaluated for analyzing the measurement data and solving the inverse scattering problem in MWR, and different levels of accuracy have been reported. In this paper, the potential of PSO in solving this problem is demonstrated at 1-2.25 GHz. Two distinct algorithms are developed for the two considered scenarios. In the first scenario, we assume no a priori knowledge of the tissue under the test, whereas, in the second scenario, a priori knowledge is assumed. It is noteworthy that, there are only a few research articles studying PSO for permittivity estimation. However, since these studies underestimate the loss encountered by the test samples, the methods are not valid for body tissue case. Here, measurement-based loss coefficients, reported in the existing literature, are included in the calculations. It is shown that the algorithm converges relatively fast, and, distinguishes between different tissues with an acceptable accuracy.

A signal processor/compressor dedicated to implantable neural recording microsystems is presented. Signal compression is performed based on Haar wavelet. It is shown in this paper that, compared to other mathematical transforms already used for this purpose, compression of neural signals using this type of wavelet transform can be of almost the same quality, while demanding less circuit complexity and smaller silicon area. Designed in a 0.13-μm standard CMOS process, the 64-channel 8-bit signal processor reported in this paper occupies 113 μm x 110 μm of silicon area. It operates under a 1.8-V supply voltage at a master clock frequency of 3.2 MHz.

This paper presents a time-varying cortical connectivity analysis for to a source localization approach within the inner layers of the newborn brain. (Additional details can be found in the comprehensive book on Time-Frequency Signal Analysis and Processing (see http://www.elsevier.com/locate/isbn/0080443354). In addition, the most recent upgrade of the original software package that calculates Time-Frequency Distributions and Instantaneous Frequency estimators can be downloaded from the web site: www.time-frequency.net. This was the first software developed in the field, and it was first released publicly in 1987 at the 1st ISSPA conference held in Brisbane, Australia, and then continuously updated). Multivariate Granger causality in the timefrequency domain as a representation of time-varying cortical connectivity in the brain has been investigated for the adult case. This is, however, not the case in newborns as the nature of the transient changes in the newborn EEG is different from that of adults. This paper aims to evaluate the performance of the time-varying versions of the two popular Granger causality measures, namely Partial Directed Coherence (PDC) and direct Directed Transfer Function (dDTF). The parameters of the time-varying AR, that models the inter-channel interactions, are estimated using Dual Extended Kalman Filter (DEKF) as it accounts for both non-stationarity and non-linearity behaviors of the EEG. Using simulated data, we show that fast changing cortical connectivity between channels can be measured more accurately using the time-varying PDC. The performance of the time-varying PDC is also tested on a neonatal EEG exhibiting seizure.

This paper introduces a kernel adaptive filter implemented with stochastic gradient on temporal differences, kernel Temporal Difference (TD)(λ), to estimate the state-action value function in reinforcement learning. The case λ=0 will be studied in this paper. Experimental results show the method's applicability for learning motor state decoding during a center-out reaching task performed by a monkey. The results are compared to the implementation of a time delay neural network (TDNN) trained with backpropagation of the temporal difference error. From the experiments, it is observed that kernel TD(0) allows faster convergence and a better solution than the neural network.

This paper describes the initial groundwork car- ried out as part of the European Commission funded Transatlantic Tumor Model Repositories project, to develop a new markup language for computational cancer modelling, TumorML. In this paper we describe the motivations for such a language, arguing that current state-of-the-art biomodelling languages are not suited to the cancer modelling domain. We go on to describe the work that needs to be done to develop TumorML, the conceptual design, and a description of what existing markup languages will be used to compose the language specification.

This paper is the continuation of a work presented at ICORR 07, in which we discussed the possibility of improving eye-hand coordination in children diagnosed with this problem, using a robotic mapping from a haptic user interface to a virtual environment. Our goal is to develop, implement and refine a system that will assess and improve the eye-hand coordination and grip strength in children diagnosed with poor graphomotor skills. A detailed analysis of patters (e.g., labyrinths, letters and angles) was conducted in order to select three very distinguishable levels of difficulty that could be included in the system, and which would yield the greatest benefit in terms of assessment of coordination and strength issues as well as in training. Support algorithms (position, force, velocity, inertia and viscosity) were also developed and incorporated into the tasks in order to introduce general computer assistance to the mapping of the user's movements to the computer screen without overriding the user's commands to the robotic device. In order to evaluate performance (given by %accuracy and time) of the executed tasks, a sophisticated evaluation function was designed based on image analysis and edge detection algorithms. This paper presents the development of the haptic tasks, the various assistance algorithms, the description of the evaluation function and the results of a study implemented at the Motor Development Clinic at Cal Poly Pomona. The results (Accuracy and Time) of this function are currently being used as inputs to an Intelligent Decision Support System (described in [5]), which in turn, suggests the next task to be executed by the subject based on his/her performance.

Different biological labs conduct similar experiments on same diseases. It is highly desirable to have a better model based on more experimental results than that on a single result. In this paper, we propose a method for integrating microarray data from multiple sources for building classification models. To test the method, we use three real world microarray data sets from different labs with different experimental devices and environments. Although microarray data is well known for its inconsistencies across labs, we demonstrate that it is possible to build consistent models using data sets from multiple labs. We report our method, experimental results and observations in the paper. Refereed/Peer-reviewed

Independent Component Analysis (ICA) is a very common instantiation of the Blind Source Separation (BSS) problem. In the context of biomedical signal analysis, ICA is generally applied to multi-channel recordings of physiological phenomena in order to de-noise and extract meaningful information underlying the recordings. This paper assesses the Spatio-Temporal ICA (ST-ICA) framework, which uses both spatial and temporal information derived from multi-channel time-series to extract underlying sources. In contrast, the standard implementation of the ICA algorithm generally uses only limited spatial information to inform the separation process. One of the major steps in the implementation of any ICA algorithm is the selection of relevant components from the many ICA usually returns. With ST-ICA there is a rich data-set of components exhibiting spatial as well as temporal/spectral information that could be used to identify the underlying process subspaces extracted by the ST-ICA algorithm. This paper highlights the methodology for performing ST-ICA and assesses the possible ways in which process subspace identification may take place.

This paper presents a new application of the Particle Swarm Optimization (PSO) algorithm to optimize Mel Frequency Cepstrum Coefficients (MFCC) parameters, in order to extract an optimal feature set for diagnosis of hypothyroidism in infants using Multi-Layer Perceptrons (MLP) neural network. MFCC features is influenced by the number of filter banks (f(b)) and the number of coefficients (n(c)) used. These parameters are critical in representation of the features as they affect the resolution and dimensionality of the features. In this paper, the PSO algorithm was used to optimize the values of f(b) and n(c). The MFCC features based on the PSO optimization were extracted from healthy and unhealthy infant cry signals and used to train MLP in the classification of hypothyroid infant cries. The results indicate that the PSO algorithm could determine the optimum combination of f(b) and n(c) that produce the best classification accuracy of the MLP.

Objective: To remove blood from an incision and find the incision spot is a key task during surgery, or else over discharge of blood will endanger a patient’s life. However, the repetitive manual blood removal involves plenty of workload contributing fatigue of surgeons. Thus, it is valuable to design a robotic system which can automatically remove blood on the incision surface. Methods: In this paper, we design a robotic system to fulfill the surgical task of the blood removal. The system consists of a pair of dual cameras, a 6-DoF robotic arm, an aspirator whose handle is fixed to a robotic arm, and a pump connected to the aspirator. Further, a path-planning algorithm is designed to generate a path, which the aspirator tip should follow to remove blood. Results: In a group of simulating bleeding experiments on ex vivo porcine tissue, the contour of the blood region is detected, and the reconstructed spatial coordinates of the detected blood contour is obtained afterward. The BRR robot cleans thoroughly the blood running out the incision. Conclusions: This study contributes the first result on designing an autonomous blood removal medical robot. The skill of the surgical blood removal operation, which is manually operated by surgeons nowadays, is alternatively grasped by the proposed BRR medical robot.

The ballistocardiogram (BCG) signal represents the movements of the body in response to cardiac ejection of blood. Recently, many groups have developed low-cost instrumentation for facilitating BCG measurement in the home. The standard method used in the literature for estimating the BCG pulse response has generally been ensemble averaging over several beats. Unfortunately, since the BCG pulse response is likely longer than a typical heartbeat interval, this standard approach does not yield a full-length estimate of the response. This paper describes a simple, novel algorithm for estimating the full-length BCG pulse response using the R-wave timing of a simultaneously acquired electrocardiogram (ECG). With this pulse response, the full signal can be reconstructed, enabling the analysis of slow transient effects in the BCG signal, and of the measurement noise. Additionally, while this paper focuses only on the BCG signal, the same algorithm could be applied to other biomedical signals such as the phonocardiogram or impedance cardiogram, particularly when the heartbeat interval is shorter than the duration of the cpulse response.

Dynamic infrared imaging is a promising technique in breast oncology. In this paper, a quantum well infrared photodetector infrared camera is used to acquire a sequence of consecutive thermal images of the patientpsilas breast for 10 s. Information on the local blood perfusion is obtained from the spectral analysis of the time series at each image pixel. Due to respiratory and motion artifacts, the direct comparison of the temperature values that a pixel assumes along the sequence becomes difficult. In fact, the small temperature changes due to blood perfusion, of the order of 10-50 mK, which constitute the signal of interest in the time domain, are superimposed onto large temperature fluctuations due to the subjectpsilas motion, which represent noise. To improve the time series S/N, and as a consequence, enhance the specificity and sensitivity of the dynamic infrared examination, it is important to realign the thermal images of the acquisition sequence, thus reducing motion artifacts. In a previous study, we demonstrated that a registration algorithm based on fiducial points is suitable to both clinical applications and research, when associated with a proper set of skin markers. In this paper, we quantitatively evaluate the performance of different marker sets by means of a model that allows for estimating the S/N increment due to registration, and we conclude that a 12-marker set is a good compromise between motion artifact reduction and the time required to prepare the patient.

Closed Angle Glaucoma causes an increase in fluid pressure in the eye due to decreased fluid flow between the iris and the cornea. To diagnose Glaucoma, ultrasound images of the eye are taken and a technician studies the images and takes measurements. This process is time consuming, and results vary among technicians. In a previous paper, we developed an image processing algorithm to extract several key landmark features in the iris. In this paper, we extend this algorithm to calculate the complete set of diagnostic parameters from ultrasound biomicroscopy. A robust algorithm is described to calculate the parameters which describe the thickness of iris; ID1, ID2, and ID3. The method is tested on 89 clinical images and compared to the values identified by the technician. Results indicate good agreement in 97% of cases.

Multi-slice Computed Tomography (MSCT) is an important medical imaging tool that provides dynamic three-dimensional (3D) volume data of the heart for diagnosis of various cardiac diseases. Due to the huge amount of data in MSCT, manual identification, segmentation and tracking of various parts of the heart are very labor intensive and inefficient. In this paper, we introduce a semi-automatic method for robustly segmenting the endocardium surface from cardiac MSCT images. A level set approach is adopted to define a flexible and powerful interface for capturing the complex anatomical structure of the heart. A novel speed function based on clustering the image intensities of the region of interest and the background is proposed for use with the level set method. The method introduced in this paper has the advantages of simple initialization and being capable of segmenting the blood pool with non-homogeneous intensities. Experiments on real data using the proposed speed function have been carried out with 2D, 3D and 4D implementations of the level sets respectively, and comparisons in terms of computational speed and segmentation results are presented.

The number of people that now go on to develop Alzheimer's disease (AD) and other types of dementia is rapidly rising. For maximum benefits from new treatments, the disease should be diagnosed as early as possible, but this is difficult with current clinical criteria. Potentially, the EEG can serve as an objective, first line of decision support tool to improve diagnosis. It is non-invasive, widely available, low-cost and could be carried out rapidly in the high-risk age group that will develop AD. Changes in the EEG due to the dementing process could be quantified as an index or marker. In this paper, we investigate two information theoretic methods (Tsallis entropy and universal compression algorithm) as a way to generate potentially robust markers from the EEG. The hypothesis is that the information theoretic makers for AD are significantly different to those of normal subjects. An attraction of the information theoretic approach is that, unlike most existing methods, there may be a natural link between the underlying ideas of information theoretic methods, the physiology of AD and its impact on brain functions. Data compression has not been investigated as a means of generating EEG markers before and is attractive because it does not require a priori knowledge of the source model. In this paper, we focus on the LZW algorithm because of its sound theoretical foundation. We used the LZW algorithm and Tsallis model to compute the markers (compression ratios and normalized entropies, respectively) from two EEG datasets. The results show that the information theoretic methods can be used to compute EEG markers for AD.

An electrocardiogram (ECG) is an important and commonly used diagnostic aid in cardiovascular disease diagnosis. Physicians routinely perform diagnosis by a simple visual examination of ECG waveform shapes. In this paper, we address the problem of shape-based retrieval of ECG recordings, both digital and scanned from paper, to infer similarity in diagnosed diseases. Specifically, we use the knowledge of ECG recording structure to segment and extract curves representing various recording channels from ECG images. We then present a method of capturing the perceptual shape similarity of ECG waveforms by combining shape matching with dynamic time warping. The shape similarity of each recording channel is combined to develop an overall shape similarity measure between ECG recordings. Results are presented that demonstrate the method on shape-based matching of various cardiovascular diseases.

The segmentation and landmark identification in infrared images of the human body are key steps in a computerized processing of large database of thermal images. The segmentation task is especially challenging due to specific characteristics of thermal images. Few papers deal with segmentation techniques for clinical infrared images and available segmentation methods (e.g. for breast or military thermal images) do not perform well on other types of images. This paper presents a few strategies for the automated segmentation and registration of anatomical landmarks on thermal images of arms and hands. The segmentation method is based on mathematical morphological operations and simple rule based processing easily available through prior knowledge about the objects of interest.

This paper outlines the development and initial testing of a new hybrid computer cursor control system based on Eye Gaze Tracking (EGT) and electromyogram (EMG) processing for hands-free control of the computer cursor. The ultimate goal of the system is to provide an efficient computer interaction mechanism for individuals with severe motor disabilities (or specialized operators whose hands are committed to other tasks, such as surgeons, pilots, etc.) The paper emphasizes the enhancements that have been made on different areas of the architecture, with respect to a previous prototype developed by our group, and demonstrates the performance improvement verified for some of the enhancements.

This paper presents a complete open software environment for control, data processing and assessment of visual experiments. Visual experiments are widely used in research on human perception physiology and the results are applicable to various visual information-based man-machine interfacing, human-emulated automatic visual systems or scanpath-based learning of perceptual habits. The toolbox is designed for Matlab platform and supports infra-red reflection-based eyetracker in calibration and scanpath analysis modes. Toolbox procedures are organized in three layers: the lower one, communicating with the eyetracker output file, the middle detecting scanpath events on a physiological background and the one upper consisting of experiment schedule scripts, statistics and summaries. Several examples of visual experiments carried out with use of the presented toolbox complete the paper.

This paper describes a novel user-guided method for grouping fibers from diffusion tensor MRI tractography into bundles. The method finds fibers, that passing through user-defined ROIs, still fit to the underlying data model given by the diffusion tensor. This is achieved by filtering the data and the ROIs with a kernel derived from a geodesic metric between tensors. A standard approach using binary decisions defining tracts passing through ROIs is critically dependent on ROIs that includes all trace lines of interest. The method described in this paper uses a softer decision mechanism through a kernel which enables grouping of bundles driven less exact, or even single point, ROIs. The method analyzes the responses obtained from the convolution with a kernel function along the fiber with the ROI data. Results in real data shows the feasibility of the approach to fiber bundling.

Functional magnetic resonance imaging (fMRI) is a\ud rapidly growing technique for studying the brain in\ud action. Since its creation [1], [2], cognitive scientists\ud have been using fMRI to understand how we remember,\ud manipulate, and act on information in our environment.\ud Working with magnetic resonance physicists, statisticians, and\ud engineers, these scientists are pushing the frontiers of knowledge\ud of how the human brain works.\ud The design and analysis of single-subject fMRI studies\ud has been well described. For example, [3], chapters 10\ud and 11 of [4], and chapters 11 and 14 of [5] all give accessible\ud overviews of fMRI methods for one subject. In contrast,\ud while the appropriate manner to analyze a group of\ud subjects has been the topic of several recent papers, we do\ud not feel it has been covered well in introductory texts and\ud review papers. Therefore, in this article, we bring together\ud old and new work on so-called group modeling of fMRI\ud data using a consistent notation to make the methods more\ud accessible and comparable.

Pain sensation is essential for survival, since it draws attention to physical threat to the body. Pain assessment is usually done through self-reports. However, self-assessment of pain is not available in the case of noncommunicative patients, and therefore, observer reports should be relied upon. Observer reports of pain could be prone to errors due to subjective biases of observers. Moreover, continuous monitoring by humans is impractical. Therefore, automatic pain detection technology could be deployed to assist human caregivers and complement their service, thereby improving the quality of pain management, especially for noncommunicative patients. Facial expressions are a reliable indicator of pain, and are used in all observer-based pain assessment tools. Following the advancements in automatic facial expression analysis, computer vision researchers have tried to use this technology for developing approaches for automatically detecting pain from facial expressions. This paper surveys the literature published in this field over the past decade, categorizes it, and identifies future research directions. The survey covers the pain datasets used in the reviewed literature, the learning tasks targeted by the approaches, the features extracted from images and image sequences to represent pain-related information, and finally, the machine learning methods used.

In this paper, we have presented Turtulebot-assisted instantaneous heart rate (HR) estimator using camera based remote photoplethysmography. We used a Turtlebot with a camera to record human face. For the face detection, we used Haar Cascade algorithm. To increase the accuracy of the HR estimation, we combined a plane-orthogonal-to-skin (POS) model with finite state machine (FSM) framework. By combining POS and FSM framework, we achieved 1.08 bpm of MAE, which is the lowest error comparing to the state-of-art methods.