JP6054244B2 - 医用画像解析装置 - Google Patents
医用画像解析装置 Download PDFInfo
- Publication number
- JP6054244B2 JP6054244B2 JP2013101689A JP2013101689A JP6054244B2 JP 6054244 B2 JP6054244 B2 JP 6054244B2 JP 2013101689 A JP2013101689 A JP 2013101689A JP 2013101689 A JP2013101689 A JP 2013101689A JP 6054244 B2 JP6054244 B2 JP 6054244B2
- Authority
- JP
- Japan
- Prior art keywords
- blood vessel
- pixel
- time
- region
- arterial
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 210000004204 blood vessel Anatomy 0.000 claims description 312
- 230000007704 transition Effects 0.000 claims description 262
- 238000010191 image analysis Methods 0.000 claims description 50
- 239000002872 contrast media Substances 0.000 claims description 30
- 238000003384 imaging method Methods 0.000 claims description 28
- 210000001367 artery Anatomy 0.000 claims description 26
- 230000002093 peripheral effect Effects 0.000 claims description 22
- 230000017531 blood circulation Effects 0.000 claims description 15
- 238000000605 extraction Methods 0.000 claims description 14
- 239000000284 extract Substances 0.000 claims description 11
- 230000002792 vascular Effects 0.000 claims description 7
- 210000001147 pulmonary artery Anatomy 0.000 description 61
- 238000004458 analytical method Methods 0.000 description 54
- 210000001710 bronchial artery Anatomy 0.000 description 20
- 230000009026 tissue transition Effects 0.000 description 20
- 238000000034 method Methods 0.000 description 19
- 230000010412 perfusion Effects 0.000 description 17
- 210000004072 lung Anatomy 0.000 description 13
- 238000010586 diagram Methods 0.000 description 12
- 238000012545 processing Methods 0.000 description 12
- 210000000709 aorta Anatomy 0.000 description 9
- 239000008280 blood Substances 0.000 description 8
- 210000004369 blood Anatomy 0.000 description 8
- 230000000694 effects Effects 0.000 description 8
- 238000012905 input function Methods 0.000 description 8
- 230000004087 circulation Effects 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 210000004556 brain Anatomy 0.000 description 3
- 238000002591 computed tomography Methods 0.000 description 3
- 210000004185 liver Anatomy 0.000 description 3
- 101150025379 RPA1 gene Proteins 0.000 description 2
- 101100360219 Schizosaccharomyces pombe (strain 972 / ATCC 24843) rpp101 gene Proteins 0.000 description 2
- 101100033861 Schizosaccharomyces pombe (strain 972 / ATCC 24843) ssb1 gene Proteins 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 238000002595 magnetic resonance imaging Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000008520 organization Effects 0.000 description 2
- 210000003240 portal vein Anatomy 0.000 description 2
- 238000005316 response function Methods 0.000 description 2
- 101100235549 Caenorhabditis elegans lin-53 gene Proteins 0.000 description 1
- 101100300924 Caenorhabditis elegans rba-1 gene Proteins 0.000 description 1
- 101150085800 RPA2 gene Proteins 0.000 description 1
- 101100251675 Schizosaccharomyces pombe (strain 972 / ATCC 24843) rpp201 gene Proteins 0.000 description 1
- 101100033868 Xenopus laevis rpa2-a gene Proteins 0.000 description 1
- 230000008321 arterial blood flow Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 210000001627 cerebral artery Anatomy 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 230000010339 dilation Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 210000002216 heart Anatomy 0.000 description 1
- 210000002767 hepatic artery Anatomy 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 210000003734 kidney Anatomy 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Description
[構成]
図1を参照して、第1の実施形態に係る医用画像解析装置1の構成を説明する。
読出部10は、被検体の複数の領域について重畳領域を有するようにそれぞれ異なる時間に被検体に造影剤を投与して撮影された複数の時系列画像を記憶部2から読み出す。また、読出部10は、さらに造影剤が投与されていない被検体の複数の領域全体を表す全体領域画像を記憶部2から読み出す。ここで、時系列画像は、被検体の一つの領域について一定の時間間隔で撮影された複数の静止画像(フレーム)を時系列的に対応付けた動画像である。被検体の複数の領域についてこの時系列画像を撮影したとき、自ずと時系列画像は複数となる。全体領域画像は、例えばX線CT装置によるヘリカルスキャンによって撮影された静止画像である。複数の時系列画像と全体領域画像とは医用画像解析装置1の内部又は外部に存在する記憶部2に記憶されている。記憶部2が医用画像解析装置1の外部に存在するとき、読出部10は一般的な通信手段を介して複数の時系列画像と全体画像とを記憶部2から読み出す。
第1の動脈推移情報部11は、読出部10が読み出した複数の時系列画像のうち一部の時系列画像中に指定された動脈領域に基づいて動脈領域の画素値の推移を表す第1の動脈推移情報を求める。換言すると、第1の動脈推移情報部11は、動脈領域が表す被検体での位置における造影剤の濃度の推移情報を求める。この処理の例として、第1の動脈推移情報部11は、各フレームにおける動脈領域に含まれる画素の画素値の平均値を算出し、その平均値が時系列的に変化した情報を推移情報として求める。動脈領域とは、解析対象の組織に対する血液流入経路が時系列画像において表された領域である。例えば、組織が脳であるとき、脳動脈が表された領域が動脈領域であり、組織が肝臓であるとき、肝動脈が表された領域と門脈が表された領域とが動脈領域である。門脈は動脈ではないがここでは便宜上動脈領域とする。また、組織が肺であるとき、肺には血液流入経路が肺動脈と気管支動脈との二つの動脈がある。通常、気管支動脈が表された領域を動脈領域として指定することは困難であるので、気管支動脈の代替として大動脈が表された領域が動脈領域として指定される。さらに、肺動脈が表された領域が動脈領域として指定される。複数の時系列画像は、組織の全体領域を複数の領域に分けて撮影された画像であるので、通常、複数の時系列画像の全てに動脈領域を指定できるわけではない。それにより、動脈領域は複数の時系列画像の一部の時系列画像すなわち組織の一部の領域を撮影した時系列画像に指定される。従って、複数の時系列画像には、動脈領域が指定された時系列画像(動脈領域を含む領域について撮影された時系列画像)と動脈領域が指定されない時系列画像(動脈領域を含まない領域について撮影された時系列画像)とが存在する。動脈領域は、例えばユーザが或るフレームを目視しながら操作部3を操作することによって指定される。また、例えば第1の動脈推移情報部11が臨床的な統計データを参照して自動的に動脈領域を指定してもよい。図2は、被検体の肺を解析対象の組織として複数の時系列画像を撮影した例を表す模式図である。この例では、肺LUについて第1の領域A1、第2の領域A2及び第3の領域A3の三つの領域に分けて撮影された複数の時系列画像のうち、第1の領域A1について撮影された時系列画像には大動脈が表された領域が大動脈の動脈領域AOとして指定され、第2の領域A2について撮影された時系列画像には肺動脈が表された領域が肺動脈の動脈領域PAとして指定されている。なお、第1の領域A1について撮影された時系列画像を第1の時系列画像、第2の領域A2について撮影された時系列画像を第2の時系列画像、第3の領域A3について撮影された画像を第3の時系列画像とする。
血管画素選択部12は、読出部10が読み出した複数の時系列画像の画素のうち血管を表す画素である血管画素を選択する。例えば血管画素選択部12は、複数の時系列画像の画素のうち血管の分岐点を表す画素を血管画素として選択する。血管画素選択部12が血管画素を選択する方法として、前述の非特許文献1に一例が開示されている。例えば、血管画素選択部12は、組織の血管をグラフ化し、そのグラフの分岐点を得る。血管画素選択部12は、得られた分岐点近傍の一定の領域における距離変換を行い、一定の領域内の枝それぞれが円柱であるとみなして積集合領域を求める。血管画素選択部12は、求めた積集合領域の重心に該当する画素を血管画素として選択する。
血管画素対応部13は、血管画素選択部12により選択された血管画素のうち、互いに同じ重畳領域に係る一方の時系列画像の血管画素と他方の時系列画像の血管画素とを対応付ける。血管対応部が血管画素を対応付ける方法として、前述の特許文献1に一例が開示されている。例えば、血管画素対応部13は、ランドマークとしての血管画素のうち何れか二つの距離や、何れか三つ以上が形成する多角形の角度に基づいて幾つかの血管画素を除外し、残った血管画素について対応付ける。すなわち、互いに同じ重畳領域に係る一方の時系列画像の血管画素と他方の時系列画像の血管画素とについて、前述の距離や角度を比較して対応付ける。
組織推移情報部14は、読出部10が読み出した複数の時系列画像を受け、被検体の組織における画素値の時系列変化を示す組織推移情報を取得する。組織推移情報部14は、時系列画像の複数フレーム間の画素を対応付けし、各画素の画素値が時系列的に変化した推移情報を取得する。つまり、被検体の組織における造影剤の濃度の推移情報を取得する。組織推移情報部14は血管推移情報部140を有する。
血管推移情報部140は、血管画素対応部13により対応付けられた血管画素の画素値の推移を表す血管推移情報を求める。換言すると、血管推移情報部140は、複数フレーム間の血管画素を対応付けし、各血管画素の画素値が時系列的に変化した情報を血管推移情報として求める。それにより、血管推移情報部140は、血管画素が表す血管における造影剤の濃度の推移情報を求める。
第2の動脈推移情報部15は、複数の時系画像のそれぞれが撮影された時間を表す時間情報と、第1の動脈推移情報部11による第1の動脈推移情報と、血管推移情報部140による血管推移情報とに基づいて、複数の時系列画像のうち動脈領域が設定された一部の時系列画像以外の時系列画像が撮影された時間における動脈領域に相当する第2の動脈推移情報を求める。第2の動脈推移情報とは、動脈領域が指定されていない領域について時系列画像が撮影された時間に、他の領域において指定された動脈領域が表す被検体での位置における造影剤の濃度の推移を表す情報である。例えば、図2に示した例の第1の領域A1について撮影した時系列画像では、第2の領域A2に指定された肺動脈の動脈領域PAが表す被検体での位置における造影剤の濃度の該時系列画像を撮影している時間における推移を表す情報が第2の動脈推移情報である。また、第2の領域A2について撮影した時系列画像では、第1の領域A1に指定された大動脈の動脈領域AOが表す被検体での位置における造影剤の濃度の該時系列画像を撮影している時間における推移を表す情報が第2の動脈推移情報である。また、第3の領域A3について撮影した時系列画像では、第1の領域A1に指定された大動脈の動脈領域AOと第2の領域A2に指定された肺動脈の動脈領域PAとのそれぞれが表す被検体での位置における造影剤の濃度の該時系列画像を撮影している時間における推移を表す情報が第2の動脈推移情報である。
解析部16は、第1の動脈推移情報部11による第1の動脈推移情報と、第2の動脈推移情報部15による第2の動脈推移情報と、組織推移情報部14による組織推移情報と、指定された解析手法とに基づいて複数の時系列画像それぞれについてパフュージョン解析を行う。解析手法には、デコンボリューション法やMaximum Slope法などがあるが、解析手法はユーザが操作部3を操作することによって指定されてもよく、解析部16にプリセットされていてもよい。例えば解析部16は、指定された解析手法に基づいて、第1の時系列画像については「C1pa(t)」と「C1ao(t)」とを入力関数とし、第2の時系列画像については、「C2pa(t)」と「C2ao(t)」とを入力関数とし、第3の時系列画像については、「C3pa(t)」と「C3ao(t)」とを入力関数としてパフュージョン解析を時系列画像の画素ごとに行う。解析部16は、各画素が表す被検体の組織の血流量や血液量などの血流動態を解析結果として画像生成部17に出力する。
画像生成部17は、血管画素対応部13によって位置合わせされた複数の時系列画像と解析部16による解析結果を受け、被検体の組織の血流動態を表すマップを生成する。マップとしては、例えば肺の組織の血液量を表す血流量マップや血液量を表す血液量マップなどが生成される。画像生成部17は生成したマップを表示部4に表示させる。
制御部は、医用画像解析装置1の各部を制御する。制御部は、例えば、処理装置と記憶装置とを含んで構成される。処理装置としては、例えば、CPU(Central Processing Unit)、GPU(Graphic Processing Unit)、又はASIC(Application Specific Integrated Circuit)が用いられる。記憶装置は、例えば、ROM(Read Only Memory)、RAM(Random Access Memory)、HDD(Hard Disc Drive)を含んで構成される。記憶装置には、医用画像解析装置1の各部の機能を実行するためのコンピュータプログラムが記憶されている。処理装置は、これらコンピュータプログラムを実行することで、上記機能を実現する。
操作部3は、ユーザによる操作を受けて、この操作の内容に応じた信号や情報を装置各部に入力する。操作部3は、例えば、キーボード、マウス、タッチパネルなどによって構成される。また、操作部3は、必ずしも医用画像解析装置1の一体として備えられる必要はなく、一般的なインターフェイスを介して信号や情報を装置各部に入力する構成でもよい。
表示部4は、例えば、CRT(Cathode Ray Tube)やLCD(Liquid Crystal Display)などの表示デバイスで構成される表示デバイスである。表示部4は、必ずしも医用画像解析装置1の一体として備えられる必要はなく、一般的なインターフェイスを介して画像を表示する構成でもよい。
この実施形態の医用画像解析装置1の動作について説明する。図5は、この実施形態の動作を表すフローチャートである。
読出部10は、被検体の複数の領域について重畳領域を有するようにそれぞれ異なる時間に被検体に造影剤を投与して撮影された複数の時系列画像を記憶部2から読み出す。また、読出部10は、さらに造影剤が投与されていない被検体の複数の領域全体を表す全体領域画像を記憶部2から読み出す。
第1の動脈推移情報部11は、読出部10が読み出した複数の時系列画像のうち一部の時系列画像中に指定された動脈領域に基づいて動脈領域の画素値の推移を表す第1の動脈推移情報を求める。
血管画素選択部12は、読出部10が読み出した複数の時系列画像の画素のうち血管を表す画素である血管画素を選択する。また、血管画素選択部12は、読出部10が読み出した全体領域画像の画素のうち血管を表す画素である全体領域血管画素を選択する。
血管画素対応部13は、血管画素選択部12により選択された血管画素と全体領域血管画素とを対応付けることによって、重畳領域に係る一方の時系列画像の血管画素と他方の時系列画像の血管画素とを対応付け、全体領域画像における重畳領域に該当する該当領域を特定する。
血管画素選択部12は、血管画素対応部13が特定した全体画像における重畳領域に該当する該当領域の画素のうち血管を表す画素である第1の重畳領域血管画素を血管画素対応部13に対応付けられた全体領域血管画素よりも高い密度で新たに選択するとともに、複数の時系列画像の重畳領域に含まれる画素のうち血管を表す画素である第2の重畳領域血管画素を血管画素よりも高い密度で新たに選択する。
血管画素対応部13は、第1の重畳領域血管画素の一部と第2の重畳領域血管画素の一部とを対応付けることによって、重畳領域に係る一方の時系列画像の血管画素と他方の時系列画像の血管画素とを対応付ける。
血管画素対応部13は、第1の重畳領域血管画素のうち第2の重畳領域血管画素に対応付けられていない第1の未対応画素を含む領域の画像である第1の周辺画像と、複数の時系列画像のうち周辺画像に相当する画素である第2の未対応画素を含む領域の画像である第2の周辺画像とを比較し、第1の周辺画像と第2の周辺画像との一致度が指定された一致度を超えたとき、第1の未対応画素と第2の未対応画素とをさらに対応付ける。
血管推移情報部140は、血管画素対応部13により対応付けられた血管画素の画素値の推移を表す血管推移情報を求める。また、組織推移情報部14は、読出部10が読み出した複数の時系列画像を受け、被検体の組織における画素値の時系列変化を示す組織推移情報を求める。
第2の動脈推移情報部15は、血管の種別である血管種別の特性を表す血管種別情報を予め記憶し、血管推移情報部140による血管推移情報と血管種別情報とに基づいて、血管画素対応部13が対応付けた血管画素が表す血管の血管種別を判断する。
第2の動脈推移情報部15は、複数の時系画像のそれぞれが撮影された時間を表す時間情報と、第1の動脈推移情報部11による動脈推移情報と、血管推移情報部140による血管推移情報とに基づいて、複数の時系列画像のうち動脈領域が設定された一部の時系列画像以外の時系列画像が撮影された時間における動脈領域に相当する第2の動脈推移情報を求める。
解析部16は、第1の動脈推移情報部11による第1の動脈推移情報と、第2の動脈推移情報部15による第2の動脈推移情報と、組織推移情報部14による組織推移情報と、指定された解析手法とに基づいて複数の時系列画像それぞれについてパフュージョン解析を行う。
画像生成部17は、血管画素対応部13によって位置合わせされた複数の時系列画像と解析部16による解析結果を受け、被検体の組織の血流動態を表すマップを生成し、表示部4に表示させる。
この実施形態の医用画像解析装置1の効果を説明する。
[構成]
図6を参照して、第2の実施形態に係る医用画像解析装置1の構成を説明する。この実施形態の医用画像解析装置1は、時相対応部19、第1の動脈推移情報部11、組織領域抽出部20、動脈画素情報部21、及び第2の動脈推移情報部15の構成が第1の実施形態の医用画像解析装置1と異なる。その他の構成要素は第1の実施形態の医用画像解析装置1と同様である。以下、第1の実施形態の医用画像解析装置1と異なる構成について特に説明する。また、図7の模式図のように、肺動脈の起始部を表す肺動脈の動脈領域PAが第1の領域A1内、かつ、第2の領域A2外に指定される場合の例を説明する。
時相対応部19は、読出部10が読み出した複数の時系列画像のそれぞれについて、フレーム同士の画素の位置合わせを行う。例えば時相対応部19は、画像相関処理によって一方のフレームと他方のフレームとの位置合わせを行う。時相対応部19は、時系列画像の最初のフレームから最後のフレームまでこの位置合わせを行う。それにより、時系列画像の撮影中に呼吸などによって組織が動いた場合において、その組織をフレーム間で対応付けることができる。時相対応部19は、時系列画像のそれぞれについてこの対応付けを行い、時相対応情報として動脈画素情報部21へ出力する。
第1の動脈推移情報部11は、読出部10が読み出した複数の時系列画像のうち一部の時系列画像中に指定された動脈領域に基づいて動脈領域の画素値の推移を表す第1の動脈推移情報を求める。例えば第1の動脈推移情報部11は、第1の領域A1に指定された肺動脈の動脈領域PAに含まれる画素の平均値を算出し、その平均値の時間変化を示す情報を第1の動脈推移情報として求める。この第1の動脈推移情報を「C1pa(t):t1s≦t≦t1e」とする。ここで時刻t1sは、第1の時系列画像の撮影開始時刻を示し、時刻t1eは、第1の時系列画像の撮影終了時刻を示す。
組織推移情報部14は、読出部10が読み出した複数の時系列画像と時相対応部19からの時相対応情報とを受け、被検体の組織における画素値の時系列変化を示す組織推移情報を取得する。
組織領域抽出部20は、読出部10が読み出した複数の時系列画像のそれぞれについて、解析対象の組織を表す画像領域を抽出する。例えば組織領域抽出部20は、解析対象の組織形状を表す形状データを予め記憶し、時系列画像に描出された画像形状と形状データとを比較し、解析対象を表す画像領域を抽出する。解析対象が肺LUであるとき、組織領域抽出部20は、時系列画像のうち、肺LUを表す画像領域を抽出する。組織領域抽出部20は、抽出した画像領域を動脈画素情報部21へ出力する。
動脈画素情報部21は、読出部10が読み出した複数の時系列画像と、時相対応部19による時相対応情報と、組織領域抽出部20による画像領域とを受け、画像領域の各画素の画素値の推移を表す動脈画素推移情報を求める。このとき、動脈画素情報部21は、複数の時系列画像のそれぞれに、グレースケールモーフォロジー処理、または一般的なその他のフィルタ処理を施して画像領域に表された組織形状をさらに明確に抽出し、各画素の動脈画素推移情報を求めてもよい。なお、動脈画素情報部21は、グレースケールモーフォロジー処理を複数の時系列画像に施す場合、グレースケールダイレーション処理を施してもよい。
第2の動脈推移情報部15は、複数の時系列画像のすべての撮影時間に亘る動脈の推移情報を第2の動脈推移情報として求める。第2の動脈推移情報部15は、例えば次式によって、各々の時系列画像について造影剤投与時刻と撮影開始時刻とのずれ時間を補正する。
解析部16は、第1の動脈推移情報部11による第1の動脈推移情報と、第2の動脈推移情報部15による第2の動脈推移情報と、組織推移情報部14による組織推移情報と、指定された解析手法とに基づいて複数の時系列画像それぞれについてパフュージョン解析を行う。このとき解析部16は、第2の動脈推移情報部15から受けた推移情報Cpa(t)のうち、時刻t1sから時刻t1eまでの時間範囲に表される推移情報を、第1の時系列画像の入力関数としてパフュージョン解析を行う。また、解析部16は、第2の動脈推移情報部15から受けた推移情報Cpa(t)のうち、時刻t2sから時刻t2eまでの時間範囲に表される推移情報を、第2の時系列画像の入力関数としてパフュージョン解析を行う。解析部16はパフュージョン解析を時系列画像の画素毎に行い、各画素が表す被検体の組織の血流量や血液量などの血流動態を解析結果として画像生成部17へ出力する。
この実施形態の医用画像解析装置1の動作について説明する。図9は、この実施形態の動作を表すフローチャートである。
読出部10は、被検体の複数の領域について重畳領域を有するようにそれぞれ異なる時間に被検体に造影剤を投与して撮影された複数の時系列画像を記憶部2から読み出す。
時相対応部19は、読出部10が読み出した複数の時系列画像のそれぞれについて、フレーム同士の画素の位置合わせを行う。
第1の動脈推移情報部11は、読出部10が読み出した複数の時系列画像のうち一部の時系列画像中に指定された動脈領域に基づいて動脈領域の画素値の推移を表す第1の動脈推移情報を求める。また、組織推移情報部14は、読出部10が読み出した複数の時系列画像を受け、被検体の組織における画素値の時系列変化を示す組織推移情報を求める。
組織領域抽出部20は、読出部10が読み出した複数の時系列画像のそれぞれについて、解析対象の組織を表す画像領域を抽出する。組織領域抽出部20は、抽出した画像領域を動脈画素情報部21へ出力する。
動脈画素情報部21は、読出部10が読み出した複数の時系列画像と、時相対応部19による時相対応情報と、組織領域抽出部20による画像領域とを受け、画像領域の各画素の画素値の推移を表す動脈画素推移情報を求める。動脈画素情報部21は、求めた推移情報を第2の動脈推移情報部15へ出力する。
第2の動脈推移情報部15は、複数の時系列画像のすべての撮影時間に亘る動脈領域の推移情報を第2の動脈推移情報として求める。第2の動脈推移情報部15は、求めた推移情報を解析部16へ出力する。
解析部16は、第2の動脈推移情報部15による第2の動脈推移情報と、組織推移情報部14による組織推移情報と、指定された解析手法とに基づいて複数の時系列画像それぞれについてパフュージョン解析を行う。
画像生成部17は、血管画素対応部13によって位置合わせされた複数の時系列画像と解析部16による解析結果を受け、被検体の組織の血流動態を表すマップを生成し、表示部4に表示させる。
第2の実施形態の医用画像解析装置1の効果について説明する。医用画像解析装置1は、時相対応部19と、第1の動脈推移情報部11と、組織領域抽出部20と、動脈画素情報部21と、第2の動脈推移情報部15とを有する。時相対応部19は、複数の時系列画像のそれぞれについて、フレーム同士の画素の位置合わせを行う。第1の動脈推移情報部11は、読出部10が読み出した複数の時系列画像のうち一部の時系列画像中に指定された動脈領域に基づいて動脈領域の画素値の推移を表す第1の動脈推移情報を求める。組織領域抽出部20は、複数の時系列画像のそれぞれから、解析対象の組織を表す画像領域を抽出する。動脈画素情報部21は、複数の時系列画像と、時相対応情報と、画像領域とを受け、画像領域の各画素の画素値の推移を表す動脈画素推移情報を求める。第2の動脈推移情報部15は複数の時系列画像のすべての撮影時間に亘る動脈領域の推移情報を第2の動脈推移情報として求める。それにより、医用画像解析装置1は、複数の領域のうち何れか一つの領域に指定された動脈領域の画素値の推移情報に基づいて、動脈領域が指定されなかった領域の時系列画像についての入力関数を求め、該時系列画像についてのパフュージョン解析を行うことができる。
2 記憶部
3 操作部
4 表示部
10 読出部
11 第1の動脈推移情報部
12 血管画素選択部
13 血管画素対応部
14 組織推移情報部
15 第2の動脈推移情報部
16 解析部
17 画像生成部
18 制御部
19 時相対応部
20 組織領域抽出部
21 動脈画素情報部
140 血管推移情報部
Claims (8)
- 被検体の時系列画像を解析して前記被検体の血流動態を求める医用画像解析装置であって、
前記被検体の複数の領域について重畳領域を有するようにそれぞれ異なる時間に前記被検体に造影剤を投与して撮影された複数の時系列画像を読み出す読出部と、
複数の前記時系列画像のうち一部の時系列画像中に指定された動脈領域に基づいて前記動脈領域の画素値の推移を表す第1の動脈推移情報を求める第1の動脈推移情報部と、
複数の前記時系列画像の画素のうち血管を表す画素である血管画素を選択する血管画素選択部と、
前記血管画素選択部により選択された前記血管画素のうち、互いに同じ前記重畳領域に係る一方の時系列画像の前記血管画素と他方の時系列画像の前記血管画素とを対応付ける血管画素対応部と、
前記血管画素対応部により対応付けられた前記血管画素の画素値の推移を表す血管推移情報を求める血管推移情報部と、
複数の前記時系列画像のそれぞれが撮影された時間を表す時間情報と、前記第1の動脈推移情報部による前記第1の動脈推移情報と、前記血管推移情報部による前記血管推移情報とに基づいて、複数の前記時系列画像のうち前記一部以外の時系列画像が撮影された時間における前記動脈領域に相当する第2の動脈推移情報を求める第2の動脈推移情報部と
を有することを特徴とする医用画像解析装置。 - 前記読出部は、さらに、造影剤が投与されていない前記被検体の複数の領域全体を表す全体領域画像を読み出し、
前記血管画素選択部は、さらに、前記全体領域画像の画素のうち血管を表す画素である全体領域血管画素を選択し、
前記血管画素対応部は、前記血管画素選択部により選択された前記血管画素と前記全体領域血管画素とを対応付けることによって、前記重畳領域に係る一方の時系列画像の前記血管画素と他方の時系列画像の前記血管画素とを対応付ける
ことを特徴とする請求項1に記載の医用画像解析装置。 - 前記血管画素対応部は、前記血管画素選択部により選択された前記全体領域血管画素と前記血管画素とを対応付けることによって、前記重畳領域に係る一方の時系列画像の前記血管画素と他方の時系列画像の前記血管画素とを対応付け、前記全体領域画像における前記重畳領域に該当する該当領域を特定し、
前記血管画素選択部は、前記血管画素対応部が特定した前記該当領域の画素のうち血管を表す画素である第1の重畳領域血管画素を前記全体領域血管画素よりも高い密度で新たに選択するとともに、複数の前記時系列画像の前記重畳領域に含まれる画素のうち血管を表す画素である第2の重畳領域血管画素を前記血管画素よりも高い密度で新たに選択し、
前記血管画素対応部は、前記血管画素選択部により新たに選択された前記第1の重畳領域血管画素の一部と前記第2の重畳領域血管画素の一部とを対応付けることによって、前記重畳領域に係る一方の時系列画像の前記血管画素と他方の時系列画像の前記血管画素とを対応付ける
ことを特徴とする請求項2に記載の医用画像解析装置。 - 前記血管画素対応部は、前記血管画素選択部により新たに選択された前記第1の重畳領域血管画素の一部と前記第2の重畳領域血管画素の一部とを対応付けたとき、前記第1の重畳領域血管画素のうち前記第2の重畳領域血管画素に対応付けられていない第1の未対応画素を含む領域の画像である第1の周辺画像と、複数の前記時系列画像のうち前記周辺画像に相当する画素である第2の未対応画素を含む領域の画像である第2の周辺画像とを比較し、前記第1の周辺画像と前記第2の周辺画像との一致度が指定された一致度を超えたとき、前記第1の未対応画素と前記第2の未対応画素とをさらに対応付けることによって、前記重畳領域に係る一方の時系列画像の前記血管画素と他方の時系列画像の前記血管画素とを対応付ける
ことを特徴とする請求項3に記載の医用画像解析装置。 - 前記第2の動脈推移情報部は、血管の種別である血管種別の特性を表す血管種別情報を予め記憶し、前記血管推移情報部による前記血管推移情報と前記血管種別情報とに基づいて、前記血管画素対応部が対応付けた前記血管画素が表す血管の前記血管種別を判断し、判断した前記血管種別を前記血管推移情報に含めて前記第2の動脈推移情報を求めることを特徴とする請求項1〜4のいずれかに記載の医用画像解析装置。
- 前記血管画素選択部は、複数の前記時系列画像の画素のうち血管の分岐点を表す画素を前記血管画素として選択することを特徴とする請求項1に記載の医用画像解析装置。
- 前記血管画素選択部は、複数の前記時系列画像の画素のうち血管の分岐点を表す画素を前記血管画素として選択するとともに、前記全体領域画像の画素のうち血管の分岐点を表す画素を前記全体領域血管画素として選択することを特徴とする請求項2に記載の医用画像解析装置。
- 被検体の時系列画像を解析して前記被検体の血流動態を求める医用画像解析装置であって、
前記被検体の複数の領域について重畳領域を有するようにそれぞれ異なる時間に前記被検体に造影剤を投与して撮影された複数の時系列画像を読み出す読出部と、
前記複数の時系列画像のうち一部の時系列画像中に指定された動脈領域に基づいて前記動脈領域の画素値の推移を表す第1の動脈推移情報を求める第1の動脈推移情報部と、
前記複数の時系列画像のそれぞれについて、フレーム同士の画素の位置合わせを行い、該位置合わせの結果を時相対応情報として求める時相対応部と、
前記複数の時系列画像のそれぞれから、解析対象の組織を表す画像領域を抽出する組織領域抽出部と、
前記複数の時系列画像と、前記時相対応情報と、前記画像領域とを受け、前記画像領域の動脈画素を抽出し、前記動脈画素それぞれの画素値の推移を表す動脈画素推移情報を求める動脈画素情報部と、
複数の前記時系列画像のそれぞれが撮影された時間を表す時間情報と、前記動脈画素推移情報と、前記第1の動脈推移情報とに基づいて、複数の時系列画像のすべての撮影時間に亘る前記動脈領域の推移情報を第2の動脈推移情報として求める第2の動脈推移情報部と
を有することを特徴とする医用画像解析装置。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013101689A JP6054244B2 (ja) | 2013-05-13 | 2013-05-13 | 医用画像解析装置 |
PCT/JP2014/062745 WO2014185424A1 (ja) | 2013-05-13 | 2014-05-13 | 医用画像解析装置 |
US14/940,260 US9734578B2 (en) | 2013-05-13 | 2015-11-13 | Medical image analyzer |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013101689A JP6054244B2 (ja) | 2013-05-13 | 2013-05-13 | 医用画像解析装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2014221120A JP2014221120A (ja) | 2014-11-27 |
JP6054244B2 true JP6054244B2 (ja) | 2016-12-27 |
Family
ID=52121050
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2013101689A Expired - Fee Related JP6054244B2 (ja) | 2013-05-13 | 2013-05-13 | 医用画像解析装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP6054244B2 (ja) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP7270453B2 (ja) * | 2019-04-26 | 2023-05-10 | キヤノン株式会社 | 画像処理装置、画像処理方法およびプログラム |
CN113616226B (zh) * | 2021-09-14 | 2023-06-23 | 上海联影医疗科技股份有限公司 | 一种血管分析方法、系统、设备及存储介质 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002095654A (ja) * | 2000-09-26 | 2002-04-02 | Toshiba Corp | X線画像診断装置 |
JP2003210456A (ja) * | 2002-01-21 | 2003-07-29 | Toshiba Corp | 時系列画像の処理装置 |
JP2004208714A (ja) * | 2002-12-26 | 2004-07-29 | Ge Medical Systems Global Technology Co Llc | X線ctシステムおよびその制御方法 |
JP5322548B2 (ja) * | 2008-09-17 | 2013-10-23 | 株式会社東芝 | X線ct装置、医用画像処理装置および医用画像処理プログラム |
JP5562610B2 (ja) * | 2009-10-22 | 2014-07-30 | 株式会社東芝 | 画像処理装置 |
JP5718576B2 (ja) * | 2010-02-23 | 2015-05-13 | 株式会社東芝 | X線画像診断装置 |
-
2013
- 2013-05-13 JP JP2013101689A patent/JP6054244B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2014221120A (ja) | 2014-11-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP1842164B1 (en) | Image processing system and method for alignment of images | |
US20100034446A1 (en) | System and Method for Coronary Digital Subtraction Angiography | |
US10052032B2 (en) | Stenosis therapy planning | |
JP6958202B2 (ja) | 動態画像処理装置及びプログラム | |
JP6381895B2 (ja) | 画像処理装置、画像処理方法およびx線診断装置 | |
US9569839B2 (en) | Image processing apparatus, method and medical image device | |
US9198603B2 (en) | Device, method and program for searching for the shortest path in a tubular structure | |
JP5215036B2 (ja) | 医用画像処理装置、及び医用画像処理プログラム | |
CN102131462A (zh) | 用于心脏减影的掩模构造 | |
US9691148B2 (en) | Medical imaging analyzer and medical image processing method | |
CN106469449A (zh) | 一种医学图像中病灶的显示方法及装置 | |
WO2004024001A1 (ja) | 血流動態解析装置とその方法、画像診断装置 | |
WO2014185424A1 (ja) | 医用画像解析装置 | |
US20170301083A1 (en) | Image registration device, method, and program | |
JP6054244B2 (ja) | 医用画像解析装置 | |
JP6054245B2 (ja) | 医用画像解析装置 | |
KR20180008134A (ko) | 분획 혈류 예비력 예측 방법 | |
US11056149B2 (en) | Medical image storage and reproduction apparatus, method, and program | |
JP5750381B2 (ja) | 領域抽出処理システム | |
US11813106B2 (en) | Image processing device, and image processing method utilizing time-series computed tomography (CT) images | |
CN113538419B (zh) | 一种图像处理方法和系统 | |
JP2006075390A (ja) | 医用画像表示装置 | |
JP4567395B2 (ja) | 血管撮影画像処理装置及び血管撮影画像処理方法及び血管撮影画像処理プログラム | |
JP2014094229A (ja) | 医用画像解析装置及び医用画像撮影装置 | |
KR20160031775A (ko) | Tof mra 혈관벽에서의 신호 강도 구배 벡터 분석을 통한 혈류 방향 추출 방법 및 이를 적용한 의료 영상 시스템 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20160215 |
|
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A711 Effective date: 20160527 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20161101 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20161130 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6054244 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
LAPS | Cancellation because of no payment of annual fees |