The application discloses a three-dimensional ionospheric disturbance extraction method, which comprises the following steps: acquiring double-frequency
observation data of a research area from a GNSS receiving
station; then, calculating STEC from the GNSS double-frequency
signal and converting the STEC into dSTEC to highlight the ionospheric disturbance
signal; secondly, constructing an
observation matrix containing prior ionospheric model information and utilizing an improved three-dimensional
tomography algorithm to iteratively correct and gradually approach the real
electron density distribution; thirdly, in the
iteration process, adopting inverse distance weighted interpolation to correct the pixels not being passed through, so as to ensure the overall smoothness; finally, inputting the dSTEC data into the improved three-dimensional
tomography model, obtaining the three-dimensional
ionospheric electron density distribution through iterative solution, and comparing the three-dimensional
ionospheric electron density distribution with the actual
observation data to verify the accuracy of the method. The method can effectively extract the three-dimensional structure of the ionospheric disturbance, improve the reconstruction accuracy and resolution, and provide an effective tool for studying the response of the
ionosphere to events such as earthquakes.