Next Article in Journal
Discrimination of Leaves in a Multi-Layered Mediterranean Forest through Machine Learning Algorithms
Next Article in Special Issue
Double-Factor Tensor Cascaded-Rank Decomposition for Hyperspectral Image Denoising
Previous Article in Journal
Unraveling Regional Patterns of Sea Level Acceleration over the China Seas
Previous Article in Special Issue
A Collaborative Superpixelwise Autoencoder for Unsupervised Dimension Reduction in Hyperspectral Images
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Method Based on GPU for Real-Time Anomaly Detection in Airborne Push-Broom Hyperspectral Sensors

1
Key Laboratory of Space Active Opto-Electronics Technology, Shanghai Institute of Technical Physics, Chinese Academy of Sciences, Shanghai 200083, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
Hangzhou Institute for Advanced Study, University of Chinese Academy of Sciences, Hangzhou 310024, China
4
Research Center for Intelligent Sensing Systems, Zhejiang Laboratory, Hangzhou 311100, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Remote Sens. 2023, 15(18), 4449; https://doi.org/10.3390/rs15184449
Submission received: 5 July 2023 / Revised: 1 September 2023 / Accepted: 5 September 2023 / Published: 10 September 2023
(This article belongs to the Special Issue Hyperspectral Remote Sensing Imaging and Processing)

Abstract

:
The airborne hyperspectral remote sensing systems (AHRSSs) acquire images with high spectral resolution, high spatial resolution, and high temporal dimension. While the AHRSS captures more detailed information from the terrain objects, the computational complexity of data processing is greatly increased. As an important application technology in the hyperspectral domain, anomaly detection (AD) processing must be real-time and high-precision in many cases, such as post-disaster rescue, military battlefield search, and natural disaster detection. In this paper, the real-time AD technology for the push-broom AHRSS is studied, the mathematical model is established, and a novel implementation framework is proposed. Firstly, the optimized kernel minimum noise fraction (OP-KMNF) transformation is employed to extract informative and discriminative features between the background and anomalies. Secondly, the Nyström method is introduced to reduce the computational complexity of OP-KMNF transformation by decomposing and extrapolating the sub-kernel matrix to estimate the eigenvector of the entire kernel matrix. Thirdly, the extracted features are transferred to hard disks for data storage. Then, taking the extracted features as input data, the background separation model-based CEM anomaly detector (BSM-CEMAD) is imported to detect anomalies. Finally, graphics processing unit (GPU) parallel computing is utilized in the Nyström-based OP-KMNF (NOP-KMNF) transformation and the BSM-CEMAD to improve the execution efficiency, and the real-time AD for the push-broom AHRSS could be realized. To test the feasibility of the implementation framework proposed in this paper, the experiment is carried out with the Airborne Multi-Modular Imaging Spectrometer (AMMIS) developed by the Shanghai Institute of Technical Physics as the data acquisition platform. The experimental results show that the proposed method outperforms many other state-of-the-art AD methods in anomalies detection and background suppression. Moreover, under the condition that the downlink data could retain most of the hyperspectral data information, the proposed method achieves real-time detection of pixel-level anomalies, with the initial delay not exceeding 1 s, the false alarm rate (FAR) less than 5%, and the true positive rate (TPR) close to 98%.

1. Introduction

With the development of remote sensing technology, the acquired images have undergone a development process from panchromatic to multispectral to hyperspectral [1]. Compared with panchromatic and multispectral images, hyperspectral images provide richer spectral information, which helps to identify the subtle features of the terrain objects [2]. Recently, hyperspectral remote sensing has been widely used in various fields, such as military reconnaissance, urban planning, disaster assessment, resource exploration, and environmental monitoring [3,4,5,6,7]. According to the platforms hyperspectral sensors carried, they can be divided into ground-based hyperspectral remote sensing systems (GHRSSs), spaceborne hyperspectral remote sensing systems (SHRSSs), and airborne hyperspectral remote sensing systems (AHRSSs). In practice, the payloads of GHRSSs could be replaced according to the requirement for data acquisition, while the GHRSSs is challenging to be widely used because of the limited data acquisition areas when they are worked. The SHRSSs could obtain data in a large area, while their payloads are tough to replace and maintain. Additionally, the area of data acquisition is limited by the satellite transit situation, and it is almost impossible to collect data in a timely manner in areas that face emergencies. The AHRSS makes up for the shortcomings of the GHRSS and the SHRSS. The payloads of AHRSSs could be replaced according to the demand in practice, allowing for capturing data rapidly in areas where the human resources are difficult to reach or involve. Moreover, the flight altitude setting of the airborne platforms is flexible, and multi-scale remote sensing data could be collected according to the requirements [8]. The basic information of the representative AHRSS is shown in Table 1, and it can be seen that push-broom scanning represents the mainstream of the AHRSS.
While the AHRSS captures more detailed information from the terrain objects, the computational complexity of data processing is greatly increased [6]. As an important application technology in the hyperspectral domain, anomaly detection (AD) processing must be real-time and high-precision in many cases, such as post-disaster rescue, military battlefield search, and natural disaster detection [9]. Methods based on deep learning, by mining the deep features of images, have emerged as a research hotspot in recent years for achieving anomaly target detection. Depending on whether annotated samples are needed, deep learning-based anomaly detection methods can be categorized into supervised and unsupervised anomaly detection. Current anomaly detection methods based on deep learning are depicted in Table 2. Over the past few years, researchers have developed numerous methods for hyperspectral anomaly detection. However, several critical technical challenges still persist in practical applications: (1) parameter tuning, (2) algorithm generality, and (3) real-time capability [10]. Furthermore, with the advancement of automated applications, conducting real-time anomaly detection and transmitting results to the ground for intelligent decision making has also become a pressing technological challenge that needs to be overcome. As mentioned above, push-broom scanning represents the mainstream of the AHRSS. Many related methods on the real-time AD for the push-broom AHRSS are proposed, mainly focusing on fast matrix computation, hardware design, and parallel computation [10]. For matrix inverse adaptive calculation and optimization, some algorithms, such as the real-time recursive causality [11], the linear algebra libraries and multicore platforms [12], and the Cholesky decomposition and linear algebra [11], are proposed to improve the computational efficiency of anomaly detection methods effectively. In terms of hardware design, a two-module field-programmable gate array (FPGA) [13] based on the Reed-Xiaoli (RX) algorithm and the constrained energy minimization (CEM) algorithm is designed to transform those two algorithms according to requirements. In parallel computation, the parallel computing version of the RX detector is studied on the multicore platform [12] and graphics processing unit (GPU) [14], which achieve a remarkable acceleration effect in execution time.
Although these studies have made some achievements, the real-time anomaly detection for the push-broom AHRSS still faces the following difficulties: (1) the nonlinear characteristics and high dimensionality of hyperspectral images are significant, and it is difficult to extract features effectively; (2) constructing a precise representation of anomalies and background information for complex scenes is challenging, and as a result, the target detection accuracy for AHRSS is limited recently; (3) the hyperspectral image has many bands, leading to high data processing calculation and low processing efficiency; (4) the AHRSS captures more detailed information from the terrain objects, and the images have high research value. Ensuring that the downlink data retain most of the hyperspectral data information while performing real-time processing is also a problem to be considered.
In this paper, the real-time AD technology for the push-broom AHRSS is studied, the mathematical model is established, and a novel implementation framework is proposed. To test the feasibility of the implementation framework, this paper designs an experiment carried out with the Airborne Multi-Modular Imaging Spectrometer (AMMIS) developed by the Shanghai Institute of Technical Physics as the data acquisition platform. The experimental results show that the proposed method outperforms many other AD methods (including RX detector [12], low-probability detector (LPD) [28], abundance- and dictionary-based low-rank decomposition (ADLR) [29], collaborative representation detector with background purification and saliency weight (CRDBPSW) [30], anomaly detection via integration of feature extraction and background purification (FEBPAD) [31], and kernel isolation forest-based hyperspectral anomaly detection method (KIFD) [32]) in anomalies detection and background suppression. Moreover, under the condition that the downlink data could retain most of the hyperspectral data information, the proposed method achieves real-time detection of pixel-level anomalies, with the initial delay not exceeding 1 s, the false alarm rate (FAR) less than 5%, and the true positive rate (TPR) close to 98%.
The main contributions of this paper are introduced as follows:
(1)
Considering the data acquisition, data transmission/storage, and data processing of the push-broom AHRSS, the real-time AD technology for the push-broom AHRSS is studied in this paper, the mathematical model of it is established, and the critical problems that need to be solved can be explained by mathematical theoretical analysis, which provides a theoretical reference for the difficulties in the push-broom AHRSS real-time AD technology.
(2)
A novel AHRSS real-time AD implementation framework is proposed based on the mathematical model. Firstly, the optimized kernel minimum noise fraction (OP-KMNF) transformation is employed to extract informative and discriminative features between the background and anomalies. Secondly, the Nyström method is introduced to reduce the computational complexity of OP-KMNF transformation by decomposing and extrapolating the sub-kernel matrix to estimate the eigenvector of the entire kernel matrix. Thirdly, the extracted features are transferred to hard disks for data storage. Then, taking the extracted features as input data, the background separation model-based CEM anomaly detector (BSM-CEMAD) is imported to detect anomalies. Finally, graphics processing unit (GPU) parallel computing is utilized in the Nyström-based OP-KMNF (NOP-KMNF) transformation and the BSM-CEMAD to improve the execution efficiency, and the real-time AD for the push-broom AHRSS could be realized.
(3)
To test the feasibility of the implementation framework, this paper designs an experiment carried out with the Airborne Multi-Modular Imaging Spectrometer (AMMIS) developed by the Shanghai Institute of Technical Physics as the data acquisition platform. The experimental results show that the proposed method outperforms many other AD methods (including RX detector, LPD, ADLR, CRDBPSW, FEBPAD, and KIFD) in anomalies detection and background suppression. Moreover, under the condition that the downlink data could retain most of the hyperspectral data information, the proposed method achieves real-time detection of pixel-level anomalies, with the initial delay not exceeding 1 s, the false alarm rate (FAR) less than 5%, and the true positive rate (TPR) close to 98%.
The remainder of the manuscript is organized as follows: In Section 2, the proposed methods are described in detail. The experimental results are shown in Section 3. Section 4 analyzes and discusses the results, and conclusions are presented in Section 5.

2. Proposed Methods

As mentioned above, the real-time AD for the push-broom AHRSS still faces many difficulties. Considering the data acquisition, data transmission/storage, and data processing of the push-broom AHRSS, the real-time AD technology based on the graphics processing unit (GPU) for the push-broom AHRSS is studied in this paper. In this section, the mathematical model of the real-time AD technology for the push-broom AHRSS is established, and the GPU-based parallel computing is described.

2.1. Mathematical Model Analysis

The realization of the push-broom AHRSS real-time AD could be abstracted as the mapping function:
y = f y x ; η ,
where y represents the expected output of the push-broom AHRSS real-time AD, x represents the image data captured by the push-broom AHRSS, η represents the parameter to be determined in the mapping function, and f y represents the abstracted mapping function dependent on the parameter η .
To solve the difficulties in the real-time AD for the push-broom AHRSS, a novel AHRSS real-time AD implementation framework is proposed in this paper. The implementation flowchart is shown in Figure 1. It can be seen that the factors affecting the expected output y include the effectiveness of the features extracted by feature extraction, the efficiency and reliability of data transmission/storage, the detection capability of the anomaly detection method, and the execution efficiency of the data processing. The determination of the parameter η in the mapping function is affected by these factors, which can be abstracted as the mapping function:
η = ε d t f η ( d r ;   a d ) ,
where d r , a d , ε , and d t represent the effects of the effectiveness of the features extracted by feature extraction, the detection capability of the anomaly detection method, the execution efficiency of the data processing and transmission/storage, and the reliability of data transmission/storage in the determination of the parameter η , respectively.
In the mathematical model, the influence of the execution efficiency of the data processing and transmission/storage and the reliability of data transmission/storage in the determination of the parameter η is represented by the perturbation factors ε and d t in the mapping function, which plays a decisive role in the whole mapping relationship. The execution efficiency of the data processing and transmission/storage, ε , is affected by two factors—(a) the execution efficiency of the data processing ε d r a d and (b) the execution efficiency of data transmission/storage, ε d t — and the mathematical relationship could be expressed as:
ε = ε d r a d + ε d t
After the mathematical model of the real-time AD technology based on the GPU for the push-broom AHRSS is established, the critical problems that need to be solved can be explained by mathematical theoretical analysis. Improving the reliability of data transmission/storage and the execution efficiency of the data processing and transmission/storage can be described as the process of optimizing the perturbation factors ε and d t , while improving the feature extraction performance and AD accuracy can be described as the process of optimizing the parameters d r   a n d   a d in the mapping function.
In previous works, feature extraction, anomaly detection, and algorithm complexity optimization are studied [33,34,35]. The optimized kernel minimum noise fraction (OP-KMNF) transformation for feature extraction [33], the background separation model-based CEM anomaly detector (BSM-CEMAD) for anomaly detection [34], and the Nyström-based OP-KMNF (NOP-KMNF) transformation for algorithm complexity optimization [35] are proposed. On the problem of the inaccurate noise estimation in kernel minimum noise fraction (KMNF) transformation, a mixed noise estimation model (MNEM) combining the Gaussian prior denoising model, the Sobel operator, and the median filter is proposed for OP-KMNF transformation. The MNEM is more robust and effective, retains more edge features and details, and is more suitable for noise estimation during KMNF transformation. Experiments using various hyperspectral datasets with different spatial and spectral resolutions were conducted. The results show that the OP-KMNF transformation overperforms other feature extraction methods (including linear discriminant analysis (LDA) [36], principal component analysis (PCA) [37], minimum noise fraction (MNF) transformation [38], optimized MNF (OMNF) [39], factor analysis (FA) [40], kernel PCA (KPCA) [41], KMNF transformation [42], optimized KMNF (OKMNF) transformation [43], and local preserving projections (LPPs) [44]) for feature extraction and has better adaptability to the changes in spatial resolution and spectral resolution [33]. To reduce the computational complexity of the KMNF transformation, a Nyström-based KMNF transformation (NKMNF) feature extraction method is proposed. The entire kernel matrix of the feature vector is estimated by decomposition and extrapolation of the sub-kernel matrix. The experimental results demonstrate that the NKMNF transformation has lower computational complexity and a better feature extraction performance than KMNF [34]. Similarly, the Nyström method is also applicable to the OP-KMNF transformation to reduce the computational complexity and improve the execution efficiency. Aiming to obtain accurate background and abnormal pixel sets, a background separation model (BSM) combining outlier removal, an iterative strategy, and an RX detector is proposed. The BSM-CEMAD takes the background and abnormal pixel sets as the input of CEMAD to improve the anomaly detection capability. The experimental results show that the BSM-CEMAD has better anti-noise performance, anomalies detection, and background suppression ability, and has better adaptability to the changes in spatial resolution and spectral resolution [35].
To overcome the difficulties in the real-time AD for the push-broom AHRSS, a novel AHRSS real-time AD implementation framework is proposed based on the mathematical model. The technical realization flowchart is shown in Figure 2. Firstly, the OP-KMNF transformation is employed to extract informative and discriminative features between the background and anomalies. Secondly, the Nyström method is introduced to reduce the computational complexity of OP-KMNF transformation by estimating the eigenvector of the entire kernel matrix by decomposing and extrapolating the sub-kernel matrix. Thirdly, the extracted features are transferred to hard disks for data storage. Then, taking the extracted features as input data, the BSM-CEMAD is imported to detect anomalies. Finally, graphics processing unit (GPU) parallel computing is utilized in the Nyström-based OP-KMNF (NOP-KMNF) transformation and the BSM-CEMAD to improve the execution efficiency, and the real-time AD for the push-broom AHRSS could be realized.

2.2. GPU-Based Parallel Computing

According to the type of hardware platforms, hyperspectral image real-time processing technology can be divided into three categories: (1) real-time hyperspectral image processing technology based on FPGA; (2) real-time hyperspectral image processing technology based on cloud computing; (3) real-time hyperspectral image processing technology based on GPU. In practice, GPU and FPGA cannot process large-scale hyperspectral images because of the limitation of their memory capacity, while the cloud computing has no such limitations. The low power consumption makes FPGA more suitable for real-time processing of AHRSS without a data downlink, while its acceleration performance is not superior to GPU. Compared to FPGA and cloud computing, GPU has lower cost and better acceleration performance, which is the most cost-effective hardware platform in high-performance computing [45]. The comparison of different hardware platforms is shown in Table 3. Recently, GPU-related technologies have developed rapidly. Taking NVIDIA GPU products as an example, its organizational structure has been upgraded from Fermi architecture to Ampere architecture, and the processing power has been significantly improved. The development of NVIDIA GPU architecture is shown in Figure 3.
Basic linear algebra subroutines (BLASs) provide a series of interface standards for basic linear algebra operations and are widely used to develop various high-quality linear algebra software, such as LAPACK and FLAME [46,47]. To perform vector and matrix operations effectively, NVIDIA computing unified device architecture (CUDA) provides a BLAS library (CUBLAS) based on the GPU parallel computing platform. Some relevant research demonstrates that CUBLAS has superior performance in processing vector and matrix operations [48,49]. An analysis of NOP-KMNF transformation and BSM-CEMAD show that voluminous vector and matrix operations exist in their implementation program. These operations mainly encompass vector multiplication, matrix multiplication, and matrix eigenvalue computation. The overall algorithmic time complexity of the processing workflow could be represented as O ( N 2 ) . In this paper, the NVIDIA CUBLAS is utilized to develop the GPU parallel computing-based NOP-KMNF (GNOP-KMNF) transformation and BSM-CEMAD (GBSM-CEMAD). The computational efficiency of GNOP-KMNF transformation and GBSM-CEMAD are evaluated in the hardware composed of an NVIDIA GeForce RTX2060 GPU card and an Intel (R) Core (TM) i7-10750H CPU. The NVIDIA GeForce RTX2060 GPU hardware specifications are shown in Table 4.

3. Results

In this section, the experimental results of the method proposed in this paper are described. In Section 3.1, the data acquisition platform AMMIS VNIR hyperspectral push-broom sensor is introduced. In Section 3.2, the real hyperspectral dataset named Vehicle dataset captured by AMMIS VNIR hyperspectral push-broom sensor and conducted in the experiments is presented. In Section 3.3, the evaluation criteria of the proposed methods are described. Two experiments are designed to assess the feasibility of the implementation framework in real-time anomaly detection for the push-broom AHRSS. In order to evaluate the effectiveness of the real-time processing, the Vehicle dataset is used as the captured image, and the runtimes of (1) feature extraction, (2) data transmission/storage, and (3) anomaly detection with different data sizes along the flight direction are tested. To ensure the reliability of the experimental results, each experiment is conducted five times, and the average values are used to analyze the results, which are shown in Section 3.4. The second experiment is designed to assess the AD capability of the proposed method, and the results in the Vehicle dataset are shown in Section 3.5.

3.1. Data Acquisition Platform

The AMMIS is an AHRSS developed by the Shanghai Institute of Technical Physics, Chinese Academy of Sciences. It is designed for various applications, such as land management, urban planning, natural disaster monitoring and assessment, agricultural and forestry surveys, water resources utilization, water quality monitoring, and so on. It can be configured in different modules according to the application requirements [50]. In this paper, the AMMIS visible and near-infrared range (VNIR) hyperspectral push-broom sensor is utilized as the data acquisition platform. The AMMIS VNIR hyperspectral push-broom sensor specifications are shown in Table 5, and its imaging diagram is shown in Figure 4.

3.2. Input Data

The Vehicle dataset used in this paper is collected by the AMMIS VNIR hyperspectral push-broom sensor in Xiong’an New Area, Hebei Province, China. The dataset consists of 2048 pixels along the across-track direction and 15,000 pixels along the flight direction. The wavelength range is 0.40~0.95 μ m , and the spatial resolution is 0.5 m. The pseudo-color and the reference map are shown in Figure 5.

3.3. Evaluation Criteria

To assess the performance of the anomaly detection capability of the proposed method, the 3D receiver operating characteristic (3D ROC) curve, the 2D ROC curve of ( P D , P F ), the 2D ROC curve of ( P D , t ), and the 2D ROC curve of ( P F , t ) are visualized in this paper. Furthermore, the area under the 2D ROC curve of ( P D , P F ) ( A U C   ( D ,   F ) ), the area under the 2D ROC curve of ( P D , t ) ( A U C   ( D ,   t ) ), the area under the 2D ROC curve of ( P F , t ) ( A U C   ( F ,   t ) ), the AUC value of target detectability ( A U C T D ), the AUC value of background suppressibility ( A U C B S ), the AUC value of target detection in the background ( A U C T D B S ), the AUC value of the overall detection probability ( A U C O D P ), the AUC value of overall detection ( A U C O D ), and the AUC value of the signal-to-noise probability ratio ( A U C S N P R ) are used as evaluation indexes to quantitatively evaluate the anomaly detection capability [51]. The calculation formulas of these evaluation criteria are shown in Table 6.
In practical applications, the specific false alarm rate (FAR) and true positive rate (TPR) give the operator a more intuitive understanding of the performance of the methods. As a matter of common knowledge, the lower the false alarm rate, the higher the detection rate, the better the method’s performance is. In this paper, a novel evaluation criterion based on the FAR and TPR in the 3D ROC curve is proposed in this section.
I n d e x = T P R 1 F A R .
As shown in Formula (4), the I n d e x values are calculated based on the FAR and TPR corresponding to each coordinate point on the 3D ROC curve. The larger the I n d e x value, the better the detection results are. Taking the F A R m i n and the T P R m a x at the maximum Index value together as the evaluation criterion of the anomaly detection methods, the results can be assessed more intuitively.

3.4. Real-Time Processing Analysis

Experiments on hyperspectral images with different data sizes are conducted in this section. To ensure the reliability of the experimental results, each experiment is conducted five times, and the average values are used to analyze the results. Figure 6 shows the computational cost comparisons for NOP-KMNF-Ratio transformation and GNOP-KMNF-Ratio transformation with different data sizes, and the computational cost comparisons for BSM-CEMAD and GBSM-CEMAD with different data sizes are shown in Figure 7.
The results shown in Figure 6 and Figure 7 indicate that (1) with the increase in data sizes, the acceleration effect of GNOP-KMNF transformation rather than that of NOP-KMNF-Ratio transformation and GBSM-CEMAD rather than that of BSM-CEMAD are more and more significant; (2) when the data size is 288 × 288 × 250, the execution efficiency of GNOP-KMNF transformation is about 336 times that of NOP-KMNF transformation, which shows the superior performance of the GNOP-KMNF transformation feature extraction algorithm in execution efficiency; (3) when the data size is 128 × 128 × 250, the execution efficiency of GBSM-CEMAD is about 223 times that of BSM-CEMAD, which indicates the superior performance of GBSM-CEMAD in execution efficiency.
As mentioned above, the data processing includes three parts: (1) feature extraction, (2) data transmission/storage, and (3) anomaly detection. The Vehicle dataset is used as the captured image, and the runtimes of these three data processing parts with different data sizes along the flight direction are tested. To ensure the reliability of the experimental results, each experiment is conducted five times, and the average values are used to analyze the results. The execution efficiency of GNOP-KMNF transformation, data transmission/storage, and GBSM-CEMAD with different data sizes are shown in Figure 8, Figure 9, and Figure 10, respectively. To quantitatively show the experimental results, the runtimes of GNOP-KMNF transformation, data transmission/storage, and GBSM-CEMAD with different data sizes are shown in Table 7.
Analyzing the data in Table 6, it can be seen that when the data size is 2048 × 250 × 160, the whole processing time of GNOP-KMNF transformation, data transmission/storage, and GBSM-CEMAD is 991.61 ms. The frame frequency of the AMMIS VNIR hyperspectral sensor is 160 Hz, the data acquisition rate is 313 MB/s, and the processing speed of the proposed method is 315.65 MB/s. In real-time processing, the hyperspectral image data of 2048 × 250 × 160 are used as a hyperspectral image cube to be processed. With the initial delay not exceeding 1 s, the real-time AD for the AMMIS VNIR hyperspectral sensor can be realized. The real-time AD processing for the AMMIS VNIR hyperspectral sensor is shown in Figure 11, and the anomaly detection results of the proposed method for the Vehicle dataset is shown in Figure 12.

3.5. Anomaly Detection Capability Evaluation

In this section, the anomaly detection capability of the proposed method and other state-of-the-art algorithms (including RX detector, LPD, ADLR, CRDBPSW, FEBPAD, and KIFD) is evaluated. Because the size of the Vehicle dataset is too large, part of the anomaly detection results using different detectors is shown in Figure 13 to display the detection results better.
Part of the anomaly detection results using different detectors for the Vehicle dataset are shown in Figure 13. By observing these result images, it can be seen that the proposed method suppresses more of the background and highlights the anomalies from the background. Compared with RX, ADLR, CRDBPSW, and FEBPAD, LPD and KIFD detect more anomaly pixels and obtain an acceptable background suppression.
To quantitatively assess the results, the 3D ROC curve and its generated three 2D ROC curves of the Vehicle dataset are demonstrated in Figure 14, and the AUC values calculated from three 2D ROC curves for the Vehicle dataset are shown in Table 8.
As shown in Figure 14, the proposed method shows better performance in the ROC curves. For the 2D ROC curves of ( P F , t ), the curve of the proposed method is closest to the bottom-right corner. Moreover, the 2D ROC curves of ( P D , P F ) and ( P D , t ) of the proposed method are closest to the top-left corner. The results demonstrate the significant performance of the proposed method.
For the AUC values calculated from three 2D ROC curves for the Vehicle dataset shown in Table 8, the A U C ( D , t ) of ADLR is better than that of the proposed method, while the A U C F , t and other evaluation criteria values of ADLR are much worse than that of the proposed method. Additionally, the A U C F , t of FEBPAD outperforms that of the proposed method, while the A U C ( D , t ) and other evaluation criteria values of FEBPAD are much worse than that of the proposed method. In summary, the results show that the proposed method has a better anomalies detection and background suppression performance.
The Index values based on the FAR and TPR in the 3D ROC curves for different detectors are calculated, which are shown in Table 9. When the I n d e x value is the largest, the F A R m i n of the proposed method is 0.045527, and the T P R m a x of the proposed method is 0.979167. The results demonstrate that the proposed method can achieve real-time detection of pixel-level anomalies, with the initial delay not exceeding 1 s, the FAR of less than 5%, and the TPR close to 98%.

4. Discussion

In this paper, two experiments are designed to assess the effectiveness of the implementation framework in real-time anomaly detection for the push-broom AHRSS.
In Section 3.4, the computational cost comparisons for NOP-KMNF-Ratio transformation and GNOP-KMNF-Ratio transformation with different data sizes are tested. The results indicate that (1) with the increase in data sizes, the acceleration effect of GNOP-KMNF transformation rather than that of NOP-KMNF-Ratio transformation and GBSM-CEMAD rather than that of BSM-CEMAD are more and more significant; (2) when the data size is 288 × 288 × 250, the execution efficiency of GNOP-KMNF transformation is about 336 times that of NOP-KMNF transformation, which shows the superior performance of the GNOP-KMNF transformation feature extraction algorithm in execution efficiency; (3) when the data size is 128 × 128 × 250, the execution efficiency of GBSM-CEMAD is about 223 times that of BSM-CEMAD, which indicates the superior performance of GBSM-CEMAD in execution efficiency. In order to evaluate the effectiveness of the real-time processing, the Vehicle dataset is used as the captured image, and the runtimes of three data processing parts with different data sizes along the flight direction are tested. It can be seen that when the data size is 2048 × 250 × 160, the whole processing time of GNOP-KMNF transformation, data transmission/storage, and GBSM-CEMAD is 991.61 ms. The frame frequency of the AMMIS VNIR hyperspectral sensor is 160 Hz, the data acquisition rate is 313 MB/s, and the processing speed of the proposed method is 315.65 MB/s. In real-time processing, the hyperspectral image data of 2048 × 250 × 160 are used as a hyperspectral image cube to be processed. With the initial delay not exceeding 1 s, the real-time AD for the AMMIS VNIR hyperspectral sensor could be realized. In Section 3.5, the experiment is designed to assess the anomaly detection capability of the proposed methods. The experimental results show that the proposed method outperforms many other AD methods (including RX detector, LPD, ADLR, CRDBPSW, FEBPAD, and KIFD) in anomalies detection and background suppression. Moreover, under the condition that the downlink data could retain most of the hyperspectral data information, the proposed method achieves real-time detection of pixel-level anomalies, with the initial delay not exceeding 1 s, the FAR less than 5%, and the TPR close to 98%.
In this paper, the real-time AD technology for the push-broom AHRSS is studied, and a novel implementation framework is proposed. In order to achieve high-precision real-time anomaly detection for the push-broom AHRSS, the whole implementation process of hyperspectral image acquisition, transmission, and processing is analyzed, and the mathematical model of the real-time AD technology for the push-broom AHRSS is established. On this basis, the solutions to the current difficulties of this research are given from the mathematical point of view, which lays a theoretical foundation for the subsequent implementation. The processing power of GPUs is measured by the ability to perform floating-point operations per second (FLOPS). In practical applications, the processing power of GPUs depends on many factors, such as hardware architecture, number of cores, frequency, and memory bandwidth. Figure 15 shows the thermal design power of these NVIDIA desktop GPUs, and the processing power of NVIDIA desktop GPUs since 2015 are shown in Figure 16.
Compared with the NVIDIA GeForce RTX 2060 used in the experiments, the processing power of the most powerful NVIDIA GeForce RTX 4090 on single-precision floating-point operations has increased by 12.80 times, and the processing power on double-precision floating-point operations has increased by 6.40 times. With the continuous development of GPU parallel computing technologies, the execution efficiency of the algorithm will be further improved, and the initial delay time will be effectively shortened. While for the thermal design power, the most powerful NVIDIA GeForce RTX 4090 has a 2.81-fold increase compared with the NVIDIA GeForce RTX 2060 used in the experiments, how to select the hardware platform based on the careful consideration of execution efficiency and power consumption in actual flight missions still needs further discussion.
Moreover, the implementation framework designed in this paper on the real-time AD technology based on the GPU for the push-broom AHRSS is guided by practical application. In practice, the data acquisition speed by the push-broom AHRSS is affected not only by the frame frequency of the sensors, but also by the flight speed of the airborne platform. Carrying out simulation experiments only considering the push-broom AHRSS specifications is limited. It is necessary to design and conduct experiments on real-time anomaly detection for the push-broom AHRSS in combination with actual flight missions.

5. Conclusions

In this paper, the real-time AD technology for the push-broom AHRSS is studied, the mathematical model is established, and a novel implementation framework is proposed. To test the feasibility of the implementation framework proposed in this paper, the experiment is carried out with the AMMIS developed by the Shanghai Institute of Technical Physics as the data acquisition platform. The experimental results show that the proposed method outperforms many other AD methods (including RX detector, LPD, ADLR, CRDBPSW, FEBPAD, and KIFD) in anomalies detection and background suppression. Moreover, under the condition that the downlink data could retain most of the hyperspectral data information, the proposed method achieves real-time detection of pixel-level anomalies, with the initial delay not exceeding 1 s, the FAR less than 5%, and the TPR close to 98%. Through the analysis of the critical issues involved in this research, a relatively complete implementation framework with an algorithm design and data processing for the real-time AD technology on the push-broom AHRSS is proposed in this paper. Furthermore, the mathematical model of the real-time AD technology for the push-broom AHRSS is constructed, the algorithm design and verification are carried out, and the test experiment validation is conducted. In this paper, the AMMIS developed by Shanghai Institute of Technical Physics, Chinese Academy of Sciences, is utilized as the data acquisition sensor to carry out test experiments. While the experiments in this paper are conducted on the obtained data, we are interested in designing and conducting the experiment under the real airborne imaging environment in the future.

Author Contributions

Conceptualization, T.X. and Y.W.; methodology, T.X.; software, T.X. and H.X.; validation, C.W. and H.X.; formal analysis, C.W. and H.X.; investigation, H.X., C.W. and T.X.; resources, Y.W.; data curation, C.W.; writing—original draft preparation, T.X. and H.X.; writing—review and editing, T.X., H.X. and C.W.; visualization, H.X.; supervision, Y.W.; project administration, Y.W.; funding acquisition, Y.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Civil Aerospace Project (No. D040102).

Data Availability Statement

Not applicable.

Acknowledgments

We acknowledge the data support of the Key Laboratory of Space Active Opto-Electronics Technology, Shanghai Institute of Technical Physics, Chinese Academy of Sciences, and the State Key Laboratory of Remote Sensing Science, Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Shaw, G.; Manolakis, D. Signal processing for hyperspectral image exploitation. IEEE Signal Process Mag. 2002, 19, 12–16. [Google Scholar] [CrossRef]
  2. Schowengerdt, R.A. Remote Sensing: Models and Methods for Image Processing; Academic Press: Cambridge, MA, USA, 1997. [Google Scholar]
  3. Kennedy, R.E.; Townsend, P.A.; Gross, J.E.; Cohen, W.B.; Bolstad, P.; Wang, Y.Q.; Adams, P. Remote sensing change detection tools for natural resource managers: Understanding concepts and tradeoffs in the design of landscape monitoring project. Remote Sens. Environ. 2009, 113, 1382–1396. [Google Scholar] [CrossRef]
  4. Kerle, N.; Oppenheimer, C. Satellite remote sensing as a tool in lahar disaster management. Disasters 2002, 26, 140–160. [Google Scholar] [CrossRef] [PubMed]
  5. Jensen, J.R.; Cowen, D.C. Remote Sensing of Urban/Suburban Infrastructure and Socio-Economic Attributes. Photogramm. Eng. Remote Sens. 1999, 65, 611–622. [Google Scholar]
  6. Lin, W. Real-Time Target and Anormy Detection Algorithm for Hyperspectral Imagery. Master’s Thesis, Hangzhou Dianzi University, Hangzhou, China, 2017. [Google Scholar]
  7. Ma, D. Hperspectral Anomaly Detection Based on Multiple Feature Low-Rank Representation. Master’s Thesis, China University of Mining and Technology, Beijing, China, 2017. [Google Scholar]
  8. Li, R.; Liu, B.; Fang, C. Aerial remote sensing in the new era. Mod Phys. 2015, 27, 54–60. [Google Scholar] [CrossRef]
  9. Fu, L.; Deng, H.; Liu, C. Novel fast real-time target detection and classification algorithms for hyperspectral imagery. Acta Opt. Sin. 2017, 37, 0230002-1. [Google Scholar] [CrossRef]
  10. Su, H.; Wu, Z.; Zhang, H.; Du, Q. Hyperspectral anomaly detection: A survey. IEEE Geosci. Remote Sens. Mag. 2022, 10, 64–90. [Google Scholar] [CrossRef]
  11. Molero, J.M.; Garzón, E.M.; García, I.; Plaza, A. Analysis and optimizations of global and local versions of the RX algorithm for anomaly detection in hyperspectral data. IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens. 2013, 6, 801–814. [Google Scholar] [CrossRef]
  12. Zhang, L.; Peng, B.; Zhang, F.; Wang, L.; Zhang, H.; Zhang, P.; Tong, Q. Fast real-time causal linewise progressive hyperspectral anomaly detection via Cholesky decomposition. IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens. 2017, 10, 4614–4629. [Google Scholar] [CrossRef]
  13. Yang, B.; Yang, M.; Plaza, A.; Gao, L.; Zhang, B. Dual-Mode FPGA implementation of target and anomaly detection algorithms for real-time hyperspectral imaging. IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens. 2015, 8, 2950–2961. [Google Scholar] [CrossRef]
  14. Paz, A.; Plaza, A. Clusters versus GPUs for parallel target and anomaly detection in hyperspectral images. EURASIP J. Adv. Signal Process. 2010, 2010, 915639. [Google Scholar] [CrossRef]
  15. Bati, E.; Çalışkan, A.; Koz, A.; Alatan, A.A. Hyperspectral anomaly detection method based on auto-encoder. SPIE Remote Sens. 2015, 9643, 220–226. [Google Scholar]
  16. Zhao, C.; Li, X.; Zhu, H. Hyperspectral anomaly detection based on stacked denoising autoencoder. J. Appl. Remote Sens. 2017, 11, 042605. [Google Scholar] [CrossRef]
  17. Zhao, C.; Zhang, L. Spectral-spatial stacked autoencoders based on low-rank and sparse matrix decomposition for hyperspectral anomaly detection. Infrared Phys. Technol. 2018, 92, 166–176. [Google Scholar] [CrossRef]
  18. Zhang, L.; Cheng, B. A stacked autoencoders-based adaptive subspace model for hyperspectral anomaly detection. Infrared Phys. Technol. 2019, 96, 52–60. [Google Scholar] [CrossRef]
  19. Xie, W.; Lei, J.; Liu, B.; Li, Y.; Jia, X. Spectral constraint adversarial autoencoders approach to feature representation in hyperspectral anomaly detection. Neural Netw. 2019, 119, 222–234. [Google Scholar] [CrossRef]
  20. Xie, W.; Liu, B.; Li, Y.; Lei, J.; Chang, C.I.; He, G. Spectral Adversarial Feature Learning for Anomaly Detection in Hyperspectral Imagery. IEEE Trans. Geosci. Remote Sens. 2020, 58, 2352–2365. [Google Scholar] [CrossRef]
  21. Xie, W.; Liu, B.; Li, Y.; Lei, J.; Du, Q. Autoencoder and Adversarial-Learning-Based Semi-supervised Background Estimation for Hyperspectral Anomaly Detection. IEEE Trans. Geosci. Remote Sens. 2020, 58, 5416–5427. [Google Scholar] [CrossRef]
  22. Jiang, T.; Li, Y.; Xie, W.; Du, Q. Discriminative Reconstruction Constrained Generative Adversarial Network for Hyperspectral Anomaly Detection. IEEE Trans. Geosci. Remote Sens. 2020, 58, 4666–4679. [Google Scholar] [CrossRef]
  23. Lu, X.; Zhang, W.; Huang, J. Exploiting Embedding Manifold of Autoencoders for Hyperspectral Anomaly Detection. IEEE Trans. Geosci. Remote Sens. 2020, 58, 1527–1537. [Google Scholar] [CrossRef]
  24. Jiang, K.; Xie, W.; Li, Y.; Lei, J.; He, G.; Du, Q. Semisupervised Spectral Learning with Generative Adversarial Network for Hyperspectral Anomaly Detection. IEEE Trans. Geosci. Remote Sens. 2020, 58, 5224–5236. [Google Scholar] [CrossRef]
  25. Li, W.; Wu, G.; Du, Q. Transferred Deep Learning for Anomaly Detection in Hyperspectral Imagery. IEEE Geosci. Remote Sens. Lett. 2017, 14, 597–601. [Google Scholar] [CrossRef]
  26. Ma, N.; Peng, Y.; Wang, S.; Liu, D. Hyperspectral image anomaly targets detection with online deep learning. In Proceedings of the 2018 IEEE International Instrumentation and Measurement Technology Conference (I2MTC), Houston, TX, USA, 14–17 May 2018. [Google Scholar]
  27. Song, S.; Zhou, H.; Yang, Y.; Song, J. Hyperspectral Anomaly Detection via Convolutional Neural Network and Low Rank with Density-Based Clustering. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2019, 12, 3637–3649. [Google Scholar] [CrossRef]
  28. Harsanyi, J.C. Detection and Classification of Subpixel Spectral Signatures in Hyperspectral Image Sequences. Ph.D. Thesis, Department of Electronic & Electrical Engineering, University of Maryland Baltimore County, Baltimore, MD, USA, 1993. [Google Scholar]
  29. Qu, Y.; Wang, W.; Guo, R.; Ayhan, B.; Kwan, C.; Vance, S.; Qi, H. Hyperspectral Anomaly Detection Through Spectral Unmixing and Dictionary-Based Low-Rank Decomposition. IEEE Trans. Geosci. Remote Sens. 2018, 56, 4391–4405. [Google Scholar] [CrossRef]
  30. Hou, Z.; Li, W.; Tao, R.; Ma, P.; Shi, W. Collaborative Representation with Background Purification and Saliency Weight for Hyperspectral Anomaly Detection. Sci. China Inf. Sci. 2022, 65, 1–12. [Google Scholar] [CrossRef]
  31. Ma, Y.; Fan, G.; Jin, Q.; Huang, J.; Mei, X.; Ma, J. Hyperspectral Anomaly Detection via Integration of Feature Extraction and Background Purification. IEEE Geosci. Remote Sens. Lett. 2021, 18, 1436–1440. [Google Scholar] [CrossRef]
  32. Li, S.; Zhang, K.; Duan, P.; Kang, X. Hyperspectral Anomaly Detection with Kernel Isolation Forest. IEEE Trans. Geosci. Remote Sens. 2020, 58, 319–329. [Google Scholar] [CrossRef]
  33. Xue, T.; Wang, Y.; Chen, Y.; Jia, J.; Wen, M.; Guo, R.; Wu, T.; Deng, X. Mixed Noise Estimation Model for Optimized Kernel Minimum Noise Fraction Transformation in Hyperspectral Image Dimensionality Reduction. Remote Sens. 2021, 13, 2607. [Google Scholar] [CrossRef]
  34. Xue, T.; Jia, J.; Xie, H.; Zhang, C.; Deng, X.; Wang, Y. Kernel Minimum Noise Fraction Transformation-Based Background Separation Model for Hyperspectral Anomaly Detection. Remote Sens. 2022, 14, 5157. [Google Scholar] [CrossRef]
  35. Xue, T.; Wang, Y.; Deng, X. A Novel Method for Fast Kernel Minimum Noise Fraction Transformation in Hyperspectral Image Dimensionality Reduction. Remote Sens. 2022, 14, 1737. [Google Scholar] [CrossRef]
  36. Bandos, T.V.; Bruzzone, L.; Camps-Valls, G. Classification of hyperspectral images with regularized linear discriminant analysis. IEEE Trans. Geosci. Remote Sens. 2009, 47, 862–873. [Google Scholar] [CrossRef]
  37. Roger, R.E. Principal Components transform with simple, automatic noise adjustment. Int. J. Remote Sens. 1996, 17, 2719–2727. [Google Scholar] [CrossRef]
  38. Green, A.; Berman, M.; Switzer, P.; Craig, M. A transformation for ordering multispectral data in terms of image quality with implications for noise removal. IEEE Trans. Geosci. Remote Sens. 1988, 26, 65–74. [Google Scholar] [CrossRef]
  39. Gao, L.; Zhang, B.; Chen, Z.; Lei, L. Study on the Issue of Noise Estimation in Dimension Reduction of Hyperspectral Images. In Proceedings of the 2011 3rd Workshop on Hyperspectral Image and Signal Processing: Evolution in Remote Sensing, Lisbon, Portugal, 6–9 June 2011; pp. 1–4. [Google Scholar]
  40. Mooi, E.; Sarstedt, M. Factor Analysis. In A Concise Guide to Market Research; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar] [CrossRef]
  41. Gu, Y.; Liu, Y.; Zhang, Y. A Selective KPCA Algorithm Based on High-Order Statistics for Anomaly Detection in Hyperspectral Imagery. IEEE Geosci. Remote Sens. Lett. 2008, 5, 43–47. [Google Scholar] [CrossRef]
  42. Bachmann, C.; Ainsworth, T.; Fusina, R. Exploiting manifold geometry in hyperspectral imagery. IEEE Trans. Geosci. Remote Sens. 2005, 43, 441–454. [Google Scholar] [CrossRef]
  43. Gao, L.; Zhao, B.; Jia, X.; Liao, W.; Zhang, B. Optimized Kernel Minimum Noise Fraction Transformation for Hyperspectral Image Classification. Remote Sens. 2017, 9, 548. [Google Scholar] [CrossRef]
  44. Feng, F.; Li, W.; Du, Q.; Zhang, B. Dimensionality Reduction of Hyperspectral Image with Graph-Based Discriminant Analysis Considering Spectral Similarity. Remote Sens. 2017, 9, 323. [Google Scholar] [CrossRef]
  45. Du, Q.; Tang, B.; Xie, W.; Li, W. Parallel and Distributed Computing for Anomaly Detection from Hyperspectral Remote Sensing Imagery. Proc. IEEE 2021, 109, 1306–1319. [Google Scholar] [CrossRef]
  46. Williams, C.; Seeger, M. Using the Nyström method to speed up kernel machines. In Proceedings of the 14th Annual Conference on Neural Information Processing Systems, Denver, CO, USA, 1 January 2001; pp. 682–688. [Google Scholar]
  47. Anderson, E.; Bai, Z.; Bischof, C.; Blackford, L.S.; Demmel, J.; Dongarra, J.; Du Croz, J.; Greenbaum, A.; Hammarling, S.; McKenney, A.; et al. LAPACK Users’ Guide; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 1999. [Google Scholar]
  48. Bientinesi, P.; Gunnels, J.A.; Myers, M.E.; Quintana-Ortí, E.S.; Van De Geijn, R.A. The science of deriving dense linear algebra algorithms. ACM Trans. Math. Softw. 2005, 31, 1–26. [Google Scholar] [CrossRef]
  49. Fujimoto, N. Faster Matrix-Vector Multiplication on GeForce 8800GTX. In Proceedings of the 2011 IEEE International Parallel & Distributed Processing Symposium, Anchorage, AK, USA, 16–20 May 2008; pp. 1–8. [Google Scholar]
  50. Barrachina, S.; Castillo, M.; Igual, F.D.; Mayo, R.; Quintana-Orti, E.S. Evaluation and Tuning of the Level 3 CUBLAS for Graphics Processors. In Proceedings of the 22nd IEEE International Parallel & Distributed Processing Symposium, Miami, FL, USA, 14–18 April 2008; pp. 1–8. [Google Scholar]
  51. Chang, C.-I. An Effective Evaluation Tool for Hyperspectral Target Detection: 3D Receiver Operating Characteristic Curve Analysis. IEEE Trans. Geosci. Remote Sens. 2021, 59, 5131–5153. [Google Scholar] [CrossRef]
Figure 1. Flowchart of the real-time AD for push-broom AHRSS.
Figure 1. Flowchart of the real-time AD for push-broom AHRSS.
Remotesensing 15 04449 g001
Figure 2. Technical realization flowchart of the real-time AD for push-broom AHRSS.
Figure 2. Technical realization flowchart of the real-time AD for push-broom AHRSS.
Remotesensing 15 04449 g002
Figure 3. The development of NVIDIA GPU architecture.
Figure 3. The development of NVIDIA GPU architecture.
Remotesensing 15 04449 g003
Figure 4. The imaging diagram of the AMMIS VNIR hyperspectral push-broom sensor.
Figure 4. The imaging diagram of the AMMIS VNIR hyperspectral push-broom sensor.
Remotesensing 15 04449 g004
Figure 5. The visualization image and reference map of Vehicle dataset.
Figure 5. The visualization image and reference map of Vehicle dataset.
Remotesensing 15 04449 g005
Figure 6. The computational cost comparisons for NOP-KMNF-Ratio transformation and GNOP-KMNF-Ratio transformation with different data sizes.
Figure 6. The computational cost comparisons for NOP-KMNF-Ratio transformation and GNOP-KMNF-Ratio transformation with different data sizes.
Remotesensing 15 04449 g006
Figure 7. The computational cost comparisons for BSM-CEMAD and GBSM-CEMAD with different data sizes.
Figure 7. The computational cost comparisons for BSM-CEMAD and GBSM-CEMAD with different data sizes.
Remotesensing 15 04449 g007
Figure 8. The execution efficiency of GNOP-KMNF transformation with different data sizes.
Figure 8. The execution efficiency of GNOP-KMNF transformation with different data sizes.
Remotesensing 15 04449 g008
Figure 9. The execution efficiency of data transmission/storage with different data sizes.
Figure 9. The execution efficiency of data transmission/storage with different data sizes.
Remotesensing 15 04449 g009
Figure 10. The execution efficiency of GBSM-CEMAD with different data sizes.
Figure 10. The execution efficiency of GBSM-CEMAD with different data sizes.
Remotesensing 15 04449 g010
Figure 11. The real-time AD processing for the AMMIS VNIR hyperspectral sensor.
Figure 11. The real-time AD processing for the AMMIS VNIR hyperspectral sensor.
Remotesensing 15 04449 g011
Figure 12. The anomaly detection result for the Vehicle dataset.
Figure 12. The anomaly detection result for the Vehicle dataset.
Remotesensing 15 04449 g012
Figure 13. Part of the anomaly detection results using different detectors for the Vehicle dataset.
Figure 13. Part of the anomaly detection results using different detectors for the Vehicle dataset.
Remotesensing 15 04449 g013
Figure 14. The 3D ROC curve and its generated three 2D ROC curves of the Vehicle dataset.
Figure 14. The 3D ROC curve and its generated three 2D ROC curves of the Vehicle dataset.
Remotesensing 15 04449 g014
Figure 15. The thermal design power of NVIDIA desktop GPUs since 2015.
Figure 15. The thermal design power of NVIDIA desktop GPUs since 2015.
Remotesensing 15 04449 g015
Figure 16. The processing power of NVIDIA desktop GPUs since 2015.
Figure 16. The processing power of NVIDIA desktop GPUs since 2015.
Remotesensing 15 04449 g016
Table 1. The basic information of the representative AHRSS.
Table 1. The basic information of the representative AHRSS.
The AHRSS NamePlace of DevelopmentProducing YearScanning Mode
AISAmerica1983push-broom
AVIRISAmerica1987whisk-broom
CASI/SASI/MASI/TASICanada1988push-broom
AISAFinland1992push-broom
HyMapAustralia1997whisk-broom
PHIChina1997push-broom
AVISGermany1997push-broom
OMISChina2000whisk-broom
AVIS-2Germany2001push-broom
APEXSwitzerland/Belgium2009push-broom
Hyper-CamCanada2009push-broom
MAKOAmerica2010push-broom
SieletersFrance2013push-broom
AISA-OWLFinland2014push-broom
AMMISChina2016push-broom
HyTESAmerica2016push-broom
AVIRIS-NGAmerica2019push-broom
Table 2. Anomaly detection methods based on deep learning.
Table 2. Anomaly detection methods based on deep learning.
CategoryRepresentative Methods
Unsupervised deep learning-based anomaly detection methodsSparse decomposition and auto-encoder (AE) [15]; stacked AE [16,17,18]; spectral constraint adversarial AE network [19,20,21]; generative adversarial network (GAN) combined with AE [22]; exploiting embedding manifold of autoencoders [23]; and semi-supervised spectral learning with GAN [24].
Supervised deep learning-based anomaly detection methodsConvolutional neural network (CNN) [25]; online CNN [26]; and CNN combined with low-rank representation (LRR) [27].
Table 3. Comparison of different hardware platforms [37].
Table 3. Comparison of different hardware platforms [37].
High Speedup PerformanceLow Power ConsumptionLow CostLarge Memory Capacity
FPGA
Cloud computing
GPU
Table 4. The NVIDIA GeForce RTX2060 GPU hardware specifications.
Table 4. The NVIDIA GeForce RTX2060 GPU hardware specifications.
ParametersPerformance
CUDA Driver Version/Runtime Version11.4/10.0
CUDA Cores1920
Total amount of global memory6144 MBytes
Memory clock rate5501 MHz
Max clock rate1200 MHz
Memory bus width192 bits
Warp size32
Total number of registers available per block65,536
Total amount of shared memory per block49,152 bytes
Total amount of constant memory65,536 bytes
Table 5. AMMIS VNIR hyperspectral push-broom sensor specifications.
Table 5. AMMIS VNIR hyperspectral push-broom sensor specifications.
ParametersPerformance
Spectral range ( μ m )0.40~0.95
Instantaneous field of view ( m r a d )0.25
Signal-to-Noise Ratio 500
Across-track pixels2048
Frame frequency (Hz)160
Effective spectral bands250
Table 6. The calculation formulas of the anomaly detection evaluation criteria.
Table 6. The calculation formulas of the anomaly detection evaluation criteria.
Evaluation CriteriaCalculation Formulas
A U C   ( D ,   F ) the area under the 2D ROC curve of ( P D , P F )
A U C   ( D ,   t ) the area under the 2D ROC curve of ( P D , t )
A U C   ( F ,   t ) the area under the 2D ROC curve of ( P F , t )
A U C T D A U C   D ,   F + A U C   ( D ,   t )
A U C B S A U C   D ,   F A U C   ( F ,   t )
A U C T D B S A U C   D ,   t A U C   ( F ,   t )
A U C O D P A U C   D ,   t + ( 1 A U C   F ,   t )
A U C O D A U C   D ,   F + A U C   ( D ,   t ) A U C   ( F ,   t )
A U C S N P R A U C   ( D ,   t ) / A U C   ( F ,   t )
Table 7. The runtimes of GNOP-KMNF transformation, data transmission/storage, and GBSM-CEMAD with different data sizes.
Table 7. The runtimes of GNOP-KMNF transformation, data transmission/storage, and GBSM-CEMAD with different data sizes.
Data Sizes Runtimes   ( m s )
GNOP-KMNFData Transmission/StorageGBSM-CEMAD
2048 × 250 × 1068.5818.2033.01
2048 × 250 × 2096.4218.4062.43
2048 × 250 × 30127.6818.8091.42
2048 × 250 × 40154.7419.10120.49
2048 × 250 × 50183.4019.20149.86
2048 × 250 × 60215.7819.50179.39
2048 × 250 × 70242.9619.70206.98
2048 × 250 × 80272.2419.80240.27
2048 × 250 × 90302.0620.10262.98
2048 × 250 × 100334.3020.30291.50
2048 × 250 × 110362.3020.40320.75
2048 × 250 × 120392.6220.50352.40
2048 × 250 × 130423.1820.90383.11
2048 × 250 × 140454.6821.10411.90
2048 × 250 × 150482.8621.40441.82
2048 × 250 × 160506.9221.60463.09
Table 8. The AUC values calculated from three 2D ROC curves for the Vehicle dataset.
Table 8. The AUC values calculated from three 2D ROC curves for the Vehicle dataset.
Evaluation CriteriaProposedRXLPDADLRCRDBPSWFEBPADKIFD
A U C ( D , F ) 0.9824760.8662530.8262300.6893210.8353420.8493850.944433
A U C ( D , t ) 0.4996930.0733890.2963740.7597430.1550150.1189040.471660
A U C ( F , t ) 0.0437130.0152640.1242300.6309170.0185770.0122520.062557
A U C T D 1.4821690.9396421.1226041.4490640.9903570.9682891.416093
A U C B S 0.9387630.8509890.7020000.0584040.8167650.8371330.881876
A U C T D B S 0.4559800.0581250.1721440.1288260.1364380.1066520.409103
A U C O D P 1.4559801.0581251.1721441.1288261.1364381.1066521.409103
A U C O D 1.4384560.9243780.9983740.8181470.9717800.9560371.353536
A U C S N P R 11.431224.8079802.3856881.2041898.3444589.7048657.539684
Table 9. The Index values based on the FAR and TPR in the 3D ROC curves for the Vehicle dataset.
Table 9. The Index values based on the FAR and TPR in the 3D ROC curves for the Vehicle dataset.
Methods I n d e x F A R m i n T P R m a x
RX0.6175820.2589010.833333
LPD0.5654610.2729790.777778
ADLR0.7395660.2604341.000000
CRDBPSW0.6259330.2162230.798611
FEBPAD0.5784590.1990560.722222
KIFD0.8097810.0362940.840278
Proposed0.9345880.0455270.979167
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xue, T.; Wang, C.; Xie, H.; Wang, Y. A Novel Method Based on GPU for Real-Time Anomaly Detection in Airborne Push-Broom Hyperspectral Sensors. Remote Sens. 2023, 15, 4449. https://doi.org/10.3390/rs15184449

AMA Style

Xue T, Wang C, Xie H, Wang Y. A Novel Method Based on GPU for Real-Time Anomaly Detection in Airborne Push-Broom Hyperspectral Sensors. Remote Sensing. 2023; 15(18):4449. https://doi.org/10.3390/rs15184449

Chicago/Turabian Style

Xue, Tianru, Chongru Wang, Hui Xie, and Yueming Wang. 2023. "A Novel Method Based on GPU for Real-Time Anomaly Detection in Airborne Push-Broom Hyperspectral Sensors" Remote Sensing 15, no. 18: 4449. https://doi.org/10.3390/rs15184449

APA Style

Xue, T., Wang, C., Xie, H., & Wang, Y. (2023). A Novel Method Based on GPU for Real-Time Anomaly Detection in Airborne Push-Broom Hyperspectral Sensors. Remote Sensing, 15(18), 4449. https://doi.org/10.3390/rs15184449

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop