<p><span style="padding: 3px 10px; border-radius: 5px; color: #ffffff; font-weight: bold; display: inline-block; background-color: #ff0000;"> External Email - Use Caution </span></p><p></p><div dir="ltr">hi,<div><br></div><div>here are the lines that generate this behavior.</div><div><br></div><div><a href="https://secure-web.cisco.com/19bmNMoYnXLOkU_iFnakp_XeteNEgsBBxfAGwS0rF2CpZNdLrEeQ2Q7gkaBm6CPlzSZf0ojCFXZ9hmC1DsZAEgb0US11kbgMbsLoCOjXYEOf-NN97XCe29r1vUCh20gEOZLm2BIVUjGq3PJnp4npGEOWbrB5lG3fytl7LKURPdRZnu695xfBaNgSLPXjQqDRJ6MW6-NU8lSF9eSsZErLm6Xs_wM9LU9iGRhhtFMZfRgwj_W4Z4Ug6Ds8xnaQsostgCuEzh10kChT37jmHmXJKlQ/https%3A%2F%2Fgithub.com%2Fmne-tools%2Fmne-python%2Fblob%2Fmaster%2Fmne%2Ftime_frequency%2F_stockwell.py%23L26"><font color="red"><b>MailScanner has detected a possible fraud attempt from "secure-web.cisco.com" claiming to be</b></font> https://github.com/mne-tools/mne-python/blob/master/mne/time_frequency/_stockwell.py#L26</a><br></div><div><br></div><div>I read that "n_fft cannot be smaller than signal size". </div><div><br></div><div>maybe simply the doc needs to be updated?</div><div><br></div><div>HTH</div><div>Alex</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Dec 20, 2020 at 6:32 PM Paul Fishback <<a href="mailto:fishbacp@mail.gvsu.edu">fishbacp@mail.gvsu.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><p><span style="padding:3px 10px;border-radius:5px;color:rgb(255,255,255);font-weight:bold;display:inline-block;background-color:rgb(255,0,0)"> External Email - Use Caution </span></p><p></p><div dir="ltr"><div dir="ltr">Alex,<div><br></div><div>I'm so sorry I left that out. Here's a snippet that includes that info as well as plots the results. I set n_fft=500, but still received the warning, <div><font color="#ff0000">RuntimeWarning: The input signal is shorter (5001) than "n_fft" (8192). Applying zero padding</font>.. <br><div><br></div><div>import numpy as np<br>from scipy.signal import chirp<br>import matplotlib.pyplot as plt<br>from mne.time_frequency import tfr_array_stockwell<br><br>Nt=5000<br>time_length=10<br><br>t = np.linspace(0, time_length, Nt+1)<br>w = chirp(t, f0=12.5, f1=2.5, t1=10, method='linear')<br>y = np.expand_dims(w, axis=[0,1])<br><br>fmin = 0 # Hz<br>fmax = 25 # Hz<br>Fs = int(np.ceil(Nt/time_length))<br><br>st_power, itc, freqs=tfr_array_stockwell(y,sfreq=Fs,fmin=fmin,fmax=fmax,n_fft=500)<br><br>extent = (t[0], t[-1], fmin, fmax)<br>fig, ax = plt.subplots(2, 1, sharex=True)<br>ax[0].plot(t, w)<br>ax[0].set(ylabel='amplitude')<br>ax[1].imshow(np.abs(st_power[0,:,:]), origin='lower', extent=extent)<br>ax[1].axis('tight')<br>ax[1].set(xlabel='time (s)', ylabel='frequency (Hz)')<br>plt.show()<br></div></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Dec 20, 2020 at 12:06 PM Alexandre Gramfort <<a href="mailto:alexandre.gramfort@inria.fr" target="_blank">alexandre.gramfort@inria.fr</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><p><span style="padding:3px 10px;border-radius:5px;color:rgb(255,255,255);font-weight:bold;display:inline-block;background-color:rgb(255,0,0)"> External Email - Use Caution </span></p><p></p><div dir="ltr">hi Paul,<div><br></div><div>can you share a full code snippet?</div><div><br></div><div>Nt and time_length not defined</div><div><br></div><div>Alex</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Dec 19, 2020 at 9:42 PM Paul Fishback <<a href="mailto:fishbacp@mail.gvsu.edu" target="_blank">fishbacp@mail.gvsu.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><p><span style="padding:3px 10px;border-radius:5px;color:rgb(255,255,255);font-weight:bold;display:inline-block;background-color:rgb(255,0,0)"> External Email - Use Caution </span></p><p></p><div dir="ltr">Hello,<div><br></div><div>I'm interested in computing the Stockwell Transform of a simple univariate time series. Because I already use MNE, I thought I'd perform this task with the aid of the </div><div>tfr_array_stockwell command. However, I'm having a problem changing the window length, n_fft. <br></div><div><br></div><div>As a simple "toy example," here's a chirp signal, where the frequency decreases linearly from 12.5 to 2.5 Hz over a 10 second period:</div><div><br></div><div>from scipy.signal import chirp<br>from mne.time_frequency import tfr_array_stockwell<br></div><div><br>t = np.linspace(0, 10, 5001)<br>w = chirp(t, f0=12.5, f1=2.5, t1=10, method='linear')<br>y = np.expand_dims(w, axis=[0,1]) # Add dimensions to signal in order to make a valid input for tfr_array_stockwell<br><br>fmin = 0 <br>fmax = 25 <br>Fs = int(np.ceil(Nt/time_length))</div><div>st_power, itc, freqs=tfr_array_stockwell(y,sfreq=Fs,fmin=fmin,fmax=fmax)<br></div><div><br></div><div>This works fine and the plotted results look great, but I do receive a warning message indicating zero padding is being used. This is because the signal length (5001)</div><div>is less than the default n_fft, which is 8192 in my case. (Since 2^(12)=4096 < 5001 < 2^(13)=8192 )</div><div><br></div><div>However, if I change the default n_fft to something else, e.g. insert n_fft=500 above, I still receive the same warning message. </div><div><br></div><div>Perhaps I'm missing something simple?</div><div><br></div><div>Thanks</div><div><br></div><div>Paul Fishback</div><div>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><font face="arial, sans-serif">Professor of Mathematics, Grand Valley State University </font></div><div dir="ltr"><font><font face="arial, sans-serif">Past-President, Pi Mu Epsilon National Honorary Mathematics Society, Inc.</font><br></font><br><div style="font-family:arial">Department of Mathematics (MAK C-2-408)</div><div style="font-family:arial">Grand Valley State University</div><div style="font-family:arial">1 Campus Dr.<br>Allendale, MI 49401<br><font color="#0000ff"><a href="mailto:fishbacp@mail.gvsu.edu" target="_blank">fishbacp@mail.gvsu.edu</a></font></div><div><font color="#0000ff"><span style="font-family:arial"><span title="Call with Google Voice"><span title="Call with Google Voice"><span title="Call with Google Voice">616.331.2040</span></span></span></span><br></font><span style="font-family:arial"><font color="#0000ff"><span title="Call with Google Voice"><span title="Call with Google Voice"><span title="Call with Google Voice">616.331.3120</span></span></span> </font>(fax)</span><br></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
_______________________________________________<br>
Mne_analysis mailing list<br>
<a href="mailto:Mne_analysis@nmr.mgh.harvard.edu" target="_blank">Mne_analysis@nmr.mgh.harvard.edu</a><br>
<a href="https://secure-web.cisco.com/1M4CbzTtCaJem1lx1U0buS9f1c2izKtbTjVFnrbjB0yPGi90XuU1Z5zbH7fAswJDcFljz4nDFhhWdKbqEcuR5SoCem1UyfXaKE-shx-JV4-lWrS1HL6eh8opMdufAMu0eO3z895xxr9vmgvNQClV-Rsi0DHYSf5fpW22GZbfJLPXPzfJjpHZB0JXycCZwjO4lHoIoqa-TwctE2IbhjzWXmpQZpxM2Yq3Ur1j_cifW1UHMNXxoIbVqgFtXv07-6Nd2_shMG3tHtbjFhbzfyYuLxA/https%3A%2F%2Fmail.nmr.mgh.harvard.edu%2Fmailman%2Flistinfo%2Fmne_analysis" rel="noreferrer" target="_blank"><font color="red"><b>MailScanner has detected a possible fraud attempt from "secure-web.cisco.com" claiming to be</b></font> https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis</a></blockquote></div>
_______________________________________________<br>
Mne_analysis mailing list<br>
<a href="mailto:Mne_analysis@nmr.mgh.harvard.edu" target="_blank">Mne_analysis@nmr.mgh.harvard.edu</a><br>
<a href="https://secure-web.cisco.com/1ctCQ2zsYnftewzKaPfkHy8euSz0Z_C9qeWYG9Vwvew04AcXo9LLpUkkieMGqt1MCu4FTZbv3dDBHdWMVCUYTWPRazAgJVP7pTD1ohCnZATpUNE83BYx3f3sX5WnGVrIDidqjndItRhUBtzXxKVVqghi8ZBZv2VbeKb3z5tCxIOCXBxM2qI_dyxZm9CP08QmgRXVcLePiwX8MPWmFP_bb024AGt5RxjY9LRXWVmEec350uCGrnjOr4MfyAHK5n0O3LC9isssbFfny5JO6oFwchA/https%3A%2F%2Fmail.nmr.mgh.harvard.edu%2Fmailman%2Flistinfo%2Fmne_analysis" rel="noreferrer" target="_blank"><font color="red"><b>MailScanner has detected a possible fraud attempt from "secure-web.cisco.com" claiming to be</b></font> https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis</a></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><font face="arial, sans-serif">Professor of Mathematics, Grand Valley State University </font></div><div dir="ltr"><font><font face="arial, sans-serif">Past-President, Pi Mu Epsilon National Honorary Mathematics Society, Inc.</font><br></font><br><div style="font-family:arial">Department of Mathematics (MAK C-2-408)</div><div style="font-family:arial">Grand Valley State University</div><div style="font-family:arial">1 Campus Dr.<br>Allendale, MI 49401<br><font color="#0000ff"><a href="mailto:fishbacp@mail.gvsu.edu" target="_blank">fishbacp@mail.gvsu.edu</a></font></div><div><font color="#0000ff"><span style="font-family:arial"><span title="Call with Google Voice"><span title="Call with Google Voice"><span title="Call with Google Voice">616.331.2040</span></span></span></span><br></font><span style="font-family:arial"><font color="#0000ff"><span title="Call with Google Voice"><span title="Call with Google Voice"><span title="Call with Google Voice">616.331.3120</span></span></span> </font>(fax)</span><br></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
_______________________________________________<br>
Mne_analysis mailing list<br>
<a href="mailto:Mne_analysis@nmr.mgh.harvard.edu" target="_blank">Mne_analysis@nmr.mgh.harvard.edu</a><br>
<a href="https://secure-web.cisco.com/1V2HoZGdNLp0IFtAYTXTEdVFL6qTpGgJRm2OCrOWu-BPKeW0GOZYTT50BA5CqqaLL4F7Cia6pZ1bgXjpEBbPf4tQcDUwPRFmQ4SO76LkQweDD5ZEX6rWY--QfBPCCK2oKkxrMznXqgSplhMudB93vKc_c1flqMbs5HYUoQlnqldVFc9weJ4urtRK4Jv-1lKttxMvCWdmSfTdYS5kYG3y7QI_RsxsiBBzL28XDuNDaMaRP8qaIuEpabwfeoA7Qspj65b9UzXKBiPlHdjeSGHbgWA/https%3A%2F%2Fmail.nmr.mgh.harvard.edu%2Fmailman%2Flistinfo%2Fmne_analysis" rel="noreferrer" target="_blank"><font color="red"><b>MailScanner has detected a possible fraud attempt from "secure-web.cisco.com" claiming to be</b></font> https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis</a></blockquote></div>