Accurately simulating the production of radiotherapy portal images using non-zero beam angles

In this study, the delivery and portal imaging of one square-field and one conformal radiotherapy treatment was simulated using the Monte Carlo codes BEAMnrc and DOSXYZnrc. The treatment fields were delivered to a humanoid phantom from different angles by a 6 MV photon beam linear accelerator, with an amorphous-silicon electronic portal imaging device (a-Si EPID) used to provide images of the phantom generated by each field. The virtual phantom preparation code CTCombine was used to combine a computed-tomography-derived model of the irradiated phantom with a simple, rectilinear model of the a-Si EPID, at each beam angle used in the treatment. Comparison of the resulting experimental and simulated a-Si EPID images showed good agreement, within \[gamma](3%, 3 mm), indicating that this method may be useful in providing accurate Monte Carlo predictions of clinical a-Si EPID images, for use in the verification of complex radiotherapy treatments.

Letter to the editor : likelihood of tetraplegia after CT clearance of the cervical spine in obtunded blunt trauma patients

...the probabilistic computer simulation study by Dunham and colleagues evaluating the impact of different cervical spine management (CSM) strategies on tetraplegia and brain injury outcomes.1 Based on literature findings, expert opinion and with use of advances programming techniques the authors conclude that early collar removal without cervical spine magnetic resonance imaging (MRI) is a preferable CSM strategy for comatose, blunt trauma patients with extremity movement and a negative cervical spine computed tomography(CT) scan. Although we do not have the required expertise to comment on the applied statistical approach, we would like to comment on one of the medical assumptions raised by the authors, namely the likelihood of tetraplegia in this specific population....

Post-operative CT assessment of interbody fusion two years after thoracoscopic scoliosis surgery

The relationship between radiologic union and clinical outcome in thoracoscopic scoliosis surgery is not clear, as apparent non-union of a spinal fusion does not always correspond to a poor clinical result. The aim of this study was to evaluate CT fusion rates 24 months after thoracoscopic anterior scoliosis surgery, and to explore the relationship between fusion scores and; (i) rod diameter, (ii) graft type, (iii) fusion level, (iv) occurrence of post-operative implant failure, and (v) lateral position of the fusion mass in the intervertebral disc space. We propose that moderate fusion scores on the Sucato scale secure successful clinical outcomes in thoracoscopic scoliosis surgery.

Postoperative low dose CT assessment of interbody fusion two years after thoracoscopic scoliosis surgery

The relationship between radiologic union and clinical outcomes in thoracoscopic scoliosis surgery is not clear, as apparent non-union of a spinal fusion does not always correspond to a poor clinical result. The aim of this study was to evaluate for the first time the interbody fusion rates using low dose CT scans at minimum 24 months after thoracoscopic scoliosis surgery, and to explore the relationship between fusion scores and; (i) rod diameter, (ii) graft type, (iii) fusion level, (iv) implant failure, and (v) lateral position in the disc space. The study found that moderate fusion scores on the Sucato scale secure successful clinical outcomes in thoracoscopic scoliosis surgery.

Post-operative CT assessment of interbody fusion two years after thoracoscopic scoliosis surgery

The relationship between radiologic union and clinical outcome in thoracoscopic scoliosis surgery is not clear, as apparent non-union of a spinal fusion does not always correspond to a poor clinical result. The aim of this study was to evaluate CT fusion rates 2yrs after thoracoscopic surgery, and to explore the relationship between fusion scores and rod diameter, graft type, fusion level, implant failure, and lateral position in the disc space. This study suggests that moderate fusion scores secure successful clinical outcomes in thoracoscopic scoliosis surgery.

A subject-specific model of human buttocks and thighs in a seated posture

Finite element analyses of the human body in seated postures requires digital models capable of providing accurate and precise prediction of the tissue-level response of the body in the seated posture. To achieve such models, the human anatomy must be represented with high fidelity. This information can readily be defined using medical imaging techniques such as Magnetic Resonance Imaging (MRI) or Computed Tomography (CT). Current practices for constructing digital human models, based on the magnetic resonance (MR) images, in a lying down (supine) posture have reduced the error in the geometric representation of human anatomy relative to reconstructions based on data from cadaveric studies. Nonetheless, the significant differences between seated and supine postures in segment orientation, soft-tissue deformation and soft tissue strain create a need for data obtained in postures more similar to the application posture. In this study, we present a novel method for creating digital human models based on seated MR data. An adult-male volunteer was scanned in a simulated driving posture using a FONAR 0.6T upright MRI scanner with a T1 scanning protocol. To compensate for unavoidable image distortion near the edges of the study, images of the same anatomical structures were obtained in transverse and sagittal planes. Combinations of transverse and sagittal images were used to reconstruct the major anatomical features from the buttocks through the knees, including bone, muscle and fat tissue perimeters, using Solidworks® software. For each MR image, B-splines were created as contours for the anatomical structures of interest, and LOFT commands were used to interpolate between the generated Bsplines. The reconstruction of the pelvis, from MR data, was enhanced by the use of a template model generated in previous work CT images. A non-rigid registration algorithm was used to fit the pelvis template into the MR data. Additionally, MR image processing was conducted to both the left and the right sides of the model due to the intended asymmetric posture of the volunteer during the MR measurements. The presented subject-specific, three-dimensional model of the buttocks and thighs will add value to optimisation cycles in automotive seat development when used in simulating human interaction with automotive seats.

3D reconstruction of long bones utilising magnetic resonance imaging (MRI)

The design of pre-contoured fracture fixation implants (plates and nails) that correctly fit the anatomy of a patient utilises 3D models of long bones with accurate geometric representation. 3D data is usually available from computed tomography (CT) scans of human cadavers that generally represent the above 60 year old age group. Thus, despite the fact that half of the seriously injured population comes from the 30 year age group and below, virtually no data exists from these younger age groups to inform the design of implants that optimally fit patients from these groups. Hence, relevant bone data from these age groups is required. The current gold standard for acquiring such data–CT–involves ionising radiation and cannot be used to scan healthy human volunteers. Magnetic resonance imaging (MRI) has been shown to be a potential alternative in the previous studies conducted using small bones (tarsal bones) and parts of the long bones. However, in order to use MRI effectively for 3D reconstruction of human long bones, further validations using long bones and appropriate reference standards are required. Accurate reconstruction of 3D models from CT or MRI data sets requires an accurate image segmentation method. Currently available sophisticated segmentation methods involve complex programming and mathematics that researchers are not trained to perform. Therefore, an accurate but relatively simple segmentation method is required for segmentation of CT and MRI data. Furthermore, some of the limitations of 1.5T MRI such as very long scanning times and poor contrast in articular regions can potentially be reduced by using higher field 3T MRI imaging. However, a quantification of the signal to noise ratio (SNR) gain at the bone - soft tissue interface should be performed; this is not reported in the literature. As MRI scanning of long bones has very long scanning times, the acquired images are more prone to motion artefacts due to random movements of the subject‟s limbs. One of the artefacts observed is the step artefact that is believed to occur from the random movements of the volunteer during a scan. This needs to be corrected before the models can be used for implant design. As the first aim, this study investigated two segmentation methods: intensity thresholding and Canny edge detection as accurate but simple segmentation methods for segmentation of MRI and CT data. The second aim was to investigate the usability of MRI as a radiation free imaging alternative to CT for reconstruction of 3D models of long bones. The third aim was to use 3T MRI to improve the poor contrast in articular regions and long scanning times of current MRI. The fourth and final aim was to minimise the step artefact using 3D modelling techniques. The segmentation methods were investigated using CT scans of five ovine femora. The single level thresholding was performed using a visually selected threshold level to segment the complete femur. For multilevel thresholding, multiple threshold levels calculated from the threshold selection method were used for the proximal, diaphyseal and distal regions of the femur. Canny edge detection was used by delineating the outer and inner contour of 2D images and then combining them to generate the 3D model. Models generated from these methods were compared to the reference standard generated using the mechanical contact scans of the denuded bone. The second aim was achieved using CT and MRI scans of five ovine femora and segmenting them using the multilevel threshold method. A surface geometric comparison was conducted between CT based, MRI based and reference models. To quantitatively compare the 1.5T images to the 3T MRI images, the right lower limbs of five healthy volunteers were scanned using scanners from the same manufacturer. The images obtained using the identical protocols were compared by means of SNR and contrast to noise ratio (CNR) of muscle, bone marrow and bone. In order to correct the step artefact in the final 3D models, the step was simulated in five ovine femora scanned with a 3T MRI scanner. The step was corrected using the iterative closest point (ICP) algorithm based aligning method. The present study demonstrated that the multi-threshold approach in combination with the threshold selection method can generate 3D models from long bones with an average deviation of 0.18 mm. The same was 0.24 mm of the single threshold method. There was a significant statistical difference between the accuracy of models generated by the two methods. In comparison, the Canny edge detection method generated average deviation of 0.20 mm. MRI based models exhibited 0.23 mm average deviation in comparison to the 0.18 mm average deviation of CT based models. The differences were not statistically significant. 3T MRI improved the contrast in the bone–muscle interfaces of most anatomical regions of femora and tibiae, potentially improving the inaccuracies conferred by poor contrast of the articular regions. Using the robust ICP algorithm to align the 3D surfaces, the step artefact that occurred by the volunteer moving the leg was corrected, generating errors of 0.32 ± 0.02 mm when compared with the reference standard. The study concludes that magnetic resonance imaging, together with simple multilevel thresholding segmentation, is able to produce 3D models of long bones with accurate geometric representations. The method is, therefore, a potential alternative to the current gold standard CT imaging.

Computer planning of stereotactic iodine-125 seed brachytherapy for recurrent malignant gliomas

At St Thomas' Hospital, we have developed a computer program on a Titan graphics supercomputer to plan the stereotactic implantation of iodine-125 seeds for the palliative treatment of recurrent malignant gliomas. Use of the Gill-Thomas-Cosman relocatable frame allows planning and surgery to be carried out at different hospitals on different days. Stereotactic computed tomography (CT) and positron emission tomography (PET) scans are performed and the images transferred to the planning computer. The head, tumour and frame fiducials are outlined on the relevant images, and a three-dimensional model generated. Structures which could interfere with the surgery or radiotherapy, such as major vessels, shunt tubing etc., can also be outlined and included in the display. Catheter target and entry points are set using a three-dimensional cursor controlled by a set of dials attached to the computer. The program calculates and displays the radiation dose distribution within the target volume for various catheter and seed arrangements. The CT co-ordinates of the fiducial rods are used to convert catheter co-ordinates from CT space to frame space and to calculate the catheter insertion angles and depths. The surgically implanted catheters are after-loaded the next day and the seeds left in place for between 4 and 6 days, giving a nominal dose of 50 Gy to the edge of the target volume. 25 patients have been treated so far.

Tooth fracture risk analysis based on a new finite element dental structure models using micro-CT data

The finite element (FE) analysis is an effective method to study the strength and predict the fracture risk of endodontically-treated teeth. This paper presents a rapid method developed to generate a comprehensive tooth FE model using data retrieved from micro-computed tomography (μCT). With this method, the inhomogeneity of material properties of teeth was included into the model without dividing the tooth model into different regions. The material properties of the tooth were assumed to be related to the mineral density. The fracture risk at different tooth portions was assessed for root canal treatments. The micro-CT images of a tooth were processed by a Matlab software programme and the CT numbers were retrieved. The tooth contours were obtained with thresholding segmentation using Amira. The inner and outer surfaces of the tooth were imported into Solidworks and a three-dimensional (3D) tooth model was constructed. An assembly of the tooth model with the periodontal ligament (PDL) layer and surrounding bone was imported into ABAQUS. The material properties of the tooth were calculated from the retrieved CT numbers via ABAQUS user's subroutines. Three root canal geometries (original and two enlargements) were investigated. The proposed method in this study can generate detailed 3D finite element models of a tooth with different root canal enlargements and filling materials, and would be very useful for the assessment of the fracture risk at different tooth portions after root canal treatments.

Towards determining soft tissue properties for modelling spine surgery : current progress and challenges

Current complication rates for adolescent scoliosis surgery necessitate the development of better surgical planning tools to improve outcomes. Here we present our approach to developing finite element models of the thoracolumbar spine for deformity surgery simulation, with patient-specific model anatomy based on low-dose pre-operative computed tomography scans. In a first step towards defining patient-specific tissue properties, an initial 'benchmark' set of properties were used to simulate a clinically performed pre-operative spinal flexibility assessment, the fulcrum bending radiograph. Clinical data for ten patients were compared with the simulated results for this assessment and in cases where these data differed by more than 10%, soft tissue properties for the costo-vertebral joint (CVJt) were altered to achieve better agreement. Results from these analyses showed that changing the CVJt stiffness resulted in acceptable agreement between clinical and simulated flexibility in two of the six cases. In light of these results and those of our previous studies in this area, it is suggested that spinal flexibility in the fulcrum bending test is not governed by any single soft tissue structure acting in isolation. More detailed biomechanical characterisation of the fulcrum bending test is required to provide better data for determination of patient-specific soft tissue properties.

Zoledronic acid preserves bone structure and increases survival but does not limit tumour incidence in a prostate cancer bone metastasis model

Background The bisphosphonate, zoledronic acid (ZOL), can inhibit osteoclasts leading to decreased osteoclastogenesis and osteoclast activity in bone. Here, we used a mixed osteolytic/osteoblastic murine model of bone-metastatic prostate cancer, RM1(BM), to determine how inhibiting osteolysis with ZOL affects the ability of these cells to establish metastases in bone, the integrity of the tumour-bearing bones and the survival of the tumour-bearing mice. Methods The model involves intracardiac injection for arterial dissemination of the RM1(BM) cells in C57BL/6 mice. ZOL treatment was given via subcutaneous injections on days 0, 4, 8 and 12, at 20 and 100 µg/kg doses. Bone integrity was assessed by micro-computed tomography and histology with comparison to untreated mice. The osteoclast and osteoblast activity was determined by measuring serum tartrate-resistant acid phosphatase 5b (TRAP 5b) and osteocalcin, respectively. Mice were euthanased according to predetermined criteria and survival was assessed using Kaplan Meier plots. Findings Micro-CT and histological analysis showed that treatment of mice with ZOL from the day of intracardiac injection of RM1(BM) cells inhibited tumour-induced bone lysis, maintained bone volume and reduced the calcification of tumour-induced endochondral osteoid material. ZOL treatment also led to a decreased serum osteocalcin and TRAP 5b levels. Additionally, treated mice showed increased survival compared to vehicle treated controls. However, ZOL treatment did not inhibit the cells ability to metastasise to bone as the number of bone-metastases was similar in both treated and untreated mice. Conclusions ZOL treatment provided significant benefits for maintaining the integrity of tumour-bearing bones and increased the survival of tumour bearing mice, though it did not prevent establishment of bone-metastases in this model. From the mechanistic view, these observations confirm that tumour-induced bone lysis is not a requirement for establishment of these bone tumours.

CT and radiographic analysis of sagittal profile changes after thoracoscopic anterior scoliosis surgery

Background. Previous studies report an increase in thoracic kyphosis after anterior approaches and a flattening of sagittal contours following posterior approaches. Difficulties with measuring sagittal parameters on radiographs are avoided with reformatted sagittal CT reconstructions due to the superior endplate clarity afforded by this imaging modality. Methods. A prospective study of 30 Lenke 1 adolescent idiopathic scoliosis (AIS) patients receiving selective thoracoscopic anterior spinal fusion (TASF) was performed. Participants had ethically approved low dose CT scans at minimum 24 months after surgery in addition to their standard care following surgery. The change in sagittal contours on supine CT was compared to standing radiographic measurements of the same patients and with previous studies. Inter-observer variability was assessed as well as whether hypokyphotic and normokyphotic patient groups responded differently to the thoracoscopic anterior approach. Results. Mean T5-12 kyphosis Cobb angle increased by 11.8 degrees and lumbar lordosis increased by 5.9 degrees on standing radiographs two years after surgery. By comparison, CT measurements of kyphosis and lordosis increased by 12.3 degrees and 7.0 degrees respectively. 95% confidence intervals for inter-observer variability of sagittal contour measurements on supine CT ranged between 5-8 degrees. TASF had a slightly greater corrective effect on patients who were hypokyphotic before surgery compared with those who were normokyphotic. Conclusions. Restoration of sagittal profile is an important goal of scoliosis surgery, but reliable measurement with radiographs suffers from poor endplate clarity. TASF significantly improves thoracic kyphosis and lumbar lordosis while preserving proximal and distal junctional alignment in thoracic AIS patients. Supine CT allows greater endplate clarity for sagittal Cobb measurements and linear relationships were found between supine CT and standing radiographic measurements. In this study, improvements in sagittal kyphosis and lordosis following surgery were in agreement with prior anterior surgery studies, and add to the current evidence suggesting that anterior correction is more capable than posterior approaches of addressing the sagittal component of both the instrumented and adjacent non instrumented segments following surgical correction of progressive Lenke 1 idiopathic scoliosis.