1. Introduction
Hyperspectral remote sensing for earth observation has attracted increasing interest in recent years. In addition to the ground cover’s spatial distribution, hyperspectral imaging sensors can synchronously measure the reflectance of materials in hundreds of narrow and contiguous bands, spanning over the visible-to-infrared and even a wider range of the electromagnetic spectrum [
1,
2]. Therefore, each pixel of hyperspectral imagery (HSI) can be presented as a continuous spectral curve. HSI has been extensively applied in various fields, such as mineral exploration, environmental protection, and military affairs [
3,
4,
5].
Target detection (TD) is one of the most important HSI processing technologies [
6,
7]. Different materials usually have distinct wavelength-dependent electromagnetic scattering and radiation characteristics, which essentially enables TD for HSI. In fact, TD is about solving a supervised dichotomous problem where pixels are labeled the target or background according to their similarity to a priori target spectra.
In recent decades, many TD algorithms have been developed. The probability and statistics-based model is one of the classic models for TD, and some traditional TD algorithms in the literature are proposed based on it. For example, constrained energy minimization (CEM) [
8] is a finite impulse response filter which minimizes the total output energy, subjected to the constraint that the output of the target’s spectral signature is one. In general, CEM can well distinguish the targets from background. The adaptive cosine/coherence estimation (ACE) [
9] and the adaptive matching filter (AMF) [
10] algorithms are based on Kelly’s generalized likelihood ratio detection operator. An assumption is made in these two detectors that a mixed background model is subject to the multivariate Gaussian distribution, and additive noise has already been included in the background. Generally, ACE and AMF detectors can be deduced from the generalized likelihood ratio test (GLRT). Closed-form solutions and a considerably low computational cost for TD can be usually achieved with these algorithms. Equipped with the characters above, CEM, ACE and AMF are commonly cited as classic TD methods in many works related to TD [
11,
12,
13].
In fact, a vast majority of HSI pixels often belong to the background, while the target only presents in a few of pixels. The aforementioned methods commonly suffer from background interference, which may give rise to false alarms. Therefore, it is very important to improve the performance of the target detection algorithm by minimizing the interference of background pixels.
Several state-of-the-art hyperspectral target detection techniques have been proposed to cope with this issue. Two types of methods, i.e., background suppression and background separation, are usually considered. The hierarchical constrained energy minimization (hCEM) method [
14] consists of different layers of classic CEM detectors, and the output of each layer is considered as the coefficient of the next layer to suppress the background. The geometric matched filter (GMF) method [
15,
16] provides an inventive way to suppress the background by estimating the background endmembers. Both hCEM and GMF detectors are adept at suppression of the undesired background, which enables them to better concentrate on the targets and acquire remarkable detection results.
Apart from background suppression, background separation also plays an important role in minimizing the background interference and has drawn increasing interest. It is thus very meaningful to develop a preprocessing method aiming at separating the background and target apart before TD is performed. Principal component analysis (PCA) [
17] can be used to address this problem and decompose the HSI into a principal component (PC) part and a residual part. The background is supposed to mainly exist in the PC part, while the target is regarded to present in the residual part. Another famous method, independent component analysis (ICA), which is an unsupervised blind source separation technique has also already been used in hyperspectral data processing [
18] and background separation [
19]. ICA can separate endmembers apart, which enables its application in hyperspectral target detection. Nevertheless, ICA only shows independent sources instead of signal intensity, which makes it difficult to discriminate the part containing the target. In this sense, the connection between ICA and our feature extraction preprocessing scheme actually seems to be considerably weak. Therefore, we consider PCA as a feasible approach for background separation.
However, PCA only takes advantage of the one-dimensional spectral features of samples, which signifies the absence of any spatial dimensions. As a result, PCA fails to employ the spatial information of HSI, which has become an inevitable drawback in HSI processing. Unfortunately, the spatial information of HSI is always of great significance. This fact prevents PCA from being a desirable instrument for HSI processing.
Mathematically, a hyperspectral data cube is actually a three-dimensional tensor [
20,
21]. A tensor can be considered as an extension of a vector and matrix, which provides a convenient way to capture inter-correlation among multiple dimensions. In this sense, both the spectral and spatial domains of hyperspectral datasets can always be better integrated jointly in such an expression than can the situations in two-dimensional matrices. Based on a tensor model, Tucker decomposition [
22,
23] can take spatial-spectral information into consideration, but the essential difference between the spatial and spectral domains of HSI is ignored. These two classic feature extraction methods perform poorly in reflecting the HSI’s intrinsic physical properties.
Despite this, algorithms related to tensor have attracted increasing attention in the field of hyperspectral image processing [
24,
25,
26] recently. Tensor principal component analysis (TPCA) has been successfully used for feature extraction in HSI classification [
27]. Furthermore, we consider that the TPCA model compared with the PCA and Tucker decomposition models should have much more potential in TD.
The TPCA method is a tensorial counterpart of PCA. Not only can the spatial and spectral information of HSI can be considered jointly, but they are also dealt with discriminatively. In this sense, TPCA can overcome the drawbacks of PCA and Tucker decomposition, leading to a more accurate description of the HSI physical property. Inspired by this, we propose a TPCA-based TD preprocessing method. This method first operates TPCA on HSI, thus the imagery is decomposed into the sum of the principal component imagery and the residual imagery. Since the target only takes up a very small number of pixels, it is often considered to appear in the residual imagery. In this way, the target detection in the residual imagery of TPCA can largely reduce the interference of the background, which belongs to the principal component of the HSI, and produce better detection results. Experimental results demonstrate that by using the proposed algorithm, the background and the target can be well separated, and better TD results can be obtained than those with other feature extraction methods and state-of-the-art methods aiming at background suppression.
The main contributions of our work can be concisely summarized as follows:
(1) A tensor principal component analysis (TPCA) based preprocessing method is proposed for hyperspectral target detection. Due to the decomposition of HSI into a principal part and a residual part for detection, the effect of background interference on detection can be well alleviated.
(2) The spectral and spatial information of HSI is not only employed jointly but is also treated differently in the proposed preprocessing method. This characteristic is especially important for revealing the HSI physical property in detection, which also makes the proposed method have superior performance compared with that of PCA and Tucker decomposition.
(3) TPCA-based detectors that combine the proposed preprocessing method with traditional detectors can significantly improve the performance of TD and can outperform some state-of-the-art TD methods.
Traditional hyperspectral target detectors suffer from background interference; thus, their detection performances are limited. Feature extraction preprocessing methods used for background separation, such as PCA and Tucker decomposition, fail to reveal the physical meaning of the HSI’s spatial or spectral domains. The precise objective of this work is to provide a TPCA-based preprocessing method for TD through which both spatial and spectral information of HSI can be reasonably exploited.
The remainder of this paper is organized as follows.
Section 2 reviews the related works. Our TPCA-based target detection algorithm is presented in
Section 3.
Section 4 shows the experimental results. The discussion is shown in
Section 5. Conclusions are drawn in
Section 6.