CN111695079A - 基于二维傅里叶变换的基尔霍夫近似下激光散射回波计算 - Google Patents
基于二维傅里叶变换的基尔霍夫近似下激光散射回波计算 Download PDFInfo
- Publication number
- CN111695079A CN111695079A CN202010528454.1A CN202010528454A CN111695079A CN 111695079 A CN111695079 A CN 111695079A CN 202010528454 A CN202010528454 A CN 202010528454A CN 111695079 A CN111695079 A CN 111695079A
- Authority
- CN
- China
- Prior art keywords
- dimensional
- fourier transform
- rough
- rough surface
- dimensional fourier
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 15
- 230000010287 polarization Effects 0.000 claims description 12
- 238000000342 Monte Carlo simulation Methods 0.000 claims description 9
- 239000000126 substance Substances 0.000 claims description 9
- 230000009467 reduction Effects 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000007689 inspection Methods 0.000 abstract description 2
- 238000000034 method Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 4
- 238000009795 derivation Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000012634 optical imaging Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
Images
Classifications
-
- 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/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
- G06F17/142—Fast Fourier transforms, e.g. using a Cooley-Tukey type algorithm
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- 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/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/145—Square transforms, e.g. Hadamard, Walsh, Haar, Hough, Slant transforms
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Remote Sensing (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Radar, Positioning & Navigation (AREA)
- Computer Networks & Wireless Communication (AREA)
- Geometry (AREA)
- Discrete Mathematics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Optical Radar Systems And Details Thereof (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明提供基于二维傅里叶变换的基尔霍夫近似下激光散射回波计算,属于粗糙目标激光散射回波分析与检查技术领域,该基于二维傅里叶变换的基尔霍夫近似下激光散射回波计算包括以下步骤;步骤S1:二维粗糙面建模;步骤S2:计算二维粗糙面高度起伏的方向斜率;步骤S3:推导基于二维傅里叶变换的基尔霍夫近似;步骤S4:粗糙凸目标(圆锥)的建模,本发明旨在解决现有技术中的基尔霍夫近似计算粗糙凸目标散射回波效率低的问题。
Description
技术领域
本发明属于粗糙目标激光散射回波分析与检查技术领域,具体涉及基于二维傅里叶变换的基尔霍夫近似下激光散射回波计算。
背景技术
在现有技术中,根据波动光学理论,系统孔径会限制光学成像系统的分辨率,合成孔径雷达利用运动模拟虚拟孔径合成,从而实现二维的高分辨率成像;激光雷达相对于微波、无线电雷达,由于其频率高波长短,在遥感探测的高度、空间分辨率、时间上的连续监测和测量精度等方面具有独特的优势。(逆)合成孔径激光雷达(SAL/ISAL)突破了衍射极限的限制,实现了超远距离上厘米级分辨率的雷达成像。但由于激光频率高、波长短的特点,对于大目标SAL/ISAL成像,计算散射回波就面临现实困难。提高目标散射回波的计算效率就十分重要。
现在计算散射回波的方法主要有积分方程方法、基尔霍夫近似、微扰法、小斜率近似和双尺度方法。其中,基尔霍分近似又称为切平面近似,将粗糙曲面用局部切平面代替,由菲涅耳反射定律获得切平面的总场,从而近似计算远区散射场。与入射激光波长相比,对应平面的水平尺寸很大,对于粗糙大目标,计算效率就很低。
发明内容
本发明的目的在于提供基于二维傅里叶变换的基尔霍夫近似下激光散射回波计算,旨在解决现有技术中的基尔霍夫近似计算粗糙凸目标散射回波效率低的问题。
为实现上述目的,本发明提供如下技术方案:
基于二维傅里叶变换的基尔霍夫近似下激光散射回波计算,包括以下步骤;
步骤S1:二维粗糙面建模;
步骤S2:计算二维粗糙面高度起伏的方向斜率;
步骤S3:推导基于二维傅里叶变换的基尔霍夫近似;
步骤S4:粗糙凸目标(圆锥)的建模。
作为本发明优选的方案,在步骤S1中包括二维粗糙面建模;
步骤S11:首先,利用蒙特卡洛方法产生了二维随机高斯粗糙面,二维高斯粗糙面的功率谱密度为:
式中,lx,ly分别为x,y方向的相关长度,δ为高度起伏的均方根,kx,ky分别为x,y方向的波数;
步骤S12:假设二维高斯粗糙面长度分别为Lx,Ly,等间隔的离散点数分别为M,N,相邻两点间的距离分别为Δx,Δy,则粗糙面上任一点处高度随机起伏为:
其中,
作为本发明优选的方案,
作为本发明优选的方案,在步骤S2中包括计算二维粗糙面高度起伏的方向斜率;
步骤S21:利用公式可以得到粗糙面上任一点的高度起伏f(xm,yn),现要求出粗糙面分别沿x,y方向的斜率zx,zy,利用斜率概念对公式求偏导;
作为本发明优选的方案,在步骤S3中包括推导基于二维傅里叶变换的基尔霍夫近似;
步骤S31:基尔霍夫近似又称切平面近似,即界面上任一点的场强由该点处的切平面反射波决定,将粗糙曲面用局部切平面代替,利用菲涅尔反射定律获得切平面的总场,从而近似得到远区的散射场,在小粗糙度的情况下,基本的散射回波表达式为:
其中,ai(i=0,1,2)为极化系数,Zx,Zy为高度起伏的方向斜率,表达式为和;
步骤S32:本发明采用VV极化,所以对应的极化系数为:
在直角坐标系下ds′=qdx′dy′/|qz|,因此散射场公式整理为:
步骤S33:二维离散Gabor变换的核函数为:
其中,式中LX和Ly分别x和y方向的窗宽度,根据傅里叶变换理论,对散射场进行二维卷积等价于对散射场二维傅里叶变换以及核函数的二维傅里叶变换在空间域求乘法,对所得到的结果进行傅里叶逆变换;
则散射场公式改写为:
作为本发明优选的方案,在步骤S4中包括粗糙凸目标(圆锥)的建模;
步骤S41:首先,利用蒙特卡洛方法生成二维随机的高斯粗糙面,然后,将光滑圆锥的锥面进行剖分,得到与已知高斯粗糙面大小相当的小面元,最后利用坐标变换将高斯粗糙面替换掉光滑小面元,从而得到粗糙圆锥。
与现有技术相比,本发明的有益效果是:
1、本方案中的基于二维傅里叶变换的基尔霍夫近似下激光散射回波计算,首先利用蒙特卡洛方法生成二维高斯粗糙面,并利用公式计算得到二维高斯粗糙面高度起伏的方向导数,然后对基尔霍夫近似公式进行整理,使其满足二维傅里叶变换关系,然后引入二维Gabor变换(即窗口选取高斯窗的窗口傅里叶变换),从而减弱边缘效应,增强边缘的图像特征,得到了基于二维傅里叶变换的基尔霍夫近似下散射回波公式,最后,对于粗糙凸目标(圆锥) 进行建模,从而利用推导公式去计算粗糙凸目标(圆锥)激光散射的回波,提高计算粗糙凸目标激光散射回波的效率,为粗糙凸目标(圆锥)的(逆) 合成孔径激光雷达(SAL/ISAL)成像提供了基础。
附图说明
附图用来提供对本发明的进一步理解,并且构成说明书的一部分,与本发明的实施例一起用于解释本发明,并不构成对本发明的限制。在附图中:
图1为本发明的高斯粗糙面斜率zx的频数直方图及概率密度函数图;
图2为本发明中的高斯粗糙面斜率zy的频数直方图及概率密度函数图;
图3为本发明中的粗糙面散射示意图图;
图4为本发明中的粗糙凸目标(圆锥)的建模过程图;
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例1
请参阅图1-4,本发明提供以下技术方案:
基于二维傅里叶变换的基尔霍夫近似下激光散射回波计算,包括以下步骤;
步骤S1:二维粗糙面建模;
步骤S2:计算二维粗糙面高度起伏的方向斜率;
步骤S3:推导基于二维傅里叶变换的基尔霍夫近似;
步骤S4:粗糙凸目标(圆锥)的建模。
在本发明的具体实施例中,首先利用蒙特卡洛方法生成二维高斯粗糙面,并利用公式计算得到二维高斯粗糙面高度起伏的方向导数,然后对基尔霍夫近似公式进行整理,使其满足二维傅里叶变换关系,然后引入二维Gabor变换(即窗口选取高斯窗的窗口傅里叶变换),从而减弱边缘效应,增强边缘的图像特征,得到了基于二维傅里叶变换的基尔霍夫近似下散射回波公式,最后,对于粗糙凸目标(圆锥)进行建模,从而利用推导公式去计算粗糙凸目标(圆锥)激光散射的回波。
具体的,请参阅图1-4,在步骤S1中包括二维粗糙面建模;
步骤S11:首先,利用蒙特卡洛方法产生了二维随机高斯粗糙面,二维高斯粗糙面的功率谱密度为:
式中,lx,ly分别为x,y方向的相关长度,δ为高度起伏的均方根,kx,ky分别为x,y方向的波数;
步骤S12:假设二维高斯粗糙面长度分别为Lx,Ly,等间隔的离散点数分别为M,N,相邻两点间的距离分别为Δx,Δy,则粗糙面上任一点处高度随机起伏为:
其中,
具体的,请参阅图1-2,在步骤S2中包括计算二维粗糙面高度起伏的方向斜率;
步骤S21:利用公式可以得到粗糙面上任一点的高度起伏f(xm,yn),现要求出粗糙面分别沿x,y方向的斜率zx,zy,利用斜率概念对公式求偏导;
在本发明的具体实施例中,以二维高斯粗糙面为例,利用公式和计算产生的斜率zx,zy,并求其频数直方图及概率密度函数,可以发现其满足高斯特性,验证其公式的正确性。
具体的,请参阅图3,在步骤S3中包括推导基于二维傅里叶变换的基尔霍夫近似;
步骤S31:基尔霍夫近似又称切平面近似,即界面上任一点的场强由该点处的切平面反射波决定,将粗糙曲面用局部切平面代替,利用菲涅尔反射定律获得切平面的总场,从而近似得到远区的散射场,在小粗糙度的情况下,基本的散射回波表达式为:
其中,ai(i=0,1,2)为极化系数,Zx,Zy为高度起伏的方向斜率,表达式为和;
步骤S32:本发明采用VV极化,所以对应的极化系数为:
在直角坐标系下ds′=qdx′dy′/|qz|,因此散射场公式整理为:
步骤S33:二维离散Gabor变换的核函数为:
其中,式中LX和Ly分别x和y方向的窗宽度,根据傅里叶变换理论,对散射场进行二维卷积等价于对散射场二维傅里叶变换以及核函数的二维傅里叶变换在空间域求乘法,对所得到的结果进行傅里叶逆变换;
则散射场公式改写为:
在本发明的具体实施例中,为了减弱边缘效应,增强边缘的图像特征,我们引入了Gabor变换(即窗口选取高斯窗的窗口傅里叶变换),选择选高斯窗的原因在于高斯窗谱的主瓣较宽,故频率分辨率低。
具体的,请参阅图4,在步骤S4中包括粗糙凸目标(圆锥)的建模;
步骤S41:首先,利用蒙特卡洛方法生成二维随机的高斯粗糙面,然后,将光滑圆锥的锥面进行剖分,得到与已知高斯粗糙面大小相当的小面元,最后利用坐标变换将高斯粗糙面替换掉光滑小面元,从而得到粗糙圆锥。
本发明的工作原理及使用流程:首先利用蒙特卡洛方法生成二维高斯粗糙面,并利用公式计算得到二维高斯粗糙面高度起伏的方向导数,然后对基尔霍夫近似公式进行整理,使其满足二维傅里叶变换关系,然后引入二维 Gabor变换(即窗口选取高斯窗的窗口傅里叶变换),从而减弱边缘效应,增强边缘的图像特征,得到了基于二维傅里叶变换的基尔霍夫近似下散射回波公式,最后,对于粗糙凸目标(圆锥)进行建模,从而利用推导公式去计算粗糙凸目标(圆锥)激光散射的回波。
最后应说明的是:以上所述仅为本发明的优选实施例而已,并不用于限制本发明,尽管参照前述实施例对本发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.基于二维傅里叶变换的基尔霍夫近似下激光散射回波计算,其特征在于:包括以下步骤;
步骤S1:二维粗糙面建模;
步骤S2:计算二维粗糙面高度起伏的方向斜率;
步骤S3:推导基于二维傅里叶变换的基尔霍夫近似;
步骤S4:粗糙凸目标(圆锥)的建模。
4.根据权利要求3所述的基于二维傅里叶变换的基尔霍夫近似下激光散射回波计算,其特征在于:在步骤S3中包括推导基于二维傅里叶变换的基尔霍夫近似;
步骤S31:基尔霍夫近似又称切平面近似,即界面上任一点的场强由该点处的切平面反射波决定,将粗糙曲面用局部切平面代替,利用菲涅尔反射定律获得切平面的总场,从而近似得到远区的散射场,在小粗糙度的情况下,基本的散射回波表达式为:
其中,ai(i=0,1,2)为极化系数,Zx,Zy为高度起伏的方向斜率,表达式为和;
步骤S32:本发明采用VV极化,所以对应的极化系数为:
在直角坐标系下ds′=qdx′dy′/|qz|,因此散射场公式整理为:
步骤S33:二维离散Gabor变换的核函数为:
其中,式中LX和Ly分别x和y方向的窗宽度,根据傅里叶变换理论,对散射场进行二维卷积等价于对散射场二维傅里叶变换以及核函数的二维傅里叶变换在空间域求乘法,对所得到的结果进行傅里叶逆变换;
则散射场公式改写为:
5.根据权利要求4所述的基于二维傅里叶变换的基尔霍夫近似下激光散射回波计算,其特征在于:在步骤S4中包括粗糙凸目标(圆锥)的建模;
步骤S41:首先,利用蒙特卡洛方法生成二维随机的高斯粗糙面,然后,将光滑圆锥的锥面进行剖分,得到与已知高斯粗糙面大小相当的小面元,最后利用坐标变换将高斯粗糙面替换掉光滑小面元,从而得到粗糙圆锥。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010528454.1A CN111695079A (zh) | 2020-06-11 | 2020-06-11 | 基于二维傅里叶变换的基尔霍夫近似下激光散射回波计算 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010528454.1A CN111695079A (zh) | 2020-06-11 | 2020-06-11 | 基于二维傅里叶变换的基尔霍夫近似下激光散射回波计算 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111695079A true CN111695079A (zh) | 2020-09-22 |
Family
ID=72480229
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010528454.1A Pending CN111695079A (zh) | 2020-06-11 | 2020-06-11 | 基于二维傅里叶变换的基尔霍夫近似下激光散射回波计算 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111695079A (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4334780A (en) * | 1979-06-29 | 1982-06-15 | Grumman Aerospace Corporation | Optical surface roughness detection method and apparatus |
CN109241556A (zh) * | 2018-07-27 | 2019-01-18 | 西安空间无线电技术研究所 | 一种随机粗糙体目标建模方法及存储介质 |
CN110489714A (zh) * | 2019-07-17 | 2019-11-22 | 西安理工大学 | 一种二维随机粗糙面散射统计矩的计算方法 |
CN111143757A (zh) * | 2019-12-10 | 2020-05-12 | 西安理工大学 | 一种二维非高斯粗糙体目标散射特性的计算方法 |
-
2020
- 2020-06-11 CN CN202010528454.1A patent/CN111695079A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4334780A (en) * | 1979-06-29 | 1982-06-15 | Grumman Aerospace Corporation | Optical surface roughness detection method and apparatus |
CN109241556A (zh) * | 2018-07-27 | 2019-01-18 | 西安空间无线电技术研究所 | 一种随机粗糙体目标建模方法及存储介质 |
CN110489714A (zh) * | 2019-07-17 | 2019-11-22 | 西安理工大学 | 一种二维随机粗糙面散射统计矩的计算方法 |
CN111143757A (zh) * | 2019-12-10 | 2020-05-12 | 西安理工大学 | 一种二维非高斯粗糙体目标散射特性的计算方法 |
Non-Patent Citations (2)
Title |
---|
任玉超等: "双频互相干函数及在粗糙面脉冲散射中的应用", 《中国科学(G辑:物理学 力学 天文学)》 * |
郭立新等: "二维带限分形粗糙面电磁散射的基尔霍夫近似", 《电子学报》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | Scattering and Doppler spectral analysis for two-dimensional linear and nonlinear sea surfaces | |
Fieber et al. | Analysis of full-waveform LiDAR data for classification of an orange orchard scene | |
Soriano et al. | A cutoff invariant two-scale model in electromagnetic scattering from sea surfaces | |
CN102402685A (zh) | 基于Gabor特征的三马尔可夫场SAR图像分割方法 | |
Guedes et al. | Observations of alongshore variability of swash motions on an intermediate beach | |
Nie et al. | Study of microwave backscattering from two-dimensional nonlinear surfaces of finite-depth seas | |
CN113376597A (zh) | 基于数字高程地图和gpu的复杂地形电磁散射快速仿真方法 | |
Wang et al. | Automatic extraction and evaluation of pavement three-dimensional surface texture using laser scanning technology | |
Wang et al. | Doppler simulation and analysis for 2-D sea surfaces up to Ku-band | |
Zhang et al. | A study on coastline extraction and its trend based on remote sensing image data mining | |
CN110208764A (zh) | 基于电磁散射计算的动态海洋大型场景回波仿真方法 | |
Wang et al. | Investigation on Doppler shift and bandwidth of backscattered echoes from a composite sea surface | |
CN111695079A (zh) | 基于二维傅里叶变换的基尔霍夫近似下激光散射回波计算 | |
Kim et al. | Radar images of rough surface scattering: Comparison of numerical and analytical models | |
Fung et al. | Kirchhoff model for a skewed random surface | |
Wang et al. | Investigation on the scattering characteristics of Gaussian beam from two dimensional dielectric rough surfaces based on the Kirchhoff approximation | |
Lyzenga et al. | Numerical calculations of radar scattering from sharply peaked ocean waves | |
Xue et al. | Simulating the Scattering Echo and Inverse Synthetic Aperture Lidar Imaging of Rough Targets | |
Di Martino et al. | Non-local means SAR despeckling based on scattering | |
Guissard | Multispectra for ocean-like random rough surface scattering | |
Shaw et al. | Green's function refinement as an approach to radar backscatter: General theory and applications to LGA scattering from the ocean | |
Schott et al. | Reflectivity of one-dimensional rough surfaces using the ray tracing technique with multiple reflections | |
Ouchi | A theory on the distribution function of backscatter radar cross section from ocean waves of individual wavelength | |
CN112764030B (zh) | 基于电磁与统计混合模型的海面sar图像杂波生成方法 | |
Shahrezaei et al. | Directional Wave Scattering Distribution Modes Analysis and Synthesis of Random Ocean Media Roughness for SAR Electromagnetic Interactions Using Feature Fusion in Dynamic Sea States: A Survey |
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 |
Application publication date: 20200922 |
|
WD01 | Invention patent application deemed withdrawn after publication |