Vessel centerline tracking and boundary sgementation in coronary MRA with minimal manual interaction.

Magnetic resonance angiography (MRA) provides a noninvasive means to detect the presence, location and severity of atherosclerosis throughout the vascular system. In such studies, and especially those in the coronary arteries, the vessel luminal area is typically measured at multiple cross-sectional locations along the course of the artery. The advent of fast volumetric imaging techniques covering proximal to mid segments of coronary arteries necessitates automatic analysis tools requiring minimal manual interactions to robustly measure cross-sectional area along the three-dimensional track of the arteries in under-sampled and non-isotropic datasets. In this work, we present a modular approach based on level set methods to track the vessel centerline, segment the vessel boundaries, and measure transversal area using two user-selected endpoints in each coronary of interest. Arterial area and vessel length are measured using our method and compared to the standard Soap-Bubble reformatting and analysis tool in in-vivo non-contrast enhanced coronary MRA images.

Self-Organising Maps and Correlation Analysis as a Tool to Explore Patterns in Excitation-Emission Matrix Data Sets and to Discriminate Dissolved Organic Matter Fluorescence Components

Dissolved organic matter (DOM) is a complex mixture of organic compounds, ubiquitous in marine and freshwater systems. Fluorescence spectroscopy, by means of Excitation-Emission Matrices (EEM), has become an indispensable tool to study DOM sources, transport and fate in aquatic ecosystems. However the statistical treatment of large and heterogeneous EEM data sets still represents an important challenge for biogeochemists. Recently, Self-Organising Maps (SOM) has been proposed as a tool to explore patterns in large EEM data sets. SOM is a pattern recognition method which clusterizes and reduces the dimensionality of input EEMs without relying on any assumption about the data structure. In this paper, we show how SOM, coupled with a correlation analysis of the component planes, can be used both to explore patterns among samples, as well as to identify individual fluorescence components. We analysed a large and heterogeneous EEM data set, including samples from a river catchment collected under a range of hydrological conditions, along a 60-km downstream gradient, and under the influence of different degrees of anthropogenic impact. According to our results, chemical industry effluents appeared to have unique and distinctive spectral characteristics. On the other hand, river samples collected under flash flood conditions showed homogeneous EEM shapes. The correlation analysis of the component planes suggested the presence of four fluorescence components, consistent with DOM components previously described in the literature. A remarkable strength of this methodology was that outlier samples appeared naturally integrated in the analysis. We conclude that SOM coupled with a correlation analysis procedure is a promising tool for studying large and heterogeneous EEM data sets.

Translation of polarity cues into asymmetric spindle positioning in "Caenorhabditis elegans"

Abstract: Asymmetric cell division is important to generate tissue diversity. The Caenorhabditis elegans embryo is well suited to study the mechanisms of asymmetric cell division. In wild type one-cell stage embryos, the spindle sets up along the anterior-posterior axis (AP). During anaphase, the spindle elongates. While the anterior spindle pole is relatively immobile, the posterior spindle pole moves towards the posterior cortex during anaphase leading to an asymmetric spindle position. As a result, the first cleavage gives rise to a large anterior blastomere and a smaller posterior one, which differs also in cell fate determinants. This posterior spindle displacement occurs in response to polarity cues set up along the AP axis by the PAR proteins and is due to imbalanced pulling forces acting on the two spindle poles, with net forces acting on the posterior spindle pole being more extensive than those at the anterior one. The project of my thesis was to characterize the involvement of two new components, gpr-1 and gpr-2, in spindle positioning. These genes encode essentially identical proteins containing a GoLoco motif characteristic of proteins interacting with α subunits of heterotrimeric G protein (Gα). In gpr-1/2(RNAi) embryos and in embryos lacking simultaneously two α subunits, goa-1 and gpa-16, (Ga(RNAi) embryos), there is a minimal posterior displacement of the spindle during anaphase, and the first division is equal. I found that the pulling forces acting on the two spindle poles is weak and equal in gpr-1/2(RNAi) and Gα (RNAi) embryos. I found that GPR-1/2 acts downstream of polarity cues for generation of pulling forces. Furthermore, I showed that GPR-1/2 distribution was enriched at the posterior cortex during metaphase whereas GOA-1 and GPA-16 were uniformly distributed at the cell cortex throughout the cell cycle. Gα subunits oscillate between GDP- and GTP-bound forms. Gα signaling is turned on by GDP/GTP exchange catalyzed by guanine nucleotide exchange factors (GEFs) and turned off by hydrolysis of GTP catalyzed by GTPase activating proteins (GAPs). A third class of proteins, the guanine dissociation inhibitors (GDIs), binds the GDP-bound form of Gα subunits and inhibits nucleotide exchange. I found that GPR-1/2 acts as a GDI for GOA-1. Taken together, my findings suggest a model in which differential activation of Gα subunits along the AP axis may translate into generation of differential pulling forces on the anterior and posterior spindle poles, and, thus, asymmetric cell division. Résumé L'embryon du nématode Caenorhabditis elegans est un modèle approprié pour étudier les mécanismes de la division asymétrique. Chez l'embryon précoce, le fuseau mitotique se forme le long de l'axe antéro-postérieur (A/P) et au centre de l'embryon, le pôle antérieur restant relativement immobile alors que le pôle postérieur du fuseau se déplace vers le cortex postérieur au cours de l'anaphase conduisant à une position excentrée du fuseau. 11 en résulte une première division qui génère un blastomère antérieur et postérieur de grande et petite taille respectivement et qui diffèrent en facteurs développementaux. Ce déplacement postérieur se produit en réponse de la polarité établie par la distribution polarisée des protéines PAR et est le résultat de la génération de forces inégales tirant sur les deux pôles du fuseau, les forces agissant sur le pôle postérieur du fuseau étant plus grandes. Le projet de ma thèse était d'identifier la fonction de deux nouveaux constituants, gpr-1 et gpr-2 dans le positionnement asymétrique du fuseau. Ces gènes codent essentiellement pour la même protéine qui contient un motif GoLoco, caractéristique des protéines interagissant avec la sous-unité alpha des protéines G hétérotrimériques. Chez l'embryon gpr-1/2(RNAi) et chez les embryons dépourvus d'activité de deux sous-unités alpha, goa-1 et gpa-16, (Gα(RNAi)), j'ai montré qu'il y avait un déplacement minimal du fuseau vers le pôle postérieur au cours de l'anaphase et la première division est symétrique en raison de forces faibles et égales agissant sur les deux pôles du fuseau. J'ai également montré que gpr-1/2 était requis en aval des signaux établissant la polarité pour générer les forces responsables du positionnement asymétrique du fuseau. De plus, j'ai montré que GPR-1/2 était enrichi au pôle postérieur lors de la métaphase alors que GOA-1 et GPA-16 étaient localisés de façon uniforme au cortex de l'embryon précoce. Gas oscillent entre une forme liée au GDP et une forme liée au GTP. La signalisation des Gas est activée par l'échange GDP/GTP qui est catalysé par des protéines GEFs. La signalisation des Gas est désactivée par l'hydrolyse du GTP qui est catalysée par des protéines GAPs. Une troisième classe de protéines, GDIs lie la forme GDP et inhibe l'échange de nucléotides. J'ai montré que GPR-1/2 agissait comme un GDI pour GOA-1. Mes résultats suggèrent un modèle dans lequel une activation différentielle des Gα le long de l'axe A/P pourrait générer des forces différentielles sur le pôle antérieur et postérieur du fuseau.

A Topological Approach to Fuzzy Sets

This thesis presents a topological approach to studying fuzzy setsby means of modifier operators. Modifier operators are mathematical models, e.g., for hedges, and we present briefly different approaches to studying modifier operators. We are interested in compositional modifier operators, modifiers for short, and these modifiers depend on binary relations. We show that if a modifier depends on a reflexive and transitive binary relation on U, then there exists a unique topology on U such that this modifier is the closure operator in that topology. Also, if U is finite then there exists a lattice isomorphism between the class of all reflexive and transitive relations and the class of all topologies on U. We define topological similarity relation "≈" between L-fuzzy sets in an universe U, and show that the class LU/ ≈ is isomorphic with the class of all topologies on U, if U is finite and L is suitable. We consider finite bitopological spaces as approximation spaces, and we show that lower and upper approximations can be computed by means of α-level sets also in the case of equivalence relations. This means that approximations in the sense of Rough Set Theory can be computed by means of α-level sets. Finally, we present and application to data analysis: we study an approach to detecting dependencies of attributes in data base-like systems, called information systems.

Distance Transforms on Gray-Level Surfaces

This thesis studies gray-level distance transforms, particularly the Distance Transform on Curved Space (DTOCS). The transform is produced by calculating distances on a gray-level surface. The DTOCS is improved by definingmore accurate local distances, and developing a faster transformation algorithm. The Optimal DTOCS enhances the locally Euclidean Weighted DTOCS (WDTOCS) with local distance coefficients, which minimize the maximum error from the Euclideandistance in the image plane, and produce more accurate global distance values.Convergence properties of the traditional mask operation, or sequential localtransformation, and the ordered propagation approach are analyzed, and compared to the new efficient priority pixel queue algorithm. The Route DTOCS algorithmdeveloped in this work can be used to find and visualize shortest routes between two points, or two point sets, along a varying height surface. In a digital image, there can be several paths sharing the same minimal length, and the Route DTOCS visualizes them all. A single optimal path can be extracted from the route set using a simple backtracking algorithm. A new extension of the priority pixel queue algorithm produces the nearest neighbor transform, or Voronoi or Dirichlet tessellation, simultaneously with the distance map. The transformation divides the image into regions so that each pixel belongs to the region surrounding the reference point, which is nearest according to the distance definition used. Applications and application ideas for the DTOCS and its extensions are presented, including obstacle avoidance, image compression and surface roughness evaluation.

A priori parameterisation of the CERES soil-crop models and tests against several European data sets

Mechanistic soil-crop models have become indispensable tools to investigate the effect of management practices on the productivity or environmental impacts of arable crops. Ideally these models may claim to be universally applicable because they simulate the major processes governing the fate of inputs such as fertiliser nitrogen or pesticides. However, because they deal with complex systems and uncertain phenomena, site-specific calibration is usually a prerequisite to ensure their predictions are realistic. This statement implies that some experimental knowledge on the system to be simulated should be available prior to any modelling attempt, and raises a tremendous limitation to practical applications of models. Because the demand for more general simulation results is high, modellers have nevertheless taken the bold step of extrapolating a model tested within a limited sample of real conditions to a much larger domain. While methodological questions are often disregarded in this extrapolation process, they are specifically addressed in this paper, and in particular the issue of models a priori parameterisation. We thus implemented and tested a standard procedure to parameterize the soil components of a modified version of the CERES models. The procedure converts routinely-available soil properties into functional characteristics by means of pedo-transfer functions. The resulting predictions of soil water and nitrogen dynamics, as well as crop biomass, nitrogen content and leaf area index were compared to observations from trials conducted in five locations across Europe (southern Italy, northern Spain, northern France and northern Germany). In three cases, the model’s performance was judged acceptable when compared to experimental errors on the measurements, based on a test of the model’s root mean squared error (RMSE). Significant deviations between observations and model outputs were however noted in all sites, and could be ascribed to various model routines. In decreasing importance, these were: water balance, the turnover of soil organic matter, and crop N uptake. A better match to field observations could therefore be achieved by visually adjusting related parameters, such as field-capacity water content or the size of soil microbial biomass. As a result, model predictions fell within the measurement errors in all sites for most variables, and the model’s RMSE was within the range of published values for similar tests. We conclude that the proposed a priori method yields acceptable simulations with only a 50% probability, a figure which may be greatly increased through a posteriori calibration. Modellers should thus exercise caution when extrapolating their models to a large sample of pedo-climatic conditions for which they have only limited information.

Transanal minimal invasive surgery for rectal lesions: should the defect be closed?

AIM: Transanal minimal invasive surgery (TAMIS) of rectal lesions is increasingly being used, but the technique is not yet standardized. The aims of this study were to evaluate peri-operative complications and long-term functional outcome of the technique and to analyse whether or not the rectal defect needs to be closed. METHOD: Consecutive patients undergoing TAMIS using the SILS port (Covidien) and standard laparoscopic instruments were studied. RESULTS: Seventy-five patients (68% male) of mean age 67 (± 15) years underwent single-port transanal surgery at three different centres for 37 benign lesions and 38 low-risk cancers located at a mean of 6.4 ± 2.3 cm from the anal verge. The median operating time was 77 (25-245) min including a median time for resection of 36 (15-75) min and for closure of the rectal defect of 38 (9-105) min. The defect was closed in 53% using interrupted (75%) or a running suture (25%). Intra-operative complications occurred in six (8%) patients and postoperative morbidity was 19% with only one patient requiring reoperation for Grade IIIb local infection. There was no difference in the incidence of complications whether the rectal defect was closed or left open. Patients were discharged after 3.4 (1-21) days. At a median follow-up of 12.8 (2-29) months, the continence was normal (Vaizey score of 1.5; 0-16). CONCLUSION: Transanal rectal resection can be safely and efficiently performed by means of a SILS port and standard laparoscopic instruments. The rectal defect may be left open and at 1 year continence is not compromised.

Minimal requirements for parameters to be used for the development of predictive models for microbial source tracking. The example of the pair somatic coliphages and phages infecting Bacteroides

Several methods and approaches for measuring parameters to determine fecal sources of pollution in water have been developed in recent years. No single microbial or chemical parameter has proved sufficient to determine the source of fecal pollution. Combinations of parameters involving at least one discriminating indicator and one universal fecal indicator offer the most promising solutions for qualitative and quantitative analyses. The universal (nondiscriminating) fecal indicator provides quantitative information regarding the fecal load. The discriminating indicator contributes to the identification of a specific source. The relative values of the parameters derived from both kinds of indicators could provide information regarding the contribution to the total fecal load from each origin. It is also essential that both parameters characteristically persist in the environment for similar periods. Numerical analysis, such as inductive learning methods, could be used to select the most suitable and the lowest number of parameters to develop predictive models. These combinations of parameters provide information on factors affecting the models, such as dilution, specific types of animal source, persistence of microbial tracers, and complex mixtures from different sources. The combined use of the enumeration of somatic coliphages and the enumeration of Bacteroides-phages using different host specific strains (one from humans and another from pigs), both selected using the suggested approach, provides a feasible model for quantitative and qualitative analyses of fecal source identification.

Self-Organising Maps and Correlation Analysis as a Tool to Explore Patterns in Excitation-Emission Matrix Data Sets and to Discriminate Dissolved Organic Matter Fluorescence Components

Dissolved organic matter (DOM) is a complex mixture of organic compounds, ubiquitous in marine and freshwater systems. Fluorescence spectroscopy, by means of Excitation-Emission Matrices (EEM), has become an indispensable tool to study DOM sources, transport and fate in aquatic ecosystems. However the statistical treatment of large and heterogeneous EEM data sets still represents an important challenge for biogeochemists. Recently, Self-Organising Maps (SOM) has been proposed as a tool to explore patterns in large EEM data sets. SOM is a pattern recognition method which clusterizes and reduces the dimensionality of input EEMs without relying on any assumption about the data structure. In this paper, we show how SOM, coupled with a correlation analysis of the component planes, can be used both to explore patterns among samples, as well as to identify individual fluorescence components. We analysed a large and heterogeneous EEM data set, including samples from a river catchment collected under a range of hydrological conditions, along a 60-km downstream gradient, and under the influence of different degrees of anthropogenic impact. According to our results, chemical industry effluents appeared to have unique and distinctive spectral characteristics. On the other hand, river samples collected under flash flood conditions showed homogeneous EEM shapes. The correlation analysis of the component planes suggested the presence of four fluorescence components, consistent with DOM components previously described in the literature. A remarkable strength of this methodology was that outlier samples appeared naturally integrated in the analysis. We conclude that SOM coupled with a correlation analysis procedure is a promising tool for studying large and heterogeneous EEM data sets.

Sierpinski curve Julia sets for quadratic rational maps

We investigate under which dynamical conditions the Julia set of a quadratic rational map is a Sierpiński curve.

Beurling-Landau densities of weighted Fekete sets and correlation kernel estimates

Let $Q$ be a suitable real function on $C$. An $n$-Fekete set corresponding to $Q$ is a subset ${Z_{n1}},\dotsb, Z_{nn}}$ of $C$ which maximizes the expression $\Pi^n_i_{sets are, in a sense, maximally spread out with respect to the equilibrium measure. In general, our results apply only to a part of the Fekete set, which is at a certain distance away from the boundary of the droplet. However, for the potential $Q=|Z|^2$ we obtain results which hold globally, and we conjecture that such global results are true for a wide range of potentials.

A minimal set of coordinates for describing humanoid shoulder motion

The kinematics of the anatomical shoulder are analysed and modelled as a parallel mechanism similar to a Stewart platform. A new method is proposed to describe the shoulder kinematics with minimal coordinates and solve the indeterminacy. The minimal coordinates are defined from bony landmarks and the scapulothoracic kinematic constraints. Independent from one another, they uniquely characterise the shoulder motion. A humanoid mechanism is then proposed with identical kinematic properties. It is then shown how minimal coordinates can be obtained for this mechanism and how the coordinates simplify both the motion-planning task and trajectory-tracking control. Lastly, the coordinates are also shown to have an application in the field of biomechanics where they can be used to model the scapulohumeral rhythm.