(57)【要約】【課題】医用画像から、組織の局所的な領域における疾患状態の変化を容易に把握できるようにすることである。【解決手段】実施形態に係る医用画像処理装置は、取得部と、解析部と、推定部とを備える。取得部は、同一被検体を撮像した画像であって、撮像時刻が異なる複数の画像を取得する。解析部は、前記複数の画像の夫々の領域の画素値に基づいて前記被検体の組織性状の解析を行うことで前記被検体の組織性状を複数の組織性状クラスに分類し、分類された前記組織性状クラスを、前記複数の画像の夫々の領域へ割り当てる。推定部は、前記複数の画像の夫々の対応領域における前記組織性状クラスの変化から、前記被検体の疾患状態の変化を推定する。【選択図】  図２
2020028706.tif 000002
【特許請求の範囲】
【請求項１】
同一被検体を撮像した画像であって、撮像時刻が異なる複数の画像を取得する取得部と、
前記複数の画像の夫々の領域の画素値に基づいて前記被検体の組織性状の解析を行うことで前記被検体の組織性状を複数の組織性状クラスに分類し、分類された前記組織性状クラスを、前記複数の画像の夫々の領域へ割り当てる解析部と、
前記複数の画像の夫々の対応領域における前記組織性状クラスの変化から、前記被検体の疾患状態の変化を推定する推定部と、
を備える医用画像処理装置。
【請求項２】
前記複数の画像のうちのいずれかにおいて、前記対応領域ごとに、前記疾患状態の変化を識別可能とする画素値を割り当てたマップ画像を生成する生成部、を更に備え、
前記推定部は、前記対応領域ごとに疾患状態の変化を推定する、
請求項１記載の医用画像処理装置。
【請求項３】
前記解析部は、前記複数の画像の夫々の領域の画素値に基づいてテクスチャ解析を行うことで、夫々の領域を複数のテクスチャパタンのいずれかに分類し、分類された前記テクスチャパタンを、前記複数の画像の夫々の領域に割り当て、
前記推定部は、前記複数の画像の夫々の対応領域における前記テクスチャパタンの変化から、前記被検体の疾患状態の変化を推定する、
請求項１又は２に記載の医用画像処理装置。
【請求項４】
前記推定部は、前記対応領域における前記テクスチャパタンの変化から、前記対応領域における前記疾患状態の変化、及び前記疾患状態の変化速度を推定する、
請求項３に記載の医用画像処理装置。
【請求項５】
前記画像は、前記被検体の肺野の画像であって、
前記テクスチャパタンは、すりガラス状、網状、結節状、蜂巣状、及び浸潤影の少なくともいずれかを含む、
請求項３又は４に記載の医用画像処理装置。
【請求項６】
前記取得部は、第１の日時に撮像した第１画像と、前記第１の日時よりも後の第２の日時に撮像した第２画像とを取得し、
前記解析部は、前記第１画像及び第２画像の夫々の領域を、前記複数のテクスチャパタンの種類に分類し、
前記推定部は、前記第１画像及び第２画像の互いに対応する局所領域における、前記テクスチャパタンの種類の変化から、前記局所領域における前記被検体の疾患状態の変化の方向が、ａ）回復方向、ｂ）増悪方向、及び、ｃ）変化なし、のいずれかであることを推定する、
請求項３乃至５のいずれか１項に記載の医用画像処理装置。
【請求項７】
前記推定部は、前記疾患状態の変化の方向を、ピクセル単位、又は、２ピクセル以上から構成されるピクセルグループ単位、で推定する、
請求項６に記載の医用画像処理装置。
【請求項８】
ディスプレイと、
前記推定部で推定される前記疾患状態の変化の方向が、ピクセル単位、又は、２ピクセル以上から構成されるピクセルグループ単位で描出された、疾患状態変化マップ、を生成する生成部と、
前記疾患状態変化マップを前記ディスプレイに表示させる表示制御部と、
をさらに備える、
請求項６に記載の医用画像処理装置。
【請求項９】
前記生成部は、前記疾患状態の変化の方向が、前記ピクセル単位、又は、前記ピクセルグループ単位ごとに、異なる色、異なる濃度、異なる数字、及び、異なる記号の少なくとも１つを含む異なる態様によって、識別可能に描出されるように、前記疾患状態変化マップを生成する、
請求項８に記載の医用画像処理装置。
【請求項１０】
前記テクスチャパタンの種類と、前記疾患状態の進行の程度とが関連付けられたルックアップテーブルを記憶する記憶部、を更に備え、
前記推定部は、前記ルックアップテーブルを参照することにより、前記被検体の疾患状態の変化の方向を推定する、
請求項６乃至９のいずれか１項に記載の医用画像処理装置。
【請求項１１】
前記推定部は、前記ルックアップテーブルを参照することにより、回復方向における回復速度、及び、増悪方向における増悪速度、の少なくとも一方を更に推定する、
請求項１０に記載の医用画像処理装置。
【請求項１２】
前記記憶部は、異なる疾患原因に夫々対応する疾患原因別のルックアップテーブルであって、疾患原因毎に、前記テクスチャパタンの種類と前記疾患状態の進行の程度とが関連付けられた複数のルックアップテーブルをさらに記憶し、
前記推定部は、疾患原因別の前記複数のルックアップテーブルの中から、ユーザによって指定された疾患原因に対応する前記疾患原因別のルックアップテーブルを参照することにより、前記被検体の疾患状態の変化の方向を推定する、
請求項１０又は１１に記載の医用画像処理装置。
【請求項１３】
前記推定部は、前記第１画像及び第２画像の互いに対応する局所領域における、前記テクスチャパタンの種類の変化から、前記被検体の疾患原因を更に推定する、
請求項１２に記載の医用画像処理装置。
【請求項１４】
前記ルックアップテーブルは、ユーザによって、その内容を作成し、編集し、或いは、変更することが可能となるように構成されている、
請求項１０乃至１３のいずれか１項に記載の医用画像処理装置。
【請求項１５】
同一被検体を撮像した画像であって、撮像時刻が異なる複数の画像を取得する取得部と、
前記複数の画像の夫々の領域の画素値に基づいて前記被検体の組織性状の解析を行うことで前記被検体の組織性状を複数の組織性状クラスに分類し、分類された前記組織性状クラスを、前記複数の画像の夫々の領域へ割り当てる解析部と、
前記複数の画像の夫々の対応領域における前記組織性状クラスの変化から、前記被検体の疾患状態の変化を推定する推定部と、
を備える医用画像処理システム。
【請求項１６】
同一被検体を撮像した画像であって、撮像時刻が異なる複数の画像を取得し、
前記複数の画像の夫々の領域の画素値に基づいて前記被検体の組織性状の解析を行うことで前記被検体の組織性状を複数の組織性状クラスに分類し、
分類された前記組織性状クラスを、前記複数の画像の夫々の領域へ割り当て、
前記複数の画像の夫々の対応領域における前記組織性状クラスの変化から、前記被検体の疾患状態の変化を推定する、
医用画像処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、医用画像処理装置、医用画像処理システム、及び、医用画像処理方法に関する。
【背景技術】
【０００２】
今日、Ｘ線ＣＴ（Computed Tomography）装置、ＭＲＩ（Magnetic Resonance Imaging）装置、或いは超音波診断装置等のモダリティで被検体を撮像して生成された画像データを解析することにより、種々の疾患の診断が行われている。
【０００３】
例えば、Ｘ線画像やＸ線ＣＴ画像を用いて、びまん性肺疾患の診断が行われている。びまん性肺疾患は、病変が左右の肺に比較的均等に広がっている肺疾患の総称であり、びまん性肺疾患には、様々の疾患が含まれる。間質性肺炎は、びまん性肺疾患の代表格であるが、この他にも、感染症や腫瘍等によるびまん性肺疾患もある。また、間質性肺炎といってもひとつの病気ではなく、最も多い特発性間質性肺炎のほか、膠原病や過敏性肺炎、じん肺・職業性肺疾患なども、広義の間質性肺炎に含まれる。このように、びまん性肺疾患は、いくつかの個別の疾患の種類に分類することができ、これらの個別の疾患種類を、以下では疾患原因と呼ぶ場合がある。
【０００４】
一方、びまん性肺疾患のＸ線ＣＴ画像は、いくつかのテクスチャパタンの種類に分類されることが知られている。また、肺のＸ線ＣＴ画像に対してテクスチャ解析を行って肺野内をいくつかのテクスチャパタンに分類し、各テクスチャパタンの肺野全体に対する体積割合や、その経時変化を表示する方法が提案されている。
【０００５】
しかしながら、これらの従来の技術では、肺等の組織の局所的な領域が、回復方向にあるのか、或いは、増悪方向にあるのかといった、組織の局所的な領域における疾患状態の変化を把握することができなかった。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】米国特許出願公開第２０１４／０１８４６０８号明細書
【非特許文献】
【０００７】
【非特許文献１】Uchiyama Y, Katsuragawa S, Abe H, et al. Quantitative computerized analysis of diffuse lung disease in high-resolution computed tomography. Med Phys. AAPM; 2003;30(9):2440−2454.
【発明の概要】
【発明が解決しようとする課題】
【０００８】
本発明が解決しようとする課題は、医用画像から、組織の局所的な領域における疾患状態の変化を容易に把握できるようにすることである。
【課題を解決するための手段】
【０００９】
実施形態に係る医用画像処理装置は、取得部と、解析部と、推定部とを備える。取得部は、同一被検体を撮像した画像であって、撮像時刻が異なる複数の画像を取得する。解析部は、前記複数の画像の夫々の領域の画素値に基づいて前記被検体の組織性状の解析を行うことで前記被検体の組織性状を複数の組織性状クラスに分類し、分類された前記組織性状クラスを、前記複数の画像の夫々の領域へ割り当てる。推定部は、前記複数の画像の夫々の対応領域における前記組織性状クラスの変化から、前記被検体の疾患状態の変化を推定する。
【図面の簡単な説明】
【００１０】
【図１】第１の実施形態に係る医用画像処理装置を含む医用画像診断システムの構成例を示す図。
【図２】第１の実施形態に係る医用画像処理装置の構成例を示すブロック図。
【図３】第１の実施形態の医用画像処理装置の動作例を示すフローチャート。
【図４】第１画像及び第２画像の一例を示す模式図。
【図５】図３のフローチャートのステップＳＴ１０２乃至ステップＳＴ１０５の処理概念を説明する図。
【図６】異なる種類のテクスチャパタンの例を示す図。
【図７】左右の肺の領域が複数のテクスチャパタンの種類に分類された例を示す図。
【図８】状態変化方向の推定用のルックアップテーブルの一例を示す図。
【図９】疾患状態変化マップの一例を示す図。
【図１０】第１の実施形態の変形例に係る医用画像処理装置の動作例を示すフローチャート。
【図１１】疾患状態の変化の方向と変化の速度を推定するためのルックアップテーブルの一例を示す図。
【図１２】疾患状態の変化の方向及び速度を描出した疾患状態変化マップの一例を示す図。
【図１３】第２の実施形態の医用画像処理装置の動作例を示すフローチャート。
【図１４】疾患原因毎の状態変化方向の推定用ルックアップテーブルの例を示す図。
【図１５】第２の実施形態の変形例に係る医用画像処理装置の動作例を示すフローチャート。
【図１６】テクスチャパタンの種類の変化から疾患原因を推定する処理の動作概念を説明する図。
【発明を実施するための形態】
【００１１】
実施形態に係る医用画像処理装置について、添付図面を参照して説明する。なお、以下の実施形態では、同一の参照符号を付した部分は同様の動作をするものとして、重複する説明を適宜省略する。
【００１２】
図１は、実施形態に係る医用画像処理装置１００を含む、医用画像処理システムの一構成例を示す図である。医用画像処理システムは、例えば、病院内において、医用画像を取得し、画像処理し、保存し、利用する、医用画像に関する一連の処理システムである。
【００１３】
医用画像処理システムは、患者等の被検体から医用画像を取得するＸ線ＣＴ装置５１１、ＭＲＩ装置５１２、超音波診断装置５１３などのモダリティ５１０（即ち、医用画像診断装置５１０）のほか、画像サーバや、医用画像処理装置１００を有している。これらの各装置は、各種のデータや医用画像を授受できるように、例えば病院内のネットワーク５００を介して互いに接続されている。
【００１４】
（第１の実施形態の構成）
図２は、第１の実施形態に係る医用画像処理装置１００の構成例を示すブロック図である。医用画像処理装置１００は、例えば、入力インターフェース回路１０、処理回路２０、記憶回路３０、入力デバイス４０、及び、ディスプレイ５０を備えて構成される。医用画像処理装置１００は、例えば、所謂ワークステーション、或いは、高性能パーソナルコンピュータとして構成される。
【００１５】
入力インターフェース回路１０は、有線又は無線のネットワーク、或いは、専用又は汎用の通信回線を介してデータを入力したり、光ディスクやＵＳＢメモリ等の記憶媒体を介してデータを入力したりするインターフェース回路である。第１の実施形態の医用画像処理装置１００では、入力インターフェース回路１０を介して、Ｘ線ＣＴ装置５１１等のモダリティ５１０で撮像された第１画像及び第２画像、或いは、画像サーバに保存されている第１画像及び第２画像を取得する。
【００１６】
ここで、第１画像及び第２画像は、同一被検体を、異なる日時に撮像した画像である。例えば、第２画像は、第１画像の撮像日よりも後の日時に撮像された画像である。第１画像及び第２画像に関しては、さらに詳しく後述する。
【００１７】
記憶回路３０は、ＲＯＭ（Read Only Memory）やＲＡＭ（Random Access Memory）の他、ＨＤＤ（Hard Disk Drive）や光ディスク装置等の外部記憶装置を含む記憶媒体である。記憶回路３０は、後述するルックアップテーブル３１を含む各種の情報やデータを記憶する他、処理回路２０が具備するプロセッサが実行する各種のプログラムを記憶する。
【００１８】
入力デバイス４０は、例えば、マウス、キーボード、トラックボール、タッチパネル等であり、各種の情報やデータを操作者が入力するための種々のデバイスを含む。
【００１９】
ディスプレイ５０は、液晶ディスプレイパネル、プラズマディスプレイパネル、有機ＥＬパネル等の表示デバイスである。
【００２０】
処理回路２０は、例えば、ＣＰＵや、専用又は汎用のプロセッサを備える回路である。プロセッサは、記憶回路３０に記憶した各種のプログラムを実行することによって、後述する各種の機能を実現する。処理回路２０は、ＦＰＧＡ（field programmable gate array）やＡＳＩＣ（application specific integrated circuit）等のハードウェアで構成してもよい。これらのハードウェアによっても後述する各種の機能を実現することができる。また、処理回路２０は、プロセッサとプログラムによるソフトウェア処理と、ハードウェア処理とを組み合わせて、各種の機能を実現することもできる。
【００２１】
また、複数の独立したプロセッサを組み合わせて処理回路２０を構成し、各プロセッサが各機能を実現してもよい。また、プロセッサが複数設けられる場合、プログラムを記憶する記憶媒体は、プロセッサごとに個別に設けられてもよいし、１つの記憶媒体が全てのプロセッサの機能に対応するプログラムを一括して記憶してもよい。
【００２２】
（第１の実施形態の動作）
第１の実施形態の処理回路２０は、図２に示す各機能、即ち、画像データ取得機能２１、テクスチャ解析機能２２、疾患状態変化推定機能２３、疾患状態変化マップ生成機能２４、及び、表示制御機能２５を実現する。これらの各機能は、例えば、処理回路２０が具備するプロセッサが、記憶回路３０に記憶される所定のプログラムを実行することによって実現される。
【００２３】
画像データ取得機能２１は、同一被検体を撮像した画像であって、撮像時刻が異なる複数の画像を取得する。例えば、第１の日時に撮像した第１画像と、第１の日時よりも後の第２の日時に撮像した第２画像とを取得する。第１画像及び第２画像は、特に限定するものではないが、例えば、肺を撮像したＸ線画像や、Ｘ線ＣＴ画像である。
【００２４】
テクスチャ解析機能２２は、画像データ取得機能２１が取得した複数の画像の夫々の領域を、複数のテクスチャパタンの種類に分類する。例えば、テクスチャ解析機能２２は、画像データ取得機能２１が取得した第１画像及び第２画像に対して、公知のテクスチャ解析を行って、第１画像及び第２画像の夫々の領域を、複数のテクスチャパタンの種類に分類する。
【００２５】
疾患状態変化推定機能２３は、テクスチャ解析機能２２によって分類された、複数の画像の互いに対応する領域間におけるテクスチャパタンの種類の変化から、被検体の疾患状態の変化を推定する。例えば、疾患状態変化推定機能２３は、第１画像及び第２画像の互いに対応する局所領域におけるテクスチャパタンの種類の変化から、局所領域における被検体の疾患状態の変化の方向が、ａ）回復方向、ｂ）増悪方向、及び、ｃ）変化なし、の少なくともいずれかであることを推定する。なお、対応する領域（即ち、対応領域）、或いは、対応する局所領域とは、例えば、複数画像間において解剖学的に対応する領域（或いは、局所領域）のことを意味している。
【００２６】
疾患状態変化マップ生成機能２４は、対応する領域ごとに、疾患状態の変化を識別可能とする画素値を割り当てたマップ画像を生成する。より具体的には、疾患状態変化マップ生成機能２４は、疾患状態変化推定機能２３によって推定された疾患状態の変化の方向が、ピクセル単位、又は、２ピクセル以上から構成されるピクセルグループ単位で描出された疾患状態変化マップを生成する。例えば、疾患状態変化マップ生成機能２４は、疾患状態の変化の方向が、ピクセル単位、又は、ピクセルグループ単位ごとに、異なる色、異なる濃度、異なる数字、及び、異なる記号の少なくとも１つを含む異なる態様によって、識別可能に描出されるように疾患状態変化マップを生成する。
【００２７】
表示制御機能２５は、生成された疾患状態変化マップを、例えば、ユーザの指示に応じて、ディスプレイ５０に表示させる。
【００２８】
図３は、第１の実施形態の医用画像処理装置１００の動作例を示すフローチャートである。また、図４乃至図９は、第１の実施形態の医用画像処理装置１００の動作概念を説明する図である。以下、図３のフローに沿って、図４乃至図９を参照しつつ、医用画像処理装置１００の動作をより具体的に説明していく。
【００２９】
まず、図３のステップＳＴ１００で第１画像を取得し、ステップＳＴ１０１で第２画像を取得する。図４は、第１画像、及び、第２画像の一例を示す模式図である。第１画像、及び、第２画像は、同一の被検体を撮像した画像であるが、互いに撮像時期が異なっている。第１画像は、例えば、2018年5月8日に行われた検査１において撮像された画像であり、第２画像は、検査１よりも後に行なわれた検査２において、例えば、2018年8月8日に行われた検査２において撮像された画像である。
【００３０】
第１画像の撮像日と、第２画像の撮像日の間隔は、上記の例では３カ月となっているが、これはあくまで一例であり、１週間、１カ月、６カ月、１２カ月、等、医師の判断によって定められるものである。また、撮像回数も２回に限定されるものではなく、３以上の複数回の撮像日に撮像された複数の画像から、撮像日の異なる所望の２つの画像を選択して、撮像日の早い一方の画像を第１画像とし、他方の画像を第２画像としてもよい。
【００３１】
なお、図４に示す第１、第２画像は、あくまでも本実施形態の医用画像処理装置１００の動作を説明するための模式図であることに留意されたい。第１、第２画像は、例えば、Ｘ線ＣＴ装置５１１で、被検体の肺野部を撮像した３次元画像を、あるコロナル断面で切り取った２次元画像である。
【００３２】
前述したように、Ｘ線画像やＸ線ＣＴ画像を用いて、例えば、びまん性肺疾患の診断が行われている。びまん性肺疾患の治療は主にステロイド薬や免疫抑制剤などの薬剤治療である。治療効果を判定するために、Ｘ線ＣＴ画像の陰影パタン、即ち、テクスチャパタンの経時変化をモニタリングすることは極めて重要である。
【００３３】
一方、びまん性肺疾患の疾患部のＸ線ＣＴ画像は、いくつかのテクスチャパタンの種類に分類されることが知られている。また、疾患が初期状態から増悪方向に進むにつれて、テクスチャパタンの種類が変化することも知られている。また逆に、治療効果が上がり、疾患が増悪状態から回復方向に向かう場合にも、テクスチャパタンの種類が変化することも知られている。
【００３４】
図４は、このようなテクスチャパタンの経時変化を、第１画像と第２画像に模式的に示したものである。第１、第２画像の夫々の肺野部の背景（周辺部を含む）の濃い領域は、疾患のない、正常な領域を示している。一方、図４では、背景の濃い領域以外に、淡い色のハッチングで示す複数の領域が図示されている。図４では、これらの、淡い色のハッチング領域が、疾患部の領域に該当するものとしている。また、図４における異なるハッチングパタンは、Ｘ線ＣＴ画像における異なる陰影パタン（テクスチャパタン）に該当するものとしている。
【００３５】
図４に示す例では、この被検体の左右の肺のうち、右肺（図４において左側の肺）上部には、第１画像と第２画像との間で淡い色の同じ種類のテクスチャパタンが観られるものの、第１画像の右肺におけるテクスチャパタンの領域よりも、第２画像の右肺におけるテクスチャパタンの領域の方が大きくなっている。このことは、第１画像を取得した2018年5月8日時点よりも、第２画像を取得した2018年8月8日時点（第１画像取得時の３カ月後）の方が、右肺における疾患部の領域が拡大しており、右肺に関しては、疾患が増悪方向にあることを意味している。
【００３６】
一方、左肺（図４において右側の肺）に関しては、第１画像において、右肺とは異なる種類のテクスチャパタンが広い範囲に亘って発生しているものの、このテクスチャパタンの種類に対応する疾患領域は、第２画像においては縮小している。その一方で、第２画像においては、左肺の下部に、上記２つの種類とは異なる種類のテクスチャパタンが発生している。このことは、左肺に関しては、回復方向にある領域と、増悪方向にある領域とが混在していることを意味している。
【００３７】
従来、このようなＸ線ＣＴ画像におけるテクスチャパタンの変化は、医師の主観的な評価に委ねられている場合が多かった。これに対して、本実施形態の医用画像処理装置１００では、以下に説明するように、テクスチャパタンの種類の局所的な変化を検出し、テクスチャパタンの種類の局所的な変化から、疾患が増悪方向に向かっているのか、或いは、回復方向に向かっているのかといった疾患状態の変化の方向を、装置が客観的に推定することを可能としている。
【００３８】
図３に戻り、ステップＳＴ１０２では、取得した第１画像のテクスチャ解析を行い、ステップＳＴ１０３では、テクスチャ解析に基づいて、第１画像の領域を、複数のテクスチャパタンの種類に分類する。同様に、ステップＳＴ１０４では、取得した第２画像のテクスチャ解析を行い、ステップＳＴ１０５では、テクスチャ解析に基づいて、第２画像の領域を、複数のテクスチャパタンの種類に分類する。
【００３９】
図５は、ステップＳＴ１０２乃至ステップＳＴ１０５の処理概念を説明する図である。図５の第２画像の中に示した四角の破線で囲んだ領域が、テクスチャ解析の判定ウィンドである。テクスチャ解析では、例えば、この判定ウィンド内の画素値（解析対象画像がＸ線ＣＴ画像の場合は、画素毎のＣＴ値でもよい）を解析することにより、当該対象領域の特徴量を算出する。例えば、画素値の平均値、標準偏差等の統計量、画素値のヒストグラム分布や、スペクトラム分布に関するパラメータ等を、特徴量として算出する。
【００４０】
判定ウィンドを、例えば、図５に例示したように、左右上下方向にシフトさせることにより、画像全体の特徴量の分布を算出することができる。シフトのステップ幅は、ピクセル単位でもよいし、複数ピクセル単位でもよい。或いは、画像全体を予め判定ウィンドの大きさの網目で分割しておき、各網目内の特徴量を順次に、或いは並行して、算出してもよい。
【００４１】
更に、ステップＳＴ１０３、ＳＴ１０５では、算出した特徴量に基づいて第１、及び第２画像の夫々の領域を複数のテクスチャパタンの種類に分類する。例えば、非特許文献１に開示されているように、びまん性肺疾患などの疾患部の医用画像（例えば、Ｘ線ＣＴ画像）を、複数種類のテクスチャパタンに分類する。
【００４２】
例えば、ステップＳＴ１０３、ＳＴ１０５では、第１画像及び第２画像を、図６に示すように、正常なテクスチャパタン（種類「Ａ」）と、５種類の異常なテクスチャパタン（種類「Ｂ」、種類「Ｃ」、種類「Ｄ」、種類「Ｅ」、及び、種類「Ｆ」）に分類する。
【００４３】
例えば、テクスチャパタンを、種類「Ａ」の「正常」パタン（normal）、種類「Ｂ」の「すりガラス状」パタン（ground-glass opacities）、種類「Ｃ」の「網状」パタン（reticular and linear opacities）、種類「Ｄ」の「結節状」パタン（nodular opacities）、種類「Ｅ」の「蜂巣状」パタン（honeycombing）、種類「Ｆ」の「浸潤影」パタン（consolidation）等の異なる種類に分類する。
【００４４】
なお、図６における夫々のテクスチャパタンは、あくまでも本実施形態の動作を説明するために例示したものであり、医学的な正しさや正確性を意図するものではない。
【００４５】
図７は、ステップＳＴ１０３及びステップＳＴ１０５の処理によって、第１画像、及び第２画像の夫々の左右の肺の領域が、複数のテクスチャパタンの種類に分類された例を示す図である。
【００４６】
第１画像に示す例では、右肺の上部が、「すりガラス状」パタンを示す種類「Ｂ」のテクスチャパタンに分類され、右肺の他の領域は、「正常」パタンを示す種類「Ａ」のテクスチャパタンに分類されている。また第１画像の左肺の大部分は、種類「Ｃ」の「網状」パタンを示すテクスチャパタンに分類されている。
【００４７】
一方、第２画像に示す例では、右肺のかなりの部分が、「すりガラス状」パタンを示す種類「Ｂ」のテクスチャパタンに分類されている。また、第２画像の左肺は、種類「Ｃ」の「網状」パタンを示すテクスチャパタンと、種類「Ｄ」の「結節状」パタンを示すテクスチャパタンとに分類されている。
【００４８】
図７に示す例では、複数のピクセルグループ（例えば、５×５ピクセルグループ）からなる正方形のメッシュによって、左右の肺の領域が、テクスチャパタンの夫々の種類に分類されている。しかしながら、前述したように、所定の大きさの判定ウィンドを、ピクセル単位で左右上下方向にシフトさせることにより、ピクセル単位で特徴量の分布を算出することも可能である。この場合、テクスチャパタンの種類の分布を、ピクセル単位の滑らかさで求めることも可能である。
【００４９】
図３に戻り、ステップＳＴ１０６では、必要に応じて第１画像と第２画像との位置合わせを行った後、第１画像と第２画像の互いに対応する局所領域における、テクスチャパタンの種類の変化から、疾患状態の変化の方向を推定する。
【００５０】
ここで、「疾患状態の変化の方向を推定する」とは、疾患状態、即ち、疾患部の病状が、ａ）回復方向、ｂ）増悪方向、及び、ｃ）変化なし、の少なくともいずれかであることを推定することである。例えば、テクスチャパタンの種類の遷移と、疾患状態の変化の方向（即ち、当該疾患部が、回復方向にあるのか、増悪方向にあるのか、或いは、変化が見られない、のいずれかの方向）とが関連付けられたルックアップテーブル３１を参照することにより、疾患状態の変化の方向を推定することができる。このルックアップテーブル３１は、例えば、記憶回路３０に保存されている。処理回路２０の疾患状態変化推定機能２３は、このルックアップテーブル３１を記憶回路３０から読み出して、ステップＳＴ１０６の処理に利用する。
【００５１】
図８の上部に、上述した状態変化方向の推定用のルックアップテーブル３１の一例を示す。また、図８の下部に、疾患状態の変化の方向と、テクスチャパタンの種類の遷移を示している。
【００５２】
図８の下部の図は、患部のテクスチャパタンの種類が、「Ａ」から「Ｂ」へ、「Ｂ」から「Ｃ」へ、「Ｃ」から「Ｄ」へと遷移するにつれて、疾患状態が、正常な状態から症状が重くなる方向へ（即ち、増悪方向へ）変化することを示している。また、逆に、図８の下部の図は、患部のテクスチャパタンの種類が、「Ｄ」から「Ｃ」へ、「Ｃ」から「Ｂ」へ、「Ｂ」から「Ａ」へと遷移するにつれて、疾患状態が、症状が重い状態から正常な状態へ（即ち、回復方向へ）変化することを示している。
【００５３】
図８の下部に示す疾患状態の変化の方向とテクスチャパタンの種類の遷移との関係は、従来からの、多数の医用画像の蓄積と、当該医用画像に基づく、医師による当該疾患の診断結果とその経過情報の蓄積等のデータベースから取得することができる。
【００５４】
図８の上部のルックアップテーブルは、図８の下部に示す疾患状態の変化の方向とテクスチャパタンの種類の遷移との関係から生成することができる。
【００５５】
例えば、図８のルックアップテーブルでは、ある被検体の検査１で撮像された画像（第１画像）のテクスチャパタンの種類「Ａ」、「Ｂ」、「Ｃ」、「Ｄ」と、同一被検体の検査１よりも後に行なわれた検査２で撮像された画像（第２画像）のテクスチャパタンの種類「Ａ」、「Ｂ」、「Ｃ」、「Ｄ」とが、疾患状態の変化を示す３つの指標（即ち、「変化なし」、「増悪」及び「回復」の３つ指標）によって、互いに関連付けられている。
【００５６】
このルックアップテーブルを参照することにより、第１画像及び第２画像の互いに対応する局所領域におけるテクスチャパタンの種類の変化から、当該局所領域における被検体の疾患状態の変化の方向が、ａ）回復方向、ｂ）増悪方向、及び、ｃ）変化なし、の少なくともいずれかであることを推定することができる。例えば、第１画像において、テクスチャパタンの種類が「Ｃ」であった局所領域が、第２画像において、テクスチャパタンの種類「Ａ」又は「Ｂ」に変化した場合は、当該局所領域は回復方向にあると推定される一方、テクスチャパタンの種類「Ｄ」に変化した場合は、当該局所領域は増悪方向にあると推定される。
【００５７】
なお、前述したように、びまん性肺疾患は複数の疾患原因を有することが知られており、疾患原因毎に、テクスチャパタンの種類の遷移と、疾患状態の変化の方向との関係が異なる可能性がある。そこで、ある特定の疾患原因（例えば、特発性間質性肺炎）によるびまん性肺疾患に対応するものであることを示すために、図８に示すルックアップテーブルには、疾患原因を特定するための付帯情報（図８の例では、疾患原因（ア））が関連付けられている。
【００５８】
また、図８のルックアップテーブルでは、「変化なし」、「増悪」及び「回復」の３つの指標の夫々に対して、３つの記号「０」、「＋」及び「−」のいずれかを併記している。これら３つの記号「０」、「＋」及び「−」は、図９に例示する「疾患状態変化マップ」との対応付けを判り易く説明するためのものであり、ルックアップテーブル自体に、必ずしもこれら３つの記号「０」、「＋」及び「−」の情報を含ませる必要はない。
【００５９】
図３に戻り、ステップＳＴ１０７では、局所領域における疾患状態の変化を描出した疾患状態変化マップを生成する。ステップＳＴ１０７の処理は、図２の疾患状態変化マップ生成機能２４に対応する処理である。
【００６０】
図９は、疾患状態変化マップの一例を示す図である。図９に示す疾患状態変化マップは、「Ａ」、「Ｂ」、「Ｃ」、「Ｄ」等のテクスチャパタンの種類に分類された後の２つの画像、即ち、図７に示す第１画像及び第２画像と、図８に示すルックアップテーブルに基づいて生成されたものである。
【００６１】
図９において、「＋」の記号で描出された局所領域は、検査１での撮像日から検査２の撮像日の期間に増悪したと推定される局所領域を示しており、「−」の記号で描出された局所領域は、逆に、検査１での撮像日から検査２の撮像日の期間に回復したと推定される局所領域を示している。一方、「０」の記号で描出された局所領域は、検査１での撮像日から検査２の撮像日の期間に、疾患状態が変化していないと推定される局所領域を示している。
【００６２】
疾患状態変化マップの解像度に関しても、図７に示したテクスチャパタンの種類分類後の第１画像及び第２画像と同様に、ピクセル単位の解像度で描出してもよいし、或いは、２ピクセル以上から構成されるピクセルグループ単位の解像度で描出してもよい。
【００６３】
また、図９に例示した疾患状態変化マップでは、疾患状態の変化の方向を、「＋」、「−」、「０」等の異なる記号で識別しているが、これに限定されるものではない。例えば、疾患状態の変化の方向を、異なる色、異なる濃度、異なる数字、及び、異なる記号の少なくとも１つを含む異なる態様によって、識別可能に描出することにより、疾患状態変化マップを生成することができる。
【００６４】
生成された疾患状態変化マップは、ステップＳＴ１０８で、本医用画像処理装置１００が具備するディスプレイ５０に表示される。或いは、生成された疾患状態変化マップを、ネットワーク５００を介してＸ線ＣＴ装置５１１等のモダリティ５１０に送信し、モダリティ５１０が有するディスプレイに表示させてもよい。
【００６５】
以上説明してきた第１の実施形態の医用画像処理装置１００によれば、肺等の組織の局所的な領域が、回復方向にあるのか、或いは、増悪方向にあるのかといった、組織の局所的な領域における疾患状態の変化を容易に把握することができる。
【００６６】
（第１の実施形態の変形例）
図１０は、第１の実施形態の変形例に係る医用画像処理装置１００の動作例を示すフローチャートである。第１の実施形態との相違点はステップＳＴ２００の処理であり、それ以外の処理は第１の実施形態と同じである。
【００６７】
ステップＳＴ２００では、第１画像と第２画像の互いに対応する局所領域におけるテクスチャパタンの種類の変化から、疾患状態の変化の方向だけでなく、疾患状態の変化の速度をさらに推定する。即ち、該当する組織の疾患が、徐々に増悪しているのか、急激に増悪しているのかといった増悪速度や、逆に、該当する組織の疾患が、徐々に回復しているのか、急激に回復しているのかといった回復速度を推定する。
【００６８】
第１の実施形態の変形例では、疾患状態の変化の速度を推定するために、状態変化方向の推定用のルックアップテーブル（図８）に変化速度に関する情報を付加したルックアップテーブル（図１１）を用いる。
【００６９】
図１１に示すルックアップテーブルでは、例えば、検査１の第１画像におけるテクスチャパタンの種類「Ａ」が、第２画像において、種類「Ａ」のままの場合は「変化なし」、種類「Ｂ」に遷移した場合は「増悪」方向でかつ増悪速度は「小」であると判定し、種類「Ｃ」に遷移した場合は「増悪」方向でかつ増悪速度は「中」であると判定し、種類「Ｄ」に遷移した場合は「増悪」方向でかつ増悪速度は「大」であると判定する。
【００７０】
また、例えば、検査１の第１画像におけるテクスチャパタンの種類「Ｄ」が、第２画像において、種類「Ｄ」のままの場合は「変化なし」、種類「Ｃ」に遷移した場合は「回復」方向でかつ回復速度は「小」であると判定し、種類「Ｂ」に遷移した場合は「回復」方向でかつ回復速度は「中」であると判定し、種類「Ａ」に遷移した場合は「回復」方向でかつ回復速度は「大」であると判定する。
【００７１】
このように、図１１に示すルックアップテーブルにおける疾患状態の変化の速度は、図１１の下部に示すテクスチャパタンの種類の互いの遷移距離に基づいて定められている。
【００７２】
なお、図１１に示すルックアップテーブルでは、増悪方向の速度を、小さい方から大きい方の順に、「＋１」、「＋２」、「＋３」の正の符号の数字で表し、回復方向の速度を、小さい方から大きい方の順に、「−１」、「−２」、「−３」の負の符号の数字で表している。
【００７３】
図９に示す疾患状態変化マップにおける「＋」、「−」、「０」といった記号を、上記の「＋１」、「＋２」、「＋３」、「−１」、「−２」、「−３」及び「０」のような符号付数字によって置き換えることによって、疾患状態の変化の方向に加えて、変化の速度も描出された疾患状態変化マップを生成することができる。
【００７４】
図１２は、記号や数字以外の識別情報を用いて、疾患状態の変化の方向及び速度を描出した疾患状態変化マップの一例を示す図である。図１２に示す疾患状態変化マップでは、疾患状態の変化の方向と速度を、ハッチングのパタンの種類とハッチングの濃淡で描出している。
【００７５】
第１の実施形態の変形例によれば、「増悪」、「回復」、「変化なし」といった疾患状態の変化の方向に加えて、増悪速度や回復速度が疾患状態変化マップに描出されるため、客観的で精度の高い診断が可能となる。
【００７６】
なお、上記の説明では、増悪速度や回復速度等の疾患状態変化の速度を、テクスチャパタンの種類の互いの遷移距離に基づいて求めたが、検査１における第１画像の撮像日と、検査２における第２画像の撮像日との間の期間の長さに基づいて、疾患状態変化の速度を推定することもできる。さらに、テクスチャパタンの種類の互いの遷移距離と、検査１における第１画像の撮像日と、検査２における第２画像の撮像日との間の期間の長さの双方に基づいて、疾患状態変化の速度を推定することもできる。
【００７７】
例えば、第１画像の撮像日と第２画像の撮像日との間隔が３カ月であったときに、検査１の第１画像におけるテクスチャパタンの種類「Ａ」が種類「Ｂ」に遷移した場合は「増悪」方向でかつ増悪速度は「小」であると判定するものの、第１画像の撮像日と第２画像の撮像日との間隔が１カ月であったときには、検査１の第１画像におけるテクスチャパタンの種類「Ａ」が種類「Ｂ」に遷移した場合であっても、増悪速度は「小」ではなく、「中」又は「大」であると判定する。
【００７８】
（第２の実施形態）
図１３は、第２の実施形態の医用画像処理装置１００の動作例を示すフローチャートである。前述したように、びまん性肺疾患は複数の疾患原因を有することが知られており、疾患原因毎に、テクスチャパタンの種類の遷移と、疾患状態の変化の方向との関係が異なる可能性がある。そこで、第２の実施形態の医用画像処理装置１００では、疾患原因毎に状態変化方向推定用のルックアップテーブルを保有する構成としている。
【００７９】
第２の実施形態では、図１３のステップＳＴ３００と、ステップＳＴ３０１が第１の実施形態と異なっている。ステップＳＴ３００では、医師等のユーザによって指定された疾患原因を取得する。例えば、マウスやキーボード等の入力デバイス４０を用いて、ユーザが、びまん性肺疾患の疾患原因として、疾患原因（ア）（例えば、特発性間質性肺炎）、疾患原因（イ）（例えば、膠原病）、疾患原因（ウ）（例えば、過敏性肺炎）、疾患原因（エ）（例えば、じん肺・職業性肺疾患）等を指定し、処理回路２０の疾患状態変化推定機能２３が、指定された疾患原因を取得する。
【００８０】
ステップＳＴ３０１では、第１画像と第２画像の互いに対応する局所領域におけるテクスチャパタンの種類の変化と、指定された疾患原因に対応する状態変化方向推定用のルックアップテーブルとから、疾患状態の変化の方向を推定する。
【００８１】
図１４は、疾患原因毎に異なる４つの状態変化方向推定用ルックアップテーブルの例を示している。図１４の左上の疾患原因（ア）に対応するルックアップテーブルから判るように、疾患原因（ア）に対応する疾患のテクスチャパタンの種類として、「Ａ」、「Ｂ」、「Ｃ」、「Ｄ」を想定している。
【００８２】
一方、図１４の左下の疾患原因（イ）に対応するルックアップテーブルから判るように、疾患原因（イ）に対応する疾患のテクスチャパタンの種類としては、「Ａ」、「Ｂ」、「Ｅ」、「Ｆ」を想定している。また、図１４の右側の各ルックアップテーブルから判るように、疾患原因（ウ）に対応する疾患のテクスチャパタンの種類としては、「Ａ」、「Ｂ」、「Ｈ」、「Ｊ」を想定し、疾患原因（エ）に対応する疾患のテクスチャパタンの種類としては、「Ａ」、「Ｂ」、「Ｋ」、「Ｆ」を想定している。
【００８３】
このように、疾患原因の種類によって、想定されるテクスチャパタンの種類も異なってくると考えられている。したがって、テクスチャパタンの種類の変化から、疾患状態の変化の方向や速度を推定する際には、疾患原因の種類を考慮するのが好ましい。
【００８４】
第２の実施形態の医用画像処理装置１００では、疾患原因別の複数のルックアップテーブルの中から、医師などのユーザによって指定された疾患原因に対応する疾患原因別のルックアップテーブルを参照することにより、被検体の疾患状態の変化の方向を推定している。この結果、高精度で、疾患状態変化の方向を推定することができ、信頼性の高い診断に資することが可能となる。
【００８５】
（第２の実施形態の変形例）
図１５は、第２の実施形態の変形例に係る医用画像処理装置１００の動作例を示すフローチャートである。第２の実施形態の変形例の動作では、第２の実施形態の処理（図１３に示すフローチャート）に対して、更に、ステップＳＴ４００を付加した処理が行われている。ステップＳＴ４００では、テクスチャパタンの種類の変化から、疾患原因を推定する。
【００８６】
図１６は、ステップＳＴ４００の処理の概念を説明する図である。図１６の中央部には、疾患原因推定用のルックアップテーブルを例示している。このルックアップテーブルは、図１４に示す４つのルックアップテーブルと実質的に同じものであり、相互の変換可能である。
【００８７】
前述したように、疾患の進行によって変化しうるテクスチャパタンの種類は、疾患原因によって異なることがあり、疾患の進行によって変化しうるテクスチャパタンの種類の全てが必ずしも同一というわけではない。例えば、図１６の疾患原因推定用のルックアップテーブルに示すように、テクスチャパタンの種類「Ａ」から「Ｂ」への変化は、４つの疾患原因（ア）乃至疾患原因（エ）において共通に生じ得るものの、テクスチャパタンの種類「Ｂ」から「Ｈ」への変化は、疾患原因（ウ）のみによって生じ得る。
【００８８】
そこで、第２の実施形態の変形例では、ステップＳＴ３０１の処理に使用されるテクスチャパタンの変化の方向と、疾患原因推定用のルックアップテーブルとから、疾患原因を推定する。例えば、図１６の左部に示すように、被検体（例えば、患者  甲）のテクスチャパタンの変化が、ある局所領域において種類「Ａ」から種類「Ｂ」へ変化し、他の局所領域において種類「Ｂ」から「種類Ｈ」に変化したとする。この場合、ステップＳＴ４００の処理は、患者  甲の疾患原因は、疾患原因（ウ）である、と推定する。
【００８９】
このように、第２の実施形態の変形例では、テクスチャパタンの種類の変化から、疾患原因を推定することが可能となるため、診断上、より有用な情報をユーザに提供することが可能となる。
【００９０】
なお、ここまで、図８、図１１、図１４、及び、図１６に、いくつかのルックアップテーブルを例示したが、これらのルックアップテーブルは、ユーザによって、その内容を作成し、編集し、或いは、変更することが可能となるように構成されている。
【００９１】
（その他の実施形態）
ここまでは、主に、複数の画像（例えば、前述した第１画像と第２画像）の夫々の領域の画素値に基づいてテクスチャ解析を行うことにより、夫々の領域を複数のテクスチャパタンのいずれかに分類する実施形態を説明してきた。しかしながら、本発明の実施形態はこれらに限定されるものではなく、例えば、ａ）複数の画像の夫々の領域の画素値に基づいて被検体の組織性状の解析を行うことでその解析結果から被検体の組織性状を複数の組織性状クラスに分類し、ｂ）分類された組織性状クラスを、複数の画像の夫々の領域へ割り当て、ｃ）複数の画像の夫々の対応領域における組織性状クラスの変化から、被検体の疾患状態の変化を推定する、実施形態を含む。
【００９２】
例えば、被検体の組織性状の解析の例として、上述したテクスチャ解析の他、ＴＩＣ（Time-Intensity Curve）解析（例えば、造影濃度の時間変化の解析）が挙げられる。この実施形態では、ＴＩＣ解析を行い、その結果から得られるＴＩＣのカーブを複数のタイプ、すなわち組織性状クラスに分類し、ＴＩＣカーブのタイプの変化から被検体の疾患状態の変化を推定する。例えば正常な肝臓内に肝細胞癌がある被検体に対して造影ＣＴ撮影を行った場合、正常な組織のＴＩＣカーブに比べ早期相（動脈相）のカーブ立ち上がりが早く、また後期相（平衡相）においてはWash outしカーブが立ち下がる。この立ち上がり・立ち下がりの程度・時期を元に組織をパタン分類することで、各組織が正常組織から癌のパタンに近づいている（増悪している）のか、あるいはその逆なのかの遷移を推定することができる。
【００９３】
また例えば、被検体の組織性状の解析の例として、画像の領域内の輝度値を解析する例が挙げられる。この実施形態では、輝度値の解析によって、例えば、輝度値のヒストグラムを算出し、ヒストグラムの形状を複数のタイプ、すなわち組織性状クラスに分類する。例えば、ヒストグラム形状のピークの鋭さを表わす尖度、ヒストグラムのピークの数、ヒストグラムのピークに対応する輝度値等のパラメータを用いてヒストグラムの形状を複数のタイプに分類し、ヒストグラムの形状のタイプの変化から被検体の疾患状態の変化を推定する。例えば肺結節を持つ被検体に対してＣＴ撮影を行った場合、肺結節となる領域では正常な領域と組織の密度が異なるため、ＣＴ撮影を行った場合の輝度値分布、すなわちヒストグラムが異なる。先の結節の例では正常組織に比べ低濃度部分の割合が減少する一方で高濃度部分の割合が増加する。このヒストグラムの形状を元に組織をパタン分類することで、各組織が正常組織か結節のパタンに近づいている（増悪している）のか否かの変化を推定することができる。
【００９４】
また例えば、被検体の組織性状の解析の例として、画像の領域内の複数種類の指標値を解析する例が挙げられる。複数種類の指標値は、例えば、領域内の輝度値と機能指標値である。また、機能指標値は、例えば、壁運動解析から取得する、或いは、心筋シンチグラフィ等の心臓核医学検査から取得する、組織の機能を示す指標値である。この実施形態では、例えば、輝度値と機能指標値の２つの値から、領域を複数のグループに分類する。例えば、輝度値が高くて機能指標値も高いグループを「グループ１」に、輝度値が高くて機能指標値が低いグループを「グループ２」に、輝度値が低く機能指標値も低いグループを「グループ３」に分類し、そして、分類したグループの変化から被検体の疾患状態の変化を推定する。
【００９５】
以上説明した少なくとも１つの実施形態によれば、医用画像から、組織の局所的な領域における疾患状態の変化を容易に把握できるようにすることができる。
【００９６】
ここまで説明してきた各実施形態や変形例では、図２に例示した各機能２１〜２５を、１つの医用画像処理装置１００で実現するものとして説明してきたが、これに限定されるものではない。例えば、複数の医用画像処理装置１００がネットワークによって連接された医用画像処理システム（図示せず）によって、図２に例示した各機能２１〜２５を分散処理するようにしてもよい。例えば、第１及び第２の医用画像処理装置１００Ａ、１００Ｂで医用画像処理システムを構成し、図２に示した各機能２１〜２５のうち、一部の機能を第１の医用画像処理装置１００Ａで処理し、残りの機能を第２の医用画像処理装置１００Ｂで処理するようにしてもよい。医用画像処理システム内の複数の医用画像処理装置１００の夫々はネットワークを介して互いに連接されているため、各医用画像処理装置１００を異なる場所に設置することができる。
【００９７】
なお、各実施形態の記載における記憶回路、画像データ取得機能、テクスチャ解析機能、疾患状態変化推定機能、疾患状態変化マップ生成機能、及び、表示制御機能は、夫々、特許請求の範囲の記載における記憶部、取得部、解析部、推定部、生成部、及び、表示制御部の一例である。
【００９８】
なお、本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれる。
【符号の説明】
【００９９】
１０  入力インターフェース回路
２０  処理回路
２１  画像データ取得機能
２２  テクスチャ解析機能
２３  疾患状態変化推定機能
２４  疾患状態変化マップ生成機能
２５  表示制御機能
３０  記憶回路
３１  ルックアップテーブル
４０  入力デバイス
５０  ディスプレイ
１００  医用画像処理装置
５００  ネットワーク
５１１  Ｘ線ＣＴ装置
５１２  ＭＲＩ装置
５１３  超音波診断装置
【図１】
2020028706.tif 000003
【図２】
2020028706.tif 000004
【図３】
2020028706.tif 000005
【図４】
2020028706.tif 000006
【図５】
2020028706.tif 000007
【図６】
2020028706.tif 000008
【図７】
2020028706.tif 000009
【図８】
2020028706.tif 000010
【図９】
2020028706.tif 000011
【図１０】
2020028706.tif 000012
【図１１】
2020028706.tif 000013
【図１２】
2020028706.tif 000014
【図１３】
2020028706.tif 000015
【図１４】
2020028706.tif 000016
【図１５】
2020028706.tif 000017
【図１６】
2020028706.tif 000018