CN110987733B - Method for judging availability of infrared detection equipment in complex environment - Google Patents
Method for judging availability of infrared detection equipment in complex environment Download PDFInfo
- Publication number
- CN110987733B CN110987733B CN201911190646.XA CN201911190646A CN110987733B CN 110987733 B CN110987733 B CN 110987733B CN 201911190646 A CN201911190646 A CN 201911190646A CN 110987733 B CN110987733 B CN 110987733B
- Authority
- CN
- China
- Prior art keywords
- infrared
- mathematical model
- infrared detection
- under
- detection equipment
- 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.)
- Active
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 88
- 238000000034 method Methods 0.000 title claims description 24
- 238000013178 mathematical model Methods 0.000 claims abstract description 35
- 239000002245 particle Substances 0.000 claims abstract description 32
- 238000002834 transmittance Methods 0.000 claims abstract description 26
- 238000011156 evaluation Methods 0.000 claims abstract description 8
- 239000000443 aerosol Substances 0.000 claims abstract description 7
- 238000010606 normalization Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 11
- 238000005315 distribution function Methods 0.000 claims description 6
- 230000003287 optical effect Effects 0.000 claims description 6
- 238000004088 simulation Methods 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 3
- 238000011056 performance test Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 4
- 230000008033 biological extinction Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 3
- 230000008447 perception Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 230000005457 Black-body radiation Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/02—Investigating particle size or size distribution
- G01N15/0205—Investigating particle size or size distribution by optical means
- G01N15/0211—Investigating a scatter or diffraction pattern
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/06—Investigating concentration of particle suspensions
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/06—Investigating concentration of particle suspensions
- G01N15/075—Investigating concentration of particle suspensions by optical means
Landscapes
- Chemical & Material Sciences (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
The invention relates to an evaluation method for availability of infrared detection equipment in a complex environment, which aims to solve the problems of poor compatibility and low evaluation accuracy in the existing performance test technology. The technical scheme adopted by the invention is as follows: firstly, sending the particle data acquired by the aerosol particle sensor under the complicated weather into a mathematical modelAs an input quantity of the mathematical model, the mathematical model outputs the atmospheric visibility L under the weather condition, then the atmospheric visibility L and the sensor data are simultaneously input into Lowtran7 simulation software, the infrared transmittance tau under the atmospheric environment can be obtained, and the infrared transmittance tau and the visibility L are input into the mathematical model togetherThe detection distance R of the infrared detection equipment at the moment can be obtained, and whether the infrared detection equipment is available or not under the current complex environment is judged according to the comparison between the calculated theoretical detection distance R and the use requirement.
Description
Technical Field
The invention relates to an availability evaluation method of infrared equipment in a complex meteorological environment, in particular to an availability evaluation method of infrared detection equipment in a complex meteorological environment.
Background
The infrared detection equipment is widely applied in various fields at present, and has more important significance in the aspects of infrared guidance and early warning. However, the usability judgment method of the infrared detection equipment in the complex meteorological environment in China still has defects at present.
The usability evaluation needs to be related to the detection accuracy, the detection distance, and the detection efficiency of the detection device. At present, the performance research of scholars at home and abroad on an infrared optical system is mainly to calculate the thermodynamic characteristics of a measured object based on the Planck black body radiation law and then calculate the detectable maximum distance. Although the method has a clear mathematical model, the method does not consider the use environment, so that the method has poor compatibility and low judgment accuracy.
Disclosure of Invention
The invention aims to provide a method for judging the availability of infrared detection equipment in a complex environment, which aims to solve the problems of poor compatibility and low judgment accuracy in the existing performance test technology.
In order to achieve the purpose, the invention adopts the technical scheme that:
a method for judging the availability of infrared detection equipment in a complex environment is characterized by comprising the following steps: firstly, sending the particle data acquired by the aerosol particle sensor under the complicated weather into a mathematical modelAs an input quantity of the mathematical model, the mathematical model outputs the atmospheric visibility L under the weather condition, then the atmospheric visibility L and the sensor data are simultaneously input into Lowtran7 simulation software, the infrared transmittance tau under the atmospheric environment can be obtained, and the infrared transmittance tau and the visibility L are input into the mathematical model togetherThe detection distance R of the infrared detection equipment at the moment can be obtained, comparison is made with the use requirement according to the calculated theoretical detection distance R, and the infrared detection equipment judges that the infrared detection equipment is under the current complex environmentIt is detected whether the device is available.
The judging method comprises the following steps:
step 1, normalization stage: according to the aerosol particle data, normalizing the different sizes of rain, snow and fog in the complex weather environment into the change condition of the atmospheric visibility L by using a mathematical model I;
step 2, simulation calculation stage of detection distance of infrared equipment: sending the infrared transmittance tau output by Lowtran7 software and the visibility L obtained by normalization by using a mathematical model I into the mathematical model II, and calculating the detection distance R of infrared detection equipment;
step 3, an infrared detection device availability evaluation stage: and comparing the detection distance R of the infrared detection equipment obtained by calculation through the mathematical model II with the required quantity of the detection distance under the current weather, and judging whether the equipment is available at the moment.
The step 1 comprises the following steps:
step 101: acquiring sensor data, namely acquiring particle concentration, particle diameter and particle size distribution functions;
step 102: and calculating the acquired data, and sending the numerical value into a mathematical model I to obtain the atmospheric visibility L at the moment.
The step 2 comprises the following steps:
step 201, simulating and calculating the infrared transmittance tau of the current environment by using simulation software;
step 202: and sending the infrared transmittance tau and the atmospheric visibility L into a mathematical model to calculate the detection distance R of the infrared detection equipment.
The step 3 comprises the following steps:
step 301: comparing the calculated detection distance R of the infrared detection equipment with the current use requirement;
step 302: if the requirement is met, the device can be used; if the requirements are not met, the number of the used devices is increased or other optical principle devices are used.
Compared with the prior art, the invention has the following advantages and effects:
1. the invention combines the infrared detection equipment with the specific atmospheric environment, thereby improving the compatibility of the method.
2. The method of the invention directly sends the data collected by the sensor into the text mathematical model and Lowtran7 software, and the accuracy of the calculated detection distance is improved.
3. The invention is based on an MIE scattering model, solves the problem that the traditional infrared equipment detection distance calculation model is too complex, and can calculate the visibility of the atmosphere and the infrared transmittance at the moment by only considering the complex weather environment as the complex attribute of aerosol suspended particles and obtaining the information of the particle size, the particle size distribution function, the particle density and the like at the moment through various particle sensors under any specific severe condition, thereby obtaining the detection distance of the infrared detection equipment under the weather state.
4. The method provides a normalization basis of a complex scene, decomposes the complex problem into simple problems, greatly simplifies a mathematical model of the infrared detection distance, provides a theoretical basis for the availability of infrared equipment in a specific environment, and provides a theoretical basis for the collocation of multiple sensing equipment in the complex scene.
Drawings
FIG. 1 is a diagram showing the arrangement of Lowtran7 software Card 1;
FIG. 2 is a diagram of the Lowtran7 software atmospheric environment setup;
FIG. 3 is a diagram of Lowtran7 software geometric path parameter setup;
FIG. 4 is a diagram of a Lowtran7 software simulation run interface;
FIG. 5 is a graph of infrared transmittance versus visibility;
fig. 6 is a graph showing the variation of infrared detection distance with transmittance.
Detailed Description
The invention relates to a method for judging the availability of infrared detection equipment in a complex environment, which specifically comprises the following steps:
step 1, normalization stage: the different degrees of rain, snow and fog in the complex weather environment are calculated according to the aerosol particle data by using a mathematical modelNormalization to the change of atmospheric visibility L;
the method specifically comprises the following steps:
step 101: acquiring sensor data, namely acquiring particle concentration, particle diameter and particle size distribution functions;
step 102: and calculating the acquired data, and sending the numerical value into a mathematical model I to obtain the atmospheric visibility L at the moment.
Step 2, simulation calculation stage of detection distance of infrared equipment: the infrared transmittance tau output by Lowtran7 software and the visibility L obtained by normalization by a mathematical model (I) are sent to the mathematical modelCalculating the detection distance R of the infrared detection equipment;
the method specifically comprises the following steps:
step 201, simulating and calculating the infrared transmittance tau of the current environment by using simulation software;
step 202: and sending the infrared transmittance tau and the atmospheric visibility L into a mathematical model to calculate the detection distance R of the infrared detection equipment.
Step 3, an infrared detection device availability evaluation stage: and comparing the detection distance R of the infrared detection equipment obtained by calculation through the mathematical model II with the required quantity of the detection distance under the current weather, and judging whether the equipment is available at the moment.
The method specifically comprises the following steps:
step 301: comparing the calculated detection distance R of the infrared detection equipment with the current use requirement;
step 302: if the requirement is met, the device can be used; if the requirements are not met, the number of the used devices is increased or other optical principle devices are used.
Example (b):
calculating visibility corresponding to different severe environments in the first step
The perception principle of the image perception equipment is similar, the outline and color information of the target to be observed are collected through an optical element, and then data format conversion and the like are carried out through a rear photoelectric conversion device, at the moment, the visibility of the environment can be used for representing the observability degree of the target to be observed, namely the detection capability influence degree of the environment on the perception equipment. The multi-dimensional environmental impact factors are normalized to visibility.
The procedure for normalization is as follows:
the light intensity I after attenuation over the transmission distance L is:
I=I0exp(-βextL) (1)
wherein I0Is the initial light intensity.
Let the atmospheric transmittance epsilon be I/I0Then, there are:
when the atmospheric transmittance ∈ is 0.02, L at this time can be considered as the atmospheric visibility in this weather condition.
Wherein beta isextAs extinction coefficients, the particle number concentration N and the particle diameter D, and the extinction efficiency factor σextFunction of the particle size distribution function f (D):
the function f (D) of the particle size distribution function can be calculated again using the following formula:
wherein DgIs the median particle diameter of the number of particles, betagAs geometric standard deviation, whengWhen 1.0, it is monodisperse.
Extinction efficiency factor sigmaextThe calculation can be made from the following equation:
wherein λ is the wavelength of light, Re (a)n+bn) Is (a)n+bn) Real part of (a)n、bnAs parameters in the MIE scattering model, calculations can be performed by the MATLAB program.
The calculation formula of the atmospheric visibility L obtained in the simultaneous steps (1) to (5) is as follows:
the accurate visibility values after normalization corresponding to the multi-dimensional weather meteorological factors of different levels can be calculated through a mathematical model (i.e. formula (6)), the complex environment type is a range value, so the complex environment type is also a range value calculated according to the formula, and for the sake of simplicity, the starting range value and the ending range value of the range are listed. And rounding the value to obtain the average visibility value of the interval.
The visibility normalization results are shown in table 1 below:
TABLE 1 visibility normalization Standard
Wherein the influence of each fog level on the detection capability of the infrared detection device is similar, and the fog levels are analyzed according to the normalized categories.
The second step is to calculate the transmittance of infrared light in a complex environment
The transmittance of infrared light in different visibility, precipitation and temperature environments can be obtained by performing simulation calculation with third-party software (LOWTRAN7), and the using process is as shown in fig. 1 to 4:
(1) setting the atmospheric mode, the geometric type of the atmospheric path, the program mode, the scattering type, the output mode and the ground reflectivity as shown in the figure 1:
(2) the environmental weather data at this time is input in the interface shown in fig. 2:
(3) the geometric path parameter data at this time is input in the interface shown in fig. 3:
(4) the upper wavelength limit and the lower wavelength limit are set at the interface shown in fig. 4:
taking a human body (with a waveband range around 9.6 micrometers) as an example, the transmittance corresponding to different visibility obtained by simulation is shown in table 2 below:
TABLE 2 different visibility corresponding transmission
The variation trend of transmittance with visibility is shown in fig. 5:
thirdly, modeling the transmission of infrared light in complex environment
For the sensing device based on the optical principle, the maximum distance (detection distance R) that an object image can distinguish can be calculated according to a mathematical model (formula (7)):
wherein L is atmospheric visibility, tau is atmospheric transmittance, G is a slope atmospheric correction coefficient, and G is 0.19. The detection distances R for the IR device with the data available are shown in Table 3 below:
TABLE 3 detection distance of infrared detection device
The relationship between the detection distance and the infrared transmittance of the infrared detection device is shown in fig. 6:
fourthly, analyzing the usability of the infrared detection equipment in the complex environment
Firstly, determining the use requirement of infrared detection equipment in the current environment, then comparing the detection distance of the infrared detection equipment with the use requirement, and if the detection distance is greater than the use requirement, directly placing the equipment in the geometric center of a detection area for use; if the requirement is not met, the distribution is carried out according to the proportion of the detection range requirement and the detection distance of the equipment.
The above-described embodiments are merely illustrative of the principles and effects of the present invention, and some embodiments may be applied, and it will be apparent to those skilled in the art that various changes and modifications may be made without departing from the inventive concept of the present invention, and these embodiments are within the scope of the present invention.
Claims (5)
1. A method for judging the availability of infrared detection equipment in a complex environment is characterized by comprising the following steps: firstly, sending the particle data acquired by the aerosol particle sensor under the complicated weather into a mathematical modelFirstly, as the input quantity of a mathematical model, the mathematical model outputs the atmospheric visibility L under the weather condition, then the atmospheric visibility L and the sensor data are simultaneously input into Lowtran7 simulation software, the infrared transmittance tau under the atmospheric environment can be obtained, and the infrared transmittance tau and the atmospheric visibility L are input into the mathematical model togetherSecondly, the theoretical detection distance R of the infrared detection equipment at the moment can be obtained, whether the infrared detection equipment is available or not under the current complex environment is judged according to the comparison between the calculated theoretical detection distance R and the use requirement,
wherein, betagFor geometric standard deviation, N is the particle number concentration, λ is the wavelength of light, N is a natural number, Re (a)n+bn) Is (a)n+bn) Real part of (a)n、bnIs a parameter in the MIE scattering model, D is the particle diameter, DgIs the number of particles and the median diameter.
2. The method for evaluating the availability of the infrared detection device under the complex environment according to claim 1, wherein: the method comprises the following steps:
step 1, normalization stage: according to the aerosol particle data, normalizing the different sizes of rain, snow and fog in the complex weather environment into the change condition of the atmospheric visibility L by using a mathematical model I;
step 2, simulation calculation stage of detection distance of infrared equipment: sending the infrared transmittance tau output by Lowtran7 software and atmospheric visibility L obtained by normalization by using a mathematical model I into the mathematical model II, and calculating the theoretical detection distance R of infrared detection equipment;
step 3, an infrared detection device availability evaluation stage: and comparing the theoretical detection distance R of the infrared detection equipment obtained by calculation through a mathematical model II with the required quantity of the detection distance under the current weather, and judging whether the equipment is available at the moment.
3. The method for evaluating the availability of the infrared detection device in the complex environment according to claim 2, wherein:
the step 1 comprises the following steps:
step 101: acquiring sensor data, namely acquiring particle concentration, particle diameter and particle size distribution functions;
step 102: and calculating the acquired data, and sending the numerical value into a mathematical model I to obtain the atmospheric visibility L at the moment.
4. The method for evaluating the availability of the infrared detection device in the complex environment according to claim 2 or 3, wherein:
the step 2 comprises the following steps:
step 201, simulating and calculating the infrared transmittance tau of the current environment by using simulation software;
step 202: and (4) sending the infrared transmittance tau and the atmospheric visibility L into a mathematical model, and calculating the theoretical detection distance R of the infrared detection equipment.
5. The method for evaluating the availability of the infrared detection device under the complex environment according to claim 4, wherein:
the step 3 comprises the following steps:
step 301: comparing the calculated theoretical detection distance R of the infrared detection equipment with the current use requirement;
step 302: if the requirement is met, the device can be used; if the requirements are not met, the number of the used devices is increased or other optical principle devices are used.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911190646.XA CN110987733B (en) | 2019-11-28 | 2019-11-28 | Method for judging availability of infrared detection equipment in complex environment |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911190646.XA CN110987733B (en) | 2019-11-28 | 2019-11-28 | Method for judging availability of infrared detection equipment in complex environment |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110987733A CN110987733A (en) | 2020-04-10 |
CN110987733B true CN110987733B (en) | 2022-02-25 |
Family
ID=70087808
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911190646.XA Active CN110987733B (en) | 2019-11-28 | 2019-11-28 | Method for judging availability of infrared detection equipment in complex environment |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110987733B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113533262B (en) * | 2021-03-24 | 2022-11-04 | 北京航空航天大学 | Atmospheric aerosol infrared scattering transmittance determination method |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102183442A (en) * | 2011-02-17 | 2011-09-14 | 中国科学院安徽光学精密机械研究所 | Visibility measuring method based on atmospheric fine particle spectrometer |
CN103942439A (en) * | 2014-04-24 | 2014-07-23 | 中国科学院遥感与数字地球研究所 | Inhalable particle concentration estimating method based on meteorological observation data |
CN107368617A (en) * | 2016-05-13 | 2017-11-21 | 南京理工大学 | Air-ground detection infrared imaging system operating distance computational methods based on Lowtran7 air softwares |
CN109060312A (en) * | 2018-08-10 | 2018-12-21 | 河南科技大学 | Visible light device operating distance test method based on controllable target |
CN109374488A (en) * | 2018-11-14 | 2019-02-22 | 国网安徽省电力有限公司电力科学研究院 | Prediction technique for photovoltaic power station haze weather visibility |
KR102043818B1 (en) * | 2019-04-10 | 2019-11-12 | 주식회사 엠에이티 | Measuring device for fog density |
-
2019
- 2019-11-28 CN CN201911190646.XA patent/CN110987733B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102183442A (en) * | 2011-02-17 | 2011-09-14 | 中国科学院安徽光学精密机械研究所 | Visibility measuring method based on atmospheric fine particle spectrometer |
CN103942439A (en) * | 2014-04-24 | 2014-07-23 | 中国科学院遥感与数字地球研究所 | Inhalable particle concentration estimating method based on meteorological observation data |
CN107368617A (en) * | 2016-05-13 | 2017-11-21 | 南京理工大学 | Air-ground detection infrared imaging system operating distance computational methods based on Lowtran7 air softwares |
CN109060312A (en) * | 2018-08-10 | 2018-12-21 | 河南科技大学 | Visible light device operating distance test method based on controllable target |
CN109374488A (en) * | 2018-11-14 | 2019-02-22 | 国网安徽省电力有限公司电力科学研究院 | Prediction technique for photovoltaic power station haze weather visibility |
KR102043818B1 (en) * | 2019-04-10 | 2019-11-12 | 주식회사 엠에이티 | Measuring device for fog density |
Non-Patent Citations (1)
Title |
---|
基于Mie理论的激光雾中传输衰减特性研究;魏海亮 等;《舰船电子工程》;20141231;第34卷(第1期);第148-151页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110987733A (en) | 2020-04-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Jiang et al. | Real‐time crack assessment using deep neural networks with wall‐climbing unmanned aerial system | |
Klose et al. | Usability study of 3D time-of-flight cameras for automatic plant phenotyping | |
CN110058258B (en) | Atmospheric boundary layer detection method based on hybrid laser radar | |
CN104634279A (en) | Vision-based automatic aviation oil mist nozzle atomization angle detection device and method | |
CN107368617B (en) | Ground-air detection infrared imaging system action distance calculation method based on Lowtran7 atmospheric software | |
CN104835124A (en) | Foam image color correction method based on spectrum information fusion and device thereof | |
Li et al. | A hybrid contextual approach to wildland fire detection using multispectral imagery | |
CN110727223A (en) | Ring-type track intelligent inspection robot based on underground coal face and application thereof | |
CN110988909A (en) | TLS-based vegetation coverage determination method for sandy land vegetation in alpine and fragile areas | |
Makky | A portable low-cost non-destructive ripeness inspection for oil palm FFB | |
CN110987733B (en) | Method for judging availability of infrared detection equipment in complex environment | |
CN111398109A (en) | Atmospheric visibility measuring method, sensor module, system and storage medium | |
CN115797760A (en) | Active and passive fusion water quality three-dimensional remote sensing inversion method and system and storage medium | |
CN114739919A (en) | Water quality detection method based on spectrum inversion analysis | |
Li et al. | Vision-based measurement of dust concentration by image transmission | |
Yang et al. | An RGB channel operation for removal of the difference of atmospheric scattering and its application on total sky cloud detection | |
JP2015032205A (en) | Image processing system and image processing method | |
CN103996266B (en) | Tailing reservoir safety monitoring device based on video analysis | |
CN116466368B (en) | Dust extinction coefficient profile estimation method based on laser radar and satellite data | |
Bhatia et al. | Sensitivity of reflectance to water vapor and aerosol optical thickness | |
CN117348574A (en) | Intelligent control system and method for paint production line | |
CN113724381A (en) | Dynamic three-dimensional scene rapid reconstruction method based on high-resolution remote sensing image | |
CN109654559B (en) | Healthy cooking system capable of identifying water vapor and oil smoke | |
Babari et al. | Computer vision for the remote sensing of atmospheric visibility | |
Deng et al. | A method for forest vegetation height modeling based on aerial digital orthophoto map and digital surface model |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |