Regularization methods for the solution of a nonlinear least-squares problem in tomography

In this work we study a polyenergetic and multimaterial model for the breast image reconstruction in Digital Tomosynthesis, taking into consideration the variety of the materials forming the object and the polyenergetic nature of the X-rays beam. The modelling of the problem leads to the resolution of a high-dimensional nonlinear least-squares problem that, due to its nature of inverse ill-posed problem, needs some kind of regularization. We test two main classes of methods: the Levenberg-Marquardt method (together with the Conjugate Gradient method for the computation of the descent direction) and two limited-memory BFGS-like methods (L-BFGS). We perform some experiments for different values of the regularization parameter (constant or varying at each iteration), tolerances and stop conditions. Finally, we analyse the performance of the several methods comparing relative errors, iterations number, times and the qualities of the reconstructed images.

Development of an LC-MS/MS method for the quantitation of 55 compounds prescribed in combined cardiovascular therapy

This paper reports an LC-MS/MS method with positive electrospray ionization for the screening of commonly prescribed cardiovascular drugs in human plasma, including compounds with antihypertensive (57), antidiabetic (12), hypolipemiant (5), anticoagulant (2) and platelet anti-aggregation (2) effects. Sample treatment consisted of a simple protein precipitation with MeOH/0.1 M ZnSO₄ (4:1, v/v) solution after the addition of internal standard, followed by evaporation and reconstitution. Analytes separation was performed on a Polar-RP column (150 m x 2 mm, 4 μm) using a gradient elution of 15 min. The MS system was operated in MRM mode, monitoring one quantitation and one confirmation transition for each analyte. The recovery of the protein precipitation step ranged from 50 to 70% for most of the compounds, while some were considerably affected by matrix effects. Since several analytes fulfilled the linearity, accuracy and precision values required by the ICH guidelines, the method proved to be suitable for their quantitative analysis. The limits of quantitation varied from 0.38 to 9.1 μg/L and the limits of detection from 0.12 to 5.34 μg/L. The method showed to be suitable for the detection of plasma samples of patients under cardiovascular treatment with the studied drugs, and for 55 compounds reliable quantitative results could be obtained.

LC-MS/MS screening method for designer amphetamines, tryptamines, and piperazines in serum

Since the late 1990s and early 2000s, derivatives of well-known designer drugs as well as new psychoactive compounds have been sold on the illicit drug market and have led to intoxications and fatalities. The LC-MS/MS screening method presented covers 31 new designer drugs as well as cathinone, methcathinone, phencyclidine, and ketamine which were included to complete the screening spectrum. All but the last two are modified molecular structures of amphetamine, tryptamine, or piperazine. Among the amphetamine derivatives are cathinone, methcathinone, 3,4-DMA, 2,5-DMA, DOB, DOET, DOM, ethylamphetamine, MDDMA, 4-MTA, PMA, PMMA, 3,4,5-TMA, TMA-6 and members of the 2C group: 2C-B, 2C-D, 2C-H, 2C-I, 2C-P, 2C-T-2, 2C-T-4, and 2C-T-7. AMT, DPT, DiPT, MiPT, DMT, and 5MeO-DMT are contained in the tryptamine group, BZP, MDBP, TFMPP, mCPP, and MeOPP in the piperazine group. Using an Applied Biosystems LC-MS/MS API 365 TurboIonSpray it is possible to identify all 35 substances. After addition of internal standards and mixed-mode solid-phase extraction the analytes are separated using a Synergi Polar RP column and gradient elution with 1 mM ammonium formate and methanol/0.1% formic acid as mobile phases A and B. Data acquisition is performed in MRM mode with positive electro spray ionization. The assay is selective for all tested substances. Limits of detection were determined by analyzing S/N-ratios and are between 1.0 and 5.0 ng/mL. Matrix effects lie between 65% and 118%, extraction efficiencies range from 72% to 90%.

Hydrogen bonding of nucleotide base pairs: Application of the PM3 method

The ability of the pm3 semiempirical quantum mechanical method to reproduce hydrogen bonding in nucleotide base pairs was assessed. Results of pm3 calculations on the nucleotides 2′-deoxyadenosine 5′-monophosphate (pdA), 2′-deoxyguanosine 5′-monophosphate (pdG), 2′-deoxycytidine 5′-monophosphate (pdC), and 2′-deoxythymidine 5′-monophosphate (pdT) and the base pairs pdA–pdT, pdG–pdC, and pdG(syn)–pdC are presented and discussed. The pm3 method is the first of the parameterized nddo quantum mechanical models with any ability to reproduce hydrogen bonding between nucleotide base pairs. Intermolecular hydrogen bond lengths between nucleotides displaying Watson–Crick base pairing are 0.1–0.2 Å less than experimental results. Nucleotide bond distances, bond angles, and torsion angles about the glycosyl bond (χ), the C4′C5′ bond (γ), and the C5′O5′ bond (β) agree with experimental results. There are many possible conformations of nucleotides. pm3 calculations reveal that many of the most stable conformations are stabilized by intramolecular CHO hydrogen bonds. These interactions disrupt the usual sugar puckering. The stacking interactions of a dT–pdA duplex are examined at different levels of gradient optimization. The intramolecular hydrogen bonds found in the nucleotide base pairs disappear in the duplex, as a result of the additional constraints on the phosphate group when part of a DNA backbone. Sugar puckering is reproduced by the pm3 method for the four bases in the dT–pdA duplex. pm3 underestimates the attractive stacking interactions of base pairs in a B-DNA helical conformation. The performance of the pm3 method implemented in SPARTAN is contrasted with that implemented in MOPAC. At present, accurate ab initio calculations are too timeconsuming to be of practical use, and molecular mechanics methods cannot be used to determine quantum mechanical properties such as reaction-path calculations, transition-state structures, and activation energies. The pm3 method should be used with extreme caution for examination of small DNA systems. Future parameterizations of semiempirical methods should incorporate base stacking interactions into the parameterization data set to enhance the ability of these methods.

Geophysical Analysis of Seasonal Montandon Gravel Ridge Water Table Fluctuation and Moisture Gradient Variation due to Storm Events

The purpose of this research project is to continue exploring the Montandon Long-Term Hydrologic Research Site(LTHR) by using multiple geophysical methods to obtain more accurate and precise information regarding subsurface hydrologic properties of a local gravel ridge,which are important to both the health of surrounding ecosystems and local agriculture. Through using non-invasive geophysical methods such as seismic refraction, Direct Current resistivity and ground penetrating radar (GPR) instead of invasive methods such as boreholedrilling which displace sediment and may alter water flow, data collection is less likely to bias the data itself. In addition to imaging the gravel ridge subsurface, another important researchpurpose is to observe how both water table elevation and the moisture gradient (moisture content of the unsaturated zone) change over a seasonal time period and directly after storm events. The combination of three types of data collection allows the strengths of each method combine together and provide a relatively strongly supported conclusions compared to previous research. Precipitation and geophysical data suggest that an overall increase in precipitation during the summer months causes a sharp decrease in subsurface resistivity within the unsaturated zone. GPR velocity data indicate significant immediate increase in moisture content within the shallow vadose zone (< 1m), suggesting that rain water was infiltrating into the shallow subsurface. Furthermore, the combination of resistivity and GPR results suggest that the decreased resistivity within the shallow layers is due to increased ion content within groundwater. This is unexpected as rainwater is assumed to have a DC resistivity value of 3.33*105 ohm-m. These results may suggest that ions within the sediment must beincorporated into the infiltrating water.

Membrane microdialysis: Evaluation of a new method to assess splanchnic tissue metabolism

OBJECTIVE: Measuring peritoneal lactate concentrations could be useful for detecting splanchnic hypoperfusion. The aims of this study were to evaluate the properties of a new membrane-based microdialyzer in vitro and to assess the ability of the dialyzer to detect a clinically relevant decrease in splanchnic blood flow in vivo. DESIGN: A membrane-based microdialyzer was first validated in vitro. The same device was tested afterward in a randomized, controlled animal experiment. SETTING: University experimental research laboratory. SUBJECTS: Twenty-four Landrace pigs of both genders. INTERVENTIONS: In vitro: Membrane microdialyzers were kept in warmed sodium lactate baths with lactate concentrations between 2 and 8 mmol/L for 10-120 mins, and microdialysis lactate concentrations were measured repeatedly (210 measurements). In vivo: An extracorporeal shunt with blood reservoir and roller pump was inserted between the proximal and distal abdominal aorta, and a microdialyzer was inserted intraperitoneally. In 12 animals, total splanchnic blood flow (measured by transit time ultrasound) was reduced by a median 43% (range, 13% to 72%) by activating the shunt; 12 animals served as controls. MEASUREMENTS AND MAIN RESULTS: In vitro: The fractional lactate recovery was 0.59 (0.32-0.83) after 60 mins and 0.82 (0.71-0.87) after 90 mins, with no further increase thereafter. At 60 and 90 mins, the fractional recovery was independent of the lactate concentration. In vivo: Abdominal blood flow reduction resulted in an increase in peritoneal microdialysis lactate concentration from 1.7 (0.3-3.8) mmol/L to 2.8 (1.3-6.2) mmol/L (p = .006). At the same time, mesenteric venous-arterial lactate gradient increased from 0.1 (-0.2-0.8) mmol/L to 0.3 (-0.3 -1.8) mmol/L (p = .032), and mesenteric venous-arterial Pco2 gradients increased from 12 (8-19) torr to 21 (11-54) torr (p = .005). CONCLUSIONS: Peritoneal membrane microdialysis provides a method for the assessment of splanchnic ischemia, with potential for clinical application.

Hardware accelerated ray cast of volume data and volume gradient for an optimized splines-based multi-resolution 2D-3D registration

This paper describes a method for DRR generation as well as for volume gradients projection using hardware accelerated 2D texture mapping and accumulation buffering and demonstrates its application in 2D-3D registration of X-ray fluoroscopy to CT images. The robustness of the present registration scheme are guaranteed by taking advantage of a coarse-to-fine processing of the volume/image pyramids based on cubic B-splines. A human cadaveric spine specimen together with its ground truth was used to compare the present scheme with a purely software-based scheme in three aspects: accuracy, speed, and capture ranges. Our experiments revealed an equivalent accuracy and capture ranges but with much shorter registration time with the present scheme. More specifically, the results showed 0.8 mm average target registration error, 55 second average execution time per registration, and 10 mm and 10° capture ranges for the present scheme when tested on a 3.0 GHz Pentium 4 computer.

Patterned-ground plant communities along a bioclimate gradient in the High Arctic, Canada

Non-sorted circles, non-sorted polygons, and earth hummocks are common ground-surface features ill arctic regions. The), are caused by a variety of physical processes that Occur in permafrost regions including contraction cracking and frost heave. Here we describe the vegetation of patterned-ground forms on zonal sites at three location!: along an N-S transect through the High Arctic of Canada. We made 75 releves on patterned-ground features (circles, polygons, earth hummocks) and adjacent tundra (Interpolygon, intercircle, interhummock areas) and identified and classified the vegetation according to the Braun-Blanquet Method. Environmental factors were correlated with the vegetation data using a nonmetric multidimensional scaling ordination (NMDS). We identified eleven commnunities: (1) Puccinellia angustata-Papaver radicalum community in xeromesic non-sorted polygons of subzone A of the Circumpolar Arctic Vegetation Map; (2) Saxifraga-Parmelia omphalodes ssp. glacialis community in hydromesic interpolygon areas of subzone A; (3) Hypogymnia subobscura-Lecanora epibryon community In xeromesic non-sorted polygons of subzone B; (4) Orthotrichum speciosum-Salix arctica community In xeromesic interpolygon areas of subzone B; (5) Cochlearia groenlandica-Luzula nivalis community in hydromesic earth Mocks Of subzone B; (6) Salix arctica-Eriophorum angustifolium ssp. triste community in hygric earth hummocks of subzone 13; (7) Puccinellia angustata-Potentilla vahliana community in xeromesic non-sorted circles and bare patches of subzone Q (8) Dryas integrifolia-Carex rupestris community in xeromesic intercircle areas and vegetated patches of subzone C; (9) Braya glabella ssp. purpurascens-Dryas integrifolia community In hydromesic non-sorted circles of subzone Q (10) Dryas integrifolia-Carex aquatilis community in hydromesic intercircle areas of subzone C; and (11) Eriophorum angustifolium ssp. triste-Carex aquatilis community ill hygric intercircle areas of subzone C. The NMDS ordination displayed the vegetation types with respect to complex environmental gradients. The first axis of the ordination corresponds to a complex soil moisture gradient and the second axis corresponds to a complex geology/elevation/climate gradient. The tundra plots have a greater moss and graminoid cover than the adjacent frost-heave communities. In general, frost-heave features have greater thaw depths, more bare ground, thinner organic horizons, and lower soil moisture than the surrounding tundra. The morphology of the investigated patterned ground forms changes along the climatic gradient, with non-sorted pollygons dominating in the northernmost sites and non-sorted circles dominating, in the southern sites.

Cosine Lobe Based Relighting from Gradient Illumination Photographs

We present an image-based method for relighting a scene by analytically fitting cosine lobes to the reflectance function at each pixel, based on gradient illumination photographs. Realistic relighting results for many materials are obtained using a single per-pixel cosine lobe obtained from just two color photographs: one under uniform white illumination and the other under colored gradient illumination. For materials with wavelength-dependent scattering, a better fit can be obtained using independent cosine lobes for the red, green, and blue channels, obtained from three achromatic gradient illumination conditions instead of the colored gradient condition. We explore two cosine lobe reflectance functions, both of which allow an analytic fit to the gradient conditions. One is non-zero over half the sphere of lighting directions, which works well for diffuse and specular materials, but fails for materials with broader scattering such as fur. The other is non-zero everywhere, which works well for broadly scattering materials and still produces visually plausible results for diffuse and specular materials. We also perform an approximate diffuse/specular separation of the reflectance, and estimate scene geometry from the recovered photometric normals to produce hard shadows cast by the geometry, while still reconstructing the input photographs exactly.

Robust Proximal Femur Segmentation in Conventional X-Ray Images via Random Forest Regression on Multi-resolution Gradient Features

In this paper, we propose a fully automatic, robust approach for segmenting proximal femur in conventional X-ray images. Our method is based on hierarchical landmark detection by random forest regression, where the detection results of 22 global landmarks are used to do the spatial normalization, and the detection results of the 59 local landmarks serve as the image cue for instantiation of a statistical shape model of the proximal femur. To detect landmarks in both levels, we use multi-resolution HoG (Histogram of Oriented Gradients) as features which can achieve better accuracy and robustness. The efficacy of the present method is demonstrated by experiments conducted on 150 clinical x-ray images. It was found that the present method could achieve an average point-to-curve error of 2.0 mm and that the present method was robust to low image contrast, noise and occlusions caused by implants.

Identification of the Causative Bacteria in Musculoskeletal Infections Using 16s rDNA - Denaturing Gradient Gel Electrophoresis Analysis

Musculoskeletal infections are infections of the bone and surrounding tissues. They are currently diagnosed based on culture analysis, which is the gold standard for pathogen identification. However, these clinical laboratory methods are frequently inadequate for the identification of the causative agents, because a large percentage (25-50%) of confirmed musculoskeletal infections are false negatives in which no pathogen is identified in culture. My data supports these results. The goal of this project was to use PCR amplification of a portion of the 16S rRNA gene to test an alternative approach for the identification of these pathogens and to assess the diversity of the bacteria involved. The advantages of this alternative method are that it should increase sample sensitivity and the speed of detection. In addition, bacteria that are non-culturable or in low abundance can be detected using this molecular technique. However, a complication of this approach is that the majority of musculoskeletal infections are polymicrobial, which prohibits direct identification from the infected tissue by DNA sequencing of the initial 16S rDNA amplification products. One way to solve this problem is to use denaturing gradient gel electrophoresis (DGGE) to separate the PCR products before DNA sequencing. Denaturing gradient gel electrophoresis (DGGE) separates DNA molecules based on their melting point, which is determined by their DNA sequence. This analytical technique allows a mixture of PCR products of the same length that electrophoreses through agarose gels as one band, to be separated into different bands and then used for DNA sequence analysis. In this way, the DGGE allows for the identification of individual bacterial species in polymicrobial-infected tissue, which is critical for improving clinical outcomes. By combining the 16S rDNA amplification and the DGGE techniques together, an alternative approach for identification has been used. The 16S rRNA gene PCR-DGGE method includes several critical steps: DNA extraction from tissue biopsies, amplification of the bacterial DNA, PCR product separation by DGGE, amplification of the gel-extracted DNA, and DNA sequencing and analysis. Each step of the method was optimized to increase its sensitivity and for rapid detection of the bacteria present in human tissue samples. The limit of detection for the DNA extraction from tissue was at least 20 Staphylococcus aureus cells and the limit of detection for PCR was at least 0.05 pg of template DNA. The conditions for DGGE electrophoreses were optimized by using a double gradient of acrylamide (6 – 10%) and denaturant (30-70%), which increased the separation between distinct PCR products. The use of GelRed (Biotium) improved the DNA visualization in the DGGE gel. To recover the DNA from the DGGE gels the gel slices were excised, shredded in a bead beater, and the DNA was allowed to diffuse into sterile water overnight. The use of primers containing specific linkers allowed the entire amplified PCR product to be sequenced and then analyzed. The optimized 16S rRNA gene PCR-DGGE method was used to analyze 50 tissue biopsy samples chosen randomly from our collection. The results were compared to those of the Memorial Hermann Hospital Clinical Microbiology Laboratory for the same samples. The molecular method was congruent for 10 of the 17 (59%) culture negative tissue samples. In 7 of the 17 (41%) culture negative the molecular method identified a bacterium. The molecular method was congruent with the culture identification for 7 of the 33 (21%) positive cultured tissue samples. However, in 8 of the 33 (24%) the molecular method identified more organisms. In 13 of the 15 (87%) polymicrobial cultured tissue samples the molecular method identified at least one organism that was also identified by culture techniques. Overall, the DGGE analysis of 16S rDNA is an effective method to identify bacteria not identified by culture analysis.

Development of a rapid column-switching LC-MS/MS method for the quantification of THCCOOH and THCCOOH-glucuronide in whole blood for assessing cannabis consumption frequency

The concentration of 11-nor-9-carboxy-Δ(9)-tetrahydrocannabinol (THCCOOH) in whole blood is used as a parameter for assessing the consumption behavior of cannabis consumers. The blood level of THCCOOH-glucuronide might provide additional information about the frequency of cannabis use. To verify this assumption, a column-switching liquid chromatography-tandem mass spectrometry (LC-MS/MS) method for the rapid and direct quantification of free and glucuronidated THCCOOH in human whole blood was newly developed. The method comprised protein precipitation, followed by injection of the processed sample onto a trapping column and subsequent gradient elution to an analytical column for separation and detection. The total LC run time was 4.5 min. Detection of the analytes was accomplished by electrospray ionization in positive ion mode and selected reaction monitoring using a triple-stage quadrupole mass spectrometer. The method was fully validated by evaluating the following parameters: linearity, lower limit of quantification, accuracy and imprecision, selectivity, extraction efficiency, matrix effect, carry-over, dilution integrity, analyte stability, and re-injection reproducibility. All acceptance criteria were analyzed and the predefined criteria met. Linearity ranged from 5.0 to 500 μg/L for both analytes. The method was successfully applied to whole blood samples from a large collective of cannabis consumers, demonstrating its applicability in the forensic field.

Multiple gradient echo propeller (MGREP) MRI: Technical development and potential applications

The PROPELLER (Periodically Rotated Overlapping Parallel Lines with Enhanced Reconstruction) magnetic resonance imaging (MRI) technique has inherent advantages over other fast imaging methods, including robust motion correction, reduced image distortion, and resistance to off-resonance effects. These features make PROPELLER highly desirable for T2*-sensitive imaging, high-resolution diffusion imaging, and many other applications. However, PROPELLER has been predominantly implemented as a fast spin-echo (FSE) technique, which is insensitive to T2* contrast, and requires time-inefficient signal averaging to achieve adequate signal-to-noise ratio (SNR) for many applications. These issues presently constrain the potential clinical utility of FSE-based PROPELLER. ^ In this research, our aim was to extend and enhance the potential applications of PROPELLER MRI by developing a novel multiple gradient echo PROPELLER (MGREP) technique that can overcome the aforementioned limitations. The MGREP pulse sequence was designed to acquire multiple gradient-echo images simultaneously, without any increase in total scan time or RF energy deposition relative to FSE-based PROPELLER. A new parameter was also introduced for direct user-control over gradient echo spacing, to allow variable sensitivity to T2* contrast. In parallel to pulse sequence development, an improved algorithm for motion correction was also developed and evaluated against the established method through extensive simulations. The potential advantages of MGREP over FSE-based PROPELLER were illustrated via three specific applications: (1) quantitative T2* measurement, (2) time-efficient signal averaging, and (3) high-resolution diffusion imaging. Relative to the FSE-PROPELLER method, the MGREP sequence was found to yield quantitative T2* values, increase SNR by ∼40% without any increase in acquisition time or RF energy deposition, and noticeably improve image quality in high-resolution diffusion maps. In addition, the new motion algorithm was found to improve the performance considerably in motion-artifact reduction. ^ Overall, this work demonstrated a number of enhancements and extensions to existing PROPELLER techniques. The new technical capabilities of PROPELLER imaging, developed in this thesis research, are expected to serve as the foundation for further expanding the scope of PROPELLER applications. ^

Hybrid topological derivative and gradient-based methods for electrical impedance tomography

We present a technique to reconstruct the electromagnetic properties of a medium or a set of objects buried inside it from boundary measurements when applying electric currents through a set of electrodes. The electromagnetic parameters may be recovered by means of a gradient method without a priori information on the background. The shape, location and size of objects, when present, are determined by a topological derivative-based iterative procedure. The combination of both strategies allows improved reconstructions of the objects and their properties, assuming a known background.