(57)【要約】【課題】Ｘ線撮影で得られる画像から散乱線成分を低減した画像を直接的に導出することができるとともに、計算時間が短い散散乱線補正装置を提供する。【解決手段】本発明に係る散乱線補正装置は、Ｘ線撮影で得られる測定画像に対して散乱線補正を行う。また、本発明に係る散乱線補正装置は、正規化された全Ｘ線による輝度値である正規化理想値を前記測定画像の各ピクセルについて算出する正規化理想値算出部と、前記測定画像のピクセル毎に、前記正規化理想値と、正規化された前記測定画像の輝度値とから、前記測定画像の散乱線成分を算出する散乱線成分算出部と、前記測定画像のピクセル毎に、前記測定画像のピクセルの輝度値から、前記散乱線成分算出処理で算出した散乱線成分を除去する補正処理とを備える。【選択図】図９
2015065991.tif 000002
【特許請求の範囲】
【請求項１】
Ｘ線撮影で得られる測定画像に対して散乱線補正を行う散乱線補正装置であって、
正規化された全Ｘ線による輝度値である正規化理想値を前記測定画像の各ピクセルについて算出する正規化理想値算出部と、
前記測定画像のピクセル毎に、前記正規化理想値と、正規化された前記測定画像の輝度値とから、前記測定画像の散乱線成分を算出する散乱線成分算出部と、
前記測定画像のピクセル毎に、前記測定画像のピクセルの輝度値から、前記散乱線成分算出処理部で算出した散乱線成分を除去する補正処理とを備えることを特徴とする散乱線補正装置。
【請求項２】
前記散乱線成分算出部が、前記測定画像のピクセル毎に、前記正規化理想値と、正規化された前記測定画像の輝度値と、前記Ｘ線撮影の投影角度とから、前記測定画像の散乱線成分を算出する請求項１に記載の散乱線補正装置。
【請求項３】
正規化された全Ｘ線による輝度値と、全Ｘ線量に対するその全Ｘ線量に含まれる散乱線量の比との関係を記憶する記憶部を備え、
前記散乱線成分算出部が、前記記憶部によって記憶されている前記関係を用いて、前記測定画像の散乱線成分を算出する請求項１または請求項２に記載の散乱線補正装置。
【請求項４】
前記記憶部が、前記関係を示す関数で用いられる定数の値を、前記正規化理想値に関連付けられたデータテーブルの形式で記憶する請求項３に記載の散乱線補正装置。
【請求項５】
前記正規化理想値算出部が、前記測定画像をそれぞれの領域が複数のピクセルの集合体である複数の領域に分割し、領域を代表する全Ｘ線による輝度値を領域毎に求め、領域を代表する全Ｘ線による輝度値を参照して前記正規化理想値を算出する請求項１〜４のいずれか１項に記載の散乱線補正装置。
【請求項６】
前記正規化理想値算出部が、シミュレーションによって求めた前記正規化理想値を、領域を代表する全Ｘ線による輝度値を参照して修正し、修正後の前記正規化理想値を算出結果とする請求項５に記載の散乱線補正装置。
【請求項７】
前記正規化理想値が、被写体の着目する部位を透過するＸ線経路上の骨あるいは金属を皮膚に置き換えた場合に得られる請求項１〜６のいずれか１項に記載の散乱線補正装置。
【請求項８】
Ｘ線撮影で得られる測定画像に対して散乱線補正を行う散乱線補正方法であって、
正規化された全Ｘ線による輝度値である正規化理想値を前記測定画像の各ピクセルについて算出する正規化理想値算出ステップと、
前記測定画像のピクセル毎に、前記正規化理想値と、正規化された前記測定画像の輝度値とから、前記測定画像の散乱線成分を算出する散乱線成分算出ステップと、
前記測定画像のピクセル毎に、前記測定画像のピクセルの輝度値から、前記散乱線成分算出処理ステップで算出した散乱線成分を除去する補正処理ステップとを備えることを特徴とする散乱線補正方法。
【請求項９】
被写体に対してＸ線を照射するＸ線照射部と、
前記被写体を透過したＸ線を検出するＸ線検出部と、
前記Ｘ線検出部の検出結果を用いて測定画像を生成する測定画像生成部と、
前記測定画像を対数変換して投影画像を生成する投影画像生成部と、
前記測定画像に対して散乱線補正を行う請求項１〜７のいずれか１項に記載の散乱線補正装置とを備えることを特徴とするＸ線撮影装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、Ｘ線撮影で得られる画像に含まれる散乱線成分を散乱線補正により低減する散乱線補正装置及び散乱線補正方法並びに散乱線補正装置を備えるＸ線撮影装置に関する。
【背景技術】
【０００２】
Ｘ線撮影で得られる画像は、被写体を透過した直接線に起因する直接線成分と、直接線成分以外の成分とを含む。直接線成分以外の成分としては、散乱線に起因する散乱線成分と、Ｘ線検出素子の光拡散およびＸ線撮影装置の系から発生する電気的ノイズ等に起因する成分とが存在するが、散乱線成分の方が支配的である。散乱線としては、被写体からの散乱線やフラットパネルディテクターの前面に設置されるカーボンからの散乱線などがある。
【０００３】
Ｘ線撮影で得られる画像に含まれる散乱線成分は、Ｘ線撮影で得られる画像の品質を大きく低下させており、例えば金属アーチファクトの除去を困難にする原因の一つとなっている。
【０００４】
そのため、Ｘ線撮影で得られる画像に含まれる散乱線成分を低減するための散乱線成分低減方法が種々提案されている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特許第４２１８９０８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
特許文献１で提案されている散乱線成分低減方法は、散乱線を含むデータから再構成ボリュームデータを得て、その再構成ボリュームデータ上で投影画像に含まれる散乱線の分布を推定し、投影画像から散乱線成分を減算した後、もう一度再構成するという処理を行っているため計算時間がかかる。
【０００７】
本発明は、上記の状況に鑑み、Ｘ線撮影で得られる画像から散乱線成分を低減した画像を直接的に導出することができるとともに、計算時間が短い散乱線補正装置、散乱線補正方法、及びＸ線撮影装置を提供することを目的とするものである。
【課題を解決するための手段】
【０００８】
上記目的を達成するために本発明に係る散乱線補正装置においては、Ｘ線撮影で得られる測定画像に対して散乱線補正を行う散乱線補正装置とし、正規化された全Ｘ線による輝度値である正規化理想値を前記測定画像の各ピクセルについて算出する正規化理想値算出部と、前記測定画像のピクセル毎に、前記正規化理想値と、正規化された前記測定画像の輝度値とから、前記測定画像の散乱線成分を算出する散乱線成分算出部と、前記測定画像のピクセル毎に、前記測定画像のピクセルの輝度値から、前記散乱線成分算出処理部で算出した散乱線成分を除去する補正処理とを備える構成（第１の構成）とする。これにより、Ｘ線撮影で得られる測定画像から散乱線成分を低減した画像を直接的に導出することができる。また、特許文献１のような二度の再構成が不要であるため、計算時間を短くすることができる。歯科用ＣＴ撮影では被写体が人間の頭部（特に顎部）に限定されているため、計算時間を短くすることが容易になる。
【０００９】
上記第１の構成の散乱線補正装置において、前記散乱線成分算出部が、前記測定画像のピクセル毎に、前記正規化理想値と、正規化された前記測定画像の輝度値と、前記Ｘ線撮影の投影角度とから、前記測定画像の散乱線成分を算出する構成（第２の構成）としてもよい。これにより、Ｘ線撮影の投影角度に応じて適切な散乱線補正を行うことができる。
【００１０】
上記第１又は第２の構成の散乱線補正装置において、正規化された全Ｘ線による輝度値と、全Ｘ線量に対するその全Ｘ線量に含まれる散乱線量の比との関係を記憶する記憶部を備え、前記散乱線成分算出部が、前記記憶部によって記憶されている前記関係を用いて、前記測定画像の散乱線成分を算出する構成（第３の構成）としてもよい。
【００１１】
上記第３の構成の散乱線補正装置において、前記記憶部が、前記関係を示す関数で用いられる定数の値を、前記正規化理想値に関連付けられたデータテーブルの形式で記憶する構成（第４の構成）としてもよい。これにより、記憶部の使用容量を節約することができる。
【００１２】
上記第１〜第４のいずれかの構成の散乱線補正装置において、前記正規化理想値算出部が、前記測定画像をそれぞれの領域が複数のピクセルの集合体である複数の領域に分割し、領域を代表する全Ｘ線による輝度値（例えば領域の中心位置の全Ｘ線による輝度値）を領域毎に求め、領域を代表する全Ｘ線による輝度値を参照して前記正規化理想値を算出する構成（第５の構成）としてもよい。これにより、各ピクセルの輝度値の誤差による影響を抑えることができる。
【００１３】
上記第５の構成の散乱線補正装置において、前記正規化理想値算出部が、シミュレーションによって求めた前記正規化理想値を、領域を代表する全Ｘ線による輝度値を参照して修正し、修正後の前記正規化理想値を算出結果とする構成（第６の構成）としてもよい。これにより、実際の被写体により適応した散乱線補正を行うことができる。
【００１４】
上記第１〜第６のいずれかの構成の散乱線補正装置において、前記正規化理想値が、被写体の着目する部位を透過するＸ線経路上の骨あるいは金属を皮膚に置き換えた場合に得られる構成（第７の構成）としてもよい。これにより、人体の構造を考慮することができる。
【００１５】
また、上記目的を達成するために本発明に係る散乱線補正方法においては、Ｘ線撮影で得られる測定画像に対して散乱線補正を行う散乱線補正方法とし、正規化された全Ｘ線による輝度値である正規化理想値を前記測定画像の各ピクセルについて算出する正規化理想値算出ステップと、前記測定画像のピクセル毎に、前記正規化理想値と、正規化された前記測定画像の輝度値とから、前記測定画像の散乱線成分を算出する散乱線成分算出ステップと、前記測定画像のピクセル毎に、前記測定画像のピクセルの輝度値から、前記散乱線成分算出処理ステップで算出した散乱線成分を除去する補正処理ステップとを備えるようにしている。
【００１６】
また、上記目的を達成するために本発明に係るＸ線撮影装置においては、被写体に対してＸ線を照射するＸ線照射部と、前記被写体を透過したＸ線を検出するＸ線検出部と、前記Ｘ線検出部の検出結果を用いて測定画像を生成する測定画像生成部と、前記測定画像を対数変換して投影画像を生成する投影画像生成部と、前記測定画像に対して散乱線補正を行う上記第１〜第７のいずれかの構成の散乱線補正装置とを備える構成とする。
【発明の効果】
【００１７】
本発明によると、Ｘ線撮影で得られる測定画像から散乱線成分を低減した画像を直接的に導出することができる。また、計算時間を短くすることができる。
【図面の簡単な説明】
【００１８】
【図１】本発明の一実施形態に係るＸ線撮影装置の本体部の外観を示す図である。
【図２】パノラマ撮影モードの標準軌道を示す図である。
【図３】局所ＣＴ撮影モードの軌道を示す図である。
【図４】局所ＣＴ撮影モードにおいて、撮影対象部位の中心を前歯の位置に設定した場合の撮影対象部位の中心及び画像再構成範囲を示す図である。
【図５】局所ＣＴ撮影モードにおいて、撮影対象部位の中心を左顎の位置に設定した場合の撮影対象部位の中心及び画像再構成範囲を示す図である。
【図６】局所ＣＴ撮影モードにおいて、撮影対象部位の中心を右第２小臼歯の位置に設定した場合の撮影対象部位の中心及び画像再構成範囲を示す図である。
【図７】全歯ＣＴ撮影モードの軌道を示す図である。
【図８】全顎ＣＴ撮影モードの軌道を示す図である。
【図９】画像処理装置の構成を示す図である。
【図１０】モンテカルロシミュレーションの計算のジオメトリーを示す図である。
【図１１】顎部、頭部、頸部を想定した円柱形の水ファントムの組み合せを示す図である。
【図１２】モンテカルロシミュレーションの計算結果を示すグラフである。
【図１３】モンテカルロシミュレーションの計算結果を示すグラフである。
【図１４】モンテカルロシミュレーションの計算結果を示すグラフである。
【図１５】モンテカルロシミュレーションの計算結果を示すグラフである。
【発明を実施するための形態】
【００１９】
本発明の実施形態について図面を参照して以下に説明する。
【００２０】
まず始めに、本発明の一実施形態に係るＸ線撮影装置の本体部１（以下、「Ｘ線撮影装置の本体部１」と称す）の構成について図１を参照して説明する。図１はＸ線撮影装置１の本体部の外観を示す図であり、図１（ａ）は上面図、図１（ｂ）は正面図、図１（ｃ）は側面図である。
【００２１】
Ｘ線撮影装置の本体部１は、歯科用あるいは耳鼻科用等のＸ線撮影装置の本体部であって、床面に載置されるベース２と、ベース２から鉛直方向に立設された下部ポール３と、鉛直方向にスライド可能に下部ポール３に接続される上部ポール４と、上部ポール４の上端部に固定されている固定アーム５と、回転可能に固定アーム５に接続される旋回アーム６と、上部ポール４の中央部に固定されており被写体（例えば歯など）を含む人体の頭部を保持する頭部保持部７とを備えている。実施形態では、固定アーム５が上部ポール４に固定されているが、例えば、Ｘ線撮影装置の本体部１を設置する部屋の壁や天井に固定アーム５が直接あるいは部屋の壁や天井との距離を調整することができる調整機構を介して取り付けられる態様であってもよい。
【００２２】
旋回アーム６は、被写体に対してＸ線を照射するＸ線照射部８と、被写体を透過したＸ線を検出するＸ線検出部９とを対向して配置している。本実施形態では、Ｘ線検出部９として、照射されたＸ線に応じて電気信号を生成する変換素子が二次元状に配置されているフラットパネルディテクターを用いる。そして、フラットパネルディテクターの前面にはカーボンが設置されている。
【００２３】
Ｘ線撮影装置の本体部１の撮影モードは特に限定されないが、例えば、パノラマ撮影モードやＣＴ撮影モードを挙げることができる。パノラマ撮影モードでは、Ｘ線照射部８及びＸ線検出部９が歯列弓の形状に沿った所定の軌跡を描くように、旋回アーム６の旋回軸を旋回軸に垂直な方向（Ｘ方向、Ｙ方向）に移動させ、旋回アーム６を旋回軸回りに旋回させながら断層撮影を行う。ＣＴ撮影モードでは、頭部の対象撮影領域（画像再構成範囲）を中心にして旋回アーム６を回転させながら、対象撮影領域（画像再構成範囲）の断層撮影を行う。
【００２４】
ここで、パノラマ撮影モードについて図２を参照してより詳細に説明する。図２はパノラマ撮影モードの標準（成人用）軌道を示している。パノラマ撮影モードの標準（成人用）軌道では、Ｘ線照射部８及びＸ線検出部９が仮想歯列弓２０１の形状に沿った所定の軌跡を描いてＸ線ビームの軌跡が包絡線状の軌跡２０２になるように、Ｘ線照射部８及びＸ線検出部９が配置されている旋回アーム６を、撮影開始位置Ｐ１から図２に示す軌道に沿って撮影終了位置Ｐ２まで移動させる。撮影開始位置Ｐ１と撮影終了位置Ｐ２との間における旋回アーム６の旋回角度は約２２０度である。なお、撮影終了位置Ｐ２を除く図２に示された旋回アーム６の位置は被写体の撮影領域の左半分における各撮影位置である。Ｘ線照射部８のＸ線焦点８Ａから射出されるＸ線は、Ｘ線照射部８に設けられているＸ線絞り８Ｂによって絞られ、Ｘ線検出部９上でのＸ線ビーム幅Ｗが調整される。
【００２５】
パノラマ撮影モードは、上述した標準（成人用）軌道の他に、小児用軌道、直行軌道、顎関節撮影軌道、上顎洞撮影軌道などを有していることが好ましい。小児用軌道は、仮想歯列弓２０１の形状が小さくなる点が標準軌道と異なっている。直行軌道は、各撮影位置でのＸ線ビームが患者歯列弓２０３の歯と歯の間を通過するようにしている点が標準軌道と異なっている。顎関節撮影軌道（側面）は、Ｘ線照射部８及びＸ線検出部９が仮想歯列弓２０１の両端部分（顎関節撮影可能部分）の形状に沿った所定の軌跡を描くように旋回アーム６を移動させる点が標準軌道と異なっている。顎関節撮影軌道（正面）は、Ｘ線照射部８及びＸ線検出部９が仮想線２０４の形状に沿った所定の軌跡を描くように旋回アーム６を移動させる点が標準軌道と異なっている。上顎洞撮影軌道は、Ｘ線照射部８及びＸ線検出部９が仮想線２０５の形状に沿った所定の軌跡を描くように旋回アーム６を移動させる点が標準軌道と異なっている。
【００２６】
続いて、ＣＴ撮影モードについて図３〜図８を参照してより詳細に説明する。なお、図３〜図８において図２と同一の部分には同一の符号を付す。
【００２７】
局所ＣＴ撮影モードは、歯顎領域内の上下歯牙領域全体よりも狭い特定の領域を撮影対象とするＣＴ撮影モードである。局所ＣＴ撮影モードの画像再構成範囲は例えば直径５１ｍｍ高さ５５ｍｍの円柱形状の空間領域である。図３は局所ＣＴ撮影モードの軌道を示している。局所ＣＴ撮影モードでは、図３に示すように、Ｘ線検出部９の中心がＸ線照射部８と旋回アーム６の旋回軸中心２０６とを結ぶラインの延長線上にくるように旋回アーム６を旋回させながら複数の撮影位置で撮影が行われる。また、局所ＣＴ撮影モードでは、通常、図３に示すように、旋回アーム６の旋回軸中心２０６は定位置になっている。なお、図３には撮影位置として４箇所が図示されているが、これはあくまで例示であり撮影位置は図示された箇所に限定されるものではない。
【００２８】
局所ＣＴ撮影モードは、後述する全歯ＣＴ撮影モードや全顎ＣＴ撮影モードに比べてＸ線検出部９上でのＸ線ビーム幅Ｗが狭いため、Ｘ線検出部９のサイズが小さくても実施可能である。
【００２９】
なお、局所ＣＴ撮影モードでは、撮影対象部位（関心領域）の中心を何処に設定するかに応じて旋回アーム６の旋回軸中心２０６の位置を変えるようにしており、通常、図３に示すように、撮影対象部位（関心領域）の中心と旋回アーム６の旋回軸中心２０６の位置とが一致するように位置調整がなされる。局所ＣＴ撮影モードにおける撮影対象部位（関心領域）の中心は任意に設定することができる。図３に示した位置設定の他にも、例えば、図４に示すように撮影対象部位（関心領域）の中心２０８を仮想歯列弓２０１上の前歯の位置に設定することもでき、図５に示すように撮影対象部位（関心領域）の中心２０８を仮想歯列弓２０１上の左顎の位置に設定することもでき、図６に示すように撮影対象部位（関心領域）の中心２０８を仮想歯列弓２０１上の右第２小臼歯の位置に設定することもでき、その他種々の位置設定が可能である。
【００３０】
全歯ＣＴ撮影モードは、上下歯牙領域全体を撮影対象とするＣＴ撮影モードである。全歯ＣＴ撮影モードの画像再構成範囲は例えば直径９７ｍｍ高さ１００ｍｍの円柱形状の空間領域である。図７は全歯ＣＴ撮影モードの軌道を示している。全歯ＣＴ撮影モードでは、図７に示すように、Ｘ線検出部９の中心がＸ線照射部８と旋回アーム６の旋回軸中心２０６とを結ぶラインの延長線上にくるように旋回アーム６を旋回させながら複数の撮影位置で撮影が行われる。また、全歯ＣＴ撮影モードでは、通常、図７に示すように、旋回アーム６の旋回軸中心２０６は定位置になっている。なお、図７には撮影位置として４箇所が図示されているが、これはあくまで例示であり撮影位置は図示された箇所に限定されるものではない。
【００３１】
全歯ＣＴ撮影モードは、上述した局所ＣＴ撮影モードに比べて撮影対象が広範囲になりＸ線検出部９上でのＸ線ビーム幅Ｗが広くなるため、その広いＸ線ビーム幅Ｗに見合ったＸ線検出部９のサイズを必要とする。
【００３２】
全顎ＣＴ撮影モードは、歯顎領域の全ての範囲を撮影対象とするＣＴ撮影モードである。全顎ＣＴ撮影モードの画像再構成範囲は例えば直径１６１ｍｍ高さ１００ｍｍの円柱形状の空間領域である。図８は全顎ＣＴ撮影モードの軌道を示している。全顎ＣＴ撮影モードでは、図８に示すように、Ｘ線検出部９の中心がＸ線照射部８と旋回アーム６の旋回軸中心２０６とを結ぶラインの延長線上からずれるように旋回アーム６を旋回させながら複数の撮影位置で撮影が行われる。また、全顎ＣＴ撮影モードでは、通常、図８に示すように、旋回アーム６の旋回軸中心２０６は定位置になっている。なお、図８には撮影位置として４箇所が図示されているが、これはあくまで例示であり撮影位置は図示された箇所に限定されるものではない。
【００３３】
全顎ＣＴ撮影モードは、Ｘ線検出部９の中心をＸ線照射部８と旋回アーム６の旋回軸中心２０６とを結ぶラインの延長線上からずらして撮影を行っているので、上述した全歯ＣＴ撮影モードよりも画像再構成範囲２０７を拡大することができる。したがって、Ｘ線検出部９のサイズアップを抑えながら歯顎領域の全ての範囲を撮影対象とすることができる
【００３４】
なお、全顎ＣＴ撮影モードにおいて、Ｘ線検出部９をサイズアップして、Ｘ線検出部９上でのＸ線ビーム幅Ｗを図８に示す場合よりも拡大し、画像再構成範囲を例えば直径２３０ｍｍ高さ１６４ｍｍの円柱形状の空間領域にすることで、歯顎領域の全ての範囲のみならず、頭頸部領域の全ての範囲を撮影対象とすることも可能である。
【００３５】
上述したパノラマ撮影モード及びＣＴ撮影モードでは、撮影時に患者歯列弓２０３が想定した位置（図２、図３、図７、図８に図示した位置）に存在することで、撮影者が意図していた通りの撮影を行うことができる。患者歯列弓２０３の想定した位置への位置合わせを容易に実現する方法としては、例えば、光ビームを利用する方法を挙げることができる。当該光ビームとしては、例えば、頭の正中線の位置を示す正中線光ビーム、眼窩下縁と外耳道を結ぶ線の位置を示す水平線光ビーム、犬歯の位置（断層撮影の基準位置）を示す断層基準線光ビームなどがあり、これらの光ビームの出力部をＸ線撮影装置に設け、これらの光ビームを参考にして患者が頭の位置を微調整するとよい。
【００３６】
また、旋回アーム６から離れた位置に設置するセファロ用ユニット（不図示）を用い、セファロ撮影モードでの撮影が行えるようにしてもよい。セファロ用ユニットは、被写体を透過したＸ線を検出して、被写体をセファロ撮影するためのセファロ用Ｘ線検出部と、頭部を固定するための頭部固定部とを備える。セファロ撮影は、歯科矯正の診断等に用いられ、頭部規格Ｘ線撮影法（セファロ撮影法）を用いて撮影する。セファロ撮影では、例えば、頭部固定部のイヤーロッドを頭部の左右の外耳孔部に挿入して固定し、旋回アーム６に設けられたＸ線照射部８からＸ線を照射して、被写体を透過したＸ線をセファロ用Ｘ線検出部で検出する。
【００３７】
本発明の一実施形態に係るＸ線撮影装置は、Ｘ線撮影装置の本体部１の他に、図９に示す画像処理装置１０も備えている。
【００３８】
画像処理装置１０は、ＲＯＭ１０２やＨＤＤ１０７に格納されているプログラムに従って画像処理装置１０全体を制御するＣＰＵ１０１と、固定的なプログラムやデータを記録するＲＯＭ１０２と、作業メモリを提供するＲＡＭ１０３と、Ｘ線撮影装置の本体部１内に格納されＸ線撮影装置の本体部１の各部を制御する制御部（不図示）との間で通信を行うための通信インターフェース部１０４と、画像データを一時的に記憶するＶＲＡＭ１０５と、ＶＲＡＭ１０５に記憶された画像データに基づいて画像を表示する表示部１０６と、前記制御部及びＣＰＵ１０１が協働してＸ線撮影動作を制御するための撮影制御プログラム、再構成画像を生成するための画像再構成処理プログラム、散乱線補正処理を行うための散乱線補正処理プログラム等の各種プログラム、各種プログラムを実行する際に用いられる各種パラメータの設定値、並びに、再構成画像データ等の各種データを記憶するＨＤＤ１０７と、キーボード、ポインティングデバイス等の入力部１０８とを備えている。
【００３９】
画像処理装置１０は、画像処理装置１０と前記制御部との通信方法は、有線通信でもよく、無線通信でもよく、有線と無線を組み合わせた通信であってもよい。画像処理装置１０としては、例えば、パーソナルコンピュータを挙げることができる。なお、画像処理装置１０は、画像処理以外に、Ｘ線撮影装置の本体部１の遠隔操作、画像表示も行う。ＨＤＤ１０７に記憶されている各プログラムは、画像処理装置１０にプリインストールされていてもよく、光ディスク等の記憶媒体に格納された形態で流通されて画像処理装置１０にインストールされてもよく、ネットワークを介して流通されて画像処理装置１０にインストールされてもよい。
【００４０】
散乱線補正処理プログラムを実行すると、画像処理装置１０は散乱線補正装置として機能する。散乱線補正処理は画像再構成処理中に割り込んで実施される。
【００４１】
＜＜散乱線補正処理の前提となる理論＞＞
散乱線補正処理の内容を説明する前に、散乱線補正処理の前提となる理論について説明する。当該理論は、測定画像のピクセルの輝度値を元に散乱線による輝度値を算出する理論であって、本発明者が独自に構築したものである。本発明者は、モンテカルロシミュレーションの結果を用いて当該理論を構築した。
【００４２】
ここで、モンテカルロシミュレーションの計算のジオメトリーを図１０に示す。図１０（ａ）は上面図であり、図１０（ｂ）は側面図である。被写体１１は、Ｘ線焦点８ＡとＸ線検出器９Ａとの間に配置される。Ｘ線検出器９Ａはシンチレーター等を備えるフラットパネルディテクターとした。また、Ｘ線検出器９Ａの手前にはカーボン９Ｂが設置されている。
【００４３】
モンテカルロシミュレーションの計算において、Ｘ線焦点８Ａで発生させるＸ線スペクトルは、実際の撮影において使用するＸ線管等の仕様に基づいてＸ線管から放出され、Ｘ線検出器９Ａの位置でＸ線検出器９Ａの有感領域となるようコリメートされた一様なＸ線束となるようにした。
【００４４】
被写体については、生体に近づけるため、顎部を想定した直径１５ｃｍの円柱形の水ファントム１２Ａ、頭部を想定した直径１８ｃｍの円柱形の水ファントム１２Ｂ、頸部を想定した直径１３ｃｍの円柱形の水ファントム１２Ｃを図１１（ａ）に示す側面図および図１１（ｂ）に示す上面図のように組み合せ、水ファントム１２Ａの内部に厚さ２ｍｍの円筒状の皮質骨を設置し、水ファントム１２Ｂの内部に厚さ１ｍｍの円筒状の皮質骨を設置し、水ファントム１２Ａおよび１２Ｃの内部に長径４ｃｍ、短径３ｃｍ、厚さ３ｍｍの楕円筒状の頸椎部分１２Ｄを設置した。ただし、図１１（ａ）に示す一点鎖線で囲った部分においては、円筒状の皮質骨および頸椎部分１２Ｄをカットした。
【００４５】
モンテカルロシミュレーションの計算においては、図１１（ａ）に示す一点鎖線で囲った部分を金属の設置場所として、金属によるＸ線の減衰および散乱線を調べることにする。これは、骨があることによる周りからの散乱線の減衰効果を維持しつつ、骨によらない金属のみによる効果を見積もるためである。
【００４６】
まず、図１１に示す水ファントムの組み合わせであって、図１１（ａ）に示す一点鎖線で囲った部分に金属を設置しないものを被写体にしてモンテカルロシミュレーションの計算を行う。その計算結果を図１２において実線で示す。図１２に示すグラフの横軸は白画像における平均輝度値に対する全Ｘ線（直接線と散乱線）による輝度値の比であり、図１２に示すグラフの縦軸はＸ線検出器９Ａの或るピクセルで検出される全Ｘ線量に対するその全Ｘ線量に含まれる散乱線量の比である。
【００４７】
白画像における平均輝度値に対する全Ｘ線による輝度値の比は被写体のＸ線照射方向の厚みが薄いほど大きくなる。図１２において実線で示されている計算結果より、被写体のＸ線照射方向の厚みが薄いほど全Ｘ線量に対する散乱線量の比が小さくなることが分かる。図１２において実線で示されている計算結果は、例えばｙ＝Ａ／（ｘ−Ｃ）＋Ｂで近似することができる。Ａ＝０．０４９５５６、Ｂ＝０．０４６５０１、Ｃ＝−０．１６２７４とすると、図１２において破線で示されている近似曲線が得られる。
【００４８】
次に、下記（１）〜（３）の場合におけるモンテカルロシミュレーションの各計算結果を図１３及び図１４に示すグラフで比較する。
【００４９】
図１３に示すグラフの横軸は図１１に示すｘ軸方向のピクセル位置を示すピクセル番号であり、図１３に示すグラフの縦軸は白画像における平均輝度値に対する全Ｘ線による輝度値の比である。なお、ピクセル番号が大きいほど、対応するＸ線検出器９Ａのピクセルに入射したＸ線が透過した被写体の厚みは厚い。図１３中の曲線Ｃ１は下記（１）の場合に得られる計算結果である。同様に、図１３中の曲線Ｃ２は下記（２）の場合に得られる計算結果であり、図１３中の曲線Ｃ３は下記（３）の場合に得られる計算結果である。
【００５０】
また、図１４に示すグラフの横軸は図１１に示すｘ軸方向のピクセル位置を示すピクセル番号であり、図１４に示すグラフの縦軸はＸ線検出器９Ａの或るピクセルで検出される全Ｘ線量に対するその全Ｘ線量に含まれる散乱線量の比である。図１４中の曲線Ｃ１１は下記（１）の場合に得られる計算結果である。同様に、図１４中の曲線Ｃ１２は下記（２）の場合に得られる計算結果であり、図１４中の曲線Ｃ１３は下記（３）の場合に得られる計算結果である。
【００５１】
（１）図１１に示す水ファントムの組み合わせであって、図１１に示す一点鎖線で囲った部分に金属を設置しないものを被写体にした場合
（２）図１１に示す水ファントムの組み合わせであって、図１１に示す一点鎖線で囲った部分に幅５ｍｍ、長さ８ｃｍ、厚さ３ｍｍのチタンを、幅方向をｚ軸、長さ方向をｙ軸、厚さ方向をｘ軸に合わせて設置したものを被写体にしてＸ線の照射方向をｘ軸に合わせた場合
（３）図１１に示す水ファントムの組み合わせであって、図１１に示す一点鎖線で囲った部分に幅５ｍｍ、長さ８ｃｍ、厚さ１ｍｍのＡｕＡｇＰｄ合金を、幅方向をｚ軸、長さ方向をｙ軸、厚さ方向をｘ軸に合わせて設置したものを被写体にしてＸ線の照射方向をｘ軸に合わせた場合
【００５２】
図１３及び図１４から、Ｘ線検出器９Ａのピクセルで検出される全Ｘ線による輝度値（測定画像のピクセルの輝度値に相当）が同程度であっても、Ｘ線が透過する金属の種類や厚さによって散乱線量が異なることが分かる。
【００５３】
例えば、白画像の平均輝度値が３万であり、Ｘ線検出器９Ａのピクセルで検出される全Ｘ線による輝度値が６００である条件すなわち白画像における平均輝度値に対する全Ｘ線による輝度値の比が０．０２である条件を満たすピクセル番号は、上記（２）の場合は２０１であり、上記（３）の場合は５である（図１３参照）。これらのピクセル番号に対応するピクセルで検出される全Ｘ線量に対するその全Ｘ線量に含まれる散乱線量の比を図１４からから求めると、それぞれ０．４９１、０．９８７となり、散乱線による輝度値に変換するとそれぞれ２９５、５９２となる。したがって、測定画像のピクセルの輝度値が同じ６００であっても、Ｘ線がどのような材質の物質を透過したかによって、すなわち、物質の減弱係数の違いによって、散乱線量（測定画像のピクセルの散乱線成分）は異なってくる。
【００５４】
図１１に示す水ファントムの組み合わせであって、図１１（ａ）に示す一点鎖線で囲った部分に金属を設置しないものを被写体にしてモンテカルロシミュレーションの計算を行って得られる白画像における平均輝度値に対する全Ｘ線による輝度値の比は、被写体のＸ線照射方向の厚みに応じて異なる（図１３中の曲線Ｃ１参照）。以下、被写体の着目する部位を透過するＸ線経路上の骨あるいは金属を皮膚に置き換えた場合にＸ線検出器９Ａのピクセルで得られる輝度値を理想値と呼び、白画像における平均輝度値に対する理想値の比を理想値比と呼ぶ。
【００５５】
理想値比を固定した状態で金属の種類や厚さを変えながら、白画像における平均輝度値に対する全Ｘ線による輝度値の比と、Ｘ線検出器９Ａの或るピクセルで検出される全Ｘ線量に対するその全Ｘ線量に含まれる散乱線量の比との関係を求めることで図１５に示す一つのグラフを得る。例えば、理想値比を０．１に固定した状態で図１３及び図１４中の白抜き丸、白抜き三角から図１５中の白抜き丸、白抜き三角の各点を得て、その各点からｙ＝Ａ／（ｘ−Ｃ）＋Ｂの曲線を近似してＡ、Ｂ、Ｃの各値を求めることで図１５中の曲線Ｃ２９を得ることができる。
【００５６】
そして、理想値比を変更することで、図１５に示すそれぞれのグラフを得る。図１５中の曲線Ｃ２１は理想値比が１のグラフであり、図１５中の曲線Ｃ２２は理想値比が０．３のグラフであり、図１５中の曲線Ｃ２３は理想値比が０．１のグラフであり、図１５中の曲線Ｃ２４は理想値比が０．０３のグラフである。
【００５７】
ここで、図１５中の各曲線はＸ線撮影の投影角度を或る値に固定した場合に対応するものである。そして、白画像における平均輝度値に対する全Ｘ線による輝度値の比と、Ｘ線検出器９Ａの或るピクセルで検出される全Ｘ線量に対するその全Ｘ線量に含まれる散乱線量の比との関係を示す曲線は、Ｘ線撮影の投影角度に応じて異なる。なぜなら、歯列がＸ線焦点８Ａに近い側に位置するような投影角度のＸ線撮影に対応する場合には、図１１（ａ）に示す一点鎖線で囲った部分内で金属をＸ線焦点８Ａに近い側に設置し、歯列がＸ線焦点８Ａに遠い側に位置するような投影角度のＸ線撮影に対応する場合には、図１１（ａ）に示す一点鎖線で囲った部分内で金属をＸ線焦点８Ａに遠い側に設置するといったように、Ｘ線撮影の投影角度によって被写体中における金属の設置位置が異なり、この金属の設置位置の違いが金属を透過した直接線を検出するＸ線検出器９Ａのピクセルで検出される散乱線量に影響を与えるからである。
【００５８】
以上により、測定画像のピクセルの輝度値を元に散乱線による輝度値を算出する方法をまとめると、以下のようになる。
［１］測定画像の全てのピクセルについて理想値比を算出する。
［２］理想値比と、白画像における平均輝度値に対する測定画像の輝度値の比と、Ｘ線撮影の投影角度とから、測定画像の全てのピクセルの散乱線成分を算出する。
【００５９】
ところで、ＣＴ撮影においては余分な被ばくを抑えるためＸ線の照射範囲がＸ線検出器９Ａの有感領域に一致するようＸ線焦点８Ａの近傍でＸ線をカットしている。よって、Ｘ線検出器９Ａの端部に位置するピクセルでは、Ｘ線検出器９Ａの中央部に位置するピクセルと比較して被写体からの散乱線による影響が少ない。その結果、端部に位置するピクセルでは散乱線成分が減少する。
【００６０】
この減少傾向は、Ｘ線検出器９Ａの端のピクセルほど強いが、Ｘ線検出器９Ａの端から離れるにつれて連続的に小さくなり、Ｘ線検出器９Ａの端から１２０ピクセル以上離れたＸ線検出器９Ａの中央部ではほぼ０となる。しかし、白画像における平均輝度値に対する全Ｘ線による輝度値の比が高いピクセルにおいては、被写体のＸ線照射方向の厚みが薄いので元々被写体からの散乱線は少ない。このため、理想値比が高いピクセルほど、ピクセルがＸ線検出器９Ａの端部に位置したことによる散乱線成分の減少は小さくなる。したがって、ピクセルがＸ線検出器９Ａの端部の方に位置し、且つ、理想値比が低い場合は、上述した散乱線成分の減少を加味して散乱線成分を求めることが望ましい。
【００６１】
＜＜散乱線補正処理の内容＞＞
散乱線補正処理は、理想値比算出処理、散乱線成分算出処理、補正処理に分けられ、理想値比算出処理、散乱線成分算出処理、補正処理の順で実行される。
【００６２】
＜理想値比算出処理＞
画像処理装置１０は、理想値比算出処理において、例えば、５１２×５１２ピクセルの測定画像を１６×１６ピクセルごとにまとめ、横方向３２個×縦方向３２個の合計１０２４個の領域に分割し、各領域の理想値比を計算した後に各ピクセルの理想値比を算出する。このような算出手順にした理由は、被写体のＸ線照射方向の厚みに応じてピクセルの輝度値は徐々に変化するものの、１ピクセルごとに調べていたのでは、各ピクセルの輝度値の誤差が大きいため、輝度値の減少または増加が被写体のＸ線照射方向の厚みが変化したためであるのか、それとも誤差によるものなのか判別するのが非常に難しくなるからである。なお、上記のピクセルサイズや分割する領域の個数はあくまで例示であり、上記に示した値に限定されない。
【００６３】
また、画像処理装置１０は、横方向３２個×縦方向３２個の領域の各中心位置の理想値比について、まず横方向に領域を移動して得られる３２個の曲線（横方向の領域をｘ、白画像における平均輝度値に対する中心位置の全Ｘ線による輝度値の比をｙとして得られる曲線）を測定画像と図１１に示す水ファントムの組み合わせであって、図１１（ａ）に示す一点鎖線で囲った部分に金属を設置しないものを被写体にしてモンテカルロシミュレーションの計算を行って得られるデータとからそれぞれ算出し、次に縦方向に領域を移動して得られる３２個の曲線（縦方向の領域をｘ、白画像における平均輝度値に対する中心位置の全Ｘ線による輝度値の比をｙとして得られる曲線）を測定画像と図１１に示す水ファントムの組み合わせであって、図１１（ａ）に示す一点鎖線で囲った部分に金属を設置しないものを被写体にしてモンテカルロシミュレーションの計算を行って得られるデータとからそれぞれ算出するという手順で理想値比算出処理を行う。なお、本実施形態では、各領域を代表する理想値比としてそれぞれの中心位置の理想値比を採用しているが、各領域を代表する理想値比はこれに限定されることはなく、例えば、各領域において領域全体の理想値比の平均値を求め、各領域を代表する理想値比として各平均値を採用してもよい。
【００６４】
測定画像を対象として横方向に輝度値の高い方から低い方へ行う方向で領域を移動したとき、白画像における平均輝度値に対する中心位置の全Ｘ線による輝度値の比は、骨部分に入るまでの皮膚のみの部分においては基本的にｙ＝Ａ／（ｘ−Ｃ）＋Ｂの曲線に従って変化し、皮膚のみの部分から骨部分に入ると当該曲線から外れて大きく落ち込む。また、当該曲線（ｘの範囲は骨部分に入るまでの皮膚のみの部分に限定）と図１３中の曲線Ｃ１とは必ずしも一致しない。
【００６５】
そこで、画像処理装置１０は、当該曲線（ｘの範囲は骨部分に入るまでの皮膚のみの部分に限定）と図１３中の曲線Ｃ１とが一致するように、図１３中の曲線Ｃ１をｘ軸方向に伸縮処理或いは縮小処理し、ｘ軸方向に伸縮処理或いは縮小処理した後の図１３中の曲線Ｃ１に基づいて領域の中心位置の理想値比を算出する。
【００６６】
しかし、Ｘ線を被写体の背面付近から照射する場合および正面付近から照射する場合については、皮膚のみの部分の範囲が極端に狭くなってしまうため、測定画像を対象として横方向に輝度値の高い方から低い方へ行う方向で領域を移動したときに求まる理想値比の算出精度が悪くなってしまう。
【００６７】
一方、縦方向に関しては、顎部や頸部等の広い範囲で皮膚のみの部分がある。そこで、Ｘ線を被写体の背面付近から照射する場合や正面付近から照射する場合に対応する投影角度においては、横方向に領域を移動させたときの理想値比の計算で皮膚のみの部分がある範囲以下となった時点で、縦方向に領域を移動させたときの理想値比の計算に切り替え、端から例えば５列分について縦方向に領域を移動させたときの理想値比の計算を行うようにする。この縦方向に領域を移動させたときの理想値比の計算結果から得られる値を、皮膚のみの部分から得られる理想値比と仮定することによって、皮膚のみの部分から理想値比を算出することが困難だった横方向に領域を移動させたときの理想値比の算出が可能となるので、その後横方向に領域を移動させたときの理想値比の計算を再開する。
【００６８】
横方向に領域を移動させたときの理想値比の計算は各行で行うため、図１３中の曲線Ｃ１に対するｘ軸方向の伸縮処理或いは縮小処理の程度も各行で異なる。このため、縦方向の理想値比の変化が滑らかではなく、凸凹になることが多い。そこで、縦方向についても理想値比の変化を滑らかにする処置を施す必要がある。ただし、横方向の場合と同様の処置を行ったのでは、次は横方向の理想値比の変化がまた凸凹になってしまうので、本実施形態では、最小二乗法を利用することにした。
【００６９】
ｙ＝Ａ／（ｘ−Ｃ）＋Ｂに最小二乗法を適用すると非常に複雑になり、Ａ、Ｂ、Ｃの値を求めることができない。そこで、理想値比に対して対数z=logyをとると、zは局所的にはxに比例することから、ｚ＝Ａｘ＋Ｂについて最小二乗法を適用し、その後にｚをｙに戻すことで、理想値比を算出するようにした。本実施形態は、最小二乗法を適用する範囲を４つに分割し、それぞれについて境界が滑らかになるようにしつつ別々に最小二乗法を適用した。なお、最小二乗法を適用する範囲の分割数は４つ以外であってもよい。
【００７０】
さらに、画像処理装置１０は、上述した手順で算出した各領域の理想値比の変化をより滑らかにするため、各領域を注目領域の対象とし、注目領域の周りの領域の理想値比の平均値を注目領域の理想値比とする平滑化処理を行う。
【００７１】
理想値比算出処理の最後において、画像処理装置１０は、各領域の中心位置の理想値比に基づいて各ピクセルの理想値比の値を算出する。例えば、理想値比の値の算出対象であるピクセルが、隣り合う領域の中心位置同士を結ぶ線分上に位置する場合は、当該線分上で理想値比が線形的に変化するものとし、理想値比の値の算出対象であるピクセルが、隣り合う領域の中心位置同士を結ぶ線分上に位置しない場合は、当該ピクセルを囲む四つの中心位置の理想値比を利用し、まず当該ピクセルと横方向の座標が同一であって横方向に隣り合う一組の領域の中心位置同士を結ぶ線分上に位置する第１の点の理想値比を当該線分上で理想値比が線形的に変化するものとして算出し、次に当該ピクセルと横方向の座標が同一であって横方向に隣り合うもう一組の領域の中心位置同士を結ぶ線分上に位置する第２の点の理想値比を当該線分上で理想値比が線形的に変化するものとし算出し、そして第１の点と第２の点とを結ぶ線分上に位置する場合は、当該線分上で理想値比が線形的に変化するものとすることで各ピクセルの理想値比の値を算出できるが、他の方法で各ピクセルの理想値比の値を求めてもよい。例えば、上記の方法において横方向を縦方向に置き換えてもよい。
【００７２】
＜散乱線成分算出処理＞
画像処理装置１０は、測定画像のピクセル毎に、上述した理想値比算出処理で算出した理想値比と、白画像における平均輝度値に対する測定画像の輝度値の比と、Ｘ線撮影の投影角度とから、図１５に示す曲線のデータを用いて測定画像の散乱線成分を算出する。図１５に示す曲線のデータは、ｙ＝Ａ／（ｘ−Ｃ）＋Ｂの定数Ａ、Ｂ、Ｃの各値が、理想値比と投影角度とに関連付けられたデータテーブルの形式でＨＤＤ１０７に記憶されている。
【００７３】
ここで、画像処理装置１０は、ピクセルがＸ線検出器９Ａの端部の方に位置し、且つ、理想値比が低い場合は、上記の通り算出した散乱線成分を減少させる修正処理を行うようにしてもよい。
【００７４】
＜補正処理＞
画像処理装置１０は、測定画像のピクセル毎に、測定画像のピクセルの輝度値から、上述した散乱線成分算出処理で算出した散乱線成分を除去する補正処理を行う。これにより、Ｘ線撮影で得られる測定画像から散乱線成分を低減した画像を直接的に導出することができる。
【００７５】
測定画像の輝度値には被写体やフラットパネルディテクターからの散乱線も反映されるので、Ａｕの含有量の多いＡｕＡｇＰｄ合金のようにＸ線を全く通さない金属が被写体であっても、必ず測定画像に輝度値が現れることになる。金属アーチファクト除去操作の過程において、散乱線成分が含まれたままであると金属部分をうまく抽出することができないため、金属アーチファクトを十分に除去することができない。しかし、散乱線成分を取り除くことによりこれら金属部分の透過Ｘ線を正確に算出することができれば、金属アーチファクトを除去できることが期待できる。
【００７６】
また、歯科用のＣＴ撮影では一般に正確なＣＴ値を算出することができない。これは、被写体やフラットパネルディテクターからの散乱線が主な原因であると考えられる。しかし、測定画像から散乱線成分を取り除くことができれば、データとして得られるのは被写体のＸ線減弱によって得られる透過Ｘ線による輝度値のみとなる。したがって、より正確なＣＴ値を算出できることが期待できる。
【００７７】
以上、本発明の一実施形態について説明したが、本発明の範囲はこれに限定されるものではなく、発明の主旨を逸脱しない範囲で種々の変更を加えて実施することができる。
【００７８】
例えば、Ｘ線撮影の撮影範囲が制限されている等の理由により、Ｘ線撮影の投影角度が変わっても白画像における平均輝度値に対する全Ｘ線による輝度値の比と、Ｘ線検出器９Ａの或るピクセルで検出される全Ｘ線量に対するその全Ｘ線量に含まれる散乱線量の比との関係を示す曲線がほとんど変化しない場合には、散乱線成分算出処理において画像処理装置１０が、測定画像のピクセル毎に、上述した理想値比算出処理で算出した理想値比と、白画像における平均輝度値に対する測定画像の輝度値の比との二つのみから、図１５に示す曲線のデータを用いて測定画像の散乱線成分を算出するようにしてもよい。
【００７９】
上述した実施形態において説明した歯科用あるいは耳鼻科用等のＸ線撮影装置では、被写体が人間の頭部に限定されている。ただし、本発明は被写体が人間の頭部である場合に限定されるものではない。被写体が人間の頭部でない場合は、モンテカルロシミュレーションにおいて使用する被写体のモデルを、図１１に示すモデルではなく、実際の被写体に適したモデルにすればよい。
【符号の説明】
【００８０】
１  本発明の一実施形態に係るアーム型Ｘ線撮影装置の本体部
２  ベース
３  下部ポール
４  上部ポール
５  固定アーム
６  旋回アーム
７  頭部保持部
８  Ｘ線照射部
８Ａ  Ｘ線焦点
８Ｂ  Ｘ線絞り
９  Ｘ線検出部
９Ａ  Ｘ線検出器
９Ｂ  カーボン
１０  画像処理装置
１１  被写体
１２Ａ〜１２Ｃ  水ファントム
１２Ｄ  頸椎部分
１０１  ＣＰＵ
１０２  ＲＯＭ
１０３  ＲＡＭ
１０４  通信インターフェース部
１０５  ＶＲＡＭ
１０６  表示部
１０７  ＨＤＤ
１０８  入力部
２０１  仮想歯列弓
２０２  包絡仮想歯列弓
２０３  患者歯列弓
２０４、２０５  仮想線
２０６  旋回アームの旋回軸中心
２０７  画像再構成範囲
Ｐ１  撮影開始位置
Ｐ２  撮影終了位置
Ｗ  Ｘ線検出部上でのＸ線ビーム幅
【図１】
2015065991.tif 000003
【図２】
2015065991.tif 000004
【図３】
2015065991.tif 000005
【図４】
2015065991.tif 000006
【図５】
2015065991.tif 000007
【図６】
2015065991.tif 000008
【図７】
2015065991.tif 000009
【図８】
2015065991.tif 000010
【図９】
2015065991.tif 000011
【図１０】
2015065991.tif 000012
【図１１】
2015065991.tif 000013
【図１２】
2015065991.tif 000014
【図１３】
2015065991.tif 000015
【図１４】
2015065991.tif 000016
【図１５】
2015065991.tif 000017