特許 7016680 生体情報処理装置、生体情報処理方法およびプログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2022-01-28

(45)【発行日】2022-02-07

(54)【発明の名称】生体情報処理装置、生体情報処理方法およびプログラム

(51)【国際特許分類】

   A61B 5/36 20210101AFI20220131BHJP

   A61B 5/355 20210101ALI20220131BHJP

   A61B 8/08 20060101ALI20220131BHJP

   A61B 8/14 20060101ALI20220131BHJP

   G06T 7/00 20170101ALI20220131BHJP

【ＦＩ】

A61B5/36 ZDM

A61B5/355

A61B8/08

A61B8/14

G06T7/00 350C

【請求項の数】 4

(21)【出願番号】P 2017228680

(22)【出願日】2017-11-29

(65)【公開番号】P2019097639

(43)【公開日】2019-06-24

【審査請求日】2020-09-30

(73)【特許権者】

【識別番号】000230962

【氏名又は名称】日本光電工業株式会社

(74)【代理人】

【識別番号】110001416

【氏名又は名称】特許業務法人 信栄特許事務所

(72)【発明者】

【氏名】高柳 恒夫

【審査官】樋口 祐介

(56)【参考文献】

【文献】特表２０１６－５３８０４３（ＪＰ，Ａ）

【文献】特開２０１７－１９２６０７（ＪＰ，Ａ）

【文献】特開２０１４－１００４７３（ＪＰ，Ａ）

【文献】特表２０１３－５２４９８４（ＪＰ，Ａ）

【文献】特表２００８－５１５４８６（ＪＰ，Ａ）

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ａ６１Ｂ５／０２－５／３９８

Ａ６１Ｂ８／００－８／１５

Ｇ０６Ｔ７／００－７／９０

(57)【特許請求の範囲】

【請求項1】

被検者の生体情報を取得する取得部と、
前記取得部で取得された前記生体情報を分類する分類部と、
前記分類部の分類結果に応じて、予め用意された複数のアルゴリズムから用いるアルゴリズムを選択し、選択した前記アルゴリズムにより前記分類された生体情報を分析する分析部と、
を有し、
前記生体情報は、心電図波形を含み、
前記分類部は、前記心電図波形を、ＱＴ間隔の長さに応じた複数の分類基準及びＴ波の振幅の大きさに応じた複数の分類基準を組み合わせて分類する、生体情報処理装置。

【請求項2】

前記分類部は、ニューラルネットワークを用いて予め生体情報を深層学習させることに基づいて構築されている、請求項１記載の生体情報処理装置。

【請求項3】

被検者の生体情報を取得する取得工程と、
前記取得工程で取得された前記生体情報を分類する分類工程と、
前記分類工程における分類結果に応じて、予め用意された複数のアルゴリズムから用いるアルゴリズムを選択する選択工程と、
前記選択工程で選択された前記アルゴリズムにより前記分類された生体情報を分析する分析工程と、
を含み、
前記生体情報は、心電図波形を含み、
前記分類工程は、前記心電図波形を、ＱＴ間隔の長さに応じた複数の分類基準及びＴ波の振幅の大きさに応じた複数の分類基準を組み合わせて分類する、生体情報処理方法。

【請求項4】

請求項３に記載の方法をコンピュータに実行させるためのプログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本開示は、生体情報処理装置、生体情報処理方法に関する。さらに、本開示は、当該生体情報処理装置をコンピュータに実行させるためのプログラムに関する。

【背景技術】

【0002】

非特許文献１には、測定した心電図波形などの生体情報を分析する方法が開示されている。

【先行技術文献】

【非特許文献】

【0003】

【文献】Ｈｅｒｍａｎｓ ＢＪＭ， Ｖｉｎｋ ＡＳ， Ｂｅｎｎｉｓ ＦＣ，Ｆｉｌｉｐｐｉｎｉ ＬＨ， Ｍｅｉｊｂｏｒｇ ＶＭＦ， Ｗｉｌｄｅ ＡＡＭ， ｅｔ ａｌ．（２０１７） 「Ｔｈｅ ｄｅｖｅｌｏｐｍｅｎｔ ａｎｄ ｖａｌｉｄａｔｉｏｎ ｏｆ ａｎ ｅａｓｙ ｔｏ ｕｓｅ ａｕｔｏｍａｔｉｃ ＱＴ－ｉｎｔｅｒｖａｌ ａｌｇｏｒｉｔｈｍ．」ＰＬｏＳ ＯＮＥ １２（９）： ｅ０１８４３５２． （ｈｔｔｐｓ：／／ｄｏｉ．ｏｒｇ／１０．１３７１／ｊｏｕｒｎａｌ．ｐｏｎｅ．０１８４３５２）

【発明の概要】

【発明が解決しようとする課題】

【0004】

測定した心電図波形におけるＴ波の終了点の検出を行う場合、時間軸上のＱＲＳ波の後方の一定の範囲において、Ｔ波のピークが検出できるはずという予測に基づいて、Ｔ波のピークを探していた。この場合、Ｔ波のピークを探す一定の範囲は、各種のデータベースを用いた一般的な値を用いていた。例えば、非特許文献１の場合、Ｔ波のピークを探す範囲はＲ＋５０ｍｓ～Ｒ＋０．７ＲＲ間隔、となっている。

【0005】

しかし、心電図波形は、被検者の各種疾患や測定時の環境などの様々な要因により、典型的な心電図波形の振幅や時間軸上の間隔から大きく外れてしまうことがある。このため、一般的な値を用いて心電図波形におけるＴ波を探しても、Ｔ波のピークを見つけることができず、Ｔ波の終了点を正しく検出することができないことがあった。

【0006】

本開示は、生体情報の分析を精度良く行うことができる生体情報処理装置、生体情報処理方法およびプログラムを提供することを目的とする。

【課題を解決するための手段】

【0007】

本開示の生体情報処理装置は、
被検者の生体情報を取得する取得部と、
前記取得部で取得された前記生体情報を分類する分類部と、
前記分類部の分類結果に応じて、予め用意された複数のアルゴリズムから用いるアルゴリズムを選択し、選択した前記アルゴリズムにより前記分類された生体情報を分析する分析部と、
を有する。

【0008】

上記構成によれば、まず、分類部によって生体情報を分類する。その分類結果に応じて、分析部は、予め用意された複数のアルゴリズムから用いるアルゴリズムを選択し、選択したアルゴリズムにより分類された生体情報を分析する。これにより、分類された生体情報に適するアルゴリズムによる分析を行うことができ、生体情報の分析を精度良く行うことができる。

【0009】

また、本開示の生体情報処理方法は、
被検者の生体情報を取得する取得工程と、
前記取得工程で取得された前記生体情報を分類する分類工程と、
前記分類工程における分類結果に応じて、予め用意された複数のアルゴリズムから用いるアルゴリズムを選択する選択工程と、
前記選択工程で選択された前記アルゴリズムにより前記分類された生体情報を分析する分析工程と、
を含む。

【0010】

上記方法によれば、まず、生体情報を分類する。この分類結果に応じて、予め用意された複数のアルゴリズムから用いるアルゴリズムを選択する。さらに、選択されたアルゴリズムにより分類された生体情報を分析する。この方法により、分類結果に適するアルゴリズムによる生体情報の分析を、より精度良く行うことができる。

【0011】

また、本開示のプログラムは、上記方法をコンピュータに実行させるためのプログラムである。

【0012】

上記プログラムによれば、分類結果に適するアルゴリズムによる分析をコンピュータに実行させることができ、生体情報の分析をより精度良く行うことができる。

【発明の効果】

【0013】

本開示によれば、生体情報の分析を精度良く行うことができる。

【図面の簡単な説明】

【0014】

【図1】本開示の実施形態に係る生体情報処理装置を示すハードウェア構成図である。

【図2】制御部の機能ブロック図である。

【図3】生体波形の分析について模式的に示す図である。

【図4】分類テーブルの一例を説明するための図である。

【図5】複数のアルゴリズムの一例を説明するための図である。

【図6】生体波形の分析について模式的に示す図である。

【図7】本実施形態に係る生体情報処理方法の一例を説明するためのフローチャートである。

【図8】従来の生体波形の分析について模式的に示す図である。

【発明を実施するための形態】

【0015】

以下、本実施形態について、図面を参照しながら説明する。図１は、本実施形態に係る生体情報処理装置１（以下「処理装置１」という。）のハードウェア構成図を示す。

【0016】

処理装置１は、入力部２と、制御部３と、記憶部４と、出力部５と、センサインターフェース６と、ネットワークインターフェース７とを有している。入力部２、制御部３、記憶部４、出力部５、センサインターフェース６、ネットワークインターフェース７は、バス８で電気的に接続されている。

【0017】

処理装置１は、生体情報の分析を行うための装置である。生体情報は、生体に関する各種の情報であり、例えば、心電図波形や脈波などの生体波形の情報、血管の断面や心臓などの生体画像の情報である。処理装置１は、心電計などの生体情報の処理を行う専用装置であってもよく、パーソナルコンピュータ、スマートフォン、タブレットなどの汎用デバイスであってもよく、Ａｐｐｌｅ Ｗａｔｃｈ（登録商標）などのウェアラブルデバイスであってもよい。

【0018】

入力部２は、処理装置１を操作する検査者などの入力操作を受け付けると共に、当該入力操作に対応する指示信号を生成するように構成されている。入力部２は、例えば、処理装置１の筐体に取り付けられた操作ボタン、筐体に接続されたマウスやキーボード、筐体に配置されたタッチパネルである。入力部２によって生成された指示信号は、バス８を介して制御部３に出力される。

【0019】

制御部３は、入力部２から入力される指示信号に応じて、処理装置１に関する制御を行うように構成されている。制御部３は、プロセッサとメモリを有している。メモリは、コンピュータ可読命令（プログラム）を記憶するように構成されており、例えば、各種プログラムなどが格納されたＲＯＭ（Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ）やプロセッサにより実行される各種プログラムなどが格納される複数ワークエリアを有するＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）などから構成される。プロセッサは、例えばＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）、ＭＰＵ（Ｍｉｃｒｏ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）及び／又はＧＰＵ（Ｇｒａｐｈｉｃｓ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）であって、ＲＯＭに組み込まれた各種プログラムから指定されたプログラムをＲＡＭ上に展開し、ＲＡＭとの協働で各種処理を実行するように構成されている。

【0020】

プロセッサが後述する生体情報処理方法を実行するための生体情報処理プログラムをＲＡＭ上に展開し、ＲＡＭとの協働で当該生体情報処理プログラムを実行することで、制御部３は、処理装置１の各種動作を制御してもよい。制御部３及び生体情報処理プログラムの詳細については後述する。

【0021】

記憶部４は、例えば、ＨＤＤ（Ｈａｒｄ Ｄｉｓｋ Ｄｒｉｖｅ）、ＳＳＤ（Ｓｏｌｉ
ｄ Ｓｔａｔｅ Ｄｒｉｖｅ）、ＵＳＢフラッシュメモリなどの記憶装置（ストレージ）であって、プログラムや各種データを格納するように構成されている。記憶部４に格納されるプログラムは、生体情報処理プログラムや、分析に用いる複数のアルゴリズムなどである。記憶部４に格納される各種データは、生体情報に関するデータなどである。生体情報に関するデータは、時間軸上に連続して発生する複数の心電図波形（生体波形の一例）を示す波形データなどである。波形データは、不図示の心電図センサによりリアルタイムで（計測と同時に）取得されて、センサインターフェース６を介して、記憶部４に取得されてもよい。波形データは、既に計測された計測結果データ（過去の計測結果）であってもよい。既に計測された計測結果データの場合、ネットワークインターフェース７を介して、記憶部４に取得されてもよく、ＵＳＢメモリなどの外部機器を介して取得されてもよい。

【0022】

出力部５は、処理装置１の分析結果を出力する機器で構成されている。出力部５は、例えば、液晶ディスプレイ、有機ＥＬディスプレイなどの表示機器や、インクジェットプリンタ、レーザプリンタなどのプリンタ装置などである。

【0023】

センサインターフェース６は、不図示のセンサを処理装置１に接続するためのインターフェースである。
ネットワークインターフェース７は、処理装置１を通信ネットワークに接続するように
構成されている。
処理装置１とセンサインターフェース６またはネットワークインターフェース７とは、有線で接続されてもよく、無線で接続されてもよい。

【0024】

制御部３は、図２に示すように、取得部３１と、分類部３２と、分析部３３とを含んでいる。

【0025】

取得部３１は、センサインターフェース６、ネットワークインターフェース７または記憶部４から、被検者の生体情報を取得するように構成されている。図３は、生体波形の分析について模式的に示す図であり、取得部３１が取得するＱＲＳ波、Ｔ波を含む心電図波形（生体波形の一例）が示されている。
取得部３１が取得した生体情報は、分類部３２に出力される。

【0026】

分類部３２は、取得部３１で取得された生体情報を分類するように構成されている。分類部３２は、生体波形などに見られるおおよその特徴に基づいて生体情報を分類する。
図４は、分類部３２による分類に用いられる分類テーブルの一例である。図４に示すように、分類部３２は、心電図波形におけるＱＴ時間（時間軸上の時間間隔の一例）及びＴ波の振幅（高さ）の違いに基づいて、生体波形を分類結果Ｃ１～Ｃ９に分類する。分類テーブルに用いられる分類基準は、例えば、ＱＴ時間について、３５０ｍｓ以下を「ＱＴ時間が短縮（以下「ＱＴ短縮」という）」、４５０ｍｓ以上を「ＱＴ時間が延長（以下「ＱＴ延長」という）」、３５０ｍｓより長く４５０ｍｓ未満を「ＱＴ時間が正常（ＱＴ時間が通常ともいう。以下「ＱＴ正常」という）と分類してもよい。
分類部３２は、一例として、ニューラルネットワークを用いて予め生体情報を深層学習させることに基づいて構築されていてもよい。
分類部３２の分類結果は、分類部３２から分析部３３に出力される。

【0027】

分析部３３は、分類部３２の分類結果に応じて、予め用意された複数のアルゴリズムから用いるアルゴリズムを選択し、選択したアルゴリズムにより分類された生体情報を分析するように構成されている。予め用意された複数のアルゴリズムは、それぞれ、分類部３２で用いたおおよその特徴に基づく分類に合致した最適な分析を行うことができるように構築されている。

【0028】

分析に用いる複数のアルゴリズムは、一例として、図５に示す第一アルゴリズム～第九アルゴリズムである。この場合、分析部３３は、分類部３２から取得した分類結果Ｃ１～Ｃ９に応じたアルゴリズムを、予め用意された複数の第一アルゴリズム～第九アルゴリズムから選択するように構成されている。分析部３３は、選択したアルゴリズムにより、分類結果に対応する心電図波形を分析する。分類結果に対応する心電図波形は、取得部３１から分析部３３に出力されてもよく、記憶部４から分析部３３に出力されても良い。

【0029】

次に、処理装置１の処理方法について、図３～図７を用いて説明する。図７は、本実施形態の生体情報処理方法（以下、「処理方法」という。）の一例を説明するためのフローチャートである。

【0030】

図７に示すように、ＳＴＥＰ１において、処理装置１の取得部３１は、被検者から取得される生体情報を取得する（取得工程）。本例では、取得部３１が、図３または図６に示す生体波形を取得した場合について説明する。取得部３１は、取得した生体波形をデジタルデータ（波形データ）にＡ／Ｄ変換し、この波形データを分類部３２に出力する。

【0031】

ＳＴＥＰ２において、分類部３２は、ＳＴＥＰ１で取得した波形データを分類する（分類工程）。分類部３２は、図３に関する波形データを取得した場合、「ＱＴ短縮」かつ「Ｔ波の振幅が正常（正常）」（分類結果Ｃ２）と分類する（図４）。また、分類部３２は、図６に関する波形データを取得した場合、「ＱＴ延長」かつ「Ｔ波の振幅が正常（正常）」（分類結果Ｃ８）と分類する（図４）。
分類部３２は、この分類結果を、分析部３３に出力する。

【0032】

ＳＴＥＰ３において、分析部３３は、分類部３２から、ＳＴＥＰ２の分類結果を取得する。ＳＴＥＰ４において、分析部３３は、ＳＴＥＰ３で取得した分類結果に応じて、この分類結果に適するアルゴリズムを、予め用意された複数のアルゴリズム（第一アルゴリズム～第九アルゴリズム）から選択する（選択工程）。分析部３３は、取得した分類結果が分類結果Ｃ２である場合、分類結果Ｃ２に適するアルゴリズムとして第二アルゴリズムを選択する。分析部３３は、取得した分類結果がＣ８である場合、分類結果Ｃ８に適するアルゴリズムとして第八アルゴリズムを選択する。
ＳＴＥＰ５において、分析部３３は、ＳＴＥＰ４で選択したアルゴリズム（第二アルゴリズム、第八アルゴリズム）を用いて、ＳＴＥＰ２で分類された波形データの分析を実行する（分析工程）。分析部３３の分析結果は、出力部５から出力されてもよく、記憶部４に記憶されてもよい。

【0033】

以上説明したように、本実施形態の処理装置１によれば、まず、分類部３２によって生体情報を分類する（図７のＳＴＥＰ２）。その分類結果Ｃ１～Ｃ９に応じて、分析部３３は、予め用意された複数のアルゴリズム（第一アルゴリズム～第九アルゴリズム）から用いるアルゴリズムを選択する（図７のＳＴＥＰ４）。そして、選択したアルゴリズムにより分類された生体情報を分析する（図７のＳＴＥＰ５）。これにより、分類された生体情報に適するアルゴリズムによる分析を行うことができ、生体情報の分析を精度良く行うことができる。

【0034】

すなわち、従来、心電図波形の分析において、例えばＴ波のピークの検出に一般的な値を用いる場合には、標準的な波形から外れる波形についてＴ波のピークを検出することができない場合があった。図８に、従来における標準的なＴ波のピークの検索範囲を示す。図８に示す例では、Ｔ波のピークの検索範囲からＴ波が外れており、Ｔ波のピークが検出できない。この場合は、ＱＴ時間の検出も行うことができなかった。
これに対して、上記した本実施形態の処理装置１によれば、分類部３２により生体情報が分類される。その分類結果に応じて、分析部３３により、用いるアルゴリズムが選択され、選択されたアルゴリズムを用いて分析が行われる。これにより、例えばＱＴ短縮（図３）やＱＴ延長（図６）のように、従来は正確に分析できなかった生体情報について、精度良く生体情報の分析を行うことができる。

【0035】

また、ニューラルネットワークを用いて予め生体情報を深層学習させることに基づいて構築されている分類部３２により、生体情報を分類するので、検査者の主観（検査者が設定した特徴値）に依拠しない高精度の分類を行うことができる。

【0036】

また、生体情報は、生体波形を含み、分類部３２は、生体波形の時間間隔に基づく分類を行うように構成されているので、生体波形をより適切に分類することができる。

【0037】

本実施形態に係る処理装置１をソフトウェアによって実現するためには、生体情報処理プログラムが制御部３のメモリ（ＲＯＭなどの記憶部）に予め組み込まれていてもよい。または、生体情報処理プログラムは、磁気ディスク（ＨＤＤ、フロッピーディスク）、光ディスク（ＣＤ－ＲＯＭ，ＤＶＤ－ＲＯＭ、Ｂｌｕ－ｒａｙ（登録商標）ディスクなど）、光磁気ディスク（ＭＯなど）、フラッシュメモリ（ＳＤカード、ＵＳＢメモリ、ＳＳＤなど）などのコンピュータ読取可能な記憶媒体に格納されていてもよい。この場合、処理装置１に設けられたディスクドライブなどによって記憶媒体に格納された生体情報処理プログラムが読み込まれることで、当該生体情報処理プログラムが、メモリに組み込まれる。そして、メモリに組み込まれた当該プログラムがＲＡＭ上にロードされて、プロセッサがＲＡＭ上にロードされた当該プログラムを実行する。
また、生体情報処理プログラムは、通信ネットワーク上のコンピュータからネットワークインターフェース７を介してダウンロードされてもよい。この場合も同様に、ダウンロードされた当該プログラムがメモリに組み込まれる。

【0038】

本開示は、上述した実施形態や変形例に限定されず、適宜、変形、改良などが自在である。その他、上述した実施形態における各構成要素の材質、形状、形態、数、配置場所などは、本開示を達成できるものであれば任意であり、限定されない。

【0039】

上記実施形態では、生体波形をＱＴ時間とＴ波の振幅の違いに基づいて分類する例を説明したが、本開示の分類の指標はＱＴ時間、Ｔ波に限らない。本開示の処理装置１及び処理方法は、例えば、ＰＲ間隔に基づく生体波形の分類にも用いることができる。
ＰＲ間隔に基づく心電図波形（生体波形の一例）を分類する場合、分類部３２は、取得した心電図波形（Ａ／Ｄ変換された波形データ）を、ＰＲ間隔に基づいて分類し、ＰＲ間隔が短縮された場合をＣ１１、ＰＲ間隔が通常の場合をＣ１２、ＰＲ間隔が延長された場合をＣ１３と分類する。この場合において、分析部３３は、分類結果Ｃ１１～Ｃ１３に応じたアルゴリズムを、予め用意された第十一アルゴリズム～第十三アルゴリズムから選択し、選択したアルゴリズムにより心電図波形の分析を行う。第十一アルゴリズムは、例えばＷＰＷ症候群に関するアルゴリズムとすることができ、分類結果Ｃ１１の場合に選択される。第十三アルゴリズムは、例えば１度房室ブロックに関するアルゴリズムとすることができ、分類結果Ｃ１３の場合に選択される。

【0040】

また、上記処理装置１及び処理方法は、生体情報が画像の場合にも用いることができる。心室内腔面積の算出に関して、超音波エコーで取得する画像を用いる例を説明する。
分類部３２は、取得部３１から心臓を含む画像を取得し、この画像について深層学習により心臓がおおよそ画像のどの位置にあるのかを分類する。この分類において、分類部３２は、分類する画像全体を複数の領域（例えば９領域）に区分し、各区分を分類結果Ｃ２１～Ｃ２９とする。分類部３２は、取得した画像において、心臓が位置するおおよその位置を、例えば分類結果Ｃ２８と分類する。
分析部３３は、分類結果Ｃ２８に応じたアルゴリズムを、予め用意された第二十一アルゴリズム～第二十九アルゴリズムから選択し、選択したアルゴリズムにより画像の分析を行う。例えば第二十八アルゴリズムは、心臓の内側の面積を算出するアルゴリズムとすることができる。分析部３３は、分類結果Ｃ２８に応じて第二十八アルゴリズムを選択し、選択した第二十八アルゴリズムを用いて、心臓の内壁を検出して心臓の内側の面積を算出する。

【0041】

上記では、予め用意されたアルゴリズムの数と分類結果の種類数とが同一の例を説明したが、予め用意するアルゴリズムの数と分類結果の種類数とは異なっていてもよい。例えば、分類結果の種類の数が、予め用意するアルゴリズムの数より少なくてもよい。

【0042】

また、分析に用いる複数のアルゴリズムは、上記した記憶部４に記憶される構成に限らない。分析に用いる複数のアルゴリズムは、制御部３に記憶されていてもよい。分析に用いる複数のアルゴリズムは、ＵＳＢメモリなどの外部機器に記憶されていてもよい。外部機器に記憶された複数のアルゴリズムは、ネットワークインターフェース７を介して分析部３３に取得されてもよい。

【0043】

また、上記処理装置１は出力部５を含んで構成されているが、本開示の処理装置１は上記構成に限らない。出力部５は、処理装置１の外部に備えられていてもよい。出力部５が外部に備えられる場合、処理装置１は、分類に関する情報や分析に関する情報を、ネットワークインターフェース７を介して出力部５に送信してもよく、ＵＳＢメモリなどの記憶媒体を介して出力部５に取得させてもよい。

【符号の説明】

【0044】

１ 生体情報処理装置、 ２ 入力部、 ３ 制御部、 ４ 記憶部、 ５ 出力部、 ６ センサインターフェース、 ７ ネットワークインターフェース、 ８ バス、 ３１ 取得部、 ３２ 分類部、 ３３ 分析部、 Ｃ１～Ｃ９ 分類結果

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】