(57)【要約】本発明の１つの態様では、医用画像データに基づく幾何学的モデルの確立方法が提供され、それは、  医用画像データを読み取るステップと、  医用画像データと組織の種類との間の変換関係によって組織の種類を定義するステップと、  組織グループ分けの数を決定するステップと、  医用画像データと濃度との変換関係によって組織の密度を定義するステップと、  組織及び密度に関する情報を携帯する３Ｄコードマトリックスを確立するステップと、  幾何学的モデルを生成するステップと、を含む。  該幾何学的モデルの確立方法は、医用画像データと組織の種類との間の変換関係によって、実際のニーズに応じて組織グループ分けの数を決定することができ、それにより組織の種類、要素の組成及び密度をより正確に提供し、確立された幾何学的モデルは、医用画像データが反映された実際の状況によりマッチングする。
2019502505.tif 000002
【特許請求の範囲】
【請求項１】
医用画像データを読み取るステップと、
前記医用画像データと組織の種類との間の変換関係によって前記組織の種類を定義するステップと、
組織グループ分けの数を決定するステップと、
前記医用画像データと濃度との変換関係によって前記組織の密度を定義するステップと、
前記組織及び前記密度に関する情報を携帯する３Ｄコードマトリックスを確立するステップと、
幾何学的モデルを生成するステップと、を含む、
医用画像データに基づく幾何学的モデルの確立方法。
【請求項２】
前記幾何学的モデルの確立方法は、中性子捕捉療法に適用され、前記幾何学的モデルの確立方法は、さらに、
Ｂ−１０濃度を与えるステップと、
Ｂ−１０濃度に関する情報を携帯する３Ｄコードマトリックスを確立するステップと、を含むことを特徴とする、
請求項１に記載の医用画像データに基づく幾何学的モデルの確立方法。
【請求項３】
前記組織グループ分けの数は、ユーザから手動で定義された組織グループ分けの数に、データベース内の既存の４種の組織グループ分けの数又は１４種の組織グループ分けの数を加えた数であることを特徴とする、
請求項１に記載の医用画像データに基づく幾何学的モデルの確立方法。
【請求項４】
前記幾何学的モデルの確立方法は、さらに、
３Ｄ組織コードマトリックスを確立するステップと、
３Ｄ密度コードマトリックスを確立するステップと、を含むことを特徴とする、
請求項３に記載の医用画像データに基づく幾何学的モデルの確立方法。
【請求項５】
前記幾何学的モデルは、ＭＣＮＰソフトウェア入力ファイルに必要な格子カード、セルカード、サーフェスカード及びマテリアルカードを含むことを特徴とする、
請求項１に記載の医用画像データに基づく幾何学的モデルの確立方法。
【請求項６】
医用画像データを読み取るステップと、
ＲＯＩ境界を定義又は読み取るステップと、
医用画像ボクセルが前記ＲＯＩ境界内にあるかどうかを判断するステップと、
前記医用画像ボクセルが前記ＲＯＩ境界内にある場合、各前記ＲＯＩ境界内のボクセルに１つの特定の組織及び密度を与えることによって、ユーザが組織の種類及び密度を手動で定義するステップに進み、又は前記医用画像データと組織の種類／密度との間の変換関係によってＲＯＩ組織の種類及び密度を自動的に定義するステップに進み、
前記医用画像ボクセルが前記ＲＯＩ境界内にない場合、前記医用画像データと組織の種類との間の変換関係によって、組織の種類を自動的に定義するステップに進み、かつ前記医用画像データと組織の種類との間の変換関係によって、組織の密度を定義するステップに進み、
前記組織及び前記密度に関する情報を携帯する３Ｄコードマトリックスを確立するステップと、
幾何学的モデルを生成するステップと、を含む、
医用画像データに基づく幾何学的モデルの確立方法。
【請求項７】
前記幾何学的モデルの確立方法は、中性子捕捉療法に適用され、前記幾何学的モデルの確立方法は、さらに、
Ｂ−１０濃度を与えるステップと、
Ｂ−１０濃度に関する情報を携帯する３Ｄコードマトリックスを確立するステップと、を含むことを特徴とする、
請求項６に記載の医用画像データに基づく幾何学的モデルの確立方法。
【請求項８】
前記組織グループ分けの数は、ユーザから手動で定義された組織グループ分けの数に、データベース内の既存の４種の組織グループ分けの数又は１４種の組織グループ分けの数を加えた数であることを特徴とする、
請求項６に記載の医用画像データに基づく幾何学的モデルの確立方法。
【請求項９】
前記幾何学的モデルの確立方法は、さらに、
３Ｄ組織コードマトリックスを確立するステップと、
３Ｄ密度コードマトリックスを確立するステップと、を含むことを特徴とする、
請求項８に記載の医用画像データに基づく幾何学的モデルの確立方法。
【請求項１０】
前記幾何学的モデルは、ＭＣＮＰソフトウェア入力ファイルに必要な格子カード、セルカード、サーフェスカード及びマテリアルカードを含むことを特徴とする、
請求項６に記載の医用画像データに基づく幾何学的モデルの確立方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は幾何学的モデルの設定方法に関し、さらに詳しくは、医用画像データに基づく幾何学的モデルの設定方法に関する。
【背景技術】
【０００２】
原子科学の発展に従って、コバルト60、線形加速器、電子ビームなどの放射線療法は、すでにがん治療の主な手段の一つとなった。しかし、従来の光子または電子療法は、放射線そのものの物理的条件の制限で腫瘍細胞を殺すとともに、ビーム経路上の数多くの正常組織に損傷を与える。また、腫瘍細胞により放射線に対する感受性の度合いが異なっており、従来の放射線療法では、放射線耐性の高い悪性腫瘍（例、多形神経膠芽腫（glioblastoma multiforme）、黒色腫（melanoma））に対する治療効果が良くない。
腫瘍の周囲の正常組織への放射線損傷を軽減するすために、化学療法（chemotherapy）における標的療法が、放射線療法に用いられている。また、放射線耐性の高い腫瘍細胞に対し、現在では生物学的効果比（relative biological effectiveness, RBE）の高い放射線源が積極的に開発されている（例えば、陽子線治療、重粒子治療、中性子捕捉療法など）。このうち、中性子捕捉療法は、上記の2つの構想を結びつけたものである。例えば、ホウ素中性子捕捉療法では、ホウ素含有薬物が腫瘍細胞に特異的に集まり、高精度な中性子ビームの制御と合わせることで、従来の放射線と比べて、より良いがん治療オプションを提供する。
ホウ素中性子捕捉療法（Boron Neutron Capture Therapy, BNCT）はホウ素（10B）含有薬物が熱中性子に対し大きい捕獲断面積を持つ特性を利用し、10B(n,α) 7Li中性子捕捉と核分裂反応により4Heと7Liという2種の重荷電粒子を生成する。図1と図2は、それぞれホウ素中性子捕捉の反応概略図と10B(n,α) 7Li中性子捕捉の原子核反応式を示す。2種の重荷電粒子は平均エネルギーが2.33MeVであり、高い線エネルギー付与(Linear Energy Transfer, LET)及び短い射程という特徴を持つ。α粒子の線エネルギー付与と射程はそれぞれ150 keV/μm、8μmであり、7Li重荷粒子の場合、それぞれ175 keV/μm、5μmである。2種の粒子の合計射程が細胞のサイズに近いので、生体への放射線損害を細胞レベルに抑えられる。ホウ素含有薬物を選択的に腫瘍細胞に集め、適切な中性子源と合わせることで、正常組織に大きな損害を与えないで腫瘍細胞を部分的に殺せる。  ３次元モデルは、科学実験分析、科学実験シミュレーション分野で広く使用されている。例えば、核放射線及び保護の分野で、特定の放射線条件下での人体の吸収線量をシミュレートするために、常にコンピュータ技術を用いて医用画像データに対して様々な処理を行ってＭＣＮＰに必要な正確な格子モデルを確立し、かつＭＣＮＰ（モンテカルロプログラム）と組み合わせてシミュレーション計算を行う。
【０００３】
モンテカルロ法は、現在で放射照射内部の３次元空間核粒子の衝突軌道及びエネルギー分布を正確にシミュレートできるツールであり、モンテカルロ法と複雑な３次元人体解剖モデルとの組み合わせは、コンピュータ技術におけるシミュレーションの飛躍を表している。診断放射線検査において、正確な人体臓器線量評価は、放射線療法に非常に有益である。現在、世界中で様々な人体モデルが成功に確立され、かつモンテカルロシミュレーションプログラムと組み合わせられ、放射線環境下での人体の吸収線量を正確に計算及び評価する。３次元人体解剖モデルがモンテカルロプログラムに正常に変換することに必要な幾何学的説明は、モンテカルロシミュレーション計算を行う前提条件であり、現時点では世界中でモンテカルロシミュレーション研究のホットスポット及び難点である。
【０００４】
磁気共鳴イメージング（Ｍａｇｎｅｔｉｃ  Ｒｅｓｏｎａｎｃｅ  Ｉｍａｇｉｎｇ、ＭＲＩ）又はコンピュータ断層撮影（Ｃｏｍｐｕｔｅｄ  Ｔｏｍｏｇｒａｐｈｙ、ＣＴ）などの医用画像データは、人体体内の特徴に対して、より詳細な組織幾何学的構造情報を提供し、人体内部構造の物理的モデル化にデータ基礎を提供することができる。中性子捕捉療法の分野で、医用画像データに基づいてＭＣＮＰに必要な幾何学的モデルをどのように確立するかが重要な問題であり、つまり、医用画像データに基づいてＭＣＮＰソフトウェアの入力に必要な格子モデルをどのように確立するか、治療計画の精度を向上させることである。
したがって、治療計画の精度を向上させる医用画像データに基づいてＭＣＮＰに必要な幾何学的モデルを確立する方法を提供する必要がある。
【発明の概要】
【０００５】
本発明の１つの態様では、医用画像データに基づく幾何学的モデルの確立方法が提供され、それは、
医用画像データを読み取るステップと、
医用画像データと組織の種類との間の変換関係によって組織の種類を定義するステップと、
組織グループ分けの数を決定するステップと、
医用画像データと濃度との変換関係によって組織の密度を定義するステップと、
組織及び密度に関する情報を携帯する３Ｄコードマトリックスを確立するステップと、
幾何学的モデルを生成するステップと、を含む。
【０００６】
該幾何学的モデルの確立方法は、医用画像データと組織の種類との間の変換関係によって、実際のニーズに応じて組織グループ分けの数を決定することができ、それにより組織の種類、要素の組成及び密度をより正確に提供し、確立された幾何学的モデルは、医用画像データが反映された実際の状況によりマッチングする。
【０００７】
好ましく、該幾何学的モデルの確立方法は、中性子捕捉療法に適用され、それは、さらに、Ｂ−１０濃度を与えるステップと、Ｂ−１０濃度に関する情報を携帯する３Ｄコードマトリックスを確立するステップと、を含む。Ｂ−１０濃度に関する情報で標識された幾何学的モデルは、各組織におけるホウ素含有薬物の濃度を明らかに了解することができ、続いて中性子捕捉療法の照射シミュレーションを行う時に、実際の状況をより現実的に反映する。
【０００８】
組織グループ分けの数は、ユーザから手動で定義された組織グループ分けの数にデータベース内の既存の４種の組織グループ分けの数又は１４種の組織グループ分けの数を加えた数である。既存のデータベースに対応する組織グループ分けの数が確立されない場合、ユーザから１つの新しい組織グループ分けの数を定義することができる。このように、既存のデータベースに対応する組織グループ分けの数に完全にマッチングしない場合、近似選択のみが行われることを回避し、それによりモデリングの精度を効果的に向上させる。
【０００９】
さらに好ましく、該幾何学的モデルの確立方法は、さらに、３Ｄ組織コードマトリックスを確立するステップと、３Ｄ密度コードマトリックスを確立するステップと、を含む。医用画像データのスライスに基づいて、対応する変換関係によって、各スライスに対応する組織コード及び密度コードを確立し、それにより３Ｄ組織コードマトリックス及び３Ｄ密度コードマトリックスを確立する。
【００１０】
幾何学的モデルは、ＭＣＮＰソフトウェア入力ファイルに必要な格子カード、セルカード、サーフェスカード及びマテリアルカードを含む。医用画像データにより、最終的にＭＣＮＰソフトウェア入力ファイルに必要な格子カード、セルカード、サーフェスカード及びマテリアルカードを得て、それによりシミュレーション計算に理論的根拠を提供し、かつ正確なシミュレーション結果を得る。
【００１１】
本発明の別の態様では、医用画像データに基づく幾何学的モデルの確立方法が提供され、それは、
医用画像データを読み取るステップと、
ＲＯＩ境界を定義又は読み取るステップと、
医用画像ボクセルがＲＯＩ境界内にあるかどうかを判断するステップと、
そうであれば、各ＲＯＩ境界内のボクセルに１つの特定の組織及び密度を与えることによって、ユーザが組織の種類及び密度を手動で定義するステップに進み、又は医用画像データと組織の種類／密度との間の変換関係によってＲＯＩ組織の種類及び密度を自動的に定義するステップに進み、
そうでなければ、医用画像データと組織の種類との間の変換関係によって、組織の種類を自動的に定義するステップに進み、かつ医用画像データと組織の種類との間の変換関係によって、組織密度を定義するステップに進み、
組織及び密度に関する情報を携帯する３Ｄコードマトリックスを確立するステップと、
幾何学的モデルを生成するステップと、を含む。
【００１２】
いわゆるＲＯＩは、関心領域（以下、ＲＯＩと総称する）であり、ユーザはＲＯＩ組織の種類、元素組成及び密度を手動で定義することができる。医用画像ボクセルがＲＯＩ境界内にない場合、医用画像データと組織の種類との間の変換関係に基づいて組織の種類の定義を行い、かつ実際のニーズに応じて組織グループ分け数を決定し、それにより組織の種類、元素組成及び密度をより正確に提供し、確立された幾何学的モデルは、医用画像データによって反映された実際の状況とより密接にマッチングする。
【００１３】
好ましく、該幾何学的モデルの確立方法は、中性子捕捉療法に適用され、幾何学的モデルの確立方法は、Ｂ−１０濃度を与えるステップと、Ｂ−１０濃度に関する情報を携帯する３Ｄコードマトリックスを確立するステップと、を含む。Ｂ−１０濃度に関する情報で標識された幾何学的モデルは、各組織におけるホウ素含有薬物の濃度を明らかに了解することができ、続いて中性子捕捉療法の照射シミュレーションを行う時に、実際の状況をより現実的に反映する。
【００１４】
組織グループ分けの数は、ユーザから手動で定義された組織グループ分けの数にデータベース内の既存の４種の組織グループ分けの数又は１４種の組織グループ分けの数を加えた数である。既存のデータベースに対応する組織グループ分けの数が確立されない場合、ユーザから１つの新しい組織グループ分けの数を定義することができる。このように、既存のデータベースに対応する組織グループ分けの数に完全にマッチングしない場合、近似選択のみが行われることを回避し、それによりモデリングの精度を効果的に向上させる。
【００１５】
さらに好ましく、該幾何学的モデルの確立方法は、さらに、３Ｄ組織コードマトリックスを確立するステップと、３Ｄ密度コードマトリックスを確立するステップと、を含む。医用画像データのスライスに基づいて、対応する変換関係によって、各スライスに対応する組織コード及び密度コードを確立し、それにより３Ｄ組織コードマトリックス及び３Ｄ密度コードマトリックスを確立する。
【００１６】
幾何学的モデルは、ＭＣＮＰソフトウェア入力ファイルに必要な格子カード、セルカード、サーフェスカード及びマテリアルカードを含む。医用画像データにより、最終的にＭＣＮＰソフトウェア入力ファイルに必要な格子カード、セルカード、サーフェスカード及びマテリアルカードを得て、それによりシミュレーション計算に理論的根拠を提供し、かつ正確なシミュレーション結果を得る。
【００１７】
医用画像データは、磁気共鳴イメージング（Ｍａｇｎｅｔｉｃ  Ｒｅｓｏｎａｎｃｅ  Ｉｍａｇｉｎｇ、ＭＲＩ）又はコンピュータ断層撮影（Ｃｏｍｐｕｔｅｄ  Ｔｏｍｏｇｒａｐｈｙ、ＣＴ）であってもよく、以下の実施例において、コンピュータ断層撮影（ＣＴ）のデータに基づいて説明され、ＣＴのファイル形式は一般的にＤＩＣＯＭである。当業者には周知のように、さらに、別の医用画像データを使用することができ、該医用画像データが組織種類、密度及びＢ−１０濃度に関する情報を携帯するＭＣＮＰ格子モデルに変換できれば、本発明の開示する医用画像データに基づく幾何学的モデルの確立方法に適用することができる。
【００１８】
本発明の実施例における有益な効果及び／又は特徴は以下のとおりである。
１、関心領域（ＲＯＩ）の組織種類、元素組成及び密度を手動で定義する。
２、非ＲＯＩのＣＴ画像ボクセルポイントについて、組織種類マッチングを自動的に行うことができ、ＣＴ値の違いにより、既存のデータベースは４種類又は１４種類の要素からなる組織を区別することができ、また実際の実験結果により他の数の異なる元素からなる組織を決定することができる。
【００１９】
３、ＣＴ値被覆範囲の広いＲＯＩについて、例えば粘膜チャンバーであり、唯一の組織密度を与えることは困難であり、ユーザは本発明の実施例が開示する方法に従ってＣＴ値と組織／密度の変換を自動的に行うことができる。
【００２０】
４、正常血液中のホウ素含有薬物の濃度、組織／腫瘍−血液ホウ素濃度比などのパラメータを入力した後、本発明の実施例が開示する方法は、Ｂ−１０元素を全てのボクセルポイントに自動的に書き込む。
【００２１】
５、最後に得られた３次元ＭＣＮＰ格子モデルは、組織種類（元素組成）、密度及びＢ−１０濃度などの情報を携帯する。
【図面の簡単な説明】
【００２２】
【図１】図１はホウ素中性子捕捉反応の概略図である。
【図２】図２は１０Ｂ（ｎ、α）７Ｌｉ中性子捕捉核反応式である。
【図３】図３は本発明の実施例における医用画像データに基づく幾何学的モデル確立方法の論理ブロック図である。
【図４】図４はＣＴ値（ＨＵ）と組織密度回帰曲線式及び相関係数図である。
【発明を実施するための形態】
【００２３】
以下は、図面を参照して本発明の実施例をさらに詳細に説明し、当業者であれば、明細書の文字を参照して実施することができる。
好ましくは、中性子捕捉療法に用いられる医用画像データに基づく幾何学的モデル確立方法を本発明の実施例とする。以下は、中性子捕捉療法、特にホウ素中性子捕捉療法を簡単に説明する。
【００２４】
中性子捕捉療法は効果的ながん治療の手段として、近年ではその適用が増加しており、そのうち、ホウ素中性子捕捉療法が最も一般的なものとなった。ホウ素中性子捕捉療法に用いられる中性子は原子炉または加速器で供給できる。本発明の実施形態は加速器ホウ素中性子捕捉療法（Accelerated-based Boron Neutron Capture Therapy）を例とする。加速器ホウ素中性子捕捉療法の基本モジュールは、一般的に荷電粒子（陽子、デューテリウム原子核など）の加速に用いられる加速器、ターゲット、熱除去システム及びビーム整形アセンブリを含む。加速後の荷電粒子と金属ターゲットとの作用により中性子が生成され、必要な中性子収率及びエネルギー、提供可能な加速荷電粒子のエネルギー及び電流、及び、金属ターゲットの物理的・化学的特性などにより、適切な原子核反応が選定される。よく検討されている原子核反応は7Li(p,n)7Be及び9Be(p,n)9Bであり、この両方はすべて吸熱反応でエネルギー閾値がそれぞれ1.881MeVと2.055MeVである。ホウ素中性子捕捉療法の理想的中性子源はkeVエネルギーレベルの熱外中性子なので、理論的には、エネルギーが閾値よりやや高い陽子によるリチウムターゲットへの衝撃で、比較的低いエネルギーの中性子が生成され、あまり多くの減速処理を要しないで臨床適用が可能になる。しかし、リチウム（Li）及びベリリウム（Be）の2種のターゲットは、閾値エネルギーの陽子と作用する断面が大きくないので、十分な中性子束を確保するために、一般的には比較的高いエネルギーを持つ陽子で原子核反応を引き起こされる。
ホウ素中性子捕捉療法（Boron Neutron Capture Therapy, BNCT）はホウ素（10B）含有薬物が熱中性子に対し大きい捕獲断面積を持つ特性を利用し、10B(n,α) 7Li中性子捕捉と核分裂反応により4Heと7Liという2種の重荷電粒子を生成する。図1と図2は、それぞれホウ素中性子捕捉の反応概略図と10B(n,α) 7Li中性子捕捉の原子核反応式を示す。2種の重荷電粒子は平均エネルギーが2.33MeVであり、高い線エネルギー付与(Linear Energy Transfer, LET)及び短い射程という特徴を持つ。α粒子の線エネルギー付与と射程はそれぞれ150 keV/μm、8μmであり、7Li重荷粒子の場合、それぞれ175 keV/μm、5μmである。2種の粒子の合計射程が細胞のサイズに近いので、生体への放射線損害を細胞レベルに抑えられる。ホウ素含有薬物を選択的に腫瘍細胞に集め、適切な中性子源と合わせることで、正常組織に大きな損害を与えないで腫瘍細胞を部分的に殺せる。  本発明の実施例の目的は、医用画像データを組織の種類、密度及びＢ−１０濃度に関する情報を携帯するＭＣＮＰ格子モデルに変換し、ホウ素中性子捕捉療法の組織線量シミュレーション計算を行うことである。医用画像データは、磁気共鳴イメージング（Ｍａｇｎｅｔｉｃ  Ｒｅｓｏｎａｎｃｅ  Ｉｍａｇｉｎｇ、ＭＲＩ）又はコンピュータ断層撮影（Ｃｏｍｐｕｔｅｄ  Ｔｏｍｏｇｒａｐｈｙ、ＣＴ）であってもよく、以下の実施例において、コンピュータ断層撮影（ＣＴ）のデータに基づいて説明され、ＣＴのファイル形式は一般的にＤＩＣＯＭである。当業者には周知のように、さらに、別の医用画像データを使用することができ、該医用画像データが組織種類、密度及びＢ−１０濃度に関する情報を携帯するＭＣＮＰ格子モデルに変換できれば、本発明の開示する医用画像データに基づく幾何学的モデルの確立方法に適用することができる。
【００２５】
簡単に言えば、本発明の実施例が開示する医用画像データに基づく幾何学的モデル確立方法は、主に以下のフローを含む。
１、コンピュータ断層撮影（ＤＩＣＯＭ形式）を入力し、ＣＴ画像が対応するインタフェース上に表示される。
【００２６】
２、ＤＩＣＯＭファイルで定義されたＲＯＩの境界を自動的に読み取り、また別のＲＯＩを追加することができる。
３、ソフトウェアは、ＣＴ画像上の各ボクセルポイント（ｖｏｘｅｌ）の位置を決定し、すべてのボクセルポイントを「ＲＯＩ境界内にある」及び「ＲＯＩ境界外にある」と分類する。
【００２７】
４、ＲＯＩ境界内のボクセルポイントについて、実際のニーズに応じて、ユーザにより各ＲＯＩに相対的唯一の組織の種類及び組織密度を手動で定義し、又はＣＴ値と組織の種類、密度を自動的にマッチングさせ、所定ＣＴ値被覆範囲の大きいＲＯＩ（粘膜チャンバーなど）唯一の組織の種類（元素組成）及び密度が、線量計算の誤差をもたらすことを回避する。
【００２８】
５、ＲＯＩ境界外のボクセルポイントについて、組織種類の定義を自動的に行い、ＣＴ値の違いにより、４種類又は１４種類の異なる元素からなる組織を区別することができ、ユーザは実際の判断に従って、ＩＣＲＵ−４６レポートに示された４種類の組織（以下に詳述する）を選択して使用し、又はＶａｎｄｅｒｓｔｒａｅｔｅｎらによって２００７年に公表された１４種類の異なる元素からなる組織を使用することができる（以下に詳述する）。
【００２９】
６、手動で定義されない密度のボクセルポイントについて、ＣＴ値の違いにより、自動的に密度を与え、合計で９６個の密度グループ分けを区別することができる。
７、ユーザは、正常血中のホウ素含有薬剤の濃度、腫瘍-血液中ホウ素濃度比、組織-血液中ホウ素濃度比などのパラメータを手動で入力し、全てのボクセルポイントにＢ−１０元素濃度を書き込む。
【００３０】
８、ソフトウェアは、組織の種類（元素組成）、組織密度、組織Ｂ−１０濃度及び他の情報を統合し、３次元ＭＣＮＰ格子モデルを生成し、かつＭＣＮＰ入力ファイルの指定形式の格子カード（ｌａｔｔｉｃｅ  ｃａｒｄ）、セルカード（ｃｅｌｌ  ｃａｒｄ）、サーフェスカード（ｓｕｒｆａｃｅ  ｃａｒｄ）、及びマテリアルカード（ｍａｔｅｒｉａｌ  ｃａｒｄ）を書き込む。
【００３１】
具体的に、図３に示すとおり、本発明の開示する医用画像データに基づく幾何学的モデル確立方法は、２つの実施例を含む。
本発明の第1実施例にでは、医用画像データに基づく幾何学的モデルの確立方法が提供され、それは、
医用画像データを読み取るステップと、
医用画像データと組織の種類との間の変換関係によって組織の種類を定義するステップと、
組織グループ分けの数を決定するステップと、
医用画像データと濃度との変換関係によって組織の密度を定義するステップと、
組織及び密度に関する情報を携帯する３Ｄコードマトリックスを確立するステップと、
幾何学的モデルを生成するステップと、を含む。
【００３２】
該幾何学的モデルの確立方法は、医用画像データと組織の種類との間の変換関係によって、実際のニーズに応じて組織グループ分けの数を決定することができ、それにより組織の種類、要素の組成及び密度をより正確に提供し、確立された幾何学的モデルは、医用画像データが反映された実際の状況によりマッチングする。
【００３３】
好ましく、該幾何学的モデルの確立方法は、中性子捕捉療法に適用され、それは、さらに、Ｂ−１０濃度を与えるステップと、Ｂ−１０濃度に関する情報を携帯する３Ｄコードマトリックスを確立するステップと、を含む。Ｂ−１０濃度に関する情報で標識された幾何学的モデルは、各組織におけるホウ素含有薬物の濃度を明らかに了解することができ、続いて中性子捕捉療法の照射シミュレーションを行う時に、実際の状況をより現実的に反映する。
【００３４】
組織グループ分けの数は、ユーザから手動で定義された組織グループ分けの数にデータベース内の既存の４種の組織グループ分けの数又は１４種の組織グループ分けの数を加えた数である。既存のデータベースに対応する組織グループ分けの数が確立されない場合、ユーザから１つの新しい組織グループ分けの数を定義することができる。このように、既存のデータベースに対応する組織グループ分けの数に完全にマッチングしない場合、近似選択のみが行われることを回避し、それによりモデリングの精度を効果的に向上させる。
【００３５】
さらに好ましく、該幾何学的モデルの確立方法は、さらに、３Ｄ組織コードマトリックスを確立するステップと、３Ｄ密度コードマトリックスを確立するステップと、を含む。医用画像データのスライスに基づいて、対応する変換関係によって、各スライスに対応する組織コード及び密度コードを確立し、それにより３Ｄ組織コードマトリックス及び３Ｄ密度コードマトリックスを確立する。
【００３６】
幾何学的モデルは、ＭＣＮＰソフトウェア入力ファイルに必要な格子カード、セルカード、サーフェスカード及びマテリアルカードを含む。医用画像データにより、最終的にＭＣＮＰソフトウェア入力ファイルに必要な格子カード、セルカード、サーフェスカード及びマテリアルカードを得て、それによりシミュレーション計算に理論的根拠を提供し、かつ正確なシミュレーション結果を得る。
本発明の第2実施例にでは、医用画像データに基づく幾何学的モデルの確立方法が提供され、それは、
医用画像データを読み取るステップと、
ＲＯＩ境界を定義又は読み取るステップと、
医用画像ボクセルがＲＯＩ境界内にあるかどうかを判断するステップと、
そうであれば、各ＲＯＩ境界内のボクセルに１つの特定の組織及び密度を与えることによって、ユーザが組織の種類及び密度を手動で定義するステップに進み、又は医用画像データと組織の種類／密度との間の変換関係によってＲＯＩ組織の種類及び密度を自動的に定義するステップに進み、
そうでなければ、医用画像データと組織の種類との間の変換関係によって、組織の種類を自動的に定義するステップに進み、かつ医用画像データと組織の種類との間の変換関係によって、組織密度を定義するステップに進み、
組織及び密度に関する情報を携帯する３Ｄコードマトリックスを確立するステップと、
幾何学的モデルを生成するステップと、を含む。
【００３７】
いわゆるＲＯＩは、関心領域（以下、ＲＯＩと総称する）であり、ユーザはＲＯＩ組織の種類、元素組成及び密度を手動で定義することができる。医用画像ボクセルがＲＯＩ境界内にない場合、医用画像データと組織の種類との間の変換関係に基づいて組織の種類の定義を行い、かつ実際のニーズに応じて組織グループ分け数を決定し、それにより組織の種類、元素組成及び密度をより正確に提供し、確立された幾何学的モデルは、医用画像データによって反映された実際の状況とより密接にマッチングする。
【００３８】
好ましく、該幾何学的モデルの確立方法は、中性子捕捉療法に適用され、幾何学的モデルの確立方法は、Ｂ−１０濃度を与えるステップと、Ｂ−１０濃度に関する情報を携帯する３Ｄコードマトリックスを確立するステップと、を含む。Ｂ−１０濃度に関する情報で標識された幾何学的モデルは、各組織におけるホウ素含有薬物の濃度を明らかに了解することができ、続いて中性子捕捉療法の照射シミュレーションを行う時に、実際の状況をより現実的に反映する。
【００３９】
組織グループ分けの数は、ユーザから手動で定義された組織グループ分けの数にデータベース内の既存の４種の組織グループ分けの数又は１４種の組織グループ分けの数を加えた数である。既存のデータベースに対応する組織グループ分けの数が確立されない場合、ユーザから１つの新しい組織グループ分けの数を定義することができる。このように、既存のデータベースに対応する組織グループ分けの数に完全にマッチングしない場合、近似選択のみが行われることを回避し、それによりモデリングの精度を効果的に向上させる。
【００４０】
さらに好ましく、該幾何学的モデルの確立方法は、さらに、３Ｄ組織コードマトリックスを確立するステップと、３Ｄ密度コードマトリックスを確立するステップと、を含む。医用画像データのスライスに基づいて、対応する変換関係によって、各スライスに対応する組織コード及び密度コードを確立し、それにより３Ｄ組織コードマトリックス及び３Ｄ密度コードマトリックスを確立する。
【００４１】
幾何学的モデルは、ＭＣＮＰソフトウェア入力ファイルに必要な格子カード、セルカード、サーフェスカード及びマテリアルカードを含む。医用画像データにより、最終的にＭＣＮＰソフトウェア入力ファイルに必要な格子カード、セルカード、サーフェスカード及びマテリアルカードを得て、それによりシミュレーション計算に理論的根拠を提供し、かつ正確なシミュレーション結果を得る。
【００４２】
医用画像データは、磁気共鳴イメージング（Ｍａｇｎｅｔｉｃ  Ｒｅｓｏｎａｎｃｅ  Ｉｍａｇｉｎｇ、ＭＲＩ）又はコンピュータ断層撮影（Ｃｏｍｐｕｔｅｄ  Ｔｏｍｏｇｒａｐｈｙ、ＣＴ）であってもよく、以下の実施例において、コンピュータ断層撮影（ＣＴ）のデータに基づいて説明され、ＣＴのファイル形式は一般的にＤＩＣＯＭである。当業者には周知のように、さらに、別の医用画像データを使用することができ、該医用画像データが組織種類、密度及びＢ−１０濃度に関する情報を携帯するＭＣＮＰ格子モデルに変換できれば、本発明の開示する医用画像データに基づく幾何学的モデルの確立方法に適用することができる。
【００４３】
図４に示すとおり、以下は、本発明の開示する医用画像データに基づく幾何学的モデル確立方法における既存のデータベース中のＣＴ値−組織種類及びＣＴ値−組織密度の変換関係図を説明する。
【００４４】
ＣＴ値は、Ｈｏｕｎｓｆｉｅｌｄ  Ｕｎｉｔ（ＨＵ）とも呼ばれ、光減衰係数を反映する単位であり、その定義は式１に示すとおり、
【数１】
2019502505.tif 000003
【００４５】
Ｖａｎｄｅｒｓｔｒａｅｔｅｎらによって２００７年に公表された文献（Ｂａｒｂａｒａ  Ｖａｎｄｅｒｓｔｒａｅｔｅｎら、「モンテカルロ線量計算のためのＣＴ値の組織パラメーターへのコンバージェンス：マルチセンター研究」、Ｐｈｙｓ．Ｍｅｄ．Ｂｉｏｌ．５２（２００７）５３９−５６２）を参照し、ＣＴ値の違いにより、１種類の空気（ａｉｒ）、１種類の肺組織（ｌｕｎｇ）、２種類の軟部組織（脂肪（ａｄｉｐｏｓｅ）及び筋肉軟部組織（ｍｕｓｃｌｅ））、１０種類の骨組織（ｂｏｎｅ）を変換することができ、換言すれば、異なるＣＴ値は、１４種類の異なる要素からなる組織に対応し、以下の表１に示すとおりである。
【００４６】
【表１】
2019502505.tif 000004
ＩＣＲＵ−４６号レポート（Ｉｎｔｅｒｎａｔｉｏｎａｌ  Ｃｏｍｍｉｓｓｉｏｎ  ｏｎ  Ｒａｄｉａｔｉｏｎ  Ｕｎｉｔｓ  ａｎｄ  Ｍｅａｓｕｒｅｍｅｎｔｓ、Ｐｈｏｔｏｎ、ｅｌｅｃｔｒｏｎ、ｐｒｏｔｏｎ  ａｎｄ  ｎｅｕｔｒｏｎ  ｉｎｔｅｒａｃｔｉｏｎ  ｄａｔａ  ｆｏｒ  ｂｏｄｙ  ｔｉｓｓｕｅｓ、ＩＣＲＵ−４６、Ｔｅｃｈ．Ｒｅｐ．１９９２．）を参照し、４種類のヒト脳部の主要な組織を取り、空気（ａｉｒ）、脳組織（ａｄｕｌｔ  ｂｒａｉｎ）、皮膚（ａｄｕｌｔ  ｓｋｉｎ）、頭蓋骨（ｃｒａｎｉｕｍ）を含み、その対応する密度及び元素組成は表２に示すとおりである。
【００４７】
【表２】
2019502505.tif 000005
同様にＶａｎｄｅｒｓｔｒａｅｔｅｎらの文献を参照し、該文献は病院の実際の実験値を統合し、ＣＴ値に対応する組織密度の関係式が整理され、図４に示すとおり、本発明の開示する医用画像データに基づく幾何学的モデルの確立方法は、図４の３つの回帰式を使用し、ＣＴ値（−１０００〜２０００）を９６組の密度グループ分けに分割し、表３に示すとおりである。
【００４８】
【表３】
2019502505.tif 000006
本発明の実施例における有益な効果及び／又は特徴は以下のとおりである。
【００４９】
１、関心領域（ＲＯＩ）の組織種類、元素組成及び密度を手動で定義する。
２、非ＲＯＩのＣＴ画像ボクセルポイントについて、組織種類マッチングを自動的に行うことができ、ＣＴ値の違いにより、既存のデータベースは４種類又は１４種類の要素からなる組織を区別することができ、また実際の実験結果により他の数の異なる元素からなる組織を決定することができる。
【００５０】
３、ＣＴ値被覆範囲の広いＲＯＩについて、例えば粘膜チャンバーであり、唯一の組織密度を与えることは困難であり、ユーザは本発明の実施例が開示する方法に従ってＣＴ値と組織／密度の変換を自動的に行うことができる。
【００５１】
４、正常血液中のホウ素含有薬物の濃度、組織／腫瘍−血液ホウ素濃度比などのパラメータを入力した後、本発明の実施例が開示する方法は、Ｂ−１０元素を全てのボクセルポイントに自動的に書き込む。
５、最後に得られた３次元ＭＣＮＰ格子モデルは、組織種類（元素組成）、密度及びＢ−１０濃度などの情報を携帯する。
【００５２】
以上で本発明の説明的な具体的な実施形態を説明し、それにより当業者が本発明を理解しやすく、しかし、本発明は具体的な実施形態の範囲に限定するものではなく、当業者にとっては、様々な変化が特許請求に限定及び確定された本発明の範囲内にあれば、これらの変化が明らかであり、いずれも本発明の保護範囲内にある。
【図１】
2019502505.tif 000007
【図２】
2019502505.tif 000008
【図３】
2019502505.tif 000009
【図４】
2019502505.tif 000010
【国際調査報告】
2019502505.tif 000011
2019502505.tif 000012
2019502505.tif 000013
2019502505.tif 000014