(57)【要約】【課題】手技用デバイス等の視認性を向上した透過像をリアルタイムに得ることが可能になる放射線画像処理システムを実現する。【解決手段】被写体に放射線を照射して透過像を撮像する放射線画像処理システムであって、第一透過像を生成する画像生成部と、第一透過像の視認性を向上させるための処理を実行する画像処理部と、を備え、画像処理部は、第一透過像から低周波成分を抽出する低周波成分抽出部と、第一透過像を低周波成分で除算する除算処理部と、第一透過像のハローを低減するハロー低減部と、を含む。【選択図】  図６
2020065614.tif 000002
【特許請求の範囲】
【請求項１】
被写体に放射線を照射して透過像を撮像する放射線画像処理システムであって、
第一透過像を生成する画像生成部と、
前記第一透過像の視認性を向上させるための処理を実行する画像処理部と、を備え、
前記画像処理部は、
前記第一透過像から低周波成分を抽出する低周波成分抽出部と、
前記第一透過像を前記低周波成分で除算する除算処理部と、
前記第一透過像のハローを低減するハロー低減部と、を含むことを特徴とする放射線画像処理システム。
【請求項２】
請求項１に記載の放射線画像処理システムであって、
画像を表示する表示部を備え、
前記第一透過像は、複数の画素の輝度値から構成される画像であり、
前記除算処理部は、前記第一透過像を前記低周波成分で除算することによって第二透過像を算出して、前記第二透過像を前記ハロー低減部に出力し、
前記ハロー低減部は、
前記第二透過像の各々の画素の輝度値を第一輝度値で減算して、第三透過像を算出し、
前記第三透過像において、輝度値が正となる画素を抽出することによって第四透過像を算出し、
前記第四透過像の低周波成分を抽出することによって第五透過像を算出し、
前記第五透過像および前記第一透過像を用いて、前記ハローが低減された第六透過像を算出し、
前記表示部に、前記第六透過像を出力することを特徴とする放射線画像処理システム。
【請求項３】
請求項２に記載の放射線画像処理システムであって、
前記ハロー低減部は、
前記第六透過像に対して、任意の特性を有する領域の輝度値を調整するための輝度値変換処理を実行することによって、第七透過像を算出し、
前記表示部に、前記第七透過像を出力することを特徴とする放射線画像処理システム。
【請求項４】
被写体に放射線を照射して透過像を生成する放射線画像処理システムが実行する画像処理方法であって、
前記放射線画像処理システムは、前記透過像を生成する画像生成部と、前記透過像の視認性を向上させるための処理を実行する画像処理部と、を有し、
画像処理方法は、
前記画像生成部が、第一透過像を生成する第１のステップと、
前記画像処理部が、前記第一透過像から低周波成分を抽出し、前記第一透過像を前記低周波成分で除算することによって第二透過像を算出する第２のステップと、
前記画像処理部が、前記第二透過像のハローを低減するための処理を実行する第３のステップと、を含むことを特徴とする画像処理方法。
【請求項５】
請求項４に記載の画像処理方法であって、
前記放射線画像処理システムは、画像を表示する表示部を有し、
前記第一透過像は、複数の画素の輝度値から構成される画像であり、
前記第３のステップは、
前記画像処理部が、前記第二透過像の各々の画素の輝度値を第一輝度値で減算して、第三透過像を算出するステップと、
前記画像処理部が、前記第三透過像において、輝度値が正となる画素を抽出することによって第四透過像を算出するステップと、
前記画像処理部が、前記第四透過像の低周波成分を抽出することによって第五透過像を算出するステップと、
前記画像処理部が、前記第五透過像および前記第一透過像を用いて、前記ハローが低減された第六透過像を算出するステップと、
前記画像処理部が、前記表示部に、前記第六透過像を出力するステップと、を含むことを特徴とする画像処理方法。
【請求項６】
請求項５に記載の画像処理方法であって、
前記第３のステップは、
前記画像処理部が、前記第六透過像に対して、任意の特性を有する領域の輝度値を調整するための輝度値変換処理を実行することによって、第七透過像を算出するステップと、
前記画像処理部が、前記表示部に、前記第七透過像を出力するステップと、を含むことを特徴とする画像処理方法。
【請求項７】
被写体に放射線を照射することによって生成された透過像を処理する計算機に実行させるためのプログラムであって、
第一透過像を受け付ける第１の手順と、
前記第一透過像から低周波成分を抽出し、前記第一透過像を前記低周波成分で除算することによって第二透過像を算出する第２の手順と、
前記第二透過像のハローを低減するための処理を実行する第３の手順と、を前記計算機に実行させることを特徴とするプログラム。
【請求項８】
請求項７に記載のプログラムであって、
前記計算機は、画像を表示する表示装置と接続し、
前記第一透過像は、複数の画素の輝度値から構成される画像であり、
前記第３の手順は、
前記第二透過像の各々の画素の輝度値を第一輝度値で減算して、第三透過像を算出する手順と、
前記第三透過像において、輝度値が正となる画素を抽出することによって第四透過像を算出する手順と、
前記第四透過像の低周波成分を抽出することによって第五透過像を算出する手順と、
前記第五透過像および前記第一透過像を用いて、前記ハローが低減された第六透過像を算出する手順と、
前記表示装置に前記第六透過像を出力する手順と、を含むことを特徴とするプログラム。
【請求項９】
請求項８に記載のプログラムであって、
前記第３の手順は、
前記第六透過像に対して、任意の特性を有する領域の輝度値を調整するための輝度値変換処理を実行することによって、第七透過像を算出する手順と、
前記表示装置に前記第七透過像を出力する手順と、を含むことを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、被写体に放射線を照射して得られた画像（透過像）を処理するシステム、方法、およびプログラムに関する。
【背景技術】
【０００２】
医療現場では、人体（患者）等の被写体に対してＸ線等の放射線を照射し、被写体を透過した放射線の強度分布を検出器で検出して得た画像（以下、透過像と記載する）をリアルタイムで表示しながら、医師や技師により治療や検査等の医療行為が行われている。
【０００３】
医療行為の一例である、内視鏡（胃カメラ）を使って胆管および膵管を造影するＥＲＣＰ（ｅｎｄｏｓｃｏｐｉｃ  ｒｅｔｒｏｇｒａｄｅ  ｃｈｏｌａｎｇｉｏｐａｎｃｒｅａｔｏｇｒａｐｈｙ：内視鏡的逆行性胆道膵管造影）検査では、口から十二指腸まで内視鏡を挿入し、内視鏡の先端から膵管および胆管等の臓器の中にカテーテル（細い管）を挿入し、そのカテーテルから造影剤を注入し、臓器のＸ線写真を撮る。このとき、内視鏡径は太く、臓器に入らないので、通過が容易なカテーテルやガイドワイヤ（金属製の細いワイヤ）等のデバイスを補助的に用いることが一般的である。本明細書では、カテーテルやガイドワイヤ等のデバイスを手技用デバイスと記載する。
【０００４】
前述した透過像のリアルタイムな表示は、前述のような医療行為において、手技用デバイスが正しい位置に挿入されているか否かを確認する場合に用いられる。
【０００５】
ＥＲＣＰ検査でガイドワイヤを膵管に挿入する場合、解剖学的に膵臓は脊椎を横断する位置にあるため、透過像においてガイドワイヤと椎体（脊椎を構成する個々の骨）とが重なりやすい。そのため、Ｘ線は椎体によって大きく減弱されるため、透過像では椎体が映り込んでいる領域が暗くなる。このとき、椎体と重なったガイドワイヤは、視認しにくくなるため、医師によるガイドワイヤの操作に支障をきたすという課題がある。
【０００６】
そこで、透過像に映り込んだ椎体等を減弱することによって、相対的に手技用デバイスの視認性を向上させる画像処理技術の実現が望まれている。
【０００７】
例えば、非特許文献１には、胸部単純Ｘ線画像上に映った肺結節病変の読影を支援するためのソフトウェアの実現を目的として、透過像に映り込んだ肋骨や鎖骨を認識し減弱することで、肺野内の視認性を向上する画像処理技術について述べられている。この技術では、（１）肺野認識処理、（２）骨認識処理、（３）骨信号減弱処理、を行い、骨に起因する信号変化のみを減弱することによって、骨に重なる異常陰影や血管等の微細構造の信号をオリジナル画像のまま残し、病変の視認性を改善している。
【先行技術文献】
【非特許文献】
【０００８】
【非特許文献１】小林剛ほか、「胸部単純Ｘ線ＣＡＤアプリケーション」  Ｂｏｎｅ  Ｓｕｐｐｒｅｓｓｉｏｎ処理の開発、KONICA MINOLTA TECHNOLOGY REPORT VOL.12 pp.71-76 (2015)
【発明の概要】
【発明が解決しようとする課題】
【０００９】
非特許文献１記載の技術では、前述したように、（１）肺野認識処理、（２）骨認識処理、（３）骨信号減弱処理、を行うことが述べられている。非特許文献１によると、これらの具体的な処理内容は、以下のとおりである。
【００１０】
（１）肺野認識処理では、辺縁の性質が異なる四つの領域（肺尖部、外胸郭部、横隔膜部、縦隔部）に肺野を行う分類し、各境界の性質に最適なエッジ情報を用いて領域を抽出している。
【００１１】
（２）骨認識処理では、大量データから構築した骨のモデル情報を用いて、示す事前情報を元にした骨構造の推定値と、対象の胸部単純Ｘ線画像から計測した被写体固有の骨構造の推定結果を合わせることにより、骨の詳細構造を精度よく抽出するロバストな骨認識を実現している。
【００１２】
（３）骨信号減弱処理では、認識された骨候補から骨の信号成分を推定し減弱を行う。推定には肋骨および鎖骨とオーバーラップしていない構造物の信号変化を利用している。
【００１３】
前述したように、非特許文献１記載の技術は、肺結節病変の読影を支援するためのソフトウェアの実現を目的としているため、リアルタイム性は必要とされておらず、高度な画像認識処理等を用いることによって、非リアルタイム処理（オフライン処理）で透過像に映り込んだ骨等を減弱している。
【００１４】
一方、前述したＥＲＣＰ検査等において、リアルタイムで手技用デバイス等の視認性を向上する場合、パーソナルコンピュータ（以下、ＰＣ）等の汎用計算機上で稼働するソフトウェアによって高度な画像認識処理等を行うと多大な処理時間がかかってしまい、手技用デバイスおよび内臓等の被写体の動きが滑らかでない不自然な動きの透過像が表示されてしまう問題がある。そのため、被写体の動きが滑らかで自然な動きの透過像を表示するためには、高速な処理を実行できる特別な計算機リソースを備えた高価な装置が必要となり、不経済である。
【００１５】
本発明は、このような状況に鑑みてなされたものであり、演算量を削減し、手技用デバイスのリアルタイムな視認性を向上できる放射線画像処理システムを提供するものである。
【課題を解決するための手段】
【００１６】
本願において開示される発明のうち、代表的なものの概要を簡単に説明すれば、次のとおりである。すなわち、被写体に放射線を照射して透過像を撮像する放射線画像処理システムであって、第一透過像を生成する画像生成部と、前記第一透過像の視認性を向上させるための処理を実行する画像処理部と、を備え、前記画像処理部は、前記第一透過像から低周波成分を抽出する低周波成分抽出部と、前記第一透過像を前記低周波成分で除算する除算処理部と、前記第一透過像のハローを低減するハロー低減部と、を含む。
【発明の効果】
【００１７】
本発明によれば、手技用デバイス等の視認性を向上した透過像をリアルタイムに得ることが可能になる。上記した以外の課題、構成および効果は、以下の実施例の説明により明らかにされる。
【図面の簡単な説明】
【００１８】
【図１】実施例１の放射線画像処理システムの構成の一例を示すブロック図である。
【図２】実施例１の画像処理部の構成の一例を示すブロック図である。
【図３】実施例１の画像処理部が実行する処理の一例を説明するフローチャートである。
【図４】実施例１の画像処理部が実行する処理の具体例を示す図である。
【図５】実施例１の効果の一例を示す図である。
【図６】実施例２の画像処理部の構成の一例を示すブロック図である。
【図７】実施例２の画像処理部が実行する処理の具体例を示す図である。
【図８】実施例２のハロー低減部の動作の説明図である。
【図９】実施例２のハロー低減部の構成の一例を示すブロック図である。
【図１０】実施例２の画像処理部が実行する処理の一例を説明するフローチャートである。
【発明を実施するための形態】
【００１９】
本発明は、手技用デバイス等の視認性を向上した透過像をリアルタイムに得ることを実現する技術を提供するものである。
【００２０】
以下、添付図面を参照して本発明の実施例について説明する。なお、添付図面は本発明の原理に則った具体的な実施例と実装例を示しているが、これらは本発明の理解のためのものであり、決して本発明を限定的に解釈するために用いられるものではない。
【００２１】
本実施例では、当業者が本発明を実施するのに十分詳細にその説明がなされているが、他の実装および形態も可能で、本発明の技術的思想の範囲と精神を逸脱することなく構成および構造の変更および多様な要素の置き換えが可能であることを理解する必要がある。したがって、以降の記述をこれに限定して解釈してはならない。
【００２２】
また、実施例を説明するための全図において、同一の部材には原則として同一の符号を付し、その繰り返しの説明は省略する。
【実施例１】
【００２３】
＜放射線画像処理システムの構成例＞
図１は、実施例１の放射線画像処理システムの構成の一例を示すブロック図である。
【００２４】
図１に示すように、放射線画像処理システム１０１は、Ｘ線管１０２、高電圧発生部１０３、Ｘ線制御部１０４、絞り１０５、Ｘ線補償フィルタ１０６、絞り／フィルタ制御部１０７、テーブル１０９、機構制御部１１０、Ｘ線検出器１１１、検出器制御部１１２、記憶部１１３、中央処理部１１４、画像処理部１１５、入力部１１６、表示部１１７を有する。
【００２５】
テーブル１０９は、人等の被写体１０８を載せる寝台である。機構制御部１１０は、テーブル１０９と電気的に接続され、被写体１０８が撮影に適した位置となるように、テーブル１０９の動きを制御する。このとき、Ｘ線検出器１１１についても、テーブル１０９と一体的に移動する構造としてもよい。
【００２６】
Ｘ線管１０２は、Ｘ線を発生させ、テーブル１０９の上に配置された被写体１０８に向けて当該Ｘ線を照射する。高電圧発生部１０３は、Ｘ線管１０２と電気的に接続され、Ｘ線管１０２に与える高電圧を発生する。Ｘ線制御部１０４は、高電圧発生部１０３と電気的に接続され、高電圧発生部１０３を制御し、Ｘ線管１０２から照射されるＸ線の線量および線質を制御する。
【００２７】
絞り１０５は、Ｘ線管１０２のＸ線照射方向に配置され、Ｘ線管１０２で発生したＸ線が照射される領域を、Ｘ線吸収率の高い金属の開閉によって制御する。Ｘ線補償フィルタ１０６は、Ｘ線吸収率の高い物質で構成され、被写体１０８のＸ線吸収率の低い部位に到達するＸ線を減衰させることによってハレーションを軽減する。
【００２８】
絞り／フィルタ制御部１０７は、絞り１０５およびＸ線補償フィルタ１０６と電気的に接続され、絞り１０５およびＸ線補償フィルタ１０６を制御する。
【００２９】
Ｘ線検出器１１１は、絞り１０５、Ｘ線補償フィルタ１０６およびテーブル１０９を介在してＸ線管１０２と対向するように配置され、画像生成部として機能する。具体的には、Ｘ線検出器１１１は、Ｘ線管１０２から照射され被写体１０８を透過したＸ線の強度分布を特徴量に変換し、画素毎の特徴量から構成される透過像のデータを出力する。特徴量は、例えば、輝度値及び分散値等である。本明細書では、輝度値を特徴量として有する画像を用いて説明する。
【００３０】
検出器制御部１１２は、Ｘ線検出器１１１と電気的に接続され、Ｘ線検出器１１１を制御することによって透過像のデータを取得し、画像処理部１１５に透過像のデータを入力する。検出器制御部１１２は、Ｘ線検出器１１１を制御することによって、透過像を静止画として生成してもよいし、時間的に異なるタイミングで撮影した複数の透過像を動画像として生成してもよい。動画像を生成するための撮影タイミングは、例えば、毎秒３０フレームおよび毎秒１５フレーム等の一定の時間間隔が考えられる。なお、本発明は、時間間隔に限定されない。
【００３１】
画像処理部１１５は、検出器制御部１１２と電気的に接続され、Ｘ線検出器１１１で撮影され、検出器制御部１１２を介して入力された透過像の補正処理を実行する。
【００３２】
中央処理部１１４は、Ｘ線制御部１０４、絞り／フィルタ制御部１０７、機構制御部１１０、検出器制御部１１２、記憶部１１３、画像処理部１１５、入力部１１６、表示部１１７と電気的に接続され、電気的に接続される各機能部を制御する。中央処理部１１４は、例えば、汎用計算機が有するＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）である。
【００３３】
記憶部１１３は、半導体メモリおよび磁気ディスク等の記録媒体を備え、画像取得条件および画像等をデータとして記憶する。なお、記録媒体の種類は、これに限定されるものではない。
【００３４】
入力部１１６は、使用者が画像取得条件等を設定するためのユーザインターフェースである。入力部１１６として、キーボード、マウス、および制御用ボタン等を有してもよいし、音声入力およびジェスチャー入力など行うためのセンサ等を有してもよい。
【００３５】
表示部１１７は、補正後の画像を表示する。表示部１１７として、ディスプレイおよびプリンタ等を有してもよい。
【００３６】
Ｘ線制御部１０４、絞り／フィルタ制御部１０７、機構制御部１１０、検出器制御部１１２、および画像処理部１１５は専用のハードウェアを用いて実現しているがこれに限定されない。例えば、各ハードウェアをソフトウェアとして実現してもよい。この場合、各ハードウェアの機能を実現するプログラムを記憶部１１３に格納し、中央処理部１１４が当該プログラムにしたがって処理を実行することによって、各ハードウェアの機能を実現する。
【００３７】
画像処理部１１５の詳細については、以下に詳しく述べる。
【００３８】
＜放射線画像処理システムの動作原理＞
まず、放射線（Ｘ線）が被写体を透過したときに得られる各画素の強度分布として出力される透過像Ｉｍ（ｘ，ｙ）、一次線の強度分布として出力される一次線像Ｉｐ（ｘ，ｙ）、および散乱線の強度分布として出力される散乱線像Ｉｓ（ｘ，ｙ）の特性について説明する。なお、以下の式中の（ｘ，ｙ）は画素の位置を示す。
【００３９】
式（１）に示すように、透過像Ｉｍ（ｘ，ｙ）は一次線像Ｉｐ（ｘ，ｙ）および散乱線像Ｉｓ（ｘ，ｙ）の和として表される。ここで、一次線像Ｉｐ（ｘ、ｙ）は式（２）で与えられ、散乱線像Ｉｓ（ｘ，ｙ）は式（３）で与えられる。なお、式（３）の記号は畳み込み演算を表す。
【００４０】
【数１】
2020065614.tif 000003
【００４１】
【数２】
2020065614.tif 000004
【００４２】
【数３】
2020065614.tif 000005
【００４３】
式（２）は、一次線像Ｉｐ（ｘ，ｙ）の強度が被写体の厚みＴ（ｘ，ｙ）に応じて指数的に減衰することを示す。式（２）におけるＩｏ（ｘ，ｙ）は、位置（ｘ，ｙ）における照射線量の強度を示し、μは単位厚みあたりの放射線減衰率（線減弱係数）を示す。式（３）は、散乱線像Ｉｓ（ｘ，ｙ）の強度が一次線像Ｉｐ（ｘ，ｙ）の強度と点拡散関数Ｓσ（Ｔ（ｘ，ｙ））の畳み込み演算で表されることを示す。点拡散関数Ｓσ（Ｔ（ｘ，ｙ））は被写体の厚みＴ（ｘ，ｙ）に応じて変化する関数である。
【００４４】
一般的に、被写体とＸ線検出器１１１の間に、グリッドと呼ばれる放射線の吸収率が高い物質と吸収率が低い物質を薄く交互に積層した板を設置し、Ｘ線検出器１１１に入射する散乱線量を低く抑える運用がなされる。したがって、透過像Ｉｍ（ｘ，ｙ）に含まれる散乱線像Ｉｓ（ｘ，ｙ）の強度は比較的弱いと考えることができる。そこで、簡単のため散乱線像Ｉｓ（ｘ，ｙ）を十分小さいものとして扱い、式（４）に示すように透過像Ｉｍ（ｘ，ｙ）は一次線像Ｉｐ（ｘ，ｙ）とほぼ等しいものと近似して、以下、本実施例の動作原理を説明する。
【００４５】
【数４】
2020065614.tif 000006
【００４６】
以下の説明では、透過像Ｉｍ（ｘ，ｙ）に厚さＴｄ（ｘ，ｙ）の手技用デバイスが映り込んでいる場合を想定する。ただし、透過像Ｉｍ（ｘ，ｙ）の手技用デバイスが存在しない位置ではＴｄ（ｘ，ｙ）＝０とする。このとき、式（４）および式（２）から、透過像Ｉｍ（ｘ，ｙ）は式（５）のように表すことができる。
【００４７】
【数５】
2020065614.tif 000007
【００４８】
また、手技用デバイスが強調表示される透過像Ｉｄ（ｘ，ｙ）は、式（６）式のように表すことができる。
【００４９】
【数６】
2020065614.tif 000008
【００５０】
ここで、手技用デバイスは細い金属で構成されることから、厚みＴｄ（ｘ，ｙ）は位置（ｘ，ｙ）に応じて急峻に変化する。そのため、透過像Ｉｄ（ｘ，ｙ）は低周波成分よりも高周波成分を多く含む。一方、手技用デバイスを除いた被写体（体内の臓器、骨、結果、および筋肉等）の厚みＴ（ｘ，ｙ）は位置（ｘ，ｙ）に応じて比較的緩やかに変化する。そのため、高周波成分よりも低周波成分を多く含む。したがって、手技用デバイスを除いた被写体の透過像Ｉｏ（ｘ，ｙ）ｅｘｐ（−μＴ（ｘ，ｙ））は、透過像Ｉｍ（ｘ，ｙ）の低周波成分像Ｉｍｌ（ｘ，ｙ）で近似することができる。具体的には、式（６）は式（７）のように近似できる。
【００５１】
【数７】
2020065614.tif 000009
【００５２】
以上で述べた放射線画像処理システム１００の動作原理をまとめると、透過像Ｉｄ（ｘ，ｙ）は、透過像Ｉｍ（ｘ，ｙ）をその低周波成分像Ｉｍｌ（ｘ，ｙ）で除算することによって、近似的に算出できる。
【００５３】
＜画像処理部の構成例＞
図２は、実施例１の画像処理部１１５の構成の一例を示すブロック図である。図３は、実施例１の画像処理部１１５が実行する処理の一例を説明するフローチャートである。
【００５４】
図２に示すように、画像処理部１１５は、検出器制御部１１２から、Ｘ線検出器１１１によって撮影された透過像Ｉｍ（ｘ，ｙ）を入力画像として受け付ける（ステップＳ３０１）。透過像Ｉｍ（ｘ，ｙ）は、低周波成分抽出部２０１および除算処理部２０２に入力される。
【００５５】
画像処理部１１５の低周波成分抽出部２０１は、低周波成分像Ｉｍｌ（ｘ，ｙ）を抽出し（ステップＳ３０２）、除算処理部２０２に低周波成分像Ｉｍｌ（ｘ，ｙ）を入力する。
【００５６】
なお、低周波成分抽出部２０１は、水平方向および垂直方向の低周波成分を通過させる２次元ローパスフィルタを用いて実現できる。２次元ローパスフィルタによる低周波成分の抽出方法は公知であるため、詳細な説明を省略する。２次元ローパスフィルタによる低周波成分の抽出処理は演算時間が短く、かつ、演算量が小さい。
【００５７】
画像処理部１１５の除算処理部２０３は、式（７）に示すように、透過像Ｉｍ（ｘ，ｙ）を低周波成分像Ｉｍｌ（ｘ，ｙ）で除算することによって、透過像Ｉｄ（ｘ，ｙ）を算出する（ステップＳ３０３）。
【００５８】
除算処理は画素毎に行われる。なお、除算処理は公知の演算方法であるため詳細な説明は省略する。除算処理は演算時間が短く、かつ、演算量が小さい。
【００５９】
画像処理部１１５の除算処理部２０３は、表示部１１７に透過像Ｉｄ（ｘ，ｙ）を出力し（ステップＳ３０４）、処理を終了する。
【００６０】
＜画像処理部の動作例＞
図４は、実施例１の画像処理部１１５が実行する処理の具体例を示す図である。
【００６１】
図４（ａ）は、透過像の一例であり、人体の胸部を模したファントム（模型）を被写体１０８として撮影することによって取得された透過像の一部である。この画像には、ガイドワイヤ４０２の陰影と、中心線４０１に沿って並ぶ複数の模擬的な椎体４０３の陰影が映っている。なお、中心線４０１は、説明のために追加した線であって、実際の透過像には映っていない。
【００６２】
図４（ｂ）は、図４（ａ）の画像から主要な輪郭線を抽出した画像である。図４（ａ）の画像におけるガイドワイヤ４０２および椎体４０３の位置関係は、図４（ｂ）のようになる。
【００６３】
図４（ｃ）は、図４（ａ）に示した画像の中心線４０１上の輝度値を描いたグラフである。グラフの縦軸は垂直方向の位置を示し、グラフの横軸はその位置における画像の輝度値を示す。なお、グラフに向かって左方向が大きな値となる。図４（ｃ）において、椎体４０３が存在する位置の信号の変化はゆるやかであり、ガイドワイヤ４０２が存在する位置の信号の変化が急峻であることがわかる。これは、ガイドワイヤが細い金属で構成されることから、ガイドワイヤが存在する位置では放射線が大きく減弱されて、ガイドワイヤの周囲よりも輝度値が急峻に小さくなる、すなわち、画像が急峻に暗くなることを示す。
【００６４】
図４（ｄ）は、図４（ｃ）に示す信号から抽出された低周波成分の輝度値を示す。椎体４０３が存在する位置の信号の変化は、図４（ｃ）に示す信号の変化よりも緩やかになる。一方、ガイドワイヤ４０２が存在する位置の信号はほぼ変化していないことがわかる。これは、ガイドワイヤ４０２は細いため、当該部分では高周波成分が主な成分であり、低周波成分抽出部２０１によってその高周波成分が強く減衰されたことを示す。
【００６５】
図４（ｅ）は、図４（ｃ）に示す信号を、図４（ｄ）に示す信号で除算することによって得られた輝度値を示す。ここで、図４（ｃ）に示す信号の輝度値が変化しない場合、図４（ｃ）に示す信号と図４（ｄ）に示した信号の各輝度値は等しくなるため、図４（ｅ）に示す輝度値は１．０となる。また、図４（ｃ）に示す、椎体４０３が存在する位置の輝度値の変化は比較的ゆるやかであるため、図４（ｅ）に示す除算後の椎体４０３が存在する位置の輝度値の変化は比較的小さい。一方、ガイドワイヤ４０２が存在する位置の輝度値の変化は比較的急峻であるため、図４（ｅ）に示す除算後のガイドワイヤ４０２が存在する位置の輝度値は急峻に小さくなる。
【００６６】
このように、透過像Ｉｍ（ｘ，ｙ）をその低周波成分像Ｉｍｌ（ｘ，ｙ）で除算した場合、椎体４０３が存在する位置の信号の変化量（振幅）は緩やかに小さくなるのに対し、ガイドワイヤ４０２が存在する位置の信号の変化量は急峻に小さく。したがって、画像におけるガイドワイヤ４０２の視認性が向上するという効果が得られる。
【００６７】
＜画像処理部における効果例＞
図５は、実施例１の効果の一例を示す図である。以下、この図を用いて、実際の透過像におけるガイドワイヤの視認性向上の効果について説明する。
【００６８】
図５（ａ１）は、人体の胸部を模したファントム（模型）の透過像Ｉｍ（ｘ，ｙ）であり、透過像Ｉｍ（ｘ，ｙ）にはガイドワイヤ５０１および椎体５０２等が映っている。
【００６９】
図５（ａ２）は、図５（ａ１）の矩形部分を拡大した画像を示す。図５（ａ２）に示すように、処理前の透過像Ｉｍ（ｘ，ｙ）では、ガイドワイヤ５０１および椎体５０２が重なっている部分では画像が暗くなっており、ガイドワイヤ５０１が視認しにくい。
【００７０】
図５（ｂ１）は、図５（ａ１）の透過像Ｉｍ（ｘ，ｙ）を処理することによって出力された透過像Ｉｄ（ｘ，ｙ）である。透過像Ｉｄ（ｘ，ｙ）では、ガイドワイヤ５０１が存在する位置が明瞭に表示され、椎体５０２が存在する位置が全体的に灰色となり、輪郭が薄く表示される。
【００７１】
図５（ｂ２）は、図５（ｂ１）の矩形部分を拡大した画像を示す。図５（ｂ２）に示すように、ガイドワイヤ５０１および椎体５０２が重なっている部分でも、ガイドワイヤ５０１が明瞭に視認できる。すなわち、透過像Ｉｄ（ｘ，ｙ）は、もとの透過像Ｉｍ（ｘ，ｙ）より、ガイドワイヤ５０１の視認性が向上している。
【００７２】
＜実施例１のまとめ＞
以上で述べたように、実施例１の放射線画像処理システム１０１は、ガイドワイヤのように金属製で輪郭が急峻な手技用デバイスが存在する位置の輝度が急峻に暗くなる画像を出力できる。これによって、椎体（骨）、内臓、血管等の背景と手技用デバイスとのコントラスト差が付くため、手技用デバイスの視認性が向上する。
【００７３】
画像処理部１１５が実行する処理は、演算時間が短いため、リアルタイムに画像を表示できる。また、画像処理部１１５が有する低周波成分抽出部２０１および除算処理部２０２は、汎用計算機を用いて実現でき、専用な計算機リソースを用意する必要がないため、コストを抑えることができる。
【実施例２】
【００７４】
次に、実施例２について説明する。実施例２では画像処理部１１５の構成および処理が一部異なる。以下、実施例１との差異を中心に実施例２について説明する。
【００７５】
実施例２の放射線画像処理システム１０１の構成は実施例１の放射線画像処理システム１０１と同一の構成である。
【００７６】
＜画像処理部の構成例＞
図６は、実施例２の画像処理部１１５の構成の一例を示すブロック図である。
【００７７】
図６に示すように、画像処理部１１５は、低周波成分抽出部２０１、除算処理部２０２、およびハロー低減部６０１を含む。低周波成分抽出部２０１および除算処理部２０２は、実施例１と同一の構成である。
【００７８】
＜画像処理部の動作例＞
図７は、実施例２の画像処理部１１５が実行する処理の具体例を示す図である。
【００７９】
図７（ａ）は、透過像Ｉｄ（ｘ，ｙ）を示す図である。椎体４０３の輪郭付近、すなわち、図４（ｃ）に示した信号の輝度値が図４（ｄ）に示した信号の輝度値よりも大きくなる位置で、除算後の輝度値が大きくなる現象、すなわち、ハロー７０１が生じる。このような、輪郭付近の輝度値が大きく（明るく）なる現象は、被写体の後ろ側から「後光」が差しているように見えることから「ハロー」と呼ばれている。
【００８０】
手技用デバイスおよび骨のように放射線（Ｘ線）を強く減弱する物体は、透過像では、当該物体の内部は暗く（黒く）映り、周辺は明るく（白く）映る。この透過像の低周波成分は、被写体の輪郭周辺の外側では入力信号よりも若干暗くなる。すなわち、低周波成分像の輝度値が透過像の輝度値よりも小さくなる。したがって、透過像の輝度値を低周波成分像の輝度値で除算すると、１．０よりも大きい値になってハローが発生する。
【００８１】
図７（ａ）に示す透過像Ｉｄ（ｘ，ｙ）において、手技用デバイスの視認性を向上するためには、低周波成分像の通過周波数帯域を狭くし、入力信号の輝度値と低周波成分の輝度値との差を大きくすることが重要である。
【００８２】
しかし、内視鏡のように比較的太く、かつ放射線の減弱量が大きい手技用デバイスが映った透過像に対して、通過周波数帯域の狭い低周波成分抽出を行った場合、ハローが太くかつ明るくなる。手技用デバイスの存在する位置の陰影の周囲に強いハローが生じた場合、不自然な画像になるばかりでなく、医師の誤診および病変部の見逃しに繋がる恐れがある。例えば、ＥＲＣＰでは必ず内視鏡を使用するため、医師が内視鏡の周辺にある臓器を視認する場合、ハローが目障りであり、また、胆管およびハローが重なって胆管を認識できなくなる。
【００８３】
同様に、造影剤を原液のまま使用した場合、造影剤の周辺に太くて明るいハローが生じるため、透過像における胆管壁が不自然に映るため、誤診に繋がる恐れがある。また、ＥＵＳ−ＢＤ（ＥＵＳ−ｇｕｉｄｅｄ  ｂｉｌｉａｒｙ  ｄｒａｉｎａｇｅ：超音波内視鏡ガイド下胆道ドレナージ）では、針（金属）を使用するためハローが発生しやすく、ハローによって針周辺の組織の認識が困難となって、医師の手技に支障をきたす恐れがある。
【００８４】
このように、強いハローの発生によって重大な支障が出てしまう恐れがある。そこで、実施例２ではハロー低減部６０１が、図７（ｂ）に示すようにハロー７０２を低減する処理を実行する。
【００８５】
＜ハロー低減部の動作例＞
図８は、実施例２のハロー低減部６０１の動作の説明図である。
【００８６】
ハロー低減部６０１は、図８（ａ）または図８（ｂ）に示すような入出力特性を有する輝度値変換器として実現することができる。
【００８７】
図８（ａ）は、入力する輝度値が１．０よりも大きい場合、傾きが１．０よりも小さくなるような折れ線状の入出力特性を持った輝度値変換器を示す。当該輝度値変換器は、入力値（ｕ）が１．０以下の場合、出力値（ｖ）はｖ＝ｕの直線８０１となる特性を持ち、入力値（ｕ）が１．０以上の場合、出力値（ｖ）がｖ＝ｋ（ｕ−１）＋１の直線８０２となる特性を持つ。ただし、ｋは０．０より大きく１．０より小さい。当該起動変換器を用いることによって、１．０よりも大きい輝度値を持つハローを目立たなくすることができる。
【００８８】
図８（ｂ）は、曲線８０３の入出力特性を持った輝度値変換を示す。当該輝度値変換器は、（ｕ，ｖ）＝（１．０，１．０）の周辺における曲線の傾きを連続的に変化させることによって出力画像が不自然になることを回避し、かつ、入力値（ｕ）が１．０より大きくなるハローに対応する出力値（ｖ）を小さく抑えることによって、ハローを目立たなくすることができる。また、入力値（ｕ）が１．０よりも小さい場合、出力値（ｖ）の値を入力値（ｕ）の値よりも意図的に小さくする変換を行うことによって、透過像Ｉｄ（ｘ，ｙ）において、希釈した造影剤等の薄い陰影が強調されるため、視認性を向上できる。
【００８９】
なお、曲線８０３の具体的な入出力特性は、椎体４０３等の濃い陰影を弱める効果と、希釈した造影剤等の薄い陰影を強調する効果との兼ね合いに基づいて決定すればよい。例えば、事前の実験等によって、透過像における各物体の視認性を確認しながら特性を決定すればよい。
【００９０】
輝度値の変換は、図８（ａ）および図８（ｂ）に示すような入出力特性を持った、一次式、高次多項式、分数式、指数式、対数式、およびこれらを組み合わせた関数式で実現してもよいし、ルックアップテーブルによって実現できる。ルックアップテーブルの入出力特性は、実験的に決定すればよい。例えば、入力値（ｕ）が１．０よりも大きい場合、ハローの見え方を確認しながら、出力値（ｖ）が入力値（ｕ）より小さくなるように設定し、入力値（ｕ）が１．０よりも小さい場合、椎体４０３等の濃い陰影を弱める効果と、希釈した造影剤等の薄い陰影を強調する効果との兼ね合いに基づいて、出力値（ｖ）を設定する。
【００９１】
＜ハロー低減部の構成例＞
図９は、実施例２のハロー低減部６０１の構成の一例を示すブロック図である。図１０は、実施例２の画像処理部１１５が実行する処理の一例を説明するフローチャートである。
【００９２】
ステップＳ１００１からステップＳ１００３までの処理は、ステップＳ３０１からステップＳ３０３までの処理と同一である。除算処理部２０２が算出した透過像Ｉｄ（ｘ，ｙ）は、ハロー低減部６０１に入力される。
【００９３】
ハロー低減部６０１は、透過像Ｉｄ（ｘ，ｙ）を所定の輝度値で減算する（ステップＳ１００４）。これは、ハローは、透過像Ｉｄ（ｘ，ｙ）において、輝度値が１．０よりも大きくなるという性質を持つためである。
【００９４】
具体的には、減算器９０１が、透過像Ｉｄ（ｘ，ｙ）の各画素の輝度値を所定の輝度値で減算する。減算する輝度値は任意に設定できる。例えば１．０とする。
【００９５】
ハロー低減部６０１は、各画素について、輝度値が正成分であるか否かを判定する（ステップＳ１００５）。
【００９６】
輝度値が正成分であると判定された場合、正成分抽出部９０２が、輝度値の大きさに応じて信号成分を抽出し、低周波成分抽出部９０３が、正成分抽出部９０２によって抽出された信号から低周波成分（低周波成分像）を抽出する（ステップＳ１００６）。
【００９７】
正成分抽出部９０２は、例えば、式（８）にしたがって、入力値（ｕ）から出力値（ｖ）を算出する。式（８）は、入力値（ｕ）および０のうち大きい値を出力する関数を表す。
【００９８】
【数８】
2020065614.tif 000010
【００９９】
医師の手技に支障をきたすような強いハローは低周波成分が主体である。したがって、ハロー低減部６０１は、正成分の低周波成分をハローとして抽出する。
【０１００】
なお、低周波成分抽出部９０３は、低周波成分抽出部２０１と同様に、２次元ローパスフィルタを用いて実現できる。なお、低周波成分抽出部９０３のタップ数（フィルタ長）および通過周波数特性は、抑えたいハローの周波数特性と一致するように決定してもよいし、ハローの見え方を確認しながら実験的に決定してもよい。
【０１０１】
次に、ハロー低減部６０１の減算器９０５は、正成分抽出部９０２が抽出した信号から低周波成分を減算する（ステップＳ１００７）。これによって、ハローを抑えることができる。全ての画素について処理が完了していない場合、ハロー低減部６０１はステップＳ１００５に戻り、同様の処理を実行する。全ての画素について処理が完了している場合、ハロー低減部６０１はステップＳ１００９に進む。
【０１０２】
ステップＳ１００５において、輝度値が負成分であると判定された場合、負成分抽出部９０４が、輝度値の大きさに応じて信号成分を抽出する（ステップＳ１００８）。全ての画素について処理が完了していない場合、ハロー低減部６０１はステップＳ１００５に戻り、同様の処理を実行する。全ての画素について処理が完了している場合、ハロー低減部６０１はステップＳ１００９に進む。
【０１０３】
負成分抽出部９０４は、例えば、式（９）にしたがって、入力値（ｕ）から出力値（ｖ）を算出する。式（９）は、入力値（ｕ）および０のうち小さい値を出力する関数を表す。
【０１０４】
【数９】
2020065614.tif 000011
【０１０５】
ステップＳ１００９では、加算器９０６が、減算器９０５および負成分抽出部９０４の各々の出力値を加算する（ステップＳ１００９）。さらに、加算器９０７は、加算器９０６の出力値に所定の輝度値を加算する（ステップＳ１０１０）。
【０１０６】
次に、ハロー低減部６０１は、ルックアップテーブル９０８を用いて、希釈した造影剤等の薄い陰影の強調等を実現するための輝度値変換を実行する（ステップＳ１０１１）。
【０１０７】
ハロー低減部６０１は、表示部１１７に、処理された透過像Ｉｄ（ｘ，ｙ）を出力し（ステップＳ１０１２）、処理を終了する。
【０１０８】
なお、低周波成分抽出部９０３および減算器９０５は、高周波成分抽出器に置き換えてもよい。また、ステップＳ１０１１の処理は省略されてもよい。
【０１０９】
＜実施例２のまとめ＞
以上で述べたように、実施例２の放射線画像処理システム１０１は、手技用デバイス、造影剤の原液、および椎体のような太い骨などの周囲に発生するハローを抑えることができる。これによって、手技用デバイスの視認性がより向上するため、医師による手技を支障なく円滑に進めることができる。
【０１１０】
また、実施例１と同様に、画像処理部１１５が実行する処理は、演算時間が短いため、リアルタイムに画像を表示できる。また、画像処理部１１５が有する低周波成分抽出部２０１、除算処理部２０２、およびハロー低減部６０１は、汎用計算機を用いて実現でき、専用な計算機リソースを用意する必要がないため、コストを抑えることができる。
【０１１１】
＜全体のまとめ＞
（ｉ）以上で説明した本発明の一形態によれば、手技用デバイス等の視認性が向上した透過像をリアルタイムに得られる。手技用デバイスと、放射線を大きく減弱する骨等の被写体が重なるような手技において、不要なハローを目立たせずに、手技用デバイスの視認性を向上することができる。
【０１１２】
（ｉｉ）本発明を用いて実現した放射線画像処理システム１００の適用先は、ＥＲＣＰおよびＥＵＳ−ＢＤに限定されない。これら以外にも、例えば、ＰＦＣ（ｐａｎｃｒｅａｔｅｉｃ  ａｎｄ  ｐｅｒｉｐａｎｃｒｅａｔｉｃ  ｆｌｕｉｄ  ｃｏｌｌｅｃｉｏｎｓ：膵炎に伴う膵および膵周囲液体貯留に対するドレナージ）、ＰＴＣＤ（ｐｅｒｃｕｔａｎｅｏｕｓ  ｔｒａｎｓｈｅｐａｔｉｃ  ｃｈｏｌａｎｇｉｏｄｒａｉｎａｇｅ：経皮経肝胆道ドレナージ）、シャント造影検査（血管造影検査）等のような幅広い適用先が考えられる。
【０１１３】
（ｉｉｉ）なお、本発明は上記した実施例に限定されるものではなく、様々な変形例が含まれる。また、例えば、上記した実施例は本発明を分かりやすく説明するために構成を詳細に説明したものであり、必ずしも説明した全ての構成を備えるものに限定されるものではない。また、各実施例の構成の一部について、他の構成に追加、削除、置換することが可能である。
【０１１４】
また、上記の各構成、機能、処理部、処理手段等は、それらの一部または全部を、例えば集積回路で設計する等によりハードウェアで実現してもよい。また、本発明は、実施例の機能を実現するソフトウェアのプログラムコードによっても実現できる。この場合、プログラムコードを記録した記憶媒体をコンピュータに提供し、そのコンピュータが備えるプロセッサが記憶媒体に格納されたプログラムコードを読み出す。この場合、記憶媒体から読み出されたプログラムコード自体が前述した実施例の機能を実現することになり、そのプログラムコード自体、およびそれを記憶した記憶媒体は本発明を構成することになる。このようなプログラムコードを供給するための記憶媒体としては、例えば、フレキシブルディスク、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、ハードディスク、ＳＳＤ（Ｓｏｌｉｄ  Ｓｔａｔｅ  Ｄｒｉｖｅ）、光ディスク、光磁気ディスク、ＣＤ−Ｒ、磁気テープ、不揮発性のメモリカード、ＲＯＭなどが用いられる。
【０１１５】
また、本実施例に記載の機能を実現するプログラムコードは、例えば、アセンブラ、Ｃ／Ｃ＋＋、ｐｅｒｌ、Ｓｈｅｌｌ、ＰＨＰ、Ｐｙｔｈｏｎ、Ｊａｖａ（登録商標）等の広範囲のプログラムまたはスクリプト言語で実装できる。
【０１１６】
さらに、実施例の機能を実現するソフトウェアのプログラムコードを、ネットワークを介して配信することによって、それをコンピュータのハードディスクやメモリ等の記憶手段またはＣＤ−ＲＷ、ＣＤ−Ｒ等の記憶媒体に格納し、コンピュータが備えるプロセッサが当該記憶手段や当該記憶媒体に格納されたプログラムコードを読み出して実行するようにしてもよい。
【０１１７】
上述の実施例において、制御線や情報線は、説明上必要と考えられるものを示しており、製品上必ずしも全ての制御線や情報線を示しているとは限らない。全ての構成が相互に接続されていてもよい。
【符号の説明】
【０１１８】
１０１  放射線画像処理システム
１０２  Ｘ線管
１０３  高電圧発生部
１０４  Ｘ線制御部
１０５  絞り
１０６  Ｘ線補償フィルタ
１０７  絞り／フィルタ制御部
１０８  被写体
１０９  テーブル
１１０  機構制御部
１１１  Ｘ線検出器
１１２  検出器制御部
１１３  記憶部
１１４  中央処理部
１１５  画像処理部
１１６  入力部
１１７  表示部
２０１、９０３  低周波成分抽出部
２０２  除算処理部
２０３  除算処理部
６０１  ハロー低減部
９０１、９０５  減算器
９０２  正成分抽出部
９０４  負成分抽出部
９０６、９０７  加算器
９０８  ルックアップテーブル
【図１】
2020065614.tif 000012
【図２】
2020065614.tif 000013
【図３】
2020065614.tif 000014
【図４】
2020065614.tif 000015
【図５】
2020065614.tif 000016
【図６】
2020065614.tif 000017
【図７】
2020065614.tif 000018
【図８】
2020065614.tif 000019
【図９】
2020065614.tif 000020
【図１０】
2020065614.tif 000021