Jinhu Bian

An Enhanced Spatial and Temporal Data Fusion Model for Fusing Landsat and MODIS Surface Reflectance to Generate High Temporal Landsat-Like Data

Remotely sensed data, with high spatial and temporal resolutions, can hardly be provided by only one sensor due to the tradeoff in sensor designs that balance spatial resolutions and temporal coverage. However, they are urgently needed for improving the ability of monitoring rapid landscape changes at fine scales (e.g., 30 m). One approach to acquire them is by fusing observations from sensors with different characteristics (e.g., Enhanced Thematic Mapper Plus (ETM+) and Moderate Resolution Imaging Spectroradiometer (MODIS)). The existing data fusion algorithms, such as the Spatial and Temporal Data Fusion Model (STDFM), have achieved some significant progress in this field. This paper puts forward an Enhanced Spatial and Temporal Data Fusion Model (ESTDFM) based on the STDFM algorithm, by introducing a patch-based ISODATA classification method, the sliding window technology, and the temporal-weight concept. Time-series ETM+ and MODIS surface reflectance are used as test data for comparing the two algorithms. Results show that the prediction ability of the ESTDFM algorithm has been significantly improved, and is even more satisfactory in the near-infrared band (the contrasting average absolute difference [AAD]: 0.0167 vs. 0.0265). The enhanced algorithm will support subsequent research on monitoring land surface dynamic changes at finer scales.

Estimating the Maximal Light Use Efficiency for Different Vegetation through the CASA Model Combined with Time-Series Remote Sensing Data and Ground Measurements

Maximal light use efficiency (LUE) is an important ecological index of a vegetation essential attribute, and a key parameter of the LUE-based model for estimating large-scale vegetation productivity by remote sensing technology. However, although currently used in different models there still exists extensive controversy. This paper takes the Zoige Plateau in China as a case area to develop a new approach for estimating the maximal LUEs for different vegetation. Based on an existing land cover map and MODIS NDVI product, the linear unmixing method with a moving window was adopted to estimate the time-series NDVI for different end members in a MODIS NDVI pixel; then Particle Swarm Optimizer (PSO) was applied to search for the optimization of LUE retrievals through the CASA (Carnegie-Ames-Stanford Approach) model combined with time-series NDVI and ground measurements. The derived maximal LUEs present significant differences among various vegetation types. These are 0.669 gC.MJ(-1), 0.450 gC.MJ(-1) and 0.126 gC.MJ(-1) for the xerophilous grasslands with high, moderate and low vegetation fraction respectively, 0.192 gC.MJ(-1) for the hygrophilous grasslands, and 0.125 gC.MJ(-1) for the helobious grasslands. The field validation shows that the estimated net primary productivity (NPP) by the derived maximal LUE is closely related to the ground references, with R-2 of 0.8698 and root-mean-square error (RMSE) of 59.37 gC.m(-2).a(-1). This indicates that the default set in the CASA model is not suitable for NPP estimation for the regional mountain area. The derived maximal LUEs can significantly improve the capability of NPP mapping, and open up the perspective for long-term monitoring of vegetation ecological health and ecosystem productivity by combining the LUE-based model with remote sensing observations.

Auto-registration and orthorecification algorithm for the time series HJ-1A/B CCD images

How to deal with geometric distortion is an open problem when using the massive amount of satellite images at a national or global scale, especially for multi-temporal image analysis. In this paper, an algorithm is proposed to automatically rectify the geometric distortion of time-series CCD multispectral data of small constellation for environmental and disaster mitigation (HJ-1A/B) which was launched by China in 2008. In this algorithm, the area-based matching method was used to automatically search tie points firstly, and then the polynomial function was introduced to correct the systematic errors caused by the satellite motion along the roll, pitch and yaw direction. The improved orthorectification method was finally used to correct pixel displacement caused by off-nadir viewing of topography, which are random errors in the images and cannot be corrected by the polynomial equation. Nine scenes of level 2 HJ CCD images from one path/row were taken as the warp images to test the algorithm. The test result showed that the overall accuracy of the proposed algorithm was within 2 pixels (the average residuals were 37.8 m, and standard deviations were 19.8 m). The accuracies of 45.96% validation points (VPs) were within 1 pixel and 90.33% VPs were within 2 pixels. The discussion showed that three main factors including the distortion patterns of HJ CCD images, percent of cloud cover and the varying altitude of the satellite orbit may affect the search of tie points and the accuracy of results. Although the influence of varying altitude of the satellite orbits is less than the other factors, it is noted that detailed satellite altitude information should be given in the future to get a more precise result. The proposed algorithm should be an efficient tool for the geo-correction of HJ CCD multi-spectral images.

Development of Dense Time Series 30-m Image Products from the Chinese HJ-1A/B Constellation: A Case Study in Zoige Plateau, China

Time series remote sensing products with both fine spatial and dense temporal resolutions are urgently needed for many earth system studies. The development of small satellite constellations with identical sensors affords novel opportunities to provide such kind of earth observations. In this paper, a new dense time series 30-m image product was proposed respectively based on an 8-day, 16-day and monthly composition. The products were composited by the Charge Coupled Device (CCD) images from the 2-day revisit small satellite constellation for environmental monitoring and disaster mitigation of China (HJ-1A/B). Taking the Zoige plateau in China as a case area where it is covered by highly heterogeneous vegetation landscapes, a detailed methodology was introduced on how to use 183 scenes of CCD images in 2010 to create composite products. The quality of the HJ CCD composites was evaluated by inter-comparison with the monthly 30-m global Web-Enabled Landsat Data (WELD), 16-day 500-m MODIS NDVI, and 8-day 500-m MODIS surface reflectance products. Results showed that the radiometric consistency between HJ and WELD composited Top Of Atmosphere (TOA) reflectance was in good agreement except for May, June, July and August when more clouds and invalid data gaps appeared in WELD. Visual assessment and temporal profile analysis also revealed that HJ possessed better visual effects and temporal coherence than that of WELD. The comparison between HJ and MODIS products indicated that HJ composites were radiometrically consistent with MODIS products over areas consisting of large patches of homogeneous surface types, but can better reflect the detailed spatial differences in regions with heterogeneous landscapes. This paper highlights the potential of compositing HJ-1A/B CCD images, allowing for providing a cloud free, time-space consistent, 30-m spatial resolution, and dense in time series image product. Meanwhile, the proposed products could also be treated as a reference to generate regional or even global composited products for the on-going satellite constellations and even for the forthcoming satellite missions such as Sentinel-2A/B.

An Improved Physics-Based Model for Topographic Correction of Landsat TM Images

Optical remotely sensed images in mountainous areas are subject to radiometric distortions induced by topographic effects, which need to be corrected before quantitative applications. Based on Li model and Sandmeier model, this paper proposed an improved physics-based model for the topographic correction of Landsat Thematic Mapper (TM) images. The model employed Normalized Difference Vegetation Index (NDVI) thresholds to approximately divide land targets into eleven groups, due to NDVIs lower sensitivity to topography and its significant role in indicating land cover type. Within each group of terrestrial targets, corresponding MODIS BRDF (Bidirectional Reflectance Distribution Function) products were used to account for land surfaces BRDF effect, and topographic effects are corrected without Lambertian assumption. The methodology was tested with two TM scenes of severely rugged mountain areas acquired under different sun elevation angles. Results demonstrated that reflectance of sun-averted slopes was evidently enhanced, and the overall quality of images was improved with topographic effect being effectively suppressed. Correlation coefficients between Near Infra-Red band reflectance and illumination condition reduced almost to zero, and coefficients of variance also showed some reduction. By comparison with the other two physics-based models (Sandmeier model and Li model), the proposed model showed favorable results on two tested Landsat scenes. With the almost half-century accumulation of Landsat data and the successive launch and operation of Landsat 8, the improved model in this paper can be potentially helpful for the topographic correction of Landsat and Landsat-like data.

Land Cover Mapping in Southwestern China Using the HC-MMK Approach

Land cover mapping in mountainous areas is a notoriously challenging task due to the rugged terrain and high spatial heterogeneity of land surfaces as well as the frequent cloud contamination of satellite imagery. Taking Southwestern China (a typical mountainous region) as an example, this paper established a new HC-MMK approach (Hierarchical Classification based on Multi-source and Multi-temporal data and geo-Knowledge), which was especially designed for land cover mapping in mountainous areas. This approach was taken in order to generate a 30 m-resolution land cover product in Southwestern China in 2010 (hereinafter referred to as CLC-SW2010). The multi-temporal native HJ (HuanJing, small satellite constellation for disaster and environmental monitoring) CCD (Charge-Coupled Device) images, Landsat TM (Thematic Mapper) images and topographical data (including elevation, aspect, slope, etc.) were taken as the main input data sources. Hierarchical classification tree construction and a five-step knowledge-based interactive quality control were the major components of this proposed approach. The CLC-SW2010 product contained six primary categories and 38 secondary categories, which covered about 2.33 million km(2) (accounting for about a quarter of the land area of China). The accuracies of primary and secondary categories for CLC-SW2010 reached 95.09% and 87.14%, respectively, which were assessed independently by a third-party group. This product has so far been used to estimate the terrestrial carbon stocks and assess the quality of the ecological environments. The proposed HC-MMK approach could be used not only in mountainous areas, but also for plains, hills and other regions. Meanwhile, this study could also be used as a reference for other land cover mapping projects over large areas or even the entire globe.

Simulation of the Grazing Effects on Grassland Aboveground Net Primary Production Using DNDC Model Combined with Time-Series Remote Sensing Data—A Case Study in Zoige Plateau, China

Measuring the impact of livestock grazing on grassland above-ground net primary production (ANPP) is essential for grass yield estimation and pasture management. However, since there is a lack of accurate and repeatable techniques to obtain the details of grazing locations and stocking rates at the regional scale, it is an extremely challenging task to study the influence of regional grazing on the grassland ANPP. Taking Zoige County as a case, this paper proposes an approach to quantify the spatial and temporal variation of grazing intensity and grazing period through time-series remote sensing data, simulated grassland ANPP through the denitrification and decomposition (DNDC) model, and then explores the impact of grazing on grassland ANPP. The result showed that the model-estimated ANPP while considering grazing had a significant relationship with the field-observed ANPP, with the coefficient of determination (R-2) of 0.75, root mean square error (RMSE) of 122.86 kgC/ha, and average relative error (RE) of 8.77%. On the contrary, if grazing activity was not considered in simulation, a large uncertainty was found when the model-estimated ANPP was compared with the field observation, showing R-2 of 0.4, RMSE of 211.51 kgC/ha, and average RE of 32.5%. For the whole area of Zoige County in 2012, the statistics of the estimation showed that the total regional ANPP was up to 3.815 x 10(5) tC, while the total regional ANPP, without considering grazing, would be overestimated by 44.4%, up to 5.51 x 10(5) tC. This indicates that the grazing parameters derived in this study could effectively improve the accuracy of ANPP simulation results. Therefore, it is feasible to combine time-series remote sensing data with the process model to simulate the grazing effects on grassland ANPP. However, some issues, such as selecting proper remote sensing data, improving the quality of model input parameters, collecting more field data, and exploring the data assimilation approaches, still should be considered in the future work.

Cloud and Snow Discrimination for CCD Images of HJ-1A/B Constellation Based on Spectral Signature and Spatio-Temporal Context

It is highly desirable to accurately detect the clouds in satellite images before any kind of applications. However, clouds and snow discrimination in remote sensing images is a challenging task because of their similar spectral signature. The shortwave infrared (SWIR, e.g., Landsat TM 1.55-1.75 mu m band) band is widely used for the separation of cloud and snow. However, for some sensors such as the CBERS-2 (China-Brazil Earth Resources Satellite), CBERS-4 and HJ-1A/B (HuanJing (HJ), which means environment in Chinese) that are designed without SWIR band, such methods are no longer practical. In this paper, a new practical method was proposed to discriminate clouds from snow through combining the spectral reflectance with the spatio-temporal contextual information. Taking the Mt. Gongga region, where there is frequent clouds and snow cover, in China as a case area, the detailed methodology was introduced on how to use the 181 scenes of HJ-1A/B CCD images in the year 2011 to discriminate clouds and snow in these images. Visual inspection revealed that clouds and snow pixels can be accurately separated by the proposed method. The pixel-level quantitative accuracy validation was conducted by comparing the detection results with the reference cloud masks generated by a random-tile validation scheme. The pixel-level validation results showed that the coefficient of determination (R-2) between the reference cloud masks and the detection results was 0.95, and the average overall accuracy, precision and recall for clouds were 91.32%, 85.33% and 81.82%, respectively. The experimental results confirmed that the proposed method was effective at providing reasonable cloud mask for the SWIR-lacking HJ-1A/B CCD images. Since HJ-1A/B have been in orbit for over seven years and these satellites still run well, the proposed method is helpful for the cloud mask generation of the historical archive HJ-1A/B images and even similar sensors.

Performance Evaluation of the Triangle-Based Empirical Soil Moisture Relationship Models Based on Landsat-5 TM Data and In Situ Measurements

Surface soil moisture (SSM) is an important parameter at the land–atmosphere interface. In past decades, passive microwave remote sensing offers a good opportunity for obtaining SSM on a global scale, and many downscaling methods have been proposed using the triangle-based empirical soil moisture relationship models to overcome the limitation of coarse spatial resolution of its SSM products for regional applications. This paper aimed to examine and compare the effectiveness of five typical triangle-based empirical soil moisture relationship models for estimating SSM with Landsat-5 data and in situ measurements from the Maqu network on the northeastern part of the Tibetan Plateau for nine cloud-free days. The results showed that the model that treats the SSM as a second-order polynomial with land surface temperature, vegetation indices (VIs), and surface albedo as inputs exhibited the best performance compared with the results of other models. The VI comparison indicated that the use of the normalized difference VI or the fractional vegetation cover in this model outperformed other VIs, with the root-mean-square deviation of approximately 0.055 m3/m3 and the coefficient of determination (

China land cover monitoring in mountainous regions by remote sensing technology — Taking the Southwestern China as a case

\text{R}^{2}