JP2004057340A - Thrombus section image extraction method, image processor and x-ray ct system - Google Patents

Thrombus section image extraction method, image processor and x-ray ct system Download PDF

Info

Publication number
JP2004057340A
JP2004057340A JP2002217726A JP2002217726A JP2004057340A JP 2004057340 A JP2004057340 A JP 2004057340A JP 2002217726 A JP2002217726 A JP 2002217726A JP 2002217726 A JP2002217726 A JP 2002217726A JP 2004057340 A JP2004057340 A JP 2004057340A
Authority
JP
Japan
Prior art keywords
region
thrombus
image
blood vessel
image information
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.)
Pending
Application number
JP2002217726A
Other languages
Japanese (ja)
Inventor
Yasuhiro Imai
今井 靖浩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
GE Medical Systems Global Technology Co LLC
Original Assignee
GE Medical Systems Global Technology Co LLC
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by GE Medical Systems Global Technology Co LLC filed Critical GE Medical Systems Global Technology Co LLC
Priority to JP2002217726A priority Critical patent/JP2004057340A/en
Publication of JP2004057340A publication Critical patent/JP2004057340A/en
Pending legal-status Critical Current

Links

Images

Abstract

<P>PROBLEM TO BE SOLVED: To provide a thrombus section image extraction method, an image processor and an X-ray CT system, for simply measuring the size of a thrombus of aneurysm with high accuracy using image processing technology. <P>SOLUTION: The thrombus section image extraction method includes the steps of: extracting vicinity region containing a blood vessel in a vertical tomographic plane intersecting perpendicularly to the blood vessel (step S404); smoothing the vicinity region (step S405); extracting an approximation region of the thrombus section from the smoothed vicinity region by region expansion processing (step S408); and extracting a circular arc region from the approximate region of the thrombus section by circular arc processing (step S409). The approximate region of the thrombus section is extracted using the region expansion processing, the circular-arc thrombus section with high accuracy is extracted from the approximate region, and the evaluation on the area of thrombus and the quantitative evaluation on a temporal change of the area of thrombus are realized. <P>COPYRIGHT: (C)2004,JPO

Description

【0001】
【発明の属する技術分野】
この発明は、造影剤を用いて撮像を行うX線CT装置の画像処理装置に関する。
【0002】
【従来の技術】
近年、動脈瘤の撮影目的で、血管造影剤を用いたX線CT撮影が行われる。動脈瘤は、血管と、この血管壁に付着し、ゼリー(jelly)状に血液が固まった血栓部分とからなる。ここで、動脈瘤の経過を予測する上で、動脈瘤部の血管径と、血栓の大きさの情報を得ることが重要である。造影剤を用いた撮影では、動脈瘤の血管部分は、血中の造影剤により鮮明に描出される。
【0003】
【発明が解決しようとする課題】
しかしながら、上記従来技術によれば血管壁に付着した血栓部分は、造影効果がなく、周辺組織からの識別に困難が伴っていた。すなわち、血栓部分には、造影剤が流入しないので造影効果は見られなかった。
【0004】
特に、識別しにくい血栓部分の情報を、画像処理技術により抽出し、計測する試みもなされている。しかし、血栓部分には、濃度斑と呼ばれる明るさの濃淡が存在し、画像処理技術を用いた血栓部分の計測を妨げる要因となっている。
【0005】
これらのことから、動脈瘤の血栓の大きさを、画像処理技術を用いて、簡易に、しかも高い確度で測定できる血栓部画像抽出方法、画像処理装置およびX線CTシステムをいかに実現するかが重要となる。
【0006】
この発明は、上述した従来技術による課題を解決するためになされたものであり、動脈瘤の血栓の大きさを、画像処理技術を用いて、簡易に、しかも高い確度で測定できる、血栓部画像抽出方法、画像処理装置およびX線CTシステムを提供することを目的とする。
【0007】
【課題を解決するための手段】
上述した課題を解決し、目的を達成するために、第1の観点の発明にかかる血栓部画像抽出方法は、X線CT装置により被検体の複数の断層画像情報を取得し、前記断層画像情報から血管中の血流方向と垂直に交わる複数の垂直断層画像情報を生成し、前記垂直断層画像情報の画像から前記血管を含む前記血管の近傍領域を抽出し、前記近傍領域の画像を平滑化し、前記平滑化された画像から前記血管に付着する血栓部の近似領域を抽出し、前記血栓領域の辺縁を円弧状に確定すること、を特徴とする。
【0008】
この第1の観点の発明によれば、X線CT装置により被検体の複数の断層画像情報を取得し、この断層画像情報から血管中の血流方向と垂直に交わる複数の垂直断層画像情報を生成し、この垂直断層画像情報の画像から血管を含む血管の近傍領域を抽出し、この近傍領域の画像を平滑化し、この平滑化された画像から血管に付着する血栓部の近似領域を抽出し、この近似領域の辺縁を円弧状に確定することとしているので、血栓部に存在する濃度斑を平滑化し、血栓部の近似領域を抽出し、この抽出部分の辺縁をさらに円弧状に確定して確度の高い血栓領域を抽出することができ、ひいては、血栓面積の評価および血栓面積の時間的変化等の定量的評価を行うことができる。
【0009】
また、第2の観点の発明にかかる画像処理装置は、X線CT装置から被検体の複数の断層画像情報を取得し、前記断層画像情報から血管中の血流方向と垂直に交わる複数の垂直断層画像情報を生成する画像処理装置であって、前記垂直断層画像情報の画像から前記血管を含む前記血管の近傍領域を抽出する抽出手段と、前記近傍領域の画像を平滑化する平滑化手段と、前記平滑化された画像から、指定範囲の画素値を有する領域を抽出する領域拡張手段と、前記近似領域の辺縁を円弧状に確定する確定手段と、を備えることを特徴とする。
【0010】
この第2の観点による発明によれば、抽出手段により、垂直断層画像情報から血管を含む血管の近傍領域を抽出し、平均化手段により、近傍領域の画像を平滑化し、領域拡張手段により、指定範囲の画素値を有する領域を抽出し、確定手段により、血栓領域の辺縁を円弧状に確定することとしているので、血栓部分に存在する濃度斑を平滑化し、血栓部の近似領域を抽出し、この近似領域の辺縁をさらに円弧状に確定して確度の高い血栓領域を抽出することができ、ひいては、血栓面積の評価および血栓面積の時間的変化等の定量的評価を行うことができる。
【0011】
また、第3の観点の発明にかかる画像処理装置は、前記抽出手段に、閾値を越える画素値、あるいは閾値を越えない画素値の画像領域を抽出する閾値手段を備えることを特徴とする。
【0012】
この第3の観点の発明によれば、抽出手段は、閾値手段により、閾値を越える画素値、あるいは閾値を越えない画素値の画像領域を抽出することとしているので、画像情報から骨部領域および背景領域を抽出して、軟部組織部分のみからなる画像情報を生成することができる。
【0013】
また、第4の観点の発明にかかる画像処理装置は、前記抽出手段に、関心領域を設定する関心領域設定手段を備えることを特徴とする。
この第4の観点の発明によれば、抽出手段は、関心領域設定手段により、関心領域を設定することとしているので、軟部組織から、血管を含む血管の近傍領域の抽出を行うことができる。
【0014】
また、第5の観点の発明にかかる画像処理装置は、前記関心領域設定手段に、前記関心領域の範囲を記憶し、前記垂直断層画像情報の断層画像に近接する複数の断層画像に前記範囲を設定する再設定手段を備えることを特徴とする。
【0015】
この第5の観点の発明によれば、関心領域設定手段は、再設定手段により、関心領域の範囲を記憶し、垂直断層画像情報の断層画像に近接する複数の断層画像に前記範囲を設定することとしているので、近接する断層画像ごとに関心領域を設定する手間を省くことができる。
【0016】
また、第6の観点の発明にかかる画像処理装置は、前記平滑化手段に、前記血栓領域の濃度斑を平滑するマスク領域の大きさを有する空間フィルタを備えること特徴とする。
【0017】
この第6の観点の発明によれば、平滑化手段は、空間フィルタに、血栓領域の濃度斑を平滑するマスク領域の大きさを有することとしているので、血栓領域の濃度斑を平滑することができる。
【0018】
また、第7の観点の発明にかかる画像処理装置は、前記空間フィルタに、マスク領域の大きさを変更する変更手段を備えることを特徴とする。
この第7の観点の発明によれば、空間フィルタは、変更手段により、マスク領域の大きさを変更することとしているので、血栓領域の濃度斑に個人差が存在する場合にも、濃度斑を平滑することができる。
【0019】
また、第8の観点の発明にかかる画像処理装置は、前記確定手段により、前記辺縁を円弧で近似することを特徴とする。
この第8の観点の発明によれば、確定手段は、辺縁を円弧で近似することとしているので、概ね円形をなす血管周囲の血栓領域を、高い確度で抽出することができる。
【0020】
また、第9の観点の発明にかかるX線CTシステムは、X線ビームを被検体に照射して、前記被検体の投影データを収集し、前記投影データの画像再構成により前記被検体の断層画像情報を取得するX線CT装置と、前記X線CT装置の前記断層画像情報を転送する転送手段と、前記転送手段により転送された前記断層画像情報から、血管中の血流方向と垂直に交わる複数の垂直断層画像情報を生成する画像処理装置と、を備えるX線CTシステムであって、前記画像処理装置は、前記垂直断層画像情報の画像から前記血管を含む前記血管の近傍領域を抽出する抽出手段と、前記近傍領域の画像を平滑化する平滑化手段と、前記平滑化された画像から、指定範囲の画素値を有する領域を抽出する領域拡張手段と、前記血栓領域の辺縁を円弧状に確定する確定手段と、を備えることを特徴とする。
【0021】
この第9の観点の発明によれば、画像処理装置は、抽出手段により、垂直断層画像情報の画像から血管を含む血管の近傍領域を抽出し、平滑化手段により、近傍領域の画像を平滑化し、領域拡張手段により、平滑化された画像から、指定範囲の画素値を有する領域を抽出し、確定手段により、血栓領域の辺縁を円弧状に確定することとしているので、血栓部分に存在する濃度斑を平滑化し、血栓部の近似領域を抽出し、この近似領域の辺縁をさらに円弧状に確定して、確度の高い血栓領域を抽出することができ、ひいては、血栓面積の評価および血栓面積の時間的変化等の定量的評価を行うことができる。
【0022】
【発明の実施の形態】
以下に添付図面を参照して、この発明にかかる血栓部画像抽出方法、画像処理装置およびX線CTシステムの好適な実施の形態について説明する。なお、これにより本発明が限定されるものではない。
【0023】
まず、本実施の形態にかかる画像処理装置およびX線CTシステムの全体構成について説明する。図1は、X線CT装置のブロック(block)図を示す。図1に示すように、本装置は、走査ガントリ(gantry)2および操作コンソール(console)6を備えている。
【0024】
走査ガントリ2はX線管20を有する。X線管20から放射された図示しないX線は、コリメータ(collimator)22により、例えば扇状のX線ビーム(beam)すなわちファンビーム(fanbeam)X線となるように成形され、検出器アレイ(array)24に照射される。
【0025】
検出器アレイ24は、ファンビームX線の広がり方向にアレイ状に配列された複数のX線検出素子を有する。検出器アレイ24は、複数のX線検出素子をアレイ状に配列した、多チャンネル(channel)のX線検出器となっている。
【0026】
複数の検出器アレイ24は、全体として、円筒凹面状に湾曲したX線入射面を形成する。検出器アレイ24は、例えばシンチレータ(scintillator)とフォトダイオード(photo diode)の組み合わせによって構成される。なお、これに限られず、例えばカドミウム・テルル(CdTe)等を利用した半導体X線検出素子またはXeガス(gas)を用いる電離箱型のX線検出素子であっても良い。X線管20、コリメータ22および検出器アレイ24は、X線照射・検出装置を構成する。
【0027】
検出器アレイ24にはデータ収集部26が接続されている。データ収集部26は、検出器アレイ24の個々のX線検出素子の検出データを収集する。X線管20からのX線の照射は、X線コントローラ(controller)28によって制御される。なお、X線管20とX線コントローラ28との接続関係およびコリメータ22とコリメータコントローラ30との接続関係については図示を省略する。コリメータ22は、コリメータコントローラ30によって制御される。
【0028】
以上の、X線管20からコリメータコントローラ30までのものが、走査ガントリ2の回転部34に搭載されている。ここで、被検体1は、回転部34の中心に位置するボア(bore)29内のクレードル(cradle)31上に、横臥状態で載置される。回転部34は、回転コントローラ36により制御されつつ回転し、X線管20からX線を爆射し、検出器アレイ24において被検体1の透過X線を検出する。なお、回転部34と回転コントローラ36との接続関係については図示を省略する。
【0029】
操作コンソール6は画像再構成部60を有する。画像再構成部60は、例えばコンピュータ等によって構成される。画像再構成部60には、制御インタフェース(interface)62が接続されている。制御インタフェース62には、走査ガントリ2が接続されている。画像再構成部60は制御インタフェース62を通じて走査ガントリ2を制御する。
【0030】
走査ガントリ2内のデータ収集部26、X線コントローラ28、コリメータコントローラ30および回転コントローラ36が制御インタフェース62を通じて制御される。なお、それら各部と制御インタフェース62との個別の接続については図示を省略する。
【0031】
画像再構成部60には、また、データ収集バッファ(buffer)64が接続されている。データ収集バッファ64には、走査ガントリ2のデータ収集部26が接続されている。データ収集部26で収集されたデータがデータ収集バッファ64を通じて画像再構成部60に入力される。
【0032】
画像再構成部60は、データ収集バッファ64を通じて収集した透過X線信号すなわち投影データを用いて画像再構成を行う。画像再構成部60には、また、記憶部66が接続されている。記憶部66は、データ収集バッファ64に収集された投影データや再構成された断層画像情報および本装置の機能を実現するためのプログラム(program)等を記憶する。
【0033】
画像再構成部60には、また、表示装置68と操作装置70がそれぞれ接続されている。表示装置68は、画像再構成部60から出力される断層画像情報やその他の情報を表示する。操作装置70は、オペレータによって操作され、各種の指示や情報等を画像再構成部60に入力する。オペレータは表示装置68および操作装置70を使用してインタラクティブ(interactive)に本装置を操作する。
【0034】
画像処理装置40は、記憶部44を介して、図示しない通信手段により、操作コンソール6の記憶部66と接続されている。画像処理装置40は、操作コンソール6とは独立したコンソール上に配設されており、記憶部66に蓄積された再構成された断層画像情報等を、記憶部44に転送し、画像処理部41で画像処理を行う。ここで、画像処理部41は、例えばコンピュータ等によって構成される。
【0035】
画像処理部41には、また、表示装置42と操作装置43がそれぞれ接続されている。表示装置42は、画像処理部41から出力される断層画像情報やその他の情報を表示する。操作装置43は、オペレータによって操作され、各種の指示や情報等を画像処理部41に入力する。オペレータは表示装置42および操作装置43を使用してインタラクティブに本装置を操作する。
【0036】
図2は、画像処理装置40の機能ブロック図である。画像処理装置40は、操作コンソール6からの断層画像情報200を有し、さらに、垂直断層画像情報生成手段210、抽出手段をなす関心領域設定手段220、閾値手段231、領域拡張手段232および確定手段240、画像処理手段をなす平滑化手段230および血栓部評価手段250、等の機能を有する。
【0037】
断層画像情報200は、記憶部66から記憶部44転送された、被検体1の連続したアキシャル(axial)断面の断層画像情報である。ここで、垂直断層画像情報生成手段210は、断層画像情報200から、断層画像中の血管の血流方向と垂直に交わる垂直断層画像を算出し、生成する。これにより、この垂直断層画像内の血管では、常に血流がこの垂直断層面と垂直をなす方向に流れ、血管断面は、概ね円形となる。
【0038】
関心領域設定手段220は、後述する画像の演算処理を行う際に、演算処理を実施する画像領域を設定するものである。この設定では、表示装置42に表示された断層画像上に関心領域(ROI;Region OF Interestとも称される)を設定し、この領域内の画素に対して各種の演算を実施する。この関心領域は、円形の領域を有し、断層画像上の位置および大きさは、円の中心および半径を、操作装置43から入力することにより決定される。なお、再設定手段は、この設定された関心領域を記憶部44に保存し、異なる断層画像上に同一の関心領域を設定する。
【0039】
閾値手段231は、画像処理装置40の操作装置43から閾値を入力し、この閾値を越えるあるいは越えない画素値を判別し、この画素値が存在する関心領域を抽出し、表示装置42に表示する。閾値手段231は、断層画面の中で画素値が大きく異なる領域を識別するのに有効に用いられる。一方、断層画面の中で画素値の異なりが小さい領域を閾値手段231で識別する際には、ノイズ(noise)その他の理由により、境界領域が確定せず、的確な領域の識別を行うには困難が伴う場合がある。なお、閾値手段231は、抽出された関心領域の境界情報を出力する。
【0040】
平滑化手段230は、指定された画像領域の画像情報に平滑化を施すもので、断層画像情報のある画素を中心にマスク(mask)と呼ばれる近傍領域を設定し、このマスク内のすべての画素値を加算あるいは重み付け加算して、新たな画素値とする。このマスクの大きさは、後述する血栓部の画像に生じる濃度斑を平滑化する大きさに実験的に決められ、概ね3×3の画素マトリクス(matrix)の大きさを有する。また、平滑化手段230は、被検体1の個人差を考慮し、オペレータが画像処理装置40の操作装置43から、画素マトリクスの大きさを変更する変更手段も有する。
【0041】
領域拡張手段232は、関心領域内にある画像の画素値が、オペレータが指定した上限および下限の範囲内にある画素を抽出し、画像領域の抽出を行う。領域拡張手段232では、画素値が、ある上限および下限を有する有限の範囲内にある、中間的なCT値を有する組織部分のみを、抽出領域が画像領域内で分離することなく抽出するのに用いられる。一方、不均一な組織、例えば多くの濃度斑が存在する血栓部の様な場合には、組織の画素値が大きな範囲に渡って点在するため、組織部分のみを抽出するには困難が伴う。
【0042】
確定手段240は、表示装置42の断層画面上に設定された閉じた境界領域に対して、円弧近似を行い、例えば境界領域との偏差が最小となる円形の境界を生成する。これにより、閉じた境界領域に最も適合する円形領域を抽出することができる。
【0043】
血栓部評価手段250は、抽出された血栓部領域に基づいて、血栓部領域の面積、あるいは血栓部領域に含まれる血管部分の面積との比率、さらには、時系列的に取得された血栓部領域の面積から、血栓部の時間的変化等の定量的情報の算出を行う。
【0044】
次に、本実施の形態にかかる画像処理部41の動作について説明する前に、操作コンソール6および走査ガントリ2による画像収集について概要を説明する。まず、図3のフローチャートを用いて、X線CTシステムによる被検体の撮像を示す。オペレータは、被検体1を図1に示した走査ガントリ2のクレードル31上に載置する(ステップS301)。そして、オペレータは、被検体1に造影剤を投与する(ステップS302)。この際、造影剤は、主としてヨード(iodin)系のものが、例えば経静脈的に投与される。
【0045】
その後、オペレータは、被検体の撮像を行い(ステップS303)、被検体の投影データを取得する。そして、画像再構成部60において、投影データから、被検体の断層画像情報を生成し、記憶部66に保存する。そして、記憶部66の断層画像情報を、図示しない転送手段により、画像処理部41の記憶部44に転送する(ステップS304)。その後、画像処理部41は、転送された断層画像情報に、血栓部抽出処理を行い(ステップS305)、本処理を終了する。
【0046】
つづいて、図4のフローチャートを用いて、画像処理部41による血栓部抽出処理を、具体的に説明する。オペレータは、操作コンソール6から転送された断層画像情報200から、断層画像内に描出された患部、すなわち動脈瘤が存在する近傍の血管と直交する垂直断層面の画像情報を垂直断層画像生成手段210により生成する(ステップS401)。図5(A)に、X線CT装置を用いて腹部を撮像した際の、腹部大動脈の垂直断層画像の一例を示した。垂直断層画像では、腹部大動脈の血管は、概ね円形に描出される。
【0047】
その後、オペレータは、垂直断層画像から、骨部領域および背景領域を、閾値手段231により選別し、血管およびその周囲に存在する血栓部を含む軟部組織のみの関心領域を抽出する(ステップS402)。図5(B)に軟部組織のみからなる関心領域を図示した。背景および脊椎の骨部分が、関心領域から除外される。
【0048】
その後、オペレータは、軟部組織のみの関心領域から血管領域を抽出する(ステップS403)。造影剤の投与により、血管領域は、高いCT値で描出されるので閾値手段231で抽出することができる。また、造影剤の効果は、一過性であるため、所定時間経過後の垂直断面の画像情報の差分を取ることにより、血管領域のみを抽出することも容易にできる。なお、この抽出により、後述する血管領域を省いた近傍領域の平滑化および血栓部のみの抽出が円滑に行われる。
【0049】
その後、オペレータは、ステップS403で抽出された血管領域を含む周囲に、近傍領域を関心領域設定手段220により設定する(ステップS404)。この設定では、血管周囲に付着した血栓部を充分に含むように、他方、ステップS402で抽出された軟部組織の関心領域に含まれるような大きさの円形の関心領域が設定される。図5(C)に軟部組織内の腹部大動脈周囲に設定された近傍領域を図示した。
【0050】
図6に、この近傍領域のみを拡大して示した。この近傍領域は、図6(A)に示すように、血管、血栓部および軟部組織よりなる。血管周囲に存在する血栓部には、図6(A)中の黒点で示したような濃度斑が存在する。また、血管内部には、血流部が存在するが、ステップS403の血管領域の抽出により、近傍領域から除外される。
【0051】
その後、オペレータは、ステッS404で設定された近傍領域に対して、平滑化手段230を用いて、平滑化を行う(ステップS405)。この平滑化により、血栓部領域に存在する濃度斑は、平滑化され、概ね均一な画像情報となる。このため、閾値手段231あるいは領域拡張手段232を実施することができる。図6(B)に平滑化された関心領域を模式的に示した。血栓部の濃度斑は、平滑化により消滅し、周囲の軟部組織に近似した一様な画素領域となる。
【0052】
その後、オペレータは、平滑化された血栓部とその周辺組織の境界が区別できるかどうかを判定し(ステップS406)、区別できる場合には(ステップS406否定)、閾値手段231の閾値を設定し、周辺組織の除去を行う(ステップS407)。また、区別できない場合には(ステップS406肯定)、ステップS407を実行せず次のステップに進む。
【0053】
その後、オペレータは、平滑化された関心領域に対して、領域拡張手段232を用いて、血栓部の近似領域のみの抽出を行う(ステップS408)。これにより、血栓部の近似領域が抽出される。図6(C)に、領域拡張手段232により抽出された血栓部の近似領域を模式的に示した。ステップS405で行った平滑化処理等の理由により、この近似領域の辺縁部は、周囲の軟部組織との境界が確定せず、ばらついた状態となる。
【0054】
その後、オペレータは、近似領域の辺縁の境界に対して、確定手段240を用いて、円弧近似を行い(ステップS409)、円形の血栓部領域を確定する。図6(D)に、確定手段240により抽出された円形の血栓部領域を模式的に示した。円弧近似を行う根拠は、血管の周囲に付着する血栓部の断面は、円形をなすという、実験的および経験的な事実に基づいている。
【0055】
その後、近傍に存在する異なる垂直断層面から血栓部の抽出を行うかどうかを判定し(ステップS410)、血栓部の抽出を行う場合には(ステップS410肯定)、ステップS401で生成された垂直断層面の中から近傍に存在する垂直断層面を選択し、表示する(ステップS411)。そして、ステップS404に移行し、表示画像に関心領域を設定する。この際、再設定手段を用いて、既に設定した関心領域をそのまま用いることも出来る。
【0056】
また、血栓部の抽出を行わない場合には(ステップS410否定)、血栓部の抽出処理を終了し、取得された血栓部情報から血栓部の評価を行う(ステップS412)。この血栓部の評価では、抽出された血栓部領域の面積あるいは血管部領域との面積比等が算出される。さらに、時系列をもって採取された血栓部領域の面積変化、手術前および手術後の血栓部領域の変化等も算出される。
【0057】
上述してきたように、本実施の形態では、血管と直交する垂直断層面内の血管を含む近傍領域を抽出し、この近傍領域に平滑化処理を行い、平滑化された関心領域から領域拡張処理により血栓部の近似領域を抽出し、この血栓部の近似領域から円弧処理により、円弧領域を抽出することとしているので、領域拡張処理を用いて血栓部の近似領域を抽出し、この近似領域から確度の高い円弧状の血栓領域を抽出し、ひいては、血栓面積の評価および血栓面積の時間的変化等の定量的評価を行うことができる。
【0058】
【発明の効果】
以上説明したように、本発明によれば、X線CT装置により被検体の複数の断層画像情報を取得し、この断層画像情報から血管中の血流方向と垂直に交わる複数の垂直断層画像情報を生成し、この垂直断層画像情報の画像から血管を含む血管の近傍領域を抽出し、この近傍領域の画像を平滑化し、領域拡張手段により、この平滑化された画像から、血管に付着する血栓部の近似領域を抽出し、この近似領域の辺縁を円弧状に確定することとしているので、血栓部に存在する濃度斑を平滑化し、血栓部の近似領域を抽出し、この抽出部分の辺縁をさらに円弧状に確定して、確度の高い血栓領域を抽出することができ、ひいては、血栓面積の評価および血栓面積の時間的変化等の定量的評価を行うことができる。
【図面の簡単な説明】
【図1】X線CTシステムの全体構成を示すブロック図である。
【図2】実施の形態の画像処理装置の機能ブロック図である。
【図3】実施の形態の断層画像情報の取得を示すフローチャートである。
【図4】実施の形態の血栓部抽出処理を示すフローチャートである。
【図5】実施の形態の垂直断層画像からの近傍領域の抽出を示す図である。
【図6】実施の形態の近傍領域の画像処理を示す模式図である。
【符号の説明】
1 被検体
2 走査ガントリ
6 操作コンソール
20 X線管
22 コリメータ
24 検出器アレイ
26 データ収集部
28 X線コントローラ
30 コリメータコントローラ
31 クレードル
34 回転部
36 回転コントローラ
40 画像処理装置
41 画像処理部
42 表示装置
43 操作装置
44 記憶部
60 画像再構成部
62 制御インタフェース
64 データ収集バッファ
66 記憶部
68 表示装置
70 操作装置
200 断層画像情報
210 断層画像生成手段
210 垂直断層画像生成手段
220 関心領域設定手段
230 平滑化手段
231 閾値手段
232 領域拡張手段
240 確定手段
250 血栓部評価手段
[0001]
TECHNICAL FIELD OF THE INVENTION
The present invention relates to an image processing apparatus of an X-ray CT apparatus that performs imaging using a contrast agent.
[0002]
[Prior art]
In recent years, X-ray CT imaging using an angiographic agent is performed for the purpose of imaging an aneurysm. An aneurysm consists of a blood vessel and a thrombus portion which adheres to the wall of the blood vessel and in which the blood solidifies in a jelly state. Here, in order to predict the progress of the aneurysm, it is important to obtain information on the blood vessel diameter of the aneurysm and the size of the thrombus. In imaging using a contrast agent, the blood vessel portion of the aneurysm is clearly depicted by the blood contrast agent.
[0003]
[Problems to be solved by the invention]
However, according to the above-mentioned prior art, the thrombus portion adhering to the blood vessel wall has no contrast effect, and it is difficult to identify the thrombus from the surrounding tissue. That is, the contrast agent did not flow into the thrombus portion, so that no contrast effect was observed.
[0004]
In particular, attempts have been made to extract and measure information on a thrombus portion that is difficult to identify using an image processing technique. However, light and shade of brightness called density unevenness exist in the thrombus portion, which is a factor that hinders measurement of the thrombus portion using image processing technology.
[0005]
From these facts, how to realize a thrombus portion image extraction method, an image processing apparatus, and an X-ray CT system that can easily and highly accurately measure the size of a thrombus of an aneurysm using image processing technology. It becomes important.
[0006]
The present invention has been made in order to solve the above-described problems of the related art, and the thrombus size of a thrombus of an aneurysm can be measured easily and with high accuracy using an image processing technique. An object is to provide an extraction method, an image processing device, and an X-ray CT system.
[0007]
[Means for Solving the Problems]
In order to solve the above-described problems and achieve the object, a thrombus portion image extracting method according to a first aspect of the present invention obtains a plurality of pieces of tomographic image information of a subject using an X-ray CT apparatus, and obtains the tomographic image information. Generates a plurality of pieces of vertical tomographic image information that intersects perpendicularly with the direction of blood flow in the blood vessel, extracts a region near the blood vessel including the blood vessel from the image of the vertical tomographic image information, and smoothes the image of the nearby region. Extracting an approximate region of a thrombus portion adhering to the blood vessel from the smoothed image, and determining an edge of the thrombus region in an arc shape.
[0008]
According to the invention of the first aspect, a plurality of tomographic image information of the subject is acquired by the X-ray CT apparatus, and a plurality of pieces of vertical tomographic image information that intersects perpendicularly with the blood flow direction in the blood vessel are obtained from the tomographic image information. Generate and extract a region near a blood vessel including a blood vessel from the image of the vertical tomographic image information, smooth the image of the nearby region, and extract an approximate region of a thrombus portion attached to the blood vessel from the smoothed image. Since the edge of this approximate area is determined to be an arc, the density unevenness existing in the thrombus is smoothed, the approximate area of the thrombus is extracted, and the edge of the extracted part is further determined to be an arc. As a result, a highly accurate thrombus region can be extracted, and as a result, a thrombus area can be evaluated and a quantitative evaluation such as a temporal change in the thrombus area can be performed.
[0009]
Further, the image processing apparatus according to the second aspect of the present invention acquires a plurality of tomographic image information of a subject from an X-ray CT apparatus, and from the tomographic image information, a plurality of perpendicular to the blood flow direction in a blood vessel. An image processing apparatus for generating tomographic image information, comprising: an extracting unit configured to extract a region near the blood vessel including the blood vessel from an image of the vertical tomographic image information; and a smoothing unit configured to smooth an image of the nearby region. A region expansion unit that extracts a region having a pixel value in a specified range from the smoothed image; and a determination unit that determines the edge of the approximate region in an arc shape.
[0010]
According to the invention according to the second aspect, the extraction unit extracts the vicinity region of the blood vessel including the blood vessel from the vertical tomographic image information, smoothes the image of the vicinity region by the averaging unit, and designates the designated region by the region expansion unit. Since the region having the pixel value of the range is extracted and the edge of the thrombus region is determined to be an arc by the determination means, the density unevenness existing in the thrombus portion is smoothed, and the approximate region of the thrombus portion is extracted. It is possible to further determine the edge of this approximate area in an arc shape and extract a highly accurate thrombus area, and thus to perform thrombus area evaluation and quantitative evaluation such as temporal change in thrombus area. .
[0011]
An image processing apparatus according to a third aspect of the present invention is characterized in that the extraction means includes a threshold means for extracting an image area having a pixel value exceeding a threshold or a pixel value not exceeding a threshold.
[0012]
According to the invention of the third aspect, the extracting means extracts the image area having the pixel value exceeding the threshold value or the pixel value not exceeding the threshold value by the threshold value means. By extracting a background region, image information consisting of only a soft tissue portion can be generated.
[0013]
Further, an image processing apparatus according to a fourth aspect of the invention is characterized in that the extraction means includes a region of interest setting means for setting a region of interest.
According to the fourth aspect of the invention, since the region of interest is set by the region of interest setting unit, the extraction unit can extract a region near a blood vessel including a blood vessel from the soft tissue.
[0014]
The image processing apparatus according to a fifth aspect of the present invention stores the range of the region of interest in the region of interest setting means, and stores the range in a plurality of tomographic images close to the tomographic image of the vertical tomographic image information. It is characterized by comprising resetting means for setting.
[0015]
According to the fifth aspect of the invention, the region of interest setting means stores the range of the region of interest by the resetting means and sets the range in a plurality of tomographic images close to the tomographic image of the vertical tomographic image information. Therefore, it is not necessary to set a region of interest for each adjacent tomographic image.
[0016]
Further, the image processing apparatus according to the invention of a sixth aspect is characterized in that the smoothing means includes a spatial filter having a size of a mask area for smoothing density unevenness of the thrombus area.
[0017]
According to the sixth aspect of the invention, since the smoothing means has the spatial filter having the size of the mask area for smoothing the density unevenness in the thrombus area, the smoothing means can smooth the density unevenness in the thrombus area. it can.
[0018]
An image processing apparatus according to a seventh aspect of the present invention is characterized in that the spatial filter includes a changing unit that changes a size of a mask area.
According to the seventh aspect of the invention, since the spatial filter changes the size of the mask area by the changing unit, even when individual differences exist in the density unevenness of the thrombus area, the spatial filter can reduce the density unevenness. Can be smoothed.
[0019]
An image processing apparatus according to an eighth aspect of the present invention is characterized in that the edge is approximated by a circular arc by the determination means.
According to the eighth aspect of the invention, since the determining means approximates the periphery with a circular arc, it is possible to extract a thrombus region around a blood vessel having a substantially circular shape with high accuracy.
[0020]
An X-ray CT system according to a ninth aspect of the present invention irradiates a subject with an X-ray beam, collects projection data of the subject, and reconstructs a tomographic image of the subject by image reconstruction of the projection data. An X-ray CT apparatus that acquires image information; a transfer unit that transfers the tomographic image information of the X-ray CT apparatus; An image processing apparatus for generating a plurality of intersecting vertical tomographic image information, wherein the image processing apparatus extracts a region near the blood vessel including the blood vessel from an image of the vertical tomographic image information. Extracting means, a smoothing means for smoothing the image of the neighboring area, an area expanding means for extracting an area having a pixel value in a specified range from the smoothed image, and an edge of the thrombus area. Arc shape And determining means for determining, characterized in that it comprises a.
[0021]
According to the ninth aspect, in the image processing apparatus, the extracting unit extracts a region near the blood vessel including the blood vessel from the image of the vertical tomographic image information, and the smoothing unit smoothes the image of the nearby region. Since the region expanding means extracts a region having a pixel value in a specified range from the smoothed image, and the deciding device decides the periphery of the thrombus region in an arc shape, the region exists in the thrombus portion. The density unevenness is smoothed, an approximate region of the thrombus portion is extracted, and the edge of the approximate region is further determined in an arc shape, so that a thrombus region with high accuracy can be extracted.Thus, evaluation of thrombus area and thrombus It is possible to perform a quantitative evaluation such as a temporal change of the area.
[0022]
BEST MODE FOR CARRYING OUT THE INVENTION
DESCRIPTION OF EMBODIMENTS Preferred embodiments of a thrombus portion image extraction method, an image processing device, and an X-ray CT system according to the present invention will be described below with reference to the accompanying drawings. Note that the present invention is not limited to this.
[0023]
First, an overall configuration of an image processing apparatus and an X-ray CT system according to the present embodiment will be described. FIG. 1 shows a block diagram of the X-ray CT apparatus. As shown in FIG. 1, the apparatus includes a scanning gantry 2 and an operation console 6.
[0024]
The scanning gantry 2 has an X-ray tube 20. An X-ray (not shown) emitted from the X-ray tube 20 is shaped by a collimator 22 into, for example, a fan-shaped X-ray beam (beam), that is, a fan beam X-ray, and a detector array (array). ) 24.
[0025]
The detector array 24 has a plurality of X-ray detection elements arranged in an array in the direction in which the fan beam X-rays spread. The detector array 24 is a multi-channel (channel) X-ray detector in which a plurality of X-ray detection elements are arranged in an array.
[0026]
The plurality of detector arrays 24 form an X-ray incident surface curved as a cylindrical concave surface as a whole. The detector array 24 is composed of, for example, a combination of a scintillator and a photodiode. The present invention is not limited to this, and may be, for example, a semiconductor X-ray detection element using cadmium tellurium (CdTe) or the like or an ionization box type X-ray detection element using Xe gas (gas). The X-ray tube 20, the collimator 22, and the detector array 24 constitute an X-ray irradiation / detection device.
[0027]
A data collection unit 26 is connected to the detector array 24. The data collection unit 26 collects detection data of individual X-ray detection elements of the detector array 24. The irradiation of X-rays from the X-ray tube 20 is controlled by an X-ray controller 28. The connection relationship between the X-ray tube 20 and the X-ray controller 28 and the connection relationship between the collimator 22 and the collimator controller 30 are not shown. The collimator 22 is controlled by a collimator controller 30.
[0028]
The components from the X-ray tube 20 to the collimator controller 30 are mounted on the rotating unit 34 of the scanning gantry 2. Here, the subject 1 is placed on a cradle 31 in a bore 29 located at the center of the rotating unit 34 in a lying state. The rotation unit 34 rotates while being controlled by the rotation controller 36, emits X-rays from the X-ray tube 20, and detects transmitted X-rays of the subject 1 in the detector array 24. The illustration of the connection relationship between the rotation unit 34 and the rotation controller 36 is omitted.
[0029]
The operation console 6 has an image reconstruction unit 60. The image reconstruction unit 60 is configured by, for example, a computer. The image reconstruction unit 60 is connected to a control interface (interface) 62. The scanning gantry 2 is connected to the control interface 62. The image reconstruction unit 60 controls the scanning gantry 2 through the control interface 62.
[0030]
The data acquisition unit 26, X-ray controller 28, collimator controller 30, and rotation controller 36 in the scanning gantry 2 are controlled through a control interface 62. It should be noted that illustration of individual connections between these units and the control interface 62 is omitted.
[0031]
The image reconstruction unit 60 is also connected to a data collection buffer (buffer) 64. The data collection buffer 26 is connected to the data collection unit 26 of the scanning gantry 2. The data collected by the data collection unit 26 is input to the image reconstruction unit 60 via the data collection buffer 64.
[0032]
The image reconstruction unit 60 performs image reconstruction using the transmitted X-ray signal, that is, the projection data collected through the data collection buffer 64. The storage unit 66 is also connected to the image reconstruction unit 60. The storage unit 66 stores the projection data collected in the data collection buffer 64, the reconstructed tomographic image information, and a program for realizing the functions of the present apparatus.
[0033]
A display device 68 and an operation device 70 are connected to the image reconstruction unit 60, respectively. The display device 68 displays tomographic image information and other information output from the image reconstruction unit 60. The operation device 70 is operated by an operator, and inputs various instructions and information to the image reconstruction unit 60. The operator uses the display device 68 and the operation device 70 to interactively operate the apparatus.
[0034]
The image processing device 40 is connected to the storage unit 66 of the operation console 6 via the storage unit 44 by communication means (not shown). The image processing device 40 is provided on a console independent of the operation console 6, transfers the reconstructed tomographic image information and the like stored in the storage unit 66 to the storage unit 44, and transmits the image processing unit 41. Performs image processing. Here, the image processing unit 41 is configured by, for example, a computer.
[0035]
A display device 42 and an operation device 43 are connected to the image processing unit 41, respectively. The display device 42 displays tomographic image information output from the image processing unit 41 and other information. The operation device 43 is operated by an operator and inputs various instructions and information to the image processing unit 41. The operator uses the display device 42 and the operation device 43 to interactively operate the apparatus.
[0036]
FIG. 2 is a functional block diagram of the image processing device 40. The image processing device 40 has tomographic image information 200 from the operation console 6, and further has a vertical tomographic image information generating means 210, a region of interest setting means 220 serving as an extracting means, a threshold value means 231, an area expanding means 232, and a fixing means. 240, a function of a smoothing means 230 serving as an image processing means, and a thrombus evaluation means 250.
[0037]
The tomographic image information 200 is tomographic image information of a continuous axial section of the subject 1 transferred from the storage unit 66 to the storage unit 44. Here, the vertical tomographic image information generating means 210 calculates and generates a vertical tomographic image that intersects perpendicularly with the blood flow direction of the blood vessel in the tomographic image from the tomographic image information 200. Thus, in the blood vessel in the vertical tomographic image, the blood flow always flows in a direction perpendicular to the vertical tomographic plane, and the blood vessel cross-section is substantially circular.
[0038]
The region-of-interest setting means 220 sets an image region to be subjected to the arithmetic processing when performing the arithmetic processing of the image described later. In this setting, a region of interest (ROI; also referred to as Region Of Interest) is set on the tomographic image displayed on the display device 42, and various calculations are performed on pixels in this region. The region of interest has a circular region, and the position and size on the tomographic image are determined by inputting the center and radius of the circle from the operation device 43. The resetting unit stores the set region of interest in the storage unit 44 and sets the same region of interest on different tomographic images.
[0039]
The threshold means 231 receives a threshold value from the operation device 43 of the image processing device 40, determines a pixel value exceeding or not exceeding the threshold value, extracts a region of interest in which the pixel value exists, and displays the region of interest on the display device 42. . The threshold means 231 is effectively used to identify a region in the tomographic screen where the pixel values differ greatly. On the other hand, when an area having a small difference in pixel value in the tomographic screen is identified by the threshold unit 231, the boundary area is not determined due to noise or other reasons, and it is necessary to identify an accurate area. Difficulties may accompany. The threshold unit 231 outputs the extracted boundary information of the region of interest.
[0040]
The smoothing means 230 performs smoothing on the image information of the designated image area. The smoothing means 230 sets a neighborhood area called a mask around a certain pixel having tomographic image information, and sets all pixels in the mask. The values are added or weighted to obtain a new pixel value. The size of this mask is experimentally determined to be a size for smoothing density unevenness generated in an image of the thrombus, which will be described later, and has a size of approximately 3 × 3 pixel matrix. The smoothing unit 230 also includes a changing unit that allows the operator to change the size of the pixel matrix from the operation device 43 of the image processing device 40 in consideration of the individual difference of the subject 1.
[0041]
The region expanding unit 232 extracts pixels whose pixel values of the image within the region of interest are within the range of the upper limit and the lower limit specified by the operator, and extracts the image region. The region expanding means 232 extracts only a tissue portion having an intermediate CT value whose pixel value is within a finite range having a certain upper limit and a lower limit without extracting the extracted region within the image region. Used. On the other hand, in the case of a non-uniform tissue, for example, a thrombus where many density spots are present, it is difficult to extract only the tissue part because the pixel values of the tissue are scattered over a large range. .
[0042]
The determination means 240 performs circular arc approximation on the closed boundary area set on the tomographic screen of the display device 42, and generates, for example, a circular boundary having a minimum deviation from the boundary area. As a result, it is possible to extract a circular region that best fits the closed boundary region.
[0043]
Based on the extracted thrombus region, the thrombus evaluation unit 250 calculates the ratio of the area of the thrombus region or the ratio of the area of the blood vessel portion included in the thrombus region to the thrombus region obtained in time series. From the area of the area, quantitative information such as a temporal change of the thrombus is calculated.
[0044]
Next, before describing the operation of the image processing unit 41 according to the present embodiment, an outline of image collection by the operation console 6 and the scanning gantry 2 will be described. First, imaging of a subject by the X-ray CT system will be described with reference to the flowchart of FIG. The operator places the subject 1 on the cradle 31 of the scanning gantry 2 shown in FIG. 1 (Step S301). Then, the operator administers the contrast agent to the subject 1 (Step S302). At this time, the iodine-based contrast agent is mainly administered, for example, intravenously.
[0045]
Thereafter, the operator performs imaging of the subject (Step S303), and acquires projection data of the subject. Then, in the image reconstruction unit 60, tomographic image information of the subject is generated from the projection data and stored in the storage unit 66. Then, the tomographic image information in the storage unit 66 is transferred to the storage unit 44 of the image processing unit 41 by a transfer unit (not shown) (step S304). After that, the image processing unit 41 performs a thrombus extraction process on the transferred tomographic image information (step S305), and ends this process.
[0046]
Next, the thrombus extraction process performed by the image processing unit 41 will be specifically described with reference to the flowchart in FIG. From the tomographic image information 200 transferred from the operation console 6, the operator converts the image information of the diseased part drawn in the tomographic image, that is, the vertical tomographic plane orthogonal to the blood vessel in the vicinity of the aneurysm, into the vertical tomographic image generating means 210. (Step S401). FIG. 5A shows an example of a vertical tomographic image of the abdominal aorta when imaging the abdomen using an X-ray CT apparatus. In the vertical tomographic image, the blood vessels of the abdominal aorta are depicted in a substantially circular shape.
[0047]
Thereafter, the operator selects a bone region and a background region from the vertical tomographic image by the threshold means 231 and extracts a region of interest only of a soft tissue including a blood vessel and a thrombus existing around the blood vessel (step S402). FIG. 5B shows a region of interest consisting only of soft tissue. The background and bone parts of the spine are excluded from the region of interest.
[0048]
Thereafter, the operator extracts a blood vessel region from the region of interest only of the soft tissue (step S403). By administering the contrast agent, the blood vessel region is drawn with a high CT value, and can be extracted by the threshold means 231. In addition, since the effect of the contrast agent is transient, it is easy to extract only the blood vessel region by taking the difference between the image information of the vertical cross section after a predetermined time has elapsed. By this extraction, smoothing of a nearby region excluding a blood vessel region described later and extraction of only the thrombus portion are smoothly performed.
[0049]
Thereafter, the operator sets a nearby region by the region-of-interest setting means 220 around the blood vessel region extracted in step S403 (step S404). In this setting, a circular region of interest is set so as to sufficiently include the thrombus attached to the periphery of the blood vessel and to be included in the region of interest of the soft tissue extracted in step S402. FIG. 5 (C) illustrates a neighboring region set around the abdominal aorta in the soft tissue.
[0050]
FIG. 6 shows an enlarged view of only this vicinity region. As shown in FIG. 6 (A), this nearby region is composed of blood vessels, thrombus, and soft tissue. In the thrombus around blood vessels, there are density spots as indicated by black dots in FIG. Although a blood flow part exists inside the blood vessel, the blood flow part is excluded from the neighboring area by extracting the blood vessel area in step S403.
[0051]
Thereafter, the operator performs smoothing on the neighboring area set in step S404 using the smoothing means 230 (step S405). By this smoothing, the density unevenness existing in the thrombus region is smoothed, and becomes substantially uniform image information. Therefore, the threshold means 231 or the area expanding means 232 can be implemented. FIG. 6B schematically shows the smoothed region of interest. The density unevenness of the thrombus disappears by the smoothing, and becomes a uniform pixel area similar to the surrounding soft tissue.
[0052]
Thereafter, the operator determines whether or not the boundary between the smoothed thrombus portion and the surrounding tissue can be distinguished (step S406). If the boundary can be distinguished (step S406 negative), the operator sets the threshold value of the threshold means 231; The surrounding tissue is removed (step S407). If it cannot be distinguished (Yes at Step S406), the process proceeds to the next step without executing Step S407.
[0053]
Thereafter, the operator extracts only the approximate region of the thrombus portion from the smoothed region of interest by using the region expanding means 232 (step S408). Thereby, the approximate region of the thrombus portion is extracted. FIG. 6C schematically shows the approximate region of the thrombus extracted by the region expanding means 232. Due to the smoothing process or the like performed in step S405, the border of this approximate region is in a state in which the boundary with the surrounding soft tissue is not determined and varies.
[0054]
Thereafter, the operator performs arc approximation on the boundary of the edge of the approximate region using the determining means 240 (step S409), and determines the circular thrombus region. FIG. 6D schematically shows the circular thrombus region extracted by the determination means 240. The basis for performing the arc approximation is based on experimental and empirical facts that the cross section of the thrombus adhering around the blood vessel is circular.
[0055]
Thereafter, it is determined whether or not to extract a thrombus from a nearby different vertical tomographic plane (step S410). If a thrombus is to be extracted (Yes at step S410), the vertical tomogram generated at step S401 is determined. A nearby vertical tomographic plane is selected from the planes and displayed (step S411). Then, the process shifts to step S404 to set a region of interest in the display image. At this time, the already set region of interest can be used as it is by using the resetting means.
[0056]
If the thrombus portion is not extracted (No at Step S410), the thrombus portion extraction process is terminated, and the thrombus portion is evaluated from the acquired thrombus portion information (Step S412). In the evaluation of the thrombus portion, the area of the extracted thrombus region or the area ratio with the blood vessel region is calculated. Further, a change in the area of the thrombus region collected in time series, a change in the thrombus region before and after the operation, and the like are also calculated.
[0057]
As described above, in the present embodiment, a nearby region including a blood vessel in a vertical tomographic plane orthogonal to the blood vessel is extracted, a smoothing process is performed on this nearby region, and a region expanding process is performed on the smoothed region of interest. The approximate region of the thrombus is extracted by the following, and the circular region is extracted by the arc processing from the approximate region of the thrombus.Therefore, the approximate region of the thrombus is extracted using the region expanding process, and the approximate region is extracted from the approximate region. It is possible to extract an arc-shaped thrombus region with high accuracy, and thus to perform thrombus area evaluation and quantitative evaluation such as temporal change in thrombus area.
[0058]
【The invention's effect】
As described above, according to the present invention, a plurality of pieces of tomographic image information of a subject are acquired by an X-ray CT apparatus, and a plurality of pieces of vertical tomographic image information that intersects perpendicularly with a blood flow direction in a blood vessel are obtained from the tomographic image information. Is extracted from the image of the vertical tomographic image information, a region near the blood vessel including the blood vessel is extracted, the image of the nearby region is smoothed, and a thrombus adhering to the blood vessel is extracted from the smoothed image by the region expanding means. Since the approximate region of the thrombus is extracted and the edge of this approximate region is determined in an arc shape, the density unevenness present in the thrombus is smoothed, and the approximate region of the thrombus is extracted. The edge is further determined to be an arc shape, and a thrombus region with high accuracy can be extracted. As a result, a thrombus area can be evaluated, and a quantitative evaluation such as a temporal change in the thrombus area can be performed.
[Brief description of the drawings]
FIG. 1 is a block diagram showing the overall configuration of an X-ray CT system.
FIG. 2 is a functional block diagram of the image processing apparatus according to the embodiment;
FIG. 3 is a flowchart illustrating acquisition of tomographic image information according to the embodiment;
FIG. 4 is a flowchart illustrating a thrombus extraction process according to the embodiment;
FIG. 5 is a diagram showing extraction of a nearby region from a vertical tomographic image according to the embodiment;
FIG. 6 is a schematic diagram illustrating image processing of a nearby area according to the embodiment;
[Explanation of symbols]
1 Subject
2 scanning gantry
6 Operation console
20 X-ray tube
22 Collimator
24 Detector array
26 Data Collection Unit
28 X-ray controller
30 Collimator controller
31 cradle
34 Rotating part
36 Rotation controller
40 Image processing device
41 Image processing unit
42 Display
43 Operating device
44 Memory
60 Image reconstruction unit
62 Control interface
64 data collection buffer
66 storage unit
68 Display device
70 Operating device
200 tomographic image information
210 tomographic image generation means
210 Vertical tomographic image generation means
220 region of interest setting means
230 smoothing means
231 threshold means
232 area expansion means
240 Confirmation means
250 thrombus evaluation means

Claims (9)

X線CT装置により被検体の複数の断層画像情報を取得し、前記断層画像情報から血管中の血流方向と垂直に交わる複数の垂直断層画像情報を生成し、
前記垂直断層画像情報の画像から前記血管を含む前記血管の近傍領域を抽出し、
前記近傍領域の画像を平滑化し、
前記平滑化された画像から前記血管に付着する血栓部の近似領域を抽出し、
前記近似領域の辺縁を円弧状に確定すること、
を特徴とする血栓部画像抽出方法。
A plurality of tomographic image information of the subject is acquired by an X-ray CT apparatus, and a plurality of vertical tomographic image information that intersects perpendicularly with a blood flow direction in a blood vessel is generated from the tomographic image information;
Extracting the vicinity region of the blood vessel including the blood vessel from the image of the vertical tomographic image information,
Smoothing the image of the neighborhood,
Extract an approximate region of a thrombus attached to the blood vessel from the smoothed image,
Determining the edge of the approximate region in an arc shape,
A method for extracting an image of a thrombus portion, characterized in that:
X線CT装置から被検体の複数の断層画像情報を取得し、前記断層画像情報から血管中の血流方向と垂直に交わる複数の垂直断層画像情報を生成する画像処理装置であって、
前記垂直断層画像情報の画像から前記血管を含む前記血管の近傍領域を抽出する抽出手段と、
前記近傍領域の画像を平滑化する平滑化手段と、
前記平滑化された画像から、指定範囲の画素値を有する領域を抽出する領域拡張手段と、
前記近似領域の辺縁を円弧状に確定する確定手段と、
を備えることを特徴とする画像処理装置。
An image processing apparatus that acquires a plurality of pieces of tomographic image information of a subject from an X-ray CT apparatus and generates a plurality of pieces of vertical tomographic image information that intersects perpendicularly with a blood flow direction in a blood vessel from the tomographic image information,
Extracting means for extracting a region near the blood vessel including the blood vessel from the image of the vertical tomographic image information,
Smoothing means for smoothing the image of the neighboring area,
Region expanding means for extracting a region having a pixel value in a specified range from the smoothed image;
Determining means for determining the edge of the approximate area in an arc shape,
An image processing apparatus comprising:
前記抽出手段は、閾値を越える画素値、あるいは閾値を越えない画素値の画像領域を抽出する閾値手段を備えることを特徴とする請求項2に記載の画像処理装置。The image processing apparatus according to claim 2, wherein the extraction unit includes a threshold unit that extracts an image area having a pixel value exceeding a threshold value or a pixel value not exceeding the threshold value. 前記抽出手段は、関心領域を設定する関心領域設定手段を備えることを特徴とする請求項2ないし3のいずれか1つに記載の画像処理装置。The image processing apparatus according to claim 2, wherein the extraction unit includes a region of interest setting unit that sets a region of interest. 前記関心領域設定手段は、前記関心領域の範囲を記憶し、前記垂直断層画像情報の断層画像に近接する複数の断層画像に前記範囲を設定する再設定手段を備えることを特徴とする請求項4に記載の画像処理装置。5. The apparatus according to claim 4, wherein the region of interest setting unit stores a range of the region of interest and resets the region to a plurality of tomographic images close to a tomographic image of the vertical tomographic image information. An image processing apparatus according to claim 1. 前記平滑化手段は、前記血栓領域の濃度斑を平滑するマスク領域の大きさを有する空間フィルタを備えること特徴とする請求項2ないし5のいずれか1つに記載の画像処理装置。The image processing apparatus according to any one of claims 2 to 5, wherein the smoothing unit includes a spatial filter having a size of a mask area for smoothing density unevenness in the thrombus area. 前記空間フィルタは、マスク領域の大きさを変更する変更手段を備えることを特徴とする請求項6に記載の画像処理装置。The image processing apparatus according to claim 6, wherein the spatial filter includes a changing unit that changes a size of a mask area. 前記確定手段は、前記辺縁を円弧で近似することを特徴とする請求項2ないし7のいずれか1つに記載の画像処理装置。The image processing apparatus according to claim 2, wherein the determination unit approximates the edge with a circular arc. X線ビームを被検体に照射して、前記被検体の投影データを収集し、前記投影データの画像再構成により前記被検体の断層画像情報を取得するX線CT装置と、
前記X線CT装置の前記断層画像情報を転送する転送手段と、
前記転送手段により転送された前記断層画像情報から、血管中の血流方向と垂直に交わる複数の垂直断層画像情報を生成する画像処理装置と、
を備えるX線CTシステムであって、
前記画像処理装置は、前記垂直断層画像情報の画像から前記血管を含む前記血管の近傍領域を抽出する抽出手段と、前記近傍領域の画像を平滑化する平滑化手段と、前記平滑化された画像から、指定範囲の画素値を有する領域を抽出する領域拡張手段と、前記血栓領域の辺縁を円弧状に確定する確定手段と、
を備えることを特徴とするX線CTシステム。
An X-ray CT apparatus that irradiates the subject with an X-ray beam, collects projection data of the subject, and obtains tomographic image information of the subject by image reconstruction of the projection data;
Transfer means for transferring the tomographic image information of the X-ray CT apparatus;
From the tomographic image information transferred by the transfer means, an image processing device that generates a plurality of vertical tomographic image information that intersects perpendicularly to the blood flow direction in the blood vessel,
An X-ray CT system comprising:
The image processing apparatus includes: an extraction unit configured to extract a region near the blood vessel including the blood vessel from an image of the vertical tomographic image information; a smoothing unit configured to smooth an image of the vicinity region; and the smoothed image. From, region expansion means for extracting a region having a pixel value of the specified range, and a determination means for determining the edge of the thrombus region in an arc,
An X-ray CT system comprising:
JP2002217726A 2002-07-26 2002-07-26 Thrombus section image extraction method, image processor and x-ray ct system Pending JP2004057340A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2002217726A JP2004057340A (en) 2002-07-26 2002-07-26 Thrombus section image extraction method, image processor and x-ray ct system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2002217726A JP2004057340A (en) 2002-07-26 2002-07-26 Thrombus section image extraction method, image processor and x-ray ct system

Publications (1)

Publication Number Publication Date
JP2004057340A true JP2004057340A (en) 2004-02-26

Family

ID=31939110

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2002217726A Pending JP2004057340A (en) 2002-07-26 2002-07-26 Thrombus section image extraction method, image processor and x-ray ct system

Country Status (1)

Country Link
JP (1) JP2004057340A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008237908A (en) * 2007-03-27 2008-10-09 General Electric Co <Ge> Detection of thrombi in ct using energy discrimination
JP2009183508A (en) * 2008-02-07 2009-08-20 Hitachi Medical Corp Diagnostic image support system
JP2010167188A (en) * 2009-01-26 2010-08-05 Toshiba Corp Medical image diagnostic apparatus, image data output device, and control program for outputting image data
JP2011255033A (en) * 2010-06-10 2011-12-22 Shimadzu Corp Radiation image region-extracting apparatus, radiation image region-extracting program, and radiographic apparatus

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000172829A (en) * 1998-12-07 2000-06-23 Hitachi Ltd Area extracting method

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000172829A (en) * 1998-12-07 2000-06-23 Hitachi Ltd Area extracting method

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008237908A (en) * 2007-03-27 2008-10-09 General Electric Co <Ge> Detection of thrombi in ct using energy discrimination
JP2009183508A (en) * 2008-02-07 2009-08-20 Hitachi Medical Corp Diagnostic image support system
JP2010167188A (en) * 2009-01-26 2010-08-05 Toshiba Corp Medical image diagnostic apparatus, image data output device, and control program for outputting image data
JP2011255033A (en) * 2010-06-10 2011-12-22 Shimadzu Corp Radiation image region-extracting apparatus, radiation image region-extracting program, and radiographic apparatus

Similar Documents

Publication Publication Date Title
JP6688557B2 (en) X-ray CT system
US9811903B2 (en) Systems and methods for identifying bone marrow in medical images
JP5588859B2 (en) X-ray CT scan simulator and X-ray CT apparatus
JP5718014B2 (en) X-ray CT system
JP2007144174A (en) Method and system for automatically determining regions in scanned object
JP2000350726A (en) Method for determining length of object and computed tomography system
WO2007114470A1 (en) X-ray ct scan simulator device, x-ray ct device and x-ray ct scan simulator program
JP2004320771A (en) Method for performing digital subtraction angiography
JP2004230154A (en) Volumetric ct system and method utilizing multiple detector panels
JP2005103263A (en) Method of operating image formation inspecting apparatus with tomographic ability, and x-ray computerized tomographic apparatus
JP4474304B2 (en) Ring artifact removal method and X-ray CT apparatus
JP2007252898A (en) Image display and x-ray ct scanner
JP2000217810A (en) Cone beam ct device
JP2007014483A (en) Medical diagnostic apparatus and diagnostic support apparatus
CN113940691A (en) System and method for patient positioning for image acquisition
JP2016168120A (en) Medical image processor and image display control method in medical image processor
JP6900218B2 (en) Image processing device and display method of radiographic image data
CN116019474B (en) Multi-source imaging device and method
KR101768520B1 (en) A method of integrated operation of chest X-ray digital radiography and chest digital tomosynthesis
JP6656883B2 (en) X-ray computed tomography apparatus and medical diagnostic system
JP2004057340A (en) Thrombus section image extraction method, image processor and x-ray ct system
EP3809376A2 (en) Systems and methods for visualizing anatomical structures
JP4305720B2 (en) X-ray CT system
KR101525040B1 (en) Method and Apparatus of Generation of reference image for determining scan range of pre-operative images
JP2004081394A (en) In-tissue fat evaluating method, image processor, and x-ray computerized tomography (ct) system

Legal Events

Date Code Title Description
A625 Written request for application examination (by other person)

Free format text: JAPANESE INTERMEDIATE CODE: A625

Effective date: 20041216

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20060411

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20070116

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20070619