A flare occurrence position determination method and device, electronic equipment and computer medium

Abstract

The application discloses a flare occurrence position determination method and device, electronic equipment and computer medium. The method comprises the following steps: when it is determined that a flare occurs, a flare peak time, a full-disk image corresponding to the flare peak time and initial disk positions of each solar active region are obtained; according to the flare peak time and the initial disk positions of each solar active region, first disk positions of each solar active region corresponding to the flare peak time are determined; the first disk positions of each solar active region are subjected to coordinate system conversion, and pixel positions of each solar active region in the full-disk image are determined; according to the pixel positions of each solar active region and the full-disk image, first pixel values of regions in which each solar active region is located in the full-disk image are respectively calculated; and the first disk position of a solar active region corresponding to a maximum value is determined as the occurrence position of the flare. The application can accurately determine the occurrence position of the flare when the flare occurs.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of space environment products, and in particular to a method and device for determining the location of a solar flare, an electronic device and a computer medium. BACKGROUND

[0002] In current space weather forecasting, timely warning of solar major flares is an important content. Solar flares are one of the important manifestations of solar eruptions, and the enhanced electromagnetic radiation can cause sudden disturbances to the Earth's ionosphere facing the sun, thereby interfering with high-frequency radio communication and satellite communication. Major flares are usually accompanied by coronal mass ejections and solar proton events, and the larger the scale of the flare, the larger the scale of other solar eruption activities that accompany it, thereby triggering a series of different forms of Earth's space environment disturbance events such as geomagnetic storms, high-energy electron storms, and ionospheric storms. Therefore, timely and accurate solar flare warning is an important means to avoid or mitigate the hazards of space weather to human production and life.

[0003] Existing flare warning information generally includes the time of flare occurrence, the end time, the event intensity, etc., and lacks information about the location of the corresponding flare. However, the location information of the flare plays a decisive role in determining the occurrence of solar proton events and whether the accompanying coronal mass ejections can reach the Earth. Therefore, in the warning of flares, it is urgent to add the location information of the flare to the warning information. SUMMARY

[0004] The present application aims to overcome the defects of the prior art and proposes a method and device for determining the location of a solar flare, an electronic device and a computer medium.

[0005] In a first aspect, the present application provides a method for determining the location of a solar flare, the method comprising:

[0006] Upon determining the occurrence of a solar flare, obtaining the peak time of the flare, a full-disk image corresponding to the peak time of the flare, and the initial heliographic position of each solar active region;

[0007] According to the peak time of the flare and the initial heliographic position of each solar active region, determining the first heliographic position of each solar active region at the peak time of the flare;

[0008] Converting the coordinate system of the first heliographic position of each solar active region to determine the pixel position of each solar active region in the full-disk image;

[0009] According to the pixel position of each solar active region and the full-disk image, calculating the first pixel value of the region in the full-disk image where each solar active region is located;

[0010] Determine a maximum pixel value from the first pixel values, and determine the first solar surface position of the solar active region corresponding to the maximum pixel value as the flare occurrence position.

[0011] In some embodiments, the first solar surface position of each solar active region at the flare peak time is determined according to the flare peak time and the initial solar surface position of each solar active region, including:

[0012] Obtain a solar active region summary file and a publishing time of the solar active region summary file;

[0013] Determine the initial solar surface position of each solar active region in the solar active region summary file;

[0014] Determine the first solar surface position of each solar active region at the flare peak time according to the flare peak time, the publishing time, the initial solar surface position of each solar active region, and a preset solar surface differential rotation formula.

[0015] In some embodiments, the sum of pixel values of each solar active region in the region of the full solar surface image is calculated according to the pixel position and the full solar surface image, including:

[0016] For each solar active region, obtain the pixel value of a target pixel point in the region of the full solar surface image;

[0017] Calculate the sum of pixel values of target pixel points based on all target pixel points in the region of the full solar surface image of the solar active region, to obtain the first pixel value of the region of the full solar surface image of the solar active region.

[0018] In some embodiments, after the target pixel point in the region of the full solar surface image of the solar active region is obtained, the method further includes:

[0019] Determine a plurality of target pixel points according to the size of the pixel value of each pixel point, wherein the pixel value of the target pixel point is greater than a preset pixel threshold value;

[0020] The first pixel value of the region of the full solar surface image of the solar active region is calculated based on the sum of pixel values of all target pixel points in the region of the full solar surface image of the solar active region, including:

[0021] Calculate the sum of pixel values of target pixel points based on the target pixel points in the region of the full solar surface image of the solar active region, to obtain the first pixel value of the region of the full solar surface image of the solar active region.

[0022] In some embodiments, the method further comprises:

[0023] obtaining four consecutive 1-minute X-ray flux values;

[0024] when the four consecutive 1-minute X-ray flux values are increasing, and the fourth 1-minute X-ray flux value is greater than 1.4 times the first 1-minute X-ray flux value, determining that a flare occurs.

[0025] In some embodiments, the method further comprises:

[0026] obtaining the first solar surface position of the solar active region where the flare occurs;

[0027] determining the active region number of the solar active region where the flare occurs according to a solar active region summary file;

[0028] generating flare alert information according to the first solar surface position and the active region number of the solar active region where the flare occurs.

[0029] In some embodiments, the method further comprises:

[0030] determining a flare occurrence time, a flare peak time, a flare end time and a flare intensity according to four consecutive 1-minute X-ray flux values;

[0031] The first solar surface position and the active region number of the solar active region where the flare occurs, generating flare alert information, comprises:

[0032] generating the flare alert information according to the flare occurrence time, the flare peak time, the flare end time, the flare intensity, the first solar surface position and the active region number of the solar active region where the flare occurs.

[0033] In a second aspect, the present application provides a flare occurrence position determination device, the device comprising:

[0034] an obtaining module, configured to, when a flare is determined to occur, obtain a flare peak time, a full-disk image corresponding to the flare peak time and initial solar surface positions of each solar active region;

[0035] a first calculating module, configured to determine first solar surface positions of each solar active region corresponding to the flare peak time according to the flare peak time and the initial solar surface positions of each solar active region;

[0036] a converting module, configured to perform coordinate system conversion on the first solar surface positions of each solar active region to determine pixel positions of each solar active region in the full-disk image;

[0037] a second calculating module, configured to calculate a first pixel value of each of the solar active regions in a region of the full-disk image according to the pixel position of each of the solar active regions and the full-disk image;

[0038] a determining module, configured to determine a maximum value in the first pixel values, and determine the first solar position of the solar active region corresponding to the maximum value of the first pixel values as the flare occurrence position.

[0039] In a third aspect, the present application provides an electronic device, comprising: a processor, a memory, a communication interface and a communication bus, the processor, the memory and the communication interface complete the communication among each other through the communication bus;

[0040] The memory is used to store at least one executable instruction, and the executable instruction makes the processor execute the operations corresponding to the flare occurrence position determination method in any one of the first aspect.

[0041] In a fourth aspect, the present application provides a computer storage medium, which stores a computer program, and the program is executed by a processor to realize the flare occurrence position determination method in any one of the first aspect.

[0042] Compared with the prior art, the present application has the following advantages:

[0043] In the flare occurrence, the flare peak time and the initial solar position of each solar active region are used to accurately determine the first solar position of each solar active region in the flare occurrence, and the coordinate system conversion of the first solar position of the solar active region is used to accurately determine the pixel position of each solar active region in the full-disk image in the flare occurrence, and the first pixel value of each solar active region in the region of the full-disk image is accurately calculated by using the pixel position of each solar active region and the full-disk image, so that the activity degree of each solar active region can be determined by using the first pixel value, and the solar active region corresponding to the maximum value in the first pixel value is determined to accurately determine the most active solar active region in the flare peak time, so that the first solar position of the most active solar active region in the flare peak time is determined as the flare occurrence position, and the flare occurrence position can be accurately determined in the flare occurrence, thereby effectively avoiding the problem of lacking the flare occurrence position in the prior art. BRIEF DESCRIPTION OF DRAWINGS

[0044] Figure 1 is a flare intensity diagram provided by the present application;

[0045] Figure 2 The present application provides a flare occurrence position determination method flowchart.

[0046] Figure 3 is a flow chart of a method for determining a location of a flare provided by the present application;

[0047] Figure 4 is a schematic diagram of a flare alarm generation provided by the present application;

[0048] Figure 5 is a schematic diagram of a device for determining a location of a flare provided by the present application;

[0049] Figure 6 is a schematic diagram of an electronic device structure provided by the present application. DETAILED DESCRIPTION

[0050] The technical solutions of the present application will be described in detail below with reference to the accompanying drawings and examples.

[0051] In current space weather service prediction, timely warning of solar major flare is an important content. Solar flare is one of the important manifestations of solar eruption, and the enhanced electromagnetic radiation can cause sudden disturbance of the Earth's ionosphere towards the sun, and then cause interference to high frequency radio communication and satellite communication. Major flare is usually accompanied by coronal mass ejection and solar proton event, and the larger the size of the flare, the larger the scale of other solar eruption activities accompanied by it, and then a series of different forms of space environment disturbance events on the earth, such as geomagnetic storm, high energy electron storm, ionospheric storm, etc. Therefore, timely and accurate solar flare warning is an important means to avoid or reduce the harm of space weather to human production and life.

[0052] In professional space environment prediction institutions, such as the Space Environmental Prediction Center (SEPC) of the Chinese Academy of Sciences, solar flare warning has realized automatic warning of X-ray flare above M-class. The data used for warning comes from real-time soft X-ray flux data provided by GOES satellite, and the warning content is the occurrence time, end time, event intensity, etc. obtained from the flux data, but lacks the information of the location of the flare. However, the location information of the flare plays a decisive role in judging the occurrence of solar proton event and whether the coronal mass ejection can reach the earth. Therefore, in the warning of solar major flare, it is urgent to add the information of the location of the flare to facilitate us to judge more quickly whether the current solar eruption is effective to the earth.

[0053] Solar flare identification based on solar observation images is a means to solve this problem. In recent decades, several algorithms for automatically detecting flares from solar images have been proposed and implemented: one algorithm segments the solar observation image into a small number of macropixels and then compares the differences between the current image and a reference image to determine whether a flare has occurred; another algorithm focuses on the intensity increase in the original image region to determine the flare location; machine learning methods such as neural networks and support vector machines have also been applied to the automatic detection of flares in solar images. Most of these algorithms are designed solely for detecting flares on solar images, and the identified flare occurrence time and number differ significantly from the flare results monitored operationally using GOES satellite X-ray flux, making location matching difficult. Therefore, there is an urgent need to design a flare identification method that can provide comprehensive information such as flare occurrence time, intensity, and source region location to meet the flare warning requirements in operational space environment forecasting.

[0054] Solar flares are primarily characterized by enhanced electromagnetic radiation, including radiation with wavelengths shorter than [missing information]. From gamma rays and X-rays to radio waves with wavelengths reaching several kilometers, the enhancement is most significant in short-wave radiation (X-rays and ultraviolet light) and the radio band. Currently, the internationally widely used... The X-ray flux value is used to classify flares. For example... Figure 1 The diagram shown is a schematic representation of the flare intensity provided in this application embodiment, based on the order of X-ray peak flux observed by the GOES satellite (10). -8 ~10 -4 W / m 2 Flares are classified into five levels: A, B, C, M, and X. The numbers following each level indicate the specific value of the peak X-ray flux.

[0055] Specifically, such as Figure 1 As shown, C2 level indicates a peak soft X-ray flux of 2 × 10⁻⁶. -6 W / m 2 M5 indicates that the peak soft X-ray flux of the flare is 5 × 10⁻⁶. -5 W / m 2 .

[0056] To address the lack of location information in flare alarms, this application provides a method for identifying the location of a flare occurrence, the specific steps of which are as follows: Figure 2 The method includes:

[0057] S201, when a flare is determined to occur, acquire the flare peak time, the full solar surface image corresponding to the flare peak time, and the initial solar surface position of each solar active region.

[0058] In the embodiments of the present application, when it is determined that a flare occurs on the sun, the flare peak time, the full-disk image corresponding to the flare peak time, and the initial heliographic position of each solar active region are obtained. Specifically, the flare peak time can be obtained by using the common flare identification algorithm. The full-disk image can be provided by the Atmospheric Imaging Assembly (AIA) of the Solar Dynamics Observatory (SDO), and the time accuracy of the full-disk image is about 3 minutes. The full-disk image can be an observation image of the 171 angstrom, 193 angstrom, 211 angstrom, 304 angstrom, or the like, and is preferably an observation image of the 193 angstrom. The initial heliographic position of each solar active region can be obtained from the Solar Region Summary (SRS) published by the Space Weather Prediction Center (SWPC) at 00:30 every day. The SRS also includes the solar active region number.

[0059] Further, when the flare peak time and the plurality of full-disk images with a time accuracy of 3 minutes are obtained, the flare peak time is compared with the time of the plurality of full-disk images with a time accuracy of 3 minutes to determine the full-disk image corresponding to the flare peak time. Specifically, in the embodiments of the present application, the full-disk image corresponding to the flare peak time is the full-disk image closest to the flare peak time.

[0060] S202, according to the flare peak time and the initial heliographic position of each solar active region, the first heliographic position of each solar active region corresponding to the flare peak time is determined.

[0061] In the embodiments of the present application, after the above data is obtained, since the Solar Region Summary is published only once a day, when it is determined that a flare occurs, the current position of each solar active region corresponding to the flare peak time, i.e., the first heliographic position of each solar active region, needs to be determined according to the flare peak time and the initial heliographic position of each solar active region.

[0062] In one example, the heliographic position in the embodiments of the present application can be represented by a heliocentric spherical coordinate system.

[0063] S203, the first heliographic position of each solar active region is converted in the coordinate system to determine the pixel position of each solar active region in the full-disk image.

[0064] In the embodiments of the present application, after the first solar disk position of each solar active region at the flare peak time is determined, the first solar disk position of each solar active region is converted in the coordinate system to convert the first solar disk position in the heliocentric spherical coordinate system into a pixel position in the pixel coordinate system of the full-disk image.

[0065] Specifically, according to the information such as the heliocentric, the solar radius, the pixel size provided by the full-disk image in the 193 angstrom band, the first solar disk position of the solar active region in the heliocentric spherical coordinate is converted into a pixel position (Px, Py) in the image pixel coordinate of the full-disk image.

[0066] S204, according to the pixel position of each solar active region and the full-disk image, the first pixel value of the region where each solar active region is located in the full-disk image is calculated respectively.

[0067] In the embodiments of the present application, after the pixel position of each solar active region is obtained, the region where each solar active region is located in the full-disk image is determined, and then the first pixel value of the region where each solar active region is located in the full-disk image is calculated by using the pixel value of the pixel point in the region where the solar active region is located in the full-disk image.

[0068] Specifically, when the region where the solar active region is located in the full-disk image is determined, a rectangular frame is pulled on the full-disk image with the pixel position (Px, Py) of the solar active region as the center.

[0069] In one example, the size of the rectangular frame can be converted according to the field of view range. For example, the image resolution in the x direction under the solar disk projection coordinate provided in the full-disk image is CDELT1=2.4 [arcsec / pixel], and then the pixel number corresponding to the length of 6' of the solar active region is about 6*60 / 2.4=150 [pixel], and the pixel number corresponding to the width of 4' is 4*60 / 2.4=100 [pixel], so that the region where the solar active region is located in the full-disk image can be determined.

[0070] S205, the maximum value is determined in the plurality of first pixel values, and the first solar disk position of the solar active region corresponding to the maximum value is determined as the occurrence position of the flare.

[0071] In the embodiments of the present application, after the first pixel value of the region where each solar active region is located in the full-disk image is obtained, the maximum value of the first pixel value is determined according to the size of the first pixel value of the region where each solar active region is located in the full-disk image, so that the target solar active region corresponding to the maximum value of the first pixel value can be determined, and then the first solar disk position of the target solar active region is determined as the occurrence position of the flare.

[0072] Specifically, the pixel value of each pixel point in the full-disk image represents the radiation intensity of the solar atmosphere at the wavelength, and thus the size of the first pixel value of the region where the solar active region is located in the full-disk reflects the activity level of the solar active region. Further, the occurrence position of the flare can be determined by using the size of the first pixel value of the region where the solar active region is located in the full-disk.

[0073] The embodiment of the present application can accurately determine the first disk position of each solar active region at the time of flare occurrence by using the flare peak time and the initial disk position of each solar active region, further accurately determine the pixel position of each solar active region in the full-disk image by performing coordinate system conversion on the first disk position of the solar active region, further accurately calculate the first pixel value of the region where each solar active region is located in the full-disk image by using the pixel position of each solar active region and the full-disk image, so as to determine the activity level of each solar active region by using the plurality of first pixel values, further accurately determine the most active solar active region at the flare peak time by determining the solar active region corresponding to the maximum first pixel value in the plurality of first pixel values, and thus accurately determine the occurrence position of the flare at the time of flare occurrence by determining the first disk position of the most active solar active region at the flare peak time as the occurrence position of the flare, thereby effectively avoiding the problem of lacking the occurrence position of the flare in the prior art.

[0074] Further, the embodiment of the present application can determine the first disk position of each solar active region corresponding to the flare peak time according to the flare peak time and the initial disk position of each solar active region, which can include the following steps:

[0075] S2021, obtaining a solar active region summary file and a publishing time of the solar active region summary file.

[0076] In the embodiment of the present application, when determining the first disk position of each solar active region corresponding to the flare peak time, first, a solar active region summary file (SRS) published by the Space Weather Prediction Center (SWPC) is obtained, and the publishing time of the solar active region summary file is also obtained.

[0077] S2022, determining the initial disk position of each solar active region in the solar active region summary file.

[0078] In the embodiment of the present application, after obtaining the solar active region summary file, the solar active region summary file is parsed to determine the initial disk position of each solar active region. Specifically, the initial disk position of the embodiment of the present application is the disk position of each solar active region at the publishing time of the solar active region summary file.

[0079] S2023, determining a first position of each solar active region corresponding to the flare peak time according to the flare peak time, the release time, initial heliographic positions of the solar active regions, and a preset heliographic differential rotation formula.

[0080] In the embodiment of the present application, after the flare peak time and the release time are determined, the time difference is calculated, and then the first heliographic position of each solar active region corresponding to the flare peak time is calculated by using the time difference, the initial heliographic positions of the solar active regions, and the preset heliographic differential rotation formula.

[0081] In one example, the preset heliographic differential rotation formula in the embodiment of the present application is the Carrington rotation period formula.

[0082] The embodiment of the present application can accurately determine the first heliographic position of the solar active region corresponding to the flare peak time by using the flare peak time, the release time, the initial heliographic positions of the solar active regions, and the preset heliographic differential rotation formula, so that the occurrence position of the flare can be determined by using the first heliographic position of the solar active region, and the problem of lacking the occurrence position of the flare in the prior art is effectively avoided.

[0083] Further, in the embodiment of the present application, the first pixel value of the region where each solar active region is located in the full-disk image is calculated according to the pixel positions of the solar active regions and the full-disk image, which can include the following steps:

[0084] S2041, for each solar active region, obtaining pixel values of all pixel points in the region where the solar active region is located in the full-disk image.

[0085] In the embodiment of the present application, for each solar active region, after the region where the solar active region is located in the full-disk image is obtained, all pixel points in the region where each solar active region is located in the full-disk image are obtained respectively, so that the pixel value of each pixel point can be obtained.

[0086] S2042, calculating the first pixel value of the region where the solar active region is located in the full-disk image based on the sum of the pixel values of all pixel points in the region.

[0087] In the embodiment of the present application, the pixel values of all pixel points in the region where the solar active region is located in the full-disk image are summed to obtain the first pixel value of the region where each solar active region is located in the full-disk image.

[0088] Specifically, the first pixel value of the region where the solar activity region is located in the full-disk image in the embodiments of the present application is used to represent the radiation intensity of the solar activity region, the radiation intensity can represent the activity of the solar activity region, the size of the first pixel value is positively correlated with the radiation intensity, and the radiation intensity is positively correlated with the activity of the solar activity region.

[0089] The embodiments of the present application can accurately calculate the first pixel value of each solar activity through the full-disk image, so that the activity of the solar activity region can be determined by using the first pixel value, and then the position where the flare occurs can be determined by using the first pixel value.

[0090] Further, when calculating the first pixel value of the region where the solar activity region is located in the full-disk image, the following steps can be further included:

[0091] S2043, determining a plurality of target pixel points according to the size of the pixel value of each pixel point.

[0092] In the embodiments of the present application, when calculating the first pixel value of the region where the solar activity region is located in the full-disk image, the preset pixel threshold can be obtained first, and the value of each pixel point in the region where the solar activity region is located in the full-disk image is compared with the preset pixel threshold to screen the target pixel point. Specifically, the value of the target pixel point is greater than the preset pixel threshold, and each solar activity region corresponds to a different preset pixel threshold.

[0093] Specifically, when determining the preset pixel threshold corresponding to each solar activity region, the pixel values of each pixel point in the region where the solar activity region is located in the full-disk image are used to sort all the pixel points in the region where the solar activity region is located in the full-disk image, and further, a target position is selected in the sorting result, and the pixel value corresponding to the target position is determined as the preset pixel threshold. Further, the target position in the embodiments of the present application is the front 5%-25% in the arrangement result, and the preferred position is 5%.

[0094] In one example, if there are 10,000 pixel points in the region where the solar activity region is located in the full-disk image, the pixel points are arranged from large to small according to the corresponding pixel values, and the pixel value corresponding to the 500th pixel point can be determined as the preset pixel threshold.

[0095] Further, when S2042 is performed, the sum of the pixel values of the target pixel points in the region where the solar activity region is located in the full-disk image is calculated to obtain the first pixel value of the region where each solar activity region is located in the full-disk image.

[0096] The embodiment of the present application can select the pixel point that can best represent the activity of the solar active region by presetting the pixel threshold to filter the target pixel point, and can obtain the first pixel value of the region of the solar active region in the full-disk image by calculating the sum of the pixel values of the target pixel point, so that the most active solar active region when the flare occurs can be determined more accurately, and the first disk position of the most active solar active region is used to determine the flare occurrence position, thereby effectively avoiding the problem of lacking the flare occurrence position in the prior art.

[0097] Further, the flare occurrence position determination method of the embodiment of the present application can further include the following steps:

[0098] S206, obtaining the flow values of X-rays of four consecutive 1 minutes.

[0099] In the embodiment of the present application, when determining whether a flare occurs, the flow values of X-rays of four consecutive 1 minutes are first obtained. Specifically, the flow values of X-rays of four consecutive 1 minutes can be determined according to the flow data of 1-8 angstroms with 1-minute accuracy provided by Geostationary Operational Environmental Satellites (GOES).

[0100] S207, when the flow values of X-rays of four consecutive 1 minutes are increasing, and the flow value of X-rays of the fourth 1 minute is greater than 1.4 times the flow value of X-rays of the first 1 minute, it is determined that a flare occurs.

[0101] In the embodiment of the present application, after obtaining the flow values of X-rays of four consecutive 1 minutes, the flow values of X-rays of four consecutive 1 minutes are compared, and after determining that the flow values of X-rays of four consecutive 1 minutes are strictly increasing, the flow value of X-rays of the fourth 1 minute is compared with the flow value of X-rays of the first 1 minute, and when the flow value of X-rays of the fourth 1 minute is greater than 1.4 times the flow value of X-rays of the first 1 minute, it is determined that a flare occurs.

[0102] The embodiment of the present application can accurately determine whether a flare occurs by using the flow values of X-rays of four consecutive 1 minutes, and then when a flare occurs, the flare occurrence position can be determined by using the flare peak time, the full-disk image corresponding to the flare peak time, and the initial disk position of each solar active region, thereby effectively avoiding the problem of lacking the flare occurrence position in the prior art.

[0103] In one example, the four continuous 1-minute X-ray flow values obtained in the embodiment S206 of the application can also be used to determine the flare occurrence time, flare peak time, flare end time, flare intensity, and other information. Specifically, the flare occurrence time is the time when the first 1-minute X-ray flow value is determined, the flare peak time is the time when the 1-minute X-ray flow value reaches the maximum value, the flare end time is the time when the 1-minute X-ray flow value decays to the middle point between the peak value and the background level before the flare occurrence, and the flare intensity is determined according to the X-ray flow value at the flare peak time.

[0104] Further, the flare occurrence position determination method of the embodiment of the application can further include the following steps:

[0105] S208, obtaining the first heliographic position of the solar active region where the flare occurs.

[0106] In the embodiment of the application, after obtaining the first heliographic position of the solar active region where the flare occurs, the position is obtained.

[0107] S209, determining the active region number of the solar active region where the flare occurs according to the solar active region summary file.

[0108] In the embodiment of the application, after obtaining the first heliographic position of the solar active region where the flare occurs, the active region number of the solar active region where the flare occurs is determined in the solar active region summary file based on the initial heliographic position of the solar active region.

[0109] S210, generating the flare warning information according to the first heliographic position and the active region number of the solar active region where the flare occurs.

[0110] In the embodiment of the application, after obtaining the first heliographic position and the active region number of the solar active region where the flare occurs, the first heliographic position and the active region number of the solar active region where the flare occurs are integrated to generate the flare warning information, thereby avoiding the limitation of the lack of flare occurrence position in the prior art.

[0111] In one example, when generating the flare warning information, the flare occurrence time, flare peak time, flare end time, and flare intensity can also be obtained, so as to use the flare occurrence time, flare peak time, flare end time, flare intensity, first heliographic position of the solar active region where the flare occurs, and the active region number to generate the flare warning information, thereby improving the completeness and practicability of the flare warning information.

[0112] Further, the embodiment of the application provides a flare occurrence position determination method, and the specific steps are as shown in Figure 3 .

[0113] Please refer to Figures 3-4As shown, the method comprises obtaining four consecutive 1-minute precision X-ray flux values at 1-8 angstrom provided by Geostationary Operational Environmental Satellites (GOES), obtaining about 3-minute precision full-disk images of the 193 angstrom band provided by the Atmospheric Imaging Assembly (AIA) of the Solar Dynamics Observatory (SDO), and obtaining a Solar Region Summary (SRS) document published by the Space Weather Prediction Center (SWPC), determining that a flare occurs when the four consecutive 1-minute X-ray flux values increase, and the fourth 1-minute X-ray flux value is greater than 1.4 times the first 1-minute X-ray flux value, the observation time corresponding to the first 1-minute X-ray flux is the flare occurrence time, and then using the 1-minute X-ray flux data to determine the flare peak time, the flare end time and the flare intensity.

[0114] Meanwhile, the flare peak time is compared with the time of a plurality of 3-minute precision full-disk images to determine the full-disk image closest to the observation of the flare peak time, and the initial solar disk position of each solar active region and the publication time of the Solar Region Summary document are obtained from the Solar Region Summary document. Further, after determining the flare peak time and the publication time, the time difference is calculated, and then the first solar disk position of each solar active region at the flare peak time is calculated using the time difference, the initial solar disk position of each solar active region and a preset solar disk differential rotation formula.

[0115] After determining the first solar disk position of each solar active region at the flare peak time, the coordinate system of the solar active region at the first solar disk position is converted to determine the pixel position of each solar active region in the pixel coordinate system of the full-disk image.

[0116] After obtaining the pixel position of the solar active region in the pixel coordinate system of the full-disk image, a rectangular frame is pulled around the pixel position to determine the area where each solar active region is located in the full-disk image.

[0117] After determining the region of each solar activity region in the full-disk image, the value of each pixel point in the region where the solar activity region is located in the full-disk image is compared with a preset pixel threshold value, a target pixel point in the region where the solar activity region is located in the full-disk image is determined, the sum of pixel values of the target pixel point in the region where the solar activity region is located in the full-disk image is calculated, and a first pixel value of the region where each solar activity region is located in the full-disk image is obtained.

[0118] After obtaining the first pixel value of the region where each solar activity region is located in the full-disk image, a first pixel value maximum is determined according to the size of the first pixel value of the region where each solar activity region is located in the full-disk image, a target solar activity region corresponding to the first pixel value maximum is determined, and the first disk position of the target solar activity region is further determined as the occurrence position of the flare.

[0119] The flare occurrence time, the flare peak time, the flare end time, and the flare intensity are obtained, and the flare warning information is generated by using the flare occurrence time, the flare peak time, the flare end time, the flare intensity, the first disk position of the solar activity region where the flare occurs, and the activity region number.

[0120] Further, the embodiment of the present application provides a flare occurrence position determination device, as shown in the figure, which comprises: Figure 5

[0121] The acquisition module 501 is configured to, when a flare occurs, acquire the flare peak time, the full-disk image corresponding to the flare peak time, and the initial disk position of each solar activity region.

[0122] The first calculation module 502 is configured to determine the first disk position of each solar activity region at the flare peak time according to the flare peak time and the initial disk position of each solar activity region.

[0123] The conversion module 503 is configured to perform coordinate system conversion on the first disk position of each solar activity region to determine the pixel position of each solar activity region in the full-disk image.

[0124] The second calculation module 504 is configured to calculate the first pixel value of the region where each solar activity region is located in the full-disk image according to the pixel position of each solar activity region and the full-disk image.

[0125] The determination module 505 is configured to determine the maximum value in the plurality of first pixel values, and determine the first disk position of the solar activity region corresponding to the first pixel value maximum as the occurrence position of the flare.

[0126] ​The flare occurrence position determination apparatus of the embodiment is used to implement the flare occurrence position determination method in the foregoing method embodiments, and has the beneficial effects of the corresponding method embodiments, which will not be repeated here. In addition, the functions of each module in the flare occurrence position determination apparatus of the embodiment can be implemented by referring to the descriptions of the corresponding parts in the foregoing method embodiments, which will not be repeated here either.

[0127] Referring to Figure 6 , a structural schematic diagram of an electronic device according to an embodiment of the present application is shown, and the specific implementation of the electronic device is not limited in the embodiments of the present application.

[0128] As Figure 6 shown, the electronic device can include a processor 602, a communications interface 604, a memory 606, and a communications bus 608.

[0129] Among them:

[0130] The processor 602, the communications interface 604, and the memory 606 complete mutual communication through the communications bus 608.

[0131] The communications interface 604 is configured to communicate with other electronic devices or servers.

[0132] The processor 602 is configured to execute the program 610, and specifically can execute the related steps in the interface display adaptation method embodiments.

[0133] Specifically, the program 610 can include program code, and the program code includes computer operation instructions.

[0134] The processor 602 can be a central processing unit CPU, or an application specific integrated circuit ASIC, or one or more integrated circuits configured to implement one or more embodiments of the present application. One or more processors included in the smart device can be the same type of processor, such as one or more CPUs; or can be different types of processors, such as one or more CPUs and one or more ASICs.

[0135] The memory 606 is configured to store the program 610. The memory 606 can include a high-speed RAM memory, and can also include a non-volatile memory, such as at least one disk memory.

[0136] The program 610 can be specifically used to make the processor 602 perform the following operations:

[0137] In an alternative implementation, the program 610 is further configured to cause the processor 602 to perform the steps of the flare occurrence location determination method. The specific implementation of the steps in the program 610 can refer to the corresponding description in the steps and units in the above-described flare occurrence location determination method embodiments, and will not be described here. It can be clearly understood by those skilled in the art that, for the convenience and brevity of description, the specific working processes of the above-described devices and modules can refer to the corresponding process description in the foregoing method embodiments, and will not be described here.

[0138] The embodiments of the present application further provide a computer storage medium, which stores a computer program. The program is executed by a processor to implement the flare occurrence location determination method described in the embodiments of the present application.

[0139] It should be noted that, according to the implementation needs, each component / step described in the embodiments of the present application can be split into more components / steps, or two or more components / steps or part of the operations of the components / steps can be combined into a new component / step, to achieve the purpose of the embodiments of the present application.

[0140] The above-described method according to the embodiments of the present application can be implemented in hardware, firmware, or as software or computer code that can be stored in a recording medium such as a CD ROM, a RAM, a floppy disk, a hard disk or a magneto-optical disk, or be downloaded through a network from a remote recording medium or a non-transitory machine-readable medium originally stored in a local recording medium and then stored in a local recording medium, so that the method described herein can be processed by such software using a general-purpose computer, a special-purpose processor or programmable or special-purpose hardware (such as ASIC or FPGA). It can be understood that the computer, processor, microprocessor controller or programmable hardware includes a storage component (for example, RAM, ROM, flash memory, etc.) that can store or receive software or computer code, when the software or computer code is accessed and executed by the computer, processor or hardware, the interface display adaptation method described herein is implemented. In addition, when a general-purpose computer accesses the code for implementing the interface display adaptation method shown herein, the execution of the code will convert the general-purpose computer into a special-purpose computer for executing the interface display adaptation method shown herein.

[0141] Those skilled in the art can appreciate that the units and method steps of the examples described in connection with the embodiments disclosed herein can be implemented in electronic hardware, or in combination of computer software and electronic hardware. Whether the functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. A skilled person can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of the embodiments of the present application.

[0142] Finally, it should be noted that the above examples are only used to illustrate the technical solutions of the present application and are not limiting. Although the present application has been described in detail with reference to the examples, those skilled in the art should understand that modifications or equivalent replacements to the technical solutions of the present application can be made without departing from the spirit and scope of the present application, and all should be covered in the scope of the claims of the present application.

Claims

1. A method for determining the location of a solar flare, characterized in that, The method includes: When a flare is determined to occur, the peak time of the flare, the full solar image corresponding to the peak time, and the initial solar position of each solar active region are obtained. Based on the flare peak time and the initial solar surface position of each of the solar active regions, determine the first solar surface position of each of the solar active regions corresponding to the flare peak time; A coordinate system transformation is performed on the first solar surface position of each solar active region to determine the pixel position of each solar active region in the full solar surface image; Based on the pixel position of each solar active region and the full solar image, calculate the first pixel value of the region in the full solar image where each solar active region is located; The maximum value among multiple first pixel values ​​is determined, and the first solar surface position of the solar active region corresponding to the maximum value of the first pixel value is determined as the location of the flare.

2. The method according to claim 1, characterized in that, The step of determining the first solar surface position of each solar active region corresponding to the flare peak time based on the flare peak time and the initial solar surface position of each solar active region includes: Obtain the solar active region overview document and the publication time of the solar active region overview document; The initial solar surface position of each of the solar active regions is determined in the solar active region overview document; Based on the flare peak time, the release time, the initial solar surface position of each of the solar active regions, and the preset solar surface difference rotation formula, the first solar surface position of each of the solar active regions corresponding to the flare peak time is determined.

3. The method according to claim 1, characterized in that, The step of calculating the first pixel value of the region in the total solar image where each of the solar active regions is located, based on the pixel position of each of the solar active regions and the total solar image, includes: For each solar active region, the values ​​of each pixel in the region where the solar active region is located in the full solar image are obtained; Based on all the pixels in the region where the solar active region is located in the full-solar image, the sum of the pixel values ​​of all the pixels is calculated to obtain the first pixel value of the region where the solar active region is located in the full-solar image.

4. The method according to claim 3, characterized in that, After obtaining the pixel values ​​of each pixel in the region where the solar active region is located in the full-plane image, the method further includes: Based on the pixel value of each pixel, multiple target pixels are determined, wherein the pixel value of each target pixel is greater than a preset pixel threshold. The step of calculating the sum of pixel values ​​of all pixels in the region where the solar active region is located in the full-plane image, to obtain the first pixel value of the region where the solar active region is located in the full-plane image, includes: Based on the target pixels in the region where the solar active region is located in the full-plane image, the sum of the pixel values ​​of the target pixels is calculated to obtain the first pixel value of the region where the solar active region is located in the full-plane image.

5. The method according to claim 1, characterized in that, The method further includes: Obtain the flow rate values ​​of four consecutive 1-minute X-rays; A flare is determined to have occurred when the flow rate values ​​of the four consecutive 1-minute X-rays increase, and the flow rate value of the fourth 1-minute X-ray is greater than 1.4 times the flow rate value of the first 1-minute X-ray.

6. The method according to claim 1, characterized in that, The method further includes: Obtain the position of the first solar surface in the solar active region where the flare occurred; Based on the solar active region overview document, determine the active region number of the solar active region where the flare occurred; A flare warning message is generated based on the first solar surface location of the solar active region where the flare occurred and the active region number.

7. The method according to claim 6, characterized in that, The method further includes: The flare occurrence time, flare peak time, flare end time, and flare intensity are determined based on four consecutive 1-minute X-ray flow values. The solar flare's location on the first solar surface and its number within the active region are used to generate a flare warning message, including: The flare warning information is generated based on the flare occurrence time, the flare peak time, the flare end time, the flare intensity, the first solar surface location of the solar active region where the flare occurred, and the active region number.

8. A device for determining the location of a solar flare, characterized in that, The device includes: The acquisition module is used to acquire the flare peak time, the full solar image corresponding to the flare peak time, and the initial solar position of each solar active region when it is determined that a flare has occurred. The first calculation module is used to determine the first solar surface position of each solar active region corresponding to the peak time of the flare, based on the peak time of the flare and the initial solar surface position of each solar active region. The conversion module is used to perform coordinate system transformation on the first solar surface position of each of the solar active regions to determine the pixel position of each of the solar active regions in the full solar surface image; The second calculation module is used to calculate the first pixel value of the region in the total solar image of each solar active region based on the pixel position of each solar active region and the total solar image. The determination module is used to determine the maximum value among multiple first pixel values, and to determine the first solar surface position of the solar active region corresponding to the maximum value of the first pixel value as the location of the flare.

9. An electronic device, comprising: The processor, memory, communication interface, and communication bus are provided, wherein the processor, memory, and communication interface communicate with each other via the communication bus. The memory is used to store at least one executable instruction that causes the processor to perform the operation corresponding to the flare location determination method as described in any one of claims 1-7.

10. A computer storage medium having a computer program stored thereon, which, when executed by a processor, implements the method for determining the location of a flare as described in any one of claims 1-7.

Images