如何使用Python计算半高宽?

2024-09-21 07:52:30 发布

您现在位置:Python中文网/ 问答频道 /正文

我正在开发一个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))

Tags: selfosnpmedminx1x2zero

热门问题