CN104502948B - 辐射剂量空间分布测量用探测器矩阵相对响应关系的检测方法 - Google Patents
辐射剂量空间分布测量用探测器矩阵相对响应关系的检测方法 Download PDFInfo
- Publication number
- CN104502948B CN104502948B CN201410771536.3A CN201410771536A CN104502948B CN 104502948 B CN104502948 B CN 104502948B CN 201410771536 A CN201410771536 A CN 201410771536A CN 104502948 B CN104502948 B CN 104502948B
- Authority
- CN
- China
- Prior art keywords
- detecting element
- measurement
- radiation
- matrix
- detector
- 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.)
- Expired - Fee Related
Links
Landscapes
- Measurement Of Radiation (AREA)
Abstract
本发明涉及一种辐射剂量空间分布测量用探测器矩阵的相对响应关系检测方法,包括以下步骤:步骤(1)射束主轴沿线矩阵中心法线从正面穿过线矩阵中心,该位置作为初始位置,待出束条件稳定后,采集一组辐射场的测量数据,记录每个探测元件初始位置的测量数据;步骤(2)整体移动探测器矩阵,使得相邻两支探测器依次在辐射野的相同位置接受同样的辐射,并记录移动后每个探测元件的测量数据;步骤(3)经过步骤(1)和步骤(2)后,在每个探测元件对应的位置分别可以获得两个相邻探测元件的测量结果,采用递归算法,计算探测元件之间的相对响应关系。该方法简便易行,对于测量所用的辐射场无特殊要求,测量过程步骤少,评价数据可靠。
Description
技术领域
本发明涉及一种辐射剂量空间分布测量用探测器矩阵相对响应关系的检测方法,放射治疗中辐射剂量的空间分布测量和剂量验证,以及其它测量辐射剂量空间分布的领域。
背景技术
探测器矩阵是一种用于测量辐射剂量空间分布的设备,该设备的探测单元由多个探测器按照一定规则排列,组成一维或二维探测器阵列。目前还没有有效评价组成探测器矩阵的探测器相对响应关系的方法,解决探测单元之间响应差异不确定对空间剂量分布测量结果的影响。
发明内容
本发明所要解决的技术问题是提供一种辐射剂量空间分布测量用探测器矩阵相对响应关系的检测方法,克服现有探测器矩阵中探测元件之间可能存在响应差异不确定导致对空间剂量分布测量结果影响的缺陷。
本发明解决上述技术问题的技术方案如下:一种辐射剂量空间分布测量用探测器矩阵相对响应关系的检测方法,包括以下步骤:
步骤(1)射束主轴沿线矩阵中心法线从正面穿过线矩阵中心,该位置作为初始位置,待出束条件稳定后,采集一组辐射场的测量数据,记录每个探测元件初始位置的测量数据;
步骤(2)整体移动探测器矩阵,使得相邻两支探测元件依次在辐射野的相同位置接受同样的辐射,并记录移动后每个探测元件的测量数据;
步骤(3)经过步骤(1)和步骤(2)后,除左边第1个或第一列探测元件对应的初始位置和平移后右边最后1个或1列探测元件对应的位置外,在每个或每列探测元件对应的位置分别可以获得两个相邻探测元件的测量结果,采用递归算法,计算探测元件之间的相对响应关系。
在上述技术方案的基础上,本发明还可以做如下改进。
本发明如上所述一种辐射剂量空间分布测量用探测器矩阵相对响应关系的检测方法,所述探测器矩阵为一维排列的探测元件,步骤(2)整体移动探测器矩阵是将一维排列的探测元件整体向一侧或另一侧移动相邻两个探测元件的中心距离。
本发明如上所述一种辐射剂量空间分布测量用探测器矩阵相对响应关系的检测方法,所述探测器矩阵为二维排列的探测元件,步骤(2)所述整体移动探测器矩阵是在初始位置的基础上,先将二维排列的探测元件整体向一侧或另一侧移动相邻两个探测元件的中心距离,记录移动后每个探测元件的测量数据;再在初始位置的基础上,将二维排列的探测元件整体向上或向下移动相邻两个探测元件的中心距离,记录移动后每个探测元件的测量数据。
本发明的有益效果是:本方法简便易行,对于测量所用的辐射场无特殊要求,测量过程步骤少,数据处理无需复杂计算,评价数据可靠,适用一维和二维探测器矩阵。
附图说明
图1为本发明一种辐射剂量空间分布测量用探测器矩阵相对响应关系的检测方法其中一种一维矩阵的实施方式的布置示意图;
图2为本发明一种辐射剂量空间分布测量用探测器矩阵相对响应关系的检测方法其中一种二维矩阵的实施方式的布置示意图;
图3为本发明一种辐射剂量空间分布测量用探测器矩阵相对响应关系的检测方法的一种一维矩阵的每一个探测器响应的相对关系的计算过程示意图;
图4为本发明一种辐射剂量空间分布测量用探测器矩阵相对响应关系的检测方法对德国IBA公司生产的MatriXX Evolution型探测器矩阵的相对响应关系在Co-60γ射线辐射场中进行了评价探测器矩阵一致性的评价结果示意图。
附图中,各标号所代表的部件列表如下:
1、初始辐照位置,2、向右侧或X轴方向移动后的位置,3、射线束,4,y轴方向移动后的位置。
具体实施方式
以下结合附图对本发明的原理和特征进行描述,所举实例只用于解释本发明,并非用于限定本发明的范围。
如图1所示一维矩阵测量方法,一维矩阵响应相对响应关系的测量过程分为两步,其中方块代表一维矩阵初始位置,方块中的三角形代表一维矩阵平移后的位置;测量时具体的操作步骤是:第一步:将被测探测器矩阵放置于稳定的辐射场中,该位置作为初始位置1,待被测探测器矩阵充分预热稳定后,采集一组辐射场的测量数据,记录每个探测元件初始位置的测量数据;
第二步:以初始位置为起点,整体移动探测器矩阵,移动距离等于相邻探测器的间距,以向右移动为例,将左边第一个探测元件移至第一步中左边第二个探测元件的位置,如图1所示移动后的位置2,其它条件不变,再采集一组辐射场的测量数据。
组成一维矩阵的探测元件从左至右可以标记为d1、d2、……、di(i=1、2、…、n,n为探测元件数量),则假设在初始位置每个探测元件的测量结果分别为do1、do2、……、doi,向右平移一个探测元件间距后,每个探测元件的测量结果为ds1、ds2、……、dsi,至此,除左边第1个探测元件对应的初始位置和平移后右边最后1个探测元件对应的位置外,在每个探测元件对应的位置分别可以获得两个相邻探测元件的测量结果,以测量过程中辐射条件没有变化为前提,两个探测元件接受的辐射相同,位置与探测元件测量结果的对应关系见表1。
表1:一维矩阵测量结果与位置对应关系
注:*位置是指线矩阵在辐射野中的初始位置,用探测元件排列顺序表示。
ds1和do2分别是探测元件d1和d2在同一位置的测量结果,两个测量结果相除就可以获得探测元件d1和d2的相对响应关系r1,2,同理,经过n-1次计算就可以获得相邻两个探测元件的相对响应关系,即探测元件d1和d2、探测元件d2和d3、探测元件d3和d4,直至探测元件di-1和di的相对响应关系r1,2、r2,3、r3,4、……、ri-1,i。
利用探测元件d1和d2相对响应关系和探测元件d2和d3的相对响应关系,采用乘法可以获得探测元件d1和d3的相对响应关系r1,3;利用探测元件d1和d3的相对响应关系和探测元件d3和d4的相对响应关系,可以获得探测元件d1和d4的相对响应关系r1,4,以此类推,经过n-2次运算,就可以获得探测元件d1和线矩阵中其余探测元件的相对响应关系r1,i,至此就可以获得组成线矩阵的每一个探测器响应的相对关系,计算过程见图3。
图2为二维矩阵检测m×n(m和n分别是探测元件的行数和列数)的面矩阵相对响应关系的思路与一维矩阵相对响应关系的思路相同,测量过程分为三步(测量过程中的位置见图2):
第一步:将被测二维探测器矩阵放置于稳定的辐射场中,该位置作为初始位置1,待出射线束3条件稳定后,采集一组辐射场的测量数据;
第二步:将二维矩阵从初始位置向左或向右横向移动一个探测元件间距,以向右(x轴方向)移动为例,将左边第一列探测元件移至第一步中左边第二列探测元件的位置2,其它条件不变,再采集一组辐射场的测量数据;
第三步:将二维矩阵从初始位置向下或向上(y轴方向)移动一个探测元件间距,以向下移动为例,将上边第一行探测元件移至第一步中上边第二列探测元件的位置4,其它条件不变,再采集一组辐射场的测量数据。
二维矩阵的三个位置测量结束后,首先采用初始位置和按行方向平移一个探测元件间距的测量结果,将测量数据按行分解,可以获得两组m行一维的测量结果,根据一维矩阵的计算方法,可以获得每一行探测元件的响应一致性结果;再采用初始位置和按列方向平移一个探测元件间距的测量结果,将测量数据按列分解,可以获得两组n列一维的测量结果,同样根据线矩阵的计算方法,可以获得每一列探测元件的相对响应关系。由于已经获得每一行探测元件的响应一致性结果,实际只需要一列探测元件的响应一致性结果,就可以获得组成二维矩阵的所有探测元件的相对响应关系。
实施例1
方法对德国IBA公司生产的MatriXX Evolution型探测器矩阵的相对响应关系在Co-60γ射线辐射场中进行了评价,并用评价结果重新修正了测量数据,同时对测量位置的剂量分布采用另外的单支探测器逐点进行测量。结果如图4所示。
由图4可以看出,未经由本评价方法获得的数据修正的测量结果和采用单独探测器逐点测量的结果存在最大±5%的差异,而经本发明评价修正后的测量结果和采用单独探测器逐点测量的结果差异小于1%,由此证明本方法可靠,能有效解决探测器中探测元件之间响应差异对空间剂量分布测量结果的影响。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种辐射剂量空间分布测量用探测器矩阵的相对响应关系检测方法,其特征在于,包括以下步骤:
步骤(1)射束主轴沿线矩阵中心法线从正面穿过线矩阵中心,该位置作为初始位置,待出束条件稳定后,采集一组辐射场的测量数据,记录每个探测元件初始位置的测量数据;
步骤(2)整体移动探测器矩阵,使得相邻两支探测元件依次在辐射野的相同位置接受同样的辐射,并记录移动后每个探测元件的测量数据;
步骤(3),经过步骤(1)和步骤(2)后,除左边第1个或第一列探测元件对应的初始位置和平移后右边最后1个或1列探测元件对应的位置外,在每个或每列探测元件对应的位置分别可以获得两个相邻探测元件的测量结果,采用递归算法,计算探测元件之间的相对响应关系。
2.根据权利要求1所述一种辐射剂量空间分布测量用探测器矩阵的相对响应关系检测方法,其特征在于,所述探测器矩阵为一维排列的探测元件,步骤(2)整体移动探测器矩阵是将一维排列的探测元件整体向一侧或另一侧移动相邻两个探测元件的中心距离。
3.根据权利要求1所述一种辐射剂量空间分布测量用探测器矩阵的相对响应关系检测方法,其特征在于,所述探测器矩阵为二维排列的探测元件,步骤(2)所述整体移动探测器矩阵是在初始位置的基础上,先将二维排列的探测元件整体向一侧或另一侧移动相邻两个探测元件的中心距离,记录移动后每个探测元件的测量数据;再在初始位置的基础上将二维排列的探测元件整体向上或向下移动相邻两个探测元件的中心距离,记录移动后每个探测元件的测量数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410771536.3A CN104502948B (zh) | 2014-12-12 | 2014-12-12 | 辐射剂量空间分布测量用探测器矩阵相对响应关系的检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410771536.3A CN104502948B (zh) | 2014-12-12 | 2014-12-12 | 辐射剂量空间分布测量用探测器矩阵相对响应关系的检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104502948A CN104502948A (zh) | 2015-04-08 |
CN104502948B true CN104502948B (zh) | 2017-06-13 |
Family
ID=52944359
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410771536.3A Expired - Fee Related CN104502948B (zh) | 2014-12-12 | 2014-12-12 | 辐射剂量空间分布测量用探测器矩阵相对响应关系的检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104502948B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111538068B (zh) * | 2020-05-07 | 2022-04-22 | 中国计量科学研究院 | 一种放射性表面污染测量装置及表面污染测量方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6429431B1 (en) * | 1999-09-24 | 2002-08-06 | Peter J. Wilk | Medical diagnostic method and apparatus utilizing radioactivity detection |
CN102656478A (zh) * | 2009-12-15 | 2012-09-05 | 皇家飞利浦电子股份有限公司 | 基于辐射剂量的成像探测器瓦片参数补偿 |
CN102725659A (zh) * | 2009-12-01 | 2012-10-10 | 原子能和能源替代品委员会 | X射线探测器的校准方法 |
CN103188996A (zh) * | 2010-10-29 | 2013-07-03 | 富士胶片株式会社 | 放射线照相相衬成像设备 |
-
2014
- 2014-12-12 CN CN201410771536.3A patent/CN104502948B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6429431B1 (en) * | 1999-09-24 | 2002-08-06 | Peter J. Wilk | Medical diagnostic method and apparatus utilizing radioactivity detection |
CN102725659A (zh) * | 2009-12-01 | 2012-10-10 | 原子能和能源替代品委员会 | X射线探测器的校准方法 |
CN102656478A (zh) * | 2009-12-15 | 2012-09-05 | 皇家飞利浦电子股份有限公司 | 基于辐射剂量的成像探测器瓦片参数补偿 |
CN103188996A (zh) * | 2010-10-29 | 2013-07-03 | 富士胶片株式会社 | 放射线照相相衬成像设备 |
Also Published As
Publication number | Publication date |
---|---|
CN104502948A (zh) | 2015-04-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101297221B (zh) | 用于谱计算机断层摄影的方法和设备 | |
CN107430201B (zh) | 放射线拍摄装置、放射线计数装置以及放射线拍摄方法 | |
JP4968384B2 (ja) | 放射線断層撮影装置 | |
CN102331586A (zh) | 电荷丢失校正 | |
CN101577284B (zh) | 用于测量辐射的半导体探测器及成像装置 | |
CN104936525B (zh) | X射线计算机断层摄影装置以及光子计数方法 | |
JP5237919B2 (ja) | 核医学診断装置 | |
JP6282435B2 (ja) | ミュオン軌跡検出器及びミュオン軌跡検出方法 | |
CN102283665A (zh) | 核医学成像装置以及核医学成像方法 | |
JP6258359B2 (ja) | 遅延を挿入する装置、核医学イメージング装置、遅延を挿入する方法および較正する方法 | |
CN103890571B (zh) | 具有偏移校正的用于探测光子的放射探测装置 | |
WO2017000711A1 (zh) | 一种探测器信号读出的通道复用方法 | |
EP2573588B1 (en) | Positron ct apparatus and timing correction method | |
US6326624B1 (en) | Device and method for determining the assumed position of a phenomenon relative to a set of photodetectors, and application to gamma-cameras | |
TWI356689B (en) | Method for determining location of gamma interacti | |
US20100084560A1 (en) | Timing response improvement in light-sharing detectors | |
CN104502948B (zh) | 辐射剂量空间分布测量用探测器矩阵相对响应关系的检测方法 | |
GB2528168A (en) | Analysis of signals from pixellated detectors of ionizing radiation | |
CN110012673A (zh) | 用于光谱分析和增益调整的系统和方法 | |
CN109521459B (zh) | 一种射线在闪烁晶体中的击中点定位方法及其系统 | |
TWI400470B (zh) | 決定晶體穿透效應與偵檢效率以及偵檢效率校正之方法 | |
JP6737154B2 (ja) | 放射線検出装置 | |
JP2012177555A (ja) | 放射線撮像装置 | |
TW201714583A (zh) | 加馬入射事件三維發生位置與能量之辨識方法及其裝置 | |
CN111624645B (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 |
Granted publication date: 20170613 Termination date: 20171212 |
|
CF01 | Termination of patent right due to non-payment of annual fee |