CN111366930A - Cloud top height verification method and device, computer equipment and storage medium - Google Patents
Cloud top height verification method and device, computer equipment and storage medium Download PDFInfo
- Publication number
- CN111366930A CN111366930A CN202010156675.0A CN202010156675A CN111366930A CN 111366930 A CN111366930 A CN 111366930A CN 202010156675 A CN202010156675 A CN 202010156675A CN 111366930 A CN111366930 A CN 111366930A
- Authority
- CN
- China
- Prior art keywords
- top height
- data
- cloud top
- cloud
- echo
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 54
- 238000012795 verification Methods 0.000 title claims abstract description 28
- 238000004422 calculation algorithm Methods 0.000 claims description 34
- 238000006243 chemical reaction Methods 0.000 claims description 34
- 238000002310 reflectometry Methods 0.000 claims description 32
- 238000004590 computer program Methods 0.000 claims description 25
- 238000010586 diagram Methods 0.000 claims description 21
- 238000005259 measurement Methods 0.000 claims description 9
- 238000005516 engineering process Methods 0.000 abstract description 3
- 238000004364 calculation method Methods 0.000 description 11
- 238000012545 processing Methods 0.000 description 6
- 238000001514 detection method Methods 0.000 description 4
- 238000002592 echocardiography Methods 0.000 description 4
- 238000001556 precipitation Methods 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- PZWWYAHWHHNCHO-FGHAYEPSSA-N (4r,7s,10s,13r,16s,19r)-n-[(2s,3r)-1-amino-3-hydroxy-1-oxobutan-2-yl]-19-[[(2r)-2-amino-3-phenylpropanoyl]amino]-10-(3-aminopropyl)-7-[(1r)-1-hydroxyethyl]-16-[(4-hydroxyphenyl)methyl]-13-(1h-indol-3-ylmethyl)-3,3-dimethyl-6,9,12,15,18-pentaoxo-1,2-dithia Chemical compound C([C@H]1C(=O)N[C@H](CC=2C3=CC=CC=C3NC=2)C(=O)N[C@@H](CCCN)C(=O)N[C@H](C(=O)N[C@@H](C(SSC[C@@H](C(=O)N1)NC(=O)[C@H](N)CC=1C=CC=CC=1)(C)C)C(=O)N[C@@H]([C@H](O)C)C(N)=O)[C@@H](C)O)C1=CC=C(O)C=C1 PZWWYAHWHHNCHO-FGHAYEPSSA-N 0.000 description 2
- 230000004075 alteration Effects 0.000 description 2
- 239000003086 colorant Substances 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000002955 isolation Methods 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 108700040302 phenylalanyl-cyclo(cysteinyltyrosyl-tryptophyl-ornithyl-threonyl-penicillamine)threoninamide Proteins 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 241000208818 Helianthus Species 0.000 description 1
- 235000003222 Helianthus annuus Nutrition 0.000 description 1
- 238000009833 condensation Methods 0.000 description 1
- 230000005494 condensation Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000010339 dilation Effects 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000003709 image segmentation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/95—Radar or analogous systems specially adapted for specific applications for meteorological use
- G01S13/958—Theoretical aspects
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/46—Measurement of colour; Colour measuring devices, e.g. colorimeters
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/02—Systems using the reflection of electromagnetic waves other than radio waves
- G01S17/06—Systems determining position data of a target
- G01S17/08—Systems determining position data of a target for measuring distance only
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
- G01S7/4052—Means for monitoring or calibrating by simulation of echoes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/90—Details of database functions independent of the retrieved data types
- G06F16/901—Indexing; Data structures therefor; Storage structures
- G06F16/9017—Indexing; Data structures therefor; Storage structures using directory or table look-up
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- Theoretical Computer Science (AREA)
- Databases & Information Systems (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The present application relates to the field of computer technologies, and in particular, to a cloud top height verification method and apparatus, a computer device, and a storage medium. The method comprises the following steps: acquiring meteorological radar data and meteorological satellite data; calculating the echo peak height according to the meteorological radar data; searching the cloud top height corresponding to the meteorological satellite data from the pre-generated incidence relation; and verifying the cloud top height by using the echo top height. By adopting the method, the accuracy of the cloud top height obtained by utilizing meteorological satellite data inversion can be accurately verified.
Description
Technical Field
The present application relates to the field of computer technologies, and in particular, to a cloud top height verification method and apparatus, a computer device, and a storage medium.
Background
The cloud top height is one of cloud physical parameters, represents the maximum height reached by condensation of water vapor in the atmosphere, can influence the balance of ground gas radiation, and has a remarkable adjusting effect on the balance of atmospheric energy balance.
In the traditional technology, satellite data and auxiliary data such as an atmospheric transmittance profile and the like need to be acquired for calculating the height of the cloud top, and more variables are introduced by the method for calculating the height of the cloud top by using the auxiliary data, so that the calculation accuracy of the height of the cloud top is low, and the efficiency of checking the height of the cloud top is reduced.
Disclosure of Invention
In view of the foregoing, it is desirable to provide a method, an apparatus, a computer device and a storage medium capable of improving cloud top height verification efficiency.
A cloud top height verification method comprises the following steps:
acquiring meteorological radar data and meteorological satellite data;
calculating the echo peak height according to the meteorological radar data;
searching the cloud top height corresponding to the meteorological satellite data from the pre-generated incidence relation;
and verifying the cloud top height by using the echo top height.
In one embodiment, the cloud top height is verified by using the echo top height, and the method comprises the following steps:
converting the data types of the cloud top height and the echo top height corresponding to the same area range into consistent data types by using a conversion algorithm;
and calculating the absolute value deviation between the echo top height and the cloud top height with the consistent data types, and verifying the cloud top height according to the absolute value deviation.
In one embodiment, the converting algorithm is used to convert the data types of the cloud top height and the echo top height corresponding to the same region range into a consistent data type, and the converting algorithm includes:
acquiring cloud top height and echo top height corresponding to the same area range;
and converting the longitude and latitude numerical value corresponding to the cloud top height into a distance numerical value according to a pre-generated conversion algorithm, wherein the conversion algorithm is specifically obtained according to the longitude and latitude numerical value corresponding to the position of the radar station, the earth radius and the equator radius.
In one embodiment, the calculating the echo peak height according to the weather radar data comprises:
acquiring the reflectivity of each data point in the meteorological radar data and echoes corresponding to the reflectivities;
extracting the highest elevation angle corresponding to the data point of the echo greater than the preset threshold value;
and calculating the echo top height of the data point corresponding to the highest elevation according to a height measurement formula.
In one embodiment, the searching for the cloud top height corresponding to the meteorological satellite data from the pre-generated association relationship comprises:
extracting brightness temperature values corresponding to two preset channels from meteorological satellite data;
calculating the difference value of the brightness temperature values in the two channels to obtain the brightness temperature difference;
and acquiring the cloud top height corresponding to the brightness temperature value and/or the brightness temperature difference according to a pre-generated incidence relation, wherein the incidence relation is specifically obtained according to the relation between the brightness temperature value and/or the brightness temperature difference and the satellite cloud top height in the cloud top height distribution diagram.
In one embodiment, the method for generating the association relationship includes:
acquiring a cloud top height distribution map corresponding to a preset channel in meteorological satellite data;
reading a color value corresponding to a data pair consisting of each bright temperature value and each bright temperature difference from the cloud top height distribution map;
searching the cloud top height corresponding to the color value;
and establishing an association relation among the bright temperature value, the bright temperature difference and the cloud top height.
In one embodiment, before reading the color value corresponding to the data pair consisting of each bright temperature value and each bright temperature difference from the cloud top height distribution map, the method further includes:
reading the corresponding pixel values of the data pairs consisting of the brightness temperature values and the brightness temperature differences from the cloud top height distribution map;
and when the pixel value corresponds to the achromatic value, acquiring a neighborhood pixel value corresponding to an adjacent region of the data pair, and assigning the pixel value of the data pair according to the neighborhood pixel value in the cloud top height distribution map to generate a cloud top height distribution map after reassignment.
A cloud top height verification device, the device comprising:
the acquisition module is used for acquiring meteorological radar data and meteorological satellite data;
the echo top height calculating module is used for calculating the echo top height according to the meteorological radar data;
the cloud top height searching module is used for searching the cloud top height corresponding to the meteorological satellite data from the pre-generated incidence relation;
and the verification module is used for verifying the cloud top height by utilizing the echo top height.
A computer device comprising a memory storing a computer program and a processor implementing the steps of the method when the processor executes the computer program.
A computer-readable storage medium, on which a computer program is stored which, when being executed by a processor, carries out the steps of the above-mentioned method.
The cloud top height verification method, the cloud top height verification device, the computer equipment and the storage medium acquire different types of meteorological data, such as meteorological radar data and meteorological satellite data; then calculating the echo peak height according to the meteorological radar data; searching the cloud top height corresponding to the meteorological satellite data from the pre-generated incidence relation; and then obtaining height values calculated by different types of meteorological data at the same position point, and then verifying the cloud top height by using the echo top height so as to verify the accuracy of the cloud top height obtained by inversion according to meteorological satellite data.
Drawings
FIG. 1 is a diagram illustrating an exemplary implementation of a cloud ceiling height verification method;
FIG. 2 is a flow diagram illustrating a cloud top height verification method according to an embodiment;
FIG. 3 is a flow chart of a method for verifying accuracy of satellite data inversion cloud top height provided in one embodiment;
fig. 4 is a schematic flow chart illustrating a process of searching for a cloud top height corresponding to weather satellite data in a pre-generated association relationship provided in an embodiment;
FIG. 5 is a schematic diagram of a cloud top height distribution map provided in one embodiment;
FIG. 6 is a block diagram of an exemplary cloud ceiling height verification apparatus;
FIG. 7 is a diagram illustrating an internal structure of a computer device according to an embodiment.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application more apparent, the present application is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the present application and are not intended to limit the present application.
The cloud top height verification method provided by the application can be applied to the application environment shown in fig. 1. The radar 102 and the satellite 103 communicate with a server 104 via a network. The server 104 acquires weather radar data from the radar 102, and the server 104 acquires weather satellite data from the satellite 103; the server 104 calculates the echo peak height according to the meteorological radar data; the server 104 searches for the cloud top height corresponding to the meteorological satellite data from the pre-generated association relation; the server 104 verifies the cloud top height using the echo top height. The server 104 may be implemented as a stand-alone server or a server cluster composed of a plurality of servers.
In one embodiment, as shown in fig. 2, a cloud top height verification method is provided, which is described by taking the method as an example applied to the server 104 in fig. 1, and in other embodiments, the method may also be applied to a terminal, and the method includes the following steps:
step 210, acquiring weather radar data and weather satellite data.
The meteorological satellite data is data obtained by meteorological observation of the earth and the atmosphere thereof by a satellite, and the meteorological satellite data can contain a plurality of channel data, wherein different channel data correspond to different central wavelengths, for example, 16 channels can be provided in total, wherein the central wavelength corresponding to the 14 th channel is 11.2 μm, and the central wavelength corresponding to the 15 th channel is 12.3 μm. In addition, the weather satellite data may also include a plurality of visible light bands, and data of different visible light bands may be synthesized into a color map, for example, the color map may include 3 visible light bands, a full-color map may be synthesized using 3 visible light bands, and further, the color map may also include 3 near-infrared bands, 10 infrared bands, and the like. The frequency of observation of the satellites may be 1 time every 10 minutes.
The radar can detect the reflectivity of precipitation particles and the speed of the radial direction of the precipitation particles, the meteorological radar data are stored in a binary mode, the meteorological radar data comprise reflectivity data, radial speed data and spectral width data, the radar detects precipitation cloud bodies from a low elevation angle to a high elevation angle, 9 elevation angle layers are provided in total, each elevation angle layer has 360 radial directions, each radial direction has 460 reflectivity values for the reflectivity, the resolution ratio is 1km, each radial direction has 920 radial speed values for the radial speed, and the resolution ratio is 0.25 km. It should be noted that the radar may be a conventional pulse radar or a doppler radar, which is not limited herein.
And step 220, calculating the echo top height according to the meteorological radar data.
The echo peak height refers to a height at which an echo with intensity greater than a certain threshold value exists in the reflectivity factor, and specifically, the echo peak height refers to a height at which an echo with intensity of 30dBz exists in the upper portion of a thunderstorm cloud. Before the echo peak height is calculated by using the meteorological radar data, the server can also preprocess the acquired meteorological radar data, and the preprocessing comprises the following steps: the server carries out isolation point removing processing and/or median filtering processing on the reflectivity data in the radar base data, wherein the isolation point removing processing refers to removing isolated echo points in the radar. By performing de-outlier processing on the reflectivity, clutter may be removed to make the reflectivity data smoother.
The server extracts the preprocessed reflectivity data from the meteorological radar data, and then calculates the echo top height according to the reflectivity data, for example, the echo top height can be calculated by using a preset calculation algorithm.
And step 230, searching the cloud top height corresponding to the meteorological satellite data from the pre-generated association relation.
The cloud top height is one of cloud physical parameters and is the maximum height to which water vapor in the atmosphere condenses. The cloud top height can influence the balance of ground gas radiation, has a remarkable adjusting effect on the balance of atmospheric energy balance, and meanwhile, the parameters of the cloud top height and the like also have important embodiments in the fields of aerial weather guarantee, numerical weather forecast and the like, and the determination of the parameters of the cloud top height and the like has important practical significance on the aspects of atmospheric physics, climate research, weather guarantee and the like.
The storage mode of the incidence relation can be a lookup table, and the incidence relation among a bright temperature value, a bright temperature difference and a cloud top height in a certain channel in the meteorological satellite data is stored in the lookup table, wherein the storage format of the incidence relation can be an excel format lookup table. Specifically, the server reads the obtained lookup table, and obtains the cloud top height corresponding to the bright temperature value and the bright temperature difference from the lookup table.
And 240, verifying the cloud top height by using the echo top height.
Because of the difference of the detection principles of the meteorological radar and the meteorological satellite, the echo top height or cloud top height calculated by the base data of the meteorological radar and the meteorological satellite according to the corresponding algorithm can have difference, but experiments show that although the height values calculated by the meteorological radar data and the meteorological satellite data have difference, certain correlation still exists between the two height values, for example, the correlation relationship between the two height values can be fitted from the statistical angle.
Specifically, the similarity between the echo top height and the cloud top height can be calculated, and the difference between the echo top height and the cloud top height can be measured according to the size of the similarity value.
Specifically, 100 samples may be selected as verification samples, and an average value, that is, an average deviation, of absolute deviations between the cloud top heights and the echo top heights corresponding to all the samples is calculated, where a calculation formula is shown in formula (1):
wherein, AD is the average deviation, x is the cloud top height value calculated according to weather satellite data, x1 is the average value of the echo height calculated according to weather radar data in the sample, and N is the total number of the sample. When the average deviation AD is larger than a preset value, the difference between the cloud top height inverted by the meteorological satellite data and the echo height calculated by the meteorological radar data is larger, and the accuracy of the cloud top height inverted by the meteorological satellite data is judged to be lower, so that the cloud top height is invalid data. On the contrary, when the average deviation AD is within the range of the preset value, the difference between the cloud top height inverted by the meteorological satellite data and the echo height calculated by the meteorological radar data is within the preset range, and the cloud top height obtained by the meteorological satellite data inversion is judged to be high in accuracy and to be effective data.
In consideration of the difference between radar and satellite detection principles, the echo top height and the cloud top height calculated according to the data of the radar and the satellite are different, the radar can detect the echo only when small water drops are formed in the cloud, the satellite can detect the echo only when the cloud exists, and a formula fitted according to statistics in advance is only a maximum approximate relation between the two, so that the difference between the cloud top height calculated by the radar and the cloud top height calculated by the satellite is within 2km, and the result that the satellite data are used for calculating the cloud top height by the method is accurate and reasonable. The average deviation obtained by verifying the method by using actual data is 1km, so that the method for efficiently calculating the height of the cloud top is accurate.
In this embodiment, in order to verify the accuracy of the cloud top height of the meteorological satellite data inversion, the echo top height is calculated by using meteorological radar data corresponding to a radar, and the difference between the cloud top height and the true value is verified by using the echo top height as the true value, so as to obtain the accuracy of the cloud top height of the meteorological satellite data inversion.
In one embodiment, the cloud top height is verified by using the echo top height, and the method comprises the following steps: and converting the data types of the cloud top height and the echo top height corresponding to the same area range into consistent data types by using a conversion algorithm.
To verify the cloud top height by using the echo top height, firstly, the consistency of the two data types needs to be ensured. The conclusion that a certain correlation exists between the echo top height and the cloud top height is obtained through data analysis of a large amount of statistical data, so that a conversion algorithm between the echo top height and the cloud top height can be obtained according to the correlation relation, and the data types of the cloud top height and the echo top height can be converted into consistent data types through the conversion algorithm. Specifically, the relationship between the echo peak height and the cloud peak height calculated by the FY2 meteorological satellite data is fitted, and if the correlation is the highest, the conversion algorithm corresponding to the fitting relationship between the echo peak height calculated by the SA model radar and the cloud peak height calculated by the FY2 satellite is shown in formula (2) when the combined reflectivity is less than 30dBz and the threshold value is 1 dBz.
1dBz ET＝0.632*FY_2C CTOP+2.053 (2)
Wherein, 1dBz ET is the echo top height calculated when the radar threshold is 1dBz, and FY _2C CTOP is the cloud top height calculated by using FY2 meteorological satellites. The echo peak height calculated by the weather radar can be converted into the cloud peak height according to equation (2). Therefore, the radar data with the combined reflectivity of less than 30dBz and the satellite data of the corresponding time are selected to verify the accuracy of the result, and specifically, 100 cases can be selected as verification samples.
In the embodiment, by converting the projection forms of the echo top height and the cloud top height into the same data type, each point of the two is the same point. And then the difference between the echo top height and the cloud top height corresponding to the same position point can be calculated, and the cloud top height can be verified by utilizing the echo top height at the same position.
In one embodiment, the converting algorithm is used to convert the data types of the cloud top height and the echo top height corresponding to the same region range into a consistent data type, and the converting algorithm includes: and acquiring the cloud top height and the echo top height corresponding to the same area range.
Because the detection principle that the radar corresponds to the satellite and the detection resolution ratio are different to a certain extent, before the cloud top height is verified by using the echo top height, the positions of the meteorological radar data and the meteorological satellite data are required to be consistent, and the fact that the cloud top height is verified by using the echo top height can be guaranteed only if the positions are consistent. Specifically, the server may verify the cloud top height of the location point according to the absolute value deviation by calculating the absolute value deviation between the echo top height and the cloud top height corresponding to the same location point.
The display range of Doppler weather radar of SA waveband of WSR-98D used is generally 920km × 920km, while the sunflower No. 8 satellite data can cover half range of the earth, and in order to verify the cloud top height calculated by the weather satellite data by using the echo top height calculated by the weather radar data, an area consistent with the coverage range of the echo top height in the weather satellite data needs to be intercepted.
Specifically, the longitude and latitude numerical value corresponding to the cloud top height is converted into a distance numerical value according to a pre-generated conversion algorithm, and the conversion algorithm is specifically obtained according to the longitude and latitude numerical value corresponding to the position of the radar station, the earth radius and the equator radius.
Since weather satellite data is stored in a projection format with equal longitude and latitude, the resolution is 0.01 longitude × 0.01.01 latitude, and the resolution of the cloud ceiling height calculated by using weather radar data is 1km × 1km, projection conversion is required, and the cloud ceiling height calculated by using the satellite is converted from 0.01 longitude × 0.01.01 latitude to the format of 1km × 1 km.
x1＝PI*eqRAD*math.cos(i*PI/180)/180 (3)
x＝x1*(i-Radarlong)+0.5 (4)
y＝(j-Radarlat)*PI*eqRAD/180+0.5 (5)
Wherein, (x, y) is the pixel point position of the cloud top height corresponding to the satellite with the same radar resolution, eqRAD is the equator radius, eqRAD is 6378.137km, and PI is 3.1415926.
In this embodiment, a conversion algorithm obtained by analyzing the data in advance converts the longitude and latitude values corresponding to each data point in the satellite into pixel point positions consistent with the radar resolution according to the conversion algorithm, thereby realizing consistent conversion of the data format.
In one embodiment, calculating the echo peak height from the weather radar data comprises: acquiring the reflectivity of each data point in the meteorological radar data and echoes corresponding to the reflectivities; extracting the highest elevation angle corresponding to the data point of the echo greater than the preset threshold value; and calculating the echo top height corresponding to the highest elevation angle of the data point according to a height measurement formula.
Specifically, the method for calculating the echo top height comprises the following steps: the server reads the meteorological radar base data, firstly determines the highest elevation angle at which the reflectivity of each data point is greater than a certain threshold echo, and then calculates the height of the grid point at the highest elevation angle layer determined in the last step by using a height measurement formula. For example, the computer device may calculate the echo height from the reflectivities in the radar-based data, and may determine the highest elevation angle at which the echo with a reflectivity greater than 18.3dBz is located at each point in the polar coordinates, e.g., in all elevation layers at a, the reflectivities of the third and fifth elevation layers are greater than 18.3dBz, and the reflectivities of the other elevation layers are less than 18.3dBz, then the highest elevation angle, which is the fifth elevation angle, may calculate the height of the store at the fifth elevation layer, i.e., the echo height at point a, using the height measurement formula. The height measurement formula is shown as formula (6).
h＝h0+Bn×sin(EImax)+Bn2/Rm (6)
Wherein h is the echo top height, h0 is the station height of the radar station, Bn is the distance between the point A and the station, EImax is the highest elevation angle of the point A with the reflectivity larger than 18.3dBz, and Rm is the radius of the earth.
As shown in fig. 3, a flow chart of a method for verifying accuracy of satellite data inversion cloud top height is provided. The method comprises the following steps:
Specifically, the server acquires weather radar data and weather satellite data.
And step 320, respectively reading the screened meteorological radar data and meteorological satellite data.
And the server screens out data corresponding to the same time and the same range from the meteorological radar data and the meteorological satellite data. Furthermore, reflectivity data are read from weather radar data, and brightness temperature value data corresponding to the preset channel are read from weather satellite data.
And step 330, calculating the echo peak height by using the screened meteorological radar data.
Specifically, the server calculates the echo peak height from the reflectivity data in the weather radar data.
Specifically, the cloud top height calculated by using the satellite data and the echo top height calculated by using the radar data have correlation, and a conversion relation is generated according to the correlation, so that the echo top height can be converted into the cloud top height by using the conversion relation.
And step 340, calculating the cloud top height by using the satellite data.
Specifically, the server searches for the cloud top height according to a statistical relationship generated in advance, such as table data.
Data calculation can be performed only between data with the same data unit, so that the cloud top height unit needs to be converted to be consistent with the data format in the radar data.
Step 350, comparing and calculating average deviation.
And calculating the deviation between the cloud top height and the echo top height with the consistent data format by using a deviation calculation formula.
And step 360, obtaining a conclusion.
And evaluating the accuracy of the cloud top height inverted according to the meteorological satellite data according to the numerical value corresponding to the deviation.
In one embodiment, as shown in fig. 4, a schematic flow chart of searching for cloud top height corresponding to weather satellite data in a pre-generated association relationship is provided, which includes:
and step 410, extracting brightness temperature values corresponding to the two preset channels from the meteorological satellite data.
The brightness temperature value is a representation of temperature characteristics, and is a comprehensive result of atmospheric temperature characteristics of various levels from the ground to the upper air. Specifically, the server may obtain the brightness temperature value corresponding to the preset channel in the weather satellite data, and in one embodiment, the channel corresponding to the center wavelength of about 11 μm and 12 μm is selected, so that the brightness temperature values corresponding to the 14 th channel (B14) and the 15 th channel (B15) are selected, and in other embodiments, the brightness temperature values corresponding to other channels may also be obtained, which is not limited in this application.
And step 420, calculating the difference value of the brightness temperature values in the two channels to obtain the brightness temperature difference.
Specifically, the geographic range covered by each channel data in the two channels acquired by the server is the same, such as range data that may be a hemisphere. And then the server calculates the difference value of the obtained brightness temperature values at the corresponding positions in the two channels to obtain the brightness temperature difference corresponding to each position. Further, the server can also perform association binding on the brightness temperature value of the preset channel and the brightness temperature difference obtained by calculating the difference value.
For example, the server can calculate the difference between B14 and B15, obtain the bright temperature difference D of the two channels, and associate and bind B14, B15 and D.
And 430, acquiring the cloud top height corresponding to the brightness temperature value and/or the brightness temperature difference according to a pre-generated incidence relation, wherein the incidence relation is specifically obtained according to the relation between the brightness temperature value and/or the brightness temperature difference and the satellite cloud top height in the cloud top height distribution diagram.
The cloud top height distribution diagram stores the corresponding relation between the brightness temperature value of a certain channel and the brightness temperature difference, wherein the brightness temperature difference is calculated according to the brightness temperature value of the channel and other channels. The cloud top height corresponding to each brightness temperature value and/or brightness temperature difference can be obtained from the cloud top height distribution diagram, and the method comprises the following steps: and acquiring the cloud top height corresponding to each brightness temperature value from the cloud top height distribution map, or acquiring the cloud top height corresponding to each brightness temperature difference from the cloud top height distribution map, and acquiring the cloud top height corresponding to each data pair consisting of each brightness temperature value and each brightness temperature difference from the cloud top height distribution map.
As shown in fig. 5, a schematic diagram of a cloud top height profile is provided. The cloud top height distribution graph can be arranged in a coordinate system, the abscissa corresponds to the brightness temperature value of a certain channel, the ordinate corresponds to the brightness temperature difference between the brightness temperature value of the channel and the brightness temperature values of other channels, each coordinate pair consisting of the brightness temperature value and the brightness temperature difference in the coordinate system also corresponds to a cloud top height, and the cloud top height corresponding to the pixel value can be obtained by reading the pixel value corresponding to the coordinate pair and searching the relation between the pixel value stored in advance and the cloud top height. Specifically, in the cloud top height distribution diagram, the abscissa is the brightness temperature value corresponding to the 11 μm wavelength channel (the 11 μm wavelength channel corresponds to 14 channels), and the ordinate is the brightness temperature difference, which is the difference between the brightness temperature values of the 11 μm wavelength channel and the 12 μm wavelength channel (the 12 μm wavelength channel corresponds to 15 channels). More specifically, in the cloud top height distribution diagram, the range of the bright temperature value of the abscissa may be 180K to 310K, and the range of the bright temperature difference of the ordinate may be-3K to 8K, so that the abscissa may be divided into 1300 parts, i.e., each part represents 0.1K, and the ordinate may be divided into 110 parts, each part represents 0.01K.
Further, the server performs data processing on the acquired cloud top height distribution map to obtain an association relation. Wherein, the incidence relation stores the corresponding relation between the brightness temperature value and/or the brightness temperature difference and the cloud top height, and comprises: the correlation relationship at least stores the corresponding relationship between one of the bright temperature value or the bright temperature difference and the cloud top height. Specifically, the method can comprise the following steps: the incidence relation stores the corresponding relation between the brightness temperature value and the cloud top height, or the incidence relation stores the corresponding relation between the brightness temperature difference and the cloud top height, or the incidence relation stores the corresponding relation between the brightness temperature value and the brightness temperature difference and the cloud top height. Further, the server may store association relations corresponding to different channels in advance, for example, the association relations between two channels B14 and B15, the association relations between two channels B16 and B17, and the like are stored in advance.
The incidence relation obtained according to the cloud top height distribution diagram is the corresponding relation between the brightness temperature value and the brightness temperature difference of different channel data of the satellite and the cloud top height, namely the corresponding relation between the brightness temperature value and the cloud top height and the corresponding relation between the brightness temperature difference and the cloud top height, and is a lookup table, so that the corresponding cloud top height can be obtained by reading the brightness temperature value and the brightness temperature difference in satellite data and then by looking up the incidence relation in the lookup table.
In this embodiment, the cloud top height can be obtained by directly reading a pre-stored association relation, such as an excel table, for the calculation of the cloud top height, and the efficiency of obtaining the cloud top height is greatly improved.
In one embodiment, the method for generating the association relationship includes:
and acquiring a cloud top height distribution map corresponding to a preset channel in the meteorological satellite data. Specifically, one channel corresponds to one cloud top height distribution map, and the server acquires the cloud top height distribution map corresponding to the channel according to the received channel.
And reading the color value corresponding to the data pair consisting of each bright temperature value and each bright temperature difference from the cloud top height distribution map. Specifically, the server may read a color value corresponding to a bright temperature value, or a color value corresponding to a bright temperature difference, or a color value corresponding to a data pair composed of each bright temperature value and each bright temperature difference from the cloud top height distribution map. And after reading the cloud top height distribution map, the server acquires the pixel value, such as the RGB value, of each pixel point in the cloud top height distribution map.
And searching the cloud top height corresponding to the color value. The server obtains a pre-stored relation table between the color value and the cloud top height, and searches the cloud top height corresponding to the color value from the relation table. For example, the color value is converted into the cloud top height by the server according to the read color value of each pixel point, for example, the cloud top height with the red color of 16km, the cloud top height with the orange color of 15km and the like.
And establishing an association relation among the bright temperature value, the bright temperature difference and the cloud top height. And the server establishes an association relation between one of the brightness temperature value or the brightness temperature difference and the cloud top height or an association relation between the brightness temperature value, the brightness temperature difference and the cloud top height according to the acquired brightness temperature value, the brightness temperature difference and the cloud top height. Specifically, the server outputs the numbers corresponding to the abscissa and the ordinate and the corresponding values of the cloud top height to an excel table to form a lookup table. The calculation of the cloud top height only needs to directly read the prestored excel table.
In this embodiment, the acquired cloud top height distribution maps of different channel data are processed to obtain the incidence relation among the brightness temperature value, the brightness temperature difference and the cloud top height, so that the subsequent cloud top height calculation process only needs to search the pre-stored incidence relation, the cloud top height distribution map does not need to be acquired again, the cloud top height data can be acquired more intuitively and accurately by searching the incidence relation, and the efficiency and the accuracy of the cloud top height calculation are improved.
In one embodiment, considering that some undesirable lines, colors or other information may exist in the acquired cloud top height distribution map, the undesirable data needs to be removed, for example, black lines, such as grid lines, coordinate scale values, dotted lines, solid lines, numbers and the like, acquired from the cloud top height distribution map are removed, so as to ensure that the bright temperature value and the color value corresponding to the bright temperature difference can be accurately acquired from the cloud top height distribution map, and further, the cloud top height data can be accurately acquired according to the color value.
And preprocessing the acquired cloud top height distribution map before reading the corresponding color value of the data pair consisting of each bright temperature value and each bright temperature difference from the cloud top height distribution map. Specifically, the step of pre-treating may comprise: the server reads pixel values corresponding to data pairs consisting of the brightness temperature values and the brightness temperature differences from the cloud top height distribution diagram, wherein the pixel values can be RGB values.
When the pixel value corresponds to an achromatic color value, such as a gray value, the server obtains a neighborhood pixel value corresponding to an adjacent region of the data pair, for example, the neighborhood pixel value may be a pixel value corresponding to a previous pixel point of a current pixel, or an average value of pixel values corresponding to three neighborhood regions, where the three neighborhood regions refer to regions composed of 3 × 3 neighborhood pixels. Or the mean value of pixel values corresponding to pixel points in other neighborhood ranges, etc., which is not limited herein.
In the cloud top height distribution map, the pixel values of the data pairs are assigned according to the neighborhood pixel values, and the cloud top height distribution map after reassignment is generated, so that the assignment of the achromatic pixel values in the cloud top height distribution map is realized, the extraction of the cloud top height information cannot be influenced by the unexpected achromatic values, and the incidence relation of the preprocessed cloud top height distribution map can be obtained.
In this embodiment, the current pixel value is assigned by using the neighborhood pixel value, so that the pixel value of the unexpected pixel point in the cloud top height distribution map is processed, the subsequent accurate data reading from the cloud top height distribution map is facilitated, and the efficiency and the accuracy of data reading are improved.
In one embodiment, when the pixel value corresponds to an achromatic value, obtaining a neighborhood pixel value corresponding to an adjacent area of the data pair, and assigning, in the cloud-top height distribution map, the pixel value of the data pair according to the neighborhood pixel value, includes: and carrying out image segmentation on the cloud top height distribution diagram to obtain a colored area and an achromatic area.
The brightness temperature value data corresponding to the colored area has cloud top height information, and the brightness temperature value data corresponding to the achromatic area does not have cloud top height information. Specifically, the server identifies the pixel value of each pixel point in the acquired cloud top height distribution map to obtain a boundary between a colored region and an achromatic region, and divides the cloud top height distribution map according to the boundary to obtain the colored region and the achromatic region. In other embodiments, the implementation manner of the algorithm for segmenting the cloud top height distribution map by the server is not limited.
And when the pixel value corresponding to the data pair in the color region is an achromatic value, acquiring a neighborhood pixel value corresponding to the data pair in the adjacent region, and when the neighborhood pixel value is the chromatic value, assigning a value to the data pair by using the chromatic value.
Specifically, in the color region, the server identifies achromatic color pixels in the color region, and assigns values by using pixel values corresponding to neighborhood pixels of the achromatic color pixels, specifically, when the pixel values of the neighborhood pixels are chromatic, assigns values by using chromatic values, and re-assigns the achromatic color pixels in the color region, so that the achromatic color pixels in the color region are free of the achromatic color pixels. Specifically, achromatic pixels such as solid lines and dotted lines in the chromatic region are assigned to chromatic RGB values of adjacent pixels.
And when the pixel value corresponding to the data in the achromatic color area is a black value in the achromatic color, assigning a value to the data pair by using the white value. And realizing the white processing of the pixel points in the achromatic region.
In this embodiment, the RGB values of the black lines in the cloud top height distribution map are all assigned to be consistent with the surrounding colors, for example, the scale lines and the grid lines in the achromatic regions are both assigned to be white RGB values, and the scale lines and the grid lines in the chromatic regions are both assigned to be chromatic RGB values.
In one embodiment, the server may further perform a filtering process on the acquired cloud top height distribution map to remove undesired pixel information, for example, remove an undesired small region pixel value in the cloud top height distribution map by using a dilation and erosion algorithm, highlight the pixel value information of the desired region by a sharpening process, and the like, which is not limited herein. Moreover, after the acquired cloud top height distribution map is preprocessed to obtain the incidence relation, the subsequent cloud top height calculation only needs to directly read the incidence relation such as an excel table, and the cloud top height distribution map does not need to be preprocessed again, so that the efficiency of calculating the cloud top height is improved.
In one embodiment, finding the cloud top height corresponding to the color value comprises: and acquiring a color target card, and extracting the cloud top height corresponding to the color value from the color target card. The color scale card represents a corresponding relationship between the cloud top height and the color value, for example, the cloud top height values corresponding to the vertical solid lines of 17, 16, 15 and 14 … … are sequentially marked from left to right, the unit km is obtained, different cloud top heights correspond to one color value respectively, and the one-to-one correspondence between the color values and the cloud top heights is realized. Therefore, the server can search the cloud top height corresponding to the color value by using the color code card.
Further, the server searches a color difference value corresponding to the color value in the cloud top height distribution graph. The color value in the color target card is a discontinuous value, and the discontinuous value is used for representing the possible error of the cloud top height. Therefore, the server is further configured to obtain cloud top height chromatic aberration corresponding to each brightness temperature value, where the chromatic aberration may be a numerical value obtained through experimental statistics. Specifically, the color difference may be marked on the cloud top height distribution map, corresponding to the bright temperature value and the bright temperature difference, for example, the dotted lines marked with 0.5, 1.0, 2.0, etc. may represent the variance on the cloud top height distribution map.
And adjusting the searched cloud top height by using the color difference value to obtain the cloud top height corresponding to the color value. If the color difference can be added or subtracted on the searched cloud top height value, the cloud top height value is adjusted to obtain the cloud top height corresponding to the color value.
In this embodiment, the cloud top height is adjusted by using the color difference value, and further, the accuracy of acquiring the cloud top height is improved.
It should be understood that although the various steps in the flow charts of fig. 2-4 are shown in order as indicated by the arrows, the steps are not necessarily performed in order as indicated by the arrows. The steps are not performed in the exact order shown and described, and may be performed in other orders, unless explicitly stated otherwise. Moreover, at least some of the steps in fig. 2-4 may include multiple sub-steps or multiple stages that are not necessarily performed at the same time, but may be performed at different times, and the order of performance of the sub-steps or stages is not necessarily sequential, but may be performed in turn or alternating with other steps or at least some of the sub-steps or stages of other steps.
In one embodiment, as shown in fig. 6, there is provided a cloud top height verification apparatus, including:
the acquiring module 610 is used for acquiring weather radar data and weather satellite data.
And an echo top height calculating module 620, configured to calculate an echo top height according to the weather radar data.
And a cloud top height searching module 630, configured to search a cloud top height corresponding to the weather satellite data from a pre-generated association relationship.
And the verification module 640 is used for verifying the cloud top height by using the echo top height.
In one embodiment, the verification module 640 includes:
and the conversion unit is used for converting the data types of the cloud top height and the echo top height corresponding to the same area range into consistent data types by using a conversion algorithm.
And the verification unit is used for calculating the absolute value deviation between the echo top height and the cloud top height with the consistent data types and verifying the cloud top height according to the absolute value deviation.
In one embodiment, the conversion unit includes:
and the acquisition subunit is used for acquiring the cloud top height and the echo top height corresponding to the same area range.
And the conversion subunit is used for converting the longitude and latitude numerical value corresponding to the cloud top height into a distance numerical value according to a pre-generated conversion algorithm, and the conversion algorithm is obtained according to the longitude and latitude numerical value corresponding to the position of the radar station, the earth radius and the equator radius.
In one embodiment, the cloud top height finding module 630 includes:
and the brightness temperature value extraction unit is used for extracting brightness temperature values corresponding to the two preset channels from the meteorological satellite data.
And the bright temperature difference calculation unit is used for calculating the difference value of the bright temperature values in the two channels to obtain the bright temperature difference.
And the cloud top height calculation unit is used for acquiring the cloud top height corresponding to the brightness temperature value and/or the brightness temperature difference according to a pre-generated incidence relation, and the incidence relation is obtained according to the relation between the brightness temperature value and/or the brightness temperature difference in the cloud top height distribution diagram and the satellite cloud top height.
In one embodiment, the apparatus further comprises:
and the height map acquisition module is used for acquiring a cloud top height distribution map corresponding to a preset channel in the meteorological satellite data.
And the color value reading module is used for reading the color value corresponding to the data pair consisting of each bright temperature value and each bright temperature difference from the cloud top height distribution map.
And the searching module is used for searching the cloud top height corresponding to the color value.
And the relation establishing module is used for establishing the incidence relation among the bright temperature value, the bright temperature difference and the cloud top height.
In one embodiment, the apparatus further comprises:
and the pixel value reading module is used for reading the pixel value corresponding to the data pair consisting of each brightness temperature value and the brightness temperature difference from the cloud top height distribution map.
And the assignment module is used for acquiring neighborhood pixel values corresponding to adjacent regions of the data pairs when the pixel values correspond to the achromatic values, assigning the pixel values of the data pairs according to the neighborhood pixel values in the cloud top height distribution map, and generating the cloud top height distribution map after reassignment.
In one embodiment, the echo top height calculating module 620 includes:
and the echo acquisition unit is used for acquiring the reflectivity of each data point in the meteorological radar data and the echo corresponding to each reflectivity.
And the elevation angle determining unit is used for extracting the highest elevation angle corresponding to the data point larger than the preset threshold echo.
And the calculating unit is used for calculating the echo peak height corresponding to the data point at the highest elevation according to the height measurement formula.
For specific limitations of the cloud top height verification apparatus, reference may be made to the above limitations of the cloud top height verification method, which are not described herein again. The modules in the cloud top height verification device can be wholly or partially realized by software, hardware and a combination thereof. The modules can be embedded in a hardware form or independent from a processor in the computer device, and can also be stored in a memory in the computer device in a software form, so that the processor can call and execute operations corresponding to the modules.
In one embodiment, a computer device is provided, which may be a server, and its internal structure diagram may be as shown in fig. 7. The computer device includes a processor, a memory, a network interface, and a database connected by a system bus. Wherein the processor of the computer device is configured to provide computing and control capabilities. The memory of the computer device comprises a nonvolatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program, and a database. The internal memory provides an environment for the operation of an operating system and computer programs in the non-volatile storage medium. The database of the computer device is used for storing weather radar data and weather satellite data. The network interface of the computer device is used for communicating with an external terminal through a network connection. The computer program is executed by a processor to implement a cloud top height verification method.
Those skilled in the art will appreciate that the architecture shown in fig. 7 is merely a block diagram of some of the structures associated with the disclosed aspects and is not intended to limit the computing devices to which the disclosed aspects apply, as particular computing devices may include more or less components than those shown, or may combine certain components, or have a different arrangement of components.
In one embodiment, a computer device is provided, comprising a memory storing a computer program and a processor implementing the following steps when the processor executes the computer program: acquiring meteorological radar data and meteorological satellite data; calculating the echo peak height according to the meteorological radar data; searching the cloud top height corresponding to the meteorological satellite data from the pre-generated incidence relation; and verifying the cloud top height by using the echo top height.
In one embodiment, the processor when executing the computer program further performs the step of verifying the cloud ceiling height using the echo ceiling height is further configured to: converting the data types of the cloud top height and the echo top height corresponding to the same area range into consistent data types by using a conversion algorithm; and calculating the absolute value deviation between the echo top height and the cloud top height with the consistent data types, and verifying the cloud top height according to the absolute value deviation.
In one embodiment, the step of converting the data types of the cloud top height and the echo top height corresponding to the same region range into the consistent data type by using the conversion algorithm is further performed when the processor executes the computer program: acquiring cloud top height and echo top height corresponding to the same area range; and converting the longitude and latitude numerical value corresponding to the cloud top height into a distance numerical value according to a pre-generated conversion algorithm, wherein the conversion algorithm is specifically obtained according to the longitude and latitude numerical value corresponding to the position of the radar station, the earth radius and the equator radius.
In one embodiment, the processor, when executing the computer program, further performs the step of searching the cloud top height corresponding to the weather satellite data from the pre-generated association relationship, and is configured to: extracting brightness temperature values corresponding to two preset channels from meteorological satellite data; calculating the difference value of the brightness temperature values in the two channels to obtain the brightness temperature difference; and acquiring the cloud top height corresponding to the brightness temperature value and/or the brightness temperature difference according to a pre-generated incidence relation, wherein the incidence relation is specifically obtained according to the relation between the brightness temperature value and/or the brightness temperature difference and the satellite cloud top height in the cloud top height distribution diagram.
In one embodiment, the processor, when executing the computer program, is further configured to: acquiring a cloud top height distribution map corresponding to a preset channel in meteorological satellite data; reading a color value corresponding to a data pair consisting of each bright temperature value and each bright temperature difference from the cloud top height distribution map; searching the cloud top height corresponding to the color value; and establishing an association relation among the bright temperature value, the bright temperature difference and the cloud top height.
In one embodiment, the step of reading the data pair consisting of each bright temperature value and bright temperature difference from the cloud top height distribution map before the processor executes the computer program is further configured to: reading the corresponding pixel values of the data pairs consisting of the brightness temperature values and the brightness temperature differences from the cloud top height distribution map; and when the pixel value corresponds to the achromatic value, acquiring a neighborhood pixel value corresponding to an adjacent region of the data pair, and assigning the pixel value of the data pair according to the neighborhood pixel value in the cloud top height distribution map to generate a cloud top height distribution map after reassignment.
In one embodiment, the processor when executing the computer program further performs the step of calculating the echo top height from the weather radar data by: acquiring the reflectivity of each data point in the meteorological radar data and echoes corresponding to the reflectivities; extracting the highest elevation angle corresponding to the data point of the echo greater than the preset threshold value; and calculating the echo top height corresponding to the highest elevation angle of the data point according to a height measurement formula.
In one embodiment, a computer-readable storage medium is provided, on which a computer program is stored, which when executed by a processor performs the steps of: acquiring meteorological radar data and meteorological satellite data; calculating the echo peak height according to the meteorological radar data; searching the cloud top height corresponding to the meteorological satellite data from the pre-generated incidence relation; and verifying the cloud top height by using the echo top height.
In one embodiment, the computer program when executed by the processor further performs the step of verifying the cloud ceiling height using the echo ceiling height is further configured to: converting the data types of the cloud top height and the echo top height corresponding to the same area range into consistent data types by using a conversion algorithm; and calculating the absolute value deviation between the echo top height and the cloud top height with the consistent data types, and verifying the cloud top height according to the absolute value deviation.
In one embodiment, the computer program when executed by the processor further performs the step of converting the data types of the cloud top height and the echo top height corresponding to the same region range into the consistent data type by using a conversion algorithm, further configured to: acquiring cloud top height and echo top height corresponding to the same area range; and converting the longitude and latitude numerical value corresponding to the cloud top height into a distance numerical value according to a pre-generated conversion algorithm, wherein the conversion algorithm is specifically obtained according to the longitude and latitude numerical value corresponding to the position of the radar station, the earth radius and the equator radius.
In one embodiment, when being executed by the processor, the computer program further performs the step of searching the cloud top height corresponding to the weather satellite data from the pre-generated association relationship, and is further configured to: extracting brightness temperature values corresponding to two preset channels from meteorological satellite data; calculating the difference value of the brightness temperature values in the two channels to obtain the brightness temperature difference; and acquiring the cloud top height corresponding to the brightness temperature value and/or the brightness temperature difference according to a pre-generated incidence relation, wherein the incidence relation is obtained according to the relation between the brightness temperature value and/or the brightness temperature difference in the cloud top height distribution diagram and the satellite cloud top height.
In one embodiment, the computer program when executed by the processor is further operable to: acquiring a cloud top height distribution map corresponding to a preset channel in meteorological satellite data; reading a color value corresponding to a data pair consisting of each bright temperature value and each bright temperature difference from the cloud top height distribution map; searching the cloud top height corresponding to the color value; and establishing an association relation among the bright temperature value, the bright temperature difference and the cloud top height.
In one embodiment, the computer program when executed by the processor performs the steps of reading the color values corresponding to the data pairs consisting of the bright temperature values and the bright temperature differences from the cloud top height distribution map, and further: reading the corresponding pixel values of the data pairs consisting of the brightness temperature values and the brightness temperature differences from the cloud top height distribution map; and when the pixel value corresponds to the achromatic value, acquiring a neighborhood pixel value corresponding to an adjacent region of the data pair, and assigning the pixel value of the data pair according to the neighborhood pixel value in the cloud top height distribution map to generate a cloud top height distribution map after reassignment.
In one embodiment, the computer program when executed by the processor performs the step of calculating the echo height from the weather radar data is further operable to: acquiring the reflectivity of each data point in the meteorological radar data and echoes corresponding to the reflectivities; extracting the highest elevation angle corresponding to the data point of the echo greater than the preset threshold value; and calculating the echo top height corresponding to the highest elevation angle of the data point according to a height measurement formula.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by hardware related to instructions of a computer program, which can be stored in a non-volatile computer-readable storage medium, and when executed, can include the processes of the embodiments of the methods described above. Any reference to memory, storage, database, or other medium used in the embodiments provided herein may include non-volatile and/or volatile memory, among others. Non-volatile memory can include read-only memory (ROM), Programmable ROM (PROM), Electrically Programmable ROM (EPROM), Electrically Erasable Programmable ROM (EEPROM), or flash memory. Volatile memory can include Random Access Memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in a variety of forms such as Static RAM (SRAM), Dynamic RAM (DRAM), Synchronous DRAM (SDRAM), Double Data Rate SDRAM (DDRSDRAM), Enhanced SDRAM (ESDRAM), Synchronous Link DRAM (SLDRAM), Rambus Direct RAM (RDRAM), direct bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM).
The technical features of the above embodiments can be arbitrarily combined, and for the sake of brevity, all possible combinations of the technical features in the above embodiments are not described, but should be considered as the scope of the present specification as long as there is no contradiction between the combinations of the technical features.
The above examples only express several embodiments of the present application, and the description thereof is more specific and detailed, but not construed as limiting the scope of the invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the concept of the present application, which falls within the scope of protection of the present application. Therefore, the protection scope of the present patent shall be subject to the appended claims.
Claims (10)
1. A cloud top height verification method, the method comprising:
acquiring meteorological radar data and meteorological satellite data;
calculating the echo peak height according to the meteorological radar data;
searching the cloud top height corresponding to the meteorological satellite data from a pre-generated incidence relation;
and verifying the cloud top height by using the echo top height.
2. The method of claim 1, wherein the verifying the cloud ceiling height using the echo ceiling height comprises:
converting the data types of the cloud top height and the echo top height corresponding to the same area range into consistent data types by using a conversion algorithm;
and calculating the absolute value deviation between the echo top height and the cloud top height with the consistent data type, and verifying the cloud top height according to the absolute value deviation.
3. The method of claim 2, wherein the converting the data types of the cloud top height and the echo top height corresponding to the same region range into the consistent data type by using a conversion algorithm comprises:
acquiring the cloud top height and the echo top height corresponding to the same area range;
and converting the longitude and latitude numerical value corresponding to the cloud top height into a distance numerical value according to a pre-generated conversion algorithm, wherein the conversion algorithm is specifically obtained according to the longitude and latitude numerical value corresponding to the position of the radar station, the earth radius and the equator radius.
4. The method of claim 1, wherein said calculating an echo ceiling height from said weather radar data comprises:
acquiring the reflectivity of each data point in the meteorological radar data and the echo corresponding to each reflectivity;
extracting the highest elevation angle corresponding to the data point of the echo greater than the preset threshold value;
and calculating the echo top height of the data point corresponding to the highest elevation according to a height measurement formula.
5. The method as claimed in claim 1, wherein the searching for the cloud top height corresponding to the weather satellite data from the pre-generated association relationship comprises:
extracting brightness temperature values corresponding to two preset channels from the meteorological satellite data;
calculating the difference value of the brightness temperature values in the two channels to obtain the brightness temperature difference;
and acquiring the cloud top height corresponding to the brightness temperature value and/or the brightness temperature difference according to a pre-generated incidence relation, wherein the incidence relation is specifically obtained according to the relation between the brightness temperature value and/or the brightness temperature difference and the satellite cloud top height in a cloud top height distribution diagram.
6. The method according to claim 5, wherein the generation method of the association relationship comprises:
acquiring a cloud top height distribution map corresponding to a preset channel in meteorological satellite data;
reading a color value corresponding to a data pair consisting of each brightness temperature value and the brightness temperature difference from the cloud top height distribution map;
searching the cloud top height corresponding to the color value;
and establishing an association relation among the brightness temperature value, the brightness temperature difference and the cloud top height.
7. The method according to claim 6, wherein before reading the color value corresponding to the data pair consisting of each of the bright temperature values and the bright temperature differences from the cloud top height distribution map, further comprising:
reading the corresponding pixel value of the data pair consisting of each bright temperature value and the bright temperature difference from the cloud top height distribution map;
and when the pixel value corresponds to an achromatic value, acquiring a neighborhood pixel value corresponding to an adjacent region of the data pair, and assigning the pixel value of the data pair according to the neighborhood pixel value in the cloud top height distribution map to generate a cloud top height distribution map after reassignment.
8. A cloud top height verification apparatus, the apparatus comprising:
the acquisition module is used for acquiring meteorological radar data and meteorological satellite data;
the echo top height calculating module is used for calculating the echo top height according to the meteorological radar data;
the cloud top height searching module is used for searching the cloud top height corresponding to the meteorological satellite data from a pre-generated incidence relation;
and the verification module is used for verifying the cloud top height by utilizing the echo top height.
9. A computer device comprising a memory and a processor, the memory storing a computer program, wherein the processor implements the steps of the method of any one of claims 1 to 7 when executing the computer program.
10. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of the method of any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010156675.0A CN111366930A (en) | 2020-03-09 | 2020-03-09 | Cloud top height verification method and device, computer equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010156675.0A CN111366930A (en) | 2020-03-09 | 2020-03-09 | Cloud top height verification method and device, computer equipment and storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111366930A true CN111366930A (en) | 2020-07-03 |
Family
ID=71208604
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010156675.0A Pending CN111366930A (en) | 2020-03-09 | 2020-03-09 | Cloud top height verification method and device, computer equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111366930A (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113238231A (en) * | 2021-05-07 | 2021-08-10 | 上海眼控科技股份有限公司 | Echo peak height forecasting method, device, equipment and storage medium |
| JP7051145B1 (en) | 2020-11-13 | 2022-04-11 | 国立研究開発法人宇宙航空研究開発機構 | Optical ground station distributed placement evaluation system, evaluation device and optical ground station distributed placement method |
| CN114677579A (en) * | 2022-03-31 | 2022-06-28 | 国家卫星气象中心(国家空间天气监测预警中心) | A real-time calculation method of cloud top height based on satellite infrared detection element parameters |
| CN114778454A (en) * | 2022-03-30 | 2022-07-22 | 中国人民解放军国防科技大学 | Multilayer cloud top height inversion method based on satellite-borne multi-channel radiation imager |
| CN117452369A (en) * | 2023-12-25 | 2024-01-26 | 江西师范大学 | Echo jacking calculation optimization method for short-time disastrous weather monitoring |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101255736B1 (en) * | 2012-09-21 | 2013-04-26 | 대한민국 | Method for classifying meteorological/non-meteorological echoes using single polarization radars |
| CN103901505A (en) * | 2014-03-11 | 2014-07-02 | 中国气象科学研究院 | Cloud top height detection method and device based on wet bulb effect |
| CN108445464A (en) * | 2018-03-12 | 2018-08-24 | 南京恩瑞特实业有限公司 | Satellite radar inverting fusion methods of the NRIET based on machine learning |
| CN110703256A (en) * | 2019-10-24 | 2020-01-17 | 上海眼控科技股份有限公司 | Radar data display method and device, computer equipment and storage medium |
-
2020
- 2020-03-09 CN CN202010156675.0A patent/CN111366930A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101255736B1 (en) * | 2012-09-21 | 2013-04-26 | 대한민국 | Method for classifying meteorological/non-meteorological echoes using single polarization radars |
| CN103901505A (en) * | 2014-03-11 | 2014-07-02 | 中国气象科学研究院 | Cloud top height detection method and device based on wet bulb effect |
| CN108445464A (en) * | 2018-03-12 | 2018-08-24 | 南京恩瑞特实业有限公司 | Satellite radar inverting fusion methods of the NRIET based on machine learning |
| CN110703256A (en) * | 2019-10-24 | 2020-01-17 | 上海眼控科技股份有限公司 | Radar data display method and device, computer equipment and storage medium |
Non-Patent Citations (3)
| Title |
|---|
| 周非非等: "FY-2卫星反演的云顶高度与多普勒雷达回波顶高的关系初探", 《气象》 * |
| 曹亚楠等: "AIRS红外高光谱卫星数据反演卷云光学厚度和云顶高度", 《光谱学与光谱分析》 * |
| 梁玉冰等: "利用分裂窗直方图法对半透明云云顶高度的反演研究", 《北京大学学报(自然科学版)》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP7051145B1 (en) | 2020-11-13 | 2022-04-11 | 国立研究開発法人宇宙航空研究開発機構 | Optical ground station distributed placement evaluation system, evaluation device and optical ground station distributed placement method |
| JP2022078679A (en) * | 2020-11-13 | 2022-05-25 | 国立研究開発法人宇宙航空研究開発機構 | Optical ground station dispersion arrangement evaluation system, evaluation device and optical ground station dispersion arrangement method |
| CN113238231A (en) * | 2021-05-07 | 2021-08-10 | 上海眼控科技股份有限公司 | Echo peak height forecasting method, device, equipment and storage medium |
| CN114778454A (en) * | 2022-03-30 | 2022-07-22 | 中国人民解放军国防科技大学 | Multilayer cloud top height inversion method based on satellite-borne multi-channel radiation imager |
| CN114677579A (en) * | 2022-03-31 | 2022-06-28 | 国家卫星气象中心(国家空间天气监测预警中心) | A real-time calculation method of cloud top height based on satellite infrared detection element parameters |
| CN114677579B (en) * | 2022-03-31 | 2025-04-08 | 国家卫星气象中心(国家空间天气监测预警中心) | Cloud top height real-time calculation method based on satellite infrared detection element parameters |
| CN117452369A (en) * | 2023-12-25 | 2024-01-26 | 江西师范大学 | Echo jacking calculation optimization method for short-time disastrous weather monitoring |
| CN117452369B (en) * | 2023-12-25 | 2024-04-05 | 江西师范大学 | An optimization method for calculating echo top height for short-term severe weather monitoring |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111366930A (en) | Cloud top height verification method and device, computer equipment and storage medium | |
| US4912770A (en) | Method of detecting change using image | |
| CN114088633A (en) | Coal mine area methane emission abnormity identification and accounting method based on satellite-ground cooperative monitoring | |
| Shaoqing et al. | The comparative study of three methods of remote sensing image change detection | |
| CN115410074B (en) | Remote sensing image cloud detection method and device | |
| CN113096129A (en) | Method and device for detecting cloud cover in hyperspectral satellite image | |
| CN110826509A (en) | A system and method for extracting grassland fence information based on high-resolution remote sensing images | |
| CN116824410A (en) | Inversion method and device for vegetation coverage aiming at hyperspectral remote sensing image data | |
| Apke et al. | Comparison of optical flow derivation techniques for retrieving tropospheric winds from satellite image sequences | |
| CN119025927A (en) | Rapid water quality inversion method, device, equipment and storage medium | |
| Korets et al. | Vegetation cover mapping based on remote sensing and digital elevation model data | |
| CN117953388B (en) | Rapid development stream identification method and device based on satellite fusion data | |
| CN114332085A (en) | Optical satellite remote sensing image detection method | |
| CN117470280B (en) | Laser SLAM real-time precision evaluation method, device, medium and equipment | |
| Li et al. | An improved single-channel polar region ice surface temperature retrieval algorithm using landsat-8 data | |
| CN113379810B (en) | Forest land point cloud data registration method, device, equipment and storage medium | |
| CN111223153A (en) | Cloud top height calculation method and device, computer equipment and storage medium | |
| CN111814692B (en) | Method for acquiring influence distance of ground object type on ground surface temperature | |
| CN115393306A (en) | Land utilization dynamic change analysis system for geographic remote sensing information service | |
| CN119247520B (en) | Compact Gao Bingyun detection method and system based on water vapor and ice crystal absorption characteristics | |
| CN111382698A (en) | Cloud-type computing method, apparatus, computer equipment and storage medium | |
| CN117190982B (en) | Geographic information mapping system based on remote sensing image | |
| JP7127927B1 (en) | Water vapor observation method | |
| Zekoll et al. | Comparison of Masking Algorithms for Sentinel-2 Imagery. Remote Sens. 2021, 13, 137 | |
| Hou et al. | A Landmark Matching Algorithm for the Geostationary Satellite Images Based on Multi-Level Grids |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| AD01 | Patent right deemed abandoned |
Effective date of abandoning: 20230324 |
|
| AD01 | Patent right deemed abandoned |