CN114460110B - Servo system error compensation method - Google Patents
Servo system error compensation method Download PDFInfo
- Publication number
- CN114460110B CN114460110B CN202210220140.4A CN202210220140A CN114460110B CN 114460110 B CN114460110 B CN 114460110B CN 202210220140 A CN202210220140 A CN 202210220140A CN 114460110 B CN114460110 B CN 114460110B
- Authority
- CN
- China
- Prior art keywords
- detector
- gamma
- fitting
- inclination angle
- rotation angle
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/03—Investigating materials by wave or particle radiation by transmission
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/10—Different kinds of radiation or particles
- G01N2223/101—Different kinds of radiation or particles electromagnetic radiation
- G01N2223/1016—X-ray
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/60—Specific applications or type of materials
- G01N2223/646—Specific applications or type of materials flaws, defects
-
- 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
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/02—Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
The invention discloses a servo system error compensation method, which belongs to the technical field of calibration of a motion mechanism, wherein X-rays are adopted to image a marker, the position of an imaged reference marker in an image is kept unchanged to serve as a constraint condition, the horizontal position of an objective table, the height of the objective table, the inclination angle and the rotation angle of a detector are counted, a five-axis motion relation is established by adopting a polynomial fitting method, the horizontal position of the objective table follows the change of the height of the objective table, the inclination angle and the rotation angle of the detector according to the motion relation, the fact that the center of a field of view is always kept unchanged under different amplification ratios and different visual angles is realized, and the region of interest is imaged; the invention can establish the absolute position relation of the point object image in space, can not only meet the requirement of keeping the center of the field of view unchanged all the time under different amplification ratios and different visual angles, but also be applied to planar CT imaging, and can realize the three-dimensional imaging of the measured object without an additional planar CT module.
Description
Technical Field
The invention relates to the technical field of calibration of a motion mechanism, in particular to a servo system error compensation method.
Background
The X-ray inspection apparatus is mainly directed to integrated circuit package reliability inspection, inspecting a circuit and its package to inspect the presence of defects in order to determine the cause of the defects. The Chinese patent application No. CN201910592325.6 discloses an X-ray detection device based on a five-axis motion platform, which can acquire a plurality of two-dimensional images or three-dimensional models of a region of interest from different view angles or projections. The servo system may allow relative movement of the detector within the sphere and linear movement of the sample with the stage in three dimensions.
Two-dimensional imaging of a region of interest is effective, but often provides insufficient information due to shadowing or the like, requiring imaging of the region of interest from different perspectives. However, in the actual operation process, due to the error sources such as the non-orthogonal axis error, the non-concentric center of the spherical motion of the detector and the source center of the radiation source, the indication error, the detector installation error, and the like, which are caused by the installation error of the radiation source, when the stage is moved in height and the detector is rotated for detection, the region of interest deviates from the center of the field of view, and the region of interest needs to be repositioned, especially when the imaging with high magnification is performed, the magnification or the viewing angle is greatly adjusted, the region of interest deviates from the field of view completely, and the repositioning is difficult.
To solve this problem, there are two general solutions for the existing X-ray detection apparatus. One is to adjust the height of the objective table and the inclination and rotation angles of the detector at a small angle, then adjust the horizontal position of the objective table to ensure that the region of interest is always positioned in the center of the field of view, and repeatedly adjust for several times to realize the imaging of the region of interest under the expected amplification ratio and inclination angle; the method is complicated in adjustment, poor in user experience, and difficult to meet the automatic detection requirement in batch processing. Another approach is to rely on a high precision mechanical system for relatively high precision movement between the X-ray source, the region of interest, and the detector. This need for high precision mechanical equipment makes the X-ray detection device expensive. For this purpose, a servo error compensation method is proposed.
Disclosure of Invention
The technical problems to be solved by the invention are as follows: how to solve the system error of the compensation servo system, the method realizes that the center of the visual field is always kept unchanged under different amplification ratios and different visual angles, images the region of interest, and provides a servo system error compensation method.
The invention solves the technical problems through the following technical proposal, and the invention comprises the following steps:
s1: acquiring an image in real time using X-rays to identify fiducial markers on the stage surface having a known geometry;
s2: moving the object stage height z, the detector inclination angle theta and the detector rotation angle gamma, and keeping the position of the characteristic of the reference mark in the image unchanged through the horizontal positions x and y of the object stage to obtain a sequence { x, y, z, theta, gamma };
s3: extracting a sequence { x, y, z, θ, γ=0 } when γ is 0, and respectively establishing the relationship between x, y, z and θ by adopting a data fitting method:
wherein p is x00 And p y00 Is constant, f x And f y A corresponding fitting function;
s4: performing circle fitting on { x, y } in the sequences of the same z, the same θ and different γ to obtain the sequences { z, θ, x } c ,y c R }, where x c And y c The values of x and y corresponding to the fitted circle center are obtained, and r is the radius of the fitted circle;
s5: establishing r and gamma by adopting data fitting method 0 Relationship with z and θ:
wherein, gamma 0 To fit the initial angle of the circle, y d Y is relative to y when γ=0 c Offset of f r And f yd A corresponding fitting function;
s6: compensating for servo errors using fitted functions, dynamically correcting p if the stage moves horizontally x00 And p y00 :
If the stage height z, the detector inclination angle θ, and the detector rotation angle γ change, then:
wherein x is T 、y T The following position is the following position after the object stage height z, the detector inclination angle theta and the detector rotation angle gamma are changed.
Still further, in said step S1, said fiducial markers of known geometry are metal balls.
Still further, in the step S2, the reference mark is characterized by a geometric center, a centroid, and a minimum circumscribing center.
Further, in the steps S3 and S5, the data fitting method is a polynomial fitting method, a neural network fitting method, or the like.
Compared with the prior art, the invention has the following advantages: the servo system error compensation method is suitable for an X-ray detection device based on a five-axis motion platform, takes metal microspheres as markers, adopts X-rays to image the markers, takes the position of the imaged circle center in the visual field as a constraint condition, counts the horizontal position, the height and the inclination angle of a detector of an objective table, and adopts a polynomial fitting method to establish a five-axis motion relation, and according to the motion relation, the horizontal position of the objective table changes along with the height, the inclination angle and the rotation angle of the objective table, so that the central of the visual field is always kept unchanged under different amplification ratios and different visual angles, and the imaging of a concerned region is realized; the method can establish the absolute position relation of the point object image in space, can meet the requirement of keeping the center of the field of view unchanged all the time under different amplification ratios and different visual angles, can be applied to planar CT imaging, and can realize three-dimensional imaging of an object to be measured without an additional planar CT module; the servo system error is compensated by adopting a software method, no additional high-precision mechanical equipment and detection equipment are needed, the cost performance is high, and the servo system error compensation method is worthy of popularization and use.
Drawings
FIG. 1 is a flow chart of data acquisition and fitting of a servo system error compensation method in a second embodiment of the present invention;
FIG. 2 is a flow chart of error compensation of a servo system error compensation method according to a second embodiment of the invention;
FIG. 3a shows f in a second embodiment of the invention x (z, θ) data fitting results schematic;
FIG. 3b is a schematic diagram of the data fitting residuals in FIG. 3 a;
FIG. 4A is a diagram of f in a second embodiment of the invention y (z, θ) data fitting results schematic;
FIG. 4A is a schematic diagram of the data fit residuals in FIG. 4A;
FIG. 4B is f in a second embodiment of the invention r (z, θ) data fitting results schematic;
FIG. 4B is a schematic diagram of the data fit residuals in FIG. 4B;
FIG. 4C is f in a second embodiment of the invention yd (z, θ) data fitting results schematic;
fig. 4C is a schematic diagram of the data fit residuals in fig. 4C.
Detailed Description
The following describes in detail the examples of the present invention, which are implemented on the premise of the technical solution of the present invention, and detailed embodiments and specific operation procedures are given, but the scope of protection of the present invention is not limited to the following examples.
Example 1
The embodiment provides a technical scheme: a method for compensating for errors in a servo system, comprising the steps of:
s1: identifying a reference mark with a known geometric shape on the surface of the object stage by utilizing X-rays, and acquiring an image in real time, wherein the image contains the reference mark;
s2: moving the object stage height z, the detector inclination angle theta and the detector rotation angle gamma, and adjusting the horizontal positions x and y of the object stage to ensure that the position of the characteristic of the reference mark in the image is unchanged, so as to obtain a sequence { x, y, z, theta, gamma };
s3: extracting a sequence { x, y, z, θ, γ=0 } when γ is 0, and respectively establishing the relationship between x, y, z and θ by adopting a data fitting method:
wherein p is x00 And p y00 Is constant, f x And f y A corresponding fitting function;
s4: performing circle fitting on { x, y } in the sequences of the same z, the same θ and different γ to obtain the sequences { z, θ, x } c ,y c R }, where x c And y c The values of x and y corresponding to the fitted circle center are obtained, and r is the radius of the fitted circle;
s5: establishing r and gamma by adopting data fitting method 0 Relationship with z and θ:
wherein, gamma 0 To fit the initial angle of the circle, y d Y is relative to y when γ=0 c Offset of f r And f yd A corresponding fitting function;
s6: compensating for servo errors using fitted functions, dynamically correcting p if the stage moves horizontally x00 And p y00 :
If the stage height z, the detector inclination angle θ, and the detector rotation angle γ change, then
Wherein x is T 、y T The position is the follow-up position after the object stage height z, the detector inclination angle theta and the detector rotation angle gamma are changed.
In this embodiment, in the steps S3 and S5, the data fitting method is a polynomial fitting method, a neural network fitting method, or the like.
In this embodiment, in said step S1, said fiducial markers of known geometry are metal balls.
In this embodiment, in the step S2, the reference mark is characterized by a minimum circumscribing circle center.
In this embodiment, the stage is movable in the horizontal direction (x-axis and y-axis, two axes being perpendicular to each other) and in the vertical direction (z-axis) under the control of the servo system.
Example two
Data acquisition and processing are carried out according to the method shown in FIG. 1, the range of the change of the stage height z is 20-50 mm, the range of the detector inclination angle theta is 0-45 degrees, the range of the detector rotation angle gamma is-180 degrees, and the sampling data are processed by adopting a polynomial fitting method and a circle fitting method to obtain f x 、f y 、f r And f yd The function, the fitting conditions of which are shown in fig. 3a, 4B and 4C, and the fitting residuals are shown in fig. 3B, 4A, 4B and 4C, the residuals are less than 50 μm, which illustrates the effectiveness of the method.
In summary, the method for compensating the error of the servo system in the embodiment is suitable for an X-ray detection device based on a five-axis motion platform, uses metal microspheres as markers, uses X-rays to image the markers, uses the position of the imaged circle center in the field of view as a constraint condition, counts the horizontal position of an objective table, the height of the objective table, the inclination angle and the rotation angle of a detector, adopts a polynomial fitting method to establish a five-axis motion relationship, and realizes that the horizontal position of the objective table changes along with the height of the objective table, the inclination angle and the rotation angle of the detector according to the motion relationship, so that the center of the field of view is always kept unchanged under different amplification ratios and different visual angles, and images a region of interest; the method can establish the absolute position relation of the point object image in space, can meet the requirement of keeping the center of the field of view unchanged all the time under different amplification ratios and different visual angles, can be applied to planar CT imaging, and can realize three-dimensional imaging of an object to be measured without an additional planar CT module; the servo system error is compensated by adopting a software method, no additional high-precision mechanical equipment and detection equipment are needed, the cost performance is high, and the servo system error compensation method is worthy of popularization and use.
While embodiments of the present invention have been shown and described above, it will be understood that the above embodiments are illustrative and not to be construed as limiting the invention, and that variations, modifications, alternatives and variations may be made to the above embodiments by one of ordinary skill in the art within the scope of the invention.
Claims (3)
1. A method for compensating for errors in a servo system, comprising the steps of:
s1: acquiring an image in real time using X-rays to identify fiducial markers on the stage surface having a known geometry;
s2: moving the object stage height z, the detector inclination angle theta and the detector rotation angle gamma, and adjusting the horizontal positions x and y of the object stage to ensure that the position of the characteristic of the reference mark in the image is unchanged, so as to obtain a sequence { x, y, z, theta, gamma };
s3: extracting a sequence { x, y, z, θ, γ=0 } when γ is 0, and respectively establishing the relationship between x, y, z and θ by adopting a data fitting mode:
wherein p is x00 And p y00 Is constant, f x And f y A corresponding fitting function;
s4: performing circle fitting on { x, y } in the sequences of the same z, the same θ and different γ to obtain the sequences { z, θ, x } c ,y c R }, where x c And y c The values of x and y corresponding to the fitted circle center are obtained, and r is the radius of the fitted circle;
s5: adopting a data fitting mode to establish r and gamma 0 Relationship with z and θ:
wherein, gamma 0 To fit the initial angle of the circle, y d Y is relative to y when γ=0 c Offset of f r Anda corresponding fitting function;
s6: compensating for servo errors using fitted functions, dynamically correcting p if the stage moves horizontally x00 And p y00 :
If the stage height z, the detector inclination angle θ, and the detector rotation angle γ change, then:
wherein x is T 、y T The position is a follow-up position after the object stage height z, the detector inclination angle theta and the detector rotation angle gamma are changed;
in the step S2, the reference mark is characterized by a geometric center, a centroid, and a minimum circumscribing circle center.
2. A method of compensating for servo errors as recited in claim 1, wherein: in the step S1, the fiducial markers having a known geometry are metal balls.
3. A method of compensating for servo errors as recited in claim 1, wherein: in the steps S3 and S5, the data fitting method is a polynomial fitting method and a neural network fitting method.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210220140.4A CN114460110B (en) | 2022-03-08 | 2022-03-08 | Servo system error compensation method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210220140.4A CN114460110B (en) | 2022-03-08 | 2022-03-08 | Servo system error compensation method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114460110A CN114460110A (en) | 2022-05-10 |
CN114460110B true CN114460110B (en) | 2023-06-06 |
Family
ID=81416582
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210220140.4A Active CN114460110B (en) | 2022-03-08 | 2022-03-08 | Servo system error compensation method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114460110B (en) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2645701A1 (en) * | 2012-03-29 | 2013-10-02 | Axis AB | Method for calibrating a camera |
CN108122203A (en) * | 2016-11-29 | 2018-06-05 | 上海东软医疗科技有限公司 | A kind of bearing calibration of geometric parameter, device, equipment and system |
Family Cites Families (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6493646B1 (en) * | 2000-02-16 | 2002-12-10 | Ge Medical Systems Global Technology Company, Llc | High order primary decay correction for CT imaging system detectors |
JP4041637B2 (en) * | 2000-03-31 | 2008-01-30 | 東芝Itコントロールシステム株式会社 | X-ray fluoroscopy system |
JP3802869B2 (en) * | 2002-12-02 | 2006-07-26 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Beam hardening post-processing method and X-ray CT apparatus |
US6848827B2 (en) * | 2003-05-13 | 2005-02-01 | General Electric Company | Method and apparatus for calibrating detector spectral response |
DE102005033187A1 (en) * | 2005-07-13 | 2007-01-25 | Carl Zeiss Industrielle Messtechnik Gmbh | Measuring arrangement e.g. computed tomography measuring arrangement, calibrating method, involves determining geometry parameter of geometrical model, which describes geometry of arrangement, by evaluation of radiographic image |
US7581884B1 (en) * | 2006-02-07 | 2009-09-01 | Barnes Gary T | Mobile radiography system and grid alignment process |
US7780672B2 (en) * | 2006-02-27 | 2010-08-24 | Biomet Manufacturing Corp. | Femoral adjustment device and associated method |
TWI574125B (en) * | 2006-09-01 | 2017-03-11 | 尼康股份有限公司 | Mobile body driving method and moving body driving system, pattern forming method and apparatus, exposure method and apparatus, and component manufacturing method |
EP2326250B1 (en) * | 2008-08-13 | 2017-01-11 | Koninklijke Philips N.V. | Calibration method for ring artifact correction in non-ideal isocentric 3d rotational x-ray scanner systems using a calibration phantom based rotation center finding algorithm |
JP2012042340A (en) * | 2010-08-19 | 2012-03-01 | Shimadzu Corp | X-ray ct equipment |
DE102013205406A1 (en) * | 2013-03-27 | 2014-10-16 | Siemens Aktiengesellschaft | X-ray imaging system for X-ray imaging at high image frequencies of an examination subject by direct measurement of the interference pattern |
GB2520711B (en) * | 2013-11-28 | 2018-06-20 | Nikon Metrology Nv | Calibration apparatus and method for computed tomography |
CN103784160B (en) * | 2014-03-06 | 2015-12-23 | 北京锐视康科技发展有限公司 | A kind of correcting unit of cone-beam CT system geometric position and bearing calibration thereof |
CN104764756B (en) * | 2015-03-30 | 2017-05-31 | 清华大学 | The scaling method of cone-beam CT imaging system |
CN105769233A (en) * | 2016-02-29 | 2016-07-20 | 江苏美伦影像系统有限公司 | Geometric correction method |
CN110220926B (en) * | 2019-07-03 | 2020-03-27 | 中国电子科技集团公司第三十八研究所 | X-ray detection device based on five-axis motion platform |
CN111678472B (en) * | 2020-06-09 | 2022-02-15 | 无锡身为度信息技术有限公司 | Error identification method for rotary table of four-axis coordinate measuring machine |
JP2022010983A (en) * | 2020-06-29 | 2022-01-17 | 株式会社ミツトヨ | Method for calibrating x-ray measuring device |
CN113902645B (en) * | 2021-10-27 | 2022-12-16 | 中国电子科技集团公司第三十八研究所 | Reverse RD positioning model-based RPC correction parameter acquisition method for satellite-borne SAR image |
-
2022
- 2022-03-08 CN CN202210220140.4A patent/CN114460110B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2645701A1 (en) * | 2012-03-29 | 2013-10-02 | Axis AB | Method for calibrating a camera |
CN108122203A (en) * | 2016-11-29 | 2018-06-05 | 上海东软医疗科技有限公司 | A kind of bearing calibration of geometric parameter, device, equipment and system |
Also Published As
Publication number | Publication date |
---|---|
CN114460110A (en) | 2022-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP0457843B1 (en) | Method and apparatus for measuring registration between layers of a semiconductor wafer | |
US5825913A (en) | System for finding the orientation of a wafer | |
US6373917B1 (en) | Z-axis elimination in an X-ray laminography system using image magnification for Z plane adjustment | |
CN109029299B (en) | Dual-camera measuring device and method for butt joint corner of cabin pin hole | |
CN105849537B (en) | Calibration device and computed tomography method | |
WO2018218611A1 (en) | Geometric parameter determination method for cone beam computed tomography system | |
TW201831858A (en) | Camera and specimen alignment to facilitate large area imaging in microscopy | |
CN108168541B (en) | Improved sub-pixel star point centroid positioning method | |
CN114695225A (en) | Wafer pre-alignment device and wafer pre-alignment method | |
CN116211328A (en) | System and method for calibrating an imaging system | |
EP2068113A1 (en) | Calibrating method of image measuring instrument | |
JP4413746B2 (en) | Charged particle beam equipment | |
CN111915685B (en) | Zoom camera calibration method | |
CN209032406U (en) | A kind of cone-beam CT system geometric calibration device | |
KR101545186B1 (en) | method of correction of defect location using predetermined wafer image targets | |
US6532310B1 (en) | Removing noise caused by artifacts from a digital image signal | |
CN111376254B (en) | Plane ranging method and system and plane adjusting method and system for mechanical arm | |
CN113781579B (en) | Geometric calibration method for panoramic infrared camera | |
CN114460110B (en) | Servo system error compensation method | |
CN109938841B (en) | Surgical instrument navigation system based on multi-view camera coordinate fusion | |
JP2021043046A (en) | Calibration method for x-ray measuring device | |
CN113793313B (en) | High-precision tool setting method and device for machining full-surface micro-pit structure of thin-wall spherical shell type micro-component | |
CN114279303B (en) | Device and method for detecting verticality of double-sided micro-cylindrical lens array | |
CN114544676B (en) | Data acquisition method for high-precision error compensation of servo system | |
CN114111578B (en) | Automatic pose determining method for large-caliber element |
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 |