Skip to main content

Implementation of Empirical Mode Decomposition (EMD) and its variations

Project description

codecov BuildStatus DocStatus

PyEMD

The project is ongoing. This is very limited part of my private collection, but before I upload everything I want to make sure it works as it should. If there is something you wish to have, do email me as there is high chance that I have already done it, but it just sits around and waits until I’ll have more time. Don’t hesitate to contact me for anything.

Introduction

This is yet another Python implementation of Empirical Mode Decomposition (EMD). The package contains many EMD variations and intends to deliver more in time.

EMD variations:
  • Ensemble EMD (EEMD),

  • Image decomposotion (EMD2D),

  • different settings and configurations of vanilla EMD.

PyEMD allows to use different splines for envelopes, stopping criteria and extrema interpolation.

Available splines:
  • Natural cubic [default]

  • Pointwise cubic

  • Akima

  • Linear

Available stopping criteria:
  • Cauchy convergence [default]

  • Fixed number of iterations

  • Number of consecutive proto-imfs

Extrema detection:
  • Discrete extrema [default]

  • Parabolic interpolation

Installation

PyPi

Packaged obtained from PyPi is/will be slightly behind this project, so some features might not be the same. However, it seems to be the easiest/nicest way of installing any Python packages, so why not this one?

$ pip install EMD-signal

Example

More detailed examples are included in documentation.

EMD

In most cases default settings are enough. Simply import EMD and pass your signal to emd() method.

from PyEMD import EMD
import numpy as np

s = np.random.random(100)
emd = EMD()
IMFs = emd.emd(s)

The Figure below was produced with input: \(S(t) = cos(22 \pi t^2) + 6t^2\)

simpleExample

EEMD

Simplest case of using Esnembld EMD (EEMD) is by importing EEMD and passing your signal to eemd() method.

from PyEMD import EEMD
import numpy as np

s = np.random.random(100)
eemd = EEMD()
eIMFs = eemd.eemd(s)

EMD2D

Simplest case is to pass image as monochromatic numpy 2D array.

from PyEMD import EMD2D
import numpy as np

x, y = np.arange(128), np.arange(128).reshape((-1,1))
img = np.sin(0.1*x)*np.cos(0.2*y)
emd2d = EMD2D()
IMFs_2D = emd2d.emd(img)

Contact

Feel free to contact me with any questions, requests or simply saying hi. It’s always nice to know that I might have contributed to saving someone’s time or that I might improve my skills/projects.

Contact me either through gmail ({my_username}@gmail) or search me favourite web search.

Project details


Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Source Distribution

EMD-signal-0.2.2.tar.gz (34.9 kB view hashes)

Uploaded Source

Built Distribution

EMD_signal-0.2.2-py2.py3-none-any.whl (36.3 kB view hashes)

Uploaded Python 2 Python 3

Supported by

AWS AWS Cloud computing and Security Sponsor Datadog Datadog Monitoring Fastly Fastly CDN Google Google Download Analytics Microsoft Microsoft PSF Sponsor Pingdom Pingdom Monitoring Sentry Sentry Error logging StatusPage StatusPage Status page