Vegetation isoline equations with first- and second-order interaction terms for modeling a canopy-soil system of layers in the red and near-infrared reflectance space

Munenori Miura; Kenta Obata; Hiroki Yoshioka

doi:10.1117/1.JRS.9.095987

18 November 2015 Vegetation isoline equations with first- and second-order interaction terms for modeling a canopy-soil system of layers in the red and near-infrared reflectance space

Munenori Miura; Kenta Obata; Hiroki Yoshioka

Author Affiliations +

J. of Applied Remote Sensing, 9(1), 095987 (2015). doi:10.1117/1.JRS.9.095987

Abstract

A previously proposed vegetation isoline equation suffers from errors if the soil background of a canopy layer is bright. These errors arise from the truncation of the second- and higher-order terms that represent photon interactions between the canopy and the soil. An isoline equation that includes a second-order interaction term is introduced. The equation was initially derived by explicitly including a second-order interaction term in both the red and near-infrared (NIR) reflectance spectra (symmetric approximation). We also examined an alternative model in which the interaction term was included only in the NIR band (asymmetric approximation). In this model, the derived isolines tend to shift upward (overcorrection effects). Numerical experiments revealed that the errors in the isoline obtained by the asymmetric approximation were reduced in magnitude to nearly one-fifth of the errors in the previously proposed method. Its accuracy was higher than that of the symmetric approximation, despite the fact that the asymmetric approximation included fewer terms than the symmetric approximation. The improved model accuracy resulted from the overcorrection effects, which compensated for the truncation error. With the simplicity and improved accuracy, the current isoline equations provide a good alternative to the previous approach.

Miura, Obata, and Yoshioka: Vegetation isoline equations with first- and second-order interaction terms for modeling a canopy-soil system of layers in the red and near-infrared reflectance space

1. Introduction

Biophysical parameter retrieval from remotely sensed reflectance spectra is a fundamental goal in the field of land remote sensing. Qin et al.¹ categorized the available retrieval algorithms into four groups based on the approaches taken: (1) techniques that relied on a spectral vegetation index (VI) and its correlation with biophysical parameters, such as the leaf area index (LAI);²3.4.5.6.^–⁷ (2) algorithms that used lookup tables;⁸^,⁹ (3) neural networks;¹⁰11.^–¹² and (4) direct inversions of numerical models [e.g., models of radiative transfer (RT)] using optimization methods.¹^,⁹^,¹³14.^–¹⁵ These approaches present advantages and disadvantages over other approaches in terms of accuracy, computational costs, complexity, and applicability. The common feature of all categories of approach is that a better model increases the accuracy of the retrieved parameters. For this reason, significant efforts have been applied toward improving the accuracy of physical and numerical models.

The concept of vegetation isolines²^,¹⁶17.18.19.^–²⁰ forms the basis of the spectral VI,²^,⁷^,¹⁸^,²⁰ which has been widely used as a proximity measure (model) of surface biophysical parameters.²¹ The isoline concept has been used as an analytical tool for investigating the influence of the soil on the retrieved parameters.²²23.24.25.26.27.^–²⁸ From this standpoint, vegetation isoline equations provide a model for the relationship between reflectances at different wavelengths. Several reports have attempted to use isolines in the analysis of the VI ²⁴^,²⁸ and in parameter retrieval.²⁹^,³⁰ Recently, these relationships were applied to intercalibration studies of the VI values obtained from different sensors.³¹32.^–³³ Thus, the accuracy of the vegetation isoline equation must be improved in order to improve parameter retrieval algorithms and better understand the factors that affect intercalibration studies.

Yoshioka et al.³⁴ found that the derived isoline by the isoline equations loses accuracy in intermediate ranges of the LAI. Loss of accuracy in the isoline equation arises from truncations of the terms that correspond to multiple interactions among the photons reflected from the canopy layer (at the bottom surface) and the soil surface. These terms are referred to as higher-order interaction terms³⁴ in this study. The isoline equations derived by Yoshioka et al. ²⁴^,³⁴^,³⁵ retain the interaction terms up to the first-order terms. The truncation, however, simplifies the derivation and yields a final form that is useful as an analytical tool.²⁴^,²⁸^,³¹32.^–³³ One drawback of this truncation, however, is the loss of accuracy, which must be improved while simultaneously retaining the simplicity of the model. This study seeks to do just this.

We conducted a series of pilot studies³⁶37.^–³⁸ to explore possible improvements to the higher-order terms in the isoline equations. Derivations were developed for use in several cases. Three model issues have yet to be clarified. First, the relationship between the previously derived isoline, which includes a first-order interaction term, and the newly derived isoline, which includes higher-order interactions, is not yet understood. Second, the mechanism by which the errors were reduced upon inclusion of the higher-order term has not yet been identified. In some cases, isolines containing fewer interaction terms to describe one of the two bands showed significantly better accuracy. This mechanism must be explored systematically. This study examines these matters from an analytical and numerical perspective. The isoline equation derived here is even simpler than the equation introduced in the pilot studies; thus, this equation may be readily applied to new analyses. Finally, numerical procedures for determining the isoline parameters were not discussed in the pilot studies. These matters require further study for application purposes.

The objective of this study is to improve the first-order isoline equation by including second-order interaction terms while maintaining the model simplicity. This objective was achieved through a novel approach. Instead of retaining the second-order interaction terms of both the red and near-infrared (NIR) bands, we retained the term only for the reflectance of the NIR band. The asymmetric treatment of the second-order interaction significantly improved the model accuracy without sacrificing the simplicity of the derived expression. This study describes the formal steps, which used to derive the improved version of the vegetation isoline equation and validate its accuracy by conducting numerical experiments based on a coupled leaf and canopy RT model, PROSAIL.³⁹

The remainder of this discussion presents a review of a previously derived vegetation isoline equation,²⁴^,³⁴ referred to as a first-order isoline equation in this study. Two forms of higher-order approximations are derived. The results of numerical experiments using an RT model are then shown to evaluate the degree to which the accuracy was improved by the introduction of the higher-order approximations. The mechanism by which the accuracy was improved by the approximations is discussed in detail. Finally, the findings of this study are summarized.

2. Analytical Model and First-Order Vegetation Isoline

2.1.

Analytical Canopy Reflectance Model

This study begins with an analytical form of the top of the canopy (TOC) reflectance model, which has been used extensively in this field of study:⁴⁰

(1)

ρ_{λ} = ω ρ_{v λ} + ω \frac{T_{λ}^{2} R_{s λ}}{1 - R_{s λ} R_{v λ}} + (1 - ω) R_{s λ},

where the variables and notations mainly follow those given in Refs. 24 and 34, except for the fraction of vegetation cover (FVC) represented by

ω

in this study. The variable

ρ_{λ}

represents the TOC reflectance at the wavelength

λ

,

ρ_{v λ}

represents the “pure” canopy reflectance, which can only be obtained by assuming perfect absorbance beneath the canopy layer,

R_{s λ}

and

R_{v_{λ}}

represent the bihemispherical reflectance of the soil and canopy layers, respectively. The variable

R_{v λ}

is somewhat special among these variables because it represents the albedo of the bottom surface of the canopy layer. The details of the model are illustrated in Fig. 2 of Ref. 34.

T_{λ}^{2}

represents the two-way transmittance of the canopy layer at the indicated wavelength,

λ

. We next defined the area-averaged two-way transmittance at a point:

(2)

\bar{T_{λ}^{2}} = ω T_{λ}^{2} + 1 - ω .

The word “area-averaged” indicates that the bare soil region,

(1 - ω)

, provides “perfect” two-way transmittance properties, and the averaged transmittance of the partially vegetated area may be modeled using a weighted average of

T_{v}^{2}

and 1.0, with weighted values of

ω

and

(1 - ω)

. In the following section, we briefly explain the vegetation isoline equation introduced previously.²⁴^,³⁴

2.2.

Vegetation Isoline Equation with a First-Order Interaction Term in the Red–Near-Infrared Reflectance Subspace

The isoline equation introduced in this subsection is referred to as the first-order isoline because only a single interaction term with the soil surface is included.

The model was derived by explicitly including the first-order interaction term as the second term of the right-handside (RHS) of Eq. (3),

(3)

ρ_{λ} = ω ρ_{v λ} + ω T_{λ}^{2} R_{s λ} + \frac{ω T_{λ}^{2} R_{s λ}^{2} R_{v λ}}{1 - R_{s λ} R_{v λ}} + (1 - ω) R_{s λ},

where the second term of the RHS in Eq. (3) indicates the first-order interaction term. The second and higher interaction terms, represented by the third term of the RHS, are further defined as a truncated-order term:

(4)

O (R_{s λ}^{2}) = \frac{T_{λ}^{2} R_{s λ}^{2} R_{v λ}}{1 - R_{s λ} R_{v λ}} .

Equation (3) can be simplified using Eq. (2):

(5)

ρ_{λ} = ω ρ_{v λ} + \bar{T_{λ}^{2}} R_{s λ} + ω O (R_{s λ}^{2}) .

Two equations were used to describe the red (denoted by the subscript $R$ ) and NIR wavelengths ( $N$ ):

(6)

ρ_{R} = ω ρ_{v R} + \bar{T_{R}^{2}} R_{s R} + ω O (R_{s R}^{2}),

(7)

ρ_{N} = ω ρ_{v N} + \bar{T_{N}^{2}} R_{s N} + ω O (R_{s N}^{2}) .

The soil line assumption of Ref. 41 was applied:

(8)

R_{s N} = a R_{s R} + b,

where

a

and

b

represent the slope and offset of the soil line, respectively. Equations (6)–(8) were used to eliminate the soil reflectances

R_{s R}

and

R_{s N}

to obtain the first-order approximated vegetation isoline equation:

(9)

ρ_{N} = a γ_{1} ρ_{R} + D_{1} + ε_{1},

where

γ

,

D

, and

ε_{1}

are defined by

(10)

γ_{1} = \bar{T_{N}^{2}} / \bar{T_{R}^{2}},

(11)

D_{1} = b \bar{T_{N}^{2}} + ω (ρ_{v N} - a γ_{1} ρ_{v R}),

(12)

ε_{1} = ω [O (R_{s N}^{2}) - a γ_{1} O (R_{s R}^{2})] .

Finally, we obtained an approximated form of the vegetation isoline equation by truncating the higher-order interaction term:

(13)

ρ_{N} \approx a γ_{1} ρ_{R} + D_{1} .

The first-order isoline model suffers from truncation errors. For example, Yoshioka et al.³⁴ indicated that the truncation error increases at higher soil reflectance values, and the relative error can reach 5%. Their findings suggest that the errors in the retrieved biophysical parameters, calculated based on the isoline formula, can reach the same magnitude. The truncation errors in the isoline formula should, therefore, be minimized to obtain more accurate parameter retrievals using the isolines.

The truncation error pattern obtained from the first-order isoline model was numerically characterized using a RT model PROSAIL.³⁹ Figure 1 shows the first-order isoline and reflectance spectra simulated using the model in the red and NIR reflectance subspaces. The empty circles denote the simulated reflectance, which was considered to be a “true” spectrum. The solid lines indicate the first-order isoline at various LAI values of the pure canopy component ranging from 0 to 4 at intervals of 0.5. The left and right figures present the results obtained from different FVC ( $ω$ ) values: full coverage ( $ω = 1$ ) or half coverage ( $ω = 0.5$ ), respectively. A detailed description of the simulation conditions is provided in the latter section. These results confirmed that the error (the discrepancy between the empty circles and the solid lines) in the first-order isoline could be reproduced numerically using the model presented in this study.

Fig. 1

First-order vegetation isolines and simulated reflectance spectra obtained using PROSAIL. The value of the fraction of vegetation cover (FVC) represented by $ω$ was assumed to be (a) unity, representing a fully covered case and (b) 0.5, representing a partially covered case.

The error trend was characterized by plotting the distance between the isolines and the true spectra, as shown in Fig. 2, at $LAI = 1.0$ (solid line) and $LAI = 4.0$ (dashed line) for both the fully covered and partially covered cases. The error was plotted as a function of the soil reflectance in the red band. The figure clearly reveals two characteristics: first, the error increased as the soil became brighter. Second, the error at $LAI = 1.0$ always exceeded that obtained for $LAI = 4.0$ . Because these trends have been analyzed previously,²⁴^,³⁴ we summarize the conclusions briefly: (1) the influence of the higher-order interaction terms increased for brighter soils and (2) this influence first increased and then decreased as the canopy thickened.

Fig. 2

Distance between the first-order isoline and the true reflectance spectrum shown in Fig. 1 for (a) $ω = 1.0$ and (b) $ω = 0.5$ .

The error could be reduced to some extent by adjusting the higher-order terms. The next section discusses the derivation steps that are used to obtain the second-order approximated reflectance and a new isoline equation referred to as the “asymmetric order isoline” in this study.

3. Second-Order Approximation and Asymmetric Order Isoline

3.1.

Parametric Form of the Second-Order Approximated Reflectance Spectra

The second-order interaction term could be explicitly separated from the higher-order term in Eq. (3) using the approach taken to separate the first-order term. Specifically, the interaction terms were retained up to the second order for both the red and NIR spectra in Eq. (3). The resulting equation then becomes

(14)

ρ_{λ} = ω ρ_{v λ} + \bar{T_{λ}^{2}} R_{s λ} + ω T_{λ}^{2} R_{v λ} R_{s λ}^{2} + ω O (R_{s λ}^{3}),

where

(15)

O (R_{s λ}^{3}) = \frac{T_{λ}^{2} R_{s λ}^{3} R_{v λ}^{2}}{1 - R_{s λ} R_{v λ}} .

We next rewrote Eq. (14) to describe the red and NIR bands using the subscripts $R$ and $N$ to obtain

(16)

ρ_{R} = ω ρ_{v R} + \bar{T_{R}^{2}} R_{s R} + ω T_{R}^{2} R_{v R} R_{s R}^{2} + ω O (R_{s R}^{3}),

(17)

ρ_{N} = ω ρ_{v N} + \bar{T_{N}^{2}} R_{s N} + ω T_{N}^{2} R_{v N} R_{s N}^{2} + ω O (R_{s N}^{3}) .

The soil line in Eq. (8) was used to obtain the reflectance spectrum to a second-order approximation, as represented by the following form with $R_{s R}$ as a parasite parameter:

(18)

[\begin{matrix} ρ_{R} \\ ρ_{N} \end{matrix}] = [\begin{matrix} \bar{T_{R}^{2}} & ω T_{R}^{2} R_{v R} \\ a \bar{T_{N}^{2}} + ω 2 a b T_{N}^{2} R_{v N} & ω a^{2} T_{N}^{2} R_{v N} \end{matrix}] [\begin{matrix} R_{s R} \\ R_{s R}^{2} \end{matrix}] + ω [\begin{matrix} ρ_{v R} \\ ρ_{v N} + b \bar{T_{N}^{2}} / ω + b^{2} T_{N}^{2} R_{v N} \end{matrix}] + ω [\begin{matrix} O (R_{s R}^{3}) \\ O (R_{s N}^{3}) \end{matrix}] .

The last term could be neglected to obtain a parametric representation of the approximated second-order spectrum. This form was used only in the numerical experiments to evaluate the isoline equation, as discussed in Sec. 3.2.

3.2.

Isoline Equation Obtained by Asymmetrically Truncating the Higher-Order Interactions

The previous subsection discussed the retention of the higher-order interaction terms up to the second-order terms for both the red and NIR bands. In this subsection, we include the second-order interaction term only in the description of the NIR band, and the red band is approximated up to the first-order interaction term.

A system of equations was obtained using Eqs. (6), (8), and (17), in which the soil reflectances were eliminated. Algebraic manipulations yielded the final results:

(19)

ρ_{N} = a^{2} ζ ρ_{R}^{2} + a (γ_{1} + δ_{1}) ρ_{R} + D_{1} + δ_{0} + ε_{2},

where

(20)

ζ = ω T_{N}^{2} R_{v N} / {(\bar{T_{R}^{2}})}^{2},

(21)

δ_{0} = ζ {(b \bar{T_{R}^{2}} - ω a ρ_{v R})}^{2},

(22)

δ_{1} = 2 ζ (b \bar{T_{R}^{2}} - ω a ρ_{v R}),

(23)

ε_{2} = ω O (R_{s N}^{3}) + a^{2} ω^{2} ζ {[O (R_{s R}^{2})]}^{2} - a ω [2 a ζ ρ_{R} + γ_{1} + 2 ζ (b \bar{T_{R}^{2}} - a ω ρ_{v R})] O (R_{s R}^{2}) .

Neglecting $ε_{2}$ from Eq. (19), we have

(24)

ρ_{N} \approx a^{2} ζ ρ_{R}^{2} + a γ_{2} ρ_{R} + D_{2},

where

(25)

γ_{2} = γ_{1} + δ_{1},

(26)

D_{2} = D_{1} + δ_{0} .

4. Results of the Numerical Simulations

4.1.

Parameter Settings Used in the Numerical Experiments

A series of numerical experiments was conducted using the canopy RT model PROSAIL.³⁹ This model consists of the leaf model PROSPECT⁴² and the canopy model SAIL;⁴³ thus, two types of input parameter were required. Numerical experiments were conducted using the set of input values provided with the code, except that the three input parameters—LAI, leaf angle distribution (LAD), and soil factor—were set as follows. The soil factor was obtained from the mixture ratio of the wet and dry soil spectra provided with the code. The parameter ranges of the three parameters were as follows. LAI was varied from 0.0 to 4.0 at 0.5 intervals (9 levels), and the soil factor was varied from 0.0 to 1.0 at 0.1 intervals (11 levels). During the numerical experiments, six LAD models (planophile, erectophile, plagiophile, extremophile, spherical, and uniform distributions) were assumed. (This section focuses on the case of a spherical LAD, as a representative case, for brevity.) Finally, the obtained reflectance spectra were linearly mixed with the pure soil reflectance spectra using the fraction of green cover ( $ω$ ) as a weight. The parameter $ω$ was varied from 0.0 to 1.0 in 0.1 intervals (11 levels). The total number of spectra used to model each LAD was 1089. In the analysis, we assumed that the reflectances at 655 and 865 nm provided representative values of the red and NIR bands in this study. This choice of wavelength pair corresponded to the center of the red and NIR bands obtained from the Landsat 8 operational land imager sensor.

The other input parameters provided with the code were fixed as follows. For the SAIL part of the code, the parameter describing the hot spot (hspot) was set to 0.01. The solar zenith, observation zenith, and relative azimuth angle were set to 30, 10, and 0 deg, respectively. For the PROSPECT part of the code, chlorophyll-a and -b, carotenoid, and the leaf mass per area were assumed to be 40, 8, and 0.009 in $g / {cm}^{2}$ , respectively. The equivalent water thickness was set to 0.01 cm, and the brown pigment content was assumed to be zero. Finally, the leaf mesophyll structure ( $N$ ) was assumed to be 1.5 (the equivalent number of layers).

4.2.

Numerical Procedure Used for the Isoline Parameter Retrieval and the Error Estimation

The isoline parameters used in Eqs. (9) and (19) were computed from $T_{v λ}^{2}$ and $ρ_{v λ}$ , which were determined using the algorithm described in Ref. 34. This algorithm required two hypothetical simulations in which the soil was assumed to be “spectrally flat,” with a zero reflectance value, or the soil was assumed to have a medium reflectance value over the entire wavelength range. In addition to $T_{v λ}^{2}$ and $ρ_{v λ}$ , the asymmetric order approximated isoline, Eq. (19), required a value of $R_{v N}$ to define $ζ$ . $R_{v N}$ was determined by conducting an additional simulation in which the soil was assumed to be spectrally flat and even brighter than was assumed in the simulation used to determine $T_{v λ}^{2}$ . The assumption of brighter soil increased the photon contributions of the higher-order interactions. The parameter $R_{v N}$ was determined by solving Eq. (17) for $R_{v N}$ and using $T_{v N}^{2}$ and $ρ_{v N}$ , which were computed prior to $R_{v N}$ . With these variables in hand, the isoline parameters $γ_{1}$ , $D_{1}$ , $γ_{2}$ , $D_{2}$ , $ζ$ , $δ_{0}$ , and $δ_{1}$ were obtained. In summary, three hypothetical simulations were conducted to determine the isoline parameters that corresponded to spectrally flat soil at three different brightness levels.

The errors in the isolines and the approximated reflectance spectra were estimated by computing the distance from the true spectra. Note that the error in the isolines should be equal to the distance between the true reflectance spectrum (which includes all higher-order interaction terms) and the “isoline” represented by Eqs. (9) and (19), corresponding to the first-order and asymmetric order isolines, respectively. Even if a spectrum approximated using a model based on truncated higher-order terms was far from the true spectrum, the error could be zero provided that the “isoline” passed through the point of the true spectrum in the red–NIR reflectance subspace. This distance was employed as a measure of the error because the error in the biophysical parameter retrieval obtained using the isoline reached zero numerically under conditions in which the true spectrum point coincided with the isoline. The goal of isoline determination is to identify the conditions under which the isoline coincides with the true spectra.

4.3.

Comparison of the Accuracy Across the Three Approximations

The performances of the two types of isoline, namely, the first-order and second-order isolines, were compared with the true and second-order reflectance spectra obtained in the red–NIR reflectance subspace. Figure 3 represents the isolines as solid lines. The first-order [Eq. (9)] and the asymmetric order [Eq. (19)] isolines are plotted as black and red lines, respectively. The empty circle represents the true spectra obtained directly from PROSAIL and corresponds to the spectra represented by Eq. (1). The crossmark was used to indicate the second-order approximated reflectance spectra [Eq. (18)]. The isolines were compared at two values of the vegetation cover: (a) $ω = 1.0$ and (b) $ω = 0.5$ .

Fig. 3

First-order isoline (black line) defined by Eq. (9), and the asymmetric isoline (red line) defined by Eq. (19), with the second-order reflectance spectra (cross mark) computed using Eq. (18) and the true spectra (empty circle) computed using PROSAIL. The value of the FVC represented by $ω$ was assumed to be (a) unity, representing a fully covered case and (b) 0.5, representing a partially covered case.

These figures reveal that the asymmetric isolines (red lines) provided much better approximates for the true spectra (empty circles) than the first-order isolines (black lines) over the entire LAI range. The error in the asymmetric order isoline was smaller than that obtained from the first-order isoline. Further analyses in this subsection (Table 1) and the next subsection (Fig. 6) indicated that the asymmetric isolines were even closer to the true spectra than the second-order approximated reflectance spectra (crossmarks) over the full spectral range. This result is, to some extent, surprising because the error in the second-order approximated spectrum was expected to be smaller than the error obtained from the asymmetric isoline, in which one of the bands (the red band in this study) was approximated to the first order instead of to the second order. This result will be further discussed later in this section.

Table 1

Average, standard deviation, and maximum differences between the true reflectance spectra and the three cases of the isoline/reflectance spectra. The differences were computed by assuming six types of leaf angle distribution (LAD): planophile, erectophile, plagiophile, extremophile, spherical, and uniform distributions.

	First order isoline	Second order spectrum	Asymmetric isoline
	LAD: planophile
Mean	1.93E−03	8.53E−04	3.46E−04
STD	2.47E−03	1.33E−03	5.06E−04
MAX	1.30E−02	8.08E−03	2.56E−03
	LAD: erectophile
Mean	2.93E−03	1.92E−03	8.44E−04
STD	3.68E−03	2.84E−03	1.16E−03
MAX	1.85E−02	1.60E−02	5.79E−03
	LAD: plagiophile
Mean	1.57E−03	5.53E−04	2.16E−04
STD	2.06E−03	8.89E−04	3.23E−04
MAX	1.10E−02	5.60E−03	1.62E−03
	LAD: extremophile
Mean	1.74E−03	6.33E−04	2.47E−04
STD	2.26E−03	1.01E−03	3.68E−04
MAX	1.18E−02	6.22E−03	1.84E−03
	LAD: spherical
Mean	1.95E−03	8.79E−04	3.57E−04
STD	2.51E−03	1.37E−03	5.21E−04
MAX	1.32E−02	8.31E−03	2.65E−03
	LAD: Uniform
Mean	1.65E−03	5.85E−04	2.28E−04
STD	2.15E−03	9.38E−04	3.41E−04
MAX	1.13E−02	5.86E−03	1.71E−03

The errors obtained from the two isolines and the second-order reflectance spectra were directly compared, and the distance between these isolines and the spectra obtained from the true values (including all higher-order terms) are plotted as a function of the soil reflectance $R_{s R}$ for the four combinations of LAI and $ω$ (Fig. 4). The errors obtained from the first-order isoline, second-order reflectance, and asymmetric order isoline are denoted using different colors. Figures 4(a) and 4(b) show the results obtained for the fully covered case, with $LAI = 1.0$ and 4.0, respectively. Figures 4(c) and 4(d) show the results obtained for the half-covered case. These figures indicated that (1) the error over $R_{s R}$ was high at high values of $R_{s R}$ and (2) in most cases, the errors of the asymmetric isoline were the smallest among the three models over the range of $R_{s R}$ . This result indicates that the accuracy of the asymmetric isoline dramatically improved over the entire range of $R_{s R}$ . Again, recall that the asymmetric isoline used a first-order approximated reflectance to model the red band, whereas a second-order reflectance was used to model the second-order interaction terms in both the red and NIR bands. The error in the second-order reflectance (blue line) was expected to be smaller than the error obtained from the asymmetric isoline (red line). These results could be explained in terms of the relationship between the overcorrection and the truncation error. This mechanism is described in detail below.

Fig. 4

Distance between the first-order isoline, second-order reflectance spectrum, or the asymmetric isoline and the corresponding true reflectance spectra, as a measure of the error shown in Fig. 3 at four combinations of leaf area index (LAI) and $ω$ : (a) $LAI = 1.0$ ; $ω = 1.0$ , (b) $LAI = 4.0$ ; $ω = 1.0$ , (c) $LAI = 1.0$ ; $ω = 0.5$ , and (d) $LAI = 4.0$ ; $ω = 0.5$ .

Before discussing these surprising results, we will analyze the error trend by investigating the error distributions obtained from the simulations conducted using combinations of the three input parameters (LAI, FVC $ω$ , and soil reflectance $R_{s R}$ ) employed in this study. Figure 5 shows a histogram of the errors obtained from the three approximations. The figure reveals that the asymmetric order isoline errors were clustered at lower errors, unlike the errors of the other two approximations, and very few simulations provided errors of $1.0 E - 3$ . On the other hand, the errors of the first- and second-order isoline reflectance simulations were uniformly distributed at distances exceeding $1.0 E - 3$ . These results confirmed that the asymmetric order isoline outperformed the other two approximations.

Fig. 5

Histogram of the errors obtained from the first-order isoline, second order reflectance, and asymmetric order isoline. The total number of simulated cases was 1089 = (9 discrete values of LAI) $\times$ (11 values of FVC, $ω$ ) $\times$ (11 values of soil reflectance $R_{s R}$ ), assuming a spherical leaf angle distribution.

The model performances were further validated by varying the LAD. The experiments compared the results obtained by assuming six different LAD models: planophile, erectophile, plagiophile, extremophile, spherical, and uniform distributions. The error values were averaged over the entire range of the three parameters (LAI, $ω$ , and $R_{s R}$ ). These results are summarized in Table 1, which lists the standard deviation and the maximum values obtained in all cases. The table indicates that the average error of the asymmetric order isoline was much smaller (by nearly one order of magnitude) than the average error of the first-order isoline for all LAD cases. The table also reveals that the performance of the asymmetric order isoline was better than that of the second-order reflectance in terms of the average, standard deviation, and the maximum value. Table 1 supports the above findings (e.g., that the asymmetric order isoline performed better than the other models), derived in Figs. 4 and 5. Section 4.4 analyzes the detailed mechanisms underlying this trend.

4.4.

Error Reduction Mechanisms Using the Second-Order Isoline

The numerical results presented above indicate that the asymmetric order isoline model was more accurate than the reflectance approximated up to the second-order interaction terms, despite the fact that the latter included a greater number of terms in the red band than did the asymmetric order isoline. This trend is not easy to understand intuitively, but it appeared to result from the tendency of the model to overcorrect by a degree that was approximately equal to the order of magnitude of the truncation error in the isoline.

We further examined the model performances by plotting the four reflectance spectra: (1) the true spectrum computed directly using PROSAIL; (2) the first-order approximated reflectance; (3) the second-order approximated reflectance; and (4) the asymmetric order approximated reflectance using two types of vegetation isoline (first-order isoline and asymmetric order isoline). Figure 6 shows a limited region of the red-NIR reflectance subspace to illustrate the differences between the four predicted spectra. This part of the subspace corresponds to the results obtained at $LAI = 2.0$ for the case of full canopy coverage.

Fig. 6

Comparison of the four approximated reflectance spectra: true spectrum calculated using PROSAIL (circle), first-order approximated spectrum (square), second-order approximated spectrum (cross), and asymmetric order approximated spectrum (triangle), with the two isolines (first-order, black line, and asymmetric order, red line.) The value of $LAI = 2.0$ and full canopy coverage ( $ω = 1.0$ ) was assumed.

Figure 6 reveals that the closest spectrum to the true spectrum (denoted by the circle) is the second-order approximated reflectance, denoted by the crosses. The asymmetric isoline represented by the red line is even closer to the true spectrum, whereas the asymmetric order approximated reflectance spectrum (triangle) is further from the true spectrum compared to the second-order approximated spectrum. Recall that the errors are measured as the distance between the true spectrum (circle) and the model spectra. Because the distance to the isoline (red line) was smaller than the distance between the true spectrum and the second-order approximated spectrum (cross), the error of the asymmetric isoline was surely smaller. Therefore, these results are consistent with the trend described in the previous subsection.

This trend could be understood as resulting from an overcorrection to the NIR band in the asymmetric approximation and the truncation error in the isoline. In the asymmetric case, the inclusion of a second-order term only in the NIR band overcorrected the spectrum upward in the direction from the first-order spectrum (square), as illustrated in Fig. 7. The overcorrection of the NIR band shifted the asymmetric order approximated spectrum (triangle) from the position of the first-order approximated spectrum (square) parallel to the NIR axis instead of toward the true reflectance spectrum (circle). Fortunately, this shift direction and distance compensated for the truncation error. As a result, the isoline of the asymmetric approximation (red line) ran through the subspace between the true (circle) and second-order approximated spectra (cross). In summary, the overcorrection of the NIR band and the truncation error in the NIR band nicely canceled each other out, thereby shifting the isoline (red) upward into the subspace to decrease the distance between the true spectrum and the isoline.

Fig. 7

Illustration of the error reduction mechanism in the asymmetric order approximated isoline.

5. Discussion and Conclusions

The truncation of higher-order interaction terms presents a major limitation to isoline models based on a first-order approximation, although this truncation can simplify the final expression. The simplicity of the analytical form is advantageous for analytical and numerical investigations of parameter retrieval algorithms, such as the LAI, FVC, as well as of proximity measures, such as the spectral VI. The accuracy of the vegetation isoline models may be increased by including the second-order interaction terms in both the red and NIR reflectances; however, the derived expression was rather complex, thereby reducing the utility of the model as a tool for analytical and numerical studies. This study took a unique approach: instead of including the second-order terms in both the red and NIR bands, this study included the term only in the NIR reflectance.

The final form of the derived asymmetric isoline is rather simple; thus, it can be easily altered, similar to the previously derived first-order isoline. This model, however, dramatically reduces the errors obtained in the first-order isoline, and the accuracy of the asymmetric case is even better than that obtained from the reflectance spectrum using the second-order term. With both simplicity and accuracy, the derived expression can contribute to a wide range of applications, from designing optimal spectral VI sets to developing inversion algorithms in which the derived expression may be used as a constraint in the optimization algorithm.

This study focused on only the relationship between the red and NIR reflectances. Although overcorrections in the NIR band nicely compensated for the truncation errors inherent in the isoline equation, this compensation mechanism may not apply to other combinations of wavelengths. The applicability of this model to other combinations, e.g., the NIR and shortwave bands, will require more thorough investigations. The findings of this study are currently limited to the combination of red and NIR wavelengths.

Overcorrection by intentionally truncating more terms in the red band than in the NIR band compensated for the truncation errors of the higher-order term in a derivation of the vegetation isoline. These findings suggest that optimal control over the overcorrection level could further reduce the errors in the asymmetric isoline. This possibility is worth exploring in future studies.

Acknowledgments

This work was supported by JSPS KAKENHI Grant No. 15H02865.

References

1.

J. Qin et al., “Development of the adjoint model of a canopy radiative transfer model for sensitivity study and inversion of leaf area index,” IEEE Trans. Geosci. Remote Sens. 46(7), 2028–2037 (2008).http://dx.doi.org/10.1109/TGRS.2008.916637Google Scholar

2.

A. R. Huete, “A soil-adjusted vegetation index (SAVI),” Remote Sens. Environ. 25(3), 295–309 (1988).http://dx.doi.org/10.1016/0034-4257(88)90106-XGoogle Scholar

3.

F. Baret, G. Guyot and D. Major, “TSAVI: a vegetation index which minimizes soil brightness effects on LAI and APAR estimation,” in 12th Canadian Symposium on Remote Sensing, Int. Geoscience and Remote Sensing Symposium, Vol. 3, pp. 1355–1358, IEEE (1989).Google Scholar

4.

D. Major, F. Baret and G. Guyot, “A ratio vegetation index adjusted for soil brightness,” Int. J. Remote Sens. 11(5), 727–740 (1990).http://dx.doi.org/10.1080/01431169008955053Google Scholar

5.

J. Qi et al., “A modified soil adjusted vegetation index,” Remote Sens. Environ. 48(2), 119–126 (1994).http://dx.doi.org/10.1016/0034-4257(94)90134-1Google Scholar

6.

M. Gilabert et al., “A generalized soil-adjusted vegetation index,” Remote Sens. Environ. 82(2), 303–310 (2002).http://dx.doi.org/10.1016/S0034-4257(02)00048-2Google Scholar

7.

Z. Jiang et al., “Development of a two-band enhanced vegetation index without a blue band,” Remote Sens. Environ. 112(10), 3833–3845 (2008).http://dx.doi.org/10.1016/j.rse.2008.06.006Google Scholar

8.

R. Myneni et al., “Global products of vegetation leaf area and fraction absorbed PAR from year one of MODIS data,” Remote Sens. Environ. 83(1), 214–231 (2002).http://dx.doi.org/10.1016/S0034-4257(02)00074-3Google Scholar

9.

R. Darvishzadeh et al., “Inversion of a radiative transfer model for estimating vegetation LAI and chlorophyll in a heterogeneous grassland,” Remote Sens. Environ. 112(5), 2592–2604 (2008).http://dx.doi.org/10.1016/j.rse.2007.12.003Google Scholar

10.

J. A. Smith, “LAI inversion using a back-propagation neural network trained with a multiple scattering model,” IEEE Trans. Geosci. Remote Sens. 31(5), 1102–1106 (1993).http://dx.doi.org/10.1109/36.263783Google Scholar

11.

C. Bacour et al., “Neural network estimation of LAI, fAPAR, fCover and LAI $\times$ C ab, from top of canopy MERIS reflectance data: principles and validation,” Remote Sens. Environ. 105(4), 313–325 (2006).http://dx.doi.org/10.1016/j.rse.2006.07.014Google Scholar

12.

F. Baret et al., “LAI, fAPAR and fCover CYCLOPES global products derived from VEGETATION: Part 1: principles of the algorithm,” Remote Sens. Environ. 110(3), 275–286 (2007).http://dx.doi.org/10.1016/j.rse.2007.02.018Google Scholar

13.

S. Jacquemoud et al., “Extraction of vegetation biophysical parameters by inversion of the PROSPECT+ SAIL models on sugar beet canopy reflectance data. Application to TM and AVIRIS sensors,” Remote Sens. Env. 52(3), 163–172 (1995).http://dx.doi.org/10.1016/0034-4257(95)00018-VGoogle Scholar

14.

R. B. Myneni et al., “Estimation of global leaf area index and absorbed PAR using radiative transfer models,” IEEE Trans. Geosci. Remote Sens. 35(6), 1380–1393 (1997).http://dx.doi.org/10.1109/36.649788Google Scholar

15.

Y. Knyazikhin et al., “Synergistic algorithm for estimating vegetation canopy leaf area index and fraction of absorbed photosynthetically active radiation from MODIS and MISR data,” J. Geophys. Res. 103(D24), 32257–32275 (1998).http://dx.doi.org/10.1029/98JD02462Google Scholar

16.

R. D. Jackson and A. R. Huete, “Interpreting vegetation indices,” Prev. Vet. Med. 11(3), 185–200 (1991).http://dx.doi.org/10.1016/S0167-5877(05)80004-2Google Scholar

17.

F. Baret and G. Guyot, “Potentials and limits of vegetation indices for LAI and APAR assessment,” Remote Sens. Environ. 35(2), 161–173 (1991).http://dx.doi.org/10.1016/0034-4257(91)90009-UGoogle Scholar

18.

M. M. Verstraete and B. Pinty, “Designing optimal spectral indexes for remote sensing applications,” IEEE Trans. Geosci. Remote Sens. 34(5), 1254–1265 (1996).http://dx.doi.org/10.1109/36.536541Google Scholar

19.

M. M. Verstraete, B. Pinty and R. B. Myneni, “Potential and limitations of information extraction on the terrestrial biosphere from satellite remote sensing,” Remote Sens. Environ. 58(2), 201–214 (1996).http://dx.doi.org/10.1016/S0034-4257(96)00069-7Google Scholar

20.

J. Hanes, Biophysical Applications of Satellite Remote Sensing, Springer Science & Business Media, Berlin Heidelberg (2013).Google Scholar

21.

R. Myneni and D. Williams, “On the relationship between FAPAR and NDVI,” Remote Sens. Environ. 49(3), 200–211 (1994).http://dx.doi.org/10.1016/0034-4257(94)90016-7Google Scholar

22.

G. Asrar et al., “Theory and applications of optical remote sensing,” Theory Appl. Opt. Remote Sens. (1989).http://dx.doi.org/10.2307/623085Google Scholar

23.

N. Gobron, B. Pinty and M. M. Verstraete, “Theoretical limits to the estimation of the leaf area index on the basis of visible and near-infrared remote sensing data,” IEEE Trans. Geosci. Remote Sens. 35(6), 1438–1445 (1997).http://dx.doi.org/10.1109/36.649798Google Scholar

24.

H. Yoshioka et al., “Analysis of vegetation isolines in red-NIR reflectance space,” Remote Sens. Environ. 74(2), 313–326 (2000).http://dx.doi.org/10.1016/S0034-4257(00)00130-9Google Scholar

25.

N. H. Broge and E. Leblanc, “Comparing prediction power and stability of broadband and hyperspectral vegetation indices for estimation of green leaf area index and canopy chlorophyll density,” Remote Sens. Environ. 76(2), 156–172 (2001).http://dx.doi.org/10.1016/S0034-4257(00)00197-8Google Scholar

26.

Z. Jiang et al., “Analysis of NDVI and scaled difference vegetation index retrievals of vegetation fraction,” Remote Sens. Environ. 101(3), 366–378 (2006).http://dx.doi.org/10.1016/j.rse.2006.01.003Google Scholar

27.

Z. Jiang et al., “Interpretation of the modified soil-adjusted vegetation index isolines in red-NIR reflectance space,” J. Appl. Remote Sens. 1(1), 013503 (2007).http://dx.doi.org/10.1117/1.2709702Google Scholar

28.

H. Yoshioka et al., “Derivation of soil line influence on two-band vegetation indices and vegetation isolines,” Remote Sens. 1(4), 842–857 (2009).http://dx.doi.org/10.3390/rs1040842Google Scholar

29.

H. Yoshioka, H. Yamamoto and T. Miura, “Use of an isoline-based inversion technique to retrieve a leaf area index for inter-sensor calibration of spectral vegetation index,” in IEEE Int. Geoscience and Remote Sensing Symposium, Vol. 3, pp. 1639–1641, IEEE (2002).Google Scholar

30.

A. Kallel et al., “Determination of vegetation cover fraction by inversion of a four-parameter model based on isoline parametrization,” Remote Sens. Environ. 111(4), 553–566 (2007).http://dx.doi.org/10.1016/j.rse.2007.04.006Google Scholar

31.

H. Yoshioka, T. Miura and A. R. Huete, “An isoline-based translation technique of spectral vegetation index using EO-1 Hyperion data,” IEEE Trans. Geosci. Remote Sens. 41(6), 1363–1372 (2003).http://dx.doi.org/10.1109/TGRS.2003.813212Google Scholar

32.

H. Yoshioka, T. Miura and K. Obata, “Derivation of relationships between spectral vegetation indices from multiple sensors based on vegetation isolines,” Remote Sensing 4(3), 583–597 (2012).RSEND3http://dx.doi.org/10.3390/rs4030583Google Scholar

33.

K. Obata et al., “Derivation of a MODIS-compatible enhanced vegetation index from visible infrared imaging radiometer suite spectral reflectances using vegetation isoline equations,” J. Appl. Remote Sens. 7(1), 073467 (2013).http://dx.doi.org/10.1117/1.JRS.7.073467Google Scholar

34.

H. Yoshioka, A. R. Huete and T. Miura, “Derivation of vegetation isoline equations in red-NIR reflectance space,” IEEE Trans. Geosci. Remote Sens. 38(2), 838–848 (2000).http://dx.doi.org/10.1109/36.842012Google Scholar

35.

H. Yoshioka, “Vegetation isoline equations for an atmosphere-canopy-soil system,” IEEE Trans. Geosci. Remote Sens. 42(1), 166–175 (2004).http://dx.doi.org/10.1109/TGRS.2003.817793Google Scholar

36.

M. Miura and H. Yoshioka, “Vegetation isoline equations with higher order interaction terms to improve accuracy in retrieval of vegetation biophysical parameters,” in IEEE Int. Geoscience and Remote Sensing Symposium, Vol. 3, p. 1056, IEEE (2008).Google Scholar

37.

M. Miura, K. Obata and H. Yoshioka, “Vegetation isoline equations for analysis of hyper-spectral data with higher order interaction terms,” in IEEE Int. Geoscience and Remote Sensing Symposium, pp. 1031–1034, IEEE (2010).Google Scholar

38.

M. Miura, K. Obata and H. Yoshioka, “Vegetation isoline equations for atmosphere-canopy-soil system of layer with second order interaction term,” Proc. SPIE 7809, 780906 (2010).http://dx.doi.org/10.1117/12.860432Google Scholar

39.

S. Jacquemoud et al., “PROSPECT+ SAIL models: a review of use for vegetation characterization,” Remote Sens. Environ. 113, S56–S66 (2009).http://dx.doi.org/10.1016/j.rse.2008.01.026Google Scholar

40.

K. Cooper, J. Smith and D. Pitts, “Reflectance of a vegetation canopy using the adding method,” Appl. Opt. 21(22), 4112–4118 (1982).http://dx.doi.org/10.1364/AO.21.004112Google Scholar

41.

F. Baret, S. Jacquemoud and J. Hanocq, “The soil line concept in remote sensing,” Remote Sens. Rev. 7(1), 65–82 (1993).http://dx.doi.org/10.1080/02757259309532166Google Scholar

42.

S. Jacquemoud and F. Baret, “PROSPECT: a model of leaf optical properties spectra,” Remote Sens. Environ. 34(2), 75–91 (1990).http://dx.doi.org/10.1016/0034-4257(90)90100-ZGoogle Scholar

43.

W. Verhoef, “Earth observation modeling based on layer scattering matrices,” Remote Sens. Environ. 17(2), 165–178 (1985).http://dx.doi.org/10.1016/0034-4257(85)90072-0Google Scholar

Biography

Munenori Miura is a PhD student at Aichi Prefectural University, Aichi, Japan. He received his BS and MS degrees in information science and technology from Aichi Prefectural University in 2008 and 2010, respectively.

Kenta Obata is a researcher in the Institute of Geology and Geoinformation at the National Institute of Advanced Industrial Science and Technology, Japan. He received his MS and PhD degrees in information science and technology from Aichi Prefectural University, Aichi, Japan, in 2010 and 2012, respectively. His research interests include intersensor calibration of biophysical retrievals, spectral mixture analysis, and radiometric calibration.

Hiroki Yoshioka is a professor in the Department of Information Science and Technology at Aichi Prefectural University, Aichi, Japan. He received his MS and PhD degrees in nuclear engineering from the University of Arizona, Tucson, in 1993 and 1999, respectively. His research interests include cross-calibrations of satellite data products for continuity and compatibility of similar data sets, and development of canopy radiative transfer models and their inversion techniques.

Citation Download Citation

Munenori Miura, Kenta Obata, Hiroki Yoshioka, "Vegetation isoline equations with first- and second-order interaction terms for modeling a canopy-soil system of layers in the red and near-infrared reflectance space," Journal of Applied Remote Sensing 9(1), 095987 (18 November 2015). http://dx.doi.org/10.1117/1.JRS.9.095987

Submission: Received ; Accepted

Access the abstract

JOURNAL ARTICLE
14 PAGES

DOWNLOAD PDF + Save for Later

SHARE

GET CITATION

< Previous Article

|

Next Article >

KEYWORDS

Reflectivity

Vegetation

Near infrared

Soil science

Spherical lenses

Systems modeling

Error analysis

Show All Keywords

Subscribe to Digital Library

Receive Journal Email Alerts

Create, Organize, Share

1.

Introduction

2.

Analytical Model and First-Order Vegetation Isoline

2.1.

Analytical Canopy Reflectance Model

(1)

(2)

2.2.

Vegetation Isoline Equation with a First-Order Interaction Term in the Red–Near-Infrared Reflectance Subspace

(3)

(4)

(5)

(6)

(7)

(8)

(9)

(10)

(11)

(12)

(13)

Fig. 1

Fig. 2

3.

Second-Order Approximation and Asymmetric Order Isoline

3.1.

Parametric Form of the Second-Order Approximated Reflectance Spectra

(14)

(15)

(16)

(17)

(18)

3.2.

Isoline Equation Obtained by Asymmetrically Truncating the Higher-Order Interactions

(19)

(20)

(21)

(22)

(23)

(24)

(25)

(26)

4.

Results of the Numerical Simulations

4.1.

Parameter Settings Used in the Numerical Experiments

4.2.

Numerical Procedure Used for the Isoline Parameter Retrieval and the Error Estimation

4.3.

Comparison of the Accuracy Across the Three Approximations

Fig. 3

Table 1

Fig. 4

Fig. 5

4.4.

Error Reduction Mechanisms Using the Second-Order Isoline

Fig. 6

Fig. 7

5.

Discussion and Conclusions

Acknowledgments

References

Biography

Create, Organize, Share

Show All Keywords

KEYWORDS/PHRASES

SEARCH IN:

PUBLICATION YEARS