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

[Osteoporosis - whom to treat? The importance of FRAX® in Switzerland]

The WHO fracture risk assessment tool FRAX® is a computer based algorithm that provides models for the assessment of fracture probability in men and women. The approach uses easily obtained clinical risk factors (CRFs) to estimate 10-year probability of a major osteoporotic fracture (hip, clinical spine, humerus or wrist fracture) and the 10-year probability of a hip fracture. The estimate can be used alone or with femoral neck bone mineral density (BMD) to enhance fracture risk prediction. FRAX® is the only risk engine which takes into account the hazard of death as well as that of fracture. Probability of fracture is calculated in men and women from age, body mass index, and dichotomized variables that comprise a prior fragility fracture, parental history of hip fracture, current tobacco smoking, ever long-term use of oral glucocorticoids, rheumatoid arthritis, other causes of secondary osteoporosis, daily alcohol consumption of 3 or more units daily. The relationship between risk factors and fracture probability was constructed using information of nine population-based cohorts from around the world. CRFs for fracture had been identified that provided independent information on fracture risk based on a series of meta-analyses. The FRAX® algorithm was validated in 11 independent cohorts with in excess of 1 million patient-years, including the Swiss SEMOF cohort. Since fracture risk varies markedly in different regions of the world, FRAX® models need to be calibrated to those countries where the epidemiology of fracture and death is known. Models are currently available for 31 countries across the world. The Swiss-specific FRAX® model was developed very soon after the first release of FRAX® in 2008 and was published in 2009, using Swiss epidemiological data, integrating fracture risk and death hazard of our country. Two FRAX®-based approaches may be used to explore intervention thresholds. They have recently been investigated in the Swiss setting. In the first approach the guideline that individuals with a fracture probability equal to or exceeding that of women with a prior fragility fracture should be considered for treatment is translated into thresholds using 10-year fracture probabilities. In that case the threshold is age-dependent and increases from 16 % at the age of 60 ys to 40 % at the age of 80 ys. The second approach is a cost-effectiveness approach. Using a FRAX®-based intervention threshold of 15 % for both, women and men 50 years and older, should permit cost-effective access to therapy to patients at high fracture probability in our country and thereby contribute to further reduce the growing burden of osteoporotic fractures.

Articulated Statistical Shape Model-Based 2D-3D Reconstruction of a Hip Joint

In this paper, reconstruction of three-dimensional (3D) patient-specific models of a hip joint from two-dimensional (2D) calibrated X-ray images is addressed. Existing 2D-3D reconstruction techniques usually reconstruct a patient-specific model of a single anatomical structure without considering the relationship to its neighboring structures. Thus, when those techniques would be applied to reconstruction of patient-specific models of a hip joint, the reconstructed models may penetrate each other due to narrowness of the hip joint space and hence do not represent a true hip joint of the patient. To address this problem we propose a novel 2D-3D reconstruction framework using an articulated statistical shape model (aSSM). Different from previous work on constructing an aSSM, where the joint posture is modeled as articulation in a training set via statistical analysis, here it is modeled as a parametrized rotation of the femur around the joint center. The exact rotation of the hip joint as well as the patient-specific models of the joint structures, i.e., the proximal femur and the pelvis, are then estimated by optimally fitting the aSSM to a limited number of calibrated X-ray images. Taking models segmented from CT data as the ground truth, we conducted validation experiments on both plastic and cadaveric bones. Qualitatively, the experimental results demonstrated that the proposed 2D-3D reconstruction framework preserved the hip joint structure and no model penetration was found. Quantitatively, average reconstruction errors of 1.9 mm and 1.1 mm were found for the pelvis and the proximal femur, respectively.

FRAX® assessment of osteoporotic fracture probability in Switzerland

A Swiss-specific FRAX model was developed. Patient profiles at increased probability of fracture beyond currently accepted reimbursement thresholds for bone mineral density (BMD) measurement by dual X-ray absorptiometry (DXA), and osteoporosis treatment were identified.

Influence of Smoothing on Voxel-Based Mesh Accuracy in Micro-Finite Element

The interest in automatic volume meshing for finite element analysis (FEA) has grown more since the appearance of microfocus CT (μCT), due to its high resolution, which allows for the assessment of mechanical behaviour at a high precision. Nevertheless, the basic meshing approach of generating one hexahedron per voxel produces jagged edges. To prevent this effect, smoothing algorithms have been introduced to enhance the topology of the mesh. However, whether smoothing also improves the accuracy of voxel-based meshes in clinical applications is still under question. There is a trade-off between smoothing and quality of elements in the mesh. Distorted elements may be produced by excessive smoothing and reduce accuracy of the mesh. In the present work, influence of smoothing on the accuracy of voxel-based meshes in micro-FE was assessed. An accurate 3D model of a trabecular structure with known apparent mechanical properties was used as a reference model. Virtual CT scans of this reference model (with resolutions of 16, 32 and 64 μm) were then created and used to build voxel-based meshes of the microarchitecture. Effects of smoothing on the apparent mechanical properties of the voxel-based meshes as compared to the reference model were evaluated. Apparent Young’s moduli of the smooth voxel-based mesh were significantly closer to those of the reference model for the 16 and 32 μm resolutions. Improvements were not significant for the 64 μm, due to loss of trabecular connectivity in the model. This study shows that smoothing offers a real benefit to voxel-based meshes used in micro-FE. It might also broaden voxel-based meshing to other biomechanical domains where it was not used previously due to lack of accuracy. As an example, this work will be used in the framework of the European project ContraCancrum, which aims at providing a patient-specific simulation of tumour development in brain and lungs for oncologists. For this type of clinical application, such a fast, automatic, and accurate generation of the mesh is of great benefit.

Multi-Scale Modeling for Image Analysis of Brain Tumor Studies

Image-based modeling of tumor growth combines methods from cancer simulation and medical imaging. In this context, we present a novel approach to adapt a healthy brain atlas to MR images of tumor patients. In order to establish correspondence between a healthy atlas and a pathologic patient image, tumor growth modeling in combination with registration algorithms is employed. In a first step, the tumor is grown in the atlas based on a new multi-scale, multi-physics model including growth simulation from the cellular level up to the biomechanical level, accounting for cell proliferation and tissue deformations. Large-scale deformations are handled with an Eulerian approach for finite element computations, which can operate directly on the image voxel mesh. Subsequently, dense correspondence between the modified atlas and patient image is established using nonrigid registration. The method offers opportunities in atlasbased segmentation of tumor-bearing brain images as well as for improved patient-specific simulation and prognosis of tumor progression.

A Portable Image Overlay Projection Device for Computer-Aided Open Liver Surgery

Image overlay projection is a form of augmented reality that allows surgeons to view underlying anatomical structures directly on the patient surface. It improves intuitiveness of computer-aided surgery by removing the need for sight diversion between the patient and a display screen and has been reported to assist in 3-D understanding of anatomical structures and the identification of target and critical structures. Challenges in the development of image overlay technologies for surgery remain in the projection setup. Calibration, patient registration, view direction, and projection obstruction remain unsolved limitations to image overlay techniques. In this paper, we propose a novel, portable, and handheld-navigated image overlay device based on miniature laser projection technology that allows images of 3-D patient-specific models to be projected directly onto the organ surface intraoperatively without the need for intrusive hardware around the surgical site. The device can be integrated into a navigation system, thereby exploiting existing patient registration and model generation solutions. The position of the device is tracked by the navigation system’s position sensor and used to project geometrically correct images from any position within the workspace of the navigation system. The projector was calibrated using modified camera calibration techniques and images for projection are rendered using a virtual camera defined by the projectors extrinsic parameters. Verification of the device’s projection accuracy concluded a mean projection error of 1.3 mm. Visibility testing of the projection performed on pig liver tissue found the device suitable for the display of anatomical structures on the organ surface. The feasibility of use within the surgical workflow was assessed during open liver surgery. We show that the device could be quickly and unobtrusively deployed within the sterile environment.

Effectiveness of the cognitive differentiation program of the integrated psychological therapy: group versus individual treatment

The aim of the current pilot study was to compare two strategies in the application of the cognitive differentiation program of Integrated Psychological Therapy for people with schizophrenia. Twenty-six outpatients were randomly assigned to the application of the program in group sessions (CDg), or to its application in individualized sessions (CDi). The program provides cognitive exercises to promote better performance in cognition, and both groups of participants completed the same number of exercises following the same number of sessions per week. Outcomes were assessed on neuropsychological measures of attention, executive functioning and everyday memory, and everyday functioning. Effect sizes showed the absence of effects in everyday memory and social functioning, higher improvements in the CDi group in attention, and a higher improvement in the CDg condition in executive functioning. The results suggest that the program application model could be individualized, depending on patient-specific cognitive deficits.

Patient-related risk factors leading to aseptic stem loosening in total hip arthroplasty: a case-control study of 5,035 patients

BACKGROUND: We hypothesized that certain patient characteristics have different effects on the risk of early stem loosening in total hip arthroplasty (THA). We therefore conducted a case-control study using register-database records with the aim of identifying patient-specific risk factors associated with radiographic signs of aseptic loosening of the femoral component in THA. METHOD: Data were derived from a multinational European registry and were collected over a period of 25 years. 725 cases with radiographic signs of stem loosening were identified and matched to 4,310 controls without any signs of loosening. Matching criteria were type of implant, size of head, date of operation, center of primary intervention, and follow-up time. The risk factors analyzed were age at operation, sex, diagnosis and previous ipsilateral operations, height, weight, body mass index and mobility based on the Charnley classification. RESULTS: Women showed significantly lower risk of radiographic loosening than men (odds ratio (OR) 0.64). Age was also a strong factor: risk decreased by 1.8% for each additional year of age at the time of surgery. Height and weight were not associated with risk of loosening. A higher body mass index, however, increased the risk of stem loosening to a significant extent (OR 1.03) per additional unit of BMI. Charnley Class B, indicating restricted mobility, was associated with lower risk of loosening (OR 0.78). INTERPRETATION: An increased activity level, as seen in younger patients and those with unrestricted mobility, is an important factor in the etiology of stem loosening. If combined with high BMI, the risk of stem loosening within 10 years is even higher. A younger person should not be denied the benefits of a total hip arthroplasty but must accept that the risk of future failure is increased.

Remaining lifetime and absolute 10-year probabilities of osteoporotic fracture in Swiss men and women

SUMMARY: Remaining lifetime and absolute 10-year probabilities for osteoporotic fractures were determined by gender, age, and BMD values. Remaining lifetime probability at age 50 years was 20.2% in men and 51.3% in women and increased with advancing age and decreasing BMD. The study validates the elements required to populate a Swiss-specific FRAX model. INTRODUCTION: Switzerland belongs to high-risk countries for osteoporosis. Based on demographic projections, burden will still increase. We assessed remaining lifetime and absolute 10-year probabilities for osteoporotic fractures by gender, age and BMD in order to populate FRAX algorithm for Switzerland. METHODS: Osteoporotic fracture incidence was determined from national epidemiological data for hospitalised fractured patients from the Swiss Federal Office of Statistics in 2000 and results of a prospective Swiss cohort with almost 5,000 fractured patients in 2006. Validated BMD-associated fracture risk was used together with national death incidence and risk tables to determine remaining lifetime and absolute 10-year fracture probabilities for hip and major osteoporotic (hip, spine, distal radius, proximal humerus) fractures. RESULTS: Major osteoporotic fractures incidence was 773 and 2,078 per 100,000 men and women aged 50 and older. Corresponding remaining lifetime probabilities at age 50 were 20.2% and 51.3%. Hospitalisation for clinical spine, distal radius, and proximal humerus fractures reached 25%, 30% and 50%, respectively. Absolute 10-year probability of osteoporotic fracture increased with advancing age and decreasing BMD and was higher in women than in men. CONCLUSION: This study validates the elements required to populate a Swiss-specific FRAX model, a country at highest risk for osteoporotic fractures.

Integrated Segmentation of Brain Tumor Images for Radiotherapy and Neurosurgery

Image-based modeling of tumor growth combines methods from cancer simulation and medical imaging. In this context, we present a novel approach to adapt a healthy brain atlas to MR images of tumor patients. In order to establish correspondence between a healthy atlas and a pathologic patient image, tumor growth modeling in combination with registration algorithms is employed. In a first step, the tumor is grown in the atlas based on a new multi-scale, multi-physics model including growth simulation from the cellular level up to the biomechanical level, accounting for cell proliferation and tissue deformations. Large-scale deformations are handled with an Eulerian approach for finite element computations, which can operate directly on the image voxel mesh. Subsequently, dense correspondence between the modified atlas and patient image is established using nonrigid registration. The method offers opportunities in atlasbased segmentation of tumor-bearing brain images as well as for improved patient-specific simulation and prognosis of tumor progression.

3D volumetric intensity reconsturction from 2D x-ray images using partial least squares regression

Reconstruction of shape and intensity from 2D x-ray images has drawn more and more attentions. Previously introduced work suffers from the long computing time due to its iterative optimization characteristics and the requirement of generating digitally reconstructed radiographs within each iteration. In this paper, we propose a novel method which uses a patient-specific 3D surface model reconstructed from 2D x-ray images as a surrogate to get a patient-specific volumetric intensity reconstruction via partial least squares regression. No DRR generation is needed. The method was validated on 20 cadaveric proximal femurs by performing a leave-one-out study. Qualitative and quantitative results demonstrated the efficacy of the present method. Compared to the existing work, the present method has the advantage of much shorter computing time and can be applied to both DXA images as well as conventional x-ray images, which may hold the potentials to be applied to clinical routine task such as total hip arthroplasty (THA).

Finite element analyses of human vertebral bodies embedded in polymethylmethalcrylate or loaded via the hyperelastic intervertebral disc models provide equivalent predictions of experimental strength

Quantitative computer tomography (QCT)-based finite element (FE) models of vertebral body provide better prediction of vertebral strength than dual energy X-ray absorptiometry. However, most models were validated against compression of vertebral bodies with endplates embedded in polymethylmethalcrylate (PMMA). Yet, loading being as important as bone density, the absence of intervertebral disc (IVD) affects the strength. Accordingly, the aim was to assess the strength predictions of the classic FE models (vertebral body embedded) against the in vitro and in silico strengths of vertebral bodies loaded via IVDs. High resolution peripheral QCT (HR-pQCT) were performed on 13 segments (T11/T12/L1). T11 and L1 were augmented with PMMA and the samples were tested under a 4° wedge compression until failure of T12. Specimen-specific model was generated for each T12 from the HR-pQCT data. Two FE sets were created: FE-PMMA refers to the classical vertebral body embedded model under axial compression; FE-IVD to their loading via hyperelastic IVD model under the wedge compression as conducted experimentally. Results showed that FE-PMMA models overestimated the experimental strength and their strength prediction was satisfactory considering the different experimental set-up. On the other hand, the FE-IVD models did not prove significantly better (Exp/FE-PMMA: R²=0.68; Exp/FE-IVD: R²=0.71, p=0.84). In conclusion, FE-PMMA correlates well with in vitro strength of human vertebral bodies loaded via real IVDs and FE-IVD with hyperelastic IVDs do not significantly improve this correlation. Therefore, it seems not worth adding the IVDs to vertebral body models until fully validated patient-specific IVD models become available.

Extraspinal ossifications after implantation of vertical expandable prosthetic titanium ribs (VEPTRs).

INTRODUCTION Though developed for thoracic insufficiency syndrome, the spinal growth-stimulating potential and the ease of placement of vertical expandable titanium ribs (VEPTRs) has resulted in their widespread use for early-onset spine deformity. Observation of implant-related ossifications warrants further assessment, since they may be detrimental to the function-preserving non-fusion strategy. PATIENTS AND METHODS Radiographs (obtained pre and post index procedure, and at 4-year follow-up) and the records of 65 VEPTR patients from four paediatric spine centres were analysed. Ossifications were classified as type I (at anchor points), type II (along the central part) or type III (re-ossification after thoracostomy). RESULTS The average age at the index procedure was 6.5 years (min 1, max 13.7). The most prevalent spine problem was congenital scoliosis (37) with rib fusions (34), followed by neuromuscular and syndromic deformities (13 and 8, respectively). Idiopathic and secondary scoliosis (e.g. after thoracotomy) were less frequent (3 and 4, respectively). Forty-two of the 65 (65 %) patients showed ossifications, half of which were around the anchors. Forty-five percent (15/33) without pre-existing rib fusions developed a type II ossification along the implant. Re-ossifications of thoracostomies were less frequent (5/34, 15 %). The occurrence of ossifications was not associated with patient-specific factors. CONCLUSIONS Implant-related ossifications around VEPTR are common. In contrast to harmless bone formation around anchors, ossifications around the telescopic part and the rod section are troublesome in view of their possible negative impact on chest cage compliance and spinal mobility. This potential side effect needs to be considered during implant selection, particularly in patients with originally normal thoracic and spinal anatomy.

Comparison of morphological and rheological conditions between conventional and eversion carotid endarterectomy using computational fluid dynamics - a pilot study

PURPOSE To compare postoperative morphological and rheological conditions after eversion carotid endarterectomy versus conventional carotid endarterectomy using computational fluid dynamics. BASIC METHODS Hemodynamic metrics (velocity, wall shear stress, time-averaged wall shear stress and temporal gradient wall shear stress) in the carotid arteries were simulated in one patient after conventional carotid endarterectomy and one patient after eversion carotid endarterectomy by computational fluid dynamics analysis based on patient specific data. PRINCIPAL FINDINGS Systolic peak of the eversion carotid endarterectomy model showed a gradually decreased pressure along the stream path, the conventional carotid endarterectomy model revealed high pressure (about 180 Pa) at the carotid bulb. Regions of low wall shear stress in the conventional carotid endarterectomy model were much larger than that in the eversion carotid endarterectomy model and with lower time-averaged wall shear stress values (conventional carotid endarterectomy: 0.03-5.46 Pa vs. eversion carotid endarterectomy: 0.12-5.22 Pa). CONCLUSIONS Computational fluid dynamics after conventional carotid endarterectomy and eversion carotid endarterectomy disclosed differences in hemodynamic patterns. Larger studies are necessary to assess whether these differences are consistent and might explain different rates of restenosis in both techniques.