Monte Carlo Dose Calculation on Deforming Anatomy

This article presents the implementation and validation of a dose calculation approach for deforming anatomical objects. Deformation is represented by deformation vector fields leading to deformed voxel grids representing the different deformation scenarios. Particle transport in the resulting deformed voxels is handled through the approximation of voxel surfaces by triangles in the geometry implementation of the Swiss Monte Carlo Plan framework. The focus lies on the validation methodology which uses computational phantoms representing the same physical object through regular and irregular voxel grids. These phantoms are chosen such that the new implementation for a deformed voxel grid can be compared directly with an established dose calculation algorithm for regular grids. Furthermore, separate validation of the aspects voxel geometry and the density changes resulting from deformation is achieved through suitable design of the validation phantom. We show that equivalent results are obtained with the proposed method and that no statistically significant errors are introduced through the implementation for irregular voxel geometries. This enables the use of the presented and validated implementation for further investigations of dose calculation on deforming anatomy.

Meta-Environment and Executable Meta-Language using Smalltalk: an Experience Report

Object-oriented modelling languages such as EMOF are often used to specify domain specific meta-models. However, these modelling languages lack the ability to describe behavior or operational semantics. Several approaches have used a subset of Java mixed with OCL as executable meta-languages. In this experience report we show how we use Smalltalk as an executable meta-language in the context of the Moose reengineering environment. We present how we implemented EMOF and its behavioral aspects. Over the last decade we validated this approach through incrementally building a meta-described reengineering environment. Such an approach bridges the gap between a code-oriented view and a meta-model driven one. It avoids the creation of yet another language and reuses the infrastructure and run-time of the underlying implementation language. It offers an uniform way of letting developers focus on their tasks while at the same time allowing them to meta-describe their domain model. The advantage of our approach is that developers use the same tools and environment they use for their regular tasks. Still the approach is not Smalltalk specific but can be applied to language offering an introspective API such as Ruby, Python, CLOS, Java and C#.

A Closed-Form, Consistent and Robust Solution to Uncalibrated Photometric Stereo Via Local Diffuse Reflectance Maxima

Images of an object under different illumination are known to provide strong cues about the object surface. A mathematical formalization of how to recover the normal map of such a surface leads to the so-called uncalibrated photometric stereo problem. In the simplest instance, this problem can be reduced to the task of identifying only three parameters: the so-called generalized bas-relief (GBR) ambiguity. The challenge is to find additional general assumptions about the object, that identify these parameters uniquely. Current approaches are not consistent, i.e., they provide different solutions when run multiple times on the same data. To address this limitation, we propose exploiting local diffuse reflectance (LDR) maxima, i.e., points in the scene where the normal vector is parallel to the illumination direction (see Fig. 1). We demonstrate several noteworthy properties of these maxima: a closed-form solution, computational efficiency and GBR consistency. An LDR maximum yields a simple closed-form solution corresponding to a semi-circle in the GBR parameters space (see Fig. 2); because as few as two diffuse maxima in different images identify a unique solution, the identification of the GBR parameters can be achieved very efficiently; finally, the algorithm is consistent as it always returns the same solution given the same data. Our algorithm is also remarkably robust: It can obtain an accurate estimate of the GBR parameters even with extremely high levels of outliers in the detected maxima (up to 80 % of the observations). The method is validated on real data and achieves state-of-the-art results.

Computational Comparison of Continuous and Discontinuous Galerkin Time-Stepping Methods for Nonlinear Initial Value Problems

This article centers on the computational performance of the continuous and discontinuous Galerkin time stepping schemes for general first-order initial value problems in R n , with continuous nonlinearities. We briefly review a recent existence result for discrete solutions from [6], and provide a numerical comparison of the two time discretization methods.

Computer Assisted Diagnosis and Treatment Planning of Femoroacetabular Impingement (FAI)

Femoroacetabular impingement (FAI) is a dynamic conflict of the hip defined by a pathological, early abutment of the proximal femur onto the acetabulum or pelvis. In the past two decades, FAI has received increasing focus in both research and clinical practice as a cause of hip pain and prearthrotic deformity. Anatomical abnormalities such as an aspherical femoral head (cam-type FAI), a focal or general overgrowth of the acetabulum (pincer-type FAI), a high riding greater or lesser trochanter (extra-articular FAI), or abnormal torsion of the femur have been identified as underlying pathomorphologies. Open and arthroscopic treatment options are available to correct the deformity and to allow impingement-free range of motion. In routine practice, diagnosis and treatment planning of FAI is based on clinical examination and conventional imaging modalities such as standard radiography, magnetic resonance arthrography (MRA), and computed tomography (CT). Modern software tools allow three-dimensional analysis of the hip joint by extracting pelvic landmarks from two-dimensional antero-posterior pelvic radiographs. An object-oriented cross-platform program (Hip2Norm) has been developed and validated to standardize pelvic rotation and tilt on conventional AP pelvis radiographs. It has been shown that Hip2Norm is an accurate, consistent, reliable and reproducible tool for the correction of selected hip parameters on conventional radiographs. In contrast to conventional imaging modalities, which provide only static visualization, novel computer assisted tools have been developed to allow the dynamic analysis of FAI pathomechanics. In this context, a validated, CT-based software package (HipMotion) has been introduced. HipMotion is based on polygonal three-dimensional models of the patient’s pelvis and femur. The software includes simulation methods for range of motion, collision detection and accurate mapping of impingement areas. A preoperative treatment plan can be created by performing a virtual resection of any mapped impingement zones both on the femoral head-neck junction, as well as the acetabular rim using the same three-dimensional models. The following book chapter provides a summarized description of current computer-assisted tools for the diagnosis and treatment planning of FAI highlighting the possibility for both static and dynamic evaluation, reliability and reproducibility, and its applicability to routine clinical use.

Reduced muscle mass and bone size in pediatric patients with inflammatory bowel disease

BACKGROUND: Decreased bone mineral density has been reported in children with inflammatory bowel disease (IBD). We used peripheral quantitative computed tomography (pQCT) to assess bone mineralization, geometry, and muscle cross-sectional area (CSA) in pediatric IBD. METHODS: In a cross-sectional study, pQCT of the forearm was applied in 143 IBD patients (mean age 13.9 +/- 3.5 years); 29% were newly diagnosed, 98 had Crohn's disease, and 45 had ulcerative colitis. Auxological data, cumulative glucocorticoid dose, disease activity indices, laboratory markers for inflammation, and bone metabolism were related to the results of pQCT. RESULTS: Patients were compromised in height (-0.82 +/- 1.1 SD), weight (-0.77 +/- 1.0 SD), muscle mass (-1.12 +/- 1.0 SD), and total bone cross-sectional area (-0.79 +/- 1.0 SD) compared to age- and sex-matched healthy controls (z-scores). In newly diagnosed patients, the ratio of bone mineral mass per muscle CSA was higher than in those with longer disease duration (1.00 versus 0.30, P = 0.007). Serum albumin level and disease activity correlated with muscle mass, accounting for 41.0% of variability in muscle mass (P < 0.01). The trabecular bone mineral density z-score was on average at the lower normal level (-0.40 +/- 1.3 SD, P < 0.05). CONCLUSIONS: Reduced bone geometry was explained only in part by reduced height. Bone disease in children with IBD seems to be secondary to muscle wasting, which is already present at diagnosis. With longer disease duration, bone adapts to the lower muscle CSA. Serum albumin concentration is a good marker for muscle wasting and abnormal bone development.

Reduced trabecular bone mineral density and cortical thickness accompanied by increased outer bone circumference in metacarpal bone of rheumatoid arthritis patients: a cross-sectional study

Introduction The objective of this study was to assess three-dimensional bone geometry and density at the epiphysis and shaft of the third meta-carpal bone of rheumatoid arthritis (RA) patients in comparison to healthy controls with the novel method of peripheral quantitative computed tomography (pQCT). Methods PQCT scans were performed in 50 female RA patients and 100 healthy female controls at the distal epiphyses and shafts of the third metacarpal bone, the radius and the tibia. Reproducibility was determined by coefficient of varia-tion. Bone densitometric and geometric parameters were compared between the two groups and correlated to disease characteristics. Results Reproducibility of different pQCT parameters was between 0.7% and 2.5%. RA patients had 12% to 19% lower trabecular bone mineral density (BMD) (P ≤ 0.001) at the distal epiphyses of radius, tibia and metacarpal bone. At the shafts of these bones RA patients had 7% to 16% thinner cortices (P ≤ 0.03). Total cross-sectional area (CSA) at the metacarpal bone shaft of pa-tients was larger (between 5% and 7%, P < 0.02), and relative cortical area was reduced by 13%. Erosiveness by Ratingen score correlated negatively with tra-becular and total BMD at the epiphyses and shaft cortical thickness of all measured bones (P < 0.04). Conclusions Reduced trabecular BMD and thinner cortices at peripheral bones, and a greater bone shaft diameter at the metacarpal bone suggest RA spe-cific bone alterations. The proposed pQCT protocol is reliable and allows measuring juxta-articular trabecular BMD and shaft geometry at the metacarpal bone.

The binding mode of side chain- and C3-modified epothilones to tubulin

The tubulin-binding mode of C3- and C15-modified analogues of epothilone A (Epo A) was determined by NMR spectroscopy and computational methods and compared with the existing structural models of tubulin-bound natural Epo A. Only minor differences were observed in the conformation of the macrocycle between Epo A and the C3-modified analogues investigated. In particular, 3-deoxy- (compound 2) and 3-deoxy-2,3-didehydro-Epo A (3) were found to adopt similar conformations in the tubulin-binding cleft as Epo A, thus indicating that the 3-OH group is not essential for epothilones to assume their bioactive conformation. None of the available models of the tubulin-epothilone complex is able to fully recapitulate the differences in tubulin-polymerizing activity and microtubule-binding affinity between C20-modified epothilones 6 (C20-propyl), 7 (C20-butyl), and 8 (C20-hydroxypropyl). Based on the results of transferred NOE experiments in the presence of tubulin, the isomeric C15 quinoline-based Epo B analogues 4 and 5 show very similar orientations of the side chain, irrespective of the position of the nitrogen atom in the quinoline ring. The quinoline side chain stacks on the imidazole moiety of beta-His227 with equal efficiency in both cases, thus suggesting that the aromatic side chain moiety in epothilones contributes to tubulin binding through strong van der Waals interactions with the protein rather than hydrogen bonding involving the heteroaromatic nitrogen atom. These conclusions are in line with existing tubulin polymerization and microtubule-binding data for 4, 5, and Epo B.

Combined bone and soft-tissue augmentation surgery in temporo-orbital contour reconstruction

Temporal hollowing due to temporal muscle atrophy after standard skull base surgery is common. Various techniques have been previously described to correct the disfiguring defect. Most often reconstruction is performed using freehand molded polymethylmethacrylate cement. This method and material are insufficient in terms of aesthetic results and implant characteristics. We herein propose reconstruction of such defects with a polyetheretherketone (PEEK)-based patient-specific implant (PSI) including soft-tissue augmentation to preserve normal facial topography. We describe a patient who presented with a large temporo-orbital hemangioma that had been repaired with polymethylmethacrylate 25 years earlier. Because of a toxic skin atrophy fistula, followed by infection and meningitis, this initial implant had to be removed. The large, disfiguring temporo-orbital defect was reconstructed with a PEEK-based PSI. The lateral orbital wall and the temporal muscle atrophy were augmented with computer-aided design and surface modeling techniques. The operative procedure to implant and adopt the reconstructed PEEK-based PSI was simple, and an excellent cosmetic outcome was achieved. The postoperative clinical course was uneventful over a 5-year follow-up period. Polyetheretherketone-based combined bony and soft contour remodeling is a feasible and effective method for cranioplasty including combined bone and soft-tissue reconstruction of temporo-orbital defects. Manual reconstruction of this cosmetically delicate area carries an exceptional risk of disfiguring results. Augmentation surgery in this anatomic location needs accurate PSIs to achieve satisfactory cosmetic results. The cosmetic outcome achieved in this case is superior compared with previously reported techniques.

Statistical finite element model for bone shape and biomechanical properties

We present a framework for statistical finite element analysis combining shape and material properties, and allowing performing statistical statements of biomechanical performance across a given population. In this paper, we focus on the design of orthopaedic implants that fit a maximum percentage of the target population, both in terms of geometry and biomechanical stability. CT scans of the bone under consideration are registered non-rigidly to obtain correspondences in position and intensity between them. A statistical model of shape and intensity (bone density) is computed by means of principal component analysis. Afterwards, finite element analysis (FEA) is performed to analyse the biomechanical performance of the bones. Realistic forces are applied on the bones and the resulting displacement and bone stress distribution are calculated. The mechanical behaviour of different PCA bone instances is compared.

Ragu: a free tool for the analysis of EEG and MEG event-related scalp field data using global randomization statistics

We present a program (Ragu; Randomization Graphical User interface) for statistical analyses of multichannel event-related EEG and MEG experiments. Based on measures of scalp field differences including all sensors, and using powerful, assumption-free randomization statistics, the program yields robust, physiologically meaningful conclusions based on the entire, untransformed, and unbiased set of measurements. Ragu accommodates up to two within-subject factors and one between-subject factor with multiple levels each. Significance is computed as function of time and can be controlled for type II errors with overall analyses. Results are displayed in an intuitive visual interface that allows further exploration of the findings. A sample analysis of an ERP experiment illustrates the different possibilities offered by Ragu. The aim of Ragu is to maximize statistical power while minimizing the need for a-priori choices of models and parameters (like inverse models or sensors of interest) that interact with and bias statistics.

Effect of lateral meniscectomy and osteochondral grafting of a lateral femoral condylar defect on contact mechanics: a cadaveric study in dogs

BACKGROUND Osteochondral autograft transfer (OAT) aims at restoring normal articular cartilage surface geometry and articular contact mechanics. To date, no studies have evaluated the contact mechanics of the canine stifle following OAT. Additionally, there are no studies that evaluated the role of the meniscus in contact mechanics following OAT in human or canine femorotibial joints. The objective of this study was to measure the changes in femorotibial contact areas (CA), mean contact pressure (MCP) and peak contact pressure (PCP) before and after osteochondral autograft transplantation (OAT) of a simulated lateral femoral condylar cartilage defect with an intact lateral meniscus and following lateral meniscectomy. RESULTS With an intact lateral meniscus, creation of an osteochondral defect caused a decrease in MCP and PCP by 11% and 30%, respectively, compared to the intact stifle (p < 0.01). With an intact meniscus, implanting an osteochondral graft restored MCP and PCP to 96% (p = 0.56) and 92% (p = 0.41) of the control values. Lateral meniscectomy with grafting decreased CA by 54% and increased PCP by 79% compared to the intact stifle (p < 0.01). CONCLUSIONS OAT restored contact pressures in stifles with a simulated lateral condylar defect when the meniscus was intact. The lateral meniscus has a significant role in maintaining normal contact pressures in both stifles with a defect or following OAT. Meniscectomy should be avoided when a femoral condylar defect is present and when performing OAT.

Two-dimensional supramolecular polymers: controlled formation and visualization on the surface

Self – assembly is a powerful tool for the construction of highly organized nanostructures. Therefore, the possibility to control and predict pathways of molecular ordering on the nanoscale level is a critical issue for the production of materials with tunable and adaptive macroscopic properties. 2D polymers are attractive objects for the field of material sciences due to their exceptional properties. [1] As shown before, amphiphilic oligopyrenotides (produced via automated solid-phase synthesis) form rod–like supramolecular polymers in water. [2] These assemblies form 1D objects. [3] By applying certain changes to the design of the oligopyrenotide units the dimensionality of the formed assemblies can be influenced. Herein, we demonstrate that Py3 (see Figure 1) forms defined supramolecular assemblies under thermodynamic conditions in water. To study Py3 self-assembly, we carried out whole set of spectroscopic (UV/vis, fluorescence, DLS) and microscopic experiments (AFM). The obtained results suggest that oligopyrenotides with the present type of geometry and linker length leads to formation of 2D supramolecular assemblies.

Cortical remodeling during menopause, rheumatoid arthritis, glucocorticoid and bisphosphonate therapy

Bone mass, bone geometry and its changes are based on trabecular and cortical bone remodeling. Whereas the effects of estrogen loss, rheumatoid arthritis (RA), glucocorticoid (GC) and bisphosphonate (BP) on trabecular bone remodeling have been well described, the effects of these conditions on the cortical bone geometry are less known. The present review will report current knowledge on the effects of RA, GC and BP on cortical bone geometry and its clinical relevance. Estrogen deficiency, RA and systemic GC lead to enhanced endosteal bone resorption. While in estrogen deficiency and under GC therapy endosteal resorption is insufficiently compensated by periosteal apposition, RA is associated with some periosteal bone apposition resulting in a maintained load-bearing capacity and stiffness. In contrast, BP treatment leads to filling of endosteal bone cavities at the epiphysis; however, periosteal apposition at the bone shaft seems to be suppressed. In summary, estrogen loss, RA and GC show similar effects on endosteal bone remodeling with an increase in bone resorption, whereas their effect on periosteal bone remodeling may differ. Despite over 50 years of GC therapy and over 25 years of PB therapy, there is still need for better understanding of the skeletal effects of these drugs as well as of inflammatory disease such as RA on cortical bone remodeling.