Надежный алгоритм определения ширины пиков

Я спросил , как программно оценивать полосы спектраи @detlyпредложил использовать FWHM (полная ширина на половине максимума) для определения ширины вершин. Я поискал и обнаружил, что FWHM можно использовать для подгонки моделей (я практически не профан в этом!), особенно моделей Гаусса. В частности, 2,354 * сигма— это ширинадля модели Гаусса.

Я рассматриваю аппроксимацию Гаусса из-заналичия плохих пиков. На этой картинке есть 4 правильно сформированных пика. Затем есть «двойной» пик (хотя оба могут быть важны) и два расплывчатых пика. Они станут невыполнимой задачей для наивного FWHM.

Можете ли вы помочь с генерацией гауссовой аппроксимации (для расчета полуширины) пика в Scip/Numpy, учитывая его приблизительное местоположениеего координаты x? Если Guassian — плохой выбор, то какая-то другая схема.