[Mne_analysis] Compute source power spectral density (PSD) in a label

[Vivek Sharma]

        External Email - Use Caution        

Hi,
In continuation with the above query.
If during the compute of psd the option db is set as false then how it can
be plotted in db.
stc = compute_source_psd(raw, inverse_operator, lambda2=1. / 9.,
method="dSPM",
                         tmin=tmin, tmax=tmax, fmin=fmin, fmax=fmax,
                         pick_ori="normal", n_fft=n_fft, label=label,
                         dB=False)

To plot power in db I did the following change in code:
plt.plot(stc.times, 10*np.log10(stc.data.T))
I need validation for the changes I have made.

On Thu, Dec 26, 2019 at 2:08 AM Vivek Sharma <vivek.sharma1510 at gmail.com>
wrote:

> Hi,
> In the following example
> (
> https://mne.tools/stable/auto_examples/time_frequency/plot_source_power_spectrum.html
> )
> why are we multiplying stc.times and stc.data.T by 1e3?
>
> Secondly, I would like to bring it to notice that in the same example the
> plot that is being obtained shows frequency as 10000, 20000, etc. and the
> unit is represented as Hz which I believe is wrong as it should be
> millihertz.
> I try to keep the frequency in Hz by making the following change in the
> code.
> plt.plot(stc.times, (1e3 * stc.data.T))
> is the change valid?
>
> --
> Vivek Sharma
>


-- 
Vivek Sharma
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.nmr.mgh.harvard.edu/pipermail/mne_analysis/attachments/20191226/1c8f7bef/attachment.html