Breathing adapted radiotherapy: a 4D gating software for lung cancer

Purpose Physiological respiratory motion of tumors growing in the lung can be corrected with respiratory gating when treated with radiotherapy (RT). The optimal respiratory phase for beam-on may be assessed with a respiratory phase optimizer (RPO), a 4D image processing software developed with this purpose. Methods and Materials Fourteen patients with lung cancer were included in the study. Every patient underwent a 4D-CT providing ten datasets of ten phases of the respiratory cycle (0-100% of the cycle). We defined two morphological parameters for comparison of 4D-CT images in different respiratory phases: tumor-volume to lung-volume ratio and tumor-to-spinal cord distance. The RPO automatized the calculations (200 per patient) of these parameters for each phase of the respiratory cycle allowing to determine the optimal interval for RT. Results Lower lobe lung tumors not attached to the diaphragm presented with the largest motion with breathing. Maximum inspiration was considered the optimal phase for treatment in 4 patients (28.6%). In 7 patients (50%), however, the RPO showed a most favorable volumetric and spatial configuration in phases other than maximum inspiration. In 2 cases (14.4%) the RPO showed no benefit from gating. This tool was not conclusive in only one case. Conclusions The RPO software presented in this study can help to determine the optimal respiratory phase for gated RT based on a few simple morphological parameters. Easy to apply in daily routine, it may be a useful tool for selecting patients who might benefit from breathing adapted RT.

Irregularities below the Enamel-Dentin Junction May Predispose for Fissure Caries

This study investigated the structure of the fissure fundus on occlusal surfaces with respect to the detection of possible irregularities below the enamel-dentin junction (EDJ). Occlusal surfaces were examined by micro-computed tomography (µCT). In total, 203 third molars with clinically sound occlusal fissures or non-cavitated lesions were selected. All specimens were scanned with µCT. Subsequently, each tooth was sectioned, and each slice was investigated by stereomicroscopy. In 7 of 203 molars (3.4%), demarcated radiolucencies below the EDJ were detected by µCT. These defects were obviously of non-carious origin, because the µCT images revealed no gradient of demineralization in the dentin. In all cases, a direct pathway between the oral cavity and the dentin was evident. The comparison of the µCT sites with conventional histological images also revealed defects in the dentin. These results demonstrate that demarcated radiolucencies below the EDJ may not necessarily be caries lesions according to µCT images and may be classified as possible developmental irregularities. To avoid misinterpreting µCT data, dental researchers should carefully consider this condition when analyzing µCT images. The clinical significance of this finding is that these defects may predispose molar teeth to early-onset caries in occlusal pits and fissures.

Simultaneous Multiscale Polyaffine Registration by Incorporating Deformation Statistics

Locally affine (polyaffine) image registration methods capture intersubject non-linear deformations with a low number of parameters, while providing an intuitive interpretation for clinicians. Considering the mandible bone, anatomical shape differences can be found at different scales, e.g. left or right side, teeth, etc. Classically, sequential coarse to fine registration are used to handle multiscale deformations, instead we propose a simultaneous optimization of all scales. To avoid local minima we incorporate a prior on the polyaffine transformations. This kind of groupwise registration approach is natural in a polyaffine context, if we assume one configuration of regions that describes an entire group of images, with varying transformations for each region. In this paper, we reformulate polyaffine deformations in a generative statistical model, which enables us to incorporate deformation statistics as a prior in a Bayesian setting. We find optimal transformations by optimizing the maximum a posteriori probability. We assume that the polyaffine transformations follow a normal distribution with mean and concentration matrix. Parameters of the prior are estimated from an initial coarse to fine registration. Knowing the region structure, we develop a blockwise pseudoinverse to obtain the concentration matrix. To our knowledge, we are the first to introduce simultaneous multiscale optimization through groupwise polyaffine registration. We show results on 42 mandible CT images.

High-accuracy patient-to-image registration for the facilitation of image-guided robotic microsurgery on the head

Image-guided microsurgery requires accuracies an order of magnitude higher than today's navigation systems provide. A critical step toward the achievement of such low-error requirements is a highly accurate and verified patient-to-image registration. With the aim of reducing target registration error to a level that would facilitate the use of image-guided robotic microsurgery on the rigid anatomy of the head, we have developed a semiautomatic fiducial detection technique. Automatic force-controlled localization of fiducials on the patient is achieved through the implementation of a robotic-controlled tactile search within the head of a standard surgical screw. Precise detection of the corresponding fiducials in the image data is realized using an automated model-based matching algorithm on high-resolution, isometric cone beam CT images. Verification of the registration technique on phantoms demonstrated that through the elimination of user variability, clinically relevant target registration errors of approximately 0.1 mm could be achieved.

In Vitro and In Vivo Imaging Characteristics Assessment of Polymeric Coils Compared with Standard Platinum Coils for the Treatment of Intracranial Aneurysms

BACKGROUND AND PURPOSE:Conventional platinum coils cause imaging artifacts that reduce imaging quality and therefore impair imaging interpretation on intraprocedural or noninvasive follow-up imaging. The purpose of this study was to evaluate imaging characteristics and artifact production of polymeric coils compared with standard platinum coils in vitro and in vivo.MATERIALS AND METHODS:Polymeric coils and standard platinum coils were evaluated in vitro with the use of 2 identical silicon aneurysm models coiled with a packing attenuation of 20% each. DSA, flat panel CT, CT, and MR imaging were performed. In vivo evaluation of imaging characteristics of polymeric coils was performed in experimentally created rabbit carotid bifurcation aneurysms. DSA, CT/CTA, and MR imaging were performed after endovascular treatment of the aneurysms. Images were evaluated regarding visibility of individual coils, coil mass, artifact production, and visibility of residual flow within the aneurysm.RESULTS:Overall, in vitro and in vivo imaging showed relevantly reduced artifact production of polymeric coils in all imaging modalities compared with standard platinum coils. Image quality of CT and MR imaging was improved with the use of polymeric coils, which permitted enhanced depiction of individual coil loops and residual aneurysm lumen as well as the peri-aneurysmal area. Remarkably, CT images demonstrated considerably improved image quality with only minor artifacts compared with standard coils. On DSA, polymeric coils showed transparency and allowed visualization of superimposed vessel structures.CONCLUSIONS:This initial experimental study showed improved imaging quality with the use of polymeric coils compared with standard platinum coils in all imaging modalities. This might be advantageous for improved intraprocedural imaging for the detection of complications and posttreatment noninvasive follow-up imaging.

Effects of anatomic conformation on three-dimensional motion of the caudal lumbar and lumbosacral portions of the vertebral column of dogs

OBJECTIVE: To determine the association between the 3-dimensional (3-D) motion pattern of the caudal lumbar and lumbosacral portions of the canine vertebral column and the morphology of vertebrae, facet joints, and intervertebral disks. SAMPLE POPULATION: Vertebral columns of 9 German Shepherd Dogs and 16 dogs of other breeds with similar body weights and body conditions. PROCEDURE: Different morphometric parameters of the vertebral column were assessed by computed tomography (CT) and magnetic resonance imaging. Anatomic conformation and the 3-D motion pattern were compared, and correlation coefficients were calculated. RESULTS: Total range of motion for flexion and extension was mainly associated with the facet joint angle, the facet joint angle difference between levels of the vertebral column in the transverse plane on CT images, disk height, and lever arm length. CONCLUSIONS AND CLINICAL RELEVANCE: Motion is a complex process that is influenced by the entire 3-D conformation of the lumbar portion of the vertebral column. In vivo dynamic measurements of the 3-D motion pattern of the lumbar and lumbosacral portions of the vertebral column will be necessary to further assess biomechanics that could lead to disk degeneration in dogs.

Effectively Selecting a Target Population for a Future Comparative Study

When comparing a new treatment with a control in a randomized clinical study, the treatment effect is generally assessed by evaluating a summary measure over a specific study population. The success of the trial heavily depends on the choice of such a population. In this paper, we show a systematic, effective way to identify a promising population, for which the new treatment is expected to have a desired benefit, using the data from a current study involving similar comparator treatments. Specifically, with the existing data we first create a parametric scoring system using multiple covariates to estimate subject-specific treatment differences. Using this system, we specify a desired level of treatment difference and create a subgroup of patients, defined as those whose estimated scores exceed this threshold. An empirically calibrated group-specific treatment difference curve across a range of threshold values is constructed. The population of patients with any desired level of treatment benefit can then be identified accordingly. To avoid any ``self-serving'' bias, we utilize a cross-training-evaluation method for implementing the above two-step procedure. Lastly, we show how to select the best scoring system among all competing models. The proposals are illustrated with the data from two clinical trials in treating AIDS and cardiovascular diseases. Note that if we are not interested in designing a new study for comparing similar treatments, the new procedure can also be quite useful for the management of future patients who would receive nontrivial benefits to compensate for the risk or cost of the new treatment.

Analyzing Panel Count Data with Informative Observation Times

In this paper, we study panel count data with informative observation times. We assume nonparametric and semiparametric proportional rate models for the underlying recurrent event process, where the form of the baseline rate function is left unspecified and a subject-specific frailty variable inflates or deflates the rate function multiplicatively. The proposed models allow the recurrent event processes and observation times to be correlated through their connections with the unobserved frailty; moreover, the distributions of both the frailty variable and observation times are considered as nuisance parameters. The baseline rate function and the regression parameters are estimated by maximizing a conditional likelihood function of observed event counts and solving estimation equations. Large sample properties of the proposed estimators are studied. Numerical studies demonstrate that the proposed estimation procedures perform well for moderate sample sizes. An application to a bladder tumor study is presented to illustrate the use of the proposed methods.

Hardware accelerated ray cast of volume data and volume gradient for an optimized splines-based multi-resolution 2D-3D registration

This paper describes a method for DRR generation as well as for volume gradients projection using hardware accelerated 2D texture mapping and accumulation buffering and demonstrates its application in 2D-3D registration of X-ray fluoroscopy to CT images. The robustness of the present registration scheme are guaranteed by taking advantage of a coarse-to-fine processing of the volume/image pyramids based on cubic B-splines. A human cadaveric spine specimen together with its ground truth was used to compare the present scheme with a purely software-based scheme in three aspects: accuracy, speed, and capture ranges. Our experiments revealed an equivalent accuracy and capture ranges but with much shorter registration time with the present scheme. More specifically, the results showed 0.8 mm average target registration error, 55 second average execution time per registration, and 10 mm and 10° capture ranges for the present scheme when tested on a 3.0 GHz Pentium 4 computer.

Pathomorphologic alterations predict presence or absence of hip osteoarthrosis

Abnormal morphology of the hip has been associated with primary osteoarthrosis. We evaluated the morphology of 464 consecutive hips contralateral to hips treated by THA. We excluded all hips with known diagnoses leading to secondary osteoarthritis and all hips with advanced arthrosis to eliminate the effect of arthritic remodeling on the morphologic measurements. Of the remaining 119 hips, 25 were in patients aged 60 years or older who had no or mild arthrosis (Tönnis Grade 0 or 1) and 94 hips had Tönnis Grade 2 osteoarthrosis. We quantified morphologic parameters on plain radiographs and CT images and simulated range of motion using virtual bone models from the CT data. The nonarthritic hips had fewer pathomorphologic findings. High alpha angles and high lateral center edge angles were strongly associated with the presence of arthritis; decreased internal and external rotation in 90 degrees flexion showed lesser correlation. The data confirm previous observations that abnormal hip morphology predates arthrosis and is not secondary to the osteoarthritic process. Hips at risk for developing arthrosis resulting from pathomorphologic changes may potentially be identified at the cessation of growth, long before the development of osteoarthrosis.

Haptics in forensics: the possibilities and advantages in using the haptic device for reconstruction approaches in forensic science

Non-invasive documentation methods such as surface scanning and radiological imaging are gaining in importance in the forensic field. These three-dimensional technologies provide digital 3D data, which are processed and handled in the computer. However, the sense of touch gets lost using the virtual approach. The haptic device enables the use of the sense of touch to handle and feel digital 3D data. The multifunctional application of a haptic device for forensic approaches is evaluated and illustrated in three different cases: the representation of bone fractures of the lower extremities, by traffic accidents, in a non-invasive manner; the comparison of bone injuries with the presumed injury-inflicting instrument; and in a gunshot case, the identification of the gun by the muzzle imprint, and the reconstruction of the holding position of the gun. The 3D models of the bones are generated from the Computed Tomography (CT) images. The 3D models of the exterior injuries, the injury-inflicting tools and the bone injuries, where a higher resolution is necessary, are created by the optical surface scan. The haptic device is used in combination with the software FreeForm Modelling Plus for touching the surface of the 3D models to feel the minute injuries and the surface of tools, to reposition displaced bone parts and to compare an injury-causing instrument with an injury. The repositioning of 3D models in a reconstruction is easier, faster and more precisely executed by means of using the sense of touch and with the user-friendly movement in the 3D space. For representation purposes, the fracture lines of bones are coloured. This work demonstrates that the haptic device is a suitable and efficient application in forensic science. The haptic device offers a new way in the handling of digital data in the virtual 3D space.

Effect of chest compressions only during experimental basic life support on alveolar collapse and recruitment

AIM: The importance of ventilatory support during cardiac arrest and basic life support is controversial. This experimental study used dynamic computed tomography (CT) to assess the effects of chest compressions only during cardiopulmonary resuscitation (CCO-CPR) on alveolar recruitment and haemodynamic parameters in porcine model of ventricular fibrillation. MATERIALS AND METHODS: Twelve anaesthetized pigs (26+/-1kg) were randomly assigned to one of the following groups: (1) intermittent positive pressure ventilation (IPPV) both during basic life support and advanced cardiac life support, or (2) CCO during basic life support and IPPV during advanced cardiac life support. Measurements were acquired at baseline prior to cardiac arrest, during basic life support, during advanced life support, and after return of spontaneous circulation (ROSC), as follows: dynamic CT series, arterial and central venous pressures, blood gases, and regional organ blood flow. The ventilated and atelectatic lung area was quantified from dynamic CT images. Differences between groups were analyzed using the Kruskal-Wallis test, and a p<0.05 was considered statistically significant. RESULTS: IPPV was associated with cyclic alveolar recruitment and de-recruitment. Compared with controls, the CCO-CPR group had a significantly larger mean fractional area of atelectasis (p=0.009), and significantly lower PaO(2) (p=0.002) and mean arterial pressure (p=0.023). The increase in mean atelectatic lung area observed during basic life support in the CCO-CPR group remained clinically relevant throughout the subsequent advanced cardiac life support period and following ROSC, and was associated with prolonged impaired haemodynamics. No inter-group differences in myocardial and cerebral blood flow were observed. CONCLUSION: A lack of ventilation during basic life support is associated with excessive atelectasis, arterial hypoxaemia and compromised CPR haemodynamics. Moreover, these detrimental effects remain evident even after restoration of IPPV.

An integral conservative gridding--algorithm using Hermitian curve interpolation

The problem of re-sampling spatially distributed data organized into regular or irregular grids to finer or coarser resolution is a common task in data processing. This procedure is known as 'gridding' or 're-binning'. Depending on the quantity the data represents, the gridding-algorithm has to meet different requirements. For example, histogrammed physical quantities such as mass or energy have to be re-binned in order to conserve the overall integral. Moreover, if the quantity is positive definite, negative sampling values should be avoided. The gridding process requires a re-distribution of the original data set to a user-requested grid according to a distribution function. The distribution function can be determined on the basis of the given data by interpolation methods. In general, accurate interpolation with respect to multiple boundary conditions of heavily fluctuating data requires polynomial interpolation functions of second or even higher order. However, this may result in unrealistic deviations (overshoots or undershoots) of the interpolation function from the data. Accordingly, the re-sampled data may overestimate or underestimate the given data by a significant amount. The gridding-algorithm presented in this work was developed in order to overcome these problems. Instead of a straightforward interpolation of the given data using high-order polynomials, a parametrized Hermitian interpolation curve was used to approximate the integrated data set. A single parameter is determined by which the user can control the behavior of the interpolation function, i.e. the amount of overshoot and undershoot. Furthermore, it is shown how the algorithm can be extended to multidimensional grids. The algorithm was compared to commonly used gridding-algorithms using linear and cubic interpolation functions. It is shown that such interpolation functions may overestimate or underestimate the source data by about 10-20%, while the new algorithm can be tuned to significantly reduce these interpolation errors. The accuracy of the new algorithm was tested on a series of x-ray CT-images (head and neck, lung, pelvis). The new algorithm significantly improves the accuracy of the sampled images in terms of the mean square error and a quality index introduced by Wang and Bovik (2002 IEEE Signal Process. Lett. 9 81-4).

Validation of a new method for determination of cup orientation in THA

Our goal was to validate accuracy, consistency, and reproducibility/reliability of a new method for determining cup orientation in total hip arthroplasty (THA). This method allows matching the 3D-model from CT images or slices with the projected pelvis on an anteroposterior pelvic radiograph using a fully automated registration procedure. Cup orientation (inclination and anteversion) is calculated relative to the anterior pelvic plane, corrected for individual malposition of the pelvis during radiograph acquisition. Measurements on blinded and randomized radiographs of 80 cadaver and 327 patient hips were investigated. The method showed a mean accuracy of 0.7 +/- 1.7 degrees (-3.7 degrees to 4.0 degrees) for inclination and 1.2 +/- 2.4 degrees (-5.3 degrees to 5.6 degrees) for anteversion in the cadaver trials and 1.7 +/- 1.7 degrees (-4.6 degrees to 5.5 degrees) for inclination and 0.9 +/- 2.8 degrees (-5.2 degrees to 5.7 degrees) for anteversion in the clinical data when compared to CT-based measurements. No systematic errors in accuracy were detected with the Bland-Altman analysis. The software consistency and the reproducibility/reliability were very good. This software is an accurate, consistent, reliable, and reproducible method to measure cup orientation in THA using a sophisticated 2D/3D-matching technique. Its robust and accurate matching algorithm can be expanded to statistical models.

A complementary radiographic projection of the equine temporomandibular joint

The complexity of the equine skull makes the temporomandibular joint a difficult area to evaluate radiographically. The goal of this study was to determine the optimal angle for a complementary radiographic projection of the equine temporomandibular joint based on a computed tomography (CT) cadaver study. CT was performed on six equine cadaver heads of horses that were euthanized for other reasons than temporomandibular joint disease. After the CT examination, 3D reconstruction of the equine skull was performed to subjectively determine the angle for a complementary radiographic projection of the temporomandibular joint. The angle was measured on the left and right temporomandibular joint of each head. Based on the measurements obtained from the CT images, a radiographic projection of the temporomandibular joint in a rostra-145 degrees ventral-caudodorsal oblique (R45 degrees V-CdDO) direction was developed by placing the X-ray unit 30 degrees laterally, maintaining at the same time the R45 degrees V-CdDO angle (R45 degrees V30 degrees L-CdDLO). This radiographic projection was applied to all cadaver heads and on six live horses. In three of the live horses abnormal findings associated with the temporomandibular joint were detected. We conclude that this new radiographic projection of the temporomandibular joint provides superior visualization of the temporomandibular joint space and the articular surface of the mandibular condyle.