Investigation on the accuracy of field widths and the quantitative relationship with energy variations for dynamic jaws helical tomotherapy

BACKGROUND AND PURPOSE

TomoEDGE is an advanced technology for TomoTherapy treatment delivery by introducing a sliding-window dynamic jaw motion. The front and back jaws move independently at the start and end of a target volume along the longitudinal couch direction to reduce the undesired dose to the normal tissues. The accuracy of field width is essential to treatment delivery in this regard. The purpose of this work was to analyze the performance of dynamic jaws on helical tomotherapy and investigate the relationship with energy variation.

METHODS

The Tomotherapy-Quality-Assurance (TQA) Dynamic Field Width procedure was performed monthly across three tomotherapy machines. All field widths were analyzed, especially the FWHM of the 10 mm field width. Field width measurements were compared with the ratio of Percentage Depth Dose at 20 and 10 cm to render the value of correlation. Changes in beam FWHM and energy were further discussed. Two-year data were collected for this purpose.

RESULTS

On average, measured field widths in each unit agreed within 1% tolerance recommendation stated. The average absolute difference between reference and measured FWs in each unit was approximately 0.07 mm. An increase of 1.5% in the FW of the 10 mm nominal beam width was correlated with a 1% increase in PDD20,10 ratio, implying a positive correlation between the two factors (p < 0.002).

CONCLUSIONS

A positive correlation between nominal 10 mm FW and PDD20,10 was observed. In the case that the PDD20,10 marginally passes the QA tests, users are recommended to consider further verification on Dynamic Jaws to ensure the smallest field width to be within tolerance, which is essential to maintain effective treatment in TomoEDGE system. Since the regression of this study was a single-factor model, other confounding factors such as the focal spot size of linear accelerator should also be considered when evaluating the machine status.

Prostate diffusion-weighted imaging (DWI) in MR-guided radiotherapy: Reproducibility assessment on 1.5 T MR-Linac and 1.5 T MR-simulator

PURPOSE

Diffusion-weighted imaging (DWI) holds promise for image-guided radiotherapy (MRgRT) in prostate cancer. However, challenges persist due to image distortion, artifacts, and apparent diffusion coefficient (ADC) reproducibility issues. This study aimed to assess DWI image quality and ADC reproducibility on both a 1.5 T MR-simulator and a 1.5 T MR-Linac, employing measurements from both an ACR MRI phantom and prostate cancer patients undergoing MRgRT.

METHODS

DW-MRI scans were conducted on 19 patients (mean age = 69 ± 8 years, with 23 MR-visible intra-prostatic lesions) and an ACR MRI phantom using a 1.5 T MR-simulator (b-values = 0, 800, 1400s/mm2) and a 1.5 T MR-Linac (b-values = 50, 500, 800 s/mm2). ADC homogeneity in the phantom was evaluated via 1D profile flatness (FL) in three directions. Image quality was assessed through qualitative 5-point Likert scale ratings and quantitative ADC and signal-to-noise ratio (SNR) measurements. Intra-observer reproducibility of image quality scores was evaluated using ICC(1, 2). Geometric distortion was measured by comparing landmark sizes on the ACR phantom against the ground truth. Mean ADC and reproducibility were assessed using Bland-Altman plots.

RESULTS

Both MR-simulator and MR-Linac demonstrated high ADC homogeneity (FL > 87.5% - MR-simulator: 97.23 ± 0.62%, MR-Linac: 94.75 ± 0.62%, p < 0.05) in the phantom. Image quality scores revealed acceptable ratings (≥3) for capsule demarcation, image artifacts, and geometric distortion in patients. However, intra-prostatic lesions were barely discernible in b800 images for both MR-simulator (average score = 2.37 ± 1.33) and MR-Linac (average score = 2.16 ± 1.28). While MR-Linac DWI scans exhibited significantly more severe geometric distortion than MR-simulator scans (p < 0.01), most phantom measurements fell within the image in-plane resolution of 3 mm. Significant differences were noted in MR-simulator ADC (CTV: 1.20 ± 0.14 × 10-3 mm2/s (MR-simulator) vs 1.06 ± 0.10 × 10-3 mm2/s (MR-Linac); GTV: 1.05 ± 0.21 × 10-3 mm2/s vs 0.91 ± 0.16 × 10 mm2/s, all p < 0.05), with a small non-zero bias observed in the Bland-Altman analysis (CTV: 12.3%; GTV: 14.5%).

CONCLUSION

The significantly larger MR-simulator ADC and the small non-zero bias hint at potential systematic differences in ADC values acquired from an MR-simulator and an MR-Linac, both at 1.5 T. Although acceptable ADC homogeneity was noted, caution is warranted in interpreting MR-Linac DWI images due to occasional severe artifacts. Further studies are essential to validate DWI and ADC as reliable imaging markers in prostate cancer MRgRT.

One-year clinical outcomes of MR-guided stereotactic body radiation therapy with rectal spacer for patients with localized prostate cancer

BACKGROUND AND PURPOSE

This prospective study aimed to investigate adaptive magnetic resonance (MR)-guided stereotactic body radiation therapy (MRgSBRT) with rectal spacer for localized prostate cancer (PC) and report 1-year clinical outcomes.

MATERIALS AND METHODS

Thirty-four consecutive patients with low- to high-risk localized PC that underwent 5-fraction adaptive MRgSBRT with rectal spacer were enrolled. The dosimetric comparison was performed on a risk- and age-matched cohort treated with MRgSBRT but without a spacer at a similar timepoint. Clinician-reported outcomes were based on Common Terminology Criteria for Adverse Events. Patient-reported outcomes were based on the Expanded Prostate Cancer Index Composite (EPIC) questionnaire at baseline, acute (1-3 months), subacute (4-12 months), and late (> 12 months) phases.

RESULTS

The median follow-up was 390 days (range 28-823) and the median age was 70 years (range 58-82). One patient experienced rectal bleeding soon after spacer insertion that subsided before MRgSBRT. The median distance between the midline of the prostate midgland and the rectum after spacer insertion measured 7.8 mm (range 2.6-15.3), and the median length of the spacer was 45.9 mm (range 16.8-62.9) based on T2-weighted MR imaging. The use of spacer resulted in significant improvements in target coverage (V100% > 95% = 98.6% [range 93.4-99.8] for spacer vs. 97.8% [range 69.6-99.7] for non-spacer) and rectal sparing (V95% < 3 cc = 0.7 cc [range 0-4.6] for spacer vs. 4.9 cc [range 0-12.5] for non-spacer). Nine patients (26.5%) experienced grade 1 gastrointestinal toxicities, and no grade ≥ 2 toxicities were observed. During the 1-year follow-up period, EPIC scores for the bowel domain remained stable and were the highest among all other domains.

CONCLUSIONS

MRgSBRT with rectal spacer for localized PC showed exceptional tolerability with minimal gastrointestinal toxicities and satisfactory patient-reported outcomes. Improvements in dosimetry, rectal sparing, and target coverage were achieved with a rectal spacer. Randomized trials are warranted for further validation.

A study of Bayesian deep network uncertainty and its application to synthetic CT generation for MR-only radiotherapy treatment planning

BACKGROUND

The use of synthetic computed tomography (CT) for radiotherapy treatment planning has received considerable attention because of the absence of ionizing radiation and close spatial correspondence to source magnetic resonance (MR) images, which have excellent tissue contrast. However, in an MR-only environment, little effort has been made to examine the quality of synthetic CT images without using the original CT images.

PURPOSE

To estimate synthetic CT quality without referring to original CT images, this study established the relationship between synthetic CT uncertainty and Bayesian uncertainty, and proposed a new Bayesian deep network for generating synthetic CT images and estimating synthetic CT uncertainty for MR-only radiotherapy treatment planning.

METHODS AND MATERIALS

A novel deep Bayesian network was formulated using probabilistic network weights. Two mathematical expressions were proposed to quantify the Bayesian uncertainty of the network and synthetic CT uncertainty, which was closely related to the mean absolute error (MAE) in Hounsfield Unit (HU) of synthetic CT. These uncertainties were examined to demonstrate the accuracy of representing the synthetic CT uncertainty using a Bayesian counterpart. We developed a hybrid Bayesian architecture and a new data normalization scheme, enabling the Bayesian network to generate both accurate synthetic CT and reliable uncertainty information when probabilistic weights were applied. The proposed method was evaluated in 59 patients (13/12/32/2 for training/validation/testing/uncertainty visualization) diagnosed with prostate cancer, who underwent same-day pelvic CT- and MR-acquisitions. To assess the relationship between Bayesian and synthetic CT uncertainties, linear and non-linear correlation coefficients were calculated on per-voxel, per-tissue, and per-patient bases. For accessing the accuracy of the CT number and dosimetric accuracy, the proposed method was compared with a commercially available atlas-based method (MRCAT) and a U-Net conditional-generative adversarial network (UcGAN).

RESULTS

The proposed model exhibited 44.33 MAE, outperforming UcGAN 52.51 and MRCAT 54.87. The gamma rate (2%/2 mm dose difference/distance to agreement) of the proposed model was 98.68%, comparable to that of UcGAN (98.60%) and MRCAT (98.56%). The per-patient and per-tissue linear correlation coefficients between the Bayesian and synthetic CT uncertainties ranged from 0.53 to 0.83, implying a moderate to strong linear correlation. Per-voxel correlation coefficients varied from -0.13 to 0.67 depending on the regions-of-interest evaluated, indicating tissue-dependent correlation. The R2 value for estimating MAE solely using Bayesian uncertainty was 0.98, suggesting that the uncertainty of the proposed model was an ideal candidate for predicting synthetic CT error, without referring to the original CT.

CONCLUSION

This study established a relationship between the Bayesian model uncertainty and synthetic CT uncertainty. A novel Bayesian deep network was proposed to generate a synthetic CT and estimate its uncertainty. Various metrics were used to thoroughly examine the relationship between the uncertainties of the proposed Bayesian model and the generated synthetic CT. Compared with existing approaches, the proposed model showed comparable CT number and dosimetric accuracies. The experiments showed that the proposed Bayesian model was capable of producing accurate synthetic CT, and was an effective indicator of the uncertainty and error associated with synthetic CT in MR-only workflows.

A pilot study of MRI radiomics for high-risk prostate cancer stratification in 1.5 T MR-guided radiotherapy

PURPOSE

To investigate the potential value of MRI radiomics obtained from a 1.5 T MRI-guided linear accelerator (MR-LINAC) for D'Amico high-risk prostate cancer (PC) classification in MR-guided radiotherapy (MRgRT).

METHODS

One hundred seventy-six consecutive PC patients underwent 1.5 T MRgRT treatment were retrospectively enrolled. Each patient received one or two pretreatment T₂ -weighted MRI scans on a 1.5 T MR-LINAC. The endpoint was to differentiate high-risk from low/intermediate-risk PC based on D'Amico criteria using MRI-radiomics. Totally 1023 features were extracted from clinical target volume (CTV) and planning target volume (PTV). Intraclass correlation coefficient of scan-rescan repeatability, feature correlation, and recursive feature elimination were used for feature dimension reduction. Least absolute shrinkage and selection operator regression was employed for model construction. Receiver operating characteristic area under the curve (AUC) analysis was used for model performance assessment in both training and testing data.

RESULTS

One hundred and eleven patients fulfilled all criteria were finally included: 76 for training and 35 for testing. The constructed MRI-radiomics models extracted from CTV and PTV achieved the AUC of 0.812 and 0.867 in the training data, without significant difference (P = 0.083). The model performances remained in the testing. The sensitivity, specificity, and accuracy were 85.71%, 64.29%, and 77.14% for the PTV-based model; and 71.43%, 71.43%, and 71.43% for the CTV-based model. The corresponding AUCs were 0.718 and 0.750 (P = 0.091) for CTV- and PTV-based models.

CONCLUSION

MRI-radiomics obtained from a 1.5 T MR-LINAC showed promising results in D'Amico high-risk PC stratification, potentially helpful for the future PC MRgRT. Prospective studies with larger sample sizes and external validation are warranted for further verification.

Magnetic field induced dose effects in radiation therapy using MR-linacs

MR-guided radiotherapy (MRgRT) is one of the most significant advances in radiotherapy in recent years. The hybrid systems were designed to visualize patient anatomical and physiological changes during the course of radiotherapy, enabling more precise treatment. However, before MR-linacs reach their full potential in delivering safe and accurate treatments to patients, the radiotherapy team must understand how a magnetic field alters the dosimetric properties of the radiation beam and its potential impact on treatment quality and clinical outcomes. This review aims to provide an in-depth description of the magnetic field induced dose effects for the two widely available systems, the 0.35 T and the 1.5 T MR-linacs. In MR-linac treatments, the primary photon beam passes through MR components that never exist in conventional linacs, which alter both in-field and out-of-field doses. More importantly, the interplay between the always-on magnetic field and the secondary electrons is not negligible. This interplay affects dose deposition in the patient, resulting in reduced in-field skin dose due to purged-out contaminant electrons, shortened build-up distance and a shifted crossline profile owing to asymmetric dose kernel. Especially two effects, namely electron return effect (ERE) and electron stream effect (ESE), are not seen in conventional radiotherapy. This review also summarizes the clinical observations on the site-specific treatments influenced mostly by the magnetic field. In MR-linac treatment, the head and neck region is one of the most challenging sites as ERE occurs at low and high density tissue interfaces and around air cavities, generating hot and cold spots. In breast cancer treatment, consideration should be given to the increased in-field skin dose induced by ERE and the increased out-of-field dose caused by ESE for regions such as the ears, chin, and neck. In lung cancer treatments, tissue inhomogeneity combined with ERE will exacerbate target dose heterogeneity and increase or decrease interface dose. Lastly, treatment in the abdomen and pelvic region will be affected by the presence of gas pockets near the target. The review provides practical recommendations to mitigate these effects. This article is protected by copyright. All rights reserved.

Determining the reliable feature change in longitudinal radiomics studies: A methodological approach using the reliable change index

PURPOSE

Determination of reliable change of radiomics feature over time is essential and vital in delta-radiomics, but has not yet been rigorously examined. This study attempts to propose a methodological approach using reliable change index (RCI), a statistical metric to determine the reliability of quantitative biomarker changes by accounting for the baseline measurement standard error, in delta-radiomics. The use of RCI was demonstrated with the MRI data acquired from a group of prostate cancer (PCa) patients treated by 1.5 T MRI-guided radiotherapy (MRgRT).

METHODS

Fifty consecutive PCa patients who underwent five-fractionated MRgRT were retrospectively included, and 1023 radiomics features were extracted from the clinical target volume (CTV) and planning target volume (PTV). The two MRI datasets acquired at the first fraction (MRI11 and MRI21) were used to calculate the baseline feature reliability against image acquisition using intraclass correlation coefficient (ICC). The RCI was constructed based on the baseline feature measurement standard deviation, ICC, and feature value differences at two time points between the fifth (MRI51) and the first fraction MRI (MRI11). The reliable change of features was determined in each patient only if the calculated RCI was over 1.96 or smaller than -1.96. The feature changes between MRI51 and MRI11 were correlated to two patient-reported quality-of-life clinical endpoints of urinary domain summary score (UDSS) and bowel domain summary score (BDSS) in 35 patients using the Spearman correlation test. Only the significant correlations between a feature that was reliably changed in ≥7 patients (20%) by RCI and an endpoint were considered as true significant correlations.

RESULTS

The 352 (34.4%) and 386 (37.7%) features among all 1023 features were determined by RCI to be reliably changed in more than five (10%) patients in the CTV and PTV, respectively. Nineteen features were found reliably changed in the CTV and 31 features in the PTV, respectively, in 10 (20%) or more patients. These features were not necessarily associated with significantly different longitudinal feature values (group p-value < 0.05). Most reliably changed features in more than 10 patients had excellent or good baseline test-retest reliability ICC, while none showed poor reliability. The RCI method ruled out the features to be reliably changed when substantial feature measurement bias was presented. After applying the RCI criterion, only four and five true significant correlations were confirmed with UDSS and BDSS in the CTV, respectively, with low true significance correlation rates of 10.8% (4/37) and 17.9% (5/28). No true significant correlations were found in the PTV.

CONCLUSIONS

The RCI method was proposed for delta-radiomics and demonstrated using PCa MRgRT data. The RCI has advantages over some other statistical metrics commonly used in the previous delta-radiomics studies, and is useful to reliably identify the longitudinal radiomics feature change on an individual basis. This proposed RCI method should be helpful for the development of essential feature selection methodology in delta-radiomics.

Analysis of online plan adaptation for 1.5T magnetic resonance-guided stereotactic body radiotherapy (MRgSBRT) of prostate cancer

PURPOSE

To analyze and characterize the online plan adaptation of 1.5T magnetic resonance-guided stereotactic body radiotherapy (MRgSBRT) of prostate cancer (PC).

METHODS

PC patients (n = 107) who received adaptive 1.5 Tesla MRgSBRT were included. Online plan adaptation was implemented by either the adapt-to-position (ATP) or adapt-to-shape (ATS) methods. Patients were assigned to the ATS group if they underwent ≥ 1 ATS fraction (n = 51); the remainder were assigned to the ATP group (n = 56). The online plan adaptation records of 535 (107 × 5) fractions were retrospectively reviewed. Rationales for ATS decision-making were determined and analyzed using predefined criteria. Statistics of ATS fractions were summarized. Associations of patient characteristics and clinical factors with ATS utilization were investigated.

RESULTS

There were 87 (16.3%) ATS fractions and 448 ATP fractions (83.7%). The numbers of ATS adoptions in fractions 1-5 were 29 (29/107, 27.1%), 18 (16.8%), 15 (14.0%), 16 (15.0%), and 9 (8.4%), respectively, with significant differences in adoption frequency between fractions (p = 0.007). Other baseline patient characteristics and clinical factors were not significantly associated with ATS classification (all p > 0.05). Underlying criteria for the determination of ATS implementation comprised anatomical changes (77 fractions in 50 patients) and discrete multiple targets (15 fractions in 3 patients). No ATS utilization was determined using dosimetric or online quality assurance criteria.

CONCLUSIONS

This study contributes to facilitating the establishment of a standardized protocol for online MR-guided adaptive radiotherapy in PC.

The impact of the ClearRT™ upgrade on target motion tracking accuracy in Radixact® Synchrony® lung treatments

Background

The objective was to investigate the change in segmentation error of Radixact^® Synchrony^® lung treatment after its kV imaging system was upgraded from Generation 1 to Generation 2 in the ClearRT™ installation.

Materials and methods

Radixact^® Lung Synchrony^® plans were created for the Model 18023 Xsight^® Lung Tracking "XLT" Phantom combined with different lung target inserts with densities of 0.280, 0.500, 0.943 and 1.093 g/cc. After Radixact^® Synchrony^® treatment delivery using the Generation 1 and Generation 2 kV systems according to each plan, the tracking performance of the two kV systems on each density insert was compared by calculating the root mean square (RMS) error (δ_RMS) between the Synchrony-predicted motion in the log file and the known phantom motion and by calculating δ95%, the maximum error within a 95% probability threshold.

Results

The δ_RMS and δ95% of Radixact^® Synchrony^® treatment for Gen1 kV systems deteriorated as the density of the target insert decreased, from 1.673 ± 0.064 mm and 3.049 ± 0.089 mm, respectively, for the 1.093 g/cc insert to 8.355 ± 5.873 mm and 15.297 ± 10.470 mm, respectively, for the 0.280 g/cc insert. In contrast, no such trend was observed in the δ_RMS or δ95% of Synchrony^® treatment using the Gen2 kV system. The δ_RMS and δ95%, respectively, fluctuated slightly from 1.586 to 1.687 mm and from 2.874 to 2.971 mm when different target inserts were tracked by the Gen2 kV system.

Conclusion

With improved image contrast in kV radiographs, the Gen2 kV imaging system can enhance the ability to track targets accurately in Radixact^® Lung Synchrony^® treatment and reduce the segmentation error. Our study showed that lung targets with density values as low as 0.280 cc/g could be tracked correctly in Synchrony treatment with the Gen2 kV imaging system.

Performance evaluation of online motion synchronizing systems for patient-specific respiratory movements with helical tomotherapy

BACKGROUND AND PURPOSE

The introduction of online motion synchronizing system on helical tomotherapy paves way for robust motion tracking. A recent upgrade launches modifications on both hardware and software of the kV tracking system. An evaluation on the kV subsystems, prior(Version1) and post upgrade(Version2), was performed to compare tracking accuracy by means of fiducial tracking error and resulted root-mean-square (RMS). Impacts influenced by various patient-specific breathing pattern regularities and target movements were also investigated to refine motion tracking error estimations upon future selection of possible candidates.

MATERIALS AND METHODS

Respiratory patterns from twenty-five lung cases were imported individually into a commercial dynamic platform model. Situating a phantom implanted with gold fiducial markers on the platform, superior-inferior(SI) movements of corresponding targets were simulated. Each case was delivered via an identical treatment plan in Version1 and was repeated in Version2. Motion tracking accuracy, by means of discrepancies between subsystem predicted model and raw data motion recorded in patient CT simulation, was analyzed statistically. Wilcoxon signed ranked test was employed to evaluate the difference in tracking error range between the two versions. Statistical model was fitted to inspect the dependence of internal target movement towards fiducial tracking errors.

RESULTS

A small difference of ±1 mm was exhibited in 99% of fiducial tracking errors for all cases experimented under both versions. RMS errors were all below 0.5 mm. Version2 demonstrated a greater extremity in fiducial tracking error (p=0.04). A positive correlation was depicted between internal target amplitudes and 95% interval of fiducial tracking errors (p< 0.02). Overall, irregular respiratory patterns tended to have greater fiducial tracking errors.

CONCLUSIONS

The excellent tracking performance in both kV subsystem versions offers motion compensations benefits, yet Version1 outperformed Version2 in fiducial tracking accuracy. It is noticeable that greater magnitude in internal target movement and irregular breathing patterns yield greater tracking error.

Clinical implementation of kVCT-guided tomotherapy with ClearRT

A helical fan-beam kilovoltage computed tomography (kVCT) was recently introduced into Tomotherapy units. This study aims to share the initial experience of kVCT in clinical workflow, compare its performance with that of the existing megavoltage computed tomography (MVCT), and explore its potential in adaptive planning. We retrospectively enrolled 23 patients who underwent both MVCT and kVCT scans. The clinical performance data regarding image acquisition time, nominal dose length product (DLP), registration time and registration corrections were extracted and compared. Image quality was scored by six experienced radiation therapists and quantified based on phantom measurements. CT number stability and the implementation of adaptive radiotherapy were dosimetrically evaluated by performing the dose recalculation on kVCT. Compared to MVCT, kVCT significantly reduced DLP (except the highest kVp protocol), image acquisition and registration time. KVCT obtained higher scores than MVCT on all criteria except artifacts. Phantom measurements also revealed a better image performance characterization of kVCT except for image uniformity. The CT number variation could lead to a dose difference of 0.5% for D95% of target and D_mean of organ-at-risk. For the treatment planning with kVCT, a systematic dose difference (> 1%) in PTV dose metrics was observed at regions with large longitudinal density discontinuities compared to the reference plans. The new kVCT imaging provides enhanced soft-tissue visualization. The improved efficiency with kVCT-guided treatment will allow more patients to be treated each day. In most cases, the dose calculation accuracy of kVCT images is acceptable except for regions with severe artifacts.

Dosimetric impact of phase shifts on Radixact Synchrony tracking system with patient-specific breathing patterns

Purpose

The Synchrony tracking system of Radixact is capable of real‐time tumor tracking by building a correlation model between external light‐emitting diodes on the patient's chest and an internal marker. A phase shift between the chest wall and a lung tumor has been reported. Hence, this study focused on evaluating the accuracy of the tracking system, especially under a patient‐specific breathing pattern with respiratory phase shifts.

Methods

A phantom containing fiducial markers was placed on a moving platform. The intrinsic delivery accuracy was verified with a patient‐specific breathing pattern. Three patient‐specific breathing patterns were then implemented, for which phase shifts, φ, were introduced. Phase shifts with +0.3 s and +1 s were tested for dosimetric aspects, whereas ±0.3, ±0.6, and ±0.8 s shifts were used for tracking accuracy. The resultant dose distributions were analyzed by γ comparison. Dose profiles in the superior‐inferior and lateral directions were compared. Logfiles of the tracking information were extracted from the system and compared with the input breathing pattern. The root mean square (RMS) difference was used to quantify the consistency.

Results

When the φ value was as large as 1 s, a severe inconsistency was observed. The target was significantly underdosed, down to 89% of the originally planned dose. γ analysis revealed that the failed portion was concentrated in the target region. The RMS of the tracking difference was close to 1 mm when φ was ±0.3 s and approximately 4 mm when φ was ±0.8 s. Tracking errors increased with an increase in the degree of phase shifts.

Conclusion

Phase shifts between the patient chest wall and the internal target may hamper treatment delivery and jeopardize treatment using Synchrony Tracking. Hence, a larger planning target volume (PTV) may be necessary if a large phase shift is observed in a patient, especially when an external surrogate shows a lag in motion when compared with the tumor.

Acquisition repeatability of MRI radiomics features in the head and neck: a dual-3D-sequence multi-scan study

Radiomics has increasingly been investigated as a potential biomarker in quantitative imaging to facilitate personalized diagnosis and treatment of head and neck cancer (HNC), a group of malignancies associated with high heterogeneity. However, the feature reliability of radiomics is a major obstacle to its broad validity and generality in application to the highly heterogeneous head and neck (HN) tissues. In particular, feature repeatability of radiomics in magnetic resonance imaging (MRI) acquisition, which is considered a crucial confounding factor of radiomics feature reliability, is still sparsely investigated. This study prospectively investigated the acquisition repeatability of 93 MRI radiomics features in ten HN tissues of 15 healthy volunteers, aiming for potential magnetic resonance-guided radiotherapy (MRgRT) treatment of HNC. Each subject underwent four MRI acquisitions with MRgRT treatment position and immobilization using two pulse sequences of 3D T1-weighed turbo spin-echo and 3D T2-weighed turbo spin-echo on a 1.5 T MRI simulator. The repeatability of radiomics feature acquisition was evaluated in terms of the intraclass correlation coefficient (ICC), whereas within-subject acquisition variability was evaluated in terms of the coefficient of variation (CV). The results showed that MRI radiomics features exhibited heterogeneous acquisition variability and uncertainty dependent on feature types, tissues, and pulse sequences. Only a small fraction of features showed excellent acquisition repeatability (ICC > 0.9) and low within-subject variability. Multiple MRI scans improved the accuracy and confidence of the identification of reliable features concerning MRI acquisition compared to simple test-retest repeated scans. This study contributes to the literature on the reliability of radiomics features with respect to MRI acquisition and the selection of reliable radiomics features for use in modeling in future HNC MRgRT applications.

A narrative review of MRI acquisition for MR-guided-radiotherapy in prostate cancer

Magnetic resonance guided radiotherapy (MRgRT), enabled by the clinical introduction of the integrated MRI and linear accelerator (MR-LINAC), is a novel technique for prostate cancer (PCa) treatment, promising to further improve clinical outcome and reduce toxicity. The role of prostate MRI has been greatly expanded from the traditional PCa diagnosis to also PCa screening, treatment and surveillance. Diagnostic prostate MRI has been relatively familiar in the community, particularly with the development of Prostate Imaging - Reporting and Data System (PI-RADS). But, on the other hand, the use of MRI in the emerging clinical practice of PCa MRgRT, which is substantially different from that in PCa diagnosis, has been so far sparsely presented in the medical literature. This review attempts to give a comprehensive overview of MRI acquisition techniques currently used in the clinical workflows of PCa MRgRT, from treatment planning to online treatment guidance, in order to promote MRI practice and research for PCa MRgRT. In particular, the major differences in the MRI acquisition of PCa MRgRT from that of diagnostic prostate MRI are demonstrated and explained. Limitations in the current MRI acquisition for PCa MRgRT are analyzed. The future developments of MRI in the PCa MRgRT are also discussed.

Phantom assessment of three-dimensional geometric distortion of a dedicated wide-bore MR-simulator for radiotherapy

This study evaluated the machine-dependent three-dimensional geometric distortion images acquired from a 1.5T 700mm-wide bore MR-simulator based on a large geometric accuracy phantom. With the consideration of radiation therapy (RT) application requirements, every sequence was examined in various combinations of acquisition-orientations and receiver-bandwidths with console-integrated distortion correction enabled. Distortion was repeatedly measured over a six-month period. The distortion measured from the images acquired at the beginning of this period was employed to retrospectively correct the distortion in the subsequent acquisitions. Geometric distortion was analyzed within the largest field-of-view allowed. Six sequences were examined for comprehensive distortion analysis - VIBE, SPACE, TSE, FLASH, BLADE and PETRA. Based on optimal acquisition parameters, their diameter-sphere-volumes (DSVs) of CT-comparable geometric fidelity (where 1mm distortion was allowed) were 333.6mm, 315.1mm, 316.0mm, 318.9mm, 306.2mm and 314.5mm respectively. This was a significant increase from 254.0mm, 245.5mm, 228.9mm, 256.6mm, 230.8mm and 254.2mm DSVs respectively, when images were acquired using un-optimized parameters. The longitudinal stability of geometric distortion and the efficacy of retrospective correction of console-corrected images, based on prior distortion measurements, were inspected using VIBE and SPACE. The retrospectively corrected images achieved over 500mm DSVs with 1mm distortion allowed. The median distortion was below 1mm after retrospective correction, proving that obtaining prior distortion map for subsequent retrospective distortion correction is beneficial. The systematic evaluation of distortion using various combinations of sequence-type, acquisition-orientation and receiver-bandwidth in a six-month time span would be a valuable guideline for optimizing sequence for various RT applications.

Out-of-field dose and its constituent components for a 1.5 T MR-Linac

This study aims to quantify the relative contributions of phantom scatter, collimator scatter and head leakage to the out-of-field doses (OFDs) of both static fields and clinical intensity-modulated radiation therapy (IMRT) treatments in a 1.5 T MR-Linac. The OFDs of static fields were measured at increasing distances from the field edge in an MR-conditional water phantom. Inline scans at depths of dmax (14 mm), 50 and 100 mm were performed for static fields of 5 × 5, 10 × 10 and 15 × 15 cm²under three different conditions: full scatter, with phantom scatter prevented, and head leakage only. Crossline scans at isocenter and offset positions were performed in full scatter condition. EBT3 radiochromic films were placed at 100 mm depth of solid water phantom to measure the OFD of clinical IMRT plans. All water tank data were normalized to Dmax of a 10 × 10 cm²field and the film results were presented as a fraction of the target mean dose.The OFD in the inline direction varied from 3.5% (15 × 15 cm², 100 mm depth, 50 mm distance) to 0.014% (5 × 5 cm², dmax, 400 mm distance). For all static fields, the collimator scatter was higher than the phantom scatter and head leakage at a distance of 100-400 mm. Head leakage remained the smallest among the three components, except at long distances (>375 mm) with small field size. Compared to the inline scans, the crossline scans at the isocenter showed higher doses at distances longer than 80 mm. All crossline profiles at longitudinal offset positions showed a cone shape with laterally shifted maxima. The OFD of IMRT deliveries varied with different target size. For prostate stereotactic body radiation therapy (SBRT) treatment, the OFD decreased from 2% to 0.03% at a distance of 50-500 mm. The OFDs have been measured for a 1.5 T MR-Linac. The presented dosimetric data are valuable for radiation safety assessments on patients treated with the MR-Linac, such as evaluating carcinogenic risk and radiation exposure to cardiac implantable electronic devices.

The association between central ambulatory blood pressure and end-organ damage

Background/Introduction

Central blood pressure (BP) measured in the clinic relates more strongly to end-organ damage and may be superior to brachial BP in predicting cardiovascular events. Measurement of central ambulatory BP (ABP) is now possible and emerging data suggest that central ABP is better correlated with left ventricular mass index (LVMI) than brachial ABP.1,2 However, the association between central ABP and other measures of end-organ damage remains unclear.

Purpose

We investigated the association between brachial or central ABP and end-organ damage, in a large community-based population of untreated individuals.

Method

1091 participants (mean age 45±18 years; 589 females) had simultaneous measurements of brachial and central ABP over 24 hours, using the Mobil-O-Graph device. Central ABP was derived using two waveform calibration methods (1: mean/diastolic BP; 2: systolic/diastolic BP). Participants also underwent measurement of aortic pulse wave velocity (aPWV; SphygmoCor device) in the clinic. In a subset of 675 individuals, LVMI was assessed by echocardiography and in 610 individuals, carotid intima-media thickness (CIMT) was measured, using ultrasound. 24-hour and daytime brachial and central ambulatory systolic BP (ASBP) and pulse pressure (APP) were considered.

Results

The Pearson's coefficient for each correlation is listed in Table 1. LVMI was most strongly correlated with 24-hour central ASBP, using calibration method 1 (MAP/DBP; r=0.403, P&lt;0.001). CIMT was most strongly correlated with daytime central ASBP, again derived from calibration method 1 (r=0.341, p&lt;0.001), whereas aPWV measured in the clinic, was most strongly correlated with 24-hour central ASBP, derived from calibration method 2 (SBP/DBP; r=0.441, P&lt;0.001). Based on z statistics, all correlations reported above were significantly stronger than equivalent correlations using brachial ASBP (P&lt;0.001 for all three comparisons of correlation coefficients).

Conclusion

Measurement of central ABP relates more closely to end-organ damage than equivalent measures based on brachial ABP. Central ABP may provide valuable additional information concerning cardiovascular risk above and beyond brachial ABP.

Funding Acknowledgement

Type of funding sources: Public Institution(s). Main funding source(s): British Heart Foundation

Reliability of MRI radiomics features in MR-guided radiotherapy for prostate cancer: Repeatability, reproducibility, and within-subject agreement

PURPOSE

The MR-guided radiotherapy (MRgRT) images on the integrated MRI and linear accelerator (MR-LINAC) might facilitate radiomics analysis for longitudinal treatment response assessment. However, the reliability of MRgRT radiomics features is largely unknown. This study aims to investigate MRgRT radiomics feature reliability acquired using a standardized 3D-T2W-TSE sequence in terms of repeatability, reproducibility, and within-subject feature agreement on a 1.5T MR-simulator and a 1.5T MR-LINAC for prostate cancer (PC).

METHODS

Twenty-six consecutive PC patients who underwent one MRI-simulator scan and two MR-LINAC scans before dose delivery were retrospectively included. The three MRI datasets were rigidly co-registered. 1023 first-order and texture radiomics features were extracted with different intensity bin widths for each scan in the manually segmented clinical target volume (CTV) and planning target volume (PTV) by an experienced radiation oncologist. Intraclass correlation coefficient (ICC) was used to evaluate feature repeatability between MR-LINAC scans and reproducibility between MRI-simulator and MR-LINAC scans. The within-subject feature value agreements were evaluated using Bland-Altman analysis. The impact of inter-observer segmentation on the radiomics feature reliability was also examined based on the second manual segmentation of CTV and PTV by an MRI researcher.

RESULTS

Based on the segmentation by the radiation oncologist and the default bin width of 25, 9.6%, 24.1%, 49.6%, and 16.8% of the total 1023 features exhibited excellent (ICC > 0.9), good (0.9 > ICC > 0.75), moderate (0.75 > ICC > 0.5), and poor (ICC < 0.5) repeatability in the CTV, and 9.2%, 26.8%, 50.5%, and 13.5% in the PTV, respectively. For reproducibility, the corresponding feature percentages were 8.9%, 19.7%, 41.9%, and 29.6% in the CTV, and 8.4%, 17.8%, 47.9%, and 26% in the PTV. Feature reliability was not notably influenced by intensity bin width for discretization. BA analysis revealed wide 95% limit-of-agreements and substantial biases of feature values between CTV and PTV and between any two MRI scans. The features even with excellent ICC were still subjected to considerable inter-scan feature variations in each individual subject. The analysis on the second segmentation by the MRI researcher showed insignificantly different feature repeatability and reproducibility in terms of ICC values.

CONCLUSIONS

Only a small proportion of features exhibited excellent/good repeatability and reproducibility, highlighting the importance of reliable MRgRT feature selection. The within-subject feature values were subjected to considerable inter-scan variations, imposing a challenge on the determination of the smallest detectable change in future MRgRT delta-radiomics studies.

Magnetic Resonance-Guided Radiation Therapy of Patients With Cardiovascular Implantable Electronic Device on a 1.5 T Magnetic Resonance-Linac

Magnetic resonance-guided radiation therapy is reported for treating patients with an insertable cardiac monitor and implantable cardiac pacemakers. All treatments were delivered using a 1.5 T MR-Linac. Among the 4 patients, 2 were treated with stereotactic body radiation therapy at a dose of 40 Gy in 5 fractions. A clinical safety protocol was developed to address the decision-making and patient selection, as well as the clarified responsibilities of different parties for management of patients with cardiovascular implantable electronic devices. Dose estimation based on out-of-field dose data are necessary for cardiovascular implantable electronic devices located outside the treatment fields.

Comparison of modeling accuracy between Radixact®and CyberKnife®Synchrony®respiratory tracking system

Synchrony Respiratory Tracking system adapted from CyberKnife has been introduced in Radixact to compensate the tumor motion caused by respiration. This study aims to compare the modeling accuracy of the Synchrony system between Radixact and CyberKnife. Two Synchrony plans based on fiducial phantoms were created for CyberKnife and Radixact, respectively. Different respiratory motion traces were used to drive a motion platform to move along the superoinferior and left-right direction. The cycle time and the amplitude of target/surrogate motion of one selected motion trace were scaled to investigate the dependence of modeling accuracy on the motion characteristic. The predicted target position, the correlation error, potential difference (Radixact only) and standard error (CyberKnife only) were extracted from raw data or log files of the two systems. The modeling accuracy was evaluated by calculating the root-mean-square (RMS) error between the predicted target positions and the input motion trace. A threshold T95 within which 95% of the potential difference or the standard error lay was defined and evaluated. Except for the motion trace with a small amplitude and a good (linear) correlation between target and surrogate motion, Radixact showed smaller RMS errors than CyberKnife. The RMS error of both systems increased with the motion amplitude and showed a decreasing trend with the increasing cycle time. No correlation was found between the RMS error and the amplitude of surrogate motion. T95 could be a good estimator of modeling accuracy for CyberKnife rather than Radixact. The correlation error defined in Radixact were largely affected by the number of fiducial markers and the setup error. In general, the modeling accuracy of the Radixact Synchrony system is better than that of the CyberKnife Synchrony system under unfavorable conditions.

3D T1-weighted turbo spin echo contrast-enhanced MRI at 1.5 T for frameless brain metastases radiotherapy

PURPOSE

Performance of 3D-T1W-TSE has been proven superior to 3D-MP-GRE at 3 T on brain metastases (BM) contrast-enhanced (CE) MRI. However, its performance at 1.5 T is largely unknown and sparsely reported. This study aims to assess image quality, lesion detectability and conspicuity of 1.5 T 3D-T1W-TSE on planning MRI of frameless BM radiotherapy.

METHODS

94 BM patients to be treated by frameless brain radiotherapy were scanned using 3D-T1W-TSE with immobilization on multi-vendor 1.5 T MRI-simulators. BMs were jointly diagnosed by 4 reviewers. Enhanced lesion conspicuity was quantitatively assessed by calculating contrast ratio (CR) and contrast-to-noise ratio (CNR). Signal-to-noise ratio (SNR) reduction of white matter due to the use of flexible coil was assessed. Lesion detectability and conspicuity were compared between 1.5 T planning MRI and 3 T diagnostic MRI by an oncologist and a radiologist in 10 patients.

RESULTS

497 BMs were jointly diagnosed. The CR and CNR were 75.2 ± 39.9% and 14.2 ± 8.1, respectively. SNR reduced considerably from 31.7 ± 8.3 to 21.9 ± 5.4 with the longer distance to coils. 3 T diagnostic MRI and 1.5 T planning MRI yielded exactly the same detection of 84 BMs. Qualitatively, lesion conspicuity at 1.5 T was not inferior to that at 3 T. Quantitatively, lower brain SNR and lesion CNR were found at 1.5 T, while lesion CR at 1.5 T was highly comparable to that at 3 T.

CONCLUSION

1.5 T 3D-T1W-TSE planning MRI of frameless BM radiotherapy was comprehensively assessed. Highly comparable BM detectability and conspicuity were achieved by 1.5 T planning MRI compared to 3 T diagnostic MRI. 1.5 T 3D-T1W-TSE should be valuable for frameless brain radiotherapy planning.

Quantitative assessment of acquisition imaging parameters on MRI radiomics features: a prospective anthropomorphic phantom study using a 3D-T2W-TSE sequence for MR-guided-radiotherapy

Background

MRI pulse sequences and imaging parameters substantially influence the variation of MRI radiomics features, thus impose a critical challenge on MRI radiomics reproducibility and reliability. This study aims to prospectively investigate the impact of various imaging parameters on MRI radiomics features in a 3D T2-weighted (T2W) turbo-spin-echo (TSE) pulse sequence for MR-guided-radiotherapy (MRgRT).

Methods

An anthropomorphic phantom was scanned using a 3D-T2W-TSE MRgRT sequence at 1.5T under a variety of acquisition imaging parameter changes. T1 and T2 relaxation times of the phantom were also measured. 93 first-order and texture radiomics features in the original and 14 transformed images, yielding 1,395 features in total, were extracted from 10 volumes-of-interest (VOIs). The percentage deviation (d%) of radiomics feature values from the baseline values and intra-class correlation coefficient (ICC) with the baseline were calculated. Robust radiomics features were identified based on the excellent agreement of radiomics feature values with the baseline, i.e., the averaged d% <5% and ICC >0.90 in all VOIs for all imaging parameter variations.

Results

The radiomics feature values changed considerably but to different degrees with different imaging parameter adjustments, in the ten VOIs. The deviation d% ranged from 0.02% to 321.3%, with a mean of 12.5% averaged for all original features in all ten VOIs. First-order and GLCM features were generally more robust to imaging parameters than other features in the original images. There were also significantly different radiomics feature values (ANOVA, P<0.001) between the original and the transformed images, exhibiting quite different robustness to imaging parameters. 330 out of 1395 features (23.7%) robust to imaging parameters were identified. GLCM and GLSZM features had the most (42.5%, 153/360) and least (3.8%, 9/240) robust features in the original and transformed images, respectively.

Conclusions

This study helps better understand the quantitative dependence of radiomics feature values on imaging parameters in a 3D-T2W-TSE sequence for MRgRT. Imaging parameter heterogeneity should be considered as a significant source of radiomics variability and uncertainty, which must be well harmonized for reliable clinical use. The identified robust features to imaging parameters are helpful for the pre-selection of radiomics features for reliable radiomics modeling.

Initial clinical experience of patient-specific QA of treatment delivery in online adaptive radiotherapy using a 1.5 T MR-Linac

Purpose. This study aims to evaluate the performance of a commercial 1.5 T MR-Linac by analyzing its patient-specific quality assurance (QA) data collected during one full year of clinical operation.Methods and Materials. The patient-specific QA system consisted of offline delivery QA (DQA) and online calculation-based QA. Offline DQA was based on ArcCHECK-MR combined with an ionization chamber. Online QA was performed using RadCalc that calculated and compared the point dose calculation with the treatment planning system (TPS). A total of 24 patients with 189 treatment fractions were enrolled in this study. Gamma analysis was performed and the threshold that encompassed 95% of QA results (T95) was reported. The plan complexity metric was calculated for each plan and compared with the dose measurements to determine whether any correlation existed.Results. All point dose measurements were within 5% deviation. The mean gamma passing rates of the group data were found to be 96.8 ± 4.0% and 99.6 ± 0.7% with criteria of 2%/2mm and 3%/3mm, respectively. T95 of 87.4% and 98.2% was reported for the overall group with the two passing criteria, respectively. No statistically significant difference was found between adaptive treatments with adapt-to-position (ATP) and adapt-to-shape (ATS), whilst the category of pelvis data showed a better passing rate than other sites. Online QA gave a mean deviation of 0.2 ± 2.2%. The plan complexity metric was positively correlated with the mean dose difference whilst the complexity of the ATS cohort had larger variations than the ATP cohort.Conclusions. A patient-specific QA system based on ArcCHECK-MR, solid phantom and ionization chamber has been well established and implemented for validation of treatment delivery of a 1.5 T MR-Linac. Our QA data obtained over one year confirms that good agreement between TPS calculation and treatment delivery was achieved.

Effects on skin dose from unwanted air gaps under bolus in an MR-guided linear accelerator (MR-linac) system

Bolus is commonly used in MV photon radiotherapy to increase superficial dose and improve dose uniformity for treating shallow lesions. However, irregular patient body contours can cause unwanted air gaps between a bolus and patient skin. The resulting dosimetric errors could be exacerbated in MR-Linac treatments, as secondary electrons generated by photons are affected by the magnetic field. This study aimed to quantify the dosimetric effect of unwanted gaps between bolus and skin surface in an MR-Linac. A parallel-plate ionization chamber and EBT3 films were utilized to evaluate the surface dose under bolus with various gantry angles, field sizes, and different air gaps. The results of surface dose measurements were then compared to Monaco 5.40 Treatment Planning System (TPS) calculations. The suitability of using a parallel-plate chamber in MR-Linac measurement was validated by benchmarking the percentage depth dose and output factors with the microDiamond detector and air-filled ionization chamber measurements in water. A non-symmetric response of the parallel-plate chamber to oblique beams in the magnetic field was characterized. Unwanted air gaps significantly reduced the skin dose. For a frontal beam, skin dose was halved when there was a 5 mm gap, a much larger difference than in a conventional linac. Skin dose manifested a non-symmetric pattern in terms of gantry angle and gap size. The TPS overestimated skin dose in general, but shared the same trend with measurement when there was no air gap, or the gap size was larger than 5 mm. However, the calculated and measured results had a large discrepancy when the bolus-skin gap was below 5 mm. When treating superficial lesions, unwanted air gaps under the bolus will compromise the dosimetric goals. Our results highlight the importance of avoiding air gaps between bolus and skin when treating superficial lesions using an MR-Linac system.

Longitudinal acquisition repeatability of MRI radiomics features: An ACR MRI phantom study on two MRI scanners using a 3D T1W TSE sequence

PURPOSE

The purpose of this study was to quantitatively assess the longitudinal acquisition repeatability of MRI radiomics features in a three-dimensional (3D) T1-weighted (T1W) TSE sequence via a well-controlled prospective phantom study.

METHODS

Thirty consecutive daily datasets of an ACR-MRI phantom were acquired on two 1.5T MRI simulators using a 3D T1W TSE sequence. Images were blindly segmented by two observers. Post-acquisition processing was minimized but an intensity discretization (fixed bin size of 25). One hundred and one radiomics features (shape n = 12; first order n = 16; texture n = 73) were extracted. Longitudinal repeatability of each feature was evaluated by Pearson correlation and coefficient of variance (CV_68% ). Interobserver feature value agreement was also quantified using intraclass correlation coefficient (ICC) and Bland-Altman analysis. A most repeatable radiomics feature set on both scanners was determined by feature coefficient of variance (CV_68% <5%), ICC (>0.75), and the ratio of the interobserver difference to the interobserver mean δ<5%.

RESULTS

No trend of radiomics feature value changed with time. Longitudinal feature repeatability CV_68% ranged 0.01-38.60% (mean/median: 12.5%/9.9%), and 0.01-40.47%, (8.49%/7.34%) on the scanners A and B. Shape features exhibited significantly better repeatability than first-order and texture features (all P < 0.01). Significant longitudinal repeatability difference was observed in texture features (P < 0.001) between the two scanners, but not in shape and first-order features (P > 0.30). First-order and texture features had smaller interobserver-dependent variation than acquisition-dependent variation. They also showed good interobserver agreement on both scanners (A:ICC = 0.80 ± 0.23; B:ICC = 0.80 ± 0.22), independent of acquisition repeatability. The repeatable radiomics features in common on both scanners, including 12 shape features, 0 first-order features, and 3 texture features, were determined as the most repeatable MRI radiomics feature set.

CONCLUSIONS

Radiomics features exhibited heterogeneous longitudinal repeatability, while the shape features were the most repeatable, in this phantom study with a 3D T1W TSE acquisition. The most repeatable radiomics feature set derived in this study should be helpful for the selection of reliable radiomics features in the future clinical use.

Reliability of radiomics features due to image reconstruction using a standardized T2 -weighted pulse sequence for MR-guided radiotherapy: An anthropomorphic phantom study

PURPOSE

To prospectively investigate the impact of image reconstruction on MRI radiomics features.

METHODS

An anthropomorphic phantom was scanned at 1.5 T using a standardized sequence for MR-guided radiotherapy under SENSE and compressed-SENSE reconstruction settings. A total of 93 first-order and texture radiomics features in 10 volumes of interest were assessed based on (1) accuracy measured by the percentage deviation from the reference, (2) robustness on reconstruction in all volumes of interest measured by the intraclass correlation coefficient, and (3) repeatability measured by the coefficient of variance over the repetitive acquisitions. Finally, reliable and unreliable radiomics features were comprehensively determined based on their accuracy, robustness, and repeatability.

RESULTS

Better accuracy and robustness of the radiomics features were achieved under SENSE than compressed-SENSE reconstruction. The feature accuracy under SENSE reconstruction was more affected by acceleration factor than direction, whereas under compressed-SENSE reconstruction, accuracy was substantially impacted by the increasing denoising levels. Feature repeatability was dependent more on feature types than on reconstruction. A total of 45 reliable features and 13 unreliable features were finally determined for SENSE, compared with 22 reliable and 26 unreliable features for compressed SENSE. First-order and gray-level co-occurrence matrix features were generally more reliable than other features.

CONCLUSION

Radiomics features could be substantially affected by MRI reconstruction, so precautions need to be taken regarding their reliability for clinical use. This study helps the guidance of the preselection of reliable radiomics features and the preclusion of unreliable features in MR-guided radiotherapy.

Deep Learning-Based Decision-Tree Classifier for COVID-19 Diagnosis From Chest X-ray Imaging

The global pandemic of coronavirus disease 2019 (COVID-19) has resulted in an increased demand for testing, diagnosis, and treatment. Reverse transcription polymerase chain reaction (RT-PCR) is the definitive test for the diagnosis of COVID-19; however, chest X-ray radiography (CXR) is a fast, effective, and affordable test that identifies the possible COVID-19-related pneumonia. This study investigates the feasibility of using a deep learning-based decision-tree classifier for detecting COVID-19 from CXR images. The proposed classifier comprises three binary decision trees, each trained by a deep learning model with convolution neural network based on the PyTorch frame. The first decision tree classifies the CXR images as normal or abnormal. The second tree identifies the abnormal images that contain signs of tuberculosis, whereas the third does the same for COVID-19. The accuracies of the first and second decision trees are 98 and 80%, respectively, whereas the average accuracy of the third decision tree is 95%. The proposed deep learning-based decision-tree classifier may be used in pre-screening patients to conduct triage and fast-track decision making before RT-PCR results are available.

A fast volumetric 4D-MRI with sub-second frame rate for abdominal motion monitoring and characterization in MRI-guided radiotherapy

Background

To propose a fast volumetric 4D-MRI based on 3D pulse sequence acquisition for abdominal motion monitoring and characterization in MRI-guided radiotherapy (MRgRT).

Methods

A 3D spoiled gradient echo sequence volumetric interpolated breath-hold examination (VIBE) [repetition time/echo time (TR/TE) =0.53/1.57 ms, flip-angle =5°, receiver bandwidth (RBW) =1,400 Hz/voxel] based 4D-MRI acquisition, accelerated by 4-fold controlled aliasing in parallel imaging results in higher acceleration (CAIPIRINHA), named CAIPIRINHA-VIBE 4D-MRI, was implemented on a 1.5T MRI simulator (MR-sim) and applied for abdominal imaging of nine healthy volunteers under free breathing. One hundred and forty-four dynamics of the entire abdomen volume (56 slices), in total 8,064 (144×56) images with a voxel size of 2.7×2.7×4.0 mm³, were acquired in 89 s for 4D-MRI. This CAIPIRINHA-VIBE 4D-MRI was qualitatively compared with a 2D half-Fourier acquisition single-shot turbo spin-echo (2D-HASTE) based 4D-MRI. The motions of liver dome, kidney and spleen were analyzed using the CAIPIRINHA-VIBE 4D-MRI data. The kidney motion was quantitatively characterized in terms of motion range and the correlations between left and right kidneys.

Results

CAIPIRINHA-VIBE 4D-MRI was successfully conducted in all subjects. CAIPIRINHA-VIBE 4D-MRI exhibited much higher effective volumetric temporal resolution (0.615 vs. ~5 s/volume) and better reconstructed volume consistency than 2D-HASTE 4D-MRI. CAIPIRINHA-VIBE 4D-MRI was able to characterize the respiratory motion of abdominal organs simultaneously in three orthogonal directions, and could potentially be used for whole abdomen deformable motion tracking. Renal motion range was most pronounced in superior-inferior (SI) direction (L: 10.03±2.65 mm; R: 10.38±2.80 mm), significantly larger (P<0.001) than that in anterior-posterior (AP) and the least in left-right (LR) directions. Right kidney had significantly larger mobility (4.18±2.19 vs. 2.32±1.34 mm, P=0.045) than left kidney in AP, but not in LR and SI directions. The Pearson correlation coefficients r between left and right kidney motion were 0.5063 (P=0.164), 0.6624 (P=0.052) and 0.5752 (P=0.105) in LR, AP and SI correspondingly. The correlation of renal motion in SI and AP was found significant in right kidney (r=0.843, P=0.004) but not in left kidney (r=0.467, P=0.205).

Conclusions

A fast volumetric 4D-MRI was implemented for abdominal motion monitoring in MRgRT. A sub-second volumetric temporal resolution of 0.615 s, covering the entire abdomen, was demonstrated for respiratory motion monitoring and characterization. This technique holds potentials for MRgRT applications.

A pilot study of highly accelerated 3D MRI in the head and neck position verification for MR-guided radiotherapy

Background

To evaluate the performance of a highly accelerated 3D MRI on inter-fractional positional measurement for MR-guided radiotherapy (MRgRT) in the head and neck (HN).

Methods

Fourteen healthy volunteers received 159 scans on a 1.5 T MR-sim to simulate MRgRT fractions. MRI acquisition included a high-resolution (HQI-MRI, voxel-size =1.05×1.05×1.05 mm³, duration =5 min) and a highly-accelerated low-resolution (true-LQI-MRI, acceleration-factor =9, voxel-size =1.4×1.4×1.4 mm³, duration =86 s) T1w spin-echo sequence (TR/TE =420/7.2 ms). The first session HQI-MRI was used as the reference to mimic planning MRI. Other HQI-MRI was also retrospectively down-sampled in K-space and GRAPPA reconstructed to generate pseudo-LQI-MRI. Inter-sessional positional shift calculated from HQI-MRI, true-LQI-MRI and pseudo-LQI-MRI rigidly registering to the reference were analyzed and compared in the overall HN and the sub-regions of brain, nasopharynx, oropharynx and hypopharynx.

Results

The calculated SD of systematic errors (Σ) from HQI-MRI/pseudo-LQI-MRI/true-LQI-MRI images for overall HN were 1.11/1.14/1.08, 0.28/0.26/0.29, 0.43/0.44/0.60, and 0.77/0.79/0.74 mm for translation in LR, AP, SI and 3D, respectively; The corresponding RMS of random errors (σ) were 0.97/0.98/0.96, 0.28/0.27/0.26, 0.77/0.77/0.72, and 0.85/0.87/0.85 mm. For all sub-regions, brain showed the smallest Σ and σ in 3D. Other sub-regions showed direction-dependent error patterns, but the positioning results were consistent, independent of the datasets used for registration.

Conclusions

A highly-accelerated 3D-MRI could be used for MR-guided HN radiotherapy without compromising position verification accuracy.

A novel method for monitoring the constancy of beam path accuracy in CyberKnife

The aim of current work was to present a novel evaluation procedure implemented for checking the constancy of beam path accuracy of a CyberKnife system based on ArcCHECK. A tailor‐made Styrofoam with four implanted fiducial markers was adopted to enable the fiducial tracking during beam deliveries. A simple two‐field plan and an isocentric plan were created for determining the density override of ArcCHECK in MultiPlan and the constancy of beam path accuracy respectively. Correlation curves for all diodes involved in the study were obtained by analyzing the dose distributions calculated by MultiPlan after introducing position shifts in anteroposterior, superoinferior, and left–right directions. The ability of detecting systematic position error was also evaluated by changing the position of alignment center intentionally. The one standard deviation (SD) result for reproducibility test showed the RMS of 0.054 mm and the maximum of 0.263 mm, which was comparable to the machine self‐test result. The mean of absolute value of position errors in the constancy test was measured to 0.091 mm with a SD of 0.035 mm, while the root‐mean‐square was 0.127 mm with a SD of 0.034 mm. All introduced systematic position errors range from 0.3 to 2 mm were detected successfully. Efficient method for evaluating the constancy of beam path accuracy of CyberKnife has been developed and proven to be sensitive enough for detecting a systematic drift of robotic manipulator. Once the workflow is streamlined, our proposed method will be an effective and easy quality assurance procedure for medical physicists.

Performance evaluation of the CyberKnife system in real-time target tracking during beam delivery using a moving phantom coupled with two-dimensional detector array

The aim of the current study was to evaluate the tracking error of the Synchrony Respiratory Tracking system by conducting beam-by-beam analyses to determine the variation in the tracking beams measured during target motion. A moving phantom of in-house design coupled with a two-dimensional (2D) detector array was used to simulate respiratory motion in the superoinferior (SI) and anteroposterior (AP) direction. A styrofoam block with four implanted fiducial markers was placed on top of the detector to enable the fiducial-based respiratory tracking. Measurements were performed with the phantom under either stationary mode or sinusoidal motion of 6-s cycle and 15/20-mm amplitude at SI and AP direction. The measurement data were saved as movie files that were used to calculate the center shift of the beam with 100-ms sampling time. The tracking accuracy of the system was defined as the targeting error, which could be tracked with probability of > 95% (Ep95). The mean ± standard deviation of Ep95 was 0.28 ± 0.08 mm under stationary condition; 0.66 ± 0.23 mm (range: 0.28-1.22 mm) under sinusoidal respiratory motion. The maximum drift of the beam center for all beam paths was 2.7 mm. The tracking accuracy of CyberKnife Synchrony system was successfully evaluated using a moving phantom and 2D detector array; the maximum tracking error was < 1.5 mm for sinusoidal motion of amplitude ≤ 20 mm.

Assessment of positional reproducibility in the head and neck on a 1.5-T MR simulator for an offline MR-guided radiotherapy solution

Background

Recently, a shuttle-based offline magnetic resonance-guided radiotherapy (MRgRT) approach was proposed. This study aims to evaluate the positional reproducibility in the immobilized head and neck using a 1.5-T MR-simulator (MR-sim) on healthy volunteers.

Methods

A total of 159 scans of 14 healthy volunteers were conducted on a 1.5-T MR-sim with thermoplastic mask immobilization. MR images with isotropic 1.05³ mm³ voxel size were rigidly registered to the first scan based on fiducial, anatomical and gross positions. Mean and standard deviation of positional displacements in translation and rotation were assessed. Systematic error and random errors of positioning in the head and neck on the MR-sim were determined in the translation of, and in the rotation of roll, pitch and yaw.

Results

The systematic error (Σ) of translation in left-right (LR), anterior-posterior (AP) and superior-inferior (SI) direction was 0.57, 0.22 and 0.26 mm for fiducial displacement, 0.28, 0.10 and 0.52 mm for anatomical displacement, and 0.53, 0.22 and 0.49 mm for gross displacement, respectively. The random error (σ) in corresponding translation direction was 2.07, 0.54 and 1.32 mm for fiducial displacement, 1.34, 0.73 and 2.04 mm for anatomical displacement, and 2.24, 0.86 and 2.61 mm for gross displacement. The systematic error and random error of rotation were generally smaller than 1°.

Conclusions

Our results suggested that high gross positional reproducibility (<1 mm translational and <1° rotational systematic error) could be achieved on an MR-sim for the proposed offline MRgRT.

Development of a novel methodology for QA of respiratory-gated and VMAT beam delivery using Octavius 4D phantom

The objective of this study was to develop and evaluate a series of quality assurance (QA) techniques based on Octavius 4D phantom for testing of respiratory-gated treatment delivery, integrity of dose rate vs gantry speed in volumetric-modulated arc therapy (VMAT) commissioning, and multileaf collimator (MLC) positioning accuracy of a linear accelerator. An Octavius 4D phantom capable of rotating with the gantry and recording the detector signal with a sampling rate of 10 Hz was isocentrally set up and an inclinometer was also installed to measure the gantry angle simultaneously. A simple arc test was created and delivered with gating function activated to measure the timing accuracy of the gating window. A tailor-made dose rate vs gantry speed plan was also designed to test the accuracy of measured dose rate, gantry speed, and actual control points. All experiments were conducted while machine log files were collected for comparison. The variations of beam flatness, symmetry, and field size were analyzed as a function of gantry angle to evaluate the influence from the modulation of dose rate and gantry speed. MLC position accuracy was evaluated based on specific garden fence plans. The time of gating window was measured to be less than 10-millisecond deviation from the log data. Gantry backlash was observed and quantified to be 1.72° with an extra stabilization time of 1.16 seconds for a gating arc with gantry speed of 6°/s. In the dose rate vs gantry speed test, the mean deviation between measured gantry angle and log data was less than 0.2° after a time delay of 0.25 second was corrected. The measured dose rate agreed with the log data very well with a mean deviation of 0.05%, and even the transit of modulation was tracked successfully. There was a statistically significant difference on the variation of beam parameters between a VMAT plan and a simple arc plan. The induced MLC position errors were detected with an accuracy of 0.05 mm. The leaf position reproducibility was found to be better than 0.02 mm, whereas the routine MLC position accuracy was better than 0.1 mm. A time-resolved method using Octavius 4D phantom has been developed and proven to be convenient for respiratory gating QA, dose rate vs gantry speed test, and MLC QA. Gating time, dose rate, and gantry speed-induced leave position error could be directly measured with high accuracy after comparison with the machine log data. This study also highlights the capability of the phantom in quantifying the variation of flatness, symmetry, and field size during gantry rotation.

Correction to: Molecular characterization of Blastocystis sp. in captive wild animals in Qinling Mountains

There were errors in Fig. 1 of the originally published article. Correct fig. 1 is presented here.

Image quality assessment of a 1.5T dedicated magnetic resonance-simulator for radiotherapy with a flexible radio frequency coil setting using the standard American College of Radiology magnetic resonance imaging phantom test

BACKGROUND

A flexible RF coil setting has to be used on an MR-simulator (MR-sim) in the head and neck simulation scan for radiotherapy (RT) purpose, while the image quality might be compromised due to the sub-optimized flexible coil compared to the normal diagnostic radiological (DR) head coil. In this study, we assessed the image quality of an MR-sim by conducting the standard American College of Radiology (ACR) MRI phantom test on a 1.5T MR-sim under RT-setting and comparing it to DR-setting.

METHODS

A large ACR MRI phantom was carefully positioned, aligned and scanned 9 times for each under RT- and DR-setting on a 1.5T MR-sim, following the ACR scanning instruction. Images were analyzed following the ACR guidance. Measurement results under two coil settings were quantitatively compared. Inter-observer disagreements under RT-setting between two physicists were compared using Bland-Altman (BA) analysis and intra-class correlation coefficient (ICC).

RESULTS

The MR-sim with RT-setting obtained sufficiently good image quality to pass all ACR recommended criteria. No significant difference was found in phantom length accuracy, high-contrast spatial resolution, slice thickness accuracy, slice position accuracy, and percent-signal ghosting. RT-setting significantly under-performed in low-contrast object detectability, while better performed in image intensity uniformity. BA analysis showed that 95% limit of agreement and biases of phantom test measurement under RT-setting between two observers were very small. Excellent inter-observer agreement (ICC >0.75) was achieved in all measurements except for slice thickness accuracy (ICC =0.42, moderate agreement) under RT-setting.

CONCLUSIONS

Very good and highly reproducible image quality could be achieved on a 1.5T MR-sim with a flexible coil setting as revealed by the standard ACR MRI phantom test. The flexible RT-setting compromised in image signal-to-noise ratio (SNR) compared to the normal DR-setting, and resulted in reduced low-contrast object detectability.

Filmless methods for quality assurance of Tomotherapy using ArcCHECK

PURPOSE

Tomotherapy delivers an intensity-modulated radiation therapy (IMRT) treatment by the synchronization of gantry rotation, multileaf collimator (MLC), and couch movement. This dynamic nature makes the quality assurance (QA) important and challenging. The purpose of this study is to develop some methodologies using an ArcCHECK for accurate QA measurements of the gantry angle and speed, MLC synchronization and leaf open time, couch translation per gantry rotation, couch speed and uniformity, and constancy of longitudinal beam profile for a Tomotherapy unit.

METHODS

Four test plans recommended by AAPM Task Group 148 (TG148) and the manufacturer were chosen for this study. Helical and static star shot tests are used for checking the leaves opened at the expected gantry angles. Another helical test is to verify the couch traveled the expected distance per gantry rotation. The final test is for checking the couch speed constancy with a static gantry. ArcCHECK can record the detector signal every 50 ms as a movie file, and has a virtual inclinometer for gantry angle measurement. These features made the measurement of gantry angle and speed, MLC synchronization and leaf open time, and longitudinal beam profile possible. A shaping parameter was defined for facilitating the location of the beam center during the plan delivery, which was thereafter used to calculate the couch translation per gantry rotation and couch speed. The full width at half maximum (FWHM) was calculated for each measured longitudinal beam profile and then used to evaluate the couch speed uniformity. Furthermore, a mean longitudinal profile was obtained for constancy check of field width. The machine trajectory log data were also collected for comparison. Inhouse programs were developed in MATLAB to process both the ArcCHECK and machine log data.

RESULTS

The deviation of our measurement results from the log data for gantry angle was calculated to be less than 0.4°. The percentage differences between measured and planned leaf open time were found to be within 0.5% in all the tests. Our results showed mean values of MLC synchronization of 0.982, 0.983, and 0.995 at static gantry angle 0°, 45°, and 135°, respectively. The mean value of measured couch translation and couch speed by ArcCHECK had less than 0.1% deviation from the planned values. The variation in the value of FWHM suggested the couch speed uniformity was better than 1%. The mean of measured longitudinal profiles was suitable for constancy check of field width.

CONCLUSION

Precise and efficient methods for measuring the gantry angle and speed, leaf open time, couch translation per gantry rotation, couch speed and uniformity, and constancy of longitudinal beam profile of Tomotherapy using ArcCHECK have been developed and proven to be accurate compared with machine log data. Estimation of the Tomotherapy binary MLC leaf open time is proven to be precise enough to verify the leaf open time as small as 277.8 ms. Our method also makes the observation and quantification of the synchronization of leaves possible.

Technical note: Correlation between TQA data trends and TomoHD functional status

TomoTherapy Quality Assurance (TQA) is a software package developed to monitor certain aspects of machine performance. In this study, the TQA quantities or data trends most effective in monitoring energy drifts and magnetron stability were determined respectively. This retrospective study used data collected from three TomoHD units. The TQA modules investigated were Step‐Wedge Helical, Step‐Wedge Static, and Basic Dosimetry. First, the TQA quantities correlated with energy changes (|r|>0.85, where r is the Pearson's correlation coefficient) were found. The corresponding sensitivities to percentage depth dose (PDD) ratio changes were then calculated and compared. Second, the pulse‐by‐pulse dose stability was compared before and after each magnetron replacement using a nonparametric comparison test (Welch's t‐test), and the raw dose profiles were surveyed. In this study, exit detector flatness obtained in Basic Dosimetry was shown to be the most sensitive (r=0.945) to energy changes, followed by the energy differences in Step‐Wedge Static (r=0.942) and Step‐Wedge Helical (r=0.898). The three quantities could detect a PDD ratio change of 5.1×10−4,5.4×10−4, and 7.1×10−4, respectively. Pulse‐by‐Pulse Dose1 from Basic Dosimetry over a one‐week period before and after a magnetron replacement showed a significant difference (p<0.05) in only three of the nine instances. On the other hand, a raw output profile free from discontinuities, frequent dropped pulses and abnormal spikes was found to indicate that the magnetron would continue to function normally for a week 89% of the time.

PACS numbers: 87.56.bd, 87.56.Fc, 84.40.Fe

Accuracy in contouring of small and low contrast lesions: comparison between diagnostic quality computed tomography scanner and computed tomography simulation scanner-A phantom study

To evaluate the accuracy in detection of small and low-contrast regions using a high-definition diagnostic computed tomography (CT) scanner compared with a radiotherapy CT simulation scanner. A custom-made phantom with cylindrical holes of diameters ranging from 2-9 mm was filled with 9 different concentrations of contrast solution. The phantom was scanned using a 16-slice multidetector CT simulation scanner (LightSpeed RT16, General Electric Healthcare, Milwaukee, WI) and a 64-slice high-definition diagnostic CT scanner (Discovery CT750 HD, General Electric Healthcare). The low-contrast regions of interest (ROIs) were delineated automatically upon their full width at half maximum of the CT number profile in Hounsfield units on a treatment planning workstation. Two conformal indexes, CI(in), and CI(out), were calculated to represent the percentage errors of underestimation and overestimation in the automated contours compared with their actual sizes. Summarizing the conformal indexes of different sizes and contrast concentration, the means of CI(in) and CI(out) for the CT simulation scanner were 33.7% and 60.9%, respectively, and 10.5% and 41.5% were found for the diagnostic CT scanner. The mean differences between the 2 scanners' CI(in) and CI(out) were shown to be significant with p < 0.001. A descending trend of the index values was observed as the ROI size increases for both scanners, which indicates an improved accuracy when the ROI size increases, whereas no observable trend was found in the contouring accuracy with respect to the contrast levels in this study. Images acquired by the diagnostic CT scanner allow higher accuracy on size estimation compared with the CT simulation scanner in this study. We recommend using a diagnostic CT scanner to scan patients with small lesions (<1 cm in diameter) for radiotherapy treatment planning, especially for those pending for stereotactic radiosurgery in which accurate delineation of small-sized, low-contrast regions is important for dose calculation.

Characterization of Baylisascaris schroederi from Qinling subspecies of giant panda in China by the first internal transcribed spacer (ITS-1) of nuclear ribosomal DNA

In the present study, a total of 20 nematode isolates, (including 10 male and 10 female worms) representing Baylisascaris schroederi from 5 Qinling subspecies of giant pandas (Ailuropoda melanoleuca) in Shaanxi Province of China, were characterized and grouped genetically by the first internal transcribed spacer (ITS-1) of nuclear ribosomal DNA (rDNA). The rDNA fragment spanning 3' end of 18S rDNA, complete ITS-1 rDNA, and 5' end of 5.8S rDNA were amplified and sequenced. The sequence variability in ITS-1 rDNA was examined within B. schroederi and among parasites in order Ascaridata available in GenBank™, and their phylogenetic relationships were also reconstructed. The sequences of ITS-1 rDNA for all the B. schroederi isolates were 427 bp in length, with no genetic variation detected among these isolates. Phylogenetic analyses based on the ITS-1 rDNA sequences revealed that all the male and female B. schroederi isolates sequenced in the present study were posited into the clade of genus Baylisascaris, sistered to zoonotic nematodes in genus Ascaris, and the ITS-1 rDNA sequence could distinguish different species in order Ascaridata. These results showed that the ITS-1 rDNA provides a suitable molecular marker for the inter-species phylogenetic analysis and differential identification of nematodes in order Ascaridata.

Rarefied-gas heat transfer in micro- and nanoscale Couette flows

The physics of the heat conduction and viscous dissipation in rarefied gases is analyzed and discussed. A heat transfer model valid for arbitrary Knudsen numbers, defined as the ratio of the molecular mean free path to the characteristic length of channels, is derived by treating the heat transfer behavior in the slip and transition regimes as an intermediate function of continuum heat transfer model and free molecular heat transfer model. Comparison studies reveal that this model not only shows good agreement with the numerical results based on the direct simulation Monte Carlo method, but also has some unique features that can overcome the deficiencies existing in the previous models. Therefore, this model is capable to study the heat transport phenomena in very dilute gas Couette flows through micro/nanochannels more accurately.

Resistance to adefovir dipivoxil in lamivudine resistant chronic hepatitis B patients treated with adefovir dipivoxil

BACKGROUND

Adefovir dipivoxil (ADV) is a potent nucleotide analogue against both the wild-type and lamivudine (LMV) resistant hepatitis B virus (HBV). The cumulative incidence of ADV resistant mutations in the nucleoside/-tide treatment naive chronic hepatitis B patient (CHB) at weeks 48, 96, and 144 was 0, 0.8-3%, and approximately 5.9%, respectively.

AIMS

The aim of this study was to characterise the genotypic and phenotypic mutation profiles to ADV in 67 LMV resistant CHB patients who were treated with ADV.

METHODS

Serum HBV DNA was quantified by real time polymerase chain reaction. The ADV mutant was detected using matrix assisted laser desorption/ionisation time of flight mass spectrometry based genotyping assays, termed restriction fragment mass polymorphism (RFMP).

RESULTS

RFMP analysis revealed that a total of 11 amino acid substitutions developed in the rt domain of the HBV polymerase in nine patients. The cumulative incidence of genotypic ADV resistance at months 12 and 24 was 6.4% and 25.4%, respectively. The rtA181V, rtN236T, and rtA181T mutations were detected in five, four, and two of the 67 patients at treatment months 12-17, 3-19, and 7-20, respectively. Serial quantification of serum HBV DNA revealed that two patients with the rtA181V mutation, with or without the rtN236T mutation, and one patient with the rtA181T mutation displayed HBV DNA rebound.

CONCLUSION

Emergence of the ADV mutation in LMV resistant patients who are treated with ADV appeared to present earlier and more frequently than was reported in previous studies on nucleoside/-tide treatment naive patients.

MAPD--an objective way to select mAs for paediatric brain CT

CT is an advanced imaging modality, but the imaging parameters are normally selected subjectively. For standard head examinations, most of the parameters used are consistent amongst different centres, with the exception of large variations in the selection of the tube current-exposure time product (mAs). As a result, CT images may contain unacceptable levels of noise, or the patient may receive excessive radiation. In this study, the maximum anteroposterior diameter (MAPD) was shown to be a good criterion for mAs selection, and could be measured in a pilot view. 200 paediatric brain CT studies were randomly selected to determine the MAPD at the mid brain level. With knowledge of MAPD distribution, a phantom study was performed to determine the relationship between MAPD and the mAs required for consistent and acceptable image noise. It was found that the required mAs increased linearly with MAPD. Assuming the manufacturer's recommended value is "appropriate" for the average MAPD, the appropriate mAs value could be estimated. Using this method, appropriate mAs values were calculated retrospectively for a group of 240 randomly selected paediatric brain CT studies and compared with the actual mAs subjectively determined by the radiographer. Although their average values were similar, the difference between the calculated and actual values deviated markedly in some cases. When the actual mAs was smaller than the calculated value, higher image noise was observed. However, reduction of image noise was barely observed when the applied mAs was larger than the calculated value. Thus, this method is more objective and appropriate for determination of the mAs value for paediatric brain CT than the traditional subjective method.