CN105242252A - 基于图像匹配的下降轨聚束sar雷达定位方法 - Google Patents
基于图像匹配的下降轨聚束sar雷达定位方法 Download PDFInfo
- Publication number
- CN105242252A CN105242252A CN201510641039.6A CN201510641039A CN105242252A CN 105242252 A CN105242252 A CN 105242252A CN 201510641039 A CN201510641039 A CN 201510641039A CN 105242252 A CN105242252 A CN 105242252A
- Authority
- CN
- China
- Prior art keywords
- point
- sar
- setting
- radar
- ground
- 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.)
- Granted
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 123
- 238000000034 method Methods 0.000 title claims abstract description 47
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 claims abstract description 14
- 239000011159 matrix material Substances 0.000 claims abstract description 14
- 230000009466 transformation Effects 0.000 claims abstract description 14
- 238000012937 correction Methods 0.000 claims abstract description 5
- 238000013507 mapping Methods 0.000 claims abstract description 5
- 230000008878 coupling Effects 0.000 claims description 31
- 238000010168 coupling process Methods 0.000 claims description 31
- 238000005859 coupling reaction Methods 0.000 claims description 31
- OQCFWECOQNPQCG-UHFFFAOYSA-N 1,3,4,8-tetrahydropyrimido[4,5-c]oxazin-7-one Chemical compound C1CONC2=C1C=NC(=O)N2 OQCFWECOQNPQCG-UHFFFAOYSA-N 0.000 claims description 10
- 230000013011 mating Effects 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 4
- 238000005070 sampling Methods 0.000 abstract 1
- 230000008859 change Effects 0.000 description 6
- 208000004350 Strabismus Diseases 0.000 description 4
- 102000010029 Homer Scaffolding Proteins Human genes 0.000 description 3
- 108010077223 Homer Scaffolding Proteins Proteins 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 230000004807 localization Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000008439 repair process Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- NQQVFXUMIDALNH-UHFFFAOYSA-N picloram Chemical compound NC1=C(Cl)C(Cl)=NC(C(O)=O)=C1Cl NQQVFXUMIDALNH-UHFFFAOYSA-N 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
Classifications
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
- G01S13/9052—Spotlight mode
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9029—SAR image post-processing techniques specially adapted for moving target detection within a single SAR image or within multiple SAR images taken at the same time
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9004—SAR image acquisition techniques
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9027—Pattern recognition for feature extraction
-
- 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
- G01S7/415—Identification of targets based on measurements of movement associated with the target
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/02—Affine transformations
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Electromagnetism (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于图像匹配的下降轨聚束SAR雷达定位方法,其思路为:以地面点目标为原点,并基于SAR雷达的下降轨建立北天东坐标系下弹体定位的几何构型,依次得到地面点目标在成像平面的SAR成像平面图像和SAR雷达的实际位置A,再利用成像平面与地平面之间的几何失真校正映射函数,得到SAR成像匹配图,然后通过图像匹配算法,得到P个匹配点,根据该P个匹配点,得到SAR成像匹配图与设定的SAR成像基准图之间的仿射变换矩阵,并利用所述仿射变换矩阵,构建欧拉四面体,再根据设定的SAR雷达的采样频率及距离单元数,得到SAR雷达分别到地面点目标T、设定的参考点C及第p个匹配点Pp的距离,计算得到SAR雷达的实际位置A的坐标值。
Description
技术领域
本发明属于雷达图像目标定位技术领域,特别涉及一种基于图像匹配的下降轨聚束SAR雷达定位方法,适用于为弹体修偏提供准确的位置信息。
背景技术
导弹在末制导阶段因惯导误差导致导弹弹体的实际位置和惯导指示位置存在较大偏差,影响了导弹导引头对目标的准确定位,无法满足传统的光学导弹导引头的末段交接班精度要求,降低了导弹的精确打击性能。弹载合成孔径雷达(SyntheticApertureRadar,SAR)导引头在导弹下降段采用聚束SAR成像模式,通过对目标场景进行高分辨率二维聚焦成像,实时获取目标场景的SAR成像,并和设定的SAR基准图像进行匹配定位,为导弹弹体修偏提供准确的目标场景与导弹弹***置信息,从而解决光学导弹导引头末端交接班精度不满足要求的问题。
燕英等人在文献“弹载合成孔径雷达成像处理及定位误差分析”中提到基于惯导数据的SAR平台定位方法,该方法利用惯导数据中的斜视角和俯仰角进行SAR平台定位,其方法虽然简单,但国内使用的惯导数据误差普遍较大,导致SAR平台定位精度较低。
李亚超等人在文献“SAR末制导中导弹定位方法及分析”中提到利用导弹在末制导中不同的三个固定目标场景以及有关该三个固定目标场景各自的斜距和地距的已知参数对导弹弹体进行定位;该方法虽然分析了导弹规则飞行和导弹不规则飞行对定位精度的影响,但是该定位方法只适用于SAR平台工作在正侧视条带模式的情况,对于大斜视聚束模式,该方法并不适用。
李亚超等人在“高精度景象匹配下的高速SAR平台定位和测速可行性研究”中提到一种基于景象匹配的斜视SAR平台定位算法,该定位算法通过景象匹配,再利用数字地图高精度配准对斜视SAR平台进行定位与测速;虽然该方法的定位精度较高,但是仍然会受到地面匹配点位置的影响,当配准误差及估计误差确定,或地面匹配点位置较近时,也会使得对斜视SAR平台的定位精度下降。
杨立波等人在“合成孔径雷达景象匹配中制导导弹定位”中建立了基于匹配点距离差和多普勒差的定位模型,该定位模型对于下降轨,弹道倾角不为零,此时下降轨的速度矢量和中心斜距矢量的夹角,与速度矢量和瞬时斜距矢量的夹角不在同一平面上,即不能利用距离差和多普勒差的定位模型进行测角定位。此外,聚束模式下采用极坐标格式算法(PolarFormatAlgorithm,PFA)进行SAR成像时,得到的斜平面SAR成像为二维波数域图,使得多普勒测角公式不能直接使用,也限制了该方法的应用范围。
发明内容
针对以上现有技术存在的不足,本发明目的在于提出一种基于图像匹配的下降轨聚束SAR雷达定位方法,该方法基于匹配图与基准图进行景象匹配得到地面点目标信息,并在基准图中重构出以地面点目标为原点的北天东坐标系定位模型,该定位模型基于欧拉四面体几何构型,使得基于欧拉四面体几何构型的定位精度能够得到较好的保证。
为达到上述技术目的,本发明采用如下技术方案予以实现。
一种基于图像匹配的下降轨聚束SAR雷达定位方法,其特征在于,包括以下步骤:
步骤1,以地面点目标T为原点,并基于SAR雷达的下降轨建立北天东坐标系下弹体定位的几何构型,依次得到地面点目标T在成像平面的SAR成像平面图像和SAR雷达的实际位置点A;
步骤2,利用成像平面与地平面之间的几何失真校正映射函数,将地面点目标T在成像平面的SAR成像平面图像投影到地平面,得到SAR成像匹配图;
步骤3,设定SAR成像基准图,并在该设定的SAR成像基准图中设定参考点C,将SAR成像匹配图和设定的SAR成像基准图进行图像匹配,得到P个匹配点,再根据该P个匹配点,得到SAR成像匹配图与设定的SAR成像基准图之间的仿射变换矩阵;其中,P表示SAR成像匹配图与设定的SAR成像基准图进行匹配的匹配点总个数;
步骤4,在设定的SAR成像基准图中分别选取地面点目标T和设定的参考点C各自对应的位置点,并利用SAR成像匹配图与设定的SAR成像基准图之间的仿射变换矩阵,得到地面点目标T与设定的参考点C在SAR成像匹配图中各自对应的位置点坐标;
然后从P个匹配点中选取第i个匹配点Pi,并利用SAR成像匹配图与设定的SAR成像基准图之间的仿射变换矩阵,构建分别以SAR雷达的实际位置点A、第i个匹配点Pi、地面点目标T和设定的参考点C为顶点的欧拉四面体,进而得到所述欧拉四面体第一底边TPi的长度S1、欧拉四面体第二底边PiC的长度S2、欧拉四面体第三底边CT的长度S3;;其中,i∈{1,2,…,P},P表示SAR成像匹配图与设定的SAR成像基准图进行匹配的匹配点总个数,TPi表示地面点目标T与第i个匹配点Pi之间的距离,PiC表示第i个匹配点Pi与设定的参考点C之间的距离,CT表示设定的参考点C与地面点目标T之间的距离;
步骤5,分别设定SAR雷达的采样频率和SAR雷达的距离单元数,根据设定的SAR雷达的采样频率及设定的SAR雷达的距离单元数,得到SAR雷达分别到地面点目标T、设定的参考点C和第i个匹配点Pi的距离,再根据所述欧拉四面体第一底边TPi的长度S1、欧拉四面体第二底边PiC的长度S2、欧拉四面体第三底边CT的长度S3,计算得到SAR雷达的实际位置点A的坐标值。
本发明的有益效果:
第一,本发明方法基于欧拉四面体几何构型,通过图像匹配算法得到的高精度仿射变换矩阵和基于图像的雷达测距保证欧拉四面体中各点的位置精度,并利用图像匹配获取地面点目标的高精度相对位置信息和SAR雷达到各点目标的高精度斜距信息,通过构建欧拉四面体几何模型,解算出以地面点目标为原点构建的北天东坐标系下的弹体实际位置信息,使得本发明具备较高的定位精度;
第二,本发明方法只需要单幅匹配图就能够实现高精度定位,不需要多幅图进行数据融合提高精度。同时,本发明方法不依赖惯导数据,因此对惯导精度没有要求,具有较好的应用价值;
第三,本发明方法只需要一个匹配点,且该匹配点可以任意选取,其误差大小不会影响定位精度,能够有效降低图像匹配的要求,进一步提高了北天东坐标系下弹体定位的几何构型的鲁棒性;
第四,本发明方法中北天东坐标系下弹体定位的几何构型中的参考点可以任意选取,也可以不是成像场景中的点,使得本发明方法适用于小场景定位,解决了雷达末制导定位中需要大场景聚束成像的要求;
第五,本发明方法中SAR雷达的实际位置点A的输出是基于地面点目标为原点建立的北天东坐标系下的雷达三维位置,有效减少了坐标转换的步骤,并可以直接传给导弹的控制***进行弹***姿调整。
附图说明
下面结合附图和具体实施方式对本发明作进一步详细说明。
图1是本发明的一种基于图像匹配的下降轨聚束SAR雷达定位方法的流程示意图;
图2是以地面点目标为原点建立的北天东坐标系下弹体定位几何构型示意图;其中,所述北天东坐标系下弹体定位几何构型为X-Y-Z三维坐标系,A表示SAR雷达的实际位置点,T表示设定的SAR成像基准图中选取的地面点目标,C表示设定的参考点,建立以地面点目标T为原点的北天东坐标系,O表示SAR雷达的实际位置点投影到地面上的位置点,且其坐标为(x,y,0),设定的SAR成像基准图的中心点坐标为(Δx,Δy,0),Pi表示从P个匹配点中选取的第i个匹配点,P表示SAR成像匹配图与设定的SAR成像基准图进行匹配的匹配点总个数,F表示O点在Y坐标轴上的投影点,G表示设定的参考点C在X轴上的投影点,H表示设定的参考点C在Y轴上的投影点,OF表示SAR雷达的实际位置点A与Y轴的距离在地面上的投影,CG表示设定的参考点C与X轴之间的距离,E表示CG的延长线与OF的交点。
具体实施方式:
参照图1,是本发明的一种基于图像匹配的下降轨聚束SAR雷达定位方法的流程示意图,该种基于图像匹配的下降轨聚束SAR雷达定位方法的流程示意图,包括以下步骤:
步骤1,以地面点目标T为原点,并基于SAR雷达的下降轨建立北天东坐标系下弹体定位的几何构型,依次得到地面点目标T在成像平面的SAR成像平面图像和SAR雷达的实际位置点A。
具体地,参照图2,是以地面点目标T为原点建立的北天东坐标系下弹体定位几何构型示意图;其中,A为SAR雷达的实际位置;通过以地面点目标T为原点的北天东坐标系定位模型可知,建立的以地面点目标T为原点的北天东坐标系是可调整的;为了降低SAR雷达的定位误差,应尽可能地保证SAR雷达与地面点目标T的斜距在地面上的投影,和所述北天东坐标系X轴的夹角尽量小,所述夹角可以为15°左右的小夹角。SAR雷达工作在聚束模式对地面点目标T进行聚焦成像,依次得到地面点目标T在成像平面的SAR成像平面图像和SAR雷达的实际位置A。
步骤2,利用SAR雷达成像算法中成像平面与地平面之间的几何失真校正映射函数,将地面点目标T在成像平面的SAR成像平面图像投影到地平面,得到SAR成像匹配图。
具体地,由于导弹使用SAR雷达得到的SAR成像平面图像在斜平面上,因此为了SAR雷达的实际位置A的定位需要,利用SAR雷达成像算法中成像平面与地平面之间的几何失真校正映射函数,将所述SAR成像平面图像投影到地平面,得到SAR成像匹配图。
步骤3,设定SAR成像基准图,并在该SAR成像基准图中设定参考点C,根据图像匹配算法将SAR成像匹配图和设定的SAR成像基准图进行图像匹配,得到P个匹配点,再根据该P个匹配点,得到SAR成像匹配图与设定的SAR成像基准图之间的仿射变换矩阵;其中,P表示SAR成像匹配图与设定的SAR成像基准图进行匹配的匹配点总个数。
步骤4,在设定的SAR成像基准图中分别选取地面点目标T和设定的参考点C各自对应的位置点,并利用SAR成像匹配图与设定的SAR成像基准图之间的仿射变换矩阵,得到地面点目标T与设定的参考点C在SAR成像匹配图中各自对应的位置点坐标;
然后从P个匹配点中选取第i个匹配点Pi,并利用SAR成像匹配图与设定的SAR成像基准图之间的仿射变换矩阵,构建分别以SAR雷达的实际位置点A、第i个匹配点Pi、地面点目标T和设定的参考点C为顶点的欧拉四面体,进而得到所述欧拉四面体第一底边TPi的长度S1、欧拉四面体第二底边PiC的长度S2、欧拉四面体第三底边CT的长度S3;其中,i∈{1,2,…,P},P表示SAR成像匹配图与设定的SAR成像基准图进行匹配的匹配点总个数,TPi表示地面点目标T与第i个匹配点Pi之间的连线,PiC表示第i个匹配点Pi与设定的参考点C之间的连线,CT表示设定的参考点C与地面点目标T之间的连线。
具体地,在设定的SAR成像基准图中分别选取地面点目标T和设定的参考点C各自对应的位置点,并利用SAR成像匹配图与设定的SAR成像基准图之间的仿射变换矩阵,分别得到地面点目标T与设定的参考点C在SAR成像匹配图中各自对应的位置点坐标;
然后从P个匹配点中选取第i个匹配点Pi,并利用SAR成像匹配图与设定的SAR成像基准图之间的仿射变换矩阵,构建分别以SAR雷达的实际位置点A、第i个匹配点Pi、地面点目标T和设定的参考点C为顶点的欧拉四面体。
由于第i个匹配点Pi、地面点目标T和设定的参考点C均能在SAR匹配成像图中找到各自对应的位置坐标;因此可以通过所述第i个匹配点Pi、地面点目标T和设定的参考点C分别在SAR成像匹配图中各自对应的位置点坐标,分别得到所述欧拉四面体第一底边TPi的长度S1、欧拉四面体第二底边PiC的长度S2、欧拉四面体第三底边CT的长度S3;其中,i∈{1,2,…,P},P表示SAR成像匹配图与设定的SAR成像基准图进行匹配的匹配点总个数,TPi表示地面点目标T与第i个匹配点Pi之间的距离,PiC表示第i个匹配点Pi与设定的参考点C之间的距离,CT表示设定的参考点C与地面点目标T之间的距离。
步骤5,分别设定SAR雷达的采样频率和SAR雷达的距离单元数,并根据设定的SAR雷达的采样频率及设定的SAR雷达的距离单元数,得到SAR雷达分别到地面点目标T、设定的参考点C和第i个匹配点Pi的距离,再根据所述欧拉四面体第一底边TPi的长度S1、欧拉四面体第二底边PiC的长度S2、欧拉四面体第三底边CT的长度S3,计算得到SAR雷达的实际位置点A的坐标值。
具体地,在实际应用过程中,SAR雷达的实际位置点A到地面点目标T的距离还需要考虑到地面点的数字高程模型(DEM)。对于下降轨,SAR雷达的波束擦地脚很大,并且在欧拉四面体构型中,SAR雷达的实际位置点A到地面点目标T的距离近似等于基于图像测距得到的距离加上地面点目标T的数字高程模型(DEM)值。
5a)分别设定SAR雷达的采样频率和SAR雷达的距离单元数,并根据设定的SAR雷达的采样频率及设定的SAR雷达的距离单元数,得到SAR雷达分别到地面点目标T、设定的参考点C和第i个匹配点Pi的距离。
具体地,SAR雷达的实际位置点A到地面点目标T的距离为R1、SAR雷达的实际位置点A到设定的参考点C的距离为R2,SAR雷达的实际位置点A到第i个匹配点Pi的距离为R3,其表达式分别为:
其中,DEM1表示设定的参考点C的当地数字高程模型值,DEM2表示第i个匹配点Pi的当地数字高程模型值,DEM3表示地面点目标T的当地数字高程模型值,R1′表示地面点目标T对应的斜距长度,R2′表示设定的参考点C对应的斜距长度,R3′表示第i个匹配点Pi对应的斜距长度。
5b)根据SAR雷达分别到地面点目标T、设定的参考点C和第i个匹配点Pi的距离,以及所述欧拉四面体第一底边TPi的长度S1、欧拉四面体第二底边PiC的长度S2、欧拉四面体第三底边CT的长度S3,计算得到SAR雷达的实际位置点A的坐标值。
具体地,所求SAR雷达的实际位置点A的坐标设为(x,y,h),
其中,V表示欧拉四面体的体积,S表示欧拉四面体的底面积。
参照图2,是图2是以地面点目标为原点建立的北天东坐标系下弹体定位几何构型示意图,图2中的ATCPp为所述欧拉四面体,通过图2中所述欧拉四面体的几何关系,所述欧拉四面体的体积V的计算公式如下所示:
利用海伦公式,所述欧拉四面体的底面三角形TCPi的面积如下式所示:
其中,R1表示SAR雷达的实际位置点A到地面点目标T的距离,R2表示SAR雷达的实际位置点A到设定的参考点C的距离,R3表示SAR雷达的实际位置点A到第i个匹配点Pi的距离,s表示所述欧拉四面体的底面三角形TCPi的周长的一半,S1表示欧拉四面体第一底边TPi的长度,S2表示欧拉四面体第二底边PiC的长度,S3表示欧拉四面体第三底边CT的长度,则SAR雷达的实际位置点A的高度h可通过所述欧拉四面体的体积公式求得,如下式所示:
然后分别计算SAR雷达的实际位置点A的长度x和SAR雷达的实际位置点A的宽度y。
由图2可知有如下关系:
OT2-OC2-CT2=(OF2+FT2)-(OE2+EC2)-(CH2+CG2)
对上式展开,可得:
(R1 2-h2)-(R2 2-h2)-S3 2=(x2+y2)-[(x-Δx)2+(y-Δy)2]-(Δx2+Δy2)
其中,设定的SAR成像基准图的中心点坐标为(Δx,Δy,0),OT表示SAR雷达的实际位置点A与地面点目标T的距离在地面上的投影,OC表示SAR雷达的实际位置点A与设定的参考点C的距离在地面上的投影,CT表示设定的参考点C到地面点目标T的距离,OF表示SAR雷达的实际位置点A与Y轴的距离在地面上的投影,FT表示SAR雷达的实际位置点A与X轴的距离在地面的投影,CH表示设定的参考点C与Y轴之间的距离,CG表示设定的参考点C与X轴之间的距离,E表示CG的延长线与OF的交点,OE表示清晰化分析过程的第一辅助线,EC表示清晰化分析过程的第二辅助线,h表示SAR雷达的实际位置点A的高度,x表示SAR雷达的实际位置点A的长度,y表示SAR雷达的实际位置点A的宽度,Δx表示设定的SAR成像基准图的中心点的长度,Δy表示设定的SAR成像基准图的中心点的宽度,R1表示SAR雷达的实际位置点A到地面点目标T的距离,R2表示SAR雷达的实际位置点A到设定的参考点C的距离,R3表示SAR雷达的实际位置点A到第i个匹配点Pi的距离。
由上式可得:
其中
u表示为了表述方便而设定的一个变量,R1g表示SAR雷达与地面点目标T的距离在地面上的投影,R2g表示SAR雷达与设定的参考点C的距离在地面上的投影。
由上述公式,得到:
进而得到式解y1和y2,然后读取惯导的弹体低精度位置坐标,粗略确定SAR雷达的实际位置点A的宽度y,即如果x≥y,则y=min(y1,y2);如果x<y,y=max(y1,y2)。
根据勾股定理,得到SAR雷达的实际位置点A坐标的长度x为:
通过上述计算,得到SAR雷达的实际位置点A的坐标值。
特别地,SAR雷达的实际位置点A的坐标值关系式表明了SAR雷达的实际位置点A的坐标与设定的SAR成像基准图的中心点之间是有函数关系的,并且,SAR雷达的实际位置点A的长度x和SAR雷达的实际位置点A的宽度y关于Δx和Δy均具有对称性,进而关于Δx和Δy的偏导也是一致的,此处求y关于Δx和Δy的偏导,所述函数关系式具体为:
其中,∠XTO表示SAR雷达的实际位置点与地面点目标T的斜距在地面的投影与X轴的夹角,x表示SAR雷达的实际位置点A的长度,y表示SAR雷达的实际位置点A的宽度,Δx表示设定的SAR成像基准图的中心点的长度,Δy表示设定的SAR成像基准图的中心点的宽度,R1g表示SAR雷达与地面点目标T的距离在地面上的投影,表示求微分。
通过所述函数关系式中SAR雷达的实际位置点A的长度x和SAR雷达的实际位置点A的宽度y分别关于Δx和Δy求偏导,得到x和y分别关于Δx和Δy的变化率,所述变化率越小,意味着当Δx和Δy的值变化时,x和y的变化情况不大。
在本发明中,所述x和y,以及Δx和Δy通过所述函数关系式进行设定的参考点C的选取;为了实现SAR雷达的实际位置点A的精确定位,所述函数关系中两个偏导数的值就需要越小越好,这样能够使得设定的参考点C的选取误差对SAR雷达的实际位置点A的坐标值带来的影响最小化;而且所述函数关系式表明,设定的参考点C的选取误差会始终存在且不受其坐标位置的影响,并且所述设定的参考点C的选取会直接影响SAR雷达的实际位置点的坐标值;因此根据惯导数据思想得出,当设定的参考点C为SAR雷达星下点,即SAR雷达的实际位置点A投影到地面上的点时,设定的参考点C的选取误差带来的影响最小,SAR雷达的实际位置点A的定位精度更好。
本发明效果可以通过下述仿真实验得到进一步说明。
为验证本发明方法对SAR雷达的三维定位精度,下面通过仿真实验进行验证。
(一)仿真条件
SAR雷达的实际位置点与地面点目标的实际距离为20984m,本发明选取的设定的参考点C相对于地面点目标T,其距离差范围为(8000m,12000m)的大间距。
表1为仿真实验参数表,其中,R1表示雷达位置点与地面点目标的距离,∠XTO表示SAR雷达的实际位置点与地面点目标T的斜距在地面的投影与X轴的夹角,Δx表示设定的参考点C在X轴上的距离,Δy表示设定的参考点C在Y轴上的距离,x表示SAR雷达的实际位置点在X轴上的距离,y表示SAR雷达的实际位置点在Y轴上的距离,z表示SAR雷达的实际位置点在Z轴上的距离。
表1
(二)仿真内容
使用传统方法得到的定位误差分析如表2所示;其中,Δx0表示使用传统方法得到的SAR雷达位置点的长度与SAR雷达位置点的实际长度的差值,Δy0表示使用传统方法得到的SAR雷达位置点的宽度与SAR雷达位置点的实际宽度的差值,Δz0表示使用传统方法得到的SAR雷达位置点的高度与SAR雷达位置点的实际高度的差值。
表2
使用本发明方法得到的定位误差分析如表3所示;其中,Δxr表示使用本发明方法得到的SAR雷达位置点的长度与SAR雷达位置点的实际长度的差值,Δyr表示使用本发明方法得到的SAR雷达位置点的宽度与SAR雷达位置点的实际宽度的差值,Δzr表示使用本发明方法得到的SAR雷达位置点的高度与SAR雷达位置点的实际高度的差值。
表3
(三)结果分析
由于第i个匹配点对于该模型测角误差没有影响,因此第i个匹配点坐标误差未计入仿真模型,该结论也通过仿真实验验证;并且,分析表3可知,设定的参考点C的坐标默认准确,未加入误差,此结果表明,SAR雷达与地面点目标的斜距在地面的投影,和所述北天东坐标系轴的夹角测角精度能够进一步提高到0.1°,甚至小于0.1°。
通过对比表2和表3,我们可以得到以下两点结论:第一,在不同分辨率和估计误差情况下,本文的定位精度均优于传统的方法;第二,本发明方法通过选取合适的参考点,使得定位精度受估计误差和图像分辨率影响较小,而传统方法的定位精度的影响较大。
综上所述,仿真实验验证了本发明的正确性,有效性和可靠性。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围;这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (5)
1.一种基于图像匹配的下降轨聚束SAR雷达定位方法,其特征在于,包括以下步骤:
步骤1,以地面点目标T为原点,并基于SAR雷达的下降轨建立北天东坐标系下弹体定位的几何构型,依次得到地面点目标T在成像平面的SAR成像平面图像和SAR雷达的实际位置点A;
步骤2,利用成像平面与地平面之间的几何失真校正映射函数,将地面点目标T在成像平面的SAR成像平面图像投影到地平面,得到SAR成像匹配图;
步骤3,设定SAR成像基准图,并在该设定的SAR成像基准图中设定参考点C,将SAR成像匹配图和设定的SAR成像基准图进行图像匹配,得到P个匹配点,再根据该P个匹配点,得到SAR成像匹配图与设定的SAR成像基准图之间的仿射变换矩阵;其中,P表示SAR成像匹配图与设定的SAR成像基准图进行匹配的匹配点总个数;
步骤4,在设定的SAR成像基准图中分别选取地面点目标T和设定的参考点C各自对应的位置点,并利用SAR成像匹配图与设定的SAR成像基准图之间的仿射变换矩阵,得到地面点目标T与设定的参考点C在SAR成像匹配图中各自对应的位置点坐标;
然后从P个匹配点中选取第i个匹配点Pi,并利用SAR成像匹配图与设定的SAR成像基准图之间的仿射变换矩阵,构建分别以SAR雷达的实际位置点A、第i个匹配点Pi、地面点目标T和设定的参考点C为顶点的欧拉四面体,进而得到所述欧拉四面体第一底边TPi的长度S1、欧拉四面体第二底边PiC的长度S2、欧拉四面体第三底边CT的长度S3;其中,i∈{1,2,…,P},P表示SAR成像匹配图与设定的SAR成像基准图进行匹配的匹配点总个数,TPi表示地面点目标T与第i个匹配点Pi之间的距离,PiC表示第i个匹配点Pi与设定的参考点C之间的距离,CT表示设定的参考点C与地面点目标T之间的距离;
步骤5,分别设定SAR雷达的采样频率和SAR雷达的距离单元数,根据设定的SAR雷达的采样频率及设定的SAR雷达的距离单元数,得到SAR雷达分别到地面点目标T、设定的参考点C和第i个匹配点Pi的距离,再根据所述欧拉四面体第一底边TPi的长度S1、欧拉四面体第二底边PiC的长度S2、欧拉四面体第三底边CT的长度S3,计算得到SAR雷达的实际位置点A的坐标值。
2.如权利要求1所述的一种基于图像匹配的下降轨聚束SAR雷达定位方法,其特征在于,在步骤5中,所述得到SAR雷达分别到地面点目标T、设定的参考点C及第i个匹配点Pi的距离,分别为SAR雷达的实际位置点A到地面点目标T的距离R1、SAR雷达的实际位置点A到设定的参考点C的距离R2,SAR雷达的实际位置点A到第i个匹配点Pi的距离为R3,其表达式分别为:
其中,DEM1表示设定的参考点C的当地数字高程模型值,DEM2表示第i个匹配点Pi的当地数字高程模型值,DEM3表示地面点目标T的当地数字高程模型值,R1′表示地面点目标T对应的斜距长度,R2′表示设定的参考点C对应的斜距长度,R3′表示第i个匹配点Pi对应的斜距长度,i∈{1,2,…,P},P表示SAR成像匹配图与设定的SAR成像基准图进行匹配的匹配点总个数。
3.如权利要求1所述的一种基于图像匹配的下降轨聚束SAR雷达定位方法,其特征在于,在步骤5中,所述得到SAR雷达的实际位置点A的坐标值,所求SAR雷达的实际位置点A的坐标设为(x,y,h),
其中,V表示欧拉四面体的体积,S表示欧拉四面体的底面积;并且,所述欧拉四面体的体积V的计算公式如下所示:
利用海伦公式,所述欧拉四面体的底面三角形TCPi的面积表达式为:
其中,R1表示SAR雷达的实际位置点A到地面点目标T的距离,R2表示SAR雷达的实际位置点A到设定的参考点C的距离,R3表示SAR雷达的实际位置点A到第i个匹配点Pi的距离,s表示所述欧拉四面体的底面三角形TCPi的周长的一半,S1表示欧拉四面体第一底边TPi的长度,S2表示欧拉四面体第二底边PiC的长度,S3表示欧拉四面体第三底边CT的长度,则SAR雷达的实际位置点A的高度h可通过所述欧拉四面体的体积公式求得,其表达式为:
然后分别计算SAR雷达的实际位置点A的长度x和SAR雷达的实际位置点A的宽度y;
由所述欧拉四面体可知有如下关系:
OT2-OC2-CT2=(OF2+FT2)-(OE2+EC2)-(CH2+CG2)
对上式展开,可得:
(R1 2-h2)-(R2 2-h2)-S3 2=(x2+y2)-[(x-△x)2+(y-△y)2]-(△x2+△y2)
其中,设定的SAR成像基准图的中心点坐标为(△x,△y,0),OT表示SAR雷达的实际位置点A与地面点目标T的距离在地面上的投影,OC表示SAR雷达的实际位置点A与设定的参考点C的距离在地面上的投影,CT表示设定的参考点C到地面点目标T的距离,OF表示SAR雷达的实际位置点A与Y轴的距离在地面上的投影,FT表示SAR雷达的实际位置点A与X轴的距离在地面的投影,CH表示设定的参考点C与Y轴之间的距离,CG表示设定的参考点C与X轴之间的距离,E表示CG的延长线与OF的交点,OE表示清晰化分析过程的第一辅助线,EC表示清晰化分析过程的第二辅助线,h表示SAR雷达的实际位置点A的高度,x表示SAR雷达的实际位置点A的长度,y表示SAR雷达的实际位置点A的宽度,R1表示SAR雷达的实际位置点A到地面点目标T的距离,R2表示SAR雷达的实际位置点A到设定的参考点C的距离,R3表示SAR雷达的实际位置点A到第i个匹配点Pi的距离;
由上式可得:
其中,u表示为了表述方便而设定的一个变量,R1g表示SAR雷达与地面点目标T的距离在地面上的投影,R2g表示SAR雷达与设定的参考点C的距离在地面上的投影;
由上述公式,得到:
进而得到式解y1和y2,然后读取惯导的弹体低精度位置坐标,粗略确定SAR雷达的实际位置点A的宽度y,即如果x≥y,则y=min(y1,y2);如果x<y,y=max(y1,y2);
根据勾股定理,得到SAR雷达的实际位置点A坐标的长度x为:
通过上述计算,得到SAR雷达的实际位置点A的坐标值。
4.如权利要求1或权利要求2或权利要求3所述的一种基于图像匹配的下降轨聚束SAR雷达定位方法,其特征在于,所述设定的参考点C,将其设定为SAR雷达星下点。
5.如权利要求4所述的一种基于图像匹配的下降轨聚束SAR雷达定位方法,其特征在于,所述SAR雷达星下点,是根据SAR雷达的实际位置点A的坐标值和设定的SAR成像基准图的中心点坐标值的函数关系进行选取的,所述函数关系为:
其中,∠XTO表示SAR雷达的实际位置点与地面点目标T的斜距在地面的投影与X轴的夹角,x表示SAR雷达的实际位置点A的长度,y表示SAR雷达的实际位置点A的宽度,△x表示设定的SAR基准图的中心点的长度,△y表示设定的SAR基准图的中心点的宽度,R1g表示SAR雷达与地面点目标T的距离在地面上的投影,表示求微分。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510641039.6A CN105242252B (zh) | 2015-09-30 | 2015-09-30 | 基于图像匹配的下降轨聚束sar雷达定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510641039.6A CN105242252B (zh) | 2015-09-30 | 2015-09-30 | 基于图像匹配的下降轨聚束sar雷达定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105242252A true CN105242252A (zh) | 2016-01-13 |
CN105242252B CN105242252B (zh) | 2017-09-12 |
Family
ID=55039959
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510641039.6A Active CN105242252B (zh) | 2015-09-30 | 2015-09-30 | 基于图像匹配的下降轨聚束sar雷达定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105242252B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106353756A (zh) * | 2016-11-25 | 2017-01-25 | 西安电子科技大学 | 基于图像匹配的下降轨聚束合成孔径雷达定位方法 |
CN108594223A (zh) * | 2018-03-26 | 2018-09-28 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 机载sar图像目标定位方法 |
CN110244267A (zh) * | 2019-05-21 | 2019-09-17 | 成都航天科工微电子***研究院有限公司 | 弹载毫米波切向前视sar三维成像模型性能分析方法 |
CN111060910A (zh) * | 2019-12-11 | 2020-04-24 | 西安电子科技大学 | 一种基于地形-影像匹配的InSAR载体反定位 |
CN115166680A (zh) * | 2022-09-07 | 2022-10-11 | 中国科学院空天信息创新研究院 | 地面地物点的几何定位方法、装置、设备及介质 |
CN115638725A (zh) * | 2022-10-26 | 2023-01-24 | 成都清正公路工程试验检测有限公司 | 一种基于自动测量***的目标点位测量方法 |
CN115932823A (zh) * | 2023-01-09 | 2023-04-07 | 中国人民解放军国防科技大学 | 基于异源区域特征匹配的飞行器对地面目标定位方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7511656B2 (en) * | 2006-02-10 | 2009-03-31 | Raytheon Company | Spotlight synthetic aperture radar (SAR) system and method for generating a SAR map in real-time using a modified polar format algorithm (PFA) approach |
CN102590797A (zh) * | 2012-02-17 | 2012-07-18 | 张朝路 | 机载聚束式sar的高效成像处理方法 |
CN104062657A (zh) * | 2014-05-30 | 2014-09-24 | 西安电子科技大学 | 一种合成孔径雷达sar的广义极坐标成像方法 |
CN104678393A (zh) * | 2015-01-30 | 2015-06-03 | 南京航空航天大学 | 斜视滑动聚束sar的子孔径波数域成像方法 |
CN104865571A (zh) * | 2015-04-02 | 2015-08-26 | 西安电子科技大学 | 多通道多子带滑动聚束模式sar成像方法 |
-
2015
- 2015-09-30 CN CN201510641039.6A patent/CN105242252B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7511656B2 (en) * | 2006-02-10 | 2009-03-31 | Raytheon Company | Spotlight synthetic aperture radar (SAR) system and method for generating a SAR map in real-time using a modified polar format algorithm (PFA) approach |
CN102590797A (zh) * | 2012-02-17 | 2012-07-18 | 张朝路 | 机载聚束式sar的高效成像处理方法 |
CN104062657A (zh) * | 2014-05-30 | 2014-09-24 | 西安电子科技大学 | 一种合成孔径雷达sar的广义极坐标成像方法 |
CN104678393A (zh) * | 2015-01-30 | 2015-06-03 | 南京航空航天大学 | 斜视滑动聚束sar的子孔径波数域成像方法 |
CN104865571A (zh) * | 2015-04-02 | 2015-08-26 | 西安电子科技大学 | 多通道多子带滑动聚束模式sar成像方法 |
Non-Patent Citations (4)
Title |
---|
ZIQIANG MENG ET AL.: "Imaging of missile-borne bistatic forward-looking SAR", 《SIGNAL AND INFORMATION PROCESSING(CHINASIP), 2014 IEEE CHINA SUMMIT&INTERNATIONAL CONFERENCE ON》 * |
周松 等: "弹载SAR下降段成像算法研究", 《西安电子科技大学学报(自然科学版)》 * |
左绍山 等: "基于几何校正的聚束SAR快速分级后投影算法", 《电子与信息学报》 * |
邵鹏 等: "一种应用于斜视聚束模式的改进极坐标格式成像算法", 《电子与信息学报》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106353756A (zh) * | 2016-11-25 | 2017-01-25 | 西安电子科技大学 | 基于图像匹配的下降轨聚束合成孔径雷达定位方法 |
CN108594223A (zh) * | 2018-03-26 | 2018-09-28 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 机载sar图像目标定位方法 |
CN110244267A (zh) * | 2019-05-21 | 2019-09-17 | 成都航天科工微电子***研究院有限公司 | 弹载毫米波切向前视sar三维成像模型性能分析方法 |
CN110244267B (zh) * | 2019-05-21 | 2020-12-29 | 航天科工微电子***研究院有限公司 | 弹载毫米波切向前视sar三维成像模型性能分析方法 |
CN111060910A (zh) * | 2019-12-11 | 2020-04-24 | 西安电子科技大学 | 一种基于地形-影像匹配的InSAR载体反定位 |
CN111060910B (zh) * | 2019-12-11 | 2023-08-29 | 西安电子科技大学 | 一种基于地形-影像匹配的InSAR载体反定位 |
CN115166680A (zh) * | 2022-09-07 | 2022-10-11 | 中国科学院空天信息创新研究院 | 地面地物点的几何定位方法、装置、设备及介质 |
CN115166680B (zh) * | 2022-09-07 | 2022-11-29 | 中国科学院空天信息创新研究院 | 地面地物点的几何定位方法、装置、设备及介质 |
CN115638725A (zh) * | 2022-10-26 | 2023-01-24 | 成都清正公路工程试验检测有限公司 | 一种基于自动测量***的目标点位测量方法 |
CN115932823A (zh) * | 2023-01-09 | 2023-04-07 | 中国人民解放军国防科技大学 | 基于异源区域特征匹配的飞行器对地面目标定位方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105242252B (zh) | 2017-09-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105242252A (zh) | 基于图像匹配的下降轨聚束sar雷达定位方法 | |
WO2018223950A1 (zh) | 卫星激光测高数据辅助提取高程控制点方法 | |
CN102654576B (zh) | 基于sar图像和dem数据的图像配准方法 | |
US10670689B2 (en) | System and method for determining geo location of a target using a cone coordinate system | |
CN100520297C (zh) | 一种基于无偏差带的星敏感器地面校准方法 | |
CN105225241A (zh) | 无人机深度图像的获取方法及无人机 | |
KR100529401B1 (ko) | 합성 개구 레이더 영상을 이용한 수치표고자료 제작 장치및 그 방법 | |
CN103235306A (zh) | 一种适用于高速机动飞行器载sar成像的运动补偿方法 | |
CN103954953B (zh) | 一种基于数据驱动的机载激光雷达盲源误差补偿方法 | |
CN104180808A (zh) | 一种用于自主空中加油的圆形锥套视觉位姿解算方法 | |
IL226941A (en) | Automatic recovery method for unmanned aerial vehicles | |
CN103869311A (zh) | 实波束扫描雷达超分辨成像方法 | |
CN106643670B (zh) | 一种无人机航摄站点坐标求解装置及方法 | |
CN107765226A (zh) | 一种sar卫星雷达回波模拟方法、***和介质 | |
CN110646782A (zh) | 一种基于波形匹配的星载激光在轨指向检校方法 | |
CN102866393B (zh) | 一种基于pos与dem数据的sar多普勒参数估计方法 | |
CN103697883B (zh) | 一种基于天际线成像的飞行器水平姿态确定方法 | |
CN104567801B (zh) | 一种基于立体视觉的高精度激光测量方法 | |
CN106353756A (zh) | 基于图像匹配的下降轨聚束合成孔径雷达定位方法 | |
CN103197291A (zh) | 一种基于非停走模型的星载sar回波信号仿真方法 | |
CN109445310A (zh) | 一种红外成像制导飞行器半实物仿真***及方法 | |
CN104482924A (zh) | 旋成体目标位姿视觉测量方法 | |
CN107991676A (zh) | 星载单航过InSAR***对流层误差校正方法 | |
CN107607948A (zh) | 一种中轨道卫星sar的运动补偿及成像方法 | |
CN103218789A (zh) | 一种基于几何形变分步解算的自动化正射纠正方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |