Kernel density correlation based non-rigid point set matching for statistical model-based 2D/3D reconstruction

This paper presents a kernel density correlation based nonrigid point set matching method and shows its application in statistical model based 2D/3D reconstruction of a scaled, patient-specific model from an un-calibrated x-ray radiograph. In this method, both the reference point set and the floating point set are first represented using kernel density estimates. A correlation measure between these two kernel density estimates is then optimized to find a displacement field such that the floating point set is moved to the reference point set. Regularizations based on the overall deformation energy and the motion smoothness energy are used to constraint the displacement field for a robust point set matching. Incorporating this non-rigid point set matching method into a statistical model based 2D/3D reconstruction framework, we can reconstruct a scaled, patient-specific model from noisy edge points that are extracted directly from the x-ray radiograph by an edge detector. Our experiment conducted on datasets of two patients and six cadavers demonstrates a mean reconstruction error of 1.9 mm

Computer algebra for x-ray spectral reconstruction between 6 and 25 MV

A previously presented algorithm for the reconstruction of bremsstrahlung spectra from transmission data has been implemented into MATHEMATICA. Spectra vectorial algebra has been used to solve the matrix system A * F = T. The new implementation has been tested by reconstructing photon spectra from transmission data acquired in narrow beam conditions, for nominal energies of 6, 15, and 25 MV. The results were in excellent agreement with the original calculations. Our implementation has the advantage to be based on a well-tested mathematical kernel. Furthermore it offers a comfortable user interface.

Accurate and robust reconstruction of proximal femur from sparse intraoperative data and dense point distribution model for surgical navigation

Constructing a 3D surface model from sparse-point data is a nontrivial task. Here, we report an accurate and robust approach for reconstructing a surface model of the proximal femur from sparse-point data and a dense-point distribution model (DPDM). The problem is formulated as a three-stage optimal estimation process. The first stage, affine registration, is to iteratively estimate a scale and a rigid transformation between the mean surface model of the DPDM and the sparse input points. The estimation results of the first stage are used to establish point correspondences for the second stage, statistical instantiation, which stably instantiates a surface model from the DPDM using a statistical approach. This surface model is then fed to the third stage, kernel-based deformation, which further refines the surface model. Handling outliers is achieved by consistently employing the least trimmed squares (LTS) approach with a roughly estimated outlier rate in all three stages. If an optimal value of the outlier rate is preferred, we propose a hypothesis testing procedure to automatically estimate it. We present here our validations using four experiments, which include 1 leave-one-out experiment, 2 experiment on evaluating the present approach for handling pathology, 3 experiment on evaluating the present approach for handling outliers, and 4 experiment on reconstructing surface models of seven dry cadaver femurs using clinically relevant data without noise and with noise added. Our validation results demonstrate the robust performance of the present approach in handling outliers, pathology, and noise. An average 95-percentile error of 1.7-2.3 mm was found when the present approach was used to reconstruct surface models of the cadaver femurs from sparse-point data with noise added.

Reducing CT radiation dose with iterative reconstruction algorithms: the influence of scan and reconstruction parameters on image quality and CTDIvol

OBJECTIVES In this phantom CT study, we investigated whether images reconstructed using filtered back projection (FBP) and iterative reconstruction (IR) with reduced tube voltage and current have equivalent quality. We evaluated the effects of different acquisition and reconstruction parameter settings on image quality and radiation doses. Additionally, patient CT studies were evaluated to confirm our phantom results. METHODS Helical and axial 256 multi-slice computed tomography scans of the phantom (Catphan(®)) were performed with varying tube voltages (80-140kV) and currents (30-200mAs). 198 phantom data sets were reconstructed applying FBP and IR with increasing iterations, and soft and sharp kernels. Further, 25 chest and abdomen CT scans, performed with high and low exposure per patient, were reconstructed with IR and FBP. Two independent observers evaluated image quality and radiation doses of both phantom and patient scans. RESULTS In phantom scans, noise reduction was significantly improved using IR with increasing iterations, independent from tissue, scan-mode, tube-voltage, current, and kernel. IR did not affect high-contrast resolution. Low-contrast resolution was also not negatively affected, but improved in scans with doses <5mGy, although object detectability generally decreased with the lowering of exposure. At comparable image quality levels, CTDIvol was reduced by 26-50% using IR. In patients, applying IR vs. FBP resulted in good to excellent image quality, while tube voltage and current settings could be significantly decreased. CONCLUSIONS Our phantom experiments demonstrate that image quality levels of FBP reconstructions can also be achieved at lower tube voltages and tube currents when applying IR. Our findings could be confirmed in patients revealing the potential of IR to significantly reduce CT radiation doses.

Feasible Dose Reduction in Routine Chest Computed Tomography Maintaining Constant Image Quality Using the Last Three Scanner Generations: From Filtered Back Projection to Sinogram-affirmed Iterative Reconstruction and Impact of the Novel Fully Integrated Detector Design Minimizing Electronic Noise

OBJECTIVE The aim of the present study was to evaluate a dose reduction in contrast-enhanced chest computed tomography (CT) by comparing the three latest generations of Siemens CT scanners used in clinical practice. We analyzed the amount of radiation used with filtered back projection (FBP) and an iterative reconstruction (IR) algorithm to yield the same image quality. Furthermore, the influence on the radiation dose of the most recent integrated circuit detector (ICD; Stellar detector, Siemens Healthcare, Erlangen, Germany) was investigated. MATERIALS AND METHODS 136 Patients were included. Scan parameters were set to a thorax routine: SOMATOM Sensation 64 (FBP), SOMATOM Definition Flash (IR), and SOMATOM Definition Edge (ICD and IR). Tube current was set constantly to the reference level of 100 mA automated tube current modulation using reference milliamperes. Care kV was used on the Flash and Edge scanner, while tube potential was individually selected between 100 and 140 kVp by the medical technologists at the SOMATOM Sensation. Quality assessment was performed on soft-tissue kernel reconstruction. Dose was represented by the dose length product. RESULTS Dose-length product (DLP) with FBP for the average chest CT was 308 mGy*cm ± 99.6. In contrast, the DLP for the chest CT with IR algorithm was 196.8 mGy*cm ± 68.8 (P = 0.0001). Further decline in dose can be noted with IR and the ICD: DLP: 166.4 mGy*cm ± 54.5 (P = 0.033). The dose reduction compared to FBP was 36.1% with IR and 45.6% with IR/ICD. Signal-to-noise ratio (SNR) was favorable in the aorta, bone, and soft tissue for IR/ICD in combination compared to FBP (the P values ranged from 0.003 to 0.048). Overall contrast-to-noise ratio (CNR) improved with declining DLP. CONCLUSION The most recent technical developments, namely IR in combination with integrated circuit detectors, can significantly lower radiation dose in chest CT examinations.

Vascularised fibula grafts for early tibia reconstruction in infants with congenital pseudarthrosis

Congenital pseudarthrosis of the tibia (CPT) is caused by an ill-defined, segmental disturbance of periosteal bone formation leading to spontaneous bowing, followed by fracture and subsequent pseudarthrosis in the first 2 years of life. The results of conventional treatment modalities (e.g., bracing, internal and external fixation and bone grafting) are associated with high failure rates in terms of persisting pseudarthrosis, malunion and impaired growth. As a more promising alternative, a more aggressive approach, including wide resection of the affected bone, reconstruction with free vascularised fibula grafts from the healthy contralateral leg and stable external fixation at a very early stage has been suggested. Between 1995 and 2007, 10 children (age 12-31 months, median 20 months) suffering from CPT were treated at our institutions according to this principle. Two patients were treated before a fracture had occurred. The length of the fibula graft was 7-9cm. End-to-end anastomoses were performed at the level of the distal tibia stump. The follow-up was 80 months (median, range 12 months to 12 years). Radiologic examination at 6 weeks postoperatively showed normal bone density and structure of the transplanted fibula in all cases and osseous consolidation at 19 of the 20 graft/tibia junctions. One nonunion was sucessfully treated with bone grafting and plate osteosynthesis. Pin-tract infection occurred in three patients. Five children sustained graft fractures that were successfully treated with internal or external fixation. Two patients developed diminished growth of the affected limb or foot; all others had equal limb length and shoe size. At long-term follow-up, tibialisation of the transplant had occurred, and normal gait and physical activities were possible in all children. We conclude that in spite of a relatively high complication rate and the reluctance to perform free flap surgery in infants at this young age, the present concept may successfully prevent the imminent severe sequelae associated with CPT.

Aesthetic outcome and oncological safety of nipple-areola complex replantation after mastectomy and immediate breast reconstruction

Immediate breast reconstruction (IBR) has become an established procedure for women necessitating mastectomy. Traditionally, the nipple-areola complex (NAC) is resected during this procedure. The NAC, in turn, is a principal factor determining aesthetic outcome after breast reconstruction, and due to its particular texture and shape, a natural-looking NAC can barely be reconstructed with other tissues. The aim of this study was to assess the oncological safety as well as morbidity and aesthetic outcome after replantation of the NAC some days after IBR. Retrospective analysis of 85 patients receiving 88 mastectomies and IBR between 1998 and 2007 was conducted. NAC (n=29) or the nipple alone (n=23) were replanted 7 days (median, range 2-10 days) after IBR in 49 patients, provided the subareolar tissue was histologically negative for tumour infiltration. Local recurrence rate was assessed after 49 months (median, range 6-120 months). Aesthetic outcome was evaluated by clinical assessment during routine follow-up at least 12 months after the last intervention. Malignant involvement of the subareolar tissue was found in eight cases (9.1%). Patients qualifying for NAC replantation were in stage 0 in 29%, stage I in 15%, stage IIa in 31%, stage IIb in 17% and stage III in 8%. Total or partial necrosis occurred in 69% and 26% if the entire NAC or only the nipple were replanted, respectively (P<0.01). Depigmentation was seen in 52% and corrective surgery was done in 11 out of 52 NAC or nipple replantations. Local recurrence and isolated regional lymph node metastasis were observed in one single case each. Another 5.8% of the patients showed distant metastases. We conclude that the replantation of the NAC in IBR is oncologically safe, provided the subareolar tissue is free of tumour. However, the long-term aesthetic outcome of NAC replantation is not satisfying, which advocates replanting the nipple alone.

Computer-assisted design and manufacture of implants in the late reconstruction of extensive orbital fractures

To evaluate the use of computer-assisted designed and manufactured (CAD/CAM) orbital wall and floor implants for late reconstruction of extensive orbital fractures.

Single fluoroscopic image based 2D/3D reconstruction of a scaled, patient-specific vertebral model

Purpose Accurate three-dimensional (3D) models of lumbar vertebrae can enable image-based 3D kinematic analysis. The common approach to derive 3D models is by direct segmentation of CT or MRI datasets. However, these have the disadvantages that they are expensive, timeconsuming and/or induce high-radiation doses to the patient. In this study, we present a technique to automatically reconstruct a scaled 3D lumbar vertebral model from a single two-dimensional (2D) lateral fluoroscopic image. Methods Our technique is based on a hybrid 2D/3D deformable registration strategy combining a landmark-to-ray registration with a statistical shape model-based 2D/3D reconstruction scheme. Fig. 1 shows different stages of the reconstruction process. Four cadaveric lumbar spine segments (total twelve lumbar vertebrae) were used to validate the technique. To evaluate the reconstruction accuracy, the surface models reconstructed from the lateral fluoroscopic images were compared to the associated ground truth data derived from a 3D CT-scan reconstruction technique. For each case, a surface-based matching was first used to recover the scale and the rigid transformation between the reconstructed surface model Results Our technique could successfully reconstruct 3D surface models of all twelve vertebrae. After recovering the scale and the rigid transformation between the reconstructed surface models and the ground truth models, the average error of the 2D/3D surface model reconstruction over the twelve lumbar vertebrae was found to be 1.0 mm. The errors of reconstructing surface models of all twelve vertebrae are shown in Fig. 2. It was found that the mean errors of the reconstructed surface models in comparison to their associated ground truths after iterative scaled rigid registrations ranged from 0.7 mm to 1.3 mm and the rootmean squared (RMS) errors ranged from 1.0 mm to 1.7 mm. The average mean reconstruction error was found to be 1.0 mm. Conclusion An accurate, scaled 3D reconstruction of the lumbar vertebra can be obtained from a single lateral fluoroscopic image using a statistical shape model based 2D/3D reconstruction technique. Future work will focus on applying the reconstructed model for 3D kinematic analysis of lumbar vertebrae, an extension of our previously-reported imagebased kinematic analysis. The developed method also has potential applications in surgical planning and navigation.

Validation of statistical shape model based reconstruction of the proximal femur - A morphology study

Seventeen bones (sixteen cadaveric bones and one plastic bone) were used to validate a method for reconstructing a surface model of the proximal femur from 2D X-ray radiographs and a statistical shape model that was constructed from thirty training surface models. Unlike previously introduced validation studies, where surface-based distance errors were used to evaluate the reconstruction accuracy, here we propose to use errors measured based on clinically relevant morphometric parameters. For this purpose, a program was developed to robustly extract those morphometric parameters from the thirty training surface models (training population), from the seventeen surface models reconstructed from X-ray radiographs, and from the seventeen ground truth surface models obtained either by a CT-scan reconstruction method or by a laser-scan reconstruction method. A statistical analysis was then performed to classify the seventeen test bones into two categories: normal cases and outliers. This classification step depends on the measured parameters of the particular test bone. In case all parameters of a test bone were covered by the training population's parameter ranges, this bone is classified as normal bone, otherwise as outlier bone. Our experimental results showed that statistically there was no significant difference between the morphometric parameters extracted from the reconstructed surface models of the normal cases and those extracted from the reconstructed surface models of the outliers. Therefore, our statistical shape model based reconstruction technique can be used to reconstruct not only the surface model of a normal bone but also that of an outlier bone.

On-table reconstruction of comminuted fractures of the radial head

The most widely accepted treatment for comminuted fractures of the radial head is either the excision or open reduction and internal fixation. The purpose of the present study is to evaluate the value of an 'on-table' reconstruction technique in severely comminuted fractures of the radial head. In this study, two patients with a Mason type-III and four patients with a Mason type-IV radial-head fracture were treated with 'on-table' reconstruction and fixation using low-profile mini-plates. After a mean follow-up of 112 months (47-154 months), the mean elbow motion was 0-6-141 degrees extension flexion with 79 degrees of pronation and 70 degrees of supination. The mean Broberg and Morrey functional rating score was 97.0 points, the Mayo Elbow Performance Index was 99.2 points and the mean Disabilities of the Arm, Shoulder, and Hand (DASH) Outcome Measure score was 1.94 points. One patient had symptoms of degenerative changes, with a slight joint-space narrowing. There were no radiographic signs of devitalisation at final examination. Comminuted fractures of the radial head, which would otherwise require excision, can be successfully treated with an 'on-table' reconstruction technique.