The use of MODIS data to derive acreage estimations for larger fields: A case study in the south-western Rostov region of Russia

Recent developments in remote sensing technology, in particular improved spatial and temporal resolution, open new possibilities for estimating crop acreage over larger areas. Remotely sensed data allow in some cases the estimation of crop acreage statistics independently of sub-national survey statistics, which are sometimes biased and incomplete. This work focuses on the use of MODIS data acquired in 2001/2002 over the Rostov Oblast in Russia, by the Azov Sea. The region is characterised by large agricultural fields of around 75 ha on average. This paper presents a methodology to estimate crop acreage using the MODIS 16-day composite NDVI product. Particular emphasis is placed on a good quality crop mask and a good quality validation dataset. In order to have a second dataset which can be used for cross-checking the MODIS classification a Landsat ETM time series for four different dates in the season of 2002 was acquired and classified. We attempted to distinguish five different crop types and achieved satisfactory and good results for winter crops. Three hundred and sixty fields were identified to be suitable for the training and validation of the MODIS classification using a maximum likelihood classification. A novel method based on a pure pixel field sampling is introduced. This novel method is compared with the traditional hard classification of mixed pixels and was found to be superior.     

A study on the polarimetric properties of various features using SIR-C data

Spaceborne Imaging Radar (SIR-C) data acquired over Gujarat, India in 1994 were processed and analysed using differnet techniques applicable to polarimetric SAR data such as polarization signatures, polarization index, decomposition of the signal and polarization phase difference and limited groundtruth data. It has been observed that multi-frequency polarimetric data enhances the potential of retrieving geo-physical parameters. The polarization signatures are found to vary with the nature of the target. Target decomposition of the returned signal will be useful for the classification of various features. Polarization Phase Difference (PPD) gives good information about the vegetation parameters.     

The Bayesian detection of discontinuities in a polynomial regression and its application to the cycle-slip problem

This paper deals with the problem of detecting and correcting cycle-slips in Global Navigation Satellite System (GNSS) phase data by exploiting the Bayesian theory. The method is here applied to undifferenced observations, because repairing cycle-slips already at this stage could be a useful pre-processing tool, especially for a network of permanent GNSS stations. If a dual frequency receiver is available, the cycle-slips can be easily detected by combining two phase observations or phase and range observations from a single satellite to a single receiver. These combinations, expressed in a distance unit form, are completely free from the geometry and depend only on the ionospheric effect, on the electronic biases and on the initial integer ambiguities; since these terms are expected to be smooth in time, at least in a short period, a cycle-slip in one or both the two carriers can be modelled as a discontinuity in a polynomial regression. The proposed method consists in applying the Bayesian theory to compute the marginal posterior distribution of the discontinuity epoch and to detect it as a maximum a posteriori (MAP) in a very accurate way. Concerning the cycle-slip correction, a couple of simultaneous integer slips in the two carriers is chosen by maximazing the conditional posterior distribution of the discontinuity amplitude given the detected epoch. Numerical experiments on simulated and real data show that the discontinuities with an amplitude 2 or 3 times larger than the noise standard deviation are successfully identified. This means that the Bayesian approach is able to detect and correct cycle-slips using undifferenced GNSS observations even if the slip occurs by one cycle. A comparison with the scientific software BERNESE 5.0 confirms the good performance of the proposed method, especially when data sampled at high frequency (e.g. every 1 s or every 5 s) are available.     

Geodetic intersection on the ellipsoid

Through each of two known points on the ellipsoid a geodesic is passing in a known azimuth. We solve the problem of intersection of the two geodesics. The solution for the latitude is obtained as a closed formula for the sphere plus a small correction, of the order of the eccentricity of the ellipsoid, which is determined by numerical integration. The solution is iterative. Once the latitude is obtained, the longitude is determined without iteration.     

Ocean tide loading (OTL) displacements from global and local grids: comparisons to GPS estimates over the shelf of Brittany, France

In this paper we examine OTL displacements detected by GPS stations of a dedicated campaign and validate ocean tide models. Our area of study is the continental shelf of Brittany and Cotentin in France. Brittany is one of the few places in the world where tides provoke loading displacements of ∼10–12 cm vertically and a few cm horizontally. Ocean tide models suffer from important discrepancies in this region. Seven global and regional ocean tide models were tested: FES2004 corrected for K2, TPXO.7.0, TPXO.6.2, GOT00.2, CSR4.0, NAO.99b and the most recent regional grids of the North East Atlantic (NEA2004). These gridded amplitudes and phases of ocean tides were convolved in order to get the predicted OTL displacements using two different algorithms. Data over a period of 3.5 months of 8 GPS campaign stations located on the north coast of Brittany are used, in order to evaluate the geographical distribution of the OTL effect. We have modified and implemented new algorithms in our GPS software, GINS 7.1. GPS OTL constituents are estimated based on 1-day batch solutions. We compare the observed GPS OTL constituents of M₂, S₂, N₂ and K₁ waves with the selected ocean tide models on global and regional grids. Large phase-lag and amplitude discrepancies over 20° and 1.5 cm in the vertical direction in the semi-diurnal band of M₂ between predictions and GPS/models are detected in the Bay of Mont St-Michel. From a least squares spectral analysis of the GPS time-series, significant harmonic peaks in the integer multiples of the orbital periods of the GPS satellites are observed, indicating the existence of multipath effects in the GPS OTL constituents. The GPS OTL observations agree best with FES2004, NEA2004, GOT00.2 and CSR4.0 tide models.     

Prediction of UT1–UTC,LOD and AAM χ<Subscript>3</Subscript> by combination of least-squares and multivariate stochastic methods

This article presents the application of a multivariate prediction technique for predicting universal time (UT1–UTC), length of day (LOD) and the axial component of atmospheric angular momentum (AAM χ ₃). The multivariate predictions of LOD and UT1–UTC are generated by means of the combination of (1) least-squares (LS) extrapolation of models for annual, semiannual, 18.6-year, 9.3-year oscillations and for the linear trend, and (2) multivariate autoregressive (MAR) stochastic prediction of LS residuals (LS + MAR). The MAR technique enables the use of the AAM χ ₃ time-series as the explanatory variable for the computation of LOD or UT1–UTC predictions. In order to evaluate the performance of this approach, two other prediction schemes are also applied: (1) LS extrapolation, (2) combination of LS extrapolation and univariate autoregressive (AR) prediction of LS residuals (LS + AR). The multivariate predictions of AAM χ ₃ data, however, are computed as a combination of the extrapolation of the LS model for annual and semiannual oscillations and the LS + MAR. The AAM χ ₃ predictions are also compared with LS extrapolation and LS + AR prediction. It is shown that the predictions of LOD and UT1–UTC based on LS + MAR taking into account the axial component of AAM are more accurate than the predictions of LOD and UT1–UTC based on LS extrapolation or on LS + AR. In particular, the UT1–UTC predictions based on LS + MAR during El Niño/La Niña events exhibit considerably smaller prediction errors than those calculated by means of LS or LS + AR. The AAM χ ₃ time-series is predicted using LS + MAR with higher accuracy than applying LS extrapolation itself in the case of medium-term predictions (up to 100 days in the future). However, the predictions of AAM χ ₃ reveal the best accuracy for LS + AR.     

Bringing GPS into harsh environments for fully automated deformation monitoring

Engineering projects that require deformation monitoring frequently utilize geodetic sensors to measure displacements of target points located in the deformation zone. In situations where control stations and targets are separated by a kilometer or more, GPS can offer higher precision position updates at more frequent intervals than can normally be achieved using total station technology. For large-scale deformation projects requiring the highest precision, it is therefore advisable to use a combination of the two sensors. In response to the need for high precision, continuous GPS position updates in harsh deformation monitoring environments, a software has been developed that employs triple-differenced carrier-phase measurements in a delayed-state Kalman filter. Two data sets were analyzed to test the capabilities of the software. In the first test, a GPS antenna was displaced using a translation stage to mimic slow deformation. In the second test, data collected at a large open pit mine were processed. It was shown that the delayed-state Kalman filter developed could detect millimeter-level displacements of a GPS antenna. The actual precision attained depends upon the amount of process noise infused at each epoch to accommodate the antenna displacements. Higher process noise values result in quicker detection times, but at the same time increase the noise in the solutions. A slow, 25 mm displacement was detected within 30 min of the full displacement with sigma values in E, N and U of ±10 mm or better. The same displacement could also be detected in less than 5 h with sigma values in E, N and U of ±5 mm or better. The software works best for detecting long period deformations (e.g., 20 mm per day or less) for which sigma values of 1–2 mm are attained in all three solution components. It was also shown that the triple-differenced carrier-phase observation can be used to significantly reduce the effects of residual tropospheric delay that would normally plague double-differenced observations in harsh GPS environments.