A contour matching approach for accurate NOAA-AVHRR image navigation

Although different methods for NOAA AVHRR image navigation have already been established, the multitemporal and multi-satellite character of most studies requires automatic and accurate methods for navigation of satellite images. In the proposed method, a simple Kepplerian orbital model for the NOAA satellites is considered as reference model, and mean orbital elements are given as input to the model from ephemeris data. In order to correct the errors caused by these simplifications, errors resulting from inaccuracies in the positioning of the satellite and failures in the satellite internal clock, an automatic global contour matching approach has been adopted. First, the sensed image is preprocessed to obtain a gradient energy map of the reliable areas (sea-land contours) using a cloud detection algorithm and a morphological gradient operator. An initial estimation of the reliable contour positions is automatically obtained. The final positions of the contours are obtained by means of an iterative local minimization procedure that allows a contour to converge on an area of high image energy (edge). Global transformation parameters are estimated based on the initial and final positions of all reliable contour points. Finally, the performance of this approach is assessed using NOAA 14 AVHRR images from different geographic areas.Postprint (published version

Corruption as a Barrier to Entry: Theory and Evidence

Conventional wisdom depicts corruption as a tax on incumbent firms. This paper challenges this view in two ways. First, by arguing that corruption matters not so much because of the value of the bribe ("tax"), but because of another less studied feature of corruption, namely bribe unavoidability. Second, we argue that the social costs of corruption arise not because corruption hurts incumbent firms, but mostly because it acts as a powerful barrier to the entry of new firms. Corruption sands and greases in tandem: it helps incumbent firms (on balance) and it hurts potential entrants. We put forward a model in which a bureaucrat chooses entry barriers to optimize bribe revenues. When the capacity to collect bribes is high, it is optimal to allow high levels of oligopoly power to incumbents. Conversely, the more avoidable are the bribes, the more firms are allowed into the market. These ideas are tested using a unique, representative sample of Brazilian manufacturing firms. Consistently with our theoretical model, we show that corruption (a) is ranked as the most important barrier to entry (above finance, taxes and regulation) and (b) while bribes' unavoidability is positively related to firm performance, the size of the bribe is not.corruption, barriers to entry, firm performance

Normal fault earthquakes or graviquakes

Earthquakes are dissipation of energy throughout elastic waves. Canonically is the elastic energy accumulated during the interseismic period. However, in crustal extensional settings, gravity is the main energy source for hangingwall fault collapsing. Gravitational potential is about 100 times larger than the observed magnitude, far more than enough to explain the earthquake. Therefore, normal faults have a different mechanism of energy accumulation and dissipation (graviquakes) with respect to other tectonic settings (strike-slip and contractional), where elastic energy allows motion even against gravity. The bigger the involved volume, the larger is their magnitude. The steeper the normal fault, the larger is the vertical displacement and the larger is the seismic energy released. Normal faults activate preferentially at about 60° but they can be shallower in low friction rocks. In low static friction rocks, the fault may partly creep dissipating gravitational energy without releasing great amount of seismic energy. The maximum volume involved by graviquakes is smaller than the other tectonic settings, being the activated fault at most about three times the hypocentre depth, explaining their higher b-value and the lower magnitude of the largest recorded events. Having different phenomenology, graviquakes show peculiar precursor

Accurate and automatic NOAA-AVHRR image navigation using a global contour matching approach

The problem of precise and automatic AVHRR image navigation is tractable in theory, but has proved to be somewhat difficult in practice. The authors' work has been motivated by the need for a fully automatic and operational navigation system capable of geo-referencing NOAA-AVHRR images with high accuracy and without operator supervision. The proposed method is based on the simultaneous use of an orbital model and a contour matching approach. This last process, relying on an affine transformation model, is used to correct the errors caused by inaccuracies in orbit modeling, nonzero value for the spacecraft's roll, pitch and yaw, errors due to inaccuracies in the satellite positioning and failures in the satellite internal clock. The automatic global contour matching process is summarized as follows: i) Estimation of the gradient energy map (edges) in the sensed image and detection of the cloudless (reliable) areas in this map. ii) Initialization of the affine model parameters by minimizing the Euclidean distance between the reference and sensed images objects. iii) Simultaneous optimization of all reference image contours on the sensed image by energy minimization in the domain of the global transformation parameters. The process is iterated in a hierarchical way, reducing the parameter searching space at each iteration. The proposed image navigation algorithm has proved to be capable of geo-referencing a satellite image within 1 pixel.Peer ReviewedPostprint (published version

About multi-resolution techniques for large eddy simulation of reactive multi-phase flows

A numerical technique for mesh refinement in the HeaRT (Heat Release and Transfer) numerical code is presented. In the CFD framework, Large Eddy Simulation (LES) approach is gaining in importance as a tool for simulating turbulent combustion pro- cesses, also if this approach has an high computational cost due to the complexity of the turbulent modeling and the high number of grid points necessary to obtain a good numerical solution. In particular, when a numerical simulation of a big domain is performed with a structured grid, the number of grid points can increase so much that the simulation becomes impossible: this problem can be overcomed with a mesh refinement technique. Mesh refinement technique developed for HeaRT numerical code (a staggered finite difference code) is based on an high order reconstruction of the variables at the grid interfaces by means of a least square quasi-eno interpolation: numerical code is written in modern Fortran (2003 standard of newer) and is parallelized using domain decomposition and message passing interface (MPI) standard

Generic features of the fluctuation dissipation relation in coarsening systems

The integrated response function in phase-ordering systems with scalar, vector, conserved and non conserved order parameter is studied at various space dimensionalities. Assuming scaling of the aging contribution $\chi_{ag} (t,t_w)= t_w ^{-a_\chi} \hat \chi (t/t_w)$ we obtain, by numerical simulations and analytical arguments, the phenomenological formula describing the dimensionality dependence of $a_\chi$ in all cases considered. The primary result is that $a_\chi$ vanishes continuously as $d$ approaches the lower critical dimensionality $d_L$. This implies that i) the existence of a non trivial fluctuation dissipation relation and ii) the failure of the connection between statics and dynamics are generic features of phase ordering at $d_L$.Comment: 6 pages, 5 figure

Scaling and universality in the aging kinetics of the two-dimensional clock model

We study numerically the aging dynamics of the two-dimensional p-state clock model after a quench from an infinite temperature to the ferromagnetic phase or to the Kosterlitz-Thouless phase. The system exhibits the general scaling behavior characteristic of non-disordered coarsening systems. For quenches to the ferromagnetic phase, the value of the dynamical exponents, suggests that the model belongs to the Ising-type universality class. Specifically, for the integrated response function $\chi (t,s)\simeq s^{-a_\chi}f(t/s)$, we find $a_\chi$ consistent with the value $a_\chi =0.28$ found in the two-dimensional Ising model.Comment: 16 pages, 14 figures (please contact the authors for figures