PySDKit 0.4.43__tar.gz → 0.4.45__tar.gz

{pysdkit-0.4.43 → pysdkit-0.4.45}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: PySDKit
-Version: 0.4.43
+Version: 0.4.45
 Summary: A Python library for signal decomposition algorithms with a unified interface.
 Author: Whenxuan Wang, RuiZhe Wang, Rongkun Zhu, Kai Wu, Lei Wang, josefinez, Deeksha Manjunath, Yuan Feng, WenTong Zhao, JacktheFowler
 Maintainer-email: Whenxuan Wang <wwhenxuan@gmail.com>, RuiZhe Wang <3133986068@qq.com>
@@ -51,10 +51,10 @@ Requires-Dist: tqdm>=4.66.5
 Requires-Dist: requests>=2.32.3
 Dynamic: license-file
 
-# PySDKit: signal decomposition in Python
-
 <div align="center">
 
+# PySDKit: signal decomposition in Python
+
 [![PyPI version](https://badge.fury.io/py/PySDKit.svg)](https://pypi.org/project/PySDKit/) 
 ![License](https://img.shields.io/github/license/wwhenxuan/PySDKit)
 [![Python](https://img.shields.io/badge/python-3.8+-blue?logo=python)](https://www.python.org/)
@@ -91,60 +91,69 @@ This project integrates simple signal processing methods, signal decomposition a
 4. Visualize and analyze the original signal and the intrinsic mode functions IMFs obtained by decomposition.
 
 ~~~python
+import numpy as np
 from pysdkit import EMD
-from pysdkit.data import test_emd
 from pysdkit.plot import plot_IMFs
 
-t, signal = test_emd()
+t = np.linspace(0, 1, 1000)
+signal = np.sin(2 * np.pi * 5 * t) + 0.7 * np.sin(2 * np.pi * 25 * t) + 0.45 * np.sin(2 * np.pi * 80 * t)
 
-# create an instance for signal decomposition
 emd = EMD()
-# implement signal decomposition
-IMFs = emd.fit_transform(signal, max_imfs=2)
-plot_IMFs(signal, IMFs)
+IMFs = emd.fit_transform(signal, max_imfs=3)
+plot_IMFs(signal, IMFs, view="2d_freq", fs=1000, freq_max=150)
 ~~~
 
-![example](https://raw.githubusercontent.com/wwhenxuan/PySDKit/main/images/example.jpg)
+![example](https://raw.githubusercontent.com/wwhenxuan/PySDKit/main/images/decomposition_demo.jpg)
 
 The EMD in the above example is the most classic [`empirical mode decomposition`](https://www.mathworks.com/help/signal/ref/emd.html) algorithm in signal decomposition. For more complex signals, you can try other algorithms such as variational mode decomposition ([`VMD`](https://ieeexplore.ieee.org/abstract/document/6655981)).
 
 ~~~python
-import numpy as np
 from pysdkit import VMD
+from pysdkit.data import test_vmd
+from pysdkit.plot import plot_IMFs
 
-# load new signal
-signal = np.load("./example/example.npy")
-
-# use variational mode decomposition
-vmd = VMD(alpha=500, K=3, tau=0.0, tol=1e-9)
-IMFs = vmd.fit_transform(signal=signal)
-print(IMFs.shape)
+t, signal, fs = test_vmd()
 
-vmd.plot_IMFs(save_figure=True)
+vmd = VMD(alpha=2000, K=4, tau=0.0, tol=1e-7)
+IMFs = vmd.fit_transform(signal)
+plot_IMFs(signal, IMFs, view="2d_freq", fs=fs, freq_max=fs / 2)
 ~~~
 
 ![vmd_example](https://raw.githubusercontent.com/wwhenxuan/PySDKit/main/images/vmd_example.jpg)
 
-Better observe the characteristics of the decomposed intrinsic mode function in the frequency domain.
+For multichannel recordings, algorithms such as multivariate VMD ([`MVMD`](https://doi.org/10.1109/TSP.2019.2951223)) keep shared oscillations **mode-aligned** across channels:
 
 ~~~python
-from pysdkit.plot import plot_IMFs_amplitude_spectra
+import numpy as np
+from pysdkit import MVMD
+from pysdkit.plot import plot_IMFs
 
-# frequency domain visualization
-plot_IMFs_amplitude_spectra(IMFs, smooth="exp")   # use exp smooth
+t = np.arange(0, 1, 0.001)
+# ch1: 2+36 Hz, ch2: 24+36 Hz, ch3: 80+36 Hz (36 Hz shared)
+signal = np.vstack([
+    np.cos(2*np.pi*2*t) + np.cos(2*np.pi*36*t),
+    np.cos(2*np.pi*24*t) + np.cos(2*np.pi*36*t),
+    np.cos(2*np.pi*80*t) + np.cos(2*np.pi*36*t),
+])
+
+mvmd = MVMD(alpha=2000, K=4, tau=0.0, init="uniform")
+IMFs = mvmd.fit_transform(signal)   # shape: (K, T, C)
+plot_IMFs(signal, IMFs)             # per-channel panels
 ~~~
 
-![frequency_example](https://raw.githubusercontent.com/wwhenxuan/PySDKit/main/images/frequency_example.jpg)
+![mvmd_example](https://raw.githubusercontent.com/wwhenxuan/PySDKit/main/images/mvmd_example.jpg)
 
 ## Target 🎯 <a id="Target"></a>
 
 `PySDKit` is still under development. We are currently working on reproducing the signal decomposition algorithms in the table below, including not only common decomposition algorithms for `univariate signals` such as EMD, VMD, OVMD, AVNCMD, ALIF and APMD, but also decomposition algorithms for `multivariate signals` such as MEMD and MVMD. We will also further reproduce the decomposition algorithms for `two-dimensional images` to make PySDKit not only suitable for signal processing, but also for image analysis and understanding. See [`Mission`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/README.md) for the reasons why we developed PySDKit.
 
+
+
 |                          Algorithm                           |                            Paper                             |                            Source                            | State |
 | :----------------------------------------------------------: | :----------------------------------------------------------: | :----------------------------------------------------------: | :---: |
 | [`EMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_emd/emd.py) (Empirical Mode Decomposition) | [[paper]](https://royalsocietypublishing.org/doi/abs/10.1098/rspa.1998.0193) | [[code]](https://www.mathworks.com/help/signal/ref/emd.html) |   ✔️   |
 | [`MEMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_emd/memd.py) (Multivariate Empirical Mode Decomposition) | [[paper]](https://royalsocietypublishing.org/doi/full/10.1098/rspa.2009.0502) | [[code]](https://www.mathworks.com/matlabcentral/fileexchange/71270-fast-and-adaptive-multivariate-and-multidimensional-emd) |   ✔️   |
-| [`BEMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_emd2d/bemd.py) (Bidimensional Empirical Mode Decomposition) | [[paper]](https://www.sciencedirect.com/science/article/abs/pii/S0262885603000945) | [[code]](https://github.com/laszukdawid/PyEMD/blob/master/PyEMD/BEMD.py) |   ✖️   |
+| [`BEMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_emd2d/bemd.py) (Bidimensional Empirical Mode Decomposition) | [[paper]](https://www.sciencedirect.com/science/article/abs/pii/S0262885603000945) | [[code]](https://github.com/laszukdawid/PyEMD/blob/master/PyEMD/BEMD.py) |   ✔️   |
 |      [`CEMD`]() (Complex Empirical Mode Decomposition)       | [[paper]](https://ieeexplore.ieee.org/abstract/document/4063369) |                          [[code]]()                          |   ✖️   |
 | [`EEMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_emd/eemd.py) (Ensemble Empirical Mode Decomposition) | [[paper]](https://www.sciencedirect.com/topics/physics-and-astronomy/ensemble-empirical-mode-decomposition) | [[code]](https://github.com/laszukdawid/PyEMD/blob/master/PyEMD/EEMD.py) |   ✔️   |
 | [`REMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_emd/remd.py) (Robust Empirical Mode Decomposition) | [[paper]](https://www.sciencedirect.com/science/article/pii/S0019057821003785) | [[code]](https://www.mathworks.com/matlabcentral/fileexchange/70032-robust-empirical-mode-decomposition-remd) |   ✔️   |
@@ -159,7 +168,8 @@ plot_IMFs_amplitude_spectra(IMFs, smooth="exp")   # use exp smooth
 | [`HVD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_tid/hvd.py) (Hilbert Vibration Decomposition) | [[paper]](https://www.sciencedirect.com/science/article/pii/S0022460X06001556) | [[code]](https://www.mathworks.com/matlabcentral/fileexchange/178804-hilbert-vibration-decomposition?s_tid=FX_rc1_behav) |   ✔️   |
 | [`ITD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_tid/itd.py) (Intrinsic Time-Scale Decomposition) | [[paper]](https://royalsocietypublishing.org/doi/10.1098/rspa.2006.1761) | [[code]](https://www.mathworks.com/matlabcentral/fileexchange/69380-intrinsic-time-scale-decomposition-itd) |   ✔️   |
 | [`ALIF`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_tid/alif.py) (Adaptive Local Iterative Filtering) |          [[paper]](https://arxiv.org/abs/1411.6051)          |           [[code]](https://github.com/Cicone/ALIF)           |   ✔️   |
-| [`APMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_apmd/apmd.py) (Adaptive Polymorphic Mode Decomposition) |     [[paper]](https://doi.org/10.1016/j.dsp.2024.104913)     | [[code]](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_apmd/apmd.py) |   ✔️   |
+| [`IMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_imd/imd.py) (Impulsive Mode Decomposition) | [[paper]](https://www.sciencedirect.com/science/article/pii/S0888327024001250) | [[code]](https://www.mathworks.com/matlabcentral/fileexchange/182635-impulsive-mode-decomposition?s_tid=FX_rc3_behav) |   ✔️   |
+| [`APMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_imd/apmd.py) (Adaptive Polymorphic Mode Decomposition) |     [[paper]](https://doi.org/10.1016/j.dsp.2024.104913)     | [[code]](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_imd/apmd.py) |   ✔️   |
 | [`LMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_lmd/lmd.py) (Local Mean Decomposition) | [[paper]](https://www.sciencedirect.com/science/article/abs/pii/S1051200418308133?via%3Dihub) | [[code]](https://github.com/shownlin/PyLMD/blob/master/PyLMD/LMD.py) |   ✔️   |
 | [`RLMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_lmd/rlmd.py) (Robust Local Mean Decomposition) | [[paper]](https://www.sciencedirect.com/science/article/abs/pii/S0888327017301619) | [[code]](https://www.mathworks.com/matlabcentral/fileexchange/66935-robust-local-mean-decomposition-rlmd) |   ✔️   |
 | [`FMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_fmd/fmd.py) (Feature Mode Decomposition) |   [[paper]](https://ieeexplore.ieee.org/document/9732251)    | [[code]](https://www.mathworks.com/matlabcentral/fileexchange/108099-feature-mode-decomposition-fmd) |   ✔️   |
@@ -188,10 +198,10 @@ plot_IMFs_amplitude_spectra(IMFs, smooth="exp")   # use exp smooth
 | [`JMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_jmd/jmd.py) (Jump Plus AM-FM Mode Decomposition) |         [[paper]](https://arxiv.org/abs/2407.07800)          | [[code]](https://www.mathworks.com/matlabcentral/fileexchange/169388-jump-plus-am-fm-mode-decomposition-jmd?s_tid=prof_contriblnk) |   ️✔️   |
 | [`MJMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_jmd/mjmd.py) (Multivariate Jump Plus AM-FM Mode Decomposition) |         [[paper]](https://arxiv.org/abs/2407.07800)          | [[code]](https://www.mathworks.com/matlabcentral/fileexchange/169393-multivariate-jump-plus-am-fm-mode-decomposition-mjmd?s_tid=prof_contriblnk) |   ✔️   |
 | [`SJMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_jmd/sjmd.py) / `SMJMD` (Successive Jump and Mode Decomposition) |         [[paper]](https://arxiv.org/abs/2504.08453)          |                          [[code]]()                          |   ✔️   |
-|   [`ESMD`]() (Extreme-Point Symmetric Mode Decomposition)    |          [[paper]](https://arxiv.org/abs/1303.6540)          |         [[code]](https://github.com/WuShichao/esmd)          |   ✔️   |
+| [`ESMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_emd/esmd.py) (Extreme-Point Symmetric Mode Decomposition) |          [[paper]](https://arxiv.org/abs/1303.6540)          |         [[code]](https://github.com/WuShichao/esmd)          |   ✔️   |
 |   [`STNBMD`]() (Short-Time Narrow-Band Mode Decomposition)   | [[paper]](https://www.sciencedirect.com/science/article/pii/S0022460X16002443?via%3Dihub) | [[code]](https://ww2.mathworks.cn/matlabcentral/fileexchange/56226-short-time-narrow-band-mode-decomposition-stnbmd-toolbox) |   ✖️   |
 | [`SWD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_osd/swd.py) (Swarm Decomposition) |   [[paper]](https://doi.org/10.1016/j.sigpro.2016.09.004)    |  [[code]](https://github.com/gkaposto/Swarm-Decomposition)   |   ✔️   |
-| [`OSD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_osd/osd.py) (Optimization-based Signal Decomposition) | [[paper]](https://web.stanford.edu/~boyd/papers/sig_decomp_mprox.html) | [[code]](https://github.com/cvxgrp/signal-decomposition) |   ✔️   |
+| [`OSD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_osd/osd.py) (Optimization-based Signal Decomposition) | [[paper]](https://web.stanford.edu/~boyd/papers/sig_decomp_mprox.html) |   [[code]](https://github.com/cvxgrp/signal-decomposition)   |   ✔️   |
 | [`STL`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/tsa/_stl.py) (Seasonal-Trend decomposition using LOESS) | [[paper]](https://www.nniiem.ru/file/news/2016/stl-statistical-model.pdf) | [[code]](https://www.statsmodels.org/stable/examples/notebooks/generated/stl_decomposition.html) |   ✔️   |
 | [`MSTL`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/tsa/_mstl.py) (Multivariate Seasonal-Trend decomposition using LOESS) |         [[paper]](https://arxiv.org/abs/2107.13462)          |         [[code]](https://github.com/KishManani/MSTL)         |   ✔️   |
 

{pysdkit-0.4.43 → pysdkit-0.4.45}/PySDKit.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: PySDKit
-Version: 0.4.43
+Version: 0.4.45
 Summary: A Python library for signal decomposition algorithms with a unified interface.
 Author: Whenxuan Wang, RuiZhe Wang, Rongkun Zhu, Kai Wu, Lei Wang, josefinez, Deeksha Manjunath, Yuan Feng, WenTong Zhao, JacktheFowler
 Maintainer-email: Whenxuan Wang <wwhenxuan@gmail.com>, RuiZhe Wang <3133986068@qq.com>
@@ -51,10 +51,10 @@ Requires-Dist: tqdm>=4.66.5
 Requires-Dist: requests>=2.32.3
 Dynamic: license-file
 
-# PySDKit: signal decomposition in Python
-
 <div align="center">
 
+# PySDKit: signal decomposition in Python
+
 [![PyPI version](https://badge.fury.io/py/PySDKit.svg)](https://pypi.org/project/PySDKit/) 
 ![License](https://img.shields.io/github/license/wwhenxuan/PySDKit)
 [![Python](https://img.shields.io/badge/python-3.8+-blue?logo=python)](https://www.python.org/)
@@ -91,60 +91,69 @@ This project integrates simple signal processing methods, signal decomposition a
 4. Visualize and analyze the original signal and the intrinsic mode functions IMFs obtained by decomposition.
 
 ~~~python
+import numpy as np
 from pysdkit import EMD
-from pysdkit.data import test_emd
 from pysdkit.plot import plot_IMFs
 
-t, signal = test_emd()
+t = np.linspace(0, 1, 1000)
+signal = np.sin(2 * np.pi * 5 * t) + 0.7 * np.sin(2 * np.pi * 25 * t) + 0.45 * np.sin(2 * np.pi * 80 * t)
 
-# create an instance for signal decomposition
 emd = EMD()
-# implement signal decomposition
-IMFs = emd.fit_transform(signal, max_imfs=2)
-plot_IMFs(signal, IMFs)
+IMFs = emd.fit_transform(signal, max_imfs=3)
+plot_IMFs(signal, IMFs, view="2d_freq", fs=1000, freq_max=150)
 ~~~
 
-![example](https://raw.githubusercontent.com/wwhenxuan/PySDKit/main/images/example.jpg)
+![example](https://raw.githubusercontent.com/wwhenxuan/PySDKit/main/images/decomposition_demo.jpg)
 
 The EMD in the above example is the most classic [`empirical mode decomposition`](https://www.mathworks.com/help/signal/ref/emd.html) algorithm in signal decomposition. For more complex signals, you can try other algorithms such as variational mode decomposition ([`VMD`](https://ieeexplore.ieee.org/abstract/document/6655981)).
 
 ~~~python
-import numpy as np
 from pysdkit import VMD
+from pysdkit.data import test_vmd
+from pysdkit.plot import plot_IMFs
 
-# load new signal
-signal = np.load("./example/example.npy")
-
-# use variational mode decomposition
-vmd = VMD(alpha=500, K=3, tau=0.0, tol=1e-9)
-IMFs = vmd.fit_transform(signal=signal)
-print(IMFs.shape)
+t, signal, fs = test_vmd()
 
-vmd.plot_IMFs(save_figure=True)
+vmd = VMD(alpha=2000, K=4, tau=0.0, tol=1e-7)
+IMFs = vmd.fit_transform(signal)
+plot_IMFs(signal, IMFs, view="2d_freq", fs=fs, freq_max=fs / 2)
 ~~~
 
 ![vmd_example](https://raw.githubusercontent.com/wwhenxuan/PySDKit/main/images/vmd_example.jpg)
 
-Better observe the characteristics of the decomposed intrinsic mode function in the frequency domain.
+For multichannel recordings, algorithms such as multivariate VMD ([`MVMD`](https://doi.org/10.1109/TSP.2019.2951223)) keep shared oscillations **mode-aligned** across channels:
 
 ~~~python
-from pysdkit.plot import plot_IMFs_amplitude_spectra
+import numpy as np
+from pysdkit import MVMD
+from pysdkit.plot import plot_IMFs
 
-# frequency domain visualization
-plot_IMFs_amplitude_spectra(IMFs, smooth="exp")   # use exp smooth
+t = np.arange(0, 1, 0.001)
+# ch1: 2+36 Hz, ch2: 24+36 Hz, ch3: 80+36 Hz (36 Hz shared)
+signal = np.vstack([
+    np.cos(2*np.pi*2*t) + np.cos(2*np.pi*36*t),
+    np.cos(2*np.pi*24*t) + np.cos(2*np.pi*36*t),
+    np.cos(2*np.pi*80*t) + np.cos(2*np.pi*36*t),
+])
+
+mvmd = MVMD(alpha=2000, K=4, tau=0.0, init="uniform")
+IMFs = mvmd.fit_transform(signal)   # shape: (K, T, C)
+plot_IMFs(signal, IMFs)             # per-channel panels
 ~~~
 
-![frequency_example](https://raw.githubusercontent.com/wwhenxuan/PySDKit/main/images/frequency_example.jpg)
+![mvmd_example](https://raw.githubusercontent.com/wwhenxuan/PySDKit/main/images/mvmd_example.jpg)
 
 ## Target 🎯 <a id="Target"></a>
 
 `PySDKit` is still under development. We are currently working on reproducing the signal decomposition algorithms in the table below, including not only common decomposition algorithms for `univariate signals` such as EMD, VMD, OVMD, AVNCMD, ALIF and APMD, but also decomposition algorithms for `multivariate signals` such as MEMD and MVMD. We will also further reproduce the decomposition algorithms for `two-dimensional images` to make PySDKit not only suitable for signal processing, but also for image analysis and understanding. See [`Mission`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/README.md) for the reasons why we developed PySDKit.
 
+
+
 |                          Algorithm                           |                            Paper                             |                            Source                            | State |
 | :----------------------------------------------------------: | :----------------------------------------------------------: | :----------------------------------------------------------: | :---: |
 | [`EMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_emd/emd.py) (Empirical Mode Decomposition) | [[paper]](https://royalsocietypublishing.org/doi/abs/10.1098/rspa.1998.0193) | [[code]](https://www.mathworks.com/help/signal/ref/emd.html) |   ✔️   |
 | [`MEMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_emd/memd.py) (Multivariate Empirical Mode Decomposition) | [[paper]](https://royalsocietypublishing.org/doi/full/10.1098/rspa.2009.0502) | [[code]](https://www.mathworks.com/matlabcentral/fileexchange/71270-fast-and-adaptive-multivariate-and-multidimensional-emd) |   ✔️   |
-| [`BEMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_emd2d/bemd.py) (Bidimensional Empirical Mode Decomposition) | [[paper]](https://www.sciencedirect.com/science/article/abs/pii/S0262885603000945) | [[code]](https://github.com/laszukdawid/PyEMD/blob/master/PyEMD/BEMD.py) |   ✖️   |
+| [`BEMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_emd2d/bemd.py) (Bidimensional Empirical Mode Decomposition) | [[paper]](https://www.sciencedirect.com/science/article/abs/pii/S0262885603000945) | [[code]](https://github.com/laszukdawid/PyEMD/blob/master/PyEMD/BEMD.py) |   ✔️   |
 |      [`CEMD`]() (Complex Empirical Mode Decomposition)       | [[paper]](https://ieeexplore.ieee.org/abstract/document/4063369) |                          [[code]]()                          |   ✖️   |
 | [`EEMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_emd/eemd.py) (Ensemble Empirical Mode Decomposition) | [[paper]](https://www.sciencedirect.com/topics/physics-and-astronomy/ensemble-empirical-mode-decomposition) | [[code]](https://github.com/laszukdawid/PyEMD/blob/master/PyEMD/EEMD.py) |   ✔️   |
 | [`REMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_emd/remd.py) (Robust Empirical Mode Decomposition) | [[paper]](https://www.sciencedirect.com/science/article/pii/S0019057821003785) | [[code]](https://www.mathworks.com/matlabcentral/fileexchange/70032-robust-empirical-mode-decomposition-remd) |   ✔️   |
@@ -159,7 +168,8 @@ plot_IMFs_amplitude_spectra(IMFs, smooth="exp")   # use exp smooth
 | [`HVD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_tid/hvd.py) (Hilbert Vibration Decomposition) | [[paper]](https://www.sciencedirect.com/science/article/pii/S0022460X06001556) | [[code]](https://www.mathworks.com/matlabcentral/fileexchange/178804-hilbert-vibration-decomposition?s_tid=FX_rc1_behav) |   ✔️   |
 | [`ITD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_tid/itd.py) (Intrinsic Time-Scale Decomposition) | [[paper]](https://royalsocietypublishing.org/doi/10.1098/rspa.2006.1761) | [[code]](https://www.mathworks.com/matlabcentral/fileexchange/69380-intrinsic-time-scale-decomposition-itd) |   ✔️   |
 | [`ALIF`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_tid/alif.py) (Adaptive Local Iterative Filtering) |          [[paper]](https://arxiv.org/abs/1411.6051)          |           [[code]](https://github.com/Cicone/ALIF)           |   ✔️   |
-| [`APMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_apmd/apmd.py) (Adaptive Polymorphic Mode Decomposition) |     [[paper]](https://doi.org/10.1016/j.dsp.2024.104913)     | [[code]](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_apmd/apmd.py) |   ✔️   |
+| [`IMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_imd/imd.py) (Impulsive Mode Decomposition) | [[paper]](https://www.sciencedirect.com/science/article/pii/S0888327024001250) | [[code]](https://www.mathworks.com/matlabcentral/fileexchange/182635-impulsive-mode-decomposition?s_tid=FX_rc3_behav) |   ✔️   |
+| [`APMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_imd/apmd.py) (Adaptive Polymorphic Mode Decomposition) |     [[paper]](https://doi.org/10.1016/j.dsp.2024.104913)     | [[code]](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_imd/apmd.py) |   ✔️   |
 | [`LMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_lmd/lmd.py) (Local Mean Decomposition) | [[paper]](https://www.sciencedirect.com/science/article/abs/pii/S1051200418308133?via%3Dihub) | [[code]](https://github.com/shownlin/PyLMD/blob/master/PyLMD/LMD.py) |   ✔️   |
 | [`RLMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_lmd/rlmd.py) (Robust Local Mean Decomposition) | [[paper]](https://www.sciencedirect.com/science/article/abs/pii/S0888327017301619) | [[code]](https://www.mathworks.com/matlabcentral/fileexchange/66935-robust-local-mean-decomposition-rlmd) |   ✔️   |
 | [`FMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_fmd/fmd.py) (Feature Mode Decomposition) |   [[paper]](https://ieeexplore.ieee.org/document/9732251)    | [[code]](https://www.mathworks.com/matlabcentral/fileexchange/108099-feature-mode-decomposition-fmd) |   ✔️   |
@@ -188,10 +198,10 @@ plot_IMFs_amplitude_spectra(IMFs, smooth="exp")   # use exp smooth
 | [`JMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_jmd/jmd.py) (Jump Plus AM-FM Mode Decomposition) |         [[paper]](https://arxiv.org/abs/2407.07800)          | [[code]](https://www.mathworks.com/matlabcentral/fileexchange/169388-jump-plus-am-fm-mode-decomposition-jmd?s_tid=prof_contriblnk) |   ️✔️   |
 | [`MJMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_jmd/mjmd.py) (Multivariate Jump Plus AM-FM Mode Decomposition) |         [[paper]](https://arxiv.org/abs/2407.07800)          | [[code]](https://www.mathworks.com/matlabcentral/fileexchange/169393-multivariate-jump-plus-am-fm-mode-decomposition-mjmd?s_tid=prof_contriblnk) |   ✔️   |
 | [`SJMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_jmd/sjmd.py) / `SMJMD` (Successive Jump and Mode Decomposition) |         [[paper]](https://arxiv.org/abs/2504.08453)          |                          [[code]]()                          |   ✔️   |
-|   [`ESMD`]() (Extreme-Point Symmetric Mode Decomposition)    |          [[paper]](https://arxiv.org/abs/1303.6540)          |         [[code]](https://github.com/WuShichao/esmd)          |   ✔️   |
+| [`ESMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_emd/esmd.py) (Extreme-Point Symmetric Mode Decomposition) |          [[paper]](https://arxiv.org/abs/1303.6540)          |         [[code]](https://github.com/WuShichao/esmd)          |   ✔️   |
 |   [`STNBMD`]() (Short-Time Narrow-Band Mode Decomposition)   | [[paper]](https://www.sciencedirect.com/science/article/pii/S0022460X16002443?via%3Dihub) | [[code]](https://ww2.mathworks.cn/matlabcentral/fileexchange/56226-short-time-narrow-band-mode-decomposition-stnbmd-toolbox) |   ✖️   |
 | [`SWD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_osd/swd.py) (Swarm Decomposition) |   [[paper]](https://doi.org/10.1016/j.sigpro.2016.09.004)    |  [[code]](https://github.com/gkaposto/Swarm-Decomposition)   |   ✔️   |
-| [`OSD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_osd/osd.py) (Optimization-based Signal Decomposition) | [[paper]](https://web.stanford.edu/~boyd/papers/sig_decomp_mprox.html) | [[code]](https://github.com/cvxgrp/signal-decomposition) |   ✔️   |
+| [`OSD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_osd/osd.py) (Optimization-based Signal Decomposition) | [[paper]](https://web.stanford.edu/~boyd/papers/sig_decomp_mprox.html) |   [[code]](https://github.com/cvxgrp/signal-decomposition)   |   ✔️   |
 | [`STL`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/tsa/_stl.py) (Seasonal-Trend decomposition using LOESS) | [[paper]](https://www.nniiem.ru/file/news/2016/stl-statistical-model.pdf) | [[code]](https://www.statsmodels.org/stable/examples/notebooks/generated/stl_decomposition.html) |   ✔️   |
 | [`MSTL`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/tsa/_mstl.py) (Multivariate Seasonal-Trend decomposition using LOESS) |         [[paper]](https://arxiv.org/abs/2107.13462)          |         [[code]](https://github.com/KishManani/MSTL)         |   ✔️   |
 

{pysdkit-0.4.43 → pysdkit-0.4.45}/PySDKit.egg-info/SOURCES.txt

@@ -7,8 +7,6 @@ PySDKit.egg-info/dependency_links.txt
 PySDKit.egg-info/requires.txt
 PySDKit.egg-info/top_level.txt
 pysdkit/__init__.py
-pysdkit/_apmd/__init__.py
-pysdkit/_apmd/apmd.py
 pysdkit/_emd/__init__.py
 pysdkit/_emd/_find_extrema.py
 pysdkit/_emd/_prepare_points.py
@@ -17,6 +15,7 @@ pysdkit/_emd/ceemdan.py
 pysdkit/_emd/eemd.py
 pysdkit/_emd/efd.py
 pysdkit/_emd/emd.py
+pysdkit/_emd/esmd.py
 pysdkit/_emd/memd.py
 pysdkit/_emd/remd.py
 pysdkit/_emd/semd.py
@@ -28,8 +27,6 @@ pysdkit/_emd2d/__init__.py
 pysdkit/_emd2d/bemd.py
 pysdkit/_emd2d/bmemd.py
 pysdkit/_emd2d/emd2d.py
-pysdkit/_esmd/__init__.py
-pysdkit/_esmd/esmd.py
 pysdkit/_ewt/__init__.py
 pysdkit/_ewt/ewt.py
 pysdkit/_ewt/ewt2d.py
@@ -47,6 +44,11 @@ pysdkit/_gdmd/agncmd.py
 pysdkit/_gdmd/gdmd.py
 pysdkit/_gdmd/ivgnmd.py
 pysdkit/_gdmd/vgnmd.py
+pysdkit/_imd/__init__.py
+pysdkit/_imd/apmd.py
+pysdkit/_imd/imd.py
+pysdkit/_imd/data/gearbox_fault_snippet.npy
+pysdkit/_imd/data/input_sig.npy
 pysdkit/_jmd/__init__.py
 pysdkit/_jmd/jmd.py
 pysdkit/_jmd/mjmd.py
@@ -93,7 +95,9 @@ pysdkit/data/_image.py
 pysdkit/data/_models.py
 pysdkit/data/_test_univariate.py
 pysdkit/data/_time_series.py
+pysdkit/data/fmd_demo.npy
 pysdkit/data/texture.txt
+pysdkit/data/vmd_example.npy
 pysdkit/entropy/__init__.py
 pysdkit/entropy/_approxiamte_entropy.py
 pysdkit/entropy/_permutation_entropy.py
@@ -126,7 +130,10 @@ pysdkit/tests/test_ewt2d.py
 pysdkit/tests/test_faemd.py
 pysdkit/tests/test_faemd2d.py
 pysdkit/tests/test_faemd3d.py
+pysdkit/tests/test_fmd.py
 pysdkit/tests/test_gdmd.py
+pysdkit/tests/test_hvd.py
+pysdkit/tests/test_imd.py
 pysdkit/tests/test_itd.py
 pysdkit/tests/test_ivgnmd.py
 pysdkit/tests/test_jmd.py

{pysdkit-0.4.43 → pysdkit-0.4.45}/README.md

@@ -1,7 +1,7 @@
-# PySDKit: signal decomposition in Python
-
 <div align="center">
 
+# PySDKit: signal decomposition in Python
+
 [![PyPI version](https://badge.fury.io/py/PySDKit.svg)](https://pypi.org/project/PySDKit/) 
 ![License](https://img.shields.io/github/license/wwhenxuan/PySDKit)
 [![Python](https://img.shields.io/badge/python-3.8+-blue?logo=python)](https://www.python.org/)
@@ -38,60 +38,69 @@ This project integrates simple signal processing methods, signal decomposition a
 4. Visualize and analyze the original signal and the intrinsic mode functions IMFs obtained by decomposition.
 
 ~~~python
+import numpy as np
 from pysdkit import EMD
-from pysdkit.data import test_emd
 from pysdkit.plot import plot_IMFs
 
-t, signal = test_emd()
+t = np.linspace(0, 1, 1000)
+signal = np.sin(2 * np.pi * 5 * t) + 0.7 * np.sin(2 * np.pi * 25 * t) + 0.45 * np.sin(2 * np.pi * 80 * t)
 
-# create an instance for signal decomposition
 emd = EMD()
-# implement signal decomposition
-IMFs = emd.fit_transform(signal, max_imfs=2)
-plot_IMFs(signal, IMFs)
+IMFs = emd.fit_transform(signal, max_imfs=3)
+plot_IMFs(signal, IMFs, view="2d_freq", fs=1000, freq_max=150)
 ~~~
 
-![example](https://raw.githubusercontent.com/wwhenxuan/PySDKit/main/images/example.jpg)
+![example](https://raw.githubusercontent.com/wwhenxuan/PySDKit/main/images/decomposition_demo.jpg)
 
 The EMD in the above example is the most classic [`empirical mode decomposition`](https://www.mathworks.com/help/signal/ref/emd.html) algorithm in signal decomposition. For more complex signals, you can try other algorithms such as variational mode decomposition ([`VMD`](https://ieeexplore.ieee.org/abstract/document/6655981)).
 
 ~~~python
-import numpy as np
 from pysdkit import VMD
+from pysdkit.data import test_vmd
+from pysdkit.plot import plot_IMFs
 
-# load new signal
-signal = np.load("./example/example.npy")
-
-# use variational mode decomposition
-vmd = VMD(alpha=500, K=3, tau=0.0, tol=1e-9)
-IMFs = vmd.fit_transform(signal=signal)
-print(IMFs.shape)
+t, signal, fs = test_vmd()
 
-vmd.plot_IMFs(save_figure=True)
+vmd = VMD(alpha=2000, K=4, tau=0.0, tol=1e-7)
+IMFs = vmd.fit_transform(signal)
+plot_IMFs(signal, IMFs, view="2d_freq", fs=fs, freq_max=fs / 2)
 ~~~
 
 ![vmd_example](https://raw.githubusercontent.com/wwhenxuan/PySDKit/main/images/vmd_example.jpg)
 
-Better observe the characteristics of the decomposed intrinsic mode function in the frequency domain.
+For multichannel recordings, algorithms such as multivariate VMD ([`MVMD`](https://doi.org/10.1109/TSP.2019.2951223)) keep shared oscillations **mode-aligned** across channels:
 
 ~~~python
-from pysdkit.plot import plot_IMFs_amplitude_spectra
+import numpy as np
+from pysdkit import MVMD
+from pysdkit.plot import plot_IMFs
 
-# frequency domain visualization
-plot_IMFs_amplitude_spectra(IMFs, smooth="exp")   # use exp smooth
+t = np.arange(0, 1, 0.001)
+# ch1: 2+36 Hz, ch2: 24+36 Hz, ch3: 80+36 Hz (36 Hz shared)
+signal = np.vstack([
+    np.cos(2*np.pi*2*t) + np.cos(2*np.pi*36*t),
+    np.cos(2*np.pi*24*t) + np.cos(2*np.pi*36*t),
+    np.cos(2*np.pi*80*t) + np.cos(2*np.pi*36*t),
+])
+
+mvmd = MVMD(alpha=2000, K=4, tau=0.0, init="uniform")
+IMFs = mvmd.fit_transform(signal)   # shape: (K, T, C)
+plot_IMFs(signal, IMFs)             # per-channel panels
 ~~~
 
-![frequency_example](https://raw.githubusercontent.com/wwhenxuan/PySDKit/main/images/frequency_example.jpg)
+![mvmd_example](https://raw.githubusercontent.com/wwhenxuan/PySDKit/main/images/mvmd_example.jpg)
 
 ## Target 🎯 <a id="Target"></a>
 
 `PySDKit` is still under development. We are currently working on reproducing the signal decomposition algorithms in the table below, including not only common decomposition algorithms for `univariate signals` such as EMD, VMD, OVMD, AVNCMD, ALIF and APMD, but also decomposition algorithms for `multivariate signals` such as MEMD and MVMD. We will also further reproduce the decomposition algorithms for `two-dimensional images` to make PySDKit not only suitable for signal processing, but also for image analysis and understanding. See [`Mission`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/README.md) for the reasons why we developed PySDKit.
 
+
+
 |                          Algorithm                           |                            Paper                             |                            Source                            | State |
 | :----------------------------------------------------------: | :----------------------------------------------------------: | :----------------------------------------------------------: | :---: |
 | [`EMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_emd/emd.py) (Empirical Mode Decomposition) | [[paper]](https://royalsocietypublishing.org/doi/abs/10.1098/rspa.1998.0193) | [[code]](https://www.mathworks.com/help/signal/ref/emd.html) |   ✔️   |
 | [`MEMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_emd/memd.py) (Multivariate Empirical Mode Decomposition) | [[paper]](https://royalsocietypublishing.org/doi/full/10.1098/rspa.2009.0502) | [[code]](https://www.mathworks.com/matlabcentral/fileexchange/71270-fast-and-adaptive-multivariate-and-multidimensional-emd) |   ✔️   |
-| [`BEMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_emd2d/bemd.py) (Bidimensional Empirical Mode Decomposition) | [[paper]](https://www.sciencedirect.com/science/article/abs/pii/S0262885603000945) | [[code]](https://github.com/laszukdawid/PyEMD/blob/master/PyEMD/BEMD.py) |   ✖️   |
+| [`BEMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_emd2d/bemd.py) (Bidimensional Empirical Mode Decomposition) | [[paper]](https://www.sciencedirect.com/science/article/abs/pii/S0262885603000945) | [[code]](https://github.com/laszukdawid/PyEMD/blob/master/PyEMD/BEMD.py) |   ✔️   |
 |      [`CEMD`]() (Complex Empirical Mode Decomposition)       | [[paper]](https://ieeexplore.ieee.org/abstract/document/4063369) |                          [[code]]()                          |   ✖️   |
 | [`EEMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_emd/eemd.py) (Ensemble Empirical Mode Decomposition) | [[paper]](https://www.sciencedirect.com/topics/physics-and-astronomy/ensemble-empirical-mode-decomposition) | [[code]](https://github.com/laszukdawid/PyEMD/blob/master/PyEMD/EEMD.py) |   ✔️   |
 | [`REMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_emd/remd.py) (Robust Empirical Mode Decomposition) | [[paper]](https://www.sciencedirect.com/science/article/pii/S0019057821003785) | [[code]](https://www.mathworks.com/matlabcentral/fileexchange/70032-robust-empirical-mode-decomposition-remd) |   ✔️   |
@@ -106,7 +115,8 @@ plot_IMFs_amplitude_spectra(IMFs, smooth="exp")   # use exp smooth
 | [`HVD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_tid/hvd.py) (Hilbert Vibration Decomposition) | [[paper]](https://www.sciencedirect.com/science/article/pii/S0022460X06001556) | [[code]](https://www.mathworks.com/matlabcentral/fileexchange/178804-hilbert-vibration-decomposition?s_tid=FX_rc1_behav) |   ✔️   |
 | [`ITD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_tid/itd.py) (Intrinsic Time-Scale Decomposition) | [[paper]](https://royalsocietypublishing.org/doi/10.1098/rspa.2006.1761) | [[code]](https://www.mathworks.com/matlabcentral/fileexchange/69380-intrinsic-time-scale-decomposition-itd) |   ✔️   |
 | [`ALIF`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_tid/alif.py) (Adaptive Local Iterative Filtering) |          [[paper]](https://arxiv.org/abs/1411.6051)          |           [[code]](https://github.com/Cicone/ALIF)           |   ✔️   |
-| [`APMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_apmd/apmd.py) (Adaptive Polymorphic Mode Decomposition) |     [[paper]](https://doi.org/10.1016/j.dsp.2024.104913)     | [[code]](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_apmd/apmd.py) |   ✔️   |
+| [`IMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_imd/imd.py) (Impulsive Mode Decomposition) | [[paper]](https://www.sciencedirect.com/science/article/pii/S0888327024001250) | [[code]](https://www.mathworks.com/matlabcentral/fileexchange/182635-impulsive-mode-decomposition?s_tid=FX_rc3_behav) |   ✔️   |
+| [`APMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_imd/apmd.py) (Adaptive Polymorphic Mode Decomposition) |     [[paper]](https://doi.org/10.1016/j.dsp.2024.104913)     | [[code]](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_imd/apmd.py) |   ✔️   |
 | [`LMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_lmd/lmd.py) (Local Mean Decomposition) | [[paper]](https://www.sciencedirect.com/science/article/abs/pii/S1051200418308133?via%3Dihub) | [[code]](https://github.com/shownlin/PyLMD/blob/master/PyLMD/LMD.py) |   ✔️   |
 | [`RLMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_lmd/rlmd.py) (Robust Local Mean Decomposition) | [[paper]](https://www.sciencedirect.com/science/article/abs/pii/S0888327017301619) | [[code]](https://www.mathworks.com/matlabcentral/fileexchange/66935-robust-local-mean-decomposition-rlmd) |   ✔️   |
 | [`FMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_fmd/fmd.py) (Feature Mode Decomposition) |   [[paper]](https://ieeexplore.ieee.org/document/9732251)    | [[code]](https://www.mathworks.com/matlabcentral/fileexchange/108099-feature-mode-decomposition-fmd) |   ✔️   |
@@ -135,10 +145,10 @@ plot_IMFs_amplitude_spectra(IMFs, smooth="exp")   # use exp smooth
 | [`JMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_jmd/jmd.py) (Jump Plus AM-FM Mode Decomposition) |         [[paper]](https://arxiv.org/abs/2407.07800)          | [[code]](https://www.mathworks.com/matlabcentral/fileexchange/169388-jump-plus-am-fm-mode-decomposition-jmd?s_tid=prof_contriblnk) |   ️✔️   |
 | [`MJMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_jmd/mjmd.py) (Multivariate Jump Plus AM-FM Mode Decomposition) |         [[paper]](https://arxiv.org/abs/2407.07800)          | [[code]](https://www.mathworks.com/matlabcentral/fileexchange/169393-multivariate-jump-plus-am-fm-mode-decomposition-mjmd?s_tid=prof_contriblnk) |   ✔️   |
 | [`SJMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_jmd/sjmd.py) / `SMJMD` (Successive Jump and Mode Decomposition) |         [[paper]](https://arxiv.org/abs/2504.08453)          |                          [[code]]()                          |   ✔️   |
-|   [`ESMD`]() (Extreme-Point Symmetric Mode Decomposition)    |          [[paper]](https://arxiv.org/abs/1303.6540)          |         [[code]](https://github.com/WuShichao/esmd)          |   ✔️   |
+| [`ESMD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_emd/esmd.py) (Extreme-Point Symmetric Mode Decomposition) |          [[paper]](https://arxiv.org/abs/1303.6540)          |         [[code]](https://github.com/WuShichao/esmd)          |   ✔️   |
 |   [`STNBMD`]() (Short-Time Narrow-Band Mode Decomposition)   | [[paper]](https://www.sciencedirect.com/science/article/pii/S0022460X16002443?via%3Dihub) | [[code]](https://ww2.mathworks.cn/matlabcentral/fileexchange/56226-short-time-narrow-band-mode-decomposition-stnbmd-toolbox) |   ✖️   |
 | [`SWD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_osd/swd.py) (Swarm Decomposition) |   [[paper]](https://doi.org/10.1016/j.sigpro.2016.09.004)    |  [[code]](https://github.com/gkaposto/Swarm-Decomposition)   |   ✔️   |
-| [`OSD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_osd/osd.py) (Optimization-based Signal Decomposition) | [[paper]](https://web.stanford.edu/~boyd/papers/sig_decomp_mprox.html) | [[code]](https://github.com/cvxgrp/signal-decomposition) |   ✔️   |
+| [`OSD`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/_osd/osd.py) (Optimization-based Signal Decomposition) | [[paper]](https://web.stanford.edu/~boyd/papers/sig_decomp_mprox.html) |   [[code]](https://github.com/cvxgrp/signal-decomposition)   |   ✔️   |
 | [`STL`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/tsa/_stl.py) (Seasonal-Trend decomposition using LOESS) | [[paper]](https://www.nniiem.ru/file/news/2016/stl-statistical-model.pdf) | [[code]](https://www.statsmodels.org/stable/examples/notebooks/generated/stl_decomposition.html) |   ✔️   |
 | [`MSTL`](https://github.com/wwhenxuan/PySDKit/blob/main/pysdkit/tsa/_mstl.py) (Multivariate Seasonal-Trend decomposition using LOESS) |         [[paper]](https://arxiv.org/abs/2107.13462)          |         [[code]](https://github.com/KishManani/MSTL)         |   ✔️   |
 

{pysdkit-0.4.43 → pysdkit-0.4.45}/pyproject.toml

@@ -68,5 +68,8 @@ include = ["pysdkit*"]
 
 [tool.setuptools.package-data]
 "*" = ["*.txt", "*.npy"]
+"pysdkit.data" = ["*.txt", "*.npy"]
 "pysdkit._tid" = ["data/*.npy"]
 "pysdkit._tid.data" = ["*.npy"]
+"pysdkit._imd" = ["data/*.npy"]
+"pysdkit._imd.data" = ["*.npy"]

{pysdkit-0.4.43 → pysdkit-0.4.45}/pysdkit/__init__.py

@@ -2,7 +2,7 @@
 A Python library for signal decomposition algorithms.
 """
 
-__version__ = "0.4.43"
+__version__ = "0.4.45"
 
 # Empirical Mode Decomposition
 from ._emd import EMD
@@ -44,7 +44,7 @@ from ._tid import HVD
 from ._tid import ITD
 
 # Extreme-Point Symmetric Mode Decomposition
-from ._esmd import ESMD, esmd
+from ._emd import ESMD, esmd
 
 # Local Mean Decomposition
 from ._lmd import LMD
@@ -112,8 +112,8 @@ from ._vncmd import AVNCMD
 # Adaptive Local Iterative Filtering (time iterative decomposition)
 from ._tid import ALIF
 
-# Adaptive Polymorphic Mode Decomposition
-from ._apmd import APMD
+# Adaptive Polymorphic Mode Decomposition / Impulsive Mode Decomposition
+from ._imd import APMD, IMD, imd
 
 # Empirical Wavelet Transform
 from ._ewt import ewt, EWT
@@ -207,6 +207,7 @@ Iterative Nonlinear Chirp Mode Decomposition    |    INCMD
 Adaptive Variational Nonlinear Chirp Mode Dec.  |    AVNCMD
 Adaptive Local Iterative Filtering              |    ALIF
 Adaptive Polymorphic Mode Decomposition         |    APMD
+Impulsive Mode Decomposition                    |    IMD
 Empirical Wavelet Transform                     |    EWT
 Empirical Wavelet Transform for 2D Image        |    EWT2D
 Jump Plus AM-FM Mode Decomposition              |    JMD
@@ -246,7 +247,8 @@ __all__ = [
     "SWD",
     "swd",
     "OSD",
-    "GDMD",    "gdmd",
+    "GDMD",
+    "gdmd",
     "VGNMD",
     "vgnmd",
     "IVGNMD",
@@ -272,6 +274,8 @@ __all__ = [
     "AVNCMD",
     "ALIF",
     "APMD",
+    "IMD",
+    "imd",
     "ewt",
     "EWT",
     "ewt2d",

{pysdkit-0.4.43 → pysdkit-0.4.45}/pysdkit/_emd/__init__.py

@@ -26,3 +26,5 @@ from .efd import EFD
 from .hht import HHT
 
 from .semd import SEMD, concatenate_signals, deconcatenate_imfs, transition_bridge
+
+from .esmd import ESMD, esmd

pysdkit-0.4.45/pysdkit/_fmd/__init__.py

@@ -0,0 +1,8 @@
+# -*- coding: utf-8 -*-
+"""
+Feature Mode Decomposition (FMD / NFMD).
+"""
+
+from .fmd import CK, FMD, TT, max_IJ
+
+__all__ = ["FMD", "TT", "CK", "max_IJ"]