Estimating areal means and variances of forest attributes using the k-Nearest Neighbors technique and satellite imagery

The k-Nearest Neighbor (k-NN) technique has become extremely popular for a variety of forest inventory mapping and estimation applications. Much of this popularity may be attributed to the non-parametric, multivariate features of the technique, its intuitiveness, and its ease of use. When used with satellite imagery and forest inventory plot data, the technique has been shown to produce useful estimates of many forest attributes including forest/non-forest, volume, and basal area. However, variance estimators for quantifying the uncertainty of means or sums of k-NN pixel-level predictions for areas of interest (AOI) consisting of multiple pixels have not been reported. The primary objectives of the study were to derive variance estimators for AOI estimates obtained from k-NN predictions and to compare precision estimates resulting from different approaches to k-NN prediction and different interpretations of those predictions. The approaches were illustrated by estimating proportion forest area, tree volume per unit area, tree basal area per unit area, and tree density per unit area for 10-km AOIs. Estimates obtained using k-NN approaches and traditional inventory approaches were compared and found to be similar. Further, variance estimates based on different interpretations of k-NN predictions were similar. The results facilitate small area estimation and simultaneous and consistent mapping and estimation of multiple forest attributes.     

Variation among and within mountain birch trees in foliage phenols,carbohydrates, and amino acids,and in growth ofEpirrita autumnata larvae

Leaf quality of the mountain birch (Betula pubescens ssp.tortuosa) for herbivores was studied at several hierarchical levels: among trees, among ramets within trees, among branches within ramets, and among short shoots within branches. The experimental units at each level were chosen randomly. The indices of leaf quality were the growth rate of the larvae of a geometrid,Epirrita autumnata, and certain biochemical traits of the leaves (total phenolics and individual phenolic compounds, total carbohydrates and individual sugars, free and protein-bound amino acids). We also discuss relationships between larval growth rate and biochemical foliage traits. Larval growth rates during two successive years correlated positively at the level of tree, the ramet, and the branch, indicating that the relationships in leaf quality remained constant between seasons both among and within trees. The distribution of variation at different hierarchical levels depended on the trait in question. In the case of larval growth rate, ramets and short shoots accounted for most of the explained variation. In the case of biochemical compounds, trees accounted for most of the variance in the content of total phenolics and individual low-molecular-weight phenolics. In the content of carbohydrates (total carbohydrates, starch, fructose, glucose, and sucrose) and amino acids, variation among branches was generally larger than variation among trees. Variation among ramets was low for most compounds. No single leaf trait played a paramount role in larval growth. Secondary compounds, represented by phenolic compounds, or primary metabolites, particularly sugars, may both be important in determining the suitability of birch leaves for larvae. If phenols are causally more important, genet-specific analyses of foliage chemistry are needed. If sugars are of primary importance, within-genet sampling and analysis of foliage chemistry are necessary.     

Using k-nn and discriminant analyses to classify rain forest types in a Landsat TM image over northern Costa Rica

Conservation and land use planning in humid tropical lowland forests urgently need accurate remote sensing techniques to distinguish among floristically different forest types. We investigated the degree to which floristically and structurally defined Costa Rican lowland rain forest types can be accurately discriminated by a non-parametric k nearest neighbors (k-nn) classifier or linear discriminant analysis. Pixel values of Landsat Thematic Mapper (TM) image and Shuttle Radar Topography Mission (SRTM) elevation model extracted from segments or from 5 × 5 pixel windows were employed in the classifications. 104 field plots were classified into three floristic and one structural type of forest (regrowth forest). Three floristically defined forest types were formed through clustering the old-growth forest plots (n = 52) by their species specific importance values. An error assessment of the image classification was conducted via cross-validation and error matrices, and overall percent accuracy and Kappa scores were used as measures of accuracy. Image classification of the four forest types did not adequately distinguish two old-growth forest classes, so they were merged into a single forest class. The resulting three forest classes were most accurately classified by the k-nn classifier using segmented image data (overall accuracy 91%). The second best method, with respect to accuracy, was the k-nn with 5 × 5 pixel windows data (89% accuracy), followed by the canonical discriminant analysis using the 5 × 5 pixel window data (86%) and the segment data (82%). We conclude the k-nn classifier can accurately distinguish floristically and structurally different rain forest types. The classification accuracies were higher for the k-nn classifier than for the canonical discriminant analysis, but the differences in Kappa scores were not statistically significant. The segmentation did not increase classification accuracy in this study.     

Results of Finnish national survey on EU legislation concerning computer work

The European Directive on computer work (VDU 90/270/EEC) is implemented in the Finnish Government Decree. The aim of the present study was to evaluate the effects of the legislation and its applications in practice. The quantitative method used an online questionnaire. The respondents were employers (N = 934), employees (N = 1872) and occupational health care (OHC) units (N = 289). The majority of all these three groups considered the provisions clear and easy to understand, comprehensive, and easy to comply with. The provisions had a great impact on preventing mental overloading, arranging sight tests at the employer's cost, and on the refunding of eyeglasses for computer work. Most employers felt they did not need any more specific provisions, whereas half of the employees and OHC professionals would have liked the provisions to be more detailed. More explicit and specific regulations were also needed for practice applications in particular, as well as for the refunding of the costs of special eyeglasses for computer work. A total of 59% of the employers reported that costs of glasses are refunded, whereas the figure for employees was clearly lower (37%). OHC professionals reported that the costs of glasses were refunded in 82% of their customer workplaces. The practical conclusion is that employees' sight examinations and compensation for eyeglasses should be promoted in computer work. Moreover, ergonomic applications are best carried out in co-operation with OHC personnel.     

An efficient deconvolution algorithm for estimating oxygen consumption during muscle activities

The reconstruction of an unknown input function from noisy measurements in a biological system is an ill-posed inverse problem. Any computational algorithm for its solution must use some kind of regularization technique to neutralize the disastrous effects of amplified noise components on the computed solution. In this paper, following a hierarchical Bayesian statistical inversion approach, we seek estimates for the input function and regularization parameter (hyperparameter) that maximize the posterior probability density function. We solve the maximization problem simultaneously for all unknowns, hyperparameter included, by a suitably chosen quasi-Newton method. The optimization approach is compared to the sampling-based Bayesian approach. We demonstrate the efficiency and robustness of the deconvolution algorithm by applying it to reconstructing the time courses of mitochondrial oxygen consumption during muscle state transitions (e.g., from resting state to contraction and recovery), from the simulated noisy output of oxygen concentration dynamics on the muscle surface. The model of oxygen transport and metabolism in skeletal muscle assumes an in vitro cylindrical structure of the muscle in which the oxygen from the surrounding oxygenated solution diffuses into the muscle and is then consumed by the muscle mitochondria. The algorithm can be applied to other deconvolution problems by suitably replacing the forward model of the system.     

Quadratic nonnegative matrix factorization

In Nonnegative Matrix Factorization (NMF), a nonnegative matrix is approximated by a product of lower-rank factorizing matrices. Most NMF methods assume that each factorizing matrix appears only once in the approximation, thus the approximation is linear in the factorizing matrices. We present a new class of approximative NMF methods, called Quadratic Nonnegative Matrix Factorization (QNMF), where some factorizing matrices occur twice in the approximation. We demonstrate QNMF solutions to four potential pattern recognition problems in graph partitioning, two-way clustering, estimating hidden Markov chains, and graph matching. We derive multiplicative algorithms that monotonically decrease the approximation error under a variety of measures. We also present extensions in which one of the factorizing matrices is constrained to be orthogonal or stochastic. Empirical studies show that for certain application scenarios, QNMF is more advantageous than other existing nonnegative matrix factorization methods.     

Easy Algorithm Animation on the Web

The animation environment Jeliot automates visualization of algorithms over the Internet. The user can visualize algorithms of his own, not merely the selection supplied by a service provider. Jeliot is based on self-animating data types: the user selects the data objects of the source code to be visualized, and Jeliot produces the animation based on operations performed on those types. We present the design of Jeliot and the principles of its implementation. In addition, we discuss the educational applications of Jeliot.     

On the existence and convergence of the solution of PML equations

In this article we study the mesh termination method in computational scattering theory known as the method of Perfectly Matched Layer (PML). This method is based on the idea of surrounding the scatterer and its immediate vicinity with a fictitious absorbing non-reflecting layer to damp the echoes coming from the mesh termination surface. The method can be formulated equivalently as a complex stretching of the exterior domain. The article is devoted to the existence and convergence questions of the solutions of the resulting equations. We show that with a special choice of the fictitious absorbing coefficient, the PML equations are solvable for all wave numbers, and as the PML layer is made thicker, the PML solution converge exponentially towards the actual scattering solution. The proofs are based on boundary integral methods and a new type of near-field version of the radiation condition, called here the double surface radiation condition. Partly supported by the Finnish Academy, project 37692.     

A hybrid approach to computing electrostatic forces in fluidized beds of charged particles

In particulate flow devices particles acquire electric charge through triboelectric charging, and resulting electrostatic forces can alter hydrodynamics. To capture this effect, the electrostatic force acting on individual particles in the device should be computed accurately. Electrostatic force is calculated using a hybrid approach consisting of: (1) long‐range contributions from an Eulerian electric field solved using the Poisson equation (2) short‐range contributions calculated using a truncated pairwise sum and (3) a correction to avoid double counting. Euler‐Lagrange simulation of flows incorporating this hybrid approach reveals that bed height oscillations in small fluidized beds of particles with monopolar charge decreases with increasing charge level, which is related to lateral segregation of particles. A ring‐like layer of particles, reported in experimental studies, forms at modestly high charge levels. Beds with equal amounts of positively and negatively charged particles are fluidized in a manner similar to uncharged particles. © 2016 American Institute of Chemical Engineers AIChE J, 62: 2282–2295, 2016