2022-04-14 13:51:53 +02:00
|
|
|
import pandas as pd
|
2022-04-20 12:44:51 +02:00
|
|
|
import numpy as np
|
2022-04-14 13:51:53 +02:00
|
|
|
|
2022-04-19 15:24:46 +02:00
|
|
|
from cr_features.helper_functions import convert_ibi_to2d_time, hrv_features, hrv_freq_features
|
2022-04-20 12:44:51 +02:00
|
|
|
from cr_features.hrv import extract_hrv_features_2d_wrapper
|
2022-04-25 15:07:03 +02:00
|
|
|
from cr_features_helper_methods import get_sample_rate, extract_second_order_features
|
2022-04-14 13:51:53 +02:00
|
|
|
|
2022-04-20 12:44:51 +02:00
|
|
|
import math
|
2022-04-14 13:51:53 +02:00
|
|
|
import sys
|
|
|
|
|
2022-04-19 15:24:46 +02:00
|
|
|
pd.set_option('display.max_rows', 1000)
|
2022-04-20 12:44:51 +02:00
|
|
|
pd.set_option('display.max_columns', None)
|
2022-04-14 13:51:53 +02:00
|
|
|
|
2022-04-19 15:24:46 +02:00
|
|
|
|
|
|
|
def extract_ibi_features_from_intraday_data(ibi_intraday_data, features, window_length, time_segment, filter_data_by_segment):
|
2022-04-14 13:51:53 +02:00
|
|
|
ibi_intraday_features = pd.DataFrame(columns=["local_segment"] + features)
|
|
|
|
|
|
|
|
if not ibi_intraday_data.empty:
|
2022-04-19 15:24:46 +02:00
|
|
|
sample_rate = get_sample_rate(ibi_intraday_data)
|
2022-04-14 13:51:53 +02:00
|
|
|
|
|
|
|
ibi_intraday_data = filter_data_by_segment(ibi_intraday_data, time_segment)
|
|
|
|
|
|
|
|
if not ibi_intraday_data.empty:
|
|
|
|
|
|
|
|
ibi_intraday_features = pd.DataFrame()
|
2022-04-25 15:07:03 +02:00
|
|
|
|
2022-04-14 13:51:53 +02:00
|
|
|
# apply methods from calculate features module
|
2022-04-19 15:24:46 +02:00
|
|
|
if window_length is None:
|
|
|
|
ibi_intraday_features = \
|
|
|
|
ibi_intraday_data.groupby('local_segment').apply(\
|
2022-04-20 12:44:51 +02:00
|
|
|
lambda x:
|
2022-04-19 15:24:46 +02:00
|
|
|
extract_hrv_features_2d_wrapper(
|
2022-04-20 12:44:51 +02:00
|
|
|
signal_2D = \
|
|
|
|
convert_ibi_to2d_time(x[['timings', 'inter_beat_interval']], math.ceil(x['timings'].iloc[-1]))[0],
|
|
|
|
ibi_timings = \
|
|
|
|
convert_ibi_to2d_time(x[['timings', 'inter_beat_interval']], math.ceil(x['timings'].iloc[-1]))[1],
|
|
|
|
sampling=None, hampel_fiter=False, median_filter=False, mod_z_score_filter=True, feature_names=features))
|
2022-04-19 15:24:46 +02:00
|
|
|
else:
|
|
|
|
ibi_intraday_features = \
|
|
|
|
ibi_intraday_data.groupby('local_segment').apply(\
|
2022-04-20 12:44:51 +02:00
|
|
|
lambda x:
|
2022-04-19 15:24:46 +02:00
|
|
|
extract_hrv_features_2d_wrapper(
|
2022-04-20 12:44:51 +02:00
|
|
|
signal_2D = convert_ibi_to2d_time(x[['timings', 'inter_beat_interval']], window_length)[0],
|
|
|
|
ibi_timings = convert_ibi_to2d_time(x[['timings', 'inter_beat_interval']], window_length)[1],
|
|
|
|
sampling=None, hampel_fiter=False, median_filter=False, mod_z_score_filter=True, feature_names=features))
|
2022-04-14 13:51:53 +02:00
|
|
|
|
|
|
|
ibi_intraday_features.reset_index(inplace=True)
|
|
|
|
|
|
|
|
return ibi_intraday_features
|
|
|
|
|
|
|
|
|
|
|
|
def cr_features(sensor_data_files, time_segment, provider, filter_data_by_segment, *args, **kwargs):
|
|
|
|
ibi_intraday_data = pd.read_csv(sensor_data_files["sensor_data"])
|
|
|
|
|
|
|
|
requested_intraday_features = provider["FEATURES"]
|
|
|
|
|
|
|
|
calc_windows = kwargs.get('calc_windows', False)
|
|
|
|
|
|
|
|
if provider["WINDOWS"]["COMPUTE"] and calc_windows:
|
|
|
|
requested_window_length = provider["WINDOWS"]["WINDOW_LENGTH"]
|
|
|
|
else:
|
|
|
|
requested_window_length = None
|
|
|
|
|
|
|
|
# name of the features this function can compute
|
2022-04-19 15:24:46 +02:00
|
|
|
base_intraday_features_names = hrv_features + hrv_freq_features
|
2022-04-14 13:51:53 +02:00
|
|
|
# the subset of requested features this function can compute
|
|
|
|
intraday_features_to_compute = list(set(requested_intraday_features) & set(base_intraday_features_names))
|
|
|
|
|
|
|
|
# extract features from intraday data
|
2022-04-19 15:24:46 +02:00
|
|
|
ibi_intraday_features = extract_ibi_features_from_intraday_data(ibi_intraday_data, intraday_features_to_compute,
|
2022-04-14 13:51:53 +02:00
|
|
|
requested_window_length, time_segment, filter_data_by_segment)
|
2022-04-25 15:07:03 +02:00
|
|
|
if calc_windows:
|
|
|
|
so_features_names = provider["WINDOWS"]["SECOND_ORDER_FEATURES"]
|
|
|
|
ibi_second_order_features = extract_second_order_features(ibi_intraday_features, so_features_names)
|
|
|
|
return ibi_intraday_features, ibi_second_order_features
|
|
|
|
|
2022-04-14 13:51:53 +02:00
|
|
|
|
|
|
|
return ibi_intraday_features
|