使用python进行半高宽计算
使用python进行半高宽计算
提问于 2018-03-05 05:38:51
我正在尝试使用python计算光谱的半高宽。光谱描述(我说的是物理方面的)对我来说有点复杂，我不能用一些简单的高斯或洛伦兹分布来拟合数据。
到目前为止，我设法处理了数据的插值，并通过半极大值绘制了一条平行于X轴的直线。
怎样才能找到山峰两侧两条线交点的坐标呢？
我知道如果我把光标放在这些点上，它会给我坐标，但我想让这个过程自动化，让它变得更加用户友好。我该怎么做呢？
回答 2
Stack Overflow用户
回答已采纳
发布于 2018-03-05 07:28:29
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)。
Stack Overflow用户
发布于 2020-08-28 08:28:35
除了前面的答案之外，如果基线不是0，那么((max- min )/2) +min。这就是我为解决我的问题所做的。Tks。
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/49100778
复制相关文章
相似问题
腾讯云开发者
