CN109839617A - A kind of non-precipitation echo removing method and system of radar network composite product - Google Patents
A kind of non-precipitation echo removing method and system of radar network composite product Download PDFInfo
- Publication number
- CN109839617A CN109839617A CN201910188842.7A CN201910188842A CN109839617A CN 109839617 A CN109839617 A CN 109839617A CN 201910188842 A CN201910188842 A CN 201910188842A CN 109839617 A CN109839617 A CN 109839617A
- Authority
- CN
- China
- Prior art keywords
- precipitation
- data
- composite product
- network composite
- radar network
- 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
Classifications
-
- 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
- Radar Systems Or Details Thereof (AREA)
- Image Processing (AREA)
Abstract
The present invention relates to the non-precipitation echo removing methods and system of a kind of radar network composite product, wherein, the non-precipitation echo removing method of radar network composite product, comprising the following steps: obtain the non-precipitation ground observation data of the original observed data of satellite sounding and the corresponding region of original observed data；Original observed data is converted into infrared bright temperature data；Non- precipitation infrared brightness temperature distributed data is generated according to infrared bright temperature data and non-precipitation ground observation data；Obtain the Satellite Observations of the corresponding region of radar network composite product and radar network composite product；Consistency treatment is carried out to radar network composite product and Satellite Observations, obtains processing result；It is identified according to non-precipitation infrared brightness temperature distributed data and rejects the non-precipitation echo in processing result.In the inventive solutions, the non-precipitation echo of weather radar can be identified, to realize the purpose for effectively eliminating the non-precipitation echo of weather radar.
Description
Technical field
The present invention relates to meteorological data monitoring management fields more particularly to a kind of non-precipitation echo of radar network composite product to disappear
Except the non-precipitation echo of method and a kind of radar network composite product eliminates system.
Background technique
When China New Generation Weather Radar is cloudless on high or has cloud without precipitation, it can often observe that the clear sky of large area returns
Wave.This kind of echo generate the main reason for be have maintain the long period, have very strong air refraction index horizontal gas-bearing formation or
It is reflected caused by the vertical gas-bearing formation of person, scattering and aerosol, insect etc. caused by refractive index caused by atmospheric turbulance rises and falls are drawn
The scattering risen.Since weather radar data is widely used in quantitative measurement of rainfall and Data Assimilation, the presence of clear air echo is to drop
Water estimates that accuracy and Data Assimilation quality generate large effect.Therefore, in China New Generation Weather Radar software systems or precipitation
In nowcasting operation system, clear air echo identification and elimination are carried out, there is important meaning to the application effect for improving Radar Data
Justice.
Currently, clear air echo recognizer has inclining test algorithm and the clear air echo recognition methods based on fuzzy logic.
According to statistics, the recognition accuracy of both methods is all lower, 40% or so.The identification elimination of clear air echo is still considered weather
The huge challenge that radar data quality control business need to face.
Summary of the invention
The present invention is directed to solve at least one of the technical problems existing in the prior art or related technologies.
For this purpose, it is an object of the present invention to provide a kind of non-precipitation echo removing method of radar network composite product,
The non-precipitation echo of weather radar can be identified, to realize the purpose for effectively eliminating the non-precipitation echo of weather radar.
It is another object of the present invention to provide a kind of non-precipitation echoes of radar network composite product to eliminate system, can
The non-precipitation echo of weather radar is identified, to realize the purpose for effectively eliminating the non-precipitation echo of weather radar.
To achieve the above object, the technical solution of first aspect present invention provides a kind of non-precipitation of radar network composite product
Method for echo cancellation, comprising the following steps: obtain the original observed data of satellite sounding and the correspondence area of original observed data
The non-precipitation ground observation data in domain；Original observed data is converted into infrared bright temperature data；According to infrared bright temperature data and non-
Precipitation ground observation data generate non-precipitation infrared brightness temperature distributed data；Obtain radar network composite product and radar network composite product
The Satellite Observations of corresponding region；Consistency treatment is carried out to radar network composite product and Satellite Observations, obtains processing knot
Fruit；It is identified according to non-precipitation infrared brightness temperature distributed data and rejects the non-precipitation echo in processing result.
In the technical scheme, the original observed data of satellite sounding and the automatic Weather Station for corresponding to original observed data are utilized
The history ground observation data of observation combine, and statistically analyze the probability distribution of precipitation and the TBB under non-precipitation condition, and use
Infrared brightness temperature threshold method and multichannel differential technique based on high-resolution Meteorological Satellites, can effectively improve non-to weather radar
The accuracy rate that precipitation echo is identified is reduced with realizing the purpose for effectively eliminating the non-precipitation echo of weather radar because clear sky returns
Influence of the presence of wave to Calculation of precipitation accuracy and Data Assimilation quality.
In the above-mentioned technical solutions, it is preferable that original observed data is converted into infrared bright temperature data, comprising: will be original
Observation data are converted to radiation value；Infrared brightness temperature is generated according to radiation value, and generates infrared bright temperature data according to longitude and latitude sequence.
In the technical scheme, in order to establish non-precipitation echo recognition function, meteorological satellite and surface observations are utilized
The satellite infrared bright temperature TBB probability distribution of automatic website corresponding region when the method statistic analysis precipitation and clear sky that combine, with
Determine that non-precipitation echo corresponds to the bright temperature TBB threshold value of satellite infrared.
In any of the above-described technical solution, it is preferable that non-precipitation ground observation data are that the rainfall of automatic Weather Station observation is small
In or equal to 0.2mm ground observation data.
In any of the above-described technical solution, it is preferable that carried out at consistency to radar network composite product and Satellite Observations
Reason, obtains processing result, comprising: counts infrared brightness temperature value according to Satellite Observations, generates statistical result；To statistical result into
Row processing, obtains statistical attribute；Radar network composite product is handled according to statistical attribute, obtains processing result.
In the technical scheme, since the factors such as satellite and radar observation mode difference, propagation time delay will cause
There are sterically defined deviations when observing the same area for satellite and radar.Using the method for simple " point-to-point " carry out satellite and
Radar Data spatial position consistency treatment is insecure.To make up deficiency existing for " point-to-point " processing method, this method
It is eliminated using the processing method of " point is to block " to avoid normal meteorological echo.After above-mentioned pre-processing of the information, identical longitude and latitude
Radar and the corresponding lattice point coordinate of satellite data are also identical.
In any of the above-described technical solution, it is preferable that identified according to non-precipitation infrared brightness temperature distributed data and reject processing
As a result the non-precipitation echo in, comprising: processing result is counted according to non-precipitation infrared brightness temperature distributed data, and is united
Evaluation；According to statistical value compared with preset threshold, when statistical value is less than preset threshold, it is identified as non-precipitation echo.
In any of the above-described technical solution, it is preferable that the corresponding region of original observed data and pair of radar network composite product
The radius for answering region is 30km；And/or preset threshold is 10%.
In any of the above-described technical solution, it is preferable that when Satellite Observations are clear sky data, statistical value are as follows:
Non- precipitation infrared brightness temperature distributed data are as follows:
When Satellite Observations are thin cloud data, statistical value are as follows:
Distribution occasion of the non-precipitation infrared brightness temperature distributed data in different interchannels are as follows:
Wherein, IR1-IR2 < 0 is expressed as Bao Yun, and R indicates the radius of corresponding region.
The non-precipitation echo that the technical solution of second aspect of the present invention provides a kind of radar network composite product eliminates system, packet
Include: first obtains module, is arranged to be used for obtaining the original observed data of satellite sounding and the correspondence of original observed data
The non-precipitation ground observation data in region；Conversion module is arranged to be used for for original observed data being converted into infrared brightness temperature number
According to；Processing module is arranged to be used for generating non-precipitation infrared brightness according to infrared bright temperature data and non-precipitation ground observation data
Warm distributed data；Second obtains module, is arranged to be used for obtaining the correspondence area of radar network composite product and radar network composite product
The Satellite Observations in domain；Consistency treatment module is arranged to be used for carrying out radar network composite product and Satellite Observations
Consistency treatment obtains processing result；Identification module is arranged to be used for being identified simultaneously according to non-precipitation infrared brightness temperature distributed data
Reject the non-precipitation echo in processing result.
In the technical scheme, the original observed data of satellite sounding and the automatic Weather Station for corresponding to original observed data are utilized
The history ground observation data of observation combine, and statistically analyze the probability distribution of precipitation and the TBB under non-precipitation condition, and use
Infrared brightness temperature threshold method and multichannel differential technique based on high-resolution Meteorological Satellites, can effectively improve non-to weather radar
The accuracy rate that precipitation echo is identified is reduced with realizing the purpose for effectively eliminating the non-precipitation echo of weather radar because clear sky returns
Influence of the presence of wave to Calculation of precipitation accuracy and Data Assimilation quality.
In the above-mentioned technical solutions, it is preferable that conversion module includes: converting unit, is arranged to be used for original observation
Data are converted to radiation value；Infrared brightness temperature processing unit is arranged to be used for generating infrared brightness temperature according to radiation value, and according to warp
Latitude sequence generates infrared bright temperature data；Consistency treatment module includes: the first statistic unit, is arranged to be used for according to satellite
Observational data statistical infrared brightness temperature value generates statistical result；Processing unit is originally provided for handling statistical result,
Obtain statistical attribute；Consistency treatment unit is arranged to be used for handling radar network composite product according to statistical attribute, obtain
To processing result；Identification module includes: the second statistic unit, is arranged to be used for according to non-precipitation infrared brightness temperature distributed data pair
Processing result is counted, and obtains statistical value；Recognition unit is arranged to be used for according to statistical value compared with preset threshold,
When statistical value is less than preset threshold, it is identified as non-precipitation echo.
In the technical scheme, in order to establish non-precipitation echo recognition function, meteorological satellite and surface observations are utilized
The satellite infrared bright temperature TBB probability distribution of automatic website corresponding region when the method statistic analysis precipitation and clear sky that combine, with
Determine that non-precipitation echo corresponds to the bright temperature TBB threshold value of satellite infrared.Due to satellite and radar observation mode difference, propagation time delay
Etc. factors will cause satellite and radar there are sterically defined deviations when observing the same area.Using simple " point-to-point "
Method carries out satellite and Radar Data spatial position consistency treatment is insecure.Exist to make up " point-to-point " processing method
Deficiency, this method is eliminated using the processing method of " point to block " to avoid normal meteorological echo.Through above-mentioned pre-processing of the information
Afterwards, identical longitude and latitude radar and the corresponding lattice point coordinate of satellite data are also identical.
In any of the above-described technical solution, it is preferable that non-precipitation ground observation data are that the rainfall of automatic Weather Station observation is small
In or equal to 0.2mm ground observation data；And/or original observed data corresponding region and radar network composite product correspondence area
The radius in domain is 30km；And/or preset threshold is 10%；And/or Satellite Observations be clear sky data when, statistical value are as follows:
Non- precipitation infrared brightness temperature distributed data are as follows:
When Satellite Observations are thin cloud data, statistical value are as follows:
Distribution occasion of the non-precipitation infrared brightness temperature distributed data in different interchannels are as follows:
Wherein, IR1-IR2 < 0 is expressed as Bao Yun, and R indicates the radius of corresponding region.
Detailed description of the invention
Above-mentioned and/or additional aspect of the invention and advantage will become from the description of the embodiment in conjunction with the following figures
Obviously and it is readily appreciated that, in which:
Fig. 1 shows the flow diagram of non-precipitation echo removing method involved by one embodiment of the invention；
Fig. 2 shows the flow diagrams of steps involved of embodiment of the present invention S600；
Fig. 3 shows the flow diagram of steps involved of embodiment of the present invention S200；
Fig. 4 shows the flow diagram of steps involved of embodiment of the present invention S500；
Fig. 5 shows the structural block diagram that non-precipitation echo involved by one embodiment of the invention eliminates system；
Fig. 6 shows the structural block diagram of identification module involved by the embodiment of the present invention；
Fig. 7 shows the structural block diagram of conversion module involved by the embodiment of the present invention；
Fig. 8 shows the structural block diagram of consistency treatment module involved by the embodiment of the present invention；
Fig. 9 shows effect picture before and after 6:00 clear air echo Quality Control on July 28th, 2016；
(a) the intensity reflectogram without Quality Control；(b) the intensity reflectogram after existing algorithm Quality Control；
(c) radar satellite unified algorithm identifies back echo figure；(d) satellite infrared cloud atlas；
(e) national automatic Weather Station hour Pluviogram；
Figure 10 shows effect picture before and after 21:50 clear air echo Quality Control on July 27th, 2016；
(a) the intensity reflectogram without Quality Control；(b) the intensity reflectogram after existing algorithm Quality Control；
(c) radar satellite unified algorithm identifies back echo figure；(d) satellite infrared cloud atlas；
(e) national automatic Weather Station hour Pluviogram；
Figure 11 shows effect picture before and after 16:20 clear air echo Quality Control on June 14th, 2016；
(a) the intensity reflectogram without Quality Control；(b) the intensity reflectogram after existing algorithm Quality Control；
(c) radar satellite unified algorithm identifies back echo figure；(d) satellite infrared cloud atlas.
Specific embodiment
To better understand the objects, features and advantages of the present invention, with reference to the accompanying drawing and specific real
Applying mode, the present invention is further described in detail.It should be noted that in the absence of conflict, the implementation of the application
Feature in example and embodiment can be combined with each other.
In the following description, numerous specific details are set forth in order to facilitate a full understanding of the present invention, still, the present invention may be used also
To be implemented using other than the one described here other modes, therefore, protection scope of the present invention is not limited to following public affairs
The limitation for the specific embodiment opened.
The non-precipitation echo elimination side of the radar network composite product of some embodiments of the invention is described referring to Fig. 1 to Fig. 8
Method and system.
As shown in Figure 1, the non-precipitation echo removing method of the radar network composite product according to one embodiment of the invention, including
Following steps:
S100 obtains the original observed data of satellite sounding and the corresponding region (radius 30km) of original observed data
Non- precipitation ground observation data (i.e. the rainfall of automatic Weather Station observation is less than or equal to the ground observation data of 0.2mm)；
Original observed data is converted into infrared bright temperature data by S200；
S300 generates non-precipitation infrared brightness temperature distributed data according to infrared bright temperature data and non-precipitation ground observation data；
S400 obtains the moonscope number of the corresponding region (radius 30km) of radar network composite product and radar network composite product
According to；
S500 carries out consistency treatment to radar network composite product and Satellite Observations, obtains processing result；
S600 is identified according to non-precipitation infrared brightness temperature distributed data and is rejected the non-precipitation echo in processing result.
In this embodiment, it is seen using the original observed data of satellite sounding with the automatic Weather Station for corresponding to original observed data
The history ground observation data of survey combine, and statistically analyze the probability distribution of precipitation and the TBB under non-precipitation condition, and use base
In the infrared brightness temperature threshold method and multichannel differential technique of high-resolution Meteorological Satellites, can effectively improve to the non-drop of weather radar
The accuracy rate that water echo is identified is reduced with realizing the purpose for effectively eliminating the non-precipitation echo of weather radar because of clear air echo
Influence of the presence to Calculation of precipitation accuracy and Data Assimilation quality.
Specifically, non-precipitation echo removing method includes but is not limited to following technical scheme:
Embodiment 1 (as depicted in figs. 1 and 2)
Non- precipitation echo removing method the following steps are included:
S100 obtains the original observed data of satellite sounding and the corresponding region (radius 30km) of original observed data
Non- precipitation ground observation data (i.e. the rainfall of automatic Weather Station observation is less than or equal to the ground observation data of 0.2mm)；
Original observed data is converted into infrared bright temperature data by S200；
S300 generates non-precipitation infrared brightness temperature distributed data according to infrared bright temperature data and non-precipitation ground observation data；
S400 obtains the moonscope number of the corresponding region (radius 30km) of radar network composite product and radar network composite product
According to；
S500 carries out consistency treatment to radar network composite product and Satellite Observations, obtains processing result；
S601 counts processing result according to non-precipitation infrared brightness temperature distributed data, and obtains statistical value；
S602 when statistical value is less than preset threshold, is identified as non-drop according to statistical value compared with preset threshold (10%)
Otherwise water echo is precipitation echo；
Wherein, when Satellite Observations are clear sky data, statistical value are as follows:
The radius of R expression corresponding region；
Non- precipitation infrared brightness temperature distributed data are as follows:
In this embodiment, if there are echoes in radar observation region, and satellite data shows that the region is clear sky area, just
Think echo in the region be it is cloudless under the conditions of non-precipitation echo.Method particularly includes: it is assumed that radar mosaic data lattice point
R(i, j)On there are Echo Ratings, find satellite data corresponding lattice point S(i, j)；To S(i, j)The clear sky feature of data in a certain range of periphery
Attribute is counted；Statistical value is more than that given threshold is considered as R(i, j)It is considered as clear air echo.
In the present embodiment, TbbI, jFor the longitudes and latitudes lattice point such as Radar Data R(i, j)Upper corresponding satellite cloud picture infrared brightness temperature value,
PI, jFor with satellite data lattice point S(i, j)The ratio of window lattice point number is accounted for for the characteristic attribute statistical value that center radius is R range.When
When greater than a certain threshold value, it is believed that the lattice point is clear sky area.If there are radar returns for corresponding radar lattice point, just it is identified as cloudless
Under the conditions of non-precipitation echo.
Embodiment 2 (as depicted in figs. 1 and 2)
Non- precipitation echo removing method the following steps are included:
S100 obtains the original observed data of satellite sounding and the corresponding region (radius 30km) of original observed data
Non- precipitation ground observation data (i.e. the rainfall of automatic Weather Station observation is less than or equal to the ground observation data of 0.2mm)；
Original observed data is converted into infrared bright temperature data by S200；
S300 generates non-precipitation infrared brightness temperature distributed data according to infrared bright temperature data and non-precipitation ground observation data；
S400 obtains the moonscope number of the corresponding region (radius 30km) of radar network composite product and radar network composite product
According to；
S500 carries out consistency treatment to radar network composite product and Satellite Observations, obtains processing result；
S601 counts processing result according to non-precipitation infrared brightness temperature distributed data, and obtains statistical value；
S602, according to statistical value compared with preset threshold (10%), when statistical value is greater than or equal to preset threshold, identification
For non-precipitation echo；
Wherein, when Satellite Observations are thin cloud data, statistical value are as follows:
Distribution occasion of the non-precipitation infrared brightness temperature distributed data in different interchannels are as follows:
Wherein, IR1-IR2 < 0 is expressed as Bao Yun, and R indicates the radius of corresponding region.
In this embodiment, since meteorological satellite infrared channel 11.2um (IR1), 12.3um (IR2) bands of a spectrum are essentially identical,
All it is the brightness temperature of displaying target object, but is influenced by water vapor absorption that there are apparent differences.It is seen according to different infrared channels
Survey same target object obtains infrared brightness temperature value and has differences this characteristic, can carry out the identification of some varieties of clouds.Multi-pass is used herein
Road interpolation method identifies the clear air echo in transition region.Method particularly includes: it is assumed that radar mosaic lattice point R(i, j)It is upper to exist back
Wave number finds satellite data corresponding lattice point S(i, j)；To the above-mentioned non-Precipitation Clouds characteristic attribute of a certain range of data in periphery into
Row statistics；Statistical value is more than that given threshold is considered as R(i, j) be thin cloud under the conditions of non-precipitation echo.With IR1-IR2 value in text
It is counted as characteristic attribute, carries out the clear air echo identification under the conditions of thin cloud.
In the present embodiment, wherein IR1-IR2 < 0 is expressed as Bao Yun, PS (i, j)For with satellite data lattice point S(i, j)For in
Heart radius is that the characteristic attribute statistical value of R range accounts for the ratio of window lattice point number.Work as PS (i, j)When greater than given threshold, it is believed that should
Lattice point overhead is thin cloud covered areas.If corresponding Radar Data lattice point thinks that the clear sky under the conditions of thin cloud is returned there are radar return
Wave.
In any of the above-described embodiment, it is preferable that as shown in figure 3, S200, is converted into infrared brightness temperature for original observed data
Data, comprising:
Original observed data is converted to radiation value by S201；
S202 generates infrared brightness temperature according to radiation value, and generates infrared bright temperature data according to longitude and latitude sequence.
In this embodiment, in order to establish non-precipitation echo recognition function, meteorological satellite and surface observations phase are utilized
In conjunction with method statistic analysis precipitation and when clear sky automatic website corresponding region the bright temperature TBB probability distribution of satellite infrared, with true
Fixed non-precipitation echo corresponds to the bright temperature TBB threshold value of satellite infrared.
In any of the above-described embodiment, it is preferable that as shown in figure 4, S500, to radar network composite product and Satellite Observations
Consistency treatment is carried out, processing result is obtained, comprising:
S501 counts infrared brightness temperature value according to Satellite Observations, generates statistical result；
S502 handles statistical result, obtains statistical attribute；
S503 handles radar network composite product according to statistical attribute, obtains processing result.
In this embodiment, it is defended since the factors such as satellite and radar observation mode difference, propagation time delay will cause
There are sterically defined deviations when observing the same area for star and radar.Satellite and thunder are carried out using the method for simple " point-to-point "
It is insecure up to data spatial position consistency treatment.To make up deficiency existing for " point-to-point " processing method, this method is adopted
It is eliminated with the processing method of " point is to block " to avoid normal meteorological echo.After above-mentioned pre-processing of the information, identical longitude and latitude thunder
It is also identical up to lattice point coordinate corresponding with satellite data.
As shown in figure 5, eliminating system according to the non-precipitation echo of the radar network composite product of another embodiment of the present invention
1000, comprising:
First obtains module 100, is arranged to be used for obtaining the original observed data of satellite sounding and original observation number
According to corresponding region non-precipitation ground observation data；
Conversion module 200 is arranged to be used for for original observed data being converted into infrared bright temperature data；
Processing module 300 is arranged to be used for generating non-drop according to infrared bright temperature data and non-precipitation ground observation data
Water infrared brightness temperature distributed data；
Second obtains module 400, is arranged to be used for obtaining the correspondence area of radar network composite product and radar network composite product
The Satellite Observations in domain；
Consistency treatment module 500 is arranged to be used for carrying out consistency to radar network composite product and Satellite Observations
Processing, obtains processing result；
Identification module 600 is arranged to be used for being identified according to non-precipitation infrared brightness temperature distributed data and rejects processing result
In non-precipitation echo.
In this embodiment, it is seen using the original observed data of satellite sounding with the automatic Weather Station for corresponding to original observed data
The history ground observation data of survey combine, and statistically analyze the probability distribution of precipitation and the TBB under non-precipitation condition, and use base
In the infrared brightness temperature threshold method and multichannel differential technique of high-resolution Meteorological Satellites, can effectively improve to the non-drop of weather radar
The accuracy rate that water echo is identified is reduced with realizing the purpose for effectively eliminating the non-precipitation echo of weather radar because of clear air echo
Influence of the presence to Calculation of precipitation accuracy and Data Assimilation quality.
Specifically, it includes but is not limited to following technical side that the non-precipitation echo of radar network composite product, which eliminates system 1000,
Case:
Embodiment 3 (as shown in Figure 5 and Figure 6)
The non-precipitation echo of radar network composite product eliminates system 1000
First obtains module 100, is arranged to be used for obtaining the original observed data of satellite sounding and original observation number
According to corresponding region (radius 30km) non-precipitation ground observation data (i.e. automatic Weather Station observation rainfall be less than or equal to
The ground observation data of 0.2mm)；
Conversion module 200 is arranged to be used for for original observed data being converted into infrared bright temperature data；
Processing module 300 is arranged to be used for generating non-drop according to infrared bright temperature data and non-precipitation ground observation data
Water infrared brightness temperature distributed data；
Second obtains module 400, is arranged to be used for obtaining the correspondence area of radar network composite product and radar network composite product
The Satellite Observations in domain (radius 30km)；
Consistency treatment module 500 is arranged to be used for carrying out consistency to radar network composite product and Satellite Observations
Processing, obtains processing result；
Second statistic unit 601 is arranged to be used for carrying out processing result according to non-precipitation infrared brightness temperature distributed data
Statistics, and obtain statistical value；
Recognition unit 602 is arranged to be used for according to statistical value compared with preset threshold, is greater than or equal in statistical value pre-
If being identified as non-precipitation echo when threshold value (10%)；
Wherein, when Satellite Observations are clear sky data, statistical value are as follows:
The radius of R expression corresponding region；
Non- precipitation infrared brightness temperature distributed data are as follows:
In this embodiment, if there are echoes in radar observation region, and satellite data shows that the region is clear sky area, just
Think echo in the region be it is cloudless under the conditions of non-precipitation echo.Method particularly includes: it is assumed that radar mosaic data lattice point
R(i, j)On there are Echo Ratings, find satellite data corresponding lattice point S(i, j)；To S(i, j)The clear sky feature of data in a certain range of periphery
Attribute is counted；Statistical value is more than that given threshold is considered as R(i, j)It is considered as clear air echo.
In the present embodiment, TbbI, jFor the longitudes and latitudes lattice point such as Radar Data R(i, j)Upper corresponding satellite cloud picture infrared brightness temperature value,
PI, jFor with satellite data lattice point S(i, j)The ratio of window lattice point number is accounted for for the characteristic attribute statistical value that center radius is R range.When
When greater than a certain threshold value, it is believed that the lattice point is clear sky area.If there are radar returns for corresponding radar lattice point, just it is identified as cloudless
Under the conditions of non-precipitation echo.
Embodiment 4 (as shown in Figure 5 and Figure 6)
The non-precipitation echo of radar network composite product eliminates system 1000
First obtains module 100, is arranged to be used for obtaining the original observed data of satellite sounding and original observation number
According to corresponding region (radius 30km) non-precipitation ground observation data (i.e. automatic Weather Station observation rainfall be less than or equal to
The ground observation data of 0.2mm)；
Conversion module 200 is arranged to be used for for original observed data being converted into infrared bright temperature data；
Processing module 300 is arranged to be used for generating non-drop according to infrared bright temperature data and non-precipitation ground observation data
Water infrared brightness temperature distributed data；
Second obtains module 400, is arranged to be used for obtaining the correspondence area of radar network composite product and radar network composite product
The Satellite Observations in domain (radius 30km)；
Consistency treatment module 500 is arranged to be used for carrying out consistency to radar network composite product and Satellite Observations
Processing, obtains processing result；
Second statistic unit 601 is arranged to be used for carrying out processing result according to non-precipitation infrared brightness temperature distributed data
Statistics, and obtain statistical value；
Recognition unit 602 is arranged to be used for according to statistical value compared with preset threshold, is greater than or equal in statistical value pre-
If being identified as non-precipitation echo when threshold value (10%)；
Wherein, when Satellite Observations are thin cloud data, statistical value are as follows:
Distribution occasion of the non-precipitation infrared brightness temperature distributed data in different interchannels are as follows:
Wherein, IR1-IR2 < 0 is expressed as Bao Yun, and R indicates the radius of corresponding region.
In this embodiment, since meteorological satellite infrared channel 11.2um (IR1), 12.3um (IR2) bands of a spectrum are essentially identical,
All it is the brightness temperature of displaying target object, but is influenced by water vapor absorption that there are apparent differences.It is seen according to different infrared channels
Survey same target object obtains infrared brightness temperature value and has differences this characteristic, can carry out the identification of some varieties of clouds.Multi-pass is used herein
Road interpolation method identifies the clear air echo in transition region.Method particularly includes: it is assumed that radar mosaic lattice point R(i, j)It is upper to exist back
Wave number finds satellite data corresponding lattice point S(i, j)；To the above-mentioned non-Precipitation Clouds characteristic attribute of a certain range of data in periphery into
Row statistics；Statistical value is more than that given threshold is considered as R(i, j) be thin cloud under the conditions of non-precipitation echo.With IR1-IR2 value in text
It is counted as characteristic attribute, carries out the clear air echo identification under the conditions of thin cloud.
In the present embodiment, wherein IR1-IR2 < 0 is expressed as Bao Yun, PS (i, j)For with satellite data lattice point S(i, j)For in
Heart radius is that the characteristic attribute statistical value of R range accounts for the ratio of window lattice point number.Work as PS (i, j)When greater than given threshold, it is believed that should
Lattice point overhead is thin cloud covered areas.If corresponding Radar Data lattice point thinks that the clear sky under the conditions of thin cloud is returned there are radar return
Wave.
In any of the above-described embodiment, it is preferable that as shown in fig. 7, conversion module 200 includes:
Converting unit 201 is arranged to be used for original observed data being converted to radiation value；
Infrared brightness temperature processing unit 202 is arranged to be used for generating infrared brightness temperature according to radiation value, and arranges according to longitude and latitude
Sequence generates infrared bright temperature data.
In this embodiment, in order to establish non-precipitation echo recognition function, meteorological satellite and surface observations phase are utilized
In conjunction with method statistic analysis precipitation and when clear sky automatic website corresponding region the bright temperature TBB probability distribution of satellite infrared, with true
Fixed non-precipitation echo corresponds to the bright temperature TBB threshold value of satellite infrared.
In any of the above-described embodiment, it is preferable that as shown in figure 8, consistency treatment module 500 includes:
First statistic unit 501 is arranged to be used for counting infrared brightness temperature value according to Satellite Observations, generates statistics knot
Fruit；
Processing unit 502 is originally provided for handling statistical result, obtains statistical attribute；
Consistency treatment unit 503 is arranged to be used for handling radar network composite product according to statistical attribute, obtain
Processing result.
In this embodiment, it is defended since the factors such as satellite and radar observation mode difference, propagation time delay will cause
There are sterically defined deviations when observing the same area for star and radar.Satellite and thunder are carried out using the method for simple " point-to-point "
It is insecure up to data spatial position consistency treatment.To make up deficiency existing for " point-to-point " processing method, this method is adopted
It is eliminated with the processing method of " point is to block " to avoid normal meteorological echo.After above-mentioned pre-processing of the information, identical longitude and latitude thunder
It is also identical up to lattice point coordinate corresponding with satellite data.
Technical effect
For assess the recognition methods of radar satellite joint effect, to number of echoes in survey region -30 days on the 1st July in 2016
According to testing.Wherein, including 30143 precipitation echoes and 294220 non-precipitation echoes.Statistical analysis discovery, radar satellite
It is 88.5% that joint recognition methods, which obtains total sample recognition correct rate, and non-precipitation echo recognition correct rate is 78.8%, and precipitation returns
Wave recognition correct rate is 98.2%.
A example 1
In morning on July 28th, 2016, by secondary high control, the Yangtze river basin central and east is large stretch of clear sky area.Fig. 9 (a) is as it can be seen that river
There is large stretch of continuous echo in the ground such as Soviet Union, Anhui, Hubei, Henan, most of intensity is between 10-30dBz.Compare the front and back 6:00
Time picture mosaic product discovery when multiple, the echo that 5:00 starts the region extend from coastal upcountry fast development, echo face more than 6 points
Product is maximum, gradually reduces later, meets the non-features of precipitation echo of clear sky.Comparison diagram 9 (a), (b) are as can be seen that existing algorithm has
Effect identifies background return, but unobvious to the recognition effect of the non-precipitation echo of clear sky.Comparison diagram 9 (b), (c) can be seen that
Radar satellite joint recognizer can identify and reject on a large scale the non-precipitation echo of clear sky, while precipitation echo completely being protected
It stays.It is tested radar satellite unified algorithm recognition effect as it can be seen that by knowing using automatic Weather Station precipitation data (Fig. 9 (e))
Not Wei the non-precipitation echo of clear sky and do not occur precipitation in the region rejected, Central Shandong Province, northeast of Jiangsu Province, Chongqing, Hubei southwest
The precipitation echo in portion is effectively retained.
A example 2
Evening on July 27th, 2016 forms shear line between secondary height and continental high almospheric pressure, and West of Hubei Province has Convective Cloud to be formed,
Such as Figure 10 (d).The original radar mosaic of that night 21:50, it can be seen that radial electricity occur in Yongchuan, Ankang, Xi'an, Luoyang radar
Magnetic disturbance, echo strength are up to 50dBZ, such as Figure 10 (a).After existing non-precipitation echo elimination algorithm control, part
Diametral interference echo is identified, but on the whole non-precipitation echo without be improved significantly, such as Figure 10 (b).Join through radar satellite
After closing recognizer Quality Control, non-precipitation echo is also significantly eliminated.Herewith, (the peace of the diametral interference echo under the conditions of clear sky
Also effectively being identified and being rejected Kang Leida), such as Figure 10 (c).Using 1 hour precipitation data of national automatic Weather Station to radar
Satellite joint recognizer effect test as it can be seen that be identified as non-precipitation echo and reject region in do not occur precipitation,
The precipitation echo in the area such as Shaanxi, Henan, Hubei can remain (Figure 10 (e)).
A example 3
In afternoon on 2016 month June 14, the Lezhou Peninsula, Northern Hainan Island are clear sky area by secondary high control, such as Figure 11 (d).When
Day 16:20 there is large stretch of weak echo region in the Lezhou Peninsula, Northern Hainan Island before quality control, there are precipitation time for Beibu Bay
Wave area, such as Figure 11 (a).After existing Quality Control algorithm, large stretch of weak echo region is identified as non-precipitation echo and rejects, red collimation mark
It is considered normal meteorological echo at knowledge to be retained, such as Figure 11 (b).After radar satellite combines recognizer Quality Control, residual
Weak echo is further identified and is rejected, while marine echo (Figure 11 b marking frame) is identified as non-precipitation echo and is rejected
(Figure 11 c).
Since region can not be tested with 1 hour precipitation of automatic Weather Station, to examine radar satellite to combine recognizer
Correctness has access to the satellite remote dada in the region and the Haikou closed on and Port of Fangcheng Weather Radar Product data respectively.Satellite
Cloud atlas data shows the region without precipitation cloud, such as Fig. 4 (d).3 elevation angle PPI data such as 0.5,1.5,2.4 ° of Haikou radar are equal
It has been shown that, the region is without precipitation echo information；Port of Fangcheng radar shows that there are echo letters in the region in 0.5 ° of elevation angle PPI data
Breath, but raise (1.5 °, 2.4 °) echo disappearance behind the elevation angle.By thinking that the radar return in the region is sea echo, and
Non- precipitation echo；It is correct that radar satellite, which combines recognizer,.
In the present invention, term " first ", " second ", " third " are only used for the purpose of description, and should not be understood as indicating
Or imply relative importance；Term " multiple " then refers to two or more, unless otherwise restricted clearly.Term " installation ",
The terms such as " connected ", " connection ", " fixation " shall be understood in a broad sense, for example, " connection " may be a fixed connection, being also possible to can
Dismantling connection, or be integrally connected；" connected " can be directly connected, can also be indirectly connected through an intermediary.For this
For the those of ordinary skill in field, the specific meanings of the above terms in the present invention can be understood according to specific conditions.
In description of the invention, it is to be understood that the instructions such as term " on ", "lower", "left", "right", "front", "rear"
Orientation or positional relationship is to be based on the orientation or positional relationship shown in the drawings, and is merely for convenience of the description present invention and simplification is retouched
It states, rather than the device or unit of indication or suggestion meaning must have specific direction, be constructed and operated in a specific orientation,
It is thus impossible to be interpreted as limitation of the present invention.
In the description of this specification, the description of term " one embodiment ", " some embodiments ", " specific embodiment " etc.
Mean that particular features, structures, materials, or characteristics described in conjunction with this embodiment or example are contained at least one reality of the invention
It applies in example or example.In the present specification, schematic expression of the above terms are not necessarily referring to identical embodiment or reality
Example.Moreover, description particular features, structures, materials, or characteristics can in any one or more of the embodiments or examples with
Suitable mode combines.
The foregoing is only a preferred embodiment of the present invention, is not intended to restrict the invention, for the skill of this field
For art personnel, the invention may be variously modified and varied.All within the spirits and principles of the present invention, made any to repair
Change, equivalent replacement, improvement etc., should all be included in the protection scope of the present invention.
Claims (10)
1. a kind of non-precipitation echo removing method of radar network composite product, which comprises the following steps:
It sees on the non-precipitation ground obtained in the original observed data of satellite sounding and the corresponding region of the original observed data
Measured data；
The original observed data is converted into infrared bright temperature data；
Non- precipitation infrared brightness temperature distributed data is generated according to the infrared bright temperature data and the non-precipitation ground observation data；
Obtain the Satellite Observations of the corresponding region of the radar network composite product and the radar network composite product；
Consistency treatment is carried out to the radar network composite product and the Satellite Observations, obtains processing result；
It is identified according to the non-precipitation infrared brightness temperature distributed data and rejects the non-precipitation echo in the processing result.
2. the non-precipitation echo removing method of radar network composite product according to claim 1, which is characterized in that by the original
Observation data conversion begin into infrared bright temperature data, comprising:
The original observed data is converted into radiation value；
Infrared brightness temperature is generated according to the radiation value, and generates the infrared bright temperature data according to longitude and latitude sequence.
3. the non-precipitation echo removing method of radar network composite product according to claim 1, which is characterized in that the non-drop
Water ground observation data are that the rainfall of automatic Weather Station observation is less than or equal to the ground observation data of 0.2mm.
4. the non-precipitation echo removing method of radar network composite product according to any one of claim 1 to 3, feature exist
In carrying out consistency treatment to the radar network composite product and the Satellite Observations, obtain processing result, comprising:
Infrared brightness temperature value is counted according to the Satellite Observations, generates statistical result；
The statistical result is handled, statistical attribute is obtained；
The radar network composite product is handled according to the statistical attribute, obtains the processing result.
5. the non-precipitation echo removing method of radar network composite product according to claim 4, which is characterized in that according to described
Non- precipitation infrared brightness temperature distributed data identifies and rejects the non-precipitation echo in the processing result, comprising:
The processing result is counted according to the non-precipitation infrared brightness temperature distributed data, and obtains statistical value；
According to the statistical value compared with preset threshold, when the statistical value is less than the preset threshold, it is identified as described non-
Precipitation echo.
6. the non-precipitation echo removing method of radar network composite product according to claim 5, it is characterised in that: described original
The radius for observing the corresponding region of data and the corresponding region of the radar network composite product is 30km；And/or the preset threshold
It is 10%.
7. the non-precipitation echo removing method of radar network composite product according to claim 5, it is characterised in that:
When the Satellite Observations are clear sky data, the statistical value are as follows:
The non-precipitation infrared brightness temperature distributed data are as follows:
When the Satellite Observations are thin cloud data, the statistical value are as follows:
Distribution occasion of the non-precipitation infrared brightness temperature distributed data in different interchannels are as follows:
Wherein, IR1-IR2 < 0 is expressed as Bao Yun, and R indicates the radius of corresponding region.
8. a kind of non-precipitation echo of radar network composite product eliminates system characterized by comprising
First obtains module, is arranged to be used for obtaining the original observed data and the original observed data of satellite sounding
The non-precipitation ground observation data of corresponding region；
Conversion module is arranged to be used for the original observed data being converted into infrared bright temperature data；
Processing module is arranged to be used for generating non-drop according to the infrared bright temperature data and the non-precipitation ground observation data
Water infrared brightness temperature distributed data；
Second obtains module, is arranged to be used for obtaining the correspondence area of the radar network composite product and the radar network composite product
The Satellite Observations in domain；
Consistency treatment module is arranged to be used for carrying out consistency to the radar network composite product and the Satellite Observations
Processing, obtains processing result；
Identification module is arranged to be used for being identified according to the non-precipitation infrared brightness temperature distributed data and rejects the processing result
In non-precipitation echo.
9. the non-precipitation echo of radar network composite product according to claim 8 eliminates system, which is characterized in that the conversion
Module includes:
Converting unit is arranged to be used for the original observed data being converted to radiation value；
Infrared brightness temperature processing unit is arranged to be used for generating infrared brightness temperature according to the radiation value, and sorts according to longitude and latitude
Generate the infrared bright temperature data；
The consistency treatment module includes:
First statistic unit is arranged to be used for counting infrared brightness temperature value according to the Satellite Observations, generates statistical result；
Processing unit is originally provided for handling the statistical result, obtains statistical attribute；
Consistency treatment unit is arranged to be used for handling the radar network composite product according to the statistical attribute, obtain
To the processing result；
The identification module includes:
Second statistic unit is arranged to be used for carrying out the processing result according to the non-precipitation infrared brightness temperature distributed data
Statistics, and obtain statistical value；
Recognition unit is arranged to be used for according to the statistical value compared with preset threshold, is less than in the statistical value described pre-
If when threshold value, being identified as the non-precipitation echo.
10. the non-precipitation echo of radar network composite product according to claim 9 eliminates system, it is characterised in that:
The non-precipitation ground observation data are that the rainfall of automatic Weather Station observation is less than or equal to the ground observation data of 0.2mm；
And/or
The radius of the corresponding region of the original observed data and the corresponding region of the radar network composite product is 30km；And/or
The preset threshold is 10%；And/or
When the Satellite Observations are clear sky data, the statistical value are as follows:
The non-precipitation infrared brightness temperature distributed data are as follows:
When the Satellite Observations are thin cloud data, the statistical value are as follows:
Distribution occasion of the non-precipitation infrared brightness temperature distributed data in different interchannels are as follows:
Wherein, IR1-IR2 < 0 is expressed as Bao Yun, and R indicates the radius of corresponding region.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910188842.7A CN109839617A (en) | 2019-03-13 | 2019-03-13 | A kind of non-precipitation echo removing method and system of radar network composite product |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910188842.7A CN109839617A (en) | 2019-03-13 | 2019-03-13 | A kind of non-precipitation echo removing method and system of radar network composite product |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109839617A true CN109839617A (en) | 2019-06-04 |
Family
ID=66885656
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910188842.7A Pending CN109839617A (en) | 2019-03-13 | 2019-03-13 | A kind of non-precipitation echo removing method and system of radar network composite product |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109839617A (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111947707A (en) * | 2020-07-03 | 2020-11-17 | 中国气象局兰州干旱气象研究所 | Arid and semi-arid region ground surface water circulation full-component monitoring and identifying method |
| CN113311416A (en) * | 2021-05-10 | 2021-08-27 | 中国科学院地理科学与资源研究所 | Mountain region small watershed radar quantitative precipitation estimation technology |
| CN113534090A (en) * | 2021-07-14 | 2021-10-22 | 中国科学院大气物理研究所 | Inversion method and device for liquid water content in cloud |
| CN117237677A (en) * | 2023-11-15 | 2023-12-15 | 南京信息工程大学 | Precipitation prediction correction method for overall similarity of strong precipitation space based on deep learning |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2011196916A (en) * | 2010-03-23 | 2011-10-06 | Mitsubishi Electric Corp | Measuring vehicle, and road feature measuring system |
| CN108761408A (en) * | 2018-06-24 | 2018-11-06 | 中国气象局上海台风研究所 | A method of the assessment non-precipitation echo recognizer effect of ground weather radar |
-
2019
- 2019-03-13 CN CN201910188842.7A patent/CN109839617A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2011196916A (en) * | 2010-03-23 | 2011-10-06 | Mitsubishi Electric Corp | Measuring vehicle, and road feature measuring system |
| CN108761408A (en) * | 2018-06-24 | 2018-11-06 | 中国气象局上海台风研究所 | A method of the assessment non-precipitation echo recognizer effect of ground weather radar |
Non-Patent Citations (3)
| Title |
|---|
| 李爽 等: "应用FY_3A_MERSI数据反演土壤水分的研究", 《现代农业科技》 * |
| 程昌玉 等: "基于气象卫星资料的天气雷达非降水回波消除方法", 《气象与减灾研究》 * |
| 肖笑: "静止气象卫星资料的强对流云团的识别与预报研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111947707A (en) * | 2020-07-03 | 2020-11-17 | 中国气象局兰州干旱气象研究所 | Arid and semi-arid region ground surface water circulation full-component monitoring and identifying method |
| CN113311416A (en) * | 2021-05-10 | 2021-08-27 | 中国科学院地理科学与资源研究所 | Mountain region small watershed radar quantitative precipitation estimation technology |
| CN113311416B (en) * | 2021-05-10 | 2024-05-28 | 中国科学院地理科学与资源研究所 | Mountain small-basin radar quantitative precipitation estimation method |
| CN113534090A (en) * | 2021-07-14 | 2021-10-22 | 中国科学院大气物理研究所 | Inversion method and device for liquid water content in cloud |
| CN113534090B (en) * | 2021-07-14 | 2024-01-30 | 中国科学院大气物理研究所 | Inversion method and device for liquid water content in cloud |
| CN117237677A (en) * | 2023-11-15 | 2023-12-15 | 南京信息工程大学 | Precipitation prediction correction method for overall similarity of strong precipitation space based on deep learning |
| CN117237677B (en) * | 2023-11-15 | 2024-03-26 | 南京信息工程大学 | Precipitation prediction correction method for overall similarity of strong precipitation space based on deep learning |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109839617A (en) | A kind of non-precipitation echo removing method and system of radar network composite product | |
| Cummins et al. | 6.1. The US National Lightning Detection Network: post-upgrade status | |
| Anagnostou | A convective/stratiform precipitation classification algorithm for volume scanning weather radar observations | |
| Robbins et al. | Changepoints in the North Atlantic tropical cyclone record | |
| KR101258668B1 (en) | Korea local radar processing system | |
| Collis et al. | Statistics of storm updraft velocities from TWP-ICE including verification with profiling measurements | |
| CN109814175B (en) | Strong convection monitoring method based on satellite and application thereof | |
| Kober et al. | Tracking and nowcasting of convective cells using remote sensing data from radar and satellite | |
| Sieglaff et al. | A satellite-based convective cloud object tracking and multipurpose data fusion tool with application to developing convection | |
| CN110261857B (en) | Spatial interpolation method for weather radar | |
| Seto et al. | Applicability of the iterative backward retrieval method for the GPM dual-frequency precipitation radar | |
| CN113642475B (en) | Atlantic hurricane strength estimation method based on convolutional neural network model | |
| CN116449331B (en) | Dust particle number concentration estimation method based on W-band radar and meteorological satellite | |
| Forsythe et al. | Ionospheric horizontal correlation distances: Estimation, analysis, and implications for ionospheric data assimilation | |
| Peterson et al. | Thunderstorm cloud-type classification from space-based lightning imagers | |
| CN115062527A (en) | Geostationary satellite sea temperature inversion method and system based on deep learning | |
| Berne et al. | Influence of the vertical profile of reflectivity on radar-estimated rain rates at short time steps | |
| Ionescu et al. | DeePS at: A deep learning model for prediction of satellite images for nowcasting purposes | |
| Chen et al. | Proactive quality control: Observing system experiments using the NCEP global forecast system | |
| CN109283600A (en) | A kind of visibility automatic Observation and artificial observation comparing appraisal procedure and system | |
| Kida et al. | Improvement of rain/no-rain classification methods for microwave radiometer observations over the ocean using a 37 GHz emission signature | |
| CN109541565A (en) | A kind of radar echo intensity homogeneity detection method and system | |
| CN115292667A (en) | Space probability analysis method and system based on rainfall forecast and remote correlation corresponding relation | |
| Zhang | Inter-Comparison of Space-and Ground-Based Observations of Lightning | |
| Bandholnopparat et al. | Optical properties of intracloud and cloud-to-ground discharges derived from JEM-GLIMS lightning observations |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190604 |