Help with ICA with the WAY-EEG-GAL dataset

Hello all I hope you are doing great.

I am trying to remove the artifacts (EOG, EMG) from this EEG dataset

https://www.nature.com/articles/sdata201447

It has EEG recordings of 12 Subjects, 9 Sessions for each, each session containing around 29-34 trials.

I filtered the data with a FIR zero phase filter between 0.1-40 Hz and re-referenced to the average, and afterwards applied ICA for each session and used the Fp1, and Fp2 channels to remove correlated components.

eog_indices, _ = ica.find_bads_eog(eeg_data1, ch_name=[“Fp1”, “Fp2”])
ica.exclude = eog_indices

is this a proper way of doing it? this is my first time working with EEG data. Any help would be appreciated

Thank you all