INTRODUCTION: EEG-correlated fMRI (EEG-fMRI) studies can reveal haemodynamic changes associated with Interictal Epileptic Discharges (IED). Methodological improvements are needed to increase sensitivity and specificity for localising the epileptogenic zone. We investigated whether the estimated EEG source activity improved models of the BOLD changes in EEG-fMRI data, compared to conventional << event-related >> designs based solely on the visual identification of IED. METHODS: Ten patients with pharmaco-resistant focal epilepsy underwent EEG-fMRI. EEG Source Imaging (ESI) was performed on intra-fMRI averaged IED to identify the irritative zone. The continuous activity of this estimated IED source (cESI) over the entire recording was used for fMRI analysis (cESI model). The maps of BOLD signal changes explained by cESI were compared to results of the conventional IED-related model. RESULTS: ESI was concordant with non-invasive data in 13/15 different types of IED. The cESI model explained significant additional BOLD variance in regions concordant with video-EEG, structural MRI or, when available, intracranial EEG in 10/15 IED. The cESI model allowed better detection of the BOLD cluster, concordant with intracranial EEG in 4/7 IED, compared to the IED model. In 4 IED types, cESI-related BOLD signal changes were diffuse with a pattern suggestive of contamination of the source signal by artefacts, notably incompletely corrected motion and pulse artefact. In one IED type, there was no significant BOLD change with either model. CONCLUSION: Continuous EEG source imaging can improve the modelling of BOLD changes related to interictal epileptic activity and this may enhance the localisation of the irritative zone.