(57)【要約】【課題】  異なる内部構造を持つ放射線検出装置が混在した状態においても構造が写り込んだ領域を適切に補正して合成画像の画質を向上させる放射線撮影システム及び放射線撮影方法を提供する。【解決手段】  放射線を検出する複数の放射線検出装置と、複数の放射線検出装置から取得される複数の放射線画像を合成して合成画像を生成する合成処理部を有した放射線撮影システムにおいて、前記合成画像における前記放射線検出装置の構造が写り込んだ領域を補正する画像補正部を備え、前記画像補正部は前記合成画像に写り込んだ放射線検出装置の構造陰影の特性に応じて補正方法を設定する。【選択図】  図１
2019068915.tif 000002
【特許請求の範囲】
【請求項１】
放射線を検出する複数の放射線検出装置と、複数の放射線検出装置から取得される複数の放射線画像を合成して合成画像を生成する合成処理部を有した放射線撮影システムにおいて、前記合成画像における前記放射線検出装置の構造が写り込んだ領域を補正する画像補正部を備え、前記画像補正部は前記合成画像に写り込んだ放射線検出装置の構造陰影の特性に応じて補正方法を設定することを特徴とする放射線撮影システム。
【請求項２】
前記合成画像に写り込む前記放射線検出装置の構造陰影を識別する識別コードを記憶する記憶部を備え、前記画像補正部は識別コードに応じて前記補正方法を設定することを特徴とする請求項１に記載の放射線撮影システム。
【請求項３】
前記識別コードは、前記放射線検出装置の種類を示す情報を含んでいることを特徴とする請求項２に記載の放射線撮影システム。
【請求項４】
前記識別コードは、前記放射線検出装置の位置情報を含んでいることを特徴とする請求項２又は３に記載の放射線撮影システム。
【請求項５】
前記画像補正部は、前記構造陰影が単純な合成画像の補正に適した補正、前記構造陰影が複雑な合成画像の補正に適した補正の一方を行うことを特徴とする請求項１に記載の放射線撮影システム。
【請求項６】
前記構造陰影が複雑な合成画像の補正に適した補正が行われた場合、前記画像補正部は、前記放射線検出装置の構造が写り込んだ領域に対して粒状性を改善する補正を行うことを特徴とする請求項５に記載の放射線撮影システム。
【請求項７】
前記画像補正部は、前記構造が写り込んだ領域を前記構造陰影の違いに応じて分割し、分割された領域毎に補正方法を設定することを特徴とする請求項１に記載の放射線撮影システム。
【請求項８】
前記画像補正部は、前記合成画像における少なくとも１つ以上の画素において定めたモデルにおけるモデルパラメータを推定し、推定したモデルパラメータを用いて前記画素の画素値を補正することを特徴とする請求項１に記載の放射線撮影システム。
【請求項９】
前記画像補正部は、前記合成画像に写り込んだ放射線検出装置の構造陰影に応じて前記モデルまたはモデルパラメータの推定方法の少なくとも何れか１つを変更することを特徴とする請求項８に記載の放射線撮影システム。
【請求項１０】
前記画像補正部は、前記合成画像に写り込んだ放射線検出装置の構造陰影に応じてノイズの補正を実行するか否かを設定することを特徴とする請求項１に記載の放射線撮影システム。
【請求項１１】
放射線を検出する複数の放射線検出装置と、複数の放射線検出装置から取得される複数の放射線画像を合成して合成画像を生成する合成処理部を有した放射線撮影システムにおいて、前記放射線検出装置の種類と位置情報を識別する識別コードを記憶する記憶部と、前記識別コードに応じて、前記合成画像における前記放射線検出装置の構造が写り込んだ領域を補正する画像補正部とを備えることを特徴とする放射線撮影システム。
【請求項１２】
放射線を検出する複数の放射線検出装置と、複数の放射線検出装置から取得される複数の放射線画像を合成して合成画像を生成する放射線撮影方法において、前記合成画像における前記放射線検出装置の構造が写り込んだ領域を補正する補正ステップを有し、前記補正ステップは前記合成画像に写り込んだ放射線検出装置の構造陰影の特性に応じて補正方法を設定することを特徴とする放射線撮影方法。
【請求項１３】
放射線を検出する複数の放射線検出装置と、複数の放射線検出装置から取得される複数の放射線画像を合成して合成画像を生成する放射線撮影方法において、前記放射線検出装置の種類と位置情報を識別する識別コードを記憶するステップと、前記識別コードに応じて、前記合成画像における前記放射線検出装置の構造が写り込んだ領域を補正するステップを有していることを特徴とする放射線撮影方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、放射線を用いて撮影を行う放射線撮影システム及び放射線撮影方法に関するものである。
【背景技術】
【０００２】
近年、例えば医療分野では被検者の脊椎や下肢の全体を撮影するといった、観察領域が広い撮影（以下、長尺撮影と呼ぶ）が行われている。特許文献１には、複数の放射線検出装置を並べて撮影することで長尺撮影を行うことができ、長尺画像（合成画像とも呼ぶ）における、該重ね合せられた放射線検出装置の構造の写り込みを補正することが開示されている。また、特許文献２では放射線検出装置の構造を工夫することで、該重ね合せられた放射線検出装置の構造の写り込みを軽減する方法が開示されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１６−１４０５１５号公報
【特許文献２】特開２０１７−９４１３１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
長尺撮影で用いられる放射線検出装置の内部構造は様々であり、特許文献２のように構造が単純化された放射線検出装置や複雑な構造を有する放射線検出装置がある。そのため、長尺画像に写り込む構造陰影に応じて合成画像の補正の種類も異なる。しかしながら、特許文献１、２には内部構造が異なる放射線検出装置が混在する場合の対策については言及されていない。
【０００５】
そこで本発明は、内部構造が異なる放射線検出装置が混在した状態においても適切に合成画像を補正することができる放射線撮影システム及び放射線撮影方法を提供することを目的とする。
【課題を解決するための手段】
【０００６】
本発明の目的を達成するために、放射線を検出する複数の放射線検出装置と、複数の放射線検出装置から取得される複数の放射線画像を合成して合成画像を生成する合成処理部を有した放射線撮影システムにおいて、前記合成画像における前記放射線検出装置の構造が写り込んだ領域を補正する画像補正部を備え、前記画像補正部は前記合成画像に写り込んだ放射線検出装置の構造陰影の特性に応じて補正方法を設定する。
【発明の効果】
【０００７】
本発明によれば、内部構造が異なる放射線検出装置が混在した状態においても適切に合成画像を補正することができる。
【図面の簡単な説明】
【０００８】
【図１】本発明の放射線撮影システムの概略構成を示す図。
【図２】本発明の放射線撮影システムの放射線検出装置と画像データの関係を示す図。
【図３】本発明の放射線撮影システム（主に画像表示制御部）の構成を示す図。
【図４】本発明の放射線撮影システムの構造データの取得方法を示す図。
【図５】本発明の実施例１の画像補正部の処理手順を示すフローチャート。
【図６】本発明の実施例１の合成画像と構造情報を示す図。
【図７】本発明の実施例１の欠損領域の位置を説明する図。
【図８】本発明の実施例２の画像補正部の処理手順を示すフローチャート。
【図９】本発明の実施例２の欠損領域の分割を説明する図。
【発明を実施するための形態】
【０００９】
以下、添付図面を参照して本発明の好適な実施形態について説明する。
【実施例１】
【００１０】
図１は、複数の放射線検出装置を並べて行われる長尺撮影に用いられる放射線撮影システムの概略構成を示す図である。
【００１１】
放射線撮影システムは、放射線を発生させる放射線発生部１１２を備えている。放射線発生部１１２は、照射範囲１１４に放射線を照射することができる。放射線発生部１１２は、床面又は天井に設置された支持部（図示しない。）を介して設置されている。放射線発生部１１２の照射面には、放射線を遮蔽する絞り（図示しない。）が設置されている。操作者は、放射線を遮蔽する絞りを制御することにより、放射線発生部１１２から照射される放射線の照射範囲１１４を設定することができる。
【００１２】
放射線撮影システムは、複数の放射線検出装置１２０、１２２、１２４を備えている。ここでは、３つの放射線検出装置１２０、１２２、１２４を備えた形態を示すが、２つの放射線検出装置、４つ以上の放射線検出装置であってもよい。複数の放射線検出装置１２０、１２２、１２４は、被検者１００を透過した放射線を検出し、放射線に応じた画像データを出力するものである。なお、画像データを放射線画像と言い換えることもできる。
【００１３】
具体的には、複数の放射線検出装置１２０、１２２、１２４は、被検者を透過した放射線を、透過放射線量に相当する電荷として検出する。例えば、複数の放射線検出装置１２０、１２２、１２４には、放射線を電荷に変換するａ−Ｓｅなどの放射線を直接的に電荷に変換する直接変換型センサや、ＣｓＩなどのシンチレータとａ−Ｓｉなどの光電変換素子を用いた間接型センサが用いられる。さらに、複数の放射線検出装置１２０、１２２、１２４は、検出された電荷をＡ／Ｄ変換することにより、画像データを生成し、画像表示制御部１３０へ出力する。
【００１４】
複数の放射線検出装置は、撮影台１１０内に収納されている。撮影台１１０は、矩形の筐体であり、筐体内は中空である。また、撮影台１１０は、複数の放射線検出装置１２０、１２２、１２４を保持する機能を有している。図１に示すように、撮影台１１０を床面に対して直立させ、撮影台１１０が設置される。被検者１００は、撮影台１１０の長手方向に沿って設置される。撮影台１１０は、被検者１００を支える支持機能を有している。
【００１５】
図１では、撮影台１１０の長手方向が鉛直方向となるように、すなわち、撮影台１１０が床面に対して直立するように撮影台１１０が設置される。なお、撮影台１１０の長手方向が水平方向となるように、すなわち、撮影台１１０が床面に対して平行となるように撮影台１１０が設置されてもよい。
【００１６】
撮影台１１０には、放射線検出装置１２０、放射線検出装置１２２、放射線検出装置１２４が撮影台１１０の長手方向に沿ってそれぞれ配置される。このとき、放射線検出装置の一部を重ねながら複数の放射線検出装置が配置される。例えば、図１に示すように、放射線検出装置１２０と放射線検出装置１２２は、一部が空間的に互いに重なるように配置される。このとき、放射線検出装置１２０と放射線検出装置１２２の撮影可能領域は互いに重なっている。同様にして、放射線検出装置１２２と放射線検出装置１２４は、一部が空間的に互いに重なるように配置される。このとき、放射線検出装置１２２と放射線検出装置１２４の撮影可能領域は互いに重なっている。また、放射線検出装置１２２は、放射線検出装置１２０と放射線検出装置１２４の背面側、つまり放射線発生部１１２から遠い位置に配置されている。
【００１７】
また、放射線撮影システムは、放射線検出装置から出力された画像データに対して画像処理を行ない、画像を生成する画像表示制御部１３０と、画像を表示する表示部１３２と、操作者から指示を行うための操作部１３４とを備えている。また、画像表示制御部１３０は、各構成要素を制御する機能を有している。
【００１８】
画像表示制御部１３０は、複数の放射線検出装置１２０、１２２、１２４に接続されている。具体的には、画像表示制御部１３０は、複数の放射線検出装置１２０、１２２、１２４と有線または無線のネットワークもしくは専用線で接続されている。複数の放射線検出装置１２０、１２２、１２４は、放射線発生部１１２で発生した放射線を撮像し、画像データを画像表示制御部１３０に出力する。画像表示制御部１３０は、コンピュータ上で動作するアプリケーション機能を有している。画像表示制御部１３０は、複数の放射線検出装置１２０、１２２、１２４の動作を制御しつつ、表示部１３２へ画像を出力したり、グラフィカルユーザーインターフェース（図示しない。）を出力したりする。
【００１９】
画像表示制御部１３０は、放射線発生部１１２の放射線を発生するタイミングと放射線の撮影条件を制御する。また、画像表示制御部１３０は、複数の放射線検出装置１２０、１２２、１２４の画像データを撮影するタイミング及び出力するタイミングを制御する。画像表示制御部１３０は、複数の放射線検出装置１２０、１２２、１２４に対して同時に撮影を行わせ、複数の放射線検出装置１２０、１２２、１２４に対して同時に画像データを出力させることができる。
【００２０】
画像表示制御部１３０は、複数の放射線検出装置１２０、１２２、１２４から出力された画像データに対して、階調変換などの画像処理を行う機能を有している。また、画像表示制御部１３０は、複数の放射線検出装置１２０、１２２、１２４から出力された画像に対してトリミングや回転といった画像処理を行なうこともできる。表示部１３２は、画像表示制御部１３０から出力される当該画像を表示させる。
【００２１】
被検者１００は、撮影台１１０に置かれた踏み台上に立ち、複数の放射線検出装置１２０、１２２、１２４および放射線発生部１１２に対して位置決めされる。本実施例では、放射線検出装置１２２の中心に垂直に放射線が入射する角度となっている。放射線発生部１１２から複数の放射線検出装置１２０、１２２、１２４に向け照射された放射線は、被検者１００を透過して複数の放射線検出装置１２０、１２２、１２４に到達して検出される。複数の放射線検出装置１２０、１２２、１２４で得られた画像データは、画像表示制御部１３０で合成処理され、被検者１００の合成画像が生成される。合成画像は、観察領域が広い長尺撮影によって取得される合成画像である。表示部１３２は、画像表示制御部１３０から出力される合成画像を表示させる。
【００２２】
本発明の放射線撮影システムでは、１回の放射線の照射によって、被検者１００の脊椎や下肢の全体を撮影する長尺撮影を行うことができる。放射線発生部１１２から照射される放射線（照射範囲１１４）が複数の放射線検出装置１２０、１２２、１２４に同時に照射される。例えば、操作者は、放射線を遮蔽する絞りを制御したり、複数の放射線検出装置１２０、１２２、１２４と放射線発生部１１２との距離を調整したりする。
【００２３】
なお、複数の放射線検出装置１２０、１２２、１２４は、放射線発生部１１２からの放射線の照射を自動検知する検知機能を有していてもよい。自動検知する検知機能は、放射線発生部１１２から放射線が照射された際、複数の放射線検出装置１２０、１２２、１２４が放射線を検知して放射線に起因する電荷を蓄積する機能である。複数の放射線検出装置１２０、１２２、１２４のいずれかに１つより放射線の照射を検知した際、複数の放射線検出装置１２０、１２２、１２４は、本読み動作を開始させて画像データを取得する。
【００２４】
上述した放射線撮影システムでは、放射線検出装置１２０、１２４の背後に放射線検出装置１２２が重なり合うように配置されている。このため、放射線検出装置１２２が出力する画像データには、放射線検出装置１２０、１２４の内部構成要素である放射線検出パネル、基板、筐体などの構造が写り込んだ領域（以下、欠損領域と呼ぶ）が生じる。
【００２５】
この欠損領域の発生メカニズムについて図２を用いて説明する。図２（ａ）は放射線検出装置１２０、１２２の内部構造の一例を示したものである。放射線検出装置１２０は、放射線入射面側から、放射線を検出する放射線検出パネル１５０、放射線検出パネル１５０を粘着してパネル基台１５８に設置させる粘着材１５６と、放射線検出パネル１５０を支持するパネル基台１５８、放射線検出パネル１５０から電気信号を出力させる制御基板１５４の順に積層された結合体が内包される。放射線検出パネル１５０と制御基板１５４は、フレキシブル基板１５２を介して接続されている。
【００２６】
また、放射線検出装置１２０の外装筺体は、金属若しくはカーボンからなる外装筺体１６０と、放射線を透過させる放射線透過部材から成る放射線透過部１６２とから構成される。放射線検出パネル１５０の放射線入射面には、放射線透過部１６２が設置され、放射線発生部１１２からの放射線の減衰を抑える。放射線検出パネル１５０は、放射線を検出可能な有効画素領域と、該有効画素領域の外周に辺縁部を有する。
【００２７】
放射線検出装置１２２は、その有効画素領域が放射線検出装置１２０の有効画素領域と一部重なるように配置され、どのラインにおいても放射線検出装置１２０、１２２のいずれかの有効画素領域が確実に画像情報を取得するように構成される。合成画像は、放射線検出装置１２０から出力される画像データ（放射線画像）と、放射線検出装置１２２から出力される画像データのうちの放射線検出装置１２０が取得していない画像領域の画像データ（放射線画像）から生成される。
【００２８】
ここで、放射線検出装置１２２から取得される画像データ３０２には、放射線検出装置１２０の構造が写り込んでいる。放射線検出装置１２２の有効画素領域の端部から放射線検出装置１２２の外装筺体の端部までの領域４１０は、放射線検出装置１２０の構造が放射線検出装置１２２に写り込んでしまう領域である。放射線検出装置１２２から取得される画像データ３０２には、放射線検出装置１２０の構造の写り込みによる欠損領域４１２が生じる。よって、合成処理部１４２において、放射線検出装置１２２から取得される画像データ３０２から生成される合成画像にも同様に欠損領域４１２が生じる。
【００２９】
放射線検出装置１２２から取得される画像データ３０２の欠損領域４１２には、放射線検出装置１２０における放射線検出パネル１５０、フレキシブル基板１５２、粘着材１５６、パネル基台１５８、金属筺体１６０の一部が構造陰影として写り込む。また、欠損領域４１２には、フレキシブル基板１５２上の基板や、ネジなどに起因する構造陰影も含まれる。
【００３０】
以上説明した通り、欠損領域は前面に配置された放射線検出装置の内部構造によって生じる構造陰影が写り込む領域であり、この欠損領域では被検者の陰影に構造陰影が重畳し診断時の妨げとなる可能性がある。
【００３１】
この構造陰影は、放射線検出装置の内部構造に応じて大きく異なる。例えば、図２（ｂ）に示す放射線検出装置１２０は内部構造を工夫することで図２（ａ）に比べ構造陰影を単純化および軽減を実現している。
【００３２】
具体的には、図２（ｂ）に示した放射線検出装置１２０では、外装筺体１６０を金属筐体から放射線の減衰を抑えた放射線透過筐体に置き換えることで構造陰影の軽減を図っている。また、粘着材１５６、パネル基台１５８が放射線検出装置１２２の有効画素領域と重ならないように内側配置することで粘着材１５６、パネル基台１５８の写り込みを防いでいる。さらに、フレキシブル基板やネジについても有効画素領域と重ならない位置に配置（図示しない。）することで基板やネジに起因する複雑な構造陰影の写り込みを防いでいる。
【００３３】
これにより、放射線検出装置１２２から取得される画像データ３０２の欠損領域４１４は、主に放射線検出パネル１５０、外装筺体１６０から構成され、その構造陰影は単純なものとなる。また、内部構造による放射線の減衰も少ないため構造陰影のコントラストも軽減される。
【００３４】
なお、放射線検出装置の内部構造は図２に示した以外にも様々なものがある。放射線検出装置の上端部と下端部で内部構造が異なる放射線検出装置も存在する。例えば、放射線検出装置の上端部によって、図２（ａ）に示すような構造陰影が写り込む領域が画像データ３０２に含まれる場合がある。つまり、放射線検出装置１２２から取得される画像データ３０２には、放射線検出装置１２０の構造の写り込みによる欠損領域４１２が生じる。また、放射線検出装置の下端部によって、図２（ｂ）に示すような構造陰影が写り込む領域が画像データ３０２に含まれる場合がある。つまり、放射線検出装置１２２から取得される画像データ３０２には、放射線検出装置１２０の構造の写り込みによる欠損領域４１４が生じる。
【００３５】
次に、図３に示す本発明の放射線撮影システムの構成図を用いて、上述した放射線検出装置の重ね合わせに起因する合成画像の欠損領域を補正し、画質を向上させる形態を説明する。
【００３６】
画像表示制御部１３０は、放射線検出装置から出力された画像データを記憶する記憶部１４０と、画像データを合成して合成画像を生成する合成処理部１４２と、合成画像に生じる欠損領域を目立たないように補正する画像補正部１４６と、画像補正部１４６によって補正された合成画像に対して階調処理を行う階調処理部１４８とを備えている。
【００３７】
記憶部１４０は、複数の放射線検出装置１２０、１２２、１２４から出力される画像データ（放射線画像）を記憶する。図３に示すように、放射線検出装置１２０、１２２、１２４は、それぞれ、放射線検出装置（Ｄ１）、放射線検出装置（Ｄ２）、放射線検出装置（Ｄ３）としている。
【００３８】
記憶部１４０は、放射線検出装置１２０、１２２、１２４から出力される画像データを時間情報とともに記憶することができる。よって、記憶部１４０は、放射線画像が取得された時間情報によって、放射線検出装置１２０、１２２、１２４から出力された放射線画像が同時に取得されたものであるかどうかを区別して記憶することができる。記憶部１４０は、被検者の画像情報が含まれた放射線画像であるのか、被検者の画像情報が含まれていない放射線画像であるのか区別して記憶することができる。
【００３９】
また記憶部１４０は、複数の放射線検出装置１２０、１２２、１２４によって同時に撮影された複数の放射線画像を、放射線検出装置の位置情報（空間的配置情報）と関連付けて記憶することができる。例えば、記憶部１４０は、放射線検出装置１２０から出力される画像データと放射線検出装置１２２から出力される画像データとが隣接していることを関連付けて記憶することができる。同様にして、記憶部１４０は、放射線検出装置１２２から出力される画像データと放射線検出装置１２４から出力される画像データとが隣接していることを関連付けて記憶することができる。さらに記憶部１４０は、放射線検出装置１２２が放射線検出装置１２０、１２４の背面側に配置されていることを関連付けて記憶することができる。記憶部１４０は、合成処理部１４２に対して、複数の画像データとその位置情報を出力することができる。
【００４０】
また、記憶部１４０は、放射線検出装置１２０、１２２、１２４を識別する識別情報（識別コード）を記憶することができる。識別コードは、放射線検出装置の種類を示す情報を含んでいる。言い換えれば、識別コードは、放射線検出装置の内部構造を示す情報であり、放射線検出装置の構造陰影（構造の写り込み）に対応する情報を含んでいる。また、識別コードは、放射線検出装置の種類に加えて、放射線検出装置の位置情報（配置情報）による情報を含んでいてもよい。放射線検出装置の位置情報（配置情報）は、例えば、放射線検出装置が撮影台１１０における位置を示す情報、放射線検出装置の上端部及び下端部を示す情報である。放射線検出装置の位置情報（配置情報）は、放射線検出装置の上下の向きの情報を含んでいてもよい。このように、識別コードによって、背面側に配置されている放射線検出装置から取得される画像データに写り込む構造陰影を識別することができる。
【００４１】
このように、識別コードによって、例えば、図２（ａ）に示す放射線検出装置１２０と、図２（ｂ）に示す放射線検出装置１２０とを識別することができる。識別コードによって、図２（ａ）に示すような構造陰影が写り込む領域が画像データ３０２に含まれるか、図２（ｂ）に示すような構造陰影が写り込む領域が画像データ３０２に含まれるかを識別することができる。
【００４２】
合成処理部１４２は、記憶部１４０に記憶された複数の画像データを合成して、合成画像を生成する。このとき、合成処理部１４２は、被検者１００の画像情報が含まれた複数の画像データについて合成して、合成画像を生成する。
【００４３】
合成処理部１４２は、放射線検出装置１２０、１２２、１２４から出力された複数の画像データとその時間情報及び位置情報に基づいて合成することにより、合成画像を生成する。具体的には、合成処理部１４２は、放射線検出装置１２０、１２２、１２４から時間情報に基づいて同時に出力された複数の画像データ（放射線画像）を合成対象と判別し、複数の画像データを合成する。合成処理部１４２は、位置情報に基づいて放射線検出装置１２０、１２２、１２４から出力された複数の画像データの位置関係を決定して合成する。ここでは、合成画像の欠損領域に存在する構造陰影を識別する識別コードが合成画像の付帯情報として記憶される。
【００４４】
例えば、図１に示す例では、放射線検出装置１２０から出力された画像データが上方に、放射線検出装置１２４から出力された画像データが下方に、放射線検出装置１２２から出力された画像データがその間に位置決めされる。さらに位置情報が示す重なり方も考慮して合成が行われる。例えば、放射線発生部１１２から遠い位置に他の放射線検出装置に重なり合うように配置された放射線検出装置１２２には、上下に欠損領域が生じる。放射線検出装置１２０、１２４には欠損領域は生じない。そこで、合成処理部１４２は、放射線検出装置が重なり合う範囲では放射線検出装置１２０、１２４が生成する画像データを用いて合成画像を生成することで合成画像に生じる欠損領域の面積を最小化することができる。このように、合成処理部１４２は、隣接する複数の撮影領域を撮影して得た複数の画像データを合成することにより、合成画像を生成することができる。
【００４５】
また、画像データの上下に生じた欠損領域の構造陰影を識別する識別コードが付帯情報として合成画像とともに記憶される。例えば、合成処理部１４２（画像表示制御部１３０）は、識別コードを合成画像とともに記憶する。例えば、図１に示した放射線検出装置１２０の内部構造が図２（ａ）で示した放射線検出装置１２０と同様のものであるならば、合成画像に存在する欠損領域に写り込む構造陰影は図２（ａ）の４１２に示したものとなる。そこで、この構造陰影を示すユニークなＩＤが付帯情報として記憶される。また、図１に示した放射線検出装置１２４の内部構造が図２（ｂ）で示した放射線検出装置１２０と同様のものであるならば、合成画像に存在する欠損領域に写り込む構造陰影は図２（ｂ）の４１４に示したものとなる。
【００４６】
画像補正部１４６は、合成処理部１４２から出力された合成画像に存在する欠損領域を構造情報に基づき補正する。ここで、構造情報とは、欠損領域に写り込む構造陰影を表す情報である。また、本実施例ではさらに構造情報として被検者が無い状態で複数の放射線検出装置を重ね合わせて撮影した構造データを利用する。なお、この構造データは放射線検出装置の内部構造の写り込みを画像として保持したデータであり、この画像の画素値は例えば放射線源弱係数が大きく厚い構造による写り込みが生じている画素では小さな値を、放射線源弱係数が小さく薄い構造による写り込みが生じている画素では大きな値となる。
【００４７】
なお、構造データは事前に取得され、記憶部１４０に記憶しておく。例えば、図１に示す放射線検出装置の配置で被検者を撮影する場合は、図４で示すように放射線検出装置１２０、１２２、１２４の配置を図１と同じようにし、被検者が無い状態で撮影を行う。このように撮影された放射線検出装置１２２から取得される画像データ３０２には、重複する放射線検出装置１２０の下端部おける内部構造の写り込み領域３０６が含まれている。また、放射線検出装置１２２から取得される画像データ３０２には、重複する放射線検出装置１２４の上端部における内部構造の写り込み領域３０８が含まれている。
【００４８】
ここで、放射線検出装置１２０から取得される画像データ（放射線画像）３００には、他の放射線検出装置の構造情報の写り込みは生じない。また、放射線検出装置１２４から取得される画像データ（放射線画像）３０４には、他の放射線検出装置の構造情報の写り込みは生じない。よって、画像データ３０２が、画像に写り込む構造陰影を表す構造データに相当し、この画像データ３０２を構造データとして記憶部１４０に保持しておけばよい。なお、上端部の写り込み領域３０６及び下端部の写り込み領域３０８を別々の構造データと見なし保持してもよい。
【００４９】
階調処理部１４８は、複数の画像データ（放射線画像）を合成して得られた合成画像に対して、階調処理を行なう。具体的には、階調処理部１４８は、放射線検出装置１２０、１２２、１２４から取得された複数の画像データを記憶部１４０から取得する。階調処理部１４８は、放射線検出装置１２０、１２２、１２４から取得された複数の画像データの特徴量をそれぞれ解析して、表示部１３２のダイナミックレンジを有効に利用することができるように、合成画像の階調変換特性を決定する。
【００５０】
そして、階調処理部１４８は、決定された階調変換特性を用いて合成画像の階調を変換する。特徴量には、各画像データのヒストグラム、最大画素値、最小画素値が含まれ、放射線検出装置１２０、１２２、１２４から取得された複数の画像データに対して解析処理を実行することにより、特徴量を算出している。
【００５１】
階調処理部１４８は、画像補正部１４６によって補正が行われた合成画像に対して、階調処理を行う。ここでは、欠損領域が低減された合成画像に対して階調処理を行うため、合成画像の階調処理を適切に行うことができる。つまり、階調処理部１４８は、放射線検出装置１２０及び放射線検出装置１２４の構造物の写り込みの影響を抑えて、合成画像の階調処理を行うことができる。表示部１３２は、階調処理部１４８によって階調処理がなされた合成画像を表示する。
【００５２】
次に、本実施の放射線撮影システムの特徴である画像補正部１４６の動作手順について図５のフローチャートを用いて詳細に説明する。
【００５３】
上述の如く、合成処理部１４２は、位置情報に基づいて放射線検出装置１２０、１２２、１２４から出力された複数の画像データの位置関係を決定して合成する。また、合成画像の欠損領域に存在する構造陰影を識別する識別情報（識別コード）が合成画像の付帯情報として記憶される。
【００５４】
画像補正部１４６は、合成処理部１４２から合成画像、および構造情報として欠損領域に写り込む構造陰影に関する識別コードを取得する（Ｓ５０１）。本実施例における識別コードは、放射線検出装置の種類と、放射線検出装置の位置情報（配置情報）による情報を含んでいるため、合成画像における欠損領域の形態（構造陰影）を特定することができる。
【００５５】
ここでは、本実施例における放射線検出装置１２０と放射線検出装置１２４の内部構造が異なるものとする。具体的には、放射線検出装置１２０は図２（ｂ）で示した放射線検出装置１２０と同様の内部構造を持つものとする。また、放射線検出装置１２４は図２（ａ）で示した放射線検出装置１２０と同様の内部構造とする。
【００５６】
よって、本実施例における合成画像には２つの異なる構造陰影が存在することとなる。具体的には、合成画像の上部に存在する欠損領域６０１には図６（ｃ）に示したような単純でかつ放射線の減衰が少ない構造が写り込む。また、合成画像の下部に存在する欠損領域６０３は図６（ｄ）に示したような複雑でかつ放射線の減衰が多い構造が写り込む。なお、以降の説明では図６（ｃ）に示した構造陰影の識別コードを０ｘ００とし、図６（ｄ）に示した構造陰影の識別コードを０ｘ０１とし説明を行う。
【００５７】
次に、画像補正部１４６は、合成画像における補正対象である欠損領域の情報を取得する（Ｓ５０２）。具体的には、画像補正部１４６は、図７に示すように合成画像上の欠損領域の開始行Ｙ０と終了行Ｙ１を欠損領域の位置情報として取得する。なお、図６（ａ）に示すように合成画像上に欠損領域が複数存在する場合は、補正が実行されていない欠損領域のうち、何れか１つの情報を任意に取得する。例えば、欠損領域６０１と欠損領域６０３の何れも補正が実行されていない場合は、どちらか一方の欠損領域の情報を任意に取得する。欠損領域６０１の補正が実行されており、欠損領域６０３の補正が実行されていない場合は、欠損領域６０３の情報を取得する。
【００５８】
さらに、画像補正部１４６は、欠損領域の情報として構造データ上の欠損領域の開始行Ｙ２と終了行Ｙ３および構造陰影の識別コードを合わせて取得する。例えば、図６の欠損領域６０１の情報を取得する場合は欠損領域６０１の開始行Ｙ０と終了行Ｙ１を取得するとともに、対応する構造データ上の欠損領域６０２の開始行Ｙ２および終了行Ｙ３を取得する。また、画像補正部１４６は、構造陰影の識別コードとして０ｘ００を取得する。
【００５９】
なお、以降の説明では補正前の合成画像をＩ、構造データをＰで表す。また、開始行Ｙ０、Ｙ２の座標を１、終了行Ｙ１、Ｙ３の座標をＨとし、ｙ番目（１≦ｙ≦Ｈ）の欠損行におけるｘ番目（１≦ｘ≦Ｗ）の座標を（ｘ，ｙ）、その座標における補正前の画素値をＩ（ｘ，ｙ）、構造データの画素値をＰ（ｘ，ｙ）と表すこととする。
【００６０】
次に、画像補正部１４６は、取得した構造陰影の識別コードをチェックし（Ｓ５０３）、識別コードに応じて補正処理を実行する。言い換えれば、画像補正部１４６は、識別コードに応じて補正方法を設定する。識別コードが０ｘ００と識別コードが０ｘ０１の場合では、画像補正部１４６による補正方法が異なる。画像補正部１４６は、識別コードが０ｘ００の場合、合成画像に対して低次の補正を行うのに対し、識別コードが０ｘ０１の場合、合成画像に対して高次の補正を行う。低次の補正とは、構造陰影が単純な画像の補正に適した補正であり、高次の補正とは、構造陰影が複雑な画像の補正に適した補正である。補正方法を区別するのは、構造陰影が単純な画像に対して、高次の補正をかけてしまうと、不必要な補正をかけてしまい、画質が劣化する可能性があるからである。例えば、低次の補正は、１次式による補正、１段階の補正などの少なくとも１つが含まれ、高次の補正は、多項式による補正、複数段階の補正などの少なくとも１つが含まれる。つまり、画像補正部１４６は、構造陰影が単純な合成画像の補正に適した補正、構造陰影が複雑な合成画像の補正に適した補正の一方を行う。
【００６１】
具体的には、識別コードが０ｘ００の場合、画像補正部１４６は、図６（ｃ）に示した単純でかつ放射線の減衰が少ない構造陰影に適した補正処理を実行する（Ｓ５０４）。具体的には、画像補正部１４６は、構造陰影が無い画像Ｏ（すなわち、理想的な補正後の画像）と補正前の画像Ｉおよび構造データＰの関係が下記式でモデル化できると仮定し補正を行う。
０（ｘ，ｙ）＝Ｉ（ｘ，ｙ）−（ａ・Ｐ（ｘ，ｙ）＋ｂ）    （数１）
【００６２】
上式において、ａ、ｂはモデルパラメータであり、画像Ｉには構造データＰの一次式で表される構造陰影が加法的に重畳していると仮定したモデルとなっている。なお、このモデルパラメータａ、ｂは未知であり、最適なモデルパラメータを推定する必要がある。そこで、本実施例では最小二乗基準の最適化問題を解くことで最適なモデルパラメータを推定する。具体的には下記式で表される誤差Ｊが最小となるモデルパラメータを算出する。
【００６３】
【数１】
2019068915.tif 000003
【００６４】
ここで、上式は各行のＷ個の画素サンプルから行毎にモデルパラメータを推定するものである。また、上式において構造陰影の無い画像Ｏは未知である。そこで、隣接する行同士の相関が高いことを利用し、構造陰影の無い隣接行のデータをＯの推定値として用いる。
【００６５】
なお、構造陰影の無い隣接行は欠損領域の開始行または終了行のみに存在する。そこで、画像補正部１４６は、モデルパラメータの推定および補正は開始行または終了行から１行ずつ補正し、補正後の行を次の行の構造陰影の無い隣接行として逐次的に補正を行っていく。
【００６６】
具体的には、画像補正部１４６は、開始行から逐次的に補正を行う場合は開始行（ｙ＝１）から、次式によるモデルパラメータの推定と補正を下の行に向かって１行ずつ順番に終了行（ｙ＝Ｈ）まで繰り返す。
【００６７】
【数２】
2019068915.tif 000004
【００６８】
また、画像補正部１４６は、終了行から逐次的に補正を行う場合は終了行（ｙ＝Ｈ）から、次式によるモデルパラメータの推定と補正を上の行に向かって１行ずつ順番に開始行（ｙ＝１）まで繰り返す。
【００６９】
【数３】
2019068915.tif 000005
【００７０】
なお、上述の補正のうち１方向の補正結果のみを用いてもよいが、上下２方向の補正結果をブレンドすることもできる。例えば、上方向から補正した結果をＯ１、下方向から補正した結果をＯ２とすれば、次式のように補正の開始行から距離に基づく重みを考慮して補正結果をブレンドすればよい。なお、この場合では補正の開始行から次の行に伝搬する誤差を抑えることができ、１方向の補正結果のみを用いるよりも高精度な補正が行える。
【００７１】
【数４】
2019068915.tif 000006
【００７２】
以上、識別コードが０ｘ００の場合、すなわち図６（ｃ）に示した構造陰影に対する補正について説明した。図６（ｃ）に示した構造陰影は単純でかつｘ方向に均一な構造であるため、一次式のような低次のモデルで十分な補正精度を実現できる。また、ｘ方向に均一な構造であるため行単位で唯一のモデルパラメータを推定すれば良く処理速度の面においても有利である。なお、単純な構造に対し高次のモデルを適用した場合や、画素単位でモデルパラメータを推定した場合ではオーバーフィッティングによる被検者の陰影劣化が起きやすいため、低次のモデルが好適である。
【００７３】
次に、識別コードが０ｘ０１の場合は、画像補正部１４６は、図６（ｄ）に示した複雑でかつ放射線の減衰が多い構造陰影に適した補正処理を実行する（Ｓ５０５）。具体的には、画像補正部１４６は、構造陰影が無い画像Ｏ（すなわち、理想的な補正後の画像）と補正前の画像Ｉおよび構造データＰの関係が下記式でモデル化できると仮定し補正を行う。
【００７４】
【数５】
2019068915.tif 000007
【００７５】
上式において、ａｉ，ｊ、ｂはモデルパラメータであり、複雑な構造陰影に適合させるために識別コード０ｘ００の場合に比べモデルパラメータの数を増やし、モデルの自由度を大きくしている。また、このモデルでは構造データＰの周辺画素を考慮することで画像Ｉと構造データＰの僅かな幾何学的ずれによる補正誤差を抑制することができる。これにより構造陰影内に多数存在するエッジの僅かな位置ずれにより生じる補正残りを抑えることができる。
【００７６】
なお、画像補正部１４６は、このモデルパラメータａｉ，ｊ、ｂは最小二乗基準の最適化問題を解くことで推定する。具体的には下記式で表される誤差Ｊが最小となるモデルパラメータを算出する。
【００７７】
【数６】
2019068915.tif 000008
【００７８】
ここで、上式においてＯの推定値として、識別コード０ｘ００と同様に補正後の隣接行データを用いるが、放射線の減衰が多い構造陰影の場合、補正後データに重畳するノイズの影響で推定精度が悪くなる可能性がある。そこで、補正後の隣接行データを平均したものを推定値として用いる。
【００７９】
具体的には、画像補正部１４６は、開始行から逐次的に補正を行う場合は開始行（ｙ＝１）から、次式によるモデルパラメータの推定と補正を下の行に向かって１行ずつ順番に終了行（ｙ＝Ｈ）まで繰り返す。
【００８０】
【数７】
2019068915.tif 000009
【００８１】
ここで、Ｎは平均する補正後の隣接行データの数を表し、例えばＮ＝５を用いる。
【００８２】
また、画像補正部１４６は、終了行から逐次的に補正を行う場合は終了行（ｙ＝Ｈ）から、次式によるモデルパラメータの推定と補正を上の行に向かって１行ずつ順番に開始行（ｙ＝１）まで繰り返す。
【００８３】
【数８】
2019068915.tif 000010
【００８４】
ここで、Ｎは平均する補正後の隣接行データの数を表し、例えばＮ＝５を用いる。
【００８５】
なお、画像補正部１４６は、上述の補正のうち１方向の補正結果のみを用いてもよいが、上下２方向の補正結果をブレンドすることもできる。この補正については、識別コード０ｘ００と同様であるため、説明は省略する。
【００８６】
また、画像補正部１４６は、上述では行単位で唯一のモデルパラメータを求めたが、行を区間に分けて区間毎にモデルパラメータを推定し、補正を行ってもよい。例えば図７に示すように行を区間Ａ〜Ｋに分けて各区間に対し上述したモデルパラメータの補正および推定を行う。また、画素毎にモデルパラメータを求めてもよい。例えば、座標（ｘ，ｙ）におけるモデルパラメータを推定する場合は、次式のように同じ行の近傍の数画素のサンプルにおける誤差Ｊが最小となるモデルパラメータを推定し、補正を行う。
【００８７】
【数９】
2019068915.tif 000011
【００８８】
ここで、Ｎはモデルパラメータを推定する際のサンプル数を規定するパラメータであり、モデルの適合度に応じて任意に設定すればよい。
【００８９】
なお、サンプル数が少ない場合は異常値の影響を受けやすい。そこで、Ｍ推定やＬＭｅｄＳ推定などのロバストな推定方法を用いてモデルパラメータを推定してもよい。
【００９０】
以上、区間または画素毎にモデルパラメータを推定する場合を説明したが、この場合では、場所依存で起こりうる誤差の影響を低減できる。そのため、場所依存で構造が大きく変化するような構造でもモデルを適合しやすく、特に複雑な構造陰影に対して好適な補正方法である。
【００９１】
次に、識別コードが０ｘ０１の場合は、画像補正部１４６は、欠損領域に対してさらなる画像補正を行う。上述の通り図６（ｄ）に示した複雑な構造陰影では放射線の減衰が多いため、ノイズによる粒状性悪化が起きる。そこで、画像補正部１４６は、欠損領域に対し粒状性を改善する補正を行う。つまり、画像補正部１４６は、合成画像に写り込んだ放射線検出装置の構造陰影に応じてノイズの補正を実行するか否かを設定することができる。具体的には特に限定するものではないが、例えば移動平均やガウシアン等の平滑化フィルタや、εフィルタやバイラテラルフィルタ等のエッジを保存する非線形なフィルタを用いたフィルタ処理を用いれば良く、詳細は公知であるため省略する。
【００９２】
次に、画像補正部１４６は、すべての欠損領域の補正が完了したかをチェックし（Ｓ５０６）、補正が完了した場合は終了する。補正が完了していない場合Ｓ５０２に戻り、画像補正部１４６は、新たに設定した欠損領域に対する補正を実行する。
【００９３】
以上、実施例１では、放射線を検出する複数の放射線検出装置と、複数の放射線検出装置から取得される複数の放射線画像を合成して合成画像を生成する合成処理部１４２を有した放射線撮影システムにおいて、合成画像における放射線検出装置の構造が写り込んだ領域を補正する画像補正部１４６を備え、画像補正部１４６は写り込んだ放射線検出装置の構造陰影の特性（識別コード）に応じて補正方法を設定する。このように、放射線検出装置の構造陰影に応じて補正方法を切り替えることで、構造陰影に応じた好適な補正を行うことができる。
【００９４】
なお、本実施例では、合成処理部１４２と画像補正部１４６を別の構成で示しているが、放射線画像の合成と同時に放射線検出装置の構造が写り込んだ領域の補正を行ってもよい。
【実施例２】
【００９５】
次に実施例２について、図８、図９を用いて説明する。実施例１と異なる点は、合成画像において放射線検出装置の構造陰影が写る欠損領域を細分化し、領域毎に補正方法を切り替える点である。ここでは、画像補正部１４６は、合成画像における構造が写り込んだ領域を構造陰影の違いに応じて分割し、分割された領域毎に補正方法を設定する。具体的には、例えば図５のＳ５０５の動作を図８のフローチャートに従った動作とする。以下より図８のフローチャートを用いて説明する。
【００９６】
画像補正部１４６は、実施例１で説明した識別コード０ｘ０１に対する補正処理において、欠損領域が写り込む構造陰影に応じて分割を行う（Ｓ８０１）。本実施例では、画像補正部１４６は、例えば構造陰影の形状の複雑さと放射線検出装置の構造による放射線の減衰に応じた分割を行う。具体的には、画像補正部１４６は、図９に示したように構造が単純でかつ放射線の減衰が少ない領域Ａ、構造が複雑でかつ放射線の減衰が多い領域Ｂ、構造が単純でかつ放射線の減衰が多い領域Ｃに分割する。なお、どのように領域を分割するかは放射線検出装置の内部構造に応じて識別コードとともに予め設定することもできる。
【００９７】
次に、画像補正部１４６は、分割した領域のうち補正が実行されていない領域を選択する（Ｓ８０２）。画像補正部１４６は、構造が単純でかつ放射線の減衰が少ない領域Ａ、構造が複雑でかつ放射線の減衰が多い領域Ｂ、構造が単純でかつ放射線の減衰が多い領域Ｃを判断する（Ｓ８０３）。そして、判断された領域の特性に応じた補正処理を実行する（Ｓ８０４〜Ｓ８０５）。
【００９８】
構造が単純でかつ放射線の減衰が少ない領域Ａでは、画像補正部１４６は、実施例１のＳ５０４で説明した補正と同じ補正を行う（Ｓ８０４）。また、構造が複雑でかつ放射線の減衰が多い領域Ｂでは、画像補正部１４６は、実施例１のＳ５０５で説明した補正と同じ補正を行う（Ｓ８０５）。すなわち、画像補正部１４６は、自由度の大きいモデルを用いた補正を行った後、さらにノイズの補正を行う。ノイズの補正とは、実施例１で説明した粒状性を改善する補正等である。次に構造が単純でかつ放射線の減衰が多い領域Ｃでは、画像補正部１４６は、実施例１のＳ５０４で説明した補正と同じ補正に加え、さらにノイズの補正を行う（Ｓ８０６）。
【００９９】
次に、画像補正部１４６は、すべての領域の補正が完了したかをチェックする（Ｓ８０７）。補正が完了した場合は終了する。補正が完了していない場合Ｓ８０２に戻り、画像補正部１４６は、未処理の領域を選択し補正を実行する。
【０１００】
以上、実施例２では構造陰影を細分化し補正方法を切り替えることで、構造陰影の変化が大きい場合であっても好適な補正を行うことができる。
【０１０１】
以上、本発明の好ましい実施例について説明したが、本発明はこれらの実施例に限定されないことはいうまでもなく、その要旨の範囲内で種々の変形及び変更が可能である。
【符号の説明】
【０１０２】
１００  被検者
１１０  撮影台
１１２  放射線発生部
１１４  照射範囲
１２０  放射線検出装置（Ｄ１）
１２２  放射線検出装置（Ｄ２）
１２４  放射線検出装置（Ｄ３）
１３０  画像表示制御部
１３２  表示部
１３４  操作部
１４０  記憶部
１４２  合成処理部
１４６  画像補正部
１４８  階調処理部
【図１】
2019068915.tif 000012
【図２】
2019068915.tif 000013
【図３】
2019068915.tif 000014
【図４】
2019068915.tif 000015
【図５】
2019068915.tif 000016
【図６】
2019068915.tif 000017
【図７】
2019068915.tif 000018
【図８】
2019068915.tif 000019
【図９】
2019068915.tif 000020