JP2019195461A - 胎児心電信号処理方法及び胎児心電信号処理装置 - Google Patents
胎児心電信号処理方法及び胎児心電信号処理装置 Download PDFInfo
- Publication number
- JP2019195461A JP2019195461A JP2018091128A JP2018091128A JP2019195461A JP 2019195461 A JP2019195461 A JP 2019195461A JP 2018091128 A JP2018091128 A JP 2018091128A JP 2018091128 A JP2018091128 A JP 2018091128A JP 2019195461 A JP2019195461 A JP 2019195461A
- Authority
- JP
- Japan
- Prior art keywords
- signal
- fetal
- signals
- electrocardiogram
- channels
- 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.)
- Granted
Links
- 238000012545 processing Methods 0.000 title claims abstract description 41
- 238000003672 processing method Methods 0.000 title claims abstract description 16
- 210000003754 fetus Anatomy 0.000 title abstract description 20
- 238000000034 method Methods 0.000 claims abstract description 61
- 238000000605 extraction Methods 0.000 claims abstract description 57
- 238000012880 independent component analysis Methods 0.000 claims abstract description 49
- 230000035935 pregnancy Effects 0.000 claims abstract description 4
- 230000001605 fetal effect Effects 0.000 claims description 202
- 230000008774 maternal effect Effects 0.000 claims description 60
- 239000000284 extract Substances 0.000 claims description 16
- 238000013144 data compression Methods 0.000 claims description 14
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 claims description 12
- 230000008707 rearrangement Effects 0.000 description 13
- 238000007781 pre-processing Methods 0.000 description 9
- 238000004458 analytical method Methods 0.000 description 5
- 238000007906 compression Methods 0.000 description 4
- 230000006835 compression Effects 0.000 description 4
- 238000005259 measurement Methods 0.000 description 3
- 210000001015 abdomen Anatomy 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000006698 induction Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000036760 body temperature Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 210000002458 fetal heart Anatomy 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012806 monitoring device Methods 0.000 description 1
- 230000035900 sweating Effects 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/344—Foetal cardiography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/339—Displays specially adapted therefor
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/346—Analysis of electrocardiograms
- A61B5/349—Detecting specific parameters of the electrocardiograph cycle
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/43—Detecting, measuring or recording for evaluating the reproductive systems
- A61B5/4306—Detecting, measuring or recording for evaluating the reproductive systems for evaluating the female reproductive systems, e.g. gynaecological evaluations
- A61B5/4343—Pregnancy and labour monitoring, e.g. for labour onset detection
- A61B5/4362—Assessing foetal parameters
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
- A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
- A61B5/6813—Specially adapted to be attached to a specific body part
- A61B5/6823—Trunk, e.g., chest, back, abdomen, hip
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7232—Signal processing specially adapted for physiological signals or for diagnostic purposes involving compression of the physiological signal, e.g. to extend the signal recording period
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/10—Pre-processing; Data cleansing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2134—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on separation criteria, e.g. independent component analysis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2503/00—Evaluating a particular growth phase or type of persons or animals
- A61B2503/02—Foetus
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/25—Bioelectric electrodes therefor
- A61B5/279—Bioelectric electrodes therefor specially adapted for particular uses
- A61B5/28—Bioelectric electrodes therefor specially adapted for particular uses for electrocardiography [ECG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/25—Bioelectric electrodes therefor
- A61B5/279—Bioelectric electrodes therefor specially adapted for particular uses
- A61B5/28—Bioelectric electrodes therefor specially adapted for particular uses for electrocardiography [ECG]
- A61B5/283—Invasive
- A61B5/288—Invasive for foetal cardiography, e.g. scalp electrodes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Pathology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Molecular Biology (AREA)
- Medical Informatics (AREA)
- Heart & Thoracic Surgery (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Cardiology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Theoretical Computer Science (AREA)
- Signal Processing (AREA)
- Data Mining & Analysis (AREA)
- Physiology (AREA)
- Psychiatry (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Physics & Mathematics (AREA)
- Pediatric Medicine (AREA)
- Pregnancy & Childbirth (AREA)
- Gynecology & Obstetrics (AREA)
- Reproductive Health (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
【課題】参照系独立成分分析法により分離された複数の信号から胎児心電成分を確実に含む信号を適切に選択できる胎児心電信号処理方法及び胎児心電信号処理装置を提供する提供すること。【解決手段】本発明の胎児心電信号処理方法は、妊娠中の母体より取得された複数チャネルの生体信号から参照系独立成分分析法により分離された信号を処理して胎児特徴表示信号を抽出する胎児特徴表示信号抽出ステップと、胎児特徴表示信号抽出ステップにより抽出された胎児特徴表示信号から母体の心電信号が出現される可能性のあるタイミングの信号を除去することで、胎児の心電信号を多く含む胎児特徴信号を得る母体心電信号除去ステップと、を備える。【選択図】 図１
Description
本発明は、母体内の胎児の心電信号を適切に取得できる胎児心電信号処理方法及び胎児心電信号処理装置に関する。
近年、非接触型（圧電素子、ミリ波レーダー、光学センサ等）や、従来とは異なる形状や素材の接触型の電極（導電繊維、腕時計型の光学センサ、小型ウェアラブル等）から、呼吸、脈、心音、体温、発汗等の生体信号を取得する機器が増加している。このような生体信号は、雑音を減衰させるため、フィルタ処理が必要であるが、バンドパスフィルタ等の周波数フィルタだけでは、目的の信号を抽出することが難しく、それぞれの信号に応じたノイズ除去処理を行う必要がある。
上述したような生体信号のうち、母体心電信号及び胎児心電信号を含む生体信号から、胎児心電信号を分離抽出する心電図信号処理方法が知られている（特許文献１参照）。
この特許文献１に記載の心電図信号処理方法では、胎児の心拍周期を検出する検出器から入力される心拍周期信号に基づいて、指定された誘導形式の胎児心電信号を分離抽出するための参照信号を生成し、妊娠中の母体より取得された複数チャネルの母体心電信号及び胎児心電信号を含む生体信号から、参照信号に基づいて、指定された誘導形式の胎児心電信号を参照系独立成分分析法により分離抽出する。そして、母体の胸部に取り付けられた電極のチャネルから母体心電信号の主な成分を、他の電極のチャネルの生体信号から取り除くことにより、胎児心電信号を分離させる。これにより、参照系独立成分分析法により分離された複数チャネルの入力信号から、胎児心電成分とその他の成分とが分離され、胎児心電成分を含む信号が生成される。
この特許文献１に記載の心電図信号処理方法では、胎児の心拍周期を検出する検出器から入力される心拍周期信号に基づいて、指定された誘導形式の胎児心電信号を分離抽出するための参照信号を生成し、妊娠中の母体より取得された複数チャネルの母体心電信号及び胎児心電信号を含む生体信号から、参照信号に基づいて、指定された誘導形式の胎児心電信号を参照系独立成分分析法により分離抽出する。そして、母体の胸部に取り付けられた電極のチャネルから母体心電信号の主な成分を、他の電極のチャネルの生体信号から取り除くことにより、胎児心電信号を分離させる。これにより、参照系独立成分分析法により分離された複数チャネルの入力信号から、胎児心電成分とその他の成分とが分離され、胎児心電成分を含む信号が生成される。
しかしながら、特許文献１の構成では、胎児心電信号は、成人の心電信号計測とは異なり、胎児に直接電極を貼ることができないので、どの部位に電極を貼るかを指定することができず、取得した信号が必ずしも計測したい向きで検出されるわけではない。また、抽出後の胎児心電成分のＲ波の向きやスペクトルの強さによって、十分な大きさの振幅が得られない場合、ノイズの影響を受けやすくなり、胎児心電成分を認識できないことがある。このような場合、解析後に誤った信号を、胎児心電成分を含む信号として選択してしまう可能性がある。
本発明は、このような事情に鑑みてなされたもので、参照系独立成分分析法により分離された複数の信号から胎児心電成分を確実に含む信号を適切に選択できる胎児心電信号処理方法及び胎児心電信号処理装置を提供することを目的とする。
本発明の胎児心電信号処理方法は、妊娠中の母体より取得された複数チャネルの生体信号から参照系独立成分分析法により分離された信号からノイズを除去して胎児特徴表示信号を抽出する胎児特徴表示信号抽出ステップと、前記胎児特徴表示信号抽出ステップにより抽出された前記胎児特徴表示信号から前記母体の心電信号が出現される可能性のあるタイミングの信号を除去することで、胎児の心電信号を多く含む胎児特徴信号を得る母体心電信号除去ステップと、を備える。
ここで、参照系独立成分分析法により分離された信号は、ノイズや母体の心電信号は、ある程度除去されているものの、その除去は完全ではない。また、胎児心電成分のＲ波の向きやスペクトルの強さによって、十分な大きさの振幅が得られない場合、ノイズの影響を受けやすくなり、胎児心電成分を認識できないことがある。このような場合、解析後に誤った信号を、胎児心電成分を含む信号として選択してしまう可能性がある。
これに対し、本発明では、参照系独立成分分析法により分離された信号からノイズを除去して胎児特徴表示信号を抽出し、抽出された胎児特徴表示信号から母体の心電信号が含まれる可能性のあるタイミングの信号を除去することで、ノイズ及び母体の心電信号を確実に除去することができる。これにより、胎児の心電信号を多く含む胎児特徴信号を得ることができ、各チャネルの胎児特徴信号を比べることによって、参照系独立成分分析法により分離された信号の中から胎児心電成分を確実に含む信号を適切に選択できる。
これに対し、本発明では、参照系独立成分分析法により分離された信号からノイズを除去して胎児特徴表示信号を抽出し、抽出された胎児特徴表示信号から母体の心電信号が含まれる可能性のあるタイミングの信号を除去することで、ノイズ及び母体の心電信号を確実に除去することができる。これにより、胎児の心電信号を多く含む胎児特徴信号を得ることができ、各チャネルの胎児特徴信号を比べることによって、参照系独立成分分析法により分離された信号の中から胎児心電成分を確実に含む信号を適切に選択できる。
本発明の胎児心電信号処理方法の好ましい態様としては、前記胎児特徴表示信号抽出ステップは、前記複数チャネルごとに前記参照系独立成分分析法により分離された信号の位相をずらして得た虚部信号を含む複数の複素信号を生成する複数信号生成ステップと、前記複数のチャネルごとの前記複数の複素信号の特定周波数帯域成分を抽出する特定周波数帯域成分抽出ステップと、を備え、前記母体心電信号除去ステップでは、前記複数の複素信号の特定周波数帯域成分から前記母体の心電信号が出現される可能性のあるタイミングの信号を除去するとよい。
なお、上記特定周波数帯域は、１０Ｈｚ以上４０Ｈｚ以下であることが望ましい。これは、参照系独立成分分析法により分離された信号に含まれる胎児のＲ波成分が１０〜４０Ｈｚを多く含んでいるからである。すなわち、特定周波数帯域成分抽出ステップでは、胎児のＲ波成分を多く含む周波数帯域成分を抽出している。
上記態様では、複数チャネルごとに参照系独立成分分析法により分離された信号の位相をずらして虚部信号を含む複数の複素信号を生成することにより、胎児特徴信号を抽出するためのサンプル数を増やすことができ、複数のチャネルごとに胎児特徴信号を複数得ることができる。このため、胎児特徴信号の比較対象が多くなるので、参照系独立成分分析法により分離された信号の中から胎児心電成分を確実に含む信号をより適切に選択できる。
上記態様では、複数チャネルごとに参照系独立成分分析法により分離された信号の位相をずらして虚部信号を含む複数の複素信号を生成することにより、胎児特徴信号を抽出するためのサンプル数を増やすことができ、複数のチャネルごとに胎児特徴信号を複数得ることができる。このため、胎児特徴信号の比較対象が多くなるので、参照系独立成分分析法により分離された信号の中から胎児心電成分を確実に含む信号をより適切に選択できる。
本発明の胎児心電信号処理方法の好ましい態様としては、前記複数チャネルごとの複数の前記胎児特徴信号から複数サンプルごとに値の大きいデータを含むサンプルからなる圧縮データを生成するデータ圧縮ステップと、前記データ圧縮ステップにより生成された圧縮データのうち、ノイズを除き最も大きい標準偏差の圧縮データを含む前記参照系独立成分分析法により分離された信号を前記複数チャネルごとに選択する信号選択ステップと、を備えるとよい。
ここで、解析処理においてサンプル点数の増加は、処理時間の増加に繋がるため、処理時間短縮のためには、できる限り冗長なデータを削除してから解析を実行することが望ましい。
これに対し、上記態様では、胎児特徴信号から複数サンプルごとに値の大きいデータを含むサンプルからなる圧縮データを用いて、胎児心電成分が多く含まれる参照系独立成分分析法により分離された信号を選択できるので、処理時間を短縮できる。また、ノイズを除いた最も大きい標準偏差の圧縮データを含むチャネルの信号を選択することで、胎児心電成分が多く含まれる参照系独立成分分析法により分離された信号をさらに確実に選択できる。
ここで、解析処理においてサンプル点数の増加は、処理時間の増加に繋がるため、処理時間短縮のためには、できる限り冗長なデータを削除してから解析を実行することが望ましい。
これに対し、上記態様では、胎児特徴信号から複数サンプルごとに値の大きいデータを含むサンプルからなる圧縮データを用いて、胎児心電成分が多く含まれる参照系独立成分分析法により分離された信号を選択できるので、処理時間を短縮できる。また、ノイズを除いた最も大きい標準偏差の圧縮データを含むチャネルの信号を選択することで、胎児心電成分が多く含まれる参照系独立成分分析法により分離された信号をさらに確実に選択できる。
本発明の胎児心電信号処理装置は、妊娠中の母体より取得された複数チャネルの生体信号から参照系独立成分分析法により分離された信号からノイズを除去して胎児特徴表示信号を抽出する胎児特徴表示信号抽出部と、前記胎児特徴表示信号抽出部により抽出された前記胎児特徴表示信号から前記母体の心電信号が出現される可能性のあるタイミングの信号を除去することで、胎児の心電信号を多く含む胎児特徴信号を得る母体心電信号除去部と、を備える。
本発明の胎児心電信号処理装置の好ましい態様としては、前記胎児特徴表示信号抽出部は、前記複数チャネルごとに前記参照系独立成分分析法により分離された信号の位相をずらして得た虚部信号を含む複数の複素信号を生成する複数信号生成部と、前記複数のチャネルごとの前記複数の複素信号の特定周波数帯域成分を抽出する特定周波数帯域成分抽出部と、を備え、前記母体心電信号除去部は、前記複数の複素信号の特定周波数帯域成分から前記母体の心電信号が出現される可能性のあるタイミングの信号を除去するとよい。
本発明の胎児心電信号処理装置の好ましい態様としては、前記複数チャネルごとの複数の前記胎児特徴信号から複数サンプルごとに値の大きいデータを含むサンプルからなる圧縮データを生成するデータ圧縮部と、前記データ圧縮部により生成された圧縮データのうち、ノイズを除き最も大きい標準偏差の圧縮データを含む前記参照系独立成分分析法により分離された信号を前記複数チャネルごとに選択する信号選択部と、を備えるとよい。
本発明によれば、参照系独立成分分析法により分離された複数の信号から胎児心電成分を確実に含む信号を適切に選択できる。
以下、本発明の胎児心電信号処理装置を備える医療機器及び胎児心電信号処理方法の実施形態について、図面を参照して説明する。
［医療機器の概略構成］
本実施形態の医療機器１は、いわゆる分娩監視装置により構成され、図２に示すように、妊娠中の母体Ｍに装着される複数の電極２と、複数の電極２により取得された母体の心電信号や胎児の心電信号を含む生体信号が送信される装置本体３と、を備えている。例えば、本実施形態では、胎児の心電信号を取得するための１１個の電極２と、母体の心電信号を取得するための１個の電極２を備え、胎児の心電信号を取得するためのデータ通信路は、１１個の電極２に対応する１１のチャネルにより構成される。
この装置本体３は、複数（１１個）のチャネルから生体信号を受信して前処理を行う前処理部４と、前処理が実行された生体信号から胎児成分（信号）を抽出（分離）する胎児成分抽出処理部５と、胎児成分抽出処理部５により抽出された信号を加工して、胎児心電成分を多く含む信号を適切に選択して並び替えるチャネル並び替え部６と、チャネル並び替え部６により並び替えられた信号の波形データ（胎児心電図に相当する信号）を表示するモニタ７と、を備えている。
［医療機器の概略構成］
本実施形態の医療機器１は、いわゆる分娩監視装置により構成され、図２に示すように、妊娠中の母体Ｍに装着される複数の電極２と、複数の電極２により取得された母体の心電信号や胎児の心電信号を含む生体信号が送信される装置本体３と、を備えている。例えば、本実施形態では、胎児の心電信号を取得するための１１個の電極２と、母体の心電信号を取得するための１個の電極２を備え、胎児の心電信号を取得するためのデータ通信路は、１１個の電極２に対応する１１のチャネルにより構成される。
この装置本体３は、複数（１１個）のチャネルから生体信号を受信して前処理を行う前処理部４と、前処理が実行された生体信号から胎児成分（信号）を抽出（分離）する胎児成分抽出処理部５と、胎児成分抽出処理部５により抽出された信号を加工して、胎児心電成分を多く含む信号を適切に選択して並び替えるチャネル並び替え部６と、チャネル並び替え部６により並び替えられた信号の波形データ（胎児心電図に相当する信号）を表示するモニタ７と、を備えている。
これらのうち、胎児成分抽出処理部５は、上記特許文献１に記載の方法を実行する。具体的には、胎児の心拍周期を検出する検出器から入力される心拍周期信号に基づいて、指定された誘導形式の胎児心電信号を分離抽出するための参照信号を生成し、妊娠中の母体より取得された複数チャネルの母体心電信号及び胎児心電信号を含む生体信号から、参照信号に基づいて、指定された誘導形式の胎児心電信号を参照系独立成分分析法により分離抽出する。そして、母体の胸部に取り付けられた電極のチャネルから母体心電信号の主な成分を、他の電極のチャネルの生体信号から取り除くことにより、胎児心電信号を分離させる。
本実施形態では、チャネル並び替え部６が胎児成分抽出処理部５により分離された信号を処理して胎児特徴表示信号を抽出し、抽出された胎児特徴表示信号から母体の心電信号が出現される可能性のあるタイミングの信号を除去することで、胎児の心電信号を多く含む胎児特徴信号を得て、これに基づいて胎児成分抽出処理部５により分離された信号の中から胎児心電成分を確実に含む信号を選択できる点に特徴を有している。以下に詳しく説明する。
本実施形態では、チャネル並び替え部６が胎児成分抽出処理部５により分離された信号を処理して胎児特徴表示信号を抽出し、抽出された胎児特徴表示信号から母体の心電信号が出現される可能性のあるタイミングの信号を除去することで、胎児の心電信号を多く含む胎児特徴信号を得て、これに基づいて胎児成分抽出処理部５により分離された信号の中から胎児心電成分を確実に含む信号を選択できる点に特徴を有している。以下に詳しく説明する。
［複数の電極の構成］
複数の電極２は、いわゆる超音波センサにより構成され、母体Ｍの胸部や腹部に装着される。例えば、１個の電極２は、母体Ｍの心電信号を検出するため、母体Ｍの胸部に装着され、その他の１１個の電極２は、母体Ｍの腹部に装着される。これら各電極２は、妊娠中の母体Ｍの心臓からの母体の心電信号及び胎児Ｍ１の心臓からの胎児の心電信号やノイズ等を含む生体信号を検出する。そして、各電極２により検出された生体信号は、ケーブル（図示省略）を介して装置本体３の前処理部４へと送信される。
例えば、複数の電極２による生体信号の検出は、略２０分間行われる。
複数の電極２は、いわゆる超音波センサにより構成され、母体Ｍの胸部や腹部に装着される。例えば、１個の電極２は、母体Ｍの心電信号を検出するため、母体Ｍの胸部に装着され、その他の１１個の電極２は、母体Ｍの腹部に装着される。これら各電極２は、妊娠中の母体Ｍの心臓からの母体の心電信号及び胎児Ｍ１の心臓からの胎児の心電信号やノイズ等を含む生体信号を検出する。そして、各電極２により検出された生体信号は、ケーブル（図示省略）を介して装置本体３の前処理部４へと送信される。
例えば、複数の電極２による生体信号の検出は、略２０分間行われる。
［前処理部の構成］
前処理部４は、各電極２から検出された生体信号を解析するために、該生体信号をデジタル信号に変換する。具体的には、前処理部４では、各電極２から入力された生体信号を計測し、生体信号に含まれる胎児の心電信号を増幅させる。また、アナログデータ（生体信号）をデジタルデータに変換するとともに、オンライン計算可能なブロック毎にデータを蓄え、胎児成分抽出処理部５に出力する。
前処理部４は、各電極２から検出された生体信号を解析するために、該生体信号をデジタル信号に変換する。具体的には、前処理部４では、各電極２から入力された生体信号を計測し、生体信号に含まれる胎児の心電信号を増幅させる。また、アナログデータ（生体信号）をデジタルデータに変換するとともに、オンライン計算可能なブロック毎にデータを蓄え、胎児成分抽出処理部５に出力する。
［胎児成分抽出処理部の構成］
胎児成分抽出処理部５は、上述したように、前処理部４から取得した複数チャネルの生体信号から参照系独立成分分析方法により胎児心電成分を分離する。そして、胎児成分抽出処理部５により抽出された複数チャネルごとの信号は、チャネル並び替え部６に出力される。
例えば、胎児成分抽出処理部５は、前処理部４から取得した１１チャネルの生体信号から胎児心電成分を分離する。また、１チャネルごと１５０００サンプル（１０００サンプル１秒であることから、１５秒分のサンプル）ごとに信号（波形データ）が分離され、チャネル並び替え部６へと送信される。
胎児成分抽出処理部５は、上述したように、前処理部４から取得した複数チャネルの生体信号から参照系独立成分分析方法により胎児心電成分を分離する。そして、胎児成分抽出処理部５により抽出された複数チャネルごとの信号は、チャネル並び替え部６に出力される。
例えば、胎児成分抽出処理部５は、前処理部４から取得した１１チャネルの生体信号から胎児心電成分を分離する。また、１チャネルごと１５０００サンプル（１０００サンプル１秒であることから、１５秒分のサンプル）ごとに信号（波形データ）が分離され、チャネル並び替え部６へと送信される。
［チャネル並び替え部の構成］
チャネル並び替え部６は、図３に示すように、胎児成分抽出処理部５から取得した参照系独立成分分析法により分離された信号からノイズを除去して胎児特徴表示信号を抽出する胎児特徴表示信号抽出部５１と、胎児特徴表示信号抽出部により抽出された胎児特徴表示信号から母体の心電信号が出現する可能性のあるタイミングの信号を除去することで、胎児の心電信号を多く含む胎児特徴信号を得る母体心電信号除去部５２と、複数チャネルごとの複数の胎児特徴信号から複数サンプルごとに値の大きいデータを含むサンプルからなる圧縮データを生成するデータ圧縮部５３と、データ圧縮部５３により生成された圧縮データのうち、ノイズを除き最も大きい標準偏差の圧縮データを含む参照系独立成分分析法により分離された信号をチャネルごとに選択する信号選択部５４と、標準偏差の大きさ順に複数チャネルの上記信号を並べて表示する画像生成部５５と、を備えている。
チャネル並び替え部６は、図３に示すように、胎児成分抽出処理部５から取得した参照系独立成分分析法により分離された信号からノイズを除去して胎児特徴表示信号を抽出する胎児特徴表示信号抽出部５１と、胎児特徴表示信号抽出部により抽出された胎児特徴表示信号から母体の心電信号が出現する可能性のあるタイミングの信号を除去することで、胎児の心電信号を多く含む胎児特徴信号を得る母体心電信号除去部５２と、複数チャネルごとの複数の胎児特徴信号から複数サンプルごとに値の大きいデータを含むサンプルからなる圧縮データを生成するデータ圧縮部５３と、データ圧縮部５３により生成された圧縮データのうち、ノイズを除き最も大きい標準偏差の圧縮データを含む参照系独立成分分析法により分離された信号をチャネルごとに選択する信号選択部５４と、標準偏差の大きさ順に複数チャネルの上記信号を並べて表示する画像生成部５５と、を備えている。
これらのうち、胎児特徴表示信号抽出部５１は、複数チャネルごとに参照系独立成分分析法により分離された信号の位相をずらして虚部信号を含む複数の複素信号を生成する複数信号生成部５１１と、複数チャネルごとの複数の複素信号の特定周波数帯域成分を抽出する特定周波数帯域成分抽出部５１２と、特定周波数帯域成分を絶対値に変換するとともに標準偏差σを算出して該特定周波数帯域成分を標準化して、胎児特徴表示信号を抽出する抽出部５１３と、を備えている。
［複数信号生成部の構成］
複数信号生成部５１１は、複数チャネルごとに参照系独立成分分析法により分離された信号の位相をずらして虚部信号を含む複数の複素信号を生成する。具体的には、複数信号生成部５１１は、複数のチャネルＤ１〜Ｄ１１の参照系独立成分分析法により分離された信号に対してヒルベルト変換を実行して、上記信号の位相を９０°遅らせる処理を実行し、上記信号の虚部（図４に示すｉｍａｇｅ）をそれぞれ生成する。
また、複数信号生成部５１１は、上記信号の虚部信号を含む複素信号（図４に示すｒｅａｌ±ｉｍａｇｅ）を生成する。この信号の虚部信号を含む複素信号は、信号の位相を４５°遅らせた信号と、信号の位相を４５°進めた信号とからなる。このため、本実施形態では、複数のチャネルＤ１〜Ｄ１１ごとに虚部信号を含む４つの複素信号が生成される。すなわち、複数信号生成部５１１は、複数のチャネルＤ１〜Ｄ１１の信号をそれぞれ４倍に増やすことにより、４つの複素信号を生成する。これにより、４４ｃｈ×１５０００サンプルの複素信号が生成される。
この虚部は、実部に対して大きなＲ波が表れている場合があり、胎児の特徴を抽出する際の判断の材料に用いるためである。すなわち、本処理では、信号のピークの高さや形状を求めることを目的としておらず、胎児心電成分を目立たせて、胎児心電成分を多く含む信号を抽出して選択することを目的としている。このため、４つの複素信号を生成している。
なお、本実施形態では、虚部信号を含む複素信号は、４倍に増やすこととしたが、これに限らず、２倍であってもよいし、実部を使わなくてもよい。また、参照系独立成分分析法により分離された信号の虚部を生成するのに、ヒルベルト変換を用いたが、これに限らず、他の方法により虚部を生成してもよい。
複数信号生成部５１１は、複数チャネルごとに参照系独立成分分析法により分離された信号の位相をずらして虚部信号を含む複数の複素信号を生成する。具体的には、複数信号生成部５１１は、複数のチャネルＤ１〜Ｄ１１の参照系独立成分分析法により分離された信号に対してヒルベルト変換を実行して、上記信号の位相を９０°遅らせる処理を実行し、上記信号の虚部（図４に示すｉｍａｇｅ）をそれぞれ生成する。
また、複数信号生成部５１１は、上記信号の虚部信号を含む複素信号（図４に示すｒｅａｌ±ｉｍａｇｅ）を生成する。この信号の虚部信号を含む複素信号は、信号の位相を４５°遅らせた信号と、信号の位相を４５°進めた信号とからなる。このため、本実施形態では、複数のチャネルＤ１〜Ｄ１１ごとに虚部信号を含む４つの複素信号が生成される。すなわち、複数信号生成部５１１は、複数のチャネルＤ１〜Ｄ１１の信号をそれぞれ４倍に増やすことにより、４つの複素信号を生成する。これにより、４４ｃｈ×１５０００サンプルの複素信号が生成される。
この虚部は、実部に対して大きなＲ波が表れている場合があり、胎児の特徴を抽出する際の判断の材料に用いるためである。すなわち、本処理では、信号のピークの高さや形状を求めることを目的としておらず、胎児心電成分を目立たせて、胎児心電成分を多く含む信号を抽出して選択することを目的としている。このため、４つの複素信号を生成している。
なお、本実施形態では、虚部信号を含む複素信号は、４倍に増やすこととしたが、これに限らず、２倍であってもよいし、実部を使わなくてもよい。また、参照系独立成分分析法により分離された信号の虚部を生成するのに、ヒルベルト変換を用いたが、これに限らず、他の方法により虚部を生成してもよい。
例えば、図５（ａ）は、参照系独立成分分析法により分離された信号の実部（図４に示すｒｅａｌ）の一部を示す波形データＨ１である。この信号の実部に対してヒルベルト変換を実行すると、図５（ｂ）に示すように、虚部の波形データＨ２が生成される。そして、上記信号の虚部信号を含む波形データＨ３，Ｈ４を合成する。この波形データＨ１〜Ｈ４が合わさると、図５（ｃ）に示すようになる。
なお、これら各波形データＨ１〜Ｈ４においては、周波数の高い成分及び周波数の低い成分のいずれも図５（ｃ）に示すように１／４周期（＝９０°）ずれる（波形データの山の形や高さなどが変わる）が、周波数の低い成分は、周波数の大きい成分に比べて１周期が長いため、ずれが大きくなる。すなわち、周波数の高い成分（胎児心電成分のＲ波成分）はほとんどピークの位置が変わらない。このため、次の特定周波数帯域成分抽出部による処理で、このピークの位置の信号を抽出し易くできる。
なお、これら各波形データＨ１〜Ｈ４においては、周波数の高い成分及び周波数の低い成分のいずれも図５（ｃ）に示すように１／４周期（＝９０°）ずれる（波形データの山の形や高さなどが変わる）が、周波数の低い成分は、周波数の大きい成分に比べて１周期が長いため、ずれが大きくなる。すなわち、周波数の高い成分（胎児心電成分のＲ波成分）はほとんどピークの位置が変わらない。このため、次の特定周波数帯域成分抽出部による処理で、このピークの位置の信号を抽出し易くできる。
［特定周波数帯域成分抽出部の構成］
特定周波数帯域成分抽出部５１２は、複数チャネルＤ１〜Ｄ１１ごとの４つの複素信号の特定周波数帯域成分を抽出する。胎児心電成分のＲ波成分は、１０Ｈｚ〜４０Ｈｚ相当の周波数を多く含むため、例えば、１０〜４０Ｈｚの周波数帯域成分を抽出する。この特定周波数帯域成分抽出部５１２による処理が実行されると、周波数の高い成分（４０Ｈｚ）は残り、周波数の低い成分は、平らになる。この際、複素信号のそれぞれの周波数の高い成分は、角度（ｒｅａｌやｉｍａｇｅ等）によって周波数の高さが若干異なる場合があるが、いずれの複素信号においても、周波数の高い成分が抽出される。
特定周波数帯域成分抽出部５１２は、複数チャネルＤ１〜Ｄ１１ごとの４つの複素信号の特定周波数帯域成分を抽出する。胎児心電成分のＲ波成分は、１０Ｈｚ〜４０Ｈｚ相当の周波数を多く含むため、例えば、１０〜４０Ｈｚの周波数帯域成分を抽出する。この特定周波数帯域成分抽出部５１２による処理が実行されると、周波数の高い成分（４０Ｈｚ）は残り、周波数の低い成分は、平らになる。この際、複素信号のそれぞれの周波数の高い成分は、角度（ｒｅａｌやｉｍａｇｅ等）によって周波数の高さが若干異なる場合があるが、いずれの複素信号においても、周波数の高い成分が抽出される。
［抽出部の構成］
抽出部５１３は、特定周波数帯域成分抽出部５１２により抽出された特定周波数帯域成分を絶対値に変換するとともに標準偏差σを算出して該特定周波数帯域成分を標準化した後、胎児特徴表示信号を抽出する。具体的には、抽出部５１３は、特定周波数帯域成分の標準偏差σを算出し、該標準偏差σで割ることにより、各チャネルのノイズの高さを揃える。これにより、ノイズと比べてどのチャネルのピークが目立つかをピークの高さを比較するだけで、容易に判断できる。ここで、特定周波数帯域成分にはノイズが含まれており、このノイズの９５％が±２σの範囲内に収まる。このため、２．５σ未満は、基線の範疇とみなしている。また、胎児成分のＲ波成分（ピーク値）は、３σ〜５σの範囲内であり、１０σを超えた範囲は、大きなノイズである可能性が高い。このため、抽出部５１３は、標準偏差σを算出した後、σで割ったのち、２．５σ〜１０σを胎児特徴表示信号として抽出する。この抽出部５１３の処理により特定周波数帯域成分からノイズのほとんどを取り除いた胎児特徴表示信号が生成される。
なお、上記実施形態では、抽出部５１３は、２．５σ〜１０σを胎児特徴表示信号として抽出することとしたが、これに限らず、この範囲は適宜設定できる。
抽出部５１３は、特定周波数帯域成分抽出部５１２により抽出された特定周波数帯域成分を絶対値に変換するとともに標準偏差σを算出して該特定周波数帯域成分を標準化した後、胎児特徴表示信号を抽出する。具体的には、抽出部５１３は、特定周波数帯域成分の標準偏差σを算出し、該標準偏差σで割ることにより、各チャネルのノイズの高さを揃える。これにより、ノイズと比べてどのチャネルのピークが目立つかをピークの高さを比較するだけで、容易に判断できる。ここで、特定周波数帯域成分にはノイズが含まれており、このノイズの９５％が±２σの範囲内に収まる。このため、２．５σ未満は、基線の範疇とみなしている。また、胎児成分のＲ波成分（ピーク値）は、３σ〜５σの範囲内であり、１０σを超えた範囲は、大きなノイズである可能性が高い。このため、抽出部５１３は、標準偏差σを算出した後、σで割ったのち、２．５σ〜１０σを胎児特徴表示信号として抽出する。この抽出部５１３の処理により特定周波数帯域成分からノイズのほとんどを取り除いた胎児特徴表示信号が生成される。
なお、上記実施形態では、抽出部５１３は、２．５σ〜１０σを胎児特徴表示信号として抽出することとしたが、これに限らず、この範囲は適宜設定できる。
［母体心電信号除去部の構成］
母体心電信号除去部５２は、胎児特徴表示信号から母体の心電信号が出現される可能性のあるタイミングの信号を除去することで、胎児の心電信号を多く含む胎児特徴信号を得る。この母体の心電信号が出現される可能性のあるタイミングは、母体の胸部に装着された電極２から取得された生体信号から胎児成分抽出処理部５により予め検出された母体の心電信号の周期（図６（ａ）の大きい波形Ｈ２１及び破線Ｈ２２で示された周期）である。
母体心電信号除去部５２は、胎児特徴表示信号から母体の心電信号が出現される可能性のあるタイミングの信号を除去することで、胎児の心電信号を多く含む胎児特徴信号を得る。この母体の心電信号が出現される可能性のあるタイミングは、母体の胸部に装着された電極２から取得された生体信号から胎児成分抽出処理部５により予め検出された母体の心電信号の周期（図６（ａ）の大きい波形Ｈ２１及び破線Ｈ２２で示された周期）である。
この母体心電信号除去部５２は、胎児特徴表示信号のデータＨ２（図６（ａ）参照）において、胎児心電成分を目立たせるため、該母体の心電信号が出現される可能性のあるタイミングの信号（例えば、１５０００サンプルのうちの上記タイミングの前後５０サンプル）を「０」にする。これにより、図６（ｂ）の破線にて囲まれた範囲Ｈ２３の信号が「０」にされて（除去されて）、胎児心電成分が目立たせられたデータＨ２が生成される。これにより、仮に胎児特徴表示信号に母体の心電信号が残っている場合であっても母体の心電信号が出現される可能性のあるタイミングの前後の信号を除去することで、胎児心電成分を目立たせた胎児特徴信号を得ることが可能となる。
なお、胎児心電成分と母体の心電信号とが重なる可能性があるが、胎児心電成分も母体の心電信号とともに除去して「０」とする。また、母体の心電信号が表れていない場合も、母体の心電信号が出現される可能性があるタイミング（例えば、上記破線Ｈ２２で表されたタイミング）の前後の信号を除去して「０」とする。すなわち、母体の心電信号が表れているかいないかの判断をすることなく、母体の心電信号が出現される可能性のあるタイミングの前後の信号を「０」とする。
なお、胎児心電成分と母体の心電信号とが重なる可能性があるが、胎児心電成分も母体の心電信号とともに除去して「０」とする。また、母体の心電信号が表れていない場合も、母体の心電信号が出現される可能性があるタイミング（例えば、上記破線Ｈ２２で表されたタイミング）の前後の信号を除去して「０」とする。すなわち、母体の心電信号が表れているかいないかの判断をすることなく、母体の心電信号が出現される可能性のあるタイミングの前後の信号を「０」とする。
［データ圧縮部の構成］
データ圧縮部５３は、複数チャネルごとの複数の胎児特徴信号を複数サンプルごとに圧縮する。具体的には、各チャネル１５０００サンプルを３００サンプルに圧縮する。この際、５０サンプルごと（図７（ａ）の各線で囲まれた範囲Ｈ２４）の最大値を抽出し、この５０サンプル分をその最大値を含む１つのサンプルとする。得られた３００サンプル（×４４チャネル）のうち隣接するサンプルのデータを比較して、最も大きいサンプル以外のサンプルのデータを「０」とした圧縮データを新たに作成する。
ここで、胎児心電成分のＲ波は、２４０ｂｐｍ程度が最大値であり、４００ｂｐｍ以上はむしろノイズを検出している可能性が高いため、４００ｂｐｍであっても１５０サンプル内に胎児心電成分のＲ波成分は、１つ以上含まれることはない。このため、圧縮後の３００サンプルのデータの１サンプル（圧縮前の５０サンプル）の最大値が両隣のサンプルの最大値より小さい場合に、これを「０」として扱うことで、３サンプルに１つの割合で胎児心電成分のＲ波成分が含まれた図７（ｂ）に示すような圧縮データが生成される。
データ圧縮部５３は、複数チャネルごとの複数の胎児特徴信号を複数サンプルごとに圧縮する。具体的には、各チャネル１５０００サンプルを３００サンプルに圧縮する。この際、５０サンプルごと（図７（ａ）の各線で囲まれた範囲Ｈ２４）の最大値を抽出し、この５０サンプル分をその最大値を含む１つのサンプルとする。得られた３００サンプル（×４４チャネル）のうち隣接するサンプルのデータを比較して、最も大きいサンプル以外のサンプルのデータを「０」とした圧縮データを新たに作成する。
ここで、胎児心電成分のＲ波は、２４０ｂｐｍ程度が最大値であり、４００ｂｐｍ以上はむしろノイズを検出している可能性が高いため、４００ｂｐｍであっても１５０サンプル内に胎児心電成分のＲ波成分は、１つ以上含まれることはない。このため、圧縮後の３００サンプルのデータの１サンプル（圧縮前の５０サンプル）の最大値が両隣のサンプルの最大値より小さい場合に、これを「０」として扱うことで、３サンプルに１つの割合で胎児心電成分のＲ波成分が含まれた図７（ｂ）に示すような圧縮データが生成される。
［信号選択部の構成］
信号選択部５４は、各チャネルの圧縮データの標準偏差を算出する。上記圧縮データにおいて、胎児心電成分のＲ波成分が検出できている場合、ピーク値が多数検出でき、標準偏差の値が大きくなる。このため、複数チャネルごとに標準偏差の大きい順に並べ、各チャネルのうち、標準偏差が２番目に大きい圧縮データを選択して、該標準偏差の大きい順に並べる。
この際、標準偏差が２番目に大きい圧縮データを選択するのは、標準偏差が最も大きい圧縮データを選択すると、突発的なノイズが含まれることにより標準偏差が最も大きくなっている圧縮データを選択しないようにするためである。具体的には、複数信号生成部５１１で複素信号が４つ生成されるが、元は同じデータ（参照系独立成分分析法により分離された信号）なので、ほとんど同じ標準偏差となるが突発的に生じたノイズにより標準偏差が大きくなったデータを外し、ほとんど同じ標準偏差となる残りのデータのうちの、２番目に大きい圧縮データを選択している。
なお、信号選択部５４は、標準偏差が２番目に大きい圧縮データを選択することとしたが、これに限らず、例えば、標準偏差が１番目に大きい圧縮データを選択してもよいし、標準偏差が３番目、若しくは４番目に大きい圧縮データを選択してもよい。
信号選択部５４は、各チャネルの圧縮データの標準偏差を算出する。上記圧縮データにおいて、胎児心電成分のＲ波成分が検出できている場合、ピーク値が多数検出でき、標準偏差の値が大きくなる。このため、複数チャネルごとに標準偏差の大きい順に並べ、各チャネルのうち、標準偏差が２番目に大きい圧縮データを選択して、該標準偏差の大きい順に並べる。
この際、標準偏差が２番目に大きい圧縮データを選択するのは、標準偏差が最も大きい圧縮データを選択すると、突発的なノイズが含まれることにより標準偏差が最も大きくなっている圧縮データを選択しないようにするためである。具体的には、複数信号生成部５１１で複素信号が４つ生成されるが、元は同じデータ（参照系独立成分分析法により分離された信号）なので、ほとんど同じ標準偏差となるが突発的に生じたノイズにより標準偏差が大きくなったデータを外し、ほとんど同じ標準偏差となる残りのデータのうちの、２番目に大きい圧縮データを選択している。
なお、信号選択部５４は、標準偏差が２番目に大きい圧縮データを選択することとしたが、これに限らず、例えば、標準偏差が１番目に大きい圧縮データを選択してもよいし、標準偏差が３番目、若しくは４番目に大きい圧縮データを選択してもよい。
［画像生成部の構成］
画像生成部５５は、標準偏差が２番目に大きい圧縮データのそれぞれを標準偏差の大きさ順で、複数のチャネルの参照系独立成分分析法により分離された信号を並べてモニタ７に表示させる。具体的には、画像生成部５５は、胎児心電成分が多く含まれる順で複数チャネルの参照系独立成分分析法により分離された信号が並び替えられた画像を生成してモニタ７へと出力する。
例えば、胎児成分抽出処理部５により分離された信号の波形データ（胎児心電図に相当する信号）が複数チャネルのチャネル順に表示された画像（図８参照）では、胎児心電成分が多く含まれる信号がどの信号であるかがわかりにくいが、本実施形態では、画像生成部５５により生成された画像は、胎児心電成分が多く含まれる信号が最も上側に配置される（図９参照）。このような図９に示すような画像がモニタ７に表示されることにより、医療従事者は、胎児心電成分を多く含む信号の波形データ（胎児心電図に相当する信号）を確実に選択できることとなる。
画像生成部５５は、標準偏差が２番目に大きい圧縮データのそれぞれを標準偏差の大きさ順で、複数のチャネルの参照系独立成分分析法により分離された信号を並べてモニタ７に表示させる。具体的には、画像生成部５５は、胎児心電成分が多く含まれる順で複数チャネルの参照系独立成分分析法により分離された信号が並び替えられた画像を生成してモニタ７へと出力する。
例えば、胎児成分抽出処理部５により分離された信号の波形データ（胎児心電図に相当する信号）が複数チャネルのチャネル順に表示された画像（図８参照）では、胎児心電成分が多く含まれる信号がどの信号であるかがわかりにくいが、本実施形態では、画像生成部５５により生成された画像は、胎児心電成分が多く含まれる信号が最も上側に配置される（図９参照）。このような図９に示すような画像がモニタ７に表示されることにより、医療従事者は、胎児心電成分を多く含む信号の波形データ（胎児心電図に相当する信号）を確実に選択できることとなる。
［胎児心電信号処理方法］
以上説明したチャネル並び替え部６（胎児心電処理装置）は、図１に示す手順に沿って胎児心電処理を実行する。
まず、複数信号生成部５１１は、複数のチャネルＤ１〜Ｄ１１の参照系独立成分分析法により分離された信号に対してヒルベルト変換を実行して、上記信号の虚部（図４に示すｉｍａｇｅ）を生成する（ステップＳ１１）。そして、複素信号生成部５１１は、ステップＳ１１にて生成した虚部信号を含む４つの複素信号を複数のチャネルＤ１〜Ｄ１１ごとに生成する波形合成を実行する（ステップＳ１２）。
ステップＳ１２を実行した後、特定周波数帯域成分抽出部５１２は、複数チャネルＤ１〜Ｄ１１ごとの４つの複素信号の特定周波数帯域成分を抽出する（ステップＳ１３）。その後、抽出部５１３は、ステップＳ１３により抽出された特定周波数帯域成分を絶対値に変換する（ステップＳ１４）とともに、標準偏差σを算出する（ステップＳ１５）。そして、抽出部５１３は、標準偏差σを算出した後、σで割ったのち、２．５σ〜１０σを胎児特徴表示信号として抽出する（ステップＳ１６）。
これにより、特定周波数帯域成分からノイズのほとんどが取り除かれた胎児特徴表示信号が抽出される。すなわち、ステップＳ１１〜Ｓ１６の処理は、本発明の胎児特徴表示信号抽出ステップに相当し、ステップＳ１１及びＳ１２が複数信号生成ステップに相当し、ステップＳ１３が特定周波数帯域成分抽出ステップに相当する。
以上説明したチャネル並び替え部６（胎児心電処理装置）は、図１に示す手順に沿って胎児心電処理を実行する。
まず、複数信号生成部５１１は、複数のチャネルＤ１〜Ｄ１１の参照系独立成分分析法により分離された信号に対してヒルベルト変換を実行して、上記信号の虚部（図４に示すｉｍａｇｅ）を生成する（ステップＳ１１）。そして、複素信号生成部５１１は、ステップＳ１１にて生成した虚部信号を含む４つの複素信号を複数のチャネルＤ１〜Ｄ１１ごとに生成する波形合成を実行する（ステップＳ１２）。
ステップＳ１２を実行した後、特定周波数帯域成分抽出部５１２は、複数チャネルＤ１〜Ｄ１１ごとの４つの複素信号の特定周波数帯域成分を抽出する（ステップＳ１３）。その後、抽出部５１３は、ステップＳ１３により抽出された特定周波数帯域成分を絶対値に変換する（ステップＳ１４）とともに、標準偏差σを算出する（ステップＳ１５）。そして、抽出部５１３は、標準偏差σを算出した後、σで割ったのち、２．５σ〜１０σを胎児特徴表示信号として抽出する（ステップＳ１６）。
これにより、特定周波数帯域成分からノイズのほとんどが取り除かれた胎児特徴表示信号が抽出される。すなわち、ステップＳ１１〜Ｓ１６の処理は、本発明の胎児特徴表示信号抽出ステップに相当し、ステップＳ１１及びＳ１２が複数信号生成ステップに相当し、ステップＳ１３が特定周波数帯域成分抽出ステップに相当する。
ノイズ除去ステップ（ステップＳ１１〜Ｓ１６）の実行後、母体心電信号除去部５２は、標準化された特定周波数帯域成分から母体の心電信号が出現される可能性のあるタイミングの信号を除去して、胎児の心電信号を多く含む胎児特徴信号を得る母体心電信号除去ステップを実行する（ステップＳ１７）。
その後、データ圧縮部５３は、複数チャネルごとの複数の胎児特徴信号から複数サンプルごとに値の大きいデータを取得するデータ圧縮ステップを実行する（ステップＳ１８）。
ステップＳ１８の実行後、信号選択部５４は、各チャネルの圧縮データの標準偏差を算出し（ステップＳ１９）、複数チャネルごとに標準偏差の大きい順に並べ（ステップＳ２０）、各チャネルのうち、標準偏差が２番目に大きい圧縮データを選択して、該標準偏差の大きい順に並べる（ステップＳ２１）。これらステップＳ１８〜Ｓ２１の処理は、本発明の信号選択ステップに相当する。
そして、画像生成部５５は、胎児心電成分が多く含まれる順で複数チャネルの参照系独立成分分析法により分離された信号が並び替えられた画像を生成してモニタ７へと出力して（ステップＳ２２）、チャネル並び替え処理（胎児心電信号処理）が終了する。
その後、データ圧縮部５３は、複数チャネルごとの複数の胎児特徴信号から複数サンプルごとに値の大きいデータを取得するデータ圧縮ステップを実行する（ステップＳ１８）。
ステップＳ１８の実行後、信号選択部５４は、各チャネルの圧縮データの標準偏差を算出し（ステップＳ１９）、複数チャネルごとに標準偏差の大きい順に並べ（ステップＳ２０）、各チャネルのうち、標準偏差が２番目に大きい圧縮データを選択して、該標準偏差の大きい順に並べる（ステップＳ２１）。これらステップＳ１８〜Ｓ２１の処理は、本発明の信号選択ステップに相当する。
そして、画像生成部５５は、胎児心電成分が多く含まれる順で複数チャネルの参照系独立成分分析法により分離された信号が並び替えられた画像を生成してモニタ７へと出力して（ステップＳ２２）、チャネル並び替え処理（胎児心電信号処理）が終了する。
本実施形態では、参照系独立成分分析法により分離された信号からノイズを除去して胎児特徴表示信号を抽出し、抽出された胎児特徴表示信号から母体の心電信号が含まれる可能性のあるタイミングの信号を除去することで、ノイズ及び母体の心電信号を確実に除去することができる。これにより、胎児の心電信号を多く含む胎児特徴信号を得ることができ、各チャネルの胎児特徴信号を比べることによって、参照系独立成分分析法により分離された信号の中から胎児心電成分を確実に含む胎児心電信号を適切に選択できる。
また、本実施形態では、複数チャネルごとに参照系独立成分分析法により分離された信号の位相をずらして虚部信号を含む複数の複素信号を生成することにより、胎児特徴信号を抽出するためのサンプル数を増やすことができ、複数のチャネルごとに胎児特徴信号を複数得ることができる。このため、胎児特徴信号の比較対象が多くなるので、参照系独立成分分析法により分離された信号の中から胎児心電成分を確実に含む信号をより適切に選択できる。
さらに、本実施形態では、胎児特徴信号から複数サンプルごとに値の大きいデータを含むサンプルからなる圧縮データを用いて、胎児心電成分が多く含まれる参照系独立成分分析法により分離された信号を選択できるので、処理時間を短縮できる。また、ノイズを除いた最も大きい標準偏差の圧縮データを含むチャネルの信号を選択することで、胎児心電成分が多く含まれる参照系独立成分分析法により分離された信号をさらに確実に選択できる。
このようにリアルタイムでより正確な胎児心電図に相当する信号をモニタ７に表示できるようになったことから、分析抽出した複数チャネルの信号から胎児心電成分の特徴を多く含むチャネルを正確に選択する確率が上昇するため、モニタ７の波形データの表示時に胎児心電成分の特徴を多く含む信号（胎児心電図に相当する信号）を確認でき、計測中は判断の一助となり、計測後の再確認の手間を減少させることができる。
また、胎児心電成分の特徴を多く含むチャネルを選択できる可能性が高まることから、１拍ごとのＲ波の位置を認識できる可能性が高まり、これをリファレンスや予測モデルとして再解析することでＲ波のピークの検出精度向上やより正確な心電図成分の推測が可能となる。また、これに加えて１拍ごとの胎児心拍数の算出が正確にできる。
また、胎児心電成分の特徴を多く含むチャネルを選択できる可能性が高まることから、１拍ごとのＲ波の位置を認識できる可能性が高まり、これをリファレンスや予測モデルとして再解析することでＲ波のピークの検出精度向上やより正確な心電図成分の推測が可能となる。また、これに加えて１拍ごとの胎児心拍数の算出が正確にできる。
また、本発明は上記実施形態に限定されるものではなく、本発明の趣旨を逸脱しない範囲において種々の変更を加えることが可能である。
１ 医療機器
２ 電極
３ 装置本体
４ 前処理部
５ 胎児成分抽出処理部
６ チャネル並び替え部（胎児心電信号処理装置）
７ モニタ
５１ 胎児特徴表示信号抽出部
５２ 母体心電信号除去部
５３ データ圧縮部
５４ 信号選択部
５５ 画像生成部
５１１ 複数信号生成部
５１２ 特定周波数帯域成分抽出部
５１３ 抽出部
２ 電極
３ 装置本体
４ 前処理部
５ 胎児成分抽出処理部
６ チャネル並び替え部（胎児心電信号処理装置）
７ モニタ
５１ 胎児特徴表示信号抽出部
５２ 母体心電信号除去部
５３ データ圧縮部
５４ 信号選択部
５５ 画像生成部
５１１ 複数信号生成部
５１２ 特定周波数帯域成分抽出部
５１３ 抽出部
本発明の胎児心電信号処理方法は、胎児の心拍周期信号に基づいて指定された誘導形式の胎児心電信号を分離抽出するための参照信号を生成し、妊娠中の母体より取得された複数チャネルの生体信号から、前記参照信号に基づいて、参照系独立成分分析法により分離された信号からノイズを除去して胎児特徴表示信号を抽出する胎児特徴表示信号抽出ステップと、前記胎児特徴表示信号抽出ステップにより抽出された前記胎児特徴表示信号から前記母体の心電信号が出現される可能性のあるタイミングの信号を除去することで、胎児の心電信号を多く含む胎児特徴信号を得る母体心電信号除去ステップと、を備える。
本発明の胎児心電信号処理装置は、胎児の心拍周期信号に基づいて指定された誘導形式の胎児心電信号を分離抽出するための参照信号を生成し、妊娠中の母体より取得された複数チャネルの生体信号から、前記参照信号に基づいて、参照系独立成分分析法により分離された信号からノイズを除去して胎児特徴表示信号を抽出する胎児特徴表示信号抽出部と、前記胎児特徴表示信号抽出部により抽出された前記胎児特徴表示信号から前記母体の心電信号が出現される可能性のあるタイミングの信号を除去することで、胎児の心電信号を多く含む胎児特徴信号を得る母体心電信号除去部と、を備える。
Claims (6)
- 妊娠中の母体より取得された複数チャネルの生体信号から参照系独立成分分析法により分離された信号からノイズを除去して胎児特徴表示信号を抽出する胎児特徴表示信号抽出ステップと、
前記胎児特徴表示信号抽出ステップにより抽出された前記胎児特徴表示信号から前記母体の心電信号が出現される可能性のあるタイミングの信号を除去することで、胎児の心電信号を多く含む胎児特徴信号を得る母体心電信号除去ステップと、を備えることを特徴とする胎児心電信号処理方法。 - 前記胎児特徴表示信号抽出ステップは、前記複数チャネルごとに前記参照系独立成分分析法により分離された信号の位相をずらして得た虚部信号を含む複数の複素信号を生成する複数信号生成ステップと、前記複数のチャネルごとの前記複数の複素信号の特定周波数帯域成分を抽出する特定周波数帯域成分抽出ステップと、を備え、
前記母体心電信号除去ステップでは、前記複数の複素信号の特定周波数帯域成分から前記母体の心電信号が出現される可能性のあるタイミングの信号を除去することを特徴とする請求項１に記載に胎児心電信号処理方法。 - 前記複数チャネルごとの複数の前記胎児特徴信号から複数サンプルごとに値の大きいデータを含むサンプルからなる圧縮データを生成するデータ圧縮ステップと、
前記データ圧縮ステップにより生成された圧縮データのうち、ノイズを除き最も大きい標準偏差の圧縮データを含む前記参照系独立成分分析法により分離された信号を前記複数チャネルごとに選択する信号選択ステップと、を備えることを特徴とする請求項１又は２に記載の胎児心電信号処理方法。 - 妊娠中の母体より取得された複数チャネルの生体信号から参照系独立成分分析法により分離された信号からノイズを除去して胎児特徴表示信号を抽出する胎児特徴表示信号抽出部と、
前記胎児特徴表示信号抽出部により抽出された前記胎児特徴表示信号から前記母体の心電信号が出現される可能性のあるタイミングの信号を除去することで、胎児の心電信号を多く含む胎児特徴信号を得る母体心電信号除去部と、を備えることを特徴とする胎児心電信号処理装置。 - 前記胎児特徴表示信号抽出部は、前記複数チャネルごとに前記参照系独立成分分析法により分離された信号の位相をずらして得た虚部信号を含む複数の複素信号を生成する複数信号生成部と、前記複数のチャネルごとの前記複数の複素信号の特定周波数帯域成分を抽出する特定周波数帯域成分抽出部と、を備え、
前記母体心電信号除去部は、前記複数の複素信号の特定周波数帯域成分から前記母体の心電信号が出現される可能性のあるタイミングの信号を除去することを特徴とする請求項４に記載に胎児心電信号処理装置。 - 前記複数チャネルごとの複数の前記胎児特徴信号から複数サンプルごとに値の大きいデータを含むサンプルからなる圧縮データを生成するデータ圧縮部と、
前記データ圧縮部により生成された圧縮データのうち、ノイズを除き最も大きい標準偏差の圧縮データを含む前記参照系独立成分分析法により分離された信号を前記複数チャネルごとに選択する信号選択部と、を備えることを特徴とする請求項４又は５に記載の胎児心電信号処理装置。
Priority Applications (6)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2018091128A JP6649982B2 (ja) | 2018-05-10 | 2018-05-10 | 胎児心電信号処理方法及び胎児心電信号処理装置 |
| KR1020207035251A KR20210008849A (ko) | 2018-05-10 | 2019-04-25 | 태아 심전 신호 처리 방법 및 태아 심전 신호 처리 장치 |
| PCT/JP2019/017759 WO2019216251A1 (ja) | 2018-05-10 | 2019-04-25 | 胎児心電信号処理方法及び胎児心電信号処理装置 |
| US17/052,424 US20210161455A1 (en) | 2018-05-10 | 2019-04-25 | Fetal electrocardiographic signal processing method and fetal electrocardiographic signal processing device |
| CN201980025110.4A CN112040862A (zh) | 2018-05-10 | 2019-04-25 | 胎儿心电信号处理方法及胎儿心电信号处理装置 |
| EP19800659.5A EP3791784A4 (en) | 2018-05-10 | 2019-04-25 | FETAL ELECTROCARDIOGRAPHIC SIGNAL PROCESSING METHOD AND FETAL ELECTROCARDIOGRAPHIC SIGNAL PROCESSING DEVICE |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2018091128A JP6649982B2 (ja) | 2018-05-10 | 2018-05-10 | 胎児心電信号処理方法及び胎児心電信号処理装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2019195461A true JP2019195461A (ja) | 2019-11-14 |
| JP6649982B2 JP6649982B2 (ja) | 2020-02-19 |
Family
ID=68468270
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2018091128A Active JP6649982B2 (ja) | 2018-05-10 | 2018-05-10 | 胎児心電信号処理方法及び胎児心電信号処理装置 |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US20210161455A1 (ja) |
| EP (1) | EP3791784A4 (ja) |
| JP (1) | JP6649982B2 (ja) |
| KR (1) | KR20210008849A (ja) |
| CN (1) | CN112040862A (ja) |
| WO (1) | WO2019216251A1 (ja) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113662557B (zh) * | 2021-09-16 | 2022-05-03 | 电子科技大学 | 一种基于因子同步性分析的多导联胎儿心电提取装置 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2006080167A1 (ja) * | 2005-01-31 | 2006-08-03 | Tohoku University | 心電図信号処理方法および心電図信号処理装置 |
| WO2009110088A1 (ja) * | 2008-03-07 | 2009-09-11 | 株式会社 東北テクノアーチ | 参照系独立部分空間抽出法を用いた胎児心電図の抽出方法およびその装置 |
| WO2009110051A1 (ja) * | 2008-03-03 | 2009-09-11 | 株式会社 東北テクノアーチ | 胎動時にも有効なオンライン生体信号処理方法およびその装置 |
| WO2011033784A1 (ja) * | 2009-09-18 | 2011-03-24 | 国立大学法人東北大学 | 信号抽出装置 |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7680531B2 (en) * | 2003-10-31 | 2010-03-16 | The Board Of Trustees Of The University Of Illinois | Separation of one or more fetal heart component signals from heart signal information obtained from a pregnant female |
| US7474915B2 (en) * | 2005-07-26 | 2009-01-06 | American University Of Sharjah And Arab Science And Technology Foundation | Separating mixed signals containing a distorted signal |
| US8064991B2 (en) * | 2008-01-08 | 2011-11-22 | The General Electric Company | Method of fetal and maternal ECG identification across multiple EPOCHS |
| WO2012102040A1 (ja) * | 2011-01-25 | 2012-08-02 | 公立大学法人奈良県立医科大学 | 胎児心電位信号抽出プログラム、胎児心電位信号判別装置及びこれを用いた妊婦見守りシステム |
| US9642544B2 (en) * | 2015-03-16 | 2017-05-09 | Nuvo Group Ltd. | Systems, apparatuses and methods for sensing fetal activity |
| JP6635607B2 (ja) | 2016-12-05 | 2020-01-29 | 株式会社Ｄｉｔ | エネルギー吸収機構及び木造建物 |
| CN107693004A (zh) * | 2017-09-05 | 2018-02-16 | 广东工业大学 | 基于hilbert变换的胎儿心电提取与胎儿心率识别方法 |
-
2018
- 2018-05-10 JP JP2018091128A patent/JP6649982B2/ja active Active
-
2019
- 2019-04-25 US US17/052,424 patent/US20210161455A1/en not_active Abandoned
- 2019-04-25 WO PCT/JP2019/017759 patent/WO2019216251A1/ja active Application Filing
- 2019-04-25 KR KR1020207035251A patent/KR20210008849A/ko active Search and Examination
- 2019-04-25 EP EP19800659.5A patent/EP3791784A4/en not_active Withdrawn
- 2019-04-25 CN CN201980025110.4A patent/CN112040862A/zh active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2006080167A1 (ja) * | 2005-01-31 | 2006-08-03 | Tohoku University | 心電図信号処理方法および心電図信号処理装置 |
| WO2009110051A1 (ja) * | 2008-03-03 | 2009-09-11 | 株式会社 東北テクノアーチ | 胎動時にも有効なオンライン生体信号処理方法およびその装置 |
| WO2009110088A1 (ja) * | 2008-03-07 | 2009-09-11 | 株式会社 東北テクノアーチ | 参照系独立部分空間抽出法を用いた胎児心電図の抽出方法およびその装置 |
| WO2011033784A1 (ja) * | 2009-09-18 | 2011-03-24 | 国立大学法人東北大学 | 信号抽出装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112040862A (zh) | 2020-12-04 |
| JP6649982B2 (ja) | 2020-02-19 |
| WO2019216251A1 (ja) | 2019-11-14 |
| EP3791784A1 (en) | 2021-03-17 |
| KR20210008849A (ko) | 2021-01-25 |
| US20210161455A1 (en) | 2021-06-03 |
| EP3791784A4 (en) | 2022-01-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3430989B1 (en) | Biological information analysis device, system, program, and biological information analysis method | |
| JP5271718B2 (ja) | 複数の時間区分にわたり胎児及び産婦のｅｃｇを特定する方法 | |
| US20120130263A1 (en) | Medical decision support system | |
| JP2004033775A (ja) | 弱い生理的信号を検出する方法及び装置 | |
| EP3170449A1 (en) | Device and method to detect diabetes in a person using pulse palpation signal | |
| Arteaga-Falconi et al. | R-peak detection algorithm based on differentiation | |
| EP3678553B1 (en) | Diagnosis of pathologies using infrasonic signatures | |
| KR101593412B1 (ko) | 파형종류 도수분포를 이용한 가속도맥파 분석 장치 및 분석 방법 | |
| Mollakazemi et al. | Fetal QRS detection in noninvasive abdominal electrocardiograms using principal component analysis and discrete wavelet transforms with signal quality estimation | |
| US7831300B2 (en) | Extrapolating ICA knowledge from one epoch to another for improved fetal ECG separation | |
| US20160287183A1 (en) | Classifying a time-series signal as ventricular premature contraction | |
| JP6649982B2 (ja) | 胎児心電信号処理方法及び胎児心電信号処理装置 | |
| Vuksanovic et al. | ECG based system for arrhythmia detection and patient identification | |
| JP6403311B2 (ja) | 心拍状態解析装置 | |
| Tanasković et al. | A new algorithm for fetal heart rate detection: Fractional order calculus approach | |
| Chanwimalueang et al. | Modelling stress in public speaking: Evolution of stress levels during conference presentations | |
| Amin et al. | Features for heartbeat sound signal normal and pathological | |
| MZ et al. | Shannon energy application for detection of ECG R-peak using bandpass filter and stockwell transform methods | |
| JP2014147486A (ja) | 血圧測定装置 | |
| RU2481060C1 (ru) | Способ обработки электрокардиосигнала | |
| Goda et al. | A phonography-based method improved by hidden Markov model for fetal breathing movement detection | |
| Javier et al. | R-wave detection algorithm using integration of virtual instruments | |
| JPH1043140A (ja) | 生体情報処理装置 | |
| US8731647B2 (en) | Method and system for analyzing the cardiac activity of a patient and uses thereof | |
| EP3967237A1 (en) | Ultrasound data processor |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20190220 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20190723 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20190806 |
|
| 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: 20200114 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20200117 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6649982 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |