(57)【特許請求の範囲】
【請求項１】
３次元画像が投影される投影面上の各画素と任意の視点とを結ぶ複数の視線に沿って、該３次元画像中に複数の探査点を設定し、該探査点の各々の画素値を該探査点の近傍画素の画素値に基づく補間演算によって算出し、該探査点の画素値に基づいて前記視線毎に前記投影面上の画素の画素値を決定し、該投影面上の各画素から構成される投影画像を生成する投影画像生成手段を備えた画像処理装置であって、
前記３次元画像は、前記探査点の設定対象となる設定対象領域および該探査点の設定対象でない設定対象外領域とに区分されたものであり、
前記投影画像生成手段は、前記探査点の近傍画素が前記設定対象外領域の画素の場合、該設定対象外領域の近傍画素の画素値が寄与しないように、前記探査点の画素値を算出するものであることを特徴とする画像処理装置。
【請求項２】
前記投影画像生成手段が、前記設定対象外領域の近傍画素の画素値を、前記投影面上の画素の画素値の決定方法に応じて定められた所定の方法を用いて前記設定対象領域中の画素の画素値から算出した値に置換し、前記補間演算において、前記設定対象外領域の近傍画素の画素値は該置換後の値を用いるものであることを特徴とする請求項１記載の画像処理装置。
【請求項３】
前記投影画像生成手段が、前記設定対象外領域の近傍画素の画素値を、前記投影面上の画素の画素値の決定方法に応じた所定の固定値に置換し、前記補間演算において、前記設定対象外領域の近傍画素の画素値は該置換後の値を用いるものであることを特徴とする請求項１記載の画像処理装置。
【請求項４】
前記投影画像生成手段が、
前記設定対象外領域の近傍画素の画素値を前記設定対象領域の画素の画素値の最大値以下の値に置換するものであり、
前記投影面上の画素の画素値を、該画素に対応する前記視線上の全ての探査点の画素値のうちの最大値に決定するものであることを特徴とする請求項２または３記載の画像処理装置。
【請求項５】
前記投影画像生成手段が、前記設定対象外領域の近傍画素の画素値を、該近傍画素に対する前記探査点における、前記設定対象領域中の近傍画素の画素値のうちの最大値以下の値に置換するものであることを特徴とする請求項４記載の画像処理装置。
【請求項６】
前記３次元画像が、前記設定対象外領域の画素の画素値が、前記設定対象領域の画素の画素値よりも大きくなるように区分されたものであることを特徴とする請求項４または５記載の画像処理装置。
【請求項７】
前記投影画像生成手段が、
前記設定対象外領域の近傍画素の画素値を前記設定対象領域の画素の画素値の最小値以上の値に置換するものであり、
前記投影面上の画素の画素値を、該画素に対応する前記視線上の全ての探査点の画素値のうちの最小値に決定するものであることを特徴とする請求項２または３記載の画像処理装置。
【請求項８】
前記投影画像生成手段が、前記設定対象外領域の近傍画素の画素値を、該近傍画素に対する探査点における、前記設定対象領域中の近傍画素の画素値のうちの最小値以上の値に置換するものであることを特徴とする請求項７記載の画像処理装置。
【請求項９】
前記３次元画像が、前記設定対象外領域の画素の画素値が、前記設定対象領域の画素の画素値よりも小さくなるように区分されたものであることを特徴とする請求項７または８記載の画像処理装置。
【請求項１０】
３次元画像が投影される投影面上の各画素と任意の視点とを結ぶ複数の視線に沿って、該３次元画像中に複数の探査点を設定し、該探査点の各々の画素値を該探査点の近傍画素の画素値に基づく補間演算によって算出し、該探査点の画素値に基づいて前記視線毎に前記投影面上の画素の画素値を決定し、該投影面上の各画素から構成される投影画像を生成する画像処理方法であって、
前記３次元画像は、前記探査点の設定対象となる設定対象領域および該探査点の設定対象でない設定対象外領域とに区分されたものであり、
前記探査点の画素値の算出の際、前記探査点の近傍画素が前記設定対象外領域の画素の場合、該設定対象外領域の近傍画素の画素値が寄与しないように、前記探査点の画素値を算出するものであることを特徴とする画像処理方法。
【請求項１１】
コンピュータに、
３次元画像が投影される投影面上の各画素と任意の視点とを結ぶ複数の視線に沿って、該３次元画像中に複数の探査点を設定させ、該探査点の各々の画素値を該探査点の近傍画素の画素値に基づく補間演算によって算出させ、該探査点の画素値に基づいて前記視線毎に前記投影面上の画素の画素値を決定させ、該投影面上の各画素から構成される投影画像を生成させる画像処理プログラムであって、
前記３次元画像は、前記探査点の設定対象となる設定対象領域および該探査点の設定対象でない設定対象外領域とに区分されたものであり、
前記探査点の画素値の算出の際、前記探査点の近傍画素が前記設定対象外領域の画素の場合、該設定対象外領域の近傍画素の画素値が寄与しないように、前記探査点の画素値を算出させることを特徴とする画像処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、３次元画像処理に関し、特に詳しくは、表示対象外領域を含む３次元画像から投影画像を生成する処理に関するものである。
【背景技術】
【０００２】
従来、医療分野では、ＣＴ装置、ＭＲＩ装置、超音波診断装置（エコー）等で得られる３次元医用画像を所望の投影面に投影した投影画像の観察、診断が行われている。このような投影画像を得る方法としては、３次元医用画像中に、任意の視点と投影面上の各投影画素とを通る視線に沿って複数の探査点を設定し、設定された複数の探査点の画素値に基づいて、視線毎に投影画素の画素値を求める画像処理が行われている。このような画像処理としては、例えば、視線毎に探査点の画素値の最大値を抽出して投影するＭＩＰ（Maximum Intensity Projection；最大値投影）処理や最小値を抽出して投影するＭｉｎＩＰ（Minimum Intensity Projection；最小値投影）処理等が知られている。
【０００３】
また、投影画像の生成処理を行う前に、３次元医用画像中に表示対象外の領域を設定し、表示対象の領域に対して画像処理を行うことが知られている。例えば、表示対象でない大動脈の血液プールと表示対象の冠状動脈を含む心臓の３次元画像から所定の画像処理によって冠状動脈を含み且つ血液プールを除いた強調画像を作成し、強調画像に対してＭＩＰ処理を行うことによって、表示対象でない血液プールの影響を受けずに冠状動脈が強調された画像を作成して表示することが提案されている（例えば、特許文献１）。
【特許文献１】特開平１１−２４２７３９号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
投影画像を生成する際、視点と投影面の位置関係によっては、視線上の各探査点の３次元画像の座標系における位置が整数座標値にならないため、３次元画像の各画素の画素値をそのまま用いることができないことがありうる。このような場合には、各探査点の近傍の画素の画素値に基づく補間演算により、各探査点の画素値を算出する必要がある。
【０００５】
ここで、３次元画像が表示対象の領域と表示対象外の領域に区分されている場合、両領域の境界付近の探査点では、探査点は当然表示対象の領域内に存在するが、その探査点の近傍の画素が表示対象外の領域に属することもありうる。図４はこのようなケースの一例を示したものであり、視線Ｅj上の探査点Ｐji（ｘ，ｙ，ｚ）において、近傍８画素のうちの１つの画素Ｑ（ｘ1，ｙ2，ｚ2）が表示対象外の領域の画素となっている。
【０００６】
このような場合に、探査点Ｐjiの画素値を近傍の補間演算によって求めると、算出された値は、本来表示対象ではない近傍画素Ｑの画素値の影響を受けたものとなってしまう。例えば、画素値の閾値100を超える画素が表示対象外（非表示）であり、近傍画素Ｑの画素値が200、他の近傍画素（表示対象領域内）の画素値が100であり、探査点Ｐjiが各近傍画素の中心にある場合、探査点Ｐjiの画素値を補間演算によって求めると、100*(1/8)*7+200*(1/8)=112.5となり、表示対象の画素値の最大値である上記閾値よりも大きい、表示対象外とすべき領域の値となってしまう。ここで、ＭＩＰ処理によって投影画像を生成するのであれば、この探査点を含む視線上では、この探査点の画素値が採用されてしまい、投影画像中ではアーチファクトとして現れてしまう。
【０００７】
本発明はこのような事情に鑑みてなされたものであり、表示対象外の領域を含む３次元画像中の画素を投影した投影画像を生成する際に、表示対象外の領域に起因するアーチファクトを抑制することを可能にする装置および方法並びにプログラムを提供することを目的とするものである。
【課題を解決するための手段】
【０００８】
本発明による画像処理装置は、３次元画像が投影される投影面上の各画素と任意の視点とを結ぶ複数の視線に沿って、該３次元画像中に複数の探査点を設定し、該探査点の各々の画素値を該探査点の近傍画素の画素値に基づく補間演算によって算出し、該探査点の画素値に基づいて前記視線毎に前記投影面上の画素の画素値を決定し、該投影面上の各画素から構成される投影画像を生成する投影画像生成手段を備えた画像処理装置であって、前記３次元画像は、前記探査点の設定対象となる設定対象領域および該探査点の設定対象でない設定対象外領域とに区分されたものであり、前記投影画像生成手段は、前記探査点の近傍画素が前記設定対象外領域の画素の場合、該設定対象外領域の近傍画素の画素値が寄与しないように、前記探査点の画素値を算出するものであることを特徴とする。
【０００９】
本発明の画像処理方法は、３次元画像が投影される投影面上の各画素と任意の視点とを結ぶ複数の視線に沿って、該３次元画像中に複数の探査点を設定し、該探査点の各々の画素値を該探査点の近傍画素の画素値に基づく補間演算によって算出し、該探査点の画素値に基づいて前記視線毎に前記投影面上の画素の画素値を決定し、該投影面上の各画素から構成される投影画像を生成する画像処理方法であって、前記３次元画像は、前記探査点の設定対象となる設定対象領域および該探査点の設定対象でない設定対象外領域とに区分されたものであり、前記探査点の画素値の算出の際、前記探査点の近傍画素が前記設定対象外領域の画素の場合、該設定対象外領域の近傍画素の画素値が寄与しないように、前記探査点の画素値を算出するものであることを特徴とする。
【００１０】
さらに、本発明による画像処理プログラムは上記の方法をコンピュータに実行させるためのものである。
【００１１】
次に、本発明による画像処理装置および方法並びにプログラムの詳細について説明する。
【００１２】
「３次元画像」は、探査点の設定対象となる設定対象領域および探査点の設定対象でない設定対象外領域とに区分されたものであるが、設定対象領域と設定対象外領域のいずれかが明示的に設定されていればよい。
【００１３】
「設定対象領域」は、本発明において３次元画像から生成された投影画像を表示する際に表示対象となりうる領域である。
【００１４】
「（設定対象領域と設定対象外領域の）区分」の具体的方法としては、所定の閾値との大小関係に基づいて両領域の区分を行う方法や、公知の画像認識処理によって抽出された領域を両領域のいずれかに設定する方法等が挙げられる。
【００１５】
「任意の視点」は１つであってもよいし、複数あってもよい、すなわち、投影画像は１つの視点からの複数の視線に沿って３次元画像中の画素を投影した中心投影法による画像であってもよいし、複数の視点からの平行な視線に沿って３次元画像中の画素を投影した平行投影法による画像であってもよい。
【００１６】
「近傍画素」の範囲は、補間演算の方法に応じて適切に設定することができる。例えば、線形補間の場合には探索点が囲まれる２×２×２＝８画素が近傍画素となり、３次補間（tricubic）の場合には探索点が囲まれる４×４×４＝６４画素が近傍画素となる。
【００１７】
「該探査点の画素値に基づいて前記視線毎に前記投影面上の画素の画素値を決定」する処理の具体例としては、視線上の全ての探査点の画素値のうちの最大値に決定するＭＩＰ（Maximum Intensity Projection；最大値投影）処理や、視線上の全ての探査点の画素値のうちの最小値に決定するＭｉｎＩＰ（Minimum Intensity Projection；最小値投影）処理、視線上の全ての探査点の画素値の平均値を求める処理、レイキャスティング法によるボリュームレンダリング処理等が挙げられる。これらの処理のうち、ＭＩＰ処理やＭｉｎＩＰ処理のように、視線上の特定の探査点の画素値を投影面上の画素の画素値に決定する処理を用いた場合に本発明の効果は顕著となる。
【００１８】
「前記探査点の近傍画素が前記設定対象外領域の画素の場合、該設定対象外領域の近傍画素の画素値が寄与しないように、前記探査点の画素値を算出する」処理の具体的実現方法としては、前記設定対象外領域の近傍画素の画素値を他の値に置換し、前記補間演算において、前記設定対象外領域の近傍画素の画素値は該置換後の値を用いることや、設定対象外領域の近傍画素の画素値だけでなく、その画素自体の存在を補間演算の対象外とし、他の設定対象領域の近傍画素のみを用いて補間演算を行うことが考えられる。
【００１９】
前者の場合、置換する他の値としては、前記投影面上の画素の画素値の決定方法に応じて定められた所定の方法を用いて前記設定対象領域中の画素の画素値から算出した値、または、前記投影面上の画素の画素値の決定方法に応じた所定の固定値が考えられる。
【００２０】
ここで、置換する「所定の方法によって算出した値」の具体例としては、設定対象領域の画素の画素値の最大値、最小値、平均値、中央値、最頻値、設定対象領域の各画素の画素値のヒストグラムにおいて最も頻度の高いクラスの階級値、その設定対象外領域の近傍画素に対する探査点における、設定対象領域中の近傍画素の画素値の最大値、最小値、平均値、中央値等が挙げられる。また、置換する「所定の固定値」の具体例としては、その３次元画像の量子化等の仕様上の最大値や最小値等が挙げられる。
【００２１】
また、前記投影面上の画素の画素値を決定する具体的方法と、それに応じて前記設定対象領域中の画素の画素値から置換後の値を算出する具体的方法の組合せの具体例としては、前記投影面上の画素の画素値を、該画素に対応する前記視線上の全ての探査点の画素値のうちの最大値／最小値に決定する方法（ＭＩＰ／ＭｉｎＩＰ処理）の場合に、前記設定対象領域の画素の画素値の最大値以下／最小値以上の値を置換後の値としたり、該近傍画素に対する前記探査点における、前記設定対象領域中の近傍画素の画素値のうちの最大値以下／最小値以上の値を置換後の値としたりすることが挙げられる。さらに、これらの場合、３次元画像が、設定対象外領域の画素の画素値が、設定対象領域の画素の画素値よりも大きく／小さくなるように区分されたものである場合に、本発明の効果は顕著となる。
【発明の効果】
【００２２】
本発明によれば、３次元画像が投影される投影面上の各画素と任意の視点とを結ぶ複数の視線に沿って、３次元画像中に複数の探査点を設定し、探査点の各々の画素値をその探査点の近傍画素の画素値に基づく補間演算によって算出し、算出された画素値に基づいて視線毎に投影面上の画素の画素値を決定し、投影面上の各画素から構成される投影画像を生成する際に、探査点の近傍画素が探索点の設定対象外の領域の画素の場合、その設定対象外の領域の近傍画素の画素値が寄与しないように、探査点の画素値を算出するようにしたので、設定対象外の領域の画素の画素値が投影面上の各画素の画素値の決定に影響を及ぼすことがなくなり、設定対象外の領域に起因するアーチファクトが抑制される。
【００２３】
また、投影面上の各画素の画素値を決定する際に、ＭＩＰ処理やＭｉｎＩＰ処理のように、視線上の特定の探査点の画素値を投影面上の画素の画素値に決定する処理を用いた場合、その特定の探索点の近傍画素に設定対象外領域のものが含まれていると、その特定の探索点を通る視線に対応する投影面上の画素の画素値は、その設定対象外領域の近傍画素の影響を大きく受けることになるので、本発明によるアーチファクト抑制効果はより顕著となる。
【発明を実施するための最良の形態】
【００２４】
以下、図面を参照して、本発明の一実施の形態について説明する。
【００２５】
図１は、３次元医用画像処理システムの概要を示すハードウェア構成図である。図に示すように、このシステムでは、モダリティ１と、画像保管サーバ２と、画像処理ワークステーション３とが、ネットワーク９を経由して通信可能な状態で接続されている。
【００２６】
モダリティ１は、被検体を表す３次元医用画像Ｖを取得するものであり、具体的には、ＣＴ装置やＭＲＩ装置、超音波診断装置等である。
【００２７】
画像保管サーバ２は、モダリティ１で取得された３次元医用画像Ｖや画像処理ワークステーション３での画像処理によって生成された医用画像を画像データベースに保存・管理するコンピュータであり、大容量外部記憶装置やデータベース管理ソフトウェア（たとえば、ＯＲＤＢ(Object Relational Database)管理ソフトウェア）を備えている。
【００２８】
画像処理ワークステーション３は、読影者からの要求に応じて、モダリティ１や画像保管サーバ２から取得した３次元医用画像Ｖに対して画像処理を行い、生成された画像を表示するコンピュータであり、特に、読影者からの要求を入力するキーボードやマウス等の入力装置と、取得した３次元医用画像Ｖを格納可能な容量の主記憶装置と、生成された画像を表示するディスプレイとを備えている。
【００２９】
画像データの格納形式やネットワーク９経由での各装置間の通信は、ＤＩＣＯＭ（Digital Imaging and Communications in Medicine）等のプロトコルに基づいている。
【００３０】
図２は、画像処理ワークステーション３の投影画像生成機能に関連する部分を示すブロック図である。図に示すように、画像処理ワークステーション３は、読影者からの要求に応じて、読影対象の患者の３次元医用画像Ｖをモダリティ１や画像保管サーバ２から取得する画像取得部１０、取得した３次元医用画像Ｖ中に表示対象領域と非表示対象領域とを設定する表示設定部２０、任意の視点と投影面上の各画素とを結ぶ複数の視線に沿って、３次元医用画像Ｖの表示対象領域中に複数の探査点を設定し、その視線上の全ての探査点の画素値に基づいて、該視線によって投影される前記投影面上の画素の画素値を決定し、投影画像を生成する投影画像生成部６０、および生成されたボリュームレンダリング画像をディスプレイに表示する画像表示部７０から構成されており、投影画像生成部６０は、探査点設定部３０、画素値算出部４０、投影値決定部５０とを備える。
【００３１】
次に、上記各処理部で行われる処理の詳細を説明しつつ、この医用画像処理システム、特に画像処理ワークステーション３によって、投影画像を生成する処理の流れについて説明する。
【００３２】
まず、画像取得部１０が、読影者からの要求に応じて、読影対象の患者の３次元医用画像Ｖをモダリティ１や画像保管サーバ２から取得する。この３次元医用画像Ｖは、マルチスライス画像を画素（ボクセル）に分割し、３次元座標空間内に配列してなるものであり、各画素の位置は、被写体の左右方向をｘ軸、前後方向をｙ軸、上下方向をｚ軸とする３次元座標系で定義され、各画素の画素値は、その画素の位置の座標と関連づけられている。
【００３３】
次に、表示設定部２０が、画像取得部１０によって取得された３次元医用画像Ｖ中に表示対象領域（探索点の設定対象領域）と非表示対象領域（探索点の設定対象外領域）とを設定する。たとえば、読影者が３次元医用画像Ｖ中の部分領域を選択的に表示又は非表示されるように指定する入力に基づいて、３次元医用画像Ｖの各画素における表示・非表示の情報を有するマスクを作成する。
【００３４】
また、閾値処理により、非表示対象領域中の画素の画素値が、表示対象領域の画素の画素値よりも大きく／小さくなるように設定したり、画素値又は不透明度の変化（エッジ情報など）に基づいて異なる生体組織間の境界を認識することにより、３次元医用画像Ｖから抽出した生体組織領域を他の領域と区別して表示又は非表示されるように設定したりし、その設定により、３次元医用画像Ｖの各画素における表示・非表示の情報を有するマスクを作成することができる。
【００３５】
他にも、抽出したい対象領域が含まれた複数のサンプル画像を予め機械学習しておき、その学習した結果に基づいて３次元医用画像Ｖからその生体組織を抽出し、抽出した領域を他の領域と区別して表示または非表示されるように設定するなど、以上に例示した方法と異なる方法で３次元医用画像Ｖ中に表示対象領域と非表示対象領域とを設定することもできる。
【００３６】
次に、投影画像生成部６０が、表示対象領域と非表示対象領域とが設定された３次元医用画像Ｖを入力として、探査点設定部３０、画素値算出部４０、投影値決定部５０により、投影画像を構成する画素の画素値（出力画素値）を求める。
【００３７】
まず、探査点設定部３０が、初期設定ファイルや読影者によるキーボードやマウス等からの入力等により設定された視点、投影面Ｆ（大きさ、位置、画素数）により、たとえば、図３に示すように、その視点Ｅと投影面Ｆ上の投影画素の各々とを結ぶ複数の視線Ｅｊ（ｊ＝１、２、…）に沿って３次元医用画像Ｖ中の表示対象領域を所定の間隔でサンプリングした複数の探査点Ｐｊｉ（ｉ＝１、２、…、ｎ；ｎは、視線Ｅｊ上の探査点の数）を設定する。
【００３８】
次に、画素値算出部４０は、表示設定部２０において作成された各画素の表示・非表示の情報を有するマスクを参照して、各探査点Ｐｊｉの近傍にある８個の画素（近傍画素）の各々が非表示と設定された画素（非表示対象領域中の画素）であるか否かを判定し、非表示と設定された画素の場合には、その近傍画素の画素値を後述の置換処理により置換し、その後、各探査点Ｐｊｉにおける画素値ｇ (Ｐｊｉ)を、その探査点の近傍画素の画素値に基づく補間演算によって算出する。
【００３９】
この置換処理は、その非表示対象領域中の近傍画素の画素値が寄与しないように、探索点Ｐｊｉの画素値を算出するためのものである。後述する投影値決定部５０においてＭＩＰ／ＭｉｎＩＰ処理によって投影値を決定する場合には、その非表示対象領域中の近傍画素の画素値を、表示対象領域中の全画素の画素値の最大値／最小値に置換したり、その近傍画素に対する探査点における、表示対象領域中の近傍画素の画素値のうちの最大値／最小値に置換したりすることができる。
【００４０】
例えば、表示設定部２０において、画素値の閾値100を超える画素が表示対象外となるように設定しており、探査点Ｐｊｉの近傍画素が、図４に示すように、表示対象の領域に属する黒色で示した画素と、表示対象外の領域に属する白色で示した画素Ｑとからなる８個の画素であり、表示対象の領域に属する画素の画素値が各々80,85,85,90,90,95,95で、表示対象外の領域の属するがその画素値が200である場合、表示対象外の領域の画素Ｑの画素値を、表示対象の領域に属する他の近傍画素の画素値のうちの最大値95に置換すればよい。ここで、画素Ｑの画素値を、表示対象の領域に属する他の近傍画素の画素値の平均値に置換してもよいし、ＭＩＰ処理の場合は0、ＭｉｎＩＰ処理の場合はとりうる最大の画素値のような固定値に置換してもよい。
【００４１】
各探査点Ｐｊｉにおける画素値ｇ (Ｐｊｉ)は、図４に示すように、その探査点Ｐｊｉが含まれる格子を構成する８個の近傍画素の画素値を線形補間して求める。ここで、探査点Ｐｊｉの座標点を（ｘ、ｙ、ｚ）とすると画素値ｇ（ｘ、ｙ、ｚ）は次の式（１）で求められる。
【数１】
0004668289.tif 000002
【００４２】
具体的には、まず、図４に示す８個の各画素の座標値（（ｘｋ、ｙｌ、ｚｍ）、ｋ、ｌ、ｍ＝１、２）およびその座標での画素値ｇ（ｘｋ、ｙｌ、ｚｍ）を上記の式（１）代入することにより、ａｔに関する８元方程式を得、ａｔ（ｔ＝０〜７）が決定したら、探査点Ｐｊｉの座標値（ｘ、ｙ、ｚ）を代入し、画素値ｇ（ｘ、ｙ、ｚ）を求める。
【００４３】
なお、上記補間演算において、非表示対象領域中の近傍画素の画素値は、前述の置換処理後の画素値を用いる。
【００４４】
また、ここでは、各探査点Ｐｊｉにおける画素値ｇ (Ｐｊｉ)を、その探査点Ｐｊｉが含まれる格子を構成する８個の画素の画素値を線形補間して求める場合について説明したが、その探査点Ｐｊｉが含まれる格子を構成する６４個の画素（ｘ軸方向に４画素×ｙ軸方向に４画素×ｚ軸方向に４画素）の画素値を３次補間(tricubic)して求めるようにしてもよい。
【００４５】
次に、投影値決定部５０が、画素値算出部４０において算出された各探査点Ｐｊｉにおける画素値ｇ (Ｐｊｉ)に基づいて、投影面上の各画素の出力画素値Ｃｊを決定する。
【００４６】
例えば、各視線Ｅｊ上の全ての探査点の画素値のうちの最大値／最小値（ＭＩＰ／ＭｉｎＩＰ処理）、あるいは視線Ｅｊ上の全ての探査点の画素値の平均値を、その視線Ｅｊによって投影される投影面上の画素の出力画素値Ｃｊとして決定したり、各視線Ｅｊに沿って、各探査点Ｐｊｉにおける輝度値と不透明度の積を加算していき、加算結果をその視線Ｅｊが通る投影面上の投影画素の出力画素値Ｃｊとして決定したりする（ボリュームレンダリング法）。
【００４７】
このような処理を各視線について行い、投影面上のすべての投影画素の出力画素値を決定し、投影画像（最大値投影画像／最小値投影画像）を生成する。生成された投影画像は、画像表示部７０によって画像処理ワークステーション３のディスプレイに表示される。
【００４８】
このように本発明による画像処理の実施形態となる３次元医用画像処理システムでは、探索点設定部３０が、３次元画像が投影される投影面上の各画素と任意の視点とを結ぶ複数の視線に沿って、３次元画像の表示対象領域中に複数の探査点を設定し、画素値算出部４０が、探査点の各々の画素値をその探査点の近傍画素の画素値に基づく補間演算によって算出し、投影値決定部５０が、算出された各探索点の画素値に基づいて視線毎に投影面上の画素の画素値を決定し、投影面上の各画素から構成される投影画像を生成する際に、画素値算出部４０が、探査点の近傍画素が表示対象外（探索点の設定対象外）の領域の画素の場合、その設定対象外の領域の近傍画素の画素値が寄与しないように、探査点の画素値を算出するようにしたので、表示対象外（探索点の設定対象外）の領域の画素の画素値が投影面上の各画素の画素値の決定に影響を及ぼすことがなくなり、設定対象外の領域に起因するアーチファクトが抑制される。
【００４９】
上記の実施形態におけるシステム構成、処理フロー、モジュール構成等に対して、本発明の趣旨から逸脱しない範囲で様々な改変を行ったものも、本発明の技術的範囲に含まれる。また、上記の実施形態はあくまでも例示であり、上記のすべての説明が本発明の技術的範囲を限定的に解釈するために利用されるべきものではない。
【００５０】
例えば、上記の実施形態では、画像処理ワークステーション３で画像処理と画像表示の両方を行うようにしたが、画像処理サーバを別途設けてネットワーク９に接続し、画像処理はこの画像処理サーバに行わせるようにしてもよい。これにより、分散処理が図られ、たとえば、画像の表示を複数の端末で行う場合には、高性能の画像処理ワークステーションを複数台設置する必要がなくなり、システム全体のコストの低減に資する。
【図面の簡単な説明】
【００５１】
【図１】本発明の実施の形態となる３次元医用画像処理システムの概略構成図
【図２】図１の画像処理ワークステーションの投影画像生成機能および処理の流れを示すブロック図
【図３】投影画像生成部による探査点の設定処理を説明するための図
【図４】画素値算出部による探査点の画素値を算出する処理を説明するための図
【符号の説明】
【００５２】
１ モダリティ
２ 画像保管サーバ
３ 画像処理ワークステーション
９ ネットワーク
２０ 表示設定部
３０ 探査点設定部
４０ 画素値算出部
５０ 投影値決定部
６０ 投影画像生成部
Ｆ 投影面
Ｅｊ  視線
Ｐｊｉ 探査点
【図１】
0004668289.tif 000003
【図２】
0004668289.tif 000004
【図３】
0004668289.tif 000005
【図４】
0004668289.tif 000006