CN109411041B - 心电信息处理方法和心电工作站系统 - Google Patents

心电信息处理方法和心电工作站系统 Download PDF

Info

Publication number
CN109411041B
CN109411041B CN201810157158.8A CN201810157158A CN109411041B CN 109411041 B CN109411041 B CN 109411041B CN 201810157158 A CN201810157158 A CN 201810157158A CN 109411041 B CN109411041 B CN 109411041B
Authority
CN
China
Prior art keywords
data
report
user
electrocardiogram
information
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.)
Active
Application number
CN201810157158.8A
Other languages
English (en)
Other versions
CN109411041A (zh
Inventor
曹君
胡传言
刘涛
姜艳
田亮
曹会丽
刘畅
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yocaly Information Science & Technology Co ltd
Original Assignee
Shanghai Lepu Yunzhi Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Lepu Yunzhi Technology Co ltd filed Critical Shanghai Lepu Yunzhi Technology Co ltd
Priority to CN201810157158.8A priority Critical patent/CN109411041B/zh
Priority to US16/766,569 priority patent/US11350868B2/en
Priority to EP18907280.4A priority patent/EP3698709A4/en
Priority to PCT/CN2018/083465 priority patent/WO2019161610A1/zh
Publication of CN109411041A publication Critical patent/CN109411041A/zh
Application granted granted Critical
Publication of CN109411041B publication Critical patent/CN109411041B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7264Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • A61B5/7267Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/024Detecting, measuring or recording pulse rate or heart rate
    • A61B5/02405Determining heart rate variability
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/024Detecting, measuring or recording pulse rate or heart rate
    • A61B5/0245Detecting, measuring or recording pulse rate or heart rate by using sensing means generating electric signals, i.e. ECG signals
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • A61B5/353Detecting P-waves
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • A61B5/355Detecting T-waves
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • A61B5/358Detecting ST segments
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • A61B5/363Detecting tachycardia or bradycardia
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • A61B5/364Detecting abnormal ECG interval, e.g. extrasystoles, ectopic heartbeats
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • A61B5/366Detecting abnormal QRS complex, e.g. widening
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7221Determining signal validity, reliability or quality
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/725Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271Specific aspects of physiological measurement analysis
    • A61B5/7282Event detection, e.g. detecting unique waveforms indicative of a medical condition
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/74Details of notification to user or communication with user or patient ; user input means
    • A61B5/742Details of notification to user or communication with user or patient ; user input means using visual displays
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H15/00ICT specially adapted for medical reports, e.g. generation or transmission thereof
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2562/00Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
    • A61B2562/08Sensors provided with means for identification, e.g. barcodes or memory chips
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/0002Remote monitoring of patients using telemetry, e.g. transmission of vital signals via a communication network
    • A61B5/0004Remote monitoring of patients using telemetry, e.g. transmission of vital signals via a communication network characterised by the type of physiological signal transmitted
    • A61B5/0006ECG or EEG signals
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • A61B5/352Detecting R peaks, e.g. for synchronising diagnostic apparatus; Estimating R-R interval
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/12Classification; Matching
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/70ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for mining of medical data, e.g. analysing previous cases of other patients

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Cardiology (AREA)
  • Physics & Mathematics (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Biomedical Technology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Artificial Intelligence (AREA)
  • Signal Processing (AREA)
  • Physiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Fuzzy Systems (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

本发明实施例涉及一种心电信息处理方法和心电工作站系统,所述方法包括:心电工作站接收多台心电监测设备输出的心电图数据,对心电图数据进行心电数据分析,生成报告数据并进行存储;接收用户输入的报告数据查询指令,根据用户的用户ID查询相应的报告数据,并生成报告数据查询结果列表数据,进行显示输出;接收用户输入的选择指令,根据选择指令获取被选择的报告数据;接收用户输入的报告数据会诊请求;根据报告数据会诊请求获取用户ID对应的关联用户的用户ID,并根据关联用户的用户ID将报告数据发送至关联用户的用户设备;接收关联用户的用户设备发送的会诊结果反馈数据,并显示。

Description

心电信息处理方法和心电工作站系统
技术领域
本发明涉及数据处理技术领域,尤其涉及一种心电信息处理方法和心电工作站系统。
背景技术
心电工作站是将人体弱小的心电信号由心电采集器采集,送入计算机处理系统放大、分析处理,从而达到诊断心脏疾病的目的,通常由计算机主机、监视器、导联线、心电信号采集盒和打印机等几部分组成,具有视窗菜单,不需要深奥的操作技术和过多的人工干预,临床医生或操作人员即可以通过它获得心电图输出报告。
心电工作站所执行的心电监测是一种常用的临床医疗监测手段。其监测的心电信号是心肌细胞的电活动在体表反映出的微弱电流,通过体表电极和放大记录系统记录下来。
在心电监测中涉及心电图数据自动分析,但由于心电图信号本身的复杂与变异性,很容易受到各种信号干扰,在记录过程中同时还会记录到其他非心源性的电信号,比如骨骼肌活动带来的肌电信号干扰等等。这些信号可能会导致不正确心搏信号检测结果的输出。
此外,心电信号是心肌电活动过程的体现,可以体现出大量的心脏状态的信息。在心脏状态出现问题的时候,心电信号会出现相应的改变,目前自动分析的准确率远远不够,致使输出的心电图检测报告的参考意义不大,还是依然依赖于医生的主观判断来形成心电图检测报告。
此外,现有心电工作站的心电图检测报告输出方式单一,无法满足不同用户或不同使用需求的差异化要求。因而如何有效提高心电图的自动分析水平,并智能化的生成和输出分析数据以及满足多方位的需求,是本发明所要解决的难题和挑战。
发明内容
本发明为实现上述目的,提出了一种心电信息处理方法和心电工作站系统,通过对心电监测设备输出的心电图数据进行完整快速的自动分析,形成报告数据;能够对报告数据进行查询、修改,按照请求或设定的输出模式和数据格式输出报告数据,并能设定关联用户,发起对报告数据的会诊请求,使得关联用户能够获取到报告数据并进行会诊处理。
本发明实施例第一方面提供了一种心电信息处理方法,包括:
心电工作站接收多台心电监测设备输出的心电图数据,所述心电图数据包括被测对象ID、心电监测设备ID和检测时间信息;
对所述心电图数据进行心电数据分析,生成报告数据并进行存储;所述报告数据包括:所述被测对象ID对应的报告结论数据、报告表项数据和/或报告图形数据;
接收用户输入的报告数据查询指令,根据所述用户的用户ID查询相应的报告数据,并生成报告数据查询结果列表数据,进行显示输出;所述报告数据查询结果列表数据至少包括查询得到的一个或多个所述报告数据的被测对象ID、心电监测设备ID和检测时间信息;
接收用户输入的选择指令,根据所述选择指令获取被选择的报告数据;
接收用户输入的报告数据会诊请求;
根据所述报告数据会诊请求获取所述用户ID对应的关联用户的用户ID,并根据所述关联用户的用户ID将所述报告数据发送至所述关联用户的用户设备;
接收所述关联用户的用户设备发送的会诊结果反馈数据,并显示。
优选的,所述根据所述用户的用户ID查询相应的报告数据具体包括:
根据所述用户ID确定所述用户的用户权限信息;
根据所述用户权限信息确定与所述用户ID相关联的一个或多个被测对象ID;
查询所述用户ID相关联的一个或多个被测对象ID的报告数据。
优选的,所述方法还包括:
接收所述用户根据所述被选择的报告数据和/或所述会诊结果反馈数据输入的报告结论修改数据和/或报告表项修改数据;
根据所述报告结论修改数据更新所述报告结论数据,根据所述报告表项修改数据更新所述报告表项数据。
优选的,所述方法还包括:
接收用户输出的报告输出指令,根据所述报告输出指令确定输出模式信息和输出格式信息;
根据所述输出模式信息选择性输出所述报告结论数据、报告表项数据和/或报告图形数据中的部分数据或全部数据;
根据所述输出格式信息,将所述部分数据或全部数据进行数据格式转换,生成报告输出数据。
优选的,在根据所述报告数据会诊请求获取所述用户ID对应的关联用户的用户ID之前,所述方法还包括:
配置各用户的级别属性信息;
根据所述用户ID的级别属性信息确定所述用户ID对应的关联用户的用户ID,并生成关联用户记录数据;
存储所述关联用户记录数据。
优选的,所述心电图数据为单导联或多导联的时间序列数据;所述对所述心电图数据进行心电数据分析,生成报告数据具体包括:
将所述心电图数据的数据格式经过重采样转换为预设标准数据格式,并对转换后的预设标准数据格式的心电图数据进行第一滤波处理;
对所述第一滤波处理后的心电图数据进行心搏检测处理,识别所述心电图数据包括的多个心搏数据;每个所述心搏数据对应一个心搏周期,包括相应的P波、QRS波群、T波的幅值数据和起止时间数据;
根据训练得到的干扰识别二分类模型对所述心搏数据进行干扰识别,确定所述心搏数据是否存在干扰噪音,以及用于判断干扰噪音的一个概率值;
根据心搏数据的导联参数和所述心搏数据,基于所述干扰识别的结果和时间规则合并生成心搏时间序列数据;根据所述心搏时间序列数据生成心搏分析数据;
对所述心搏分析数据进行信号质量分析,评估所述心搏分析数据的信号质量,得到所述心搏分析数据的信号质量评价参数;
根据训练得到的心搏分类模型对所述心搏分析数据进行幅值和时间表征数据的特征提取和分析,得到所述心搏分析数据的一次分类信息;
对所述一次分类信息结果中的特定心搏的心搏分析数据输入到训练好的ST段和T波改变模型进行识别,确定ST段和T波评价信息;
根据心搏时间序列数据,对心搏分析数据进行P波和T波特征检测,确定每个心搏中P波和T波的详细特征信息;
对所述心搏分析数据在所述一次分类信息下根据所述心电图基本规律参考数据、P波和T波的详细特征信息以及所述ST段和T波评价信息进行二次分类处理,得到心搏分类信息;
将所述心搏分析数据根据所述心搏分类信息和心电图基本规律参考数据生成心电图事件数据;根据所述信号质量评价参数对心电图事件数据进行筛选,得到所述报告结论数据和报告表项数据;根据每种心电图事件数据中的典型数据片段,生成所述报告图形数据。
进一步优选的,对所述第一滤波处理后的心电图数据进行心搏检测处理还包括:
根据所述QRS波群确定RR间期并计算所述RR间期内噪声的估计值;
根据所述噪声的估计值和每个QRS波群中的最大幅值确定每个QRS波群的检测置信度。
优选的,所述根据训练得到的干扰识别二分类模型对所述心搏数据进行干扰识别,具体包括:
以第一数据量对所述心搏数据进行切割采样,并将切割采样得到的数据输入到干扰识别二分类模型进行干扰识别;
识别心搏数据中,心搏间期大于等于预设间期判定阈值的数据片段;
对所述心搏间期大于等于预设间期判定阈值的数据片段进行信号异常判断,确定是否为异常信号;
如果不是异常信号,则以预设时间宽度,根据设定时值确定所述数据片段中滑动取样的起始数据点和终止数据点,并由所述起始数据点开始对所述数据片段进行滑动取样,至所述终止数据点为止,得到多个取样数据段;
对每个所述取样数据段进行所述干扰识别。
进一步优选的,所述对所述一次分类信息结果中的特定心搏的心搏分析数据输入到训练好的ST段和T波改变模型进行识别,确定ST段和T波评价信息,具体包括:
将所述一次分类信息中的特定心搏数据,按导联依次输入到根据训练得到的ST段和T波改变模型,对各导联的特定心搏数据进行幅值和时间表征数据的特征提取和分析,得到各导联的ST段和T波改变信息,确定ST段和T波评价信息;所述ST段和T波评价信息具体为所述心搏片段数据对应的ST段和T波发生改变的导联位置信息。
进一步优选的,所述根据每种心电图事件数据中的典型数据片段,生成报告图形数据,具体包括:
根据所述信号质量评价参数,对每种心电图事件的数据片段进行评估,选择信号质量评价参数最高的数据片段为所述心电图事件中的典型数据片段。
本发明实施例提供的心电信息处理方法,通过对心电监测设备输出的心电图数据进行完整快速的自动分析,形成报告数据;能够对报告数据进行查询、修改,按照请求或设定的输出模式和数据格式输出报告数据,并能设定关联用户,发起对报告数据的会诊请求,使得关联用户能够获取到报告数据并进行会诊处理。
本发明实施例第二方面提供了一种心电工作站系统,该系统包括第一方面所述的心电工作站和多台心电监测设备;
其中,所述心电工作站包括存储器和处理器,存储器用于存储程序,处理器用于执行第一方面及第一方面的各实现方式中的方法。
本发明实施例第三方面提供了一种包含指令的计算机程序产品,当计算机程序产品在计算机上运行时,使得计算机执行第一方面及第一方面的各实现方式中的方法。
本发明实施例第四方面提供了一种计算机可读存储介质,计算机可读存储介质上存储有计算机程序,计算机程序被处理器执行时实现第一方面及第一方面的各实现方式中的方法。
附图说明
图1为本发明实施例提供的心电信息处理方法的流程图;
图2为本发明实施例提供的心电图数据的处理方法的流程图;
图3为本发明实施例提供的干扰识别二分类模型的示意图;
图4为本发明实施例提供的心搏分类模型的示意图;
图5为本发明实施例提供的ST段和T波改变模型的示意图;
图6为本发明实施例提供的心电工作站系统的结构示意图;
图7a-图7c为本发明实施例提供的一种静态心电监测数据的报告输出模式;
图8a-图8g为本发明实施例提供的一种动态心电监测数据的报告输出模式。
具体实施方式
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
心电工作站是将采集到的人体弱小的心电信号送入计算机处理系统放大、分析处理,从而达到心电检查、监护的目的。心电信号是心肌细胞的电活动在体表反映出的微弱电流,通过体表电极和放大记录系统记录下来。
本发明就是基于记录的心电信号,进行完整快速的自动分析形成报告数据;并建立了对报告数据进行查询、按请求或设定的输出模式和数据格式输出、对报告数据的更新、修改,以及请求报告数据的会诊处理的机制。
下面首先结合图1所示的心电信息处理方法的流程图,对本发明的心电信息处理方法进行详述。
如图1所示,本发明的心电信息处理方法主要包括如下步骤:
步骤110,心电工作站接收多台心电监测设备输出的心电图数据;
具体的,心电监测设备具体可以是单导联或多导联的心电监测仪,每台心电检测设备都有唯一的设备ID,从而能够对测得的数据进行标记,以表明该数据是由哪台设备获得的。
心电工作站在本发明中是指具有心电图数据分析处理功能的且能够提供用户进行信息交互的设备,能够实现数据传输、用户管理和数据管理以及心电图数据的报告输出管理等功能。
它与心电监测设备可以集成于一体,也可以是二者通过有线方式或无线网络进行连接进行数据传输。其中无线网络包括无线但不限于基于IEEE 802.11b标准的无线局域网(WIFI),蓝牙,3G/4G/5G移动通信网络,物联网等方式。在本实施例中的心电监测设备,侧重于心电信号采集的功能实现,输出心电图数据由心电工作站进行处理。
心电监测设备通过无创心电图检查的方式对心脏细胞电生理活动产生的信号以单导联或多导联的形式进行采集记录,得到心电图数据发送到心电工作站。心电图数据中包括心电图数据包括被测对象ID、心电监测设备ID和检测时间信息;
步骤120,对心电图数据进行心电数据分析,生成报告数据并进行存储;
考虑到在心电信号的记录过程中同时还会记录到其他非心源性的电信号,比如骨骼肌活动带来的肌电信号干扰等等,因此我们认为需要对心电信号进行有效的干扰识别和排除,才能够有效降低因为干扰信号造成的误报。
此外,心电信号是心肌电活动过程的体现,因此心电信号除了可以用来检测心率以外,还可以体现出大量的心脏状态的信息。在心脏状态出现问题的时候,心电信号会出现相应的改变。在对业内现有的心电信号处理方法的研究中我们发现,目前只对心电信号进行了非常有限的分析和报警。对此,除了对心电信号进行有效的干扰识别和排除,以降低因为干扰信号造成的误报之外,我们认为还可以从以下几点进行改进:
第一,在心搏特征提取中需要对P波、T波进行准确识别,可以避免心搏检测的多检和漏检,比如对一些特殊心电图信号,例如心律比较缓慢患者的高大T波,或者T波肥大的信号的多检。
第二,对心搏的分类进行更加细致的划分,而不能仅停留在窦性、室上性和室性这三种分类,从而满足临床心电图医生复杂全面的分析要求。
第三,准确识别房扑房颤和ST-T改变,从而能够有助于提供对ST段和T波改变对心肌缺血分析的帮助。
第四,对心搏和心电事件的准确识别。
在本发明中,我们针对上述几点,通过对心电图数据的分析计算,特别是引入人工智能(AI)技术,对采集的数字信号进行心律失常分析、长间歇停搏,扑动和颤动,传导阻滞,早搏和逸搏,心动过缓,心动过快,ST-T改变检测、心电事件的分析与归类,以达到产生准确报警信号的目的,从而有效的进行病人生命体征的监护。
基于上述几点,本发明的心电图数据的处理过程,采用了基于人工智能自学习的心电图自动分析方法,是基于人工智能卷积神经网络(CNN)模型来实现的。CNN模型是深度学习中的监督学习方法,就是一个模拟神经网络的多层次网络(隐藏层hidden layer)连接结构,输入信号依次通过每个隐藏层,在其中进行一系列复杂的数学处理(Convolution卷积、Pooling池化、Regularization正则化、防止过拟合、Dropout暂时丢弃、Activation激活、一般使用Relu激活函数),逐层自动地抽象出待识别物体的一些特征,然后把这些特征作为输入再传递到高一级的隐藏层进行计算,直到最后几层的全连接层(FullConnection)重构整个信号,使用Softmax函数进行逻辑(logistics)回归,达到多目标的分类。
CNN属于人工智能中的监督学习方法,在训练阶段,输入信号经过多个的隐藏层处理到达最后的全连接层,softmax逻辑回归得到的分类结果,与已知的分类结果(label标签)之间会有一个误差,深度学习的一个核心思想就是通过大量的样本迭代来不断地极小化这个误差,从而计算得到连接各隐藏层神经元的参数。这个过程一般需要构造一个特别的损失函数(costfunction),利用非线性优化的梯度下降算法和误差反向传播算法(backpropagation algorithm,BP),快速有效地极小化整个深度(隐藏层的层数)和广度(特征的维数)都十分复杂的神经网络结构中所有连接参数。
深度学习把需要识别的数据输入到训练模型,经过第一隐藏层、第二隐藏层、第三隐藏层,最后是输出识别结果。
在本发明中,对心电图数据进行波群特征识别、干扰识别、心搏分类等都是基于人工智能自学习的训练模型来得到输出结果,分析速度快,准确程度高。
具体的,该步骤通过对心电图数据进行波群特征识别,得到心电图数据的特征信号,根据特征信号对心电图数据进行心搏分类,结合心电图基本规律参考数据得到心搏分类信息,根据心搏分类信息处理生成心电图事件数据,并最终生成报告数据来实现。进一步的,该过程具体可以通过如图2所示的如下步骤来实现。
步骤1201,将心电图数据的数据格式经过重采样转换为预设标准数据格式,并对转换后的预设标准数据格式的心电图数据进行第一滤波处理;
具体的,心电图数据的格式适配读取,对不同的设备有不同的读取实现,读取后,需要调整基线、根据增益转换成毫伏数据。经过数据重采样,把数据转换成全流程能够处理的采样频率。然后通过滤波去除高频,低频的噪音干扰和基线漂移,提高人工智能分析准确率。将处理后的心电图数据以预设标准数据格式保存。
通过本步骤解决不同在使用不同导联,采样频率和传输数据格式的差异,以及通过数字信号滤波去除高频,低频的噪音干扰和基线漂移。
数字信号滤波可以分别采用高通滤波器,低通滤波器和中值滤波,把工频干扰、肌电干扰和基线漂移干扰消除,避免对后续分析的影响。
更具体的,可以采用低通、高通巴特沃斯滤波器进行零相移滤波,以去除基线漂移和高频干扰,保留有效的心电信号;中值滤波则可以利用预设时长的滑动窗口内数据点电压幅值的中位数替代窗口中心序列的幅值。可以去除低频的基线漂移。
步骤1202,对第一滤波处理后的心电图数据进行心搏检测处理,识别心电图数据包括的多个心搏数据;
每个心搏数据对应一个心搏周期,包括相应的P波、QRS波群、T波的幅值和起止时间数据。本步骤中的心搏检测由两个过程构成,一是信号处理过程,从所述第一滤波处理后的心电图数据中提取QRS波群的特征频段;二是通过设置合理的阈值确定QRS波群的发生时间。在心电图中,一般会包含P波、QRS波群、T波成分以及噪声成分。一般QRS波群的频率范围在5到20Hz之间,可以通过一个在此范围内的带通滤波器提出QRS波群信号。然而P波、T波的频段以及噪声的频段和QRS波群频段有部分重叠,因此通过信号处理的方法并不能完全去除非QRS波群的信号。因此需要通过设置合理的阈值来从信号包络中提取QRS波群位置。具体的检测过程是一种基于峰值检测的过程。针对信号中每一个峰值顺序进行阈值判断,超过阈值时进入QRS波群判断流程,进行更多特征的检测,比如RR间期、形态等。
在心电信息的记录过程中心搏信号的幅度和频率时时刻刻都在变化,并且在疾病状态下,这种特性会表现的更强。在进行阈值设定时,需要根据数据特征在时域的变化情况动态的进行阈值调整。为了提高检测的准确率和阳性率,QRS波群检测大多采用双幅度阈值结合时间阈值的方式进行,高阈值具有更高的阳性率,低阈值具有更高的敏感率,在RR间期超过一定时间阈值,使用低阈值进行检测,减少漏检情况。而低阈值由于阈值较低,容易受到T波、肌电噪声的影响,容易造成多检,因此优先使用高阈值进行检测。
对于不同导联的心搏数据都具有导联参数,用以表征该心搏数据为哪个导联的心搏数据。因此在得到心电图数据的同时也就可以根据其传输来源确定了其导联的信息,将此信息作为心搏数据的导联参数。
步骤1203,根据心搏数据确定每个心搏的检测置信度;
具体的,置信度计算模块在心搏检测的过程中,根据QRS波群的幅度以及RR间期内噪声信号的幅度比例可以提供针对QRS波群检测置信度的估计值。
步骤1204,根据干扰识别二分类模型对心搏数据进行干扰识别,得到心搏数据是否存在干扰噪音,以及用于判断干扰噪音的一个概率值;
因为在长时间记录过程中易受多种影响出现干扰现象,导致获取的心搏数据无效或不准确,不能正确反映受测者的状况,同时也增加医生诊断难度及工作量;而且干扰数据也是导致智能分析工具无法有效工作的主要因素。因此,将外界信号干扰降到最低显得尤为重要。
本步骤基于以深度学习算法为核心的端到端二分类识别模型,具有精度高,泛化性能强的特点,可有效地解决电极片脱落、运动干扰和静电干扰等主要干扰来源产生的扰动问题,克服了传统算法因干扰数据变化多样无规律而导致的识别效果差的问题。
具体可以通过如下方法来实现:
步骤A,对心搏数据使用干扰识别二分类模型进行干扰识别;
步骤B,识别心搏数据中,心搏间期大于等于预设间期判定阈值的数据片段;
步骤C,对心搏间期大于等于预设间期判定阈值的数据片段进行信号异常判断,确定是否为异常信号;
其中,异常信号的识别主要包括是否为电极片脱落、低电压等情况。
步骤D,如果不是异常信号,则以预设时间宽度,根据设定时值确定数据片段中滑动取样的起始数据点和终止数据点,并由起始数据点开始对数据片段进行滑动取样,至终止数据点为止,得到多个取样数据段;
步骤E,对每个取样数据段进行干扰识别。
以一个具体的例子对上述步骤A-E进行说明。对每个导联的心搏数据以设定的第一数据量进行切割采样,然后分别输入到干扰识别二分类模型进行分类,获得干扰识别结果和对应结果的一个概率值;对心搏间期大于等于2秒的心搏数据,先判断是否是信号溢出,低电压,电极脱落;如果不是上述情况,就按照第一数据量,从左边心搏开始,向右连续以第一数据量不重叠滑动取样,进行识别。
输入可以是任一导联的第一数据量心搏数据,然后采用干扰识别二分类模型进行分类,直接输出是否为干扰的分类结果,获得结果快,精确度高,稳定性好,可为后续分析提供更有效优质的数据。
因为干扰数据往往是由外界扰动因素的作用而引起的,主要有电极片脱落、低电压、静电干扰和运动干扰等情况,不但不同扰动源产生的干扰数据不同,而且相同扰动源产生的干扰数据也是多种多样;同时考虑到干扰数据虽然多样性布较广,但与正常数据的差异很大,所以在收集干扰的训练数据时也是尽可能的保证多样性,同时采取移动窗口滑动采样,尽可能增加干扰数据的多样性,以使模型对干扰数据更加鲁棒,即使未来的干扰数据不同于以往任何的干扰,但相比于正常数据,其与干扰的相似度也会大于正常数据,从而使模型识别干扰数据的能力增强。
本步骤中采用的干扰识别二分类模型可以如图3所示,网络首先使用2层卷积层,卷积核大小是1x5,每层后加上一个最大值池化。卷积核数目从128开始,每经过一次最大池化层,卷积核数目翻倍。卷积层之后是两个全连接层和一个softmax分类器。由于该模型的分类数为2,所以softmax有两个输出单元,依次对应相应类别,采用交叉熵做为损失函数。
对于该模型的训练,我们采用了来源于30万病人近400万精确标注的数据片段。标注分为两类:正常心电图信号或者是有明显干扰的心电图信号片段。我们通过定制开发的工具进行片段标注,然后以自定义标准数据格式保存干扰片段信息。
在训练过程,使用两台GPU服务器进行几十次轮循训练。在一个具体的例子中,采样率是200Hz,数据长度是300个心电图电压值(毫伏)的一个片段D[300],输入数据是:InputData(i,j),其中,i是第i个导联,j是导联i第j个片段D。输入数据全部经过随机打散才开始训练,保证了训练过程收敛,同时,控制从同一个病人的心电图数据中收集太多的样本,提高模型的泛化能力,既真实场景下的准确率。训练收敛后,使用100万独立的测试数据进行测试,准确率可以到达99.3%。另有具体测试数据如下表1。
干扰 正常
敏感率(Sensitivity) 99.14% 99.32%
阳性预测率(Positive Predicitivity) 96.44% 99.84%
表1
步骤1205,根据检测置信度确定心搏数据的有效性,并且,根据确定有效的心搏数据的导联参数和心搏数据,基于干扰识别的结果和时间规则合并生成心搏时间序列数据,并根据心搏时间序列数据生成心搏分析数据;
具体的,由于心电图信号的复杂性以及每个导联可能受到不同程度的干扰影响,依靠单个导联检测心搏会存在多检和漏检的情况,不同导联检测到心搏结果的时间表征数据没有对齐,所以需要对所有导联的心搏数据根据干扰识别结果和时间规则进行合并,生成一个完整的心搏时间序列数据,统一所有导联心搏数据的时间表征数据。其中,时间表征数据用于表示每个数据点在心电图数据信号时间轴上的时间信息。根据这个统一的心搏时间序列数据,在后续的分析计算时,可以使用预先设置好的阀值,对各导联心搏数据进行切割,从而生成具体分析需要的各导联的心搏分析数据。
上述每个导联的心搏数据在合并前,需要根据步骤1203中获得的检测置信度确定心搏数据的有效性。
具体的,导联心搏合并模块执行的心搏数据合并过程如下:根据心电图基本规律参考数据的不应期获取不同导联心搏数据的时间表征数据组合,丢弃其中偏差较大的心搏数据,对上述时间表征数据组合投票产生合并心搏位置,将合并心搏位置加入合并心搏时间序列,移动到下一组待处理的心搏数据,循环执行直至完成所有心搏数据的合并。
其中,心电图活动不应期可以优选在200毫秒至280毫秒之间。获取的不同导联心搏数据的时间表征数据组合应满足以下条件:心搏数据的时间表征数据组合中每个导联最多包含一个心搏数据的时间表征数据。在对心搏数据的时间表征数据组合进行投票时,使用检出心搏数据的导联数占有效导联数的百分比来决定;若心搏数据的时间表征数据对应导联的位置为低电压段、干扰段以及电极脱落时认为该导联对此心搏数据为无效导联。在计算合并心搏具体位置时,可以采用心搏数据的时间表征数据平均值得到。在合并过程中,本方法设置了一个不应期来避免错误合并。
在本步骤中,通过合并操作输出一个统一的心搏时间序列数据。该步骤同时能够降低心搏的多检率和漏检率,有效的提高心搏检测的敏感度和阳性预测率。
步骤1206,对心搏分析数据进行信号质量分析,评估心搏分析数据的信号质量,得到心搏分析数据的信号质量评价参数。
如果监测的心电图数据是静态心电图数据,那么信号质量会相对较好,有效信号的比例比较高,这种情况下可以不用进行信号质量评估。而如果监测的是动态心电图数据,则在监测过程中会发生各种扰动,导致无效信号的产生会相对较多,为了避免无效信号影响最终的结果计算,需要进行信号质量评估,降低无效信号对最终结果的影响。
具体的,通过提取心搏分析数据中的RR间期信号,对RR间期信号进行第二滤波处理和包络计算,确定噪声信号强度,并计算相对于相应的心搏时间序列数据的最大振幅的信号噪声比,确定信号质量评价参数。其具体步骤可以包括:
步骤A,提取心搏分析数据中QRS波群的位置信息和宽度信息;
步骤B,提取相邻两个QRS波群信号之间的RR间期的信号;
步骤C,对RR间期的信号进行滤波处理,并对滤波处理后的信号进行包络计算,得到RR间期的噪声信号的平均功率;其中,RR间期的噪声信号的平均功率为RR间期内包络幅度的均值。
步骤D,根据噪声信号的平均功率和QRS波群信号的功率得到信号质量评价参数;
进一步具体的,信号质量评价参数可以由RR间期内相对于QRS波群的噪声水平来表示。具体通过根据QRS波群的功率和噪声信号的平均功率计算得到。
以公式表示为:
Figure BDA0001581773380000151
其中,Si为第i个QRS波群信号的幅度,Ni,t为第i个RR间期区间内第t个样本点的幅值,T为RR间期长度。
步骤1207,根据心搏分类模型对心搏分析数据进行幅值和时间表征数据的特征提取和分析,得到心搏分析数据的一次分类信息;
具体的,该步骤在步骤1205之后执行。不同心电监测设备在信号测量、采集或者输出的导联数据等方面存在的差异,因此可以根据具体情况,使用简单的单导联分类方法,或者是多导联分类方法。多导联分类方法又包括导联投票决策分类方法和导联同步关联分类方法两种。导联投票决策分类方法是基于各导联的心搏分析数据进行导联独立分类,再把结果投票融合确定分类结果的投票决策方法;导联同步关联分类方法则采用对各导联的心搏分析数据进行同步关联分析的方法。单导联分类方法就是对单导联设备的心搏分析数据,直接使用对应导联模型进行分类,没有投票决策过程。下面对以上所述几种分类方法分别进行说明。
单导联分类方法包括:
根据心搏时间序列数据,将单导联心搏数据进行切割生成单导联的心搏分析数据,并输入到训练得到的对应该导联的心搏分类模型进行幅值和时间表征数据的特征提取和分析,得到单导联的一次分类信息。
导联投票决策分类方法可以具体包括:
第一步、根据心搏时间序列数据,对各导联心搏数据进行切割,从而生成各导联的心搏分析数据;
第二步、根据训练得到的各导联对应的心搏分类模型对各导联的心搏分析数据进行幅值和时间表征数据的特征提取和分析,得到各导联的分类信息;
第三步、根据各导联的分类信息和导联权重值参考系数进行分类投票决策计算,得到所述一次分类信息。具体的,导联权重值参考系数是基于心电数据贝叶斯统计分析得到各导联对不同心搏分类的投票权重系数。
导联同步关联分类方法可以具体包括:
根据心搏时间序列数据,对各导联心搏数据进行切割,从而生成各导联的心搏分析数据;然后根据训练得到的多导联同步关联分类模型对各导联的心搏分析数据进行同步幅值和时间表征数据的特征提取和分析,得到心搏分析数据的一次分类信息。
心搏数据的同步关联分类方法输入是动态心电图设备所有导联数据,按照心搏分析数据统一的心搏位点,截取各导联上相同位置和一定长度的数据点,同步输送给经过训练的人工智能深度学习模型进行计算分析,输出是每个心搏位置点综合考虑了所有导联心电图信号特征,以及心搏在时间上前后关联的心律特征的准确心搏分类。
本方法充分考虑了心电图不同导联数据实际上就是测量了心脏电信号在不同的心电轴向量方向传递的信息流,把心电图信号在时间和空间上传递的多维度数字特征进行综合分析,极大地改进了传统方法仅仅依靠单个导联独立分析,然后把结果汇总进行一些统计学的投票方式而比较容易得出的分类错误的缺陷,极大地提高了心搏分类的准确率。
本步骤中采用的心搏分类模型可以如图4所示,具体可以为基于人工智能深度学习的卷积神经网络AlexNet,VGG16,Inception等模型启发的端对端多标签分类模型。具体的讲,该模型的网络是一个7层的卷积网络,每个卷积之后紧跟一个激活函数。第一层是两个不同尺度的卷积层,之后是六个卷积层。七层卷积的卷积核分别是96,256,256,384,384,384,256。除第一层卷积核有两个尺度分别是5和11外,其他层卷积核尺度为5。第三、五、六、七层卷积层后是池化层。最后跟着两个全连接层。
本步骤中的心搏分类模型,我们采用了训练集包含30万病人的1700万数据样本进行训练。这些样本是根据动态心电图分析诊断的要求对数据进行准确的标注产生的,标注主要是针对常见心律失常,传导阻滞以及ST段和T波改变,可满足不同应用场景的模型训练。具体以预设标准数据格式保存标注的信息。在训练数据的预处理上,为增加模型的泛化能力,对于样本量较少的分类做了小幅的滑动来扩增数据,具体的说,就是以每个心搏为基础,按照一定步长(比如10-50个数据点)移动2次,这样就可以增加2倍的数据,提高了对这些数据量比较少的分类样本的识别准确率。经过实际结果验证,泛化能力也得到了改善。
在一个实际训练过程使用了两台GPU服务器进行几十次轮循训练,训练收敛后,使用500万独立的测试数据进行测试,准确率可以到达91.92%。
其中,训练数据的截取的长度,可以是1秒到10秒。比如采样率是200Hz,以2.5s为采样长度,取得的数据长度是500个心电图电压值(毫伏)的一个片段D[500],输入数据是:InputData(i,j),其中,i是第i个导联,j是导联i第j个片段D。输入数据全部经过随机打散才开始训练,保证了训练过程收敛,同时,控制从同一个病人的心电图数据中收集太多的样本,提高模型的泛化能力,既真实场景下的准确率。训练时候,同步输入所有导联对应的片段数据D,按照图像分析的多通道分析方法,对每个时间位置的多个空间维度(不同心电轴向量)的导联数据进行同步学习,从而得到一个比常规算法更准确的分类结果。
步骤1208,对一次分类信息结果中的特定心搏的心搏分析数据输入到ST段和T波改变模型进行识别,确定ST段和T波评价信息;
ST段和T波评价信息具体为心搏分析数据对应的ST段和T波发生改变的导联位置信息。因为临床诊断要求对于ST段和T波的改变定位到具体的导联。
其中,一次分类信息的特定心搏数据是指包含窦性心搏(N)和其它可能包含ST改变的心搏类型的心搏分析数据。
ST段和T波改变导联定位模块将一次分类信息的特定心搏数据,按照每个导联依次输入到一个为识别ST段和T波改变的人工智能深度学习训练模型,进行计算分析,输出的结果说明导联片段的特征是否符合ST段和T波改变的结论,这样就可以确定ST段和T波改变发生的在具体那些导联的信息,即ST段和T波评价信息。具体方法可以是:把一次分类信息中结果是窦性心搏的各导联心搏分析数据,输入给ST段和T波改变模型,对窦性心搏分析数据进行逐一识别判断,以确定窦性心搏分析数据是否存在ST段和T波改变特征以及发生的具体导联位置信息,确定ST段和T波评价信息。
本步骤中采用的ST段和T波改变模型可以如图5所示,具体可以为基于人工智能深度学习的卷积神经网络AlexNet和VGG16等模型启发的端对端分类模型。具体的讲,该模型是一个7层的网络,模型包含了7个卷积,5个池化和2个全连接。卷积使用的卷积核均为1x5,每层卷积的滤波器个数各不相同。第1层卷积滤波器个数为96;第2层卷积和第3层卷积连用,滤波器个数为256;第4层卷积和第5层卷积连用,滤波器个数为384;第6层卷积滤波器个数为384;第7层卷积滤波器个数为256;第1、3、5、6、7层卷积层后是池化。随后是两个全连接,最后还采用Softmax分类器将结果分为两类。为了增加模型的非线性,提取数据更高维度的特征,故采用两个卷积连用的模式。
因为带有ST段和T波改变的心搏在所有心搏中的占比较低,为了兼顾训练数据的多样性及各个类别数据量的均衡性,选取无ST段和T波改变以及有ST段和T波改变的训练数据比例约为2:1,保证了模型在分类过程中良好的泛化能力且不出现对训练数据占比较多一类的倾向性。由于心搏的形态多种多样,不同个体表现的形态不尽相同,因此,为了模型更好估计各分类的分布,能有效提取特征,训练样本从不同年龄,体重,性别和居住地区的个体收集;另外,因为单个个体在同一时间段内的心电图数据往往是高度相似的,所以为了避免过度学习,在获取单个个体的数据时,从所有数据中随机选取不同时间段的少量样本;最后,由于患者的心搏形态存在个体间差异大,而个体内相似度高的特点,因而在划分训练、测试集时,把不同的患者分到不同的数据集,避免同一个体的数据同时出现在训练集与测试集中,由此,所得模型测试结果最接近真实应用场景,保证了模型的可靠性和普适性。
步骤1209,根据心搏时间序列数据,对心搏分析数据进行P波和T波特征检测,确定每个心搏中P波和T波的详细特征信息;
具体的,详细特征信息包括幅值、方向、形态和起止时间的数据;在对心搏信号的分析中,P波、T波以及QRS波中的各项特征也是心电图分析中的重要依据。
在P波和T波特征检测模块中,通过计算QRS波群中切分点位置,以及P波和T波的切分点位置,来提取P波、T波以及QRS波群中的各项特征。可以分别通过QRS波群切分点检测、单导联PT检测算法和多导联PT检测算法来实现。
QRS波群切分点检测:根据QRS波群检测算法提供的QRS波群段功率最大点以及起止点,寻找单个导联中QRS波群的R点,R’点,S点以及S’点。在存在多导联数据时,计算各个切分点的中位数作为最后的切分点位置。
单导联P波、T波检测算法:P波和T波相对QRS波群幅度低、信号平缓,容易淹没在低频噪声中,是检测中的难点。本方法依据QRS波群检测的结果,在消除QRS波群对低频频段的影响后,使用低通滤波器对信号进行第三滤波,使PT波相对幅度增高。之后通过峰值检测的方法在两个QRS波群之间寻找T波。因为T波是心室复极产生的波群,因此T波和QRS波群之间有明确的锁时关系。以检测到的QRS波群为基准,在每个QRS波群到下一个QRS波群间期取中点(比如限制在第一个QRS波群后400ms到600ms之间的范围)作为T波检测结束点,在此区间内选取最大的峰作为T波。再在剩余的峰值内选择幅度最大的峰为P波。同时也根据P波和T波的峰值与位置数据,确定P波和T波的方向和形态特征。优选的,低通滤波的截止频率设置为10-30Hz之间。
多导联P波、T波检测算法:在多导联的情况中,由于心搏中各个波的产生时间相同,空间分布不同,而噪声的时间空间分布不同,可以通过溯源算法来进行P、T波的检测。首先对信号进行QRS波群消除处理并使用低通滤波器对信号进行第三滤波以去除干扰。之后通过独立成分分析算法计算原始波形中的各个独立成分。在分离出的各个独立成分中,依据其峰值的分布特征以及QRS波群位置,选取相应的成分作为P波和T波信号,同时确定P波和T波的方向和形态特征。
步骤1210,对心搏分析数据在一次分类信息下根据心电图基本规律参考数据、P波和T波的详细特征信息以及ST段和T波评价信息进行二次分类处理,得到心搏分类信息;
具体的,心电图基本规律参考数据是遵循权威心电图教科书中对心肌细胞电生理活动和心电图临床诊断的基本规则描述生成的,比如两个心搏之间最小的时间间隔,P波与R波的最小间隔等等,用于将心搏分类后的一次分类信息再进行细分;主要根据是心搏间RR间期以及不同心搏信号在各导联上的医学显著性;心搏审核模块依据心电图基本规律参考数据结合一定连续多个心搏分析数据的分类识别,以及P波和T波的详细特征信息将室性心搏分类拆分更细的心搏分类,包括:室性早搏(V)、室性逸搏(VE)、加速性室性早搏(VT),将室上性类心搏细分为室上性早搏(S)、房性逸搏(SE)、交界性逸搏(JE)和房性加速性早搏(AT)等等。
此外,通过二次分类处理,还可以纠正一次分类中发生的不符合心电图基本规律参考数据的错误分类识别。将细分后的心搏分类按照心电图基本规律参考数据进行模式匹配,找到不符合心电图基本规律参考数据的分类识别,根据RR间期及前后分类标识纠正为合理的分类。
具体的,经过二次分类处理,可以输出多种心搏分类,比如:正常窦性心搏(N)、完全性右束支阻滞(N_CRB)、完全性左束支阻滞(N_CLB)、室内阻滞(N_VB)、一度房室传导阻滞(N_B1)、预激(N_PS)、室性早搏(V)、室性逸搏(VE)、加速性室性早搏(VT)、室上性早搏(S)、房性逸搏(SE)、交界性逸搏(JE)、加速性房性早搏(AT)、房扑房颤(AF)、伪差(A)等分类结果。
通过本步骤,还可以完成基础心率参数的计算。其中基础计算的心率参数包括:RR间期、心率、QT时间、QTc时间等参数。
随后,根据心搏二次分类结果,按照心电图基本规律参考数据进行模式匹配,可以得到分类对应于心电图事件数据的以下这些典型的心电图事件,包括但不限于:
室上性早搏
室上性早搏成对
室上性早搏二联律
室上性早搏三联律
房性逸搏
房性逸搏心律
交界性逸搏
交界性逸搏心律
非阵发性室上性心动过速
最快室上性心动过速
最长室上性心动过速
室上性心动过速
短阵室上性心动过速
心房扑动-心房颤动
室性早搏
室性早搏成对
室性早搏二联律
室性早搏三联律
室性逸搏
室性逸搏心律
加速性室性自主心律
最快室性心动过速
最长室性心动过速
室性心动过速
短阵室性心动过速
二度I型窦房传导阻滞
二度II型窦房传导阻滞
一度房室传导阻滞
二度I型房室传导阻滞
二度II型房室传导阻滞
二度II型(2:1)房室传导阻滞
高度房室传导阻滞
完全性左束支阻滞
完全性右束支阻滞
室内阻滞
预激综合症
ST段和T波改变
最长RR间期
步骤1211,将心搏分析数据根据心搏分类信息和心电图基本规律参考数据生成心电图事件数据;根据信号质量评价参数对心电图事件数据进行筛选,得到报告结论数据和报告表项数据;根据每种心电图事件数据中的典型数据片段,生成报告图形数据。
根据参考信号质量评价参数,剔除心电图事件数据中不可信事件,计算各种心率参数,统计事件的数量和发生时间等,得到报告结论数据和报告表项数据。
具体的,心率参数计算包括计算平均心率、及最大、最小心率等。在计算最大和最小心率时,以固定长度片段为统计单位,逐个心搏进行全程扫描和统计比较。片段的长度一般为8-10秒,也可以根据需要自由设定。计算心率时,针对窦性心率占主体和异位心率占主体采用不同的心搏类型统计计算方法。在计算最大最小心率时:针对窦性心率占主体的心电图,仅计算窦性分类心搏;针对房扑房颤类占主体的心电图,仅计算房扑房颤类心搏及窦性心搏;针对非房扑房颤的其他异位心搏占主体的心电图,除伪差外所有类型心搏均参与计算。
对心电图事件数据,根据信号质量评价参数进行质量评估,选取数据信号质量最高的事件片段,同时分析片段中的包含的事件种类数量,优先选取最具代表性的片段,优选片段的开始位置,尽量保证事件心搏位于所选取片段的中部,生成报告图形数据。
在优选的方案中,心电图事件的选择规则可以具体如下例子所述:
对于一般心电图事件,选取单个片段。当片段第一个心搏心率大于等于100时,片段起始点位置距离第一个心搏距离为0.3秒。当片段第一个心搏心率小于等于45时,片段起始点位置距离第一个心搏距离为0.37秒。
对于有多个片段的心电图事件,需要对多个片段进行质量评估,依据信号质量评价参数计算截取片段的无干扰信号占比,当无干扰信号占比达到设定的阈值(优选的在60%-95%范围中确定阈值),片段符合筛选条件,从中获得最优典型数据段。
根据被测对象ID,报告结论数据、报告表项数据和报告图形数据相关联,存储为报告数据。因此在查询调用时,可以根据被测对象ID直接获取到该被测对象ID相关的各种报告数据。报告数据中还相应的包括有报告数据的检测时间信息。这个信息应当至少包括心电图数据的获取时间/时间段的信息,还可以包括报告数据最后更新的时间。因为考虑到可能对报告数据进行修改,所以也可以对该时间进行记录。
步骤130,接收用户输入的报告数据查询指令,根据用户的用户ID查询相应的报告数据,并生成报告数据查询结果列表数据,进行显示输出;
具体的,心电工作站能够向用户提供人机交互界面,通过人机交互界面,用户可以进行登录,并查询报告数据。
在进行报告数据查询时,可以根据用户ID确定用户的用户权限信息;然后根据用户权限信息确定与用户ID相关联的一个或多个被测对象ID;查询用户ID相关联的一个或多个被测对象ID的报告数据。
心电工作站根据用户权限信息对报告数据进行数据提取,得到报告数据查询结果列表数据进行输出。报告数据查询结果列表数据至少包括查询得到的一个或多个报告数据的被测对象ID和检测时间信息。
下面以具体的例子来进行说明,此处仅为举例,并不限定其具体实现方式。
比如,用户为某科室主治医生,在其登录心电工作站后,进行患者报告数据查询操作时,心电工作站根据其权限获取该科室的患者的ID,根据患者的ID去查询匹配所有被测对象ID,得到该主治医生有权限查阅的患者的报告数据,优选的可以是设定时间段内的报告数据。然后对报告数据进行数据提取,至少提取其中的被测对象ID和检测时间信息,生成报告数据查询结果列表数据,优选的以列表方式在心电工作站的人机交互界面进行显示。其中,被测对象ID具体可以是姓名或患者档案号等。列表的展示顺序可以按照被测对象ID顺序排列,或者按照检测时间信息顺序/逆序排列。
又如,用户为某患者,在其登录心电工作站后,进行患者报告数据查询操作时,心电工作站根据其权限获取该患者ID相关的所有或设定时间段内的报告数据。同样的,对报告数据进行数据提取,至少提取其中的被测对象ID和检测时间信息,生成报告数据查询结果列表数据,并以列表方式在心电工作站的人机交互界面进行显示。
其中,报告数据查询结果列表数据至少包括查询得到的一个或多个报告数据的被测对象ID、心电监测设备ID和检测时间信息。
步骤140,接收用户输入的选择指令,根据选择指令获取被选择的报告数据;
具体的,用户通过人机交互界面进行选择输入,从而确定索要选择的报告数据。
基于不同的操作目的,支持对报告数据的单个处理操作或批处理操作。
比如需要对报告数据进行输出查看、编辑、修改等操作,可以从报告数据查询结果列表数据选择单个查询结果以获取相应的一个报告数据。
如果是对报告数据进行打印、传送、下载等操作时,可以从报告数据查询结果列表数据同时选择多个查询结果以一次性的获取相应的多个报告数据。这样可以简化操作。
通过心电工作站,可以对用户进行操作权限管理,比如可以设定用户的用户类别信息,例如为医生、普通工作人员、患者等,对于不同类别用户的权限设定为不同。比如用户类别信息为医生,则可以对报告数据进行查询、修改、下载、转发等操作,如果为普通工作人员或患者,只能对报告数据进行查询和打印。
步骤150,接收用户输入的报告数据会诊请求;
具体的,为了最大化的发挥心电工作站的作用,在心电工作站中还包括了收发服务,能够实现报告数据的发送和接收。
心电工作站通过收发服务接入局域网或广域网,心电工作站能够支持以太网、通用串行总线(USB)和RS232串口通信。
心电工作站的人机交互界面提供了会诊服务的服务接口,用户通过选择会诊服务的操作图标,以发起报告数据的会诊请求。
步骤160,根据报告数据会诊请求获取用户ID对应的关联用户的用户ID,并根据关联用户的用户ID将报告数据发送至关联用户的用户设备;
具体的,关联用户的信息可以在心电工作站中进行配置。配置方法可以包括人工配置或自动配置。
人工配置可以是具有相应用户管理权限的用户获取各用户的用户信息,并相应的配置各用户的级别属性信息,然后根据级别属性信息再设定各用户之间的关联信息,并保存为关联用户记录数据。
自动配置的方法可以是可以根据用户ID所对应的用户的级别属性信息确定用户ID对应的关联用户的用户ID,并生成关联用户记录数据。在这里,用户级别属性信息可以以预先设定的规则而定,比如,用户级别属性信息包括用户所在医疗机构的级别信息、用户所在科室信息和用户职务或专业等级的信息,因此根据用户信息中记录的用户的医疗机构的信息、科室信息、职务或专业登记信息等,就能够自动生成用户级别属性信息,并能够根据设定的关联规则,通过用户级别属性信息自动产生关联关系,将相互关联的用户ID进行记录,生成关联用户记录数据。
在心电工作站中对关联用户记录数据进行存储。
为便于理解本方案的意图,下面以一个具体的例子为例对发出会诊请求的过程进行说明。
心电工作站的操作用户为某县级医院的一个医生,该医生操作心电工作站输入报告数据会诊请求。心电工作站接收到报告数据后,根据医生的用户ID查询关联用户信息。在本例中,关联用户优选为该医生所在医院的上级医院的同科室的主治医生或该医生所在医院及上级医院的相关科室的主治医生。关联用户信息即为上述主治医生的信息。
根据关联用户信息,心电工作站进行报告数据的发送,使得关联用户能够接收到该报告数据。
优选的,在输入报告数据会诊请求时,还可以输入会诊请求信息,在进行报告数据发送时,一并发送会诊请求信息。通过会诊请求信息可以描述所需进行会诊的主要项目或问题等等,以使得关联用户能够更清晰的了解会诊请求发起者的需求。
需要说明的是,关联用户用以处理报告数据会诊请求的设备同样的可以是心电工作站。因此,用户也可以通过心电工作站接收来自其他用户的发起的报告数据会诊请求,并输出接收到的报告数据,通过心电工作站的用户界面输入会诊结果反馈数据,并发送给提出报告数据会诊请求的心电工作站。
当然接收报告数据会诊请求的用户设备并不限于心电工作站,还可以是加载运行有相应服务软件的智能手机、平板电脑,也可以是通过相应定制软件或者浏览器获取报告数据的台式机、笔记本电脑等。
步骤170,接收关联用户的用户设备发送的会诊结果反馈数据,并显示。
具体的,在根据报告数据会诊请求将报告数据发出之后,心电工作站可以接收关联用户的用户设备返回的报告数据的会诊结果数据,心电工作站对会诊结果数据进行存储,并自动生成与报告数据的关联信息。
优选的,在接收到会诊结果反馈数据时,心电工作站还可以生成提示信息,以提示用户进行查阅。
当然,在步骤140之后,还可以对报告数据进行本地输出,具体包括如下步骤:
步骤180,接收用户输出的报告输出指令,根据报告输出指令确定输出模式信息和输出格式信息;
具体的,当用户在报告数据查询结果列表数据中选择报告数据并请求对报告数据执行输出时,心电工作站根据用户的选择指令获取被选择的报告数据,并根据报告输出模式信息和报告输出格式信息响应的生成报告输出数据。
同时,心电工作站获取报告输出模式信息和输出格式信息。报告输出模式信息是指报告以哪种格式输出,可以根据用户权限的不同对应有不同的报告输出模式,比如用户类别信息为患者的用户,仅输出报告数据中的设定的部分数据。报告输出模式还包括输出的数据展示方式。在图7a-图7c和图8a-8g中分别展示了不同的报告输出模式,其中图7a-图7c为一种静态心电监测数据的报告输出模式,图8a-8g为一种动态心电检测数据的报告输出模式,用户可以根据需要,输入或选择相应的报告输出模式信息,以使心电工作站输出相应输出模式的报告数据。
报告的输出格式信息可以是根据不同用途而进行选择的输出格式,用以对报告数据进行格式转换,转换为所需要的数据格式再进行输出。具体格式可以包括:可扩展标记语言(XML)、心电图数据记录与传输的标准通信协议(SCP)、便携式文档格式(PDF)、联合图像专家组(JPEG)、医学数字成像和通信(DICOM)3.0、HL7组织制定的心电图数据注释格式(HI7)等多种格式。
心电工作站中预设有默认的报告输出模式信息和输出格式信息,在用户不进行指定时,心电工作站可以按照默认格式进行将报告数据转换为报告输出数据进行输出。
步骤190,根据输出模式信息选择性输出报告结论数据、报告表项数据和/或报告图形数据中的部分数据或全部数据;并且,根据所述输出格式信息,将所述部分数据或全部数据进行数据格式转换,生成报告输出数据。
具体的,报告输出数据的输出包括在心电工作站的人机交互界面上输出,以及在心电工作站中集成的或连接的打印机进行打印输出等。
此外,在步骤170之后,还可以根据会诊结果反馈数据对报告数据进行修改的操作。
具体包括:接收用户根据被选择的报告数据和/或会诊结果反馈数据输入的报告结论修改数据和/或报告表项修改数据;根据报告结论修改数据更新报告结论数据,根据报告表项修改数据更新报告表项数据。
也就是说,在用户查阅了会诊结果反馈数据之后,可以根据会诊反馈的结果进行报告结论数据、报告表项数据的主动修改。
优选的,在报告数据中,对修改时间和执行修改操作的用户ID的信息进行记录。
图6为本发明实施例提供的一种心电工作站系统的结构示意图,该心电工作站系统包括:心电工作站1和多台心电监测设备2。
心电工作站1包括处理器11和存储器12。存储器12可通过总线13与处理器11连接。存储器13可以是非易失存储器,例如硬盘驱动器和闪存,存储器中存储有软件程序和设备驱动程序。软件程序能够执行本发明实施例提供的上述方法的各种功能;设备驱动程序可以是网络和接口驱动程序。处理器11用于执行软件程序,该软件程序被执行时,能够实现本发明实施例提供的方法。
需要说明的是,本发明实施例还提供了一种计算机可读存储介质。该计算机可读存储介质上存储有计算机程序,该计算机程序被处理器执行时,能够实现本发明实施例提供的方法。
本发明实施例还提供了一种包含指令的计算机程序产品。当该计算机程序产品在计算机上运行时,使得处理器执行上述方法。
本发明实施例提供的心电信息处理方法和心电工作站系统,通过对心电监测设备输出的心电图数据进行完整快速的自动分析,形成报告数据;能够对报告数据进行查询,按照请求或设定的输出模式和数据格式输出报告数据,并能设定关联用户,发起对报告数据的会诊请求,使得关联用户能够获取到报告数据并进行会诊处理。
专业人员应该还可以进一步意识到,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、计算机软件或者二者的结合来实现,为了清楚地说明硬件和软件的可互换性,在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本发明的范围。
结合本文中所公开的实施例描述的方法或算法的步骤可以用硬件、处理器执行的软件模块,或者二者的结合来实施。软件模块可以置于随机存储器(RAM)、内存、只读存储器(ROM)、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其它形式的存储介质中。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (11)

1.一种心电信息处理方法,其特征在于,所述方法包括:
心电工作站接收多台心电监测设备输出的心电图数据,所述心电图数据包括被测对象ID、心电监测设备ID和检测时间信息;
对所述心电图数据进行心电数据分析,生成报告数据并进行存储;所述报告数据包括:所述被测对象ID对应的报告结论数据、报告表项数据和/或报告图形数据;
接收用户输入的报告数据查询指令,根据所述用户的用户ID查询相应的报告数据,并生成报告数据查询结果列表数据,进行显示输出;所述报告数据查询结果列表数据至少包括查询得到的一个或多个所述报告数据的被测对象ID、心电监测设备ID和检测时间信息;
接收用户输入的选择指令,根据所述选择指令获取被选择的报告数据;
接收用户输入的报告数据会诊请求;
根据所述报告数据会诊请求获取所述用户ID对应的关联用户的用户ID,并根据所述关联用户的用户ID将所述报告数据发送至所述关联用户的用户设备;
接收所述关联用户的用户设备发送的会诊结果反馈数据,并显示;
其中,所述心电图数据为单导联或多导联的时间序列数据;所述对所述心电图数据进行心电数据分析,生成报告数据具体包括:
将所述心电图数据的数据格式经过重采样转换为预设标准数据格式,并对转换后的预设标准数据格式的心电图数据进行第一滤波处理;
对所述第一滤波处理后的心电图数据进行心搏检测处理,识别所述心电图数据包括的多个心搏数据;每个所述心搏数据对应一个心搏周期,包括相应的P波、QRS波群、T波的幅值数据和起止时间数据;
根据训练得到的干扰识别二分类模型对所述心搏数据进行干扰识别,确定所述心搏数据是否存在干扰噪音,以及用于判断干扰噪音的一个概率值;
根据心搏数据的导联参数和所述心搏数据,基于所述干扰识别的结果和时间规则合并生成心搏时间序列数据;根据所述心搏时间序列数据生成心搏分析数据;
对所述心搏分析数据进行信号质量分析,评估所述心搏分析数据的信号质量,得到所述心搏分析数据的信号质量评价参数;
根据训练得到的心搏分类模型对所述心搏分析数据进行幅值和时间表征数据的特征提取和分析,得到所述心搏分析数据的一次分类信息;
对所述一次分类信息结果中的特定心搏的心搏分析数据输入到训练好的ST段和T波改变模型进行识别,确定ST段和T波评价信息;
根据心搏时间序列数据,对心搏分析数据进行P波和T波特征检测,确定每个心搏中P波和T波的详细特征信息;
对所述心搏分析数据在所述一次分类信息下根据所述心电图基本规律参考数据、P波和T波的详细特征信息以及所述ST段和T波评价信息进行二次分类处理,得到心搏分类信息;
将所述心搏分析数据根据所述心搏分类信息和心电图基本规律参考数据生成心电图事件数据;根据所述信号质量评价参数对心电图事件数据进行筛选,得到所述报告结论数据和报告表项数据;根据每种心电图事件数据中的典型数据片段,生成所述报告图形数据;
所述方法还包括:
接收用户输出的报告输出指令,根据所述报告输出指令确定输出模式信息和输出格式信息;
根据所述输出模式信息选择性输出所述报告结论数据、报告表项数据和/或报告图形数据中的部分数据或全部数据;
根据所述输出格式信息,将所述部分数据或全部数据进行数据格式转换,生成报告输出数据;
所述干扰识别二分类模型为基于人工智能卷积神经网络CNN构建的分类模型;所述干扰识别二分类模型的网络结构包括第一卷积层、第一最大池化层、第二卷积层、第二最大池化层、第一全连接层、第二全连接层和softmax分类器;所述第一卷积层的卷积核大小为1x5,卷积核数目为128;所述第二卷积层的卷积核大小为1x5,卷积核数目为256;所述softmax分类器包括两个分类输出;
所述根据训练得到的心搏分类模型对所述心搏分析数据进行幅值和时间表征数据的特征提取和分析,得到所述心搏分析数据的一次分类信息,具体包括:若所述心搏分析数据为单导联的分析数据,则将所述心搏分析数据输入训练得到的对应导联的心搏分类模型进行幅值和时间表征数据的特征提取和分析,得到所述一次分类信息;若所述心搏分析数据为多导联的分析数据,则对多导联分析方式进行确认;若所述多导联分析方式为投票决策分类方式,则根据训练得到的各导联对应的心搏分类模型对各导联的心搏分析数据进行幅值和时间表征数据的特征提取和分析,得到各导联的分类信息;根据所述各导联的分类信息和导联权重值参考系数进行分类投票决策计算,得到所述一次分类信息;所述导联权重值参考系数为基于心电数据贝叶斯统计分析得到各导联对不同心搏分类的投票权重系数;若所述多导联分析方式为同步关联分类方式,则根据训练得到的多导联同步关联分类模型对各导联的心搏分析数据进行同步幅值和时间表征数据的特征提取和分析,得到所述一次分类信息;
所述报告输出数据的数据格式至少包括:可扩展标记语言XML数据格式、心电图数据记录与传输的标准通信协议SCP数据格式、便携式文档PDF数据格式、联合图像专家组JPEG数据格式、医学数字成像和通信DICOM3.0数据格式和HL7组织制定的心电图数据注释HI7数据格式。
2.根据权利要求1所述的心电信息处理方法,其特征在于,所述根据所述用户的用户ID查询相应的报告数据具体包括:
根据所述用户ID确定所述用户的用户权限信息;
根据所述用户权限信息确定与所述用户ID相关联的一个或多个被测对象ID;
查询所述用户ID相关联的一个或多个被测对象ID的报告数据。
3.根据权利要求1所述的心电信息处理方法,其特征在于,所述方法还包括:
接收所述用户根据所述被选择的报告数据和/或所述会诊结果反馈数据输入的报告结论修改数据和/或报告表项修改数据;
根据所述报告结论修改数据更新所述报告结论数据,根据所述报告表项修改数据更新所述报告表项数据。
4.根据权利要求1所述的心电信息处理方法,其特征在于,在根据所述报告数据会诊请求获取所述用户ID对应的关联用户的用户ID之前,所述方法还包括:
配置各用户的级别属性信息;
根据所述用户ID的级别属性信息确定所述用户ID对应的关联用户的用户ID,并生成关联用户记录数据;
存储所述关联用户记录数据。
5.根据权利要求1所述的心电信息处理方法,其特征在于,对所述第一滤波处理后的心电图数据进行心搏检测处理还包括:
根据所述QRS波群确定RR间期并计算所述RR间期内噪声的估计值;
根据所述噪声的估计值和每个QRS波群中的最大幅值确定每个QRS波群的检测置信度。
6.根据权利要求1所述的心电信息处理方法,其特征在于,所述根据训练得到的干扰识别二分类模型对所述心搏数据进行干扰识别,具体包括:
以第一数据量对所述心搏数据进行切割采样,并将切割采样得到的数据输入到干扰识别二分类模型进行干扰识别;
识别心搏数据中,心搏间期大于等于预设间期判定阈值的数据片段;
对所述心搏间期大于等于预设间期判定阈值的数据片段进行信号异常判断,确定是否为异常信号;
如果不是异常信号,则以预设时间宽度,根据设定时值确定所述数据片段中滑动取样的起始数据点和终止数据点,并由所述起始数据点开始对所述数据片段进行滑动取样,至所述终止数据点为止,得到多个取样数据段;
对每个所述取样数据段进行所述干扰识别。
7.根据权利要求1所述的心电信息处理方法,其特征在于,所述对所述一次分类信息结果中的特定心搏的心搏分析数据输入到训练好的ST段和T波改变模型进行识别,确定ST段和T波评价信息,具体包括:
将所述一次分类信息中的特定心搏数据,按导联依次输入到根据训练得到的ST段和T波改变模型,对各导联的特定心搏数据进行幅值和时间表征数据的特征提取和分析,得到各导联的ST段和T波改变信息,确定ST段和T波评价信息;所述ST段和T波评价信息具体为所述心搏片段数据对应的ST段和T波发生改变的导联位置信息。
8.根据权利要求1所述的心电信息处理方法,其特征在于,所述根据每种心电图事件数据中的典型数据片段,生成所述报告图形数据,具体包括:
根据所述信号质量评价参数,对每种心电图事件的数据片段进行评估,选择信号质量评价参数最高的数据片段为所述心电图事件中的典型数据片段。
9.一种心电工作站系统,其特征在于,所述心电工作站系统包括上述权利要求1-8任一项所述的心电工作站和多台心电监测设备;
其中,所述心电工作站包括存储器和处理器,所述存储器用于存储程序,所述处理器用于执行如权利要求1至8任一项所述的方法。
10.一种包含指令的计算机程序产品,当其在计算机上运行时,使得计算机执行如权利要求1至8任一项所述的方法。
11.一种计算机可读存储介质,包括指令,当所述指令在计算机上运行时,使所述计算机执行根据权利要求1至8任一项所述的方法。
CN201810157158.8A 2018-02-24 2018-02-24 心电信息处理方法和心电工作站系统 Active CN109411041B (zh)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN201810157158.8A CN109411041B (zh) 2018-02-24 2018-02-24 心电信息处理方法和心电工作站系统
US16/766,569 US11350868B2 (en) 2018-02-24 2018-04-18 Electrocardiogram information processing method and electrocardiogram workstation system
EP18907280.4A EP3698709A4 (en) 2018-02-24 2018-04-18 ELECTROCARDIOGRAM INFORMATION PROCESSING AND ELECTROCARDIOGRAM WORKPLACE SYSTEM
PCT/CN2018/083465 WO2019161610A1 (zh) 2018-02-24 2018-04-18 心电信息处理方法和心电工作站系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810157158.8A CN109411041B (zh) 2018-02-24 2018-02-24 心电信息处理方法和心电工作站系统

Publications (2)

Publication Number Publication Date
CN109411041A CN109411041A (zh) 2019-03-01
CN109411041B true CN109411041B (zh) 2021-08-03

Family

ID=65464004

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810157158.8A Active CN109411041B (zh) 2018-02-24 2018-02-24 心电信息处理方法和心电工作站系统

Country Status (4)

Country Link
US (1) US11350868B2 (zh)
EP (1) EP3698709A4 (zh)
CN (1) CN109411041B (zh)
WO (1) WO2019161610A1 (zh)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110164522A (zh) * 2019-04-26 2019-08-23 广州优视一号生物科技有限公司 一种视力检测报告系统
US11751772B2 (en) * 2019-10-31 2023-09-12 Turtle Shell Technologies Private Limited System and a method for myocardial performance determination
CN113208603A (zh) * 2020-01-20 2021-08-06 深圳市理邦精密仪器股份有限公司 心电图机及其报告生成方法以及相关装置
CN113223694A (zh) * 2020-01-20 2021-08-06 深圳市理邦精密仪器股份有限公司 心电图机的心电分析方法、心电图机和存储装置
CN113160923A (zh) * 2020-01-20 2021-07-23 深圳市理邦精密仪器股份有限公司 病患资料的获取方法、心电图机及计算机可读存储介质
CN113712525A (zh) * 2020-05-21 2021-11-30 深圳市理邦精密仪器股份有限公司 一种生理参数处理方法、装置及医疗设备
CN111881158B (zh) * 2020-07-31 2024-06-18 平安国际融资租赁有限公司 一种管理报表数据处理方法、装置、计算机系统及可读存储介质
CN111882005B (zh) * 2020-09-28 2020-12-15 平安科技(深圳)有限公司 数据类别确定方法及装置、终端设备、存储介质
CN112401903B (zh) * 2020-11-03 2023-12-22 沈阳东软智能医疗科技研究院有限公司 心电数据的识别方法、装置、存储介质和电子设备
CN112373215B (zh) * 2020-11-13 2021-10-26 深圳市凯沃尔电子有限公司 一种心电信号的热敏打印方法
CN113010722B (zh) * 2021-04-18 2023-10-03 南通大学 一种融合眼底图像的慢病临床研究队列查询系统及方法
CN113349791A (zh) * 2021-05-31 2021-09-07 平安科技(深圳)有限公司 异常心电信号的检测方法、装置、设备及介质
CN113470804B (zh) * 2021-06-30 2023-07-28 上海联影医疗科技股份有限公司 医学扫描成像系统、权限管理方法、电子设备及存储介质
CN117281529B (zh) * 2023-11-23 2024-01-30 西南医科大学附属医院 一种心电图波形数据智能处理方法

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006296888A (ja) * 2005-04-22 2006-11-02 Fukuda Denshi Co Ltd 生体情報出力装置及び方法並びに生体情報レポート
JP4788927B2 (ja) * 2007-12-26 2011-10-05 日本光電工業株式会社 モニタリングネットワークシステム
KR101002020B1 (ko) * 2008-03-27 2010-12-16 계명대학교 산학협력단 실시간 심전도 모니터링 시스템 및 방법, 패치형 심전도측정장치, 통신장치
CN102028459B (zh) * 2010-12-02 2013-01-02 广东宝莱特医用科技股份有限公司 一种心电图机通道起搏信号检测方法
CN203539346U (zh) * 2013-10-17 2014-04-16 苏州百慧华业精密仪器有限公司 一种支持无线网络及无线打印的心电工作站
CN104042210B (zh) * 2014-07-04 2015-10-28 深圳市凯沃尔电子有限公司 一种数字式心电图机及其系统
US20160045117A1 (en) * 2014-08-14 2016-02-18 Nehemiah T. Liu Peak Detection System and Method for Calculation of Signal-Derived Metrics
JP2016043041A (ja) * 2014-08-22 2016-04-04 セイコーエプソン株式会社 生体情報検出装置及び生体情報検出方法
CN106821366A (zh) 2015-10-14 2017-06-13 张胜国 智能心电监护系统
CN205433651U (zh) 2016-02-23 2016-08-10 济宁中科大象医疗电子科技有限公司 可穿戴式运动及心电信号实时采集及监护设备
CN106027566B (zh) 2016-07-13 2018-10-30 上海数创医疗科技有限公司 基于sip传输协议的移动心电采集及诊断系统
CN107016248A (zh) 2017-04-18 2017-08-04 成都琅瑞医疗技术股份有限公司 一种心电数据分析系统及分析方法
CN107714023B (zh) * 2017-11-27 2020-09-01 上海优加利健康管理有限公司 基于人工智能自学习的静态心电图分析方法和装置

Also Published As

Publication number Publication date
EP3698709A1 (en) 2020-08-26
CN109411041A (zh) 2019-03-01
US20210007618A1 (en) 2021-01-14
EP3698709A4 (en) 2021-08-18
WO2019161610A1 (zh) 2019-08-29
US11350868B2 (en) 2022-06-07

Similar Documents

Publication Publication Date Title
CN109411041B (zh) 心电信息处理方法和心电工作站系统
CN109411042B (zh) 心电信息处理方法和心电工作站
CN108478209B (zh) 心电信息动态监护方法和动态监护系统
CN107837082B (zh) 基于人工智能自学习的心电图自动分析方法和装置
CN107714023B (zh) 基于人工智能自学习的静态心电图分析方法和装置
WO2019100565A1 (zh) 基于人工智能自学习的动态心电图分析方法和装置
CN107358196B (zh) 一种心搏类型的分类方法、装置及心电仪
Mei et al. Automatic atrial fibrillation detection based on heart rate variability and spectral features
CN106815570B (zh) 一种基于动态模式识别的心电信号st-t段识别方法
Sumathi et al. A wavelet transform based feature extraction and classification of cardiac disorder
KR101366101B1 (ko) 개인별 ecg 신호의 정상 신호 분류 시스템 및 방법
US11672464B2 (en) Electrocardiogram processing system for delineation and classification
Rashkovska et al. Clustering of heartbeats from ECG recordings obtained with wireless body sensors
Degirmenci et al. ECG Heartbeat Classification based on Signal-to-Image Transformation Using Convolutional Neural Networks
Singh et al. Signal Quality Evaluation and Processing for QRS Detection in ECG based Smart Healthcare Systems
Banjarnahor et al. Application of Decision Tree Method in ECG Signal Classification For Heart Disorder Detection
CN117390391A (zh) 一种针对动态心电图的数据处理系统
CN117503155A (zh) 一种基于Lorenz散点图的心房颤动诊断系统
Radha et al. Hrudyalysis: A Novel Cloud-Based ECG Analytics Method

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right

Effective date of registration: 20190617

Address after: 102200 Building No. 37 Qianqian Road, Changping Science and Technology Park, Changping District, Beijing

Applicant after: Lepu (Beijing) Medical Equipment Co.,Ltd.

Address before: 518000 North Side of 9 Buildings, 4 Floors, Baiwangxin Industrial Zone, Songbai Road, Xili Street, Nanshan District, Shenzhen City, Guangdong Province

Applicant before: Kai Woer Electronics Co., Ltd. of Shenzhen

TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20200511

Address after: 16th Floor, Block A, 668 Shenzhuan Highway, Songjiang District, Shanghai, 20112

Applicant after: YOCALY INFORMATION SCIENCE & TECHNOLOGY Co.,Ltd.

Address before: 102200, Beijing Changping Changping District science and Technology Park, super Road, No. 3, building 37

Applicant before: Lepu Medical Technology (Beijing) Co.,Ltd.

TA01 Transfer of patent application right
CB02 Change of applicant information

Address after: 201612 16 / F, block a, no.668, Xinzhuan Road, Songjiang District, Shanghai

Applicant after: Shanghai Lepu Yunzhi Technology Co.,Ltd.

Address before: 201612 16 / F, block a, no.668, Xinzhuan Road, Songjiang District, Shanghai

Applicant before: YOCALY INFORMATION SCIENCE & TECHNOLOGY Co.,Ltd.

CB02 Change of applicant information
GR01 Patent grant
GR01 Patent grant