Homec4science
Diffusion iCAPs public (master) Browse

iCAPs/packages/spm12/external/mne/mne_ex_read_write_raw.m
2cb18aae5832master

Last Change
Raw File

mne_ex_read_write_raw.m
View Options

function mne_ex_read_write_raw(infile,outfile)
%
% function mne_ex_read_write_raw(infile,outfile);
%
% Read and write raw data in 60-sec blocks
%
%
% Author : Matti Hamalainen, MGH Martinos Center
% License : BSD 3-clause
%
global FIFF;
if isempty(FIFF)
FIFF = fiff_define_constants();
end
%
me = 'MNE:mne_ex_read_write_raw';
%
if nargin ~= 2
error(me,'Incorrect number of arguments');
end
%
% Setup for reading the raw data
%
try
raw = fiff_setup_read_raw(infile);
catch
error(me,'%s',mne_omit_first_line(lasterr));
end
%raw.info.projs = [];
%
% Set up pick list: MEG + STI 014 - bad channels
%
%
want_meg = true;
want_eeg = false;
want_stim = false;
include{1} = 'STI 014';
try
picks = fiff_pick_types(raw.info,want_meg,want_eeg,want_stim,include,raw.info.bads);
catch
%
% Failure: Try MEG + STI101 + STI201 + STI301 - bad channels instead
%
include{1} = 'STI101';
include{2} = 'STI201';
include{3} = 'STI301';
try
picks = fiff_pick_types(raw.info,want_meg,want_eeg,want_stim,include,raw.info.bads);
catch
error(me,'%s (channel list may need modification)',mne_omit_first_line(lasterr));
end
end
%
[outfid,cals] = fiff_start_writing_raw(outfile,raw.info,picks);
%
% Set up the reading parameters
%
from = raw.first_samp;
to = raw.last_samp;
quantum_sec = 10;
quantum = ceil(quantum_sec*raw.info.sfreq);
%
% To read the whole file at once set
%
%quantum = to - from + 1;
%
%
% Read and write all the data
%
first_buffer = true;
for first = from:quantum:to
last = first+quantum-1;
if last > to
last = to;
end
try
[ data, times ] = fiff_read_raw_segment(raw,first,last,picks);
catch
fclose(raw.fid);
fclose(outfid);
error(me,'%s',mne_omit_first_line(lasterr));
end
%
% You can add your own miracle here
%
fprintf(1,'Writing...');
if first_buffer
if first > 0
fiff_write_int(outfid,FIFF.FIFF_FIRST_SAMPLE,first);
end
first_buffer = false;
end
fiff_write_raw_buffer(outfid,data,cals);
fprintf(1,'[done]\n');
end
fiff_finish_writing_raw(outfid);
fclose(raw.fid);
c4science · Help