Assessing quality assurance of multi-leaf collimator using the structural similarity index

PURPOSE

This study aims to evaluate the viability of utilizing the Structural Similarity Index (SSI*) as an innovative imaging metric for quality assurance (QA) of the multi-leaf collimator (MLC). Additionally, we compared the results obtained through SSI* with those derived from a conventional Gamma index test for three types of Varian machines (Trilogy, Truebeam, and Edge) over a 12-week period of MLC QA in our clinic.

METHOD

To assess sensitivity to MLC positioning errors, we designed a 1 cm slit on the reference MLC, subsequently shifted by 0.5-5 mm on the target MLC. For evaluating sensitivity to output error, we irradiated five 25 cm × 25 cm open fields on the portal image with varying Monitor Units (MUs) of 96-100. We compared SSI* and Gamma index tests using three linear accelerator (LINAC) machines: Varian Trilogy, Truebeam, and Edge, with MLC leaf widths of 1, 0.5, and 0.25 mm. Weekly QA included VMAT and static field modes, with Picket fence test images acquired. Mechanical uncertainties related to the LINAC head, electronic portal imaging device (EPID), and MLC during gantry rotation and leaf motion were monitored.

RESULTS

The Gamma index test started detecting the MLC shift at a threshold of 4 mm, whereas the SSI* metric showed sensitivity to shifts as small as 2 mm. Moreover, the Gamma index test identified dose changes at 95MUs, indicating a 5% dose difference based on the distance to agreement (DTA)/dose difference (DD) criteria of 1 mm/3%. In contrast, the SSI* metric alerted to dose differences starting from 97MUs, corresponding to a 3% dose difference. The Gamma index test passed all measurements conducted on each machine. However, the SSI* metric rejected all measurements from the Edge and Trilogy machines and two from the Truebeam.

CONCLUSIONS

Our findings demonstrate that the SSI* exhibits greater sensitivity than the Gamma index test in detecting MLC positioning errors and dose changes between static and VMAT modes. The SSI* metric outperformed the Gamma index test regarding sensitivity across these parameters.

Enhancing liver tumor localization accuracy by prior-knowledge-guided motion modeling and a biomechanical model

BACKGROUND

Pre-treatment liver tumor localization remains a challenging task for radiation therapy, mostly due to the limited tumor contrast against normal liver tissues, and the respiration-induced liver tumor motion. Recently, we developed a biomechanical modeling-based, deformation-driven cone-beam CT estimation technique (Bio-CBCT), which achieved substantially improved accuracy on low-contrast liver tumor localization. However, the accuracy of Bio-CBCT is still affected by the limited tissue contrast around the caudal liver boundary, which reduces the accuracy of the boundary condition that is fed into the biomechanical modeling process. In this study, we developed a motion modeling and biomechanical modeling-guided CBCT estimation technique (MM-Bio-CBCT), to further improve the liver tumor localization accuracy by incorporating a motion model into the CBCT estimation process.

METHODS

MM-Bio-CBCT estimates new CBCT images through deforming a prior high-quality CT or CBCT volume. The deformation vector field (DVF) is solved by iteratively matching the digitally-reconstructed-radiographs (DRRs) of the deformed prior image to the acquired 2D cone-beam projections. Using the same solved DVF, the liver tumor volume contoured on the prior image can be transferred onto the new CBCT image for automatic tumor localization. To maximize the accuracy of the solved DVF, MM-Bio-CBCT employs two strategies for additional DVF optimization: (I) prior-knowledge-guided liver boundary motion modeling with motion patterns extracted from a prior 4D imaging set like 4D-CTs/4D-CBCTs, to improve the liver boundary DVF accuracy; and (II) finite-element-analysis-based biomechanical modeling of the liver volume to improve the intra-liver DVF accuracy. We evaluated the accuracy of MM-Bio-CBCT on both the digital extended-cardiac-torso (XCAT) phantom images and real liver patient images. The liver tumor localization accuracy of MM-Bio-CBCT was evaluated and compared with that of the purely intensity-driven 2D-3D deformation technique, the 2D-3D deformation technique with motion modeling, and the Bio-CBCT technique. Metrics including the DICE coefficient and the center-of-mass-error (COME) were assessed for quantitative evaluation.

RESULTS

Using limited-view 20 projections for CBCT estimation, the average (± SD) DICE coefficients between the estimated and the 'gold-standard' liver tumors of the XCAT study were 0.57±0.31, 0.78±0.26, 0.83±0.21, and 0.89±0.11 for 2D-3D deformation, 2D-3D deformation with motion modeling, Bio-CBCT and MM-Bio-CBCT techniques, respectively. Using 20 projections for estimation, the patient study yielded average DICE results of 0.63±0.21, 0.73±0.13 and 0.78±0.12, and 0.83±0.09, correspondingly. The MM-Bio-CBCT localized the liver tumor to an average COME of ~2 mm for both the XCAT and the liver patient studies.

CONCLUSIONS

Compared to Bio-CBCT, MM-Bio-CBCT further improves the accuracy of liver tumor localization. MM-Bio-CBCT can potentially be used towards pre-treatment liver tumor localization and intra-treatment liver tumor location verification to achieve substantial radiotherapy margin reduction.

4D liver tumor localization using cone-beam projections and a biomechanical model

PURPOSE

To improve the accuracy of liver tumor localization, this study tests a biomechanical modeling-guided liver cone-beam CT (CBCT) estimation (Bio-CBCT-est) technique, which generates new CBCTs by deforming a prior high-quality CT or CBCT image using deformation vector fields (DVFs). The DVFs can be used to propagate tumor contours from the prior image to new CBCTs for automatic 4D tumor localization.

METHODS/MATERIALS

To solve the DVFs, the Bio-CBCT-est technique employs an iterative scheme that alternates between intensity-driven 2D-3D deformation and biomechanical modeling-guided DVF regularization and optimization. The 2D-3D deformation step solves DVFs by matching digitally reconstructed radiographs of the 3D deformed prior image to 2D phase-sorted on-board projections according to imaging intensities. This step's accuracy is limited at low-contrast intra-liver regions without sufficient intensity variations. To boost the DVF accuracy in these regions, we use the intensity-driven DVFs solved at higher-contrast liver boundaries to fine-tune the intra-liver DVFs by finite element analysis-based biomechanical modeling. We evaluated Bio-CBCT-est's accuracy with seven liver cancer patient cases. For each patient, we simulated 4D cone-beam projections from 4D-CT images, and used these projections for Bio-CBCT-est based image estimations. After Bio-CBCT-est, the DVF-propagated liver tumor/cyst contours were quantitatively compared with the manual contours on the original 4D-CT 'reference' images, using the DICE similarity index, the center-of-mass-error (COME), the Hausdorff distance (HD) and the voxel-wise cross-correlation (CC) metrics. In addition to simulation, we also performed a preliminary study to qualitatively evaluate the Bio-CBCT-est technique via clinically acquired cone beam projections. A quantitative study using an in-house deformable liver phantom was also performed.

RESULTS

Using 20 projections for image estimation, the average (±s.d.) DICE index increased from 0.48 ± 0.13 (by 2D-3D deformation) to 0.77 ± 0.08 (by Bio-CBCT-est), the average COME decreased from 7.7 ± 1.5 mm to 2.2 ± 1.2 mm, the average HD decreased from 10.6 ± 2.2 mm to 5.9 ± 2.0 mm, and the average CC increased from -0.004 ± 0.216 to 0.422 ± 0.206. The tumor/cyst trajectory solved by Bio-CBCT-est matched well with that manually obtained from 4D-CT reference images.

CONCLUSIONS

Bio-CBCT-est substantially improves the accuracy of 4D liver tumor localization via cone-beam projections and a biomechanical model.

Lung surface deformation prediction from spirometry measurement and chest wall surface motion

PURPOSE

The authors have developed and evaluated a method to predict lung surface motion based on spirometry measurements, and chest and abdomen motion at selected locations.

METHODS

A patient-specific 3D triangular surface mesh of the lung region was obtained at the end expiratory phase by the threshold-based segmentation method. Lung flow volume changes were recorded with a spirometer for each patient. A total of 192 selected points at a regular spacing of 2 × 2 cm matrix points were used to detect chest wall motion over a total area of 32 × 24 cm covering the chest and abdomen surfaces. QR factorization with column pivoting was employed to remove redundant observations of the chest and abdominal areas. To create a statistical model between the lung surface and the corresponding surrogate signals, the authors developed a predictive model based on canonical ridge regression. Two unique weighting vectors were selected for each vertex on the lung surface; they were optimized during the training process using all other 4D-CT phases except for the test inspiration phase. These parameters were employed to predict the vertex locations of a testing data set.

RESULTS

The position of each lung surface mesh vertex was estimated from the motion at selected positions within the chest wall surface and from spirometry measurements in ten lung cancer patients. The average estimation of the 98th error percentile for the end inspiration phase was less than 1 mm (AP = 0.9 mm, RL = 0.6 mm, and SI = 0.8 mm). The vertices located at the lower region of the lung had a larger estimation error as compared with those within the upper region of the lung. The average landmark motion errors, derived from the biomechanical modeling using real surface deformation vector fields (SDVFs), and the predicted SDVFs were 3.0 and 3.1 mm, respectively.

CONCLUSIONS

Our newly developed predictive model provides a noninvasive approach to derive lung boundary conditions. The proposed system can be used with personalized biomechanical respiration modeling to derive lung tumor motion during radiation therapy from noninvasive measurements.

A Biomechanical Modeling Guided CBCT Estimation Technique

Two-dimensional-to-three-dimensional (2D-3D) deformation has emerged as a new technique to estimate cone-beam computed tomography (CBCT) images. The technique is based on deforming a prior high-quality 3D CT/CBCT image to form a new CBCT image, guided by limited-view 2D projections. The accuracy of this intensity-based technique, however, is often limited in low-contrast image regions with subtle intensity differences. The solved deformation vector fields (DVFs) can also be biomechanically unrealistic. To address these problems, we have developed a biomechanical modeling guided CBCT estimation technique (Bio-CBCT-est) by combining 2D-3D deformation with finite element analysis (FEA)-based biomechanical modeling of anatomical structures. Specifically, Bio-CBCT-est first extracts the 2D-3D deformation-generated displacement vectors at the high-contrast anatomical structure boundaries. The extracted surface deformation fields are subsequently used as the boundary conditions to drive structure-based FEA to correct and fine-tune the overall deformation fields, especially those at low-contrast regions within the structure. The resulting FEA-corrected deformation fields are then fed back into 2D-3D deformation to form an iterative loop, combining the benefits of intensity-based deformation and biomechanical modeling for CBCT estimation. Using eleven lung cancer patient cases, the accuracy of the Bio-CBCT-est technique has been compared to that of the 2D-3D deformation technique and the traditional CBCT reconstruction techniques. The accuracy was evaluated in the image domain, and also in the DVF domain through clinician-tracked lung landmarks.

Sensitivity of tumor motion simulation accuracy to lung biomechanical modeling approaches and parameters

Finite element analysis (FEA)-based biomechanical modeling can be used to predict lung respiratory motion. In this technique, elastic models and biomechanical parameters are two important factors that determine modeling accuracy. We systematically evaluated the effects of lung and lung tumor biomechanical modeling approaches and related parameters to improve the accuracy of motion simulation of lung tumor center of mass (TCM) displacements. Experiments were conducted with four-dimensional computed tomography (4D-CT). A Quasi-Newton FEA was performed to simulate lung and related tumor displacements between end-expiration (phase 50%) and other respiration phases (0%, 10%, 20%, 30%, and 40%). Both linear isotropic and non-linear hyperelastic materials, including the neo-Hookean compressible and uncoupled Mooney-Rivlin models, were used to create a finite element model (FEM) of lung and tumors. Lung surface displacement vector fields (SDVFs) were obtained by registering the 50% phase CT to other respiration phases, using the non-rigid demons registration algorithm. The obtained SDVFs were used as lung surface displacement boundary conditions in FEM. The sensitivity of TCM displacement to lung and tumor biomechanical parameters was assessed in eight patients for all three models. Patient-specific optimal parameters were estimated by minimizing the TCM motion simulation errors between phase 50% and phase 0%. The uncoupled Mooney-Rivlin material model showed the highest TCM motion simulation accuracy. The average TCM motion simulation absolute errors for the Mooney-Rivlin material model along left-right, anterior-posterior, and superior-inferior directions were 0.80 mm, 0.86 mm, and 1.51 mm, respectively. The proposed strategy provides a reliable method to estimate patient-specific biomechanical parameters in FEM for lung tumor motion simulation.

Improvements and artifact analysis in conductivity images using multiple internal electrodes

Electrical impedance tomography is an attractive functional imaging method. It is currently limited in resolution and sensitivity due to the complexity of the inverse problem and the safety limits of introducing current. Recently, internal electrodes have been proposed for some clinical situations such as intensive care or RF ablation. This paper addresses the research question related to the benefit of one or more internal electrodes usage since these are invasive. Internal electrodes would be able to reduce the effect of insulating boundaries such as fat and bone and provide improved internal sensitivity. We found there was a measurable benefit with increased numbers of internal electrodes in saline tanks of a cylindrical and complex shape with up to two insulating boundary gel layers modeling fat and muscle. The internal electrodes provide increased sensitivity to internal changes, thereby increasing the amplitude response and improving resolution. However, they also present an additional challenge of increasing sensitivity to position and modeling errors. In comparison with previous work that used point sources for the internal electrodes, we found that it is important to use a detailed mesh of the internal electrodes with these voxels assigned to the conductivity of the internal electrode and its associated holder. A study of different internal electrode materials found that it is optimal to use a conductivity similar to the background. In the tank with a complex shape, the additional internal electrodes provided more robustness in a ventilation model of the lungs via air filled balloons.

Real-time estimation of prostate tumor rotation and translation with a kV imaging system based on an iterative closest point algorithm

Previous studies have shown that during cancer radiotherapy a small translation or rotation of the tumor can lead to errors in dose delivery. Current best practice in radiotherapy accounts for tumor translations, but is unable to address rotation due to a lack of a reliable real-time estimate. We have developed a method based on the iterative closest point (ICP) algorithm that can compute rotation from kilovoltage x-ray images acquired during radiation treatment delivery. A total of 11 748 kilovoltage (kV) images acquired from ten patients (one fraction for each patient) were used to evaluate our tumor rotation algorithm. For each kV image, the three dimensional coordinates of three fiducial markers inside the prostate were calculated. The three dimensional coordinates were used as input to the ICP algorithm to calculate the real-time tumor rotation and translation around three axes. The results show that the root mean square error was improved for real-time calculation of tumor displacement from a mean of 0.97 mm with the stand alone translation to a mean of 0.16 mm by adding real-time rotation and translation displacement with the ICP algorithm. The standard deviation (SD) of rotation for the ten patients was 2.3°, 0.89° and 0.72° for rotation around the right-left (RL), anterior-posterior (AP) and superior-inferior (SI) directions respectively. The correlation between all six degrees of freedom showed that the highest correlation belonged to the AP and SI translation with a correlation of 0.67. The second highest correlation in our study was between the rotation around RL and rotation around AP, with a correlation of -0.33. Our real-time algorithm for calculation of rotation also confirms previous studies that have shown the maximum SD belongs to AP translation and rotation around RL. ICP is a reliable and fast algorithm for estimating real-time tumor rotation which could create a pathway to investigational clinical treatment studies requiring real-time measurement and adaptation to tumor rotation.

Measurement of retinal arteriolar diameters from auto scale phase congruency with fuzzy weighting and L1 regularization

Manual measurements of small changes in retinal vascular diameter are slow and may be subject to considerable observer-related biases. Among the conventional automatic methods the sliding linear regression filter (SLRF) demonstrates the least scattered and most repeatable coefficients. For optimal performance it relies on the choice of the correct filter scale for different vessel sizes. A small scale extracts fine details at the expense noise sensitivity, while large scales have poor edge localization. Here we use auto scale phase congruency to select the filter scales with fuzzy weighting to reduce noise, and L1 regularization for edge smoothing. Our method uses a one dimensional analysis normal to the vessel and so is faster than the 2D phase congruency. In 65 vessels randomly selected from 20 images the proposed method showed better repeatability and over three times less scattering than conventional SLRF.

Dry electrode bio-potential recordings

As wireless bio-medical long term monitoring moves towards personal monitoring it demands very high input impedance systems capable to extend the reading of bio-signal during the daily activities offering a kind of "stress free", convenient connection, with no need for skin preparation. In particular we highlight the development and broad applications of our own circuits for wearable bio-potential sensor systems enabled by the use of an FET based amplifier circuit with sufficiently high impedance to allow the use of passive dry electrodes which overcome the significant barrier of gel based contacts. In this paper we present the ability of dry electrodes in long term monitoring of ECG, EEG and fetal ECG.

A comparison between compressed sensing algorithms in electrical impedance tomography

Electrical Impedance Tomography (EIT) calculates the internal conductivity distribution within a body using electrical contact measurements. Conventional EIT reconstruction methods solve a linear model by minimizing the least squares error, i.e., the Euclidian or L2-norm, with regularization. Compressed sensing provides unique advantages in Magnetic Resonance Imaging (MRI) [1] when the images are transformed to a sparse basis. EIT images are generally sparser than MRI images due to their lower spatial resolution. This leads us to investigate ability of compressed sensing algorithms currently applied to MRI in EIT without transformation to a new basis. In particular, we examine four new iterative algorithms for L1 and L0 minimization with applications to compressed sensing and compare these with current EIT inverse L1-norm regularization methods. The four compressed sensing methods are as follows: (1) an interior point method for solving L1-regularized least squares problems (L1-LS); (2) total variation using a Lagrangian multiplier method (TVAL3); (3) a two-step iterative shrinkage / thresholding method (TWIST) for solving the L0-regularized least squares problem; (4) The Least Absolute Shrinkage and Selection Operator (LASSO) with tracing the Pareto curve, which estimates the least squares parameters subject to a L1-norm constraint. In our investigation, using 1600 elements, we found all four CS algorithms provided an improvement over the best conventional EIT reconstruction method, Total Variation, in three important areas: robustness to noise, increased computational speed of at least 40x and a visually apparent improvement in spatial resolution. Out of the four CS algorithms we found TWIST was the fastest with at least a 100x speed increase.