(57)【要約】【課題】Ｘ線ＣＴ装置による再構成画像において、アーチファクトを抑えつつ、軟部組織領域のノイズを増加させることなく、高コントラスト領域の空間分解能を向上させる。【解決手段】  投影データから画像を生成する過程において、投影データ空間にて、データ値のプロファイルにおけるエッジ領域を強調する処理を行うステップと、投影データ空間または画像空間にて、データ値のコントラストが一定レベル以下の領域に対して、強調を緩和する処理を行うステップとを実行する。例えば、投影データに対して、そのプロファイルにおけるエッジ領域を強調する強調処理と、そのプロファイルにおけるコントラストが一定レベル以下の領域を平滑化する平滑化処理とを行い、上記強調処理および平滑化処理が行われた投影データに基づいて画像を再構成する。【選択図】  図４
2013027520.tif 000002
【特許請求の範囲】
【請求項１】
Ｘ線ＣＴ撮影により得られた撮影対象の投影データから画像を生成する過程において、
投影データ空間にて、データ値のプロファイルにおけるエッジ領域を強調する処理を行うステップと、
前記投影データ空間または画像空間にて、データ値のコントラストが一定レベル以下の領域に対して、強調を緩和する処理を行うステップとを実行する画像生成方法。
【請求項２】
Ｘ線ＣＴ撮影により得られた撮影対象の投影データから画像を生成する過程において、投影データ空間にて、データ値のプロファイルにおけるエッジ領域を強調する処理を行う手段と、
前記過程において、前記投影データ空間または画像空間にて、データ値のコントラストが一定レベル以下の領域に対して、強調を緩和する処理を行う手段とを備えた画像生成装置。
【請求項３】
コンピュータを、
Ｘ線ＣＴ撮影により得られた撮影対象の投影データから画像を生成する過程において、投影データ空間にて、データ値のプロファイルにおけるエッジ領域を強調する処理を行う手段と、
前記過程において、前記投影データ空間または画像空間にて、データ値のコントラストが一定レベル以下の領域に対して、強調を緩和する処理を行う手段として機能させるためのプログラム。
【請求項４】
Ｘ線ＣＴ撮影を行って撮影対象の投影データを収集し、前記投影データから画像を生成するＸ線ＣＴ装置であって、
前記画像を生成する過程において、投影データ空間にて、データ値のプロファイルにおけるエッジ領域を強調する手段と、
前記過程において、前記投影データ空間または画像空間にて、データ値のコントラストが一定レベル以下の領域に対して、強調を緩和する処理を行う手段とを備えているＸ線ＣＴ装置。
【請求項５】
前記投影データのプロファイルにおけるエッジ領域を強調する強調処理を行う強調手段と、
前記プロファイルにおけるコントラストが一定レベル以下の領域を平滑化する平滑化処理を行う平滑化手段と、
前記強調処理および平滑化処理が行われた前記投影データに基づいて画像を再構成する再構成手段とを備えた請求項４に記載のＸ線ＣＴ装置。
【請求項６】
前記強調処理は、エッジの鋭さが所定範囲内であるエッジ領域のみを強調する処理である請求項５に記載のＸ線ＣＴ装置。
【請求項７】
前記強調処理は、前記投影データのプロファイルの高周波成分を表すプロファイルにおけるデータ値が正の第１の閾値以上であり所望の上限値以下である領域と、前記データ値が負の第２の閾値以下であり所望の下限値以上である領域とを強調する処理である請求項６に記載のＸ線ＣＴ装置。
【請求項８】
前記強調処理は、前記投影データのチャネル方向、列方向およびビュー方向の少なくとも一つのプロファイルにおいて行われる請求項５から請求項７のいずれか一項に記載のＸ線ＣＴ装置。
【請求項９】
前記平滑化処理は、前記投影データのプロファイルの高周波成分を表すプロファイルにおけるデータ値が正の第３の閾値以下であり負の第４の閾値以上である領域を平滑化する処理である請求項５から請求項８のいずれか一項に記載のＸ線ＣＴ装置。
【請求項１０】
前記平滑化処理は、前記投影データのチャネル方向、列方向およびビュー方向の少なくとも一つのプロファイルにおいて行われる請求項７から請求項１１のいずれか一項に記載のＸ線ＣＴ装置。
【請求項１１】
前記投影データのプロファイルにおけるエッジ領域を強調する強調処理を行う強調手段と、
前記強調処理が行われた前記投影データに基づいて第１の画像を再構成する第１の再構成手段と、
前記強調処理が行われていない前記投影データに基づいて第２の画像を再構成する第２の再構成手段と、
前記第１の画像または第２の画像において、コントラストが一定レベルを超える第１の領域に前記第１の画像の対応領域の画像を充てており、コントラストが該一定レベル以下である第２の領域に前記第２の画像の対応領域の画像を充てた画像を生成する生成手段とを備えた請求項４に記載のＸ線ＣＴ装置。
【請求項１２】
前記強調処理は、前記投影データの高周波成分を表すプロファイルにおけるデータ値が正の第１の閾値以上である領域と、前記データ値が負の第２の閾値以下である領域とを強調する処理である請求項１１に記載のＸ線ＣＴ装置。
【請求項１３】
前記強調処理は、前記投影データのチャネル方向、列方向およびビュー方向の少なくとも一つのプロファイルにおいて行われる請求項１１または請求項１２に記載のＸ線ＣＴ装置。
【請求項１４】
前記第１の領域は、前記第１の画像または第２の画像における、画素値のバラツキの程度が一定レベルを超える領域であり、
前記第２の領域は、該バラツキの程度が該一定レベル以下である領域である請求項１１から請求項１３のいずれか一項に記載のＸ線ＣＴ装置。
【請求項１５】
前記第１の領域は、前記第１の画像または第２の画像の所望方向の高周波成分を表すプロファイルにおけるデータ値が一定レベルを超える領域であり、
前記第２の領域は、該データ値が該一定レベル以下である領域である請求項１１から請求項１３のいずれか一項に記載のＸ線ＣＴ装置。
【請求項１６】
前記第１および第２の閾値は、前記第１の再構成手段が画像の再構成に用いる再構成関数の種類に応じて変化する請求項７または請求項１２に記載のＸ線ＣＴ装置。
【請求項１７】
前記強調処理は、先鋭化処理である請求項５から請求項１６のいずれか一項に記載のＸ線ＣＴ装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、Ｘ線ＣＴ（Computed Tomography）画像の画質の最適化技術に関する。
【背景技術】
【０００２】
Ｘ線ＣＴ画像の空間分解能を決定する主要因の一つとして、Ｘ線ＣＴ装置のＸ線検出器を構成する検出素子の配列間隔が挙げられる。そのナイキスト（Nyquist）周波数を最大限に引き出して空間分解能を高めようとすると、その高周波成分により、ノイズ（noise）やアーチファクト（artifact）の増大を伴ってしまう。逆に、ノイズやアーチファクトを抑えようとすると、空間分解能が低下してしまう。このように、空間分解能とノイズやアーチファクトの程度とは、通常、トレードオフ（trade off）の関係にある。
【０００３】
そのため、従来のＸ線ＣＴ装置では、再構成される画像の高周波成分を調整することにより、空間分解能とノイズやアーチファクトとのバランス（balance）を取り、実用に耐えるようにしている。すなわち、多くのＸ線ＣＴ装置において、軟部組織領域用である低周波強調の再構成関数や、高コントラスト（contrast）領域用である高周波強調の再構成関数など、種々の周波数強調による再構成関数を提供している（例えば特許文献１，段落［０００５］等参照）。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００６−０６８２２９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
ところで、観察者は、状況により、軟部組織領域と高コントラスト領域の両方を観察したい場合がある。
【０００６】
このような場合には、例えば、軟部組織領域用の再構成関数による画像と、高コントラスト領域用の再構成関数による画像とをそれぞれ再構成する。そして、これら２種類の画像を、モニタ（monitor）に同時に表示させたり、切り換えて表示させたりする。
【０００７】
しかし、この方法では、１つの画像のサイズ（size）が小さくなって見づらくなったり、視線の移動が大きくなったりして、観察作業が煩雑になり、診断効率が悪くなる。
【０００８】
また、さらにこの問題を解決するため、例えば、まず高コントラスト領域用である高周波強調の再構成関数で画像を再構成し、次に、この画像の軟部組織領域だけにノイズ低減処理を行って、軟部組織領域と高コントラスト領域の両方がそれぞれ適当な画質となっている１画像を得る方法も考えられる。
【０００９】
しかし、この方法では、得られた画像において、軟部組織領域のアーチファクトが強調されるなど、副作用が生じやすい。
【００１０】
このような事情により、Ｘ線ＣＴ装置により得られる画像において、アーチファクトを抑えつつ、軟部組織領域のノイズを増加させることなく、高コントラスト領域の空間分解能を向上させることが望まれている。
【課題を解決するための手段】
【００１１】
第１の観点の発明は、Ｘ線ＣＴ撮影により得られた撮影対象の投影データ（data）から画像を生成する過程において、投影データ空間にて、データ値のプロファイル（profile）におけるエッジ（edge）領域を強調する処理を行うステップ（step）と、前記投影データ空間または画像空間にて、データ値のコントラストが一定レベル（level）以下の領域に対して、強調を緩和する処理を行うステップとを実行する画像生成方法を提供する。
【００１２】
第２の観点の発明は、Ｘ線ＣＴ撮影により得られた撮影対象の投影データから画像を生成する過程において、投影データ空間にて、データ値のプロファイルにおけるエッジ領域を強調する処理を行う手段と、前記過程において、前記投影データ空間または画像空間にて、データ値のコントラストが一定レベル以下の領域に対して、強調を緩和する処理を行う手段とを備えた画像生成装置を提供する。
【００１３】
第３の観点の発明は、コンピュータ（computer）を、Ｘ線ＣＴ撮影により得られた撮影対象の投影データから画像を生成する過程において、投影データ空間にて、データ値のプロファイルにおけるエッジ領域を強調する処理を行う手段と、前記過程において、前記投影データ空間または画像空間にて、データ値のコントラストが一定レベル以下の領域に対して、強調を緩和する処理を行う手段として機能させるためのプログラム（program）を提供する。
【００１４】
第４の観点の発明は、Ｘ線ＣＴ撮影を行って撮影対象の投影データを収集し、前記投影データから画像を生成するＸ線ＣＴ装置であって、前記画像を生成する過程において、投影データ空間にて、データ値のプロファイルにおけるエッジ領域を強調する手段と、前記過程において、前記投影データ空間または画像空間にて、データ値のコントラストが一定レベル以下の領域に対して、強調を緩和する処理を行う手段とを備えているＸ線ＣＴ装置を提供する。
【００１５】
第５の観点の発明は、前記投影データのプロファイルにおけるエッジ領域を強調する強調処理を行う強調手段と、前記プロファイルにおけるコントラストが一定レベル以下の領域を平滑化する平滑化処理を行う平滑化手段と、前記強調処理および平滑化処理が行われた前記投影データに基づいて画像を再構成する再構成手段とを備えた上記第４の観点のＸ線ＣＴ装置を提供する。
【００１６】
ここで、「平滑化手段」は、「強調を緩和する処理を行う手段」の一例である。
【００１７】
第６の観点の発明は、前記強調処理が、エッジの鋭さが所定範囲内であるエッジ領域のみを強調する処理である上記第５の観点のＸ線ＣＴ装置を提供する。
【００１８】
第７の観点の発明は、前記強調処理が、前記投影データのプロファイルの高周波成分を表すプロファイルにおけるデータ値が正の第１の閾値以上であり所望の上限値以下である領域と、前記データ値が負の第２の閾値以下であり所望の下限値以上である領域とを強調する処理である上記第６の観点のＸ線ＣＴ装置を提供する。
【００１９】
第８の観点の発明は、前記強調処理が、前記投影データのチャネル（channel）方向、列方向およびビュー（view）方向の少なくとも一つのプロファイルにおいて行われる上記第５の観点から第７の観点のいずれか一つの観点のＸ線ＣＴ装置を提供する。
【００２０】
第９の観点の発明は、前記平滑化処理が、前記投影データのプロファイルの高周波成分を表すプロファイルにおけるデータ値が正の第３の閾値以下であり負の第４の閾値以上である領域を平滑化する処理である上記第５の観点から第８の観点のいずれか一つの観点のＸ線ＣＴ装置を提供する。
【００２１】
第１０の観点の発明は、前記平滑化処理が、前記投影データのチャネル方向、列方向およびビュー方向の少なくとも一つのプロファイルにおいて行われる上記第７の観点から第１１の観点のいずれか一つの観点のＸ線ＣＴ装置を提供する。
【００２２】
第１１の観点の発明は、前記投影データのプロファイルにおけるエッジ領域を強調する強調処理を行う強調手段と、前記強調処理が行われた前記投影データに基づいて第１の画像を再構成する第１の再構成手段と、前記強調処理が行われていない前記投影データに基づいて第２の画像を再構成する第２の再構成手段と、前記第１の画像または第２の画像において、コントラストが一定レベルを超える第１の領域に前記第１の画像の対応領域の画像を充てており、コントラストが該一定レベル以下である第２の領域に前記第２の画像の対応領域の画像を充てた画像を生成する生成手段とを備えた上記第４の観点のＸ線ＣＴ装置を提供する。
【００２３】
ここで、「第１の再構成手段」、「第２の再構成手段」および「生成手段」は、「強調を緩和する処理を行う手段」の一例である。
【００２４】
第１２の観点の発明は、前記強調処理が、前記投影データの高周波成分を表すプロファイルにおけるデータ値が正の第１の閾値以上である領域と、前記データ値が負の第２の閾値以下である領域とを強調する処理である上記第１１の観点のＸ線ＣＴ装置を提供する。
【００２５】
第１３の観点の発明は、前記強調処理が、前記投影データのチャネル方向、列方向およびビュー方向の少なくとも一つのプロファイルにおいて行われる上記第１１の観点または第１２の観点のＸ線ＣＴ装置を提供する。
【００２６】
第１４の観点の発明は、前記第１の領域が、前記第１の画像または第２の画像における、画素値のバラツキの程度が一定レベルを超える領域であり、前記第２の領域が、該バラツキの程度が該一定レベル以下である領域である上記第１１の観点から第１３の観点のいずれか一つの観点のＸ線ＣＴ装置を提供する。
【００２７】
第１５の観点の発明は、前記第１の領域が、前記第１の画像または第２の画像の所望方向の高周波成分を表すプロファイルにおけるデータ値が一定レベルを超える領域であり、前記第２の領域が、該データ値が該一定レベル以下である領域である上記第１１の観点から第１３の観点のいずれか一つの観点のＸ線ＣＴ装置を提供する。
【００２８】
第１６の観点の発明は、前記第１および第２の閾値が、前記第１の再構成手段が画像の再構成に用いる再構成関数の種類に応じて変化する上記第７の観点または第１２の観点のＸ線ＣＴ装置を提供する。
【００２９】
第１７の観点の発明は、前記強調処理が、先鋭化処理である上記第５の観点から第１６の観点のいずれか一つの観点のＸ線ＣＴ装置を提供する。
【発明の効果】
【００３０】
上記観点の発明によれば、エッジ領域の強調を、アーチファクトの発生要因である過補正の調整が難しい再構成関数や画像空間で行うのではなく、過補正の調整が容易な投影データ空間で行い、これに伴う画像上での非高コントラスト領域の不要な強調や、本来的に発生するノイズを抑える処理を、投影データ空間または画像空間で行うことができ、再構成画像において、アーチファクトを抑えつつ、軟部組織領域のノイズを増加させることなく、高コントラスト領域の空間分解能を向上させることができる。
【図面の簡単な説明】
【００３１】
【図１】第一実施形態に係るＸ線ＣＴ装置の構成を概略的に示す図である。
【図２】第一実施形態における中央処理装置の機能ブロック（block）図である。
【図３】第一実施形態に係るＸ線ＣＴ装置における処理の流れを示すフロー（flow）図である。
【図４】第一実施形態に係るＸ線ＣＴ装置における処理の要部を概念的に示す図である。
【図５】投影データ空間におけるエッジ領域強調処理の一例を示す図である。
【図６】投影データ空間における非高コントラスト領域平滑化処理の一例を示す図である。
【図７】第一実施形態による画像と他の画像とにおけるＭＴＦ（modulation transfer function）曲線の比較結果を示す図である。
【図８】第一実施形態による画像と他の画像におけるＭＴＦ値の比較結果を示す図である。
【図９】第一実施形態による画像と他の画像におけるノイズレベルの比較結果を示す図である。
【図１０】第一実施形態による画像と他の画像におけるエリアシングアーチファクト（aliasing artifact）の比較結果を示す図である。
【図１１】第二実施形態における中央処理装置の機能ブロック図である。
【図１２】第一実施形態に係るＸ線ＣＴ装置における処理の流れを示すフロー図である。
【図１３】第二実施形態に係るＸ線ＣＴ装置における処理の要部を概念的に示す図である。
【図１４】第１の画像と第２の画像との合成処理の例を示す図である。
【図１５】第二実施形態による画像と他の画像とにおけるＭＴＦ曲線の比較結果を示す図である。
【図１６】第二実施形態による画像と他の画像におけるＭＴＦ値の比較結果を示す図である。
【図１７】第二実施形態による画像と他の画像におけるノイズレベルの比較結果を示す図である。
【図１８】第二実施形態による画像と他の画像におけるエリアシングアーチファクトの比較結果を示す図である。
【発明を実施するための形態】
【００３２】
以下、発明の実施形態について説明する。なお、これにより発明が限定されるものではない。
【００３３】
（第一実施形態）
図１は、第一実施形態に係るＸ線ＣＴ装置の構成を概略的に示す図である。
【００３４】
Ｘ線ＣＴ装置１００は、操作コンソール（console）１と、撮影テーブル（table）１０と、走査ガントリ（gantry）２０とを具備している。
【００３５】
操作コンソール１は、操作者からの入力を受け付ける入力装置２と、被検体の撮影を行うための各部の制御や画像を生成するためのデータ処理などを行う中央処理装置３と、走査ガントリ２０で取得したデータを収集するデータ収集バッファ（buffer）５と、画像を表示するモニタ（monitor）６と、プログラムやデータなどを記憶する記憶装置７とを具備している。
【００３６】
撮影テーブル１０は、被検体４０を載せて走査ガントリ２０の開口部Ｂに入れ出しするクレードル（cradle）１２を具備している。クレードル１２は、撮影テーブル１０に内蔵するモータ（motor）で昇降および水平直線移動される。なお、ここでは、被検体４０の体軸方向すなわちクレードル１２の水平直線移動方向をｚ方向、鉛直方向をｙ方向、ｚ方向およびｙ方向に垂直な水平方向をｘ方向とする。
【００３７】
走査ガントリ２０は、回転部１５と、回転部１５を回転可能に支持する本体部２０ａとを有する。回転部１５には、Ｘ線管２１と、Ｘ線管２１を制御するＸ線コントローラ（controller）２２と、Ｘ線管２１から発生したＸ線８１をファンビーム（fan
beam）或いはコーンビーム（cone beam）に整形するアパーチャ（aperture）２３と、被検体４０を透過したＸ線８１を検出するＸ線検出器２４と、Ｘ線検出器２４の出力をＸ線投影データに変換して収集するＤＡＳ（Data Acquisition System）（データ収集装置ともいう）２５と、Ｘ線コントローラ２２，アパーチャ２３，ＤＡＳ２５の制御を行う回転部コントローラ２６とが搭載されている。本体部２０ａは、制御信号などを操作コンソール１や撮影テーブル１０と通信する制御コントローラ２９を具備する。回転部１５と本体部２０ａとは、スリップリング（slip ring）３０を介して電気的に接続されている。
【００３８】
Ｘ線管２１およびＸ線検出器２４は、被検体４０が載置される撮影空間、すなわち走査ガントリ２０の空洞部Ｂを挟んで互いに対向して配置されている。回転部１５が回転すると、Ｘ線管２１およびＸ線検出器２４は、その位置関係を維持したまま、被検体４０の周りを回転する。Ｘ線管２１から放射されアパーチャ２３で整形されたファンビーム或いはコーンビームのＸ線８１は、被検体４０を透過し、Ｘ線検出器２４の検出面に照射される。このファンビーム或いはコーンビームのＸ線８１のｘｙ平面における広がり方向をチャネル方向、ｚ方向における広がり方向もしくはｚ方向そのものを列方向という。
【００３９】
Ｘ線検出器２４は、チャネル方向および列方向に配列された複数の検出素子を有している。ＤＡＳ２５は、個々の検出素子で得られた検出信号を個々のチャネルデータに変換し、これらをＸ線投影データとして収集する。
【００４０】
図２に、第一実施形態における中央処理装置の機能ブロック図を示す。図２に示すように、第一実施形態における中央処理装置３は、機能的には、エッジ領域強調部３０１、非高コントラスト領域平滑化部３０２、および画像再構成部３０３を備えている。
【００４１】
エッジ領域強調部３０１は、Ｘ線投影データＰ１に対して、そのプロファイルにおけるエッジ領域を強調するエッジ領域強調処理を行う。
【００４２】
非高コントラスト領域平滑化部３０２は、Ｘ線投影データＰ１に対して、そのプロファイルにおける非高コントラスト領域を平滑化する非高コントラスト領域平滑化処理を行う。
【００４３】
画像再構成部３０３は、エッジ領域強調処理および非高コントラスト領域平滑化処理が行われた処理済Ｘ線投影データＰ２に基づいて、画像ＳＧ１を再構成する。
【００４４】
なお、エッジ領域強調部３０１は、発明における強調手段の一例であり、非高コントラスト領域平滑化部３０２は、発明における平滑化手段の一例であり、画像再構成部３０３は、発明における再構成手段の一例である。
【００４５】
これより、第一実施形態に係るＸ線ＣＴ装置における処理の流れについて説明する。
【００４６】
図３は、第一実施形態に係るＸ線ＣＴ装置における処理の流れを示すフロー図である。また、図４は、当該処理の要部である、投影データ空間でのエッジ領域強調処理および非高コントラス領域平滑化処理を概念的に示す図である。
【００４７】
ステップ（step）Ｓ１では、被検体のスキャン（scan）を行って、複数ビューのＸ線投影データＰ０を収集する。各ビューのＸ線投影データＰ０は、それぞれ、Ｘ線検出器２４を構成する検出素子、すなわちチャネル方向および列方向に配列された個々の検出素子により得られたチャネルデータにより構成されている。なお、収集されたＸ線投影データＰ０には、通常、ビームハードニング補正、検出器チャネル間特性補正、Ｘ線強度補正などを含む前処理が行われるが、ここでは、説明を省略する。
【００４８】
ステップＳ２では、エッジ領域強調部３０１が、収集された各ビューのＸ線投影データに対して、図４（ａ）に示すように、投影データ空間におけるエッジ領域強調処理を行う。エッジ領域強調処理は、例えば、Ｘ線投影データのビュー方向、チャネル方向および列方向の少なくとも一つのプロファイルにおいて行う。
【００４９】
図５に、投影データ空間におけるエッジ領域強調処理の一例を示す。
【００５０】
図５の例では、まず、図５（ａ）に示すような、Ｘ線投影データのプロファイルｆ０にハイパスフィルタ（highpass filter）を適用するなどして、図５（ｂ）に示すような、プロファイルｆ０の高周波成分を表す高周波成分プロファイルｆ１を得る。
【００５１】
次に、図５（ｃ）に示すように、この高周波成分プロファイルｆ１において、データ値が、所望の第１の閾値Ｔｈ１（＞０）以上となる領域と、所望の第２の閾値Ｔｈ２（＜０）以下となる領域とを、強調処理の対象となるエッジ領域ＥＰとして検出する。
【００５２】
ただし、ここでは、高周波成分プロファイルｆ１におけるデータ値が、所望の上限値Ｌｔ１（＞Ｔｈ１＞０）以上である領域と、所望の下限値Ｌｔ２（＜Ｔｈ２＜０）以下である領域とは、強調処理の対象としてのエッジ領域ＥＰからは除外する。このようにする理由は、既に急峻な傾きを持っている鋭いエッジに対してさらに強調処理を行うと、再構成画像におけるＣＴ値のオーバーシュート（overshoot）やアンダーシュート（undershoot）などの過補正を生じさせ、それに伴うアーチファクトを発生させる要因となるからである。なお、上記の上限値Ｌｔ１および下限値Ｌｔ２の適正値は、例えば、シミュレーション（simulation）の結果や、実際の試行結果により求める。
【００５３】
エッジ領域ＥＰが検出されたら、図５（ｄ）に示すように、Ｘ線投影データのプロファイルｆ０において、その検出されたエッジ領域ＥＰを強調する強調処理を行う。この強調処理としては、例えば、ラプラシアンフィルタ（Laplacian filter）やアンシャープマスク（unsharp mask）による先鋭化（鮮鋭化）処理、エッジ領域ＥＰのチャネルデータ値に１より大きい所定の係数を乗算する処理などが考えられる。
【００５４】
このように、投影データ空間でエッジ領域ＥＰの強調を行うと、再構成画像において、高コントラスト領域がより強調され、空間分解能を向上させることができる。また、エッジ領域ＥＰの強調を再構成関数や画像空間で行う場合と比較して、再構成画像におけるＣＴ値のオーバーシュートやアンダーシュートなどの過補正の調整が容易であり、その過補正に伴うアーチファクトの発生を抑えることができる。また、Ｘ線投影データのプロファイルｆ０において、既に鋭いエッジを表すエッジ領域に対しては強調処理を行わないため、再構成画像におけるＣＴ値の過補正を防ぐことができ、それに伴うアーチファクトの発生をさらに抑えることができる。
【００５５】
本例では、このようなエッジ領域ＥＰの検出および強調処理を、Ｘ線投影データのビュー方向、チャネル方向および列方向のそれぞれのプロファイルについて行うものとする。
【００５６】
ところで、このように、画像空間ではなく投影データ空間でエッジ領域ＥＰの強調を行うと、上述の利点がある一方、逆投影などの画像再構成法の特性により、再構成画像において本来高コントラストでない領域さえもが不必要に強調される。すなわち、再構成画像における非高コントラスト領域つまり軟部組織領域においても、ノイズが増大してしまう。
【００５７】
これは、図４（ａ）を参照すると分かるように、エッジ領域ＥＰが強調されたＸ線投影データＰ１を逆投影すると、強調されたエッジ領域ＥＰのデータが、その逆投影されるパス（path）と重なるすべての画像領域に反映され、その画像領域のコントラストが上がり、ノイズが増大するためである。
【００５８】
そこで、次ステップにおいては、このノイズの増大を抑えるため、図４（ｂ）に示すように、Ｘ線投影データのプロファイルｆ０において、エッジ領域強調処理に加えて、非高コントラスト領域平滑化処理を行う。非高コントラスト領域ＮＰが平滑化されたＸ線投影データＰ２を逆投影すると、平滑化された非高コントラスト領域ＮＰのデータが、その逆投影されるパスと重なるすべての画像領域に反映され、その画像領域のコントラストが下がり、ノイズが減少する。多数のビュー方向から同様に逆投影すれば、エッジ領域ＥＰのデータが１回または少ない回数しか逆投影されていない画像領域については、他の方向から逆投影される投影データによって平滑化される。ただし、再構成画像においての本来のエッジ領域は、多数の方向からエッジ領域ＥＰが強調されたデータが逆投影されるので、平滑化された非高コントラスト領域ＮＰのデータに打ち消されることなく、強調されたまま残る。
【００５９】
このように、投影データ空間でエッジ領域ＥＰの強調と非高コントラスト領域ＮＰの平滑化とを組み合わせて行うことにより、再構成画像に展開した場合に、高コントラスト領域での空間分解能を向上させつつ、軟部組織領域でのノイズの増大を防ぐことができる。ちなみに、投影データ空間におけるエッジ領域ＥＰの強調と非高コントラスト領域ＮＰの平滑化の調整次第では、逆にノイズを低下させることも可能であることが分かっている。
【００６０】
ステップＳ３では、非高コントラスト領域平滑化部３０２が、エッジ領域強調処理が行われた各ビューのＸ線投影データＰ１に対して、投影データ空間における非高コントラスト領域平滑化処理を行う。非高コントラスト領域平滑化処理は、例えば、Ｘ線投影データＰ０のビュー方向、チャネル方向および列方向の少なくとも一つのプロファイルにおいて行う。
【００６１】
図６に、投影データ空間における非高コントラスト領域平滑化処理の一例を示す。
【００６２】
図６の例では、まず、図６（ａ）に示すような、Ｘ線投影データのプロファイルｆ０において、図６（ｂ）に示すように、チャネルデータ値のバラツキ程度を表すバラツキ指標値σ１、例えばプロファイル方向の単位幅当たりにおけるチャネルデータ値の標準偏差を算出する。そして、図６（ｃ）に示すように、バラツキ指標値σ１が所望の第３の閾値Ｔｈ３（＞０）以下となる領域を非高コントラスト領域ＮＰとして検出する。
【００６３】
非高コントラスト領域ＮＰが検出されたら、図６（ｄ）に示すように、その検出された非高コントラスト領域ＮＰを平滑化する平滑化処理を行う。この平滑化処理としては、例えば、移動平均フィルタ（平均化フィルタ）やガウシアンフィルタ（Gaussian filter）を用いた平滑化処理が考えられる。
【００６４】
このように、投影データ空間でエッジ領域強調処理と非高コントラスト平滑化処理とを組み合わせて行うことにより、再構成画像に展開した場合の軟部組織領域でのノイズの増大を防ぐことができ、かつ、調整次第でノイズを低下させることも可能になる。
【００６５】
ステップＳ４では、画像再構成部３０３が、エッジ領域強調処理および非高コントラスト領域平滑化処理が行われた処理済Ｘ線投影データＰ２に基づいて画像ＳＧ１を再構成する。再構成には、例えば、三次元フィルタ逆投影法を用いる。
【００６６】
ステップＳ５では、再構成画像ＳＧ１を表示する。
【００６７】
これより、第一実施形態の方法により得られた画像と他の方法により得られた画像との比較結果について説明する。
【００６８】
図７に、ＭＴＦ曲線の比較結果を示す。比較対象は、エッジ領域強調処理および非コントラスト領域平滑化処理を行わないオリジナル（original）画像ＯＧ、単純に再構成関数（recon kernel）をブースト（boost）して高周波成分を強調したカーネルブースト画像ＢＧ、および、エッジ領域強調処理および非高コントラスト領域平滑化処理を行った第一実施形態の方法による画像ＳＧ１である。これらの画像は、いずれも同一のＸ線投影データに基づいて生成されたものである。本比較結果における撮影対象は、空間分解能測定用のファントム（phantom）である。図７に示すグラフの横軸は、空間分解能（line-pair/cm）を表している。同グラフの縦軸は、ＭＴＦ値を表しており、オリジナル画像ＯＧのＭＴＦ値を１００％とした場合の相対値である。
【００６９】
第一実施形態による画像ＳＧ１のＭＴＦ曲線は、カーネルブースト画像ＢＧと同様に、オリジナル画像ＯＧのＭＴＦ曲線より周波数特性がよくなっており、空間分解能が向上しているのが分かる。
【００７０】
図８に、上記の各画像におけるＭＴＦ値の比較結果を示す。本比較結果における撮影対象は、空間分解能測定用のファントムである。図８に示すグラフの縦軸は、ＭＴＦ値の相対値がそれぞれ５％，１０％，５０％になるときの空間分解能（lp/cm）を表しており、オリジナル画像ＯＧにおける上記の各空間分解能を１００％とした場合の相対値である。
【００７１】
この結果から、第一実施形態による画像ＳＧ１のＭＴＦ曲線は、カーネルブースト画像ＢＧのものと同様に、オリジナル画像ＯＧよりもＭＴＦ値を向上させることができることが分かる。
【００７２】
図９に、上記の各画像におけるノイズレベル（ＣＴ値（ＨＵ）のＳＤ値）の比較結果を示す。本比較結果における撮影対象は、空間分解能測定用のファントムである。図９に示すグラフの縦軸は、互いに異なる第１〜第３の関心領域ＲＯＩ１〜ＲＯＩ３についてのノイズレベル（ＣＴ値のＳＤ値）を表しており、オリジナル画像ＯＧにおける上記の各ノイズレベルを１００％とした場合の相対値である。
【００７３】
一般的に、ＭＴＦ値などで表される空間分解能を向上させた場合、その画像のノイズレベルは上昇し、悪化するものである。再構成関数でブーストしたカーネルブースト画像ＢＧの場合は、この例に漏れず、そのノイズレベルが上昇し、悪化している。これに対して、第一実施形態による画像ＳＧ１では、空間分解能を示すＭＴＦ値が向上しているにもかかわらず、そのノイズレベルは小さくなり、改善されている。
【００７４】
このように、第一実施形態によれば、一般的にはトレードオフの関係になる、空間分解能の改善とノイズレベルの低減といった背反する項目を両立させることができるという非常に優れた効果を奏する。
【００７５】
図１０に、上記の各画像におけるエリアシングアーチファクトの比較結果を示す。本比較結果における撮影対象は、人体の頭部である。単純に再構成関数でブーストしたカーネルブースト画像ＢＧの場合には、そのエリアシングアーチファクトが増加してしまうが、第一実施形態の場合には、そのエリアシングアーチファクトは、オリジナル画像と同等か、場合によっては低下して改善されるケースもある。
【００７６】
以上、第一実施形態によれば、エッジ領域の強調を、アーチファクトの発生要因である過補正の調整が難しい再構成関数や画像空間で行うのではなく、過補正の調整が容易な投影データ空間で行い、これに伴う画像上での非高コントラスト領域の不要な強調や、本来的に発生するノイズを抑える処理として、投影データ空間で非高コントラスト領域の平滑化を行うことができ、その結果、Ｘ線ＣＴ装置により得られる画像において、アーチファクトを抑えつつ、軟部組織領域のノイズを増加させることなく、高コントラスト領域の空間分解能を向上させることができる。
【００７７】
また、第一実施形態によれば、強調処理の対象となるエッジ領域の検出に際し、上限値Ｌｔ１，下限値Ｌｔ２を設けているので、エッジ領域の強調による過補正およびそれに伴って発生するアーチファクトの抑制効果を高めることができる。
【００７８】
（第二実施形態）
図１１に、第二実施形態における中央処理装置の機能ブロック図を示す。図１１に示すように、第二実施形態において、中央処理装置３は、機能的には、エッジ領域強調部３１１、第１の画像再構成部３１２、第２の画像再構成部３１３、および画像合成部３１４を備えている。
【００７９】
エッジ領域強調部３１１は、第一実施形態におけるエッジ領域強調部３０１とほぼ同様に、Ｘ線投影データに対して、そのプロファイルにおけるエッジ領域を強調するエッジ領域強調処理を行う。ただし、エッジ領域の検出に際し、上限値Ｌｔ１、下限値Ｌｔ２は、必ずしも設ける必要はない。
【００８０】
第１の画像再構成部３１２は、エッジ領域強調処理が行われた処理済Ｘ線投影データＰ１に基づいて、第１の画像Ｇ１を再構成する。
【００８１】
第２の画像再構成部３１３は、エッジ領域強調処理が行われる前のＸ線投影データＰ０に基づいて、第２の画像Ｇ０を再構成する。
【００８２】
画像合成部３１４は、第１の画像Ｇ１または第２の画像Ｇ０の画像領域において、高コントラスト領域には第１の画像Ｇ１の対応領域の画像を充て、非高コントラスト領域には、第２の画像Ｇ０の対応領域の画像を充てる合成処理を行い、合成画像ＳＧ２を得る。
【００８３】
なお、エッジ領域強調部３１１は、発明における強調手段の一例であり、第１の画像再構成部３１２は、発明における第１の再構成手段の一例であり、第２の画像再構成部３１３は、発明における第２の再構成手段の一例であり、画像合成部３１４は、発明における生成手段の一例である。
【００８４】
これより、第二実施形態に係るＸ線ＣＴ装置における処理の流れについて説明する。
【００８５】
図１２は、第二実施形態に係るＸ線ＣＴ装置における処理の流れを示すフロー図である。また、図１３は、当該処理の要部である、投影データ空間でのエッジ領域強調処理および画像空間での画像合成処理（非高コントラス領域の強調緩和）を概念的に示す図である。
【００８６】
ステップ（step）Ｓ１１では、被検体のスキャン（scan）を行って、複数ビューのＸ線投影データＰ０を収集する。
【００８７】
ステップＳ１２では、エッジ領域強調部３１１が、収集された各ビューのＸ線投影データＰ０に対して、図１３（ａ）に示すように、Ｘ線投影データのプロファイルｆ０におけるエッジ領域ＥＰに対して強調処理を行う。ここでのエッジ領域の強調処理は、第一実施形態と同様に、例えば、Ｘ線投影データのビュー方向、チャネル方向および列方向の少なくとも一つのプロファイルにおいて行う。ただし、第一実施形態では、エッジ領域ＥＰの検出に際しては、上限値Ｌｔ１および下限値Ｌｔ２を設けるなどして、急峻な傾きを持つエッジを強調処理の対象から除外しているが、ここではこのような除外をしなくてもよい。
【００８８】
第一実施形態と同様、このように、投影データ空間でエッジ領域の強調を行うと、再構成画像における高コントラスト領域がより強調され、空間分解能が向上する。また、エッジ領域の強調を再構成関数や画像空間で行う場合と比較して、過補正の調整が容易であり、その過補正に伴うアーチファクトの発生をより抑えることができる。
【００８９】
本例では、このようなエッジ領域の検出および強調処理を、Ｘ線投影データのビュー方向、チャネル方向および列方向のそれぞれのプロファイルについて行うものとする。
【００９０】
ところで、第一実施形態でも説明したが、このように、画像空間ではなく投影データ空間でエッジ領域ＥＰに強調処理を行うと、上述の利点がある一方、逆投影などの画像再構成法の特性により、再構成画像に展開した場合の軟部組織領域、つまり高コントラストではない領域さえも不必要に強調されて、ノイズが増大してしまう。
【００９１】
これは、図１３（ａ）を参照すると分かるように、エッジ領域ＥＰが強調されたＸ線投影データを逆投影すると、強調されたエッジ領域ＥＰのデータが、その逆投影されるパスと重なるすべての画像領域に反映され、その画像領域のコントラストが上がり、ノイズが増大するためである。
【００９２】
そこで、次ステップ以降においては、このノイズの増大を抑えるため、図１３（ｂ）に示すように、エッジ領域の強調処理が行われたＸ線投影データＰ１に基づいて再構成された第１の画像Ｇ１、またはエッジ領域の強調処理が行われていないＸ線投影データＰ０に基づいて再構成された第２の画像Ｇ０において、ＣＴ値のバラツキ等により、高コントラスト領域ＨＭと非高コントラスト領域ＮＭとを求める。そして、高コントラスト領域ＨＭには、エッジ領域の強調が行われた第１の画像Ｇ１の対応領域の画像を充て、それ以外の非高コントラスト領域ＮＭに対しては、エッジ領域強調処理が行われていない第２の画像Ｇ０の対応領域の画像を充てる。これにより、再構成画像に展開した場合の軟部組織領域でのノイズの増大を防ぐことができ、かつ、調整次第でノイズを低下させることも可能になる。
【００９３】
ステップＳ１３では、第１の画像再構成部３１２が、エッジ領域強調処理が行われた処理済Ｘ線投影データＰ１に基づいて第１の画像Ｇ１を再構成する。再構成には、例えば、３次元フィルタ逆投影法を用いる。
【００９４】
ステップＳ１４では、第２の画像再構成部３１３が、エッジ領域強調処理が行われる前の（行われていない）Ｘ線投影データＰ０に基づいて第２の画像Ｇ０を再構成する。再構成には、例えば、３次元フィルタ逆投影法を用いる。
【００９５】
ステップＳ１５では、画像合成部３１４が、第１の画像Ｇ１と第２の画像Ｇ０とを、第１の画像Ｇ１または第２の画像Ｇ０の画像領域において、高コントラスト領域であるか、非高コントラスト領域であるかに応じて合成する。高コントラスト領域ＨＭには、第１の画像Ｇ１の対応領域の画像を充て、非高コントラスト領域ＮＭには、第２の画像Ｇ０の対応領域の画像を充てて、画像の合成処理を行い、合成画像ＳＧ２を得る。
【００９６】
図１４に、第１の画像と第２の画像との合成処理の例を示す。
【００９７】
図１４の例では、まず、第１の画像Ｇ１または第２の画像Ｇ０において、各局所領域ごとに、ＣＴ値のバラツキ程度を示す指標値σ２、例えば単位面積当たりにおけるＣＴ値の標準偏差を算出する。そして、このバラツキ指標値σ２が所望の第４の閾値Ｔｈ４（＞０）以下である領域を非高コントラスト領域ＮＭ、それ以外の領域を高コントラスト領域ＨＭとして検出する。
【００９８】
なお、第１の画像Ｇ１または第２の画像Ｇ０上の一方向におけるＣＴ値のプロファイルにおいて、ＣＴ値のバラツキ程度が一定レベル以下の領域、例えばプロファイルの単位幅当たりにおけるＣＴ値の標準偏差が所望の閾値Ｔｈ４′以下である領域を非高コントラスト領域ＮＭ、それ以外を高コントラスト領域ＨＭとして検出し、さらに、このような検出を、第１の画像Ｇ１または第２の画像Ｇ０上の多方向において行い、高コントラスト領域および非高コントラスト領域を２次元的に検出するようにしてもよい。
【００９９】
高コントラスト領域ＨＭおよび非高コントラスト領域ＮＭが検出されたら、高コントラスト領域ＨＭには、エッジ領域強調処理が行われた第１の画像Ｇ１の対応領域の画像を充て、非高コントラスト領域ＮＭには、エッジ領域強調処理が行われていない第２の画像Ｇ０の対応領域の画像を充てることにより、第１の画像Ｇ１と第２の画像Ｇ０との合成画像ＳＧ２を生成する。これにより、エッジ領域強調処理により画像上の非高コントラスト領域ＮＭに対して現れる不本意な強調を抑制することができる。
【０１００】
このように、投影データ空間でエッジ領域強調処理をすることと、画像空間で非高コントラスト領域の画像を、エッジ領域強調処理を行っていないＸ線投影データに基づく画像に代えることとを組み合わせることにより、再構成画像に展開した場合の軟部組織領域でのノイズの増大を防ぐことができ、かつ、調整次第でノイズを低下させることも可能になる。
【０１０１】
ステップＳ１６では、その生成された合成画像ＳＧ２を表示する。
【０１０２】
これより、第二実施形態の方法により得られた画像と他の方法により得られた画像との比較結果について説明する。
【０１０３】
図１５に、ＭＴＦ曲線の比較結果を示す。比較対象は、第一実施形態と同様に、オリジナル画像ＯＧ、カーネルブースト画像ＢＧ、および、第二実施形態の方法による画像ＳＧ２である。これらの画像は、いずれも同一のＸ線投影データに基づいて生成されたものである。本比較結果における撮影対象は、空間分解能測定用のファントムである。図１５に示すグラフの横軸は、空間分解能（lp/cm）を表している。同グラフの縦軸は、ＭＴＦ値を表しており、オリジナル画像ＯＧのＭＴＦ値を１００％とした場合の相対値である。
【０１０４】
第二実施形態による画像ＳＧ２のＭＴＦ曲線は、カーネルブースト画像ＢＧと同様に、オリジナル画像ＯＧのＭＴＦ曲線より周波数特性がよくなっており、空間分解能が向上しているのが分かる。
【０１０５】
図１６に、上記の各画像におけるＭＴＦ値の比較結果を示す。本比較結果における撮影対象は、空間分解能測定用のファントムである。図１６に示すグラフの縦軸は、ＭＴＦ値の相対値がそれぞれ５％，１０％，５０％になるときの空間分解能（lp/cm）を表しており、オリジナル画像ＯＧにおける上記の各空間分解能を１００％とした場合の相対値である。
【０１０６】
この結果から、第二実施形態による画像ＳＧ２のＭＴＦ曲線は、カーネルブースト画像ＢＧのものと同様に、オリジナル画像ＯＧよりもＭＴＦ値を向上させることができることが分かる。
【０１０７】
図１７に、上記の各画像におけるノイズレベル（ＣＴ値（ＨＵ）のＳＤ値）の比較結果を示す。本比較結果における撮影対象は、空間分解能測定用のファントムである。図１７に示すグラフの縦軸は、互いに異なる第１〜第３の関心領域ＲＯＩ１〜ＲＯＩ３についてのノイズレベル（ＣＴ値のＳＤ値）を表しており、オリジナル画像ＯＧにおける上記の各ノイズレベルを１００％とした場合の相対値である。
【０１０８】
カーネルブースト画像ＢＧの場合は、ノイズレベルが上昇し、悪化している。これに対して、第二実施形態による画像ＳＧ２では、空間分解能を示すＭＴＦ値が向上しているにもかかわらず、そのノイズレベルは維持されており、カーネルブースト画像ＢＧとの比較では、より小さくなっており、改善されている。
【０１０９】
このように、第二実施形態によれば、第一実施形態と同様に、一般的にはトレードオフの関係になる、空間分解能の改善とノイズレベルの低減といった背反する項目を両立させることができる。
【０１１０】
図１８に、上記の各画像におけるエリアシングアーチファクトの比較結果を示す。本比較結果における撮影対象は、人体の頭部である。単純に再構成関数でブーストしたカーネルブースト画像ＢＧの場合には、そのエリアシングアーチファクトが増加してしまうが、第二実施形態による画像ＳＧ２の場合には、第一実施形態と同様に、そのエリアシングアーチファクトは、オリジナル画像ＯＧと同等である。
【０１１１】
以上、第二実施形態によれば、エッジ領域の強調を、アーチファクトの発生要因である過補正の調整が難しい再構成関数や画像空間で行うのではなく、過補正の調整が容易な投影データ空間で行い、これに伴う画像上での非高コントラスト領域の不要な強調や、本来的に発生するノイズを抑える処理として、画像上での非高コントラスト領域に、エッジ領域の強調を行っていない画像を充てて、画像空間での平滑化に相当する処理を行うことができ、その結果、Ｘ線ＣＴ装置により得られる画像において、アーチファクトを抑えつつ、軟部組織領域のノイズを増加させることなく、高コントラスト領域の空間分解能を向上させることができる。
【０１１２】
なお、発明は、これらの実施形態に限定されず、その趣旨を逸脱しない範囲において種々の変更が可能である。
【０１１３】
例えば、第一実施形態において、強調処理の対象となるエッジ領域の検出に際し、閾値Ｔｈ１，Ｔｈ２，Ｔｈ３、上限値Ｌｔ１、下限値Ｌｔ２の少なくともいずれかを、画像再構成に用いる再構成関数（kernel）に応じて調整してもよい。これにより、再構成関数で強調あるいは抑制される、画像の周波数成分を考慮して、ノイズ低減やアーチファクト抑制、空間分解能の向上といった複数のファクタの総合的なバランスを最適化することができる。
【０１１４】
また例えば、第二実施形態において、強調処理の対象となるエッジ領域の検出に際し、第一実施形態と同様、上限値Ｌｔ１，下限値Ｌｔ２を設けてもよい。この場合、エッジ領域の強調による過補正およびそれに伴って発生するアーチファクトの更なる抑制が期待できる。また、閾値Ｔｈ１，Ｔｈ２，Ｔｈ４（Ｔｈ４′），上限値Ｌｔ１，下限値Ｌｔ２の少なくともいずれかを、画像再構成に用いる再構成関数に応じて調整してもよい。これにより、再構成関数で強調あるいは抑制される、画像の周波数成分を考慮して、ノイズ低減やアーチファクト抑制、空間分解能の向上といった複数のファクタ（factor）の総合的なバランスを最適化することができる。
【０１１５】
また例えば、第一および第二実施形態は、Ｘ線ＣＴ装置であるが、発明は、Ｘ線ＣＴ装置とＰＥＴまたはＳＰＥＣＴとを組み合わせたＰＥＴ−ＣＴ装置やＳＰＥＣＴ−ＣＴ装置などの他、撮影手段を有していない画像生成装置などにも適用可能である。また、コンピュータを、上述の画像生成機能を備えた画像生成装置として機能させるためのプログラム、当該プログラムが記憶された記憶媒体もまた、発明の実施形態の一例である。
【符号の説明】
【０１１６】
１  操作コンソール
２  入力装置
３  中央処理装置
５  データ収集バッファ
６  モニタ
７  記憶装置
１０  撮影テーブル
１２  クレードル
１５  回転部
２０  走査ガントリ
２１  Ｘ線管
２２  Ｘ線コントローラ
２３  アパーチャ
２４  Ｘ線検出器
２５  ＤＡＳ
２６  回転部コントローラ
２９  制御コントローラ
３０  スリップリング
４０  被検体
８１  Ｘ線
１００  Ｘ線ＣＴ装置
３０１  エッジ領域強調部
３０２  非高コントラスト領域平滑化部
３０３  画像再構成部
３１１  エッジ領域強調部
３１２  第１の画像再構成部
３１３  第２の画像再構成部
３１４  画像合成部
【図１】
2013027520.tif 000003
【図２】
2013027520.tif 000004
【図３】
2013027520.tif 000005
【図４】
2013027520.tif 000006
【図５】
2013027520.tif 000007
【図６】
2013027520.tif 000008
【図７】
2013027520.tif 000009
【図８】
2013027520.tif 000010
【図１１】
2013027520.tif 000011
【図１２】
2013027520.tif 000012
【図１３】
2013027520.tif 000013
【図１４】
2013027520.tif 000014
【図１５】
2013027520.tif 000015
【図１６】
2013027520.tif 000016
【図９】
2013027520.tif 000017
【図１０】
2013027520.tif 000018
【図１７】
2013027520.tif 000019
【図１８】
2013027520.tif 000020