特許 6491823 医用画像診断支援装置、方法およびプログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6491823

(24)【登録日】2019年3月8日

(45)【発行日】2019年3月27日

(54)【発明の名称】医用画像診断支援装置、方法およびプログラム

(51)【国際特許分類】

   A61B 6/03 20060101AFI20190318BHJP

   A61B 5/055 20060101ALI20190318BHJP

【ＦＩ】

   A61B6/03 360G

   A61B6/03 360B

   A61B5/055 380

【請求項の数】8

【全頁数】15

(21)【出願番号】特願2014-93581(P2014-93581)

(22)【出願日】2014年4月30日

(65)【公開番号】特開2015-208599(P2015-208599A)

(43)【公開日】2015年11月24日

【審査請求日】2017年4月25日

(73)【特許権者】

【識別番号】503313373

【氏名又は名称】株式会社ＡＺＥ

(74)【代理人】

【識別番号】110001771

【氏名又は名称】特許業務法人虎ノ門知的財産事務所

(72)【発明者】

【氏名】阪本 剛

(72)【発明者】

【氏名】浦山 博昭

(72)【発明者】

【氏名】畦元 将吾

【審査官】 松岡 智也

(56)【参考文献】

【文献】 特開２００５−３２３９２６（ＪＰ，Ａ）

【文献】 特開平０３−００９７３０（ＪＰ，Ａ）

【文献】 米国特許出願公開第２０１２／００３９５１８（ＵＳ，Ａ１）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ａ６１Ｂ ６／００−６／１４

(57)【特許請求の範囲】

【請求項1】

３次元医用画像データにより生成される被検体の複数の観察用画像に対して、ノイズを低減させるためのノイズ低減処理を実行する医用画像診断支援装置であって、
前記観察用画像を構成する画像構成要素に対応する信号値の所定の範囲ごとに、前記画像構成要素の個数を求める計算手段と、
前記計算手段により求められた前記範囲ごとの個数を、前記範囲ごとに補正する補正手段と、
前記補正手段により補正された前記個数に基づき、前記観察用画像ごとの前記ノイズ低減処理の実行強度を決定する決定手段と
を備えることを特徴とする医用画像診断支援装置。

【請求項2】

前記補正手段により補正された前記範囲ごとの前記個数の積算値を求める積算手段を更に備え、
前記決定手段は、前記積算手段で求められた前記積算値に基づいて、前記ノイズ低減処理の実行強度を決定すること
を特徴とする請求項１に記載の医用画像診断支援装置。

【請求項3】

前記計算手段、前記補正手段、および前記積算手段による各処理を、前記被検体の複数の観察用画像に対し実行し、当該複数の観察用画像の各々の積算値を求める複数積算値算出手段を備え、
前記決定手段は、前記複数積算値算出手段により求められた前記複数の観察用画像の各々の積算値に基づいて、前記ノイズ低減処理の実行強度を決定すること
を特徴とする請求項２に記載の医用画像診断支援装置。

【請求項4】

前記決定手段は、前記複数の観察用画像の各々について求められた前記積算値の中から基準となる基準積算値を選出し、当該基準積算値に対する、所定の観察用画像における前記積算値の比率を求め、当該比率に応じて、前記所定の観察用画像に対する前記ノイズ低減処理の実行強度を決定すること
を特徴とする請求項３に記載の医用画像診断支援装置。

【請求項5】

前記複数の観察用画像の各々について求められた各前記積算値のグラフを作成するグラフ作成手段を備え、
前記決定手段は、前記グラフ作成手段により作成された前記グラフを用いて、前記ノイズ低減処理の実行強度を決定すること
を特徴とする請求項３または４に記載の医用画像診断支援装置。

【請求項6】

信号値の所定の範囲ごとに、前記補正手段での補正のための重み付け係数を記憶する記憶手段を更に備え、
前記補正手段は前記記憶手段で記憶される前記重み付け係数に基づいて、前記範囲ごとに補正を行うこと
を特徴とする請求項１乃至５のいずれか１項に記載の医用画像診断支援装置。

【請求項7】

３次元医用画像データにより生成される被検体の複数の観察用画像に対して、ノイズを低減させるためのノイズ低減処理を実行する医用画像診断支援装置の制御方法であって、 前記観察用画像を構成する画像構成要素に対応する信号値の所定の範囲ごとに、前記画像構成要素の個数を求める計算ステップと、
前記計算ステップにより求められた前記範囲ごとの個数を、前記範囲ごとに補正する補正ステップと、
前記補正ステップにより補正された前記個数に基づき、前記観察用画像ごとの前記ノイズ低減処理の実行強度を決定する決定ステップと
を含むことを特徴とする医用画像診断支援装置の制御方法。

【請求項8】

３次元医用画像データにより生成される被検体の複数の観察用画像に対して、ノイズを低減させるためのノイズ低減処理を実行する医用画像診断支援装置で実行可能なプログラムであって、
前記医用画像診断支援装置を、
前記観察用画像を構成する画像構成要素に対応する信号値の所定の範囲ごとに、前記画像構成要素の個数を求める計算手段と、
前記計算手段により求められた前記範囲ごとの個数を、前記範囲ごとに補正する補正手段と、
前記補正手段により補正された前記個数に基づき、前記観察用画像ごとの前記ノイズ低減処理の実行強度を決定する決定手段
として機能させることを特徴とするプログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、医療分野等における画像診断を支援するための画像処理技術に関するものであり、特に、３次元医用画像データにより画像を生成する際に生じるノイズの低減処理を行う医用画像診断支援装置、方法およびプログラムに関する。

【背景技術】

【0002】

従来、Ｘ線ＣＴ（Computed Tomography）やＭＲＩ（Magnetic Resonance Imaging）等の画像診断装置により得られた３次元医用画像データ（３次元的な情報を担持したボリュームデータ）を、コンピュータを用いて２次元的に再構成することにより、被検体（患者）の任意の領域（例えば、胸部や腹部における体軸断面）に対応した観察用画像を作成し、その観察用画像を利用して、患者に病状を説明したり手術計画を立てたりすることが行われている。

【0003】

このような観察用画像には、様々な要因によりノイズが発生する。例えば、被検体の断面の観察用画像を作成する場合には、被検体の断面積（ボクセル換算の体積と考えてもよい）が大きいほどノイズ量が増える傾向がある。また、Ｘ線ＣＴによる観察用画像では、Ｘ線が通過し難い、高いＣＴ値を有する骨等の組織の割合が大きいほどノイズが増える傾向があることや、コンプトン散乱等の物理的現象の影響でノイズが生じることも知られている。

【0004】

Ｘ線ＣＴでは、Ｘ線量を増やすとノイズが低減する（例えば、Ｘ線量を２倍に増やすとノイズが２の平方根（２分の１乗）分の１に減少する）ことが知られているが、Ｘ線量を増やすと被検体の被曝線量も増えることになる。近年のＸ線ＣＴ装置においては、Ｘ線量を抑制しつつノイズの少ない観察用画像を作成することを目標としている。そのために、観察用画像の構成手法として、演算量は少ないが低線量下でのノイズが大きくなるフィルタ補正逆投影（FBP：Filtered Back Projection）法に替えて、演算量は多いが低線量下でもノイズを抑えることができる逐次近似再構成（IR：Iterative Reconstruction）法を応用した手法を用いるなどの対策を講じることにより、低線量下における画像ノイズの低減を図っている。

【0005】

一方、画像診断装置により構成された観察用画像（ボクセルデータの集合体としての画像データ）に対し、画像処理技術を用いた後処理としてのノイズ低減処理を行うことも近年一般的に行われている。このようなノイズ低減処理としては、移動平均フィルタや加重平均フィルタ等の線形フィルタや、メディアンフィルタやエッジ保存フィルタ等の非線形フィルタによる平滑化処理が一般的であるが、他にも種々の手法が知られている。

【0006】

例えば、被検体の体軸方向に沿った、互いに隣接する複数の観察用画像（体軸断面画像データ）がある場合に、その中の一の観察用画像のノイズパターンを、隣接する他の観察用画像を用いて抽出し、そのノイズパターンを一の観察用画像から除去する方法が知られている。また、Ｘ線ＣＴ装置により得られた複数の薄いスライス画像に対してそれぞれ３次元の画像フィルタ処理を施し、その後に複数の薄いスライス画像を束ねることによりノイズ低減を図る手法も提案されている（例えば、下記特許文献１を参照）。

【先行技術文献】

【特許文献】

【0007】

【特許文献1】特開２００５−１３１２８７号公報

【発明の概要】

【発明が解決しようとする課題】

【0008】

画像処理により良好なノイズ低減効果を得るためには、ノイズ低減処理を施す観察用画像に生じているノイズ量を的確に把握し、そのノイズ量に応じて、観察用画像に対しどの程度ノイズ低減処理を施すのか（以下、適宜「ノイズ低減処理の実行強度」と称する）を適正に決定することが重要とされている。

【0009】

しかしながら、観察用画像に生じているノイズ量を正確に把握することは難しく時間を要する。このため、ノイズ低減処理の実行強度を適正に決定しようとするあまり、画像ノイズ量を厳密に把握することに拘泥されると、ノイズ低減処理を実行するのに多大な時間を要することとなる。

【0010】

本発明は、３次元医用画像データにより生成される観察用画像に対し実行するノイズ低減処理の実行強度を、適切に決定することが可能な仕組みを提供することを目的とする。

【課題を解決するための手段】

【0011】

上記目的を達成するため、本発明に係る医用画像診断支援装置は、以下の特徴を備えている。

【0012】

本発明に係る医用画像診断支援装置は、３次元医用画像データにより生成される被検体の複数の観察用画像に対して、ノイズを低減させるためのノイズ低減処理を実行する医用画像診断支援装置であって、前記観察用画像を構成する画像構成要素に対応する信号値の所定の範囲ごとに、前記画像構成要素の個数を求める計算手段と、前記計算手段により求められた前記範囲ごとの個数を、前記範囲ごとに補正する補正手段と、前記補正手段により補正された前記個数に基づき、前記観察用画像ごとの前記ノイズ低減処理の実行強度を決定する決定手段とを備えることを特徴とする。

【0023】

上記「観察用画像」とは、紙面や画像モニタ面等の画像表示面に出力された画像の他、画像表示面に出力される前の、仮想の画像データ空間上に構築された、ボクセルやピクセル等の画像構成要素（画像データ）の集合体を含む概念である。

【0024】

また、上記「観察用画像を構成する複数の画像構成要素」とは、観察用画像の全領域を構成する全ての画像構成要素を意味する場合の他、観察用画像における特定の領域（例えば、観察対象となる領域や画像ノイズの特に多い領域）を構成する一部の画像構成要素を意味する場合を含む。

【0025】

また、上記「ノイズ低減処理」とは、画像処理技術を用いた後処理としてのノイズ低減処理を意味している。具体的には、空間領域における、移動平均フィルタや加重平均フィルタ（ガウシアンフィルタを含む）等の線形フィルタや、メディアンフィルタやエッジ保存フィルタ等の非線形フィルタによる平滑化処理が挙げられる。また、周波数領域において高周波成分を除去するためのローパスフィルタによる処理も他の例として挙げられる。

【0026】

また、上記「ノイズ低減処理の実行強度」とは、観察用画像に対し実行するノイズ低減処理の程度を示すもので、具体的には例えば、平滑化処理に用いられるフィルタにおける各係数の値の比率やフィルタの大きさ、ローパスフィルタにおける通過域や阻止域の大きさ、あるいは一の観察用画像に対し実行するノイズ低減処理の回数等が挙げられる。

【発明の効果】

【0027】

本発明によれば、３次元医用画像データにより生成される観察用画像に対し実行するノイズ低減処理の実行強度を、適切に決定することが可能となる。

【図面の簡単な説明】

【0028】

【図1】本発明の一実施形態に係る医用画像診断支援装置の概要を示す図（（Ａ）は概略構成図、（Ｂ）は演算処理装置の内部構成を示すブロック図）である。

【図2】図１に示す演算処理装置の機能的構成を示すブロック図である。

【図3】ＣＴ値に応じて設定する各群と群別ボクセル数と重みづけ係数との対応関係を示す図である。

【図4】群別ボクセル数の一例を示す図（（Ａ）は変更前の群別ボクセル数の度数分布曲線、（Ｂ）は変更後の群別ボクセル数の度数分布曲線）である。

【図5】積算値比較グラフの一例を示す図である。

【図6】正規化した積算値比較グラフの一例を示す図である。

【図7】ノイズ低減処理に使用するガウシアンフィルタ（（Ａ）は３×３、（Ｂ）は５×５、（Ｃ）７×７）を例示する図である。

【図8】ノイズ低減処理を実行した観察用画像（（Ａ）は胸部の体軸断面、（Ｂ）は大腿部の体軸断面）を例示する図である。

【図9】他の観察用画像（骨盤付近の体軸断面）を示す図（（Ａ）はノイズ低減処理を実行したもの、（Ｂ）はノイズ低減処理を実行していないもの）である。

【図10】ノイズ低減処理を実行していない観察用画像（（Ａ）は肩付近の体軸断面、（Ｂ）は腸骨付近の体軸断面）を例示する図である。

【図11】ノイズ低減処理を実行していない他の観察用画像（（Ａ）は膝関節付近の体軸断面、（Ｂ）は肝臓付近の体軸断面）を例示する図である。

【発明を実施するための形態】

【0029】

以下、本発明の実施形態について、上述の図面を参照しつつ詳細に説明する。まず、図１、２を参照しながら本発明の一実施形態に係る医用画像診断支援装置の構成について説明する。

【0030】

本実施形態の医用画像診断支援装置１は、図１（Ａ）に示すように、コンピュータ等からなる演算処理装置１０と、キーボード２１やマウス２２等からなる入力装置２０と、液晶表示パネル等からなる画像表示部３１を有する画像表示装置３０とを備えてなる。この医用画像診断支援装置１は、Ｘ線ＣＴ装置やＭＲＩ装置等の画像診断装置により得られた３次元医用画像データにより生成される、被検体の体軸断面等の観察用画像（画像診断装置に付属する画像処理部により構成されたものを含む）に対し、画像処理技術を用いた後処理としてのノイズ低減処理を行うワークステーションとして構成されている。

【0031】

上記演算処理装置１０は、図１（Ｂ）に示すように、画像処理等における各種演算を実行するＣＰＵ（Central Processing Unit）１０１と、処理プログラムや処理すべき画像データ等が置かれる、ＤＲＡＭ（Dynamic Random Access Memory）等からなるメインメモリ１０２と、ＣＰＵ１０１やメインメモリ１０２等へのデータの入出力を制御するチップセット１０３とを備えている。また、ＣＰＵ１０１から送られた画像データや処理命令に基づき表示用画像を構成するＧＰＵ（Graphic Processing Unit）１０５および該ＧＰＵ１０５により構成された表示用画像を記憶するＶＲＡＭ（Video Random Access Memory）１０６を有するグラフィックスボード１０４と、ハードディスク等からなる記憶装置１０７と、これら各構成要素間におけるデータ等のやり取りを仲介するデータバス１０８とを備えている。さらに、当該演算処理装置１０と、入力装置２０および画像表示装置３０との間におけるデータ等のやり取りを仲介するインターフェース１０９を備えている。

【0032】

また、この演算処理装置１０は、機能的な観点から説明すれば、図２に示すように、画像データ取得部１１、画像データ解析部１２および画像生成部１３を備えている。なお、図２に示す各ブロックは、図１（Ｂ）に示すＣＰＵ１０１、メインメモリ１０２、グラフィックスボード１０４、記憶装置１０７等のハードウェアや、ＣＰＵ１０１やＧＰＵ１０５が実行する各種プログラム（本発明に係る医用画像診断支援プログラムを含む）等により実現される各機能を具象化し、構成要素として示したものである。

【0033】

上記画像データ取得部１１は、Ｘ線ＣＴ装置やＭＲＩ装置等の各種の医用画像診断装置により得られた種々の３次元医用画像データを取得するものである。なお、以下では、画像データ取得部１１が、マルチスライスＸ線ＣＴ装置により得られた、被検体の全身（頭部から足部まで）のＣＴ画像データを取得する場合を例にとって説明する。

【0034】

上記画像データ解析部１２は、画像データ取得部１１により取得されたＣＴ画像データに基づき、各種の解析や計測を行うものであり、画像構成要素分別部１２１、画像構成要素計数部１２２、画像構成要素数変更部１２３、画像構成要素数積算部１２４、複数画像積算値算出部１２５、ノイズ低減処理実行強度決定部１２６およびノイズ低減処理実行部１２７を備えている。

【0035】

画像構成要素分別部１２１（画像構成要素分別手段に相当）は、取得されたＣＴ画像データにより構築される画像データ空間（以下、適宜「ボクセル空間」と称する）において、被検体の所定の領域（本実施形態では、被検体の頭部から足部までの各体軸断面のうちの一断面）に対応する所定の観察用画像（一断面画像）を構成する（本実施形態では、観察用画像中の被検体領域を構成する）複数のボクセルを、各ボクセルが担持するＣＴ値（信号値に相当）の大きさに応じて複数の群（詳しくは後述する）に分別するように構成されている。

【0036】

画像構成要素計数部１２２（画像構成要素計数手段に相当）は、画像構成要素分別部１２１により上記複数の群に分別されたボクセルの、群別の個数（各群に属するボクセルの個数、以下、適宜「群別ボクセル数」と称する）をそれぞれ求めるように構成されている。

【0037】

画像構成要素数変更部１２３（画像構成要素数変更手段に相当）は、画像構成要素計数部１２２により求められた上記群別ボクセル数を、複数の群の各々に属する各ボクセルが担持するＣＴ値の大きさに応じて変更する（変更方法の詳細については後述する）ように構成されている。

【0038】

画像構成要素数積算部１２４（画像構成要素数積算手段に相当）は、画像構成要素数変更部１２３による変更後の群別ボクセル数の積算値（変更後の各群別ボクセル数を合算した値、以下、適宜「変更後ボクセル数積算値」と称する）を求めるように構成されている。

【0039】

複数画像積算値算出部１２５（複数画像積算値算出手段に相当）は、画像構成要素分別部１２１、画像構成要素計数部１２２、画像構成要素数変更部１２３および画像構成要素数積算部１２４による各処理を、被検体の複数の領域（本実施形態では、被検体の頭部から足部までの各体軸断面）の各々に対応した複数の観察用画像に対し実行させ、その複数の観察用画像の各々について上記変更後ボクセル数積算値を求めるように構成されている。

【0040】

ノイズ低減処理実行強度決定部１２６（ノイズ低減処理実行強度決定手段に相当）は、複数画像積算値算出部１２５の指示の下、画像構成要素数積算部１２４により求められた各観察画像の変更後ボクセル数積算値に基づき、各観察用画像に施すノイズ低減処理の実行強度を画像ごとに決定する（実行強度の決定方法の詳細については後述する）ように構成されている。

【0041】

ノイズ低減処理実行部１２７（ノイズ低減処理実行手段に相当）は、ノイズ低減処理実行強度決定部１２６により画像ごとに決定されたノイズ低減処理の実行強度に基づき、ノイズ低減処理の実行対象となる観察用画像に対しノイズ低減処理を実行するように構成されている。

【0042】

一方、上記画像生成部１３は、画像データ取得部１１により取得されたＣＴ画像データや、画像データ解析部１２により求められた各種の解析または計測結果に基づき、画像診断支援に供する種々の画像を生成するものであり、積算値グラフ作成部１３１および観察用画像作成部１３２を備えている。

【0043】

積算値グラフ作成部１３１（積算値グラフ作成手段に相当）は、被検体の頭部から足部までの各体軸断面に対応した複数の観察用画像の各々について求められた変更後ボクセル数積算値の大小関係を示すグラフ（以下、適宜「積算値比較グラフ」と称する）を作成するように構成されている。

【0044】

観察用画像出力部１３２は、ノイズ低減処理が施された観察用画像を出力用の画像に変換し、画像表示装置３０の画像表示部３１に表示するように構成されている。

【0045】

次に、本発明に係る医用画像診断支援方法の一実施形態について、主に図３〜図１１を参照しながら説明する。この医用画像診断支援方法は、本発明の医用画像診断支援プログラムに基づき、上記医用画像診断支援装置１により実行されるものであり、大別すると以下の〈１〉〜〈９〉の手順からなる。

【0046】

〈１〉 被検体の頭部から足部までの各体軸断面のうちの一断面（例えば、胸部の一断面）に対応する観察用画像（以下、適宜「第１観察用画像」と称する）を構成する（詳しくは、第１観察用画像中の被検体領域を構成する）複数のボクセル（以下、適宜「対象ボクセル」と称する）を、各ボクセルが担持するＣＴ値の大きさに応じて複数の群に分別する（画像構成要素分別ステップ）。この複数のボクセルの分別は、詳細には例えば、以下の手順で行われる。

【0047】

まず、ＣＴ値に応じて複数の群を設定する。本実施形態では、図３に示すように、第１群、第２群および第３群の３つの群を設定する。第１群は、空気に該当する領域のボクセルが属するようにＣＴ値の範囲を−５００未満に設定している。第２群は、脂肪や水、軟組織に該当するボクセルが属するようにＣＴ値の範囲を−５００以上１００未満に設定しており、第３群は、骨に該当する領域のボクセルが属するようにＣＴ値の範囲を１００以上に設定している。

【0048】

次に、第１観察用画像における上記対象ボクセルの全てを、各ボクセルが担持するＣＴ値の大きさに応じて、第１群、第２群および第３群の各々に分別する。例えば、ＣＴ値が４０３のボクセルは第３群に分別し、ＣＴ値が−８４２のボクセルは第１群に分別し、ＣＴ値が２８のボクセルは第２群に分別するという具合である。

【0049】

〈２〉 第１群〜第３群に分別されたボクセルの群別の個数（群別ボクセル数）を求める（画像構成要素計数ステップ）。ここでは、図３に示すように、第１群の群別ボクセル数がＮ₁、第２群の群別ボクセル数がＮ₂、第３群の群別ボクセル数がＮ₃であるとして説明を続ける。なお、対象ボクセルの分別時に、各群の群別ボクセル数を計数するようにしてもよい。

【0050】

〈３〉 上記〈２〉の手順により求められた各群の群別ボクセル数Ｎ₁、Ｎ₂、Ｎ₃を、各群に属するボクセルが担持するＣＴ値の大きさに応じて変更する（画像構成要素数変更ステップ）。このＣＴ値の変更は、詳細には例えば、以下の手順で行われる。なお、以下の手順による各処理は、上述の画像構成要素数変更部１２３により実行される。

【0051】

まず、上記〈１〉の手順により設定された第１群、第２群および第３群の各々に対応する重みづけ係数Ｋ₁、Ｋ₂、Ｋ₃を設定する。この重みづけ係数Ｋ₁、Ｋ₂、Ｋ₃は、上記〈２〉の手順により求められた各群の群別ボクセル数Ｎ₁、Ｎ₂、Ｎ₃に乗じて変更するためのもので、各群に対応づけたＣＴ値の範囲と画像ノイズの発生度合に応じて値が決定される。

【0052】

例えば、本実施形態では、第２群のＣＴ値範囲に属する脂肪や水、軟部組織等の占める割合が大きい画像中ではノイズ発生量が比較的少ないことを考慮し、第２群の重みづけ係数Ｋ₂の値を１に設定する。一方、第１群のＣＴ値範囲に属する空気等の占める割合が増えるとノイズ量がかなり増えることを考慮し、第１群の重みづけ係数Ｋ₁の値を２に設定する。また、第３群のＣＴ値範囲に属する骨や歯、金属等の占める割合が増えるとノイズ量がさらに大きく増えることを考慮し、第３群の重みづけ係数Ｋ₃の値を３に設定する。

【0053】

次に、下式（１）〜（３）に示すように、群別ボクセル数Ｎ₁、Ｎ₂、Ｎ₃に重みづけ係数Ｋ₁、Ｋ₂、Ｋ₃をそれぞれ乗じ、変更後の群別ボクセル数Ｎ₁´、Ｎ₂´、Ｎ₃´を算出する。

【0054】

Ｎ₁´＝Ｋ₁・Ｎ₁（＝２Ｎ₁） ・・・ （１）
Ｎ₂´＝Ｋ₂・Ｎ₂（＝Ｎ₂） ・・・ （２）
Ｎ₃´＝Ｋ₃・Ｎ₃（＝３Ｎ₃） ・・・ （３）

【0055】

なお、上述の変更後の群別ボクセル数Ｎ₁´、Ｎ₂´、Ｎ₃´を算出する過程を、第１観察画像における対象ボクセルの度数分布曲線（階級幅＝１）を用いて表したものが、図４である。図４（Ａ）は、重みづけ係数Ｋ₁、Ｋ₂、Ｋ₃を乗じる前、すなわち変更前の群別ボクセル数Ｎ₁、Ｎ₂、Ｎ₃の度数分布を示し、図４（Ｂ）は、重みづけ係数Ｋ₁、Ｋ₂、Ｋ₃を乗じた場合の、変更後の群別ボクセル数Ｎ₁´、Ｎ₂´、Ｎ₃´の度数分布を示している。このような度数分布曲線を用いた表示を行うことは必須ではないが、医師等に変更後の群別ボクセル数Ｎ₁´、Ｎ₂´、Ｎ₃´を算出する過程を分かり易く示すために行ってもよい。また、表示された度数分布を見た医師等が、重みづけ係数Ｋ₁、Ｋ₂、Ｋ₃を補正し得るようにしてもよい。

【0056】

〈４〉 上記〈３〉の手順により算出された、変更後の群別ボクセル数Ｎ₁´、Ｎ₂´、Ｎ₃´の積算値（変更後ボクセル数積算値）ＴＮを、下式（４）により求める（画像構成要素数積算ステップ）。

【0057】

ＴＮ＝Ｎ₁´＋Ｎ₂´＋Ｎ₃´ ・・・ （４）

【0058】

〈５〉 上記〈１〉、〈２〉、〈３〉、〈４〉の手順による各処理を、第１観察用画像以外の、被検体の頭部から足部までの各体軸断面に対応した複数の観察用画像（例えば、Ｘ線ＣＴ装置の各スライス位置に対応した観察用画像）に対し実行し、その複数の観察用画像の各々について上記変更後ボクセル数積算値を求める（複数画像積算値算出ステップ）。

【0059】

〈６〉 上記〈５〉の手順により複数の観察用画像の各々について求められた各変更後ボクセル数積算値の大小関係を示すグラフ（積算値比較グラフ）を作成する（積算値グラフ作成ステップ）。

【0060】

図５に、このような積算値比較グラフの一例を示す。図５に示す積算値比較グラフは、各観察用画像の断面位置（体軸方向の位置）を横軸にとり、各観察用画像における変更後ボクセル数積算値を縦軸にとって、体軸方向における変更後ボクセル数積算値の変化状態を曲線（実線による曲線）で表したものである。

【0061】

なお、図５では、各観察用画像における上記対象ボクセル（被検体領域を構成するボクセル）の個数の変化状態を、２点鎖線による曲線（以下、適宜「参考曲線」と称する）で表示している。この参考曲線は、体軸方向の各位置における被検体領域全体の断面積値をそのまま示したものに相当する。一方、本実施形態における変更後ボクセル数積算値は、観察用画像における、画像ノイズが発生し易い空気や骨等の被検体領域の断面積値を、画像ノイズの発生量において基準となる脂肪や水、軟部組織等の被検体領域の断面積値に換算した場合の被検体領域全体の断面積値に該当する。したがって、図５における実線による曲線（積算値比較グラフ）は、このような断面積値の換算をした場合の、体軸方向の各位置における被検体領域全体の断面積値を示したものに相当する。

【0062】

〈７〉 上記〈５〉の手順により求められた各変更後ボクセル数積算値に基づき、ノイズ低減処理の実行強度を決定する（ノイズ低減処理実行強度決定ステップ）。この実行強度の決定は、上記〈６〉の手順により作成された積算値比較グラフを用いて、詳細には例えば、以下の手順で行われる。なお、以下の手順による各処理は、上述のノイズ低減処理実行強度決定部１２６により実行される。

【0063】

まず、上記積算値比較グラフにおいて、被検体の頭部から足部までの各体軸断面に対応した複数の観察用画像の各々について求められた各変更後ボクセル数積算値の中から基準となる基準積算値を選出する。本実施形態では、各変更後ボクセル数積算値のうち最大値となるものを、基準積算値として選出する（図５を参照）。

【0064】

次いで、上記積算値比較グラフを正規化する。具体的には、変更後ボクセル数積算値の比較を容易とするために、基準積算値に「１」を対応づける正規化を行う（基準積算値以外の変更後ボクセル数積算値は、基準積算値との比率に応じた数値に対応づける）。このような正規化を行った積算値比較グラフ（以下「正規化した積算値比較グラフ」と称する）を図６に示す。

【0065】

次に、ノイズ低減処理を施す観察用画像（ここでは、第１観察用画像とする）における変更後ボクセル数積算値の上記基準積算値に対する比率を、上記正規化した積算値比較グラフから求め（ここでは、その比率を０．８５とする）、その比率に応じて、第１観察用画像に対するノイズ低減処理の実行強度を決定する。本実施形態では、ノイズ低減処理の実行強度を、例えば、上記比率が０以上０．４未満の範囲に属する場合は強度１に設定し、０．４以上０．７未満の範囲に属する場合は強度２に設定し、０．７以上１．０以下の範囲に属する場合は強度３（数値が大きいほど実行強度が大きい）に設定する（強度設定の詳細については後述する）。

【0066】

なお、このような積算値比較グラフや正規化した積算値比較グラフを作成し表示することは必須ではないが、医師等にノイズ低減処理の実行強度を決定する過程を分かり易く示すために行ってもよい。また、医師等が、決定されたノイズ低減処理の実行強度を補正し得るようにしてもよい。

【0067】

〈８〉 上記〈７〉の手順により決定された実行強度に基づき、第１観察用画像に対しノイズ低減処理を実行する（ノイズ低減処理実行ステップ）。本実施形態では、ノイズ低減処理を、例えば、２次元のガウシアンフィルタを用いた平滑化処理により行う。ガウシアンフィルタは、カーネル（オペレータ、マトリックス等とも称される）における重み係数の分布を、下式（５）で表される正規分布に近似させた加重平均フィルタである。

【0068】

ｆ（ｘ，ｙ）＝（１／２πσ²）exp｛−（ｘ²＋ｙ²）／２σ²｝ ・・・ （５）
ここで、ｘ，ｙは座標値、σは分散を表す。

【0069】

本実施形態においては、ノイズ低減処理の実行強度（ガウシアンフィルタによる平滑化の度合）を、例えば、上記（５）中のσの値を変えてカーネルの重み係数を変更することにより設定する。具体的には例えば、ノイズ低減処理の実行強度が、強度１の場合はσの値を０．８に設定し、強度２の場合はσの値を１．０に設定し、強度３の場合はσの値を１．３に設定して、カーネルの重み係数を決定する。なお、σの値に応じてカーネルの大きさを変更してもよいし、σの値によらずカーネルの大きさを一定（例えば、５×５）としてもよい。

【0070】

図７に、実行強度に応じてカーネルの大きさを変更する場合のガウシアンフィルタの例を示す。図７（Ａ）に示す３×３のガウシアンフィルタは強度１のものでσ≒０．８となるように重み係数が設定されている。同様に、図７（Ｂ）に示す５×５のガウシアンフィルタは強度２のものでσ≒１．０となるように重み係数が設定され、図７（Ｃ）に示す７×７のガウシアンフィルタは強度３のものでσ≒１．３となるように重み係数が設定されている。なお、図７に示す各カーネル中の数値は、実際の重み係数の値を分数で表記する場合の分子の値のみを示しており、実際の係数値は、カーネル中の各数値に各カーネルの左に記す括弧内の数値を乗じた値となる。

【0071】

上記第１観察用画像の場合、上記〈７〉の手順により、変更後ボクセル数積算値の基準積算値に対する比率（０．８５）からノイズ低減処理の実行強度が強度３に設定されるので、強度３のガウシアンフィルタを用いてノイズ低減処理（平滑化処理）が実行される。

【0072】

〈９〉 上記〈８〉の手順によりノイズ低減処理が施された観察用画像を、出力環境に応じた諧調処理等を経て出力用の画像に変換し、画像表示装置３０の画像表示部３１に表示する（観察用画像出力ステップ）。

【0073】

図８に、本実施形態によるノイズ低減処理を実行した観察用画像（出力画像）の例を示す。図８（Ａ）に示すのは、ノイズ低減処理を実行した胸部の観察用画像（上述の第１観察用画像に相当する体軸断面画像）である。このような胸部の観察用画像（ノイズ低減処理前）では、肺野部の空気領域の占める割合が大きく、被検体領域のＣＴ値の平均値は低くなる。このため、画像ノイズの発生量は少なくなるとも考えられるが、実際には、ノイズ発生量は多くなる。また、肺野の場合、末梢血管や気管支などの構造を読影するためにエッジ強調フィルタが用いられることが多く、その場合ノイズ成分も強調されて画像に表示されることになる。このため、ノイズ量の十分な低減が必要になる。本実施形態では上述したように、空気領域のＣＴ値範囲を含む上記第１群の重みづけ係数の値を高く設定することにより、変更後の群別ボクセル数が多くなるようにし、ノイズ低減処理の実行強度が高レベル（強度３）に設定されるようにしている。このため、本実施形態によるノイズ低減処理を施すことにより良好なノイズ低減効果が得られる。

【0074】

図８（Ｂ）に示すのは、本実施形態によるノイズ低減処理を実行した大腿部の観察用画像（体軸断面画像）である。このような大腿部の観察用画像（ノイズ低減処理前）では、上記胸部のものと比べ被検体領域の断面積が小さく骨部の領域も多くないため、画像ノイズの発生量は比較的少なくなる。本実施形態を適用すると、ノイズ低減処理の実行強度が中間レベル（強度２）に設定され、これにより良好なノイズ低減効果が得られる。

【0075】

図９に、本実施形態によるノイズ低減処理を実行した観察用画像（出力画像）の他の例を示す。図９（Ａ）に示すのは、ノイズ低減処理を実行した骨盤付近の観察用画像（体軸断面画像）であり、図９（Ｂ）に示すのは、ノイズ低減処理を実行していない骨盤付近の観察用画像である。このような骨盤付近の観察用画像（ノイズ低減処理前）では、被検体領域の断面積が大きく骨部の領域も多くなるため、画像ノイズの発生量は多くなる。本実施形態を適用すると、ノイズ低減処理の実行強度が高レベル（強度３）に設定され、これにより良好なノイズ低減効果が得られる。

【0076】

なお、以下、体軸断面位置による画像ノイズの発生量の傾向について説明する。図１０は、ノイズ低減処理を実行していない肩付近の観察用画像（同図（Ａ））とノイズ低減処理を実行していない腸骨付近の観察用画像（同図（Ｂ））を示している。これらの観察用画像では、被検体領域の断面積が大きく骨部の領域も多くなるため、画像ノイズの発生量は多くなる。本実施形態を適用する場合、ノイズ低減処理の実行強度が高レベル（強度３）に設定される。

【0077】

図１１は、ノイズ低減処理を実行していない膝関節付近の観察用画像（同図（Ａ））とノイズ低減処理を実行していない肝臓付近の観察用画像（同図（Ｂ））を示している。膝関節付近の観察用画像では、被検体領域の全体に占める骨部の割合が比較的高いが被検体領域の断面積が小さいため、画像ノイズの発生量は少なくなる。本実施形態を適用する場合、ノイズ低減処理の実行強度が低レベル（強度１）に設定される。肝臓付近の観察用画像では、被検体領域の全体に占める骨部の割合は少ないが被検体領域の断面積が大きいため、画像ノイズの発生量は比較的多くなる。本実施形態を適用する場合、ノイズ低減処理の実行強度が中レベル（強度２）から高レベル（強度３）に設定される。

【0078】

以上、本発明の実施形態について説明したが、本発明は上記実施形態に限定されるものではなく、種々に実施の態様を変更することが可能である。

【0079】

例えば、上述の実施形態では、観察用画像における対象ボクセルを分別する群の数を３としているが、２以上の任意の数に設定可能である。例えば、上述の第２群をさらに３つの群（ＣＴ値の範囲が−５００〜０、０〜５０、５０〜１００の各群）に分け、第３群をさらに２つの群（ＣＴ値の範囲が１００〜１０００の群と１０００以上の群）に分けることが挙げられる。また、各群に対応づけるＣＴ値の範囲や各群に適用する重みづけ係数の値についても上記実施形態とは異なるように適宜設定可能である。

【0080】

また、上述の実施形態では、ノイズ低減処理の実行強度を３段階（強度１〜３）に設定しているが、例えば、強度１〜１０の１０段階を設定したり、強度１〜５の５段階を設定したりするなど、２以上の任意の数の段階を設定可能である。

【0081】

さらに、上述の実施形態では、複数の観察用画像の各々について求めた各変更後ボクセル数積算値の中から、最大値となるものを基準積算値として選出しているが、別の条件で基準積算値を選出してもよい。例えば、変更後ボクセル数積算値が所定の数値範囲に属するもの（複数ある場合はそのうちの１つ）を、基準積算値として選出することが挙げられる。また、基準積算値との比率に基づき設定するノイズ低減処理の実行強度についても、基準積算値の設定の仕方に応じて適宜変更することが可能である（例えば、基準積算値に対応する実行強度が低レベルや中レベルになることもある）。

【0082】

また、上述の実施形態では、ガウシアンフィルタを用いた平滑化処理をノイズ低減処理として実行しているが、移動平均フィルタ等の別の線形フィルタやメディアンフィルタ等の非線形フィルタによる平滑化処理をノイズ低減処理として実行することや、周波数領域においてローパスフィルタにより高周波成分を除去する処理をノイズ低減処理として実行することも可能である。

【0083】

また、上記実施形態では、Ｘ線ＣＴ装置により取得されたＣＴ画像データを用いて、体軸断面画像を観察用画像として作成する場合を例にとって説明しているが、本発明は他の場合にも適用可能である。例えば、ＤＳＡ（Digital Subtraction Angiography）装置により取得された画像データを用いて血管画像を観察用画像として作成する場合や、ＭＲＩ装置により取得された画像データを用いて任意断面の画像を観察用画像として作成する場合が挙げられる。

【符号の説明】

【0084】

１ 医用画像診断支援装置
１０ 演算処理装置
２０ 入力装置
２１ キーボード
２２ マウス
３０ 画像表示装置
３１ 画像表示部
Ｎ₁、Ｎ₂、Ｎ₃ 群別ボクセル数
Ｎ₁´、Ｎ₂´、Ｎ₃´ 変更後の群別ボクセル数
Ｋ₁、Ｋ₂、Ｋ₃ 重みづけ係数

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】