CN104569975A - 微波凝视关联成像系统中时空随机辐射场的远场测量方法 - Google Patents

微波凝视关联成像系统中时空随机辐射场的远场测量方法 Download PDF

Info

Publication number
CN104569975A
CN104569975A CN201510067466.8A CN201510067466A CN104569975A CN 104569975 A CN104569975 A CN 104569975A CN 201510067466 A CN201510067466 A CN 201510067466A CN 104569975 A CN104569975 A CN 104569975A
Authority
CN
China
Prior art keywords
radiation field
measurement
radiation
field
measurement point
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
Application number
CN201510067466.8A
Other languages
English (en)
Other versions
CN104569975B (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.)
University of Science and Technology of China USTC
Original Assignee
University of Science and Technology of China USTC
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 University of Science and Technology of China USTC filed Critical University of Science and Technology of China USTC
Priority to CN201510067466.8A priority Critical patent/CN104569975B/zh
Publication of CN104569975A publication Critical patent/CN104569975A/zh
Application granted granted Critical
Publication of CN104569975B publication Critical patent/CN104569975B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种微波凝视关联成像系统中时空随机辐射场的远场测量方法,其包括:在成像区域中选取M个测量点以及一参考点,在相同的辐照条件下对第一个测量点处的辐射场与参考点处的辐射场进行同步连续的无扰测量及记录;在不改变辐射系统构型与激励信号的情况下,更换下一个测量点,并对该测量点处的辐射场与参考点处的辐射场进行同步连续的无扰测量及记录,直到遍历所有测量点;以参考点处记录的辐射场变化为时间维匹配基准,将所有测量点处的记录的辐射场在时间维对齐;提取不同时间片段测量点处的辐射场分布,获得辐射场随时间与空间两维变化的信息。该方法能够准确测量辐射场时间维和空间维的变化,可用于辐射场标定或者关联成像。

Description

微波凝视关联成像系统中时空随机辐射场的远场测量方法
技术领域
本发明涉及雷达技术领域,尤其涉及一种微波凝视关联成像系统中时空随机辐射场的远场测量方法。
背景技术
现有技术中的微波凝视成像,其通过构建具备时、空两维随机分布特征的辐射场,其与被探测目标相互作用后,将接收的散射回波与已知辐射场进行关联处理,获得目标反演像,该方法克服了传统实孔径雷达空间分辨率取决于天线孔径的缺陷,能实现高分辨的微波凝视成像;时空两维随机辐射场是上述微波凝视成像的关键,其随机特征表现为时空两维的非相关特性。
传统雷达成像的分析多从信号层面进行,而微波凝视关联成像完全从场的角度分析成像过程。凝视关联成像系统中各阵元天线辐射非相干场,以多相位中心的形式在远区形成时空两维随机辐射场,同时通过回波信号反演回波场的分布,将随机辐射场与回波场通过场与场关联的方法得到目标区域的图像。因此,微波凝视成像的方法需得到与回波场关联的辐射场,辐射场的准确性影响了成像的质量与分辨率。
传统单天线近场分布通过测量或者电磁计算方法得到,再根据惠更斯原理进行近远场外推,得到天线的远区场的幅相分布,单天线的远场分布通常使用方向图描述,主要关心远区功率的分布情况,而凝视成像系统关联处理时需要随机辐射雷达阵列远区辐射场的幅度及相位分布。
微波凝视关联成像系统中各阵元天线的激励信号具有较大带宽,且具有不同的调制方式,各天线的类型、形状及尺寸可不一致,极化方式也可能不相同,此外系统多天线辐射存在天线近场的耦合现象,以上因素导致微波凝视关联成像系统构造的远区辐射场具有时空两维随机的特性,远区场的获取存在较大难度。
发明内容
本发明的目的是提供一种微波凝视关联成像系统中时空随机辐射场的远场测量方 法,可以测量得到微波凝视关联成像中远区辐射场的时空变化,可用于辐射场标定或者关联成像。
本发明的目的是通过以下技术方案实现的:
一种微波凝视关联成像系统中时空随机辐射场的远场测量方法,该方法包括:
在成像区域中选取M个测量点以及一参考点,在相同的辐照条件下对第一个测量点处的辐射场与参考点处的辐射场进行同步连续的无扰测量及记录;
在不改变辐射系统构型与激励信号的情况下,更换下一个测量点,并对该测量点处的辐射场与参考点处的辐射场进行同步连续的无扰测量及记录,直到遍历所有测量点;
以参考点处记录的辐射场变化为时间维匹配基准,将所有测量点处的记录的辐射场在时间维对齐;
提取不同时间片段测量点处的辐射场分布,获得辐射场随时间与空间两维变化的信息。
进一步的,在相同的辐照条件下对第一个测量点处的辐射场与参考点处的辐射场进行同步连续的无扰测量及记录;在不改变辐射系统构型与激励信号的情况下,更换下一个测量点,并对该测量点处的辐射场与参考点处的辐射场进行同步连续的无扰测量及记录,直到遍历所有测量点的步骤包括:
搭建微波随机辐射系统,在随机辐射系统的辐照场景中放置一弱散射平板,板上铺吸波材料并划分网格;选取的M个测量点的位置为r1,...ri,...,rM
选取一参考点,位置为r0,同时在第一个测量点r1和参考点r0处放置一标准测量小天线接收直达的辐射场信号,并且输入接收机同步测量两点辐射场的变化并且记录,得到一组两点的辐射场变化测量数据;
在不改变辐射系统构型与激励信号的情况下,更换下一个测量点,并重复上述测量,直到遍历所有M个测量点;第i(i=1,...,M)次测量得到一组离散数据包括xi(n),xoi(n),n=1,...,N,其中xi(n)表示测量得到ri处的辐射场变化,xoi(n)表示对应的参考点处辐射场的变化,N表示一次测量的采样点数。
进一步的,所述以参考点处记录的辐射场变化为时间维匹配基准,将所有测量点处的记录的辐射场在时间维对齐包括:
将xo1(n)记为xo(n),计算关于k变量的函数取得最大值处的自变量值ki,此时xoi(n-ki)与xo(n)变化一致,参考信号时间维对齐,而 xi(n)是与xoi(n)同步测量得到,则xi(n-ki),i=1,...,M为各测量点处在时间维对齐后辐射场随时间的变化。
由上述本发明提供的技术方案可以看出,以连续测量记录辐射场时间维的变化,成像区域多点测量记录其空间维的变化,并且选择一参考点,记录参考点处场的变化作为时间维匹配的基准,用于各测量点辐射场时间维的对齐,通过上述方式克服了微波关联成像系统成像区域辐射场时空两维随机特性对测量带来的困难,能够准确测量辐射场时间维和空间维的变化。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域的普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他附图。
图1为本发明实施例微波关联成像的场景示意图;
图2为本发明实施例提供的一种微波凝视关联成像系统中时空随机辐射场的远场测量方法的流程图;
图3为本发明实施例测量场景示意图;
图4为本发明实施例测量辐射场变化示意图;
图5为本发明实施例各测量点辐射场的对齐;
图6为本发明实施例一次均匀致密测量时测量点的位置;
图7为本发明实施例获取不同时刻各测量点辐射场分布示意图。
具体实施方式
下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明的保护范围。
微波关联成像通过在成像区域构建具备时、空两维随机分布特征的辐射场,其与被探测目标相互作用后,根据接收到的散射回波信号反演回波场,并将回波场与已知辐射场进行关联处理,实现了场与场的关联处理,最终获得目标反演像。该方法克服了传统 实孔径雷达空间分辨率取决于天线孔径的缺陷,能实现高分辨的微波凝视成像;时空两维随机辐射场是上述微波凝视成像的关键,而如何得到准确的随机阵列远区辐射场为非常重要一环,直接关系到成像结果的可靠性。
图1给出了微波关联成像的场景示意图。随机辐射源阵列D对凝视成像区域S进行辐照,成像区域辐射场的分布函数为Einc(t,r),辐射场与目标作用后的散射场被接收天线接收,为Esca(t,r0′),将Einc(t,r)与Esca(t,r0′)进行关联处理,可得到成像区域目标的后向散射系数分布。随机辐射源阵列远区辐射场分布Einc(t,r)直接影响到成像结果的可靠性及质量,由于各阵元天线的激励信号具有较大带宽,且具有不同的调制方式,各天线的类型、形状及尺寸可不一致,极化方式也可能不相同,此外系统多天线辐射存在天线近场的耦合现象,导致了Einc(t,r)具有时空两维随机的特性,给测量带来了较大困难。
本发明实施例,可以克服微波关联成像系统成像区域辐射场时空两维随机特性对测量带来的困难,能够准确测量辐射场时间维和空间维的变化;具体如下:
实施例
图2为本发明实施例提供的一种微波凝视关联成像系统中时空随机辐射场的远场测量方法的流程图。如图2所示,该方法主要包括如下步骤:
步骤21、在成像区域中选取M个测量点以及一参考点,在相同的辐照条件下对第一个测量点处的辐射场与参考点处的辐射场进行同步连续的无扰测量及记录。
步骤22、在不改变辐射系统构型与激励信号的情况下,更换下一个测量点,并对该测量点处的辐射场与参考点处的辐射场进行同步连续的无扰测量及记录,直到遍历所有测量点。
具体来说,步骤21~步骤22可以通过如下方式来实现:
搭建微波随机辐射系统,在随机辐射系统的辐照场景中放置一弱散射平板,板上铺吸波材料并划分网格,所形成的格点有助于记录测量点位置;选取的M个测量点的位置为r1,...ri,...,rM
选取一参考点,位置为r0,同时在第一个测量点r1和参考点r0处放置一标准测量小天线接收直达的辐射场信号,并且输入接收机同步测量两点辐射场的变化并且记录,得到一组两点的辐射场变化测量数据;x1(n),xo1(n)n=1,...,N,N为一次测量中记录的点数,其中x1(n)表示r1处辐射场变化,xo1(n)表示与其对应的r0处辐射场变化。
图3给出了测量场景示意图,标准接收天线接收的直达辐射场信号经低噪声放大器后 接入高采样率存储深度大的示波器,观察并存储参考点及测量点处的接收波形,测量点处的示波器采样波形反应了各观察点处直达辐射场的变化规律。为最大限度地避免接收测量器件对辐射场分布的扰动,将它们置于平板背面(如图3(a)所示),同时,这大大缩短了天线到低噪声放大器(LNA)及LNA到高采样率示波器之间同轴线的长度,降低了信号的衰减,但标定验证场景受到实际条件的限制不能将接收测量器件置于平板背面时,需采用图3(b)所示的方法,采用吸波材料包裹的同轴线将天线接收到的信号引到较远的地方放大接收,该方法中,选取小尺寸LNA,将其紧邻天线放置而不是放在示波器处,是为了改善示波器采样信号的信噪比。高采样率大存储示波器直接对标准测量天线接收到的射频信号进行采样,保留了辐射场变化的细节,为后续的信号处理提供了更多的信息。
在不改变辐射系统构型与激励信号的情况下,更换下一个测量点,并重复上述测量,直到遍历所有M个测量点;第i(i=1,...,M)次测量得到一组离散数据包括xi(n),xoi(n),n=1,...,N,其中xi(n)表示测量得到ri处的辐射场变化,xoi(n)表示对应的参考点处辐射场的变化,N表示一次测量的采样点数。如图4为测量的波形示意图。
步骤23、以参考点处记录的辐射场变化为时间维匹配基准,将所有测量点处的记录的辐射场在时间维对齐。
图4中每组测量数据在相同的辐射系统构型和激励信号下获得,参考点处辐射场变化只有细微的差别,由辐照条件不可能完全一致及测量误差引起,可以将其作为时间维匹配基准;将所有测量点处的辐射场变化数据在时间维对齐的步骤如下:将xo1(n)记为xo(n),求关于k变量的函数取得最大值处的自变量值ki(其中xoi(n)的自变量取n-k时得到xoi(n-k),表示xoi(n-k)的共轭复数),此时xoi(n-ki)与xo(n)(即xo1(n))变化一致,参考信号时间维对齐,如图5(a)所示。而xi(n)与xoi(n)同步测量得到,则xi(n-ki),i=1,...,M为各测量点处在时间维对齐后辐射场随时间的变化,如图5(b)所示。
步骤24、提取不同时间片段测量点处的辐射场分布,获得辐射场随时间与空间两维变化的信息。
各测量点辐射场时间维对齐后,得到各测量点辐射场变化xi(n-ki),i=1,...,M,当n取相同的值时得到同一时刻各测量点上辐射场的分布。如图6为在一次均匀致密的测量条件下测量点的位置分布,测量点数目M=100,使用如上的方法测量,一次测量的采样点数为N=16000,各测量点处辐射场变化时间维对齐后如图7(a)所示,提取n=T0,T1两时刻的辐射场分布,分布图如7(b)所示。
本发明实施例中,以连续测量记录辐射场时间维的变化,成像区域多点测量记录其空间维的变化,并且选择一参考点,记录参考点处场的变化作为时间维匹配的基准,用于各测量点辐射场时间维的对齐,通过上述方式克服了微波关联成像系统成像区域辐射场时空两维随机特性对测量带来的困难,能够准确测量辐射场时间维和空间维的变化。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明披露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求书的保护范围为准。

Claims (3)

1.一种微波凝视关联成像系统中时空随机辐射场的远场测量方法,其特征在于,该方法包括:
在成像区域中选取M个测量点以及一参考点,在相同的辐照条件下对第一个测量点处的辐射场与参考点处的辐射场进行同步连续的无扰测量及记录;
在不改变辐射系统构型与激励信号的情况下,更换下一个测量点,并对该测量点处的辐射场与参考点处的辐射场进行同步连续的无扰测量及记录,直到遍历所有测量点;
以参考点处记录的辐射场变化为时间维匹配基准,将所有测量点处的记录的辐射场在时间维对齐;
提取不同时间片段测量点处的辐射场分布,获得辐射场随时间与空间两维变化的信息。
2.根据权利要求1所述的方法,其特征在于,在相同的辐照条件下对第一个测量点处的辐射场与参考点处的辐射场进行同步连续的无扰测量及记录;在不改变辐射系统构型与激励信号的情况下,更换下一个测量点,并对该测量点处的辐射场与参考点处的辐射场进行同步连续的无扰测量及记录,直到遍历所有测量点的步骤包括:
搭建微波随机辐射系统,在随机辐射系统的辐照场景中放置一弱散射平板,板上铺吸波材料并划分网格;选取的M个测量点的位置为r1,...ri,...,rM
选取一参考点,位置为r0,同时在第一个测量点r1和参考点r0处放置一标准测量小天线接收直达的辐射场信号,并且输入接收机同步测量两点辐射场的变化并且记录,得到一组两点的辐射场变化测量数据;
在不改变辐射系统构型与激励信号的情况下,更换下一个测量点,并重复上述测量,直到遍历所有M个测量点;第i(i=1,...,M)次测量得到一组离散数据包括xi(n),xoi(n),n=1,...,N,其中xi(n)表示测量得到ri处的辐射场变化,xoi(n)表示对应的参考点处辐射场的变化,N表示一次测量的采样点数。
3.根据权利要求1所述的方法,其特征在于,所述以参考点处记录的辐射场变化为时间维匹配基准,将所有测量点处的记录的辐射场在时间维对齐包括:
将xo1(n)记为xo(n),计算关于k变量的函数取得最大值处的自变量值ki,此时xoi(n-ki)与xo(n)变化一致,参考信号时间维对齐,而xi(n)是与xoi(n)同步测量得到,则xi(n-ki),i=1,...,M为各测量点处在时间维对齐后辐射场随时间的变化。2 -->
CN201510067466.8A 2015-02-09 2015-02-09 微波凝视关联成像系统中时空随机辐射场的远场测量方法 Expired - Fee Related CN104569975B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510067466.8A CN104569975B (zh) 2015-02-09 2015-02-09 微波凝视关联成像系统中时空随机辐射场的远场测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510067466.8A CN104569975B (zh) 2015-02-09 2015-02-09 微波凝视关联成像系统中时空随机辐射场的远场测量方法

Publications (2)

Publication Number Publication Date
CN104569975A true CN104569975A (zh) 2015-04-29
CN104569975B CN104569975B (zh) 2017-04-19

Family

ID=53086493

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510067466.8A Expired - Fee Related CN104569975B (zh) 2015-02-09 2015-02-09 微波凝视关联成像系统中时空随机辐射场的远场测量方法

Country Status (1)

Country Link
CN (1) CN104569975B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106980109A (zh) * 2017-05-03 2017-07-25 四川九洲电器集团有限责任公司 一种多子阵低空雷达
CN111257871A (zh) * 2020-03-09 2020-06-09 中国科学技术大学 可用于微波凝视关联成像的单天线辐射源设计方法
CN113835088A (zh) * 2021-09-24 2021-12-24 电子科技大学 一种自适应步进频积累的随机辐射雷达伪影抑制方法
CN117555197A (zh) * 2024-01-10 2024-02-13 湖南工商大学 一种基于弱光探测的扫描全景关联成像系统及其抑噪方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090237092A1 (en) * 2008-03-20 2009-09-24 The Curators Of The University Of Missouri Microwave and millimeter wave imaging system
CN102141617A (zh) * 2011-01-04 2011-08-03 中国科学技术大学 微波凝视成像关联方法
CN102141618A (zh) * 2011-01-04 2011-08-03 中国科学技术大学 微波凝视成像的方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090237092A1 (en) * 2008-03-20 2009-09-24 The Curators Of The University Of Missouri Microwave and millimeter wave imaging system
CN102141617A (zh) * 2011-01-04 2011-08-03 中国科学技术大学 微波凝视成像关联方法
CN102141618A (zh) * 2011-01-04 2011-08-03 中国科学技术大学 微波凝视成像的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
YUANYUE GUO ET AL.: "Research on Sensing Matrix Characteristics in Microwave Staring Correlated Imaging Based on Compressed Sensing", 《2014 IEEE INTERNATIONAL CONFERENCE ON IMAGING SYSTEMS AND TECHNIQUES(IST) PROCEEDINGS》 *
孙慧峰 等: "SAR天线相位中心的分析及测量", 《SAR天线相位中心的分析及测量 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106980109A (zh) * 2017-05-03 2017-07-25 四川九洲电器集团有限责任公司 一种多子阵低空雷达
CN106980109B (zh) * 2017-05-03 2020-11-03 四川九洲防控科技有限责任公司 一种多子阵低空雷达
CN111257871A (zh) * 2020-03-09 2020-06-09 中国科学技术大学 可用于微波凝视关联成像的单天线辐射源设计方法
CN111257871B (zh) * 2020-03-09 2023-06-16 中国科学技术大学 可用于微波凝视关联成像的单天线辐射源设计方法
CN113835088A (zh) * 2021-09-24 2021-12-24 电子科技大学 一种自适应步进频积累的随机辐射雷达伪影抑制方法
CN117555197A (zh) * 2024-01-10 2024-02-13 湖南工商大学 一种基于弱光探测的扫描全景关联成像系统及其抑噪方法
CN117555197B (zh) * 2024-01-10 2024-03-26 湖南工商大学 一种基于弱光探测的扫描全景关联成像系统及其抑噪方法

Also Published As

Publication number Publication date
CN104569975B (zh) 2017-04-19

Similar Documents

Publication Publication Date Title
Appel-Hansen Reflectivity level of radio anechoic chambers
CN112098961B (zh) 一种多极化sar载荷内定标模式设计及数据分析方法
US12117371B2 (en) Method and apparatus of full-field vibration measurement via microwave sensing
Zhang et al. A generalized approach for SAR and MIMO radar imaging of building interior targets with compressive sensing
Wang et al. High-resolution imaging using a wideband MIMO radar system with two distributed arrays
CN105572487B (zh) 一种阵列天线近场幅相测量方法及幅相测量器
WO2021143306A1 (zh) 一种基于时间反转和压缩感知的电磁非视线成像方法
US10101489B2 (en) System for exploring underground geophysical properties and method for analyzing underground geophysical properties using the same
US9400325B2 (en) Method and apparatus for increasing angular resolution in an automotive radar system
CN104569975A (zh) 微波凝视关联成像系统中时空随机辐射场的远场测量方法
CN104833863B (zh) 一种高频相控阵天线的远场暗室测试系统及方法
CN107003388A (zh) 用于增加汽车雷达系统中的角分辨率的方法和设备
KR101777381B1 (ko) Pcl 수신기에서 적응 필터를 이용한 표적반사 신호의 도래각 추정 장치 및 방법
CN104991249A (zh) 一种山体滑坡的mimo雷达监测系统及监测方法
CN106556783A (zh) 一种变电站内基于特高频相控阵原理的局部放电测向方法
CN111413666A (zh) 一种阵列测频测向联合接收机设计方法
US8570207B1 (en) Method, technique, and system for detecting Brillouin precursors at microwave frequencies for enhanced performance in various applications
KR101564729B1 (ko) 위상배열 레이더의 배열면 정렬 방법
KR100962419B1 (ko) 다배열 지하레이더를 이용한 3차원 시뮬레이션 시스템 및그 지하정보 영상구현방법
Ma et al. System error analysis and calibration methods for multi-channel SAR
Guo et al. Super-resolution staring imaging radar based on stochastic radiation fields
CN114325700A (zh) 一种星载多通道sar动目标成像方法
Li et al. Time reversal imaging based on synchronism
CN111708023B (zh) 毫米波多发多收线阵的相位不一致性校正方法及成像系统
CN110267862A (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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170419