本文整理汇总了Python中mne.io.Raw.fix_mag_coil_types方法的典型用法代码示例。如果您正苦于以下问题:Python Raw.fix_mag_coil_types方法的具体用法?Python Raw.fix_mag_coil_types怎么用?Python Raw.fix_mag_coil_types使用的例子?那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类mne.io.Raw
的用法示例。
在下文中一共展示了Raw.fix_mag_coil_types方法的2个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: test_fix_types
# 需要导入模块: from mne.io import Raw [as 别名]
# 或者: from mne.io.Raw import fix_mag_coil_types [as 别名]
def test_fix_types():
"""Test fixing of channel types
"""
for fname, change in ((hp_fif_fname, True), (test_fif_fname, False), (ctf_fname, False)):
raw = Raw(fname)
mag_picks = pick_types(raw.info, meg="mag")
other_picks = np.setdiff1d(np.arange(len(raw.ch_names)), mag_picks)
# we don't actually have any files suffering from this problem, so
# fake it
if change:
for ii in mag_picks:
raw.info["chs"][ii]["coil_type"] = FIFF.FIFFV_COIL_VV_MAG_T2
orig_types = np.array([ch["coil_type"] for ch in raw.info["chs"]])
raw.fix_mag_coil_types()
new_types = np.array([ch["coil_type"] for ch in raw.info["chs"]])
if not change:
assert_array_equal(orig_types, new_types)
else:
assert_array_equal(orig_types[other_picks], new_types[other_picks])
assert_true((orig_types[mag_picks] != new_types[mag_picks]).all())
assert_true((new_types[mag_picks] == FIFF.FIFFV_COIL_VV_MAG_T3).all())
示例2: test_triux
# 需要导入模块: from mne.io import Raw [as 别名]
# 或者: from mne.io.Raw import fix_mag_coil_types [as 别名]
def test_triux():
"""Test TRIUX system support"""
raw = Raw(tri_fname).crop(0, 0.999)
raw.fix_mag_coil_types()
# standard
sss_py = maxwell_filter(raw, coord_frame='meg', regularize=None)
assert_meg_snr(sss_py, Raw(tri_sss_fname), 37, 700)
# cross-talk
sss_py = maxwell_filter(raw, coord_frame='meg', regularize=None,
cross_talk=tri_ctc_fname)
assert_meg_snr(sss_py, Raw(tri_sss_ctc_fname), 35, 700)
# fine cal
sss_py = maxwell_filter(raw, coord_frame='meg', regularize=None,
calibration=tri_cal_fname)
assert_meg_snr(sss_py, Raw(tri_sss_cal_fname), 31, 360)
# ctc+cal
sss_py = maxwell_filter(raw, coord_frame='meg', regularize=None,
calibration=tri_cal_fname,
cross_talk=tri_ctc_fname)
assert_meg_snr(sss_py, Raw(tri_sss_ctc_cal_fname), 31, 350)
# regularization
sss_py = maxwell_filter(raw, coord_frame='meg', regularize='in')
sss_mf = Raw(tri_sss_reg_fname)
assert_meg_snr(sss_py, sss_mf, 0.6, 9)
_check_reg_match(sss_py, sss_mf, 1)
# all three
sss_py = maxwell_filter(raw, coord_frame='meg', regularize='in',
calibration=tri_cal_fname,
cross_talk=tri_ctc_fname)
sss_mf = Raw(tri_sss_ctc_cal_reg_in_fname)
assert_meg_snr(sss_py, sss_mf, 0.6, 9)
_check_reg_match(sss_py, sss_mf, 1)
# tSSS
raw = Raw(tri_fname).fix_mag_coil_types()
sss_py = maxwell_filter(raw, coord_frame='meg', regularize=None,
st_duration=4., verbose=True)
assert_meg_snr(sss_py, Raw(tri_sss_st4_fname), 700., 1600)