我正在开发一个python系统,在这里我需要计算图形的半高宽,用户在其中表示X轴和Y轴的txt文件(我将轴保存在2个列表xList和yList2中)。我运行了两个图形,结果是:30.918nm(正确值:34.786nm)和76.601nm(正确值:74.710nm)。我正在使用这个解决方案:FWHM calculation using python。我的代码是:
import numpy as np
def lin_interp(x, y, i, half):
return x[i] + (x[i+1] - x[i]) * ((half - y[i]) / (y[i+1] - y[i]))
class FWHM:
def __init__(self, xList, yList2):
self.xList = xList
self.yList2 = yList2
#pega os valores máximo e mínimo do eixo Y e divide por 2
y_min = (min(yList2))
y_max = (max(yList2))
med_y = ((y_max-y_min)/2)
#cálculo para descobrir os pontos (x1 e x2) que estão à meia altura
signs = np.sign(np.add(yList2, -med_y))
zero_crossings = (signs[0:-2] != signs[1:-1])
zero_crossings_i = np.where(zero_crossings)[0]
x1 = lin_interp(xList, yList2, zero_crossings_i[1], med_y)
x2 = lin_interp(xList, yList2, zero_crossings_i[0], med_y)
#printa os pontos do eixo X e o FWHM
print(x1, x2, abs(x1-x2))
目前没有回答
相关问题 更多 >
编程相关推荐