Skip to content

Instantly share code, notes, and snippets.

@jhouck
jhouck / sample_phantom_distance_to_helmet.py
Created January 18, 2022 23:47
Test distance between phantom_otaniemi and sensor helmet
#!/usr/bin/env python3
import os
from mne.io import read_info
from mne.datasets.brainstorm import bst_phantom_elekta
from mne import get_meg_helmet_surf
from mne.datasets import fetch_phantom
from mne.io.constants import FIFF
from mne.viz import plot_alignment
from mne.transforms import Transform
@jhouck
jhouck / mne.sample.helmet-distance.py
Created January 3, 2022 20:14
Quick little script to get median distance between sensor helmet and head surface points (excluding the face)
import mne
from mne.io.constants import FIFF
import os.path as op
import numpy as np
from scipy.stats import rankdata
#data_path = mne.datasets.sample.data_path(update_path=False)
data_path = '/export/research/analysis/human/jhouck/shared/tools/MNE-sample-data/MNE-sample-data'
subjects_dir = op.join(data_path, 'subjects')
subject='sample'
@jhouck
jhouck / mne_quantify_motion.py
Last active July 27, 2019 17:23
Quantify motion from head position file
#!/usr/bin/env python3
# Rationale for using degrees (pitch/roll/yaw) interchangeably with mm is "borrowed" from afni/3dvolreg. Basically, at a radius of about
# 57 mm (57.2598), circumference=360 mm, and 1 degree of rotation == 1 mm shift on the spherical surface.
# See https://afni.nimh.nih.gov/afni/community/board/read.php?1,82641,82645#msg-82645
import os
import mne
from nibabel.eulerangles import mat2euler
from mne.fixes import einsum
import numpy as np
@jhouck
jhouck / check_coreg.py
Created June 30, 2019 14:58
Borrow a bunch of private functions from mne.gui to extract coreg QA measures
#!/usr/bin/env python3
# So basically what we want to do is:
# - load up dig source, MRISubjectSource, and create a MRIHeadWithFiducialsModel
# - load coreg transformation
# - apply transform
# - Compare distance
import scipy as sp
import numpy as np
from mne import coreg
@jhouck
jhouck / mne_sample_inverse_psd.py
Last active March 31, 2019 19:42
Show power spectra after inverse using norm vs. using radial component
#!/usr/bin/env python3
import numpy as np
import mne
from mne.datasets import sample
from mne import setup_volume_source_space
from mne.minimum_norm import make_inverse_operator, apply_inverse_raw
import matplotlib.pyplot as plt
import os.path as op
data_path = sample.data_path()
@jhouck
jhouck / mne_sample_test_complex.py
Last active February 8, 2019 23:57
Test output of apply_inverse in surface and volume using MNE sample data
#!/usr/bin/env python3
import mne
from mne import setup_volume_source_space
from mne.io import Raw
from mne.datasets import sample
from mne.minimum_norm import make_inverse_operator, apply_inverse_epochs, apply_inverse_raw
import numpy as np
data_path = sample.data_path()
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.