Next Article in Journal
Soil Moisture and Water Transport through the Vadose Zone and into the Shallow Aquifer: Field Observations in Irrigated and Non-Irrigated Pasture Fields
Next Article in Special Issue
The Role of Proximity in Transformational Development: The Case of Resource-Based Cities in China
Previous Article in Journal
City Size and Household Consumption in China
Previous Article in Special Issue
Function Evaluation and Coordination Analysis of Production–Living–Ecological Space Based on the Perspective of Type–Intensity–Connection: A Case Study of Suzhou, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Superpixel-Based Long-Range Dependent Network for High-Resolution Remote-Sensing Image Classification

1
College of Geological Engineering and Geomatics, Chang’an University, Xi’an 710064, China
2
School of Land Engineering, Chang’an University, Xi’an 710064, China
3
Academy of Social Governance, Laboratory of Social Entrepreneurship, Center of Social Welfare and Governance, School of Public Affairs, Zhejiang University, Hangzhou 310058, China
4
School of Economics and Management, Xi’an Aeronautical Institute, Xi’an 710077, China
5
Business School, NingboTech University, Ningbo 315100, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Land 2022, 11(11), 2028; https://doi.org/10.3390/land11112028
Submission received: 1 September 2022 / Revised: 4 November 2022 / Accepted: 7 November 2022 / Published: 12 November 2022
(This article belongs to the Special Issue Territory Spatial Planning toward High-Quality Development in China)

Abstract

:
Data-driven deep neural networks have demonstrated their superiority in high-resolution remote-sensing image (HRSI) classification based on superpixel-based objects. Currently, most HRSI classification methods that combine deep learning and superpixel object segmentation use multiple scales of stacking to satisfy the contextual semantic-information extraction of one analyzed object. However, this approach does not consider the long-distance dependencies between objects, which not only weakens the representation of feature information but also increases computational redundancy. To solve this problem, a superpixel-based long-range dependent network is proposed for HRSI classification. First, a superpixel segmentation algorithm is used to segment HRSI into homogeneous analysis objects as input. Secondly, a multi-channel deep convolutional neural network is proposed for the feature mapping of the analysis objects. Finally, we design a long-range dependent framework based on a long short-term memory (LSTM) network for obtaining contextual relationships and outputting classes of analysis objects. Additionally, we define the semantic range and investigate how it affects classification accuracy. A test is conducted by using two HRSI with overall accuracy (0.79, 0.76) and kappa coefficients ( κ ) (0.92, 0.89). Both qualitative and quantitative comparisons are adopted to test the proposed method’s efficacy. Findings concluded that the proposed method is competitive and consistently superior to the benchmark comparison method.

1. Introduction

Remote-sensing images are generally ground images taken by sensing equipment mounted on an air or space vehicle [1,2]. The spatial resolution of the images has increased steadily as remote-sensing sensors have advanced. These images greatly benefit a wide range of fields, including urban spatial planning, land surveying, environmental monitoring, etc. [3,4]. However, remote-sensing images with high spatial resolution usually contain complex ground targets and objects of different scales [5,6,7]. They have scale dependency, which further presents a challenge to HRSI classification. According to classification strategies, they can be classified into pixel-based and object-based classification methods [8].
The pixel-based classification methods usually determine the category to which each pixel belongs based on the grayscale information of image elements [9]. Depending on the level of automation, they can be categorized into supervised and unsupervised classifications. Unsupervised methods include the iterative self-organizing analysis algorithm [10], K-means [11,12], and fuzzy clustering [13]. Zhang et al. [14] proposed a K-means-based framework to extract effective feature-learning representations for remote-sensing image classification. Common methods for supervised classification include the maximum likelihood method [15], multilayer perceptron [16], SVM [17] and random forest classifier [18]. Dong et al. [19] combined conditional random fields, SVM and random forests to propose a method based on multi-model fusion for high-resolution remote-sensing applications. Though the aforementioned pixel-based classification approaches are more appropriate for high-spectral resolution remote-sensing images with a greater number of wavebands, it is still a challenge for these methods to capture the long-range contextual relationship between pixels.
Further, with the development of deep learning, semantic segmentation has been successfully applied in HRSI classification, such as in fully connected network [20], U-Net [21], DeepLab [22]. These methods are the data-driven approach to learning high-level features of pixels and establish contextual relationships between pixels, which compensate for the limitations of pixel classification based on manual features. Theoretically, these methods can fit all the pixel features to be classified. However, these methods rely heavily on the richness of the training data. For HRSI, the detailed information (texture features, geometric features and spatial features) is much richer due to the richness [23]. These pixel-based classification methods are prone to the “pretzel phenomenon” and generate misclassification. To solve these problems, many research works have been carried out to reduce the pretzel effect by increasing the network depth, redesigning the network structure and adding Markov random fields for post-processing [24].
Compared with the pixel-oriented approach, the object-oriented method is mostly applied for HRSI classification [25]. They apply the strategy of segmentation and then classification for HRSI classification, which compensates for the pretzel effect caused by pixel-based classification methods. Kim et al. [26] analyzed the role of texture, scale and objects in image classification of high-resolution aerial images, and the results showed that the object-based multi-scale classification algorithm had the highest accuracy. Ma et al. [27] analyzed the factors affecting the accuracy of object-oriented classification, and found that the optimal segmentation scale was related to the image spatial resolution and the study area, and that random forests performed better in object-based classification. Zheng et al. [28] proposed an object-oriented Markov random field model. Specifically, this method built the weighted-region neighborhood graph based on region size and edge information as feature information, and then achieved semantic segmentation by probabilistic inference from random fields. However, the above object-based approaches still hardly meet the accuracy needs of HRSI classification, despite these methods improving the classification accuracy.
Currently, many studies use object-based segmentation integrating a deep neural network approach for HRSI classification. These methods avoid complex artificially designed features and improve the classification accuracy. Hong et al. [29] proposed a depth feature-based remote-sensing image classification method using multi-scale object analysis and a convolutional neural network (CNN) model. Zhou et al. [30] proposed a fine-grained functional-area classification method based on segmented objects and convolutional neural networks, and combined frequency statistics methods to identify the functional classes of basic units. Although the above object-based classification methods can obtain higher accuracy using deep learning networks, they fail to determine segmentation scale due to the network output size, which can easily cause over-segmentation or under-segmentation problems.
Superpixel segmentation, which situates neighboring pixels into uniformly distributed irregular pixel blocks, has been successfully applied to HRSI classification. Lv et al. [31] proposed a deep-learning method based on CNN and energy-driven sampling for HRSI classification. Li et al. [32] adopted a deep neural network method for objects’ standardized segmentation for HRSI classification. These superpixel-based methods can effectively delineate and map the features of high spatial resolution images. However, these methods only perform feature extraction at multiple scales, which not only increases the information redundancy, but also increases the non-separability between features. Figure 1 shows the current deep neural network approach using superpixels, s 1 , s 2 , . . . , s i scales are sampled and these data are simultaneously used as the input to feature extraction for this superpixel block classification. This approach only obtains a category, while not obtaining the long-range dependencies of superpixel blocks.
Although some progress has been made in deep learning for HRSI classification based on superpixel segmentation, it is still worth exploring. The main problems faced are as follows. (1) Semantic scale. HRSI usually contains features at different scales, and using images of fixed scale range as input will increase the burden of heterogeneous feature representation of superpixel objects at different scales by the network. (2) Long-distance dependence. The superpixel objects only contain pixels with homogeneity in a smaller range, and they fail to determine the class to which the object belongs using only one superpixel object because of the homospectral and heterospectral, and the long-distance dependency between surrounding objects.
To tackle the above problems, this paper proposes a long-distance dependent deep neural network structure for HRSI classification. The main contributions are: (1) For semantic scaling, we propose a multiple-channel all-inclusive shared deep neural network considering the multiple scales between different superpixel objects. A larger range of superpixel objects is used as input, and each object is used as a feature extraction unit, which enhances the feature contribution of each segmented object to feature classification. (2) For long-distance dependence, we design a deep neural network with long-range dependencies. A mesh of contextual correspondences between input objects is established, and contextual dependencies between surrounding distant objects are incrementally enhanced, while the class of each superpixel object is determined.

2. Method

2.1. Methodological Framework

Figure 2 illustrates the technical flowchart of HRSI classification based on a long-range dependent deep neural network, and the process is divided into three main parts: (1) Superpixel segmentation. Superpixel segmentation is performed by simple linear iterative clustering (SLIC) on HRSI to obtain superpixel segmentation objects. The centroid of objects is defined as the center of the cropped patches to generate the data as the network input. (2) Model training and classification. A deep neural network with long-range dependencies is designed. It can extract the contextual relationships between objects for the category output as well as the features of the input data and their long-range dependencies. We divide the objects into training and test data, which are used for the training and prediction of the model, respectively. (3) Classification accuracy analysis and semantic range selection. For HRSI, the trained model’s prediction accuracy is examined. The feedback-designed network structure is utilized to determine the semantic range of the network’s input superpixel objects (the number of input objects) based on accuracy analysis.

2.2. Superpixel Segmentation

Simple linear iterative clustering (SLIC) is a superpixel segmentation algorithm. First, the color space of the image is transformed into CIELab, and the initial clustering centers are obtained by sampling on S pixels apart. To produce superpixels of approximately the same size, the distance of the clustering points is set to S = N k , where k is the desired number of superpixel sizes. The clustering centers need to be moved to the lowest gradient in the 3 × 3 domain. The specific computational steps are followings.
(1)
The seed points (clustering centers) are initialized according to the set number of superpixels, and the seed points are distributed evenly within the image. Suppose the image has a total of N pixel points, pre-segmented into k superpixels of the same size; then, the size of each superpixel is N k , and the distance (step) of adjacent seed points is approximated as S = N k .
(2)
The seed points are reselected in the n × n neighborhood, the gradient values of all pixel points in that neighborhood are calculated, and the seed points are moved to the place with the smallest gradient in that neighborhood.
(3)
Each pixel point is assigned class labels within the neighborhood around each seed point. The search range is limited to 2 S × 2 S , which can accelerate the convergence of the algorithm. The desired superpixel size is S × S , and the search range is 2 S × 2 S .
(4)
Distance metric. The distance metric includes color and spatial distance. For each searched pixel point, they are calculated as the distance to that seed point. The distance is calculated as follows.
d c = l j l i 2 + a j a i 2 + b j b i 2 d s = x j x i 2 + y j y i 2 D = d c N c 2 + d s N s 2
where d c represents the color distance, d s represents the spatial distance, and N s is the maximum spatial distance within a class, defined as N s = S = N k , which is applicable to each cluster. The maximum color distance N c varies both with pictures and with clusters, so we take a fixed constant m instead. The final distance metric D is as follows.
D = d c m 2 + d s S 2
Since each pixel point is searched by multiple seed points, each pixel point is given a distance from the surrounding seed points, and the seed point corresponding to the minimum value is taken as the clustering center of that pixel point.
(5)
Iterative optimization. The above steps are iterated continuously until the error converges.
Given the limitations of CNN, the input data must be rectangular image patches. Although objects have homogeneous regions, they are irregularly polymorphic in size. To unify the structure of the input data, we crop patches containing superpixel objects as input. First, the position of the centroid of the irregular graph of each superpixel partition object is calculated. Second, the image blocks are cropped at equidistant positions above and below the center of mass according to the size of the cropped image blocks. The detailed formula for calculating the centroid of the irregular image is as follows.
( x ¯ , y ¯ ) = i = 1 n A i x ¯ i i = 1 n A i , i = 1 n A i y ¯ i i = 1 n A i
where x and y are the position coordinates of each pixel within the superpixel object and A is the number of pixels of the irregular graph. Figure 3 depicts the flowchart of cropping.

2.3. Long-Range Dependent Network

Figure 4 depicts the framework of a deep neural network with long-range dependencies. The network has two components: feature extraction based on a multichannel CNN and long-distance dependent mapping based on LSTM. The feature-extraction structure of the multichannel-based CNN takes superpixel objects as input, and the network structure of each of their channels is the same. Weight sharing is used to reduce the number of parameters in the training phase. The network is first operated by a convolution and activation function, and then output to Stage-1, Stage-2, Stage-3 and Stage-4 for feature extraction, respectively. Feature reuse is performed between each stage, i.e., the input and output of each stage are summed as the final output. Each stage is composed of CNN, activation function, normalization, and pooling layers. The size of the convolution kernel is set to 3 × 3 , and the step size of the convolution operation is designed to be 2 × 2 to reduce the size of the feature map in Stage-2 and Stage-4, respectively. Finally, 2 × 2 pooling is used to reduce the size of the feature map to 1/8 of the original size, and reshape it to 1 × 1 / 8 as the input of the LSTM. The detailed parameters of the feature-extraction structure of the multichannel-based CNN are shown in Table 1.
LSTM is used in the long-distance dependent mapping structure to handle long-distance relationships in the sequence of segmented objects, remembering the contextual information between each superpixel-segmented object feature information. The LSTM “gate” mechanism is designed as a way to allow information to be selected for passage. It consists of an oblivion gate, an input gate, and an output gate, which determines which information needs to be discarded from the cell state and which new information can be stored in the state as the final output. Suppose the input i t 1 of the previous layer and the sequence data i t to be input in this layer; the output is obtained by an activation function. The output values range from [0, 1] from the forgetting gate, which indicates the probability that the information state of the previous layer is forgotten, with 1 being completely retained and 0 being completely discarded.
o t = σ W f · i t 1 , x t + b f
The input gate contains two parts, each consisting of two activation functions, and the outputs are i t , i t , where i t are taken in the interval [0, 1], indicating the degree to which the information in i t is retained. i t × i t denotes the new information that is retained in this layer, and the new information can be updated to the cell state i t in this layer.
i t = σ W f · i t 1 , x t + b i
i t = tanh W c · i t 1 , x t + b c
C t = f t i t 1 + i t i t
o t = σ W o h t 1 , x t + b o
h t = o t tanh C t
The output gate is used to control the amount of information being filtered in this layer. An activation function is first used to obtain a o t with values in the interval [0, 1], followed by multiplying the cell state with o t after processing it through another activation function, which is the output h t of this layer.

3. Data and Parameter Settings

3.1. Data

WHDLD [33] is a densely labeled dataset which can be used for remote-sensing image retrieval and pixel-based tasks such as remote-sensing classification. We use the pixels of each image with the following five category labels, namely, buildings, roads, bare soil, vegetation and water.
Zurich Summer dataset [34] is a densely labeled dataset of high-resolution remote-sensing image classification which can be used for the semantic segmentation of remote-sensing images. We choose the image file (named mediumresidential) with four category labels, namely, building, road, woodland and vegetation, for training.
Figure 5 shows the two test HRSI used in the experiment. Both images consist of four multispectral bands (including red, green, blue and near infrared). The two images were acquired by two different sensors with rich coverage of feature classes. Therefore, it is suitable to test the generalization ability of our proposed method for classification.
Figure 5a is a Gaofen-2 image (GF) with a resolution of 1 m and a size of 1200 × 1120, located in Wuhan, China. Figure 5c is a QuickBird satellite image (QB) with the location of a Zurich suburb, resolution of 0.61 m and size of 1100 × 1320. Figure 5b,d represent the data used for training and testing. The main feature classes covered on QB are buildings, roads, vegetation, woodland and water. Since these features constitute different proportions of the whole image, the same number of samples cannot be selected for training according to each category. Therefore, the training data is selected based on the proportion of each category in the whole image. In QB, the proportion of buildings is high, and the proportion of water bodies and woodland is low. The proportion of the training sample for each feature in the whole training sample should be selected to be similar to the proportion on the original image. Similarly, in the GF, the samples are selected according to the same strategy.

3.2. Parameter Settings

The accuracy of classification is significantly influenced by the settings of the superpixel-segmentation parameters. Pixels within superpixel-segmented objects act as homogeneous classes; therefore, it is necessary to experiment with the superpixel parameters of the images. Increasing the number of superpixel objects can increase the homogeneity of pixels within each object, which obviously requires a long time for image classification. Therefore, a balance between the number of objects and the tightness is needed to reduce the number of segmented objects while ensuring the homogeneity of each segmented object. Due to the limitation of the deep neural network, it must be ensured that the cropped image patches contain the segmented superpixel objects. Therefore, this requires increasing the tightness to ensure that each superpixel object has approximately the same shape.
Figure 6 shows the partial segmentation effect of the two images with different segmentation parameters. The number of segmented superpixels is determined according to the size of the two images, and the tightness parameter is set to 1–60 for comparison. Figure 6a–e show the segmentation results for tightness of 1, 15, 30, 45, and 60. For GF, the tightness coverage parameter is set between 15–30 to obtain a better segmentation effect. This is due to GF covering complex features with more mixed pixels, and the reduced tightness can ensure the heterogeneity of each segmented object. For QB, the covered features are mainly artificial buildings with distinct geometries, and increasing the tightness can obviously obtain a higher homogeneity of segmented objects. The results show that the tightness parameter of QB image when set between 30–45 is more suitable.

3.2.1. Semantic-Range Selection

The choice of semantic range is the basis for establishing the superpixel distant dependencies. A circle with different radius sizes is defined as the semantic range, i.e., superpixel objects contained within the circle as the input to the network. Figure 7 shows the superpixel inputs to the network for classification contained in circles of different radii. Theoretically, the smallest semantic range should include one superpixel object. The maximum semantic range should be smaller than the size of the remote-sensing image. A larger semantic range can obtain the superpixel-object long-range dependency. However, a larger number of inputs will inevitably increase the difficulty of training and the complexity of computation. Therefore, choosing the appropriate semantic range is very important for remote-sensing image classification. In Section 4.2, we analyze in detail the effect of semantic range on classification accuracy for selecting the appropriate semantic range.

3.2.2. Model-Training Parameters Setting

Different semantic ranges correspond to the number of inputs; therefore, a network model with multiple channels is required to be designed for feature extraction and the classification of superpixel objects. The network parameters for each channel are the same as in Table 1. The parameters are shared between multiple channels during the training and prediction phases, which can reduce the number of parameters of the model and improve computational efficiency. The model is implemented by the Pytorch library and all the experiments are implemented on DUAL-RTX 2070-O8G-EVO. The optimization model of the network is used with Adam and the learning rate, epoch time and batch size were set to 0.0001, 4000 and 32.

4. Experiment and Analysis

In this section, we first perform a qualitative and quantitative comparison with competitive methods (OBIA-SVM, Superpixel-DCNN and DeepLab v3 [35]) on two test images to validate the performance of our proposed method. OBIA-SVM is implemented using the eCognition Developer 9.0. In GF and QB, the optimal scale, shape and compactness obtained by human visual judgment are set to the parameters of 60, 0.6, 0.8 and 60, 0.85, 0.9, respectively. For Superpixel-DCNN and the proposed, they both use SLIC for superpixel segmentation. Superpixel-DCNN uses the network structure in Table 2, for comparison. DeepLab v3 follows the parameter settings of [35] for training and testing on the dataset.
The classification accuracy of each category in the image, overall accuracy (OA) and statistical Kappa ( κ ) coefficients are used as evaluation metrics. Taking the binary classification as an example, the detailed calculation is given below.
O A = TP + TN T
where T is the total number of pixels in the accuracy assessment.
κ = T ( TP + TN ) T 2
where ∑ is the chance accuracy represented by (TP + FP)(TP + FN) + (FN + TN)(FP + TN).

4.1. Classification Results

Classification results.The quantitative comparison results of our proposed method and competitive methods on two images are provided in Table 3. For GF, the proposed method obtains the highest classification accuracy with OA of 0.79 and κ of 0.76. The OA and κ for OBIA-SVM, Superpixel-DCNN and DeepLab v3 are (0.57, 0.51), (0.65, 0.62), and (0.70, 0.68), respectively. For each category of classification accuracy, the effectiveness of the proposal is further demonstrated. The accuracy of buildings, roads, vegetation, water, and bare soil increased by 0.15, 0.41, 1.05 and 0.34 compared to OBIA-SVM, with especially vegetation and water increasing more significantly. The accuracy of road increased dramatically by 0.15 in comparison to Superpixel-DCNN, and the accuracy of buildings, vegetation, water, and bare soil increased by 0.17, 0.25, and 0.23, respectively. The proposal increased considerably by 0.22 for buildings and slightly by 0.08, 0.08, 0.13 and 0.13 for roads, vegetation, water, and bare soil as compared to DeepLab v3.
The proposal obtained the highest classification accuracy on QB with OA and κ higher than GF, as shown in Table 4. On QB, the OA of the proposed method is 0.92 and κ is 0.89. The OA and κ of the proposed method increased by (0.48, 0.62), (0.26, 0.21) and (0.19, 0.20), respectively, in comparison to OBIA-SVM, Superpixel-DCNN, and DeepLab v3. The overall performance ranking of classification accuracy is OBIA-SVM < Superpixel-DCNN < DeepLab V3 < the proposed method. The classification performance of each category on QB proved the effectiveness. The accuracy of each category was much higher than OBIA-SVM, at 0.91, 0.90, 0.89, 0.93, and 0.94, respectively. The proposal improved more considerably than Superpixel-DCNN for roads and water, by 0.20 and 0.13, respectively. However, for buildings and vegetation, their classification accuracy increased less. For buildings, roads, woodland, vegetation and water, the accuracy increased by 0.21, 0.20, 0.30 and 0.08 relative to DeepLab v3, respectively.
Discussion and analysis. The proposed method had a higher classification performance than OBIA-SVM, Superpixel-DCNN and DeepLab v3 on two images. OBIA-SVM method used the strategy of object-based segmentation before classification for remote-sensing image. However, it was limited by the selection of segmentation parameters and the feature representation of the examined objects, which resulted in classification outcomes which were noticeably inferior to those of other approaches, as shown in Figure 8 and Figure 9. The influence of the SVM features allowed confusing each category, despite the fact that OBIA-SVM could obtain fine boundaries on specific feature categories, such as buildings and roads. For example, roads on QB were misclassified as buildings, with water bodies misclassified into vegetation. Superpixel-DCNN classified each segmented object using a deep neural network after applying a superpixel-segmentation approach to obtain a precise boundary. Compared to OBIA-SVM, Superpixel-DCNN was more accurate in classifying objects. Since Superpixel-DCNN only performed feature extraction on a segmented object, lacking contextual information between objects, it achieved lower classification accuracy obtained on the water in GF and roads in QB. This was owing to the large scale of these objects, the difficulty of characterizing the linear properties of roads using only one segmentation object, and the lack of information transmission between upper and lower objects in Superpixel-DCNN. DeepLab v3 was an end-to-end classification method based on deep neural networks, which was a semantic segmentation algorithm. However, DeepLab v3 had lower classification accuracy than our proposed method in the classification of small samples. This is a result of the DeepLab v3 algorithm’s ability to classify each individual pixel and need for a wide range of training data. Deeplab v3 provided a lower classification accuracy due to the homospectral and heterospectral on high-resolution remote-sensing images.
The superpixel-based proposal broadened the semantic range of the data input and could identify long-distance semantic relationships between numerous segmented objects, which enabled our proposed method to correctly classify large-scale features such as water in GF and roads in QB. Additionally, objects were input into CNN and LSTM networks, which not only gathered high-level features of objects, but also increased the effectiveness of classification.
Finally, the speed of classification is also important. Since the comparison methods OBIA-SVM, Superpixel-DCNN and the proposed method all used a segmentation first and then classification strategy, the time of the optimal segmentation parameters was difficult to quantify. For a fair comparison, we tested the processing speed (sec/sample) of the classification models on the test images. All evaluations were based on an Intel(R) i7-7700 CPU and a DUAL-RTX 2070-O8G-EVO. Table 5 compared the estimated time complexity of the different methods. Although Superpixel-DCNN, DeepLab v3 and the proposed method spent more time on the CPU, the compatibility with the GPU made these methods the fastest in processing speed.

4.2. The Effect of Semantic Range on Classification Accuracy

The semantic range was quantified as the number of input segmentation objects. Superpixel segmentation objects were introduced in Section 3.2, and they were encircled by various-sized circles. The number of input objects (N) was chosen as 9, 12, 16, 20, 25, and 36 for measuring the classification performance to investigate the impact of various semantic ranges on classification accuracy. The overall classification maps for GF and QB with various semantic ranges are displayed in Figure 10. The classification maps for roads and water bodies in GF exhibited discontinuity, or misclassification, at smaller semantic ranges. With no essential change in classification accuracy, plane buildings and vegetation display extremely good color continuity. On GF, bare soil had mixed classifications with buildings and an overall linear and plane shape. Bare soil exhibited an increase in correct classification on GF as the semantic range widens. For QB, the classification map exhibited a similar phenomenon to GF. The QB contained a large number of homogeneous regions such as vegetation, water, and textures on QB. They displayed less complexity while maintaining excellent classification accuracy at a constrained semantic range. Since roads and buildings were linear types on QB, the classification map showed the same class color constancy as the semantic range expands. On QB, woodland was plane and point-like, and the classification differed dramatically from the maps for smaller semantic ranges.
Figure 11 depicted the accuracy of each category in the two test images. Overall, as the semantic range widens, the classification accuracy of each category tended to steadily increase, with various categories showing different patterns in classification accuracy. The overall results were similar to conclusions obtained above. As the semantic range widens, the classification accuracy for roads and water in GF increased from 0.38 to a maximum value of 0.79, with a more significant growth between 9 and 16. For buildings and vegetation, their classification accuracy increased insignificantly. For bare soil, the classification accuracy trend was between linear features (roads, water) and plane features (buildings, vegetation). Both the roads in QB and GF followed the same pattern, whereby the classification accuracy increased as the semantic range raised, peaking at values between 20 and 25. While most of the buildings in QB exhibited a linear type, its classification accuracy was significantly affected by the semantic range. Since the accuracy change curves for vegetation and water were nearly constant, it was possible to classify these objects accurately within a limited semantic range. The classification accuracy for vegetation on QB reached its highest level in the semantic range of 16–25.
Result analysis. Objects with linear shapes are subject to larger semantic range images. This might be because there were not as many objects with the linear category within a certain semantic range, which caused a lack of long-distance dependencies between the obtained objects and misclassification. We designed the semantic range with plane circular-selection segmentation objects as input, which had a major impact on why the classification progress of faceted objects was less impacted by the semantic range. Strong dependencies between plane segmentation objects allowed the network to construct a heterogeneous representation between them. Additionally, the homogeneity of plane objects was more influenced than that of the semantic range, as was the case for vegetation and water with high homogeneity in QB. Contrarily, since distinct buildings in GF were impacted by mixed image elements, they were heterogeneous and required a wider semantic range for classification.

4.3. Ablation Studies for Network Configuration

An ablation study of the long-range dependency network was conducted to evaluate the efficiency of joining various network operations to assist in the design of classification networks. A multichannel convolutional network and an LSTM constituted the classification network, whose primary function was to identify the long-distance dependencies among the segmented input objects. Consequently, the multichannel convolutional network serves as the foundational network (MCCB). To evaluate the performance of the classification, MCCB and LSTM networks were combined, either with or without configured longitudinal connections (NL, YL). Additionally, we contrasted the impact of various LSTM layer counts (l = 3, l = 4, l = 5) on the dependencies between the input objects. The detailed network combinations and classification results were shown in Table 6. The networks for each combination allowed the random training variables (epoch time, learning rate, batch size) to be deterministic. The effectiveness of the configured networks was evaluated using the classification accuracy of each category on the evaluation test images.
The experimental results showed that MCCB performed significantly worse than MCCB+NL, MCCB+YL, MCCB + YL (l = 3), MCCB + YL (l = 4) and MCCB + YL (l = 5) for classification on all features. Despite the fact that the input objects in MCCB shared the same parameters, they lacked the contextual relationship between the objects, which led to a lower level of classification accuracy. The classification performance of MCCB + NL was slightly higher than that of MCCB overall. This was because MCCB+NL merely expanded the underlying network by a laterally connected LSTM and did not perform information dependency extraction between the input objects. Therefore, the overall classification performance was lower than that of MCCB + YL.
MCCB + YL (l = 3), MCCB + YL (l = 4) and MCCB + YL (l = 5) established long-range dependencies for the input objects and had higher classification accuracy on the whole. From the classification results of MCCB + YL (l = 4) and MCCB + YL (l = 5), increasing the number of layers of the network did not significantly increase the classification accuracy. This indicated that the LSTM’s layer number had less of an impact on classification accuracy. The accuracy affecting the classification of features was mainly dominated by the longitudinally connected LSTM.

5. Conclusions

In this paper, a long-range-dependent deep neural network is constructed for high-resolution remote-sensing image classification. A weight-sharing multi-channel feature-extraction network was proposed to alleviate the network’s training burden. An LSTM-based mesh mapping network was developed to obtain the long-range dependencies between the input objects. The experimental findings demonstrate that the proposed long-range dependency network structure can obtain contextual relationships between the input objects. The proposed network provides a heterogeneous representation for segmented objects obtaining high classification accuracy. It was tested using two high-resolution remote-sensing images, and their OA and κ are 0.79, 0.76 and 0.92, 0.89, respectively. In addition, we defined the concept of semantic range and explored the effects of different semantic ranges for different categories. The results show that linear features and non-homogeneous plane objects require a larger semantic range. This opens up additional solutions for the joint use of superpixel segmentation and deep neural networks.

Author Contributions

Conceptualization, L.L. and Q.M.; methodology, L.L.; software, Y.J.; validation, Y.J., L.L. and L.H.; formal analysis, Y.J. and L.H.; resources, L.L.; writing—original draft preparation, L.L.; writing—review and editing, Y.J., Y.Z. and Y.Y; funding acquisition, L.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by Grant information: This paper was supported by the National Social Science Foundation of China(21&ZD184).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We are grateful to those involved in data processing and manuscript writing revision.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yao, H.; Qin, R.; Chen, X. Unmanned aerial vehicle for remote sensing applications—A review. Remote Sens. 2019, 11, 1443. [Google Scholar] [CrossRef] [Green Version]
  2. Lv, Z.; Liu, T.; Benediktsson, J.A.; Falco, N. Land cover change detection techniques: Very-high-resolution optical images: A review. IEEE Geosci. Remote Sens. Mag. 2021, 10, 44–63. [Google Scholar] [CrossRef]
  3. Zhang, C.; Pan, X.; Li, H.; Gardiner, A.; Sargent, I.; Hare, J.; Atkinson, P.M. A hybrid MLP-CNN classifier for very fine resolution remotely sensed image classification. ISPRS J. Photogramm. Remote Sens. 2018, 140, 133–144. [Google Scholar] [CrossRef] [Green Version]
  4. Zhang, C.; Yue, P.; Tapete, D.; Jiang, L.; Shangguan, B.; Huang, L.; Liu, G. A deeply supervised image fusion network for change detection in high resolution bi-temporal remote sensing images. ISPRS J. Photogramm. Remote Sens. 2020, 166, 183–200. [Google Scholar] [CrossRef]
  5. Yao, X.; Feng, X.; Han, J.; Cheng, G.; Guo, L. Automatic weakly supervised object detection from high spatial resolution remote sensing images via dynamic curriculum learning. IEEE Trans. Geosci. Remote Sens. 2020, 59, 675–685. [Google Scholar] [CrossRef]
  6. Ghaffarian, S.; Valente, J.; Van Der Voort, M.; Tekinerdogan, B. Effect of attention mechanism in deep learning-based remote sensing image processing: A systematic literature review. Remote Sens. 2021, 13, 2965. [Google Scholar] [CrossRef]
  7. Kotaridis, I.; Lazaridou, M. Remote sensing image segmentation advances: A meta-analysis. ISPRS J. Photogramm. Remote Sens. 2021, 173, 309–322. [Google Scholar] [CrossRef]
  8. Pan, X.; Zhang, C.; Xu, J.; Zhao, J. Simplified object-based deep neural network for very high resolution remote sensing image classification. ISPRS J. Photogramm. Remote Sens. 2021, 181, 218–237. [Google Scholar] [CrossRef]
  9. Fang, Y.; Xu, L.; Peng, J.; Yang, H.; Wong, A.; Clausi, D.A. Unsupervised Bayesian classification of a hyperspectral image based on the spectral mixture model and Markov random field. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2018, 11, 3325–3337. [Google Scholar] [CrossRef]
  10. Hoque, M.E.; Kipli, K.; Zulcaffle, T.M.A.; Mat, D.A.A.; Joseph, A.; Zamhari, N.; Sapawi, R.; Arafat, M.Y. Segmentation of retinal microvasculature based on iterative self-organizing data analysis technique (ISODATA). In Proceedings of the 2019 International UNIMAS STEM 12th Engineering Conference (EnCon), Kuching, Malaysia, 28–29 August 2019; IEEE: Piscataway, NJ, USA, 2019; pp. 59–64. [Google Scholar]
  11. Ahmed, M.; Seraj, R.; Islam, S.M.S. The k-means algorithm: A comprehensive survey and performance evaluation. Electronics 2020, 9, 1295. [Google Scholar] [CrossRef]
  12. Navin, M.S.; Agilandeeswari, L. Land use land cover change detection using k-means clustering and maximum likelihood classification method in the javadi hills, Tamil Nadu, India. Int. J. Eng. Adv. Technol. (IJEAT) 2019, 9, 51–56. [Google Scholar] [CrossRef]
  13. Ruspini, E.H.; Bezdek, J.C.; Keller, J.M. Fuzzy clustering: A historical perspective. IEEE Comput. Intell. Mag. 2019, 14, 45–55. [Google Scholar] [CrossRef]
  14. Zhang, F.; Du, B.; Zhang, L.; Zhang, L. Hierarchical feature learning with dropout k-means for hyperspectral image classification. Neurocomputing 2016, 187, 75–82. [Google Scholar] [CrossRef]
  15. Peng, J.; Li, L.; Tang, Y.Y. Maximum likelihood estimation-based joint sparse representation for the classification of hyperspectral remote sensing images. IEEE Trans. Neural Netw. Learn. Syst. 2018, 30, 1790–1802. [Google Scholar] [CrossRef] [PubMed]
  16. Zhang, B.; Zhang, M.; Kang, J.; Hong, D.; Xu, J.; Zhu, X. Estimation of pmx concentrations from landsat 8 oli images based on a multilayer perceptron neural network. Remote Sens. 2019, 11, 646. [Google Scholar] [CrossRef] [Green Version]
  17. Sheykhmousa, M.; Mahdianpari, M.; Ghanbari, H.; Mohammadimanesh, F.; Ghamisi, P.; Homayouni, S. Support vector machine versus random forest for remote sensing image classification: A meta-analysis and systematic review. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2020, 13, 6308–6325. [Google Scholar] [CrossRef]
  18. Izquierdo-Verdiguier, E.; Zurita-Milla, R. An evaluation of Guided Regularized Random Forest for classification and regression tasks in remote sensing. Int. J. Appl. Earth Obs. Geoinf. 2020, 88, 102051. [Google Scholar] [CrossRef]
  19. Dong, L.; Du, H.; Mao, F.; Han, N.; Li, X.; Zhou, G.; Zhu, D.; Zheng, J.; Zhang, M.; Xing, L.; et al. Very High Resolution Remote Sensing Imagery Classification Using a Fusion of Random Forest and Deep Learning Technique—Subtropical Area for Example. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2020, 13, 113–128. [Google Scholar] [CrossRef]
  20. Fu, G.; Liu, C.; Zhou, R.; Sun, T.; Zhang, Q. Classification for high resolution remote sensing imagery using a fully convolutional network. Remote Sens. 2017, 9, 498. [Google Scholar] [CrossRef] [Green Version]
  21. Ren, Y.; Li, X.; Yang, X.; Xu, H. Development of a dual-attention U-Net model for sea ice and open water classification on SAR images. IEEE Geosci. Remote Sens. Lett. 2021, 19, 1–5. [Google Scholar] [CrossRef]
  22. Niu, Z.; Liu, W.; Zhao, J.; Jiang, G. DeepLab-based spatial feature extraction for hyperspectral image classification. IEEE Geosci. Remote Sens. Lett. 2018, 16, 251–255. [Google Scholar] [CrossRef]
  23. Sandborn, A.; Engstrom, R.N. Determining the relationship between census data and spatial features derived from high-resolution imagery in Accra, Ghana. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2016, 9, 1970–1977. [Google Scholar] [CrossRef]
  24. Ghamisi, P.; Maggiori, E.; Li, S.; Souza, R.; Tarablaka, Y.; Moser, G.; De Giorgi, A.; Fang, L.; Chen, Y.; Chi, M.; et al. New frontiers in spectral-spatial hyperspectral image classification: The latest advances based on mathematical morphology, Markov random fields, segmentation, sparse representation, and deep learning. IEEE Geosci. Remote Sens. Mag. 2018, 6, 10–43. [Google Scholar] [CrossRef]
  25. Li, Z.; Li, E.; Samat, A.; Xu, T.; Liu, W.; Zhu, Y. An Object-Oriented CNN Model Based on Improved Superpixel Segmentation for High-Resolution Remote Sensing Image Classification. IEEE J. Sel. Top. Appl. Earth Obs. Remote. Sens. 2022, 15, 4782–4796. [Google Scholar] [CrossRef]
  26. Kim, M.; Warner, T.A.; Madden, M.; Atkinson, D.S. Multi-scale GEOBIA with very high spatial resolution digital aerial imagery: Scale, texture and image objects. Int. J. Remote Sens. 2011, 32, 2825–2850. [Google Scholar] [CrossRef]
  27. Ma, L.; Li, M.; Ma, X.; Cheng, L.; Du, P.; Liu, Y. A review of supervised object-based land-cover image classification. ISPRS J. Photogramm. Remote Sens. 2017, 130, 277–293. [Google Scholar] [CrossRef]
  28. Zheng, C.; Wang, L. Semantic segmentation of remote sensing imagery using object-based Markov random field model with regional penalties. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2014, 8, 1924–1935. [Google Scholar] [CrossRef]
  29. Hong, L.; Zhang, M. Object-oriented multiscale deep features for hyperspectral image classification. Int. J. Remote Sens. 2020, 41, 5549–5572. [Google Scholar] [CrossRef]
  30. Zhou, W.; Ming, D.; Lv, X.; Zhou, K.; Bao, H.; Hong, Z. SO–CNN based urban functional zone fine division with VHR remote sensing image. Remote Sens. Environ. 2020, 236, 111458. [Google Scholar] [CrossRef]
  31. Lv, X.; Ming, D.; Chen, Y.; Wang, M. Very high resolution remote sensing image classification with SEEDS-CNN and scale effect analysis for superpixel CNN classification. Int. J. Remote Sens. 2019, 40, 506–531. [Google Scholar] [CrossRef]
  32. Li, L.; Han, L.; Hu, H.; Liu, Z.; Cao, H. Standardized object-based dual CNNs for very high-resolution remote sensing image classification and standardization combination effect analysis. Int. J. Remote Sens. 2020, 41, 6635–6663. [Google Scholar] [CrossRef]
  33. Shao, Z.; Zhou, W.; Deng, X.; Zhang, M.; Cheng, Q. Multilabel remote sensing image retrieval based on fully convolutional network. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2020, 13, 318–328. [Google Scholar] [CrossRef]
  34. Volpi, M.; Ferrari, V. Semantic segmentation of urban scenes by learning local class interactions. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition Workshops, Boston, MA, USA, 7–12 June 2015; pp. 1–9. [Google Scholar]
  35. Chen, L.C.; Zhu, Y.; Papandreou, G.; Schroff, F.; Adam, H. Encoder-decoder with atrous separable convolution for semantic image segmentation. In Proceedings of the European Conference on Computer Vision (ECCV), Munich, Germany, 8–14 September 2018; pp. 801–818. [Google Scholar]
Figure 1. Illustration of a superpixel segmented object by multi-scale input to deep neural networks for classification, where C i denotes the class of a superpixel object.
Figure 1. Illustration of a superpixel segmented object by multi-scale input to deep neural networks for classification, where C i denotes the class of a superpixel object.
Land 11 02028 g001
Figure 2. Methodology flowchart.
Figure 2. Methodology flowchart.
Land 11 02028 g002
Figure 3. Intuitive description of the superpixel object cropping process. (a) Superpixel segmentation. (b) Cropping of image patches containing superpixel objects.
Figure 3. Intuitive description of the superpixel object cropping process. (a) Superpixel segmentation. (b) Cropping of image patches containing superpixel objects.
Land 11 02028 g003
Figure 4. The proposed network structure.
Figure 4. The proposed network structure.
Land 11 02028 g004
Figure 5. Original images and corresponding samples. (a,c) are the Gaofen-2 and QuickBird images, respectively. (b,d) are the samples selected from the original images for testing.
Figure 5. Original images and corresponding samples. (a,c) are the Gaofen-2 and QuickBird images, respectively. (b,d) are the samples selected from the original images for testing.
Land 11 02028 g005
Figure 6. Superpixel-segmentation results on partial GF and QB by different segmentation parameters, where the tightness parameters of (ae) are 1, 15, 30, 45, and 60, respectively.
Figure 6. Superpixel-segmentation results on partial GF and QB by different segmentation parameters, where the tightness parameters of (ae) are 1, 15, 30, 45, and 60, respectively.
Land 11 02028 g006
Figure 7. Semantic range. Circles of different radii are used to define the semantic range, which contains segmented objects for input.
Figure 7. Semantic range. Circles of different radii are used to define the semantic range, which contains segmented objects for input.
Land 11 02028 g007
Figure 8. Qualitative classification of the competitive and proposed methods on GF.
Figure 8. Qualitative classification of the competitive and proposed methods on GF.
Land 11 02028 g008
Figure 9. Qualitative classification of the competitive and proposed methods on QB.
Figure 9. Qualitative classification of the competitive and proposed methods on QB.
Land 11 02028 g009
Figure 10. Classification of GF and QB with different semantic ranges, where (ae) are under 9, 16, 20, 25 and 36 semantic ranges, respectively.
Figure 10. Classification of GF and QB with different semantic ranges, where (ae) are under 9, 16, 20, 25 and 36 semantic ranges, respectively.
Land 11 02028 g010
Figure 11. Classification accuracy on the two test images for each category.
Figure 11. Classification accuracy on the two test images for each category.
Land 11 02028 g011
Table 1. Network structure parameters.
Table 1. Network structure parameters.
LayerOperationFeature MapSizeKernel SizeStrideActivate
Input-- 32 × 32 --
2× CNN64 32 × 32 3 × 3 1 × 1 ReLU
Stage-12× CNN128 32 × 32 3 × 3 1 × 1 ReLU
Stage-1Max Pooling128 32 × 32 3 × 3 1 × 1 ReLU
Stage-22× CNN256 32 × 32 3 × 3 1 × 1 ReLU
Stage-2Max Pooling256 16 × 16 3 × 3 2 × 2 ReLU
Stage-32× CNN512 16 × 16 3 × 3 1 × 1 ReLU
Stage-3Max Pooling512 16 × 16 3 × 3 2 × 2 ReLU
Stage-42× CNN512 16 × 16 3 × 3 1 × 1 ReLU
Stage-4Max Pooling512 8 × 8 3 × 3 2 × 2 ReLU
FCNFC-256---
Table 2. A confusion matrix. TP is the number of correctly classification C 1 pixels. FN is the number of undetected C 1 pixels. FP is the number of incorrectly extracted C 1 pixels. TN is the number of correctly rejected C 2 pixels.
Table 2. A confusion matrix. TP is the number of correctly classification C 1 pixels. FN is the number of undetected C 1 pixels. FP is the number of incorrectly extracted C 1 pixels. TN is the number of correctly rejected C 2 pixels.
Reference Data
C 1 C 2
Classified data C 1 TPFP
C 2 FNTN
Table 3. Classification results on GF.
Table 3. Classification results on GF.
MethodBuildingRoadBare GroundVegetationWater
Accuracy OA κ
OBIA-SVM 0.65 ± 0.16 0.55 ± 0.18 0.52 ± 0.07 0.39 ± 0.24 0.63 ± 0.06 0.57 0.51
Superpixel-DCNN 0.64 ± 0.17 0.66 ± 0.17 0.62 ± 0.06 0.64 ± 0.21 0.69 ± 0.07 0.65 0.62
Deeplab V3 0.61 ± 0.15 0.72 ± 0.10 0.66 ± 0.08 0.74 ± 0.05 0.75 ± 0.08 0.70 0.68
Proposed 0.75 ± 0.09 0.78 ± 0.18 0.75 ± 0.04 0.80 ± 0.19 0.85 ± 0.07 0.79 0.76
Table 4. Classification results on QB.
Table 4. Classification results on QB.
MethodBuildingRoadWoodlandVegetationWater
Accuracy OA κ
OBIA-SVM 0.63 ± 0.07 0.66 ± 0.14 0.61 ± 0.06 0.46 ± 0.21 0.88 ± 0.05 0.62 0.61
Superpixel-DCNN 0.71 ± 0.06 0.74 ± 0.13 0.70 ± 0.05 0.72 ± 0.24 0.83 ± 0.08 0.73 0.73
Deeplab v3 0.75 ± 0.06 0.75 ± 0.11 0.73 ± 0.08 0.71 ± 0.26 0.87 ± 0.06 0.77 0.74
Proposed 0.91 ± 0.03 0.90 ± 0.13 0.89 ± 0.07 0.93 ± 0.15 0.94 ± 0.06 0.92 0.89
Table 5. Time complexity comparisons on GF and QB. The unit is second (s).
Table 5. Time complexity comparisons on GF and QB. The unit is second (s).
OBIA-SVMSuperpixel-DCNNDeepLab v3ProposedSuperpixel-DCNNDeepLab v3Proposed
BackendCPUGPU
GF 10.27 15.25 9.46 12.85 4.64 3.15 4.06
QB 10.59 15.77 10.19 13.25 4.91 3.53 4.62
Table 6. Ablation study classification results on GF and QB.
Table 6. Ablation study classification results on GF and QB.
DataObjectMCCBMCCB + NLMCCB + YLMCCB + YL (l = 3)MCCB + YL (l = 4)MCCB + YL (l = 5)
Building 0.46 ± 0.06 0.54 ± 0.08 0.75 ± 0.08 0.74 ± 0.07 0.73 ± 0.08 0.76 ± 0.03
Road 0.41 ± 0.06 0.53 ± 0.01 0.78 ± 0.06 0.71 ± 0.08 0.75 ± 0.06 0.75 ± 0.07
GFBare ground 0.56 ± 0.07 0.67 ± 0.03 0.77 ± 0.05 0.70 ± 0.06 0.77 ± 0.07 0.79 ± 0.06
Vegetation 0.62 ± 0.04 0.61 ± 0.01 0.79 ± 0.07 0.77 ± 0.07 0.76 ± 0.08 0.73 ± 0.07
Water 0.68 ± 0.07 0.68 ± 0.00 0.83 ± 0.06 0.78 ± 0.04 0.71 ± 0.08 0.77 ± 0.06
Building 0.59 ± 0.07 0.60 ± 0.06 0.91 ± 0.04 0.92 ± 0.03 0.87 ± 0.07 0.82 ± 0.03
Road 0.57 ± 0.05 0.64 ± 0.06 0.89 ± 0.07 0.87 ± 0.03 0.81 ± 0.07 0.88 ± 0.07
QBWoodland 0.54 ± 0.08 0.60 ± 0.02 0.92 ± 0.05 0.84 ± 0.06 0.82 ± 0.05 0.87 ± 0.04
Vegetation 0.63 ± 0.03 0.66 ± 0.04 0.93 ± 0.06 0.85 ± 0.08 0.84 ± 0.05 0.86 ± 0.01
Water 0.69 ± 0.04 0.70 ± 0.01 0.94 ± 0.05 0.90 ± 0.06 0.89 ± 0.06 0.89 ± 0.02
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, L.; Han, L.; Miao, Q.; Zhang, Y.; Jing, Y. Superpixel-Based Long-Range Dependent Network for High-Resolution Remote-Sensing Image Classification. Land 2022, 11, 2028. https://doi.org/10.3390/land11112028

AMA Style

Li L, Han L, Miao Q, Zhang Y, Jing Y. Superpixel-Based Long-Range Dependent Network for High-Resolution Remote-Sensing Image Classification. Land. 2022; 11(11):2028. https://doi.org/10.3390/land11112028

Chicago/Turabian Style

Li, Liangzhi, Ling Han, Qing Miao, Yang Zhang, and Ying Jing. 2022. "Superpixel-Based Long-Range Dependent Network for High-Resolution Remote-Sensing Image Classification" Land 11, no. 11: 2028. https://doi.org/10.3390/land11112028

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