add option in "mne report" to get sensitivity maps of forward models

[agramfort]

it's a good quality insurance plot for forward models.

[agramfort]

something along these lines (copied from a script of @dengemann)

mag_map = mne.sensitivity_map(fwd, ch_type=fwd_ch_type, mode=fwd_ori)
brain = mag_map.plot(subject=subject, time_label='Magnetometer sensitivity',
                     hemi='rh', subjects_dir=subjects_dir, **fwd_plot_args)
brain.save_image('sensitivity_mag.png')
report.add_images_to_section('sensitivity_mag.png', 'sensitivity map',
                             'forward')
os.remove('sensitivity_mag.png')

fig = plt.figure()
plt.hist(mag_map.data.ravel(),
         bins=20, label=['Magnetometers'],
         color=['b'])
plt.title('Normal orientation sensitivity')
plt.xlabel('sensitivity')
plt.ylabel('count')
plt.legend()

[jasmainak]

+1

[larsoner]

@jasmainak do you have time for this in the next two weeks?

[jasmainak]

okay I'll give some time for fixing MNE stuff on Monday. I've a couple of things to do.

[dengemann]

+1

On Fri, Sep 9, 2016 at 2:30 PM Mainak Jas notifications@github.com wrote:

okay I'll give some time for fixing MNE stuff on Monday. I've a couple of
things to do.

—
You are receiving this because you were mentioned.
Reply to this email directly, view it on GitHub
#2828 (comment),
or mute the thread
https://github.com/notifications/unsubscribe-auth/AB0fiiaVcWckZu6BFXgoJmmjaOHSIgVkks5qoVFOgaJpZM4HPzH3
.

[larsoner]

@jasmainak any time for this one in the next couple of days? If not I'll remove the milestone (or you can)

[jasmainak]

I have a local branch @Eric89GXL but I fear it's going to take more time. I'll remove it from the milestone

[parulsethi]

hi,
I’m interested in solving this issue, but will need some help through it.

If I understand correctly, a Report class method needs to be added to append 'sensitivity map' to the report? and is ‘sensitivity/count histogram’ mentioned above also required to be appended?

and if a method, will the input parameters to it would be same as that of 'mne.sensitivity_map()+SourceEstimate.plot()' or something else?

[jasmainak]

hi @parulsethi I'd suggest that you start from my work in progress commit here. I don't think I got it working, but it might be easier than starting from scratch.

[jasmainak]

Also, make sure you understand how to generate a report. You can read the tutorial on the mne report here

[parulsethi]

@jasmainak thanks for the suggestion, I've forked from there.