Technical note: GAMMORA,a free,open-source,and validated GATE-based model for Monte-Carlo simulations of the Varian TrueBeam

PurposeMonte Carlo (MC) is the reference computation method for medical physics. In radiotherapy, MC computations are necessary for some issues (such as assessing figures of merit, double checks, and dose conversions). A tool based on GATE is proposed to easily create full MC simulations of the Varian TrueBeam STx.MethodsGAMMORA is a package that contains photon phase spaces as a pre-trained generative adversarial network (GAN) and the TrueBeam’s full geometry. It allows users to easily create MC simulations for simple or complex radiotherapy plans such as VMAT. To validate the model, the characteristics of generated photons are first compared to those provided by Varian (IAEA format). Simulated data are also compared to measurements in water and heterogeneous media. Simulations of 8 SBRT plans are compared to measurements (in a phantom). Two examples of applications (a second check and interplay effect assessment) are presented.ResultsThe simulated photons generated by the GAN have the same characteristics (energy, position, and direction) as the IAEA data. Computed dose distributions of simple cases (in water) and complex plans delivered in a phantom are compared to measurements, and the Gamma index (3%/3mm) was always superior to 98%. The feasibility of both clinical applications is shown.ConclusionsThis model is now shared as a free and open-source tool that generates radiotherapy MC simulations. It has been validated and used for five years. Several applications can be envisaged for research and clinical purposes.     

Test of Accuracy of LAI Estimation by LAI-2000 under Artificially Changed Leaf to Wood Area Proportions

The accuracy of LAI-2000 Plant Canopy Analyzer for leaf (LAI) and plant (PAI) area indexes measurements was tested in 20-year-old Norway spruce stand using the reduction of canopy biomass. Needle and branch areas were reduced progressively upward every one meter. Values of effective leaf area index (LAI_e), as an uncorrected product of LAI-2000, were compared with directly estimated LAI and PAI values after each reduction step. LAI-2000 underestimates PAI and LAI values according to LAI-2000 rings readings, and varied proportions between leaf and wood areas. The values of LAI_c have been increased with decreasing of the view angle of the relevant LAI-2000 rings. Therefore, the underestimation of LAI becomes smaller when the readings near the horizon are masked. More accurate results, for projected LAI (LAI_p) calculation, are produced by LAI-2000 when some dense grids of measurement points and the most vertical ring readings (0 –13 °) are used. Correction factor 1.6 is possible to use for unreduced canopy hemi-surface LAI estimation, when the last rings (i.e. 5^th and 4^th rings, 47 –74 °) are excluded. Correction factor of 1.25 can be used to compute LAI_p if the angle readings under 43 °are also masked.     

Are the numbers adding up? Exploiting discrepancies among complementary population models

Large carnivores are difficult to monitor because they tend to be sparsely distributed, sensitive to human activity, and associated with complex life histories. Consequently, understanding population trend and viability requires conservationists to cope with uncertainty and bias in population data. Joint analysis of combined data sets using multiple models (i.e., integrated population model) can improve inference about mechanisms (e.g., habitat heterogeneity and food distribution) affecting population dynamics. However, unobserved or unobservable processes can also introduce bias and can be difficult to quantify. We developed a Bayesian hierarchical modeling approach for inference on an integrated population model that reconciles annual population counts with recruitment and survival data (i.e., demographic processes). Our modeling framework is flexible and enables a realistic form of population dynamics by fitting separate density-dependent responses for each demographic process. Discrepancies estimated from shared parameters among different model components represent unobserved additions (i.e., recruitment or immigration) or removals (i.e., death or emigration) when annual population counts are reliable. In a case study of gray wolves in Wisconsin (1980–2011), concordant with policy changes, we estimated that a discrepancy of 0% (1980–1995), −2% (1996–2002), and 4% (2003–2011) in the annual mortality rate was needed to explain annual growth rate. Additional mortality in 2003–2011 may reflect density-dependent mechanisms, changes in illegal killing with shifts in wolf management, and nonindependent censoring in survival data. Integrated population models provide insights into unobserved or unobservable processes by quantifying discrepancies among data sets. Our modeling approach is generalizable to many population analysis needs and allows for identifying dynamic differences due to external drivers, such as management or policy changes.     

Finite element simulation of spinal deformities correction by in situ contouring technique

Biomechanical models have been proposed in order to simulate the surgical correction of spinal deformities. With these models, different surgical correction techniques have been examined: distraction and rod rotation. The purpose of this study was to simulate another surgical correction technique: the in situ contouring technique. In this way, a comprehensive three-dimensional Finite Element (FE) model with patient-specific geometry and patient-specific mechanical properties was used. The simulation of the surgery took into account elasto–plastic behavior of the rod and multiple moments loading and unloading representing the surgical maneuvers. The simulations of two clinical cases of hyperkyphosis and scoliosis were coherent with the surgeon's experience. Moreover, the results of simulation were compared to post-operative 3D measurements. The mean differences were under 5° for vertebral rotations and 5 mm for spinal lines. These simulations open the way for future predictive tools for surgical planning.     

Quantification of blockiness in pectins—A comparative study using vibrational spectroscopy and chemometrics

The gelling properties of pectins are related not only to the degree of esterification (DE), but also to the distribution of the ester groups. In this study, we have examined an experimentally designed series of 31 pectins originating from the same mother pectin and de-esterified using combinations of two different enzymatic mechanisms. The potential of using infrared (IR), Raman, and near infrared (NIR) spectroscopies combined with chemometrics for reliable and rapid determination of the DE and distribution patterns of methyl ester groups in a designed set of pectin powders was investigated. Quantitative calibration models using partial least squares (PLS) regression were developed and compared. The calibration models for prediction of DE obtained on extended inverse signal correction (EISC)-treated spectra of all three spectroscopic methods yielded models with cross-validated prediction errors (RMSECV) between 1.1%p and 1.6%p DE and correlation coefficients of 0.99. A calibration model predicting degree of random de-esterification (R) and block de-esterification (B) was developed for each spectroscopic method, yielding RMSECV values between 4.4 and 6.7 and correlation coefficients (r) between 0.79 and 0.92. Variable selection using interval PLS (iPLS) significantly improved the prediction of R for IR spectroscopy, yielding RMSECV of 3.5 and correlation coefficients of 0.95. All three spectroscopic methods were able to distinguish the spectral patterns of pectins with different enzyme treatments in simple classification models by principal component analysis (PCA). Extended canonical variate analysis revealed one specific signal in the Raman (1045 cm⁻¹) spectrum and one significant area (1250-1400 cm⁻¹) in the IR spectrum which are able to classify the pectin samples according to the four different enzyme treatments. In both Raman and IR spectra, the signal intensity decreased in the sequence R-B > B > B-R > R > re-methylated pectin.     

Enhancement of soft-tissue contrast in cone-beam CT using an anti-scatter grid with a sparse sampling approach

PurposeAnti-scatter grids suppress the scatter substantially thus improving image contrast in radiography. However, its active use in cone-beam CT for the purpose of improving contrast-to-noise ratio (CNR) has not been successful mainly due to the increased noise related to Poisson statistics of photons. This paper proposes a sparse-view scanning approach to address the above issue.MethodCompared to the conventional cone-beam CT imaging framework, the proposed method reduces the number of projections and increases exposure in each projection to enhance image quality without an additional cost of radiation dose to patients. For image reconstruction from sparse-view data, an adaptive-steepest-descent projection-onto-convex-sets (ASD POCS) algorithm regularized by total-variation (TV) minimization was adopted. Contrast and CNR with various scattering conditions were evaluated in projection domain by a simulation study using GATE. Then we evaluated contrast, resolution, and image uniformity in CT image domain with Catphan phantom. A head phantom with soft-tissue structures was also employed for demonstrating a realistic application. A virtual grid-based estimation and reduction of scatter has also been implemented for comparison with the real anti-scatter grid.ResultsIn the projection domain evaluation, contrast and CNR enhancement was observed when using an anti-scatter grid compared to the virtual grid. In the CT image domain, the proposed method produced substantially higher contrast and CNR of the low-contrast structures with much improved image uniformity.ConclusionWe have shown that the proposed method can provide high-quality CBCT images particularly with an increased contrast of soft-tissue at a neutral dose for image-guidance.     

Small Sample Properties of the Mantel-Haenszel Test Statistic for Density Follow-Up Studies

For the analysis of combinations of 2×2 non-contingency tables as obtained from density follow-up studies (relating a number of events to a number of person-years of follow-up) an analogue of the Mantel-Haenszel test for 2×2 contingency tables is widely used. In this paper the small sample properties of this test, both with and without continuity correction, are evaluated. Also the improvement of the test-statistic by using the first four cumulants via the Edgeworth expansion was studied. Results on continuity correction agree with similar studies on the Mantel-Haenszel statistic for 2×2 contingency tables: Continuity correction gives a p-value which approximates the exact p-value better than the p-value obtained without this correction; both the exact test and its approximations show considerable conservatism in small samples; the uncorrected Mantel-Haenszel test statistic gives a p-value that agrees more with the nominal significance level, but can be anti-conservative. The p-value based on the first four cumulants gives a better approximation of the exact p-value than the continuity corrected test, especially when the distribution has marked skewness.     

Evaluation of three scatter correction methods based on estimation of photopeak scatter spectrum in SPECT imaging: A simulation study

Three practical methods for scatter correction of Tc-99m SPECT images are evaluated. Among these, two methods, three-energy window (TEW) methods using the trapezoidal and triangular approximations, have been described previously by investigators, and a new approximation is offered in this work. The SIMIND (SIMulation of Imaging Nuclear Detectors) Monte Carlo program is used to simulate a line source placed at on-axis and 5 cm off-axis locations, a cold-sphere/hot-background phantom, a hot-sphere/cold-background phantom, and a more clinically realistic NCAT (Nonuniform Rational B-spline-based CArdiac-Torso) phantom. For evaluation of these methods, the scatter line-spread functions and scatter fractions for the on- and off-axis line source, image contrast, signal-to-noise ratio and relative noise for the cold spheres, and recovery coefficient for the hot spheres of different diameters are compared. For the NCAT phantom, a line profile through a slice of the reconstructed image is considered before and after scatter correction, and also image contrast defined by this profile is used to compare the correction methods. The results of this study indicate that for the line source simulation the scatter fractions obtained from the proposed method are a better estimation of true scatter fractions. Also, for both the sphere simulation and NCAT simulation, the proposed method improves the image contrast as compared to the two other methods.     

On corrected score approach for proportional hazards model with covariate measurement error

In the presence of covariate measurement error with the proportional hazards model, several functional modeling methods have been proposed. These include the conditional score estimator (Tsiatis and Davidian, 2001, Biometrika 88, 447-458), the parametric correction estimator (Nakamura, 1992, Biometrics 48, 829-838), and the nonparametric correction estimator (Huang and Wang, 2000, Journal of the American Statistical Association 95, 1209-1219) in the order of weaker assumptions on the error. Although they are all consistent, each suffers from potential difficulties with small samples and substantial measurement error. In this article, upon noting that the conditional score and parametric correction estimators are asymptotically equivalent in the case of normal error, we investigate their relative finite sample performance and discover that the former is superior. This finding motivates a general refinement approach to parametric and nonparametric correction methods. The refined correction estimators are asymptotically equivalent to their standard counterparts, but have improved numerical properties and perform better when the standard estimates do not exist or are outliers. Simulation results and application to an HIV clinical trial are presented.