Method and system for determining azimuth of marsquake/moonquake event

A method for determining an azimuth of a marsquake or moonquake event, comprising:

• collecting original waveform data of a marsquake or moonquake by a plurality of seismometers in a single seismic station, wherein when the single seismic station is situated on Mars, the original waveform data collected by the plurality of the seismometers in the single seismic station is original waveform data of a marsquake;
• and when the single seismic station is situated on Moon, the original waveform data collected by the plurality of the seismometers in the single seismic station is original waveform data of a moonquake;
• acquiring, by a computer, the original waveform data of the marsquake or moonquake from the plurality of the seismometers in the single seismic station;
• filtering, by the computer, the original waveform data to obtain filtered waveform data;
• intercepting, by the computer, filtered waveform data of a set time length to obtain first arrival waveform data of a P wave, wherein a starting time of the set time length is an arrival time of the P wave;
• performing, by the computer, data analysis on the first arrival waveform data of the P wave, and determining the first arrival waveform data of the P wave of an east (E) component obtained after the data analysis as X-column data and the first arrival waveform data of the P wave of a north (N) component obtained after the data analysis as Y-column data;
• and
• performing, by the computer, optimal fitting on the X-column data and the Y-column data using a direct least square ellipse fitting algorithm to obtain a first arrival elliptical locus of the P wave obtained after the optimal fitting, and determining an azimuth of the marsquake or moonquake relative to the single seismic station according to the first arrival elliptical locus of the P wave obtained after the optimal fitting.