(57)【要約】【課題】マーカを用いた治療において、治療直前の患者体内の構造を従来に比べてより効率良く確認することができ、患者を計画位置に設置するのに要する時間を短縮することが可能な放射線治療装置を提供する。【解決手段】粒子線治療装置１は、治療計画の作成時に撮像した治療計画ＣＴ画像の中のマーカ６１の位置と粒子線の照射直前ＣＴ画像の中のマーカ６１の位置との位置関係、および治療計画ＣＴ画像の中の標的４７の位置と照射直前ＣＴ画像の中の標的４７の位置との位置関係を算出する算出装置３１６Ｂと、算出装置３１６Ｂで算出された治療計画ＣＴ画像と照射直前ＣＴ画像における標的４７およびマーカ６１の位置関係を表示する表示装置３１５と、を備えた。【選択図】  図９
2019072393.tif 000002
【特許請求の範囲】
【請求項１】
標的に対して放射線を照射する放射線治療装置であって、
治療計画の作成時に撮像した治療計画ＣＴ画像の中のマーカの位置と前記放射線の照射直前ＣＴ画像の中のマーカの位置との位置関係、および前記治療計画ＣＴ画像の中の標的の位置と前記照射直前ＣＴ画像の中の標的の位置との位置関係を算出する算出装置と、
前記算出装置で算出された前記治療計画ＣＴ画像と前記照射直前ＣＴ画像における標的およびマーカの位置関係を表示する表示装置と、を備えた
ことを特徴とする放射線治療装置。
【請求項２】
請求項１に記載の放射線治療装置において、
前記算出装置は、前記マーカを基準として、前記標的の重心の位置関係を算出する
ことを特徴とする放射線治療装置。
【請求項３】
請求項１に記載の放射線治療装置において、
前記位置関係の表示では、前記治療計画ＣＴ画像と前記照射直前ＣＴ画像におけるマーカの位置を一致させて、前記治療計画ＣＴ画像と前記照射直前ＣＴ画像のそれぞれに対応する前記標的の形状を重ねて表示する
ことを特徴とする放射線治療装置。
【請求項４】
請求項３に記載の放射線治療装置において、
前記位置関係の表示では、前記標的の形状の差分を強調表示する
ことを特徴とする放射線治療装置。
【請求項５】
請求項３に記載の放射線治療装置において、
前記位置関係の表示では、前記標的の形状の差分を予め定めた閾値と比較し、前記差分が前記閾値を超えると判定されたときは、前記閾値を超えた差分を強調表示する
ことを特徴とする放射線治療装置。
【請求項６】
請求項４または５に記載の放射線治療装置において、
前記強調表示では、画面上に文字を表示するとともに、色を変えて表示する
ことを特徴とする放射線治療装置。
【請求項７】
請求項１に記載の放射線治療装置において、
前記算出装置は、前記治療計画ＣＴ画像と前記照射直前ＣＴ画像におけるマーカの位置を一致させて、前記治療計画ＣＴ画像と前記照射直前ＣＴ画像のそれぞれにおける前記標的や正常臓器の照射方向の手前側および／または奥側の境界までの水等価厚とその差分を求め、
前記表示装置では、前記治療計画ＣＴ画像と前記照射直前ＣＴ画像の前記水等価厚の差分を表示する
ことを特徴とする放射線治療装置。
【請求項８】
請求項１に記載の放射線治療装置において、
前記算出装置は、前記マーカと前記標的との位置関係に基づいて、前記標的を有する患者を載せるベッドの移動量を算出する
ことを特徴とする放射線治療装置。
【請求項９】
請求項１に記載の放射線治療装置において、
前記マーカおよび前記標的をＸ線により複数の方向から撮像するＸ線透視装置を更に備え、
前記算出装置は、前記Ｘ線透視装置によって取得したＸ線画像から画像再構成を行い、前記照射直前ＣＴ画像としてコーンビームＣＴ画像を生成する
ことを特徴とする放射線治療装置。
【請求項１０】
請求項９に記載の放射線治療装置において、
前記Ｘ線透視装置によって取得したＸ線画像から患者内の前記標的の３次元位置を認識する認識装置を更に備え、
前記算出装置は、前記Ｘ線透視装置によって取得したＸ線画像の中から、前記認識装置によって認識した標的の位置が動体追跡治療の治療放射線照射条件を満たす画像から画像再構成を行う
ことを特徴とする放射線治療装置。
【請求項１１】
請求項９に記載の放射線治療装置において、
患者内の標的の位置が動体追跡治療の治療放射線照射条件を満たすときに前記標的に放射線を照射するよう制御する治療放射線制御装置を更に備えた
ことを特徴とする放射線治療装置。
【請求項１２】
請求項１１に記載の放射線治療装置において、
前記治療放射線制御装置は、前記マーカと前記標的との位置関係に基づいて、前記治療放射線照射条件を変更する
ことを特徴とする放射線治療装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、陽子線や炭素線等の粒子線またはＸ線等の放射線を患部に照射する放射線治療装置に関する。
【背景技術】
【０００２】
放射線治療の患者位置決めにおいて必要とされる、治療放射線照射状態での標的や放射線通過領域、及び危険臓器の位置関係を把握可能な画像を取得することができる放射線治療システムの一例として、特許文献１には、回転しながら複数の方向から被検体にＸ線を照射しＸ線撮影するＸ線撮影装置と、Ｘ線撮影装置によって取得したＸ線画像から被検体内の標的の３次元位置を認識する標的認識装置と、Ｘ線撮影装置によって取得したＸ線画像の中から、標的認識装置によって認識した標的の位置が動体追跡治療の治療放射線照射条件を満たす画像を選択して画像再構成を行い、コーンビームＣＴ画像を生成するＣＴ画像生成装置と、を備える放射線治療システムが記載されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１５−２９７９３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
癌などの患者の患部に陽子線などの放射線を照射する治療方法が知られている。この治療方法に用いる放射線治療装置は、例えば陽子線の場合、加速器及びビーム輸送系からなるビーム発生装置、照射野形成装置、及びＸ線透視装置と患者ベッドからなる位置決め装置を備えている。
【０００５】
加速器で加速された粒子線は、ビーム輸送系を経て照射野形成装置に達し、この照射野形成装置によってモニタされ、かつ患者の患部の形に合うよう整形される。
【０００６】
陽子線や炭素線等の粒子線は、停止する直前にエネルギーの大部分を放出する特性を有しており、その結果得られる線量分布の形はブラッグピークと呼ばれている。粒子線治療装置は、この特性を利用し、粒子線のエネルギーを選択することで粒子線を照射目標で停止させてエネルギーの大部分を患部に放出する。
【０００７】
放射線治療装置を用いて、患者の患部に向けて陽子線を照射し、患部に合致した形状の線量分布を形成するためには、患者を計画通りの位置に設置することが重要である。
【０００８】
患者の位置を計画した位置に精度良く設置するため、治療直前の標的と標的周囲の位置や形状を確認する方法として、上述の特許文献１に記載されたようなコーンビームＣＴ装置により３次元Ｘ線ＣＴ画像（以下、コーンビームＣＴ画像と記載）を取得する方法がある。
【０００９】
これは、計画用に取得した３次元Ｘ線ＣＴ画像（以下、計画ＣＴ画像と記載）上には、アイソセンタと呼ばれる照射の基準点が指定されており、コーンビームＣＴ画像と計画ＣＴ画像が一致するように患者の位置を動かすことで、患者を位置精度良く設置することができるものである。
【００１０】
放射線治療では、呼吸などで動く臓器を照射する場合に標的に集中して陽子線を照射するため、患部などの標的近傍に設置したマーカが、計画した位置に来たときのみ陽子線を照射する方法がある。
【００１１】
コーンビームＣＴ画像の撮影では、患者の周りを１ｒｐｍ程度の速度でＸ線透視装置が回転する。ＣＴ画像の再構成には１８０度分以上の方向からのＸ線画像が必要となるため、撮影時間は３０秒以上が必要である。
【００１２】
特に、呼吸などで動く臓器の撮像では、標的が撮像中に動くため、特許文献１では、標的またはマーカの位置を計測しながらＸ線画像を取得し、標的またはマーカが計画位置にあるときのＸ線画像のみを使用してＣＴ画像を作成する方法が開示されている。
【００１３】
以上のようにして取得したコーンビームＣＴ画像は、計画ＣＴ画像と比較して表示することで患者の設置位置を決めるために使用されている。
【００１４】
標的近傍に設置したマーカが計画位置にきたときのみ照射する場合、コーンビームＣＴ画像を使用することにより、マーカの位置と標的の位置の両方を確認することができる。
【００１５】
ここで、計画ＣＴ画像におけるマーカと標的の位置とコーンビームＣＴにおけるマーカと標的の位置の比較，確認は、従来は、コーンビームＣＴ画像と計画ＣＴ画像を重ね合わせて表示することで、操作者がその違いを目視により確認していた。
【００１６】
本発明は、マーカを用いた治療において、治療直前の患者体内の構造を従来に比べて効率良く確認することができ、患者を計画位置に設置するのに要する時間を短縮することが可能な放射線治療装置を提供する。
【課題を解決するための手段】
【００１７】
本発明は、上記課題を解決する手段を複数含んでいるが、その一例を挙げるならば、標的に対して放射線を照射する放射線治療装置であって、治療計画の作成時に撮像した治療計画ＣＴ画像の中のマーカの位置と前記放射線の照射直前ＣＴ画像の中のマーカの位置との位置関係、および前記治療計画ＣＴ画像の中の標的の位置と前記照射直前ＣＴ画像の中の標的の位置との位置関係を算出する算出装置と、前記算出装置で算出された前記治療計画ＣＴ画像と前記照射直前ＣＴ画像における標的およびマーカの位置関係を表示する表示装置と、を備えたことを特徴とする。
【発明の効果】
【００１８】
本発明によれば、マーカと標的の位置の確認における操作者の操作量を従来に比べて削減することができるため、治療直前の患者体内の構造を従来に比べて効率良く確認することができ、患者を計画位置に設置するのに要する時間を短縮することができる。
【図面の簡単な説明】
【００１９】
【図１】本発明の一実施例の粒子線治療装置の全体構成の概略を示す図である。
【図２】本実施例の粒子線治療装置における照射野形成装置の概略を示す図である。
【図３】本実施例の粒子線治療装置による粒子線照射での最大エネルギーのスポット配置を示す図である。
【図４】本実施例の粒子線治療装置による粒子線照射での２番目のエネルギーのスポット配置を示す図である。
【図５】本実施例の粒子線治療装置における位置決めシステムの概略を示す図である。
【図６】本実施例の粒子線治療装置による治療計画を立案する流れを説明するフロー図である。
【図７】本実施例の粒子線治療装置による粒子線を照射する流れを説明するフロー図である。
【図８】本実施例の粒子線治療装置における位置決めの流れを説明するフロー図である。
【図９】本実施例の粒子線治療装置におけるベクトルの表示画面の一例を示す図である。
【図１０】本実施例の粒子線治療装置における計画ＣＴとコーンビーム画像での標的とを重ねて表示する表示画面の一例を示す図である。
【図１１】本実施例の粒子線治療装置における計画ＣＴとコーンビーム画像での標的との差分を強調表示する表示画面の一例を示す図である。
【図１２】本実施例の粒子線治療装置における水等価厚を表示する表示画面の一例を示す図である。
【発明を実施するための形態】
【００２０】
本発明の好適な一実施例である粒子線治療装置について、図１乃至図１２を用いて説明する。
【００２１】
本実施例では、放射線治療装置の一種として陽子線や炭素線等の粒子線を用いる粒子線治療装置を例に説明するが、Ｘ線を用いるＸ線治療装置に対しても本発明を適用することで同様の効果が得られる。
【００２２】
最初に、図１を用いて粒子線治療装置の全体構成について説明する。図１は粒子線治療装置の全体構成の概略を示す図である。
【００２３】
図１において、標的４７に対して粒子線を照射する粒子線治療装置１は、荷電粒子ビーム発生装置１１、高エネルギービーム輸送系２０、回転装置２１、中央制御装置３１２、メモリ３１３、照射制御システム３１４、表示装置３１５、位置決め制御システム３１６、照射野形成装置（照射装置）４０、ベッド５０、Ｘ線透視装置５１、治療計画装置５０１を備えている。
【００２４】
荷電粒子ビーム発生装置１１は、イオン源１２、前段加速器１３、粒子ビーム加速装置１４から構成される。本実施例では、粒子ビーム加速装置１４としてシンクロトロン型の粒子ビーム加速装置を想定しているが、粒子ビーム加速装置１４としてサイクロトロン等、他の様々な粒子ビーム加速装置を用いることができる。
【００２５】
シンクロトロン型の粒子ビーム加速装置１４は、図１に示すように、その周回軌道上に偏向電磁石１５、加速装置１６、出射用の高周波印加装置１７、出射用デフレクタ１８、および４極電磁石（図示省略）を備える。
【００２６】
図１を用いて、粒子ビームが、シンクロトロン型の粒子ビーム加速装置１４を利用した荷電粒子ビーム発生装置１１から発生し、患者へ向けて出射されるまでの経過を説明する。
【００２７】
イオン源１２より供給された粒子は、前段加速器１３にて加速され、ビーム加速装置であるシンクロトロンへと送られる。シンクロトロン内を周回する粒子ビームは、加速装置１６を通過する周期に同期させて加速装置１６に設けられた高周波加速空胴（図示省略）から高周波が印加され、加速される。このようにして粒子ビームが所定のエネルギーに達するまで加速される。
【００２８】
所定のエネルギー（例えば７０〜２５０ＭｅＶ）まで粒子ビームが加速された後、中央制御装置３１２より、照射制御システム３１４を介して出射開始信号が出力されると、高周波印加装置１７に設置された高周波印加電極から高周波電源１９からの高周波電力がシンクロトロン内を周回している粒子ビームに印加され、粒子ビームがシンクロトロンから出射される。
【００２９】
高エネルギービーム輸送系２０は、シンクロトロンと照射野形成装置４０とを連絡しており、取り出された粒子ビームは、高エネルギービーム輸送系２０を介して回転装置２１に設置された照射野形成装置４０まで導かれる。
【００３０】
回転装置２１は、患者３０の任意の方向からビームを照射するためにあって、ガントリーと呼ばれる筒状の構造物に搭載させており、装置全体が回転することで患者３０の設置されたベッド５０の周囲どの方向へも回転することができる。
【００３１】
照射野形成装置４０は、最終的に患者３０へ照射する粒子ビームの形状を整形する装置であり、その構造は照射方式により異なる。散乱体法とスキャニング法が、代表的な照射方式であり、本発明はどちらの照射方式でも有効である。本実施例はスキャニング法を用いて説明する。
【００３２】
スキャニング法は、高エネルギービーム輸送系２０から輸送された細いビームをそのまま標的へ照射し、これを３次元的に走査することで、最終的に標的のみに高線量領域を形成することができる。
【００３３】
Ｘ線透視装置５１とベッド５０は位置決め制御システム３１６を介して表示装置３１５と中央制御装置３１２と繋がっている。
【００３４】
Ｘ線透視装置５１は、それぞれ２台のＸ線発生装置とＸ線検出装置から構成されており、患者３０の体内を透視したＸ線透視画像を二つの直交する方向から撮影することができる。このＸ線透視画像は、患者３０内のマーカ６１（図２参照）の位置を検出するためにマーカ６１および標的４７（図２参照）をＸ線により複数の方向から撮像するものであり、粒子線照射前の位置決めと、粒子線照射中の標的位置の計測に用いられる。
【００３５】
位置決め制御システム３１６は、Ｘ線透視装置５１によって撮像されたＸ線透視画像から、粒子ビーム照射時のベッド５０の位置を決定するための処理を実行するシステムであり、認識装置３１６Ａおよび算出装置３１６Ｂを備えている。
【００３６】
認識装置３１６Ａは、Ｘ線透視装置５１によって取得したＸ線画像から患者３０内の標的４７の３次元位置を認識する。
【００３７】
算出装置３１６Ｂは、治療計画の作成時に３次元撮像装置７００によって撮像され、データサーバ６００に保存されている治療計画ＣＴ画像の中のマーカ６１の位置と粒子線の照射直前ＣＴ画像の中のマーカ６１の位置との位置関係を算出する。また、算出装置３１６Ｂは、治療計画ＣＴ画像の中の標的４７の位置と照射直前ＣＴ画像の中の標的４７の位置との位置関係を算出する。この際、マーカ６１の位置を基準として、標的４７の重心の位置関係を算出する。
【００３８】
算出装置３１６Ｂは、算出した位置関係に基づいて、標的４７を有する患者３０を載せるベッド５０の移動量を算出する。
【００３９】
マーカ６１の位置関係や標的４７の位置関係を算出するために、算出装置３１６Ｂは、Ｘ線透視装置５１によって取得したＸ線画像から画像再構成を行い、照射直前ＣＴ画像としてコーンビームＣＴ画像を生成する。この際、算出装置３１６Ｂは、Ｘ線透視装置５１によって取得したＸ線画像の中から、認識装置３１６Ａによって認識した標的４７の位置が所定の範囲内（照射許可領域、動体追跡治療の治療粒子線照射条件）にあるときに撮像されたＸ線画像から画像再構成を行う。
【００４０】
更には、算出装置３１６Ｂは、治療計画ＣＴ画像と照射直前ＣＴ画像におけるマーカ６１の位置を一致させて、その際の治療計画ＣＴ画像と照射直前ＣＴ画像のそれぞれにおける標的４７や正常臓器の照射方向の手前側および／または奥側の境界までの水等価厚とその差分を求める。
【００４１】
表示装置３１５は、位置決め制御システム３１６の算出装置３１６Ｂで算出された治療計画ＣＴ画像と照射直前ＣＴ画像における標的４７およびマーカ６１の位置関係を表示する。この表示装置３１５では、算出装置３１６Ｂで標的４７や正常臓器の照射方向の手前側および／または奥側の境界までの水等価厚とその差分を求めた場合には、その水等価厚とその差分を表示することが可能である。
【００４２】
以上の認識装置３１６Ａおよび算出装置３１６Ｂの処理の詳細や、表示装置３１５で表示される表示画面の詳細は後述する。
【００４３】
図２を使って、照射野形成装置４０内の機器のそれぞれの役割と機能とを簡単に述べる。図２は、スキャニング法に対応した照射野形成装置４０の構成を示す図である。
【００４４】
照射野形成装置４０は、上流側から二つの走査電磁石４１，４２、線量モニタ４３、ビーム位置モニタ４４を備える。線量モニタ４３はモニタを通過した粒子ビームの量を計測する。ビーム位置モニタ４４は、粒子ビームが通過した位置を計測することができる。これら線量モニタ４３およびビーム位置モニタ４４からの情報により、計画通りの位置に、計画通りの量のビームが照射されていることを、照射制御システム３１４が管理することが可能となる。
【００４５】
荷電粒子ビーム発生装置１１から高エネルギービーム輸送系２０を経て輸送された細い粒子ビームは、走査電磁石４１，４２によりその進行方向が偏向される。これらの走査電磁石４１，４２は、ビーム進行方向と垂直な方向に磁力線が生じるように設けられており、例えば図２では、走査電磁石４１は走査方向４５の方向にビームを偏向させ、走査電磁石４２は走査方向４５に垂直な方向にビームを偏向させる。これらの二つの電磁石を利用することで、ビーム進行方向と垂直な面内において任意の位置にビームを移動させることができ、標的４７へのビーム照射が可能となる。
【００４６】
ここで、照射制御システム３１４は、患者３０内の標的４７の位置が照射許可領域内に存在するときに標的４７に粒子線を照射する制御、いわゆる動体追跡制御を実行する。この照射制御システム３１４は、マーカ６１と標的４７との位置関係から、照射許可領域を変更することが可能である。
【００４７】
また、照射制御システム３１４は、走査電磁石磁場強度制御装置４９を介して、走査電磁石４１，４２に流す電流の量を制御する。走査電磁石４１，４２には、走査電磁石用電源４８より電流が供給され、電流量に応じた磁場が励起されることでビームの偏向量を自由に設定できる。粒子ビームの偏向量と電流量との関係は、あらかじめテーブルとして中央制御装置３１２の中のメモリ３１３に保持されており、それを参照する。
【００４８】
スキャニング法のビームの走査方式は二通りある。一つは照射位置の移動と停止を繰り返す離散的な方式、もう一つは連続的に照射位置を変化させる方式である。
【００４９】
離散的な方式では、まず、照射位置をある点に留めたまま、規定量のビームが照射される。この点のことをスポットと呼ぶ。規定量のビームがスポットへ照射されたら、続いて、一時的にビームの照射を停止させた後、次の位置へ照射できるように走査電磁石の電流量が変更される。電流量が変更され、次の照射位置に移動後、再びビームを照射させる。
【００５０】
照射位置を連続的に移動させる方式は、ビームを照射したまま照射位置を変化させる。すなわち、電磁石の励磁量を連続的に変化させながら、照射野内全体を通過するようにビームを照射しながら移動させる。この方法における照射位置ごとの照射量の変化は、走査速度かビームの電流量、あるいはその両方を変調させることで実現する。
【００５１】
離散的な方式による照射の概念図を図３に示す。図３は、立方体の標的４７を照射する例である。粒子ビームは、進行方向におけるある位置で停止し、その停止位置にエネルギーの大部分を付与するため、ビームの停止する深さが標的４７内または標的４７近傍となるようにエネルギーが調整される。図３では、同一エネルギーで照射される面８２付近で停止するエネルギーのビームが選ばれている。この面上に、離散的なビーム照射位置（スポット）がスポット間隔８３で配置されている。一つのスポットで規定量を照射すると、次のスポットへ移動する。スポット８４は、スポット８４を照射するビームの軌跡８５を通るビームで照射される。標的４７内に配置された同一エネルギーのスポットを順次照射し終わると、標的４７内の他の深さ位置を照射するために、ビームを停止させる深さが変更される。
【００５２】
ビームの停止する深さを変化させるためには、患者３０に照射するビームのエネルギーを変化させる。エネルギーを変化させる方法の一つは、粒子ビーム加速装置、すなわち本実施例においてはシンクロトロンの設定を変更することである。粒子はシンクロトロンにおいて設定されたエネルギーになるまで加速されるが、この設定値を変更することで患者３０に入射するエネルギーを変更することができる。この場合、シンクロトロンから取り出されるエネルギーが変化するため、高エネルギービーム輸送系２０を通過する際のエネルギーも変化し、高エネルギービーム輸送系２０の設定変更も必要になる。
【００５３】
図３の例では、同一エネルギーで照射される面８２に相当する領域に主にエネルギーを付与していた。エネルギーを変更することで、例えば図４のような状況となる。図４では、図３で使用したエネルギーよりも低いエネルギーのビームが照射される。そのため、ビームはより浅い位置で停止する。この面を同一エネルギーで照射される面９１で表わす。このエネルギーのビームに対応するスポットの一つであるスポット９２は、スポット９２を照射するビームの軌跡９３を通るビームで照射される。
【００５４】
このように、スポット毎のエネルギー、照射位置、照射量の情報が入力されると、粒子線治療装置１は、標的４７に向けて粒子線を照射し、線量分布を形成する。このスポット毎のエネルギー、照射位置、照射量などの情報を決定するのが治療計画装置５０１である。
【００５５】
図５に示すように照射装置の先端部には、Ｘ線検出器５３が２台配置されている。また、回転装置内には、２台のＸ線発生器５２が配置されている。Ｘ線発生器５２からのＸ線で患者の体内を直交する二つの方向から透視可能である。なお、２方向のＸ線は直行している必要はない。
【００５６】
この２組のＸ線発生器５２とＸ線検出器５３からなるＸ線透視装置５１を用いることで、粒子線照射前にコーンビームＣＴ画像の撮影ができ、また粒子線照射中に標的４７近傍に留置したマーカ６１の位置を計測することができる。
【００５７】
次に、操作者による治療計画立案の流れと治療計画装置５０１の操作の流れを、図６に沿って説明する。図６は治療計画装置５０１の計算内容の流れを示したものである。
【００５８】
治療に先立ち、治療計画用の計画ＣＴ画像が撮像される。治療計画用の計画ＣＴ画像として最も一般的に利用されるのは３次元Ｘ線ＣＴ画像である。３次元Ｘ線ＣＴ画像は、ＣＴ装置などの３次元撮像装置７００により患者の複数の方向から取得した透視画像から、３次元のデータを再構成することで生成する。撮像された計画ＣＴ画像は、データサーバ６００に保存されている。治療計画装置５０１は、この計画ＣＴ画像を利用する。
【００５９】
図６において、治療計画装置５０１は、まず、治療計画の立案が開始されたことを認識（ステップＳ１０１）すると、治療計画装置５０１の操作者である技師（または医師）は、マウス等の入力装置を用いて、データサーバ６００から対象となる計画ＣＴデータを読み込む。治療計画装置５０１は、入力装置の操作により、データサーバ６００からＣＴ画像をメモリ上にコピーする（ステップＳ１０２）。
【００６０】
データサーバ６００からメモリへの計画ＣＴ画像の読み込みが完了し、計画ＣＴ画像が表示装置３１５に表示されると、操作者は表示装置３１５に表示された計画ＣＴ画像を確認しながら、入力装置を用いて、計画ＣＴ画像のスライス、すなわち２次元ＣＴ画像ごとに標的４７として指定すべき領域を入力する（ステップＳ１０３）。ここで入力すべき標的４７は、腫瘍細胞が存在する、あるいは存在する可能性があるために十分な量の放射線を照射すべきと判断された領域である。これを標的４７と呼ぶ。照射線量を極力抑えるべき正常組織が標的４７の近傍に存在するなど、他に評価、制御を必要とする領域がある場合、操作者はそれら正常組織等の領域も同様に指定する。
【００６１】
このステップＳ１０３での領域描画は、過去に描画した領域描画データをもとに治療計画装置５０１が自動で領域指定してもよい。また、領域描画データをもとに治療計画装置５０１が自動で描画したものを、操作者が修正して指定することもできる。他にも、ＭＲＩに代表される異なるモダリティの画像上で実行されてもよい。
【００６２】
次に操作者は、登録された標的４７に対して照射すべきビームの位置やエネルギーの情報を含む治療計画を作成するため、必要な照射パラメータを設定する（ステップＳ１０４）。操作者は、まず、照射方向を設定する。本実施例を適用した粒子線治療装置では、回転装置２１とベッド５０の角度を選択することで、患者の任意の方向からビームの照射を行うことができる。照射方向は一つの標的４７に対して複数設定することが可能である。ある方向からビームが照射される場合、照射時には標的４７の重心位置がアイソセンタ（回転装置２１の回転中心位置）に一致するように位置決めされることが想定される。
【００６３】
他に操作者が決定すべき照射のためのパラメータとしては、ステップＳ１０３で登録した領域に照射すべき線量値（処方線量）と照射回数がある。処方線量は標的４７に照射すべき線量や、正常組織に許容される最大線量が含まれる。
【００６４】
以上のパラメータが決まった後、操作者の指示に従って治療計画装置５０１が自動で計算を行う（ステップＳ１０５）。以下で、治療計画装置５０１が行う線量計算に係わる内容の詳細に関して説明する。
【００６５】
初めに、治療計画装置５０１は、ビーム照射位置を決定する。本実施例で扱う前述した離散的な走査方式であれば、離散的なスポット位置を算出する。照射位置は標的４７を覆うように設定される。照射方向（回転装置２１とベッド５０の角度）として複数の方向が指定されている場合は、各方向に関して同じ操作を行う。
【００６６】
全ての照射位置が決定されると、治療計画装置５０１は照射量の最適化計算を開始する。各スポットへの照射量が、ステップＳ１０４で設定された目標の処方線量に近づくように決定される。この計算では、スポットごとの照射量をパラメータとした目標線量からのずれを数値化した目的関数を用いる方法が広く採用されている。目的関数は線量分布が目標とする線量を満たすほど小さな値となるように定義されており、目標関数が最小となるような照射量を反復計算により探索することで、最適とされる照射量を算出する。
【００６７】
反復計算が終了すると、最終的に各スポットに必要な照射量が定まる。複数あるスポットの照射順序も、この段階で定まる。通常は図３の走査経路８６で示したようにジグザグな経路が設定されるが、走査時間や走査方向を考慮した基準に従って照射順序を並び替えることも可能である。
【００６８】
次に、治療計画装置５０１は演算処理装置により、最終的に得られたスポット位置とスポット照射量を用いて、線量分布を計算する。計算した線量分布の結果は、表示装置３１５に表示される（ステップＳ１０６）。
【００６９】
操作者は、表示装置３１５に表示された線量分布を確認し、線量分布が目標通りの線量分布となっていなければ、ステップＳ１０４に戻り、新たな領域に線量の制約を加えるなどの設定を変更して、目標通りの線量分布が得られるまでステップＳ１０４からステップＳ１０７を繰り返す（ステップＳ１０７）。
【００７０】
ステップＳ１０７において目標通りの線量分布が得られた場合、操作完了となる（ステップＳ１０８）。
【００７１】
次に、図７のフロー図を用いて粒子線照射の全体フローを説明する。図７に粒子線治療装置による粒子線を照射する流れを説明するフロー図を示す。
【００７２】
図７において、最初に、患者３０が治療室に入室する（ステップＳ２０１）。
【００７３】
次いで、患者３０がベッド５０に横になり、本発明の特徴である位置決めを実施する（ステップＳ２０２）。
【００７４】
本ステップＳ２０２では、治療計画装置５０１が決定した情報を中央制御装置３１２が読み込む。位置決めに必要な情報は、中央制御装置３１２から位置決め制御システム３１６へ送信される。位置決めの詳細なフローについては図８を用いて後述する。
【００７５】
次いで、中央制御装置３１２は、粒子線を照射する回転装置の角度を設定する（ステップＳ２０３）。また、照射に必要なスポット毎のエネルギーや照射量などの情報を中央制御装置３１２から照射制御システム３１４に送信する。
【００７６】
次いで、中央制御装置３１２は荷電粒子ビーム発生装置１１と高エネルギービーム輸送系２０と照射野形成装置４０を照射制御システム３１４と共に制御して粒子線を患者へ向けて照射する（ステップＳ２０４）。
【００７７】
本ステップＳ２０４における粒子線の照射中は、標的４７近傍に留置したマーカ６１の位置を連続的に計測する。マーカ６１の計測には、２方向のＸ線透視画像を用いる。２台のＸ線透視装置５１におけるＸ線を発生する位置と、Ｘ線透視画像上に写ったマーカ６１の位置を結ぶ２本の直線が交わる点、又は最も接近する点をマーカ６１の計測位置とする。
【００７８】
本ステップでは、そのマーカ６１の計測位置が照射許可領域内にあるときのみ、粒子線を標的４７に向けて照射する。照射許可領域は、マーカ６１の計画位置に対して例えば±２ｍｍの領域として設定される。マーカ６１の計画位置は、計画ＣＴ画像におけるマーカ６１の位置である。
【００７９】
粒子線の照射が完了すると、ステップＳ１０５において、ガントリーを回転して粒子線を照射する。
【００８０】
次いで、中央制御装置３１２は、計画したすべての照射方向からの照射が完了したか否かを判定し（ステップＳ２０５）、完了したと判定されたときは処理をステップＳ２０６に進めて、患者３０はベッド５０から降り治療室から退室する（ステップＳ２０６）。これに対し、完了していないと判定されたときは処理をステップＳ２０３に戻し、全ての照射方向からの照射を完了させる。
【００８１】
次に、図８のフロー図を用いて図７のステップＳ２０２の位置決め処理の詳細について説明する。図８は本実施例の粒子線治療装置における位置決めの流れを説明するフロー図である。
【００８２】
図８において、患者３０がベッド５０上で横になり、基準ＣＴ画像を撮像したときと同じ姿勢を再現される（ステップＳ３０１）。両腕の位置や膝を曲げる角度などを、サポート具を用いて再現する。標的４７の位置がアイソセンタに一致するようなおおよその位置にベッド５０を移動する。
【００８３】
おおよその位置へのベッド５０の移動後、例えば患者３０が息を吐いた状態で２台のＸ線透視装置５１を用いて２つの方向のＸ線透視画像が取得される（ステップＳ３０２）。取得したＸ線透視画像と、計画ＣＴ画像から計算により得られる投影画像を比較し、患者の位置が計画位置に設置されるようにベッド５０を移動させる。
【００８４】
ステップＳ３０２でのベッド５０の移動後、２台のＸ線透視装置５１によってＸ線撮像を連続的に実行し、標的４７の近傍にあるマーカ６１が計画位置の近傍にあることを確認する（ステップＳ３０３）。例えば、マーカ６１が計画位置を中心とした±２ｍｍの領域にあるときに粒子線を照射することから、その領域を照射許可領域と呼ぶ。ここでは、マーカ６１が照射許可領域内を通過する時間が十分に長いことを確認する。
【００８５】
次いで、２台のＸ線透視装置５１を用いてコーンビームＣＴ画像を撮像する（ステップＳ３０４）。コーンビームＣＴ画像を取得するため、二つの方向のＸ線透視画像を取得しながらガントリーを回転させる。ガントリーは少なくとも９０度分回転させる必要があり、回転角度が多いほど高画質のコーンビームＣＴ画像を取得することができる。
【００８６】
Ｘ線透視画像は、二つの方向から取得したものをペアとして記録する。認識装置３１６Ａは、二つの画像のそれぞれの中で、マーカ６１が写っている位置をテンプレートマッチングなどにより認識する。Ｘ線発生器５２からＸ線が発生する位置とマーカ６１が写った画像上の位置を結ぶ２本の直線が交わる、又は最も接近する位置をマーカ６１の位置とする。以上により、ペアで記録されたＸ線透視画像毎にマーカ６１の位置を記録することができる。認識装置３１６Ａはこのマーカ６１の位置が照射許可領域内にあるＸ線透視画像のみを抽出し、算出装置３１６Ｂは抽出したＸ線透視画像からコーンビームＣＴ画像を再構成する。
【００８７】
次いで、算出装置３１６Ｂによって、計画ＣＴ上で指定されている標的４７や臓器の情報を取得したコーンビームＣＴ画像上に転記する（ステップＳ３０５）。この転記は、操作者が画面上でコーンビームＣＴ画像上の標的４７と臓器を指定することもできるが、非剛体レジストレーションと呼ばれる手法を用いることもできる。非剛体レジストレーションを用いた手法では、計画ＣＴからコーンビームＣＴへの歪みベクトルを算出し、その歪みベクトルに合わせて計画ＣＴ上の標的４７と臓器を歪めることで、コーンビームＣＴ上の標的４７と臓器が指定される。
【００８８】
次いで、算出装置３１６Ｂによって、計画ＣＴ上とコーンビームＣＴ上のそれぞれで、マーカ６１の座標と、標的４７の重心座標を算出し、マーカ６１の座標を起点とした標的重心座標へのベクトルを算出する（ステップＳ３０６）。
【００８９】
次いで、算出装置３１６Ｂによって、算出されたベクトルのｘ，ｙ，ｚのそれぞれの成分を、計画ＣＴとコーンビームＣＴのそれぞれに分けて操作画面上に表示する（ステップＳ３０７）。
【００９０】
このステップＳ３０７では、ステップＳ３０６において算出したベクトルを表示装置３１５の表示画面３１５Ａに表示する。この表示画面３１５Ａの一例を図９に示す。図９はベクトルの表示画面の一例を示す図である。
【００９１】
図９に示す表示画面３１５Ａにおいて、算出された標的４７の重心座標が現在値と記載されている領域３１５Ａ１に表示され、Ｘ，Ｙ，Ｚのそれぞれが表示される。また、その下の領域３１５Ａ２には、計画ＣＴ上における計画時の標的４７の重心座標の計画値が表示される。３段目の領域３１５Ａ３には、軸毎の差分が表示される。
【００９２】
操作者は、図９に示すような表示画面３１５Ａの領域３１５Ａ１，３１５Ａ２，３１５Ａ３に表示された値が許容範囲内であり、問題ないことを確認した後には、位置決めが完了したと判断し、粒子線の照射開始を中央制御装置３１２に対して指示する。この指示の入力を受けて、図７のステップＳ２０４に示すように、中央制御装置３１２は粒子線の照射を開始する。
【００９３】
これに対し、より計画時の位置に一致させたいと判断した時は、再びステップＳ３０２でのベッド５０の移動を実施し、その後の位置決め処理を再度実行する。
【００９４】
なお、ステップＳ３０６において算出するのは重心座標である必要はなく、たとえば標的４７の３軸それぞれの中点などとすることができる。
【００９５】
また、ステップＳ３０６において用いる標的４７の情報は点である必要はない。
【００９６】
また、ステップＳ３０７では、コーンビームＣＴ画像または計画ＣＴ画像上にマーカ６１の位置を一致させたうえで、標的４７の形状全体を重ねて描き、差分を強調表示することができる。以下、図１０および図１１を用いてその表示画面３１５Ｂについて説明する。図１０に計画ＣＴとコーンビーム画像での標的４７を重ねて表示する例を示す。図１１に差分を強調表示する際の画面の一例を示す。
【００９７】
そのためには、ステップＳ３０６において、算出装置３１６Ｂは、治療計画ＣＴ画像と照射直前ＣＴ画像におけるマーカ６１の３次元位置を一致させて、治療計画ＣＴ画像と照射直前ＣＴ画像のそれぞれに対応する標的４７の形状の３次元座標位置を演算する。
【００９８】
図１０に示すように、算出装置３１６Ｂによる演算処理の結果、表示装置３１５では、コーンビームＣＴ画像上に治療直前または治療時のマーカ６１と標的輪郭６２が表示される。また、計画ＣＴ画像上の標的輪郭６０が、マーカ６１の位置がコーンビームＣＴ画像のマーカ６１の位置に一致するように、マーカ６１の位置と標的輪郭の相対関係を保って転記した状態で表示される。
【００９９】
さらに、表示画面３１５Ｂでは、コーンビームＣＴでの標的輪郭６２と計画ＣＴでの標的輪郭６０の重なり合わない差分６４が、図１０の中で塗潰された部分として強調表示される。
【０１００】
また、標的４７以外の正常臓器の形状についても重ねて表示させて、その差分を表示することで、正常臓器が計画時と同じ位置にあることを確認することができる。
【０１０１】
更に、計画ＣＴ画像とコーンビームＣＴ画像の値の差分について予め複数の閾値を設けておき、求めた差分が閾値を超える場合には、図１１に示すように、操作画面上に文字と色の違いで差異が大きいことを表示することができる。例えば、表示画面３１５Ｃでは、差分が閾値を下回っている部分については差分６４Ａとして表示し、閾値を超えている部分については差分６４Ｂとして表示するとともに、警告表示６５を表示する。文字と色は閾値毎に異なるものを用いることが好ましい。
【０１０２】
また、ステップＳ３０７では、マーカ６１の位置を基準として一致させた場合に、照射方向から見た標的４７や正常臓器の手前側または奥側の境界までの水等価厚について、コーンビームＣＴ画像と計画ＣＴ画像の差分を表示することも可能である。粒子線では特に水等価厚の変化量が重要である。
【０１０３】
図１２に、水等価厚の差について表示する例を示す。図１２は紙面手前方向から標的４７に向けて粒子線を照射するときの例である。図１２は水等価厚を表示する表示画面の一例を示す図である。
【０１０４】
そのためには、ステップＳ３０６において、算出装置３１６Ｂは、コーンビームＣＴでの標的輪郭と計画ＣＴでの標的輪郭の差を、粒子線を照射する方向の水等価に換算しておく。その上で、ステップＳ３０７において、換算された水等価のうち、差が大きい場合には濃い色で塗潰して領域３１５Ｄ１のように表示し、差が小さい場合には薄い色で塗潰して領域３１５Ｄ２のように表示画面３１５Ｄに表示する。また、差がない場合には塗潰さない。これを、照射方向に対して手前側（浅い側）と奥側（深い側）の両方に対して計算し、それぞれを切り替えて表示する。
【０１０５】
また、ステップＳ３０６において標的４７の位置が大きく変化したと演算される場合には、算出装置３１６Ｂは、例えば標的４７が一致するようにベッド５０の移動量を算出し、その値に基づいてベッド５０を移動させることもできる。ベッド５０を移動させた場合に、照射許可領域とマーカ６１の位置が離れる場合には、照射許可領域の位置についてもベッド５０の移動量に合わせて移動させることが可能である。
【０１０６】
また、ステップＳ３０４では、コーンビームＣＴ画像を取得する際に得られるマーカ６１の軌跡を用いて、認識装置３１６Ａは、マーカ６１の軌跡を複数の領域毎に分割し、算出装置３１６Ｂは、その分割領域毎にコーンビームＣＴ画像を再構成することができる。
【０１０７】
また、算出装置３１６Ｂは、複数の領域ごとに分割した画像を用いてステップＳ３０７においてマーカ６１基準の標的４７の位置を計算することで、分割領域毎のマーカ６１と標的４７の関係を算出することが可能である。これによれば、算出した結果を基に、マーカ６１と標的４７の関係が計画時と最も近い分割領域の位置を表示することが可能である。また、その分割領域の中心が、照射許可領域の中心と重なるようにベッド５０の移動量を算出し、ベッドを移動させることができる。
【０１０８】
次に、本実施例の効果について説明する。
【０１０９】
上述した本実施例の粒子線治療装置１は、治療計画の作成時に撮像した治療計画ＣＴ画像の中のマーカ６１の位置と粒子線の照射直前ＣＴ画像の中のマーカ６１の位置との位置関係、および治療計画ＣＴ画像の中の標的４７の位置と照射直前ＣＴ画像の中の標的４７の位置との位置関係を算出する算出装置３１６Ｂと、算出装置３１６Ｂで算出された治療計画ＣＴ画像と照射直前ＣＴ画像における標的４７およびマーカ６１の位置関係を表示する表示装置３１５と、を備えている。
【０１１０】
これによって、マーカ６１の位置と標的４７の位置が計画時と一致しているか否かを判断するために有意義な情報を操作者が容易に把握することができる形態で提供することができる。このため、標的４７の位置を計画時と同等とするために要する時間を従来に比べて短縮することができる。特に動く標的４７を照射する際の準備時間を短縮することができる。従って、標的４７に対して粒子線を照射する際に、従来に比べてより精度良く短時間で粒子線を標的４７に照射することが可能となる。
【０１１１】
また、算出装置３１６Ｂは、マーカ６１を基準として、標的４７の重心の位置関係を算出するため、計画時と照射直前の標的４７の位置が一致しているか否かをより正確に求めることができ、標的４７の位置をより正確に計画時と一致させることが可能となる。
【０１１２】
更に、位置関係の表示では、治療計画ＣＴ画像と照射直前ＣＴ画像におけるマーカ６１の位置を一致させて、治療計画ＣＴ画像と照射直前ＣＴ画像のそれぞれに対応する標的４７の形状を重ねて表示することで、標的４７のうち、計画時と一致しない領域がある場合には容易にそのズレを確認することができ、標的４７の位置を計画時と同等とするために要する時間をより短縮することができる。
【０１１３】
また、位置関係の表示では、標的４７の形状の差分を強調表示することにより、治療計画時と治療直前での標的４７の位置関係のズレをより明確に把握することができる。
【０１１４】
更に、位置関係の表示では、標的４７の形状の差分を予め定めた閾値と比較し、差分が閾値を超えると判定されたときは、閾値を超えた差分を強調表示することで、治療計画時と治療直前での標的４７の位置関係のズレを更に明確に把握することができる。
【０１１５】
また、強調表示では、画面上に文字を表示するとともに、色を変えて表示することによっても、治療計画時と治療直前での標的４７の位置関係のズレを更に明確に把握することができる。
【０１１６】
更に、算出装置３１６Ｂは、治療計画ＣＴ画像と照射直前ＣＴ画像におけるマーカ６１の位置を一致させて、治療計画ＣＴ画像と照射直前ＣＴ画像のそれぞれにおける標的４７や正常臓器の照射方向の手前側および／または奥側の境界までの水等価厚とその差分を求め、表示装置３１５では、治療計画ＣＴ画像と照射直前ＣＴ画像の水等価厚の差分を表示することで、放射線として粒子線を照射する際に重要な指標となる水等価厚の治療計画時に対する照射直前のズレについても容易に把握することができる。すなわち、マーカ６１の位置と標的４７の位置が計画時と一致しているか否かを水等価厚のズレに応じても判断することができる。例えば、標的４７の位置が計画時と少々ずれが生じている場合であっても、水等価厚が同じであると判断される場合は標的４７の位置は計画時と同等であると判断できる。このため、標的４７の位置が計画時と一致しているか否かの判断材料をより多く提供することができ、標的４７を照射する際の準備時間をより短縮することができる。
【０１１７】
また、算出装置３１６Ｂは、マーカ６１と標的４７との位置関係に基づいて、標的４７を有する患者３０を載せるベッド５０の移動量を算出することにより、標的４７の位置を計画時と同等の位置に容易に移動させることができ、標的４７を照射する際の準備時間をより短縮することができる。
【０１１８】
更に、マーカ６１および標的４７をＸ線により複数の方向から撮像するＸ線透視装置５１を更に備え、算出装置３１６Ｂは、Ｘ線透視装置５１によって取得したＸ線画像から画像再構成を行い、照射直前ＣＴ画像としてコーンビームＣＴ画像を生成することで、いわゆる動体追跡で利用される装置を用いてマーカ６１の位置と標的４７の位置が計画時と一致しているか否かを判断することができ、新たな装置を追加することなく標的４７を照射する際の準備時間を短縮することが可能となる。
【０１１９】
また、Ｘ線透視装置５１によって取得したＸ線画像から患者３０内の標的４７の３次元位置を認識する認識装置３１６Ａを更に備え、算出装置３１６Ｂは、Ｘ線透視装置５１によって取得したＸ線画像の中から、認識装置３１６Ａによって認識した標的４７の位置が動体追跡治療の治療粒子線照射条件を満たす画像から画像再構成を行うことにより、粒子線を照射するタイミングでの体内構造を反映したコーンビームＣＴ画像を取得することができる。これにより、治療粒子線を照射する状態の患者３０により近い状態を３次元ＣＴ画像化（可視化）することができ、放射線治療の患者位置決めにおいて必要とされる、治療放射線照射状態での標的４７等の位置関係を把握可能な画像を取得することができる。その結果、治療計画通りの精密な患者位置決めが可能となる。
【０１２０】
更に、患者３０内の標的４７の位置が動体追跡治療の照射許可領域を満たすときに標的４７に粒子線を照射するよう制御する照射制御システム３１４を更に備えたことで、標的４７に対する治療粒子線の照射精度をより向上させることができる。
【０１２１】
また、照射制御システム３１４は、マーカ６１と標的４７との位置関係に基づいて、照射許可領域を変更することにより、照射許可領域を照射直前の標的４７の位置に応じて設定することが可能となり、標的４７に対する治療粒子線の照射精度をより向上させることができる。
【０１２２】
＜その他＞  
なお、本発明は上記の実施例に限られず、種々の変形、応用が可能なものである。上述した実施例は本発明を分かりやすく説明するために詳細に説明したものであり、必ずしも説明した全ての構成を備えるものに限定されない。
【０１２３】
例えば、上述の実施例では、最初にステップＳ３０２でＸ線透視画像を用いてベッドを移動した後にコーンビームＣＴ画像を撮像する例を示したが、Ｘ線透視画像の位置合わせをすることなくコーンビームＣＴ画像を取得することもできる。
【０１２４】
この場合は、コーンビームＣＴ画像取得時に得られるマーカ６１の軌跡から、照射許可領域となるべき位置を選択し、選択した位置を基準のマーカ６１の位置として、計画時のものとマーカ６１と標的４７の位置関係を比較する。比較した結果、標的４７の位置が一致していれば、選択した位置が照射許可領域と一致するようなベッド５０の移動量を算出し、算出した結果に基づいてベッド５０を移動させることができる。
【０１２５】
また、治療直前のＸ線ＣＴ画像を取得する手段としてコーンビームＣＴ画像を用いる手法について述べたが、Ｘ線ＣＴ装置を治療室内に別に置き、その装置にて撮像したデータをコーンビームＣＴ画像の代わりに用いることができる。
【０１２６】
更に、照射許可領域の大きさと、コーンビームＣＴ画像を再構成するためのマーカ６１の領域の大きさとは一致している必要はない。コーンビームＣＴ画像を再構成するためのマーカ６１の領域は、マーカ６１の動きの大きさと必要なコーンビームＣＴ画像の画質から適宜調整することが有効である。
【符号の説明】
【０１２７】
１…粒子線治療装置
３１２…中央制御装置
３１４…照射制御システム（治療粒子線制御装置）
３１５…表示装置
３１５Ａ，３１５Ｂ，３１５Ｃ，３１５Ｄ…表示画面
３１６…位置決め制御システム
３１６Ａ…認識装置
３１６Ｂ…算出装置
４０…照射野形成装置
４６…患者
４７…標的
５１…Ｘ線透視装置
５２…Ｘ線発生器
５３…Ｘ線検出器
６１…マーカ
６２…標的輪郭
６４，６４Ａ，６４Ｂ…差分
６５…警告表示
５０１…治療計画装置
７００…３次元撮像装置
【図１】
2019072393.tif 000003
【図２】
2019072393.tif 000004
【図３】
2019072393.tif 000005
【図４】
2019072393.tif 000006
【図５】
2019072393.tif 000007
【図６】
2019072393.tif 000008
【図７】
2019072393.tif 000009
【図８】
2019072393.tif 000010
【図９】
2019072393.tif 000011
【図１０】
2019072393.tif 000012
【図１１】
2019072393.tif 000013
【図１２】
2019072393.tif 000014