Hello,
the way this is usually treated is that you combine the channels by subtraction, such that:
VEOG = VEOG_lower - VEOG_upper
HEOG = HEOG_left - HEOG_right
You would roughly do something like the following:
veog_ch_name = 'VEOG' # the new channel that will be created
mne.set_bipolar_reference(
inst=raw,
anode='VEOG_lower',
cathode='VEOG_upper',
ch_name=veog_ch_name,
drop_refs=True, # drop anode and cathode from the data
copy=False # modify in-place
)
If not both input channels were already of type eog
, you need to set the type of the created channel accordingly, too:
raw.set_channel_types({veog_ch_name: 'eog'})