Subjectively adapted high capacity lossless image data hiding based on prediction errors

This article reports on a lossless data hiding scheme for digital images where the data hiding capacity is either determined by minimum acceptable subjective quality or by the demanded capacity. In the proposed method data is hidden within the image prediction errors, where the most well-known prediction algorithms such as the median edge detector (MED), gradient adjacent prediction (GAP) and Jiang prediction are tested for this purpose. In this method, first the histogram of the prediction errors of images are computed and then based on the required capacity or desired image quality, the prediction error values of frequencies larger than this capacity are shifted. The empty space created by such a shift is used for embedding the data. Experimental results show distinct superiority of the image prediction error histogram over the conventional image histogram itself, due to much narrower spectrum of the former over the latter. We have also devised an adaptive method for hiding data, where subjective quality is traded for data hiding capacity. Here the positive and negative error values are chosen such that the sum of their frequencies on the histogram is just above the given capacity or above a certain quality.

Reversible and high capacity data hiding in medical images

In this paper we introduce a highly efficient reversible data hiding system. It is based on dividing the image into tiles and shifting the histograms of each image tile between its minimum and maximum frequency. Data are then inserted at the pixel level with the largest frequency to maximize data hiding capacity. It exploits the special properties of medical images, where the histogram of their nonoverlapping image tiles mostly peak around some gray values and the rest of the spectrum is mainlyempty. The zeros (or minima) and peaks (maxima) of the histograms of the image tiles are then relocated to embed the data. The grey values of some pixels are therefore modified.High capacity, high fidelity, reversibility and multiple data insertions are the key requirements of data hiding in medical images. We show how histograms of image tiles of medical images can be exploited to achieve these requirements. Compared with data hiding method applied to the whole image, our scheme can result in 30%-200% capacity improvement and still with better image quality, depending on the medical image content. Additional advantages of the proposed method include hiding data in the regions of non-interest and better exploitation of spatial masking.

Mating system parameters in a high density population of andirobas in the Amazon forest

The objective of this work was to estimate the mating system parameters of a andiroba (Carapa guianensis) population using microsatellite markers and the mixed and correlated mating models. Twelve open‑pollinated progeny arrays of 15 individuals were sampled in an area with C. guianensis estimated density of 25.7 trees per hectare. Overall, the species has a mixed reproductive system, with a predominance of outcrossing. The multilocus outcrossing rate (t m = 0.862) was significantly lower than the unity, indicating that self‑pollination occurred. The rate of biparental inbreeding was substantial (t m ‑ t s = 0.134) and significantly different from zero. The correlation of selfing within progenies was high (r s = 0.635), indicating variation in the individual outcrossing rate. Consistent with this result, the estimate of the individual outcrossing rate ranged from 0.598 to 0.978. The multilocus correlation of paternity was low (r p(m) = 0.081), but significantly different from zero, suggesting that the progenies contain full‑sibs. The coancestry within progenies (Θ = 0.185) was higher and the variance effective size (Ne(v) = 2.7) was lower than expected for true half‑sib progenies (Θ = 0.125; Ne(v) = 4). These results suggest that, in order to maintain a minimum effective size of 150 individuals for breeding, genetic conservation, and environmental reforestation programs, seeds from at least 56 trees must be collected.

Carbon-isotope stratigraphy and ammonite faunal turnover for the Middle Jurassic in the Southern Iberian palaeomargin

Variations in the stable carbon-isotope ratio of marine and continental sediments can reflect changes in sink and flux modifications of the palaeocarbon cycle. Here we report carbon-isotope compositions of Middle Jurassic marine carbonates from the Betic Cordillera (southern Spain), which represents an ideal region to link the stable carbon-isotope curves directly to ammonite zones and subzones, and thereby for the first time achieve an accurate chronostratigraphic calibration. The five sections studied represent basin and high swell deposits of the Southern Iberian palaeomargin. We find a similar delta C-13 of carbonates between different oceanic areas, suggesting a homogeneous carbon-isotope oceanic reservoir through the Middle Jurassic. The Aalenian-Bajocian transition is a critical period in ammonite evolution; hence the Early Jurassic fauna are replaced by new ammonite families which become dominant throughout the Middle and Late Jurassic. For this reason, we compared the delta C-13 values of carbonates with ammonite diversity and extinction rates at different taxonomical levels in order to explore the possible relationship between the carbon cycle and ammonite evolution. The carbon-isotope values of carbonates are not exactly linearly correlated with the extinction rate and ammonite diversity, but the main faunal turnovers follow minimum delta C-13 values, where extinct taxa are replaced by new ones. Likewise, radiation episodes are associated with increasing delta C-13 values and with transgressive sea-level rise. All these data support the idea that perturbations in the global carbon cycle reflect rapid palaeoenvironmental changes. We made detailed analyses of these faunal turnovers, using them as a proxy to identify major palaeoenvironmental crises in their ecosystems forced by modification in the carbon cycle. (c) 2006 Elsevier B.V All rights reserved.

A novel semi-fragile forensic watermarking scheme for remote sensing images

Peer-reviewed

High capacity FFT-based audio watermarking

Peer-reviewed

A robust audio watermarking scheme based on MPEG 1 Layer 3 Compression

This paper describes an audio watermarking scheme based on lossy compression. The main idea is taken from an image watermarking approach where the JPEG compression algorithm is used to determine where and how the mark should be placed. Similarly, in the audio scheme suggested in this paper, an MPEG 1 Layer 3 algorithm is chosen for compression to determine the position of the mark bits and, thus, the psychoacoustic masking of the MPEG 1 Layer 3compression is implicitly used. This methodology provides with a high robustness degree against compression attacks. The suggested scheme is also shown to succeed against most of the StirMark benchmark attacks for audio.

The Gaussian assumption in second-order estimation problems in digital communications

This paper deals with the goodness of the Gaussian assumption when designing second-order blind estimationmethods in the context of digital communications. The low- andhigh-signal-to-noise ratio (SNR) asymptotic performance of the maximum likelihood estimator—derived assuming Gaussiantransmitted symbols—is compared with the performance of the optimal second-order estimator, which exploits the actualdistribution of the discrete constellation. The asymptotic study concludes that the Gaussian assumption leads to the optimalsecond-order solution if the SNR is very low or if the symbols belong to a multilevel constellation such as quadrature-amplitudemodulation (QAM) or amplitude-phase-shift keying (APSK). On the other hand, the Gaussian assumption can yield importantlosses at high SNR if the transmitted symbols are drawn from a constant modulus constellation such as phase-shift keying (PSK)or continuous-phase modulations (CPM). These conclusions are illustrated for the problem of direction-of-arrival (DOA) estimation of multiple digitally-modulated signals.

Quantitative modelling of transcription factor binding sites

Abstract One of the most important issues in molecular biology is to understand regulatory mechanisms that control gene expression. Gene expression is often regulated by proteins, called transcription factors which bind to short (5 to 20 base pairs),degenerate segments of DNA. Experimental efforts towards understanding the sequence specificity of transcription factors is laborious and expensive, but can be substantially accelerated with the use of computational predictions. This thesis describes the use of algorithms and resources for transcriptionfactor binding site analysis in addressing quantitative modelling, where probabilitic models are built to represent binding properties of a transcription factor and can be used to find new functional binding sites in genomes. Initially, an open-access database(HTPSELEX) was created, holding high quality binding sequences for two eukaryotic families of transcription factors namely CTF/NF1 and LEFT/TCF. The binding sequences were elucidated using a recently described experimental procedure called HTP-SELEX, that allows generation of large number (> 1000) of binding sites using mass sequencing technology. For each HTP-SELEX experiments we also provide accurate primary experimental information about the protein material used, details of the wet lab protocol, an archive of sequencing trace files, and assembled clone sequences of binding sequences. The database also offers reasonably large SELEX libraries obtained with conventional low-throughput protocols.The database is available at http://wwwisrec.isb-sib.ch/htpselex/ and and ftp://ftp.isrec.isb-sib.ch/pub/databases/htpselex. The Expectation-Maximisation(EM) algorithm is one the frequently used methods to estimate probabilistic models to represent the sequence specificity of transcription factors. We present computer simulations in order to estimate the precision of EM estimated models as a function of data set parameters(like length of initial sequences, number of initial sequences, percentage of nonbinding sequences). We observed a remarkable robustness of the EM algorithm with regard to length of training sequences and the degree of contamination. The HTPSELEX database and the benchmarked results of the EM algorithm formed part of the foundation for the subsequent project, where a statistical framework called hidden Markov model has been developed to represent sequence specificity of the transcription factors CTF/NF1 and LEF1/TCF using the HTP-SELEX experiment data. The hidden Markov model framework is capable of both predicting and classifying CTF/NF1 and LEF1/TCF binding sites. A covariance analysis of the binding sites revealed non-independent base preferences at different nucleotide positions, providing insight into the binding mechanism. We next tested the LEF1/TCF model by computing binding scores for a set of LEF1/TCF binding sequences for which relative affinities were determined experimentally using non-linear regression. The predicted and experimentally determined binding affinities were in good correlation.

Second-order parameter estimation

This work provides a general framework for the design of second-order blind estimators without adopting anyapproximation about the observation statistics or the a prioridistribution of the parameters. The proposed solution is obtainedminimizing the estimator variance subject to some constraints onthe estimator bias. The resulting optimal estimator is found todepend on the observation fourth-order moments that can be calculatedanalytically from the known signal model. Unfortunately,in most cases, the performance of this estimator is severely limitedby the residual bias inherent to nonlinear estimation problems.To overcome this limitation, the second-order minimum varianceunbiased estimator is deduced from the general solution by assumingaccurate prior information on the vector of parameters.This small-error approximation is adopted to design iterativeestimators or trackers. It is shown that the associated varianceconstitutes the lower bound for the variance of any unbiasedestimator based on the sample covariance matrix.The paper formulation is then applied to track the angle-of-arrival(AoA) of multiple digitally-modulated sources by means ofa uniform linear array. The optimal second-order tracker is comparedwith the classical maximum likelihood (ML) blind methodsthat are shown to be quadratic in the observed data as well. Simulationshave confirmed that the discrete nature of the transmittedsymbols can be exploited to improve considerably the discriminationof near sources in medium-to-high SNR scenarios.

Computational and Experimental Analysis of Flow Field in the Diffusers of Centrifugal Compressors

Centrifugal compressors are widely used for example in process industry, oil and gas industry, in small gas turbines and turbochargers. In order to achieve lower consumption of energy and operation costs the efficiency of the compressor needs to be improve. In the present work different pinches and low solidity vaned diffusers were utilized in order to improve the efficiency of a medium size centrifugal compressor. In this study, pinch means the decrement of the diffuser flow passage height. First different geometries were analyzed using computational fluid dynamics. The flow solver Finflo was used to solve the flow field. Finflo is a Navier-Stokes solver. The solver is capable to solve compressible, incompressible, steady and unsteady flow fields. Chien's k-e turbulence model was used. One of the numerically investigated pinched diffuser and one low solidity vaned diffuser were studied experimentally. The overall performance of the compressor and the static pressure distribution before and after the diffuser were measured. The flow entering and leaving the diffuser was measured using a three-hole Cobra-probe and Kiel-probes. The pinch and the low solidity vaned diffuser increased the efficiency of the compressor. Highest isentropic efficiency increment obtained was 3\% of the design isentropic efficiency of the original geometry. It was noticed in the numerical results that the pinch made to the hub and the shroud wall was most beneficial to the operation of the compressor. Also the pinch made to the hub was better than the pinchmade to the shroud. The pinch did not affect the operation range of the compressor, but the low solidity vaned diffuser slightly decreased the operation range.The unsteady phenomena in the vaneless diffuser were studied experimentally andnumerically. The unsteady static pressure was measured at the diffuser inlet and outlet, and time-accurate numerical simulation was conducted. The unsteady static pressure showed that most of the pressure variations lay at the passing frequency of every second blade. The pressure variations did not vanish in the diffuser and were visible at the diffuser outlet. However, the amplitude of the pressure variations decreased in the diffuser. The time-accurate calculations showed quite a good agreement with the measured data. Agreement was very good at the design operation point, even though the computational grid was not dense enough inthe volute and in the exit cone. The time-accurate calculation over-predicted the amplitude of the pressure variations at high flow.

Ground-state order and spin-lattice coupling in tetrahedral spin systems Cu2Te2O5X2 (X = Br or Cl)

High-resolution ac susceptibility and thermal conductivity measurement on Cu2Te2O5X2 (X=Br,Cl) single crystals are reported. For Br-sample, sample dependence prevents one from distinguishing between possibilities of magnetically ordered and spin-singlet ground states. In Cl-sample a three-dimensional transition at 18.5 K is accompanied by almost isotropic behavior of susceptibility and almost switching behavior of thermal conductivity. Thermal conductivity studies suggest the presence of a tremendous spin-lattice coupling characterizing Cl- but not Br-sample. Below the transition Cl-sample is in a complex magnetic state involving AF order but also the elements consistent with the presence of a gap in the excitation spectrum.

Comparison of plantar pressure distribution in adolescent runners at low vs. high running velocity

This study aimed to compare foot plantar pressure distribution while jogging and running in highly trained adolescent runners. Eleven participants performed two constant-velocity running trials either at jogging (11.2 ± 0.9 km/h) or running (17.8 ± 1.4 km/h) pace on a treadmill. Contact area (CA in cm(2)), maximum force (F(max) in N), peak pressure (PP in kPa), contact time (CT in ms), and relative load (force time integral in each individual region divided by the force time integral for the total plantar foot surface, in %) were measured in nine regions of the right foot using an in-shoe plantar pressure device. Under the whole foot, CA, F(max) and PP were lower in jogging than in running (-1.2% [p<0.05], -12.3% [p<0.001] and -15.1% [p<0.01] respectively) whereas CT was higher (+20.1%; p<0.001). Interestingly, we found an increase in relative load under the medial and central forefoot regions while jogging (+6.7% and +3.7%, respectively; [p<0.05]), while the relative load under the lesser toes (-8.4%; p<0.05) was reduced. In order to prevent overloading of the metatarsals in adolescent runners, excessive mileage at jogging pace should be avoided.

Streptabody, a high avidity molecule made by tetramerization of in vivo biotinylated, phage display-selected scFv fragments on streptavidin.

Phage display is a powerful method of isolating of antibody fragments from highly diverse naive human antibody repertoires. However, the affinity of the selected antibodies is usually low and current methods of affinity maturation are complex and time-consuming. In this paper, we describe an easy way to increase the functional affinity (avidity) of single chain variable fragments (scFvs) by tetramerization on streptavidin, following their site-specific biotinylation by the enzyme BirA. Expression vectors have been constructed that enable addition of the 15 amino acid biotin acceptor domain (BAD) on selected scFvs. Different domains were cloned at the C-terminus of scFv in the following order: a semi-rigid hinge region (of 16 residues), the BAD, and a histidine tail. Two such recombinant scFvs directed against the carcinoembryonic antigen (CEA) were previously selected from human non-immune and murine immune phage display libraries. The scFvs were first synthesized in Escherichia coli carrying the plasmid encoding the BirA enzyme, and then purified from the cytoplasmic extracts by Ni-NTA affinity chromatography. Purified biotinylated scFvs were tetramerized on the streptavidin molecule to create a streptabody (StAb). The avidity of various forms of anti-CEA StAbs, tested on purified CEA by competitive assays and surface plasmon resonance showed an increase of more than one log, as compared with the scFv monomer counterparts. Furthermore, the percentage of direct binding of 125I-labeled StAb or monomeric scFv on CEA-Sepharose beads and on CEA-expressing cells showed a dramatic increase for the tetramerized scFv (>80%), as compared with the monomeric scFv (<20%). Interestingly, the percentage binding of 125I-labeled anti-CEA StAbs to CEA-expressing colon carcinoma cells was definitely higher (>80%) than that obtained with a reference high affinity murine anti-CEA mAb (30%). Another advantage of using scFvs in a StAb format was demonstrated by Western blot analysis, where tetramerized anti-CEA scFv could detect a small quantity of CEA at a concentration 100-fold lower than the monomeric scFv.

Buprenorphine and norbuprenorphine quantification in human plasma by simple protein precipitation and ultra-high performance liquid chromatography tandem mass spectrometry.

A highly sensitive ultra-high performance liquid chromatography tandem mass spectrometry (UHPLC-MS/MS) method was developed for the quantification of buprenorphine and its major metabolite norbuprenorphine in human plasma. In order to speed up the process and decrease costs, sample preparation was performed by simple protein precipitation with acetonitrile. To the best of our knowledge, this is the first application of this extraction technique for the quantification of buprenorphine in plasma. Matrix effects were strongly reduced and selectivity increased by using an efficient chromatographic separation on a sub-2μm column (Acquity UPLC BEH C18 1.7μm, 2.1×50mm) in 5min with a gradient of ammonium formate 20mM pH 3.05 and acetonitrile as mobile phase at a flow rate of 0.4ml/min. Detection was made using a tandem quadrupole mass spectrometer operating in positive electrospray ionization mode, using multiple reaction monitoring. The procedure was fully validated according to the latest Food and Drug Administration guidelines and the Société Française des Sciences et Techniques Pharmaceutiques. Very good results were obtained by using a stable isotope-labeled internal standard for each analyte, to compensate for the variability due to the extraction and ionization steps. The method was very sensitive with lower limits of quantification of 0.1ng/ml for buprenorphine and 0.25ng/ml for norbuprenorphine. The upper limit of quantification was 250ng/ml for both drugs. Trueness (98.4-113.7%), repeatability (1.9-7.7%), intermediate precision (2.6-7.9%) and internal standard-normalized matrix effects (94-101%) were in accordance with international recommendations. The procedure was successfully used to quantify plasma samples from patients included in a clinical pharmacogenetic study and can be transferred for routine therapeutic drug monitoring in clinical laboratories without further development.