A red edge band is a sensitive spectral band of crops, which helps to improve the accuracy of crop classification. In view of the characteristics of GF-6 WFV data with multiple red edge bands, this paper took Hengshui City, Hebei Province, China, as the study area to carry out red edge feature analysis and crop classification, and analyzed the influence of different red edge features on crop classification. On the basis of GF-6 WFV red edge band spectral analysis, different red edge feature extraction and red edge indices feature importance evaluation, 12 classification schemes were designed based on GF-6 WFV of four bands (only including red, green, blue and near-infrared bands), stepwise discriminant analysis (SDA) and random forest (RF) method were used for feature selection and importance evaluation, and RF classification algorithm was used for crop classification. The results show the following: (1) The red edge 750 band of GF-6 WFV data contains more information content than the red edge 710 band. Compared with the red edge 750 band, the red edge 710 band is more conducive to improving the separability between different crops, which can improve the classification accuracy; (2) According to the classification results of different red edge indices, compared with the SDA method, the RF method is more accurate in the feature importance evaluation; (3) Red edge spectral features, red edge texture features and red edge indices can improve the accuracy of crop classification in different degrees, and the red edge features based on red edge 710 band can improve the accuracy of crop classification more effectively. This study improves the accuracy of remote sensing classification of crops, and can provide reference for the application of GF-6 WFV data and its red edge bands in agricultural remote sensing.

High-intensity coal mining (large mining height, shallow mining depth, and rapid advancing) frequently causes large-scale ground damage within a short period of time. Understanding mining subsidence under high-intensity mining can provide a basis for mining-induced damage assessment, land remediation in a subsidence area, and ecological reconstruction in vulnerable ecological regions in Western China. In this study, the mining subsidence status of Shendong Coalfield was investigated and analyzed using two-pass differential interferometric synthetic aperture radar (DInSAR) technology based on high-resolution synthetic aperture radar data (RADARSAT-2 precise orbit, multilook fine, 5 m) collected from 20 January 2012 to June 2013. Surface damages in Shendong Coalfield over a period of 504 days under open-pit mining and underground mining were observed. Ground deformation of the high-intensity mining working faces 22201-1/2 in Bu’ertai Mine, Shendong Coalfield was monitored using small baseline subset (SBAS) InSAR technology. (1) DInSAR detected and located 85 ground deformation areas (including ground deformations associated with past-mining activity). The extent of subsidence in Shendong Coalfield presented a progressive increase at an average monthly rate of 13.09 km² from the initial 54.98 km² to 225.20 km², approximately, which accounted for 7% of the total area of Shendong Coalfield; (2) SBAS-InSAR reported that the maximum cumulative subsidence area reached 5.58 km² above the working faces 22201-1/2. The advance speed of ground destruction (7.9 m/day) was nearly equal to that of underground mining (8.1 m/day).

Time series of vegetation indices can be utilized to capture crop phenology information, and have been widely used in land cover and crop classification, phenological feature extraction, and planting structure monitoring. This is of great significance for guiding agricultural production and formulating agricultural policies. According to the characteristics of the GF-6 satellite’s newly-added red edge bands, wide field view and high-frequency imaging, the time series of vegetation indices about multi-temporal GF-6 WFV data are used for the study of land cover and crop classification. In this study, eight time steps of GF-6 WFV data were selected from March to October 2019 in Hengshui City. The normalized difference vegetation index (NDVI) time series and 10 different red edge spectral indices time series were constructed. Then, based on principal component analysis (PCA), using two feature selection and evaluation methods, stepwise discriminant analysis (SDA) and random forest (RF), the red edge vegetation index of normalized difference red edge (NDRE) was selected. Seven different lengths of NDVI, NDRE and NDVI&NDRE time series were reconstructed by the Savizky-Golay (S-G) smoothing algorithm. Finally, an RF classification algorithm was used to analyze the influence of time series length and red edge indices features on land cover and crop classification, and the planting structure and distribution of crops in the study area were obtained. The results show that: (1) Compared with the NDRE red edge time series, the NDVI time series is more conducive to the improvement of the overall classification accuracy of crops, and NDRE can assist NDVI in improving the crop classification accuracy; (2) With the shortening of NDVI and NDRE time series, the accuracy of crop classification is gradually decreased, and the decline is gradually accelerated; and (3) Through the combination of the NDVI and NDRE time series, the accuracy of crop classification with different time series lengths can be improved compared with the single NDVI time series, which is conducive to improving the classification accuracy and timeliness of crops. This study has fully tapped the application potential of the new red edge bands of GF-6 WFV time series data, which can provide references for crop identification and classification of time series data such as NDVI and red edge vegetation index of different lengths. At the same time, it promotes the application of optical satellite data with red edge bands in the field of agricultural remote sensing.

An understanding of the response of interannual vegetation variations to climate change is critical for the future projection of ecosystem processes and developing effective coping strategies. In this study, the spatial pattern of interannual variability in the growing season normalized difference vegetation index (NDVI) for different biomes and its relationships with climate variables were investigated in Inner Mongolia during 1982–2015 by jointly using linear regression, geographical detector, and geographically weighted regression methodologies. The result showed that the greatest variability of the growing season NDVI occurred in typical steppe and desert steppe, with forest and desert most stable. The interannual variability of NDVI differed monthly among biomes, showing a time gradient of the largest variation from northeast to southwest. NDVI interannual variability was significantly related to that of the corresponding temperature and precipitation for each biome, characterized by an obvious spatial heterogeneity and time lag effect marked in the later period of the growing season. Additionally, the large slope of NDVI variation to temperature for desert implied that desert tended to amplify temperature variations, whereas other biomes displayed a capacity to buffer climate fluctuations. These findings highlight the relationships between vegetation variability and climate variability, which could be used to support the adaptive management of vegetation resources in the context of climate change.

High-resolution crop type mapping is of importance for site-specific agricultural management and food security in smallholder farming regions, but is challenging due to limited data availability and the need for image-based algorithms. In this paper, we developed an efficient object- and pixel-based mapping algorithm to generate a 10 m resolution crop type map over large spatial domains by integrating time series optical images (Sentinel-2) and synthetic aperture radar (SAR) images (Sentinel-1) using the Google Earth Engine (GEE) platform. The results showed that the proposed method was reliable for crop type mapping in the study area with an overall accuracy (OA) of 93.22% and a kappa coefficient (KC) of 0.89. Through experiments, we also found that the monthly median values of the vertical transmit/vertical receive (VV) and vertical transmit/horizontal receive (VH) bands were insensitive to crop type mapping itself, but adding this information to supplement the optical images improved the classification accuracy, with an OA increase of 0.09–2.98%. Adding the slope of vegetation index change (VI_slope) at the critical period to crop type classification was obviously better than that of relative change ratio of vegetation index (VI_ratio), both of which could make an OA improvement of 2.58%. These findings not only highlighted the potential of the VI_slope and VI_ratio indices during the critical period for crop type mapping in small plots, but suggested that SAR images could be included to supplement optical images for crop type classification.

The geological environment damage caused by coal mining subsidence has become an important factor affecting the sustainable development of mining areas. Reconstruction of the Coal Mining Subsidence Field (CMSF) is the key to preventing geological disasters, and the needs of CMSF reconstruction cannot be met by solely relying on a single remote sensing technology. The combination of Unmanned Aerial Vehicle (UAV) and Synthetic Aperture Radar (SAR) has complementary advantages; however, the data fusion strategy by refining the SAR deformation field through UAV still needs to be updated constantly. This paper proposed a Prior Weighting (PW) method based on Satellite Aerial (SA) heterogeneous remote sensing. The method can be used to fuse SAR and UAV Light Detection and Ranging (LiDAR) data for ground subsidence parameter inversion. Firstly, the subsidence boundary of Differential Interferometric SAR (DInSAR) combined with the large gradient subsidence of Pixel Offset Tracking (POT) was developed to initialize the SAR preliminary CMSF. Secondly, the SAR preliminary CMSF was refined by UAV LiDAR data; the weights of SAR and UAV LiDAR data are 0.4 and 0.6 iteratively. After the data fusion, the subsidence field was reconstructed. The results showed that the overall CMSF accuracy improved from ±144 mm to ±51 mm. The relative errors of the surface subsidence factor and main influence angle tangent calculated by the physical model and in situ measured data are 1.3% and 1.7%. It shows that the proposed SAR/UAV fusion method has significant advantages in the reconstruction of CMSF, and the PW method contributes to the prevention and control of mining subsidence.