CN106052867A - 基于辐亮度参数的多通道led照明系统的光谱匹配方法 - Google Patents

基于辐亮度参数的多通道led照明系统的光谱匹配方法 Download PDF

Info

Publication number
CN106052867A
CN106052867A CN201610318759.3A CN201610318759A CN106052867A CN 106052867 A CN106052867 A CN 106052867A CN 201610318759 A CN201610318759 A CN 201610318759A CN 106052867 A CN106052867 A CN 106052867A
Authority
CN
China
Prior art keywords
led
channel
value
spectral
spoke
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
CN201610318759.3A
Other languages
English (en)
Other versions
CN106052867B (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.)
Qilu University of Technology
Original Assignee
Qilu University of 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 Qilu University of Technology filed Critical Qilu University of Technology
Priority to CN201610318759.3A priority Critical patent/CN106052867B/zh
Publication of CN106052867A publication Critical patent/CN106052867A/zh
Application granted granted Critical
Publication of CN106052867B publication Critical patent/CN106052867B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/28—Investigating the spectrum

Abstract

本发明公开了一种基于辐亮度参数的多通道LED照明系统的光谱匹配方法，属于LED照明技术领域。本发明首先利用光谱辐射计测量出每一LED通道最大数字驱动值时的光谱功率分布数据以及以若干数字驱动等级条件下的辐亮度值；然后利用最大数字驱动值时的光谱功率分布数据和目标光谱功率分布数据，经牛顿迭代算法计算出用以模拟目标光谱功率分布的每种LED通道之间的辐亮度比例值；最后通过三次样条插值算法计算出每种LED通道所需辐亮度比例条件下的数字驱动值，从而最终实现目标光谱。本发明解决了如何利用多通道照明系统模拟任意光谱时各LED通道所需的数字驱动值问题，为高精度图像再现、颜色复制等领域提供了必要的观察照明条件。

Description

基于辐亮度参数的多通道LED照明系统的光谱匹配方法
技术领域
本发明涉及LED照明技术领域，特别涉及一种基于辐亮度参数的多通道LED照明系统的光谱匹配方法。
背景技术
近年来，随着国家对节能减排的重视，LED照明设备已越来越多的应用在室外照明领域，而LED室内照明将成为下一个照明领域研究的热点。
LED除了低耗能、寿命长的特点外，其窄带光谱以及瞬间响应的发光特点使光谱可调成为可能。传统光源的光谱分布是确定的，但LED光源可以通过组合不同色光LED来实现任意光谱组成的光环境。
LED光谱可调光源是利用迭代方法优化由多种不同峰值波长的LED辐射叠加所形成的模拟光谱，使其与目标光谱差最小，最终实现预定的光源光谱功率分布。
由于各项光度学和色度学参数都是基于辐射度学中的光谱功率分布或光谱反射比的基础上推演而来，因此当光源的光谱功率分布确定之后，光源的辐亮度值以及光源色品坐标、色温等一系列色度学量也就相应确定下来。也即是说，如能很好的模拟任意光源或者发光体的光谱功率分布，即可以对其色品坐标、颜色三刺激值等一系列色度学参数进行匹配和复制，在实际应用中具有非常大的应用前景，尤其是对光学遥感器的辐射定标、光源显色性研究、以及光生物学等领域提供可靠的技术支持。
发明内容
为了解决现有技术中的问题，本发明提供了一种基于辐亮度参数的多通道LED照明系统的光谱匹配方法。
本发明的技术方案为：
一种基于辐亮度参数的多通道LED照明系统的光谱匹配方法，包括步骤：
1)测量标准白板或标准灰板的光谱反射比P白(λ)；
2)测量经所述标准白板或标准灰板反射的每种LED通道在最大数字驱动值255等级下的绝对光谱功率分布数据P测(λ)；
3)通过公式(1)计算多通道LED照明系统的每种LED通道在最大数字驱动值255等级下的绝对光谱功率分布数据公式(1)中，i＝1～n，其中n为LED照明系统通道数；
4)对每种LED通道的数字驱动值从0～255等级范围内采样，采样数量不小于3，且包括最大数字驱动值255；然后采用步骤2)和步骤3)的方法测量出所采样本驱动值等级下每种LED通道的辐亮度数据，将数字驱动值为Gi的样本测量所得的辐亮度值记作
5)将步骤4)所得每种LED通道在各采样数字驱动值下的辐亮度值按照公式(2)进行归一化处理；其中，将255等级数字驱动值下的辐亮度值定义为1，按照各数字驱动值下的辐亮度相对于255等级数字驱动值下的辐亮度的比例进行调整，得到在0～1之间的一系列归一化后的相对辐亮度值
R G i = I G i / I 255 i - - - ( 2 )
6)预定义可以实现目标光谱的每种LED通道辐亮度的比例值ki(0)，i＝1～n，其中n为LED照明系统通道数，然后用迭代公式(3)计算第j次迭代后ki(j)的值，式中，a为经验值比例系数，取100，为每种LED通道最大数字驱动值时的光谱功率分布，PS(λ)为目标光谱；
k i ( j ) = k i ( j - 1 ) - a × ∂ Σ λ = 380 780 [ Σ i = 1 n k i ( j - 1 ) × P LED i ( λ ) - P s ( λ ) ] 2 ∂ k i ( j - 1 ) - - - ( 3 )
7)根据判别公式(4)判断j次迭代后输出光谱与目标光谱差是否小于预设的迭代阈值，如果小于迭代阈值，迭代过程结束；如果大于迭代阈值，则进行下一次迭代，其中，ε为迭代阈值；
Σ λ = 380 780 | Σ i = 1 n k i ( j ) × P LED i ( λ ) - P s ( λ ) | ≤ ϵ - - - ( 4 )
8)判断利用步骤6)和7)的迭代算法所获得的ki(j)数组中的每一个元素是否都小于1，如果都小于1执行下一步，如果有任何一个元素大于1，降低测量面照度后，并要求返回到步骤1重新执行；
9)根据数字驱动值Gi和归一化后的相对辐亮度值的比例关系，利用插值算法计算出预测所需的辐亮度I′i所对应的每种LED通道的数字驱动值G′i。
本发明首先利用光谱辐射计测量出每一LED通道最大数字驱动值时的光谱功率分布数据以及以若干数字驱动等级条件下的辐亮度值；然后利用最大数字驱动值时的光谱功率分布数据和目标光谱功率分布数据，经牛顿迭代算法计算出用以模拟目标光谱功率分布的每种LED通道之间的辐亮度比例值；最后通过三次样条插值算法计算出每种LED通道所需辐亮度比例条件下的数字驱动值，从而最终实现目标光谱。
作为优选方案，步骤1)中，采用分光光度计测量标准白板或标准灰板的光谱反射比P白(λ)。
作为优选方案，步骤2)中，采用光谱辐射计测量经所述标准白板或标准灰板反射的每种LED通道在最大数字驱动值255等级下的绝对光谱功率分布数据P测(λ)。
进一步地，步骤4)中，采用光谱辐射计测量所采样本驱动值等级下每种LED通道的辐亮度数据
作为优选方案，步骤4)中对每种LED通道的数字驱动值从0～255等级范围内均匀采样，采用数目为9-12个。
作为优选方案，步骤7)中，迭代阈值取值为0.0001。
作为优选方案，步骤9)中所述插值算法为三次样条插值算法。
本发明的有益效果为：
本发明的光谱匹配方法适用于解决LED多通道照明系统的光谱匹配问题，利用辐亮度作为重要参数指标，对辐射光谱可覆盖或基本覆盖可见光范围且大于3通道的任意通道数的LED照明系统均可由本发明方法进行光谱匹配，精度高，且简便实用。
本发明解决了如何利用多通道照明系统模拟任意光谱时各LED通道所需的数字驱动值问题，此问题的解决为高精度图像再现、颜色复制等科技领域提供了必要的观察照明条件。
具体实施方式
以11通道LED照明系统模拟测量面照度1000lx条件下的标准照明体A的光谱功率分布为例，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例1
一种基于辐亮度参数的多通道LED照明系统的光谱匹配方法，包括步骤：
1)利用X-Rite SpectralEye分光光度计测量X-Rite ColorChecker Passport标准灰板的光谱反射比P白(λ)，结果如表1所示；
2)利用Topcon SR-3A光谱辐射计测量经X-Rite ColorChecker Passport标准灰板反射的11种LED通道在最大数字驱动值255等级下的绝对光谱功率分布数据P测(λ)，11种LED通道的光谱数据如表2-1、表2-2、表2-3、表2-4所示；
3)通过公式(1)计算多通道LED照明系统的每种LED通道在最大数字驱动值255等级下的绝对光谱功率分布数据公式(1)中，i＝1～n，其中n为LED照明系统通道数
绝对光谱功率分布数据根据表1与表2-1、表2-2、表2-3、表2-4中的数据计算即可；
4)对每种LED通道的数字驱动值进行采样，采样值分别为255、235、215、195、175、150、100、80、60、40、25；采用Topcon SR-3A光谱辐射计以步骤2)和步骤3)的方法测量出所采样本驱动值等级下11种LED通道的辐亮度数据，将数字驱动值为Gi的样本测量所得的辐亮度值记作11种LED通道的辐亮度值如表3-1、表3-2所示；
5)将步骤4)所得11种LED通道在11个采样数字驱动值下的辐亮度值按照公式(2)进行归一化处理；其中，将255等级数字驱动值下的辐亮度值定义为1，按照各数字驱动值下的辐亮度相对于255等级数字驱动值下的辐亮度的比例进行调整，得到在0～1之间的一系列归一化后的相对辐亮度值相对辐亮度值如表4-1、表4-2所示；
R G i = I G i / I 255 i - - - ( 2 )
6)预定义11种LED通道辐亮度的比例值ki(0)为1:1:1:1:1:1:1:1:1:1:1，用以实现标准照明体A的光谱功率分布，然后用迭代公式(3)计算第j次迭代后的ki(j)值，本例共迭代20次，ki的20次迭代取值如表5-1、表5-2所示，a为经验系数，取值100；为每种LED通道最大数字驱动值时的光谱功率分布，PS(λ)为目标光谱；
k i ( j ) = k i ( j - 1 ) - a × ∂ Σ λ = 380 780 [ Σ i = 1 n k i ( j - 1 ) × P LED i ( λ ) - P s ( λ ) ] 2 ∂ k i ( j - 1 ) - - - ( 3 )
7)根据判别公式(4)判断20次迭代后输出光谱与目标光谱差是否小于预设的迭代阈值ε＝0.0001；本实施例中，19次迭代后的光谱差大于阈值，20次迭代后光谱差为4.6903e-05，即，4.6903×10-5；
Σ λ = 380 780 | Σ i = 1 n k i ( j ) × P LED i ( λ ) - P s ( λ ) | ≤ ϵ - - - ( 4 )
8)获得的ki(20)中代表11个LED通道的元素值全部小于1，继续执行下一步；
9)根据数字驱动值Gi和归一化后的相对辐亮度值RGi的比例关系，利用三次样条插值算法和取整算法计算出预测所需11种LED通道的数字驱动值G′i，分别为156、139、16、183、186、128、259、164、39、72、16。
表1X-Rite ColorChecker Passport标准灰板的光谱反射比
表2-1 11种LED通道经灰板反射后的绝对光谱功率分布
表2-2 11种LED通道经灰板反射后的绝对光谱功率分布
表2-3 11种LED通道经灰板反射后的绝对光谱功率分布
表2-4 11种LED通道经灰板反射后的绝对光谱功率分布
表3-1 11种LED通道的辐亮度值
表3-2 11种LED通道的辐亮度值
表4-1 11种LED通道归一化的辐亮度
表4-2 11种LED通道归一化的辐亮度
表5-1 11种LED通道20次迭代值
迭代次数 K1 K2 K3 K4 K5 K6
1 0.614551408 0.572259283 0.405591542 0.777636015 0.95004936 0.894271622
2 0.472454137 0.419577868 0.236633741 0.689876084 0.915882229 0.817482005
3 0.390842598 0.330835873 0.137896099 0.638034341 0.891472045 0.758692518
4 0.344294431 0.279279622 0.079173289 0.607400209 0.873317636 0.711562044
5 0.318077638 0.249421599 0.043038898 0.589468612 0.859333696 0.672306037
6 0.304823758 0.233384619 0.021225257 0.579582903 0.848175733 0.638481439
7 0.2995051 0.225832711 0.007981182 0.57468302 0.838994688 0.608541038
8 0.299088343 0.223538311 0 0.57288638 0.831242593 0.581482574
9 0.301718612 0.224520057 0 0.572643947 0.824492492 0.556521084
10 0.305887367 0.227142299 0 0.573462613 0.818557701 0.533309326
11 0.311103081 0.230880013 0 0.574986066 0.813297419 0.511583489
12 0.31701172 0.235353838 0 0.576955745 0.808604559 0.491141152
13 0.323357544 0.240288924 0 0.579183948 0.804396946 0.47182534
14 0.32995595 0.245486324 0 0.581534871 0.800610888 0.453512868
15 0.336674254 0.250802668 0 0.58391098 0.797196445 0.436105763
16 0.34341786 0.25613549 0 0.586243072 0.794113936 0.419524919
17 0.350120181 0.261412504 0 0.588482948 0.791331344 0.403705399
18 0.356735232 0.266583699 0 0.59059798 0.788822402 0.388592944
19 0.363232146 0.271615469 0 0.592567045 0.786565163 0.374141385
20 0.369591116 0.276486257 0 0.594377495 0.784540945 0.360310716
表5-2 11种LED通道20次迭代值

Claims (7)

1.一种基于辐亮度参数的多通道LED照明系统的光谱匹配方法，其特征在于，包括步骤：
1)测量标准白板或标准灰板的光谱反射比P白(λ)；
2)测量经所述标准白板或标准灰板反射的每种LED通道在最大数字驱动值255等级下的绝对光谱功率分布数据P测(λ)；
3)通过公式(1)计算多通道LED照明系统的每种LED通道在最大数字驱动值255等级下的绝对光谱功率分布数据公式(1)中，i＝1～n，其中n为LED照明系统通道数；
4)对每种LED通道的数字驱动值从0～255等级范围内采样，采样数量不小于3，且包括最大数字驱动值255；然后采用步骤2)和步骤3)的方法测量出所采样本驱动值等级下每种LED通道的辐亮度数据，将数字驱动值为Gi的样本测量所得的辐亮度值记作
5)将步骤4)所得每种LED通道在各采样数字驱动值下的辐亮度值按照公式(2)进行归一化处理；其中，将255等级数字驱动值下的辐亮度值定义为1，按照各数字驱动值下的辐亮度相对于255等级数字驱动值下的辐亮度的比例进行调整，得到在0～1之间的一系列归一化后的相对辐亮度值
R G i = I G i / I 255 i - - - ( 2 )
6)预定义可以实现目标光谱的每种LED通道辐亮度的比例值ki(0)，i＝1～n，其中n为LED照明系统通道数，然后用迭代公式(3)计算第j次迭代后ki(j)的值，式中，a为经验值比例系数，取100，为每种LED通道最大数字驱动值时的光谱功率分布，PS(λ)为目标光谱；
k i ( j ) = k i ( j - 1 ) - a × ∂ Σ λ = 380 780 [ Σ i = 1 n k i ( j - 1 ) × P LED i ( λ ) - P s ( λ ) ] 2 ∂ k i ( j - 1 ) - - - ( 3 )
7)根据判别公式(4)判断j次迭代后输出光谱与目标光谱差是否小于预设的迭代阈值，如果小于迭代阈值，迭代过程结束；如果大于迭代阈值，则进行下一次迭代，其中，ε为迭代阈值；
Σ λ = 380 780 | Σ i = 1 n k i ( j ) × P LED i ( λ ) - P s ( λ ) | ≤ ϵ - - - ( 4 )
8)判断利用步骤6)和7)的迭代算法所获得的ki(j)数组中的每一个元素是否都小于1，如果都小于1执行下一步，如果有任何一个元素大于1，降低测量面照度后，并要求返回到步骤1重新执行；
9)根据数字驱动值Gi和归一化后的相对辐亮度值的比例关系，利用插值算法计算出预测所需的辐亮度I′i所对应的每种LED通道的数字驱动值G′i。
2.如权利要求1所述基于辐亮度参数的多通道LED照明系统的光谱匹配方法，其特征在于：步骤1)中，采用分光光度计测量标准白板或标准灰板的光谱反射比P白(λ)。
3.如权利要求1或2所述基于辐亮度参数的多通道LED照明系统的光谱匹配方法，其特征在于：步骤2)中，采用光谱辐射计测量经所述标准白板或标准灰板反射的每种LED通道在最大数字驱动值255等级下的绝对光谱功率分布数据P测(λ)。
4.如权利要求3所述基于辐亮度参数的多通道LED照明系统的光谱匹配方法，其特征在于：步骤4)中，采用光谱辐射计测量所采样本驱动值等级下每种LED通道的辐亮度数据
5.如权利要求1所述基于辐亮度参数的多通道LED照明系统的光谱匹配方法，其特征在于：步骤4)中对每种LED通道的数字驱动值从0～255等级范围内均匀采样，采用数目为9-12个。
6.如权利要求1所述基于辐亮度参数的多通道LED照明系统的光谱匹配方法，其特征在于：步骤7)中，迭代阈值取值为0.0001。
7.如权利要求1所述基于辐亮度参数的多通道LED照明系统的光谱匹配方法，其特征在于：步骤9)中所述插值算法为三次样条插值算法。
CN201610318759.3A 2016-05-12 2016-05-12 基于辐亮度参数的多通道led照明系统的光谱匹配方法 Expired - Fee Related CN106052867B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610318759.3A CN106052867B (zh) 2016-05-12 2016-05-12 基于辐亮度参数的多通道led照明系统的光谱匹配方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610318759.3A CN106052867B (zh) 2016-05-12 2016-05-12 基于辐亮度参数的多通道led照明系统的光谱匹配方法

Publications (2)

Publication Number Publication Date
CN106052867A true CN106052867A (zh) 2016-10-26
CN106052867B CN106052867B (zh) 2017-12-22

Family

ID=57176289

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610318759.3A Expired - Fee Related CN106052867B (zh) 2016-05-12 2016-05-12 基于辐亮度参数的多通道led照明系统的光谱匹配方法

Country Status (1)

Country Link
CN (1) CN106052867B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110708797A (zh) * 2019-10-22 2020-01-17 滨州学院 一种用于三色led光源的精确调光方法及系统
CN111006762A (zh) * 2019-12-10 2020-04-14 东华理工大学 一种基于led分段拟合太阳光谱的空间遥感仪器星上定标方法
CN111654941A (zh) * 2020-06-23 2020-09-11 滨州学院 一种飞机客舱模拟日光照明的方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101988894A (zh) * 2009-07-31 2011-03-23 恩德莱斯和豪瑟尔测量及调节技术分析仪表两合公司 用于确定液体样品的参数的方法
US20120098924A1 (en) * 2008-11-26 2012-04-26 Honeywell International Inc. Signal spectra detection system
CN103208102A (zh) * 2013-03-29 2013-07-17 上海交通大学 一种基于稀疏表示的遥感图像融合方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120098924A1 (en) * 2008-11-26 2012-04-26 Honeywell International Inc. Signal spectra detection system
CN101988894A (zh) * 2009-07-31 2011-03-23 恩德莱斯和豪瑟尔测量及调节技术分析仪表两合公司 用于确定液体样品的参数的方法
CN103208102A (zh) * 2013-03-29 2013-07-17 上海交通大学 一种基于稀疏表示的遥感图像融合方法

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110708797A (zh) * 2019-10-22 2020-01-17 滨州学院 一种用于三色led光源的精确调光方法及系统
CN110708797B (zh) * 2019-10-22 2021-04-02 滨州学院 一种用于三色led光源的精确调光方法及系统
CN111006762A (zh) * 2019-12-10 2020-04-14 东华理工大学 一种基于led分段拟合太阳光谱的空间遥感仪器星上定标方法
CN111006762B (zh) * 2019-12-10 2021-09-24 东华理工大学 一种基于led分段拟合太阳光谱的空间遥感仪器星上定标方法
CN111654941A (zh) * 2020-06-23 2020-09-11 滨州学院 一种飞机客舱模拟日光照明的方法及系统
CN111654941B (zh) * 2020-06-23 2022-04-26 滨州学院 一种飞机客舱模拟日光照明的方法及系统

Also Published As

Publication number Publication date
CN106052867B (zh) 2017-12-22

Similar Documents

Publication Publication Date Title
Lynn et al. Color rendering properties of semi-transparent thin-film PV modules
CN105136432B (zh) 基于主客观实验数据的led照明质量评价方法及系统
CN105973470B (zh) 一种多色led实现色度限制的光谱匹配方法
CN106017868B (zh) 基于亮度参数的多通道led照明系统的光谱匹配方法
CN105102946B (zh) 用于多流颜色匹配的系统和方法
EP3324160B1 (en) Standard illuminant color matching box
CN106053024B (zh) 一种面向单色系物体的led光源喜好度预测方法
CN113673389B (zh) 光源光谱功率分布相关的绘画照明视觉评价方法
CN104991988B (zh) 基于多颗单色大功率led实现类日光光源的方法
CN102507151A (zh) 一种基于工效学实验的阅读台灯评价方法
CN106052867B (zh) 基于辐亮度参数的多通道led照明系统的光谱匹配方法
CN109511195A (zh) 适用于中国脆弱文物照明的led光谱调节方法
US20090103075A1 (en) Method For Measuring Chroma Value By a Colorimeter
CN105934020A (zh) 一种多色led匹配光谱和照度的方法
CN110677954B (zh) 一种超三色led光源的精确调光方法及系统
CN106353264A (zh) 适用于绘画色彩保护照明的白光led光谱获取方法
CN103196559A (zh) 电子罗维朋测色法及其检测系统
CN104062010B (zh) 一种优化定标算法的分光光源颜色照度测量仪器
CN109870447A (zh) 判定光源对中国脆弱文物照明损伤度的方法
CN106061000B (zh) 基于光谱峰值辐亮度参数的多通道照明系统光谱匹配方法
Burgos-Fernandez et al. Spectrally tunable light source based on light-emitting diodes for custom lighting solutions
Ajrina et al. Spectral reflectance and chromaticity differences of various colors of interior finishing material samples under tunable LED lamps
CN106323849A (zh) 衡量照明光源对传统绘画色彩损伤的检测及评价方法
CN106908150A (zh) 嵌入式系统相关色温快速计算方法
Dang et al. Illumination in museums: Four-primary white LEDs to optimize the protective effect and color quality

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

Granted publication date: 20171222

CF01 Termination of patent right due to non-payment of annual fee