Fully automatic localization and segmentation of 3D vertebral bodies from CT/MR images via a learning-based method

In this paper, we address the problems of fully automatic localization and segmentation of 3D vertebral bodies from CT/MR images. We propose a learning-based, unified random forest regression and classification framework to tackle these two problems. More specifically, in the first stage, the localization of 3D vertebral bodies is solved with random forest regression where we aggregate the votes from a set of randomly sampled image patches to get a probability map of the center of a target vertebral body in a given image. The resultant probability map is then further regularized by Hidden Markov Model (HMM) to eliminate potential ambiguity caused by the neighboring vertebral bodies. The output from the first stage allows us to define a region of interest (ROI) for the segmentation step, where we use random forest classification to estimate the likelihood of a voxel in the ROI being foreground or background. The estimated likelihood is combined with the prior probability, which is learned from a set of training data, to get the posterior probability of the voxel. The segmentation of the target vertebral body is then done by a binary thresholding of the estimated probability. We evaluated the present approach on two openly available datasets: 1) 3D T2-weighted spine MR images from 23 patients and 2) 3D spine CT images from 10 patients. Taking manual segmentation as the ground truth (each MR image contains at least 7 vertebral bodies from T11 to L5 and each CT image contains 5 vertebral bodies from L1 to L5), we evaluated the present approach with leave-one-out experiments. Specifically, for the T2-weighted MR images, we achieved for localization a mean error of 1.6 mm, and for segmentation a mean Dice metric of 88.7% and a mean surface distance of 1.5 mm, respectively. For the CT images we achieved for localization a mean error of 1.9 mm, and for segmentation a mean Dice metric of 91.0% and a mean surface distance of 0.9 mm, respectively.

Follicular diameter range based on morphological features in Synbranchus marmoratus (Bloch, 1795) (Teleostei, Synbranchiformes, Synbranchidae) from the South-central region of Brazil

The morphological characteristics of Synbranchus marmoratus female germ cells in various development stages were described in details; then measurements of ovarian follicle diameters were taken from primary and secondary growth as during these development stages the oocyte size varied considerably along the fish growth. The results were correlated to total fish length, using the individuals division in six size classes. It was possible to group oocytes by stages according to histological characteristics but not according to morphometric diameter, as there was a wide variation in diameter in each stage and overlap between different maturation stages. These data make available new information on the reproductive biology of Synbranchidae. (c) 2004 Published by Elsevier Ltd.

Morphometry of the crab Hexapanopeus schmitti (Decapoda : Xanthoidea) on the northern coast of the state of São Paulo, Brazil

A morphometric study of the xanthoid crab Hexapanopeus schmitti was carried out, using the allometric method. Samples were taken monthly for two years (1998-1999) in the Ubatuba region, northern coast of São Paulo, Brazil. Sex and size were assessed for each specimen, and all crabs were measured to obtain their carapace width (CW) and length (CL), abdomen width (AW) of females, major cheliped propodus length and height (PL and PH), and gonopod length (GL) of males. A total of 301 crabs were analyzed, 209 males and 92 females. The CWs of the crabs ranged from 2.5 to 9.8 mm for males and from 2.8 to 9.4 mm for females. The relative growth equation (y=ax(b)) based on the relationship between GL and CW suggested that males reach their morphological sexual maturity near 6.1 mm CW. In females, the estimated size at 50 % maturity was 4.8 mm CW, based on the relationship of AW vs. CW. Males reach larger sizes than females, which probably favors their ability to guard the females during courtship. In approximately 83 % of the crabs (n= 371), disregarding sex, the right cheliped was larger.

Morphometry and histochemistry of the Rectus abdominis muscle fibers of tufted capuchin monkeys (Cebus apella Linnaeus, 1758)

The phylogenetic proximity of primates to humans, along with their behavioral, biochemical, and anatomical similarities, make such animals more interesting experimental models for biomedical researches, as compared to classical laboratory animals. Another aspect that has called the attention of researchers is the differentiated quadrupedalism present in some primates. The tufted capuchin monkey uses the ground and tree branches as its support for locomotion, showing various postures while performing this task. On the basis of this information, we have decided to study the rectus abdominis muscle of the tufted capuchin monkey, with the following goals: the frequency and area of fiber types; its possible compartmentalization; and identify if this muscle is better adapted to phasic or postural activities. To do this, samples were removed from 4 regions of the rectus abdominis muscle of 6 adult male tufted capuchin monkeys, and were submitted to reaction with m-ATPase, (with alkaline and acid pre-incubation), NADH, and H.E.. Results showed: a statistically significant difference (P<0.05) for both frequency and area, between fiber types FG and FOG and FG and SO, but did not show a statistically significant difference between fibers FOG and SO, in all studied regions; similarity in frequency and area of a same fiber type (FG, FOG, and SO) among the studied regions. Based on these data, it was concluded that: the rectus abdominis muscle of the tufted capuchin monkey does not show fiber compartmentalization, since the distribution and size patterns of the different fiber types are similar in the studied regions; there is a predominance of fast twitch fibers (FG + FOG) over slow twitch fibers (SO), for frequency and area, which characterizes the muscle as being more dedicated to phasic than to postural activities. © 2006 Sociedad Chilena de Anatom.

Differentiation of deer species of the genus Mazama by track morphometry

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Inter- and intra-specific diversity of Cuban Pinguicula (Lentibulariaceae) based on morphometric analyses and its relation with geographical distribution

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Objective CT-Based Quantification of Lung Sequelae in Treated Patients With Paracoccidioidomycosis

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Multicentre study for a robust protocol in single-voxel spectroscopy: quantification of MRS signals by time-domain fitting algorithms

Magnetic Resonance Spectroscopy (MRS) is an advanced clinical and research application which guarantees a specific biochemical and metabolic characterization of tissues by the detection and quantification of key metabolites for diagnosis and disease staging. The "Associazione Italiana di Fisica Medica (AIFM)" has promoted the activity of the "Interconfronto di spettroscopia in RM" working group. The purpose of the study is to compare and analyze results obtained by perfoming MRS on scanners of different manufacturing in order to compile a robust protocol for spectroscopic examinations in clinical routines. This thesis takes part into this project by using the GE Signa HDxt 1.5 T at the Pavillion no. 11 of the S.Orsola-Malpighi hospital in Bologna. The spectral analyses have been performed with the jMRUI package, which includes a wide range of preprocessing and quantification algorithms for signal analysis in the time domain. After the quality assurance on the scanner with standard and innovative methods, both spectra with and without suppression of the water peak have been acquired on the GE test phantom. The comparison of the ratios of the metabolite amplitudes over Creatine computed by the workstation software, which works on the frequencies, and jMRUI shows good agreement, suggesting that quantifications in both domains may lead to consistent results. The characterization of an in-house phantom provided by the working group has achieved its goal of assessing the solution content and the metabolite concentrations with good accuracy. The goodness of the experimental procedure and data analysis has been demonstrated by the correct estimation of the T2 of water, the observed biexponential relaxation curve of Creatine and the correct TE value at which the modulation by J coupling causes the Lactate doublet to be inverted in the spectrum. The work of this thesis has demonstrated that it is possible to perform measurements and establish protocols for data analysis, based on the physical principles of NMR, which are able to provide robust values for the spectral parameters of clinical use.

How much is there really? Why stereology is essential in lung morphometry

Quantitative data on lung structure are essential to set up structure-function models for assessing the functional performance of the lung or to make statistically valid comparisons in experimental morphology, physiology, or pathology. The methods of choice for microscopy-based lung morphometry are those of stereology, the science of quantitative characterization of irregular three-dimensional objects on the basis of measurements made on two-dimensional sections. From a practical perspective, stereology is an assumption-free set of methods of unbiased sampling with geometric probes, based on a solid mathematical foundation. Here, we discuss the pitfalls of lung morphometry and present solutions, from specimen preparation to the sampling scheme in multiple stages, for obtaining unbiased estimates of morphometric parameters such as volumes, surfaces, lengths, and numbers. This is demonstrated on various examples. Stereological methods are accurate, efficient, simple, and transparent; the precision of the estimates depends on the size and distribution of the sample. For obtaining quantitative data on lung structure at all microscopic levels, state-of-the-art stereology is the gold standard.

A nonlinear QCT-based finite element model validation study for the human femur tested in two configurations in vitro

Purpose Femoral fracture is a common medical problem in osteoporotic individuals. Bone mineral density (BMD) is the gold standard measure to evaluate fracture risk in vivo. Quantitative computed tomography (QCT)-based homogenized voxel finite element (hvFE) models have been proved to be more accurate predictors of femoral strength than BMD by adding geometrical and material properties. The aim of this study was to evaluate the ability of hvFE models in predicting femoral stiffness, strength and failure location for a large number of pairs of human femora tested in two different loading scenarios. Methods Thirty-six pairs of femora were scanned with QCT and total proximal BMD and BMC were evaluated. For each pair, one femur was positioned in one-legged stance configuration (STANCE) and the other in a sideways configuration (SIDE). Nonlinear hvFE models were generated from QCT images by reproducing the same loading configurations imposed in the experiments. For experiments and models, the structural properties (stiffness and ultimate load), the failure location and the motion of the femoral head were computed and compared. Results In both configurations, hvFE models predicted both stiffness (R2=0.82 for STANCE and R2=0.74 for SIDE) and femoral ultimate load (R2=0.80 for STANCE and R2=0.85 for SIDE) better than BMD and BMC. Moreover, the models predicted qualitatively well the failure location (66% of cases) and the motion of the femoral head. Conclusions The subject specific QCT-based nonlinear hvFE model cannot only predict femoral apparent mechanical properties better than densitometric measures, but can additionally provide useful qualitative information about failure location.

3D morphometry using automated aortic segmentation in native MR angiography: an alternative to contrast enhanced MRA?

INTRODUCTION Native-MR angiography (N-MRA) is considered an imaging alternative to contrast enhanced MR angiography (CE-MRA) for patients with renal insufficiency. Lower intraluminal contrast in N-MRA often leads to failure of the segmentation process in commercial algorithms. This study introduces an in-house 3D model-based segmentation approach used to compare both sequences by automatic 3D lumen segmentation, allowing for evaluation of differences of aortic lumen diameters as well as differences in length comparing both acquisition techniques at every possible location. METHODS AND MATERIALS Sixteen healthy volunteers underwent 1.5-T-MR Angiography (MRA). For each volunteer, two different MR sequences were performed, CE-MRA: gradient echo Turbo FLASH sequence and N-MRA: respiratory-and-cardiac-gated, T2-weighted 3D SSFP. Datasets were segmented using a 3D model-based ellipse-fitting approach with a single seed point placed manually above the celiac trunk. The segmented volumes were manually cropped from left subclavian artery to celiac trunk to avoid error due to side branches. Diameters, volumes and centerline length were computed for intraindividual comparison. For statistical analysis the Wilcoxon-Signed-Ranked-Test was used. RESULTS Average centerline length obtained based on N-MRA was 239.0±23.4 mm compared to 238.6±23.5 mm for CE-MRA without significant difference (P=0.877). Average maximum diameter obtained based on N-MRA was 25.7±3.3 mm compared to 24.1±3.2 mm for CE-MRA (P<0.001). In agreement with the difference in diameters, volumes obtained based on N-MRA (100.1±35.4 cm(3)) were consistently and significantly larger compared to CE-MRA (89.2±30.0 cm(3)) (P<0.001). CONCLUSIONS 3D morphometry shows highly similar centerline lengths for N-MRA and CE-MRA, but systematically higher diameters and volumes for N-MRA.

Aortic morphometry at endograft position as assessed by 3D image analysis affects risk of type I endoleak formation after TEVAR

PURPOSE The purpose of this study was to identify morphologic factors affecting type I endoleak formation and bird-beak configuration after thoracic endovascular aortic repair (TEVAR). METHODS Computed tomography (CT) data of 57 patients (40 males; median age, 66 years) undergoing TEVAR for thoracic aortic aneurysm (34 TAA, 19 TAAA) or penetrating aortic ulcer (n = 4) between 2001 and 2010 were retrospectively reviewed. In 28 patients, the Gore TAG® stent-graft was used, followed by the Medtronic Valiant® in 16 cases, the Medtronic Talent® in 8, and the Cook Zenith® in 5 cases. Proximal landing zone (PLZ) was in zone 1 in 13, zone 2 in 13, zone 3 in 23, and zone 4 in 8 patients. In 14 patients (25%), the procedure was urgent or emergent. In each case, pre- and postoperative CT angiography was analyzed using a dedicated image processing workstation and complimentary in-house developed software based on a 3D cylindrical intensity model to calculate aortic arch angulation and conicity of the landing zones (LZ). RESULTS Primary type Ia endoleak rate was 12% (7/57) and subsequent re-intervention rate was 86% (6/7). Left subclavian artery (LSA) coverage (p = 0.036) and conicity of the PLZ (5.9 vs. 2.6 mm; p = 0.016) were significantly associated with an increased type Ia endoleak rate. Bird-beak configuration was observed in 16 patients (28%) and was associated with a smaller radius of the aortic arch curvature (42 vs. 65 mm; p = 0.049). Type Ia endoleak was not associated with a bird-beak configuration (p = 0.388). Primary type Ib endoleak rate was 7% (4/57) and subsequent re-intervention rate was 100%. Conicity of the distal LZ was associated with an increased type Ib endoleak rate (8.3 vs. 2.6 mm; p = 0.038). CONCLUSIONS CT-based 3D aortic morphometry helps to identify risk factors of type I endoleak formation and bird-beak configuration during TEVAR. These factors were LSA coverage and conicity within the landing zones for type I endoleak formation and steep aortic angulation for bird-beak configuration.

Motion-insensitive determination of B1+ amplitudes based on the bloch-siegert shift in single voxels of moving organs including the human heart.

PURPOSE To reliably determine the amplitude of the transmit radiofrequency ( B1+) field in moving organs like the liver and heart, where most current techniques are usually not feasible. METHODS B1+ field measurement based on the Bloch-Siegert shift induced by a pair of Fermi pulses in a double-triggered modified Point RESolved Spectroscopy (PRESS) sequence with motion-compensated crusher gradients has been developed. Performance of the sequence was tested in moving phantoms and in muscle, liver, and heart of six healthy volunteers each, using different arrangements of transmit/receive coils. RESULTS B1+ determination in a moving phantom was almost independent of type and amplitude of the motion and agreed well with theory. In vivo, repeated measurements led to very small coefficients of variance (CV) if the amplitude of the Fermi pulse was chosen above an appropriate level (CV in muscle 0.6%, liver 1.6%, heart 2.3% with moderate amplitude of the Fermi pulses and 1.2% with stronger Fermi pulses). CONCLUSION The proposed sequence shows a very robust determination of B1+ in a single voxel even under challenging conditions (transmission with a surface coil or measurements in the heart without breath-hold). Magn Reson Med, 2015. © 2015 Wiley Periodicals, Inc.

MorphCol 2004-2013. A collection of Fortran 77 programs for morphometry

This is a technical description in html format of simple fortran programs for Macintosh for the morphometric analysis of tests planktonic foraminifera under reflected light, with special focus on the Neogene group of Globorotalia menardii. The second part of this report gives information and performance tests about the development of AMOR (Automated Measurement system for the mORphometry of microfossils). AMOR is Windows based and helps to orientate and collect digital images of menardiform globorotalids. The above fortran programs may be useful to extract and analyze some morphometric parameters from images collected with AMOR. After unzipping the archive file please open the Start.html file using a common web browser like firefox. In case of any questions or problems, please contact Michael W. Knappertsbusch (mailto:michael.knappertsbusch@unibas.ch).

Geochemistry and morphometry on planktonic foraminifera

About one third of the anthropogenic carbon dioxide (CO2) released into the atmosphere in the past two centuries has been taken up by the ocean. As CO2 invades the surface ocean, carbonate ion concentrations and pH are lowered. Laboratory studies indicate that this reduces the calcification rates of marine calcifying organisms, including planktic foraminifera. Such a reduction in calcification resulting from anthropogenic CO2 emissions has not been observed, or quantified in the field yet. Here we present the findings of a study in the Western Arabian Sea that uses shells of the surface water dwelling planktic foraminifer Globigerinoides ruber in order to test the hypothesis that anthropogenically induced acidification has reduced shell calcification of this species. We found that light, thin-walled shells from the surface sediment are younger (based on 14C and d13C measurements) than the heavier, thicker-walled shells. Shells in the upper, bioturbated, sediment layer were significantly lighter compared to shells found below this layer. These observations are consistent with a scenario where anthropogenically induced ocean acidification reduced the rate at which foraminifera calcify, resulting in lighter shells. On the other hand, we show that seasonal upwelling in the area also influences their calcification and the stable isotope (d13C and d18O) signatures recorded by the foraminifera shells. Plankton tow and sediment trap data show that lighter shells were produced during upwelling and heavier ones during non-upwelling periods. Seasonality alone, however, cannot explain the 14C results, or the increase in shell weight below the bioturbated sediment layer. We therefore must conclude that probably both the processes of acidification and seasonal upwelling are responsible for the presence of light shells in the top of the sediment and the age difference between thick and thin specimens.