CN101797423A - 自校正、多曲线拟合数字式剂量监控方法和系统 - Google Patents
自校正、多曲线拟合数字式剂量监控方法和系统 Download PDFInfo
- Publication number
- CN101797423A CN101797423A CN 201010143793 CN201010143793A CN101797423A CN 101797423 A CN101797423 A CN 101797423A CN 201010143793 CN201010143793 CN 201010143793 CN 201010143793 A CN201010143793 A CN 201010143793A CN 101797423 A CN101797423 A CN 101797423A
- Authority
- CN
- China
- Prior art keywords
- ionization chamber
- signal
- output signal
- current
- fitting function
- 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
Links
Images
Landscapes
- Measurement Of Radiation (AREA)
Abstract
本发明公开了一种自校正、多曲线拟合数字式剂量监控方法和系统,所述方法包括:预先根据采集值进行数据拟合,拟合函数为分段函数,从而建立起电离室输出信号与实际剂量率的实际非线性对应关系,并非简单的线性关系,因此本发明的监控系统能够获得精确的实际剂量率,即便在剂量率有较大变化的情况下,也不会产生监控到的剂量与实际剂量之间的大偏差。
Description
技术领域
本发明涉及放射治疗器械的技术领域,具体涉及自校正、多曲线拟合数字式剂量监控方法和系统。
背景技术
剂量监控系统是将医用电子加速器的X辐射或电子辐射应用于放射治疗的基础之一。安全、准确、可靠地传送所需的剂量正是剂量监控系统的主要功能。
在放射治疗过程中,剂量监控系统采用透射型电离室采集电子加速器的实际剂量率,并发送给监控主机,由监控主机根据实际剂量率调整后续给予剂量率的大小,以保证向病人安全、准确、可靠地传送所需剂量。
但是,电离室输出信号与实际剂量率并非简单的线性关系,如果不对电离室的输出信号进行校正,那么在剂量率有较大变化的情况下,监控剂量与实际剂量就会有较大的偏差。在这种情况下,监控主机根据电离室采集的实际剂量率信号进行剂量率调节不够精确,从而导致给予的剂量不符合病人需求,危害病人健康。
发明内容
本发明提供了一种自校正、多曲线拟合数字式剂量监控方法,该方法能够在剂量率有较大变化的情况下,减小监控到的剂量与实际剂量的偏差。
该方法包括:
预先采集不同剂量率射束照射下对应的电离室输出信号的大小,将采集的电离室输出信号进行分段,根据采集数据进行分段拟合,得到各电离室输出信号分段的拟合函数;
记录电离室输出信号分段与拟合函数之间的对应关系;
在实际监控过程中,采集当前电离室的输出信号,查找建立的对应关系,查找到当前电离室输出信号所在的分段,以及该所在分段对应的拟合函数;采用得到的拟合函数计算当前电离室的输出信号对应的当前剂量率;再采用当前工作温度和大气压力对当前剂量率进行再次校正,得到精确的当前剂量率。
为了减少存储量,建立对应关系时,记录电离室输出信号与拟合函数的参数之间的对应关系,约定各分段的函数类型,例如是一次函数还是二次函数,那么在根据对应关系获取分段函数时,根据当前电离室的输出信号查找所述对应关系,获得当前电离室输出信号所在分段以及该所在分段对应的拟合函数的参数,再根据获得的拟合函数的参数以及预先约定的各分段的函数类型,组成拟合函数,将当前电离室的输出信号值代入拟合函数得到当前剂量率。
本发明还提供了一种自校正、多曲线拟合数字式剂量监控系统,该系统包括两路监控通道,每一路监控通道都包括电离室、前置放大板和剂量监控板;
对于每一路监控通道:
电离室被分为两部分,输出第一信号和第二信号;其中,其中一个电离室被分为左右两部分,输出左信号和右信号,另一个电离室被分为前后两部分,输出前信号和后信号;
在电离室安装的温度传感器输出电离室的温度信号、气压传感器输出电离室的气压信号;第一信号、第二信号、温度信号和气压信号经过前置放大板的放大输入给剂量监控板;
剂量监控板包括箝位模块、脉冲积分模块、模数转换模块、同步脉冲差分模块以及单片机;所述同步脉冲差分模块将同步脉冲信号进行差分处理后,送入脉冲积分模块和单片机;所述第一信号和第二信号分别通过箝位模块的箝位处理和脉冲积分模块的脉冲积分处理输入模数转换模块,温度信号和气压信号输入模数转换模块;
模数转换模块对第一信号、第二信号、温度信号和气压信号进行模数转换输入单片机;
单片机记录有电离室输出信号分段与拟合函数之间的对应关系,该对应关系的确定方式为:预先采集不同剂量率射束照射下对应的电离室输出信号的大小,将采集的电离室输出信号进行分段,根据采集数据进行分段拟合,得到各电离室输出信号分段的拟合函数,之后建立电离室输出信号分段与拟合函数之间的对应关系;
在实际监控过程中,所述单片机接收来自模数转换模块的第一信号、第二信号、温度信号和气压信号,将第一信号和第二信号相加,根据相加结果在自身记录的所述对应关系中查找到当前电离室输出信号所在分段,以及该所在分段对应的拟合函数;采用得到的拟合函数计算电离室的当前输出信号对应的当前剂量率,再采用电离室的当前工作温度和大气压力对当前剂量率进行再次校正,得到精确的当前剂量率。
由以上方案可以看出,采用该方法预先根据采集值进行数据拟合,拟合函数为分段函数,从而建立起电离室输出信号与实际剂量率的实际非线性对应关系,并非简单的线性关系,因此本发明的监控系统能够获得精确的实际剂量率,即便在剂量率有较大变化的情况下,也不会产生监控到的剂量与实际剂量之间的大偏差。而且本发明采用了单片机,采用数字器件实现根据对应关系的自校正,系统搭建方便,以较小的代价实现了减小监控到的剂量与实际剂量之间偏差的效果。
附图说明
图1为本发明自校正、多曲线拟合数字式剂量监控方法的流程图。
图2为本发明自校正、多曲线拟合数字式剂量监控系统的结构示意图。
图3为图2中剂量控制板的结构示意图。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
本发明提供了一种自校正、多曲线拟合数字式剂量监控方案,其基本思想为:预先根据采集值进行数据拟合,拟合函数为分段函数,从而建立起电离室输出信号与实际剂量率的实际非线性对应关系,并非简单的线性关系,因此本发明的监控系统能够获得精确的实际剂量率,即便在剂量率有较大变化的情况下,也不会产生监控到的剂量与实际剂量之间的大偏差。
图1为自校正、多曲线拟合数字式剂量监控方法的流程图。如图1所示,该流程包括:
步骤1:数据采集。
本步骤中,将电离室安装在直线加速器样机上,用不同剂量率的射束透过电离室,记录在各个剂量率时电离室的输出信号的大小,将采集的电离室输出信号进行分段,根据这些采集数值进行分段拟合,得到各电离室输出信号分段的拟合函数。
具体来说,将剂量率[m,n]区域分为N个子区域,记为[mi,ni](i=1、2…N),针对每个子区域进行拟合,得到子区域1的拟合函数y1=k1x+b1,y2=k2x2+b2x+c2,等等。其中,y为电离室输出信号值,x为实际剂量率。
例如,在250-360伦琴/分钟的区域,x=100时,y=251;x=145时,y=342,把这两组数据分别代入250-360伦琴/分钟的区域对应的拟合函数,求出k≈2.02,b≈24.12。
步骤2、建立数据库。
本步骤中,记录电离室输出信号分段与拟合函数的参数之间的对应关系,且记录约定各分段的函数类型。在实际中,也可以直接记录电离室输出信号分段与拟合函数之间的对应关系。
具体来说,对于[mi,ni],i=1时在数据库记录[m1,n1]与k1、b1的对应关系,i=2时,在数据库记录[m2,n2]与k2、b2、c2的对应关系。表1示出了一种实际数据库的内容。
表1
步骤3:剂量监控过程。
该步骤3包括以下子步骤:
步骤31、采集当前电离室的输出信号,查找记录在数据库中的对应关系,查找到当前电离室输出信号所在的分段,以及该所在分段对应的拟合函数的参数,从而得到相应的拟合函数。例如,采集到的电离室输出信号y为300,查上表,得到它所在分段的参数k=2.02,b=24.12,则拟合函数为
步骤32、采用得到的拟合函数计算当前电离室的输出信号对应的当前剂量率。例如,将y=300、k=2.02、b=24.12代入步骤31得到的拟合函数,即可算出当前电离室的输出信号对应的当前剂量率为124.4伦琴/分钟。
步骤33、再采用电离室的当前工作温度和大气压力对当前剂量率进行校正,得到精确的当前剂量率。本步骤的校正可以采用常规的根据温度和压力对剂量率进行校正的方法实现。
例如,将采用温度传感器获取的当前工作温度和采用大气压力传感器获取的当前大气压力信号经A/D转换后送入系统中的计算模块,采用如下公式对当前剂量率进行校正,得到精确的当前剂量率。假设当前的工作温度t为25℃、大气压力P为1.000×105Pa,代入下式即可算出当前的校正因子J≈1.03;其中,P0=1.013×105Pa。
最后采用校正因子J对当前剂量率进行校正,得到准确剂量率为124.4×1.03≈128.1伦琴/分钟。
至此,本流程结束。
为了实现本发明的自校正、多曲线拟合数字式剂量监控方法,本发明还提供了一种如图2所示的自校正、多曲线拟合数字式剂量监控系统。该数字式剂量监控系统包括两组监控通道,每一路监控通道都包括安装在直线加速器上的电离室、前置放大板和剂量监控板,各路监控通道的处理过程相通。
下面针对一路监控通道进行描述。
如图2所示,电离室安装在直线加速器上,且被分为两部分,分别输出第一信号和第二信号;其中,电离室1被分为左右两部分,输出左(L)信号和右(R)信号,电离室2被分为前后两部分,输出前(G)信号和后(T)信号。
并且在电离室附近安装的温度传感器输出温度信号,在电离室附近安装的气压传感器输出气压信号,第一信号、第二信号、温度信号和气压信号通过前置放大板的放大,输出给剂量监控板。
如图3所示,剂量监控板包括箝位模块、脉冲积分模块、模数(AD)转换模块、同步脉冲差分模块以及单片机(例如89C52)。其中,
同步脉冲差分模块将同步脉冲信号进行差分处理后,送入脉冲积分模块和单片机。
第一信号和第二信号分别通过箝位模块的箝位处理和脉冲积分模块的脉冲积分处理后输入AD转换模块,温度信号和气压信号直接输入AD转换模块。
AD转换模块对第一信号、第二信号、温度信号和气压信号进行模数转换输入单片机。
单片机中记录有电离室输出信号分段与拟合函数之间的对应关系,该对应关系的确定方式为:预先采集不同剂量率射束照射下对应的电离室输出信号的大小,将采集的电离室输出信号进行分段,根据采集数据进行分段拟合,得到各电离室输出信号分段的拟合函数,之后建立电离室输出信号分段与拟合函数之间的对应关系。其中,记录的对应关系可以具体为:电离室输出信号分段与拟合函数的参数之间的对应关系,并约定分段与函数类型之间的关系。
在实际监控过程中,所述单片机接收来自模数转换模块的第一信号、第二信号、温度信号和气压信号,将第一信号和第二信号相加,根据相加结果在自身记录的所述对应关系中查找到当前电离室输出信号所在分段,以及该所在分段对应的拟合函数;采用得到的拟合函数计算电离室的当前输出信号对应的当前剂量率,再采用电离室的当前工作温度和大气压力对当前剂量率进行再次校正,得到精确的当前剂量率。
如图2所示,剂量监控板还根据第一信号和第二信号进行比较,得到第一信号和第二信号之间的误差,输出到相应的输出导向电源。例如,L和R之间的误差输入到X输出导向电源,X输出导向电源根据LR之间的误差控制束流照射电离室的位置,用以控制电离室1的左右对称性。同理,T和G信号之间的误差输入到Y输出导向电源,Y输出导向电源根据TG之间的误差控制束流照射电离室的位置,用以控制电离室2的左右对称性。
此外,剂量控制板进一步向前置放大板发送能档控制信号(能档控制1和能档控制0),用于控制前置放大的倍数,不同能档对应不同的放大倍数。剂量监控板处理后的剂量率值通过485总线输入给后级的控制台,控制台上的工控机根据剂量率进行控制处理。由于剂量率是根据多段拟合的函数计算得到,因此工控机可以根据精确的剂量率进行控制,从而避免了控制上的大偏差。
Claims (4)
1.一种自校正、多曲线拟合数字式剂量监控方法,其特征在于,该方法包括:
预先采集不同剂量率射束照射下对应的电离室输出信号的大小,将采集的电离室输出信号进行分段,根据采集数据进行分段拟合,得到各电离室输出信号分段的拟合函数;
记录电离室输出信号分段与拟合函数之间的对应关系;
在实际监控过程中,采集当前电离室的输出信号,根据所述对应关系,查找当前电离室输出信号所在分段,以及该所在分段对应的拟合函数;采用得到的拟合函数计算当前电离室的输出信号对应的当前剂量率。
2.如权利要求1所述的方法,其特征在于,所述采用得到的拟合函数计算当前电离室的输出信号对应的当前剂量率之后,该方法进一步包括:采用电离室的当前工作温度和大气压力对当前剂量率进行再次校正,得到精确的当前剂量率。
3.如权利要求1所述的方法,其特征在于,所述记录电离室输出信号分段与拟合函数之间的对应关系为:记录电离室输出信号与拟合函数的参数之间的对应关系,约定各分段的函数类型;
所述查找所述对应关系,查找到当前电离室输出信号所在分段,以及该所在分段对应的拟合函数为:根据当前电离室的输出信号查找所述对应关系,获得当前电离室输出信号所在分段以及该所在分段对应的拟合函数的参数,再根据获得的拟合函数的参数以及预先约定的各分段的函数类型,组成拟合函数,将当前电离室的输出信号值代入拟合函数得到当前剂量率。
4.一种自校正、多曲线拟合数字式剂量监控系统,其特征在于,该系统包括两路监控通道,每一路监控通道都包括电离室、前置放大板和剂量监控板;
对于每一路监控通道:
电离室被分为两部分,分别输出第一信号和第二信号;其中,其中一个电离室被分为左右两部分,分别输出左信号和右信号,另一个电离室被分为前后两部分,输出前信号和后信号;
在电离室附近安装的温度传感器输出电离室的温度信号、气压传感器输出电离室的气压信号;第一信号、第二信号、温度信号和气压信号经过前置放大板的放大输入给剂量监控板;
剂量监控板包括箝位模块、脉冲积分模块、模数转换模块、同步脉冲差分模块以及单片机;所述同步脉冲差分模块将同步脉冲信号进行差分处理后,送入脉冲积分模块和单片机;所述第一信号和第二信号分别通过箝位模块的箝位处理和脉冲积分模块的脉冲积分处理输入模数转换模块,温度信号和气压信号输入模数转换模块;
模数转换模块对第一信号、第二信号、温度信号和气压信号进行模数转换输入单片机;
单片机记录有电离室输出信号分段与拟合函数之间的对应关系,该对应关系的确定方式为:预先采集不同剂量率射束照射下对应的电离室输出信号的大小,将采集的电离室输出信号进行分段,根据采集数据进行分段拟合,得到各电离室输出信号分段的拟合函数,之后建立电离室输出信号分段与拟合函数之间的对应关系;
在实际监控过程中,所述单片机接收来自模数转换模块的第一信号、第二信号、温度信号和气压信号,将第一信号和第二信号相加,根据相加结果在自身记录的所述对应关系中查找到当前电离室输出信号所在分段,以及该所在分段对应的拟合函数;采用得到的拟合函数计算电离室的当前输出信号对应的当前剂量率,再采用电离室的当前工作温度和大气压力对当前剂量率进行再次校正,得到精确的当前剂量率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101437934A CN101797423B (zh) | 2010-04-12 | 2010-04-12 | 自校正、多曲线拟合数字式剂量监控方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101437934A CN101797423B (zh) | 2010-04-12 | 2010-04-12 | 自校正、多曲线拟合数字式剂量监控方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101797423A true CN101797423A (zh) | 2010-08-11 |
CN101797423B CN101797423B (zh) | 2011-06-15 |
Family
ID=42593360
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010101437934A Active CN101797423B (zh) | 2010-04-12 | 2010-04-12 | 自校正、多曲线拟合数字式剂量监控方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101797423B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103175547A (zh) * | 2011-12-21 | 2013-06-26 | 北京普源精电科技有限公司 | 一种数据采集装置的参数拟合方法 |
CN109603016A (zh) * | 2018-12-14 | 2019-04-12 | 江苏海明医疗器械有限公司 | 一种医用加速器剂量率稳定控制系统 |
CN112082627A (zh) * | 2020-08-24 | 2020-12-15 | 浙江威星智能仪表股份有限公司 | 一种基于光信号采样的软件计量校准处理方法 |
CN112232381A (zh) * | 2020-09-25 | 2021-01-15 | 国网上海市电力公司 | 主导负荷参数类噪声辨识的模型参数后处理方法和装置 |
CN113032732A (zh) * | 2021-03-24 | 2021-06-25 | 中国兵器装备集团自动化研究所有限公司 | 一种基于相对误差分段的剂量率拟合方法及系统 |
CN115779281A (zh) * | 2022-10-28 | 2023-03-14 | 中子高新技术产业发展(重庆)有限公司 | 一种混合辐射场的测量装置和方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060273746A1 (en) * | 2005-05-19 | 2006-12-07 | Marc Delaperriere | Method for function monitoring in medical accelerator systems |
CN101290354A (zh) * | 2008-06-12 | 2008-10-22 | 中国测试技术研究院电离辐射研究所 | 用于放射治疗设备临床剂量分布的交换探测器测量方法 |
-
2010
- 2010-04-12 CN CN2010101437934A patent/CN101797423B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060273746A1 (en) * | 2005-05-19 | 2006-12-07 | Marc Delaperriere | Method for function monitoring in medical accelerator systems |
CN101290354A (zh) * | 2008-06-12 | 2008-10-22 | 中国测试技术研究院电离辐射研究所 | 用于放射治疗设备临床剂量分布的交换探测器测量方法 |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103175547A (zh) * | 2011-12-21 | 2013-06-26 | 北京普源精电科技有限公司 | 一种数据采集装置的参数拟合方法 |
CN109603016A (zh) * | 2018-12-14 | 2019-04-12 | 江苏海明医疗器械有限公司 | 一种医用加速器剂量率稳定控制系统 |
CN112082627A (zh) * | 2020-08-24 | 2020-12-15 | 浙江威星智能仪表股份有限公司 | 一种基于光信号采样的软件计量校准处理方法 |
CN112232381A (zh) * | 2020-09-25 | 2021-01-15 | 国网上海市电力公司 | 主导负荷参数类噪声辨识的模型参数后处理方法和装置 |
CN112232381B (zh) * | 2020-09-25 | 2024-03-01 | 国网上海市电力公司 | 主导负荷参数类噪声辨识的模型参数后处理方法和装置 |
CN113032732A (zh) * | 2021-03-24 | 2021-06-25 | 中国兵器装备集团自动化研究所有限公司 | 一种基于相对误差分段的剂量率拟合方法及系统 |
CN115779281A (zh) * | 2022-10-28 | 2023-03-14 | 中子高新技术产业发展(重庆)有限公司 | 一种混合辐射场的测量装置和方法 |
CN115779281B (zh) * | 2022-10-28 | 2023-10-31 | 中子高新技术产业发展(重庆)有限公司 | 一种混合辐射场的测量装置和方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101797423B (zh) | 2011-06-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101797423B (zh) | 自校正、多曲线拟合数字式剂量监控方法和系统 | |
CN108986939A (zh) | 用于验证核反应堆功率量程功率系数Gk标定值的方法 | |
CN101934930A (zh) | 无定量喂料管控制物料流量均衡的方法 | |
CN105634077A (zh) | 一种充电方法和电子设备 | |
CN103061828B (zh) | 涡轮控制方法、装置和系统以及麻醉机和呼吸机 | |
CN203414450U (zh) | 气体传感器自动校准系统 | |
Bolt et al. | A multi-centre analysis of radiotherapy beam output measurement | |
CN112034699A (zh) | 一种冗余系统高精度散热控制方法 | |
CN107096135A (zh) | 移动检测系统和方法 | |
CN107529270A (zh) | 高能医用电子加速器用的具有自稳幅技术的微波激励源 | |
CN116756969A (zh) | 一种磨煤机的煤粉出力在线计算方法和系统 | |
CN111987720A (zh) | 多通道限额约束下受端电网受电及供电裕度区间评估方法 | |
CN111721181A (zh) | 螺栓群预紧数据测量和无线采集装置、系统及方法 | |
CN103211607A (zh) | X 光机及其控制方法和装置 | |
CN201765087U (zh) | 压力变送器在线校准装置 | |
CN101728829B (zh) | 电网有功和无功综合安全控制系统 | |
CN102551758B (zh) | 放射线检测器、其制造方法和放射线摄影图像拍摄设备 | |
CN117045898B (zh) | 一种智能推麻仪流量控制系统 | |
CN206573950U (zh) | 一种基于dcs的硫酸钾反应室内微负压自动控制装置 | |
CN103976755A (zh) | 一种时间补偿方法及装置 | |
CN203705635U (zh) | 一种新型tvlf探水雷达同步系统 | |
CN115167327B (zh) | 一种新能源汽车零件生产设备的控制方法及系统 | |
CN113917238B (zh) | 基于电磁辐射信号时域频域特征的煤岩动力灾害预警方法 | |
CN111489797A (zh) | 一种病患信息随访系统 | |
CN114706357A (zh) | 基于5g通信的工控设备信息化控制指令传输系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |