How to take the fourth derivative of a reflectance graph

Question

Erin Browne 2019-3-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/451331-how-to-take-the-fourth-derivative-of-a-reflectance-graph

评论： Torsten 2019-3-27

Hi I am looking to calculate and plot the fourth derivtative of a reflectance plot for coral spectra - commonly done in papers by the savistky and golay, 1964 method.

I have used the smoothing function on my coral mean spectra (x3) (1st image)

S2 = smooth(S1);

plot(wv1, S2);

xlim([400 750])

Then went to take my 1st derivative (2nd image)

This turned out ok using this coding:

dy = diff(S2)/diff(wv1);

plot(wv1(2:end),dy)

xlim([400 750])

But I believe this is the wrong method for taking the derivative of the reflectance spectra as I loose data going fromm 1446x1 to 1445x1 each time i derive i.e.

taking the second derivative like so:

d2ydx2 = diff(dy)./diff(wv1);

plot(wv1(1:1444),d2ydx2)

xlim([400 750])

producing a graph like image 3 - which when it comes to the 4th i have lost nearly all the data (image 4)

Does anyone know where I am going wrong, as this must be much more straight forward given it has been used so frequently.

Thanks

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Torsten 2019-3-21

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/451331-how-to-take-the-fourth-derivative-of-a-reflectance-graph#answer_366534

If you have the signal processing toolbox licenced, use "sgolay" and "sgolayfilt".

Or download from the file exchange:

https://de.mathworks.com/matlabcentral/fileexchange/30299-savitzky-golay-smooth-differentiation-filters-and-filter-application

4 个评论
显示 2更早的评论隐藏 2更早的评论

Erin Browne 2019-3-26

thanks for the help!

I've found that this method isn't working as my data is in 1nm and due to the fact it is going up and down so close together when taking any avergaing method it loses the quite obvious peaks, if you understand what im trying to say.

I believe binning into coarser resolutions may help?