Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “image segmentation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,423 records · Page 79

SU-E-T-148: Efficient Verification Method for Modulated Electron Radiotherapy Treatment Plans

Purpose: For shallow tumors, modulated electron radiotherapy (MERT) promises a reduction of dose to distal organs at risk. At our institution a framework was developed in order to create treatment plans for MERT employing both forward and inverse optimization. In this work, an efficient quality assurance (QA) process is established. Methods: Treatment plans for three different tumor sites were created using an inverse optimization. These plans consist of 6–12 segments and energies between 6 and 18 MeV. An already established QA process for photon IMRT plans is now extended to additionally handle MERT plans. First, the dose distributions are calculated in a homogenous water phantom. For this task a dedicated Monte Carlo (MC) framework for MERT is used. Second, the segments are applied on a stand-alone amorphous silicon electronic portal imaging device (EPID) using a source-to-surface distance of 70 cm. This device was calibrated for electron beams in a previous work. An in-house developed analysis software, is then utilized for comparisons and evaluation of the measured and calculated dose distributions. Results: For all three plans the calculated dose distributions agree well with the measured ones. Using a 2D gamma comparison (2% of dose max/2 mm and 10% dose threshold) passing rates >98% are achieved.The dose calculation for each plan on the water phantom, using voxels of 0.2×0.2×0.2 cm{sup 3}, takes at maximum 30 min on a single core Pentium 2.66 GHz system with 6 GB RAM, to reach a statistical uncertainty of 2% (1 std. dev.). Conclusion: An already established QA procedure for IMRT photon plans was applied for MERT. The dedicated MC framework and the use of EPID measurements allow an efficient QA procedure in a clinical environment. This work was supported by Varian Medical Systems.

60 APPLIED LIFE SCIENCES↗

WE-H-207A-05: Spatial Co-Localization of F-18 NaF Vs. F-18 FDG Defined Disease Volumes

Purpose: Both [F-18]NaF and [F-18]FDG show promise for quantitative PET/CT assessment in metastatic prostate cancer to bone. Broad agreement between the tracers has been shown but voxel-wise correspondence has not been explored in depth. This study evaluates the spatial co-localization of [F-18]NaF PET and [F-18]FDG PET in bone lesions. Methods: Seventy-three lesion contours were identified in six patients receiving dynamic NaF PET/CT and FDG PET/CT scans two hours apart using identical fields-of-view. Tracer uptake (SUV) reflecting 60 minutes post-injection was modeled from kinetic parameters. Lesions were segmented by a physician separately on NaF PET and FDG PET. PET images were rigidly aligned using skeletal references on CT images. Lesion size, degree of overlap, voxel-wise tracer uptake values (SUV), and CT density distributions were compared using Dice coefficient, Positive Predictive Value (PPV), and Spearman rank correlation tests. Results: Across all patients, 42 lesions were identified on NaF PET (median 1.4 cm{sup 3}, range <1–204 cm{sup 3}) compared to 31 using FDG PET (median 1.8 cm{sup 3}, range <1–244 cm{sup 3}). Spatial cooccurrence was found in 25 lesion pairs. Lesions on NaF PET had PPV of 0.91 and on FDG a PPV of 0.65. Overall, NaF-defined lesions were 47% (±24%) larger by volume with moderate overlap to FDG, resulting in mean Dice coefficient of 34% (±22%). In areas of overlap, voxel-wise correlation of NaF and FDG SUV was moderate (ρ=0.56). Expanding to regions of non-spatial overlap, voxels contained in FDG-only contours were almost exclusively low HU (median 118), compared to dense regions of NaF-only voxels (median 250). In sclerotic sub-volumes (HU > 300) NaF-defined contours encompassed 83% of total FDG volume. Conclusion: Moderate voxel-wise correlation of FDG and NaF PET/CT uptake was observed. Spatial discrepancies in FDG and NaF PET/CT imaging of boney metastases could be influenced by poor sensitivity of FDG PET/CT in sclerotic regions. Funded by Prostate Cancer Foundation.

60 APPLIED LIFE SCIENCES↗

An Overview of the Thematic Mapper Geometric Correction System

The processing concepts which form the basis of the Thematic Mapper (TM) Geometric Correction System are examined. The principle flight and ground segment subsystems are discussed. Correction data is generated and compared to TM image data. Geometric accuracy is defined and investigated along with a general system overview.

Beyer, E. P.↗

SU-F-J-171: Robust Atlas Based Segmentation of the Prostate and Peripheral Zone Regions On MRI Utilizing Multiple MRI System Vendors

Purpose: Automatically generated prostate MRI contours can be used to aid in image registration with CT or ultrasound and to reduce the burden of contouring for radiation treatment planning. In addition, prostate and zonal contours can assist to automate quantitative imaging features extraction and the analyses of longitudinal MRI studies. These potential gains are limited if the solutions are not compatible across different MRI vendors. The goal of this study is to characterize an atlas based automatic segmentation procedure of the prostate collected on MRI systems from multiple vendors. Methods: The prostate and peripheral zone (PZ) were manually contoured by an expert radiation oncologist on T2-weighted scans acquired on both GE (n=31) and Siemens (n=33) 3T MRI systems. A leave-one-out approach was utilized where the target subject is removed from the atlas before the segmentation algorithm is initiated. The atlas-segmentation method finds the best nine matched atlas subjects and then performs a normalized intensity-based free-form deformable registration of these subjects to the target subject. These nine contours are then merged into a single contour using Simultaneous Truth and Performance Level Estimation (STAPLE). Contour comparisons were made using Dice similarity coefficients (DSC) and Hausdorff distances. Results: Using the T2 FatSat (FS) GE datasets the atlas generated contours resulted in an average DSC of 0.83±0.06 for prostate, 0.57±0.12 for PZ and 0.75±0.09 for CG. Similar results were found when using the Siemens data with a DSC of 0.79±0.14 for prostate, 0.54±0.16 and 0.70±0.9. Contrast between prostate and surrounding anatomy and between the PZ and CG contours for both vendors demonstrated superior contrast separation; significance was found for all comparisons p-value < 0.0001. Conclusion: Atlas-based segmentation yielded promising results for all contours compared to expertly defined contours in both Siemens and GE 3T systems providing fast and automatic segmentation of the prostate. Funding Support, Disclosures, and Conflict of Interest: AS Nelson is a partial owner of MIM Software, Inc. AS Nelson, and A Swallen are current employees at MIM Software, Inc.

60 APPLIED LIFE SCIENCES↗

Anatomy structure creation and editing using 3D implicit surfaces

Purpose: To accurately reconstruct, and interactively reshape 3D anatomy structures' surfaces using small numbers of 2D contours drawn in the most visually informative views of 3D imagery. The innovation of this program is that the number of 2D contours can be very much smaller than the number of transverse sections, even for anatomy structures spanning many sections. This program can edit 3D structures from prior segmentations, including those from autosegmentation programs. The reconstruction and surface editing works with any image modality. Methods: Structures are represented by variational implicit surfaces defined by weighted sums of radial basis functions (RBFs). Such surfaces are smooth, continuous, and closed and can be reconstructed with RBFs optimally located to efficiently capture shape in any combination of transverse (T), sagittal (S), and coronal (C) views. The accuracy of implicit surface reconstructions was measured by comparisons with the corresponding expert-contoured surfaces in 103 prostate cancer radiotherapy plans. Editing a pre-existing surface is done by overdrawing its profiles in image views spanning the affected part of the structure, deleting an appropriate set of prior RBFs, and merging the remainder with the new edit contour RBFs. Two methods were devised to identify RBFs to be deleted based only on the geometry of the initial surface and the locations of the new RBFs. Results: Expert-contoured surfaces were compared with implicit surfaces reconstructed from them over varying numbers and combinations of T/S/C planes. Studies revealed that surface-surface agreement increases monotonically with increasing RBF-sample density, and that the rate of increase declines over the same range. These trends were observed for all surface agreement metrics and for all the organs studied--prostate, bladder, and rectum. In addition, S and C contours may convey more shape information than T views for CT studies in which the axial slice thickness is greater than the pixel size. Surface editing accuracy likewise improves with larger sampling densities, and the rate of improvement similarly declines over the same conditions. Conclusions: Implicit surfaces based on RBFs are accurate representations of anatomic structures and can be interactively generated or modified to correct segmentation errors. The number of input contours is typically smaller than the number of T contours spanned by the structure.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Laser Guide Star for Large Segmented-aperture Space Telescopes. I. Implications for Terrestrial Exoplanet Detection and Observatory Stability

Precision wavefront control on future segmented-aperture space telescopes presents significant challenges, particularly in the context of high-contrast exoplanet direct imaging. We present a new wavefront control architecture that translates the ground-based artificial guide star concept to space with a laser source on board a second spacecraft, formation flying within the telescope’s field of view. We describe the motivating problem of mirror segment motion and develop wavefront sensing requirements as a function of guide star magnitude and segment motion power spectrum. Several sample cases with different values for transmitter power, pointing jitter, and wavelength are presented to illustrate the advantages and challenges of having a non-stellar-magnitude noise limited wavefront sensor for space telescopes. These notional designs allow increased control authority, potentially relaxing spacecraft stability requirements by two orders of magnitude and increasing terrestrial exoplanet discovery space by allowing high-contrast observations of stars of arbitrary brightness.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluation of magnetic shear in off-disk center active regions

The changes that projection effects produce in the evaluation of magnetic shear in off-disk center active regions by comparing angular shear calculated in image plane and heliographic coordinates are analyzed, and the procedure for properly evaluating magnetic shear by transforming the observed vector magnetic field into the heliographic system is described. This procedure is then used to evaluate magnetic shear along the magnetic neutral line in an active region that was observed on April 24, 1984 at a longitude offset of -45 deg. In particular, the number of 'critically sheared' pixels along an east-west directed segment of the neutral line in the leader sunspot group changes from 16 in the image plane magnetogram to 14 in the heliographic magnetogram. The critical shear as calculated in the image plane served as a good predictor for the location of flaring activity since the flare ribbons of the great flare of April 24 bracketed the inversion line where the critical shear was located. These results indicate that for this particular region, projection effects did not significantly affect the evaluation of critical shear.

Venkatakrishnan, P.↗

Method for implementation of recursive hierarchical segmentation on parallel computers

A method, computer readable storage, and apparatus for implementing a recursive hierarchical segmentation algorithm on a parallel computing platform. The method includes setting a bottom level of recursion that defines where a recursive division of an image into sections stops dividing, and setting an intermediate level of recursion where the recursive division changes from a parallel implementation into a serial implementation. The segmentation algorithm is implemented according to the set levels. The method can also include setting a convergence check level of recursion with which the first level of recursion communicates with when performing a convergence check.

Tilton, James C.↗

Estimating nonrigid motion from inconsistent intensity with robust shape features

Purpose: To develop a nonrigid motion estimation method that is robust to heterogeneous intensity inconsistencies amongst the image pairs or image sequence. Methods: Intensity and contrast variations, as in dynamic contrast enhanced magnetic resonance imaging, present a considerable challenge to registration methods based on general discrepancy metrics. In this study, the authors propose and validate a novel method that is robust to such variations by utilizing shape features. The geometry of interest (GOI) is represented with a flexible zero level set, segmented via well-behaved regularized optimization. The optimization energy drives the zero level set to high image gradient regions, and regularizes it with area and curvature priors. The resulting shape exhibits high consistency even in the presence of intensity or contrast variations. Subsequently, a multiscale nonrigid registration is performed to seek a regular deformation field that minimizes shape discrepancy in the vicinity of GOIs. Results: To establish the working principle, realistic 2D and 3D images were subject to simulated nonrigid motion and synthetic intensity variations, so as to enable quantitative evaluation of registration performance. The proposed method was benchmarked against three alternative registration approaches, specifically, optical flow, B-spline based mutual information, and multimodality demons. When intensity consistency was satisfied, all methods had comparable registration accuracy for the GOIs. When intensities among registration pairs were inconsistent, however, the proposed method yielded pronounced improvement in registration accuracy, with an approximate fivefold reduction in mean absolute error (MAE = 2.25 mm, SD = 0.98 mm), compared to optical flow (MAE = 9.23 mm, SD = 5.36 mm), B-spline based mutual information (MAE = 9.57 mm, SD = 8.74 mm) and mutimodality demons (MAE = 10.07 mm, SD = 4.03 mm). Applying the proposed method on a real MR image sequence also provided qualitatively appealing results, demonstrating good feasibility and applicability of the proposed method. Conclusions: The authors have developed a novel method to estimate the nonrigid motion of GOIs in the presence of spatial intensity and contrast variations, taking advantage of robust shape features. Quantitative analysis and qualitative evaluation demonstrated good promise of the proposed method. Further clinical assessment and validation is being performed.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Oriented Connectivity-Based Method for Segmenting Solar Loops

A method based on oriented connectivity that can automatically segnient arc-like structures (solar loops) from intensity images of the Sun's corona is introduced. The method is a constructive approach that uses model-guided processing to enable extraction of credible loop structures. Since the solar loops are vestiges of the solar magnetic field, the model-guided processing exploits external estimates of this field s local orientations that are derived from a physical magnetic field model. Empirical studies of the method s effectiveness are also presented. The Oriented Connectivity- Based Method is the first automatic method for the segmentation of solar loops.

Lee, J. K.↗

SU-E-J-209: Verification of 3D Surface Registration Between Stereograms and CT Images

Purpose: Stereography can provide a visualization of the skin surface for radiation therapy patients. The aim of this study was to verify the registration algorithm in a commercial image analysis software, 3dMDVultus, for the fusion of stereograms and CT images. Methods: CT and stereographic scans were acquired of a head phantom and a deformable phantom. CT images were imported in 3dMDVultus and the surface contours were generated by threshold segmentation. Stereograms were reconstructed in 3dMDVultus. The resulting surfaces were registered with Vultus algorithm and then exported to in-house registration software and compared with four algorithms: rigid, affine, non-rigid iterative closest point (ICP) and b-spline algorithm. RMS (root-mean-square residuals of the surface point distances) error between the registered CT and stereogram surfaces was calculated and analyzed. Results: For the head phantom, the maximum RMS error between registered CT surfaces to stereogram was 6.6 mm for Vultus algorithm, whereas the mean RMS error was 0.7 mm. For the deformable phantom, the maximum RMS error was 16.2 mm for Vultus algorithm, whereas the mean RMS error was 4.4 mm. Non-rigid ICP demonstrated the best registration accuracy, as the mean of RMS errors were both within 1 mm. Conclusion: The accuracy of registration algorithm in 3dMDVultus was verified and exceeded RMS of 2 mm for deformable cases. Non-rigid ICP and b-spline algorithms improve the registration accuracy for both phantoms, especially in deformable one. For those patients whose body habitus deforms during radiation therapy, more advanced nonrigid algorithms need to be used.

60 APPLIED LIFE SCIENCES↗

Pulmonary lobe segmentation based on ridge surface sampling and shape model fitting

Purpose: Performing lobe-based quantitative analysis of the lung in computed tomography (CT) scans can assist in efforts to better characterize complex diseases such as chronic obstructive pulmonary disease (COPD). While airways and vessels can help to indicate the location of lobe boundaries, segmentations of these structures are not always available, so methods to define the lobes in the absence of these structures are desirable. Methods: The authors present a fully automatic lung lobe segmentation algorithm that is effective in volumetric inspiratory and expiratory computed tomography (CT) datasets. The authors rely on ridge surface image features indicating fissure locations and a novel approach to modeling shape variation in the surfaces defining the lobe boundaries. The authors employ a particle system that efficiently samples ridge surfaces in the image domain and provides a set of candidate fissure locations based on the Hessian matrix. Following this, lobe boundary shape models generated from principal component analysis (PCA) are fit to the particles data to discriminate between fissure and nonfissure candidates. The resulting set of particle points are used to fit thin plate spline (TPS) interpolating surfaces to form the final boundaries between the lung lobes. Results: The authors tested algorithm performance on 50 inspiratory and 50 expiratory CT scans taken from the COPDGene study. Results indicate that the authors' algorithm performs comparably to pulmonologist-generated lung lobe segmentations and can produce good results in cases with accessory fissures, incomplete fissures, advanced emphysema, and low dose acquisition protocols. Dice scores indicate that only 29 out of 500 (5.85%) lobes showed Dice scores lower than 0.9. Two different approaches for evaluating lobe boundary surface discrepancies were applied and indicate that algorithm boundary identification is most accurate in the vicinity of fissures detectable on CT. Conclusions: The proposed algorithm is effective for lung lobe segmentation in absence of auxiliary structures such as vessels and airways. The most challenging cases are those with mostly incomplete, absent, or near-absent fissures and in cases with poorly revealed fissures due to high image noise. However, the authors observe good performance even in the majority of these cases.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Three-dimensional super line-localization in low signal-to-noise microscope images via prior-apprised unsupervised learning (PAUL)

Biological processes such as processive enzyme turnover and intracellular cargo tracking involve the dynamic motion of a small "article" along a curvilinear biopolymer track. To understand these processes that occur across multiple length and time scales, one must acquire both the trajectory of the particle and the position of the track along which it moves, possibly by combining high-resolution single-particle tracking with conventional microscopy. Yet, usually there is a significant resolution mismatch between these modalities: while the tracked particle is localized with a precision of 10 nm, the image of the surroundings is limited by optical difraction, with 200 nm lateral and 500 nm axial resolutions. Compared to the particle's trajectory, the surrounding curvilinear structure appears as a blurred and noisy image. This disparity in the spatial resolutions of the particle trajectory and the surrounding curvilinear structure image makes data reconstruction, as well as interpretation, particularly challenging. Analysis is further complicated when the curvilinear structures are oriented arbitrarily in 3D space. Here, we present a prior-apprised unsupervised learning (PAUL) approach to extract information from 3D images where the underlying features resemble a curved line such as a filament or microtubule. This three-stage framework starts with a Hessian-based feature enhancement, which is followed by feature registration, where local line segments are detected on repetitively sampled subimage tiles. In the final stage, statistical learning, segments are clustered based on their geometric relationships. Principal curves are then approximated from each segment group via statistical tools including principal component analysis, bootstrap and kernel transformation. This procedure is characterized on simulated images, where sub-voxel medium deviations from true curves have been achieved. The 3D PAUL approach has also been implemented for successful line localization in experimental 3D images of gold nanowires obtained using a multifocal microscope. Lastly, this work not only bridges the resolution gap between two microscopy modalities, but also allows us to conduct 3D super line-localization imaging experiments, without using super-resolution techniques.

42 ENGINEERING↗

SU-F-J-158: Respiratory Motion Resolved, Self-Gated 4D-MRI Using Rotating Cartesian K-Space Sampling

Purpose: Dynamic MRI has been used to quantify respiratory motion of abdominal organs in radiation treatment planning. Many existing 4D-MRI methods based on 2D acquisitions suffer from limited slice resolution and additional stitching artifacts when evaluated in 3D{sup 1}. To address these issues, we developed a 4D-MRI (3D dynamic) technique with true 3D k-space encoding and respiratory motion self-gating. Methods: The 3D k-space was acquired using a Rotating Cartesian K-space (ROCK) pattern, where the Cartesian grid was reordered in a quasi-spiral fashion with each spiral arm rotated using golden angle{sup 2}. Each quasi-spiral arm started with the k-space center-line, which were used as self-gating{sup 3} signal for respiratory motion estimation. The acquired k-space data was then binned into 8 respiratory phases and the golden angle ensures a near-uniform k-space sampling in each phase. Finally, dynamic 3D images were reconstructed using the ESPIRiT technique{sup 4}. 4D-MRI was performed on 6 healthy volunteers, using the following parameters (bSSFP, Fat-Sat, TE/TR=2ms/4ms, matrix size=500×350×120, resolution=1×1×1.2mm, TA=5min, 8 respiratory phases). Supplemental 2D real-time images were acquired in 9 different planes. Dynamic locations of the diaphragm dome and left kidney were measured from both 4D and 2D images. The same protocol was also performed on a MRI-compatible motion phantom where the motion was programmed with different amplitude (10–30mm) and frequency (3–10/min). Results: High resolution 4D-MRI were obtained successfully in 5 minutes. Quantitative motion measurements from 4D-MRI agree with the ones from 2D CINE (<5% error). The 4D images are free of the stitching artifacts and their near-isotropic resolution facilitates 3D visualization and segmentation of abdominal organs such as the liver, kidney and pancreas. Conclusion: Our preliminary studies demonstrated a novel ROCK 4D-MRI technique with true 3D k-space encoding and respiratory motion self-gating. The technique leads to high-resolution and artifacts-free 4D images for improved abdominal organ motion studies. K.S acknowledges funding support from NIH R01CA188300.

60 APPLIED LIFE SCIENCES↗

Technical Note: Suppression of artifacts arising from simultaneous cone-beam imaging and RF transponder tracking in prostate radiotherapy

Purpose: Artifacts in treatment-room cone-beam reconstructions have been observed at the authors' center when cone-beam acquisition is simultaneous with radio frequency (RF) transponder tracking using the Calypso 4D system (Calypso Medical, Seattle, WA). These artifacts manifest as CT-number modulations and increased CT-noise. The authors present a method for the suppression of the artifacts. Methods: The authors propose a three-stage postprocessing technique that can be applied to image volumes previously reconstructed by a cone-beam system. The stages are (1) segmentation of voxels into air, soft-tissue, and bone; (2) application of a 2D spatial-filter in the axial plane to the soft-tissue voxels; and (3) normalization to remove streaking along the axial-direction. The algorithm was tested on patient data acquired with Synergy XVI cone-beam CT systems (Elekta, Crawley, United Kingdom). Results: The computational demands of the suggested correction are small, taking less than 15 s per cone-beam reconstruction on a desktop PC. For a moderate loss of spatial-resolution, the artifacts are strongly suppressed and low-contrast visibility is improved. Conclusions: The correction technique proposed is fast and effective in removing the artifacts caused by simultaneous cone-beam imaging and RF-transponder tracking.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

WE-AB-207A-11: Respiratory Motion Guided 4DCBCT - A Step Towards Controlling 4DCBCT Image Quality

Purpose: We have developed a method, called respiratory motion guided 4DCBCT (RMG-4DCBCT), in which the gantry speed and projection frequency are varied in response to the patient’s real-time respiratory signal to eliminate streaking artifacts and to suppress duplicate projections in 4DCBCT images. In 2015, we realized RMG-4DCBCT on an Elekta Synergy linear accelerator with a mechanical relay to suppress projections and a potentiometer to adjust the gantry speed in response to the patient’s real-time respiratory signal. The aim of this study was to analyse the image quality to determine what can and cannot be controlled. Methods: Using RMG-4DCBCT, we acquired 40 (RMG-4DCBCT-40) and 60 (RMG-4DCBCT-60) equally spaced projections per respiratory phase of the CIRS dynamic thorax phantom with breathing periods from 2s to 8s and two breathing traces from lung cancer patients. The contrast to noise ratio (CNR) and edge response width (ERW) were used to compare image quality between RMG-4DCBCT and conventional 4DCBCT. Results: Regardless of the breathing period, for RMG-4DCBCT, the CNR is approximately 7 and 9 with RMG-4DCBCT-40 and RMG-4DCBCT-60 respectively. Conventional 4DCBCT has a CNR dropping from 20 down to 6 as the breathing period drops from 2s to 8s. With RMG-4DCBCT, the ERW, in the direction of phantom motion, ranges from 2.1mm to 2.5mm as the breathing period drops from 2s to 8s which compares to a higher range of 2.0mm to 2.5mm with conventional 4DCBCT. Images with similar quality to conventional 4DCBCT can be acquired with RMG-4DCBCT-40 which has a 70% reduction in imaging dose. Conclusion: The image contrast can be controlled with RMG-4DCBCT regardless of the patients breathing rate. However, although the image sharpness is better with RMG-4DCBCT, image sharpness has a small dependence on the breathing period; the accuracy of registration and segmentation will therefore vary with the patient’s breathing period. This project was supported by a National Health and Medical Research Council (NHMRC) project grant 1034060 and Cancer Australia grant number 1084566.

60 APPLIED LIFE SCIENCES↗