Method and device for the automatic identification of phase-locked high-frequency oscillations (PLHFO) to localize epileptogenic brain for neurosurgical intervention, including filtering brain signals into low frequency and high frequency oscillation (HFO) data streams. Applying ICA to the HFO data stream, transforming the data streams to produce an HFO instantaneous amplitude (HFOIA) and a low-frequency instantaneous phase (LFIP) data stream. Transforming the normalized HFOIA to produce an instantaneous phase of the normalized HFOIA. Determining a continuous or discrete PLHFO calculation that measures cross frequency coupling between the instantaneous phase of the low frequency data stream, and the instantaneous amplitude of the HFO data stream based at least in part on LFIP, raw or normalized HFOIA, and may include the instantaneous phase of normalized or raw HFOIA. Determining that at least a portion of the electrical signals from the brain display PLHFO if the PLHFO calculation is above a statistical threshold.
