实时调整注释点的制作方法

1.本发明整体涉及体内医学规程和器械，并且具体地涉及体内心脏心电图(ecg)感测和可视化。


背景技术：

2.当测量和注释由大量电极生成的内部心电图(iecg)信号时，可能期望处理这些信号(例如，通过计算机)，以便减少嵌入的噪声。
3.存在用于此类iecg信号处理的各种方法。例如，美国专利申请公布2016/0089048描述了一种确定四个或更多个多通道心脏电描记图信号的局部激动时间(lat)的自动方法，该多通道心脏电描记图信号包括心室通道、映射通道和多个参考通道。
4.另一个示例是美国专利申请公布2017/0311833，该专利申请公布描述了一种用于通过对身体内的空间电生理(ep)模式进行测量、分类、分析和标测来诊断心律失常并引导导管治疗的系统。
5.又一个示例是美国专利申请公布2017/0042436，该专利申请公布描述了将对多次心跳进行的测量自动整合到单个心脏标测图中的系统和方法。


技术实现要素：

6.本文描述的本发明的实施方案提供了一种包括信号采集电路和处理单元的系统。信号采集电路被配置为接收由患者的心脏中的心脏内探头的多个电极采集的多个心脏内信号。处理单元被配置为选择一组心脏内信号，根据似然准则从该组中的心脏内信号中的每个心脏内信号中提取相应的最可能注释值，在该组中识别其最可能注释值在该组中统计上偏差超过预定义的偏差的量度的心脏内信号，根据似然准则从具有统计上偏差的注释值的心脏内信号中提取至少第二最可能注释值，并且响应于该第二最可能注释值的统计偏差，为对应的心脏内信号选择有效注释值。
7.在一些实施方案中，处理单元被配置为根据注释值的标准分数来定义偏差的量度。在所公开的实施方案中，处理单元被配置为计算心脏内信号上的注释值的偏差，该心脏内信号由在心脏中彼此相距不超过预定义的距离的空间相关电极的选择的子集采集。
8.在示例性实施方案中，处理单元被进一步配置为：对于至少统计上偏差的最可能注释值，识别具有减小的似然秩的一组另选注释值；以及响应于另选注释值的统计偏差和似然秩来选择有效注释值。在另一个实施方案中，注释值包括局部激动时间(lat)。
9.在一些实施方案中，处理单元被配置为通过在心动周期中找到给定心脏内信号的极值来提取给定心脏内信号中的最可能注释值，并且通过找到心脏内信号的第二最高局部极值来提取第二最可能注释值。在其他实施方案中，处理单元被配置为通过在心动周期中找到给定心脏内信号的极值导数来提取给定心脏内信号中的最可能注释值，并且通过找到该导数的第二最高局部极值来提取第二最可能注释值。
10.根据本发明的一个实施方案，另外提供了一种方法，包括接收由患者的心脏中的
心脏内探头的多个电极采集的多个心脏内信号。选择一组心脏内信号。根据似然准则，从该组中的心脏内信号中的每个心脏内信号中提取相应的最可能注释值。在该组中识别心脏内信号，该心脏内信号的最可能注释值在该组中统计上偏差超过预定义的偏差的量度。根据似然准则，从具有统计上偏差的最可能注释值的心脏内信号中提取至少第二最可能注释值。响应于第二最可能注释值的统计偏差，为对应的心脏内信号选择有效注释值。
11.根据本发明的实施方案，还提供了一种获得心脏内心电图信号的有效局部激动时间的方法。该方法包括采集表示心脏内心电图(ecg)信号的一组数字化信号，以及从该组ecg信号中提取有效局部激动时间的第一最佳估计。计算该组ecg信号的统计特性。基于组统计特性计算该组ecg信号中的每个信号的标准分数。将每个信号的标准分数与预设限值进行比较。有效局部激动时间的第一最佳估计被具有在预设限值内的标准分数的信号的局部激动时间替换。
12.结合附图，通过以下对本发明的实施方案的详细描述，将更全面地理解本发明，其中：
附图说明
13.图1为根据本发明的实施方案的用于心脏内ecg信号的多通道测量的电解剖系统的示意性图解；
14.图2为根据本发明的实施方案的示意性地示出一组空间相关的电极在单个心动周期中采集信号的图；
15.图3为根据本发明的实施方案的示意性地示出用于增强注释值的可靠性的方法的流程图；并且
16.图4为根据本发明的实施方案的示意性地示出用于增强一组空间相关的lat值的注释值的可靠性的方法的流程图。
具体实施方式
17.概述
18.基于心脏内探头的(例如，基于导管的)心脏诊断和治疗系统可在侵入式规程期间测量多个心脏内信号，诸如心电图(ecg)。此类系统可使用装配在探头的远侧端部处的电极(在下文中也称为“远侧电极”)来采集多个心脏内信号。所测量的信号通常被分析，并且局部激动信号(lat)被注释，它们可用于向医师提供视觉心脏信息，诸如患者的心脏内病理电图源的3
‑
d标测，并支持矫正医疗规程，诸如消融。
19.所测量的信号通常较弱，信噪比(snr)较低。此外，病理心电图诸如由心房扑动或心房纤颤引起的那些病理心电图可在心动周期中表现出多个峰，这使lat值的测量复杂化。另一方面，使用许多电极，因此系统从电极接收的数据可能存在一些冗余。
20.本文所公开的本发明的实施方案提供了基于心脏内探头的电解剖测量和分析系统以及使用远侧电极收集的信号的统计特性的方法，以改善所收集的数据的质量和可靠性。这些方法是快速的，并且因此可实时完成；例如，在侵入式规程期间。
21.尽管下文的描述是指局部激动时间(lat)的注释值，但在本发明的各种实施方案中，可以必要的变更使用其他合适的信号参数。因此，术语“注释值”包括其他参数以及与
lat相关的那些参数。
22.根据实施方案，可通过使用预设准则分析心脏内信号来估计所测量的心脏内信号的最可能lat注释值。可使用若干lat估计准则，包括(但不限于)最大信号电压事件、最小信号电压事件、正电压斜率事件的最大变化率和负电压斜率事件的最大变化率。然后根据选择的准则将lat估计设定为事件的发生时间。事件处的电压或斜率的值将在本文中称为“y值”。
23.如上所述，由于噪声和/或心脏病理，此类lat估计有时可能错误地识别局部极值(根据准则，最大值或最小值)而不是期望的极值。在那些情况下，如果处理器使用同一准则重新分析信号，但搜索第二极值、第三极值、第四极值或其他极值，则可实现更准确的lat估计。(为了避免最小/最大/高/低之间的混淆，我们将在下文使用术语诸如“最佳lat估计”、“第二最佳”、“第三最佳”等。对于最大电压或最大dv/dt事件，最佳估计是最高电压或最高dv/dt；第二最佳是第二最高局部最大值，等等。类似地，对于最小电压或最小dv/dt事件，最佳估计是最小电压或最小dv/dt；第二最佳估计是第二最低局部最小值；等等。)
24.上述lat估计有时被称为“lat注释值”。
25.本文所公开的技术假设，没有噪声和不规则的电流连接，物理上彼此靠近的电极(“相邻电极”)和/或在时间上相邻的心跳处采集的信号表现出类似的注释值。从相邻电极提取的信号和从此类信号注释的lat值被称为“空间相关”，而从时间相邻的心跳提取的信号和从此类信号注释的lat值被称为“时间相关”。空间和/或时间相关的信号和注释值将统称为“相关”。
26.根据本发明的实施方案利用相关信号的预期类似性来提高可视化lat值的可靠性。
27.在一些实施方案中，处理器例如使用上述四个准则中的一个准则来注释一组相关的心脏内电生理信号的lat值。然后处理器计算该组的统计特性，并且使用响应于统计特性的偏差的量度来确定信号中的哪些信号基本上偏离该组的lat值(例如，lat值不在距组平均值的预设距离内的信号)。
28.根据实施方案，对于基本上偏离组值的信号，处理器决定使用第二最佳估计、第三最佳估计等重新评估lat值，直到处理器找到基本上不偏离组值的lat值。在实施方案中，统计特性包括该组的lat值的平均值(例如，)和标准偏差(例如，)。偏差的量度是lat值的标准分数(被定义为该值与平均值之间的差值除以标准偏差)，将其与预设限值进行比较。例如，基本上偏离的lat值可被定义为比平均值大超过3.5个标准偏差(标准分数＝3.5)，或比平均值低超过1.5个标准偏差(标准分数＝1.5)的任何值。
29.概括地说，根据本发明的实施方案，通过计算注释值的统计特性，并且使用接下来的最佳估计重新计算基本上偏离组值的lat值，可实时改善为用户可视化的空间相关的心脏内信号的一组注释值的质量和可靠性。
30.系统描述
31.图1为根据本发明的实施方案的用于心脏内ecg信号的多通道测量的电解剖系统21的示意性图解。在一些实施方案中，系统21用于心脏的电解剖标测。
32.图1描绘了医师22使用电解剖导管23来执行患者25的心脏24的电解剖标测。导管23在其远侧端部处包括可以是机械柔性的一个或多个臂26，一个或多个远侧电极27耦接到
所述一个或多个臂中的每一个臂。应当理解，尽管图1描绘了具有五个臂的导管，但在根据本发明的另选实施方案中可使用其他类型的导管。电极通过接口32耦接到处理器34。
33.在电解剖标测规程期间，跟踪系统用于跟踪远侧电极27的心脏内位置，使得所采集的电生理信号中的每一个电生理信号可与已知的心脏内位置相关联。跟踪系统的一个示例是在美国专利8456182中描述的有效电流位置(acl)系统。在acl系统中，处理器基于在远侧电极27中的每一个远侧电极与耦接到患者25的皮肤的多个表面电极28之间测量的阻抗来估计远侧电极的相应位置(为了便于说明，图1中仅示出一个表面电极)。然后，处理器可将从远侧电极27接收到的任何电生理信号与采集信号的位置相关联。
34.在一些实施方案中，多个远侧电极27从心脏24的心腔的组织采集心脏内ecg信号。处理器包括被耦接以从接口32接收心脏内信号的信号采集电路36、存储数据和/或指令的存储器38以及处理单元42(例如，cpu或其他处理器)。
35.信号采集电路36使心脏内信号数字化，以便产生多个数字信号。采集电路然后将数字化信号传送到包括在处理器28中的处理单元42。
36.在其他任务中，处理单元42被配置为根据来自一组lat估计准则的选择的准则从信号中提取注释参数，诸如局部激动时间(“lat”)，包括(但不限于)最大ecg电压、最小ecg电压、ecg电压的最大正变化率和ecg电压的最大负变化率。当处理单元根据例如最大ecg电压准则估计lat时，lat值将等于ecg信号的y值处于最大值(对于每个心动周期)时的时间。类似地，估计的lat可等于当y值处于最小值时、当y值的一阶导数处于其最大值时，或者当一阶导数处于其最小值时的时间。
37.处理单元42还被配置为计算统计特性，诸如可能类似的相邻信号的组的注释参数的平均值(在当前上下文中，相邻信号是指来自彼此靠近定位的电极的信号(“空间相关”))。
38.根据实施方案，由于差的电流连接、来自各种源的感应噪声、信号采集电路36中的噪声或来自任何其他源的噪声，所采集的ecg信号可能是嘈杂的。此外，病理ecg信号(诸如与心房扑动和右束分支阻塞相关联的那些病理ecg信号)可在每个心动周期表现出多个电压峰值和/或下降，以及多个dv/dt峰值和下降。通过比较该组的信号中的每个信号的lat值，处理单元可确定显著偏离该组的值的lat值可能是错误的。不太可能是错误的lat值(例如，最佳估计lat值)在下文被称为有效lat值。
39.在实施方案中，对于lat估计准则中的任一个lat估计准则，除了最可能值(其对应于绝对最大值或最小值)之外，处理单元还可提取具有下降似然性的一系列另选lat值。例如，对于最大电压准则，处理单元可被配置为找到ecg信号的局部最大值，然后从第二最高最大值注释第二最佳lat估计，从第三最高最大值注释第三最佳估计，等等。以类似的方式，处理单元可通过找到电压的局部最大值或最小值以及ecg信号的电压的一阶导数来构造用于四个准则中的任一个准则的一系列下降似然性的另选lat值。
40.如果处理单元42确定lat值无效，则处理单元可尝试另选lat值，从第二最佳估计开始，并且前进至不太可能的lat值，直到找到有效值(例如，不太可能错误的lat值)(如果处理单元未找到任何有效值，则处理单元可例如丢弃信号，或使用可能错误的最佳估计)。
41.在一些实施方案中，处理单元42例如使用有效注释值以构造心脏的电解剖标测图50上并在屏幕52上向医师22实时显示该标测图。另选地，处理单元42可以任何其他合适的
方式呈现有效注释值。
42.图1所示的示例性图示完全是为了概念清晰而选择的。在本发明的另选实施方案中，例如，也可以通过在多对表面电极28之间施加电压梯度并利用远侧电极27测量所得电势(例如，使用由加利福尼亚州欧文市(irvine，california)的biosense
‑
webster生产的4技术)来进行位置测量。因此，本发明的实施方案适用于任何位置感测方法。
43.可等效地采用其他类型的导管，诸如导管(由biosense
‑
webster公司生产)或篮状导管。接触传感器可装配在电解剖导管23的远侧端部处。可以类似的方式在远侧电极27上利用其他类型的电极(诸如用于消融的电极)来采集心脏内电生理信号。
44.图1主要示出了与本发明的实施方案相关的部分。省略了其他系统元件，诸如外部ecg记录电极及其连接部。省略了各种ecg记录系统元件，以及用于对电路进行滤波、数字化、保护等的元件。
45.在可选的实施方案中，读出专用集成电路(asic)用于测量心脏内ecg信号。用于路由信号采集电路36的各种元件可以硬件实现，例如，使用一个或多个分立部件，诸如现场可编程门阵列(fpga)或asic。在一些实施方案中，信号采集电路36和/或处理单元42中的一些元件可以软件实现，或者通过使用软件元件和硬件元件的组合来实现。
46.处理单元42通常包括通用处理器，该通用处理器具有经编程以执行本文所述功能的软件。例如，软件可以电子形式通过网络下载，或者另选地或除此之外，该软件可被提供和/或存储在非临时性有形介质诸如磁存储器、光学存储器或电子存储器上。
47.图2为根据本发明的实施方案的示意性地示出一组空间相关的电极在单个心动周期中采集信号的图200。水平轴线表示时间，并且竖直轴线表示电压电平。示出了四个信号——第一电极信号202、第二电极信号204、第三电极信号206和第四电极信号208。每个信号表现出r峰值，然后是t峰值。
48.处理单元42(图1)采用最大电压准则来注释信号的lat值——对应于信号208的r峰值的值210；对应于信号202的r峰值的值212；以及对应于信号206的r峰值的值214。由于噪声和/或病理，对于信号204，t峰值的电压高于r峰值的电压，并且因此信号204的初始lat注释(最佳lat估计)是对应于t信号的峰值的lat值216。然而，信号204lat的第二最佳估计(即，第二最高局部最大值)是由lat值218表示的r峰值。
49.根据图2所示的示例性实施方案，lat值216基本上偏离该组的lat值。因此，处理单元42将尝试替代值，从第二最佳估计开始。由于lat值218不偏离该组的lat值并因此不偏离最可能注释值，因此lat值218是有效的lat值。
50.因此，使用可实时完成的相对快速的方法，处理单元可通过在最佳lat估计显著偏离该组的lat值的情况下检查与低于绝对最大值的局部最大值相关联的lat值来增强注释的lat值的可靠性。
51.图2所示的示例性例证完全是为了概念清晰而选择的。例如，在本发明的另选实施方案中，可使用更空间相关的信号；局部最大值的数量可以多于两个，并且可以使用其他lat估计准则。
52.图3为根据本发明的实施方案的示意性地示出用于增强注释值的可靠性的方法的流程图300。该方法由处理单元42执行(图1)。
53.该方法在“采集信号”步骤302处开始，其中处理单元采集多个心脏内ecg信号。然
后，在“选择组”步骤304中，处理单元从所采集的信号中选择一组相关信号。此种选择可例如使用4技术或通过其他合适的技术来完成。
54.然后，在“选择准则”步骤304中，处理单元选择lat估计准则。选择的准则可为例如最大电压、最小电压、最大dv/dt和最小(最负)dv/dt中的一者。该选择可由用户指示或通过任何其他方式进行选择。
55.使用在步骤304中选择的准则，处理单元在“提取组注释值”步骤306处，通过找到对应于选择的准则的最佳lat估计来提取所有组信号的lat值。例如，在图2中，处理器可查看具有lat值210、lat值212、lat值214、lat值216和lat值218的该组信号(210至218)，并且针对选择的似然准则考虑这些值以选择lat值218(例如，从图2)为步骤306(在图3中)中的“最佳lat估计”(或不太可能是错误的lat值)成为可以图形方式显示的“有效注释值”。
56.应当注意，根据循环310至314、循环316至320或循环324至328中使用的统计测试的结果，处理器可将来自步骤306的最佳lat估计替换为相应统计测试循环316至320的最佳估计(即，将步骤306的第一最佳估计替换为步骤322处的第二最佳估计)，或替换为循环324至328的最佳估计(以将步骤306的第一最佳估计替换为步骤330处的第三最佳估计)。
57.接下来，在“计算统计特性”步骤308中，处理单元计算该组lat值的一个或多个参数，诸如平均值、标准偏差等。
58.在步骤308中计算组统计值之后，处理单元开始循环(步骤310至步骤314)，其中对照该组的统计特性(来自步骤308)检查每个信号的lat注释，并且在lat值显著偏离组lat值的情况下重新计算。即，处理器检查相对于在步骤308中获得的组统计特性的“统计上偏差”超过预定义的统计上偏差的量度的lat值。循环(步骤310至步骤314)从“选择下一个信号”步骤310开始，其中处理单元从组信号(例如，图2中的信号210至信号218)的列表中选择下一个信号，并且在“计算标准分数”步骤312中，计算信号相对于组统计特性的标准分数(即，处理器计算信号的lat值与组平均值之间的标准偏差数)。
59.接下来，在“与限制值进行比较”步骤314中，处理单元检查lat信号的标准分数是否在预设限值内。如果lat值在预设限值(诸如，例如，0.25)内，则处理单元将进入“检查考虑的所有信号”步骤315。如果在步骤315中要考虑更多信号，则处理单元将重新进入步骤310，以检查该组中的下一个信号；然而，如果不再考虑信号，则流程图将结束(或另选地，可针对另一组信号重新运行)。
60.步骤310至步骤314允许处理器识别该组心脏内信号中具有“统计上偏差”注释值的任何心脏内信号；例如——偏离组平均值超过统计标准偏差的预定义量度的注释值。
61.如果在步骤314处，信号的lat值不在预设限值内，则处理单元进入“提取第二最佳lat”步骤316并计算第二最佳lat值(如上文关于步骤312所述)，在“计算标准分数”步骤318中计算第二最佳lat的标准分数，并且在“与限制值进行比较”步骤320中，将标准分数与预设限值进行比较。
62.如果在步骤320中，标准分数在预设限值内，则第二最佳估计优于第一估计；处理单元将进入“通过第二最佳替换最佳”步骤322，并且用第二最佳替换该组信号中的最佳估计lat(来自步骤306)，然后重新进入步骤310以检查下一个信号。然而，如果在步骤320中，标准分数不在预设限值内，则处理单元假设第二最佳估计是错误的并且尝试第三最佳估计——在“提取第三最佳lat”步骤324中，处理单元提取第三最佳值；在“计算标准分数”步
骤326中，处理单元计算第三最佳lat估计的标准分数，并且检查标准分数是否在“与限制值进行比较”步骤328中的预设限值内。
63.如果在步骤328中，标准分数在预设限值内，则第三最佳估计优于步骤306中的第一最佳估计(以及步骤316中的第二最佳估计)；处理单元将进入“通过第三最佳替换最佳”步骤330，用第三最佳替换该组信号中的最佳估计lat(来自步骤306)，然后重新进入步骤310以检查下一个信号。然而，如果在步骤328中，第三最佳估计的标准分数不在限制值内，则不应改变最佳估计lat，并且处理单元重新进入步骤310以检查下一个信号。
64.如上所述，当不再有要考虑的当前组的信号时，流程图在步骤315处结束。在该阶段，处理单元将为组中的所有信号具有注释的lat值。处理单元然后可重新输入用于另一组相关信号的流程图，或执行超出本公开的范围的其他任务。
65.图3中所示的示例性流程图完全是为了概念清晰而选择的。在另选的实施方案中，例如，流程图可包括检查第四和/或更不太可能的最佳估计；在其他实施方案中，仅检查第二最佳估计。在一些实施方案中，计算并存储不太可能的估计作为用于整个组的步骤306的一部分，并且因此步骤312、步骤318和步骤326是冗余的。在一些实施方案中，一些或所有步骤可同时执行。将参考图4描述一些典型的另选的流程图。
66.图4为根据本发明的实施方案的示意性地示出用于增强一组空间相关的lat值的注释值的可靠性的方法的流程图400。该方法由处理单元42(图1)在针对该组的每个成员计算该组的统计特性之后(例如，在图3的步骤308之后；然而，应当注意，图4替换了图3的从步骤310至步骤330的另外步骤)。
67.与图3的示例性实施方案不同，根据图4所示的方法，处理单元在检查该组信号的同时创建一组新的有效lat信号，然后可将其可视化给用户。
68.流程图从“获取最佳lat”步骤402开始，其中处理单元存储来自关于图3的流程图所述的例程的最佳lat和对应的y值。我们可以假设，当处理器例如在图3的步骤306中计算组注释值时，提取了最佳lat和对应的y值。
69.接下来，在“计算标准分数”步骤404中，处理单元计算信号的lat值与组平均值之间的标准偏差数，并且前进至“与限制值进行比较”步骤406，以检查lat信号的标准分数是否在预设限值内。
70.如果lat值的标准分数不在预设限值内，则处理单元将进入“检查是否允许下一次迭代”步骤408。当处理单元第一次进入步骤408时，该处理单元将开始第一迭代，并且每当处理单元重新进入步骤408时，迭代次数将增加。在实施方案中，迭代次数可存在限制
‑
例如，处理器在第三最佳估计之后不太可能跳过估计(并且因此，在步骤408中，在三次迭代之后将不允许下一次迭代)。在另一个示例中，如果不再存在局部最大值或最小值(根据选择的准则)，则处理单元可跳过进一步的迭代。
71.如果在步骤408中允许下一次迭代，则处理单元将进入“计算下一个最佳lat”步骤410，其中处理单元计算下一个最佳lat和对应的y值。
72.根据图4所示的示例性实施方案，为了符合有效lat值的资格，局部最大值(或最小值，根据选择的准则)应当相对于最佳估计是实质性的。例如，如果准则为最大dv/dt；绝对最大值为25mv/sec并且第二局部最大值(第二最佳估计)为1mv/sec，但处理单元可拒绝第二最佳估计，尽管最佳估计可能基本上偏离组平均值。
73.处理单元在步骤410之后执行的“检查y值比率”步骤412中执行该测试。处理单元将第二最佳估计(在上述示例中为1mv/sec)除以最佳估计的y值(25mv/sec)，并且将结果(0.04)与预设阈值进行比较。如果该比率大于阈值，则处理单元前进至“计算标准分数”步骤414并计算相对于组统计值的当前lat估计的标准分数，并且进入“与限制值进行比较”步骤416，以检查下一个最佳lat信号的标准分数是否在预设限值内。如果标准分数不在限制值内，则处理单元将重新进入步骤408，并尝试下一个最佳估计。
74.根据图4所示的示例性实施方案，处理单元可处于丢弃模式，其中无效lat值不被可视化(而不是可视化最佳估计，尽管最可能是错误的)。如果在步骤408中，不再允许迭代，或者如果在步骤412中，比率太低，认为当前最佳lat估计无效，则不存在信号的有效lat值。处理单元进入“检查丢弃模式”步骤418；如果丢弃模式打开，则该过程结束(并且没有有效值将进入该组有效lat值，并且因此当前信号的lat将不被可视化)。如果在步骤418中，丢弃模式关闭，则处理单元将进入“恢复到最佳估计”步骤420，并且将当前lat估计改变为在步骤402中存储的最佳lat。
75.如果在步骤406中，最佳lat在限制值内；或者，如果在步骤416中，较低的最佳估计在限制值内；并且，在步骤420之后，处理单元进入“将lat添加到有效lat”步骤422，并且将当前lat估计添加到该组有效lat估计，并且流程图结束(或者，通常针对该组的下一个信号重新运行)。
76.应当理解，图4中所示的示例性流程图纯粹是为了概念清晰而选择的。在本发明的另选实施方案中，例如，不存在丢弃模式——在实施方案中，无效值从未被可视化，以及在另一个实施方案中，无效值总是被可视化。在又一个实施方案中，处理单元可根据约束较小的测试来决定是否可视化或隐藏非有效lat值；例如——如果值的标准分数偏离组平均值的阈值大于步骤406和步骤416中使用的阈值，则该值将被可视化。在一些实施方案中，无效值被可视化但被清楚地标记。
77.应当理解，上述实施方案以举例的方式被引用，并且本发明不限于上文具体示出和描述的内容。相反，本发明的范围包括上文描述的各种特征的组合和子组合以及它们的变型和修改，本领域的技术人员在阅读上述描述时将会想到该变型和修改，并且该变型和修改并未在现有技术中公开。