GF-5 Hyperspectral Data for Species Mapping of Mangrove in Mai Po, Hong Kong

: Hyperspectral data has been widely used in species discrimination of plants with rich spectral information in hundreds of spectral bands, while the availability of hyperspectral data has hindered its applications in many speciﬁc cases. The successful operation of the Chinese satellite, Gaofen-5 (GF-5), provides potentially promising new hyperspectral dataset with 330 spectral bands in visible and near infrared range. Therefore, there is much demand for assessing the e ﬀ ectiveness and superiority of GF-5 hyperspectral data in plants species mapping, particularly mangrove species mapping, to better support the e ﬃ cient mangrove management. In this study, mangrove forest in Mai Po Nature Reserve (MPNR), Hong Kong was selected as the study area. Four dominant native mangrove species were investigated in this study according to the ﬁeld surveys. Two machine learning methods, Random Forests and Support Vector Machines, were employed to classify mangrove species with Landsat 8, Simulated Hyperion and GF-5 data sets. The results showed that 97 more bands of GF-5 over Hyperion brought a higher over accuracy of 87.12%, in comparison with 86.82% from Hyperion and 73.89% from Landsat 8. The higher spectral resolution of 5 nm in GF-5 was identiﬁed as making the major contribution, especially for the mapping of Aegiceras corniculatum . Therefore, GF-5 is likely to improve the classiﬁcation accuracy of mangrove species mapping via enhancing spectral resolution and thus has promising potential to improve mangrove monitoring at species level to support mangrove management.


Introduction
Mangrove forests are tropical and subtropical ecosystems, growing in inter-tidal areas, being the interface of land and oceans [1]. They can provide coastal protection, economic benefit from aquaculture, and significant eco-services, such as habitat provision and carbon sequestration [2,3]. However, mangrove forests have been suffering great habitat loss due to human activities and global climate change [4], which urge us to monitor mangrove at various scales. Monitoring mangrove at species level can provide essential information on biodiversity, which may support mangrove management and protection in terms of biodiversity conservation, ecological succession analysis, biomass and carbon estimation etc.

Study Site
The Mai Po Nature Reserve (MPNR, 113°59'E-114°03´E, 22°28´N-22°32´N) is located in the mouth of Shenzhen River, the northwest of Hong Kong, opposite to the Futian Nature Reserve, Shenzhen (Error! Reference source not found.a). The reserve comprises of Geiwai, freshwater ponds, inter-tidal mudflats, mangroves, reed beds and fishponds, providing habitats for various wildlife (World Wild Fund for Nature -Hong Kong, WWF-HK 2016). Being a key station of the East Asian-Australasian Flyway (EAAF), it services million migratory birds every year. Due to the special location and significant eco-services provision, the reserve has been listed as restricted area via Wild Animals Protection Ordinance in 1975 and designated as a Site of Special Scientific Importance (SSSI) in 1976. In 1983, the nature reserve was developed and run by WWF-HK; it was further designated as a Ramsar Site in 1995.
The mangrove forest in the Mai Po is the largest one in Hong Kong with an area of around 319 ha [14], mainly occupying the core zone of the reserve. Six native and two exotic mangrove species have been found in the reserve. Specifically, Kandelia obovate (KO), Avicennia marina (AM), Aegiceras corniculatum (AC), and Acanthus ilicifolius (AI) are dominant, while Bruguiera gymnorrhiza and Excoecaria agallocha are less seen. It is believed that two exotic species including Sonneratia caseolaris and Sonneratia apetala are floating from the Futian Nature Reserve, and they are treated as invasive species and annually removed by Agriculture, Fisheries and Conservation Department (AFCD) due to potential threats to local species.  The mangrove forest in the Mai Po is the largest one in Hong Kong with an area of around 319 ha [14], mainly occupying the core zone of the reserve. Six native and two exotic mangrove species have been found in the reserve. Specifically, Kandelia obovate (KO), Avicennia marina (AM), Aegiceras corniculatum (AC), and Acanthus ilicifolius (AI) are dominant, while Bruguiera gymnorrhiza and Excoecaria agallocha are less seen. It is believed that two exotic species including Sonneratia caseolaris and Sonneratia apetala are floating from the Futian Nature Reserve, and they are treated as invasive species and annually removed by Agriculture, Fisheries and Conservation Department (AFCD) due to potential threats to local species.

Hyperspectral Data and Preprocessing
The hyperspectral data adopted in this study was acquired by Advanced Hyperspectral Imager (AHSI) loaded by GF-5 satellite at 1:40 pm on 5 October, 2018. It provides Visible/Near-infrared (VNIR) and Shortwave infrared (SWIR) data with spectral range from 400-2500nm. Similar to Hyperion, the spatial resolution of GF-5 hyperspectral data is 30 m and the spectral resolution in SWIR is 10 nm. Due to the higher spectral resolution of 5 nm in VNIR, GF-5 hyperspectral data has double bands in VNIR and provides 330 bands in total (Table 1). The data was delivered as a Level 1 product, which has been preprocessed with radiometric calibration and noise removal by the China Centre for Resources Satellite Data and Application (CRESDA), China (http://www.cresda.com/). The radiometric calibration was conducted with a linear model with absolute coefficients calculated via laboratory experiments by CRESDA, China. After removing the bands without calibration or those contaminated by strips (Table 2), we conducted atmospheric correction using FLAASH module in ENVI 5.3.1 [28]. For the spectral overlap between VNIR (1006.68-1028.98nm) and SWIR (1004..57-1029.85nm), we kept the overlapped data in VNIR rather than SWIR to conserve accurate spectral information. Finally, we obtained 250 bands for GF-5 hyperspectral data (the details can be referred to the Appendix A). Table 2. Excluded bands of hyperspectral data from GF-5 satellite (* the absorption was referred to [28] To save the costs on computation, a region covering the MPNR maximumly was cropped from the Gaofen-5 hyperspectral image (Figure 1b), and fifteen ground control points were picked up from the cropped image for geometric correction with a result of Root Mean Square (RMS) of 0.048068. Because we focus on mangroves, the remaining objects like urban areas were firstly masked by manually outlining the initial region containing mangrove forests. The mudflat and water body were then removed from the initial region by thresholding Normalized Difference Vegetation Index (NDVI), which is calculated as a ratio between a red band around 660 nm and a near infrared band around 860 nm and widely used for primary separation vegetation from non-vegetation or mangrove from non-mangrove [32].

Field Survey and Sample Data
To investigate the real distribution of mangrove species in MPNR, a field survey was conducted along the floating bridge in the Mai Po on November 10, 2015. Based on the field survey, samples were manually identified with visual interpretation of a very high resolution WorldView-3 imagery Remote Sens. 2020, 12, 656 5 of 16 with 1.6 m resolution in multispectral bands and 0.4 m resolution in the panchromatic band. Initially, 416 samples of six mangrove species containing the information on mangrove species and Global Positioning System (GPS) location were obtained [18]. To reduce the effect of the difference of date acquisition between hyperspectral data and samples on results, we identified the altered areas by comparing two very high spatial resolution images in two periods and removed the samples located in the changed areas. Finally, we added extra samples and obtained 293 samples of six classes (Table 3). In particular, AI was further divided into two subclasses (AI 1 and AI 2 ) according to whether leaves had smooth or serrated edges, and KO was also refined with two subsets of KO 1 and KO 2 according to the difference of tree height and leaf density [6,18,28], which may be related to the seaward or landward location [33]. Sonneratia was not selected as samples because of difficult percept of small patches from hyperspectral data after removal by AFCD. These samples were divided into training and testing samples with proportion of 70% and 30% respectively using the strategy of stratified random sampling. Based on the samples, the mean reflectance spectra of six mangrove species can be presented in Figure 2. in the changed areas. Finally, we added extra samples and obtained 293 samples of six classes (Error! Reference source not found.). In particular, AI was further divided into two subclasses (AI1 and AI2) according to whether leaves had smooth or serrated edges, and KO was also refined with two subsets of KO1 and KO2 according to the difference of tree height and leaf density [6,18,28], which may be related to the seaward or landward location [33]. Sonneratia was not selected as samples because of difficult percept of small patches from hyperspectral data after removal by AFCD. These samples were divided into training and testing samples with proportion of 70% and 30% respectively using the strategy of stratified random sampling. Based on the samples, the mean reflectance spectra of six mangrove species can be presented in Error! Reference source not found.. Aegiceras corniculatum AC usually found as undergrowth 19 Figure 2. The mean reflectance spectra of six mangrove species (the scale was set to be 10 for demonstration).

Methods
We designed a comparison experiment to assess the performance of GF-5 hyperspectral data in mangrove species mapping. Hyperion from EO-1 may be for benchmark since its global coverage with high quality and free charge. However, the available Hyperion covering the Mai Po is in 2008 and no more Hyperion data is provided after EO-1 decommission in 2017. The date difference of available Hyperion and Gaofen-5 hyperspectral data makes the comparison between them impossible. To overcome the limitation of available hyperspectral data for comparison with GF-5 hyperspectral data, we Simulated Hyperion (SH) data according to the spectral positions and bandwidth of Hyperion satellite based on GF-5 hyperspectral data with a simulation model. Besides, we also added another benchmark of Landsat 8 to be representative of multispectral data. Specifically, the Landsat 8 image closed to the date when GF-5 hyperspectral data was acquired and without cloud coverage was selected for comparison. The panchromatic band, cirrus band and thermal infrared bands were removed due to irrelevance to vegetation. Detailed experiment design is listed in Error! Reference source not found..  The mean reflectance spectra of six mangrove species (the scale was set to be 10 for demonstration).

Methods
We designed a comparison experiment to assess the performance of GF-5 hyperspectral data in mangrove species mapping. Hyperion from EO-1 may be for benchmark since its global coverage with high quality and free charge. However, the available Hyperion covering the Mai Po is in 2008 and no more Hyperion data is provided after EO-1 decommission in 2017. The date difference of available Hyperion and Gaofen-5 hyperspectral data makes the comparison between them impossible. To overcome the limitation of available hyperspectral data for comparison with GF-5 hyperspectral data, we Simulated Hyperion (SH) data according to the spectral positions and bandwidth of Hyperion satellite based on GF-5 hyperspectral data with a simulation model. Besides, we also added another benchmark of Landsat 8 to be representative of multispectral data. Specifically, the Landsat 8 image closed to the date when GF-5 hyperspectral data was acquired and without cloud coverage was selected for comparison. The panchromatic band, cirrus band and thermal infrared bands were removed due to irrelevance to vegetation. Detailed experiment design is listed in Table 4.

Generation of Simulated Hyperion Data
EO-1 Hyperion data contains 242 bands, including VNIR of 70 bands and SWIR 172 bands, with spectral resolution of 10 nm and spatial resolution of 30 m. With same range from 400 to 2500 nm, GF-5 hyperspectral data has more bands since the spectral resolution was improved from 10 nm to 5 nm in VNIR. Regardless of the response difference for two spectra with an interval of less 5 nm, Hyperion can be viewed as the subset of GF-5 hyperspectral data. Consequently, simulating Hyperion from GF-5 can be regarded as the subset construction. Namely, the band of simulated Hyperion b SH i with wavelength of i can be described with the band from GF-5 hyperspectral data (b GF5 j ) with wavelength of j if the difference between their wavelength is minimum and less than 5 nm, which can be formulated in Equation (1) where F is the set of wavelengths of the valid bands from GF-5 hyperspectral data.

Mangrove Species Classification with Machine Learning Methods
More spectral features provided by Hyperspectral data make better performance possible. Meanwhile, they also promote a high requirement of classifiers due to high dimension [34,35]. Since this study aims to assess the efficiency of spectra, common approaches of dimension reduction, such as Principal Component Analysis (PCA) and Minimum Noise Fraction (MNF) will, because of their spectral characteristics, not apply to our purposes, because the contribution of the original spectral features cannot be tracked after applying the transformation from PCA or MNF. Therefore, the approaches which can inherently handle high-dimension features will be considered.
Random Forests adopts the principle of ensemble learning to vote the results from multiple weaker decision trees [36]. On the one hand, Random Forests randomly select a subset of features from all the features to split tree nodes, so they can deal with high-dimension features as well as maintain the robustness to noises. On the other hand, Random Forests randomly select a subset of samples from all the training samples to construct individual decision tree. When constructing a decision tree, around 63% of the samples are selected for training the tree while the rest (also named out-of-bag, OOB) for calculating an OOB error; and m features are randomly picked up from p features (total number of features). These m most distinguished features are selected for node splitting based on the measurement of Gini coefficient [37]. This procedure iterates until it reaches one of the stopping criteria, such as approaching maximum depth, minimum error threshold and maximum member of decision trees. Therefore, the tree number and feature number m are important for Random Forests construction. In principle, m was set to be √ p [38]. Support Vector Machine (SVM) is another approach which is not sensitive to feature dimension and well adapted to process hyperspectral data [39,40]. Based on structural risk minimization, SVM searches a hyperplane which can maximally separate two parts by projecting the data into a high-dimension linear space. For the case of nonlinearity, a kernel function is used to project the data into a higher Remote Sens. 2020, 12, 656 7 of 16 dimensional space where the projected data become linearly distributed. Moreover, slack variable and penalty factors are introduced to control the bias and make them approximately separated. For the multiclass tasks of mapping multiple mangrove species, the "one-against-one" approach was used to construct many binary SVMs for each possible pair of two classes for training, while a voting strategy selected the class with maximum number of votes for classifying the new data. The details can be referred to in [41]. Kernel function has significant influence on the performance of SVM. In general, linear kernel and Radial Basis Function (RBF) kernel are widely used due to them having fewer parameter settings. We chose RBF for this study in consideration that linear kernel is a special case of RBF [42], for which the Gamma coefficient was set to be the reciprocal of dimension of input feature (HARRIS 2019).

Accuracy Assessment
Confusion matrix is employed to analyze the details on the classification of different classes over the samples, from which overall accuracy, Kappa coefficient, user accuracy and producer accuracy are derived for accuracy assessment in different aspects. In this study, we also adopted the overall accuracy and Kappa coefficient for accuracy assessment [43].

Comparison of Simulated Hyprion and GF-5
To better simulate the Hyperion, we adopted the valid bands from the available data acquired in 2008 [28] and got 153 bands for Hyperion. Due to sensors differences, nine bands with wavelength near 1336.87nm, 1346.97nm, 1488.29nm, 1498.4nm, 2002.8nm, 2012.9nm, 2375.98nm, 2386.08nm, 2396.18 nm are not available in GF-5 hyperspectral data. Therefore, 144 bands including 52 VNIR bands and 92 SWIR bands were simulated. Compared to simulated Hyperion, GF-5 hyperspectral data has nearly three times valid bands (142) in VNIR, while has 108 valid bands in SWIR. The details of simulated Hyperion can be seen in to the Appendix A.

Mangrove Species Mapping with Accuracy Assessment
According to the data availability of GF-5 hyperspectral data and simulated Hyperion, the bands of experiment data including SH, SH-GF5VN, SH-GF5SW, GF-5 are 144, 234, 160 and 250 individually. Together with Landsat 8, these data are fed into both Random Forests and SVM for comparison. The number of decision trees for random forests was set 100, and the penalty factor for SVM with RBF kernel was set 100. To reduce the bias of sample division for training, we generated ten groups of training and testing samples through random selection with proportion of 70% and 30% respectively for ten times. As a result, we can obtain ten overall accuracies and mean overall accuracy for each dataset. The results using SVM and Random Forests can be seen from Figure 3.
Compared to multispectral data of Landsat 8, hyperspectral data including GF-5 and simulated Hyperion demonstrated significant advantage in mangrove species mapping with improvement in mean overall accuracy by 13%-17% and mean Kappa coefficient by 0.16-0.24 using two classifiers. Against simulated Hyperion, the increase in spectral resolution in VNIR (namely SH-GF5VN) or appended data in SWIR (SH-GF5SW) can make contribution to mangrove species mapping, showing growth of mean overall accuracy by 0.19% and 0.07% with Random Forests, while 0.47% and 0.31% with SVM. The mean Kappa coefficient of 0.835 over 0.834 (Random Forests) and 0.74 over 0.73 (SVM) for SH-GF5VN and SH-GF5SW also shows that the improvement from the increase in spectral resolution in VNIR is more significant than that from addictive bands in SWIR. Surprisingly, GF-5 hyperspectral, namely simultaneous enhancement in VNIR and SWIR based on simulated Hyperion, boosts further growth only with Random Forests, while it lowers the overall accuracy with SVM. The influence from classifiers on mangrove species mapping with remote sensing data can also be revealed from the comparison between Random Forests and SVM with same input data. In general, the overall accuracy Remote Sens. 2020, 12, 656 8 of 16 resulted from Random Forests (Figure 3a) is higher than that resulted from SVM (Figure 3b)  (a) (b) Figure 3. The overall accuracy of mangrove species mapping using Random Forests (a) and Support Vector Machine (SVM) (b) with different datasets and ten different groups of training and testing samples (number 1-10) through stratified random sampling with proportion of 70% and 30% respectively; the mean Overall Accuracy (OA) and mean Kappa coefficient were also shown in lines.
The average confusion matrixes for the classification with Landsat 8, SH, and SH-GF5VN were calculated for insight into influence of spectral resolution increasement on mangrove species mapping (Error! Reference source not found.). From the producer accuracy derived from Landsat 8, we can observe that KO1 is difficult to be classify correctly with accuracy of 68.08%, and the confusion between AM and other species is serious with user accuracy of 68.04%, especially with KO1. This situation will be improved with hyperspectral data (simulated Hyperion and SH-GF5VN), showing growth of around 17% and 10% in producer accuracy for KO1 and user accuracy for AM. For the classification with simulated Hyperion and SH-GF5VN, the major difference comes from AC mapping. With the help of high spectral resolution in VNIR (SH-GF5VN), AC can be classified with higher accuracy of 84.00% against 78.00% from SH.
Visually, the misclassification and confusion can be also observed from the mapping (Error! Reference source not found.). In the south part of study area (area A in Error! Reference source not found.a), the majority of KO1 was misclassified into AM, which is consistence with low producer accuracy of 68.08% in KO1. In the area B where narrow zonation of AM was between AI1 and KO1, the confusion happens when using Landsat 8 (low user accuracy of 68.04%); while for area C, AI1 and AM are confused. For the distribution of mangrove species resulting from simulated Hyperion and GF-5 hyperspectral data, they present a similar pattern with a minor difference in area D and E, where a small amount of AM was misclassified into AI1 and more AC were correctly identified. Similarly, the results from SVM with datasets from Landsat 8 (Error! Reference source not found.d), SH (Error! Reference source not found.e) to SH-GF5VN (Error! Reference source not found.f) show the major changes in AI1 and AI2. More specifically, the classification of AI1 and AI2 were refined with fine particles when using hyperspectral data (SH and SH-GF5VN). In contrast, fewer ACs were found from the mangrove species mapping using SVM, and KO located at the seaward side cannot be identified unless SH and SH-GF5VN are used.
Finally, we calculated the areas for each mangrove species in the study area using Random Forests with SH-GF5VN. The results show that the most dominant species is KO1 with area of 103.95±3.56 ha (46% of the mangrove forests in the study area), followed by AM, AI2, AI1, and KO2. The AC species are rare with an area of 6.93±1.07 ha, occupying only 3% of the mangrove forests in the study area. The details can be found in Error! Reference source not found..   The average confusion matrixes for the classification with Landsat 8, SH, and SH-GF5VN were calculated for insight into influence of spectral resolution increasement on mangrove species mapping ( Table 6). From the producer accuracy derived from Landsat 8, we can observe that KO 1 is difficult to be classify correctly with accuracy of 68.08%, and the confusion between AM and other species is serious with user accuracy of 68.04%, especially with KO 1 . This situation will be improved with hyperspectral data (simulated Hyperion and SH-GF5VN), showing growth of around 17% and 10% in producer accuracy for KO 1 and user accuracy for AM. For the classification with simulated Hyperion and SH-GF5VN, the major difference comes from AC mapping. With the help of high spectral resolution in VNIR (SH-GF5VN), AC can be classified with higher accuracy of 84.00% against 78.00% from SH.
Visually, the misclassification and confusion can be also observed from the mapping (Figure 4). In the south part of study area (area A in Figure 4a), the majority of KO 1 was misclassified into AM, which is consistence with low producer accuracy of 68.08% in KO 1 . In the area B where narrow zonation of AM was between AI 1 and KO 1 , the confusion happens when using Landsat 8 (low user accuracy of 68.04%); while for area C, AI 1 and AM are confused. For the distribution of mangrove species resulting from simulated Hyperion and GF-5 hyperspectral data, they present a similar pattern with a minor difference in area D and E, where a small amount of AM was misclassified into AI 1 and more AC were correctly identified. Similarly, the results from SVM with datasets from Landsat 8 (Figure 4d), SH (Figure 4e) to SH-GF5VN (Figure 4f) show the major changes in AI 1 and AI 2 . More specifically, the classification of AI 1 and AI 2 were refined with fine particles when using hyperspectral data (SH and SH-GF5VN). In contrast, fewer ACs were found from the mangrove species mapping using SVM, and KO located at the seaward side cannot be identified unless SH and SH-GF5VN are used.    Finally, we calculated the areas for each mangrove species in the study area using Random Forests with SH-GF5VN. The results show that the most dominant species is KO 1 with area of 103.95±3.56 ha (46% of the mangrove forests in the study area), followed by AM, AI 2 , AI 1 , and KO 2 . The AC species are rare with an area of 6.93±1.07 ha, occupying only 3% of the mangrove forests in the study area. The details can be found in Table 7.

Quality Assessment of Simulated Hyprion
Simulated Hyperion based on GF-5 hyperspectral data is formulated by subset extraction, which hypothesizes that the deviation from two different sensors is minor for mangrove species mapping. To assess the quality of simulated Hyperion, we collected the spectral characteristics of six mangrove species based on shared samples which were picked up from the stable area without human interference in simulated Hyperion and true Hyperion data in 2008. Consequently, the spectral characteristics from two data sets should be similar with only a difference in growth, which can be viewed as intra-class variation. The correlation analysis of two spectra of six mangrove species indicates that it is reasonable to simulate Hyperion based on GF-5 hyperspectral data ( Figure 5). Except for AC species, the Simulated Hyperion has low correlation with true Hyperion, which is possibly because AC grows at the outmost of mangrove forests where mixed pixels with mudflats are easily encountered. The reflectance of the rest simulated from GF-5 data is highly correlated with that of true Hyperion. Therefore, the Simulated Hyperion can be an alternative for comparison if Hyperion is not available.
Remote Sens. 2020, 12, x FOR PEER REVIEW 11 of 20 indicates that it is reasonable to simulate Hyperion based on GF-5 hyperspectral data (Error! Reference source not found.). Except for AC species, the Simulated Hyperion has low correlation with true Hyperion, which is possibly because AC grows at the outmost of mangrove forests where mixed pixels with mudflats are easily encountered. The reflectance of the rest simulated from GF-5 data is highly correlated with that of true Hyperion. Therefore, the Simulated Hyperion can be an alternative for comparison if Hyperion is not available.

Classifier Selection
According to the aforementioned findings, we learned that for the same input Random Forests outperform SVM if there is no feature selection. The contribution from extra minor information may be masked by abundant information due to high-dimension features with dependency when using SVM, which may explain the reason why further improvement was not expected with SVM but did

Classifier Selection
According to the aforementioned findings, we learned that for the same input Random Forests outperform SVM if there is no feature selection. The contribution from extra minor information may be masked by abundant information due to high-dimension features with dependency when using SVM, which may explain the reason why further improvement was not expected with SVM but did with Random Forests. In another words, SVM may perform better with feature selection, which was also indicated in [39]. It should be noted that no feature selection was conducted in this study, because we did not aim for classifier optimization but for the contribution from the extra spectra on mangrove species mapping, and feature selection may lead to information loss.

Contribution of the Increase in Spectral Resolution in VNIR
We have learned from the confusion matrixes that the advantage of simulated Hyperion over Landsat 8 is to improve the accuracy of identification of many species including KO, AM, and AI, while, in contrast, simulated Hyperion GF-5 hyperspectral data shows metrics on the AC identification. From the Figure 2, we can see that the spectra of different species ranging from 403 nm to 740 nm (where most spectra of Landsat 8 lie) are similar. In contrast, most spectra from 780 nm to 1300 nm show distinguished differences. That is why the hyperspectral data (GF-5 and simulated Hyperion) can outperform Landsat 8 in mangrove species mapping with great improvement. In visible spectra, KO, AI and AC are with slight differences in reflectance, and the reflectance of AM lying between KO and AI are more similar to KO 1 and AI 1 . This may explain the confusion of AM, KO 1 and AI 1 resulting in low accuracies in classification. In VNIR, KO has distinct features from others, and AM can be easily separated from one other with spectra ranging from 750 to 960nm, while AC and AI have similar spectra. This reveals that the increase in spectral resolution (from Landsat 8 to simulated Hyperion) can better discriminate mangrove species, and increasing the spectral resolution for hyperspectral data (from simulated Hyperion to GF-5 hyperspectral data) can better discriminate species with minor differences, which does make sense. The side effect of accuracy drop in KO and AM, which should have been easily separated with distinct features difference, is possibly because the introduced high-resolution spectra enlarge the intra-class variation for them. The similar problem of intra-class enlargement can be seen when increasing the spatial resolution. This may be resolved through individual processing by category, which was also used to exclude inter-class variation first to improve the results of mangrove biomass estimation [44]. This means that the effect of the additional spectra due to increase in spectral resolution should be explored in a new way.

Limitation of the Study
Due to a lack of true Hyperion data for comparison, a simulated Hyperion in this study was generated based on GF-5 hyperspectral in formulation of subset extraction. However, there is not always a solution for Equation (1). For example, the spectra in 1340-1520 nm for Hyperion cannot be simulated due to the differences of invalid bands between two sensors. The spectral loss (see Section 3.1 from simulated Hyperion should be analyzed, although there are only nine bands and they are close to the absorption windows of water vapor. Similarly, the work to which spectral features contribute most for one certain mangrove species should be analyzed. This may be useful for mangrove species discrimination in other sites. Another limitation is the sample number, especially for AC with small patches. More data will make the results more robust and convincing. Due to the small and narrow areas of each mangrove species, the GF-5 hyperspectral data with spatial resolution of 30 m cannot capture many pixels for each species, which makes it hard to collect more samples from the GF-5 hyperspectral data in the study area. Compared to the samples used in [6,14,18], we accepted the repeated use of our samples.

Conclusions
To date, GF-5 as the latest satellite which, in commission with hyperspectral data provision, offers potential opportunities for accurate observation due to higher spectral resolution in VNIR. In our study, we used GF-5 hyperspectral data for accurate mapping of mangroves in the Mai Po at species level. Our aim was to assess the performance of GF-5 hyperspectral data on mangrove species discrimination. We simulated Hyperion for the control hyperspectral data. Together with Landsat 8, various data with different spectral resolution is compared using Random Forests and SVM. According to the aforementioned results, we can draw the conclusions that follow.
1. GF-5 hyperspectral data can be used for accurately mapping six mangrove species in the Mai Po with accuracy of 87.12%; 2. GF-5 shows an advantage over Landsat 8 and Hyperion in mangrove species mapping due to an increase in spectral resolution in VNIR. Therefore, it is recommended to include more specific bands in VNIR for future satellites with the aims of mangrove species discrimination.
3. Extra spectra from GF-5 hyperspectral data with higher spectral resolution in VNIR acts on AC mapping; 4. Random Forests are preferred to SVM for mapping mangrove species when using GF-5 hyperspectral data without feature selection for spectral features.

Acknowledgments:
The authors would like to thank the Shanghai Institute of Technical Physics, CAS for providing the GF-5 AHSI data, and also thank the editor and two anonymous reviewers for their critical comments and suggestions to improve the original manuscript.

Conflicts of Interest:
The authors declare no conflict of interest.