hiro/bandfit
1
0
Fork 0
mirror of https://github.com/vale981/bandfit synced 2025-03-05 09:31:42 -05:00
Code Issues Projects Releases Packages Wiki Activity Actions

improve continuity slightly

Browse source
This commit is contained in:
Valentin Boettcher 2023-03-21 10:24:56 -04:00
parent ad106facea
commit 4cbaf6f5cd
1 changed files with 17 additions and 4 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

19
bandfit/bandfit.py
View file

@ -170,14 +170,16 @@ def refine_band_fit(k, e0, data):
pass
def detect_bands_fixed_k(k, data, γ, last_separation=0, min_height=0.5):
def detect_bands_fixed_k(
k, data, γ, last_separation=0, min_height=0.5, separation_continuity=1 / 2
):
col = data[:, k].copy()
col -= col.min()
col /= col.max()
e_axis = np.arange(col.size)
guesses = np.array(sc.signal.find_peaks(col, distance=2, height=min_height)[0])
guesses, props = sc.signal.find_peaks(col, distance=2, height=min_height)
means = guesses[:, None] + guesses[None, :]
guess_i_1, guess_i_2 = np.unravel_index(
@ -186,9 +188,20 @@ def detect_bands_fixed_k(k, data, γ, last_separation=0, min_height=0.5):
guess_1, guess_2 = np.sort((guesses[guess_i_1], guesses[guess_i_2]))
if last_separation > 0 and abs(guess_2 - guess_1) < last_separation / 4:
if (
last_separation > 0
and abs((abs(guess_2 - guess_1) - last_separation)) / last_separation
> separation_continuity
):
guess_idx = (guesses[:, None] == [guess_1, guess_2]).argmax(axis=0)
heights = props["peak_heights"][guess_idx]
if heights[0] > heights[1]:
guess_2 = col.size - guess_1
else:
guess_1 = col.size - guess_2
col /= col[guess_1]
(e_1, e_2, γ, _), cov = sc.optimize.curve_fit(