# Hyperspectral Image Classification Model Using Squeeze and Excitation Network with Deep Learning

Authors: Rajendran T; Prajoona Valsalan; Amutharaj J; Jenifer M; Rinesh S; Charlyn Pushpa Latha G; Anitha T
Journal: Computational Intelligence and Neuroscience
 (2022)

---

## Abstract

In the domain of remote sensing, the classification of hyperspectral image (HSI) has become a popular topic. In general, the complicated features of hyperspectral data cause the precise classification difficult for standard machine learning approaches. Deep learning-based HSI classification has lately received a lot of interest in the field of remote sensing and has shown promising results. As opposed to conventional hand-crafted feature-based classification approaches, deep learning can automatically learn complicated features of HSIs with a greater number of hierarchical layers. Because HSI’s data structure is complicated, applying deep learning to it is difficult. The primary objective of this research is to propose a deep feature extraction model for HSI classification. Deep networks can extricate features of spatial and spectral from HSI data simultaneously, which is advantageous for increasing the performances of the proposed system. The squeeze and excitation (SE) network is combined with convolutional neural networks (SE-CNN) in this work to increase its performance in extracting features and classifying HSI. The squeeze and excitation block is designed to improve the representation quality of a CNN. Three benchmark datasets are utilized in the experiment to evaluate the proposed model: Pavia Centre, Pavia University, and Salinas. The proposed model’s performance is validated by a performance comparison with current deep transfer learning approaches such as VGG-16, Inception-v3, and ResNet-50. In terms of accuracy on each class of datasets and overall accuracy, the proposed SE-CNN model outperforms the compared models. The proposed model achieved an overall accuracy of 96.05% for Pavia University, 98.94% for Pavia Centre dataset, and 96.33% for Salinas dataset.


---

## Body

## 1. Introduction

With the advancements of remote sensing technology, the use of hyperspectral imaging is becoming increasingly common. The precise classifications of ground features using HSI is a significant research topic that has received a lot of interest. Because of its high resolving powers for good spectra, hyperspectral images have a wide variety of uses in the environmental, medical, defense, and mining. The collection of HSI is dependent on imaging spectrometer deployed in various locations. In the 1980s, the imaging spectrum was created. It was utilized to image electromagnetic waves in the ultraviolet, near-infrared, visible, and midinfrared ranges. The imaging spectrometer can photograph in a variety of continuous and extremely narrow band, allowing every pixel to obtain a completely emitted or reflected spectrum in the wavelength range of interest. As a result, hyperspectral images feature great spectral resolution, multiple bands, and a lot of information. Image corrections, noise reductions, transformations, dimensionality reductions, and classification are the most common methods for processing hyperspectral remote sensing images. As opposed to conventional images, hyperspectral images include a wealth of spectral data. These spectral data could represent the physical structures and chemical structures of the target objects, which supports in classification of image [1].The purpose of computer remote sensing image classification was to detect and classify data from the environment and surface of Earth on remote sensing image in order for determining the feature data matching to the image data and extricate the valuable feature data. The specific use of automated pattern identification technology in the remote sensing domain is computer classification of remote sensing images. The count of imaging bands in HSI was more than in multispectral image, and the capacity to resolve objects was solid, that is, the higher the spectral resolutions. Thus, because of the higher-dimensional features of HSI, as well as the similarity among the mixed pixels and the spectra, the HSI classification methodology still confronts a number of difficulties; the most significant of those are listed below.(i)
Hyperspectral image data have a higher dimensionality. Because HSIs are produced by combining values of spectral reflectance acquired by space-borne or airborne imaging spectrometer in hundreds of bands, the related spectrum data dimension of HSIs can likewise be hundreds of dimensions.(ii)
Variability of spectral information in space. The HSIs are affected by atmospheric conditions, structures, and distributions of ground features, sensors, and the surrounding environments.(iii)
The interference of noise and background variables during the capture of hyperspectral images has a significant impact on the quality of the obtained data. The classification accuracy of hyperspectral images is directly affected by the quality of image [1].Deep neural network’s work is based on extracting features from raw input data via processing layer-by-layer of input data. Although the use of spectral characteristics in isolation can be beneficial, it might not be sufficient in many instances. When two diverse items have the similar spectral signatures, their forms and textures can be used to distinguish them. As a result, spectral and spatial data could be merged to enhance HSI classification. Deep networks could extricate spatial and spectral characteristics from HSI data at the same time. However, traditional deep neural networks have significant limitations. They need a large number of training samples and a significant amount of work to tune hyperparameters [2].Deep learning networks are classified into many types. The convolutional neural networks is one of the most well-known networks for extracting and classifying HSI features [3]. The squeeze and excitation (SE) network was combined with CNN (SE-CNN) in this work to increase its performance in extracting features and classifying HSI images. The squeeze and excitation block is designed to improve the representation quality of a convolutional neural network. In the experiment, three HSI benchmark datasets were utilized for evaluating the proposed model. The proposed model is verified based on the obtained performance by comparing it to the existing approaches used for HSI classification.


## 2. Related Works

Bandar and Munif proposed the deep hypernetwork architecture capable of learning the deep features of HSIs and provided fine performances without the need for dataset augmentation or extensive preprocessing. The primary aim of this HSI classification algorithm was to anticipate the type of land cover by allocating and labelling separate pixels with multiple frequency bands into individual classes. By integrating and merging the basic principles of both Inception and ResNet into a single model, a novel CNN model was presented that used both the deep residual network and inception models. The performance was not improved over conventional techniques because of the minimal volume of labelled training sample [4]. Madhumitha et al. proposed the HSI classification model using an edge detection, optimization, dimensionality reduction, and classification procedures. The image denoising with filter functionality and edge identification method were used to implement the edge detection procedure. The optimization procedure was carried out by first producing the initial populations, then computing fitness, Cauchy mutation, crossover, and lastly determining the decreasing rate with a value of optimal threshold. The input HSI with L-bands was segmented utilizing auto encoder after the denoising was removed, and dimensionality reduction was performed, yielding the output images with M components. CNN was used to classify the HSI, and the M components were categorized using two stages of convolution and max-pooling. The dimensionality reduction technique in HSI improved overall performance while reducing computing complexity [5].Xin et al. presented the multisource deep transfer learning framework (MS-DTL). This multisource transfer learning approach was utilized to classify HSIs in order to get more information and ease the problem of limited samples. The ResNet-like basic model was created to address the needs of many sources and a single target. A loss that included the cross-entropy losses from each source HSI was used to train this model. The shallow layers were moved, the number of classes was raised by utilizing various sources to extract global features, and this method's performance was enhanced [6]. Yahya et al. used CNNs to automatically detect haploid and diploid maize seeds using transfer learning-based technique. For this work, AlexNet, VGGNet, GoogLeNet, and ResNet were used particularly. In contrast, the VGG-19 outperformed the others [7]. Fuding et al. proposed a spectral–spatial HSI classification technique based on superpixel pooling CNN and transfer learning (SP–CNN). The model was divided into three phases. The first section comprised of a convolution and pooling operation, which was a down-sampling method used to extract the key spectral features of an HSI. The next section used up-sampling and super-pixel pooling to investigate an HSI’s spatial structure information. Finally, the hyperspectral data were loaded into a fully connected neural network using each superpixel as a basic input rather than a pixel [8]. Yakoub et al. proposed a classification of remote sensing scene model based on CNN and other deep transfer learning algorithms. The detailed performance analysis was discussed in [9].For extracting meaningful features of spectral–spatial from HSI images, Yushi et al. presented the 3-D CNN-based features extraction approach with combined regularization. A regularized deep features extraction technique for HSI classification using CNN was presented. This model extracted deep features from HSIs using multiple convolutional and pooling layers that were discriminant, nonlinear, and invariant. These features were beneficial for classification of image and identification of target. This model has an overfitting problem [10].Shengjie et al. presented MDL4OW, a multitask deep learning technique for HSI classification with unknown classes, in which a multitask network was used to do classification and reconstruction at the same time. The classification supplied the likelihood of known classes, but the reconstruction estimated the unknown value. According to the findings, the integrity of a classification system was significant in HSI classification with unknown classes [11]. To solve the overfitting problem, Haokui et al. presented AINet, a 3D asymmetric inception network for HSI classification. Initially, AINet utilized a lightweight 3D CNN that was nevertheless highly deep, allowing deep learning to be used to extract representative features and ease the difficulty caused by the restricted annotation datasets. Second, while assessing the HSI’s property, spectral signatures were prioritized above geographical contexts. Furthermore, a data fusion transfer learning approach was used to improve model initialization and reduce training time [12].Xiangyong et al. presented an active deep learning technique for HSI classification that combined active and deep learning into a single framework. First, the CNN was trained with a small quantity of labelled pixels. Following that, the most informative pixels from the candidate pool were actively chosen for labelling. The CNN was then fine-tuned using the newly generated training set, which included the newly labelled pixels. Finally, for improving the classification performance, the Markov random field was used to enforce class label smoothness [13]. From the review of related works, it can be seen that most of the researches used the deep learning techniques for the classification of HSI. Deep learning proved as the best approach for the classification, especially CNN when it comes to image processing. Although, the CNN has some drawbacks of degraded performance due to the overfitting and underfitting issues. To solve these issues, the deep transfer learning techniques can be used, which is a pretrained model that supports to overcome these issues and improves the performance of the model in classification. Different DTL techniques have been proposed earlier, from that the squeeze and excitation network was a recently proposed model that is used in this proposed research with the CNN for the classification of HSI. The SE-Net has an improved performance comparing to the other DTL techniques like Inception, VGG, ResNet, ResNext, and so on.


## 3. Proposed Methodology

The HSI classification model developed in this work is based on a feature extraction method that combines a squeeze and excitation network with a deep convolutional neural network. Deep learning is difficult to apply to HSI because HSI’s data structure is complicated. In general, the neural network has a strong representation capability and a larger volume of training samples. The primary objective of this work was to create a deep feature extraction model for HSI classification. Deep networks are capable of extracting spatial and spectral characteristics from HSI data simultaneously, which is advantageous for increasing the performances of the presented model. The SE network was combined with CNN (SE-CNN) in this research to increase its performance in extracting features and classifying HSI.


### 3.1. Convolutional Neural Network

CNN is a feed forward neural network with multiple convolutions and pooling operations, which has the benefits of automated spatial data learning and overfitting issue management. It performs admirably in image recognition, tracking target, and natural language processing. In this work, the CNN method was used for extracting deep HSI classification features. The CNN training method consists mostly of forward propagations and backward parameter updates, and every network layer was described in detail in the following sections.


#### 3.1.1. Forward Propagation Process

CNN convolves input image to the C-layer by using varying widths of the convolution kernel. After applying the bias, the features of abstract textures of the input image were extricated using the activation function to accomplish feature improvement. Convolution could be represented as follows:(1)xjl=f∑i∈Mjxil−1∗wijl+bjl.Here,xil was jth factor of lth layer; Mj was jth convolution segment of l-1 layered feature maps; xil−1 was the factor in that; wijl was the lth layer’s weight matrix; bjl was the bias; f was commonly the function of nonlinear ReLU activations, and it was represented by,(2)fx=max0,x.


#### 3.1.2. Pooling Layer

The pooling layer, otherwise called as down-sampling layer, was often found beyond the C-layer. It could reduce unnecessary features through employing a down-sampling mechanism for avoiding additional over fitting and parameters of lower networks. Assuming thatl-1th layer was C-layer, the following lth P-layer may be written like,(3)xjl=fβjl∗downxjl−1+bjl,where βjl indicates the weight and bjl indicates the bias of the lth layer’s jth feature map, correspondingly. The down (∗) symbol represents the function of down-sampling, which usually includes mean, max, and stochastic pooling. To decrease the dimension of output feature maps, the max-pooling with shift-invariances was utilized.


#### 3.1.3. FC Layer

The last FC layer can acquire the output class or the input sample probability after numerous alternative convolution and pooling processes. All neuron nodes in thekth layer are linked to every output node in the upper k-1th layer in the FC-layer, and its mathematical expression is as follows:(4)yk=fwkxk−1+bk,where k denotes the network layer’s numerical order; yk was the FC-layer outcome; xk−1 was the unfolded single dimension eigen vector; wk was the weighted coefficient; bk indicates the bias; and f (∗), activation functions of the last layer in a FC-layer was Softmax functions for classifications operations. Furthermore, after every C-layer or FC-layer, a dropout function called as a regularization approach could be applied to improve the generalization capabilities of proposed CNN approach and minimize overfitting while training the approach [14].


#### 3.1.4. Back Propagation for Parameters Update

The input samples target could be obtained after forward propagation. The network parameters are then returned to their original state by reducing the loss functions of the actual and target outputs. HSI classification may be viewed as an image multiclassification process in this work. As a result, the categorical cross-entropy loss function was used, which can be expressed as:(5)E=1n∑k=1nyklntk+1−ykln1−tk,where n was the sample amount, yk was the actual target, and tk was the kth sample’s predictive value, accordingly. The gradient descent approach could minimize the loss function during model training. The CNN model’s adaptive parameters w and b are gradually updated by computing the partial derivatives of (5). The following are the expressions:(6)wijl=wijl−α∂E∂wijl,(7)bjl=bjl−α∂E∂bjl.In (6) and (7), α is the learning rate that controls the parameter updates stride. A suitable α may increase network convergence speed and prevent the network from declining into the local optimal solution. As a result, the time-based learning rate schedule used in this study is as follows:(8)α=α∗11+decay∗epoch,where decay represents the reduction in learning value from the earlier or prior epoch through the particular fixed values; and epoch was present training epoch. At last, the aforementioned back-propagation method completes the whole CNN’s training process [15].


### 3.2. Channel Domain Attention

To improve the model’s feature extraction performance and therefore to achieve precise classification, the channel attention in Figure1, SENet, will be used in this work for recalibrating deep featured maps produced by CNNs. SENet could learn global data automatically and apply varying weight ratios for filtering channels to choose target on key features and eliminate unnecessary data using the squeeze and excitation procedures. A series of convolutional transformations Ftr can convert the provided image X, as input with the dimension of (W′, H′, C′), it could be mapped to the feature map U while U ∈ RH×W×C. The output U = [u1, u2, …, uC] may be written as follows:(9)uc=vc∗X=∑s=1C′vCs∗Xs.where ∗ represents the convolution operation; V =  v1,v2…vC denotes the learnt convolution kernels; and vC  =  vc1,vC2…vCc is the Cth 2D spatial filter kernels; X = [x1, x2, ..., xC′] and uc∈RH×W. The squeeze transform Fsq then translates feature mappings U to the global spatial single dimension feature vectors, and a statistic z ∈ RC produced by compacting U with spatial dimension H×W may be expressed as follows:(10)zc=Fsquc=1H×W∑i=1H∑j=1Wuci,j.Figure 1
Squeeze and excitation block [16].According to (10), the above procedure was a global average pooling for obtaining the statistic z. Hence, using the self-gating method with dual FC-layers, the excitation operation was developed for performing evaluation of weight on all channels for adaptive featured recalibrations. It may be written as follows:(11)s=Fexz,W=σgz,W=σW2δW1,z.where δ was the ReLU activation function, W1∈RC/r×C and W2∈RC×C/r, and r represents a ratio of dimensionality reduction. At last, the SE block’s output was obtained by rescaling U with the activation s.(12)X˜c=Fscaleuc,sc=sc⋅uc,where X˜=X˜1,X˜2,…,X˜c and Fscaleuc,sc signify channel-wise scalar multiplication sc and feature map uc∈RH×W [17].The descriptive realization procedure of SE-CNN model was on the basis of theoretical analysis discussed above, where the deep featured maps were extricated initially by a CNN approach, and single dimension weighted vectors of filter channel were acquired by the group of consecutive performances comprising FC-layer, global pooling, ReLU, and Sigmoid activations in Figure2. Thus, for precise image identification, the feature maps generated by the CNN were recalibrated utilizing the aforementioned channel weight ratio. As a result, this research will integrate CNN and SENet for creating a SE-CNN model for HSI classification.Figure 2
Proposed SE-CNN model for HSI classification.To train the CNN model more easily, the image size was typically 2n, such as 32, 64, 128, 256, and so on. In addition, based on the efficiency of detection and HSI classification accuracy, this model resizes the image as 64 × 64 pixels in size. Following that, the SE-CNN model was created for extracting deep features in comparison to other CNN models. For input images, the SE-CNN model has several interchanging P-layers and C-layers, as well as one or few FC-layers.Each of the aforementioned CNN models has 32, 32, 64, or 128 filters in its C-layer, and strides of P-layers and C-layer performances were fixed to one. In addition, a zero-padding approach with the same padding was used on CNN models for keeping feature maps the same size. The SE module recalibrates the CNNs’ output feature maps, including squeeze, excitation, and scale procedures. Finally, the FC-layers flatten the calibrated feature maps, allowing a Softmax classifier to classify the images. Table1 represents the particular parameters of the proposed deep learning model, where C1∼C4 denotes the convolutional layers, P1∼P4 denotes the max-pooling layers, and FC-layer finally for HSI classification. Furthermore, the models described above all were applied with batch normalization following the initial C-layer and trained with a batch size of 32. Simultaneously, to prevent the network from overfitting, a dropout operation was presented after the FC-layer. Because the SE module’s dimensionality reduction ratio r impacts classification accuracy, relevant tests and analysis are carried out in this research.Table 1
Parameters of CNN model.
LayerParametersInputInput image with a size of (64 × 64, 3)C1Conv 2D (32, 3 × 3)BN1Batch normalizationP1Maxpooling 2D (2, 2)C2Conv 2D (32, 3 × 3)P2Maxpooling 2D (2, 2)C3Conv 2D (64, 3 × 3) drop out (0.35) ReLUP3Maxpooling 2D (2, 2)C4Conv 2D (128, 3 × 3)P4Maxpooling 2D (2, 2)FC11024FC2256OutputClassification of images


## 3.1. Convolutional Neural Network

CNN is a feed forward neural network with multiple convolutions and pooling operations, which has the benefits of automated spatial data learning and overfitting issue management. It performs admirably in image recognition, tracking target, and natural language processing. In this work, the CNN method was used for extracting deep HSI classification features. The CNN training method consists mostly of forward propagations and backward parameter updates, and every network layer was described in detail in the following sections.


### 3.1.1. Forward Propagation Process

CNN convolves input image to the C-layer by using varying widths of the convolution kernel. After applying the bias, the features of abstract textures of the input image were extricated using the activation function to accomplish feature improvement. Convolution could be represented as follows:(1)xjl=f∑i∈Mjxil−1∗wijl+bjl.Here,xil was jth factor of lth layer; Mj was jth convolution segment of l-1 layered feature maps; xil−1 was the factor in that; wijl was the lth layer’s weight matrix; bjl was the bias; f was commonly the function of nonlinear ReLU activations, and it was represented by,(2)fx=max0,x.


### 3.1.2. Pooling Layer

The pooling layer, otherwise called as down-sampling layer, was often found beyond the C-layer. It could reduce unnecessary features through employing a down-sampling mechanism for avoiding additional over fitting and parameters of lower networks. Assuming thatl-1th layer was C-layer, the following lth P-layer may be written like,(3)xjl=fβjl∗downxjl−1+bjl,where βjl indicates the weight and bjl indicates the bias of the lth layer’s jth feature map, correspondingly. The down (∗) symbol represents the function of down-sampling, which usually includes mean, max, and stochastic pooling. To decrease the dimension of output feature maps, the max-pooling with shift-invariances was utilized.


### 3.1.3. FC Layer

The last FC layer can acquire the output class or the input sample probability after numerous alternative convolution and pooling processes. All neuron nodes in thekth layer are linked to every output node in the upper k-1th layer in the FC-layer, and its mathematical expression is as follows:(4)yk=fwkxk−1+bk,where k denotes the network layer’s numerical order; yk was the FC-layer outcome; xk−1 was the unfolded single dimension eigen vector; wk was the weighted coefficient; bk indicates the bias; and f (∗), activation functions of the last layer in a FC-layer was Softmax functions for classifications operations. Furthermore, after every C-layer or FC-layer, a dropout function called as a regularization approach could be applied to improve the generalization capabilities of proposed CNN approach and minimize overfitting while training the approach [14].


### 3.1.4. Back Propagation for Parameters Update

The input samples target could be obtained after forward propagation. The network parameters are then returned to their original state by reducing the loss functions of the actual and target outputs. HSI classification may be viewed as an image multiclassification process in this work. As a result, the categorical cross-entropy loss function was used, which can be expressed as:(5)E=1n∑k=1nyklntk+1−ykln1−tk,where n was the sample amount, yk was the actual target, and tk was the kth sample’s predictive value, accordingly. The gradient descent approach could minimize the loss function during model training. The CNN model’s adaptive parameters w and b are gradually updated by computing the partial derivatives of (5). The following are the expressions:(6)wijl=wijl−α∂E∂wijl,(7)bjl=bjl−α∂E∂bjl.In (6) and (7), α is the learning rate that controls the parameter updates stride. A suitable α may increase network convergence speed and prevent the network from declining into the local optimal solution. As a result, the time-based learning rate schedule used in this study is as follows:(8)α=α∗11+decay∗epoch,where decay represents the reduction in learning value from the earlier or prior epoch through the particular fixed values; and epoch was present training epoch. At last, the aforementioned back-propagation method completes the whole CNN’s training process [15].


## 3.1.1. Forward Propagation Process

CNN convolves input image to the C-layer by using varying widths of the convolution kernel. After applying the bias, the features of abstract textures of the input image were extricated using the activation function to accomplish feature improvement. Convolution could be represented as follows:(1)xjl=f∑i∈Mjxil−1∗wijl+bjl.Here,xil was jth factor of lth layer; Mj was jth convolution segment of l-1 layered feature maps; xil−1 was the factor in that; wijl was the lth layer’s weight matrix; bjl was the bias; f was commonly the function of nonlinear ReLU activations, and it was represented by,(2)fx=max0,x.


## 3.1.2. Pooling Layer

The pooling layer, otherwise called as down-sampling layer, was often found beyond the C-layer. It could reduce unnecessary features through employing a down-sampling mechanism for avoiding additional over fitting and parameters of lower networks. Assuming thatl-1th layer was C-layer, the following lth P-layer may be written like,(3)xjl=fβjl∗downxjl−1+bjl,where βjl indicates the weight and bjl indicates the bias of the lth layer’s jth feature map, correspondingly. The down (∗) symbol represents the function of down-sampling, which usually includes mean, max, and stochastic pooling. To decrease the dimension of output feature maps, the max-pooling with shift-invariances was utilized.


## 3.1.3. FC Layer

The last FC layer can acquire the output class or the input sample probability after numerous alternative convolution and pooling processes. All neuron nodes in thekth layer are linked to every output node in the upper k-1th layer in the FC-layer, and its mathematical expression is as follows:(4)yk=fwkxk−1+bk,where k denotes the network layer’s numerical order; yk was the FC-layer outcome; xk−1 was the unfolded single dimension eigen vector; wk was the weighted coefficient; bk indicates the bias; and f (∗), activation functions of the last layer in a FC-layer was Softmax functions for classifications operations. Furthermore, after every C-layer or FC-layer, a dropout function called as a regularization approach could be applied to improve the generalization capabilities of proposed CNN approach and minimize overfitting while training the approach [14].


## 3.1.4. Back Propagation for Parameters Update

The input samples target could be obtained after forward propagation. The network parameters are then returned to their original state by reducing the loss functions of the actual and target outputs. HSI classification may be viewed as an image multiclassification process in this work. As a result, the categorical cross-entropy loss function was used, which can be expressed as:(5)E=1n∑k=1nyklntk+1−ykln1−tk,where n was the sample amount, yk was the actual target, and tk was the kth sample’s predictive value, accordingly. The gradient descent approach could minimize the loss function during model training. The CNN model’s adaptive parameters w and b are gradually updated by computing the partial derivatives of (5). The following are the expressions:(6)wijl=wijl−α∂E∂wijl,(7)bjl=bjl−α∂E∂bjl.In (6) and (7), α is the learning rate that controls the parameter updates stride. A suitable α may increase network convergence speed and prevent the network from declining into the local optimal solution. As a result, the time-based learning rate schedule used in this study is as follows:(8)α=α∗11+decay∗epoch,where decay represents the reduction in learning value from the earlier or prior epoch through the particular fixed values; and epoch was present training epoch. At last, the aforementioned back-propagation method completes the whole CNN’s training process [15].


## 3.2. Channel Domain Attention

To improve the model’s feature extraction performance and therefore to achieve precise classification, the channel attention in Figure1, SENet, will be used in this work for recalibrating deep featured maps produced by CNNs. SENet could learn global data automatically and apply varying weight ratios for filtering channels to choose target on key features and eliminate unnecessary data using the squeeze and excitation procedures. A series of convolutional transformations Ftr can convert the provided image X, as input with the dimension of (W′, H′, C′), it could be mapped to the feature map U while U ∈ RH×W×C. The output U = [u1, u2, …, uC] may be written as follows:(9)uc=vc∗X=∑s=1C′vCs∗Xs.where ∗ represents the convolution operation; V =  v1,v2…vC denotes the learnt convolution kernels; and vC  =  vc1,vC2…vCc is the Cth 2D spatial filter kernels; X = [x1, x2, ..., xC′] and uc∈RH×W. The squeeze transform Fsq then translates feature mappings U to the global spatial single dimension feature vectors, and a statistic z ∈ RC produced by compacting U with spatial dimension H×W may be expressed as follows:(10)zc=Fsquc=1H×W∑i=1H∑j=1Wuci,j.Figure 1
Squeeze and excitation block [16].According to (10), the above procedure was a global average pooling for obtaining the statistic z. Hence, using the self-gating method with dual FC-layers, the excitation operation was developed for performing evaluation of weight on all channels for adaptive featured recalibrations. It may be written as follows:(11)s=Fexz,W=σgz,W=σW2δW1,z.where δ was the ReLU activation function, W1∈RC/r×C and W2∈RC×C/r, and r represents a ratio of dimensionality reduction. At last, the SE block’s output was obtained by rescaling U with the activation s.(12)X˜c=Fscaleuc,sc=sc⋅uc,where X˜=X˜1,X˜2,…,X˜c and Fscaleuc,sc signify channel-wise scalar multiplication sc and feature map uc∈RH×W [17].The descriptive realization procedure of SE-CNN model was on the basis of theoretical analysis discussed above, where the deep featured maps were extricated initially by a CNN approach, and single dimension weighted vectors of filter channel were acquired by the group of consecutive performances comprising FC-layer, global pooling, ReLU, and Sigmoid activations in Figure2. Thus, for precise image identification, the feature maps generated by the CNN were recalibrated utilizing the aforementioned channel weight ratio. As a result, this research will integrate CNN and SENet for creating a SE-CNN model for HSI classification.Figure 2
Proposed SE-CNN model for HSI classification.To train the CNN model more easily, the image size was typically 2n, such as 32, 64, 128, 256, and so on. In addition, based on the efficiency of detection and HSI classification accuracy, this model resizes the image as 64 × 64 pixels in size. Following that, the SE-CNN model was created for extracting deep features in comparison to other CNN models. For input images, the SE-CNN model has several interchanging P-layers and C-layers, as well as one or few FC-layers.Each of the aforementioned CNN models has 32, 32, 64, or 128 filters in its C-layer, and strides of P-layers and C-layer performances were fixed to one. In addition, a zero-padding approach with the same padding was used on CNN models for keeping feature maps the same size. The SE module recalibrates the CNNs’ output feature maps, including squeeze, excitation, and scale procedures. Finally, the FC-layers flatten the calibrated feature maps, allowing a Softmax classifier to classify the images. Table1 represents the particular parameters of the proposed deep learning model, where C1∼C4 denotes the convolutional layers, P1∼P4 denotes the max-pooling layers, and FC-layer finally for HSI classification. Furthermore, the models described above all were applied with batch normalization following the initial C-layer and trained with a batch size of 32. Simultaneously, to prevent the network from overfitting, a dropout operation was presented after the FC-layer. Because the SE module’s dimensionality reduction ratio r impacts classification accuracy, relevant tests and analysis are carried out in this research.Table 1
Parameters of CNN model.
LayerParametersInputInput image with a size of (64 × 64, 3)C1Conv 2D (32, 3 × 3)BN1Batch normalizationP1Maxpooling 2D (2, 2)C2Conv 2D (32, 3 × 3)P2Maxpooling 2D (2, 2)C3Conv 2D (64, 3 × 3) drop out (0.35) ReLUP3Maxpooling 2D (2, 2)C4Conv 2D (128, 3 × 3)P4Maxpooling 2D (2, 2)FC11024FC2256OutputClassification of images


## 4. Experiment Analysis

This section focuses on and presents a detailed review of experiments based on HSI classification utilizing three HSI datasets. Two evaluation measures were utilized to assess the proposed method’s efficiency in distinguishing different classes of hyperspectral images. The accuracy of all classes in the testing set was determined as the first measure, and the overall accuracy was the second measure, which was computed as the total of properly classified hyperspectral samples split by the overall count of HSI samples in testing set. To demonstrate the validity of the proposed technique, different existing classification methods for hyperspectral imagery were used for the comparison, namely, VGG-16, Inception-v3, and ResNet-50. MATLAB 2017a was used to run the tests using a system with Intel core i5 CPU with 2.9 GHz speed, RAM with 8 GB, and a 64 bit Windows 10-OS [18–25].


### 4.1. Dataset Description

In this research, three HSI benchmark datasets such as Pavia Centre, Pavia University, and Salinas, were utilized for evaluating the proposed approach, and their descriptions are as follows.


#### 4.1.1. Pavia Centre and University Datasets

Pavia Centre and University are two images captured by the ROSIS sensor during the fly above Pavia, Italy. Pavia University has 103 spectral bands, where Pavia Centre has 102. Pavia Centre has the 1096 × 1096 pixels images, while Pavia University has the 610 × 610 pixels images although some samples in both images are insufficient in information and deleted before processing. The spatial resolution was 1.3 meters. Both image ground truths distinguish nine groups. Table2 represents the description of the Pavia Centre dataset, and Table 3 represents the description of the Pavia University dataset.Table 2
Pavia Centre dataset description.
ClassSampleWater824Trees820Asphalt816Self-blocking bricks808Bitumen808Tiles1260Shadows476Meadows824Bare soil820Total7456Table 3
Pavia University dataset description.
ClassSampleAsphalt6631Meadows18649Gravel2099Trees3064Painted metal sheets1345Bare soil5029Bitumen1330Self-blocking bricks3682Shadows947Total42776


#### 4.1.2. Salinas

This dataset was collected by the airborne visible/infrared imaging spectrometer (AVIRIS) sensor above the city of Salinas Valley, California, USA. This collection incorporates images with the size of 512 × 217 pixels, a spatial resolution of 3.7 meters, and 204 bands after eliminating 20 water absorption bands. The accessible ground references map includes 16 different classes, such as vegetables, vineyards, and bare soil as shown in Table4.Table 4
Salinas scene dataset description.
ClassSampleCorn_senesced_green_weeds3278Celery3579Brocoli_green_weeds_12009Fallow_smooth2678Fallow1976Brocoli_green_weeds_23726Fallow_rough_plow1394Soil_vineyard_develop6203Stubble3959Vineyard_untrained7268Grapes_untrained11271Lettuce_romaine_6 wk916Lettuce_romaine_4 wk1068Vineyard_vertical_trellis1807Lettuce_romaine_7 wk1070Lettuce_romaine_5 wk1927Total54129Figure3 represents the ground truth images of the datasets utilized in this research for evaluation. It is worth noting that the zero value in the ground truth maps displayed in Figures 3(a)–3(c) indicates clutter pixels, which could not be classified as a new class.Figure 3
Ground Truth images: (a) Pavia University; (b) Pavia Centre; (c) Salinas.
(a)(b)(c)


### 4.2. Performance Analysis

For the performance evaluation, the datasets are divided into 50% to train and 50% to test. As shown in the following tabulations, Table5 represents the accuracy of every class and overall accuracy of total classes based on the classification performance of the proposed model using the Pavia Centre dataset. This Pavia Centre dataset is divided into 3728 samples for training and 3728 samples for testing for the experiment. The classification performance of the proposed model was compared with other existing deep transfer learning approaches like Inception-v3, VGG-16, and ResNet-50.Table 5
Evaluation of accuracy on each class of Pavia Centre dataset with overall accuracy.
ClassVGG-16Inception-v3ResNet-50ProposedWater99.89100100100Trees94.8595.7895.1095.43Asphalt93.9096.3996.0295.18Self-blocking bricks87.5689.0890.1790.38Bitumen95.4496.7196.5097.84Tiles96.3498.5898.4998.95Shadows95.0495.2094.9995.46Meadows97.6398.0598.5799.65Bare soil96.5096.8998.1799.68Overall accuracy96.9397.9097.5598.94As shown in Table5, the proposed model has obtained better classification results in most of the classes and obtained an overall accuracy of 98.94%, which is 1.04–2.01% improved than the other compared models. The Inception-v3 and ResNet-50 models have achieved some best results in few classes compared to the proposed model, and in most classes, both these models achieved some close results. Figure 4 represents the graphical plot of overall accuracy obtained by the models.Figure 4
Graphical plot of overall accuracy on Pavia Centre dataset.Table6 represents the evaluation of the classification performances of the proposed model based on the accuracy achieved on every class and the overall accuracy of the proposed model using the Pavia University dataset. The Pavia University dataset is divided into 21388 samples for training and 21388 samples for testing in this performance evaluation. As noted in Table 6, the proposed model has obtained better classification results in majority of classes and obtained an average accuracy of 96.05%, which is 0.48–1.2% improved than the other compared models. The Inception-v3 and ResNet-50 models have achieved some best results in few classes compared to the proposed model, and in most classes, both these models achieved some close results. VGG-16 has the least performance in this evaluation. Figure 5 represents the graphical plot of accuracy acquired by the models with Pavia University dataset.Table 6
Evaluation of accuracy on each class of Pavia University dataset with overall accuracy.
ClassVGG-16Inception-v3ResNet-50ProposedAsphalt93.8592.6196.2096.89Meadows96.0496.3597.5298.74Gravel78.3381.2680.0981.63Trees90.1296.7996.6296.08Painted metal sheets99.9010099.85100Bare soil89.4491.7394.2693.38Bitumen85.6090.8386.9089.16Self-blocking bricks86.7589.4892.0792.94Shadows10099.7399.97100Overall accuracy94.8595.1495.5796.05Figure 5
Graphical plot of overall accuracy on Pavia University dataset.Table7 represents the classifications performance evaluation of the proposed system based on the accuracy obtained in every class and the overall accuracy of the proposed system using the Salinas dataset. Salinas dataset is divided into 27065 samples for training and 27064 samples for testing in this performance evaluation.Table 7
Evaluation of accuracy on each class of Salinas dataset with overall accuracy.
ClassVGG-16Inception-v3ResNet-50ProposedBrocoli_green_weeds_198.7499.5699.3899.65Brocoli_green_weeds_297.4999.8999.6099.96Fallow98.5099.9310099.92Fallow_rough_plow98.8799.5099.1099.44Fallow_smooth98.7998.2799.7499.50Stubble96.0096.8099.2899.81Celery98.9499.1399.4199.90Grapes_untrained86.5291.9593.8294.14Soil_vinyard_develop99.87100100100Corn_senesced_green_weeds94.3897.5898.2597.40Lettuce_romaine_4 wk93.8297.7599.0198.57Lettuce_romaine_5 wk97.0699.4999.7299.80Lettuce_romaine_6 wk98.3199.2498.9099.35Lettuce_romaine_7 wk91.9598.8999.0599.26Vinyard_untrained65.1578.3681.2982.90Vinyard_vertical_trellis97.8798.2199.0099.28Overall accuracy90.4392.8995.1596.33The above table shows that the proposed model has obtained better classification results in many classes and obtained an overall accuracy of 96.33%, which is 1.18–5.9% improved than the other compared models. The Inception-v3 and ResNet-50 models have obtained better results in one or two classes than the proposed model. In most classes, both these models obtained some close results. VGG-16 has the least performance in all these experiment evaluations. Figure6 represents the graphical plot of accuracy acquired by the models with the Salinas dataset. Figure 7 represents the graphical plot of accuracy acquired by the proposed model based on the three datasets used in this research.Figure 6
Graphical plot of overall accuracy on Salinas dataset.Figure 7
Proposed model’s accuracy on each dataset.


## 4.1. Dataset Description

In this research, three HSI benchmark datasets such as Pavia Centre, Pavia University, and Salinas, were utilized for evaluating the proposed approach, and their descriptions are as follows.


### 4.1.1. Pavia Centre and University Datasets

Pavia Centre and University are two images captured by the ROSIS sensor during the fly above Pavia, Italy. Pavia University has 103 spectral bands, where Pavia Centre has 102. Pavia Centre has the 1096 × 1096 pixels images, while Pavia University has the 610 × 610 pixels images although some samples in both images are insufficient in information and deleted before processing. The spatial resolution was 1.3 meters. Both image ground truths distinguish nine groups. Table2 represents the description of the Pavia Centre dataset, and Table 3 represents the description of the Pavia University dataset.Table 2
Pavia Centre dataset description.
ClassSampleWater824Trees820Asphalt816Self-blocking bricks808Bitumen808Tiles1260Shadows476Meadows824Bare soil820Total7456Table 3
Pavia University dataset description.
ClassSampleAsphalt6631Meadows18649Gravel2099Trees3064Painted metal sheets1345Bare soil5029Bitumen1330Self-blocking bricks3682Shadows947Total42776


### 4.1.2. Salinas

This dataset was collected by the airborne visible/infrared imaging spectrometer (AVIRIS) sensor above the city of Salinas Valley, California, USA. This collection incorporates images with the size of 512 × 217 pixels, a spatial resolution of 3.7 meters, and 204 bands after eliminating 20 water absorption bands. The accessible ground references map includes 16 different classes, such as vegetables, vineyards, and bare soil as shown in Table4.Table 4
Salinas scene dataset description.
ClassSampleCorn_senesced_green_weeds3278Celery3579Brocoli_green_weeds_12009Fallow_smooth2678Fallow1976Brocoli_green_weeds_23726Fallow_rough_plow1394Soil_vineyard_develop6203Stubble3959Vineyard_untrained7268Grapes_untrained11271Lettuce_romaine_6 wk916Lettuce_romaine_4 wk1068Vineyard_vertical_trellis1807Lettuce_romaine_7 wk1070Lettuce_romaine_5 wk1927Total54129Figure3 represents the ground truth images of the datasets utilized in this research for evaluation. It is worth noting that the zero value in the ground truth maps displayed in Figures 3(a)–3(c) indicates clutter pixels, which could not be classified as a new class.Figure 3
Ground Truth images: (a) Pavia University; (b) Pavia Centre; (c) Salinas.
(a)(b)(c)


## 4.1.1. Pavia Centre and University Datasets

Pavia Centre and University are two images captured by the ROSIS sensor during the fly above Pavia, Italy. Pavia University has 103 spectral bands, where Pavia Centre has 102. Pavia Centre has the 1096 × 1096 pixels images, while Pavia University has the 610 × 610 pixels images although some samples in both images are insufficient in information and deleted before processing. The spatial resolution was 1.3 meters. Both image ground truths distinguish nine groups. Table2 represents the description of the Pavia Centre dataset, and Table 3 represents the description of the Pavia University dataset.Table 2
Pavia Centre dataset description.
ClassSampleWater824Trees820Asphalt816Self-blocking bricks808Bitumen808Tiles1260Shadows476Meadows824Bare soil820Total7456Table 3
Pavia University dataset description.
ClassSampleAsphalt6631Meadows18649Gravel2099Trees3064Painted metal sheets1345Bare soil5029Bitumen1330Self-blocking bricks3682Shadows947Total42776


## 4.1.2. Salinas

This dataset was collected by the airborne visible/infrared imaging spectrometer (AVIRIS) sensor above the city of Salinas Valley, California, USA. This collection incorporates images with the size of 512 × 217 pixels, a spatial resolution of 3.7 meters, and 204 bands after eliminating 20 water absorption bands. The accessible ground references map includes 16 different classes, such as vegetables, vineyards, and bare soil as shown in Table4.Table 4
Salinas scene dataset description.
ClassSampleCorn_senesced_green_weeds3278Celery3579Brocoli_green_weeds_12009Fallow_smooth2678Fallow1976Brocoli_green_weeds_23726Fallow_rough_plow1394Soil_vineyard_develop6203Stubble3959Vineyard_untrained7268Grapes_untrained11271Lettuce_romaine_6 wk916Lettuce_romaine_4 wk1068Vineyard_vertical_trellis1807Lettuce_romaine_7 wk1070Lettuce_romaine_5 wk1927Total54129Figure3 represents the ground truth images of the datasets utilized in this research for evaluation. It is worth noting that the zero value in the ground truth maps displayed in Figures 3(a)–3(c) indicates clutter pixels, which could not be classified as a new class.Figure 3
Ground Truth images: (a) Pavia University; (b) Pavia Centre; (c) Salinas.
(a)(b)(c)


## 4.2. Performance Analysis

For the performance evaluation, the datasets are divided into 50% to train and 50% to test. As shown in the following tabulations, Table5 represents the accuracy of every class and overall accuracy of total classes based on the classification performance of the proposed model using the Pavia Centre dataset. This Pavia Centre dataset is divided into 3728 samples for training and 3728 samples for testing for the experiment. The classification performance of the proposed model was compared with other existing deep transfer learning approaches like Inception-v3, VGG-16, and ResNet-50.Table 5
Evaluation of accuracy on each class of Pavia Centre dataset with overall accuracy.
ClassVGG-16Inception-v3ResNet-50ProposedWater99.89100100100Trees94.8595.7895.1095.43Asphalt93.9096.3996.0295.18Self-blocking bricks87.5689.0890.1790.38Bitumen95.4496.7196.5097.84Tiles96.3498.5898.4998.95Shadows95.0495.2094.9995.46Meadows97.6398.0598.5799.65Bare soil96.5096.8998.1799.68Overall accuracy96.9397.9097.5598.94As shown in Table5, the proposed model has obtained better classification results in most of the classes and obtained an overall accuracy of 98.94%, which is 1.04–2.01% improved than the other compared models. The Inception-v3 and ResNet-50 models have achieved some best results in few classes compared to the proposed model, and in most classes, both these models achieved some close results. Figure 4 represents the graphical plot of overall accuracy obtained by the models.Figure 4
Graphical plot of overall accuracy on Pavia Centre dataset.Table6 represents the evaluation of the classification performances of the proposed model based on the accuracy achieved on every class and the overall accuracy of the proposed model using the Pavia University dataset. The Pavia University dataset is divided into 21388 samples for training and 21388 samples for testing in this performance evaluation. As noted in Table 6, the proposed model has obtained better classification results in majority of classes and obtained an average accuracy of 96.05%, which is 0.48–1.2% improved than the other compared models. The Inception-v3 and ResNet-50 models have achieved some best results in few classes compared to the proposed model, and in most classes, both these models achieved some close results. VGG-16 has the least performance in this evaluation. Figure 5 represents the graphical plot of accuracy acquired by the models with Pavia University dataset.Table 6
Evaluation of accuracy on each class of Pavia University dataset with overall accuracy.
ClassVGG-16Inception-v3ResNet-50ProposedAsphalt93.8592.6196.2096.89Meadows96.0496.3597.5298.74Gravel78.3381.2680.0981.63Trees90.1296.7996.6296.08Painted metal sheets99.9010099.85100Bare soil89.4491.7394.2693.38Bitumen85.6090.8386.9089.16Self-blocking bricks86.7589.4892.0792.94Shadows10099.7399.97100Overall accuracy94.8595.1495.5796.05Figure 5
Graphical plot of overall accuracy on Pavia University dataset.Table7 represents the classifications performance evaluation of the proposed system based on the accuracy obtained in every class and the overall accuracy of the proposed system using the Salinas dataset. Salinas dataset is divided into 27065 samples for training and 27064 samples for testing in this performance evaluation.Table 7
Evaluation of accuracy on each class of Salinas dataset with overall accuracy.
ClassVGG-16Inception-v3ResNet-50ProposedBrocoli_green_weeds_198.7499.5699.3899.65Brocoli_green_weeds_297.4999.8999.6099.96Fallow98.5099.9310099.92Fallow_rough_plow98.8799.5099.1099.44Fallow_smooth98.7998.2799.7499.50Stubble96.0096.8099.2899.81Celery98.9499.1399.4199.90Grapes_untrained86.5291.9593.8294.14Soil_vinyard_develop99.87100100100Corn_senesced_green_weeds94.3897.5898.2597.40Lettuce_romaine_4 wk93.8297.7599.0198.57Lettuce_romaine_5 wk97.0699.4999.7299.80Lettuce_romaine_6 wk98.3199.2498.9099.35Lettuce_romaine_7 wk91.9598.8999.0599.26Vinyard_untrained65.1578.3681.2982.90Vinyard_vertical_trellis97.8798.2199.0099.28Overall accuracy90.4392.8995.1596.33The above table shows that the proposed model has obtained better classification results in many classes and obtained an overall accuracy of 96.33%, which is 1.18–5.9% improved than the other compared models. The Inception-v3 and ResNet-50 models have obtained better results in one or two classes than the proposed model. In most classes, both these models obtained some close results. VGG-16 has the least performance in all these experiment evaluations. Figure6 represents the graphical plot of accuracy acquired by the models with the Salinas dataset. Figure 7 represents the graphical plot of accuracy acquired by the proposed model based on the three datasets used in this research.Figure 6
Graphical plot of overall accuracy on Salinas dataset.Figure 7
Proposed model’s accuracy on each dataset.


## 5. Conclusion

In this research, the HSI classification model was proposed based on a feature extraction method that combined a squeeze and excitation network with a deep convolutional neural network. The primary objective of this research was to create a deep feature extraction model for HSI classification. Deep networks were capable of extracting spatial and spectral features from HSI data simultaneously, which is a benefit for enhancing the performance of the proposed model. The squeeze and excitation network was proposed in this work for this feature extraction model utilizing CNN. The squeeze and excitation block was designed to improve the representation quality of a convolutional neural network. Three benchmark datasets, such as Pavia Centre, Pavia University, and Salinas, were utilized in the experiment to evaluate the proposed model. The proposed model’s performance was validated by a performance comparison with current deep transfer learning approaches such as Inception-v3, VGG-16, and ResNet-50. The proposed SE-CNN model outperforms the compared models in terms of accuracy across all dataset classes as well as overall accuracy. The proposed model achieved an overall accuracy of 98.94% for Pavia Centre dataset, 96.05% for Pavia University dataset, and 96.33% for the Salinas dataset. The time required to train the proposed model is the primary drawback of this study. In future, this limitation will be solved by improving the preprocessing function and adjusting the number of layers in the network and finding a solution for solving the data imbalance in HSI classification.



---

*Source: 9430779-2022-08-04.xml*