Py6S is a Python interface to the 6S Radiative Transfer Model. It allows you to run many 6S simulations using a simple Python syntax, rather than dealing with the rather cryptic 6S input and output files. As well as generally making it easier to use 6S, Py6S adds some new features:
- The ability to run many simulations easily and quickly, with no manual editing of input files
- The ability to run for many wavelengths and/or angles and easily plot the results
- The ability to import real-world data to parameterise 6S, from radiosonde measurements and AERONET sun photometer measurements
Py6S has been designed to be easy to use, and to work on the ‘principle of least surprise’. Far more details are available in the rest of this documentation, but a quick code example should give you an idea of what Py6S can do:
# Import the Py6S module from Py6S import * # Create a SixS object s = SixS() # Set the wavelength to 0.675um s.wavelength = Wavelength(0.675) # Set the aerosol profile to Maritime s.aero_profile = AeroProfile.PredefinedType(AeroProfile.Maritime) # Run the model s.run() # Print some outputs print(s.outputs.pixel_reflectance, s.outputs.pixel_radiance, s.outputs.direct_solar_irradiance) # Run the model across the VNIR wavelengths, and plot the result wavelengths, results = SixSHelpers.Wavelengths.run_vnir(s, output_name='pixel_radiance') SixSHelpers.Wavelengths.plot_wavelengths(wavelengths, results, "Pixel radiance ($W/m^2$)")
This will produce the results shown below:
0.283 112.095 667.589
Py6S was described in a journal article in Computers and Geosciences which must be cited if Py6S is used for producing outputs for a scientific report/publication - see the Publications page for more information.
I’m very interested in receiving feedback, bug reports and feature suggestions - please see the Support page to find out how to get in touch with me.