python - 使用 python 计算 FWHM

Question

我正在尝试使用 python 计算光谱的 FWHM。光谱描述（我说的是物理学）对我来说有点复杂，我无法使用一些简单的高斯或洛伦齐分布来拟合数据。

到目前为止，我设法管理数据的插值并通过半最大值绘制一条平行于 X 轴的直线。

如何找到峰两侧两条线的交点坐标？

我知道如果我将光标放在这些点上，它会给我坐标，但我想自动化这个过程，以便它变得更加用户友好。我怎样才能做到这一点？

Answer 1

from matplotlib import pyplot as mp
import numpy as np

def peak(x, c):
    return np.exp(-np.power(x - c, 2) / 16.0)

def lin_interp(x, y, i, half):
    return x[i] + (x[i+1] - x[i]) * ((half - y[i]) / (y[i+1] - y[i]))

def half_max_x(x, y):
    half = max(y)/2.0
    signs = np.sign(np.add(y, -half))
    zero_crossings = (signs[0:-2] != signs[1:-1])
    zero_crossings_i = np.where(zero_crossings)[0]
    return [lin_interp(x, y, zero_crossings_i[0], half),
            lin_interp(x, y, zero_crossings_i[1], half)]

# make some fake data
x=np.linspace(0,20,21)
y=peak(x,10)

# find the two crossing points
hmx = half_max_x(x,y)

# print the answer
fwhm = hmx[1] - hmx[0]
print("FWHM:{:.3f}".format(fwhm))

# a convincing plot
half = max(y)/2.0
mp.plot(x,y)
mp.plot(hmx, [half, half])
mp.show()

(x, y)两点坐标为(hmx[0], half)和(hmx[1], half)。

Answer 2

除了前面的答案，如果基线不为 0，则 ((max-min)/2) + min。这就是我为解决我的问题所做的。谢了。