Total electron content (TEC) changes as an earthquake develops, so the seismic anomaly information can be analyzed through TEC data. However, TEC data is more sensitive to solar activity than seismic events, so previous studies can only choose the time period with relatively quiet solar ac- tivity. Affected by the magnetic storm, TEC data cannot normally reflect the seismic anomaly infor- mation in Dali Prefecture, Yunnan Province in 2021. Therefore, it is particularly important to sup- press the magnetic storm interference. Based on the sliding quartile method, a magnetic storm suppres- sion method for TEC data is proposed using the ARMA model. By refitting the TEC data before and after the magnetic storm instead of using the original data, the influence of magnetic storm on seismic information can be reduced to a certain extent. Then, a model of seismic features and their correlation is proposed to analyze the correlation between TEC anomalies and seismic events. Based on the global ionospheric TEC data provided by the IGS data analysis center, the 2021 Ms6.4 earthquake in Yanghe County, Dali Prefecture, Yunnan Province and the Ms7.4 earthquake in Maduo County, Guoluo Pre- fecture, Qinghai Province are simulated and analyzed respectively. The traditional sliding quartile method and sliding window method are used respectively, and the features and event correlation re- sults before and after magnetic storm suppression are compared through the proposed correlation mod- el. The results show that the correlation between earthquakes and features after removing magnetic storms increases by about 13~18%.