Due to the complex processes of earthquake preparation, the observations and studies associated with earthquakes have attracted the attention of geophysicists for many years. The CSES was successfully launched on 2 February 2018. This satellite can provide global data of the electromagnetic field, plasma, and energetic particles in the ionosphere to monitor and study the ionospheric perturbations associated with earthquakes. Focusing on the characteristics of CSES, a spatial analysis method was proposed to extract the disturbances of electron density prior to earthquakes. Taking Indonesia Mw6.9 earthquake that occurred on 5 August 2018 as an example, the spatial method was illustrated and verified by another analysis method also using the data of electron density and GPS TEC data with the same analysis method. Based on the electron density of CSES for more than 2 years, this method was applied to carry out the statistical study prior to Mw ≥ 6.0 global earthquakes using the superposed epoch and space approach (SESA) method. It was found that 1) relative to the epicenters, seismo-ionospheric disturbances are more obvious in the equator direction than those in the polar direction; 2) the anomalies within 300 km distance from the epicenter are significant 11, 3, and 2 days prior to Mw ≥ 6.0 earthquakes; 3) the influence region of perturbances associated with earthquakes enlarges with the magnitude increase, and the stronger magnitude is the earlier disturbance appears. These statistical characteristics were not detected for the random earthquakes. Comparing the statistical result with the simulation output, the electric field pathway could be considered as the main channel of lithosphere–atmosphere–ionosphere coupling.