JP2020027081A - Gamma camera - Google Patents

Gamma camera Download PDF

Info

Publication number
JP2020027081A
JP2020027081A JP2018153467A JP2018153467A JP2020027081A JP 2020027081 A JP2020027081 A JP 2020027081A JP 2018153467 A JP2018153467 A JP 2018153467A JP 2018153467 A JP2018153467 A JP 2018153467A JP 2020027081 A JP2020027081 A JP 2020027081A
Authority
JP
Japan
Prior art keywords
information
gamma
measurement target
energy
measurement
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.)
Granted
Application number
JP2018153467A
Other languages
Japanese (ja)
Other versions
JP6983129B2 (en
Inventor
西沢 博志
Hiroshi Nishizawa
博志 西沢
中西 正一
Shoichi Nakanishi
正一 中西
理 笹野
Makoto Sasano
理 笹野
真照 林
Masateru Hayashi
真照 林
哲史 東
Tetsushi Azuma
哲史 東
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

Abstract

To achieve gamma cameras that can highly sensitively measure a three-dimensional concentration distribution of a radioactive substance in a short time even when a measurement object is a high density substance.SOLUTION: A gamma camera 100 is configured to: calculate an energy change and direction change of a gamma ray scattered inside a measurement object on the basis of incident direction information and energy information on a gamma ray acquired by a radiation detector, physical information on a measurement object acquired by a physical information acquisition device, and measurement condition information on the radiation detector acquired by a measurement condition acquisition device; restrict a distribution area of a radioactive substance; and further restrict the distribution area of the radioactive substance due to weighting by a scattering probability and absorption probability calculated based on the physical information of the measurement object. Accordingly, even when the measurement object is a high density substance, the gamma camera can highly sensitively measure a three-dimensional concentration distribution of the radioactive substance in a short time.SELECTED DRAWING: Figure 7

Description

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

従来、ガンマカメラは、放射性物質の濃度分布の測定に用いられるが、奥行き方向の分布を測定することが難しく、1台のガンマカメラで三次元濃度分布を測定する方法が検討されている。例えば特許文献1では、1台のガンマカメラを用い、測定対象の正面に対する垂直線から左右にそれぞれ同一の角度傾いた撮像軸上で2つの画像を撮像し、得られた2つの画像から立体画像を得ている。   Conventionally, a gamma camera is used for measuring the concentration distribution of a radioactive substance, 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 Literature 1, a single gamma camera is used to capture two images on imaging axes inclined at the same angle to the left and right from a vertical line with respect to the front of the measurement target, and a stereoscopic image is obtained from the obtained two images. Have gained.

また、ガンマ線源の三次元分布を画像化する方法として、コンプトンカメラを利用する方法が知られている。典型的なコンプトンカメラは、2つの放射線検出器を備えており、ガンマ線源からのガンマ線が前段検出器にてコンプトン散乱された後、後段検出器にて光電吸収される事象に注目し、各検出器におけるガンマ線の検出位置及び検出エネルギーに基づいてガンマ線源の三次元分布を画像化する。   As a method for imaging the 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, and focuses on events in which gamma rays from a gamma ray source are Compton-scattered by the first detector and then photoelectrically absorbed by the second detector. The three-dimensional distribution of the gamma ray source is imaged based on the detection position and the detected energy of the gamma ray in the detector.

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

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

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.W. H. Clinthorne, J .; A. Fessler, and W.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.

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

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

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

本願に開示されるガンマカメラは、測定対象から飛来するガンマ線の入射方向とエネルギーを検出する放射線検出器と、放射線検出器により取得された入射方向情報とエネルギー情報を記憶するガンマ線データ記憶部と、少なくとも測定対象の成分及び密度を含む物性情報を取得する物性情報取得装置と、物性情報取得装置により取得された測定対象の物性情報を記憶する物性情報記憶部と、放射線検出器の測定条件情報を取得する測定条件取得装置と、測定条件取得装置により取得された測定条件情報を記憶する測定条件記憶部と、ガンマ線データ記憶部に記憶された入射方向情報とエネルギー情報、物性情報記憶部に記憶された物性情報、及び測定条件記憶部に記憶された測定条件情報に基づいて、測定対象内部で散乱されたガンマ線のエネルギー変化と方向変化を計算し、放射性物質の三次元濃度分布の演算を行う三次元分布演算装置を備えたものである。   The gamma camera disclosed in the present application, 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 the incident direction information and energy information acquired by the radiation detector, Physical property information acquisition device that acquires physical property information including at least the component and density 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 the energy information stored in the gamma ray data storage unit, and the physical property information storage unit Based on the physical property information and the measurement condition information stored in the measurement condition storage unit. Calculate the energy change and direction change, those having a three-dimensional distribution calculating device for performing an operation of a three-dimensional concentration distribution of radioactive material.

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

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

実施の形態1.
以下に、実施の形態1によるガンマカメラについて、図面に基づいて説明する。図1は、実施の形態1に係るガンマカメラの構成図、図2は、実施の形態1に係るガンマカメラによる放射性物質の三次元濃度分布の測定方法を説明する図、図3は、実施の形態1に係るガンマカメラにおける三次元分布演算装置のハードウェア構成図である。なお、各図において、同一、相当部分には同一符号を付している。
Embodiment 1 FIG.
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 Embodiment 1, 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 Embodiment 1, and FIG. FIG. 3 is a hardware configuration diagram of a three-dimensional distribution calculation device in the gamma camera according to Embodiment 1. In each of the drawings, the same or corresponding parts have the same reference characters allotted.

実施の形態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 a human body, such as concrete, metal, and the like. The radioactive material present in the measurement object contains radionuclides that emit radiation and decay to be converted into other nuclei. Gamma rays (hereinafter referred to as γ-rays) are emitted from γ-ray nuclides. As shown in FIG. 2, γ-rays coming from a radioactive substance 11 present inside the measurement target 10 are divided into a direct component 12 directly incident on the gamma camera 100 and a 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 a measurement condition acquisition. The apparatus 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 γ-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 scatter component 13.

放射線検出器1は、γ線の入射によりγ線が反応した位置と、エネルギーに応じて決まる電圧もしくは電流強度を出力するγ線センサを有している。γ線センサとしては、半導体式検出器を複数配置したもの、半導体式検出器において信号を取り出す電極を複数有するもの、γ線を紫外光、可視光、及び赤外光に変換するシンチレータと光検出部を含むユニットを複数配置したもの等が用いられる。   The radiation detector 1 has a γ-ray sensor that outputs a voltage or current intensity determined according to the position where the γ-ray has reacted due to the incidence of the γ-ray. A gamma ray sensor includes a plurality of semiconductor detectors, a plurality of electrodes for extracting signals from the semiconductor detector, a scintillator for converting gamma rays into ultraviolet light, visible light, and infrared light, and light detection. 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 detector, a detector using Si, CdTe, CdZnTe, or the like is used. In the case of a semiconductor detector, the position where the γ-ray has reacted is detected from the position where the semiconductor detector to which the γ-ray has reacted or the position of the electrode from which the signal has been extracted. On the other hand, the scintillator includes a scintillation material such as NaI (Tl), CsI (Tl), BGO, LSO, and the like, and the light from the scintillator is detected by a light detection unit. In the case of a scintillator, the position where the γ-ray has reacted is detected from the position where 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 combining a γ-ray sensor and a housing having one or more holes in a shielding body that blocks γ-rays, or a configuration combining a γ-ray sensor and a collimator that limits the incident direction of γ-rays Or a configuration in which two or more γ-ray sensors are arranged. The γ-ray data storage unit 2 includes a detected energy storage unit 21 and an incident direction storage unit 22, and stores energy information and incident direction information of γ-rays detected by the radiation detector 1, respectively.

物性情報取得装置3は、測定対象の物性情報として、測定対象の密度、密度分布、成分(構成元素)、構成元素の混合比、及び厚さ寸法等を測定する。物性情報取得装置3により取得される物性情報には、少なくとも測定対象の成分及び密度が含まれる。物性情報取得装置3は、エックス線検査装置、γ線検査装置、中性子線検査装置、超音波検査装置、及び音波検査装置の中から選ばれた複数の装置を含む。なお、上記以外の装置を用いることも可能である。   The physical property information acquisition device 3 measures the density, the density distribution, the components (constituent elements), the mixing ratio of the constituent elements, the thickness, and the like of the measurement target as the physical property information of the measurement target. The physical property information acquired by the physical property information acquiring device 3 includes at least a component to be measured and a density. The physical property information acquiring device 3 includes a plurality of devices selected from an X-ray inspection device, a γ-ray inspection device, a neutron inspection device, an ultrasonic inspection device, and an ultrasonic inspection device. In addition, it is also possible to use apparatuses 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 object is measured by an ultrasonic inspection device or an ultrasonic inspection device. The physical property information storage unit 4 stores the physical property information acquired by the physical property information acquiring device 3. The measurement of the physical property information by the physical property information acquisition device 3 does not need to be performed at the same time as the measurement of the three-dimensional concentration distribution of the radioactive substance, and a part of the measurement target is collected in advance and the physical property information is measured. It is also possible to record it in the storage unit 4.

測定条件取得装置5は、放射線検出器1の測定条件情報を取得する。測定条件情報には、放射線検出器1と測定対象の距離、及び放射線検出器1の視野角が含まれる。測定条件取得装置5は、レーザ距離計、電波距離計、超音波距離計、及び巻尺のいずれかを含む。ただし、上記以外の装置を用いることも可能である。測定条件記憶部6は、測定条件取得装置5により取得された測定条件情報を記憶する。   The measurement condition acquisition device 5 acquires 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 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 them, the energy change and direction change of the γ-ray scattered inside the measurement object are calculated, and the three-dimensional concentration distribution of the radioactive material is calculated. Specifically, the radiation energy emitted from the radionuclide estimated based on the physical property information of the measurement target is calculated, and the γ-ray is calculated 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 θ at 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 calculated. Based on this, the distribution area of the radioactive material is limited. Further, the scattering position of the γ-ray 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 based on the estimated scattering position. As one of the calculation methods at this time, it is possible to apply the ML-EM method, the maximum entropy method, or the like that estimates the three-dimensional concentration distribution of the radioactive substance so that the distribution becomes the most reliable.

三次元分布演算装置7による上記の演算処理は、放射線検出器1のγ線センサが検出した一つのγ線に対して行われる。放射線検出器1が複数のγ線を検出した場合、三次元分布演算装置7は、複数のγ線それぞれの入射方向情報とエネルギー情報に基づいて放射性物質の分布領域を推定し、さらに推定した複数の分布領域が重なり合った部分を放射性物質の分布領域とする。三次元分布演算装置7による演算結果は、表示部8に表示される。   The above calculation processing by the three-dimensional distribution calculation device 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 area 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. A portion where the distribution regions overlap each other is defined as a radioactive material distribution region. 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にプログラムが入力される。   Note that the three-dimensional distribution calculation device 7 has a processing circuit 70 including a processor 71 and a 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, an auxiliary storage device of a hard disk may be provided instead of the flash memory. The processor 71 executes a 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, a calculation method by the three-dimensional distribution calculation device 7 will be described in detail. Γ-rays emitted from radioactive materials undergo energy changes and direction changes due to scattering. The three-dimensional distribution calculation device 7 calculates an energy change and a direction change due to γ-ray scattering, and limits the distribution region of the radioactive substance based on the calculation result.

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

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

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

測定対象が複数の元素の集合体で構成される場合には、それぞれの元素の集合体での密度ρ、γ線がそれぞれの元素の集合体を進んだ距離L、及び散乱断面積σによる下の式2によって、測定対象内部の箇所毎に散乱される確率pが算出される。 When the measurement object is composed of an aggregate of a plurality of elements, the density ρ i of the aggregate of each element, the distance L i that the γ-ray has traveled through the aggregate of each element, and the scattering cross section σ The probability p that is scattered for each location inside the measurement target is calculated by the following Expression 2 based on i .

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

なお、γ線が散乱される確率pの算出方法は、式1及び式2を用いた上記の方法に限定されるものではない。別の算出方法として、物性情報取得装置3で取得した測定対象の物性情報を用い、乱数計算によって測定対象の物性情報を再現した条件でシミュレーションを行い、γ線が散乱される確率pを算出する方法がある。   Note that the method of calculating the probability p that γ rays are scattered is not limited to the above method using Expressions 1 and 2. As another calculation method, the physical property information of the measurement target acquired by the physical property information acquisition device 3 is used, a simulation is performed under conditions in which 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, it is possible to calculate the probability that γ-rays are absorbed at each location inside the measurement target based on the physical property information obtained from the physical property information storage unit 4 by the above method using Equations 1 and 2. It is. In that case, the cross-sectional area where the element constituting the measurement object absorbs γ-rays is applied to Equations 1 and 2. Similarly, the simulation of calculating the probability p that γ rays are scattered can also calculate the probability that γ rays are absorbed.

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

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

測定対象内部のある箇所において、その位置でγ線が散乱される確率は、散乱確率と吸収確率の差(図4中、Dで示す)のような確率分布となる。従って、測定対象内部でγ線が散乱された箇所の位置情報は、散乱確率と吸収確率の差の確率分布を用いて推定することができる。測定対象が複数の元素の集合体で構成される場合には、各々の元素に対応する散乱確率と吸収確率の差の確率分布が用いられる。散乱確率と吸収確率の差の確率分布により推定された位置情報は、放射性物質の分布領域を限定するための重み付けに用いられる。   At a certain location inside the measurement target, the probability that the γ-ray is scattered at that location has a probability distribution like 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 using the probability distribution of the difference between the scattering probability and the absorption probability. When the measurement object 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 for limiting the distribution region of the radioactive material.

また、測定対象内部でγ線が散乱された場合には、コンプトン散乱によるγ線のエネルギー損失が起き、エネルギー損失後のγ線のエネルギーE´は、散乱される前のエネルギーE、電子の静止エネルギーE、及び散乱された角度θにより、下の式3で示される。式3において、Eは散乱される前のγ線のエネルギー、E´はコンプトン散乱によるエネルギー損失後のγ線のエネルギー、θは放射性物質から放出されたγ線が測定対象内部で散乱された角度である。 When γ rays are scattered inside the measurement object, energy loss of γ rays occurs due to Compton scattering, and the energy E ′ of the γ rays after the energy loss becomes the energy E before being scattered and the static energy of electrons. The energy E 0 and the scattered angle θ are represented by the following Equation 3. In Equation 3, E is the energy of the γ-ray before being scattered, E ′ is the energy of the γ-ray after the energy loss due to Compton scattering, and θ is the angle at which the γ-ray emitted from the radioactive material is scattered inside the measurement object. It 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. The energy E ′ of the γ-ray after the energy loss is the energy of the γ-ray detected by the radiation detector 1 and is stored in the detected energy storage unit 21 as energy information. The angle θ at which the γ-rays emitted from the radioactive substance are scattered inside the measurement target is obtained from the equation (3).

実施の形態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 Embodiment 1 will be described with reference to the flowchart in FIG. First, in step S <b> 1, energy information of γ-rays detected by the radiation detector 1 is obtained from the detected energy storage unit 21. Subsequently, in step S2, a 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 object using the γ-ray energy information obtained 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-described 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 operation for limiting the distribution region of the radioactive material is performed using the angle θ calculated in step S4 and the γ-ray incident direction information obtained in step S5. The distribution region of the radioactive substance obtained from the angle θ calculated by Expression 3 and the incident direction information of the γ-ray 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.

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

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

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

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

三次元分布演算装置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 first γ-ray energy change and direction change (step S6 in FIG. 6), and further weights the scattering probability and the absorption probability. (Step S9 in FIG. 6) to estimate the distribution region A1 of the radioactive material. Similarly, the distribution area A2 is estimated by the second γ-ray, and the distribution area A3 of the radioactive substance is estimated by superimposing these distribution areas A1 and A2. In the region where the radioactive substance is truly present, the overlapping of the distribution regions of the radioactive substance calculated from the respective γ-rays increases. In other words, the greater the overlap, the higher the probability of being a true radioactive substance distribution area, and the more accurate the measurement of the three-dimensional concentration distribution.

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

また、測定対象の物性情報に基づいて算出した散乱確率と吸収確率による重み付けにより、放射性物質の分布領域をさらに限定するようにしたので、測定精度が高くなる。また、複数のγ線を検出した場合には、それぞれのγ線から推定した放射性物質の分布領域を重ね合わせることにより、さらに測定精度が向上する。   In addition, the distribution area of the radioactive substance is further limited by weighting based on the scattering probability and the absorption probability calculated based on the physical property information of the measurement object, so that the measurement accuracy is improved. When a plurality of γ-rays are detected, the measurement accuracy is further improved by overlapping 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.
FIG. 8 shows a configuration diagram of the gamma camera according to the second embodiment. The gamma camera 100A according to the second embodiment stores, in the gamma camera 100 according to the first embodiment, a shape information acquisition device 14 that acquires a shape of a measurement target and shape information acquired by the shape information acquisition device 14. And a shape information storage unit 15 to be added. As the shape information acquisition device 14, a three-dimensional shape measurement device, a contour shape measurement device, or the like is used. The other configuration is the same as that of the first embodiment, and the description is omitted here.

上記実施の形態1に係るガンマカメラ100は、測定対象以外の領域からのγ線を含めて演算を行っている。これに対し、実施の形態2に係るガンマカメラ100Aの三次元分布演算装置7Aは、形状情報記憶部15から取得した測定対象の形状情報に基づいて、放射性物質の分布領域を補正することができる。具体的には、演算により得られた放射性物質の分布領域と、形状情報取得装置14により取得した測定対象の形状情報とを比較し、測定対象以外の領域について演算した結果を取り除く。   The gamma camera 100 according to Embodiment 1 performs the 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 material 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 result of the calculation 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 Embodiment 2 will be described with reference to the flowchart in FIG. Note that the processing from step S11 to step S18 in FIG. 9 is the same as the processing from step S1 to step S8 in the flowchart in FIG. 6 described in the first embodiment, and a description thereof will not be repeated. In step S19 of FIG. 9, the γ in the measurement target is determined based on the incident direction information of the γ-ray acquired in step S15, the scattering probability of the γ-ray calculated in step S17, and the absorption probability of the γ-ray 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 area of the radioactive material defined in step S16 is weighted by the γ-ray scattering position estimated in step S19 in the same manner as in the first embodiment, and the distribution area of the radioactive substance is determined. Is further limited. Next, in step S21, the shape information of the measurement target stored in the shape information storage unit 15 is obtained. 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 defined in step S20, and the distribution region is corrected.

なお、ステップS11からステップS22までの処理の順番は、適宜入れ替えが可能である。ステップS11からステップS15の処理は、ステップS16の処理の前に行われていれば良く、ステップS11の処理をステップS12とステップS13の処理の後に行っても良い。また、ステップS21の形状情報取得処理をステップS16の処理の前に行い、測定対象以外の領域での放射性物質の三次元濃度分布の演算を行わないようにしてもよい。   Note that the order of the processes 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 processing in step S21 may be performed before the processing in step S16, and 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 effects as in the first embodiment, the distribution region of the radioactive substance obtained by the calculation is compared with the shape information of the measurement object acquired by the shape information acquiring device 14 to obtain the distribution. Since the correction of the region is performed, the calculation result for the region 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 from being performed, the time required for the calculation can be reduced as compared with the first embodiment.

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

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

1 放射線検出器、2 ガンマ線データ記憶部、3 物性情報取得装置、4 物性情報記憶部、5 測定条件取得装置、6 測定条件記憶部、7、7A 三次元分布演算装置、8 表示部、10 測定対象、11 放射性物質、12 直接成分、13 散乱成分、14 形状情報取得装置、15 形状情報記憶部、21 検出エネルギー記憶部、22 入射方向記憶部、70 処理回路、71 プロセッサ、72 メモリ、100、100A ガンマカメラ   Reference Signs List 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 Object, 11 radioactive material, 12 direct component, 13 scattering component, 14 shape information acquisition device, 15 shape information storage unit, 21 detected energy storage unit, 22 incident direction storage unit, 70 processing circuit, 71 processor, 72 memory, 100, 100A gamma camera

Claims (11)

測定対象から飛来するガンマ線の入射方向とエネルギーを検出する放射線検出器と、
前記放射線検出器により取得された入射方向情報とエネルギー情報を記憶するガンマ線データ記憶部と、
少なくとも測定対象の成分及び密度を含む物性情報を取得する物性情報取得装置と、
前記物性情報取得装置により取得された測定対象の物性情報を記憶する物性情報記憶部と、
前記放射線検出器の測定条件情報を取得する測定条件取得装置と、
前記測定条件取得装置により取得された測定条件情報を記憶する測定条件記憶部と、
前記ガンマ線データ記憶部に記憶された入射方向情報とエネルギー情報、前記物性情報記憶部に記憶された物性情報、及び前記測定条件記憶部に記憶された測定条件情報に基づいて、測定対象内部で散乱されたガンマ線のエネルギー変化と方向変化を計算し、放射性物質の三次元濃度分布の演算を行う三次元分布演算装置を備えたことを特徴とするガンマカメラ。
A radiation detector for detecting the incident direction and energy of gamma rays coming from the measurement object,
A gamma ray data storage unit for storing incident direction information and energy information obtained by the radiation detector,
Physical property information acquisition device that acquires physical property information including at least the component and density of the measurement target,
Physical property information storage unit that stores the physical property information of the measurement target acquired by the physical property information acquisition device,
A measurement condition acquisition device for acquiring measurement condition information of the radiation detector,
A measurement condition storage unit that stores measurement condition information acquired by the measurement condition acquisition device,
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 scattering inside the measurement target based on the measurement condition information stored in the measurement condition storage unit. A gamma camera comprising a three-dimensional distribution calculation device for calculating a change in energy and a change in direction of the gamma ray and calculating a three-dimensional concentration distribution of the radioactive material.
測定対象の形状情報を取得する形状情報取得装置と、前記形状情報取得装置により取得された形状情報を記憶する形状情報記憶部とを備えたことを特徴とする請求項1記載のガンマカメラ。   The gamma camera according to claim 1, further comprising: a shape information acquisition device that acquires shape information of a measurement target; and a shape information storage unit that stores 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 the radioactive material based on shape information of the measurement target acquired from the shape information storage unit. 前記物性情報取得装置は、測定対象の物性情報として、測定対象の密度分布、構成元素の混合比、及び厚さ寸法を取得することを特徴とする請求項1から請求項3のいずれか一項に記載のガンマカメラ。   4. The physical property information acquisition device according to claim 1, wherein the physical property information of the measurement target acquires a density distribution, a mixing ratio of constituent elements, and a thickness dimension of the measurement target. The gamma camera according to the above. 前記物性情報取得装置は、エックス線検査装置、ガンマ線検査装置、中性子線検査装置、超音波検査装置、及び音波検査装置から選ばれた複数の装置を含むことを特徴とする請求項1から請求項4のいずれか一項に記載のガンマカメラ。   The said physical-property information acquisition apparatus contains several apparatuses selected from the X-ray inspection apparatus, the gamma-ray inspection apparatus, the neutron beam inspection apparatus, the ultrasonic inspection apparatus, and the sonic inspection apparatus. The gamma camera according to any one of the above. 前記測定条件取得装置は、測定条件情報として、前記放射線検出器と測定対象の距離と、前記放射線検出器の視野角とを取得することを特徴とする請求項1から請求項5のいずれか一項に記載のガンマカメラ。   6. The measurement condition acquisition apparatus according to claim 1, wherein the measurement condition information acquires a distance between the radiation detector and a measurement target and a viewing angle of the radiation detector as measurement condition information. The gamma camera according to the item. 前記測定条件取得装置は、レーザ距離計、電波距離計、超音波距離計、及び巻尺のいずれかを含むことを特徴とする請求項1から請求項6のいずれか一項に記載のガンマカメラ。   The gamma camera according to any one of claims 1 to 6, wherein the measurement condition acquisition device includes 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 radionuclide estimated based on the physical property information of the measurement target, and calculates the gamma ray based on the calculated radiation energy and the energy information obtained from the gamma ray data storage unit. The gamma camera according to any one of claims 1 to 7, wherein an energy change of the gamma camera 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 calculates the calculated angle and the incident direction information obtained from the gamma ray data storage unit. The gamma camera according to any one of claims 1 to 8, wherein a distribution region of the radioactive substance is limited on the basis of: 前記三次元分布演算装置は、測定対象の物性情報に基づいて測定対象内部におけるガンマ線の散乱位置を推定し、推定した散乱位置による重み付けにより放射性物質の分布領域を限定することを特徴とする請求項1から請求項9のいずれか一項に記載のガンマカメラ。   The three-dimensional distribution calculation device estimates a gamma ray scattering position inside the measurement target based on physical property information of the measurement target, and limits a distribution region of the radioactive material by weighting with 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 the distribution area of the radioactive material based on the incident direction information and the energy information of each of the plurality of gamma rays, and estimates The gamma camera according to any one of claims 1 to 10, wherein a portion where the plurality of distribution regions overlap each other is defined as a distribution region of a 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 true JP2020027081A (en) 2020-02-20
JP6983129B2 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)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20230204798A1 (en) * 2021-12-28 2023-06-29 Nuctech Company Limited Gamma camera imaging method and gamma camera imaging device

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120053847A1 (en) * 2005-04-28 2012-03-01 Utah State University Suppressed correlation method for identifying radioactive sources
JP2014122793A (en) * 2012-12-20 2014-07-03 Hitachi-Ge Nuclear Energy Ltd Radiation measuring device
JP2016138753A (en) * 2015-01-26 2016-08-04 学校法人早稲田大学 Method and device for measuring three-dimensional distribution of radioactive substance

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120053847A1 (en) * 2005-04-28 2012-03-01 Utah State University Suppressed correlation method for identifying radioactive sources
JP2014122793A (en) * 2012-12-20 2014-07-03 Hitachi-Ge Nuclear Energy Ltd Radiation measuring device
JP2016138753A (en) * 2015-01-26 2016-08-04 学校法人早稲田大学 Method and device for measuring three-dimensional distribution of radioactive substance

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20230204798A1 (en) * 2021-12-28 2023-06-29 Nuctech Company Limited Gamma camera imaging method and gamma camera imaging device

Also Published As

Publication number Publication date
JP6983129B2 (en) 2021-12-17

Similar Documents

Publication Publication Date Title
US9867580B2 (en) X-ray imaging based on image data from a photon-counting multi-bin X-ray detector
US8471210B2 (en) Radiation imaging method with individual signal resolution
US8354648B2 (en) Device for locating and imaging gamma or X-radiation sources
Turecek et al. Compton camera based on Timepix3 technology
US20120104262A1 (en) Evaluation of measurements from a pixelated detector
JP2007212242A (en) Nuclear medicine system
González-Montoro et al. Novel method to measure the intrinsic spatial resolution in PET detectors based on monolithic crystals
JP6990412B2 (en) Gamma ray image acquisition device and gamma ray image acquisition method
JP7317586B2 (en) MEDICAL IMAGE PROCESSING APPARATUS, METHOD AND PROGRAM
Baek et al. Optimization of large-angle pinhole collimator for environmental monitoring system
CN114252900A (en) Counting activity meter for measuring activity of radioactive source
Chen et al. Calibration of scintillation cameras and pinhole SPECT imaging systems
Ueno et al. Development of a high sensitivity pinhole type gamma camera using semiconductors for low dose rate fields
JP4893950B2 (en) Radioactivity absolute measurement method, radiation detector assembly detection efficiency determination method, and radiation measurement apparatus calibration method
JP6983129B2 (en) Gamma camera
JPH07113873A (en) Scattering simultaneous counting measuring method by gamma-ray absorber of positron tomograph and the same tomograph
Huber et al. Development of the LBNL positron emission mammography camera
Peterson et al. Multi-pinhole SPECT imaging with silicon strip detectors
Iltis et al. Temporal imaging CeBr3 Compton camera: A new concept for nuclear decommissioning and nuclear waste management
Ghelman et al. Design of 4π High-Efficiency Directional Radiation Detector Based on Compton Scattering
JP2007510150A (en) Calibration method and apparatus for pixelated solid state detector
US20070200066A1 (en) Emission-data-based photon scatter correction in computed nuclear imaging technology
JP2007218769A (en) Nuclear medicine imaging device
Lee et al. Portable active collimation imager using URA patterned scintillator
Lee et al. A Dual Modality Gamma Camera Using ${\rm LaCl} _ {3}({\rm Ce}) $ Scintillator

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