Next Article in Journal
Analysis of the Recent Agricultural Situation of Dakhla Oasis, Egypt, Using Meteorological and Satellite Data
Previous Article in Journal
Improved SRGAN for Remote Sensing Image Super-Resolution Across Locations and Sensors
Previous Article in Special Issue
Urban Nighttime Leisure Space Mapping with Nighttime Light Images and POI Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New GPU Implementation of Support Vector Machines for Fast Hyperspectral Image Classification

Hyperspectral Computing Laboratory (HyperComp), Department of Computer Technology and Communications. Escuela Politecnica de Caceres, University of Extremadura, Avenida de la Universidad sn, E-10002 Caceres, Spain
*
Author to whom correspondence should be addressed.
Remote Sens. 2020, 12(8), 1257; https://doi.org/10.3390/rs12081257
Submission received: 29 February 2020 / Revised: 7 April 2020 / Accepted: 15 April 2020 / Published: 16 April 2020
(This article belongs to the Special Issue Big Data in Remote Sensing for Urban Mapping)

Abstract

:
The storage and processing of remotely sensed hyperspectral images (HSIs) is facing important challenges due to the computational requirements involved in the analysis of these images, characterized by continuous and narrow spectral channels. Although HSIs offer many opportunities for accurately modeling and mapping the surface of the Earth in a wide range of applications, they comprise massive data cubes. These huge amounts of data impose important requirements from the storage and processing points of view. The support vector machine (SVM) has been one of the most powerful machine learning classifiers, able to process HSI data without applying previous feature extraction steps, exhibiting a robust behaviour with high dimensional data and obtaining high classification accuracies. Nevertheless, the training and prediction stages of this supervised classifier are very time-consuming, especially for large and complex problems that require an intensive use of memory and computational resources. This paper develops a new, highly efficient implementation of SVMs that exploits the high computational power of graphics processing units (GPUs) to reduce the execution time by massively parallelizing the operations of the algorithm while performing efficient memory management during data-reading and writing instructions. Our experiments, conducted over different HSI benchmarks, demonstrate the efficiency of our GPU implementation.

Graphical Abstract

1. Introduction

Recent advances in computer technology allowed for the development of powerful instruments for remote sensed data acquisition, lowering both the cost of their production and the cost of launching new Earth Observation (EO) missions. In particular, imaging spectroscopy (also known as hyperspectral imaging) [1] has attracted the attention of many researchers because of the great potential of hyperspectral images (HSIs) in characterizing the surface of the Earth by covering the visible, near-infrared and shortwave infrared regions of the electromagnetic spectrum. To exploit these data, multiple EO missions are now using imaging spectrometers, such as the Environmental Mapping and Analysis Programme (EnMAP) [2,3], or the Hyperspectral Precursor of the Application Mission (PRISMA) [4]. In this context, a high-dimensional stream of remotely sensed HSI data is now being generated, which can be applied to multiple tasks after being properly processed, such as managing of natural and environmental resources, urban planning and monitoring of agricultural fields, risk prevention and natural/human disaster management, among others.
However, HSI data processing faces important memory and computational requirements, due the large amount of information provided by EO airborne/satellite instruments. In fact, imaging spectrometers are able to collect hundreds of images over large areas along the Earth’s surface, gathering hundreds of narrow and continuous spectral bands along different wavelengths. As a result, each pixel in a HSI cube measures the reflection and absorption of electromagnetic radiation from ground objects into several spectral channels, creating a unique spectral signature for each material optically detected by the spectrometer [5]. This allows for a very accurate characterization of land cover surface, which is useful for the modeling and mapping of materials of interest but presents significant computational requirements. In particular, spectral-based classification algorithms need to process large amounts of spectral information in order to correctly assign a label (which corresponds to a fixed land cover class) to each pixel in the HSI scene, facing heavy computation burdens [6]. In this sense, the development of parallel implementations of traditional classification algorithms for fast and accurate processing is a requirement in order to efficiently process large HSI data repositories through proper management of computational resources.
With this aim, high performance computing (HPC) has become an efficient tool, not only for managing and analyzing the increasing amount of remotely sensed data that are currently being collected [7,8,9,10], but also for efficiently dealing with the high dimensionality of HSI data cubes [11,12]. From a computational point of view, most spectral-based HSI data classification algorithms exhibit inherent parallelism at three different levels [13]:
  • through HSI pixel vectors (coarse-grained pixel level parallelism),
  • through spectral-band information (fine-grained spectral level parallelism), and
  • through tasks (task-level parallelism).
HSI classification algorithms generally map nicely on HPC systems [14], with several HPC approaches (from coarse-grained and fine-grained parallelism techniques to complex distributed environments) already successfully implemented. For instance, distributed approaches based on commodity (homogeneous and heterogeneous) clusters [15,16], grid computing [17,18], and cloud computing techniques [19,20,21,22] provided good results in terms of reducing execution times, enabling an efficient processing of massive data sets on the ground segment. However, dedicated resources such as massively parallel clusters and networks of computers are generally expensive and hard to maintain, being cloud computing a more appropriate and cheaper approach (in addition to a fully distributed solution) due to its “service-oriented” computing and “pay-per-use” model [23]. Neither cluster nor cloud computing models allow on-board processing.
On the other hand, parallel solutions based on multicore CPUs [10,24,25], graphics processing units (GPUs) [14,26,27,28,29], and field programmable gate arrays (FPGAs) [30,31,32,33,34] were also successful, leading to near real-time performance in many HSI applications [35]. These platforms generally provide a cheaper solution when compared to large-scale distributed systems. Also, these approaches consume much less energy in comparison with cluster-based or distributed systems, being GPU systems (in particular, those with low-power consumption [36]) and FPGAs the ones currently offering the best performance/consumption ratio. Nevertheless, although both kinds of devices are appealing for on-board processing, FPGAs (despite their reconfigurability) are still more difficult to program than GPUs. In fact, GPUs are massively parallel programmable systems that can satisfy the extremely high computational requirements by HSI-based applications at low cost, with very good programmability. As a result, many HSI data processing algorithms were implemented on GPUs, including target and anomaly detection methods [37,38], and techniques for image compression [39], scene classification [40], and spectral unmixing [41].
In this paper, we develop a new GPU implementation of the traditional support vector machine (SVM) algorithm [42]. Traditionally considered to be a powerful classifier in a wide variety of fields, such as medical computing or text and image processing (among others), the SVM has been succefully adopted in HSI classification problems. It separates two classes of samples with the maximal margin by exploiting an optimal decision boundary [43]. As a result, SVMs have often been found to provide higher classification accuracies than other widely used pattern recognition techniques, such as the maximum likelihood, the random forest (RF) or the multi-layer perceptron (MLP) classifiers. Furthermore, SVMs appear to be particularly advantageous in heterogeneous environments, for which only a few training samples are available per class. The SVM has demonstrated its success in HSI data classification in a multitude of scientific works within the remote sensing field [44] due to its robust behaviour when dealing with high-dimensional data and its great ability to generalize to unseen data. However, SVMs exhibit a high computational cost, which strongly discourages its use when facing large-scale, real-time HSI classification problems.
To mitigate the aforementioned problem, several techniques were adopted to reduce the execution time of SVMs. For instance, in [45], Osuna et al. presented a decomposition algorithm by addressing sub-problems iteratively to enable tackling larger problems. In [46], Platt proposed the sequential minimal optimization (SMO) algorithm by decomposing the quadratic programming (QP) problem into a series of smaller QP subproblems that can be solved analytically, without the need for a time-consuming QP optimization. Ease of implementation is a typical feature of the SMO algorithm. In [47], Fan et al. developed a series of working set selection (WSS) heuristics, based on second order information, to achieve a faster convergence of the SMO. Recently, some researchers used the portability and excellent computing performance of GPUs to address the computational requirements of SVMs. In particular, Tan et al. [48] proposed a novel two-level parallel computing framework based on NVIDIA’s compute device unified architecture (CUDA) and OpenMP to accelerate SVM-based classification, while Li et al. [49] used GPU devices to improve the speed performance of SVM training and prediction stages.
In this work, we present a novel GPU implementation of SVMs that does not only reduces the execution time by massively parallelizing the operations of the algorithm on the GPU, but also performing efficient memory management during data-reading and writing operations [50]. We empirically evaluate the efficiency of our method considering different HSI real scenes, which comprise urban and agricultural land-cover information, and observed that it is able to maintain the precision in the results while significantly improving the computational performance.
The remainder of the paper is organized as follows. Section 2 presents the fundamentals of the SVM method adapted to HSI data. Section 3 describes the GPU implementation of the considered classifier. Section 4 validates the GPU implementation by comparing it with other existing implementations. Finally, Section 5 concludes the paper with some remarks and hints at plausible future research lines.

2. Support Vector Machines (SVMs): A Review

2.1. Linear SVM

Any HSI classification method y = f ( x ) can be described as the mapping function f : N n b N , where the domain is composed by those spectral pixels x i N n b = [ x i 1 , x i 2 , , x i n b ] that encode the reflectance values of the HSI cube X N n h × n w × n b , which is composed by n h · n w pixels (note that i = { 1 , , n h · n w } ) with n b spectral bands, while the co-domain is composed by the corresponding label-vectors y i [ 1 , 2 , , n c ] , being n c the number of different land cover categories. The final goal is to obtain the sample-label pairs for each HSI image element, { x i , y i } i = 1 n h · n w , obtaining as a result the final classification map Y N n h × n w .
In this context, the original SVM method [51] was proposed for binary classification problems, i.e., n c = 2 so y i { 1 , 1 } , i { 1 , , n h · n w } , where the two considered classes are linearly separable, with at least one hyperplane separating them (see Figure 1). This hyperplane is defined by its normal vector w R n b and the bias b R that defines the distance between the samples and the hyperplane. In this regard, the classifier can be defined as the discriminant function:
f ( x ) = w x + b ,
where f ( x ) = 0 is the hyperplane and the samples lie on the planes f ( x ) = ± 1 depending on the class they belong to. Both parameters w and b are estimated by the SVM during the training stage in order to maximize the distance between the closest sample and the hyperplane, i.e., with the aim of maximizing the margin m = 2 / | | w | | [42]. This maximization problem should be conveniently replaced by the equivalent minimization problem 1 2 | | w | | 2 . In the end, the SVM finds the optimal hyperplane along the space of samples by minimizing the quadratic optimization problem given by Equation (2):
min w , b 1 2 | | w | | 2 , subject   to y i w · x i + b 1 0 i { 1 , , n t }
where n t is the number of training samples, x i is the i-th sample, and y i is the correct output provided by the SVM, i.e., +1 for positive samples and -1 for negative ones. Considering the Lagrangian formulation [52], Equation (2) can be simplified by replacing the inequality ≥through a dual problem. In particular, employing the Lagrange multipliers α i 0 i { 1 , , n t } , we define the optimization problem given by Equation (3), which only depends on the set of Lagrange multipliers α = { α 1 , α 2 , , α n t } :
min α 1 2 i = 1 n t j = 1 n t y i y j x i · x j α i α j i = 1 n t α i , subject   to   α i 0 and i = 1 n t y i α i = 0 i { 1 , , n t } ,
where x i and x j are two samples of the training set with labels y i and y j and Lagrange multipliers α i and α j , respectively. Equation (3) can be computed using QP methods, such as the SMO [53], which breaks the problem into a series of smallest and equivalent QP problems, solving each one analytically.
With this problem representation, the normal vector w can be derived as the sum of all the samples that belong to the training set and modified by their corresponding Lagrange multiplier together with the associated class of the sample, w = i = 1 n t α i y i x i , while the bias is obtained for some multiplier α j > 0 as b = w x j y j . Moreover, the discriminant function f ( x ) provided by Equation (1) can be solely described in terms of Lagrange multipliers and samples as:
f ( x ) = i = 1 n t α i y i ( x i · x ) + b
It is worth noting that each Lagrange multiplier α i measures the relevance of the sample x i when calculating the hyperplane and classifying the data, so that samples with α i = 0 will be ignored, while samples with α i > 0 will be considered to be support vectors. Also, from Equation (4), we can observe that, to perform the classification of any sample x , a significant volume of memory will be required, consuming also a large amount of computational resources to compute the dot products needed.

2.2. Linear SVM for Linearly Nonseparable Data Classification

The linear SVM has a limitation: with linearly nonseparable data classification there may not be a hyperplane, so the solution becomes infinite. To overcome this drawback, the linear soft-margin SVM implements a new cost function, which takes into account the original margin maximization and including a loss term that penalizes misclassified data, as it can be observed in Equation (5):
min w , b 1 2 | | w | | 2 + C i = 1 n t δ i , subject   to y i w · x i + b 1 + δ i 0 and δ i 0 ,
where δ i 0 are some slack variables that are introduced into the inequality restriction y i w · x i + b 1 + δ i 0 to relax the condition that all the samples that belongs to the same class lie on the same side of the hyperplane. Thus, i = 1 n t δ i can be interpreted as some measure of the amount of misclassifications. Also, a regularization constant C [ 0 , ) is included to control the weight of minimizing 1 2 | | w | | 2 while penalizing solutions with very large δ i . In fact, C can be considered to be a parameter that adjusts the robustness or flexibility of the model, depending its value.
Once more, the Lagrangian formulation can be applied to translate the primal problem defined by Equation (5) into the dual problem of Equation (3), by simply introducing the box constraint C α i 0 . The slack variables δ i do not affect the optimization problem.

2.3. Kernel SVM for Non-Linear Data Classification

As an alternative to the soft-margin implementation, the SVM method can be adapted to perform non-linear classification by including a kernel [54,55] within the classifier. The main idea is to map the original feature vectors x i N n b into a higher dimensional Euclidean space, denoted as H , by using a non-linear vector function Φ : N n b H , as we can observe in Figure 2. In this context, the optimal margin problem can be posed in the space H by replacing the original inner product between the samples x i · x with the transformed vectors Φ ( x i ) · Φ ( x ) , so the discriminant function will be defined as:
f ( x ) = i = 1 n t α i y i Φ ( x i ) · Φ ( x ) + b Replacing K ( x i , x ) = Φ ( x i ) · Φ ( x ) f ( x ) = i = 1 n t α i y i K ( x i , x ) + b
Equation (6) can be easily simplified by assuming there is such a kernel function that K ( x i , x ) = Φ ( x i ) · Φ ( x ) . This assumption allows avoiding the computation of the inner product between transformed vectors, significantly reducing the computational consumption of the algorithm. To ensure the existence of the underlying map Φ , the kernel function is positively defined [51] by satisfying Mercer’s conditions [56].
Finally, the solution of the dual problem can be also meaningfully simplified as:
min α 1 2 i = 1 n t j = 1 n t y i y j K x i x j α i α j i = 1 n t α i , subject   to C α i 0   and i = 1 n t y i α i = 0
In fact, the kernel function usually measures the distance between the input sample x i and the other training samples x j . Table 1 shows some of the most widely used kernels. In particular, the RBF kernel has quite interesting properties. It can be simplified as exp γ | | x i x j | | 2 , where γ controls the model’s performance, in the sense that a small value of γ will provide classification solutions with low bias and high variance, while a high value of γ will give solutions with high bias and low variance.

2.4. Multi-Class SVM

The aforementioned SVM implementations can be generalized to develop multi-class classification by combining several binary classifiers [42]. In this regard, there are two main strategies that can be adopted to perform multi-class classification [57]: (i) constructing and combining several binary classifiers, and (ii) considering all data in one optimization formulation.
Regarding the first approach, it is usual to find one-against-all, one-against-one and decision directed acyclic graph (DDAG)-based approaches for multi-class SVM in the available literature [58,59,60]. The first approach constructs n c pairwise classifiers, where the k-th model is trained with those samples that belong to the k-th class by assigning them a positive label (while the rest of the samples are paired with negative labels). The classification procedure applies the n c SVM models to each sample x i , selecting the label of the classifier f k ( x ) with the largest margin:
y i = arg max k f k ( x i )
As a result, n c n t -variable QP-problems have to be solved, which implies that training time scales linearly with n c . Instead of developing n c classifiers, the one-against-one approach adopts a pairwise decomposition strategy, implementing n c ( n c 1 ) / 2 binary models on all pairs of training data. Each model f k t is trained on data that belong to two different classes, k and t, ignoring all the other classes. In this way, those samples that belong to class k are labeled as positive, while the samples that belong to class t correspond to the negative label (note that f t k = f k t ). The final label is selected by majority voting strategy:
y i = arg max k t f k t ( x i )
It is noteworthy that the one-against-one approach exhibits a much more robust behaviour when classifying imbalanced datasets; however, it is also more computationally expensive than the one-against-all for simpler problems, such as those with fewer classes.
Finally, the DDAGSVM also solves n c ( n c 1 ) / 2 binary models during the training stage, implementing a rooted binary directed acyclic graph with also n c ( n c 1 ) / 2 internal nodes (each one is a binary SVM) and n c leaves (i.e., the predicted class) during the testing stage.

3. GPU-Accelerated SVM for HSI Data Classification

3.1. Previous Works and Proposal Overview

Despite obtaining great accuracy results, the training and prediction steps of the SVM algorithm are very expensive from a computational point of view, especially for large and complex problems. Furthermore, it is noteworthy that kernel-based techniques can greatly simplify the computations, thus reducing the execution times. However, these methods need to use all the data to properly calculate the kernel, which generally requires large memory consumption [61]. In particular, HSI data offers a great challenge to the SVM classifier, mainly due to (both) the large spectral dimensionality and volume of data to be processed. As noted in Section 1, there were some previous works in the literature focused on optimizing and parallelizing the execution of the SVM. However, there are few efforts in the recent literature oriented at the exploitation of GPU-based implementations for HSI data classification. From the current state-of-the-art in the literature, we can highlight two works. Tan et al. [48] combined the NVIDIA Compute Unified Device Architecture (CUDA) [50] and OpenMP [62] to optimize the classification problem, and implemented as a multi-class SVM following one-against-one strategy. In this sense, they were able to run the one-against-one strategy in parallel by implementing an OpenMP-based approach, while each single WSS-based binary SVM classifier was parallelized into GPUs. In addition, the calculation of every kernel element K ( x i , x j ) was done by the GPU, being the constant values stored into the shared memory of the GPU device, following column-major order, and sent (when needed) through a broadcasting mechanism. Also, Wu et al. [63] developed a GPU-parallelized SVM method based on composite kernels (SVMCK) [64] and able to integrate spatial and spectral information, with the aim of improving the accuracy in HSI data classification tasks. In this context, they carried the calculations of the composite kernel matrix in the GPU. The resulting K ( x i , x j ) were copied to the host (CPU), and then stored in the available memory of the host. The rest of the computations were carried out in the CPU.
In this context, we propose a new parallel version of the SVM algorithm for high dimensionality data. Our proposal can be efficiently generalized to be applied in other fields, being our parallel SVM adapted to the specific case of remote sensing HSI data classification. In particular, we propose a parallel implementation of a multi-class SVM, following one-against-one strategy, with the RBF kernel. It must be noted that each single binary SVM has been solved through the SMO method. Furthermore, the SMO solver (during the training stage) and the full inference stage were paralleled into a general purpose GPU (GPGPU) framework, employing the NVIDIA CUDA platform. Our goal is not only to achieve a good speed up ratio in comparison with other traditional SVM implementations, but also to develop a scalable implementation with regards to both the amount of data and the spectral dimensionality of the data. In this regard, the implemented approach performs the algorithmic parallelization of the SVM operations during the training and inference stages, implementing at the same time a series of optimizations in the memory read and write operations, with the aim of reducing the number of communications between the host and the device, and to reduce the memory latency. Also, it should be noted that a data-oriented parallelization strategy has been adopted, parallelizing the operations related to data management and reusing as much information as possible between different iterations. Regarding to this, a batch-based approach has been implemented within the SMO with the aim of solving multiple sub-problems in parallel. Particular attention has been paid to the calculation of the kernel matrix, as it consumes a massive amount of computational and storage resources. In the following sections we will provide details about our new implementation.

3.2. CUDA Platform

To obtain an accelerated version of the SVM for HSI data classification, a GPGPU-based implementation has been employed using NVIDIA CUDA platform [50]. The main goal is to exploit the full computing power of GPU devices (using the general-purpose parallel computing platform and the programming model offered by CUDA) to solve many complex computational problems in a more efficient way. Specifically, we solve n c binary SVM instances (with the RBF kernel) that compose our multi-class SVM algorithm for remote sensing HSI data classification. In this context, we consider the GPU device as a set of stream multiprocessors, composed by several cores with a 3-level hierarchical memory system: (i) several fast registers that are accessible by each stream processor, (ii) a multiprocessor-level memory that is shared between all the cores that compose the multiprocessor, and (iii) a global memory shared between all the multiprocessors.
As we can observe in Figure 3, CUDA maintains the hierarchy of memories at the same time that provides a programming model which abstracts the system of cores and multiprocessors into a three-level logical representation of threads, blocks and grids, where the thread is the minimum processing unit (that can be grouped into 1D/2D/3D blocks of warps), and the blocks can be organized into a 1D/2D grid environment. In this sense, the programmer determines the logical representation by organizing the structure of the threads in the grid that best fits the data structure handled by the problem to be parallelized. Underneath, CUDA maps this representation to the actual hardware resources. This distinction between physical and logical representations allows CUDA-based algorithms to be ported to any NVIDIA GPU platform, improving their scalability and computing performance as new and higher-capacity GPUs are released to the market. In addition, it should be highlighted that each thread has access to three different memory spaces: its local memory, the memory shared by all the threads of the same block, and the global memory shared by all the threads of the grid. Therefore the correct handling of the threads and their accesses to memory are fundamental for the correct and optimal parallelization of the problem.
Through the organization and synchronization of the threads, we can develop the parallel processing of the operations carried out by the SVM method. In particular, we can differentiate between those operations performed during the training stage of the SVM and those operations performed during the inference step.

3.3. Parallel SMO during the Training Stage

3.3.1. Previous Concepts about the SMO Algorithm

As pointed out before, we adopted the one-against-one approach for the multi-class SVM to perform remote sensing HSI data classification. Therefore n c ( n c 1 ) / 2 binary models must be properly trained to extract the support vectors and derive the corresponding Lagrange multipliers for each classifier. In particular, the binary models were trained using a decomposition method to solve the convex optimization problem defined by Equation (7). It is noticeable that the main difficulty when solving Equation (7) is how to calculate the kernel matrix K R n t × n t that stores the kernel values K ( x i , x j ) , i [ 1 , n t ] and j [ 1 , n t ] , as:
K = K ( x 1 , x 1 ) K ( x 1 , x 2 ) K ( x 1 , x n t ) K ( x 2 , x 1 ) K ( x 2 , x 2 ) K ( x 2 , x n t ) K ( x n t , x 1 ) K ( x n t , x 2 ) K ( x n t , x n t )
Usually, K is a dense matrix and may be too difficult to be efficiently stored and handled if n t is too large. To deal with this limitation, some decomposition methods were designed [46,47,65,66,67,68,69,70,71,72,73,74] to break down the problem into several smaller and easier to handle sub-problems, where small subsets of Lagrange variables α are modified by employing some columns of K instead of the entire matrix. In particular, the SMO [46] algorithm has been considered in this work.
The SMO is a simple and iterative algorithm that allows for the quick and effective resolution of very large QP problems (such as those involved in the SVM calculations) by decomposing such overall QP problem into smaller QP sub-problems [45], which are solved analytically without the need for numerical optimization. It solves the optimization problem described by Equations (3) and (7), by solving the smallest possible optimization problem at every step until the optimal condition of the SVM classifier is reached. As the linear equality constraint i = 1 n t y i α i = 0 involves the Lagrange multipliers α i , the smallest possible optimization problem will involve two such multipliers, in the sense that, if we change one α t by an amount in either direction, then the same change must be applied to another α l in the opposite direction. That is, α t and α l should be on the same line in order to maintain the constraint (see Figure 4):
y t α t + y l α l = y t α ^ t + y l α ^ l , subject   to α t 0 and α l C
In this way, the SMO algorithm comprises two main components: (i) a heuristic for selecting the pair of Lagrange multipliers to be optimized, and (ii) an analytic method for solving those multipliers. In this sense, in each iteration the SMO algorithm heuristically chooses two Lagrange multipliers α t and α l at every step to jointly optimize, then it analytically obtains the new optimal values α ^ t and α ^ l , and finally it updates the SVM to reflect the new values.
Focusing on the heuristic procedure, the SMO applies two heuristic searches, one for each Lagrange multiplier. The first multiplier α t is chosen by iterating over the entire training set, looking for those samples that violate the Karush-Kuhn-Tucker (KKT) conditions [75] that help to find an optimal separating hyperplane. In particular, the KKT conditions for Equation (7) are:
α i = 0 y i ( w x i + b ) 1 , C > α i < 0 y i ( w x i + b ) = 1 , α i = C y i ( w x i + b ) 1
where y i is the correct SVM output and ( w x i + b ) is the current output of the SVM for the i-th sample x i , i [ 1 , n t ] . Any α i that satisfies the KKT conditions will be an optimal solution for the QP optimization problem defined by Equation (7). On the contrary, any α i that violates the KKT conditions will be eligible for optimization, so the SMO’s goal is to iterate until all these conditions are satisfied within a tolerance threshold (in our case, this tolerance has been set to 0.001). Once the first α t has been chosen, the second Lagrange multiplier α l is selected in order to maximize the size of the step taken during joint optimization. To do this, the SMO method implements an optimality indicator vector E = [ E 1 , E 2 , , E n t ] , where each E j is the optimality indicator of the j-th training sample, i.e., the classification errors on the j-th sample:
E j = f ( x j ) y j = i = 1 n t α i y i K ( x i , x j ) + b y j
Related to this, α t and α l can be selected by looking for those samples x t and x l that have the maximum and minimum optimality indicators that maximize | E t E l | , so if E t is positive, the SMO will choose a sample with minimum E l , while if E t is negative, the SMO will select a sample with maximum E l . The desired indexes t and l can be directly obtained by computing Equation (13) [67]:
t = arg min i { E i | x i X u p p e r } l = arg max i ( E t E i ) 2 μ i | E t < E i , x i X l o w e r
where μ i = K ( x t , x t ) + K ( x i , x i ) 2 K ( x t , x i ) , E t and E i are the optimality indicator of samples x t and x i , respectively, and X u p p e r = X 1 X 2 X 3 and X l o w e r = X 1 X 4 X 5 are two data subsets, where each X is composed by the following training samples:
X 1 = { x i | x i X , 0 < α i < C } X 2 = { x i | x i X , y i = + 1 , α i = 0 } X 3 = { x i | x i X , y i = 1 , α i = C } X 4 = { x i | x i X , y i = + 1 , α i = C } X 5 = { x i | x i X , y i = 1 , α i = 0 }
Once both Lagrange multipliers have been obtained, the SMO method will compute their optimal values α ^ t and α ^ l with the aim of obtaining the optimal class-separating hyperplane. In particular, it begins by calculating α ^ l , whose feasible values are framed into the constraint U α ^ l V to meet the original C α l 0 , being U and V two boundaries defined as:
I f y t y l U = max { 0 , α l α t } V = min { C , C α t + α l }
I f y t = y l U = max { 0 , α t + α l C } V = min { C , α t + α l }
Besides, the optimal α l ^ value within the range [ U , V ] will be obtained as:
α ^ l = V , if α ˜ l > V α ˜ j , if V α ˜ l U U , if α ˜ l < U being α ˜ l = α l + y l ( E t E l ) μ
where μ = K ( x t , x t ) + K ( x l , x l ) 2 K ( x t , x l ) and E t and E l are the classification errors on the t-th and l-th training samples respectively. Once α ^ l has been obtained, an optimal α ^ t is easily calculated as:
α ^ t = α t + y t y l α l α ^ l
Once α ^ t and α ^ l have been obtained, the SMO updates the bias threshold b such that the KKT conditions are satisfied for the t-th and l-th samples. In this sense, three cases can occur, as we can observe in Equation (18):
b ^ = b ^ 1 if C > α t > 0 b ^ 2 if C > α l > 0 ( b ^ 1 + b ^ 2 ) 2 otherwise
where b ^ 1 and b ^ 2 are defined as follows:
b ^ 1 = b E t y t ( α ^ t α t ) K ( x t , x t ) y l ( α ^ l α l ) K ( x t , x l ) b ^ 2 = b E l y t ( α ^ t α t ) K ( x t , x l ) y l ( α ^ l α l ) K ( x l , x l )
Finally, the SMO updates the SVM. For each training sample x i , the SMO updates its E i using the following Equation (19):
E i = E i + ( α ^ t α t ) y t K ( x t , x i ) + ( α ^ l ) y l K ( x l , x i )
The full procedure is repeated until the optimal condition is reached, i.e., E t E m a x , where E m a x acts as a threshold, E m a x = max { E i | x i X l o w e r } .

3.3.2. CUDA Optimization of SMO Algorithm

The SVM starts by dividing the HSI scene into training and inference subsets. We can consider the training set as a collection of instances and their associated labels, i.e., D t r a i n = { X , Y } . The training instances are represented by a 2D-matrix X N n t × n b composed by n t training vectors, where each x i N n b = [ x i , 1 , x i , 2 , , x i , n b ] comprises n b spectral bands, while the training labels are stored into the matrix Y N n t × n c , being y i = [ y i , 1 , y i , 2 , , y i , n c ] the corresponding label of sample x i in one-hot encoding, where n c indicates the number of different land cover classes. The training stage starts by creating the different single binary SVMs in a sequential way, confronting each class i to each different class j, i , j [ 1 , n c ] . Each binary SVM was optimized by employing a parallel SMO solver.
The parallel SMO solver begins by creating its working set B. Traditionally, B is of size two; however, our proposal implements a bigger working set to solve in parallel multiple subproblems of SMO in a batch. Moreover, the proposed implementation precomputes all the kernel values for the current working set, storing them as a data buffer on the device’s global memory, in order to reduce the high latency memory accesses and also to avoid a large number of small read/write operations in the CPU memory. This implies that in each iteration of the SMO algorithm, the current working batch B will be updated, being n B Lagrange multipliers optimized. This allows us to compute (at once) n B rows of the kernel matrix K , performing a more efficient use of the GPU by reducing the number of accesses to the device’s global memory (which is much slower than the shared memory) and enabling the re-use of kernel information (which in turn, reduces repeated kernel value computations). Taking this into account, our parallel SMO solver can be divided into three different steps.
During the step 1, the parallel SMO solver looks for the n B extreme training instances which can potentially improve the SVM the most according to Equation (9). This is parallelized by applying consecutive parallel reductions [76] over the training samples. For each working set, the optimality indicators of the training samples are sorted in ascending order, selecting the first n B / 2 and the last n B / 2 training samples to optimize the corresponding n B Lagrangre multipliers. During the reduction, one thread per sample takes the data from the global device memory to the block shared memory (in a coalesced way) to perform a fast execution. It must be noted that shared memory is around 7 times faster than global memory. Once synchronized, the threads operates over the data, where each one takes and compares the optimality indicators of the samples from the batch, choosing the smaller or the larger one depending on the desired Lagrange multiplier through the application of the corresponding heuristics given by Equation (13).
Once the n B Lagrange multipliers have been selected, the improvements of the Lagrangian multiplier pair α t and α l are sequentially obtained by one GPU thread per pair (step 2). It is worth noting from previous Equations (13), (16) and (19) that, during training stage, the same kernel values may be used in indifferent iterations. This implies that some kernel values can be stored into the GPU memory with the aim of reducing the high latency memory and avoiding a large number of small read/write operations from the CPU memory. In this sense, all the kernel values related to the n B Lagrange multipliers are obtained, computing n B rows of the kernel matrix K through parallel matrix multiplications [77,78]:
K B = K ( x 1 , x 1 ) K ( x 1 , x 2 ) K ( x 1 , x n t ) K ( x 2 , x 1 ) K ( x 2 , x 2 ) K ( x 2 , x n t ) K ( x n B , x 1 ) K ( x n B , x 2 ) K ( x n B , x n t )
In particular, the RBF-based K B matrix is computed in parallel by the GPU and stored into the device’s global memory. Algorithm 1 provides the pseudo-code of the parallel kernel RBF function, which were implemented following the approximate RBF kernel form:
K ( x i , x j ) = exp γ x i x j 2 = exp γ x i 2 + x j 2 2 x i x j
In this context, the input parameters x x and x 2 x contain the | | x i | | 2 , i [ 1 , , n B ] and x i x j , i , j [ 1 , , n B ] values respectively, which were previously computed through parallel matrix multiplications [78]. Then, n B threads compute the desired n B rows of the kernel matrix K B :
Algorithm 1 Parallel Kernel RBF for HSI classification
Require: x x matrix of | | x i | | 2 values,
x 2 x matrix of x i x j values,
K B resulting kernel matrix,
n B number of rows,
n t number of training samples.
 
i = b l o c k I d x . x b l o c k D i m . x + t h r e a d I d x . x
if i d x < n B then
   j = 0
  while j < n t do
    K B [ i , j ] = exp ( γ ( x x [ i ] + x x [ j ] 2 x 2 x [ i , j ] ) )
    j + +
  end while
end if
During the training stage, the kernel values are extracted from the global memory and gathered into the shared memory as they are needed, avoiding repeated computations.
Finally, the parallel SMO solver updates the optimality indicator vector E of the training instances by launching n t GPU threads that compute Equation (14) in parallel (step 3).
These training steps are sequentially repeated until the optimality condition is met or when the SVM classifier is not able to improve.

3.4. Parallel Classification during the Inference Stage

The proposed one-against-all SVM for HSI remote sensing data multi-class classification provides the corresponding label y i of a given test sample x i by applying Equation (8) during the inference stage. In this sense, Equation (8) is computed in parallel, making use of the kernel values in K to reduce the latency memory and avoid the repeated computations. Also, from Equation (6) we can observe that K ( x i , x j ) is the same that K ( x j , x i ) , where i and j are related to the number of test samples and support vectors, respectively. In this sense, we can reduce the number or read/copy operations into the GPU memory following two strategies. On the one hand, if the number of test samples is larger than the number of support vectors, we read all the rows of the kernel matrix K that correspond to the support vectors by reading the indexes respected to j. On the other hand, if the number of support vectors is larger than the number of test samples, those rows that correspond to the test samples will be read by reading the indexes respected to i.

4. Experimental Results

4.1. Experimental Environment

To evaluate the performance and the benefits of the proposed parallel SVM for HSI remote sensing data classification, several implementations of the proposed classifier were developed and tested over two different hardware platforms:
  • Platform 1: it is composed by an Intel Core Coffee Lake Refresh i7-9750H processor, 32 GB of DDR4 RAM with 2667 MHz, and an NVIDIA GeForce RTX 2070 with 8 GB of RAM, graphic clock at 2100 MHz and 14,000 MHZ of memory transfer rate. It is equipped with 2304 CUDA cores. These processors were named CPU0 and GPU0.
  • Platform 2: Intel i9-9940X processor, 128 GB of DDR4 RAM with 2100 MHz, and an NVIDIA GTX 1080Ti with 11 GB of RAM, 2037 MHz of graphic clock and 11,232 MHz of memory transfer rate. It is equipped with 3584 CUDA cores. These processors were named CPU1 and GPU1.
CPU1 will serve as the baseline for the performance comparisons that were conducted during the experimentation due to its characteristics. Moreover, both environments use Ubuntu 18.04.3 x64 as operating system.

4.2. Hyperspectral Datasets

Our experiments were conducted over six different and well-known HSI scenes: Indian Pines and Big Indian Pines, Pavia University, Salinas Valley and University of Houston. Figure 5 shows the ground truth and the number of pixels per class for each image. Below, the characteristics of each scene are provided.
  • The first dataset is known as Indian Pines, which was collected by the Airborne Visible Infra-Red Imaging Spectrometer (AVIRIS) [79] over the Indian Pines test site in North-western Indiana, which is characterized by several agricultural crops and irregular forest and pasture areas. It has 145 × 145 pixels, each of which has 224 spectral reflectance bands covering the wavelengths from 400 nm to 2500 nm. We remove the bands 104–108, 150–163 and 220 (water absorption and null bands), and keep 200 bands in our experiments. This scene has 16 different ground-truth classes (see Figure 5).
  • Big Indian Pines is a larger version of the first dataset, which has 2678 × 614 pixels with wavelengths ranging from 400 nm to 2500 nm. We also remove the water absorption and null bands, retaining 200 spectral bands in our experiments. This scene has 58 ground-truth classes (see Figure 5).
  • The third dataset is the Pavia University scene, which was collected by the Reflective Optics Spectrographic Imaging System (ROSIS) [80] during a flight campaign over Pavia, nothern Italy. In this sense, it is characterized by being an urban area, with areas of buildings, roads and parking lots. In particular, the Pavia University scene has 610 × 340 pixels, and its spatial resolution is 1.3 m. The original pavia dataset contains 115 bands in the spectral region of 0.43–0.86 μ m. We remove the water absorption bands, and retain 103 bands in our experiments. The number of classes in this scene is 9 (see Figure 5).
  • The fourth dataset is Pavia Centre and was also gathered by ROSIS sensor. It is composed by 1096 × 1096 pixels and 102 spectral bands. This scene also has 9 ground-truth classes from an urban area (see Figure 5).
  • The fifth dataset is Houston University [81], which was acquired by the Compact Airborne Spectrographic Imager (CASI) sensor [82] over the Houston University campus in June 2012, collecting spectral information from an urban area. This scene has 114 bands and 349 × 1905 pixels with wavelengths ranging from 380nm to 1050nm. It comprises 15 ground-truth classes (see Figure 5).
  • Finally, the sixth dataset is Salinas Valley, which was also acquired by AVIRIS sensor over an agricultural area. It has 512 ×217 pixels and covers Salinas Valley in California. We remove the water absorption bands 108–112, 154–167 and 224, and keep 204 bands in our experiments. This scene contains 16 classes (see Figure 5).

4.3. Performance Evaluation

With the aim of evaluating the computational performance obtained by the proposed GPU implementation of SVMs for HSI data classification, and making a thorough analysis of the implemented classifiers also in terms of accuracy, several experiments were carried out:
  • The first experiment focuses on the classification accuracy obtained by our GPU implementation as compared to a standard (CPU) implementation in LibSVM [83]. In particular, proposed GPU-SVM was compared with its CPU counterpart, the random forest (RF) [40] and the multinomial logistic regression (MLR) [40].
  • The second experiment focuses on the scalability and speedups achieved by the GPU implementation with regards to the CPU implementation, from a global perspective. As we pointed before, CPU1 will be considered to be the baseline due its characteristics that make it the slowest device.
  • The third and last experiment focuses on some specific aspects of the GPU implementation, including data-transfer times.
In the following, we describe in detail the results obtained in each of the aforementioned experiments.

4.3.1. Experiment 1: Accuracy Performance

In this experiment, we use the fixed training and test data for the AVIRIS Indian Pines image [displayed in Figure 6c,d], the ROSIS Pavia University image [displayed in Figure 7c,d] and the University of Houston image [Figure 8c,d]. These fixed training and test sets are publicly available online from the IEEE Geoscience and Remote Sensing Society (GRSS) data and algorithm evaluation website (DASE) at http://dase.grss-ieee.org. The main rationale for using these fixed training sets is that all algorithms can be compared in a fair way, regardless of random splits of training and test sets.
The obtained classification results can be graphically observed in Figure 6 and Figure 8. In particular, Figure 6a, Figure 7a and Figure 8a respectively show false color compositions of the aforementioned three images, while Figure 6b, Figure 7b and Figure 8b show the entire ground-truth data for each image. The results obtained by two well-known classifiers: RF [see Figure 6e, Figure 7e and Figure 8e] and the multinomial logistic regression (MLR) [see Figure 6f, Figure 7f and Figure 8f] are reported, together with the classification maps obtained by the CPU (LibSVM) implementation of the SVM [see Figure 6g, Figure 7g and Figure 8g] and our GPU implementation of SVM [see Figure 6h, Figure 7h and Figure 8h]. Although all classification maps have the typical “salt and pepper” noise of spectral classifiers, we can see that SVM classifiers achieve a higher quality map by classifying certain regions of the scenes more precisely, for instance the Soybeans-notill and Oats land-covers in Indian Pines scene or the the tree-line areas of Pavia University. Moreover, those classification maps obtained by CPU-based SVM and GPU-based SVM are quite similar.
The individual class accuracies are reported on Table 2. Also the overall (OA), average (AA) accuracies and kappa coefficient are shown for each HSI dataset. We can observe that the CPU and GPU implementations of the SVM classifier reach the best OA and AA percentages in all datasets, exhibiting also the best kappa value. In particular, as shown in all cases, the OA obtained by our GPU implementation of SVM is very similar to that obtained by the corresponding CPU implementation, and superior to that obtained by other well-known classifiers such as RF and MLR.

4.3.2. Experiment 2: Scalability and Speedup

In this experiment, we used randomly selected training/test sets to evaluate the scalability of our GPU implementation as the amount of training becomes larger. In particular, 1%, 3%, 5%, 10%, 15%, 20%, 25%, 30%, 40%, 60% and 80% of training data were considered for the University of Pavia, Pavia Center, Salinas Valley and Big Indian Pines datasets.
Figure 9 reports the processing times (and speedups) measured in the considered CPUs (CPU0 in the first hardware platform and CPU1 in the second platform) and GPUs (GPU0 in the platform 1 and GPU1 in the platform 2) as the percentage of training increases. As we can observe in Figure 9, for all the considered images in this experiment (University of Pavia, Pavia Center, Salinas Valley and Big Indian Pines), the processing times increase as the percentage of training samples increases. This is expected, as the complexity of the classification problem becomes larger. However, we can also observe that the speedups obtained become more significant with the training size, particularly for the larger images. For instance, with University of Pavia (about 34,220 pixels when training with 80% of data) the implementation is able to reach a speedup of x6 with the fastest GPU, while with Big Indian Pines (around 1,315,433 samples when considering 80% of training data) is x140 faster than the baseline. This indicates that our newly developed GPU implementation of SVM scales with complexity and problem size, which is a highly desirable feature in parallel implementations.

4.3.3. Experiment 3: GPU Transfer-Memory and Kernel Runtimes

In this experiment, we fix the training percentage to 40% and study the memory transfers required by the proposed GPU implementation for different HSI datasets over GPU0 and GPU1 devices. Here, we specifically focus on the two largest datasets (Pavia Center and Big Indian Pines) and report the memory times consumed by the considered GPU platforms. Table 3 provides a breakdown of the obtained processing times for different images and training percentages, indicating the transfer times from host to device (HtoD), device to host (DtoH), device to device (DtoD) and the time taken by the kernels. The times reported on Table 3 indicate that most of the time spent by our parallel implementation is consumed by kernels and not by data transfers, despite these transfers become more significant as the training percentage becomes larger (this does not affect scalability, as reported in Figure 9. We can graphically observe these results in Figure 10, where the amount of GPU video memory required by our implementation is very small in all cases (below 5.9%), which indicates that our implementation is efficient in terms of memory usage, even with large percentages of training.

5. Conclusions and Future Lines

This work presented a new GPU implementation of the well-known SVM technique in the context of HSI remote sensing data classification, with the goal of reducing computation times by massively parallelizing operations across GPU threads, and reducing memory latencies by efficiently handling data read/write operations. Three experiments were conducted over six widely used and real HSI datasets, providing very heterogeneous information in terms of content (different land-cover classes extracted from agricultural and urban environments) and amount of data (different spatial sizes with also different spatial resolutions). These experiments show the effectiveness of our accelerated SVM implementation, which not only scales with data volume, but also with training percentage. This reveals that the acceleration is more effective when the classification problem is more complex. Moreover, the study of data transfer times and algorithmic computation times shows an efficient use of memory, where the kernel computations dominate the parallel processing times. In addition, the proposed GPU implementation achieves good performance in comparison with other popular SVM implementations, achieving similar accuracy results but with faster performance.
As with any new approach, there are some unresolved issues that may present challenges over time. In the future, we will improve our SVM implementation in order to reduce the accesses to global device memory and make a more efficient use of CUDA streams. In addition, we will use our accelerated SVM in conjunction with other techniques for HSI processing and classification (e.g., supervised and semi-supervised techniques) with the aim of improving even more the obtained classification results.

Author Contributions

The authors have contributed as equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

This work has been supported by the Spanish Ministry (FPU15/02090), Junta de Extremadura, Ref. GR18060 and the European Union’s Horizon 2020 research and innovation programme under grant agreement No. 734541 (EOXPOSURE).

Acknowledgments

We gratefully thank the Associate Editor and the five Anonymous Reviewers for theiroutstanding comments and suggestions, which greatly helped us to improve the technical quality and presentationof our work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Goetz, A.F.; Vane, G.; Solomon, J.E.; Rock, B.N. Imaging spectrometry for earth remote sensing. Science 1985, 228, 1147–1153. [Google Scholar] [CrossRef] [PubMed]
  2. Heldens, W.; Heiden, U.; Esch, T.; Stein, E.; Müller, A. Can the future EnMAP mission contribute to urban applications? A literature survey. Remote Sens. 2011, 3, 1817–1846. [Google Scholar] [CrossRef] [Green Version]
  3. Guanter, L.; Kaufmann, H.; Segl, K.; Foerster, S.; Rogass, C.; Chabrillat, S.; Kuester, T.; Hollstein, A.; Rossner, G.; Chlebek, C.; et al. The EnMAP spaceborne imaging spectroscopy mission for earth observation. Remote Sens. 2015, 7, 8830–8857. [Google Scholar] [CrossRef] [Green Version]
  4. Pignatti, S.; Palombo, A.; Pascucci, S.; Romano, F.; Santini, F.; Simoniello, T.; Umberto, A.; Vincenzo, C.; Acito, N.; Diani, M.; et al. The PRISMA hyperspectral mission: Science activities and opportunities for agriculture and land monitoring. In Proceedings of the 2013 IEEE International Geoscience and Remote Sensing Symposium-IGARSS, Melbourne, Australia, 21–26 July 2013; pp. 4558–4561. [Google Scholar]
  5. Chang, C.I. Hyperspectral Imaging: Techniques for Spectral Detection and Classification; Springer: Berlin/Heidelberg, Germany, 2003. [Google Scholar] [CrossRef]
  6. Christophe, E.; Michel, J.; Inglada, J. Remote sensing processing: From multicore to GPU. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2011, 4, 643–652. [Google Scholar] [CrossRef]
  7. Messmer, P. High-Performance Computing in Earth- and Space-Science: An Introduction. In Applied Parallel Computing. State of the Art in Scientific Computing, Proceedings of the 7th International Workshop, PARA 2004, Lyngby, Denmark, 20–23 June 2004; Revised Selected Papers; Dongarra, J., Madsen, K., Waśniewski, J., Eds.; Springer: Berlin/Heidelberg, Germany, 2006; pp. 527–529. [Google Scholar] [CrossRef]
  8. Plaza, A.J.; Chang, C.I. High Performance Computing in Remote Sensing; CRC Press: Boca Raton, FL, USA, 2007. [Google Scholar]
  9. Lee, C.A.; Gasster, S.D.; Plaza, A.; Chang, C.I.; Huang, B. Recent developments in high performance computing for remote sensing: A review. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2011, 4, 508–527. [Google Scholar] [CrossRef]
  10. León, G.; Molero, J.M.; Garzón, E.M.; García, I.; Plaza, A.; Quintana-Ortí, E.S. Exploring the performance-power-energy balance of low-power multicore and manycore architectures for anomaly detection in remote sensing. J. Supercomput. 2015, 71, 1893–1906. [Google Scholar] [CrossRef] [Green Version]
  11. Plaza, A.; Benediktsson, J.A.; Boardman, J.W.; Brazile, J.; Bruzzone, L.; Camps-Valls, G.; Chanussot, J.; Fauvel, M.; Gamba, P.; Gualtieri, A.; et al. Recent advances in techniques for hyperspectral image processing. Remote Sens. Environ. 2009, 113, S110–S122. [Google Scholar] [CrossRef]
  12. Plaza, A.; Du, Q.; Chang, Y.L. High Performance Computing for Hyperspectral Image Analysis: Perspective and State-of-the-art. In Proceedings of the 2009 IEEE International Geoscience and Remote Sensing Symposium, Cape Town, South Africa, 12–17 July 2009; pp. V-72–V-75. [Google Scholar] [CrossRef]
  13. Setoain, J.; Tenllado, C.; Prieto, M.; Valencia, D.; Plaza, A.; Plaza, J. Parallel Hyperspectral Image Processing on Commodity Graphics Hardware. In Proceedings of the 2006 International Conference on Parallel Processing Workshops (ICPPW’06), Columbus, OH, USA, 14–18 August 2006. [Google Scholar] [CrossRef]
  14. Setoain, J.; Prieto, M.; Tenllado, C.; Tirado, F. GPU for Parallel On-Board Hyperspectral Image Processing. Int. J. High Perform. Comput. Appl. 2008. [Google Scholar] [CrossRef]
  15. Plaza, A.; Valencia, D.; Plaza, J.; Martinez, P. Commodity cluster-based parallel processing of hyperspectral imagery. J. Parallel Distrib. Comput. 2006, 66, 345–358. [Google Scholar] [CrossRef]
  16. Plaza, J.; Pérez, R.; Plaza, A.; Martínez, P.; Valencia, D. Parallel Morphological/Neural Classification of Remote Sensing Images Using Fully Heterogeneous and Homogeneous Commodity Clusters. In Proceedings of the 2006 IEEE International Conference on Cluster Computing, Barcelona, Spain, 25–28 September 2006; pp. 1–10. [Google Scholar] [CrossRef] [Green Version]
  17. Aloisio, G.; Cafaro, M. A dynamic earth observation system. Parallel Comput. 2003, 29, 1357–1362. [Google Scholar] [CrossRef]
  18. Gorgan, D.; Bacu, V.; Stefanut, T.; Rodila, D.; Mihon, D. Grid based Satellite Image Processing Platform for Earth Observation Application Development. In Proceedings of the 2009 IEEE International Workshop on Intelligent Data Acquisition and Advanced Computing Systems: Technology and Applications, Rende, Italy, 21–23 September 2009; Volume 21, pp. 247–252. [Google Scholar] [CrossRef]
  19. Chen, Z.; Chen, N.; Yang, C.; Di, L. Cloud computing enabled web processing service for earth observation data processing. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2012, 5, 1637–1649. [Google Scholar] [CrossRef]
  20. Wu, Z.; Li, Y.; Plaza, A.; Li, J.; Xiao, F.; Wei, Z. Parallel and Distributed Dimensionality Reduction of Hyperspectral Data on Cloud Computing Architectures. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2016, 9, 2270–2278. [Google Scholar] [CrossRef]
  21. Haut, J.; Paoletti, M.; Plaza, J.; Plaza, A. Cloud implementation of the K-means algorithm for hyperspectral image analysis. J. Supercomput. 2017, 73. [Google Scholar] [CrossRef]
  22. Haut, J.M.; Gallardo, J.A.; Paoletti, M.E.; Cavallaro, G.; Plaza, J.; Plaza, A.; Riedel, M. Cloud deep networks for hyperspectral image analysis. IEEE Trans. Geosci. Remote Sens. 2019, 57, 9832–9848. [Google Scholar] [CrossRef]
  23. Quirita, V.A.A.; Da Costa, G.A.O.P.; Happ, P.N.; Feitosa, R.Q.; Da Silva Ferreira, R.; Oliveira, D.A.B.; Plaza, A. A New Cloud Computing Architecture for the Classification of Remote Sensing Data. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2017, 14, 409–416. [Google Scholar] [CrossRef]
  24. Zheng, X.; Xue, Y.; Guang, J.; Liu, J. Remote sensing data processing acceleration based on multi-core processors. In Proceedings of the 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Beijing, China, 10–15 July 2016; pp. 641–644. [Google Scholar] [CrossRef]
  25. Paoletti, M.E.; Haut, J.M.; Plaza, J.; Plaza, A.; Liu, Q.; Hang, R. Multicore Implementation of the Multi-Scale Adaptive Deep Pyramid Matching Model for Remotely Sensed Image Classification. In Proceedings of the 2017 IEEE International Geoscience and Remote Sensing Symposium, Fort Worth, TX, USA, 23–28 July 2017; pp. 2247–2250. [Google Scholar]
  26. Sevilla, J.; Bernabe, S.; Plaza, A. Unmixing-based content retrieval system for remotely sensed hyperspectral imagery on GPUs. J. Supercomput. 2014, 70, 588–599. [Google Scholar] [CrossRef]
  27. Wu, Z.; Wang, Q.; Plaza, A.; Li, J.; Wei, Z. Real-Time Implementation of the Sparse Multinomial Logistic Regression for Hyperspectral Image Classification on GPUs. IEEE Geosci. Remote Sens. Lett. 2015, 12, 1456–1460. [Google Scholar] [CrossRef]
  28. Paoletti, M.E.; Haut, J.M.; Plaza, J.; Plaza, A. Scalable recurrent neural network for hyperspectral image classification. J. Supercomput. 2020, 1–17. [Google Scholar] [CrossRef]
  29. Jaramago, J.A.G.; Paoletti, M.E.; Haut, J.M.; Fernandez-Beltran, R.; Plaza, A.; Plaza, J. GPU Parallel Implementation of Dual-Depth Sparse Probabilistic Latent Semantic Analysis for Hyperspectral Unmixing. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2019, 12, 3156–3167. [Google Scholar] [CrossRef]
  30. Leeser, M.; Belanovic, P.; Estlick, M.; Gokhale, M.; Szymanski, J.J.; Theiler, J. Applying Reconfigurable Hardware to the Analysis of Multispectral and Hyperspectral Imagery. In Proceedings of the International Symposium on Optical Science and Technology, San Diego, CA, USA, 17 January 2002; pp. 100–107. [Google Scholar]
  31. Williams, J.A.; Dawood, A.S.; Visser, S.J. FPGA-based cloud detection for real-time onboard remote sensing. In Proceedings of the 2002 IEEE International Conference on FieId-Programmable Technology (FPT 2002), Hong Kong, China, 16–18 December 2002; pp. 110–116. [Google Scholar] [CrossRef]
  32. Nie, Z.; Zhang, X.; Yang, Z. An FPGA Implementation of Multi-Class Support Vector Machine Classifier Based on Posterior Probability. Int. Proc. Comput. Sci. Inf. Technol. 2012. [Google Scholar] [CrossRef]
  33. Gonzalez, C.; Sánchez, S.; Paz, A.; Resano, J.; Mozos, D.; Plaza, A. Use of FPGA or GPU-based architectures for remotely sensed hyperspectral image processing. Integration 2013, 46, 89–103. [Google Scholar] [CrossRef]
  34. González, C.; Bernabé, S.; Mozos, D.; Plaza, A. FPGA implementation of an algorithm for automatically detecting targets in remotely sensed hyperspectral images. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2016, 9, 4334–4343. [Google Scholar] [CrossRef]
  35. Torti, E.; Acquistapace, M.; Danese, G.; Leporati, F.; Plaza, A. Real-time identification of hyperspectral subspaces. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2014, 7, 2680–2687. [Google Scholar] [CrossRef]
  36. Haut, J.M.; Bernabé, S.; Paoletti, M.E.; Fernandez-Beltran, R.; Plaza, A.; Plaza, J. Low–high-power consumption architectures for deep-learning models applied to hyperspectral image classification. IEEE Geosci. Remote Sens. Lett. 2018, 16, 776–780. [Google Scholar] [CrossRef]
  37. Paz, A.; Plaza, A. Clusters versus GPUs for parallel target and anomaly detection in hyperspectral images. EURASIP J. Adv. Signal Process. 2010, 2010, 1–18. [Google Scholar] [CrossRef] [Green Version]
  38. Bernabe, S.; López, S.; Plaza, A.; Sarmiento, R. GPU implementation of an automatic target detection and classification algorithm for hyperspectral image analysis. IEEE Geosci. Remote Sens. Lett. 2012, 10, 221–225. [Google Scholar] [CrossRef]
  39. Santos, L.; Magli, E.; Vitulli, R.; López, J.F.; Sarmiento, R. Highly-parallel GPU architecture for lossy hyperspectral image compression. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2013, 6, 670–681. [Google Scholar] [CrossRef]
  40. Paoletti, M.; Haut, J.; Plaza, J.; Plaza, A. Deep learning classifiers for hyperspectral imaging: A review. ISPRS J. Photogramm. Remote Sens. 2019, 158, 279–317. [Google Scholar] [CrossRef]
  41. Agathos, A.; Li, J.; Petcu, D.; Plaza, A. Multi-GPU implementation of the minimum volume simplex analysis algorithm for hyperspectral unmixing. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2014, 7, 2281–2296. [Google Scholar] [CrossRef]
  42. Melgani, F.; Bruzzone, L. Classification of hyperspectral remote sensing images with support vector machines. IEEE Trans. Geosci. Remote Sens. 2004, 42, 1778–1790. [Google Scholar] [CrossRef] [Green Version]
  43. Cristianini, N.; Shawe-Taylor, J. An Introduction to Support Vector Machines and other Kernel-Based Learning Methods; Cambridge University Press: Cambridge, UK, 2000. [Google Scholar]
  44. Mountrakis, G.; Im, J.; Ogole, C. Support vector machines in remote sensing: A review. ISPRS J. Photogramm. Remote Sens. 2011, 66, 247–259. [Google Scholar] [CrossRef]
  45. Osuna, E.; Freund, R.; Girosi, F. An improved training algorithm for support vector machines. In Proceedings of the Neural Networks for Signal Processing VII. 1997 IEEE Signal Processing Society Workshop, Amelia Island, FL, USA, 24–26 September 1997; pp. 276–285. [Google Scholar]
  46. Platt, J. Fast Training of Support Vector Machines Using Sequential Minimal Optimization; Advances in Kernel Methods—Support Vector Learning; AJ, MIT Press: Cambridge, MA, USA, 1999; pp. 185–208. [Google Scholar]
  47. Fan, R.E.; Chen, P.H.; Lin, C.J. Working set selection using second order information for training support vector machines. J. Mach. Learn. Res. 2005, 6, 1889–1918. [Google Scholar]
  48. Tan, K.; Zhang, J.; Du, Q.; Wang, X. GPU parallel implementation of support vector machines for hyperspectral image classification. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2015, 8, 4647–4656. [Google Scholar] [CrossRef]
  49. Li, Q.; Salman, R.; Kecman, V. An intelligent system for accelerating parallel SVM classification problems on large datasets using GPU. In Proceedings of the 2010 10th International Conference on Intelligent Systems Design and Applications, Cairo, Egypt, 29 November–1 December 2010; pp. 1131–1135. [Google Scholar]
  50. Wen, Z.; Shi, J.; Li, Q.; He, B.; Chen, J. ThunderSVM: A fast SVM library on GPUs and CPUs. J. Mach. Learn. Res. 2018, 19, 797–801. [Google Scholar]
  51. Cortes, C.; Vapnik, V. Support-vector networks. Mach. Learn. 1995, 20, 273–297. [Google Scholar] [CrossRef]
  52. Gale, D.; Kuhn, H.W.; Tucker, A.W. Linear programming and the theory of games. Act. Anal. Prod. Alloc. 1951, 13, 317–335. [Google Scholar]
  53. Platt, J. Sequential Minimal Optimization: A Fast Algorithm for Training Support Vector Machines; Technical Report MSR-TR-98-14; Microsoft Research: Redmond, WA, USA, 1998. [Google Scholar]
  54. Aiserman, M.; Braverman, E.M.; Rozonoer, L. Theoretical foundations of the potential function method in pattern recognition. Avtomat. I Telemeh 1964, 25, 917–936. [Google Scholar]
  55. Moreno, P.J.; Ho, P.P.; Vasconcelos, N. A Kullback-Leibler divergence based kernel for SVM classification in multimedia applications. In Proceedings of the Advances in Neural Information Processing Systems; 2004; pp. 1385–1392. Available online: https://dl.acm.org/doi/10.5555/2981345.2981516 (accessed on 16 April 2020).
  56. Mercer, J. Functions ofpositive and negativetypeand theircommection with the theory ofintegral equations. Philos. Trinsdictions Rogyal Soc. 1909, 209, 4–415. [Google Scholar]
  57. Hsu, C.W.; Lin, C.J. A comparison of methods for multiclass support vector machines. IEEE Trans. Neural Netw. 2002, 13, 415–425. [Google Scholar]
  58. Knerr, S.; Personnaz, L.; Dreyfus, G. Single-layer learning revisited: A stepwise procedure for building and training a neural network. In Neurocomputing; Springer: Berlin/Heidelberg, Germany, 1990; pp. 41–50. [Google Scholar]
  59. Platt, J.C.; Cristianini, N.; Shawe-Taylor, J. Large margin DAGs for multiclass classification. In Proceedings of the Advances in Neural Information Processing Systems; 2000; pp. 547–553. Available online: https://papers.nips.cc/paper/1773-large-margin-dags-for-multiclass-classification (accessed on 16 April 2020).
  60. Duan, K.B.; Rajapakse, J.C.; Nguyen, M.N. One-versus-one and one-versus-all multiclass SVM-RFE for gene selection in cancer classification. In Proceedings of the European Conference on Evolutionary Computation, Machine Learning and Data Mining in Bioinformatics, Valencia, Spain, 11–13 April 2007; pp. 47–56. [Google Scholar]
  61. Boser, B.E.; Guyon, I.M.; Vapnik, V.N. A training algorithm for optimal margin classifiers. In Proceedings of the Fifth Annual Workshop on Computational Learning Theory; 1992; pp. 144–152. Available online: https://dl.acm.org/doi/10.1145/130385.130401 (accessed on 16 April 2020).
  62. Chandra, R.; Dagum, L.; Kohr, D.; Menon, R.; Maydan, D.; McDonald, J. Parallel Programming in OpenMP; Morgan Kaufmann: San Mateo, CA, USA, 2001. [Google Scholar]
  63. Wu, Z.; Liu, J.; Plaza, A.; Li, J.; Wei, Z. GPU implementation of composite kernels for hyperspectral image classification. IEEE Geosci. Remote Sens. Lett. 2015, 12, 1973–1977. [Google Scholar]
  64. Camps-Valls, G.; Gomez-Chova, L.; Muñoz-Marí, J.; Vila-Francés, J.; Calpe-Maravilla, J. Composite kernels for hyperspectral image classification. IEEE Geosci. Remote Sens. Lett. 2006, 3, 93–97. [Google Scholar] [CrossRef]
  65. Osuna, E.; Freund, R.; Girosit, F. Training support vector machines: An application to face detection. In Proceedings of the IEEE Computer Society Conference on Computer Vision and Pattern Recognition, San Juan, PR, USA, 17–19 June 1997; pp. 130–136. [Google Scholar]
  66. Joachims, T. Making Large-Scale SVM Learning Practical; Technical Report; MIT Press: Cambridge, MA, USA, 1998. [Google Scholar]
  67. Keerthi, S.S.; Shevade, S.K.; Bhattacharyya, C.; Murthy, K.R.K. Improvements to Platt’s SMO algorithm for SVM classifier design. Neural Comput. 2001, 13, 637–649. [Google Scholar] [CrossRef]
  68. Hsu, C.W.; Lin, C.J. A simple decomposition method for support vector machines. Mach. Learn. 2002, 46, 291–314. [Google Scholar] [CrossRef] [Green Version]
  69. Palagi, L.; Sciandrone, M. On the convergence of a modified version of SVM light algorithm. Optim. Methods Softw. 2005, 20, 317–334. [Google Scholar] [CrossRef] [Green Version]
  70. Glasmachers, T.; Igel, C. Maximum-gain working set selection for SVMs. J. Mach. Learn. Res. 2006, 7, 1437–1466. [Google Scholar]
  71. Bo, L.; Jiao, L.; Wang, L. Working set selection using functional gain for LS-SVM. IEEE Trans. Neural Netw. 2007, 18, 1541–1544. [Google Scholar] [CrossRef]
  72. Glasmachers, T.; Igel, C. Second-order SMO improves SVM online and active learning. Neural Comput. 2008, 20, 374–382. [Google Scholar] [CrossRef]
  73. Dogan, U.; Glasmachers, T.; Igel, C. Fast Training of Multi-Class Support Vector Machines; Faculty of Science, University of Copenhagen: Copenhagen, Denmark, 2011. [Google Scholar]
  74. Tuma, M.; Igel, C. Improved working set selection for LaRank. In Proceedings of the International Conference on Computer Analysis of Images and Patterns, Seville, Spain, 29–31 August 2011; pp. 327–334. [Google Scholar]
  75. Wu, H.C. The Karush–Kuhn–Tucker optimality conditions in an optimization problem with interval-valued objective function. Eur. J. Oper. Res. 2007, 176, 46–59. [Google Scholar] [CrossRef]
  76. Merrill, D. CUB v1. 5.3: CUDA Unbound, a library of warp-wide, blockwide, and device-wide GPU parallel primitives. NVIDIA Res. 2015. [Google Scholar]
  77. Nvidia, C. Cublas library. NVIDIA Corp. Santa Clara Calif. 2008, 15, 31. [Google Scholar]
  78. Naumov, M.; Chien, L.; Vandermersch, P.; Kapasi, U. Cusparse library. In Proceedings of the GPU Technology Conference, San Jose, CA, USA, 20–23 September 2010. [Google Scholar]
  79. Vane, G.; Green, R.O.; Chrien, T.G.; Enmark, H.T.; Hansen, E.G.; Porter, W.M. The airborne visible/infrared imaging spectrometer (AVIRIS). Remote Sens. Environ. 1993, 44, 127–143. [Google Scholar] [CrossRef]
  80. Kunkel, B.; Blechinger, F.; Lutz, R.; Doerffer, R.; van der Piepen, H.; Schroder, M. ROSIS (Reflective Optics System Imaging Spectrometer)—A candidate instrument for polar platform missions. In Proceedings of the SPIE 0868 Optoelectronic technologies for remote sensing from space, Cannes, France, 17–20 November 1988; p. 8. [Google Scholar] [CrossRef]
  81. Xu, X.; Lil, f.; Plaza, A. Fusion of hyperspectral and LiDAR data using morphological component analysis. In Proceedings of the 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Beijing, China, 10–15 July 2016; pp. 3575–3578. [Google Scholar] [CrossRef]
  82. Babey, S.; Anger, C. A compact airborne spectrographic imager (CASI). In Quantitative Remote Sensing: An Economic Tool for the Nineties, Proceedings of the 12th IGARSS ’89 and Canadian Symposium on Remote Sensing, Vancouver, CO, Canada, 10–14 July 1989; Institute of Electrical and Electronics Engineers: New York, NY, USA, 1989; Volume 1, pp. 1028–1031. [Google Scholar]
  83. Chang, C.C.; Lin, C.J. LIBSVM: A Library for Support Vector Machines. ACM Trans. Intell. Syst. Technol. 2011, 2. [Google Scholar] [CrossRef]
Figure 1. General SVM scheme. The final goal is to approximate parameters w and b with the aim of defining the optimal hyperplane that maximizes the margin between the two classes.
Figure 1. General SVM scheme. The final goal is to approximate parameters w and b with the aim of defining the optimal hyperplane that maximizes the margin between the two classes.
Remotesensing 12 01257 g001
Figure 2. In order to deal with non-linearly separable data, the SVM applies a transformation Φ to translate the data from the original non-linear representation in space N n b , to a linear representation in space H , i.e., Φ : N n b H .
Figure 2. In order to deal with non-linearly separable data, the SVM applies a transformation Φ to translate the data from the original non-linear representation in space N n b , to a linear representation in space H , i.e., Φ : N n b H .
Remotesensing 12 01257 g002
Figure 3. Programming and memory model provided by NVIDIA CUDA for a GPU device.
Figure 3. Programming and memory model provided by NVIDIA CUDA for a GPU device.
Remotesensing 12 01257 g003
Figure 4. The inequality constraint C α i 0 forces the Lagrange multipliers to be into a box while the linear equality constraint i = 1 n t y i α i = 0 forces them to lie on a diagonal line. Therefore, one step of the SMO algorithm should find the optimum values α ^ t and α ^ l on a diagonal line segment.
Figure 4. The inequality constraint C α i 0 forces the Lagrange multipliers to be into a box while the linear equality constraint i = 1 n t y i α i = 0 forces them to lie on a diagonal line. Therefore, one step of the SMO algorithm should find the optimum values α ^ t and α ^ l on a diagonal line segment.
Remotesensing 12 01257 g004
Figure 5. Number of available samples in the considered HSI datasets.
Figure 5. Number of available samples in the considered HSI datasets.
Remotesensing 12 01257 g005
Figure 6. Classification maps for the Indian Pines dataset with the fixed training test sets in http://dase.grss-ieee.org. (a) False color composition. (b) Ground-truth (available labeled samples). (c) Fixed training set. (d) Fixed test set. (e) Classification map obtained by the RF classifier. (f) Classification map obtained by the MLR classifier. (g) Classification map obtained by the SVM classifier (LibSVM implementation). (h) Classification map obtained by the proposed GPU implementation. The corresponding overall classification accuracies are shown in brackets.
Figure 6. Classification maps for the Indian Pines dataset with the fixed training test sets in http://dase.grss-ieee.org. (a) False color composition. (b) Ground-truth (available labeled samples). (c) Fixed training set. (d) Fixed test set. (e) Classification map obtained by the RF classifier. (f) Classification map obtained by the MLR classifier. (g) Classification map obtained by the SVM classifier (LibSVM implementation). (h) Classification map obtained by the proposed GPU implementation. The corresponding overall classification accuracies are shown in brackets.
Remotesensing 12 01257 g006
Figure 7. Classification maps for the Pavia University dataset with the fixed training test sets in http://dase.grss-ieee.org. (a) False color composition. (b) Ground-truth (available labeled samples). (c) Fixed training set. (d) Fixed test set. (e) Classification map obtained by the RF classifier. (f) Classification map obtained by the MLR classifier. (g) Classification map obtained by the SVM classifier (LibSVM implementation). (h) Classification map obtained by the proposed GPU implementation. The corresponding overall classification accuracies are shown in brackets.
Figure 7. Classification maps for the Pavia University dataset with the fixed training test sets in http://dase.grss-ieee.org. (a) False color composition. (b) Ground-truth (available labeled samples). (c) Fixed training set. (d) Fixed test set. (e) Classification map obtained by the RF classifier. (f) Classification map obtained by the MLR classifier. (g) Classification map obtained by the SVM classifier (LibSVM implementation). (h) Classification map obtained by the proposed GPU implementation. The corresponding overall classification accuracies are shown in brackets.
Remotesensing 12 01257 g007
Figure 8. Classification maps for the University of Houston dataset with the fixed training test sets in http://dase.grss-ieee.org. (a) False color composition. (b) Ground-truth (available labeled samples). (c) Fixed training set. (d) Fixed test set. (e) Classification map obtained by the RF classifier. (f) Classification map obtained by the MLR classifier. (g) Classification map obtained by the SVM classifier (LibSVM implementation). (h) Classification map obtained by the proposed GPU implementation. The corresponding overall classification accuracies are shown in brackets.
Figure 8. Classification maps for the University of Houston dataset with the fixed training test sets in http://dase.grss-ieee.org. (a) False color composition. (b) Ground-truth (available labeled samples). (c) Fixed training set. (d) Fixed test set. (e) Classification map obtained by the RF classifier. (f) Classification map obtained by the MLR classifier. (g) Classification map obtained by the SVM classifier (LibSVM implementation). (h) Classification map obtained by the proposed GPU implementation. The corresponding overall classification accuracies are shown in brackets.
Remotesensing 12 01257 g008
Figure 9. Processing times (and speedups) measured in the considered CPUs (CPU0 in the platform 1 and CPU1 in the platform 2) and GPUs (GPU0 in the platform 1 and GPU1 in the platform 2) as the percentage of training samples increases.
Figure 9. Processing times (and speedups) measured in the considered CPUs (CPU0 in the platform 1 and CPU1 in the platform 2) and GPUs (GPU0 in the platform 1 and GPU1 in the platform 2) as the percentage of training samples increases.
Remotesensing 12 01257 g009
Figure 10. Memory-transfer times required by the proposed GPU implementation for the Pavia Center and Big Indian Pines datasets. In all cases, the training percentage was fixed to 40% of the available samples.
Figure 10. Memory-transfer times required by the proposed GPU implementation for the Pavia Center and Big Indian Pines datasets. In all cases, the training percentage was fixed to 40% of the available samples.
Remotesensing 12 01257 g010
Table 1. Most common kernels in the SVM method, where x i and x j are two samples of dataset X .
Table 1. Most common kernels in the SVM method, where x i and x j are two samples of dataset X .
Linear x i , x j
Polynomial ( γ x i , x j + c ) n
Gaussian exp ( γ | | x i x j | | 2 ) γ
Sigmoid tanh ( γ x i , x j + c ) n
Radial basis function (RBF) exp | | x i x j | | 2 2 σ 2
Table 2. Classification results obtained by different techniques for the Indian Pines, University of Pavia and University of Houston scenes, using the fixed training and test sets available for these datasets in http://dase.grss-ieee.org.
Table 2. Classification results obtained by different techniques for the Indian Pines, University of Pavia and University of Houston scenes, using the fixed training and test sets available for these datasets in http://dase.grss-ieee.org.
Indian PinesUniversity of PaviaHouston University
ClassRFMLR RFMLR RFMLR
SVM SVM SVM
[40][40]CPUGPU[40][40]CPUGPU[40][40]CPUGPU
132.8068.088.088.079.5977.783.0881.8582.5582.2482.3482.15
251.5678.0780.7479.3255.258.7867.4567.283.582.583.3683.36
344.4159.4167.5767.4345.4267.2264.8565.097.9499.899.899.8
426.4625.2551.5249.998.7374.2998.2898.3591.4698.398.9697.86
579.3488.3287.2386.2899.1498.8899.2899.396.6997.4498.7798.6
695.7196.8996.6196.6778.7793.5292.0692.399.1694.4197.996.78
720.050.0100.0100.080.4185.1288.8988.9575.3573.4177.4375.95
8100.099.298.898.890.9687.5892.1292.3133.0363.8260.357.78
916.040.060.050.097.6999.2296.3596.669.270.2576.7777.88
108.4756.1481.7182.03 43.955.661.2961.56
1189.6381.6486.9587.77 69.7974.1980.5581.2
1226.668.4477.6677.94 54.1270.4179.9280.65
1389.2596.2593.7594.25 59.8667.7270.8872.56
1492.089.9590.6490.83 99.3598.79100.0100.0
1538.7982.8376.7781.21 97.4295.5696.4197.42
1693.6493.1888.6488.64
OA65.6978.1684.2184.2570.1572.2378.9178.6772.9778.9881.8681.69
AA56.5473.3582.9182.4480.6682.4886.9386.8776.8981.6384.3184.24
K(x100)59.8774.9981.9882.063.0165.4573.3773.0970.9677.3180.4380.25
Table 3. Processing times (ms) obtained by our GPU implementation of the SVM algorithm on different platforms, using the datasets: Pavia Center and Big Indian Pines. The times indicate the transfers from host to device (HtoD), device to host (DtoH), device to device (DtoD) and the total processing times consumed by the kernels.
Table 3. Processing times (ms) obtained by our GPU implementation of the SVM algorithm on different platforms, using the datasets: Pavia Center and Big Indian Pines. The times indicate the transfers from host to device (HtoD), device to host (DtoH), device to device (DtoD) and the total processing times consumed by the kernels.
PAVIA CENTER
NVIDIA GeForce RTX 2070 (Laptop)NVIDIA GeForce GTX 1080Ti (Desktop)
Tr PercentHtoDDtoHDtoDKernelsHtoDDtoHDtoDKernels
11.9008001.0084500.40128089.6835102.0655400.9096630.40182453.969180
55.1964101.2119100.389169342.752095.9283301.3088900.445162192.82042
109.6503601.7107200.404029588.1151411.103641.7984700.448121292.20041
2018.907922.6587500.405595845.3827322.659272.7310100.465384503.46422
4037.788424.7610700.4447971369.810045.996254.8565800.505738816.56403
6056.839776.9371000.4837091924.370071.392926.9302100.5069871180.4300
8077.182979.0775300.4763792425.940095.875009.0327000.5476601448.0500
BIG INDIAN PINES
NVIDIA GeForce RTX 2070 (Laptop)NVIDIA GeForce GTX 1080Ti (Desktop)
Tr PercentHtoDDtoHDtoDKernelsHtoDDtoHDtoDKernels
172.27061044.04178020.125482167.880086.87498045.89442023.893302445.0100
5189.4537953.01270021.6584515000.320215.9661856.35232025.6550210679.120
10345.0400068.47627023.6376034717.630395.1083471.71834027.3472922486.610
20683.13051118.7884429.1764177072.650813.34057123.3149432.4095346963.590
401439.0210296.1476643.88958184683.751703.9300301.1484046.36005112824.33
602325.2100579.4322560.25476335392.102868.1800581.1937461.65440255449.66
803296.8700969.3961068.03797519363.014077.2100972.8647477.46072307944.21

Share and Cite

MDPI and ACS Style

Paoletti, M.E.; Haut, J.M.; Tao, X.; Miguel, J.P.; Plaza, A. A New GPU Implementation of Support Vector Machines for Fast Hyperspectral Image Classification. Remote Sens. 2020, 12, 1257. https://doi.org/10.3390/rs12081257

AMA Style

Paoletti ME, Haut JM, Tao X, Miguel JP, Plaza A. A New GPU Implementation of Support Vector Machines for Fast Hyperspectral Image Classification. Remote Sensing. 2020; 12(8):1257. https://doi.org/10.3390/rs12081257

Chicago/Turabian Style

Paoletti, Mercedes E., Juan M. Haut, Xuanwen Tao, Javier Plaza Miguel, and Antonio Plaza. 2020. "A New GPU Implementation of Support Vector Machines for Fast Hyperspectral Image Classification" Remote Sensing 12, no. 8: 1257. https://doi.org/10.3390/rs12081257

APA Style

Paoletti, M. E., Haut, J. M., Tao, X., Miguel, J. P., & Plaza, A. (2020). A New GPU Implementation of Support Vector Machines for Fast Hyperspectral Image Classification. Remote Sensing, 12(8), 1257. https://doi.org/10.3390/rs12081257

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