R Package for harmonic modelling of time-series data

R Package for harmonic modelling of time-series data

April 23, 2020

Sentinel-2 NDVI time-series over the Steigerwald. Left: Original satellite scenes after cloud, cloud shadow and snow maksing. Right: Interpolated time-series using a harmonic modelling.

In order to fully exploit the monitoring potential of the satellite systems, challenges such as noise and data gaps must be effectively addressed. These quality losses are mainly caused by sensor artifacts, clouds, cloud shadows and other weather conditions (Verbesselt et al., 2012). In the context of time series analysis, harmonic modeling is a powerful tool to fill these gaps and reduce noise by smoothing the original signal (de Jong et al., 2011).

For that purpose, I created an R-package rHarmonics which enables the user to perform a harmonic analysis on a given time-series data set.

To calculate the harmonic fitted curve of a periodic signal, ordinary least squares regressions are computed using coupled sine and cosine curves on time-series data. The underlying algorithm which is based on Shumway & Stoffer (2017) equations 4.1 – 4.2 can be seen below:

MODIS NDVI time-series data together with a harmonic fitted curve using 3 cylces per year.



Literature:

de Jong, R., de Bruin, S., de Wit, A., Schaepman, M. E., & Dent, D. L. (2011). Analysis of monotonic greening and browning trends from global ndvi time-series. Remote Sensing of Environment, 115 (2), 692–702.

Shumway, R. H., & Stoffer, D. S. (2017). Time series analysis and its applications: with r examples. Springer.

Verbesselt, J., Zeileis, A., & Herold, M. (2012). Near real-time disturbance detection using satellite image time series. Remote Sensing of Environment, 123 , 98–108.

you may also like:

PhD defense by Thilo Erbertseder

PhD defense by Thilo Erbertseder

Thilo Erbertseder will defend his PhD thesis "Satellite-based analysis of NO2 air pollution: from global to urban aspect" on Wednesday 23rd of July at 3pm in John-Skilton Str. 4a, seminar room 1. All interested staff, students, family and friends are cordially invited...

upcoming PhD defense by Adomas Liepa

upcoming PhD defense by Adomas Liepa

Our PhD student Adomas Liepa will defend his Phd "Potential of Satellite Earth Observation in seasonal monitoring of complex agricultural environments of East Africa" on Thursday 24th of July at 11am. The defense will take place at John Skilton Str. 4a, seminarroom 1....

EORC Staff and EAGLE Students at ESA Living Planet Symposium 2025

EORC Staff and EAGLE Students at ESA Living Planet Symposium 2025

This week, our EORC team and EAGLE MSc students are joining the global Earth observation community at the ESA Living Planet Symposium (LPS) 2025 — one of the most important gatherings for Earth system scientists, remote sensing experts, and space agencies worldwide....

Our Contributions to the ESA Living Planet Symposium 2025

Our Contributions to the ESA Living Planet Symposium 2025

This week, the global Earth observation community gathered in Vienna for the ESA Living Planet Symposium 2025 — one of the most anticipated events for anyone passionate about understanding our planet through remote sensing. Our team was proud to contribute with an...