This study uses a network that is comprised of 10 total intensity magnetometers to detect azimuthal propagation of seismo-magnetic emission waves during 26 earthquakes that occurred between July 2007 and December 2008 in Taiwan. The propagation azimuth and phase velocity of the seismo-magnetic waves are calculated using frequency wavenumber analysis at the ultra low frequency of 0.05 Hz every 30 min. We superimpose the derived azimuths within a moving window of 30 days as the monitored distributions, and the entire dataset as the background distribution. We also find the propagation azimuths of the seismomagnetic anomalies of each earthquake by subtracting the background from the monitored distributions. The results show that frequency wavenumber analysis can be applied to evaluate azimuthal propagation of seismo-magnetic emission waves using a scalar of geomagnetic total intensity fields. The success detection rate of seismo-magnetic anomalies increases from 62% of the 26 earthquakes to 77% using the surface magnetic anomalous reference tip (SMART) to substitute the epicenters. Meanwhile, the odds proportions between the azimuths of the seismomagnetic emission waves towards and away from SMART reveal the associated anomalous propagation.