WO2006129497A1

WO2006129497A1 - 画像処理方法および画像処理装置

Info

Publication number: WO2006129497A1
Application number: PCT/JP2006/309999
Authority: WO
Inventors: Daisuke Kaji
Original assignee: Konica Minolta Medical & Graphic, Inc.
Priority date: 2005-05-31
Filing date: 2006-05-19
Publication date: 2006-12-07
Also published as: JP4844560B2; US20090136110A1; US8036443B2; CN101184441A; JPWO2006129497A1; CN100562290C

Abstract

　適切な条件で画像処理を実行する。　放射線画像に含まれる骨部あるいは軟部を重要領域として検出する重要領域検出手段１２０と、重要領域の画素に所定の重みを持たせる重み画像を生成する重み画像生成手段１３０と、重み画像の重みと放射線画像の画素値とを乗算し、この乗算結果により重み付きヒストグラムを算出するヒストグラム算出手段１４０と、重み付きヒストグラムを所定の評価関数により評価し、評価結果が最大となるシフト値を算出する評価手段１５０と、シフト値を得る評価関数の最大値に対応する放射線画像の画素値で所定の処理結果が得られるように画像処理条件を決定する画像処理条件決定手段１６０と、画像処理条件に基づいて画像処理を実行する画像処理手段１７０と、を有する。

Description

明細書

画像処理方法および画像処理装置

技術分野

[0001] 本発明は放射線画像を処理する際の画像処理方法および画像処理装置に関し、さら〖こ詳しくは、診断等に適した放射線画像を得ることができる画像処理方法および画像処理装置に関する。

背景技術

[0002] 近年、放射線画像を直接デジタル画像として撮影できる装置が開発されてヽる。たとえば、被写体に照射された放射線量を検出し、その検出量に対応して形成される放射線画像を電気信号として得る装置としては、輝尽性蛍光体を用いたディテクタを用いる方法が特開昭 55-12429号公報、特開昭 63-189853号公報など、多数開示されている。

[0003] このような装置では、シート状の基板に輝尽性蛍光体を塗布、あるいは蒸着等によつて固着したディテクタに、 Vヽつたん被写体を透過した放射線を照射して輝尽性蛍光体に放射線を吸収させる。

[0004] その後、この輝尽性蛍光体を光または熱エネルギーで励起することにより、この輝尽性蛍光体が上記吸収によって蓄積している放射線エネルギ—を蛍光として放射させ、この蛍光を光電変換して画像信号を得るようにしてヽる。

[0005] 一方、照射された放射線の強度に応じた電荷を光導電層に生成し、生成された電荷を二次元的に配列された複数のコンデンサに蓄積し、それら蓄積された電荷を取り出すことにより得られる放射線画像検出装置が提案されて、る。

[0006] このような放射線画像検出装置では、フラットパネルディテクタ (FPD)と呼ばれるものを使用している。この種の FPDは、特開平 9-90048号公報に記載されているよう〖こ、照射された放射線強度に応じた蛍光を発する蛍光体と、蛍光体から発する蛍光を直接または縮小光学系を介して受光して光電変換を行うフォトダイオードや CCDのような光電変換素子の組み合わせによって実現されるものが知られている。

[0007] また特開平 6-342098号公報に記載されているように、照射された放射を直接電荷に変換するものも知られて、る。

[0008] これらの放射線画像検出装置では、取得した画像に対して、診断に適した画像となるよう、階調変換処理やエッジ強調処理などの各種の画像処理を行うことが一般的に行われている。

[0009] また、このようにして得られた画像データに基づく放射線画像を表示したり出力する場合、撮影条件の変動を受けることなく見やすヽ放射線画像となるような画像処理も行われている。

[0010] このため、放射線画像の所望の領域内の画像データの累積ヒストグラムを求めたのち、この累積ヒストグラムが所定の割合となるデータレベルを基準信号値と設定して画像処理を行うことが、以下の特許文献 1に記載されてヽる。

[0011] また、高信号値領域と低信号値領域の分布状態に基づいて、画像処理条件を定めて、適切な画像処理を行うことが、以下の特許文献 2に記載されている。

[0012] さらに、同様な処理が以下の特許文献 3や特許文献 4にも記載されている。

特許文献 1 :特開平 6— 61325号公報 (第 1頁、図 1)

特許文献 2 :特開 2000— 157518号公報 (第 1〜5頁、図 4)

特許文献 3：特開 2001— 238868号公報 (第 1頁、図 1)

特許文献 4:特開平 7— 271972号公報 (第 1頁、図 1)

発明の開示

発明が解決しょうとする課題

[0013] ところで、放射線画像では、被写体の撮影部位によって高濃度領域 (放射線透過量の大き!、領域)と低濃度領域 (放射線透過量の小さ!、領域)の割合が大きく異なる。また、被写体の状態、例えば患者の呼吸状態によって肺野内の濃度が著しく異なるちのとされる。

[0014] このため、放射線画像の関心領域 (ROI)などの所望の領域内の画像データの累積ヒストグラムを求めたのち、この累積ヒストグラムが所定の割合となるデータレベルを基準信号値と設定して画像処理を行う場合 (特開平 6— 61325号公報）には、高濃度領域と低濃度領域の割合の違いによっては、診断等に適した見やすい放射線画像を得ることができな、恐れがある。 [0015] 例えば、基準信号値が所定の濃度となるように階調処理を行った場合、低濃度領域が支配的であるときには基準信号値が低い値となるので全体的に高濃度な画像となり、高濃度領域が支配的であるときには基準信号値が高い値となるので全体的に低濃度な画像となる。

[0016] しかし、例えば肋骨のような部位では上部肋骨部（高濃度領域の多い肺野部分が含まれる)や下部肋骨部 (高濃度領域の少な!、腹部部分が含まれる)を撮影して診断を行う場合でも、診断で必要とされる部分は高濃度領域と低濃度領域の両方にまたがっており、画像全体の濃度が低濃度あるいは高濃度に偏ることは、診断を行う上で好ましくない。

[0017] すなわち、従来は、画像処理のパラメータを決定するために様々な手法が提案されて!、るが、その手法は画像内の特定領域（関心領域 (ROI)など）の信号値やヒストグラム力決定される信号値あるいはパラメータとして一意的に与えるなどのように決定されているため、得られる領域やパラメータが不適切な場合、診断を行うには不適切な画像となり、再調整が必要となる場合があった。

[0018] すなわち、従来技術においては、各被写体あるいは被写体の部位毎に、専用の処理ゃ専用の関心領域を設定する必要があり、処理拡張性に難点があった。例えば、部位情報などの画像処理のノラメータが間違った場合には、この影響を受けてしま

Vヽ、診断を行うには不適切な画像となってしまう問題を有して!/ヽた。

[0019] 又、画像内における関心領域を自動認識する場合に於いては、関心領域の認識精度に、診断対象部分の仕上がり濃度 (処理後の濃度)が大きく依存する結果、誤つた認識を行った時には、診断に使用できず、再調整が必要となる場合があった。

[0020] 一般に、比較的規模の大きな施設では、専任の画像処理担当者をおき、当該担当者力 CR (Computed Radiography)撮影装置や FPD等で生成された各放射線画像に対し、一定の画像品質を満足しているか否か、すなわち診断に提供できる画像である力否かのチェックを行っている。したがって、設定不良に起因する画質不良が発生していても、濃度'コントラスト等の修正処理を加えて上記一定の画像品質を満足するものとして力ゝら、次工程である読影 (診断）に提供するシステムが採用されている [0021] 一方、開業医やクリック等の小規模施設では、医師一人で、問診 ·撮影条件設定 •画像処理 ·診断 (読影)等の全てを行うことが多ぐ設定不良や自動部位認識不良等に起因する画質不良 (診断に使用できない画像)が発生した場合には、上記大規模施設のような専任の画像処理担当者と同レベルの修正処理が行われることは期待できない。

[0022] 従って、専任の画像処理担当者のいない小規模施設向けの画像処理としては、診断を行う上で高画質の画像を提供することも可能だが、画質不良が発生して再処理する確率が高い処理よりも、中庸品質ではあるが再処理が不要となるロバスト性の高 V、画像処理の方が好ま、。

[0023] 本発明は以上のような課題に鑑みてなされたものであって、関心領域設定方式などにおけるロバスト性の低下を防止し、適切な条件で画像処理することが可能な画像処理方法および画像処理装置を実現することを目的とする。

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

[0024] すなわち、前記した課題の解決は、以下に列記する発明により解決される。

[0025] 請求の範囲第 1項に記載の発明は、被写体を透過した放射線の照射線量に応じた信号を有する放射線画像に対し、診断に適した画像を得るための画像処理を行う画像処理方法であって、前記放射線画像に含まれる骨部あるいは軟部を重要領域として検出する重要領域検出ステップと、前記重要領域の画素に所定の重みを持たせる重み画像を生成する重み画像生成ステップと、前記重み画像の重みと前記放射線画像の画素値とを乗算し、この乗算結果により重み付きヒストグラムを算出するヒストグラム算出ステップと、前記重み付きヒストグラムを所定の評価関数により評価し、評価結果が最大となるシフト値を算出する評価ステップと、前記シフト値を得る前記評価関数の最大値に対応する前記放射線画像の画素値で所定の処理結果が得られるように画像処理条件を決定する画像処理条件決定ステップと、前記画像処理条件に基づ!/ヽて画像処理を実行する画像処理ステップと、を有することを特徴とする画像処理方法である。

[0026] 請求の範囲第 2項に記載の発明は、前記画像処理ステップは、階調処理である、ことを特徴とする請求の範囲第 1項に記載の画像処理方法である。 [0027] 請求の範囲第 3項に記載の発明は、前記評価関数の評価特性の最大値を、前記放射線画像に含まれる骨部の画素値の範囲内に有する、ことを特徴とする請求の範囲第 1項または第 2項に記載の画像処理方法である。

[0028] 請求の範囲第 4項に記載の発明は、前記評価関数の評価特性の最大値を、前記放射線画像に含まれる骨部の画素値の範囲内に有すると共に、該評価特性の別の山となる部分を前記骨部の画素値より大きい画素値に有する、ことを特徴とする請求の範囲第 1項または第 2項に記載の画像処理方法である。

[0029] 請求の範囲第 5項に記載の発明は、前記評価関数の評価特性の最大値を、前記放射線画像に含まれる骨部の画素値の範囲内に有すると共に、前記評価特性の最大値となる画素値の近傍であって該画素値よりも大きい画素値の領域に比較して、前記評価特性の最大値となる画素値の近傍であって該画素値よりも小さい画素値の領域に前記評価特性に緩やかな傾きを有する、ことを特徴とする請求の範囲第 1項または第 2項に記載の画像処理方法である。

[0030] 請求の範囲第 6項に記載の発明は、前記評価関数の評価特性の最大値を、前記放射線画像に含まれる骨部の画素値の範囲内に有し、該評価特性の別の山となる部分を前記骨部の画素値より大きい画素値に有し、前記評価特性の最大値となる画素値の近傍であって該画素値よりも大き!、画素値の領域に比較して、前記評価特性の最大値となる画素値の近傍であって該画素値よりも小さい画素値の領域に前記評価特性に緩やかな傾きを有する、ことを特徴とする請求の範囲第 1項または第 2項に記載の画像処理方法である。

[0031] 請求の範囲第 7項に記載の発明は、前記評価関数の評価特性の最大値を、前記放射線画像に含まれる軟部の画素値の範囲内に有する、ことを特徴とする請求の範囲第 1項または第 2項に記載の画像処理方法である。

[0032] 請求の範囲第 8項に記載の発明は、被写体を透過した放射線の照射線量に応じた信号を有する放射線画像に対し、診断に適した画像を得るための画像処理を行う画像処理装置であって、前記放射線画像に含まれる骨部あるいは軟部を重要領域として検出する重要領域検出手段と、前記重要領域の画素に所定の重みを持たせる重み画像を生成する重み画像生成手段と、前記重み画像の重みと前記放射線画像の画素値とを乗算し、この乗算結果により重み付きヒストグラムを算出するヒストグラム算出手段と、前記重み付きヒストグラムを所定の評価関数により評価し、評価結果が最大となるシフト値を算出する評価手段と、前記シフト値を得る前記評価関数の最大値に対応する前記放射線画像の画素値で所定の処理結果が得られるように画像処理条件を決定する画像処理条件決定手段と、前記画像処理条件に基づ!、て画像処理を実行する画像処理手段と、を有することを特徴とする画像処理装置である。

[0033] 請求の範囲第 9項に記載の発明は、前記画像処理手段は、階調処理を画像処理として実行する、ことを特徴とする請求の範囲第 8項に記載の画像処理装置である。

[0034] 請求の範囲第 10項に記載の発明は、前記評価関数の評価特性の最大値を、前記放射線画像に含まれる骨部の画素値の範囲内に有する、ことを特徴とする請求の範囲第 8項または第 9項に記載の画像処理装置である。

[0035] 請求の範囲第 11項に記載の発明は、前記評価関数の評価特性の最大値を、前記放射線画像に含まれる骨部の画素値の範囲内に有すると共に、該評価特性の別の山となる部分を前記骨部の画素値より大きい画素値に有する、ことを特徴とする請求の範囲第 8項または第 9項に記載の画像処理装置である。

[0036] 請求の範囲第 12項に記載の発明は、前記評価関数の評価特性の最大値を、前記放射線画像に含まれる骨部の画素値の範囲内に有すると共に、前記評価特性の最大値となる画素値の近傍であって該画素値よりも大きい画素値の領域に比較して、前記評価特性の最大値となる画素値の近傍であって該画素値よりも小さい画素値の領域に前記評価特性に緩やかな傾きを有する、ことを特徴とする請求の範囲第 8項または第 9項に記載の画像処理装置である。

[0037] 請求の範囲第 13項に記載の発明は、前記評価関数の評価特性の最大値を、前記放射線画像に含まれる骨部の画素値の範囲内に有し、該評価特性の別の山となる部分を前記骨部の画素値より大きい画素値に有し、前記評価特性の最大値となる画素値の近傍であって該画素値よりも大き!、画素値の領域に比較して、前記評価特性の最大値となる画素値の近傍であって該画素値よりも小さい画素値の領域に前記評価特性に緩やかな傾きを有する、ことを特徴とする請求の範囲第 8項または第 9項に記載の画像処理装置である。 [0038] 請求の範囲第 14項に記載の発明は、前記評価関数の評価特性の最大値を、前記放射線画像に含まれる軟部の画素値の範囲内に有する、ことを特徴とする請求の範囲第 8項または第 9項に記載の画像処理装置である。

発明の効果

[0039] この発明は、被写体を透過した放射線の照射線量に応じた信号を有する放射線画像に対し、診断に適した画像を得るための画像処理を行う際に、放射線画像に含まれる骨部あるいは軟部を重要領域として検出し、この重要領域の画素に所定の重みを持たせる重み画像を生成し、この重み画像の重みと放射線画像の画素値とを乗算して重み付きヒストグラムを算出し、この重み付きヒストグラムを所定の評価関数により評価し、評価結果が最大となるシフト値を算出し、このようなシフト値を得る評価関数の最大値に対応する放射線画像の画素値で所定の処理結果が得られるように画像処理条件を決定し、画像処理条件に基づ!ヽて画像処理を実行する。

[0040] ここで、画像処理は、階調処理である。また、評価関数の評価特性の最大値を、放射線画像に含まれる骨部の画素値の範囲内に有する。

[0041] また、ここで、評価関数の評価特性の最大値を、放射線画像に含まれる骨部の画素値の範囲内に有すると共に、該評価特性の別の山となる部分を骨部の画素値より大きい画素値に有する。または、評価関数の評価特性の最大値を、放射線画像に含まれる骨部の画素値の範囲内に有すると共に、評価特性の最大値となる画素値の近傍であって該画素値よりも大き!/、画素値の領域に比較して、評価特性の最大値となる画素値の近傍であって該画素値よりも小さい画素値の領域に評価特性に緩やかな傾きを有する。または、評価関数の評価特性の最大値を、放射線画像に含まれる骨部の画素値の範囲内に有し、該評価特性の別の山となる部分を骨部の画素値より大き、画素値に有し、評価特性の最大値となる画素値の近傍であって該画素値よりも大き、画素値の領域に比較して、評価特性の最大値となる画素値の近傍であって該画素値よりも小さい画素値の領域に評価特性に緩やかな傾きを有する。

[0042] このように、放射線画像に含まれる重要領域の画素に所定の重みを持たせて重み付きヒストグラムを算出し、この重み付きヒストグラムを所定の評価関数により評価した評価結果が最大となるシフト値を算出して、このようなシフト値を得る評価関数の最大値に対応する放射線画像の画素値で所定の処理結果が得られるように画像処理条件を決定しているため、重要領域で所定の処理結果が得られるようになり、関心領域設定方式などにおける専用の関心領域 (ROI)を設定する必要がなくなり、ロバスト性の高い、適切な条件の画像処理を実行することが可能になる。

図面の簡単な説明

[0043] [図 1]本発明の一実施の形態例の全体構成を機能的に示す機能ブロック図である。

[図 2]本発明の一実施の形態例の全体処理の流れを示すフローチャート図である。

[図 3]本発明の実施の形態例における処理の様子を示す説明図である。

[図 4]本発明の実施の形態例における処理の様子を示す説明図である。

[図 5]本発明の実施の形態例における処理の様子を示す説明図である。

[図 6]本発明の実施の形態例における処理の様子を示す説明図である。

[図 7]本発明の実施の形態例における処理の様子を示す説明図である。

[図 8]本発明の実施の形態例における処理の様子を示す説明図である。

[図 9]本発明の実施の形態例における処理の様子を示す説明図である。

[図 10]本発明の実施の形態例における評価関数の一例を示す図である。

符号の説明

[0044] 5 被写体

30 放射線発生装置

40 放射線画像読取器

100 画像処理装置

101 制御部

110 画像データ生成部

120 重要領域検出部

130 重み画像生成部

140 ヒストグラム算出部

150 評価部

160 画像処理条件決定部

170 画像処理部発明を実施するための最良の形態

[0045] 以下、図面を参照して本発明の実施をするための最良の形態を詳細に説明する。

[0046] 本発明の実施をするための最良の形態の画像処理方法および画像処理装置の好適な実施の形態について説明する。なお、これにより本発明が限定されるものではない。

[0047] なお、本実施例の形態例の各手段は、ハードウェアやファームウェア、またはソフトウェアで構成することが可能である。このため、画像処理方法の各ステップ，画像処理装置の各手段の処理手順に沿った機能ブロック図としての図 1を示す。

[0048] 以下、本実施の形態例の構成および動作について、図 1のブロック図、図 2のフロ一チャート、および、その他の説明図を参照しつつ詳細に説明する。なお、図 1における各手段や各部は、画像形成装置の各手段や各部を示すだけでなぐ画像処理方法の各ステップ、を示すものでもある。

[0049] [第一実施形態]

〈全体構成および処理の流れ〉

システム全体構成：

放射線発生装置 30、放射線画像読取器 40、画像処理装置 100が、図 1に示されるように構成されている。

[0050] また、画像処理装置 100内部には、制御手段を構成する制御部 101、画像データを生成する画像データ生成部 110、画像データに含まれる骨部ある!/、は軟部を重要領域として検出する重要領域検出部 (重要領域検出手段） 120、重要領域の画素に所定の重みを持たせる重み画像を生成する重み画像生成部（重み画像生成手段） 1 30、重み画像の重みと放射線画像の画素値とを乗算し、この乗算結果により重み付きヒストグラムを算出するヒストグラム算出部 (ヒストグラム算出手段） 140、重み付きヒストグラムを所定の評価関数により評価し、評価結果が最大となるシフト値を算出する評価部 (評価手段） 150、以上のシフト値を得る評価関数の最大値に対応する放射線画像の画素値で所定の処理結果が得られるように画像処理条件を決定する画像処理条件決定部（画像処理条件決定手段） 160、画像処理条件に基づいて画像処理を実行する画像処理部（画像処理手段） 170が、図 1に示すように構成される。処理の流れ：

図 1に示した画像処理システムで実行される処理の流れにつ!、て、図 2を用いて説明する。なお、以下の処理の説明では、重要領域として骨部を基準にして処理を進める場合を具体例として説明を行う。

(1)制御部 101は放射線画像撮影'読み取り、および、画像処理に伴う各種処理についての制御を行う（図示せず)。

(2)放射線発生装置 30からの放射線は被写体 5を透過し、この被写体 5を透過した放射線は放射線画像読取器 40により読み取られる。放射線画像読取器 40により読み取られた信号は、画像データ生成部 110で画像データに変換される (ステップ S1)

(3)重要領域検出部 120は、画像データに含まれる骨部エッジにより骨部を重要領域として検出する (ステップ S 2)。

(4)重み画像生成部 130は、骨部が検出された画像データについて二値ィヒ処理を行い、所定の重み、たとえば、骨部の画素で値を 1、骨部以外の画素では値を 0とする重み画像を生成する (ステップ S3)。

(5)ヒストグラム算出部 140は、各画素につ、て重み画像の重みと放射線画像の画素値とを乗算した上で、この乗算結果により重み付きヒストグラムを算出する (ステップ S4)。

(6)評価部 150は、重み付きヒストグラムを所定の評価関数 (特徴量評価関数）により評価して評価結果を求め、この評価結果が最大となるシフト値 (S値)を算出する (ステツプ S5)。

(7)画像処理条件決定部 160は、以上のシフト値を得る評価関数の最大値に対応する画像データの画素値で所定の処理結果が得られるように階調処理の画像処理条件を決定する。ここで、画像処理が階調処理であれば、決定した画像処理条件に応じて階調処理用の LUTを定める（ステップ S6)。なお、画像処理が他の処理であれば、その処理に合わせた画像処理条件を定める。

(8)画像処理部 170は、決定された画像処理条件、階調処理であれば LUTに基づ V、て画像処理を実行する (ステップ S7)。 [0052] 〈各手段、各ステップ、各ルーチンの詳細内容〉

以下、各手段、各ステップ、各ルーチンにおける詳細な処理手順 '内容について説明する。

[0053] 撮影 ·放射線画像入力：

制御部 101では、まず、撮影部位または撮影方向等の情報をユーザインターフエースなど力も取得することとしても良い。これらの情報は、ユーザが撮影部位等を特定することで行う。例えば、表示部とタツチパネルとを兼ね備えた当該画像処理装置 100のユーザインタフェース（図示せず)から、撮影部位が表示されたボタンを押すこと等により入力したり、骨部あるいは軟部を重要領域として検出するかを指定しても良い。磁気カード、バーコード、 HIS (病院内情報システム:ネットワークによる情報管理 )等を利用しても行われる。

[0054] 放射線発生装置 30は、画像処理装置 100内の制御部 101と各種データの授受を行っている。なお、放射線発生装置 30が制御部 101によって制御される場合と、放射線発生装置 30内の制御部（図示せず）と画像処理装置 100内の制御部 101とが通信を行う場合とがある。

[0055] 放射線発生装置 30から放射された放射線は、被写体 5を通して放射線画像読取器 40の前面に装着されている撮像パネルに照射される。そして、放射線画像読取器 40では、被写体 5を透過した放射線を検知して、画像信号として取得する。

[0056] 具体的な構成例としては、輝尽性蛍光体プレートを用いたものとして、特開平 11-1 42998号公報ゃ特開 2002- 156716号公報に記載されたものがある。また、フラットパネルディテクタ（FPD)を入力装置として用いるものには、特開平 6-342098号公報に記載された、検出した X線を直接電荷に変換し、画像信号として取得するものや、特開平 9-90048号公報に記載された、検出した X線をー且光に変換した後、その光を受光して電荷に変換する、間接方式のものがある。

[0057] なお、放射線画像読取器 40は、放射線画像を記録した銀塩フィルムにレーザ、蛍光灯等の光源からの光を照射し、この銀塩フィルムの透過光を光電変換して画像データを生成してもよい。また、放射線量子計数型検出器を用いて放射線エネルギーを直接電気信号に変換して画像データを生成する構成であってもよい。 [0058] 被写体 5の放射線画像を得る際には、放射線発生装置 30と放射線画像読取器 40 の撮像パネルの間に被写体 5が位置するものとされて、放射線発生装置 30から放射された放射線が被写体 5に照射されると共に、被写体 5を透過した放射線が撮像パネルに入射される。

[0059] そして、放射線画像読取器 40により読み取られた画像信号は画像データ生成部 1

10送られ、画像データ生成部 110で所定の画素数や階調ビット数の画像データに変換される。

[0060] (2)重要領域検出：

図 2のステップ S2で実行される重要領域検出処理について説明する。

[0061] 重要領域検出部 120は、重要領域が骨部である場合には、トップハット処理ゃスケルトン処理などの骨格ィ匕を行う各種画像処理によって、画像データに含まれる骨部エッジなどにより骨部を重要領域として検出する。

[0062] この場合、トップハット処理やスケルトン処理のほかに、ラプラシアンフィルタゃソーベルフィルタなどによるエッジ検出手法を用いることも可能である。

[0063] なお、重要領域が軟部である場合には、初めに画像端部にあるエッジを検出することで照射野を認識し、その後に二値化と判別分析とで素抜け領域を抽出し、残りの領域から上述したように抽出した骨部領域を除くような画像処理によって、画像データに含まれる軟部を重要領域として検出する。

[0064] また、複数の手法を用いた場合には、複数の手法による結果を、ある重みに従って結合する。この場合、ニューラルネットワークに従って足し合わせ方を決定する手法も存在している。

[0065] なお、本実施の形態にぉ、て、骨部とは、被写体 5の骨を透過した放射線を検知した画像の領域のことをいい、軟部とは、被写体 5の骨以外の組織や臓器を透過した放射線を検知した画像の領域のことを、う。一般に骨部は放射線の透過量が少なく、軟部は放射線の透過量が多いが、被写体の存在しない領域 (いわゆる素抜け領域 )よりは放射線の透過量が少な、。

[0066] (3)重み画像生成：

図 2のステップ S3で実行される重み画像生成処理について説明する。 [0067] 重み画像生成部 130は、重要領域検出部 120によって骨部が検出された画像データについて二値ィ匕処理を行い、所定の重み、たとえば、骨部の各画素で値を 1、骨部以外の各画素では値を 0とする重み画像を生成する。なお、重要領域が軟部の場合には、軟部の各画素で値を 1、軟部以外の各画素では値を 0とする重み画像を生成する。

[0068] 図 3 (a)は、図 2のステップ S2で入力された元となる放射線画像である。図 3 (b)は、図 2のステップ S3において、重要領域検出部 120によって骨部が検出された画像データについて、重み画像生成部 130が二値ィ匕処理として、骨部の画素で値を 1、骨部以外の画素では値を 0とする重み画像を生成した様子を示している。

[0069] なお、複数の特徴量力も最終的な重み画像を決定するには、ニューラルネットヮークによる手法を用いることが可能である。この場合、教師画像として予め骨部，軟部，素抜け，照射野などに分割した画像を用意しておき、逆誤差伝播法などで最適化を行うことになる。

[0070] (4)ヒストグラム算出：

図 2のステップ S4で実行されるヒストグラム算出処理について説明する。

[0071] ヒストグラム算出部 140は、各画素について、ステップ S3で生成された重み画像の重みと、ステップ S2で読み取られた、重み画像と対応する放射線画像の各画素の画素値とを乗算した乗算結果を得て、この乗算結果の値と出現頻度との重み付きヒストグラムを算出する。なお、この重み付きヒストグラムで、所定の出現頻度以下の部分はノイズなどの可能性があるので除去してもよ!/、。

[0072] 図 4 (a)は、図 3 (b)に示したような骨部を重要領域とした重み画像を用いた場合の重み付きヒストグラムを生成した場合の一例である。ここでは、骨部に相当する画素値の範囲で出現頻度が大きくなつて、る。

[0073] (5)重要度算出：

図 2のステップ S5で実行される重要度算出処理について説明する。

[0074] 評価部 150は、以上の重み付きヒストグラムを所定の評価関数 (特徴量評価関数）により評価して、評価結果及び評価結果が最大となるシフト値 (S値)を算出する。

[0075] ここで、評価関数について説明する。 [0076] 一般に、放射線画像は骨部領域の画像濃度を 1. 0付近となるように出力すると良好な画像となる場合が多い。これを達成するには例えば、画像内より骨部領域を抽出し、骨部領域を構成する画素値の平均を上記の 1. 0付近の濃度で出力するように画像処理するなどの方法が考えられる。

[0077] しかし、画像力骨部領域のみを正確に抽出しようとしても困難なことが多ぐ実際には骨部以外の他の領域もある程度含まれるのが普通である。また、他の領域を除外しようとすれば、画像処理に必要な骨部領域まで除外される可能性も増加し、画像処理結果が所望の処理結果から大きく離れる可能性もある。

[0078] そこで、画像処理のロバスト性を重視するために本実施の形態では、画像処理の基準となる画像の領域に、ある程度骨部領域以外、すなわち重要領域以外の領域を含むことを許容し、骨部領域を重点的に抽出するものである。そして、抽出した画像力も得られたヒストグラムと評価関数を用いて、ロバスト性の高、画像処理条件 (例えば、 LUT (ルックアップテーブル)スライド量）の決定方法を与えるものである。

[0079] 評価関数は、重要領域の抽出方法 (画像内から骨部領域を抽出する方法)や、診断対象とする部位などにより、さまざまな形態が考えられる力例えば、次のようにして与えることができる。

[0080] 図 10は、画素値 Xを横軸としたときの、各種評価関数により算出される評価値を示した図である。シグモイド関数 fl (図 10 (a) )、正規分布 f 2 (図 10 (b) )及び評価関数 f 3 (図 10 (c) )はそれぞれ下記式で表せる。

[0081] [数 1]

1+ e

[0082] [数 2]

f 2

[0083] [数 3]

f 3 = u · e ^{t 2}

[0084] ここで、 a, b, c, s, t, u, kは定数 (パラメータ）である。上記 fl、 f2、 f3を用いて、 [0085] [数 4] f = a1 · f 1 (a2 - f2 + a3 - f3)

[0086] といった評価関数を考える。

[0087] 例えば、骨部領域のみを比較的容易に正確に抽出することが可能な場合は、 a3 = 0とすること〖こより、単峰性の正規分布のみになるので、図 10 (b)に示した評価関数のパラメータを変更した評価関数を得ることが可能である。

[0088] また、 al = a2 = a3 = lとして適したパラメータ、たとえば、下記のようなパラメータを用いることで図 10 (d)に示したような評価関数の評価特性を得ることができる。

[0089] a= 1500

b =400

c = 0. 5

s= 1000

t= 2000

u=0. 9

k= 3000

これらのパラメータの決定は、例えば、ステップ S2で入力された画像を観察しながら、試行錯誤的に決定することができる。例えば、低濃度な (透過放射線量が少ない）画像の比率が高い場合は flの切片 51を調整し、低画素値での評価値を非常に小さい値にすることで階調処理後の画像濃度が低濃度になることを抑制することが可能となるようにパラメータを決定することができる。

[0090] また、パラメータの決定は、予め定められた目標となる画像に近づけるように決定することもできる。目標画像に近づける場合、診断に好適な画像濃度範囲に調整されたサンプル画像を用意し、 1.予め定められた範囲内でパラメータ値を変更し、 2.変更した各パラメータ値における処理済画像とサンプル画像との平均濃度差を算出し、 3 .当該平均濃度差が最小になるようなパラメータ値を最適なものとする方法でもよい。ここで、平均濃度差とはサンプル画像と各パラメータ値で画像処理された処理済画像、それぞれにつ！、画像全体の信号値若しくは信号値に対応した画像濃度の平均値を算出し、それぞれ算出された平均値の差の絶対値を取ったものである。

[0091] ここで、図 4 (b)は中央重点型の評価関数の特性の一例を示している。この中央重点型の評価関数は、中央部分 Cの評価値の最大値が 1. 0であり、左右の周辺にいくにしたがって評価値が徐々に低下する特性になっている。そして、この評価関数は、評価関数の中心値 Bをどの画素値の位置に置くか、左右にシフトさせることが可能に構成されている。

[0092] ここで、評価部 150は、評価関数の中心値をある画素値の位置に置いて、評価関数の評価値と重み付きヒストグラムの出現頻度とを乗じることで評価結果を得る。そして、評価部 150は、評価関数を所定量シフトさせ、評価関数の中心値を別な画素値の位置に置いて、評価関数の評価値と重み付きヒストグラムの出現頻度とを乗じることで評価結果を得る。このようにして、評価部 150は、図 4 (c)の fa、 fb、 fc、 fdに示した如ぐ評価関数の中心値をシフトさせて、複数の評価結果を得る。

[0093] ここで、評価値が最大となるシフトの量を、所定の画素値 (例えば、画素値 0)を基準としたシフト値 (S値)として算出する。なお、このような評価関数による評価の場合、重み付きヒストグラムの重要領域の頻度の分布の形状にもよる力重み付きヒストグラムの出現頻度のピークの部分と評価関数の中央部分とがほぼ一致した状態で評価値が最大となる。すなわち、重み付きヒストグラムの出現頻度のピークの画素値 E力画素値 0からのシフト値にほぼ該当する。

[0094] なお、骨部を評価関数により評価する場合には、濃度 0. 5〜1. 5に相当する画素値の範囲内で評価関数をシフトさせることになる。よって、この範囲を外れた場合には何らかの異常があるとして処理を中断する。

[0095] また、軟部を評価関数により評価する場合には、濃度 1. 5〜2. 5に相当する画素値の範囲内で評価関数をシフトさせることになる。よって、この範囲を外れた場合には何らかの異常があるとして処理を中断する。

[0096] (6)画像処理条件設定：

図 2のステップ S6で実行される画像処理条件の設定処理について説明する。

[0097] 画像処理条件決定部 160は、以上のシフト値を得る評価関数の最大値に対応する画像データの画素値を画像処理することで所定の処理結果が得られるように、階調処理の画像処理条件を決定する。

[0098] 画像処理が階調処理であれば、上述したシフト値に該当する画素値を画像処理することで所望の信号値 (濃度値、輝度値)を得られるような画像処理条件を決定し、その画像処理条件に応じて階調処理用の LUTを定める。

[0099] たとえば、図 4 (d)に示すように、ステップ S5の処理で算出されたシフト値に該当する画素値 (骨部領域の出現頻度ピークの画素値)で濃度 Dl ( = l. 0)となるような画像処理条件を決定し、その画像処理条件に合致する特性の LUTを定める。

[0100] 骨部領域を重要領域として抽出'出力する場合、最大濃度が 3.0で出力されるときにシフト値に該当する画素値では 0.7〜1.3の濃度となるように画像処理するのがよい

。また、ディスプレイなどの輝度表示に関する場合は、シフト値に該当する画素値では最大輝度の 60〜80%で出力するように画像処理するのがよ、。

[0101] また、軟部領域を抽出'出力する場合、最大濃度が 3.0で出力されるときにシフト値に該当する画素値では 1.5〜2.5とするように画像処理するのがよい。また、ディスプレィなどの輝度表示に関する場合は、シフト値に該当する画素値では最大輝度の 15〜

50%で出力するように画像処理するのがよい。

[0102] (7)画像処理：

図 2のステップ S7で実行される画像処理について説明する。

[0103] 画像処理部 170は、画像データ生成部 110からの画像データについて、ステップ S

1乃至ステップ S6の処理によって画像処理条件決定部 160で決定された画像処理条件 (パラメータ）にしたがって、画像処理を実行する。

[0104] 具体的には、例えば階調処理においては、上記の評価関数により決定されたシフト値に従って定められた LUTを使用して、画像データの各画素を LUTによって変換する。このようにすることで、処理済み画像は与えられた評価関数、すなわち上記例では変換後の重要領域が所望の濃度値や輝度値になるように変換された画像となる。また、周波数強調処理ゃィコライゼーシヨン処理においても同様に上記のような評価関数により決定した強調パラメータあるいは補正パラメータにより、処理が行われる。この処理済み画像データは、図示されない画像表示手段、各種外部機器、画像出力装置などに供給される。

[0105] 〈実施形態により得られる効果〉

以上のように、この実施形態では、放射線画像に含まれる重要領域の画素に所定の重みを持たせて重み付きヒストグラムを算出し、この重み付きヒストグラムを所定の評価関数により評価した評価結果が最大となるシフト値を算出して、このようなシフト値を得る評価関数の最大値に対応する放射線画像の画素値で所定の処理結果が得られるように画像処理条件を決定して、るため、重要領域で所定の処理結果が得られるようになり、関心領域設定方式などにおける専用の関心領域 (ROI)を設定する必要がなくなり、ロノスト性の高い、適切な条件の画像処理を実行することが可能になる。

[0106] [第二実施形態]

図 4 (a)に示した重み付きヒストグラムでは、図 3 (b)などの骨部を重要領域とした重み画像について重み付きヒストグラムを生成した場合の一例であつたが、実際には照射野の影響が重み付きヒストグラムに現れることもある。

[0107] 図 5 (a)に示した重み付きヒストグラムは、骨部を重要領域としつつ、照射野の影響の現れた重み画像にっ、て重み付きヒストグラムを生成した場合の一例を示す。ここでは、骨部より小さ、画素値に照射野の画素が大きな頻度で出現して、る。

[0108] このような場合に、図 4 (b)と同じ中央重点型の評価関数によって評価値が最大になるような評価を行うと、図 5 (b)のように、照射野のピーク P1と骨部のピーク P2の中間部分に、両方のピークが評価関数 fの比較的山になっている部分にほぼ均等にか力るような状態になる。

[0109] このため、この状態で評価値が最大となるシフトの量をシフト値 (S値）として算出すると、重み付きヒストグラムの骨部の出現頻度のピーク P2の画素値 E1が、評価関数のシフト値とは一致しなくなり、照射野の出現頻度に影響されて若干小さいシフト値となる (E2)。

[0110] この結果、画像処理条件決定部 160は、以上のシフト値で所定の濃度 Dl (= 1. 0 )となるような画像処理条件を決定してしまい、その画像処理条件に合致する特性の LUTを定める（図 5 (c) )。この場合、シフト値が小さいため LUTの特性が左に寄っており、骨部の画素値は所定の濃度 D1よりも大きな濃度 D1 'に変換されてしまう。すなわち、本来の重要領域で所定の処理結果が得らな、状態になる。

[0111] 一方、図 6 (a)に示した重み付きヒストグラムは、骨部を重要領域としつつ、放射線の素抜けの影響の現れた重み画像について重み付きヒストグラムを生成した場合の一例を示す。ここでは、骨部より大きい画素値に素抜けの画素が大きな頻度で出現している。

[0112] このような場合に、図 4 (b)と同じ中央重点型の評価関数によって評価値が最大になるような評価を行うと、図 6 (b)のように、素抜けのピーク P3と骨部のピーク P4の中間部分に、両方のピークが評価関数 fの比較的山になっている部分にほぼ均等にか力るような状態になる。

[0113] このため、この状態で評価値が最大となるシフトの量をシフト値 (S値）として算出すると、重み付きヒストグラムの骨部の出現頻度のピーク P4の画素値 E3が、評価関数 f のシフト値 E4とは一致しなくなり、素抜けの出現頻度に影響されて大きいシフト値となる（図 6 (b) )。

[0114] この結果、画像処理条件決定部 160は、以上のシフト値で所定の濃度 Dl (= 1. 0 )となるような画像処理条件を決定してしまい、その画像処理条件に合致する特性の LUTを定める（図 6 (c) )。この場合、シフト値が大きいため LUTの特性が右に寄っており、骨部の画素値は所定の濃度 D1よりも小さい濃度 D1"に変換されてしまう。すなわち、本来の重要領域で所定の処理結果が得らな、状態になる。

[0115] そこで、図 4に示した中央重点型の評価関数ではなぐ以下のように変形させた評価関数を用いることが望まし、。

[0116] (_a)骨部と照射野：

評価関数の評価特性の最大値 C2を、放射線画像に含まれる重要領域 (骨部)の画素値の範囲内に有すると共に、評価特性の最大値となる画素値 E5の近傍であつて該画素値よりも大き!、画素値の領域に比較して、評価特性の最大値となる画素値の近傍であって該画素値よりも小さ!、画素値の領域 (照射野による画素値近傍) E7 に評価特性に緩やかな傾きを有する形状の特性にする。

[0117] この評価関数の形状を図 7 (b)に示す。

[0118] このような形状の評価関数 feとすることで、照射野によって評価関数のシフト値 E5 がずれることがなくなる。また、ヒストグラムに照射野が存在しなくても評価関数のシフト値がずれることはな、（図 7 (c) )。

[0119] (b)骨部と素抜け：

評価関数の評価特性の最大値 C3を、放射線画像に含まれる重要領域 (骨部)の画素値の範囲内に有すると共に、該評価特性の別の山 C4となる部分を骨部の画素値より大きい画素値の領域 (素抜けによる画素値近傍) E9に有する形状の特性にする。

[0120] この評価関数の形状を図 8 (b)に示す。

[0121] このような形状の評価関数 fgとすることで、素抜けによって評価関数のシフト値 E8 がずれることがなくなる。また、ヒストグラムに素抜けが存在しなくても評価関数のシフト値がずれることはな、（図 8 (c) )。

[0122] (c)骨部と照射野と素抜け：

評価関数の評価特性の最大値 C5を、放射線画像に含まれる重要領域 (骨部）の画素値の範囲内に有すると共に、評価特性の最大値となる画素値 E 11の近傍であつて該画素値よりも大き!、画素値の領域に比較して、評価特性の最大値となる画素値の近傍であって該画素値よりも小さい画素値の領域 (照射野による画素値近傍) E12 に評価特性に緩やかな傾きを有する形状の特性にし、さらに、該評価特性の別の山 C6となる部分を骨部の画素値より大き、画素値の領域 (素抜けによる画素値近傍) E 13に有する形状の特性にする。

[0123] この評価関数の形状を図 9 (b)に示す。

[0124] このような形状の評価関数 fhとすることで、照射野によって評価関数のシフト値 E11 がずれることがなくなる。また、ヒストグラムに照射野が存在しなくても評価関数のシフト値がずれることはない（図 9 (c) )。そして、このような形状の評価関数とすることで、素抜けによって評価関数のシフト値がずれることがなくなる。また、ヒストグラムに素抜けが存在しなくても評価関数のシフト値がずれることはな、（図 9 (c) )。

[0125] すなわち、重要領域 (骨部)で所定の処理結果が得られるようになり、関心領域設定方式などにおける専用の関心領域 (ROI)を設定する必要がなくなり、ロバスト (堅牢)性の高！ヽ、適切な条件の画像処理を実行することが可能になる。

[0126] (d)骨部と、照射野あるいは素抜け以外：

以上の照射野あるいは素抜け以外の何らかの影響によりヒストグラムに大きな出現頻度の領域が生じる場合には、その画素値に合わせた形状の評価関数とすることで、以上の（a) (b) (c)の場合と同様に良好な結果を得ることが可能になる。

[0127] [その他の実施形態]

以上の実施形態では、骨部を重要領域として具体例を用いて説明してきたが、本実施形態を適用することで、軟部組織など他の重要領域で所望の画像処理を実行することも可會である。

Claims

請求の範囲

[1] 被写体を透過した放射線の照射線量に応じた信号を有する放射線画像に対し、診断に適した画像を得るための画像処理を行う画像処理方法であって、

前記放射線画像に含まれる骨部あるいは軟部を重要領域として検出する重要領域検出ステップと、

前記重要領域の画素に所定の重みを持たせる重み画像を生成する重み画像生成ステップと、

前記重み画像の重みと前記放射線画像の画素値とを乗算し、この乗算結果により重み付きヒストグラムを算出するヒストグラム算出ステップと、

前記重み付きヒストグラムを所定の評価関数により評価し、評価結果が最大となるシフト値を算出する評価ステップと、

前記シフト値を得る前記評価関数の最大値に対応する前記放射線画像の画素値で所定の処理結果が得られるように画像処理条件を決定する画像処理条件決定ステツプと、

前記画像処理条件に基づ！/ヽて画像処理を実行する画像処理ステップと、を有することを特徴とする画像処理方法。

[2] 前記画像処理ステップは、階調処理を画像処理として実行する、ことを特徴とする請求の範囲第 1項に記載の画像処理方法。

[3] 前記評価関数の評価特性の最大値を、前記放射線画像に含まれる骨部の画素値の範囲内に有する、ことを特徴とする請求の範囲第 1項または第 2項に記載の画像処理方法。

[4] 前記評価関数の評価特性の最大値を、前記放射線画像に含まれる骨部の画素値の範囲内に有すると共に、該評価特性の別の山となる部分を前記骨部の画素値より大きい画素値に有する、ことを特徴とする請求の範囲第 1項または第 2項に記載の画像処理方法。

[5] 前記評価関数の評価特性の最大値を、前記放射線画像に含まれる骨部の画素値の範囲内に有すると共に、

前記評価特性の最大値となる画素値の近傍であって該画素値よりも大きい画素値の領域に比較して、前記評価特性の最大値となる画素値の近傍であって該画素値よりも小さヽ画素値の領域に前記評価特性に緩やかな傾きを有する、ことを特徴とする請求の範囲第 1項または第 2項に記載の画像処理方法。

[6] 前記評価関数の評価特性の最大値を、前記放射線画像に含まれる骨部の画素値の範囲内に有し、

該評価特性の別の山となる部分を前記骨部の画素値より大き、画素値に有し、前記評価特性の最大値となる画素値の近傍であって該画素値よりも大きい画素値の領域に比較して、前記評価特性の最大値となる画素値の近傍であって該画素値よりも小さヽ画素値の領域に前記評価特性に緩やかな傾きを有する、ことを特徴とする請求の範囲第 1項または第 2項に記載の画像処理方法。

[7] 前記評価関数の評価特性の最大値を、前記放射線画像に含まれる軟部の画素値の範囲内に有する、ことを特徴とする請求の範囲第 1項または第 2項に記載の画像処理方法。

[8] 被写体を透過した放射線の照射線量に応じた信号を有する放射線画像に対し、診断に適した画像を得るための画像処理を行う画像処理装置であって、

前記放射線画像に含まれる骨部あるいは軟部を重要領域として検出する重要領域検出手段と、

前記重要領域の画素に所定の重みを持たせる重み画像を生成する重み画像生成手段と、

前記重み画像の重みと前記放射線画像の画素値とを乗算し、この乗算結果により重み付きヒストグラムを算出するヒストグラム算出手段と、

前記重み付きヒストグラムを所定の評価関数により評価し、評価結果が最大となるシフト値を算出する評価手段と、

前記シフト値を得る前記評価関数の最大値に対応する前記放射線画像の画素値で所定の処理結果が得られるように画像処理条件を決定する画像処理条件決定手段と、

前記画像処理条件に基づ！/ヽて画像処理を実行する画像処理手段と、を有することを特徴とする画像処理装置。

[9] 前記画像処理手段は、階調処理を画像処理として実行する、ことを特徴とする請求の範囲第 8項に記載の画像処理装置。

[10] 前記評価関数の評価特性の最大値を、前記放射線画像に含まれる骨部の画素値の範囲内に有する、ことを特徴とする請求の範囲第 8項または第 9項に記載の画像処理装置。

[11] 前記評価関数の評価特性の最大値を、前記放射線画像に含まれる骨部の画素値の範囲内に有すると共に、該評価特性の別の山となる部分を前記骨部の画素値より大きい画素値に有する、ことを特徴とする請求の範囲第 8項または第 9項に記載の画像処理装置。

[12] 前記評価関数の評価特性の最大値を、前記放射線画像に含まれる骨部の画素値の範囲内に有すると共に、

前記評価特性の最大値となる画素値の近傍であって該画素値よりも大きい画素値の領域に比較して、前記評価特性の最大値となる画素値の近傍であって該画素値よりも小さヽ画素値の領域に前記評価特性に緩やかな傾きを有する、ことを特徴とする請求の範囲第 8項または第 9項に記載の画像処理装置。

[13] 前記評価関数の評価特性の最大値を、前記放射線画像に含まれる骨部の画素値の範囲内に有し、

該評価特性の別の山となる部分を前記骨部の画素値より大き、画素値に有し、前記評価特性の最大値となる画素値の近傍であって該画素値よりも大きい画素値の領域に比較して、前記評価特性の最大値となる画素値の近傍であって該画素値よりも小さヽ画素値の領域に前記評価特性に緩やかな傾きを有する、ことを特徴とする請求の範囲第 8項または第 9項に記載の画像処理装置。

[14] 前記評価関数の評価特性の最大値を、前記放射線画像に含まれる軟部の画素値の範囲内に有する、ことを特徴とする請求の範囲第 8項または第 9項に記載の画像処理装置。