research/labsw
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

put logn ans short SG

Browse source
This commit is contained in:
François Boulogne 2025-04-10 15:33:52 +02:00
parent 4108e25bcb
commit d46897a920
2 changed files with 6 additions and 4 deletions
Download patch file Download diff file Expand all files Collapse all files

2
labsw/__init__.py
View file

@ -1 +1 @@
__version__ = '0.1.7.1'
__version__ = '0.1.7.2'

8
labsw/io.py
View file

@ -42,7 +42,7 @@ def load_data_RH_logger(filepath, every=1):
return df.iloc[::every]
def process_data_RH_logger(filepath, every, diff_period=1_000, SG_window=10_000):
def process_data_RH_logger(filepath, every, diff_period=1_000, SG_window_short=100, SG_window_long=10_000):
df = load_data_RH_logger(filepath, every=1)
df = df.drop(columns=['time'])
@ -56,10 +56,12 @@ def process_data_RH_logger(filepath, every, diff_period=1_000, SG_window=10_000)
# Derivative
delta = np.mean(df['time'].diff())
df['dMdt_SG'] = savgol_filter(df['M'], window_length=SG_window, polyorder=1, deriv=1, delta=delta)
df['dMdt_SG_short'] = savgol_filter(df['M'], window_length=SG_window_short, polyorder=1, deriv=1, delta=delta)
df['dMdt_SG_long'] = savgol_filter(df['M'], window_length=SG_window_long, polyorder=1, deriv=1, delta=delta)
df['dMdt_diff'] = df['M'].diff(periods=diff_period) / df['time'].diff(periods=diff_period)
df['dmdt_SG'] = savgol_filter(df['m'], window_length=SG_window, polyorder=1, deriv=1, delta=delta)
df['dmdt_SG_short'] = savgol_filter(df['m'], window_length=SG_window_short, polyorder=1, deriv=1, delta=delta)
df['dmdt_SG_long'] = savgol_filter(df['m'], window_length=SG_window_long, polyorder=1, deriv=1, delta=delta)
df['dmdt_diff'] = df['m'].diff(periods=diff_period) / df['time'].diff(periods=diff_period)
h5path = os.path.splitext(filepath)[0]