.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "emd_tutorials/04_utilities/emd_tutorial_04_utils_02_logger.py" .. LINE NUMBERS ARE GIVEN BELOW. .. only:: html .. note:: :class: sphx-glr-download-link-note :ref:`Go to the end ` to download the full example code .. rst-class:: sphx-glr-example-title .. _sphx_glr_emd_tutorials_04_utilities_emd_tutorial_04_utils_02_logger.py: Using the logger ================ EMD has a built in logger which can be tuned to print out the progress of an analysis to the console, to a file or both. .. GENERATED FROM PYTHON SOURCE LINES 10-12 Logger Basics ^^^^^^^^^^^^^^ .. GENERATED FROM PYTHON SOURCE LINES 14-17 The logger must be initialised by calling ``emd.logger.set_up``. All subsequent calls to functions in the EMD library will then print output messages according to the logger specification. .. GENERATED FROM PYTHON SOURCE LINES 17-27 .. code-block:: Python # sphinx_gallery_thumbnail_path = '_static/emd_logger_thumb.png' # Import numpy for later import numpy as np # Import EMD and initialise the logger import emd emd.logger.set_up() .. rst-class:: sphx-glr-script-out .. code-block:: none EMD Logger Started EMD v installed in /home/docs/checkouts/readthedocs.org/user_builds/emd/checkouts/stable/emd .. GENERATED FROM PYTHON SOURCE LINES 28-38 The detail of logger output can be tuned by changing the logger level. The available levels are ``CRITICAL`` (only print output when the program is about to crash), ``WARNING`` (only print output when something unusual is happening or an analysis is potentially wrong), ``INFO`` (print general statements about which processes are running) and ``DEBUG`` (print loads of info including details of computations). The default level is ``INFO``, so if we re-initialise the logger on ``DEBUG`` we get more detailed outputs .. GENERATED FROM PYTHON SOURCE LINES 38-41 .. code-block:: Python emd.logger.set_up(level='DEBUG') .. rst-class:: sphx-glr-script-out .. code-block:: none EMD logger: handler 'console' level set to 'DEBUG' EMD Logger Started EMD v installed in /home/docs/checkouts/readthedocs.org/user_builds/emd/checkouts/stable/emd .. GENERATED FROM PYTHON SOURCE LINES 42-44 Lets explore the logger by running some a few sifts. Here we create a simple simulated oscillation and run a standard sift with default options. .. GENERATED FROM PYTHON SOURCE LINES 44-58 .. code-block:: Python # Initialise the logger with default settings (level=INFO) emd.logger.set_up(level='INFO') # Generate a simulation peak_freq = 12 sample_rate = 512 seconds = 10 noise_std = .5 x = emd.simulate.ar_oscillator(peak_freq, sample_rate, seconds, noise_std=noise_std, random_seed=42, r=.99) # Run a standard sift imf = emd.sift.sift(x) .. rst-class:: sphx-glr-script-out .. code-block:: none EMD logger: handler 'console' level set to 'INFO' EMD Logger Started STARTED: sift sifting IMF : 0 sifting IMF : 1 sifting IMF : 2 sifting IMF : 3 sifting IMF : 4 sifting IMF : 5 Finishing sift: reached energy ratio 56.149185420729 COMPLETED: sift .. GENERATED FROM PYTHON SOURCE LINES 59-64 With the level on ``INFO`` the logger tells us that the sift is running but not much else. If we change the logger level to ``DEBUG`` we get more output about how the sift is performing. The level of an initialised logger can be changed using ``emd.logger.set_level``. .. GENERATED FROM PYTHON SOURCE LINES 64-70 .. code-block:: Python emd.logger.set_level('DEBUG') # Run a standard sift imf = emd.sift.sift(x) .. rst-class:: sphx-glr-script-out .. code-block:: none EMD logger: handler 'console' level set to 'DEBUG' STARTED: sift Input data size: 5120 Input Sift Args: {} Checking sift inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) sifting IMF : 0 Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) Interpolating 'both' with method 'splrep' found 1644 minima and 1644 maxima on mode both Padding 2 extrema in signal X (5120,) using method 'rilling' L: max earlier than min, first val smaller than first min R: max earlier than min, last val larger than last max SD-thresh stop metric evaluated at iter 1 is : 5.35109985804001 Interpolating 'both' with method 'splrep' found 1743 minima and 1743 maxima on mode both Padding 2 extrema in signal X (5120,) using method 'rilling' L: max earlier than min, first val smaller than first min R: max earlier than min, last val smaller than last max Sift stopped by SD-thresh in 2 iters with sd 0.06412282310135851 sifting IMF : 1 Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) Interpolating 'both' with method 'splrep' found 528 minima and 529 maxima on mode both Padding 2 extrema in signal X (5120,) using method 'rilling' L: max earlier than min, first val larger than first min R: max later than min, last val larger than last min SD-thresh stop metric evaluated at iter 1 is : 2.4941806363513344 Interpolating 'both' with method 'splrep' found 650 minima and 651 maxima on mode both Padding 2 extrema in signal X (5120,) using method 'rilling' L: max earlier than min, first val larger than first min R: max later than min, last val larger than last min SD-thresh stop metric evaluated at iter 2 is : 0.8465302375518771 Interpolating 'both' with method 'splrep' found 680 minima and 681 maxima on mode both Padding 2 extrema in signal X (5120,) using method 'rilling' L: max earlier than min, first val larger than first min R: max later than min, last val larger than last min SD-thresh stop metric evaluated at iter 3 is : 0.526217083548086 Interpolating 'both' with method 'splrep' found 689 minima and 690 maxima on mode both Padding 2 extrema in signal X (5120,) using method 'rilling' L: max earlier than min, first val larger than first min R: max later than min, last val larger than last min SD-thresh stop metric evaluated at iter 4 is : 0.5224466697290732 Interpolating 'both' with method 'splrep' found 698 minima and 698 maxima on mode both Padding 2 extrema in signal X (5120,) using method 'rilling' L: max earlier than min, first val larger than first min R: max earlier than min, last val larger than last max SD-thresh stop metric evaluated at iter 5 is : 0.11871687916738188 Interpolating 'both' with method 'splrep' found 699 minima and 700 maxima on mode both Padding 2 extrema in signal X (5120,) using method 'rilling' L: max earlier than min, first val larger than first min R: max later than min, last val larger than last min Sift stopped by SD-thresh in 6 iters with sd 0.021359130568091945 sifting IMF : 2 Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) Interpolating 'both' with method 'splrep' found 169 minima and 168 maxima on mode both Padding 2 extrema in signal X (5120,) using method 'rilling' L: max later than min, first val smaller than first max R: max earlier than min, last val smaller than last max Sift stopped by SD-thresh in 1 iters with sd 0.08971365742623837 sifting IMF : 3 Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) Interpolating 'both' with method 'splrep' found 66 minima and 67 maxima on mode both Padding 2 extrema in signal X (5120,) using method 'rilling' L: max earlier than min, first val smaller than first min R: max later than min, last val larger than last min SD-thresh stop metric evaluated at iter 1 is : 0.33577493501437566 Interpolating 'both' with method 'splrep' found 67 minima and 68 maxima on mode both Padding 2 extrema in signal X (5120,) using method 'rilling' L: max earlier than min, first val smaller than first min R: max later than min, last val larger than last min Sift stopped by SD-thresh in 2 iters with sd 0.022552967284660053 sifting IMF : 4 Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) Interpolating 'both' with method 'splrep' found 27 minima and 28 maxima on mode both Padding 2 extrema in signal X (5120,) using method 'rilling' L: max earlier than min, first val smaller than first min R: max later than min, last val larger than last min SD-thresh stop metric evaluated at iter 1 is : 0.44204086673456716 Interpolating 'both' with method 'splrep' found 27 minima and 28 maxima on mode both Padding 2 extrema in signal X (5120,) using method 'rilling' L: max earlier than min, first val smaller than first min R: max later than min, last val larger than last min Sift stopped by SD-thresh in 2 iters with sd 0.016991665806757508 sifting IMF : 5 Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) Interpolating 'both' with method 'splrep' found 11 minima and 11 maxima on mode both Padding 2 extrema in signal X (5120,) using method 'rilling' L: max earlier than min, first val smaller than first min R: max earlier than min, last val larger than last max SD-thresh stop metric evaluated at iter 1 is : 0.12271568367694073 Interpolating 'both' with method 'splrep' found 11 minima and 11 maxima on mode both Padding 2 extrema in signal X (5120,) using method 'rilling' L: max earlier than min, first val smaller than first min R: max earlier than min, last val larger than last max Sift stopped by SD-thresh in 2 iters with sd 0.015311478985488196 Finishing sift: reached energy ratio 56.149185420729 Returning 7 imfs COMPLETED: sift .. GENERATED FROM PYTHON SOURCE LINES 71-77 If we don't want to change the logger level for our whole script, some functions allow us to override the logger level for a single call to the function. Any functions with the ``verbose`` option in 'Other Parameters' can override the logger level. Here we run a sift on ``WARNING`` and should see no outputs. .. GENERATED FROM PYTHON SOURCE LINES 77-80 .. code-block:: Python imf = emd.sift.sift(x, verbose='WARNING') .. GENERATED FROM PYTHON SOURCE LINES 81-85 We can also disable the logger altogether for a section of code using ``emd.logger.disable`` and then restart it using ``emd.logger.enable``. Here we disable logging for a sift call and renable it afterwards. .. GENERATED FROM PYTHON SOURCE LINES 85-90 .. code-block:: Python emd.logger.disable() imf = emd.sift.sift(x) emd.logger.enable() .. rst-class:: sphx-glr-script-out .. code-block:: none EMD logging disabled EMD logging enabled .. GENERATED FROM PYTHON SOURCE LINES 91-93 Advanced Logging ^^^^^^^^^^^^^^^^^ .. GENERATED FROM PYTHON SOURCE LINES 95-101 This section contains some logger optionality for more advanced use cases. Firstly, we can supplement the EMD logger output from a script by loading the EMD logger into a script and adding our own logger calls. In this example, we load the logger and add some custom updates. To print some output at the ``INFO`` level through the logger, we call ``logger.info``. .. GENERATED FROM PYTHON SOURCE LINES 101-124 .. code-block:: Python # Initialised EMD logger emd.logger.set_up(level='INFO') # Load logger into this script import logging logger = logging.getLogger('emd') # Check the time import time tic = time.perf_counter() # Run a sift logger.info('Now starting my new analysis') imf = emd.sift.sift(x) # Check time again toc = time.perf_counter() # Print sift run-time elapsed = toc - tic logger.info('My new analysis finished in {0:4f} seconds'.format(elapsed)) .. rst-class:: sphx-glr-script-out .. code-block:: none EMD logger: handler 'console' level set to 'INFO' EMD Logger Started Now starting my new analysis STARTED: sift sifting IMF : 0 sifting IMF : 1 sifting IMF : 2 sifting IMF : 3 sifting IMF : 4 sifting IMF : 5 Finishing sift: reached energy ratio 56.149185420729 COMPLETED: sift My new analysis finished in 0.033055 seconds .. GENERATED FROM PYTHON SOURCE LINES 125-131 This output respects the overall logger level, so info statements will be printed at levels ``INFO`` and ``DEBUG`` but suppressed if the overall logger level is at ``WARNING`` or ``CRITICAL``. Try changing the logger level in the example above to see the effect on the following output. .. GENERATED FROM PYTHON SOURCE LINES 133-139 Next, we define an analysis function which runs a sift followed by a frequency transform and simple cycle analysis. We'll run this example a few times in the next sections. Note that we've included some custom logger calls and a temporary logger override on the ``mask_sift`` .. GENERATED FROM PYTHON SOURCE LINES 139-178 .. code-block:: Python def my_analysis(x): # Check time, and load the logger into the session tic = time.perf_counter() import logging logger = logging.getLogger('emd') # Print start-up message logger.info('Starting my example analysis') logger.info('----------------------------') # Run a mask-sift with detailed logging imf = emd.sift.mask_sift(x, verbose='INFO') # Compute frequency stats IP, IF, IA = emd.spectra.frequency_transform(imf, sample_rate, 'nht') logger.info('Avg frequency of IMF-2 is {0:2f}Hz'.format(np.average(IF[:, 2], weights=IA[:, 2]))) # Find cycles in IMF-2 mask = IA[:, 2] > .05 cycles = emd.cycles.get_cycle_vector(IP, return_good=True, mask=mask) # Compute cycle stats cycle_freq = emd.cycles.get_cycle_stat(cycles[:, 2], IF[:, 2], func=np.mean)[1:] cycle_amp = emd.cycles.get_cycle_stat(cycles[:, 2], IA[:, 2], func=np.mean)[1:] # Print some cycle correlations logger.info('Freq-Amp correlation: r={0:2f}'.format(np.corrcoef(cycle_freq, cycle_amp)[0, 1])) # Print the elapsed time of the analysis toc = time.perf_counter() elapsed = toc - tic logger.info('My new analysis finished in {0:4f} seconds'.format(elapsed)) return cycle_freq, cycle_amp .. GENERATED FROM PYTHON SOURCE LINES 179-181 We can run this function as normal and inspect the logger outputs as they appear in the console. .. GENERATED FROM PYTHON SOURCE LINES 181-185 .. code-block:: Python # Run the analysis freq, amp = my_analysis(x) .. rst-class:: sphx-glr-script-out .. code-block:: none Starting my example analysis ---------------------------- EMD logger: handler 'console' level set to 'INFO' STARTED: mask_sift Computing first mask frequency with method zc Getting first IMF with no mask Found first mask frequency of 0.33955078125 Sifting IMF-0 Defining masks with freq 0.33955078125 and amp 8236.609552003989 at 4 phases Sifting IMF-1 Defining masks with freq 0.169775390625 and amp 8236.609552003989 at 4 phases Sifting IMF-2 Defining masks with freq 0.0848876953125 and amp 8236.609552003989 at 4 phases Sifting IMF-3 Defining masks with freq 0.04244384765625 and amp 8236.609552003989 at 4 phases Sifting IMF-4 Defining masks with freq 0.021221923828125 and amp 8236.609552003989 at 4 phases Sifting IMF-5 Defining masks with freq 0.0106109619140625 and amp 8236.609552003989 at 4 phases Sifting IMF-6 Defining masks with freq 0.00530548095703125 and amp 8236.609552003989 at 4 phases Sifting IMF-7 Defining masks with freq 0.002652740478515625 and amp 8236.609552003989 at 4 phases Sifting IMF-8 Defining masks with freq 0.0013263702392578125 and amp 8236.609552003989 at 4 phases Finishing sift: reached max number of imfs (9) COMPLETED: mask_sift EMD logger: handler 'console' level set to 'INFO' STARTED: compute frequency stats Using Amplitude-Normalised Hilbert transform STARTED: Amplitude-Normalise Normalise of IMF-0-0 complete in 3 iters (val=0.0) Normalise of IMF-1-0 complete in 3 iters (val=0.0) Normalise of IMF-2-0 complete in 3 iters (val=0.0) Normalise of IMF-3-0 complete in 2 iters (val=0.0) Normalise of IMF-4-0 complete in 3 iters (val=0.0) Normalise of IMF-5-0 complete in 3 iters (val=0.0) Normalise of IMF-6-0 complete in 3 iters (val=0.0) Normalise of IMF-7-0 complete in 3 iters (val=0.0) Normalise of IMF-8-0 complete in 2 iters (val=0.0) COMPLETED: Amplitude-Normalise COMPLETED: compute frequency stats. Returning 9 imfs Avg frequency of IMF-2 is 36.310268Hz STARTED: get cycle indices found 1 cycles in IMF-0 found 53 cycles in IMF-1 found 105 cycles in IMF-2 found 100 cycles in IMF-3 found 92 cycles in IMF-4 found 27 cycles in IMF-5 found 7 cycles in IMF-6 found 1 cycles in IMF-7 found 0 cycles in IMF-8 COMPLETED: get cycle indices STARTED: get_cycle_stat STARTED: get_cycle_stat Freq-Amp correlation: r=-0.366304 My new analysis finished in 0.205080 seconds .. GENERATED FROM PYTHON SOURCE LINES 186-195 We can see a lot of information about which processes were running as the logger is set to ``INFO``. If we configure the logger to run on ``WARNING`` level, we should only see output about potential errors. This can be useful when running familiar code where you only need output when something potentially strange is going on. So in this case, after changing to ``WARNING`` we should only see output from the ``mask_sift`` call (as this has a logger override to ``DEBUG``). All other output is suppressed. .. GENERATED FROM PYTHON SOURCE LINES 195-202 .. code-block:: Python # Change logger level emd.logger.set_level('WARNING') # Run the analysis freq, amp = my_analysis(x) .. rst-class:: sphx-glr-script-out .. code-block:: none STARTED: mask_sift Computing first mask frequency with method zc Getting first IMF with no mask Found first mask frequency of 0.33955078125 Sifting IMF-0 Defining masks with freq 0.33955078125 and amp 8236.609552003989 at 4 phases Sifting IMF-1 Defining masks with freq 0.169775390625 and amp 8236.609552003989 at 4 phases Sifting IMF-2 Defining masks with freq 0.0848876953125 and amp 8236.609552003989 at 4 phases Sifting IMF-3 Defining masks with freq 0.04244384765625 and amp 8236.609552003989 at 4 phases Sifting IMF-4 Defining masks with freq 0.021221923828125 and amp 8236.609552003989 at 4 phases Sifting IMF-5 Defining masks with freq 0.0106109619140625 and amp 8236.609552003989 at 4 phases Sifting IMF-6 Defining masks with freq 0.00530548095703125 and amp 8236.609552003989 at 4 phases Sifting IMF-7 Defining masks with freq 0.002652740478515625 and amp 8236.609552003989 at 4 phases Sifting IMF-8 Defining masks with freq 0.0013263702392578125 and amp 8236.609552003989 at 4 phases Finishing sift: reached max number of imfs (9) COMPLETED: mask_sift .. GENERATED FROM PYTHON SOURCE LINES 203-210 If we're dealing with logger output from multiple sources or perhaps from multiple analyses running in parallel, then we can add a prefix to the logger to help distinguish the output coming from each. This can be specified when initialising the logger. For example, here we reun the analysis function with a logger prefix indicating that we're processing Subject 001. .. GENERATED FROM PYTHON SOURCE LINES 210-217 .. code-block:: Python # Initialise logger with a prefix emd.logger.set_up(level='INFO', prefix='Subj001') # Run the analysis freq, amp = my_analysis(x) .. rst-class:: sphx-glr-script-out .. code-block:: none Subj001 : EMD logger: handler 'console' level set to 'INFO' Subj001 : EMD Logger Started Subj001 : Starting my example analysis Subj001 : ---------------------------- Subj001 : EMD logger: handler 'console' level set to 'INFO' Subj001 : STARTED: mask_sift Subj001 : Computing first mask frequency with method zc Subj001 : Getting first IMF with no mask Subj001 : Found first mask frequency of 0.33955078125 Subj001 : Sifting IMF-0 Subj001 : Defining masks with freq 0.33955078125 and amp 8236.609552003989 at 4 phases Subj001 : Sifting IMF-1 Subj001 : Defining masks with freq 0.169775390625 and amp 8236.609552003989 at 4 phases Subj001 : Sifting IMF-2 Subj001 : Defining masks with freq 0.0848876953125 and amp 8236.609552003989 at 4 phases Subj001 : Sifting IMF-3 Subj001 : Defining masks with freq 0.04244384765625 and amp 8236.609552003989 at 4 phases Subj001 : Sifting IMF-4 Subj001 : Defining masks with freq 0.021221923828125 and amp 8236.609552003989 at 4 phases Subj001 : Sifting IMF-5 Subj001 : Defining masks with freq 0.0106109619140625 and amp 8236.609552003989 at 4 phases Subj001 : Sifting IMF-6 Subj001 : Defining masks with freq 0.00530548095703125 and amp 8236.609552003989 at 4 phases Subj001 : Sifting IMF-7 Subj001 : Defining masks with freq 0.002652740478515625 and amp 8236.609552003989 at 4 phases Subj001 : Sifting IMF-8 Subj001 : Defining masks with freq 0.0013263702392578125 and amp 8236.609552003989 at 4 phases Subj001 : Finishing sift: reached max number of imfs (9) Subj001 : COMPLETED: mask_sift Subj001 : EMD logger: handler 'console' level set to 'INFO' Subj001 : STARTED: compute frequency stats Subj001 : Using Amplitude-Normalised Hilbert transform Subj001 : STARTED: Amplitude-Normalise Subj001 : Normalise of IMF-0-0 complete in 3 iters (val=0.0) Subj001 : Normalise of IMF-1-0 complete in 3 iters (val=0.0) Subj001 : Normalise of IMF-2-0 complete in 3 iters (val=0.0) Subj001 : Normalise of IMF-3-0 complete in 2 iters (val=0.0) Subj001 : Normalise of IMF-4-0 complete in 3 iters (val=0.0) Subj001 : Normalise of IMF-5-0 complete in 3 iters (val=0.0) Subj001 : Normalise of IMF-6-0 complete in 3 iters (val=0.0) Subj001 : Normalise of IMF-7-0 complete in 3 iters (val=0.0) Subj001 : Normalise of IMF-8-0 complete in 2 iters (val=0.0) Subj001 : COMPLETED: Amplitude-Normalise Subj001 : COMPLETED: compute frequency stats. Returning 9 imfs Subj001 : Avg frequency of IMF-2 is 36.310268Hz Subj001 : STARTED: get cycle indices Subj001 : found 1 cycles in IMF-0 Subj001 : found 53 cycles in IMF-1 Subj001 : found 105 cycles in IMF-2 Subj001 : found 100 cycles in IMF-3 Subj001 : found 92 cycles in IMF-4 Subj001 : found 27 cycles in IMF-5 Subj001 : found 7 cycles in IMF-6 Subj001 : found 1 cycles in IMF-7 Subj001 : found 0 cycles in IMF-8 Subj001 : COMPLETED: get cycle indices Subj001 : STARTED: get_cycle_stat Subj001 : STARTED: get_cycle_stat Subj001 : Freq-Amp correlation: r=-0.366304 Subj001 : My new analysis finished in 0.206710 seconds .. GENERATED FROM PYTHON SOURCE LINES 218-219 Finally, we can direct the logger output into a text file as well as the console. .. GENERATED FROM PYTHON SOURCE LINES 219-232 .. code-block:: Python # Define a temporary file import tempfile log_file = tempfile.NamedTemporaryFile(prefix="ExampleEMDLog", suffix='.log').name # OR uncomment this line and define your own filepath # log_file = /path/to/my/log_file.log # Initialise the logger with a prefix and a text file emd.logger.set_up(level='INFO', prefix='Subj001', log_file=log_file) # Run the analysis freq, amp = my_analysis(x) .. rst-class:: sphx-glr-script-out .. code-block:: none Subj001 : EMD logger: handler 'console' level set to 'INFO' Subj001 : EMD Logger Started Subj001 : logging to file: /tmp/ExampleEMDLogbfmc_xf7.log Subj001 : Starting my example analysis Subj001 : ---------------------------- Subj001 : EMD logger: handler 'console' level set to 'INFO' Subj001 : STARTED: mask_sift Subj001 : Computing first mask frequency with method zc Subj001 : Getting first IMF with no mask Subj001 : Found first mask frequency of 0.33955078125 Subj001 : Sifting IMF-0 Subj001 : Defining masks with freq 0.33955078125 and amp 8236.609552003989 at 4 phases Subj001 : Sifting IMF-1 Subj001 : Defining masks with freq 0.169775390625 and amp 8236.609552003989 at 4 phases Subj001 : Sifting IMF-2 Subj001 : Defining masks with freq 0.0848876953125 and amp 8236.609552003989 at 4 phases Subj001 : Sifting IMF-3 Subj001 : Defining masks with freq 0.04244384765625 and amp 8236.609552003989 at 4 phases Subj001 : Sifting IMF-4 Subj001 : Defining masks with freq 0.021221923828125 and amp 8236.609552003989 at 4 phases Subj001 : Sifting IMF-5 Subj001 : Defining masks with freq 0.0106109619140625 and amp 8236.609552003989 at 4 phases Subj001 : Sifting IMF-6 Subj001 : Defining masks with freq 0.00530548095703125 and amp 8236.609552003989 at 4 phases Subj001 : Sifting IMF-7 Subj001 : Defining masks with freq 0.002652740478515625 and amp 8236.609552003989 at 4 phases Subj001 : Sifting IMF-8 Subj001 : Defining masks with freq 0.0013263702392578125 and amp 8236.609552003989 at 4 phases Subj001 : Finishing sift: reached max number of imfs (9) Subj001 : COMPLETED: mask_sift Subj001 : EMD logger: handler 'console' level set to 'INFO' Subj001 : STARTED: compute frequency stats Subj001 : Using Amplitude-Normalised Hilbert transform Subj001 : STARTED: Amplitude-Normalise Subj001 : Normalise of IMF-0-0 complete in 3 iters (val=0.0) Subj001 : Normalise of IMF-1-0 complete in 3 iters (val=0.0) Subj001 : Normalise of IMF-2-0 complete in 3 iters (val=0.0) Subj001 : Normalise of IMF-3-0 complete in 2 iters (val=0.0) Subj001 : Normalise of IMF-4-0 complete in 3 iters (val=0.0) Subj001 : Normalise of IMF-5-0 complete in 3 iters (val=0.0) Subj001 : Normalise of IMF-6-0 complete in 3 iters (val=0.0) Subj001 : Normalise of IMF-7-0 complete in 3 iters (val=0.0) Subj001 : Normalise of IMF-8-0 complete in 2 iters (val=0.0) Subj001 : COMPLETED: Amplitude-Normalise Subj001 : COMPLETED: compute frequency stats. Returning 9 imfs Subj001 : Avg frequency of IMF-2 is 36.310268Hz Subj001 : STARTED: get cycle indices Subj001 : found 1 cycles in IMF-0 Subj001 : found 53 cycles in IMF-1 Subj001 : found 105 cycles in IMF-2 Subj001 : found 100 cycles in IMF-3 Subj001 : found 92 cycles in IMF-4 Subj001 : found 27 cycles in IMF-5 Subj001 : found 7 cycles in IMF-6 Subj001 : found 1 cycles in IMF-7 Subj001 : found 0 cycles in IMF-8 Subj001 : COMPLETED: get cycle indices Subj001 : STARTED: get_cycle_stat Subj001 : STARTED: get_cycle_stat Subj001 : Freq-Amp correlation: r=-0.366304 Subj001 : My new analysis finished in 0.241116 seconds .. GENERATED FROM PYTHON SOURCE LINES 233-238 The log file is a simple text file containing very detailed outputs of which functions were executed and when. Here we read the log file and print its contents to the console. Note that the log file contains a much more detailed output that the console! .. GENERATED FROM PYTHON SOURCE LINES 238-244 .. code-block:: Python # Open the text file and print its contents with open(log_file, 'r') as f: txt = f.read() print(txt) .. rst-class:: sphx-glr-script-out .. code-block:: none [2024-03-14 14:13:09] Subj001 : - INFO - emd.logger:195 - set_level() : EMD logger: handler 'console' level set to 'INFO' [2024-03-14 14:13:09] Subj001 : - INFO - emd.logger:180 - set_up() : EMD Logger Started [2024-03-14 14:13:09] Subj001 : - INFO - emd.logger:184 - set_up() : logging to file: /tmp/ExampleEMDLogbfmc_xf7.log [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : EMD v installed in /home/docs/checkouts/readthedocs.org/user_builds/emd/checkouts/stable/emd [2024-03-14 14:13:09] Subj001 : - INFO - emd.emd_tutorial_04_utils_02_logger:149 - my_analysis() : Starting my example analysis [2024-03-14 14:13:09] Subj001 : - INFO - emd.emd_tutorial_04_utils_02_logger:150 - my_analysis() : ---------------------------- [2024-03-14 14:13:09] Subj001 : - INFO - emd.logger:195 - set_level() : EMD logger: handler 'console' level set to 'INFO' [2024-03-14 14:13:09] Subj001 : - INFO - emd.logger:261 - sift_logger() : STARTED: mask_sift [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.logger:271 - sift_logger() : Input data size: 5120 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.logger:274 - sift_logger() : Input Sift Args: {'verbose': 'INFO'} [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking sift inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - INFO - emd.sift:1084 - get_mask_freqs() : Computing first mask frequency with method zc [2024-03-14 14:13:09] Subj001 : - INFO - emd.sift:1085 - get_mask_freqs() : Getting first IMF with no mask [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 1644 minima and 1644 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:256 - _pad_extrema_rilling() : L: max earlier than min, first val smaller than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:287 - _pad_extrema_rilling() : R: max earlier than min, last val larger than last max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.sift:287 - stop_imf_sd() : SD-thresh stop metric evaluated at iter 1 is : 5.35109985804001 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 1743 minima and 1743 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:256 - _pad_extrema_rilling() : L: max earlier than min, first val smaller than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:281 - _pad_extrema_rilling() : R: max earlier than min, last val smaller than last max [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 2 iters with sd 0.06412282310135851 [2024-03-14 14:13:09] Subj001 : - INFO - emd.sift:1093 - get_mask_freqs() : Found first mask frequency of 0.33955078125 [2024-03-14 14:13:09] Subj001 : - INFO - emd.sift:1259 - mask_sift() : Sifting IMF-0 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf_mask inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - INFO - emd.sift:1030 - get_next_imf_mask() : Defining masks with freq 0.33955078125 and amp 8236.609552003989 at 4 phases [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 1735 minima and 1734 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:271 - _pad_extrema_rilling() : L: max later than min, first val smaller than first max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:281 - _pad_extrema_rilling() : R: max earlier than min, last val smaller than last max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.sift:287 - stop_imf_sd() : SD-thresh stop metric evaluated at iter 1 is : 1.3789635105451155 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 1737 minima and 1736 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:271 - _pad_extrema_rilling() : L: max later than min, first val smaller than first max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:287 - _pad_extrema_rilling() : R: max earlier than min, last val larger than last max [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 2 iters with sd 0.004623719861115919 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 1735 minima and 1735 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:271 - _pad_extrema_rilling() : L: max later than min, first val smaller than first max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.sift:287 - stop_imf_sd() : SD-thresh stop metric evaluated at iter 1 is : 1.3147285371322173 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 1739 minima and 1739 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:271 - _pad_extrema_rilling() : L: max later than min, first val smaller than first max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 2 iters with sd 0.004133391824521104 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 1737 minima and 1737 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:256 - _pad_extrema_rilling() : L: max earlier than min, first val smaller than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:281 - _pad_extrema_rilling() : R: max earlier than min, last val smaller than last max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.sift:287 - stop_imf_sd() : SD-thresh stop metric evaluated at iter 1 is : 1.3114641484437155 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 1737 minima and 1737 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:281 - _pad_extrema_rilling() : R: max earlier than min, last val smaller than last max [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 2 iters with sd 0.0039783220679981525 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 1735 minima and 1735 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:287 - _pad_extrema_rilling() : R: max earlier than min, last val larger than last max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.sift:287 - stop_imf_sd() : SD-thresh stop metric evaluated at iter 1 is : 1.3704622921607965 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 1738 minima and 1738 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:287 - _pad_extrema_rilling() : R: max earlier than min, last val larger than last max [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 2 iters with sd 0.0043234337469501985 [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Averaging across 4 proto IMFs [2024-03-14 14:13:09] Subj001 : - INFO - emd.sift:1259 - mask_sift() : Sifting IMF-1 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf_mask inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - INFO - emd.sift:1030 - get_next_imf_mask() : Defining masks with freq 0.169775390625 and amp 8236.609552003989 at 4 phases [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 869 minima and 869 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:265 - _pad_extrema_rilling() : L: max later than min, first val larger than first max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.sift:287 - stop_imf_sd() : SD-thresh stop metric evaluated at iter 1 is : 1.5390230090961547 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 869 minima and 869 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:265 - _pad_extrema_rilling() : L: max later than min, first val larger than first max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 2 iters with sd 0.002497476774192765 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 869 minima and 869 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:271 - _pad_extrema_rilling() : L: max later than min, first val smaller than first max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.sift:287 - stop_imf_sd() : SD-thresh stop metric evaluated at iter 1 is : 1.5756192467542398 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 869 minima and 869 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:271 - _pad_extrema_rilling() : L: max later than min, first val smaller than first max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 2 iters with sd 0.0026827942833537004 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 868 minima and 869 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.sift:287 - stop_imf_sd() : SD-thresh stop metric evaluated at iter 1 is : 1.560900501833889 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 868 minima and 869 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 2 iters with sd 0.0023174964672515817 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 869 minima and 869 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:281 - _pad_extrema_rilling() : R: max earlier than min, last val smaller than last max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.sift:287 - stop_imf_sd() : SD-thresh stop metric evaluated at iter 1 is : 1.529833317604314 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 869 minima and 869 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:281 - _pad_extrema_rilling() : R: max earlier than min, last val smaller than last max [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 2 iters with sd 0.0024305461759021965 [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Averaging across 4 proto IMFs [2024-03-14 14:13:09] Subj001 : - INFO - emd.sift:1259 - mask_sift() : Sifting IMF-2 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf_mask inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - INFO - emd.sift:1030 - get_next_imf_mask() : Defining masks with freq 0.0848876953125 and amp 8236.609552003989 at 4 phases [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 434 minima and 434 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:265 - _pad_extrema_rilling() : L: max later than min, first val larger than first max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.sift:287 - stop_imf_sd() : SD-thresh stop metric evaluated at iter 1 is : 1.531132019094807 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 434 minima and 434 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:265 - _pad_extrema_rilling() : L: max later than min, first val larger than first max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 2 iters with sd 0.0030009458306292235 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 435 minima and 434 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:271 - _pad_extrema_rilling() : L: max later than min, first val smaller than first max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:281 - _pad_extrema_rilling() : R: max earlier than min, last val smaller than last max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.sift:287 - stop_imf_sd() : SD-thresh stop metric evaluated at iter 1 is : 1.5103862889963942 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 435 minima and 434 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:271 - _pad_extrema_rilling() : L: max later than min, first val smaller than first max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:281 - _pad_extrema_rilling() : R: max earlier than min, last val smaller than last max [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 2 iters with sd 0.002603505258214974 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 434 minima and 435 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.sift:287 - stop_imf_sd() : SD-thresh stop metric evaluated at iter 1 is : 1.4665697234971053 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 434 minima and 435 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 2 iters with sd 0.0028835702286500705 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 434 minima and 435 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.sift:287 - stop_imf_sd() : SD-thresh stop metric evaluated at iter 1 is : 1.4868036775239224 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 434 minima and 435 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 2 iters with sd 0.0029152251830059255 [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Averaging across 4 proto IMFs [2024-03-14 14:13:09] Subj001 : - INFO - emd.sift:1259 - mask_sift() : Sifting IMF-3 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf_mask inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - INFO - emd.sift:1030 - get_next_imf_mask() : Defining masks with freq 0.04244384765625 and amp 8236.609552003989 at 4 phases [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 210 minima and 211 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.sift:287 - stop_imf_sd() : SD-thresh stop metric evaluated at iter 1 is : 0.13110876324665582 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 210 minima and 211 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 2 iters with sd 0.028557455670537003 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 206 minima and 206 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:271 - _pad_extrema_rilling() : L: max later than min, first val smaller than first max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.sift:287 - stop_imf_sd() : SD-thresh stop metric evaluated at iter 1 is : 0.12747281555032836 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 206 minima and 206 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:265 - _pad_extrema_rilling() : L: max later than min, first val larger than first max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 2 iters with sd 0.02367094727973033 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 207 minima and 207 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:256 - _pad_extrema_rilling() : L: max earlier than min, first val smaller than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:281 - _pad_extrema_rilling() : R: max earlier than min, last val smaller than last max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.sift:287 - stop_imf_sd() : SD-thresh stop metric evaluated at iter 1 is : 0.11731517843127388 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 207 minima and 207 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:281 - _pad_extrema_rilling() : R: max earlier than min, last val smaller than last max [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 2 iters with sd 0.026190556382307998 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 211 minima and 212 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.sift:287 - stop_imf_sd() : SD-thresh stop metric evaluated at iter 1 is : 0.11469751415340777 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 211 minima and 212 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 2 iters with sd 0.032779118918049986 [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Averaging across 4 proto IMFs [2024-03-14 14:13:09] Subj001 : - INFO - emd.sift:1259 - mask_sift() : Sifting IMF-4 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf_mask inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - INFO - emd.sift:1030 - get_next_imf_mask() : Defining masks with freq 0.021221923828125 and amp 8236.609552003989 at 4 phases [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 109 minima and 108 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:271 - _pad_extrema_rilling() : L: max later than min, first val smaller than first max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:281 - _pad_extrema_rilling() : R: max earlier than min, last val smaller than last max [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 1 iters with sd 0.014292054510223652 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 109 minima and 108 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:271 - _pad_extrema_rilling() : L: max later than min, first val smaller than first max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:281 - _pad_extrema_rilling() : R: max earlier than min, last val smaller than last max [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 1 iters with sd 0.02269985577307988 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 109 minima and 109 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:271 - _pad_extrema_rilling() : L: max later than min, first val smaller than first max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 1 iters with sd 0.031891690424448234 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 108 minima and 109 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:256 - _pad_extrema_rilling() : L: max earlier than min, first val smaller than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:301 - _pad_extrema_rilling() : R: max later than min, last val smaller than last min [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 1 iters with sd 0.014882633541519674 [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Averaging across 4 proto IMFs [2024-03-14 14:13:09] Subj001 : - INFO - emd.sift:1259 - mask_sift() : Sifting IMF-5 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf_mask inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - INFO - emd.sift:1030 - get_next_imf_mask() : Defining masks with freq 0.0106109619140625 and amp 8236.609552003989 at 4 phases [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 54 minima and 55 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 1 iters with sd 0.012477217831270666 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 55 minima and 54 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:271 - _pad_extrema_rilling() : L: max later than min, first val smaller than first max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:281 - _pad_extrema_rilling() : R: max earlier than min, last val smaller than last max [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 1 iters with sd 0.03173757419559748 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 54 minima and 54 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:256 - _pad_extrema_rilling() : L: max earlier than min, first val smaller than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:281 - _pad_extrema_rilling() : R: max earlier than min, last val smaller than last max [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 1 iters with sd 0.011561903979950975 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 54 minima and 55 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:256 - _pad_extrema_rilling() : L: max earlier than min, first val smaller than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 1 iters with sd 0.008042658099794611 [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Averaging across 4 proto IMFs [2024-03-14 14:13:09] Subj001 : - INFO - emd.sift:1259 - mask_sift() : Sifting IMF-6 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf_mask inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - INFO - emd.sift:1030 - get_next_imf_mask() : Defining masks with freq 0.00530548095703125 and amp 8236.609552003989 at 4 phases [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 27 minima and 28 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 1 iters with sd 0.005571997537588026 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 27 minima and 27 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:271 - _pad_extrema_rilling() : L: max later than min, first val smaller than first max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 1 iters with sd 0.027507814058890287 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 27 minima and 27 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:256 - _pad_extrema_rilling() : L: max earlier than min, first val smaller than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:281 - _pad_extrema_rilling() : R: max earlier than min, last val smaller than last max [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 1 iters with sd 0.004292108979085951 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 27 minima and 27 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:281 - _pad_extrema_rilling() : R: max earlier than min, last val smaller than last max [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 1 iters with sd 0.0021921019611720553 [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Averaging across 4 proto IMFs [2024-03-14 14:13:09] Subj001 : - INFO - emd.sift:1259 - mask_sift() : Sifting IMF-7 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf_mask inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - INFO - emd.sift:1030 - get_next_imf_mask() : Defining masks with freq 0.002652740478515625 and amp 8236.609552003989 at 4 phases [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 14 minima and 14 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:281 - _pad_extrema_rilling() : R: max earlier than min, last val smaller than last max [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 1 iters with sd 0.002349210487310889 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 14 minima and 13 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:271 - _pad_extrema_rilling() : L: max later than min, first val smaller than first max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:281 - _pad_extrema_rilling() : R: max earlier than min, last val smaller than last max [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 1 iters with sd 0.028984048367158653 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 13 minima and 14 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:256 - _pad_extrema_rilling() : L: max earlier than min, first val smaller than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 1 iters with sd 0.0017924683005907703 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 13 minima and 14 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 1 iters with sd 0.0007426849451273461 [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Averaging across 4 proto IMFs [2024-03-14 14:13:09] Subj001 : - INFO - emd.sift:1259 - mask_sift() : Sifting IMF-8 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf_mask inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - INFO - emd.sift:1030 - get_next_imf_mask() : Defining masks with freq 0.0013263702392578125 and amp 8236.609552003989 at 4 phases [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 7 minima and 7 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:281 - _pad_extrema_rilling() : R: max earlier than min, last val smaller than last max [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 1 iters with sd 0.0030141254472833526 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 7 minima and 7 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:271 - _pad_extrema_rilling() : L: max later than min, first val smaller than first max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 1 iters with sd 0.04499571071584374 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 6 minima and 7 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:256 - _pad_extrema_rilling() : L: max earlier than min, first val smaller than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 1 iters with sd 0.0023666068736690115 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:217 - ensure_1d_with_singleton() : Checking get_next_imf inputs - Trimming trailing singletons from input 'X' (input size (5120, 1)) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'both' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 7 minima and 7 maxima on mode both [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:281 - _pad_extrema_rilling() : R: max earlier than min, last val smaller than last max [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Sift stopped by SD-thresh in 1 iters with sd 0.0009500624908985781 [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Averaging across 4 proto IMFs [2024-03-14 14:13:09] Subj001 : - INFO - emd.sift:1276 - mask_sift() : Finishing sift: reached max number of imfs (9) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.logger:282 - sift_logger() : Returning 9 imfs [2024-03-14 14:13:09] Subj001 : - INFO - emd.logger:287 - sift_logger() : COMPLETED: mask_sift [2024-03-14 14:13:09] Subj001 : - INFO - emd.logger:195 - set_level() : EMD logger: handler 'console' level set to 'INFO' [2024-03-14 14:13:09] Subj001 : - INFO - emd.spectra:84 - frequency_transform() : STARTED: compute frequency stats [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.spectra:87 - frequency_transform() : computing on 5120 samples over 9 imfs at sample rate 512 [2024-03-14 14:13:09] Subj001 : - INFO - emd.spectra:102 - frequency_transform() : Using Amplitude-Normalised Hilbert transform [2024-03-14 14:13:09] Subj001 : - INFO - emd.imftools:58 - amplitude_normalise() : STARTED: Amplitude-Normalise [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.imftools:61 - amplitude_normalise() : Normalising 5120 samples across 9 IMFs [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.imftools:64 - amplitude_normalise() : Using pchip interpolation with threshold of 1e-10 and max_iters 3 [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 1611 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 1652 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 1653 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 1653 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - INFO - emd.imftools:106 - amplitude_normalise() : Normalise of IMF-0-0 complete in 3 iters (val=0.0) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 1437 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 1442 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 1442 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 1442 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - INFO - emd.imftools:106 - amplitude_normalise() : Normalise of IMF-1-0 complete in 3 iters (val=0.0) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 763 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 764 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 764 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 764 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - INFO - emd.imftools:106 - amplitude_normalise() : Normalise of IMF-2-0 complete in 3 iters (val=0.0) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 260 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 260 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 260 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - INFO - emd.imftools:106 - amplitude_normalise() : Normalise of IMF-3-0 complete in 2 iters (val=0.0) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 228 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 228 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 228 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 228 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - INFO - emd.imftools:106 - amplitude_normalise() : Normalise of IMF-4-0 complete in 3 iters (val=0.0) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 87 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 87 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 87 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 87 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - INFO - emd.imftools:106 - amplitude_normalise() : Normalise of IMF-5-0 complete in 3 iters (val=0.0) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 44 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 44 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 44 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 44 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - INFO - emd.imftools:106 - amplitude_normalise() : Normalise of IMF-6-0 complete in 3 iters (val=0.0) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 22 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 22 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 22 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 22 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - INFO - emd.imftools:106 - amplitude_normalise() : Normalise of IMF-7-0 complete in 3 iters (val=0.0) [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 10 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 10 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 10 extrema on mode abs_peaks [2024-03-14 14:13:09] Subj001 : - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:09] Subj001 : - INFO - emd.imftools:106 - amplitude_normalise() : Normalise of IMF-8-0 complete in 2 iters (val=0.0) [2024-03-14 14:13:09] Subj001 : - INFO - emd.imftools:119 - amplitude_normalise() : COMPLETED: Amplitude-Normalise [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'upper' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 1826 minima and 1826 maxima on mode peaks [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:287 - _pad_extrema_rilling() : R: max earlier than min, last val larger than last max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'upper' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 777 minima and 777 maxima on mode peaks [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:271 - _pad_extrema_rilling() : L: max later than min, first val smaller than first max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'upper' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 395 minima and 396 maxima on mode peaks [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'upper' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 132 minima and 133 maxima on mode peaks [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'upper' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 116 minima and 115 maxima on mode peaks [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:271 - _pad_extrema_rilling() : L: max later than min, first val smaller than first max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:281 - _pad_extrema_rilling() : R: max earlier than min, last val smaller than last max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'upper' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 45 minima and 45 maxima on mode peaks [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:256 - _pad_extrema_rilling() : L: max earlier than min, first val smaller than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:281 - _pad_extrema_rilling() : R: max earlier than min, last val smaller than last max [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'upper' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 23 minima and 24 maxima on mode peaks [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:256 - _pad_extrema_rilling() : L: max earlier than min, first val smaller than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'upper' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 11 minima and 12 maxima on mode peaks [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:256 - _pad_extrema_rilling() : L: max earlier than min, first val smaller than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'upper' with method 'splrep' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 6 minima and 6 maxima on mode peaks [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:09] Subj001 : - DEBUG - emd._sift_core:281 - _pad_extrema_rilling() : R: max earlier than min, last val smaller than last max [2024-03-14 14:13:09] Subj001 : - INFO - emd.spectra:193 - frequency_transform() : COMPLETED: compute frequency stats. Returning 9 imfs [2024-03-14 14:13:09] Subj001 : - INFO - emd.emd_tutorial_04_utils_02_logger:157 - my_analysis() : Avg frequency of IMF-2 is 36.310268Hz [2024-03-14 14:13:09] Subj001 : - INFO - emd.cycles:131 - get_cycle_vector() : STARTED: get cycle indices [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.support:261 - ensure_2d() : Checking get_cycle_vector inputs - Adding dummy dimension to input 'mask' [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.cycles:138 - get_cycle_vector() : computing on 5120 samples over 9 IMFs [2024-03-14 14:13:09] Subj001 : - DEBUG - emd.cycles:141 - get_cycle_vector() : 5027 (98.18%) samples masked out [2024-03-14 14:13:09] Subj001 : - INFO - emd.cycles:191 - get_cycle_vector() : found 1 cycles in IMF-0 [2024-03-14 14:13:09] Subj001 : - INFO - emd.cycles:191 - get_cycle_vector() : found 53 cycles in IMF-1 [2024-03-14 14:13:09] Subj001 : - INFO - emd.cycles:191 - get_cycle_vector() : found 105 cycles in IMF-2 [2024-03-14 14:13:09] Subj001 : - INFO - emd.cycles:191 - get_cycle_vector() : found 100 cycles in IMF-3 [2024-03-14 14:13:09] Subj001 : - INFO - emd.cycles:191 - get_cycle_vector() : found 92 cycles in IMF-4 [2024-03-14 14:13:09] Subj001 : - INFO - emd.cycles:191 - get_cycle_vector() : found 27 cycles in IMF-5 [2024-03-14 14:13:09] Subj001 : - INFO - emd.cycles:191 - get_cycle_vector() : found 7 cycles in IMF-6 [2024-03-14 14:13:09] Subj001 : - INFO - emd.cycles:191 - get_cycle_vector() : found 1 cycles in IMF-7 [2024-03-14 14:13:09] Subj001 : - INFO - emd.cycles:191 - get_cycle_vector() : found 0 cycles in IMF-8 [2024-03-14 14:13:09] Subj001 : - INFO - emd.cycles:193 - get_cycle_vector() : COMPLETED: get cycle indices [2024-03-14 14:13:09] Subj001 : - INFO - emd.cycles:423 - get_cycle_stat() : STARTED: get_cycle_stat [2024-03-14 14:13:09] Subj001 : - INFO - emd.cycles:423 - get_cycle_stat() : STARTED: get_cycle_stat [2024-03-14 14:13:09] Subj001 : - INFO - emd.emd_tutorial_04_utils_02_logger:168 - my_analysis() : Freq-Amp correlation: r=-0.366304 [2024-03-14 14:13:09] Subj001 : - INFO - emd.emd_tutorial_04_utils_02_logger:173 - my_analysis() : My new analysis finished in 0.241116 seconds .. GENERATED FROM PYTHON SOURCE LINES 245-247 If we want this extra detailed output in the console as well, we can specify the ``console_format`` when setting up the logger. .. GENERATED FROM PYTHON SOURCE LINES 247-253 .. code-block:: Python # Initialise logger with a verbose console format emd.logger.set_up(level='DEBUG', prefix='Subj001', console_format='verbose') # Run the analysis freq, amp = my_analysis(x) .. rst-class:: sphx-glr-script-out .. code-block:: none Subj001 EMD logger: handler 'console' level set to 'DEBUG' [2024-03-14 14:13:09,948] Subj001 - INFO - emd.logger:220 - set_format() : EMD logger: handler console format changed to verbose [2024-03-14 14:13:09,948] Subj001 - INFO - emd.logger:180 - set_up() : EMD Logger Started [2024-03-14 14:13:09,948] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : EMD v installed in /home/docs/checkouts/readthedocs.org/user_builds/emd/checkouts/stable/emd [2024-03-14 14:13:09,948] Subj001 - INFO - emd.emd_tutorial_04_utils_02_logger:149 - my_analysis() : Starting my example analysis [2024-03-14 14:13:09,948] Subj001 - INFO - emd.emd_tutorial_04_utils_02_logger:150 - my_analysis() : ---------------------------- [2024-03-14 14:13:09,948] Subj001 - INFO - emd.logger:195 - set_level() : EMD logger: handler 'console' level set to 'INFO' [2024-03-14 14:13:09,948] Subj001 - INFO - emd.logger:261 - sift_logger() : STARTED: mask_sift [2024-03-14 14:13:09,948] Subj001 - INFO - emd.sift:1084 - get_mask_freqs() : Computing first mask frequency with method zc [2024-03-14 14:13:09,948] Subj001 - INFO - emd.sift:1085 - get_mask_freqs() : Getting first IMF with no mask [2024-03-14 14:13:09,954] Subj001 - INFO - emd.sift:1093 - get_mask_freqs() : Found first mask frequency of 0.33955078125 [2024-03-14 14:13:09,954] Subj001 - INFO - emd.sift:1259 - mask_sift() : Sifting IMF-0 [2024-03-14 14:13:09,954] Subj001 - INFO - emd.sift:1030 - get_next_imf_mask() : Defining masks with freq 0.33955078125 and amp 8236.609552003989 at 4 phases [2024-03-14 14:13:09,975] Subj001 - INFO - emd.sift:1259 - mask_sift() : Sifting IMF-1 [2024-03-14 14:13:09,976] Subj001 - INFO - emd.sift:1030 - get_next_imf_mask() : Defining masks with freq 0.169775390625 and amp 8236.609552003989 at 4 phases [2024-03-14 14:13:09,993] Subj001 - INFO - emd.sift:1259 - mask_sift() : Sifting IMF-2 [2024-03-14 14:13:09,993] Subj001 - INFO - emd.sift:1030 - get_next_imf_mask() : Defining masks with freq 0.0848876953125 and amp 8236.609552003989 at 4 phases [2024-03-14 14:13:10,008] Subj001 - INFO - emd.sift:1259 - mask_sift() : Sifting IMF-3 [2024-03-14 14:13:10,008] Subj001 - INFO - emd.sift:1030 - get_next_imf_mask() : Defining masks with freq 0.04244384765625 and amp 8236.609552003989 at 4 phases [2024-03-14 14:13:10,023] Subj001 - INFO - emd.sift:1259 - mask_sift() : Sifting IMF-4 [2024-03-14 14:13:10,023] Subj001 - INFO - emd.sift:1030 - get_next_imf_mask() : Defining masks with freq 0.021221923828125 and amp 8236.609552003989 at 4 phases [2024-03-14 14:13:10,032] Subj001 - INFO - emd.sift:1259 - mask_sift() : Sifting IMF-5 [2024-03-14 14:13:10,032] Subj001 - INFO - emd.sift:1030 - get_next_imf_mask() : Defining masks with freq 0.0106109619140625 and amp 8236.609552003989 at 4 phases [2024-03-14 14:13:10,040] Subj001 - INFO - emd.sift:1259 - mask_sift() : Sifting IMF-6 [2024-03-14 14:13:10,040] Subj001 - INFO - emd.sift:1030 - get_next_imf_mask() : Defining masks with freq 0.00530548095703125 and amp 8236.609552003989 at 4 phases [2024-03-14 14:13:10,048] Subj001 - INFO - emd.sift:1259 - mask_sift() : Sifting IMF-7 [2024-03-14 14:13:10,048] Subj001 - INFO - emd.sift:1030 - get_next_imf_mask() : Defining masks with freq 0.002652740478515625 and amp 8236.609552003989 at 4 phases [2024-03-14 14:13:10,056] Subj001 - INFO - emd.sift:1259 - mask_sift() : Sifting IMF-8 [2024-03-14 14:13:10,056] Subj001 - INFO - emd.sift:1030 - get_next_imf_mask() : Defining masks with freq 0.0013263702392578125 and amp 8236.609552003989 at 4 phases [2024-03-14 14:13:10,065] Subj001 - INFO - emd.sift:1276 - mask_sift() : Finishing sift: reached max number of imfs (9) [2024-03-14 14:13:10,065] Subj001 - INFO - emd.logger:287 - sift_logger() : COMPLETED: mask_sift [2024-03-14 14:13:10,065] Subj001 - INFO - emd.logger:195 - set_level() : EMD logger: handler 'console' level set to 'DEBUG' [2024-03-14 14:13:10,065] Subj001 - INFO - emd.spectra:84 - frequency_transform() : STARTED: compute frequency stats [2024-03-14 14:13:10,065] Subj001 - DEBUG - emd.spectra:87 - frequency_transform() : computing on 5120 samples over 9 imfs at sample rate 512 [2024-03-14 14:13:10,065] Subj001 - INFO - emd.spectra:102 - frequency_transform() : Using Amplitude-Normalised Hilbert transform [2024-03-14 14:13:10,066] Subj001 - INFO - emd.imftools:58 - amplitude_normalise() : STARTED: Amplitude-Normalise [2024-03-14 14:13:10,066] Subj001 - DEBUG - emd.imftools:61 - amplitude_normalise() : Normalising 5120 samples across 9 IMFs [2024-03-14 14:13:10,066] Subj001 - DEBUG - emd.imftools:64 - amplitude_normalise() : Using pchip interpolation with threshold of 1e-10 and max_iters 3 [2024-03-14 14:13:10,066] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,066] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 1611 extrema on mode abs_peaks [2024-03-14 14:13:10,066] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,067] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,067] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 1652 extrema on mode abs_peaks [2024-03-14 14:13:10,067] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,068] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,068] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 1653 extrema on mode abs_peaks [2024-03-14 14:13:10,068] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,069] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,069] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 1653 extrema on mode abs_peaks [2024-03-14 14:13:10,069] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,070] Subj001 - INFO - emd.imftools:106 - amplitude_normalise() : Normalise of IMF-0-0 complete in 3 iters (val=0.0) [2024-03-14 14:13:10,070] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,070] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 1437 extrema on mode abs_peaks [2024-03-14 14:13:10,071] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,071] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,072] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 1442 extrema on mode abs_peaks [2024-03-14 14:13:10,072] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,073] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,073] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 1442 extrema on mode abs_peaks [2024-03-14 14:13:10,073] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,074] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,074] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 1442 extrema on mode abs_peaks [2024-03-14 14:13:10,074] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,075] Subj001 - INFO - emd.imftools:106 - amplitude_normalise() : Normalise of IMF-1-0 complete in 3 iters (val=0.0) [2024-03-14 14:13:10,075] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,075] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 763 extrema on mode abs_peaks [2024-03-14 14:13:10,075] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,076] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,076] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 764 extrema on mode abs_peaks [2024-03-14 14:13:10,076] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,077] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,077] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 764 extrema on mode abs_peaks [2024-03-14 14:13:10,077] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,078] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,078] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 764 extrema on mode abs_peaks [2024-03-14 14:13:10,078] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,079] Subj001 - INFO - emd.imftools:106 - amplitude_normalise() : Normalise of IMF-2-0 complete in 3 iters (val=0.0) [2024-03-14 14:13:10,079] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,079] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 260 extrema on mode abs_peaks [2024-03-14 14:13:10,079] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,080] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,080] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 260 extrema on mode abs_peaks [2024-03-14 14:13:10,080] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,081] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,081] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 260 extrema on mode abs_peaks [2024-03-14 14:13:10,081] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,082] Subj001 - INFO - emd.imftools:106 - amplitude_normalise() : Normalise of IMF-3-0 complete in 2 iters (val=0.0) [2024-03-14 14:13:10,082] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,082] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 228 extrema on mode abs_peaks [2024-03-14 14:13:10,082] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,083] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,083] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 228 extrema on mode abs_peaks [2024-03-14 14:13:10,083] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,084] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,084] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 228 extrema on mode abs_peaks [2024-03-14 14:13:10,084] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,085] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,085] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 228 extrema on mode abs_peaks [2024-03-14 14:13:10,085] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,086] Subj001 - INFO - emd.imftools:106 - amplitude_normalise() : Normalise of IMF-4-0 complete in 3 iters (val=0.0) [2024-03-14 14:13:10,086] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,086] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 87 extrema on mode abs_peaks [2024-03-14 14:13:10,086] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,087] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,087] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 87 extrema on mode abs_peaks [2024-03-14 14:13:10,087] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,088] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,088] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 87 extrema on mode abs_peaks [2024-03-14 14:13:10,088] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,089] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,089] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 87 extrema on mode abs_peaks [2024-03-14 14:13:10,089] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,090] Subj001 - INFO - emd.imftools:106 - amplitude_normalise() : Normalise of IMF-5-0 complete in 3 iters (val=0.0) [2024-03-14 14:13:10,090] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,090] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 44 extrema on mode abs_peaks [2024-03-14 14:13:10,090] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,091] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,091] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 44 extrema on mode abs_peaks [2024-03-14 14:13:10,091] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,092] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,092] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 44 extrema on mode abs_peaks [2024-03-14 14:13:10,092] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,093] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,093] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 44 extrema on mode abs_peaks [2024-03-14 14:13:10,093] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,094] Subj001 - INFO - emd.imftools:106 - amplitude_normalise() : Normalise of IMF-6-0 complete in 3 iters (val=0.0) [2024-03-14 14:13:10,094] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,094] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 22 extrema on mode abs_peaks [2024-03-14 14:13:10,094] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,095] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,095] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 22 extrema on mode abs_peaks [2024-03-14 14:13:10,095] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,096] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,096] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 22 extrema on mode abs_peaks [2024-03-14 14:13:10,096] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,097] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,097] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 22 extrema on mode abs_peaks [2024-03-14 14:13:10,097] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,098] Subj001 - INFO - emd.imftools:106 - amplitude_normalise() : Normalise of IMF-7-0 complete in 3 iters (val=0.0) [2024-03-14 14:13:10,098] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,098] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 10 extrema on mode abs_peaks [2024-03-14 14:13:10,098] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,099] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,099] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 10 extrema on mode abs_peaks [2024-03-14 14:13:10,099] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,100] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'combined' with method 'pchip' [2024-03-14 14:13:10,100] Subj001 - DEBUG - emd._sift_core:92 - get_padded_extrema() : found 10 extrema on mode abs_peaks [2024-03-14 14:13:10,100] Subj001 - VERBOSE - emd.logger:91 - log_for_level() : Padding 2 extrema in signal X 5120 using method 'numpypad' [2024-03-14 14:13:10,101] Subj001 - INFO - emd.imftools:106 - amplitude_normalise() : Normalise of IMF-8-0 complete in 2 iters (val=0.0) [2024-03-14 14:13:10,101] Subj001 - INFO - emd.imftools:119 - amplitude_normalise() : COMPLETED: Amplitude-Normalise [2024-03-14 14:13:10,102] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'upper' with method 'splrep' [2024-03-14 14:13:10,103] Subj001 - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 1826 minima and 1826 maxima on mode peaks [2024-03-14 14:13:10,103] Subj001 - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:10,103] Subj001 - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:10,103] Subj001 - DEBUG - emd._sift_core:287 - _pad_extrema_rilling() : R: max earlier than min, last val larger than last max [2024-03-14 14:13:10,105] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'upper' with method 'splrep' [2024-03-14 14:13:10,105] Subj001 - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 777 minima and 777 maxima on mode peaks [2024-03-14 14:13:10,105] Subj001 - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:10,105] Subj001 - DEBUG - emd._sift_core:271 - _pad_extrema_rilling() : L: max later than min, first val smaller than first max [2024-03-14 14:13:10,105] Subj001 - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:10,106] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'upper' with method 'splrep' [2024-03-14 14:13:10,106] Subj001 - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 395 minima and 396 maxima on mode peaks [2024-03-14 14:13:10,106] Subj001 - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:10,106] Subj001 - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:10,106] Subj001 - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:10,107] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'upper' with method 'splrep' [2024-03-14 14:13:10,107] Subj001 - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 132 minima and 133 maxima on mode peaks [2024-03-14 14:13:10,107] Subj001 - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:10,107] Subj001 - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:10,108] Subj001 - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:10,108] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'upper' with method 'splrep' [2024-03-14 14:13:10,108] Subj001 - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 116 minima and 115 maxima on mode peaks [2024-03-14 14:13:10,109] Subj001 - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:10,109] Subj001 - DEBUG - emd._sift_core:271 - _pad_extrema_rilling() : L: max later than min, first val smaller than first max [2024-03-14 14:13:10,109] Subj001 - DEBUG - emd._sift_core:281 - _pad_extrema_rilling() : R: max earlier than min, last val smaller than last max [2024-03-14 14:13:10,109] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'upper' with method 'splrep' [2024-03-14 14:13:10,110] Subj001 - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 45 minima and 45 maxima on mode peaks [2024-03-14 14:13:10,110] Subj001 - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:10,110] Subj001 - DEBUG - emd._sift_core:256 - _pad_extrema_rilling() : L: max earlier than min, first val smaller than first min [2024-03-14 14:13:10,110] Subj001 - DEBUG - emd._sift_core:281 - _pad_extrema_rilling() : R: max earlier than min, last val smaller than last max [2024-03-14 14:13:10,110] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'upper' with method 'splrep' [2024-03-14 14:13:10,111] Subj001 - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 23 minima and 24 maxima on mode peaks [2024-03-14 14:13:10,111] Subj001 - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:10,111] Subj001 - DEBUG - emd._sift_core:256 - _pad_extrema_rilling() : L: max earlier than min, first val smaller than first min [2024-03-14 14:13:10,111] Subj001 - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:10,112] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'upper' with method 'splrep' [2024-03-14 14:13:10,112] Subj001 - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 11 minima and 12 maxima on mode peaks [2024-03-14 14:13:10,112] Subj001 - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:10,112] Subj001 - DEBUG - emd._sift_core:256 - _pad_extrema_rilling() : L: max earlier than min, first val smaller than first min [2024-03-14 14:13:10,112] Subj001 - DEBUG - emd._sift_core:295 - _pad_extrema_rilling() : R: max later than min, last val larger than last min [2024-03-14 14:13:10,113] Subj001 - DEBUG - emd._sift_core:484 - interp_envelope() : Interpolating 'upper' with method 'splrep' [2024-03-14 14:13:10,113] Subj001 - DEBUG - emd._sift_core:78 - get_padded_extrema() : found 6 minima and 6 maxima on mode peaks [2024-03-14 14:13:10,113] Subj001 - DEBUG - emd._sift_core:239 - _pad_extrema_rilling() : Padding 2 extrema in signal X (5120,) using method 'rilling' [2024-03-14 14:13:10,113] Subj001 - DEBUG - emd._sift_core:250 - _pad_extrema_rilling() : L: max earlier than min, first val larger than first min [2024-03-14 14:13:10,113] Subj001 - DEBUG - emd._sift_core:281 - _pad_extrema_rilling() : R: max earlier than min, last val smaller than last max [2024-03-14 14:13:10,120] Subj001 - INFO - emd.spectra:193 - frequency_transform() : COMPLETED: compute frequency stats. Returning 9 imfs [2024-03-14 14:13:10,120] Subj001 - INFO - emd.emd_tutorial_04_utils_02_logger:157 - my_analysis() : Avg frequency of IMF-2 is 36.310268Hz [2024-03-14 14:13:10,120] Subj001 - INFO - emd.cycles:131 - get_cycle_vector() : STARTED: get cycle indices [2024-03-14 14:13:10,120] Subj001 - DEBUG - emd.support:261 - ensure_2d() : Checking get_cycle_vector inputs - Adding dummy dimension to input 'mask' [2024-03-14 14:13:10,120] Subj001 - DEBUG - emd.cycles:138 - get_cycle_vector() : computing on 5120 samples over 9 IMFs [2024-03-14 14:13:10,120] Subj001 - DEBUG - emd.cycles:141 - get_cycle_vector() : 5027 (98.18%) samples masked out [2024-03-14 14:13:10,125] Subj001 - INFO - emd.cycles:191 - get_cycle_vector() : found 1 cycles in IMF-0 [2024-03-14 14:13:10,138] Subj001 - INFO - emd.cycles:191 - get_cycle_vector() : found 53 cycles in IMF-1 [2024-03-14 14:13:10,146] Subj001 - INFO - emd.cycles:191 - get_cycle_vector() : found 105 cycles in IMF-2 [2024-03-14 14:13:10,149] Subj001 - INFO - emd.cycles:191 - get_cycle_vector() : found 100 cycles in IMF-3 [2024-03-14 14:13:10,152] Subj001 - INFO - emd.cycles:191 - get_cycle_vector() : found 92 cycles in IMF-4 [2024-03-14 14:13:10,153] Subj001 - INFO - emd.cycles:191 - get_cycle_vector() : found 27 cycles in IMF-5 [2024-03-14 14:13:10,154] Subj001 - INFO - emd.cycles:191 - get_cycle_vector() : found 7 cycles in IMF-6 [2024-03-14 14:13:10,154] Subj001 - INFO - emd.cycles:191 - get_cycle_vector() : found 1 cycles in IMF-7 [2024-03-14 14:13:10,155] Subj001 - INFO - emd.cycles:191 - get_cycle_vector() : found 0 cycles in IMF-8 [2024-03-14 14:13:10,155] Subj001 - INFO - emd.cycles:193 - get_cycle_vector() : COMPLETED: get cycle indices [2024-03-14 14:13:10,155] Subj001 - INFO - emd.cycles:423 - get_cycle_stat() : STARTED: get_cycle_stat [2024-03-14 14:13:10,158] Subj001 - INFO - emd.cycles:423 - get_cycle_stat() : STARTED: get_cycle_stat [2024-03-14 14:13:10,161] Subj001 - INFO - emd.emd_tutorial_04_utils_02_logger:168 - my_analysis() : Freq-Amp correlation: r=-0.366304 [2024-03-14 14:13:10,161] Subj001 - INFO - emd.emd_tutorial_04_utils_02_logger:173 - my_analysis() : My new analysis finished in 0.212678 seconds .. rst-class:: sphx-glr-timing **Total running time of the script:** (0 minutes 1.273 seconds) .. _sphx_glr_download_emd_tutorials_04_utilities_emd_tutorial_04_utils_02_logger.py: .. only:: html .. container:: sphx-glr-footer sphx-glr-footer-example .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: emd_tutorial_04_utils_02_logger.ipynb ` .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: emd_tutorial_04_utils_02_logger.py ` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_