(57)【特許請求の範囲】
【請求項１】
心電図データを取得するステップと、
前記心電図データを処理して複数の心拍データを取得するステップと、
前記複数の心拍データのうち少なくとも２つの心拍データに基づいて、機械学習における学習用データの拡張データを生成するステップとを含む、ことを特徴とする心電図データの拡張方法。
【請求項２】
前記心電図データを処理して複数の心拍データを取得する前記ステップは、
前記心電図データを前処理するステップと、
前処理された心電図データに対して心拍識別を行って、複数の心拍データを取得するステップとを含む、ことを特徴とする請求項１に記載の方法。
【請求項３】
前記前処理は、
離散ウェーブレット変換によって前記心電図データに対してマルチレベルウェーブレット分解を行って、最低レベルの近似値をゼロに設定し、離散ウェーブレット再構成によってベースライン較正後の心電図データを取得するステップを含む、ことを特徴とする請求項２に記載の方法。
【請求項４】
前記前処理は、
定常ウェーブレット変換によって前記心電図データに対してマルチレベルウェーブレット分解を行って、詳細値に対して閾値フィルタリングを行い、定常ウェーブレット再構成によってノイズ低減後の心電図データを取得するステップを含む、ことを特徴とする請求項２または３に記載の方法。
【請求項５】
前処理された心電図データに対して心拍識別を行って複数の心拍データを取得する前記ステップは、
前処理された心電図データにおける第１誘導の誘導データに対して、定常ウェーブレット変換によってマルチレベルウェーブレット分解を行って、分解結果を取得するステップと、
前記分解結果に対して特徴識別を行って、前記第１誘導の誘導データにおけるＲ波の波ピークの位置を確定するステップと、
前記第１誘導の誘導データにおけるＲ波の波ピークの位置に基づいて、前記第１誘導及び／又は他の誘導の誘導データを複数の心拍データに分割するステップとを含む、ことを特徴とする請求項２に記載の方法。
【請求項６】
前記複数の心拍データのうち少なくとも２つの心拍データに基づいて、機械学習における学習用データの拡張データを生成する前記ステップは、
前記複数の心拍データから第１心拍データ及び第２心拍データを確定するステップと、
前記第１心拍データと第２心拍データとの少なくとも１種の誘導データが接合可能である場合、第３心拍データが生成されるように前記第１心拍データと第２心拍データとを接合するステップと、
前記第１心拍データ及び／又は第２心拍データのラベルを取得するステップと、
前記第１心拍データ及び／又は第２心拍データのラベルに基づいて、前記第３心拍データのラベルを確定し、ラベル付き第３心拍データを拡張データとして確定するステップとを含む、ことを特徴とする請求項１に記載の方法。
【請求項７】
前記第１心拍データと第２心拍データとの少なくとも１種の誘導データが接合可能である場合、第３心拍データが生成されるように前記第１心拍データと第２心拍データとを接合する前記ステップは、
誘導データの種類毎に対して、前記第１心拍データと第２心拍データとにおける当該種類の誘導データが、同じＲ波方向を有し且つＲ波後の予め設定された時間内にゼロクロス点がある条件を満たす場合、前記ゼロクロス点の位置に基づいて、前記第１心拍データと第２心拍データとにおける前記種類の誘導データを、仮心拍データの当該種類の誘導データとして接合するステップと、前記条件が満たされない場合は、前記第２心拍データの前記種類の誘導データを、仮心拍データの当該種類の誘導データとして確定するステップと、
前記仮心拍データと第２心拍データとが完全同一でない場合、前記仮心拍データを第３心拍データとして確定するステップとを含む、ことを特徴とする請求項６に記載の方法。
【請求項８】
心電図データを取得するように構成されている取得モジュールと、
前記心電図データを処理して複数の心拍データを取得するように構成されている処理モジュールと、
前記複数の心拍データのうち少なくとも２つの心拍データに基づいて、機械学習における学習用データの拡張データを生成するように構成されている生成モジュールとを含む、ことを特徴とする心電図データの拡張装置。
【請求項９】
電子機器であって、メモリとプロセッサーを含み、前記メモリは、１つ又は複数のコンピュータ命令を記憶するために使用され、前記１つ又は複数のコンピュータ命令は、請求項１〜７のいずれか一項に記載の方法を実現するように前記プロセッサーによって実行される、ことを特徴とする電子機器。
【請求項１０】
コンピュータ命令が記憶されている読み取り可能な記憶媒体であって、このコンピュータ命令は、請求項１〜７のいずれか一項に記載の方法を実現するようにプロセッサーによって実行される、ことを特徴とする読み取り可能な記憶媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は医療技術の分野に関し、具体的には、心電図データの拡張方法並びにその装置、電子機器及び媒体に関する。
【背景技術】
【０００２】
データ拡張は、深層学習によく使用されている技巧の１つであり、主に学習用データセットを増加させ、データセットをできるだけ多様にして、訓練済みモデルの汎化能力をより強くするために使用される。実際においては、すべての拡張方式が現実の学習用データに適するわけではない。データセットの特徴に応じて、どのデータ拡張方式を使用するかを決定する必要がある。現在のデータ拡張は主に水平／垂直反転、回転、拡大縮小、クロッピング、カット、平行移動、コントラスト、色ジッター、ノイズなどを含む。
【０００３】
心電図は、心電図測定装置を使用して、体表から心臓の各動きのサイクルで発生する電気活動の変化が記録されたパターンを示すものであり、各セグメントの信号には、特定の医学的意義がある。しかし、反転、回転、クロッピングなどの従来のデータ拡張方法の使用は、心電図の医学的意義を破壊し、機械学習モデルの訓練に有効な効果を果たすことができない。つまり、一般的な画像処理のデータ拡張手段は、心電信号のデータ拡張に適用することができない。
【発明の概要】
【０００４】
本開示の実施例は、関連技術における問題を解決するための心電図データ拡張方法、装置、電子機器及び媒体を提供する。
【０００５】
第１態様では、本開示の実施例は、心電図データの拡張方法を提供する。この方法は、
心電図データを取得するステップと、
前記心電図データを処理して複数の心拍データを取得するステップと、
前記複数の心拍データのうち少なくとも２つの心拍データに基づいて、機械学習における学習用データの拡張データを生成するステップとを含む。
【０００６】
第１態様を参照して、本開示の第１態様の第１実施形態では、前記心電図データを処理して複数の心拍データを取得する前記ステップは、
前記心電図データを前処理するステップと、
前処理された心電図データに対して心拍識別を行って、複数の心拍データを取得するステップとを含む。
【０００７】
第１態様の第１実施形態を参照して、本開示の第１態様の第２実施形態では、前記前処理は、
離散ウェーブレット変換によって前記心電図データに対してマルチレベルウェーブレット分解を行って、最低レベルの近似値をゼロに設定し、離散ウェーブレット再構成によってベースライン較正後の心電図データを取得するステップを含む。
【０００８】
第１態様の第１実施形態又は第２実施形態を参照して、本開示の第１態様の第３実施形態では、前記前処理は、
定常ウェーブレット変換によって前記心電図データに対してマルチレベルウェーブレット分解を行って、詳細値に対して閾値フィルタリングを行い、定常ウェーブレット再構成によってノイズ低減後の心電図データを取得するステップを含む。
【０００９】
第１態様の第１実施形態を参照して、本開示の第１態様の第４実施形態では、前処理された心電図データに対して心拍識別を行って複数の心拍データを取得する前記ステップは、
前処理された心電図データにおける第１誘導の誘導データに対して、定常ウェーブレット変換によってマルチレベルウェーブレット分解を行って、分解結果を取得するステップと、
前記分解結果に対して特徴識別を行って、前記第１誘導の誘導データにおけるＲ波の波ピークの位置を確定するステップと、
前記第１誘導の誘導データにおけるＲ波の波ピークの位置に基づいて、前記第１誘導及び／又は他の誘導の誘導データを複数の心拍データに分割するステップとを含む。
【００１０】
第１態様を参照して、本開示の第１態様の第５実施形態では、前記複数の心拍データのうち少なくとも２つの心拍データに基づいて、拡張データを生成する前記ステップは、
前記複数の心拍データから第１心拍データ及び第２心拍データを確定するステップと、
前記第１心拍データと第２心拍データとのうち少なくとも１種の誘導データが接合可能である場合、第３心拍データを生成されるように前記第１心拍データと第２心拍データとを接合するステップと、
前記第１心拍データ及び／又は第２心拍データのラベルを取得するステップと、
前記第１心拍データ及び／又は第２心拍データのラベルに基づいて、前記第３心拍データのラベルを確定し、ラベル付き第３心拍データを拡張データとして確定するステップとを含む。
【００１１】
第１態様の第５実施形態を参照して、本開示の第１態様の第６実施形態では、前記第１心拍データと第２心拍データとのうち少なくとも１種の誘導データが接合可能である場合、第３心拍データを生成されるように前記第１心拍データと第２心拍データとを接合する前記ステップは、
誘導データの種類毎に対して、前記第１心拍データと第２心拍データとにおける当該種類の誘導データが、同じＲ波方向を有し且つＲ波後の予め設定された時間内にゼロクロス点がある条件を満たす場合、前記ゼロクロス点の位置に基づいて、前記第１心拍データと第２心拍データとにおける前記種類の誘導データを、仮心拍データの当該種類の誘導データとして接合するステップと、前記条件が満たされない場合は、前記第２心拍データの前記種類の誘導データを、仮心拍データの当該種類の誘導データとして確定するステップと、
前記仮心拍データと第２心拍データとが完全同一でない場合、前記仮心拍データを第３心拍データとして確定するステップとを含む。
【００１２】
第２態様では、本開示の実施例は、心電図データ拡張装置を提供する。この装置は、
心電図データを取得するように構成されている取得モジュールと、
前記心電図データを処理して複数の心拍データを取得するように構成されている処理モジュールと、
前記複数の心拍データのうち少なくとも２つの心拍データに基づいて、機械学習における学習用データの拡張データを生成するように構成されている生成モジュールとを含む。
【００１３】
第３態様では、本開示の実施例は電子機器を提供する。この電子機器は、メモリとプロセッサーを含む。前記メモリは、１つ又は複数のコンピュータ命令を記憶するために使用され、前記１つ又は複数のコンピュータ命令は、第１態様及び第１態様の第１〜第６実施形態のいずれかに記載の方法を実現するように前記プロセッサーによって実行される。
【００１４】
第４態様では、本開示の実施例は、コンピュータ命令が記憶されているコンピュータ読み取り可能な記憶媒体を提供する。このコンピュータ命令は、第１態様及び第１態様の第１〜第６実施形態のいずれかに記載の方法を実現するようにプロセッサーによって実行される。
【００１５】
本開示の実施例によって提供される技術的解決策によれば、心電図データを取得し、前記心電図データを処理して複数の心拍データを取得し、前記複数の心拍データのうち少なくとも２つの心拍データに基づいて、機械学習における学習用データの拡張データを生成することによって、心電図データに対してデータ拡張を行うことができ、学習用データセットのデータ量を増加させることができ、そして、機械学習のモデルを十分に訓練して、訓練済みモデルのパフォーマンスを高めることができる。
【００１６】
上記の一般的な説明および以下の詳細な説明は、例示的かつ説明的なものにすぎず、本開示を限定するものではないことを理解されたい。
【図面の簡単な説明】
【００１７】
添付の図面を参照しながら、以下の非限定的な実施形態を詳細に説明することにより、本開示の他の特徴、目的及び利点はより明らかになる。
【図１】本開示の実施例に係る心電図データの拡張方法のフローチャートを示す。
【図２】本開示の実施例に係る心電図データの模式図を示す。
【図３】本開示の実施例に係るベースライン較正後の心電図データの模式図を示す。
【図４】本開示の実施例に係るノイズ低減処理後の心電図データの模式図を示す。
【図５】本開示の実施例に係る心拍識別のフローチャートを示す。
【図６Ａ】本開示の実施例に係る逆方向のＲ波の模式図を示す。
【図６Ｂ】本開示の実施例に係る逆方向のＲ波の模式図を示す。
【図７】本開示の実施例に係る心拍識別の模式図を示す。
【図８】本開示の実施例に係る拡張データを生成するフローチャートを示す。
【図９】本開示の別の実施例に係る拡張データを生成するフローチャートを示す。
【図１０】本開示の実施例に係る、心拍Ａと心拍Ｂに基づいて生成された心拍Ｃの模式図を示す。
【図１１Ａ】本開示の実施例に係る心拍Ｂと心拍Ｃの対照図を示す。
【図１１Ｂ】本開示の実施例に係る心拍Ａと心拍Ｃの対照図を示す。
【図１２】本開示の実施例に係る心電図データ拡張装置のブロック図を示す。
【図１３】本開示の実施例に係る電子機器のブロック図を示す。
【図１４】本開示の実施例に係る、心電図データの拡張方法の実現に適したコンピューターシステムの構造模式図を示す。
【発明を実施するための形態】
【００１８】
以下、本開示の例示的な実施例は、当業者がそれらを容易に実施できるように、添付の図面を参照して詳細に説明する。また、明確にするために、例示的な実施例の説明に関係のない部分は、添付の図面で省略されている。
【００１９】
本開示では、「含む」や「有する」などの用語は、本明細書に開示される特徴、数字、ステップ、行為、部材、部分又はそれらの組み合わせの存在を示すことを意図し、１つ又は複数の他の特徴、数字、ステップ、行為、部材、部分又はそれらの組み合わせの存在の可能性を排除することを意図しないことを理解されたい。
【００２０】
また、本開示の実施例と実施例における特徴は、矛盾がない場合、互いに組み合わせることができることに留意されたい。以下、本開示は、添付の図面を参照し、実施例と併せて詳細に説明する。
【００２１】
ラベル付きデータセットの構築は、膨大な仕事量が必要となる。また、人々の異なる疾患の発症率が異なり、臨床医学で直接収集されたデータセットには、カテゴリーデータ量が不均衡であるという問題が多い。データセットのデータ量が不十分であるか、異なるカテゴリーのデータ量が不均衡であると、例えば分類の正確率が低くなり、再現率が低くなるなど、機械学習アルゴリズムのパフォーマンスが低下する。
【００２２】
心電図は、心電図測定装置を使用して、体表から心臓の各動きのサイクルで発生する電気活動の変化が記録されたパターンを示すものであり、各セグメントの信号には、特定の医学的意義がある。しかし、反転、回転、クロッピングなどの従来のデータ拡張方法の使用は、心電図の医学的意義を破壊し、機械学習モデルの学習に有効な効果を果たすことができない。つまり、一般的な画像処理のデータ拡張手段は、心電信号のデータ拡張に適用することができない。
【００２３】
本開示の実施例によって提供される技術的解決策によれば、心電図データを取得し、前記心電図データを処理して複数の心拍データを取得し、前記複数の心拍データの少なくとも２つの心拍データに基づいて、機械学習における学習用データの拡張データを生成することによって、心電図データに対してデータ拡張を行うことができ、学習用データセットのデータ量を増加させることができ、そして、機械学習のモデルを十分に訓練して、訓練済みモデルのパフォーマンスを高めることができる。
【００２４】
図１は、本開示の実施例に係る心電図データの拡張方法のフローチャートを示す。
図１に示すように、この方法はステップＳ１０１〜Ｓ１０３を含む。
【００２５】
ステップＳ１０１では、心電図データを取得する。
ステップＳ１０２では、前記心電図データを処理して複数の心拍データを取得する。
ステップＳ１０３では、前記複数の心拍データの少なくとも２つの心拍データに基づいて、機械学習における学習用データの拡張データを生成する。
【００２６】
本開示の実施例によって提供される技術的解決策によれば、心電図データを取得し、前記心電図データを処理して複数の心拍データを取得し、前記複数の心拍データの少なくとも２つの心拍データに基づいて、機械学習における学習用データの拡張データを生成することによって、心電図データに対してデータ拡張を行うことができ、学習用データセットのデータ量を増加させることができ、そして、機械学習のモデルを十分に訓練して、訓練済みモデルのパフォーマンスを高めることができる。
【００２７】
図２は、本開示の実施例に係る心電図データの模式図を示す。図２に示すように、心電図データは連続的に収集され、一般的に複数の心拍を含むので、一定の周期性の特徴を有する。各心拍にある最大の波ピークはＲ波と呼ばれる。Ｒ波は顕著な特徴を有するため、心拍を分割する基準として使用され得る。
【００２８】
心電図の専門用語では、心電図を記録する際に人体の体表に取り付ける電極の配置位置、電極と増幅器の接続方式は心電図の誘導と呼ばれる。一般的に使用される誘導は１２種類あり、それぞれ第I誘導、第II誘導、第III誘導、ａＶＲ誘導、ａＶＬ誘導、ａＶＦ誘導、Ｖ１誘導、Ｖ２誘導、Ｖ３誘導、Ｖ４誘導、Ｖ５誘導及びＶ６誘導である。すべての誘導の心電信号は、ほぼ同時に検出される。それらは、Ｒ波の分布密度や位置などに関してほぼ同じである。図２で例示される心電図データは１つの誘導データのみであり、他の誘導データもそれと類似する可能性がある。
【００２９】
本開示の実施例によれば、前記心電図データを処理して複数の心拍データを取得する前記ステップは、前記心電図データを前処理し、前処理された心電図データに対して心拍識別を行って複数の心拍データを取得するステップを含む。
【００３０】
本開示の実施例によって提供される技術的解決策によれば、前記心電図データを前処理し、前処理された心電図データに対して心拍識別を行って複数の心拍データを取得することによって、心電図データに対してデータ拡張を行うことができ、学習用データのデータ量を増加させることができ、そして、機械学習のモデルを十分に訓練して、モデルの表現を高めることができる。
【００３１】
本開示の実施例によれば、前記前処理は、離散ウェーブレット変換によって前記心電図データに対してマルチレベルウェーブレット分解を行って、最低レベルの近似値をゼロに設定し、離散ウェーブレット再構成によってベースライン較正後の心電図データを取得するステップを含む。
【００３２】
ウェーブレット変換は、信号の時間周波数分析及び処理を実行する理想的なツールである。コンピュータで一般的に使用されるウェーブレット変換アルゴリズムはｍａｌｌａｔアルゴリズムである。このアルゴリズムの中心的なアイデアは、フィルタを使用してウェーブレット変換を実行することである。
【００３３】
元の入力信号Ｓは、２つの相補的フィルタ（ローパスフィルタとハイパスフィルタ）を通過してＡとＤの２つの信号を生成する。Ａは信号の近似値（低周波フィルタで取得された値）であり、Ｄは信号の詳細値（高周波フィルタで取得された値）である。このような一対のフィルタによる原信号の分解は、１レベル分解と呼ばれる。信号の分解過程を繰り返すことができる。即ち、マルチレベル分解を実行することができる。高周波数成分をさらに分解せず、低周波数成分を継続的に分解すると、より低い周波数の成分を取得し、ウェーブレット分解ツリーを形成することができる。必要に応じて、分解レベル数を決定することができる。
【００３４】
本開示の実施例では、離散ウェーブレット変換（ＤＷＴ）を採用して、各レベルの分解後にデータをダウンサンプリングするため、データ量が比較的小さく、計算が比較的速い。レベル数が増え続けると、得られる低周波数成分の解像度は徐々に低下する。対照的に、離散ウェーブレット再構成アルゴリズム（ＩＤＷＴ）は、アップサンプリングとフィルタリングの２つの過程を含む。アップサンプリングとは、ダウンサンプリングデータの間に０を挿入することを意味する。
【００３５】
ベースライン変動の除去中に使用されるウェーブレットは、ｄｂ５ウェーブレットである。図３に示すように、ＤＷＴを使用して信号に対して９レベルのウェーブレット分解を実行した後、最低レベルの近似値をゼロに設定してから、ＩＤＷＴを実行すると、ベースライン変動を除去した後の心電信号を取得することができる。
【００３６】
本開示の実施例によって提供される技術的解決策によれば、離散ウェーブレット変換によって前記心電図データに対してマルチレベルウェーブレット分解を行って、最低レベルの近似値をゼロに設定し、離散ウェーブレット再構成によってベースライン較正後の心電図データを取得することで、心電図データに対してベースライン較正をよりよく行うことができる。
【００３７】
本開示の実施例によれば、前処理過程では、ベースライン較正に加えて、ノイズ低減処理を実行することもできる。心電図データセットは、筋電干渉及び５０／６０Ｈｚの商用周波数ノイズを含む。例えば、ウェーブレット変換閾値法を使用して、心電信号の高周波ノイズを除去することができる。
【００３８】
本開示の実施例によれば、前記前処理は、定常ウェーブレット変換によって前記心電図データに対してマルチレベルウェーブレット分解を行って、詳細値に対して閾値フィルタリングを行い、定常ウェーブレット再構成によってノイズ低減後の心電図データを取得するステップを含む。
【００３９】
本開示の実施例によれば、定常ウェーブレット変換（ＳＷＴ）を採用してノイズを除去する。ＳＷＴとＤＷＴの最大の違いは、ＳＷＴは各レベル分解の後にダウンサンプリングを実行しないことである。そのため、データ量が比較的大きいが、各レベルの解像度が変化しないままである。高周波ノイズの除去は、単にあるレベルの信号を完全に削除するのではなく、閾値削除を行うことであるため、解像度は特に重要である。同じ理由で、定常ウェーブレット再構成（ＩＳＷＴ）の過程でも、アップサンプリングを実行する必要はない。
【００４０】
高周波ノイズの除去中に使用されるウェーブレットは、ｂｉｏｒ２．６ウェーブレットである。ＳＷＴを使用して信号に対して６レベルウェーブレット分解を実行した後、まず、最初の２つのレベルでの周波数の最も高い詳細値をゼロに設定する。これらの詳細値で表わされる周波数は高すぎて、ほとんどの心電信号情報を含まない。次に、３〜６レベルでの詳細値に対して閾値フィルタリングを実行する。閾値は、次の式を使用して選択される。
【数１】
0006831944.tif 000002
ｄは詳細値、Ｎはその長さ、ｍｉｄは中央値である。使用される閾値方法は、ハーフソフト閾値法であり、ハード閾値法のようにウェーブレット係数の突然変化を生成することも、ソフト閾値法のように偏差を生成することもない。ハーフソフト閾値法の式は、以下のとおりである。
【数２】
0006831944.tif 000003
式では、λ１＝λｔｈ、λ２＝１．２５λｔｈとする。閾値フィルタリングの後にＩＳＷＴを行うと、高周波ノイズを除去することができる。ノイズ低減処理後の心電図データは図４に示される。
【００４１】
本開示の実施例によって提供される技術的解決策によれば、定常ウェーブレット変換によって前記心電図データに対してマルチレベルウェーブレット分解を行い、詳細値に対して閾値フィルタリングを行い、定常ウェーブレット再構成によってノイズ低減後の心電図データを取得することで、心電図データのノイズ低減効果を向上させることができる。
【００４２】
図５は、本開示の実施例に係る心拍識別のフローチャートを示す。
【００４３】
図５に示すように、この方法はステップＳ５０１〜Ｓ５０３を含む。
【００４４】
ステップＳ５０１では、前処理された心電図データにおける第１誘導の誘導データに対して、定常ウェーブレット変換によってマルチレベルウェーブレット分解を行って、分解結果を取得する。
ステップＳ５０２では、前記分解結果に対して特徴識別を行って、前記第１誘導の誘導データにおけるＲ波の波ピークの位置を確定する。
ステップＳ５０３では、前記第１誘導の誘導データにおけるＲ波の波ピークの位置に基づいて、前記第１誘導及び／又は他の誘導の誘導データを複数の心拍データに分割する。
【００４５】
本開示の実施例によれば、ステップＳ５０１では、ＳＷＴを使用して信号に対してウェーブレット変換を行うことができる。使用されるウェーブレットは、例えば、
【数３】
0006831944.tif 000004
であってもよい。
【００４６】
本開示の実施例によれば、ステップＳ５０２では、変換された第５レベルの詳細値を選択することができる。このレベルのＲ波の特徴が最も明白である。上向きのＲ波は、このレベルで負−正の極大値ペアの形で現れる。下向きのＲ波は、このレベルで正−負の極大値ペアの形で現れる。これらの極大値ペアを見つけて、それらのゼロクロス点の位置を見つける。対応するＲ波ピークは、ゼロクロス点±０．０５ｓ内の極大／極小値である。
【００４７】
健常者の心電図においても、各種類誘導のＲ波のピーク方向は一定とは限らない。病変後の心電図では、Ｒ波のピーク方向は通常の方向と反対になる状況もある。図６Ａ及び図６Ｂは、本開示の実施例に係る逆方向のＲ波の模式図を示す。図６Ａは、病理学的に引き起こされるＲ波反転の状況の模式図を示す。図６Ｂは、逆に配置された電極によるＲ波反転の模式図を示す。
【００４８】
実際の操作では、同じ人の同時に検出された心電信号の１２種類の誘導（これらの１２種類の誘導はそれぞれ第Ｉ誘導、第II誘導、第III誘導、ａＶＲ誘導、ａＶＬ誘導、ａＶＦ誘導、Ｖ１誘導、Ｖ２誘導、Ｖ３誘導、Ｖ４誘導、Ｖ５誘導、Ｖ６誘導である。以下、上記の１２種類の誘導の言及順序に基づき、各種類の誘導を指すことができる。例えば、第１誘導は第Ｉ誘導を指し、第７誘導はＶ１誘導を指すことができる）に対して、それぞれ単独のベースライン除去及びフィルタリング操作を実行する。その後、第１誘導のみに対してＲ波識別を行い、第１誘導のＲ波識別の結果を使用して他のすべての誘導を分割する。この理由は、すべての誘導の心電信号がほぼ同時に検出されるため、それらがＲ波の分布密度や位置などに関してほぼ同じであることである。
【００４９】
分割後に、第１誘導以外の１１種類の誘導のＲ波位置を較正する必要がある。人体の隣接する細胞の間の電位伝導に一定の時間がかかるため、これらの１２種類の誘導のＲ波ピークの位置は、最大で±０．０５ｓの違いがある可能性がある。他の１１種類の誘導に対してそれぞれ次の操作を実行する。第１誘導のＲ波ピークの位置の近くの±０．０５ｓの範囲内で最大値と最小値を見つける。最大値が最小値の絶対値の三分の一よりも小さい場合、最小値のある位置はこの誘導のＲ波ピークであると認め、そうでない場合、最大値のある位置はこの誘導のＲ波ピークであると認める。従って、本開示の実施例の方法は、前記第１誘導のデータにおけるＲ波の波ピークの位置に基づいて、他の誘導のデータにおけるＲ波の波ピークの位置を較正することをさらに含み得る。
【００５０】
本開示の実施例によれば、分割時に、最初と最後の心拍信号は不完全な可能性があるため、最初と最後の心拍信号を棄却し、隣接する２つの心拍のＲピークの中間点を分割点として使用することができる。この方法により、完全な心電信号は、一つ一つの心拍に分割することができる。分割後の結果は図７に示される。
【００５１】
本開示の実施例によって提供される技術的解決策によれば、前処理された心電図データにおける第１誘導のデータに対して、定常ウェーブレット変換によってマルチレベルウェーブレット分解を行って、分解結果を取得するステップと、前記分解結果に対して特徴識別を行って、前記第１誘導のデータにおけるＲ波の波ピークの位置を確定するステップと、前記第１誘導のデータにおけるＲ波の波ピークの位置に基づいて、前記第１誘導及び／又は他の誘導のデータを複数の心拍データに分割するステップとを実行することによって、心拍データを正確に識別することができ、心電図データに対してデータ拡張を行うことができ、学習用データセットのデータ量を増加させることができ、そして、機械学習のモデルを十分に訓練して、訓練済みモデルのパフォーマンスを高めることができる。
【００５２】
図８は、本開示の実施例に係る拡張データを生成するフローチャートを示す。
【００５３】
図８に示すように、この方法はステップＳ８０１〜Ｓ８０４を含む。
【００５４】
ステップＳ８０１では、前記複数の心拍データから第１心拍データ及び第２心拍データを確定する。
ステップＳ８０２では、前記第１心拍データと第２心拍データのうち少なくとも１種類の誘導データが接合可能である場合、前記第１心拍データと第２心拍データとを接合して第３心拍データを生成する。
ステップＳ８０３では、前記第１心拍データ及び／又は第２心拍データのラベルを取得する。
ステップＳ８０４では、前記第１心拍データ及び／又は第２心拍データのラベルに基づいて、前記第３心拍データのラベルを確定し、ラベル付き第３心拍データを拡張データとして確定する。
【００５５】
本開示の実施例によれば、２つの分割された心拍信号（各誘導）を、それぞれ心拍Ａと心拍Ｂとしてランダムに選択することができる。このランダム選択は、次の３つの方法の１つに基づいたものである。
（１）同じ患者の異なる時間での２つの心拍信号をＡとＢとしてランダムに選択する。
（２）同じ種類の疾患に罹患した２人の異なる患者の２つの心拍信号をＡとＢとしてランダムに選択する。
（３）１人の健常者の心拍信号をＡとして、１人の患者の心拍信号をＢとしてランダムに選択する。
【００５６】
本開示の実施例によれば、３つの方法を組み合わせて同時に使用することによって、データの汎化程度を増加させることができる。
【００５７】
また、選択した２つの信号が接合可能であるか否かを事前に検出することもできる。例えば、心拍Ａと心拍Ｂとにおける各種誘導のＲ波のピーク方向をそれぞれ比較して、ＡとＢとにおけるある種類の誘導データのＲ波のピーク方向が同じである場合、この誘導データを接合可能としてマークし、当該ピーク方向が同じではない場合、接合不可としてマークする。すべての誘導データを接合不可としてマークする場合、この接合が失敗し、心拍を再度選択する。接合可能としてマークされた誘導データが少なくとも１種類がある場合、接合を実行することができる。
【００５８】
本開示の実施例によれば、上記の３つの方法については、病変を検出するための主な特徴データの位置を確定し、第３心拍データのこの位置でのデータが心拍Ａ又は心拍Ｂに由来するかを判断する。心拍Ａに由来する場合、心拍Ａのラベルに基づいて第３心拍データのラベルを配置する。心拍Ｂに由来する場合、心拍Ｂのラベルに基づいて第３心拍データのラベルを配置する。例えば、心筋梗塞病変の主な検出位置は、ＳＴセグメント（Ｒ波の後のセグメント）であるため、拡張された心拍ＣのＳＴセグメントが心拍Ｂに由来する場合、心拍Ｃのラベルは心拍Ｂと同じであると認める。例えば、心拍Ｂが心筋梗塞患者の心電図データに由来する場合、心拍Ｃのラベルを心筋梗塞疾患罹患として設定する。心拍Ｂが健康な人の心電図データに由来する場合、心拍Ｃのラベルを健康として設定する。
【００５９】
本開示の実施例によって提供される技術的解決策によれば、前記複数の心拍データから第１心拍データ及び第２心拍データを確定するステップと、前記第１心拍データと第２心拍データにおける少なくとも１種類の誘導データが接合可能である場合、前記第１心拍データと第２心拍データとを接合して第３心拍データを生成するステップと、前記第１心拍データ及び／又は第２心拍データのラベルを取得するステップと、前記第１心拍データ及び／又は第２心拍データのラベルに基づいて、前記第３心拍データのラベルを確定し、ラベル付き第３心拍データを拡張データとして確定するステップとを実行することによって、心電図データに対してデータ拡張を行うことができ、学習用データセットのデータ量を増加させることができ、そして、機械学習のモデルを十分に訓練して、訓練済みモデルの表現を高めることができる。
【００６０】
本開示の実施例によれば、前記第１心拍データと第２心拍データとのうち少なくとも１種類の誘導が接合可能である場合、前記第１心拍データと第２心拍データとを接合して第３心拍データを生成する前記ステップは、
誘導データの種類毎に対して、前記第１心拍データと第２心拍データとにおける当該種類の誘導データが、同じＲ波方向を有し且つＲ波後の予め設定された時間内にゼロクロス点がある条件を満たす場合、前記ゼロクロス点の位置に基づいて、前記第１心拍データと第２心拍データとにおける前記種類の誘導データを、仮心拍データの当該種類の誘導データとして接合するステップと、前記条件が満たされない場合は、前記第２心拍データの前記種類の誘導データを、仮心拍データの当該種類の誘導データとして確定するステップと、
前記仮心拍データと第２心拍データとが完全同一でない場合、前記仮心拍データを第３心拍データとして確定するステップとを含む。
【００６１】
以下、図９、図１０、図１１Ａ及び図１１Ｂを参照して、本開示の実施例の方法を例示として説明する。
【００６２】
図９は、本開示の別の実施例に係る拡張データを生成するフローチャートを示す。
【００６３】
図９に示すように、この方法はステップＳ９０１〜Ｓ９１０を含む。
【００６４】
ステップＳ９０１では、心拍Ａと心拍Ｂを選択する。上記の図８の説明を参照することができるので、ここでは繰り返さない。
【００６５】
ステップＳ９０２では、未処理の一種類の誘導を選択する。
【００６６】
ステップＳ９０３では、この誘導に対し、心拍Ａ及び心拍ＢとにおけるＲ波方向が同じか否かを判断し、同じである場合、ステップＳ９０４を継続して実行し、同じではない場合、ステップＳ９０７を実行する。
【００６７】
ステップＳ９０４では、心拍Ａ及び心拍Ｂとにおける当該種類の誘導データのＲ波ピークの後の１００ｍｓ内の１番目のゼロクロス点の位置を見つける。
【００６８】
ステップＳ９０５では、この誘導に対し、Ａ及びＢの両方にゼロクロス点を見つけたか否かを判断し、見つけた場合、ステップＳ９０６を継続して実行し、見つけない場合、ステップＳ９０７を実行する。
【００６９】
実際には、ゼロクロス点を見つけることができない（例えば、ＳＴセグメント上昇）特別な波形の病変がたまにある。このような波形の場合、それと接合できる他の波形を見つけることさえできないため、接合前にこのような波形を回避することは、比較的良好な処理方式である。
【００７０】
ステップＳ９０６では、心拍Ａの当該種類の誘導データの前半部分を再サンプリングして心拍Ｃの前半部分に割り当て、心拍Ｂの当該種類の誘導データの後半部分を心拍Ｃの後半部分に割り当て、心拍Ｃの当該種類の誘導データを取得する。
【００７１】
本開示の実施例によれば、心拍Ａの当該種類の誘導データの前半部分を再サンプリングして心拍Ｃの前半部分に割り当てるのは、例えば、水平方向の圧縮又は伸張により、心拍Ａの前半部分の長さを心拍Ｂの前半部分の長さと一致させることで実行され得る。前半部分とは、心拍の開始位置から上記のゼロクロス点の位置までのデータを指し、後半部分とは、ゼロクロス点の位置から終了位置までのデータを指す。
【００７２】
図１０は、本開示の実施例に係る、心拍Ａと心拍Ｂに基づいて生成された心拍Ｃの模式図を示す。図１０に示すように、２つの心拍データの心拍Ａと心拍Ｂにより、１つの新しい心拍データの心拍Ｃを生成することができる。心拍Ｃは、心拍Ａ及び心拍Ｂの両方とは異なるものである。図１１Ａは、本開示の実施例に係る心拍Ｂと心拍Ｃの対照図を示す。図１１Ａに示すように、心拍Ｃの後半部分は、心拍Ｂと一致している。図１１Ｂは、本開示の実施例に係る心拍Ａと心拍Ｃの対照図を示す。図１１Ｂに示すように、心拍Ｃの前半部分は、心拍Ａの前半部分を再サンプリングした（伸張）結果である。図１０、図１１Ａ及び図１１Ｂは、１種類の誘導データのみを示す。本実施例の方法で処理された他の誘導のデータは、類似の特徴を有し得る。
【００７３】
図９に戻って参照されたい。ステップＳ９０７では、心拍Ｂの当該種類の誘導データを心拍Ｃに直接割り当てる。
【００７４】
ステップＳ９０８では、すべての誘導の処理が完了したか否かを判断し、完了した場合、ステップＳ９０９を実行し、完了していない場合、ステップＳ９０２に戻って、未処理の一種類の誘導を継続的に選択する。
【００７５】
ステップＳ９０９では、心拍Ｃの各誘導が心拍Ｂと同じか否かを判断し、同じである場合、心拍Ｃを棄却し、Ｓ９０１に戻って心拍を選択し直し、同じではない場合、ステップＳ９１０を実行する。
【００７６】
ステップＳ９１０では、心拍Ｃを１つの拡張データとして確定する。
【００７７】
ステップＳ９１０を実行した後、ステップＳ９０１に戻って、新しい拡張データを継続的に生成することができる。
【００７８】
本開示の実施例によって提供される技術的解決策によれば、種類毎の誘導データに対して、前記第１心拍データと第２心拍データとにおける当該種類の誘導データが、同じＲ波方向を有し且つＲ波後の予め設定された時間内にゼロクロス点がある条件を満たす場合、前記ゼロクロス点の位置に基づいて、前記第１心拍データと第２心拍データとにおける前記種類の誘導データを、仮心拍データの当該種類の誘導データとして接合するステップと、前記条件が満たされない場合は、前記第２心拍データの前記種類の誘導データを、仮心拍データの当該種類の誘導データとして確定するステップと、前記仮心拍データと第２心拍データとが完全同一でない場合、前記仮心拍データを第３心拍データとして確定するステップとを実行することによって、心電図データに対してデータ拡張を行うことができ、学習用データセットのデータ量を増加させることができ、そして、機械学習のモデルを十分に訓練して、訓練済みモデルのパフォーマンスを高めることができる。
【００７９】
本開示の実施例の方法は、既存の深層学習方法に基づいて識別精度を改善することができる。心電信号を拡張した後、学習用データセットの不十分なデータ又は不均衡なデータの問題を解決し、更に学習用データセットの不十分なデータ又は不均衡なデータにより引き起こされる識別精度が低いという問題を解決するため、このようなデータの識別精度を大幅に向上させる。全体的な有効データの量を増加させるため、不均衡なデータの問題を解決するとき、個別的な少数のデータの識別精度を向上させるとともに、すべてのデータに対する全体的な識別精度も向上させる。
【００８０】
図１２は、本開示の実施例に係る心電図データ拡張装置のブロック図を示す。この装置は、ソフトウェア、ハードウェア、又は両方の組み合わせによって、電子機器の一部又は全部になることができる。
【００８１】
図１２に示すように、前記心電図データ拡張装置１２００は、取得モジュール１２１０、処理モジュール１２２０及び生成モジュール１２３０を含む。
【００８２】
取得モジュール１２１０は、心電図データを取得するように構成されている。
処理モジュール１２２０は、前記心電図データを処理して複数の心拍データを取得するように構成されている。
生成モジュール１２３０は、前記複数の心拍データのうち少なくとも２つの心拍データに基づいて、機械学習における学習用データの拡張データを生成するように構成されている。
【００８３】
本開示の実施例によって提供された技術的解決策によれば、心電図データを取得するように構成されている取得モジュール、前記心電図データを処理して複数の心拍データを取得するように構成されている処理モジュール、及び前記複数の心拍データの少なくとも２つの心拍データに基づいて、機械学習における学習用データの拡張データを生成するように構成されている生成モジュールによって、心電図データに対してデータ拡張を行うことができ、学習用データセットのデータ量を増加させることができ、そして、機械学習のモデルを十分に訓練して、訓練済みモデルのパフォーマンスを高めることができる。
【００８４】
本開示は電子機器をさらに開示する。図１３は、本開示の実施例に係る電子機器のブロック図を示す。
【００８５】
図１３に示すように、前記電子機器１３００はメモリ１３０１とプロセッサー１３０２を含む。前記メモリ１３０１は、電子機器が上記の実施例のいずれかの心電図データの拡張方法又はコード生成方法を実行するのを支援するプログラムを記憶するために使用される。前記プロセッサー１３０２は、前記メモリ１３０１に記憶されているプログラムを実行するように構成されている。
【００８６】
本開示の実施例によれば、前記メモリ１３０１は、１つ又は複数のコンピュータ命令を記憶するために使用され、前記１つ又は複数のコンピュータ命令は、前記プロセッサー１３０２によって実行されると、：
心電図データを取得するステップと、
前記心電図データを処理して複数の心拍データを取得するステップと、
前記複数の心拍データのうち少なくとも２つの心拍データに基づいて、機械学習における学習用データの拡張データを生成するステップとが実現される。
【００８７】
本開示の実施例によれば、前記心電図データを処理して複数の心拍データを取得する前記ステップは、
前記心電図データを前処理するステップと、
前処理された心電図データに対して心拍識別を行って、複数の心拍データを取得するステップとを含む。
本開示の実施例によれば、前記前処理は、
離散ウェーブレット変換によって前記心電図データに対してマルチレベルウェーブレット分解を行って、最低レベルの近似値をゼロに設定し、離散ウェーブレット再構成によってベースライン較正後の心電図データを取得するステップを含む。
【００８８】
本開示の実施例によれば、前記前処理は、
定常ウェーブレット変換によって前記心電図データに対してマルチレベルウェーブレット分解を行って、詳細値に対して閾値フィルタリングを行い、定常ウェーブレット再構成によってノイズ低減後の心電図データを取得するステップを含む。
【００８９】
本開示の実施例によれば、前処理された心電図データに対して心拍識別を行って、複数の心拍データを取得する前記ステップは、
前処理された心電図データにおける第１誘導の誘導データに対して、定常ウェーブレット変換によってマルチレベルウェーブレット分解を行って、分解結果を取得するステップと、
前記分解結果に対して特徴識別を行って、前記第１誘導の誘導データにおけるＲ波の波ピークの位置を確定するステップと、
前記第１誘導の誘導データにおけるＲ波の波ピークの位置に基づいて、前記第１誘導及び／又は他の誘導の誘導データを複数の心拍データに分割するステップとを含む。
【００９０】
本開示の実施例によれば、前記複数の心拍データの少なくとも２つの心拍データに基づいて、拡張データを生成する前記ステップは、
前記複数の心拍データから第１心拍データ及び第２心拍データを確定するステップと、
前記第１心拍データと第２心拍データにおける少なくとも１種類の誘導データが接合可能である場合、前記第１心拍データと第２心拍データを接合して第３心拍データを生成するステップと、
前記第１心拍データ及び／又は第２心拍データのラベルを取得するステップと、
前記第１心拍データ及び／又は第２心拍データのラベルに基づいて、前記第３心拍データのラベルを確定し、ラベル付き第３心拍データを拡張データとして確定するステップとを含む。
【００９１】
本開示の実施例によれば、前記第１心拍データと第２心拍データの少なくとも１種類の誘導データが接合可能である場合、前記第１心拍データと第２心拍データを接合して第３心拍データを生成する前記ステップは、
種類毎の誘導データに対して、前記第１心拍データと第２心拍データとにおける当該種類の誘導データが、同じＲ波方向を有し且つＲ波後の予め設定された時間内にゼロクロス点がある条件を満たす場合、前記ゼロクロス点の位置に基づいて、前記第１心拍データと第２心拍データとにおける前記種類の誘導データを、仮心拍データの当該種類の誘導データとして接合するステップと、前記条件が満たされない場合は、前記第２心拍データの当該種類の誘導データを、仮心拍データの当該種類の誘導データとして確定するステップと、
前記仮心拍データと第２心拍データとが完全同一でない場合、前記仮心拍データを第３心拍データとして確定するステップとを含む。
【００９２】
図１４は、本開示の実施例に係る、心電図データの拡張方法の実現に適したコンピューターシステムの構造模式図を示す。
【００９３】
図１４に示すように、コンピューターシステム１４００は処理ユニット１４０１を含む。処理ユニット１４０１は、読取り専用メモリ（ＲＯＭ）１４０２に記憶されているプログラム、又は記憶部分１４０８からランダムアクセスメモリ（ＲＡＭ）１４０３にロードされたプログラムに従って、上記実施例における各種の処理を実行することができる。ＲＡＭ  １４０３には、システム１４００の操作に必要となる様々なプログラム及びデータが記憶されている。処理ユニット  １４０１、ＲＯＭ  １４０２及びＲＡＭ  １４０３は、バス１４０４を介して互いに接続される。入力／出力（Ｉ／Ｏ）インタフェース１４０５もバス１４０４に接続される。
【００９４】
Ｉ／Ｏインタフェース１４０５には、キーボードやマウスなどを含む入力部分１４０６と、陰極線管（ＣＲＴ）や液晶ディスプレー（ＬＣＤ）など、及びスピーカーなどを含む出力部分１４０７と、ハードディスクなどを含む記憶部分１４０８と、ＬＡＮカード、モデムなどを含むネットワークインタフェースカードの通信部分１４０９とが接続される。通信部分１４０９は、インターネットなどのネットワークを介して通信処理を実行する。ドライバ１４１０も必要に応じてＩ／Ｏインタフェース１４０５に接続される。磁気ディスク、光ディスク、光磁気ディスク、半導体メモリなどのリムーバブルメディア１４１１は、必要に応じてドライバ１４１０に取り付けられる。それにより、リムーバブルメディア１４１１から読み取られたコンピュータプログラムは、必要に応じて記憶部分１４０８にインストールされる。前記処理ユニット１４０１は、ＣＰＵ、ＧＰＵ、ＴＰＵ、ＦＰＧＡ、ＮＰＵなどの処理ユニットとして実現され得る。
【００９５】
特に、本開示の実施例によれば、上述した方法は、コンピュータソフトウェアプログラムとして実現され得る。例えば、本開示の実施例は、コンピュータプログラム製品を含む。このコンピュータプログラム製品は、読み取り可能な媒体に物理的に含まれているコンピュータプログラムを含む。前記コンピュータプログラムは、上記方法を実行するためのプログラムコードを含む。このような実施例では、このコンピュータプログラムは、通信部分１４０９を介してネットワークからダウンロードしてインストールすることができ、及び／又はリムーバブルメディア１４１１からインストールすることができる。
【００９６】
添付の図面におけるフローチャート及びブロック図は、本開示の各種の実施例のシステム、方法及びコンピュータプログラム製品による実現可能なシステムアーキテクチャ、機能及び操作を図示する。これに関して、フローチャート又はブロック図中の各ブロックは、１つのモジュール、プログラムセグメント又はコードの一部を表すことができる。前記モジュール、プログラムセグメント又はコードの一部は、所定の論理機能を実装するための１つ又は複数の実行可能命令を含む。また、一部の代替実装では、ブロックにマークされた機能が、添付の図面にマークされた順序とは異なる順序で発生する場合があることにも注意されたい。例えば、次々に表示される２つのブロックは、実際には基本的に並行して実行することも、関連する機能に応じて逆の順序で実行することもできる。また、ブロック図及び／又はフローチャート中の各ブロック、ならびにブロック図及び／又はフローチャート中のブロックの組み合わせは、所定の機能又は操作を実行するための専用のハードウェアに基づいたシステムによって実装することができ、又は専用のハードウェアとコンピュータ命令の組み合わせによって実装することができることにも注意されたい。
【００９７】
本開示の実施例で言及されているユニット又はモジュールは、ソフトウェアの形態で実現することができ、プログラマブル・ハードウェアの形態で実現することもできる。説明されているユニット又はモジュールは、プロセッサーに配置することもできる。これらのユニット又はモジュールの名称は、特定の状況で該ユニット又はモジュール自体に対する制限を構成するものではない。
【００９８】
別の態様として、本開示は、コンピュータ読み取り可能な記憶媒体をさらに提供する。このコンピュータ読み取り可能な記憶媒体は、上記実施例における電子機器又はコンピューターシステムに含まれるコンピュータ読み取り可能な記憶媒体であってもよいし、機器にインストールされていない単独で存在しているコンピュータ読み取り可能な記憶媒体であってもよい。コンピュータ読み取り可能な記憶媒体記憶には、１つ又は複数のプログラムが記憶されている。前記プログラムは、本開示に記載されている方法を実行するために、１つ又は複数のプロセッサーによって使用される。
【００９９】
上記の説明は、本開示の好ましい実施例、及び使用される技術原理の説明にすぎない。当業者であれば、本開示に係る発明の範囲は、上記の技術的特徴の特定の組み合わせによる技術的解決策に限定されるものではなく、前記発明の構想から逸脱しない前提で、上記の技術的特徴又はそれらの同等の特徴の任意の組み合わせによる他の技術的解決策も含むことを理解すべきである。例えば、上記の特徴と、本開示で開示される（これらに限定されない）類似の機能を有する技術的特徴の相互交換によって形成される技術的解決策が挙げられる。
【要約】
本開示の実施例は、心電図データの拡張方法並びにその装置、電子機器及び媒体を開示する。この心電図データ拡張方法は、心電図データを取得するステップと、前記心電図データを処理して複数の心拍データを取得するステップと、前記複数の心拍データのうち少なくとも２つの心拍データに基づいて、拡張データを生成するステップとを含む。それにより、心電図データに対してデータ拡張を行うことができ、学習用データセットのデータ量を増加させることができ、そして、機械学習のモデルを十分に訓練して、訓練済みモデルの表現を高めることができる。
【選択図】なし【図１】
0006831944.tif 000005
【図２】
0006831944.tif 000006
【図３】
0006831944.tif 000007
【図４】
0006831944.tif 000008
【図５】
0006831944.tif 000009
【図６Ａ】
0006831944.tif 000010
【図６Ｂ】
0006831944.tif 000011
【図７】
0006831944.tif 000012
【図８】
0006831944.tif 000013
【図９】
0006831944.tif 000014
【図１０】
0006831944.tif 000015
【図１１Ａ】
0006831944.tif 000016
【図１１Ｂ】
0006831944.tif 000017
【図１２】
0006831944.tif 000018
【図１３】
0006831944.tif 000019
【図１４】
0006831944.tif 000020