Epilepsy surgery is the most efficient treatment option for patients with refractory epilepsy. Before surgery, it is of utmost importance to accurately delineate the seizure onset zone (SOZ). Non-invasive EEG is the most used neuroimaging technique to diagnose epilepsy, but it is hard to localize the SOZ from EEG due to its low spatial resolution and because epilepsy is a network disease, with several brain regions becoming active during a seizure. In this work, we propose and validate an approach based on EEG source imaging (ESI) combined with functional connectivity analysis to overcome these problems. We considered both simulations and real data of patients. Ictal epochs of 204-channel EEG and subsets down to 32 channels were analyzed. ESI was done using realistic head models and LORETA was used as inverse technique. The connectivity pattern between the reconstructed sources was calculated, and the source with the highest number of outgoing connections was selected as SOZ. We compared this algorithm with a more straightforward approach, i.e. selecting the source with the highest power after ESI as the SOZ. We found that functional connectivity analysis estimated the SOZ consistently closer to the simulated EZ/RZ than localization based on maximal power. Performance, however, decreased when 128 electrodes or less were used, especially in the realistic data. The results show the added value of functional connectivity analysis for SOZ localization, when the EEG is obtained with a high-density setup. Next to this, the method can potentially be used as objective tool in clinical settings.