Active deformation fields: dense deformation field estimation for atlas-based segmentation using the active contour framework.

This paper presents a new and original variational framework for atlas-based segmentation. The proposed framework integrates both the active contour framework, and the dense deformation fields of optical flow framework. This framework is quite general and encompasses many of the state-of-the-art atlas-based segmentation methods. It also allows to perform the registration of atlas and target images based on only selected structures of interest. The versatility and potentiality of the proposed framework are demonstrated by presenting three diverse applications: In the first application, we show how the proposed framework can be used to simulate the growth of inconsistent structures like a tumor in an atlas. In the second application, we estimate the position of nonvisible brain structures based on the surrounding structures and validate the results by comparing with other methods. In the final application, we present the segmentation of lymph nodes in the Head and Neck CT images, and demonstrate how multiple registration forces can be used in this framework in an hierarchical manner.

Comparing interval breast cancer rates in Norway and North Carolina: results and challenges.

OBJECTIVE: To compare interval breast cancer rates (ICR) between a biennial organized screening programme in Norway and annual opportunistic screening in North Carolina (NC) for different conceptualizations of interval cancer. SETTING: Two regions with different screening practices and performance. METHODS: 620,145 subsequent screens (1996-2002) performed in women aged 50-69 and 1280 interval cancers were analysed. Various definitions and quantification methods for interval cancers were compared. RESULTS: ICR for one year follow-up were lower in Norway compared with NC both when the rate was based on all screens (0.54 versus 1.29 per 1000 screens), negative final assessments (0.54 versus 1.29 per 1000 screens), and negative screening assessments (0.53 versus 1.28 per 1000 screens). The rate of ductal carcinoma in situ was significantly lower in Norway than in NC for cases diagnosed in both the first and second year after screening. The distributions of histopathological tumour size and lymph node involvement in invasive cases did not differ between the two regions for interval cancers diagnosed during the first year after screening. In contrast, in the second year after screening, tumour characteristics remained stable in Norway but became prognostically more favorable in NC. CONCLUSION: Even when applying a common set of definitions of interval cancer, the ICR was lower in Norway than in NC. Different definitions of interval cancer did not influence the ICR within Norway or NC. Organization of screening and screening performance might be major contributors to the differences in ICR between Norway and NC.

Grade and Cross Slope Estimation from LIDAR-based Surface Models, 2003

Many transportation agencies maintain grade as an attribute in roadway inventory databases; however, the information is often in an aggregated format. Cross slope is rarely included in large roadway inventories. Accurate methods available to collect grade and cross slope include global positioning systems, traditional surveying, and mobile mapping systems. However, most agencies do not have the resources to utilize these methods to collect grade and cross slope on a large scale. This report discusses the use of LIDAR to extract roadway grade and cross slope for large-scale inventories. Current data collection methods and their advantages and disadvantages are discussed. A pilot study to extract grade and cross slope from a LIDAR data set, including methodology, results, and conclusions, is presented. This report describes the regression methodology used to extract and evaluate the accuracy of grade and cross slope from three dimensional surfaces created from LIDAR data. The use of LIDAR data to extract grade and cross slope on tangent highway segments was evaluated and compared against grade and cross slope collected using an automatic level for 10 test segments along Iowa Highway 1. Grade and cross slope were measured from a surface model created from LIDAR data points collected for the study area. While grade could be estimated to within 1%, study results indicate that cross slope cannot practically be estimated using a LIDAR derived surface model.

Estimation of forest canopy cover : a comparison of field measurement techniques

Abstract

Intermittent atrial tachycardia facilitates atrial fibrillation by a shortening of activation recovery interval

Introduction: We recently observed in a chronic ovine model that a shortening of action potential duration (APD) as assessed by the activation recovery interval (ARI) may be a mechanism whereby pacing-induced atrial tachycardia (PIAT) facilitates atrial fibrillation (AF), mediated by a return to 1:1 atrial capture after the effective refractory period has been reached. The aim of the present study is to evaluate the effect of long term intermittent burst pacing on ARI before induction of AF.Methods: We specifically developed a chronic ovine model of PIAT using two pacemakers (PM) each with a right atrial (RA) lead separated by ∼2cm. The 1st PM (Vitatron T70) was used to record a broadband unipolar RA EGM (800 Hz, 0.4 Hz high pass filter). The 2nd was used to deliver PIAT during electrophysiological protocols at decremental pacing CL (400 beats, from 400 to 110ms) and long term intermittent RA burst pacing to promote electrical remodeling (5s of burst followed by 2s of sinus rhythm) until onset of sustained AF. ARI was defined as the time difference between the peak of the atrial repolarization wave and the first atrial depolarization. The mean ARIs of paired sequences (before and after remodeling), each consisting of 20 beats were compared.Results: As shown in the figure, ARIs (n=4 sheep, 46 recordings) decreased post remodeling compared to baseline (86±19 vs 103±12 ms, p<0.05). There was no difference in atrial structure as assessed by light microscopy between control and remodeled sheep.Conclusions: Using standard pacemaker technology, atrial ARIs as a surrogate of APDs were successfully measured in vivo during the electrical remodeling process leading to AF. The facilitation of AF by PIAT mimicking salvos from pulmonary veins is heralded by a significant shortening of ARI.

Software for interpolation of vegetative growth of yerba mate plants in 3D

The objective of this work was to build mock-ups of complete yerba mate plants in several stages of development, using the InterpolMate software, and to compute photosynthesis on the interpolated structure. The mock-ups of yerba-mate were first built in the VPlants software for three growth stages. Male and female plants grown in two contrasting environments (monoculture and forest understory) were considered. To model the dynamic 3D architecture of yerba-mate plants during the biennial growth interval between two subsequent prunings, data sets of branch development collected in 38 dates were used. The estimated values obtained from the mock-ups, including leaf photosynthesis and sexual dimorphism, are very close to those observed in the field. However, this similarity was limited to reconstructions that included growth units from original data sets. The modeling of growth dynamics enables the estimation of photosynthesis for the entire yerba mate plant, which is not easily measurable in the field. The InterpolMate software is efficient for building yerba mate mock-ups.

Effects of computing parameters and measurement locations on the estimation of 3D NPS in non-stationary MDCT images.

The goal of this study was to investigate the impact of computing parameters and the location of volumes of interest (VOI) on the calculation of 3D noise power spectrum (NPS) in order to determine an optimal set of computing parameters and propose a robust method for evaluating the noise properties of imaging systems. Noise stationarity in noise volumes acquired with a water phantom on a 128-MDCT and a 320-MDCT scanner were analyzed in the spatial domain in order to define locally stationary VOIs. The influence of the computing parameters in the 3D NPS measurement: the sampling distances bx,y,z and the VOI lengths Lx,y,z, the number of VOIs NVOI and the structured noise were investigated to minimize measurement errors. The effect of the VOI locations on the NPS was also investigated. Results showed that the noise (standard deviation) varies more in the r-direction (phantom radius) than z-direction plane. A 25 × 25 × 40 mm(3) VOI associated with DFOV = 200 mm (Lx,y,z = 64, bx,y = 0.391 mm with 512 × 512 matrix) and a first-order detrending method to reduce structured noise led to an accurate NPS estimation. NPS estimated from off centered small VOIs had a directional dependency contrary to NPS obtained from large VOIs located in the center of the volume or from small VOIs located on a concentric circle. This showed that the VOI size and location play a major role in the determination of NPS when images are not stationary. This study emphasizes the need for consistent measurement methods to assess and compare image quality in CT.

The role of interval nodes in sentinel lymph node mapping and dissection for melanoma patients.

In sentinel node (SN) biopsy, an interval SN is defined as a lymph node or group of lymph nodes located between the primary melanoma and an anatomically well-defined lymph node group directly draining the skin. As shown in previous reports, these interval SNs seem to be at the same metastatic risk as are SNs in the usual, classic areas. This study aimed to review the incidence, lymphatic anatomy, and metastatic risk of interval SNs. METHODS: SN biopsy was performed at a tertiary center by a single surgical team on a cohort of 402 consecutive patients with primary melanoma. The triple technique of localization was used-that is, lymphoscintigraphy, blue dye, and gamma-probe. Otolaryngologic melanoma and mucosal melanoma were excluded from this analysis. SNs were examined by serial sectioning and immunohistochemistry. All patients with metastatic SNs were recommended to undergo a radical selective lymph node dissection. RESULTS: The primary locations of the melanomas included the trunk (188), an upper limb (67), or a lower limb (147). Overall, 97 (24.1%) of the 402 SNs were metastatic. Interval SNs were observed in 18 patients, in all but 2 of whom classic SNs were also found. The location of the primary was truncal in 11 (61%) of the 18, upper limb in 5, and lower limb in 2. One patient with a dorsal melanoma had drainage exclusively in a cervicoscapular area that was shown on removal to contain not lymph node tissue but only a blue lymph channel without tumor cells. Apart from the interval SN, 13 patients had 1 classic SN area and 3 patients 2 classic SN areas. Of the 18 patients, 2 had at least 1 metastatic interval SN and 2 had a classic SN that was metastatic; overall, 4 (22.2%) of 18 patients were node-positive. CONCLUSION: We found that 2 of 18 interval SNs were metastatic: This study showed that preoperative lymphoscintigraphy must review all known lymphatic areas in order to exclude an interval SN.

Estimating parameters for generalized mass action models using constraint propagation.

As modern molecular biology moves towards the analysis of biological systems as opposed to their individual components, the need for appropriate mathematical and computational techniques for understanding the dynamics and structure of such systems is becoming more pressing. For example, the modeling of biochemical systems using ordinary differential equations (ODEs) based on high-throughput, time-dense profiles is becoming more common-place, which is necessitating the development of improved techniques to estimate model parameters from such data. Due to the high dimensionality of this estimation problem, straight-forward optimization strategies rarely produce correct parameter values, and hence current methods tend to utilize genetic/evolutionary algorithms to perform non-linear parameter fitting. Here, we describe a completely deterministic approach, which is based on interval analysis. This allows us to examine entire sets of parameters, and thus to exhaust the global search within a finite number of steps. In particular, we show how our method may be applied to a generic class of ODEs used for modeling biochemical systems called Generalized Mass Action Models (GMAs). In addition, we show that for GMAs our method is amenable to the technique in interval arithmetic called constraint propagation, which allows great improvement of its efficiency. To illustrate the applicability of our method we apply it to some networks of biochemical reactions appearing in the literature, showing in particular that, in addition to estimating system parameters in the absence of noise, our method may also be used to recover the topology of these networks.

Soybean crop area estimation by Modis/Evi data

The objective of this work was to develop a procedure to estimate soybean crop areas in Rio Grande do Sul state, Brazil. Estimations were made based on the temporal profiles of the enhanced vegetation index (Evi) calculated from moderate resolution imaging spectroradiometer (Modis) images. The methodology developed for soybean classification was named Modis crop detection algorithm (MCDA). The MCDA provides soybean area estimates in December (first forecast), using images from the sowing period, and March (second forecast), using images from the sowing and maximum crop development periods. The results obtained by the MCDA were compared with the official estimates on soybean area of the Instituto Brasileiro de Geografia e Estatística. The coefficients of determination ranged from 0.91 to 0.95, indicating good agreement between the estimates. For the 2000/2001 crop year, the MCDA soybean crop map was evaluated using a soybean crop map derived from Landsat images, and the overall map accuracy was approximately 82%, with similar commission and omission errors. The MCDA was able to estimate soybean crop areas in Rio Grande do Sul State and to generate an annual thematic map with the geographic position of the soybean fields. The soybean crop area estimates by the MCDA are in good agreement with the official agricultural statistics.

Stratified estimation of forest inventory variables using spatially summarized stratifications

Abstract

Bone mineral density (BMD), micro-architecture estimation (TBS) and vertebral fracture assessment (VFA) extracted from a single DXA device in combination with clinical risk factors improve significantly the identification of women at high risk of fracture.

Osteoporosis (OP) is a systemic skeletal disease characterized by a low bone mineral density (BMD) and a micro-architectural (MA) deterioration. Clinical risk factors (CRF) are often used as a MA approximation. MA is yet evaluable in daily practice by the trabecular bone score (TBS) measure. TBS is very simple to obtain, by reanalyzing a lumbar DXA-scan. TBS has proven to have diagnosis and prognosis values, partially independent of CRF and BMD. The aim of the OsteoLaus cohort is to combine in daily practice the CRF and the information given by DXA (BMD, TBS and vertebral fracture assessment (VFA)) to better identify women at high fracture risk. The OsteoLaus cohort (1400 women 50 to 80 years living in Lausanne, Switzerland) started in 2010. This study is derived from the cohort COLAUS who started in Lausanne in 2003. The main goal of COLAUS is to obtain information on the epidemiology and genetic determinants of cardiovascular risk in 6700 men and women. CRF for OP, bone ultrasound of the heel, lumbar spine and hip BMD, VFA by DXA and MA evaluation by TBS are recorded in OsteoLaus. Preliminary results are reported. We included 631 women: mean age 67.4 ± 6.7 years, BMI 26.1 ± 4.6, mean lumbar spine BMD 0.943 ± 0.168 (T-score − 1.4 SD), and TBS 1.271 ± 0.103. As expected, correlation between BMD and site matched TBS is low (r2 = 0.16). Prevalence of VFx grade 2/3, major OP Fx and all OP Fx is 8.4%, 17.0% and 26.0% respectively. Age- and BMI-adjusted ORs (per SD decrease) are 1.8 (1.2-2.5), 1.6 (1.2-2.1), and 1.3 (1.1-1.6) for BMD for the different categories of fractures and 2.0 (1.4-3.0), 1.9 (1.4-2.5), and 1.4 (1.1-1.7) for TBS respectively. Only 32 to 37% of women with OP Fx have a BMD < − 2.5 SD or a TBS < 1.200. If we combine a BMD < − 2.5 SD or a TBS < 1.200, 54 to 60% of women with an osteoporotic Fx are identified. As in the already published studies, these preliminary results confirm the partial independence between BMD and TBS. More importantly, a combination of TBS subsequent to BMD increases significantly the identification of women with prevalent OP Fx which would have been misclassified by BMD alone. For the first time we are able to have complementary information about fracture (VFA), density (BMD), micro- and macro architecture (TBS and HAS) from a simple, low ionizing radiation and cheap device: DXA. Such complementary information is very useful for the patient in the daily practice and moreover will likely have an impact on cost effectiveness analysis.