ANNA: A new prediction method for bioassessment programs

1. Cluster analysis of reference sites with similar biota is the initial step in creating River Invertebrate Prediction and Classification System (RIVPACS) and similar river bioassessment models such as Australian River Assessment System (AUSRIVAS). This paper describes and tests an alternative prediction method, Assessment by Nearest Neighbour Analysis (ANNA), based on the same philosophy as RIVPACS and AUSRIVAS but without the grouping step that some people view as artificial. 2. The steps in creating ANNA models are: (i) weighting the predictor variables using a multivariate approach analogous to principal axis correlations, (ii) calculating the weighted Euclidian distance from a test site to the reference sites based on the environmental predictors, (iii) predicting the faunal composition based on the nearest reference sites and (iv) calculating an observed/expected (O/E) analogous to RIVPACS/AUSRIVAS. 3. The paper compares AUSRIVAS and ANNA models on 17 datasets representing a variety of habitats and seasons. First, it examines each model's regressions for Observed versus Expected number of taxa, including the r(2), intercept and slope. Second, the two models' assessments of 79 test sites in New Zealand are compared. Third, the models are compared on test and presumed reference sites along a known trace metal gradient. Fourth, ANNA models are evaluated for western Australia, a geographically distinct region of Australia. The comparisons demonstrate that ANNA and AUSRIVAS are generally equivalent in performance, although ANNA turns out to be potentially more robust for the O versus E regressions and is potentially more accurate on the trace metal gradient sites. 4. The ANNA method is recommended for use in bioassessment of rivers, at least for corroborating the results of the well established AUSRIVAS- and RIVPACS-type models, if not to replace them.

Measuring macromolecular diffusion using heteronuclear multiple-quantum pulsed-field-gradient NMR

We have previously shown that H-1 pulsed-field-gradient (PFG) NMR spectroscopy provides a facile method for monitoring protein self-association and can be used, albeit with some caveats, to measure the apparent molecular mass of the diffusant [Dingley et al. (1995) J. Biomol. NMR, 6, 321-328]. In this paper we show that, for N-15-labelled proteins, selection of H-1-N-15 multiple-quantum (MQ) coherences in PFG diffusion experiments provides several advantages over monitoring H-1 single-quantum (SQ) magnetization. First, the use of a gradient-selected MQ filter provides a convenient means of suppressing resonances from both the solvent and unlabelled solutes. Second, H-1-N-15 zero-quantum coherence dephases more rapidly than H-1 SQ coherence under the influence of a PFG. This allows the diffusion coefficients of larger proteins to be measured more readily. Alternatively, the gradient length and/or the diffusion delay may be decreased, thereby reducing signal losses from relaxation. In order to extend the size of macromolecules to which these experiments can be applied, we have developed a new MQ PFG diffusion experiment in which the magnetization is stored as longitudinal two-spin order for most of the diffusion period, thus minimizing sensitivity losses due to transverse relaxation and J-coupling evolution.

Validation of a Cutoff Value on Echo Doppler Analysis to Replace Right Heart Catheterization During Pulmonary Hypertension Evaluation in Heart Transplant Candidates

Background. Heart transplantation (OHT) has traditionally been contraindicated in the presence of severe pulmonary hypertension (PH), as detected by right heart catheterization. Noninvasive methods are still not reliably accurate to make this evaluation. Objectives. Determine the efficacy of echo Doppler analysis for the diagnosis of severe PH. Methods. One hundred thirty patients (mean age = 42 +/- 15 years, 82 men) showed severe left ventricular dysfunction (mean ejection fraction = 29 +/- 12%; functional class III-IV). We excluded patients with atrial fibrillation, heart failure secondary to congenital disease, and valvulopathy. The pulmonary parameters defined as severe PH were: systolic pulmonary artery pressure (sPAP) >= 60 mm Hg; a mean transpulmonary gradient >= 15; or pulmonary vascular resistance >= 5 Wood units. Patients underwent a right heart catheterization using a Swan-Ganz catheter to measure hemodynamic parameters and to noninvasively estimate right-sided pressures from spectral Doppler recordings of tricuspid regurgitation velocity (right ventricular systolic pressure [RVsP]). A Pearson correlation of sPAP was obtained with RVsP by; the sensitivity of RVsP for the diagnosis of PH was determined by a receiver operating characteristic (ROC) curve. Results. A good correlation between sPAP and RVsP was obtained by Pearson correlation analysis (r = 0.64; 95% confidence interval [CI] 0.50-0.75; P < .001). The ROC curve analysis showed a sensitivity of 100%, a specificity of 37.2%, (95% CI 0.69-0.83, P < .0001) of a RVsP < 45 mm Hg (cutoff) on the exclusion of severe PH. Conclusions. The cutoff of RVsP < 45 mm Hg, on noninvasive echo Doppler evaluation of PH is an efficient method to replace invasive heart catheterization in OHT candidates.

A comparison of low-storage strategies for spectral analysis in dissipative systems: filter diagonalisation in the Lanczos representation and harmonic inversion of the Chebychev-order-domain autocorrelation function

We compare the performance of two different low-storage filter diagonalisation (LSFD) strategies in the calculation of complex resonance energies of the HO2, radical. The first is carried out within a complex-symmetric Lanczos subspace representation [H. Zhang, S.C. Smith, Phys. Chem. Chem. Phys. 3 (2001) 2281]. The second involves harmonic inversion of a real autocorrelation function obtained via a damped Chebychev recursion [V.A. Mandelshtam, H.S. Taylor, J. Chem. Phys. 107 (1997) 6756]. We find that while the Chebychev approach has the advantage of utilizing real algebra in the time-consuming process of generating the vector recursion, the Lanczos, method (using complex vectors) requires fewer iterations, especially for low-energy part of the spectrum. The overall efficiency in calculating resonances for these two methods is comparable for this challenging system. (C) 2001 Elsevier Science B.V. All rights reserved.

Time evolution in the unimolecular master equation at low temperatures: full spectral solution with scalable iterative methods and high precision

A new method is presented to determine an accurate eigendecomposition of difficult low temperature unimolecular master equation problems. Based on a generalisation of the Nesbet method, the new method is capable of achieving complete spectral resolution of the master equation matrix with relative accuracy in the eigenvectors. The method is applied to a test case of the decomposition of ethane at 300 K from a microcanonical initial population with energy transfer modelled by both Ergodic Collision Theory and the exponential-down model. The fact that quadruple precision (16-byte) arithmetic is required irrespective of the eigensolution method used is demonstrated. (C) 2001 Elsevier Science B.V. All rights reserved.

A novel target-field method for finite-length magnetic resonance shim coils: I. Zonal shims

This paper presents a new approach for the design of genuinely finite-length shim and gradient coils, intended for use in magnetic resonance imaging equipment. A cylindrical target region is located asymmetrically, at an arbitrary position within a coil of finite length. A desired target field is specified on the surface of that region, and a method is given that enables winding patterns on the surface of the coil to be designed, to produce the desired field at the inner target region. The method uses a minimization technique combined with regularization, to find the current density on the surface of the coil. The method is illustrated for linear, quadratic and cubic magnetic target fields located asymmetrically within a finite-length coil.

Application of the maximum entropy method to the (2 + 1)D four-fermion model

We investigate spectral functions extracted using the maximum entropy method from correlators measured in lattice simulations of the (2+1)-dimensional four-fermion model. This model is particularly interesting because it has both a chirally broken phase with a rich spectrum of mesonic bound states and a symmetric phase where there are only resonances. In the broken phase we study the elementary fermion, pion, sigma, and massive pseudoscalar meson; our results confirm the Goldstone nature of the π and permit an estimate of the meson binding energy. We have, however, seen no signal of σ→ππ decay as the chiral limit is approached. In the symmetric phase we observe a resonance of nonzero width in qualitative agreement with analytic expectations; in addition the ultraviolet behavior of the spectral functions is consistent with the large nonperturbative anomalous dimension for fermion composite operators expected in this model.

Finite difference time domain (FDTD) method for modeling the effect of switched gradients on the human body in MRI

Numerical modeling of the eddy currents induced in the human body by the pulsed field gradients in MRI presents a difficult computational problem. It requires an efficient and accurate computational method for high spatial resolution analyses with a relatively low input frequency. In this article, a new technique is described which allows the finite difference time domain (FDTD) method to be efficiently applied over a very large frequency range, including low frequencies. This is not the case in conventional FDTD-based methods. A method of implementing streamline gradients in FDTD is presented, as well as comparative analyses which show that the correct source injection in the FDTD simulation plays a crucial rule in obtaining accurate solutions. In particular, making use of the derivative of the input source waveform is shown to provide distinct benefits in accuracy over direct source injection. In the method, no alterations to the properties of either the source or the transmission media are required. The method is essentially frequency independent and the source injection method has been verified against examples with analytical solutions. Results are presented showing the spatial distribution of gradient-induced electric fields and eddy currents in a complete body model.

On the induced electric field gradients in the human body for magnetic stimulation by gradient coils in MRI

Prior theoretical studies indicate that the negative spatial derivative of the electric field induced by magnetic stimulation may he one of the main factors contributing to depolarization of the nerve fiber. This paper studies this parameter for peripheral nerve stimulation (PNS) induced by time.-varying gradient fields during MRI scans. The numerical calculations are based on an efficient, quasi-static, finite-difference scheme and an anatomically realistic human, full-body model. Whole-body cylindrical and planar gradient sets in MRI systems and various input signals have been explored. The spatial distributions of the induced electric field and their gradients are calculated and attempts are made to correlate these areas with reported experimental stimulation data. The induced electrical field pattern is similar for both the planar coils and cylindrical coils. This study provides some insight into the spatial characteristics of the induced field gradients for PNS in MRI, which may be used to further evaluate the sites where magnetic stimulation is likely to occur and to optimize gradient coil design.

A novel target-field method for magnetic resonance shim coils: III. Shielded zonal and tesseral coils

This paper continues the development of a new approach for the design of shim and gradient coils, used in magnetic resonance imaging (MRI) applications. A cylindrical primary coil of radius a and length 2L is placed inside a co-axial shield cylinder of radius b. An active shielding strategy is used to create a desired target field at an arbitrarily specified (cylindrical) location within the primary coil, and to annul the field at a certain radius outside the shield. The form of the interior target field may be chosen arbitrarily by the designer, although zonal and tesseral harmonics are typically used in MRI applications. The method presented here designs coil windings on both the primary and shielding cylinders, to produce fields that conform to the specified interior target field and the annulled field exterior to the shield. An additional feature of the method presented here is that the target field inside the primary coil is matched at two different radii, to improve overall accuracy. The method is illustrated by designing several shielded shim coils, for creating higher order tesseral fields located asymmetrically within the coil. The simpler case of pure zonal fields is discussed separately and applied to the design of some higher order shielded coils.

Implementação em hardware reconfigurável de método de separação de dados hiperespetrais

Relatório do Trabalho Final de Mestrado para obtenção do grau de Mestre em Engenharia de Electrónica e Telecomunicações

Wetland Habitat Studies using various Classification Techniques on Multi-Spectral Landsat Imagery: Case study: Tram chim National Park, Dong Thap Vietnam

Dissertation submitted in partial fulfilment of the requirements for the Degree of Master of Science in Geospatial Technologies

Spectral solution for the air stripping pollutants removal dynamic model with non linear steady state conditions

This work deals with the numerical simulation of air stripping process for the pre-treatment of groundwater used in human consumption. The model established in steady state presents an exponential solution that is used, together with the Tau Method, to get a spectral approach of the solution of the system of partial differential equations associated to the model in transient state.

A comprehensive high-throughput FTIR spectroscopy-based method for evaluating the transfection event: estimating the transfection efficiency and extracting associated metabolic responses

Reporter genes are routinely used in every laboratory for molecular and cellular biology for studying heterologous gene expression and general cellular biological mechanisms, such as transfection processes. Although well characterized and broadly implemented, reporter genes present serious limitations, either by involving time-consuming procedures or by presenting possible side effects on the expression of the heterologous gene or even in the general cellular metabolism. Fourier transform mid-infrared (FT-MIR) spectroscopy was evaluated to simultaneously analyze in a rapid (minutes) and high-throughput mode (using 96-wells microplates), the transfection efficiency, and the effect of the transfection process on the host cell biochemical composition and metabolism. Semi-adherent HEK and adherent AGS cell lines, transfected with the plasmid pVAX-GFP using Lipofectamine, were used as model systems. Good partial least squares (PLS) models were built to estimate the transfection efficiency, either considering each cell line independently (R 2 ≥ 0.92; RMSECV ≤ 2 %) or simultaneously considering both cell lines (R 2 = 0.90; RMSECV = 2 %). Additionally, the effect of the transfection process on the HEK cell biochemical and metabolic features could be evaluated directly from the FT-IR spectra. Due to the high sensitivity of the technique, it was also possible to discriminate the effect of the transfection process from the transfection reagent on KEK cells, e.g., by the analysis of spectral biomarkers and biochemical and metabolic features. The present results are far beyond what any reporter gene assay or other specific probe can offer for these purposes.

Parallel GPU architecture for hyperspectral unmixing based on augmented Lagrangian method

Hyperspectral imaging has become one of the main topics in remote sensing applications, which comprise hundreds of spectral bands at different (almost contiguous) wavelength channels over the same area generating large data volumes comprising several GBs per flight. This high spectral resolution can be used for object detection and for discriminate between different objects based on their spectral characteristics. One of the main problems involved in hyperspectral analysis is the presence of mixed pixels, which arise when the spacial resolution of the sensor is not able to separate spectrally distinct materials. Spectral unmixing is one of the most important task for hyperspectral data exploitation. However, the unmixing algorithms can be computationally very expensive, and even high power consuming, which compromises the use in applications under on-board constraints. In recent years, graphics processing units (GPUs) have evolved into highly parallel and programmable systems. Specifically, several hyperspectral imaging algorithms have shown to be able to benefit from this hardware taking advantage of the extremely high floating-point processing performance, compact size, huge memory bandwidth, and relatively low cost of these units, which make them appealing for onboard data processing. In this paper, we propose a parallel implementation of an augmented Lagragian based method for unsupervised hyperspectral linear unmixing on GPUs using CUDA. The method called simplex identification via split augmented Lagrangian (SISAL) aims to identify the endmembers of a scene, i.e., is able to unmix hyperspectral data sets in which the pure pixel assumption is violated. The efficient implementation of SISAL method presented in this work exploits the GPU architecture at low level, using shared memory and coalesced accesses to memory.