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.

Is noncoplanar plan more robust to inter-fractional variations than coplanar plan in treating bilateral HN tumors with pencil-beam scanning proton beams?

PURPOSE: Noncoplanar plans (NCPs) are commonly used for proton treatment of bilateral head and neck (HN) malignancies. NCP requires additional verification setup imaging between beams to correct residual errors of robotic couch motion, which increases imaging dose and total treatment time. This study compared the quality and robustness of NCPs with those of coplanar plans (CPs). METHODS AND MATERIALS: Under an IRB-approved study, CPs were created retrospectively for 10 bilateral HN patients previously treated with NCPs maintaining identical beam geometry of the original plan but excluding couch rotations. Plan robustness to the inter-fractional variation (IV) of both plans was evaluated through the Dose Volume Histograms (DVH) of weekly quality assurance CT (QACT) sets (39 total). In addition, delivery efficiency for both plans was compared using total treatment time (TTT) and beam-on time (BOT). RESULTS: No significant differences in plan quality were observed in terms of clinical target volume (CTV) coverage (D95) or organ-at-risk (OAR) doses (p > 0.4 for all CTVs and OARs). No significant advantage of NCPs in the robustness to IV was found over CP, either. Changes in D95 of QA plans showed a linear correlation (slope = 1.006, R2  > 0.99) between NCP and CP for three CTV data points (CTV1, CTV2, and CTV3) in each QA plan (117 data points for 39 QA plans). NCPs showed significantly higher beam delivery time than CPs for TTT (539 ± 50 vs. 897 ± 142 s; p < 0.001); however, no significant differences were observed for BOT. CONCLUSION: NCPs are not more robust to IV than CPs when treating bilateral HN tumors with pencil-beam scanning proton beams. CPs showed plan quality and robustness similar to NCPs while reduced treatment time (∼6 min). This suggests that CPs may be a more efficient planning technique for bilateral HN cancer proton therapy.

A study of skin marker alignment using different diamond-shaped light fields for prone breast external-beam radiation therapy.

For breast cancer patients treated in the prone position with tangential fields, a diamond-shaped light field (DSLF) can be used to align with corresponding skin markers for image-guided radiation therapy (IGRT). This study evaluates and compares the benefits of different DSLF setups. Seventy-one patients who underwent daily tangential kilovoltage (kV) IGRT were categorized retrospectively into four groups: (1) DSLF field size (FS) = 10 × 10 cm2 , gantry angle = 90° (right breast)/270° (left breast), with the same isocenter as treatment tangential beams; (2) same as group 1, except DSLF FS = 4 × 4 cm2 ; (3) DSLF FS = 4 × 4-6 × 8 cm2 , gantry angle = tangential treatment beam, off-isocenter so that the DSLF was at the approximate breast center; and (4) No-DSLF. We compared their total setup time (including any DSLF/marker-based alignment and IGRT) and relative kV-based couch shift corrections. For groups 1-3, DSLF-only dose distributions (excluding kV-based correction) were simulated by reversely shifting the couch positions from the computed tomography plans, which were assumed equivalent to the delivered dose when both DSLF and IGRT were used. For patient groups 1-4, the average daily setup time was 2.6, 2.5, 5.0, and 8.3 min, respectively. Their mean and standard deviations of daily kV-based couch shifts were 0.64 ± 0.4, 0.68 ± 0.3, 0.8 ± 0.6, and 1.0 ± 0.6 cm. The average target dose changes after excluding kV-IGRT for groups 1-3 were-0.2%, -0.1%, and +0.4%, respectively, whereas DSLF-1 was most efficient in sparing heart and chest wall, DSLF-2 had lowest lung Dmax ; and DSLF-3 maintained the highest target coverage at the cost of highest OAR dose. In general, the use of DSLF greatly reduces patient setup time and may result in smaller IGRT corrections. If IGRT is limited, different DSLF setups yield different target coverage and OAR dose sparing. Our findings will help DSLF setup optimization in the prone breast treatment setting.

Is a weekly qualitative picket fence test sufficient? A proposed alternate EPID-based weekly MLC QA program.

PURPOSE: Well-designed routine multileaf collimator (MLC) quality assurance (QA) is important to assure external-beam radiation treatment delivery accuracy. This study evaluates the clinical necessity of a comprehensive weekly (C-Weekly) MLC QA program compared to the American Association of Physics in Medicinerecommended weekly picket fence test (PF-Weekly), based on our seven-year experience with weekly MLC QA. METHODS: The C-Weekly MLC QA program used in this study includes 5 tests to analyze: (1) absolute MLC leaf position; (2) interdigitation MLC leaf position; (3) picket fence MLC leaf positions at static gantry angle; (4) minimum leaf-gap setting; and (5) volumetric-modulated arc therapy delivery. A total of 20,226 QA images from 16,855 tests (3,371 tests × 5) for 11 linacs at 5 photon clinical sites from May 2014 to June 2021 were analyzed. Failure mode and effects analysis was performed with 5 failure modes related to the 5 tests. For each failure mode, a risk probability number (RPN) was calculated for a C-Weekly and a PF-Weekly MLC QA program. The probability of occurrence was evaluated from statistical analyses of the C-Weekly MLC QA. RESULTS: The total number of failures for these 16,855 tests was 143 (0.9%): 39 (27.3%) for absolute MLC leaf position, 13 (9.1%) for interdigitation position, 9 (6.3%) for static gantry picket fence, 2 (1.4%) for minimum leaf-gap setting, and 80 (55.9%) for VMAT delivery. RPN scores for PF-Weekly MLC QA ranged from 60 to 192 and from 48 to 96 for C-Weekly MLC QA. CONCLUSION: RPNs for the 5 failure modes of MLC QA tests were quantitatively determined and analyzed. A comprehensive weekly MLC QA is imperative to lower the RPNs of the 5 failure modes to the desired level (<125); those from the PF-Weekly MLC QA program were found to be higher (>125). This supports the clinical necessity for comprehensive weekly MLC QA.

Toward automation of initial chart check for photon/electron EBRT: the clinical implementation of new AAPM task group reports and automation techniques.

PURPOSE: The recently published AAPM TG-275 and the public review version of TG-315 list new recommendations for comprehensive and minimum physics initial chart checks, respectively. This article addresses the potential development and benefit of initial chart check automation when these recommendations are implemented for clinical photon/electron EBRT. METHODS: Eight board-certified physicists with 2-20 years of clinical experience performed initial chart checks using checklists from TG-275 and TG-315. Manual check times were estimated for three types of plans (IMRT/VMAT, 3D, and 2D) and for prostate, whole pelvis, lung, breast, head and neck, and brain cancers. An expert development team of three physicists re-evaluated the automation feasibility of TG-275 checklist based on their experience of developing and implementing the in-house and the commercial automation tools in our institution. Three levels of initial chart check automation were simulated: (1) Auto_UMMS_tool (which consists of in-house program and commercially available software); (2) Auto_TG275 (with full and partial automation as indicated in TG-275); and (3) Auto_UMMS_exp (with full and partial automation as determined by our experts' re-evaluation). RESULTS: With no automation of initial chart checks, the ranges of manual check times were 29-56 min (full TG-315 list) and 102-163 min (full TG-275 list), which varied significantly with physicists but varied little at different tumor sites. The 69 of 71 checks which were considered as "not fully automated" in TG-275 were re-evaluated with more automation feasibility. Compared to no automation, the higher levels of automation yielded a great reduction in both manual check times (by 44%-98%) and potentially residual detectable errors (by 15-85%). CONCLUSION: The initial chart check automation greatly improves the practicality and efficiency of implementing the new TG recommendations. Revisiting the TG reports with new technology/practice updates may help develop and utilize more automation clinically.

A practical cyberattack contingency plan for radiation oncology.

PURPOSE: This article presents a solution for continuing radiation therapy without interruption in the event of a cyberattack to the radiation oncology information systems (ROIS). This process could be easily deployed to any radiation oncology practice, with little clinical overhead or burden. METHODS AND MATERIALS: The solution automatically retrieves all essential information from the clinical ROIS for each under-treatment patient and periodically (e.g., daily) saves these data to a dedicated secure server for recovery. In the event that the clinical ROIS is not functioning as a result of a cyberattack, this essential information is used to build a new secondary ROIS server to continue radiotherapy treatments until the main ROIS is recovered. Once the cyberattack threat is cleared, the clinical ROIS server is rebuilt from the institution's enterprise backup. The newly accumulated treatment information for each patient is then exported from the secondary ROIS to bring the clinical ROIS up to date. RESULTS: The Department of Radiation Oncology at the University of Maryland Medical System implemented this solution for clinical use with the Varian ARIA ROIS in the management of ~250 daily radiotherapy treatments, inclusive of a proton center. This solution was determined to be a feasible and affordable business continuity plan for the radiation oncology practice by minimizing radiation treatment downtime to a couple of hours in a simulated cyberattack drill. CONCLUSIONS: The proposed solution can achieve continuation of radiation therapy treatment without treatment breaks in the event of a cyberattack. It also provides cushion time for radiation oncology departments to rebuild their clinical ROIS systems from the enterprise data backup.

Clinical Implementation and Dosimetric Evaluation of a Robust Proton Lattice Planning Strategy Using Primary and Robust Complementary Beams.

PURPOSE: Pencil-beam scanning proton therapy has been considered as a potential modality for the 3D form of spatially-fractionated-radiation-therapy called lattice therapy. However, few practical solutions have been introduced in the clinic. Existing limitations include degradation in plan quality and robustness when using single-field versus multifield lattice plans, respectively. We propose a practical and robust proton lattice (RPL) planning method using multifield and evaluate its dosimetric characteristics compared to clinically acceptable photon lattice plans. METHODS: Seven cases previously treated with photon lattice therapy were used to evaluate a novel RPL planning technique using two-orthogonal beams: a primary beam (PB) and a robust complementary beam (RCB) that deliver 67% and 33%, respectively, of the prescribed dose to vertices inside the gross-target-volume (GTV). Only RCB is robustly optimized for setup and range uncertainties. The number and volume of vertices, peak-to-valley dose ratios (PVDRs), and volume of low dose to GTV of proton and photon plans were compared. The RPL technique was then used in treatment of two patients and their dosimetric parameters are reported. RESULTS: The RPL strategy was able to achieve the clinical planning goals. Compared to previously-treated photon plans, the average number of vertices increased by 30%, average vertex volume by 49% (18.2±25.9cc vs. 12.2±14.5cc, P=0.21), and higher PVDR (10.5±4.8 vs. 2.5±0.9, P<0.005) was achieved. In addition, RPL plans show more conformal dose with less low-dose to GTV (V30%: 60.9±7.2% vs. 81.6±13.9% and V10%: 88.3±4.5% vs. 98.6±3.6% [P<0.01]). The RPL plan for two treated patients showed PVDRs of 4.61 and 14.85 with vertices-to-GTV ratios of 1.52% and 1.30%, respectively. CONCLUSION: A novel RPL planning strategy using a pair of orthogonal beams was developed and successfully translated to the clinic. The proposed method can generate better quality plans, a higher number of vertices, and higher PVDRs than currently used photon lattice plans.

Readiness for Radiation Treatment Continuity: Survey on Contingency Plans Against Cyberattacks.

Purpose: Cyberattacks on health care systems have been on the rise over the past 5 years. Formulation and implementation of a robust postattack business continuity plan and/or contingency plan (CP) is essential for minimal disruption to patient care. The level of awareness and planning within the radiation oncology community for cyberattacks is not clear. This study was undertaken to survey and assess cyberattack CP awareness and preparedness. Methods and Materials: A survey instrument comprising 5 questions on awareness and preparedness of cyberattack CPs was e-mailed to 150 radiation oncology departments. Recipients included 105 institutions with residency programs in therapeutic medical physics, as listed by the Commission on Accreditation of Medical Physics Education Program (usually either school-based or large institutional settings), and 45 additional smaller settings within the United States, representing community practices. Results: Forty-three responses were deemed evaluable for analysis. Forty-two percent (18 respondents) of respondents responded that they are well-aware of the concept of a cyberattack CP. A large discrepancy in awareness exists between larger hospitals (LH) that have 5 or more treatment machines and smaller hospitals (SH) that have 4 or fewer, 54% versus 24 % (P < .05). Fifty-eight percent of respondents considered it "essential" to have such a plan in place, and 28% considered it "desirable" to do so but not practical. Nine percent regarded a cyberattack CP as unnecessary. No significant differences in responses were noted among different types or sizes of institutions on this issue. Sixty-two percent of LH responded that they were either preparing or evaluating a CP, compared with only 29% of SH (P = .03). However, no respondents explicitly replied that they already had a CP in place in their practices. Conclusions: The importance of cyberattack preparedness and implementation does not seem to be well-recognized in radiation oncology. Both the awareness and the preparedness of SH are substantially less than those of LH. Specific and ongoing education efforts in parallel with development of appropriate programs are needed to counter the increasingly pervasive and complex threat of cyberattacks.

Use of CBCT plus plan robustness for reducing QACT frequency in intensity-modulated proton therapy: Head-and-neck cases.

PURPOSE: Anatomic variation has a significant dosimetric impact in intensity-modulated proton therapy. Weekly or biweekly computed tomography (CT) scans, called quality assurance CTs (QACTs), are used to monitor anatomic and resultant dose changes to determine whether adaptive plans are needed. Frequent CT scans result in unwanted QACT dose and increased clinical workloads. This study proposed utilizing patient setup cone-beam CTs (CBCTs) and treatment plan robustness to reduce the frequency of QACTs. METHODS: We retrospectively analyzed data from 27 patients with head-and-neck cancer, including 594 CBCTs, 136 QACTs, and 19 adaptive plans. For each CBCT, water-equivalent thickness (WET) along the pencil-beam path was calculated. For each treatment plan, the WET of the first-day CBCT was used as the reference, and the mean WET changes (ΔWET) in each following CBCT was used as the surrogate of proton range change. Using CBCTs acquired prior to a QACT, we predicted the ΔWET on the QACT day by a linear regression model. The impact of range change on target dose was calculated as the predicted ΔWET weighted by the monitor units of each field. In addition, plan robustness was estimated from the robust dose-volume histograms (DVHs) and combined with ΔWET to reduce QACT frequency. Robustness was estimated from the distance between the DVH curves of the nominal and worst scenarios. RESULTS: When the estimated mean ΔWET was <6.5 mm (or <7.5 mm if the robustness was >95%), the QACT could be skipped without missing any adaptive planning; otherwise a QACT was required. Overall, 41% of QACTs could be eliminated when ΔWET was <6.5 mm and 56% when ΔWET was <7.5 mm, and robustness was >95%. At least one QACT could have been omitted in 25 of the 27 cases under skipping thresholds at ΔWETs <7.5 mm and R > 95%. CONCLUSION: This study suggests that the number of QACTs can be greatly reduced by calculating range change in patient setup CBCTs and can be further reduced by combining this information with analyses of plan robustness.

Semi-Supervised Segmentation of Radiation-Induced Pulmonary Fibrosis From Lung CT Scans With Multi-Scale Guided Dense Attention.

Computed Tomography (CT) plays an important role in monitoring radiation-induced Pulmonary Fibrosis (PF), where accurate segmentation of the PF lesions is highly desired for diagnosis and treatment follow-up. However, the task is challenged by ambiguous boundary, irregular shape, various position and size of the lesions, as well as the difficulty in acquiring a large set of annotated volumetric images for training. To overcome these problems, we propose a novel convolutional neural network called PF-Net and incorporate it into a semi-supervised learning framework based on Iterative Confidence-based Refinement And Weighting of pseudo Labels (I-CRAWL). Our PF-Net combines 2D and 3D convolutions to deal with CT volumes with large inter-slice spacing, and uses multi-scale guided dense attention to segment complex PF lesions. For semi-supervised learning, our I-CRAWL employs pixel-level uncertainty-based confidence-aware refinement to improve the accuracy of pseudo labels of unannotated images, and uses image-level uncertainty for confidence-based image weighting to suppress low-quality pseudo labels in an iterative training process. Extensive experiments with CT scans of Rhesus Macaques with radiation-induced PF showed that: 1) PF-Net achieved higher segmentation accuracy than existing 2D, 3D and 2.5D neural networks, and 2) I-CRAWL outperformed state-of-the-art semi-supervised learning methods for the PF lesion segmentation task. Our method has a potential to improve the diagnosis of PF and clinical assessment of side effects of radiotherapy for lung cancers.

Assessing the Need for Adjusted Organ-at-Risk Planning Goals for Patients Undergoing Adjuvant Radiation Therapy for Locally Advanced Breast Cancer with Proton Radiation.

PURPOSE: Locally advanced breast cancer requires surgical management via lumpectomy or mastectomy with or without systemic therapy followed by chest wall or breast (CW) and comprehensive nodal irradiation (CNI). Radiation (RT) dose constraints for the heart and ipsilateral lung have been developed based on photon RT. Proton therapy (PBT) can deliver significantly lower doses of RT to these organs-at-risk (OARs) and may warrant adjustments to OAR planning goals. METHODS AND MATERIALS: The RT plans of consecutive patients undergoing adjuvant CW-CNI RT with PBT within a single center were reviewed. A inital treatment volume, comprised of CW/intact breast + CNI (CTV_init) structure, including the CW and CNI but excluding any boost plans was analyzed. Frequency distributions were generated based on doses received by the heart, lungs, and esophagus for validated dosimetric parameters. Frequency distributions were generated and then stratified by laterality and compared using the Kruskal-Wallis H test. The 75th, 85th, and 95th percentiles for each dosimetric parameter were calculated, overall and by laterality. The 75th percentile (Q3), was used as a suggested primary goal, and the 95th percentile was used as a suggested secondary goal. RESULTS: One hundred and seventy-two plans were analyzed. Forty-nine plans were right-sided, 107 were left-sided, and 16 were bilateral. The overall Q3 of the mean and V25 of the heart were 1.5 Gy and 1.7%, respectively. The mean and V25 to the heart differed significantly by laterality. Pulmonary values were similar to current recommendations. For all lateralites, the median volume of the esophagus receiving 70% prescription dose was ≤1 cm3. CONCLUSIONS: We present the first dosimetric study providing complete OAR dose-volume histograms data for patients undergoing adjuvant pencil-beam scanning-PBT for locally advanced breast cancer, with detailed information on central tendencies, ranges and distributions of data. We have provided suggested planning goals and metrics for the lungs, heart, and esophagus; the latter 2 differing significantly from current Quantitative Analysis of Normal Tissue Effects in the Clinic (QUANTEC) constraints and classical photon goals.

A systematic quality assurance framework for the upgrade of radiation oncology information systems.

In spite of its importance, no systematic and comprehensive quality assurance (QA) program for radiation oncology information systems (ROIS) to verify clinical and treatment data integrity and mitigate against data errors/corruption and/or data loss risks is available. Based on data organization, format and purpose, data in ROISs falls into five different categories: (1) the ROIS relational database and associated files; (2) the ROIS DICOM data stream; (3) treatment machine beam data and machine configuration data; (4) electronic medical record (EMR) documents; and (5) user-generated clinical and treatment reports from the ROIS. For each data category, this framework proposes a corresponding data QA strategy to very data integrity. This approach verified every bit of data in the ROIS, including billions of data records in the ROIS SQL database, tens of millions of ROIS database-associated files, tens of thousands of DICOM data files for a group of selected patients, almost half a million EMR documents, and tens of thousands of machine configuration files and beam data files. The framework has been validated through intentional modifications with test patient data. Despite the 'big data' nature of ROIS, the multiprocess and multithread nature of our QA tools enabled the whole ROIS data QA process to be completed within hours without clinical interruptions. The QA framework suggested in this study proved to be robust, efficient and comprehensive without labor-intensive manual checks and has been implemented for our routine ROIS QA and ROIS upgrades.

Spatially Fractionated Radiotherapy (GRID) Prior to Standard Neoadjuvant Conventionally Fractionated Radiotherapy for Bulky, High-Risk Soft Tissue and Osteosarcomas: Feasibility, Safety, and Promising Pathologic Response Rates.

Spatially fractionated radiotherapy (GRID) has been utilized primarily in the palliative and definitive treatment of bulky tumors. Delivered in the modern era primarily with megavoltage photon therapy, this technique offers the promise of safe dose escalation with potential immunogenic, bystander and microvasculature effects that can augment a conventionally fractionated course of radiotherapy. At the University of Maryland, an institutional standard has arisen to incorporate a single fraction of GRID radiation in large (>8 cm), high-risk soft tissue and osteosarcomas prior to a standard fractionated course. Herein, we report on the excellent pathologic responses and apparent safety of this regimen in 26 consecutive patients.

Action Levels on Dose and Anatomic Variation for Adaptive Radiation Therapy Using Daily Offline Plan Evaluation: Preliminary Results.

PURPOSE: This study aimed to develop action levels for replanning to accommodate dosimetric variations resulting from anatomic changes during the course of treatments, using daily cone beam computed tomography (CBCT). METHODS AND MATERIALS: Daily or weekly CBCT images of 20 patients (10 head and neck, 5 lung, and 5 prostate cancers) who underwent resimulation per physicians' clinical decisions, mainly from the comparison of CBCT scans, were used to determine action levels. The first CBCT image acquired before the first treatment was used as the reference image to rule out effects of dose inaccuracy from the CBCT. The Pearson correlation of clinical target volume (CTV) was used as a parameter of anatomic variation. Parameters for action levels on dose and anatomic variation were deduced by comparing the parameters and clinical decisions made for replanning. A software tool was developed to automatically perform all procedures, including dose calculations, using the CBCT and plan evaluations. RESULTS: Replans were clinically decided based on either significant dose or anatomic changes in 13 cases. The 7 cases that did not require replanning showed dose differences <5%, and the Pearson correlation of the CTV was >75% for all fractions. A difference in planning target volume dose >5% or a difference in the image correlation coefficient of the CTV <0.75 proved to be indicators for replanning. Once the results of the CBCT plan met the replanning criteria, the software tool automatically alerted the attending physician and physicist by both e-mail and pager so that the case could be examined closely. CONCLUSIONS: Our study shows that a dose difference of 5% and/or anatomy variation at 0.75 Pearson correlations are practical action levels on dose and anatomic variation for replanning for the given data sets.

Automated compromised right lung segmentation method using a robust atlas-based active volume model with sparse shape composition prior in CT.

To resolve challenges in image segmentation in oncologic patients with severely compromised lung, we propose an automated right lung segmentation framework that uses a robust, atlas-based active volume model with a sparse shape composition prior. The robust atlas is achieved by combining the atlas with the output of sparse shape composition. Thoracic computed tomography images (n=38) from patients with lung tumors were collected. The right lung in each scan was manually segmented to build a reference training dataset against which the performance of the automated segmentation method was assessed. The quantitative results of this proposed segmentation method with sparse shape composition achieved mean Dice similarity coefficient (DSC) of (0.72, 0.81) with 95% CI, mean accuracy (ACC) of (0.97, 0.98) with 95% CI, and mean relative error (RE) of (0.46, 0.74) with 95% CI. Both qualitative and quantitative comparisons suggest that this proposed method can achieve better segmentation accuracy with less variance than other atlas-based segmentation methods in the compromised lung segmentation.

A Voxel-by-Voxel Comparison of Deformable Vector Fields Obtained by Three Deformable Image Registration Algorithms Applied to 4DCT Lung Studies.

BACKGROUND: Commonly used methods of assessing the accuracy of deformable image registration (DIR) rely on image segmentation or landmark selection. These methods are very labor intensive and thus limited to relatively small number of image pairs. The direct voxel-by-voxel comparison can be automated to examine fluctuations in DIR quality on a long series of image pairs. METHODS: A voxel-by-voxel comparison of three DIR algorithms applied to lung patients is presented. Registrations are compared by comparing volume histograms formed both with individual DIR maps and with a voxel-by-voxel subtraction of the two maps. When two DIR maps agree one concludes that both maps are interchangeable in treatment planning applications, though one cannot conclude that either one agrees with the ground truth. If two DIR maps significantly disagree one concludes that at least one of the maps deviates from the ground truth. We use the method to compare 3 DIR algorithms applied to peak inhale-peak exhale registrations of 4DFBCT data obtained from 13 patients. RESULTS: All three algorithms appear to be nearly equivalent when compared using DICE similarity coefficients. A comparison based on Jacobian volume histograms shows that all three algorithms measure changes in total volume of the lungs with reasonable accuracy, but show large differences in the variance of Jacobian distribution on contoured structures. Analysis of voxel-by-voxel subtraction of DIR maps shows differences between algorithms that exceed a centimeter for some registrations. CONCLUSION: Deformation maps produced by DIR algorithms must be treated as mathematical approximations of physical tissue deformation that are not self-consistent and may thus be useful only in applications for which they have been specifically validated. The three algorithms tested in this work perform fairly robustly for the task of contour propagation, but produce potentially unreliable results for the task of DVH accumulation or measurement of local volume change. Performance of DIR algorithms varies significantly from one image pair to the next hence validation efforts, which are exhaustive but performed on a small number of image pairs may not reflect the performance of the same algorithm in practical clinical situations. Such efforts should be supplemented by validation based on a longer series of images of clinical quality.

Optical measurement of azimuthal anchoring strength in nematic liquid crystals.

We have observed the azimuthal switching at the interface of nematic liquid crystals (LCs) and the aligning substrate, induced by the planar electric field of a fine comb electrode. Optical transmittance as a function of applied voltage was modeled by both the elastic theory with rigid anchoring and the Landau-de Gennes theory with the interfacial energy expression 1/2Wtr[(Q two down arrows -Q two down arrows (0))(2)], where Q two down arrows is the liquid crystal order parameter and Q two down arrows (0) is the surface order parameter induced by the aligning substrate. Optical data on the in-plane switching LC cells were found to differ qualitatively with the predictions of the rigid anchoring model but to agree well with those of the Landau-de Gennes theory. We obtain not only the strength W of the azimuthal anchoring, but also find the surface order parameter S(surface) to be 20-30 % less than that of the bulk. The optically measured azimuthal anchoring strength is in good agreement with the literature values determined through other means.

First-order liquid crystal orientation transition on inhomogeneous substrates.

In a recent experiment, we uncovered an unconventional liquid crystal (LC) orientation transition on microtextured substrates consisting of alternating horizontal and vertical corrugations. When the period of alternation was decreased toward approximately 1 microm, the LC alignment underwent an abrupt transition from inhomogeneous planar to a more uniform configuration with a large pretilt angle ( approximately 40 degrees ). With the aid of a model based on the competition between the Frank-Oseen elastic energy and a phenomenological surface potential of the form W(theta,phi)=(1/2)W((2))(theta) sin(2) theta+(1/4)W((4))(theta) sin(4) theta+(1/2)W(phi) cos(2) theta sin(2) phi(x,y) (where theta and phi are, respectively, the pretilt and azimuthal angles of the LC director and W((2))(theta), W((4))(theta), and W(phi) are constants) that demonstrated good agreement with experiment, we investigated the microscopic origin of the observed transition. It was found that this transition comprises two steps. First, the LC director homogenizes toward the phi=45 degrees azimuthal direction in the plane to relax the elastic energy. The resulting rise in azimuthal anchoring energy subsequently drives the LC to adopt a finite pretilt. The values of the W's deduced from the model reveal that the polar anchoring energy is about approximately 1/10 of the typical values, with the sin(4) theta term dominating the sin(2) theta term. We present a possible explanation for this unexpected finding.

Liquid crystal orientation transition on microtextured substrates.

A uniform alignment of liquid crystal (LC) with finite pretilt was observed on microtextured substrates that were lithographically fabricated with alternating horizontal and vertical corrugations. As the period of alternation was decreased toward 0.8 microm, the nematic LC alignment on these substrates changed from inhomogeneous in plane, copying the substrate corrugations, to a uniform configuration with a large pretilt of approximately 40 degrees. This transition is pertinent to a frustrated boundary wherein a lowering in the LC elastic energy due to spatial variation in the LC orientation compromises an increase in the surface anchoring energy. A model based on this idea demonstrates good agreement with the experiment. This result may open up a new arena for tailoring substrate characteristics for LC alignment.