JP7412223B2 - 画像処理装置、医用画像診断装置、画像処理方法、プログラム、および学習装置 - Google Patents

Description

本発明は、画像処理装置、医用画像診断装置、画像処理方法、プログラム、および学習装置に関する。

医用の分野では、種々の医用画像診断装置（モダリティ）によって取得される断面画像（断層像とも称す）を用いた診断が行われている。この診断の中では、断面画像の種別（被検体のどの断面を撮像した画像であるかによる分類。以下、断面種別と称す）を識別し、断面種別を表す情報を断面画像に関連づけて表示または保存することが行われている。また、識別した断面種別に応じた画像処理を当該断面画像に施すことや、識別した情報を利用した解析処理を実行することが行われている。この断面種別を識別する作業は一般に、人手で行う場合が多い。

特許文献１では、心臓を撮像した２次元超音波画像を対象として、その断面種別を自動的に識別する技術が開示されている。特許文献１では、一例として、心尖部アプローチ像の断面種別を左室長軸断面、四腔断面、五腔断面、二腔断面のなかから識別する例が記載されている。この技術では、各断面種別の代表例となる画像をテンプレートとして保持し、入力画像がどのテンプレートに近いかを算出することで断面種別を識別する。

国際公開第２００７／０５８１９５号

しかしながら、特許文献１では単一の症例（断面画像）のみをテンプレートとして用いているため、断面画像の多様なバリエーションに対応するには不十分であるという課題がある。

例えば、被検体（患者）によって臓器等の形状が異なるため、同じ断面種別の画像であっても、実際の断面画像は被検体によって異なっている。また、医用画像診断装置の機種・性能や撮像パラメータ（撮像条件）によっても、得られる断面画像に違いがある。それゆえ、断面種別の高精度な識別を実現するためには、このような断面画像の多様なバリエーション（ばらつき）への十分な対応が重要となる。

本発明は上記課題に鑑みてなされたものであり、断面画像の種別を精度よく識別するための技術を提供することを目的とする。

なお、前記目的に限らず、後述する発明を実施するための形態に示す各構成により導かれる作用効果であって、従来の技術によっては得られない作用効果を奏することも本明細書の開示の他の目的の１つとして位置付けることができる。

本発明の第一態様は、被検体の断面画像から当該断面画像の種別である断面種別を識別する画像処理装置であって、処理対象となる入力画像を取得する画像取得部と、断面種別が同じ断面画像を主成分分析することによって生成された、複数の断面種別に対応する複数の統計モデルを用いて、前記入力画像の断面種別を識別する識別部と、を有し、前記統
計モデルは、断面種別が同じ複数の断面画像の特徴を表す部分空間の情報であることを特徴とする画像処理装置を提供する。

本発明の第二態様は、被検体の断面画像から当該断面画像の種別である断面種別を識別する医用画像診断装置であって、処理対象となる医用画像を取得する画像取得部と、断面種別が同じ断面画像を主成分分析することによって生成された、複数の断面種別に対応する複数の統計モデルを用いて、前記医用画像の断面種別を識別する識別部と、を有し、前記統計モデルは、断面種別が同じ複数の断面画像の特徴を表す部分空間の情報であることを特徴とする医用画像診断装置を提供する。

本発明の第三態様は、被検体の断面画像から当該断面画像の種別である断面種別を識別する画像処理方法であって、処理対象となる入力画像を取得するステップと、断面種別が同じ断面画像を主成分分析することによって生成された、複数の断面種別に対応する複数の統計モデルを用いて、前記入力画像の断面種別を識別するステップと、を有し、前記統計モデルは、断面種別が同じ複数の断面画像の特徴を表す部分空間の情報であることを特徴とする画像処理方法を提供する。

本発明の第四態様は、上記画像処理方法の各ステップをコンピュータに実行させるためのプログラムを提供する。

本発明の第五態様は、被検体の断面画像から当該断面画像の種別である断面種別を識別する処理に用いられるモデルを学習するための学習装置であって、複数の断面種別のそれぞれについて、複数の断面画像を学習データとして準備する学習データ準備部と、断面種別が同じ断面画像を主成分分析することによって、複数の断面種別に対応する複数の統計モデルを生成するモデル生成部と、を備え、前記統計モデルは、断面種別が同じ複数の断面画像の特徴を表す部分空間の情報であることを特徴とする学習装置を提供する。

本発明によれば、断面画像の種別を精度よく識別することが可能となる。

第１の実施形態に係る画像処理システムの構成例を示す図。 第１の実施形態に係る画像処理装置の処理手順の例を示すフローチャート。 心尖部二腔像、心尖部三腔像、心尖部四腔像、空中放置像の例を示す図。 学習データから統計モデルを生成する学習処理の例を示すフローチャート。 固有ベクトルの数と累積寄与率の関係の例を示すグラフ。 断面種別の識別処理の例を示すフローチャート。 第２の実施形態に係る画像処理システムの構成例を示す図。 第２の実施形態に係る画像処理装置の処理手順の例を示すフローチャート。

以下、図面を参照して、本発明の実施形態を例示的に詳しく説明する。ただし、この実施形態に記載されている構成要素はあくまで例示であり、本発明の技術的範囲は、特許請求の範囲によって確定されるのであって、以下の個別の実施形態によって限定されるわけではない。

本発明の実施形態に係る画像処理装置は、被検体の断面画像から当該断面画像の種別である断面種別を推定ないし識別する機能を提供する。断面画像は、モダリティと呼ばれる撮像装置によって取得された、被検体（人体など）の内部の構造を表す医用画像であり、断層像、再構成画像などとも呼ばれる。例えば、超音波診断装置によって得られる超音波画像、Ｘ線ＣＴ装置によって得られるＸ線ＣＴ画像、ＭＲＩ装置によって得られるＭＲＩ画像などが断面画像の典型例である。このような断面画像は、例えば、医用分野において、診断、検査、研究などに利用される。断面種別とは、被検体のどの断面を撮像した画像
であるかという観点による分類である。ここでは、画像処理装置を例にして説明するが、医用画像診断装置に本発明を適用してもよい。

本発明の実施形態に係る画像処理装置は、断面種別が既知の多数の断面画像（学習データ）を統計解析した結果に基づいて、処理対象となる入力画像（断面種別が未知の断面画像）の断面種別の推定ないし識別を行う点に特徴の一つを有する。具体的には、画像処理装置は、複数の断面種別にそれぞれ対応する複数の統計モデルを用い、複数の統計モデルのいずれが入力画像に良く当てはまるかを評価することによって入力画像の断面種別を識別する。例えば、夫々の統計モデルを用いて入力画像を再構築し、再構築された画像と入力画像の差異が最も小さい、言い換えると、入力画像を最も良く再現する統計モデルと対応付けられた断面種別を識別結果とする。

ここで、統計モデルとは、断面種別が同じ複数の断面画像を統計解析することによって生成されたモデルであり、その断面種別の画像の特徴の分布（統計的傾向）を表現するものである。統計モデルの生成（学習）は、例えば、学習装置によって、断面種別が同じ画像群を主成分分析し、その断面種別の画像群の特徴を表す部分空間を求めることにより行ってもよい。この場合、部分空間の情報が、当該断面種別に対応する統計モデルに該当する。なお、統計モデルは学習済モデルとも呼ばれる。

従来技術では、単一の症例（画像）のみから生成したテンプレートを用いていたため、被検体、撮像装置、撮像条件等による断面画像のばらつきに起因して、入力画像とテンプレートの間の特徴が乖離し、識別精度が低下するという課題があった。これに対し本発明の実施形態に係る画像処理装置では、複数の断面画像から生成した統計モデルを用いることで、多数症例の統計的傾向に基づく識別が可能であり、従来技術に比べて識別精度の向上を図ることができる。

＜第１の実施形態＞
第１の実施形態に係る画像処理装置は、入力画像の断面種別を識別する装置である。本実施形態では、入力画像および学習データとして心臓の２次元超音波画像を用い、入力画像が心尖部二腔像、心尖部三腔像、心尖部四腔像、および空中放置像の４種類のいずれであるかを識別する、というケースを例示する。心尖部二腔像は、左心房・左心室の２腔が写った画像であり、心尖部三腔像は、左心房・左心室・右心室の３腔が写った画像であり、心尖部四腔像は、左右心房・左右心室の計４腔が写った画像である。また、空中放置像は、プローブが被検体に接しておらずいずれの腔も写っていない画像である。なお、これらは断面種別の一例であり、その他の種別の画像についても、同様の処理で識別可能である。

以下、図１を用いて本実施形態の画像処理装置の構成および処理を説明する。図１は、本実施形態の画像処理装置を含む画像処理システム（医用画像処理システムともいう）の構成例を示すブロック図である。画像処理システムは、画像処理装置１０およびデータベース２２を備える。画像処理装置１０は、ネットワーク２１を介してデータベース２２に通信可能な状態で接続されている。ネットワーク２１は、例えばＬＡＮ（Ｌｏｃａｌ Ａｒｅａ Ｎｅｔｗｏｒｋ）やＷＡＮ（Ｗｉｄｅ Ａｒｅａ Ｎｅｔｗｏｒｋ）を含む。

データベース２２は、多数の画像と各画像に関連付けられた情報とを保持し、管理する。データベース２２に保持される画像は、断面種別が未知の画像と、断面種別が既知の画像とを含む。前者の画像は、画像処理装置１０による断面種別の識別処理に供される。後者の画像は、少なくとも断面種別の情報（正解データ）が関連付けられており、統計モデルを生成するための学習データとして利用される。各画像には、超音波診断装置が当該画像を撮像した際の撮像パラメータ（撮像条件）などの情報が関連付けられていてもよい。
画像処理装置１０は、ネットワーク２１を介して、データベース２２が保持しているデータを取得することが可能である。

画像処理装置１０は、通信ＩＦ（Ｉｎｔｅｒｆａｃｅ）３１（通信部）、ＲＯＭ（Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ）３２、ＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）３３、記憶部３４、操作部３５、表示部３６、および制御部３７を備える。

通信ＩＦ３１（通信部）はＬＡＮカードなどにより構成され、外部装置（例えばデータベース２２など）と画像処理装置１０との通信を実現する。ＲＯＭ３２は、不揮発性のメモリなどにより構成され、各種プログラムや各種データを記憶する。ＲＡＭ３３は揮発性のメモリなどにより構成され、実行中のプログラムや作業中のデータを一時的に記憶する。記憶部３４はＨＤＤ（Ｈａｒｄ Ｄｉｓｋ Ｄｒｉｖｅ）などにより構成され、プログラムやデータを記憶する。操作部３５は、キーボードやマウス、タッチパネルなどにより構成され、ユーザ（例えば医師や検査技師）からの指示を各種装置に入力する。

表示部３６はディスプレイなどにより構成され、各種情報をユーザに提示する。制御部３７は、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）などにより構成され、画像処理装置１０における処理を統括制御する。制御部３７は、その機能的な構成として、画像取得部５１、統計モデル取得部５２、識別部５３、表示処理部５４を備える。制御部３７は、ＧＰＵ（Ｇｒａｐｈｉｃｓ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）やＤＳＰ（Ｄｉｇｉｔａｌ Ｓｉｇｎａｌ Ｐｒｏｃｅｓｓｏｒ）やＦＰＧＡ（Ｆｉｅｌｄ－Ｐｒｏｇｒａｍｍａｂｌｅ Ｇａｔｅ Ａｒｒａｙ）などを備えてもよい。

画像取得部５１は、処理対象となる入力画像をデータベース２２から取得する。入力画像は、モダリティにより撮像された被検体の断面画像である。なお、画像取得部５１は、入力画像をモダリティから直接取得してもよい。この場合、画像処理装置１０はモダリティに接続された装置であってもよいし、モダリティの中（例えばコンソール）に実装された機能であってもよい。

また、画像取得部５１は、学習データをデータベース２２から取得する。学習データは、入力画像とは異なる複数の断面画像（被検体画像）からなる画像群であり、識別対象とする断面種別の夫々に属する断面画像を複数枚含んでいる。学習データにおける夫々の画像は、その断面種別の情報（心尖部二腔像・心尖部三腔像・心尖部四腔像・空中放置像のいずれかを表す識別子）を既知の値（付帯情報）として保持している。なお、統計モデルのロバスト性を高めるために、学習データとして用いる複数の断面画像は、互いに異なる被検体（患者）の画像であることが望ましい。ただし、同一の被検体（患者）を異なる条件（時刻、時期、モダリティの機種、撮像パラメータ等）で撮影した画像が学習データに含まれていてもよい。また、１心拍周期にわたって撮像された２次元動画像の夫々のフレームを学習データとして用いてもよい。学習データとして用いる断面画像の数は少なくとも数十から数百程度あるとよい。なお、学習データにおける断面種別情報の保持の仕方は、夫々の断面画像に付帯情報として断面種別情報を保持させる方法以外の方法でもよい。例えば、断面画像の識別子（画像ＩＤ、ファイル名など）と断面種別とが対応付けられたリストを用いたり、断面種別ごとに断面画像を格納するフォルダを分けたり、断面画像のファイル名やプロパティに断面種別情報を埋め込むなど、いかなる方法でもよい。

本実施形態では入力画像および学習データが心臓領域を撮像した１枚の２次元超音波画像である場合を例に挙げて説明するが、他の種類（他の臓器、他の断面種別）の画像であってもよい。本発明は、複数の２次元画像（例えば、複数の時相の画像）や２次元の動画像に対しても適用可能である。また、モダリティの種類によらず適用可能である。

統計モデル取得部５２は、画像取得部５１が取得した学習データに含まれる断面画像の夫々について画素値情報と断面種別情報を取得し、断面種別が一致する画像群毎に（すなわち、夫々の断面種別毎に）統計解析を行う。すなわち、本実施形態では、統計モデル取得部５２は、学習データに含まれる心尖部二腔像の画像群・心尖部三腔像の画像群・心尖部四腔像の画像群・空中放置像の画像群に対して、それぞれ統計解析を行う。そして、統計モデル取得部５２は、該統計解析の結果から、夫々の断面種別の統計モデルとして、夫々の断面種別の部分空間情報（部分空間を構成する基底の情報）を取得する。なお、あらかじめ生成された統計モデルがデータベース２２又は記憶部３４に保持されている場合は、統計モデル取得部５２は、データベース２２又は記憶部３４から統計モデルを取得してもよい。すなわち、統計モデル取得部５２は、学習データから統計モデルを生成する学習装置として機能してもよいし、あらかじめ生成され記憶装置に格納されている統計モデルを読み込む機能であってもよいし、両方の機能を有していてもよい。

識別部５３は、画像取得部５１が取得した入力画像と、統計モデル取得部５２で取得した統計モデルである部分空間情報に基づいて、入力画像の断面種別を識別する。表示処理部５４は、入力画像と、該入力画像の断面種別の識別結果を、表示部３６の画像表示領域内に表示させる。表示部３６には、前記表示内容のほかにも、表示されている入力画像の縮尺情報や心電図の波形情報など、２次元超音波画像の観察に必要な種々の情報を表示することができる。

上記画像処理装置１０の各構成要素は、コンピュータプログラムに従って機能する。例えば、制御部３７（ＣＰＵ）がＲＡＭ３３を作業領域としてＲＯＭ３２又は記憶部３４などに記憶されたコンピュータプログラムを読み込み、実行することで、各構成要素の機能が実現される。なお、画像処理装置１０の構成要素の一部又はすべての機能が専用の回路を用いることで実現されてもよい。また、制御部３７の構成要素の一部の機能が、クラウドコンピューティングを用いることで実現されてもよい。例えば、画像処理装置１０とは異なる場所にある演算装置がネットワーク２１を介して画像処理装置１０と通信可能に接続され、画像処理装置１０と演算装置がデータの送受信を行うことで、画像処理装置１０又は制御部３７の構成要素の機能が実現されてもよい。

次に、図２のフローチャートを用いて、図１の画像処理装置１０の処理の例について説明する。

（ステップＳ１０１：統計モデルの取得）
ステップＳ１０１において、統計モデル取得部５２は、識別対象とする複数の断面種別にそれぞれ対応する複数の統計モデルを取得する。

図３Ａ～図３Ｄに、本実施形態において識別対象として例示する４種類の断面種別の断面画像を模式的に示す。図３Ａが心尖部二腔像、図３Ｂが心尖部三腔像、図３Ｃが心尖部四腔像、そして図３Ｄが空中放置像をそれぞれ表している。心尖部二腔像には、左心室３０３と左心房３０４の２つの領域、心尖部三腔像には、左心室３０５と左心房３０６と右心室３０７の３つの領域、心尖部四腔像には、左心室３１０と左心房３１１と右心室３０８と右心房３０９の４つの領域が夫々写っている。空中放置像（図３Ｄ）はプローブが被検体に接していない状態の画像であるため、いずれの領域も写っていない。

図４に、統計モデル取得部５２によって実行される、学習データから統計モデルを生成する処理（学習処理）の一例を示す。

ステップＳ４０１において、統計モデル取得部５２は、画像取得部５１を介してデータベース２２から学習データを取得する。そして、統計モデル取得部５２は、学習データと
して取得した複数の断面画像を、断面種別ごとにグループ分けする。本実施形態の場合は、心尖部二腔像・心尖部三腔像・心尖部四腔像・空中放置像という４つの画像群が生成される。以降のステップＳ４０２～Ｓ４０６の処理は、グループ分けされた画像群ごとに行われる。

ステップＳ４０２において、統計モデル取得部５２は、画像群に含まれる各画像に対して空間的正規化を施す。空間的正規化は、画像同士のスケールを揃えるための処理である。学習データには、異なる機種で撮像された画像や、異なる撮像パラメータ（撮像条件）で撮像された画像などが混在している。それゆえ、画像内の撮像領域（画像において信号が存在する領域（図３Ａ～図３Ｄの扇形の領域））の位置やサイズなどが、画像によって異なり得る。したがって、学習データとして集められた画像群をそのまま統計解析にかけると、撮像領域の位置やサイズのばらつきが統計モデルの表現力低下を招く可能性がある。そこで本実施形態では、統計解析に先立って空間的正規化を実行して、各画像における撮像領域の位置およびサイズを一致させる。

具体的には、まず統計モデル取得部５２は、学習データの画像の夫々について、プローブ位置座標と有効画像表示サイズを、空間的正規化に用いる情報として取得する。プローブ位置座標は、扇形の撮像領域の頂点（図３Ａの符号３０１）の画像座標であり、有効画像表示サイズは、扇形の撮像領域の半径（図３Ａの符号３０２）のピクセル数である。有効画像表示サイズは、超音波画像の視野深度（デプス）に対応している。プローブ位置座標および有効画像表示サイズの情報は、例えば、学習データ内に保持されている既知の情報を読み込むことで取得できる。この既知の情報は、学習データを準備する者が作成した情報であってもよいし、画像を撮像したモダリティ（超音波診断装置）から得られた情報であってもよい。あるいは、統計モデル取得部５２が、夫々の画像から画像処理で撮像領域（扇形の輪郭）を抽出することで、プローブ位置座標および有効画像表示サイズを取得してもよい。次に、統計モデル取得部５２は、各画像の有効画像表示サイズが所定のピクセル数となるように、画像全体を拡大又は縮小する。所定のピクセル数、すなわち正規化処理後の有効画像表示サイズは、あらかじめ定義されている（例えば５１２ピクセル）。その後、統計モデル取得部５２は、各画像のプローブ位置座標が全て一致するように画像を平行移動させる。

なお、空間的正規化処理は必須ではない。例えば、撮像領域の位置およびサイズのばらつきが十分に小さいことがあらかじめわかっている場合には、空間的正規化処理を省略してもよい。また、本実施形態では有効画像表示サイズを用いて画像のスケーリングを行う場合を例に挙げて説明したが、他の情報を用いてもよい。例えば、画像サイズ（画像全体のｘ方向およびｙ方向のピクセル数）が画像ごとに異なる場合に、画像サイズを一致させるように正規化してもよい。また、ピクセル数ではなく、１ピクセル当たりの物理サイズを考慮して、１ピクセル当たりの物理サイズが一致するように正規化してもよい。これによって、有効画像表示サイズの情報を取得できない場合であっても、正規化処理を実行することができる。

ステップＳ４０３において、統計モデル取得部５２は、空間的正規化を施した各画像のピクセル数を一致させるため、クロップ処理を行う。具体的には、統計モデル取得部５２は、扇形の撮像領域を包含する所定サイズの矩形領域を定義し、各画像から該矩形領域を切り出す。このとき、該矩形領域のｘ方向とｙ方向のピクセル数をそれぞれＮｘ、Ｎｙとする。なお、本実施形態では空間的正規化処理の後にクロップ処理を行ったが、最初に画像から撮像領域を切り出し、切り出された画像のサイズがＮｘ×Ｎｙとなるように拡大又は縮小することでも同様の結果を得ることができる。

ステップＳ４０４において、統計モデル取得部５２は、各画像の画素値情報を列ベクト
ルに変換する。クロップ処理後の画像はＮｘ×Ｎｙ個の画素で構成されている。ある画像Ｉの画素値情報の列ベクトルａは、画像左上を原点（０，０）として、画素位置（ｘ，ｙ）における画素値をＩ（ｘ，ｙ）で表すものとして、以下の通り定義される。

なお、処理対象の画像が１枚の２次元画像ではなく複数枚の２次元画像の場合には、それらの画素値情報を直列に並べたものを列ベクトルａとすればよい。

ステップＳ４０５において、統計モデル取得部５２は、画像群に含まれる全画像に対応する列ベクトル群を統計解析する。該統計解析には、主成分分析（ＰＣＡ：Ｐｒｉｎｃｉｐａｌ Ｃｏｍｐｏｎｅｎｔ Ａｎａｌｙｓｉｓ）などの既知の手法を用いればよく、Ｋｅｒｎｅｌ ＰＣＡやＷｅｉｇｈｔｅｄ ＰＣＡなどの手法を用いてもよい。また独立成分分析（ＩＣＡ：Ｉｎｄｅｐｅｎｄｅｎｔ Ｃｏｍｐｏｎｅｎｔ Ａｎａｌｙｓｉｓ）といった主成分分析以外の統計解析手法を用いてもよい。以下、ＰＣＡを用いる場合を例に説明する。

ＰＣＡを用いることで、列ベクトル群の平均ベクトルと固有ベクトル、そして各固有ベクトルに対応する固有値が計算される。これらを用いることで、同じ断面種別の画像群の特徴を表す部分空間を下記式のように表現できる。ｄは部分空間内に存在する点である。

ここで、ａバーが平均ベクトル、ｅ_ｉがｉ番目の基底における固有ベクトル、ｇ_ｉが固有ベクトルｅ_ｉに対応する係数である。また、Ｌは計算に用いる固有ベクトルの数を表している。

Ｌの値について述べる。Ｌの値は、学習データを必要十分に表現できる数（例えば、累積寄与率が９５％を超える数）よりも大幅に小さい値であることが望ましい。具体的には、学習データの累積寄与率が５０％以上７０％以下となるようにＬの値を設定するとよい。本発明者らの実験では、累積寄与率が５０％より小さい場合に、断面種別の識別率の低下が認められた。これは、累積寄与率が小さすぎる場合は、統計モデルの表現力が低すぎて、同じ断面種別内での画像の多様性に十分に対応できないからと考えられる。また、累積寄与率が７０％より大きい場合にも、断面種別の識別率の低下が認められた。これは、累積寄与率が大きすぎる場合は、統計モデルの表現力が高すぎ、全ての断面種別の統計モデルの当てはまりが良くなってしまうために、結果として識別力が低下するからと考えられる。

また、累積寄与率が５０％以上７０％以内に収まるＬの値の範囲が断面種別ごとに異なる場合であっても、Ｌの値を断面種別ごとに変化させるのではなく、同一の値に揃えることが望ましい。図５を例に挙げて説明する。図５は、固有ベクトルの数（Ｌの値）と累積寄与率の関係を示すグラフである。図５によると、空中放置像の方が他の種類（二腔像・三腔像・四腔像）に比べ、より小さい基底数で高い累積寄与率を示す。この例の場合、空中放置像のみＬの値を他の３種類と比べて小さく設定するのではなく、すべての断面種別においてＬ＝７～２１と設定することが望ましい。なぜなら、発明者らの実験によると、後段のステップＳ６０４において再構築誤差を算出する際、他の種別よりもＬの値が小さい種別は、目視で観察される違い以上に再構築誤差が大きく算出される場合があり、識別性能が低下する場合があるためである。このことから、本実施形態では、二腔像・三腔像
・四腔像における累積寄与率をもとに、Ｌ＝１０ないし２０という単一の値を採用する。

ステップＳ４０６において、統計モデル取得部５２は、統計解析の結果から、断面種別ごとに、統計モデルである部分空間情報（部分空間を構成する平均ベクトル（ａバー）とＬ個の固有ベクトル（ｅ_１～ｅ_Ｌ））を生成して、ＲＡＭ３３に格納する。以上で、統計モデルの学習処理は終了である。

なお、統計モデル取得部５２によるステップＳ４０１～Ｓ４０６の機能が本発明の「学習装置」の一例である。また、ステップＳ４０１の機能、ステップＳ４０２の機能、ステップＳ４０４～Ｓ４０６の機能がそれぞれ、本発明の「学習データ準備部」、「前処理部」、「モデル生成部」の一例である。

なお、ステップＳ４０１～Ｓ４０６における統計モデルの生成処理（学習処理）は、後述する入力画像の断面種別の識別処理とは独立したデータ処理である。そのため、統計モデルの生成処理を事前に実施し、生成された統計モデルをあらかじめデータベース２２又は記憶部３４に保存しておいてもよい。この場合、図２のステップＳ１０１では、統計モデル取得部５２は、統計モデルをデータベース２２又は記憶部３４から読み込んでＲＡＭ３３に格納すればよい。あらかじめ統計モデルを算出しておくことで、入力画像における断面種別の識別処理を行う際の処理時間を短縮できる効果がある。なお、統計モデルの生成処理は、画像処理装置１０の統計モデル取得部５２により行ってもよいし、画像処理装置１０とは異なる他の学習装置で行ってもよい。他の学習装置で統計モデルを生成する場合の処理も、図４に示したものと同様である。

なお、空間的正規化処理（ステップＳ４０２）とクロップ処理（ステップＳ４０３）の間に、学習データ内の画像の水増し処理を行ってもよい。すなわち、学習データ内の各画像を例えば各軸方向に±１０ピクセルずつ平行移動させ、そのように変位させた画像も学習データとして用いても良い。水増しは、平行移動だけでなく回転によって行ってもよい。また、両者を組み合わせてもよい。また、ステップＳ４０５で統計解析を行う際は、オリジナルの画像と水増しした画像をまとめて統計解析してもよいし、変位量（平行移動・回転量）が同一である画像ごとにグループ分けして、それぞれ統計解析を行ってもよい。変位量が同一である画像ごとにグループ分けして統計解析する場合、断面種別の数×変位量の組み合わせ数の統計モデルが生成されることになる。このようにすることで、入力画像の位置や角度のバリエーションに対する統計モデルの頑健性を向上させることができる。また、本実施形態では画素値情報として画素値そのものを利用したが、他の画像特徴（例えば、画像のテクスチャに関する特徴など）を画素値情報として利用してもよい。

（ステップＳ１０２：入力画像の取得）
ステップＳ１０２において、ユーザが操作部３５を介して画像の取得を指示すると、画像取得部５１は、ユーザが指定した入力画像をデータベース２２から取得し、ＲＡＭ３３に格納する。このとき、表示処理部５４は、入力画像を表示部３６の画像表示領域内に表示させてもよい。

（ステップＳ１０３：断面種別の識別）
ステップＳ１０３において、識別部５３は、ステップＳ１０１で取得した複数の統計モデルそれぞれの入力画像への当てはまりを評価することによって、入力画像の断面種別を識別する。本実施形態では、統計モデルとして、断面種別ごとの部分空間情報が用いられる。すなわち、識別すべき断面種別の数が４種類の場合、４つの部分空間情報が取得される。この場合、識別部５３は、断面種別ごとに入力画像の特徴が部分空間に属する可能性を表すスコアを算出し、断面種別ごとのスコアに基づいて入力画像の断面種別を識別するとよい。スコアの算出方法は問わないが、例えば、入力画像の特徴を部分空間へ投影した
のち元の画像空間に逆投影することで得られる再構築画像を断面種別ごとに生成し、各再構築画像と入力画像の類似度に基づいて断面種別ごとのスコアを算出するとよい。再構築画像と入力画像の類似度が高いほど、すなわち、再構築画像と入力画像の差異が小さいほど、入力画像がその部分空間に属する可能性が高いと評価することができる。したがって、断面種別ごとのスコアを比較することにより、入力画像が属する部分空間、すなわち、入力画像の断面種別を推定することができる。なお、スコアは、再構築画像と入力画像の類似度と正の相関をもつスコアでもよいし、負の相関をもつスコアでもよい。例えば、再構築画像と入力画像のあいだの差異（再構築誤差と呼ぶ）は、類似度と負の相関をもつスコアの一例である。

図６に、識別部５３によって実行される、断面種別の識別処理の具体例を示す。

ステップＳ６０１において、識別部５３は、画像取得部５１から入力画像の画素値情報を取得する。また、識別部５３は、入力画像の空間的正規化に用いる情報として、プローブ位置情報、有効画像表示サイズ情報を取得する。これらの情報は、例えば、超音波診断装置が入力画像を撮像する際に用いた撮像パラメータを入力画像に付帯情報として保持させておき、既知の情報として入力画像と共に読み込むことで取得できる。あるいは、入力画像が当該情報を保持していない場合には、入力画像から撮像領域を画像処理で抽出することで取得してもよい。

ステップＳ６０２において、識別部５３は、入力画像に対し、ステップＳ４０２、Ｓ４０３と同様の空間的正規化処理とクロップ処理を施すことで、入力画像のｘ方向とｙ方向のピクセル数を学習データと同じＮｘ・Ｎｙに揃える。そして、識別部５３は、式（１）に従い、入力画像を列ベクトルａ_ｔに変換する。

なお、入力画像に対する空間的正規化処理は、ステップＳ４０２において学習データに対して施した空間的正規化処理と同じ処理であることが望ましい。学習データと入力画像のスケールが揃っている方がより高い識別精度を期待できるからである。しかしながら、空間的正規化処理は必須ではなく、省略してもよい。また、有効画像表示サイズを基準とした正規化ではなく、画像サイズ（画像全体のｘ方向およびｙ方向のピクセル数）を基準とした正規化や、１ピクセル当たりの物理サイズを基準とした正規化を行ってもよい。なお、統計モデルである部分空間情報を算出するときは学習データに対し空間的正規化を行うが、入力画像に対しては空間的正規化を行わない、という構成にすることも可能である。この方法は、学習時はオフラインのマニュアル作業によって有効画像表示サイズを既知の情報として取得できるが、識別時は有効画像表示サイズが不明という場合に対応できる。あるいは、入力画像の有効画像表示サイズに合わせて学習データの空間的正規化を行うことで、入力画像の空間的正規化を省略してもよい。この方法は、入力画像の有効画像表示サイズが既知かつ固定の場合に特に有効である。

ステップＳ６０３において、識別部５３は、入力画像を夫々の部分空間で再構築する。いま、夫々の部分空間について、すべての固有ベクトルｅ_１～ｅ_Ｌを並べた行列Ｅを以下の式で定義する。

識別部５３は、以下の式を用いて入力画像ａ_ｔの固有空間への投影点ｐ_ｔを計算する。

ここでＥ^ＴはＥの転置行列を表す。次に、以下の式を用いて、固有空間への投影点ｐ_ｔを今度は逆投影する。

左辺（ａハット_ｔ）は、行列Ｅで表現される部分空間の表現力の範囲内で入力画像ａ_ｔを表現した、再構築画像である。識別部５３は、この再構築画像を、断面種別ごとのＥを用いてそれぞれ算出する。

ステップＳ６０４において、識別部５３は、入力画像と、断面種別ごとに算出された再構築画像との差異（すなわち、再構築誤差）を表すスコアを算出し、断面種別ごとのスコアを比較することで断面種別を識別する。本実施形態では、公知の画像類似度評価指標であるＳＳＤ（Ｓｑｕａｒｅｄ Ｓｕｍ Ｄｉｆｆｅｒｅｎｃｅ）を、入力画像の特徴を表す列ベクトルａ_ｔと、各再構築画像を表す列ベクトルａハット_ｔの差異を示すスコアとして用いる。そして、スコアが最も小さい再構築画像をもたらす部分空間と対応付いた断面種別を識別結果とする。なお、再構築誤差の算出には、上記したＳＳＤ以外にも、どのような画像類似評価尺度を用いてもよい。

なお、部分空間夫々に対して入力画像のスコアを算出するものであれば、上記以外の手法を用いてもよい。例えば、式（４）の投影処理のみを行い、最も平均ベクトルに近い、すなわち投影点ｐ_ｔのノルムが小さいものを識別結果としてもよい。あるいは、カーネル非線型部分空間法、相互部分空間法といった種々の公知の部分空間法を用いてもよい。

なお、識別部５３は、夫々の断面種別について算出されたスコアが所定の条件を満たす場合に、「入力画像はいずれの断面種別にも該当しない」という判定を行ってもよい。所定の条件は、例えば、所定の閾値より小さいスコアが存在しない（すなわち、どの断面種別についても再構築誤差が大きい）、断面種別間のスコアの差が所定値より小さい（すなわち、どちらの断面種別であるか明確に識別できない）、などが考えられる。これにより、例えば心臓領域以外を撮像した（あるいは、識別対象でない断面種別の）２次元超音波画像が入力された場合に、誤った識別結果が表示されてユーザに混乱を生じさせる可能性を低減させることができる。

ステップＳ６０５において、識別部５３は、ステップＳ６０４の識別結果をＲＡＭ３３に格納する。また、識別部５３は、データベース２２に識別結果を出力し、データベース２２が保持する入力画像の付帯情報として保存させるようにしてもよい。なお、１つの断面種別を特定することは行わずに、夫々の部分空間に対して算出されたスコアを識別結果としてＲＡＭ３３に格納する構成でもよい。なお、前記スコアは、再構築誤差の値そのものでもよいし、例えば０．０～１．０の範囲内に収まるように正規化した値でもよい。

（ステップＳ１０４：識別結果の表示）
ステップＳ１０４において、表示処理部５４は、入力画像とその種別識別結果を、それらが容易に視認できるような表示形態で、表示部３６の画像表示領域内に表示する。なお、断面種別の記録を目的とする場合には、本ステップの表示処理は必ずしも必要ではなく、識別した断面種別情報をデータベース２２又は記憶部３４に保存するだけの構成であってもよい。

なお、表示処理部５４は、断面種別だけでなく、ステップＳ１０３で算出されたスコアを表示してもよい。その場合、ユーザは、断面種別間のスコアの差を観察することで、どの程度識別結果が確からしいか、という判断をする材料を得ることができる。また、スコ
アだけを表示し、ユーザに断面種別を判断させる構成であってもよい。

なお、断面種別に応じた入力画像の解析や計測を目的とする場合には、本ステップのあと、画像処理装置１０は、断面種別に応じた処理を実行する。例えば、入力画像の断面種別を心尖部四腔像と判定した場合に左心室・左心房・右心室・右心房の抽出処理を行い、心尖部二腔像と判定した場合に左心室と左心房の抽出処理を行うような分岐処理を行う。この場合、断面種別の保存や表示は必ずしも行わなくてよい。

なお、本実施形態では４種類（心尖部二腔像・心尖部三腔像・心尖部四腔像・空中放置像）の断面種別の識別を行う場合を例に挙げて説明したが、識別する断面種別の数が異なる場合であっても同様に処理を行うことができる。例えば、心尖部四腔像・心尖部二腔像の２種類を識別してもよいし、短軸像や傍胸骨長軸像を加えた５種類以上の識別を行ってもよい。すなわち、断面種別は、心尖部二腔像、心尖部三腔像、心尖部四腔像、短軸像、傍胸骨長軸像、および、空中放置像のうちの、少なくとも２種類以上を含んでいればよい。

なお、画像取得部５１が取得する入力画像が動画像の各フレームの場合に、当該動画像の各フレームに対して、ステップＳ１０２、Ｓ１０３、Ｓ１０４の処理を逐次実施する構成であってもよい。

また、本実施形態では、入力画像と統計解析に用いる画像がともに２次元超音波画像である場合、すなわちモダリティが一致している場合を例に挙げて説明したが、両者のモダリティは異なっていてもよい。例えば、統計解析に用いる画像は３次元ＣＴ画像から切り出された２次元断面画像で入力画像は２次元超音波画像であってもよい。この場合、統計モデルで表現される画像と入力画像とでは、画素値分布や大小関係が異なる可能性がある。そのため、ステップＳ１０３では、正規化相関係数（ＮＣＣ：Ｎｏｒｍａｌｉｚｅｄ Ｃｏｒｒｅｌａｔｉｏｎ Ｃｏｅｆｆｉｃｉｅｎｔ）や相互情報量（ＭＩ：Ｍｕｔｕａｌ
Ｉｎｆｏｒｍａｔｉｏｎ）といった、画素値分布の違いに不偏な類似度評価尺度を用いて再構築誤差を算出することが望ましい。

本実施形態によれば、学習データから抽出された統計的傾向に基づいて、画像の断面種別の識別をより高い精度で実行できるという効果がある。

＜第２の実施形態＞
第１の実施形態では、部分空間を算出したり（ステップＳ１０１）、入力画像を投影・逆投影したり（ステップＳ１０３）する際、画像中の矩形領域を列ベクトルに変換していた。しかしながら、識別の対象がセクタ型プローブで撮像した２次元超音波画像の場合、その撮像領域は扇形であるため、範囲を矩形で定義すると、信号が存在しない撮像領域の外側を含むことになり、無駄が生じる。そこで、信号が存在する扇形の撮像領域のみを処理対象とするとよい。この方法によれば、統計モデルの表現能力を劣化させずに、統計モデルである部分空間情報の保持に必要なデータサイズを削減できる。本実施形態では、この方法を実現する具体的な方法を述べる。

図７は第２の実施形態に係る画像処理システムの構成例を示している。本実施形態に係る画像処理装置１０の構成は、第１の実施形態の構成（図１参照）と概ね同じである。ただし、機能的な構成として、制御部３７に計算対象範囲取得部６１が追加されており、統計モデル取得部５２および識別部５３にも機能が追加されている。

第２の実施形態において、計算対象範囲取得部６１は、データベース２２から取得した学習データに対して、部分空間の計算対象となる範囲（計算対象範囲）を取得する。詳細
は、ステップＳ８０１の説明で後述する。

統計モデル取得部５２および識別部５３は、計算対象範囲取得部６１で定義された計算対象範囲のみを対象として、部分空間の算出や識別処理を行う。

本実施形態においても第１の実施形態と同様に、入力画像および学習データが、心臓を撮像した２次元超音波画像である場合を例に挙げて説明する。そして、本実施形態では、入力画像として与えられた２次元超音波画像が心尖部二腔像・心尖部三腔像・心尖部四腔像・空中放置像の４種類のいずれかであるかを識別する。ただし、その他の種別の画像についても、同様の処理で識別可能である。

次に、図８を用いて、本実施形態の画像処理装置１０の処理の例について説明する。図８は第２の実施形態に係る画像処理装置１０の処理手順の例を示すフローチャートである。このフローチャートで示したステップにおいて、ステップＳ８０３およびステップＳ８０５は、図２に示した第１の実施形態におけるステップＳ１０２およびステップＳ１０４と同様の処理を行う。以下、追加された処理と第１の実施形態との相違部分を主に説明する。

（ステップＳ８０１：計算対象範囲の取得）
ステップＳ８０１において、計算対象範囲取得部６１は、学習データに含まれる夫々の画像の画素値情報（画像データ）をデータベース２２から取得する。そして、学習データ内の複数の画像の画素値情報を用いて、以降のステップで計算対象とする空間的範囲（計算対象範囲）を決定する。

計算対象範囲の決定方法を具体的に述べる。２次元超音波画像では、図３Ａ～図３Ｄで示しているとおり、矩形で表現された画像領域のうち、信号が存在するのは、扇形の撮像領域（以後、扇形領域と称す）だけである。そのため、以降のステップＳ８０２やＳ８０４を実行する際には、前記扇形領域のみを計算の対象とすればよいことが分かる。しかしながら、単純に、各画像における扇形領域をそれぞれ計算対象範囲とすることはできない。なぜなら、以降のステップＳ８０２やＳ８０４では、画像を列ベクトル化したときのベクトルの大きさ（式（１）のベクトルａの要素数）が、全ての学習データと入力画像の間で同一である必要があるからである。症例が異なったり、モダリティの機種・性能や撮像パラメータ（撮像条件）が異なったりすると、扇形領域の中心角の大きさが違うことがあるため、何らかの方法で、学習データと入力画像全てにおいて不変な計算対象範囲を定義する必要がある。

そこで本実施形態では、計算対象範囲取得部６１は、全ての学習データの（空間的正規化後の）撮像領域を全てカバーする扇形領域を算出して、これを計算対象範囲として取得する。ここで、全ての画像の撮像領域をカバーすることは必須ではなく、所定の割合（例えば９８％）の画像の撮像領域をカバーする領域を計算対象範囲としてもよい。なお、計算対象範囲の算出方法はこれに限らず、例えば、学習データの夫々から画像処理で心室・心房を抽出し、当該領域を全てカバーする扇形領域を計算対象範囲として取得してもよい。このとき、心室・心房の領域は、画像処理で算出した情報はなく、人手で与えた情報であってもよい。計算対象範囲取得部６１は、決定した計算対象範囲を以降のステップで使用するために、ＲＡＭ３３および記憶部３４に格納する。

なお、計算対象範囲は、操作者が手動で定義した範囲を用いるようにしてもよい。例えば、操作者が全学習データの心房・心室領域を包含する領域を目視で設定し、これを計算対象範囲としてデータベース２２又は記憶部３４に事前に保存しておいてもよい。この場合、計算対象範囲取得部６１は、データベース２２又は記憶部３４に保存されている計算
対象範囲を読み出すことで、計算対象範囲を取得する。

また、計算対象範囲は、信号が存在する扇形領域の全てを包含する領域でなくともよい。例えば、画素値が一定値以上の領域を包含する扇形領域を計算対象範囲とすることにより、扇形領域辺縁の画素値が低い（暗い）領域を除外しても良い。また、計算対象範囲は扇形に限らず、円形・三角形を含む任意の形状を用いることができる。

定義された計算対象範囲は、具体的には、下式で示す通り、ピクセルの番号と座標値の対応関係を表した対応表（ルックアップテーブル）として表現される。

ここで、例えばｍ［１］は計算対象範囲における１番目のピクセルを表す記号で、ｘ１・ｙ１はそれぞれ該ピクセルのｘ座標値・ｙ座標値である。そして、Ｋは計算対象範囲を構成するピクセルの数を表す。

なお、計算対象範囲の定義の保持方法は上記の例に限られず、例えば、計算対象範囲のピクセルに１、それ以外のピクセルに０が設定されたＮｘ×Ｎｙピクセルのマスク画像を利用してもよい。ただし、本実施形態のように計算対象範囲を連番（１～Ｋ）と座標値の対応表として保持する方法は、上で述べた、マスク画像を用いる一般的な方法に比べて次のような利点がある。マスク画像を用いる方法では、計算対象範囲に限定して何らかの処理を行う場合、Ｎｘ×Ｎｙピクセル全てに対して、該ピクセルが計算対象範囲に含まれているか否か検査する必要がある。すなわち、Ｎｘ×Ｎｙ個のピクセルを走査する必要がある。一方、本実施形態の保持方法であれば、計算対象範囲に含まれるＫ個（Ｎｘ×Ｎｙ個よりも少ない）のピクセルを走査すればよい。そのため、処理時間を短縮することが可能である。

（ステップＳ８０２：統計モデルの取得）
ステップＳ８０２において、統計モデル取得部５２は、学習データに含まれる夫々の画像の画素値情報（画像データ）と該画像の断面種別情報（心尖部二腔像・心尖部三腔像・心尖部四腔像・空中放置像のいずれか）をデータベース２２から取得する。また、統計モデル取得部５２は、ステップＳ８０１で決定した計算対象範囲をＲＡＭ３３又は記憶部３４から取得する。第１の実施形態におけるステップＳ１０１で述べた通り、前記画素値情報は、統計解析によって部分空間の情報を算出するために用いる学習データである。統計モデル取得部５２は、学習データを統計解析して、断面種別ごとに、統計モデルである部分空間情報を取得する。なお、学習処理が事前に行われている場合には、統計モデル取得部５２は、ステップＳ８０２において、データベース２２又は記憶部３４から、学習済みの統計モデルである部分空間情報を取得すればよい。

本ステップにおいて統計モデルである部分空間情報を取得するための方法は、第１の実施形態におけるステップＳ１０１と同様である。ただし、各画像の画素値情報を列ベクトルに変換する際、計算対象範囲内の画素値情報のみを用いる点で異なる。本実施形態においてある画像Ｉの画素値情報の列ベクトルａは、以下の式で定義される。

すなわち、第１の実施形態では列ベクトルａのサイズはＮｘ×Ｎｙであったのに対し、本実施形態ではＫ（Ｋ＜Ｎｘ×Ｎｙ）である。これ以降の処理は、第１の実施形態におけるステップＳ１０１と同一である。

このようにして、統計モデル取得部５２は、計算対象範囲に限定して統計解析を行う。そして、前記統計解析の結果から、断面種別ごとに、統計モデルである部分空間情報（部分空間を構成する固有ベクトルと平均ベクトルの情報）を取得して、ＲＡＭ３３に格納する。

（ステップＳ８０４：断面種別の識別）
ステップＳ８０４において、識別部５３は、ステップＳ８０３で取得された入力画像と、ステップＳ８０２で取得した統計モデルを取得する。また、識別部５３は、ステップＳ８０１で決定された計算対象範囲をＲＡＭ３３又は記憶部３４から読み込む。そして、識別部５３は、断面種別ごとの統計モデルの入力画像への当てはまりを評価することによって、入力画像の断面種別を識別する。

本ステップにおいて断面種別を識別するための方法は、第１の実施形態におけるステップＳ１０３と同じである。ただし、入力画像を列ベクトルａ_ｔに変換する際、計算対象範囲内の画素値情報のみを用いる点で異なる。これは、ステップＳ８０２における式（７）と同じである。したがって、入力画像と断面種別ごとに算出された再構築画像との間の差異（再構築誤差）を算出するときも、計算対象範囲に基づいて算出された、サイズＫの列ベクトルを用いる。すなわち、入力画像を表す列ベクトルａ_ｔと、各再構築画像を表す列ベクトルａハット_ｔの差異を求める。このようにすることで、サイズＮｘ×Ｎｙの列ベクトルを用いる第１の実施形態の方法よりも、計算コストを削減でき、処理を高速化することが可能である。

本実施形態によれば、断面種別識別に必要な統計モデル（部分空間情報）のデータサイズを削減し、識別処理に要する時間を短縮することができるという効果がある。

＜その他＞
以上、本発明の実施形態について説明したが、本発明はこれらに限定されるものではなく、請求項に記載された範囲内において、変更・変形することが可能である。

また、開示の技術は例えば、システム、装置、方法、プログラム若しくは記録媒体（記憶媒体）等としての実施態様をとることが可能である。具体的には、複数の機器（例えば、ホストコンピュータ、インターフェイス機器、撮像装置、ｗｅｂアプリケーション等）から構成されるシステムに適用してもよいし、また、１つの機器からなる装置に適用してもよい。

また、本発明の目的は、以下のようにすることによって達成されることはいうまでもない。すなわち、前述した実施形態の機能を実現するソフトウェアのプログラムコード（コンピュータプログラム）を記録した記録媒体（または記憶媒体）を、システムあるいは装置に供給する。係る記憶媒体は言うまでもなく、コンピュータ読み取り可能な記憶媒体である。そして、そのシステムあるいは装置のコンピュータ（またはＣＰＵやＭＰＵ）が記録媒体に格納されたプログラムコードを読み出し実行する。この場合、記録媒体から読み出されたプログラムコード自体が前述した実施形態の機能を実現することになり、そのプログラムコードを記録した記録媒体は本発明を構成することになる。

１０ 画像処理装置
３７ 制御部
５１ 画像取得部
５２ 統計モデル取得部
５３ 識別部

Claims (21)

1. 被検体の断面画像から当該断面画像の種別である断面種別を識別する画像処理装置であって、
   処理対象となる入力画像を取得する画像取得部と、
   断面種別が同じ断面画像を主成分分析することによって生成された、複数の断面種別に対応する複数の統計モデルを用いて、前記入力画像の断面種別を識別する識別部と、
   を有し、
   前記統計モデルは、断面種別が同じ複数の断面画像の特徴を表す部分空間の情報であることを特徴とする画像処理装置。
2. 前記識別部は、前記複数の統計モデルのいずれが前記入力画像に当てはまるかを評価することによって、前記入力画像の断面種別を識別する
   ことを特徴とする請求項１に記載の画像処理装置。
3. 断面種別が同じ断面画像を主成分分析することによって、複数の断面種別に対応する前記複数の統計モデルを生成するモデル生成部を有する
   ことを特徴とする請求項１又は２に記載の画像処理装置。
4. 前記部分空間の情報は、複数の固有ベクトルを含み、
   前記複数の固有ベクトルの数は、前記複数の断面画像の累積寄与率が５０％以上７０％以下となるように設定される
   ことを特徴とする請求項１～３のうちいずれか１項に記載の画像処理装置。
5. 前記識別部は、断面種別ごとの統計モデルを用いて、断面種別ごとに前記入力画像の特徴が部分空間に属する可能性を表すスコアを算出し、断面種別ごとのスコアに基づいて前記入力画像の断面種別を識別する
   ことを特徴とする請求項１～４のうちいずれか１項に記載の画像処理装置。
6. 前記識別部は、前記入力画像の特徴を部分空間へ投影したのち元の画像空間に逆投影することで得られる再構築画像を断面種別ごとに生成し、断面種別ごとの再構築画像と前記入力画像との類似度に基づいて断面種別ごとのスコアを算出する
   ことを特徴とする請求項５に記載の画像処理装置。
7. 前記識別部は、前記複数の断面種別のそれぞれについて算出されたスコアが所定の条件を満たす場合に、前記入力画像が前記複数の断面種別のいずれにも該当しないと判定することを特徴とする請求項５又は６に記載の画像処理装置。
8. 前記断面画像および前記入力画像は、心臓の２次元超音波画像である
   ことを特徴とする請求項１～７のうちいずれか１項に記載の画像処理装置。
9. 前記断面種別は、心尖部二腔像、心尖部三腔像、心尖部四腔像、短軸像、傍胸骨長軸像、および、プローブが被検体に接していない状態に取得される空中放置像のうちの、少なくとも２種類以上を含む
   ことを特徴とする請求項１～８のうちいずれか１項に記載の画像処理装置。
10. 前記統計モデルは、スケールを揃えるための空間的正規化が施された複数の断面画像を用いて生成される
    ことを特徴とする請求項１～９のうちいずれか１項に記載の画像処理装置。
11. 前記識別部は、前記入力画像に対し前記複数の断面画像と同じ空間的正規化を施したのちに、前記入力画像の断面種別の識別を行う
    ことを特徴とする請求項１０に記載の画像処理装置。
12. 前記断面画像および前記入力画像は、２次元超音波画像であり、
    前記空間的正規化は、２次元超音波画像の視野深度に対応する有効画像表示サイズが所定のピクセル数となるように画像を拡大又は縮小する処理を含む
    ことを特徴とする請求項１０又は１１に記載の画像処理装置。
13. 前記画像取得部は、超音波診断装置が前記入力画像を撮像した際の撮像パラメータの情報を前記入力画像とともに取得し、
    前記識別部は、前記撮像パラメータの情報に含まれる前記入力画像の視野深度に関する情報に基づいて前記空間的正規化を行う
    ことを特徴とする請求項１２に記載の画像処理装置。
14. 被検体の断面画像から当該断面画像の種別である断面種別を識別する医用画像診断装置であって、
    処理対象となる医用画像を取得する画像取得部と、
    断面種別が同じ断面画像を主成分分析することによって生成された、複数の断面種別に対応する複数の統計モデルを用いて、前記医用画像の断面種別を識別する識別部と、
    を有し、
    前記統計モデルは、断面種別が同じ複数の断面画像の特徴を表す部分空間の情報であることを特徴とする医用画像診断装置。
15. 被検体の断面画像から当該断面画像の種別である断面種別を識別する画像処理方法であって、
    処理対象となる入力画像を取得するステップと、
    断面種別が同じ断面画像を主成分分析することによって生成された、複数の断面種別に対応する複数の統計モデルを用いて、前記入力画像の断面種別を識別するステップと、
    を有し、
    前記統計モデルは、断面種別が同じ複数の断面画像の特徴を表す部分空間の情報であることを特徴とする画像処理方法。
16. 請求項１５に記載の画像処理方法の各ステップをコンピュータに実行させるためのプログラム。
17. 被検体の断面画像から当該断面画像の種別である断面種別を識別する処理に用いられるモデルを学習するための学習装置であって、
    複数の断面種別のそれぞれについて、複数の断面画像を学習データとして準備する学習データ準備部と、
    断面種別が同じ断面画像を主成分分析することによって、複数の断面種別に対応する複数の統計モデルを生成するモデル生成部と、
    を備え、
    前記統計モデルは、断面種別が同じ複数の断面画像の特徴を表す部分空間の情報であることを特徴とする学習装置。
18. 前記部分空間を表す情報は、複数の固有ベクトルを含み、
    前記モデル生成部は、前記複数の固有ベクトルの数を、累積寄与率に基づいて決定することを特徴とする請求項１７に記載の学習装置。
19. 前記モデル生成部は、前記複数の固有ベクトルの数を、前記複数の断面画像の累積寄与率が５０％以上７０％以下となるように決定する
    ことを特徴とする請求項１８に記載の学習装置。
20. 前記断面画像は、２次元超音波画像であり、
    前記学習装置は、前記複数の断面画像に対し、２次元超音波画像の視野深度に対応する有効画像表示サイズが所定のピクセル数となるように画像を拡大又は縮小する処理を含む空間的正規化を施す前処理部を有する
    ことを特徴とする請求項１７～１９のうちいずれか１項に記載の学習装置。
21. 前記断面画像は、心臓の２次元超音波画像であり、
    前記断面種別は、心尖部二腔像、心尖部三腔像、心尖部四腔像、短軸像、傍胸骨長軸像、および、プローブが被検体に接していない状態に取得される空中放置像のうちの、少なくとも２種類以上を含む
    ことを特徴とする請求項１７～２０のうちいずれか１項に記載の学習装置。
    

Images