CN104266762A - 一种基于环境辐射改变的现场目标发射率测量***及方法 - Google Patents

一种基于环境辐射改变的现场目标发射率测量***及方法 Download PDF

Info

Publication number
CN104266762A
CN104266762A CN201410484662.0A CN201410484662A CN104266762A CN 104266762 A CN104266762 A CN 104266762A CN 201410484662 A CN201410484662 A CN 201410484662A CN 104266762 A CN104266762 A CN 104266762A
Authority
CN
China
Prior art keywords
radiation
target
infrared eye
emissivity
spectroscope
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
CN201410484662.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.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and Technology
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 Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN201410484662.0A priority Critical patent/CN104266762A/zh
Publication of CN104266762A publication Critical patent/CN104266762A/zh
Pending legal-status Critical Current

Links

Landscapes

  • Radiation Pyrometers (AREA)

Abstract

本发明公开了一种基于环境辐射改变的现场目标发射率测量***及方法。该***包括沿光路方向依次共光轴设置的目标、分光镜和红外探测器;该***还包括黑体和聚光镜,所述黑体的辐射光经聚光镜汇聚于分光镜上,并由分光镜反射到达目标,黑体的辐射光通过与目标作用然后反射至分光镜并到达探测器。首先对红外探测器进行辐射定标;采用黑体辐射为基准建立***响应函数;测量在环境辐射强度为Vs1的条件下,红外探测器的输出响应值V1;改变环境辐射强度为Vs2的条件下,测量红外探测器的输出响应值V2;根据测得的两种不同的辐射强度条件下红外探测器的输出响应值,算出目标的反射率并求得目标的发射率。本发明为红外辐射非接触精确测温提供了必要参量。

Description

一种基于环境辐射改变的现场目标发射率测量***及方法
技术领域
本发明属于现场目标发射率的测量技术领域,特别是一种基于环境辐射改变的现场目标发射率测量***及方法。
背景技术
目标表面真实温度的辐射测量在许多科学研究中是一项重要的任务,尤其是表面真实温度的非接触精确测量更为困难和关键。红外成像测温技术具有非接触(不会破坏被测物体的温度场)、直观、灵敏度高(可分辨0.01℃的温度差)、快速(几毫秒内测出目标温度)、测温范围广(从-170℃到3200℃以上)、检测距离可近可远、可实现夜视、安全等优点,不但在高压电线巡检、工业生产等民用领域得到广泛的应用;而且在侦察与制导、伪装设计与检测等军用领域中也得到了广泛应用。红外探测器测温主要受被测物体表面发射率的影响,但反射率、环境温度、大气温度、测量距离和大气衰减等因素的影响也不容忽视。这些影响因素会导致热像仪的测温不准,影响其应用领域,尤其是对物体表面发射率估计的不准确,更影响温度测量的精确性。发射率指物体通过表面向外辐射的电磁能与同温度的黑体在相同条件下所辐射的电磁能的比值,是衡量物体辐射能力强弱的参数,其数值在0与1之间变化,表示实际物体的热辐射与黑体辐射的接近程度,根据黑体辐射定律,只要知道材料的发射率,就知道任何物体的红外辐射特性。理想黑体的发射率为1,实际物体的发射率小于1。因为发射率是辐射波长、温度、方向以及表面状态的函数,并在很大程度上取决于物体的表面状态。因此在现场测量目标的红外辐射特性中,手册中的数据就不可靠了。这就需要在现场对发射率进行实际测量。本文提出一种现场测量常温物体发射率的方法,进行了理论分析和实验测量,给出了实验结果。
目前,测量目标发射率主要采用直接测量法。直接测量法是将目标与黑体置于同一条件同一温度下,目标输出值与黑体输出值之比即为目标发射率。直接测量法需要在无大气、无环境影响和已知物体本身温度的情况下进行,只适合于实验室测量。这种测量方法对于野外现场的目标发射率测量毫无意义。
发明内容
本发明的目的在于提供一种适应范围广泛的基于环境辐射改变的现场目标发射率测量***及方法,能够无接触、精确地测量目标的发射率。
实现本发明目的的技术解决方案为:一种基于环境辐射改变的现场目标发射率测量***,包括沿光路方向依次共光轴设置的目标、分光镜和红外探测器;该***还包括黑体和聚光镜,所述黑体的辐射光经聚光镜汇聚于分光镜上,并由分光镜反射到达目标,黑体的辐射光通过与目标作用然后反射至分光镜并到达探测器。
一种基于环境辐射改变的现场目标发射率测量方法,步骤如下:
步骤1,在利用红外探测器测试之前,先对红外探测器进行辐射定标;
步骤2,采用黑体辐射为基准建立***响应函数;
步骤3,根据步骤2建立的***响应函数,测量在环境辐射强度为Vs1的条件下,红外探测器的输出响应值V1
步骤4,根据步骤2建立的***响应函数,改变环境辐射强度为Vs2的条件下,测量红外探测器的输出响应值V2
步骤5,根据步骤3和步骤4测得的两种不同的辐射强度条件下红外探测器的输出响应值,算出目标的反射率,再根据目标反射率和发射率的关系求得目标的发射率。
本发明与现有技术相比,其显著优点为:(1)非接触测量:采用红外探测器作为探测器来测量目标和环境的辐射特性,不需要接触被测目标,因此可以实现非接触测量;(2)测量精确度高:采用非接触式测量方法,所以测量过程中不会改变目标和环境的温度场,因此也就不会改变目标和环境的辐射特性,这样测量结果的精确度非常高,误差最高为1.7%;(3)适用范围广:测量目标发射率的方法不仅适用于实验室条件下的测量,而且在野外条件下测量时也能得到比较精确的结果。
附图说明
图1是本发明基于环境辐射改变的现场目标发射率测量***的结构示意图。
图2是黑体辐射亮度与输出值之间的关系。
图3是黑体辐射到瓷板的图像。
图4是瓷板上黑体辐射部分的图像。
具体实施方式
下面结合附图及具体实施例对本发明作进一步详细说明。
结合图1,本发明基于环境辐射改变的现场目标发射率测量***,包括沿光路方向依次共光轴设置的目标、分光镜和红外探测器,所述分光镜倾斜放置且分光镜与光路方向之间的夹角为135°;该***还包括黑体和聚光镜,所述黑体的辐射光经聚光镜汇聚于分光镜上,并由分光镜反射到达目标,黑体的辐射光通过与目标作用然后反射至分光镜并到达探测器,聚光镜的作用是为了使到达目标的辐射强烈而均匀。
本发明基于环境辐射改变的现场目标发射率测量方法,步骤如下:
步骤1,在利用红外探测器测试之前,先对红外探测器进行辐射定标;红外探测器是以黑体辐射源为基准,其理论基础是黑体辐射理论,通过与黑体辐射源的辐射量进行对比测量,因此在利用红外探测器测试之前,先要进行定标。
物体在波段λ1~λ2内的辐射亮度,由普朗克公式可推得:
L = ϵ π ∫ λ 1 λ 2 C 1 λ - 5 ( e C 2 λT - 1 ) dλ - - - ( 1 )
式中,ε为物体发射率,C1为第一辐射常量,C1=2πhc2=3.7415×108W·μm4/m2,C2为第二辐射常量,C2=hc/k=1.43879×104μm·K,λ为波长,T为开氏温度;
红外探测器的标定模型为:
V=αL+Lf  (2)
式中,V为红外探测器的输出值,Lf为红外探测器本身杂散能量,α为红外探测器的亮度响应度。
实验室中利用黑体在不同温度下测量,给出一系列的输出值,画出目标亮度与输出值的关系曲线,通过最小二乘法对输出值进行线性拟合,得出红外探测器的辐射亮度响应关系。针对Y=kX+b形式的最小二乘法拟合公式如下:
k = ( xy ‾ - x ‾ × y ‾ ) ÷ ( x 2 - x ‾ 2 ) - - - ( 3 )
b = y ‾ - k x ‾ - - - ( 4 )
其中, x ‾ = 1 n Σ x i ; xy ‾ = 1 n Σ x i y i .
根据实验测得的黑体温度对应的灰度和辐射亮度数据表1,可以绘制出定标曲线如图2所示,该曲线表示的是黑体辐射亮度与输出值之间的关系。
表1.黑体温度对应的灰度和辐射亮度数据
通过实验测量由最小二乘法拟合的公式(3)和(4)得到该红外探测器的辐射亮度响应关系为:
y=20.9672x-981.3774  (5)
步骤2,建立***响应函数,采用黑体辐射为基准建立***响应函数,***响应函数是反映测量所得输出信号与输入辐射之间函数关系的,它由光学***、电子线路和探测器的响应率等因素决定。
在实际测量时,红外探测器接收到的有效辐射包括三个部分:目标自身辐射、目标反射周围环境辐射和大气辐射,表示如下式(6)所示:
L=L0+ρLs+Lp  (6)
其中,L为到达探测器表面的辐射亮度,L0为目标自身辐射亮度,ρ为目标反射率,Ls为环境辐射,Lp为大气辐射;
设大气透过率为1,将式(6)代入式(2)中得:
V=α(L0+ρLs+Lp)+Lf  (7)
对式(7)进行处理可得:
V=(αL0+Lf)+ρ(αLs+Lf)-ρLf+αLp  (8)
观察式(8),其中令
V0=αL0+Lf  (9)
Vs=αLs+Lf  (10)
其中,V0为只有目标的辐射作用到探测器上时红外探测器输出值,Vs为只有环境辐射作用到探测器上时红外探测器的输出值。
由上两式可将式(8)整理为:
V=V0+ρVs-ρLf+αLp  (11)
步骤3,根据步骤2建立的***响应函数,根据式(11),测量在环境辐射强度为Vs1的条件下,红外探测器的输出响应值V1
V1=V0+ρVs1-ρLf+αLp  (12)
步骤4,根据步骤2建立的***响应函数,改变环境辐射强度为Vs2的条件下,红外探测器的输出响应值V2。根据式(11)可通过改变环境辐射Vs,在同一传输路径上大气辐射没有改变,得到另一组输出:
V2=V0+ρVs2-ρLf+αLp  (13)
步骤5,根据步骤3和步骤4测得的两种不同的辐射强度条件下红外探测器的输出响应值,算出目标的反射率,再根据目标反射率和发射率的关系求得目标的发射率。
将将步骤4所得红外探测器的输出响应值V2和步骤3所得红外探测器的输出响应值V1相减,得:
V2-V1=ρ(Vs2-Vs1)  (14)
整理得
ρ = V 2 - V 1 V s 2 - V s 1 - - - ( 15 )
由基尔霍夫定律,当辐射能入射到物体表面时,包括三个过程:吸收、反射、透射。对于不透明表面,只包括吸收、反射两个过程,其发射率等于热吸收率:
ρ+ε=1  (16)
则根据式(15)和式(16)可得
ϵ = 1 - ρ = 1 - V 2 - V 1 V s 2 - V s 1 - - - ( 17 )
其中,在环境辐射强度为Vs1的条件下,红外探测器的输出响应值为V1;改变环境辐射强度为Vs2的条件下,测量红外探测器的输出响应值为V2。通过改变黑体温度达到改变环境辐射的目的。
实施例1
选取陶瓷、白纸和塑料板三种材料作为测量目标,将目标与黑体置于同一条件同一温度下,黑体和目标与探测器的距离0.5m,每种材料重复测量5次,取平均值。试验测试结果如表2所示。
表2材料发射率测量值
为了更能准确的评价测量结果的不确定程度,分别计算三组测量值的标准不确定度,利用公式(18),它反映了数值相对于平均值的离散程度。不确定度小于2%,由此可以说明该方法测量的稳定性较好。
S = Σ ( x i - x ‾ ) 2 N - 1 - - - ( 18 )
图3和图4为瓷板测试图,图中较亮部分为瓷板表面反射黑体所辐射的(8~14um)红外波段,提取探测器的输出值。改变黑体温度,获取不同环境辐射下另一输出值,并通过式(12)和式(13)计算两次环境辐射值,最后通过式(17)计算目标发射率。表3是三组材料经过多次测量得到的发射率值。
表3目标发射率测量数据
由表3可知,本文方法得到的物体的发射率实验结果与直接测量法数值的最大误差为0.017(1.7%)。

Claims (6)

1.一种基于环境辐射改变的现场目标发射率测量***,其特征在于,包括沿光路方向依次共光轴设置的目标、分光镜和红外探测器;该***还包括黑体和聚光镜,所述黑体的辐射光经聚光镜汇聚于分光镜上,并由分光镜反射到达目标,黑体的辐射光通过与目标作用然后反射至分光镜并到达探测器。
2.根据权利要求1所述的基于环境辐射改变的现场目标发射率测量***,其特征在于,所述分光镜倾斜放置且分光镜与光路方向之间的夹角为135°。
3.一种基于环境辐射改变的现场目标发射率测量方法,其特征在于,步骤如下:
步骤1,在利用红外探测器测试之前,先对红外探测器进行辐射定标;
步骤2,采用黑体辐射为基准建立***响应函数;
步骤3,根据步骤2建立的***响应函数,测量在环境辐射强度为Vs1的条件下,红外探测器的输出响应值V1
步骤4,根据步骤2建立的***响应函数,改变环境辐射强度为Vs2的条件下,测量红外探测器的输出响应值V2
步骤5,根据步骤3和步骤4测得的两种不同的辐射强度条件下红外探测器的输出响应值,算出目标的反射率,再根据目标反射率和发射率的关系求得目标的发射率。
4.根据权利要求3所述的基于环境辐射改变的现场目标发射率测量方法,其特征在于,步骤1所述的辐射标定,具体如下:
物体在波段λ1~λ2内的辐射亮度L,由普朗克公式推得:
L = ϵ π ∫ λ 1 λ 2 C 1 λ - 5 ( e C 2 λT - 1 ) dλ
式中,ε为物体发射率,C1为第一辐射常量,C1=2πhc2=3.7415×108W·μm4/m2,C2为第二辐射常量,C2=hc/k=1.43879×104μm·K,λ为波长,T为开氏温度;
红外探测器的标定模型为:
V=αL+Lf
式中,V为红外探测器的输出值,Lf为红外探测器本身杂散能量,α为红外探测器的亮度响应度。
5.根据权利要求3所述的基于环境辐射改变的现场目标发射率测量方法,其特征在于,步骤2所述采用黑体辐射为基准建立***响应函数,具体如下:
在实际测量时,红外探测器接收到的有效辐射包括三个部分:目标自身辐射、目标反射周围环境辐射和大气辐射,如下式所示:
L=L0+ρLs+Lp
其中,L为到达探测器表面的辐射亮度,L0为目标自身辐射亮度,ρ为目标反射率,Ls为环境辐射,Lp为大气辐射;
设大气透过率为1,将L代入红外探测器的标定模型得:
V=α(L0+ρLs+Lp)+Lf
对上式进行处理得:
V=(αL0+Lf)+ρ(αLs+Lf)-ρLf+αLp
令上式中,
V0=αL0+Lf
Vs=αLs+Lf
因此得:
V=V0+ρVs-ρLf+αLp
其中,V0为只有目标的辐射作用到探测器上时红外探测器输出值,Vs为只有环境辐射作用到探测器上时红外探测器的输出值。
6.根据权利要求3所述的基于环境辐射改变的现场目标发射率测量方法,其特征在于,步骤5所述根据目标反射率和发射率的关系求得目标的发射率,具体步骤如下:
将步骤4所得红外探测器的输出响应值V2和步骤3所得红外探测器的输出响应值V1相减,得:
V2-V1=ρ(Vs2-VS1)
整理得:
ρ = V 2 - V 1 V s 2 - V s 1
由基尔霍夫定律,当辐射能入射到物体表面时,对于不透明表面,只包括吸收、反射两个过程,其发射率等于热吸收率,即:
ρ+ε=1
因此:
ϵ = 1 - ρ = 1 - V 2 - V 1 V s 2 - V s 1
其中,在环境辐射强度为Vs1的条件下,红外探测器的输出响应值为V1;改变环境辐射强度为Vs2的条件下,测量红外探测器的输出响应值为V2
CN201410484662.0A 2014-09-19 2014-09-19 一种基于环境辐射改变的现场目标发射率测量***及方法 Pending CN104266762A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410484662.0A CN104266762A (zh) 2014-09-19 2014-09-19 一种基于环境辐射改变的现场目标发射率测量***及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410484662.0A CN104266762A (zh) 2014-09-19 2014-09-19 一种基于环境辐射改变的现场目标发射率测量***及方法

Publications (1)

Publication Number Publication Date
CN104266762A true CN104266762A (zh) 2015-01-07

Family

ID=52158303

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410484662.0A Pending CN104266762A (zh) 2014-09-19 2014-09-19 一种基于环境辐射改变的现场目标发射率测量***及方法

Country Status (1)

Country Link
CN (1) CN104266762A (zh)

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104634765A (zh) * 2015-03-05 2015-05-20 姜志富 基于光学辐射测量仪测量大气透过率的装置和方法
CN104833429A (zh) * 2015-03-27 2015-08-12 中国计量科学研究院 一种基于控制背景辐射的黑体发射率测量装置及方法
CN105628215A (zh) * 2016-02-17 2016-06-01 无锡元创华芯微机电有限公司 一种红外探测器单黑体响应率测试方法
CN105784631A (zh) * 2015-12-10 2016-07-20 中国航空工业集团公司北京长城计量测试技术研究所 一种负温及常温不透明材料发射率测量装置
CN105784333A (zh) * 2016-03-29 2016-07-20 渤海大学 红外窗口的高温光谱发射率测量方法及装置
CN107478548A (zh) * 2017-08-16 2017-12-15 清华大学 一种确定运动中高温颗粒表观发射率的装置和方法
CN108760770A (zh) * 2018-05-24 2018-11-06 上海航天测控通信研究所 一种热控涂层损耗测试装置及方法
CN108896187A (zh) * 2018-05-15 2018-11-27 重庆邮电大学 一种回转窑表面红外双波段测温方法
CN111579488A (zh) * 2020-06-16 2020-08-25 烟台艾睿光电科技有限公司 物体表面发射率的测量方法、装置、设备及存储介质
CN112504471A (zh) * 2020-11-25 2021-03-16 天津津航技术物理研究所 一种应用于智能监控***的实时红外测温方法
CN113324663A (zh) * 2021-07-09 2021-08-31 中国人民解放军32801部队 一种外场目标发射率的测量方法及***
CN114235690A (zh) * 2021-11-25 2022-03-25 中国人民解放军空军工程大学 航空器涂层的面红外发射率测量方法及装置
CN114441045A (zh) * 2022-01-28 2022-05-06 电子科技大学 一种准确测量辐射温度的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
邓明德等: "黑体辐射公式的积分解及应用", 《遥感信息》 *

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104634765B (zh) * 2015-03-05 2017-06-23 姜志富 基于光学辐射测量仪测量大气透过率的装置和方法
CN104634765A (zh) * 2015-03-05 2015-05-20 姜志富 基于光学辐射测量仪测量大气透过率的装置和方法
CN104833429A (zh) * 2015-03-27 2015-08-12 中国计量科学研究院 一种基于控制背景辐射的黑体发射率测量装置及方法
CN105784631A (zh) * 2015-12-10 2016-07-20 中国航空工业集团公司北京长城计量测试技术研究所 一种负温及常温不透明材料发射率测量装置
CN105784631B (zh) * 2015-12-10 2018-10-30 中国航空工业集团公司北京长城计量测试技术研究所 一种负温及常温不透明材料发射率测量装置
CN105628215B (zh) * 2016-02-17 2019-02-19 无锡元创华芯微机电有限公司 一种红外探测器单黑体响应率测试方法
CN105628215A (zh) * 2016-02-17 2016-06-01 无锡元创华芯微机电有限公司 一种红外探测器单黑体响应率测试方法
CN105784333A (zh) * 2016-03-29 2016-07-20 渤海大学 红外窗口的高温光谱发射率测量方法及装置
CN107478548B (zh) * 2017-08-16 2019-10-25 清华大学 一种确定运动中高温颗粒表观发射率的装置和方法
CN107478548A (zh) * 2017-08-16 2017-12-15 清华大学 一种确定运动中高温颗粒表观发射率的装置和方法
CN108896187A (zh) * 2018-05-15 2018-11-27 重庆邮电大学 一种回转窑表面红外双波段测温方法
CN108760770B (zh) * 2018-05-24 2021-05-11 上海航天测控通信研究所 一种热控涂层损耗测试装置及方法
CN108760770A (zh) * 2018-05-24 2018-11-06 上海航天测控通信研究所 一种热控涂层损耗测试装置及方法
CN111579488A (zh) * 2020-06-16 2020-08-25 烟台艾睿光电科技有限公司 物体表面发射率的测量方法、装置、设备及存储介质
CN112504471A (zh) * 2020-11-25 2021-03-16 天津津航技术物理研究所 一种应用于智能监控***的实时红外测温方法
CN113324663A (zh) * 2021-07-09 2021-08-31 中国人民解放军32801部队 一种外场目标发射率的测量方法及***
CN113324663B (zh) * 2021-07-09 2022-02-25 中国人民解放军32801部队 一种外场目标发射率的测量方法及***
CN114235690A (zh) * 2021-11-25 2022-03-25 中国人民解放军空军工程大学 航空器涂层的面红外发射率测量方法及装置
CN114235690B (zh) * 2021-11-25 2023-11-21 中国人民解放军空军工程大学 航空器涂层的面红外发射率测量方法及装置
CN114441045A (zh) * 2022-01-28 2022-05-06 电子科技大学 一种准确测量辐射温度的方法
CN114441045B (zh) * 2022-01-28 2023-07-04 电子科技大学 一种准确测量辐射温度的方法

Similar Documents

Publication Publication Date Title
CN104266762A (zh) 一种基于环境辐射改变的现场目标发射率测量***及方法
Zhang et al. A method for reducing the influence of measuring distance on infrared thermal imager temperature measurement accuracy
CN102661799B (zh) 故障定位方法与***
CN103076101B (zh) 一种红外热像仪像元点的标定方法
RU2523775C2 (ru) Способ и система коррекции на основе квантовой теории для повышения точности радиационного термометра
CN109655161B (zh) 基于红外热像仪的目标红外积分辐射强度测试方法及设备
CN104155008A (zh) 一种红外温度监测***测量误差修正方法
CN108072459A (zh) 一种测量钢板温度场及计算其辐射强度的方法
CN103604504A (zh) 一种红外辐***确测温方法
CN107588854B (zh) 基于内置参考体的高精度测温方法
CN103063312A (zh) 一种测量物体发射率的测量***及方法
CN102889934A (zh) 实时测量温度的方法
CN105628208B (zh) 一种基于红外成像***的温度测量方法
CN105784184A (zh) 高温拉伸试验中一种温度-应变的一体化测量方法
EP2653843A1 (en) Method and system for measuring surface temperature
CN105352988A (zh) 一种建筑物外墙保温性能评估***及方法
CN103335717B (zh) 一种基于变积分模式的红外热像仪高精度抗温漂测温方法
CN103675019A (zh) 一种红外热像仪快速测量材料表面发射率的方法
CN105004754A (zh) 一种发射率的测量方法
CN103557965B (zh) 水泥回转窑温度测定及其温度场在线检测方法、装置
CN111024237A (zh) 非接触式宽温差红外温度测量方法
CN214471418U (zh) 基于双脉冲激光的温度传感器动态校准装置
CN105223230A (zh) 一种红外透波材料辐射传输特性测量方法
CN201892573U (zh) 一种近红外辐射温度计
CN104833419A (zh) 一种1-3μm准直光源辐射照度测量仪

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20150107