Inverse Modeling of GPS Multipath for Snow Depth Estimation-Part I: Formulation and Simulations

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Inverse Modeling of GPS Multipath for Snow Depth Estimation-Part II: Application and Validation

GPS multipath reflectometry (GPS-MR) is a technique that uses geodetic quality GPS receivers to estimate snow depth. The accuracy and precision of GPS-MR retrievals are evaluated at three different sites: grasslands, alpine, and forested. The assessment yields a correlation of 0.98 and an rms error of 6-8 cm for observed snow depths of up to 2.5 m. GPS-MR underestimates in situ snow depth by 10%-15% at these three sites, although the validation methods do not measure the same footprint as GPS-MR.

Comparison between multiplanar and rendering modes in the assessment of fetal atrioventricular valve areas by 3D/4D ultrasonography

Objective: To compare the agreement of multiplanar and rendering modes in the assessment fetal atrioventricular valves (mitral and tricuspid) areas by three-dimensional (3D) ultrasonography using the software spatio-temporal image correlation (STIC). Methods: We conducted a prospective cross-sectional study with normal pregnant women, with single fetuses, between 18-33 weeks. To measure the areas, we used the plan of four-chamber view. In the case of multiplanar, the plane was rotated on the axis "Z" form the heart to position at 9h. For rendering, the green line (region of interest - ROI) was placed from the atria of the heart perpendicular to the crux. The agreement was assessed by a Bland-Altman (limits of agreement) using the relative difference between the measures: ((rendering mode) - (multiplanar mode)) / (average). Results: 328 fetuses were evaluated. We have not identified the occurrence of systematic error between methods: the average relative difference was 1.62% (-2.07% to 5.32%, confidence interval 95%) in the mitral and 1.77% (- 1.08% to 4.62%) in the tricuspid valve. The limits of agreement between methods were -65.26% to 68.51% for the mitral and -49.91% to 53.45% for the tricuspid. Conclusions: There was no systematic error between modes and thus the observed values for the area of fetal atrioventricular valves can be used for comparisons needs to be corrected. However, relatively large variations may be observed when repeating the measurement area by different modes.

An inverse method to estimate the moving heat source in machining process

The present work propounds an inverse method to estimate the heat sources in the transient two-dimensional heat conduction problem in a rectangular domain with convective bounders. The non homogeneous partial differential equation (PDE) is solved using the Integral Transform Method. The test function for the heat generation term is obtained by the chip geometry and thermomechanical cutting. Then the heat generation term is estimated by the conjugated gradient method (CGM) with adjoint problem for parameter estimation. The experimental trials were organized to perform six different conditions to provide heat sources of different intensities. This method was compared with others in the literature and advantages are discussed. (C) 2012 Elsevier Ltd. All rights reserved.

Post-breakup tectonics in southeast Brazil from thermochronological data and combined inverse-forward thermal history modeling

The continental margin of southeast Brazil is elevated. Onshore Tertiary basins and Late Cretaceous/Paleogene intrusions are good evidence for post breakup tectono-magmatic activity. To constrain the impact of post-rift reactivation on the geological history of the area, we carried out a new thermochronological study. Apatite fission track ages range from 60.7 +/- 1.9 Ma to 129.3 +/- 4.3 Ma, mean track lengths from 11.41 +/- 0.23 mu m to 14.31 +/- 0.24 mu m and a subset of the (U-Th)/He ages range from 45.1 +/- 1.5 to 122.4 +/- 2.5 Ma. Results of inverse thermal history modeling generally support the conclusions from an earlier study for a Late Cretaceous phase of cooling. Around the onshore Taubate Basin, for a limited number of samples, the first detectable period of cooling occurred during the Early Tertiary. The inferred thermal histories for many samples also imply subsequent reheating followed by Neogene cooling. Given the uncertainty of the inversion results, we did deterministic forward modeling to assess the range of possibilities of this Tertiary part of the thermal history. The evidence for reheating seems to be robust around the Taubate Basin, but elsewhere the data cannot discriminate between this and a less complex thermal history. However, forward modeling results and geological information support the conclusion that the whole area underwent cooling during the Neogene. The synchronicity of the cooling phases with Andean tectonics and those in NE Brazil leads us to assume a plate-wide compressional stress that reactivated inherited structures. The present-day topographic relief of the margin reflects a contribution from post-breakup reactivation and uplift.

Pimavanserin, a selective serotonin (5-HT)2A-inverse agonist, enhances the efficacy and safety of risperidone, 2 mg/day, but does not enhance efficacy of haloperidol, 2 mg/day: Comparison with reference dose risperidone, 6 mg/day

Most atypical antipsychotic drugs (APDs), e. g. risperidone (RIS), produce more extensive blockade of brain serotonin (5-HT)(2A) than dopamine (DA) D-2 receptors. This distinguishes them from typical APDs, e.g. haloperidol (HAL). Our objective was to test the hypothesis that augmentation of low doses of RIS or HAL (2 mg/day) with pimavanserin (PIM), a selective 5-HT2A inverse agonist, to enhance 5-HT2A receptor blockade, can achieve efficacy comparable to RIS, 6 mg/day, but with lesser side effects. In a multi-center, randomized, double-blind, 6 week trial, 423 patients with chronic schizophrenia experiencing a recent exacerbation of psychotic symptoms were randomized to RIS2mg + placebo (RIS2PBO), RIS2mg + PIM20mg (RIS2PIM), RIS6mg + PBO (RIS6PBO), HAL2mg + PBO (HAL2PBO), or HAL2mg + PIM20mg (HAL2PIM). Improvement in psychopathology was measured by the PANSS and CGI-S. The reduction in PANSS Total Score with RIS2PIM at endpoint was significantly greater than RIS2PBO: -23.0 vs. -16.3 (p = 0.007), and not significantly different from the RIS6PBO group: -23.2 points. The percentage of patients with >= 20% improvement at day 15 in the RIS2PIM group was 62.3%, significantly greater than the RIS6PBO (42.1%; p = 0.01) and the RIS2PBO groups (37.7%; p = 0.002). Weight gain and hyperprolactinemia were greater in the RIS6PBO group than the RIS2PIM group but there was no difference in extrapyramidal side effects (EPS). HAL2PBO and HAL2PIM were not significantly different from each other in efficacy but HAL2PIM had less EPS at end point. Both HAL groups and RIS6PBO showed equal improvement in psychopathology at endpoint, indicating HAL 2 mg/day is effective to treat an acute exacerbation in chronic schizophrenia patients. In conclusion, a sub-effective RIS dose combined with PIM to enhance 5-HT2A receptor blockade provided faster onset of action, and at endpoint, equal efficacy and better safety, compared to standard dose RIS. These results support the conclusion that 5-HT2A receptor blockade is a key component of the action of some atypical APDs and can reduce EPS due to a typical APD. (C) 2012 Elsevier B.V. All rights reserved.

A graded subring of an inverse limit of polynomial rings

We study the power series ring R= K[[x1,x2,x3,...]]on countably infinitely many variables, over a field K, and two particular K-subalgebras of it: the ring S, which is isomorphic to an inverse limit of the polynomial rings in finitely many variables over K, and the ring R', which is the largest graded subalgebra of R. Of particular interest are the homogeneous, finitely generated ideals in R', among them the generic ideals. The definition of S as an inverse limit yields a set of truncation homomorphisms from S to K[x1,...,xn] which restrict to R'. We have that the truncation of a generic I in R' is a generic ideal in K[x1,...,xn]. It is shown in Initial ideals of Truncated Homogeneous Ideals that the initial ideal of such an ideal converge to the initial ideal of the corresponding ideal in R'. This initial ideal need no longer be finitely generated, but it is always locally finitely generated: this is proved in Gröbner Bases in R'. We show in Reverse lexicographic initial ideals of generic ideals are finitely generated that the initial ideal of a generic ideal in R' is finitely generated. This contrast to the lexicographic term order. If I in R' is a homogeneous, locally finitely generated ideal, and if we write the Hilbert series of the truncated algebras K[x1,...,xn] module the truncation of I as qn(t)/(1-t)n, then we show in Generalized Hilbert Numerators that the qn's converge to a power series in t which we call the generalized Hilbert numerator of the algebra R'/I. In Gröbner bases for non-homogeneous ideals in R' we show that the calculations of Gröbner bases and initial ideals in R' can be done also for some non-homogeneous ideals, namely those which have an associated homogeneous ideal which is locally finitely generated. The fact that S is an inverse limit of polynomial rings, which are naturally endowed with the discrete topology, provides S with a topology which makes it into a complete Hausdorff topological ring. The ring R', with the subspace topology, is dense in R, and the latter ring is the Cauchy completion of the former. In Topological properties of R' we show that with respect to this topology, locally finitely generated ideals in R'are closed.

Topics in geometry, analysis and inverse problems

The thesis consists of three independent parts. Part I: Polynomial amoebas We study the amoeba of a polynomial, as de ned by Gelfand, Kapranov and Zelevinsky. A central role in the treatment is played by a certain convex function which is linear in each complement component of the amoeba, which we call the Ronkin function. This function is used in two di erent ways. First, we use it to construct a polyhedral complex, which we call a spine, approximating the amoeba. Second, the Monge-Ampere measure of the Ronkin function has interesting properties which we explore. This measure can be used to derive an upper bound on the area of an amoeba in two dimensions. We also obtain results on the number of complement components of an amoeba, and consider possible extensions of the theory to varieties of codimension higher than 1. Part II: Differential equations in the complex plane We consider polynomials in one complex variable arising as eigenfunctions of certain differential operators, and obtain results on the distribution of their zeros. We show that in the limit when the degree of the polynomial approaches innity, its zeros are distributed according to a certain probability measure. This measure has its support on the union of nitely many curve segments, and can be characterized by a simple condition on its Cauchy transform. Part III: Radon transforms and tomography This part is concerned with different weighted Radon transforms in two dimensions, in particular the problem of inverting such transforms. We obtain stability results of this inverse problem for rather general classes of weights, including weights of attenuation type with data acquisition limited to a 180 degrees range of angles. We also derive an inversion formula for the exponential Radon transform, with the same restriction on the angle.

Computing inverse optical flow

[EN] We propose four algorithms for computing the inverse optical flow between two images. We assume that the forward optical flow has already been obtained and we need to estimate the flow in the backward direction. The forward and backward flows can be related through a warping formula, which allows us to propose very efficient algorithms. These are presented in increasing order of complexity. The proposed methods provide high accuracy with low memory requirements and low running times.In general, the processing reduces to one or two image passes. Typically, when objects move in a sequence, some regions may appear or disappear. Finding the inverse flows in these situations is difficult and, in some cases, it is not possible to obtain a correct solution. Our algorithms deal with occlusions very easy and reliably. On the other hand, disocclusions have to be overcome as a post-processing step. We propose three approaches for filling disocclusions. In the experimental results, we use standard synthetic sequences to study the performance of the proposed methods, and show that they yield very accurate solutions. We also analyze the performance of the filling strategies. 

A generalization of the Moore-Penrose inverse related to matrix subspaces of C(nxm)

[EN]A natural generalization of the classical Moore-Penrose inverse is presented. The so-called S-Moore-Penrose inverse of a m x n complex matrix A, denoted by As, is defined for any linear subspace S of the matrix vector space Cnxm. The S-Moore-Penrose inverse As is characterized using either the singular value decomposition or (for the nonsingular square case) the orthogonal complements with respect to the Frobenius inner product. These results are applied to the preconditioning of linear systems based on Frobenius norm minimization and to the linearly constrained linear least squares problem.

A generalization of the optimal diagonal approximate inverse preconditioner

[EN ]The classical optimal (in the Frobenius sense) diagonal preconditioner for large sparse linear systems Ax = b is generalized and improved. The new proposed approximate inverse preconditioner N is based on the minimization of the Frobenius norm of the residual matrix AM − I, where M runs over a certain linear subspace of n × n real matrices, defined by a prescribed sparsity pattern. The number of nonzero entries of the n×n preconditioning matrix N is less than or equal to 2n, and n of them are selected as the optimal positions in each of the n columns of matrix N. All theoretical results are justified in detail…

Bayesian Analysis of Linear Inverse Problems with Applications in Economics and Finance

In my PhD thesis I propose a Bayesian nonparametric estimation method for structural econometric models where the functional parameter of interest describes the economic agent's behavior. The structural parameter is characterized as the solution of a functional equation, or by using more technical words, as the solution of an inverse problem that can be either ill-posed or well-posed. From a Bayesian point of view, the parameter of interest is a random function and the solution to the inference problem is the posterior distribution of this parameter. A regular version of the posterior distribution in functional spaces is characterized. However, the infinite dimension of the considered spaces causes a problem of non continuity of the solution and then a problem of inconsistency, from a frequentist point of view, of the posterior distribution (i.e. problem of ill-posedness). The contribution of this essay is to propose new methods to deal with this problem of ill-posedness. The first one consists in adopting a Tikhonov regularization scheme in the construction of the posterior distribution so that I end up with a new object that I call regularized posterior distribution and that I guess it is solution of the inverse problem. The second approach consists in specifying a prior distribution on the parameter of interest of the g-prior type. Then, I detect a class of models for which the prior distribution is able to correct for the ill-posedness also in infinite dimensional problems. I study asymptotic properties of these proposed solutions and I prove that, under some regularity condition satisfied by the true value of the parameter of interest, they are consistent in a "frequentist" sense. Once I have set the general theory, I apply my bayesian nonparametric methodology to different estimation problems. First, I apply this estimator to deconvolution and to hazard rate, density and regression estimation. Then, I consider the estimation of an Instrumental Regression that is useful in micro-econometrics when we have to deal with problems of endogeneity. Finally, I develop an application in finance: I get the bayesian estimator for the equilibrium asset pricing functional by using the Euler equation defined in the Lucas'(1978) tree-type models.