Phytoplankton Photophysiology Utils

This is a tool to read and process active chlorophyll fluorescence data from raw format and apply the biophysical model of Kolber et al. (1998). For more information see the documentation and demo file, below is a short example of how to use the data to read in and process variables.

EXAMPLE USAGE

This package is meant to be used in an interactive environment - ideally JuPyTer Notebook

import phyto_photo_utils as ppu

fname = '/path_to_data/data'
output = '/output_path'

# Load all variables needed for fitting saturation and relaxation models
df = ppu.load_FASTTrackaI_files(fname, append=False, save_files=True, res_path=output, seq_len=120, irrad=545.62e10)

# Perform a calculate ro saturation model fit on the data
sat = ppu.fit_saturation(pfd, flevel, seq, datetime, blank=0, sat_len=100, skip=0, model_type='calc_ro', bounds=True, ro_lims=[0.01,1.0], sig_lims =[100,2200])

# Perform a single decay relaxation model fit on the data
rel = ppu.fit_relaxation(flevel, seq_time, seq, datetime, blank=0, sat_len=100, rel_len=40, model_type='single', bounds=True, tau1_lims=[100, 50000])

# Perform time averaging (5 minute averages) on raw transients, including the removal of outliers (mean + stdev * 3)
dfm = ppu.remove_outlier_from_time_average(df, time=5, multiplier=3)

# Correct for FIRe instrument detector bias
dfb = ppu.correct_fire_instrument_bias(df, sat=False, pos=1, sat_len=100)

# See the demo file for more info

ABOUT

This work was funded by the CSIR. This research was partially supported by the Australian Government through the Australian Research Council's Discovery Projects funding scheme (DP160103387).

• Version: 1.6.4
• Author: Thomas Ryan-Keogh, Charlotte Robinson
• Email: tjryankeogh(at)gmail.com
• Date: 2018-12-06
• Institution: Council for Scientific and Industrial Research, Curtin University
• Research group: Southern Ocean Carbon - Climate Observatory (SOCCO), Remote Sensing and Satellite Research Group

Please use the guidlines given on https://integrity.mit.edu/handbook/writing-code to cite this code.

Example citation: Source: phyto_photo_utils [https://gitlab.com/tjryankeogh/phytophotoutils] retrieved on 30 May 2019.

PACKAGE STRUCTURE

NOTE: This package structure is defined by the __init__.py file

• load
  • load_FIRe_files
  • load_FASTTrackaI_files
  • load_FastOcean_files
  • load_LIFT_FRR_files
• saturation
  • fit_saturation
• relaxation
  • fit_relaxation
• tools
  • remove_outlier_from_time_average
  • correct_fire_instrument_bias
  • calculate_blank_FastOcean _ calculate_blank_FIRe
• spectral_correction
  • calculate_chl_specific_absorption
  • calculate_instrument_led_correction
• etr
  • calculate_amplitude_etr
• plot
  • plot_saturation_data
  • plot_relaxation_data
  • plot_fluorescence_light_curve
• equations
  • fit_kolber_nop
  • calculatate_residual_saturation_nop
  • fit_kolber_p
  • calculate_residual_saturation_p
  • fit_single_relaxation
  • calculate_residual_single_relaxation
  • fit_triple_relaxation
  • calculate_residual_triple_relaxation
  • calculate_alpha_model
  • calculate_residual_etr
  • calculate_modified_alpha_model
  • calculate_residual_phi
  • calculate_beta_model
  • calculate_residual_beta
  • calculate_modified_beta_model
  • calculate_residual_mbeta
  • calculate_bias
  • calculate_rmse
  • calculate_nrmse
  • calculate_fit_errors
• fitting
  • fit_kpf_saturation
  • fit_kpf_relaxation

ACKNOWLEDGEMENTS

• Kevin Oxborough (Chelsea Technology Groups) for linear methods to estimate Fo and Fm

TO DO

• Add a_PSII calculation. PFD (micromol quanta m-2 s-1) * sigmaPSII (m-2 PSII-1).
• Add in function to read FastTracka I binary files
• Add deconvolution method for triple decay relaxation
• Add additional load method for CTG LabSTAF
• Add function to fit saturation and relaxation together
• Check fit definitions match SOP terminology.