Open access peer-reviewed chapter

Establishing the Downscaling and Spatiotemporal Scale Conversion Models of NDVI Based on Fractal Methodology

Written By

Haijun Luan

Submitted: May 30th, 2019 Reviewed: January 26th, 2020 Published: March 25th, 2020

DOI: 10.5772/intechopen.91359

Chapter metrics overview

625 Chapter Downloads

View Full Metrics


Scale effect is a crucial scientific problem in quantitative remote sensing (RS), and scholars attempt to solve it with scale conversion models, which can characterize the numerical relationship of RS land surface parameters at different resolutions (scales). As a significant land surface parameter, scale conversion of normalized difference vegetation index (NDVI) has been studied for a long time. Therefore, taking NDVI as an example, the development of scaling research is described and analyzed in the paper, and based on fractal theory, the development trends are discussed for land surface parameters in quantitative remote sensing. These are our conclusions: it will be the new trend to establish downscaling models based on fractal theory for land surface parameters in quantitative remote sensing; additionally, it still is the hotspot to establish spatiotemporal scale conversion models for land surface parameters in quantitative remote sensing in the future, and addressed on that, the multi-fractal scaling methodology is proposed, and its availability is analyzed in the paper, which presents significant potential.


  • remote sensing (RS)
  • normalized difference vegetation index (NDVI)
  • scaling
  • fractal
  • iterated function system (IFS)
  • multi-fractal

1. Introduction

The scale problem is one of the important and fundamental problems of quantitative remote sensing [1, 2, 3]. Scholars have studied the scale effects of different remote sensing (RS) land surface parameters. The study of scale effect is conducive to the synergistic use of RS data of different spatial and temporal resolutions (scales) to solve the problem that “massive” RS images cannot be fully utilized and has important application potential and scientific research value [3]. In view of the spatiotemporal characteristics of the ground objects, the RS land surface parameters not only have spatial scale effects but also have temporal scale effects. Scholars have conducted extensive and in-depth research on the scale effect of land surface parameters, which includes the mechanism, manifestation, effect analysis, and solution of scale effects. The author has previously discussed it in detail [4]. Based on the above research aspects, scale conversion as a solution to scale effects has received attention. The scale conversion model can characterize the numerical or physical relationship of RS land surface parameter images at different resolutions (scales) and can quantitatively describe scale effects. This paper will also focus on the research progress of the spatial down-scaling and the spatiotemporal scaling.


2. Downscaling of NDVI based on fractal IFS

2.1 Review of downscaling of RS land surface parameters

Liang [1] has reviewed several current downscaling methods, including linear decomposition methods and nonlinear statistical decomposition methods, methods for generating continuous regions, normalized difference vegetation index (NDVI) time series decomposition, multi-resolution data fusion, the statistical downscaling method of global climate model products (GCM), etc. Further, Gao et al. [5], Zhu et al. [6], and Huang et al. [7, 8] have done systematic and effective work in the spatiotemporal fusion downscaling of land surface reflectance, which has become a research hot topic. The spectral-spatial feature fusion by Wang et al. [9, 10, 11, 12] and Shi and Wang [13] also achieved good results for subpixel mapping. These studies, however, scarcely considered the scale conversion process from the perspective of dynamics, which studies of surface parameter downscaling based on the fractal iterated function system (IFS) have paid attention to.

As a fractal branch of mathematics, because of its complete and rigorous theoretical system, it can systematically study the performance, nature, and causes of multi-scale characteristics of natural phenomena. In the fractal geometry theory system, in addition to the familiar fractal phenomenon description and fractal measurement, the internal causes or dynamic processes of mathematical fractals (interaction, feedback, and iteration, represented by IFS-iteration function system) and the physical causes of statistical fractals (such as critical or abrupt changes) are also important research contents of fractal geometry, and fractal geometry has become a part of nonlinear dynamics research [14]. Although the current research on fractal dynamics has just started, there are still many problems waiting to be solved, but its potential value and significance in dynamics research cannot be denied.

In quantitative remote sensing research, fractal methods are mostly used in the mapping of surface morphology (spatial structure) such as active radar imagery and snow and ocean imagery [15], but it also has important applications in scale conversion research and is further deepened and expanded. The use of fractals for surface parameter scale conversion modeling usually contains two important research components:

  1. The performance of fractal features, that is, fractal metrics, and also the fractal dimension of the research object. For example, Zhang et al. [16, 17] used the information dimension method to describe the fractal dimension of leaf area index (LAI) scale conversion. Luan et al. [18, 19] and Wu et al. [20] used the similar dimension method to measure the fractal dimension of NDVI and LAI scale-up conversion, respectively.

  2. The intrinsic nature of the fractal phenomenon, that is, the dynamics produced, which is the combined effect of multifactor surface effects.

The mathematical basis of fractal generation is IFS. Kim and Barros [21] first constructed the r function from the dynamic factors (soil sediment content, vegetation water content) of soil moisture scale conversion and then established the IFS to describe the soil moisture downscaling, and the conversion effect was good. The model can describe the dynamic process of soil moisture scale conversion, which has physical significance and demonstrates the advantages of downscaling surface parameters based on fractal IFS. In general, there is currently little research into the causes of fractal dynamics. In mathematics, the fractal IFS is a continuously iterative calculation based on the whole research object [14], and the RS land surface parameter image is created in units of local pixels. This ensures that the mathematical IFS vertical conversion factor (r function) is usually constant [21], while the vertical conversion factor of RS land surface parameters (such as soil moisture) is based on the physical elements of each pixel (such as sandy soil). The amount of space and the vegetation water content varies dynamically and temporally [21]. This is why the IFS function can describe the scale switching dynamics of surface parameters and why the model has certain physical meanings. The vertical conversion factor is used to describe the interscale conversion of surface parameter values and is the key to determining the IFS function. Different surface parameters have different values due to the spatial distribution and scale conversion factors (or dynamic factors), and the vertical conversion factor (r function) contains different types of variables and function forms. How to determine the r function is the difficulty in determining the IFS function, which is also an important reason why the latter is less frequently applied in descriptions of quantitative RS land surface parameter scale conversion. Therefore, the NDVI downscaling model based on the fractal IFS function can be considered to describe the dynamic process of scale conversion. This research covers a wide area and is of great significance. The following is a description of a preliminary implementation [22].

2.2 Methodology

How does one build an NDVI downscaling model based on the fractal IFS function? The following points need to be considered: first, how to identify the sensitive factors affecting the spatial distribution and scale effect of NDVI for NDVI; second, how to use this sensitive factor to establish the vertical scale conversion factor r function in the IFS and then determine the IFS function to achieve NDVI downscaling; and finally, how to evaluate the downscaled conversion results. The solution incorporating these considerations is described below.

2.2.1 Identify sensitive factors

According to the above description, water body is an important parameter affecting the spatial distribution and scale effect of NDVI; thus it can be determined that the pixel water parameter is one of the important dynamic factors of NDVI scale conversion. In addition, Wen et al. [23] gave a method for albedo conversion from small-scale to large-scale images and used the pixel topographical influencing factors to correct the converted results, which demonstrated that the method was effective for albedo scale conversion of rugged terrain. Considering the close relationship between the surface reflectivity and the surface albedo, and that the surface reflectance is the basic parameter for calculating NDVI, the topographic factor parameter can be determined as one of the important kinetic factors for NDVI scale conversion. Therefore, the important dynamic factors in NDVI spatial distribution and scale conversion are determined to be the pixel water parameters and topographic factors.

2.2.2 Determine the vertical conversion factor r function and establish the IFS function

Referring to Kim [21], IFS formula (1), horizontal transformation formula (2), and vertical transformation formula (3) for large-scale surface parameter pixel downscaling are obtained as follows. The IFS formula is calculated by pixel-by-pixel sliding. Get the full image downscaling results:


where IFSi,jn,mxiyjsij represents the surface parameter of the pixel at the ij location when the large-scale pixel of the surface parameter is downscaled to the small-scale image of the n×m dimension; xi, yj, and sij correspond, respectively, to the x-direction coordinate pnxi, the y-direction coordinate qmyj, and the surface parameter values In,mxiyjsij of the three-dimensional data of the pixel; xn1i and x0i represent, respectively, the x-direction starting coordinate of the ij pixel in the n×m dimensional small-scale image and the x-direction starting coordinate of the large-scale pixel; α represents the downscaling ratio (small-scale/large-scale, which is less than or equal to 1); en,m, fn,m, gn,m, and kn,m are, respectively, functions of the x and y coordinates of the lower left corner and upper right corner of the large-scale pixel, the downscaled surface parameter data, and the vertical scale conversion surface function; and r1xiyj and r2xiyj represent, respectively, the two different vertical conversion factors in the vertical scale conversion surface function. Reference [21] should be consulted for the parameters or factors not represented in the formula, which will not be explained here. Generally, the pnxi and qmyj coordinates of the ij pixel are obtained by dividing the large-scale pixel equally into 1/α parts, and the In,mxiyjsij calculation is the key. In formula (3), r2xiyj is the same as the r1xiyj function, but their argument coefficients are different.

For NDVI, gn,m, en,m, fn,m, and kn,m represent the functions of thenm pixel of the downscaled NDVI image. Based on the special downscaled NDVI 3D values of the four corner pixels, gn,m, en,m, fn,m, and kn,m can be calculated as formulas (4)(11):




Therefore, the calculation of the r1(xi,yj) function is significant, and r1nm (0r11) is used to adjust the NDVI surface roughness. The following treatment focuses on establishing the vertical transformation formula for NDVI, that is, the determination of the r function (containing r1xiyj and r2xiyj).

Based on the above sensitivity factors, a vertical conversion factor r function can be constructed:


where Swater represents the pixel water parameter; s represents the topographic information, taking into account the magnitude of the r function; the normalized difference water index (NDWI) and slope (calculated from the digital elevation model (DEM) image) represent, respectively, the water body effect and the topographic influence in the pixel; γ and β are the coefficients of the two parameters, respectively; and δ represents the adjustment constant. Two different orders of magnitude of r are calculated as follows:


For NDVI, the γ, β, and δ coefficients can be calculated by linear regression between the high-resolution NDVI image and its NDWI/slope images.

Following construction of the r function, formulas (1)(3) can be solved in combination with other known conditions, and NDVI downscaling can be achieved.

2.2.3 Evaluation of downscaling results

In order to obtain more accurate downscaling results, if the resolution of the low resolution image is too different from the resolution of the target resolution image (such as downscaling from 250 m MODIS NDVI to 30 m NDVI), a hierarchical downscaling method will be adopted. First, the low-resolution surface parameter image is downscaled to an intermediate resolution image, and then the intermediate resolution image is further downscaled to the target resolution image, which can largely guarantee the accuracy of the result.

Referring to the study by Kim and Barros [21], the accuracy of the downscaled results can be evaluated using statistical indicators such as the maximum, minimum, variance, and standard deviation (compared to high-resolution NDVI images). Moreover, the histograms of the downscaled NDVI and true NDVI images were drawn and compared, and their correlation coefficient was calculated. With those indexes, the accuracy of the downscaled images and methodology could be validated.

2.3 Experiment and result analysis

2.3.1 Experiment

As the best indicator of the status of vegetation growth and vegetation coverage, the normalized difference vegetation index is widely used in the study of environmental (climate) changes, crop yield estimation, and other fields. Among existing vegetation index products, the moderate resolution imaging spectroradiometer (MODIS) vegetation index products are highly valued for their ease of use, ready availability, global coverage, and continuous phase. They have been widely used in studies of forest fires [24, 25], grassland vegetation growth [26, 27], drought [28, 29], land desertification [30], and other studies involving ecological environment monitoring. The maximum spatial resolution of MODIS vegetation index products, however, is only 250 m. The validation of this remote sensing land surface parameter is an important issue that cannot be avoided [31, 32, 33] and needs to be carried out by means of scale conversion. The most representative MODIS NDVI product, namely, MOD13 Q1, will be studied in this paper, which will also focus on establishing a downscaled model of NDVI and validating the MOD13 Q1 product based on it.

This is the experiment. A Landsat8 OLI NDVI image (Figure 1) was utilized to validate a MODIS NDVI image (MOD13 Q1, Figure 2) with nearest imaging time in Xiamen, China. Based on the downscaling formulas in Section 2.2, the MOD13 Q1 image of Xiamen was directly downscaled by ⅛ multiples, and the 30 m downscaled NDVI was obtained as Figure 3. The histograms of the original and processed NDVI images are drawn as Figure 4, and the statistics and correlation coefficients of the NDVI images are presented in Table 1. Based on these data, the downscaled results were evaluated and the MOD13 Q1 image was validated.

Figure 1.

30 m OLI NDVI image.

Figure 2.

240 m MOD13 Q1 image.

Figure 3.

30 m downscaled image of MOD13 Q.

Figure 4.

Histograms of original and processed NDVI images. (a) MOD13 Q1, (b) 30 m downscaled MOD13 Q1, (c) OLI NDVI, (d) difference image between Figures 3 and 1 (Figure 3Figure 1).

NDVI imagesMaximumMinimumMeanVarianceCorrelation coefficient
MOD13 Q10.999100−0.5615000.3974000.001100
30 m downscaled MOD13 Q10.999950−0.9997930.3264120.2064950.937449
OLI NDVI0.984436−0.9612100.2994180.227313
Difference image (Figure 3Figure 1)1.672050−1.088930−0.0295940.027916

Table 1.

Statistics of original and processed NDVI images.

Note: The correlation coefficient is the one between the downscaled MOD13 Q1 and the OLI NDVI.

2.3.2 Result analysis

By analyzing Figures 14 and Table 1, it is found that:

  1. Compared with the real 30 m OLI NDVI image, the 30 m downscaled MOD13 Q1 image has smaller differences in maximum value, minimum value, mean value, and variance. The correlation coefficient between the two images is 0.93, which is highly correlated. The overall quality of the NDVI image obtained by downscaling the MOD13 Q1 image is considered to be good, indicating that the overall quality of MOD13 Q1 is good. Comparing Figure 4(a) and (b), there is a certain similarity between the distribution patterns of the two images, which indicates that the downscaled image retains the spatial distribution structure of the original image to a good degree, which proves to some extent that the original MOD13 Q1 image is of good quality. In addition, comparing Figure 4(b) and (c), it is found that the downscaled NDVI image has a higher proportion in the vicinity of the zero value (mainly artificial features) than the real image. In the range of 0.2–0.6, the difference is greater. The downscaled image generally has a higher proportion in this range of values, and the histogram is smoother, indicating that the image recognition of the NDVI difference is not high. Referring to the correlation between Figure 4(b) and (a), it is known that MOD13 Q1 also has these problems within the abovementioned range of values. Analysis of the original MOD13 Q1 image shows that it is a 16-day NDVI composite product, and each pixel takes the maximum value of NDVI within 16 days as the result of the product release. Therefore, the histogram is reasonable to a certain degree in the larger value area. At the same time, the histogram distribution of the difference image indicates that the pixel values are distributed in the range [−1, 1] and the distribution pattern is low on both sides and high in the middle (approximately a value of 0), which also indicates that the downscaled image and the real image are highly consistent.

  2. Further, the analysis of Table 1 shows that the maximum value of the difference image exceeds the range of [−1, 1]. This may be due to a certain error which is caused by the MOD13 Q1 and OLI image during preprocessing process (atmospheric correction, geometric correction, etc.), which causes a large difference in pixel values between the MOD13 Q1 downscale image and the OLI NDVI image. However, the analysis of Table 1 shows that the mean and variance of the difference image are small, so the above abnormal situation only occupies a small space and does not affect the overall evaluation conclusion.

According to the above analysis, the overall quality of the MOD13 Q1 downscaled image is good, indicating that the overall quality of MOD13 Q1 is good. In the NDVI range of values from 0.2 to 0.6, MOD13 Q1 is overestimated, and its discrimination ability of NDVI difference is low, which should be taken into account in practical applications.

2.4 Discussion

Based on the fractal iterated function system, downscaling models of remote sensing land surface parameters can be established. The models can then be merged with more ancillary data, which relate to the scale effects of land surface parameters. Therefore, the models are of benefit for obtaining accurate downscaled results.

In summary, although the breadth and depth of the fractal IFS application in establishing RS land surface parameters downscaling models is still insufficient, the inherent physical meaning and the advantages of the dynamic process expression of this method confer great potential on it, which needs further investigation. It is expected to become a new universal method for quantitative downscaling of RS land surface parameters and lead to the discovery of new research methods.


3. Establishing spatiotemporal scale conversion models of RS land surface parameters based on multi-fractal theory and method

3.1 Review of establishing spatiotemporal scale conversion models of RS land surface parameters

The phase is an important feature of RS images. When the phase changes, the spectrum of the objects in the image changes accordingly. Then, the parameters calculated based on the spectral information will also change, such as surface reflectivity, NDVI, and so on. The temporal response of RS land surface parameters will be further reflected in the variation of its spatial scale conversion model, i.e., the phase characteristics of spatial scale effects.

In order to quantitatively characterize the phase characteristics of spatial scale effects, that is, to establish a spatiotemporal scale conversion model (also called a spatiotemporal scaling fusion model), scholars combined the advantages of higher temporal-resolution feature of low spatial-resolution images and higher spatial-resolution feature of medium spatial-resolution images, and carried out a series of studies on spatiotemporal fusion of remotely sensed surface parameters such as surface reflectance [34], land surface temperature [35, 36], vegetation indexes [37], leaf area index [38], and so on. And then Huang et al. [39] reviewed this and presented the systematic achievements in theory and application. From the theoretical basis of spatiotemporal scale conversion fusion (the spatial scale consistency of the time-phase variation model and the time consistency of the spatial downscaling model) to the type division of the spatiotemporal scale conversion fusion algorithms (the algorithms based on features’ components, the algorithms based on surface spatial information, the algorithms based on features’ temporal change, and the combination algorithms of the ones above), and then to the key problems and challenges encountered in existing research (the imaging geometry and radiation characteristics, differences between multi-source RS images, the complexity of subpixel unmixing models, the complexity of features’ temporal change models, etc.), and the possible development trend in future (improvement in the versatility and robustness of the algorithms), he made a detailed and in-depth explanation, so that we have a more comprehensive understanding of the development of spatiotemporal scale conversion fusion research. In fact, in addition to this method, the multi-fractal method has important potential to solve the above problems [21, 40]. The following is an example of NDVI analysis and how to establish a spatiotemporal scale conversion model (or spatiotemporal scaling fusion model) based on multi-fractal theory and method.

3.2 Spatiotemporal scale conversion models of NDVI based on multi-fractal theory and method

As the best indicator of vegetation growth status and vegetation coverage, NDVI has typical phenological characteristics. This means that in the same area where the surface cover type is unchanged, the physiological characteristics and external forms of the plant can change significantly in different growth stages, and this change will be directly reflected in the changes in image spectrum and NDVI of the surface. Furthermore, the NDVI spatial scale conversion model based on RS images of different growth periods (i.e., different phases) will also change. How to effectively reflect the influence of the phase characteristics of RS images on the construction of this model and then construct a more universal NDVI spatial scale conversion model that can be integrated with surface phenological features, namely, NDVI spatiotemporal scale conversion model? This issue has important research value. Kim and Barros [21] proposed the idea of multi-fractal method for multi-temporal remote sensing soil moisture spatial down-scaling model to describe the phase characteristics of soil moisture spatial down-scaling, but did not do specific research.

Referring to the existing knowledge, the specific method of establishing the NDVI spatiotemporal scale conversion model is given here: first, analyze the surface condition of the study area, determine the type of the main cover of the study area, and based on its phenological knowledge, select enough low and medium-high spatial resolution images finely corresponding to important “nodes” of vegetation throughout the growing season; secondly, the NDVI spatial downscaling models for different growing stages “nodes” are constructed based on the down-scaling methods such as fractal IFS; third, according to multi-fractal theory and method, using the time phase as a factor in the fractal dimension calculation method, the models corresponding to each growth stages are “fused” to obtain a unified and full growth period NDVI scale conversion model (i.e., NDVI spatiotemporal scale conversion model). At this time, the time phase (i.e., different growth stages) has been embodied as a parameter in the model. This model is more universal than the downscaling model based on the single phase image. To obtain a medium to high spatial resolution NDVI image of a certain phase during vegetation growth, the corresponding phase and the low spatial resolution NDVI image of the phase are brought into the model calculation. Of course, this method requires the research object to have a more significant phase or time periodicity, and the established spatiotemporal scale conversion model is more accurate.

Besides, there is another method of multi-fractal modeling of NDVI spatiotemporal scaling. The implementation idea is similar with Section 2.2, while the r function changes. The r functional parameters may need to be recalibrated when the spatial distribution of vegetation cover changes obviously with time (e.g., sowing stage, heading stage, maturity stage, etc.). Therefore, r function will be merged with temporal parameters of NDVI distribution, such as LAI. And the multi-fractal model of NDVI spatiotemporal scaling should be a function of NDVI to capture temporal changes in relation to ancillary data such as LAI.

Although the multi-fractal theory and method has advantages in constructing a spatiotemporal scale conversion model of RS land surface parameters, the theory and implementation of this method are more complicated, and few research cases are currently seen. However, this method is expected to become a new method for the construction of spatiotemporal scale conversion model of RS land surface parameters, which is worthy of further study.


4. Conclusions

Taking normalized difference vegetation index (NDVI) as an example, the establishment of scaling models based on fractal theory was described and analyzed in the paper. It was concluded that fractal iterated function system was an effective methodology to establish downscaling models for remote sensing land surface parameters such as NDVI and multi-fractal modeling may be a novel methodology to establish spatiotemporal scale conversion models for land surface parameters such as NDVI in the future.



This work was supported by the National Natural Science Foundation of China “Coupling of NDVI’s up-scaling and down-scaling fusing with ground objects classification” (No. 41601350), the Natural Science Foundation of Fujian Province, China “Research on NDVI’s scaling fusing with ground objects classification” (No. 2017J05069), and the “Scientific Research Climbing Plan” Project from Xiamen University of Technology “Spatial Distribution Estimation and Dynamic Monitoring of Soil Organic Matter Based on Multi-source and Heterogeneous Data” (No. XPDKT19010).


  1. 1. Liang S. Quantitative Remote Sensing. Translated by Fan WJ. Beijing: Science Press; 2009. pp. 180-183. (in Chinese)
  2. 2. Mohanty BP. Soil hydraulic property estimation using remote sensing: A review. Vadose Zone Journal. 2013;12(4). DOI: 10.2136/vzj2013.06.0100. Available from:
  3. 3. Li XW, Wang YT. Prospects on future developments of quantitative remote sensing. Acta Geographica Sinica. 2013;68(9):1163-1169. (in Chinese). DOI: 10.11821/dlxb201309001
  4. 4. Luan HJ, Tian QJ, Yu T, et al. Review of up-scaling of quantitative remote sensing. Advances in Earth Science. 2013;28(6):657-664. (in Chinese) 10.11867/j.issn.1001-8166.2013.06.0657
  5. 5. Gao F, Masek J, Schwaller M, et al. On the blending of the Landsat and MODIS surface reflectance: Predicting daily Landsat surface reflectance. IEEE Transactions on Geoscience and Remote Sensing. 2006;44(8):2207-2218. DOI: 10.1109/TGRS.2006.872081
  6. 6. Zhu XL, Chen JM, Gao F, et al. An enhanced spatial and temporal adaptive reflectance fusion model for complex heterogeneous regions. Remote Sensing of Environment. 2010;114(11):2610-2623. DOI: 10.1016/j.rse.2010.05.032
  7. 7. Huang B, Zhang HK, Song HH, et al. Unified fusion of remote sensing imagery: Generating simultaneously high-resolution synthetic spatial-temporal spectral earth observations. Remote Sensing Letters. 2013;4:561-569. DOI: 10.1080/2150704X.2013.769283
  8. 8. Huang B, Zhang HK. Spatio-temporal reflectance fusion via unmixing: Accounting for both phenological and land-cover changes. International Journal of Remote Sensing. 2014;35(16):6213-6233. DOI: 10.1080/01431161.2014.951097
  9. 9. Wang Q, Shi W, Wang L. Allocating classes for soft-then-hard subpixel mapping algorithms in units of class. IEEE Transactions on Geoscience and Remote Sensing. 2014;52(5):2940-2959. DOI: 10.1109/TGRS.2013.2267802
  10. 10. Wang Q, Shi W, Atkinson PM, et al. Downscaling MODIS images with area-to-point regression kriging. Remote Sensing of Environment. 2015;166:191-204. DOI: 10.1016/j.rse.2015.06.003
  11. 11. Wang Q, Atkinson PM, Shi W. Indicator cokriging-based subpixel mapping without prior spatial structure information. IEEE Transactions on Geoscience and Remote Sensing. 2015;53(1):309-323. DOI: 10.1109/TGRS.2014.2321834
  12. 12. Wang Q, Atkinson PM, Shi W. Fast sub-pixel mapping algorithms for sub-pixel resolution change detection. IEEE Transactions on Geoscience and Remote Sensing. 2015;53(4):1692-1706. DOI: 10.1109/TGRS.2014.2346535
  13. 13. Shi W, Wang Q. Soft-then-hard sub-pixel mapping with multiple shifted images. International Journal of Remote Sensing. 2015;36(5):1329-1348. DOI: 10.1080/01431161.2015.1009650
  14. 14. Chen Y, Chen L. Fractal Geometry. 2nd ed. Beijing: Earthquake Press; 2005. pp. 49-51, 95-98. (in Chinese)
  15. 15. Riccio D, Ruello G. Synthesis of fractal surfaces for remote-sensing applications. IEEE Transactions on Geoscience and Remote Sensing. 2015;53(7):3803-3814. DOI: 10.1109/TGRS.2014.2384595
  16. 16. Zhang RH, Tian J, Li ZL, et al. Spatial scaling and information fractal dimension of surface parameters used in quantitative remote sensing. International Journal of Remote Sensing. 2008;29:5145-5159. DOI: 10.1080/01431160802036581
  17. 17. Zhang RH, Tian J, Li ZL, et al. Principles and methods for the validation of quantitative remote sensing products. Science China Earth Sciences. 2010;53:741-751. DOI: 10.1007/s11430-010-0021-3
  18. 18. Luan HJ, Tian QJ, Gu XF, et al. Establishing continuous scaling of NDVI based on fractal theory and GEOEYE-1 image. Journal of Infrared and Millimeter Waves. 2013;32(6):538-544, 549. (in Chinese). DOI: 10.3724/SP.J.1010.2013.00538
  19. 19. Luan HJ, Tian QJ, Yu T, et al. Establishing continuous spatial scaling model of NDVI on fractal theory and five-index estimation system. Journal of Remote Sensing. 2015;19(1):116-125. (in Chinese). DOI: 1007-4619(2015)01-0116-10
  20. 20. Wu L, Qin Q, Liu X, et al. Spatial up-scaling correction for leaf area index based on the fractal theory. Remote Sensing. 2016;8(3):197. DOI: 10.3390/rs8030197
  21. 21. Kim G, Barros AP. Downscaling of remotely sensed soil moisture with a modified fractal interpolation method using contraction mapping and ancillary data. Remote Sensing of Environment. 2002;83:400-413. DOI: 10.1016/S0034-4257(02)00044-5
  22. 22. Luan HJ, Zhang M, Wan YY, et al. Establishing the downscaling model of NDVI based on the iterated function system. In: Lecture Notes in Electrical Engineering. In press
  23. 23. Wen JG, Liu Q, Liu QH, et al. Scale effect and scale correction of land-surface albedo in rugged terrain. International Journal of Remote Sensing. 2009;30(20):5397-5420. DOI: 10.1080/01431160903130903
  24. 24. Huesca M, Litago J, Palaciosorueta A, et al. Assessment of forest fire seasonality using MODIS fire potential: A time series approach. Agricultural and Forest Meteorology. 2009;149(11):1946-1955. DOI: 10.1016/j.agrformet.2009.06.022
  25. 25. Leon JR, Van Leeuwen WJ, Casady GM, et al. Using MODIS-NDVI for the modeling of post-wildfire vegetation response as a function of environmental conditions and pre-fire restoration treatments. Remote Sensing. 2012;4(3):598-621. DOI: 10.3390/rs4030598
  26. 26. Zhang B, Zhang L, Xie D, et al. Application of synthetic NDVI time series blended from Landsat and MODIS data for grassland biomass estimation. Remote Sensing. 2015;8(1):10. DOI: 10.3390/rs8010010
  27. 27. Nestola E, Calfapietra C, Emmerton CA, et al. Monitoring grassland seasonal carbon dynamics, by integrating MODIS NDVI, proximal optical sampling, and Eddy covariance measurements. Remote Sensing. 2016;8(3):260. DOI: 10.3390/rs8030260
  28. 28. Du L, Tian Q, Yu T, et al. A comprehensive drought monitoring method integrating MODIS and TRMM data. International Journal of Applied Earth Observation and Geoinformation. 2013;23(1):245-253. DOI: 10.1016/j.jag.2012.09.010
  29. 29. Kim Y. Drought and elevation effects on MODIS vegetation indices in northern Arizona ecosystems. International Journal of Remote Sensing. 2013;34(14):4889-4899. DOI: 10.1080/2150704X.2013.781700
  30. 30. Veron SR, Paruelo JM. Desertification alters the response of vegetation to changes in precipitation. Journal of Applied Ecology. 2010;47(6):1233-1241. DOI: 10.1111/j.1365-2664.2010.01883.x
  31. 31. Gao X, Huete AR, Didan K, et al. Multisensor comparisons and validation of MODIS vegetation indices at the semiarid Jornada experimental range. IEEE Transactions on Geoscience and Remote Sensing. 2003;41(10):2368-2381. DOI: 10.1109/TGRS.2003.813840
  32. 32. Fensholt R, Sandholt I, Stisen S, et al. Evaluating MODIS, MERIS, and VEGETATION vegetation indices using in situ measurements in a semiarid environment. IEEE Transactions on Geoscience and Remote Sensing. 2006;44(7):1774-1786. DOI: 10.1109/TGRS.2006.875940
  33. 33. Geng L, Ma M, Yu W, et al. Validation of the MODIS NDVI products in different land-use types using In situ measurements in the Heihe River basin. IEEE Geoscience and Remote Sensing Letters. 2014;11(9):1649-1653. DOI: 10.1109/LGRS.2014.2314134
  34. 34. Emelyanova IV, Mcvicar TR, Van Niel TG, et al. Assessing the accuracy of blending Landsat-MODIS surface reflectance in two landscapes with contrasting spatial and temporal dynamics: A framework for algorithm selection. Remote Sensing of Environment. 2013;133:193-209. DOI: 10.1016/j.rse.2013.02.007
  35. 35. Huang B, Wang J, Song H, et al. Generating high spatiotemporal resolution land surface temperature for urban Heat Island monitoring. IEEE Geoscience and Remote Sensing Letters. 2013;10(5):1011-1015. DOI: 10.1109/LGRS.2012.2227930
  36. 36. Kim J, Hogue TS. Evaluation and sensitivity testing of a coupled Landsat-MODIS downscaling method for land surface temperature and vegetation indices in semi-arid regions. Journal of Applied Remote Sensing. 2012;6(1):063569-1-063569-17. DOI: 10.1117/1.JRS.6.063569
  37. 37. Ouyang W, Hao F, Skidmore AK, et al. Integration of multi-sensor data to assess grassland dynamics in a Yellow River sub-watershed. Ecological Indicators. 2012;18(1):163-170. DOI: 10.1016/j.ecolind.2011.11.013
  38. 38. Zhang HK, Chen JM, Huang B, et al. Reconstructing seasonal variation of Landsat vegetation index related to leaf area index by fusing with MODIS data. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing. 2014;7(3):950-960. DOI: 10.1109/JSTARS.2013.2284528
  39. 39. Huang B, Zhao YQ. Research status and Prospect of spatiotemporal fusion of multi-source satellite remote sensign imagery. Acta Geodaetica et Cartographica Sinica. 2017;46(10):1492-1499. (in Chinese). DOI: 10.11947/j.AGCS.2017.20170376
  40. 40. Chen JM. Spatial scaling of a remotely sensed surface parameter by contexture. Remote Sensing of Environment. 1999;69:30-42. DOI: 10.1016/S0034-4257(99)00006-1

Written By

Haijun Luan

Submitted: May 30th, 2019 Reviewed: January 26th, 2020 Published: March 25th, 2020