CN114820966A - Method and device for reconstructing three-dimensional cloud through satellite-ground cloud fusion, electronic equipment and medium - Google Patents
Method and device for reconstructing three-dimensional cloud through satellite-ground cloud fusion, electronic equipment and medium Download PDFInfo
- Publication number
- CN114820966A CN114820966A CN202210457626.XA CN202210457626A CN114820966A CN 114820966 A CN114820966 A CN 114820966A CN 202210457626 A CN202210457626 A CN 202210457626A CN 114820966 A CN114820966 A CN 114820966A
- Authority
- CN
- China
- Prior art keywords
- cloud
- data
- satellite
- radar
- height
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
-
- 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/95—Retrieval from the web
- G06F16/953—Querying, e.g. by the use of web search engines
- G06F16/9537—Spatial or temporal dependent retrieval, e.g. spatiotemporal queries
-
- 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/95—Retrieval from the web
- G06F16/953—Querying, e.g. by the use of web search engines
- G06F16/9538—Presentation of query results
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Databases & Information Systems (AREA)
- Theoretical Computer Science (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Data Mining & Analysis (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Environmental & Geological Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Electromagnetism (AREA)
- Computer Graphics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Atmospheric Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Ecology (AREA)
- Environmental Sciences (AREA)
- Image Processing (AREA)
Abstract
The embodiment of the disclosure discloses a method, a device, electronic equipment and a medium for reconstructing a three-dimensional cloud body through satellite-ground cloud fusion, wherein the method for reconstructing the three-dimensional cloud body through satellite-ground cloud fusion comprises the following steps: acquiring foundation cloud observation data of at least three cloud radars; acquiring data of a stationary orbit meteorological satellite multichannel scanning imaging radiometer and data of a cloud detection product; performing space-time matching on the foundation cloud observation data, the multi-channel scanning imaging radiometer data and the cloud detection product data to obtain satellite-ground matching data; based on the satellite-ground matching data, cloud top height inversion of the satellite in a preset area around each cloud radar installation point is achieved, and cloud top height of a cloud radar networking area is obtained; obtaining the cloud base height of a cloud radar networking area based on the foundation cloud observation data; obtaining a three-dimensional cloud body boundary based on the cloud top height and the cloud bottom height; and obtaining a reconstructed three-dimensional cloud body based on the reflectivity factor profile and the radial velocity profile in the foundation cloud observation data.
Description
Technical Field
The disclosure relates to the field of atmospheric detection and atmospheric remote sensing, in particular to a method and a device for reconstructing a three-dimensional cloud body through satellite-ground cloud fusion, electronic equipment and a medium.
Background
One way of obtaining the cloud observation data is to perform large-area remote sensing observation through a meteorological satellite, particularly to realize high-frequency continuous observation on the same place of the ground by keeping the relative position between a static orbit meteorological satellite and the earth unchanged. The meteorological satellite can acquire information such as large-scale cloud distribution, cloud amount and cloud shape, and further invert to obtain cloud top height information based on the satellite infrared cloud chart by utilizing the characteristics of infrared radiation and temperature correlation of the cloud. However, since the satellite observation data is influenced by various factors such as the earth atmosphere, the earth surface, the sun position, the working state of an on-satellite instrument and the like, and the infrared radiation intensity of the cloud is related to factors such as the density of the cloud, the height of the cloud top, the thickness of the cloud body and the like, the cloud top height directly inverted by the satellite sometimes has a large error.
Another way to obtain cloud observation data is through single-point observation of foundation clouds, such as laser ceilometers, all-sky imagers, millimeter wave cloud radars, and the like. The millimeter wave cloud radar works in a Ka (35GHz) waveband, can penetrate through a thicker cloud layer by utilizing the scattering characteristic of a water condensate on millimeter waves, obtains basic information such as echo intensity, radial velocity, velocity spectrum width and the like in a cloud body, and further processes the basic information to obtain cloud top height and cloud bottom height products. However, data obtained by single-point observation of foundation cloud has the characteristics of more 'point' products, less 'surface' products and single vertical product, and the characteristics limit the exertion of comprehensive observation benefits.
Disclosure of Invention
The embodiment of the disclosure provides a method and a device for reconstructing a three-dimensional cloud body through satellite-ground cloud fusion, electronic equipment and a medium.
In a first aspect, a method for reconstructing a three-dimensional cloud body through satellite-ground cloud fusion is provided in the embodiments of the present disclosure.
The method for reconstructing the three-dimensional cloud body through satellite-ground cloud fusion comprises the following steps:
acquiring ground cloud observation data of at least three cloud radars, wherein the at least three cloud radars are networked at preset station distances;
based on the obtained foundation cloud observation data, obtaining the cloud top height, the cloud bottom height, the reflectivity factor profile and the radial speed profile of each cloud radar installation point vertical headspace;
acquiring multi-channel scanning imaging radiometer data and cloud detection product data of a stationary orbit meteorological satellite;
performing space-time matching on the foundation cloud observation data, the multi-channel scanning imaging radiometer data and the cloud detection product data to obtain satellite-ground matching data;
based on the satellite-ground matching data, cloud top height inversion of the satellite in a preset area around each cloud radar installation point is achieved, the cloud top height of the preset area around each cloud radar installation point is obtained, and further the cloud top height of a cloud radar networking area is obtained;
interpolating the cloud base heights of the cloud radars to obtain the cloud base heights of the cloud radar networking areas;
matching the cloud bottom height and the cloud top height of the cloud radar networking area to obtain a three-dimensional cloud body boundary;
performing on-plane interpolation on each vertical height layer based on the cloud radar reflectivity factor profile and the radial velocity profile to fill the three-dimensional cloud body, so as to obtain a reconstructed three-dimensional cloud body, wherein the height of each vertical height layer is determined by the vertical distance resolution of the cloud radar.
According to the embodiment of the disclosure, the obtaining of the cloud top height, the cloud bottom height, the reflectivity factor profile and the radial velocity profile of the vertical headspace of each cloud radar installation point based on the obtained foundation cloud observation data includes:
acquiring cloud radar power spectrum data in the foundation cloud observation data;
obtaining a main peak of a power spectrum according to the position, the width and the amplitude of a peak body in the power spectrum data, and obtaining a reflectivity factor and a radial velocity according to the main peak of the power spectrum;
acquiring power spectrum data on each vertical height layer, and obtaining a reflectivity factor profile and a radial velocity profile in the vertical direction according to the power spectrum data on each vertical height layer;
and obtaining the cloud top height and the cloud bottom height in the vertical direction according to the reflectivity factor profile and the radial speed profile.
According to the embodiment of the present disclosure, the performing space-time matching on the ground-based cloud observation data, the multi-channel scanning imaging radiometer data, and the cloud detection product data to obtain satellite-ground matching data includes:
taking M minutes before and after satellite observation as a foundation cloud observation data time window;
determining a cloud free time fraction within the time window based on the cloud detection product data;
when the cloud-free time ratio in the time window is smaller than a first preset threshold value, searching a pixel point which is closest to a satellite according to the position information of a cloud radar installation point for matching to obtain space matching data;
and respectively carrying out normalization processing on the cloud top height and the cloud bottom height of the cloud radar in the time window, and then matching the data after the normalization processing with the space matching data to obtain the satellite-ground matching data, wherein M is a positive integer.
According to the embodiment of the disclosure, the method for achieving cloud top height inversion of the satellite in the preset area around each cloud radar installation point based on the satellite-ground matching data to obtain the cloud top height of the preset area of each cloud radar installation point includes:
normalizing the cloud top height of the cloud radar in the time window to obtain the normalized cloud top height H of the matched pixel points a ；
Obtaining thermal infrared channel data Z in multi-channel scanning imaging radiometer data matched with image element points a ；
Based on the formula
Obtaining a satellite-ground cloud fusion factor of the cloud radar at the time of observation corresponding to the time window
Based on the formula
Obtaining the cloud top height H of other pixels around the secondary cloud radar mounting point during observation of the satellite ij Wherein i and j are the number of rows and columns of positions of other pixel points of the satellite;
acquiring the cloud top height of the preset area of each cloud radar mounting point based on the cloud top height of the pixels around all the cloud radar mounting points during observation of the satellite;
the preset area is a circular area which takes a cloud radar installation point as a center and takes a first preset distance as a radius.
According to the embodiment of the disclosure, the combination of the cloud top heights of the preset areas of the cloud radar installation points to obtain the cloud top height of the cloud radar networking area comprises the following steps:
and when the preset areas of the cloud radar installation points are overlapped, determining the average value of the cloud top heights of the preset areas of the cloud radar installation points forming the overlap as the cloud top height of the overlap area.
According to the embodiment of the disclosure, when multiple layers of clouds exist in the vertical observation direction of the cloud radar, the cloud top height, the cloud bottom height, the reflectivity factor profile and the radial velocity profile of the vertical headspace of the cloud radar installation point are determined based on the observation data of the topmost cloud with the thickness larger than the second preset threshold.
According to the embodiment of the disclosure, the interpolation refers to interpolation with the satellite spatial resolution as a grid point.
According to the embodiment of the disclosure, the method further comprises the step of taking extreme values of the reflectivity factor and the radial speed of the three-dimensional cloud body to obtain a combined reflectivity factor and a combined radial speed of the cloud radar networking area, wherein the extreme values of the reflectivity factor and the radial speed of the three-dimensional cloud body are taken as follows, the maximum value of the reflectivity factor is taken, and the minimum value of the radial speed is taken.
In a second aspect, an embodiment of the present disclosure provides an apparatus for reconstructing a three-dimensional cloud through satellite-ground cloud fusion.
The device for reconstructing the three-dimensional cloud body through satellite-ground cloud fusion comprises:
the system comprises a first acquisition module, a second acquisition module and a third acquisition module, wherein the first acquisition module is configured to acquire foundation cloud observation data of at least three cloud radars, and based on the acquired foundation cloud observation data, the cloud top height, the cloud bottom height, the reflectivity factor and the radial speed profile of the vertical headspace of each cloud radar installation point are acquired, and the at least three cloud radars are networked at a preset station distance;
a second acquisition module configured to acquire multi-channel scanning imaging radiometer data and cloud detection product data of a stationary orbit meteorological satellite;
a space-time matching module configured to perform space-time matching on the foundation cloud observation data, the multi-channel scanning imaging radiometer data and the cloud detection product data to obtain satellite-ground matching data;
the inversion module is configured to realize cloud top height inversion of the satellite in a preset area around each cloud radar installation point based on the satellite-ground matching data, obtain cloud top heights of the preset areas of the cloud radar installation points, and further obtain cloud top heights of a cloud radar networking area;
the cloud base height determining module is configured to interpolate the cloud base heights of the cloud radars to obtain the cloud base heights of the cloud radar networking areas;
the reconstruction module is configured to match the cloud bottom height and the cloud top height of the cloud radar networking area to obtain a three-dimensional cloud body boundary; performing on-plane interpolation on each vertical height layer based on the cloud radar reflectivity factor profile and the radial velocity profile to fill the three-dimensional cloud body, so as to obtain a reconstructed three-dimensional cloud body, wherein the height of each vertical height layer is determined by the vertical distance resolution of the cloud radar.
In a third aspect, an embodiment of the present disclosure provides an electronic device, including a memory and a processor, where the memory is used to store one or more computer instructions for supporting a data processing apparatus to execute the above method for reconstructing a three-dimensional cloud through satellite-ground cloud fusion, and the processor is configured to execute the computer instructions stored in the memory. The electronic device may also include a communication interface for the electronic device to communicate with other devices or a communication network.
In a fourth aspect, an embodiment of the present invention provides a computer-readable storage medium for storing computer instructions for a data processing apparatus, where the computer instructions include computer instructions for executing the method for reconstructing a three-dimensional cloud through satellite-ground cloud fusion as described above.
The technical scheme provided by the embodiment of the disclosure can have the following beneficial effects:
the method for reconstructing the three-dimensional cloud body through satellite-ground cloud fusion effectively combines the characteristics of large observation scale of the geostationary orbit meteorological satellite and high observation precision of the foundation cloud radar, performs space-time matching and data fusion on the geostationary orbit meteorological satellite and the ground cloud radar, can quickly and accurately realize three-dimensional reconstruction of cloud in a cloud radar networking area, and improves the speed and precision of three-dimensional cloud body reconstruction; the reconstructed three-dimensional cloud body can provide various products such as cloud top height, cloud bottom height and three-dimensional cloud fusion products on the cloud radar networking area surface, and can be widely applied to cloud visual monitoring and application.
It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the disclosure.
Drawings
Other features, objects, and advantages of the present disclosure will become more apparent from the following detailed description of non-limiting embodiments when taken in conjunction with the accompanying drawings. In the drawings:
fig. 1 shows a flowchart of a method for reconstructing a three-dimensional cloud through satellite-ground cloud fusion according to an embodiment of the present disclosure.
Fig. 2 shows a schematic diagram of a cloud radar networking area acquisition process according to an embodiment of the present disclosure.
Fig. 3 shows a structure diagram of a device for reconstructing a three-dimensional cloud through satellite-ground cloud fusion according to an embodiment of the present disclosure.
Fig. 4 shows a block diagram of an electronic device according to an embodiment of the present disclosure.
FIG. 5 is a schematic block diagram of a computer system suitable for use in implementing a data transfer method according to an embodiment of the present disclosure.
Detailed Description
Hereinafter, exemplary embodiments of the present disclosure will be described in detail with reference to the accompanying drawings so that those skilled in the art can easily implement them. Also, for the sake of clarity, parts not relevant to the description of the exemplary embodiments are omitted in the drawings.
In the present disclosure, it is to be understood that terms such as "including" or "having," etc., are intended to indicate the presence of the disclosed features, numbers, steps, actions, components, parts, or combinations thereof, and do not preclude the possibility that one or more other features, numbers, steps, actions, components, parts, or combinations thereof are present or added.
It should be further noted that the embodiments and features of the embodiments in the present disclosure may be combined with each other without conflict. The present disclosure will be described in detail below with reference to the accompanying drawings in conjunction with embodiments.
The details of the embodiments of the present disclosure are described in detail below with reference to specific embodiments.
As mentioned above, since the satellite observation data is influenced by various factors such as the earth atmosphere, the earth surface, the sun position, and the working state of the satellite instrument, and the infrared radiation intensity of the cloud is related to the density of the cloud, the height of the cloud top, the thickness of the cloud, and the like, the error of the cloud top height directly inverted by the satellite is sometimes large. Meanwhile, data obtained by single-point observation of foundation cloud has the characteristics of more 'point' products, less 'surface' products and single vertical product, and the characteristics limit the exertion of comprehensive observation benefits. Therefore, the conventional three-dimensional cloud body reconstruction method cannot rapidly and accurately reconstruct the three-dimensional cloud body in a larger range.
In view of the above, the method for reconstructing the three-dimensional cloud body through satellite-ground cloud fusion is provided, and by effectively combining the characteristics of large observation scale of the geostationary orbit meteorological satellite and high observation precision of the foundation cloud radar, the space-time matching and the data fusion are carried out on the satellite-ground cloud satellite and the foundation cloud radar, so that the three-dimensional stereo reconstruction of the cloud in the cloud radar networking area can be rapidly and accurately realized, and the speed and the precision of the three-dimensional cloud body reconstruction are improved; the reconstructed three-dimensional cloud body can provide various products such as cloud top height, cloud bottom height, combined reflectivity factor, combined radial velocity products, three-dimensional cloud fusion products and the like on the surface of a cloud radar networking area, and can be widely applied to cloud visual monitoring and application.
Fig. 1 shows a flowchart of a method for reconstructing a three-dimensional cloud through satellite-ground cloud fusion according to an embodiment of the present disclosure. As shown in fig. 1, the method for reconstructing a three-dimensional cloud through satellite-ground cloud fusion includes the following steps S101 to S107:
in step S101, acquiring ground-based cloud observation data of at least three cloud radars, wherein the at least three cloud radars are networked at a preset station distance;
in step S102, based on the obtained ground cloud observation data, a cloud top height, a cloud bottom height, a reflectivity factor profile and a radial velocity profile of the vertical headspace of each cloud radar installation point are obtained;
in step S103, acquiring multichannel scanning imaging radiometer data and cloud detection product data of the stationary orbit meteorological satellite;
in step S104, performing space-time matching on the ground-based cloud observation data, the multi-channel scanning imaging radiometer data, and the cloud detection product data to obtain satellite-ground matching data;
in step S105, based on the satellite-ground matching data, cloud top height inversion of the satellite in a preset area around each cloud radar installation point is realized, so as to obtain the cloud top height of the preset area of each cloud radar installation point, and further obtain the cloud top height of a cloud radar networking area;
in step S106, interpolating the cloud base heights of the cloud radars to obtain the cloud base heights of the cloud radar networking areas;
in step S107, matching the cloud bottom height and the cloud top height of the cloud radar networking area to obtain a three-dimensional cloud body boundary;
in step S108, performing an above-surface interpolation on each vertical height layer based on the cloud radar reflectivity factor profile and the radial velocity profile to fill the three-dimensional cloud, so as to obtain a reconstructed three-dimensional cloud, wherein the height of each vertical height layer is determined by a vertical distance resolution of the cloud radar.
In the embodiment of the present disclosure, the cloud radar may be a millimeter wave cloud radar, and the at least three cloud radars may be arranged in a polygonal manner for networking, for example, when the cloud radar is three, a triangular layout is recommended, when the cloud radar is four, a square or rectangular layout is recommended, and when the cloud radar is five, a pentagonal layout is recommended. In addition, when the number of the cloud radars is large, besides forming a polygon, the layout needs to be performed in the middle area of the polygon, so as to realize the overall coverage of the inner area of the polygon. And a certain station distance is formed between the cloud radars, and the certain station distance can be the farthest 100 kilometers of two adjacent stations.
In this disclosure, the obtaining of the cloud top height, the cloud bottom height, the reflectivity factor profile, and the radial velocity profile of the vertical headspace of each cloud radar installation point based on the obtained ground cloud observation data includes: acquiring cloud radar power spectrum data in the foundation cloud observation data; obtaining a main peak of a power spectrum according to the position, the width and the amplitude of a peak body in the power spectrum data, and obtaining a reflectivity factor and a radial velocity according to the main peak of the power spectrum; acquiring power spectrum data on each vertical height layer, and obtaining a reflectivity factor profile and a radial velocity profile in the vertical direction according to the power spectrum data on each vertical height layer; and obtaining the cloud top height and the cloud bottom height in the vertical direction according to the reflectivity factor profile and the radial speed profile. Wherein the height of each vertical height layer is determined by the vertical range resolution of the cloud radar, and may be, for example, one vertical height layer every 30 meters.
In the embodiment of the disclosure, the geostationary orbit meteorological satellite may be a wind cloud four meteorological satellite, and the multi-channel scanning imaging radiometric data includes thermal infrared channel data which is most sensitive to the height of a cloud top in a plurality of channels. The multichannel scanning imaging radiometer data is first-level data observed by a meteorological satellite, is generally represented by pixel brightness values of remote sensing images, records gray values of target objects, is unit-free, and has the size related to the radiation resolution of a sensor, the emissivity of the target objects, the atmospheric transmittance, the scattering rate and the like. In the embodiment of the disclosure, the cloud detection product data is secondary data obtained by processing data of a meteorological satellite, and can represent whether each pixel point observed by the satellite is a cloud point. The cloud detection product can be used for quality control of three-dimensional clouds, for example, if the satellite does not monitor the cloud, the cloud top height and the cloud bottom height of a cloud radar networking area and the three-dimensional clouds should not exist.
In the embodiment of the disclosure, in order to eliminate the uncertainty error introduced by the radiometric calibration conversion, when the geostationary orbit meteorological satellite observes in the cloud radar networking area, the atmospheric angle, the earth surface angle, the solar altitude angle and the working state of the on-satellite instrument in the networking area can be approximately considered to be kept unchanged in a short time, and the only change is the radiation capability of the cloud body, so that the change of satellite channel data on each pixel can directly reflect the change of cloud information when the same satellite in the networking area observes. Furthermore, considering that the thermal infrared channel data is related to the infrared radiation of a target object, the size of the infrared radiation energy directly reflects the height of the cloud top temperature, and the cloud top temperature determines the cloud top height, so that space-time matching is firstly carried out on the cloud top height obtained from the cloud radar observation data and the thermal infrared channel data of the meteorological satellite to obtain satellite-ground matching data, then a satellite-ground cloud fusion factor is obtained based on the satellite-ground matching data, and finally the cloud top height inversion of the satellite in a preset area around each cloud radar installation point is realized based on the satellite-ground cloud fusion factor, so that the cloud top height inversion of the satellite in a cloud radar networking area can be effectively realized.
In an embodiment of the present disclosure, the performing space-time matching on the ground-based cloud observation data, the multi-channel scanning imaging radiometer data, and the cloud detection product data to obtain satellite-ground matching data includes: taking M minutes before and after satellite observation as a foundation cloud observation data time window; determining a cloud free time fraction within the time window based on the cloud detection product data; when the cloud-free time ratio in the time window is smaller than a first preset threshold value, searching a pixel point which is closest to a satellite according to the position information of a cloud radar installation point for matching to obtain space matching data; and respectively carrying out normalization processing on the cloud top height and the cloud bottom height in the time window, and then matching the data after the normalization processing with the space matching data to obtain the satellite-ground matching data, wherein M is a positive integer. And when the cloud-free time ratio in the time window is greater than or equal to the first preset threshold, determining that the time matching fails, and waiting for the next time window without performing space matching.
Specifically, the spatial resolution of the thermal infrared channel of the wind and cloud four meteorological satellite is 4km × 4km, the observation horizontal range of the millimeter wave cloud radar at the height of 10km is hundreds of meters, and the meteorological satellite data is acquired only at the preset observation time, so that the two meteorological satellite data have certain difference in space and time, and space-time matching needs to be performed on the multichannel scanning imaging radiometer data of the meteorological satellite and the ground-based cloud observation data of the cloud radar before data fusion. The space-time matching specifically comprises: firstly, determining a foundation cloud observation data time window according to the observation time of the satellite, for example, M minutes before and after the observation time of the satellite can be used as the foundation cloud observation data time window, wherein M is a positive integer; then determining a cloud-free time fraction based on the cloud detection product data within a determined time window; when the cloud-free time ratio in the time window is smaller than a first preset threshold, searching a satellite and a corresponding pixel point thereof for matching according to cloud radar mounting point position information to obtain space matching data, wherein the first preset threshold can be comprehensively determined according to historical data and actual conditions, and can be 30% for example; and finally, respectively carrying out normalization processing on the cloud top height and the cloud bottom height in the time window, and matching the data after the normalization processing with the space matching data to obtain the satellite-ground matching data.
In the embodiment of the present disclosure, when the cloud radar data time window is determined, the time of satellite observation may be respectively pushed forward and backward for N minutes, and data of 2N +1 minutes is obtained altogether, so as to further expand data samples and improve processing accuracy. In the embodiment of the present disclosure, the normalization process may be performed by taking an average, a mode, a median, or the like.
In this disclosure, the implementing, based on the satellite-ground matching data, cloud top height inversion of the satellite in a preset area around each cloud radar installation point to obtain a cloud top height of the preset area of each cloud radar installation point includes: normalizing the cloud top height in the time window to obtain the normalized cloud top height H of the matched pixel point a (ii) a Thermal infrared channel data Z in meteorological satellite multi-channel scanning imaging radiometer data for acquiring matched image element points a (ii) a Based on the formula
Obtaining a satellite-ground cloud fusion factor of the cloud radar at the time of observation corresponding to the time window
Based on the formula
Obtaining the cloud top height H of other pixels around the secondary cloud radar mounting point during observation of the satellite ij Wherein i and j are the number of rows and columns of positions of other pixel points of the satellite; the cloud top of the preset area of each cloud radar mounting point is obtained based on the cloud top height of the pixels around all the cloud radar mounting points during observation of the satelliteHeight.
In this disclosure, the further obtaining the cloud top height of the cloud radar networking area may be to combine the cloud top heights of the preset areas of the cloud radar installation points to obtain the cloud top height of the cloud radar networking area, and includes: and when the preset areas of the cloud radar installation points are overlapped, determining the average value of the cloud top heights of the preset areas of the cloud radar installation points forming the overlap as the cloud top height of the overlap area. Specifically, the cloud radar networking area is formed by combining preset areas of cloud radar installation points, and therefore the cloud top height of the cloud radar networking area can be obtained by combining the cloud top heights of the preset areas of the cloud radar installation points. However, due to many factors, there may be overlap between the preset areas of the cloud radars, and at this time, the cloud top heights of the preset areas of the cloud radar installation points that form the overlap may be averaged to serve as the cloud top height of the overlap area.
In the embodiment of the disclosure, after the cloud top height of the cloud radar networking area is determined, the cloud bottom height of the vertical headspace of each cloud radar mounting point in the foundation cloud observation data can be interpolated to obtain the cloud bottom height of the cloud radar networking area, and then the cloud bottom height of the cloud radar networking area is matched with the cloud top height to obtain the three-dimensional cloud body boundary.
In the embodiment of the disclosure, when multiple layers of clouds exist in the vertical observation direction of the cloud radar, the cloud top height, the cloud bottom height, the reflectivity factor profile and the radial velocity profile of the vertical headspace of the cloud radar installation point may be determined based on observation data of the topmost cloud with the thickness greater than the second preset threshold. Specifically, sky-ground matching needs to consider the sky state, and generally, the matching effect is the best when the zenith is a large-area layered cloud. Meanwhile, because the existence of multiple layers of clouds in the vertical direction is also a common condition, when multiple layers of clouds exist in the vertical observation direction of the cloud radar, the stationary orbit meteorological satellite can only acquire data of the uppermost layer or the next layer of clouds when the uppermost layer or the layers of clouds are thin, so that the cloud top height, the cloud bottom height, the reflectivity factor profile and the radial velocity profile of the vertical headspace of the cloud radar installation point can be determined based on the observation data of the uppermost layer of clouds with the thickness larger than a second preset threshold value. The second preset threshold may be set as required, for example, 1 km.
In the embodiment of the present disclosure, as shown in fig. 2, the cloud radar networking area may be determined in the following manner: firstly, determining a cloud radar observation area needing to be observed; then networking at least three foundation millimeter wave cloud radars in the cloud radar observation area at a certain station distance, wherein the networking cloud radars are in irregular shape layout with different station distances due to the influence of factors such as equipment erection conditions and the like, and the optimal satellite-ground matching radius can be selected according to the actual station distance of the adjacent cloud radars to obtain a preset area Q1 which can be effectively covered by each cloud radar; secondly, clustering and analyzing the multichannel scanning imaging radiometric data of the meteorological satellite by using a satellite remote sensing technology to obtain an identification result of whether each pixel is a cloud at each observation time, further obtaining cloud detection product data, and identifying in the cloud radar networking observation area based on the cloud detection product data to obtain a cloud range Q2; and finally, acquiring the coverage range Q1 and the overlapping area Q12 of the cloud range Q2, so that the cloud radar networking area can be obtained.
In the embodiment of the present disclosure, the interpolation is performed on each vertical height layer and the cloud base height is interpolated on the basis of the cloud radar reflectivity factor profile and the radial velocity profile, which may be performed by using a satellite spatial resolution as a lattice point.
According to the technical scheme provided by the embodiment of the disclosure, the characteristics of large observation scale of the meteorological satellite and high observation precision of the foundation cloud radar are effectively combined, and the meteorological satellite and the foundation cloud radar are subjected to time-space matching and data fusion, so that the three-dimensional reconstruction of cloud in a cloud radar networking area can be rapidly and accurately realized, and the speed and precision of the three-dimensional cloud reconstruction are improved; the reconstructed three-dimensional cloud body can provide various products such as cloud top height, cloud bottom height and three-dimensional cloud fusion products on the cloud radar networking area surface, and can be widely applied to cloud visual monitoring and application. In this disclosure, after the reconstructed three-dimensional cloud object is obtained, extreme values may be further taken for the reflectivity factor and the radial velocity of the three-dimensional cloud object to obtain a combined reflectivity factor and a combined radial velocity of the cloud radar networking area, where the extreme values taken for the reflectivity factor and the radial velocity of the three-dimensional cloud object are that the maximum value is taken for the reflectivity factor, and the minimum value is taken for the radial velocity.
Fig. 3 is a structural diagram of an apparatus for reconstructing a three-dimensional cloud through satellite-earth cloud fusion according to an embodiment of the present disclosure, and as shown in fig. 3, the apparatus for reconstructing a three-dimensional cloud through satellite-earth cloud fusion includes:
the first acquisition module 301 is configured to acquire foundation cloud observation data of at least three cloud radars, and obtain a cloud top height, a cloud bottom height, a reflectivity factor profile and a radial velocity profile of a vertical headspace of each cloud radar installation point based on the acquired foundation cloud observation data, wherein the at least three cloud radars are networked at a preset station distance;
a second acquisition module 302 configured to acquire multi-channel scanning imaging radiometer data of a stationary orbit meteorological satellite, the multi-channel scanning imaging radiometer data including thermal infrared channel data, and cloud detection products;
a space-time matching module 303 configured to perform space-time matching on the ground-based cloud observation data, the multi-channel scanning imaging radiometer data, and the cloud detection product data to obtain satellite-ground matching data;
the inversion module 304 is configured to implement cloud top height inversion of the satellite in a preset area around each cloud radar installation point based on the satellite-ground matching data, obtain cloud top heights of the preset areas of the cloud radar installation points, and further obtain cloud top heights of the cloud radar networking areas;
a determining module 305 configured to interpolate the cloud base heights of the cloud radars to obtain cloud base heights of the cloud radar networking areas;
a reconstruction module 306 configured to match the cloud bottom height and the cloud top height of the cloud radar networking area to obtain a three-dimensional cloud body boundary; and performing surface interpolation on each vertical height layer based on the cloud radar reflectivity factor profile and the radial velocity profile to fill the three-dimensional cloud body to obtain a reconstructed three-dimensional cloud body, wherein the height of each vertical height layer is determined by the vertical distance resolution of the cloud radar.
In the embodiment of the present disclosure, the cloud radar may be a millimeter wave cloud radar, and the at least three cloud radars may be arranged in a polygonal manner for networking, for example, when the cloud radar is three, a triangular layout is recommended, when the cloud radar is four, a square or rectangular layout is recommended, and when the cloud radar is five, a pentagonal layout is recommended. In addition, when the number of the cloud radars is large, besides forming a polygon, the layout needs to be performed in the middle area of the polygon, so as to realize the overall coverage of the inner area of the polygon. And a certain station distance is formed between the cloud radars, and the certain station distance can be the farthest 100 kilometers of two adjacent stations.
In this disclosure, the obtaining of the cloud top height, the cloud bottom height, the reflectivity factor, and the radial velocity profile of the vertical headspace of each cloud radar installation point based on the obtained ground cloud observation data includes: acquiring cloud radar power spectrum data in the foundation cloud observation data; obtaining a main peak of a power spectrum according to the position, the width and the amplitude of a peak body in the power spectrum data, and obtaining a reflectivity factor and a radial velocity according to the main peak of the power spectrum; acquiring power spectrum data on each vertical height layer, and obtaining a reflectivity factor profile and a radial velocity profile in the vertical direction according to the power spectrum data on each vertical height layer; and obtaining the cloud top height and the cloud bottom height in the vertical direction according to the reflectivity factor profile and the radial speed profile. Wherein the height of each vertical height layer is determined by the vertical range resolution of the cloud radar, and may be, for example, one vertical height layer every 30 meters.
In the embodiment of the disclosure, the geostationary orbit meteorological satellite may be a wind cloud four meteorological satellite, and the multi-channel scanning imaging radiometric data includes thermal infrared channel data which is most sensitive to the height of a cloud top in 14 channels. The multichannel scanning imaging radiometer data is primary data observed by a satellite, is generally represented by pixel brightness values of remote sensing images, records gray values of target objects, is unit-free, and has the size related to the radiation resolution of a sensor, the emissivity of the target objects, the atmospheric transmittance, the scattering rate and the like. In the embodiment of the present disclosure, the cloud detection product data is secondary data obtained by processing data of a meteorological satellite, and can represent whether each pixel point observed by the satellite is a cloud point. The cloud detection product can be used for quality control of three-dimensional clouds, for example, if the satellite does not monitor the cloud, the cloud top height and the cloud bottom height of a cloud radar networking area and the three-dimensional clouds should not exist.
In the embodiment of the disclosure, in order to eliminate the uncertainty error introduced by the radiometric calibration conversion, when the geostationary orbit meteorological satellite observes in the cloud radar networking area, the atmospheric angle, the earth surface angle, the solar altitude angle and the working state of the on-satellite instrument in the networking area can be approximately considered to be kept unchanged in a short time, and the only change is the radiation capability of the cloud body, so that the change of satellite channel data on each pixel can directly reflect the change of cloud information when the same satellite in the networking area observes. Furthermore, considering that the thermal infrared channel data is related to the infrared radiation of a target object, the cloud top temperature determines the cloud top height, and the temperature directly reflects the size of the infrared radiation energy, so that the cloud top height obtained from the cloud radar observation data is subjected to space-time matching with the thermal infrared channel data of the meteorological satellite to obtain satellite-ground matching data, a satellite-ground cloud fusion factor is obtained based on the satellite-ground matching data, and finally the cloud top height of the satellite in a preset area around each cloud radar installation point is inverted based on the satellite-ground cloud fusion factor, so that the inversion of the cloud top height of the satellite in a cloud radar networking area can be effectively realized. The preset area is a circular area with the cloud radar installation point as the center and the first preset distance as the radius, and the first preset distance can be 80 kilometers.
In an embodiment of the present disclosure, the performing space-time matching on the ground-based cloud observation data, the multi-channel scanning imaging radiometer data, and the cloud detection product data to obtain satellite-ground matching data includes: taking M minutes before and after satellite observation as a foundation cloud observation data time window; determining a cloud free time fraction within the time window based on the cloud detection product data; when the cloud-free time ratio in the time window is smaller than a first preset threshold value, searching a pixel point which is closest to a satellite according to the position information of a cloud radar installation point for matching to obtain space matching data; and respectively carrying out normalization processing on the cloud top height and the cloud bottom height in the time window, and then matching the data after the normalization processing with the space matching data to obtain the satellite-ground matching data, wherein M is a positive integer. And when the cloud-free time ratio in the time window is greater than or equal to the first preset threshold, determining that the time matching fails, and waiting for the next time window without performing space matching.
Specifically, the spatial resolution of the thermal infrared channel of the wind and cloud four meteorological satellite is 4km × 4km, the observation horizontal range of the millimeter wave cloud radar at the height of 10km is hundreds of meters, and the meteorological satellite data is acquired only at the preset observation time, so that the two meteorological satellite data have certain difference in space and time, and space-time matching needs to be performed on the multichannel scanning imaging radiometer data of the meteorological satellite and the ground-based cloud observation data of the cloud radar before data fusion. The space-time matching specifically comprises: firstly, determining a foundation cloud observation data time window according to the observation time of the satellite, for example, M minutes before and after the observation time of the satellite can be used as the foundation cloud observation data time window, wherein M is a positive integer; then determining a cloud-free time fraction based on the cloud detection product data within a determined time window; when the cloud-free time ratio in the time window is smaller than a first preset threshold, searching a satellite and a corresponding pixel point thereof for matching according to cloud radar mounting point position information to obtain space matching data, wherein the first preset threshold can be comprehensively determined according to historical data and actual conditions, and can be 30% for example; and finally, respectively carrying out normalization processing on the cloud top height and the cloud bottom height in the time window, and matching the data after the normalization processing with the space matching data to obtain the satellite-ground matching data.
In the embodiment of the present disclosure, when the cloud radar data time window is determined, the time of satellite observation may be respectively pushed forward and backward for N minutes, and data of 2N +1 minutes is obtained altogether, so as to further expand data samples and improve processing accuracy. In the embodiment of the present disclosure, the normalization process may be performed by taking an average, a mode, a median, and the like.
In this disclosure, the obtaining of the cloud top height of the preset area of each cloud radar installation point by implementing the cloud top height inversion of the satellite in the preset area around each cloud radar installation point based on the satellite-ground matching data includes: normalizing the cloud top height in the time window to obtain the normalized cloud top height H of the matched pixel point a (ii) a Thermal infrared channel data Z in meteorological satellite multi-channel scanning imaging radiometer data for acquiring matched image element points a (ii) a Based on the formula
Obtaining a satellite-ground cloud fusion factor of the cloud radar at the time of observation corresponding to the time window
Based on the formula
Obtaining the cloud top height Hij of other pixels around the cloud radar mounting point during observation of the satellite, wherein i and j are the number of rows and columns of the positions of other pixel points of the satellite; and obtaining the cloud top height of the preset area of each cloud radar mounting point based on the cloud top height of the pixels around all the cloud radar mounting points during observation of the satellite.
In this disclosure, the further obtaining the cloud top height of the cloud radar networking area may be to combine the cloud top heights of the preset areas of the cloud radar installation points to obtain the cloud top height of the cloud radar networking area, and includes: and when the preset areas of the cloud radar installation points are overlapped, determining the average value of the cloud top heights of the preset areas of the cloud radar installation points forming the overlap as the cloud top height of the overlap area. Specifically, the cloud radar networking area is formed by combining preset areas of cloud radar installation points, and therefore the cloud top height of the cloud radar networking area can be obtained by combining the cloud top heights of the preset areas of the cloud radar installation points. However, due to many factors, there may be overlap between the preset areas of the cloud radars, and at this time, the cloud top heights of the preset areas of the cloud radar installation points that form the overlap may be averaged to serve as the cloud top height of the overlap area.
In the embodiment of the disclosure, after the cloud top height of the cloud radar networking area is determined, interpolation may be performed on the cloud bottom height of the vertical headspace of each cloud radar mounting point in the foundation cloud observation data to obtain the cloud bottom height of the cloud radar networking area, and then the cloud bottom height of the cloud radar networking area is matched with the cloud top height to obtain the three-dimensional cloud body boundary.
In the embodiment of the disclosure, when multiple layers of clouds exist in the vertical observation direction of the cloud radar, the cloud top height, the cloud bottom height, the reflectivity factor profile and the radial velocity profile of the vertical headspace of the cloud radar installation point may be determined based on observation data of the topmost cloud with the thickness greater than the second preset threshold. Specifically, the sky-ground matching needs to consider the sky state, and generally, the matching effect is the best when the sky top is a layered cloud with a continuous large area. Meanwhile, because the existence of multiple layers of clouds in the vertical direction is also a common condition, when multiple layers of clouds exist in the vertical observation direction of the cloud radar, the meteorological satellite can only acquire data of the uppermost layer or the next layer of clouds when the uppermost layer or the layers of clouds are thin, so that the cloud top height, the cloud bottom height, the reflectivity factor profile and the radial velocity profile of the vertical headspace of the cloud radar installation point can be determined based on the observation data of the uppermost layer of clouds with the thickness larger than a second preset threshold value. The second preset threshold may be set as required, for example, 1 km.
In the embodiment of the present disclosure, the interpolation is performed on each vertical height layer and the cloud base height is interpolated on the basis of the cloud radar reflectivity factor profile and the radial velocity profile, which may be performed by using a satellite spatial resolution as a lattice point.
According to the technical scheme provided by the embodiment of the disclosure, the characteristics of large observation scale of the meteorological satellite and high observation precision of the foundation cloud radar are effectively combined, and the meteorological satellite and the foundation cloud radar are subjected to time-space matching and data fusion, so that the three-dimensional reconstruction of cloud in a cloud radar networking area can be rapidly and accurately realized, and the speed and precision of the three-dimensional cloud reconstruction are improved; the reconstructed three-dimensional cloud body can provide various products such as cloud top height, cloud bottom height and three-dimensional cloud fusion products on the cloud radar networking area surface, and can be widely applied to cloud visual monitoring and application.
In this disclosure, after the reconstructed three-dimensional cloud is obtained, extreme values may be further obtained for the reflectivity factor and the radial velocity of the three-dimensional cloud to obtain a combined reflectivity factor and a combined radial velocity of the cloud radar networking area, where the extreme values for the reflectivity factor and the radial velocity of the three-dimensional cloud are that the maximum value is obtained for the reflectivity factor, and the minimum value is obtained for the radial velocity.
An embodiment of the present disclosure also discloses an electronic device, fig. 4 shows a block diagram of a structure of an electronic device according to an embodiment of the present disclosure, as shown in fig. 4, the electronic device 400 includes a memory 401 and a processor 402; wherein the content of the first and second substances,
the memory 401 is used to store one or more computer instructions, which are executed by the processor 402 to implement any of the method steps described above.
Fig. 5 is a schematic structural diagram of a computer system suitable for implementing a data transmission method according to an embodiment of the present disclosure.
As shown in fig. 5, the computer system 500 includes a processing unit 501 that can execute various processes in the above-described embodiments according to a program stored in a Read Only Memory (ROM)502 or a program loaded from a storage section 508 into a Random Access Memory (RAM) 503. In the RAM503, various programs and data necessary for the operation of the computer system 500 are also stored. The processing unit 501, the ROM502, and the RAM503 are connected to each other by a bus 504. An input/output (I/O) interface 505 is also connected to bus 504.
The following components are connected to the I/O interface 505: an input portion 506 including a keyboard, a mouse, and the like; an output portion 507 including a display such as a Cathode Ray Tube (CRT), a Liquid Crystal Display (LCD), and the like, and a speaker; a storage portion 508 including a hard disk and the like; and a communication section 509 including a network interface card such as a LAN card, a modem, or the like. The communication section 509 performs communication processing via a network such as the internet. A drive 510 is also connected to the I/O interface 505 as needed. A removable medium 511 such as a magnetic disk, an optical disk, a magneto-optical disk, a semiconductor memory, or the like is mounted on the drive 510 as necessary, so that a computer program read out therefrom is mounted into the storage section 508 as necessary. The processing unit 501 may be implemented as a CPU, a GPU, a TPU, an FPGA, an NPU, or other processing units.
In particular, the above described method may be implemented as a computer software program according to an embodiment of the present invention. For example, embodiments of the present invention include a computer program product comprising a computer program tangibly embodied on a medium readable thereby, the computer program comprising program code for performing the data transmission method. In such an embodiment, the computer program may be downloaded and installed from a network through the communication section 509, and/or installed from the removable medium 511.
The flowchart and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of systems, methods and computer program products according to various embodiments of the present invention. In this regard, each block in the flowcharts or block diagrams may represent a module, a program segment, or a portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that, in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems which perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.
The units or modules described in the embodiments of the present invention may be implemented by software, or may be implemented by hardware. The units or modules described may also be provided in a processor, and the names of the units or modules do not in some cases constitute a limitation of the units or modules themselves.
As another aspect, an embodiment of the present invention further provides a computer-readable storage medium, where the computer-readable storage medium may be a computer-readable storage medium included in the apparatus in the foregoing embodiment; or it may be a separate computer readable storage medium not incorporated into the device. The computer readable storage medium stores one or more programs for use by one or more processors in performing the methods described in embodiments of the invention.
The foregoing description is only exemplary of the preferred embodiments of the invention and is illustrative of the principles of the technology employed. It will be appreciated by those skilled in the art that the scope of the invention according to the embodiments of the present invention is not limited to the specific combination of the above-mentioned features, but also encompasses other embodiments in which any combination of the above-mentioned features or their equivalents is made without departing from the inventive concept. For example, the above features and (but not limited to) the features with similar functions disclosed in the embodiments of the present invention are mutually replaced to form the technical solution.
Claims (11)
1. A method for reconstructing a three-dimensional cloud body through satellite-ground cloud fusion is characterized by comprising the following steps:
acquiring ground cloud observation data of at least three cloud radars, wherein the at least three cloud radars are networked at preset station distances;
based on the obtained foundation cloud observation data, obtaining the cloud top height, the cloud bottom height, the reflectivity factor profile and the radial speed profile of each cloud radar installation point vertical headspace;
acquiring multi-channel scanning imaging radiometer data and cloud detection product data of a stationary orbit meteorological satellite;
performing space-time matching on the foundation cloud observation data, the multi-channel scanning imaging radiometer data and the cloud detection product data to obtain satellite-ground matching data;
based on the satellite-ground matching data, cloud top height inversion of the satellite in a preset area around each cloud radar installation point is achieved, the cloud top height of the preset area around each cloud radar installation point is obtained, and further the cloud top height of a cloud radar networking area is obtained;
interpolating the cloud base heights of the cloud radars to obtain the cloud base heights of the cloud radar networking areas;
matching the cloud bottom height and the cloud top height of the cloud radar networking area to obtain a three-dimensional cloud body boundary;
and performing surface interpolation on each vertical height layer based on the cloud radar reflectivity factor profile and the radial velocity profile to fill the three-dimensional cloud body to obtain a reconstructed three-dimensional cloud body, wherein the height of each vertical height layer is determined by the vertical distance resolution of the cloud radar.
2. The method of claim 1, wherein the obtaining of the cloud top height, the cloud bottom height, the reflectivity factor profile, and the radial velocity profile of each cloud radar installation point vertical headspace based on the obtained ground based cloud observation data comprises:
acquiring cloud radar power spectrum data in the foundation cloud observation data;
obtaining a main peak of a power spectrum according to the position, the width and the amplitude of a peak body in the power spectrum data, and obtaining a reflectivity factor and a radial velocity according to the main peak of the power spectrum;
acquiring power spectrum data on each vertical height layer, and obtaining a reflectivity factor profile and a radial velocity profile in the vertical direction according to the power spectrum data on each vertical height layer;
and obtaining the cloud top height and the cloud bottom height in the vertical direction according to the reflectivity factor profile and the radial speed profile.
3. The method of claim 1, the spatio-temporal matching of the ground based cloud observation data, the multi-channel scanning imaging radiometer data, and the cloud detection product data to obtain satellite-to-ground matching data, comprising:
taking M minutes before and after satellite observation as a foundation cloud observation data time window;
determining a cloud free time fraction within the time window based on the cloud detection product data;
when the cloud-free time ratio in the time window is smaller than a first preset threshold value, searching a pixel point which is closest to a satellite according to the position information of a cloud radar installation point for matching to obtain space matching data;
and respectively carrying out normalization processing on the cloud top height and the cloud bottom height of the cloud radar in the time window, and then matching the data after the normalization processing with the space matching data to obtain the satellite-ground matching data, wherein M is a positive integer.
4. The method according to claim 3, wherein the step of performing cloud top height inversion of the satellite in a preset area around each cloud radar installation point based on the satellite-ground matching data to obtain the cloud top height of each preset area of the cloud radar installation point comprises the steps of:
normalizing the cloud top height of the cloud radar in the time window to obtain the normalized cloud top height H of the matched pixel points a ；
Obtaining thermal infrared channel data Z in multi-channel scanning imaging radiometer data matched with image element points a ；
Based on the formula
Obtaining the cloud radar in the time windowCorresponding satellite-ground cloud fusion factor of observation time
Based on the formula
Obtaining the cloud top height H of other pixels around the secondary cloud radar mounting point during observation of the satellite ij Wherein i and j are the number of rows and columns of positions of other pixel points of the satellite;
acquiring the cloud top height of the preset area of each cloud radar mounting point based on the cloud top height of the pixels around all the cloud radar mounting points during observation of the satellite;
the preset area is a circular area which takes a cloud radar installation point as a center and takes a first preset distance as a radius.
5. The method according to any one of claims 1 to 4, wherein the step of combining the cloud top heights of the preset areas of the cloud radar installation points to obtain the cloud top height of the cloud radar networking area comprises the following steps:
and when the preset areas of the cloud radar installation points are overlapped, determining the average value of the cloud top heights of the preset areas of the overlapped cloud radar installation points as the cloud top height of the overlapped area.
6. The method according to any one of claims 1 to 4, wherein, when there are multiple layers of clouds in the vertical observation direction of the cloud radar, the cloud top height, the cloud bottom height, the reflectivity factor profile and the radial velocity profile of the vertical headspace of the cloud radar installation point are determined based on the observation data of the topmost cloud with the thickness greater than a second preset threshold.
7. The method of claim 1, wherein the interpolation is performed with a satellite spatial resolution as a grid point.
8. The method of claim 1, further comprising:
and taking extreme values of the reflectivity factor and the radial speed of the three-dimensional cloud body to obtain a combined reflectivity factor and a combined radial speed of the cloud radar networking area, wherein the extreme values of the reflectivity factor and the radial speed of the three-dimensional cloud body are taken as follows, the maximum value of the reflectivity factor is taken, and the minimum value of the radial speed is taken.
9. A device for reconstructing a three-dimensional cloud body through satellite-ground cloud fusion is characterized by comprising:
the system comprises a first acquisition module, a second acquisition module and a third acquisition module, wherein the first acquisition module is configured to acquire foundation cloud observation data of at least three cloud radars, and based on the acquired foundation cloud observation data, the cloud top height, the cloud bottom height, a reflectivity factor profile and a radial speed profile of each cloud radar installation point vertical headspace are obtained, and the at least three cloud radars are networked at a preset station distance;
a second acquisition module configured to acquire multi-channel scanning imaging radiometer data and cloud detection product data of a stationary orbit meteorological satellite;
a space-time matching module configured to perform space-time matching on the foundation cloud observation data, the multi-channel scanning imaging radiometer data and the cloud detection product data to obtain satellite-ground matching data;
the inversion module is configured to realize cloud top height inversion of the satellite in a preset area around each cloud radar installation point based on the satellite-ground matching data, obtain cloud top heights of the preset areas of the cloud radar installation points, and further obtain cloud top heights of a cloud radar networking area;
the cloud base height determining module is configured to interpolate the cloud base heights of the cloud radars to obtain the cloud base heights of the cloud radar networking areas;
the reconstruction module is configured to match the cloud bottom height and the cloud top height of the cloud radar networking area to obtain a three-dimensional cloud body boundary; performing on-plane interpolation on each vertical height layer based on the cloud radar reflectivity factor profile and the radial velocity profile to fill the three-dimensional cloud body, so as to obtain a reconstructed three-dimensional cloud body, wherein the height of each vertical height layer is determined by the vertical distance resolution of the cloud radar.
10. An electronic device comprising a memory and a processor; wherein the content of the first and second substances,
the memory is configured to store one or more computer instructions, wherein the one or more computer instructions are executed by the processor to implement the method steps of any of claims 1-8.
11. A computer readable storage medium having stored thereon computer instructions which, when executed by a processor, carry out the method steps of any of claims 1-8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210457626.XA CN114820966A (en) | 2022-04-27 | 2022-04-27 | Method and device for reconstructing three-dimensional cloud through satellite-ground cloud fusion, electronic equipment and medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210457626.XA CN114820966A (en) | 2022-04-27 | 2022-04-27 | Method and device for reconstructing three-dimensional cloud through satellite-ground cloud fusion, electronic equipment and medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114820966A true CN114820966A (en) | 2022-07-29 |
Family
ID=82508805
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210457626.XA Pending CN114820966A (en) | 2022-04-27 | 2022-04-27 | Method and device for reconstructing three-dimensional cloud through satellite-ground cloud fusion, electronic equipment and medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114820966A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115238514A (en) * | 2022-08-02 | 2022-10-25 | 北京华云星地通科技有限公司 | Method and system for calculating satellite load observation simulation data |
| CN117452360A (en) * | 2023-09-26 | 2024-01-26 | 北京华云星地通科技有限公司 | Matching method and system for satellite-borne radar and ground-based radar data |
-
2022
- 2022-04-27 CN CN202210457626.XA patent/CN114820966A/en active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115238514A (en) * | 2022-08-02 | 2022-10-25 | 北京华云星地通科技有限公司 | Method and system for calculating satellite load observation simulation data |
| CN117452360A (en) * | 2023-09-26 | 2024-01-26 | 北京华云星地通科技有限公司 | Matching method and system for satellite-borne radar and ground-based radar data |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114820966A (en) | Method and device for reconstructing three-dimensional cloud through satellite-ground cloud fusion, electronic equipment and medium | |
| US6340946B1 (en) | Method for determining storm predictability | |
| Nguyen et al. | Stereographic methods for cloud base height determination using two sky imagers | |
| Darden et al. | Utilizing total lightning information to diagnose convective trends | |
| CN109725317B (en) | Sea surface bright temperature imaging simulation method based on one-dimensional synthetic aperture microwave radiometer | |
| JPH07502591A (en) | Marine and meteorological data | |
| KR102151362B1 (en) | Image decoding apparatus based on airborn using polar coordinates transformation and method of decoding image using the same | |
| CN114371182A (en) | Satellite-borne GNSS-R high-precision soil moisture estimation method based on CYGNSS data | |
| CN111404593B (en) | Method and device for processing satellite remote sensing data | |
| Henning et al. | The arecibo l-band feed array zone of avoidance survey. i. precursor observations through the inner and outer galaxy | |
| EP3239735B1 (en) | Methods and apparatus for adaptive multisensor analysis and aggregation | |
| Miller et al. | Coordinated study of coherent radar backscatter and optical airglow depletions in the central Pacific | |
| CN112213721A (en) | Millimeter wave three-dimensional imaging method for scanning outer or inner scenes of cylinder for security inspection | |
| CN109945969B (en) | Method and device for determining earth radiation balance based on meteorological satellite observation | |
| Jing et al. | Sub-pixel accuracy evaluation of FY-3D MERSI-2 geolocation based on OLI reference imagery | |
| CN111487621A (en) | Sea surface flow field inversion method based on radar image and electronic equipment | |
| Suwarsono et al. | Detecting the lava flow deposits from 2018 Anak Krakatau eruption using data fusion Landsat-8 optic and Sentinel-1 SAR | |
| Narayanan et al. | Imaging observations of upper mesospheric nightglow emissions from Tirunelveli (8.7 o N) | |
| JP5035782B2 (en) | Split beam synthetic aperture radar | |
| Yuan et al. | Exploring the influence of various factors on microwave radiation image simulation for Moon-based Earth observation | |
| CN113936009B (en) | Cloud shadow removing method, device and equipment for meteorological satellite flood monitoring | |
| CN113850514A (en) | Oil-gas pipeline inspection method and system based on constellation remote sensing technology | |
| Ma et al. | Optical observation of plasma bubbles and comparative study of multiple methods of observing the ionosphere over China | |
| Patel et al. | Geometric Properties of Water-ice Clouds as Observed from Jezero Crater in the First 600 sols with the NavCam Instrument On Board the Mars2020 Rover, Perseverance | |
| Ushio et al. | Development and observation of the phased array radar at X band |
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 |