(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2024139509
(43)【公開日】2024-10-09
(54)【発明の名称】プログラム、画像処理方法、画像処理装置及びモデル生成方法
(51)【国際特許分類】
A61B 8/12 20060101AFI20241002BHJP
A61B 1/00 20060101ALI20241002BHJP
A61B 1/045 20060101ALI20241002BHJP
A61B 8/06 20060101ALN20241002BHJP
【FI】
A61B8/12
A61B1/00 526
A61B1/00 530
A61B1/045 614
A61B1/045 618
A61B8/06
【審査請求】未請求
【請求項の数】10
【出願形態】OL
(21)【出願番号】P 2023050481
(22)【出願日】2023-03-27
(71)【出願人】
【識別番号】000109543
【氏名又は名称】テルモ株式会社
(74)【代理人】
【識別番号】100114557
【弁理士】
【氏名又は名称】河野 英仁
(74)【代理人】
【識別番号】100078868
【弁理士】
【氏名又は名称】河野 登夫
(72)【発明者】
【氏名】坂口 雄紀
(72)【発明者】
【氏名】井上 耕一
(72)【発明者】
【氏名】松下 祐次郎
【テーマコード(参考)】
4C161
4C601
【Fターム(参考)】
4C161AA22
4C161BB08
4C161CC07
4C161FF40
4C161FF46
4C161MM10
4C161NN01
4C161QQ09
4C161RR01
4C161RR18
4C161WW02
4C161WW08
4C161WW16
4C601BB03
4C601BB13
4C601BB14
4C601BB24
4C601DD14
4C601EE07
4C601EE09
4C601EE14
4C601FE01
4C601FE04
4C601JB34
4C601JC06
4C601KK25
4C601KK31
4C601LL33
(57)【要約】
【課題】処理負荷の増加を抑制しつつ、術者にとって有用な情報の提供が可能なプログラム等を提供する。
【解決手段】コンピュータは、プログラムに従って、カテーテルを用いて管腔器官をスキャンしたスキャンデータを取得し、取得したスキャンデータのうちで、前記管腔器官の径方向の複数角度でスキャンしたスキャンデータに基づいて、前記管腔器官の縦断層画像を複数生成する。そして、コンピュータは、生成した複数の縦断層画像に基づいて、前記管腔器官の状態を検出する第1処理を実行し、第1処理の検出結果を出力する。
【選択図】
図1
【特許請求の範囲】
【請求項1】
カテーテルを用いて管腔器官をスキャンしたスキャンデータを取得し、
前記スキャンデータのうちで、前記管腔器官の径方向の複数角度でスキャンしたスキャンデータに基づいて、前記管腔器官の縦断層画像を複数生成し、
生成した複数の前記縦断層画像に基づいて、前記管腔器官の状態を検出する第1処理を実行し、
前記第1処理の検出結果を出力する
処理をコンピュータに実行させるプログラム。
【請求項2】
前記第1処理の検出結果に基づいて、前記管腔器官の横断層画像に基づいて前記管腔器官の状態を検出する第2処理を優先的に実行すべき優先領域を特定し、
特定した優先領域を撮影位置とする前記スキャンデータに基づいて、前記管腔器官の横断層画像を複数生成し、
生成した前記横断層画像に基づいて前記第2処理を実行し、
前記第2処理の検出結果を出力する
処理を前記コンピュータに実行させる請求項1に記載のプログラム。
【請求項3】
前記管腔器官の横断層画像に基づいて前記管腔器官の状態を検出する第2処理を優先的に実行すべき優先領域の指定を受け付け、
指定された優先領域を撮影位置とする前記スキャンデータに基づいて、前記管腔器官の横断層画像を複数生成し、
生成した前記横断層画像に基づいて前記第2処理を実行し、
前記第2処理の検出結果を出力する
処理を前記コンピュータに実行させる請求項1に記載のプログラム。
【請求項4】
前記優先領域以外の領域を撮影位置とする前記スキャンデータに基づいて、前記管腔器官の横断層画像を複数生成し、
前記優先領域の横断層画像に実行した前記第2処理の処理間隔よりも大きい処理間隔で、前記優先領域以外の領域の横断層画像に基づいて前記第2処理を実行する
処理を前記コンピュータに実行させる請求項2又は3に記載のプログラム。
【請求項5】
管腔器官の縦断層画像を入力した場合に前記縦断層画像中の管腔器官の状態に関する情報を出力するように学習された第1学習モデルに、生成した前記縦断層画像を入力して、前記縦断層画像中の管腔器官の状態に関する情報を取得する前記第1処理を実行する
処理を前記コンピュータに実行させる請求項1~3のいずれか1つに記載のプログラム。
【請求項6】
管腔器官の横断層画像を入力した場合に前記横断層画像中の管腔器官の状態に関する情報を出力するように学習された第2学習モデルに、生成した前記横断層画像を入力して、前記横断層画像中の管腔器官の状態に関する情報を取得する前記第2処理を実行する
処理を前記コンピュータに実行させる請求項2又は3に記載のプログラム。
【請求項7】
管腔器官の縦断層画像を入力した場合に前記縦断層画像中の管腔器官の状態に関する情報を出力するように学習された第1学習モデルに、生成した前記縦断層画像を入力して、前記縦断層画像中の管腔器官の状態に関する情報を取得し、
取得した前記縦断層画像中の管腔器官の状態に関する情報に基づいて、前記管腔器官の横断層画像中の管腔器官の状態に関する情報を推定した推定横断層画像を取得し、
前記第2学習モデルは、前記管腔器官の横断層画像と、前記縦断層画像中の管腔器官の状態に関する情報に基づいて推定された推定横断層画像とを入力した場合に、前記横断層画像中の管腔器官の状態に関する情報を出力するように学習されており、
前記第2学習モデルに、生成した前記横断層画像と、取得した前記推定横断層画像とを入力して、生成した前記横断層画像中の管腔器官の状態に関する情報を取得する
処理を前記コンピュータに実行させる請求項6に記載のプログラム。
【請求項8】
カテーテルを用いて管腔器官をスキャンしたスキャンデータを取得し、
前記スキャンデータのうちで、前記管腔器官の径方向の複数角度でスキャンしたスキャンデータに基づいて、前記管腔器官の縦断層画像を複数生成し、
生成した複数の前記縦断層画像に基づいて、前記管腔器官の状態を検出する第1処理を実行し、
前記第1処理の検出結果を出力する
処理をコンピュータが実行する画像処理方法。
【請求項9】
制御部を備える画像処理装置であって、
前記制御部は、
カテーテルを用いて管腔器官をスキャンしたスキャンデータを取得し、
前記スキャンデータのうちで、前記管腔器官の径方向の複数角度でスキャンしたスキャンデータに基づいて、前記管腔器官の縦断層画像を複数生成し、
生成した複数の前記縦断層画像に基づいて、前記管腔器官の状態を検出する第1処理を実行し、
前記第1処理の検出結果を出力する
画像処理装置。
【請求項10】
カテーテルを用いて管腔器官をスキャンしたスキャンデータを取得し、
前記スキャンデータのうちで、前記管腔器官の径方向の所定角度でスキャンしたスキャンデータに基づいて、前記管腔器官の縦断層画像を生成し、
前記スキャンデータに基づいて、前記管腔器官の長軸方向の複数位置で撮影した横断層画像を複数生成し、
管腔器官の横断層画像を入力した場合に前記横断層画像中の管腔器官の状態に関する情報を出力するように学習された第2学習モデルに、生成した複数の横断層画像を入力して、それぞれの横断層画像中の管腔器官の状態に関する情報を取得し、
それぞれの横断層画像中の管腔器官の状態に基づいて、前記所定角度でスキャンしたスキャンデータに基づく縦断層画像中の管腔器官の状態に関する情報を生成し、
生成した前記管腔器官の縦断層画像と、生成した前記縦断層画像中の管腔器官の状態に関する情報とを関連付けた訓練データを取得し、
前記訓練データを用いて、前記管腔器官の縦断層画像を入力した場合に前記縦断層画像中の管腔器官の状態に関する情報を出力する第1学習モデルを生成する
処理をコンピュータが実行するモデル生成方法。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、プログラム、画像処理方法、画像処理装置及びモデル生成方法に関する。
【背景技術】
【0002】
冠動脈における狭窄部の治療法の一つに経皮的冠状動脈インターベンション(Percutaneous Coronary Intervention:PCI)と称される方法がある。PCIは、狭窄した病変部をバルーンカテーテルで拡張し、ステントを留置して血管を再建する低侵襲治療である。PCIの術中又は術後において、画像診断用カテーテルを用いた血管内超音波診断法(IVUS:Intra Vascular Ultra Sound)又は光干渉断層診断法(OCT:Optical Coherence Tomography)により、血管内の状態を観察することができる。特許文献1には、超音波診断のためのイメージングコアを備える超音波カテーテルが開示されている。
【先行技術文献】
【特許文献】
【0003】
【発明の概要】
【発明が解決しようとする課題】
【0004】
PCIの術中に、術者は、IVUS又はOCTによって血管内の状態を観察し、ステントのサイズ及び留置位置等を決定しており、術者の意思決定を支援するために、血管内を撮影した画像に基づいて、血管の状態を検出し、検出結果を術者に提示することが行われている。例えば、血管の横断面を撮影したIVUS画像又はOCT画像に基づいて、撮影された血管の内腔及び血管壁等の領域が検出されて提示される。このような検出処理は、ステントを留置する可能性のある治療対象領域では高精度に行う必要があるが、それ以外の領域では高精度に行う必要は必ずしもない。撮影範囲内で撮影された全ての画像に対して高精度の処理を行う場合、処理負荷及び処理時間が増加し、また、高精度の処理を行うために高性能な装置を採用する場合にはコストアップが必要となる。
【0005】
一つの側面では、処理負荷の増加を抑制しつつ、術者にとって有用な情報の提供が可能なプログラム等を提供することを目的とする。
【課題を解決するための手段】
【0006】
(1)本開示は、カテーテルを用いて管腔器官をスキャンしたスキャンデータを取得し、前記スキャンデータのうちで、前記管腔器官の径方向の複数角度でスキャンしたスキャンデータに基づいて、前記管腔器官の縦断層画像を複数生成し、生成した複数の前記縦断層画像に基づいて、前記管腔器官の状態を検出する第1処理を実行し、前記第1処理の検出結果を出力する処理をコンピュータに実行させるプログラムである。
【0007】
(2)上記(1)のプログラムにおいて、前記第1処理の検出結果に基づいて、前記管腔器官の横断層画像に基づいて前記管腔器官の状態を検出する第2処理を優先的に実行すべき優先領域を特定し、特定した優先領域を撮影位置とする前記スキャンデータに基づいて、前記管腔器官の横断層画像を複数生成し、生成した前記横断層画像に基づいて前記第2処理を実行し、前記第2処理の検出結果を出力する処理を更に前記コンピュータに実行させることが好ましい。
【0008】
(3)上記(1)又は(2)のプログラムにおいて、前記管腔器官の横断層画像に基づいて前記管腔器官の状態を検出する第2処理を優先的に実行すべき優先領域の指定を受け付け、指定された優先領域を撮影位置とする前記スキャンデータに基づいて、前記管腔器官の横断層画像を複数生成し、生成した前記横断層画像に基づいて前記第2処理を実行し、前記第2処理の検出結果を出力する処理を更に前記コンピュータに実行させることが好ましい。
【0009】
(4)上記(2)又は(3)のプログラムにおいて、前記優先領域以外の領域を撮影位置とする前記スキャンデータに基づいて、前記管腔器官の横断層画像を複数生成し、前記優先領域の横断層画像に実行した前記第2処理の処理間隔よりも大きい処理間隔で、前記優先領域以外の領域の横断層画像に基づいて前記第2処理を実行する処理を更に前記コンピュータに実行させることが好ましい。
【0010】
(5)上記(1)~(4)のいずれか1つに記載のプログラムにおいて、管腔器官の縦断層画像を入力した場合に前記縦断層画像中の管腔器官の状態に関する情報を出力するように学習された第1学習モデルに、生成した前記縦断層画像を入力して、前記縦断層画像中の管腔器官の状態に関する情報を取得する前記第1処理を実行する処理を更に前記コンピュータに実行させることが好ましい。
【0011】
(6)上記(2)~(5)のいずれか1つに記載のプログラムにおいて、管腔器官の横断層画像を入力した場合に前記横断層画像中の管腔器官の状態に関する情報を出力するように学習された第2学習モデルに、生成した前記横断層画像を入力して、前記横断層画像中の管腔器官の状態に関する情報を取得する前記第2処理を実行する処理を更に前記コンピュータに実行させることが好ましい。
【0012】
(7)上記(6)のプログラムにおいて、管腔器官の縦断層画像を入力した場合に前記縦断層画像中の管腔器官の状態に関する情報を出力するように学習された第1学習モデルに、生成した前記縦断層画像を入力して、前記縦断層画像中の管腔器官の状態に関する情報を取得し、取得した前記縦断層画像中の管腔器官の状態に関する情報に基づいて、前記管腔器官の横断層画像中の管腔器官の状態に関する情報を推定した推定横断層画像を取得し、前記第2学習モデルは、前記管腔器官の横断層画像と、前記縦断層画像中の管腔器官の状態に関する情報に基づいて推定された推定横断層画像とを入力した場合に、前記横断層画像中の管腔器官の状態に関する情報を出力するように学習されており、前記第2学習モデルに、生成した前記横断層画像と、取得した前記推定横断層画像とを入力して、生成した前記横断層画像中の管腔器官の状態に関する情報を取得する処理を更に前記コンピュータに実行させることが好ましい。
【0013】
(8)また、本開示は、カテーテルを用いて管腔器官をスキャンしたスキャンデータを取得し、前記スキャンデータのうちで、前記管腔器官の径方向の複数角度でスキャンしたスキャンデータに基づいて、前記管腔器官の縦断層画像を複数生成し、生成した複数の前記縦断層画像に基づいて、前記管腔器官の状態を検出する第1処理を実行し、前記第1処理の検出結果を出力する処理をコンピュータが実行する画像処理方法である。
【0014】
(9)また、本開示は、制御部を備える画像処理装置であって、前記制御部は、カテーテルを用いて管腔器官をスキャンしたスキャンデータを取得し、前記スキャンデータのうちで、前記管腔器官の径方向の複数角度でスキャンしたスキャンデータに基づいて、前記管腔器官の縦断層画像を複数生成し、生成した複数の前記縦断層画像に基づいて、前記管腔器官の状態を検出する第1処理を実行し、前記第1処理の検出結果を出力する画像処理装置である。
【0015】
(10)また、本開示は、カテーテルを用いて管腔器官をスキャンしたスキャンデータを取得し、前記スキャンデータのうちで、前記管腔器官の径方向の所定角度でスキャンしたスキャンデータに基づいて、前記管腔器官の縦断層画像を生成し、前記スキャンデータに基づいて、前記管腔器官の長軸方向の複数位置で撮影した横断層画像を複数生成し、管腔器官の横断層画像を入力した場合に前記横断層画像中の管腔器官の状態に関する情報を出力するように学習された第2学習モデルに、生成した複数の横断層画像を入力して、それぞれの横断層画像中の管腔器官の状態に関する情報を取得し、それぞれの横断層画像中の管腔器官の状態に基づいて、前記所定角度でスキャンしたスキャンデータに基づく縦断層画像中の管腔器官の状態に関する情報を生成し、生成した前記管腔器官の縦断層画像と、生成した前記縦断層画像中の管腔器官の状態に関する情報とを関連付けた訓練データを取得し、前記訓練データを用いて、前記管腔器官の縦断層画像を入力した場合に前記縦断層画像中の管腔器官の状態に関する情報を出力する第1学習モデルを生成する処理をコンピュータが実行するモデル生成方法である。
【発明の効果】
【0016】
一つの側面では、処理負荷の増加を抑制しつつ、術者にとって有用な情報を提供することができる。
【図面の簡単な説明】
【0017】
【
図2】画像診断用カテーテルの概要を説明する説明図である。
【
図3】センサ部を挿通させた血管の断面を示す説明図である。
【
図5】画像処理装置の構成例を示すブロック図である。
【
図7】縦断層用モデルの生成処理手順の一例を示すフローチャートである。
【
図8】縦断層用モデルの生成に用いる訓練データの説明図である。
【
図9】PCIで撮影された画像の表示処理手順の一例を示すフローチャートである。
【
図12】実施形態2の画像の表示処理手順の一例を示すフローチャートである。
【
図13】横断層用モデルの変形例を示す説明図である。
【発明を実施するための形態】
【0018】
以下、本開示のプログラム、画像処理方法、画像処理装置及びモデル生成方法について、その実施形態を示す図面に基づいて詳述する。以下の各実施形態では、血管内治療である心臓カテーテル治療を一例に説明するが、カテーテル治療の対象とする管腔器官は血管に限定されず、例えば胆管、膵管、気管支、腸等の他の管腔器官であってもよい。
【0019】
(実施形態1)
図1は、画像診断装置100の構成例を示す説明図である。本実施形態では、血管内超音波診断法(IVUS)及び光干渉断層診断法(OCT)の両方の機能を備えるデュアルタイプのカテーテルを用いた画像診断装置について説明する。デュアルタイプのカテーテルでは、IVUSのみによって超音波断層画像を取得するモードと、OCTのみによって光干渉断層画像を取得するモードと、IVUS及びOCTによって両方の断層画像を取得するモードとが設けられており、これらのモードを切り替えて使用することができる。以下、超音波断層画像及び光干渉断層画像をそれぞれIVUS画像及びOCT画像と記載する。IVUS画像及びOCT画像は血管の断層画像の一例であり、IVUS画像及びOCT画像はそれぞれ、血管の径方向の断面画像である横断層画像と、血管の長軸方向の断面画像である縦断層画像とを含む。
【0020】
本実施形態の画像診断装置100は、血管内検査装置101と、血管造影装置102と、画像処理装置3と、表示装置4と、入力装置5とを備える。血管内検査装置101は、画像診断用カテーテル1及びMDU(Motor Drive Unit)2を備える。画像診断用カテーテル1は、MDU2を介して画像処理装置3に接続されている。画像処理装置3には、表示装置4及び入力装置5が接続されている。表示装置4は、例えば液晶ディスプレイ又は有機ELディスプレイ等であり、入力装置5は、例えばキーボード、マウス、タッチパネル又はマイク等である。入力装置5と画像処理装置3とは、一体に構成されていてもよい。更に入力装置5は、ジェスチャ入力又は視線入力等を受け付けるセンサであってもよい。
【0021】
血管造影装置102は画像処理装置3に接続されている。血管造影装置102は、患者の血管に造影剤を注入しながら、患者の生体外からX線を用いて血管を撮影し、当該血管の透視画像であるアンギオ画像を得るためのアンギオグラフィ装置である。血管造影装置102は、X線源及びX線センサを備え、X線源から照射されたX線をX線センサが受信することにより、患者のX線透視画像をイメージングする。なお、画像診断用カテーテル1にはX線を透過しないマーカが設けられており、アンギオ画像において画像診断用カテーテル1(マーカ)の位置が可視化される。血管造影装置102は、撮影して得られたアンギオ画像を画像処理装置3へ出力し、画像処理装置3を介して当該アンギオ画像を表示装置4に表示させる。なお、表示装置4には、アンギオ画像と、画像診断用カテーテル1を用いて撮影された断層画像とが表示される。
【0022】
図2は画像診断用カテーテル1の概要を説明する説明図である。なお、
図2中の上側の一点鎖線の領域は、下側の一点鎖線の領域を拡大したものである。画像診断用カテーテル1は、プローブ11と、プローブ11の端部に配置されたコネクタ部15とを有する。プローブ11は、コネクタ部15を介してMDU2に接続される。以下の説明では画像診断用カテーテル1のコネクタ部15から遠い側を先端側と記載し、コネクタ部15側を基端側と記載する。プローブ11は、カテーテルシース11aを備え、その先端部には、ガイドワイヤGWが挿通可能なガイドワイヤ挿通部14が設けられている。ガイドワイヤ挿通部14はガイドワイヤルーメンを構成し、予め血管内に挿入されたガイドワイヤGWを受け入れ、ガイドワイヤGWによってプローブ11を患部まで導くのに使用される。カテーテルシース11aは、ガイドワイヤ挿通部14との接続部分からコネクタ部15との接続部分に亘って連続する管部を形成している。カテーテルシース11aの内部にはシャフト13が挿通されており、シャフト13の先端側にはセンサ部12が接続されている。
【0023】
センサ部12は、ハウジング12dを有し、ハウジング12dの先端側は、カテーテルシース11aの内面との摩擦や引っ掛かりを抑制するために半球状に形成されている。ハウジング12d内には、超音波を血管内に送信すると共に血管内からの反射波を受信する超音波送受信部12a(以下ではIVUSセンサ12aという)と、近赤外光を血管内に送信すると共に血管内からの反射光を受信する光送受信部12b(以下ではOCTセンサ12bという)とが配置されている。
図2に示す例では、プローブ11の先端側にIVUSセンサ12aが設けられており、基端側にOCTセンサ12bが設けられている。画像診断用カテーテル1において、IVUSセンサ12a及びOCTセンサ12bは、シャフト13の軸方向に対して略90度となる方向(シャフト13の径方向)を超音波又は近赤外光の送受信方向として取り付けられている。なお、IVUSセンサ12a及びOCTセンサ12bは、カテーテルシース11aの内面での反射波又は反射光を受信しないように、径方向よりややずらして取り付けられることが望ましい。
【0024】
シャフト13には、IVUSセンサ12aに接続された電気信号ケーブル(図示せず)と、OCTセンサ12bに接続された光ファイバケーブル(図示せず)とが内挿されている。プローブ11は、先端側から血管内に挿入される。センサ部12及びシャフト13は、カテーテルシース11aの内部で進退可能であり、また、周方向に回転することができる。センサ部12及びシャフト13は、シャフト13の中心軸を回転軸として回転する。画像診断装置100では、センサ部12及びシャフト13によって構成されるイメージングコアを用いることにより、血管の内側から撮影されたIVUS画像、及び/又は、血管の内側から撮影されたOCT画像によって血管内部の状態を測定する。
【0025】
MDU2は、コネクタ部15によってプローブ11(画像診断用カテーテル1)が着脱可能に取り付けられる駆動装置であり、医療従事者の操作に応じて内蔵モータを駆動することにより、血管内に挿入された画像診断用カテーテル1の動作を制御する。例えばMDU2は、プローブ11に内挿されたセンサ部12及びシャフト13を一定の速度でMDU2側に向けて引っ張りながら周方向に回転させるプルバック操作を行う。センサ部12は、プルバック操作によって先端側から基端側に移動しながら回転しつつ、所定の時間間隔で連続的に血管内を走査することにより、IVUSセンサ12aが送信した超音波の血管内からの反射波を受信し、OCTセンサ12bが送信した光の血管内からの反射光を受信する。MDU2は、IVUSセンサ12aが受信した超音波の反射波データと、OCTセンサ12bが受信した反射光データとを画像処理装置3へ出力する。
【0026】
画像処理装置3は、MDU2を介してIVUSセンサ12aが受信した超音波の反射波データである信号データセットと、OCTセンサ12bが受信した反射光データである信号データセットとを取得する。画像処理装置3は、超音波の信号データセットから超音波ラインデータを生成し、生成した超音波ラインデータに基づいて血管の横断層(横断面)を撮影したIVUS横断層画像と、血管の縦断層(縦断面)を撮影したIVUS縦断層画像とを構築する。また、画像処理装置3は、反射光の信号データセットから光ラインデータを生成し、生成した光ラインデータに基づいて血管の横断層を撮影したOCT横断層画像と、血管の縦断層を撮影したOCT縦断層画像とを構築する。なお、超音波の信号データセットから超音波ラインデータを生成する処理、及び、反射光の信号データセットから光ラインデータを生成する処理は、画像処理装置3で実行されるほかに、MDU2で実行されてもよい。この場合、画像処理装置3は、MDU2から超音波ラインデータ及び光ラインデータを取得する構成となる。
【0027】
ここで、IVUSセンサ12a及びOCTセンサ12bが取得する信号データセットと、信号データセットから構築される断層画像とについて説明する。
図3は、センサ部12を挿通させた血管の断面を示す説明図であり、
図4は断層画像を説明する説明図である。まず、
図3を用いて、血管内におけるIVUSセンサ12a及びOCTセンサ12bの動作と、IVUSセンサ12a及びOCTセンサ12bによって取得される信号データセット(超音波ラインデータ及び光ラインデータ)について説明する。イメージングコアが血管内に挿通された状態で断層画像の撮影が開始されると、イメージングコアが矢符で示す方向に、シャフト13の中心軸を回転中心として回転する。このとき、IVUSセンサ12aは、各回転角度において超音波の送信及び受信を行う。ライン1,2,…512は各回転角度における超音波の送受信方向を示している。本実施形態では、IVUSセンサ12aは、血管内において360度回動(1回転)する間に512回の超音波の送信及び受信を断続的に行う。IVUSセンサ12aは、1回の超音波の送受信により、送受信方向の1ラインのデータを取得するので、1回転の間に、回転中心から放射線状に延びる512本の超音波ラインデータを得ることができる。512本の超音波ラインデータは、回転中心の近傍では密であるが、回転中心から離れるにつれて互いに疎になっていく。そこで、画像処理装置3は、各ラインの空いた空間における画素を周知の補間処理によって生成することにより、
図4A左側に示すような2次元の超音波断層画像を構築することができる。このように512本のラインデータから構築される2次元の超音波断層画像を1フレーム(1枚)のIVUS横断層画像という。なお、センサ部12は血管内を移動しながら走査するため、移動範囲内において1回転した各位置で1フレームのIVUS横断層画像が取得される。即ち、移動範囲においてプローブ11の先端側から基端側への各位置で1フレームのIVUS横断層画像が取得されるので、
図4A右側に示すように、移動範囲内で複数フレームのIVUS横断層画像が取得される。
【0028】
また、画像処理装置3は、移動範囲内で取得した超音波ラインデータのうちで、同一の回転角度で受信した超音波ラインデータを、各ラインデータの取得位置(血管の長軸方向の位置)に応じて並べることにより、
図4Bに示すような2次元の超音波断層画像を構築することができる。具体的には、画像処理装置3は、任意の回転角度で受信したラインデータ(同じライン番号のラインデータ)と、この回転角度に180度を加算した回転角度で受信したラインデータ(前記ライン番号に256を加算したライン番号のラインデータ)とから2次元の超音波断層画像を構築し、このような2次元の超音波断層画像をIVUS縦断層画像という。
【0029】
同様に、OCTセンサ12bも、各回転角度において近赤外光(測定光)の送信及び受信を行う。OCTセンサ12bも血管内において360度回動する間に512回の測定光の送信及び受信を行うので、1回転の間に、回転中心から放射線状に延びる512本の光ラインデータを得ることができる。光ラインデータについても、画像処理装置3は、各ラインの空いた空間における画素を周知の補間処理によって生成することにより、
図4Aに示すIVUS横断層画像と同様の2次元のOCT横断層画像を構築することができる。また、光ラインデータについても、画像処理装置3は、任意の回転角度で受信した光ラインデータと、この回転角度に180度を加算した回転角度で受信した光ラインデータとから、
図4Bに示すIVUS縦断層画像と同様の2次元のOCT縦断層画像を構築することができる。
【0030】
画像診断用カテーテル1は、IVUSセンサ12aによって得られるIVUS画像及び/又はOCTセンサ12bによって得られるOCT画像と、血管造影装置102によって得られるアンギオ画像との位置関係を確認するために、X線を透過しないマーカを有する。
図2に示す例では、カテーテルシース11aの先端部、例えばガイドワイヤ挿通部14にマーカ14aが設けられており、センサ部12のシャフト13側にマーカ12cが設けられている。このように構成された画像診断用カテーテル1をX線で撮影すると、マーカ14a,12cが可視化されたアンギオ画像が得られる。マーカ14a,12cを設ける位置は一例であり、マーカ12cはセンサ部12ではなくシャフト13に設けてもよく、マーカ14aはカテーテルシース11aの先端部以外の箇所に設けてもよい。また、本実施形態では、IVUSセンサ12a及びOCTセンサ12bはそれぞれ、512本のラインデータを取得する構成とするが、IVUSセンサ12a及びOCTセンサ12bが取得するラインデータの数は512本に限定されない。
【0031】
図5は画像処理装置3の構成例を示すブロック図である。画像処理装置3はコンピュータ(情報処理装置)であり、制御部31、主記憶部32、入出力部33、通信部34、補助記憶部35、読取部36を備える。制御部31は、一又は複数のCPU(Central Processing Unit)、MPU(Micro-Processing Unit)、GPU(Graphics Processing Unit)、GPGPU(General-purpose computing on graphics processing units)、TPU(Tensor Processing Unit)等の演算処理装置を用いて構成されている。制御部31は、バスを介して画像処理装置3を構成するハードウェア各部と接続されている。
【0032】
主記憶部32は、SRAM(Static Random Access Memory)、DRAM(Dynamic Random Access Memory)、フラッシュメモリ等の一時記憶領域であり、制御部31が演算処理を実行するために必要なデータを一時的に記憶する。
【0033】
入出力部33は、血管内検査装置101、血管造影装置102、表示装置4、入力装置5等の外部装置が接続されるインタフェースを備える。制御部31は、入出力部33を介して、血管内検査装置101から超音波の反射波データ及び測定光の反射光データを取得し、血管造影装置102からアンギオ画像を取得する。なお、制御部31は、血管内検査装置101から取得した反射波データから超音波ラインデータを生成し、更にIVUS画像を構築する。また制御部31は、血管内検査装置101から取得した反射光データから光ラインデータを生成し、更にOCT画像を構築する。また、制御部31は、入出力部33を介して、IVUS画像、OCT画像、又はアンギオ画像の医用画像信号を表示装置4へ出力することによって、表示装置4に医用画像を表示する。更に、制御部31は、入出力部33を介して、入力装置5に入力された情報を受け付ける。
【0034】
通信部34は、例えば、4G、5G、WiFi等の通信規格に準拠した通信インタフェースを備える。画像処理装置3は、通信部34を介して、インターネット等の外部ネットワークに接続されるクラウドサーバ等の外部サーバと通信を行う。制御部31は、通信部34を介して、外部サーバにアクセスし、当該外部サーバのストレージに記憶されている各種のデータを参照するものであってもよい。また、制御部31は、当該外部サーバと、例えばプロセス間通信を行うことにより、本実施形態における処理を協働して行うものであってもよい。
【0035】
補助記憶部35は、ハードディスク、SSD(Solid State Drive)等の記憶装置である。補助記憶部35は、制御部31が実行するプログラムP、制御部31の処理に必要な各種データを記憶する。また補助記憶部35は、例えば機械学習によって訓練データを学習済みの学習モデルが記憶されている。学習モデルは、人工知能ソフトウェアを構成するプログラムモジュールとしての利用が想定される。学習モデルは、入力値に対して所定の演算を行い、演算結果を出力するものであり、補助記憶部35には、この演算を規定する関数の係数や閾値等のデータが学習モデルとして記憶される。本実施形態では、学習モデルとして、後述する横断層用モデルM1と、縦断層用モデルM2とを記憶している。なお、補助記憶部35は画像処理装置3に接続された外部記憶装置であってもよい。プログラムP及びモデルM1,M2の一部又は全部は、画像処理装置3の製造段階において補助記憶部35に書き込まれてもよいし、遠隔のサーバ装置が配信するものを画像処理装置3が通信にて取得して補助記憶部35に記憶させてもよい。また、プログラムP及びモデルM1,M2の一部又は全部は、磁気ディスク、光ディスク、半導体メモリ等の記録媒体30に読み出し可能に記録された態様であってもよく、読取部36が記録媒体30から読み出して補助記憶部35に記憶させてもよい。
【0036】
画像処理装置3は、単一のコンピュータに限らず、複数のコンピュータを含んで構成されるマルチコンピュータであってよい。また、画像処理装置3は、サーバクライアントシステムや、クラウドサーバ、ソフトウェアによって仮想的に構築された仮想マシンであってもよい。以下の説明では、画像処理装置3が1台のコンピュータであるものとして説明する。本実施形態では、画像処理装置3に、2次元のアンギオ画像を撮影する血管造影装置102が接続されているが、生体外の複数の方向から患者の管腔器官及び画像診断用カテーテル1を撮影する装置であれば、血管造影装置102に限定されない。
【0037】
図6はモデルM1,M2の構成例を示す説明図である。
図6Aは横断層用モデルM1の構成例を示す。横断層用モデルM1(第2学習モデル)は、IVUS横断層画像又はOCT横断層画像を入力とし、入力された横断層画像中の血管内の血管内腔の領域及び血管壁の領域を認識するモデルであり、例えばセマンティックセグメンテーションにより、画像中のオブジェクト(血管内腔及び血管壁の領域)を画素単位で分類することができる。なお、血管壁は、血管内腔側から内膜、内弾性板、中膜、外弾性板、外膜で構成されている。横断層用モデルM1は、IVUS横断層画像用のモデルと、OCT横断層画像用のモデルとを含むが、以下ではまとめて横断層用モデルM1と記載する。横断層用モデルM1は、例えばU-Net、FCN、SegNet、PSPNet等の画像セグメンテーション用のアルゴリズムを用いて構成することができ、複数のアルゴリズムを組み合わせて構成してもよい。また、横断層用モデルM1は、YOLO(You Only Look Once)、SSD(Single Shot Multi-Box Detector)、ViT(Vision Transformer)等の物体検出用の学習モデルであってもよい。更に、横断層用モデルM1は、IVUS横断層画像及びOCT横断層画像の両方を入力とし、入力された横断層画像中の血管内の血管内腔の領域及び血管壁の領域を認識するモデルであってもよい。
【0038】
横断層用モデルM1は、1枚の横断層画像を入力とし、入力された横断層画像に基づいて、当該横断層画像に含まれる血管内腔の領域及び血管壁の領域を認識する演算を行い、認識結果を示す情報を出力するように学習される。具体的には、横断層用モデルM1は、入力された横断層画像の各画素を、血管内腔の領域と血管壁の領域とその他の領域とに分類し、各画素に領域毎のラベルを対応付けた分類済みの画像(以下ではラベル画像という)を出力する。
図6Aに示す横断層用モデルM1は、横断層画像中の各画素が、血管内腔の領域と、血管壁の領域と、その他の領域とに分類され、各領域に応じた画素値が割り当てられたラベル画像(血管の状態に関する情報)を出力する。
図6Aに示すラベル画像では、血管内腔の領域に分類された画素を白で示し、血管壁の領域に分類された画素をグレーで示し、その他の領域に分類された画素を黒で示している。例えば、IVUS画像においては、画像セグメンテーションで、血管内腔の輪郭線(血管内腔と血管壁との境界線)として血管内膜が認識され、血管壁の外側の輪郭線として外弾性板(EEM:External Elastic Membrane)が認識される。
【0039】
横断層用モデルM1は、訓練用の横断層画像と、この横断層画像中の各画素に対して判別すべきオブジェクト(ここでは血管内膜及び血管壁の領域)を示すデータがラベリングされた正解のラベル画像もしくはラベル配列(ラベルデータ、以下ではまとめてラベル画像と記載する)とを含む訓練データを用いて機械学習することにより生成される。正解のラベル画像では、訓練用の横断層画像に対して、各オブジェクトの領域に対応する座標範囲と、各オブジェクトの種類とを表すラベルが付与されている。正解のラベル画像は、医師等のアノテーション担当者が横断層画像中の各画素を各領域に分類し、各画素に対して、分類した領域を示す正解ラベルを割り当てることにより生成される。横断層用モデルM1は、訓練データに含まれる少なくとも1フレームの横断層画像が入力された場合に、訓練データに含まれる正解のラベル画像を出力するように学習する。具体的には、横断層用モデルM1は、入力された横断層画像に基づいて演算を行い、画像中の各オブジェクト(ここでは血管内腔及び血管壁の領域)を検出した検出結果を取得する。より具体的には、横断層用モデルM1は、横断層画像中の各画素に対して、分類されたオブジェクトの種類を示す値がラベリングされたラベル画像を出力として取得する。そして、横断層用モデルM1は、検出結果のラベル画像を、正解のラベル画像におけるオブジェクト領域の座標範囲及びオブジェクトの種類と比較し、両者が近似するように、最急降下法、誤差逆伝播法等を用いて、ニューロン間の重み等のパラメータを最適化する。これにより、横断層画像が入力された場合に、入力画像中の血管内腔の領域及び血管壁の領域を示すラベル画像を出力する横断層用モデルM1が得られる。なお、横断層用モデルM1は、横断層画像から、血管内腔及び血管壁の領域に加えて、ガイドワイヤ及びカテーテルの領域を認識する構成でもよく、血管壁中のプラークの領域を認識する構成でもよい。
【0040】
図6Bは縦断層用モデルM2の構成例を示す。縦断層用モデルM2(第1学習モデル)は、IVUS縦断層画像又はOCT縦断層画像を入力とし、入力された縦断層画像中の血管内の血管内腔の領域及び血管壁の領域を認識するモデルである。縦断層用モデルM2も、例えばセマンティックセグメンテーションにより、画像中のオブジェクト(血管内腔及び血管壁の領域)を画素単位で分類することができる。縦断層用モデルM2は、IVUS縦断層画像用のモデルと、OCT縦断層画像用のモデルとを含むが、以下ではまとめて縦断層用モデルM2と記載する。縦断層用モデルM2は、例えばU-Net、FCN、SegNet、PSPNet等の画像セグメンテーション用のアルゴリズムを用いて構成することができ、複数のアルゴリズムを組み合わせて構成してもよく、YOLO、SSD、ViT等の物体検出用の学習モデルであってもよい。また、縦断層用モデルM2は、IVUS縦断層画像及びOCT縦断層画像の両方を入力とし、入力された縦断層画像中の血管内の血管内腔の領域及び血管壁の領域を認識するモデルであってもよい。
【0041】
縦断層用モデルM2は、少なくとも1フレームの縦断層画像を入力とし、入力された縦断層画像に基づいて、当該縦断層画像に含まれる血管内腔の領域及び血管壁の領域を認識する演算を行い、認識結果を示す情報を出力するように学習される。具体的には、縦断層用モデルM2は、入力された縦断層画像の各画素を、血管内腔の領域と血管壁の領域とその他の領域とに分類し、各画素に領域毎のラベルを対応付けた分類済みの画像(以下ではラベル画像という)を出力する。
図6Bに示す縦断層用モデルM2は、縦断層画像中の各画素が、血管内腔の領域と、血管壁の領域と、その他の領域とに分類され、各領域に応じた画素値が割り当てられたラベル画像を出力する。
図6Bに示すラベル画像では、血管内腔の領域に分類された画素を白で示し、血管壁の領域に分類された画素をグレーで示し、その他の領域に分類された画素を黒で示している。
【0042】
縦断層用モデルM2は、訓練用の縦断層画像と、この縦断層画像中の各画素に対して判別すべきオブジェクト(ここでは血管内膜及び血管壁の領域)を示すデータがラベリングされた正解のラベル画像とを含む訓練データを用いて機械学習することにより生成される。正解のラベル画像では、訓練用の縦断層画像に対して、各オブジェクトの領域に対応する座標範囲と、各オブジェクトの種類とを表すラベルが付与されている。ここでも、正解のラベル画像は、医師等のアノテーション担当者が縦断層画像中の各画素を各領域に分類し、各画素に対して、分類した領域を示す正解ラベルを割り当てることにより生成される。縦断層用モデルM2は、訓練データに含まれる縦断層画像が入力された場合に、訓練データに含まれる正解のラベル画像を出力するように学習する。具体的には、縦断層用モデルM2は、入力された縦断層画像に基づいて演算を行い、画像中の各オブジェクト(ここでは血管内腔及び血管壁の領域)を検出した検出結果を取得する。より具体的には、縦断層用モデルM2は、縦断層画像中の各画素に対して、分類されたオブジェクトの種類を示す値がラベリングされたラベル画像を出力として取得する。そして、縦断層用モデルM2は、検出結果のラベル画像を、正解のラベル画像におけるオブジェクト領域の座標範囲及びオブジェクトの種類と比較し、両者が近似するように、最急降下法、誤差逆伝播法等を用いて、ニューロン間の重み等のパラメータを最適化する。これにより、縦断層画像が入力された場合に、入力画像中の血管内腔の領域及び血管壁の領域を示すラベル画像を出力する縦断層用モデルM2が得られる。なお、縦断層用モデルM2は、縦断層画像から、血管内腔及び血管壁の領域に加えて、ガイドワイヤ及びカテーテルの領域を認識する構成でもよく、血管壁中のプラークの領域を認識する構成でもよい。
【0043】
上述したモデルM1,M2の学習は、画像処理装置3で行われてもよく、他の学習装置で行われてもよい。他の学習装置で学習が行われて生成された学習済みのモデルM1,M2は、例えばネットワーク経由又は記録媒体30経由で学習装置から画像処理装置3にダウンロードされて補助記憶部35に記憶される。
【0044】
以下に、縦断層用モデルM2の生成処理について説明する。
図7は縦断層用モデルM2の生成処理手順の一例を示すフローチャート、
図8は縦断層用モデルM2の生成に用いる訓練データの説明図である。以下の処理は、画像処理装置3の制御部31が、補助記憶部35に記憶してあるプログラムPに従って行うが、他の学習装置で行われてもよい。以下の処理では、制御部31はまず、PCIの術中等に血管内検査装置101によって取得されたラインデータに基づいて訓練データを生成し、その後、生成した訓練データを用いて縦断層用モデルM2の学習を行う。なお、IVUS用の縦断層用モデルM2及びOCT用の縦断層用モデルM2は共に以下の処理によって生成可能であるので以下ではまとめて説明する。よって、以下の説明において、IVUS用の縦断層用モデルM2の生成処理では、ラインデータは超音波ラインデータ、横断層画像はIVUS横断層画像、横断層用モデルM1はIVUS用の横断層用モデルM1、縦断層画像はIVUS縦断層画像とそれぞれ読み替える。また、OCT用の縦断層用モデルM2の生成処理では、ラインデータは光ラインデータ、横断層画像はOCT横断層画像、横断層用モデルM1はOCT用の横断層用モデルM1、縦断層画像はOCT縦断層画像とそれぞれ読み替える。
【0045】
画像処理装置3の制御部31は、血管内検査装置101によって撮影されて、例えば補助記憶部35に記憶してあるラインデータを取得する(S11)。例えば制御部31は、血管内検査装置101がプルバック操作を行いつつ撮影した一連のラインデータ(血管内をスキャンしたスキャンデータ)を取得する。制御部31は、ラインデータに基づいて横断層画像を構築する(S12)。具体的には、制御部31は、センサ部12が1回転する間に取得された512本のラインデータ毎に、512本のラインデータに対して補間処理によって画素を補間して2次元の横断層画像を複数構築する。これにより、制御部31は、血管の長軸方向における撮影範囲内の複数の撮影位置で撮影した横断層画像を生成する。なお、ラインデータから構築された横断層画像が補助記憶部35に記憶されている場合、制御部31は、ステップS11~S12の代わりに、補助記憶部35から横断層画像を取得する処理を行う。
【0046】
制御部31は、取得した横断層画像のそれぞれに対してセグメンテーションを実行する(S13)。具体的には、制御部31は、それぞれの横断層画像を横断層用モデルM1に入力し、横断層用モデルM1から出力されたラベル画像を取得する。これにより、
図8上側に示すように、一連の横断層画像から、各横断層画像中の画素が血管内腔の領域、血管壁の領域、又は他の領域に分類されたラベル画像が生成される。よって、制御部31は、それぞれの横断層画像中の血管の状態(血管内腔及び血管壁の領域)に関する情報を取得できる。
【0047】
制御部31は、処理対象のライン番号に1を設定する(S14)。そして、制御部31は、横断層画像のラベル画像から、設定したライン番号(ここでは1)の縦断層画像のラベル画像を生成する(S15)。なお、横断層画像のラベル画像において、画像中心から放射線状に延びる各ライン上の画素は、横断層画像における各ライン上の画素に対応している。よって、制御部31は、横断層画像の各ラベル画像から、処理対象のライン番号、及びこのライン番号に256を加算したライン番号のライン上の画素を抽出する。例えば制御部31は、
図8上側に示す横断層画像のラベル画像中に実線で示すライン上の画素を抽出する。そして、制御部31は、横断層画像の各ラベル画像から抽出したライン上の画素を、横断層画像の撮影位置の順に並べることにより、
図8下側に示すような縦断層画像のラベル画像を生成する。これにより、制御部31は、一連の横断層画像のラベル画像(血管の状態に関する情報)に基づいて、所定角度(ここでは、ライン番号1,257の撮影方向)でスキャンしたスキャンデータに基づくラベル画像(血管の状態に関する情報)を生成する。
【0048】
次に制御部31は、ステップS11で取得したラインデータから、ここでのライン番号の縦断層画像を構築する(S16)。具体的には、制御部31は、ここでのライン番号のラインデータと、このライン番号に256を加算したライン番号のラインデータとを抽出して、撮影位置の順に並べて2次元の縦断層画像を構築する。これにより、血管の径方向の所定角度(ここでは、ライン番号1,257の撮影方向)でスキャンしたスキャンデータに基づいて、血管の縦断層画像が生成される。制御部31は、ステップS15で生成した縦断層画像のラベル画像(血管の状態に関する情報)を正解のラベル画像として、ステップS16で生成した縦断層画像と関連付けて訓練データを生成し、主記憶部32又は補助記憶部35に記憶する(S17)。制御部31は、生成した訓練データを、主記憶部32又は補助記憶部35に用意された訓練DB(図示せず)に記憶する。
【0049】
制御部31は、処理対象のライン番号に所定数を加算する(S18)。例えば、ライン番号1,2,3…のラインデータに基づいて訓練データを生成する場合、制御部31は、ライン番号に1を加算する。また、ライン番号1,3,5,7…のラインデータに基づいて訓練データを生成する場合、制御部31は、ライン番号に2を加算する。そして、制御部31は、加算後のライン番号が256よりも大きいか否かを判断し(S19)、256以下であると判断した場合(S19:NO)、ステップS15に戻り、加算後のライン番号について、ステップS15~S17の処理を実行する。これにより、所定間隔のライン番号毎に各ラインデータに基づく訓練データが生成されて記憶される。
【0050】
加算後のライン番号が256よりも大きいと判断した場合(S19:YES)、制御部31は、ステップS11で取得したラインデータに基づく訓練データの生成処理を終了したと判断し、訓練データの処理に用いる他のラインデータ(処理対象のラインデータ)があるか否かを判断する(S20)。処理対象のラインデータがあると判断した場合(S20:YES)、制御部31は、ステップS11に戻り、処理対象のラインデータについて、ステップS11~S19の処理を実行する。制御部31は、訓練データの処理に用いるラインデータがないと判断するまでステップS11~S20の処理を繰り返す。これにより、プルバック操作を行いつつ撮影した一連のラインデータを複数用いて訓練データが生成されて記憶される。なお、訓練データの生成に用いるラインデータは、補助記憶部35に記憶してある構成に限定されず、制御部31は、例えばネットワーク経由で他の装置からラインデータを取得する構成でもよい。
【0051】
処理対象のラインデータがないと判断した場合(S20:NO)、制御部31は、上述したように生成して訓練DBに蓄積した訓練データを用いて、縦断層用モデルM2の学習を行う。制御部31は、上述した処理によって訓練DBに蓄積した訓練データのうちの1つを読み出す(S21)。そして、制御部31は、読み出した訓練データに基づいて、縦断層用モデルM2の学習処理を行う(S22)。ここでは、制御部31は、訓練データに含まれる縦断層画像を縦断層用モデルM2に入力し、当該縦断層画像が入力されることによって縦断層用モデルM2から出力される出力値(縦断層画像のラベル画像)を取得する。制御部31は、縦断層用モデルM2から出力された出力値と、訓練データに含まれる正解のラベル画像とを比較し、両者が近似するように縦断層用モデルM2を学習させる。なお、学習処理において、縦断層用モデルM2は、例えば、出力層側から入力層側に向かって順次更新する誤差逆伝播法を用いて、演算処理に用いるパラメータを最適化する。
【0052】
制御部31は、訓練DBに記憶してある訓練データのうちで、学習処理が行われていない未処理の訓練データがあるか否かを判断する(S23)。未処理の訓練データがあると判断した場合(S23:YES)、制御部31は、ステップS21の処理に戻り、学習処理が未処理の訓練データに基づいてステップS21~S22の処理を行う。未処理の訓練データがないと判断した場合(S23:NO)、制御部31は、一連の処理を終了する。
【0053】
上述した学習処理により、ラインデータから構築された縦断層画像が入力された場合に、入力画像中の各画素が血管内腔の領域、血管壁の領域、及びその他の領域に分類されたラベル画像を出力する縦断層用モデルM2が生成される。なお、上述した処理において、ステップS11~S20による訓練データの生成処理と、ステップS21~S23による縦断層用モデルM2の生成処理とは、各別の装置で行われてもよい。縦断層用モデルM2は、上述したような訓練データを用いた学習処理を繰り返し行うことにより更に最適化することが可能である。また、既に学習済みの縦断層用モデルM2についても、上述した学習処理で再学習させることにより、判別精度が更に向上した縦断層用モデルM2を生成できる。
【0054】
上述した処理では、ラインデータから構築された縦断層画像と、横断層画像のラベル画像に基づいて生成された縦断層画像のラベル画像とを関連付けた訓練データを、ライン番号が1,2,3…の順に生成したが、この構成に限定されない。例えば、訓練データの生成処理を開始するライン番号は1に限定されず、任意のライン番号から処理を行ってもよい。また、訓練データの生成に用いるラインデータは、ライン番号1,2,3…の全てのデータを用いる必要はなく、例えば所定数間隔のライン番号のデータを用いて訓練データを生成してもよい。
【0055】
以下に、横断層用モデルM1及び縦断層用モデルM2を用いて画像処理装置3が行う処理について説明する。
図9は、PCIで撮影された画像の表示処理手順の一例を示すフローチャート、
図10は画像の表示処理を示す説明図、
図11は画面例を示す説明図である。本実施形態の画像処理装置3において、制御部31は、補助記憶部35に記憶されたプログラムPを読み出して実行することにより以下の処理を行う。なお、OCT画像に対して行う処理と、IVUS画像に対して行う処理とは同様の処理であり、以下ではまとめて説明する。よって、以下の説明において、OCT画像に基づく処理では、ラインデータは光ラインデータ、縦断層画像はOCT縦断層画像、横断層画像はOCT横断層画像とそれぞれ読み替える。また、IVUS画像に基づく処理では、ラインデータは超音波ラインデータ、縦断層画像はIVUS縦断層画像、横断層画像はIVUS横断層画像とそれぞれ読み替える。
【0056】
PCIの術中に、術者は、バルーンカテーテルによって血管を拡張する前、血管を拡張した後(ステントの留置前)、ステントを留置した後、留置したステントをバルーンカテーテルで圧接(後拡張)した後等、適宜のタイミングでIVUSセンサ12a及び/又はOCTセンサ12bによる撮影を行い、治療対象領域を血管内から撮影した断層画像(IVUS画像及び/又はOCT画像)を取得する。そして、術者は、取得した断層画像によって治療対象領域を観察する。以下の処理は、上述したタイミングのいずれで実行されてもよい。
【0057】
画像処理装置3の制御部31は、血管内検査装置101による血管内の撮影処理を行い、血管内検査装置101がプルバック操作を行いつつ撮影した一連のラインデータ(血管内をスキャンしたスキャンデータ)を取得する(S31)。ここでは、血管内検査装置101は、画像診断用カテーテル1のセンサ部12を先端側(遠位側)から基端側(近位側)に移動させながら血管内を走査して、IVUSセンサ12aによる撮影で得られた超音波ラインデータと、OCTセンサ12bによる撮影で得られた光ラインデータとを取得する。画像処理装置3の制御部31は、入出力部33を通じて、血管内検査装置101が取得した一連のラインデータを取得する。なお、制御部31は、MDU2を介してIVUSセンサ12aから超音波の反射波データを取得する場合、取得した超音波の反射波データから超音波ラインデータを生成し、MDU2を介してOCTセンサ12bから反射光データを取得する場合、取得した反射光データから光ラインデータを生成する。
【0058】
制御部31は、ステップS31で取得したラインデータに基づいて、所定枚数の縦断層画像を構築する(S32)。例えば制御部31は、4枚の縦断層画像を生成する場合、ライン番号1及び257のラインデータに基づく縦断層画像と、ライン番号65及び321のラインデータに基づく縦断層画像と、ライン番号129及び385のラインデータに基づく縦断層画像と、ライン番号193及び449のラインデータに基づく縦断層画像とを構築する。具体的には、制御部31は、取得したラインデータのうちで、上述した各ライン番号のラインデータを抽出し、各ラインデータを、長軸方向における撮影位置の順に並べて2次元の縦断層画像を生成する。これにより、血管の径方向の複数角度(ここでは、64毎のライン番号の撮影方向)でスキャンしたスキャンデータに基づいて4枚の縦断層画像が生成される。制御部31は、構築した縦断層画像を主記憶部32又は補助記憶部35に記憶する。なお、ここで構築する縦断層画像の枚数は4枚に限定されず、2枚、6枚、8枚等であってもよい。
【0059】
制御部31は、生成した縦断層画像のそれぞれに対してセグメンテーションを実行する(S33)。具体的には、制御部31は、それぞれの縦断層画像を縦断層用モデルM2に入力し、縦断層用モデルM2から出力されたラベル画像を取得する。これにより、制御部31は、それぞれの縦断層画像中の画素が血管内腔の領域、血管壁の領域、又は他の領域に分類されたラベル画像を取得する。制御部31は、それぞれの縦断層画像から取得したラベル画像から、血管の長軸方向の各撮影位置における血管内腔及び血管壁の領域を特定し、特定した血管内腔及び血管壁の領域に基づいて仮想横断層画像を生成する(S34)。具体的には、制御部31は、
図10A左側に示すように、4枚の縦断層画像において、長軸方向のある撮影位置における血管内腔の輪郭線(具体的には血管内膜)及び血管壁の外側の輪郭線(具体的にはEEM)を検出する。
図10Aでは血管内膜を白丸で示し、EEMを黒丸で示す。そして、制御部31は、4枚の縦断層画像から検出した血管内膜上の8点に対して、スプライン補間等によって血管内膜を補間し、4枚の縦断層画像から検出したEEM上の8点に対して、スプライン補間等によってEEMを補間する。これにより、
図10A右側に示すように、補間された血管内腔の輪郭線と血管壁の輪郭線とを含む仮想横断層画像(縦断層画像のラベル画像から推定された推定横断層画像)が生成される。制御部31は、このような仮想横断層画像を、センサ部12が移動した移動範囲内において1回転した位置毎に生成してもよく、所定間隔毎に生成してもよい。このように生成された仮想横断層画像は、横断層用モデルM1の学習に用いてもよい。即ち、各撮影位置での横断層画像に対して、各撮影位置について生成された仮想横断層画像を正解のラベル画像とする訓練データを、横断層用モデルM1の訓練データに用いてもよい。
【0060】
制御部31は、ステップS32で生成した縦断層画像と、上述のように生成した仮想横断層画像とに基づいて、画像中の血管の状態を検出する第1処理を実行する(S35)。例えば制御部31は、各撮影位置における仮想横断層画像について、画像中の血管内のプラーク領域を特定する。例えば制御部31は、各仮想横断層画像中の血管内腔の面積と、血管内腔及び血管壁を合わせた血管全体の面積とを算出し、例えば以下の式を用いてプラークバーデンを算出し、プラークバーデンが所定値(例えば50%)以上である撮影位置を、プラーク領域がある位置に特定する。なお、制御部31は、仮想横断層画像中の各画素の画素値に基づいて、プラークの領域を特定してもよい。また、OCT画像では脂質性組織で形成されて線維性被膜が薄い脂質性プラークを高精度で撮影できるため、制御部31は、OCT画像から生成した仮想横断層画像に基づいて脂質性プラークの領域を特定してもよい。IVUS画像では超音波の減衰を伴う減衰性プラーク(attenuated plaque)を撮影できるため、制御部31は、IVUS画像から生成した仮想横断層画像に基づいて減衰性プラークの領域を特定してもよい。また、制御部31は、仮想横断層画像においてプラーク領域の血管周方向の広がりを計測してもよい。
【0061】
プラークバーデン=(血管全体の面積-血管内腔の面積)/血管全体の面積 …(式)
【0062】
また制御部31は、仮想横断層画像に基づいて、画像中の血管内の石灰化領域を特定してもよい。また制御部31は、例えばIVUS画像から生成した仮想横断層画像に基づいて、画像診断用カテーテル1が挿入されている血管(本幹)から分岐して伸びる血管(側枝)の位置を検出してもよい。また、ステップS31で取得したラインデータが、ステント留置後に取得したラインデータである場合、制御部31は、ステップS32で生成した縦断層画像に基づいて、画像中の血管内のステントの留置範囲を検出してもよい。なお、制御部31は、IVUS及びOCTの各モダリティが得意とする検出対象について、IVUS縦断層画像及びOCT縦断層画像のそれぞれに対するセグメンテーションの結果から、IVUS横断層画像及びOCT横断層画像のそれぞれに対するセグメンテーションを実行するか否かを判断してもよい。また、横断層画像(IVUS横断層画像及び/又はOCT横断層画像)に対するセグメンテーションを実行する場合、縦断層画像に対するセグメンテーションの結果と、横断層画像に対するセグメンテーションの結果との間にズレが生じていれば、横断層画像及び縦断層画像に基づいてセグメンテーションの結果を補正してもよい。例えば、縦断層画像に対するセグメンテーションの結果から生成された仮想横断層画像と、横断層画像に対するセグメンテーションの結果のラベル画像とにおいて、各画素の平均値を補正後のラベル画像の画素値としてもよい。また、撮影位置の前後方向のズレも加味して、縦断層画像に対するセグメンテーションの結果と、横断層画像に対するセグメンテーションの結果との間のズレを補正してもよい。
【0063】
制御部31は、ステップS32で構築して主記憶部32又は補助記憶部35に記憶した縦断層画像を表示装置4に表示する(S36)。ここでは、制御部31は、
図11Aに示すような画面を表示する。
図11Aの画面は、OCT縦断層画像及びIVUS縦断層画像を並べて表示しており、この2つの縦断層画像は、血管の径方向の同じ方向を撮影方向とした縦断層画像である。また
図11Aの画面では、OCT縦断層画像の上側に、撮影対象の血管の長軸方向において、プラーク領域(例えば脂質性プラーク)が検出された仮想横断層画像の撮影位置を示すマークL1が表示されており、IVUS縦断層画像の上側に、プラーク領域(例えば減衰性プラーク)が検出された仮想横断層画像の撮影位置を示すマークL2が表示されている。なお、マークL1は、例えばプラークの血管周方向の広がりが180度以上であると判定された仮想横断層画像の撮影位置を示している。また、マークL2は、例えばプラークバーデンが50%以上であると判定された仮想横断層画像の撮影位置、プラークの血管周方向の広がりが180度以上であると判定された仮想横断層画像の撮影位置、石灰化領域の血管周方向の広がりが180度以上であると判定された仮想横断層画像の撮影位置等を示しており、血管の状態に応じてマークL2内のハッチングを異ならせてもよい。また、
図11Aの画面では、IVUS縦断層画像上に、側枝が分岐している位置を示すマークL3が表示されている。
【0064】
ステップS31で取得したラインデータがステント留置後に取得したラインデータである場合、制御部31は、
図11Bに示すような画面を表示する。
図11Bの画面は、
図11Aの画面と同様に、OCT縦断層画像及びIVUS縦断層画像を並べて表示している。また
図11Bの画面では、OCT縦断層画像及びIVUS縦断層画像の上側に、血管の長軸方向におけるステントの留置範囲を示すマークL4と、MSA(Minimum Stent Area:最小ステント拡張面積)の位置を示すマークL6とが表示されている。更に
図11Bの画面では、ステントと血管内腔面(内膜)との隙間(マルアポジション)が所定値(例えば0.4mm)以上である位置を示すマークL5が表示されている。このような画面を表示装置4に表示することにより、縦断層画像及び仮想横断層画像に基づく第1処理によって検出した血管の状態を術者に提示できる。ここで、画像処理装置3は、横断層画像として例えば512画素×512画素の画像を処理しており、血管内検査装置101が1回のプルバック操作によって1024回転する場合、縦断層画像として例えば512画素×1024画素の画像を処理する。よって、縦断層画像は横断層画像の2枚分のデータ量であり、このような少ないデータ量に基づいて、血管の長軸方向の概要を検出することにより、センサ部12の撮影終了から早期に血管の概要を術者に提示できる。なお、縦断層画像は、512画素×1024画素の画像のほかに、血管の長軸方向の画素を1画素毎に間引くことにより、512画素×512画素の画像としてもよく、この場合、更にデータ量を削減することができる。なお、512画素×1024画素の縦断層画像を512画素×512画素の画像2枚で生成してもよく、この場合、2枚の画像のつなぎ目部分を含む画像を更に1枚追加して3枚の画像で縦断層画像を生成してもよい。例えば1フレーム~1024フレームの各撮影位置において、1フレーム目~512フレーム目のラインデータに基づく縦断層画像と、256フレーム目~768フレーム目のラインデータに基づく縦断層画像と、512フレーム目~1024フレーム目のラインデータに基づく縦断層画像とを生成してもよい。
【0065】
次に制御部31は、縦断層画像及び仮想横断層画像に基づく第1処理の結果から、詳細に血管状態を確認すべき優先領域を特定する(S37)。例えば制御部31は、
図10Bに示すように、上述した第1処理によって検出した、プラークバーデンが所定値以上の位置、血管内腔の面積又は直径が所定値未満の位置(狭窄部の位置)、及び、これらの位置を含む近傍の領域等を優先領域に特定する。また、ステップS31で取得したラインデータがステント留置後に取得したラインデータである場合、制御部31は、
図10Cに示すように、ステントの留置範囲を優先領域に特定してもよい。本実施形態では、第1処理の結果から優先領域を特定するが、優先領域は、縦断層画像(IVUS縦断層画像及び/又はOCT縦断層画像)に対する、縦断層用モデルM2を用いたセグメンテーションの結果から特定されればよい。例えば、縦断層画像(IVUS縦断層画像及び/又はOCT縦断層画像)を入力した場合に、入力された縦断層画像中の優先領域を出力するモデルを用いて、優先領域を特定してもよい。
【0066】
制御部31は、ステップS31で取得したラインデータのうちで、特定した優先領域内を撮影位置として取得したラインデータに基づいて、当該優先領域を撮影位置とする横断層画像を構築する(S38)。ここでは制御部31は、優先領域内でセンサ部12が一回転する間に取得された512本のラインデータ毎に、512本のラインデータに対して補間処理によって画素を補間して2次元の横断層画像を構築する。制御部31は、構築した横断層画像を主記憶部32又は補助記憶部35に記憶する。
【0067】
制御部31は、ステップS38で構築して主記憶部32又は補助記憶部35に記憶した優先領域内の横断層画像のそれぞれに対してセグメンテーションを実行する(S39)。ここでは制御部31は、それぞれの横断層画像を横断層用モデルM1に入力し、横断層用モデルM1から出力されたラベル画像に基づいて、それぞれの横断層画像中の血管内腔及び血管壁の領域を特定する。そして、制御部31は、セグメンテーションによって得られた各横断層画像中の血管内腔及び血管壁の領域に基づいて、画像中の血管の状態を検出する第2処理を実行する(S40)。第2処理は、処理対象が横断層画像であるが、第1処理と同じ処理である。よって、制御部31は、各横断層画像について、画像中の血管内のプラーク領域及び石灰化領域を特定し、プラーク領域及び石灰化領域の血管周方向の広がりを計測し、側枝の位置を検出する。また制御部31は、各横断層画像中の血管内腔面(ルーメンの輪郭)、EEM(External Elastic Membrane:外弾性板)(ヴェッセルの輪郭)を検出してもよく、ステントの留置後のラインデータに基づく横断層画像については画像中の血管内のステントを検出してもよい。
【0068】
制御部31は、構築して主記憶部32又は補助記憶部35に記憶した横断層画像を表示装置4に表示する(S41)。ここでは、制御部31は、
図11Cに示すような画面を表示する。
図11Cの画面は、ステップS32で生成したOCT縦断層画像及びIVUS縦断層画像と、ステップS38で生成したOCT横断層画像及びIVUS横断層画像とを並べて表示している。なお、OCT横断層画像及びIVUS横断層画像は、血管の長軸方向の同じ位置で撮影された横断層画像であり、この横断層画像の撮影位置は、OCT縦断層画像及びIVUS縦断層画像においてマークL7で示されている。
図11Cに示す画面において、例えば入力装置5を介して一方の縦断層画像上のマークL7が左右に移動された場合、制御部31は、他方の縦断層画像上のマークL7も連動して移動し、また、表示中の横断層画像を、移動後のマークL7の位置を撮影位置とする横断層画像に切り替える。これにより、術者は、横断層画像を、観察位置(撮影位置)を適宜変更しながら確認できる。
【0069】
また
図11Cの画面では、
図11Aの画面中のマークL1~L3に加えて、OCT横断層画像及びIVUS横断層画像上に、表示中の横断層画像中に検出されたプラーク(例えば脂質性プラーク)の領域を示すマークL8,L9が表示されている。マークL8,L9は、表示中の横断層画像中のプラーク領域の血管周方向の広がりを示している。これにより、術者は、脂質性プラーク等のプラークの位置及び広がり範囲を把握できる。また、
図11Cの画面において、横断層画像上に、石灰化領域の位置及び広がり範囲を示すマークが表示されてもよい。更に、横断層画像上に、横断層画像から検出された血管内腔面(ルーメンの輪郭)、EEM(ヴェッセルの輪郭)、ステントが表示されてもよい。このような画面を表示装置4に表示することにより、横断層画像に基づく第2処理によって検出した血管の状態を術者に提示できる。
【0070】
制御部31は、優先領域を撮影位置とする横断層画像に対してステップS39~S41の処理を実行した後、優先領域以外の領域を撮影位置として取得したラインデータに基づいて、ステップS38~S40と同様の処理を行う(S42~S44)。具体的には、制御部31は、優先領域以外の領域を撮影位置とする横断層画像を構築し(S42)、構築した横断層画像のそれぞれに対してセグメンテーションを実行し(S43)、セグメンテーションによって得られた各横断層画像中の血管内腔及び血管壁の領域に基づいて、画像中の血管の状態を検出する第2処理を実行する(S44)。なお、制御部31は、ステップS38~S39では、センサ部12が一回転する間に取得された512本のラインデータ毎に、即ち1フレーム毎に横断層画像の生成及び各横断層画像のセグメンテーションを行うのに対して、ステップS42~S43では、複数フレーム毎(例えば0.5mm毎又は1mm毎)に横断層画像の生成及びセグメンテーションを行ってもよい。この場合、詳細に血管状態を観察すべき優先領域では、横断層画像の各フレームに対して第2処理が行われて血管の状態が詳細に検出されるのに対して、優先領域以外の領域では、第2処理を行うフレームを間引いて複数フレーム毎に第2処理が行われる。即ち、優先領域以外の領域では、第2処理を実行する横断層画像の撮影位置の間隔(第2処理の処理間隔)が複数フレーム毎であり、優先領域よりも大きく設定されている。これにより、優先領域では1フレーム毎の検出処理によって血管の状態を高精度に検出すると共に、優先領域以外の領域では処理対象のフレーム数を削減することにより、演算量を削減することができる。
【0071】
制御部31は、ステップS44の処理後、ステップS42で生成した横断層画像と、ステップS44で検出した血管の状態を示す各マークとを表示装置4に表示してもよい。なお、術者は、ステップS41で表示された
図11Cに示す画面によって血管の状態を確認し、例えばバルーンカテーテルによる血管拡張前又はステントの留置前であれば、拡張すべき位置及びステントの留置位置を決定し、ステントの留置後であればステントの留置位置が適切であるか否かを判断する。よって、ステップS41の処理後、術者は、決定した処置内容に基づいて次の処置を開始してもよく、この場合、制御部31は、ステップS42~S44の処理を行わずに一連の処理を終了してもよい。
【0072】
上述した処理により、本実施形態では、プルバック操作によって一連のラインデータを取得した後、縦断層画像に基づいて血管の長軸方向の概要が検出され、血管全体の状態をおおまかに把握するための情報が術者に提供される。縦断層画像はラインデータの一部(任意のライン番号のラインデータ及びこのライン番号に256を加算したライン番号のラインデータ)から構築されるので、全てのラインデータを処理する必要がなく、早期に血管状態の検出が可能となる。よって、ラインデータの取得後すぐに、縦断層画像に基づく演算量の少ない処理によって、血管全体の状態を大まかに把握できる情報を提供できる。よって、術者は、ラインデータの取得後、早期に提供された情報に基づいて、注目すべき領域(狭窄部、病変部等)を早急に特定することが可能となる。また、本実施形態では、注目すべき領域として、縦断層画像に基づいて優先領域が設定され、優先領域を撮影位置とする横断層画像に対する血管状態の検出処理を優先的に行うことにより、注目すべき領域についての検査結果が早期に術者に提供される。このように縦断層画像から検出された血管の状態に基づいて、優先的に検出処理を行うべき領域を決定できるので、長軸方向の各領域に対して、横断層画像に基づく詳細な検出処理を行う際の優先順位を設定できる。また、優先領域については1フレーム毎に検出処理を実行することにより、詳細な検出結果が得られる。一方、注目すべき領域ではない領域については、検出処理の対象を削減することにより、処理負荷の増加を抑制できる。よって、本実施形態では、ラインデータの取得直後に、注目すべき領域の特定に利用できる情報として、縦断層画像に基づいて検出された血管の情報が提示され、注目すべき領域が設定された後は、注目すべき領域で撮影された横断層画像に基づいて血管の詳細な状態が検出される。このように、術者が意思決定を行う際に有益な情報を適切なタイミングで術者に提供できる。
【0073】
図9に示す処理は、ステップS31で超音波ラインデータ及び光ラインデータを取得し、超音波ラインデータ及び光ラインデータのそれぞれについてステップS32~S44の処理を行うものである。本実施形態では、超音波ラインデータに基づく検出処理の結果と、光ラインデータに基づく検出処理の結果とから、最終的な検出結果を判断してもよい。例えばステップS35において、IVUS画像に基づく検出結果と、OCT画像に基づく検出結果とに基づいて、最終的な血管の状態を判断してもよい。また、ステップS37において、IVUS画像に基づいて特定した優先領域と、OCT画像に基づいて特定した優先領域とに基づいて、最終的な優先領域を特定してもよい。また、ステップS40,S44において、IVUS画像に基づく検出結果と、OCT画像に基づく検出結果とに基づいて、最終的な血管の状態を判断してもよい。このような構成とした場合、より高精度の検出処理を実現できる。
【0074】
(実施形態2)
横断層画像に基づく血管状態の検出処理(第2処理)の処理対象の領域を術者が指定できる画像診断装置について説明する。本実施形態の画像診断装置100は、実施形態1の画像診断装置100における各装置と同様の装置によって実現できるので、同様の構成については説明を省略する。
【0075】
図12は、実施形態2の画像の表示処理手順の一例を示すフローチャートである。
図12に示す処理は、
図9に示す処理において、ステップS36,S37の間にステップS51~S55を追加したものである。
図9と同じステップについては説明を省略する。
【0076】
本実施形態の画像診断装置100では、画像処理装置3の制御部31は、
図9中のステップS31~S36と同じ処理を実行する。これにより、例えば
図11A又は
図11Bに示すような画面が表示装置4に表示される。本実施形態の画像処理装置3は、例えば入力装置5を介して、横断層画像に基づく血管状態の検出処理(第2処理)の実行を優先的に行う処理領域の指定を受け付けるように構成されている。処理領域の指定は、例えば
図11A又は
図11Bの画面中の縦断層画像に対して、血管の長軸方向の任意の領域が選択されることによって行われる。制御部31は、処理領域の指定を受け付けたか否かを判断する(S51)。処理領域の指定を受け付けていないと判断した場合(S51:NO)、制御部31は、ステップS37の処理に移行する。この場合、制御部31は、
図9に示す実施形態1の処理と同じ処理を行う。
【0077】
処理領域の指定を受け付けたと判断した場合(S51:YES)、制御部31は、指定された処理領域(指定領域)内を撮影位置として取得したラインデータに基づいて、ステップS38~S41と同様の処理を行う(S52~S55)。ここでは制御部31は、指定領域を撮影位置とする横断層画像を構築し(S52)、構築した横断層画像のそれぞれに対してセグメンテーションを実行し(S53)、セグメンテーションによって得られた各横断層画像中の血管内腔及び血管壁の領域に基づいて、画像中の血管の状態を検出する第2処理を実行する(S54)。なお、ステップS52~S53では、制御部31は、1フレーム毎に横断層画像の生成及び各横断層画像のセグメンテーションを行う。これにより、術者に指定された領域に対して血管状態を高精度に検出することができる。
【0078】
そして、制御部31は、
図11Cに示すような画面によって、ステップS52で生成した横断層画像を表示装置4に表示する(S55)。その後、制御部31は、ステップS42の処理に移行する。なお、制御部31は、ステップS55の処理後、ステップS37の処理に移行してもよい。この場合、術者が指定した領域に対する検出処理が終了した後に、仮想横断層画像に基づく血管状態の検出結果から優先領域が特定され、特定した優先領域に対して横断層画像に基づく血管状態の検出処理が行われる構成となる。
【0079】
上述した処理により、本実施形態では、横断層画像に基づく血管状態の検出処理(第2処理)を優先的に行う領域の指定を術者から受け付けることができる。よって、術者が指定した領域に対して、横断層画像に基づく血管状態の検出処理を優先的に実行することができる。これにより、縦断層画像に基づく血管状態の検出処理から設定される優先領域よりも、術者による指定領域に対する処理を割込処理として優先的に行うことができる。よって、術者が観察したい領域に対する検出結果を早期に提示することができる。本実施形態においても、上述した実施形態1で適宜説明した変形例の適用が可能である。
【0080】
(実施形態3)
横断層用モデルM1の変形例について説明する。
図13は横断層用モデルM1の変形例を示す説明図である。
図13に示す横断層用モデルM1は、
図6Aに示す実施形態1の横断層用モデルM1の構成に加えて、
図9に示す処理中のステップS34で生成される仮想横断層画像(IVUS仮想横断層画像又はOCT仮想横断層画像)が入力される構成を有する。具体的には、横断層用モデルM1は、横断層画像(IVUS横断層画像又はOCT横断層画像)と、仮想横断層画像(IVUS仮想横断層画像又はOCT仮想横断層画像)とを入力とし、入力された横断層画像中の血管内の血管内腔の領域及び血管壁の領域を認識するモデルである。
【0081】
本実施形態の横断層用モデルM1は、1枚の横断層画像と、この横断層画像と同じ撮影位置であり、
図9中のステップS34で生成された仮想横断層画像とを入力とし、入力された横断層画像及び仮想横断層画像に基づいて、当該横断層画像に含まれる血管内腔の領域及び血管壁の領域を認識する演算を行い、認識結果を示す情報を出力するように学習される。具体的には、横断層用モデルM1は、入力された横断層画像の各画素を、血管内腔の領域と血管壁の領域とその他の領域とに分類し、各画素に領域毎のラベルを対応付けた分類済みの画像(ラベル画像)を出力する。
【0082】
本実施形態の横断層用モデルM1は、訓練用の横断層画像及び仮想横断層画像と、この横断層画像中の各画素に対して判別すべきオブジェクト(ここでは血管内膜及び血管壁の領域)を示すデータがラベリングされた正解のラベル画像とを含む訓練データを用いて機械学習することにより生成される。本実施形態の横断層用モデルM1は、訓練データに含まれる横断層画像及び仮想横断層画像が入力された場合に、訓練データに含まれる正解のラベル画像を出力するように学習する。具体的には、横断層用モデルM1は、入力された横断層画像中の各画素に対して、分類されたオブジェクトの種類を示す値がラベリングされたラベル画像を出力として取得する。そして、横断層用モデルM1は、検出結果のラベル画像を、正解のラベル画像におけるオブジェクト領域の座標範囲及びオブジェクトの種類と比較し、両者が近似するように、最急降下法、誤差逆伝播法等を用いて、ニューロン間の重み等のパラメータを最適化する。これにより、横断層画像及び仮想横断層画像が入力された場合に、入力画像中の血管内腔の領域及び血管壁の領域を示すラベル画像を出力する横断層用モデルM1が得られる。本実施形態の横断層用モデルM1においても、横断層画像及び仮想横断層画像から、血管内腔及び血管壁の領域に加えて、ガイドワイヤ及びカテーテルの領域を認識する構成でもよく、血管壁中のプラークの領域を認識する構成でもよい。
【0083】
上述した構成の横断層用モデルM1を用いる場合であっても、画像処理装置3は、
図9の処理の実行が可能である。なお、ステップS39において、制御部31は、ステップS38で構築した優先領域内の横断層画像のそれぞれと、各横断層画像と同じ撮影位置の仮想横断層画像とを横断層用モデルM1に入力し、横断層用モデルM1から出力されたラベル画像に基づいて、それぞれの横断層画像中の血管内腔及び血管壁の領域を特定する。また、ステップS43において、制御部31は、ステップS42で構築した優先領域以外の領域内の横断層画像のそれぞれと、各横断層画像と同じ撮影位置の仮想横断層画像とを横断層用モデルM1に入力し、横断層用モデルM1から出力されたラベル画像に基づいて、それぞれの横断層画像中の血管内腔及び血管壁の領域を特定する。
【0084】
本実施形態の横断層用モデルM1では、横断層画像と仮想横断層画像とが入力される構成であるが、例えば、横断層画像と仮想横断層画像とを重ねた状態で入力される構成でもよい。例えば、横断層画像に対して、仮想横断層画像中の血管内腔の輪郭線と血管壁の輪郭線とが付加された画像を入力画像としてもよく、横断層画像中の各画素に対して、各画素が対応する仮想横断層画像中の画素が血管内腔の領域であるか血管壁の領域であるかを示すラベルが付加された画像を入力画像としてもよい。
【0085】
本実施形態の横断層用モデルM1の構成は、上述した実施形態1~2の横断層用モデルM1に適用可能であり、実施形態1~2の横断層用モデルM1に適用した場合であっても同様の効果が得られる。なお、本実施形態の横断層用モデルM1のように、横断層画像に加えて、当該横断層画像と同じ撮影位置について生成された仮想横断層画像に基づいて、当該横断層画像中の血管内腔及び血管壁の領域を分類する構成により、より分類精度の高いモデルを実現できる。また、本実施形態においても、上述した実施形態1~2で適宜説明した変形例の適用が可能である。
【0086】
上述した各実施形態では、超音波を用いて血管内の断層画像を撮影するIVUSセンサ12aと、近赤外光を用いて血管内の断層画像を撮影するOCTセンサ12bとを使用する構成であるが、このような構成に限定されない。例えばIVUSセンサ12a又はOCTセンサ12bの代わりに、血管内からのラマン散乱光を受光して血管内の断層画像を撮影するセンサ、血管内からの励起光を受光して血管内の断層画像を撮影するセンサ等、血管の状態を観察できる各種のセンサを使用する構成とすることもできる。
【0087】
今回開示された実施の形態はすべての点で例示であって、制限的なものではないと考えられるべきである。本発明の範囲は、上記した意味ではなく、特許請求の範囲によって示され、特許請求の範囲と均等の意味及び範囲内でのすべての変更が含まれることが意図される。
【符号の説明】
【0088】
1 画像診断用カテーテル
2 MDU
3 画像処理装置
4 表示装置
5 入力装置
31 制御部
32 主記憶部
33 入出力部
35 補助記憶部
101 血管内検査装置
102 血管造影装置
P プログラム