(57)【要約】【課題】観察対象の動き伝搬が異常な領域を検出することが可能な画像処理装置、画像処理方法及び画像処理プログラムを提供する。【解決手段】本技術の一形態に係る画像処理装置は、動きベクトル算出部と、伝搬異常領域抽出部とを具備する。  上記動きベクトル算出部は、観察対象を経時的に撮像した対象画像から、上記対象画像上に設定された各部分領域に含まれる観察対象の動きベクトルを算出する。  上記伝搬異常領域抽出部は、上記動きベクトルに基づいて、上記観察対象の動き伝搬が異常な部分領域を伝搬異常領域として抽出する。【選択図】図１
2016041023.tif 000002
【特許請求の範囲】
【請求項１】
観察対象を経時的に撮像した対象画像から、前記対象画像上に設定された各部分領域に含まれる観察対象の動きベクトルを算出する動きベクトル算出部と、
前記動きベクトルに基づいて、前記観察対象の動き伝搬が異常な部分領域を伝搬異常領域として抽出する伝搬異常領域抽出部と
を具備する画像処理装置。
【請求項２】
請求項１に記載の画像処理装置であって、
前記伝搬異常領域抽出部は、前記動きベクトルから検出された動きの向きに基づいて、前記伝搬異常領域を抽出する
画像処理装置。
【請求項３】
請求項２に記載の画像処理装置であって、
前記伝搬異常領域抽出部は、前記動きベクトルから正常伝搬方向を検出し、前記正常伝搬方向とは異なる異常伝搬方向を示す部分領域を伝搬異常領域として抽出する
画像処理装置。
【請求項４】
請求項１に記載の画像処理装置であって、
前記観察対象は、心筋細胞であり、
前記動き伝搬は、拍動の伝搬である
画像処理装置。
【請求項５】
請求項４に記載の画像処理装置であって、
前記動きベクトルから検出された動きの速度に基づいて、速度変化特性の異常な部分領域を拍動異常領域として抽出する拍動異常領域抽出部
をさらに具備する画像処理装置。
【請求項６】
請求項５に記載の画像処理領域であって、
前記速度変化特性は、前記速度の周期的な変化特性である
画像処理装置。
【請求項７】
請求項１に記載の画像処理装置であって、
前記動きベクトル算出部は、前記対象画像を分割することで、前記部分領域を設定する
画像処理装置。
【請求項８】
請求項１に記載の画像処理装置であって、
前記伝搬異常領域抽出部は、前記動きベクトルから検出された動きの速度に基づいて、前記動き伝搬の速度が異常な部分領域を伝搬異常領域として抽出する
画像処理装置。
【請求項９】
動きベクトル算出部が、観察対象を経時的に撮像した対象画像から、前記対象画像上に設定された各部分領域に含まれる観察対象の動きベクトルを算出し、
伝搬異常領域抽出部が、前記動きベクトルに基づいて、前記観察対象の動き伝搬が異常な部分領域を伝搬異常領域として抽出する
画像処理方法。
【請求項１０】
観察対象を経時的に撮像した対象画像から、前記対象画像上に設定された各部分領域に含まれる観察対象の動きベクトルを算出する動きベクトル算出部と、
前記動きベクトルに基づいて、前記観察対象の動き伝搬が異常な部分領域を伝搬異常領域として抽出する伝搬異常領域抽出部と
として情報処理装置を動作させる画像処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本技術は、観察対象の動き伝搬が異常な領域を検出することが可能な画像処理装置、画像処理方法及び画像処理プログラムに関する。
【背景技術】
【０００２】
心筋細胞は、拍動の生成、伝達が可能な刺激伝導系と呼ばれる細胞群を含む。刺激伝導系では、自律的に拍動するペースメーカー細胞から下流の心筋細胞へ拍動が伝搬し、各細胞がリズミックに弛緩と収縮を行う。
【０００３】
このような心筋細胞を培養し、再生医療や創薬等の分野に役立てようという試みがされている。例えば、創薬の分野においては、培養心筋細胞に薬剤を投与し、拍動の伝搬等に変化が見られるか否かについて評価することが可能である。
【０００４】
そこで、例えば特許文献１には、心筋細胞等を撮像した画像を分割し、同一分割内で経時的な変化を評価する技術が記載されている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特表２０１０−５３８６０３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、特許文献１に記載の技術では、拍動の伝搬についての評価はできなかった。さらに、分割された領域間の比較はユーザによる作業に頼らざるを得なかったため、ある時点においてどの領域に異常な部位があるか、評価することが難しかった。
【０００７】
以上のような事情に鑑み、本技術の目的は、観察対象の動き伝搬が異常な領域を検出することが可能な画像処理装置、画像処理方法及び画像処理プログラムを提供することにある。
【課題を解決するための手段】
【０００８】
上記目的を達成するため、本技術の一形態に係る画像処理装置は、動きベクトル算出部と、伝搬異常領域抽出部とを具備する。
上記動きベクトル算出部は、観察対象を経時的に撮像した対象画像から、上記対象画像上に設定された各部分領域に含まれる観察対象の動きベクトルを算出する。
上記伝搬異常領域抽出部は、上記動きベクトルに基づいて、上記観察対象の動き伝搬が異常な部分領域を伝搬異常領域として抽出する。
【０００９】
この構成によれば、各部分領域に含まれる観察対象について算出された動きベクトルに基づいて、伝搬異常領域を抽出することができる。したがって、観察対象の動き伝搬が異常な領域を検出することが可能になる。さらに、部分領域間において動き伝搬や拍動プロファイル等を比較することが容易になり、対象画像全体の観察対象の動きからは検出することができない異常な領域を検出することができる。
【００１０】
また、上記伝搬異常領域抽出部は、上記動きベクトルから検出された動きの向きに基づいて、上記伝搬異常領域を抽出してもよい。
【００１１】
これにより、伝搬方向が異常となった伝播異常領域を抽出し、ユーザに通知することが可能となる。
【００１２】
例えば、上記伝搬異常領域抽出部は、上記動きベクトルから正常伝搬方向を検出し、上記正常伝搬方向とは異なる異常伝搬方向を示す部分領域を伝搬異常領域として抽出してもよい。
【００１３】
また、上記観察対象は、心筋細胞であり、
上記動き伝搬は、拍動の伝搬であってもよい。
【００１４】
あるいは、上記画像処理装置は、
上記動きベクトルから検出された動きの速度に基づいて、速度変化特性の異常な部分領域を拍動異常領域として抽出する拍動異常領域抽出部
をさらに具備してもよい。
【００１５】
これにより、速度変化特性によって示される心筋細胞の拍動プロファイルが異常な拍動異常領域を抽出し、ユーザに通知することができる。
【００１６】
上記速度変化特性は、上記速度の周期的な変化特性であってもよい。
【００１７】
これにより、心筋細胞の拍動プロファイルの周期性が異常な領域を抽出することが可能となる。
【００１８】
また、上記動きベクトル算出部は、上記対象画像を分割することで、上記部分領域を設定してもよい。
【００１９】
これにより、略同一の大きさの部分領域について評価することができるとともに、対象画像内を漏れなく評価することが可能となる。
【００２０】
また、上記伝搬異常領域抽出部は、上記動きベクトルから検出された動きの速度に基づいて、上記動き伝搬の速度が異常な部分領域を伝搬異常領域として抽出してもよい。
【００２１】
これにより、伝搬速度が変化した領域を伝搬異常領域として抽出し、ユーザに通知することができる。
【００２２】
本技術の一形態に係る画像処理方法は、動きベクトル算出部が、観察対象を経時的に撮像した対象画像から、上記対象画像上に設定された各部分領域に含まれる観察対象の動きベクトルを算出する。
伝搬異常領域抽出部が、上記動きベクトルに基づいて、上記観察対象の動き伝搬が異常な部分領域を伝搬異常領域として抽出する。
【００２３】
本技術の一形態に係る画像処理プログラムは、動きベクトル算出部と、伝搬異常領域抽出部と、として情報処理装置を動作させる。
上記上記動きベクトル算出部は、観察対象を経時的に撮像した対象画像から、上記対象画像上に設定された各部分領域に含まれる観察対象の動きベクトルを算出する。
上記伝搬異常領域抽出部は、上記動きベクトルに基づいて、上記観察対象の動き伝搬が異常な部分領域を伝搬異常領域として抽出する。
【発明の効果】
【００２４】
以上のように、本技術によれば、観察対象の動き伝搬が異常な領域を検出することが可能な画像処理装置、画像処理方法及び画像処理プログラムを提供することができる。
なお、ここに記載された効果は必ずしも限定されるものではなく、本開示中に記載されたいずれかの効果であってもよい。
【図面の簡単な説明】
【００２５】
【図１】本技術の一実施形態に係る画像処理装置を示すブロック図である。
【図２】上記画像処理装置の処理に供される対象画像の例を示す図である。
【図３】上記画像処理装置の動きベクトル算出部による部分領域毎の動きベクトルの算出例を示す図である。
【図４】上記画像処理装置の伝搬異常領域抽出部による処理について説明する図である。
【図５】上記画像処理装置の観察対象である心筋細胞の速度変化特性の検出例を示す図である。
【図６】上記画像処理装置の異常通知部による異常な領域の通知例を示す図である。
【図７】上記動きベクトル算出部による部分領域の設定例を示す図である。
【図８】図７Ａ，Ｂの部分領域各々に含まれる心筋細胞の動きベクトルの解析例を示すグラフである。
【図９】上記異常通知部よる異常な領域の通知例を示す図である。
【図１０】上記伝搬異常領域抽出部による正常伝搬方向及び異常伝搬方向の検出例を示す図である。
【図１１】上記伝搬異常領域抽出部による拍動伝搬の速度が異常な部分領域の抽出例を示す図である。
【図１２】Ａは、観察対象全体の動きベクトルの解析例を示すグラフであり、Ｂ〜Ｅは、対象画像中に設定された各部分領域の例と、各部分領域の動きベクトルの解析例を示すグラフである。
【図１３】上記画像処理装置の拍動異常領域抽出部の処理の詳細を説明する図である。
【図１４】Ａは、正常な拍動プロファイルに基づく速度変化特性を示すグラフであり、Ｂは、上記動きベクトル算出部による部分領域の動きベクトルの解析に基づく速度変化特性を示すグラフである。
【図１５】同画像処理装置の機能的構成を実現するためのハードウェア構成を示す模式図である。
【発明を実施するための形態】
【００２６】
以下、本技術に係る実施形態を、図面を参照しながら説明する。
【００２７】
［画像処理装置の構成］
図１は、本実施形態に係る画像処理装置１００を示すブロック図である。同図に示すように、画像処理装置１００は、画像取得部１０１と、動き検出部１０２と、動きベクトル算出部１０３と、伝搬異常領域抽出部１０４と、拍動異常領域抽出部１０５と、異常通知部１０６とを有する。
【００２８】
画像取得部１０１は、対象画像を取得する。対象画像は、観察対象を経時的に撮像した動画や連続的に撮像された複数の静止画であるものとすることができる。具体的には、対象画像は、明視野撮像、暗視野撮像、位相差撮像、蛍光撮像、共焦点撮像、多光子励起蛍光撮像、吸収光撮像、乱光撮像等の各種の光学的な画像撮像方法を利用して撮像された画像であるものとすることができる。
【００２９】
観察対象は、例えば心筋細胞とすることができる。あるいは、観察対象は、動きが伝搬する特性を有するものであればよく、心筋細胞の他、骨格筋細胞、平滑筋細胞、神経細胞やその他の細胞、生体組織、又は生物等とすることができる。
【００３０】
図２は、対象画像１０の例であり、複数の心筋細胞を含む画像を示す。画像取得部１０１は、図示しない撮像装置（顕微鏡撮像装置等）から対象画像１０を取得してもよく、ストレージ１３（図１５参照）に格納されている画像やネットワークから供給された画像を対象画像１０として取得するものとすることも可能である。画像取得部１０１は、取得した対象画像１０を動き検出部１０２へ供給する。
【００３１】
動き検出部１０２は、対象画像から観察対象の動きを検出する。動き検出部１０２は、対象画像全体から動きを検出することもできるし、後述する各部分領域毎に検出することもできる。また、観察対象の動きは、対象画像に対するブロックマッチング等の画像処理によって検出することができる。
【００３２】
動き検出部１０２は、具体的には、観察対象毎の動きベクトルを検出することができる。動きベクトルは、例えば、動きの向き及び大きさを表すものとすることができる。これにより、動き検出部１０２は、観察対象全体から検出された動きの向き及び大きさの代表値を算出し、当該代表値を観察対象全体の動きベクトルとして表すことができる。当該代表値は、動きの向き及び大きさの平均値とすることもできるし、最も多くの観察対象が示す動きの向き及び大きさの値等とすることもできる。
【００３３】
動きベクトル算出部１０３は、対象画像から、対象画像上に設定された各部分領域に含まれる観察対象の動きベクトルを算出する。動きベクトル算出部１０３は、具体的には、部分領域毎に、動き検出部１０２によって検出された各観察対象の動きの向き及び大きさの上記代表値を算出し、当該代表値を、各部分領域に含まれる観察対象の動きベクトルとして表すことができる。
【００３４】
動きベクトル算出部１０３が算出した動きベクトルの向きは、各部分領域での観察対象の動き伝搬の向きとみなすことができる。またこの動きベクトルの大きさは、各部分領域での観察対象の動き伝搬の大きさとみなすことができる。さらに、当該動きベクトルの大きさを、動きベクトルの検出元である対象画像の時間差分で除することにより、動きの速度を算出することができる。
【００３５】
図３は、動きベクトル算出部１０３による部分領域毎の動きベクトルの算出例を示す図である。動きベクトル算出部１０３は、対象画像を分割することで、部分領域を設定することができる。部分領域１１は、同図に示すように、対象画像１０を（ａ×ｂ）個に分割した領域とすることができる。これにより、略同一の大きさの部分領域について評価することができるとともに、対象画像１０内を漏れなく評価することが可能となる。
【００３６】
部分領域１１は、予め決められた所定の分割数で分割するものとすることもできるし、ユーザによって指定された分割数で分割するものとすることもできる。あるいは、動きベクトル算出部１０３は、動き検出部１０２により検出された各観察対象の動きに基づいて、各観察対象の動きベクトルの向きや大きさが異なる領域毎に部分領域を設定することができる。また、動きベクトル算出部１０３は、対象画像から検出された観察対象の密度に基づいて、部分領域を設定することもできる。
【００３７】
伝搬異常領域抽出部１０４は、動きベクトル算出部１０３により算出された動きベクトルに基づいて、観察対象の動き伝搬が異常な部分領域を伝搬異常領域として抽出する。動き伝搬は、一つの観察対象の動きが他の観察対象の動きに伝搬する現象であり、観察対象が心筋細胞の場合は、拍動の伝搬とすることができる。
【００３８】
ここで、心筋細胞は、生体内において常時収縮と弛緩を繰り返しながら拍動する。培養心筋細胞も、その全体が収縮と弛緩を繰り返すように、例えば各細胞が所定の方向に運動することで、拍動が伝搬する。心筋細胞に薬剤等を投与した場合、細胞の一部で他の細胞とは異なる方向に拍動の伝搬が行われたり、拍動の伝搬方向が旋回するパターンを示したり、拍動の伝搬速度が遅くなったりといった変化が生じる可能性がある。本実施形態によれば、このような拍動の伝搬の異常な領域を検出することが可能となる。
【００３９】
伝搬異常領域抽出部１０４は、本実施形態において、動きベクトルから検出された動きの向きに基づいて、伝搬異常領域を抽出するものとすることができる。より具体的には、伝搬異常領域抽出部１０４は、部分領域毎に算出された動きベクトルから正常伝搬方向を検出し、正常伝搬方向とは異なる異常伝搬方向を示す部分領域を伝搬異常領域として抽出するものとすることができる。正常伝搬方向は、例えば、検出された動きベクトルの向きのうち、最も多数の観察対象から検出される向きとすることができる。あるいは、正常伝搬方向は、観察対象の種類や配置から予め予想される動きベクトルの向き等とすることもできる。
【００４０】
図４は、伝搬異常領域抽出部１０４による処理について説明する図であり、図４Ａは拍動伝搬の向きの検出例を示し、図４Ｂは伝搬異常領域の抽出例を示す。同図Ａに示すように、黒矢印で示す正常伝搬方向は、一方向を示すのに対し、白矢印で示す異常伝搬方向は、様々な向きを示す。したがって、同図Ｂに示すように、伝搬異常領域抽出部１０４は、異常伝搬方向が検出された領域を伝搬異常領域１２として抽出することができる。
【００４１】
拍動異常領域抽出部１０５は、動きベクトルから検出された動きの速度に基づいて、速度変化特性の異常な部分領域を拍動異常領域として抽出する。速度変化特性は、動きの速度の周期的な変化特性であるものとすることができる。また、例えば心筋細胞は、周期的に速度のピークを有する速度変化特性を示す。すなわち心筋細胞における速度変化特性は、心筋細胞の拍動プロファイルを示すものとすることができる。
【００４２】
図５は、観察対象である心筋細胞の速度変化特性の検出例を示す図であり、いずれも縦軸は動きの速度、横軸は時間を示す。また同図Ａは図３に示す部分領域１１Ａの結果を示し、同図Ｂは図３に示す部分領域１１Ｂの結果を示し、同図Ｃは対象画像１０全体の代表値の結果を示す。
【００４３】
図５Ｂ及び図５Ｃでは、２つの速度ピークＰ１，Ｐ２が検出されているのに対し、同図Ａでは、２つの速度ピークＰ１，Ｐ２の他、小さな速度ピークＰ３が検出されている。したがって、拍動異常領域抽出部１０５は、同図Ａの検出された部分領域１１Ａにおける速度変化特性が異常であると判断し、部分領域１１Ａを拍動異常領域として抽出することができる。
【００４４】
なお、拍動異常領域抽出部１０５は、伝搬異常領域抽出部１０４の伝搬異常領域の抽出と並行して拍動異常領域の抽出を行うこともできるし、伝搬異常領域の抽出の前後に拍動異常領域の抽出を行うこともできる。
【００４５】
異常通知部１０６は、伝搬異常領域抽出部１０４及び拍動異常領域抽出部１０５によって抽出された異常な領域を通知する。当該異常な領域は、伝搬異常領域、又は伝搬異常領域及び拍動異常領域を含む領域であるものとすることができる。
【００４６】
図６は、異常な領域１３Ａ，１３Ｂの通知例を示す図である。同図に示すように、異常な領域１３Ａ，１３Ｂは、実際に抽出された伝搬異常領域及び拍動異常領域よりも大きい領域として示してもよい。なお、異常な領域を示す方法は同図に示す方法に限定されない。例えば、異常通知部１０６は、伝搬異常領域及び拍動異常領域の中心上に表示されたマーク等によって異常な領域を示すこともできるし、図４Ｂに示すように網掛けや色を変化させることにより異常な領域を示すこともできる。
【００４７】
画像処理装置１００は、以上のように構成されている。画像処理装置１００は、上述のように、部分領域毎に算出される動きベクトルに基づいて、伝搬異常領域や拍動異常領域を抽出することが可能である。したがって、観察対象の動き伝搬が異常な領域を検出することが可能となる。
【００４８】
また、画像処理装置１００によれば、部分領域間において動き伝搬や拍動プロファイル等を比較することが容易になり、対象画像全体の観察対象の動きからは検出することができない異常な領域を検出することが可能となる。すなわち、従来ユーザの作業に頼らざるを得なかった部分領域間の細胞の反応性や挙動等の評価を容易に行うことが可能となる。
【００４９】
さらに、画像処理装置１００によれば、対象画像の全体にわたり異常をきたした心筋細胞等が存在する領域について自動で解析を行い、評価することができる。これにより、ユーザの作業による見落としを防止することができる。
【００５０】
以上のように、画像処理装置１００は、心筋細胞のように領域依存的に異常が生じ、不整脈の原因となる心筋梗塞等が発生する観察対象に対して、特に有効に用いることができる。以下、培養心筋細胞に対する動作例を示し、本実施形態についてより詳細に説明する。
【００５１】
［動作例１］
動作例１として、培養心筋細胞に所定の薬剤を投与し、その薬剤の効果を各部分領域の動きベクトルに基づいて評価した。以下、図１に沿って説明する。
【００５２】
まず、画像取得部１０１が、対象画像を取得する。ここでは、対象画像は、薬剤が投与された培養心筋細胞を位相差撮像した画像であるものとすることができる。
【００５３】
続いて、動き検出部１０２が、対象画像全体の各観察対象の動きベクトルを検出する。
【００５４】
続いて、動きベクトル算出部１０３が、部分領域を設定する。
【００５５】
図７Ａ，Ｂは、動きベクトル算出部１０３による部分領域の設定例を示す図である。部分領域１１Ａ，１１Ｂは、例えば、ユーザの入力操作等により指定された大きさ、位置に設定されてもよい。またこれらの図に示すように、同一対象画像１０内の部分領域１１Ａ，１１Ｂの大きさは異なっていてもよい。
【００５６】
続いて、動きベクトル算出部１０３が、各部分領域における動きベクトルや動きの速度等を算出する。
【００５７】
図８Ａ，Ｂは、図７Ａ，Ｂの部分領域１１Ａ，１１Ｂのそれぞれに含まれる心筋細胞の動きベクトルの解析例を示すグラフであり、いずれも縦軸は動きベクトルの大きさから算出された動きの速度、横軸は時間を示す。また、グラフ中の実線の矢印は、算出された動きベクトルの向きを模式的に示す。図８Ａに示すように、部分領域１１Ａに含まれる観察対象の動きベクトルの向きは、いずれも、同一の向きである。一方、図８Ｂに示すように、部分領域１１Ｂに含まれる観察対象の動きベクトルの向きは、部分領域１１Ａの動きベクトルと異なる向きのものを含む。
【００５８】
続いて、伝搬異常領域抽出部１０４が、部分領域１１Ａの動きベクトルの向き（図８Ａ参照）を正常伝搬方向と決定し、部分領域１１Ｂの動きベクトルの向き（図８Ｂ参照）をいずれも異常伝搬方向と決定する。これにより、伝搬異常領域抽出部１０４が、部分領域１１Ｂを、伝搬異常領域として抽出する。
【００５９】
さらに、異常通知部１０６が、図９に示すように、抽出した伝搬異常領域を通知する。図９では、対象画像１０上に、黒矢印で正常伝搬方向を示し、白矢印で異常伝搬方向を示し、かつ異常伝搬方向が検出された伝搬異常領域を含む異常な領域１３を矩形の枠で示している。
【００６０】
また、上記の解析によって、部分領域１１Ｂにおいて異常伝搬方向の動きが生じることで、図８Ａ，Ｂの破線矢印に示す異常な速度変化が生じている可能性があることがわかる。これにより、対象画像１０内の領域全体で異常な動きを示している可能性についても検討することができる。
【００６１】
以上の動作により、ユーザが、伝搬異常領域を容易に判別し、薬剤の効果や異常の生じた心筋細胞の解析を詳細かつ効率的に進めることが可能となる。また以下に、拍動伝搬の解析方法の他の例を説明する。
【００６２】
［動作例２］
動作例２として、動作例１と同様の培養心筋細胞に所定の薬剤を投与した例を示す。なお、本動作例において、画像取得部１０１及び動き検出部１０２の動作は動作例１と同様であるため、その説明を省略する。
【００６３】
動きベクトル算出部１０３が、対象画像を（ａ×ｂ）個に分割し、部分領域を設定する（図３参照）。そして動きベクトル算出部１０３が、各部分領域１１についての動きベクトルを算出する。
【００６４】
続いて、伝搬異常領域抽出部１０４が、算出された各部分領域１１の動きベクトルの向き及び大きさに基づいて、正常伝搬方向及び異常伝搬方向を検出する。
【００６５】
図１０Ａ，Ｂは、伝搬異常領域抽出部１０４による正常伝搬方向及び異常伝搬方向の検出例を示す図である。図１０Ａに示す例では、対象画像１０の右下部に正常伝搬方向（黒矢印）が検出されているのに対し、対象画像１０の中央部から左上部にかけて、正常伝搬方向とは異なる方向の異常伝搬方向（白矢印）が検出されている。したがって、伝搬異常領域検出部１０４は、対象画像１０の中央部及び左上部の部分領域を、伝搬異常領域１２と抽出することができる。
【００６６】
また図１０Ｂに示す例では、対象画像１０の下部に正常伝搬方向（黒矢印）が検出されているが、そこから異常伝搬方向（白矢印）が旋回するパターンが検出されている。なお、同図において、説明のため旋回方向をグレーの矢印で示している。したがって、伝搬異常領域検出部１０４は、旋回パターンを示す部分領域を、伝搬異常領域１２と抽出することができる。
【００６７】
これにより、異常通知部１０６が、伝搬異常領域抽出部１０４によって抽出された伝搬異常領域１２に基づいて異常な領域を通知することができる。
【００６８】
［動作例３］
動作例３として、伝搬異常領域抽出部１０４が、拍動伝搬速度が変化した領域を伝搬異常領域として抽出する例を示す。本動作例において、伝搬異常領域抽出部１０４以外の各要素の動作は動作例２と同様であるため、その説明を省略する。
【００６９】
伝搬異常領域抽出部１０４は、例えば、動きベクトルから検出された動きの速度に基づいて、動き伝搬（拍動伝搬）の速度が異常な部分領域を伝搬異常領域として抽出することができる。
【００７０】
図１１は、伝搬異常領域抽出部１０４による拍動伝搬の速度が異常な部分領域の抽出例を示す図である。同図に示すように、対象画像１０全体において拍動伝搬の向きは同一であるが、拍動伝搬の速度は徐々に遅くなっており、特に、対象画像１０の左上部において、速度の変化が顕著に見られる。このため、伝搬異常領域抽出部１０４は、速度の変化が顕著に見られる部分領域を、伝搬異常領域１２として抽出することができる。
【００７１】
拍動伝搬の速度の異常の判定基準は特に限定されないが、例えば、所定の閾値を設定し、その速度以下の場合は速度が異常であると判定するものとすることができる。
【００７２】
以上の動作例のように、拍動伝搬の速さが異常な領域を検出することでも、薬剤の効果や異常の生じた心筋細胞の解析の効率化に寄与することが可能となる。
【００７３】
［動作例４］
動作例４として、伝搬異常領域抽出部１０４による処理と並行して、又はその前後に、拍動異常領域抽出部１０５が速度変化特性の異常な領域を抽出する例を示す。なお、本動作例において、画像取得部１０１、動き検出部１０２及び異常通知部１０６の動作は動作例１と同様であるため、その説明を省略する。
【００７４】
動きベクトル算出部１０３は、部分領域を設定し、各部分領域についての動きベクトルを算出する。
【００７５】
図１２Ａは、観察対象全体の動きベクトルの解析例を示すグラフでであり、図１２Ｂ〜Ｅは、対象画像１０中に設定された部分領域１１Ｂ，１１Ｃ，１１Ｄ，１１Ｅの例と、各部分領域１１Ｂ〜１１Ｅの動きベクトルの解析例を示すグラフである。また、いずれのグラフも縦軸は動きベクトルの大きさから算出された動きの速度、横軸は時間を示す。上述のように、同図のグラフに示す速度変化特性は、心筋細胞の拍動プロファイルを示す。
【００７６】
図１２Ａに示す観察対象全体の動きベクトルの解析例からは、正常な拍動プロファイルに見られる速度ピーク以外に、目立って異常な速度ピークは検出されないため、全体として正常な拍動プロファイルの対象画像１０と判断される可能性がある。一方、図１２Ｂ及び図１２Ｄに示す部分領域１１Ｂ，１１Ｄの動きベクトルの解析例からは、図１２Ａではほとんど検出されなかった異常な速度ピークＰＢ，ＰＤがそれぞれ検出されている。
【００７７】
そこで、拍動異常領域抽出部１０５が、正常な拍動プロファイルに基づく速度変化特性においては検出されない異常な速度ピークが検出した部分領域１１Ｂ，１１Ｄを、拍動異常領域として抽出する。以下、当該処理についてより詳細に説明する。
【００７８】
図１３は、拍動異常領域抽出部１０５の処理の詳細を説明する図であり、Ａは正常な拍動プロファイルに基づく速度変化特性を示すグラフ、Ｂは、動きベクトル算出部１０３による部分領域の動きベクトルの解析に基づく速度変化特性を示すグラフである。また、同図Ａ，Ｂに示す符号Ｐ１〜Ｐ７は、それぞれ動きの速度の主なピークとボトムを示す。
【００７９】
図１３Ａに示すように、正常な拍動プロファイルは、Ｐ１〜Ｐ５の５点で一周期（一拍）を構成する。一方、早期後脱分極（early afterdepolarization, ＥＡＤ）や遅延後脱分極 （delayed afterdepolarization, ＤＡＤ）等の不整脈が生じると、図１３Ｂに示すように、Ｐ５の後にＰ６，Ｐ７等に示される異常な拍動が見られる。
【００８０】
このため、例えば拍動異常領域抽出部１０５は、図１３Ｂに示すように、正常な拍動プロファイルに見られる拍動Ｐ１〜Ｐ５の他Ｐ６，Ｐ７のような拍動が見られ、かつＰ７の後にピークの見られない時間が所定以上ある場合、Ｐ６，Ｐ７を異常な速度変化であると判定する。そして拍動異常領域抽出部１０５は、このような異常な速度変化が検出された部分領域を拍動異常領域と抽出することができる。これにより、精度よく異常な拍動プロファイルを判定することができる。
【００８１】
以上のように、本動作例によれば、部分領域毎に心筋細胞の速度変化特性を確認することができ、拍動プロファイルが異常な領域を容易に検出することができる。
【００８２】
［動作例５］
動作例５として、拍動異常領域抽出部１０５の他の動作例について説明する。
【００８３】
図１４Ａは、正常な拍動プロファイルに基づく速度変化特性を示すグラフであり、図１４Ｂは、動きベクトル算出部１０３による部分領域の動きベクトルの解析に基づく速度変化特性を示すグラフである。また、同図に示す符号ａ，ｂ，ｃ，ｄ及びａ'，ｂ'，ｃ'，ｄ'は、動きの速度の主なピークとボトムの時間を示す。
【００８４】
拍動異常領域抽出部１０５は、例えば、図１４Ａのａ−ｄ間に示される正常な拍動プロファイルの一周期の時間ｔ１に対し、図１４Ｂのａ'−ｄ' 間の時間ｔ'１が変化していた場合、そのような時間変化が検出された部分領域１１を拍動異常領域と抽出することができる。図１４Ｂに示す例では、図１４Ａに示す時間ｔ１よりも、時間ｔ'１が延長している。したがって、拍動異常領域抽出部１０５は、この部分領域を拍動異常領域と抽出することができる。なお、時間ｔ１の具体的な設定値や、時間ｔ'１が時間ｔ１に対して変化しているか否かの判断基準は、一般的な拍動プロファイルに鑑み適宜設定することができる。
【００８５】
あるいは、拍動異常領域抽出部１０５は、ｂ−ｃ間に示される正常な拍動プロファイルの一部の時間ｔ２に対し、ｂ'−ｃ' 間の時間ｔ'２が変化していた場合、その部分領域１１を拍動異常領域と抽出することができる。例えばｂ−ｃ間は、動きの速度ピーク間の時間とすることができる。図１４Ｂに示す例では、ｂ−ｃ間の時間よりも、ｂ'−ｃ' 間の時間が延長している。したがって、やはりこの部分領域を拍動異常領域と抽出することができる。
【００８６】
［動作例７］
他の動作例として、例えば、動きベクトル算出部１０３が部分領域を設定し、各部分領域毎の動きベクトルを算出した後、設定された部分領域の一部を新たな部分領域として設定してもよい。この場合、初めに設定された部分領域について伝搬異常領域抽出部１０４や拍動異常領域抽出部１０５によって異常な領域が抽出されなかった場合、動きベクトル算出部１０３が新たな部分領域を設定し、再度処理を繰り返してもよい。あるいは、初めに設定された部分領域について異常通知部１０６が結果を通知した後、ユーザの操作等により動きベクトル算出部１０３が新たな部分領域を設定し、再度処理を繰り返してもよい。
【００８７】
[画像処理装置のハードウェア構成]
上記のような画像処理装置１００の機能的構成は、以下に示すハードウェア構成によって実現することが可能である。
【００８８】
図１５は、画像処理装置１００のハードウェア構成を示す模式図である。同図に示すように画像処理装置１００はハードウェア構成として、ＣＰＵ１１、メモリ１２、ストレージ１３及び入出力部（Ｉ／Ｏ）１４を有する。これらはバス１５によって互いに接続されている。
【００８９】
ＣＰＵ（Central Processing Unit）１１は、メモリ１２に格納されたプログラムに従って他の構成を制御すると共に、プログラムに従ってデータ処理を行い、処理結果をメモリ１０２に格納する。ＣＰＵ１１はマイクロプロセッサであるものとすることができる。
【００９０】
メモリ１２はＣＰＵ１１によって実行されるプログラム及びデータを格納する。メモリ１０２はＲＡＭ（Random Access Memory）であるものとすることができる。
【００９１】
ストレージ１３は、プログラムやデータを格納する。ストレージ１３はＨＤＤ（Hard disk drive）やＳＳＤ（solid state drive）であるものとすることができる。
【００９２】
入出力部１４は画像処理装置１００に対する入力を受け付け、また画像処理装置１００の出力を外部に供給する。入出力部１４は、キーボードやマウス等の入力機器やディスプレイ等の出力機器、ネットワーク等の接続インターフェイスを含む。
【００９３】
画像処理装置１００のハードウェア構成はここに示すものに限られず、画像処理装置１００の機能的構成を実現できるものであればよい。また、上記ハードウェア構成の一部又は全部はネットワーク上に存在していてもよい。
【００９４】
以上、本技術の実施形態について説明したが、本技術は上述の実施形態にのみ限定されるものではなく、本開示の要旨を逸脱しない範囲内において種々変更され得る。
【００９５】
例えば、画像処理装置１００は、画像処理システムの一部を構成し、少なくとも、動きベクトル算出部１０３と、伝搬異常領域抽出部１０４とを有するものとすることができる。この場合、画像取得部１０１と、動き検出部１０２と、異常通知部１０６とは、画像処理装置１００と有線又は無線により接続された他の装置等に含まれるものとすることができる。また、上述の実施形態で示したように、画像処理装置１００は、拍動異常領域抽出部１０５を有しないものとすることもできる。
【００９６】
なお、本技術は以下のような構成もとることができる。
（１）観察対象を経時的に撮像した対象画像から、上記対象画像上に設定された各部分領域に含まれる観察対象の動きベクトルを算出する動きベクトル算出部と、
上記動きベクトルに基づいて、上記観察対象の動き伝搬が異常な部分領域を伝搬異常領域として抽出する伝搬異常領域抽出部と
を具備する画像処理装置。
（２）上記（１）に記載の画像処理装置であって、
上記伝搬異常領域抽出部は、上記動きベクトルから検出された動きの向きに基づいて、上記伝搬異常領域を抽出する
画像処理装置。
（３）上記（２）に記載の画像処理装置であって、
上記伝搬異常領域抽出部は、上記動きベクトルから正常伝搬方向を検出し、上記正常伝搬方向とは異なる異常伝搬方向を示す部分領域を伝搬異常領域として抽出する
画像処理装置。
（４）上記（１）から（３）のうちいずれか１つに記載の画像処理装置であって、
上記観察対象は、心筋細胞であり、
上記動き伝搬は、拍動の伝搬である
画像処理装置。
（５）上記（４）に記載の画像処理装置であって、
上記動きベクトルから検出された動きの速度に基づいて、速度変化特性の異常な部分領域を拍動異常領域として抽出する拍動異常領域抽出部
をさらに具備する画像処理装置。
（６）上記（５）に記載の画像処理装置であって、
上記速度変化特性は、上記速度の周期的な変化特性である
画像処理装置。
（７）上記（１）から（６）のうちいずれか１つに記載の画像処理装置であって、
上記動きベクトル算出部は、上記対象画像を分割することで、上記部分領域を設定する
画像処理装置。
（８）上記（１）から（７）のうちいずれか１つに記載の画像処理装置であって、
上記伝搬異常領域抽出部は、上記動きベクトルから検出された動きの速度に基づいて、上記動き伝搬の速度が異常な部分領域を伝搬異常領域として抽出する
画像処理装置。
（９）動きベクトル算出部が、観察対象を経時的に撮像した対象画像から、上記対象画像上に設定された各部分領域に含まれる観察対象の動きベクトルを算出し、
伝搬異常領域抽出部が、上記動きベクトルに基づいて、上記観察対象の動き伝搬が異常な部分領域を伝搬異常領域として抽出する
画像処理方法。
（１０）観察対象を経時的に撮像した対象画像から、上記対象画像上に設定された各部分領域に含まれる観察対象の動きベクトルを算出する動きベクトル算出部と、
上記動きベクトルに基づいて、上記観察対象の動き伝搬が異常な部分領域を伝搬異常領域として抽出する伝搬異常領域抽出部と
として情報処理装置を動作させる画像処理プログラム。
【符号の説明】
【００９７】
１０…対象画像
１１…部分領域
１２…伝搬異常領域
１００…画像処理装置
１０３…動きベクトル算出部
１０４…伝搬異常領域抽出部
１０５…拍動異常領域抽出部
【図１】
2016041023.tif 000003
【図３】
2016041023.tif 000004
【図５】
2016041023.tif 000005
【図６】
2016041023.tif 000006
【図１３】
2016041023.tif 000007
【図１４】
2016041023.tif 000008
【図１５】
2016041023.tif 000009
【図２】
2016041023.tif 000010
【図４】
2016041023.tif 000011
【図７】
2016041023.tif 000012
【図８】
2016041023.tif 000013
【図９】
2016041023.tif 000014
【図１０】
2016041023.tif 000015
【図１１】
2016041023.tif 000016
【図１２】
2016041023.tif 000017