In jbferet/prospect: PROSPECT leaf radiative transfer model and inversion routines

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>",
  eval=FALSE
)

Brief introduction & useful links to go further

This tutorial summarizes the different functionalities of the R package prospect based on the eponym leaf physical model.

This tutorial does not aim at detailling the principles of PROSPECT or comparing it to other models, as:

• The comprehensive description of the model and the physical principes it relies on can be found in the original paper from Jacquemoud & Baret (1990)
• this website and this one provide a lot of information on this topic
• this website gathers references and important links for anyone interested in leaf optical properties and physical modeling

PROSPECT aims at simulating leaf optical properties in the optical domain from 400 nm to 2500 nm based on their biophysical properties, including a limited number of biochemical constituents and a unique structure parameter, N. PROSPECT is based on a simple representation derived from the extended plate model (Allen et al., 1970).

Leaf optical properties correspond to conical-hemispherical reflectance and transmittance, and are typically measured with an integrating sphere. These optical properties are often described as directional-hemispherical reflectance and transmittance, although such directional quantities are conceptual quantities (see Schaepman-Strub et al. (2006) for a comprehensive definition of optical measurements).

Leaf reflectance measured with a leaf clip / contact probe (such as the device provided with the ASD FieldSpec spectroradiometer) does not correspond to conical-hemispherical reflectance. Therefore comparing such reflectance with PROSPECT simulations (in forward and inverse mode) may lead to biased or uncertain results. See Li et al. (2018) for alternative approaches taking advantage ofphysical modeling for the estimation of leaf chemistry from such reflectance measurements.

 

 

Fig. 1. Representation of a leaf according to PROSPECT and corresponding leaf optical properties obtained as output

 

Note that the alpha parameter is also available as input, and corresponds to the maximum incidence angle relative to the normal defining the solid angle of incident light at the surface of leaf (aims at including surface roughness, default value = 40 degrees).

Many simplications can be identified, and alternative versions of the model have been developed through the years in order to increase the realism of the model. Among the many simplifications:

• The refractive index identifcal for all leaves. It should theoretically change with leaf properties

• The leaf anatomy is not differentiated between an adaxial and an abaxial face

• the leaf surface, which is defined my many properties in real life (presence of waxes, hairs, ...) is identical among leaves in PROSPECT

• ...

The current version of the model implemented in prospect is PROSPECT-PRO, which includes the following biochemical constituents (defined as C_i in Fig. 1):

• Chlorophyll a + b CHL
• Carotenoids CAR
• Anthocyanins ANT
• Prown pigments BROWN
• Equivalent water thickness EWT
• proteins PROT
• carbon-based consituents CBC (constituents of dry matter other than proteins)

If you want to use PROSPECT-D instead of PROSPECT-PRO, please define a value for LMA (Leaf mass per area) and set PROT and CBC to 0, or leave no value (default value = 0)

The specific absoprtion coefficients corresponding to these constituents are recorded in the variable SpecPROSPECT available when loading the package prospect.

prospect: forward and inverse mode

In the forward mode, PROSPECT simulates leaf optical propertes based on a set of input parameters corresponding to the biochemical constituents and the N structure parameter.

The performances of PROSPECT for the simulation of the leaf optical properties are based on the proper calibration of the optical constants, more particularly the specific absorption coefficients corresponding to each biochemical constituent.

In the inverse mode, an algorithm is used to derive the input parameters from the leaf optical properties.

See the following pages of the tutorial for an illustration of how to use prospect in forward and inverse mode.

 

 

Fig. 2. Forward and inverse modes of PROSPECT

 

prospect: Miscellaneous

Fluorescence & extended infrared domain

The package currently only includes a version of the model covering the domain from 400 nm to 2500 nm. However, several valuable alternative versions exist, but are not available in the current version of the package:

• PROSPECT-VISIR, modeling directional–hemispherical reflectance and transmittance of fresh and dry leaves from 0.4 μm to 5.7 μm.
• FLUSPECT, radiative transfer model for leaf chlorophyll fluorescence.

These models may be added in a future version of the package.

Alternative distributions and coupling with SAIL

A version of the R package prospect coupled with the SAIL canopy model is available in the R package prosail.

Alternative implementations of PROSPECT and PROSAIL can be found in various languages. A non-exhaustive list of distributions of these models is available at this webpage. It includes Matlab, Fortran and R implementations of versions 4, 5, D and PRO of PROSPECT.

We strongly recommend using version D or PRO instead of versions 4 and 5.

Note that additional versions of PROSPECT and PROSAIL model are also available in packages written in python, Julia and R.

jbferet/prospect documentation built on Feb. 10, 2025, 9:35 a.m.