JP2003325503A

JP2003325503A - 画像表示装置及びその制御方法、並びにプログラムコード、記憶媒体

Info

Publication number: JP2003325503A
Application number: JP2002318004A
Authority: JP
Inventors: Akihiko Nishide; 明彦西出; Hirofumi Yanagida; 弘文柳田; Koji Segawa; 晃司瀬川
Original assignee: GE Medical Systems Global Technology Co LLC
Current assignee: GE Medical Systems Global Technology Co LLC
Priority date: 2002-03-06
Filing date: 2002-10-31
Publication date: 2003-11-18
Anticipated expiration: 2022-10-31
Also published as: JP4130118B2

Abstract

(57)【要約】【課題】針の先が表示された画面を操作者に通知する
こと。また、針の先と関心領域との位置関係に応じて、
操作者にメッセージを通知すること。【解決手段】３次元画像に３次元ラベリング処理を施
し（Ｓ１２０３）、針の候補と推定される領域を楕円体
近似し、主軸方向や各径、各径の比や円形度、球体度な
どを求め（Ｓ１２０４）、求めた結果に応じて針ではな
いと判断された場合、針の候補から落とすためにラベル
番号を0にクリアする（Ｓ１２０６）。針の画像と見な
された場合、針の画像に対して細線化処理、端点処理を
施す（Ｓ１２０７）。次に、針の画像が表示されている
画面を特定し（Ｓ１２０８）、特定された画面のフレー
ムの色、針の先の位置（もしくはその近傍の位置）の画
素の色を変更する（Ｓ１２０９）。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】本発明は、針を穿刺された被
検体に対してフルオロスキャンを行うことで得られる複
数のＸ線断層像を表示する表示部を複数備える画像表示
装置及びその制御方法並びにプログラムコード、記憶媒
体に関するものである。

【０００２】

【従来の技術】Ｘ線ＣＴシステムの放射線断層撮像シス
テムによる実時間（リアルタイム）透視検査（fluorosc
opy）とよばれる診断においては、患者（被検体）の患
部等の興味のある部位（関心領域；ＲＯＩ、Region of
Ineterestとよばれる）に薬を注入したり組織の一部を
取り出したりする目的で針を穿刺する。操作者は、針を
穿刺した後、被検体のＸ線断層像により、針の位置を確
認している。

【０００３】上述の目的で針を用いる場合、操作者は得
られるＸ線断層像から針の位置、特に針の先の位置を確
認する必要がある。得られるＸ線断層像の表示例を図３
に示す。

【０００４】図３（ａ）において３０１はＸ線管で、３
０２は４列の検出器列である。図３（ｂ）は図３（ａ）
に示した４列の検出器列３０２を被検体の搬送方向（以
下、ｚ軸方向）にＡ、Ｂ、Ｃの３つの検出器群に分けた
ときに、夫々の検出器群で得られたＸ線断層像の表示例
を示す。検出器群Ｂには２列の検出器列が含まれてお
り、この検出器群からは２つの画像が得られるが、この
２つの画像を合成して図３（ｂ）の画面に表示する。

【０００５】図３（ｃ）は図３（ａ）に示した４列の検
出器列３０２を被検体の搬送方向（以下、ｚ軸方向）に
Ｄ、Ｅの２つの検出器群に分けたときに、夫々の検出器
群で得られたＸ線断層像の表示例を示す。尚、検出器群
Ｄ、Ｅには夫々２つの検出器列が含まれているが、上述
の検出器群Ｂと同様に、含まれている２列の検出器によ
る画像を合成してから表示する。例えば検出器群Ｄに属
する２列の検出器から得られる２つの画像を合成して図
３（ｃ）に示す画面に表示する。

【０００６】

【発明が解決しようとする課題】しかし従来では、被検
者のＸ線断層像は複数得られるので、特に実時間性（リ
アルタイム性）を要求されるＣＴフルオロスキャンで
は、図３（ｂ）や図３（ｃ）に示す表示例では、どの検
出器群による画面に針の画像が表示されているかを瞬間
的に探すのが手間であった。また、表示されている画面
が分かっても、その断層像内に針の先の位置を特定する
ことが困難であった。

【０００７】また上述の通り、針は関心領域に対して穿
刺するので、関心領域に対して正確に、且つ安全に針を
穿刺するためには、操作者は針の先と関心領域との位置
関係を常に把握する必要がある。

【０００８】本発明はかかる問題点に鑑みてなされたも
のであり、針の先が表示された画面を操作者に通知する
ことを目的とする。

【０００９】また本発明の別の目的としては、針の先と
関心領域との位置関係に応じて、操作者にメッセージを
通知することを目的とする。

【００１０】

【課題を解決するための手段】本発明の目的を達成する
ために、例えば本発明の画像表示装置は以下の構成を備
える。

【００１１】すなわち、針を穿刺された被検体に対して
フルオロスキャンを行うことで得られるＸ線断層像を表
示する表示部を複数備える画像表示装置であって、前記
針の画像を含むＸ線断層像から、前記針の端点を検出す
る検出手段と、前記複数の表示部において、前記検出手
段で検出された針の端点が表示されている表示部を特定
し、当該表示部を通知する通知手段とを備えることを特
徴とする。

【００１２】本発明の目的を達成するために、例えば本
発明の画像表示装置は以下の構成を備える。

【００１３】すなわち、針を穿刺された被検体に対して
フルオロスキャンを行うことで得られるＸ線断層像を表
示する表示部を複数備える画像表示装置であって、前記
針の画像を含むＸ線断層像であって、当該Ｘ線断層像の
集合である３次元画像を２値化して３次元ラベリング処
理を施し、３次元連続領域を抽出する３次元連続領域抽
出手段と、前記３次元ラベリング処理の結果に従って、
前記３次元画像から前記針の３次元連続領域と推定され
る領域を特定し、当該領域を包含する楕円体近似情報、
もしくは直方体近似情報に基づいて前記画像が針の画像
であるか否かを判断する判断手段と、前記判断手段で針
の画像と判断された画像から針の端点を検出する検出手
段と、前記複数の表示部において、前記検出手段で検出
された針の端点が表示されている表示部を特定し、当該
表示部を通知する通知手段とを備えることを特徴とす
る。

【００１４】

【発明の実施の形態】以下添付した図面を参照して、本
発明の画像表示装置を操作コンソールに適用した好適な
実施形態に従って、詳細に説明する。

【００１５】［第１の実施形態］図１は、本実施形態の
Ｘ線ＣＴシステムのブロック構成図である。図示のよう
に本システムは、被検体へのＸ線照射と被検体を透過し
たＸ線を検出するためのＸ線検出機構を一体的に取り付
けるガントリ装置１００と、ガントリ装置１００に対し
て各種動作設定を行うとともに、ガントリ装置１００か
ら出力されたデータに基づいてＸ線断層像を再構成し、
表示する操作コンソール２００により構成されている。

【００１６】ガントリ装置１００は、その全体の制御を
つかさどるメインコントローラ１をはじめ、以下の構成
を備える。

【００１７】２は操作コンソール２００との通信を行う
ためのインタフェース、３はテーブル１２上に横たえた
被検体（患者）を搬送（図面に垂直な方向で以下、ｚ軸
ともいう）するための空洞部を有するガントリであり、
内部には、Ｘ線発生源であって、Ｘ線管コントローラ５
により駆動制御されるＸ線管４、Ｘ線の照射範囲を制御
するためのスリットを有するコリメータ６、コリメータ
６のＸ線照射範囲を制御するスリット幅の調整用モータ
であるモータ７ａが設けられている。このモータ７ａの
駆動はコリメータコントローラ７により制御される。

【００１８】また、ガントリ３は、被検体を透過したＸ
線を検出する、複数列の検出器から成るマルチ検出器で
あるＸ線検出部８、およびＸ線検出部８で得た透過Ｘ線
より得られる投影データを収集するデータ収集部９も備
える。Ｘ線管４及びコリメータ６と、Ｘ線検出部８は互
いに空洞部分を挟んで、すなわち、被検体を挟んで対向
する位置に設けられ、その関係が維持された状態でガン
トリ３のまわりを回転運動するようになっている。この
回動運動は、モータコントローラ１１からの駆動信号に
より駆動される回転モータ１０によって行われる。ま
た、被検体を乗せるテーブル１２は、ｚ軸方向への搬送
がなされるが、その駆動はテーブルモータ１３によって
行われる。

【００１９】メインコントローラ１は、インタフェース
２を介して受信した各種コマンドの解析を行い、それに
基づいて上記のＸ線管コントローラ５、コリメータコン
トローラ７、モータコントローラ１１、テーブルモータ
コントローラ１４、そして、データ収集部９に対し、各
種制御信号を出力することになる。また、データ収集部
９で収集された投影データは、データ入力インタフェー
ス６０を介して操作コンソール２００に送出される。

【００２０】操作コンソール２００は、いわゆるワーク
ステーションであり、図示するように、装置全体の制御
をつかさどるＣＰＵ５１、ブートプログラム等を記憶し
ているＲＯＭ５２、主記憶装置として機能するＲＡＭ５
３をはじめ、以下の構成を備える。

【００２１】ＨＤＤ５４は、ハードディスク装置であっ
て、ここにＯＳ、後述するスキャン制御プログラムのほ
か、ガントリ装置１００に各種指示を与えたり、ガント
リ装置１００より受信したデータに基づいてＸ線断層像
を再構成するためのプログラムが格納されている。ま
た、後述する画像表示のためのプログラムコードも格納
されている。ＶＲＡＭ５５は表示しようとするイメージ
データを展開するメモリであり、ここにイメージデータ
等を展開することでＣＲＴ５６やＣＲＴ６１に表示させ
ることができる。なお、ＣＲＴ６１は針の穿刺を行う操
作者がその表示画面が目視可能な位置に設けられてい
る。

【００２２】５７及び５８は、各種設定を行うためのキ
ーボードおよびマウスである。また、５９はガントリ装
置１００と通信を行うためのインタフェースである。

【００２３】次に、以上の構成を備えるＸ線ＣＴシステ
ムを用いて、被検体のＸ線断層像を得、このＸ線断層像
から針の画像を特定する方法について説明する。なお、
Ｘ線断層像を得る際のスキャン方法は、実時間で断層像
を表示できるフルオロスキャンを用いる。

【００２４】図２に上述のＸ線ＣＴシステムにより得ら
れる被検体の複数のＸ線断層像から成る３次元画像を示
す。２０１はｚ軸方向に検出器の列の数だけ得られる被
検体のＸ線断層像の集合である３次元画像であって、各
画素が３次元的に配置された画像である。２０２は針の
画像（３次元連続領域）である。以下では、３次元画像
２０１から針の画像２０２を得るために、この３次元画
像２０１に対して適当なＣＴ値範囲で２値化した後、３
次元ラベリング処理を施す。

【００２５】図４に本実施形態における３次元ラベリン
グ装置の機能構成を示す。２値画像である３次元画像４
０１は読み出し部４０２において、２次元平面（２次元
平面の定義は上述の通り）毎に読み出される。セレクタ
４０４は後述のラベリング部４０５において用いられる
画素データを読み出し部４０２から選択し、ラベリング
部４０５に出力する。ラベリング部４０５はセレクタ４
０４が選択した画素データと後述の近傍マスクを用いて
ラベリングを行う。そして各２次元平面毎の３次元ラベ
リング情報を生成し、又各ラベルの接続情報（詳細は後
述）を生成し、ラベル接続情報格納部４０６に出力す
る。ラベル接続情報格納部４０６はラベリング部４０５
からのラベル接続情報を格納し、後述のリ・ラベリング
部４０７に出力する。リ・ラベリング部４０７は、ラベ
ル接続情報を用いて、ラベリングされた３次元画像を構
成する。

【００２６】以上の機能構成を有する本実施形態の３次
元ラベリング装置が行うラベリング処理について、以下
説明する。

【００２７】本実施形態におけるラベリング処理に用い
る３次元近傍マスクの構成を図６（ａ），（ｂ）に示
す。本実施形態における３次元近傍マスクは２層を成
し、図６（ｂ）に示す近傍マスクの上部（同図において
矢印の方向）に図６（ａ）に示す近傍マスクを重ねた構
成を備える。上述の２つの近傍マスクを重ねる際には、
図６（ａ）に示す６０２の画素と、図６（ｂ）に示す６
０３の画素との位置をあわせるようにして重ねる。本実
施形態では図６（ｃ）に示すように、注目画素を中心と
して上９画素、周囲８画素、下９画素の計２６画素を参
照する２６近傍の３次元近傍マスクを用いる。尚、図６
（ａ）に示す６０１の画素を以下注目画素と呼称する
（詳細は後述する）。

【００２８】図８に上述の本実施形態の３次元近傍マス
クを用いて３次元ラベリング処理を行う際の、各２次元
平面と３次元近傍マスクを示す。同図において８０１
ａ、８０１ｂ、８０１ｃは２次元平面（２値画像）、８
０２は図６（ｂ）に示した近傍マスクの上部（走査がｚ
軸座標の大きい方向に進む場合に、ｚ軸座標の小さい方
を上部と呼ぶ）、８０３は図６（ａ）に示した近傍マス
クの下部、８０３ａは注目画素（図６での画素６０１）
である。

【００２９】以下、２次元平面８０１ｂに含まれる各画
素に対してラベリング処理を行う場合について説明する
が、ここで説明するラベリング処理を３次元画像を構成
するすべての２次元平面（例えば２次元平面８０１ａ，
８０１ｂ，８０１ｃ）に対して行うことが３次元画像に
対してラベリング処理を行うことを包含していることは
明白である。つまりこの近傍マスクはｘ、ｙ、ｚ軸の順
に３次元的に走査される。本実施形態ではまずｘ軸方向
に走査され、次にｙ座標を進めてｘ軸方向に走査され、
そのｘ−ｙ平面の走査完了後に次のｚ軸座標に進めてい
る。即ち図８に示すように走査１−１，１−２，１−
３、、２−１，２−２，２−３、、という具合に走査さ
れる。

【００３０】まず、平面８０１ｂの左上隅のｘ＝０，ｙ
＝０，ｚ＝０の画素からｘ軸方向に１画素毎にスキャン
を行い、画素値が１である画素を検索する。右端までス
キャンを行ったら、ｙ座標を進めた１段下の行の左端
（ｘ＝０）から再度スキャンを行う。尚、スキャンを行
う最初の画素の位置やスキャンの方向はこれに限定され
るものではない。そして最初に見つかった（注目画素の
画素値が１である）画素のラベル番号を１とする。そし
て以後、画素値が１である画素が見つかったら（注目画
素の画素値が１であったら）、注目画素の周辺の近傍マ
スク内（８０２，８０３内）の既に走査済みの各画素
（各マスク画素と各画素は１対１で対応する）のラベル
番号を参照する。

【００３１】近傍マスク内の画素であって、注目画素以
外の画素に対してラベル番号ｉがついていた場合、注目
画素のラベル番号をｉとする。

【００３２】一方、近傍マスク内の画素であって、注目
画素以外の画素に対して複数のラベル番号ｉ、ｋ、ｍが
ついていた場合、本実施形態ではラベル番号ｉ、ｋ、ｍ
のうち最も小さい番号を注目画素のラベル番号とする。
この場合、ラベル番号ｉ、ｋ、ｍの画素（注目画素も含
む）はリ・ラベリングのために全て３次元接続している
ラベルであることを示す接続情報（ラベル接続情報）を
生成する。

【００３３】この接続情報について図１０を用いて説明
する。図１０は、２次元平面１０００内で同じ画像領域
であるにもかかわらず、３次元走査する方向の優先度に
より異なるラベル番号が付けられた場合に、それらの画
像領域を接続し、リ・ラベリング処理において同じラベ
ル番号を付ける処理を説明するための図である。図１０
において左の図より、画像領域１００１，１００２は同
じ画像領域であるにもかかわらず、異なるラベル番号
（１，３）が付けられている。本実施形態におけるラベ
リング処理では、このような状態が発生することがあ
る。上述の接続情報は領域１００１と１００２とが同じ
画像領域に含まれるものであることを示す情報であっ
て、その情報の記述方法は特には限定しない。この接続
情報を参照することで、領域１００１と領域１００２と
を接続して、リ・ラベリング処理で同じラベル番号（例
えばお互いのラベル番号のうちで最小の番号）を付ける
ことで、一つの領域１００３（図１０右の図参照）とし
て扱うことができる。

【００３４】また、近傍マスク内の画素であって、注目
画素以外の画素の値が全て０である場合、つまり、近傍
マスク内の（注目画素以外の）画素のどれもラベル番号
がついていない場合、既にそのラベリング処理の中で付
けたラベル番号で最も大きい番号に１を加算したラベル
番号を注目画素のラベル番号とする。

【００３５】以降、平面８０１ｂ内の全ての画素に対し
てスキャンを行いながら、上述の注目画素に対するラベ
リング処理を行う。その結果、平面８０１ｂに対するラ
ベリング情報を作成することができる。

【００３６】そして平面８０１ｂ内の全ての画素に対し
て３次元ラベリング処理が終わったら、次の平面８０１
ｃに移行し、平面８０１ｃに対する上述の３次元ラベリ
ング処理を行う。その場合にはセレクタ４０４は平面８
０１ｂの次の平面８０１ｃを読み出し部４０２から読み
出し、ラベリング処理部４０５に出力することで、ラベ
リング処理部４０５は平面８０１ｂ、８０１ｃを用いて
平面８０１ｃに含まれる画素のラベリング処理を行う。

【００３７】以上のラベリング処理によって、注目平面
（上述の例では平面８０１ｂ）に含まれる注目画素（上
述の例は注目画素）のラベル番号は、その上部の平面
（上述の例では平面８０１ａ）内で、近傍マスクに含ま
れる画素の値、もしくはラベル番号によって決まるの
で、全ての平面に対してラベリング処理を終えた時点
で、連続した３次元領域の”１”であった画像領域は平
面間で接続されて同じラベル番号になっていることにな
る。その結果、各平面内で独立したラベリング処理を行
い、後に平面間で接続して3次元領域として同じ領域に
含まれる画像領域の接続を行う処理方法に比べて、少な
くとも平面間の画像領域の接続を調べる処理時間の分だ
け、本実施形態におけるラベリング処理は高速に行うこ
とができる。これはリ・ラベリングも含めて２回の３次
元走査により３次元的に画像領域の接続関係を参照する
ことができることにも起因している。

【００３８】なお、最上面（ｚ＝0）の平面８０１ａを
注目平面とするときは上部の平面がいるため、上部平面
を構成する全てのデータを０として扱い、処理を進めて
いく。最下面も同様の処理を行う。

【００３９】以上の３次元ラベリング処理のフローチャ
ートを図９に示す。

【００４０】まず、ラベル番号の変数ｉを１に初期化す
る（ステップＳ９０１）。そして平面内の各画素を３次
元走査し、画素値が１である画素を検索する（ステップ
Ｓ９０２）。画素値が１である画素を最初に見つけた場
合（ステップＳ９０３）、その画素のラベル番号をｉ、
すなわち１とする（ステップＳ９０４）。そして以降
も、画素値が１である画素を検索する。注目画素の画素
値が０の場合は、注目画素のラベル番号は０とする（ス
テップＳ９１５）。

【００４１】画素値が１である画素（注目画素が１であ
る画素）を見つけた場合（ステップＳ９０５）、図７に
示した近傍マスクを用いて、近傍マスク内で、ラベル番
号がついた画素を検索する（ステップＳ９０６）。ステ
ップＳ９０６における検索の結果、複数の画素（ここで
は例えば３つとし、夫々ラベル番号がｊ、ｋ、ｌ、ｊ＜
ｋ＜ｌ）が見つかった場合、処理をステップＳ９０７に
移行し、ｊ、ｋ、ｌのうちの最小のラベル番号であるｊ
を注目画素のラベル番号とする（ステップＳ９０７）。
又、ラベル番号がｊ、ｋ、ｌの画素が３次元接続してい
ることを示すラベル接続情報を生成する（ステップＳ９
０８）。

【００４２】一方、ステップＳ９０６における検索の結
果、単数の画素（ラベル番号ｊ）が見つかった場合、処
理をステップＳ９０９に移行し、このラベル番号ｊを注
目画素のラベル番号とする（ステップＳ９０９）。ま
た、ステップＳ９０６における検索の結果、近傍マスク
内の（注目画素以外の）画素値が全て０である場合、処
理をステップＳ９１０に移行し、今まで付けたラベル番
号で最も大きい番号（ｉ）に１を加算した番号（ｉ＋
１）を注目画素のラベル番号とする（ステップＳ９１
０）。

【００４３】そして再び、３次元走査を行う（ステップ
Ｓ９１１）。そしてステップＳ９１２における判断によ
り、未だ全ての画素についてスキャンを行っていない場
合には処理をステップＳ９０５に移行し、上述の処理を
繰り返す。一方、平面内の画素全てのスキャンが終わっ
た場合には、３次元画像内の全ての平面を走査し終える
まで繰り返す（ステップＳ９１３）。又その後、走査途
中で見つかった３次元接続情報に基づき、リ・ラベリン
グ処理を行う（ステップＳ９１４）。具体的には上述の
接続情報を元に、各平面内で各画像領域の接続を行い、
接続を行った各画像領域で同じラベル番号を付ける。

【００４４】なお、上述の例では３次元画像に対して直
接２値化処理を行った後、３次元ラベリング処理を行っ
たが、まず３次元画像に対して３次元空間フィルタ（ロ
ーパスフィルタや中間値フィルタなど）を用いて３次元
画像雑音除去を行って、また必要に応じて針を強調する
３次元空間フィルタリングを行って、その後にラベリン
グ処理を行った方がより良い精度でラベリング結果を得
ることができる。

【００４５】以上のラベリング処理の結果、例えば図２
に示した３次元画像２０１において、針の画像２０２の
領域とそれ以外のノイズの領域とで異なるラベルが付け
られたことになる。しかし一般にラベリングされた領域
が針の画像であるか否かはこの時点では判断できない。
よって次に、得られた領域が針の画像であるか否かの判
断処理について説明する。

【００４６】３次元画像に含まれる針の画像は大まかに
見て、楕円球体もしくは直方体の画像であると考えられ
る。このことを利用して、上述の３次元ラベリング処理
により得られた針と推定される候補の領域を本実施形態
では楕円球体の３次元領域と見なし、その特徴パラメー
タを抽出して、この画像が針であるか否かを判断する。

【００４７】図５（ａ）に針と推定される領域と、この
領域を近似した楕円体を示す。５０１は針と推定される
領域で、５０２はこの領域５０１を近似する楕円体であ
る。この楕円を示す特徴パラメータとして、重心周りの
２次モーメントを求め、それから主軸方向と３つの径Ｌ
１，Ｌ２，Ｌ３を図５（ｂ）に示すように求める。実際
にはその長さは画素数として求められることになる。そ
してＬ１，Ｌ２，Ｌ３の夫々の値、Ｌｉ／Ｌｊ（ｉ、ｊ
＝１，２，３）の比が所定の基準を満たしているかの判
断を行う。針であれば、図５によると、Ｌ１，Ｌ３は小
さな値、Ｌ２は大きな値つまり、Ｌ１／Ｌ３≒１Ｌ１／Ｌ２≪１Ｌ３／Ｌ２≪１となるはずである。又、それ以外にもこの楕円の平均Ｃ
Ｔ値、円形度、球体度などを求めて、それらが針である
所定の基準を満たしているか否かの判断を行っても良
い。そしてこれらの判断から、針と推定される領域（図
５（ａ）では画像５０１）が針の画像であるか否かの判
断を行う。そしてこの判断により、針の画像でないと判
断された画像領域に対しては、ラベル番号を”０”にし
て針の候補から落とす。以上の処理により、３次元画像
における針の３次元領域を特定することができる。

【００４８】また、ラベリングされた３次元画像を各平
面に分解した場合、針の領域もしくは針の先端を含んで
いる平面をラベルを参照することで特定することができ
る。また、上記各平面のうちいくつか（例えば図３
（ｂ）では３つ）はＣＲＴ６１の表示画面上に表示され
ており、例えば図３（ｂ）に示したように各検出器群か
らのＸ線断層像を表示する場合、針の領域を含んでいる
平面を特定し、特定された画面のフレームの色を変えて
表示することで、針を含む画面を見つけやすくすること
ができる。なお、図３（ｂ）に示した３つの画像は１つ
のＣＲＴ６１の表示画面上に表示しても良いし、ＣＲＴ
６１が複数のＣＲＴで構成されている場合、夫々のＣＲ
Ｔの表示画面上に夫々の画像を表示しても良い。

【００４９】図７にＣＲＴ６１の表示画面上に表示され
る、針の画像が表示されている画面を見やすくした表示
例を示す。７０１，７０２，７０３は各検出器群からの
画像を表示する画面であって、図３（ｂ）に示した例で
は、それそれ検出器群Ａ，Ｂ，Ｃから得られた画像を表
示する画面である。同図において７０５が針の画像であ
るので、この画像を含む画面７０２のフレーム７０２ａ
の色を変化させることで、操作者は針の画像が画面７０
２に表示されていることが認識しやすくなる。

【００５０】また、３次元画像における針の先を得るた
めに、３次元画像における針の画像に対して、例えば後
述の３次元論理フィルタを用いて細線化処理を行った画
像を生成する。そして細線化処理を施された針の画像に
対して端点検出を行う。端点検出処理について具体的に
は、端点検出の３次元論理フィルタで検出された端点候
補の中のうちで、楕円体近似したときの主軸方向のベク
トルにほぼ合致したものを選ぶ。

【００５１】その結果、３次元画像における針の先端が
特定できるので、そこから針の先端の存在する平面を特
定することができる。そして平面内で針の先端の位置
（もしくはその近傍でも良い）に存在する画素の値を操
作者が認識しやすい色（例えば赤や黄色など）に変更す
る（図１１のポインタ１１０１を参照）。その結果、操
作者は、近傍に設けられたＣＲＴ６１の上記表示画面を
参照することで針の先端の位置を瞬間的に認識しやすく
なり、穿刺の作業に専念できる。なお、上述のフレーム
の色を変える処理は針の先を表示している画面のフレー
ムに対してのみ行っても良い。

【００５２】次に、上記の３次元論理フィルタについて
簡単に説明する。本実施形態で細線化処理に用いる３次
元論理フィルタは、３次元画像の注目画素と注目画素の
２６近傍の画素とを用いて、細線化処理の結果として”
１”を出力すべき場合のデータパターンを論理結合する
条件部と、その条件部に記述されたデータパターンに対
する細線化処理結果に相当する値”１”及び、条件部に
記述されたデータパターン以外のデータパターンに対す
る細線化処理結果結果に相当する値”０”を出力する出
力部となら構成される。

【００５３】図１２に上述の処理のフローチャートを示
す。３次元画像が入力されると（ステップＳ１２０
１）、入力された３次元画像に対して上述の３次元空間
フィルタリング処理（ローパスフィルタや中間値フィル
タなど）を施す（ステップＳ１２０２）。ステップＳ１
２０２における処理は必要に応じて省いても良い。次に
（必要に応じてフィルタ処理が施された）３次元画像に
対して３次元ラベリング処理が施される（ステップＳ１
２０３）。ステップＳ１２０３における３次元ラベリン
グ処理の詳細は前述のとおりである。

【００５４】次に、針と推定される３次元領域を近似し
た楕円を求め、この楕円の主軸方向と３つの径や夫々の
径の比や、楕円の円形度、球体度などを求める（ステッ
プＳ１２０４）。そして、求めた特徴パラメータの結果
が所定の条件を満たしているかの判断を行い（ステップ
Ｓ１２０５）、針ではないと判断された場合、処理をス
テップＳ１２０６にすすめ、この領域に対してラベル番
号を０にして、針の候補から外す（ステップＳ１２０
６）。全てのラベル番号について上記の処理を終えた後
に（ステップＳ１２０７，Ｓ１２１１）処理をステップ
Ｓ１２０８に進め、針と判断されたラベル番号の領域に
対して細線化処理を施す（ステップＳ１２０８）。又、
その後に端点処理も行う。次に、画面（図３（ｂ）では
３つ、図３（ｃ）では２つ）に表示されている画面で、
針が表示されている画面を特定する（ステップＳ１２０
９）。そして特定された画面のフレームの色を変更す
る。また、針の先の位置（もしくはその近傍の位置）の
画素の色を変更する（ステップＳ１２１０）。

【００５５】また、針の先が表示された画面のフレーム
の色を変更する場合、ステップＳ１２０９において針の
先が表示されている画面を特定し、ステップＳ１２１０
で特定された画面のフレームの色を変更する。

【００５６】又、本実施形態では針の画像を表示してい
る画面のフレーム、もしくは針の先を表示している画面
のフレームを操作者が認識しやすいようにその色を変更
していたが、これに限定されるものではなく、他にも例
えばフレームを所定の色で点滅させたり、針の先端近辺
を指すポインタを針の先端近傍に表示してもよい。

【００５７】図１３にガントリ装置１００と操作コンソ
ール２００の全体的な処理のフローチャートを示す。

【００５８】まず、操作コンソール２００においてキー
ボード５７，マウス５８を用いてスキャン計画が入力さ
れる（ステップＳ１３５１）。そしてスキャン計画を操
作コンソール２００に入力し終えると（ステップＳ１３
５２）、入力したスキャン計画をガントリ装置１００に
送信する（ステップＳ１３５３）。ガントリ装置１００
は操作コンソール２００からスキャン計画を受信すると
（ステップＳ１３０１）、受信したスキャン計画に基づ
いてスキャンの準備を行う（ステップＳ１３０２）。そ
して、ガントリ装置１００は受信したスキャン計画に基
づいてスキャンを行う（ステップＳ１３０３）。そして
スキャンした結果（投影データ）を操作コンソール２０
０に対して送信する（ステップＳ１３０４）。投影デー
タを操作コンソール２００に対して送信するタイミング
などは公知の方法に従う。

【００５９】操作コンソール２００は投影データを受信
し（ステップＳ１３５４）、受信した投影データに基づ
いて公知の技術によりＸ線断層像を再構成する（ステッ
プＳ１３５５）。そしてＸ線検出部８の各列を複数の検
出器群に分割した場合の各検出器群に対応した画面に断
層像を表示する（ステップＳ１３５６）。なお本ステッ
プでは、複数の検出器からの投影データの合成処理を行
い、その結果を表示する処理も含まれる。

【００６０】そして、ステップＳ１３５７において、針
の画像、もしくは針の先が表示されている画面を特定
し、特定した画面のフレームの色を変更する。また、本
ステップにおいて、針の位置、もしくはその近傍の画素
の色を変更する。本ステップにおける処理の詳細は図１
２に示すとおりである。

【００６１】以上のスキャン（ステップＳ１３０２）か
ら針先表示（ステップＳ１３５７）までの処理はスキャ
ンが連続する限り繰り返し行われる。

【００６２】尚、本実施形態では針の画像、もしくは針
の先が表示された画面のフレームの色変更処理と、針の
先の位置、もしくはその近傍の画素の値を変更する処理
の両方を行っているが、いずれか１つの処理のみを実行
しても良い。

【００６３】［第２の実施形態］本実施形態では、第１
の実施形態で説明した３次元画像における針先の位置の
検出方法に加えて、検出した針先の位置をチェックする
方法について説明する。

【００６４】一般に、針の先にはアーチファクト（偽
像）が発生し易い。図１４に針と推定される画像とアー
チファクトとの関係を示す。１４０１は針と推定される
画像で、１４０２はアーチファクトである。そこで、３
次元画像（２次元平面でも良いが）において、アーチフ
ァクト１４０２を針の先端を中心に極座標変換、もしく
はＨｏｕｇｈ変換する。その結果は図１５のようにな
り、ある径の範囲、もしくは径方向に加算したプロフィ
ールを求め、その標準偏差値（θ方向の頻度のばらつ
き）が大きければ針の先端であると判断する。その判断
結果によって、このアーチファクトが針先に特有のもの
であるか否かを判断することができ、同時に、このアー
チファクトから針の先端か否かを特定することができ
る。

【００６５】［第３の実施形態］また上記実施形態での
説明による、Ｘ線断層像中における針の先の位置を特定
する方法を用いれば、Ｘ線断層像中の病変部等の検査対
象（関心領域）と針の先との位置関係に応じたメッセー
ジを操作者に対して通知することが出来る。針の穿刺は
検査対象に対して行われるものであるので、例えば現在
の針の位置がこの検査対象に近い、遠い等のように、針
の先の位置が検査対象に対してどの位置にあるのかを操
作者に対して通知することは重要なことである。

【００６６】以下、このように針の先と検査対象との位
置関係に応じたメッセージを操作者に対して通知する場
合について説明する。なお本実施形態ではこのメッセー
ジを音声で通知する。図１６は本実施形態に係る本実施
形態のＸ線ＣＴシステムのブロック構成図である。図１
と同じ部分については同じ番号を付けており、その説明
を省略する。６２はスピーカで、ＣＰＵ５１からの指示
に応じて、ＲＡＭ５３に記憶されている音のデータをＤ
／Ａ変換し、音として出力する。

【００６７】まず同図のＸ線ＣＴシステムを用いて、ヘ
リカルスキャン、もしくはアキシャルスキャンを行う範
囲を決めるためのスカウトスキャンを行う。ここでヘリ
カルスキャンとはスキャン位置をｚ軸方向に移動させつ
つ、スキャンを連続的に行うスキャンを指す。またアキ
シャルスキャンとは、１つのスキャン位置におけるスキ
ャンが終了してから次のスキャン位置におけるスキャン
を行う作業を全てのスキャン位置に対して行うスキャン
を指す。またスカウトスキャンとは、Ｘ線管の位置を固
定したまま、被検体をｚ軸方向にスキャンするスキャン
を指す。

【００６８】操作者はスカウトスキャンにより得られた
２次元像に基づいて、ヘリカルスキャン、もしくはアキ
シャルスキャンを行う範囲を決め、決めた範囲内で更に
ヘリカル、もしくはアキシャルスキャンを行う。そして
ヘリカル、もしくはアキシャルスキャンを行った結果得
られる複数のＸ線断層像から検査対象を包含する領域で
あるＲＯＩをキーボード５７やマウス５８を用いて操作
者が指定する。

【００６９】図１８に検査対象に対してＲＯＩを設定し
た様子を示す。同図において１８００はＸ線断層像で、
１８０１が検査対象部で、１８０２がＲＯＩである。同
図では検査対象部１８０１を包含する矩形領域をＲＯＩ
として設定している。なおＲＯＩの形状は矩形に限定さ
れるものではない。なおＲＯＩを指定するＸ線断層像は
１つでも良いし、複数であっても良い。

【００７０】上記実施形態による説明の通り、Ｘ線断層
像における針の先の位置を求めることが出来るので、求
めたＸ線断層像における針の先の位置とＲＯＩの位置
（例えばＲＯＩの重心の位置）との距離ｄを求める。そ
して求めた距離ｄに応じたメッセージをスピーカ６２か
ら出力する。距離ｄとメッセージとの関係を示すテーブ
ルの一例を図１７に示す。

【００７１】ＣＰＵ５１は同図のテーブルを参照して、
距離ｄが０≦ｄ＜ｄ１であれば、ＨＤＤ５４に保存され
ているメッセージデータ「ｍｓ１．ｄａｔ」をＲＡＭ５
３に読み出し、このデータに基づいた音をスピーカ６２
から出力されるようにＣＰＵ５１がスピーカ６２を制御
する。距離ｄがｄ１≦ｄ＜ｄ２であれば、ＨＤＤ５４に
保存されているメッセージデータ「ｍｓ２．ｄａｔ」を
ＲＡＭ５３に読み出し、このデータに基づいた音をスピ
ーカ６２から出力されるようにＣＰＵ５１がスピーカ６
２を制御する。上記メッセージの出力は一定時間毎（例
えば５秒毎）に行われる。

【００７２】距離ｄの値が小さければ小さいほど、針の
先がＲＯＩに近い（検査対象に近い）ことを示している
ので、「ｍｓ２．ｄａｔ」は例えば「まだＲＯＩには遠
いです」という音のデータ、「ｍｓ１．ｄａｔ」は例え
ば「ＲＯＩに入りました」という音のデータ、というよ
うに、距離ｄが小さいほど、針の先がＲＯＩに近いこと
を操作者に通知するようにメッセージを設定する。な
お、上記ｄ１は、距離ｄが０≦ｄ＜ｄ１である場合に、
針の先がＲＯＩ内であるように動的に求めても良い。ま
た、上記テーブルのデータ、各メッセージの音のデータ
はＨＤＤ５４に保存されており、必要に応じてＲＡＭ５
３にロードされる。

【００７３】なお、複数のＸ線断層像に対してＲＯＩを
設定した場合には、針の先が存在するＸ線断層像におい
て上記距離ｄを求める処理、求めた距離ｄに基づいてメ
ッセージを出力する処理を行う。

【００７４】なお、操作者に対して通知するメッセージ
の内容は単に針の先がＲＯＩの位置に対して遠い、近い
のみを示すものに限定されることなく、他にも例えば、
ＲＯＩの位置に対して針の先が上にあるのか下にあるの
か、右にあるのか、左にあるのかを示すものでも良い。
また、針の先がＲＯＩの位置に対してずれている場合、
どの方向に修正すればよいのかを示すものでも良い。そ
の場合、針の先がどの平面に存在しているのかを求め、
存在している平面（即ち針の先の位置）とＲＯＩの位置
とを比較することで、その位置関係を通知することが出
来る。もちろんその位置関係に応じたメッセージの内容
をテーブルとして予めＨＤＤ５４に保存しておく必要が
ある。

【００７５】またメッセージの内容は言葉の発声に限定
されるものではなく、他にも例えば単調な音（例えば
「ピッ」という音）を繰り返し出力し、この繰り返し周
期を距離ｄに応じて変化させても良い（例えば距離ｄが
小さいほど繰り返し周期を短くする）。

【００７６】以上の処理を行うことで、操作者は上記メ
ッセージを聞くことで針の先がＲＯＩに対してどの位置
にあるのかを把握しながら針の穿刺を行うことができ、
より正確に且つ安全に関心領域に対して針の穿刺を行う
ことが出来る。

【００７７】［第４の実施形態］上記第３の実施形態に
おいて用いる針の先の位置を求める方法は第１，２の実
施形態で説明した方法に限定されるものではなく、より
単純な方法を用いても良い。すなわち、３次元画像を二
値化し、二値化した結果をラベリングする。ラベリング
処理以降は第１の実施形態で説明した方法を用いる。

【００７８】［その他の実施形態］また、本発明の目的
は、上述した実施形態の機能を実現するソフトウェアの
プログラムコードを記録した記憶媒体（または記録媒
体）を、システムあるいは装置に供給し、そのシステム
あるいは装置のコンピュータ（またはＣＰＵやＭＰＵ）
が記憶媒体に格納されたプログラムコードを読み出し実
行することによっても実現できるものである。この場
合、記憶媒体から読み出されたプログラムコード自体が
前述した実施形態の機能を実現することになり、そのプ
ログラムコードを記憶した記憶媒体は本発明を構成する
ことになる。また、コンピュータ（操作コンソール）が
読み出したプログラムコードを実行することにより、前
述した実施形態の機能が実現されるだけでなく、そのプ
ログラムコードの指示に基づき、コンピュータ上で稼働
しているオペレーティングシステム(ＯＳ)などが実際の
処理の一部または全部を行い、その処理によって前述し
た実施形態の機能が実現される場合も含まれる。

【００７９】本発明を上記記憶媒体に適用する場合、そ
の記憶媒体には、先に説明した（図９、及び／又は図１
２、及び／又は図１３に示すフローチャートの一部もし
くは全部）に対応するプログラムコードが格納されるこ
とになる。

【００８０】このようなプログラムコードを格納する記
憶媒体としては、例えばフロッピー（登録商標）ディス
ク、ハードディスク、光ディスク、光磁気ディスク、Ｃ
Ｄ−ＲＯＭ、磁気テープ、不揮発性のメモリカード、Ｒ
ＯＭ等を用いることができる。更には、ネットワーク
（例えばインターネット）という媒体を介してダウンロ
ードしても良いであろう。

【００８１】

【発明の効果】以上の説明により、本発明によって、針
の先が表示された画面を操作者に通知することができ
る。

【図面の簡単な説明】

【図１】本発明の第１の実施形態におけるＸ線ＣＴシス
テムのブロック構成図である。

【図２】図１に示したＸ線ＣＴシステムにより得られる
被検体の複数のＸ線断層像から成る３次元画像を示す図
である。

【図３】各検出器から得られるＸ線断層像の表示例を示
す図である。

【図４】３次元ラベリング処理の構成を示す図である。

【図５】針と推定される画像と、それを包含する楕円体
を示す図である。

【図６】近傍マスクの構成を示す図である。

【図７】ＣＲＴ６１の表示画面上に表示されている、針
が表示されている画面を見やすくした表示例を示す図で
ある。

【図８】平面と近傍マスクを示す図である。

【図９】３次元ラベリング処理のフローチャートであ
る。

【図１０】接続情報を説明する図である。

【図１１】針の先端の位置、もしくはその近傍の画素の
色を変更した画面を示す図である。

【図１２】針が表示されている画面のフレームの色を変
更する、または、針の先の位置（もしくはその近傍の位
置）の画素の色を変更する処理のフローチャートであ
る。

【図１３】ガントリ装置１００と操作コンソール２００
の全体的な処理のフローチャートである。

【図１４】針と推定される領域とアーチファクトとの関
係を示す図である。

【図１５】３次元画像において、アーチファクトを針の
先端を中心に極座標変換、もしくはＨｏｕｇｈ変換した
結果と、ある径の範囲、もしくは径方向に加算して求ま
るプロフィールを示す図である。

【図１６】本発明の第３の実施形態に係る本実施形態の
Ｘ線ＣＴシステムのブロック構成図である。

【図１７】距離ｄとメッセージとの関係を示すテーブル
の一例を示す図である。

【図１８】検査対象に対してＲＯＩを設定した様子を示
す図である。

─────────────────────────────────────────────────────

【手続補正書】

【提出日】平成１４年１１月７日（２００２．１１．
７）

【手続補正１】

【補正対象書類名】図面

【補正対象項目名】図６

【補正方法】変更

【補正内容】

【図６】

フロントページの続き (72)発明者西出明彦東京都日野市旭が丘４丁目７番地の127 ジーイー横河メディカルシステム株式会社内 (72)発明者柳田弘文東京都日野市旭が丘４丁目７番地の127 ジーイー横河メディカルシステム株式会社内 (72)発明者瀬川晃司東京都日野市旭が丘４丁目７番地の127 ジーイー横河メディカルシステム株式会社内Ｆターム(参考） 4C093 AA22 BA09 CA23 EB18 FD03 FD05 FD08 FD09 FF03 FF15 FF19 FF20 FF21 FF23 FF28 FF42 FG01 FG08 FG13 FG19

Claims

【特許請求の範囲】

【請求項１】針を穿刺された被検体に対してフルオロ
スキャンを行うことで得られるＸ線断層像を表示する表
示部を複数備える画像表示装置であって、前記針の画像を含むＸ線断層像から、前記針の端点を検
出する検出手段と、前記複数の表示部において、前記検出手段で検出された
針の端点が表示されている表示部を特定し、当該表示部
を通知する通知手段とを備えることを特徴とする画像表
示装置。
【請求項２】針を穿刺された被検体に対してフルオロ
スキャンを行うことで得られるＸ線断層像を表示する表
示部を複数備える画像表示装置であって、前記針の画像を含むＸ線断層像であって、当該Ｘ線断層
像の集合である３次元画像を２値化して３次元ラベリン
グ処理を施し、３次元連続領域を抽出する３次元連続領
域抽出手段と、前記３次元ラベリング処理の結果に従って、前記３次元
画像から前記針の３次元連続領域と推定される領域を特
定し、当該領域を包含する楕円体近似情報、もしくは直
方体近似情報に基づいて、各Ｘ線断層像における針の画
像を判断する判断手段と、前記判断手段で針の画像と判断された画像から針の端点
を検出する検出手段と、前記複数の表示部において、前記検出手段で検出された
針の端点が表示されている表示部を特定し、当該表示部
を通知する通知手段とを備えることを特徴とする画像表
示装置。
【請求項３】更に、前記３次元画像に対して画像ノイ
ズ除去を行う３次元空間フィルタリング処理を施す３次
元空間フィルタリング処理手段を備えることを特徴とす
る請求項２に記載の画像表示装置。
【請求項４】前記３次元空間フィルタリング処理手段
が前記フィルタリング処理に用いるフィルタとしては、
ローパスフィルタ、中間値フィルタ、最大値フィルタ、
最小値フィルタを含むことを特徴とする請求項３に記載
の画像表示装置。
【請求項５】前記３次元連続領域抽出手段は、注目平
面に含まれる注目画素に対してラベリングを行う場合、
前記注目画素を含む平面、及び当該平面に隣接する平面
において、前記注目画素近傍の画素群を前記近傍マスク
を用いて参照し、前記注目画素の値が１で、前記近傍マスク内でラベル番
号が付けられている画素の数が複数種の場合、前記複数
の画素のラベル番号のうち、最小のラベル番号を前記注
目画素のラベル番号とし、前記注目画素の値が１で、前記近傍マスク内でラベル番
号が付けられている画素の数が１種の場合、前記画素の
ラベル番号を前記注目画素のラベル番号とし、前記注目画素の値が１で、前記近傍マスク内でラベル番
号が付けられている画素が存在しない場合、既に用いた
ラベル番号で最も大きいラベル番号に１を加算した数を
前記注目画素のラベル番号とすることを特徴とする請求
項１乃至３のいずれか１項に記載の画像表示装置。
【請求項６】前記３次元連続領域抽出手段は、前記近
傍マスク内でラベル番号が付けられている画素の数が複
数の場合、前記複数の画素が接続していることを示す接
続情報を生成することを特徴とする請求項５に記載の画
像表示装置。
【請求項７】更に、前記接続情報に基づいて、前記複
数の画素を接続し、前記複数の画素のラベル番号を同じ
番号にする接続手段を備えることを特徴とする請求項６
に記載の画像表示装置。
【請求項８】前記近傍マスクは注目画素を含むＸ線断
層像、及び当該Ｘ線断層像に隣接するＸ線断層像に渡る
形状であって、前記注目画素近傍の画素群を参照可能な
３次元形状の近傍マスクであることを特徴とする請求項
５乃至７のいずれか１項に記載の画像表示装置。
【請求項９】前記判断手段は前記楕円体近似情報の長
径と短径に基づいて、針の３次元連続領域と推定される
領域が針であるか否かの判断を行うことを特徴とする請
求項２に記載の画像表示装置。
【請求項１０】前記判断手段は、前記長径、短径、短
径と長径の比、円形度、球体度のうち、少なくとも一つ
を求め、求めた結果を用いて針の３次元連続領域と推定
される領域が針の領域であるか否かの判断を行うことを
特徴とする請求項９に記載の画像表示装置。
【請求項１１】前記判断手段により、針であると判断
された３次元連続領域に対して、対応するラベル番号を
針の表示色にデータ変換することを特徴とする請求項２
に記載の画像表示装置。
【請求項１２】前記検出手段は、前記判断手段で針の
画像と判断された画像に対して、細線化処理を施した
後、針の端点を検出することを特徴とする請求項２に記
載の画像表示装置。
【請求項１３】前記通知手段は、前記検出手段で検出
された針の端点が表示されている表示部のフレームの色
を変更することを特徴とする請求項２に記載の画像表示
装置。
【請求項１４】前記通知手段は、前記検出手段で検出
された針の端点が表示されている表示部のフレームの色
を所定の色で点滅させることを特徴とする請求項２に記
載の画像表示装置。
【請求項１５】前記通知手段は更に、前記検出手段で
検出された針の端点が表示されている表示部において、
前記針の端点の位置、もしくはその近傍の画素の色を変
更することを特徴とする請求項２乃至１４のいずれか１
項に記載の画像表示装置。
【請求項１６】更に、前記３次元画像、もしくはＸ線
断層像において、前記検出手段で検出された針の端点近
傍に生じるアーチファクトの有無を用いて、当該針の端
点か否かをチェックするチェック手段を備えることを特
徴とする請求項２乃至１５のいずれか１項に記載の画像
表示装置。
【請求項１７】前記チェック手段は、前記Ｘ線断層像
の針の端点候補を中心に極座標変換、もしくはＨｏｕｇ
ｈ変換し、各針の端点候補まわりのＣＴ値の標準偏差を
求め、当該標準偏差を用いて前記針の端点の位置をチェ
ックすることを特徴とする請求項１６に記載の画像表示
装置。
【請求項１８】更に、Ｘ線断層像に設けられた注目領
域と前記検出手段により検出された針の端点との位置関
係に応じた音声メッセージを通知する音声通知手段を備
えることを特徴とする請求項１乃至１７の何れか１項に
記載の画像表示装置。
【請求項１９】針を穿刺された被検体に対してフルオ
ロスキャンを行うことで得られるＸ線断層像を表示する
表示部を複数備える画像表示装置が行う画像表示方法で
あって、前記針の画像を含むＸ線断層像から、前記針の端点を検
出する検出工程と、前記複数の表示部において、前記検出工程で検出された
針の端点が表示されている表示部を特定し、当該表示部
を通知する通知工程とを備えることを特徴とする画像表
示方法。
【請求項２０】針を穿刺された被検体に対してフルオ
ロスキャンを行うことで得られるＸ線断層像を表示する
表示部を複数備える画像表示装置が行う画像表示方法で
あって、前記針の画像を含むＸ線断層像であって、当該Ｘ線断層
像の集合である３次元画像を２値化して３次元ラベリン
グ処理を施し、３次元連続領域を抽出する３次元連続領
域抽出工程と、前記３次元ラベリング処理の結果に従って、前記３次元
画像から前記針の３次元連続領域と推定される領域を特
定し、当該領域を包含する楕円体近似情報、もしくは直
方体近似情報に基づいて、各Ｘ線断層像における針の画
像を判断する判断工程と、前記判断工程で針の画像と判断された画像から針の端点
を検出する検出工程と、前記複数の表示部において、前記検出工程で検出された
針の端点が表示されている表示部を特定し、当該表示部
を通知する通知工程とを備えることを特徴とする画像表
示方法。
【請求項２１】更に、前記３次元画像に対して画像ノ
イズ除去を行う３次元空間フィルタリング処理を施す３
次元空間フィルタリング処理工程を備えることを特徴と
する請求項２０に記載の画像表示方法。
【請求項２２】前記３次元空間フィルタリング処理工
程で前記フィルタリング処理に用いるフィルタとして
は、ローパスフィルタ、中間値フィルタ、最大値フィル
タ、最小値フィルタを含むことを特徴とする請求項２１
に記載の画像表示方法。
【請求項２３】前記３次元連続領域抽出工程では、注
目平面に含まれる注目画素に対してラベリングを行う場
合、前記注目画素を含む平面、及び当該平面に隣接する
平面において、前記注目画素近傍の画素群を前記近傍マ
スクを用いて参照し、前記注目画素の値が１で、前記近傍マスク内でラベル番
号が付けられている画素の数が複数種の場合、前記複数
の画素のラベル番号のうち、最小のラベル番号を前記注
目画素のラベル番号とし、前記注目画素の値が１で、前記近傍マスク内でラベル番
号が付けられている画素の数が１種の場合、前記画素の
ラベル番号を前記注目画素のラベル番号とし、前記注目画素の値が１で、前記近傍マスク内でラベル番
号が付けられている画素が存在しない場合、既に用いた
ラベル番号で最も大きいラベル番号に１を加算した数を
前記注目画素のラベル番号とすることを特徴とする請求
項１９乃至２１のいずれか１項に記載の画像表示方法。
【請求項２４】前記３次元連続領域抽出工程では、前
記近傍マスク内でラベル番号が付けられている画素の数
が複数の場合、前記複数の画素が接続していることを示
す接続情報を生成することを特徴とする請求項２３に記
載の画像表示方法。
【請求項２５】更に、前記接続情報に基づいて、前記
複数の画素を接続し、前記複数の画素のラベル番号を同
じ番号にする接続工程を備えることを特徴とする請求項
２４に記載の画像表示方法。
【請求項２６】前記近傍マスクは注目画素を含むＸ線
断層像、及び当該Ｘ線断層像に隣接するＸ線断層像に渡
る形状であって、前記注目画素近傍の画素群を参照可能
な３次元形状の近傍マスクであることを特徴とする請求
項２３乃至２５のいずれか１項に記載の画像表示方法。
【請求項２７】前記判断工程では前記楕円体近似情報
の長径と短径に基づいて、針の３次元連続領域と推定さ
れる領域が針であるか否かの判断を行うことを特徴とす
る請求項２０に記載の画像表示方法。
【請求項２８】前記判断工程では、前記長径、短径、
短径と長径の比、円形度、球体度のうち、少なくとも一
つを求め、求めた結果を用いて針の３次元連続領域と推
定される領域が針の領域であるか否かの判断を行うこと
を特徴とする請求項２７に記載の画像表示方法。
【請求項２９】前記判断工程で、針であると判断され
た３次元連続領域に対して、対応するラベル番号を針の
表示色にデータ変換することを特徴とする請求項２０に
記載の画像表示方法。
【請求項３０】前記検出工程では、前記判断工程で針
の画像と判断された画像に対して、細線化処理を施した
後、針の端点を検出することを特徴とする請求項２０に
記載の画像表示方法。
【請求項３１】前記通知工程では、前記検出工程で検
出された針の端点が表示されている表示部のフレームの
色を変更することを特徴とする請求項２０に記載の画像
表示方法。
【請求項３２】前記通知工程では、前記検出工程で検
出された針の端点が表示されている表示部のフレームの
色を所定の色で点滅させることを特徴とする請求項２０
に記載の画像表示方法。
【請求項３３】前記通知工程は更に、前記検出工程で
検出された針の端点が表示されている表示部において、
前記針の端点の位置、もしくはその近傍の画素の色を変
更することを特徴とする請求項２０乃至３２のいずれか
１項に記載の画像表示方法。
【請求項３４】更に、前記３次元画像、もしくはＸ線
断層像において、前記検出工程で検出された針の端点近
傍に生じるアーチファクトの有無を用いて、当該針の端
点か否かをチェックするチェック工程を備えることを特
徴とする請求項２０乃至３３のいずれか１項に記載の画
像表示方法。
【請求項３５】前記チェック工程では、前記Ｘ線断層
像の針の端点候補を中心に極座標変換、もしくはＨｏｕ
ｇｈ変換し、各針の端点候補まわりのＣＴ値の標準偏差
を求め、当該標準偏差を用いて前記針の端点の位置をチ
ェックすることを特徴とする請求項３４に記載の画像表
示方法。
【請求項３６】更に、Ｘ線断層像に設けられた注目領
域と前記検出工程で検出された針の端点との位置関係に
応じた音声メッセージを通知する音声通知工程を備える
ことを特徴とする請求項１９乃至３５の何れか１項に記
載の画像表示方法。
【請求項３７】コンピュータを請求項１乃至１８のい
ずれか１項に記載の画像表示装置として機能させること
を特徴とするプログラム。
【請求項３８】コンピュータに請求項１９乃至３６の
いずれか１項に記載の画像表示装置の制御方法を実行さ
せることを特徴とするプログラム。
【請求項３９】請求項３７又は３８に記載のプログラ
ムを格納することを特徴とするコンピュータ読み取り可
能な記憶媒体。