Anatomical pulmonary magnetic resonance imaging segmentation for regional structure-function measurements of asthma

PURPOSE

Pulmonary magnetic-resonance-imaging (MRI) and x-ray computed-tomography have provided strong evidence of spatially and temporally persistent lung structure-function abnormalities in asthmatics. This has generated a shift in their understanding of lung disease and supports the use of imaging biomarkers as intermediate endpoints of asthma severity and control. In particular, pulmonary (1)H MRI can be used to provide quantitative lung structure-function measurements longitudinally and in response to treatment. However, to translate such biomarkers of asthma, robust methods are required to segment the lung from pulmonary (1)H MRI. Therefore, their objective was to develop a pulmonary (1)H MRI segmentation algorithm to provide regional measurements with the precision and speed required to support clinical studies.

METHODS

The authors developed a method to segment the left and right lung from (1)H MRI acquired in 20 asthmatics including five well-controlled and 15 severe poorly controlled participants who provided written informed consent to a study protocol approved by Health Canada. Same-day spirometry and plethysmography measurements of lung function and volume were acquired as well as (1)H MRI using a whole-body radiofrequency coil and fast spoiled gradient-recalled echo sequence at a fixed lung volume (functional residual capacity + 1 l). We incorporated the left-to-right lung volume proportion prior based on the Potts model and derived a volume-proportion preserved Potts model, which was approximated through convex relaxation and further represented by a dual volume-proportion preserved max-flow model. The max-flow model led to a linear problem with convex and linear equality constraints that implicitly encoded the proportion prior. To implement the algorithm, (1)H MRI was resampled into ∼3 × 3 × 3 mm(3) isotropic voxel space. Two observers placed seeds on each lung and on the background of 20 pulmonary (1)H MR images in a randomized dataset, on five occasions, five consecutive days in a row. Segmentation accuracy was evaluated using the Dice-similarity-coefficient (DSC) of the segmented thoracic cavity with comparison to five-rounds of manual segmentation by an expert observer. The authors also evaluated the root-mean-squared-error (RMSE) of the Euclidean distance between lung surfaces, the absolute, and percent volume error. Reproducibility was measured using the coefficient of variation (CoV) and intraclass correlation coefficient (ICC) for two observers who repeated segmentation measurements five-times.

RESULTS

For five well-controlled asthmatics, forced expiratory volume in 1 s (FEV1)/forced vital capacity (FVC) was 83% ± 7% and FEV1 was 86 ± 9%pred. For 15 severe, poorly controlled asthmatics, FEV1/FV C = 66% ± 17% and FEV1 = 72 ± 27%pred. The DSC for algorithm and manual segmentation was 91% ± 3%, 92% ± 2% and 91% ± 2% for the left, right, and whole lung, respectively. RMSE was 4.0 ± 1.0 mm for each of the left, right, and whole lung. The absolute (percent) volume errors were 0.1 l (∼6%) for each of right and left lung and ∼0.2 l (∼6%) for whole lung. Intra- and inter-CoV (ICC) were <0.5% (>0.91%) for DSC and <4.5% (>0.93%) for RMSE. While segmentation required 10 s including ∼6 s for user interaction, the smallest detectable difference was 0.24 l for algorithm measurements which was similar to manual measurements.

CONCLUSIONS

This lung segmentation approach provided the necessary and sufficient precision and accuracy required for research and clinical studies.

Quantitative 3D and 2D Head Ultrasound to Determine Thresholds for Intervention In Preterm Neonates with Posthemorrhagic Ventricular Dilation

BACKGROUND: Preterm neonates with intraventricular hemorrhage (IVH) often acquire post hemorrhagic ventricle dilation (PHVD), which, when severe, can lead to neurological impairment. Cranial 2D ultrasound (US) images are used for the diagnosis and monitoring of PHVD; however, there is no consensus on the use of 2D US images to guide treatment. This can lead to delays in interventions, and the potential for brain injury. We have developed a 3D US system that has been shown to accurately detect changes ventricle volumes (VV).

OBJECTIVES: We investigate the utility of using 3D and 2D US measurements to determine thresholds for treatment of neonates with PHVD and to predict the need for further treatments.

DESIGN/METHODS: Neonates were imaged twice weekly in accordance to a protocol approved by the research ethics board. 3D US images were manually segmented to obtain VV. 2D measurements included ventricle index, anterior horn widths, third ventricle width, and largest thalamo-occipital distance. The rate of change for each measurement was calculated. Decisions to perform ventricular taps (VTs) to relieve intracranial pressure were made independently by neurosurgeons who were blinded to study images. Receiver operator curves (ROC) were generated using the sensitivity and specificity of the rates of change of sonographic parameters in predicting the need for V T. For each parameter optimal threshold for intervention was estimated by the area under ROC; and positive and negative predictive values (PPV, NPV) were calculated. Additionally, we investigated whether US measurements predicted the need for multiple interventions.

RESULTS: 23 neonates with PHVD were enrolled, 8 required interventions. The best predictor to determine initial intervention was the rate of change in VV when a threshold of &gt;2.04 cm3/day was used within the first three weeks of life (NPV and PPV of 1) and, this measurement was able to determine if then a patient would require further interventions when a threshold of -0.04 cm3/day was used looking at imaging time points after the first intervention (NPV and PPV of 1). 2D measurements were less sensitive and/or less specific (sensitivity of 88-57%, specificity of 100-79%, PPV of 0.88-0.57 and NPV of 0.93-0.79).

CONCLUSION: 3D US VV can predict the requirement for interven-tional ventricular tap in neonates with IVH, and can identify patients that have resolving PHVD following initial intervention, with higher sensitivity and specificity than 2D US measurements. These findings show promise for early classification of neonates using 3D US for prediction of interven-tional therapy, potentiallyaiding in timely management of these patients.

3D prostate histology reconstruction: an evaluation of image-based and fiducial-based algorithms

PURPOSE

Evaluation of in vivo prostate imaging modalities for determining the spatial distribution and aggressiveness of prostate cancer ideally requires accurate registration of images to an accepted reference standard, such as histopathological examination of radical prostatectomy specimens. Three-dimensional (3D) reconstruction of prostate histology facilitates these registration-based evaluations by reintroducing 3D spatial information lost during histology processing. Because the reconstruction accuracy may constrain the clinical questions that can be answered with these data, it is important to assess the tradeoffs between minimally disruptive methods based on intrinsic image information and potentially more robust methods based on extrinsic fiducial markers.

METHODS

Ex vivo magnetic resonance (MR) images and digitized whole-mount histology images from 12 radical prostatectomy specimens were used to evaluate four 3D histology reconstruction algorithms. 3D reconstructions were computed by registering each histology image to the corresponding ex vivo MR image using one of two similarity metrics (mutual information or fiducial registration error) and one of two search domains (affine transformations or a constrained subset thereof). The algorithms were evaluated for accuracy using the mean target registration error (TRE) computed from homologous intrinsic point landmarks (3-16 per histology section; 232 total) identified on histology and MR images, and for the sensitivity of TRE to rotational, translational, and scaling initialization errors.

RESULTS

The algorithms using fiducial registration error and mutual information had mean ± standard deviation TREs of 0.7 ± 0.4 and 1.2 ± 0.7 mm, respectively, and one algorithm using fiducial registration error and affine transforms had negligible sensitivities to initialization errors. The postoptimization values of the mutual information-based metric showed evidence of errors due to both the optimizer and the similarity metric, and variation of parameters of the mutual information-based metric did not improve its performance.

CONCLUSIONS

The extrinsic fiducial-based algorithm had lower mean TRE and lower sensitivity to initialization than the intrinsic intensity-based algorithm using mutual information. A model relating statistical power to registration error for certain imaging validation study designs estimated that a reconstruction algorithm with a mean TRE of 0.7 mm would require 27% fewer subjects than the method used to initialize the algorithms (mean TRE 1.3 ± 0.7 mm), suggesting the choice of reconstruction technique can have a substantial impact on the design of imaging validation studies, and on their overall cost.

Three-dimensional segmentation of three-dimensional ultrasound carotid atherosclerosis using sparse field level sets

PURPOSE

Three-dimensional ultrasound (3DUS) vessel wall volume (VWV) provides a 3D measurement of carotid artery wall remodeling and atherosclerotic plaque and is sensitive to temporal changes of carotid plaque burden. Unfortunately, although 3DUS VWV provides many advantages compared to measurements of arterial wall thickening or plaque alone, it is still not widely used in research or clinical practice because of the inordinate amount of time required to train observers and to generate 3DUS VWV measurements. In this regard, semiautomated methods for segmentation of the carotid media-adventitia boundary (MAB) and the lumen-intima boundary (LIB) would greatly improve the time to train observers and for them to generate 3DUS VWV measurements with high reproducibility.

METHODS

The authors describe a 3D algorithm based on a modified sparse field level set method for segmenting the MAB and LIB of the common carotid artery (CCA) from 3DUS images. To the authors' knowledge, the proposed algorithm is the first direct 3D segmentation method, which has been validated for segmenting both the carotid MAB and the LIB from 3DUS images for the purpose of computing VWV. Initialization of the algorithm requires the observer to choose anchor points on each boundary on a set of transverse slices with a user-specified interslice distance (ISD), in which larger ISD requires fewer user interactions than smaller ISD. To address the challenges of the MAB and LIB segmentations from 3DUS images, the authors integrated regional- and boundary-based image statistics, expert initializations, and anatomically motivated boundary separation into the segmentation. The MAB is segmented by incorporating local region-based image information, image gradients, and the anchor points provided by the observer. Moreover, a local smoothness term is utilized to maintain the smooth surface of the MAB. The LIB is segmented by constraining its evolution using the already segmented surface of the MAB, in addition to the global region-based information and the anchor points. The algorithm-generated surfaces were sliced and evaluated with respect to manual segmentations on a slice-by-slice basis using 21 3DUS images.

RESULTS

The authors used ISD of 1, 2, 3, 4, and 10 mm for algorithm initialization to generate segmentation results. The algorithm-generated accuracy and intraobserver variability results are comparable to the previous methods, but with fewer user interactions. For example, for the ISD of 3 mm, the algorithm yielded an average Dice coefficient of 94.4% ± 2.2% and 90.6% ± 5.0% for the MAB and LIB and the coefficient of variation of 6.8% for computing the VWV of the CCA, while requiring only 1.72 min (vs 8.3 min for manual segmentation) for a 3DUS image.

CONCLUSIONS

The proposed 3D semiautomated segmentation algorithm yielded high-accuracy and high-repeatability, while reducing the expert interaction required for initializing the algorithm than the previous 2D methods.

3D ultrasound system to investigate intraventricular hemorrhage in preterm neonates

Intraventricular hemorrhage (IVH) is a common disorder among preterm neonates that is routinely diagnosed and monitored by 2D cranial ultrasound (US). The cerebral ventricles of patients with IVH often have a period of ventricular dilation (ventriculomegaly). This initial increase in ventricle size can either spontaneously resolve, which often shows clinically as a period of stabilization in ventricle size and eventual decline back towards a more normal size, or progressive ventricular dilation that does not stabilize and which may require interventional therapy to reduce symptoms relating to increased intracranial pressure. To improve the characterization of ventricle dilation, we developed a 3D US imaging system that can be used with a conventional clinical US scanner to image the ventricular system of preterm neonates at risk of ventriculomegaly. A motorized transducer housing was designed specifically for hand-held use inside an incubator using a transducer commonly used for cranial 2D US scans. This system was validated using geometric phantoms, US/MRI compatible ventricle volume phantoms, and patient images to determine 3D reconstruction accuracy and inter- and intra-observer volume estimation variability. 3D US geometric reconstruction was found to be accurate with an error of <0.2%. Measured volumes of a US/MRI compatible ventricle-like phantom were within 5% of gold standard water displacement measurements. Intra-class correlation for the three observers was 0.97, showing very high agreement between observers. The coefficient of variation was between 1.8-6.3% for repeated segmentations of the same patient. The minimum detectable difference was calculated to be 0.63 cm(3) for a single observer. Results from ANOVA for three observers segmenting three patients of IVH grade II did not show any significant differences (p > 0.05) for the measured ventricle volumes between observers. This 3D US system can reliably produce 3D US images of the neonatal ventricular system. There is the potential to use this system to monitor the progression of ventriculomegaly over time in patients with IVH.

3D ultrasound system to investigate intraventricular hemorrhage in preterm neonates

Intraventricular hemorrhage (IVH) is a common disorder among preterm neonates that is routinely diagnosed and monitored by 2D cranial ultrasound (US). The cerebral ventricles of patients with IVH often have a period of ventricular dilation (ventriculomegaly). This initial increase in ventricle size can either spontaneously resolve, which often shows clinically as a period of stabilization in ventricle size and eventual decline back towards a more normal size, or progressive ventricular dilation that does not stabilize and which may require interventional therapy to reduce symptoms relating to increased intracranial pressure. To improve the characterization of ventricle dilation, we developed a 3D US imaging system that can be used with a conventional clinical US scanner to image the ventricular system of preterm neonates at risk of ventriculomegaly. A motorized transducer housing was designed specifically for hand-held use inside an incubator using a transducer commonly used for cranial 2D US scans. This system was validated using geometric phantoms, US/MRI compatible ventricle volume phantoms, and patient images to determine 3D reconstruction accuracy and inter- and intra-observer volume estimation variability. 3D US geometric reconstruction was found to be accurate with an error of <0.2%. Measured volumes of a US/MRI compatible ventricle-like phantom were within 5% of gold standard water displacement measurements. Intra-class correlation for the three observers was 0.97, showing very high agreement between observers. The coefficient of variation was between 1.8-6.3% for repeated segmentations of the same patient. The minimum detectable difference was calculated to be 0.63 cm(3) for a single observer. Results from ANOVA for three observers segmenting three patients of IVH grade II did not show any significant differences (p > 0.05) for the measured ventricle volumes between observers. This 3D US system can reliably produce 3D US images of the neonatal ventricular system. There is the potential to use this system to monitor the progression of ventriculomegaly over time in patients with IVH.

Feasibility of 3D ultrasound to evaluate upper extremity nerves

PURPOSE

This study investigates the performance of a 3 D Ultrasound (US) system in imaging elbow and wrist nerves.

MATERIALS AND METHODS

Twenty healthy volunteers with asymptomatic median, ulnar and radial nerves were prospectively investigated. Bilateral 3DUS scans of the elbows and wrists were acquired by using a commercially available US scanner (18 MHz, AplioXG, Toshiba) and stored as a 3 D volume by a dedicated software (CURE, Robarts Research Institute). Retrospectively, qualitative (image quality, atypical nerve location, findings potentially associated with compression neuropathy) and quantitative (cross-sectional area measurements) evaluations were performed.

RESULTS

In all 200 nerves 3DUS was feasible (100%). Image quality was insufficient in 13.5% (25 ulnar nerve elbow, 2 radial nerve) and sonomorphology was not assessable in those nerves. Measurement of cross sectional areas was feasible in all nerves (100%). Median cross-sectional area (range) were: median nerve elbow 7 mm2 (6-9), radial nerve 3 mm2 (1-4), ulnar nerve elbow 8 mm2 (5-11), median nerve wrist 8 mm2 (5-10), and ulnar nerve wrist 4 mm2 (2-6). No significant changes in nerve cross-sectional area along each nerve was found. Ulnar nerve subluxation was found in 2 nerves (6.7%). No anconeus epitrochlearis muscle or osteophytes were found.

CONCLUSION

3DUS is a feasible method for assessing nerves of the upper extremity and has been shown to provide a good overview of the median, ulnar and radial nerve at the elbow and wrist, but is limited for evaluation of the ulnar nerve in the cubital tunnel. This technique enables reliable measurements at different locations along the nerve.

An automated neural-fuzzy approach to malignant tumor localization in 2D ultrasonic images of the prostate

In this paper, a new neural-fuzzy approach is proposed for automated region segmentation in transrectal ultrasound images of the prostate. The goal of region segmentation is to identify suspicious regions in the prostate in order to provide decision support for the diagnosis of prostate cancer. The new automated region segmentation system uses expert knowledge as well as both textural and spatial features in the image to accomplish the segmentation. The textural information is extracted by two recurrent random pulsed neural networks trained by two sets of data (a suspicious tissues' data set and a normal tissues' data set). Spatial information is captured by the atlas-based reference approach and is represented as fuzzy membership functions. The textural and spatial features are synthesized by a fuzzy inference system, which provides a binary classification of the region to be evaluated.

A new lateral guidance device for stereotactic breast biopsy using an add-on unit to an upright mammography system

Stereotactic breast biopsy (SBB) is the gold standard for noninvasive breast cancer diagnosis. Current systems rely on one of two methods for needle insertion: a vertical-approach (from above the breast compression plate) or a lateral-approach (parallel to the compression plate). While the vertical-approach is more commonly used, it is not feasible in patients with thin breasts (less than 3 cm thickness after compression) or with superficial lesions. We present a novel design of lateral guidance device for SBB which addresses these limitations of the vertical-approach, and provides improvements over existing lateral guidance hardware. This device incorporates spherical linkages to allow two degrees of rotational freedom in the needle trajectory for increased targeting flexibility, as well as an adjustable rigid needle support to minimize needle deflection within the tissue. Needle placement error in SBB experiments is compared using both the new lateral guidance device and a commercial lateral guidance device in agar phantoms. The effect of elevation angle on needle placement accuracy using the new lateral guidance device is also assessed. Finally, a biopsy accuracy experiment is presented using a certified SBB phantom to compare the new design and the commercial lateral guidance device. In these experiments, SBB performed using the new lateral guidance device resulted in improved needle placement error and biopsy accuracy, while increasing targeting flexibility and maintaining procedural workflow.

Poster - Thurs Eve-02: 3 dimensional ultrasound-guided breast brachytherapy

Breast cancer is one of Canada's leading causes of death, taking the lives of approximately 5000 people annually. Breast-conserving tumour excision, or lumpectomy followed by radiation therapy is becoming an increasingly common treatment method for smaller tumours. High dose rate (HDR) brachytherapy is a precise form of radiation delivery following surgery involving the delivery of radiation dose through an HDR afterloader attached to catheters inserted into the breast. Currently, a CT scan of the patient is taken to properly reconstruct the tumour and guide catheter insertion. We propose to use our three-dimensional ultrasound (3DUS) scanner as the primary treatment planning device, eliminating the need for a CT scan. This would greatly increase patient comfort along with saving time and money. We have designed and constructed a 3DUS scanner specifically to be used in breast brachytherapy. It attaches to a Kuske breast application kit already used in the clinic. Software to view 3DUS images produced by the device is already being used, and needle guidance software is currently being developed. Laboratory tests on agar phantoms are set to begin shortly to evaluate the precision of the device and perform brachytherapy catheter insertion simulations. When the results of these tests are satisfactory, a full brachytherapy procedure will be performed in Quebec City using the 3DUS scanner.

Sci-Fri AM: YIS-05: A new guidance device for lateral-approach stereotactic breast biopsy

Stereotactic breast biopsy (SBB) is the gold standard for noninvasive breast cancer diagnosis. Current systems rely on one of two methods for needle insertion: a top-approach (from above the breast compression plate) or a lateral-approach (parallel to the compression plate). While the top-approach is more commonly used, it is not feasible in patients with thin breasts (less than 2.5 cm thickness after compression), or with superficial lesions. We present a novel design of lateral guidance support for SBB, which addresses these limitations of the top-approach, and provides improvements over existing lateral support hardware. This device incorporates spherical linkages to allow two degrees of rotational freedom in the needle trajectory for increased targeting flexibility, as well as an adjustable rigid needle support to minimize needle deflection within the tissue. Needle placement error in SBB experiments is compared using both the new lateral guidance device and a commercial lateral guidance device in agar phantoms. The effect of elevation angle on needle placement accuracy using the new lateral guidance device is also assessed. Finally, a biopsy accuracy experiment is presented using a certified SBB phantom to compare the new design and the commercial lateral guidance device. In these experiments, SBB performed using the new lateral guidance device results in improved needle placement error and biopsy accuracy, while increasing targeting flexibility and maintaining procedural workflow.

Stereotactic mammography imaging combined with 3D US imaging for image guided breast biopsy

Stereotactic X-ray mammography (SM) and ultrasound (US) guidance are both commonly used for breast biopsy. While SM provides three-dimensional (3D) targeting information and US provides real-time guidance, both have limitations. SM is a long and uncomfortable procedure and the US guided procedure is inherently two dimensional (2D), requiring a skilled physician for both safety and accuracy. The authors developed a 3D US-guided biopsy system to be integrated with, and to supplement SM imaging. Their goal is to be able to biopsy a larger percentage of suspicious masses using US, by clarifying ambiguous structures with SM imaging. Features from SM and US guided biopsy were combined, including breast stabilization, a confined needle trajectory, and dual modality imaging. The 3D US guided biopsy system uses a 7.5 MHz breast probe and is mounted on an upright SM machine for preprocedural imaging. Intraprocedural targeting and guidance was achieved with real-time 2D and near real-time 3D US imaging. Postbiopsy 3D US imaging allowed for confirmation that the needle was penetrating the target. The authors evaluated 3D US-guided biopsy accuracy of their system using test phantoms. To use mammographic imaging information, they registered the SM and 3D US coordinate systems. The 3D positions of targets identified in the SM images were determined with a target localization error (TLE) of 0.49 mm. The z component (x-ray tube to image) of the TLE dominated with a TLEz of 0.47 mm. The SM system was then registered to 3D US, with a fiducial registration error (FRE) and target registration error (TRE) of 0.82 and 0.92 mm, respectively. Analysis of the FRE and TRE components showed that these errors were dominated by inaccuracies in the z component with a FREz of 0.76 mm and a TREz of 0.85 mm. A stereotactic mammography and 3D US guided breast biopsy system should include breast compression for stability and safety and dual modality imaging for target localization. The system will provide preprocedural x-ray mammography information in the form of SM imaging along with real-time US imaging for needle guidance to a target. 3D US imaging will also be available for targeting, guidance, and biopsy verification immediately postbiopsy.

3D thoracoscopic ultrasound volume measurement validation in an ex vivo and in vivo porcine model of lung tumours

The purpose of this study was to validate the accuracy and reliability of volume measurements obtained using three-dimensional (3D) thoracoscopic ultrasound (US) imaging. Artificial "tumours" were created by injecting a liquid agar mixture into spherical moulds of known volume. Once solidified, the "tumours" were implanted into the lung tissue in both a porcine lung sample ex vivo and a surgical porcine model in vivo. 3D US images were created by mechanically rotating the thoracoscopic ultrasound probe about its long axis while the transducer was maintained in close contact with the tissue. Volume measurements were made by one observer using the ultrasound images and a manual-radial segmentation technique and these were compared with the known volumes of the agar. In vitro measurements had average accuracy and precision of 4.76% and 1.77%, respectively; in vivo measurements had average accuracy and precision of 8.18% and 1.75%, respectively. The 3D thoracoscopic ultrasound can be used to accurately and reproducibly measure "tumour" volumes both in vivo and ex vivo.