raw_empty.interpolate_bads()# pick MEG channels onlypicks = pick_types(raw_empty.info, meg=True, ref_meg=False, eeg=False, stim=False, eog=False, exclude='bads')# calculate noise-covariance matrixnoise_cov_mat = cp_covariance(raw_empty, picks=picks, verbose=verbose)# write noise-covarian...