CN110988393A - 基于超声波风速风向仪的无人机风速风向测量及校正算法 - Google Patents
基于超声波风速风向仪的无人机风速风向测量及校正算法 Download PDFInfo
- Publication number
- CN110988393A CN110988393A CN201911278177.7A CN201911278177A CN110988393A CN 110988393 A CN110988393 A CN 110988393A CN 201911278177 A CN201911278177 A CN 201911278177A CN 110988393 A CN110988393 A CN 110988393A
- Authority
- CN
- China
- Prior art keywords
- unmanned aerial
- aerial vehicle
- wind
- wind speed
- 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.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
- G01P5/24—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting acoustical wave
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P13/00—Indicating or recording presence, absence, or direction, of movement
- G01P13/02—Indicating direction only, e.g. by weather vane
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P21/00—Testing or calibrating of apparatus or devices covered by the preceding groups
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P21/00—Testing or calibrating of apparatus or devices covered by the preceding groups
- G01P21/02—Testing or calibrating of apparatus or devices covered by the preceding groups of speedometers
- G01P21/025—Testing or calibrating of apparatus or devices covered by the preceding groups of speedometers for measuring speed of fluids; for measuring speed of bodies relative to fluids
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Aviation & Aerospace Engineering (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开了一种基于超声波风速风向仪的无人机风速风向测量及校正算法,包括两种状态的风速风向测量及校正算法,一种是无人机悬停状态下的风速风向测量及校正算法,另一种是无人机前进状态下的风速风向测量及校正算法,两种风速风向测量及校正算法都是利用四元数与姿态阵之间的关系来得出无人机的三个欧拉角,进而得到无人机的倾角,通过倾角的大小来对超声波风速风向仪测得的风速风向结果进行校正。与现有技术相比,本发明的有益效果是:本发明的测量及其校正算法克服了因无人机前进或者有风吹来时机身平面倾斜导致的测量结果误差,具有测量精度高、使用寿命长的优点。
Description
技术领域
本发明涉及一种测风校正算法,尤其涉及一种基于超声波风速风向仪的无人机风速风向测量及校正算法。
背景技术
目前,无人飞行器已在野外施工、勘探、交通运输、旅游、救援,特别是江河湖泊上工作等领域中得到广泛应用。随着无人机应用的深入,采用多旋翼无人机进行气象监测也正成为一个研究热点。
但是在现有的技术中,使用多旋翼无人机测风主要是通过直接在机体上部加装塔状超声波风速风向仪来实现的。现有的技术不仅会因为改变重心、增大风阻而导致无人机飞行稳定性的降低,而且还会因为超声波风速风向仪在工作时受到旋翼自身产生风场的影响而极大降低测量精度。故现有技术无法实现在多旋翼无人机上对风速、风向的精确监测。
考虑到能装在无人机上,这里我们选用了原理简单的超声波风速风向仪,发明了一种基于超声波风速风向仪的无人机风速风向测量的校正算法。利用超声波的传播特性来测量风速风向,具有测量精度高、使用寿命长的优点。它适用于多种环境,可靠性高。所以我们把这种测风方法与无人机相结合,可以充分发挥其优势。但因为传感器放置在无人机上,无人机前进或者有风吹来时机身平面会倾斜,超声波风速风向仪也随之倾斜,由于超声波风速风向仪仅测得与超声波风速风向仪平面平行的风,这势必造成测量结果的误差,因此要对其结果进行校正。
发明内容
本发明针对上述现有技术的不足,而提供一种基于超声波风速风向仪的无人机风速风向测量及校正算法。
为解决上述技术问题,本发明采用的技术方案是:
一种基于超声波风速风向仪的无人机风速风向测量及校正算法,包括两种状态的风速风向测量及校正算法,一种是无人机悬停状态下的风速风向测量及校正算法,另一种是无人机前进状态下的风速风向测量及校正算法,两种风速风向测量及校正算法都是利用四元数与姿态阵之间的关系来得出无人机的三个欧拉角,进而得到无人机的倾角,通过倾角的大小来对超声波风速风向仪测得的风速风向结果进行校正;
设有一参考坐标系R,简称R系,一刚体相对于R系作定点转动,定点为O;选取坐标系b与刚体固联,坐标系b简称b系,假设初始时刻b系与R系重合,设OA=r为起始位置向量,OA'=r'为旋转后的向量;根据欧拉定理,仅仅考虑初始时刻与最终时刻的位置,刚体从A位置转至A'位置等效成绕单位瞬轴转过θ1角度一次完成;将向量分解、旋转、合成可以得到:
r'=rcosθ1+(1-cosθ1)(u·r)u+u×rsinθ1; (1)
由三重矢积公式变换得:
r'=r+u×rsinθ1+(1-cosθ1)u×(u×r); (2)
则:
u×r=Ur;
u×(u×r)=U·Ur;
所以有:
令:
则式(3)可以写成:
r'=Dr; (5)
记初始时刻的刚体固联坐标系为b0,由于初始时刻刚体固联坐标系b与参考坐标系R重合,所以有:
在转动过程中位置向量和b系都和刚体固联,所以位置向量和b系的相对角位置始终不变,即有:
所以得到:
r=r'b; (8)
将式(8)带入式(5)得:
r'=Dr'b; (9)
该式说明D为b系到R系的坐标变换矩阵;
即:
令:
以q0,q1,q2,q3构造四元数:
设参考坐标系为导航坐标系n,与刚体固联的坐标系为机体坐标系,那么坐标变换矩阵就是姿态矩阵在机体轴坐标系下,无人机先绕x轴转动横滚角α,再绕y轴转动俯仰角β,最后再绕z轴转动航向角γ;则在此坐标系下的变换矩阵为:
对比式(13)和式(14)可得到三个姿态角:
β=-arcsinT31=-arcsin2(q1q3-q0q2);
其中,T11为cosβcosγ,T21为cosβsinγ,T31为-sinβ,T32为sinαcosβ,T33为cosαcosβ;
所以无人机的倾角θ为:
θ=arccos(cosαcosβ); (17)
由超声波测风原理可以得到风速ν为:
其中,t1、t2为x轴方向传播时间,t3、t4为y轴方向传播时间,
其中,k为整数,公式后面标了各种情况下的k取什么值。
当无人机处于悬停状态下时,无人机内内置的GPS测得的无人机速度大小为0,此时超声波风速风向仪所测得的风向与实际风向一致;将风分解为与无人机平面平行的分量和垂直的分量,垂直分量对超声波风速风向仪无影响,所以超声波风速风向仪测得的风为平行分量,测得的风向角与原风向角一致,所以修正的结果为:
当无人机处于前进状态时,无人机内内置的GPS测得无人机前进的速度大小为v1,无人机以一定的速度前进,等效于无人机悬停时,同样大小的风沿着无人机前进方向的反方向吹来;超声波风速风向仪测量的速度大小为v,此时的v合为无人机前进状态下等效的风速v1与实际要测的风速v2的合成风,则:
在机体坐标系下,无论无人机平面如何旋转都会过原点,所以无人机平面的方程为:
Ax+By+Cz=0
又因为其法向量为:
求得:
令z=0可得交线方程为:
Ax+By=0; (25)
实际风速大小为:
在这种情况下求出的实际风的风向角θ2为:
当B<0时的平面合风向分解示意图,此时要求的风向角为负,同理可得:
将θ2转化为[0,2π)内的角度,即:
其中,θ′2为求出的风向角转化为[0,2π)内的角;
随着时间的变化,某时刻的风向为:
与现有技术相比,本发明的有益效果在于:本发明的测量及其校正算法克服了因无人机前进或者有风吹来时机身平面倾斜导致的测量结果误差,具有测量精度高、使用寿命长的优点。
附图说明
图1为刚体旋转示意图;
图2为风向分解示意图;
图3为B>0时平面合风向分解示意图;
图4为B<0时平面合风向分解示意图。
其中,1-无人机平面,2-水平面,3-风向,4-风向沿无人机平面的分量,5-风向垂直无人机平面的分量,6-无人机平面与水平面的交线。
具体实施方式
现在结合附图对本发明作进一步详细的说明。这些附图均为简化的示意图,仅以示意方式说明本发明的基本结构,因此其仅显示与本发明有关的构成。
一种基于超声波风速风向仪的无人机风速风向测量及校正算法,包括两种状态的风速风向测量及校正算法,一种是无人机悬停状态下的风速风向测量及校正算法,另一种是无人机前进状态下的风速风向测量及校正算法,两种风速风向测量及校正算法都是利用四元数与姿态阵之间的关系来得出无人机的三个欧拉角,进而得到无人机的倾角,通过倾角的大小来对超声波风速风向仪测得的风速风向结果进行校正;
图1是刚体旋转示意图,设有一参考坐标系R,简称R系,一刚体相对于R系作定点转动,定点为O;选取坐标系b与刚体固联,坐标系b简称b系,假设初始时刻b系与R系重合,设OA=r为起始位置向量,OA'=r'为旋转后的向量;根据欧拉定理,仅仅考虑初始时刻与最终时刻的位置,刚体从A位置转至A'位置等效成绕单位瞬轴转过θ1角度一次完成;将向量分解、旋转、合成可以得到:
r'=rcosθ1+(1-cosθ1)(u·r)u+u×rsinθ1; (1)
由三重矢积公式变换得:
r'=r+u×rsinθ1+(1-cosθ1)u×(u×r); (2)
令:
则:
u×r=Ur;
u×(u×r)=U·Ur;
所以有:
令:
则式(3)可以写成:
r'=Dr; (5)
记初始时刻的刚体固联坐标系为b0,由于初始时刻刚体固联坐标系b与参考坐标系R重合,所以有:
在转动过程中位置向量和b系都和刚体固联,所以位置向量和b系的相对角位置始终不变,即有:
所以得到:
r=r'b; (8)
将式(8)带入式(5)得:
r'=Dr'b; (9)
该式说明D为b系到R系的坐标变换矩阵;
即:
令:
以q0,q1,q2,q3构造四元数:
设参考坐标系为导航坐标系n,与刚体固联的坐标系为机体坐标系,那么坐标变换矩阵就是姿态矩阵在机体轴坐标系下,无人机先绕x轴转动横滚角α,再绕y轴转动俯仰角β,最后再绕z轴转动航向角γ;则在此坐标系下的变换矩阵为:
对比式(13)和式(14)可得到三个姿态角:
β=-arcsinT31=-arcsin2(q1q3-q0q2);
其中,T11为cosβcosγ,T21为cosβsinγ,T31为-sinβ,T32为sinαcosβ,T33为cosαcosβ;
所以无人机的倾角θ为:
θ=arccos(cosαcosβ); (17)
由超声波测风原理可以得到风速ν为:
其中,t1、t2为x轴方向传播时间,t3、t4为y轴方向传播时间,
其中,k为整数,公式后面标了各种情况下的k取什么值。
图2为风向分解示意图,1为无人机平面,2为水平面,3为风向,4为风向沿无人机平面的分量,5为风向垂直无人机平面的分量,6为无人机平面与水平面的交线。
当无人机处于悬停状态下时,GPS测得的无人机速度大小为0,此时风向3即为实际的风向。将风分解为与无人机平面平行的分量和垂直的分量,垂直分量对超声波风速风向仪无影响,所以超声波风速风向仪测得的风为平行分量,测得的风向角与原风向角一致,所以修正的结果为:
当无人机处于前进状态时,GPS测得无人机前进的速度大小为v1,无人机以一定的速度前进,等效于无人机悬停时,同样大小的风沿着无人机前进方向的反方向吹来。超声波风速风向仪测量的速度大小为v,此时的v合为无人机前进状态下等效的风速v1与实际要测的风速v2的合成风,则:
在机体坐标系下,无论无人机平面如何旋转都会过原点,所以无人机平面的方程为:
Ax+By+Cz=0
又因为其法向量为:
求得:
令z=0可得交线方程为:
Ax+By=0; (25)
实际风速大小为:
在这种情况下求出的风向角为:
图4为B<0时的平面合风向分解示意图,此时要求的风向角为负,同理可得:
将θ2转化为[0,2π)内的角度,即
其中,θ′2为求出的风向角转化为[0,2π)内的角;
随着时间的变化,某时刻的风向为:
以上述依据本发明的理想实施例为启示,通过上述的说明内容,相关工作人员完全可以在不偏离本项发明技术思想的范围内,进行多样的变更以及修改。本项发明的技术性范围并不局限于说明书上的内容,必须要根据权利要求范围来确定其技术性范围。
Claims (3)
1.一种基于超声波风速风向仪的无人机风速风向测量及校正算法,包括两种状态的风速风向测量及校正算法,一种是无人机悬停状态下的风速风向测量及校正算法,另一种是无人机前进状态下的风速风向测量及校正算法,其特征在于:两种风速风向测量及校正算法都是利用四元数与姿态阵之间的关系来得出无人机的三个欧拉角,进而得到无人机的倾角,通过倾角的大小来对超声波风速风向仪测得的风速风向结果进行校正;
设有一参考坐标系R,简称R系,一刚体相对于R系作定点转动,定点为O;选取坐标系b与刚体固联,坐标系b简称b系,假设初始时刻b系与R系重合,设OA=r为起始位置向量,OA'=r'为旋转后的向量;根据欧拉定理,仅仅考虑初始时刻与最终时刻的位置,刚体从A位置转至A'位置等效成绕单位瞬轴转过θ1角度一次完成;将向量分解、旋转、合成可以得到:
r'=rcosθ1+(1-cosθ1)(u·r)u+u×rsinθ1;(1)
由三重矢积公式变换得:
r'=r+u×rsinθ1+(1-cosθ1)u×(u×r);(2)
令:
则:
u×r=Ur;
u×(u×r)=U·Ur;
所以有:
令:
则式(3)可以写成:
r'=Dr; (5)
记初始时刻的刚体固联坐标系为b0,由于初始时刻刚体固联坐标系b与参考坐标系R重合,所以有:
在转动过程中位置向量和b系都和刚体固联,所以位置向量和b系的相对角位置始终不变,即有:
所以得到:
r=r'b; (8)
将式(8)带入式(5)得:
r'=Dr'b; (9)
该式说明D为b系到R系的坐标变换矩阵;
即:
令:
以q0,q1,q2,q3构造四元数:
将q0,q1,q2,q3代入式(10)中进一步化简得:
设参考坐标系为导航坐标系n,与刚体固联的坐标系为机体坐标系,那么坐标变换矩阵就是姿态矩阵在机体轴坐标系下,无人机先绕x轴转动横滚角α,再绕y轴转动俯仰角β,最后再绕z轴转动航向角γ;则在此坐标系下的变换矩阵为:
对比式(13)和式(14)可得到三个姿态角:
β=-arcsinT31=-arcsin2(q1q3-q0q2);
其中,T11为cosβcosγ,T21为cosβsinγ,T31为-sinβ,T32为sinαcosβ,T33为cosαcosβ;
所以无人机的倾角θ为:
θ=arccos(cosαcosβ);(17)
由超声波测风原理可以得到风速v为:
其中,t1、t2为x轴方向传播时间,t3、t4为y轴方向传播时间,
其中,k为整数,公式后面标了各种情况下的k取什么值。
3.根据权利要求1所述的基于超声波风速风向仪的无人机风速风向测量及校正算法,其特征在于:当无人机处于前进状态时,无人机内内置的GPS测得无人机前进的速度大小为v1,无人机以一定的速度前进,等效于无人机悬停时,同样大小的风沿着无人机前进方向的反方向吹来;超声波风速风向仪测量的速度大小为v,此时的v合为无人机前进状态下等效的风速v1与实际要测的风速v2的合成风,则:
在机体坐标系下,无论无人机平面如何旋转都会过原点,所以无人机平面的方程为:
Ax+By+Cz=0
又因为其法向量为:
求得:
令z=0可得交线方程为:
Ax+By=0; (25)
实际风速大小为:
在这种情况下求出的实际风的风向角θ2为:
当B<0时的平面合风向分解示意图,此时要求的风向角为负,同理可得:
故综上所述,可得:
将θ2转化为[0,2π)内的角度,即:
其中,θ′2为求出的风向角转化为[0,2π)内的角;
随着时间的变化,某时刻的风向为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911278177.7A CN110988393A (zh) | 2019-12-12 | 2019-12-12 | 基于超声波风速风向仪的无人机风速风向测量及校正算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911278177.7A CN110988393A (zh) | 2019-12-12 | 2019-12-12 | 基于超声波风速风向仪的无人机风速风向测量及校正算法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110988393A true CN110988393A (zh) | 2020-04-10 |
Family
ID=70092984
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911278177.7A Pending CN110988393A (zh) | 2019-12-12 | 2019-12-12 | 基于超声波风速风向仪的无人机风速风向测量及校正算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110988393A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113092813A (zh) * | 2021-03-15 | 2021-07-09 | 南京航空航天大学 | 一种翼伞空降情境下自适应风速风向检测***和方法 |
CN113252294A (zh) * | 2021-06-16 | 2021-08-13 | 西南交通大学 | 一种跨海大桥空间风速风向测试***及监测方法 |
CN118169427A (zh) * | 2024-05-14 | 2024-06-11 | 浙江大学 | 飞行器环境测风装置、方法及温室气体检测装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160109475A1 (en) * | 2014-10-16 | 2016-04-21 | The Boeing Company | Methods and systems for airborne vehicle wind depiction |
CN105954819A (zh) * | 2016-05-24 | 2016-09-21 | 南京信息工程大学 | 基于无人机倾角检测的风速测量装置及操作方法 |
CN108152529A (zh) * | 2017-11-02 | 2018-06-12 | 成都飞机工业(集团)有限责任公司 | 一种基于飞行参数计算风速及风向的方法 |
CN109374924A (zh) * | 2018-10-25 | 2019-02-22 | 西北工业大学 | 一种基于四旋翼无人机倾斜角的横侧向风场估计方法 |
CN110134134A (zh) * | 2019-05-24 | 2019-08-16 | 南京信息工程大学 | 一种无人机悬停状态下的测风方法 |
-
2019
- 2019-12-12 CN CN201911278177.7A patent/CN110988393A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160109475A1 (en) * | 2014-10-16 | 2016-04-21 | The Boeing Company | Methods and systems for airborne vehicle wind depiction |
CN105954819A (zh) * | 2016-05-24 | 2016-09-21 | 南京信息工程大学 | 基于无人机倾角检测的风速测量装置及操作方法 |
CN108152529A (zh) * | 2017-11-02 | 2018-06-12 | 成都飞机工业(集团)有限责任公司 | 一种基于飞行参数计算风速及风向的方法 |
CN109374924A (zh) * | 2018-10-25 | 2019-02-22 | 西北工业大学 | 一种基于四旋翼无人机倾斜角的横侧向风场估计方法 |
CN110134134A (zh) * | 2019-05-24 | 2019-08-16 | 南京信息工程大学 | 一种无人机悬停状态下的测风方法 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113092813A (zh) * | 2021-03-15 | 2021-07-09 | 南京航空航天大学 | 一种翼伞空降情境下自适应风速风向检测***和方法 |
CN113092813B (zh) * | 2021-03-15 | 2022-01-04 | 南京航空航天大学 | 一种翼伞空降情境下自适应风速风向检测***和方法 |
CN113252294A (zh) * | 2021-06-16 | 2021-08-13 | 西南交通大学 | 一种跨海大桥空间风速风向测试***及监测方法 |
CN118169427A (zh) * | 2024-05-14 | 2024-06-11 | 浙江大学 | 飞行器环境测风装置、方法及温室气体检测装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110988393A (zh) | 基于超声波风速风向仪的无人机风速风向测量及校正算法 | |
CN110134134B (zh) | 一种无人机悬停状态下的测风方法 | |
CN101813465B (zh) | 非接触精密测转角的单目视觉测量方法 | |
CN101067554A (zh) | 补偿地磁传感器姿态误差的设备和方法 | |
CN107831515B (zh) | 水下定位方法和*** | |
CN104459728B (zh) | 一种基于gnss定位的磁偏角校准方法 | |
CN110926468A (zh) | 基于传递对准的动中通天线多平台航姿确定方法 | |
CN104504255B (zh) | 一种螺旋翼升力和阻力力矩的确定方法 | |
CN105841698A (zh) | 一种无需调零的auv舵角精确实时测量*** | |
CN103439727A (zh) | 一种地面坐标的测量方法 | |
CN107063244A (zh) | 一种飞行器飞行过程模拟方法 | |
CN105892498A (zh) | 基于三轴云台的目标凝视与扫视控制*** | |
CN110779514B (zh) | 面向仿生偏振导航辅助定姿的分级卡尔曼融合方法及装置 | |
CN106123917A (zh) | 考虑外杆臂效应的捷联惯导***罗经对准方法 | |
CN101793517A (zh) | 一种提高机载平台定姿精度的在线快速方法 | |
CN104614746B (zh) | 一种基于卫星接收机的航姿仪校准装置和方法 | |
CN105321186B (zh) | 一种基于正射投影的封闭空间地图采集方法及装置 | |
CN105953797A (zh) | 利用单轴陀螺仪、倾角仪和里程计的组合导航装置及方法 | |
CN109470274A (zh) | 一种车载光电经纬仪载车平台变形测量***及方法 | |
CN106968665A (zh) | 一种利用惯导***进行石油井测斜的方法 | |
CN111649738B (zh) | 微重力场下的加速度计初始姿态解算方法 | |
CN103557875B (zh) | 一种姿态测量***的fpga硬件实现方法 | |
CN112013716A (zh) | 一种电视导引头模拟***捕捉快速移动目标的方法 | |
DE102015121517A1 (de) | Verfahren und Vorrichtung zur Bestimmung eines Geschwindigkeitsvektors eines in der Umgegebung eines Fluggeräts herrschenden Windes, und Fluggerät | |
CN114322970B (zh) | 一种双陀螺仪寻北方法、***及存储介质 |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20200410 |