Biomedical subjects

Steven Haker

Publications and source records attributed to Steven Haker.

12 recordsLinked to original sources

PubMed2006

Deformable structure registration of bladder through surface mapping.

Cumulative dose distributions in fractionated radiation therapy depict the dose to normal tissues and therefore may permit an estimation of the risk of normal tissue complications. However, calculation of these distributions is highly challenging because of interfractional changes in the geometry of patient anatomy. This work presents an algorithm for deformable structure registration of the bladder and the verification of the accuracy of the algorithm using phantom and patient data. In this algorithm, the registration process involves conformal mapping of genus zero surfaces using finite element analysis, and guided by three control landmarks. The registration produces a correspondence between fractions of the triangular meshes used to describe the bladder surface. For validation of the algorithm, two types of balloons were inflated gradually to three times their original size, and several computerized tomography (CT) scans were taken during the process. The registration algorithm yielded a local accuracy of 4 mm along the balloon surface. The algorithm was then applied to CT data of patients receiving fractionated high-dose-rate brachytherapy to the vaginal cuff, with the vaginal cylinder in situ. The patients' bladder filling status was intentionally different for each fraction. The three required control landmark points were identified for the bladder based on anatomy. Out of an Institutional Review Board (IRB) approved study of 20 patients, 3 had radiographically identifiable points near the bladder surface that were used for verification of the accuracy of the registration. The verification point as seen in each fraction was compared with its predicted location based on affine as well as deformable registration. Despite the variation in bladder shape and volume, the deformable registration was accurate to 5 mm, consistently outperforming the affine registration. We conclude that the structure registration algorithm presented works with reasonable accuracy and provides a means of calculating cumulative dose distributions.

Algorithms↗

PubMed2005

Automated atlas-based clustering of white matter fiber tracts from DTMRI.

A new framework is presented for clustering fiber tracts into anatomically known bundles. This work is motivated by medical applications in which variation analysis of known bundles of fiber tracts in the human brain is desired. To include the anatomical knowledge in the clustering, we invoke an atlas of fiber tracts, labeled by the number of bundles of interest. In this work, we construct such an atlas and use it to cluster all fiber tracts in the white matter. To build the atlas, we start with a set of labeled ROIs specified by an expert and extract the fiber tracts initiating from each ROI. Affine registration is used to project the extracted fiber tracts of each subject to the atlas, whereas their B-spline representation is used to efficiently compare them to the fiber tracts in the atlas and assign cluster labels. Expert visual inspection of the result confirms that the proposed method is very promising and efficient in clustering of the known bundles of fiber tracts.

Algorithms↗

PubMed2005

Harmonic skeleton guided evaluation of stenoses in human coronary arteries.

This paper presents a novel approach that three-dimensionally visualizes and evaluates stenoses in human coronary arteries by using harmonic skeletons. A harmonic skeleton is the center line of a multi-branched tubular surface extracted based on a harmonic function, which is the solution of the Laplace equation. This skeletonization method guarantees smoothness and connectivity and provides a fast and straightforward way to calculate local cross-sectional areas of the arteries, and thus provides the possibility to localize and evaluate coronary artery stenosis, which is a commonly seen pathology in coronary artery disease.

Algorithms↗

PubMed2005

Combining classifiers using their receiver operating characteristics and maximum likelihood estimation.

In any medical domain, it is common to have more than one test (classifier) to diagnose a disease. In image analysis, for example, there is often more than one reader or more than one algorithm applied to a certain data set. Combining of classifiers is often helpful, but determining the way in which classifiers should be combined is not trivial. Standard strategies are based on learning classifier combination functions from data. We describe a simple strategy to combine results from classifiers that have not been applied to a common data set, and therefore can not undergo this type of joint training. The strategy, which assumes conditional independence of classifiers, is based on the calculation of a combined Receiver Operating Characteristic (ROC) curve, using maximum likelihood analysis to determine a combination rule for each ROC operating point. We offer some insights into the use of ROC analysis in the field of medical imaging.

Artificial Intelligence↗

PubMed2005

Mass preserving registration for heart MR images.

This paper presents a new algorithm for non-rigid registration between two doubly-connected regions. Our algorithm is based on harmonic analysis and the theory of optimal mass transport. It assumes an underlining continuum model, in which the total amount of mass is exactly preserved during the transformation of tissues. We use a finite element approach to numerically implement the algorithm.

Algorithms↗

PubMed2005

Multiscale 3D shape analysis using spherical wavelets.

Shape priors attempt to represent biological variations within a population. When variations are global, Principal Component Analysis (PCA) can be used to learn major modes of variation, even from a limited training set. However, when significant local variations exist, PCA typically cannot represent such variations from a small training set. To address this issue, we present a novel algorithm that learns shape variations from data at multiple scales and locations using spherical wavelets and spectral graph partitioning. Our results show that when the training set is small, our algorithm significantly improves the approximation of shapes in a testing set over PCA, which tends to oversmooth data.

Algorithms↗

PubMed2005

Flattening maps for the visualization of multibranched vessels.

In this paper, we present two novel algorithms which produce flattened visualizations of branched physiological surfaces, such as vessels. The first approach is a conformal mapping algorithm based on the minimization of two Dirichlet functionals. From a triangulated representation of vessel surfaces, we show how the algorithm can be implemented using a finite element technique. The second method is an algorithm which adjusts the conformal mapping to produce a flattened representation of the original surface while preserving areas. This approach employs the theory of optimal mass transport. Furthermore, a new way of extracting center lines for vessel fly-throughs is provided.

Algorithms↗

PubMed2004

The tale of Phineas Gage, digitally remastered.

The injury of Phineas Gage has fueled research on and fascination with the localization of cerebral functions in the past century and a half. Most physicians and anatomists believed that Gage sustained a largely bilateral injury to the frontal lobes. However, previous studies seem to have overlooked a few less obvious, but essential details. This has led us to reanalyze the injury using three-dimensional reconstruction and quantitative computer-aided techniques and to propose a new biomechanical model, in order to determine the location and extent of the injury and explain Gage's improbable survival. Unlike previous studies on this subject, our findings are based on computer-generated three-dimensional reconstructions of a thin-slice computed tomography scan (CAT) of Phineas Gage's skull. The results of our image analysis were corroborated with the clinical findings, thoroughly recorded by Dr. Harlow in 1848, as well as with a systematic examination of the original skull specimen. Our results show that the cerebral injury was limited to the left frontal lobe, did not extend to the contralateral side, did not affect the ventricular system, and did not involve vital intracranial vascular structures. Although modern neuroscience has perhaps outgrown the speculations prompted by this famous case, it is still a living part of the medical folklore and education. Setting the record straight based on clinical reasoning, observation of the physical evidence, and sound quantitative computational methods is more than mere minutia and of interest for the broad medical community.

Adult↗

PubMed2004

Workflow modeling and analysis of computer guided prostate brachytherapy under MR imaging control.

We demonstrate that classical Business Process Reengineering (BPR) methods can be successfully applied to Computer Aided Surgery while increasing safety and efficiency of the overall procedure through an integrated Workflow Management System. Computer guided Prostate Brachytherapy, as a sophisticated treatment by an interdisciplinary team, is perfectly suited to apply our method. Detailed suggestions for improvement of the whole procedure could be derived by our modified BPR method.

Brachytherapy↗

PubMed2003

A unified variational approach to denoising and bias correction in MR.

We propose a novel bias correction method for magnetic resonance (MR) imaging that uses complementary body coil and surface coil images. The former are spatially homogeneous but have low signal intensity; the latter provide excellent signal response but have large bias fields. We present a variational framework where we optimize an energy functional to estimate the bias field and the underlying image using both observed images. The energy functional contains smoothness-enforcing regularization for both the image and the bias field. We present extensions of our basic framework to a variety of imaging protocols. We solve the optimization problem using a computationally efficient numerical algorithm based on coordinate descent, preconditioned conjugate gradient, half-quadratic regularization, and multigrid techniques. We show qualitative and quantitative results demonstrating the effectiveness of the proposed method in producing debiased and denoised MR images.

Algorithms↗

PubMed2003

Multislice CT colonography: current status and limitations.

CT colonography (CTC) is a promising method for colorectal screening providing a full structural evaluation of the entire colon and gaining in popularity due to a superior safety profile, a low rate of complications, and high patient acceptance. Multislice CT (MSCT) has further improved the diagnostic potential of CTC by generating high-resolution CT images of the abdomen and pelvis in shorter acquisition times than was previously possible. Over the past year, multiple studies have been published on every aspect of CTC including techniques, image display, image reconstruction, and clinical trial results assessing the feasibility of CTC as a screening tool. Yet despite increasing clinical use, the appropriate role of CTC in colorectal cancer screening remains undefined and barriers to widespread adoption remain. In particular, though the test is generally regarded as easy to perform, accurate interpretation requires a steep learning curve. While several large studies have found high sensitivity and specificity, the accuracy of CTC in a screening population has yet to be verified and almost no health insurance plans reimburse for its use in colorectal screening. Ongoing research in computer-aided detection and new software tools, however, have the potential to increase accuracy and ease of interpretation significantly, accelerating its acceptance as a colorectal screening tool.

Colon↗

PubMed2003

Detection of prostate cancer by integration of line-scan diffusion, T2-mapping and T2-weighted magnetic resonance imaging; a multichannel statistical classifier.

A multichannel statistical classifier for detecting prostate cancer was developed and validated by combining information from three different magnetic resonance (MR) methodologies: T2-weighted, T2-mapping, and line scan diffusion imaging (LSDI). From these MR sequences, four different sets of image intensities were obtained: T2-weighted (T2W) from T2-weighted imaging, Apparent Diffusion Coefficient (ADC) from LSDI, and proton density (PD) and T2 (T2 Map) from T2-mapping imaging. Manually segmented tumor labels from a radiologist, which were validated by biopsy results, served as tumor "ground truth." Textural features were extracted from the images using co-occurrence matrix (CM) and discrete cosine transform (DCT). Anatomical location of voxels was described by a cylindrical coordinate system. A statistical jack-knife approach was used to evaluate our classifiers. Single-channel maximum likelihood (ML) classifiers were based on 1 of the 4 basic image intensities. Our multichannel classifiers: support vector machine (SVM) and Fisher linear discriminant (FLD), utilized five different sets of derived features. Each classifier generated a summary statistical map that indicated tumor likelihood in the peripheral zone (PZ) of the prostate gland. To assess classifier accuracy, the average areas under the receiver operator characteristic (ROC) curves over all subjects were compared. Our best FLD classifier achieved an average ROC area of 0.839(+/-0.064), and our best SVM classifier achieved an average ROC area of 0.761(+/-0.043). The T2W ML classifier, our best single-channel classifier, only achieved an average ROC area of 0.599(+/-0.146). Compared to the best single-channel ML classifier, our best multichannel FLD and SVM classifiers have statistically superior ROC performance (P=0.0003 and 0.0017, respectively) from pairwise two-sided t-test. By integrating the information from multiple images and capturing the textural and anatomical features in tumor areas, summary statistical maps can potentially aid in image-guided prostate biopsy and assist in guiding and controlling delivery of localized therapy under image guidance.

Algorithms↗