JP6983129B2 - Gamma camera - Google Patents

Gamma camera Download PDF

Info

Publication number
JP6983129B2
JP6983129B2 JP2018153467A JP2018153467A JP6983129B2 JP 6983129 B2 JP6983129 B2 JP 6983129B2 JP 2018153467 A JP2018153467 A JP 2018153467A JP 2018153467 A JP2018153467 A JP 2018153467A JP 6983129 B2 JP6983129 B2 JP 6983129B2
Authority
JP
Japan
Prior art keywords
measurement target
gamma
physical property
storage unit
information
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2018153467A
Other languages
Japanese (ja)
Other versions
JP2020027081A (en
Inventor
博志 西沢
正一 中西
理 笹野
真照 林
哲史 東
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Mitsubishi Electric Corp
Original Assignee
Mitsubishi Electric Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Mitsubishi Electric Corp filed Critical Mitsubishi Electric Corp
Priority to JP2018153467A priority Critical patent/JP6983129B2/en
Publication of JP2020027081A publication Critical patent/JP2020027081A/en
Application granted granted Critical
Publication of JP6983129B2 publication Critical patent/JP6983129B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Description

本願は、放射性物質の三次元濃度分布を測定するガンマカメラに関する。 The present application relates to a gamma camera that measures a three-dimensional concentration distribution of radioactive substances.

従来、ガンマカメラは、放射性物質の濃度分布の測定に用いられるが、奥行き方向の分布を測定することが難しく、1台のガンマカメラで三次元濃度分布を測定する方法が検討されている。例えば特許文献1では、1台のガンマカメラを用い、測定対象の正面に対する垂直線から左右にそれぞれ同一の角度傾いた撮像軸上で2つの画像を撮像し、得られた2つの画像から立体画像を得ている。 Conventionally, a gamma camera has been used to measure the concentration distribution of radioactive substances, but it is difficult to measure the distribution in the depth direction, and a method of measuring a three-dimensional concentration distribution with one gamma camera has been studied. For example, in Patent Document 1, one gamma camera is used to capture two images on an imaging axis tilted at the same angle to the left and right from a vertical line with respect to the front surface of the measurement target, and a stereoscopic image is obtained from the two obtained images. Is getting.

また、ガンマ線源の三次元分布を画像化する方法として、コンプトンカメラを利用する方法が知られている。典型的なコンプトンカメラは、2つの放射線検出器を備えており、ガンマ線源からのガンマ線が前段検出器にてコンプトン散乱された後、後段検出器にて光電吸収される事象に注目し、各検出器におけるガンマ線の検出位置及び検出エネルギーに基づいてガンマ線源の三次元分布を画像化する。 Further, as a method of imaging a three-dimensional distribution of a gamma ray source, a method using a Compton camera is known. A typical Compton camera is equipped with two radiation detectors, focusing on the phenomenon that gamma rays from a gamma ray source are scattered by the pre-stage detector and then photoelectrically absorbed by the post-stage detector, and each detection is performed. The three-dimensional distribution of the gamma ray source is imaged based on the detection position and the detection energy of the gamma ray in the detector.

このコンプトンカメラを利用した測定方法として、例えば特許文献2では、放射線発生側から見て前後に並ぶ2つの検出ユニットを備えた放射線測定装置により、散乱点の位置、吸収点の位置、及び放射線エネルギーから、散乱点を頂点とする円錐面(コンプトンコーン)として放射線源位置を特定している。 As a measurement method using this Compton camera, for example, in Patent Document 2, a radiation measuring device provided with two detection units arranged in front and behind when viewed from the radiation generation side, the position of a scattering point, the position of an absorption point, and the radiation energy are used. Therefore, the position of the radiation source is specified as a conical surface (Compton cone) whose apex is the scattering point.

また、特許文献3には、コンプトンカメラによるガンマ線源の三次元分布画像の適正化を図る方法として、統計学の手法である最尤推定期待値最大化法(ML−EM法)に基づいた画像再構成法が示されている。さらに、非特許文献1には、核医学におけるコンプトン散乱カメラ画像のリストモード最尤再構成法(LM−ML−EM法)について示されている。LM−ML−EM法によれば、観測されたガンマ線の測定データから空間中のガンマ線源分布を最尤推定により求めることが困難な場合でも、反復計算により尤度の期待値を最大化するようにガンマ線源分布を求めることが可能である。 Further, Patent Document 3 describes an image based on the maximum likelihood estimation expected value maximization method (ML-EM method), which is a statistical method, as a method for optimizing the three-dimensional distribution image of a gamma ray source by a Compton camera. The reconstruction method is shown. Further, Non-Patent Document 1 describes a list mode maximum likelihood reconstruction method (LM-ML-EM method) of a Compton scattering camera image in nuclear medicine. According to the LM-ML-EM method, even if it is difficult to obtain the gamma ray source distribution in space from the observed gamma ray measurement data by maximum likelihood estimation, the expected value of likelihood is maximized by iterative calculation. It is possible to obtain the gamma ray source distribution.

特開平9−218268号公報Japanese Unexamined Patent Publication No. 9-218268 特開2011−85418号公報Japanese Unexamined Patent Publication No. 2011-85418 WO2016−035706号公報WO2016-035706 Gazette

S.J.Wilderman,N.H.Clinthorne,J.A.Fessler,and W.L.Rogers, “List−mode maximum likelihood reconstruction of Compton scatter camera image in nuclear medicine”,IEEE Nucl. Sci. Symp., 1998, vol.3,pp.1716−1720.S. J. Wilderman, N.M. H. Clinthorne, J. et al. A. Fessler, and W. L. Rogers, "List-mode maximum likelihood reconstruction of Compton scatter camera camera image in nuclear medicine", IEEE Nucle. Sci. Symp. , 1998, vol. 3, pp. 1716-1720.

従来のガンマカメラでは、放射性物質の三次元濃度分布を測定するために、放射性物質から直接入射するガンマ線のみをデータとして用い、放射性物質と放射線検出器との間で散乱されたガンマ線を検出してもデータとして用いていない。一方、測定対象が人体等に比べて高密度な物質、例えばコンクリートまたは金属等の場合、大部分のガンマ線が高密度物質の物体内部で散乱される。このため、従来のガンマカメラでは、測定対象が高密度物質である場合、放射性物質の三次元濃度分布の測定が難しいという課題があった。 In a conventional gamma camera, in order to measure the three-dimensional concentration distribution of a radioactive substance, only the gamma rays directly incident from the radioactive substance are used as data, and the gamma rays scattered between the radioactive substance and the radiation detector are detected. Is not used as data. On the other hand, when the measurement target is a substance having a higher density than that of the human body, for example, concrete or metal, most of the gamma rays are scattered inside the object of the high density substance. For this reason, the conventional gamma camera has a problem that it is difficult to measure the three-dimensional concentration distribution of a radioactive substance when the measurement target is a high-density substance.

測定対象における放射性物質の三次元濃度分布を測定する方法として、特許文献1のように、1台のガンマカメラを移動させて複数回の測定を行う方法では、ガンマカメラを駆動部と連結させ移動させる機構が必要であり、複数回の測定を行うため測定時間が長くなるという課題があった。また、特許文献2、特許文献3、及び非特許文献1の方法では、1台のガンマカメラで1回の測定によって放射性物質の三次元濃度分布を測定しているが、コンプトン散乱を使った計測方法、いわゆるコンプトンカメラにのみに適用される方法である。 As a method of measuring the three-dimensional concentration distribution of a radioactive substance in a measurement target, in a method of moving one gamma camera to perform multiple measurements as in Patent Document 1, the gamma camera is connected to a drive unit and moved. There is a problem that the measurement time becomes long because a mechanism for making the measurement is required and the measurement is performed a plurality of times. Further, in the methods of Patent Document 2, Patent Document 3, and Non-Patent Document 1, the three-dimensional concentration distribution of the radioactive substance is measured by one measurement with one gamma camera, but the measurement using Compton scattering is used. The method, a method that applies only to so-called Compton cameras.

本願は、上記のような課題を解決するための技術を開示するものであり、測定対象が高密度物質である場合も、放射性物質の三次元濃度分布を短時間で高感度に測定することが可能なガンマカメラを得ることを目的とする。 This application discloses a technique for solving the above-mentioned problems, and it is possible to measure the three-dimensional concentration distribution of a radioactive substance with high sensitivity in a short time even when the measurement target is a high-density substance. The aim is to get a possible gamma camera.

本願に開示されるガンマカメラは、測定対象から飛来するガンマ線の入射方向とエネルギーを検出する放射線検出器と、放射線検出器により取得された入射方向情報とエネルギー情報を記憶するガンマ線データ記憶部と、少なくとも測定対象の成分及び密度を含む物性情報を取得する物性情報取得装置と、物性情報取得装置により取得された測定対象の物性情報を記憶する物性情報記憶部と、放射線検出器の測定条件情報を取得する測定条件取得装置と、測定条件取得装置により取得された測定条件情報を記憶する測定条件記憶部と、ガンマ線データ記憶部に記憶された入射方向情報とエネルギー情報、物性情報記憶部に記憶された物性情報、及び測定条件記憶部に記憶された測定条件情報に基づいて、測定対象内部で散乱されたガンマ線のエネルギー変化と方向変化を計算し、放射性物質の三次元濃度分布の演算を行う三次元分布演算装置を備えたものである。 The gamma camera disclosed in the present application includes a radiation detector that detects the incident direction and energy of gamma rays coming from a measurement target, a gamma ray data storage unit that stores incident direction information and energy information acquired by the radiation detector, and a gamma ray data storage unit. A physical property information acquisition device that acquires physical property information including at least the components and densities of the measurement target, a physical property information storage unit that stores the physical property information of the measurement target acquired by the physical property information acquisition device, and measurement condition information of the radiation detector. The measurement condition acquisition device to be acquired, the measurement condition storage unit that stores the measurement condition information acquired by the measurement condition acquisition device, the incident direction information and energy information stored in the gamma ray data storage unit, and the physical property information storage unit are stored. Based on the physical property information and the measurement condition information stored in the measurement condition storage unit, the energy change and direction change of the gamma rays scattered inside the measurement target are calculated, and the three-dimensional concentration distribution of the radioactive substance is calculated. It is equipped with an elemental distribution calculation device.

本願に開示されるガンマカメラによれば、ガンマ線の入射方向情報とエネルギー情報、測定対象の物性情報、及び放射線検出器の測定条件情報に基づいて、測定対象内部で散乱されたガンマ線のエネルギー変化と方向変化を計算し、放射性物質の分布領域を限定するようにしたので、物体内部で大部分のガンマ線が散乱される高密度物質が測定対象である場合も、放射性物質の三次元濃度分布を短時間で高感度に測定することが可能である。 According to the gamma camera disclosed in the present application, the energy change of the gamma ray scattered inside the measurement target is based on the incident direction information and energy information of the gamma ray, the physical property information of the measurement target, and the measurement condition information of the radiation detector. Since the direction change is calculated and the distribution area of the radioactive substance is limited, the three-dimensional concentration distribution of the radioactive substance is shortened even when the measurement target is a high-density substance in which most of the gamma rays are scattered inside the object. It is possible to measure with high sensitivity in time.

実施の形態1に係るガンマカメラの構成図である。It is a block diagram of the gamma camera which concerns on Embodiment 1. FIG. 実施の形態1に係るガンマカメラによる放射性物質の三次元濃度分布の測定方法を説明する図である。It is a figure explaining the method of measuring the three-dimensional concentration distribution of a radioactive substance by the gamma camera which concerns on Embodiment 1. FIG. 実施の形態1に係るガンマカメラにおける三次元分布演算装置のハードウェア構成図である。It is a hardware block diagram of the 3D distribution arithmetic unit in the gamma camera which concerns on Embodiment 1. FIG. 放射性物質からのガンマ線の直接成分と散乱成分のエネルギー分布を説明する図である。It is a figure explaining the energy distribution of the direct component and the scattering component of the gamma ray from a radioactive substance. 測定対象内部におけるガンマ線の散乱確率と吸収確率の差の確率分布を説明する図である。It is a figure explaining the probability distribution of the difference between the scattering probability and the absorption probability of gamma rays in the measurement target. 実施の形態1に係るガンマカメラの三次元分布演算装置における処理の流れを説明する図である。It is a figure explaining the flow of processing in the 3D distribution arithmetic unit of the gamma camera which concerns on Embodiment 1. FIG. 実施の形態1に係るガンマカメラにおける放射性物質の三次元濃度分布の推定方法を説明する図である。It is a figure explaining the method of estimating the three-dimensional concentration distribution of a radioactive substance in the gamma camera which concerns on Embodiment 1. FIG. 実施の形態2に係るガンマカメラの構成図である。It is a block diagram of the gamma camera which concerns on Embodiment 2. FIG. 実施の形態2に係るガンマカメラの三次元分布演算装置における処理の流れを説明する図である。It is a figure explaining the flow of processing in the 3D distribution arithmetic unit of the gamma camera which concerns on Embodiment 2. FIG.

実施の形態1.
以下に、実施の形態1によるガンマカメラについて、図面に基づいて説明する。図1は、実施の形態1に係るガンマカメラの構成図、図2は、実施の形態1に係るガンマカメラによる放射性物質の三次元濃度分布の測定方法を説明する図、図3は、実施の形態1に係るガンマカメラにおける三次元分布演算装置のハードウェア構成図である。なお、各図において、同一、相当部分には同一符号を付している。
Embodiment 1.
Hereinafter, the gamma camera according to the first embodiment will be described with reference to the drawings. FIG. 1 is a configuration diagram of a gamma camera according to the first embodiment, FIG. 2 is a diagram illustrating a method of measuring a three-dimensional concentration distribution of a radioactive substance by the gamma camera according to the first embodiment, and FIG. 3 is an embodiment. It is a hardware block diagram of the 3D distribution arithmetic unit in the gamma camera which concerns on Embodiment 1. FIG. In each figure, the same and corresponding parts are designated by the same reference numerals.

実施の形態1に係るガンマカメラ100の測定対象は主に、人体等に比べて高密度な物質、例えばコンクリート、金属等である。測定対象に存在する放射性物質は、放射線を放出して崩壊し他の原子核に変わる放射性核種を含んでいる。ガンマ線(以下、γ線と記す)は、γ線核種から放出される。図2に示すように、測定対象10の内部に存在する放射性物質11から飛来するγ線は、ガンマカメラ100に直接入射する直接成分12と、測定対象10の内部で散乱された散乱成分13に分けられる。測定対象10が高密度物質の場合、大部分のγ線が測定対象内部で散乱される。 The measurement target of the gamma camera 100 according to the first embodiment is mainly a substance having a higher density than that of the human body, for example, concrete, metal, or the like. The radioactive material present in the measurement target contains radionuclides that emit radiation and decay to change into other nuclei. Gamma rays (hereinafter referred to as γ rays) are emitted from gamma ray nuclides. As shown in FIG. 2, the γ-rays coming from the radioactive substance 11 existing inside the measurement target 10 are directed to the direct component 12 directly incident on the gamma camera 100 and the scattered component 13 scattered inside the measurement target 10. Divided. When the measurement target 10 is a high-density substance, most of the γ-rays are scattered inside the measurement target.

ガンマカメラ100は、図1に示すように、放射線検出器1、ガンマ線データ記憶部2(以下、γ線データ記憶部2と記す)、物性情報取得装置3、物性情報記憶部4、測定条件取得装置5、測定条件記憶部6、三次元分布演算装置7、及び表示部8を備えている。放射線検出器1は、測定対象に存在する放射性物質11から飛来するγ線の入射方向とエネルギーを測定する。この時、放射線検出器1が検出するγ線は、直接成分12と散乱成分13の両方を含む。 As shown in FIG. 1, the gamma camera 100 includes a radiation detector 1, a gamma ray data storage unit 2 (hereinafter referred to as a gamma ray data storage unit 2), a physical property information acquisition device 3, a physical property information storage unit 4, and measurement condition acquisition. It includes a device 5, a measurement condition storage unit 6, a three-dimensional distribution calculation device 7, and a display unit 8. The radiation detector 1 measures the incident direction and energy of the γ-rays coming from the radioactive substance 11 existing in the measurement target. At this time, the γ-ray detected by the radiation detector 1 includes both the direct component 12 and the scattering component 13.

放射線検出器1は、γ線の入射によりγ線が反応した位置と、エネルギーに応じて決まる電圧もしくは電流強度を出力するγ線センサを有している。γ線センサとしては、半導体式検出器を複数配置したもの、半導体式検出器において信号を取り出す電極を複数有するもの、γ線を紫外光、可視光、及び赤外光に変換するシンチレータと光検出部を含むユニットを複数配置したもの等が用いられる。 The radiation detector 1 has a γ-ray sensor that outputs a voltage or current intensity determined by energy and a position where the γ-ray reacts due to the incident of the γ-ray. Gamma-ray sensors include those with multiple semiconductor detectors, those with multiple electrodes that extract signals from semiconductor detectors, and scintillators and photodetectors that convert γ-rays into ultraviolet light, visible light, and infrared light. A unit in which a plurality of units including a unit are arranged is used.

半導体式検出器としては、Si、CdTe、CdZnTe等を用いた検出器が用いられる。半導体式検出器の場合、γ線が反応した半導体式検出器が配置された位置、または信号が取り出された電極の位置から、γ線が反応した位置が検出される。一方、シンチレータは、NaI(Tl)、CsI(Tl)、BGO、LSO等のシンチレーション材料を含み、シンチレータからの光を光検出部で検出する。シンチレータの場合、シンチレーション材料からの光を検出した位置から、γ線が反応した位置が検出される。 As the semiconductor type detector, a detector using Si, CdTe, CdZnTe or the like is used. In the case of a semiconductor detector, the position where the γ-ray reacts is detected from the position where the semiconductor detector with which the γ-ray reacts is arranged or the position of the electrode from which the signal is taken out. On the other hand, the scintillator contains a scintillation material such as NaI (Tl), CsI (Tl), BGO, LSO, etc., and the light from the scintillator is detected by the photodetector. In the case of a scintillator, the position where the γ-ray reacts is detected from the position where the light from the scintillation material is detected.

また、放射線検出器1は、γ線の入射方向を特定するための構成を備えている。例えばγ線センサと、γ線を遮蔽する遮蔽体に一つ以上の穴を空けた筐体とを組み合わせた構成、またはγ線センサと、γ線の入射方向を制限するコリメータとを組み合わせた構成、またはγ線センサを二つ以上並べた構成等がこれに該当する。γ線データ記憶部2は、検出エネルギー記憶部21と入射方向記憶部22を有し、放射線検出器1により検出されたγ線のエネルギー情報と入射方向情報をそれぞれ記憶する。 Further, the radiation detector 1 has a configuration for specifying the incident direction of the γ-ray. For example, a configuration in which a γ-ray sensor is combined with a housing having one or more holes in a shield that shields γ-rays, or a configuration in which a γ-ray sensor and a collimeter that limits the incident direction of γ-rays are combined. , Or a configuration in which two or more γ-ray sensors are arranged side by side corresponds to this. The γ-ray data storage unit 2 has a detection energy storage unit 21 and an incident direction storage unit 22, and stores γ-ray energy information and incident direction information detected by the radiation detector 1, respectively.

物性情報取得装置3は、測定対象の物性情報として、測定対象の密度、密度分布、成分(構成元素)、構成元素の混合比、及び厚さ寸法等を測定する。物性情報取得装置3により取得される物性情報には、少なくとも測定対象の成分及び密度が含まれる。物性情報取得装置3は、エックス線検査装置、γ線検査装置、中性子線検査装置、超音波検査装置、及び音波検査装置の中から選ばれた複数の装置を含む。なお、上記以外の装置を用いることも可能である。 The physical property information acquisition device 3 measures the density, density distribution, component (constituent element), mixing ratio of the constituent elements, thickness dimension, etc. of the measurement target as the physical property information of the measurement target. The physical property information acquired by the physical property information acquisition device 3 includes at least the components and densities to be measured. The physical property information acquisition device 3 includes a plurality of devices selected from an X-ray inspection device, a γ-ray inspection device, a neutron beam inspection device, an ultrasonic inspection device, and a sound wave inspection device. It is also possible to use a device other than the above.

測定対象の構成元素及び混合比は、エックス線検査装置またはγ線検査装置により測定される。測定対象の密度及び密度分布は中性子線検査装置により測定され、測定対象の厚さ寸法は超音波検査装置または音波検査装置により測定される。物性情報記憶部4は、物性情報取得装置3により取得された物性情報を記憶する。なお、物性情報取得装置3による物性情報の測定は、放射性物質の三次元濃度分布の測定と同時である必要はなく、事前に測定対象の一部を採取して物性情報を測定し、物性情報記憶部4に記録しておくことも可能である。 The constituent elements to be measured and the mixing ratio are measured by an X-ray inspection device or a γ-ray inspection device. The density and density distribution of the measurement target are measured by a neutron beam inspection device, and the thickness dimension of the measurement target is measured by an ultrasonic inspection device or a sound wave inspection device. The physical property information storage unit 4 stores the physical property information acquired by the physical property information acquisition device 3. The measurement of the physical property information by the physical property information acquisition device 3 does not have to be performed at the same time as the measurement of the three-dimensional concentration distribution of the radioactive substance. It is also possible to record in the storage unit 4.

測定条件取得装置5は、放射線検出器1の測定条件情報を取得する。測定条件情報には、放射線検出器1と測定対象の距離、及び放射線検出器1の視野角が含まれる。測定条件取得装置5は、レーザ距離計、電波距離計、超音波距離計、及び巻尺のいずれかを含む。ただし、上記以外の装置を用いることも可能である。測定条件記憶部6は、測定条件取得装置5により取得された測定条件情報を記憶する。 The measurement condition acquisition device 5 acquires the measurement condition information of the radiation detector 1. The measurement condition information includes the distance between the radiation detector 1 and the measurement target, and the viewing angle of the radiation detector 1. The measurement condition acquisition device 5 includes any one of a laser range finder, a radio wave range finder, an ultrasonic range finder, and a tape measure. However, it is also possible to use a device other than the above. The measurement condition storage unit 6 stores the measurement condition information acquired by the measurement condition acquisition device 5.

三次元分布演算装置7は、γ線データ記憶部2に記憶されたγ線データ、物性情報記憶部4に記憶された物性情報、及び測定条件記憶部6に記憶された測定条件情報を取得し、それらに基づいて測定対象の内部で散乱されたγ線のエネルギー変化と方向変化を計算し、放射性物質の三次元濃度分布の演算を行う。具体的には、測定対象の物性情報に基づいて推定した放射性核種から放出される放射線エネルギーを算出し、算出した放射線エネルギーと、検出エネルギー記憶部21から取得したエネルギー情報とに基づいて、γ線のエネルギー変化を計算する。 The three-dimensional distribution calculation device 7 acquires the γ-ray data stored in the γ-ray data storage unit 2, the physical property information stored in the physical property information storage unit 4, and the measurement condition information stored in the measurement condition storage unit 6. Based on these, the energy change and direction change of the γ-rays scattered inside the measurement target are calculated, and the three-dimensional concentration distribution of the radioactive substance is calculated. Specifically, the radiation energy emitted from the radioactive nuclei estimated based on the physical property information of the measurement target is calculated, and the γ-ray is based on the calculated radiation energy and the energy information acquired from the detected energy storage unit 21. Calculate the energy change of.

また、測定対象内部で散乱されたγ線のエネルギー変化からγ線が測定対象内部で散乱された角度θを算出し、算出した角度θと、入射方向記憶部22から取得した入射方向情報とに基づいて、放射性物質の分布領域を限定する。また、測定対象の物性情報に基づいて測定対象内部におけるγ線の散乱位置を推定し、推定した散乱位置による重み付けにより放射性物質の分布領域をさらに限定する。このときの演算方法の一つとして、分布が最も確からしくなるように放射性物質の三次元濃度分布を推定するML−EM法または最大エントロピー法等を適用することが可能である。 Further, the angle θ in which the γ-rays are scattered inside the measurement target is calculated from the energy change of the γ-rays scattered inside the measurement target, and the calculated angle θ and the incident direction information acquired from the incident direction storage unit 22 are used. Based on this, the distribution area of radioactive substances is limited. In addition, the scattering position of γ-rays inside the measurement target is estimated based on the physical property information of the measurement target, and the distribution region of the radioactive substance is further limited by weighting according to the estimated scattering position. As one of the calculation methods at this time, it is possible to apply the ML-EM method or the maximum entropy method for estimating the three-dimensional concentration distribution of the radioactive substance so that the distribution is most reliable.

三次元分布演算装置7による上記の演算処理は、放射線検出器1のγ線センサが検出した一つのγ線に対して行われる。放射線検出器1が複数のγ線を検出した場合、三次元分布演算装置7は、複数のγ線それぞれの入射方向情報とエネルギー情報に基づいて放射性物質の分布領域を推定し、さらに推定した複数の分布領域が重なり合った部分を放射性物質の分布領域とする。三次元分布演算装置7による演算結果は、表示部8に表示される。 The above arithmetic processing by the three-dimensional distribution arithmetic unit 7 is performed on one γ-ray detected by the γ-ray sensor of the radiation detector 1. When the radiation detector 1 detects a plurality of γ-rays, the three-dimensional distribution calculation device 7 estimates the distribution region of the radioactive substance based on the incident direction information and the energy information of each of the plurality of γ-rays, and further estimates the plurality of γ-rays. The part where the distribution areas of the above overlap is defined as the distribution area of radioactive substances. The calculation result by the three-dimensional distribution calculation device 7 is displayed on the display unit 8.

なお、三次元分布演算装置7は、図3に示すプロセッサ71とメモリ72を含む処理回路70を有している。メモリ72は、ランダムアクセスメモリ等の揮発性記憶装置と、フラッシュメモリ等の不揮発性の補助記憶装置とを具備する。また、フラッシュメモリの代わりにハードディスクの補助記憶装置を具備してもよい。プロセッサ71は、メモリ72から入力されたプログラムを実行する。この場合、補助記憶装置から揮発性記憶装置を介してプロセッサ71にプログラムが入力される。 The three-dimensional distribution arithmetic unit 7 has a processing circuit 70 including the processor 71 and the memory 72 shown in FIG. The memory 72 includes a volatile storage device such as a random access memory and a non-volatile auxiliary storage device such as a flash memory. Further, the auxiliary storage device of the hard disk may be provided instead of the flash memory. The processor 71 executes the program input from the memory 72. In this case, a program is input from the auxiliary storage device to the processor 71 via the volatile storage device.

次に、三次元分布演算装置7による演算方法について詳細に説明する。放射性物質から放出されたγ線は、散乱によりエネルギー変化と方向変化を起こしている。三次元分布演算装置7は、γ線の散乱によるエネルギー変化と方向変化を計算し、その計算結果に基づいて放射性物質の分布領域を限定する。 Next, the calculation method by the three-dimensional distribution calculation device 7 will be described in detail. Gamma rays emitted from radioactive materials cause energy changes and direction changes due to scattering. The three-dimensional distribution calculation device 7 calculates the energy change and the direction change due to the scattering of γ-rays, and limits the distribution area of the radioactive substance based on the calculation result.

図4は、放射性物質からのγ線の直接成分と散乱成分のエネルギー分布を示している。図4において、縦軸はカウント数、横軸はエネルギーであり、実線Aは直接成分のエネルギー分布、点線Bは散乱成分のエネルギー分布をそれぞれ示している。図4に示すように、γ線の直接成分は所定のエネルギーを有するが、散乱成分のエネルギーは直接成分のエネルギーよりも小さく且つ一定でない。このようなエネルギー分布の違いから、放射線検出器1で検出されたγ線が直接成分であるか散乱成分であるかを判別することができる。 FIG. 4 shows the energy distribution of the direct component and the scattered component of γ-rays from a radioactive substance. In FIG. 4, the vertical axis represents the number of counts, the horizontal axis represents the energy, the solid line A shows the energy distribution of the direct component, and the dotted line B shows the energy distribution of the scattering component. As shown in FIG. 4, the direct component of the γ-ray has a predetermined energy, but the energy of the scattering component is smaller and not constant than the energy of the direct component. From such a difference in energy distribution, it is possible to determine whether the γ-ray detected by the radiation detector 1 is a direct component or a scattering component.

三次元分布演算装置7は、物性情報記憶部4に記憶された物性情報に基づいて、測定対象内部の箇所毎にγ線が散乱される確率pを算出する。確率pの算出方法として一般的に用いられるのは、測定対象の密度ρ、γ線が測定対象内部を進んだ距離L、及び測定対象を構成する元素がγ線を散乱する散乱断面積σによる式1の関数である。 The three-dimensional distribution calculation device 7 calculates the probability p that γ-rays are scattered at each location inside the measurement target based on the physical property information stored in the physical property information storage unit 4. Generally used as a method for calculating the probability p is the density ρ of the measurement target, the distance L in which the γ-ray travels inside the measurement target, and the scattering cross section σ in which the elements constituting the measurement target scatter the γ-ray. It is a function of Equation 1.

p=1−exp(−ρσL) (式1) p = 1-exp (-ρσL) (Equation 1)

測定対象が複数の元素の集合体で構成される場合には、それぞれの元素の集合体での密度ρ、γ線がそれぞれの元素の集合体を進んだ距離L、及び散乱断面積σによる下の式2によって、測定対象内部の箇所毎に散乱される確率pが算出される。 If the measurement object is composed of a set of a plurality of elements, the density [rho i of a collection of the respective elements, the distance γ rays advances the assembly of the respective elements L i, and the scattering cross section σ The probability p of being scattered at each location inside the measurement target is calculated by the following equation 2 according to i.

p=1−exp(Σ(−ρσ)) (式2) p = 1-exp (Σ i (−ρ i σ i L i )) (Equation 2)

なお、γ線が散乱される確率pの算出方法は、式1及び式2を用いた上記の方法に限定されるものではない。別の算出方法として、物性情報取得装置3で取得した測定対象の物性情報を用い、乱数計算によって測定対象の物性情報を再現した条件でシミュレーションを行い、γ線が散乱される確率pを算出する方法がある。 The method for calculating the probability p at which γ-rays are scattered is not limited to the above method using Equations 1 and 2. As another calculation method, using the physical property information of the measurement target acquired by the physical property information acquisition device 3, a simulation is performed under the condition that the physical property information of the measurement target is reproduced by random number calculation, and the probability p that γ-rays are scattered is calculated. There is a way.

同様に、式1及び式2を用いた上記の方法により、物性情報記憶部4から取得した物性情報に基づいて、測定対象内部の箇所毎にγ線が吸収される確率を算出することが可能である。その際は、測定対象を構成する元素がγ線を吸収する断面積を式1及び式2に適用する。同様に、γ線が散乱される確率pを算出するシミュレーションによっても、γ線が吸収される確率を算出することができる。 Similarly, by the above method using the formulas 1 and 2, it is possible to calculate the probability that γ-rays are absorbed at each location inside the measurement target based on the physical property information acquired from the physical property information storage unit 4. Is. In that case, the cross-sectional area in which the element constituting the measurement target absorbs γ-rays is applied to the formulas 1 and 2. Similarly, the probability that γ-rays are absorbed can also be calculated by a simulation that calculates the probability p that γ-rays are scattered.

図5は、測定対象内部におけるγ線の散乱確率と吸収確率を示している。ここで、散乱確率とは、測定対象内部の放射性物質から放射されるγ線が測定対象内部で散乱される確率であり、吸収確率とは、測定対象内部の放射性物質から放射されるγ線が測定対象内部で吸収される確率である。図5において、縦軸は確率、横軸は放射線(γ線)が測定対象内部を進んだ距離(L)であり、実線Sは散乱確率、点線Aは吸収確率、一点鎖線Dは散乱確率と吸収確率の差をそれぞれ示している。 FIG. 5 shows the scattering probability and the absorption probability of γ-rays inside the measurement target. Here, the scattering probability is the probability that γ-rays emitted from the radioactive substance inside the measurement target are scattered inside the measurement target, and the absorption probability is the γ-ray emitted from the radioactive material inside the measurement target. It is the probability of being absorbed inside the measurement target. In FIG. 5, the vertical axis is the probability, the horizontal axis is the distance (L) that the radiation (γ ray) has traveled inside the measurement target, the solid line S is the scattering probability, the dotted line A is the absorption probability, and the alternate long and short dash line D is the scattering probability. The difference in absorption probability is shown respectively.

図5に示すように、測定対象内部においてγ線が進んだ距離が長くなると、γ線は100%散乱され、さらに距離が長くなるとγ線は100%吸収される。このため、測定対象内部においてγ線が進んだ距離が長くなると、散乱確率と吸収確率はいずれも1になる。測定対象内部でγ線が100%吸収されると、放射線検出器1はγ線を検出することができない。 As shown in FIG. 5, when the distance traveled by the γ-rays is long inside the measurement target, 100% of the γ-rays are scattered, and when the distance is further long, 100% of the γ-rays are absorbed. Therefore, when the distance traveled by the γ-rays becomes long inside the measurement target, both the scattering probability and the absorption probability become 1. If 100% of the γ-rays are absorbed inside the measurement target, the radiation detector 1 cannot detect the γ-rays.

測定対象内部のある箇所において、その位置でγ線が散乱される確率は、散乱確率と吸収確率の差(図4中、Dで示す)のような確率分布となる。従って、測定対象内部でγ線が散乱された箇所の位置情報は、散乱確率と吸収確率の差の確率分布を用いて推定することができる。測定対象が複数の元素の集合体で構成される場合には、各々の元素に対応する散乱確率と吸収確率の差の確率分布が用いられる。散乱確率と吸収確率の差の確率分布により推定された位置情報は、放射性物質の分布領域を限定するための重み付けに用いられる。 The probability that γ-rays are scattered at a certain place inside the measurement target has a probability distribution such as the difference between the scattering probability and the absorption probability (indicated by D in FIG. 4). Therefore, the position information of the location where the γ-rays are scattered inside the measurement target can be estimated by using the probability distribution of the difference between the scattering probability and the absorption probability. When the measurement target is composed of an aggregate of a plurality of elements, the probability distribution of the difference between the scattering probability and the absorption probability corresponding to each element is used. The position information estimated by the probability distribution of the difference between the scattering probability and the absorption probability is used for weighting to limit the distribution region of the radioactive material.

また、測定対象内部でγ線が散乱された場合には、コンプトン散乱によるγ線のエネルギー損失が起き、エネルギー損失後のγ線のエネルギーE´は、散乱される前のエネルギーE、電子の静止エネルギーE、及び散乱された角度θにより、下の式3で示される。式3において、Eは散乱される前のγ線のエネルギー、E´はコンプトン散乱によるエネルギー損失後のγ線のエネルギー、θは放射性物質から放出されたγ線が測定対象内部で散乱された角度である。 Further, when γ-rays are scattered inside the measurement target, energy loss of γ-rays occurs due to Compton scattering, and the energy E'of γ-rays after the energy loss is the energy E before being scattered and the electron stationary. It is represented by the following equation 3 by the energy E 0 and the scattered angle θ. In Equation 3, E is the energy of the γ-rays before being scattered, E'is the energy of the γ-rays after the energy loss due to Compton scattering, and θ is the angle at which the γ-rays emitted from the radioactive substance are scattered inside the measurement target. Is.

E´=E/(1+E×(1−cosθ)/E) (式3) E'= E / (1 + E × (1-cosθ) / E 0 ) (Equation 3)

式3において、散乱される前のエネルギーEは、測定対象の物性情報から推定したγ線核種から放出される放射線エネルギーから求められる。また、エネルギー損失後のγ線のエネルギーE´は、放射線検出器1により検出されたγ線のエネルギーであり、エネルギー情報として検出エネルギー記憶部21に記憶されている。この式3により、放射性物質から放出されたγ線が測定対象内部で散乱された角度θが求められる。 In Equation 3, the energy E before being scattered is obtained from the radiation energy emitted from the γ-ray nuclide estimated from the physical property information of the measurement target. Further, the γ-ray energy E'after the energy loss is the γ-ray energy detected by the radiation detector 1, and is stored in the detection energy storage unit 21 as energy information. From this equation 3, the angle θ at which the γ-rays emitted from the radioactive substance are scattered inside the measurement target can be obtained.

実施の形態1に係るガンマカメラ100の三次元分布演算装置7における処理の流れについて、図6のフローチャートを用いて説明する。まず、ステップS1において、放射線検出器1により検出されたγ線のエネルギー情報を、検出エネルギー記憶部21から取得する。続いてステップS2において、物性情報記憶部4から取得した測定対象の物性情報に基づいて、測定対象に存在する放射性核種(γ線核種)を推定する。さらにステップS3において、推定したγ線核種から放出される放射線エネルギーを算出する。 The flow of processing in the three-dimensional distribution calculation device 7 of the gamma camera 100 according to the first embodiment will be described with reference to the flowchart of FIG. First, in step S1, the energy information of the γ-rays detected by the radiation detector 1 is acquired from the detection energy storage unit 21. Subsequently, in step S2, the radionuclide (γ-ray nuclide) existing in the measurement target is estimated based on the physical property information of the measurement target acquired from the physical property information storage unit 4. Further, in step S3, the radiation energy emitted from the estimated γ-ray nuclide is calculated.

続いてステップS4において、ステップS1で取得したγ線のエネルギー情報と、ステップS3で算出した推定γ線核種から放出される放射線エネルギーを用い、放射性物質から放出されたγ線が測定対象内部で散乱された角度θを前述の方法で算出する。次にステップS5において、放射線検出器1により検出された一つのγ線の入射方向情報を、入射方向記憶部22から取得する。 Subsequently, in step S4, the γ-rays emitted from the radioactive material are scattered inside the measurement target using the γ-ray energy information acquired in step S1 and the radiation energy emitted from the estimated γ-ray nuclide calculated in step S3. The calculated angle θ is calculated by the above method. Next, in step S5, the incident direction information of one γ-ray detected by the radiation detector 1 is acquired from the incident direction storage unit 22.

さらにステップS6において、ステップS4で算出された角度θと、ステップS5で取得したγ線の入射方向情報を用いて、放射性物質の分布領域を限定するための演算を行う。式3により計算された角度θと、放射線検出器1により取得したγ線の入射方向情報から求めた放射性物質の分布領域は、円錐形状となる。すなわち、ステップS6において、放射性物質の分布領域は、円錐形状の内部の領域に限定される。 Further, in step S6, an calculation for limiting the distribution region of the radioactive substance is performed using the angle θ calculated in step S4 and the incident direction information of the γ-rays acquired in step S5. The distribution region of the radioactive substance obtained from the angle θ calculated by the equation 3 and the incident direction information of the γ-rays acquired by the radiation detector 1 has a conical shape. That is, in step S6, the distribution region of the radioactive substance is limited to the region inside the cone shape.

続いてステップS7及びステップS8において、物性情報記憶部4から取得した測定対象の物性情報に基づいて、γ線の散乱確率と吸収確率を、前述の方法でそれぞれ算出する。最後にステップS9において、ステップS6において限定された放射性物質の分布領域に対し、ステップS7及びステップS8で算出した散乱確率と吸収確率による重み付けを前述の方法で行い、放射性物質の分布領域をさらに限定する。 Subsequently, in steps S7 and S8, the scattering probability and the absorption probability of γ-rays are calculated by the above-mentioned methods, respectively, based on the physical property information of the measurement target acquired from the physical property information storage unit 4. Finally, in step S9, the distribution region of the radioactive material limited in step S6 is weighted by the scattering probability and the absorption probability calculated in steps S7 and S8 by the above-mentioned method, and the distribution region of the radioactive material is further limited. do.

なお、図6のフローチャートにおいて、ステップS1からステップS9までの処理の順番は、適宜入れ替えが可能である。ステップS1からステップS5の処理は、ステップS6の処理の前に行われていれば良く、ステップS1の処理をステップS2とステップS3の処理の後に行っても良い。同様に、ステップS7及びステップS8の処理は、ステップS9の処理の前に行われていれば良い。 In the flowchart of FIG. 6, the order of processing from step S1 to step S9 can be changed as appropriate. The processing of steps S1 to S5 may be performed before the processing of step S6, and the processing of step S1 may be performed after the processing of steps S2 and S3. Similarly, the processing of step S7 and step S8 may be performed before the processing of step S9.

以上のステップS1からステップS9の処理は、放射線検出器1のγ線センサが検出した一つのγ線に対して行われる。γ線センサが二つ以上のγ線を検出した場合は、それぞれのγ線によって演算された放射性物質の三次元濃度分布を重ね合わせることが可能となる。ガンマカメラ100が複数のγ線を検出したときの放射性物質の三次元濃度分布の推定方法について、図7を用いて説明する。図7(a)は、ガンマカメラと測定対象が紙面に対し垂直方向に配置された状態を示し、図7(b)は、ガンマカメラと測定対象が紙面に対し平行に配置された状態を示している。 The processing of steps S1 to S9 is performed on one γ-ray detected by the γ-ray sensor of the radiation detector 1. When the γ-ray sensor detects two or more γ-rays, it is possible to superimpose the three-dimensional concentration distribution of the radioactive substance calculated by each γ-ray. A method of estimating the three-dimensional concentration distribution of the radioactive substance when the gamma camera 100 detects a plurality of gamma rays will be described with reference to FIG. 7. FIG. 7A shows a state in which the gamma camera and the measurement target are arranged in the direction perpendicular to the paper surface, and FIG. 7B shows a state in which the gamma camera and the measurement target are arranged parallel to the paper surface. ing.

図7(a)及び図7(b)において、二点鎖線T1はガンマカメラ100に検出された1つ目のγ線の軌跡、一点鎖線T2はガンマカメラ100に検出された2つ目のγ線の軌跡を示し、A1は1つ目のγ線から推定された放射性物質の分布領域、A2は2つ目のγ線から推定された放射性物質の分布領域、さらに、A3は2つのγ線から推定された放射性物質の分布領域を示している。 In FIGS. 7 (a) and 7 (b), the two-point chain line T1 is the locus of the first γ-ray detected by the gamma camera 100, and the one-point chain line T2 is the second γ detected by the gamma camera 100. The locus of the line is shown, A1 is the distribution area of the radioactive substance estimated from the first γ-ray, A2 is the distribution area of the radioactive substance estimated from the second γ-ray, and A3 is the distribution area of the two γ-rays. The distribution area of radioactive substances estimated from is shown.

三次元分布演算装置7は、1つ目のγ線のエネルギー変化と方向変化の演算結果に基づいて放射性物質の分布領域を限定し(図6のステップS6)、さらに散乱確率と吸収確率による重み付けを行い(図6のステップS9)、放射性物質の分布領域A1を推定する。同様に、2つ目のγ線によって分布領域A2を推定し、これらの分布領域A1、A2を重ね合わせることにより、放射性物質の分布領域A3を推定する。真に放射性物質が存在する領域では、それぞれのγ線から演算された放射性物質の分布領域の重なりが増える。すなわち、この重なりが多いほど、真の放射性物質の分布領域である確率が高くなり、三次元濃度分布の測定精度が向上する。 The three-dimensional distribution calculation device 7 limits the distribution region of the radioactive material based on the calculation results of the energy change and the direction change of the first γ-ray (step S6 in FIG. 6), and further weights by the scattering probability and the absorption probability. (Step S9 in FIG. 6) to estimate the distribution region A1 of the radioactive substance. Similarly, the distribution region A2 is estimated by the second γ-ray, and the distribution region A3 of the radioactive substance is estimated by superimposing these distribution regions A1 and A2. In the region where the radioactive material is truly present, the overlap of the distribution region of the radioactive material calculated from each γ-ray increases. That is, the greater the overlap, the higher the probability that the region is a true radioactive substance distribution region, and the higher the measurement accuracy of the three-dimensional concentration distribution.

以上のように、実施の形態1に係るガンマカメラ100によれば、γ線の入射方向情報とエネルギー情報、測定対象の物性情報、及び放射線検出器1の測定条件情報に基づいて、測定対象内部で散乱されたγ線のエネルギー変化と方向変化を計算し、放射性物質の分布領域を限定するようにしたので、物体内部で大部分のγ線が散乱される高密度物質が測定対象である場合も、放射性物質の三次元濃度分布を短時間で高感度に測定することが可能である。 As described above, according to the gamma camera 100 according to the first embodiment, the inside of the measurement target is based on the incident direction information and energy information of the γ-ray, the physical property information of the measurement target, and the measurement condition information of the radiation detector 1. Since the energy change and direction change of the γ-rays scattered in the above are calculated and the distribution area of the radioactive material is limited, the measurement target is a high-density material in which most of the γ-rays are scattered inside the object. However, it is possible to measure the three-dimensional concentration distribution of radioactive substances in a short time with high sensitivity.

また、測定対象の物性情報に基づいて算出した散乱確率と吸収確率による重み付けにより、放射性物質の分布領域をさらに限定するようにしたので、測定精度が高くなる。また、複数のγ線を検出した場合には、それぞれのγ線から推定した放射性物質の分布領域を重ね合わせることにより、さらに測定精度が向上する。 Further, since the distribution region of the radioactive substance is further limited by the weighting based on the scattering probability and the absorption probability calculated based on the physical property information of the measurement target, the measurement accuracy is improved. Further, when a plurality of γ-rays are detected, the measurement accuracy is further improved by superimposing the distribution regions of the radioactive substances estimated from the respective γ-rays.

実施の形態2.
図8は、実施の形態2に係るガンマカメラの構成図を示している。実施の形態2に係るガンマカメラ100Aは、上記実施の形態1に係るガンマカメラ100に、測定対象の形状を取得する形状情報取得装置14と、形状情報取得装置14により取得された形状情報を記憶する形状情報記憶部15とを付加したものである。形状情報取得装置14としては、三次元形状測定装置、または輪郭形状測定装置等が用いられる。その他の構成については上記実施の形態1と同様であるため、ここでは説明を省略する。
Embodiment 2.
FIG. 8 shows a configuration diagram of a gamma camera according to the second embodiment. The gamma camera 100A according to the second embodiment stores the shape information acquisition device 14 for acquiring the shape of the measurement target and the shape information acquired by the shape information acquisition device 14 in the gamma camera 100 according to the first embodiment. The shape information storage unit 15 is added. As the shape information acquisition device 14, a three-dimensional shape measuring device, a contour shape measuring device, or the like is used. Since other configurations are the same as those in the first embodiment, the description thereof will be omitted here.

上記実施の形態1に係るガンマカメラ100は、測定対象以外の領域からのγ線を含めて演算を行っている。これに対し、実施の形態2に係るガンマカメラ100Aの三次元分布演算装置7Aは、形状情報記憶部15から取得した測定対象の形状情報に基づいて、放射性物質の分布領域を補正することができる。具体的には、演算により得られた放射性物質の分布領域と、形状情報取得装置14により取得した測定対象の形状情報とを比較し、測定対象以外の領域について演算した結果を取り除く。 The gamma camera 100 according to the first embodiment performs a calculation including γ-rays from a region other than the measurement target. On the other hand, the three-dimensional distribution calculation device 7A of the gamma camera 100A according to the second embodiment can correct the distribution region of the radioactive substance based on the shape information of the measurement target acquired from the shape information storage unit 15. .. Specifically, the distribution region of the radioactive substance obtained by the calculation is compared with the shape information of the measurement target acquired by the shape information acquisition device 14, and the calculation result for the region other than the measurement target is removed.

実施の形態2に係るガンマカメラ100Aの三次元分布演算装置7Aにおける処理の流れについて、図9のフローチャートを用いて説明する。なお、図9のステップS11からステップS18の処理は、上記実施の形態1で説明した図6のフローチャートのステップS1からステップS8までの処理と同様であるので、ここでは説明を省略する。図9のステップS19では、ステップS15で取得したγ線の入射方向情報、ステップS17で算出したγ線の散乱確率、及びステップS18で算出したγ線の吸収確率に基づいて、測定対象内部におけるγ線の散乱位置を推定する。 The flow of processing in the three-dimensional distribution calculation device 7A of the gamma camera 100A according to the second embodiment will be described with reference to the flowchart of FIG. Since the processes of steps S11 to S18 of FIG. 9 are the same as the processes of steps S1 to S8 of the flowchart of FIG. 6 described in the first embodiment, the description thereof will be omitted here. In step S19 of FIG. 9, γ inside the measurement target is based on the incident direction information of γ-rays acquired in step S15, the scattering probability of γ-rays calculated in step S17, and the absorption probability of γ-rays calculated in step S18. Estimate the scattering position of the line.

続いてステップS20において、ステップS16において限定された放射性物質の分布領域に対し、ステップS19で推定したγ線の散乱位置による重み付けを上記実施の形態1と同様の方法で行い、放射性物質の分布領域をさらに限定する。次にステップS21において、形状情報記憶部15に記憶された測定対象の形状情報を取得する。最後にステップS22において、ステップS21で取得した測定対象の形状情報と、ステップS20で限定した放射性物質の分布領域とを比較し、分布領域を補正する。 Subsequently, in step S20, the distribution region of the radioactive substance limited in step S16 is weighted by the scattering position of the γ-rays estimated in step S19 by the same method as in the first embodiment, and the distribution region of the radioactive substance is performed. Is further limited. Next, in step S21, the shape information of the measurement target stored in the shape information storage unit 15 is acquired. Finally, in step S22, the shape information of the measurement target acquired in step S21 is compared with the distribution region of the radioactive substance limited in step S20, and the distribution region is corrected.

なお、ステップS11からステップS22までの処理の順番は、適宜入れ替えが可能である。ステップS11からステップS15の処理は、ステップS16の処理の前に行われていれば良く、ステップS11の処理をステップS12とステップS13の処理の後に行っても良い。また、ステップS21の形状情報取得処理をステップS16の処理の前に行い、測定対象以外の領域での放射性物質の三次元濃度分布の演算を行わないようにしてもよい。 The order of processing from step S11 to step S22 can be changed as appropriate. The processing of steps S11 to S15 may be performed before the processing of step S16, and the processing of step S11 may be performed after the processing of steps S12 and S13. Further, the shape information acquisition process in step S21 may be performed before the process in step S16 so that the calculation of the three-dimensional concentration distribution of the radioactive substance in the region other than the measurement target may not be performed.

実施の形態2によれば、上記実施の形態1と同様の効果に加え、演算により得られた放射性物質の分布領域と、形状情報取得装置14により取得した測定対象の形状情報を比較して分布領域の補正を行うようにしたので、測定対象以外の領域について演算した結果を取り除くことができ、測定精度が向上する。また、測定対象以外の領域での放射性物質の三次元濃度分布の演算を行わないようにできるため、上記実施の形態1よりも演算に要する時間を短縮することが可能となる。 According to the second embodiment, in addition to the same effect as that of the first embodiment, the distribution region of the radioactive substance obtained by the calculation and the shape information of the measurement target acquired by the shape information acquisition device 14 are compared and distributed. Since the area is corrected, the result of calculation for the area other than the measurement target can be removed, and the measurement accuracy is improved. Further, since the calculation of the three-dimensional concentration distribution of the radioactive substance in the region other than the measurement target can be prevented, the time required for the calculation can be shortened as compared with the first embodiment.

本開示は、様々な例示的な実施の形態及び実施例が記載されているが、1つ、または複数の実施の形態に記載された様々な特徴、態様、及び機能は特定の実施の形態の適用に限られるのではなく、単独で、または様々な組み合わせで実施の形態に適用可能である。従って、例示されていない無数の変形例が、本願明細書に開示される技術の範囲内において想定される。例えば、少なくとも1つの構成要素を変形する場合、追加する場合または省略する場合、さらには、少なくとも1つの構成要素を抽出し、他の実施の形態の構成要素と組み合わせる場合が含まれるものとする。 The present disclosure describes various exemplary embodiments and examples, although the various features, embodiments, and functions described in one or more embodiments are those of a particular embodiment. It is not limited to application, but can be applied to embodiments alone or in various combinations. Therefore, innumerable variations not exemplified are envisioned within the scope of the techniques disclosed herein. For example, it is assumed that at least one component is modified, added or omitted, and further, at least one component is extracted and combined with the components of other embodiments.

本願は、放射性物質の三次元濃度分布を測定するガンマカメラとして利用することができる。 The present application can be used as a gamma camera for measuring the three-dimensional concentration distribution of radioactive substances.

1 放射線検出器、2 ガンマ線データ記憶部、3 物性情報取得装置、4 物性情報記憶部、5 測定条件取得装置、6 測定条件記憶部、7、7A 三次元分布演算装置、8 表示部、10 測定対象、11 放射性物質、12 直接成分、13 散乱成分、14 形状情報取得装置、15 形状情報記憶部、21 検出エネルギー記憶部、22 入射方向記憶部、70 処理回路、71 プロセッサ、72 メモリ、100、100A ガンマカメラ 1 radiation detector, 2 gamma ray data storage unit, 3 physical property information acquisition device, 4 physical property information storage unit, 5 measurement condition acquisition device, 6 measurement condition storage unit, 7, 7A three-dimensional distribution calculation device, 8 display unit, 10 measurement Target, 11 radioactive material, 12 direct component, 13 scattered component, 14 shape information acquisition device, 15 shape information storage unit, 21 detection energy storage unit, 22 incident direction storage unit, 70 processing circuit, 71 processor, 72 memory, 100, 100A gamma camera

Claims (11)

測定対象から飛来するガンマ線の入射方向とエネルギーを検出する放射線検出器と、
前記放射線検出器により取得された入射方向情報とエネルギー情報を記憶するガンマ線データ記憶部と、
少なくとも測定対象の成分及び密度を含む物性情報を取得する物性情報取得装置と、
前記物性情報取得装置により取得された測定対象の物性情報を記憶する物性情報記憶部と、
前記放射線検出器の測定条件情報を取得する測定条件取得装置と、
前記測定条件取得装置により取得された測定条件情報を記憶する測定条件記憶部と、
前記ガンマ線データ記憶部に記憶された入射方向情報とエネルギー情報、前記物性情報記憶部に記憶された物性情報、及び前記測定条件記憶部に記憶された測定条件情報に基づいて、測定対象内部で散乱されたガンマ線のエネルギー変化と方向変化を計算し、放射性物質の三次元濃度分布の演算を行う三次元分布演算装置を備えたことを特徴とするガンマカメラ。
A radiation detector that detects the incident direction and energy of gamma rays coming from the measurement target,
A gamma-ray data storage unit that stores incident direction information and energy information acquired by the radiation detector, and
A physical property information acquisition device that acquires physical property information including at least the components and densities to be measured,
A physical property information storage unit that stores physical property information of a measurement target acquired by the physical property information acquisition device, and a physical property information storage unit.
A measurement condition acquisition device that acquires measurement condition information of the radiation detector,
A measurement condition storage unit that stores measurement condition information acquired by the measurement condition acquisition device, and a measurement condition storage unit.
Scattered inside the measurement target based on the incident direction information and energy information stored in the gamma ray data storage unit, the physical property information stored in the physical property information storage unit, and the measurement condition information stored in the measurement condition storage unit. A gamma camera characterized by being equipped with a three-dimensional distribution calculation device that calculates the energy change and direction change of the gamma rays and calculates the three-dimensional concentration distribution of radioactive substances.
測定対象の形状情報を取得する形状情報取得装置と、前記形状情報取得装置により取得された形状情報を記憶する形状情報記憶部とを備えたことを特徴とする請求項1記載のガンマカメラ。 The gamma camera according to claim 1, further comprising a shape information acquisition device for acquiring shape information of a measurement target and a shape information storage unit for storing the shape information acquired by the shape information acquisition device. 前記三次元分布演算装置は、前記形状情報記憶部から取得した測定対象の形状情報に基づいて、放射性物質の分布領域を補正することを特徴とする請求項2記載のガンマカメラ。 The gamma camera according to claim 2, wherein the three-dimensional distribution calculation device corrects a distribution region of a radioactive substance based on the shape information of a measurement target acquired from the shape information storage unit. 前記物性情報取得装置は、測定対象の物性情報として、測定対象の密度分布、構成元素の混合比、及び厚さ寸法を取得することを特徴とする請求項1から請求項3のいずれか一項に記載のガンマカメラ。 One of claims 1 to 3, wherein the physical property information acquisition device acquires the density distribution of the measurement target, the mixing ratio of the constituent elements, and the thickness dimension as the physical property information of the measurement target. Gamma camera described in. 前記物性情報取得装置は、エックス線検査装置、ガンマ線検査装置、中性子線検査装置、超音波検査装置、及び音波検査装置から選ばれた複数の装置を含むことを特徴とする請求項1から請求項4のいずれか一項に記載のガンマカメラ。 Claims 1 to 4 are characterized in that the physical property information acquisition device includes a plurality of devices selected from an X-ray inspection device, a gamma ray inspection device, a neutron beam inspection device, an ultrasonic inspection device, and a sound wave inspection device. The gamma camera described in any one of the items. 前記測定条件取得装置は、測定条件情報として、前記放射線検出器と測定対象の距離と、前記放射線検出器の視野角とを取得することを特徴とする請求項1から請求項5のいずれか一項に記載のガンマカメラ。 One of claims 1 to 5, wherein the measurement condition acquisition device acquires the distance between the radiation detector and the measurement target and the viewing angle of the radiation detector as measurement condition information. The gamma camera described in the section. 前記測定条件取得装置は、レーザ距離計、電波距離計、超音波距離計、及び巻尺のいずれかを含むことを特徴とする請求項1から請求項6のいずれか一項に記載のガンマカメラ。 The gamma camera according to any one of claims 1 to 6, wherein the measurement condition acquisition device includes any one of a laser range finder, a radio range finder, an ultrasonic range finder, and a tape measure. 前記三次元分布演算装置は、測定対象の物性情報に基づいて推定した放射性核種から放出される放射線エネルギーを算出し、算出した放射線エネルギーと前記ガンマ線データ記憶部から取得したエネルギー情報とに基づいてガンマ線のエネルギー変化を計算することを特徴とする請求項1から請求項7のいずれか一項に記載のガンマカメラ。 The three-dimensional distribution calculation device calculates the radiation energy emitted from the radioactive nuclei estimated based on the physical property information of the measurement target, and gamma rays are based on the calculated radiation energy and the energy information acquired from the gamma ray data storage unit. The gamma camera according to any one of claims 1 to 7, wherein the energy change of the above is calculated. 前記三次元分布演算装置は、測定対象内部で散乱されたガンマ線のエネルギー変化からガンマ線が測定対象内部で散乱された角度を算出し、算出した角度と前記ガンマ線データ記憶部から取得した入射方向情報とに基づいて放射性物質の分布領域を限定することを特徴とする請求項1から請求項8記載のいずれか一項に記載のガンマカメラ。 The three-dimensional distribution calculation device calculates the angle at which gamma rays are scattered inside the measurement target from the energy change of the gamma rays scattered inside the measurement target, and the calculated angle and the incident direction information acquired from the gamma ray data storage unit are used. The gamma camera according to any one of claims 1 to 8, wherein the distribution region of the radioactive substance is limited based on the above. 前記三次元分布演算装置は、測定対象の物性情報に基づいて測定対象内部におけるガンマ線の散乱位置を推定し、推定した散乱位置による重み付けにより放射性物質の分布領域を限定することを特徴とする請求項1から請求項9のいずれか一項に記載のガンマカメラ。 The claim is characterized in that the three-dimensional distribution calculation device estimates the scattering position of gamma rays inside the measurement target based on the physical property information of the measurement target, and limits the distribution region of the radioactive substance by weighting by the estimated scattering position. The gamma camera according to any one of claims 1 to 9. 前記放射線検出器が複数のガンマ線の入射方向とエネルギーを検出した場合、前記三次元分布演算装置は、複数のガンマ線それぞれの入射方向情報とエネルギー情報に基づいて放射性物質の分布領域を推定し、推定した複数の分布領域が重なり合った部分を放射性物質の分布領域とすることを特徴とする請求項1から請求項10のいずれか一項に記載のガンマカメラ。 When the radiation detector detects the incident direction and energy of a plurality of gamma rays, the three-dimensional distribution calculation device estimates and estimates the distribution region of the radioactive substance based on the incident direction information and energy information of each of the plurality of gamma rays. The gamma camera according to any one of claims 1 to 10, wherein the portion where the plurality of distribution regions overlap each other is used as the distribution region of the radioactive substance.
JP2018153467A 2018-08-17 2018-08-17 Gamma camera Active JP6983129B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2018153467A JP6983129B2 (en) 2018-08-17 2018-08-17 Gamma camera

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2018153467A JP6983129B2 (en) 2018-08-17 2018-08-17 Gamma camera

Publications (2)

Publication Number Publication Date
JP2020027081A JP2020027081A (en) 2020-02-20
JP6983129B2 true JP6983129B2 (en) 2021-12-17

Family

ID=69619958

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018153467A Active JP6983129B2 (en) 2018-08-17 2018-08-17 Gamma camera

Country Status (1)

Country Link
JP (1) JP6983129B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116363024A (en) * 2021-12-28 2023-06-30 同方威视技术股份有限公司 Gamma camera imaging method and device

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8461547B2 (en) * 2005-04-28 2013-06-11 Utah State University Suppressed correlation method for identifying radioactive sources
JP5955760B2 (en) * 2012-12-20 2016-07-20 日立Geニュークリア・エナジー株式会社 Radiation measurement equipment
JP2016138753A (en) * 2015-01-26 2016-08-04 学校法人早稲田大学 Method and device for measuring three-dimensional distribution of radioactive substance

Also Published As

Publication number Publication date
JP2020027081A (en) 2020-02-20

Similar Documents

Publication Publication Date Title
US9867580B2 (en) X-ray imaging based on image data from a photon-counting multi-bin X-ray detector
US8188437B2 (en) Radiation detecting method utilizing energy information and positional information and equipment thereof
US8354648B2 (en) Device for locating and imaging gamma or X-radiation sources
JP4836094B2 (en) Nuclear medicine imaging apparatus and initial scattering position determination method thereof
US20110147594A1 (en) Radiation imaging method with individual signal resolution
US20110309252A1 (en) Nuclear medicine imaging apparatus, and nuclear medicine imaging method
EP2702449B1 (en) System and method for correction of geometric distortion of multi-camera flat panel x-ray detectors
Turecek et al. Single layer Compton camera based on Timepix3 technology
US20120104262A1 (en) Evaluation of measurements from a pixelated detector
Turecek et al. Compton camera based on Timepix3 technology
González-Montoro et al. Novel method to measure the intrinsic spatial resolution in PET detectors based on monolithic crystals
US20160370474A1 (en) System and method to improve image quality of emission tomography when using advanced radionuclides
US11703605B2 (en) System and method for imaging by gamma radiation detection
Vaquero et al. Performance characteristics of a compact position-sensitive LSO detector module
JP7317586B2 (en) MEDICAL IMAGE PROCESSING APPARATUS, METHOD AND PROGRAM
Shimazoe et al. Double photon emission coincidence imaging using GAGG-SiPM pixel detectors
JP2535762B2 (en) Simultaneous Scattering Counting Method with Gamma Absorber in Positron Tomography Equipment and Positron Tomography Equipment
Chen et al. Calibration of scintillation cameras and pinhole SPECT imaging systems
JP6983129B2 (en) Gamma camera
Ueno et al. Development of a high sensitivity pinhole type gamma camera using semiconductors for low dose rate fields
JPH1172564A (en) Gamma camera system
Iltis et al. Temporal imaging CeBr3 Compton camera: A new concept for nuclear decommissioning and nuclear waste management
JP4984963B2 (en) Nuclear medicine diagnostic equipment
US7569827B2 (en) Emission-data-based photon scatter correction in computed nuclear imaging technology
JP2007510150A (en) Calibration method and apparatus for pixelated solid state detector

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20201127

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20201127

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20211019

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20211026

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20211122

R151 Written notification of patent or utility model registration

Ref document number: 6983129

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151