JP6214987B2 - 医用画像表示装置および医用画像表示システム - Google Patents

医用画像表示装置および医用画像表示システム Download PDF

Info

Publication number
JP6214987B2
JP6214987B2 JP2013198857A JP2013198857A JP6214987B2 JP 6214987 B2 JP6214987 B2 JP 6214987B2 JP 2013198857 A JP2013198857 A JP 2013198857A JP 2013198857 A JP2013198857 A JP 2013198857A JP 6214987 B2 JP6214987 B2 JP 6214987B2
Authority
JP
Japan
Prior art keywords
image
partial
display
generation processing
image generation
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.)
Active
Application number
JP2013198857A
Other languages
English (en)
Other versions
JP2015062584A (ja
Inventor
恵太 三森
恵太 三森
泰輔 岩村
泰輔 岩村
哲也 横田
哲也 横田
真人 渋谷
真人 渋谷
後藤 巌
巌 後藤
宏城 藤井
宏城 藤井
文徳 沢田
文徳 沢田
陽介 大久保
陽介 大久保
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.)
Canon Medical Systems Corp
Original Assignee
Toshiba Medical Systems Corp
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 Toshiba Medical Systems Corp filed Critical Toshiba Medical Systems Corp
Priority to JP2013198857A priority Critical patent/JP6214987B2/ja
Publication of JP2015062584A publication Critical patent/JP2015062584A/ja
Application granted granted Critical
Publication of JP6214987B2 publication Critical patent/JP6214987B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Magnetic Resonance Imaging Apparatus (AREA)

Description

本発明の実施形態は、医用画像表示装置および医用画像表示システムに関する。
近年、医用画像撮像機器(以下、モダリティ装置とする)は被検体の多様な情報を収集できるようになり、取得できる画像データも高解像度化している。それに伴って医用画像のデータサイズは増大し、医用画像表示装置における表示速度を低下させる原因となっている。特に、画像診断において大量の画像を読影して診断を下すユーザにとって、画像表示速度が低速なことは、読影に対するストレスとなる。早期診断および早期治療が求められる現在の医療においては、大量の画像データを円滑に処理することは必然である。このように、医用画像診断装置は画像表示速度の向上が求められるようになった(たとえば、特許文献1)。
特開2006−338293号公報
CTやMRIなどのモダリティ装置で取得される医用画像データは、1回の撮像で異なる撮像位置において撮像した多数のフレーム画像を取得できる3次元画像データが多くなってきている。たとえば、X線CT装置では高速回転ヘリカルスキャンと8列、16列、64列などの並列検出器との組合せシステムによる3次元撮像法が確立している。また、磁気共鳴イメージング(MRI:Magnetic Resonance Imaging)装置においては、傾斜磁場系、高周波磁場系、RFコイル系の性能向上に伴う高速撮影法の性能アップにより息止め下の3次元撮像法が確立している。3次元画像データは従来の2次元画像データより、大量の情報を取得することができる。このように大量の情報を取得することができるようになったことで、医用画像データを用いた画像診断の精度は向上してきている。
特に、MRIを用いた拡散メージングは神経線維の異方性の検出や、異方性を用いた神経線維の描出などの脳神経域において拡散テンソルイメージング(DTI:Diffusion Tensor Imaging)として発展している。DTIは、神経線維内の水分子などの粒子が熱によるブラウン運動により散らばっていく拡散効果を強調した拡散強調画像(DWI:Diffusion Weighted Image)をテンソル解析することにより得られる画像である。生体内の水分子は神経の長軸方向に沿って大きく拡散するため、DWIでは神経の走行方向に関する情報を取得することができる。多方向から撮像した複数のDWIを用いることで、神経の走行方向に対応する固有のベクトル成分をカラーRGB成分で表現したcolor encoded DTI(CEDTI)を得ることもできる。
さらに、DTIの有する神経線維の走行方向の情報から画像処理により、各ベクトルをトラッキングして神経を連続的な3次元像として表示する拡散テンソルトラクトグラフィ(DTT:Diffusion Tensor Tractgraphy)が提供されている。このDTTは脳の神経線維がどのような神経活動領域と連絡しているかを知る手掛かりとなる。したがって、DTTは脳の高次機能解析や手術支援システムへの応用が期待されている。
このようなDTIもしくはDTTは、脳に存在する何万本もの神経線維の3次元的な走行方向を表現した画像である。1本の神経線維に対する画像データは多数の点データから構成されており、神経線維の数は何万本にも及ぶ。神経線維の3次元的な走行方向を表すデータの集合を走行データと呼ぶこととする。このような神経線維の走行状態を正しく把握するためには、奥行情報が表現された投影画像の表示が不可欠となるが、多数の神経線維それぞれについて、奥行情報を求める必要があるため、投影画像の描画処理には時間がかかる。さらに、医師等のユーザは複雑に絡み合った神経の走行方向を画面上で確認するには、表示された奥行情報が付加された投影画像(以下、表示画像という)を拡大し、あるいは回転させる操作を行う必要がある。このような医師等ユーザが行う表示の拡大、回転といった操作に合わせて、医用画像表示装置は表示画像の描画処理を逐次的に行い、画像を表示する必要がある。
したがって、表示画像の描画処理を高速化し、表示速度を速めることで、医師等ユーザの拡大または回転といった操作に対応して円滑な画面表示を行う医用画像表示装置および医用画像表示システムが要望されている。
本実施形態に係る医用画像表示装置は、神経線維の3次元的な走行方向を表すデータの集合である走行データを取得するデータ取得部と、表示画像の拡大縮小率および回転角度の少なくともいずれか一方の操作指示を取得する画面情報取得部と、前記画面情報取得部によって取得された前記操作指示に基づいて、前記走行データの各点データから奥行情報を除外した投影画像を生成し、当該投影画像を複数の2次元描画領域に分けて描画処理を行う部分画像生成処理部と、前記部分画像生成処理部で前記描画領域毎に生成した画像を合成して、前記表示画像を生成する表示画像生成部と、を備え、前記部分画像生成処理部は、前記投影画像の点データの分布に基づいて前記複数の2次元描画領域を設定し、前記投影画像を当該設定した描画領域で分割し、前記描画領域毎に前記描画処理を行うことを特徴とする。
実施形態に係る医用画像表示装置の一例を示す概念的な構成図。 実施形態に係る医用画像表示装置の機能構成例を示す機能ブロック図。 脳神経線維の表示画像を模式的に表した図。 脳神経線維の表示画像を模式的に説明した図。 実施形態に係る医用画像表示装置の動作の一例を示すフローチャート。 実施形態に係る医用画像表示装置の走行データから単純投影画像を生成する方法および点データ数の算出方法を説明する図。 実施形態に係る医用画像表示装置の部分画像生成処理部が処理する点データ数が均等になるように2次元描画領域を分割する方法を説明する図。 実施形態に係る医用画像表示装置の部分画像生成処理部が2次元描画領域および走行データを分割する方法を説明する図。 実施形態に係る医用画像表示装置の各部分画像生成処理部が部分走行データから部分表示画像を生成する方法を説明する図。 実施形態に係る医用画像表示装置の表示画像生成部がそれぞれの部分表示画像を結合する方法を説明する図。 実施形態に係る医用画像表示装置の部分画像生成処理部の処理時間が均等になるように2次元描画領域を分割する方法を説明する図。 実施形態に係る医用画像表示システムの一例を示す概念的な構成図。
以下、本発明の実施の形態について添付図面を参照して説明する。
(1)構成
図1は、実施形態に係る医用画像表示装置の一例を示す概念的な構成図である。図1が示すように、医用画像表示装置100は、通信制御装置10、主記憶装置20、メインプロセッサ30、複数のサブプロセッサ40、表示部50、入力部60を備えている。
医用画像表示装置100は、通信制御装置10を介して電子ネットワーク経由で医用画像一元管理サーバ200、モダリティ装置300と接続している。医用画像表示装置100は医用画像一元管理サーバ200またはモダリティ装置300から医用画像データを電子ネットワーク経由で取得する。
実施形態に係る医用画像表示装置100で取得される医用画像データは、脳神経線維の走行方向の情報を有するDTIもしくはDTTといった奥行情報を有する3次元画像データである。脳神経線維の走行方向の情報を有する3次元画像データを走行データと呼ぶこととする。
通信制御装置10は、ネットワーク形態に応じた種々の通信プロトコルを実装する。この各種通信プロトコルに従って、医用画像表示装置100は医用画像一元管理サーバ200、モダリティ装置300と電子ネットワークを介して接続する。ここで、電子ネットワークとは、電気通信技術を利用した情報通信網全体を意味し、病院基幹LAN、無線/有線LANやインターネット網のほか、電話通信回線網、光ファイバー通信ネットワーク、ケーブル通信ネットワークおよび衛星通信ネットワークなどを含む。
取得した走行データをもとに、主記憶装置20に格納されたプログラムが、複数のサブプロセッサ40によって実行されることで、奥行情報が付加された表示画像が描画される。走行データをもとに生成された奥行情報が付加された表示画像を、単に表示画像と呼ぶこととする。
複数のサブプロセッサ40への描画処理の分配は、主記憶装置20に格納されたプログラムをメインプロセッサ30が実行することで処理される。このように医用画像表示装置100は複数のサブプロセッサ40が描画処理を分散して処理することにより、描画処理を高速化している。サブプロセッサ40は、内部に複数のプロセッサコアを実装したマルチコアプロセッサで構成されてもよいし、複数のプロセッサから成るマルチプロセッサで構成されてもよい。
主記憶装置20は、RAMとROMをはじめとする記憶媒体などにより構成され、磁気的もしくは光学的記憶媒体または半導体メモリなどの、メインプロセッサ30およびサブプロセッサ40により読み取り可能な記憶媒体を含んだ構成を有し、これらの記憶媒体内のプログラムおよびデータの一部または全部は電子ネットワークを介してダウンロードされるように構成してもよい。
表示部50は、たとえば液晶ディスプレイやOLED(Organic Light Emitting Diode)ディスプレイなどの一般的な表示装置により構成されるほか、メインプロセッサ30およびサブプロセッサ40の制御に従って表示画像をディスプレイ表示する。
入力部60は、たとえばキーボード、タッチパネル、テンキー、マウスなどの一般的な入力装置により構成され、ユーザの画像の移動、回転または拡大といった操作指示に対応した入力信号をメインプロセッサ30に出力する。
図2は、実施形態に係る医用画像表示装置100の機能構成例を示す機能ブロック図である。図2が示すように、医用画像表示装置100は、データ取得部31、画面情報取得部33、表示画像生成部35、部分画像生成処理部41、表示部50、入力部60を有する。上記構成のうち、データ取得部31、画面情報取得部33、および表示画像生成部35は、メインプロセッサ30が主記憶装置20に格納されたプログラムを実行することで実現する機能である。また、部分画像生成処理部41の各機能は主記憶装置20に格納されたプログラムを複数のサブプロセッサ40が実行することによって実現される機能である。
図2が示すように、医用画像表示装置100の部分画像生成処理部41は、部分画像生成処理部A41a、部分画像生成処理部B41b、部分画像生成処理部C41c、部分画像生成処理部D41dの4つを備えている。これらの部分画像生成処理部A、B、C、Dは、複数のサブプロセッサ40でそれぞれ並列処理される機能である。また、本実施例に係る医用画像表示装置100では4つの部分画像生成処理部A、B、C、Dから成る部分画像生成処理部41を示しているが、部分画像生成処理部41が備える部分画像生成処理部の数は2つ以上であれば4つに限らなくてよい。
データ取得部31は、医用画像一元管理サーバ200またはモダリティ装置300から電子ネットワークを介してDTTまたはDTIなどの走行データを取得する。
画面情報取得部33は、表示部50に表示画像を表示する際に、ユーザが入力部60より入力した画像の拡大縮小率および回転角度の少なくともいずれか1つを取得する。ここで回転角度には回転方向が含まれるものとする。この表示画像の拡大縮小率および回転角度は、表示部50に最初に画像を表示する際にユーザが設定してもよいし、初期値をあらかじめ設定しておき、初期値に対する増減分をユーザが設定するようにしてもよい。表示画像がいったん表示部50に表示されたら、ユーザは表示画像を拡大もしくは回転させることによって、観察したい神経線維が見えるように表示を操作する。この操作によって取得される拡大縮小率および回転角度を画面情報取得部33は取得する。
部分画像生成処理部41は、画面情報取得部33が取得した拡大縮小率および回転角度に従って、データ取得部31で取得した走行データから、単純投影画像を生成する。単純投影画像は、走行データを画面情報取得部33が取得した拡大縮小率および回転角度で変換し、奥行情報を排除して投影処理することで生成することができる。画面情報取得部33が取得した拡大縮小率および回転角度に従って走行データから単純投影画像を生成する方法は後述する。
さらに、部分画像生成処理部41は、単純投影画像の2次元描画領域を複数ある部分画像生成処理部41の数にあわせて分割する。2次元描画領域の分割は、それぞれの部分画像生成処理部41の処理速度が均一になるように行われる。図2に示した医用画像表示装置100の実施形態では、部分画像生成処理部A41a、部分画像生成処理部B41b、部分画像生成処理部C41c、部分画像生成処理部D41dの4つの部分画像生成処理部を有するため、2次元描画領域を4つに分割する。続いて、部分画像生成処理部41は、4つに分割された部分2次元描画領域のそれぞれに対応するように走行データを分割し、4つの部分走行データを生成する。
部分画像生成処理部A41a、部分画像生成処理部B41b、部分画像生成処理部C41c、部分画像生成処理部D41dは、分割された4つの部分走行データに対応する部分表示画像を生成する。部分画像生成処理部41での部分表示画像の生成については後述する。
部分画像生成処理部41が生成した部分表示画像は表示画像生成部35によって結合され、表示画像が生成される。
表示部50は、上述したように一般的な表示装置としての機能を持ち、生成された表示画像を表示する。
図3は、脳神経線維の表示画像を模式的に表した図である。DTTは神経線維一本一本の走行方向を3次元的に表現した3次元画像である。神経線維一本一本は直線ではなく、脳内を蛇行して走行し、らせん状に他の神経線維と複雑に絡み合っている。DTTではその複雑な神経線維網を画面上に表現する際に、ユーザの手前に存在する神経線維と奥に存在する神経線維とを奥行情報を使って区別して表示している。このようなDTT画像を拡大したり、回転させたりすることによって、ユーザは観察面を様々に変更し、網目状に張り巡らされた神経線維の走行を把握する。
脳内の神経線維の数は何万本とあり、それぞれの神経線維一本一本の空間的な配置を再現したDTT画像を表示するには、描画処理に大きな負荷が生じる。また、前述したように、神経線維は直線ではなく、前後左右にある他の神経線維と複雑に絡み合っている。そのため、ユーザはDTT画像を拡大したり回転させたりして観察することによってはじめて神経線維の態様を知ることができる。しかしながら、このようなDTT画像の拡大や回転といった操作毎に描画処理は実行されるため、この描画処理速度が障壁となり、円滑に画像を表示することが困難である。
DTIもDTTと同様に神経線維の奥行き情報を有している。この奥行き情報をもとに神経線維に着色を行うCEDTIでは、DTTの描画処理と同様に、描画処理速度が障壁となり、表示速度が低下するという問題がある。
図4は、脳神経線維の表示画像を模式的に説明した図である。図4(a)は図3に示した脳神経線維の表示画像と同じ模式図である。図4(b)は、図4(a)右上の白抜き矢印で示した破線円の中央にある立方体を拡大したものである。図4(b)には3本の神経線維が立方体内に示されている。図4(c)は、図4(b)をさらに拡大した図である。表示画像は、神経線維一本一本を点データの集合で表現している。このように図4(b)では線状に示されている3本の神経線維が、これを拡大した図4(c)では、それぞれの神経線維は点が集合した線として表現されている。
このようにDTTおよびDTIによる脳神経線維の表示画像の描画は、この点データのもつ奥行情報をもとに、神経線維一本一本の奥行きを考慮して行われる。そのため、点データ数が多いほど表示画像の描画処理は遅くなる。1本の神経線維は多数の点データから構成されており、脳内の神経線維の数は何万本とあるため、これらすべての点データについて表示画像を生成しようとすると、描画処理に多大な時間がかかる。本実施形態に係る医用画像表示装置100は、表示画像の表示において障壁となっている描画処理を、複数の部分画像生成処理部41が分散して実行し、個々の部分画像生成処理部41が実行する描画処理の量を減らすことで個々の処理速度を速め、全体として表示速度を高速化するというものである。
(2)動作
図5は、実施形態に係る医用画像表示装置100の動作の一例を示すフローチャートである。
図5のST101では、データ取得部31が医用画像一元管理サーバ200もしくは、モダリティ装置300から走行データを取得する。
図5のST103では、画面情報取得部33が表示画像を表示する際にユーザが入力部60より入力した、画像の拡大縮小率および回転角度の少なくともいずれか1つを取得する。
図5のST105では、部分画像生成処理部41が走行データから単純投影画像を生成する。
図5のST107では、部分画像生成処理部41が単純投影画像の2次元描画領域に含まれる点データの総数が算出する。この点データの総数が単純投影画像のデータ量を示している。単純投影画像生成部41における単純投影画像の生成処理と2次元描画領域の点データ数の算出方法について、図6を用いて詳細に説明する。
図6は、実施形態に係る医用画像表示装置100の走行データから単純投影画像を生成する方法および点データ数の算出方法を説明する図である。図6(a)は、X軸の上部の表示画像を白抜き矢印Aの方向から観察しようとする場合を説明している。白抜き矢印Aの方向から観察する場合、ユーザは表示画像を、X軸を中心として右回りにθ度回転させる必要がある。表示画像をX軸方向にθ度回転させると、図6(a)のX軸の下方に示されたように、表示画像をユーザの視線Bの方向から観察することができるようになる。
実施形態に係る医用画像表示装置100は、このように画像を回転、拡大縮小した後の表示画像を描画する処理を分散して行うことで、表示速度を高速化するものである。走行データに含まれる多数の点データから、奥行情報を表現した表示画像の描画処理を行うには時間がかかる。その描画処理を分散するために、画面情報取得部33で取得した拡大縮小率および回転角度で投影した、奥行情報を持たない単純投影画像を生成し、その単純投影画像の2次元描画領域が有するデータ量に基づき、描画処理の速度が均一になるように処理を分散させている。単純投影画像は、奥行情報を持たない画像であるため、表示画像の描画処理と比較してはるかに高速に生成することができる。単純投影画像は、走行データを画面情報取得部33で取得した拡大縮小率および回転角度で投影する際に、奥行情報を排除して生成される。
図6(b)および図6(c)は、単純投影画像の生成方法を模式的に説明する図である。図6(b)は5つの点データから成る走行データを示している。図6(b)の右側は、この走行データを視線Bから観察する例を示している。図6(b)の右下には、視線Bと平行にZ軸が示されており、Z軸と垂直方向に直行するY軸と、水平方向に直行するX軸とがそれぞれ示されている。図6(b)左側には、各点データの座標が示されており、各点データの座標はX軸、Y軸およびZ軸で表される。たとえば、1番の点データの座標は「(X1,Y1,Z1)」、2番の点データの座標は、「(X2,Y2,Z2)」と示されている。
図6(b)の右側に示した走行データから表示画像を描画するにあたり、その描画処理を分割するために図6(c)右側に示す単純投影画像を生成する。単純投影画像は前述したとおり、走行データから奥行情報を取り除いて視線Bの方向に投影した画像である。図6(b)左側に示した各点データの座標において、奥行情報は視線Bと平行なZ軸の情報である。したがって、図6(c)左側に示すように、単純投影画像は、(b)左側に示した各点データの座標からZ座標を取り除いたデータとなる。すなわち、1番の点データの座標は「(X1,Y1)」、2番の点データの座標は、「(X2,Y2)」となる。このように生成された単純投影画像は、描画処理の障壁となっている奥行情報を持たない画像であるため、生成が容易である。
図6(b)右側に示した走行データを、ユーザが視線Bの方向から観察した場合、2番と3番の点データは4番の点データと重なり合う。したがって、図6(c)右側に示した単純投影画像では、1番、4番および5番の点データしか表示されない。しかし、単純投影画像は、図6(b)左側に示した各点データの座標からZ軸の情報を取り除いた図6(c)左側に示すデータから生成されている。したがって、走行データにおいて存在した2番と3番の点データの情報が消去されたわけではない。図6(c)左側に示すように、単純投影画像の2次元描画領域には、図6(b)左側に示したもとの走行データに存在する点データと同じ数の点データが存在する。したがって、図6(c)右側に示す単純投影画像では完全に重なって表示されない点データの場合でも、それぞれの点を1個として数えることができる。したがって、図6(c)の単純投影画像の2次元描画領域に存在する総点データ数は5個と算出することができる。
図5のST109では、部分画像生成処理部41が単純投影画像の2次元描画領域を分割するにあたり、各部分画像生成処理部41に負荷の偏りや、性能の差があるかどうか判定する。各部分画像生成処理部41に負荷の偏りや、性能の差がない場合は、図5のST111で総点データ数を部分画像生成処理部41の数で割った値が算出される。
図5のST113では、単純投影画像の2次元描画領域が部分画像生成処理部41の数に応じて分割される。それぞれの2次元描画領域が有する点データ数が均一になるように、ST111で算出した点データ数を有する部分2次元描画領域に分割する。
さらに、図5のST115では、走行データを部分画像生成処理部41の数に応じて部分走行データに分割する。走行データが分割されることで、部分2次元描画領域のそれぞれ対応する部分走行データが生成される。
なお、ST109の処理は省略可能であり、各部分画像生成処理部41に負荷の偏りや、性能の差があるかどうか判定せずにST107からST111の処理に進んでもよい。
図7は、実施形態に係る医用画像表示装置100の部分画像生成処理部41が処理する点データ数が均等になるように2次元描画領域を分割する方法を説明する図である。図2で示されたように、4つの、部分画像生成処理部A41a、部分画像生成処理部B41b、部分画像生成処理部C41c、部分画像生成処理部D41dで処理される例を説明する。
図7(a)に示した走行データを視線Bの方向から観察した表示画像を生成しようとするとき、図7(b)左側に示した単純投影画像が取得される。図7(b)右側に示した表は、単純投影画像の2次元描画領域に含まれる総点データ数、医用画像表示装置100の部分画像生成処理部41の数、および部分2次元描画領域あたりの点データ数を示す。図7(b)右表にあるように、たとえば、図7(b)左側の単純投影画像の2次元描画領域には80個の点データが含まれるものとする。本実施形態に係る医用画像表示装置100の部分画像生成処理部41の数は4つであるから、1つの部分画像生成処理部41が処理を担当する点データの数は、総点データ数を部分画像生成処理部41の数で割った数である。したがって、部分画像生成処理部41はそれぞれの部分2次元描画領域に20個の点データが含まれるように2次元描画領域を分割する。
図8は、実施形態に係る医用画像表示装置100の部分画像生成処理部41が2次元描画領域および走行データを分割する方法を説明する図である。図8(a)は2次元描画領域の点データ数が均等になるように分割する図を示している。図8(a)左側に示すように、2次元描画領域を破線で示した分割線で分割する。分割線で分けられたそれぞれA、B、C、Dの領域には、点データが20個ずつ含まれる。図8(a)右側は、2次元描画領域をA、B、C、Dの各領域を実際に分割し、部分2次元描画領域A、部分2次元描画領域B、部分2次元描画領域C、部分2次元描画領域Dに分割した図を示している。
図8(b)は走行データをそれぞれの部分2次元描画領域A乃至Dに対応する部分走行データに分割する図を示している。図8(b)左側に示すように、走行データも2次元描画領域と同じ分割線で分割する。したがって、分割線で分けられたそれぞれA、B、C、Dの領域には、点データが20個ずつ含まれる。それぞれの点データは2次元描画領域の点データと対応する点データである。図8(b)右側は、走行データのA、B、C、Dの各領域を実際に分割し、部分走行データA、部分走行データB、部分走行データC、部分走行データDに分割した図を示している。
図5のST117では、各部分画像生成処理部が部分表示画像を投影する。
図9は、実施形態に係る医用画像表示装置100の各部分画像生成処理部が部分走行データから部分表示画像を生成する方法を説明する図である。図9(a)は部分画像生成処理部A41a、図9(b)は部分画像生成処理部B41b、図9(c)は部分画像生成処理部C41c、図9(d)は部分画像生成処理部D41dを示している。各部分画像生成処理部41を示す枠内の左側には、部分画像生成処理部41で生成された部分走行データA、部分走行データB、部分走行データC、部分走行データDがそれぞれ示されている。各部分画像生成処理部41の右側には、それぞれが担当する部分走行データを視線Bの方向から投影した場合に生成される部分表示画像が示されている。各部分画像生成処理部41は、担当する部分走行データの各点データの奥行情報を計算し、部分表示画像を描画する。
図5のST119では、表示画像生成部35が各部分画像生成処理部41によって生成された部分表示画像を結合して1つの表示画像を生成する。
図5のST121では、表示画像生成部35が生成した表示画像を表示部50が表示する。
図10は、実施形態に係る医用画像表示装置100の表示画像生成部35がそれぞれの部分表示画像を結合し、結合した表示画像を表示部50が表示する方法を説明する図である。図10の表示画像生成部35の左側には、各部分画像生成処理部41が生成した部分表示画像A、部分表示画像B、部分表示画像C、部分表示画像Dが示されている。図10の表示画像生成部35の右側は、これらの部分表示画像A、B、C、Dを結合して生成された一つの表示画像が示されている。図10右側には、表示部50が示されており、表示画像生成部35が部分表示画像を結合した表示画像を表示している。
このように描画処理を分散させることによって、医用画像表示装置100の描画処理にかかる負荷が軽減され、表示画像の表示に必要な時間を短縮することができる。図7が示すように、部分2次元描画領域に含まれる点データの数が均一になるように2次元描画領域を振り分けると、各部分画像生成処理部が行う処理が均一になるため、それぞれの応答時間が均一化される。したがって、処理の終わっていない部分画像生成処理部の応答を待つ時間が発生しにくくなり、描画処理全体の処理時間を短縮し、より円滑に表示画像を表示することができる。
図5のST131では、各部分画像生成処理部41に負荷の偏りや、性能の差がある場合に、各部分画像生成処理部の処理速度が均等になる点データ数が算出される。
図11は実施形態に係る医用画像表示装置100の部分画像生成処理部41の処理時間が均等になるように2次元描画領域を分割する方法を説明する図である。図11(a)左表は、図7(b)と同様の単純投影画像の2次元描画領域について算出された総点データ数と、実施形態に係る医用画像表示装置100の部分画像生成処理部41の数が示されている。図11(a)右表は、各部分画像生成処理部41の処理時間が均一になるように部分2次元描画領域あたりの点データ数を算出した表である。この表における処理速度は、部分画像生成処理部A41の性能を1.0として算出している。処理速度を求める基準は、一定の処理に必要な時間、クロック数、コア数、スレッド数、メモリ量、データI/O量などのプロセッサの性能を示す一般的な数値に基づいて決定することができる。また、ある部分画像生成処理部41の処理速度を基準する方法以外にも、基準となる値を設定し、その基準値をもとに処理速度を算出してもよい。さらに、表示画像の拡大もしくは回転といった操作が逐次行われている場合は、部分画像生成処理部41での処理が完了しないうちに、新たな描画処理命令が到達する場合も考えられる。このような場合には、各部分画像生成処理部41の負荷に軽重が生じてしまうことも考えられる。そこで、負荷の軽重によって生ずる処理速度を動的に算出することにより、点データ数を分配してもよい。
図11(a)右表の各部分画像生成処理部41の処理速度の列から、部分画像生成処理部B41bは同じ時間で処理できる量が部分画像生成処理部A41の半分である。図11(a)右表に示されるように、処理速度はそれぞれ、部分画像生成処理部A41が1.0、部分画像生成処理部B41bが0.5、部分画像生成処理部C43が1.5、部分画像生成処理部D43が2.0となる。総点データ数80を全体の処理速度5.0で割った値に、各部分画像生成処理部41の処理速度を掛けると、同じ処理時間で各部分画像生成処理部41が処理できる点データ数を算出することができる。このように、図11(a)右表の2次元描画領域あたりの点データ数を算出することができる。
図5のST131の処理により算出された点データ数に基づき、前述した図5のST115において2次元描画領域が分割される。図11(b)左図は、図11(a)右表に示した2次元描画領域あたりの点データ数にもとづいて分割した図を示している。図8(a)左図と同様に、2次元描画領域を破線で示した分割線で分割する。図11(b)の例では、分割線で分けられたそれぞれA、B、C、Dの領域には、点データがそれぞれ16、8、24、32個含まれる。図11(b)右図は、2次元描画領域のA、B、C、Dの各領域を、部分2次元描画領域A、部分2次元描画領域B、部分2次元描画領域C、部分2次元描画領域Dを分割した図を示している。図11(c)は、図8(b)と同様に走行データのA、B、C、Dの各領域を、部分走行データA、部分走行データB、部分走行データC、部分走行データDに分割した図を示している。
図11で分割された各部分2次元描画領域および部分走行データは、それぞれの部分画像生成処理部A、B、C、Dに振り分けられ、部分表示画像が描画される。さらに、前述の図5のST119の処理により、各部分画像生成処理部41が生成した部分表示画像を合成し、表示画像を生成する。図5のST121で示すように、表示部50が表示画像をディスプレイに表示する。
図11で示したように、各部分画像生成処理部41の処理速度に合わせて2次元描画領域を分割することで、医用画像表示装置100を構成する複数の部分画像生成処理部41を最大限に活用することができる。各部分画像生成処理部41の処理速度を加味して描画処理が分配されれば、個々の部分画像生成処理部41の処理時間を均一にすることができ、処理速度の遅い部分画像生成処理部41の応答を待つ時間が発生しにくくなる。また、描画命令が発生する都度、各部分画像生成処理部41の処理速度を把握し、動的に描画処理を分配することが可能となるため、より効率的に部分画像生成処理部41を稼働させ、表示速度を高速化することができる。
(3)その他の実施形態
上述した実施形態にかかる医用画像表示装置100において、マルチプロセッサまたはマルチコアプロセッサによる描画処理の分散処理について説明したが、実施形態はこれに限られるものではない。たとえば、表示装置を備えたクライアント端末により、電子ネットワークを介して制御された複数のサーバによって、描画処理を分散して行うクライアントサーバシステムで構成されてもよい。このように複数のコンピュータによって描画処理が分散処理される医用画像表示システム1の実施形態について、以下で説明する。
図12は実施形態に係る医用画像表示システムの一例を示す概念的な構成図である。図1に示すように、医用画像表示システム1は医用画像一元管理サーバ200、医用画像撮像システム500、医用画像診断装置600、複数のコンピュータ700を備えている。医用画像撮像システム500は医用画像を実際に撮像するモダリティ装置300と、そのモダリティ装置を制御するモダリティ制御装置400から構成されている。
モダリティ装置300は、医用画像を撮像する装置であり、X線CT装置、MRIなどの各種医用画像撮像装置を含む。
モダリティ制御装置400は、モダリティ装置300で医用画像を撮像する際の撮影条件や造影処理の入力や、取得した医用画像を出力することができるモダリティ装置300のコントロール用コンピュータである。モダリティ制御装置400は、モダリティ装置300に組み込まれた構成であってもよいし、モダリティ装置300とネットワークを介して接続された独立した構成であって、一般的なパーソナルコンピュータまたは、ワークステーションのようなパーソナルコンピュータより高性能なコンピュータで構成されていてもよい。
医用画像診断装置600は、医用画像の診断に使用される端末であって、医用画像診断に使用されるソフトウェアなどが搭載されたコンピュータである。医用画像診断装置600は、一般的なパーソナルコンピュータで構成されていてもよいし、ワークステーションで構成されていてもよい。
モダリティ制御装置400および医用画像診断装置600は、その内部構成として主記憶装置20、プロセッサ30、表示手段50、入力手段60などを有している。医用画像表示システム1で描画処理された表示画像はモダリティ制御装置400および医用画像診断装置600の表示手段50を用いて表示される。画像の拡大縮小率および回転角度といった画像の操作は、入力手段60によって取得される。
医用画像一元管理サーバ200は、医用画像撮像システム500で取得した医用画像データを蓄積したデータベースである。医用画像撮像システム500で取得された医用画像データは電子ネットワーク経由で医用画像一元管理サーバ200に送信される。また、医用画像診断装置600も電子ネットワークを介して医用画像一元管理サーバ200と接続しており、医用画像診断装置600は電子ネットワーク経由で医用画像一元管理サーバ200から医用画像データを取得する。
複数のコンピュータ700は電子ネットワークを介してお互いに接続しており、さらに電子ネットワーク経由で医用画像一元管理サーバ200、医用画像撮像システム500、医用画像診断装置600から医用画像データや描画処理命令などを取得する。図12に示した医用画像表示システム1では、4つのコンピュータ700が示されているが、コンピュータ700の数は4つに限らない。
本実施形態に係る医用画像表示システム1は、モダリティ制御装置400または医用画像診断装置600が有する表示手段50に表示画像を表示する際の描画処理を、複数のコンピュータ700に分散処理させることで、表示画像の表示速度を高速化している。
このような分散制御を行う医用画像表示システム1では、モダリティ制御装置400または医用画像診断装置600がクライアント、複数のコンピュータ700がサーバとして機能するクライアントサーバシステムの形式をとる。このクライアントサーバシステムにおける分散処理の制御は、クライアント側がどのサーバにどの描画処理を実行させるかを制御する場合と、クライアント側では分散処理命令の生成しか行わず、すべての処理をサーバ側で行う場合とがある。
クライアント側がどのサーバにどの描画処理を実行させるかを制御するクライアントサーバシステムの場合、サーバ側のコンピュータ700は部分表示画像の描画処理のみを行う。したがって、クライアントであるモダリティ制御装置400または医用画像診断装置600は、単純投影画像の生成、2次元描画領域および走行データの分割を行う。クライアントは、サーバである複数のコンピュータ700に電子ネットワーク経由で、分割された部分2次元描画領域と部分走行データと描画処理命令とを送信し、各コンピュータ700は部分表示画像を生成し、クライアントに返却する。クライアントであるモダリティ制御装置400または医用画像診断装置600は、コンピュータ700から取得した、部分表示画像を結合し、表示画像を表示手段50に表示する。
一方、クライアント側が分散処理命令の生成しか行わず、すべての処理をサーバ側で行うクライアントサーバシステムの場合、クライアントであるモダリティ制御装置400または医用画像診断装置600は分散処理命令を生成し、サーバであるコンピュータ700に送信する。モダリティ制御装置400または医用画像診断装置600の主記憶装置20に格納されたプログラムがプロセッサ30によって実行されることで分散処理命令が生成される。生成された分散処理命令は、ネットワークを介してサーバである複数のコンピュータ700に送信される。複数のコンピュータ700が連携して、単純投影画像の生成、2次元描画領域および走行データの分割、部分表示画像の生成を行う。サーバである複数のコンピュータ700によって生成された部分表示画像は、電子ネットワーク経由でクライアントであるモダリティ制御装置400または医用画像診断装置600に送信された後、結合されて表示画像として表示手段50に表示される。
前述の単純投影画像の生成、2次元描画領域の分割および走行データの分割、部分表示画像の生成、部分表示画像の結合といった各処理は、医用画像表示装置100で説明した方法と同様の方法で行われる。
表示画像は何万本もの神経線維それぞれについて奥行情報を加味した画像である。神経線維一本一本は、直線ではなく、蛇行していたり、らせん状に他の神経線維と絡み合っているため、奥行情報を計算しながら描画する処理には多大な負荷がかかる。表示画像を表示する医用画像表示装置100は、このような画像を拡大し、回転させる操作をするたびに画像を表示する必要がある。医用画像表示システム1では、特にこの描画処理をサーバである複数のコンピュータ700に分散させて処理させることで、描画処理時間を短縮し、表示速度を高速化するものである。
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
1 医用画像表示システム
10 通信制御装置
20 主記憶装置
30 メインプロセッサ、プロセッサ
31 データ取得部
33 画面情報取得部
35 表示画像生成部
40 サブプロセッサ
41 部分画像生成処理部
50 表示部、表示手段
60 入力部、入力手段
100 医用画像表示装置
200 医用画像一元管理サーバ
300 モダリティ装置
400 モダリティ制御装置
500 医用画像撮像システム
600 医用画像診断装置
700 コンピュータ

Claims (9)

  1. 神経線維の3次元的な走行方向を表すデータの集合である走行データを取得するデータ取得部と、
    表示画像の拡大縮小率および回転角度の少なくともいずれか一方の操作指示を取得する画面情報取得部と、
    前記画面情報取得部によって取得された前記操作指示に基づいて、前記走行データの各点データから奥行情報を除外した投影画像を生成し、当該投影画像を複数の2次元描画領域に分けて描画処理を行う部分画像生成処理部と、
    前記部分画像生成処理部で前記描画領域毎に生成した画像を合成して、前記表示画像を生成する表示画像生成部と、を備え、
    前記部分画像生成処理部は、前記投影画像の点データの分布に基づいて前記複数の2次元描画領域を設定し、前記投影画像を当該設定した描画領域で分割し、前記描画領域毎に前記描画処理を行うこと、
    を特徴とする医用画像表示装置。
  2. 前記部分画像生成処理部は、前記複数の2次元描画領域それぞれの描画処理の処理時間が均一になるように前記複数の2次元描画領域を設定すること、
    を特徴とする請求項1に記載の医用画像表示装置。
  3. 前記部分画像生成処理部は、前記複数の2次元描画領域それぞれに含まれる点データの数が均一になるように前記複数の2次元描画領域を設定すること、
    を特徴とする請求項1に記載の医用画像表示装置。
  4. 前記部分画像生成処理部は、描画処理の処理速度が異なる複数の部分画像生成処理部を備え、
    当該複数の部分画像生成処理部の処理速度と前記投影画像の点データの分布とに応じて、前記複数の部分画像生成処理部でそれぞれ行われる描画処理の処理時間が均一になるよう前記複数の2次元描画領域を設定すること、
    を特徴とする請求項2に記載の医用画像表示装置。
  5. 前記走行データはDTI(Diffusion Tensor Imaging)およびDTT(Diffusion Tensor Tractgraphy)のいずれか一方の画像データである、
    ことを特徴とする請求項1に記載の医用画像表示装置。
  6. ネットワークを介して接続した複数のコンピュータによって分散処理を行う医用画像表示システムであって、
    神経線維の3次元的な走行方向を表すデータの集合である走行データを取得するデータ取得手段と、
    表示画像の拡大縮小率および回転角度の少なくともいずれか一方の操作指示を取得する画面情報取得手段と、
    前記画面情報取得手段によって取得された前記操作指示に基づいて、前記走行データの各点データから奥行情報を除外した投影画像を生成し、当該投影画像を複数の2次元描画領域に分けて描画処理を行う部分画像生成処理手段と、
    前記部分画像生成処理手段前記描画領域毎に生成した画像を合成して、前記表示画像を生成する表示画像生成手段と、を備え、
    前記部分画像生成処理手段は、前記投影画像の点データの分布に基づいて前記複数の2次元描画領域を設定し、前記投影画像を当該設定した描画領域で分割し、前記描画領域毎に前記描画処理を行うこと、
    を特徴とする医用画像表示システム。
  7. 前記部分画像生成処理手段は、前記複数の2次元描画領域のそれぞれの描画処理の処理時間が均一になるように前記複数の2次元描画領域を設定すること、
    を特徴とする請求項6に記載の医用画像表示システム。
  8. 前記部分画像生成処理手段は、前記複数の2次元描画領域それぞれに含まれる点データの数が均一になるように前記複数の2次元描画領域を設定すること、
    を特徴とする請求項6に記載の医用画像表示システム。
  9. 前記部分画像生成処理手段は、描画処理の処理速度が異なる複数の部分画像生成処理手段を備え、当該複数の部分画像生成処理手段の処理速度と前記投影画像の点データの分布とに応じて、前記複数の部分画像生成処理手段でそれぞれ行われる描画処理の処理時間が均一になるように前記複数の2次元描画領域を設定すること、
    を特徴とする請求項7に記載の医用画像表示システム。
JP2013198857A 2013-09-25 2013-09-25 医用画像表示装置および医用画像表示システム Active JP6214987B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2013198857A JP6214987B2 (ja) 2013-09-25 2013-09-25 医用画像表示装置および医用画像表示システム

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2013198857A JP6214987B2 (ja) 2013-09-25 2013-09-25 医用画像表示装置および医用画像表示システム

Publications (2)

Publication Number Publication Date
JP2015062584A JP2015062584A (ja) 2015-04-09
JP6214987B2 true JP6214987B2 (ja) 2017-10-18

Family

ID=52831056

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2013198857A Active JP6214987B2 (ja) 2013-09-25 2013-09-25 医用画像表示装置および医用画像表示システム

Country Status (1)

Country Link
JP (1) JP6214987B2 (ja)

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05324583A (ja) * 1992-05-26 1993-12-07 Dainippon Screen Mfg Co Ltd 画像データ処理装置
US6642716B1 (en) * 2002-05-15 2003-11-04 Koninklijke Philips Electronics, N.V. Diffusion tensor magnetic resonance imaging including fiber rendering using hyperstreamlines
JP2006338293A (ja) * 2005-06-01 2006-12-14 Sony Corp 画像データ割当て方法、画像処理制御装置および画像処理システム
WO2010122916A1 (ja) * 2009-04-22 2010-10-28 株式会社 日立メディコ 磁気共鳴イメージング装置及び繊維状組織の走行方向表示方法

Also Published As

Publication number Publication date
JP2015062584A (ja) 2015-04-09

Similar Documents

Publication Publication Date Title
US20110227934A1 (en) Architecture for Volume Rendering
US9808213B2 (en) Image processing apparatus, image processing method, medical image diagnostic system, and storage medium
JP5992819B2 (ja) 画像処理装置、医用画像診断装置、画像処理方法、及び画像処理プログラム
JP5661382B2 (ja) 画像表示装置
Jang et al. Multi GPU implementation of iterative tomographic reconstruction algorithms
JP5822554B2 (ja) 画像処理装置、画像処理方法、撮影システム及びプログラム
Dehal et al. Gpu computing revolution: Cuda
CN102835974A (zh) 基于并行计算机的医学超声三维成像方法
JP7083015B2 (ja) 超解像処理装置、方法及びプログラム
US7647593B2 (en) Image processing system for volume rendering
JP2024010085A (ja) 情報処理装置、情報処理方法およびプログラム
JP4668289B2 (ja) 画像処理装置および方法並びにプログラム
JP6491471B2 (ja) 画像処理装置、画像処理方法およびプログラム
JP6214987B2 (ja) 医用画像表示装置および医用画像表示システム
Ciżnicki et al. Efficient isosurface extraction using marching tetrahedra and histogram pyramids on multiple GPUs
JP6399554B2 (ja) Mr画像の超解像化システム
Zhang et al. MR and ultrasound cardiac image dynamic visualization and synchronization over Internet for distributed heart function diagnosis
JP2008067915A (ja) 医用画像表示装置
JP2020062382A (ja) 医用画像表示装置、医用画像表示方法、およびプログラム
Xu et al. Probabilistic brain fiber tractography on GPUs
JP2006000126A (ja) 画像処理方法および装置並びにプログラム
Bagarinao et al. Enabling on-demand real-time functional MRI analysis using grid technology
JP2020191061A (ja) 画像処理装置及び画像生成方法
Wu et al. Toward a multimodal diagnostic exploratory visualization of focal cortical dysplasia
Yu et al. 3D Reconstruction of Medical Image Based on Improved Ray Casting Algorithm

Legal Events

Date Code Title Description
A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A711

Effective date: 20160511

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20160610

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20170530

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20170531

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20170728

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: 20170822

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20170920

R150 Certificate of patent or registration of utility model

Ref document number: 6214987

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