Influence of Atmospheric Water Vapour on IRS NIR Measurements for Detecting Vegetation Signal. Part I: A Simulation Study Using MODTRAN

Quantitative remote sensing involving accurate estimation of vegetation properties relies greatly on the measurements of the near infrared (NIR) channel because of unique interaction property between light and leaf. It is generally assumed that the NIR measurements are made in the absence of atmospheric absorption. However, relatively weak water vapour absorption features still persist in the NIR channel, which has bearing on the quantitative estimates of the vegetation properties and long-term data series. This paper reports the results of a study that was carried out to infer the possible influence of the atmospheric water vapour (WV) on the NIR measurements (0.77–0.86 μm) of Indian Remote Sensing (IRS) satellite sensors through radiative transfer simulations using MODTRAN model. The study also suggests and evaluates the alternate band-positions for the NIR channel to improve the IRS NIR measurements. It was observed that the water absorption features present around 0.810 μm reduces the WV transmission of IRS NIR channel from 1 to 0.91 when atmospheric WV content increased from 0 to 6 g/cm² and thus hampered the NIR reflectance by 14% as compared to reference signal. A significant improvement of the order of 6.5 to 12% in the NIR reflectance and 4.2 to 7% in NDVI was observed, when IRS NIR channel was split into NIR1 (0.775–0.805 μm) and NIR2 (0.845–0.875 μm) channels by avoiding the WV absorption features. The companion paper in this issue (Pandya et al. 2011) will support results of this simulation study through the EO1-Hyperion data analysis.     

Comparison of UAV and WorldView-2 imagery for mapping leaf area index of mangrove forest

Unmanned Aerial Vehicle (UAV) remote sensing has opened the door to new sources of data to effectively characterize vegetation metrics at very high spatial resolution and at flexible revisit frequencies. Successful estimation of the leaf area index (LAI) in precision agriculture with a UAV image has been reported in several studies. However, in most forests, the challenges associated with the interference from a complex background and a variety of vegetation species have hindered research using UAV images. To the best of our knowledge, very few studies have mapped the forest LAI with a UAV image. In addition, the drawbacks and advantages of estimating the forest LAI with UAV and satellite images at high spatial resolution remain a knowledge gap in existing literature. Therefore, this paper aims to map LAI in a mangrove forest with a complex background and a variety of vegetation species using a UAV image and compare it with a WorldView-2 image (WV2).In this study, three representative NDVIs, average NDVI (AvNDVI), vegetated specific NDVI (VsNDVI), and scaled NDVI (ScNDVI), were acquired with UAV and WV2 to predict the plot level (10 × 10 m) LAI. The results showed that AvNDVI achieved the highest accuracy for WV2 (R² = 0.778, RMSE = 0.424), whereas ScNDVI obtained the optimal accuracy for UAV (R² = 0.817, RMSE = 0.423). In addition, an overall comparison results of the WV2 and UAV derived LAIs indicated that UAV obtained a better accuracy than WV2 in the plots that were covered with homogeneous mangrove species or in the low LAI plots, which was because UAV can effectively eliminate the influence from the background and the vegetation species owing to its high spatial resolution. However, WV2 obtained a slightly higher accuracy than UAV in the plots covered with a variety of mangrove species, which was because the UAV sensor provides a negative spectral response function(SRF) than WV2 in terms of the mangrove LAI estimation.     

Estimation of LAI and above-ground biomass in deciduous forests: Western Ghats of Karnataka, India

This study demonstrates the potentials of IRS P6 LISS-IV high-resolution multispectral sensor (IGFOV  6 m)-based estimation of biomass in the deciduous forests in the Western Ghats of Karnataka, India. Regression equations describing the relationship between IRS P6 LISS-IV data-based vegetation index (NDVI) and field measured leaf area index (ELAI) and estimated above-ground biomass (EAGB) were derived. Remote sensing (RS) data-based leaf area index (PLAI) image is generated using regression equation based on NDVI and ELAI (r² = 0.68, p ≤ 0.05). RS-based above-ground biomass (PAGB) image was generated based on regression equation developed between PLAI and EAGB (r² = 0.63, p ≤ 0.05). The mean value of estimated above-ground biomass and RS-based above-ground biomass in the study area are 280(±72.5) and 297.6(±55.2) Mg ha⁻¹, respectively. The regression models generated in the study between NDVI and LAI; LAI and biomass can also help in generating spatial biomass map using RS data alone. LISS-IV-based estimation of biophysical parameters can also be used for the validation of various coarse resolution satellite products derived from the ground-based measurements alone.     

Measurement and Scaling of Carbon Dioxide (CO<Subscript>2</Subscript>) Exchanges in Wheat Using Flux-Tower and Remote Sensing

The present study investigates the characteristics of CO₂ exchange (photosynthesis and respiration) over agricultural site dominated by wheat crop and their relationship with ecosystem parameters derived from MODIS. Eddy covariance measurement of CO₂ and H₂O exchanges was carried out at 10 Hz interval and fluxes of CO₂ were computed at half-hourly time steps. The net ecosystem exchange (NEE) was partitioned into gross primary productivity (GPP) and ecosystem respiration (R _e) by taking difference between day-time NEE and respiration. Time-series of daily reflectance and surface temperature products at varying resolution (250–1000 m) were used to derive ecosystem variables (EVI, NDVI, LST). Diurnal pattern in Net ecosystem exchange reveals negative NEE during day-time representing CO₂ uptake and positive during night as release of CO₂. The amplitude of the diurnal variation in NEE increased as LAI crop growth advances and reached its peak around the anthesis stage. The mid-day uptake during this stage was around 1.15 mg CO₂ m⁻² s⁻¹ and night-time release was around 0.15 mg CO₂ m⁻² s⁻¹. Linear and non-linear least square regression procedures were employed to develop phenomenological models and empirical fits between flux tower based GPP and NEE with satellite derived variables and environmental parameters. Enhanced vegetation index was found significantly related to both GPP and NEE. However, NDVI showed little less significant relationship with both GPP and NEE. Furthemore, temperature-greenness (TG) model combining scaled EVI and LST was parameterized to estimate daily GPP over dominantly wheat crop site. (R ² = 0.77). Multi-variate analysis shows that inclusion of LST or air temperature with EVI marginally improves variance explained in daily NEE and GPP.     

Artificial Neural Network (ANN) Based Inversion of Benthic Substrate Bottom Type and Bathymetry in Optically Shallow Waters - Initial Model Results

Ocean-colour remote sensing in optically shallow waters is influenced by contribution from the water column depth as well as by the substrate type. Therefore, it is required to include the contribution from the water column and substrate bottom type for bathymetry estimation. In this report we demonstrate the use of Artificial Neural Network (ANN) based approach to spectrally distinguish various benthic bottom types and estimate depth of substrate bottom simultaneously in optically shallow waters. We have used in-water radiative transfer simulation modeling to generate simulated top-of-the-water column reflectance the four major benthic bottom types viz. sea grass, coral sand, green algae and red algae using Hydrolight simulation model. The simulated remote sensing reflectance, for the four benthic bottom types having benthic bottom depth up to 30 m were generated for moderately clear waters. A multi-layer perceptron (MLP) type neural network was trained using the simulated data. ANN based approach was used for classification of the benthic bottom type and simultaneous inversion of bathymetry. Simulated data was inverted to yield benthic bottom type classification with an accuracy of ~98% for the four benthic substrate types and the substrate depth were estimated with an error of 0% for sea grass, 1% for coral sand and 1–3% for green and red algae up to 25 m, whereas for substrate bottom deeper than 25 m depth the classification errors increased by 2–5% for three substrate bottom types except sea grass bottom type. The initial results are promising which needs validation using the in-situ measured remote sensing reflectance spectra for implementing further on satellite data.     

Utility of Hyperspectral Data for Potato Late Blight Disease Detection

The study was carried out to investigate the utility of hyperspectral reflectance data for potato late blight disease detection. The hyperspectral data was collected for potato crop at different level of disease infestation using hand-held spectroradiometer over the spectral range of 325–1075 nm. The data was averaged into 10-nm wide wavebands, resulting in 75 narrowbands. The reflectance curve was partitioned into five regions, viz. 400–500 nm, 520–590 nm, 620–680 nm, 770–860 nm and 920–1050 nm. The notable differences in healthy and diseased potato plants were noticed in 770–860 nm and 920–1050 nm range. Vegetation indices, namely NDVI, SR, SAVI and red edge were calculated using reflectance values. The differences between the vegetation indices for plants at different levels of disease infestation were found highly significant. The optimal hyperspectral wavebands to discriminate the healthy plants from disease infested plants were 540, 610, 620, 700, 710, 730, 780 and 1040 nm whereas upto 25% infestation could be discriminated using reflectance at 710, 720 and 750 nm.     

Formulation of Time Series Vegetation Index from Indian Geostationary Satellite and Comparison with Global Product

To study impact of climate change on vegetation time series vegetation index has a vital role to know the behaviour of vegetation dynamics over a time period. INSAT 3A CCD (Charged Couple Device) is the only geostationary sensor to acquire regular coverage of Asia continent at 1 km × 1 km spatial resolution with high temporal frequency (half-an-hour). A formulation of surface reflectances in red, near infrared (NIR), short wave infrared (SWIR) and NDVI from INSAT 3A CCD has been defined and integrated in the operational chain. The atmospheric correction of at-sensor reflectances using SMAC (Simple Model for Atmospheric Correction) model improved the NDVI by 5–40% and also increased its dynamic range. The temporal dynamics of 16-day NDVI composite at 0500 GMT for a growing year (June 2008–March 2009) showed matching profiles with reference to global products (MODIS TERRA) over known land targets. The root mean square deviation (RMSD) between the two was 0.14 with correlation coefficient (r) 0.84 from 200 paired datasets. This inter-sensor cross-correlation would help in NDVI calibration to add continuity in long term NDVI database for climate change studies.     

Study of the optimum time span for distinguishing <Emphasis Type="Italic">Avena ludoviciana</Emphasis> from wheat crop based on their spectral reflectance characteristics

The study to establish the optimum time span for distinguishing Avena ludoviciana from wheat crop based on their spectral signatures was carried out at Student’s Research Farm, Department of Agronomy during 2006–07 and 2007–08. The experimental sites during both the seasons were sandy loam in texture, with normal soil reaction and electrical conductivity, low in organic carbon and available nitrogen and medium in available phosphorus and potassium. The experiment was laid out in randomized block design with four replications and consisting of twelve treatments comprising 0, 10, 15, 25, 50, 75, 100, 125, 150, 200, 250 plants m⁻² and a pure Avena ludoviciana plot (T_max). The results revealed that in all the treatments irrespective of wheat and weeds, the red reflectance (%) value decreased from 34 to 95 DAS (days after sowing) in 2006–07 and 45 DAS to 100 DAS during 2007–08, and thereafter a sharp increase was observed in all the treatments. This trend might be due to increased chlorophyll index after 34 DAS as red reflectance was reduced by chlorophyll absorption. Among all the treatments, T_max (Pure Avena ludoviciana plot) had the highest red reflectance and T₀ (Pure wheat plot) had a lowest value of red reflectance during both the years. The highest value of IR reflectance was obtained at 95 DAS (2006–07) and 70 DAS (2007–08) in all the treatments. IR reflectance of wheat crop ranged between 24.61 and 61.21 per cent during 2006–07 and 27.33 and 67.3 per cent during 2007–08. However, IR reflectance values declined after 95 DAS and 70 DAS up to harvesting during 2006–07 and 2007–08. This lower reflectance may have been due to the onset of senescence. The highest RR and NDVI values were recorded under pure wheat treatment and minimum under pure weed plots. This may be due to dark green colour and better vigor of the wheat as compared to Avena ludoviciana. It was observed that by using RR and NDVI, pure wheat can be distinguished from pure populations of Avena ludoviciana after 34 DAS and different levels of weed populations can be discriminated amongst themselves from 68 DAS up to 107 DAS during both the years of investigation.     

Retrieval of wheat leaf area index from AWiFS multispectral data using canopy radiative transfer simulation

Accurate representation of leaf area index (LAI) from high resolution satellite observations is obligatory for various modelling exercises and predicting the precise farm productivity. Present study compared the two retrieval approach based on canopy radiative transfer (CRT) method and empirical method using four vegetation indices (VI) (e.g. NDVI, NDWI, RVI and GNDVI) to estimate the wheat LAI. Reflectance observations available at very high (56 m) spatial resolution from Advanced Wide-Field Sensor (AWiFS) sensor onboard Indian Remote Sensing (IRS) P6, Resourcesat-1 satellite was used in this study. This study was performed over two different wheat growing regions, situated in different agro-climatic settings/environments: Trans-Gangetic Plain Region (TGPR) and Central Plateau and Hill Region (CPHR). Forward simulation of canopy reflectances in four AWiFS bands viz. green (0.52–0.59 μm), red (0.62–0.68 μm), NIR (0.77–0.86 μm) and SWIR (1.55–1.70 μm) were carried out to generate the look up table (LUT) using CRT model PROSAIL from all combinations of canopy intrinsic variables. An inversion technique based on minimization of cost function was used to retrieve LAI from LUT and observed AWiFS surface reflectances. Two consecutive wheat growing seasons (November 2005–March 2006 and November 2006–March 2007) datasets were used in this study. The empirical models were developed from first season data and second growing season data used for validation. Among all the models, LAI-NDVI empirical model showed the least RMSE (root mean square error) of 0.54 and 0.51 in both agro-climatic regions respectively. The comparison of PROSAIL retrieved LAI with in situ measurements of 2006–2007 over the two agro-climatic regions produced substantially less RMSE of 0.34 and 0.41 having more R² of 0.91 and 0.95 for TGPR and CPHR respectively in comparison to empirical models. Moreover, CRT retrieved LAI had less value of errors in all the LAI classes contrary to empirical estimates. The PROSAIL based retrieval has potential for operational implementation to determine the regional crop LAI and can be extendible to other regions after rigorous validation exercise.     

Assessment of RapidEye vegetation indices for estimation of leaf area index and biomass in corn and soybean crops

Leaf area index (LAI) and biomass are important indicators of crop development and the availability of this information during the growing season can support farmer decision making processes. This study demonstrates the applicability of RapidEye multi-spectral data for estimation of LAI and biomass of two crop types (corn and soybean) with different canopy structure, leaf structure and photosynthetic pathways. The advantages of Rapid Eye in terms of increased temporal resolution (∼daily), high spatial resolution (∼5 m) and enhanced spectral information (includes red-edge band) are explored as an individual sensor and as part of a multi-sensor constellation. Seven vegetation indices based on combinations of reflectance in green, red, red-edge and near infrared bands were derived from RapidEye imagery between 2011 and 2013. LAI and biomass data were collected during the same period for calibration and validation of the relationships between vegetation indices and LAI and dry above-ground biomass. Most indices showed sensitivity to LAI from emergence to 8 m²/m². The normalized difference vegetation index (NDVI), the red-edge NDVI and the green NDVI were insensitive to crop type and had coefficients of variations (CV) ranging between 19 and 27%; and coefficients of determination ranging between 86 and 88%. The NDVI performed best for the estimation of dry leaf biomass (CV = 27% and r² = 090) and was also insensitive to crop type. The red-edge indices did not show any significant improvement in LAI and biomass estimation over traditional multispectral indices. Cumulative vegetation indices showed strong performance for estimation of total dry above-ground biomass, especially for corn (CV ≤ 20%). This study demonstrated that continuous crop LAI monitoring over time and space at the field level can be achieved using a combination of RapidEye, Landsat and SPOT data and sensor-dependant best-fit functions. This approach eliminates/reduces the need for reflectance resampling, VIs inter-calibration and spatial resampling.     

Remote Sensing Study of Sittampundi Anorthosite Complex,India

In the present study, The Landsat 7 ETM satellite data was collected for the Sittampundi anorthosites complex and digital image analysis was carried out. The anorthositic rocks available at Sittampundi complex is considered as an equivalent of lunar highland rocks. Hence, a remote sensing study comprises of image analysis and spectral profile analysis was carried out. The satellite data was digitally processed and generated various outputs like band combinations, color composites, stretched outputs, and PCA. The suitable processed outputs were identified for delineating the anorthosite complex. The diagnostic absorption features of reflectance spectra are the sensitive indicators of mineralogy and chemical composition of rocks, which are interest to the planetary scientists. The spectral profile of Landsat ETM plotted for pure and mixed anorthosite pixels and compared with the field and lab reflectance spectra. The percentages of image spectra vary from 30% to 60% for Sittampundi anorthosite. The spectral bands 2, 4 and 6 have low reflectance and bands 3 and 5 have high reflectance. The spectral range of bands 2,3,4,5 and 6 are 525 nm–605 nm, 630 nm–690 nm, 750 nm–900 nm, 1550 nm–1750 nm and 10400 nm–12500 nm respectively. The field spectral curve has weak absorptions at 650 nm and 1000 nm due to the iron transition absorption and low ca- pyroxene respectively available in the anorthosite, matching with the image spectra. However, hyperspectal image with narrow bandwidth could be more useful in selecting the suitable spectrum for remotely mapping the anorthosite region, as equivalent test site for lunar highland region.     

Spectral wheat growth profile in Punjab using IRS WiFS data

A functional form of crop spectral profile suggested by Badhwar was applied to district-wise wheat Normalised Difference Vegetation Index (NDVI) values relatively normalised by Pseudo-Invariant Feature (urban and built-up) NDVI values, derived from Wide Field Sensor (WiFS) onboard Indian Remote Sensing Satellites (IRS) for 17 dates during 1999–2000 rabi season. The goodness of overall profile fitting and the three basic parameters i.e., crop emergence date (T_o), and crop specific parameters (a and P) was found to be statistically significant. While a corresponds to profile progressive growth rate, β corresponds to profile decay rate. A comparison with earlier studies in Punjab using NOAA-AVHRR indicated improvement in relation between peak NDVI and wheat yield. The estimated time of spectral emergence and profile-derived peak NDVI follow the observed behaviour of shortened crop pre-anthesis period with delayed sowing.     

Retrieval of leaf area index in alpine wetlands using a two-layer canopy reflectance model

In this paper, we focused on the retrieval of the LAI in an alpine wetland located in western part of China in late August and early July 2011. A two-layer canopy reflectance model (ACRM) was used to establish the relationships between the LAI and the reflectance of near-infrared (NIR) and red (RED) wavebands. The reflectance data were derived from Landsat TM L1T product and the Terra and Aqua MODIS 16-day and 8-day composite reflectance products (MOD/MYD09) at 250 m resolution. Due to the lack of the information about some major input parameters for ACRM, which are sensitive to model outputs in the reflectance of NIR and RED wavebands, the inverse problem was ill-posed. To overcome this problem, a method of increasing the sensitivity of the LAI while reducing the influence of other model free parameters based on the study of free parameters’ sensitivity to the ACRM outputs and the region’s features was studied. The area of interest was divided into two parts using the approximately statistic normalized difference vegetation index (NDVI) value around 0.5. One part was sparse vegetation (0.1 < NDVI < 0.5), which is more sensitive to soil background effects and less sensitive to the canopy biophysical and biochemical variables. The other part was dense vegetation (0.5 ≤ NDVI < 1.0), which is less sensitive to soil background effects and more sensitive to plant canopies and leaf parameters. Then, the relationships of ρ_nir–LAI and ρ_red–LAI were established using a look-up table algorithm for the two parts. Furthermore, a regularization technique for fast pixel-wise retrieval was introduced to reduce the elements of LUT sets while maintaining a relatively high accuracy. The results were very promising compared to the field measured LAI values that the correlation (R²) of the measured LAI values and retrieved LAI values reached 0.95, and the root-mean-square deviation (RMSD) was 0.33 for late August, 2011, while the R² reached 0.82 and RMSD was 0.25 for early July 2011.     

Atmospheric Correction Over Coastal Turbid Waters of Bay of Bengal Using OCEANSAT-I Ocean Colour Monitor (OCM) Data

Space-borne ocean-colour remote sensor-detected radiance is heavily contaminated by solar radiation backscattered by the atmospheric air molecules and aerosols. Hence, the first step in ocean-colour data processing is the removal of this atmospheric contribution from the sensor-detected radiance to enable detection of optically active oceanic constituents e.g. chlorophyll-a, suspended sediment etc. In standard atmospheric correction procedure for OCEANSAT-1 Ocean Colour Monitor (OCM) data, NIR bands centered at 765 and 865 nm wavelengths were used for aerosol characterization. Due to high absorption by water molecules, ocean surface in these two wavelengths acts as dark background, therefore, sensor detected radiance can be assumed to have major contribution from atmospheric scattering. For coastal turbid waters this assumption of dark surface fails due to the presence of highly scattering sediments which causes sufficient water-leaving radiance in NIR bands and lead to over-estimation of aerosol radiance resulting in negative water leaving radiance for λ < 700 nm. In the present study, for the turbid coastal waters in the northern Bay of Bengal, the concept of spatial homogeneity of aerosol and water leaving reflectance has been applied to perform atmospheric correction of OCAEANSAT-1 OCM data. The results of the turbid water atmospheric correction have also been validated using in-situ measured water-leaving radiance. Comparison of satellite derived water-leaving radiance for five coastal stations with in-situ measured radiance spectra, indicates an improvement over the standard atmospheric correction algorithm giving physically realistic and positive values. Root Mean Square Error (RMSE) between the in-situ measured and satellite derived water leaving radiance for wavelengths 412 nm, 443 nm, 490 nm, 512 nm and 555 nm was found to be 1.11, 0.718, 0.575, 0.611 and 0.651%, respectively, using standard atmospheric correction procedure. By the use of spatial homogeneity concept, this error was reduced to 0.125, 0.173, 0.176, 0.225, and 0.290 and the correlation coefficient arrived at 0.945, which is an improvement over the standard atmospheric correction procedure.     

Land surface temperature retrieval and its validation using NOAA AVHRR thermal data

The retrieval of land (soil-vegetation complex) surface temperature (LST) was carried out over semi-arid mixed agriculture landscape of Gujarat using thermal bands (channel 4 and 5) and ground emissivity from atmospherically corrected NDVI of NOAA AVHRR LAC images. The atmospheric correction of Visible and NIR band reflectance was done using SMAC model. The LST computed from split-window method and subsequently corrected with fractional vegetation cover were then compared with near synchronous ground observations of soil and air temperatures made during 13–17 January and April, 1997 at five Land Surface Processes Experiment (LASPEX) sites of Anand, Sanand, Derol, Arnej and Khandha covering 100 km x 100 km. The fractional vegetation cover corrected LST at noon hrs. varied from 301.6 – 311.9K in January and from 315.8 – 325.6K in April. The LST_corr were found to lie in the mid way between AT and ST during January. But in April, LST were found to be more close to ST which may be due to relatively poor vegetation growth as indicated by lower NDVI values in April indicating more contribution to LST from exposed soil surface.     

Inversion of PROSAIL Model for Retrieval of Plant Biophysical Parameters

The current development of satellite technology particularly in the sensors like POLDER and MISR, has emphasized more on directional reflectance measurements (i.e. spectral reflectance of the target measured from different view zenith and azimuth angles) of the earth surface features mainly the vegetation for retrieval of biophysical parameters at regional scale using radiative transfer models. This approach being physical process based and uses directional reflectance measurement has been found to better and more reliable compared to the conventional statistical approach used till date and takes care of anisotropic nature (i.e. reflectance from the target is different if measured from different view angles) of the target. Keeping this in view a field experiment was conducted in mustard crop to evaluate the radiative transfer model for biophysical parameter retrieval through its inversion with the objectives set as (i) to relate canopy biophysical parameters and geometry to its bidirectional reflectance, (ii) to evaluate a canopy reflectance model to best represent the radiative transfer within the canopy for its inversion and (iii) to retrieve crop biophysical parameters through inversion of the model. Two varieties of the mustard crop (Brassica juncea L) were grown with two nitrogen treatments. The bidirectional reflectance data obtained at 5 nm interval for a range of 400–1100 nm were integrated to IRS LISS–II sensor’s four band values using Newton Cotes Integration technique. Biophysical parameters like leaf area index, leaf chlorophyll content, leaf length, plant height and average leaf inclination angle, biomass etc were estimated synchronizing with the bi-directional reflectance measurements. Radiative transfer model PROSAIL model was validated and its inversion was done to retrieve LAI and ALA. Look Up Table (LUT) of Bidirectional reflectance distribution function (BRDF) was prepared simulating through PROSAIL model varying only LAI (0.2 interval from 1.2 to 5.4 ) and ALA (5° interval from 40° to 55°) parameters and inversion was done using a merit function and numerical optimization technique given by Press et al. (1986). The derived LAI and ALA values from inversion were well matched with observed one with RMSE 0.521 and 5.57, respectively.     

Sensitivity of spectral reflectance values to different burn and vegetation ratios: A multi-scale approach applied in a fire affected area

The aim of our study was to explore the spectral properties of fire-scorched (burned) and non fire-scorched (vegetation) areas, as well as areas with different burn/vegetation ratios, using a multisource multiresolution satellite data set. A case study was undertaken following a very destructive wildfire that occurred in Parnitha, Greece, July 2007, for which we acquired satellite images from LANDSAT, ASTER, and IKONOS. Additionally, we created spatially degraded satellite data over a range of coarser resolutions using resampling techniques. The panchromatic (1 m) and multispectral component (4 m) of IKONOS were merged using the Gram-Schmidt spectral sharpening method. This very high-resolution imagery served as the basis to estimate the cover percentage of burned areas, bare land and vegetation at pixel level, by applying the maximum likelihood classification algorithm. Finally, multiple linear regression models were fit to estimate each land-cover fraction as a function of surface reflectance values of the original and the spatially degraded satellite images.The main findings of our research were: (a) the Near Infrared (NIR) and Short-wave Infrared (SWIR) are the most important channels to estimate the percentage of burned area, whereas the NIR and red channels are the most important to estimate the percentage of vegetation in fire-affected areas; (b) when the bi-spectral space consists only of NIR and SWIR, then the NIR ground reflectance value plays a more significant role in estimating the percent of burned areas, and the SWIR appears to be more important in estimating the percent of vegetation; and (c) semi-burned areas comprising 45–55% burned area and 45–55% vegetation are spectrally closer to burned areas in the NIR channel, whereas those areas are spectrally closer to vegetation in the SWIR channel. These findings, at least partially, are attributed to the fact that: (i) completely burned pixels present low variance in the NIR and high variance in the SWIR, whereas the opposite is observed in completely vegetated areas where higher variance is observed in the NIR and lower variance in the SWIR, and (ii) bare land modifies the spectral signal of burned areas more than the spectral signal of vegetated areas in the NIR, while the opposite is observed in SWIR region of the spectrum where the bare land modifies the spectral signal of vegetation more than the burned areas because the bare land and the vegetation are spectrally more similar in the NIR, and the bare land and burned areas are spectrally more similar in the SWIR.     

Testing the discrimination and detection limits of WorldView-2 imagery on a challenging invasive plant target

Invasive plants pose significant threats to biodiversity and ecosystem function globally, leading to costly monitoring and management effort. While remote sensing promises cost-effective, robust and repeatable monitoring tools to support intervention, it has been largely restricted to airborne platforms that have higher spatial and spectral resolutions, but which lack the coverage and versatility of satellite-based platforms. This study tests the ability of the WorldView-2 (WV2) eight-band satellite sensor for detecting the invasive shrub mesquite (Prosopis spp.) in the north-west Pilbara region of Australia. Detectability was challenged by the target taxa being largely defoliated by a leaf-tying biological control agent (Gelechiidae: Evippe sp. #1) and the presence of other shrubs and trees. Variable importance in the projection (VIP) scores identified bands offering greatest capacity for discrimination were those covering the near-infrared, red, and red-edge wavelengths. Wavelengths between 400 nm and 630 nm (coastal blue, blue, green, yellow) were not useful for species level discrimination in this case. Classification accuracy was tested on three band sets (simulated standard multispectral, all bands, and bands with VIP scores ≥1). Overall accuracies were comparable amongst all band-sets (Kappa = 0.71–0.77). However, mesquite omission rates were unacceptably high (21.3%) when using all eight bands relative to the simulated standard multispectral band-set (9.5%) and the band-set informed by VIP scores (11.9%). An incremental cover evaluation on the latter identified most omissions to be for objects <16 m². Mesquite omissions reduced to 2.6% and overall accuracy significantly improved (Kappa = 0.88) when these objects were left out of the confusion matrix calculations. Very high mapping accuracy of objects >16 m² allows application for mapping mesquite shrubs and coalesced stands, the former not previously possible, even with 3 m resolution hyperspectral imagery. WV2 imagery offers excellent portability potential for detecting other species where spectral/spatial resolution or coverage has been an impediment. New generation satellite sensors are removing barriers previously preventing widespread adoption of remote sensing technologies in natural resource management.     

The total optimal search criterion in solving the mixed integer linear model with GNSS carrier phase observations

Existing algorithms for GPS ambiguity determination can be classified into three categories, i.e. ambiguity resolution in the measurement domain, the coordinate domain and the ambiguity domain. There are many techniques available for searching the ambiguity domain, such as FARA (Frei and Beutler in Manuscr Geod 15(4):325–356, 1990), LSAST (Hatch in Proceedings of KIS’90, Banff, Canada, pp 299–308, 1990), the modified Cholesky decomposition method (Euler and Landau in Proceedings of the sixth international geodetic symposium on satellite positioning, Columbus, Ohio, pp 650–659, 1992), LAMBDA (Teunissen in Invited lecture, section IV theory and methodology, IAG general meeting, Beijing, China, 1993), FASF (Chen and Lachapelle in J Inst Navig 42(2):371–390, 1995) and modified LLL Algorithm (Grafarend in GPS Solut 4(2):31–44, 2000; Lou and Grafarend in Zeitschrift für Vermessungswesen 3:203–210, 2003). The widely applied LAMBDA method is based on the Least Squares Ambiguity Search (LSAS) criterion and employs an effective decorrelation technique in addition. G. Xu (J Glob Position Syst 1(2):121–131, 2002) proposed also a new general criterion together with its equivalent objective function for ambiguity searching that can be carried out in the coordinate domain, the ambiguity domain or both. Xu’s objective function differs from the LSAS function, leading to different numerical results. The cause of this difference is identified in this contribution and corrected. After correction, the Xu’s approach and the one implied in LAMBDA are identical. We have developed a total optimal search criterion for the mixed integer linear model resolving integer ambiguities in both coordinate and ambiguity domain, and derived the orthogonal decomposition of the objective function and the related minimum expressions algebraically and geometrically. This criterion is verified with real GPS phase data. The theoretical and numerical results show that (1) the LSAS objective function can be derived from the total optimal search criterion with the constraint on the fixed integer ambiguity parameters, and (2) Xu’s derivation of the equivalent objective function was incorrect, leading to an incorrect search procedure. The effects of the total optimal criterion on GPS carrier phase data processing are discussed and its practical implementation is also proposed.     

Estimating chlorophyll with thermal and broadband multispectral high resolution imagery from an unmanned aerial system using relevance vector machines for precision agriculture

Precision agriculture requires high-resolution information to enable greater precision in the management of inputs to production. Actionable information about crop and field status must be acquired at high spatial resolution and at a temporal frequency appropriate for timely responses. In this study, high spatial resolution imagery was obtained through the use of a small, unmanned aerial system called AggieAir^TM. Simultaneously with the AggieAir flights, intensive ground sampling for plant chlorophyll was conducted at precisely determined locations. This study reports the application of a relevance vector machine coupled with cross validation and backward elimination to a dataset composed of reflectance from high-resolution multi-spectral imagery (VIS–NIR), thermal infrared imagery, and vegetative indices, in conjunction with in situ SPAD measurements from which chlorophyll concentrations were derived, to estimate chlorophyll concentration from remotely sensed data at 15-cm resolution. The results indicate that a relevance vector machine with a thin plate spline kernel type and kernel width of 5.4, having LAI, NDVI, thermal and red bands as the selected set of inputs, can be used to spatially estimate chlorophyll concentration with a root-mean-squared-error of 5.31 μg cm⁻², efficiency of 0.76, and 9 relevance vectors.