A nonlinear template-based approach for BCG artifact removal in EEG-fMRI recordings at high fields