def test_io_cov(): """Test IO for noise covariance matrices """ cov = read_cov(cov_fname) cov.save(op.join(tempdir, 'cov.fif')) cov2 = read_cov(op.join(tempdir, 'cov.fif')) assert_array_almost_equal(cov.data, cov2.data) cov2 = read_cov(cov_gz_fname) assert_array_almost_equal(cov.data, cov2.data) cov2.save(op.join(tempdir, 'cov.fif.gz')) cov2 = read_cov(op.join(tempdir, 'cov.fif.gz')) assert_array_almost_equal(cov.data, cov2.data) cov['bads'] = ['EEG 039'] cov_sel = pick_channels_cov(cov, exclude=cov['bads']) assert_true(cov_sel['dim'] == (len(cov['data']) - len(cov['bads']))) assert_true(cov_sel['data'].shape == (cov_sel['dim'], cov_sel['dim'])) cov_sel.save(op.join(tempdir, 'cov.fif')) cov2 = read_cov(cov_gz_fname) assert_array_almost_equal(cov.data, cov2.data) cov2.save(op.join(tempdir, 'cov.fif.gz')) cov2 = read_cov(op.join(tempdir, 'cov.fif.gz')) assert_array_almost_equal(cov.data, cov2.data)