CN113343435A - Long-wave emission radiation calculation method suitable for AGRI instrument on FY4A satellite - Google Patents

Long-wave emission radiation calculation method suitable for AGRI instrument on FY4A satellite Download PDF

Info

Publication number
CN113343435A
CN113343435A CN202110550963.9A CN202110550963A CN113343435A CN 113343435 A CN113343435 A CN 113343435A CN 202110550963 A CN202110550963 A CN 202110550963A CN 113343435 A CN113343435 A CN 113343435A
Authority
CN
China
Prior art keywords
olr
channel
satellite
long
agri
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
Application number
CN202110550963.9A
Other languages
Chinese (zh)
Inventor
张婉春
吴晓
白文广
张鹏
孙凌
刘健
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
National Satellite Meteorological Center
Original Assignee
National Satellite Meteorological Center
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by National Satellite Meteorological Center filed Critical National Satellite Meteorological Center
Priority to CN202110550963.9A priority Critical patent/CN113343435A/en
Publication of CN113343435A publication Critical patent/CN113343435A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Algebra (AREA)
  • Probability & Statistics with Applications (AREA)
  • Operations Research (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Image Processing (AREA)

Abstract

The invention discloses a calculation method for long-wave radiation emitted by an AGRI instrument on an FY4A satellite, which comprises the following steps: reading IMAGER observation data; FY4-RTM mode simulation; carrying out statistical regression analysis to derive a regression coefficient; the method is suitable for the technical field of meteorological satellite OLR calculation, adopts multi-channel data to carry out OLR inversion, can better reflect the influence of near-ground gas and water vapor on the radiation quantity, and has higher inversion precision; meanwhile, the channel information used by the method has extremely high correlation with the OLR, and the atmospheric top-emitted long-wave radiant flux-OLR is calculated by the satellite channel radiance by establishing the regression relationship between satellite channel observation and OLR simulation, so that the method overcomes the limitation of the existing method, obtains higher precision, has a larger application range and is strong in popularization.

Description

Long-wave emission radiation calculation method suitable for AGRI instrument on FY4A satellite
Technical Field
The invention belongs to the technical field of OLR calculation of meteorological satellites, and particularly relates to a calculation method for long-wave radiation emitted by an AGRI instrument on an FY4A satellite.
Background
The atmospheric top emits long-wave radiation, called OLR for short, is the thermal radiation flux density radiated to the outside air by the earth atmospheric system, is one of important parameters of earth atmospheric radiation energy balance, the OLR space-time distribution observed by a satellite can reflect the thermal condition of an underlying surface, and the size of the OLR is closely related to the emission of the earth surface in clear areas; in cloudy areas, the OLR is related to the cloud top temperature, and high cloud amount is usually corresponding to a lower OLR value, so that OLR data obtained by satellite inversion can be used for researching tropical convection and precipitation estimation, ENSO, ITCZ, monsoon and the like;
however, in the existing satellite inversion algorithm, different band combinations are considered to establish the inversion algorithm, and because each channel reflects the characteristics of different properties of the earth surface or the atmosphere, only a proper infrared channel combination is selected to obtain a reasonable calculation result, so that the method has great limitation, and is difficult to reflect the influence of near-ground gas and water vapor on the radiation quantity, and the inversion accuracy is not high.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provides a calculation method for long-wave radiation emitted by an AGRI instrument on an FY4A satellite.
In order to achieve the purpose, the invention adopts the following technical scheme:
a method for calculating the emission long-wave radiation of an AGRI instrument on an adaptive FY4A satellite comprises the following steps:
reading IMAGER observation data;
FY4-RTM mode simulation;
carrying out statistical regression analysis to derive a regression coefficient;
and calculating the long-wave radiation OLR emitted by the atmospheric top.
Preferably, the IMAGER observation data includes AGRI channel radiance and satellite zenith angle.
Preferably, the AGRI channels include channel 9, channel 10, channel 12, channel 14, and a water vapor channel.
Preferably, the FY4-RTM mode simulation includes: based on the atmospheric profile, the AGRI channel spectral response function and the satellite zenith angle, calculating the long-wave radiance emitted by the simulated atmospheric dome by using FY4-RTM software, wherein the calculation formula is as follows:
Figure BDA0003075354990000021
Figure BDA0003075354990000022
Figure BDA0003075354990000023
wherein the content of the first and second substances,
Figure BDA0003075354990000024
for simulating the atmospheric top-emitted long-wave radiance, epsilonvIs the surface spectral radiance, B is the Planck function, c1，c2Is the radiation constant, v is the wavenumber, τvIs monochromatic transmittance, theta is the local zenith angle,
Figure BDA0003075354990000025
is an azimuth angle, ztIs the top height of the atmosphere, z' is any height in the atmosphere, TsIs the ground temperature, n is the atmosphere stratification, and the formula (3) is a series form of the formula (1).
Preferably, the statistical regression analysis for deriving the regression coefficient includes: based on FY-4AGRIOLR inversion mode:
Figure BDA0003075354990000031
based on ECMWF2010 reanalysis data, selecting a plurality of atmosphere contour lines of simulated atmospheric dome emissivity and AGRI channel emissivity, inputting the selected values into a multiple linear regression software, and counting a regression coefficient a0、ai、bi、ci、di、ei。
Preferably, 2521 is selected as the optimal choice of the simulated atmospheric dome emission long-wave radiance of a plurality of atmospheric profiles based on the ECMWF2010 reanalysis data.
Preferably, the statistical regression analysis derives regression coefficients for satellite zenith angles at 0 °, 10 °, 20 °, 30 °, 40 °, 50 °, 60 °, 70 °, and 80 °.
Preferably, the calculating of the atmospheric top emission long-wave radiation OLR includes the following steps:
based on the statistical regression analysis result, finding out two adjacent angle value ranges to which the satellite zenith angles belong;
searching and acquiring a corresponding regression coefficient according to the value range;
calculating OLR of two adjacent angles;
and performing interpolation processing to obtain the OLR of the satellite zenith angle.
Preferably, the interpolation processing adopts a pixel-by-pixel processing mode.
In summary, due to the adoption of the technical scheme, the invention has the beneficial effects that:
according to the invention, the OLR inversion is carried out by adopting multi-channel data, so that the influence of near-ground gas and steam on the radiation quantity can be better reflected, and the inversion precision is higher;
the channel information used by the method has extremely high correlation with the OLR, and the atmospheric top-emitted long-wave radiant flux-OLR is calculated by the satellite channel radiance by establishing the regression relationship between the satellite channel observation and the OLR, so that the method overcomes the limitation of the existing method, obtains higher precision, has a larger application range and is strong in popularization.
Drawings
FIG. 1 is a flow chart of a method of calculating the long-wave emission of an AGRI apparatus adapted to an FY4A satellite according to the present invention.
Detailed Description
The following further describes a specific embodiment of the method for calculating the long-wave radiation emitted by the AGRI instrument on the FY4A satellite according to the present invention with reference to fig. 1. The method for calculating the long-wave radiation emitted by the AGRI apparatus adapted to the FY4A satellite is not limited to the following description of the embodiment.
Example 1:
the embodiment provides a specific structure of a calculation method for long-wave radiation emitted by an AGRI instrument adapted to an FY4A satellite, as shown in fig. 1, which includes the following steps:
reading IMAGER observation data;
FY4-RTM mode simulation;
carrying out statistical regression analysis to derive a regression coefficient;
and calculating the long-wave radiation OLR emitted by the atmospheric top.
In particular, the IMAGER observation data comprises AGRI channel radiance and satellite zenith angle.
Specifically, the AGRI channels include channel 9, channel 10, channel 12, channel 14, and a moisture channel.
Specifically, the simulation of FY4-RTM mode comprises the following steps: based on the atmospheric profile, the AGRI channel spectral response function and the satellite zenith angle, calculating the long-wave radiance emitted by the simulated atmospheric dome by using FY4-RTM software, wherein the calculation formula is as follows:
Figure BDA0003075354990000051
Figure BDA0003075354990000052
Figure BDA0003075354990000053
wherein the content of the first and second substances,
Figure BDA0003075354990000054
for simulating the atmospheric top-emitted long-wave radiance, epsilonvIs the surface spectral radiance, B is the Planck function, c1，c2Is the radiation constant, v is the wavenumber, τvThe single color transmittance is shown, and theta is the local dayThe top angle is a vertical angle,
Figure BDA0003075354990000055
is an azimuth angle, ztIs the top height of the atmosphere, z' is any height in the atmosphere, TsIs the ground temperature, n is the atmosphere stratification, and the formula (3) is a series form of the formula (1).
Specifically, statistical regression analysis, deriving regression coefficients, includes: based on FY-4AGRIOLR inversion mode:
Figure BDA0003075354990000056
based on ECMWF2010 reanalysis data, selecting a plurality of atmosphere contour lines of simulated atmospheric dome emissivity and AGRI channel emissivity, inputting the selected values into a multiple linear regression software, and counting a regression coefficient a0、ai、bi、ci、di、ei。
Further, based on the ECMWF2010 reanalysis data, 2521 is selected as the optimal selection of the simulated atmosphere ejection long-wave radiance of a plurality of atmosphere profiles.
Further, statistical regression analysis is carried out to derive the regression coefficients of the satellite zenith angles at 0 °, 10 °, 20 °, 30 °, 40 °, 50 °, 60 °, 70 °, and 80 °.
Further, calculating the long-wave radiation OLR emitted by the atmospheric top comprises the following steps:
based on the statistical regression analysis result, finding out two adjacent angle value ranges to which the satellite zenith angles belong;
searching and acquiring a corresponding regression coefficient according to the value range;
calculating OLR of two adjacent angles;
and performing interpolation processing to obtain the OLR of the satellite zenith angle.
Furthermore, the interpolation process adopts a pixel-by-pixel processing mode.
Example 2:
the embodiment provides a verification method of an emitting long-wave radiation calculation method suitable for an AGRI instrument on an FY4A satellite, generally, the precision of product quality is illustrated by comparing with high-precision satellite data similar to those in foreign countries, the OLR product is obtained by processing grade 1 data of FY4A/AGRI by using the technology, the OLR product is compared with the instantaneous OLR of an Aqua satellite CERES instrument, RMSE is 8-11W/m2, under the condition of uniform underlying surface, the OLR product is 3-6W/m2, FY4A/AGRI is compared with HIRSCDR daily average OLR, and RMSE is 6-8W/m 2; FY4A in comparison to the HIRSCDR mean OLR, RMSE between 4 and 6W/m 2; FY4A compared with HIRSCDR average OLR, RMSE between 4-5W/m 2; the comparison of FY4A with HIRSCDR monthly average OLR, RMSE between 3 and 4W/m2, and these verification results prove that the method is suitable for calculating the FY4A/AGRI data to emit long-wave radiation products.
The working principle is as follows: as shown in fig. 1, the physical basis for inversion mode establishment is: the long-wave radiation emitted from the atmosphere roof is determined by the temperature of the emitting surface, the air temperature and the content of the absorbed gas, and the CO in the atmosphere2、O3、CH4The content is fixed, the emission energy mainly changes with the surface temperature and the water vapor, an FY-4AGRI channel 12(10.8 μm) provides surface temperature information in an atmospheric window area, channels 9(6.3 μm) and 10(7.1 μm) detect the water vapor of 500 and 700, the channel information has extremely high correlation with OLR, and the atmospheric top emission long wave radiant flux-OLR is calculated by establishing a regression relation between satellite channel observation and OLR and the satellite channel radiance;
firstly, reading in an atmospheric profile, radiation rates of AGRI channels 9, 10, 12 and 14 and a satellite zenith angle, and inquiring spectral response functions of the AGRI channels 9, 10, 12 and 14;
secondly, calculating the long-wave radiance emitted by the simulated atmospheric dome by using FY4-RTM software based on the atmospheric profile, the AGRI channel spectral response function and the satellite zenith angle;
then, based on the FY-4AGRIOLR inversion mode (4): inputting the simulated atmospheric top emission long-wave radiance and the AGRI channel radiance of 2512 atmospheric profiles into a multiple linear regression software, and counting a regression coefficient a0、ai、bi、ci、di、ei；
Finally, based on the statistical regression analysis result, finding out two adjacent angle value ranges to which the satellite zenith angle belongs, searching and obtaining corresponding regression coefficients according to the value ranges, calculating the OLR of the two adjacent angles, and performing interpolation processing to obtain the OLR of the satellite zenith angle;
in the OLR inversion processing process, input data are only the radiation rate of a satellite channel and the satellite zenith angle, an OLR inversion mode provides regression coefficients of the satellite zenith angle at 0 degrees, 10 degrees, 20 degrees, 30 degrees, 40 degrees, 50 degrees, 60 degrees, 70 degrees and 80 degrees, two adjacent angle value ranges (such as 0 degrees and 10 degrees) to which the satellite zenith angle theta belongs are found out firstly during calculation, corresponding regression coefficients are searched according to the value ranges, OLRs of two angles are calculated, and the OLRs of the theta angle are worked out through interpolation processing; because the OLR inversion processing relates to the OLR interpolation processing problem of zenith angles of 2 adjacent satellites, the pixel-by-pixel processing is most safe and reasonable;
according to the invention, the OLR inversion is carried out by adopting multi-channel data, so that the influence of near-ground gas and steam on the radiation quantity can be better reflected, and the inversion precision is higher;
the channel information used by the method has extremely high correlation with the OLR, and the atmospheric top-emitted long-wave radiant flux-OLR is calculated by the satellite channel radiance by establishing the regression relationship between the satellite channel observation and the OLR, so that the method overcomes the limitation of the existing method, obtains higher precision, has a larger application range and is strong in popularization.
The foregoing is a more detailed description of the invention in connection with specific preferred embodiments and it is not intended that the invention be limited to these specific details. For those skilled in the art to which the invention pertains, several simple deductions or substitutions can be made without departing from the spirit of the invention, and all shall be considered as belonging to the protection scope of the invention.

Claims (9)

1. A method for calculating the emission long-wave radiation of an AGRI instrument on an FY4A satellite is characterized by comprising the following steps:
reading IMAGER observation data;
FY4-RTM mode simulation;
carrying out statistical regression analysis to derive a regression coefficient;
and calculating the long-wave radiation OLR emitted by the atmospheric top.
2. The method of claim 1, wherein the method comprises the steps of: the IMAGER observation data comprise AGRI channel radiance and satellite zenith angles.
3. The method of claim 2, wherein the method comprises the steps of: the AGRI channel includes a channel 9, a channel 10, a channel 12, a channel 14, and a water vapor channel.
4. The method of claim 3, wherein the FY4-RTM model simulation comprises: based on the atmospheric profile, the AGRI channel spectral response function and the satellite zenith angle, calculating the long-wave radiance emitted by the simulated atmospheric dome by using FY4-RTM software, wherein the calculation formula is as follows:
Figure FDA0003075354980000011
Figure FDA0003075354980000012
Figure FDA0003075354980000021
wherein the content of the first and second substances,
Figure FDA0003075354980000022
for simulating the atmospheric top-emitted long-wave radiance, epsilonvIs the surface spectral radiance, B is the Planck function, c1，c2Is the radiation constant, v is the wavenumber, τvIs monochromatic transmittance, theta is the local zenith angle,
Figure FDA0003075354980000023
is an azimuth angle, ztIs the top height of the atmosphere, z' is any height in the atmosphere, TsIs the ground temperature, n is the atmosphere stratification, and the formula (3) is a series form of the formula (1).
5. The method of claim 4, wherein the statistical regression analysis to derive the regression coefficients comprises: based on FY-4AGRIOLR inversion mode:
Figure FDA0003075354980000024
based on ECMWF2010 reanalysis data, selecting a plurality of atmosphere contour lines of simulated atmospheric dome emissivity and AGRI channel emissivity, inputting the selected values into a multiple linear regression software, and counting a regression coefficient a0、ai、bi、ci、di、ei。
6. The method of claim 5, wherein the method comprises the steps of: based on ECMWF2010 reanalysis data, 2521 optimal choices of the simulated atmospheric dome long-wave radiance of a plurality of atmospheric profiles are selected.
7. The method of claim 6, wherein the method comprises the steps of: the statistical regression analysis derives the regression coefficients of the satellite zenith angles of 0 degrees, 10 degrees, 20 degrees, 30 degrees, 40 degrees, 50 degrees, 60 degrees, 70 degrees and 80 degrees.
8. The method of claim 7, wherein the method comprises the steps of: the method for calculating the long-wave radiation OLR emitted by the atmospheric top comprises the following steps:
based on the statistical regression analysis result, finding out two adjacent angle value ranges to which the satellite zenith angles belong;
searching and acquiring a corresponding regression coefficient according to the value range;
calculating OLR of two adjacent angles;
and performing interpolation processing to obtain the OLR of the satellite zenith angle.
9. The method of claim 8, wherein the method comprises the steps of: the interpolation processing adopts a pixel-by-pixel processing mode.
CN202110550963.9A 2021-05-20 2021-05-20 Long-wave emission radiation calculation method suitable for AGRI instrument on FY4A satellite Pending CN113343435A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110550963.9A CN113343435A (en) 2021-05-20 2021-05-20 Long-wave emission radiation calculation method suitable for AGRI instrument on FY4A satellite

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110550963.9A CN113343435A (en) 2021-05-20 2021-05-20 Long-wave emission radiation calculation method suitable for AGRI instrument on FY4A satellite

Publications (1)

Publication Number Publication Date
CN113343435A true CN113343435A (en) 2021-09-03

Family

ID=77470165

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110550963.9A Pending CN113343435A (en) 2021-05-20 2021-05-20 Long-wave emission radiation calculation method suitable for AGRI instrument on FY4A satellite

Country Status (1)

Country Link
CN (1) CN113343435A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114325129A (en) * 2021-12-18 2022-04-12 国家卫星气象中心(国家空间天气监测预警中心) Emission long wave radiation value attenuation correction method suitable for FY3BVIRR satellite

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20140115777A (en) * 2013-03-22 2014-10-01 이화여자대학교 산학협력단 Method for estimating longwave climate feedback using sea surface temperature
CN107727244A (en) * 2017-11-23 2018-02-23 中国科学院地理科学与资源研究所 A kind of contactless earth's surface temperature-indicating instrument and method
US20190018918A1 (en) * 2017-07-11 2019-01-17 The Trustees Of Princeton University System and method for performing accurate hydrologic determination using disparate weather data sources

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20140115777A (en) * 2013-03-22 2014-10-01 이화여자대학교 산학협력단 Method for estimating longwave climate feedback using sea surface temperature
US20190018918A1 (en) * 2017-07-11 2019-01-17 The Trustees Of Princeton University System and method for performing accurate hydrologic determination using disparate weather data sources
CN107727244A (en) * 2017-11-23 2018-02-23 中国科学院地理科学与资源研究所 A kind of contactless earth's surface temperature-indicating instrument and method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
吴晓等: "基于FY-3/IRAS利用非线性模式反演OLR", 应用气象学报, vol. 28, no. 2, 31 March 2017 (2017-03-31), pages 189 - 199 *
方红亮: "基于地球静止气象卫星的地表参数遥感研究进展", 遥感学报, vol. 25, no. 1, 31 January 2021 (2021-01-31), pages 109 - 125 *
闵敏;吴晓;: "从FY-4A卫星遥感数据和GFS资料估算全天空状况下的地表长波辐射通量", 气象, vol. 46, no. 03, 21 March 2020 (2020-03-21), pages 336 - 345 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114325129A (en) * 2021-12-18 2022-04-12 国家卫星气象中心(国家空间天气监测预警中心) Emission long wave radiation value attenuation correction method suitable for FY3BVIRR satellite

Similar Documents

Publication Publication Date Title
CN109581372A (en) A kind of Remote Sensing Monitoring of Ecological Environment method
CN109492563A (en) A kind of tree species classification method based on unmanned plane Hyperspectral imaging and LiDAR point cloud
CN108020322B (en) Airborne thermal infrared hyperspectral remote sensing quantitative detection method for coal field fire area
Yu et al. Kriging interpolation method and its application in retrieval of MODIS aerosol optical depth
CN109297605B (en) Surface temperature inversion method based on mid-infrared and thermal infrared data
CN108985959A (en) Wheat powdery mildew remote sensing monitoring method based on surface temperature inversion technology
CN107219497A (en) Cross-radiometric calibration method and system for grand sight angular measurement sensor
CN108827878B (en) Passive microwave remote sensing detection method for surface air pressure
CN111323129A (en) Earth surface temperature inversion method based on broadband thermal infrared image
CN109946235A (en) The multi layer cloud inversion method of wind and cloud 4A meteorological satellite Multichannel Scan Imagery Radiometer
CN112749494B (en) Method for acquiring dynamic accumulated snow depth
CN110388986B (en) Land surface temperature inversion method based on TASI data
CN111737913A (en) MWHTS clear sky observation bright temperature selection method based on cloud water content inversion
CN111104888A (en) Automatic generation technology of cloud detection algorithm supported by AVIRIS high-resolution data
CN102901563B (en) Method and device for determining land surface emissivity of narrow band and broad band simultaneously
Maddy et al. Using Metop-A AVHRR clear-sky measurements to cloud-clear Metop-A IASI column radiances
CN113343435A (en) Long-wave emission radiation calculation method suitable for AGRI instrument on FY4A satellite
CN113076645B (en) Cosmic ray neutron instrument space characteristic simulation method combining remote sensing data
Filioglou et al. Profiling water vapor mixing ratios in Finland by means of a Raman lidar, a satellite and a model
CN117523406A (en) Clamping correction method for ocean water color and temperature scanner and computer readable medium
de Lange et al. Methane profiles from GOSAT thermal infrared spectra
CN116465502A (en) Atmospheric correction method for ground station radiation temperature measurement
CN105259145B (en) A kind of while remote sensing islands and reefs underwater topography and the method for atural object
CN113740263B (en) Aerosol optical thickness inversion method and atmospheric particulate remote sensing inversion method
CN109738372A (en) A kind of Airborne Hyperspectral rock mine detection multi-data processing method

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