眼科超声测量设备中应用的自适应时间-增益补偿方法
技术领域
本发明涉及一种眼科A型眼轴长度超声测量技术。特别是涉及一种能够提高眼轴长度测量精度的眼科超声测量设备中应用的自适应时间-增益补偿方法。
背景技术
目前医学超声诊断、测量设备普遍采用的技术是“脉冲-回波”探测法,即超声换能器在电脉冲的激励下辐射超声波进入人体,再接受人体组织的反射波转换成电信号进行显示,形成“声像图”。由于不同人体组织对超声波的反射特性不同,正常组织和病变组织的反射波存在差异,医生即可依据“声像图”反映的组织结构特性为临床诊断提供依据。
“时间-增益补偿(TGC)”是超声诊断设备中普遍应用的技术。由于人体组织对超声波的衰减特性,超声波在人体内的传播随深度的增加而衰减,因此对反射波的放大器增益需要随之提升,对衰减进行补偿。
目前在超声诊断设备中普遍应用的“时间-增益补偿(TGC)”技术,是依据超声波在软组织中的平均衰减系数斜率(0.6~0.7)dB/(cm·MHz),以线性或其他曲线方式对增益进行补偿,并附之手动分段调节。
A型超声是以反射波的脉冲幅度表示界面反射波的强度,其最主要的临床应用领域是眼科用于眼轴长度AL的测量,为人工晶体植入术的实施提供数据。依据大量临床实践的相关文献报道,在人工晶体植入术中,每1mm眼轴长度的测量误差将产生屈光度约300度的偏差,因此提高眼轴长度AL的测量精度具有重要意义。
现有A型超声眼轴长度AL的测量技术是:通过对眼球前房深度ACD、晶体厚度LENS、玻璃体长度VITR的测量,计算出眼轴长度AL,即:
AL=ACD+LENS+VITR=v1t1+v2t2+v3t3
式中v1、v2、v3分别为声波在前房、晶体、玻璃体中传播的速度,均为已知量。
A型超声眼轴长度测量即依据该三种不同组织的界面反射,测量出t1、t2、t3,计算出各段距离及眼轴长度AL。
由于A型超声眼轴长度测量仅需要获得超声脉冲在前房、晶体、玻璃体和视网膜之间3个声学界面的反射波,而该三种组织由于声学特性的较大差异导致界面反射较强,因此在仪器设计中设定一个适当的“阈值”,幅度超过该阈值的脉冲被识别并被作为时间计算的依据,如图1所示。
由于界面反射波上升的前沿具有一定的斜率,到达触发阈值的时间存在延迟,反射波的幅度越低、延迟时间越长;因此界面反射波的幅度差异引起的延迟时间差异是造成测量系统误差的重要原因。
为此,眼科超声测量设备一般也采用现有的“时间-增益补偿”技术,以图减小测量误差。
现有的“时间-增益补偿”技术用于眼科的不足之处及其原因是:
1)眼球组织结构声学特性的特殊性
在眼轴长度测量中需要测量的前房、晶体、玻璃体三种组织的声衰减系数斜率存在很大差异,前房和玻璃体的声衰减系数斜率极低,仅为0.10dB/(cm·MHz),而晶体的声衰减系数斜率达到2.0dB/(cm·MHz),与前者差异达20倍。现有的“时间-增益补偿”技术所采用的补偿曲线对前房和玻璃体会造成过补偿、而对晶体则会有较严重的欠补偿,造成界面反射波幅度的很大差异而加大了测量误差。
2)眼球组织结构的个体差异
据“眼科影像诊断学”等资料报导,对于正常人群,眼球前房深度的典型值为2.27±0.541mm,晶体厚度的典型值为4.07±0.5mm,玻璃体长度的典型值为15.7±1.1mm,具有较大的离散型;近视眼的眼轴长度比正常人长,而远视眼的眼轴长度则偏短;何况患有白内障的患者晶体混浊程度不同所造成的声衰减系数斜率的变化。现有的“时间-增益补偿”技术无论采用那种近似眼球声学特性的补偿曲线,也不可能对每个患者个体都达到最佳的补偿效果。
所以,现有的“时间-增益补偿”技术不能准确地适应眼球结构的生理特性和个体差异。
按目前该类产品普遍采用的探头频率为10MHz计算,其上升时间约为25ns;现有的“时间-增益补偿”技术,由于界面反射波幅度差异而达到触发阈值的时间差可达20ns(如图1所示),按超声波在眼组织的平均声速1600m/s计算,产生的测量系统误差可达0.032mm。
我国现行有效的医药行业标准YY0107-2005《眼科A型超声测量仪》对眼轴长度AL的测量误差要求为不大于±0.1mm,现有的“时间-增益补偿”技术引起的测量系统误差,已成为影响测量精度的重要因素。
发明内容
本发明所要解决的技术问题是,提供一种能够针对受检者眼球生理结构的具体情况,简便的自动调用符合其眼球结构特性的“时间-增益补偿”曲线,大大降低现有技术的测量系统误差,提高眼轴长度测量的精度,为临床诊断和人工晶体植入术的实施提供更准确数据的眼科超声测量的自适应时间-增益补偿方法。
本发明所采用的技术方案是:一种眼科超声测量设备中应用的自适应时间-增益补偿方法,包括如下阶段:
1)分别建立眼轴长度测量中需要测量的前房、晶体、玻璃体三种组织的标准时间-增益补偿曲线,存贮TGC存贮器中;
2)依据临床检查结果,在眼科A型超声测量设备的人机对话界面输入患者相关信息;
3)进行预测量;
4)TGC存贮器依据获得的数据,自动进行自适应TGC曲线拟合,完成符合受检者个体情况的“时间-增益补偿”曲线,用于正式测量;
5)进行正式测量,预测量完成之后,在完成“自适应TGC曲线”拟合的同时,自动开始使用上述适应受检者具体情况的“时间-增益补偿”曲线进行正式测量,获得眼轴长度AL及其前房深度AC、晶体厚度LEN和玻璃体长度VITR各段的数据并进行显示,作为测量的结果输出。
所述的阶段1包括如下步骤:
a)前房、晶体、玻璃体的声衰减系数斜率数据如下:
前房、玻璃体的声衰减系数斜率取0.10dB/(cm·MHz);
晶体的声衰减系数斜率:正常眼取2.0dB/(cm·MHz);对于患有白内障的眼,按常规公知的临床诊断晶体混浊程度4级分类法,分4级对上述声衰减系数斜率进行校正;
求出不同混浊程度白内障眼的声衰减系数斜率的平均值;与正常眼的声衰减系数斜率相比较,求出校正因子α1~α4。
b)对眼组织各部分的超声频率修正:
前房部分由于声衰减系数斜率极低,不需要校正,该超声频率取自超声探头的标称工作频率f0。
晶体部分由于衰减系数斜率较大,需要对超声频率进行修正,校正方法为:
以“脉冲-回波法”测量多个超声换能器的反射波频谱,取其中反射波频谱的中心频率等于换能器标称中心频率f0作为典型频谱;
从典型频谱中读出不同频率对应的回波幅值dB,将连续的频谱简化为一组“频率-幅值”数据;
将上述不同频率对应的回波幅值用声衰减系数斜率2.0dB/(cm·MHz)校正,得到超声波通过单位晶体厚度后的一组新的“幅值-频率”关系数据,从该组数据中可获得超声波传播单位路径衰减后的中心频率值;
重复上述步骤,直至覆盖晶体厚度可能的最大值,得到中心频率的偏移与晶体厚度的关系曲线,从该曲线得出超声波通过一定厚度的晶体后修正的中心频率f1。
玻璃体部分由于声衰减系数斜率极低,不需要校正,该超声频率取自超声波通过晶体后修正的中心频率f1。
c)将上述求出的声衰减随距离变化的时间-增益补偿曲线存入TGC存贮器中。
所述的阶段2包括如下步骤:
a)选择“正常眼”:声衰减系数校正因子α自动选择α=1,即不需校正;
b)如果患者属于白内障眼,按临床诊断的晶体混浊程度选择1、2、3或4,自动调用已存入TGC存贮器中的声衰减系数斜率校正因子α1、α2、α3或α4。
所述的阶段3包括如下步骤:
a)完成患者信息输入后,眼科A型超声测量设备自动按常规方法、使用现有时间-增益补偿技术对眼轴长度进行一组含有多次发射的预测量,获得受检者前房深度AC、晶体厚度LEN和玻璃体长度VITR的一组数据,通过对该组数据的统计学处理,求出该组数据的平均值和标准差;
b)将该组数据的平均值和标准差与依据临床经验设定的判据进行比较,如果该组数据的重复性不符合判据要求,给出提示并自动进行下一次预测量;如果预测量数据的重复性符合判据要求,将该组数据的平均值AC、LEN、VITR输出至TGC存贮器,作为TGC存贮器中声衰减随距离变化函数曲线的距离值。
本发明的眼科超声测量的自适应时间-增益补偿方法,具有如下特点:
1、对A型超声眼轴长度测量中应用的“时间-增益补偿”曲线的变量“频率”,进行通过晶体后的偏移校正。
2、对A型超声眼轴长度测量中应用的“时间-增益补偿”曲线的变量“距离”,自动调用对受检者预测量获得的实际数据。
本发明与现有的、由生产厂商预置的“时间-增益补偿”曲线相比,可以方便的依据受检者个体的实际情况,自动采用适应受检者自身特点的“时间-增益补偿”曲线,有效降低了界面反射波的幅度差异引起的延迟时间差异,提高了眼轴长度AL的测量精度,为临床诊断和人工晶体植入术的实施提供更准确的数据。
附图说明
图1是A型超声眼轴长度测量工作原理示意图;
图2是本发明方法的流程图。
其中:
1:眼科A超探头 2:前房
3:晶体 4:玻璃体
5:各声学界面的反射波 6:触发阈值
7:触发脉冲——其脉冲间隔时间为计算眼轴长度的基准
具体实施方式
下面结合实施例和附图对本发明的眼科超声测量设备中应用的自适应时间-增益补偿方法做出详细说明。
本发明需要解决的技术问题是,针对受检者眼球生理结构的具体情况,简便的自动调用符合其眼球结构特性的“时间-增益补偿”曲线。
“眼科超声测量设备”在临床应用中的直接目的不是获得诊断结果或健康状况,而只是从活的人体或动物体获取作为中间结果的“眼轴长度”这一生理参数;该信息本身不能够直接得出疾病的诊断结果或健康状况,属于为实施人工晶体植入术而需要测量的中间数据。
本发明的眼科超声测量设备中应用的自适应时间-增益补偿方法的特点是:
1、对A型超声眼轴长度测量中应用的“时间-增益补偿”曲线的变量“频率”,进行通过晶体后的偏移校正。
2、对A型超声眼轴长度测量中应用的“时间-增益补偿”曲线的变量“距离”,自动调用对受检者预测量获得的实际数据。
如图2所示,本发明的眼科超声测量设备中应用的自适应时间-增益补偿方法,具体包括如下阶段:
1)分别建立眼轴长度测量中需要测量的前房、晶体、玻璃体三种组织的标准时间-增益补偿曲线,存贮TGC存贮器中。包括如下步骤:
a)前房、晶体、玻璃体的声衰减系数斜率数据取自权威文献报道的、公知公认数据,如下:
前房、玻璃体的声衰减系数斜率取0.10dB/(cm·MHz);
晶体的声衰减系数斜率:正常眼取2.0dB/(cm·MHz);对于患有白内障的眼,按常规公知的临床诊断晶体混浊程度4级分类法,分4级对上述声衰减系数斜率进行校正;
为获取不同混浊程度白内障眼的声衰减系数斜率的校正因子,可通过大样本临床测量超声波在不同混浊程度的晶体的前、后沿反射波幅度的差异,经过公知的统计学处理方法,求出不同混浊程度白内障眼的声衰减系数斜率的平均值;与正常眼的声衰减系数斜率相比较,求出校正因子α1~α4。
b)由于声衰减与超声频率有关,超声波在组织传播过程中,超声频谱的高频成分由于衰减量较大而使频谱发生频移,因此需要对眼组织各部分的超声频率修正:
前房部分由于声衰减系数斜率极低,不需要校正,该超声频率取自超声探头的标称工作频率f0(中心频率):
晶体部分由于衰减系数斜率较大,需要对超声频率进行修正,校正方法为:
以“脉冲-回波法”测量多个超声换能器的反射波频谱,取其中反射波频谱的中心频率等于换能器标称中心频率f0的作为典型频谱;
从典型频谱中读出不同频率对应的回波幅值dB,将连续的频谱简化为一组“频率-幅值”数据;
将上述不同频率对应的回波幅值用声衰减系数斜率2.0dB/(cm·MHz)校正,得到超声波通过单位晶体厚度后的一组新的幅值-频率关系数据,从该组数据中可获得超声波传播单位路径衰减后的中心频率值;
重复上述步骤,直至覆盖晶体厚度可能的最大值,可得到中心频率的偏移与晶体厚度的关系曲线,从该曲线可以得出超声波通过一定厚度的晶体后修正的中心频率f1。
玻璃体部分由于声衰减系数斜率极低,不需要校正,该超声频率取自超声波通过晶体后修正的中心频率f1。
c)将上述求出的声衰减随距离变化的时间增益补偿曲线存入TGC存贮器中。
2)依据临床检查结果,在眼科A型超声测量设备的人机对话界面输入患者相关信息。包括如下步骤:
a)选择“正常眼”:声衰减系数校正因子α自动选择α=1,即不需校正;
b)如果患者属于白内障眼,按临床诊断的晶体混浊程度选择1、2、3或4,自动调用已存入TGC存贮器中的声衰减系数斜率α1、α2、α3或α4。
3)进行预测量。包括如下步骤:
a)完成患者信息输入后,眼科A型超声测量设备自动按常规方法、使用现有时间-增益补偿技术对眼轴长度进行含有多次发射的一组预测量,约在100ms时间内即可获得受检者前房深度AC、晶体厚度LEN和玻璃体长度VITR的一组数据,通过对该组数据的统计学处理,求出该组数据的平均值和标准差;
b)将该组数据的平均值和标准差与依据临床经验设定的判据进行比较,如果该组数据的重复性不符合判据要求,给出提示并自动进行下一次预测量;如果预测量数据的重复性符合判据要求,将该组数据的平均值AC、LEN、VITR输出至TGC存贮器,作为TGC存贮器中声衰减随距离变化函数曲线的距离值。
4)TGC存贮器依据获得的数据,自动进行自适应TGC曲线拟合,完成符合受检者个体情况的“时间-增益补偿”曲线,用于正式测量;
5)进行正式测量。预测量完成之后,在完成“自适应TGC曲线”拟合的同时,自动开始使用上述适应受检者具体情况的“时间-增益补偿”曲线进行正式测量,获得眼轴长度AL及其前房深度AC、晶体厚度LEN和玻璃体长度VITR各段的数据并进行显示,作为测量的结果输出。
本发明的眼科超声测量设备中应用的自适应时间-增益补偿方法具体应用如下:
1、预测量的符合性判据
与当前产品应用的现有技术相同。即将测量值与预设的AC、LEN、VITR合理范围比较,多次重复测量数据的统计学处理等;以排除超声波束未对准眼轴、探头接触角膜压力过大等医生操作手法不当产生的错误数据。
2、声衰减系数校正因子α
正常眼晶体的声衰减系数斜率为2.0dB/(cm·MHz),对于白内障眼,由于晶体混浊的散射会使该数值变化。眼科临床一般将晶体混浊程度分为4级。
衰减系数校正因子确定的最佳方案是:
临床医生可以通过临床B型图像扫描及其他光学检查手段等现有技术,对晶体混浊程度进行级别分类。
依据临床晶体混浊程度级别分类的结果,分别对不同级别的患者采用眼轴方向的A
超扫描,记录晶体前、后沿反射波的幅度差异数据,生产商依据仪器内置的TGC曲线模式,获得其衰减数据;再通过大样本的统计数据处理,求出不同混浊程度白内障眼的声衰减系数斜率的平均值;与正常眼的声衰减系数斜率相比较,求出不同晶体混浊级别的校正因子α1~α4。
3、中心频率校正
由于晶体的声衰减系数斜率远高于其他眼组织,超声波通过晶体后的频谱偏移不能被忽略。
中心频率校正的最佳实施方案是:
1)以“脉冲-回波法”测量多个超声换能器的反射波频谱,取其中反射波频谱的中心频率等于换能器标称中心频率(f0)的作为典型频谱。
2)以中心频率的10%为起点,以1MHz为增量,从典型频谱中读出不同频率对应的回波幅值(dB),直至频率达到中心频率的2倍为止;将连续的频谱简化为一组“频率-幅值”数据。
3)由于反射波通过晶体的路程是晶体厚度的2倍,将步骤2)得到的不同频率对应的回波幅值以0.4mm为单位,用声衰减系数斜率2.0dB/(cm·MHz)校正,即可得到超声波通过晶体厚度为0.2mm后的一组“幅值-频率”关系数据,从该组数据中即可获得超声波传播单位路径衰减后的中心频率值;
4)以0.4mm为增量,重复步骤3),直至覆盖2倍晶体厚度可能的最大值(12mm),即可得到晶体厚度每增加0.2mm对应的中心频率值;
5)获得晶体厚度的预测量数据后,采用线性内插法即可获得超声波通过晶体后修正的中心频率f1。