CN107290745A - 一种电离层探测方法及系统 - Google Patents

一种电离层探测方法及系统 Download PDF

Info

Publication number
CN107290745A
CN107290745A CN201710383277.0A CN201710383277A CN107290745A CN 107290745 A CN107290745 A CN 107290745A CN 201710383277 A CN201710383277 A CN 201710383277A CN 107290745 A CN107290745 A CN 107290745A
Authority
CN
China
Prior art keywords
electromagnetic wave
ionospheric
usrp
ionosphere
frequency
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
Application number
CN201710383277.0A
Other languages
English (en)
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.)
Nanchang University
Original Assignee
Nanchang University
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 Nanchang University filed Critical Nanchang University
Priority to CN201710383277.0A priority Critical patent/CN107290745A/zh
Publication of CN107290745A publication Critical patent/CN107290745A/zh
Pending legal-status Critical Current

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
    • 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/882Radar or analogous systems specially adapted for specific applications for altimeters
    • 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
    • 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
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/12Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electromagnetic waves

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种电离层探测方法及系统,其中,方法包括:使用USRP发射设定频率的电磁波和接收电离层反射的电磁波,并记录发射和接收的电磁波信息;USRP通过千兆网口与计算机通信,计算机根据发射和接收的电磁波信息计算设定电磁波发射频率下的电离层高度、电离层电子密度、电离层电子径向速度;改变USRP发射的电磁波频率,获取不同电磁波发射频率下的电离层高度、电离层电子密度及电离层电子径向速度并绘制电离层分布图。有益效果:USRP通过千兆网口与计算机通信,USRP实现电磁波的收发及初步的信号处理,计算机进行较为复杂的信号处理和相关探测目标的计算,利用USRP和已编程的计算机即可实现电离层探测,成本低廉,灵活、便利、开放程度高,易于维护。

Description

一种电离层探测方法及系统
技术领域
本发明涉及雷达探测领域,尤其是涉及一种电离层探测方法及系统。
背景技术
电离层探测雷达是对电离层进行探测的重要工具,传统电离层探测雷达中对信号的处理都是依赖硬件实现的,虽然可以精确地测量出电离层的相关数据,但是其价格一般很高,硬件升级和维护很不方便,需要变更探测需求时,往往因为探测硬件固化而不够灵活。利用传统的电离层探测雷达研究、探测电离层,成本太高,且不够灵活、不够开放、不够方便。
发明内容
本发明的目的在于克服上述技术不足,提出一种电离层探测方法及系统,解决现有技术中的上述技术问题。
为达到上述技术目的,本发明的技术方案提供一种电离层探测方法,包括:
S1、使用USRP发射设定频率的电磁波,并记录发射的电磁波的信息为第一电磁波信息;
S2、使用USRP接收经电离层反射的电磁波,获取接收的电磁波的信息为第二电磁波信息;
S3、USRP通过千兆网口与计算机通信,计算机根据第一电磁波信息、第二电磁波信息计算设定电磁波发射频率下的电离层高度、电离层电子密度、电离层电子径向速度;
S4、通过计算机向USRP发送指令,改变USRP发射的电磁波频率,再次执行步骤S1-S3,获取不同电磁波发射频率下的电离层高度、电离层电子密度及电离层电子径向速度并绘制电离层分布图。
本发明还提供一种电离层探测系统,包括:
发射模块:使用USRP发射设定频率的电磁波,并记录发射的电磁波的信息为第一电磁波信息;
接收模块:使用USRP接收经电离层反射的电磁波,获取接收的电磁波的信息为第二电磁波信息;
计算模块:USRP通过千兆网口与计算机通信,计算机根据第一电磁波信息、第二电磁波信息计算设定电磁波发射频率下的电离层高度、电离层电子密度、电离层电子径向速度;
循环计算与绘图模块:通过计算机向USRP发送指令,改变USRP发射的电磁波频率,再次执行发射模块、接收模块、计算模块的操作,获取不同电磁波发射频率下的电离层高度、电离层电子密度及电离层电子径向速度并绘制电离层分布图。
与现有技术相比,本发明的有益效果包括:USRP通过千兆网口与计算机通信,USRP实现电磁波的收发及初步的信号处理,计算机进行较为复杂的信号处理和电离层高度、电离层电子密度、电离层电子径向速度的计算,利用USRP和已编程的计算机即可实现电离层探测,成本低廉,灵活、便利、开放程度高,易于维护。
附图说明
图1是本发明提供的一种电离层探测方法流程图;
图2是本发明提供的一种电离层探测系统结构框图。
附图中:1、电离层探测系统,11、发射模块,12、接收模块,13、计算模块,14、循环计算与绘图模块。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
USRP全称是通用软件无线电外围设备,实现了射频信号发射和接收。它由一个母板和不同的子板组成。本发明的实施例中,母板包含4个12bit/64M模/数转换器(ADC)、4个14bit/128M数/模转换器(DAC)、一个可重复编程的FPGA芯片。每个USRP母板上可安装四个子板,其中两个用于接收,另外两个用于发射。通过连接不同的子板,USRP可以处理不同频段的信号。
本发明的实施例中,将USRP和计算机通过千兆网口连接,USRP通过千兆网口与计算机通信,USRP中的FPGA进行信号的上下变频、抽样和内插等数字信号简单操作,USRP发射或者接收的电磁波信号经过FPGA处理后传递给计算机,计算机进行FFT处理、调制解调等较为复杂的数字信号处理,且计算机负责计算需要探测的电离层高度、电离层电子密度、电离层电子径向速度,并绘制电离层分布图。
本发明的实施例中,与USRP连接的计算机上的软件模块编程通过GNU Radio平台实现,通过对软件模块的编程实现对USRP的控制、计算机上的信号处理以及相关的计算任务。GNU是一个由FSF(Free Software Foundation)提供支持的开放源代码著名项目。GNURadio运行于Linux系统,其软件代码完全公开,其设计的核心思想是:基带信号处理过程由通用处理器而不是DSP(DSP芯片即指能够实现数字信号处理技术的芯片)实现;基带信号通过千兆网口与射频前端USRP相连,由USRP上的FPGA实现数字上下变频(DUC/DDC),然后通过ADC/DAC芯片与射频子板连接;覆盖不同频段的射频子板将模拟信号变换为射频进行信号收与发。GNU Radio提供了一个信号处理模块的库,并且有个机制可以把单个的处理模块连接在一起形成一个系统。编程者通过建立一个流向图(flow graph)就可以快速搭建一个无线电通信系统。GNU Radio的编程实现主要基于Python脚本语言和C++语言。C++语言被用于编写各种信号处理模块,具有较高的执行效率。GNU Radio提供了超过100个现成的信号处理模块,并且编写新的处理模块也是非常容易。而软件图形化接口和信号处理模块的链接机制则是通过Python脚本语言实现的。
本发明提供了一种电离层探测方法,包括:
S1、使用USRP发射设定频率的电磁波,并记录发射的电磁波的信息为第一电磁波信息;
S2、使用USRP接收经电离层反射的电磁波,获取接收的电磁波的信息为第二电磁波信息;
S3、USRP通过千兆网口与计算机通信,计算机根据第一电磁波信息、第二电磁波信息计算设定电磁波发射频率下的电离层高度、电离层电子密度、电离层电子径向速度;
S4、通过计算机向USRP发送指令,改变USRP发射的电磁波频率,再次执行步骤S1-S3,获取不同电磁波发射频率下的电离层高度、电离层电子密度及电离层电子径向速度并绘制电离层分布图。
本发明所述的电离层探测方法,步骤S1中:
第一电磁波信息包括:发射电磁波的时间t1、发射的电磁波的频率f1
本发明所述的电离层探测方法,步骤S2中:
第二电磁波信息包括:接收到经电离层反射的电磁波的时间t2、接收的经电离层反射的电磁波的频率f2
本发明所述的电离层探测方法,步骤S3中:
根据计算电离层高度,其中,R表示电离层高度,C表示电磁波传播速度(光速),t0表示发射电磁波与接收到反射的电磁波之间的时间间隔,t0=t2-t1
根据计算电离层电子密度,其中,Ne表示电离层电子密度;
根据计算电离层电子径向速度,其中,Vr表示电离层电子径向速度,fd表示多普勒频移,fd=f1-f2
本发明所述的电离层探测方法,步骤S4中:
由于不同频率电磁波对于电离层的穿透作用不同,所以不同发射频率电磁波在电离层的不同高度被反射,由公式计算出的电离层高度是该频率电磁波的电离层反射高度,改变电磁波的发射频率,则可以测量该区域电离层的区域范围,并且可将电磁波发射频率和电磁波发射频率对应的电离层高度输入MATLAB中,以电磁波发射频率、电离层高度分别为x、y轴绘制电离层分布图。
本发明还提供一种电离层探测系统1,包括:
发射模块11:使用USRP发射设定频率的电磁波,并记录发射的电磁波的信息为第一电磁波信息;
接收模块12:使用USRP接收经电离层反射的电磁波,获取接收的电磁波的信息为第二电磁波信息;
计算模块13:USRP通过千兆网口与计算机通信,计算机根据第一电磁波信息、第二电磁波信息计算设定电磁波发射频率下的电离层高度、电离层电子密度、电离层电子径向速度;
循环计算与绘图模块14:通过计算机向USRP发送指令,改变USRP发射的电磁波频率,再次执行发射模块11、接收模块12、计算模块的操作13,获取不同电磁波发射频率下的电离层高度、电离层电子密度及电离层电子径向速度并绘制电离层分布图。
本发明所述的电离层探测系统1,发射模块11中:
第一电磁波信息包括:发射电磁波的时间t1、发射的电磁波的频率f1
本发明所述的电离层探测系统1,接收模块12中:
第二电磁波信息包括:接收到经电离层反射的电磁波的时间t2、接收的经电离层反射的电磁波的频率f2
本发明所述的电离层探测系统1,计算模块13中:
根据计算电离层高度,其中,R表示电离层高度,C表示电磁波传播速度,t0表示发射电磁波与接收到反射的电磁波之间的时间间隔,t0=t2-t1
根据计算电离层电子密度,其中,Ne表示电离层电子密度;
根据计算电离层电子径向速度,其中,Vr表示电离层电子径向速度,fd表示多普勒频移,fd=f1-f2
本发明所述的电离层探测系统1,循环计算与绘图模块14中:
获取不同电磁波发射频率下的电离层高度后,将电磁波发射频率和电磁波发射频率对应的电离层高度输入MATLAB中,从而绘制电离层分布图。
与现有技术相比,本发明的有益效果包括:USRP通过千兆网口与计算机通信,USRP实现电磁波的收发及初步的信号处理,计算机进行较为复杂的信号处理和电离层高度、电离层电子密度、电离层电子径向速度的计算,利用USRP和已编程的计算机即可实现电离层探测,成本低廉,灵活、便利、开放程度高,易于维护。
以上所述本发明的具体实施方式,并不构成对本发明保护范围的限定。任何根据本发明的技术构思所做出的各种其他相应的改变与变形,均应包含在本发明权利要求的保护范围内。

Claims (10)

1.一种电离层探测方法,其特征在于,包括步骤:
S1、使用USRP发射设定频率的电磁波,并记录发射的电磁波的信息为第一电磁波信息;
S2、使用所述USRP接收经电离层反射的电磁波,获取接收的电磁波的信息为第二电磁波信息;
S3、所述USRP通过千兆网口与计算机通信,所述计算机根据所述第一电磁波信息、所述第二电磁波信息计算设定电磁波发射频率下的电离层高度、电离层电子密度、电离层电子径向速度;
S4、通过所述计算机向所述USRP发送指令,改变所述USRP发射的电磁波频率,再次执行步骤S1-S3,获取不同电磁波发射频率下的所述电离层高度、所述电离层电子密度及所述电离层电子径向速度并绘制电离层分布图。
2.如权利要求1所述的电离层探测方法,其特征在于,步骤S1中:
所述第一电磁波信息包括:发射电磁波的时间t1、发射的电磁波的频率f1
3.如权利要求2所述的电离层探测方法,其特征在于,步骤S2中:
所述第二电磁波信息包括:接收到经电离层反射的电磁波的时间t2、接收的经电离层反射的电磁波的频率f2
4.如权利要求3所述的电离层探测方法,其特征在于,步骤S3中:
根据计算所述电离层高度,其中,R表示电离层高度,C表示电磁波传播速度,t0表示发射电磁波与接收到反射的电磁波之间的时间间隔,t0=t2-t1
根据计算所述电离层电子密度,其中,Ne表示电离层电子密度;
根据计算所述电离层电子径向速度,其中,Vr表示电离层电子径向速度,fd表示多普勒频移,fd=f1-f2
5.如权利要求1所述的电离层探测方法,其特征在于,步骤S4中:
获取不同电磁波发射频率下的所述电离层高度后,将电磁波发射频率和电磁波发射频率对应的所述电离层高度输入MATLAB中,从而绘制电离层分布图。
6.一种电离层探测系统,其特征在于,包括:
发射模块:使用所述USRP发射设定频率的电磁波,并记录发射的电磁波的信息为所述第一电磁波信息;
接收模块:使用所述USRP接收经电离层反射的电磁波,获取接收的电磁波的信息为所述第二电磁波信息;
计算模块:USRP通过千兆网口与计算机通信,计算机根据所述第一电磁波信息、所述第二电磁波信息计算设定电磁波发射频率下的所述电离层高度、所述电离层电子密度、所述电离层电子径向速度;
循环计算与绘图模块:通过所述计算机向所述USRP发送指令,改变所述USRP发射的电磁波频率,再次执行发射模块、接收模块、计算模块的操作,获取不同电磁波发射频率下的所述电离层高度、所述电离层电子密度及所述电离层电子径向速度并绘制电离层分布图。
7.如权利要求6所述的电离层探测系统,其特征在于,发射模块中:
所述第一电磁波信息包括:发射电磁波的时间t1、发射的电磁波的频率f1
8.如权利要求7所述的电离层探测系统,其特征在于,接收模块中:
所述第二电磁波信息包括:接收到经电离层反射的电磁波的时间t2、接收的经电离层反射的电磁波的频率f2
9.如权利要求8所述的电离层探测系统,其特征在于,计算模块中:
根据计算所述电离层高度,其中,R表示电离层高度,C表示电磁波传播速度,t0表示发射电磁波与接收到反射的电磁波之间的时间间隔,t0=t2-t1
根据计算所述电离层电子密度,其中,Ne表示电离层电子密度;
根据计算所述电离层电子径向速度,其中,Vr表示电离层电子径向速度,fd表示多普勒频移,fd=f1-f2
10.如权利要求6所述的电离层探测系统,其特征在于,循环计算与绘图模块中:
获取不同电磁波发射频率下的所述电离层高度后,将电磁波发射频率和电磁波发射频率对应的所述电离层高度输入MATLAB中,从而绘制所述电离层分布图。
CN201710383277.0A 2017-05-26 2017-05-26 一种电离层探测方法及系统 Pending CN107290745A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710383277.0A CN107290745A (zh) 2017-05-26 2017-05-26 一种电离层探测方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710383277.0A CN107290745A (zh) 2017-05-26 2017-05-26 一种电离层探测方法及系统

Publications (1)

Publication Number Publication Date
CN107290745A true CN107290745A (zh) 2017-10-24

Family

ID=60095200

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710383277.0A Pending CN107290745A (zh) 2017-05-26 2017-05-26 一种电离层探测方法及系统

Country Status (1)

Country Link
CN (1) CN107290745A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110988884A (zh) * 2019-12-30 2020-04-10 陇东学院 一种基于高频地波雷达的中纬度电离层探测方法
CN113376631A (zh) * 2021-08-11 2021-09-10 武汉大学 一种电离层Es层的超分辨率垂直探测方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1998038532A2 (en) * 1997-02-26 1998-09-03 Nathan Blaunstein Method of earthquake prediction
CN106134521B (zh) * 2011-05-19 2014-08-27 中国电子科技集团公司第二十二研究所 一种小型化电离层垂直探测装置
CN104898101A (zh) * 2015-05-11 2015-09-09 南昌大学 电离层非相干散射雷达探测威力仿真方法
KR101605385B1 (ko) * 2016-01-06 2016-03-22 (주)에스이랩 전리권의 전자 밀도 분포 예측 장치 및 그 방법

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1998038532A2 (en) * 1997-02-26 1998-09-03 Nathan Blaunstein Method of earthquake prediction
CN106134521B (zh) * 2011-05-19 2014-08-27 中国电子科技集团公司第二十二研究所 一种小型化电离层垂直探测装置
CN104898101A (zh) * 2015-05-11 2015-09-09 南昌大学 电离层非相干散射雷达探测威力仿真方法
KR101605385B1 (ko) * 2016-01-06 2016-03-22 (주)에스이랩 전리권의 전자 밀도 분포 예측 장치 및 그 방법

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
汪正: ""软件雷达电离层探测技术研究"", 《中国优秀硕士学位论文全文数据库信息科技辑》 *
白波 等: ""一种新型星载无线电等离子体成像仪"", 《遥测遥控》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110988884A (zh) * 2019-12-30 2020-04-10 陇东学院 一种基于高频地波雷达的中纬度电离层探测方法
CN113376631A (zh) * 2021-08-11 2021-09-10 武汉大学 一种电离层Es层的超分辨率垂直探测方法

Similar Documents

Publication Publication Date Title
CN105676198B (zh) 一种用于脉冲式雷达测试的回波脉冲延迟产生装置
CN201765324U (zh) 一种相控阵天线自动测试仪
CN111693978B (zh) 基于mimo毫米波雷达的散点检测方法
CN102520388A (zh) 一种利用相控阵原理并结合rssi测距的定位装置
CN105005040B (zh) 一种雷达测角的方法
CN104714209A (zh) 一种基于uwb与激光测距组合的动态定位方法及装置
CN107607797A (zh) 基于无人机的天线性能测量方法及装置
CN110441732B (zh) 一种陌生散射环境中的单站多目标定位方法
CN104391183A (zh) 近场测量天线远区场特性快速计算方法
CN101776746A (zh) 信号到达角度定位与信号波程差测量装置和方法
CN105333887A (zh) 无线电高度表信号模拟系统
CN105426602B (zh) 一种基于多核dsp实现脉内调制特性分析的装置及方法
CN109001723A (zh) 雷达流速仪
CN107290745A (zh) 一种电离层探测方法及系统
RU2003120811A (ru) Способ и система для радиолокационного измерения скоростей и координат объектов (варианты)
CN108490412A (zh) 一种双基地雷达海杂波测量试验装置
CN112014812A (zh) 一种相控阵风廓线雷达标定系统及方法
CN105911556B (zh) 一种带接收模拟器的多波束测深系统
CN109788550A (zh) 一种基于超宽带技术的室内田径场定位系统及方法
CN206132864U (zh) 一种超高频雷达天线阵列方向图测量装置
CN111175734B (zh) 一种基于虚拟等效合成孔径的辐射源定位测速方法及系统
CN108594172A (zh) 一种射频信号测距的方法、存储介质、系统及中央处理器
CN104749391A (zh) 一种基于三轴加速度传感器的速度测量系统及方法
CN109828276B (zh) 一种基于干涉算法的风廓线雷达
CN107707315B (zh) 一种基于IQxel的无线产品边缘信号量测方法

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20171024