Burst error correction using partial fourier matrices and block sparse representation

There is a strong relation between sparse signal recovery and error control coding. It is known that burst errors are block sparse in nature. So, here we attempt to solve burst error correction problem using block sparse signal recovery methods. We construct partial Fourier based encoding and decoding matrices using results on difference sets. These constructions offer guaranteed and efficient error correction when used in conjunction with reconstruction algorithms which exploit block sparsity.

A MatLAB Based Virtual Phantom for 2D Electrical Impedance Tomography (MatVP2DEIT): Studying the Medical Electrical Impedance Tomography Reconstruction in Computer

Electrical Impedance Tomography (EIT) is a computerized medical imaging technique which reconstructs the electrical impedance images of a domain under test from the boundary voltage-current data measured by an EIT electronic instrumentation using an image reconstruction algorithm. Being a computed tomography technique, EIT injects a constant current to the patient's body through the surface electrodes surrounding the domain to be imaged (Omega) and tries to calculate the spatial distribution of electrical conductivity or resistivity of the closed conducting domain using the potentials developed at the domain boundary (partial derivative Omega). Practical phantoms are essentially required to study, test and calibrate a medical EIT system for certifying the system before applying it on patients for diagnostic imaging. Therefore, the EIT phantoms are essentially required to generate boundary data for studying and assessing the instrumentation and inverse solvers a in EIT. For proper assessment of an inverse solver of a 2D EIT system, a perfect 2D practical phantom is required. As the practical phantoms are the assemblies of the objects with 3D geometries, the developing of a practical 2D-phantom is a great challenge and therefore, the boundary data generated from the practical phantoms with 3D geometry are found inappropriate for assessing a 2D inverse solver. Furthermore, the boundary data errors contributed by the instrumentation are also difficult to separate from the errors developed by the 3D phantoms. Hence, the errorless boundary data are found essential to assess the inverse solver in 2D EIT. In this direction, a MatLAB-based Virtual Phantom for 2D EIT (MatVP2DEIT) is developed to generate accurate boundary data for assessing the 2D-EIT inverse solvers and the image reconstruction accuracy. MatVP2DEIT is a MatLAB-based computer program which simulates a phantom in computer and generates the boundary potential data as the outputs by using the combinations of different phantom parameters as the inputs to the program. Phantom diameter, inhomogeneity geometry (shape, size and position), number of inhomogeneities, applied current magnitude, background resistivity, inhomogeneity resistivity all are set as the phantom variables which are provided as the input parameters to the MatVP2DEIT for simulating different phantom configurations. A constant current injection is simulated at the phantom boundary with different current injection protocols and boundary potential data are calculated. Boundary data sets are generated with different phantom configurations obtained with the different combinations of the phantom variables and the resistivity images are reconstructed using EIDORS. Boundary data of the virtual phantoms, containing inhomogeneities with complex geometries, are also generated for different current injection patterns using MatVP2DEIT and the resistivity imaging is studied. The effect of regularization method on the image reconstruction is also studied with the data generated by MatVP2DEIT. Resistivity images are evaluated by studying the resistivity parameters and contrast parameters estimated from the elemental resistivity profiles of the reconstructed phantom domain. Results show that the MatVP2DEIT generates accurate boundary data for different types of single or multiple objects which are efficient and accurate enough to reconstruct the resistivity images in EIDORS. The spatial resolution studies show that, the resistivity imaging conducted with the boundary data generated by MatVP2DEIT with 2048 elements, can reconstruct two circular inhomogeneities placed with a minimum distance (boundary to boundary) of 2 mm. It is also observed that, in MatVP2DEIT with 2048 elements, the boundary data generated for a phantom with a circular inhomogeneity of a diameter less than 7% of that of the phantom domain can produce resistivity images in EIDORS with a 1968 element mesh. Results also show that the MatVP2DEIT accurately generates the boundary data for neighbouring, opposite reference and trigonometric current patterns which are very suitable for resistivity reconstruction studies. MatVP2DEIT generated data are also found suitable for studying the effect of the different regularization methods on reconstruction process. Comparing the reconstructed image with an original geometry made in MatVP2DEIT, it would be easier to study the resistivity imaging procedures as well as the inverse solver performance. Using the proposed MatVP2DEIT software with modified domains, the cross sectional anatomy of a number of body parts can be simulated in PC and the impedance image reconstruction of human anatomy can be studied.

Network reconstruction using knock-out and over-expression data

This paper outlines necessary and sufficient conditions for network reconstruction of linear, time-invariant systems using data from either knock-out or over-expression experiments. These structural system perturbations, which are common in biological experiments, can be formulated as unknown system inputs, allowing the network topology and dynamics to be found. We assume that only partial state measurements are available and propose an algorithm that can reconstruct the network at the level of the measured states using either time-series or steady-state data. A simulated example illustrates how the algorithm successfully reconstructs a network from data. © 2013 EUCA.

Phase reconstruction from three interferograms based on integral of phase gradient

A novel algorithm of phase reconstruction based on the integral of phase gradient is presented. The algorithm directly derives two real-valued partial derivatives from three phase-shifted interferograms. Through integrating the phase derivatives, the desired phase is reconstructed. During the phase reconstruction process, there is no need for an extra rewrapping manipulation to ensure values of the phase derivatives lie in the interval [-pi, pi] as before, thus this algorithm can prevent error or distortion brought about by the phase unwrapping operation. Additionally, this algorithm is fast and easy to implement, and insensitive to the nonuniformity of the intensity distribution of the interferogram. The feasibility of the algorithm is demonstrated by both computer simulation and experiment.

ENERGY-LEVEL CALCULATIONS OF THE RECONSTRUCTED 90-DEGREES PARTIAL DISLOCATION IN SILICON

An LCAO scheme taking into account 10 atomic orbitals (s-, p-, and d-type) applied to a supercell containing 256 atoms is used to calculate the bound states of the reconstructed 90-degrees partial dislocation in Si. The results differ significantly from our earlier calculations on the unreconstructed 90-degrees partial using the same method. We find two bands separate from each other in the entire Brillouin zone and the upper band penetrates deep into the indirect band gap which is in contradiction with the general opinion that core reconstruction clears the band gap of dislocation states.

The local and global dimensions of metalliferous pollution derived from a reconstruction of an eight thousand year record of copper smelting and mining at a desert-mountain frontier in southern Jordan

Grattan, J.P., Gilbertson, D.D., Hunt, C.O. (2007). The local and global dimensions of metaliferrous air pollution derived from a reconstruction of an 8 thousand year record of copper smelting and mining at a desert-mountain frontier in southern Jordan. Journal of Archaeological Science 34, 83-110

Analysis of Extended Partial Least Squares for Monitoring Large-Scale Processes

This paper theoretically analysis the recently proposed "Extended Partial Least Squares" (EPLS) algorithm. After pointing out some conceptual deficiencies, a revised algorithm is introduced that covers the middle ground between Partial Least Squares and Principal Component Analysis. It maximises a covariance criterion between a cause and an effect variable set (partial least squares) and allows a complete reconstruction of the recorded data (principal component analysis). The new and conceptually simpler EPLS algorithm has successfully been applied in detecting and diagnosing various fault conditions, where the original EPLS algorithm did only offer fault detection.

Creating replicas of scultural heritage through 3D reconstruction and low cost 3D printer in Education

The process of making replicas of heritage has traditionally been developed by public agencies, corporations and museums and is not commonly used in schools. Currently there are technologies that allow creating cheap replicas. The new 3D reconstruction software, based on photographs and low cost 3D printers allow to make replicas at a cost much lower than traditional. This article describes the process of creating replicas of the sculpture Goslar Warrior of artist Henry Moore, located in Santa Cruz de Tenerife. To make this process, first, a digital model have been created using Autodesk Recap 360, Autodesk 123D Catch and Autodesk Meshmixer MarkerBot MakerWare applications. Physical replication, has been reproduced in polylactic acid (PLA) by MakerBot Replicator 2 3D printer. In addition, a cost analysis using, in one hand, the printer mentioned, and in the other hand, 3D printing services both online and local, is included. Finally, there has been a specific action with 141 students and 12 high school teachers, who filled a questionnary about the use of sculptural replicas in education.

Extracting S-matrix poles for resonances from numerical scattering data: Type-II Pade reconstruction

We present a FORTRAN 77 code for evaluation of resonance pole positions and residues of a numerical scattering matrix element in the complex energy (CE) as well as in the complex angular momentum (CAM) planes. Analytical continuation of the S-matrix element is performed by constructing a type-II Pade approximant from given physical values (Bessis et al. (1994) [421: Vrinceanu et al. (2000) [24]; Sokolovski and Msezane (2004) [23]). The algorithm involves iterative 'preconditioning' of the numerical data by extracting its rapidly oscillating potential phase component. The code has the capability of adding non-analytical noise to the numerical data in order to select 'true' physical poles, investigate their stability and evaluate the accuracy of the reconstruction. It has an option of employing multiple-precision (MPFUN) package (Bailey (1993) [451) developed by D.H. Bailey wherever double precision calculations fail due to a large number of input partial waves (energies) involved. The code has been successfully tested on several models, as well as the F + H-2 -> HE + H, F + HD : HE + D, Cl + HCI CIH + Cl and H + D-2 -> HD + D reactions. Some detailed examples are given in the text.

Post mastectomy radiotherapy in breast cancer: a survey of current United Kingdom practice

Purpose: Radiotherapy (RT) is increasingly used following mastectomy for breast cancer While indications for post-mastectomy radiotherapy (PMRT) are clear in patient groups at high risk of local recurrence, guidelines are less clear in intermediate-risk patients and patients with ductal carcinoma in situ (DCIS). This study aimed to determine variations in the use of PMRT in the United Kingdom (UK).

Methods: A postal survey of all consultant breast surgeon members of the Association of Breast Surgery in the UK.

Results: Tumour size and nodal status were confirmed as the most important indications for PMRT There was significant variation in the influence of other factors such as tumour grade, lymphovascular invasion and margin status. Nineteen per cent of respondents stated that they would consider the use of PMRT in cases of DCIS alone.

Conclusions: There is significant variation in practice across the UK with regard to the use of PMRT in intermediate risk breast cancer and patients with DCIS. Further work is required to determine which patients in these groups are likely to benefit from the use of PMRT.

Accelerated partial breast irradiation (APBI): Are breath-hold and volumetric radiation therapy techniques useful?

Background: In a selective group of patients accelerated partial breast irradiation (APBI) might be applied after conservative breast surgery to reduce the amount of irradiated healthy tissue. The role of volumetric modulated arc therapy (VMAT) and voluntary moderately deep inspiration breath-hold (vmDIBH) techniques in further reducing irradiated healthy – especially heart – tissue is investigated.

Material and methods: For 37 partial breast planning target volumes (PTVs), three-dimensional conformal radiotherapy (3D-CRT) (3 – 5 coplanar or non-coplanar 6 and/or 10 MV beams) and VMAT (two partial 6 MV arcs) plans were made on CTs acquired in free-breathing (FB) and/or in vmDIBH. Dose-volume parameters for the PTV, heart, lungs, and breasts were compared. 

Results: Better dose conformity was achieved with VMAT compared to 3D-CRT (conformity index 1.24 0.09 vs. 1.49 0.20). Non-PTV ipsilateral breast receiving 50% of the prescribed dose was on average reduced by 28% in VMAT plans compared to 3D-CRT plans. Mean heart dose (MHD) reduced from 2.0 (0.1 – 5.1) Gy in 3D-CRT(FB) to 0.6 (0.1 – 1.6) Gy in VMAT(vmDIBH). VMAT is benefi cial for MHD reduction if MHD with 3D-CRT exceeds 0.5Gy. Cardiac dose reduction as a result of VMAT increases with increasing initial MHD, and adding vmDIBH reduces the cardiac dose further. Mean dose to the ipsilateral lung decreased from 3.7 (0.7 – 8.7) to 1.8 (0.5 – 4.0) Gy with VMAT(vmDIBH) compared to 3D-CRT(FB). VMAT resulted in a slight increase in the contralateral breast dose (DMean ) always remaining 1.9 Gy). 

Conclusions: For APBI patients, VMAT improves PTV dose conformity and delivers lower doses to the ipsilateral breast and lung compared to 3D-CRT. This goes at the cost of a slight but acceptable increase of the contralateral breast dose. VMAT reduces cardiac dose if MHD exceeds 0.5 Gy for 3D-CRT. Adding vmDIBH results in a further reduction of heart and ipsilateral lung dose. 

Ecology of testate amoebae in an Amazonian peatland and development of a transfer function for palaeohydrological reconstruction

Tropical peatlands represent globally important carbon sinks with a unique biodiversity and are currently threatened by climate change and human activities. It is now imperative that proxy methods are developed to understand the ecohydrological dynamics of these systems and for testing peatland development models. Testate amoebae have been used as environmental indicators in ecological and palaeoecological studies of peatlands, primarily in ombrotrophic Sphagnum-dominated peatlands in the mid- and high-latitudes. We present the first ecological analysis of testate amoebae in a tropical peatland, a nutrient-poor domed bog in western (Peruvian) Amazonia. Litter samples were collected from different hydrological microforms (hummock to pool) along a transect from the edge to the interior of the peatland. We recorded 47 taxa from 21 genera. The most common taxa are Cryptodifflugia oviformis, Euglypha rotunda type, Phryganella acropodia, Pseudodifflugia fulva type and Trinema lineare. One species found only in the southern hemisphere, Argynnia spicata, is present. Arcella spp., Centropyxis aculeata and Lesqueresia spiralis are indicators of pools containing standing water. Canonical correspondence analysis and non-metric multidimensional scaling illustrate that water table depth is a significant control on the distribution of testate amoebae, similar to the results from mid- and high-latitude peatlands. A transfer function model for water table based on weighted averaging partial least-squares (WAPLS) regression is presented and performs well under cross-validation (r 2apparent=0.76,RMSE=4.29;r2jack=0.68,RMSEP=5.18. The transfer function was applied to a 1-m peat core, and sample-specific reconstruction errors were generated using bootstrapping. The reconstruction generally suggests near-surface water tables over the last 3,000 years, with a shift to drier conditions at c. cal. 1218-1273 AD

Molecular pedigree reconstruction and estimation of evolutionary parameters in a wild Atlantic salmon river system with incomplete sampling: a power analysis

Background: Pedigree reconstruction using genetic analysis provides a useful means to estimate fundamental population biology parameters relating to population demography, trait heritability and individual fitness when combined with other sources of data. However, there remain limitations to pedigree reconstruction in wild populations, particularly in systems where parent-offspring relationships cannot be directly observed, there is incomplete sampling of individuals, or molecular parentage inference relies on low quality DNA from archived material. While much can still be inferred from incomplete or sparse pedigrees, it is crucial to evaluate the quality and power of available genetic information a priori to testing specific biological hypotheses. Here, we used microsatellite markers to reconstruct a multi-generation pedigree of wild Atlantic salmon (Salmo salar L.) using archived scale samples collected with a total trapping system within a river over a 10 year period. Using a simulation-based approach, we determined the optimal microsatellite marker number for accurate parentage assignment, and evaluated the power of the resulting partial pedigree to investigate important evolutionary and quantitative genetic characteristics of salmon in the system.

Results: We show that at least 20 microsatellites (ave. 12 alleles/locus) are required to maximise parentage assignment and to improve the power to estimate reproductive success and heritability in this study system. We also show that 1.5 fold differences can be detected between groups simulated to have differing reproductive success, and that it is possible to detect moderate heritability values for continuous traits (h(2) similar to 0.40) with more than 80% power when using 28 moderately to highly polymorphic markers.

Conclusion: The methodologies and work flow described provide a robust approach for evaluating archived samples for pedigree-based research, even where only a proportion of the total population is sampled. The results demonstrate the feasibility of pedigree-based studies to address challenging ecological and evolutionary questions in free-living populations, where genealogies can be traced only using molecular tools, and that significant increases in pedigree assignment power can be achieved by using higher numbers of markers.

Molecular reconstruction of a genetic code alteration

The genetic code establishes the rules that govern gene translation into proteins. It was established more than 3.5 billion years ago and it is one of the most conserved features of life. Despite this, several alterations to the standard genetic code have been discovered in both prokaryotes and eukaryotes, namely in the fungal CTG clade where a unique seryl transfer RNA (tRNACAG Ser) decodes leucine CUG codons as serine. This tRNACAG Ser appeared 272±25 million years ago through insertion of an adenosine in the middle position of the anticodon of a tRNACGA Ser gene, which changed its anticodon from 5´-CGA-3´ to 5´-CAG-3´. This most dramatic genetic event restructured the proteome of the CTG clade species, but it is not yet clear how and why such deleterious genetic event was selected and became fixed in those fungal genomes. In this study we have attempted to shed new light on the evolution of this fungal genetic code alteration by reconstructing its evolutionary pathway in vivo in the yeast Saccharomyces cerevisiae. For this, we have expressed wild type and mutant versions of the C. albicans tRNACGA Ser gene into S. cerevisiae and evaluated the impact of the mutant tRNACGA Ser on fitness, tRNA stability, translation efficiency and aminoacylation kinetics. Our data demonstrate that these mutants are expressed and misincorporate Ser at CUGs, but their expression is repressed through an unknown molecular mechanism. We further demonstrate, using in vivo forced evolution methodologies, that the tRNACAG Ser can be easily inactivated through natural mutations that prevent its recognition by the seryl-tRNA synthetase. The overall data show that repression of expression of the mistranslating tRNACAG Ser played a critical role on the evolution of CUG reassignment from Leu to Ser. In order to better understand the evolution of natural genetic code alterations, we have also engineered partial reassignment of various codons in yeast. The data confirmed that genetic code ambiguity affects fitness, induces protein aggregation, interferes with the cell cycle and results in nuclear and morphologic alterations, genome instability and gene expression deregulation. Interestingly, it also generates phenotypic variability and phenotypes that confer growth advantages in certain environmental conditions. This study provides strong evidence for direct and critical roles of the environment on the evolution of genetic code alterations.

Reconstruction tridimensionnelle pour projection sur surfaces arbitraires.

Ce mémoire s'inscrit dans le domaine de la vision par ordinateur. Elle s'intéresse à la calibration de systèmes de caméras stéréoscopiques, à la mise en correspondance caméra-projecteur, à la reconstruction 3D, à l'alignement photométrique de projecteurs, au maillage de nuages de points, ainsi qu'au paramétrage de surfaces. Réalisé dans le cadre du projet LightTwist du laboratoire Vision3D, elle vise à permettre la projection sur grandes surfaces arbitraires à l'aide de plusieurs projecteurs. Ce genre de projection est souvent utilisé en arts technologiques, en théâtre et en projection architecturale. Dans ce mémoire, on procède au calibrage des caméras, suivi d'une reconstruction 3D par morceaux basée sur une méthode active de mise en correspondance, la lumière non structurée. Après un alignement et un maillage automatisés, on dispose d'un modèle 3D complet de la surface de projection. Ce mémoire introduit ensuite une nouvelle approche pour le paramétrage de modèles 3D basée sur le calcul efficace de distances géodésiques sur des maillages. L'usager n'a qu'à délimiter manuellement le contour de la zone de projection sur le modèle. Le paramétrage final est calculé en utilisant les distances obtenues pour chaque point du modèle. Jusqu'à maintenant, les méthodes existante ne permettaient pas de paramétrer des modèles ayant plus d'un million de points.