特許 7589363 放射線源から得られたα線スペクトルを分析する方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2024-11-15

(45)【発行日】2024-11-25

(54)【発明の名称】放射線源から得られたα線スペクトルを分析する方法

(51)【国際特許分類】

   G01T 1/167 20060101AFI20241118BHJP

   G01T 1/36 20060101ALI20241118BHJP

   A61N 5/10 20060101ALN20241118BHJP

【ＦＩ】

G01T1/167 C

G01T1/36 B

A61N5/10

【請求項の数】 9

(21)【出願番号】P 2023546119

(86)(22)【出願日】2021-12-09

(65)【公表番号】

(43)【公表日】2024-02-01

(86)【国際出願番号】 KR2021018656

(87)【国際公開番号】W WO2022164011

(87)【国際公開日】2022-08-04

【審査請求日】2023-07-28

(31)【優先権主張番号】10-2021-0013059

(32)【優先日】2021-01-29

(33)【優先権主張国・地域又は機関】KR

(73)【特許権者】

【識別番号】518460417

【氏名又は名称】コリア・インスティテュート・オブ・ラディオロジカル・アンド・メディカル・サイエンシーズ

(74)【代理人】

【識別番号】100108453

【弁理士】

【氏名又は名称】村山 靖彦

(74)【代理人】

【識別番号】100110364

【弁理士】

【氏名又は名称】実広 信哉

(74)【代理人】

【識別番号】100133400

【弁理士】

【氏名又は名称】阿部 達彦

(72)【発明者】

【氏名】グン・ア・キム

(72)【発明者】

【氏名】ジョン・グク・キム

(72)【発明者】

【氏名】イル・ハン・リム

(72)【発明者】

【氏名】ビョン・ス・キム

(72)【発明者】

【氏名】カン・ヒョン・ソン

【審査官】中尾 太郎

(56)【参考文献】

【文献】特開平０５－２０９９６５（ＪＰ，Ａ）

【文献】特開平０６－２６５６３９（ＪＰ，Ａ）

【文献】特開平０８－１３６６５８（ＪＰ，Ａ）

【文献】特開２００７－０９３４７１（ＪＰ，Ａ）

【文献】特開２０１９－２１９３４０（ＪＰ，Ａ）

【文献】特表２０２０－５２８１４９（ＪＰ，Ａ）

【文献】米国特許出願公開第２００７／００６９１４７（ＵＳ，Ａ１）

【文献】Teemu SIISKONEN et al.，New Approach to Alpha Spectrum Analysis: Iterative Monte Carlo Simulations and Fitting，Progress in Nuclear Science and Technology，日本，Atomic Energy Society of Japan，2011年10月，Vol. 2，pp. 437-441

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ｇ０１Ｔ １／１６７

Ｇ０１Ｔ １／３６

Ａ６１Ｎ ５／１０

(57)【特許請求の範囲】

【請求項1】

放射線源から得られたα線スペクトルを分析する方法であって、
前記放射線源は、互いに異なる第１核種および第２核種を含み、
前記第１核種および前記第２核種の放射能関数を利用したシミュレーションによって、前記第１核種および前記第２核種から放出されるα線に対する模倣スペクトルを得るステップ；
前記放射線源から放出されるα線を測定して、前記α線スペクトルを得るステップ；および
前記模倣スペクトルを基底（ｂａｓｉｓ）として、前記α線スペクトルを前記第１核種から放出された第１α線スペクトルと前記第２核種から放出された第２α線スペクトルとに分解（ｄｅｃｏｍｐｏｓｉｔｉｏｎ）するステップ；
を含み、
前記測定は、
体の一部を模倣した模倣体を対象として行われる、分析方法。

【請求項2】

前記第２核種は、前記第１核種の娘核種である、請求項１に記載の分析方法。

【請求項3】

前記第１核種および前記第２核種は、［^２２５Ａｃ，^２１３Ｂｉ］、［^２１１Ａｔ，^２１２Ｐｏ］、［^２１２Ｂｉ，^２１２Ｐｏ］、［^２１３Ｂｉ，^２１３Ｐｏ］、［^２１２Ｐｂ，^２１２Ｐｏ］、［^２２３Ｒａ，^２１２Ｂｉ］、［^２２７Ｔｈ，^２１１Ｂｉ］の組み合わせから選択されるいずれか一つである、請求項１に記載の分析方法。

【請求項4】

前記放射線源は、第１放射線源および第２放射線源を含み、
前記第１放射線源および前記第２放射線源は、前記第１核種および前記第２核種の分布が相異する、請求項２に記載の分析方法。

【請求項5】

前記測定で、
前記第１放射線源および前記第２放射線源が離隔して配置されている、請求項４に記載の分析方法。

【請求項6】

分離以後に、
特定測定位置での前記第１放射線源からのα線と前記第２放射線源からのα線との比率を導出するステップをさらに含む、請求項５に記載の分析方法。

【請求項7】

前記放射能関数は、前記第１核種および前記第２核種の半減期情報および放出確率情報を利用して得られる、請求項１に記載の分析方法。

【請求項8】

測定スペクトルの分離は、
特定測定位置での前記第１α線スペクトルと前記第２α線スペクトルとの係数比を得るステップを含む、請求項１に記載の分析方法。

【請求項9】

放射線源から得られたα線スペクトルを分析する方法において、
前記放射線源は、
親核種を含む第１放射線源；および
前記第１放射線源から離隔されており、前記親核種の娘核種を含む第２放射線源を含み、
前記親核種および前記娘核種の放射能関数を利用したシミュレーションによって、前記親核種および前記娘核種から放出されるα線に対する模倣スペクトルを得るステップ；
前記放射線源から放出されるα線を測定して、前記α線スペクトルを得るステップ；
前記模倣スペクトルを基底として、前記α線スペクトルを前記親核種から放出された第１α線スペクトルと前記娘核種から放出された第２α線スペクトルとに分解するステップ；および
特定測定位置で測定されたα線に対して、前記親核種からのα線と前記娘核種からのα線に区分するステップ；を含み、
前記測定は、
体の一部を模倣した模倣体を対象として行われる、分析方法。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、放射線源から得られたα線スペクトルを分析する方法に関する。

【背景技術】

【0002】

α線放出核種を利用した放射線治療（Ｔａｒｇｅｔｅｄ ａｌｐｈａ ｔｈｅｒａｐｙ、ＴＡＴ）が最近注目されている。

【0003】

ＴＡＴにおいて、当該核種（親核種）が崩壊して娘核種を生成する場合、治療をしようとする位置に娘核種も存在するようになる。

【0004】

計画された量のα線放出娘核種は、治療効果を極大化することができるが、過剰な娘核種は、毒性（ｔｏｘｉｃｉｔｙ）を増加させる原因になる。

【0005】

したがって、当該娘核種が治療をしようとする位置の親核種によって生成された娘核種であるのか、または他の位置の親核種によって生成された娘核種が生理作用によって当該位置に流入したものであるのかに対する分析が必要である。

【発明の概要】

【発明が解決しようとする課題】

【0006】

したがって、本発明の目的は、放射線源から得られたα線スペクトルを分析する方法を提供することである。

【課題を解決するための手段】

【0007】

上記本発明の目的は、放射線源から得られたα線スペクトルを分析する方法において、前記放射線源は、互いに異なる第１核種および第２核種を含み、前記第１核種および前記第２核種の放射能関数を利用したシミュレーションによって、前記第１核種および前記第２核種から放出されるα線に対する模倣スペクトルを得るステップ；前記放射線源から放出されるα線を測定して、前記α線スペクトルを得るステップ；および前記模倣スペクトルを基底（ｂａｓｉｓ）として、前記α線スペクトルを前記第１核種から放出された第１α線スペクトルと前記第２核種から放出された第２α線スペクトルとに分解（ｄｅｃｏｍｐｏｓｉｔｉｏｎ）するステップ；を含むことによって達成される。

【0008】

前記第２核種は、前記第１核種の娘核種であってもよい。

【0009】

前記第１核種と前記第２核種は、［^２２５Ａｃ，^２１３Ｂｉ］、［^２１１Ａｔ，^２１２Ｐｏ］、［^２１２Ｂｉ，^２１２Ｐｏ］、［^２１３Ｂｉ，^２１３Ｐｏ］、［^２１２Ｐｂ，^２１２Ｐｏ］、［^２２３Ｒａ，^２１２Ｂｉ］、［^２２７Ｔｈ，^２１１Ｂｉ］の組み合わせから選択されるいずれか一つであってもよい。

【0010】

前記放射線源は、第１放射線源および第２放射線源を含み、前記第１放射線源および前記第２放射線源は、前記第１核種と前記第２核種の分布が相異してもよい。

【0011】

前記測定で、前記第１放射線源および前記第２放射線源が離隔して配置されてもよい。

【0012】

前記分離以後に、特定測定位置での前記第１放射線源からのα線と前記第２放射線源からのα線との比率を導出するステップをさらに含んでもよい。

【0013】

前記測定は、体の一部を模倣した模倣体を対象として行われてもよい。

【0014】

前記放射能関数は、前記第１核種および前記第２核種の半減期情報および放出確率情報を利用して得られてもよい。

【0015】

前記測定スペクトルの分離は、特定測定位置での前記第１α線スペクトルと前記第２α線スペクトルとの係数比を得るステップを含んでもよい。

【0016】

上記本発明の目的は、放射線源から得られたα線スペクトルを分析する方法において、前記放射線源は、親核種を含む第１放射線源；および前記第１放射線源から離隔されており、前記親核種の娘核種を含む第２放射線源を含み、前記親核種および前記娘核種の放射能関数を利用したシミュレーションによって、前記親核種および前記娘核種から放出されるα線に対する模倣スペクトルを得るステップ；前記放射線源から放出されるα線を測定して、前記α線スペクトルを得るステップ；前記模倣スペクトルを基底として、前記α線スペクトルを前記親核種から放出された第１α線スペクトルと前記娘核種から放出された第２α線スペクトルとに分解するステップ；および特定測定位置で測定されたα線に対して、前記親核種からのα線と前記娘核種からのα線に区分するステップ；を含むことによって達成される。

【0017】

前記模倣スペクトルは、前記測定を行う前に、ライブラリ形態で保存されてもよく、前記測定の時には、既に完成された模倣スペクトルライブラリを利用して、前記スペクトル分解ステップを達成することができる。

【0018】

前記測定は、体の一部を模倣した模倣体を対象として行われてもよい。

【発明の効果】

【0019】

本発明に係ると、放射線源から得られたα線スペクトルを分析する方法が提供される。

【図面の簡単な説明】

【0020】

【図1】本発明の一実施形態によるα線スペクトルの分析方法に対する順序図である。

【図2】本発明の実験例で得た模倣スペクトルである。

【図3】本発明の実験例で用いた模倣体と放射線源を示した図である。

【図4】本発明の実験例で用いた模倣体と放射線源を示した図である。

【図5】本発明の実験例で得られた２次元α線分布画像である。

【図6】本発明の実験例で領域２でのα線スペクトルを示した図である。

【図7a】本発明の実験例で各領域での第１α線スペクトルと第２α線スペクトルを示した図である。

【図7b】本発明の実験例で各領域での第１α線スペクトルと第２α線スペクトルを示した図である。

【図7c】本発明の実験例で各領域での第１α線スペクトルと第２α線スペクトルを示した図である。

【図7d】本発明の実験例で各領域での第１α線スペクトルと第２α線スペクトルを示した図である。

【図8a】本発明の実験例で^２２５Ａｃからのα線画像を示した図である。

【図8b】本発明の実験例で^２１３Ｂｉからのα線画像を示した図である。

【図9】本発明の実験例で^２２５Ａｃからのα線と^２１３Ｂｉからのα線画像を合算して示した図である。

【図10】本発明の実験例で実際に注入された核種の放射能と推定された核種の放射能との比較を示した図である。

【発明を実施するための形態】

【0021】

以下、図面を参照して本発明をより詳しく説明する。

【0022】

添付図面は、本発明の技術的思想をより具体的に説明するために示した一例に過ぎないので、本発明の思想が添付図面に限定されるものではない。また、添付図面は各構成要素の間の関係を説明するために、大きさと間隔などが実際と異なり、誇張されていてもよい。

【0023】

以下の説明において、α線放出核種を利用した放射線治療での親核種と娘核種の区分を例示して説明するが、本発明はこれに限定されない。本発明は、原子力発電所の安全に必要なα線核種監視またはα線放出核種ミルキングシステム（親核種から娘核種を分離して、放射性同位元素を生成するシステム）での品質管理などに用いられてもよい。

【0024】

図１を参照して、本発明に係るα線スペクトルの分析方法について説明する。

【0025】

まず、分析しようとする放射線源を決める（Ｓ１００）。

【0026】

放射線源は、第１核種および第２核種を含む。

【0027】

第１核種および第２核種は、親核種と娘核種の関係であってもよく、これに限定されないが、［^２２５Ａｃ，^２１３Ｂｉ］、［^２１１Ａｔ，^２１２Ｐｏ］、［^２１２Ｂｉ，^２１２Ｐｏ］、［^２１３Ｂｉ，^２１３Ｐｏ］、［^２１２Ｐｂ，^２１２Ｐｏ］、［^２２３Ｒａ，^２１２Ｂｉ］、［^２２７Ｔｈ，^２１１Ｂｉ］の組み合わせから選択されてもよい。

【0028】

次に、模倣スペクトルを得る（Ｓ２００）。

【0029】

模倣スペクトルは、第１核種と第２核種に対する放射能関数を利用して得る。放射能関数は、各核種に対する半減期情報と放出情報から得られ、これに限定されないが、下記の数式１のベイトマン方程式（Ｂａｔｅｍａｎ ｅｑｕａｔｉｏｎ）に情報を入力して得られる。第１核種に対する放射能関数は、第１核種とその娘核種の情報が反映されており、第２核種に対する放射能関数は、第２核種とその娘核種の情報が反映されている。

【0030】

【数1】

【0031】

数式１のｔは、特定時間であり、λ_ｉおよびＮ_ｉは、それぞれ選定した核種の崩壊曲線にある核種のうち、ｉ番目核種の崩壊定数と量であり、Ｎ_ｎ（ｔ）は、特定時間ｔで存在するｎ番目核種の量である。

【0032】

得られた放射能関数、第１核種の初期放射能情報（例えば、１０ｋＢｑ）と第２核種の初期放射能情報（例えば、１０ｋＢｑ）を利用して、シミュレーションによって、第１核種と第２核種から放出されるα線に関する模倣スペクトル（模倣エネルギースペクトル）を得る。

【0033】

シミュレーションでは、これに限定されないが、ＧＡＴＥ（Ｇｅａｎｔ４ ａｐｐｌｉｃａｔｉｏｎ ｆｏｒ ｔｏｍｏｇｒａｐｈｉｃ ｅｍｉｓｓｉｏｎ）コードを用いたモンテカルロシミュレーションを利用してもよい。他の実施形態では、他のコードとして、ＦＬＵＫＡコード、ＭＣＮＰＸコードまたはＧｅａｎｔ４コードを用いてもよい。

【0034】

模倣では、第１核種および第２核種のエネルギー別係数情報を得るようになり、得られたエネルギー別係数情報から模倣スペクトルを得る。

【0035】

次に、放射線源から放出されるα線を測定して、α線スペクトルを得る（Ｓ３００）。

【0036】

測定では、試験片に放射線源を注入し、α線画像化装置で撮影して、２次元ピクセルごとにα線スペクトルを得ることができる。

【0037】

α線画像化装置は、シンチレータと光センサなどを組み合わせて設けられてもよい。

【0038】

試験片は、小動物の薄片や体の一部を模倣した模倣体であってもよい。模倣体は、例えば、腎臓などのような臓器を模倣したものであってもよいが、これに限定されない。

【0039】

測定で、放射線源は互いに離隔されている複数で設けられ、各放射線源の核種分布は相異してもよい。例えば、第１放射線源は第１核種のみ含み、第２放射線源は第２核種のみ含むか、または第１放射線源は第１核種および第２核種を含み、第２放射線源は第２核種のみ含んでもよい。または、第１放射線源および第２放射線源のそれぞれは、第１核種および第２核種を両方含むが、第１核種と第２核種との比率が相異してもよい。

【0040】

その後、α線スペクトルを分離する（Ｓ４００）。

【0041】

α線スペクトルは、分離によって第１核種から放出されたα線に対する第１α線スペクトルと第２核種から放出されたα線に対する第２α線スペクトルとに分離される。

【0042】

分離過程を詳しく説明すれば、次の通りである。

【0043】

分離は、模倣スペクトルを基礎ないし基底として行われ、具体的には、基底分解（ｂａｓｉｓ ｄｅｃｏｍｐｏｓｉｔｉｏｎ）方式で行われてもよい。

【0044】

ｘ、ｙ位置でのα線スペクトル（ｂ_ｘ，ｙ）を模倣スペクトルで得た第１核種と第２核種の模倣スペクトル（Ａ）との係数比に分離する。分離のためには、これに限定されないが、以下の数式２の公式を用いてもよい。

【0045】

【数2】

【0046】

第１核種として、１０ｋＢｑの^２２５Ａｃを用い、第２核種として、１０ｋＢｑの^２１３Ｂｉを用いた場合、撮影されたα線スペクトル（ｂ_ｘ，ｙ）を上記公式に入れた結果は、ｗ^＊ _ｘ，ｙである。すなわち、ｗ^＊ _ｘ，ｙは、ｂ_ｘ，ｙをＡで表現した時の係数である。例えば、Ａが［３，１］として撮影されたｂ_ｘ，ｙが１９である場合、上記公式は、ｂ_ｘ，ｙ＝ｗ１＊３＋ｗ２＊１を満たし、特定制約条件があるｗ^＊ _ｘ，ｙを捜し、その答えがｗ^＊ _ｘ，ｙ＝［５，４］とすると、特定位置のｘ，ｙでの^２２５Ａｃは、５＊１０ｋＢｑである５０ｋＢｑであり、ｘ，ｙでの^２１３Ｂｉは、４＊１０ｋＢｑである４０ｋＢｑになる。

【0047】

すなわち、分離ステップは、特定測定位置での前記第１α線スペクトルと前記第２α線スペクトルとの係数比を得るステップを含む。

【0048】

最後に、α線の核種比率を導出する（Ｓ５００）。

【0049】

このステップでは、特定位置でのα線が第１核種から放出されたのか、または第２核種から放出されたのか、およびその比率を導出する。すなわちα線が放出された線源を区分するのである。

【0050】

前で得たｗ^＊ _ｘ，ｙは、２つの値で、１つは、第１α線スペクトルに対する係数比、２番目の値は、第２α線スペクトルに対する係数比である。これによって、特定ｘ，ｙ位置での第１核種と第２核種との存在比を知ることができる。

【0051】

第１核種と第２核種の存在比の把握を二次元（２Ｄ）領域全体で行って、核種別放射能分布画像をそれぞれ獲得することができる。これによって追加的に第１核種と第２核種の色を異なるように示した後に融合して、フュージョン（Ｆｕｓｉｏｎ）画像を獲得し、これを利用して核種別α線分布を可視化することができる。

【0052】

他の実施形態では、別途の画像化作業なしに位置別核種情報を提供することもできる。

【0053】

以上で説明した本発明は、エネルギー分解能が相対的に低い設備を利用した測定でも放出されたα線の放射線源を区分することができるようにする。

【0054】

本発明に係る方法を利用して、画像化システムで核種別分布画像を獲得可能である。

【0055】

また、ＴＡＴのための非臨床実験で、当該分布画像を利用して、娘核種の生成機構を追跡し、過剰な娘核種による毒性防止、および治療効果を極大化するのに寄与することができる。

【0056】

以下、実験例によって本発明を詳しく説明する。

【0057】

以下の実験例では、分離しようとするα線スペクトルを実際の測定ではないシミュレーションによって得た。

【0058】

実験例での核種は、親核種である^２２５Ａｃと娘核種である^２１３Ｂｉである。

【0059】

親核種である^２２５Ａｃの崩壊曲線は次の通りである。

【0060】

親核種に対する模倣スペクトルでは、^２２５Ａｃを含んで、^２２１Ｆｒ、^２１３Ｂｉなど^２２５Ａｃの全ての娘核種が放出するエネルギースペクトルを示し、娘核種に対する模倣スペクトルでは、^２１３Ｂｉを含んで、^２１３Ｐｏ、^２０９Ｔｌなど^２１３Ｂｉの全ての娘核種が放出するエネルギースペクトルのみを含む。

【0061】

【化1】

【0062】

親核種である^２２５Ａｃと娘核種である^２１３Ｂｉの両方に対して、１０ｋＢｑの放射能での模倣スペクトルを図２のように得た。図２は、数式２のＡに対応する。

【0063】

模倣スペクトルは、ベイトマン方程式（Ｂａｔｅｍａｎ ｅｑｕａｔｉｏｎ）によって獲得した時間別核種別放射能情報およびα線画像化装置をＧＡＴＥ（Ｇｅａｎｔ４ ａｐｐｌｉｃａｔｉｏｎ ｆｏｒ ｔｏｍｏｇｒａｐｈｉｃ ｅｍｉｓｓｉｏｎ）コードに適用したモンテカルロシミュレーションによって得た。

【0064】

α線スペクトルを得るためのシミュレーションで用いた模倣体および放射線源は、それぞれ図３および図４の通りである。

【0065】

３００×３００×８次元の腎臓（ｋｉｄｎｅｙ）形態の薄片ファントムを製作し、横×縦は、２９．４３×１８．１１ｍｍ^２である。ファントムの厚さは、８つのスライスを合わせて２８μｍであった。ファントムの材質は、図４に表示した通りである。領域２と領域３は、密度１．０５を有する腎臓（ｋｉｄｎｅｙ）成分であり、領域４および領域５は、密度１．０６を有する血管または血液（ｂｌｏｏｄ）成分である。当該ファントムをα線画像化装置の上に置いた状況をＧＡＴＥを用いて模倣した。

【0066】

シミュレーションで、α線放出核種である^２２５Ａｃおよび^２１３Ｂｉを、図４のように、領域１～領域５に注入した。領域２には^２２５Ａｃと^２１３Ｂｉを両方、領域３および領域５には^２２５Ａｃのみを、領域４には^２１３Ｂｉのみを注入した状況を模倣した。ベイトマン方程式（Ｂａｔｅｍａｎ ｅｑｕａｔｉｏｎ）を考慮して、放射平衡を成した状態の^２２５Ａｃおよび^２１３Ｂｉを注入したので、^２２５Ａｃのみを注入した場合にも^２１３Ｂｉを含んだ娘核種を含んでいる。^２１３Ｂｉのみを注入した場合にも^２１３Ｂｉの娘核種を含んでいる。

【0067】

注入直後の測定状況を模倣しており、総測定時間は、１５０秒であった。^２２５Ａｃの崩壊曲線にある核種がβ線とγ線を放出する場合もあるため、これらの影響を減らすために、５００ｋｅＶ以上のエネルギーを有する場合のみを計数して（＝閾値（Ｔｈｒｅｓｈｏｌｄ）を５００ｋｅＶに設定して）、エネルギースペクトルをピクセルごとに獲得した。１ピクセルのスペクトルのＹ値（係数値）を全部加えて１つの値で表現すると、１ピクセルの強度（ｉｎｔｅｎｓｉｔｙ）を示すことができる。

【0068】

獲得したピクセルの強度を２次元で獲得して画像を獲得した。当該画像は３００×３００次元で、各ピクセルの面積は１００μｍ^２であるので、全体画像の大きさは３０×３０ｍｍ^２である。得られた画像は、図５の通りであり、黒白イメージで、この画像ではピクセルごとの放射能を知ることができるだけで、核種別分布情報は分からない。

【0069】

上記のような過程を通じて得られた領域２でのα線スペクトルは、図６の通りである。図６は、数式２でのｂ_ｘ，ｙにあたる。

【0070】

各位置で、式２を利用した係数導出によって得られた親核種に対する第１α線スペクトル（赤色）と娘核種に対する第２α線スペクトル（青色）は、図７ａ～図７ｄの通りである。

【0071】

ここで、親核種は、親核種および当該領域の親核種の崩壊で生成された娘核種を指し、娘核種は、他の領域で生成された後に流入された娘核種を指す。

【0072】

図８ａおよび図８ｂは、それぞれ本発明の実験例で、^２２５Ａｃからのα線と^２１３Ｂｉからのα線画像を示したものである。

【0073】

分離結果を領域別に合算して領域別分離スペクトルを得、これを図９のように画像で示した。図９は、カラーイメージで得られる。

【0074】

図９の画像を得る方法を説明すれば次の通りである。まず、親核種と娘核種に対する全体２次元画像（図８）を得る。

【0075】

^２２５Ａｃ画像をＲＧＢ画像に変え、Ｒ：Ｇ：Ｂ＝１：０：０を適用する。値の高いピクセルは濃い赤色で、値の低いピクセルは薄い赤色で表現される。同様に、^２１３Ｂｉ画像をＲＧＢ画像に変えて、Ｒ：Ｇ：Ｂ＝０：１：０を適用する。値の高いピクセルは濃い緑色で、値の低いピクセルは薄い緑色で表現される。この２つの画像を加えると、^２２５Ａｃのみ存在する所は赤色、^２１３Ｂｉのみ存在する所は緑色で表現され、この２つが一緒に存在する領域は赤色と緑色の合成値（例えば黄色）で表現される。

【0076】

図９を利用すると、図５とは異なり、関心領域でのα線がどの放射線源から発生したものであるのかを確認することができる。

【0077】

図９を参照すると、２番領域の^２２５Ａｃと^２１３Ｂｉが全部存在する領域が赤色と緑色が混じっている黄色で表現されていることから、注入した放射線源の分布をよく推定していることを定性的に確認することができる。

【0078】

次の数式３を用いて、実際に注入された核種放射能と本発明によって推定された核種の放射能を比較した。

【0079】

【数3】

【0080】

比較の結果、図１０のように、領域別に５％未満の相対誤差のみを示した。

【0081】

上述した実施形態は、本発明を説明するための例示であって、本発明がこれに限定されるものではない。本発明が属する技術分野において通常の知識を有する者であれば、これから多様に変形して本発明を実施することが可能であるので、本発明の技術的保護範囲は、特許請求の範囲によって定められるべきである。

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7a】

【図7b】

【図7c】

【図7d】

【図8a】

【図8b】

【図9】

【図10】