File Exchange

image thumbnail

Aggregate timeseries data to solstices and equinoxes.

version 1.0 (2.47 KB) by

Aggregates serial datetime data to summer/winter solstices, and spring/fall equinoxes.

1 Download

Updated

View License

accumSeasonTS- aggregate timeseries to solstices and equinoxes

Using G. Cardillo function to compute solstices and equinoxes, this function will aggregate time series data (e.g. Precipitation) to seasonally summarized results. To accomplish this three core functions are used: eqnsol, unique, and accumarray. The rest of the function is manipulating data to pass around in between the functions.

This function will output cumulative seasonal totals of data. It's
application is limited default parameterization of ACCUMARRAY matlab function (i.e. sum totals). It's certainly possible to modify this function to allow the user to apply different statistics in the
ACCUMARRAY function, but not done here.

Additionally, this function uses another function (EQNSOL) to calculate the summer/winter solstices and spring/fall equinoxes per year, then sum up the data within each season.

Function does not discriminate on partial seasons. Partial seasons that are either because of beginning or ending dates, or because of missing data are not censored. Partial seasons will be output along with all complete seasons. It's up to the user to either pre-process the data or post process the results.

Expected inputs is a timeseries of data in two columns: a decimal date, and the values to be accumulated.

datain = n x 2
  datain(:,1) = decimal date/time
  datain(:,2) = values to sum by season
fn = supplied string specifying a valid summary statistic function.
Example: @sum, @max, @min, etc.

dataout = n x 2
  dataout(:,1) = decimal date/time of the start of the season
  dataout(:,2) = values summed for the season

Acknowledgements:

Cardillo G. (2007) Equinoxes and Solstices: compute the date and time of equinoxes and solstices.
http://www.mathworks.com/matlabcentral/fileexchange/17977
  
Written by:
  Jeff Burkey
  King County- DNRP
  email: jeff.burkey@kingcounty.gov
  May 30, 2009

Syntax:
  [dataout] = accumSeasonTS(datain,fn)

Comments and Ratings (0)

MATLAB Release
MATLAB 7.8 (R2009a)
Acknowledgements

Inspired by: Equinoxes and Solstices

Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.

» Watch video