File Exchange

image thumbnail


version (5.54 MB) by Chad Greene
Calculate monthly means, annual medians, and other statistics for 1D or 3D datasets.


Updated 30 Dec 2014

View Version History

View License

This function downsamples 1D or 3D data to monthly, yearly, hourly, minutely, or secondly data. This function was originally designed to create monthly mean time series from daily geospatial climate data.
Z_downsamp = downsample_ts(Z,t)
Z_downsamp = downsample_ts(...,'DownsamplingPeriod')
Z_downsamp = downsample_ts(...,'function')
[Z_downsamp,t_downsamp] = downsample_ts(...)

Z_downsamp = downsample_ts(Z,t) downsamples Z, which must be provided with a corresponding time vector t. Z can be 1D if its length matches the length of t. If Z is three-dimensional, the length of its third dimension must match the length of t. For geospatial climate data arrays, dimensions of Z might correspond to lat x lon x time or lon x lat x time.

Z_downsamp = downsample_ts(...,'DownsamplingPeriod') specifies a downsampling period as

'month' (default)

Z_downsamp = downsample_ts(...,'function') specifies a function to perform on the data. By default, monthly averages are taken, but you may wish to return the monthly median or monthly standard deviation or any of the functions listed below.

A note on functions which ignore NaNs: To get the monthly means of data while ignoring NaN values, you can use the 'nanmean' option. The nanmean function is part of the Statistics Toolbox, but may also be found as part of the NaN Suite on File Exchange. However, the File Exchange versions mix up the order of dimensions and flags for nanstd, nanvar, nanmin, and nanmax, so you will need the Statistics Toolbox for those particular functions. In all, the following functions are available:

'mean' (default)
'nanmean' ignores NaN values in Z. Requires Statistics toolbox or NaN Suite.
'nanmedian' ignores NaN values in Z. Requires Statistics toolbox or NaN Suite.
'nanmin' ignores NaN values in Z. Requires Statistics toolbox.
'nanmax' ignores NaN values in Z. Requires Statistics toolbox.
'std' standard deviation.
'nanstd' ignores NaN values in Z. Requires Statistics toolbox.
'var' variance.
'nanvar' ignores NaN values in Z. Requires Statistics toolbox.

[Z_downsamp,t_downsamp] = downsample_ts(...) also returns a time array corresponding to Z_downsamp. If Z is 3D or, t_downsamp corresponds to the third dimension of Z_downsamp. Each value in t_downsamp represents the mean time of all data contributing to that slice of Z_downsamp.

Cite As

Chad Greene (2020). downsample_ts (, MATLAB Central File Exchange. Retrieved .

Comments and Ratings (5)

Veronica Bozzini

thanks, exactly what I was looking for!

Peter Mills

My data has lots of missing timestamps. I get the message
"Warning: Some time steps in input t are more than a 2 hours apart. This may cause strange results when taking hourly
> In downsample_ts (line 207) "
Does that mean I can't trust the result?
How can I convert 1 min wind speeds to hourly wind speeds when the data has lots of missing values. So far I have rounded the timestamps to the hour and taken an average of the values with the same hourly timestamp. Like this:

% Convert the time stamps to 60 minutes invervals rounding down to the
% last 60 min

%pre define

for i=1:length(dateTime60min)
if sum(isnan(windSpeed(dateTime60==dateTime60min(i))))==0

mohammed ali

very nice script,thanks for your excellent work

pedro dollar

rajitha reddy putta

MATLAB Release Compatibility
Created with R2014b
Compatible with any release
Platform Compatibility
Windows macOS Linux

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!