Scaling Effect of Area-Averaged NDVI : Monotonicity along the Spatial Resolution

Changes in the spatial distributions of vegetation across the globe are routinely monitored by satellite remote sensing, in which the reflectance spectra over land surface areas are measured with spatial and temporal resolutions that depend on the satellite instrumentation. The use of multiple synchronized satellite sensors permits long-term monitoring with high spatial and temporal resolutions. However, differences in the spatial resolution of images collected by different sensors can introduce systematic biases, called scaling effects, into the biophysical retrievals. This study investigates the mechanism by which the scaling effects distort normalized difference vegetation index (NDVI). This study focused on the monotonicity of the area-averaged NDVI as a function of the spatial resolution. A monotonic relationship was proved analytically by using the resolution transform model proposed in this study in combination with a two-endmember linear mixture model. The monotonicity allowed the inherent uncertainties introduced by the scaling effects (error bounds) to be explicitly determined by averaging the retrievals at the extrema of the resolutions. Error bounds could not be estimated, on the other hand, for non-monotonic relationships. Numerical simulations were conducted to demonstrate the monotonicity of the averaged NDVI along spatial resolution. This study provides a theoretical basis for the scaling effects and develops techniques for rectifying the scaling effects in biophysical Remote Sens. 2012, 4 161 retrievals to facilitate cross-sensor calibration for the long-term monitoring of vegetation dynamics.


Introduction
Satellite remote sensing provides a historical record of the biophysical parameters that may be used to model the global vegetation dynamics [1,2], and the data also provides input variables for climate and surface process models [3][4][5].The cloud-free image series, which has been collected by multiple sensors, can improve the quality of the spatial and temporal data [6][7][8][9][10][11][12][13][14].For example, MODIS, which is on board the Terra satellite, covers the entire earth in about two days with a spatial resolution that is coarser than the resolution of images gathered by middle-resolution satellites (e.g., Terra-ASTER and SPOT-HRG).The low resolution of the MODIS data introduces uncertainties into the spatial map.The simultaneous use of satellite sensors (e.g., MODIS and ASTER) can compensate for the spatio-temporal gaps, producing a reliable parameterization of the terrestrial vegetation [15].
Differences in the sensor characteristics, including the spatial, radiometric, or spectral resolution, often introduce systematic biases in the arithmetically averaged (area-averaged) biophysical parameters [16][17][18][19], even though the data are gathered over a single field.The biases impose limitations on the product continuity [20], and they degrade the accuracy of the retrievals and the climate change predictions [21].In view of the vast archives containing past and future satellite data, cross-sensor calibration is required for the accurate implementation of long-term coherent monitoring using multiple sensors [22,23].
This study focuses on the systematic bias present in retrievals (parameters) due to differences in the spatial resolution of the data [24][25][26], called 'scaling effects'.The term 'scaling effects' has several meanings in remote sensing, depending on the issues to be addressed.Chen [27] categorized the scaling effects into three types: (1) the effects of the area-averaging operation and the point spread function [28,29]; (2) the effects of surface heterogeneity and the derivation algorithms on the area-averaged retrievals [24,25]; and (3) the effects of surface heterogeneity and the correlations between the surface and atmospheric variables involved in the estimation processes [30,31].In this study, we focused on the second type of scaling effects, which have been widely discussed in remote sensing studies.
In practice, biophysical parameter such as leaf area index, fraction of vegetation cover, and biomass is often derived from spectral vegetation index with empirical regression models.Therefore, the bias error in the NDVI caused by the scaling effect will be propagated into such parameters.For instance, a practical case study can be seen in [27,32] for the case of leaf area index and [33] for fraction of vegetation cover.
A key approach to improving our understanding of the scaling effects in the calculation of NDVI has been to analytically clarify the mechanisms underlying the scaling effects.Raffy proved that the error bounds associated with the scaling effects were associated with changes in the spectra at any point within a target pixel (for constant pixel-scale reflectance spectra) [56,57].As a result, the error bounds could be determined using the functions R ∧ and R ∨ , which are obtained using retrieval algorithms and pixel-scale reflectances [56,57].Calculated in this way, the error bounds are sensitive to the degree of nonlinearity in the retrieval algorithm and surface heterogeneity over the target field.Raffy remarked that maximal error bounds could be obtained for fields in which the number of land cover types (land cover heterogeneity) were equal to or fewer than (n + 2), where n represents the dimensions of an input parameter.Their study also examined the scaling effects associated in the calculation of the NDVI [57].Hu et al. investigated the scaling effects in the NDVI using a Taylor series expansion.They concluded that the scaling effects depended on the nonlinearity of the algorithm and the variances of the reflectance spectra within a target pixel [24].However, those previous works did not consider how the dependence of the spectra on the land cover type affected the scaling effects.Jiang et al. addressed the scaling effects in the NDVI using spectral mixture analysis [25].They found that if the values of the level-1 norm of the vegetation and non-vegetation spectra were identical, NDVI was scale-independent, which implicitly agreed with the findings reported in [57].They also stated that the spectral contrast was central to determining the scaling effects of the NDVI [25].Huete et al. conducted a numerical simulation of the dependence of the vegetation and non-vegetation spectra on the NDVI scaling effects [39].However, biases in the area-averaged NDVI as a function of the spatial resolution were never investigated.These biases should be studied to determine how the error bounds are affected by changes in the spatial resolution.Moreover, the spectral features of a land cover type might influence the scaling effects, and this relationship has not previously been investigated analytically.For example, a land cover type can generate either positive or negative biases.
To investigate the behavior of the area-averaged NDVI as a function of the spatial resolution, we developed a theoretical framework [58] that incorporated a linear mixture model (LMM) [59][60][61][62][63] in combination with resolution transform theory.We elucidated the error bounds of the scaling effects of the NDVI by proving the monotonicity of the area-averaged VIs as a function of the spatial resolution under certain conditions.Note that we assumed that the target spectrum consisted of vegetation and non-vegetation endmember spectra that represented pure spectra reflected from homogeneous surfaces under the two-endmember LMM.The point spread function is assumed to be uniform over the xand y-axes, in other words, the sensor response is assumed to be perfect for the purposes of an analytic discussion.

Background
Several studies have implicitly or explicitly examined the monotonic behavior of the area-averaged NDVI over a given area as a function of the spatial resolution [24,25].The monotonicity is central to understanding the error bounds imposed by the scaling effects.Jiang et al. [25] suggested that the area-averaged NDVI should shift monotonically in moving from coarser to finer resolution because the land surface heterogeneity within any given pixel should decrease as the spatial resolution increases.A somewhat controversial conclusion was drawn based on the findings reported by several studies.For example, Hu et al. [24] approximated the differences in the area-averaged NDVI calculated for two cases defined by extreme resolutions using a polynomial parameterized by the reflectance variance and covariance.Their findings predicted that the area-averaged NDVI would not necessarily be monotonic because the variance and covariance of the reflectance are not always monotonic in moving from coarser to finer resolutions.
The non-monotonic behavior of the average NDVI can be easily demonstrated using numerical simulations under two-endmember assumptions (vegetation and non-vegetation endmembers), as shown in Figure 1.Spectral data for a fixed area observed over one pixel (Figure 1  Although the monotonicity of the area-averaged NDVI is central to determining the error bounds, it is not always guaranteed.Hence, monotonic and non-monotonic behaviors are expected to arise in the calculation of an area-averaged NDVI as a function of the spatial resolution.In this study, we try to clarify the mechanisms underlying the scaling effects in terms of the monotonicity of the average NDVI, particularly under conditions in which the average NDVI changes monotonically or non-monotonically.Models of the transformation of the spatial resolution and target spectra will be discussed further in the next section.

Endmember Spectra and Their NDVI
In this study, a target spectrum was modeled using LMM under the constraint that the sum of the endmember spectral weights is unity [62].The endmember spectra for the vegetation and non-vegetation backgrounds are represented by the vectors ρ ρ ρ v = (ρ v,r , ρ v,n ) and ρ ρ ρ s = (ρ s,r , ρ s,n ), where the subscripts v and s represent vegetation and non-vegetation, respectively, and r and n represent the red and NIR bands, respectively.The NDVI value for each endmember is represented by v v and v s , so that v v larger than v s is guaranteed throughout the study, (3)

Resolution Transform Model and Area-Averaged NDVI
To proof the monotonicity, we should focuses on a comparison among the area-averaged NDVI values obtained at different spatial resolutions.For this purpose, the area of a target region was fixed at a given size selected to be larger than or equal to the size of the lowest sensor spatial resolution.The number of pixels contained in a target region depended on the spatial resolution of each sensor, that is, the 'resolution level' was given by the number of pixels within a fixed area, indicated by the index j (j ≥ 1).
The resolution transformation from one level to another was modeled by first defining a simple rule for partitioning pixels, that is, each pixel was partitioned into two at a given time point, as illustrated in Figure 2. We then transformed one spatial resolution into another by applying this partitioning rule to all pixels in the image at the original resolution.A finer resolution could be achieved by repeating this process.Examples of the resolution transform are illustrated in Figure 3.The upper and lower sequences show the resolution transform from resolutions 1 to 4, and 1 to 9, respectively.A fixed target region contained a total of j pixels at resolution j, and the subscript 'k' indicates an individual pixel of a certain resolution level (1 ≤ k ≤ j).The reflectance at each pixel for an image of resolution j, ρ R,j,k and ρ N,j,k for the red and NIR bands, respectively, is then modeled under the normalized (to unity) weight constraints according to Because the target region was divided into j (pixels) at the j-th resolution level, the relative abundance of each endmember must be defined for each pixel k.The endmember abundances corresponding to the vegetation and non-vegetation backgrounds at pixel k under the j-th resolution are indicated by ω j,k and (1 − ω j,k ), as illustrated in Figure 4 (for vegetation fraction).Because the size of the target region is fixed at a given resolution, the total abundance for each endmember species is conserved over the entire region, leading to the following constraints with respect to ω j,k .
or simply, Note that areas of all pixels are equal in these cases, whereas a set of appropriate weights should be used for the area averaging process when the pixel size is not identical.The NDVI can be obtained for each pixel, v j,k , by The variables, parameters, and indices are illustrated in Figure 5 for the two-endmember LMM.Finally, the area-averaged NDVI at the j-th resolution is defined as Figure 4.The endmember abundances of the vegetated surface (fraction of vegetation cover) for pixel k under the j-th resolution is represented by ω j,k .
Figure 5. Illustration of the variables used in this study.v j indicates the area-averaged NDVI at the j-th resolution.v j,k indicates the NDVI calculated for pixel k under the j-th resolution.ρ ρ ρ j,k is the reflectance spectrum of pixel k at a resolution of j and is a linear combination of the vegetation and non-vegetation endmember spectra, ρ ρ ρ v and ρ ρ ρ s with weights ω j,k .

v j Calculated at Various Resolutions
Consider a sequence of images obtained by repeating the partitioning process described in the previous subsection (Figure 2).In such a sequence, an image at resolution j can be expressed in terms of the image at the previous resolution (v j−1 ) and the difference ∆v, where Because ∆v indicates the difference between the results calculated in the forward and reverse steps of the partitioning process applied to one pixel, it can be expressed by where α indicates a fraction of a pixel after the partitioning process (Figure 6).The above equation suggests that, for a resolution sequence, the area-averaged NDVI should change monotonically if the sign of ∆v is the same for all j.This is the case if the number of endmember spectra is limited to two, as further explained in the next section.In this section, we analyze ∆v by focusing on its sign as a function of ω 2,1 , which indicates the FVC for a pixel after the partitioning process.The partial derivative of ∆v = v 2 − v 1 with respect to ω 2,1 is equal to the partial derivative of v 2 because v 1 is independent of ω 2,1 .
The value of the VI obtained by averaging over an area after the partitioning process can be expressed as where the reflectance spectra for the two finer pixels (ρ )) are modeled by a two-endmember LMM, Using Equations ( 13) and ( 14), the VI values for the two pixels obtained after the partitioning process can be expressed as for k = 1, 2.
Because the value of the FVC in the original target pixel (ω 1,1 ) must be conserved after the partitioning process, the following relationships among the FVCs holds: Solving for ω 2,2 and substituting the result into Equation ( 14) yields an expression for ρ ρ ρ 2,2 in terms of Equations ( 13) and (17) indicate that ρ ρ ρ 2,1 and ρ ρ ρ 2,2 can be expressed in terms of a single parameter, ω 2,1 for a given pair of endmember spectra and for a fixed value of ω 1,1 (the total fraction of vegetation over a given target region).Consequently, v 2,1 and v 2,2 can be expressed as functions of ω 2,1 .
Although ω 2,1 is independent of ω 1,1 , its range depends on ω 1,1 such that with the following definitions, It will be useful to calculate the partial derivative of v 2 with respect to ω 2,1 .From Equation ( 12), the derivative can be expressed as a weighted sum of partial derivatives of v 2,1 and v 2,2 .
The partial derivative of v 2,k with respect to ω 2,k becomes where The first term of the right-hand side of Equation ( 20) can be obtained directly from Equation (21).The partial derivative of v 2,2 with respect to ω 2,1 can be obtained from Again, from Equation ( 14) and because ω 2,2 can be written by ω 2,1 , the partial derivative of v 2,2 with respect to ω 2,1 , then, becomes Finally, the partial derivative of v 2 with respect to ω 2,1 becomes where S 1 and S 2 are defined by Note that to prove the monotonicity of the average NDVI, it is important that the sign of Equation ( 27) varies.The first parenthesis in Equation ( 27) can be expressed as The parenthetical equation given above becomes positive, as indicated in following equation derived from Equation (3), This indicates that the sign of Equation ( 30) is necessarily positive.Hence, the sign of ∂v 2 /∂ω 2,1 (Equation ( 27)) depends on the relative differences between ||ρ ρ ρ v || 1 and ||ρ ρ ρ s || 1 , and ω 2,1 .We define the parameter η to discriminate the following cases in this study.
Using the above definition, the sign of the partial derivative behaves as follows: Additionally, when η = 1, then ∂v 2 /∂ω 2,1 = 0, indicating that the area-averaged NDVI does not depend on the spatial resolution [25].
Using Equations ( 33) and (34), the sign of ∆v can be determined by the value of η as follows: Note that the value of ∆v is equal to zero for ω 2,1 = ω 1,1 .Therefore, from Equation (35), the difference between the two average NDVIs can be expressed as These results indicate that the sign of the difference between the area-averaged NDVI values at resolutions 1 and 2 is determined by η, regardless of the spatial parameters such as the endmember abundances ω j,k , or the dividing proportion α under the two-endmember LMM.

The Monotonicity of the NDVI within a Resolution Class
The results of the previous section were applied to determine the spatial resolution at which the average NDVI transformed monotonically with the spatial resolution.
According to Equation (36), which describes the area-averaged NDVI for the 2nd resolution, v 2 is smaller than v 1 , the value before partitioning, if the level-1 norm of the vegetation endmember in a given area exceeds that of the non-vegetation endmember (η > 1).Similarly, v 2 is larger than v 1 if the level-1 norm of the vegetation endmember is less than that of the non-vegetation endmember (η < 1).Therefore, the average NDVI calculated for a resolution sequence generated by the repeated application of the partitioning process is expected to shift monotonically in moving from coarser to finer resolution, as shown in Figure 7.In this figure, for η > 1, v 2 is less than v 1 .When the resolution is transformed from level 2 to level 3, one of the two rectangular pixels (of level 2) is divided into two.The average NDVI after partitioning the pixel, v ′ 2 , is then smaller than the average NDVI before partitioning, v ′ 1 .As a consequence, v 3 is certainly smaller than v 2 .Therefore, the average NDVI decreases monotonically over a resolution sequence as the spatial resolution increases.The resolution sequences define a 'resolution class'.These results are summarized in the following theorems.This condition is sufficient for NDVI to vary monotonically with the resolution.The following theorem relating to the trends in the averaged NDVI can be proven based on Equation (36).
Theorem 2 The trend in the average NDVI over a resolution class (decreasing or increasing) can be determined from η.If η exceeds unity, the average NDVI over a resolution class will decrease monotonically.Similarly, if η is less than unity, the average NDVI will increase monotonically.
Two example resolution classes are illustrated in Figure 8(a,b).These sequences represent different resolution classes, class A and class B, respectively.Figure 8(c) represents a resolution sequence that includes elements of class A and B. v j is the area-averaged NDVI for each resolution level (j).v j is constant or decreases monotonically as a function of the number of pixels within a fixed area in class A and class B (Figure 8(a,b)).On the other hand, v j varies non-monotonically as a function of the spatial resolution along the resolution sequence that contains elements of class A and class B (Figure 8(c)).These elements are not related under application of the partitioning rule in the forward and reverse directions.Therefore, v j does not always vary monotonically if elements from a different resolution class are included in a resolution sequence, although v j certainly varies monotonically within a resolution class.
As a result, VI at the extreme resolutions (the coarsest and finest resolutions) should reach a maximum or minimum because (1) the average NDVI varies monotonically within a resolution class, (2) extreme resolutions belong to the same resolution class, and (3) any type of resolution certainly belongs to a resolution class.These results indicate that the error bounds on the average NDVI associated with a given spatial resolution are determined by the two extreme resolution cases.

Numerical Simulation
Numerical simulations were conducted to verify the monotonic behavior of the area-averaged NDVI.To this end, we considered two areas of the same size that included both vegetation and non-vegetation endmembers, as shown in Figure 9(a,e).The two areas were identical except that the endmember spectra differed.The endmember spectra in each area were determined to follow decreasing or increasing trends for the average NDVI in Figure 9(a,e), respectively.In other words, η was larger than 1 in Figure 9(a) and smaller than 1 in Figure 9(e).The target field was assumed to be acquired at the N -th resolution.Note that the pixel size was significantly smaller than any scene object at this resolution.For the simulation, we generated all possible resolutions of the spectral data based on the original data (at the N -th resolution) by arithmetic averaging using a moving block window under the following conditions: (1) the number of rows and columns in the resolution data set remains constant, and (2) each pixel is an element of a square grid.We define a set B that includes elements of all divisors of n (the square root of N ), l, B = {l; l ∈ Z Z Z + , l|n}.(37) where Z Z Z + represents the set of positive integers.Then all possible members of the N -th resolution under the limited conditions are represented by a set A N as Members of the resolution class are selected from the set A N in the simulations.To define a general representation of the set of members in a resolution class (C ⊂ A N ), the members of C are indicated as where z ∈ Z Z Z + and z ̸ = 1.
The monotonic trend within the resolution class was verified by selecting two examples of the resolution class that satisfied Equation ( 39), containing five elements.For example, the subsets C 1 ⊂ A 40000 and C 2 ⊂ A 40000 were assumed to be C 2 = {1, 25, 625, 2500, 10000, 40000}.
The subsets (resolution class) C 1 and C 2 can be referred to as 'class-1' and 'class-2', respectively, as illustrated in Figure 10.According to the findings discussed in the previous section, the area-averaged NDVI did not vary monotonically over the resolution sequence comprising set A 40000 .On the other hand, the area-averaged NDVI varied monotonically over the resolution class sets (C 1 and C 2 ). Figure 9(b,f) shows that the average NDVI varied non-monotonically if the resolution sequence included two classes (class-1 and class-2).The average NDVI certainly behaved monotonically within a single class.The average NDVI within a single class (indicated by the solid line) decreased monotonically, as plotted in Figure 9(c,d) for class-1.The same observation held for the increasing case (class-2) shown in Figure 9(g,h) for class-2.

Discussion
A resolution transform model with a certain partitioning rule was proposed to analyze the area-averaged NDVI for a fixed area described under different spatial resolutions.To elucidate the error bounds on the calculated NDVI, the resolution class (condition for monotonicity) was introduced.The resolution class indicates the values that are interconnected under the forward and reverse application of the partitioning rule.Within a class, the area-averaged NDVI certainly vary monotonically with the spatial resolution.The trends in the average NDVI (decreasing or increasing) within a resolution class can be determined based on the factor η calculated for a given NDVI, in agreement with the scaling effects predicted in a numerical study [39].As a consequence, the values of the VI at extreme resolutions (at the coarsest and finest resolutions) should reach maximum or minimum values because (1) the average NDVI varied monotonically within a resolution class, (2) extreme resolutions belong to the same resolution class, and (3) any resolution certainly belongs to one resolution class.These results indicate that the error bounds on the inherent uncertainties due to the scaling effects can be specified by the maximum and minimum values.
A major limitation of this work is the number of endmember spectra assumed in the LMM.This limitation may cause a loss of practicality in some extent, since in general one encounters more numbers of distinct surfaces within a satellite imagery.This assumption, however, is reasonable at this stage of investigation for the following reason.The monotonicity of NDVI in the framework of scaling issue has not been fully and thoroughly investigated to date.As an initial step of tackling the theme, it is appropriate to start with the simplest case.Moreover, this assumption enables us to analyze the influence of spatial resolution analytically from the beginning to the end.And, owing to this simplicity, one can even reach the understanding of deep insight of the complex phenomena with no ambiguity.Therefore, the assumption brings us both advantage and disadvantage.
The assumption is indeed not practical knowing the fact that pixel size of sensors with middle to lower resolution (with wide field-of-view) is most likely large enough to include three or more numbers of endmember species.Nevertheless, considering the trend of sensor technology, spatial resolution of future sensor will tend to be higher.As a consequence, the pixel size eventually becomes fine enough to the level such that region of interest can be considered as a two-endmember case.In this sense, the results and findings from this study would serve as theory that reasonably explains behavior of scaling effect with some degrees of practicality.

Conclusions
We analytically investigated the mechanisms underlying the scaling effects present in the calculation of NDVI under a two-endmember LMM, and we related the scaling effects to the monotonic behavior as a function of the spatial resolution.The scaling effects depended on the spectral features, and this relationship was clarified to show that the trend (increasing or decreasing) depended on the vegetation and non-vegetation endmember spectra.The proof of the monotonicity indicates that the error bounds may be derived deterministically using the area-averaged values calculated from the data at the extreme resolutions.
This study establishes a theoretical framework for describing the scaling effects in the calculation of biophysical or climatological parameters retrieved from remotely sensed data.In this framework, remotely sensed data is modeled using a linear mixture of endmember spectra.The applicability of the results and the findings of this study are mostly restricted by the assumption of two endmembers.Further investigations are required to expand the discussion to (1) the analysis of a greater number of endmembers, and (2) the analysis of other biophysical parameters, such as LAI or FAPAR.Nevertheless, the fundamental behavior of the scaling effects form a theoretical basis for similar investigations.
(a)), four pixels (Figure1(b)), or nine pixels (Figure1(c)) were prepared.The NDVI values for each pixel at each resolution were obtained, and the area-averaged NDVI was calculated at each resolution.The average NDVI values as a function of the spatial resolution are plotted in Figure1(d).It is clear that the average NDVI based on two endmember LMMs did not always vary monotonically as a function of the spatial resolution.

Figure 1 .
Figure1.Non-monotonic behavior on an area-averaged NDVI as a function of the number of pixels (spatial resolution) over a fixed area.(a), (b), or (c) illustrate the hypothetical image obtained for one, four, or nine pixels within a fixed area that includes vegetation and soil backgrounds.The vegetation spectrum (covering both red and NIR reflectance spectra) is assumed to be (0.05, 0.4), and the soil background spectrum is (0.15, 0.15).(d) The area-averaged NDVI as a function of the number of pixels.

Figure 2 .
Figure 2. Elements of the resolution transformation process based on a simple partitioning procedure.

Figure 3 .
Figure 3. Illustration of the resolution transform using the partitioning rule.The upper panel shows an example of the resolution transform from levels 1 to 4, and the lower panel shows the transform from levels 1 to 9.

Figure 6 .
Figure 6.The fraction of a pixel after the partitioning process, represented by α.

Figure 7 .Theorem 1
Figure 7. Monotonicity of the area-averaged NDVI in consecutive resolution sequences.A trend increases or decreases depending on the factor η.

Figure 8 .
Figure 8. Illustration of the resolution transform from resolution levels 1 to 4. (a) and (b) describe the elements of the resolution sequences class A and class B, respectively.(c) depicts the elements of a resolution sequence in which elements from both class A and class B are included.The average NDVI varies monotonically in (a) and (b).In contrast, non-monotonic behavior is observed in (c) at the resolution level 3.

Figure 9 .
Figure 9. Numerical demonstration of the scaling effects in the calculation of the area-averaged NDVI.(a) Hypothetical field over which was observed a vegetation endmember spectrum of (ρ r , ρ n ) = (0.05, 0.4) and a non-vegetation endmember spectrum of (0.15,0.15).(b-d) present the average NDVI calculated from the field (a).(b) Variations in the average NDVI, including the two resolution classes (set A 40000 ); (c) Variations within class-1 (subset C 1 ) only, or (d) for class-2 (subset C 2 ) only.(e) A hypothetical field over which was observed a vegetation endmember spectrum of (0.05,0.4) and a non-vegetation endmember spectrum of (0.3,0.3).(f-h) present the average NDVI calculated from the field (e).(f) Variations in the average NDVI over the two resolution classes; (g) Variations within class-1 only, or (h) for class-2 only.

Figure 10 .
Figure 10.Two resolution classes used in the experiments.The upper sequence (class-1, C 1 ) presents a different class than the lower sequence (class-2, C 2 ).