A system using a wide-slot beam and simple antiscatter grids or slots has been designed to provide a localized map of tissue type that could be overlaid on the simultaneous conventional transmission image to provide an inexpensive, low dose adjunct to conventional screening mammography. Depth information is obtainable from the stereoscopic viewing angles. The system was demonstrated to produce observable contrast between adipose tissue and a phantom chosen to mimic carcinoma at an exposure comparable with screening mammography. Imaging data was collected over a range of system parameters to optimize contrast and to allow verification of simulation modeling.

The purpose of this study was to use three-dimensional (3-D) printing techniques to construct liver and brain phantoms having realistic pathologies, anatomic structures, and heterogeneous backgrounds. Patient liver and head computed tomography (CT) images were segmented into tissue, vessels, liver lesion, white and gray matter, and cerebrospinal fluid (CSF). Stereolithography files of each object were created and imported into a commercial 3-D printer. Printing materials were assigned to each object after test scans, which showed that the printing materials had CT numbers ranging from 70 to 121 HU at 120 kV. Printed phantoms were scanned on a CT scanner and images were evaluated. CT images of the liver phantom had measured CT numbers of 77.8 and 96.6 HU for the lesion and background, and 137.5 to 428.4 HU for the vessels channels, which were filled with iodine solutions. The difference in CT numbers between lesions and background (18.8 HU) was representative of the low-contrast values needed for optimization tasks. The liver phantom background was evaluated with Haralick features and showed similar texture between patient and phantom images. CT images of the brain phantom had CT numbers of 125, 134, and 108 HU for white matter, gray matter, and CSF, respectively. The CT number differences were similar to those in patient images.

A microwave imaging system has been developed as a clinical diagnostic tool operating in the 3- to 8-GHz region using multistatic data collection. A total of 86 patients recruited from a symptomatic breast care clinic were scanned with a prototype design. The resultant three-dimensional images have been compared “blind” with available ultrasound and mammogram images to determine the detection rate. Images show the location of the strongest signal, and this corresponded in both older and younger women, with sensitivity of <74%, which was found to be maintained in dense breasts. The pathway from clinical prototype to clinical evaluation is outlined.

Dynamic perfusion imaging can provide the morphologic details of the scanned organs as well as the dynamic information of blood perfusion. However, due to the polyenergetic property of the x-ray spectra, beam hardening effect results in undesirable artifacts and inaccurate CT values. To address this problem, this study proposes a segmentation-free polyenergetic dynamic perfusion imaging algorithm (pDP) to provide superior perfusion imaging. Dynamic perfusion usually is composed of two phases, i.e., a precontrast phase and a postcontrast phase. In the precontrast phase, the attenuation properties of diverse base materials (e.g., in a thorax perfusion exam, base materials can include lung, fat, breast, soft tissue, bone, and metal implants) can be incorporated to reconstruct artifact-free precontrast images. If patient motions are negligible or can be corrected by registration, the precontrast images can then be employed as a priori information to derive linearized iodine projections from the postcontrast images. With the linearized iodine projections, iodine perfusion maps can be reconstructed directly without the influence of various influential factors, such as iodine location, patient size, x-ray spectrum, and background tissue type. A series of simulations were conducted on a dynamic iodine calibration phantom and a dynamic anthropomorphic thorax phantom to validate the proposed algorithm. The simulations with the dynamic iodine calibration phantom showed that the proposed algorithm could effectively eliminate the beam hardening effect and enable quantitative iodine map reconstruction across various influential factors. The error range of the iodine concentration factors (HU·ml/mg) was reduced from [−8.01,14.77] for filtered back-projection (FBP) to [−0.05,0.08] for pDP. The quantitative results of the simulations with the dynamic anthropomorphic thorax phantom indicated that the maximum error of iodine concentrations can be reduced from 1.1  mg/cc for FBP to less than 0.1  mg/cc for pDP, which suggested that the proposed algorithm could not only effectively eliminate beam hardening artifacts but also significantly reduce the influence of the metal artifacts and accurately reconstruct the iodine map regardless of the influential factors. A segmentation-free polyenergetic dynamic perfusion imaging algorithm was proposed and validated via simulations. This method can accurately reconstruct artifact-free iodine maps for quantitative analyses.

Conventional mammography can suffer from poor contrast between healthy and cancerous tissues due to the small difference in attenuation properties. Coherent scatter slot scan imaging is an imaging technique which provides additional information and is compatible with conventional mammography. A Monte Carlo simulation of coherent scatter slot scan imaging was performed to assess its performance and provide system optimization. Coherent scatter could be exploited using a system similar to conventional slot scan mammography system with antiscatter grids tilted at the characteristic angle of cancerous tissues. System optimization was performed across several parameters, including source voltage, tilt angle, grid distances, grid ratio, and shielding geometry. The simulated carcinomas were detectable for tumors as small as 5 mm in diameter, so coherent scatter analysis using a wide-slot setup could be promising as an enhancement for screening mammography. Employing coherent scatter information simultaneously with conventional mammography could yield a conventional high spatial resolution image with additional coherent scatter information.

Quantification of sun-related changes in conjunctival ultraviolet autofluorescence (CUVAF) images is a subjective and tedious task, in which reproducibility of results is difficult. Thus, we have developed a semiautomatic method in MATLAB^® to analyze CUVAF images retrospectively. The algorithm was validated on 200 images from 50 randomly selected participants from the Western Australian Pregnancy Cohort (Raine) study 20-year follow-up assessment, in which CUVAF area measurements were available from previous manual analysis. Algorithm performance was compared to manual measurements and yielded better than 95% correspondence in both intra- and interobserver agreement. Furthermore, the semiautomatic method significantly reduced analysis time by 50%.

Detecting the position of retinal structures, including the fovea center and macula, in retinal images plays a key role in diagnosing eye diseases such as optic nerve hypoplasia, amblyopia, diabetic retinopathy, and macular edema. However, current detection methods are unreliable for infants or certain ethnic populations. Thus, a methodology is proposed here that may be useful for infants and across ethnicities that automatically localizes the fovea center and segments the macula on digital fundus images. First, dark structures and bright artifacts are removed from the input image using preprocessing operations, and the resulting image is transformed to polar space. Second, the fovea center is identified, and the macula region is segmented using the proposed dynamic identification and classification of edges (DICE) model. The performance of the method was evaluated using 1200 fundus images obtained from the relatively large, diverse, and publicly available Messidor database. In 96.1% of these 1200 cases, the distance between the fovea center identified manually by ophthalmologists and automatically using the proposed method remained within 0 to 8 pixels. The dice similarity index comparing the manually obtained results with those of the model for macula segmentation was 96.12% for these 1200 cases. Thus, the proposed method displayed a high degree of accuracy. The methodology using the DICE model is unique and advantageous over previously reported methods because it simultaneously determines the fovea center and segments the macula region without using any structural information, such as optic disc or blood vessel location, and it may prove useful for all populations, including infants.

Recent findings indicate a strong correlation between the risk of future heart disease and the volume of adipose tissue inside of the pericardium. So far, large-scale studies have been hindered by the fact that manual delineation of the pericardium is extremely time-consuming and that existing methods for automatic delineation lack accuracy. An efficient and fully automatic approach to pericardium segmentation and epicardial fat volume (EFV) estimation is presented, based on a variant of multi-atlas segmentation for spatial initialization and a random forest classifier for accurate pericardium detection. Experimental validation on a set of 30 manually delineated computer tomography angiography volumes shows a significant improvement on state-of-the-art in terms of EFV estimation [mean absolute EFV difference: 3.8 ml (4.7%), Pearson correlation: 0.99] with run times suitable for large-scale studies (52 s). Further, the results compare favorably with interobserver variability measured on 10 volumes.

Identification of the basal slice in cardiac imaging is a key step to measuring the ejection fraction of the left ventricle. Despite all the effort placed on automatic cardiac segmentation, basal slice identification is routinely performed manually. Manual identification, however, suffers from high interobserver variability. As a result, an automatic algorithm for basal slice identification is required. Guidelines published in 2013 identify the basal slice based on the percentage of myocardium surrounding the blood cavity in the short-axis view. Existing methods, however, assume that the basal slice is the first short-axis view slice below the mitral valve and are consequently at times identifying the incorrect short-axis slice. Correct identification of the basal slice under the Society for Cardiovascular Magnetic Resonance guidelines is challenging due to the poor image quality and blood movement during image acquisition. This paper proposes an automatic tool that utilizes the two-chamber view to determine the basal slice while following the guidelines. To this end, an active shape model is trained to segment the two-chamber view and create temporal binary profiles from which the basal slice is identified. From the 51 tested cases, our method obtains 92% and 84% accurate basal slice detection for the end-systole and the end-diastole, respectively.

Extraction of bone contours from x-ray radiographs plays an important role in joint space width assessment, preoperative planning, and kinematics analysis. We present a robust segmentation method to accurately extract the distal femur and proximal tibia in knee radiographs of varying image quality. A spectral clustering method based on the eigensolution of an affinity matrix is utilized for x-ray image denoising. An active shape model-based segmentation method is employed for robust and accurate segmentation of the denoised x-ray images. The performance of the proposed method is evaluated with x-ray images from the public-use dataset(s), the osteoarthritis initiative, achieving a root mean square error of 0.48±0.13  mm for femur and 0.53±0.18  mm for tibia. The results demonstrate that this method outperforms previous segmentation methods in capturing anatomical shape variations, accounting for image quality differences and guiding accurate segmentation.

Convolutional neural networks (CNNs) show potential for computer-aided diagnosis (CADx) by learning features directly from the image data instead of using analytically extracted features. However, CNNs are difficult to train from scratch for medical images due to small sample sizes and variations in tumor presentations. Instead, transfer learning can be used to extract tumor information from medical images via CNNs originally pretrained for nonmedical tasks, alleviating the need for large datasets. Our database includes 219 breast lesions (607 full-field digital mammographic images). We compared support vector machine classifiers based on the CNN-extracted image features and our prior computer-extracted tumor features in the task of distinguishing between benign and malignant breast lesions. Five-fold cross validation (by lesion) was conducted with the area under the receiver operating characteristic (ROC) curve as the performance metric. Results show that classifiers based on CNN-extracted features (with transfer learning) perform comparably to those using analytically extracted features [area under the ROC curve (AUC)=0.81]. Further, the performance of ensemble classifiers based on both types was significantly better than that of either classifier type alone (AUC=0.86 versus 0.81, p=0.022). We conclude that transfer learning can improve current CADx methods while also providing standalone classifiers without large datasets, facilitating machine-learning methods in radiomics and precision medicine.

Cochlear implants (CIs) are neural prostheses that restore hearing by stimulating auditory nerve pathways within the cochlea using an implanted electrode array. Research has shown when multiple electrodes stimulate the same nerve pathways, competing stimulation occurs and hearing outcomes decline. Recent clinical studies have indicated that hearing outcomes can be significantly improved by using an image-guided active electrode set selection technique we have designed, in which electrodes that cause competing stimulation are identified and deactivated. In tests done to date, an expert is needed to perform the electrode selection step with the assistance of a method to visualize the spatial relationship between electrodes and neural sites determined using image analysis techniques. We propose to automate the electrode selection step by optimizing a cost function that captures the heuristics used by the expert. Further, we propose an approach to estimate the values of parameters used in the cost function using an existing database of expert electrode selections. We test this method with different electrode array models from three manufacturers. Our automatic approach generates acceptable active electrode sets in 98.3% of the subjects tested. This approach represents a crucial step toward clinical translation of our image-guided CI programming system.

Spectral computed tomography (SCT) generates better image quality than conventional computed tomography (CT). It has overcome several limitations for imaging atherosclerotic plaque. However, the literature evaluating the performance of SCT based on objective image assessment is very limited for the task of discriminating plaques. We developed a numerical-observer method and used it to assess performance on discrimination vulnerable-plaque features and compared the performance among multienergy CT (MECT), dual-energy CT (DECT), and conventional CT methods. Our numerical observer was designed to incorporate all spectral information and comprised two-processing stages. First, each energy-window domain was preprocessed by a set of localized channelized Hotelling observers (CHO). In this step, the spectral image in each energy bin was decorrelated using localized prewhitening and matched filtering with a set of Laguerre–Gaussian channel functions. Second, the series of the intermediate scores computed from all the CHOs were integrated by a Hotelling observer with an additional prewhitening and matched filter. The overall signal-to-noise ratio (SNR) and the area under the receiver operating characteristic curve (AUC) were obtained, yielding an overall discrimination performance metric. The performance of our new observer was evaluated for the particular binary classification task of differentiating between alternative plaque characterizations in carotid arteries. A clinically realistic model of signal variability was also included in our simulation of the discrimination tasks. The inclusion of signal variation is a key to applying the proposed observer method to spectral CT data. Hence, the task-based approaches based on the signal-known-exactly/background-known-exactly (SKE/BKE) framework and the clinical-relevant signal-known-statistically/background-known-exactly (SKS/BKE) framework were applied for analytical computation of figures of merit (FOM). Simulated data of a carotid-atherosclerosis patient were used to validate our methods. We used an extended cardiac-torso anthropomorphic digital phantom and three simulated plaque types (i.e., calcified plaque, fatty-mixed plaque, and iodine-mixed blood). The images were reconstructed using a standard filtered backprojection (FBP) algorithm for all the acquisition methods and were applied to perform two different discrimination tasks of: (1) calcified plaque versus fatty-mixed plaque and (2) calcified plaque versus iodine-mixed blood. MECT outperformed DECT and conventional CT systems for all cases of the SKE/BKE and SKS/BKE tasks (all p<0.01). On average of signal variability, MECT yielded the SNR improvements over other acquisition methods in the range of 46.8% to 65.3% (all p<0.01) for FBP-Ramp images and 53.2% to 67.7% (all p<0.01) for FBP-Hanning images for both identification tasks. This proposed numerical observer combined with our signal variability framework is promising for assessing material characterization obtained through the additional energy-dependent attenuation information of SCT. These methods can be further extended to other clinical tasks such as kidney or urinary stone identification applications.

Chronic liver disease is a worldwide health problem, and hepatic fibrosis (HF) is one of the hallmarks of the disease. The current reference standard for diagnosing HF is biopsy followed by pathologist examination; however, this is limited by sampling error and carries a risk of complications. Pathology diagnosis of HF is based on textural change in the liver as a lobular collagen network that develops within portal triads. The scale of collagen lobules is characteristically in the order of 1 to 5 mm, which approximates the resolution limit of in vivo gadolinium-enhanced magnetic resonance imaging in the delayed phase. We use MRI of formalin-fixed human ex vivo liver samples as phantoms that mimic the textural contrast of in vivo Gd-MRI. We have developed a local texture analysis that is applied to phantom images, and the results are used to train model observers to detect HF. The performance of the observer is assessed with the area-under-the-receiver–operator-characteristic curve (AUROC) as the figure-of-merit. To optimize the MRI pulse sequence, phantoms were scanned with multiple times at a range of flip angles. The flip angle that was associated with the highest AUROC was chosen as optimal for the task of detecting HF.

The use of a channelization mechanism on model observers not only makes mimicking human visual behavior possible, but also reduces the amount of image data needed to estimate the model observer parameters. The channelized Hotelling observer (CHO) and channelized scanning linear observer (CSLO) have recently been used to assess CT image quality for detection tasks and combined detection/estimation tasks, respectively. Although the use of channels substantially reduces the amount of data required to compute image quality, the number of scans required for CT imaging is still not practical for routine use. It is our desire to further reduce the number of scans required to make CHO or CSLO an image quality tool for routine and frequent system validations and evaluations. This work explores different data-reduction schemes and designs an approach that requires only a few CT scans. Three different kinds of approaches are included in this study: a conventional CHO/CSLO technique with a large sample size, a conventional CHO/CSLO technique with fewer samples, and an approach that we will show requires fewer samples to mimic conventional performance with a large sample size. The mean value and standard deviation of areas under ROC/EROC curve were estimated using the well-validated shuffle approach. The results indicate that an 80% data reduction can be achieved without loss of accuracy. This substantial data reduction is a step toward a practical tool for routine-task-based QA/QC CT system assessment.

This study aims to characterize the effect of background tissue density and heterogeneity on the detection of irregular masses in breast tomosynthesis, while demonstrating the capability of the sophisticated tools that can be used in the design, implementation, and performance analysis of virtual clinical trials (VCTs). Twenty breast phantoms from the extended cardiac-torso (XCAT) family, generated based on dedicated breast computed tomography of human subjects, were used to extract a total of 2173 volumes of interest (VOIs) from simulated tomosynthesis images. Five different lesions, modeled after human subject tomosynthesis images, were embedded in the breasts and combined with the lesion absent condition yielded a total of 6×2173 VOIs. Effects of background tissue density and heterogeneity on the detection of the lesions were studied by implementing a composite hypothesis signal detection paradigm with location known exactly, lesion known exactly or statistically, and background known statistically. Using the area under the receiver operating characteristic curve, detection performance deteriorated as density was increased, yielding findings consistent with clinical studies. A human observer study was performed on a subset of the simulated tomosynthesis images, confirming the detection performance trends with respect to density and serving as a validation of the implemented detector. Performance of the implemented detector varied substantially across the 20 breasts. Furthermore, background tissue density and heterogeneity affected the log-likelihood ratio test statistic differently under lesion absent and lesion present conditions. Therefore, considering background tissue variability in tissue models can change the outcomes of a VCT and is hence of crucial importance. The XCAT breast phantoms have the potential to address this concern by offering realistic modeling of background tissue variability based on a wide range of human subjects, comprising various breast shapes, sizes, and densities.

To address the error introduced by computed tomography (CT) scanners when assessing volume and unidimensional measurement of solid tumors, we scanned a precision manufactured pocket phantom simultaneously with patients enrolled in a lung cancer clinical trial. Dedicated software quantified bias and random error in the X,Y, and Z dimensions of a Teflon sphere and also quantified response evaluation criteria in solid tumors and volume measurements using both constant and adaptive thresholding. We found that underestimation bias was essentially the same for X,Y, and Z dimensions using constant thresholding and had similar values for adaptive thresholding. The random error of these length measurements as measured by the standard deviation and coefficient of variation was 0.10 mm (0.65), 0.11 mm (0.71), and 0.59 mm (3.75) for constant thresholding and 0.08 mm (0.51), 0.09 mm (0.56), and 0.58 mm (3.68) for adaptive thresholding, respectively. For random error, however, Z lengths had at least a fivefold higher standard deviation and coefficient of variation than X and Y. Observed Z-dimension error was especially high for some 8 and 16 slice CT models. Error in CT image formation, in particular, for models with low numbers of detector rows, may be large enough to be misinterpreted as representing either treatment response or disease progression.

The purpose of this study was to compare computed tomography (CT) low-contrast detectability from human readers with observer model-based surrogates of image quality. A phantom with a range of low-contrast signals (five contrasts, three sizes) was imaged on a state-of-the-art CT scanner (Siemens’ force). Images were reconstructed using filtered back projection and advanced modeled iterative reconstruction and were assessed by 11 readers using a two alternative forced choice method. Concurrently, contrast-to-noise ratio (CNR), area-weighted CNR (CNRA), and observer model-based metrics were estimated, including nonprewhitening (NPW) matched filter, NPW with eye filter (NPWE), NPW with internal noise, NPW with an eye filter and internal noise (NPWEi), channelized Hotelling observer (CHO), and CHO with internal noise (CHOi). The correlation coefficients (Pearson and Spearman), linear discriminator error, E, and magnitude of confidence intervals, |CI_95%|, were used to determine correlation, proper characterization of the reconstruction algorithms, and model precision, respectively. Pearson (Spearman) correlation was 0.36 (0.33), 0.83 (0.84), 0.84 (0.86), 0.86 (0.88), 0.86 (0.91), 0.88 (0.90), 0.85 (0.89), and 0.87 (0.84), E was 0.25, 0.15, 0.2, 0.25, 0.3, 0.25, 0.4, and 0.45, and |CI_95%| was 2.84×10⁻³, 5.29×10⁻³, 4.91×10⁻³, 4.55×10⁻³, 2.16×10⁻³, 1.24×10⁻³, 4.58×10⁻², and 7.95×10⁻² for CNR, CNRA, NPW, NPWE, NPWi, NPWEi, CHO, and CHOi, respectively.

The aim of this study was to evaluate the signal intensity and signal attenuation of swept source optical coherence tomography (SS-OCT) for dental caries in relation to the variation of mineral density. SS-OCT observation was performed on the enamel and dentin artificial demineralization and on natural caries. The artificial caries model on enamel and dentin surfaces was created using Streptococcus mutans biofilms incubated in an oral biofilm reactor. The lesions were centrally cross sectioned and SS-OCT scans were obtained in two directions to construct a three-dimensional data set, from the lesion surface (sagittal scan) and parallel to the lesion surface (horizontal scan). The integrated signal up to 200  μm in depth (IS200) and the attenuation coefficient (μ) of the enamel and dentin lesions were calculated from the SS-OCT signal in horizontal scans at five locations of lesion depth. The values were compared with the mineral density obtained from transverse microradiography. Both enamel and dentin demineralization showed significantly higher IS200 and μ than the sound tooth substrate from the sagittal scan. Enamel demineralization showed significantly higher IS200 than sound enamel, even with low levels of demineralization. In demineralized dentin, the μ from the horizontal scan consistently trended downward compared to the sound dentin.

Silicon-based nanoparticles are ideally suited for use as biomedical imaging agents due to their biocompatibility, biodegradability, and simple surface chemistry that facilitates drug loading and targeting. A method of hyperpolarizing silicon particles using dynamic nuclear polarization, which increases magnetic resonance imaging signals by several orders-of-magnitude through enhanced nuclear spin alignment, has recently been developed to allow silicon particles to function as contrast agents for in vivo magnetic resonance imaging. The enhanced spin polarization of silicon lasts significantly longer than other hyperpolarized agents (tens of minutes, whereas <1  min for other species at room temperature), allowing a wide range of potential applications. We report our recent characterizations of hyperpolarized silicon particles, with the ultimate goal of targeted, noninvasive, and nonradioactive molecular imaging of various cancer systems. A variety of particle sizes (20 nm to 2  μm) were found to have hyperpolarized relaxation times ranging from ∼10 to 50 min. The addition of various functional groups to the particle surface had no effect on the hyperpolarization buildup or decay rates and allowed in vivo imaging over long time scales. Additional in vivo studies examined a variety of particle administration routes in mice, including intraperitoneal injection, rectal enema, and oral gavage.

Active shape models (ASMs) have been widely used for extracting human anatomies in medical images given their capability for shape regularization of topology preservation. However, sensitivity to model initialization and local correspondence search often undermines their performances, especially around highly variable contexts in computed-tomography (CT) and magnetic resonance (MR) images. In this study, we propose an augmented ASM (AASM) by integrating the multiatlas label fusion (MALF) and level set (LS) techniques into the traditional ASM framework. Using AASM, landmark updates are optimized globally via a region-based LS evolution applied on the probability map generated from MALF. This augmentation effectively extends the searching range of correspondent landmarks while reducing sensitivity to the image contexts and improves the segmentation robustness. We propose the AASM framework as a two-dimensional segmentation technique targeting structures with one axis of regularity. We apply AASM approach to abdomen CT and spinal cord (SC) MR segmentation challenges. On 20 CT scans, the AASM segmentation of the whole abdominal wall enables the subcutaneous/visceral fat measurement, with high correlation to the measurement derived from manual segmentation. On 28 3T MR scans, AASM yields better performances than other state-of-the-art approaches in segmenting white/gray matter in SC.

The use of three-dimensional (3-D) echocardiography is limited by signal dropouts and narrow field of view. Data compounding is proposed as a solution to overcome these limitations by combining multiple 3-D recordings to form a wide field of view. The first step of the solution requires registration between the recordings both in the spatial and temporal dimension for dynamic organs such as the heart. Accurate registration between the individual echo recordings is crucial for the quality of compounded volumes. A temporal registration method based on a piecewise one-dimensional cubic B-spline in combination with multiscale iterative Farnebäck optic flow method for spatial registration was described. The temporal registration method was validated on in vivo data sets with annotated timing of mitral valve opening. The spatial registration method was validated using in vivo data and compared to registration with Procrustes analysis using manual contouring as a benchmark. The spatial accuracy was assessed in terms of mean of absolute distance and Hausdorff distance between the left ventricular contours. The results showed that the temporal registration accuracy is in the range of half the time resolution of the echo recordings and the achieved spatial accuracy of the proposed method is comparable to manual registration.