Jani Keyriläinen

A dual model HU conversion from MRI intensity values within and outside of bone segment for MRI‐based radiotherapy treatment planning of prostate cancer

PURPOSE The lack of electron density information in magnetic resonance images (MRI) poses a major challenge for MRI-based radiotherapy treatment planning (RTP). In this study the authors convert MRI intensity values into Hounsfield units (HUs) in the male pelvis and thus enable accurate MRI-based RTP for prostate cancer patients with varying tissue anatomy and body fat contents. METHODS T1/T2*-weighted MRI intensity values and standard computed tomography (CT) image HUs in the male pelvis were analyzed using image data of 10 prostate cancer patients. The collected data were utilized to generate a dual model HU conversion technique from MRI intensity values of the single image set separately within and outside of contoured pelvic bones. Within the bone segment local MRI intensity values were converted to HUs by applying a second-order polynomial model. This model was tuned for each patient by two patient-specific adjustments: MR signal normalization to correct shifts in absolute intensity level and application of a cutoff value to accurately represent low density bony tissue HUs. For soft tissues, such as fat and muscle, located outside of the bone contours, a threshold-based segmentation method without requirements for any patient-specific adjustments was introduced to convert MRI intensity values into HUs. The dual model HU conversion technique was implemented by constructing pseudo-CT images for 10 other prostate cancer patients. The feasibility of these images for RTP was evaluated by comparing HUs in the generated pseudo-CT images with those in standard CT images, and by determining deviations in MRI-based dose distributions compared to those in CT images with 7-field intensity modulated radiation therapy (IMRT) with the anisotropic analytical algorithm and 360° volumetric-modulated arc therapy (VMAT) with the Voxel Monte Carlo algorithm. RESULTS The average HU differences between the constructed pseudo-CT images and standard CT images of each test patient ranged from -2 to 5 HUs and from 22 to 78 HUs in soft and bony tissues, respectively. The average local absolute value differences were 11 HUs in soft tissues and 99 HUs in bones. The planning target volume doses (volumes 95%, 50%, 5%) in the pseudo-CT images were within 0.8% compared to those in CT images in all of the 20 treatment plans. The average deviation was 0.3%. With all the test patients over 94% (IMRT) and 92% (VMAT) of dose points within body (lower than 10% of maximum dose suppressed) passed the 1 mm and 1% 2D gamma index criterion. The statistical tests (t- and F-tests) showed significantly improved (p ≤ 0.05) HU and dose calculation accuracies with the soft tissue conversion method instead of homogeneous representation of these tissues in MRI-based RTP images. CONCLUSIONS This study indicates that it is possible to construct high quality pseudo-CT images by converting the intensity values of a single MRI series into HUs in the male pelvis, and to use these images for accurate MRI-based prostate RTP dose calculations.

Human breast cancer in vitro: matching histo-pathology with small-angle x-ray scattering and diffraction enhanced x-ray imaging

Twenty-eight human breast tumour specimens were studied with small-angle x-ray scattering (SAXS), and 10 of those were imaged by the diffraction enhanced x-ray imaging (DEI) technique. The sample diameter was 20 mm and the thickness 1 mm. Two examples of ductal carcinoma are illustrated by histology images, DEI, and maps of the collagen d-spacing and scattered intensity in the Porod regime, which characterize the SAXS patterns from collagen-rich regions of the samples. Histo-pathology reveals the cancer-invaded regions, and the maps of the SAXS parameters show that in these regions the scattering signal differs significantly from scattering by the surrounding tissue, indicating a degradation of the collagen structure in the invaded regions. The DEI images show the borders between collagen and adipose tissue and provide a co-ordinate system for tissue mapping by SAXS. In addition, degradation of the collagen structure in an invaded region is revealed by fading contrast of the DEI refraction image. The 28 samples include fresh, defrosted tissue and formalin-fixed tissue. The d-values with their standard deviations are given. In the fresh samples there is a systematic 0.76% increase of the d-value in the invaded regions, averaged over 11 samples. Only intra-sample comparisons are made for the formalin-fixed samples, and with a long fixation time, the difference in the d-value stabilizes at about 0.7%. The correspondence between the DEI images, the SAXS maps and the histo-pathology suggests that definitive information on tumour growth and malignancy is obtained by combining these x-ray methods.

Phase-contrast X-ray imaging of breast

When an X-ray wave traverses an object, its amplitude and phase change, resulting in attenuation, interference, and refraction, and in phase-contrast X-ray imaging (PCI) these are converted to intensity changes. The relative change of the X-ray phase per unit path length is even orders of magnitude larger than that of the X-ray amplitude, so that the image contrast based on variation of the X-ray phase is potentially much stronger than the contrast based on X-ray amplitude (absorption contrast). An important medical application of PCI methods is soft-tissue imaging, where the absorption contrast is inherently weak. It is shown by in vitro examples that signs of malignant human breast tumor are enhanced in PCI images. Owing to the strong contrast, the radiation dose can be greatly reduced, so that a high-resolution phase-contrast X-ray tomography of the breast is possible with about 1 mGy mean glandular dose. Scattered radiation carries essential information on the atomic and molecular structure of the object, and particularly small-angle X-ray scattering can be used to trace cancer. The imaging methods developed at the synchrotron radiation facilities will become available in the clinical environment with the ongoing development of compact radiation sources, which produce intense X-ray beams of sufficient coherence. Several developments that are under way are described here.

Pulmonary toxicity after radiotherapy in primary breast cancer patients: Results from a randomized chemotherapy study

PURPOSE Pulmonary toxicity was prospectively evaluated within a randomized trial for breast cancer patients at high risk for relapse, who postoperatively received as adjuvant therapy either 9 cycles of tailored chemotherapy (20 patients) (cyclophosphamide, epirubicin, 5-fluorouracil [FEC]) or standard FEC x 3 followed by high-dose chemotherapy (cyclophosphamide, thiotepa, carboplatin [CTCb]) supported by peripheral blood stem cell transplantation (14 patients). After high-dose chemotherapy or tailored FEC, all patients received locoregional radiotherapy (50 Gy/5 weeks), plus tamoxifen for 5 years. METHODS AND MATERIALS Lung function tests (FVC, FEV1, and DL(CO)) were performed before chemotherapy and 9 months after radiotherapy. Computed tomography of the lungs was performed before radiotherapy and 6 weeks, 3 months, and 9 months after radiotherapy. RESULTS Clinical signs of suspected pneumonitis were noted in 29% of patients, but only 1 patient needed symptomatic therapy. Radiologic changes were detected in 68% of patients, and they were most frequent at 3 months after radiotherapy. FVC decreased in both groups (tailored FEC: mean difference, -6.5%, p = 0.0005; CTCb: -2.0%, p = 0.21; tailored FEC vs. CTCb: -4.5%, p = 0.05). DL(CO) decreased significantly in both groups (tailored FEC: mean difference, -11.2%, p < 0.0001; CTCb: -5.6%, p = 0.02; tailored FEC vs. CTCb: -5.6%, p = 0.07). FEV1 decreased by 7.3% in patients treated with tailored FEC (p < 0.0001) and by 2.5% in patients treated with CTCb (p = 0.03) (tailored FEC vs. CTCb: 3.7%, p = 0.08). CONCLUSIONS Changes in pulmonary function were thus detected in both groups, although to a greater extent in the tailored FEC group. The clinical significance of these findings should be balanced carefully against the improved, statistically significant relapse-free survival achieved with the tailored FEC regimen compared to high-dose CTCb + peripheral blood stem cell transplantation (PSCT).

Enhancement of survival of 9L gliosarcoma bearing rats following intracerebral delivery of drugs in combination with microbeam radiation therapy.

Microbeam radiation therapy (MRT) is a form of radiosurgery first dedicated to the treatment of brain tumors. It uses arrays of synchrotron generated X-rays microbeams of very high doses (typically 625 Gy). Microbeams are typically few micrometers large (25 microm) and few hundred micrometers spaced (200 microm). Previous experiments have shown that despite a good tumor eradication rate (5/11), a 100-microm spacing unidirectional irradiation (skin dose 625 Gy, width 25 microm) was too invasive for normal tissue. On the contrary, a 200-microm spacing unidirectional irradiation preserved healthy tissue with a low tumor eradication rate (2/32). The purpose of this study was to enhance the potential of the 200 microm spacing irradiation protocol. After diagnosis of the tumor by MRI, 9L tumor-bearing rats were laterally irradiated with 51 microbeams (625 Gy, 25 microm, 200 microm) 14 days after implantation. Three drugs (Gd-DTPA, CisPt, temozolomide) were tested, after intratumoral injection at the theoretical center of the tumor. Control rats displayed a median survival time of 19 days. There was no significant difference between drug-treated rats and control group. Irradiated animals showed an increase in life span (ILS) of 60.5%. Interestingly, the ILS increased to 131.6% and 1/6 rat survived more than 1 year in case of MRT combined with gadolinium injection. These results showed that the synergy between gadolinium injection (acting as a dose enhancer) and MRT improved significantly the life span of tumor bearing rats (more than a factor 2).

Feasibility of MRI-based reference images for image-guided radiotherapy of the pelvis with either cone-beam computed tomography or planar localization images

Abstract Purpose. This study introduces methods to conduct image-guided radiotherapy (IGRT) of the pelvis with either cone-beam computed tomography (CBCT) or planar localization images by relying solely on magnetic resonance imaging (MRI)-based reference images. Material and methods. Feasibility of MRI-based reference images for IGRT was evaluated against kV CBCT (50 scans, 5 prostate cancer patients) and kV & MV planar (5 & 5 image pairs and patients) localization images by comparing the achieved patient position corrections to those obtained by standard CT-based reference images. T1/T2*-weighted in-phase MRI, Hounsfield unit conversion-based heterogeneous pseudo-CT, and bulk pseudo-CT images were applied for reference against localization CBCTs, and patient position corrections were obtained by automatic image registration. IGRT with planar localization images was performed manually by 10 observers using reference digitally reconstructed radiographs (DRRs) reconstructed from the pseudo-CTs and standard CTs. Quality of pseudo-DRRs against CT-DRRs was evaluated with image similarity metrics. Results. The SDs of differences between CBCT-to-MRI and CBCT-to-CT automatic gray-value registrations were ≤ 1.0 mm & ≤ 0.8° and ≤ 2.5 mm & ≤ 3.6° with 10 cm diameter cubic VOI and prostate-shaped VOI, respectively. The corresponding values for reference heterogeneous pseudo-CT were ≤ 1.0 mm & ≤ 0.7° and ≤ 2.2 mm & ≤ 3.3°, respectively. Heterogeneous pseudo-CT was the only type of MRI-based reference image working reliably with automatic bone registration (SDs were ≤ 0.9 mm & ≤ 0.7°). The differences include possible residual errors from planning CT to MRI registration. The image similarity metrics were significantly (p ≤ 0.01) better in agreement between heterogeneous pseudo-DRRs and CT-DRRs than between bulk pseudo-DRRs and CT-DRRs. The SDs of differences in manual registrations (3D) with planar kV and MV localization images were ≤ 1.0 mm and ≤ 1.7 mm, respectively, between heterogeneous pseudo-DRRs and CT-DRRs, and ≤ 1.4 mm and ≤ 2.1 mm between bulk pseudo-DRRs and CT-DRRs. Conclusion. This study demonstrated that it is feasible to conduct IGRT of the pelvis with MRI-based reference images.

Absorbed doses behind bones with MR image‐based dose calculations for radiotherapy treatment planning

PURPOSE Magnetic resonance (MR) images are used increasingly in external radiotherapy target delineation because of their superior soft tissue contrast compared to computed tomography (CT) images. Nevertheless, radiotherapy treatment planning has traditionally been based on the use of CT images, due to the restrictive features of MR images such as lack of electron density information. This research aimed to measure absorbed radiation doses in material behind different bone parts, and to evaluate dose calculation errors in two pseudo-CT images; first, by assuming a single electron density value for the bones, and second, by converting the electron density values inside bones from T(1)∕T(2)∗-weighted MR image intensity values. METHODS A dedicated phantom was constructed using fresh deer bones and gelatine. The effect of different bone parts to the absorbed dose behind them was investigated with a single open field at 6 and 15 MV, and measuring clinically detectable dose deviations by an ionization chamber matrix. Dose calculation deviations in a conversion-based pseudo-CT image and in a bulk density pseudo-CT image, where the relative electron density to water for the bones was set as 1.3, were quantified by comparing the calculation results with those obtained in a standard CT image by superposition and Monte Carlo algorithms. RESULTS The calculations revealed that the applied bulk density pseudo-CT image causes deviations up to 2.7% (6 MV) and 2.0% (15 MV) to the dose behind the examined bones. The corresponding values in the conversion-based pseudo-CT image were 1.3% (6 MV) and 1.0% (15 MV). The examinations illustrated that the representation of the heterogeneous femoral bone (cortex denser compared to core) by using a bulk density for the whole bone causes dose deviations up to 2% both behind the bone edge and the middle part of the bone (diameter <2.5 cm), but in the opposite directions. The measured doses and the calculated ones in the standard CT image were within 0.4% (through gelatine only) and 0.9% (behind bones). CONCLUSIONS This study indicates that the decrease in absorbed dose is not dependent on the bone diameter with all types of bones. Thus, performing dose calculation in a pseudo-CT image by assuming a single electron density value for the bones can lead to a substantial misrepresentation of the dose distribution profile. This work showed that dose calculation accuracy can be improved by using a pseudo-CT image in which the electron density values have been converted from the MR image intensity values inside bones.

Analyser-based x-ray imaging for biomedical research

Analyser-based imaging (ABI) is one of the several phase-contrast x-ray imaging techniques being pursued at synchrotron radiation facilities. With advancements in compact source technology, there is a possibility that ABI will become a clinical imaging modality. This paper presents the history of ABI as it has developed from its laboratory source to synchrotron imaging. The fundamental physics of phase-contrast imaging is presented both in a general sense and specifically for ABI. The technology is dependent on the use of perfect crystal monochromator optics. The theory of the x-ray optics is developed and presented in a way that will allow optimization of the imaging for specific biomedical systems. The advancement of analytical algorithms to produce separate images of the sample absorption, refraction angle map and small-angle x-ray scattering is detailed. Several detailed applications to biomedical imaging are presented to illustrate the broad range of systems and body sites studied preclinically to date: breast, cartilage and bone, soft tissue and organs. Ultimately, the application of ABI in clinical imaging will depend partly on the availability of compact sources with sufficient x-ray intensity comparable with that of the current synchrotron environment.

Implementation of adaptive radiation therapy for urinary bladder carcinoma: Imaging, planning and image guidance

Abstract Background. Adaptive radiation therapy (ART) for urinary bladder cancer has emerged as a promising alternative to conventional RT with potential to minimize radiation-induced toxicity to healthy tissues. In this work we have studied bladder volume variations and their effect on healthy bladder dose sparing and intrafractional margins, in order to refine our ART strategy. Material and methods. An online ART treatment strategy was followed for five patients with urinary bladder cancer with the tumors demarcated using Lipiodol®. A library of 3–4 predefined treatment plans for each patient was created based on four successive computed tomography (CT) scans. Cone beam CT (CBCT) images were acquired before each treatment fraction and after the treatment at least weekly. In partial bladder treatment the sparing of the healthy part of the bladder was investigated. The bladder wall displacements due to bladder filling were determined in three orthogonal directions (CC, AP, DEX-SIN) using the treatment planning CT scans. An ellipsoidal model was applied in order to find the theoretical maximum values for the bladder wall displacements. Moreover, the actual bladder filling rate during treatment was evaluated using the CBCT images. Results. In partial bladder treatment the volume of the bladder receiving high absorbed doses was generally smaller with a full than empty bladder. The estimation of the bladder volume and the upper limit for the intrafractional movement of the bladder wall could be represented with an ellipsoidal model with a reasonable accuracy. Observed maximum growth of bladder dimensions was less than 10 mm in all three orthogonal directions during 15 minute interval. Conclusion. The use of Lipiodol contrast agent enables partial bladder treatment with reduced irradiation of the healthy bladder volume. The ellipsoidal bladder model can be used for the estimation of the bladder volume changes and the upper limit of the bladder wall movement during the treatment fraction.

Influence of MRI-based bone outline definition errors on external radiotherapy dose calculation accuracy in heterogeneous pseudo-CT images of prostate cancer patients

Abstract Background. This work evaluates influences of susceptibility-induced bone outline shift and perturbations, and bone segmentation errors on external radiotherapy dose calculation accuracy in magnetic resonance imaging (MRI)-based pseudo-computed tomography (CT) images of the male pelvis. Material and methods. T1/T2*-weighted fast gradient echo, T1-weighted spin echo and T2-weighted fast spin echo images were used in bone detection investigation. Bone edge location and bone diameter in MRI were evaluated by comparing those in the images with actual physical measurements of fresh deer bones positioned in a gelatine phantom. Dose calculation accuracy in pseudo-CT images was investigated for 15 prostate cancer patients. Bone outlines in T1/T2*-weighted images were contoured and additional segmentation errors were simulated by expanding and contracting the bone contours with 1 mm spacing. Heterogeneous pseudo-CT images were constructed by adopting a technique transforming the MRI intensity values into Hounsfield units with separate conversion models within and outside of bone segment. Results. Bone edges and diameter in the phantom were illustrated correctly within a 1 mm-pixel size in MRI. Each 1 mm-sized systematic error in bone segment resulted in roughly 0.4% change to the prostate dose level in the pseudo-CT images. The prostate average (range) dose levels in pseudo-CT images with additional systematic bone segmentation errors of −2 mm, 0 mm and 2 mm were 0.5% (−0.5–1.4%), −0.2% (−1.0–0.7%), and −0.9% (−1.8–0.0%) compared to those in CT images, respectively, in volumetric modulated arc therapy treatment plans calculated by Monte Carlo algorithm. Conclusions. Susceptibility-induced bone outline shift and perturbations do not result in substantial uncertainty for MRI-based dose calculation. Dose consistency of 2% can be achieved reliably for the prostate if heterogeneous pseudo-CT images are constructed with ≤± 2 mm systematic error in bone segment.