Optical remote sensing indices play an important role in vegetation information extraction and have been widely serving ecology, agriculture and forestry, urban monitoring, and other communities. Remote sensing indices are constructed from individual bands depending on special characteristics to enhance the typical spectral features for the identification or distinction of surface land covers. With the development of quantitative remote sensing, there is a rapid increasing requirement for accurate data processing and modeling. It is well known that the geometry-induced variation observed in surface reflectance is not ignorable, but the situation of uncertainty thereby introduced into these indices still needs further detailed understanding. We adopted the ground multi-angle hyperspectrum, spectral response function (SRF) of Thematic Mapper (TM), Enhanced Thematic Mapper (ETM+), Operational Land Imager (OLI), Moderate-Resolution Imaging Spectroradiometer (MODIS), and Multi-Spectral Instrument (MSI) optical sensors and simulated their sensor-like spectral reflectance; then, we investigated the potential angle effect uncertainty on optical indices that have been frequently involved in vegetation monitoring and examined the forward/backward effect over both the ground-based level and the actual Landsat TM/ETM+ overlapped region. Our results on the discussed indices and sensors show as following: (1) Identifiable angle effects exist with a more elevated influence than that introduced by band difference among sensors; (2) The absolute difference between forward and backward direction can reach up to −0.03 to 0.1 within bands of the TM/ETM+ overlapped region; (3) The investigation at ground level indicates that there are different variations of angle effect transmitted to each remote sensing index. Regarding cases of crop canopy at various growth phases, most of the discussed indices have more than a 20% relative difference to nadir value except Normalized Difference Vegetation Index (NDVI) and Normalized Difference Water Index (NDWI) with the magnitude lower than 10%, and less than 16% of Normalized Burn Ratio (NBR). For the case of wax maturity stage, the relative difference to nadir value of Enhanced Vegetation Index (EVI), Soil-Adjusted Vegetation Index (SAVI), Ratio Vegetation Index (RVI), Char Soil Index (CSI), NBR, Normalized Difference Moisture Index (NDMI), and SWIR2/NIR exceeded 50%, while the values for NBR and NDMI can reach up to 115.8% and 206.7%, respectively; (4) Various schemes of index construction imply different propagation of angle effect uncertainty. The “difference” indices can partially suppress the directional influence, while the “ratio” indices show high potential to amplify the angle effect. This study reveals that the angle-induced uncertainty of these indices is greater than that induced by the spectrum mismatch among sensors, especially under the case of senescence. In addition, based on this work, indices with a suppressed potential of angle effect are recommended for vegetation monitoring or information retrieval to avoid unexpected effects.

Automatic determination of quantitative parameters describing the pattern of urban expansion is extremely important for urban planning, urban management and civic resource configuration. Though the widely adopted LEI (landscape expansion index) has exhibited the potential to capture the evolution of urban landscape patterns using multi-temporal remote sensing data, erroneous determination still exists, especially for patches with special shapes due to the limited consideration of spatial relationships among patches. In this paper, we improve the identification of urban landscape expansion patterns with an enhancement of the topological relationship. We propose MCI (Mean patch Compactness Index) and AWCI (Area-Weighted Compactness Index) in terms of the moment of inertia shape index. The effectiveness of the improved approach in identifying types of expansion patches is theoretically demonstrated with a series of designed experiments. Further, we apply the proposed approaches to the analysis of urban expansion features and dynamics of urban compactness over Urumqi at various 5-year stages using available SUBAD-China data from 1990–2015. The results achieved by the theoretical experiments and case application show our approach effectively suppressed the effects induced by shapes of patches and buffer or envelope box parameters for the accurate identification of patch type. Moreover, the modified MCI and AWCI exhibited an improved potential in capturing the landscape scale compactness of urban dynamics. The investigated 25-year urban expansion of Urumqi is dominated by edge-expansion patches and supplemented by outlying growth, with opposite trends of increasing and decreasing, with a gradual decrease in landscape fragmentation. Our examination using the proposed MCI and AWCI indicates Urumqi was growing more compact in latter 15-year period compared with the first 10 years studied, with the primary urban patches tending to be compacted earlier than the entire urban setting. The historical transformation trajectories based on remote sensing data show a significant construction land gain—from 1.06% in 1990 to 6.96% in 2015—due to 289.16 km² of recently established construction, accompanied by fast expansion northward, less dynamic expansion southward, and earlier extension in the westward direction than eastward. This work provides a possible means to improve the identification of patch expansion type and further understand the compactness of urban dynamics.

The spatial distribution of cotton fields is primary information for national farm management, the agricultural economy and the textile industry. Therefore, accurate cotton information at the regional scale is required with a rapid increase due to the chance provided by the huge amounts of satellite images accumulated in recent decades. Research has started to introduce the phenology characteristics shown at special growth phases of cotton but frequently focuses on limited vegetation indices with less consideration on the whole growth period. In this paper, we investigated a set of phenological and time-series features with optimization depending on each feature permutation’s importance and redundancy, followed by its performance evaluation through the cotton extraction using the Random Forest (RF) classifier. Three sets of 31 features are involved: (1) phenological features were determined by the biophysical and biochemical characteristics in the spectral space of cotton during each of its five distinctive phenological stages, which were identified from 2307 representative cotton samples using 21,237 Sentinel-2 images; (2) three typical vegetation indices were functionalized into time-series features by harmonic analysis; (3) three terrain factors were derived from the digital elevation model. Our analysis of feature determination revealed that the most valuable discriminators for cotton involve the boll opening stage and harmonic coefficients. Moreover, both qualitative and quantitative validation were performed to evaluate the retrieval of the optimized features-based cotton information. Visual examination of the map exhibited high spatial consistency and accurate delineation of the cotton field. Quantitative comparison indicates that classification of RF-coupled optimized features achieves improved overall accuracy 5.53% higher than that which works with either the limited vegetation indices. Compared with all 31 features, the optimized features realized greater identification accuracy while using only about half the number of features. Compared with test samples, the cotton map achieved an overall accuracy greater than 98% and a kappa more than 0.96. Further comparison of the cotton map area at the county-level showed a high level of consistency with the National Bureau of Statistics data from 2020, with R² over 0.96, RMSE no more than 14.62 Kha and RRMSE less than 17.78%.

Floral phenology as a special indicator of climate change and vegetation dynamics is drawing more attention. The long-term observations of flowering events collected at scattered ground sites have accumulated valuable priority on the understanding of floral phenology, but with insufficient investigation on the spatio-temporal dynamics at regional scale, which is mainly induced by the lack of effective ways to capture the pixel-based flower events from remote sensing images. The existing yellowness indices are constructed for rape (Brassica napus L.) with less suppression to the bright background and dark green vegetation, and further with inadequate consideration on physiological characteristics and the temporal spectral signature of investigated vegetation. In this paper, we examined rape and several other representative vegetation types to determine spectral features of yellow-flower period within the growing season, then selected the visible and near-infrared bands to construct a Novel Yellowness Index (NYI) with an enhancement on the physiological mechanism of plants. The proposed NYI were discussed on the variation of mathematical properties with representative instances, cross-compared with three typical yellowness indices—Ratio Yellowness Index (RYI), Normalized Difference Yellowness Index (NDYI), and Ashourloo Canola Index (ACI) —over various yellow-flowering vegetation species at multiple scales, and validated with ground observations of three available PhenoCam network stations and field phenological observations at Görlitz, Sachsen, and Germany. In addition, we applied NYI to detect the rape field using Sentinel-2 image at Görlitz with typical rape area as a case study. Results show that the proposed NYI exhibits the potential to capture yellow-flowering events with increased sensitivity to the variation of flower density, and reduction of noise introduced by bright background or dark green vegetation of multiple vegetation species at different scales. As the flower density increases from 33% to 78%, the relative differences of NYI captured can reach up to 74%, compared with other three indices which have the relative differences no more than 57%. The cross-comparison indicates NYI performs better with higher consistent with PhenoCam observation and Deutscher Wetterdienst phenological station than other yellowness indices in capturing the variation of yellow flower density. The case study of NYI application in the identification of rape field exhibits good accuracy with the overall accuracy up to 97.5%, the Kappa coefficient of 0.94, and F score of 0.96. Consequently, the satellite-derived yellowness index will be a potential means to investigate the flowering dynamics and planting range of yellow-flowering vegetation such as rape.

Surface albedo, as an important parameter for land surface geo-biophysical and geo-biochemical processes, has been widely used in the research communities involved in surface energy balance, weather forecasting, atmospheric circulation, and land surface process models. In recent years, operational products using satellite-based surface albedo have, from time to time, been rapidly developed, contributing significantly to the estimation of energy balance at regional or global scales. The increasing number of research topics on dynamic monitoring at a decades-long scale requires a combination of albedo products generated from various sensors or programs, while the quantitative assessment of agreement or divergence among different surface albedo products still needs further understanding. In this paper, we investigated the consistency of three classical operational surface albedo products that have been frequently used by researchers globally via the official issued datasets-MODIS, GLASS (Global LAnd Surface Satellite), and CGLS (Copernicus Global Land Service). The cross-comparison was performed on all the identical dates available during 2000–2017 to represent four season-phases. We investigated the pixel-based validity of each product, consistency of global annual mean, spatial distribution and different temporal dynamics among the discussed products in white-sky (WSA) and black-sky (BSA) albedo at visible (VIS), near-infrared (NIR), and shortwave (SW) regimes. Further, varying features along with the change of seasons was also examined. In addition, the variation in accuracy of shortwave albedo magnitude was explored using ground measurements collected by the Baseline Surface Radiation Network (BSRN) and the Surface Radiation Budget Network (SUFRAD). Results show that: (1) All three products can provide valid long-term albedo for dominant land surface, while GLASS can provide additional estimation over sea surfaces, with the highest percentage of valid land surface pixels, at up to 93% in 24 October. The invalid pixels mainly existed in the 50°N–60°N latitude belt in December for GLASS, Central Africa in April and August for MODIS, and northern high latitudes for CGLS. (2) The global mean albedo of CGLS at the investigated bands has significantly higher values than those of MODIS and GLASS, with a relative difference of ~20% among the three products. The global mean albedo of MODIS and GLASS show a generally increasing trend from April to December, with an abrupt rise at NIR and SW of CGLS in June of 2014. Compared with SW and VIS bands, the linear temporal trend of the NIR global albedo mean in three products continues to increase, but the slope of CGLS is 10–100 times greater than that of the other two products. (3) The differences in albedo, which are higher in April, October, and December than in August, exhibit a small variation over the main global land surface regions, except for Central Eurasia, North Africa, and middle North America. The magnitude of global absolute difference among the three products usually varies within 0.02–0.06, but with the largest value occasionally exceeding 0.1. The relative difference is mainly within 10–20%, and can deviate more than 40% away from the baseline. In addition, CGLS has a greater opportunity to achieve the largest difference compared with MODIS and GLASS. (4) The comparison with ground measurements indicates that MODIS generally performs better than GLASS and CGLS at the sites discussed. This study demonstrates that apparent differences exist among the three investigated albedo products due to the ingested source data, algorithm, atmosphere correction etc., and also points at caution regarding data fusion when multiple albedo products were organized to serve the following applications.

Offshore drilling rigs are the foundation of oil and gas exploitation in water areas. Their spatial and temporal distribution, state attributes and other information directly reflect the exploitation level of regional oil and gas resources. Therefore, it is very important to build an automatic detecting method for offshore drilling rigs with good performance to accurately capture the temporal and spatial distribution and state of oil and gas exploitation activities. At present, there are two main groups of methods for offshore drilling rigs: invariant feature-based methods and nighttime firelight-based methods. Methods based on invariant location are more subjective in terms of their parameter settings and require intensive computation. Nighttime light-based methods are largely unable to identify offshore drilling rigs without associated waste gas ignition. Furthermore, multiple offshore drilling rigs in close proximity to one another cannot be effectively distinguished with low spatial resolution imagery. To address these shortcomings, we propose a new method for the automatic identification of offshore drilling rigs based on Landsat-7 ETM+ images from 2018 to 2019, taking the Caspian Sea as the research area. We build a nominal annual cloud and cloud shadow-free Normalized Difference Water Index (NDWI) composite by designing an optimal NDWI compositing method based of the influence of cloud and cloud shadow on the NDWI values of water, bare land (island) and offshore drilling rigs. The classification of these objects is simultaneously done during the compositing process, with the following rules: water body (Max_NDWI > 0.55), bare land (island) (Min_NDWI < −0.05) and offshore drilling rig (0 < Mean_NDWI < 0.4). A threshold segmentation and postprocessing were carried out to further refine the results. Using this method, 497 offshore platforms were automatically identified using a nominal annual cloud and cloud shadow-free NDWI composite image and Google Earth Engine. Validation using Sentinel-2 Multispectral Imager (MSI) and Google Earth images demonstrated that the correct rate of offshore drilling rig detection in the Caspian Sea is 90.2%, the missing judgment rate is 5.3% and the wrong judgment rate is 4.5%, proving the performance of the proposed method. This method can be used to identify offshore drilling rigs within a large water surface area relatively quickly, which is of great significance for exploring the exploitation status of offshore oil and gas resources. It can also be extended to finer spatial resolution optical remote sensing images; thus small-size drilling rigs can be effectively detected.