Coltop3D: A new software for structural analysis with high resolution 3D point clouds and DEM

Coltop3D is a software that performs structural analysis by using digital elevation model (DEM) and 3D point clouds acquired with terrestrial laser scanners. A color representation merging slope aspect and slope angle is used in order to obtain a unique code of color for each orientation of a local slope. Thus a continuous planar structure appears in a unique color. Several tools are included to create stereonets, to draw traces of discontinuities, or to compute automatically density stereonet. Examples are shown to demonstrate the efficiency of the method.

Stochastic time-concentration activity models for cytotoxicity in 3D brain cell cultures.

BACKGROUND: In vitro aggregating brain cell cultures containing all types of brain cells have been shown to be useful for neurotoxicological investigations. The cultures are used for the detection of nervous system-specific effects of compounds by measuring multiple endpoints, including changes in enzyme activities. Concentration-dependent neurotoxicity is determined at several time points. METHODS: A Markov model was set up to describe the dynamics of brain cell populations exposed to potentially neurotoxic compounds. Brain cells were assumed to be either in a healthy or stressed state, with only stressed cells being susceptible to cell death. Cells may have switched between these states or died with concentration-dependent transition rates. Since cell numbers were not directly measurable, intracellular lactate dehydrogenase (LDH) activity was used as a surrogate. Assuming that changes in cell numbers are proportional to changes in intracellular LDH activity, stochastic enzyme activity models were derived. Maximum likelihood and least squares regression techniques were applied for estimation of the transition rates. Likelihood ratio tests were performed to test hypotheses about the transition rates. Simulation studies were used to investigate the performance of the transition rate estimators and to analyze the error rates of the likelihood ratio tests. The stochastic time-concentration activity model was applied to intracellular LDH activity measurements after 7 and 14 days of continuous exposure to propofol. The model describes transitions from healthy to stressed cells and from stressed cells to death. RESULTS: The model predicted that propofol would affect stressed cells more than healthy cells. Increasing propofol concentration from 10 to 100 μM reduced the mean waiting time for transition to the stressed state by 50%, from 14 to 7 days, whereas the mean duration to cellular death reduced more dramatically from 2.7 days to 6.5 hours. CONCLUSION: The proposed stochastic modeling approach can be used to discriminate between different biological hypotheses regarding the effect of a compound on the transition rates. The effects of different compounds on the transition rate estimates can be quantitatively compared. Data can be extrapolated at late measurement time points to investigate whether costs and time-consuming long-term experiments could possibly be eliminated.

Comparability of clinical wear measurements by optical 3D laser scanning in two different centers.

OBJECTIVE: The purpose of this study was to compare the use of different variables to measure the clinical wear of two denture tooth materials in two analysis centers. METHODS: Twelve edentulous patients were provided with full dentures. Two different denture tooth materials (experimental material and control) were placed randomly in accordance with the split-mouth design. For wear measurements, impressions were made after an adjustment phase of 1-2 weeks and after 6, 12, 18, and 24 months. The occlusal wear of the posterior denture teeth of 11 subjects was assessed in two study centers by use of plaster replicas and 3D laser-scanning methods. In both centers sequential scans of the occlusal surfaces were digitized and superimposed. Wear was described by use of four different variables. Statistical analysis was performed after log-transformation of the wear data by use of the Pearson and Lin correlation and by use of a mixed linear model. RESULTS: Mean occlusal vertical wear of the denture teeth after 24 months was between 120μm and 212μm, depending on wear variable and material. For three of the four variables, wear of the experimental material was statistically significantly less than that of the control. Comparison of the two study centers, however, revealed correlation of the wear variables was only moderate whereas strong correlation was observed among the different wear variables evaluated by each center. SIGNIFICANCE: Moderate correlation was observed for clinical wear measurements by optical 3D laser scanning in two different study centers. For the two denture tooth materials, wear measurements limited to the attrition zones led to the same qualitative assessment.

Imaging and quantifying salt-tracer transport in a riparian groundwater system by means of 3D ERT monitoring

Determining groundwater flow paths of infiltrated river water is necessary for studying biochemical processes in the riparian zone, but their characterization is complicated by strong temporal and spatial heterogeneity. We investigated to what extent repeat 3D surface electrical resistance tomography (ERT) can be used to monitor transport of a salt-tracer plume under close to natural gradient conditions. The aim is to estimate groundwater flow velocities and pathways at a site located within a riparian groundwater system adjacent to the perialpine Thur River in northeastern Switzerland. Our ERT time-lapse images provide constraints on the plume's shape, flow direction, and velocity. These images allow the movement of the plume to be followed for 35 m. Although the hydraulic gradient is only 1.43 parts per thousand, the ERT time-lapse images demonstrate that the plume's center of mass and its front propagate with velocities of 2x10(-4) m/s and 5x10(-4) m/s, respectively. These velocities are compatible with groundwater resistivity monitoring data in two observation wells 5 m from the injection well. Five additional sensors in the 5-30 m distance range did not detect the plume. Comparison of the ERT time-lapse images with a groundwater transport model and time-lapse inversions of synthetic ERT data indicate that the movement of the plume can be described for the first 6 h after injection by a uniform transport model. Subsurface heterogeneity causes a change of the plume's direction and velocity at later times. Our results demonstrate the effectiveness of using time-lapse 3D surface ERT to monitor flow pathways in a challenging perialpine environment over larger scales than is practically possible with crosshole 3D ERT.

Flow-R, a model for susceptibility mapping of debris flows and other gravitational hazards at a regional scale: some case studies

Every year, debris flows cause huge damage in mountainous areas. Due to population pressure in hazardous zones, the socio-economic impact is much higher than in the past. Therefore, the development of indicative susceptibility hazard maps is of primary importance, particularly in developing countries. However, the complexity of the phenomenon and the variability of local controlling factors limit the use of processbased models for a first assessment. A debris flow model has been developed for regional susceptibility assessments using digital elevation model (DEM) with a GIS-based approach.. The automatic identification of source areas and the estimation of debris flow spreading, based on GIS tools, provide a substantial basis for a preliminary susceptibility assessment at a regional scale. One of the main advantages of this model is its workability. In fact, everything is open to the user, from the data choice to the selection of the algorithms and their parameters. The Flow-R model was tested in three different contexts: two in Switzerland and one in Pakistan, for indicative susceptibility hazard mapping. It was shown that the quality of the DEM is the most important parameter to obtain reliable results for propagation, but also to identify the potential debris flows sources.

Brain damage in methylmalonic aciduria: 2-methylcitrate induces cerebral ammonium accumulation and apoptosis in 3D organotypic brain cell cultures.

BACKGROUND: Methylmalonic aciduria is an inborn error of metabolism characterized by accumulation of methylmalonate (MMA), propionate and 2-methylcitrate (2-MCA) in body fluids. Early diagnosis and current treatment strategies aimed at limiting the production of these metabolites are only partially effective in preventing neurological damage. METHODS: To explore the metabolic consequences of methylmalonic aciduria on the brain, we used 3D organotypic brain cell cultures from rat embryos. We challenged the cultures at two different developmental stages with 1 mM MMA, propionate or 2-MCA applied 6 times every 12 h. In a dose-response experiment cultures were challenged with 0.01, 0.1, 0.33 and 1 mM 2-MCA. Immunohistochemical staining for different brain cell markers were used to assess cell viability, morphology and differentiation. Significant changes were validated by western blot analysis. Biochemical markers were analyzed in culture media. Apoptosis was studied by immunofluorescence staining and western blots for activated caspase-3. RESULTS: Among the three metabolites tested, 2-MCA consistently produced the most pronounced effects. Exposure to 2-MCA caused morphological changes in neuronal and glial cells already at 0.01 mM. At the biochemical level the most striking result was a significant ammonium increase in culture media with a concomitant glutamine decrease. Dose-response studies showed significant and parallel changes of ammonium and glutamine starting from 0.1 mM 2-MCA. An increased apoptosis rate was observed by activation of caspase-3 after exposure to at least 0.1 mM 2-MCA. CONCLUSION: Surprisingly, 2-MCA, and not MMA, seems to be the most toxic metabolite in our in vitro model leading to delayed axonal growth, apoptosis of glial cells and to unexpected ammonium increase. Morphological changes were already observed at 2-MCA concentrations as low as 0.01 mM. Increased apoptosis and ammonium accumulation started at 0.1 mM thus suggesting that ammonium accumulation is secondary to cell suffering and/or cell death. Local accumulation of ammonium in CNS, that may remain undetected in plasma and urine, may therefore play a key role in the neuropathogenesis of methylmalonic aciduria both during acute decompensations and in chronic phases. If confirmed in vivo, this finding might shift the current paradigm and result in novel therapeutic strategies.

How Wages and Employment Adjust to Trade Liberalization: Quasi-Experimental Evidence from Austria

We study the response of regional employment and nominal wages to trade liberalization, exploiting the natural experiment provided by the opening of Central and Eastern European markets after the fall of the Iron Curtain in 1990. Using data for Austrian municipalities, we examine differential pre- and post-1990 wage and employment growth rates between regions bordering the formerly communist economies and interior regions. If the "border regions" are defined narrowly, within a band of less than 50 kilometers, we can identify statistically significant liberalization effects on both employment and wages. While wages responded earlier than employment, the employment effect over the entire adjustment period is estimated to be around three times as large as the wage effect. The implied slope of the regional labor supply curve can be replicated in an economic geography model that features obstacles to labor migration due to immobile housing and to heterogeneous locational preferences.

Mòdul de visualització 3D d’ ajuda al diagnòstic de les neoplàsies rectals

El diagnòstic mitjançant la imatge mèdica s’ha convertit en una eina fonamental en la pràctica clínica, permet entre altres coses, reconstruir a partir d’un conjunt d’imatges 2D, obtingudes a partir d’aparells de captació, qualsevol part de l’organisme d’un pacient i representar-lo en un model 3D. Sobre aquest model 3D poden realitzar-se diferents operacions que faciliten el diagnòstic i la presa de decisions als especialistes. El projecte que es presenta forma part del desenvolupament de la plataforma informàtica de visualització i tractament de dades mèdiques, anomenada Starviewer, que desenvolupen conjuntament el laboratori de Gràfics i Imatge (GiLab) de la Universitat de Girona i l’ Institut de Diagnòstic per la Imatge (IDI) de l’Hospital Josep Trueta de Girona. En particular, en aquest projecte es centra en el diagnòstic del càncer colorectal i el desenvolupament de mètodes i tècniques de suport al seu diagnòstic. Els dos punts claus en el tractament d’aqueta patologia són: la detecció de les lesions I l’estudi de l’evolució d’aquestes lesions, una vegada s’ha iniciat el tractament tumoral. L’objectiu principal d’aquest projecte és implementar i integrar en la plataforma Starviewer les tècniques de visualització i processament de dades necessàries per donar suport als especialistes en el diagnòstic de les lesions del colon. Donada la dificultat en el processament de les dades reals del budell ens proposem: dissenyar i implementar un sistema per crear models sintètics del budell; estudiar, implementar i avaluar les tècniques de processament d’imatge que calen per segmentar lesions de budell; dissenyar i implementar un sistema d’exploració del budell iintegrar de tots els mòduls implementats en la plataforma starviewer

High throughput 3D super-resolution microscopy reveals Caulobacter crescentus in vivo Z-ring organization.

We created a high-throughput modality of photoactivated localization microscopy (PALM) that enables automated 3D PALM imaging of hundreds of synchronized bacteria during all stages of the cell cycle. We used high-throughput PALM to investigate the nanoscale organization of the bacterial cell division protein FtsZ in live Caulobacter crescentus. We observed that FtsZ predominantly localizes as a patchy midcell band, and only rarely as a continuous ring, supporting a model of "Z-ring" organization whereby FtsZ protofilaments are randomly distributed within the band and interact only weakly. We found evidence for a previously unidentified period of rapid ring contraction in the final stages of the cell cycle. We also found that DNA damage resulted in production of high-density continuous Z-rings, which may obstruct cytokinesis. Our results provide a detailed quantitative picture of in vivo Z-ring organization.

Evolution of a debris flow channel monitored using a 3D terrestrial laser scanner

Like numerous torrents in mountainous regions, the Illgraben creek (canton of Wallis, SW Switzerland) produces almost every year several debris flows. The total area of the active catchment is only 4.7 km², but large events ranging from 50'000 to 400'000 m³ are common (Zimmermann 2000). Consequently, the pathway of the main channel often changes suddenly. One single event can for instance fill the whole river bed and dig new several-meters-deep channels somewhere else (Bardou et al. 2003). The quantification of both, the rhythm and the magnitude of these changes, is very important to assess the variability of the bed's cross section and long profile. These parameters are indispensable for numerical modelling, as they should be considered as initial conditions. To monitor the channel evolution an Optech ILRIS 3D terrestrial laser scanner (LIDAR) was used. LIDAR permits to make a complete high precision 3D model of the channel and its surroundings by scanning it from different view points. The 3D data are treated and interpreted with the software Polyworks from Innovmetric Software Inc. Sequential 3D models allow for the determination of the variation in the bed's cross section and long profile. These data will afterwards be used to quantify the erosion and the deposition in the torrent reaches. To complete the chronological evolution of the landforms, precise digital terrain models, obtained by high resolution photogrammetry based on old aerial photographs, will be used. A 500 m long section of the Illgraben channel was scanned on 18th of August 2005 and on 7th of April 2006. These two data sets permit identifying the changes of the channel that occurred during the winter season. An upcoming scanning campaign in September 2006 will allow for the determination of the changes during this summer. Preliminary results show huge variations in the pathway of the Illgraben channel, as well as important vertical and lateral erosion of the river bed. Here we present the results of a river bank on the left (north-western) flank of the channel (Figure 1). For the August 2005 model the scans from 3 viewpoints were superposed, whereas the April 2006 3D image was obtained by combining 5 separate scans. The bank was eroded. The bank got eroded essentially on its left part (up to 6.3 m), where it is hit by the river and the debris flows (Figures 2 and 3). A debris cone has also formed (Figure 3), which suggests that a part of the bank erosion is due to shallow landslides. They probably occur when the river erosion creates an undercut slope. These geometrical data allow for the monitoring of the alluvial dynamics (i.e. aggradation and degradation) on different time scales and the influence of debris flows occurrence on these changes. Finally, the resistance against erosion of the bed's cross section and long profile will be analysed to assess the variability of these two key parameters. This information may then be used in debris flow simulation.

Weighted Shape-based Averaging with Neighborhood Prior Model for Multiple Atlas Fusion-based Medical Image Segmentation

In medical imaging, merging automated segmentations obtained from multiple atlases has become a standard practice for improving the accuracy. In this letter, we propose two new fusion methods: "Global Weighted Shape-Based Averaging" (GWSBA) and "Local Weighted Shape-Based Averaging" (LWSBA). These methods extend the well known Shape-Based Averaging (SBA) by additionally incorporating the similarity information between the reference (i.e., atlas) images and the target image to be segmented. We also propose a new spatially-varying similarity-weighted neighborhood prior model, and an edge-preserving smoothness term that can be used with many of the existing fusion methods. We first present our new Markov Random Field (MRF) based fusion framework that models the above mentioned information. The proposed methods are evaluated in the context of segmentation of lymph nodes in the head and neck 3D CT images, and they resulted in more accurate segmentations compared to the existing SBA.

Assessment of distribution and evolution of mechanical dyssynchrony in a porcine model of myocardial infarction by cardiovascular magnetic resonance.

BACKGROUND: We sought to investigate the relationship between infarct and dyssynchrony post- myocardial infarct (MI), in a porcine model. Mechanical dyssynchrony post-MI is associated with left ventricular (LV) remodeling and increased mortality. METHODS: Cine, gadolinium-contrast, and tagged cardiovascular magnetic resonance (CMR) were performed pre-MI, 9 ± 2 days (early post-MI), and 33 ± 10 days (late post-MI) post-MI in 6 pigs to characterize cardiac morphology, location and extent of MI, and regional mechanics. LV mechanics were assessed by circumferential strain (eC). Electro-anatomic mapping (EAM) was performed within 24 hrs of CMR and prior to sacrifice. RESULTS: Mean infarct size was 21 ± 4% of LV volume with evidence of post-MI remodeling. Global eC significantly decreased post MI (-27 ± 1.6% vs. -18 ± 2.5% (early) and -17 ± 2.7% (late), p < 0.0001) with no significant change in peri-MI and MI segments between early and late time-points. Time to peak strain (TTP) was significantly longer in MI, compared to normal and peri-MI segments, both early (440 ± 40 ms vs. 329 ± 40 ms and 332 ± 36 ms, respectively; p = 0.0002) and late post-MI (442 ± 63 ms vs. 321 ± 40 ms and 355 ± 61 ms, respectively; p = 0.012). The standard deviation of TTP in 16 segments (SD16) significantly increased post-MI: 28 ± 7 ms to 50 ± 10 ms (early, p = 0.012) to 54 ± 19 ms (late, p = 0.004), with no change between early and late post-MI time-points (p = 0.56). TTP was not related to reduction of segmental contractility. EAM revealed late electrical activation and greatly diminished conduction velocity in the infarct (5.7 ± 2.4 cm/s), when compared to peri-infarct (18.7 ± 10.3 cm/s) and remote myocardium (39 ± 20.5 cm/s). CONCLUSIONS: Mechanical dyssynchrony occurs early after MI and is the result of delayed electrical and mechanical activation in the infarct.

Anàlisi de l'aplicació GeoMedia 3D

Treball de final de carrera sobre els sistemes d'informació geogràfica on s'analitza l'eina GeoMedia 3D i la seva aplicació als SIG 3D. L'anàlisi està fet desenvolupant un model 3D de la ciutat de Tarragona.

A 2D/3D image analysis system to track fluorescently labeled structures in rod-shaped cells: application to measure spindle pole asymmetry during mitosis.

BACKGROUND: The yeast Schizosaccharomyces pombe is frequently used as a model for studying the cell cycle. The cells are rod-shaped and divide by medial fission. The process of cell division, or cytokinesis, is controlled by a network of signaling proteins called the Septation Initiation Network (SIN); SIN proteins associate with the SPBs during nuclear division (mitosis). Some SIN proteins associate with both SPBs early in mitosis, and then display strongly asymmetric signal intensity at the SPBs in late mitosis, just before cytokinesis. This asymmetry is thought to be important for correct regulation of SIN signaling, and coordination of cytokinesis and mitosis. In order to study the dynamics of organelles or large protein complexes such as the spindle pole body (SPB), which have been labeled with a fluorescent protein tag in living cells, a number of the image analysis problems must be solved; the cell outline must be detected automatically, and the position and signal intensity associated with the structures of interest within the cell must be determined. RESULTS: We present a new 2D and 3D image analysis system that permits versatile and robust analysis of motile, fluorescently labeled structures in rod-shaped cells. We have designed an image analysis system that we have implemented as a user-friendly software package allowing the fast and robust image-analysis of large numbers of rod-shaped cells. We have developed new robust algorithms, which we combined with existing methodologies to facilitate fast and accurate analysis. Our software permits the detection and segmentation of rod-shaped cells in either static or dynamic (i.e. time lapse) multi-channel images. It enables tracking of two structures (for example SPBs) in two different image channels. For 2D or 3D static images, the locations of the structures are identified, and then intensity values are extracted together with several quantitative parameters, such as length, width, cell orientation, background fluorescence and the distance between the structures of interest. Furthermore, two kinds of kymographs of the tracked structures can be established, one representing the migration with respect to their relative position, the other representing their individual trajectories inside the cell. This software package, called "RodCellJ", allowed us to analyze a large number of S. pombe cells to understand the rules that govern SIN protein asymmetry. CONCLUSIONS: "RodCell" is freely available to the community as a package of several ImageJ plugins to simultaneously analyze the behavior of a large number of rod-shaped cells in an extensive manner. The integration of different image-processing techniques in a single package, as well as the development of novel algorithms does not only allow to speed up the analysis with respect to the usage of existing tools, but also accounts for higher accuracy. Its utility was demonstrated on both 2D and 3D static and dynamic images to study the septation initiation network of the yeast Schizosaccharomyces pombe. More generally, it can be used in any kind of biological context where fluorescent-protein labeled structures need to be analyzed in rod-shaped cells. AVAILABILITY: RodCellJ is freely available under http://bigwww.epfl.ch/algorithms.html, (after acceptance of the publication).