Alexandr Malusek

Use of energy deposition spectrometer Liulin for individual monitoring of aircrew

Silicon energy deposition spectrometer Liulin was primarily developed for cosmic radiation monitoring onboard spacecrafts. Nowadays, Liulin type detectors are also used to characterise radiation field on board aircraft, at alpine observatories and behind the shielding of heavy ion accelerators. In this work, experiments and calibrations performed in these radiation fields are presented and the method developed for calculation of ambient dose equivalent H*(10) on board aircraft is described. Since 2001, a simple method employing the energy deposition spectra had been used to determine H*(10) on board aircraft but, in 2004, it became clear that the resulting values were strongly biased at locations close to Earths equator. An improved method for the determination of H*(10) on board aircraft using the Liulin detector was developed. It took into account the composition of the radiation field via the ratio of absorbed doses D(low) and D(neut) reflecting the contributions from low-LET particles and neutrons, respectively. It resulted in much better agreement with the EPCARD computer code for all aircraft locations; relative differences were within 11 % for low-LET and 20 % for neutron components of H*(10).

CTmod-A toolkit for Monte Carlo simulation of projections including scatter in computed tomography

The CTmod toolkit is a set of C++ class libraries based on the CERNs application development framework ROOT. It uses the Monte Carlo method to simulate energy imparted to a CT-scanner detector array. Photons with a given angle-energy distribution are emitted from the X-ray tube approximated by a point source, transported through a phantom, and their contribution to the energy imparted per unit surface area of each detector element is scored. Alternatively, the scored quantity may be the fluence, energy fluence, plane fluence, plane energy fluence, or kerma to air in the center of each detector element. Phantoms are constructed from homogenous solids or voxel arrays via overlapping. Implemented photon interactions (photoelectric effect, coherent scattering, and incoherent scattering) are restricted to the energy range from 10 to 200keV. Variance reduction techniques include the collision density estimator and survival biasing combined with the Russian roulette. The toolkit has been used to estimate the amount of scatter in cone beam computed tomography and planar radiography.

Iterative reconstruction for quantitative tissue decomposition in dual-energy CT

Quantitative tissue classification using dual-energy CT has the potential to improve accuracy in radiation therapy dose planning as it provides more information about material composition of scanned objects than the currently used methods based on single-energy CT. One problem that hinders successful application of both single- and dualenergy CT is the presence of beam hardening and scatter artifacts in reconstructed data. Current pre- and post-correction methods used for image reconstruction often bias CT numbers and thus limit their applicability for quantitative tissue classification. Here we demonstrate simulation studies with a novel iterative algorithm that decomposes every soft tissue voxel into three base materials: water, protein and adipose. The results demonstrate that beam hardening artifacts can effectively be removed and accurate estimation of mass fractions of all base materials can be achieved. In the future, the algorithm may be developed further to include segmentation of soft and bone tissue and subsequent bone decomposition, extension from 2-D to 3-D and scatter correction.

Monte Carlo study of the dependence of the KAP-meter calibration coefficient on beam aperture, x-ray tube voltage and reference plane

The Monte Carlo method was used to study the dependence of the calibration coefficient on the tube voltage, beam aperture and reference plane in simplified over-couch geometries modelling VacuTecs type 70157 KAP-meter both with and without an additional filter. The MCNP5 code was used to calculate (i) energy imparted to air cavities of the KAP-meter and (ii) spatial distribution of air collision kerma at entrance and exit planes of the KAP-meter and at a plane close to the patient. From these data, the air kerma area product and calibration coefficient were calculated and their dependence on the tube voltage and beam aperture was analysed. It was found that the variation of the calibration coefficient as a function of tube voltage was up to 40% when the additional filter was used. The additional filter placed closely in front of the KAP-meter decreased the calibration coefficient for the patient plane by about 10% compared to the ideal additional filter. The effect of the beam aperture was small at the patient plane and negligible for the exit plane.

Liulin silicon semiconductor spectrometers as cosmic ray monitors at the high mountain observatories Jungfraujoch and Lomnický štít

Currently, most cosmic ray data are obtained by detectors on satellites, aircraft, high-altitude balloons and ground (neutron monitors). In our work, we examined whether Liulin semiconductor spectrometers (simple silicon planar diode detectors with spectrometric properties) located at high mountain observatories could contribute new information to the monitoring of cosmic rays by analyzing data from selected solar events between 2005 and 2013. The decision thresholds and detection limits of these detectors placed at Jungfraujoch (Switzerland; 3475 m a.s.l.; vertical cut-off rigidity 4.5 GV) and Lomnický stit (Slovakia; 2633 m a.s.l.; vertical cut-off rigidity 3.84 GV) high-mountain observatories were determined. The data showed that only the strongest variations of the cosmic ray flux in this period were detectable. The main limitation in the performance of these detectors is their small sensitive volume and low sensitivity of the PIN photodiode to neutrons.

AUTOMATIC SEGMENTATION OF PELVIS FOR BRACHYTHERAPY OF PROSTATE

Advanced model-based iterative reconstruction algorithms in quantitative computed tomography (CT) perform automatic segmentation of tissues to estimate material properties of the imaged object. Compared with conventional methods, these algorithms may improve quality of reconstructed images and accuracy of radiation treatment planning. Automatic segmentation of tissues is, however, a difficult task. The aim of this work was to develop and evaluate an algorithm that automatically segments tissues in CT images of the male pelvis. The newly developed algorithm (MK2014) combines histogram matching, thresholding, region growing, deformable model and atlas-based registration techniques for the segmentation of bones, adipose tissue, prostate and muscles in CT images. Visual inspection of segmented images showed that the algorithm performed well for the five analysed images. The tissues were identified and outlined with accuracy sufficient for the dual-energy iterative reconstruction algorithm whose aim is to improve the accuracy of radiation treatment planning in brachytherapy of the prostate.

PARALLELISATION OF THE MODEL-BASED ITERATIVE RECONSTRUCTION ALGORITHM DIRA

New paradigms for parallel programming have been devised to simplify software development on multi-core processors and many-core graphical processing units (GPU). Despite their obvious benefits, the parallelisation of existing computer programs is not an easy task. In this work, the use of the Open Multiprocessing (OpenMP) and Open Computing Language (OpenCL) frameworks is considered for the parallelisation of the model-based iterative reconstruction algorithm DIRA with the aim to significantly shorten the codes execution time. Selected routines were parallelised using OpenMP and OpenCL libraries; some routines were converted from MATLAB to C and optimised. Parallelisation of the code with the OpenMP was easy and resulted in an overall speedup of 15 on a 16-core computer. Parallelisation with OpenCL was more difficult owing to differences between the central processing unit and GPU architectures. The resulting speedup was substantially lower than the theoretical peak performance of the GPU; the cause was explained.

In-situ calibration of clinical built-in KAP meters with traceability to a primary standard using a reference KAP meter

The air kerma-area product (KAP) is used for settings of diagnostic reference levels. The International Atomic Energy Agency (IAEA) recommends that doses in diagnostic radiology (including the KAP values) be estimated with an accuracy of at least ± 7% (k = 2). Industry standards defined by the International Electrotechnical Commission (IEC) specify that the uncertainty of KAP meter measurements should be less than ± 25% (k = 2). Medical physicists willing to comply with the IAEAs recommendation need to apply correction factors to KAP values reported by x-ray units. The aim of this work is to present and evaluate a calibration method for built-in KAP meters on clinical x-ray units. The method is based on (i) a tandem calibration method, which uses a reference KAP meter calibrated to measure the incident radiation, (ii) measurements using an energy-independent ionization chamber to correct for the energy dependence of the reference KAP meter, and (iii) Monte Carlo simulations of the beam quality correction factors that correct for differences between beam qualities at a standard laboratory and the clinic. The method was applied to the KAP meter in a Siemens Aristos FX plus unit. It was found that values reported by the built-in KAP meter differed from the more accurate values measured by the reference KAP meter by more than 25% for high tube voltages (more than 140 kV) and heavily filtered beams (0.3 mm Cu). Associated uncertainties were too high to claim that the IECs limit of 25% was exceeded. Nevertheless the differences were high enough to justify the need for a more accurate calibration of built-in KAP meters.

ACCURATE KAP METER CALIBRATION AS A PREREQUISITE FOR OPTIMISATION IN PROJECTION RADIOGRAPHY

Modern X-ray units register the air kerma-area product, PKA, with a built-in KAP meter. Some KAP meters show an energy-dependent bias comparable with the maximum uncertainty articulated by the IEC (25 %), adversely affecting dose-optimisation processes. To correct for the bias, a reference KAP meter calibrated at a standards laboratory and two calibration methods described here can be used to achieve an uncertainty of <7 % as recommended by IAEA. A computational model of the reference KAP meter is used to calculate beam quality correction factors for transfer of the calibration coefficient at the standards laboratory, Q0, to any beam quality, Q, in the clinic. Alternatively, beam quality corrections are measured with an energy-independent dosemeter via a reference beam quality in the clinic, Q1, to beam quality, Q Biases up to 35 % of built-in KAP meter readings were noted. Energy-dependent calibration factors are needed for unbiased PKA Accurate KAP meter calibration as a prerequisite for optimisation in projection radiography.

Prostate tissue decomposition via DECT using the model based iterative image reconstruction algorithm DIRA

Better knowledge of elemental composition of patient tissues may improve the accuracy of absorbed dose delivery in brachytherapy. Deficiencies of water-based protocols have been recognized and work is ongoing to implement patient-specific radiation treatment protocols. A model based iterative image reconstruction algorithm DIRA has been developed by the authors to automatically decompose patient tissues to two or three base components via dual-energy computed tomography. Performance of an updated version of DIRA was evaluated for the determination of prostate calcification. A computer simulation using an anthropomorphic phantom showed that the mass fraction of calcium in the prostate tissue was determined with accuracy better than 9%. The calculated mass fraction was little affected by the choice of the material triplet for the surrounding soft tissue. Relative differences between true and approximated values of linear attenuation coefficient and mass energy absorption coefficient for the prostate tissue were less than 6% for photon energies from 1 keV to 2 MeV. The results indicate that DIRA has the potential to improve the accuracy of dose delivery in brachytherapy despite the fact that base material triplets only approximate surrounding soft tissues.