Automated 3D Lumbar Intervertebral Disc Segmentation from MRI Data Sets

This paper proposed an automated 3D lumbar intervertebral disc (IVD) segmentation strategy from MRI data. Starting from two user supplied landmarks, the geometrical parameters of all lumbar vertebral bodies and intervertebral discs are automatically extracted from a mid-sagittal slice using a graphical model based approach. After that, a three-dimensional (3D) variable-radius soft tube model of the lumbar spine column is built to guide the 3D disc segmentation. The disc segmentation is achieved as a multi-kernel diffeomorphic registration between a 3D template of the disc and the observed MRI data. Experiments on 15 patient data sets showed the robustness and the accuracy of the proposed algorithm.

SALSA: a tool to estimate the stray light contamination for low-Earth orbit observatories

Stray light contamination reduces considerably the precision of photometric of faint stars for low altitude spaceborne observatories. When measuring faint objects, the necessity of coping with stray light contamination arises in order to avoid systematic impacts on low signal-to-noise images. Stray light contamination can be represented by a flat offset in CCD data. Mitigation techniques begin by a comprehensive study during the design phase, followed by the use of target pointing optimisation and post-processing methods. We present a code that aims at simulating the stray-light contamination in low-Earth orbit coming from reflexion of solar light by the Earth. StrAy Light SimulAtor (SALSA) is a tool intended to be used at an early stage as a tool to evaluate the effective visible region in the sky and, therefore to optimise the observation sequence. SALSA can compute Earth stray light contamination for significant periods of time allowing missionwide parameters to be optimised (e.g. impose constraints on the point source transmission function (PST) and/or on the altitude of the satellite). It can also be used to study the behaviour of the stray light at different seasons or latitudes. Given the position of the satellite with respect to the Earth and the Sun, SALSA computes the stray light at the entrance of the telescope following a geometrical technique. After characterising the illuminated region of the Earth, the portion of illuminated Earth that affects the satellite is calculated. Then, the flux of reflected solar photons is evaluated at the entrance of the telescope. Using the PST of the instrument, the final stray light contamination at the detector is calculated. The analysis tools include time series analysis of the contamination, evaluation of the sky coverage and an objects visibility predictor. Effects of the South Atlantic Anomaly and of any shutdown periods of the instrument can be added. Several designs or mission concepts can be easily tested and compared. The code is not thought as a stand-alone mission designer. Its mandatory inputs are a time series describing the trajectory of the satellite and the characteristics of the instrument. This software suite has been applied to the design and analysis of CHEOPS (CHaracterizing ExOPlanet Satellite). This mission requires very high precision photometry to detect very shallow transits of exoplanets. Different altitudes and characteristics of the detector have been studied in order to find the best parameters, that reduce the effect of contamination. © (2014) COPYRIGHT Society of Photo-Optical Instrumentation Engineers (SPIE). Downloading of the abstract is permitted for personal use only.

A steerable UV laser system for the calibration of liquid argon time projection chambers

A number of liquid argon time projection chambers (LAr TPCs) are being built or are proposed for neutrino experiments on long- and short baseline beams. For these detectors, a distortion in the drift field due to geometrical or physics reasons can affect the reconstruction of the events. Depending on the TPC geometry and electric drift field intensity, this distortion could be of the same magnitude as the drift field itself. Recently, we presented a method to calibrate the drift field and correct for these possible distortions. While straight cosmic ray muon tracks could be used for calibration, multiple coulomb scattering and momentum uncertainties allow only a limited resolution. A UV laser instead can create straight ionization tracks in liquid argon, and allows one to map the drift field along different paths in the TPC inner volume. Here we present a UV laser feed-through design with a steerable UV mirror immersed in liquid argon that can point the laser beam at many locations through the TPC. The straight ionization paths are sensitive to drift field distortions, a fit of these distortion to the linear optical path allows to extract the drift field, by using these laser tracks along the whole TPC volume one can obtain a 3D drift field map. The UV laser feed-through assembly is a prototype of the system that will be used for the MicroBooNE experiment at the Fermi National Accelerator Laboratory (FNAL).

MorphoGraphX: A platform for quantifying morphogenesis in 4D

Morphogenesis emerges from complex multiscale interactions between genetic and mechanical processes. To understand these processes, the evolution of cell shape, proliferation and gene expression must be quantified. This quantification is usually performed either in full 3D, which is computationally expensive and technically challenging, or on 2D planar projections, which introduces geometrical artifacts on highly curved organs. Here we present MorphoGraphX (www.MorphoGraphX.org), a software that bridges this gap by working directly with curved surface images extracted from 3D data. In addition to traditional 3D image analysis, we have developed algorithms to operate on curved surfaces, such as cell segmentation, lineage tracking and fluorescence signal quantification. The software’s modular design makes it easy to include existing libraries, or to implement new algorithms. Cell geometries extracted with MorphoGraphX can be exported and used as templates for simulation models, providing a powerful platform to investigate the interactions between shape, genes and growth.DOI: http://dx.doi.org/10.7554/eLife.05864.001Author keywordsResearch organism

Enabling the Next Generation of Spaceborne Quadrupole Mass Spectrometers

The quadrupole mass spectrometer (QMS) has over 30 years of spaceflight heritage in making important neutral gas and low energy ion observations. Given their geometrical constraints, these instruments are currently operated at the extreme limit of their capabilities. However, a technique called higher order auxiliary excitation provides a set of novel, robust, electronics-based solutions for improving the performance of these sensors. By driving the quadrupole rods with an additional frequency nearly twice that of the normal RF operating frequency, substantially increased abundance sensitivity, maximum attainable mass resolution, and peak stability can be achieved through operation of voltage scan lines through the center of formed upper stability islands. Such improvements are modeled using numerical simulations of ion trajectories in a quadrupole field with and without applied higher order auxiliary excitation. When compared to a traditional QMS with a mass range up to 500Da, sensors can be designed with the same precision electronics to have expected mass ranges beyond 1500Da with a power increase of less than twice that of its heritage implementations.

Streak detection algorithm for space debris detection on optical images.

Any image processing object detection algorithm somehow tries to integrate the object light (Recognition Step) and applies statistical criteria to distinguish objects of interest from other objects or from pure background (Decision Step). There are various possibilities how these two basic steps can be realized, as can be seen in the different proposed detection methods in the literature. An ideal detection algorithm should provide high recognition sensitiv ity with high decision accuracy and require a reasonable computation effort . In reality, a gain in sensitivity is usually only possible with a loss in decision accuracy and with a higher computational effort. So, automatic detection of faint streaks is still a challenge. This paper presents a detection algorithm using spatial filters simulating the geometrical form of possible streaks on a CCD image. This is realized by image convolution. The goal of this method is to generate a more or less perfect match between a streak and a filter by varying the length and orientation of the filters. The convolution answers are accepted or rejected according to an overall threshold given by the ackground statistics. This approach yields as a first result a huge amount of accepted answers due to filters partially covering streaks or remaining stars. To avoid this, a set of additional acceptance criteria has been included in the detection method. All criteria parameters are justified by background and streak statistics and they affect the detection sensitivity only marginally. Tests on images containing simulated streaks and on real images containing satellite streaks show a very promising sensitivity, reliability and running speed for this detection method. Since all method parameters are based on statistics, the true alarm, as well as the false alarm probability, are well controllable. Moreover, the proposed method does not pose any extraordinary demands on the computer hardware and on the image acquisition process.

Automated 3D Lumbar Intervertebral Disc Segmentation from MRI Data Sets

This paper proposed an automated three-dimensional (3D) lumbar intervertebral disc (IVD) segmentation strategy from Magnetic Resonance Imaging (MRI) data. Starting from two user supplied landmarks, the geometrical parameters of all lumbar vertebral bodies and intervertebral discs are automatically extracted from a mid-sagittal slice using a graphical model based template matching approach. Based on the estimated two-dimensional (2D) geometrical parameters, a 3D variable-radius soft tube model of the lumbar spine column is built by model fitting to the 3D data volume. Taking the geometrical information from the 3D lumbar spine column as constraints and segmentation initialization, the disc segmentation is achieved by a multi-kernel diffeomorphic registration between a 3D template of the disc and the observed MRI data. Experiments on 15 patient data sets showed the robustness and the accuracy of the proposed algorithm.

Quantitative multi-slice computed tomography assessment of the mitral valvular complex for transcatheter mitral valve interventions part 1: systematic measurement methodology and inter-observer variability

AIMS Transcatheter mitral valve replacement (TMVR) is an emerging technology with the potential to treat patients with severe mitral regurgitation at excessive risk for surgical mitral valve surgery. Multimodal imaging of the mitral valvular complex and surrounding structures will be an important component for patient selection for TMVR. Our aim was to describe and evaluate a systematic multi-slice computed tomography (MSCT) image analysis methodology that provides measurements relevant for transcatheter mitral valve replacement. METHODS AND RESULTS A systematic step-by-step measurement methodology is described for structures of the mitral valvular complex including: the mitral valve annulus, left ventricle, left atrium, papillary muscles and left ventricular outflow tract. To evaluate reproducibility, two observers applied this methodology to a retrospective series of 49 cardiac MSCT scans in patients with heart failure and significant mitral regurgitation. For each of 25 geometrical metrics, we evaluated inter-observer difference and intra-class correlation. The inter-observer difference was below 10% and the intra-class correlation was above 0.81 for measurements of critical importance in the sizing of TMVR devices: the mitral valve annulus diameters, area, perimeter, the inter-trigone distance, and the aorto-mitral angle. CONCLUSIONS MSCT can provide measurements that are important for patient selection and sizing of TMVR devices. These measurements have excellent inter-observer reproducibility in patients with functional mitral regurgitation.

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.

Dimensions of projections of sets on Riemannian surfaces of constant curvature

We apply the theory of Peres and Schlag to obtain generic lower bounds for Hausdorff dimension of images of sets by orthogonal projections on simply connected two-dimensional Riemannian manifolds of constant curvature. As a conclusion we obtain appropriate versions of Marstrand's theorem, Kaufman's theorem, and Falconer's theorem in the above geometrical settings.

Large-scale dust jets in the coma of 67P/Churyumov-Gerasimenko as seen by the OSIRIS instrument onboard Rosetta

Context. During the most recent perihelion passage in 2009 of comet 67P/Churyumov-Gerasimenko (67P), ground-based observations showed an anisotropic dust coma where jet-like features were detected at similar to 1.3 AU from the Sun. The current perihelion passage is exceptional as the Rosetta spacecraft is monitoring the nucleus activity since March 2014, when a clear dust coma was already surrounding the nucleus at 4.3 AU from the Sun. Subsequently, the OSIRIS camera also witnessed an outburst in activity between April 27 and 30, and since mid-July, the dust coma at rh similar to 3.7-3.6 AU preperihelion is clearly non-isotropic, pointing to the existence of dust jet-like features. Aims. We aim to ascertain on the nucleus surface the origin of the dust jet-like features detected as early as in mid-July 2014. This will help to establish how the localized comet nucleus activity compares with that seen in previous apparitions and will also help following its evolution as the comet approaches its perihelion, at which phase most of the jets were detected from ground-based observations. Determining these areas also allows locating them in regions on the nucleus with spectroscopic or geomorphological distinct characteristics. Methods. Three series of dust images of comet 67P obtained with the Wide Angle Camera (WAC) of the OSIRIS instrument onboard the Rosetta spacecraft were processed with different enhancement techniques. This was made to clearly show the existence of jet-like features in the dust coma, whose appearance toward the observer changed as a result of the rotation of the comet nucleus and of the changing observing geometry from the spacecraft. The position angles of these features in the coma together with information on the observing geometry, nucleus shape, and rotation, allowed us to determine the most likely locations on the nucleus surface where the jets originate from. Results. Geometrical tracing of jet sources indicates that the activity of the nucleus of 67P gave rise during July and August 2014 to large-scale jet-like features from the Hapi, Hathor, Anuket, and Aten regions, confirming that active regions may be present on the nucleus localized at 60. northern latitude as deduced from previous comet apparitions. There are also hints that large-scale jets observed from the ground are possibly composed, at their place of origin on the nucleus surface, of numerous small-scale features.

NONLINEAR TIME SERIES/SYSTEM ANALYSIS (STATE-SPACE, DYNAMICS, INTERVENTION, RESILIENCE, KALMAN)

A discussion of nonlinear dynamics, demonstrated by the familiar automobile, is followed by the development of a systematic method of analysis of a possibly nonlinear time series using difference equations in the general state-space format. This format allows recursive state-dependent parameter estimation after each observation thereby revealing the dynamics inherent in the system in combination with random external perturbations.^ The one-step ahead prediction errors at each time period, transformed to have constant variance, and the estimated parametric sequences provide the information to (1) formally test whether time series observations y(,t) are some linear function of random errors (ELEM)(,s), for some t and s, or whether the series would more appropriately be described by a nonlinear model such as bilinear, exponential, threshold, etc., (2) formally test whether a statistically significant change has occurred in structure/level either historically or as it occurs, (3) forecast nonlinear system with a new and innovative (but very old numerical) technique utilizing rational functions to extrapolate individual parameters as smooth functions of time which are then combined to obtain the forecast of y and (4) suggest a measure of resilience, i.e. how much perturbation a structure/level can tolerate, whether internal or external to the system, and remain statistically unchanged. Although similar to one-step control, this provides a less rigid way to think about changes affecting social systems.^ Applications consisting of the analysis of some familiar and some simulated series demonstrate the procedure. Empirical results suggest that this state-space or modified augmented Kalman filter may provide interesting ways to identify particular kinds of nonlinearities as they occur in structural change via the state trajectory.^ A computational flow-chart detailing computations and software input and output is provided in the body of the text. IBM Advanced BASIC program listings to accomplish most of the analysis are provided in the appendix. ^

Ordenaciones de material genético a partir de información multidimensional

Nuevas biotecnologías permiten obtener información para caracterizar materiales genéticos a partir de múltiples marcadores, ya sean éstos moleculares y/o morfológicos. La ordenación del material genético a través de la exploración de patrones de variabilidad multidimensionales se aborda mediante diversas técnicas de análisis multivariado. Las técnicas multivariadas de reducción de dimensión (TRD) y la representación gráfica de las mismas cobran sustancial importancia en la visualización de datos multivariados en espacios de baja dimensión ya que facilitan la interpretación de interrelaciones entre las variables (marcadores) y entre los casos u observaciones bajo análisis. Tanto el Análisis de Componentes Principales, como el Análisis de Coordenadas Principales y el Análisis de Procrustes Generalizado son TRD aplicables a datos provenientes de marcadores moleculares y/o morfológicos. Los Árboles de Mínimo Recorrido y los biplots constituyen técnicas para lograr representaciones geométricas de resultados provenientes de TRD. En este trabajo se describen estas técnicas multivariadas y se ilustran sus aplicaciones sobre dos conjuntos de datos, moleculares y morfológicos, usados para caracterizar material genético fúngico.

Paleomagnetic of igneous rocks of ODP Leg 114

A paleomagnetic study was made of 12 samples of trachytic basalt from the base of ODP Hole 698A on the Northeast Georgia Rise (southwest Atlantic) and four samples of andesitic basalt and nine samples of volcanic breccia from the base of ODP Hole 703A on the Meteor Rise (southeast Atlantic). The magnetic intensities of the Hole 703A samples are anomalously low, possibly reflecting alteration effects. The mean magnetic intensity of the Hole 698A samples is high, and compatible with the model of Bleil and Petersen (1983) for the variation of magnetic intensity with age in oceanic basalts, involving progressive low-temperature oxidation of titanomagnetite to titanomaghemite for some 20 m.y. followed by inversion to intergrowths of magnetite and other Fe-Ti oxides during the subsequent 100 m.y. These results support the interpretation of the Hole 698A basalts as true oceanic basement of Late Cretaceous age rather than a younger intrusion. Well-defined stable components of magnetization were identified from AF and thermal demagnetization of the Hole 698A basalts, and less well-defined components were identified for the Hole 703A samples. Studies of the magnetic homogeneity of the Hole 698A basalts, involving harmonic analysis of the spinner magnetometer output, indicate the presence of an unevenly distributed low-coercivity component superimposed on the more homogeneous high-coercivity characteristic magnetization. The former component is believed to reside in irregularly distributed multidomain magnetite grains formed along cracks within the basalt, whilst the latter resides in more uniformly distributed finer magnetic grains. The inclination values for the high-coercivity magnetization of five Hole 698A basalt samples form an internally consistent set with a mean value of 59° ± 5°. The corresponding Late Cretaceous paleolatitude of 40° ± 5° is shallower than expected for this site but is broadly compatible with models for the opening of the South Atlantic involving pivoting of South America away from Africa since the Early Cretaceous. The polarity of the stable characteristic magnetization of the Site 698 basalts is normal. This is consistent with their emplacement during the long Campanian to Maestrichtian normal polarity Chron C33N.

Geochemistry of dark-light rhythm sediments of the Japan Sea

Ocean Drilling Program Legs 127 and 128 in the Japan Sea have revealed the existence of numerous dark-light rhythms of remarkable consistency in sediments of late Miocene, latest Pliocene, and especially Pleistocene age. Light-colored units within these rhythms are massive or bioturbated, consist of diatomaceous clays, silty clays, or nannofossil-rich clays, and are generally poor in organic matter. Dark-colored units are homogeneous, laminated, or thinly bedded and include substantial amounts of biogenic material such as well-preserved diatoms, planktonic foraminifers, calcareous nannofossils, and organic matter (maximum 7.4 wt%). The dark-light rhythms show a similar geometrical pattern on three different scales: First-order rhythms consist of a cluster dominated by dark-colored units followed by a cluster dominated by light-colored units (3-5 m). Spectral analysis of a gray-value time series suggests that the frequencies of the first-order rhythms in sediments of latest Pliocene and Pleistocene age correlate to the obliquity and the eccentricity cycles. The second-order dark-light rhythms include a light and a dark-colored unit (10-160 cm). They were formed in time spans of several hundred to several ten thousand years, with variance centering around 10,500 yr. This frequency may correspond to half the precessional cycle. Third-order rhythms appear as laminated or thinly bedded dark-light couplets (2-15 mm) within the dark-colored units of the second-order rhythms and may represent annual frequencies. In interpreting the rhythms, we have to take into account that (1) the occurrence of the first- and second-order rhythms is not necessarily restricted to glacial or interglacial periods as is shown by preliminary stable-isotope analysis and comparison with the published d18O record; (2) they appear to be Milankovitch-controlled; and (3) a significant number of the rhythms are sharply bounded. The origin of the dark-light rhythms is probably related to variations in monsoonal activity in the Japan Sea, which show annual frequencies, but also operates in phase with the orbital cycles.