Automatic multi-resolution shape modeling of multi-organ structures.

Point Distribution Models (PDM) are among the most popular shape description techniques and their usefulness has been demonstrated in a wide variety of medical imaging applications. However, to adequately characterize the underlying modeled population it is essential to have a representative number of training samples, which is not always possible. This problem is especially relevant as the complexity of the modeled structure increases, being the modeling of ensembles of multiple 3D organs one of the most challenging cases. In this paper, we introduce a new GEneralized Multi-resolution PDM (GEM-PDM) in the context of multi-organ analysis able to efficiently characterize the different inter-object relations, as well as the particular locality of each object separately. Importantly, unlike previous approaches, the configuration of the algorithm is automated thanks to a new agglomerative landmark clustering method proposed here, which equally allows us to identify smaller anatomically significant regions within organs. The significant advantage of the GEM-PDM method over two previous approaches (PDM and hierarchical PDM) in terms of shape modeling accuracy and robustness to noise, has been successfully verified for two different databases of sets of multiple organs: six subcortical brain structures, and seven abdominal organs. Finally, we propose the integration of the new shape modeling framework into an active shape-model-based segmentation algorithm. The resulting algorithm, named GEMA, provides a better overall performance than the two classical approaches tested, ASM, and hierarchical ASM, when applied to the segmentation of 3D brain MRI.

Fast Prediction of Femoral Biomechanics Using Supervised Machine Learning and Statistical Shape Modeling

Finite element (FE) analysis is an important computational tool in biomechanics. However, its adoption into clinical practice has been hampered by its computational complexity and required high technical competences for clinicians. In this paper we propose a supervised learning approach to predict the outcome of the FE analysis. We demonstrate our approach on clinical CT and X-ray femur images for FE predictions ( FEP), with features extracted, respectively, from a statistical shape model and from 2D-based morphometric and density information. Using leave-one-out experiments and sensitivity analysis, comprising a database of 89 clinical cases, our method is capable of predicting the distribution of stress values for a walking loading condition with an average correlation coefficient of 0.984 and 0.976, for CT and X-ray images, respectively. These findings suggest that supervised learning approaches have the potential to leverage the clinical integration of mechanical simulations for the treatment of musculoskeletal conditions.

Shape model, reference system definition, and cartographic mapping standards for comet 67P/Churyumov-Gerasimenko Stereo-photogrammetric analysis of Rosetta/OSIRIS image data

We analyzed more than 200 OSIRIS NAC images with a pixel scale of 0.9-2.4 m/pixel of comet 67P/Churyumov-Gerasimenko (67P) that have been acquired from onboard the Rosetta spacecraft in August and September 2014 using stereo-photogrammetric methods (SPG). We derived improved spacecraft position and pointing data for the OSIRIS images and a high-resolution shape model that consists of about 16 million facets (2 m horizontal sampling) and a typical vertical accuracy at the decimeter scale. From this model, we derive a volume for the northern hemisphere of 9.35 km(3) +/- 0.1 km(3). With the assumption of a homogeneous density distribution and taking into account the current uncertainty of the position of the comet's center-of-mass, we extrapolated this value to an overall volume of 18.7 km(3) +/- 1.2 km(3), and, with a current best estimate of 1.0 X 10(13) kg for the mass, we derive a bulk density of 535 kg/m(3) +/- 35 kg/m(3). Furthermore, we used SPG methods to analyze the rotational elements of 67P. The rotational period for August and September 2014 was determined to be 12.4041 +/- 0.0004 h. For the orientation of the rotational axis (z-axis of the body-fixed reference frame) we derived a precession model with a half-cone angle of 0.14 degrees, a cone center position at 69.54 degrees/64.11 degrees (RA/Dec J2000 equatorial coordinates), and a precession period of 10.7 days. For the definition of zero longitude (x-axis orientation), we finally selected the boulder-like Cheops feature on the big lobe of 67P and fixed its spherical coordinates to 142.35 degrees right-hand-rule eastern longitude and -0.28 degrees latitude. This completes the definition of the new Cheops reference frame for 67P. Finally, we defined cartographic mapping standards for common use and combined analyses of scientific results that have been obtained not only within the OSIRIS team, but also within other groups of the Rosetta mission.

Facies distribution in CRP sediment cores from the Ross Sea, Antarctica

Clasts from the Cape Roberts Project cores CRP-2/2A and CRP-3 provide indications of glacially influenced depositional environments in Oligocene and Miocene strata in the western Victoria Land Basin, Antarctica. CRP-2/2A is interpreted to represent strongly glacially influenced, unconformity bound depositional sequences produced by repeated advance and retreat of floating and grounded ice across the shelf. A similar interpretation is extended to the upper 330 meters of the CRP-3 core, but the lower part of the core records shallow marine deposition with significantly less glacial influence. Clast shape analysis from selected coarse-grained facies throughout the cored interval indicates that most clasts are glacially sourced, with little distinction between diamictite and conglomeratic facies. Three dimensional clast fabric analysis from units immediately above sequence boundaries generally display weak or random fabrics and do not suggest that grounded ice actually reached the drillsite at these intervals. Striated and outsized clasts present in fine-grained lithofacies throughout the cores provide further evidence of sub-glacially transported sediment and iceberg rafting. The distribution of these striated and out-sized clasts indicate that a significant glacial influence persisted through most of the time represented by the cores with glaciers actively calving at sea-level introducing ice-berg rafted glacial debris even in the earliest Oligocene.

Distribution of phytolithes in CRP sediment cores from the Ross Sea, Antarctica

Phytoliths (siliceous plant microfossils) have been recovered from Cenozoic sediments (c. 34 to 17 Ma) in the CRP-2/2A and CRP-3 drillholes cored off Cape Roberts, Victoria Land Basin, Antarctica. The phytolith assemblages are sparse, but well-preserved and dominated by spherical forms similar to those of modern trees or shrubs. Rare phytoliths comparable to modern grass forms are also present. However, due to the paucity of phytolith data, any interpretations made are necessarily tentative. The assemblages of CRP-2/2A and the upper c. 250 m of CRP-3 are interpreted as representing a predominantly woody vegetation, including Nothofagus and Libocedrus with local areas of grass in the more exposed locations. A cool climate is interpreted to have prevailed throughout both cores. However, beneath c. 250 metres below sea floor in CRP-3, the dominant woody vegetation is supplemented by pockets of Palmae, ?Proteaceae and 'warm' climate grasses. This association represents vegetation growth in sheltered, moist sites - possibly north-facing mid-slopes or the coastal fringe. It may also represent remnant vegetation that grew in moist, temperate conditions during the Middle to Late Eocene, previously interpreted from the Southern McMurdo Sound erratics and lower part of the CIROS-1 drillhole. The phytolith analysis compares well to the terrestrial palynomorph record from both cores and provides additional independent taxonomic and climatic interpretations.

Maps and acoustic profiles of free gas distribution, Western Baltic Sea and Kattgat-Skagerrak area

A shallow gas depth-contour map covering the Skagerrak-western Baltic Sea region has been constructed using a relatively dense grid of existing shallow seismic lines. The digital map is stored as an ESRI shape file in order to facilitate comparison with other data from the region. Free gas usually occurs in mud and sandy mud but is observed only when sediment thickness exceeds a certain threshold value, depending on the water depth of the area in question. Gassy sediments exist at all water depths from approx. 20 m in the coastal waters of the Kattegat to 360 m in the Skagerrak. In spite of the large difference in water depths, the depth of free gas below seabed varies only little within the region, indicating a relatively fast movement of methane in the gas phase towards the seabed compared to the rate of diffusion of dissolved methane. Seeps of old microbial methane occur in the northern Kattegat where a relatively thin cover of sandy sediments exists over shallow, glacially deformed Pleistocene marine sediments. Previous estimates of total methane escape from the area may be correct but the extrapolation of local methane seepage rate data to much larger areas on the continental shelf is probably not justified. Preliminary data on porewater chemistry were compared with the free gas depth contours in the Aarhus Bay area, which occasionally suffers from oxygen deficiency, in order to examine if acoustic gas mapping may be used for monitoring the condition of the bay.

Ice thickness distribution and glacier bed of Hallstätter Gletscher

Hallstätter Glacier is the northernmost glacier of Austria. Appendant to the northern Limestone Alps, the glacier is located at 47°28'50'' N, 13°36'50'' E in the Dachstein-region. At the same time with its advance linked to the Little Ice Age (LIA), research on changes in size and mass of Hallstätter glacier was started in 1842 by Friedrich Simony. He observed and documented the glacier retreat related to its last maximum extension in 1856. In addition, Hallstätter Glacier is a subject to scientific research to date. In this thesis methods and results of ongoing mass balance measurements are presented and compared to long term volume changes and meteorological observations. The current mass balance monitoring programm using the direct glaciological method was started 2006. In this context, 2009 the ice thickness was measured with ground penetrating radar. The result are used with digital elevation models reconstucted from historical maps and recent digital elevation models to calculate changes in shape and volume of Hallstätter Glacier. Based on current meteorological measurements near the glacier and longtime homogenized climate data provided by HISTALP, time series of precipitation and temperature beginning at the LIA are produced. These monthly precipitation and monthly mean temperature data are used to compare results of a simple degree day model with the volume change calculated from the difference of the digital elevation models. The two years of direct mass balance measurements are used to calibrate the degree day model. A number of possible future scenarios are produced to indicate prospective changes. Within the 150-year-period between 1856 and 2007 the Hallstätter Glacier lost 1940 meters of its length and 2.23 km**2 in area. 37% of the initial volume of 1856 remained. This retreat came along with a change in climate. The application of a running avarage of 30 years shows an increase in precipitation of 18.5% and a warming of 1.3°C near the glacier between 1866 and 1993. The mass loss was continued in the hydrological years 2006/2007 and 2007/2008 showing mean specific mass balance of -376 mm and -700 mm, respectively. Applying a temperature correction for the different minimum elevations of the glacier, the degree day approach based on the two measured mass balances can reproduce sign and order of magnitude of the volume change of Hallstätter Glacier since 1856. Nevertheless, the relative deviation is significant. Future scenarios show, that 30% of the entire glacier volume remains after subtracting the elevation changes between the digital elevation models of 2002 and 2007 ten times from the surface of 2007. The past and present mass changes of Hallstätter Glacier are showing a retreating glacier as a consequence of rising temperatures. Due to high precepitation, increased with previous warming, the Hallstätter Glacier can and will exist in lower elevation compared to inner alpine glaciers.

Climate sensitivity across marine domains of life: limits to evolutionary adaptation shape species interactions, supplementary material

Organisms in all domains, Archaea, Bacteria, and Eukarya will respond to climate change with differential vulnerabilities resulting in shifts in species distribution, coexistence, and interactions. The identification of unifying principles of organism functioning across all domains would facilitate a cause and effect understanding of such changes and their implications for ecosystem shifts. For example, the functional specialization of all organisms in limited temperature ranges leads us to ask for unifying functional reasons. Organisms also specialize in either anoxic or various oxygen ranges, with animals and plants depending on high oxygen levels. Here, we identify thermal ranges, heat limits of growth, and critically low (hypoxic) oxygen concentrations as proxies of tolerance in a meta-analysis of data available for marine organisms, with special reference to domain-specific limits. For an explanation of the patterns and differences observed, we define and quantify a proxy for organismic complexity across species from all domains. Rising complexity causes heat (and hypoxia) tolerances to decrease from Archaea to Bacteria to uni- and then multicellular Eukarya. Within and across domains, taxon-specific tolerance limits likely reflect ultimate evolutionary limits of its species to acclimatization and adaptation. We hypothesize that rising taxon-specific complexities in structure and function constrain organisms to narrower environmental ranges. Low complexity as in Archaea and some Bacteria provide life options in extreme environments. In the warmest oceans, temperature maxima reach and will surpass the permanent limits to the existence of multicellular animals, plants and unicellular phytoplankter. Smaller, less complex unicellular Eukarya, Bacteria, and Archaea will thus benefit and predominate even more in a future, warmer, and hypoxic ocean.

Relative Distribution Methods in Stata

The concept of the relative density seems like a fruitful nonparametric approach to studying distributional differences between groups (Handcock and Morris 1999), yet it appears that the technique has gone more or less unnoticed in applied social science research. A scarcity of canned software might be one of the reasons the method is underutilized. Therefore, I present a new Stata command called reldist to plot the relative density, decompose distributional differences into location and shape effects, and compute relative distribution summary measures. The command is illustrated by an application comparing earnings by sex.

Application of active heat pulse method with fiber optic temperature sensing for estimation of wetting bulbs and water distribution in drip emitters.

Through the use of the Distributed Fiber Optic Temperature Measurement (DFOT) method, it is possible to measure the temperature in small intervals (on the order of centimeters) for long distances (on the order of kilometers) with a high temporal frequency and great accuracy. The heat pulse method consists of applying a known amount of heat to the soil and monitoring the temperature evolution, which is primarily dependent on the soil moisture content. The use of both methods, which is called the active heat pulse method with fiber optic temperature sensing (AHFO), allows accurate soil moisture content measurements. In order to experimentally study the wetting patterns, i.e. shape, size, and the water distribution, from a drip irrigation emitter, a soil column of 0.5 m of diameter and 0.6 m high was built. Inside the column, a fiber optic cable with a stainless steel sheath was placed forming three concentric helixes of diameters 0.2 m, 0.4 m and 0.6 m, leading to a 148 measurement point network. Before, during, and after the irrigation event, heat pulses were performed supplying electrical power of 20 W/m to the steel. The soil moisture content was measured with a capacitive sensor in one location at depths of 0.1 m, 0.2 m, 0.3 m and 0.4 m during the irrigation. It was also determined by the gravimetric method in several locations and depths before and right after the irrigation. The emitter bulb dimensions and shape evolution was satisfactorily measured during infiltration. Furthermore, some bulb's characteristics difficult to predict (e.g. preferential flow) were detected. The results point out that the AHFO is a useful tool to estimate the wetting pattern of drip irrigation emitters in soil columns and show a high potential for its use in the field.

Optimization of curing cycle in carbon fiber-reinforced laminates: Void distribution and mechanical properties

A strategy is presented to optimize out-of-autoclave processing of quasi-isotropic carbon fiber-reinforced laminates. Square panels of 4.6 mm nominal thickness with very low porosity ð6 0:2%Þ were manufactured by compression molding at low pressure (0.2 MPa) by careful design of the temperature cycle to maximize the processing window. The mechanisms of void migration during processing were ascertained by means of X-ray microtomography and the effect of ply clustering on porosity and on void shape was explained. Finally, the effect of porosity and ply clustering on the compressive strength before and after impact was studied.

Computer simulation of random parckings for self-similar particle size distributions in soil and granular materials: porosity and pore size distribution

A 2D computer simulation method of random packings is applied to sets of particles generated by a self-similar uniparametric model for particle size distributions (PSDs) in granular media. The parameter p which controls the model is the proportion of mass of particles corresponding to the left half of the normalized size interval [0,1]. First the influence on the total porosity of the parameter p is analyzed and interpreted. It is shown that such parameter, and the fractal exponent of the associated power scaling, are efficient packing parameters, but this last one is not in the way predicted in a former published work addressing an analogous research in artificial granular materials. The total porosity reaches the minimum value for p = 0.6. Limited information on the pore size distribution is obtained from the packing simulations and by means of morphological analysis methods. Results show that the range of pore sizes increases for decreasing values of p showing also different shape in the volume pore size distribution. Further research including simulations with a greater number of particles and image resolution are required to obtain finer results on the hierarchical structure of pore space.

Nuclear Pore Complex Number and Distribution throughout the Saccharomyces cerevisiae Cell Cycle by Three-Dimensional Reconstruction from Electron Micrographs of Nuclear Envelopes

The number of nuclear pore complexes (NPCs) in individual nuclei of the yeast Saccharomyces cerevisiae was determined by computer-aided reconstruction of entire nuclei from electron micrographs of serially sectioned cells. Nuclei of 32 haploid cells at various points in the cell cycle were modeled and found to contain between 65 and 182 NPCs. Morphological markers, such as cell shape and nuclear shape, were used to determine the cell cycle stage of the cell being examined. NPC number was correlated with cell cycle stage to reveal that the number of NPCs increases steadily, beginning in G1-phase, suggesting that NPC assembly occurs continuously throughout the cell cycle. However, the accumulation of nuclear envelope observed during the cell cycle, indicated by nuclear surface area, is not continuous at the same rate, such that the density of NPCs per unit area of nuclear envelope peaks in apparent S-phase cells. Analysis of the nuclear envelope reconstructions also revealed no preferred NPC-to-NPC distance. However, NPCs were found in large clusters over regions of the nuclear envelope. Interestingly, clusters of NPCs were most pronounced in early mitotic nuclei and were found to be associated with the spindle pole bodies, but the functional significance of this association is unknown.

Cultured melanocytes from dilute mutant mice exhibit dendritic morphology and altered melanosome distribution

Mutant alleles at the dilute unconventional myosin heavy chain locus cause diluted coat color, opisthotonic seizures, and death. The dilute coat color phenotype is caused by irregular clumping of pigment in the hair, but amounts of melanin are unchanged from wild-type controls. The melanocyte phenotype has been described as adendritic, since hair bulb and Harderian gland melanocytes appear to be rounded in tissue sections. These observations do not exclude the possibility that the processes lack pigment, since the melanocyte shape was judged by the distribution of melanin. We have tested this hypothesis by culturing primary melanocytes from dilute mutant and wild-type mice. The mutant melanocytes do not lack processes; instead, they exhibit a concentrated perinuclear distribution of melanosomes, while wild-type melanocytes have a very uniform cytoplasmic distribution of melanosomes. Electron micrographs show no detectable differences in melanosome morphology or maturation between dilute and wild-type melanocytes. Immunofluorescence experiments indicate that the dilute protein is concentrated in regions of the cytoplasm that contain melanosomes. These experiments show that the dilute myosin is necessary for the localization of melanosomes, either by active transport or tethering.