tyler/flares
2
0
Fork 0
Code Issues 12 Pull Requests Actions Packages Projects Releases 10 Wiki Activity

fix to bandpass filter

Browse Source
This commit is contained in:
Tyler
2025-10-20 16:07:18 -07:00
parent 20b255321b
commit 953ea90c67
1 changed files with 3 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

5
flares.py
View File

@@ -1070,7 +1070,8 @@ def filter_the_data(raw_haemo):
average=True, xscale="log", color="r", show=False, amplitude=False average=True, xscale="log", color="r", show=False, amplitude=False
) )
raw_haemo = raw_haemo.filter(l_freq=None, h_freq=0.4, h_trans_bandwidth=0.2) #raw_haemo = raw_haemo.filter(l_freq=None, h_freq=0.4, h_trans_bandwidth=0.2)
raw_haemo = raw_haemo.filter(0.05, 0.7, h_trans_bandwidth=0.2, l_trans_bandwidth=0.02)
raw_haemo.compute_psd(fmax=2).plot( raw_haemo.compute_psd(fmax=2).plot(
average=True, xscale="log", axes=fig_filter.axes, color="g", amplitude=False, show=False average=True, xscale="log", axes=fig_filter.axes, color="g", amplitude=False, show=False
@@ -2827,7 +2828,7 @@ def calculate_dpf(file_path):
# order is hbo / hbr # order is hbo / hbr
with h5py.File(file_path, 'r') as f: with h5py.File(file_path, 'r') as f:
wavelengths = f['/nirs/probe/wavelengths'][:] wavelengths = f['/nirs/probe/wavelengths'][:]
logger.info("Wavelengths (nm):", wavelengths) logger.info(f"Wavelengths (nm): {wavelengths}")
wavelengths = sorted(wavelengths, reverse=True) wavelengths = sorted(wavelengths, reverse=True)
age = float(AGE) age = float(AGE)
logger.info(f"Their age was {AGE}") logger.info(f"Their age was {AGE}")