Polarimetric SAR characterization of man-made structures in urban areas using normalized circular-pol correlation coefficients

Polarimetric Synthetic Aperture Radar (SAR) backscatter from man-made structures in urban areas is quite different than backscatter from predominantly natural areas. Backscatter from natural areas is often reflection symmetric; i.e., characterized by near zero values for covariance matrix off-diagonal terms of the form 〈S_HVS_HH^?〉, 〈S_HVS_VV^?〉 and their conjugates. A new approach is proposed to detect scattering from non-reflection symmetric structures using circular-pol, RR-LL, correlation coefficients, |ρ|. This method creates a normalization term, |ρ₀|, and then forms a ratio, |ρ|/|ρ₀|. The normalization term, |ρ₀|, contains the same diagonal terms of the covariance matrix. The 〈S_HVS_HH^?〉 and 〈S_HVS_VV^?〉 off-diagonal terms and their conjugates are purposely set to zero. The ratio, |ρ|/|ρ₀|, is rewritten as a product of separable helicity (τ) and orientation angle (θ) dependencies. The mathematical form of the τ dependence is a resonant singularity, or pole, term. This pole significantly enhances returns from man-made, high helicity, non-reflection symmetric structures. These structures have values of τ near the resonance value at τ = ± 1. Natural scatterers possess very strong RR / LL symmetry (τ ≈ 0) and the pole response for them is correspondingly weak. The dependence of |ρ|/|ρ₀| on the orientation angle (θ) is known from previous studies to be useful for measuring urban building alignments (relative to the azimuth direction) and measuring surface topography. The ratio |ρ|/|ρ₀| reduces much of the un-needed image detail of backscatter variations from natural areas of different surface roughness. This image simplification further facilitates detection of localized man-made targets.     

Setting due dates in a stochastic single machine environment

A set of n   jobs with statistically independent random processing times has to be processed on a single machine without idling between jobs and without preemption. It is required to set due dates and promise them to customers. During the production stage, earliness and tardiness against the promised due dates will be penalized. The goal is to minimize the total expected penalties. We consider two due date setting procedures with optimum customer service level, and an O(nlogn)$\mathrm{O}(n\log n)$ time complexity. We show that one is asymptotically optimal but the other is not. Both heuristics include safety time and the sequence remains the same regardless of disruptions, so the result is robust. For the normal distribution we provide sufficient optimality conditions, precedence relationships that the optimal sequence must obey, and tight bounds.     

Second order (F, α, ρ, d)-convexity and duality in multiobjective programming

A class of second order (F, α, ρ, d)-convex functions and their generalizations is introduced. Using the assumptions on the functions involved, weak, strong and strict converse duality theorems are established for a second order Mond-Weir type multiobjective dual.     

Influence of doping BiYO3 and Nb2O5 on PTCR characteristics of BaTiO3 thermistor ceramics

Nb₂O₅-doped (1 − x)Ba_0.96Ca_0.04TiO₃-xBiYO₃ (where x = 0.01, 0.02, 0.03 and 0.04) lead-free PTC thermistor ceramics were prepared by a conventional solid state reaction method. X-ray diffraction, scanning electron microscope, Agilent E4980A and resistivity-temperature measurement instrument, were used to characteristic the lattice distortion, microstructure, temperature dependence of permittivity and resitivity-temperature dependence. It was revealed that the tetragonality c/a of the perovskite lattice, the microstructure and the Curie temperature changed with the BiYO₃ content. In order to decrease the room temperature resistivity, the effect of Nb₂O₅ on the room temperature resistivity was also studied, and its optimal doping content was finally chosen as 0.2 mol%. The 0.97Ba_0.96Ca_0.04TiO₃-0.03BiYO₃-0.002Nb₂O₅ thermistor ceramic exhibited a low ρ_RT of 3.98 × 10³ Ω cm, a typical PTCR effect of ρ_max/ρ_min > 10³ and a T_c of 153 °C.     

Scheduling projects with limited number of preemptions

A recent paper (Ballestín F, Valls V, Quintanilla S. Preemption in resource-constrained project scheduling. European Journal of Operational Research 2008;189:1136–152) revealed the potential benefits of allowing one interruption when scheduling activities in a resource-constrained project. This paper further investigates the effect of interruption on project length in more general cases. To achieve this, a new model that covers most practical applications of discrete activity preemption is proposed and a metaheuristic algorithm for the problem is developed. Computational experiments on the standard j120 and j30 sets generated using ProGen study the difference in makespan between allowing m$m$ interruptions per activity, m=0,1,2,…$m=0,1,2,\dots$ . In the second part of the paper the usefulness of preemption in the presence of due dates is analysed.     

Online interval scheduling on a single machine with finite lookahead

We study an online weighted interval scheduling problem on a single machine, where all intervals have unit length and the objective is to maximize the total weight of all completed intervals. We investigate how the function of finite lookahead improves the competitivities of deterministic online heuristics, under both preemptive and non-preemptive models. The lookahead model studied in this paper is that an online heuristic is said to have a lookahead ability of LD if at any time point it is able to foresee all the intervals to be released within the next LD   units of time. We investigate both competitive online heuristics and lower bounds on the competitive ratio, with lookahead 0≤LD≤1$0\le \mathit{LD}\le 1$ under the preemptive model, and lookahead 0≤LD≤2$0\le \mathit{LD}\le 2$ under the non-preemptive model. A method to transform a preemptive lookahead online algorithm to a non-preemptive online algorithm with enhanced lookahead ability is also given. Computational tests are performed to compare the practical competitivities of the online heuristics with different lookahead abilities.     

Impact of sea ice on the retrieval of water-leaving reflectance, chlorophyll a concentration and inherent optical properties from satellite ocean color data

Two physical phenomena by which satellite remotely sensed ocean color data are contaminated by sea ice at high latitudes are described through simulations and observations: (1) the adjacency effect that occurs along sea ice margins and (2) the sub-pixel contamination by a small amount of sea ice within an ocean pixel. The signal at the top of the atmosphere (TOA) was simulated using the 6S radiative transfer code that allows modeling of the adjacency effect for various types of sea ice surrounding an open water area. In situ sea ice reflectance spectra used in the simulations were measured prior to and during the melt period as part of the 2004 Canadian Arctic Shelf Exchange Study (CASES). For sub-pixel contamination, the TOA signal was simulated for various surface reflectances obtained by linear mixture of both sea ice and water-leaving reflectances (ρ_w). The standard atmospheric correction algorithm was then applied to the simulated TOA spectra to retrieve ρ_w spectra from which chlorophyll a concentrations (CHL) and inherent optical properties (IOPs) were derived. The adjacency effect was associated with large errors (> 0.002) in the retrieval of ρ_w as far as 24 km from an ice edge in the blue part of the spectrum (443 nm). Therefore, for moderate to high CHL (> 0.5 mg m^− 3), any pixel located within a distance of ∼ 10-20 km from the ice edge were unreliable. It was also found necessary to consider the adjacency effect when the total absorption coefficient (a_t) was to be retrieved using a semi-analytical algorithm. a_t(443) was underestimated by more than 35% at a distance of 20 km from an ice edge for CHL > 0.5 mg m^− 3. The effect on the retrieval of the particle backscattering coefficient (b_bp) was important only for clear waters (CHL ∼ 0.05 mg m^− 3). In contrast, sub-pixel contamination by a small amount of sea ice produced systematic underestimation of ρ_w in the blue because of incorrect interpretation of enhanced reflectance in the near infrared that is attributed to higher concentrations of atmospheric aerosols. In general, sub-pixel contamination was found to result in overestimations of CHL and a_t, and underestimations of b_bp. A simple method was proposed to flag pixels contaminated by adjacency effect.     

Mathematical models for preemptive shop scheduling problems

More than half a century has passed since Bowman and Dantzig (1959) [13] and [14] introduced their models for preemptive shop scheduling problems. A more efficient model seems to be needed to address all the aspects involved in the problem. We introduce a new Integer Linear Programming (ILP) formulation as a new method for solving the preemptive Job Shop Scheduling Problem (pJSSP). The dimension of the new model, unlike those of the existing ones, depends solely on the number of jobs and machines irrespective of processing times. The proposed model is used as an optimal, two-phase approach. In phase one, the model is solved to obtain the start and completion times of each operation on each machine. In phase two, a simple algorithm in O(mn log n) steps is used to turn these times into a complete and optimal schedule. Different preemptive flow shop problems are studied as special cases of the pJSSP while some related properties are also discussed. Finally, the higher efficiency of the proposed model is verified both theoretically and computationally through its comparison with conventional methods commonly in use.     

Remote sensing of woody shrub cover in desert grasslands using MISR with a geometric-optical canopy reflectance model

A new method is described for the retrieval of fractional cover of large woody plants (shrubs) at the landscape scale using moderate resolution multi-angle remote sensing data from the Multiangle Imaging SpectroRadiometer (MISR) and a hybrid geometric-optical (GO) canopy reflectance model. Remote sensing from space is the only feasible method for regularly mapping woody shrub cover over large areas, an important application because extensive woody shrub encroachment into former grasslands has been seen in arid and semi-arid grasslands around the world during the last 150 years. The major difficulty in applying GO models in desert grasslands is the spatially dynamic nature of the combined soil and understory background reflectance: the background is important and cannot be modeled as either a Lambertian scatterer or by using a fixed bidirectional reflectance distribution function (BRDF). Candidate predictors of the background BRDF at the Sun-target-MISR angular sampling configurations included the volume scattering kernel weight from a Li-Ross BRDF model; diffuse brightness (ρ0) from the Modified Rahman-Pinty-Verstraete (MRPV) BRDF model; other Li-Ross kernel weights (isotropic, geometric); and MISR near-nadir bidirectional reflectance factors (BRFs) in the blue, green, and near infra-red bands. The best method was multiple regression on the weights of a kernel-driven model and MISR nadir camera blue, green, and near infra-red bidirectional reflectance factors. The results of forward modeling BRFs for a 5.25 km² area in the USDA, ARS Jornada Experimental Range using the Simple Geometric Model (SGM) with this background showed good agreement with the MISR data in both shape and magnitude, with only minor spatial discrepancies. The simulations were shown to be accurate in terms of both absolute value and reflectance anisotropy over all 9 MISR views and for a wide range of canopy configurations (r² = 0.78, RMSE = 0.013, N = 3969). Inversion of the SGM allowed estimation of fractional shrub cover with a root mean square error (RMSE) of 0.03 but a relatively weak correlation (r² = 0.19) with the reference data (shrub cover estimated from high resolution IKONOS panchromatic imagery). The map of retrieved fractional shrub cover was an approximate spatial match to the reference map. Deviations reflect the first-order approximation of the understory BRDF in the MISR viewing plane; errors in the shrub statistics; and the 12 month lag between the two data sets.     

Minimizing Total Completion Time Subject to Job Release Dates and Preemption Penalties

Extensive research has been devoted to preemptive scheduling. However, little attention has been paid to problems where a certain time penalty must be incurred if preemption is allowed. In this paper, we consider the single-machine scheduling problem of minimizing the total completion time subject to job release dates and preemption penalties, where each time a job is started, whether initially or after being preempted, a job-independent setup must take place. The problem is proved to be strongly NP-hard even if the setup time is one unit. We also study a natural extension of a greedy algorithm, which is optimal in the absence of preemption penalty. It is proved that the algorithm has a worst-case performance bound of 25/16, even when the maximum completion time, i.e., makespan, criterion is considered simultaneously.     

Mapping leaf chlorophyll and leaf area index using inverse and forward canopy reflectance modeling and SPOT reflectance data

Reflectance data in the green, red and near-infrared wavelength region were acquired by the SPOT high resolution visible and geometric imaging instruments for an agricultural area in Denmark (56°N, 9°E) for the purpose of estimating leaf chlorophyll content (C_ab) and green leaf area index (LAI). SPOT reflectance observations were atmospherically corrected using aerosol data from MODIS and profiles of air temperature, humidity and ozone from the Atmospheric Infrared Sounder (AIRS), and used as input for the inversion of a canopy reflectance model. Computationally efficient inversion schemes were developed for the retrieval of soil and land cover-specific parameters which were used to build multiple species and site dependent formulations relating the two biophysical properties of interest to vegetation indices or single spectral band reflectances. Subsequently, the family of model generated relationships, each a function of soil background and canopy characteristics, was employed for a fast pixel-wise mapping of C_ab and LAI.The biophysical parameter retrieval scheme is completely automated and image-based and solves for the soil background reflectance signal, leaf mesophyll structure, specific dry matter content, Markov clumping characteristics, C_ab and LAI without utilizing calibration measurements.Despite the high vulnerability of near-infrared reflectances (ρ_nir) to variations in background properties, an efficient correction for background influences and a strong sensitivity of ρ_nir to LAI, caused LAI-ρ_nir relationships to be very useful and preferable over LAI-NDVI relationships for LAI prediction when LAI > 2. Reflectances in the green waveband (ρ_green) were chosen for producing maps of C_ab.The application of LAI-NDVI, LAI-ρ_nir and C_ab-ρ_green relationships provided reliable quantitative estimates of C_ab and LAI for agricultural crops characterized by contrasting architectures and leaf biochemical constituents with overall root mean square deviations between estimates and in-situ measurements of 0.74 for LAI and 5.0 μg cm^− 2 for C_ab.The results of this study illustrate the non-uniqueness of spectral reflectance relationships and the potential of physically-based inverse and forward canopy reflectance modeling techniques for a reasonably fast and accurate retrieval of key biophysical parameters at regional scales.     

Enhancing a leaf radiative transfer model to estimate concentrations and in vivo specific absorption coefficients of total carotenoids and chlorophylls a and b from single-needle reflectance and transmittance

The relative concentrations of different pigments within a leaf have significant physiological and spectral consequences. Photosynthesis, light use efficiency, mass and energy exchange, and stress response are dependent on relationships among an ensemble of pigments. This ensemble also determines the visible characteristics of a leaf, which can be measured remotely and used to quantify leaf biochemistry and structure. But current remote sensing approaches are limited in their ability to resolve individual pigments. This paper focuses on the incorporation of three pigments—chlorophyll a, chlorophyll b, and total carotenoids—into the LIBERTY leaf radiative transfer model to better understand relationships between leaf biochemical, biophysical, and spectral properties.Pinus ponderosa and Pinus jeffreyi needles were collected from three sites in the California Sierra Nevada. Hemispheric single-leaf visible reflectance and transmittance and concentrations of chlorophylls a and b and total carotenoids of fresh needles were measured. These data were input to the enhanced LIBERTY model to estimate optical and biochemical properties of pine needles. The enhanced model successfully estimated reflectance (RMSE = 0.0255, BIAS = 0.00477, RMS%E = 16.7%), had variable success estimating transmittance (RMSE = 0.0442, BIAS = 0.0294, RMS%E = 181%), and generated very good estimates of carotenoid concentrations (RMSE = 2.48 µg/cm², BIAS = 0.143 µg/cm², RMS%E = 20.4%), good estimates of chlorophyll a concentrations (RMSE = 10.7 µg/cm², BIAS = − 0.992 µg/cm², RMS%E = 21.1%), and fair estimates of chlorophyll b concentrations (RMSE = 7.49 µg/cm², BIAS = − 2.12 µg/cm², RMS%E = 43.7%). Overall root mean squared errors of reflectance, transmittance, and pigment concentration estimates were lower for the three-pigment model than for the single-pigment model. The algorithm to estimate three in vivo specific absorption coefficients is robust, although estimated values are distorted by inconsistencies in model biophysics. The capacity to invert the model from single-leaf reflectance and transmittance was added to the model so it could be coupled with vegetation canopy models to estimate canopy biochemistry from remotely sensed data.     

Testing the potential of multi-spectral remote sensing for retrospectively estimating fire severity in African Savannahs

The remote sensing of fire severity is a noted goal in studies of forest and grassland wildfires. Experiments were conducted to discover and evaluate potential relationships between the characteristics of African savannah fires and post-fire surface spectral reflectance in the visible to shortwave infrared spectral region. Nine instrumented experimental fires were conducted in semi-arid woodland savannah of Chobe National Park (Botswana), where fire temperature (T_max) and duration (dt) were recorded using thermocouples positioned at different heights and locations. These variables, along with measures of fireline intensity (FLI), integrated temperature with time (T_sum) and biomass (and carbon/nitrogen) volatilised were compared to post-fire surface spectral reflectance. Statistically significant relationships were observed between (i) the fireline intensity and total nitrogen volatilised (r² = 0.54, n = 36, p < 0.001), (ii) integrated temperature (T_sum−μ) and total biomass combusted (r² = 0.72, n = 32, p < 0.001), and (iii) fire duration as measured at the top-of-grass sward thermocouple (dt_T) and total biomass combusted (r² = 0.74, n = 34, p < 0.001) and total nitrogen volatilised (r² = 0.73, n = 34, p < 0.001). The post-fire surface spectral reflectance was found to be related to dt and T_sum via a quadratic relationship that varied with wavelength. The use of visible and shortwave infrared band ratios produced statistically significant linear relationships with fire duration as measured by the top thermocouple (dt_T) (r² = 0.76, n = 34, p < 0.001) and the mean of T_sum (r² = 0.82, n = 34, p < 0.001). The results identify fire duration as a versatile measure that relates directly to the fire severity, and also illustrate the potential of spectrally-based fire severity measures. However, the results also point to difficulties when applying such spectrally-based techniques to Earth Observation satellite imagery, due to the small-scale variability noted on the ground. Results also indicate the potential for surface spectral reflectance to increase following higher severity fires, due to the laying down of high albedo white mineral ash. Most current techniques for mapping burned area rely on the general assumption that surface albedo decreases following a fire, and so if the image spatial resolution was high enough such methods may fail. Determination of the effect of spatial resolution on a sensor's ability to detect white ash was investigated using a validated optical mixture modelling approach. The most appropriate mixing model to use (linear or non-linear) was assessed using laboratory experiments. A linear mixing model was shown most appropriate, with results suggesting that sensors having spatial resolutions significantly higher than those of Landsat ETM+ will be required if patches of white ash are to be used to provide EO-derived information on the spatial variation of fire severity.     

Hyperspectral remote sensing of cyanobacterial pigments as indicators for cell populations and toxins in eutrophic lakes

The growth of mass populations of toxin-producing cyanobacteria is a serious concern for the ecological status of inland waterbodies and for human and animal health. In this study we examined the performance of four semi-analytical algorithms for the retrieval of chlorophyll a (Chl a) and phycocyanin (C-PC) from data acquired by the Compact Airborne Spectrographic Imager-2 (CASI-2) and the Airborne Imaging Spectrometer for Applications (AISA) Eagle sensor. The retrieval accuracies of the semi-analytical models were compared to those returned by optimally calibrated empirical band-ratio algorithms. The best-performing algorithm for the retrieval of Chl a was an empirical band-ratio model based on a quadratic function of the ratio of reflectance at 710 and 670 nm (R² = 0.832; RMSE = 29.8%). However, this model only provided a marginally better retrieval than the best semi-analytical algorithm. The best-performing model for the retrieval of C-PC was a semi-analytical nested band-ratio model (R² = 0.984; RMSE = 3.98 mg m⁻³). The concentrations of C-PC retrieved using the semi-analytical model were correlated with cyanobacterial cell numbers (R² = 0.380) and the particulate and total (particulate plus dissolved) pools of microcystins (R² = 0.858 and 0.896 respectively). Importantly, both the empirical and semi-analytical algorithms were able to retrieve the concentration of C-PC at cyanobacterial cell concentrations below current warning thresholds for cyanobacteria in waterbodies. This demonstrates the potential of remote sensing to contribute to early-warning detection and monitoring of cyanobacterial blooms for human health protection at regional and global scales.     

GNet: A generalized network model and its applications in qualitative spatial reasoning

A data model, named generalized network (GNet), is proposed to perform various network-tracing tasks, especially tracing conceptual proposition networks in qualitative spatial reasoning (QSR). The GNet model can be defined as a 6-tuple: (V, A, q, ⊕, ∼, L). By specifying each element in the 6-tuple, a GNet can function as a conventional network, or an activity on edge (AOE) network, etc. The algorithm for searching for the generalized optimum path weight (GOPW) between two vertices in a GNet is developed by extending the Bellman-Ford algorithm (EBFA). Based on the GNet model, this paper focuses on representing spatial knowledge, which consists of a set of binary relations. We present two applications of GNets, namely the RCC8 network and the hybrid RCC8 network involving cardinal direction relations. Both can be traced to infer new spatial knowledge using EBFA. The applications demonstrate that the GNet model provides a promising approach to dealing with proposition-based geospatial knowledge based on weak composition. We also point out that EBFA can check whether a network is algebraically closed, or path-consistent when the corresponding composition table is extensional.     

Minimization of earliness, tardiness and due date penalties on uniform parallel machines with identical jobs

We consider a problem of scheduling n identical nonpreemptive jobs with a common due date on m uniform parallel machines. The objective is to determine an optimal value of the due date and an optimal allocation of jobs onto machines so as to minimize a total cost function, which is the function of earliness, tardiness and due date values. For the problem under study, we establish a set of properties of an optimal solution and suggest a two-phase algorithm to tackle the problem. First, we limit the number of due dates one needs to consider in pursuit of optimality. Next, we provide a polynomial-time algorithm to build an optimal schedule for a fixed due date. The key result is an O(m² log m) algorithm that solves the main problem to optimality.Scope and purpose: To extend the existing research on cost minimization with earliness, tardiness and due date penalties to the case of uniform parallel machines.     

Capturing tree crown formation through implicit surface reconstruction using airborne lidar data

Forest structure data derived from lidar is being used in forest science and management for inventory analysis, biomass estimation, and wildlife habitat analysis. Regression analysis dominated previous approaches to the derivation of tree stem and crown parameters from lidar. The regression model for tree parameters is locally applied based on vertical lidar point density, the tree species involved, and stand structure in the specific research area. The results of this approach, therefore, are location-specific, limiting its applicability to other areas. For a more widely applicable approach to derive tree parameters, we developed an innovative method called ‘wrapped surface reconstruction’ that employs radial basis functions and an isosurface. Utilizing computer graphics, we capture the exact shape of an irregular tree crown of various tree species based on the lidar point cloud and visualize their exact crown formation in three-dimensional space. To validate the tree parameters given by our wrapped surface approach, survey-grade equipment (a total station) was used to measure the crown shape. Four vantage points were established for each of 55 trees to capture whole-tree crown profiles georeferenced with post-processed differential GPS points. The observed tree profiles were linearly interpolated to estimate crown volume. These fieldwork-generated profiles were compared with the wrapped surface to assess goodness of fit. For coniferous trees, the following tree crown parameters derived by the wrapped surface method were highly correlated (p < 0.05) with the total station-derived measurements: tree height (R² = 0.95), crown width (R² = 0.80), live crown base (R² = 0.92), height of the lowest branch (R² = 0.72), and crown volume (R² = 0.84). For deciduous trees, wrapped surface-derived parameters of tree height (R² = 0.96), crown width (R² = 0.75), live crown base (R² = 0.53), height of the lowest branch (R² = 0.51), and crown volume (R² = 0.89) were correlated with the total station-derived measurements. The wrapped surface technique is less susceptible to errors in estimation of tree parameters because of exact interpolation using the radial basis functions. The effect of diminished energy return causes the low correlation for lowest branches in deciduous trees (R² = 0.51), even though leaf-off lidar data was used. The wrapped surface provides fast and automated detection of micro-scale tree parameters for specific applications in areas such as tree physiology, fire modeling, and forest inventory.     

Innovations in science teacher education: Effects of integrating technology and team-teaching strategies

The purpose of this study was to integrate technology and team-teaching techniques into science teacher education method courses in order to explore the effects of such integration on preservice teachers. The participants included one instructor and a total of 42 preservice teachers. A technology team-teaching model (TTT) was designed in this study to restructure science method courses with technology. This study used a mixed-method design, incorporating both quantitative and qualitative techniques. The results revealed that there were significant differences in “designing an appropriate science topic to be taught with technology” and “integrating computer activities with appropriate pedagogy in classroom instruction” (F = 5.260, p < 0.05, and F = 10.260, p < 0.01, respectively). The results also showed that the TTT model could enhance the integration of science teaching theories and practice. Team-teaching technique facilitated the integration of technology in science lesson design and teaching practice, and enhanced friendship through interaction. The TTT model could better the science learning experience of preservice teachers and serve as useful reference for other teacher education institutes.     

Remote sensing of water quality in an Australian tropical freshwater impoundment using matrix inversion and MERIS images

The purpose of this study was to investigate how semi-analytical inversion techniques developed for the remote sensing of water quality parameters (chlorophyll a, tripton and coloured dissolved organic matter (CDOM)) in inland waters could be adapted or improved for application to Australian tropical and sub-tropical water bodies. The Matrix Inversion Method (MIM) with a semi-analytic model of the anisotropy of the in-water light field was applied to MERIS images of Burdekin Falls Dam, Australia, a tropical freshwater impoundment. Specific attention was required to improve the atmospheric correction of the MERIS data. The performance of the conventional three band exact solution of the MIM was compared to that of over-determined solutions that used constant and differential weighting for each sensor band.The results of the application of the MIM algorithm showed that the best weighting scheme had a mean chlorophyll a retrieval difference of 1.0 μgl^− 1, the three band direct matrix inversion scheme had a mean difference of 4.2 μgl^− 1 and the constant weight scheme had a mean difference of 5.5 μgl^− 1. For tripton, the best performed weighting scheme had a mean difference of 1.2 mgl^− 1, the three band scheme had a mean difference of 3.4 mgl^− 1 and the constant weight scheme had a mean difference of 1.8 mgl^− 1. For the CDOM retrieval, the mean difference was found to be 0.12 m^− 1 for the best performed weighting scheme, 0.25 m^− 1 for the three band scheme and 0.52 m^− 1 for the constant weight scheme. It was found that significant improvements in the accuracy and precision of retrieved water quality parameter values can be obtained by using differentially weighted, over-determined systems of equations, rather than exact solutions. These more reliable estimates of water quality parameters will allow water resource managers to improve their monitoring regimes.     

Brain mass estimation by head circumference and body mass methods in neonatal glycaemic modelling and control

Introduction

Hyperglycaemia is a common complication of stress and prematurity in extremely low-birth-weight infants. Model-based insulin therapy protocols have the ability to safely improve glycaemic control for this group. Estimating non-insulin-mediated brain glucose uptake by the central nervous system in these models is typically done using population-based body weight models, which may not be ideal.

Method

A head circumference-based model that separately treats small-for-gestational-age (SGA) and appropriate-for-gestational-age (AGA) infants is compared to a body weight model in a retrospective analysis of 48 patients with a median birth weight of 750 g and median gestational age of 25 weeks. Estimated brain mass, model-based insulin sensitivity (S_I) profiles, and projected glycaemic control outcomes are investigated. SGA infants (5) are also analyzed as a separate cohort.

Results

Across the entire cohort, estimated brain mass deviated by a median 10% between models, with a per-patient median difference in S_I of 3.5%. For the SGA group, brain mass deviation was 42%, and per-patient S_I deviation 13.7%. In virtual trials, 87–93% of recommended insulin rates were equal or slightly reduced (Δ < 0.16 mU/h) under the head circumference method, while glycaemic control outcomes showed little change.

Conclusion

The results suggest that body weight methods are not as accurate as head circumference methods. Head circumference-based estimates may offer improved modelling accuracy and a small reduction in insulin administration, particularly for SGA infants.