Automated Target Detection and Discrimination Using Constrained Kurtosis Maximization

Exploiting hyperspectral imagery without prior information is a challenge. Under this circumstance, unsupervised target detection becomes an anomaly detection problem. We propose an effective algorithm for target detection and discrimination based on the normalized fourth central moment named kurtosis, which can measure the flatness of a distribution. Small targets in hyperspectral imagery contribute to the tail of a distribution, thus making it heavier. The Gaussian distribution is completely determined by the first two order statistics and has zero kurtosis. Consequently, kurtosis measures the deviation of a distribution from the background and is suitable for anomaly/target detection. When imposing appropriate inequality constraints on the kurtosis to be maximized, the resulting constrained kurtosis maximization (CKM) algorithm will be able to quickly detect small targets with several projections. Compared to the widely used unconstrained kurtosis maximization algorithm, i.e., fast independent component analysis, the CKM algorithm may detect small targets with fewer projections and yield a slightly higher detection rate.     

A Microwave Tomographic Approach for Nondestructive Testing of Dielectric Coated Metallic Surfaces

A microwave imaging method for nondestructive testing of perfectly conducting surfaces beyond a layered media is presented. The method is an adaptation of the surface reconstruction approach by Yapar et al. to the present problem. It is based on the analytical continuation of the measured data to the surface under test through a special representation of the scattered field in terms of Fourier transform and Taylor expansion. Then the problem is reduced to the solution of a nonlinear equation which is solved iteratively via the Newton method and regularization in the least squares sense. Numerical simulations show that defects as small as lambda/500 can be recovered through the presented algorithm.     

A data-driven approach to local gravity field modelling using spherical radial basis functions

We propose a methodology for local gravity field modelling from gravity data using spherical radial basis functions. The methodology comprises two steps: in step 1, gravity data (gravity anomalies and/or gravity disturbances) are used to estimate the disturbing potential using least-squares techniques. The latter is represented as a linear combination of spherical radial basis functions (SRBFs). A data-adaptive strategy is used to select the optimal number, location, and depths of the SRBFs using generalized cross validation. Variance component estimation is used to determine the optimal regularization parameter and to properly weight the different data sets. In the second step, the gravimetric height anomalies are combined with observed differences between global positioning system (GPS) ellipsoidal heights and normal heights. The data combination is written as the solution of a Cauchy boundary-value problem for the Laplace equation. This allows removal of the non-uniqueness of the problem of local gravity field modelling from terrestrial gravity data. At the same time, existing systematic distortions in the gravimetric and geometric height anomalies are also absorbed into the combination. The approach is used to compute a height reference surface for the Netherlands. The solution is compared with NLGEO2004, the official Dutch height reference surface, which has been computed using the same data but a Stokes-based approach with kernel modification and a geometric six-parameter “corrector surface” to fit the gravimetric solution to the GPS-levelling points. A direct comparison of both height reference surfaces shows an RMS difference of 0.6 cm; the maximum difference is 2.1 cm. A test at independent GPS-levelling control points, confirms that our solution is in no way inferior to NLGEO2004.     

Feature Extraction in Remote Sensing High-Dimensional Image Data

High-dimensional image data open new possibilities in remote sensing digital image classification, particularly when dealing with classes that are spectrally very similar. The main problem refers to the estimation of a large number of classifier's parameters. One possible solution to this problem consists in reducing the dimensionality of the original data without a significant loss of information. In this letter, a new approach to reduce data dimensionality is proposed. In the proposed methodology, each pixel's curve of spectral response is initially segmented, and the digital numbers (DNs) at each segment are replaced by a smaller number of statistics. In this letter, the proposed statistics are the mean and variance of the segment's DNs, which are supposed to carry information about the segment's position and shape, respectively. Tests were performed by using Airborne Visible/Infrared Imaging Spectrometer hyperspectral image data. The experiments have shown that this methodology is capable of providing very acceptable results, in addition of being computationally efficient     

High-Resolution Forecast Models of Water Vapor Over Mountains: Comparison With MERIS and Meteosat Data

Propagation delay due to variable tropospheric water vapor (WV) is one of the most intractable problems for radar interferometry, particularly over mountains. The WV field can be simulated by an atmospheric model, and the difference between the two fields is used to correct the radar interferogram. Here, we report our use of the U.K. Met Office Unified Model in a nested mode to produce high-resolution forecast fields for the 3-km-high Mount Etna volcano. The simulated precipitable-water field is validated against that retrieved from the Medium-Resolution Imaging Spectrometer (MERIS) radiometer on the Envisat satellite, which has a resolution of 300 m. Two case studies, one from winter (November 24, 2004) and one from summer (June 25, 2005), show that the mismatch between the model and the MERIS fields ( rms = 1.1 and 1.6 mm, respectively) is small. One of the main potential sources of error in the models is the timing of the WV field simulation. We show that long-wavelength upper tropospheric troughs of low WV could be identified in both the model output and Meteosat WV imagery for the November 24, 2004 case and used to choose the best time of model output.     

Testing satellite and ground thermal imaging of low-temperature fumarolic fields: The dormant Nisyros Volcano (Greece)

The Nisyros Volcano (Greece) was monitored by satellite and ground thermal imaging during the period 2000–2002. Three night-scheduled Landsat-7 ETM⁺ thermal (band 6) images of Nisyros Island were processed to obtain land surface temperature. Ground temperature data were also collected during one of the satellite overpasses. Processed results involving orthorectification and 3-D atmospheric correction clearly show the existence of a thermal anomaly inside the Nisyros Caldera. This anomaly is associated mainly with the largest hydrothermal craters and has land surface temperatures 5–10 °C warmer than its surroundings. The ground temperature generally increased by about 4 °C inside the main crater over the period 2000–2002. Ground thermal images of the hydrothermal Stephanos Crater were also collected in 2002 using a portable thermal infrared camera. These images were calibrated to ground temperature data and orthorectified. A difference of about 0–2 °C was observed between the ground thermal images and the ground temperature data. The overall study demonstrates that satellite remote sensing of low-temperature fumarolic fields within calderas can provide a reliable long-term monitoring tool of dormant volcanoes that have the potential to reactivate. Similarly, a portable thermo-imager can easily be deployed for real-time monitoring using telemetric data transfer. The operational costs for both systems are relatively low for an early warning system.     

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.     

Hyperspectral anomaly detection within the signal subspace

This letter describes the extension of signal subspace processing (SSP) to the arena of anomaly detection. In particular, we develop an SSP-based, local anomaly detector that exploits the rich information available in the multiple bands of a hyperspectral (HS) image. This SSP approach is based on signal processing considerations, and its entire formulation reduces to a straightforward (and intuitively pleasing) geometric and algebraic development. We extend the basic SSP concepts to the HS anomaly detection problem, develop an SSP HS anomaly detector, and evaluate this algorithm using multiple HS data files.