ID identification
limited feature extraction / anomaly detection
There are couple of ways to identify the ID.
Variance method (Liu et al. 2022) : Large variance in the magnetic field compared with neighboring intervals (see notebook)
Partial variance increment (PVI) method :
- Vasko et al. (2022)
B-criterion (Burlaga and Ness 1969) : a directional change of the magnetic field larger than 30° during 60 s
TS-criterion (Tsurutani and Smith 1979) : \(|ΔB|/|B| \geq 0.5\) within 3 minutes
Traditional methods (B-criterion and TS-criterion) rely on magnetic field variations with a certain time lag. B-criterion has, as its main condition. In their methods, the IDs below the thresholds are artificially abandoned. Therefore, identification criteria may affect the statistical results, and there is likely to be a discrepancy between the findings via B-criterion and TS- criterion.
02-Jun-24 09:00:27: UserWarning: Traceback (most recent call last):
File "/Users/zijin/projects/ids_finder/.pixi/envs/default/lib/python3.11/site-packages/pdpipe/__init__.py", line 85, in <module>
from . import skintegrate
File "/Users/zijin/projects/ids_finder/.pixi/envs/default/lib/python3.11/site-packages/pdpipe/skintegrate.py", line 20, in <module>
from sklearn.base import BaseEstimator
ModuleNotFoundError: No module named 'sklearn'
02-Jun-24 09:00:27: UserWarning: pdpipe: Scikit-learn or skutil import failed. Scikit-learn-dependent pipeline stages will not be loaded.
02-Jun-24 09:00:27: UserWarning: Traceback (most recent call last):
File "/Users/zijin/projects/ids_finder/.pixi/envs/default/lib/python3.11/site-packages/pdpipe/__init__.py", line 105, in <module>
from . import nltk_stages
File "/Users/zijin/projects/ids_finder/.pixi/envs/default/lib/python3.11/site-packages/pdpipe/nltk_stages.py", line 19, in <module>
import nltk
ModuleNotFoundError: No module named 'nltk'
02-Jun-24 09:00:27: UserWarning: pdpipe: nltk import failed. nltk-dependent pipeline stages will not be loaded.
pl_format_time
pl_format_time (df:polars.lazyframe.frame.LazyFrame, tau:datetime.timedelta)
Pipelines
detect_events
detect_events (data:polars.lazyframe.frame.LazyFrame, tau:datetime.timedelta, ts:datetime.timedelta, bcols, sparse_num=None, method='liu', **kwargs)