A method for acquiring horizontal propagation phase velocity of gravitational wave based on networking of three same series micro-light satellites
By acquiring low-light nighttime data from three satellites, performing preprocessing and parallax correction model calibration, and calculating the spacing of gravity wave fringes, the problem of the inability to calculate the horizontal propagation phase velocity of gravity waves in existing technologies has been solved, achieving higher precision gravity wave velocity measurement.
Patent Information
- Application Number
- CN202411436416.8
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2024-10-15
- Publication Date
- 2025-12-19
- Estimated Expiration
- 2044-10-15
AI Technical Summary
There is no existing technology that uses networked data from three satellites of the same series to calculate the horizontal propagation phase velocity of gravity waves in the mesosphere.
By acquiring nighttime low-light data from three low-light satellites of the same series, preprocessing was performed to extract gravity wave fringes. The position of the gravity wave fringes was corrected using a parallax correction model, and the spacing between the gravity wave fringes was calculated to obtain the horizontal propagation phase velocity.
It improves the temporal resolution of low-light data, enabling accurate calculation of the horizontal propagation phase velocity of gravity waves, thus filling a gap in existing technology.
Smart Images
Figure CN119291797B_ABST
Abstract
Description
TECHNICAL FIELD
[0001] The application belongs to the field of remote sensing technology analysis, and particularly relates to a method for acquiring gravity wave horizontal propagation phase velocity based on three satellites in the same series of micro-light satellites. BACKGROUND
[0002] Visible Infrared Imaging Radiometer Suite (VIIRS) can detect the light radiance information of the earth surface at night, and then form light data transmitted back to the ground by satellite. In the past two decades, night micro-light detection sensors have been launched one after another, and the research and application technology related to such light data has also developed rapidly. In 2011, the S-NPP (Suomi-National Polar-orbiting Partnership) satellite launched by the United States carried the VIIRS (Visible Infrared Imaging Radiometer Suite) / DNB (Day Night Band) micro-light load with strong night detection capability. The spaceborne micro-light imager is originally used for imaging the atmosphere and the ground environment under low illumination conditions at night. Miller et al. (2012) first found that the existence of the middle-top gravity wave usually makes the airglow radiation near it present obvious strip image characteristics, providing a new method for the observation of the middle-top gravity wave. At present, spaceborne micro-light imagers are entering an era of rapid development, and satellite micro-light data is gradually becoming abundant. After S-NPP, the United States launched JPSS-1 and JPSS-2 satellites carrying the same type of micro-light load in November 2017 and 2022 respectively. The three satellites run on the same orbit, so that the time resolution of micro-light observation data is significantly improved. In about 50 minutes, three gravity wave observations will be performed in succession, and the gravity wave propagation phase velocity parameter can be more accurately extracted. The SDR data set developed based on the VIIRS observation data of the S-NPP, JPSS-1 and JPSS-2 satellites can provide daily satellite light data at night. This product is publicly released by NOAA, and the website is:
[0003] https: / / www.aev.class.noaa.gov / saa / products / search?sub_id=0&datatype_famil y=VIIRS_SDR&submit.x=20&submit.y=5, this paper mainly uses the SVDNB and GDNBO products in the SDR data set composed of 3072x4064 pixels. Among them, the SVDNB product provides the entrance pupil radiance which has been subjected to business radiation calibration; the GDNBO product provides the observation geometry information such as longitude and latitude and satellite zenith angle corresponding to the entrance pupil radiance.
[0004] In the prior art, no relevant research is found on network research on faint light data observed by three satellites and calculation of horizontal propagation phase velocity of middle layer top gravity wave through the faint light data. Therefore, the application provides a method for obtaining horizontal propagation phase velocity of gravity wave based on network of three faint light satellites in the same series. SUMMARY
[0005] In view of the above problems, the application aims to fill the blank of the prior art in calculating the horizontal propagation phase velocity of gravity wave by using satellite faint light data, and provides a method for obtaining horizontal propagation phase velocity of gravity wave based on network of three faint light satellites in the same series, which can obtain the horizontal propagation phase velocity of middle layer top gravity wave based on night faint light data.
[0006] In order to achieve the above technical effects, in order to achieve the above purpose, the application provides a method for obtaining horizontal propagation phase velocity of gravity wave based on network of three faint light satellites in the same series, comprising:
[0007] S1: obtaining night faint light data of three satellites, the night faint light data comprising SVDNB and GDNBO;
[0008] S2: pre-processing the night faint light data, and extracting gravity wave stripes in the pre-processed night faint light image;
[0009] S3: correcting the extracted gravity wave stripes by a parallax correction model, drawing the corrected gravity wave stripes together for pairing, and obtaining a corresponding relationship of the gravity wave stripes;
[0010] S4: based on the corresponding relationship of the gravity wave stripes in S3, extracting interval distances between the gravity wave stripes, and then calculating the horizontal propagation phase velocity of the gravity wave.
[0011] Preferably, the night faint light data pre-processing in S2 comprises:
[0012] S2.1: based on SVDNB, obtaining an entrance pupil radiance;
[0013] S2.2: based on GDNBO, obtaining observation geometry information of longitude and latitude and satellite zenith angle corresponding to the entrance pupil radiance;
[0014] S2.3: based on geographical information contained in GDNBO, constructing a GLT file in ENVI software, and performing geometric correction on the entrance pupil radiance in S2.1 to extract the gravity wave stripes after geometric correction.
[0015] Preferably, the correction by the parallax correction model in S3 comprises:
[0016] S3.1: the height of satellite from the ground H; the radius of the earth R; the height of the middle layer top from the ground h, the pixel resolution d of the airglow image, O is the earth center; S is the satellite; A is the subsatellite point position; W is the real position of the gravity wave; T is the corresponding position of the gravity wave on the ground; E is the observed position of the gravity wave in the airglow image, the corresponding row number is r and the column number is c in the observed data E E ; the length of the circular arc ET is the distance deviation of the gravity wave between the observed position and the real corresponding position on the ground;
[0017] When the satellite zenith angle θ is known, then:
[0018] ∠OES = π - θ
[0019] ∠OES is the angle formed by the earth center and the satellite with the observed position of the gravity wave in the airglow image as the vertex;
[0020] S3.2: in ΔOES formed by the earth center, the observed position of the gravity wave in the airglow image and the satellite, and in ΔOWS formed by the earth center, the real position of the gravity wave and the satellite, according to the sine theorem, we have:
[0021]
[0022]
[0023] Therefore:
[0024] (R + H) sin ∠OSE = R sin (π - θ)
[0025]
[0026] Then ET is:
[0027] ET = R · ∠WOE
[0028] Wherein, ∠OSE is the angle formed by the satellite with the earth center and the observed position of the gravity wave in the airglow image as the vertex; ∠OWS is the angle formed by the real position of the gravity wave with the earth center and the satellite as the vertex; ∠WOE is the angle formed by the real position of the gravity wave and the observed position of the gravity wave in the airglow image with the earth center as the vertex;
[0029] S3.3: on the basis of S3.2, the pixel number N of the deviation is,
[0030]
[0031] Wherein, d represents the pixel resolution of the airglow image;
[0032] S3.4: the row number r of the gravity wave at the corresponding position T on the ground surface after the processing of S3.1-S3.3 T and the column number c T are respectively:
[0033] r T = r E
[0034]
[0035] S3.5: the row number r of the gravity wave at the corresponding position T on the ground surface obtained according to S3.4 T and the column number c T , that is, the longitude and latitude of the gravity wave at the corresponding position T on the ground surface can be found according to the GDNBO.
[0036] Preferably, step 4 is specifically:
[0037] Based on the distance D between the longitude and latitude of the gravity wave at the corresponding position T on the ground surface obtained according to S3.5 and the satellite scanning interval time T, the horizontal propagation phase velocity V is defined, and then:
[0038]
[0039] Compared with the prior art, the present application has the beneficial effects that:
[0040] The present application uses a different satellite data source from the past, and the satellites for obtaining micro-light data are increased from NPP to NPP, JPSS-1 and JPSS-2, the time resolution of the night micro-light data of the polar orbit satellite is improved, a method for observing the gravity wave based on three polar orbit satellites in the same orbit is proposed, a parallax correction method for the position of the gravity wave is proposed, and after the position of the gravity wave stripe after the geometric correction is extracted and the parallax correction is performed, the horizontal propagation phase velocity of the gravity wave is calculated. BRIEF DESCRIPTION OF DRAWINGS
[0041] The accompanying drawings are included to provide a further understanding of the present application, and constitute a part of the specification, and are used together with embodiments of the present application to explain the present application, and do not constitute a limitation on the present application.
[0042] In the drawings:
[0043] Figure 1 The flowchart of the embodiment of the present application;
[0044] Figure 2 The parallax correction model of the embodiment of the present application;
[0045] Figure 3 The position of the gravity wave after the parallax correction of the embodiment of the present application;
[0046] Figure 4A schematic diagram of the calculation of the horizontal propagation phase velocity of gravity waves according to an embodiment of the present application. DETAILED DESCRIPTION
[0047] The preferred embodiments of the present application are described below in conjunction with the accompanying Figures 1-4 The preferred embodiments described herein are intended to be illustrative only and are not intended to limit the scope of the present application.
[0048] Embodiments
[0049] A method for obtaining the horizontal propagation phase velocity of gravity waves based on three satellites in the same series of micro-satellite networking, comprising:
[0050] S1: Obtain night micro-light data of three satellites, wherein the night micro-light data comprises SVDNB and GDNBO.
[0051] S2: Pre-process the night micro-light data, and extract the gravity wave stripes in the pre-processed night micro-light image (under the condition that the moonlight is weak and there is no cloud, the gravity wave in the micro-light image will present the characteristics of alternating bright and dark stripes, and the clearly identifiable stripes are drawn out). The pre-processing of the night micro-light data in S2 comprises:
[0052] S2.1: Obtain the entrance pupil radiance based on the SVDNB.
[0053] S2.2: Obtain the observation geometry information of the longitude, latitude and satellite zenith angle corresponding to the entrance pupil radiance based on the GDNBO.
[0054] S2.3: Based on the geographical information contained in the GDNBO, construct a GLT file in the ENVI software, and perform geometric correction on the entrance pupil radiance of S2.1 to extract the gravity wave stripes after geometric correction.
[0055] S3: Correct the extracted gravity wave stripes through a parallax correction model, draw the corrected gravity wave stripes together for pairing, and obtain the corresponding relationship of the gravity wave stripes.
[0056] A schematic diagram of the parallax correction model of an embodiment of the present application is shown in Figure 2 , a schematic diagram of the position of the gravity wave after parallax correction is shown in Figure 3 , and a schematic diagram of the calculation of the horizontal propagation phase velocity of gravity waves is shown in Figure 4 It can be seen that the position of the gravity wave stripe after parallax correction has a large change.
[0057] The correction by the parallax correction model in S3 comprises:
[0058] S3.1: the height H of satellite from the ground, the radius R of the earth, the height h of the middle layer top from the ground, the pixel resolution d of the airglow image, O is the earth center, S is the satellite, A is the subsatellite point position, W is the real position of the gravity wave, T is the corresponding position of the gravity wave on the ground, E is the observed position of the gravity wave in the airglow image, the corresponding row number is r and the column number is c in the observation data E E The length of the circular arc ET is the distance deviation of the gravity wave between the observed position and the real corresponding position on the ground.
[0059] When the satellite zenith angle θ is known, then:
[0060] ∠OES = π - θ
[0061] ∠OES is the angle formed by the earth center and the satellite with the observed position of the gravity wave in the airglow image as the vertex.
[0062] S3.2: In ΔOES formed by the earth center, the observed position of the gravity wave in the airglow image and the satellite, and in ΔOWS formed by the earth center, the real position of the gravity wave and the satellite, according to the sine theorem, we have:
[0063]
[0064] Therefore:
[0065] (R + H) sin ∠OSE = R sin (π - θ)
[0066]
[0067] Then ET is:
[0068] ET = R ∠WOE
[0069] Where ∠OSE is the angle formed by the satellite and the earth center with the observed position of the gravity wave in the airglow image as the vertex, ∠OWS is the angle formed by the real position of the gravity wave and the earth center with the satellite as the vertex, and ∠WOE is the angle formed by the earth center with the real position of the gravity wave and the observed position of the gravity wave in the airglow image as the vertex.
[0070] S3.3: Based on S3.2, the number of pixel deviation N is:
[0071]
[0072] Where d represents the pixel resolution of the airglow image.
[0073] S3.4: The row number r and column number c of the gravity wave on the ground corresponding position T after S3.1-S3.3 processing T T respectively,
[0074] r T = r E
[0075]
[0076] S3.5: the row number r of the gravity wave at the corresponding position T on the ground surface obtained according to S3.4 T and the column number c T , that is, the longitude and latitude of the gravity wave at the corresponding position T on the ground surface can be found according to the GDNBO.
[0077] S4: based on the corresponding relationship of the gravity wave stripes in S3, the interval distance between the gravity wave stripes is extracted, and then the horizontal propagation phase velocity of the gravity wave is calculated. Specifically,
[0078] based on the distance D between the longitude and latitude of the gravity wave at the corresponding position T on the ground surface in S3.5 and the satellite scanning interval time T, the horizontal propagation phase velocity V is defined as:
[0079]
[0080] As can be seen from the summary and the examples, the present application obtains night micro-light data through a visible light / infrared radiation imager, extracts the positions of the gravity wave stripes through geometric correction, obtains the real positions of the gravity wave through a parallax correction model, thereby studying the propagation process of the gravity wave in the micro-light images of three satellites, and then extracts the corresponding distance in the stripes and calculates the horizontal propagation phase velocity of the gravity wave, thereby making up for the defect that the horizontal propagation phase velocity of the gravity wave cannot be calculated through night satellite light data at present.
[0081] The basic principles, main features and advantages of the present application are shown and described above. It should be understood by those skilled in the art that the present application is not limited by the above examples, and the above examples and descriptions in the specification are only to illustrate the principles of the present application. Without departing from the spirit and scope of the present application, various changes and improvements can be made to the present application, and these changes and improvements all fall within the scope of the present application. The scope of protection of the present application is defined by the appended claims and their equivalents.
Claims
1. A method for obtaining horizontal propagation phase velocity of gravitational waves based on networking of three microsatellites in the same series, characterized in that: Comprise: S1: Obtain three satellite nightglow data, nightglow data includes SVDNB and GDNBO; S2: Preprocess the nightglow data, extract the gravity wave stripes in the preprocessed nightglow image; S3: The extracted gravity wave stripes are corrected by the parallax correction model, and the corrected gravity wave stripes are drawn together for pairing to obtain the corresponding relationship of the gravity wave stripes; S3 parallax correction model for correction includes: S3.1: Height of satellite from the ground H ; Earth radius R ; Height of the top of the middle layer from the ground where the gravity wave is located h ; Pixel resolution of the glimmer image d ; O is the Earth's center; S is the satellite; A is the sub-satellite point position; W is the true position of the gravity wave; T is the corresponding position on the ground of the gravity wave; E is the observed position of the gravity wave in the glimmer image, the corresponding row number in the observation data is r E , and the column number is c E ; Arc length is the distance deviation of the gravity wave at the observed position and the true corresponding position on the ground; Satellite zenith angle θ It is known that: ； ∠OES is the angle formed by the observation position of the gravity wave in the nightglow image and the satellite and the center of the earth; S3.2: The three points of the Earth's center, the observed position of the gravitational wave in the micro-light image and the satellite form a OES, and the three points of the Earth's center, the real position of the gravitational wave and the satellite form a OWS, by the sine theorem, we have: ； ； Therefore: ； ； ； Then is: ； Wherein, ∠OSE is the angle formed by the satellite and the center of the earth and the observation position of the gravity wave in the nightglow image; ∠OWS is the angle formed by the real position of the gravity wave and the satellite and the center of the earth; ∠WOE is the angle formed by the real position of the gravity wave and the observation position of the gravity wave in the nightglow image with the center of the earth as the vertex; S3.3: On the basis of S3.2, the number of pixels of deviation N is, ； wherein, d represents the pixel resolution of the low-light image; S3.4: the row number of the gravity wave on the ground surface corresponding to the position T after the processing of S3.1-S3.3 r T and the column number c T are respectively: ； ； S3.5: the row number of the gravity wave at the corresponding position T on the ground surface according to S3.4 r T and the column number c T that is, the longitude and latitude of the gravity wave at the corresponding position T on the ground surface can be found according to the GDNBO; S4: Based on the corresponding relationship of the gravity wave stripes in S3, the interval distance between the gravity wave stripes is extracted, and then the horizontal propagation phase velocity of the gravity wave is calculated.
2. The method for obtaining the horizontal propagation phase velocity of the gravity wave based on three microglow satellites in the same series according to claim 1, wherein the nightglow data preprocessing comprises: S2.1: Based on SVDNB, obtain the entrance pupil radiance; S2.2: Based on GDNBO, obtain the observation geometry information of the entrance pupil radiance corresponding to the longitude and latitude and satellite zenith angle; S2.3: Based on the geographical information contained in GDNBO, construct GLT file in ENVI software, and perform geometric correction on the entrance pupil radiance of S2.1 to extract the geometrically corrected gravity wave stripes.
3. The method of claim 2, wherein the three satellites are of the same series. Step 4 is specifically: the distance between the longitude and latitude of the corresponding position T on the ground surface based on S3.5 gravitational waves D , and the satellite scanning interval time T , the horizontal propagation phase velocity is defined as V then: 。