(57)【要約】【課題】散乱線推定処理に要する演算負荷を低減し、処理時間の短縮化を図ること。【解決手段】  複数の画素を含む放射線検出部から出力される放射線画像を処理する画像処理装置は、放射線発生部から放射線検出部に放射線が直接到達する直接線領域を放射線画像に基づいて判定する領域判定部と、直接線領域と判定された領域の画素値を修正し、放射線画像に含まれる散乱線を推定する散乱線推定部と、を備える。【選択図】図２
2020089599.tif 000002
【特許請求の範囲】
【請求項１】
複数の画素を含む放射線検出手段から出力される放射線画像を処理する画像処理装置であって、
放射線発生手段から前記放射線検出手段に放射線が直接到達する直接線領域を前記放射線画像に基づいて判定する領域判定手段と、
前記直接線領域と判定された領域の画素値を修正し、前記放射線画像に含まれる散乱線を推定する散乱線推定手段と、
を備えることを特徴とする画像処理装置。
【請求項２】
前記領域判定手段は、前記放射線画像における画素値に基づいて、前記直接線領域と、当該直接線領域以外の領域とを判定することを特徴とする請求項１に記載の画像処理装置。
【請求項３】
前記領域判定手段は、前記判定の結果に基づいて、前記直接線領域の分布を示す直接線領域マップを取得することを特徴とする請求項１または２に記載の画像処理装置。
【請求項４】
前記直接線領域の画素値を修正した画像を取得する画素値修正手段を更に備え、前記画素値修正手段は、前記直接線領域の画素値を、前記散乱線が発生し得ない画素値に修正することを特徴とする請求項１乃至３のいずれか１項に記載の画像処理装置。
【請求項５】
前記散乱線が発生し得ない画素値には、前記放射線発生手段から被写体を介さず前記放射線検出手段に到達した到達線量相当の画素値が含まれ、
前記画素値修正手段は、前記直接線領域の画素値を、前記到達線量相当の画素値に修正することを特徴とする請求項４に記載の画像処理装置。
【請求項６】
前記画素値修正手段は、前記到達線量相当の画素値を、前記放射線発生手段に設定されている撮影条件、または、前記放射線画像の画像処理の結果に基づいて取得することを特徴とする請求項５に記載の画像処理装置。
【請求項７】
前記画素値修正手段は、前記直接線領域の画素値を修正した前記画像を取得することを特徴とする請求項４乃至６のいずれか１項に記載の画像処理装置。
【請求項８】
前記散乱線推定手段は、前記放射線画像に基づいて、被写体厚の情報を示す体厚画像を取得し、
散乱線推定手段は、前記領域判定手段の判定結果に基づいて、前記直接線領域では散乱線が発生しないように前記体厚画像の画素値を修正することを特徴とする請求項１に記載の画像処理装置。
【請求項９】
前記散乱線推定手段は、前記直接線領域では前記散乱線が発生しないように前記体厚画像の画素値をゼロに修正することを特徴とする請求項８に記載の画像処理装置。
【請求項１０】
前記散乱線推定手段は、前記画素値を修正した前記体厚画像に基づいて、前記放射線画像に含まれる散乱線を推定することを特徴とする請求項８または９に記載の画像処理装置。
【請求項１１】
前記散乱線の散乱線推定処理により取得された散乱線推定画像に対して画素値の調整処理を行う調整処理手段を更に備え、
前記調整処理手段は、前記散乱線推定画像と、前記直接線領域の分布を示す情報に基づいて、前記散乱線推定画像における直接線領域の画素値を調整することを特徴とする請求項１乃至１０のいずれか１項に記載の画像処理装置。
【請求項１２】
前記調整処理手段は、前記直接線領域以外の領域との境となる直接線領域の散乱線推定画像の画素値を、直接線領域以外の領域の画素値に近づけるように滑らかに減少させて前記画素値の調整を行うことを特徴とする請求項１１に記載の画像処理装置。
【請求項１３】
前記調整処理手段は、前記直接線領域から前記直接線領域以外の領域への距離が近づくにつれて、前記直接線領域の画素値の減少の度合いを大きくして前記画素値の調整を行うことを特徴とする請求項１２に記載の画像処理装置。
【請求項１４】
前記調整処理手段は、ローパスフィルタ処理により前記散乱線推定画像から取得した低周波成分の散乱線推定画像と、前記ローパスフィルタ処理により前記直接線領域の分布を示す情報から取得した、直接線の比率を示す直接線比率画像とに基づいて、前記画素値の調整を行うことを特徴とする請求項１１乃至１３のいずれか１項に記載の画像処理装置。
【請求項１５】
複数の画素を含む放射線検出手段から出力される放射線画像を処理する画像処理方法であって、
放射線発生手段から前記放射線検出手段に放射線が直接到達する直接線領域を前記放射線画像に基づいて判定する工程と、
前記直接線領域と判定された領域の画素値を修正し、前記放射線画像に含まれる散乱線を推定する工程と、
を有することを特徴とする画像処理方法。
【請求項１６】
コンピュータを、請求項１乃至１４のいずれか１項に記載の画像処理装置の各手段として機能させるプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像処理装置、画像処理方法およびプログラムに関するものである。
【背景技術】
【０００２】
近年ではFlat Panel Detector(以下、「ＦＰＤ」と略す)と呼ばれる、放射線を電気信号に変換するための半導体素子を２次元行列状に多数配設させたものを用いた放射線撮影装置が広く普及している。
【０００３】
放射線撮影装置で被写体を撮影した際、ＦＰＤに入る放射線は、主に放射線源からＦＰＤまで直進して到達する１次放射線と、コンプトン効果により被写体内で放射線の方向が変化した後、ＦＰＤに到達する２次放射線（以降、「散乱線」と記す）と、主に２種類の成分に分けられる。
【０００４】
画像処理で放射線画像から散乱線成分を推定し、推定した散乱線成分を放射線画像から減算することで散乱線を低減する方法(以降、「散乱線低減処理」)では、放射線発生部から放射線検出部に放射線が直接到達する領域（以降、「直接線領域」）といった、散乱線の特性が被写体領域と異なる領域があると散乱線低減処理における推定精度に影響が生じ得る。
【０００５】
特許文献１には、直接線領域の影響を取り除く方法として、直接線領域以外の領域から散乱線成分を算出するフィルタ処理が開示されている。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開平１-２５８０７４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、散乱線推定処理において、直接線領域と直接線領域以外の領域とで条件分岐した処理を実行する場合、演算負荷が増大することにより散乱線推定処理に長時間を要する場合が生じ得る。
【０００８】
本発明は上記の課題に鑑み、散乱線推定処理に要する演算負荷を低減し、処理時間の短縮化を図ることが可能な画像処理技術を提供する。
【課題を解決するための手段】
【０００９】
本発明の一態様に係る画像処理装置は、複数の画素を含む放射線検出手段から出力される放射線画像を処理する画像処理装置であって、放射線発生手段から前記放射線検出手段に放射線が直接到達する直接線領域を前記放射線画像に基づいて判定する領域判定手段と、
前記直接線領域と判定された領域の画素値を修正し、前記放射線画像に含まれる散乱線を推定する散乱線推定手段と、を備えることを特徴とする。
【００１０】
本発明の他の態様に係る画像処理方法は、複数の画素を含む放射線検出手段から出力される放射線画像を処理する画像処理方法であって、放射線発生手段から前記放射線検出手段に放射線が直接到達する直接線領域を前記放射線画像に基づいて判定する工程と、
前記直接線領域と判定された領域の画素値を修正し、前記放射線画像に含まれる散乱線を推定する工程と、を有することを特徴とする。
【発明の効果】
【００１１】
散乱線低減処理において、散乱線推定処理に要する演算負荷を低減し、処理時間の短縮化を図ることが可能になる。
【図面の簡単な説明】
【００１２】
【図１】第１実施形態に係る放射線撮影装置の構成例を示す図。
【図２】第１実施形態に係る画像処理部の各機能構成の処理を説明する図。
【図３Ａ】直接線領域判定処理により、直接線領域と、それ以外の領域とに分けた結果を模式的に示す図。
【図３Ｂ】画素値が修正された結果を模式的に示す図。
【図３Ｃ】散乱線の分布率を取得する処理を説明する図。
【図４】第２実施形態に係る画像処理部の各機能構成の処理を説明する図。
【図５】第２実施形態に係る散乱線推定処理の流れを説明する図。
【図６】第３実施形態に係る放射線撮影装置の構成例を示す図。
【図７】第３実施形態に係る画像処理部の各機能構成の処理を説明する図。
【図８Ａ】被写体領域及び直接線領域の画素値の分布を模式的に示す図。
【図８Ｂ】画素値が修正された結果を模式的に示す図。
【図９】散乱線調整処理の流れを説明する図。
【図１０】調整処理部の内部処理の流れを示す図。
【図１１】直接線領域における画素値の調整結果を模式的に示す図。
【発明を実施するための形態】
【００１３】
以下、図面を参照して、本発明の実施形態を例示的に詳しく説明する。ただし、この実施形態に記載されている構成要素はあくまで例示であり、本発明の技術的範囲は、特許請求の範囲によって確定されるのであって、以下の個別の実施形態によって限定されるわけではない。
【００１４】
［第１実施形態］
本実施形態では、放射線撮影装置の構成例として散乱線推定処理の前段で、直接線領域の画素値を修正する構成について説明する。
【００１５】
図1は、実施形態に係る放射線撮影装置１０の構成例を示す図である。まず、放射線管１００から被写体１５の延長線上に配置されているＦＰＤ２００（放射線検出部）に放射線を照射する。放射線を照射されたＦＰＤ２００は、放射線を撮影画像に変換し、変換した撮影画像を画像処理部３００のＩ／Ｏ部３０１へ送信する。この時、複数の画素を含むＦＰＤ２００（放射線検出部）から出力される撮影画像（以降、放射線画像ともいう）を処理する画像処理部３００は、放射線管１００から、線量や管電圧などの画像撮影時の撮影条件等に関する撮影情報を取得することが可能である。操作部５００は、画像処理部３００の操作や撮影画像の入力等に使用される。
【００１６】
Ｉ／Ｏ部３０１（取得部）は、ＦＰＤ２００から送信された撮影画像の取得と、放射線管１００から送信された撮影情報の取得を行い、記憶媒体３０２に保存する。記憶媒体３０２には、画像処理プログラムの実行により各種の処理を実行する機能構成として、放射線管１００（放射線発生部）からＦＰＤ２００（放射線検出部）に放射線が直接到達する直接線領域を放射線画像に基づいて判定する領域判定部３０３、画素値の修正処理を行う画素値修正部３０４、散乱線推定部３０５及び散乱線減算部３０６が保存されている。制御部として機能するＣＰＵ３０９の制御の下に、各機能構成（領域判定部３０３、画素値修正部３０４、散乱線推定部３０５及び散乱線減算部３０６）は、記憶媒体３０２に保存されている撮影情報および撮影画像を、作業領域として機能するメモリ３０８に一時的に読み出し、各種の演算処理を実行することが可能である。
【００１７】
ＣＰＵ３０９（制御部）は画像処理の結果を表示部４００に表示するための表示制御を行うこが可能である。尚、ＣＰＵ３０９の代わりにＧＰＵや画像処理用チップ等の演算装置を用いることも可能である。
【００１８】
次に、画像処理部３００の各機能構成の処理について、図２に示したブロック図を用いて説明する。放射線管１００から照射された放射線が被写体１５を透過してＦＰＤ２００に到達すると、ＦＰＤ２００は、放射線を撮影画像に変換し、変換した撮影画像を画像処理部３００へ転送する（Ｓ１０１）。画像処理部３００において、領域判定部３０３は、撮影画像（Ｓ１０１）から、直接線領域の情報（直接線領域情報）を取得する直接線領域の判定処理を行う（Ｓ１０２）。直接線領域の判定処理により、撮影画像を直接線領域と、それ以外の領域とに分けることができる。
【００１９】
次に、画素値修正部３０４は、直接線領域判定処理（Ｓ１０２）で取得した直接線領域情報に基づいて、撮影画像（Ｓ１０１）における直接線領域の画素値を修正した画像（修正画像）を取得する画素値修正処理を行う（Ｓ１０３）。画素値修正処理の詳細については後に説明する。
【００２０】
散乱線推定部３０５は、直接線領域と判定された領域の画素値を修正し、放射線画像に含まれる散乱線を推定する。本実施形態では、散乱線推定部３０５は、画素値修正処理（Ｓ１０３）で生成された修正画像に基づいて、散乱線を推定する散乱線推定処理を行う（Ｓ１０４）。そして、散乱線減算部３０６は、撮影画像（Ｓ１０１）と、散乱線推定処理（Ｓ１０４）で取得した散乱線推定画像とを用いて、撮影画像から散乱線推定画像を除いた散乱線低減画像を取得する散乱線減算処理を行う（Ｓ１０６）。
【００２１】
（直接線領域判定処理：Ｓ１０２）
次に、画像処理プログラムの各処理について個々の具体例を説明する。まず、領域判定部３０３が、Ｓ１０１で取得した撮影画像から直接線領域情報を取得する直接線領域判定処理（Ｓ１０２）を実行する。領域判定部３０３は、撮影画像（放射線画像）における画素値に基づいて、撮影画像から直接線領域と、当該直接線領域以外の領域とを判定する。例えば、撮影画像における画像値の上位から一定比率の値を閾値とし、閾値を超えた画素を直接線領域とすることが可能である。あるいは、画素値のヒストグラムから直接線領域とそれ以外の画素値をｋ-means法によって、直接線領域と、直接線領域以外の領域とを分別してもよい。
【００２２】
領域判定部３０３は、判定の結果に基づいて、直接線領域の分布を示す直接線領域情報（直接線領域マップ）を取得する。直接線領域情報（直接線領域マップ）は、例えば直接線領域を示す画素を「１」、それ以外の領域を「０」とした画像として区別できるようにすることで実現できる。直接線領域を示す画素を「１」、それ以外の領域を「０」とした画像は、直接線領域マップとして取得され、直接線領域判定処理により、撮影画像を直接線領域と、それ以外の領域に分けることができる。
【００２３】
図３Ａは、領域判定部３０３の直接線領域判定処理により、直接線領域と、それ以外の領域（被写体領域）とに分けた結果を模式的に示す図である。図３Ａでは、放射線管１００からＦＰＤ２００まで直進して到達する１次放射線３１に基づいた画素値は、例えば、ＦＰＤ２００の飽和画素値３２であり、ＦＰＤ２００に到達した放射線量（到達線量）相当の画素値３３に比べて低い画素値となる。
【００２４】
被写体１５に到達した放射線３４は、被写体１５内で放射線の方向が変化した散乱線３５となりＦＰＤ２００に到達する。被写体領域の画素値３６は被写体１５を透過した画素値を示す。
【００２５】
直接線領域における画素値（例えば、飽和画素値３２）は、到達線量相当の画素値３３に比べて低いため、散乱線推定処理ではこの領域に被写体が存在するかのような演算処理が実行され得る。直接線領域と被写体領域とを区別しやすくするため、本実施形態では、次の画素値修正処理により直接線領域の画素値（例えば、飽和画素値３２）を修正する処理を行う。
【００２６】
（画素値修正処理：Ｓ１０３）
画素値修正部３０４は、直接線領域判定処理（Ｓ１０２）で取得した直接線領域情報に基づいて、撮影画像（Ｓ１０１）における直接線領域の画素値を修正した画像（修正画像）を取得する画素値修正処理（Ｓ１０３）を実行する。ここで、画素値修正処理（Ｓ１０３）において、画素値の修正とは、直接線領域判定処理Ｓ１０２で取得した直接線領域の画素値を、散乱線が発生し得ない画素値に修正する処理（変更する処理）をいう。散乱線が発生し得ない画素値には、例えば、放射線管１００（放射線発生部）から被写体１５を介さずＦＰＤ２００（放射線検出部）に到達した到達線量相当の画素値が含まれ、画素値修正部３０４は、直接線領域の画素値を、到達線量相当の画素値に修正する。画素値修正部３０４は、直接線領域の画素値を修正した画像を取得する。
【００２７】
画素値修正部３０４は、到達線量相当の画素値を、放射線管１００（放射線発生部）に設定されている撮影条件や、Ｓ１０１で取得した撮影画像（放射線画像）の画像処理の結果に基づいて取得することが可能であり、画素値修正部３０４は、取得した到達線量相当の画素値に基づいて画素値修正処理（Ｓ１０３）を実行する。
【００２８】
図３Ｂは、画素値修正部３０４により、画素値が修正された結果を模式的に示す図である。画素値修正処理は、直接線領域の画素値を、散乱線推定処理Ｓ１０４における演算処理で散乱線が発生し得ない画素値に修正することで実現できる。図３Ｂの例では、直接線領域における画素値（飽和画素値３２）が、到達線量相当の画素値３３に変更されている例を示している。図３Ｂに示す画素値の修正において、画素値修正部３０４は、飽和画素値３２に対して所定値を加算する加算処理ではなく、直接線領域の画素値（飽和画素値３２）自体を、散乱線が発生し得ない画素値として到達線量相当の画素値３３に変更している。すなわち、画素値修正部３０４は、直接線領域の画素値を書き換えている。
【００２９】
散乱線が発生し得ない画素値は、散乱線推定処理（Ｓ１０４）によって異なるため、散乱線推定処理（Ｓ１０４）の具体的処理を示してから、その後に散乱線推定処理に適用する画素値修正処理Ｓ１０３の具体的な処理を説明する。
【００３０】
（散乱線推定処理：Ｓ１０４）
まず、散乱線推定部３０５が実行する散乱線推定処理について、具体例として（１）式を用いた散乱線推定処理（Ｓ１０４）を説明する。
【００３１】
【数１】
2020089599.tif 000003
【００３２】
（１）式において、Ｓ（ｘ,ｙ）は散乱線推定画像の画素値であり、Ｖcは修正画像の画素値である。ｘ,ｙはそれぞれＸ軸とＹ軸の注目画素の座標を示し、iとjはそれぞれ修正画像のＸ軸とＹ軸の座標を示す。
【００３３】
ｎとｍはそれぞれ修正画像のＸ軸とＹ軸の画像サイズであり、ｋは撮影時の管電圧、ｇは撮影時の線量、ｄは撮影距離である。また、ＶＴは撮影条件と撮影画像の画素値に依存した散乱線総発生量であり、ＷはＸ軸とＹ軸の距離に対する散乱線の分布率を示す。
【００３４】
（１）式において、散乱線推定部３０５は、注目画素の座標（ｘ,ｙ）の散乱線推定画像の画素値Ｓ（ｘ,ｙ）を、散乱線総発生量ＶＴと、散乱線の分布率Ｗの畳み込み積分に基づいて取得することができる。（１）式を詳しく説明すると、座標（ｉ,ｊ）から発生する散乱線の総量である散乱線総発生量ＶＴと、散乱線の分布率Ｗの積算が、座標（ｉ,ｊ）から座標（ｘ,ｙ）へ到達した散乱線量を示しており、散乱線推定部３０５は、それらの総和から散乱線推定画像の画素値Ｓ（ｘ,ｙ）を取得する。ここで、管電圧ｋ、線量ｇ、撮影距離ｄの撮影条件は、操作部５００からユーザに入力されたものを用いてもよいし、事前に設定した撮影条件の代表値を用いることが可能である。
【００３５】
なお、公知の手法では、直接線領域の影響を受けない様にするため、（１）式の座標（ｉ,ｊ）が直接線領域であるとき、この領域を計算から除外する必要があった。そのため、プログラムとしては畳み込み積分を行う際に除外判定のための条件分岐が必要となり処理時間がかかるという課題があった。そこで、本発明の実施形態では画素値修正処理（Ｓ１０３）にて予め直接線領域の影響を受けない画素値に修正することで、条件分岐が不要な構成としている。画素値修正処理（Ｓ１０３）による画素値の修正により、散乱線推定処理に要する演算負荷を低減することが可能になり、処理時間の短縮化を図ることが可能になる。
【００３６】
次に、（１）式における散乱線発生量ＶＴと散乱線の分布率Ｗを具体的に説明する。散乱線総発生量ＶＴは、例えば、事前の実験で均一な厚さの人体に近い放射線減衰率の被写体に対し、散乱線除去率の高いグリッドを挿入したグリッド撮影画像と、グリッドを挿入しない非グリッド撮影画像と、グリッドの１次放射線透過率とに基づいて、取得することができる。
【００３７】
均一な厚さの被写体を撮影した場合、散乱線総発生量ＶＴは、周囲から受けた散乱線の総量である散乱線の画素値ＳＴと、おおよそ等しいと言える。これを利用して、散乱線の画素値ＳＴを求めることにより、散乱線総発生量ＶＴを導出する手法の説明をする。散乱線の画素値ＳＴは、以下の（２）式、（３）式の様に求めることができる。
【００３８】
【数２】
2020089599.tif 000004
【００３９】
（２）式において、ｋ,ｇ,ｄの定義は（１）式と同様である。ＡＮは撮影条件が管電圧ｋ、線量ｇ、撮影距離ｄのときの非グリッド撮影画像の平均画素値であり、ＳＴは、撮影条件が管電圧ｋ、線量ｇ、撮影距離ｄのときの非グリッド撮影画像の平均画素値がＡＮのときの散乱線の画素値であり、Ｐは撮影条件が管電圧ｋ、線量ｇ、撮影距離ｄのときの１次放射線画像の平均画素値を示す。
【００４０】
（２）式で表される様に、散乱線の画素値ＳＴは、非グリッド撮影画像の平均画素値から１次放射線画像の平均画素値Ｐ（k,ｇ,ｄ）を減算することで得られる。ここで、（２）式の１次放射線画像の平均画素値Ｐ（k,ｇ,ｄ）は、以下（３）式で取得することができる。
【００４１】
【数３】
2020089599.tif 000005
【００４２】
（３）式において、ＡＧは撮影条件が電圧ｋ、線量ｇ、撮影距離ｄのときのグリッドを挿入した撮影画像（グリッド撮影画像）の平均画素値を示し、ＴＰはグリッドの１次放射線透過率を示す。
【００４３】
（３）式に示すように、１次放射線画像の平均画素値Ｐは、散乱線が除去されたグリッド撮影画像の平均画素値ＡＧに、グリッド特性である１次放射線透過率ＴＰの逆数を積算し、グリッドによる放射線の減衰を補正して取得することができる。尚、（３）式において、１次放射線透過率ＴＰはグリッドのメーカー公表値を用いてもよいし、ＪＩＳ規格（ＪＩＳ  Ｚ  ４９１０：２０００）で示されている測定方法によって取得することも可能である。
【００４４】
上記の（２）式で取得した散乱線の画素値ＳＴ(ＡＮ、ｋ、ｇ、ｄ)を、散乱線総発生量ＶＴ(ＡＮ、ｋ、ｇ、ｄ)と等しいものとして代用する。具体的には、実験で得た条件を様々に変更したときの画像から、散乱線総発生量ＶＴと等しい散乱線の画素値ＳＴを算出し、テーブル化しておくことも可能である。実施方法の例としては、例えば、テーブルをそのまま用いてもよいし、テーブルを撮影条件と撮影画像の画素値を変数とした多変量多項式によりテーブルの値を近似してもよい。また、テーブルを用いる場合は、離散的な撮影条件データの間を線形補間するなどしてテーブルの値を補ってもよい。
【００４５】
散乱線の分布率Ｗは、例えば、図３ＣのＳＴ１に示すように均一な被写体１５の前にピンホールの空いた放射線を遮蔽する遮蔽板３１０を配置して撮影した、ＦＰＤ２００の画像から取得することができる。
【００４６】
具体的な例としては、例えば、図３ＣのＳＴ２に示すような画像のプロファイルのように、撮影画像のピンホール部分は１次放射線の画素値なので除外した散乱線の分布画像を作り、更に、図３ＣのＳＴに示すような散乱線の分布画像の画素値の総和で正規化することで散乱線の分布率Ｗを取得することができる。散乱線の分布率Ｗは２次元のテーブルとして用いても良いし、例えば式（４）のように近似式を用いてもよい。
【００４７】
【数４】
2020089599.tif 000006
【００４８】
ここで、ｒｘとｒｙはそれぞれＸ軸とＹ軸の散乱線の発生画素から到達画素までの距離を示し、βは散乱線の分布率のテーブルを近似するための係数を示す。
【００４９】
（４）式の近似式の例は、散乱線の発生画素から到達画素までの距離に対する散乱線の分布率を、ガウシアン関数を基に近似したものである。係数βの取得方法は、例えば最小二乗法や最尤法を用いて取得することが可能である。
【００５０】
次に、上記の散乱線推定処理（Ｓ１０４）の演算処理に適用するように、直接線領域の画素値を変更する画素値修正処理を説明する。画素値修正部３０４が実行する画素値修正処理（Ｓ１０３）により、直接線領域判定処理（Ｓ１０２）で取得した直接線領域情報に基づいて、Ｓ１０１で取得した撮影画像における直接線領域の画素値を修正した修正画像が生成される。修正画像では、直接線領域の画素は、散乱線が発生し得ない画素値に修正されている。
【００５１】
具体的には、(１)式の散乱線推定画像の画素値Ｓ（ｘ,ｙ）を求める際、座標（ｉ,ｊ）が直接線領域であった場合に、散乱線が発生しないような画素値として、散乱線総発生量ＶＴが０になるような修正画像の画素値Ｖｃに直接線領域の画素値を修正すればよい。（２）式に示すように、撮影条件によらず必ず散乱線総発生量ＶＴが０になる条件は、非グリッド撮影画像の平均画素値ＡＮが０でかつ、１次放射線画像の平均画素値Ｐが０の場合である。
【００５２】
但し、１次放射線画像の平均画素値Ｐは非グリッド撮影画像の平均画素値ＡＮの散乱線が入っていない画素値を示すので必ずＡＮ以下になり、非グリッド撮影画像の平均画素値ＡＮが０であれば必ず散乱線総発生量ＶＴは０になる。つまり、（１）式において、散乱線総発生量ＶＴを０にするには、（２）式の散乱線総発生量ＶＴが０になる非グリッド撮影画像の平均画素値ＡＮの値０を、修正画像の画素値ＶＣに代入すればよい。
【００５３】
上記の関係を式にすると、修正画像の画素値は、（５）式のように、撮影画像に対して直接線領域の画素値を修正することで取得することができる。散乱線推定処理（Ｓ１０４）では、（５）式の関係を（１）式の修正画像の画素値ＶＣに適用することで、Ｓ１０３における画素値修正処理の結果を散乱線推定処理に反映することができる。
【００５４】
【数５】
2020089599.tif 000007
【００５５】
（５）式において、ｉとｊはそれぞれ画像のＸ軸とＹ軸の座標を示し、Ｖｃは修正画像のｉ，ｊ座標の画素値を示し、Ｖinは、Ｓ１０１で取得した撮影画像のｉ,ｊ座標の画素値を示す。
【００５６】
(５)式のように、Ｓ１０１で取得した撮影画像に画素値修正処理（Ｓ１０３）を行った修正画像に基づいて、散乱線推定処理（Ｓ１０４）を行うことで、散乱線推定処理における条件分岐の演算処理は低減される。これにより、散乱線推定画像の取得において、散乱線推定処理に要する演算負荷を低減することが可能になり、処理時間の短縮化を図ることが可能になる。
【００５７】
（散乱線減算処理：Ｓ１０５）
次に、散乱線減算部３０６は、散乱線推定処理（Ｓ１０４）で取得した散乱線推定画像に基づいて、撮影画像（Ｓ１０１）に含まれる散乱線成分を減算（低減）する散乱線減算処理（Ｓ１０５）を実行する。具体的な例としては、（６）式のような様な計算式に基づいて、散乱線減算部３０６は、散乱線を低減した散乱線低減画像を取得することができる（Ｓ１０６）。
【００５８】
【数６】
2020089599.tif 000008
【００５９】
ここで、Ｖｏｕｔはｘ，ｙ座標の散乱線低減画像の画素値であり、Ｖｉｎは撮影画像の画素値であり、αは散乱線透過率、Ｓはｘ，ｙ座標の散乱線推定画像の画素値を示す。
【００６０】
（６）式では、撮影画像の画素値Ｖｉｎから散乱線推定画像Ｓの画素値を減算し１次放射線画像を求めた１次放射線の項と、散乱線推定画像Ｓに散乱線透過率αを掛けた散乱線の項と、の和から散乱線低減画像を取得している（Ｓ１０６）。ここで、散乱線透過率αは、事前にユーザの好みに合わせて設定してもよいし、グリッドのメーカーが公表している散乱線透過率に合わせて設定してもよい。
【００６１】
散乱線減算部３０６は、取得した散乱線低減画像（Ｓ１０６）をＩ／Ｏ部３０１（入出力部）を通して表示部４００に送り、ＣＰＵ３０９（制御部）は散乱線低減画像を表示部４００に表示するための表示制御を行う。ここで、散乱線減算部３０６は、Ｉ／Ｏ部３０１に送る前に階調処理や周波数強調処理等の画像処理を散乱線低減画像に対して行いユーザに見やすい画像に調整することが可能である。
【００６２】
以上の処理によって、直接線領域の画素値を散乱線推定処理（Ｓ１０４）の前段階で画素値を修正することで、散乱線推定画像の取得において、散乱線推定処理に要する演算負荷を低減することが可能になり、処理時間の短縮化を図ることが可能になる。
【００６３】
［第２実施形態］
第１実施形態では、散乱線推定処理（Ｓ１０４）の前段階で直接線領域の画素値を修正する構成について説明したが、第２実施形態では、散乱線推定処理に、直接線領域の画素値を修正する構成を組み込んだ処理の構成について説明する。
【００６４】
本実施形態における放射線撮影装置１０の構成例は図１と同様であり、画像処理部３００の各機能構成の処理について、図４に示したブロック図を用いて説明する。放射線管１００から照射された放射線が被写体１５を透過してＦＰＤ２００に到達すると、ＦＰＤ２００は、放射線を撮影画像に変換し、変換した撮影画像を画像処理部３００へ転送する（Ｓ１０１）。
【００６５】
画像処理部３００において、領域判定部３０３は、撮影画像（Ｓ１０１）から、直接線領域の情報（直接線領域情報）を取得する直接線領域の判定処理を行う（Ｓ１０２）。直接線領域の判定処理により、撮影画像を直接線領域と、それ以外の領域とに分けることができる。
【００６６】
次に、散乱線推定部３０５は、Ｓ１０１で取得した撮影画像と、直接線領域判定処理（Ｓ１０２）で取得した直接線領域の情報（直接線領域情報）と、を用いて散乱線を推定する散乱線推定処理を行う（Ｓ１０４）。そして、散乱線減算部３０６は、撮影画像（Ｓ１０１）と、散乱線推定処理（Ｓ１０４）で取得した散乱線推定画像とを用いて、撮影画像から散乱線推定画像を除いた散乱線低減画像を取得する散乱線減算処理を行う（Ｓ１０６）。
【００６７】
第２実施形態の構成と第１実施形態との構成の違いは、画像処理プログラムにおける散乱線推定処理（Ｓ１０４）であるので、以下に散乱線推定処理Ｓ１０４の具体的な処理を説明する。
【００６８】
（散乱線推定処理：Ｓ１０４）
図５は、第２実施形態に係る散乱線推定処理の流れを説明する図であり、散乱線推定処理（Ｓ１０４）では、Ｓ１０１で取得した撮影画像と、直接線領域判定処理（Ｓ１０２）で取得した直接線領域情報と、を用いて散乱線推定画像を取得する。
【００６９】
まず、Ｓ２０１において、散乱線推定部３０５は、撮影画像（Ｓ１０１）から仮の１次放射線画像を取得するために１次放射線画像取得処理を行う。
【００７０】
次に、Ｓ２０２において、散乱線推定部３０５は、１次放射線画像取得処理（Ｓ２０１）で取得した１次放射線画像と、Ｓ１０２で取得した直接線領域情報から体厚画像を取得する体厚画像取得処理を行う。
【００７１】
次に、Ｓ２０３において、散乱線推定部３０５は、体厚画像取得処理（Ｓ２０２）により取得した体厚画像を用いて、仮の散乱線推定画像を取得する仮散乱線推定画像取得処理を行う。
【００７２】
次に、Ｓ２０４において、散乱線推定部３０５は、撮影画像（Ｓ１０１）と、１次放射線画像（Ｓ２０１）と、Ｓ２０３で取得した仮の散乱線推定画像とに基づいて、１次放射線画像を更新する１次放射線画像更新処理を実行する。
【００７３】
Ｓ２０５において、散乱線推定部３０５は、１次放射線画像更新処理（Ｓ２０４）で更新される前の仮の１次放射線画像と更新後の１次放射線画像とに基づいて、ループ終了判定処理を行う。ループ終了判定処理Ｓ２０５でループ終了の判定となった場合は（Ｓ２０５−Ｙｅｓ）、散乱線推定処理Ｓ１０４を終了し、ループ未終了の判定となった場合は（Ｓ２０５−Ｎｏ）、散乱線推定部３０５は、処理をステップＳ２０２（体厚画像取得処理）に戻し、ループ終了判定になるまで同様の処理を繰り返す。
【００７４】
そして、Ｓ２０６において、散乱線推定部３０５は、ループの最後に更新された１次放射線画像と撮影画像の画素値を用いて散乱線推定画像を取得する。
【００７５】
次に、図５で説明したフローチャートの各ステップにおける具体的な処理方法を説明する。
【００７６】
図５の１次放射線画像取得処理（Ｓ２０１）では、被写体１５の被写体厚の情報を示す体厚画像を取得するために仮の１次放射線画像を取得する。Ｓ２０１で取得する仮の１次放射線画像は、後の１次放射線画像更新処理（Ｓ２０４）で真値に近づけるため、あくまで初期値としての１次放射線画像を設定する処理となる。
【００７７】
具体的な例としては、１次放射線画像として、Ｓ１０１で取得した撮影画像をそのまま用いてもよいし、１次放射線画像の真値に近づけるために事前に１次放射線含有率を設定しておき仮の１次放射線画像取得時において撮影画像（Ｓ１０１）に１次放射線含有率を積算して、仮の１次放射線画像を取得してもよい。
【００７８】
次に、図５の体厚画像取得処理（Ｓ２０２）では、仮散乱線推定画像取得処理（Ｓ２０３）で散乱線の算出に用いられる被写体厚の情報を示す体厚画像を取得する。
【００７９】
散乱線推定部３０５は、直接線領域と判定された領域の画素値を修正し、放射線画像に含まれる散乱線を推定する。本実施形態では、散乱線推定部３０５は、撮影画像（放射線画像）に基づいて、被写体厚の情報を示す体厚画像を取得し、散乱線推定部３０５は、領域判定部３０３の判定結果に基づいて、直接線領域では散乱線が発生しないように体厚画像の画素値を修正する。すなわち、散乱線推定部３０５は、直接線領域判定処理（Ｓ１０２）で取得した直接線領域情報（領域判定部３０３の判定結果）に基づいて、直接線領域では散乱線が発生しないように体厚画像の画素値を修正する。具体的な例としては、体厚画像は次の（７）式で取得することができる。（７）式に示すように、散乱線推定部３０５は、直接線領域では散乱線が発生しないように体厚画像の画素値をゼロに修正する。体厚画像における画素値の修正は、散乱線推定処理において、直接線領域の画素値を修正する構成を組み込んだ処理となり、散乱線推定部３０５は、画素値を修正した体厚画像に基づいて、放射線画像に含まれる散乱線を推定する。
【００８０】
【数７】
2020089599.tif 000009
【００８１】
ここで、ｘとｙはそれぞれ画像のＸ軸とＹ軸の座標を示し、Ｂ（ｘ,ｙ）は体厚画像の画素値を示す。また、μは放射線減衰係数を示し、Ｐ（ｘ,ｙ）はｘ,ｙ座標の１次放射線画像の画素値を示し、Ｅは被写体を介さずＦＰＤ２００に到達した直接到達線量値（到達線量相当の画素値）を示す。
【００８２】
（７）式に示したように、１次放射線画像と体厚画像の画素値はランベルトの法則を用いて取得することができる。（７）式の直接到達線量値Ｅは、被写体を介さずＦＰＤ２００に到達した放射線の画素値を示す。具体的な例としては、直接線領域判定処理（Ｓ１０２）で取得した直接線領域情報に基づいて、撮影画像（Ｓ１０１）における直接線領域の画素値の平均値を使ってもよいし、ＮＤＤ法を使っての撮影条件等から到達画素値を推定する等してもよい。
【００８３】
（７）式の放射線減衰係数μは、被写体の放射線減衰係数の代表値を示す。例えば、事前に実験で得られた脂肪の放射線減衰係数と筋肉の放射線減衰係数の平均値を用いることが可能である。
【００８４】
図５の仮散乱線推定画像取得処理（Ｓ２０３）では、例えば、（８）式に示すような散乱線の発生に関する近似式を用いて仮の散乱線推定画像を取得する。
【００８５】
【数８】
2020089599.tif 000010
【００８６】
ここで、（８）式において、Ｓ（ｘ,ｙ）は散乱線推定画像の画素値であり、Ｂは体厚画像の画素値である。ｘ,ｙはＸ軸とＹ軸のそれぞれの座標であり、iとjはそれぞれ修正画像のＸ軸とＹ軸の座標を示す。また、ｎとｍはそれぞれ画像のＸ軸とＹ軸の画像サイズを示し、αは係数、Ｅは被写体を介さずＦＰＤ２００に到達した直接到達線量値（到達線量相当の画素値）、Ｗは散乱線の分布率を示す。
【００８７】
（８）式に示したように、ある注目画素の散乱線量Ｓ（ｘ,ｙ）は、体厚画像の画素値Ｂと直接到達線量値Ｅとその２つの変数から散乱線の発生量の画素値に換算する係数αからなる散乱線総発生量と、散乱線の分布を示す散乱線分布率Ｗの畳み込み積分で取得することができる。
【００８８】
従来技術では（８）式の座標（ｉ,ｊ）が直接線領域であるとき、直接線領域の影響を受けないようにするため、計算から除外される条件分岐処理が必要とされる。本発明の実施形態では、体厚画像取得処理（Ｓ２０２）にて直接線領域の影響を受けないように直接線領域の画素値が修正されるため、条件分岐が必要なくなっている。
【００８９】
（８）式の係数αは、具体的には(２)式と同じく事前に均一な被写体をグリッド挿入時と非グリッド挿入時で撮影した画像から散乱線総発生量を計算し、散乱線総発生量から（９）式のように取得することができる。
【００９０】
【数９】
2020089599.tif 000011
【００９１】
ここで、αは係数であり、ＶＴは散乱線総発生量を示す。また、Ｅは被写体を介さずＦＰＤ２００に到達した直接到達線量値（到達線量相当の画素値）を示し、Ｂ（ｉ,ｊ）はｉ,ｊ座標における体厚画像の画素値を示す。
【００９２】
また、（８）式の散乱線の分布率Ｗは、第１実施形態で説明した（４）式と同様の近似式である。
【００９３】
次に、図５の１次放射線画像更新処理（Ｓ２０４）では、撮影画像（Ｓ１０１）と、１次放射線画像（Ｓ２０１）と、仮散乱線推定画像取得処理（Ｓ２０３）で取得した仮の散乱線推定画像とに基づいて、１次放射線画像が真値に近づくように更新する。
【００９４】
１次放射線画像更新処理（Ｓ２０４）は、例えば、最尤法や最小二乗法等を用いることにより実行することが可能である。具体的な例として、（１０）式は最尤法を用いた場合の１次放射線画像更新処理を例示するものである。
【００９５】
【数１０】
2020089599.tif 000012
【００９６】
ここで、Ｐ（ｘ,ｙ）は更新前の１次放射線画像の画素値を示し、Ｐ'（ｘ,ｙ）は更新後の１次放射線画像の画素値を示す。Ｓ（ｘ,ｙ）は散乱線推定画像の画素値を示し、Ｖｉｎは撮影画像（Ｓ１０１）の画素値を示す。
【００９７】
（１０）式に示したように、更新後の１次放射線画像の画素値Ｐ'（ｘ,ｙ）は、撮影画像の画素値Ｖｉｎを、更新前の１次放射線画像Ｐ（ｘ,ｙ）と散乱線推定画像Ｓの画素値とを加算した加算値で除算し、除算結果に１次放射線画像Ｐ（ｘ,ｙ）を乗算することにより取得される。
【００９８】
更新後の１次放射線画像の画素値Ｐ'（ｘ,ｙ）は、（１０）式に含まれる（１１）式の因子によって、（１２）式の関係式を満たすように更新することで真値に近づく。
【００９９】
【数１１】
2020089599.tif 000013
【０１００】
【数１２】
2020089599.tif 000014
【０１０１】
次に、図５のループ終了判定処理（Ｓ２０５）では、例えば、（１０）式に示したような最尤法により真値に近い値に収束させていくための繰り返し演算を継続するか否かを判定する。具体的な例としては、１次放射線画像更新処理（Ｓ２０４）の更新前の１次放射線画像Ｐ（ｘ,ｙ）の画素値と更新後の１次放射線画像Ｐ'（ｘ,ｙ）の画素値とを比較し、両者の差分が所定の閾値より小さければループを終了するようにしてもよいし、処理時間が長くなり過ぎない様に、ループ回数を事前に設定する等してもよい。
【０１０２】
次に、図５の散乱線推定画像取得処理（Ｓ２０６）では、ステップＳ２０５において、ループの最後に更新された１次放射線画像と撮影画像の画素値を用いて散乱線推定画像を取得する。
【０１０３】
（１３）式は、散乱線推定画像取得処理（Ｓ２０６）における具体的な処理の例を示す演算式である。
【０１０４】
【数１３】
2020089599.tif 000015
【０１０５】
ここで、Ｓ（ｘ,ｙ）は散乱線推定画像の画素値を示し、Ｖｉｎは撮影画像（Ｓ１０１）の画素値を示す。Ｐ（ｉ,ｊ）はループの最後に更新された１次放射線画像の画素値を示す。
【０１０６】
（１３）式に示すように、散乱線推定画像Ｓ（ｘ,ｙ）の画素値は、撮影画像の画素値Ｖｉｎから１次放射線画像の画素値を減算することにより取得することができる。以上により散乱線推定処理（Ｓ１０４）は終了する。
【０１０７】
本実施形態によれば、散乱線の発生しない直接線領域の影響を排除する散乱線低減処理において、画素間の演算負荷の多い仮散乱線推定画像取得処理（Ｓ２０３）の前段の処理（Ｓ２０２）で、体厚画像における直接線領域の画像値を修正することで、散乱線推定処理における条件分岐の演算処理は低減される。これにより、散乱線推定画像の取得において、散乱線推定処理に要する演算負荷を低減することが可能になり、処理時間の短縮化を図ることが可能になる。
【０１０８】
［第３実施形態］
本実施形態では、第１及び第２実施形態で説明したように直接線領域の画素値を修正する構成に加えて、直接線領域における散乱線の値を調整する構成について説明する。
【０１０９】
図６は本実施形態における放射線撮影装置３０の構成例を示す図であり、画像処理部３００の機能構成として、調整処理部３０７が追加されている点で図１の放射線撮影装置１０の構成と相違する。
【０１１０】
次に画像処理部３００の各機能構成の処理について、図７に示したブロック図を用いて説明する。放射線管１００から照射された放射線が被写体１５を透過してＦＰＤ２００に到達すると、ＦＰＤ２００は、放射線を撮影画像に変換し、変換した撮影画像を画像処理部３００へ転送する（Ｓ１０１）。画像処理部３００において、領域判定部３０３は、撮影画像（Ｓ１０１）から、直接線領域の情報（直接線領域情報）を取得する直接線領域の判定処理を行う（Ｓ１０２）。直接線領域の判定処理により、撮影画像を直接線領域と、それ以外の領域とに分けることができる。
【０１１１】
画素値修正部３０４は、直接線領域判定処理（Ｓ１０２）で取得した直接線領域情報に基づいて、撮影画像（Ｓ１０１）における直接線領域の画素値を修正した画像（修正画像）を取得する画素値修正処理を行う（Ｓ１０３）。この画素値修正処理は第１実施形態で説明した処理と同様の処理である。図７では、第１実施形態の構成に基づいて説明しているが、第２実施形態の図４で説明したように、散乱線推定処理（Ｓ１０４）に画素値の修正処理を組み込んでもよい。
【０１１２】
散乱線推定部３０５は、画素値修正処理（Ｓ１０３）で生成された修正画像に基づいて、散乱線を推定する散乱線推定処理を行う（Ｓ１０４）。
【０１１３】
調整処理部３０７は、散乱線の散乱線推定処理により取得された散乱線推定画像に対して画素値の調整処理を行う。調整処理部３０７は、散乱線推定画像と、直接線領域の分布を示す情報に基づいて、散乱線推定画像における直接線領域の画素値を調整する。ここで、直接線領域の分布を示す情報とは、領域判定部３０３の判定の結果により取得された直接線領域情報（直接線領域マップ）である。すなわち、調整処理部３０７は、散乱線推定処理（Ｓ１０４）により取得された散乱線推定画像と、直接線領域判定処理（Ｓ１０２）で取得された、直接線領域の分布を示す直接線領域情報（直接線領域マップ）に基づいて、散乱線推定画像における直接線領域の画素値を調整する散乱線調整処理を行う（Ｓ１０７）。散乱線調整処理として、調整処理部３０７は、散乱線推定画像における直接線領域と当該直接線領域以外の領域（被写体領域）との間における画素値の差（段差）が滑らかになるように、直接線領域判定処理により判定された直接線領域の画素値を調整する。この散乱線調整処理の詳細については後に説明する。
【０１１４】
そして、散乱線減算部３０６は、撮影画像（Ｓ１０１）と、散乱線調整処理に基づいて、被写体領域との境となる直接線領域の画素値が調整された散乱線推定画像とを用いて、撮影画像から調整された散乱線推定画像を除いた散乱線低減画像を取得する散乱線減算処理を行う（Ｓ１０６）。
【０１１５】
図７に示したブロック図において、直接線領域判定処理（Ｓ１０２）、画素値修正処理（Ｓ１０３）、散乱線推定処理（Ｓ１０４）、散乱線減算処理（Ｓ１０５）の具体的な処理内容は、第１実施形態と同様の処理となる。第１実施形態の構成と相違する散乱線調整処理（Ｓ１０７）について、以下、図８Ａ〜図１１を参照して具体的な処理の内容を説明する。
【０１１６】
図８Ａは、被写体領域及び直接線領域における画素値の分布を模式的に示す図である。図８Ａにおいて、放射線管１００からＦＰＤ２００まで直進して到達する１次放射線３１に基づいた画素値はＦＰＤ２００の飽和画素値３２であり、ＦＰＤ２００に到達した放射線量（到達線量）相当の画素値３３に比べて低い画素値となる。
【０１１７】
被写体１５に到達した放射線３４は、被写体１５内で放射線の方向が変化した散乱線３５となりＦＰＤ２００に到達する。被写体領域の画素値３６は被写体１５を透過した画素値を示す。
【０１１８】
被写体領域の一部が飽和している画像や手など厚さの薄い被写体でシェーディングが発生している画像では、被写体領域及び直接線領域において画素値の分布が同等の分布となる場合が生じ得る。例えば、被写体領域の端部の領域８０１では、直接線領域における飽和画素値３２と同等に被写体領域の画素値が高くなるため、画素値を用いて、被写体領域と直接線領域とを分別することが困難となる場合が生じ得る。
【０１１９】
図８Ｂは、画素値修正部３０４により、画素値が修正された結果を模式的に示す図である。直接線領域における画素値（飽和画素値３２）と同等の画素値を有する被写体領域の端部の領域８０１は、被写体領域に含まれる領域であるが、直接線領域と誤判定され、到達線量相当の画素値３３に変更されている。すなわち、直接線領域における画素値（飽和画素値３２）は、到達線量相当の画素値３３に書き換えられている。誤判定された画素と正しく判定された画素との境界で推定される散乱線量に連続性が無くなり画素値の差（段差）が発生する。この段差がアーチファクトとして見えてしまう場合が生じ得る。
【０１２０】
本実施形態の調整処理部３０７は、散乱線推定画像における直接線領域と被写体領域との間における画素値の差（段差）が滑らかになるように、直接線領域判定処理により判定された直接線領域の画素値を調整する。
【０１２１】
図９は、調整処理部３０７が実行する散乱線調整処理（Ｓ１０７）の流れを説明する図である。Ｓ９０１において、調整処理部３０７は、散乱線推定部３０５から取得した散乱線推定画像にローパスフィルタ（ＬＰＦ）を掛けて低周波散乱線推定画像を取得する（散乱線推定画像の低周波成分の抽出処理）。
【０１２２】
Ｓ９０２において、調整処理部３０７は、Ｓ９０１で用いたローパスフィルタ（ＬＰＦ）と同じＬＰＦを、領域判定部３０３から取得した直接線領域の分布を示す直接線領域情報（直接線マップ）（例えば、直接線領域＝１[LSB],その他の領域＝０）に掛けて、０〜１の値で直接線の比率を示す直接線マップの低周波画像（直接線比率画像）を取得する（直接線マップの低周波成分の抽出処理）。
【０１２３】
そして、Ｓ９０３において、調整処理部３０７は、散乱線推定部３０５から取得した散乱線推定画像と領域判定部３０３から取得した直接線マップと、Ｓ９０１で取得した低周波散乱線推定画像と、Ｓ９０２で取得した直接線比率画像とに基づいて、直接線領域の画素値を調整した散乱線推定画像（調整後散乱線推定画像）を生成する。
【０１２４】
図１０は、Ｓ９０３の処理（散乱線調整処理）を実行する際の調整処理部３０７の内部処理の流れを示すブロック図である。調整処理部３０７は、散乱線推定部３０５から取得した散乱線推定画像１０１０にローパスフィルタ（ＬＰＦ１０１２）を掛けて低周波成分の散乱線推定画像（低周波散乱線推定画像１０１４）を取得する。
【０１２５】
また、調整処理部３０７は、領域判定部３０３から取得した直接線領域マップ１０２０にローパスフィルタ（ＬＰＦ１０２２）を掛けて直接線比率画像１０２４を取得する。ここで、ＬＰＦ１０２２とＬＰＦ１０１２とは同一のローパスフィルタである。
【０１２６】
図１０の画素値調整処理１０３０において、調整処理部３０７は、直接線領域マップ１０２０を飽和画素判定用（領域判定用）として用いる。すなわち、調整処理部３０７は、直接線領域マップ１０２０に基づいて、散乱線推定画像１０１０における直接線領域と、それ以外の領域に分けることができる。
【０１２７】
調整処理部３０７は、画素値調整処理１０３０において直接線領域における画素値を以下の（１４）式に基づいて取得し、直接線領域以外のその他の領域の画素値を（１５）式に基づいて取得する。
【０１２８】
調整処理部３０７は、ローパスフィルタ処理（ＬＰＦ１０１２）により散乱線推定画像１０１０から取得した低周波成分の散乱線推定画像（低周波散乱線推定画像１０１４）と、ローパスフィルタ処理（ＬＰＦ１０２２）により直接線領域の分布を示す情報（直接線領域マップ１０２０）から取得した、直接線の比率を示す直接線比率画像１０２４とに基づいて、画素値の調整を行う。この画素値の調整は以下の（１４）式、（１５）式に基づいて実行される。
【０１２９】
【数１４】
2020089599.tif 000016
【０１３０】
【数１５】
2020089599.tif 000017
【０１３１】
（１４）及び（１５）式において、ｘ,ｙはＸ軸とＹ軸のそれぞれの座標であり、Ｓ（ｘ,ｙ）は散乱線推定画像の画素値を示し、Ｓ'（ｘ,ｙ）は調整後の散乱線推定画像の画素値を示している。また、（１４）式において、ＳLPF（ｘ,ｙ）は、低周波散乱線推定画像を示し、ＤLPF（ｘ,ｙ）は直接線比率画像を示している。
【０１３２】
（１４）式の上段の式は、被写体領域に対して近い側に位置する直接線領域の内側における画素値の分布を示す演算式であり、（１４）式の下段の式は、被写体領域に対して遠い側に位置する直接線領域の外側における画素値の分布を示す演算式である。直接線比率画像ＤLPF（ｘ,ｙ）は、直接領域から被写体領域に近づくに従い、直接線の比率が漸減する傾向を示す画像であり、ＤLPF（ｘ,ｙ）が０．５となるとき、（１４）式における上段と下段の式は一致する。
【０１３３】
（１５）式に示すように、直接線領域以外のその他の領域、すなわち領域判定処理で誤判定されずに正常に被写体領域と判定され領域の画素値は、散乱線調整処理において変更されず、調整後の散乱線推定画像の画素値Ｓ'（ｘ,ｙ）は、散乱線推定画像における画素値Ｓ（ｘ,ｙ）と等しくなるように設定される。
【０１３４】
図１０の画素値調整処理１０３０において、調整処理部３０７は、（１４）式及び（１５）式に基づいて、散乱線推定画像における直接線領域の画素値を調整した調整後散乱線推定画像１０４０を生成する。
【０１３５】
図１１は、調整処理部３０７による直接線領域における画素値の調整結果を模式的に示す図である。調整処理部３０７は、散乱線推定画像において、散乱線推定画像における直接線領域と被写体領域との間における画素値の差（段差）が滑らかになるように、直接線領域判定処理により判定された直接線領域の画素値を調整する。すなわち、図１１に示すように、調整処理部３０７は、直接線領域以外の領域（被写体領域）との境となる直接線領域の画素値（散乱線推定画像の画素値）を、直接線領域以外の領域（被写体領域）の画素値に近づけるように滑らかに減少させて画素値の調整を行う。すなわち、調整処理部３０７は、直接線領域以外の領域（被写体領域）との境となる直接線領域の画素値（散乱線推定画像の画素値）を、直接線領域以外の領域（被写体領域）の画素値に近づけるように画素値を連続的に減少させて画素値の調整を行う。調整処理部３０７は、直接線領域から直接線領域以外の領域（被写体領域）への距離が近づくにつれて、直接線領域の画素値の減少の度合いを大きくして画素値の調整を行う。
【０１３６】
誤判定されずに正常に被写体領域と判定され領域の画素値３６は、散乱線調整処理（Ｓ１０７）において変更されず、散乱線推定画像における画素値Ｓ（ｘ,ｙ）の値が調整後の散乱線推定画像の画素値として設定される。
【０１３７】
被写体領域に含まれる領域であるが、直接線領域と誤判定された領域８０１を含む直接線領域においては、散乱線調整処理（Ｓ１０７）により、直接線領域における修正後の画素値と、誤判定されずに正常に被写体領域と判定された画素値３６とが滑らかに接続される。
【０１３８】
本実施形態によれば、被写体領域が領域判定処理で誤判定され、直接線領域に含まれてしまった場合でも、直接線領域の修正後の画素値と被写体領域の画素値との間における画素値の差（段差）が滑らかになるように、直接線領域判定処理により判定された直接線領域の画素値を調整することができ、これにより、直接線領域と被写体領域との境界に生じ得るアーチファクトを低減することが可能になる。
【０１３９】
以上、本発明の実施形態（第１実施形態〜第３実施形態）を説明したが、本発明はこれらの実施形態に限定されないことはいうまでもなく、その要旨の範囲内で種々の変形及び変更が可能である。散乱線推定処理における具体的な計算式は例示的なものであり、本発明の適用例はこれに限るものではない。
【０１４０】
（その他の実施形態）
本発明は、上述の実施形態の１以上の機能を実現するプログラムを、ネットワーク又は記憶媒体を介してシステム又は装置に供給し、そのシステム又は装置のコンピュータにおける１つ以上のプロセッサーがプログラムを読出し実行する処理でも実現可能である。また、１以上の機能を実現する回路（例えば、ＡＳＩＣ）によっても実現可能である。
【符号の説明】
【０１４１】
１５：被写体、１００：放射線管、２００：ＦＰＤ、３００：画像処理部、３０１：Ｉ／Ｏ部、３０２：記憶媒体、３０３：領域判定部、３０４：画素値修正部、３０５：散乱線推定部、３０６：散乱線減算部、３０７：調整処理部、３０８：メモリ、３０９：ＣＰＵ、４００：表示装置、５００：操作装置
【図１】
2020089599.tif 000018
【図２】
2020089599.tif 000019
【図３Ａ】
2020089599.tif 000020
【図３Ｂ】
2020089599.tif 000021
【図３Ｃ】
2020089599.tif 000022
【図４】
2020089599.tif 000023
【図５】
2020089599.tif 000024
【図６】
2020089599.tif 000025
【図７】
2020089599.tif 000026
【図８Ａ】
2020089599.tif 000027
【図８Ｂ】
2020089599.tif 000028
【図９】
2020089599.tif 000029
【図１０】
2020089599.tif 000030
【図１１】
2020089599.tif 000031