(57)【要約】【課題】画像のノイズを低減させることができる医用画像処理装置、医用画像診断装置及び画像処理方法を提供すること。【解決手段】本発明の実施形態は、医用画像処理装置、方法及び医用画像診断装置に関する。医用画像処理装置は、変形場計算部と、変形場重み付け部及び画像変形部を含む。変形計算部は、対象物の第１画像の第２画像に対する変形場を計算する。変形場重み付け部は、対象物の各部位の動き量に応じて変形場における重み付けを行う。画像変形部は、前記重み付けが行われた変形場を用いて前記第１画像を変形する。【選択図】図１Ａ
2014174174.tif 000002
【特許請求の範囲】
【請求項１】
対象物の第１画像の第２画像に対する変形場を計算する変形場計算部と、
前記対象物の各部位の動き量に応じて前記変形場における重み付けを行う変形場重み付け部と、
前記重み付けが行われた変形場を用いて前記第１画像を変形する画像変形部と
を備えたことを特徴とする医用画像処理装置。
【請求項２】
前記変形場重み付け部は、前記第２画像と複数の前記第１画像の前記対象物の各部位の動き量が反映される分散画像に基づいて、前記変形場における重み付けを行うことを特徴とする請求項１に記載の医用画像処理装置。
【請求項３】
前記変形場重み付け部は、
アトラスに基づいて、前記第１画像における前記対象物の各部位を識別するアトラス識別部と、
識別された各部位の動き量の空間分布に基づいて、前記変形場における重み付けを行う変形場重み付け部と
を有することを特徴とする請求項１に記載の医用画像処理装置。
【請求項４】
前記変形場重み付け部は、前記対象物の運動部位の変形場の重みを１とし、非運動部位の変形場の重みを１より小さくすることを特徴とする請求項１、２又は３に記載の医用画像処理装置。
【請求項５】
前記変形場重み付け部は、前記対象物の運動部位と非運動部位との間の移行部位の変形場の重みを、前記運動部位の変形場の重みと前記非運動部位の変形場の重みとの間の値とすることを特徴とする請求項４に記載の医用画像処理装置。
【請求項６】
前記第１画像のヒストグラムと前記第２画像のヒストグラムとの間でヒストグラムマッチングを行い、前記第２画像とグレースケール範囲が近い第１画像を取得するヒストグラムマッチング部をさらに備え、
前記変形場計算部は、前記第２画像とグレースケール範囲が近い第１画像の前記第２画像に対する変形場を計算することを特徴とする請求項１〜５のいずれか一つに記載の医用画像処理装置。
【請求項７】
前記変形場計算部は、
前記第１画像の前記第２画像に対する力場を計算する力場計算部と、
前記力場に基づいて変形場を計算する変形場計算部と、
前記変形場を用いて前記第１画像を変形する画像変形部と、
前記力場計算部、前記変形場計算部及び前記画像変形部を制御し、変形前の前記第１画像に替えて変形後の前記第１画像を用いて、所定の条件を満足するまで処理を繰り返し行って、最終的な変形場を取得する制御部と
を有することを特徴とする請求項１〜６のいずれか一つに記載の医用画像処理装置。
【請求項８】
計算して得られた変形場に対して弾性制限を行う弾性制限部、及び／又は、計算して得られた変形場に対して粘性制限を行う粘性制限部をさらに備えたことを特徴とする請求項７に記載の医用画像処理装置。
【請求項９】
前記重み付けが行われた変形場に対して弾性制限を行う弾性制限部、及び／又は、前記重み付が行われた変形場に対して粘性制限を行う粘性制限部をさらに備えたことを特徴とする請求項７に記載の医用画像処理装置。
【請求項１０】
前記弾性制限は、ガウスフィルタであることを特徴とする請求項８又は９に記載の医用画像処理装置。
【請求項１１】
前記粘性制限は、メディアンフィルタであることを特徴とする請求項８又は９に記載の医用画像処理装置。
【請求項１２】
複数の変形後の前記第１画像と少なくとも２つの前記第２画像とを統合して、前記対象物の最終画像を取得する画像統合部をさらに備えたことを特徴とする請求項１〜１１のいずれか一つに記載の医用画像処理装置。
【請求項１３】
前記対象物は心臓と同一位置であることを特徴とする請求項１〜１２のいずれか一つに記載の医用画像処理装置。
【請求項１４】
前記第１画像及び前記第２画像は、医用画像診断装置により得られたデータに基づいて生成された医用画像であることを特徴とする請求項１〜１３のいずれか一つに記載の医用画像処理装置。
【請求項１５】
請求項１〜１４のいずれか一つに記載の医用画像処理装置を含むことを特徴とする医用画像診断装置。
【請求項１６】
対象物の第１画像の第２画像に対する変形場を計算するステップと、
前記対象物の各部位の動き量に応じて前記変形場における重み付けを行うステップと、
前記重み付けが行われた変形場を用いて前記第１画像を変形するステップと
を含んだことを特徴とする画像処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は情報処理分野に関し、特に、医用画像処理装置、医用画像診断装置及び画像処理方法に関する。
【背景技術】
【０００２】
同一の対象物について連続して複数の画像が得られている場合に、全体的な画像の情報を得ることが求められる場合がある。例えば、全体的な画像の情報は、複数の画像を統合することによって得られる。
【０００３】
しかし、そのような画像を得るには一定の時間が必要であるため、対象物が各種の運動をしてしまう可能性があり、全体的な画像を得るために複数の画像を統合することが難しくなる。
【０００４】
したがって、上述した問題を解決するための技術が必要である。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】米国特許出願公開第２０１１／０１０３６６４号明細書
【特許文献２】米国特許出願公開第２０１１／０１４８９２８号明細書
【発明の概要】
【発明が解決しようとする課題】
【０００６】
本発明が解決しようとする課題は、画像のノイズを低減させることができる医用画像処理装置、医用画像診断装置及び画像処理方法を提供することである。
【課題を解決するための手段】
【０００７】
実施形態に係る医用画像処理装置は、変形場計算部と、変形場重み付け部と、画像変形部とを備える。変形場計算部は、対象物の第１画像の第２画像に対する変形場を計算する。変形場重み付け部は、前記対象物の各部位の動き量に応じて前記変形場における重み付けを行う。画像変形部は、前記重み付けが行われた変形場を用いて前記第１画像を変形する。
【図面の簡単な説明】
【０００８】
【図１Ａ】図１Ａは、本発明の一実施形態に基づく医用画像処理装置を示すブロック図である。
【図１Ｂ】図１Ｂは、変形場における重み付けを行う際に適用する重みの特定の動向を示す曲線図である。
【図１Ｃ】図１Ｃは、変形場における重み付けを行う際に適用する重みの他の特定の動向を示す曲線図である。
【図２】図２は、本発明の他の一実施形態に基づく医用画像処理装置のブロック図である。
【図３】図３は、本発明の他の一実施形態に基づく医用画像処理装置のブロック図である。
【図４】図４は、本発明の他の一実施形態に基づく医用画像処理装置のブロック図である。
【図５】図５は、本発明の他の一実施形態に基づく医用画像処理装置のブロック図である。
【図６】図６は、本発明の他の一実施形態に基づく医用画像処理装置のブロック図である。
【図７】図７は、本発明の他の一実施形態に基づく医用画像処理装置のブロック図である。
【図８】図８は、本発明の一実施形態に基づく医用画像診断装置のブロック図である。
【図９】図９は、本発明の一実施形態に基づく画像処理方法のフロー図である。
【図１０】図１０は、本発明の一実施形態に基づいて変形場を計算する流れを示すフロー図である。
【図１１】図１１は、本発明の実施形態における画像処理及び方法を実施するためのコンピュータの構成を示す構成図である。
【発明を実施するための形態】
【０００９】
まず、本発明の基本的な理解のために、本発明の幾つかの実施形態についての概要を紹介する。この概要は、本発明のキーになる領域又は重要な領域を確定するものではなく、また本発明の範囲を限定するものでもない。その目的は、ただ簡略的に説明することによって、その後に説明する更に詳細な記述のイントロダクションのために説明するものである。
【００１０】
本発明の目的は、画像の運動部位をフリーズさせて画像ノイズを減少させることができる医用画像処理装置、医用画像診断装置及び画像処理方法を提供することである。
【００１１】
本出願の一実施形態に係る医用画像処理装置は、変形場計算部と、変形場重み付け部と、画像変形部とを備える。変形場計算部は、対象物の第１画像の第２画像に対する変形場を計算する。変形場重み付け部は、上記対象物の各部位の動き量に応じて上記変形場における重み付けを行う。画像変形部は、重み付けが行われた変形場を用いて上記第１画像を変形する。
【００１２】
本出願の一実施形態に係る画像処理方法は、対象物の第１画像の上記基準画像に対する変形場を計算するステップと、上記対象物の各部位の動き量に応じて上記変形場における重み付けを行うステップと、重み付けが行われた変形場を用いて上記第１画像を変形するステップとを含む。
【００１３】
また、本出願の他の一実施形態によれば、上記画像処理方法を実現するためのコンピュータプログラムを提供する。
【００１４】
さらに、本出願の他の一実施形態によれば、上記画像処理方法を実現するためのコンピュータプログラムのコードを記憶した少なくともコンピュータ読み出し可能な記録媒体の形でプログラム製品を提供する。
【００１５】
以下、図面を参照して、本発明の実施形態について詳細に説明する。
【００１６】
以下、図面を参照しながら本実施形態の説明することにより、さらに本実施形態の目的、特徴、メリットを理解し易くすることができる。図面中の構成は、ただ本実施形態の原理を示すためのものである。図面において、同じ又は類似の技術的特徴又は構成は、同様の又は類似の図面表記を用いて表現することとする。
【００１７】
以下、図面を参照しながら本実施形態について説明をする。なお、説明において、一つの図面又は一つの実施形態において記載した構成や特徴は、一つ又は複数の他の図面又は実施形態において示した構成や特徴の組み合わせることができる。さらに、明瞭にするため、図面や説明において本実施形態と無関係な内容や、当業者にとって周知の構成や処理については、表示や記載を省略する。
【００１８】
図１Ａは、本発明の一実施形態に基づく医用画像処理装置１００を示すブロック図である。
【００１９】
図１Ａに示すように、医用画像処理装置１００は、変形場計算部１１０、変形場重み付け部１２０及び画像変形部１３０を含む。
【００２０】
変形場計算部１１０は、対象物の第１画像の第２画像に対する変形場を計算する。例えば、変形場計算部１１０は、対象物の第２画像を基準画像とし、非剛性レジストレーションに基づいて該対象物の第１画像の基準画像に対する変形場を計算する。
【００２１】
例えば、対象物は、心臓と同一位置で良い。言い換えれば、第１画像と第２画像とは、ある物体の同一位置に対して得られる画像で良い。例えば、ある物体に対して同一方向の連続したマルチスライス画像でも良い。この他、ある物体の異なる方向に連続したマルチスライス画像でも良く、この場合、本文中の医用画像処理装置で処理する前に、同一方向へ揃えるようないくつかの画像変換が必要である。ただし、物体は心臓に限定されず、人体のその他の部位でも適用可能である。なお、ここで言う「同一位置」は、二次元又は三次元で表される同一の領域で良く、二次元又は三次元で表される心臓全体であって良い。
【００２２】
例えば、第１画像と第２画像とは、医用画像診断装置により得られたデータに基づいて生成された医用画像であって良い。また、各画像は二次元画像に限定されず、三次元画像であっても良い。
【００２３】
画像レジストレーションは、空間変換関係に基づいて、剛体レジストレーション（rigid  registration）及び非剛体レジストレーション（non-rigid  registration）の２種類に大別できる。剛体レジストレーションは、画像に対する全体的な変形であり、例えば画像全体にわたる回転、平行移動、拡大等であって、剛体レジストレーションには画像における局所的な変形によるアライメント処理は含まれない。一方、非剛体レジストレーションは、画像に対する各種局所的な変形であって、非剛体レジストレーションは、画像の局所的な変形（たとえば、画像の局所的なねじれ）を利用してアライメント処理を行うものである。
【００２４】
この他、例えば、Ｄｅｍｏｎｓ計算法を採用して非剛性レジストレーションを行うことができる。ただし、Ｄｅｍｏｎｓ計算法に限定されず、適宜その他の非剛体レジストレーション法を適用することができる。例えば、基本関数に基づいて変形領域を描出させてレジストレーションを行う方法でも良く、薄板スプラインやＢ−スプラインのレジストレーションや、弾性変形に基づく物理的模型のレジストレーション等の非剛体レジストレーション法を適用することができる。ここでは逐一列挙はしない。
【００２５】
この他、第２画像を基準画像として使用することには限定されず、第１画像を基準画像とした非剛性レジストレーションを適用して、該対象物の第２画像の基準画像としての第１画像に対する変形場を計算しても良い。
【００２６】
心臓で言えば、取得された心臓の画像シリーズにおいて、拡張末期（end  diastolic）の心臓の画像を選んで選択する。なぜなら、拡張末期の心臓の画像はその他のタイミングの心臓の画像よりも明瞭だからである。人体の心臓以外の部位で言えば、必要に応じて部位を指定して、取得されたあるタイミングの画像を基準画像とする、又は、任意の指定されたあるタイミングの画像を基準画像として、その他の画像をその画像とレジストレーションするようにすれば良い。
【００２７】
変形場重み付け部１２０は、対象物の各部位の動き量に応じて変形場における重み付けを行う。
【００２８】
言い換えれば、対象物の運動部位をフリーズさせて画像のノイズを減少させるため、区別して変形場における運動部位に対応する部分と変形場における非運動部位に対応する部分を処理する。例えば、変形場重み付け部１２０は、対象物の運動部位の変形場の重みを１とし、非運動部位の変形場の重みを１より小さくする（例えば、非運動部位の変形場の重みを０と設定する）。
【００２９】
例えば、動き量が第１閾値よりも大きい部位を運動部位とし、動き量が第２閾値よりも小さい部位を非運動部位とする。ここで、第１閾値は第２閾値値よりも大きいとする。このような状況において、動き量が第１閾値と第２閾値との間の値である部位を移行部位とすることができる。例えば、変形場重み付け部１２０は、対象物の運動部位と非運動部位との間の移行部位の変形場の重みを、運動部位の変形場の重みと非運動部位の変形場の重みとの間の値とすることができる。
【００３０】
当然ながら、運動部位と非運動部位との境界はこれに限らない、例えば、第１閾値が第２閾値と等しい状況下で、動き量が該特定の閾値よりも大きい部位を運動部位とし、動き量が該特定の閾値よりも小さい部位を非運動部位としても良い。変形場重み付け部１２０は、同じように、対象物の運動部位の変形場の重みを１とし、非運動部位の変形場の重みを１よりも小さく（例えば、非運動部位の変形場の重みを０とする）とするようにできる。
【００３１】
ここで、閾値に基づいて運動部位と非運動部位とを区別し、各部位の動き量の値と比例関係を成すように変形場の対応する部分の重みを設定する以外の方法でも良いが、重みの最大値は１である。
【００３２】
また、例えば、変形場重み付け部１２０は、第２画像と複数の第１画像における各部位の動き量が反映される分散画像に基づいて、運動部位と非運動部位とを定め、変形場における重み付けを行ってもよい。
【００３３】
なお、ここで用いられる分散画像は、例えば、第２画像及び複数の第１画像のグレースケールの分散画像でも良い。しかし、ここで用いられる分散画像はこれに限られず、対象物の各部位の動き量を反映できるものであれば良い。
【００３４】
以下の数式（１）により、分散画像の各点における分散値が計算できる：
【００３５】
【数１】
2014174174.tif 000003
【００３６】
ここで、（ｉ，ｊ，ｋ）は位置を表し，Ｄ（ｉ，ｊ，ｋ）は分散値を表し，Ｎは画像の総数を表す。
【００３７】
Ｒn（ｉ，ｊ，ｋ）は、第ｎ番目の画像における（ｉ，ｊ，ｋ）でのグレースケール値を表し、Ｒ（￣付き）（ｉ，ｊ，ｋ）は、（ｉ，ｊ，ｋ）でのグレースケールの平均値を表す。以下の数式（２）により、平均値Ｒ（￣付き）（ｉ，ｊ，ｋ）が計算できる。
【００３８】
【数２】
2014174174.tif 000004
【００３９】
分散値を算出した後に、以下の数式（３）及び（４）により、重み付けが行われた変形場が計算できる。
【００４０】
【数３】
2014174174.tif 000005
【００４１】
【数４】
2014174174.tif 000006
【００４２】
ここで、（ｉ，ｊ，ｋ）は、位置を表し、ＷＦ（ｉ，ｊ，ｋ）は、運動重み付けの前の変形場を表し、ＷＦnew（ｉ，ｊ，ｋ）は、運動重み付けの後の変形場を表す。また、ＭＷ（ｉ，ｊ，ｋ）は、分散値Ｄ（ｉ，ｊ，ｋ）及び動き量を示す運動重みを表す。ｆは、運動重みＭＷが特定の動向を示すようにすることができる関数（例えば、分散画像の正規化関数又は運動／非運動領域の非線形関数）である。
【００４３】
図１Ｂは、変形場における重み付けを行う際に適用する重みの特定の動向を示す曲線図である。例えば、図１Ｂに示すように、運動部位の重みを１とし、非運動部位の重みを０とし、移行部位の重みを０と１との間の線形値とする。
【００４４】
図１Ｃは、変形場における重み付けを行う際に適用する重みの他の特定の動向を示す曲線図である。図１Ｃに示すように、運動部位の重みを１とし、非運動部位の重みを０とし、移行部位の重みを０と１との間の非線形値とする。図１Ｃに示すような重みを適用した場合、得られる第３画像では、運動部位は移行部位との間で接続され、かつ、移行部位は非運動部位との間で接続されるので、さらに滑らかになる。
【００４５】
しかし、ｆは、図１Ｂに示す特定の動向に限らず、動き量が大きくなるほど重みが大きくなるようになっていれば、任意の適切な動向で良いが、重みの最大値は１である。例えば、運動部位の重みを１とし、非運動部位の重みを０．５とし、移行部位の重みを０．５と１との間（例えば、移行部位の重みを０．５と１との間の線形値又は非線形値）にするようにしても良い。
【００４６】
画像変形部１３０は、重み付けが行われた変形場を用いて第１画像を変形する。例えば、画像変形部１３０は、変形場重み付け部１２０によって重み付けが行われた変形場を用いて第１画像を変形して、第３画像を得る。
【００４７】
変形場の計算において、第１画像に基づいて得られる変形場は、第１画像におけるノイズの影響又は変形場の計算方法そのものの欠陥によって、対象物の実際の運動部位で変形が生じるだけでなく、対象物の実際の非運動部位でも変形が生じる可能性がある。例えば、変形場の全ての部分を区別せずに第１画像に適用すれば、得られる第３画像では、実際には運動していない部位にノイズ（変形すべきでない部分にまで変形をしてしまう）が混入する、又は、もともとあるノイズを強調させてしまう。本出願で開示する重み付けにより、すなわち、変形場において運動部位に対応する部分の重みを変形場における非運動部位に対応する部分の重みよりも大きくし、非運動部位と移行部位に対して受け得る変形場の作用を消し去ることによって、変形すべきでない部分まで変形が及ぶ、又は、ノイズが強調されてしまうといった事態をある程度回避することが可能である。言い換えれば、このようにして得られる各第３画像において、運動部位はフリーズされ（変形により、第１画像における運動部位を第２画像において相対する位置又は相対位置の付近まで変形する）、非運動部位と移行部位に対して、新たなノイズの発生又はもともと存在するノイズの強調を抑制することができる。
【００４８】
その上さらに、以下に記載するように、複数の第３画像を統合する際に、非運動部位に対して小さい重みを設定することによってノイズに対する変形効果を抑制することができ、変形により同じ又は近い位置まで変形されることによるノイズを回避し（複数の第１画像における類似のノイズは変形場計算法によりノイズの運動と誤って認識される）、さらに、複数の第１画像におけるノイズの相乗強調も回避できる。また、第１画像においてもともと存在したノイズのランダム性は維持されるので、複数の第３画像を統合するとノイズ低減の効果が期待できる。
【００４９】
以下、心臓のＰＥＴ（Positron  Emission  Tomography）及び／又はＳＰＥＣＴ（Single  Photon  Emission  Computed  Tomography）によるイメージングの場合を例にして説明する。
【００５０】
心臓に対してＰＥＴ及び／又はＳＰＥＣＴイメージングする場合、通常は一定の時間が必要であり、画像を得るまでに、例えば数分間、長い場合は十数分間かかってしまう。
【００５１】
この時間に呼吸や心臓の動きによって患者が動いてしまい、動いた部分の画像が不鮮明になってしまう。
【００５２】
この問題を解決するために、心臓運動の一周を複数の時間間隔で分割し、各時間間隔に対応する画像データを取得するという技術がある。しかし、各時間間隔において取得した画像におけるＳＮ比は小さくなり、フォトンカウント数も低い。
【００５３】
この状況において、変形場計算部１１０により、変形場重み付け部１２０及び画像変形部１３０は、取得された複数の画像に対して上記の処理を行う。具体的には、例えば、運動部位（例えば、心壁等）の変形場の重みを１とし、非運動部位（その他の部位）の変形場の重みを１より小さい値（例えば、０）として運動部位をフリーズさせ、非運動部位と移行部位とは部分的にフリーズ又はフリーズしないようにする。これにより、重み付けが行われた変形場を用いて変形後に得られた画像において、運動部位（例えば、心壁等）の空間上のそれぞれのつなぎ目をつないで、非運動部位と移行部位は空間上で交差する。これらの画像を統合すれば、運動部位（例えば、心壁等）は鮮明で、非運動部位と移行部位とは平均値によりノイズが除去された画像が得られる。
【００５４】
なお、この処理における第１画像及び第２画像はフォトンカウント数が低く、ＳＮ比が低い画像であったが、他の画像であっても、運動をフリーズしてノイズ低減の作用を引き起こすことができる。
【００５５】
なお、この処理における画像は心臓の画像に限らず、その他の動く臓器に対して応用してもよく、例えば、肝臓、肺、胆嚢等の器官に関する画像でも良い。
【００５６】
また、この処理における画像は医用画像に限らず、類似の運動をする物体に対する非医用画像であっても適用できる。
【００５７】
実施形態において、第２画像と複数の第１画像とを統合することができる。画像統合の実施形態については、後述にて図７を参照しながら説明する。
【００５８】
図２は、本発明の他の一実施形態に基づく医用画像処理装置２００のブロック図である。
【００５９】
図２に示すように、医用画像処理装置２００は、変形場計算部２１０、変形場重み付け部２２０及び画像変形部２３０を含む。
【００６０】
図２に示す医用画像処理装置２００と図１Ａに示す医用画像処理装置１００との違いは、図２に示す変形場重み付け部２２０が、アトラス識別部２２１と変形場重み付け部２２２とを含む点である。
【００６１】
アトラス識別部２２１は、アトラスに基づいて第１画像における対象物の各部位を識別する。ここで、アトラスとは、人間の臓器の標準的な形状を表すデータである。
【００６２】
対象物の各部位を識別する場合、例えば、各部位に対応する動き量情報を含むルックアップテーブルを準備し、そのルックアップテーブルに基づいてこれらの部位の運動情報を探し出すことができる。
【００６３】
又は、アトラス中に各部位に対応する動き量情報を含ませておき、識別した対象物の各部位に対応する動き量情報を取得することができる。
【００６４】
変形場重み付け部２２２は、アトラス識別部２２１によって識別された各部位の動き量の空間分布に基づいて、変形場における重み付けを行う。
【００６５】
図３は、本発明の他の一実施形態に基づく医用画像処理装置のブロック図である。
【００６６】
図３に示す通り、医用画像処理装置３００は、変形場計算部３１０、変形場重み付け部３２０、画像変形部３３０及びヒストグラムマッチング部３４０を含む。
【００６７】
図３に示す医用画像処理装置３００と図１Ａに示す医用画像処理装置１００との違いは、医用画像処理装置３００はさらにヒストグラムマッチング部３４０を含む点である。
【００６８】
ある物体（例えば、心壁）に基づいて取得した画像において、グレースケールは大きく変化する。この状況において、変形場を計算する前に第１画像と第２画像との間でヒストグラムマッチングを行い、第１画像と第２画像とのグレースケール範囲を近づける。
【００６９】
ヒストグラムマッチング部３４０は、第１画像のヒストグラムを第２画像のヒストグラムとの間でヒストグラムマッチングを行い、第２画像とグレースケール範囲が近い第１画像を取得する。その状況において、変形場計算部３１０は、第２画像とグレースケール範囲が近い第１画像の第２画像に対する変形場を計算できる。なお、グレースケール範囲が近いということは、言い換えると、グレースケール範囲の差が小さいことである。
【００７０】
図４は、本発明の他の一実施形態に基づく医用画像処理装置４００のブロック図である。
【００７１】
図４に示す通り、医用画像処理装置４００は変形場計算部４１０、変形場重み付け部４２０及び画像変形部４３０を含み、変形場計算部４１０は力場計算部４１１、変形場計算部４１２、画像変形部４１３及び制御部４１４を含む。
【００７２】
力場計算部４１１は、第１画像の基準画像としての第２画像に対する力場を計算する。
【００７３】
Ｄｅｍｏｎｓ計算法を適用した場合、以下の数式（５）により力場ｖ（￣付き）を計算する。
【００７４】
【数５】
2014174174.tif 000007
【００７５】
ここで、定数Ｋは、勾配下降速度を制御するための量であり、ｒは、基準画像としての第一画像のグレースケールである。また、ｆは、第１画像のグレースケールであり、ｓは、前回の反復計算時の計算によって出された変形場である。
【００７６】
変形場計算部４１２は、力場計算部４１１によって算出された力場に基づいて変形場を計算する。
【００７７】
第ｎ回の変形場をＷＦn、第ｎ−１回の変形場をＷＦn-1とすると、それらは以下の数式（６）を満足する。
【００７８】
【数６】
2014174174.tif 000008
【００７９】
ここで、Δnは、第ｎ回の反復時の変形場の増加量を表す。例えば、増加量Δnは、群最適化（group  optimization）を用いて、力場に基づいて算出することができる。
【００８０】
画像変形部４１３は、変形場計算部４１２によって算出された変形場を用いて第１画像を変形する。
【００８１】
制御部４１４は、力場計算部４１１、変形場計算部４１２及び画像変形部４１３を制御し、変形前の第１画像に替えて変形後の第１画像を用いて、所定の条件を満足するまで上記処理を繰り返し行って、最終的な変形場を取得する。
【００８２】
実際に、非剛体レジストレーションの過程において、反復回数の増加にしたがって、第１画像と第２画像の対応する部分の空間位置はどんどん近くなり、両者の画像間の相似度（例えば、正規化に関する係数）はどんどん増加するが、増加の動向は緩やかになっていく。しかし、時には、相似度は反復回数の増加に伴ってわずかに減少することもある。そこで、以下の３つの所定の条件の中の少なくとも一つを採用する。
【００８３】
第１の条件は、反復回数である。例えば、反復回数は４〜２０の間で選択できる。ただし、反復回数はこれに限らず、その他の適切な反復回数を選択することはできる。
【００８４】
第２の条件は、相似度の減少量である。例えば、正規化の関係係数の減少量である。正規化の関係係数が素早く減少しているか否かを判断する。例えば、その減少量が閾値より大きいかどうかを判断する。もし、減少量が該閾値より大きい場合は、その条件を満足していると判断する。
【００８５】
第３の条件は、最近の何回かの反復（例えば、４回）における相似度の増加量である。例えば、正規化の関係係数の増加量である。もし、最近の何回かの反復の正規化の関係係数の平均増加量が所定の閾値（例えば、０．００１）よりも小さい場合、その条件は満足していると判断する。
【００８６】
この他、制御部４１４は、力場計算部４１１を制御して力場に対して粘性制限を行わないようにする。
【００８７】
制御部４１４は、変形場計算部４１２を制御して変形場に対する弾性制限を行わないようにする。すなわち、変形場計算部４１２が最終的な変形場を取得する前に、変形場に対して弾性制限を行わない。
【００８８】
言い換えれば、変形場計算部４１０が行うレジストレーションの反復処理において、弾性制限及び／又は粘性制限を除去する。
【００８９】
取得した画像のフォトンカウント数が低い、及び／又は、画像が小さい状況において、弾性制限及び／又は粘性制限を除去することにより、変形が弱くなるのを回避することができる。
【００９０】
例えば、心臓のＰＥＴ及び／又はＳＰＥＣＴ画像について言えば、二次元的に観察し、心臓の収縮末期の画像における心壁（ｍｙｏｃａｒｄｉｕｍ）及び血液プール（ｂｌｏｏｄ  ｐｏｏｌ）等の部分はいくつかの画素の厚みを有する。同時に、心臓の収縮拡張運動により隣接画素の運動も異なる方向を向いている可能性がある。このような状況において、弾性制限及び／又は粘性制限を除去することにより、血液プール等の部分の変形が弱くなるのを回避することができる。
【００９１】
図５は、本発明の他の一実施形態に基づく医用画像処理装置５００のブロック図である。
【００９２】
図５に示すように、医用画像処理装置５００は、変形場計算部５１０、弾性制限部５５０、粘性制限部５６０、変形場重み付け部５２０及び画像変形部５３０を含む。
【００９３】
図５に示す医用画像処理装置５００と図１Ａに示す医用画像処理装置１００との違いは、医用画像処理装置５００は弾性制限部５５０と粘性制限部５６０を含む点である。
【００９４】
弾性制限部５５０は、変形場計算部５１０が算出した変形場に対して弾性制限を行う。例えば、弾性制限はガウスフィルタで良い。ただし、これに限らず、任意のその他の適切な弾性制限であれば良い。
【００９５】
粘性制限部５６０は、変形場計算部５１０が算出した変形場に対して粘性制限を行う。
【００９６】
例えば、粘性制限はメディアンフィルタで良い。ただし、これに限定されず、任意の適切な粘性制限であれば良い。
【００９７】
変形場計算部５１０におけるレジストレーションの反復処理において力場に対して粘性制限を行わず、かつ、変形場に対して弾性制限を行わずに、レジストレーションの反復処理の後に変形場に対して弾性制限及び粘性制限を行うことにより、変形場における異常点を除去又は減少させることができる。
【００９８】
実施形態において、弾性制限部５５０と粘性制限部５６０は異なるタイミングで機能させることができる。必要に応じて、弾性制限部５５０又は粘性制限部５６０だけを機能させることができる。
【００９９】
この他、実施形態において、必要に応じて、先に粘性制限部５６０により粘性制限を行い、そして弾性制限部５５０により弾性制限を行うことができる。
【０１００】
図６は、本発明の他の一実施形態に基づく医用画像処理装置６００のブロック図である。
【０１０１】
医用画像処理装置６００は、変形場計算部６１０、変形場重み付け部６２０、弾性制限部６５０、粘性制限部６６０及び画像変形部６３０を含む。
【０１０２】
図６に示す医用画像処理装置６００と図５に示す医用画像処理装置５００との違いは、弾性制限部６５０は変形場重み付け部６２０が重み付けを行った変形場に対して弾性制限を行い、粘性制限部６６０は変形場重み付け部６２０が重み付けを行った変形場に対して粘性制限を行い、弾性制限部５５０は変形場計算部５１０が計算して出した変形場に対して弾性制限を行い、粘性制限部５６０は変形場計算部５１０が計算して出した変形場に対して粘性制限を行う点である。
【０１０３】
変形場計算部６１０が行うレジストレーションの反復処理において力場に対して粘性制限を行わず、かつ、変形場に対して弾性制限を行わず、レジストレーションの反復処理の後に重み付けが行われた変形場に対して弾性制限及び粘性制限を行うことにより、重み付けが行われた変形場における異常点を除去又は減少させることができる。
【０１０４】
実施形態において、弾性制限部６５０と粘性制限部６６０は異なるタイミングで機能させることができる。必要に応じて、弾性制限部６５０又は粘性制限部６６０だけを機能させることができる。
【０１０５】
この他、実施形態において、必要に応じて、先に粘性制限部６６０により粘性制限を行い、そして弾性制限部６５０により弾性制限を行うことができる。
【０１０６】
例えば、弾性制限は、ガウスフィルタであり、粘性制限は、メディアンフィルタであっても良い。
【０１０７】
図７は、本発明の他の一実施形態に基づく医用画像処理装置７００のブロック図である。
【０１０８】
医用画像処理装置７００は、変形場計算部７１０、変形場重み付け部７２０、画像変形部７３０及び画像統合部７７０を含む。
【０１０９】
図７に示す医用画像処理装置７００と図１Ａに示す医用画像処理装置１００の差は、医用画像処理装置７００が画像統合部７７０を含む点である。
【０１１０】
医用画像処理装置１００と同様に、変形場計算部７１０は、対象物の第２画像を基準画像とし、非剛性レジストレーションにより該対象物の第１画像の上記基準画像に対する変形場を計算し、変形場重み付け部７２０は、上記対象物の各部位の動き量に応じて上記変形場における重み付けを行い、画像変形部７３０は、重み付けが行われた変形場を用いて上記第１画像を変形して、第３画像を取得する。
【０１１１】
画像統合部７７０は、画像変形部７３０によって得られた複数の第３画像と、基準画像としての少なくとも２つの第２画像とを統合して、対象物の最終画像を取得する。
【０１１２】
言い換えれば、少なくとも２つの第２画像と変形により得られた第３画像とを統合して対象物の最終画像を取得する。
【０１１３】
重み付けが行われた変形場を用いて第１画像を変形して第３画像を取得する中で、運動部分はフリーズされ、かつ、ノイズは低減される。運動部分がフリーズされ、かつ、ノイズが低減された第３画像と第２画像とを統合し、さらに高画質の最終画像を得ることができる。例えば、医学分野において、統合して得られた最終画像において、さらに、特定部位の観察が容易となり、たとえば病変部の観察が容易になる。
【０１１４】
ここでの統合は、複数の画像に対して平均値フィルタリングを実行するものであるが、それに限定されることはなく、適宜任意の適切な方法によって統合することができる。
【０１１５】
図８は、本発明の一実施形態に基づく医用画像診断装置８００のブロック図である。
【０１１６】
本発明の主旨と範囲を明瞭させるために、図８において医用画像診断装置８００のそのほかの可能な部材を省略する。医用画像診断装置８００は医用画像処理装置８１０を含み、画像に対して処理を行う。
【０１１７】
医用画像処理装置８１０は医用画像処理装置１００〜７００の中の任意の一つとすることができる。医用画像診断装置８００は、例えば、Ｘ線イメージング診断装置、超音波（ＵＬ）診断イメージング装置、コンピュータ断層スキャン（ＣＴ）装置、磁気共鳴イメージング（ＭＲＩ）診断装置、又は陽電子放出断層スキャン（Ｐｏｓｉｔｒｏｎ  Ｅｍｉｓｓｉｏｎ  Ｔｏｍｏｇｒａｐｈｙ，ＰＥＴ）装置などであるが、特に限定はされない。
【０１１８】
上記医用画像処理装置を医用画像診断装置中に含める場合、用いられる具体的な手段又は方法は当業者にとって周知のものであり、ここでは重複して説明はしない。
【０１１９】
上文中の実施形態における医用画像処理装置の過程において、明らかな処理や方法を公開した。下文中では、上文中で既に議論された詳細な状況下における方法については重複して説明しない。ただし、記載された医用画像処理装置の過程において方法を開示したが、これらの方法は上記の部分に採用され実行されるとは限らない。例えば、部分的に又は全体的にハードウェア及び／、又はソフトウェアを用いて医用画像処理装置の実施方法を実現させ、以下の議論の画像処理方法も全体的にコンピュータを用いて実行できるプログラムにより実現でき、それらの方法も医用画像処理装置のハードウェア及び／、又はソフトウェアを採用することができる。
【０１２０】
図９は、本発明の一実施形態に基づく画像処理方法のフロー図である。
【０１２１】
図９に示す通り、ステップＳ９１０において、対象物の第２画像を基準画像とし、非剛性レジストレーションに基づいて該対象物の第１画像の基準画像に対する変形場を計算する。例えば、ステップＳ９１０は、変形場計算部１１０、２１０、３１０、４１０、５１０、６１０及び７１０のいずれか一つによって実行される。
【０１２２】
ステップＳ９２０において、対象物の各部位の動き量に応じて変形場における重み付けを行う。例えば、ステップＳ９２０は、変形場重み付け部１２０、２２０、３２０、４２０、５２０、６２０及び７２０のいずれか一つによって実行される。
【０１２３】
例えば、第２画像と複数の第１画像の対象物の各部位の動き量が反映される分散画像に基づいて、変形場における重み付けを行う。
【０１２４】
また、アトラスに基づいて、第１画像において対象物の各部位を識別した後、識別された各部位の動き量の空間分布に基づいて変形場における重み付けを行うことができる。
【０１２５】
付加される重みに関して、例えば、対象物の運動部位の変形場の重みを１とし、非運動部位の変形場の重みを１より小さく（例えば、０）とする。この他、対象物の運動部位と非運動部位との間の移行部位の変形場の重みを運動部位の変形場の重みと非運動部位の変形場の重みとの間の値にすれば良い。
【０１２６】
ステップＳ９３０において、ステップＳ９２０において重み付けが行われた変形場を用いて第１画像を変形する。例えば、ステップＳ９３０は、画像変形部１３０、２３０、３３０、４３０、５３０、６３０及び７３０のいずれか一つによって実行される。
【０１２７】
なお、ステップＳ９１０における変形場の計算の前に、第１画像のヒストグラムと第２画像のヒストグラムとの間でヒストグラムマッチングを行っても良く、第２画像とグレースケール範囲が近い第１画像を得ることができる。このような状況において、ステップＳ９１０において、第２画像とグレースケール範囲が近い第１画像の第２画像に対する変形場を計算することができる。
【０１２８】
実施形態において、該対象物の第１画像の基準画像に対する変形場を計算した後、最終的に得られる変形場に対して弾性制限及び／又は粘性制限を行うことができる。例えば、図５を参照して記述した弾性制限部５５０と粘性制限部５６０によりそれぞれ弾性制限と粘性制限が行うことができる。
【０１２９】
実施形態において、変形場における重み付けを行った後、重み付けが行われた変形場に対して弾性制限及び／又は粘性制限を行う。例えば、図６を参照して記述した弾性制限部６５０と粘性制限部６６０によりそれぞれ弾性制限と粘性制限を行うことが可能である。
【０１３０】
例えば、弾性制限はガウスフィルタで、粘性制限はメディアンフィルタとすることができる。
【０１３１】
この他、複数の第３画像を得た後、複数の第３画像と基準画像のうちの少なくとも２つを統合して対象物の最終画像を得ることができる。例えば、図７を参照して記述した画像統合部７７０により該統合の処理が可能である。
【０１３２】
例えば、対象物は心臓の同一位置である。この他、第２画像と上記第１画像は医用画像診断装置により取得されたデータにより生成された医用画像である。
【０１３３】
図１０は、本発明の一実施形態に基づいて変形場を計算する流れＳ９１０を示すフロー図である。
【０１３４】
図１０に示す通り、ステップＳ９１１において、第１画像の基準画像としての第２画像に対する力場を計算する。例えば、ステップＳ９１１は、図４を参照して説明した力場計算部４１１によって実行される。
【０１３５】
ステップＳ９１２において、ステップＳ９１１で算出された力場に基づいて変形場を計算する。例えば、ステップＳ９１２は、図４を参照して説明した変形場計算部４１２によって実行される。
【０１３６】
ステップＳ９１３において、ステップＳ９１２で算出された変形場を用いて第１画像を変形する。例えば、ステップＳ９１３は、図４を参照して説明した画像変形部４１３によって実行される。
【０１３７】
ステップＳ９１４において、所定の条件を満足するか否かを判定する。例えば、所定の条件を満足しない場合は（ステップＳ９１４における「Ｎｏ」）、変形前の第１画像に替えて変形後の第１画像を用いて、所定の条件を満足するまで上記ステップＳ９１１〜Ｓ９１３の処理を繰り返し、最終的な変形場を取得する。一方、所定の条件を満足する場合は（ステップＳ９１４における「Ｙｅｓ」）、最終的な変形場が取得され、フローは終了となる。例えば、ステップＳ９１４は、図４を参照して説明した制御部４１４によって実行される。
【０１３８】
なお、ステップＳ９１１で計算して出された力場に対しては粘性制限を行わない、及び／又は、最終的な変形場を得る前に変形場に対して弾性制限を行わない。
【０１３９】
上記画像処理方法の各ステップのより具体的な詳細及び多くの可能なステップに関しては、以上の本発明の実施形態による医用画像処理装置における各部材に対する関連記述を参考することができ、ここでは重複して説明はしない。
【０１４０】
上記医用画像処理装置が医用画像診断装置に配置される場合、そこで用いられる具体的な手段又は方法は当業者にとって周知のものであり、ここでは重複して説明はしない。
【０１４１】
一例として、上記画像処理方法の各ステップ及び上記医用画像処理装置の各構成及び／又は部分はソフトウェア、ファームウェア、ハードウェア又はそれらの組み合わせとして実施しても良い。ソフトウェア又はファームウェアを介して実現した場合、上記方法のソフトウェアアプログラムを実施するため、メモリ媒体から又はネットワークを介して専用のハードウェア構造のコンピュータ（例えば、図１１に示す汎用コンピュータ１１００）へダウンロードして構成することができ、該コンピュータに各種プログラムがダウンロードされた状態で、各種機能等を実施することができる。
【０１４２】
図１１において、演算処理部（即ち、ＣＰＵ）１１０１は、読み取り専用メモリ（ＲＯＭ）１１０２の中に記憶されているプログラム、又は、メモリ部１１０８から読み書き兼用メモリ（ＲＡＭ）１１０３へ書き込まれたプログラムに基づいて、各種処理を実施する。ＲＡＭ１１０３では、必要に応じて、ＣＰＵ１１０１が各種処理等を実施するときに必要なデータも記憶しておく。ＣＰＵ１１０１、ＲＯＭ１１０２及びＲＡＭ１１０３は、綜合ライン１１０４を経由してそれぞれ接続されている。入力／出力インターフェース１１０５も、綜合ライン１１０４につながっている。
【０１４３】
下記の各部は、入力／出力インターフェース１１０５に接続されている：入力部１１０６（キーボード、マウス等を含む）、出力部１１０７（モニタ、例えば、ブラウン管（ＣＲＴ）、液晶モニタ（ＬＣＤ）等や、スピーカ等を含む）、メモリ部１１０８（キーボードを含む）、通信部１１０９（ネットワークインターフェースカード、例えば、ＬＡＮカード、モデム等）。通信部１１０９は、ネットワーク（例えば、インターネット）を介して通信処理を実施する。必要に応じて、駆動器１１１０も入力／出力インターフェース１１０５に接続可能である。取り外し可能な媒体１１１１は、例えば、磁気ディスク、光ディスク、ＭＯ、半導体メモリ等であって、必要に応じて駆動器１１１０に装着され、必要に応じてコンピュータプログラムを読み出して、メモリ部１１０８へダウンロードされる。
【０１４４】
ソフトウェアを介して上記システム処理を実施する場合、ネットワーク（例えば、インターネット又は記憶媒体（例えば、取外し可能な媒体１１１１））からプログラムをダウンロードしても良い。
【０１４５】
当業者においては、このような記憶媒体は図１１に示すようなプログラムを記憶した記憶媒体は、装置とは離れたところからユーザにプログラムを提供する取り外し可能な媒体１１１１に限らない。取り外し可能な媒体１１１１の例としては、磁気ディスク（フロッピー（登録商標）ディスク、光ディスク（ＣＤ−ＲＯＭやＤＶＤを含む）、磁気光ディスク（ＭｉｎｉＤｉｓｃ（ＭＤ、登録商標）を含む）らを含む。また、記憶媒体はＲＯＭ１１０２であっても良く、メモリ部１１０８に含まれるハードディスク等、その中にプログラムが記憶され、それらを含む装置からユーザへプログラムが送られる形態でも良い。
【０１４６】
本発明では、更に、メモリとして、機器読み取り可能なコマンドコードを記憶しているプログラム製品でも応用でき、上記コマンドコードが機器を介して読み取られると、本発明の実施形態における画像分割法が実施される。
【０１４７】
上記機器読み取り可能なコマンドコードを記憶しているプログラム製品を受け入れるための記憶媒体も本発明に適用できる。その記憶媒体は、ハードディスク、光ディスク、磁気光ディスク、メモリカード、メモリスティックには限定されない。
【０１４８】
上記の具体的実施形態においては、一つの実施方法に示す特徴について、同様の方法を一つ又は複数の他の実施方法の中で適用したり、その他の実施方法と組み合わせたり、又はその他の実施方法における特徴に替えるといったことも可能である。
【０１４９】
さらに、“包含する／含む”といった用語を使用したときは、特徴、構成、ステップ又は構造の存在を指し示す。ただし、その他の特徴、構成、ステップ又は構造の存在や付加の排除を意味するものではない。
【０１５０】
上記実施形態においては、数字構成の図番記号を用いて各ステップや構成を表記している。ただし、これらの図番記号は単なる説明や画図の都合への考慮によるものであって、その順序やいかなるほかの限定を表すものではない、と当業者は理解すべきである。
【０１５１】
このほか、本実施形態の方法は、詳細な説明の欄において説明された時間順序に沿って実施されるものに限らず、その他の時間順序に沿って、同時に、又は独立して実施されても良い。それゆえ、本願の詳細な説明において説明された方法の実施順序は、本実施形態の技術範囲に対する構成を制限するものではない。
【０１５２】
上記では、既に、本実施形態の具体的実施形態の説明をもって、本実施形態の説明を行っているものの、上記の全ての実施形態は全て単なる例示に過ぎず、限定するものではない。当業者は、特許請求の主旨や範囲において、本実施形態の各種手直し、改良又は同等物の設計を行うことが可能である。これらの手直し、改良又は同等物は、本実施形態の保護範囲内に含まれるものである。
【符号の説明】
【０１５３】
１００  医用画像処理装置
１１０  変形場計算部
１２０  変形場重み付け部
１３０  画像変形部
【図１Ａ】
2014174174.tif 000009
【図１Ｂ】
2014174174.tif 000010
【図１Ｃ】
2014174174.tif 000011
【図２】
2014174174.tif 000012
【図３】
2014174174.tif 000013
【図４】
2014174174.tif 000014
【図５】
2014174174.tif 000015
【図６】
2014174174.tif 000016
【図７】
2014174174.tif 000017
【図８】
2014174174.tif 000018
【図９】
2014174174.tif 000019
【図１０】
2014174174.tif 000020
【図１１】
2014174174.tif 000021