Automatic extraction of the mid-facial plane for cranio-maxillofacial surgery planning

Recently developed computer applications provide tools for planning cranio-maxillofacial interventions based on 3-dimensional (3D) virtual models of the patient's skull obtained from computed-tomography (CT) scans. Precise knowledge of the location of the mid-facial plane is important for the assessment of deformities and for planning reconstructive procedures. In this work, a new method is presented to automatically compute the mid-facial plane on the basis of a surface model of the facial skeleton obtained from CT. The method matches homologous surface areas selected by the user on the left and right facial side using an iterative closest point optimization. The symmetry plane which best approximates this matching transformation is then computed. This new automatic method was evaluated in an experimental study. The study included experienced and inexperienced clinicians defining the symmetry plane by a selection of landmarks. This manual definition was systematically compared with the definition resulting from the new automatic method: Quality of the symmetry planes was evaluated by their ability to match homologous areas of the face. Results show that the new automatic method is reliable and leads to significantly higher accuracy than the manual method when performed by inexperienced clinicians. In addition, the method performs equally well in difficult trauma situations, where key landmarks are unreliable or absent.

Sonographic visualization and ultrasound-guided block of the third occipital nerve: prospective for a new method to diagnose C2-C3 zygapophysial joint pain

BACKGROUND: Chronic neck pain after whiplash injury is caused by cervical zygapophysial joints in 50% of patients. Diagnostic blocks of nerves supplying the joints are performed using fluoroscopy. The authors' hypothesis was that the third occipital nerve can be visualized and blocked with use of an ultrasound-guided technique. METHODS: In 14 volunteers, the authors placed a needle ultrasound-guided to the third occipital nerve on both sides of the neck. They punctured caudal and perpendicular to the 14-MHz transducer. In 11 volunteers, 0.9 ml of either local anesthetic or normal saline was applied in a randomized, double-blind, crossover manner. Anesthesia was controlled in the corresponding skin area by pinprick and cold testing. The position of the needle was controlled by fluoroscopy. RESULTS: The third occipital nerve could be visualized in all subjects and showed a median diameter of 2.0 mm. Anesthesia was missing after local anesthetic in only one case. There was neither anesthesia nor hyposensitivity after any of the saline injections. The C2-C3 joint, in a transversal plane visualized as a convex density, was identified correctly by ultrasound in 27 of 28 cases, and 23 needles were placed correctly into the target zone. CONCLUSIONS: The third occipital nerve can be visualized and blocked with use of an ultrasound-guided technique. The needles were positioned accurately in 82% of cases as confirmed by fluoroscopy; the nerve was blocked in 90% of cases. Because ultrasound is the only available technique today to visualize this nerve, it seems to be a promising new method for block guidance instead of fluoroscopy.

RF-sputter fabrication of magnetic garnet thin films and simulation modeling for 1-D magnetic photonic crystal waveguide devices

One dimensional magnetic photonic crystals (1D-MPC) are promising structures for integrated optical isolator applications. Rare earth substituted garnet thin films with proper Faraday rotation are required to fabricate planar 1D-MPCs. In this thesis, flat-top response 1D-MPC was proposed and spectral responses and Faraday rotation were modeled. Bismuth substituted iron garnet films were fabricated by RF magnetron sputtering and structures, compositions, birefringence and magnetooptical properties were studied. Double layer structures for single mode propagation were also fabricated by sputtering for the first time. Multilayer stacks with multiple defects (phase shift) composed of Ce-YIG and GGG quarter-wave plates were simulated by the transfer matrix method. The transmission and Faraday rotation characteristics were theoretically studied. It is found that flat-top response, with 100% transmission and near 45o rotation is achievable by adjusting the inter-defect spacing, for film structures as thin as 30 to 35 μm. This is better than 3-fold reduction in length compared to the best Ce-YIG films for comparable rotations, thus allows a considerable reduction in size in manufactured optical isolators. Transmission bands as wide as 7nm were predicted, which is considerable improvement over 2 defects structure. Effect of repetition number and ratio factor on transmission and Faraday rotation ripple factors for the case of 3 and 4 defects structure has been discussed. Diffraction across the structure corresponds to a longer optical path length. Thus the use of guided optics is required to minimize the insertion losses in integrated devices. This part is discussed in chapter 2 in this thesis. Bismuth substituted iron garnet thin films were prepared by RF magnetron sputtering. We investigated or measured the deposition parameters optimization, crystallinity, surface morphologies, composition, magnetic and magnetooptical properties. A very high crystalline quality garnet film with smooth surface has been heteroepitaxially grown on (111) GGG substrate for films less than 1μm. Dual layer structures with two distinct XRD peaks (within a single sputtered film) start to develop when films exceed this thickness. The development of dual layer structure was explained by compositional gradient across film thickness, rather than strain gradient proposed by other authors. Lower DC self bias or higher substrate temperature is found to help to delay the appearance of the 2nd layer. The deposited films show in-plane magnetization, which is advantageous for waveguide devices application. Propagation losses of fabricated waveguides can be decreased by annealing in an oxygen atmosphere from 25dB/cm to 10dB/cm. The Faraday rotation at λ=1.55μm were also measured for the waveguides. FR is small (10° for a 3mm long waveguide), due to the presence of linear birefringence. This part is covered in chapter 4. We also investigated the elimination of linear birefringence by thickness tuning method for our sputtered films. We examined the compressively and tensilely strained films and analyze the photoelastic response of the sputter deposited garnet films. It has been found that the net birefringence can be eliminated under planar compressive strain conditions by sputtering. Bi-layer GGG on garnet thin film yields a reduced birefringence. Temperature control during the sputter deposition of GGG cover layer is critical and strongly influences the magnetization and birefringence level in the waveguide. High temperature deposition lowers the magnetization and increases the linear birefringence in the garnet films. Double layer single mode structures fabricated by sputtering were also studied. The double layer, which shows an in-plane magnetization, has an increased RMS roughness upon upper layer deposition. The single mode characteristic was confirmed by prism coupler measurement. This part is discussed in chapter 5.

Muffler characterization with implementation of the finite element method and experimental techniques

Determining how an exhaust system will perform acoustically before a prototype muffler is built can save the designer both a substantial amount of time and resources. In order to effectively use the simulation tools available it is important to understand what is the most effective tool for the intended purpose of analysis as well as how typical elements in an exhaust system affect muffler performance. An in-depth look at the available tools and their most beneficial uses are presented in this thesis. A full parametric study was conducted using the FEM method for typical muffler elements which was also correlated to experimental results. This thesis lays out the overall ground work on how to accurately predict sound pressure levels in the free field for an exhaust system with the engine properties included. The accuracy of the model is heavily dependent on the correct temperature profile of the model in addition to the accuracy of the source properties. These factors will be discussed in detail and methods for determining them will be presented. The secondary effects of mean flow, which affects both the acoustical wave propagation and the flow noise generation, will be discussed. Effective ways for predicting these secondary effects will be described. Experimental models will be tested on a flow rig that showcases these phenomena.

Non-invasive method to predict viscosity and size using total internal reflection fluorescence microscopy (TIRFM) system

This study develops an automated analysis tool by combining total internal reflection fluorescence microscopy (TIRFM), an evanescent wave microscopic imaging technique to capture time-sequential images and the corresponding image processing Matlab code to identify movements of single individual particles. The developed code will enable us to examine two dimensional hindered tangential Brownian motion of nanoparticles with a sub-pixel resolution (nanoscale). The measured mean square displacements of nanoparticles are compared with theoretical predictions to estimate particle diameters and fluid viscosity using a nonlinear regression technique. These estimated values will be confirmed by the diameters and viscosities given by manufacturers to validate this analysis tool. Nano-particles used in these experiments are yellow-green polystyrene fluorescent nanospheres (200 nm, 500 nm and 1000 nm in diameter (nominal); 505 nm excitation and 515 nm emission wavelengths). Solutions used in this experiment are de-ionized (DI) water, 10% d-glucose and 10% glycerol. Mean square displacements obtained near the surface shows significant deviation from theoretical predictions which are attributed to DLVO forces in the region but it conforms to theoretical predictions after ~125 nm onwards. The proposed automation analysis tool will be powerfully employed in the bio-application fields needed for examination of single protein (DNA and/or vesicle) tracking, drug delivery, and cyto-toxicity unlike the traditional measurement techniques that require fixing the cells. Furthermore, this tool can be also usefully applied for the microfluidic areas of non-invasive thermometry, particle tracking velocimetry (PTV), and non-invasive viscometry.

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.

Null-wave giant gravitons from thermal spinning brane probes

We construct and analyze thermal spinning giant gravitons in type II/M-theory based on spherically wrapped black branes, using the method of thermal probe branes originating from the blackfold approach. These solutions generalize in different directions recent work in which the case of thermal (non-spinning) D3-brane giant gravitons was considered, and reveal a rich phase structure with various new properties. First of all, we extend the construction to M-theory, by constructing thermal giant graviton solutions using spherically wrapped M2- and M5-branes. More importantly, we switch on new quantum numbers, namely internal spins on the sphere, which are not present in the usual extremal limit for which the brane world volume stress tensor is Lorentz invariant. We examine the effect of this new type of excitation and in particular analyze the physical quantities in various regimes, including that of small temperatures as well as low/high spin. As a byproduct we find new stationary dipole-charged black hole solutions in AdS m × S n backgrounds of type II/M-theory. We finally show, via a double scaling extremal limit, that our spinning thermal giant graviton solutions lead to a novel null-wave zero-temperature giant graviton solution with a BPS spectrum, which does not have an analogue in terms of the conventional weakly coupled world volume theory.

Effects of Thoratec pulsatile ventricular assist device (PVAD) timing on the abdominal aortic wave intensity pattern

Arterial waves are seen as possible independent mediators of cardiovascular risks, and the wave intensity analysis (WIA) has therefore been proposed as a method for patient selection for ventricular assist device (VAD) implantation. Interpreting measured wave intensity (WI) is challenging and complexity is increased by the implantation of a VAD. The waves generated by the VAD interact with the waves generated by the native heart, and this interaction varies with changing VAD settings. Eight sheep were implanted with a pulsatile VAD (PVAD) through ventriculo-aortic cannulation. The start of PVAD ejection was synchronized to the native R-wave and delayed between 0 % - 90 % of the cardiac cycle in 10 % steps or phase shifts (PS). Pressure and velocity signals were registered, using a combined Doppler and pressure wire positioned in the abdominal aorta, and used to calculate the WI. Depending on the PS, different wave interference phenomena occurred. Maximum unloading of the left ventricle (LV) coincided with constructive interference and maximum blood flow pulsatility, and maximum loading of the LV coincided with destructive interference and minimum blood flow pulsatility. We believe, that non-invasive WIA could potentially be used clinically to assess the mechanical load of the LV, and to monitor the peripheral hemodynamics such as blood flow pulsatility and risk of intestinal bleeding.

Combining Monte Carlo methods with coherent wave optics for the simulation of phase-sensitive X-ray imaging

Phase-sensitive X-ray imaging shows a high sensitivity towards electron density variations, making it well suited for imaging of soft tissue matter. However, there are still open questions about the details of the image formation process. Here, a framework for numerical simulations of phase-sensitive X-ray imaging is presented, which takes both particle- and wave-like properties of X-rays into consideration. A split approach is presented where we combine a Monte Carlo method (MC) based sample part with a wave optics simulation based propagation part, leading to a framework that takes both particle- and wave-like properties into account. The framework can be adapted to different phase-sensitive imaging methods and has been validated through comparisons with experiments for grating interferometry and propagation-based imaging. The validation of the framework shows that the combination of wave optics and MC has been successfully implemented and yields good agreement between measurements and simulations. This demonstrates that the physical processes relevant for developing a deeper understanding of scattering in the context of phase-sensitive imaging are modelled in a sufficiently accurate manner. The framework can be used for the simulation of phase-sensitive X-ray imaging, for instance for the simulation of grating interferometry or propagation-based imaging.

Correction of dephasing oscillations in matter-wave interferometry

Vibrations, electromagnetic oscillations, and temperature drifts are among the main reasons for dephasing in matter-wave interferometry. Sophisticated interferometry experiments, e.g., with ions or heavy molecules, often require integration times of several minutes due to the low source intensity or the high velocity selection. Here we present a scheme to suppress the influence of such dephasing mechanisms—especially in the low-frequency regime—by analyzing temporal and spatial particle correlations available in modern detectors. Such correlations can reveal interference properties that would otherwise be washed out due to dephasing by external oscillating signals. The method is shown experimentally in a biprism electron interferometer where a perturbing oscillation is artificially introduced by a periodically varying magnetic field. We provide a full theoretical description of the particle correlations where the perturbing frequency and amplitude can be revealed from the disturbed interferogram. The original spatial fringe pattern without the perturbation can thereby be restored. The technique can be applied to lower the general noise requirements in matter-wave interferometers. It allows for the optimization of electromagnetic shielding and decreases the efforts for vibrational or temperature stabilization.

Which side of the focal plane are you on?

Defocus blur is an indicator for the depth structure of a scene. However, given a single input image from a conventional camera one cannot distinguish between blurred objects lying in front or behind the focal plane, as they may be subject to exactly the same amount of blur. In this paper we address this limitation by exploiting coded apertures. Previous work in this area focuses on setups where the scene is placed either entirely in front or entirely behind the focal plane. We demonstrate that asymmetric apertures result in unique blurs for all distances from the camera. To exploit asymmetric apertures we propose an algorithm that can unambiguously estimate scene depth and texture from a single input image. One of the main advantages of our method is that, within the same depth range, we can work with less blurred data than in other methods. The technique is tested on both synthetic and real images.

Measurement of event-plane correlations in √s_NN=2.76 TeV lead-lead collisions with the ATLAS detector

A measurement of event-plane correlations involving two or three event planes of different order is presented as a function of centrality for 7 μb −1 Pb+Pb collision data at √s NN =2.76 TeV, recorded by the ATLAS experiment at the Large Hadron Collider. Fourteen correlators are measured using a standard event-plane method and a scalar-product method, and the latter method is found to give a systematically larger correlation signal. Several different trends in the centrality dependence of these correlators are observed. These trends are not reproduced by predictions based on the Glauber model, which includes only the correlations from the collision geometry in the initial state. Calculations that include the final-state collective dynamics are able to describe qualitatively, and in some cases also quantitatively, the centrality dependence of the measured correlators. These observations suggest that both the fluctuations in the initial geometry and the nonlinear mixing between different harmonics in the final state are important for creating these correlations in momentum space.

(Table 2) Compressional and shear wave velocites from basalts and sediments of the Atlantic and Pacific Ocean

Compressional (Vp) and shear (Vs) wave velocities have been measured to 1.0 kbar for 14 cores of well-consolidated sedimentary rock from Atlantic and Pacific sites of the Deep Sea Drilling Project. The range of VP (2.05-5.38 km/sec at 0.5 kbar) shows significant overlap with the range of oceanic layer-2 seismic velocities determined by marine refraction surveys, suggesting that sedimentary rocks may, in some regions, constitute the upper portion of layer 2. Differing linear relationships between VP and Vs for basalts and sedimentary rocks, however, may provide a method of resolving layer-2 composition. This is illustra ted for a refraction survey site on the flank of the Mid-Atlantic Ridge where layer-2 velocities agree with basalt, and two sites on the Saya de Malha Bank in the Indian Ocean where layer-2 velocities appear to represent sedimentary rock.