-
Notifications
You must be signed in to change notification settings - Fork 3
/
median_filter.asv
49 lines (34 loc) · 1.08 KB
/
median_filter.asv
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
%% load raw data, this time load all channels
set(0,'DefaultFigureWindowStyle','docked'); % fix matlab's figure positioning bug
% raw data available on
% usb sticks
%
datapath='C:\Users\Jakob Voigts\Desktop\pondababa_2021-09-18_09-33-02\Record Node 111\'
data_raw=[];
for ch=[1:4*6] % grab a few channels of raw data from one tetrode
fname=sprintf('100_%d.continuous',ch)
[data, timestamps, info]=load_open_ephys_data_faster(fullfile(datapath,fname));
data_raw(:,end+1) = data;
end;
disp('done loading');
%
data_raw=data_raw.*info.header.bitVolts;
fs = info.header.sampleRate;
%%data_raw=data_raw(1:30000,:); % cut away some data for faster testing
%% plot raw data
clf; hold on;
% data_raw(50000+[0:2],:)=-80; % add fake spike?
plotlim=200000;
spacing=200;
figure(1);
clf;
hold on;
median_data_window = median(data_raw([1:plotlim],:)')';
for i=1:size(data_raw,2)
plot(data_raw([1:plotlim],i)+(i-1)*spacing,'r');
plot(data_raw([1:plotlim],i)-median_data_window+(i-1)*spacing,'k');
end;
%% make median for entire recording
median_all_data=zeros(1,T);
stepsize=10000;
for t=1: