Digitoitujen väri-infrailmakuvien numeerinen monikuvatulkinta metsäninventoinnissa

The use of remote sensing imagery as auxiliary data in forest inventory is based on the correlation between features extracted from the images and the ground truth. The bidirectional reflectance and radial displacement cause variation in image features located in different segments of the image but forest characteristics remaining the same. The variation has so far been diminished by different radiometric corrections. In this study the use of sun azimuth based converted image co-ordinates was examined to supplement auxiliary data extracted from digitised aerial photographs. The method was considered as an alternative for radiometric corrections. Additionally, the usefulness of multi-image interpretation of digitised aerial photographs in regression estimation of forest characteristics was studied. The state owned study area located in Leivonmäki, Central Finland and the study material consisted of five digitised and ortho-rectified colour-infrared (CIR) aerial photographs and field measurements of 388 plots, out of which 194 were relascope (Bitterlich) plots and 194 were concentric circular plots. Both the image data and the field measurements were from the year 1999. When examining the effect of the location of the image point on pixel values and texture features of Finnish forest plots in digitised CIR photographs the clearest differences were found between front-and back-lighted image halves. Inside the image half the differences between different blocks were clearly bigger on the front-lighted half than on the back-lighted half. The strength of the phenomenon varied by forest category. The differences between pixel values extracted from different image blocks were greatest in developed and mature stands and smallest in young stands. The differences between texture features were greatest in developing stands and smallest in young and mature stands. The logarithm of timber volume per hectare and the angular transformation of the proportion of broadleaved trees of the total volume were used as dependent variables in regression models. Five different converted image co-ordinates based trend surfaces were used in models in order to diminish the effect of the bidirectional reflectance. The reference model of total volume, in which the location of the image point had been ignored, resulted in RMSE of 1,268 calculated from test material. The best of the trend surfaces was the complete third order surface, which resulted in RMSE of 1,107. The reference model of the proportion of broadleaved trees resulted in RMSE of 0,4292 and the second order trend surface was the best, resulting in RMSE of 0,4270. The trend surface method is applicable, but it has to be applied by forest category and by variable. The usefulness of multi-image interpretation of digitised aerial photographs was studied by building comparable regression models using either the front-lighted image features, back-lighted image features or both. The two-image model turned out to be slightly better than the one-image models in total volume estimation. The best one-image model resulted in RMSE of 1,098 and the two-image model resulted in RMSE of 1,090. The homologous features did not improve the models of the proportion of broadleaved trees. The overall result gives motivation for further research of multi-image interpretation. The focus may be improving regression estimation and feature selection or examination of stratification used in two-phase sampling inventory techniques. Keywords: forest inventory, digitised aerial photograph, bidirectional reflectance, converted image co­ordinates, regression estimation, multi-image interpretation, pixel value, texture, trend surface

Individual-based modelling of the efficacy of fumigation tactics to control lesser grain borer (Rhyzopertha dominica) in stored grain

Increasing resistance to phosphine (PH 3) in insect pests, including lesser grain borer (Rhyzopertha dominica) has become a critical issue, and development of effective and sustainable strategies to manage resistance is crucial. In practice, the same grain store may be fumigated multiple times, but usually for the same exposure period and concentration. Simulating a single fumigation allows us to look more closely at the effects of this standard treatment.We used an individual-based, two-locus model to investigate three key questions about the use of phosphine fumigant in relation to the development of PH 3 resistance. First, which is more effective for insect control; long exposure time with a low concentration or short exposure period with a high concentration? Our results showed that extending exposure duration is a much more efficient control tactic than increasing the phosphine concentration. Second, how long should the fumigation period be extended to deal with higher frequencies of resistant insects in the grain? Our results indicated that if the original frequency of resistant insects is increased n times, then the fumigation needs to be extended, at most, n days to achieve the same level of insect control. The third question is how does the presence of varying numbers of insects inside grain storages impact the effectiveness of phosphine fumigation? We found that, for a given fumigation, as the initial population number was increased, the final survival of resistant insects increased proportionally. To control initial populations of insects that were n times larger, it was necessary to increase the fumigation time by about n days. Our results indicate that, in a 2-gene mediated resistance where dilution of resistance gene frequencies through immigration of susceptibles has greater effect, extending fumigation times to reduce survival of homozygous resistant insects will have a significant impact on delaying the development of resistance. © 2012 Elsevier Ltd.

Dosage consistency is the key factor in avoiding evolution of resistance to phosphine and population increase in stored-grain pests

BACKGROUND Control of pests in stored grain and the evolution of resistance to pesticides are serious problems worldwide. A stochastic individual-based two-locus model was used to investigate the impact of two important issues, the consistency of pesticide dosage through the storage facility and the immigration rate of the adult pest, on overall population control and avoidance of evolution of resistance to the fumigant phosphine in an important pest of stored grain, the lesser grain borer. RESULTS A very consistent dosage maintained good control for all immigration rates, while an inconsistent dosage failed to maintain control in all cases. At intermediate dosage consistency, immigration rate became a critical factor in whether control was maintained or resistance emerged. CONCLUSION Achieving a consistent fumigant dosage is a key factor in avoiding evolution of resistance to phosphine and maintaining control of populations of stored-grain pests; when the dosage achieved is very inconsistent, there is likely to be a problem regardless of immigration rate. © 2012 Society of Chemical Industry

A Probabilistic Approach to the Primary Visual Cortex

What can the statistical structure of natural images teach us about the human brain? Even though the visual cortex is one of the most studied parts of the brain, surprisingly little is known about how exactly images are processed to leave us with a coherent percept of the world around us, so we can recognize a friend or drive on a crowded street without any effort. By constructing probabilistic models of natural images, the goal of this thesis is to understand the structure of the stimulus that is the raison d etre for the visual system. Following the hypothesis that the optimal processing has to be matched to the structure of that stimulus, we attempt to derive computational principles, features that the visual system should compute, and properties that cells in the visual system should have. Starting from machine learning techniques such as principal component analysis and independent component analysis we construct a variety of sta- tistical models to discover structure in natural images that can be linked to receptive field properties of neurons in primary visual cortex such as simple and complex cells. We show that by representing images with phase invariant, complex cell-like units, a better statistical description of the vi- sual environment is obtained than with linear simple cell units, and that complex cell pooling can be learned by estimating both layers of a two-layer model of natural images. We investigate how a simplified model of the processing in the retina, where adaptation and contrast normalization take place, is connected to the nat- ural stimulus statistics. Analyzing the effect that retinal gain control has on later cortical processing, we propose a novel method to perform gain control in a data-driven way. Finally we show how models like those pre- sented here can be extended to capture whole visual scenes rather than just small image patches. By using a Markov random field approach we can model images of arbitrary size, while still being able to estimate the model parameters from the data.

Determination of optimal trajectory under design constraints for a satellite launch vehicle

An optimal pitch steering programme of a solid-fuel satellite launch vehicle to maximize either (1) the injection velocity at a given altitude, or (2) the size of circular orbit, for a given payload is presented. The two-dimensional model includes the rotation of atmosphere with the Earth, the vehicle's lift and drag, variation of thrust with time and altitude, inverse-square gravitational field, and the specified initial vertical take-off. The inequality constraints on the aerodynamic load, control force, and turning rates are also imposed. Using the properties of the central force motion the terminal constraint conditions at coast apogee are transferred to the penultimate stage burnout. Such a transformation converts a time-free problem into a time-fixed one, reduces the number of terminal constraints, improves accuracy, besides demanding less computer memory and time. The adjoint equations are developed in a compact matrix form. The problem is solved on an IBM 360/44 computer using a steepest ascent algorithm. An illustrative analysis of a typical launch vehicle establishes the speed of convergence, and accuracy and applicability of the algorithm.

The Macroeconomic Effects of Fiscal Policy in a Small Open Economy under a Flexible Exchange Rate Regime

This licentiate's thesis analyzes the macroeconomic effects of fiscal policy in a small open economy under a flexible exchange rate regime, assuming that the government spends exclusively on domestically produced goods. The motivation for this research comes from the observation that the literature on the new open economy macroeconomics (NOEM) has focused almost exclusively on two-country global models and the analyses of the effects of fiscal policy on small economies are almost completely ignored. This thesis aims at filling in the gap in the NOEM literature and illustrates how the macroeconomic effects of fiscal policy in a small open economy depend on the specification of preferences. The research method is to present two theoretical model that are extensions to the model contained in the Appendix to Obstfeld and Rogoff (1995). The first model analyzes the macroeconomic effects of fiscal policy, making use of a model that exploits the idea of modelling private and government consumption as substitutes in private utility. The model offers intuitive predictions on how the effects of fiscal policy depend on the marginal rate of substitution between private and government consumption. The findings illustrate that the higher the substitutability between private and government consumption, (i) the bigger is the crowding out effect on private consumption (ii) and the smaller is the positive effect on output. The welfare analysis shows that the less fiscal policy decreases welfare the higher is the marginal rate of substitution between private and government consumption. The second model of this thesis studies how the macroeconomic effects of fiscal policy depend on the elasticity of substitution between traded and nontraded goods. This model reveals that this elasticity a key variable to explain the exchange rate, current account and output response to a permanent rise in government spending. Finally, the model demonstrates that temporary changes in government spending are an effective stabilization tool when used wisely and timely in response to undesired fluctuations in output. Undesired fluctuations in output can be perfectly offset by an opposite change in government spending without causing any side-effects.

Analysis of human chorionic gonadotropin-monoclonal antibody interaction in BIAcore

Kinetic studies of macromolecular ligand-ligate interaction have generated ample interest since the advent of plasmon resonance based instruments like BIAcore. Most of the studies reported in literature assume a simple 1 : 1 Langmuir binding and complete reversibility of the system. However we observed that in a high affinity antigen-antibody system [human chorionic gonadotropin-monoclonal antibody (hCG-mAb)] dissociation is insignificant and the sensogram data cannot be used to measure the equilibrium and kinetic parameters. At low concentrations of mAb the complete sensogram could be fitted to a single exponential. Interestingly we found that at higher mAb concentrations, the binding data did not conform to a simple bimolecular model. Instead, the data fitted a two-step model, which may be because of surface heterogeneity of affinity sites. In this paper, we report on the global fit of the sensograms. We have developed a method by which a single two-minute sensogram can be used in high affinity systems to measure the association rate constant of the reaction and the functional capacity of the ligand (hCG) immobilized on the chip. We provide a rational explanation for the discrepancies generally observed in most of the BIAcore sensograms

Foreign aid, domestic investment and welfare

This paper examines the welfare implications of temporary foreign aid in a simple two-period, two-country model of trade. Domestic investment is endogenous, providing an important link between aid in period one and the terms of trade in periods one and two. Transfer-induced changes in the terms of trade redistribute present and future income between the donor and the recipient. In the presence of barriers to international borrowing and lending, such redistribution gives rise to the possibility of temporary aid being both potentially and strictly Pareto improving.

Logic of aid in an intertemporal setting

The paper studies the welfare implications of temporary foreign aid in the context of a simple two-country model of trade. In addition to its usual effects, a transfer of income in one period is assumed to influence the preferences of the recipient country in the following period. The implied changes in the terms of trade over the two periods are consistent with a number of possible outcomes with respect to the intertemporal welfare of the donor, the recipient, and the world as a whole. Particular attention is devoted to the conditions for strict Pareto improvement and the circumstances under which temporary aid transactions are likely to occur.

The transfer problem and the intertemporal terms of trade

In this paper the effects of a transfer on the intertemporal terms of trade are examined in the context of a simple two-country, two-period model. When intertemporal trade occurs because the two economies have different rates of time preference, a transfer improves the terms of trade of the paying country. Alternatively, when trade occurs owing to international differences in the endowments of goods over the two periods, the effect of a transfer depends on (a) the relationship between the interest rate and the rates of time preference of the two countries and (b) the relationship between their elasticities of intertemporal consumption substitution.

Phases in the (1-x)Na0.5Bi0.5TiO3-(x)CaTiO3 system

The structures of (1 - x) Na0.5Bi0.5TiO3-(x) CaTiO3 at room temperature have been investigated using neutron powder diffraction and dielectric studies. The system exhibits an orthorhombic (Pbnm) structure for x >= 0.15 and rhombohedral (R3c) for x <= 0.05. For x = 0.10, though the neutron diffraction pattern shows features of the orthorhombic (Pbnm) structure, Rietveld refinement using this structure shows a drastic reduction in the in-phase tilt angle (similar to 4 degrees) as compared to the corresponding value (similar to 8 degrees) for a neighbouring composition x = 0.15. The neutron diffraction pattern of x = 0.10 could be fitted equally well using a two-phase model (R3c + Pbnm) with orthorhombic as the minor phase (22%), without the need for a drastic decrease in the in-phase tilt angle. The dielectric studies of x = 0.10 revealed the presence of the polar R3c phase, thereby favouring the phase coexistence model, instead of a single-phase Pbnm structure, for this composition.

Quasiequilibrium optical nonlinearities from spin-polarized carriers in GaAs

Semiconductor Bloch equations, which microscopically describe the dynamics of a Coulomb interacting, spin-unpolarized electron-hole plasma, can be solved in two limits: the coherent and the quasiequilibrium regimes. These equations have been recently extended to include the spin degree of freedom and used to explain spin dynamics in the coherent regime. In the quasiequilibrium limit, one solves the Bethe-Salpeter equation in a two-band model to describe how optical absorption is affected by Coulomb interactions within a spin unpolarized plasma of arbitrary density. In this work, we modified the solution of the Bethe-Salpeter equation to include spin polarization and light holes in a three-band model, which allowed us to account for spin-polarized versions of many-body effects in absorption. The calculated absorption reproduced the spin-dependent, density-dependent, and spectral trends observed in bulk GaAs at room temperature, in a recent pump-probe experiment with circularly polarized light. Hence, our results may be useful in the microscopic modeling of density-dependent optical nonlinearities due to spin-polarized carriers in semiconductors.

Effects of climate and morphology on temperature conditions of lakes

To a large extent, lakes can be described with a one-dimensional approach, as their main features can be characterized by the vertical temperature profile of the water. The development of the profiles during the year follows the seasonal climate variations. Depending on conditions, lakes become stratified during the warm summer. After cooling, overturn occurs, water cools and an ice cover forms. Typically, water is inversely stratified under the ice, and another overturn occurs in spring after the ice has melted. Features of this circulation have been used in studies to distinguish between lakes in different areas, as basis for observation systems and even as climate indicators. Numerical models can be used to calculate temperature in the lake, on the basis of the meteorological input at the surface. The simple form is to solve the surface temperature. The depth of the lake affects heat transfer, together with other morphological features, the shape and size of the lake. Also the surrounding landscape affects the formation of the meteorological fields over the lake and the energy input. For small lakes the shading by the shores affects both over the lake and inside the water body bringing limitations for the one-dimensional approach. A two-layer model gives an approximation for the basic stratification in the lake. A turbulence model can simulate vertical temperature profile in a more detailed way. If the shape of the temperature profile is very abrupt, vertical transfer is hindered, having many important consequences for lake biology. One-dimensional modelling approach was successfully studied comparing a one-layer model, a two-layer model and a turbulence model. The turbulence model was applied to lakes with different sizes, shapes and locations. Lake models need data from the lakes for model adjustment. The use of the meteorological input data on different scales was analysed, ranging from momentary turbulent changes over the lake to the use of the synoptical data with three hour intervals. Data over about 100 past years were used on the mesoscale at the range of about 100 km and climate change scenarios for future changes. Increasing air temperature typically increases water temperature in epilimnion and decreases ice cover. Lake ice data were used for modelling different kinds of lakes. They were also analyzed statistically in global context. The results were also compared with results of a hydrological watershed model and data from very small lakes for seasonal development.

Hydrogen soil deposition and atmospheric variations in the boreal zone

This research has been prompted by an interest in the atmospheric processes of hydrogen. The sources and sinks of hydrogen are important to know, particularly if hydrogen becomes more common as a replacement for fossil fuel in combustion. Hydrogen deposition velocities (vd) were estimated by applying chamber measurements, a radon tracer method and a two-dimensional model. These three approaches were compared with each other to discover the factors affecting the soil uptake rate. A static-closed chamber technique was introduced to determine the hydrogen deposition velocity values in an urban park in Helsinki, and at a rural site at Loppi. A three-day chamber campaign to carry out soil uptake estimation was held at a remote site at Pallas in 2007 and 2008. The atmospheric mixing ratio of molecular hydrogen has also been measured by a continuous method in Helsinki in 2007 - 2008 and at Pallas from 2006 onwards. The mean vd values measured in the chamber experiments in Helsinki and Loppi were between 0.0 and 0.7 mm s-1. The ranges of the results with the radon tracer method and the two-dimensional model were 0.13 - 0.93 mm s-1 and 0.12 - 0.61 mm s-1, respectively, in Helsinki. The vd values in the three-day campaign at Pallas were 0.06 - 0.52 mm s-1 (chamber) and 0.18 - 0.52 mm s-1 (radon tracer method and two-dimensional model). At Kumpula, the radon tracer method and the chamber measurements produced higher vd values than the two-dimensional model. The results of all three methods were close to each other between November and April, except for the chamber results from January to March, while the soil was frozen. The hydrogen deposition velocity values of all three methods were compared with one-week cumulative rain sums. Precipitation increases the soil moisture, which decreases the soil uptake rate. The measurements made in snow seasons showed that a thick snow layer also hindered gas diffusion, lowering the vd values. The H2 vd values were compared to the snow depth. A decaying exponential fit was obtained as a result. During a prolonged drought in summer 2006, soil moisture values were lower than in other summer months between 2005 and 2008. Such conditions were prevailing in summer 2006 when high chamber vd values were measured. The mixing ratio of molecular hydrogen has a seasonal variation. The lowest atmospheric mixing ratios were found in the late autumn when high deposition velocity values were still being measured. The carbon monoxide (CO) mixing ratio was also measured. Hydrogen and carbon monoxide are highly correlated in an urban environment, due to the emissions originating from traffic. After correction for the soil deposition of H2, the slope was 0.49±0.07 ppb (H2) / ppb (CO). Using the corrected hydrogen-to-carbon-monoxide ratio, the total hydrogen load emitted by Helsinki traffic in 2007 was 261 t (H2) a-1. Hydrogen, methane and carbon monoxide are connected with each other through the atmospheric methane oxidation process, in which formaldehyde is produced as an important intermediate. The photochemical degradation of formaldehyde produces hydrogen and carbon monoxide as end products. Examination of back-trajectories revealed long-range transportation of carbon monoxide and methane. The trajectories can be grouped by applying cluster and source analysis methods. Thus natural and anthropogenic emission sources can be separated by analyzing trajectory clusters.

Studies in bubble formation — III

Bubble formation from single horizontal orifices submerged in Newtonian liquids has been investigated for such chamber volumes that both the pressure inside the chamber and flow rate into the bubble are time dependent. The data collected show that under these conditions the bubble volume decreases exponentially with increase in orifice submergence. The equations for the generalized two stage model of bubble formation, taking the variation of gas flow rate with time into account, have been derived. These equations reduce to the cases of constant gas flow rate and constant pressure when adequate constraints are imposed. The results obtained under intermediate conditions have been quantitatively explained on the basis of these equations.