(57)【要約】      （修正有）【課題】医用画像のノイズ成分を低減する処理を行う画像処理装置において、画像のエッジが不鮮明になるのを防ぎつつ、エッジのない部分のノイズ成分を十分に低減できるようにする。【解決手段】画像処理装置は、画像にエッジの延長方向に沿った平滑化処理を施して第１平滑化画像を生成する第１平滑化処理と、画像又は第１平滑化画像の全体に平滑化処理を施して第２平滑化画像を生成する第２平滑化手段と、画像及び第１平滑化画像から第１特徴量を抽出する第１特徴量抽出手段と、第１平滑化画像及び第２平滑化画像から第２特徴量を抽出する第２特徴量抽出手段と、第１特徴量及び第２特徴量から所定の指標値を算出する指標値算出手段と、指標値に基づいて、画像、第１平滑化画像及び第２平滑化画像のうちの少なくともいずれかの画像の画素の信号値を用いた新たな信号値を算出し、新たな信号値の画素からなる合成画像を生成する合成画像生成手段と、を備える。【選択図】図２
2018151715.tif 000002
【特許請求の範囲】
【請求項１】
画像に、当該画像に映る物体のエッジの延長方向に沿った平滑化処理を施して第１平滑化画像を生成する第１平滑化手段と、
前記画像又は前記第１平滑化画像の全体に平滑化処理を施して第２平滑化画像を生成する第２平滑化手段と、
前記画像の各画素の信号値及び前記第１平滑化画像における対応する画素の信号値から第１特徴量を抽出する第１特徴量抽出手段と、
前記第１平滑化画像の各画素の信号値及び前記第２平滑化画像における対応する画素の信号値から第２特徴量を抽出する第２特徴量抽出手段と、
抽出された前記第１特徴量及び前記第２特徴量から、所定の指標値を画素ごとに算出する指標値算出手段と、
算出された指標値に基づいて、前記画像、前記第１平滑化画像及び前記第２平滑化画像のうちの少なくともいずれかの画像の画素の信号値を用いた新たな信号値を画素ごとに算出し、算出した前記新たな信号値の画素からなる合成画像を生成する合成画像生成手段と、を備えることを特徴とする画像処理装置。
【請求項２】
前記合成画像生成手段は、
前記決定された指標値に基づいて、前記画像、前記第１平滑化画像及び前記第２平滑化画像のうちの少なくとも二つの画像の画素の信号値の加重平均を前記新たな信号値とする重み付け手段を備えることを特徴とする請求項１に記載の画像処理装置。
【請求項３】
前記第１平滑化手段は、画素ごとにエッジの延長方向を判定し、判定した延長方向に沿って画像を平滑化することを特徴とする請求項１又は２に記載の画像処理装置。
【請求項４】
前記第１特徴量抽出手段は、前記画像における各画素の信号値と前記第１平滑化画像における対応する画素の信号値との差分を前記第１特徴量として抽出することを特徴とする請求項１から３のいずれか一項に記載の画像処理装置。
【請求項５】
前記第２特徴量抽出手段は、前記第１平滑化画像における各画素の信号値と前記第２平滑化画像における対応する画素の信号値との差分を前記第２特徴量として抽出することを特徴とする請求項１から４のいずれか一項に記載の画像処理装置。
【請求項６】
前記指標値算出手段は、前記第１特徴量と前記第２特徴量の比の値を前記指標値として算出することを特徴とする請求項１から５のいずれか一項に記載の画像処理装置。
【請求項７】
前記画像は、放射線画像撮影装置を用いて撮影された放射線画像であることを特徴とする請求項１から６のいずれか一項に記載の画像処理装置。
【請求項８】
前記画像は、放射線画像撮影装置が出力したローデータに基づく原画像、前記原画像から一部の周波数成分が除かれた帯域制限画像又は前記原画像を構成する複数の画素のうちの一部からなる低解像度近似画像であることを特徴とする請求項７に記載の画像処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、医療用のモダリティーで撮影された医用画像を処理する画像処理装置に関し、特に、画像データに混入したノイズを低減する処理を行う画像処理装置に関する。
【背景技術】
【０００２】
医療用のモダリティー（例えば、ＣＲ（Computed Radiography）やＣＴ（Computed Tomography）、ＭＲＩ（Magnetic Resonance Imaging）等）を用いて医用画像を撮影する際、電磁波等の影響により、画像データにノイズが混入し、表示される医用画像が不明瞭になってしまうことがある。
こうしたノイズ成分を除去するための技術は、これまで多数提案されてきているが、従来のノイズ除去技術では、場合によって画像のエッジが不鮮明になってしまうことがあった。こうした事態を防ぐため、近年では、画素ごとに画素ベクトルを算出し、算出した画素ベクトルに応じてノイズ成分とエッジ成分とを分離し、ノイズ成分に平滑化処理を施すとともにエッジ成分に強調処理を施す画像処理方法が提案されている（特許文献１参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０００−３０６０８９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、特許文献１に記載された画像処理方法では、エッジの方向に沿って平滑化を行うため、エッジのない平坦な部分のノイズ成分が十分に低減できない場合がある。
本発明は、このような問題点に鑑みてなされたものであり、医用画像のノイズ成分を低減する処理を行う画像処理装置において、画像のエッジが不鮮明になるのを防ぎつつ、エッジのない部分のノイズ成分を十分に低減できるようにすることを課題とする。
【課題を解決するための手段】
【０００５】
上記課題を解決するため、本発明は、
画像にエッジの延長方向に沿った平滑化処理を施して第１平滑化画像を生成する第１平滑化手段と、
前記画像又は前記第１平滑化画像にエッジをぼかす平滑化処理を施して第２平滑化画像を生成する第２平滑化手段と、
前記画像の各画素の信号値及び前記第１平滑化画像における対応する画素の信号値から第１特徴量を抽出する第１特徴量抽出手段と、
前記第１平滑化画像の各画素の信号値及び前記第２平滑化画像における対応する画素の信号値から第２特徴量を抽出する第２特徴量抽出手段と、
抽出された前記第１特徴量及び前記第２特徴量から、所定の指標値を画素ごとに算出する指標値算出手段と、
算出された指標値に基づいて、前記画像、前記第１平滑化画像及び前記第２平滑化画像のうちの少なくともいずれかの画像の画素の信号値を用いた新たな信号値を画素ごとに算出し、算出した前記新たな信号値の画素からなる合成画像を生成する合成画像生成手段と、を備えることを特徴とする。
【発明の効果】
【０００６】
本発明によれば、画像のエッジが不鮮明になるのを防ぎつつ、エッジのない部分のノイズ成分を十分に低減することができる。
【図面の簡単な説明】
【０００７】
【図１】本発明の第１（第２）実施形態に係る画像処理装置の機能的構成の一例を示すブロック図である。
【図２】図１の画像処理装置が実行するノイズ低減処理のフローチャートである。
【図３】図２ノイズ低減処理における第１平滑化処理で用いるフィルターの例を示した模式図である。
【図４】第１平滑化処理におけるフィルター選択の概念図である。
【図５】図１の画像処理装置が実行するノイズ低減処理の概念図である。
【図６】本発明の第２実施形態に係る画像処理装置が実行するノイズ低減処理のフローチャートである。
【発明を実施するための形態】
【０００８】
＜第１実施形態＞
以下、図を参照して本発明の第１実施形態について詳細に説明する。ただし、発明の範囲は、図示例に限定されるものではない。
【０００９】
〔画像処理装置の構成〕
まず、本実施形態における画像処理装置１の構成について説明する。図１は、本実施の形態における画像処理装置１の機能的構成の一例を示すブロック図である。
【００１０】
本実施形態の画像処理装置１は、医用画像（特に放射線画像）に混入したノイズ成分を低減するためのものである。画像処理装置１は、図示しない医療用のモダリティー（例えば、ＣＲ（Computed Radiography）やＣＴ（Computed Tomography）、ＭＲＩ（Magnetic Resonance Imaging）等）と直接もしくはサーバー（例えば、医療用画像管理システム（Picture Archiving and Communication System：ＰＡＣＳ）を構成するもの）を介して接続され、モダリティー又はサーバーとの間で医用画像の画像データの送受信が可能となっている。
【００１１】
画像処理装置１は、図１に示すように、制御部１１やＲＡＭ１２、通信部１３、記憶部１４、操作部１５、表示部１６等を備えて構成されており、各部１１〜２１はバス１７によって接続されている。
【００１２】
制御部１１は、ＣＰＵ（Central Processing Unit）等から構成されており、画像処理装置１の各部の処理動作を統括的に制御する。具体的には、記憶部１４に記憶されている各種処理プログラムを読み出してＲＡＭ１２に展開し、当該プログラムとの協働により、後述するノイズ低減処理等の各種処理を行う。
【００１３】
ＲＡＭ１２は、制御部１１により実行制御される各種処理において、記憶部１４から読み出された各種プログラム、入力若しくは出力データ及びパラメーター等を一時的に記憶するワークエリアを形成する。
【００１４】
通信部１３は、ネットワークインターフェース等により構成され、ＬＡＮ（Local Area Network）、ＷＡＮ（Wide Area Network）、インターネット等の通信ネットワークを介して接続されたモダリティー又はサーバーとの間で医用画像の画像データ等の送受信を行うことが可能となっている。
【００１５】
記憶部１４は、ＨＤＤ（Hard Disk Drive）や半導体の不揮発性メモリー等により構成
され、各種処理プログラム、当等を記憶することが可能となっている。
また、記憶部１４は、モダリティー等から送信されてきた医用画像の画像データを管理するためのデータベースを有している。画像データには、例えば、患者名、検査種類、日付等の所定の付帯情報が付帯され、各画像データの識別を行うことが可能となっている。
また、記憶部１４は、各種画像処理を施した画像データや、画像処理を行うためのツール（フィルター等）を記憶することも可能となっている。
【００１６】
操作部１５は、処理する対象となる画像データを指定したり、処理の開始を指示したりするためのもので、キーボードやマウス、表示部１６に設けられたタッチパネル等で構成されている。
【００１７】
表示部１６は、画像データに基づく医用画像を表示するためのもので、液晶ディスプレイ（Liquid Crystal Display、ＬＣＤ）やエレクトロルルミネッセンス（Electroluminescence：ＥＬ）を利用した表示装置等で構成されている。
【００１８】
（画像処理装置１の動作）
次に、上記画像処理装置１の動作について説明する。図２は画像処理装置１が実行するノイズ低減処理のフローチャート、図３はノイズ低減処理における第１平滑化処理で用いるフィルターの例を示した模式図、図４は第１平滑化処理におけるフィルター選択の概念図、図５は画像処理装置１が実行するノイズ低減処理の概念図である。
【００１９】
画像処理装置１の制御部１１は、操作部１５に対し処理対象の画像データの選択を指示する操作及び当該画像データに対する画像処理の開始を指示する操作がなされると、図２に示したノイズ低減処理を実行する。
ノイズ低減処理では、まず、処理対象となる画像Ｉ０の画像データを準備する（ステップＳ１）。具体的には、モダリティー又はサーバーから画像データを受信する又は記憶部１４に記憶されている画像データを読み込む。
【００２０】
なお、画像Ｉ０には、モダリティーが出力したローデータに基づく原画像や、原画像に多重解像度解析（例えばウェーブレット変換等）を施すことにより得られた帯域の異なる複数の帯域制限画像のうちのいずれか（すなわち、一部の周波数成分が除かれた画像）、原画像から所定間隔ごとに抽出した一部の画素で構成した低解像度近似画像等が含まれる。
画像Ｉ０として帯域制限画像を用いる場合には、中周波帯域のものを用いるのが好ましい。
【００２１】
画像データの準備を終えた後は、画像データに第１平滑化処理を施す（ステップＳ２）。具体的には、画像Ｉ０に、当該画像Ｉ０に映る物体（患者の撮影対象部位）のエッジ（輪郭や表面の急峻な構造変化）の延長方向に沿った平滑化処理（選択的平滑化処理）を施し、選択的平滑化画像Ｉ１を生成する。すなわち、制御部１１は、本発明における第１平滑化手段として機能するものである。
また、制御部１１は、選択的平滑化画像Ｉ１を生成したら、その画像データを記憶部１４に記憶する。
【００２２】
本実施形態においては、構造パターンが異なる複数のフィルターを順次当てはめてその結果を評価することで、エッジの延長方向（平滑化する方向）を絞り込むようにしている。具体的には、まず、図３に示したような複数種類のフィルターを用意しておく。各フィルターは、縦横とも複数の画素Ｐ（例えば縦７×横７ピクセル）に跨るようになっており、破線に沿って構造を近似するようになっている。例えば、図３のＦ１で示されたフィルターを画素Ｐの上に当てはめると、図４（ａ）に示したように、中心の画素を含む横一列の画素（太線で囲まれた領域）からなる領域の構造（複数の信号値の傾向）を直線Ｌで近似する。また、図３のＦ２で示されたフィルターを画素の上に当てはめると、図４（ｂ）に示したように、中心部の縦３×横３ピクセルの画素（太線で囲まれた領域）からなる領域の構造を二次曲面Ｃで近似する。このようなフィルターを用いるのは、構造の多くは、局所的にみると、直線、各、円のいずれかでできているとみなすことができるためである。
【００２３】
このような構造の近似を、一の画素に対し複数のフィルターを順次当てはめて行い、フィルターごとに実際の構造と直線又は曲面との誤差の大きさを評価する。なお、同じ画素に対し、直線近似用のフィルターＦ１と曲面近似用のフィルターＦ２のいずれを当てはめても、測定する誤差は、直線又は曲面からの離れ方の程度であるため、同じ基準で誤差を評価することができる。そして、当てはめた複数のフィルターの中から誤差が最も小さかったフィルターを用いて平滑化を行う。
このように、構造をフィルターの範囲内だけで評価するため、フィルターの範囲外の構造やノイズ成分の影響を受けることなく、適切なフィルターの選択を行うことができる。
【００２４】
第１平滑化処理を終えたら、画像Ｉ０又は選択的平滑化画像Ｉ１に第２平滑化処理を施す（ステップＳ３）。具体的には、画像Ｉ０又は選択的平滑化画像Ｉ１に対し、画像全体に平滑化を施すことで、ノイズ成分は十分に除去されているがエッジが不鮮明な非鮮鋭画像を生成する。すなわち、制御部１１は、本発明における第２平滑化手段として機能するものである。平滑化には、従来画像の平滑化に用いられてきた各種フィルター（例えば移動平均フィルターやガウシアンフィルター等）を用いることで施すことが可能である。
また、制御部１１は、非鮮鋭画像Ｉ２を生成したら、その画像データを記憶部１４に記憶する。
なお、第２平滑化処理は、第１平滑化処理と並行して実行してもよい。また、画像Ｉ０に対して第２平滑化処理を施す場合には、第１平滑化処理に先行して実行することも可能である。
【００２５】
第１，第２平滑化処理を終えた後は、画素ごとに第１差分処理を行う（ステップＳ４）。具体的には、画像Ｉ０における一の画素の信号値及び選択的平滑化画像Ｉ１における対応する画素の信号値から、第１特徴量Ｎを抽出する。すなわち、制御部１１は、本発明における第１特徴量抽出手段として機能するものである。
この動作は、画像Ｉ０や選択的平滑化画像Ｉ１を構成する各画素について繰り返し行う。
【００２６】
本実施形態においては、図５に示したように、画像Ｉ０における一の画素の信号値と選択的平滑化画像Ｉ１における対応する画素の信号値との差分の信号値を第１特徴量Ｎとしている。画像Ｉ０にはノイズ成分とエッジ成分の両方が含まれており、選択的平滑化画像Ｉ１には多くの構造成分と若干のノイズ成分が含まれているので、両者の差分、すなわち第１特徴量Ｎは、ノイズ成分を主とする信号値となる。
各画素の第１特徴量が得られると、図５に示したようなノイズ成分画像ＩＮを表示部１６に表示することができる。
【００２７】
第１差分処理を終えた後は、画素ごとに第２差分処理を行う（ステップＳ５）。具体的には、選択的平滑化画像Ｉ１における一の画素の信号値及び非鮮鋭画像Ｉ２における対応する画素の信号値から、第２特徴量Ｃを抽出する。すなわち、制御部１１は、本発明における第２特徴量抽出手段として機能するものである。
この動作は、選択的平滑化画像Ｉ１や非鮮鋭画像Ｉ２を構成する各画素について繰り返し行う。
なお、第２特徴量Ｃの抽出は、第１特徴量の抽出と並行して実行してもよい。
【００２８】
本実施形態においては、図５に示したように、選択的平滑化画像Ｉ１における一の画素の信号値と非鮮鋭画像Ｉ２における対応する画素の信号値との差分の信号値を第２特徴量Ｃとしている。選択的平滑化画像Ｉ１には多くの構造成分と若干のノイズ成分とが含まれており、非鮮鋭画像Ｉ２には構造成分及びノイズ成分が少しずつ含まれているので、両者の差分、すなわち第２特徴量Ｃは、構造（コントラスト）成分を主とする信号値となる。
各画素の第２特徴量が得られると、図５に示したような構造成分画像ＩＣを表示部１６に表示することができる。
なお、画像Ｉ０の画素の信号値と、非鮮鋭画像Ｉ２の画素の信号値との差分の信号値を第２特徴量Ｃとしてもよい。
【００２９】
図２の説明に戻る。第１，第２差分処理を終えた後は、得られた第１特徴量Ｎ及び第２特徴量Ｃから重み付け指標値Ｔを画素ごとに算出する（ステップＳ６）。すなわち、制御部１１は、本発明における指標値算出手段として機能するものである。重み付け指標値Ｔとは、対象の画素には、ノイズ成分と構造成分のどちらが相対的に多いかを示す値で、例えば、第１特徴量Ｎに対する第２特徴量Ｃの比の値（Ｃ／Ｎ）や、対象の画素を中心とする所定範囲（例えば縦９×横９ピクセル）における各画素の第１特徴量Ｎの二乗の合計に対する、同範囲における各画素の第２特徴量Ｃの二乗の合計の比の値の平方根（sqrt(sum(C^2)/sum(N^2))等を用いて算出する。
重み付け指標値Ｔを算出したら、重み付け指標値Ｔと予め記憶部に記憶されている所定の閾値Ｔ_thとを比較する（ステップＳ７）。閾値Ｔ_thは、各画素に共通の定数（例えば１．５）である。
【００３０】
ステップＳ７において、重み付け指標値Ｔの値が閾値Ｔ_thよりも大きい場合、すなわち、対象の画素が、ノイズ成分の少ないあるいは構造成分の多い（又はその両方の）画素である場合には（ステップＳ７；Ｙｅｓ）、重み付け指標値Ｔに基づいて、画像Ｉ０及び選択的平滑化画像Ｉ１のうちの少なくともいずれかの画像の画素の信号値を用いた新たな信号値を算出する（ステップＳ８）。具体的には、画像Ｉ０における対応する画素の信号値と同じ信号値にする、あるいは画像Ｉ０と選択的平滑化画像Ｉ１との間でより画像Ｉ０の信号値に近い信号値を設定する、といった方法が挙げられる。
【００３１】
画像Ｉ０の画素の信号値に近い新たな信号値は、例えば、画像Ｉ０と選択的平滑化画像Ｉ１がＴ−Ｔ_th：Ｔ_thの比率となるように、すなわち、下記（１）で示した式を用いて原画像Ｉ０の画素の信号値と選択的平滑化画像Ｉ１の画素の信号値との加重平均を算出することで得ることができる。
新たな信号値＝（画像Ｉ０の画素の信号値×（Ｔ−Ｔ_th）＋選択的平滑化画像Ｉ１の画素の信号値×Ｔ_th）／Ｔ・・（１）
このように、エッジ成分が相対的に多いと判断された画素については、ノイズ成分の低減量が少ない画像に基づく信号値とすることで、エッジが不鮮明になるのを防ぐことができる。
【００３２】
一方、ステップＳ７において、重み付け指標値Ｔが閾値Ｔ_thよりも小さい場合、すなわち、対象の画素が、ノイズ成分が多いあるいはエッジ成分が少ない（又はその両方の）画素である場合には（ステップＳ７；Ｎｏ）、重み付け指標値Ｔに基づいて、非鮮鋭画像Ｉ２及び選択的平滑化画像Ｉ１のうちの少なくともいずれかの画像の画素の信号値を用いた新たな信号値を算出する（ステップＳ９）。具体的には、非鮮鋭画像Ｉ２における対応する画素の信号値と同じ値にする、あるいは非鮮鋭画像Ｉ２と選択的平滑化画像Ｉ１との間でより非鮮鋭画像Ｉ２の信号値に近い信号値を設定する、といった方法が挙げられる。
【００３３】
非鮮鋭画像Ｉ２の画素の信号値に近い新たな信号値は、例えば、非鮮鋭画像Ｉ２と選択的平滑化画像Ｉ１がＴ_th−Ｔ：Ｔの比率となるように、すなわち、下記（２）で示した式を用いて非鮮鋭画像Ｉ２の画素の信号値と選択的平滑化画像Ｉ１の画素の信号値との加重平均を算出することで得ることができる。
新たな信号値＝（非鮮鋭画像Ｉ２の画素の信号値×（Ｔ−Ｔ_th）＋選択的平滑化画像Ｉ１の画素の信号値×Ｔ）／（Ｔ_th）・・（２）
このように、構造成分が相対的に多いと判断された画素については、ノイズ成分の低減量が多い画像に基づく信号値とすることで、ノイズ成分を大幅に低下させることができる。
【００３４】
ステップＳ８又はステップＳ９の処理の後は、直前に行った処理が対象の画素のうち最後の画素について行ったものであるか否か、すなわち、対象となる全ての画素についてステップＳ８又はステップＳ９の処理を行ったか否かを判定する（ステップＳ１０）。
ここで、最後の画素について行った処理ではないと判定した場合には（ステップＳ１０；Ｎｏ）、画素を特定するカウント番号を＋１更新する（処理対象の画素を一つ隣に移す）処理を行って（ステップＳ１１）、ステップＳ７の処理に戻る。
【００３５】
一方、ステップＳ１０の処理において、最後の画素について行ったものであると判定した場合には（ステップＳ１０；Ｙｅｓ）、得られた新たな信号値の画素からなる合成画像を生成する（ステップＳ１２）。すなわち、制御部１１は、本発明における合成画像生成手段をなす。上述したように、構造成分の多い画素については、画像Ｉ０又は選択的平滑化画像Ｉ１に近い信号値とし、構造成分の少ない画素については、非鮮鋭画像Ｉ２又は選択的平滑化画像Ｉ１に近い信号値としているので、合成画像は、エッジがはっきりとしつつも、ノイズ成分が十分に低減されたものとなる。
【００３６】
このように、本実施形態に係る画像処理装置１は、画像Ｉ０に、当該画像に映る物体のエッジの延長方向に沿った平滑化処理を施して選択的平滑化画像Ｉ１（第１平滑化画像）を生成する第１平滑化手段と、画像Ｉ０又は選択的平滑化画像Ｉ１の全体にエッジをぼかす平滑化処理を施して非鮮鋭画像Ｉ２（第２平滑化画像）を生成する第２平滑化手段と、画像Ｉ０の各画素の信号値及び選択的平滑化画像Ｉ１における対応する画素の信号値から第１特徴量（両者の差分）を抽出する第１特徴量抽出手段と、選択的平滑化画像Ｉ１の各画素の信号値及び非鮮鋭画像Ｉ２における対応する画素の信号値から第２特徴量（両者の差分）を抽出する第２特徴量抽出手段と、抽出された第１特徴量及び第２特徴量から、所定の重み付け指標値Ｔを画素ごとに算出する指標値算出手段と、算出された重み付け指標値Ｔに基づいて、画像Ｉ０、選択的平滑化画像Ｉ１及び非鮮鋭画像Ｉ２のうちの少なくともいずれかの画像の画素の信号値を用いた新たな信号値を画素ごとに算出し、算出した新たな信号値の画素からなる合成画像を生成する合成画像生成手段と、を備えたものとなっている。
【００３７】
選択的平滑化画像Ｉ１は、構造は比較的鮮明であるが、ノイズ成分がわずかに残ったものとなってしまう。一方、非鮮鋭画像Ｉ２は、ノイズ成分はほぼ除去された状態となっているが、エッジが不鮮明な状態となってしまう。しかし、本実施形態のようにすることで、画素ごとに、構造が比較的鮮明な画像（画像Ｉ０又は選択的平滑化画像Ｉ１）又はノイズ成分が比較的少ない画像（選択的平滑化画像Ｉ１又は非鮮鋭画像Ｉ２）のいずれかの画素の信号値を利用し合成画像を生成するので、画像のエッジが不鮮明になるのを防ぎつつ、エッジのない部分のノイズ成分を十分に低減することができる。
【００３８】
また、本実施形態に係る画像処理装置１は、合成画像生成手段が、重み付け指標値Ｔに基づいて、所定の重み付け比率を画素ごとに決定する比率決定手段と、決定された重み付け比率に基づいて、画像Ｉ０、選択的平滑化画像Ｉ１及び非鮮鋭画像Ｉ２のうちの少なくとも二つの画像の画素の信号値の加重平均を新たな信号値とする重み付け手段と、を備えたものとなっている。
合成画像の画素の信号値として、画像Ｉ０又は非鮮鋭画像Ｉ２の画素の信号値をそのまま用いると、その画素にはノイズ成分が低減されることなく混入されてしまうか、エッジが不鮮明なものとなってしまう。一方、選択的平滑化画像Ｉ１の画素の信号値をそのまま用いると、ノイズ成分が若干残りつつエッジも不鮮明な中途半端な状態となってしまう。しかし、本実施形態のようにすることで、より適切な信号値を設定することができる。
【００３９】
また、本実施形態に係る画像処理装置１は、指標値算出手段が、第１特徴量Ｎと第２特徴量Ｃの比の値（Ｃ／Ｎ）等を重み付け指標値Ｔとして算出するものとなっている。
これにより、重み付け指標値Ｔが大きい画素は構造成分が相対的に多く、小さい画素はノイズ成分が相対的に多いことになるので、重み付け指標値Ｔが大きい画素については、非鮮鋭画像Ｉ２よりも構造成分が多い画像（画像Ｉ０と選択的平滑化画像Ｉ１の少なくとも一方の画像）の画素の信号値を用いて新たな信号値を算出し、重み付け指標値Ｔが小さい画素については、画像Ｉ０よりもノイズ成分が少ない画像（選択的平滑化画像Ｉ１と非鮮鋭画像Ｉ２の少なくとも一方の画像）の画素の信号値を用いて新たな信号値を算出する、といったように、新たな信号値の算出に用いる画像の選択を容易に行うことができる。
【００４０】
＜第２実施形態＞
次に、本発明の第２実施形態について説明する。なお、第１実施形態と共通する構成については説明を省略する。
【００４１】
本実施形態に係る画像処理装置１は、重み付け指標値Ｔに基づいて、画像Ｉ０、選択的平滑化画像Ｉ１及び非鮮鋭画像Ｉ２のうちの少なくともいずれかの画像の画素の信号値を用いた新たな信号値を画素ごとに算出し、算出した新たな信号値の画素からなる合成画像を生成していたが、本実施形態に係る画像処理装置１Ａは、重み付け指標値Ｔに基づいて、選択的平滑化画像Ｉ１及び非鮮鋭画像Ｉ２のうちの少なくとも一方の画像の画素の信号値を用いた新たな信号値を画素ごとに算出し、画像Ｉ０の信号値と新たな信号値との差分をとることによって第１実施形態の合成画像と同等の差分画像を得るようにしている。
このため、本実施形態の画像処理装置１Ａは、ノイズ低減処理の流れ（記憶部１４Ａが記憶している処理プログラム）が第１実施形態と異なる。
【００４２】
図６は本実施形態に係る画像処理装置１Ａが実行するノイズ低減処理のフローチャートである。
本実施形態のノイズ低減処理は、ステップＳ７（重み付け指標値Ｔと閾値Ｔ_thとの比較）までの処理が第１実施形態と共通している。ステップＳ７において、重み付け指標値Ｔの値が閾値Ｔ_thよりも大きい場合には（ステップＳ７；Ｙｅｓ）、選択的平滑化画像Ｉ１における対応する画素の信号値と同じ信号値にする、あるいは選択的平滑化画像Ｉ１と非鮮鋭画像Ｉ２の間でより選択的平滑化画像Ｉ１の信号値に近い信号値を設定する（ステップＳ８Ａ）。選択的平滑化画像Ｉ１の画素の信号値に近い新たな信号値は、第１実施形態の説明で挙げたような形で加重平均を算出することにより得ることができる。
【００４３】
一方、ステップＳ７において、重み付け指標値Ｔが閾値Ｔ_thよりも小さい場合には（ステップＳ７；Ｎｏ）、非鮮鋭画像Ｉ２における対応する画素の信号値と同じ値にする、あるいは非鮮鋭画像Ｉ２と選択的平滑化画像Ｉ１との間でより非鮮鋭画像Ｉ２の信号値に近い信号値を設定する（ステップＳ９Ａ）。非鮮鋭画像Ｉ２の画素の信号値に近い新たな信号値は、第１実施形態の説明で挙げたような形で加重平均を算出することにより得ることができる。
【００４４】
ステップＳ８Ａ又はステップＳ９Ａの処理の後は、第１実施形態と同様に、直前に行った処理が対象の画素のうち最後の画素について行ったものであるか否か、すなわち、対象となる全ての画素についてステップＳ８Ａ又はステップＳ９Ａの処理を行ったか否かを判定する（ステップＳ１０）。
ここで、最後の画素について行った処理ではないと判定した場合には（ステップＳ１０；Ｎｏ）、画素を特定するカウント番号を＋１更新する（処理対象の画素を一つ隣に移す）処理を行って（ステップＳ１１）、ステップＳ７の処理に戻る。
【００４５】
一方、ステップＳ１０の処理において、最後の画素について行ったものであると判定した場合には、画像Ｉ０の各画素の信号値と、得られた新たな信号値との差分を取る処理を行い（ステップＳ１２Ａ）、得られた差分の信号値の画素からなる差分画像を生成する（ステップＳ１３Ａ）。
このような処理を行うことによって得られた差分画像は、第１実施形態の合成画像と同様、エッジがはっきりとしつつも、ノイズ成分が十分に低減されたものとなる。
【符号の説明】
【００４６】
１  画像処理装置
１１  制御部（第１，第２平滑化手段、第１，第２特徴量抽出手段、指標値算出手段、合成画像生成手段、重み付け手段）
１２  ＲＡＭ
１３  通信部
１４  記憶部
１５  操作部
１６  表示部
１７  バス
Ｉ０  画像（放射線画像、原画像、帯域制限画像、低解像度近似画像）
Ｉ１  選択的平滑化画像（第１平滑化画像）
Ｉ２  非鮮鋭画像（第２平滑化画像）
Ｐ  画素
【図１】
2018151715.tif 000003
【図２】
2018151715.tif 000004
【図３】
2018151715.tif 000005
【図４】
2018151715.tif 000006
【図５】
2018151715.tif 000007
【図６】
2018151715.tif 000008