CN108324295B - 能谱寻峰方法、装置和计算机存储介质 - Google Patents
能谱寻峰方法、装置和计算机存储介质 Download PDFInfo
- Publication number
- CN108324295B CN108324295B CN201711363000.8A CN201711363000A CN108324295B CN 108324295 B CN108324295 B CN 108324295B CN 201711363000 A CN201711363000 A CN 201711363000A CN 108324295 B CN108324295 B CN 108324295B
- Authority
- CN
- China
- Prior art keywords
- energy spectrum
- target
- address
- determining
- target energy
- 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
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 312
- 238000000034 method Methods 0.000 title claims abstract description 49
- 230000002285 radioactive effect Effects 0.000 claims abstract description 37
- 230000006870 function Effects 0.000 claims description 22
- 230000005484 gravity Effects 0.000 claims description 20
- 238000001514 detection method Methods 0.000 claims description 19
- 239000013078 crystal Substances 0.000 claims description 17
- 238000004422 calculation algorithm Methods 0.000 claims description 16
- 230000005855 radiation Effects 0.000 claims description 10
- 238000001914 filtration Methods 0.000 claims description 9
- 238000004590 computer program Methods 0.000 claims description 2
- 239000003814 drug Substances 0.000 abstract description 2
- 238000002600 positron emission tomography Methods 0.000 description 16
- 230000006854 communication Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 230000003595 spectral effect Effects 0.000 description 6
- 238000004891 communication Methods 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 4
- 238000001816 cooling Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000009499 grossing Methods 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 230000003902 lesion Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000000700 radioactive tracer Substances 0.000 description 2
- 230000005236 sound signal Effects 0.000 description 2
- 206010028980 Neoplasm Diseases 0.000 description 1
- 230000005255 beta decay Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000005251 gamma ray Effects 0.000 description 1
- 230000005253 gamme decay Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 238000010079 rubber tapping Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000010897 surface acoustic wave method Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/037—Emission tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating thereof
- A61B6/587—Alignment of source unit to detector unit
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Engineering & Computer Science (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Measurement Of Radiation (AREA)
Abstract
本发明公开了一种能谱寻峰方法、装置和计算机存储介质,属于医学技术领域。方法包括:通过探测器探测放射源,得到所述放射源的目标能谱,所述目标能谱包括多个能谱点,每个能谱点用于指示道址和计数的对应关系;从所述目标能谱中确定所述目标能谱的最大计数对应的第一道址;根据所述第一道址,从所述目标能谱中确定包括所述第一道址的目标区域;从所述目标区域中,确定所述目标能谱的第二道址;根据所述第二道址,确定所述目标能谱的参考峰值。在发明实施例中,寻找峰值对应的第一道址,然后基于第一道址确定目标区域,从目标区域中确定出该目标能谱的参考峰值,从而提高了寻峰准确性。
Description
技术领域
本发明涉及医学技术领域,特别涉及一种能谱寻峰方法、装置和计算机存储介质。
背景技术
随着医学技术的发展,PET(Positron Emission Tomography,正电子发射断层显像)探测器能够检测人体的病灶,例如肿瘤。然而PET探测器可能有一定的探测误差,因此,在使用PET探测器检测人体病灶之前,需要先通过该PET探测器探测放射源,获取该放射源对应的能谱的参考峰值,基于该参考峰值和该放射源的基准峰值,确定该PET探测器的漂移值,根据该漂移值对PET探测器进行调整,以补偿该探测误差。
由于PET探测器所在环境的温度对PET探测器的漂移影响较大;因此,在现有技术中在对获取该放射源对应的能谱的参考峰值时,一般使用水冷或者风冷系统对该PET探测器进行温控,然后从该能谱中找到最大计数对应的道址,将该道址作为该能谱的参考峰值。
在实现本发明的过程中,发明人发现现有技术至少存在以下问题:
由于PET探测器产生的热量并不能完全被有效带走,不同探测器本底形状和能峰形状并不完全一致。在探测器尺寸较大,性能发生非均匀变化的情况下,随着温度的漂移,会造成探测器探测到的峰宽和形状发生一定的变化,因此,上述方法中得到的参考峰值不准确。
发明内容
本发明提供了一种能谱寻峰方法、装置和计算机存储介质,可以提高在探测到的核事件计数率稳定的情况下确定出的参考峰值的准确性。技术方案如下:
第一方面,本发明提供了一种能谱寻峰方法,所述方法包括:
通过探测器探测放射源,得到所述放射源的目标能谱,所述目标能谱包括多个能谱点,每个能谱点用于指示道址和计数的对应关系;
从所述目标能谱中确定所述目标能谱的最大计数对应的第一道址;
根据所述第一道址,从所述目标能谱中确定包括所述第一道址的目标区域;
从所述目标区域中,确定所述目标能谱的第二道址;
根据所述第二道址,确定所述目标能谱的参考峰值。
在一个可能的实现方式中,所述根据所述第一道址,从所述目标能谱中确定包括所述第一道址的目标区域,包括:
从所述第一道址开始向两侧进行积分,寻找积分值不低于第一预设阈值的第一能谱点和第二能谱点;
将所述目标能谱中所述第一能谱点和所述第二能谱点之间的能谱组成所述目标区域。
在一个可能的实现方式中,所述从所述目标区域中,确定所述目标能谱的第二道址,包括:
确定所述目标区域的重心;
将所述重心对应的道址确定为所述目标能谱的第二道址。
在一个可能的实现方式中,所述根据所述第二道址,确定所述目标能谱的参考峰值之后,所述方法还包括:
获取所述目标能谱的基准峰值;
根据所述基准峰值和所述参考峰值,确定所述探测器的漂移值;
根据所述漂移值,对所述探测器进行调整。
在一个可能的实现方式中,所述第二道址包括至少一个第二子道址,所述根据所述第二道址,确定所述目标能谱的参考峰值,包括:
当所述第二道址包括一个第二子道址时,将所述第二子道址作为所述目标能谱的参考峰值;
当所述第二道址包括多个第二子道址时,根据每个第二子道址,通过预设迭代算法进行迭代运算,输出所述目标能谱的参考峰值。
在一个可能的实现方式中,所述从所述目标能谱中确定所述目标能谱的最大计数对应的第一道址之前,所述方法还包括:
通过数字滤波器,将所述目标能谱中的高频成分过滤掉。
在一个可能的实现方式中,所述放射源为所述探测器自身具有衰变功能的晶体。
第二方面,本发明提供了一种能谱寻峰装置,所述装置包括:
探测模块,用于通过探测器探测放射源,得到所述放射源的目标能谱,所述目标能谱包括多个能谱点,每个能谱点用于指示道址和计数的对应关系;
第一确定模块,用于从所述目标能谱中确定所述目标能谱的最大计数对应的第一道址;
第二确定模块,用于根据所述第一道址,从所述目标能谱中确定包括所述第一道址的目标区域;
第三确定模块,用于从所述目标区域中,确定所述目标能谱的第二道址;
第四确定模块,用于根据所述第二道址,确定所述目标能谱的参考峰值。
在一个可能的实现方式中,所述第二确定模块,包括:
寻找单元,用于从所述第一道址开始向两侧进行积分,寻找积分值不低于第一预设阈值的第一能谱点和第二能谱点;
组成单元,用于将所述目标能谱中所述第一能谱点和所述第二能谱点之间的能谱组成所述目标区域。
在一个可能的实现方式中,所述第三确定模块,还用于确定所述目标区域的重心,将所述重心对应的道址确定为所述目标能谱的第二道址。
在一个可能的实现方式中,所述装置还包括:
调整模块,用于获取所述目标能谱的基准峰值,根据所述基准峰值和所述参考峰值,确定所述探测器的漂移值,根据所述漂移值,对所述探测器进行调整。
在一个可能的实现方式中,所述第二道址包括至少一个第二子道址,
所述第四确定模块,还用于当所述第二道址包括一个第二子道址时,将所述第二子道址作为所述目标能谱的参考峰值;当所述第二道址包括多个第二子道址时,根据每个第二子道址,通过预设迭代算法进行迭代运算,输出所述目标能谱的参考峰值。
在一个可能的实现方式中,所述装置还包括:
滤波模块,用于通过数字滤波器,将所述目标能谱中的高频成分过滤掉。
在一个可能的实现方式中,所述放射源为所述探测器自身具有衰变功能的晶体。
第三方面,本发明提供了一种能谱寻峰装置,所述装置包括:
处理器;
用于存储处理器可执行指令的存储器;
其中,所述处理器被配置为:
通过探测器探测放射源,得到所述放射源的目标能谱,所述目标能谱包括多个能谱点,每个能谱点用于指示道址和计数的对应关系;
从所述目标能谱中确定所述目标能谱的最大计数对应的第一道址;
根据所述第一道址,从所述目标能谱中确定包括所述第一道址的目标区域;
从所述目标区域中,确定所述目标能谱的第二道址;
根据所述第二道址,确定所述目标能谱的参考峰值。
第四方面,本发明提供了一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序,所述程序被处理器执行时实现以下步骤:
通过探测器探测放射源,得到所述放射源的目标能谱,所述目标能谱包括多个能谱点,每个能谱点用于指示道址和计数的对应关系;
从所述目标能谱中确定所述目标能谱的最大计数对应的第一道址;
根据所述第一道址,从所述目标能谱中确定包括所述第一道址的目标区域;
从所述目标区域中,确定所述目标能谱的第二道址;
根据所述第二道址,确定所述目标能谱的参考峰值。
在发明实施例中,寻找峰值对应的第一道址,然后基于第一道址确定目标区域,从目标区域中确定出该目标能谱的参考峰值,从而提高了寻峰准确性。
附图说明
图1是本发明实施例提供的一种能谱寻峰的方法流程图;
图2是本发明实施例提供的一种能谱寻峰的方法流程图;
图3是本发明实施例提供的一种数字滤波的示意图;
图4是本发明实施例提供的一种目标区域的示意图;
图5是本发明实施例提供的一种包含多个目标子区域的示意图;
图6是本发明实施例提供的一种能谱寻峰装置装置结构示意图;
图7是本发明实施例提供的一种能谱寻峰装置装置结构示意图;
图8是本发明实施例提供的一种能谱寻峰装置装置结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
本发明实施例提供了一种能谱寻峰方法,该方法的执行主体可以为控制终端;参见图1,该方法包括:
步骤101:通过探测器探测放射源,得到放射源的目标能谱,目标能谱包括多个能谱点,每个能谱点用于指示道址和计数的对应关系。
步骤102:从目标能谱中确定目标能谱的最大计数对应的第一道址。
步骤103:根据第一道址,从目标能谱中确定包括第一道址的目标区域。
步骤104:从目标区域中,确定目标能谱的第二道址。
步骤105:根据第二道址,确定目标能谱的参考峰值。
在一个可能的实现方式中,根据第一道址,从目标能谱中确定包括第一道址的目标区域,包括:
从第一道址开始向两侧进行积分,寻找积分值不低于第一预设阈值的第一能谱点和第二能谱点;
将目标能谱中第一能谱点和第二能谱点之间的能谱组成目标区域。
在一个可能的实现方式中,从目标区域中,确定目标能谱的第二道址,包括:
确定目标区域的重心;
将重心对应的道址确定为目标能谱的第二道址。
在一个可能的实现方式中,根据第二道址,确定目标能谱的参考峰值之后,方法还包括:
获取目标能谱的基准峰值;
根据基准峰值和参考峰值,确定探测器的漂移值;
根据漂移值,对探测器进行调整。
在一个可能的实现方式中,第二道址包括至少一个第二子道址,根据第二道址,确定目标能谱的参考峰值,包括:
当第二道址包括一个第二子道址时,将第二子道址作为目标能谱的参考峰值;
当第二道址包括多个第二子道址时,根据每个第二子道址,通过预设迭代算法进行迭代运算,输出目标能谱的参考峰值。
在一个可能的实现方式中,从目标能谱中确定目标能谱的最大计数对应的第一道址之前,方法还包括:
通过数字滤波器,将目标能谱中的高频成分过滤掉。
在一个可能的实现方式中,放射源为探测器自身具有衰变功能的晶体。
在发明实施例中,寻找峰值对应的第一道址,然后基于第一道址确定目标区域,从目标区域中确定出该目标能谱的参考峰值,从而提高了寻峰准确性。
本发明实施例提供了一种能谱寻峰方法,该方法的执行主体可以为控制终端;参见图2,该方法包括:
步骤201:控制终端通过探测器探测放射源,得到该放射源的目标能谱,目标能谱包括多个能谱点,每个能谱点用于指示道址和计数的对应关系。
探测器可以为PET探测器或者CT(Computed Tomography,计算机断层扫描)探测器;放射源可以为放射性示踪剂或者探测器自身具有衰变功能的晶体。放射性示踪剂可以为18F-FDG(2-Fluorine-18-Fluoro-2-deeoxy-D-glucose,2-氟-18-氟-2-脱氧-D-葡萄糖)溶液。自身具有衰变功能的晶体可以为含176LU的LYSO晶体或者含138La的晶体选自LaBr晶体。
本步骤可以通过以下步骤(1)至(2)实现,包括:
(1):控制终端控制探测器连续探测放射源达到指定时长,得到探测能谱。
控制终端确定目标能量的能量标识,向探测器发送探测指令,该探测指令包括指定时长和该目标能量的能量标识。探测器接收探测命令,根据该目标能量的能量标识,连续探测该目标能量达到指定时长,得到该目标能量的探测能谱。
放射源能够放射至少一种能量,在本步骤中控制终端可以控制探测器采集放射源全部或者部分能量对应的探测能谱。因此,目标能量可以为一个或者多个能量。并且,控制终端也可以从放射源放射的能量中选择便于探测能峰的能量标识。相应的,控制终端确定目标能量的能量标识的步骤可以为:
控制终端根据该放射源的标识,从放射源标识和能量标识集合中获取该放射源对应的能量标识集合,将该能量集合标识中的能量标识确定为目标能量的能量标识。其中,该能量标识集合中包括便于探测能峰的能量标识。能量标识可以为目标能量的名称或者大小。
例如,当该放射源为含176LU的LYSO晶体时,该晶体具有自衰变性能,其主要衰变方式为β衰变瞬态伴随γ衰变,发出γ射线能量有三种,分别为88keV、202keV、307keV。其中,202keV和307keV为便于探测能峰的能量,则176LU的LYSO晶体对应的能量标识集合中包括202keV和307keV。当然如果仅探测一个能量的能谱时,该176LU的LYSO晶体对应的能量标识集合中可以包括307keV。
在本发明实施例中,控制终端可以探测放射源发出的所有能量的能谱,从而基于所有能量的能谱确定射线的峰值,从而提高准确性。控制终端也可以不探测放射源发出的所有能量的能谱,只探测能峰便于区分的部分能量的能谱,从而提高效率。
指定时长可以根据需要进行设置并更改,在本发明实施例中不对指定时长做具体限定;例如,指定时长为3分钟或者5分钟。为了提高后续对探测器进行调整的效率,在本发明实施例中,可以通过低计数率进行寻峰,从而提高寻峰效率。因此,指定时长小于10分钟。
需要说明的是,如果目标能量包括多个能量时,该探测能谱中包括多个能量的能谱。
(2):控制终端从该探测能谱中提取目标能谱。
目标能谱可以为一个,也即以下第一种方式,目标能谱也可以为多个,也即以下第二种方式。相应的,对于第一种方式,本步骤可以为:
控制终端将探测能谱进行等份数据分割,得到多个能谱,从多个能谱中选择一个能谱作为目标能谱。
对于第二种方式,本步骤可以为:
控制终端将探测能谱进行等份数据分割,得到多个能谱,将多个能谱分别作为目标能谱。
例如,终端连续采集了200s的数据,在本步骤中,控制终端将200s数据等份划分为4分,则得到4个能谱。控制终端可以从4个能谱中选择一个能谱作为目标能谱。控制终端也可以将4个能谱分别作为目标能谱。
优选的,在本发明实施例中,以探测器自身具有衰变功能的晶体作为放射源。这样可以无需购买放射性示踪剂,也无需对放射源进行每日质控,较少成本以及额外操作。并且合理利用了PET探测器自身晶体衰变环境的数据收集,可以实现短时间内对PET探测器即时状态的检查并且有较高的计算精度。
在本发明实施例中,由于采用低计数率进行寻峰,因此,控制终端得到目标能谱之后,可以通过步骤202对能谱的谱线进行平滑处理,从而提高寻峰准确性,也即执行完步骤201之后,执行步骤202。当然,控制终端得到目标能谱之后,不执行步骤202,直接执行步骤203。
步骤202:控制终端通过数字滤波器,将该目标能谱中的高频成分过滤掉。
在本发明实施例中,通过数字滤波器,对目标能谱进行平滑处理,从而有效去除目标能谱中的高频成分,并且还能实现谱线形状的跟随。并且在本领域中对谱线进行平滑处理时,一般使用多点平滑处理,然而多点平滑处理,不能实现谱线形状的跟随。而本发明实施例中的通过数字滤波器进行平滑处理,不仅可以有效滤除高频成分,还能实现谱线形状的跟随,例如参见图3。
步骤203:控制终端从该目标能谱中确定该目标能谱的最大计数对应的第一道址。
目标能谱中包括多个能谱点,每个能谱点用于指示计数和道址的对应关系。控制终端确定出目标能谱之后,从该目标能谱中确定最大计数,将该最大计数对应的道址确定为第一道址。
如果该目标能谱包括多个目标能量的能谱,则在本步骤中,控制终端从该目标能谱中确定每个目标能谱的目标子能谱,对于每个目标能量的目标子能谱,控制终端确定该目标子能谱中的最大计数对应的第一子道址,将每个目标子能谱中的第一子道址组成第一道址。
例如,该目标能谱包括202keV的子能谱和307keV的子能谱;则在本步骤中,控制终端从该目标能谱中分别确定出202keV对应的第一子道址,以及307keV对应的第一子道址,将202keV对应的第一子道址,以及307keV对应的第一子道址组成第一道址。
步骤204:控制终端根据第一道址,从该目标能谱中确定包括第一道址的目标区域。
在本发明实施例中,控制终端可以通过最大值法、对称零面积法、求导法、高斯拟合法或者重心法进行寻峰。鉴于目标能谱的谱线的良好跟随,采用重心法能够保持单峰的时效性。因此,在本发明实施例中,以通过重心法进行寻峰为例进行说明。相应的,本步骤可以通过以下步骤(1)至(4)实现,包括:
(1):控制终端从第一道址开始向两侧进行积分,寻找积分值不低于第一预设阈值的第一能谱点和第二能谱点。
控制终端两侧为第一能谱点和第二能谱点,对两个能谱点中间的计数减掉基线后以指定计数率进行积分,如果积分值低于预设阈值,那么计数较高的能谱点向外侧移动,重新计算积分值,如果积分值高于预设阈值,那么当前的能谱点为第一能谱点和第二能谱点。
指定计数率可以根据需要进行设置并更改,在本发明实施例中,对指定计数率不作具体限定;例如,指定计数率可以为5万或者10万等。预设阈值也可以根据需要进行设置并更改,在本发明实施例中,对预设阈值不作具体限定,例如预设阈值可以为10或者20等。
(2):控制终端将目标能谱中第一能谱点和第二能谱点之间的能谱组成目标区域。
例如,参见图4,控制终端将第一能谱点和第二能谱点之间的能谱组成目标区域。
需要说明的是,如果该目标能谱中包括多个目标能量的能谱,则第一道址包括多个第一子道址,一个目标能量对应一个第一子道址。则在本步骤中,控制终端根据每个第一子道址,从该目标能谱中确定多个目标能量对应的多个目标区域。
例如,参见图5,目标区域中包括202keV对应的目标子区域和307keV对应的目标子区域;左侧的目标子区域为202keV对应的目标子区域,右侧的目标子区域为307keV对应的目标子区域。
步骤205:控制终端从目标区域中,确定目标能谱的第二道址。
控制终端确定目标区域的重心;将该重心对应的道址确定为该目标能谱的第二道址。
当目标区域中包括多个目标子区域时,对于每个目标子区域,控制终端确定该目标子区域的重心,将该目标子区域的重心对应的道址确定为该目标子区域的第二子道址,将每个第二子道址组成第二道址。
例如,目标区域中包括202keV对应的目标子能谱和307keV对应的目标子能谱,则控制终端分别确定202keV对应的第二子道址,307keV对应的第二子道址,将202keV对应的第二子道址和307keV对应的第二子道址组成第二道址。
步骤206:控制终端根据第二道址,确定目标能谱的参考峰值。
第二道址包括至少一个第二子道址。当该目标能谱包括一个目标能量的子能谱时,第二道址包括一个第二子道址。当该目标能谱包括多个目标能量的子能谱时,第二道址包括多个第二子道址。
当第二道址包括一个第二子道址时,将第二子道址作为该目标能谱的参考峰值。当第二道址包括多个第二子道址时,根据每个第二子道址,通过预设迭代算法进行迭代运算,输出该目标能谱的参考峰值。
预设迭代算法可以根据需要进行设置并更改,在本发明实施例中,对预设迭代算法不作具体限定。例如,预设迭代算法可以为卡尔曼滤波算法或者权重加权算法。
例如,当预设迭代算法为卡尔曼滤波算法,则控制终端根据每个第二子道址,通过预设迭代算法进行迭代运算,输出该目标能谱的参考峰值的步骤可以为:
控制终端获取每个目标能量对应的基准子道址;根据每个目标能量对应的基准子道址和每个能量的第二子道址,确定每个能量的基准子道址与第二子道址之间的差值,利用高斯函数拟合每个能量对应的差值的分布,根据每个能量对应的差值的分布确定两次寻峰之间的差值分布的分辨率,然后通过卡尔曼滤波算法对两次寻峰之间的差值分布的分辨率进行方差模型的建立,然后更新迭代,得到分辨率小于预设分辨率的参考峰值。
在本发明实施例中,在PET系统上进行无源状态下的连续采集、寻峰、调整过程的迭代测试,能够将峰位有效控制在8‰(1024道分辨)以内。
在本发明实施例中,可以结合水冷或者风冷系统对该PET探测器进行温控,从而进一步提高寻峰准确性,并且对于维持系统性能稳定性更加有效。
步骤207:控制终端获取该目标能谱的基准峰值,根据该基准峰值和该参考峰值,确定该探测器的漂移值。
控制终端中存储该目标能谱的基准能谱,在本步骤中控制终端直接获取已存储的目标能谱的基准峰值,确定该参考峰值与该基准峰值之间的第一差值,将第一差值确定为该探测器的漂移值。当然,控制终端也可以确定该基准峰值与该参考峰值之间的第二差值,将第二差值确定为该探测器的偏移值。
例如,该目标能谱的基准峰值为850,该参考峰值为800,则该参考峰值与该基准峰值之间的第一差值为-50,则该探测器的漂移值为-50,则说明该探测器的峰值向左漂移了50道址。
步骤208:控制终端根据该漂移值,对该探测器进行调整。
控制终端根据该漂移值,对该探测器进行补偿调整。例如,当该漂移值用于指示探测器向左漂移了该漂移值时,则控制终端将该探测器向右调整该漂移值。当该漂移值用于指示该探测器向右漂移了该漂移值时,则控制终端将该探测器向左调整该漂移值。
在发明实施例中,寻找峰值对应的第一道址,然后基于第一道址确定目标区域,从目标区域中确定出该目标能谱的参考峰值,从而提高了寻峰准确性。
本发明实施例提供了一种能谱寻峰装置,该装置应用在控制终端中,用于执行上述能谱寻峰方法中的控制终端执行的步骤,参见图6,该装置包括:
探测模块601,用于通过探测器探测放射源,得到放射源的目标能谱,目标能谱包括多个能谱点,每个能谱点用于指示道址和计数的对应关系;
第一确定模块602,用于从目标能谱中确定目标能谱的最大计数对应的第一道址;
第二确定模块603,用于根据第一道址,从目标能谱中确定包括第一道址的目标区域;
第三确定模块604,用于从目标区域中,确定目标能谱的第二道址;
第四确定模块605,用于根据第二道址,确定目标能谱的参考峰值。
在一个可能的实现方式中,第二确定模块603,包括:
寻找单元,用于从第一道址开始向两侧进行积分,寻找积分值不低于第一预设阈值的第一能谱点和第二能谱点;
组成单元,用于将目标能谱中第一能谱点和第二能谱点之间的能谱组成目标区域。
在一个可能的实现方式中,第三确定模块604,还用于确定目标区域的重心,将重心对应的道址确定为目标能谱的第二道址。
在一个可能的实现方式中,装置还包括:
调整模块,用于获取目标能谱的基准峰值,根据基准峰值和参考峰值,确定探测器的漂移值,根据漂移值,对探测器进行调整。
在一个可能的实现方式中,第二道址包括至少一个第二子道址,
第四确定模块605,还用于当第二道址包括一个第二子道址时,将第二子道址作为目标能谱的参考峰值;当第二道址包括多个第二子道址时,根据每个第二子道址,通过预设迭代算法进行迭代运算,输出目标能谱的参考峰值。
在一个可能的实现方式中,参见图7,该装置还包括:
滤波模块606,用于通过数字滤波器,将目标能谱中的高频成分过滤掉。
在一个可能的实现方式中,放射源为探测器自身具有衰变功能的晶体。
在发明实施例中,寻找峰值对应的第一道址,然后基于第一道址确定目标区域,从目标区域中确定出该目标能谱的参考峰值,从而提高了寻峰准确性。
需要说明的是:上述实施例提供的能谱寻峰装置在能谱寻峰时,仅以上述各功能模块的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的功能模块完成,即将装置的内部结构划分成不同的功能模块,以完成以上描述的全部或者部分功能。另外,上述实施例提供的能谱寻峰装置与能谱寻峰方法实施例属于同一构思,其具体实现过程详见方法实施例,这里不再赘述。
图8是本发明实施例提供的一种控制终端的结构示意图。该控制终端可以用于实施上述实施例所示出的获取网络接入信息的方法中的控制终端所执行的功能。具体来讲:
控制终端800可以包括RF(Radio Frequency,射频)电路810、包括有一个或一个以上计算机可读存储介质的存储器820、输入单元830、显示单元840、传感器850、音频电路860、传输模块870、包括有一个或者一个以上处理核心的处理器890、以及电源880等部件。本领域技术人员可以理解,图8中示出的控制终端结构并不构成对控制终端的限定,可以包括比图示更多或更少的部件,或者组合某些部件,或者不同的部件布置。其中:
RF电路810可用于收发信息或通话过程中,信号的接收和发送,特别地,将基站的下行信息接收后,交由一个或者一个以上处理器890处理;另外,将涉及上行的数据发送给基站。通常,RF电路810包括但不限于天线、至少一个放大器、调谐器、一个或多个振荡器、用户身份模块(SIM)卡、收发信机、耦合器、LNA(Low Noise Amplifier,低噪声放大器)、双工器等。此外,RF电路810还可以通过无线通信与网络和其他控制终端通信。所述无线通信可以使用任一通信标准或协议,包括但不限于GSM(Global System of Mobilecommunication,全球移动通讯系统)、GPRS(General Packet Radio Service,通用分组无线服务)、CDMA(Code Division Multiple Access,码分多址)、WCDMA(Wideband CodeDivision Multiple Access,宽带码分多址)、LTE(Long Term Evolution,长期演进)、电子邮件、SMS(Short Messaging Service,短消息服务)等。
存储器820可用于存储软件程序以及模块,如上述示例性实施例所示出的控制终端所对应的软件程序以及模块,处理器890通过运行存储在存储器820的软件程序以及模块,从而执行各种功能应用以及数据处理,如实现基于视频的交互等。存储器820可主要包括存储程序区和存储数据区,其中,存储程序区可存储操作系统、至少一个功能所需的应用程序(比如声音播放功能、图像播放功能等)等;存储数据区可存储根据控制终端800的使用所创建的数据(比如音频数据、电话本等)等。此外,存储器820可以包括高速随机存取存储器,还可以包括非易失性存储器,例如至少一个磁盘存储器件、闪存器件、或其他易失性固态存储器件。相应地,存储器820还可以包括存储器控制器,以提供处理器890和输入单元830对存储器820的访问。
输入单元830可用于接收输入的数字或字符信息,以及产生与用户设置以及功能控制有关的键盘、鼠标、操作杆、光学或者轨迹球信号输入。具体地,输入单元830可包括触敏表面831以及其他输入控制终端832。触敏表面831,也称为触摸显示屏或者触控板,可收集用户在其上或附近的触摸操作(比如用户使用手指、触笔等任何适合的物体或附件在触敏表面831上或在触敏表面831附近的操作),并根据预先设定的程式驱动相应的链接装置。可选的,触敏表面831可包括触摸检测装置和触摸控制器两个部分。其中,触摸检测装置检测用户的触摸方位,并检测触摸操作带来的信号,将信号传送给触摸控制器;触摸控制器从触摸检测装置上接收触摸信息,并将它转换成触点坐标,再送给处理器890,并能接收处理器890发来的命令并加以执行。此外,可以采用电阻式、电容式、红外线以及表面声波等多种类型实现触敏表面831。除了触敏表面831,输入单元830还可以包括其他输入控制终端832。具体地,其他输入控制终端832可以包括但不限于物理键盘、功能键(比如音量控制按键、开关按键等)、轨迹球、鼠标、操作杆等中的一种或多种。
显示单元840可用于显示由用户输入的信息或提供给用户的信息以及控制终端800的各种图形用户接口,这些图形用户接口可以由图形、文本、图标、视频和其任意组合来构成。显示单元840可包括显示面板841,可选的,可以采用LCD(Liquid Crystal Display,液晶显示器)、OLED(Organic Light-Emitting Diode,有机发光二极管)等形式来配置显示面板841。进一步的,触敏表面831可覆盖显示面板841,当触敏表面831检测到在其上或附近的触摸操作后,传送给处理器890以确定触摸事件的类型,随后处理器890根据触摸事件的类型在显示面板841上提供相应的视觉输出。虽然在图8中,触敏表面831与显示面板841是作为两个独立的部件来实现输入和输入功能,但是在某些实施例中,可以将触敏表面831与显示面板841集成而实现输入和输出功能。
控制终端800还可包括至少一种传感器850,比如光传感器、运动传感器以及其他传感器。具体地,光传感器可包括环境光传感器及接近传感器,其中,环境光传感器可根据环境光线的明暗来调节显示面板841的亮度,接近传感器可在控制终端800移动到耳边时,关闭显示面板841和/或背光。作为运动传感器的一种,重力加速度传感器可检测各个方向上(一般为三轴)加速度的大小,静止时可检测出重力的大小及方向,可用于识别手机姿态的应用(比如横竖屏切换、相关游戏、磁力计姿态校准)、振动识别相关功能(比如计步器、敲击)等;至于控制终端800还可配置的陀螺仪、气压计、湿度计、温度计、红外线传感器等其他传感器,在此不再赘述。
音频电路860、扬声器861,传声器862可提供用户与控制终端800之间的音频接口。音频电路860可将接收到的音频数据转换后的电信号,传输到扬声器861,由扬声器861转换为声音信号输出;另一方面,传声器862将收集的声音信号转换为电信号,由音频电路860接收后转换为音频数据,再将音频数据输出处理器890处理后,经RF电路810以发送给比如另一控制终端,或者将音频数据输出至存储器820以便进一步处理。音频电路860还可能包括耳塞插孔,以提供外设耳机与控制终端800的通信。
控制终端800通过传输模块870可以帮助用户收发电子邮件、浏览网页和访问流式媒体等,它为用户提供了无线或有线的宽带互联网访问。虽然图8示出了传输模块870,但是可以理解的是,其并不属于控制终端800的必须构成,完全可以根据需要在不改变发明的本质的范围内而省略。
处理器890是控制终端800的控制中心,利用各种接口和线路链接整个手机的各个部分,通过运行或执行存储在存储器820内的软件程序和/或模块,以及调用存储在存储器820内的数据,执行控制终端800的各种功能和处理数据,从而对手机进行整体监控。可选的,处理器890可包括一个或多个处理核心;优选的,处理器890可集成应用处理器和调制解调处理器,其中,应用处理器主要处理操作系统、用户界面和应用程序等,调制解调处理器主要处理无线通信。可以理解的是,上述调制解调处理器也可以不集成到处理器890中。
控制终端800还包括给各个部件供电的电源880(比如电池),优选的,电源可以通过电源管理系统与处理器890逻辑相连,从而通过电源管理系统实现管理充电、放电、以及功耗管理等功能。电源880还可以包括一个或一个以上的直流或交流电源、再充电系统、电源故障检测电路、电源转换器或者逆变器、电源状态指示器等任意组件。
尽管未示出,控制终端800还可以包括摄像头、蓝牙模块等,在此不再赘述。具体在本实施例中,控制终端800的显示单元是触摸屏显示器,控制终端800还包括有存储器以及至少一条指令、至少一段程序、代码集或指令集,其中至少一条指令、至少一段程序、代码集或指令集存储于存储器中,且经配置以由一个或者一个以上处理器加载并执行,以实现上述实施例中的能谱寻峰方法中所执行的操作。
本发明实施例还提供了一种计算机可读存储介质,该计算机可读存储介质应用于终端,该计算机可读存储介质中存储有至少一条指令、至少一段程序、代码集或指令集,该指令、该程序、该代码集或该指令集由处理器加载并执行以实现上述实施例的能谱寻峰方法中控制终端所执行的操作。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成,也可以通过程序来指令相关的硬件完成,所述的程序可以存储于一种计算机可读存储介质中,上述提到的存储介质可以是只读存储器,磁盘或光盘等。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (16)
1.一种能谱寻峰方法,其特征在于,所述方法包括:
通过探测器探测放射源,得到所述放射源的目标能谱,所述目标能谱包括多个能谱点,每个能谱点用于指示道址和计数的对应关系;
从所述目标能谱中确定所述目标能谱的最大计数对应的第一道址;
根据所述第一道址,从所述目标能谱中确定包括所述第一道址的目标区域;
从所述目标区域中,确定所述目标能谱的第二道址;
根据所述第二道址,确定所述目标能谱的参考峰值。
2.根据权利要求1所述的方法,其特征在于,所述根据所述第一道址,从所述目标能谱中确定包括所述第一道址的目标区域,包括:
从所述第一道址开始向两侧进行积分,寻找积分值不低于第一预设阈值的第一能谱点和第二能谱点;
将所述目标能谱中所述第一能谱点和所述第二能谱点之间的能谱组成所述目标区域。
3.根据权利要求1所述的方法,其特征在于,所述从所述目标区域中,确定所述目标能谱的第二道址,包括:
确定所述目标区域的重心;
将所述重心对应的道址确定为所述目标能谱的第二道址。
4.根据权利要求1所述的方法,其特征在于,所述根据所述第二道址,确定所述目标能谱的参考峰值之后,所述方法还包括:
获取所述目标能谱的基准峰值;
根据所述基准峰值和所述参考峰值,确定所述探测器的漂移值;
根据所述漂移值,对所述探测器进行调整。
5.根据权利要求1所述的方法,其特征在于,所述第二道址包括至少一个第二子道址,所述根据所述第二道址,确定所述目标能谱的参考峰值,包括:
当所述第二道址包括一个第二子道址时,将所述第二子道址作为所述目标能谱的参考峰值;
当所述第二道址包括多个第二子道址时,根据每个第二子道址,通过预设迭代算法进行迭代运算,输出所述目标能谱的参考峰值。
6.根据权利要求1所述的方法,其特征在于,所述从所述目标能谱中确定所述目标能谱的最大计数对应的第一道址之前,所述方法还包括:
通过数字滤波器,将所述目标能谱中的高频成分过滤掉。
7.根据权利要求1-6任一所述的方法,其特征在于,所述放射源为所述探测器自身具有衰变功能的晶体。
8.一种能谱寻峰装置,其特征在于,所述装置包括:
探测模块,用于通过探测器探测放射源,得到所述放射源的目标能谱,所述目标能谱包括多个能谱点,每个能谱点用于指示道址和计数的对应关系;
第一确定模块,用于从所述目标能谱中确定所述目标能谱的最大计数对应的第一道址;
第二确定模块,用于根据所述第一道址,从所述目标能谱中确定包括所述第一道址的目标区域;
第三确定模块,用于从所述目标区域中,确定所述目标能谱的第二道址;
第四确定模块,用于根据所述第二道址,确定所述目标能谱的参考峰值。
9.根据权利要求8所述的装置,其特征在于,所述第二确定模块,包括:
寻找单元,用于从所述第一道址开始向两侧进行积分,寻找积分值不低于第一预设阈值的第一能谱点和第二能谱点;
组成单元,用于将所述目标能谱中所述第一能谱点和所述第二能谱点之间的能谱组成所述目标区域。
10.根据权利要求8所述的装置,其特征在于,
所述第三确定模块,还用于确定所述目标区域的重心,将所述重心对应的道址确定为所述目标能谱的第二道址。
11.根据权利要求8所述的装置,其特征在于,所述装置还包括:
调整模块,用于获取所述目标能谱的基准峰值,根据所述基准峰值和所述参考峰值,确定所述探测器的漂移值,根据所述漂移值,对所述探测器进行调整。
12.根据权利要求8所述的装置,其特征在于,所述第二道址包括至少一个第二子道址,
所述第四确定模块,还用于当所述第二道址包括一个第二子道址时,将所述第二子道址作为所述目标能谱的参考峰值;当所述第二道址包括多个第二子道址时,根据每个第二子道址,通过预设迭代算法进行迭代运算,输出所述目标能谱的参考峰值。
13.根据权利要求8所述的装置,其特征在于,所述装置还包括:
滤波模块,用于通过数字滤波器,将所述目标能谱中的高频成分过滤掉。
14.根据权利要求8-13任一所述的装置,其特征在于,所述放射源为所述探测器自身具有衰变功能的晶体。
15.一种能谱寻峰装置,其特征在于,所述装置包括:
处理器;
用于存储处理器可执行指令的存储器;
其中,所述处理器被配置为:
通过探测器探测放射源,得到所述放射源的目标能谱,所述目标能谱包括多个能谱点,每个能谱点用于指示道址和计数的对应关系;
从所述目标能谱中确定所述目标能谱的最大计数对应的第一道址;
根据所述第一道址,从所述目标能谱中确定包括所述第一道址的目标区域;
从所述目标区域中,确定所述目标能谱的第二道址;
根据所述第二道址,确定所述目标能谱的参考峰值。
16.一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序,其特征在于,所述程序被处理器执行时实现以下步骤:
通过探测器探测放射源,得到所述放射源的目标能谱,所述目标能谱包括多个能谱点,每个能谱点用于指示道址和计数的对应关系;
从所述目标能谱中确定所述目标能谱的最大计数对应的第一道址;
根据所述第一道址,从所述目标能谱中确定包括所述第一道址的目标区域;
从所述目标区域中,确定所述目标能谱的第二道址;
根据所述第二道址,确定所述目标能谱的参考峰值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711363000.8A CN108324295B (zh) | 2017-12-18 | 2017-12-18 | 能谱寻峰方法、装置和计算机存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711363000.8A CN108324295B (zh) | 2017-12-18 | 2017-12-18 | 能谱寻峰方法、装置和计算机存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108324295A CN108324295A (zh) | 2018-07-27 |
CN108324295B true CN108324295B (zh) | 2021-03-23 |
Family
ID=62921840
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711363000.8A Active CN108324295B (zh) | 2017-12-18 | 2017-12-18 | 能谱寻峰方法、装置和计算机存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108324295B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111006873B (zh) * | 2019-12-05 | 2022-02-01 | 中国航发四川燃气涡轮研究院 | 叶尖间隙信号处理过程中的峰峰值获取方法和装置 |
CN112684490B (zh) * | 2020-12-22 | 2024-05-24 | 上海工物高技术产业发展有限公司 | 能谱分区方法、装置、计算机设备和存储介质 |
CN113917517B (zh) * | 2021-09-03 | 2023-12-05 | 上海工物高技术产业发展有限公司 | 能谱数据处理方法、装置、计算机设备和存储介质 |
CN115393303B (zh) * | 2022-08-17 | 2023-11-17 | 上海精积微半导体技术有限公司 | 周期性结构的周期提取方法及晶圆缺陷检测方法 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6635878B2 (en) * | 2001-06-13 | 2003-10-21 | Koninklijke Philips Electronics N.V. | Gamma camera with automatic adjustment of the energy spectrum |
US6624422B2 (en) * | 2001-09-25 | 2003-09-23 | Ge Medical Systems Global Technology Company, Llc | Method for dynamic stabilization of PET detector gains |
US7479639B1 (en) * | 2007-08-30 | 2009-01-20 | Orbotech Medical Solutions Ltd. | Apparatus, method and system for determining energy windows in medical-imaging systems |
CN101713829B (zh) * | 2009-12-22 | 2012-11-14 | 清华大学 | 一种射线能谱仪的快速稳谱方法和装置 |
US8907270B2 (en) * | 2010-06-30 | 2014-12-09 | Schlumberger Technology Corporation | Method and apparatus for gain regulation in a gamma detector |
CN102608649B (zh) * | 2012-03-02 | 2013-11-27 | 成都理工大学 | 统计分布式γ或X射线能谱解谱方法 |
US20150153461A1 (en) * | 2012-06-07 | 2015-06-04 | National University Corporation Shizuoka University | Radiation dosimeter and radiation dose calculation method |
CN104337531B (zh) * | 2013-07-25 | 2016-12-28 | 苏州瑞派宁科技有限公司 | 针对全数字pet系统的在线能量符合方法及系统 |
CN105030263B (zh) * | 2015-07-22 | 2018-06-12 | 湖北锐世数字医学影像科技有限公司 | 一种数字pet的能量反馈校正方法及系统 |
CN105115994A (zh) * | 2015-07-22 | 2015-12-02 | 武汉数字派特科技有限公司 | 一种数字pet能量参数化校正方法及系统 |
CN105182402B (zh) * | 2015-09-29 | 2018-08-03 | 沈阳东软医疗系统有限公司 | 一种闪烁晶体探测器增益的校正方法和装置 |
-
2017
- 2017-12-18 CN CN201711363000.8A patent/CN108324295B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108324295A (zh) | 2018-07-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108324295B (zh) | 能谱寻峰方法、装置和计算机存储介质 | |
CN106778585B (zh) | 一种人脸关键点跟踪方法和装置 | |
CN105933904B (zh) | 网络连接方法及装置 | |
KR102229699B1 (ko) | 정보 운용 방법 및 전자 장치 | |
US20170094652A1 (en) | Apparatus and method for configuring wifi channel of mobile device | |
CN107038681A (zh) | 图像虚化方法、装置、计算机可读存储介质和计算机设备 | |
KR20150108502A (ko) | 전자 장치 및 그 동작 방법 | |
CN105095432A (zh) | 网页批注显示方法及装置 | |
KR101848696B1 (ko) | 위치 정보를 콜라주에 중첩시키는 방법, 단말기 및 서버 | |
CN104360307A (zh) | 一种定位方法和装置 | |
CN107124556A (zh) | 对焦方法、装置、计算机可读存储介质和移动终端 | |
CN104133761A (zh) | 一种移动终端的内存占用分析方法、装置和系统 | |
CN108151716A (zh) | 飞行装置测绘作业区域规划方法、装置和终端 | |
CN105653112B (zh) | 一种显示浮层的方法及装置 | |
JP2022515481A (ja) | ネットワーク探索リスト生成方法、ネットワーク探索方法及び対応する装置 | |
CN110413639A (zh) | 数据核查方法及装置、电子设备和计算机可读存储介质 | |
EP2998843A1 (en) | Intelligent terminal and method for displaying input operation interface thereof | |
US10869293B2 (en) | Positioning cycle adjustment method and apparatus | |
CN107153715B (zh) | 在页面上添加文件的方法及装置 | |
CN104239005A (zh) | 图形对齐方法及装置 | |
CN109427101A (zh) | 一种获取增强现实图像的方法和系统 | |
CN107741839B (zh) | 一种基于文本阅读器的文本显示方法和装置 | |
CN104135464A (zh) | 一种联运中网页应用的帐号信息处理方法、装置和系统 | |
CN104134044B (zh) | 一种数据安全性的检测方法、装置和系统 | |
CN105486302A (zh) | 一种地图显示方法及终端 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |