rtsBaselineCorrect.m
1.33 KB
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
48
49
50
51
52
53
54
55
56
57
58
59
60
61
function R = rtsBaselineCorrect(A, b, dim)
%This function baselines an n-dimensional array given a series of points b
%B = baseline-corrected data
%A = raw data
%b = baseline points in the format [b0 b1 ... bn]
%dim = dimension along which the baseline correction will occur
%allocate space for the sub-component indices
sz = size(A);
inds = repmat({1},1,ndims(A));
%allocate space for the result
R = zeros(size(A));
%for each baseline point
for bp = 1:length(b) - 1
%extract the data at the first baseline point
n = b(bp);
inds = getIndices(A, n, dim);
BP1 = A(inds{:});
%extract the data at the second baseline point
n = b(bp+1);
inds = getIndices(A, n, dim);
BP2 = A(inds{:});
%for each point in between the baseline points
for p = b(bp):b(bp+1)
%compute the weighting of the linear function between BP1 and BP2
a2 = (p - b(bp))/(b(bp+1) - b(bp));
a1 = 1.0 - a2;
n = p;
inds = getIndices(A, n, dim);
baseFunc = a1 * BP1 + a2 * BP2;
R(inds{:}) = A(inds{:}) - baseFunc;
end
end
function I = getIndices(A, n, dim)
%this function returns indices for all elements in S at
%point n along dimension dim
sz = size(A);
I = repmat({1},1,ndims(A));
for d = 1:ndims(A)
if d ~= dim
I{d} = 1:sz(d);
else
I{d} = n;
end
end