(57)【要約】 【課題】操作者が意図した基準でより適切な線量管理を行えるようにする。【解決手段】画像処理装置は、被検体に対する放射線撮影により得られた放射線画像を複数の解剖学的領域に分割し、該複数の解剖学的領域のうちの少なくとも一つの領域を抽出し、該抽出された領域における画素値に基づいて、該抽出された領域に対する該放射線撮影の線量指標値を算出する。【選択図】図１
2020025730.tif 000002
【特許請求の範囲】
【請求項１】
画像処理装置であって、
被検体に対する放射線撮影により得られた放射線画像を複数の解剖学的領域に分割する分割手段と、
前記複数の解剖学的領域のうちの少なくとも一つの領域を抽出する抽出手段と、
前記抽出手段により抽出された領域における画素値に基づいて、前記抽出された領域に対する前記放射線撮影の線量指標値を算出する算出手段と、
を有することを特徴とする画像処理装置。
【請求項２】
前記分割手段は、前記複数の解剖学的領域に異なるラベル番号を付与することにより、放射線画像を複数の解剖学的領域に分割することを特徴とする請求項１に記載の画像処理装置。
【請求項３】
前記分割手段は、予め機械学習によって生成されたパラメータを用いて、前記放射線画像を複数の解剖学的領域に分割することを特徴とする請求項１または２に記載の画像処理装置。
【請求項４】
機械学習により生成されたパラメータを記憶する記憶手段を更に有することを特徴とする請求項３に記載の画像処理装置。
【請求項５】
前記記憶手段は、機械学習されたパラメータを部位と関連付けて記憶することを特徴とする請求項４に記載の画像処理装置。
【請求項６】
前記抽出手段は、予め設定された、複数の部位と当該複数の部位のそれぞれに対応するラベル番号との対応情報を用いて、操作者の指示に従って、前記複数の解剖学的領域のうちの少なくとも一つの領域を抽出することを特徴とする請求項３から５のいずれか１項に記載の画像処理装置。
【請求項７】
予め設定された、所定の画像データと当該所定の画像データに対応する正解の分割割り当てデータに基づいて機械学習を行うことにより前記パラメータを更新する機械学習手段と更に有し、
前記分割手段は、前記更新されたパラメータを用いて、前記放射線画像を複数の解剖学的領域に分割することを特徴とする請求項６に記載の画像処理装置。
【請求項８】
前記分割手段により前記更新されたパラメータを用いて分割された複数の解剖学的領域に応じて、前記対応情報を更新して設定する設定手段を更に有することを特徴とする請求項７に記載の画像処理装置。
【請求項９】
前記パラメータが更新された場合、または、前記対応情報が更新された場合に、当該更新の前後に前記算出手段により算出されたそれぞれの前記線量指標値を用いて、前記線量指標値の目標値としての線量目標値を更新する更新手段を更に有することを特徴とする請求項８に記載の画像処理装置。
【請求項１０】
前記機械学習は、ＣＮＮ（Convolutional Neural Network）、ＦＣＮ（Fully Convolutional Networks）、ＳｅｇＮｅｔ、Ｕ−ｎｅｔのいずれかを用いた機械学習であることを特徴とする請求項３から９のいずれか１項に記載の画像処理装置。
【請求項１１】
前記算出手段は、前記放射線画像における、前記抽出手段により抽出された領域の中心傾向を示す値を代表値として算出し、前記代表値を用いて前記線量指標値を算出することを特徴とする請求項１から１０のいずれか１項に記載の画像処理装置。
【請求項１２】
前記抽出手段により、前記複数の解剖学的領域のうち複数の領域が抽出された場合、
前記算出手段は、前記抽出手段により抽出された複数の領域のそれぞれに対して前記代表値を算出し、複数の前記代表値から、前記抽出された複数の領域に対する複数の前記放射線撮影の線量指標値を算出することを特徴とする請求項１１に記載の画像処理装置。
【請求項１３】
画像処理方法であって、
被検体に対する放射線撮影により得られた放射線画像を複数の解剖学的領域に分割する分割工程と、
前記複数の解剖学的領域のうちの少なくとも一つの領域を抽出する抽出工程と、
前記抽出工程において抽出された領域における画素値に基づいて、前記抽出された領域に対する前記放射線撮影の線量指標値を算出する算出工程と、
を備えることを特徴とする画像処理方法。
【請求項１４】
コンピュータを、請求項１から１２のいずれか１項に記載の画像処理装置として機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、放射線撮影により得られた画像から線量指標値を出力する画像処理技術に関する。
【背景技術】
【０００２】
近年、医療分野ではデジタル画像の利用が進んでいる。デジタル画像を生成するために、放射線（Ｘ線）を間接的または直接的に電気信号に変換するセンサを用いたデジタルラジオグラフィ装置（放射線撮影装置）を用いることが主流となっている。
【０００３】
この放射線撮影装置は、放射線の線量に対して極めて広いダイナミックレンジを有しており、また画像処理による自動濃度補正によって、従来のアナログ放射線撮影に比べ線量不足や線量過多の状態でも安定した濃度の出力が得られるというメリットがある。その一方で、撮影技師が不適切な線量で撮影した場合も、撮影技師はそのことに気づきにくく、特に、線量過多の場合では、患者への被曝量が増加してしまうという課題がある。
【０００４】
そこで、この課題を解決するためにデジタル放射線画像の撮影線量の目安となる値（以下、線量指標値と呼ぶ）を撮影画像とともに表示することが一般的に行われている。また、線量指標値の算出方法についても様々な方法が提案されている。近年では、ＩＥＣ（国際電気標準会議）により国際規格ＩＥＣ６２４９４−１が発行され、標準化された線量指標値としてＥＩ（Exposure Index）が定義された。また、この国際規格では、目標とすべき線量の値（以下、線量目標値と呼ぶ）として、ＥＩＴ（Target Exposure Index）を定め、線量指標値ＥＩと線量目標値ＥＩＴとのずれ量を表すＤＩ（Deviation Index）を用いた線量管理の運用方法も与えられている。各メーカーは、この国際規格に準拠した線量管理機能を提供している。例えば特許文献１や特許文献２の方法など、メーカーから様々な算出方法が提案されている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１４−１５８５８０号公報
【特許文献２】特開２０１５−２１３５４６号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
線量指標値ＥＩの算出方法についてはブラックボックスの場合が多い。そのため、操作者にとっては線量指標値ＥＩの数値的意味合いが明確ではなく、線量管理の基準値として使い勝手の悪いものとなっていた。
【０００７】
そこで本発明は、操作者が意図した基準でより適切な線量管理を行うことを目的とする。
【課題を解決するための手段】
【０００８】
上記目的を達成するための一手段として、本発明による画像処理装置は以下の構成を有する。すなわち、画像処理装置であって、被検体に対する放射線撮影により得られた放射線画像を複数の解剖学的領域に分割する分割手段と、前記複数の解剖学的領域のうちの少なくとも一つの領域を抽出する抽出手段と、前記抽出手段により抽出された領域における画素値に基づいて、前記抽出された領域に対する前記放射線撮影の線量指標値を算出する算出手段と、を有する。
【発明の効果】
【０００９】
本発明によれば、操作者が意図した基準でより適切な線量管理を行うことが可能となる。
【図面の簡単な説明】
【００１０】
【図１】実施形態１による放射線撮影装置の構成例。
【図２】実施形態１による放射線撮影装置１００の処理手順を示すフローチャート。
【図３】実施形態１による線量指標値を算出するための領域を変更する処理手順を示すフローチャート。
【図４】実施形態１による放射線撮影装置の構成例。
【図５】実施形態２による線量指標値を算出するための領域の分割構成を変更する処理手順を示すフローチャート。
【図６】実施形態３による放射線撮影装置の構成例。
【図７】実施形態３による線量目標値ＥＩＴを自動で更新する処理手順を示すフローチャート。
【図８】（ａ）実施形態１によるセグメンテーションマップの一例を示す図、（ｂ）実施形態２による変更した正解のセグメンテーションマップの一例を示す図。
【図９】複数の撮影部位とラベル番号との対応表の一例を示す図。
【発明を実施するための形態】
【００１１】
以下、添付の図面を参照して、本発明をその実施形態の一例に基づいて詳細に説明する。なお、以下の実施形態において示す構成は一例に過ぎず、本発明は図示された構成に限定されるものではない。
【００１２】
［実施形態１］
（放射線撮影装置の構成）
実施形態１による放射線撮影装置１００の構成例を図１に示す。放射線撮影装置１００は、放射線撮影により得られた画像から線量指標値を出力する画像処理機能を有する放射線撮影装置である。すなわち、放射線撮影装置１００は、画像処理装置としても機能する。放射線撮影装置１００は、放射線発生部１０１、放射線検出器１０４、データ収集部１０５、前処理部１０６、ＣＰＵ（Central Processing Unit）１０８、記憶部１０９、操作部１１０、表示部１１１、画像処理部１１２を備えており、これらはＣＰＵバス１０７を介して互いにデータ授受が可能に接続されている。また、画像処理部１１２は、分割部１１３、抽出部１１４、算出部１１５、設定部１１６を備えている。記憶部１０９は、ＣＰＵ１０８での処理に必要な各種のデータや制御プログラムを記憶すると共に、ＣＰＵ１０８のワーキング・メモリとして機能する。ＣＰＵ１０８は、記憶部１０９を用いて、操作部１１０からの操作に従った装置全体の動作制御等を行う。これにより放射線撮影装置１００は、以下のように動作する。
【００１３】
まず、操作部１１０を介して操作者から撮影指示が入力されると、この撮影指示はＣＰＵ１０８によりデータ収集部１０５を介して放射線発生部１０１および放射線検出器１０４に伝えられる。続いて、ＣＰＵ１０８は、放射線発生部１０１及び放射線検出器１０４を制御して放射線撮影を実行させる。放射線撮影では、まず放射線発生部１０１が、被検体１０３に対して放射線ビーム１０２を照射する。放射線発生部１０１から照射された放射線ビーム１０２は、被検体１０３を透過して、放射線検出器１０４に到達する。そして、放射線検出器１０４は、到達した放射線ビーム１０２の強度（放射線強度）に応じた信号を出力する。なお、本実施形態では被検体１０３を人体とする。よって、放射線検出器１０４から出力される信号は、人体を撮影したデータとなる。
【００１４】
データ収集部１０５は、放射線検出器１０４から出力された信号を所定のデジタル信号に変換して画像データとして前処理部１０６に供給する。前処理部１０６は、データ収集部１０５から供給された画像データに対して、オフセット補正やゲイン補正等の前処理を行う。前処理部１０６で前処理が行われた画像データは、ＣＰＵ１０８の制御により、ＣＰＵバス１０７を介して、記憶部１０９、画像処理部１１２に順次転送される。
【００１５】
画像処理部１１２は、前処理部１０６から得られた画像データ（以下、放射線画像）から線量指標値を算出するための処理を行う。線量指標値とは、上述したように撮影線量の目安となる値である。分割部１１３は、入力された放射線画像（被検体１０３に対する放射線撮影により得られた放射線画像）を複数の解剖学的領域に分割する。本実施形態では、分割部１１３は、後述するセグメンテーションマップ（多値画像）を作成する。抽出部１１４は、操作者による操作部１１０への操作に基づいて、分割部１１３により分割された複数の解剖学的領域のうちの少なくとも一つの領域を、線量指標値を算出するための領域として、抽出する。算出部１１５は、抽出部１１４により抽出された領域における画素値に基づいて、抽出された領域に対する放射線撮影の線量指標値を算出する。設定部１１６は、後述する教師データや撮影部位とラベル番号との対応情報の設定や管理を行う。
【００１６】
画像処理部１１２により算出された線量指標値は、ＣＰＵ１０８の制御により、前処理部１０６から得られた放射線画像とともに表示部１１１にて表示される。操作者が表示部１１１に表示された線量指標値および放射線画像を確認後、一連の撮影動作が終了する。表示部１１１に表示された線量指標値および放射線画像は、図示しないプリンタ等に出力されても良い。
【００１７】
（放射線撮影装置の動作）
続いて、本実施形態による放射線撮影装置１００の動作に関して、図２のフローチャートを用いて具体的に説明する。図２は、本実施形態による放射線撮影装置１００の処理手順を示すフローチャートである。図２に示すフローチャートは、ＣＰＵ１０８が記憶部１０９に記憶されている制御プログラムを実行し、情報の演算および加工並びに各ハードウェアの制御を実行することにより実現され得る。
【００１８】
上述したように、前処理部１０６によって得られた放射線画像は、ＣＰＵバス１０７を介して画像処理部１１２に転送される。転送された放射線画像は、分割部１１３に入力される。分割部１１３は、入力された放射線画像からセグメンテーションマップ（多値画像）を作成する（Ｓ２０１）。具体的には、分割部１１３は、放射線画像の各画素に対し、その画素が属する解剖学的領域を示すラベルを付与する。このような、画像を任意の解剖学的領域に分割することは、セマンティックセグメンテーション（意味的領域分割）と呼ばれる。なお、本実施形態では、当該ラベルは、画素値によって識別可能なラベル番号とする。すなわち、分割部１１３は、複数の解剖学的領域に異なるラベル番号を付与することにより、放射線画像を複数の解剖学的領域に分割する。
【００１９】
図８（ａ）に、Ｓ２０１で作成されるセグメンテーションマップの例を示す。図８（ａ）は、入力された放射線画像が胸部を写す画像の場合のセグメンテーションマップの一例を示す図である。図８（ａ）のように、分割部１１３は、肺野８０１に属する領域における画素に同じ値（例えば画素値０）を与え、脊椎８０２に属する領域における画素に、肺野とは異なる同一値（例えば画素値１）を与える。また、肺野８０１、脊椎８０２の何れにも属さない被写体構造８０３に属する画素には、分割部１１３は、さらに別の同一値（例えば画素値２）を与える。
【００２０】
なお、図８（ａ）で示した分割は一例であり、放射線画像をどのような解剖学的領域の粒度で分割するかは特に限定されない。例えば、分割部１１３は、操作者により操作部１１０を介して指示された、線量指標値を算出するための所望の領域に応じて、分割の粒度を適宜決定しても良い。また、分割部１１３は、被写体構造以外の領域についても同様にラベル番号（画素値）を付与することによりセグメンテーションマップを作成してもよい。例えば、図８（ａ）において、放射線が直接に放射線検出器１０４に到達する領域８０４や放射線がコリメータ（不図示）により遮蔽された領域８０５に関しても其々別のラベル番号を付与したセグメンテーションマップを作成することも可能である。
【００２１】
本実施形態では、分割部１１３は、Ｓ２０１におけるセグメンテーションマップの作成を、公知の方法を用いて行うこととする。公知の方法の一例として、ここではＣＮＮ（Convolutional Neural Network）を用いる。ＣＮＮとは、畳込み層、プーリング層、全結合層などから構成されるニューラルネットワークであり、解決する問題に応じて各層を適当に組み合わせることで実現される。また、ＣＮＮは、機械学習アルゴリズムの一種であり、事前学習を必要とする。具体的には、畳み込み層で用いられるフィルタ係数や、各層の重みやバイアス値などのパラメータ（変数）を、多数の学習データを用いた、いわゆる教師あり学習によって調整（最適化）する必要がある。
【００２２】
ここで、教師あり学習は、ＣＮＮに入力される画像（入力画像）と当該入力画像が与えられたときに期待される出力結果（正解）の組み合わせのサンプルを多数用意し、期待する結果が出力されるようにパラメータを繰り返し調整することを含む。この調整には、一般には、誤差逆伝搬法（バックプロバケーション）が用いられる。すなわち、正解と実際の出力結果の差（損失関数で定義された誤差）が小さくなる方向に、ＣＮＮのパラメータが繰り返し調整される。
【００２３】
なお、本実施形態では、ＣＮＮへの入力画像は、前処理部１０６によって得られた放射線画像とし、期待される出力結果は、正解のセグメンテーションマップ（例えば図８（ａ））となる。このような正解のセグメンテーションマップは、予め操作者により、所望の解剖学的領域の粒度に応じて作成され得る。ＣＮＮのパラメータ（学習されたパラメータ２０２）は、入力画像と期待される出力結果との組み合わせの複数のサンプルを用いて、機械学習により予め生成される。学習されたパラメータ２０２は予め記憶部１０９に記憶され、分割部１１３がＳ２０１でセグメンテーションマップを作成する際に、学習されたパラメータ２０２を呼び出し、ＣＮＮによるセマンティックセグメンテーションを行う。すなわち、分割部１１３は、予め機械学習によって生成されたパラメータを用いて、放射線画像を複数の解剖学的領域に分割する。また、操作者は、所定の画像データと当該所定の画像データに対応する正解のセグメンテ―ションマップ（分割割り当てデータ）を、教師データとして決定することができる。当該教師データは、設定部１１６により管理される。
【００２４】
なお、全部位を合わせたデータを用いて唯一の学習されたパラメータ２０２が生成されてもよいが、部位毎（例えば、頭部、胸部、腹部、四肢など）に、学習されたパラメータ２０２が生成されても良い。すなわち、部位毎の、入力画像と期待される出力結果の組み合わせの複数のサンプルを用いて、別々に学習を行って、各部位に対する、複数のセットの学習されたパラメータ２０２が生成されても良い。複数のセットの学習されたパラメータ２０２が生成された場合は、それぞれのセットの学習されたパラメータを部位情報と関連付けて記憶部１０９に予め記憶し、分割部１１３は、入力画像の部位に応じて記憶部１０９から撮影部位に対応する学習されたパラメータ２０２を呼び出し、ＣＮＮによるセマンティックセグメンテーションを行えば良い。
【００２５】
なお、ＣＮＮのネットワーク構造については特に限定するものではなく、一般的に知られたものを用いれば良い。具体的には、機械学習として、ＦＣＮ（Fully Convolutional Networks）、ＳｅｇＮｅｔ、Ｕ−ｎｅｔ等を用いることができる。また、本実施形態では、ＣＮＮへの入力画像を前処理部１０６によって得られた放射線画像としたが、当該放射線画像を縮小した放射線画像をＣＮＮへの入力画像として用いても良い。ＣＮＮによるセマンティックセグメンテーションは、計算量が多く計算時間が掛かるため、縮小した画像データを用いることにより、計算時間の削減につながり得る。
【００２６】
次に、抽出部１１４は、線量指標値を算出するための領域を抽出する（Ｓ２０３）。具体的には、抽出部１１４は、操作者により操作部１１０を介して指定された領域を、線量指標値を算出するための領域として抽出する。本実施形態では、当該抽出処理のために、領域情報として、撮影部位とラベル番号の対応情報２０４を用いる。すなわち、抽出部１１４は、予め設定された、複数の部位と当該複数の部位のそれぞれに対応するラベル番号との対応情報２０４を用いて、操作者の指示に従って、複数の解剖学的領域のうちの少なくとも一つの領域を抽出する。図９は、対応情報２０４の一例としての、複数の撮影部位と当該複数の部位の夫々に対応するラベル番号との対応表の一例を示す。図９に示す対応表９００は、操作者等により予め作成され、記憶部１０９に格納されているものとする。抽出部１１４は、対応表９００を参照し、操作者により操作部１１０を介して指定（指示）された領域の部位に対応するラベル番号を取得する。なお、抽出部１１４は、操作者により操作部１１０を介して指定されたラベル番号を直接取得しても良い。続いて、抽出部１１４は、取得したラベル番号に該当する画素の値を１としたマスク画像（Ｍａｓｋ）を、（数１）にて生成する。
（数１）
2020025730.tif 000003
ただし、Ｍａｐはセグメンテーションマップを表し、（ｘ、ｙ）は画像の座標を表す。また、Ｌは取得されたラベル番号を表す。
【００２７】
なお、操作者が指定する領域は１つに限るのもではなく、複数であっても良い。操作者が複数の領域を指定した場合は、抽出部１１４は、複数の領域に対応する複数のラベル番号の何れかに該当する画素の値を１、それ以外を０としてマスク画像を生成すれば良い。なお、操作者が領域を設定する方法の詳細に関しては、図３のフローチャートを用いて後述する。
【００２８】
次に、算出部１１５は、Ｓ２０３で抽出された領域（すなわち、マスク画像のうち、画素値が１の領域）における代表値として、抽出された領域の中心傾向を示す値Ｖを算出する（Ｓ２０５）。本実施形態では、値Ｖを、（数２）のように算出する。
（数２）
2020025730.tif 000004
ただし、Ｏｒｇは入力画像（本実施形態では前処理部１０６によって得られた放射線画像）、Ｍａｓｋはマスク画像を表し、（ｘ、ｙ）は画像の座標を表し、Ｏｒｇ（ｘ、ｙ）は、入力画像における座標（ｘ、ｙ）における画素値を表す。
【００２９】
次に、算出部１１５は、求めた値Ｖを、線量指標値ＥＩに変換する（Ｓ２０６）。具体的には、算出部１１５は、国際規格ＩＥＣ６２４９４−１の定義に従って、（数３）のように値Ｖを線量指標値ＥＩに変換する。
（数３）
2020025730.tif 000005
ただし、関数ｇは値Ｖを空気カーマに変換するための関数であり、規定された条件下における値Ｖと空気カーマの関係に従って予め定められる。なお、関数ｇは、放射線検出器１０４の特性により異なる。そのため、操作者は、使用し得る放射線検出器１０４に応じて複数の関数ｇを予め記憶部１０９に記憶しておき、算出部１１５は、実際に使用する放射線検出器１０４に対応する関数ｇを用いて変換を行い得る。
【００３０】
また、算出部１１５にて、線量目標値ＥＩＴと線量指標値ＥＩとのずれ量（偏差）ＤＩを（数４）にて算出する（Ｓ２０７）。
（数４）
2020025730.tif 000006
ただし、線量目標値ＥＩＴは、操作者が定める線量管理基準に従って部位毎に予め設定されるものとする。偏差ＤＩは、線量目標値ＥＩＴと線量指標値ＥＩとの偏差を表す数値であるから、線量目標値ＥＩＴと線量指標値ＥＩが同じであれば０となる。また、線量指標値ＥＩが大きくなるほど、すなわち、撮影された画像の線量が線量目標値ＥＩＴよりも多くなるほど、偏差ＤＩは大きくなる。例えば、撮影された画像の線量が線量目標値ＥＩＴの２倍であれば偏差ＤＩは約３となる。また、撮影された画像の線量が線量目標値ＥＩＴよりも少ないほど偏差ＤＩは小さくなり、例えば、撮影された画像の線量が線量目標値の半分であれば偏差ＤＩは約−３となる。これにより、操作者は、撮影された画像の線量が、基準である線量目標値ＥＩＴに対して不足しているのか、逆に線量過多の状態なのかを、瞬時に把握することができる。
【００３１】
次に、ＣＰＵ１０８は、求めた線量指標値ＥＩおよび偏差ＤＩを、表示部１１１に表示する（Ｓ２０８）。ここで、表示の方法は特に限定するものではなく、例えば、ＣＰＵ１０８は、表示部１１１において、前処理部１０６により得られた放射線画像と合わせて表示することができる、また、ＣＰＵ１０８は、表示部１１１による表示エリアの一部に、求めた線量指標値ＥＩおよび偏差ＤＩをオーバーレイ表示するように制御しても良い。
【００３２】
なお、本実施形態では、操作者が指定した領域から、唯一の線量指標値ＥＩおよび偏差ＤＩを算出したが、例えば、複数の線量指標値ＥＩおよび偏差ＤＩを求める構成としても良い。具体的には、操作者が領域を複数指定することにより、複数の領域が抽出された場合は、放射線撮影装置１００は、複数の領域それぞれに対して値Ｄを算出し、当該値Ｄからそれぞれの線量指標値ＥＩおよび偏差ＤＩを算出し、表示部１１１に表示しても良い。
【００３３】
以上、撮影により得られた放射線画像から線量指標値を算出する動作について説明した。次に、操作者が線量指標値を算出するための領域を変更する場合の画像処理部１１２の動作に関して、図３のフローチャートを用いて説明する。図３は、線量指標値を算出するための領域を変更する処理手順を示すフローチャートである。図３に示すフローチャートは、ＣＰＵ１０８が記憶部１０９に記憶されている制御プログラムを実行し、情報の演算および加工並びに各ハードウェアの制御を実行することにより実現され得る。なお、図３のフローチャートにおいて、図２に示したステップと同じステップには同じ符号を付し、説明を省略する。
【００３４】
まず、操作者は、操作部１１０を介して、変更しようとする領域の部位（当該領域が属する部位）を選択する（Ｓ３０１）。次に、設定部１１６は、選択された部位に対応する教師データがあるか否かを判定する（Ｓ３０２）。ここで、教師データとは、上述した通り、操作者により決定された、正解のセグメンテーションマップ（分割割り当てデータ）を指す。選択された部位に対応する教師データがある場合（Ｓ３０２でＹｅｓ）、設定部１１６は、教師データを記憶部１０９から取得する（Ｓ３０３）。選択された部位に対応する教師データがない場合（Ｓ３０２でＮｏ）、設定部１１６は、記憶部１０９から過去の撮影により得られた選択された部位の画像データ（放射線画像）を取得し（Ｓ３０４）、分割部１１３が取得した画像データからセグメンテーションマップを作成する（Ｓ２０１）。セグメンテーションマップの作成方法は、図２のＳ２０１において説明した通りである。
【００３５】
次に、ＣＰＵ１０８は、図８（ａ）のようなセグメンテーションマップを表示部１１１に表示する（Ｓ３０５）。操作者は、変更する領域である、線量指標値を計算するための任意の領域を、図示しないマウス等を使って指定する（Ｓ３０６）。次に、設定部１１６は、指定された領域に対応するラベル番号を取得し、対応情報２０４（例えば図９に示した撮影部位とラベル番号の対応表９００）を更新し、記憶部１０９に保存（設定）する（Ｓ３０７）。
【００３６】
以上、実施形態１では、放射線画像において分割された領域の中から、線量指標値を算出するための領域を自由に変更することが可能となり、操作者が意図した基準で適切な線量管理を行うことが可能となる。
【００３７】
［実施形態２］
（放射線撮影装置の構成）
実施形態２による放射線撮影装置４００の構成例を図４に示す。放射線撮影装置４００は、図１に示した放射線撮影装置１００に対し、機械学習部４０１を備える構成である。機械学習部４０１は、線量指標値を算出するための領域の分割構成を変更するための学習を行う機能を有する。具体的には、機械学習部４０１は、教師データ（所定の画像データと当該所定の画像データに対応する正解のセグメンテ―ションマップ（分割割り当てデータ））に基づいて機械学習（ＣＮＮの再学習）を行う。
【００３８】
（放射線撮影装置の動作）
以下、実施形態１と異なる動作である、線量指標値を算出するための領域の分割構成を変更処理について、図５を用いて説明する。図５は、線量指標値を算出するための領域の分割構成を変更する処理手順を示すフローチャートである。図５に示すフローチャートは、ＣＰＵ１０８が記憶部１０９に記憶されている制御プログラムを実行し、情報の演算および加工並びに各ハードウェアの制御を実行することにより実現され得る。
【００３９】
まず、機械学習部４０１は、教師データ５０２に基づいてＣＮＮの再学習を行う（Ｓ５０１）。この再学習は、予め準備した教師データ５０２を用いた学習である。なお、学習の具体的な方法については、実施形態１で説明したと同様に、誤差逆伝搬法（バックプロバケーション）を用い、正解と実際の出力結果の差（損失関数で定義された誤差）が小さくなる方向にＣＮＮのパラメータを繰り返し調整することで行われる。
【００４０】
ここで、設定部１１６は、機械学習部４０１による再学習のための教師データを以下のように設定することができる。すなわち、設定部１１６は、教師データである、正解のセグメンテーションマップを変更し、機械学習部４０１に設定することができる。図８（ｂ）に変更した正解のセグメンテーションマップの一例を示す。例えば、実施形態１では、分割部１１３は、図８（ａ）に示すように、肺野８０１を一つの領域ととみなして同一のラベルを付与した。一方、本実施形態では、設定部１１６は、肺野を右肺野と左肺野で異なる領域として分割する場合を考慮して、図８（ｂ）のように右肺野８０１ａと左肺野８０１ｂに対し異なるラベルを付与した正解のセグメンテーションマップを、教師データ５０２として準備する。また、実施形態１では、分割部１１３は、図８（ａ）に示すように、脊椎８０２を１つの領域とみなして同一のラベルを付与した。一方、本実施形態では、設定部１１６は、脊椎を胸椎と腰椎で分割する場合を考慮して、図８（ｂ）のように胸椎８０２ａ、腰椎８０２ｂに対し異なるラベルを付与した正解のセグメンテーションマップを準備する。
【００４１】
次に、機械学習部４０１は、再学習により得られたパラメータをＣＮＮの新たなパラメータとして、記憶部１０９において更新（保存）する（Ｓ５０３）。続いて、画像処理部１１２は、線量指標値を計算する領域を再設定する（Ｓ５０４）。ここで、再設定の方法は図３のフローチャートの動作と同様であるため説明を省略する。なお、図３では、上記処理で新たな教師データを用いることにより、新たなセグメンテーションマップが作成され、Ｓ３０７の処理では、対応情報２０４（例えば図９に示した撮影部位とラベル番号の対応表９００）が更新される。すなわち、設定部１１６は、更新されたパラメータを用いて分割された複数の解剖学的領域に応じて、対応情報２０４を更新して設定する。これにより、例えば次回以降の撮影（図２）では、学習されたパラメータ２０２（図２）をＳ５０３で更新されたパラメータに替えることで、分割部１１３は、更新されたパラメータを用いて、放射線画像を複数の解剖学的領域に分割することができる。さらに、Ｓ２０３で用いる対応情報２０４を、本実施形態により更新されたものに替えて用いることで、新たに定義された領域での線量指標値の算出が可能となる。
【００４２】
以上、実施形態２では、線量指標値を算出するための領域の分割構成を変更することが可能となり、操作者が線量指標値を算出する領域を自由に変更できるという効果がある。
【００４３】
［実施形態３］
（放射線撮影装置の構成）
実施形態３による放射線撮影装置６００の構成例を図６に示す。放射線撮影装置６００は、図４に示した放射線撮影装置４００に対し、目標値更新部６０１を備える構成である。目標値更新部６０１は、線量目標値ＥＩＴを自動で設定する機能を有する。
【００４４】
線量目標値ＥＩＴは，線量指標値ＥＩの基準となる値であり、線量指標値ＥＩを算出するための領域が変われば線量目標値ＥＩＴも変更する必要がある。この変更は線量管理基準に合わせて操作者が手動で設定するものであるが、線量指標値ＥＩを算出するための領域が変更されるたびに一から設定するのは非常に煩雑である。そこで、目標値更新部６０１は、線量指標値ＥＩを算出するための領域の変更前後の差に基づいて、変更前の線量目標値ＥＩＴと略等価の値を用いて自動で更新する機能を有する。
【００４５】
（放射線撮影装置の動作）
以下、線量目標値ＥＩＴを自動で更新する方法について、図７を用いて説明する。図７は、線量目標値ＥＩＴを自動で更新する処理手順を示すフローチャートである。図７に示すフローチャートは、ＣＰＵ１０８が記憶部１０９に記憶されている制御プログラムを実行し、情報の演算および加工並びに各ハードウェアの制御を実行することにより実現され得る。なお、このフローチャートは、線量目標値を算出するための領域が変更されたタイミングで動作するものである。具体的には、線量指標値を算出するための領域を変更した場合（図３のフローチャートに従った動作）や、線量指標値を算出するための領域の分割構成を変更した場合（図５のフローチャートに従った動作）の後に実行される。
【００４６】
まず、目標値更新部６０１は、線量目標値を算出するための領域に現在設定されているＥＩＴを取得する（Ｓ７０１）。次に、算出部１１５は、記憶部１０９から、過去の撮影よる複数の放射線画像をロードし各々の放射線画像に対する線量指標値ＥＩを算出する（Ｓ７０２）。ここで、線量指標値ＥＩの算出方法は図２のフローチャートで説明した通りである。なお、線量指標値ＥＩを算出する際に用いる学習されたパラメータ２０２および、対応情報２０４は、設定変更（線量指標値を算出するための領域の変更（図３）や、線量指標値を算出するための領域の分割構成の変更（図５））の前のものが利用されるものとする。次に、算出部１１５は、設定変更の後の線量指標値ＥＩの算出をＳ７０２と同様に実行する（Ｓ７０３）。Ｓ７０２との違いは、設定変更の後の学習されたパラメータ２０２および対応情報２０４を利用する点である。
【００４７】
次に、設定変更前後でのＥＩの誤差Ｅｒｒを（数５）にて算出する（Ｓ７０４）。
（数５）
2020025730.tif 000007
ただし、ｋは画像番号を表し、ＮはＥＩを算出した画像の総数を表す。また、ＥＩ１（ｋ）、ＥＩ２（ｋ）は画像番号ｋの画像から算出したＥＩであり、ＥＩ１が設定変更前、ＥＩ２が設定変更後のＥＩを表している。
【００４８】
次に求めた誤差Ｅｒｒにて線量目標値ＥＩＴを更新する（Ｓ７０５）。
（数６）
2020025730.tif 000008
ただし、ＥＩＴ１が更新前の線量目標値であり、ＥＩＴ２が更新後の線量目標値を表す。
【００４９】
このように、目標値更新部６０１は、学習されたパラメータ２０２が更新された場合、または、対応情報２０４が更新された場合に、当該更新の前後に算出部１１５により算出されたそれぞれの線量指標値を用いて、線量指標値の目標値としての線量目標値を更新する。以上、実施形態３では設定変更時にＥＩＴの値を自動で更新することで、操作者の手間を低減できる効果がある。
【００５０】
以上、いくつかの実施形態について説明したが、本発明はこれらの実施形態に限定されないことはいうまでもなく、その要旨の範囲内で種々の変形及び変更が可能である。
【００５１】
［その他の実施形態］
本発明は、上述の実施形態の１以上の機能を実現するプログラムを、ネットワーク又は記憶媒体を介してシステム又は装置に供給し、そのシステム又は装置のコンピュータにおける１つ以上のプロセッサーがプログラムを読出し実行する処理でも実現可能である。また、１以上の機能を実現する回路（例えば、ＡＳＩＣ）によっても実現可能である。
【符号の説明】
【００５２】
１００  放射線撮影装置、１０１  放射線発生部、１０２  放射線ビーム、１０３  被検体、１０４  放射線検出器、１０５  データ収集部、１０６  前処理部、１０７  ＣＰＵバス、１０８  ＣＰＵ、１０９  記憶部、１１０  操作部、１１１  表示部、１１２  画像処理部、１１３  分割部、１１４  抽出部、１１５  線量指標値算出部、１１６  設定部、４０１  機械学習部、６０１  線量目標値更新部
【図１】
2020025730.tif 000009
【図２】
2020025730.tif 000010
【図３】
2020025730.tif 000011
【図４】
2020025730.tif 000012
【図５】
2020025730.tif 000013
【図６】
2020025730.tif 000014
【図７】
2020025730.tif 000015
【図８】
2020025730.tif 000016
【図９】
2020025730.tif 000017