JP6091852B2 - Nuclear medicine diagnostic apparatus and image processing program - Google Patents
Nuclear medicine diagnostic apparatus and image processing program Download PDFInfo
- Publication number
- JP6091852B2 JP6091852B2 JP2012249536A JP2012249536A JP6091852B2 JP 6091852 B2 JP6091852 B2 JP 6091852B2 JP 2012249536 A JP2012249536 A JP 2012249536A JP 2012249536 A JP2012249536 A JP 2012249536A JP 6091852 B2 JP6091852 B2 JP 6091852B2
- Authority
- JP
- Japan
- Prior art keywords
- data
- detectors
- period
- output data
- heart
- 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
Links
Images
Landscapes
- Nuclear Medicine (AREA)
Description
本発明の実施形態は、核医学診断装置および画像処理プログラムに関する。 Embodiments described herein relate generally to a nuclear medicine diagnostic apparatus and an image processing program.
核医学診断装置は、放射性同位元素(Radio Isotope、以下RIという)を含む薬品(血流マーカ、トレーサ)が生体内の特定組織や臓器に選択的に取り込まれる性質を利用して、生体内に分布したRIから放射されるガンマ線を生体外に配設されたガンマ線の検出器で検出するようになっている。ガンマ線の検出結果は、ガンマ線の線量分布を画像化することによる核医学画像の生成や、体内臓器等の機能の診断などに利用される。 A nuclear medicine diagnostic apparatus uses a property that a medicine (blood flow marker, tracer) containing a radioisotope (hereinafter referred to as RI) is selectively taken into a specific tissue or organ in the living body. Gamma rays emitted from the distributed RI are detected by a gamma ray detector disposed outside the living body. The detection result of gamma rays is used for generation of nuclear medicine images by imaging the dose distribution of gamma rays, diagnosis of functions of internal organs and the like.
胸部や腹部など、肺に近い部位の撮像画像は、呼吸動による影響を受ける。このため、胸部や腹部などを撮影する場合は、呼吸動の影響を低減させることが好ましい。 A captured image of a portion close to the lung, such as the chest and abdomen, is affected by respiratory motion. For this reason, when photographing the chest or abdomen, it is preferable to reduce the influence of respiratory motion.
しかし、呼吸動によって撮像画像が影響を受ける部位には、心臓の拍動によっても影響を受けるものがある。このため、この種の部位を撮影する場合には、呼吸動に加え拍動による画像に対する影響を低減させることが好ましい。 However, some parts of the captured image that are affected by respiratory motion are also affected by the heartbeat. For this reason, when this type of part is imaged, it is preferable to reduce the influence of the pulsation on the image in addition to the respiratory movement.
本発明の一実施形態に係る核医学診断装置は、上述した課題を解決するために、被検体に投与された放射性同位元素から放射されるガンマ線を検出する複数の検出器と、被検体に取り付けられた呼吸センサおよび心電センサの出力データと、複数の検出器の出力データと、を収集するデータ収集部と、データ収集部により収集されたデータにもとづいて、呼吸動の同位相期間であるとともに心臓の動きが閾値より小さい期間である第1の期間における複数の検出器の出力データにもとづく画像を生成する画像生成部と、複数の検出器の出力データにもとづいて、対消滅ガンマ線対の検出時間差から心臓の領域に由来するデータを特定する心臓領域特定部と、を備えたものである。画像生成部は、第1の期間における複数の検出器の出力データと、呼吸動の同位相期間であるとともに心臓の動きが閾値より大きい期間である複数の第2の期間のうち所定数の期間における複数の検出器の出力データから心臓の領域に由来するデータを削除することにより、複数の第2の期間における複数の検出器の出力データから心臓の領域に由来するデータを間引いたデータと、にもとづいて画像を生成する。 In order to solve the above-described problem, a nuclear medicine diagnostic apparatus according to an embodiment of the present invention is provided with a plurality of detectors that detect gamma rays emitted from a radioisotope administered to a subject, and attached to the subject. A data collection unit for collecting the output data of the respiratory sensor and the electrocardiogram sensor and the output data of a plurality of detectors, and the in-phase period of respiratory motion based on the data collected by the data collection unit In addition, an image generation unit that generates images based on the output data of the plurality of detectors in the first period, in which the heart motion is smaller than the threshold , and the pair of annihilation gamma ray pairs based on the output data of the plurality of detectors And a heart region specifying unit that specifies data derived from the heart region from the detection time difference . The image generation unit includes a predetermined number of periods among a plurality of second periods in which the output data of the plurality of detectors in the first period and the in-phase period of the respiratory motion are periods in which the heart motion is greater than the threshold value By deleting data derived from the heart region from the output data of the plurality of detectors in, data obtained by thinning out the data derived from the heart region from the output data of the plurality of detectors in the plurality of second periods, An image is generated based on this.
本発明に係る核医学診断装置および画像処理プログラムの実施の形態について、添付図面を参照して説明する。
(第1の実施形態)
Embodiments of a nuclear medicine diagnosis apparatus and an image processing program according to the present invention will be described with reference to the accompanying drawings.
(First embodiment)
図1は、本発明の第1実施形態に係る核医学診断装置の一例を示すブロック図である。なお、以下の説明では、本発明に係る核医学診断装置としてPET(Positron Emission Tomography)装置を用いる場合の一例について示す。 FIG. 1 is a block diagram showing an example of a nuclear medicine diagnostic apparatus according to the first embodiment of the present invention. In the following description, an example in which a PET (Positron Emission Tomography) apparatus is used as the nuclear medicine diagnostic apparatus according to the present invention will be described.
核医学診断装置10は、スキャナ装置11および画像処理装置12を有する。スキャナ装置11は、天板21、天板駆動装置22、複数の検出器23、検出器カバー24、データ収集部25を有する。
The nuclear medicine
天板21は、患者(被検体)Oを載置可能に構成される。天板駆動装置22は、画像処理装置12に制御されて、天板21を昇降動させる。また、天板駆動装置22は、画像処理装置12に制御されて、検出器カバー24の中央部分の開口部へ天板21を天板21の長軸方向に沿って移送する。
The
検出器23は、FDG(フルオロデオキシグルコース)などの薬品に含まれて患者Oに投与された放射性同位元素から放射されるガンマ線を検出する検出器である。検出器23としては、シンチレータ型検出器を用いてもよいし、半導体型検出器を用いてもよい。
The
シンチレータ型検出器を用いる場合は、検出器23は、ガンマ線の入射角度を規定するためのコリメータ、コリメートされたガンマ線が入射すると瞬間的な閃光を発するシンチレータ、シンチレータから射出された光を検出するための2次元に配列された複数の光電子増倍管、およびシンチレータ用電子回路などを有する。
When a scintillator type detector is used, the
シンチレータは、たとえばタリウム活性化ヨウ化ナトリウムNaI(Tl)により構成される。シンチレータ用電子回路は、ガンマ線が入射する事象(イベント)が発生するごとに、複数の光電子増倍管の出力にもとづいて複数の光電子増倍管により構成される検出面内におけるガンマ線の入射位置情報(位置情報)および強度情報を生成しデータ収集部25に出力する。
The scintillator is made of, for example, thallium activated sodium iodide NaI (Tl). The scintillator electronic circuit has information on the incident position of gamma rays in a detection plane constituted by a plurality of photomultiplier tubes based on the output of the plurality of photomultiplier tubes every time an event in which gamma rays are incident occurs. (Position information) and intensity information are generated and output to the
一方、半導体型検出器を用いる場合は、検出器23は、コリメータ、コリメートされたガンマ線を検出するための2次元に配列された複数のガンマ線検出用半導体素子(以下、半導体素子という)および半導体用電子回路を有する。
On the other hand, when a semiconductor detector is used, the
半導体素子は、たとえばCdTeやCdZnTe(CZT)により構成される。半導体用電子回路は、ガンマ線が入射する事象(イベント)が発生するごとに、半導体素子の出力にもとづいて位置情報および強度情報を生成しデータ収集部25に出力する。
The semiconductor element is made of, for example, CdTe or CdZnTe (CZT). The semiconductor electronic circuit generates position information and intensity information based on the output of the semiconductor element and outputs them to the
複数の検出器23は、たとえば患者Oの周囲を囲むように、六角形または円形に検出器カバー24内に配置される。なお、複数の検出器23の配置態様はリング配列型に限られず、たとえば平板上に配列された複数の検出器23が2つ患者Oを挟んで対向配置されつつ患者Oの周りに回転可能に保持される2検出器対向型に配置されてもよい。また、複数の検出器23は多層のリングに配列されて隣接する層間の画像を取得可能に構成されてもよい。
The plurality of
患者Oには、呼吸センサ26および心電センサ27が取り付けられる。呼吸センサ26が出力する呼吸センサ信号および心電センサ27が出力する心電センサ信号(ECG(Electro CardioGram)信号)は、データ収集部25に与えられる。
A
データ収集部25は、複数の検出器23の出力をリストモードで収集する。リストモードでは、ガンマ線の検出位置情報、強度(エネルギー)情報、検出器23と患者Oとの相対位置を示す情報(検出器23の位置や角度など)、およびガンマ線の検出時刻がガンマ線の入射イベントごとに収集される。また、データ収集部25は、呼吸センサ26および心電センサ27からそれぞれ受けた患者Oの呼吸センサ信号およびECG信号の情報をリストモードデータに関連付けて、データ収集部25に内包するRAMなどの記憶媒体に記憶させる。
The
画像処理装置12は、図1に示すように、制御部31、表示部32、入力部33および記憶部34を有する。
As illustrated in FIG. 1, the
制御部31は、CPU、RAMおよびROMをはじめとする記憶媒体などにより構成され、この記憶媒体に記憶されたプログラムに従って画像処理装置12の処理動作を制御する。
The
制御部31のCPUは、ROMをはじめとする記憶媒体に記憶された画像処理プログラムおよびこのプログラムの実行のために必要なデータをRAMへロードし、このプログラムに従って呼吸動および拍動による撮像画像に対する影響を低減するための処理を実行する。
The CPU of the
制御部31のRAMは、CPUが実行するプログラムおよびデータを一時的に格納するワークエリアを提供する。制御部31のROMをはじめとする記憶媒体は、画像処理装置12の起動プログラム、画像処理プログラムや、これらのプログラムを実行するために必要な各種データを記憶する。
The RAM of the
なお、ROMをはじめとする記憶媒体は、磁気的もしくは光学的記録媒体または半導体メモリなどの、CPUにより読み取り可能な記録媒体を含んだ構成を有し、これら記憶媒体内のプログラムおよびデータの一部または全部は電子ネットワークを介してダウンロードされるように構成してもよい。 A storage medium such as a ROM has a configuration including a recording medium readable by a CPU, such as a magnetic or optical recording medium or a semiconductor memory, and a part of programs and data in the storage medium. Or you may comprise so that all may be downloaded via an electronic network.
表示部32は、たとえば液晶ディスプレイやOLED(Organic Light Emitting Diode)ディスプレイなどの一般的な表示出力装置により構成され、制御部31の制御に従って核医学診断画像などの各種情報を表示する。
The
入力部33は、たとえばキーボード、タッチパネル、テンキーなどの一般的な入力装置により構成され、ユーザの操作に対応した操作入力信号を制御部31に出力する。
The
記憶部34は、磁気的もしくは光学的記録媒体または半導体メモリなどの、CPUにより読み取り可能な記録媒体を含んだ構成を有し、これら記憶媒体内のプログラムおよびデータの一部または全部は電子ネットワークを介してダウンロードされるように構成してもよい。記憶部34は、制御部31により制御されて同期計数情報や核医学画像などを記憶する。
The storage unit 34 has a configuration that includes a CPU-readable recording medium such as a magnetic or optical recording medium or a semiconductor memory, and some or all of the programs and data in the storage medium are stored in an electronic network. You may comprise so that it may be downloaded via. The storage unit 34 is controlled by the
図2は、第1実施形態に係る制御部31のCPUによる機能実現部の構成例を示す概略的なブロック図である。なお、この機能実現部は、CPUを用いることなく回路などのハードウエアロジックによって構成してもよい。
FIG. 2 is a schematic block diagram illustrating a configuration example of the function realization unit by the CPU of the
図2に示すように、制御部31のCPUは、画像処理プログラムによって、少なくともスキャン制御部41、同時計数部42、心肺同期データ抽出部43および画像生成部44として機能する。この各部41〜44は、RAMの所要のワークエリアをデータの一時的な格納場所として利用する。
As illustrated in FIG. 2, the CPU of the
スキャン制御部41は、ユーザから入力部33を介してスキャン計画の実行指示を受けて、スキャン計画にもとづいてスキャナ装置11を制御することにより、スキャンを実行する。この結果、患者Oから放出されたガンマ線の情報、患者Oの呼吸センサ信号およびECG信号の情報がスキャナ装置11からデータ収集部25を介して同時計数部42に与えられる。
The
図3は、同時計数部42により扱われるリストモードデータの一例を示す説明図である。図3には、同時計数した対消滅ガンマ線に関するデータのみを示した。なお、図3において、Detector#1、#2は同時計数した対消滅ガンマ線対が入射した検出器23の位置を示す番号を、Tはデータ収集開始からガンマ線が入射するまでの時間(対消滅ガンマ線2つの検出時間の平均)を、Δtは対消滅ガンマ線の検出時間の差(以下、入射時間差という)を、Energy#1、#2はそれぞれDetector#1、#2で検出されたガンマ線の入射エネルギーを、それぞれ示す。
FIG. 3 is an explanatory diagram showing an example of list mode data handled by the
同時計数部42は、データ収集部25から呼吸センサ信号の情報(呼吸動データ)およびECG信号の情報(心電波形データ)が関連付けられたリストモードデータを受ける。
The
同時計数部42は、リストモードデータのうち、ガンマ線の入射時間差(対消滅ガンマ線の検出時間の差)が所定の時間ウインドウ幅(たとえば1ns以内など)にあり、かつ対消滅ガンマ線2つのそれぞれの入射エネルギーがともに所定のエネルギーウインドウ幅内にある組み合わせを抽出する。同時計数部42は、この抽出した組み合わせのリストモードデータ(以下、同時計数情報という)を心肺同期データ抽出部43および画像生成部44に与える。
The
なお、本実施形態では画像処理装置12の制御部31が同時計数情報を生成する場合の例について示したが、同時計数情報はたとえばスキャナ装置11のデータ収集部25によって生成されて画像処理装置12に与えられてもよい。
In the present embodiment, an example in which the
心肺同期データ抽出部43は、同時計数部42から受けた同時計数情報から呼吸動の同位相期間に属する同時計数情報を抽出する。また、心肺同期データ抽出部43は、抽出した呼吸動の同位相期間に属するデータについて、心臓の動きが所定の閾値より小さい期間(以下、第1の期間という)と閾値以上の期間(以下、第2の期間という)のいずれに属するデータであるかを特定する。また、心肺同期データ抽出部43は、抽出した呼吸動の同位相データから第2の期間に属するデータを除外して、すなわち抽出した同位相期間であるとともに心臓の動きが閾値より小さい期間(第1の期間)に属する同時計数情報を抽出して、画像生成部44に与える。
The cardiopulmonary synchronization
図4(a)は呼吸動の同位相期間に属する同時計数情報から肺全体51の画像を生成する従来の画像生成方法の一例を示す説明図であり、(b)は第1の期間に属する同時計数情報から肺全体51の画像を生成する本実施形態に係る画像生成方法の一例を示す説明図である。
FIG. 4A is an explanatory diagram showing an example of a conventional image generation method for generating an image of the
呼吸動の同位相期間に属する同時計数情報を抽出することにより、呼吸動による画像への影響を低減することができる。しかし、呼吸動の同位相期間には、心臓52の動き(拍動)が小さい第1の期間と心臓の動きが大きい期間(第2の期間)がある。このため、呼吸動の同位相期間に属する同時計数情報から肺全体51の画像を生成すると、図4(a)に示すように、心臓52の動きの影響を受けてしまい、心臓52の画像がぶれてしまうほか、心臓52の動きの影響を受ける気管支53などの肺野内の臓器の画像もぶれてしまう。
By extracting the coincidence counting information belonging to the same phase period of the respiratory motion, it is possible to reduce the influence of the respiratory motion on the image. However, in the same phase period of respiratory motion, there are a first period in which the movement (beat) of the
そこで、本実施形態に係る核医学診断装置10は、呼吸動および拍動による画像に対する影響を低減するため、図4(b)に示すように第1の期間に属する同時計数情報から肺全体51の画像を生成する。
Therefore, the nuclear
画像生成部44は、心肺同期データ抽出部43により抽出された第1の期間に属する同時計数情報にもとづいて同期再構成を行い、画像を生成して表示部32に表示させる。
The
次に、本実施形態に係る核医学診断装置10の動作の一例について説明する。
Next, an example of the operation of the nuclear medicine
図5は、図1に示す制御部31のCPUにより、呼吸動および拍動による撮像画像に対する影響を低減する際の手順を示すフローチャートである。図5において、Sに数字を付した符号は、フローチャートの各ステップを示す。
FIG. 5 is a flowchart showing a procedure when the influence of the respiratory motion and pulsation on the captured image is reduced by the CPU of the
まず、ステップS1において、FDGなどの薬剤を投与した患者Oを天板21に載置する。また、患者Oに呼吸センサ26および心電センサ27が取り付けられ、患者Oの呼吸センサ信号およびECG信号の情報がデータ収集部25に与えられる。
First, in step S <b> 1, the patient O who has received a drug such as FDG is placed on the
次に、ステップS2において、スキャン制御部41はユーザから入力部33を介してスキャン計画の実行指示を受けて、スキャン計画にもとづいてスキャナ装置11を制御することにより、スキャンを実行する。
Next, in step S2, the
次に、ステップS3において、データ収集部25は呼吸動データおよび心電波形データが関連付けられたリストモードデータを生成する。同時計数部42は、この呼吸動データおよび心電波形データが関連付けられたリストモードデータを取得する。
Next, in step S3, the
次に、ステップS4において、心肺同期データ抽出部43は、同時計数部42から受けた同時計数情報から、呼吸動の同位相期間に属するデータを抽出する。
Next, in step S <b> 4, the cardiopulmonary synchronization
次に、ステップS5において、心肺同期データ抽出部43は、抽出した呼吸動の同位相期間に属するデータについて、心臓の動きが所定の閾値より小さい期間(第1の期間)と閾値以上の期間(第2の期間)のいずれに属するデータであるかを特定する。
Next, in step S5, the cardiopulmonary synchronization
次に、ステップS6において、心肺同期データ抽出部43は、抽出した呼吸動の同位相データから第2の期間に属するデータを除外して画像生成部44に与える。
Next, in step S <b> 6, the cardiopulmonary synchronization
次に、ステップS7において、画像生成部44は心肺同期データ抽出部43から受けた第1の期間に属するデータ(同時計数情報)にもとづいて同期再構成を行い、画像を生成して表示部32に表示させる。
Next, in step S7, the
以上の手順により、呼吸動および拍動による撮像画像に対する影響を低減することができる。 By the above procedure, it is possible to reduce the influence on the captured image due to respiratory motion and pulsation.
本実施形態に係る核医学診断装置10は、呼吸動の同位相期間のうち心臓の動きがしきい値より小さい期間(第1の期間)における同時計数情報にもとづいて画像を生成する事ができる。このため、呼吸動の同位相期間に属するデータをすべて用いて画像を生成する場合に比べ、呼吸動および拍動による影響を低減した画像を生成することができる。このため、心臓周辺の肺野の画像の分解能を向上させることができ、ユーザはより正確に定量的な診断を行うことができる。
The nuclear
なお、本実施形態は核医学診断装置10としてPET装置を用いる場合について説明したが、本実施形態に係る技術はSPECT装置に適用することが可能である。SPECT装置に適用する場合、検出器23の構成が一般的なSPECT装置に用いられる構成とすればよい。また、SPECT装置を用いる場合、同時計数部42が不要であり、呼吸動データおよび心電波形データが関連付けられたリストモードデータは心肺同期データ抽出部43がデータ収集部25から取得すればよい。
(第2の実施形態)
In addition, although this embodiment demonstrated the case where a PET apparatus was used as the nuclear medicine
(Second Embodiment)
次に、本発明に係る核医学診断装置10および画像処理プログラムの第2実施形態について説明する。
Next, a second embodiment of the nuclear
本実施形態に係る核医学診断装置10は、第2の期間における同時計数情報のうち、心臓52から遠い部分のデータを同期再構成に用いるものである。
The nuclear medicine
図6は、第2実施形態に係る制御部31AのCPUによる機能実現部の構成例を示す概略的なブロック図である。
FIG. 6 is a schematic block diagram illustrating a configuration example of a function implementing unit by the CPU of the
この第2実施形態に示す核医学診断装置10は、制御部31Aの構成が第1実施形態に示す核医学診断装置10の制御部31と異なる。他の構成および作用については図1に示す核医学診断装置10と実質的に異ならないため、同じ構成には同一符号を付して説明を省略する。
The nuclear
図6に示すように、制御部31AのCPUは、画像処理プログラムによって、少なくともスキャン制御部41、同時計数部42、心肺同期データ抽出部43A、心臓領域特定部61および画像生成部44Aとして機能する。この各部は、RAMの所要のワークエリアをデータの一時的な格納場所として利用する。
As shown in FIG. 6, the CPU of the
心肺同期データ抽出部43Aは、同時計数部42から受けた同時計数情報から呼吸動の同位相期間に属する同時計数情報を抽出する。また、心肺同期データ抽出部43は、抽出した呼吸動の同位相期間に属するデータについて、第1の期間と第2の期間のいずれに属するデータであるかを特定する。
The cardiopulmonary synchronization
心臓領域特定部61は、同時計数情報の入射時間差にもとづいて対消滅ガンマ線対の発生位置を特定することにより、同時計数情報のうち心臓52の領域(心臓領域)に由来するデータを特定する。ここで、心臓領域の外郭は、心臓52の輪郭に設定されてもよいし、心臓52を含み心臓52の輪郭から所定の範囲に設定されてもよい。後者の場合、心臓領域は心臓52の近辺を含む。
The heart
画像生成部44Aは、第1の期間に属する同時計数情報については、そのまま同期再構成に用いる一方、第2の期間に属する同時計数情報については、心臓領域に由来するデータを間引くか、重みを小さくするか、または除外して、同期再構成に用いて画像を生成し、表示部32に表示させる。
The
ここで、「第2の期間に属する同時計数情報について心臓領域に由来するデータを間引く」とは、複数の第2の期間のうち、所定の期間おきのデータ(たとえば5期間で1つ)については心臓領域を含めたすべてのデータを同期再構成に用い、その他のデータからは心臓領域に由来するデータを削除することをいうものとする。
Here, “thinning the data derived from the heart region for the coincidence counting information belonging to the second period” means the data every predetermined period (for example, one in 5 periods) among the plurality of second periods. is used for synchronization reconstruct all of the data, including the heart臓領region, from other data shall refer to delete data from the cardiac region.
図7は、図6に示す制御部31AのCPUにより、呼吸動および拍動による撮像画像に対する影響を低減する際の手順を示すフローチャートである。図7において、Sに数字を付した符号は、フローチャートの各ステップを示す。図5と同等のステップには同一符号を付し、重複する説明を省略する。
FIG. 7 is a flowchart illustrating a procedure when the influence of the respiratory motion and the pulsation on the captured image is reduced by the CPU of the
心肺同期データ抽出部43Aにより、同時計数情報から呼吸動の同位相期間に属するデータが抽出されるとともに(ステップS4)、この呼吸動の同位相期間に属するデータについて第1の期間第2の期間のいずれに属するデータであるかが特定される(ステップS5)。
The cardiopulmonary synchronization
次に、ステップS21において、心臓領域特定部61は、同時計数情報の入射時間差にもとづいて対消滅ガンマ線対の発生位置を特定することにより、同時計数情報のうち心臓52の領域に由来するデータを特定する。
Next, in step S21, the heart
次に、ステップS22において、画像生成部44Aは、第1の期間に属する同時計数情報と、第2の期間に属する同時計数情報から心臓領域に由来するデータを間引くか、重みを小さくするか、または除外したものと、を同期再構成に用いて画像を生成し、表示部32に表示させる。
Next, in step S22, the
以上の手順によっても、呼吸動および拍動による撮像画像に対する影響を低減することができる。 Also by the above procedure, the influence on the captured image by respiratory motion and pulsation can be reduced.
第2実施形態に係る核医学診断装置10は、第1の期間における同時計数情報については第1実施形態と同様に同期再構成に利用するほか、第2の期間における同時計数情報については、全て除外してしまうのではなく心臓領域に係るデータを間引くか、重みを小さくするか、または除外して、同期再構成に用いる。
The nuclear
肺野領域のうち、心臓52から遠い部分については心臓の動きの影響を受けづらい。このため、第2の期間における同時計数情報のうち、心臓52から遠い部分のデータを同期再構成に用いることにより、第1実施形態に係る核医学診断装置10に比べ、肺野領域のうち心臓52から遠い部分の画像をより鮮明にすることができる。
(第3の実施形態)
Of the lung field region, the part far from the
(Third embodiment)
次に、本発明に係る核医学診断装置および画像処理プログラムの第3実施形態について説明する。 Next, a nuclear medicine diagnostic apparatus and an image processing program according to a third embodiment of the present invention will be described.
本実施形態に係る核医学診断装置10は、第1の期間と第2の期間において心臓領域に対して異なる強度の位置フィルタを適用して画像を生成するものである。
The nuclear
図8は、第3実施形態に係る制御部31BのCPUによる機能実現部の構成例を示す概略的なブロック図である。 FIG. 8 is a schematic block diagram illustrating a configuration example of a function realization unit by the CPU of the control unit 31B according to the third embodiment.
この第3実施形態に示す核医学診断装置10は、制御部31Bの構成が第1実施形態に示す核医学診断装置10の制御部31と異なる。他の構成および作用については図1に示す核医学診断装置10と実質的に異ならないため、同じ構成には同一符号を付して説明を省略する。また、心肺同期データ抽出部43Aおよび心臓領域特定部61の機能については第2実施形態に係る構成と実質的に異ならないため、説明を省略する。
In the nuclear medicine
一般に、核医学画像を生成する際には、同時計数情報から得られる位置情報に対してガウシアンフィルタなどの位置フィルタが施される。これは、飛来時間(TOF(Time Of Flight))の誤差を考慮した措置である。同時計数情報から得られる位置情報に対して位置フィルタを施して位置をぼかすことにより、位置情報をTOFの計測誤差を考慮したものに修正することができる。 In general, when a nuclear medicine image is generated, a position filter such as a Gaussian filter is applied to position information obtained from coincidence count information. This is a measure that takes into account an error in time of flight (TOF). By applying a position filter to the position information obtained from the coincidence information to blur the position, the position information can be corrected to take account of the TOF measurement error.
本実施形態に係る画像生成部44Bは、この位置フィルタを利用した位置のぼかしかた、すなわち位置フィルタの強度について、動きのある心臓領域と他の領域(たとえば肺野領域)とで異ならせる。 The image generation unit 44B according to the present embodiment varies the position blurring using this position filter, that is, the strength of the position filter, between a heart region in motion and another region (for example, a lung field region).
具体的には、画像生成部44Bは、第1の期間に属する同時計数情報から得られる位置情報については所定のカーネル幅で位置フィルタを施し、第2の期間に属する同時計数情報から得られる位置情報については、肺野領域に対しては所定のカーネル幅を用いる一方心臓領域に対しては所定のカーネル幅とは異なるカーネル幅で位置フィルタを施す。位置フィルタとしてガウシアンフィルタを用いる場合、ガウス関数の分散値を変更することによりカーネル幅(フィルタ強度)を変更することができる。 Specifically, the image generation unit 44B performs position filtering on the position information obtained from the coincidence counting information belonging to the first period with a predetermined kernel width, and the position obtained from the coincidence counting information belonging to the second period. For the information, a predetermined kernel width is used for the lung field region, while a position filter is applied to the heart region with a kernel width different from the predetermined kernel width. When a Gaussian filter is used as the position filter, the kernel width (filter strength) can be changed by changing the variance value of the Gaussian function.
図9は、図8に示す制御部31BのCPUにより、呼吸動および拍動による撮像画像に対する影響を低減する際の手順を示すフローチャートである。図9において、Sに数字を付した符号は、フローチャートの各ステップを示す。図5および図7と同等のステップには同一符号を付し、重複する説明を省略する。 FIG. 9 is a flowchart illustrating a procedure when the influence of the respiratory motion and the pulsation on the captured image is reduced by the CPU of the control unit 31B illustrated in FIG. In FIG. 9, the code | symbol which attached | subjected the number to S shows each step of a flowchart. Steps equivalent to those in FIGS. 5 and 7 are denoted by the same reference numerals, and redundant description is omitted.
ステップS31において、画像生成部44Bは、第1の期間に属する同時計数情報から得られる位置情報については所定のカーネル幅で位置フィルタを施し、第2の期間に属する同時計数情報から得られる位置情報については、肺野領域に対しては所定のカーネル幅を用いる一方心臓領域に対しては所定のカーネル幅とは異なるカーネル幅で位置フィルタを施す。そして、画像生成部44Bは、位置フィルタを施したデータを用いて画像を生成し表示部32に表示させる。
In step S31, the image generation unit 44B performs position filtering on the position information obtained from the coincidence counting information belonging to the first period with a predetermined kernel width, and the position information obtained from the coincidence counting information belonging to the second period. For the lung field region, a predetermined kernel width is used, while for the heart region, a position filter is applied with a kernel width different from the predetermined kernel width. Then, the image generation unit 44B generates an image using the data on which the position filter has been applied, and causes the
以上の手順によっても、呼吸動および拍動による撮像画像に対する影響を低減することができる。 Also by the above procedure, the influence on the captured image by respiratory motion and pulsation can be reduced.
第3実施形態に係る核医学診断装置10は、第2実施形態に係る核医学診断装置10と同様に、第1の期間および第2の期間における同時計数情報を用いる。このため、第1実施形態に係る核医学診断装置10に比べ、肺野領域のうち心臓52から遠い部分の画像をより鮮明にすることができる。また、第2の期間に属する同時計数情報から得られる位置情報については、肺野領域に対しては所定のカーネル幅を用いる一方心臓領域に対しては所定のカーネル幅とは異なるカーネル幅で位置フィルタを施す。このため、単に呼吸動の同位相期間に属するデータをすべて用いて画像を生成する場合に比べ、画質を向上させることができる。
The nuclear
たとえば、第2の期間(心臓の動きが大きい期間)において、心臓領域のカーネル幅を所定のカーネル幅よりもシャープにする場合、拍動のためにぼやけがちとなる心臓領域の画像をシャープにすることができ、心臓領域と肺野領域とで均一な画像にすることができる。 For example, when the kernel width of the heart region is sharpened more than a predetermined kernel width in the second period (a period in which the heart motion is large), an image of the heart region that tends to be blurred due to pulsation is sharpened. It is possible to obtain a uniform image in the heart region and the lung region.
一方、第2の期間(心臓の動きが大きい期間)において、心臓領域のカーネル幅を所定のカーネル幅よりもブロードにする場合、心臓領域の画像はぼやけるが、心臓52の動き全体を捉えた画像(心臓の動きの誤差を含めた存在確率を反映した画像)を生成することができる。
(第4の実施形態)
On the other hand, when the kernel width of the heart region is broader than the predetermined kernel width in the second period (period in which the heart motion is large), the image of the heart region is blurred, but the image capturing the entire motion of the heart 52 (An image reflecting the existence probability including an error of heart motion) can be generated.
(Fourth embodiment)
次に、本発明に係る核医学診断装置および画像処理プログラムの第4実施形態について説明する。 Next, a nuclear medicine diagnosis apparatus and an image processing program according to a fourth embodiment of the present invention will be described.
本実施形態に係る核医学診断装置10は、第2実施形態に係る核医学診断装置10と第3実施形態に係る核医学診断装置10とを組み合わせたものである。
A nuclear
図10は、第3実施形態に係る制御部31BのCPUによる機能実現部の構成例を示す概略的なブロック図である。 FIG. 10 is a schematic block diagram illustrating a configuration example of the function realization unit by the CPU of the control unit 31B according to the third embodiment.
この第3実施形態に示す核医学診断装置10は、制御部31Cの構成が第1実施形態に示す核医学診断装置10の制御部31と異なる。他の構成および作用については図1に示す核医学診断装置10と実質的に異ならないため、同じ構成には同一符号を付して説明を省略する。また、心肺同期データ抽出部43Aおよび心臓領域特定部61の機能については第2実施形態および第3実施形態に係る構成と実質的に異ならないため、説明を省略する。
In the nuclear medicine
また、図11は、図10に示す制御部31CのCPUにより、呼吸動および拍動による撮像画像に対する影響を低減する際の手順を示すフローチャートである。図11において、Sに数字を付した符号は、フローチャートの各ステップを示す。図5、図7および図9と同等のステップには同一符号を付し、重複する説明を省略する。
FIG. 11 is a flowchart illustrating a procedure when the influence of the respiratory motion and the pulsation on the captured image is reduced by the CPU of the
画像生成部44Cは、第1の期間に属する同時計数情報については、そのまま同期再構成に用いる一方、第2の期間に属する同時計数情報については、心臓領域に由来するデータを間引くか、重みを小さくする(図11のステップS41)。第2実施形態に係る画像生成部44A都の相違点は、心臓領域に由来するデータを削除することがない点である。
The image generation unit 44C uses the coincidence counting information belonging to the first period as it is for the synchronization reconstruction, while the coincidence counting information belonging to the second period is thinned or weighted for data derived from the heart region. It is made smaller (step S41 in FIG. 11). The difference between the
また、画像生成部44Cは、第3実施形態に係る画像生成部44Bと同様に、第1の期間に属する同時計数情報から得られる位置情報については所定のカーネル幅で位置フィルタを施し、第2の期間に属する同時計数情報から得られる位置情報については、肺野領域に対しては所定のカーネル幅を用いる一方心臓領域に対しては所定のカーネル幅とは異なるカーネル幅で位置フィルタを施して、画像を生成し表示部32に表示させる(図11のステップS42)。 Similarly to the image generation unit 44B according to the third embodiment, the image generation unit 44C applies a position filter with a predetermined kernel width to the position information obtained from the coincidence counting information belonging to the first period, For the position information obtained from the coincidence count information belonging to this period, a predetermined kernel width is used for the lung field region, while a position filter is applied to the heart region with a kernel width different from the predetermined kernel width. Then, an image is generated and displayed on the display unit 32 (step S42 in FIG. 11).
第4実施形態に係る核医学診断装置10は、第1の期間および第2の期間における同時計数情報を用いる。このため、第1実施形態に係る核医学診断装置10に比べ、肺野領域のうち心臓52から遠い部分の画像をより鮮明にすることができる。また、第2の期間の心臓領域のデータについて、間引くまたは重み付けを小さくするだけでなく、位置フィルタの強度を肺野領域と異ならせることにより、心臓領域の画質をより向上させることができる。
The nuclear
なお、本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれる。 In addition, although some embodiment of this invention was described, these embodiment is shown as an example and is not intending limiting the range of invention. These novel embodiments can be implemented in various other forms, and various omissions, replacements, and changes can be made without departing from the scope of the invention. These embodiments and modifications thereof are included in the scope and gist of the invention, and are included in the invention described in the claims and the equivalents thereof.
また、本発明の実施形態では、フローチャートの各ステップは、記載された順序に沿って時系列的に行われる処理の例を示したが、必ずしも時系列的に処理されなくとも、並列的あるいは個別実行される処理をも含むものである。 Further, in the embodiment of the present invention, each step of the flowchart shows an example of processing that is performed in time series in the order described. The process to be executed is also included.
10 核医学診断装置
23 検出器
25 データ収集部
26 呼吸センサ
27 心電センサ
31、31A、31B、31C 制御部
32 表示部
42 同時計数部
43、43A 心肺同期データ抽出部
44、44A、44B、44C 画像生成部
52 心臓
61 心臓領域特定部
DESCRIPTION OF
Claims (6)
前記被検体に取り付けられた呼吸センサおよび心電センサの出力データと、前記複数の検出器の出力データと、を収集するデータ収集部と、
前記データ収集部により収集されたデータにもとづいて、呼吸動の同位相期間であるとともに心臓の動きが閾値より小さい期間である第1の期間における前記複数の検出器の出力データにもとづく画像を生成する画像生成部と、
前記複数の検出器の出力データにもとづいて、対消滅ガンマ線対の検出時間差から心臓の領域に由来するデータを特定する心臓領域特定部と、
を備え、
前記画像生成部は、
前記第1の期間における前記複数の検出器の出力データと、前記呼吸動の同位相期間であるとともに心臓の動きが前記閾値より大きい期間である複数の第2の期間のうち所定数の期間における前記複数の検出器の出力データから前記心臓の領域に由来するデータを削除することにより、複数の前記第2の期間における前記複数の検出器の出力データから前記心臓の領域に由来するデータを間引いたデータと、にもとづいて画像を生成する、
核医学診断装置。 A plurality of detectors for detecting gamma rays emitted from the radioisotope administered to the subject;
A data collection unit that collects output data of a respiratory sensor and an electrocardiographic sensor attached to the subject, and output data of the plurality of detectors;
Based on the data collected by the data collection unit, an image is generated based on the output data of the plurality of detectors in the first period in which the cardiac motion is in the same phase period and the heart motion is smaller than the threshold value. An image generation unit to
Based on the output data of the plurality of detectors, a heart region specifying unit for specifying data derived from the heart region from the detection time difference of the pair annihilation gamma ray pair;
Equipped with a,
The image generation unit
The output data of the plurality of detectors in the first period and the in-phase period of the respiratory motion and in a predetermined number of periods among a plurality of second periods in which the heart motion is greater than the threshold By deleting data derived from the heart region from the output data of the plurality of detectors, data derived from the heart region is thinned out from the output data of the plurality of detectors in the plurality of second periods. Image based on the collected data,
Nuclear medicine diagnostic equipment.
少なくとも前記被検体の肺野領域および前記心臓の領域に由来する対消滅ガンマ線を検出し、
前記画像生成部は、
前記第1の期間における前記複数の検出器の出力データに対して所定の強度で位置フィルタを施したデータと、前記第2の期間における前記複数の検出器の出力データのうち前記肺野領域に由来するデータに対して前記所定の強度で位置フィルタを施す一方前記心臓の領域に由来するデータに対しては前記所定の強度とは異なる強度で位置フィルタを施したデータと、にもとづいて画像を生成する、
請求項1記載の核医学診断装置。 The plurality of detectors are:
Detecting pair annihilation gamma rays originating from at least the lung field region of the subject and the heart region;
The image generation unit
In the lung field region of the output data of the plurality of detectors in the first period, the data subjected to position filtering with a predetermined intensity, and the output data of the plurality of detectors in the second period A position filter is applied to the derived data at the predetermined intensity, while an image is generated based on the data obtained by applying the position filter at an intensity different from the predetermined intensity for the data derived from the heart region. Generate,
The nuclear medicine diagnostic apparatus according to claim 1 .
をさらに備え、
前記複数の検出器は、
少なくとも前記被検体の肺野領域および前記心臓の領域に由来する対消滅ガンマ線を検出し、
前記画像生成部は、
前記第1の期間における前記複数の検出器の出力データに対して所定の強度で位置フィルタを施したデータと、前記呼吸動の同位相期間であるとともに心臓の動きが前記閾値より大きい期間である第2の期間における前記複数の検出器の出力データのうち前記肺野領域に由来するデータに対して前記所定の強度で位置フィルタを施す一方前記心臓の領域に由来するデータに対しては前記所定の強度とは異なる強度で位置フィルタを施したデータと、にもとづいて画像を生成する、
請求項1記載の核医学診断装置。 Based on the output data of the plurality of detectors, a heart region specifying unit that specifies data derived from the heart region from the detection time difference of the pair annihilation gamma ray pair;
Further comprising
The plurality of detectors are:
Detecting pair annihilation gamma rays originating from at least the lung field region of the subject and the heart region;
The image generation unit
The data obtained by applying a position filter to the output data of the plurality of detectors at a predetermined intensity in the first period and the in-phase period of the respiratory motion and the period in which the heart motion is greater than the threshold value Among the output data of the plurality of detectors in the second period, the position filter is applied to the data derived from the lung field region at the predetermined intensity, while the data derived from the heart region is applied to the predetermined data. Generate an image based on the data that has been subjected to position filtering at an intensity different from the intensity of
The nuclear medicine diagnostic apparatus according to claim 1.
請求項2または3に記載の核医学診断装置。 The position filter is a Gaussian filter, and the intensity is a variance value of a Gaussian function used in the Gaussian filter.
The nuclear medicine diagnostic apparatus according to claim 2 or 3 .
前記複数の検出器の出力データにもとづいて、前記検出器へガンマ線が入射するごとに入射したガンマ線のエネルギーおよび入射時間の情報を取得するリストモードでデータを収集する、
請求項1ないし4のいずれか1項に記載の核医学診断装置。 The data collection unit
Based on the output data of the plurality of detectors, data is collected in a list mode for acquiring information on the energy and incident time of the incident gamma rays every time gamma rays are incident on the detectors.
The nuclear medicine diagnostic apparatus according to any one of claims 1 to 4 .
被検体に投与された放射性同位元素から放射されるガンマ線を検出する複数の検出器の出力データと前記被検体に取り付けられた呼吸センサおよび心電センサの出力データとにもとづいて、呼吸動の同位相期間であるとともに心臓の動きが閾値より小さい期間である第1の期間における前記複数の検出器の出力データを抽出するステップと、
前記複数の検出器の出力データにもとづいて、対消滅ガンマ線対の検出時間差から心臓の領域に由来するデータを特定するステップと、
前記呼吸動の同位相期間であるとともに心臓の動きが前記閾値より大きい期間である複数の第2の期間のうち、所定数の期間における前記複数の検出器の出力データから前記心臓の領域に由来するデータを削除することにより、複数の前記第2の期間における前記複数の検出器の出力データから前記心臓の領域に由来するデータを間引いたデータを生成するステップと、
前記第1の期間における前記複数の検出器の出力データと、複数の前記第2の期間における前記複数の検出器の出力データから前記心臓の領域に由来するデータを間引いたデータと、にもとづいて画像を生成するステップと、
を実行させる画像処理プログラム。 On the computer,
Based on the output data of a plurality of detectors that detect gamma rays emitted from the radioisotope administered to the subject and the output data of the respiratory sensor and the electrocardiographic sensor attached to the subject, the respiratory motion is the same. Extracting output data of the plurality of detectors in a first period that is a phase period and wherein the heart motion is less than a threshold;
Identifying data derived from the heart region from the detection time difference of the pair of annihilation gamma rays based on the output data of the plurality of detectors;
From the output data of the plurality of detectors in a predetermined number of periods among the plurality of second periods in which the respiratory motion is the same phase period and the heart motion is greater than the threshold value, derived from the heart region Generating data obtained by thinning out data derived from the heart region from output data of the plurality of detectors in the plurality of second periods by deleting data to be
Based on output data of the plurality of detectors in the first period and data obtained by thinning out data derived from the heart region from output data of the plurality of detectors in the plurality of second periods. Generating an image; and
An image processing program for executing
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012249536A JP6091852B2 (en) | 2012-11-13 | 2012-11-13 | Nuclear medicine diagnostic apparatus and image processing program |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012249536A JP6091852B2 (en) | 2012-11-13 | 2012-11-13 | Nuclear medicine diagnostic apparatus and image processing program |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2014098583A JP2014098583A (en) | 2014-05-29 |
JP6091852B2 true JP6091852B2 (en) | 2017-03-08 |
Family
ID=50940711
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2012249536A Active JP6091852B2 (en) | 2012-11-13 | 2012-11-13 | Nuclear medicine diagnostic apparatus and image processing program |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP6091852B2 (en) |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5526456A (en) * | 1978-08-17 | 1980-02-25 | Hitachi Medical Corp | Scinti-camera control unit |
JPH0627834B2 (en) * | 1989-10-30 | 1994-04-13 | 株式会社島津製作所 | Biometric data collection device |
US5632276A (en) * | 1995-01-27 | 1997-05-27 | Eidelberg; David | Markers for use in screening patients for nervous system dysfunction and a method and apparatus for using same |
WO2008127368A2 (en) * | 2007-02-07 | 2008-10-23 | Koninklijke Philips Electronics, N.V. | Motion compensation in quantitative data analysis and therapy |
JP6243121B2 (en) * | 2009-12-10 | 2017-12-06 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | Method and apparatus for motion detection and correction in imaging scans using time-of-flight information |
-
2012
- 2012-11-13 JP JP2012249536A patent/JP6091852B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
JP2014098583A (en) | 2014-05-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5782020B2 (en) | Functional imaging | |
WO2014189134A1 (en) | Nuclear medicine diagnosis device and nuclear medicine image formation method | |
CN104331914B (en) | The method and apparatus for rebuilding cardiac phase image | |
US20190133542A1 (en) | Systems and methods for data-driven respiratory gating in positron emission tomography | |
US8515148B2 (en) | Geometrical transformations preserving list-mode format | |
JP2020076584A (en) | Medical image processing device | |
JP6668085B2 (en) | Medical image diagnostic apparatus and medical image processing apparatus | |
JP6526428B2 (en) | Medical image processing apparatus, medical image processing method and medical image diagnostic apparatus | |
CN102247156A (en) | Image processing apparatus and X-ray CT system | |
JP6091852B2 (en) | Nuclear medicine diagnostic apparatus and image processing program | |
JP6682243B2 (en) | Medical image diagnostic apparatus, image processing apparatus, and image processing method | |
US9239393B2 (en) | Positron emission computed tomography apparatus and image processing apparatus | |
US20230022425A1 (en) | Apparatus, system, method and computer probram for providing a nuclear image of a region of interest of a patient | |
Mohammadi et al. | Motion in nuclear cardiology imaging: types, artifacts, detection and correction techniques | |
JP4640143B2 (en) | Diagnostic imaging support device | |
CN106548464B (en) | A kind of method and device constructing phase image | |
JP6425885B2 (en) | Nuclear medicine diagnostic device, image processing device and image processing program | |
JP6249317B2 (en) | Tomographic image display device | |
US20200029928A1 (en) | Systems and methods for improved motion correction | |
JP2016020866A (en) | Computer program, image processing device, and method | |
US20160203621A1 (en) | Medical diagnostic imaging apparatus, image processing apparatus and image processing method | |
JP6147512B2 (en) | Nuclear medicine diagnostic apparatus, image processing apparatus, and image reconstruction program | |
US20240233211A1 (en) | Method and apparatus for cardiac signal removal in data-driven respiratory gating in positron emission tomography systems | |
JP6682373B2 (en) | SPECT equipment | |
Buechel et al. | Single-photon emission computed tomography |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20151008 |
|
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A711 Effective date: 20160506 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20160914 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20160920 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20161118 |
|
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: 20170110 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20170208 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6091852 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |