Motion Analysis in SAR Images of Unfocused Objects Using Time–Frequency Methods

Synthetic aperture radar (SAR) image formation processing assumes that the scene is stationary, and to focus an object, one coherently sums a large number of independent returns. Any target motion introduces phases that distort and/or translate the target's image. Target motion produces a smear primarily in the azimuth direction of the SAR image. Time-frequency (TF) modeling is used to analyze and correct the residual phase distortions. An interactive focusing algorithm based on TF modeling demonstrates how to correct the phase and to rapidly focus the mover. This is demonstrated on two watercraft observed in a SAR image. Then, two time-frequency representations (TFRs) are applied to estimate the motion parameters of the movers or refocus them or both. The first is the short-time Fourier transform, from which a velocity profile is constructed based on the length of the smear. The second TFR is the time-frequency distribution series, which is a robust derivative of the Wigner-Ville distribution that works well in this SAR environment. The smear is a modulated chirp, from which a velocity profile is plotted and the phase corrections are integrated to focus the movers. The relationship between these two methods is discussed. Both methods show good agreement on the example.     

Best prediction in linear models with mixed integer/real unknowns: theory and application

In this contribution, we extend the existing theory of minimum mean squared error prediction (best prediction). This extention is motivated by the desire to be able to deal with models in which the parameter vectors have real-valued and/or integer-valued entries. New classes of predictors are introduced, based on the principle of equivariance. Equivariant prediction is developed for the real-parameter case, the integer-parameter case, and for the mixed integer/real case. The best predictors within these classes are identified, and they are shown to have a better performance than best linear (unbiased) prediction. This holds true for the mean squared error performance, as well as for the error variance performance. We show that, in the context of linear model prediction, best predictors and best estimators come in pairs. We take advantage of this property by also identifying the corresponding best estimators. All of the best equivariant estimators are shown to have a better precision than the best linear unbiased estimator. Although no restrictions are placed on the probability distributions of the random vectors, the Gaussian case is derived separately. The best predictors are also compared with least-squares predictors, in particular with the integer-based least-squares predictor introduced in Teunissen (J Geodesy, in press, 2006).     

Evaluating a thermal image sharpening model over a mixed agricultural landscape in India

Fine spatial resolution (e.g., <300 m) thermal data are needed regularly to characterise the temporal pattern of surface moisture status, water stress, and to forecast agriculture drought and famine. However, current optical sensors do not provide frequent thermal data at a fine spatial resolution. The TsHARP model provides a possibility to generate fine spatial resolution thermal data from coarse spatial resolution (≥1 km) data on the basis of an anticipated inverse linear relationship between the normalised difference vegetation index (NDVI) at fine spatial resolution and land surface temperature at coarse spatial resolution. The current study utilised the TsHARP model over a mixed agricultural landscape in the northern part of India. Five variants of the model were analysed, including the original model, for their efficiency. Those five variants were the global model (original); the resolution-adjusted global model; the piecewise regression model; the stratified model; and the local model. The models were first evaluated using Advanced Space-borne Thermal Emission Reflection Radiometer (ASTER) thermal data (90 m) aggregated to the following spatial resolutions: 180 m, 270 m, 450 m, 630 m, 810 m and 990 m. Although sharpening was undertaken for spatial resolutions from 990 m to 90 m, root mean square error (RMSE) of <2 K could, on average, be achieved only for 990–270 m in the ASTER data. The RMSE of the sharpened images at 270 m, using ASTER data, from the global, resolution-adjusted global, piecewise regression, stratification and local models were 1.91, 1.89, 1.96, 1.91, 1.70 K, respectively. The global model, resolution-adjusted global model and local model yielded higher accuracy, and were applied to sharpen MODIS thermal data (1 km) to the target spatial resolutions. Aggregated ASTER thermal data were considered as a reference at the respective target spatial resolutions to assess the prediction results from MODIS data. The RMSE of the predicted sharpened image from MODIS using the global, resolution-adjusted global and local models at 250 m were 3.08, 2.92 and 1.98 K, respectively. The local model consistently led to more accurate sharpened predictions by comparison to other variants.     

Interferometry with ENVISAT wide swath ScanSAR data

The possibility to get efficient topographic mapping and monitoring of large-scale motions with ScanSAR interferometry has been demonstrated with the Shuttle Radar Topography Mission and RADARSAT mission. The Environmental Satellite Advanced Synthetic Aperture Radar (ASAR) sensor has been designed to provide enhanced capabilities for interferometric applications. Different types of interferometric products can be obtained by combining the various ASAR modes as stripmap synthetic aperture radar [image mode (IM)] and ScanSAR [wide swath (WS) mode]. This letter deals with the possibility to use WS data to get either mixed-mode (IM/WS) or ScanSAR mode (WS/WS) differential interferograms. The impact of digital elevation model localization errors on IM/WS interferograms and of scan pattern synchronization on WS/WS interferograms is investigated. Experimental results are encouraging and show that ASAR ScanSAR data can be routinely used for interferometric applications in both cases.     

The Global Impact of Clouds on the Production of MODIS Bidirectional Reflectance Model-Based Composites for Terrestrial Monitoring

A global data set of cloud occurrence probability derived from Moderate Resolution Imaging Spectroradiometer (MODIS) Terra and Aqua gridded daily data is analyzed to investigate the probability of obtaining at least a minimum number of cloud-free observations within various compositing periods. The probabilities derived from Terra and Aqua, with morning and afternoon overpass times, respectively, are similar and increase with compositing period. Compositing both Terra and Aqua observations results in considerably higher probabilities of obtaining a sufficient number of observations for bidirectional reflectance model-based compositing. Given that the only alternative to obtaining sufficient samples is to extend the observation period, which can cause significant problems when the surface state changes, it is concluded that using data from the two MODIS sensors provides the most effective way of generating composited products. Findings with respect to the availability of cloud-free composites when n-day composites are generated on a temporally overlapping daily rolling basis, i.e., every day, rather than every n-days, are also discussed for regional and global applications     

The Potential of the Maximum Cross-Correlation Technique to Estimate Surface Currents From Thermal AVHRR Global Area Coverage Data

Having already shown its potential of deriving the vector fields representing the ocean-surface advection from sequential 1.1-km-resolution local area coverage (LAC) Advanced Very High Resolution Radiometer (AVHRR) images, the maximum cross-correlation (MCC) technique here is applied to four 4.4-km-resolution global area coverage (GAC) AVHRR images. The resulting three vector fields are compared to the vector fields obtained from the LAC imagery corresponding to the same satellite passages. To quantify the reduction in accuracy inevitable when applying the method to the lower resolution imagery, the LAC vector fields were assumed to be error free. The deviation of the GAC vectors from the LAC vectors is expressed as percentage errors of the signal variance of meridional u and zonal v velocity components, and they are 16%/30%, respectively, for the best case and 62%/117% and 92%/111% for the other two cases. These results indicate that, in its present state, the GAC data do not allow the MCC technique to extract reliable current-vector information from it     

Map-Based Exploratory Evaluation of Non-Medical Determinants of Population Health

Multi‐criteria evaluation (MCE) and decision‐making are increasingly combined with interactive tools to assist users with visual thinking and exploring decision strategies. The interactive control of criterion combination rules and the simultaneous observation of geographic space and criterion space provide a means of investigating the sensitivity of the decision outcome to the decision‐maker's preferences. The Analytic Hierarchy Process (AHP) is an MCE method that has been successfully implemented in management processes including those addressed by Geographic Information Systems. In this paper, we present a map‐based, interactive AHP implementation, which provides a link between a well‐understood decision support method and exploratory geographic visualization. Using a case study with public health data for the Province of Ontario, Canada, we demonstrate that exploratory map use increases the effectiveness of the AHP‐based evaluation of population health.     

Software for validating parameters retrieved from satellite

A PC-based interactive software has been developed and presented here for validating geophysical data retrieved from satellite mounted sensors operating in visible, infrared and microwave frequencies. The program, coded in Visual Basic, is user interactive and runs on Windows-98 or higher platforms. The system prepares the database on a pre-selected Microsoft platform to enhance processing efficiency. Sub-setting option is also provided to reduce the processing time. Data retrieved from ‘Multi-channel Scanning Microwave Radiometer (MSMR) onboard the Indian satellites Oceansat-1 during 1999–2001 were validated using this software as a case study. The program has several added advantages over the conventional method of validation that involves strenuous efforts to incorporate subroutines to meet every minute requirement. Satellite-sea truth relationships on various space-time window combinations are determined and exhibited in matrix form to visualize the nature of correlation. User has the option to visualize the satellite-sea truth relationship through graphical representations before selecting optimum relationship for prediction.     

Plate kinematics of Nubia–Somalia using a combined DORIS and GPS solution

We have used up to 12 years of data to assess DORIS performance for geodynamics applications. We first examine the noise characteristics of the DORIS time-series of weekly station coordinates to derive realistic estimates of velocity uncertainties. We find that a combination of white and flicker noise best explains the DORIS time-series noise characteristics. Second, weekly solutions produced by the Institut Géographique National/Jet Propulsion Laboratory (IGN/JPL) DORIS Analysis Centre are combined to derive a global velocity field. This solution is combined with two independent GPS solutions, including 11 sites on Nubia and 5 on the Somalia plate. The combination indicates that DORIS horizontal velocities have an average accuracy of 3 mm/year, with best-determined sites having velocity accuracy better than 1 mm/year (one-sigma levels). Using our combined velocity field, we derive an updated plate kinematics model with a focus on the Nubia–Somalia area. Including DORIS data improves the precision of the angular velocity vector for Nubia by 15%. Our proposed model provides robust bounds on the maximum opening rates along the East African Rift (4.7–6.7 mm/year). It indicates opening rates 15 and 7% slower than values predicted by NUVEL-1A for the southern Atlantic Ocean and Indian Ocean, respectively. These differences are likely to arise from the fact that NUVEL-1A considered Africa as a single non-deforming plate, while here we use a more refined approach.