论文部分内容阅读
摘 要：医用诊断X射线机已广泛应用于各级医疗机构。广大计量检定人员在检定X射线机时需要进行大量数据处理工作，既繁杂，又容易出错。编制一份电子原始记录模板可以大大减轻检定人员的计算量，从而提高工作效率。
关键词：电子原始记录；X射线机；数据处理
0 引言
传统的工作模式中，检定员需要将现场检测数据记录在纸质原始记录上，进行数据处理后，将结果录入电脑，制作检定证书。大量枯燥重复的数据处理工作往往伴随着不小的差错率，直接影响工作效率和证书质量。为了解决这一问题，本所从2012年开始实行原始记录电子化，利用电子原始记录系统减少了大量繁复的计算工作，降低了数据的差错率，提高了工作效率和证书质量[1]。该系统中嵌入了Excel电子表格软件，利用其强大的数据处理和函数运算功能，可以方便的编制电子记录模板，自动对测量数据进行运算分析，生成检定证书。下面，本文以医用诊断X射线机的电子模板为例，介绍一下电子原始记录系统在计量检定中的应用。
1 检定空气比释动能率
根据JJG744-2004《医用诊断X射线辐射源检定规程》，辐射输出的空气比释动能率K=MKtpNK，其中，
使用纸质记录时，检定员需要记录三次剂量仪的示值，算出平均值M，然后根据检定时的温度、气压计算KTP，再乘以校正因子NK，才能得出最后的检定结果，计算量相当之大。有了电子记录模板，我们只要填入三次剂量仪的示值以及当时的温度、气压值，模板便会自动算出空气比释动能率的结果，既方便又快捷。检定员只需要输入温度、气压和3次测量值，电子模板便会自动算出空气比释动能率的检定结果。用到的公式如下：
计算KTP： 在单元格中输入 =（（273.15+C1）/293.15）*（101.3/C2）
计算M： 在单元格中输入 =AVERAGE（F1：F3）
计算K： 在单元格中输入 =F4*C3*C4
也可以直接输入：=AVERAGE（F1：F3）*C3*（（273.15+C1）/293.15）*（101.3/C2），直接计算出空气比释动能率。
公式中统计函数AVERAGE可以求出指定数组的平均值。其他类似的函数还有标准差函数STDEV，最大值函数MAX，最小值函数MIN，求和函数SUM，绝对值函数ABS等。将这些函数与四则运算符有机结合起来便可以方便快捷地处理测量数据，得出辐射输出的重复性、线性、X射线管焦点、光野辐射野一致性等项目的检定结果。这些项目模板的制作过程类似，只需要更改相应的计算公式即可，本文不再赘述。
2 检定辐射输出的质
对于可以直接读出HVL的仪器，直接在模板中记录数据即可。对于使用吸收片测量HVL的，使用电子模板可以减轻很大的工作量。在实际检定中，检定员往往只测量无吸收片和加一次吸收片时的剂量率，利用两点线性插值公式计算HVL。由于采用一阶线性插算法，插值余项较大，给测量带来了较大误差，而增加测量点又会大大增加检定员的计算量。利用电子模板提供的TREND函数，可以解决这一难题。
TREND函数用来返回一条线性回归拟合线的值。利用最小二乘法找到适合已知数组Yn和Xn的直线，并返回指定数组Xi在直线上对应的Yi值。其语法为：
TREND（Yn， Xn， Xi， [const]）
其中，Yn、Xn为已知数组，Xi为所求Yi对应的值，const参数可选，一般不用设置。
我们只需输入铝片厚度和相应的剂量率，模板会自动算出HVL。用到的公式如下：
计算K0/2： 在单元格中输入 =C14/2
计算HVL： 在单元格中输入 =TREND（C13：E13，C14：E14，C15）
如果要增加一组测量点，只需将公式改为=TREND（C13：F13，C14：F14，C15）即可，以此类推。
3 检定X射线管的电压
该项目的检定并不复杂，只需要按照规程要求编制公式即可。难点在于在三个测量点中要挑出一个最大偏差作为检定结果。我们需要用到逻辑函数IF。
IF函数可以对指定的条件进行判断，返回不同的结果。其语法为：
IF（logical_test，value_if_true，value_if_false）
即对逻辑表达式logical_test进行判断，如果结果为true，则返回表达式value_if_true的值；如果结果为false，则返回表达式value_if_false的值。表达式也可以是单元格甚至常数。如图3，利用函数可以很方便的在G21～G23单元格算出三个测量点的相对偏差，接下来我们需要找出其中偏差最大的一个。
由于偏差有正有负，我们不能简单的求最大值或者最小值，而要对其进行分析。在单元格中输入=IF（ABS（MAX（G21：G23））>ABS（MIN（G21：G23）），MAX（G21：G23），MIN（G21：G23）），逻辑表达式ABS（MAX（G21：G23））>ABS（MIN（G21：G23））的作用是判断出绝对值最大的偏差是正偏差还是负偏差，也可以写成MAX（G21：G23）+MIN（G21：G23）>0，结果是一样的。模板中X射线管的电流、加载时间等项目与电压类似，不再赘述。
4 检定结果的判定
对每个单项算出测量结果之后，要依据该项目的技术要求进行判定，得出单项结论并填写到证书的结果页中。我们需要用到IF函数。以X射线管的电压为例，规程要求管电压的相对误差不超过±10%为合格，那么我们可以在模板中输入：=IF（ABS（D24）>10%，"不合格"，"合格"）。这样，只要我们输入测量数据，模板会自动算出测量结果并下结论。
但是，并不是所有的项目判定依据都是单一的。比如辐射输出的质的判定依据，规程的表1列出了多种测量条件下的最小HVL。如果要让电子模板进行判定，输入的判定条件必须是单一的。要实现这一目的，我们需要用到查询函数HLOOKUP。
首先，我们需要在模板中，将几种常用的测量条件及技术要求做成表格。然后，输入判定公式=IF（C16 HLOOKUP函数可以在表格或数值数组的首行查找指定的数值，并在表格或数组中指定行的同一列中返回一个数值。其语法为：
HLOOKUP（lookup_value， table_array， row_index_num， [range_lookup]）
HLOOKUP（C12，D10：H11，2，0）的作用是，在技术要求表格D10～H11的第一行，查找测量条件C12的值70kV所在的列，并返回表格中该列指定的第2行的值2.1。当C12中输入70时，D16等效于IF（C16<2.1，"不合格"，"合格"），結果为合格；如果在C12中输入90，D16等效于IF（C16<2.5，"不合格"，"合格"）结果为不合格。
HLOOKUP函数是查找某列的值，如果要查找某一行的值，需要用VLOOKUP函数，语法类似。
5结语
综上所述，电子原始记录系统有着强大的数据处理和运算能力。我们仅需要将原始测量数据输入电子模板，系统就可以自动算出测量结果，加以判断，并生成证书。这大大减少了检定员的工作量，显著提高了工作效率，值得推广。
参考文献：
[1]张俊峰，宫宁生，计量检测原始记录电子化系统设计与实现[J]，中国计量，2013（03）：101-103.
[2]马明涛，徐海鹏，许秋茹，医用诊断X射线机半值层计算方法探析[J]，计量技术，2002（03）：44-45.
关键词：电子原始记录；X射线机；数据处理
0 引言
传统的工作模式中，检定员需要将现场检测数据记录在纸质原始记录上，进行数据处理后，将结果录入电脑，制作检定证书。大量枯燥重复的数据处理工作往往伴随着不小的差错率，直接影响工作效率和证书质量。为了解决这一问题，本所从2012年开始实行原始记录电子化，利用电子原始记录系统减少了大量繁复的计算工作，降低了数据的差错率，提高了工作效率和证书质量[1]。该系统中嵌入了Excel电子表格软件，利用其强大的数据处理和函数运算功能，可以方便的编制电子记录模板，自动对测量数据进行运算分析，生成检定证书。下面，本文以医用诊断X射线机的电子模板为例，介绍一下电子原始记录系统在计量检定中的应用。
1 检定空气比释动能率
根据JJG744-2004《医用诊断X射线辐射源检定规程》，辐射输出的空气比释动能率K=MKtpNK，其中，
使用纸质记录时，检定员需要记录三次剂量仪的示值，算出平均值M，然后根据检定时的温度、气压计算KTP，再乘以校正因子NK，才能得出最后的检定结果，计算量相当之大。有了电子记录模板，我们只要填入三次剂量仪的示值以及当时的温度、气压值，模板便会自动算出空气比释动能率的结果，既方便又快捷。检定员只需要输入温度、气压和3次测量值，电子模板便会自动算出空气比释动能率的检定结果。用到的公式如下：
计算KTP： 在单元格中输入 =（（273.15+C1）/293.15）*（101.3/C2）
计算M： 在单元格中输入 =AVERAGE（F1：F3）
计算K： 在单元格中输入 =F4*C3*C4
也可以直接输入：=AVERAGE（F1：F3）*C3*（（273.15+C1）/293.15）*（101.3/C2），直接计算出空气比释动能率。
公式中统计函数AVERAGE可以求出指定数组的平均值。其他类似的函数还有标准差函数STDEV，最大值函数MAX，最小值函数MIN，求和函数SUM，绝对值函数ABS等。将这些函数与四则运算符有机结合起来便可以方便快捷地处理测量数据，得出辐射输出的重复性、线性、X射线管焦点、光野辐射野一致性等项目的检定结果。这些项目模板的制作过程类似，只需要更改相应的计算公式即可，本文不再赘述。
2 检定辐射输出的质
对于可以直接读出HVL的仪器，直接在模板中记录数据即可。对于使用吸收片测量HVL的，使用电子模板可以减轻很大的工作量。在实际检定中，检定员往往只测量无吸收片和加一次吸收片时的剂量率，利用两点线性插值公式计算HVL。由于采用一阶线性插算法，插值余项较大，给测量带来了较大误差，而增加测量点又会大大增加检定员的计算量。利用电子模板提供的TREND函数，可以解决这一难题。
TREND函数用来返回一条线性回归拟合线的值。利用最小二乘法找到适合已知数组Yn和Xn的直线，并返回指定数组Xi在直线上对应的Yi值。其语法为：
TREND（Yn， Xn， Xi， [const]）
其中，Yn、Xn为已知数组，Xi为所求Yi对应的值，const参数可选，一般不用设置。
我们只需输入铝片厚度和相应的剂量率，模板会自动算出HVL。用到的公式如下：
计算K0/2： 在单元格中输入 =C14/2
计算HVL： 在单元格中输入 =TREND（C13：E13，C14：E14，C15）
如果要增加一组测量点，只需将公式改为=TREND（C13：F13，C14：F14，C15）即可，以此类推。
3 检定X射线管的电压
该项目的检定并不复杂，只需要按照规程要求编制公式即可。难点在于在三个测量点中要挑出一个最大偏差作为检定结果。我们需要用到逻辑函数IF。
IF函数可以对指定的条件进行判断，返回不同的结果。其语法为：
IF（logical_test，value_if_true，value_if_false）
即对逻辑表达式logical_test进行判断，如果结果为true，则返回表达式value_if_true的值；如果结果为false，则返回表达式value_if_false的值。表达式也可以是单元格甚至常数。如图3，利用函数可以很方便的在G21～G23单元格算出三个测量点的相对偏差，接下来我们需要找出其中偏差最大的一个。
由于偏差有正有负，我们不能简单的求最大值或者最小值，而要对其进行分析。在单元格中输入=IF（ABS（MAX（G21：G23））>ABS（MIN（G21：G23）），MAX（G21：G23），MIN（G21：G23）），逻辑表达式ABS（MAX（G21：G23））>ABS（MIN（G21：G23））的作用是判断出绝对值最大的偏差是正偏差还是负偏差，也可以写成MAX（G21：G23）+MIN（G21：G23）>0，结果是一样的。模板中X射线管的电流、加载时间等项目与电压类似，不再赘述。
4 检定结果的判定
对每个单项算出测量结果之后，要依据该项目的技术要求进行判定，得出单项结论并填写到证书的结果页中。我们需要用到IF函数。以X射线管的电压为例，规程要求管电压的相对误差不超过±10%为合格，那么我们可以在模板中输入：=IF（ABS（D24）>10%，"不合格"，"合格"）。这样，只要我们输入测量数据，模板会自动算出测量结果并下结论。
但是，并不是所有的项目判定依据都是单一的。比如辐射输出的质的判定依据，规程的表1列出了多种测量条件下的最小HVL。如果要让电子模板进行判定，输入的判定条件必须是单一的。要实现这一目的，我们需要用到查询函数HLOOKUP。
首先，我们需要在模板中，将几种常用的测量条件及技术要求做成表格。然后，输入判定公式=IF（C16
HLOOKUP（lookup_value， table_array， row_index_num， [range_lookup]）
HLOOKUP（C12，D10：H11，2，0）的作用是，在技术要求表格D10～H11的第一行，查找测量条件C12的值70kV所在的列，并返回表格中该列指定的第2行的值2.1。当C12中输入70时，D16等效于IF（C16<2.1，"不合格"，"合格"），結果为合格；如果在C12中输入90，D16等效于IF（C16<2.5，"不合格"，"合格"）结果为不合格。
HLOOKUP函数是查找某列的值，如果要查找某一行的值，需要用VLOOKUP函数，语法类似。
5结语
综上所述，电子原始记录系统有着强大的数据处理和运算能力。我们仅需要将原始测量数据输入电子模板，系统就可以自动算出测量结果，加以判断，并生成证书。这大大减少了检定员的工作量，显著提高了工作效率，值得推广。
参考文献：
[1]张俊峰，宫宁生，计量检测原始记录电子化系统设计与实现[J]，中国计量，2013（03）：101-103.
[2]马明涛，徐海鹏，许秋茹，医用诊断X射线机半值层计算方法探析[J]，计量技术，2002（03）：44-45.