Relative Radiometric Normalization Research Articles

LIRRN: Location-Independent Relative Radiometric Normalization of Bitemporal Remote-Sensing Images.

Relative radiometric normalization (RRN) is a critical pre-processing step that enables accurate comparisons of multitemporal remote-sensing (RS) images through unsupervised change detection. Although existing RRN methods generally have promising results in most cases, their effectiveness depends on specific conditions, especially in scenarios with land cover/land use (LULC) in image pairs in different locations. These methods often overlook these complexities, potentially introducing biases to RRN results, mainly because of the use of spatially aligned pseudo-invariant features (PIFs) for modeling. To address this, we introduce a location-independent RRN (LIRRN) method in this study that can automatically identify non-spatially matched PIFs based on brightness characteristics. Additionally, as a fast and coregistration-free model, LIRRN complements keypoint-based RRN for more accurate results in applications where coregistration is crucial. The LIRRN process starts with segmenting reference and subject images into dark, gray, and bright zones using the multi-Otsu threshold technique. PIFs are then efficiently extracted from each zone using nearest-distance-based image content matching without any spatial constraints. These PIFs construct a linear model during subject-image calibration on a band-by-band basis. The performance evaluation involved tests on five registered/unregistered bitemporal satellite images, comparing results from three conventional methods: histogram matching (HM), blockwise KAZE, and keypoint-based RRN algorithms. Experimental results consistently demonstrated LIRRN's superior performance, particularly in handling unregistered datasets. LIRRN also exhibited faster execution times than blockwise KAZE and keypoint-based approaches while yielding results comparable to those of HM in estimating normalization coefficients. Combining LIRRN and keypoint-based RRN models resulted in even more accurate and reliable results, albeit with a slight lengthening of the computational time. To investigate and further develop LIRRN, its code, and some sample datasets are available at link in Data Availability Statement.

Radiometric Normalization Using a Pseudo−Invariant Polygon Features−Based Algorithm with Contemporaneous Sentinel−2A and Landsat−8 OLI Imagery

As sensor parameters and atmospheric conditions create uncertainties for at−sensor radiation detection, radiometric consistency among satellite images is difficult to achieve. Relative radiometric normalization is a method that can improve multi−image consistency with accurate pseudo−invariant features (PIFs), especially over large areas or long time series satellite images. Although there are algorithms that manually or automatically select PIFs, the spatial mismatch of satellite images can affect PIF extraction, particularly with artificial pixels. To alleviate this problem, we proposed to use Landsat−8 OLI as the reference image and Sentinel−2A as the subject image, to apply pseudo−invariant features−based algorithms with polygon features through the single−band and multiple−band regression. Compared to pseudo−invariant point features, hyperspectral library, and histogram matching approaches, the results demonstrate the superiority of the proposed algorithms with correlation coefficients of 0.9948 and 0.9945, and an RMSE of 0.0097 and 0.0095 with multiple− and single−band regression, respectively. We also found more accurate linear fitting and better shape matching through band scattering and reflectance frequency analysis. The proposed algorithms are a significant improvement in radiometric normalization, within artificial pixels, achieving spectral signature consistency.

Robust Radiometric Normalization of the Near Equatorial Satellite Images Using Feature Extraction and Remote Sensing Analysis

Relative radiometric normalization (RRN) minimizes radiometric differences among images caused by inconsistencies of acquisition conditions rather than changes in surface. Scale invariant feature transform (SIFT) has the ability to automatically extract control points (CPs) and is commonly used for remote sensing images. However, its results are mostly inaccurate and sometimes contain incorrect matching caused by generating a small number of false CP pairs. These CP pairs have high false alarm matching. This paper presents a modified method to improve the performance of SIFT CPs matching by applying sum of absolute difference (SAD) in a different manner for the new optical satellite generation called near-equatorial orbit satellite and multi-sensor images. The proposed method, which has a significantly high rate of correct matches, improves CP matching. The data in this study were obtained from the RazakSAT satellite a new near equatorial satellite system. The proposed method involves six steps: 1) data reduction, 2) applying the SIFT to automatically extract CPs, 3) refining CPs matching by using SAD algorithm with empirical threshold, and 4) calculation of true CPs intensity values over all image’ bands, 5) preforming a linear regression model between the intensity values of CPs locate in reverence and sensed image’ bands, 6) Relative radiometric normalization conducting using regression transformation functions. Different thresholds have experimentally tested and used in conducting this study (50 and 70), by followed the proposed method, and it removed the false extracted SIFT CPs to be from 775, 1125, 883, 804, 883 and 681 false pairs to 342, 424, 547, 706, 547, and 469 corrected and matched pairs, respectively.

A Multirule-Based Relative Radiometric Normalization for Multisensor Satellite Images

A Multirule-Based Relative Radiometric Normalization for Multisensor Satellite Images

Tensor-based keypoint detection and switching regression model for relative radiometric normalization of bitemporal multispectral images

ABSTRACT In some remote sensing applications, such as unsupervised change detection, bitemporal multispectral images must be first aligned/harmonized radiometrically. For doing so, Many Relative Radiometric Normalization (RRN) algorithms exist; however, most suffer from misregistration problems and can only operate on geo/co-registered image pairs, while unregistered multispectral pairs are required. To tackle this situation, keypoint-based RRN methods were introduced, which can radiometrically calibrate unregistered/registered image pairs using keypoint matching algorithms. However, they ignore the spatial and spectral characteristics of spectral bands of input images, resulting in potential RRN errors. They also employ a linear mapping function for RRN modelling, which can not handle non-linear radiometric distortions. To address these limitations, this paper proposes a robust algorithm for RRN of bitemporal multispectral images, using a new extension of SURF detector for multispectral images, namely the Weighted Spectral Structure Tensor SURF (WSST-SURF), and a flexible Switching Regression (SR) model. Taking advantage of the tensor theory, WSST-SURF efficiently preserves both spatial and spectral information distributed over all bands of multispectral images for keypoint detection, resulting in extracting reliable inliers (or keypoints) for RRN. An adaptive SR model is introduced based on the normalized mutual information, accurately approximating a linear/non-linear relationship between inliers in multispectral images. Six unregistered multispectral image pairs captured by inter/intra remote sensing sensors were employed to validate the efficacy of the proposed method. The results indicate that adopting spectral tensor-based SURF methods in the RRN process exhibits better local and global performance than using the original SURF. Furthermore, the proposed method outperforms the existing conventional RRN methods in terms of accuracy and visual quality, indicating its competence for RRN of bitemporal multispectral images with high illumination, viewpoint, and scale differences.

GISD30: global 30 m impervious-surface dynamic dataset from 1985 to 2020 using time-series Landsat imagery on the Google Earth Engine platform

Abstract. Accurately mapping impervious-surface dynamics has great scientific significance and application value for research on urban sustainable development, the assessment of anthropogenic carbon emissions and global ecological-environment modeling. In this study, a novel and automatic method of combining the advantages of spectral-generalization and automatic-sample-extraction strategies was proposed, and then an accurate global 30 m impervious-surface dynamic dataset (GISD30) for 1985 to 2020 was produced using time-series Landsat imagery on the Google Earth Engine cloud computing platform. Firstly, the global training samples and corresponding reflectance spectra were automatically derived from prior global 30 m land-cover products after employing the multitemporal compositing method and relative radiometric normalization. Then, spatiotemporal adaptive classification models, trained with the migrated reflectance spectra of impervious surfaces from 2020 and transferred pervious-surface samples in each epoch for every 5∘×5∘ geographical tile, were applied to map the impervious surface in each period. Furthermore, a spatiotemporal-consistency correction method was presented to minimize the effects of independent classification errors and improve the spatiotemporal consistency of impervious-surface dynamics. Our global 30 m impervious-surface dynamic model achieved an overall accuracy of 90.1 % and a kappa coefficient of 0.865 using 23 322 global time-series validation samples. Cross-comparisons with five existing global 30 m impervious-surface products further indicated that our GISD30 dynamic product achieved the best performance in capturing the spatial distributions and spatiotemporal dynamics of impervious surfaces in various impervious landscapes. The statistical results indicated that the global impervious surface has doubled in the past 35 years, from 5.116×105 km2 in 1985 to 10.871×105 km2 in 2020, and Asia saw the largest increase in impervious surface area compared to other continents, with a total increase of 2.946×105 km2. Therefore, it was concluded that our global 30 m impervious-surface dynamic dataset is an accurate and promising product and could provide vital support in monitoring regional or global urbanization as well as in related applications. The global 30 m impervious-surface dynamic dataset from 1985 to 2020 generated in this paper is free to access at https://doi.org/10.5281/zenodo.5220816 (Liu et al., 2021b).

Automatic Relative Radiometric Normalization of Bi-Temporal Satellite Images Using a Coarse-to-Fine Pseudo-Invariant Features Selection and Fuzzy Integral Fusion Strategies

Relative radiometric normalization (RRN) is important for pre-processing and analyzing multitemporal remote sensing (RS) images. Multitemporal RS images usually include different land use/land cover (LULC) types; therefore, considering an identical linear relationship during RRN modeling may result in potential errors in the RRN results. To resolve this issue, we proposed a new automatic RRN technique that efficiently selects the clustered pseudo-invariant features (PIFs) through a coarse-to-fine strategy and uses them in a fusion-based RRN modeling approach. In the coarse stage, an efficient difference index was first generated from the down-sampled reference and target images by combining the spectral correlation, spectral angle mapper (SAM), and Chebyshev distance. This index was then categorized into three groups of changed, unchanged, and uncertain classes using a fast multiple thresholding technique. In the fine stage, the subject image was first segmented into different clusters by the histogram-based fuzzy c-means (HFCM) algorithm. The optimal PIFs were then selected from unchanged and uncertain regions using each cluster’s bivariate joint distribution analysis. In the RRN modeling step, two normalized subject images were first produced using the robust linear regression (RLR) and cluster-wise-RLR (CRLR) methods based on the clustered PIFs. Finally, the normalized images were fused using the Choquet fuzzy integral fusion strategy for overwhelming the discontinuity between clusters in the final results and keeping the radiometric rectification optimal. Several experiments were implemented on four different bi-temporal satellite images and a simulated dataset to demonstrate the efficiency of the proposed method. The results showed that the proposed method yielded superior RRN results and outperformed other considered well-known RRN algorithms in terms of both accuracy level and execution time.

Relative radiometric normalization for mosaicking IRS CartoSat-2 panchromatic images using genetic algorithm

Urban features can be analysed using high-resolution satellite images. Satellites procure panchromatic images at a spatial resolution even better than 1 m. However, it cannot cover the large area of a city in one scene. Therefore, satellite images from different dates and times must be mosaicked to study the urban region as a whole. Efficient and seamless mosaicking requires relative radiometric normalization (RRN), which generally is an iterative process. In this paper, an approach using the genetic algorithm has been developed for normalization in a single step. The algorithm has been tested on Jaipur city in India. Results produced for the mosaic using genetic algorithm provide visually seamless mosaic over the study region with significant improvement of the RMSE among the images’ overlap region. The algorithm has been compared against the mean-standard deviation (MS) and histogram matching (HM) methods of RRN with average RMSE improvement of 40% and 38%, respectively.

Integrated Preprocessing of Multitemporal Very-High-Resolution Satellite Images via Conjugate Points-Based Pseudo-Invariant Feature Extraction

Multitemporal very-high-resolution (VHR) satellite images are used as core data in the field of remote sensing because they express the topography and features of the region of interest in detail. However, geometric misalignment and radiometric dissimilarity occur when acquiring multitemporal VHR satellite images owing to external environmental factors, and these errors cause various inaccuracies, thereby hindering the effective use of multitemporal VHR satellite images. Such errors can be minimized by applying preprocessing methods such as image registration and relative radiometric normalization (RRN). However, as the data used in image registration and RRN differ, data consistency and computational efficiency are impaired, particularly when processing large amounts of data, such as a large volume of multitemporal VHR satellite images. To resolve these issues, we proposed an integrated preprocessing method by extracting pseudo-invariant features (PIFs), used for RRN, based on the conjugate points (CPs) extracted for image registration. To this end, the image registration was performed using CPs extracted using the speeded-up robust feature algorithm. Then, PIFs were extracted based on the CPs by removing vegetation areas followed by application of the region growing algorithm. Experiments were conducted on two sites constructed under different acquisition conditions to confirm the robustness of the proposed method. Various analyses based on visual and quantitative evaluation of the experimental results were performed from geometric and radiometric perspectives. The results evidence the successful integration of the image registration and RRN preprocessing steps by achieving a reasonable and stable performance.

Radiometric Normalization for Cross-Sensor Optical Gaofen Images with Change Detection and Chi-Square Test

As the number of cross-sensor images increases continuously, the surface reflectance of these images is inconsistent at the same ground objects due to different revisit periods and swaths. The surface reflectance consistency between cross-sensor images determines the accuracy of change detection, classification, and land surface parameter inversion, which is the most widespread application. We proposed a relative radiometric normalization (RRN) method to improve the surface reflectance consistency based on the change detection and chi-square test. The main contribution was that a novel chi-square test automatically extracts the stably unchanged samples between the reference and subject images from the unchanged regions detected by the change-detection method. We used the cross-senor optical images of Gaofen-1 and Gaofen-2 to test this method and four metrics to quantitatively evaluate the RRN performance, including the Root Mean Square Error (RMSE), spectral angle cosine, structural similarity, and CIEDE2000 color difference. Four metrics demonstrate the effectiveness of our proposed RRN method, especially the reduced percentage of RMSE after normalization was more than 80%. Comparing the radiometric differences of five ground features, the surface reflectance curve of two Gaofen images showed more minor differences after normalization, and the RMSE was smaller than 50 with the reduced percentages of about 50–80%. Moreover, the unchanged feature regions are detected by the change-detection method from the bitemporal Sentinel-2 images, which can be used for RRN without detecting changes in subject images. In addition, extracting samples with the chi-square test can effectively improve the surface reflectance consistency.

A novel automatic method on pseudo-invariant features extraction for enhancing the relative radiometric normalization of high-resolution images

ABSTRACT Relative radiometric normalization (RRN) is a critical preprocessing step that is widely applied to remote sensing data. Essential to RRN are the pseudo-invariant features (PIFs) – ground objects with invariant reflectance across the acquisition dates of multiple images. To meet the requirements of RRN for change detection in multi-temporal images, this paper proposed a novel automatic method for PIF extraction: first, we calculated spectral, spatial, and visual features for the images. Next, a set of change vectors were generated to quantify the features’ differences. Then, the initial PIF was obtained by synthesizing similar pixels from each change vector. Finally, an iterative linear regression was performed to generate the final PIFs. Furthermore, we employed the random sample consensus (RANSAC) algorithm to establish the RRN model for normalizing the target image. The methods and workflow of PIF extraction and the RRN model establishment were validated with two TripleSat-2 multi-spectral images with a spatial resolution of 3.2 metres. The applicability of the proposed PIF extraction method was further validated by Landsat-8 OLI, GF-1 WFV, Sentinel-2 MSI, and WorldView-2 images. The results showed that: (1) The proposed PIF extraction method can automatically obtain the PIFs with high precision. Based on the extracted PIFs, the RANSAC algorithm is capable of building the RRN model with a significant linear relationship. (2) The PIF extraction method is applicable to multi-spectral images with different spatial resolutions. (3) Compared with the multivariate alteration detection (MAD) and the iteratively reweighted MAD (IR-MAD), the proposed PIF extraction method achieved a better performance in selecting the pixels with smaller radiometric differences, which is helpful to reach a higher accuracy of the RRN model. (4) Differences in image brightness, spectral domain, and spatial domain are considered in the PIF extraction method with fewer parameters and high portability. (5) Spectral features are dominant to obtain ideal PIFs, and feature combinations can get the advantages to enhance the effectiveness of the PIFs.

Distortion Robust Relative Radiometric Normalization of Multitemporal and Multisensor Remote Sensing Images Using Image Features

In this article, we propose a novel framework to radiometrically correct unregistered multisensor image pairs based on the extracted feature points with the KAZE detector and the conditional probability (CP) process in the linear model fitting. In this method, the scale, rotation, and illumination invariant radiometric control set samples (SRII-RCSS) are first extracted by the blockwise KAZE strategy. They are then distributed uniformly over both textured and texture-less land use/land cover (LULC) using grid interpolation and a set of nearest-neighbors. Subsequently, SRII-RCSS are scored by a similarity measure, and the histogram of the scores is then used to refine SRII-RCSS. The normalized subject image is produced by adjusting the subject image to the reference image using the CP-based linear regression (CPLR) based on the optimal SRII-RCSS. The registered normalized image is finally generated by registration of the normalized subject image to the reference image through a two-pass registration method, namely affine-B-spline and, then, it is enhanced by updating the normalization coefficient of CPLR based on the SRII-RCSS. In this study, eight multitemporal data sets acquired by inter/intra satellite sensors were used in tests to comprehensively assess the efficiency of the proposed method. Experimental results show that the proposed method outperforms the existing state-of-the-art relative radiometric normalization (RRN) methods both qualitatively and quantitatively, indicating its capability for RRN of unregistered multisensor image pairs.

Radiometric Performance Evaluation of FY-4A/AGRI Based on Aqua/MODIS.

Fengyun-4A (FY-4A) is the first satellite of the Chinese second-generation geostationary orbit meteorological satellites (FY-4). The Advanced Geostationary Radiation Imager (AGRI), onboard FY-4A does not load with high-precision calibration facility in visible and near infrared (VNIR) channel. As a consequence, it is necessary to comprehensively evaluate its radiometric performance and quantitatively describe the attenuation while using its VNIR data. In this paper, the radiometric performance at VNIR channels of FY-4A/AGRI is evaluated based on Aqua/MODIS data using the deep convective cloud (DCC) target. In order to reduce the influence of view angle and spectral response difference, the bi-directional reflectance distribution function (BRDF) correction and spectral matching have been performed. The evaluation result shows the radiometric performance of FY-4A/AGRI: (1) is less stable and with obvious fluctuations; (2) has a lower radiation level because of 24.99% lower compared with Aqua/MODIS; 3) has a high attenuation with 9.11% total attenuation over 2 years and 4.0% average annual attenuation rate. After the evaluation, relative radiometric normalization between AGRI and MODIS in VNIR channel is performed and the procedure is proved effective. This paper proposed a more reliable reference for the quantitative applications of FY-4A data.

A Nonlinear Radiometric Normalization Model for Satellite Imgaes Time Series Based on Artificial Neural Networks and Greedy Algroithm

Satellite Image Time Series (SITS) is a data set that includes satellite images across several years with a high acquisition rate. Radiometric normalization is a fundamental and important preprocessing method for remote sensing applications using SITS due to the radiometric distortion caused by noise between images. Normalizing the subject image based on the reference image is a general strategy when using traditional radiometric normalization methods to normalize multi-temporal imagery (usually two or three scenes in different time phases). However, these methods are unsuitable for calibrating SITS because they cannot minimize the radiometric distortion between any pair of images in SITS. The existing relative radiometric normalization methods for SITS are based on linear assumptions, which cannot effectively reduce nonlinear radiometric distortion caused by continuously changing noise in SITS. To overcome this problem and obtain a more accurate SITS, we propose a nonlinear radiometric normalization model (NMAG) for SITS based on Artificial Neural Networks (ANN) and Greedy Algorithm (GA). In this method, GA is used to determine the correction order of SITS and calculate the error between the image to be corrected and normalized images, which avoids the selection of a single reference image. ANN is used to obtain the optimal solution of error function, which minimizes the radiometric distortion between different images in SITS. The SITS composed of 21 Landsat-8 images in Tianjin, China, from October 2017 to January 2019 was selected to test the method. We compared NMAG with other two contrast methods (Contrast Method 1 (CM1) and Contrast Method 2 (CM2)), and found that the average root mean square error (μRMSE) of NMAG (497.22) is significantly smaller than those of CM1 (641.39) and CM2 (543.47), and the accuracy of normalized SITS obtained using NMAG increases by 22.4% and 8.5% compared with CM1 and CM2, respectively. These experimental results confirm the effectiveness of NMAG in reducing radiometric distortion caused by continuously changing noise between images in SITS.

A Hybrid Deep Learning-Based Spatiotemporal Fusion Method for Combining Satellite Images with Different Resolutions

Spatiotemporal fusion (STF) is considered a feasible and cost-effective way to deal with the trade-off between the spatial and temporal resolution of satellite sensors, and to generate satellite images with high spatial and high temporal resolutions. This is achieved by fusing two types of satellite images, i.e., images with fine temporal but rough spatial resolution, and images with fine spatial but rough temporal resolution. Numerous STF methods have been proposed, however, it is still a challenge to predict both abrupt landcover change, and phenological change, accurately. Meanwhile, robustness to radiation differences between multi-source satellite images is crucial for the effective application of STF methods. Aiming to solve the abovementioned problems, in this paper we propose a hybrid deep learning-based STF method (HDLSFM). The method formulates a hybrid framework for robust fusion with phenological and landcover change information with minimal input requirements, and in which a nonlinear deep learning-based relative radiometric normalization, a deep learning-based superresolution, and a linear-based fusion are combined to address radiation differences between different types of satellite images, landcover, and phenological change prediction. Four comparative experiments using three popular STF methods, i.e., spatial and temporal adaptive reflectance fusion model (STARFM), flexible spatiotemporal data fusion (FSDAF), and Fit-FC, as benchmarks demonstrated the effectiveness of the HDLSFM in predicting phenological and landcover change. Meanwhile, HDLSFM is robust for radiation differences between different types of satellite images and the time interval between the prediction and base dates, which ensures its effectiveness in the generation of fused time-series data.

Random Sampling-Based Relative Radiometric Normalization

Relative radiometric normalization (RRN) is widely used for radiometric calibration of bitemporal multispectral images prior to any temporal analysis such as change detection. However, standard RRN methods are not robust against anomalous (or changed) pixels, which warp the calibration and decrease the spectral similarity of processed images. This letter proposes a novel random sample consensus-based RRN method, which only uses small pixel subsets to implement the linear mapping relationship for RRN, and does not require the calibration of its parameters. The experimental results show that the proposed method performs favorably against the widely used RRN methods in all metrics considered in this letter.

Comparison of Keypoint Detectors and Descriptors for Relative Radiometric Normalization of Bitemporal Remote Sensing Images

This article compares the performances of the most commonly used keypoint detectors and descriptors (SIFT, SURF, KAZE, AKAZE, ORB, and BRISK) in keypoint-based relative radiometric normalization (RRN) of unregistered bitemporal multispectral images. The keypoints matched between subject and reference images represent possible unchanged regions and form a radiometric control set (RCS). The initial RCS is further refined by removing the matched keypoints with a low cross-correlation. The final RCS is used to approximate a linear mapping between the corresponding bands of the subject and reference images. This procedure is validated on five datasets of unregistered multispectral image pairs acquired by inter/intra sensors in terms of RRN accuracy, visual quality, quality, and quantity of the samples in the RCS, and computational time. The experimental results show that keypoint-based RRN is robust against variations in spatial-resolution, illumination, and sensors. The blob detectors (SURF, SIFT, KAZE, and AKAZE) are more accurate on average than the corner detectors (ORB and BRISK) in RRN, with an expense of higher computational cost. The source code and samples of datasets used in this study are made available at https://github.com/ArminMoghimi/keypoint-based-RRN to support reproducible research in remote sensing.

Radiometric Normalization of Multitemporal Landsat and Sentinel-2 Images Using a Reference MODIS Product Through Spatiotemporal Filtering

Radiometric normalization is an essential preprocessing step for almost all remote sensing applications such as change detection, image mosaic, and 3-D reconstruction. This article proposes a novel radiometric normalizing method based on spatiotemporal filtering using a reference moderate resolution imaging spectroradiometer (MODIS) product. This differs from traditional relative radiometric normalization (RRN) methods in two folds: first, the number of reference images is more than one, which introduces more complexities than RRN with a single reference image; second, the resolution of MODIS product is significantly lower, thus requiring the algorithms to accommodate scale differences. To address, our approach extends the traditional spatiotemporal filtering method with per image bias that represents both internal (e.g., sensor characteristics) and external (e.g., atmosphere and topography) against the reference data. In addition, we use the Kullback-Leibler divergence metric to statistically determine the resemblance degree between the temporal images for weighting. We applied our proposed method to normalize Landsat Operational Land Imager, Enhanced Thematic Mapper Plus +, and Sentinel MSI using MODIS Nadir BRDF-adjusted reflectance product, covering two study areas of 30 × 15 km2 and 32 × 52 km2, respectively, and we show a notable radiometric consistency over both temporal and spatial dimension after the processing through three comparative experiments with state-of-the-art methods. 1) 3–7% improvement in the contexts of transfer learning, which favors only images with consistent radiometric properties and 2) Mosaic results using our processed images show no apparent seamlines as compared with images processed by other methods.

Automatically Monitoring Impervious Surfaces Using Spectral Generalization and Time Series Landsat Imagery from 1985 to 2020 in the Yangtze River Delta

Accurately monitoring the spatiotemporal dynamics of impervious surfaces is very important for understanding the process of urbanization. However, the complicated makeup and spectral heterogeneity of impervious surfaces create difficulties for impervious surface monitoring. In this study, we propose an automatic method to capture the spatiotemporal expansion of impervious surfaces using spectral generalization and time series Landsat imagery. First, the multitemporal compositing and relative radiometric normalization methods were used to extract phenological information and ensure spectral consistency between reference imagery and monitored imagery. Second, we automatically derived training samples from the prior MSMT_IS30-2020 impervious surface products and migrated the surface reflectance of impervious surfaces in the reference period of 2020 to other periods (1985–2015). Third, the random forest classification method, trained using the migrated surface reflectance of impervious surfaces and pervious surface training samples at each period, was employed to extract temporally independent impervious surfaces. Further, a temporal consistency check method was applied to ensure the consistency and reliability of the monitoring results. According to qualitative and quantitative validation results, the method achieved an overall accuracy of 90.9% and kappa coefficient of 0.859 in identifying the spatiotemporal expansion of impervious surfaces and performed better in capturing the impervious surface dynamics when compared with other impervious surface datasets. Lastly, our results indicate that a rapid increase of impervious surfaces was observed in the Yangtze River Delta, and the area of impervious surfaces in 2000 and 2020 was 1.86 times and 4.76 times that of 1985, respectively. Therefore, it could be concluded that the proposed method offered a novel perspective for providing timely and accurate impervious surface dynamics.

Radiometric normalization of Landsat thermal imagery for detection of tundra land cover changes: experience from West Siberia

ABSTRACT The potentiality and specificity of using multitemporal stacks of thermal infrared (TIR) satellite scenes for modelling changes in tundra land cover and predicting geohazards in areas under development are studied using the example of two areas in northern West Siberia. The study aims at developing and testing a method for relative radiometric normalization of Land Remote-Sensing Satellite (Landsat) TIR images applied to estimate changes in land surface temperature (LST) in the tundra with regard to the high-latitude climate conditions. The suggested normalization method (RNUC) differs from other similar approaches in the use of unsupervised classification to determine pseudo-invariant features where master and slave images are spectrally identical. The error in LST contrasts is minimized iteratively by choosing the optimal number of land-cover classes and linear regression coefficients by cross-validation. The normalization procedure ensures error reduction by a factor of two to four, which is evidence of its high performance. The results of processing and analysis of 19 pairs of Landsat thermal images at the sites of the Marre-Sale and Urengoy oil and gas-condensate fields using the RNUC method showed that the relative error of LST contrast normalization could be in the range of approximately 3% to 10%. Based on the time series of normalized values of LST and the normalized difference vegetation index (NDVI) of sites with pyrogenic damage in the area of the Urengoy oil and gas-condensate field, it is shown that the process of stabilization of the tundra cover is characterized by a negative correlation of these parameters and lasts 25 to 28 years. A positive correlation between LST and NDVI characterizes a different type of landscape transformation; in particular, it may be an increase in the thickness of secondary vegetation cover or a change in the hydrological conditions.