Hi, I’m getting different unit ranges when doing psd.
When I epoch raw data with equal length segments and do a plot_psd_topomap per frequency bands, the range is different from the plot_psd one.
See last figures on this notebook.
I marked dB=True and Normalize=False. Any other thing that should be taking into account?
The topomaps show power within specific frequency bands, whereas the PSD shows density. I’m thinking that maybe the units should not be per Hz in the topomaps because you’re basically summing power at specific frequencies.
In addition, if the unit is dB it cannot also be e.g. fT2/Hz, so we should probably fix the labels.
@cbrnr After bringing this up in a 10-second convo with @agramfort it appears that the actual difference is that in the topomaps, the unit is scaled differently (it’s actually in T^2/Hz and not in fT^2/Hz), so we should fix this or at least correct the label. Note that so far this is only a suspicion and I have yet to verify this.
This doesn’t sound right to me. If it really were T^2/Hz, you would see nothing in the maps because the actual signal is orders of magnitudes smaller. If you consider what I suspected and add up the power e.g. within 0-4Hz, you’d actually get the values in the topomap. But I’m also not 100% and would need to confirm.
True, I didn’t see the minus. But we’re talking about dB so it looks like the denominator in the two plot is different. But the other comment is still valid, the unit is dB and cannot be fT^2/Hz (or T^2/Hz) at the same time.
But why would you want to do this? dB is a ratio of two power values on a log scale. It doesn’t matter if the original unit was fT/Hz or T/Hz, because the ratio will be the same.
To make it obvious what’s the “power” that’s considered here in the first place. I guess I’m more thinking about how to help beginners understand what they’re seeing here…
Hi! I think labels help a lot for units etc, so I wouldn’t removed them.
The issue you opened Richard dealing with the psd_plot, it’s instructive, I reproduced it as well using psd welch and multitapper funtions (old posts), but still I miss the relation with the plot_psd_topomap
plot_psd_topomap doesn’t apply the unit scaling. If you run my example code from the GH issue and drop the line that scales the values, you’ll get results that resemble what you have in the topomaps.
.