CN103489176B - 一种对于严重几何畸变的sar图像进行同名点提取的方法 - Google Patents

一种对于严重几何畸变的sar图像进行同名点提取的方法 Download PDF

Info

Publication number
CN103489176B
CN103489176B CN201210194637.XA CN201210194637A CN103489176B CN 103489176 B CN103489176 B CN 103489176B CN 201210194637 A CN201210194637 A CN 201210194637A CN 103489176 B CN103489176 B CN 103489176B
Authority
CN
China
Prior art keywords
image
same place
obtains
coherence
sar
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
Application number
CN201210194637.XA
Other languages
English (en)
Other versions
CN103489176A (zh
Inventor
李银伟
向茂生
毛永飞
韦立登
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Electronics of CAS
Original Assignee
Institute of Electronics of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Institute of Electronics of CAS filed Critical Institute of Electronics of CAS
Priority to CN201210194637.XA priority Critical patent/CN103489176B/zh
Publication of CN103489176A publication Critical patent/CN103489176A/zh
Application granted granted Critical
Publication of CN103489176B publication Critical patent/CN103489176B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Image Processing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种对于严重几何畸变的SAR图像进行同名点提取的方法。该方法包括:确定同航向相邻航带的照射区域位置信息;根据所述位置信息计算所述SAR图像的重叠位置,并裁剪得到图像I,即整个成像场景的远距端图像,和图像II,即整个成像场景的近距端图像;根据图像I和图像II的构像几何信息及斜距关系,将所述图像II仿射变换到图像I的斜距坐标系下得到图像IV;对图像I和图像IV进行候选同名点的提取;以相干系数图为导引对候选同名点进行筛选剔除;将得到的同名点在图像IV上的坐标逆映射至图像II上,得到同名点坐标。本发明方法能够提取由于重叠区域较小导致几何畸变严重不一致的同航向相邻条带间的同名点。

Description

一种对于严重几何畸变的SAR图像进行同名点提取的方法
技术领域
本发明涉及信号处理技术领域,尤其是一种对于严重几何畸变的SAR图像进行同名点提取的方法。
背景技术
合成孔径雷达(SyntheticApertureRadar,SAR)是一种主动式二维高分辨成像雷达,其在距离向(波束照射方向)通过发射大时间带宽积的线性调频信号,采用脉冲压缩技术来获取高分辨率,在方位向(平台运动方向,通常与距离向垂直)利用目标和雷达的相对运动形成的轨迹来构成一个合成孔径来取代庞大的阵列实孔径来获取高分辨率。干涉SAR(InterferometricSAR,InSAR)是指利用两部或多部不同位置处的SAR对同一场景进行观测,并通过数据的后处理获取场景的高程信息。由InSAR系统的回波数据获取场景数学高程模型(DigitalElevationModel,DEM)的主要步骤包括运动补偿、成像、配准、滤波、相位解缠、干涉定标、高程反演。
InSAR数据处理中,干涉参数(基线长度、基线角和干涉相位偏置)直接影响最终获取DEM的精度。这需要在各景中实地布放足量且分布合理的地面控制点(GroundControlPoints,GCPs)来对它们进行标定。而在大区域地形测绘时,机载InSAR系统通常采取分若干条带获取数据,每条带分若干场景进行处理的方式来作业。这使得在作业过程中实地布放足量且分布合理的GCPs工作量大、作业效率低;而且野外采集GCPs受地形条件限制,存在某些测区如荒山、沼泽等难以实现GCPs布放。因此需要充分利用相邻影像之间具有同一地理位置特征的同名点(TiePoints,TPs),并联合少量GCPs,对测绘区域内所有景数据的干涉参数进行联合标定。在这过程中,TPs的自动提取显得尤为重要。
SAR的侧视成像模式和相干成像机理使得自动提取TPs较为困难。首先,相邻场景重叠区域的成像视角差异,造成重叠区影像存在灰度不一致、旋转变形、尺度缩放以及仿射变形;其次,同一分辨单元内各散射点随机散射信号相互叠加产生的相干斑噪声,降低了图像的质量,掩盖了图像的细节结构,很大程度上影响SAR影像的特征提取。
TPs的提取分为两类:同航带相邻场景和同航向相邻航带。其中,同航带相邻场景属于同一飞行航迹,成像几何如图1(a)所示。相邻场景的重叠区域是同一数据经过成像处理获得,只是处理时的理想运动轨迹不一致造成其存在较小的尺度缩放、角度缩放及仿射变形等,因此同名点提取相对来说较为容易。
同航向相邻航带属于两次飞行航迹,成像几何如图1(b)所示。重叠区位于航带I的远距端,航带Ⅱ的近距端。SAR侧视成像几何和斜距投影使得重叠区在航带Ⅱ的地距采样间隔大于航带I,例如重叠度50%,近、远距和中心视角分别为30°,60°和45°,此时重叠测区在两航带的地距分辨单元满足关系,进而造成相邻影像重叠区的几何畸变差异不一致(如图1(c)所示)。图1(c)中,SAR影像的近距压缩效应使得重叠测区等地距间隔摆放的9个圆点在两幅影像的几何畸变(距离向)不一致;当两航带的重叠范围非常小(小于15%)时,该局部几何畸变将更加显著,这使得难以在原始图像上直接提取同名点。
发明内容
本发明的目的是解决现有技术中存在的问题,提供一种对于严重几何畸变的SAR图像进行同名点提取的方法。该方法能够提取由于重叠区域较小导致几何畸变严重不一致的同航向相邻条带问的同名点。
本发明所提出的一种对于严重几何畸变的SAR图像进行同名点提取的方法,其特征在于,该方法包括以下步骤:
步骤S1,根据IMU/GPS记录的载机航迹信息及飞行前测量的天线初始安装姿态信息,确定同航向相邻航带的照射区域位置信息;
步骤S2,对于重叠区域较小的两幅同航向相邻航带SAR图像,根据所述位置信息计算具有重叠区域的所述同航向相邻航带SAR图像的重叠位置,并对所述同航向相邻航带SAR图像进行裁剪,得到对应的图像I,即整个成像场景的远距端图像,和图像II,即整个成像场景的近距端图像;
步骤S3,根据图像I和图像II的构像几何信息及斜距关系,将所述图像II仿射变换到图像I的斜距坐标系下得到图像IV;
步骤S4,对图像I和图像IV进行候选同名点的提取;
步骤S5,以相干系数图为导引对候选同名点进行筛选剔除;
步骤S6,将经筛选剔除后得到的同名点在图像IV上的坐标根据步骤S3的映射关系逆映射至图像II上,得到重叠区域图像I和图像II的同名点坐标。
本发明方法能够提取由于重叠区域较小(小于15%)导致几何畸变严重不一致的同航向相邻条带间的同名点。
附图说明
图1是相邻场景成像的几何关系图;
图2是本发明方法的流程图;
图3是重叠区域较小的同航向相邻航带SAR图像;
图4是图像II变换后的图像III和图像IV;
图5是同航向相邻航带SAR的斜距-地距剖面图;
图6是在图像I和图像IV上提取的同名点示意图;
图7是同名点在图像I和图像III(II旋转后)上的示意图;
图8是同名点在原始图像(I和II)上的示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。
图2是本发明方法的流程图,如图2所示,本发明所提出的一种对于严重几何畸变的SAR图像进行同名点提取的方法包括以下步骤:
步骤S1,根据IMU(InertialMeasurementUnit,惯性测量装置)/GPS记录的载机航迹信息及飞行前测量的天线初始安装姿态信息,确定同航向相邻航带的照射区域位置信息;
所述位置信息比如可以是照射区域的四个顶点的经纬度坐标。
步骤S2,对于重叠区域较小的两幅同航向相邻航带SAR图像(如图3所示),根据所述位置信息计算具有重叠区域的所述同航向相邻航带SAR图像的重叠位置,并将所述重叠位置在所述同航向相邻航带SAR图像中对应的图像区域裁剪下来,得到对应的图像I,即整个成像场景的远距端图像,和图像II,即整个成像场景的近距端图像(如图3所示);
步骤S3,根据图像I和图像II的构像几何信息及斜距关系,将所述图像II仿射变换到图像I的斜距坐标系下得到图像IV;
所述步骤S3进一步包括以下步骤:
步骤S31,根据对同航向相邻航带原始回波数据在成像处理时进行运动补偿后得到的理想运动轨迹的航向角、斜视角等构像几何信息,计算图像II的理想运动轨迹相对于图像I的旋转角度θ;
理论上,SAR原理要求天线相位中心(AntennaPhaseCenter,APC)做匀速直线运动;实际中,由于受大气湍流和人为控制的影响,APC总是偏离理想运动状态而产生运动误差。这需要在成像处理时根据IMU/GPS记录的载机姿态信息进行理想运动轨迹的拟合,然后通过计算APC偏离理想运动轨迹的大小对回波数据进行运动补偿。
步骤S32,根据所述旋转角度θ对图像II进行旋转变换,得到与图像I的理想运动轨迹相互平行的图像III(如图4(a)所示);
步骤S33,根据所述两幅同航向相邻航带SAR图像的重叠区域的雷达平台高度、最近斜距、采样间隔等成像信息,建立图像I和图像III之间对应的斜距关系;
同航向相邻航带SAR的斜距-地距剖面图如图5所示,其中,A1、A2分别表示获取图像I和图像III时的APC位置,h1、h2分别表示A1、A2的高度,R1、R2分别表示APC位置A1、A2与重叠区域左边缘P0的最近斜距,dR1表示图像I的斜距采样间隔。对于图像I上与点P0在斜距向上相隔n点的点Pn,可以通过式(1)计算点P0和点Pn的地距距离
| P 0 P n | = ( R 1 + n · dR 1 ) 2 - h 1 2 - R 1 2 - h 1 2 - - - ( 1 )
然后由式(2)可以计算在图像III中点Pn与点P0的斜距差ΔR2
Δ R 2 = ( R 2 2 - h 2 2 + | P 0 P n | ) 2 + h 2 2 - R 2 - - - ( 2 )
最后根据图像III的斜距采样间隔dR2,可以获得Pn在图像III中与点P0相隔的点数,即Pn在图像III中相对于点P0的斜距向位置ΔR2/dR2,从而得到两图像之间对应的斜距关系。
步骤S34,根据步骤S33得到的斜距关系,将旋转后的图像III重采样到图像I的斜距坐标系下得到图像IV(如图4(b)所示)。
步骤S4,采用高斯滤波、SIFT(ScaleInvariantFeatureTransform,尺度不变特征变换)或SURF(Speeded-UpRobustFeature,快速鲁棒特征)方法、Euclidean相似性度量距离、BBF算法、NN/SN准则、双向匹配及RANSAC鲁棒算法对图像I和图像IV进行“候选TPs”的提取;
所述步骤S4进一步包括以下步骤:
步骤S41,对图像I和图像IV进行高斯滤波,以降低图像中的斑点噪声;
步骤S42,利用SIFT或SURF方法对高斯滤波后的两图像检测多尺度特征并建立特征描述符,得到图像I的特征集合A={aj}(j=1,2,…)和图像IV的特征集合B={bi}(i=1,2,…);
步骤S43,利用两图像多尺度特征之间的Euclidean相似性度量距离,使用BBF(BestBinFirst,最优节点优先)算法,按照NN/SN准则(如果特征点与其最近邻和次近邻的距离之比小于指定阈值,则判定该特征点和其最近邻为匹配点)对两图像的特征点进行快速的双向匹配:先在图像IV的特征集合B={bi}(i=1,2,…)中寻找与图像I特征集合A={aj}(j=1,2,…)中的特征点相匹配的点集,然后对于B中已匹配的特征点,在A中寻找与其匹配的点,从而建立两图像特征点的初始匹配对集合;
步骤S44,利用RANSAC(RandomSampleConsensus,随机抽样一致性)鲁棒算法估计两幅图像的单应矩阵H,初始匹配对集合中满足|bi-Haj|<ε的内点称为“候选TPs”。其中,ε为距离阈值,取值越小说明两个兴趣点的匹配精度越高。
步骤S5,以相干系数图为导引对“候选TPs”进行筛选剔除,剔除位于相干性低的阴影、水域及相位噪声大的人工目标邻域的匹配点对;
在干涉定标中需要利用TPs的相位信息,而相干系数图是相位数据质量好坏的评价标准。采用最大似然相干估算器分别计算每幅图像的相干系数,得到每幅图像的相干系数图,若候选TPs在图像I中的相干系数γ1在图像II中的相干系数γ2及其相应邻域内的相干系数均值都分别满足:γi>0.95,(i=1,2),则满足上述条件的候选TP为所需要的TP(如图6所示)。
步骤S6,将经筛选剔除后得到的TPs在图像IV上的坐标根据步骤S3的映射关系逆映射至图像II上(如图7和8所示),得到图像I和图像II重叠区域的TPs坐标。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (8)

1.一种对于严重几何畸变的SAR图像进行同名点提取的方法,其特征在于,该方法包括以下步骤:
步骤S1,根据IMU/GPS记录的载机航迹信息及飞行前测量的天线初始安装姿态信息,确定同航向相邻航带的照射区域位置信息;
步骤S2,对于重叠区域较小的两幅同航向相邻航带SAR图像,根据所述位置信息计算具有重叠区域的所述同航向相邻航带SAR图像的重叠位置,并对所述同航向相邻航带SAR图像进行裁剪,得到对应的图像I,即整个成像场景的远距端图像,和图像II,即整个成像场景的近距端图像;
步骤S3,根据图像I和图像II的构像几何信息及斜距关系,将所述图像II仿射变换到图像I的斜距坐标系下得到图像IV;
步骤S4,对图像I和图像IV进行候选同名点的提取;
步骤S5,以相干系数图为导引对候选同名点进行筛选剔除;
步骤S6,将经筛选剔除后得到的同名点在图像IV上的坐标根据步骤S3的映射关系逆映射至图像II上,得到重叠区域图像I和图像II的同名点坐标;
所述步骤S4进一步包括以下步骤:
步骤S41,对图像I和图像IV进行高斯滤波,以降低图像中的斑点噪声;
步骤S42,利用SIFT或SURF方法对高斯滤波后的两图像检测多尺度特征并建立特征描述符,得到图像I的特征集合A={aj},j=1,2,…,和图像IV的特征集合B={bi},i=1,2,…;
步骤S43,利用两图像多尺度特征之间的Euclidean相似性度量距离,使用最优节点优先BBF算法,按照NN/SN准则对两图像的特征点进行快速的双向匹配,从而建立两图像特征点的初始匹配对集合;
步骤S44,利用随机抽样一致性RANSAC鲁棒算法估计两幅图像的单应矩阵H,初始匹配对集合中满足|bi-Haj|<ε的内点称为候选同名点,其中,ε为距离阈值。
2.根据权利要求1所述的方法,其特征在于,所述步骤S3进一步包括以下步骤:
步骤S31,根据对同航向相邻航带原始回波数据在成像处理时进行运动补偿后得到的理想运动轨迹的航向角、斜视角构像几何信息,计算图像II的理想运动轨迹相对于图像I的旋转角度θ;
步骤S32,根据所述旋转角度θ对图像II进行旋转变换,得到与图像I的理想运动轨迹相互平行的图像III;
步骤S33,根据重叠区域的成像信息,建立图像I和图像III之间对应的斜距关系;
步骤S34,根据步骤S33得到的斜距关系,将旋转后的图像III重采样到图像I的斜距坐标系下得到图像IV。
3.根据权利要求2所述的方法,其特征在于,所述步骤S33中重叠区域的成像信息包括雷达平台高度、最近斜距和采样间隔。
4.根据权利要求1所述的方法,其特征在于,所述NN/SN准则为:如果特征点与其最近邻和次近邻的距离之比小于指定阈值,则判定该特征点和其最近邻为匹配点。
5.根据权利要求1所述的方法,其特征在于,所述步骤S43中的双向匹配进一步为:先在图像IV的特征集合B={bi},i=1,2,…,中寻找与图像I特征集合A={aj},j=1,2,…,中的特征点相匹配的点集,然后对于B中已匹配的特征点,在A中寻找与其匹配的点。
6.根据权利要求1所述的方法,其特征在于,所述步骤S5进一步包括采用最大似然相干估算器分别计算每幅图像的相干系数,得到每幅图像的相干系数图。
7.根据权利要求6所述的方法,其特征在于,所述步骤S5中,若候选同名点在图像I中的相干系数γ1、在图像II中的相干系数γ2及其相应邻域内的相干系数均值都分别满足一特定条件,则该候选同名点为所需要的同名点。
8.根据权利要求7所述的方法,其特征在于,所述特定条件为:γi>0.95,i=1,2。
CN201210194637.XA 2012-06-13 2012-06-13 一种对于严重几何畸变的sar图像进行同名点提取的方法 Active CN103489176B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210194637.XA CN103489176B (zh) 2012-06-13 2012-06-13 一种对于严重几何畸变的sar图像进行同名点提取的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210194637.XA CN103489176B (zh) 2012-06-13 2012-06-13 一种对于严重几何畸变的sar图像进行同名点提取的方法

Publications (2)

Publication Number Publication Date
CN103489176A CN103489176A (zh) 2014-01-01
CN103489176B true CN103489176B (zh) 2016-02-03

Family

ID=49829374

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210194637.XA Active CN103489176B (zh) 2012-06-13 2012-06-13 一种对于严重几何畸变的sar图像进行同名点提取的方法

Country Status (1)

Country Link
CN (1) CN103489176B (zh)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104978742B (zh) * 2015-06-12 2017-12-05 北京邮电大学 基于级联结构的图像配准方法及装置
CN105842691B (zh) * 2016-03-17 2018-01-12 中国测绘科学研究院 大区域InSAR形变结果的自动化平差及镶嵌方法及装置
CN106529587B (zh) * 2016-08-26 2019-06-25 北京航空航天大学 基于目标点识别的视觉航向识别方法
CN106778587B (zh) * 2016-12-09 2019-09-27 自然资源部四川测绘产品质量监督检验站(四川省测绘产品质量监督检验站) 基于无人机影像的飞行状态检测方法及装置
CN106844733B (zh) * 2017-02-13 2020-04-03 哈尔滨理工大学 基于词汇树信息融合与豪斯多夫距离结合的图像检索方法
CN108335262A (zh) * 2017-12-16 2018-07-27 中煤航测遥感集团有限公司 一种基于物方反算的dem拼接方法及dsm拼接方法
CN108304883B (zh) * 2018-02-12 2022-09-16 西安电子科技大学 基于改进sift的sar图像匹配方法
CN108896009A (zh) * 2018-04-29 2018-11-27 天津大学 一种基于有效相干性的大型线状结构体沉降监测方法
CN109188433B (zh) * 2018-08-20 2022-11-04 南京理工大学 基于无控制点的双机载sar图像目标定位的方法
CN109493274B (zh) * 2018-10-18 2020-12-11 华南农业大学 一种基于坐标转换的图像扭曲处理方法
CN109903352B (zh) * 2018-12-24 2021-03-30 中国科学院遥感与数字地球研究所 一种卫星遥感影像大区域无缝正射影像制作方法
WO2021051439A1 (zh) * 2019-09-26 2021-03-25 深圳市速腾聚创科技有限公司 校正方法、装置、存储介质及多通道激光雷达
CN111968035B (zh) * 2020-08-05 2023-06-20 成都圭目机器人有限公司 基于损失函数的图像相对旋转角计算方法
CN112379377B (zh) * 2020-10-30 2022-08-12 上海卫星工程研究所 分布式InSAR卫星长条带测绘优化SAR任务规划方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1426021A (zh) * 2002-12-19 2003-06-25 上海交通大学 遥感图像的非线性配准方法
WO2005008281A1 (en) * 2003-07-19 2005-01-27 Gamma Remote Sensing Research And Consulting Ag Method to improve interferometric signatures by coherent point scatterers
CN1987896A (zh) * 2005-12-23 2007-06-27 中国科学院中国遥感卫星地面站 高分辨率sar影像配准处理方法及系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1426021A (zh) * 2002-12-19 2003-06-25 上海交通大学 遥感图像的非线性配准方法
WO2005008281A1 (en) * 2003-07-19 2005-01-27 Gamma Remote Sensing Research And Consulting Ag Method to improve interferometric signatures by coherent point scatterers
CN1987896A (zh) * 2005-12-23 2007-06-27 中国科学院中国遥感卫星地面站 高分辨率sar影像配准处理方法及系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
MULTISENSOR MAPPING USING SAR IN CONJUNCTION WITH OPTICAL DATA;J. Raggam 等;《isprs.org》;19931231;556-563 *
基于大尺度双边SIFT的SAR图像同名点自动提取方法;王山虎 等;《电子与信息学报》;20120229;第34卷(第2期);287-293 *
面向InSAR稀疏控制点测图的同名点提取方法;姜丽敏 等;《电子与信息学报》;20111231;第33卷(第12期);2837-2845 *

Also Published As

Publication number Publication date
CN103489176A (zh) 2014-01-01

Similar Documents

Publication Publication Date Title
CN103489176B (zh) 一种对于严重几何畸变的sar图像进行同名点提取的方法
CN113359097B (zh) 一种毫米波雷达和相机联合标定的方法
CN110244302B (zh) 地基合成孔径雷达影像像元坐标三维变换方法
CN103487809B (zh) 一种基于BP算法和时变基线的机载InSAR数据处理方法
KR20130004227A (ko) Sar 영상 내의 픽셀의 지리 좌표를 결정하기 위한 방법
CN102866393B (zh) 一种基于pos与dem数据的sar多普勒参数估计方法
CN113093184B (zh) 一种基于视频合成孔径雷达的干涉测量方法
CN101620272B (zh) 一种逆合成孔径雷达的目标转速估计方法
CN104007439A (zh) 一种干涉圆迹sar高程估计处理方法
CN103308031B (zh) 一种基于卫星三线阵ccd影像的云顶高度反演方法
WO2020151213A1 (zh) 一种空地结合的潮间带一体化测绘方法
CN112882030B (zh) InSAR成像干涉一体化处理方法
CN103472450A (zh) 基于压缩感知的非均匀空间构形分布式sar动目标三维成像方法
CN109633639B (zh) Topsar干涉数据的高精度快速配准方法
CN108983231B (zh) 一种基于视频合成孔径雷达的干涉视频测量方法
CN104361563A (zh) 基于gps的高光谱遥感图像几何精校正方法
CN110703252B (zh) 干涉合成孔径雷达阴影区域数字高程模型修正方法
CN103454633A (zh) 一种基于后向投影算法的干涉sar动基线处理方法
Feng et al. A hierarchical network densification approach for reconstruction of historical ice velocity fields in East Antarctica
KR101873732B1 (ko) 위성레이더 오프셋트래킹 지표변위 지도의 왜곡 보정 장치 및 그 방법
CN102798380B (zh) 线阵图像中目标运动参数的测量方法
CN109946682B (zh) 基于ICESat/GLAS的GF3数据基线估计方法
CN105403886A (zh) 一种机载sar定标器图像位置自动提取方法
CN116202410A (zh) 地质灾害监测方法、装置、电子设备和存储介质
Jaud et al. Method for orthorectification of terrestrial radar maps

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant