Skip to content

Python code for running ProspectD and 4SAIL Rasiative Transfer Models

License

Notifications You must be signed in to change notification settings

hectornieto/pypro4sail

Repository files navigation

pyPro4Sail

Prospect5 and 4SAIL Radiative Transfer Models for simulating the transmission of radiation in leaves and canopies.

Synopsis

This project contains Python code for Prospect and 4SAIL Radiative Transfer Models (RTM) for simulating the transmission of optical and thermal electromagnetic radiation through leaves and vegetated canopies.

The project consists of:

  1. lower-level modules with the basic functions needed in Prospec5 and 4SAIL RTMs.

  2. higher-level scripts for easily running ProSAIL in both forward and inverse mode.

Installation

Download the project to your local system, enter the download directory and then type.

python setup.py install

if you want to install pyTSEB and its low-level modules in your Python distribution.

The following Python library ir required for running Prospect and 4SAIL:

  • Numpy

In addition, the inversion of both RTMS requires.

  • Scipy

  • cma [Optional]

Code Example

High-level example

You can automatically run the coupled leaf+canopyt Prospect5+4SAIL RTM with pyPro4Sail.py module.

[N, chloro, caroten, brown, EWT, LMA, LAI, hot_spot, solar_zenith, solar_azimuth, view_zenith, view_azimuth, LIDF]=[1.5,40,8,0.0,0.01,0.009,3,0.01,30,180,10,180,(-0.35,-0.15)]
import pyPro4SAIL
wl,rho=pyPro4SAIL.run(N, chloro, caroten, brown, EWT, LMA, LAI, hot_spot, solar_zenith, solar_azimuth, view_zenith, view_azimuth, LIDF, skyl=0.2, soilType=pyPro4SAIL.DEFAULT_SOIL)

Also it is possible to simulate the surface land-leaving thermal radiance with the function run_TIR.

Low-level example

Prospect5 RTM

You can run Prospect by importing the module Prospect5.py and then either calling the function Prospect5 for simulating the full optical spectrum (400-2500nm), or the function Prospect5_wl for simulating the leaf reflectance and transmittance for a given wavelength.

# Running Prospect5
import Prospect5
# Simulate leaf full optical spectrum (400-2500nm) 
wl, rho_leaf, tau_leaf = Prospect5.Prospect5(N, chloro, caroten, brown, EWT, LMA)

You can type help(Prospect5.Prospect5) to understand better the inputs needed and the outputs returned

4SAIL RTM

You can run 4SAIL by importing the module FourSAIL.py and then either calling the function FourSAIL for simulating the reflectance and transmittance factor of a given canopy given a list of leaf reflectances and trasmittances, or you can call the function FourSAIL_wl for simulating the leaf reflectance and transmittance factor of a given canopy at for a single wavelenght.

# Running the coupled Prospect and 4SAIL
import Prospect5, FourSAIL
# Simulate leaf full optical spectrum (400-2500nm) 
wl, rho_leaf, tau_leaf = Prospect5.Prospect5(N, chloro, caroten, brown, EWT, LMA)
# Estimate the Leaf Inclination Distribution Function of a canopy
LIDF = FourSAIL.CalcLIDF_Campbell(alpha)
# Simulate leaf reflectance and transmittance factors of a canopy 
tss,too,tsstoo,rdd,tdd,rsd,tsd,rdo,tdo,rso,rsos,rsod,rddt,rsdt,rdot,rsodt,rsost,rsot,gammasdf,gammasdb,gammasowl = FourSAIL.FourSAIL(LAI,hot_spot,LIDF,solar_zenith,view_zenith,solar_azimuth-view_azimuth,rho_leaf,tau_leaf,rho_soil)
# Simulate the canopy reflectance factor for a given difuse/total radiation condition (skyl)
rho_canopy = rdot*skyl+rsot*(1-skyl)

You can type help(FourSAIL.FourSAIL) to understand better the inputs needed and the outputs returned

Basic Contents

High-level modules

  • .src/pyPro4SAIL.py

Runs the coupled Prospect5+4SAIL to estimate the canopy directional reflectance factor and 4SAIL to estimate the land-leaving broadband thermal radiance.

Low-level modules

The low-level modules in this project are aimed at providing customisation and more flexibility in running TSEB. The following modules are included.

  • .src/Prospect5.py

core functions for running Prospect5 Leaf Radiative Transfer Model.

  • .src/Prospect5Jacobian.py

core functions for computing the Jacobian of Prospect5 Leaf Radiative Transfer Model.

  • .src/FourSAIL.py

core functions for running 4SAIL Canopy Radiative Transfer Model.

  • .src/FourSAILJacobian.py

core functions for computing the Jacobian of 4SAIL Canopy Radiative Transfer Model.

  • .src/CostFunctionsPROSPECT4SAIL.py

merit functions used to invert Prospect and/or 4SAIL from a given spectrum

  • .src/cma.py

Covariance Matrix Adaptation Evolution Strategy optimization method for inverting Prospect5 and/or 4SAIL.

API Reference

http://pyPro4Sail.readthedocs.org/en/latest/index.html

Main Scientific References

  • S. Jacquemoud, F. Baret, PROSPECT: A model of leaf optical properties spectra, Remote Sensing of Environment, Volume 34, Issue 2, November 1990, Pages 75-91, ISSN 0034-4257, http://dx.doi.org/10.1016/0034-4257(90)90100-Z.
  • Jean-Baptiste Feret, Christophe François, Gregory P. Asner, Anatoly A. Gitelson, Roberta E. Martin, Luc P.R. Bidel, Susan L. Ustin, Guerric le Maire, Stéphane Jacquemoud, PROSPECT-4 and 5: Advances in the leaf optical properties model separating photosynthetic pigments, Remote Sensing of Environment, Volume 112, Issue 6, 16 June 2008, Pages 3030-3043, ISSN 0034-4257, http://dx.doi.org/10.1016/j.rse.2008.02.012.
  • W. Verhoef, Light scattering by leaf layers with application to canopy reflectance modeling: The SAIL model, Remote Sensing of Environment, Volume 16, Issue 2, October 1984, Pages 125-141, ISSN 0034-4257, http://dx.doi.org/10.1016/0034-4257(84)90057-9.
  • W. Verhoef, L. Jia, Q. Xiao and Z. Su, Unified Optical-Thermal Four-Stream Radiative Transfer Theory for Homogeneous Vegetation Canopies, IEEE Transactions on Geoscience and Remote Sensing, vol. 45, no. 6, pp. 1808-1822, June 2007. http://dx.doi.org/10.1109/TGRS.2007.895844.

Tests

to be included

Contributors

License

pyPro4Sail: a Python Two Source Energy Balance Model

Copyright 2016 Hector Nieto and contributors.

This program is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version.

This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.

You should have received a copy of the GNU General Public License along with this program. If not, see http://www.gnu.org/licenses/.

About

Python code for running ProspectD and 4SAIL Rasiative Transfer Models

Resources

License

Stars

Watchers

Forks

Packages

No packages published

Languages