Computational Tools for Stereo and Light Field Photography

This thesis covers a broad part of the field of computational photography, including video stabilization and image warping techniques, introductions to light field photography and the conversion of monocular images and videos into stereoscopic 3D content. We present a user assisted technique for stereoscopic 3D conversion from 2D images. Our approach exploits the geometric structure of perspective images including vanishing points. We allow a user to indicate lines, planes, and vanishing points in the input image, and directly employ these as guides of an image warp that produces a stereo image pair. Our method is most suitable for scenes with large scale structures such as buildings and is able to skip the step of constructing a depth map. Further, we propose a method to acquire 3D light fields using a hand-held camera, and describe several computational photography applications facilitated by our approach. As the input we take an image sequence from a camera translating along an approximately linear path with limited camera rotations. Users can acquire such data easily in a few seconds by moving a hand-held camera. We convert the input into a regularly sampled 3D light field by resampling and aligning them in the spatio-temporal domain. We also present a novel technique for high-quality disparity estimation from light fields. Finally, we show applications including digital refocusing and synthetic aperture blur, foreground removal, selective colorization, and others.

Preoperative Planning of Periacetabular Osteotomy (PAO)

Pelvic osteotomies improve containment of the femoral head in cases of developmental dysplasia of the hip or in femoroacetabular impingement due to acetabular retroversion. In the evolution of osteotomies, the Ganz Periacetabular Osteotomy (PAO) is among the complex reorientation osteotomies and allows for complete mobilization of the acetabulum without compromising the integrity of the pelvic ring. For the complex reorientation osteotomies, preoperative planning of the required acetabular correction is an important step, due to the need to comprehend the three-dimensional (3D) relationship between acetabulum and femur. Traditionally, planning was performed using conventional radiographs in different projections, reducing the 3D problem to a two-dimensional one. Known disturbance variables, mainly tilt and rotation of the pelvis make assessment by these means approximate at the most. The advent of modern enhanced computation skills and new imaging techniques gave room for more sophisticated means of preoperative planning. Apart from analysis of acetabular geometry on conventional x-rays by sophisticated software applications, more accurate assessment of coverage and congruency and thus amount of correction necessary can be performed on multiplanar CT images. With further evolution of computer-assisted orthopaedic surgery, especially the ability to generate 3D models from the CT data, examiners were enabled to simulate the in vivo situation in a virtual in vitro setting. Based on this ability, different techniques have been described. They basically all employ virtual definition of an acetabular fragment. Subsequently reorientation can be simulated using either 3D calculation of standard parameters of femoroacetabular morphology, or joint contact pressures, or a combination of both. Other techniques employ patient specific implants, templates or cutting guides to achieve the goal of safe periacetabular osteotomies. This chapter will give an overview of the available techniques for planning of periacetabular osteotomy.

Computational fluid dynamics modelling in cardiovascular medicine.

This paper reviews the methods, benefits and challenges associated with the adoption and translation of computational fluid dynamics (CFD) modelling within cardiovascular medicine. CFD, a specialist area of mathematics and a branch of fluid mechanics, is used routinely in a diverse range of safety-critical engineering systems, which increasingly is being applied to the cardiovascular system. By facilitating rapid, economical, low-risk prototyping, CFD modelling has already revolutionised research and development of devices such as stents, valve prostheses, and ventricular assist devices. Combined with cardiovascular imaging, CFD simulation enables detailed characterisation of complex physiological pressure and flow fields and the computation of metrics which cannot be directly measured, for example, wall shear stress. CFD models are now being translated into clinical tools for physicians to use across the spectrum of coronary, valvular, congenital, myocardial and peripheral vascular diseases. CFD modelling is apposite for minimally-invasive patient assessment. Patient-specific (incorporating data unique to the individual) and multi-scale (combining models of different length- and time-scales) modelling enables individualised risk prediction and virtual treatment planning. This represents a significant departure from traditional dependence upon registry-based, population-averaged data. Model integration is progressively moving towards 'digital patient' or 'virtual physiological human' representations. When combined with population-scale numerical models, these models have the potential to reduce the cost, time and risk associated with clinical trials. The adoption of CFD modelling signals a new era in cardiovascular medicine. While potentially highly beneficial, a number of academic and commercial groups are addressing the associated methodological, regulatory, education- and service-related challenges.

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.

Outcomes of Active Surveillance for Ductal Carcinoma in Situ: A Computational Risk Analysis.

BACKGROUND Ductal carcinoma in situ (DCIS) is a noninvasive breast lesion with uncertain risk for invasive progression. Usual care (UC) for DCIS consists of treatment upon diagnosis, thus potentially overtreating patients with low propensity for progression. One strategy to reduce overtreatment is active surveillance (AS), whereby DCIS is treated only upon detection of invasive disease. Our goal was to perform a quantitative evaluation of outcomes following an AS strategy for DCIS. METHODS Age-stratified, 10-year disease-specific cumulative mortality (DSCM) for AS was calculated using a computational risk projection model based upon published estimates for natural history parameters, and Surveillance, Epidemiology, and End Results data for outcomes. AS projections were compared with the DSCM for patients who received UC. To quantify the propagation of parameter uncertainty, a 95% projection range (PR) was computed, and sensitivity analyses were performed. RESULTS Under the assumption that AS cannot outperform UC, the projected median differences in 10-year DSCM between AS and UC when diagnosed at ages 40, 55, and 70 years were 2.6% (PR = 1.4%-5.1%), 1.5% (PR = 0.5%-3.5%), and 0.6% (PR = 0.0%-2.4), respectively. Corresponding median numbers of patients needed to treat to avert one breast cancer death were 38.3 (PR = 19.7-69.9), 67.3 (PR = 28.7-211.4), and 157.2 (PR = 41.1-3872.8), respectively. Sensitivity analyses showed that the parameter with greatest impact on DSCM was the probability of understaging invasive cancer at diagnosis. CONCLUSION AS could be a viable management strategy for carefully selected DCIS patients, particularly among older age groups and those with substantial competing mortality risks. The effectiveness of AS could be markedly improved by reducing the rate of understaging.

SrGAP2-Dependent Integration of Membrane Geometry and Slit-Robo-Repulsive Cues Regulates Fibroblast Contact Inhibition of Locomotion

Migrating fibroblasts undergo contact inhibition of locomotion (CIL), a process that was discovered five decades ago and still is not fully understood at the molecular level. We identify the Slit2-Robo4-srGAP2 signaling network as a key regulator of CIL in fibroblasts. CIL involves highly dynamic contact protrusions with a specialized actin cytoskeleton that stochastically explore cell-cell overlaps between colliding fibroblasts. A membrane curvature-sensing F-BAR domain pre-localizes srGAP2 to protruding edges and terminates their extension phase in response to cell collision. A FRET-based biosensor reveals that Rac1 activity is focused in a band at the tip of contact protrusions, in contrast to the broad activation gradient in contact-free protrusions. SrGAP2 specifically controls the duration of Rac1 activity in contact protrusions, but not in contact-free protrusions. We propose that srGAP2 integrates cell edge curvature and Slit-Robo-mediated repulsive cues to fine-tune Rac1 activation dynamics in contact protrusions to spatiotemporally coordinate CIL.

Effects of Strand Geometry on Selected Properties of Long-Stand Structural Composite Lumber Made from Northeastern Hardwoods

Structural composite lumber (SCL) products often possess significantly higher design values than the top grades of solid lumber, making it a popular choice for both residential and commercial applications. The enhanced mechanical properties of SCL are mainly due to defect randomization and densification of the wood fiber, both largely functions of the size, shape and composition (species) of the wood element. Traditionally, SCL manufacturers have used thin, rectangular elements produced from either moderate density softwoods or low density hardwoods. Higher density hardwood species have been avoided, as they require higher pressures to adequately densify and consolidate the wood furnish. These higher pressures can lead to increased manufacturing costs, damage to the wood fiber and/or a product that is too dense, making it heavy and unreceptive to common mechanical fastening techniques. In the northeastern United States high density, diffuse-porous hardwoods (such as maple, beech and birch) are abundant. Use of these species as primary furnish for a SCL product may allow for a competitive advantage in terms of resource cost against products that rely on veneer grade logs. Proximity to this abundant and relatively inexpensive resource may facilitate entry of SCL production facilities in the northeastern United States, where currently none exist. However, modifications to current strand sizes, geometries or production techniques will likely be required to allow for use of these species. A new SCL product concept has been invented allowing for use of these high density hardwoods. The product, referred to as long-strand structural composite lumber (LSSCL), uses strands of significantly larger cross sectional areas and volumes than existing SCL products. In spite of the large strand size, satisfactory consolidation is achieved without excessive densification of the wood fiber through use of a symmetrical strand geometric cross-section. LSSCL density is similar to that of existing SCL products, but is due mainly to the inherent density of the species, rather than through densification. An experiment was designed and conducted producing LSSCL from both large (7/16”) and small (1/4”) strands, of both square and triangular geometric cross sections. Testing results indicate that the large, triangular strands produce LSSCL beams with projected design values of: Modulus of elasticity (MOEapp) – 1,750,000 psi; Allowable bending stress (Fb) – 2750 psi; Allowable shear stress (Fv) – 260 psi. Several modifications are recommended which may lead to improvement of these values, likely allowing for competition against existing SCL products.

The Effects of Credit Risk on Dynamic Portfolio Management: A New Computational Approach

The study investigates the role of credit risk in a continuous time stochastic asset allocation model, since the traditional dynamic framework does not provide credit risk flexibility. The general model of the study extends the traditional dynamic efficiency framework by explicitly deriving the optimal value function for the infinite horizon stochastic control problem via a weighted volatility measure of market and credit risk. The model's optimal strategy was then compared to that obtained from a benchmark Markowitz-type dynamic optimization framework to determine which specification adequately reflects the optimal terminal investment returns and strategy under credit and market risks. The paper shows that an investor's optimal terminal return is lower than typically indicated under the traditional mean-variance framework during periods of elevated credit risk. Hence I conclude that, while the traditional dynamic mean-variance approach may indicate the ideal, in the presence of credit-risk it does not accurately reflect the observed optimal returns, terminal wealth and portfolio selection strategies.

An integrated molecular biology and computational approach to CYP1A1 expression in human brain

The cytochromes P450 comprise a superfamily of heme-containing mono-oxygenases. These enzymes metabolize numerous xenobiotics, but also play a role in metabolism of endogenous compounds. The P450 1A1 enzyme generally metabolizes polycyclic aromatic hydrocarbons, and its expression can be induced by aryl hydrocarbon receptor (AhR) activation. CYP1A1 is an exception to the generality that the majority of CYPs demonstrate highest expression in liver; CYP1Al is present in numerous extrahepatic tissues, including brain. This P450 has been observed in two forms, wildtype (WT) and brain variant (BV), arising from alternatively spliced mRNA transcripts. The CYP1A1 BV mRNA presented an exon deletion and was detected in human brain but not liver tissue of the same individuals. ^ Quantitative PCR analyses were performed to determine CYP1A1 WT and BV transcript expression levels in normal, bipolar disorder or schizophrenic groups. In our samples, we show that CYP1A1 BV mRNA, when present, is found alongside the full-length form. Furthermore, we demonstrate a significant decrease in expression of CYP1A1 in patients with bipolar disorder or schizophrenia. The expression level was not influenced by post-mortem interval, tissue pH, age, tobacco use, or lifetime antipsychotic medication load. ^ There is no indication of increased brain CYP1A1 expression in normal smokers versus non-smokers in these samples. We observed slightly increased CYP1A1 expression only in bipolar and schizophrenic smokers versus non-smokers. This may be indicative of complex interactions between neuronal chemical environments and AhR-mediated CYP1A1 induction in brain. ^ Structural homology modeling demonstrated that P450 1A1 BV has several alterations to positions/orientations of substrate recognition site residues compared to the WT isoform. Automated substrate docking was employed to investigate the potential binding of neurological signaling molecules and neurotropic drugs, as well as to differentiate specificities of the two P450 1A1 isoforms. We consistently observed that the BV isoform produced energetically favorable substrate dockings in orientations not observed for the same substrate in the WT isoform. These results demonstrated that structural differences, namely an expanded substrate access channel and active site, confer greater capacity for unique compound docking positions suggesting a metabolic profile distinct from the wildtype form for these test compounds. ^

Present-day King George Island ice cap geometry measured by ground penetrating radar in 2008/09

We present a consistent data set for the ice thickness, the bedrock topography and the ice surface topography of the King George Island ice cap (Arctowski Icefield and the adjacent central part). The data set is composed of groundbased and airborne Ground Penetrating Radar (GPR) and differential GPS (DGPS) measurements, obtained during several field campaigns. The data set incorporates groundbased measurements in the safely accessible inner parts and airborne measurements in the heavily crevassed coastal areas of the ice cap. In particular, the inclusion of airborne GPR measurements with the 30MHz BGR-P30-System developed at the Institute of Geophysics (University of Münster) completes the picture of the ice geometry substantially. The compiled digital elevation model of the bedrock shows a rough, highly variable topography with pronounced valleys, ridges, and troughs. The mean ice thickness is approx. 238m, with a maximum value of approx. 400m in the surveyed area. Noticeable are bounded areas in the bedrock topography below sea level where marine based ice exists.