Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP7709232B2 - 心臓疾患予測モデルの生成方法、サーバー及びコンピュータプログラム - Google Patents
[go: Go Back, main page]

JP7709232B2 - 心臓疾患予測モデルの生成方法、サーバー及びコンピュータプログラム - Google Patents

心臓疾患予測モデルの生成方法、サーバー及びコンピュータプログラム

Info

Publication number
JP7709232B2
JP7709232B2 JP2024140097A JP2024140097A JP7709232B2 JP 7709232 B2 JP7709232 B2 JP 7709232B2 JP 2024140097 A JP2024140097 A JP 2024140097A JP 2024140097 A JP2024140097 A JP 2024140097A JP 7709232 B2 JP7709232 B2 JP 7709232B2
Authority
JP
Japan
Prior art keywords
model
data
electrocardiogram data
electrocardiogram
generating
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
JP2024140097A
Other languages
English (en)
Other versions
JP2025031652A (ja
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.)
Synergy AI Co Ltd
Original Assignee
Synergy AI 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
Priority claimed from KR1020230110112A external-priority patent/KR102733909B1/ko
Priority claimed from KR1020240104439A external-priority patent/KR102780189B1/ko
Application filed by Synergy AI Co Ltd filed Critical Synergy AI Co Ltd
Publication of JP2025031652A publication Critical patent/JP2025031652A/ja
Priority to JP2025107995A priority Critical patent/JP7837100B2/ja
Application granted granted Critical
Publication of JP7709232B2 publication Critical patent/JP7709232B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/024Measuring 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 for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/024Measuring pulse rate or heart rate
    • A61B5/0245Measuring 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
    • 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
    • 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/361Detecting fibrillation
    • 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/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/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271Specific aspects of physiological measurement analysis
    • A61B5/7275Determining trends in physiological measurement data; Predicting development of a medical condition based on physiological measurements, e.g. determining a risk factor
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N20/00Machine learning
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N20/00Machine learning
    • G06N20/10Machine learning using kernel methods, e.g. support vector machines [SVM]
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N20/00Machine learning
    • G06N20/20Ensemble learning
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/044Recurrent networks, e.g. Hopfield networks
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • G06N3/0455Auto-encoder networks; Encoder-decoder networks
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/0464Convolutional networks [CNN, ConvNet]
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/047Probabilistic or stochastic networks
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/0475Generative networks
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/048Activation functions
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • G06N3/084Backpropagation, e.g. using gradient descent
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • G06N3/088Non-supervised learning, e.g. competitive learning
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • G06N3/09Supervised learning
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • G06N3/094Adversarial learning
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N5/00Computing arrangements using knowledge-based models
    • G06N5/01Dynamic search techniques; Heuristics; Dynamic trees; Branch-and-bound
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N7/00Computing arrangements based on specific mathematical models
    • G06N7/01Probabilistic graphical models, e.g. probabilistic networks
    • 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/20ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
    • 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/30ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • 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/50ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
    • 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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2560/00Constructional details of operational features of apparatus; Accessories for medical measuring apparatus
    • A61B2560/02Operational features
    • A61B2560/0223Operational features of calibration, e.g. protocols for calibrating sensors

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Biomedical Technology (AREA)
  • General Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Data Mining & Analysis (AREA)
  • Public Health (AREA)
  • Cardiology (AREA)
  • Software Systems (AREA)
  • Mathematical Physics (AREA)
  • Evolutionary Computation (AREA)
  • General Physics & Mathematics (AREA)
  • Pathology (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Veterinary Medicine (AREA)
  • Animal Behavior & Ethology (AREA)
  • Surgery (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Computational Linguistics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Signal Processing (AREA)
  • Psychiatry (AREA)
  • Primary Health Care (AREA)
  • Epidemiology (AREA)
  • Databases & Information Systems (AREA)
  • Fuzzy Systems (AREA)
  • Probability & Statistics with Applications (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)

Description

本発明の様々な実施形態は、将来発生しうる患者の心臓疾患を予測するニューラルネットワークモデルを提供するための方法に係り、より具体的には、ディープラーニングモデルを活用して心電図データに基づいて、心房細動を含む不整脈及びその他の心臓疾患のリスクを評価及び予測する技術に関する。
心電図(Electrocardiogram、ECG)データは、心臓の電気的活動を記録して心臓状態を診断するのに重要なツールとして広く使用されている。既存の多誘導心電図計は、心臓のリズム及び電気的伝導異常を探知して不整脈、特に心房細動などの心臓疾患の診断に核心的な役割を果たしてきた。心房細動を含む不整脈は、不規則又は異常な心臓拍動を特徴とする一般的な疾患であって、無症状又は軽度の症状を示して適切な検診なしでは発見し難い場合が多い。このような不整脈は、治療せずに放置すると、脳卒中、心不全、心臓突然死などの重篤な健康合併症につながる可能性があり、特に、心房細動は、脳卒中及び心不全のリスク増加に関連するため、重要な不整脈タイプの1つである。したがって、心房細動と不整脈の早期発見と介入は、これらの疾患を効果的に管理するために非常に重要である。
一方、既存の心電図計は、主に診断に焦点を置いており、将来の心臓疾患のリスクを予測又は階層化する上で限界がある。これは、既存の心電図分析方法が主に統計的分析に依存するためである。統計的方法論は、非線形的で複雑なパターンを十分に反映せず、患者の個々の心臓状態と様々な不整脈のリスクを精密に評価することが困難である。
さらに、統計ベースのモデルは、大量のデータを処理する際にディープラーニングベースの確率モデルに比べて性能が劣る。心電図データは、筋肉の動き、呼吸、他の装備から発生する電気的ノイズなどの信号の重なりに影響される可能性がある。これらの外部要因は、統計的方法論の精度を低下させるおそれがあり、結果として心臓疾患のリスクを正確に階層化及び定量化するのに大きな障害となる。
韓国公開特許公報第10-2022-0104583号
本発明が解決しようとする課題は、前述した背景技術に対応して案出されたものであり、心電図データに基づいて、将来の心臓疾患のリスクを予測し、階層化することができる技術を提供するためである。
本発明が解決しようとする課題は、上述した課題に限定されず、上述していない別の課題は、以降の記載から通常の技術者に明確に理解できるであろう。
上記の課題を解決するための本発明の一実施形態に係る心臓疾患予測モデルの生成方法が開示される。前記方法は、複数の心電図データを取得するステップと、前記複数の心電図データに基づいて心臓疾患予測のためのモデルを生成するステップと、を含むことができる。
代案的な実施形態において、前記心臓疾患予測のためのモデルを生成するステップは、前記複数の心電図データに基づいて学習データセットを構築するステップと、前記学習データセットに基づいて一つ以上のネットワーク関数に対する学習を行うことにより、心臓疾患リスク予測モデルを生成するステップと、を含み、前記学習データセットは、互いに異なる学習目的に分類された第1学習データセット及び第2学習データセットを含み、前記心臓疾患リスク予測モデルは、患者の心電図データに基づいて、将来の心臓疾患のリスクに関する予測情報を階層化して提供することを特徴とすることができる。
代案的な実施形態において、前記第1学習データセットは、心電図データの特性を特性空間に変換する過程に対応する学習のためのデータを含み、前記第2学習データセットは、心臓リスク予測に関連するキャリブレーションのためのデータを含むことができる。
代案的な実施形態において、前記学習データセットを構築するステップは、前記複数の心電図データに対する前処理を行うステップを含み、前記前処理を行うステップは、前記複数の心電図データに対するノイズ前処理を行うステップと、前記ノイズ前処理された複数の心電図データを所定のウィンドウサイズに分割して複数のROI信号を生成するステップと、前記複数の心電図データのリード単位でHRV特性情報を抽出するステップと、を含み、前記HRV特性情報は、心拍間隔間の変動性に関する指標情報を含むことができる。
代案的な実施形態において、前記心臓疾患リスク予測モデルを生成するステップは、前記第1学習データセットを利用したマスキングベースの自己教師あり学習を介して埋め込みモデルを生成するステップと、前記埋め込みモデルを利用して、前記第2学習データセットに対応する複数のROI信号のそれぞれに対応する潜在ベクトル(latent vector)を抽出する埋め込み実行ステップと、クラスタリングモデルを利用して前記潜在ベクトルに対するクラスタリングを行うステップと、前記潜在ベクトル、及び各潜在ベクトルに対応するそれぞれのクラスターに対する情報に基づいて、ベクトルデータベース(DB、DataBase)を構築するステップと、を含むことができる。
代案的な実施形態において、前記埋め込みモデルを生成するステップは、自己再構成モデルに、前記第1学習データセットに含まれたデータを入力として処理して、前記自己再構成モデルが、前記入力されたデータと類似した出力を生成するように、自己教師あり学習(self-supervised learning)を誘導するステップと、前記学習の完了した自己再構成モデルから、エンコーダを抽出して前記埋め込みモデルを生成するステップと、を含み、前記自己再構成モデルは、入力されるデータの一部をマスキングし、マスキングされた部分を復元するニューラルネットワークモデルであって、次元減少ネットワーク関数及び次元復元ネットワーク関数を含むことができる。
代案的な実施形態において、前記クラスタリングモデルは、前記複数のROI信号のそれぞれに対応する潜在ベクトル間の類似度距離に基づくクラスタリングを行うことを特徴とし、前記クラスタリングを行うステップは、前記クラスタリング結果に対応する複数のクラスターのそれぞれに対する代表ベクトルを抽出するステップを含むことができる。
代案的な実施形態において、前記方法は、予測対象者の心電図データを取得するステップと、前記予測対象者の心電図データに対する前処理を行うステップと、前記埋め込みモデルを利用して、前記前処理の完了した心電図データに対応する複数の潜在ベクトルを生成するステップと、前記複数の潜在ベクトルのそれぞれを前記複数のクラスターのそれぞれに対応する代表ベクトルのそれぞれと比較して、前記複数の潜在ベクトルのそれぞれを前記複数のクラスターのうちの1つに分類するステップと、前記複数の潜在ベクトルのそれぞれが前記複数のクラスターに分類される分類結果に基づいて、心臓疾患のリスクに関する階層化情報を生成するステップと、を含むことができる。
代案的な実施形態において、前記心臓疾患リスク予測モデルを生成するステップは、前記複数の心電図データのそれぞれに対応するユーザーメタ情報を取得するステップと、前記ベクトルデータベースから前記複数の潜在ベクトル、及び前記複数の潜在ベクトルに対応するクラスター情報を取得するステップと、前記HRV特性情報、前記ユーザーメタ情報、前記複数の潜在ベクトル、及び前記潜在ベクトルに対応するクラスター情報に基づいて複数のツリーモデルに対する学習を行い、各ツリーモデルの出力を統合するアンサンブル学習を行って前記心臓疾患リスク予測モデルを生成するステップと、を含むことができる。
代案的な実施形態において、前記心臓疾患予測のためのモデルを生成するステップは、取得された前記複数の心電図データを訓練データセット、検証データセット及び試験データセットに分類するステップと、前記訓練データセット、検証データセット及び試験データセットに分類された心電図データのそれぞれから所定のウィンドウサイズに分割されたROI心電図データを取得するステップと、前記訓練データセットに含まれた心電図データから取得された第1ROI心電図データを患者の心臓疾患を予測するためのディープラーニングモデルに入力することに対応して、前記第1ROI心電図データのそれぞれの心臓疾患クラスを予測するように前記ディープラーニングモデルを学習させるステップと、前記検証データセットに含まれた心電図データから取得された第2ROI心電図データを前記学習されたディープラーニングモデルに適用して前記心臓疾患クラスを区分するための閾値を決定するステップと、をさらに含むことができる。
代案的な実施形態において、前記決定するステップは、前記第2ROI心電図データのそれぞれに対して同一の心電図データから分離されたROI心電図データのそれぞれに対する確率スコアを収集するステップと、前記収集された第2ROI心電図データのそれぞれに対する確率スコアを平均して最終確率スコアを導出するステップと、前記導出された最終確率スコアに基づいて前記心臓疾患クラスを区分するための閾値を微調整するステップと、を含むことができる。
代案的な実施形態において、前記決定された閾値は、前記ディープラーニングモデルの重みと共に保存されることができる。
本発明の他の実施形態によれば、心臓疾患予測モデルの生成方法を行うサーバーが開示される。前記サーバーは、1つ以上のインストラクションを保存するメモリーと、前記メモリーに保存された1つ以上のインストラクションを実行するプロセッサと、を含み、前記プロセッサは、前記1つ以上のインストラクションを実行することにより、前述した心臓疾患予測モデルの生成方法を行うことができる。
本発明の別の実施形態によれば、コンピュータから読み出し可能な記録媒体に保存されたコンピュータプログラムが開示される。前記コンピュータプログラムは、ハードウェアであるコンピュータに結合され、心臓疾患予測モデルの生成方法を行うことができる。
本発明のその他の具体的な事項は、詳細な説明及び図面に含まれている。
本発明の様々な実施形態によれば、患者の心電図データからROI心電図データを抽出し、抽出されたROI心電図データをディープラーニングモデルに適用することにより、患者の心臓疾患予測精度を向上させることができる。
また、本発明は、複数の機械学習モデルを利用して将来の心臓疾患のリスクを予測することにより、患者が本人の健康状態を把握し、心臓疾患を予防することができる効果を提供することができる。これにより、心臓疾患の早期発見と迅速な対応を可能にして、患者の生命を保護し、治療結果を改善することができる機会を提供する。
また、本発明は、心電図データを一々分析しなければならない医療スタッフの物理的時間を自動化されたAIベースの処理過程を介して大幅に減少させることにより、スクリーニングの効率性を高め、医療現場の機会費用を減少させることができる効果を提供することができる。
また、本発明は、HRV特性情報と、埋め込みモデルから抽出された潜在ベクトルと、患者の年齢などのメタデータとをアンサンブルするために、複数のツリー系列機械学習モデルを用いてアンサンブルを行うことにより、心臓疾患リスクの予測精度をさらに向上させることができる。
本発明の効果は、上述した効果に限定されず、上述していない別の効果は、以降の記載から通常の技術者に明確に理解できるであろう。
本発明の一実施形態に係る心臓疾患予測モデルの生成方法を行うコンピューティング装置の構成図である。 本発明の一実施形態に係る心臓疾患を予測するためのディープラーニングモデルの学習方法を示すフローチャートである。 本発明の一実施形態に係る心電図データの分類方法を示す図である。 本発明の一実施形態に係る心電図データを個別心臓拍動に細分化する方法を示す図である。 本発明の一実施形態に係るディープラーニングモデルの種類を示す図である。 本発明の一実施形態に係るディープラーニングモデルの種類を示す図である。 本発明の一実施形態に係るディープラーニングモデルの種類を示す図である。 本発明の一実施形態に係るディープラーニングモデルの学習ステップを示す図である。 本発明の一実施形態に係る学習されたディープラーニングモデルを用いて心臓疾患を予測する方法を示すフローチャートである。 本発明の一実施形態に係るディープラーニングモデルを用いた心臓疾患予測ステップを示す図である。 本発明の一実施形態に係る心臓疾患予測モデルの生成方法に関する例示的なフローチャートである。 本発明の一実施形態に係る複数の心電図データに対する前処理過程を例示的に示すフローチャートである。 本発明の一実施形態に係る潜在ベクトルに対する埋め込みを行い、その遂行結果に基づいてベクトルデータベースを構築する過程を例示的に示すフローチャートである。 本発明の一実施形態に係る埋め込みモデルを生成する過程を例示的に示すフローチャートである。 本発明の一実施形態に係る自己再構成モデル、及び学習された自己再構成モデルから埋め込みモデルを生成する過程を説明するための例示図である。 本発明の一実施形態に係る複数のツリー系列モデルのアンサンブル学習を介して心臓疾患リスク予測モデルを生成する過程を説明するための例示的なフローチャートである。
実施形態に対する特定の構造的又は機能的説明は、例示のための目的のみで開示されたものであって、様々な形態に変更されて実現できる。よって、実際に実現される形態は、開示された特定の実施形態のみに限定されるものではなく、本明細書の範囲は、実施形態で説明した技術的思想に含まれる変更、均等物、又は代替物を含む。
第1又は第2などの用語を、様々な構成要素を説明するために使用できるが、これらの用語は、ある構成要素を他の構成要素から区別する目的でのみ解釈されるべきである。例えば、第1構成要素は第2構成要素と命名でき、同様に、第2構成要素は第1構成要素と命名できる。
ある構成要素が他の構成要素に「連結されて」いると記載された場合には、他の構成要素に直接連結又は接続されていることも、それらの間に別の構成要素が介在することもあると理解されるべきである。
単数の表現は、文脈上明らかに異なる意味を有しない限り、複数の表現を含む。本文書において、「A又はB」、「A及びBのうちの少なくとも1つ」、「A又はBのうちの少なくとも1つ」、「A、B又はC」、「A、B及びCのうちの少なくとも1つ」、及び「A、B又はCのうちの少なくとも1つ」などのフレーズのそれぞれは、そのフレーズのうちの該当するフレーズに一緒に列挙された項目のいずれか、又はそれらの可能なすべての組み合わせを含むことができる。本明細書において、「含む」又は「有する」などの用語は、説明された特徴、数字、ステップ、動作、構成要素、部品、又はこれらの組み合わせが存在することを指定するものであり、1つ又はそれ以上の他の特徴、数字、ステップ、動作、構成要素、部品、又はこれらの組み合わせの存在又は付加の可能性を予め排除しないと理解されるべきである。
他に定義されない限り、技術的又は科学的用語を含めてここで使用されるすべての用語は、当技術分野における通常の知識を有する者によって一般に理解されるのと同じ意味を有する。一般的に使用される辞書で定義されている用語は、関連技術の文脈上の意味と一致する意味を有すると解釈されるべきであり、本明細書で明白に定義しない限り、理想的又は過度に形式的な意味で解釈されない。
以下、実施形態を添付図面を参照して詳細に説明する。添付図面を参照して説明するにあたり、図面符号に関係なく、同一の構成要素は同一の参照符号を付し、これに対する重複説明は省略する。
本発明は、複数の患者に対応する複数の心電図データを取得し、取得された複数の心電図データに基づいて心臓疾患予測のためのモデルを生成することができる。実施形態において、心臓疾患を予測するためのモデルは、心電図データから重要なパターンと特徴を抽出し、これに基づいて心臓疾患の可能性を予測するニューラルネットワークモデルであり得る。本発明のコンピューティング装置は、複数の心電図データを学習に適した形態に加工し、加工されたデータを利用して心臓疾患予測のためのモデルを学習させる。
一実施形態によれば、本発明は、心臓疾患予測のためのモデルとして、ディープラーニングモデル及び心臓疾患リスク予測モデルを生成して提供することができる。実施形態において、本発明は、複数の心電図データを多様な形態に加工し、それぞれ異なる方式でニューラルネットワークを学習させ、様々な予測モデルを生成して提供することができる。以下では、様々な図面を参照して、心臓疾患予測のためのモデルの生成方法を詳細に説明する。
図1は、本発明の一実施形態に係る心臓疾患予測モデルの生成方法を行うコンピューティング装置の構成図である。
図1に示すように、コンピューティング装置100は、1つ以上のプロセッサ110と、プロセッサ110によって行われるプログラム130をロード(load)又は保存するメモリー120と、を含むことができる。図1のコンピューティング装置100に含まれた構成要素は、一例に過ぎず、本発明の属する技術分野における通常の技術者であれば、図1に示された構成要素に加えて他の汎用的な構成要素がさらに含まれることができることが分かる。
プロセッサ110は、コンピューティング装置100の各構成の全般的な動作を制御する。プロセッサ110は、CPU(Central Processing Unit)、MPU(Micro Processor Unit)、MCU(Micro Controller Unit)、GPU(Graphic Processing Unit)、NPU(Neural Processing Unit)、DSP(Digital Signal Processor)、又は本発明の技術分野によく知られている任意の形態のプロセッサのうちの少なくとも1つを含めで構成できる。また、プロセッサ110は、本発明の様々な実施形態に係る方法/動作を実行するための少なくとも1つのアプリケーション又はプログラムに対する演算を行うことができる。コンピューティング装置100は、1つ以上のプロセッサを備えることができる。
メモリー120は、コンピューティング装置100に含まれた構成要素(例えば、プロセッサ110)によって使用される様々なデータ、命令及び情報のうちの1つ又は2つ以上の組み合わせを保存する。メモリー120は、揮発性メモリー及び/又は不揮発性メモリーを含むことができる。
プログラム130は、本発明の様々な実施形態に係る方法/動作が実現された1つ以上の動作(action)を含むことができ、メモリー120にソフトウェア形態で保存されることができる。ここで、動作は、プログラム130で実現される命令語に対応する。例えば、プログラム130は、複数の患者から取得された心電図データを訓練データセット、検証データセット及び試験データセットに分類する動作、訓練データセット、検証データセット及び試験データセットに分類された心電図データのそれぞれから所定のウィンドウサイズに分割されたROI心電図データを取得する動作、訓練データセットに含まれた心電図データから取得された第1ROI心電図データを、患者の心臓疾患を予測するためのディープラーニングモデルに入力することに対応して、第1ROI心電図データのそれぞれの心臓疾患クラスを予測するようにディープラーニングモデルを学習させる動作、検証データセットに含まれた心電図データから取得された第2ROI心電図データを学習されたディープラーニングモデルに適用して心臓疾患クラスを区分するための閾値を決定する動作を行うようにするインストラクションを含むことができる。
実施形態によって、所定のウィンドウサイズに分割されたROI心電図データは、個別心臓拍動を意味することができるが、これに限定されない。以下では、実施形態によって説明の便宜のために所定のウィンドウサイズに分割されたROI心電図データに対して個別心臓拍動を例示として説明することができるが、各実施形態に適用されるROI心電図データの種類は、これに限定されない。
プログラム130がメモリー120にロードされると、プロセッサ110は、プログラム130を実現するための複数の動作を実行させることにより、本発明の様々な実施形態に係る方法/動作を行うことができる。
プログラム130の実行画面は、ディスプレイ140を介して表示できる。図1の場合、ディスプレイ140は、コンピューティング装置100に連結される別途の装置として表されるが、スマートフォン、タブレットなど、ユーザーが携帯しうる端末などのコンピューティング装置100の場合、ディスプレイ140がコンピューティング装置100の構成要素になることができる。ディスプレイ140に表現される画面は、プログラムに情報を入力する前又はプログラムの実行結果であり得る。
図2は、本発明の一実施形態に係る心臓疾患を予測するためのディープラーニングモデルの学習方法を示すフローチャートである。
図2に示されたディープラーニングモデルの学習方法は、図1に示されたコンピューティング装置のプロセッサによって行われる。ステップ210で、プロセッサは、複数の患者から取得された心電図データを訓練データセット、検証データセット及び試験データセットに分類することができる。まず、複数の患者から取得された心電図データは、正常洞リズム(True-Normal Sinus Rhythm、T-NSR)、心房細動-正常洞リズム(Afib-Normal Sinus Rhythm、AF-NSR)及び臨床的に重要な不整脈-正常洞リズム(Clinically Important Arrhythmia-Normal Sinus Rhythm、CIA-NSR)のグループに区分されることができる。このとき、複数の患者から取得された心電図データは、10秒12リード(lead)心電図データであり得るが、このような心電図データの種類は、一つの例示に過ぎず、上記の例に限定されない。
一方、臨床的に重要な不整脈(Clinically important arrhythmia、CIA)は、Atrial arrhythmia(心房性不整脈)、Ventricular arrhythmia(心室性不整脈)、Atrial fibrillation(心房細動)及びBBB(脚ブロック)を含むことができる。Atrial arrhythmiaは、心房に起因する不整脈であり、30秒以上の持続的又は30秒以内の非持続的な不整脈を意味し、Atrial premature complex、or Sustained/non-sustained Atrial rhythmなどが含まれることができる。Ventricular arrhythmiaは、心室に起因する不整脈であって、30秒以上の持続的又は30秒以内の非持続的な不整脈を意味し、Ventricular premature complex、or Sustained/Non-sustained Ventricular arrhythmiaなどが含まれることができる。Atrial fibrillationは、心房の規則的な電気信号と収縮がなく、これによる不規則な心室の収縮を示す不整脈を意味することができる。最後に、BBBは、心室を介した心臓の信号を伝達する右脚(Right bundle)又は左脚(Left bundel)の信号伝達が遮断されて特徴的な心電図上のパターンを示す不整脈を意味することができる。
T-NSRグループは、心房細動や不整脈の病歴がなく、1年に正常洞リズム心電図が3回以上ある患者の心電図で構成できる。AF-NSRグループは、正常洞リズム心電図と、当該正常洞リズム心電図から14日以内に発生した心房細動又は心房粗動心電図とが対をなす患者の心電図で構成できる。同様に、CIA-NSRグループは、正常洞リズム心電図と、当該正常洞リズム心電図から14日以内に発生した不整脈心電図とが対をなす患者の心電図で構成できる。
プロセッサは、このように区分されたT-NSRグループ、AF-NSRグループ及びCIA-NSRグループに属する心電図データを、心臓疾患予測のためのディープラーニングモデルに適用可能な訓練データセット、検証データセット及び試験データセットに分類することができる。
一例として、図3は、本発明の一実施形態に係る心電図データの分類方法を示す図である。プロセッサは、T-NSRグループ、AF-NSRグループ、及びCIA-NSRグループに属する心電図データ310を予測しようとする心臓疾患別に区分することができる。一例として、プロセッサは、心電図データ310を、心房細動を予測するためにT-NSRグループ及びAF-NSRグループに属する心電図データ320と、不整脈を予測するためにT-NSRグループ及びCIA-NSRグループに属する心電図データ330とに区分することができる。
その後、プロセッサは、それぞれの心臓疾患別に区分された心電図データ320、330を、当該心電図データ320、330が生成された任意の日付を基準に訓練データセット、検証データセット、及び試験データセットに分類することができる。
一例として、図3に示すように、心電図データ320、330が2017年5月23日~2022年5月23日の間に取得された場合、プロセッサは、心電図データ320、330を任意の日付(例えば、2021年6月11日)に基づいて、該当日より前の心電図データを一定の割合によって訓練データセット及び検証データセットに分類し、該当日以後(該当日を含む)の心電図データを試験データセットに分類することができる。図3の例では、心電図データ320、330が訓練データセット(60%)、検証データセット(20%)及び試験データセット(20%)の割合に分類されたが、このような分類割合は、1つの例示に過ぎず、上記の例に限定されない。
ステップ220で、プロセッサは、訓練データセット、検証データセット及び試験データセットに分類された心電図データのそれぞれから個別心臓拍動を取得することができる。プロセッサは、正確で信頼できるデータを取得するために、訓練データセット、検証データセット及び試験データセットに分類された10秒12リード心電図データに対して前処理を行うことができる。
より具体的には、プロセッサは、訓練データセット、検証データセット及び試験データセットに分類された心電図データの入力をXML(eXtensible Markup Language)ファイル形式で受けることができる。プロセッサは、入力データから、患者の名前、年齢、性別などの定型データ部分と連続的な信号からなる非定型データ部分をパーシング(parsing)することができる。その後、プロセッサは、連続的な信号からなる非定型データ部分に対してノイズ除去前処理を行い、心臓拍動のキーマーカーを区分することができる。プロセッサは、このように区分された心臓拍動のキーマーカーを介して、連続的な信号からなる非定型データから不連続的な個別心臓拍動を取得することができる。
一例として、図4を参照すると、プロセッサは、10秒12リード心電図データをBase64暗号化でデコードした後、ノイズ除去及びクレンジングのために移動平均カーネルのあるIIRバタワースSOSフィルターと電力線ノイズフィルターを通過させることができる。次に、プロセッサは、ノイズ除去された10秒12リード心電図データを、QRSピーク感知アルゴリズムを用いて個別心臓拍動に細分化することができる。
その結果、プロセッサは、図3に示すように、単一の10秒12リード心電図データから複数の不連続的な個別心臓拍動を取得することができ、このように取得された個別心臓拍動は、より正確な心臓疾患予測のためにディープラーニングモデルを学習するのに使用できる。
ステップ230で、プロセッサは、訓練データセットに含まれた心電図データから取得された第1個別心臓拍動を、患者の心臓疾患を予測するためのディープラーニングモデルに入力することに対応して、第1個別心臓拍動のそれぞれの心臓疾患クラスを予測するようにディープラーニングモデルを学習させることができる。
このとき、プロセッサは、ResNet-18、LSTM(Long Short-Term Memory)が含まれたConv1D、及びトランスフォーマが含まれたConv1Dのうちのいずれかのディープラーニングモデルを用いて、第1個別心臓拍動のそれぞれの心臓疾患クラスを予測するように学習させることができる。
一例として、図5aのようなResNet-18は、様々なコンボリューションニューラルネットワーク(Convolutional Neural Networks、CNN)などの畳み込み演算を用いて入力の必須特徴を抽出することができるディープラーニングモデルである。CNNアーキテクチャの勾配消失(Vanishing Gradient)問題を解決するために、ResNet-18は、図5aに示すように、入力データがネットワーク上の複数の層をスキップして、出力層に直接連結される方式であるスキップ連結を介して残余学習を行うことができる。
ResNet-18のディープラーニングモデルは、固定長の入力を必要とするため、プロセッサは、個別心臓拍動の長さをすべての個別心臓拍動の平均長さに固定することができる。一例として、全ての個別心臓拍動の平均長さが700である場合、プロセッサは、700より長い心臓拍動長さを有する個別心臓拍動に対してはスライシング(Slicing)を行い、700より短い心臓拍動長さを有する個別心臓拍動には、ゼロパディングを行うことで、長さが700となるように固定することができる。
他の一例として、図5bのようなLSTMが含まれたConv1Dは、順次データからローカル時間パターン及び長距離時間パターンを全てキャプチャすることができる。このとき、Conv1D層は、ローカル時間パターンの感知に優れており、LSTM層は、長期的な従属性のモデル化に優れている。
別の一例として、図5cのようなトランスフォーマが含まれたConv1Dは、入力データのローカルパターン及びグローバル従属性を全てキャプチャすることができる。このとき、トランスフォーマ層は、グローバル従属性のモデル化に適しており、Conv1D層は、ローカルパターンの感知に効果的であり得る。入力長が固定されているResNet-18とは異なり、トランスフォーマが含まれたConv1Dは、様々な入力サイズを収容することができるという利点がある。
一方、図6は、本発明の一実施形態に係るディープラーニングモデルの学習ステップを示す図である。図6を参照すると、プロセッサは、対数損失のある2進交差エントロピーと初期学習率0.0001のAdamWオプティマイザーを用いて、ディープラーニングモデルのパラメータを最適化することができる。このとき、2進交差エントロピーは、ディープラーニングモデルの学習における予測結果と実際の正解との差を減らす損失関数であり、AdamWオプティマイザーは、このような損失関数に基づいて実際のディープラーニングモデルの更新に関与するアルゴリズムである。プロセッサは、このような2進交差エントロピー及びAdamWオプティマイザーに最適化されたディープラーニングモデルの出力に対してシグモイド(Sigmoid)関数を適用することにより、第1個別心臓拍動のそれぞれの心臓疾患クラスに対する確率値を0~1の範囲で取得することができる。
ステップ240で、プロセッサは、検証データセットに含まれた心電図データから取得された第2個別心臓拍動率を、学習されたディープラーニングモデルに適用して、心臓疾患クラスを区分するための閾値を決定することができる。
より具体的には、プロセッサは、図6に示すように、第2個別心臓拍動のそれぞれに対して同一の心電図データから分離された個別心臓拍動のそれぞれに対する確率値を収集し、収集されたすべての第2個別心臓拍動のそれぞれに対する確率値を収集し、収集された全ての第2個別心臓拍動に対する確率値を平均して最終確率スコアを導出した後、導出された最終確率スコアに基づいて心臓疾患クラスを区分するための最適な閾値を微調整することができる。
このとき、最適な閾値は、0~1の間の閾値を0.01単位で適用して心臓疾患クラス、すなわちT-NSR及びAF-NSRクラス又はT-NSR及びCIA-NSRクラスに対する検証データセットにおける最高のF1スコアを達成することにより取得できる。
このとき、F1スコアは、下記式1のように定義できる。
ここで、
は、正しく分類された予測陽性の割合を示し、検出された結果がどれほど正確であるか、すなわち、検出結果のうち実際の物体がどれほど含まれているかを示すことができる。
は、正しく分類された真陽性の割合を示し、実際の対象物体を見逃すことなくどれほどよく捕捉するかを示す指標である。また、
は、それぞれ実陽性、真陰性に分類されたサンプルを意味し、
は、それぞれ偽陽性、偽陰性に分類されたサンプルを意味する。
プロセッサは、このような心臓疾患クラスに対する最適な閾値をディープラーニングモデルの重みと共に保存することができる。
図7は、本発明の一実施形態に係る学習されたディープラーニングモデルを用いて心臓疾患を予測する方法を示すフローチャートである。
図7に示されたディープラーニングモデルを用いた心臓疾患予測方法は、図1に示されたコンピューティング装置のプロセッサによって行われる。プロセッサは、ステップ710のように心臓疾患クラスを予測しようとする患者の心電図データを受信することで、ステップ720のように受信された心電図データから個別心臓拍動を取得することができる。このとき、受信される心電図データ及び取得される個別心臓拍動は、図3に開示された試験データセットに対応することができる。
ステップ730で、プロセッサは、取得された個別心臓拍動を、患者の心臓疾患予測のためのディープラーニングモデルに入力して、個別心臓拍動のそれぞれに対して特定の心臓疾患クラスとして予測される確率値を導出することができる。
ステップ740で、プロセッサは、個別心臓拍動のそれぞれに対して導出された確率値を用いて患者の心臓疾患を予測することができる。より具体的には、プロセッサは、すべての個別心臓拍動のそれぞれに対して導出された確率値の平均を閾値と比較することにより、患者の心臓疾患クラスを区分することができる。
一例として、図8は、本発明の一実施形態に係るディープラーニングモデルを用いた心臓疾患予測ステップを示す図である。図8を参照すると、プロセッサは、学習されたディープラーニングモデルの重み及び閾値をロード(Load)した後、個別心臓拍動のそれぞれに対する確率値を導出することができる。プロセッサは、このように導出されたすべての個別心臓拍動の確率値の平均値を求めることができる。
一例として、プロセッサは、T-NSRと判断する確率値の平均であるT-NSR Logit値が第1閾値θよりも大きく、且つCIA-NSRと判断する確率値の平均であるCA-NSR Logit値が第2閾値θより小さい場合、当該患者の心臓疾患クラスを正常状態、すなわちT-NSRと判断することができる。
又は、プロセッサは、T-NSR Logit値が第1閾値θより大きく、且つCIA-NSR Logit値が第2閾値θより大きい場合、T-NSR Logit値とCIA-NSR Logit値とを比較してより大きい値を最終予測として選択することができる。一例として、CIA-NSR Logit値がT-NSR Logit値より大きい場合、当該患者の心臓疾患クラスを、不整脈が発生しうるCIA-NSRと判断することができる。
また、プロセッサは、T-NSR Logit値が第1閾値θより小さく、且つCIA-NSR Logit値が第2閾値θより小さい場合、T-NSR Logit値とCIA-NSR Logit値とを比較してより大きい値を最終予測として選択することができる。
最後に、プロセッサは、T-NSR Logit値が第1閾値θより小さく、且つCIA-NSR Logit値が第2閾値θより大きい場合、該当患者の心臓疾患クラスを不整脈が発生しうるCIA-NSRと判断することができる。
このような心臓疾患クラスを区分する4つの方法は、T-NSRとAF-NSRの区分にも同様に適用できる。
以下では、図9~図14を参照して、本発明の一実施形態に係るHRV特性情報、潜在ベクトル、潜在ベクトルに対応するクラスター情報、及びメタデータを利用する方法と、複数のツリーモデル出力をアンサンブルして心臓疾患を予測する方法について具体的に説明する。
図9は、本発明の一実施形態に係る心臓疾患予測モデルの生成方法に関する例示的なフローチャートを示す。
本発明の一実施形態によれば、心臓疾患予測モデルの生成方法は、複数の心電図データを取得するステップ(S1000)を含むことができる。
実施形態によれば、複数の心電図データは、複数の患者に対応して取得されるものであり、様々な方式で収集できる。一実施形態において、本発明における心電図データは、心臓の電気的活動を記録して心臓状態を評価するのに関連する多誘導心電図(ECG、Electrocardiogram)データを意味することができる。多誘導心電図データは、様々な誘導(すなわち、リード)を介して得られた心臓電気信号を含み、これにより心臓のリズム、電気伝導状態及びその他の心臓機能に対する総合的な情報を提供する。
実施形態において、複数の患者から取得される複数の心電図データは、正常心電図データ、不整脈に関連する心電図データ、及びその他の心臓状態に対する心電図データを含むことができる。具体的な実施形態において、複数の患者から取得された心電図データは、正常洞リズム(True-Normal Sinus Rhythm、T-NSR)、心房細動-正常洞リズム(Afib-Normal Sinus Rhythm、AF-NSR)、及び臨床的に重要な不整脈-正常洞リズム(Clinically Important Arrhythmia-Normal Sinus Rhythm、CIA-NSR)のグループに区分されることができる。このとき、複数の患者から取得された心電図データは、10秒12リード心電図データであり得るが、このような心電図データの種類は、一例に限定されない。
一実施形態において、複数の心電図データの取得は、メモリー120に保存されたデータを受信又はロード(load)することであり得る。複数の心電図データの取得は、有線/無線通信手段に基づいて他のストレージメディアに、他のコンピューティング装置、同一のコンピューティング装置内の別途の処理モジュールから複数の心電図データを受信又はロードすることであり得る。一例として、ユーザー(例えば、患者、又は医療提供者)は、ユーザー端末を介してコンピューティング装置100に接続してコンピューティング装置100から心臓疾患予測のためのユーザーインタフェースの提供を受け、提供されたユーザーインタフェースに心電図データをドラッグアンドドロップ(drag and drop)する方式でコンピューティング装置100へ心電図データを伝達することができる。
様々な実施形態によれば、本発明のコンピューティング装置100は、取得される複数の心電図データに対するデータ増強を行うことができる。具体的には、コンピューティング装置100は、リード対を利用したデータ増強を行うことにより、ニューラルネットワークの学習のための学習データを増強させることができる。
詳細に説明すると、心電図データのそれぞれの場合、リード別に取得できる。例えば、1-lead、3-lead、12-leadなどの様々なリード構成で心電図データが収集されることができ、これは、心臓の様々な部分における電気的活動を記録してより総合的な心臓状態を評価するのに使用される。
例えば、12-leadシグナルの場合、比較的短い時間の間、多くの信号を取得することができるが、総データ量が制限的であり得る。一方、1-lead(ワンリード)シグナルは、さらに長い時間の間記録されるが、データの多様性が制限されている。これらの2類型のデータが一緒に使用されると、より正確な心臓状態評価が可能であるが、実際に、このようなデータペアの数は不十分であることがある。
この問題を解決するために、本発明のコンピューティング装置100は、1つのリードデータに基づいて他のリードデータを生成する再構築モデルを生成することができる。再構築モデルは、ワンリードデータを入力として受け、これに基づいてマルチリードデータを生成する役割を果たす。この過程は、ワンリードデータを3-lead、又はそれ以上のリード(12-lead)に拡張して、より多くのリードデータがあるかのように学習データを増強する効果を持つ。
再構築モデルの学習は、ワンリードとマルチリードデータとの相関関係を学習する過程を含む。このために、ワンリードとマルチリードデータの両方が存在するペアデータを用いて再構築モデルを学習させる。学習過程で、モデルは、ワンリードデータからマルチリードデータを再現する方法を学び、このとき、ワンリードデータの特徴とこれから類推することが可能なマルチリード信号のパターンを学習する。
モデル学習が完了した後、再構築モデルは、ワンリードデータのみでもマルチリード信号を生成することができ、これにより実際のマルチリードデータの不足を補うことができる。これにより、データの多様性が増加し、モデルの一般化性能が向上する。このように生成されたマルチリードデータ(すなわち、補強されたマルチリードデータ)は、既存のマルチリードデータと共にニューラルネットワーク学習に使用されることにより、モデルがさらに様々なリード構成を学習し、様々な心臓状態に対する予測性能を向上させることができる。
結果的に、再構築モデルによるデータ増強は、学習データの量を実質的に増加させて、モデルの予測精度を高め、心臓疾患予測モデルの信頼性を強化させる。これにより、医療スタッフが患者の健康状態をさらに正確に診断し、より効果的な治療及び予防計画を樹立することができるように役立つ。
本発明の一実施形態によれば、心臓疾患予測モデルの生成方法は、複数の心電図データに基づいて学習データセットを構築するステップ(S2000)を含むことができる。
一実施形態において、複数の心電図データを介して学習データセットを構築することは、各心電図データの重要な特徴とパターンを学習するために、様々なデータを収集及び構成する過程である。これにより、心臓疾患の予測に必要な様々な事例とパターンを含むデータセットが生成される。
実施形態によれば、本発明の学習データセットは、心臓疾患の予測に関連する様々なデータを含むものであり、ニューラルネットワークモデルを訓練するための訓練データセット、モデルの性能を評価し最適化するための検証データセット、及びモデルの一般化性能を評価するための試験データを含むことができる。
一実施形態において、学習データセットは、互いに異なる学習目的に分類された第1学習データセット及び第2学習データを含むことができる。第1学習データセット及び第2学習データセットは、互いに異なるニューラルネットワークモデルを学習させるために分類されたデータであり得る。第1学習データセット及び第2学習データセットは、コンピューティング装置100によって分類できる。
例えば、最初の学習データセットは、埋め込みモデルの学習過程で利用され、当該データセットは、心電図データの主要特徴を抽出し表現するのに重点を置く。第2学習データセットは、心臓疾患のリスクを予測するモデルの学習に使用され、当該データセットは、予測モデルの精度を高めるためのキャリブレーション及び性能評価に使用され得る。
具体的には、第1学習データセットは、心電図データの特性を特性空間に変換する過程に対応する学習のためのデータを含むことができる。第1学習データセットは、ECG表現学習(ECG Representation Learning)のためのデータであり、心電図データ信号の多様な特性とパターンを抽出し、埋め込むのに利用できる。また、第2学習データセットは、心臓のリスクの予測に関連するキャリブレーションのためのデータを含むことができる。すなわち、第2学習データセットは、心臓のリスクを補正するためのデータであり、モデルが予測したリスクをより正確に当てるために利用できる。
実施形態によれば、学習データセットを構築するステップは、複数の心電図データに対する前処理を行うステップを含むことができる。
図10は、本発明の一実施形態に係る複数の心電図データに対する前処理過程を例示的に示すフローチャートである。
図10を参照すると、複数の心電図データに対する前処理方法は、複数の心電図データに対するノイズ前処理を行うステップ(S2100)、ノイズ前処理された複数の心電図データを所定のウィンドウサイズに分割して複数のROI信号を取得するステップ(S2200)と、複数の心電図データのリード単位でHRV(Heart Rate Variability)特性情報を抽出するステップ(S2300)と、を含むことができる。一実施形態において、HRV特性情報は、心拍間隔間の変動性に関する指標情報を含むことができる。
より詳細に説明すると、まず、ステップS2100で取得された複数の心電図データからノイズ除去のための前処理が行われる。ここで、ノイズ除去のための前処理は、ローパスフィルタリング(Lowpass Filtering)及びノイズ除去(Noise Removal)を含み、これは、心電図信号のベースラインノイズを除去するためのものである。具体的には、ローパスフィルタリングは、心電図データから周波数帯域制限方式を介して高周波成分を除去して信号の純度を高め、ノイズ除去は、時系列分析技法を活用して心電図データから不要な雑音を減らして重要な信号成分をより明確にする。このような前処理過程は、心電図信号の歪みを最小限に抑え、心臓活動をより正確に評価するための信号品質の向上に寄与する。
また、実施形態によれば、コンピューティング装置100は、モデル学習の多様性とロバスト性を高めるためのデータ増強を行うことができる。コンピューティング装置100は、様々なノイズパターンをシミュレートする方式によってデータ増強を行い、これによりモデルが実際の環境で発生しうる様々なノイズ状況に対して学習し適応することができるようにする。
具体的な実施形態において、心電図データの前処理過程で、ガウスノイズ(Gaussian noise)、ジッターノイズ(Jitter noise)、ランダムタイムウォーピング(Random Time Warping)、スケーリング(Scaling)、アンプリチュードモジュレーション(Amplitude Modulation)、及びフェーズシフティング(Phase Shifting)などの実際の心電図測定環境で発生しうる様々なノイズを人為的に追加してデータセットを増加することにより、モデルが様々な環境で発生しうるノイズに対して学習し適応することができるようにする。このようなノイズの追加及び除去過程は、データの品質を向上させ、モデルの一般化性能を向上させるという利点がある。
その後、ステップS2200でノイズ除去された心電図データは、所定のウィンドウサイズに分割されて複数のROI信号に整理される。実施形態によって、ROI信号は、心電図データのそれぞれから意味のある心臓拍動が含まれた変形可能なウィンドウサイズ単位で、QRS Complex(心臓拍動)を少なくとも3~5個含む範囲として定義される。一実施形態において、QRS Complexは、心臓の電気的活動中の心室の脱分極を示し、心電図における主要な特徴として心拍動の開始と終了を識別することに寄与することができる。例えば、各ROI信号は、30秒~1分単位に分割でき、これは、心電図信号における重要なパターンを分析し、心臓拍動の規則性及び異常を識別するのに有用である。
コンピューティング装置100は、ROI信号を抽出して、心臓拍動の主要特性をより正確に分析することができるデータに整理する。これは、心電図信号の特定のパターンと異常の兆候を識別し、心臓疾患のリスクを評価するのに重要な役割を果たす。
また、ステップ(S2300)で、HRV特性情報抽出過程が行われる。HRV特性情報は、時間の経過による心拍の周期的な変化を測定する指標であり、心臓の活動を反映する重要な信号特性である。例えば、コンピューティング装置100は、心電図データを30秒~1分単位に分割して複数のROI信号を取得し、各ROI信号に対応してHRV特性を抽出する。これは、既存のビット単位分析よりもさらに長い時間範囲でデータを分析することにより、心臓の全般的な活動性をより正確に評価するためである。
コンピューティング装置100は、周波数分析又は時間領域分析などの方式によって各ROI信号からHRV特性を導出する。具体的な実施形態において、コンピューティング装置100は、様々な分析技法を適用して各ROI信号に対応するHRV特性情報を抽出することができる。具体的な例として、コンピューティング装置は、スペクトル分析による高周波及び低周波成分の評価、及びNN間隔の統計的分析を介してHRV特性情報を抽出することができる。
スペクトル分析は、周波数領域で心拍変異度を評価して、高周波成分は主に副交感神経系の活動を反映し、低周波成分は交感及び副交感神経系の混合活動を示す。これは、自律神経系のバランス状態を把握するのに有用である。NN間隔の統計的分析は、時間領域における心拍間隔の変動性を評価して心臓の自律神経系活動性を定量化する。これは、心拍の不規則性、ストレス反応、心臓の健康状態などを評価する上で重要な指標になる可能性がある。
一実施形態において、リード(lead)単位で抽出されたHRV特性情報には、RMSSD(Root Mean Square of Successive Differences)、SDNN(Standard Deviation of NN intervals)、SDANN(Standard Deviation of Average NN intervals)などが含まれることができる。具体的に、RMSSDは、連続的な心拍間隔の差の二乗平均の平方根を意味し、主に副交感神経系の影響を反映する。SDNNは、すべての正常心拍間隔の標準偏差を示し、自律神経系の全体的な変動性を示す。SDANNは、一定期間の平均心拍間隔の標準偏差を測定して長期的な心拍変動性を評価する。これらの指標は、心臓の自律神経系の活動性を評価し、心臓の健康状態を総合的に分析する上で重要な役割を果たす。すなわち、HRV特性情報は、それぞれ心拍変異度の様々な側面を反映し、心拍動の不規則性、自律神経系活動性などを評価するのに使用され、心臓疾患のリスクの予測に基盤となり得る。
本発明の一実施形態によれば、心臓疾患予測モデルの生成方法は、学習データセットに基づいて1つ以上のネットワーク関数に対する学習を行うことにより、心臓疾患リスク予測モデルを生成するステップ(S3000)を含むことができる。
図11は、本発明の一実施形態に係る潜在ベクトルに対する埋め込みを行い、その遂行結果に基づいてベクトルデータベースを構築する過程を例示的に示すフローチャートである。
図11を参照すると、心臓疾患リスク予測モデルを生成するステップは、第1学習データセットを利用したマスキングベースの自己教師あり学習を介して埋め込みモデルを生成するステップ(S3100)と、埋め込みモデルを利用して第2学習データセットに対応する複数のROI信号のそれぞれに対応する潜在ベクトル(latent vector)を抽出する埋め込み実行ステップ(S3200)と、クラスタリングモデルを利用して潜在ベクトルに対するクラスタリングを行うステップ(S3300)と、潜在ベクトル、及び各潜在ベクトルに対応するそれぞれのクラスターに対する情報に基づいてベクトルデータベース(DB、DataBase)を構築するステップ(S3400)と、を含むことができる。
より詳細に説明すると、ステップS3100で、コンピューティング装置100は、第1学習データセットに基づいて埋め込みモデルを生成することができる。
埋め込みモデルは、心電図データの複雑なパターンを低次元空間で表現することができるように学習されたモデルであり得る。埋め込みモデルは、心電図データに対応する複数のROI信号が入力される場合、各ROI信号に対応する潜在ベクトルを生成することができる。
ここで、潜在ベクトルは、心電図データの重要な特徴を圧縮的に表現した低次元ベクトルであってもよい。これは、高次元データである心電図信号を理解し分析しやすくするために、信号の主要情報を含んでいる形態に変換されたデータ表現である。潜在ベクトルは、心電図データのパターンに基づいて心臓の健康状態を評価するか、或いは疾患の兆候を感知するのに有用な情報を提供することができる。
埋め込みモデルを生成する方法についての具体的な説明は、以下で図12及び図13を参照して後述する。
図12は、本発明の一実施形態に係る埋め込みモデルを生成する過程を例示的に示すフローチャートである。図13は、本発明の一実施形態に係る自己再構成モデル及び学習された自己再構成モデルから埋め込みモデルを生成する過程を説明するための例示図である。
図12を参照すると、埋め込みモデルを生成するステップは、自己再構成モデルに、第1学習データセットに含まれたデータを入力として処理して、自己再構成モデルが、前記入力されたデータと同様の出力を生成するように、自己教師あり学習(self-supervised learning)を誘導するステップ(S3110)と、学習の完了した自己再構成モデルから、エンコーダを抽出して埋め込みモデルを生成するステップ(S3120)と、を含むことができる。
自己再構成モデルは、入力と出力を介して学習する過程でデータの重要な特徴を内在的に学習するように設計されたモデルであり得る。
一実施形態において、自己再構成モデルは、入力されるデータの一部をマスキングし、マスキングされた部分を復元するニューラルネットワークモデルであり、エンコーダ(又は、次元減少ネットワーク関数)及びデコーダ(又は、次元復元ネットワーク関数)を含むことができる。一実施形態において、本発明の自己再構成モデルは、心電図データで重要なパターンを学習して、マスキングされた信号を正確に復元するMAE(Masked AutoEncoder)であり得るが、これに限定されない。
自己再構成モデルは、心電図信号(具体的には、各ROI信号)の一部を意図的に隠し、残りの部分に基づいて信号全体を再構成する過程を介してデータの重要な特徴を効果的に学習することができる。
自己再構成モデルは、入力と出力との間の復元誤りを最小限に抑える方式で学習を行い、これにより信号の重要な特性を効果的に抽出し表現することができる。自己再構成モデルは、特にデータが複雑又は不完全であるとき、データの内在された構造を理解し活用することに強みがある。
より詳細に説明すると、図13に示すように、自己再構成モデル800は、入力データを次元縮小させた後、再び復元する過程によって重要な特徴を学習することができる。自己再構成モデルは、エンコーダ(又は次元減少ネットワーク関数)810及びデコーダ(又は次元復元ネットワーク関数)820を含むことができる。自己再構成モデル800は、エンコーダ810を介して高次元データを低次元潜在空間に変換し、この過程でデータの核心情報を圧縮し、ノイズを低減する。その後、デコーダ820を用いて、圧縮された低次元表現を元の高次元空間に復元しようとし、これによりデータの重要なパターンを学習し再現することができる。
実施形態によれば、コンピューティング装置100は、自己再構成モデル800を、マスキング学習方式を介して学習されるように誘導することができる。この方式では、入力データの一部を意図的にマスキングし、モデルがマスキングされた部分を復元するように学習する。例えば、マスクを被せなければ、自己再構成モデルは、QRS Complexではなく、ノイズに偏って学習できるため、マスキングを介して重要な特徴に集中して学習することができるようにする。これにより、モデルは、入力データの残りの部分から重要なパターンと特徴を抽出して、マスキングされた情報を正確に復元する能力を開発する。これは、データのノイズや不完全な部分にも頑強に作動することができるという利点がある。
また、マスキングされた部分を復元する過程でデータの脈絡を理解し、重要な特徴を学習することができるため、心電図データなどの複雑な医療データを分析し解釈するのに非常に有用である。マスクを適用することにより、モデルは、QRS Complexなどの主要信号に集中して、心臓疾患の正確な診断と評価のための高信頼性のデータ分析を可能にする。
すなわち、本発明は、心電図データの分析と解釈において信頼性を高めるために、自己再構成モデルをマスキング学習方式で訓練して重要な心臓信号に集中するようにする。これは、モデルがデータの脈絡を理解し、主要特徴を学習することができるようにして、QRS Complexなどの核心信号を正確に把握して復元することができるようにする。
自己再構成モデルは、信号全体を再現することを目的とし、特定の有意味な範囲の信号をマスキングした後、そのマスキングされた部分を復元することに焦点を合わせる。具体的な例として、心電図信号の全体範囲で有意な部分をマスキングし、このマスキングされた領域内でQRS Complexが3つ程度含まれるようにする。自己再構成モデルは、この領域を集中的に学習し、損失関数は、この領域に対してのみ計算される。このような方式で、モデルは、重要な信号を中心にデータの特徴を学習し、これによりデータの全体的な脈絡をより正確に理解し再現することができる。
その結果、心電図データから発生しうるノイズや不完全な部分にも頑強に作動し、心臓疾患の早期診断と評価に必要な高信頼性のデータ分析を可能にする。このような接近法は、心電図信号の分析における重要な情報の損失を最小限に抑え、正確な心臓状態評価を支援して医療データの解釈を大幅に向上させることができる。
一実施形態によれば、自己再構成モデル800は、少なくとも1つの隠し層(hidden layer)を含むことができ、奇数の隠し層が入出力層の間に配置されることができる。それぞれの層のノードの数は、入力層のノードの数からボトルネック層(エンコード)と呼ばれる中間層まで縮小され、その後、ボトルネック層から出力層(入力層と対称)への縮小と対称的に拡張されてもよい。この場合、図13では、エンコーダ810の層とデコーダ820の層とが対称であることが示されているが、本発明は、これに限定されず、エンコーダ810の層とデコーダ820の層のノードは対称であってもそうでなくてもよい。自己再構成モデル800は、非線形次元減少を行うことができる。入力層及び出力層の数は、入力データの前処理後に残った入力データの項目の数と対応することができる。自己再構成モデル800の構造でエンコーダ810に含まれた隠れ層のノードの数は、入力層から離れるほど減少する構造を有することができる。ボトルネック層(エンコーダとデコーダとの間に位置する最も少ないノードを有する層)のノードの数は、小さすぎる場合、十分な量の情報が伝達されない可能性があるので、特定の数以上(例えば、入力層の半分以上など)に維持されることもできる。
実施形態において、コンピューティング装置100は、第1学習データセットに対応する複数のROI信号を介して自己再構成モデルを学習させ、これにより、モデルは心電図データの重要な特徴及びパターンを学習する。学習された自己再構成モデルは、複数の心電図データのそれぞれに適用され、入力されたデータの一部をマスキングし、これを復元する過程を経て、心電図信号の主な特徴を抽出し分析する。
まとめると、コンピューティング装置100は、心電図データの分析及び解釈において高い信頼性を確保するために自己再構成モデルのマスキング学習方式で学習を誘導する。この過程で、自己再構成モデルは、QRS Complexなどの主要信号に集中して、心臓拍動の固有の分布を学習する。自己再構成モデルが心電図データの原本信号を再現する間、心電図データ自体が持っている重要な特徴、例えばQRS Complexのピークなどを効果的に学習する。
次のステップS3120で、学習が完了した自己再構成モデルからエンコーダを抽出して埋め込みモデルを生成する。学習された自己再構成モデルにおけるエンコーダ部分は、高次元データの重要な特徴を圧縮し、ノイズを除去するように学習(すなわち、第1学習データセットによる事前学習)されているので、入力されるデータ(例えば、第2学習データセットに含まれた複数のROI信号)に対応して低次元潜在ベクトルを出力することができる。
実施形態によれば、生成された埋め込みモデルの場合、データの特徴を圧縮して表現するため、入力されるデータ間の類似性を潜在空間で効果的に反映することができることを特徴とすることができる。例えば、類似した入力データに対応して、埋め込みモデルは、類似した潜在空間上の潜在ベクトルを出力する。例えば、第1ROI信号と類似した第2ROI信号の第2潜在ベクトルは、第1ROI信号に対応する第1潜在ベクトルと類似した潜在空間上に位置する。すなわち、潜在空間上で近い位置に配置される。
これとは逆に、入力データが互いに異なる場合(すなわち、大きな差がある場合)、埋め込みモデルは、これらのデータを潜在空間内で互いに離れた位置に配置する。これにより、心電図データの特徴的なパターンに基づいてデータ間の類似性と差異を明確に区分することができ、これは、心臓疾患の診断及び評価過程で利用できる。
すなわち、本発明は、第1学習データセットを利用して自己再構成モデルが自己学習されるように誘導し、学習の完了した自己再構成モデルからエンコーダを抽出して埋め込みモデルを形成する。
コンピューティング装置100は、第1学習データセットを利用して埋め込みモデルを生成した後、埋め込みモデルを利用して第2学習データセットに対応する複数のROI信号のそれぞれに対応する潜在ベクトルを抽出する埋め込みを行うことができる。
より具体的には、第2学習データセットに含まれた各ROI信号は、コンピューティング装置100によって埋め込みモデルに入力として提供され、埋め込みモデルは、当該信号の主要特徴を圧縮して低次元潜在ベクトルを出力する。潜在ベクトルは、ROI信号の複雑なパターンを簡潔に表現したものであり、各ROI信号間の類似性と差異を潜在空間で反映する。例えば、第2学習データセットに含まれた心電図データのROI信号のうち、類似したパターンの信号は、潜在空間上で互いに近い位置の潜在ベクトルとして出力される。
その後、ステップ(S3300)で、コンピューティング装置100は、クラスタリングモデルを利用して潜在ベクトルに対するクラスタリングを行うことができる。
実施形態によれば、コンピューティング装置100は、クラスタリングモデルを利用して前記潜在ベクトルに対するクラスタリングを行うことができる。
一例として、クラスタリングモデルは、埋め込みモデルが抽出した潜在ベクトルを利用してクラスタリングを行うように学習されたモデルであり、KNN(K-Nearest Neighbors)モデルであり得る。
一実施形態において、クラスタリングモデルは、複数のROI信号のそれぞれに対応する潜在ベクトル間の類似度距離によるクラスタリングを行うことを特徴とすることができる。
より詳細に説明すると、クラスタリング過程では、各潜在ベクトル間の距離又は類似度を計算し、これらのベクトルがどれほど近いか或いは類似するかを判断する。ここで、距離は、一般にユークリッド距離又はコサイン類似度などの測定方法を用いて計算される。
例えば、特定のROI信号から抽出された潜在ベクトルAと、他のROI信号から抽出された潜在ベクトルBとの距離を計算して、これらの2つのベクトルがどれほど類似するかを評価する。もしAとBの距離が近い場合、これらの2つのベクトルは、同じ又は類似のパターンを含む心電図信号を表すものと見なすことができる。クラスタリングモデルは、このような距離情報に基づいて潜在ベクトルを複数のクラスターにグループ化する。
KNNモデルの場合、特定のベクトルのK個の最も近い隣接ベクトルを見つけ、それらの隣接ベクトルの属するクラスターに当該ベクトルを含ませる方式でクラスタリングを行う。例えば、K=3の場合、潜在ベクトルAの最も近い3つの近隣ベクトルを見つけ、それらの隣接ベクトルの属するクラスターのうち、多数の属するクラスターにAを割り当てる。このような過程は、すべての潜在ベクトルに対して繰り返し行われ、最終的に類似した特徴を共有するベクトルが同じクラスターにグループ化される。
この過程は、心電図データの様々なパターンを分類し、各パターンが心臓疾患の特定の類型又はリスクレベルとどのように関連するかを分析するのに有用である。例えば、特定のクラスターが心臓拍動の不規則性に関連する信号を含んでいる場合、当該クラスターは、不整脈などの疾患の指標として利用できる。
より具体的な例として、クラスタリングモデルを利用して各潜在ベクトルに対するクラスタリングを行う場合、4つのクラスターが形成できる。このクラスタリング過程では、各潜在ベクトル間の類似度に基づいてベクトルをグループ化する。
第1クラスター、第2クラスター、第3クラスター、第4クラスターの4つのクラスターがあると仮定すると、これらのクラスターは、それぞれ異なる心電図パターンを代表することができる。例えば、第1クラスターは、正常な心拍パターンを有する信号を含むことができ、第2クラスターは、不整脈などの特定の異常を有する信号を含むことができる。第3クラスターと第4クラスターは、それぞれ異なる類型の心臓疾患又は状態を反映する信号を含むことができる。前述した説明では、クラスター(又は群)が4つ形成できることを説明するが、これは、例示に過ぎず、実際のクラスターの数は、データの特性及びクラスタリングモデルの設定に応じて変わり得る。
このようなクラスタリングは、心電図データにおいて特定のパターン又は異常を自動的に識別し、各クラスターの特徴を分析して心臓疾患の早期診断及び管理に有用な情報を提供することができる。クラスタリング結果は、医療専門家が患者の状態を総合的に評価し、カスタマイズされた治療計画を樹立するのに重要なデータを提供する。
また、ステップS3400で、コンピューティング装置100は、潜在ベクトル、及び各潜在ベクトルに対応するそれぞれのクラスターに対する情報に基づいてベクトルデータベースを構築することができる。
本発明のベクトルデータベースは、埋め込みベクトルと所属クラスター情報、並びに各ベクトル間の距離に対する情報を保存し管理するためのNoSQLデータベースサーバーであり得る。
本発明のベクトルデータベースは、様々な心電図データから抽出された潜在ベクトルの分析及び管理に最適化されている。これにより、心電図信号の特徴的なパターンを体系的に保存し、必要に応じてこれを検索して活用することができる。
一例として、本発明のベクトルデータベースは、患者の心電図データに基づくリスク評価及び状態モニタリングを支援することができる。特定の患者の新しく収集された心電図データが入ると、当該データから抽出された潜在ベクトルを、ベクトルデータベースに格納された既存のベクトルと比較して類似度を評価することができる。この類似度評価を介して、患者の心電図データがどのクラスターに属するのか、及び当該クラスターが以前に観察されたパターンのどれと類似するのかを把握することができる。
また、実施形態において、コンピューティング装置100は、クラスタリング結果に対応する複数のクラスターのそれぞれに対する代表ベクトル、すなわちprincipal vectorを抽出することができる。このようなprincipal vectorは、各クラスターを代表するベクトルであり、クラスター内の潜在ベクトルの中心を表すことができる。
より具体的には、ベクトルデータベースに格納されたデータが膨大である場合、例えば、30万件の心電図データを検索するのに多くの時間がかかることがある。一般に、特定の患者の新しい心電図データが入力される場合、前処理過程で複数のROI信号に分割され、各ROI信号に対する埋め込み過程を介して80~120個の埋め込みベクトルが生成される。例えば、新しい心電図データが、データベースに格納された全てのベクトルと個別的に比較される場合、計算量が多いため非効率的であり得る。これは、多くの時間がかかり、リアルタイム分析が必要な状況では適さない。
これにより、探索時間を短縮し且つ効率化を高めるために、患者の心電図データに基づいて生成された潜在ベクトルを各クラスターのprincipal vectorと比較して類似のクラスターを迅速に識別することができる。この過程で各埋め込みベクトルとprincipal vectorとの内積を計算して類似度を評価し、この情報に基づいて当該心電図データがどのクラスターに属するかを判断する。
この方法は、データ全体を一々比較する代わりに、各クラスターを代表するprincipal vectorと比較することにより、迅速に関連クラスターを識別することを可能にする。その結果、心電図データの迅速な分析と効率的な管理が可能であり、心臓疾患のリスクを迅速に評価し、患者の状態をリアルタイムでモニタリングする上で重要な役割を果たす。
また、実施形態において、以下の過程を介して、予測対象者の心電図データに対する心臓疾患予測情報が生成できる。
実施形態において、予測対象者の心電図データに対する心臓疾患予測情報の生成方法は、予測対象者の心電図データを取得するステップと、予測対象者の心電図データに対する前処理を行うステップと、埋め込みモデルを利用して前処理の完了した心電図データに対応する複数の潜在ベクトルを生成するステップと、複数の潜在ベクトルのそれぞれを複数のクラスターのそれぞれに対応する代表ベクトルのそれぞれと比較して、複数の潜在ベクトルのそれぞれを前記複数のクラスターのうちの1つに分類するステップと、複数の潜在ベクトルのそれぞれが複数のクラスターに分類される分類結果に基づいて、心臓疾患のリスクに関する階層化情報を生成するステップと、を含むことができる。
実施形態において、予測対象者は、心臓疾患のリスクを評価し管理するための目的で心電図検査を受けた個人であり得る。これは、一般的な健康診断を受ける人々、心臓疾患の家族歴がある人々、又は以前に心臓疾患と診断されたか或いは関連症状を示した患者であり得る。予測対象者は、定期的な心電図モニタリングを介して心臓状態を継続的に評価し、心臓疾患の早期発見及び予防のために本発明のコンピューティング装置100を利用することができる。
より詳細に説明すると、予測対象者の心電図データを取得した後、このデータを前処理してノイズを除去し、信号を精製する。ここで、前処理は、前述した説明でのように、心電図データの品質を向上させるために不要な信号成分を除去し、重要な信号特徴を強調する過程を含むことができ、各心電図信号を適切な単位(所定のウィンドウサイズ)に分割して処理する過程を含むことができる。所定のウィンドウサイズは、一般に、心電図データにおいて重要な情報が含まれた区間、例えば30秒~1分程度の長さとして定義できる。この過程で、心電図データは、一定の時間間隔で分割され、それぞれの区間が独立的に分析できる。
前処理された心電図データは、埋め込みモデル(例えば、第1学習データセットを介して事前に学習された自己再構成モデルから抽出されたエンコーダ)に入力され、それぞれのデータから潜在ベクトル(又は埋め込みベクトル)が生成される。生成された潜在ベクトルは、ベクトルデータベースに格納された既存のベクトルと比較され、類似度に応じて各ベクトルが特定のクラスターに分類される。各クラスターは、心臓状態を反映する特定のパターンを持っているので、潜在ベクトルがどのクラスターに属するかによって心臓疾患のリスクが評価される。
また、コンピューティング装置100は、複数の潜在ベクトルが複数のクラスターに分類される分類結果に基づいて心臓疾患のリスクに関する階層化情報を生成することができる。より具体的には、予測対象者の心電図データの場合、前処理過程で特定の時間単位に分割されて複数のROI信号として取得され、各ROI信号が埋め込みモデルの入力として処理されることにより、複数のROI信号のそれぞれに対する潜在ベクトルが生成される。この場合、各潜在ベクトルがどのクラスターにどれほど分類されるかに対する分類結果に応じて、心臓疾患のリスクに関する予測情報が生成できる。
より具体的な例として、予測対象者の心電図データが10個の30秒のROI信号で構成された場合、当該10個のROI信号のそれぞれは、埋め込みモデルを介して10個の潜在ベクトルに変換され、当該10個のベクトルは、所定の4つのクラスターのうちの一つに分類できる。ここで、所定の各クラスターは、特定の心臓状態を示し、第1クラスターが最も高いリスクを示す群集であり、第4クラスターに行くほどリスクが低い群集であり得る。
コンピューティング装置100は、10個の入力信号が各クラスターにどのように分類されるかによって全体リスクを予測することができる。
例えば、8つの信号が第1クラスターに、2つの信号が第2クラスターに属する場合、患者の心臓疾患のリスクが非常に高いと評価できる。これに対し、5つの信号が第2クラスター、3つの信号が第3クラスター、2つの信号が第4クラスターに属する場合には、リスクが中程度と評価されることができる。もし第3クラスターに2つ、第4クラスターに8つが属する場合には、リスクが低いと評価され、第4クラスターに10個がすべて属する場合には、患者の心臓疾患のリスクは非常に低いと評価され得る。
すなわち、本発明のコンピューティング装置100は、埋め込みモデルを利用して心電図データのROI信号を分析して潜在ベクトルを取得し、各潜在ベクトルが学習データセットに基づいてクラスタリングモデルを介して所定のクラスターのうちどのクラスターに属するかを評価することにより、患者の心臓疾患のリスクを階層化された形態で予測することができる。言い換えれば、コンピューティング装置100は、埋め込みモデルによって生成された潜在ベクトルを所定のクラスターと比較して、各ベクトルの属するクラスターを把握し、この情報に基づいて患者の心臓状態を多次元的に評価する。
このような階層化されたリスク予測方式は、心電図データ分析の精度を高め、患者の心臓状態に対する総合的な理解を助ける。したがって、本発明のコンピューティング装置は、患者の心臓疾患のリスクを体系的かつ明確に階層化して提供することにより、医療専門家が適切な予防措置又は治療計画を樹立する上で重要な情報を提供することができる。これにより、心臓疾患の早期発見と適切な対応を可能にして、患者の健康管理に寄与することができる。
図14は、本発明の一実施形態に係る複数のツリー系列モデルのアンサンブル学習を介して心臓疾患リスク予測モデルを生成する過程を説明するための例示的なフローチャートを示す。
図14を参照すると、心臓疾患リスク予測モデルを生成するステップは、複数の心電図データのそれぞれに対応するユーザーメタ情報を取得するステップ(S3500)、ベクトルデータベースから複数の潜在ベクトル及び複数の潜在ベクトルに対応するクラスター情報を取得するステップ(S3600)と、HRV特性情報、ユーザーメタ情報、複数の潜在ベクトル、及び潜在ベクトルに対応するクラスター情報に基づいて複数のツリーモデルに対する学習を行い、各ツリーモデルの出力を統合するアンサンブル学習を行って心臓疾患リスク予測モデルを生成するステップ(S3700)と、を含むことができる。
より具体的には、ステップS3500で、ユーザーメタ情報を取得するステップは、患者の年齢、性別、病歴、生活習慣などの個人健康情報を収集する過程を含むことができる。このような情報は、個々人の健康状態と特性をよく反映し、心電図データの解釈に重要な影響を及ぼす。特に、このようなメタ情報は、心臓疾患のリスクをより正確に評価するために必要な重要な変数として利用される。これにより、患者の個人的な健康要因を考慮して、カスタマイズされた予測と評価が可能になる。
ステップS3600では、ベクトルデータベースから複数の潜在ベクトルと各ベクトルに対応するクラスター情報を取得する。潜在ベクトルは、心電図データの特徴を要約したものであり、クラスター情報は、当該ベクトルの属する特定の心臓状態又は疾患のパターンを示す。これにより、患者の心電図データがどのクラスターに属するかを確認し、そのクラスターが示す心臓疾患のリスクを評価することができる。
ステップS3700では、HRV特性情報、ユーザーメタ情報、潜在ベクトル及びクラスター情報に基づいて複数のツリーモデルを学習させる。各ツリーモデルは、決定木、ランダムフォレスト、勾配ブースティング木など、様々なツリーベースの機械学習アルゴリズムを用いて心臓疾患のリスクを予測する。
ここで、クラスター情報は、ベクトルデータベースから取得でき、各ベクトルの所属クラスター情報、及びベクトル間の距離に対する情報を含むことができる。例えば、クラスター情報は、各潜在ベクトルがどのクラスターに属するのかと、クラスター内でベクトル間の類似度又は差異を示す距離情報を提供する。このような情報は、特定の心電図パターンがどのクラスターに属するのか、及びこのパターンが他の心電図データと比較してどれほど類似するか或いは異なるかを評価する上で重要な役割を果たす。これにより、心電図データの特性を分析し、当該データが特定の疾患又は状態に関連しているかをより正確に把握することができる。
より具体的には、HRV特性情報、ユーザーメタ情報、潜在ベクトル及びクラスター情報に基づいて、複数のツリー系列モデルを学習させる過程は、次の通りである。
まず、HRV特性情報は、心電図データから抽出された心拍変異度指標であり、心臓自律神経系の活動性を反映する。ユーザーメタ情報は、患者の年齢、性別、病歴、生活習慣などの個人別の健康状態を示す情報である。潜在ベクトルは、心電図データの特徴を低次元空間において表現した値であり、クラスター情報は、これらのベクトルがどの心臓状態又は疾患に関連するかを示す。
これらの様々な情報は、複数のツリーベースのモデルを学習させるための入力変数として使用される。例えば、ランダムフォレストモデルは、これらの入力変数を用いて複数の決定木を構成し、各ツリーの結果を総合して最終予測を提供する。勾配ブースティング木は、連続的なモデル学習を介して各ステップで以前モデルの誤差を補正し、より精密な予測を目標とする。
この場合、クラスター情報は、特に重要であり得る。各潜在ベクトルの属するクラスターと他のベクトルとの距離情報は、心電図データの特性をより深く理解し、心臓疾患のリスクを評価するのに使用される。このような情報は、複数のツリーモデルにおいて重要な決定要素として作用し、各モデルがこの情報に基づいてデータのパターンを学習させる。
最終アンサンブル過程では、個別ツリーモデルの予測結果を統合して一つの総合的な予測結果を導出する。この過程で、各モデルの重みが調整されて最終予測の信頼度と精度を高める。アンサンブル過程は、個別モデルの弱点を補完し、強みを結合して、より頑強で一貫した予測を提供する。
結果的に、このような総合的学習とアンサンブル過程は、心電図データの多次元な特性を反映して心臓疾患のリスクを正確かつ精密に予測することを可能にする。
その後、各ツリーモデルの出力を統合するアンサンブル学習が行われ、これは、個別モデルの予測結果を結合して最終予測の精度と信頼性を高める。アンサンブル学習を介して、心臓疾患のリスクをより精密に予測することができ、これにより患者の健康状態を総合的に評価し、早期診断と予防的治療計画を樹立することができる。
すなわち、本発明は、複数のツリー系列モデルに対するアンサンブル学習を介して様々なモデルの強みを結合して単一モデルよりもさらに高い予測性能を提供し、データのノイズや変動に対するロバスト性を向上させることができる。これにより、モデルが様々な患者データに対して一貫した予測を提供することができ、予測の不確実性を減らす。
まとめると、本発明のコンピューティング装置は、心臓疾患のリスクをより正確に予測するために3つのニューラルネットワークモデルを利用することを特徴とする。第一に、埋め込みモデルは、心電図データを低次元潜在ベクトルに変換して重要な信号特徴を要約する。第二に、クラスタリングモデルは、潜在ベクトルをクラスタリングして各ベクトルの所属クラスターを決定する。最後に、心臓疾患リスク予測モデルは、HRV特性情報、ユーザーメタ情報、潜在ベクトル、クラスター情報を結合して心臓疾患リスクを予測する。このような統合された接近方式によって、コンピューティング装置は、心臓疾患のリスクをより精密かつ一貫して予測することができる。
特に、最終心臓疾患リスク予測過程で、HRV特性情報、ユーザーメタ情報、潜在ベクトル、クラスター情報の4つの変数を利用して複数のツリー系列モデルが使用され、これらのモデルは、それぞれの出力を生成する。アンサンブル学習を介して、この出力を統合して最終出力、すなわち心臓疾患のリスクを予測する結果を生成する。このような構成により、心臓疾患リスクの予測精度が向上し、様々な変数の組み合わせを介して、より多次元的かつ包括的な評価が可能となる。これは、患者の個別的な健康状態をより正確に反映し、医療スタッフが患者のリスクを明確に理解し、最適な治療計画を樹立する上で重要な役割を果たす。
これにより、最終学習完了を介して生成される心臓疾患リスク予測モデルは、患者の心電図データに基づいて将来の心臓疾患のリスクに関する予測情報を階層化して提供することができる。階層化された情報は、患者の心臓の健康状態を多次元的に評価し、個別的な治療計画を樹立する上で重要な役割を果たす。これは、医療スタッフが患者のリスクをより明確に理解し、より精密な治療及び予防的措置を講じることを可能にすることにより、心臓疾患の早期発見及び管理に寄与する。
具体的には、このシステムは、心電図データと患者のメタ情報を総合的に分析して心臓疾患の発生可能性を様々なリスクレベルに区分する。例えば、低、中、高などとリスクを階層化して提供することにより、医療スタッフが特定の患者に対してどの程度の注意を払うべきかを明確に分かることができる。これにより、高リスク群患者は、早期診断と治療を受ける可能性が高くなり、低リスク群患者は、不要な検査を減らし、効率的な健康管理が可能となる。その結果、医療資源の効率的な配分を支援し、心臓疾患による合併症を予防し、患者の生活の質を向上させる上で重要な役割を果たすことができる。
以上で説明された実施形態は、ハードウェア構成要素、ソフトウェア構成要素、及び/又はハードウェア構成要素及びソフトウェア構成要素の組み合わせで実現できる。例えば、実施形態で説明された装置、方法及び構成要素は、例えば、プロセッサ、コントローラ、ALU(arithmetic logic unit)、デジタル信号プロセッサ(digital signal processor)、マイクロコンピュータ、FPGA(field programmable gate array)、PLU(programmable logic unit)、マイクロプロセッサ、又は命令(instruction)を実行し応答することができる他の任意の装置などのように、汎用コンピュータ又は特殊目的コンピュータを用いて実現できる。処理装置は、オペレーティングシステム(OS)、及び前記オペレーティングシステム上で行われるソフトウェアアプリケーションを行うことができる。また、処理装置は、ソフトウェアの実行に応答して、データをアクセス、保存、操作、処理及び生成することもできる。理解の便宜のために、処理装置は、1つが使用されると説明された場合もあるが、当該技術分野における通常の知識を有する者は、処理装置が複数の処理要素(processing element)及び/又は複数類型の処理要素を含むことができることが分かる。例えば、処理装置は、複数のプロセッサ又は1つのプロセッサ及び1つのコントローラを含むことができる。また、並列プロセッサ(parallel processor)などの、他の処理構成(processing configuration)も可能である。
ソフトウェアは、コンピュータプログラム(computer program)、コード(code)、命令(istruction)、又はこれらのうちの1つ以上の組み合わせを含むことができ、所望の通りに動作するように処理装置を構成するか、或いは独立的又は結合的に(collectively)処理装置を命令することができる。ソフトウェア及び/又はデータは、処理装置によって解釈されるか、或いは処理装置に命令又はデータを提供するために、任意のタイプの機械、構成要素(component)、物理的装置、仮想装置(virtual equipment)、コンピュータストレージメディア又は装置に保存できる。ソフトウェアは、ネットワークに接続されたコンピュータシステム上に分散され、分散された方法で保存又は実行されることもできる。ソフトウェア及びデータは、コンピュータ可読記録媒体に保存できる。
実施形態に係る方法は、様々なコンピュータ手段を介して行われるプログラム命令の形で実現されてコンピュータ可読媒体に記録できる。コンピュータ可読媒体は、プログラム命令、データファイル、データ構造などを単独で又は組み合わせて保存することができ、媒体に記録されるプログラム命令は、実施形態のために特別に設計及び構成されたものであってもよく、コンピュータソフトウェアの当業者に公知されて使用可能なものであってもよい。コンピュータ可読記録媒体の例には、ハードディスク、フロッピーディスク、磁気テープなどの磁気媒体(magnetic media)、CD-ROM、DVDなどの光記録媒体(optical media)、フロプティカルディスク(floptical disk)などの磁気光媒体(magneto-optical media)、及びROM、RAM、フラッシュメモリーなどのプログラム命令を保存及び実行するように特別に構成されたハードウェア装置が含まれる。プログラム命令の例には、コンパイラによって作成されるような機械語コードだけでなく、インタープリターなどを用いてコンピュータによって実行できる高級言語コードを含む。
上述したハードウェア装置は、実施形態の動作を行うために、1つ又は複数のソフトウェアモジュールとして作動するように構成されてもよく、その逆も同様である。
以上のように、実施形態はたとえ限定された図面によって説明されたが、当該技術分野における通常の知識を有する者であれば、これに基づいて様々な技術的修正及び変形を適用することができる。例えば、説明された技術が、説明された方法とは異なる順序で行われるか、及び/又は、説明されたシステム、構造、装置、回路などの構成要素が、説明された方法とは異なる形態で組み合わせ又は結合されるか、或いは他の構成要素又は均等物によって対置又は置換されても、適切な結果が達成され得る。
したがって、他の実現、他の実施形態及び特許請求の範囲と均等なものも、後述する特許請求の範囲の範囲に属する。
100 コンピューティング装置
110 プロセッサ
120 メモリー
130 プログラム
140 ディスプレイ

Claims (12)

  1. コンピューティング装置の一つ以上のプロセッサで行われる方法において、
    複数の心電図データを取得するステップと、
    前記複数の心電図データに基づいて心臓疾患予測のためのモデルを生成するステップと、を含み、
    前記心臓疾患予測のためのモデルを生成するステップは、
    前記複数の心電図データに基づいて学習データセットを構築するステップと、
    前記学習データセットに基づいて一つ以上のネットワーク関数に対する学習を行うことにより、心臓疾患リスク予測モデルを生成するステップと、を含み、
    前記学習データセットは、互いに異なる学習目的に分類された第1学習データセット及び第2学習データセットを含み、
    前記心臓疾患リスク予測モデルは、患者の心電図データに基づいて、将来の心臓疾患のリスクに関する予測情報を階層化して提供し、
    前記学習データセットを構築するステップは、
    前記複数の心電図データに対する前処理を行うステップを含み、
    前記前処理を行うステップは、
    前記複数の心電図データに対するノイズ前処理を行うステップと、
    前記ノイズ前処理された複数の心電図データを所定のウィンドウサイズに分割して複数のROI信号を生成するステップと、
    前記複数の心電図データのリード単位でHRV特性情報を抽出するステップと、を含み、
    前記HRV特性情報は、心拍間隔間の変動性に関する指標情報を含む、
    心臓疾患予測モデルの生成方法。
  2. 前記第1学習データセットは、心電図データの特性を特性空間に変換する過程に対応する学習のためのデータを含み、
    前記第2学習データセットは、心臓リスク予測に関連するキャリブレーションのためのデータを含む、請求項に記載の心臓疾患予測モデルの生成方法。
  3. 前記心臓疾患リスク予測モデルを生成するステップは、
    前記第1学習データセットを利用したマスキングベースの自己教師あり学習を介して埋め込みモデルを生成するステップと、
    前記埋め込みモデルを利用して、前記第2学習データセットに対応する複数のROI信号のそれぞれに対応する潜在ベクトル(latent vector)を抽出する埋め込み実行ステップと、
    クラスタリングモデルを利用して前記潜在ベクトルに対するクラスタリングを行うステップと、
    前記潜在ベクトル、及び各潜在ベクトルに対応するそれぞれのクラスターに対する情報に基づいて、ベクトルデータベース(DB、DataBase)を構築するステップと、を含む、請求項に記載の心臓疾患予測モデルの生成方法。
  4. 前記埋め込みモデルを生成するステップは、
    自己再構成モデルに、前記第1学習データセットに含まれたデータを入力として処理して、前記自己再構成モデルが、前記入力されたデータと類似した出力を生成するように、自己教師あり学習(self-supervised learning)を誘導するステップと、
    前記学習の完了した自己再構成モデルから、エンコーダを抽出して前記埋め込みモデルを生成するステップと、を含み、
    前記自己再構成モデルは、
    入力されるデータの一部をマスキングし、マスキングされた部分を復元するニューラルネットワークモデルであって、次元減少ネットワーク関数及び次元復元ネットワーク関数を含む、請求項に記載の心臓疾患予測モデルの生成方法。
  5. 前記クラスタリングモデルは、
    前記複数のROI信号のそれぞれに対応する潜在ベクトル間の類似度距離によるクラスタリングを行うことを特徴とし、
    前記クラスタリングを行うステップは、
    前記クラスタリング結果に対応する複数のクラスターのそれぞれに対する代表ベクトルを抽出するステップを含む、請求項に記載心臓疾患予測モデルの生成方法。
  6. 前記方法は、
    予測対象者の心電図データを取得するステップと、
    前記予測対象者の心電図データに対する前処理を行うステップと、
    前記埋め込みモデルを利用して、前記前処理の完了した心電図データに対応する複数の潜在ベクトルを生成するステップと、
    前記複数の潜在ベクトルのそれぞれを前記複数のクラスターのそれぞれに対応する代表ベクトルのそれぞれと比較して、前記複数の潜在ベクトルのそれぞれを前記複数のクラスターのうちの1つに分類するステップと、
    前記複数の潜在ベクトルのそれぞれが前記複数のクラスターに分類される分類結果に基づいて、心臓疾患のリスクに関する階層化情報を生成するステップと、を含む、請求項に記載の心臓疾患予測モデルの生成方法。
  7. 前記心臓疾患リスク予測モデルを生成するステップは、
    前記複数の心電図データのそれぞれに対応するユーザーメタ情報を取得するステップと、
    前記ベクトルデータベースから前記複数の潜在ベクトル、及び前記複数の潜在ベクトルに対応するクラスター情報を取得するステップと、
    前記HRV特性情報、前記ユーザーメタ情報、前記複数の潜在ベクトル、及び前記潜在ベクトルに対応するクラスター情報に基づいて複数のツリーモデルに対する学習を行い、各ツリーモデルの出力を統合するアンサンブル学習を行って前記心臓疾患リスク予測モデルを生成するステップと、を含む、請求項に記載の心臓疾患予測モデルの生成方法。
  8. コンピューティング装置の一つ以上のプロセッサで行われる方法において、
    複数の心電図データを取得するステップと、
    前記複数の心電図データに基づいて心臓疾患予測のためのモデルを生成するステップと、を含み、
    前記心臓疾患予測のためのモデルを生成するステップは、
    取得された前記複数の心電図データを訓練データセット及び検証データセットに分類するステップと、
    前記訓練データセット及び検証データセットに分類された心電図データのそれぞれから所定のウィンドウサイズに分割されたROI心電図データを取得するステップと、
    前記訓練データセットに含まれた心電図データから取得された第1ROI心電図データを、患者の心臓疾患を予測するためのディープラーニングモデルに入力することに対応して、前記第1ROI心電図データのそれぞれの心臓疾患クラスを予測するように前記ディープラーニングモデルを学習させるステップと、
    前記検証データセットに含まれた心電図データから取得された第2ROI心電図データを前記学習されたディープラーニングモデルに適用して前記心臓疾患クラスを区分するための閾値を決定するステップと、をさらに含む、心臓疾患予測モデルの生成方法。
  9. 前記決定するステップは、
    前記第2ROI心電図データのそれぞれに対して同一の心電図データから分離されたROI心電図データのそれぞれに対する確率スコアを収集するステップと、
    前記収集された第2ROI心電図データのそれぞれに対する確率スコアを平均して最終確率スコアを導出するステップと、
    前記導出された最終確率スコアに基づいて前記心臓疾患クラスを区分するための閾値を微調整するステップと、を含む、請求項に記載の心臓疾患予測モデルの生成方法。
  10. 前記決定された閾値は、前記ディープラーニングモデルの重みと共に保存される、請求項に記載の心臓疾患予測モデルの生成方法。
  11. 一つ以上のインストラクションを保存するメモリーと、
    前記メモリーに保存された前記一つ以上のインストラクションを実行するプロセッサと、を含み、
    前記プロセッサは、前記1つ以上のインストラクションを実行することにより、請求項1又は8に記載の方法を行う、サーバー。
  12. ハードウェアであるコンピュータに結合され、請求項1又は8に記載の方法を行うことができるようにコンピュータ可読記録媒体に保存されたコンピュータプログラム。
JP2024140097A 2023-08-22 2024-08-21 心臓疾患予測モデルの生成方法、サーバー及びコンピュータプログラム Active JP7709232B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2025107995A JP7837100B2 (ja) 2023-08-22 2025-06-26 心臓疾患予測のためのディープラーニングモデル学習方法

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
KR10-2023-0110112 2023-08-22
KR1020230110112A KR102733909B1 (ko) 2023-08-22 2023-08-22 심장 질환 예측 방법 및 상기 방법을 수행하는 컴퓨팅 장치
KR1020240104439A KR102780189B1 (ko) 2024-08-06 2024-08-06 심장 질환 예측 모델 생성 방법, 서버 및 컴퓨터 프로그램
KR10-2024-0104439 2024-08-06

Related Child Applications (1)

Application Number Title Priority Date Filing Date
JP2025107995A Division JP7837100B2 (ja) 2023-08-22 2025-06-26 心臓疾患予測のためのディープラーニングモデル学習方法

Publications (2)

Publication Number Publication Date
JP2025031652A JP2025031652A (ja) 2025-03-07
JP7709232B2 true JP7709232B2 (ja) 2025-07-16

Family

ID=92503712

Family Applications (2)

Application Number Title Priority Date Filing Date
JP2024140097A Active JP7709232B2 (ja) 2023-08-22 2024-08-21 心臓疾患予測モデルの生成方法、サーバー及びコンピュータプログラム
JP2025107995A Active JP7837100B2 (ja) 2023-08-22 2025-06-26 心臓疾患予測のためのディープラーニングモデル学習方法

Family Applications After (1)

Application Number Title Priority Date Filing Date
JP2025107995A Active JP7837100B2 (ja) 2023-08-22 2025-06-26 心臓疾患予測のためのディープラーニングモデル学習方法

Country Status (4)

Country Link
US (2) US20250064406A1 (ja)
EP (2) EP4661028A1 (ja)
JP (2) JP7709232B2 (ja)
CN (2) CN119517423A (ja)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20240420841A1 (en) * 2023-06-15 2024-12-19 Case Western Reserve University Multi-site cross-organ calibrated deep learning (muscid) method and apparatus

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20220084679A1 (en) 2020-09-15 2022-03-17 Vuno Inc. Deep neural network pre-training method for classifying electrocardiogram (ecg) data
JP2022148631A (ja) 2021-03-24 2022-10-06 アステラス製薬株式会社 心電図解析支援装置、プログラム、心電図解析支援方法、及び心電図解析支援システム
JP2022550705A (ja) 2019-09-18 2022-12-05 テンパス・ラボズ・インコーポレイテッド Ecgベースの将来の心房細動予測システムおよび方法
JP2023025436A (ja) 2021-08-10 2023-02-22 日本光電工業株式会社 学習済みモデルの生成方法、教師データの生成装置、学習モデルの生成装置、および疾病発症リスクの予測装置
WO2023022516A1 (ko) 2021-08-17 2023-02-23 주식회사 메디컬에이아이 딥러닝 알고리즘을 기반으로 복수개의 표준 심전도 데이터를 생성하는 방법

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10602942B2 (en) 2017-08-25 2020-03-31 Cambridge Heartwear Limited Method of detecting abnormalities in ECG signals
US11147463B2 (en) 2017-11-03 2021-10-19 Samsung Electronics Co., Ltd. Method and apparatus for high accuracy photoplethysmogram based atrial fibrillation detection using wearable device
CN107837082B (zh) 2017-11-27 2020-04-24 乐普(北京)医疗器械股份有限公司 基于人工智能自学习的心电图自动分析方法和装置
EP3748647A1 (en) 2019-06-05 2020-12-09 Assistance Publique Hôpitaux de Paris Method for detecting risk of torsades de pointes
US11657921B2 (en) * 2019-09-18 2023-05-23 Tempus Labs, Inc. Artificial intelligence based cardiac event predictor systems and methods
CN112022144B (zh) 2020-09-09 2022-05-24 生物岛实验室 心电异常检测网络训练方法、心电异常预警方法及装置
CA3197070A1 (en) * 2020-09-30 2022-04-07 Cardiologs Technologies Sas Electrocardiogram processing system for detecting and/or predicting cardiac events
KR102737360B1 (ko) 2021-01-18 2024-12-04 주식회사 메디컬에이아이 인공지능을 통한 2차원 심전도의 3차원 심전도로의 변환 전처리 및 심장질환 예측 시스템
KR102798691B1 (ko) 2021-08-17 2025-04-23 주식회사 메디컬에이아이 딥러닝기반 모델 및 원칙기반 모델 통합 심전도 판독 시스템
CN113951893B (zh) 2021-12-02 2023-02-03 清华大学 一种结合深度学习和电生理知识的多导联心电信号特征点提取方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2022550705A (ja) 2019-09-18 2022-12-05 テンパス・ラボズ・インコーポレイテッド Ecgベースの将来の心房細動予測システムおよび方法
US20220084679A1 (en) 2020-09-15 2022-03-17 Vuno Inc. Deep neural network pre-training method for classifying electrocardiogram (ecg) data
JP2022148631A (ja) 2021-03-24 2022-10-06 アステラス製薬株式会社 心電図解析支援装置、プログラム、心電図解析支援方法、及び心電図解析支援システム
JP2023025436A (ja) 2021-08-10 2023-02-22 日本光電工業株式会社 学習済みモデルの生成方法、教師データの生成装置、学習モデルの生成装置、および疾病発症リスクの予測装置
WO2023022516A1 (ko) 2021-08-17 2023-02-23 주식회사 메디컬에이아이 딥러닝 알고리즘을 기반으로 복수개의 표준 심전도 데이터를 생성하는 방법

Also Published As

Publication number Publication date
EP4661028A1 (en) 2025-12-10
CN119517423A (zh) 2025-02-25
US20250064406A1 (en) 2025-02-27
JP7837100B2 (ja) 2026-03-30
US20250295365A1 (en) 2025-09-25
JP2025126288A (ja) 2025-08-28
EP4512322A1 (en) 2025-02-26
JP2025031652A (ja) 2025-03-07
CN120565094A (zh) 2025-08-29

Similar Documents

Publication Publication Date Title
Xia et al. Generative adversarial network with transformer generator for boosting ECG classification
US10869610B2 (en) System and method for identifying cardiac arrhythmias with deep neural networks
Dhyani et al. Arrhythmia disease classification utilizing ResRNN
Parvaneh et al. Analyzing single-lead short ECG recordings using dense convolutional neural networks and feature-based post-processing to detect atrial fibrillation
CN113995419A (zh) 一种基于心跳节律信号的房颤发生风险预测系统及其应用
CN116504398A (zh) 用于使用基于变换器的神经网络来进行心律失常预测的方法和系统
Prakash et al. A system for automatic cardiac arrhythmia recognition using electrocardiogram signal
KR20220077583A (ko) 설명 가능한 다중 심전도 부정맥 진단 장치 및 방법
Aziz et al. Pulse plethysmograph signal analysis method for classification of heart diseases using novel local spectral ternary patterns
JP7837100B2 (ja) 心臓疾患予測のためのディープラーニングモデル学習方法
CN120203520A (zh) 基于多尺度卷积神经网络的阻塞性睡眠呼吸暂停检测方法、系统、存储介质和电子设备
CN118902426A (zh) 心律失常的个体化风险预测方法以及相关装置
US20260047794A1 (en) Methods and systems for analyzing ecg signals using neural networks
Liu et al. Integrating ECG and PCG Signals through a Dual-Modal ViT for Coronary Artery Disease Detection
KR102780189B1 (ko) 심장 질환 예측 모델 생성 방법, 서버 및 컴퓨터 프로그램
Wen et al. A scalable hybrid model for atrial fibrillation detection
KR20260021200A (ko) 심장 질환 예측 모델 생성 방법, 서버 및 컴퓨터 프로그램
Ghosh et al. A precise deep learning-based ECG arrhythmia classification scheme using deep bidirectional capsule network classifier
KR102733909B1 (ko) 심장 질환 예측 방법 및 상기 방법을 수행하는 컴퓨팅 장치
KR102804742B1 (ko) 심전도 데이터를 이용하여 심장 질환을 예측하는 방법 및 장치
Mathur Detecting myocardial infarctions using machine learning methods
Prudvi et al. Applying machine learning techniques to find important attributes for heart failure severity assessment
Irshad et al. Meta-Learning and Deep Convolutional Recurrent Neural Network with Foraging Phero Trap Optimization for Heart Disease Prediction
Guo Deep Learning for Healthcare: Empirical Studies and Methodological Advances in Physiological Time Series Data
Hayek et al. Hybrid Deep Learning Model for QRS Detection and Fourfold Heartbeat Classification

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20240821

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20250121

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20250416

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20250523

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20250528

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20250627

R150 Certificate of patent or registration of utility model

Ref document number: 7709232

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150