Area between baseline and data peak
显示 更早的评论
Hi,
I am trying to determine the area under a peak curve, see image. The baseline illustrated in image is not zero, and I would like to subtract this from the area under the peak. The data is from a picture and been summed to obtain this profile. I want to determine the intensity by summing the area under the peak.
1. Is there a way to determine the baseline (y- value) automatically by code.
2. How do I determine the area minus the baseline? I use sum to determine the area under the peak curve.
采纳的回答
Star Strider
2014-9-22
编辑：Star Strider
2014-9-22
Here is how I would do it:
% Create Data
x = 0:480;
y = 1E4*(0.5+rand(size(x)));
y(400:450) = 1.3E+5*exp(-(x(400:450)-425).^2/50)+y(400:450);
% Generate Statistics & Identify Peak
ysts = [mean(y) median(y); mean(y)-1.96*std(y)/sqrt(length(y)) mean(y)+1.96*std(y)/sqrt(length(y))];
[ypk, pki] = max(y);
pkiv = find(y >= ysts(2,2)); % Find Peak Indices
Iy = cumtrapz(x, y-ysts(1,1)); % Integrate Entire Record
Iydb = [ones(pkiv(1),1) x(1:pkiv(1))']\Iy(1:pkiv(1))';
Iyd = Iy' - [ones(size(x))' x']*Iydb; % Detrend Using Baseline
Ipk = Iy(max(pkiv)) - Iy(min(pkiv)); % Find Peak Area
% Plot Data & Detrended Integral & Selected Statistics
figure(1)
plot(x, y)
hold on
plot(x, Iyd, '-g')
hold off
legend('Data', 'Detrended Integral', 'Location', 'NW')
text(100, 1E+6, sprintf('Peak Area = %23.15E', Ipk))
The idea is reasonably straightforward:
- Identify the peak as y-values greater than the upper 95% confidence interval for the mean of all the data;
- Integrate the entire record using cumtrapz;
- Do a linear regression on all the integrated data up to the first index of the peak;
- Use that regression to detrend the integrated data;
- Use the identified indices of the peak to calculate the approximate integral of the peak.
This approach may not generalise well if most of the baseline information is not at x-values less than the peak, because it depends on that information to detrend the integral.
The plot produces:
7 个评论
Muthu Annamalai
2014-9-22
I was mentally thinking do a cumsum and then diff the result to find the peak. Nice solution there!
Star Strider
2014-9-22
Thank you!
Star Strider
2014-9-23
Those two lines calculate the trend line ( ‘Iydb’ ) and detrend ( ‘Iyd’ ) the integral ‘Iy’. (The integral of a constant creates a linear trend. Subtracting that trend straightens out the integral line. It may not strictly be necessary, but it improves the accuracy of the peak integral calculation.)
I calculated a few statistics in ‘ysts’ and then referred to them rather than recalculating them. The value of ‘ysts(1,1)’ is the mean. I subtracted the mean from the signal before integrating it because it essentially zeros the baseline, improving the accuracy of the integral and peak area calculation. This does not affect the integral of the peak, since that is measured from the chosen threshold. (The other values in ‘ysts’ are the median and the upper and lower 95% confidence intervals on the mean. I later used the upper limit as the threshold to define the peak. I experimented with the other values in zeroing the initial baseline and defining the integral before deciding on the best option, that I posted.)
Lizan
2014-9-23
Star Strider
2014-9-23
My pleasure!
Image Analyst
2017-7-28
When you do \Iy you're doing matrix division. Did you really want element-by-element division? If so, use dot slash instead of slash.
Also ones() is returning a column vector, so make sure x is a row vector, since you transpose x to be a column vector and want to stitch it to the right of ones.
Also be away of automatic expansion of your ly vector since you're dividing a pkiv(1)-by-2 vector by a vector array.
更多回答（2 个）
Andrew Reibold
2014-9-22
编辑：Andrew Reibold
2014-9-22
0 个投票
The area under a curve is the definition of an integral. You will want to take the integral and subtract the baseline * baseline_width to find the area.
I'm not sure if its really scientific to determine the baseline just from raw data. Maybe this entire unique data set has some factor that shifts every point in some direction which makes you perceive the baseline somewhere that it isn't. I'm sure you could take the mean of non-peaked data or something, but normally the baseline is determined from the original function, historical data, or whatever you are stepping off from.
3 个评论
Image Analyst
2014-9-23
Lizan, you could speed along a solution to your problem by just posting your data.
Lizan
2014-10-21
Image Analyst
2014-9-22
0 个投票
You can find the peak (at least for this plot) very easily with max(). Then you can just "fall down" the peak on each side (with a for loop) and find the place (index) where the data is no longer decreasing. Then you can just fit a line between those two valley locations, say with polyfit() or whatever, to get the baseline.
2 个评论
Lizan
2014-10-20
Image Analyst
2014-10-20
signal(i) will be greater than or equal to signal(i-1).
类别
在 帮助中心 和 File Exchange 中查找有关 Parametric Modeling 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
