JPS6346138A - 同期画像再構成装置 - Google Patents
同期画像再構成装置Info
- Publication number
- JPS6346138A JPS6346138A JP61189173A JP18917386A JPS6346138A JP S6346138 A JPS6346138 A JP S6346138A JP 61189173 A JP61189173 A JP 61189173A JP 18917386 A JP18917386 A JP 18917386A JP S6346138 A JPS6346138 A JP S6346138A
- Authority
- JP
- Japan
- Prior art keywords
- phase
- projection data
- image
- data
- scanner
- 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
Links
- 230000001360 synchronised effect Effects 0.000 title claims description 12
- 230000000737 periodic effect Effects 0.000 claims description 4
- 238000013480 data collection Methods 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 16
- 238000001514 detection method Methods 0.000 description 4
- 239000002872 contrast media Substances 0.000 description 3
- 230000003205 diastolic effect Effects 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000002595 magnetic resonance imaging Methods 0.000 description 2
- 238000000034 method Methods 0.000 description 2
- 230000029058 respiratory gaseous exchange Effects 0.000 description 2
- MXBCYQUALCBQIJ-RYVPXURESA-N (8s,9s,10r,13s,14s,17r)-13-ethyl-17-ethynyl-11-methylidene-1,2,3,6,7,8,9,10,12,14,15,16-dodecahydrocyclopenta[a]phenanthren-17-ol;(8r,9s,13s,14s,17r)-17-ethynyl-13-methyl-7,8,9,11,12,14,15,16-octahydro-6h-cyclopenta[a]phenanthrene-3,17-diol Chemical compound OC1=CC=C2[C@H]3CC[C@](C)([C@](CC4)(O)C#C)[C@@H]4[C@@H]3CCC2=C1.C1CC[C@@H]2[C@H]3C(=C)C[C@](CC)([C@](CC4)(O)C#C)[C@@H]4[C@@H]3CCC2=C1 MXBCYQUALCBQIJ-RYVPXURESA-N 0.000 description 1
- 241000257465 Echinoidea Species 0.000 description 1
- AQSMZTIEJMZQEC-DCAQKATOSA-N Pro-His-Ser Chemical compound C1C[C@H](NC1)C(=O)N[C@@H](CC2=CN=CN2)C(=O)N[C@@H](CO)C(=O)O AQSMZTIEJMZQEC-DCAQKATOSA-N 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000010894 electron beam technology Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000000004 hemodynamic effect Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 210000003127 knee Anatomy 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000004165 myocardium Anatomy 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 230000000241 respiratory effect Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/54—Control of apparatus or devices for radiation diagnosis
- A61B6/541—Control of apparatus or devices for radiation diagnosis involving acquisition triggered by a physiological signal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/412—Dynamic
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
- Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10S378/00—X-ray or gamma ray systems or devices
- Y10S378/901—Computer tomography program or processor
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Theoretical Computer Science (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Pulmonology (AREA)
- Physiology (AREA)
- General Physics & Mathematics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
［発明の目的］
（産業上の利用分野）
本発明は、生体等の周期的な動きのある物体について得
られた画像データに基づき、時間軸方向で複数枚の同期
画像を再構成する同期画像再構成装置に関する。
られた画像データに基づき、時間軸方向で複数枚の同期
画像を再構成する同期画像再構成装置に関する。
（従来の技術）
例えば、人体の断層面を画像化する医用診断装置として
Ｘ線ＣＴ装置、磁気共鳴イメージング装置（〜ＩＲＩ
）を挙げることができる。
Ｘ線ＣＴ装置、磁気共鳴イメージング装置（〜ＩＲＩ
）を挙げることができる。
以下、Ｘ線ＣＴ装置を例に挙げ、心拍によって周期的に
動く部位を画像化する場合について説明する。尚、ＭＲ
Ｉについても本質的に同様であり、動きのある部位とし
ては心臓の動きに基因するものの他呼吸によって動く腹
部等についても同様と考えられる。
動く部位を画像化する場合について説明する。尚、ＭＲ
Ｉについても本質的に同様であり、動きのある部位とし
ては心臓の動きに基因するものの他呼吸によって動く腹
部等についても同様と考えられる。
造影剤を注入し、それによる造影状態の変化あるいは造
影剤の流れを経時的に追跡するダイナミックスキャンは
Ｘ線ＣＴ装置でよく行われる検査法である。
影剤の流れを経時的に追跡するダイナミックスキャンは
Ｘ線ＣＴ装置でよく行われる検査法である。
ところで、心臓は動く部位でおるので、スキャン時間の
遅い一般の装置では、動きによるアーチファクトやボケ
を避けて診断価値のある画像を得るために心電同期スキ
ャンを行っている。
遅い一般の装置では、動きによるアーチファクトやボケ
を避けて診断価値のある画像を得るために心電同期スキ
ャンを行っている。
これは、スキャンを例えば１０回位繰り返し、心電波形
も同時に記録し、その後、心電波形の特定位相状態（あ
る狭い幅を持つ）、に一致するタイミングで得られた投
影データを複数回のスキャンで得られた投影データ群の
中から寄せ集めて特定位相状態の画像を形成するもので
ある。
も同時に記録し、その後、心電波形の特定位相状態（あ
る狭い幅を持つ）、に一致するタイミングで得られた投
影データを複数回のスキャンで得られた投影データ群の
中から寄せ集めて特定位相状態の画像を形成するもので
ある。
しかし、上記のようにして得られた一枚の画像は、１０
回のスキャンに亘って得られた投影データから作られた
ものであるため、数１０秒〜１０Ｏ秒の時間幅を待ち、
速い血行動態等は追跡不能となってしまう。勿論、従来
のＸ線ＣＴ装置のスキャン時間は数秒であるから、この
ような心電同期スキャンなくしてはアーチファクトがひ
どく、ダイナミックスキャン以前に診断価値のある画像
が得られない。
回のスキャンに亘って得られた投影データから作られた
ものであるため、数１０秒〜１０Ｏ秒の時間幅を待ち、
速い血行動態等は追跡不能となってしまう。勿論、従来
のＸ線ＣＴ装置のスキャン時間は数秒であるから、この
ような心電同期スキャンなくしてはアーチファクトがひ
どく、ダイナミックスキャン以前に診断価値のある画像
が得られない。
心電同期に因らずに心臓の診断価値のめる両膝を得るス
キャナとしてスリップリングを有するＸ線ＣＴ装置では
、１秒程度の高速スキャンを数１０秒に亘って繰り返し
行うことができる。
キャナとしてスリップリングを有するＸ線ＣＴ装置では
、１秒程度の高速スキャンを数１０秒に亘って繰り返し
行うことができる。
しかし、このようなＸ線ＣＴ装置にあっても下記のよう
な問題がある。
な問題がある。
１秒でスキャンが完了し、これを継続できるものとし、
一つの画像は２／３回転分のデータから再構成するもの
とし、心拍は１秒よりヤヤ短いものとした場合、画像の
収集タイミング、心電Ｒ波のタイミングは第１１図のよ
うになる。
一つの画像は２／３回転分のデータから再構成するもの
とし、心拍は１秒よりヤヤ短いものとした場合、画像の
収集タイミング、心電Ｒ波のタイミングは第１１図のよ
うになる。
同図において、Ｒ波タイミングは略々収縮期であり、Ｒ
波とＲ波との中間は略々拡張期である。
波とＲ波との中間は略々拡張期である。
従って、同図における画像■、■は心臓の像はやヤ小さ
く、画像■、■、■、■、０は心臓の像はやや大きく観
測される。このため、心臓の動きはもとの正しい動きと
は全く異なるように見え、ざらに悪いことには、ダイナ
ミックスタデイでは関心領域（ＲＯＩ）を数ケ所に設定
し、各ＲＯＩ内のＣＴ値の変化を追跡するのであるが、
ＲＯＩは画像マトリックス内の固有の座標に対して設定
するものであるから、画像が動いてしまえば特定場所の
ＣＴ値の変化を観測することができなくなる。
く、画像■、■、■、■、０は心臓の像はやや大きく観
測される。このため、心臓の動きはもとの正しい動きと
は全く異なるように見え、ざらに悪いことには、ダイナ
ミックスタデイでは関心領域（ＲＯＩ）を数ケ所に設定
し、各ＲＯＩ内のＣＴ値の変化を追跡するのであるが、
ＲＯＩは画像マトリックス内の固有の座標に対して設定
するものであるから、画像が動いてしまえば特定場所の
ＣＴ値の変化を観測することができなくなる。
（発明が解決しようとする問題点）
上述したように、近年、１回のスキャン時間が短く、し
かも複数回に亘って継続してスキャンを行うことができ
る装置が開発されていながらも、動きのある部位に対し
て位相的に一致した複数枚の同期画像を形成することが
できなかった。
かも複数回に亘って継続してスキャンを行うことができ
る装置が開発されていながらも、動きのある部位に対し
て位相的に一致した複数枚の同期画像を形成することが
できなかった。
そこで、本発明の目的とするところは上述した従来の欠
点を解消し、動きのある部位について位相的に一致した
複数枚の同期画像を再構成することができ、特にダイナ
ミックスタデイに最適な同期画像再構成装置を提供する
ことにある。
点を解消し、動きのある部位について位相的に一致した
複数枚の同期画像を再構成することができ、特にダイナ
ミックスタデイに最適な同期画像再構成装置を提供する
ことにある。
［発明の構成］
（問題点を解決するための手段）
本発明は、所定の周期で変動する部位について、複数周
期に亘って前記部位をスキャンして投影データを収集す
るスキャナと、前記データ収集中に亘って前記部位の周
期的な動きと対応する信号を検出する信号検出手段と、
前記スキャナからのデータ及び前記信号検出手段からの
信号を関連付けて記憶する記憶部と、前記変動周期中の
特定位相を設定する位相設定部と、設定された特定（ｌ
相に対応するデータを含んで少なくとも一枚の画像を再
構成するのに必要な投影データを各画像毎に前記記憶部
より読み出し、時間軸方向に複数枚の同期画像を再構成
する画像再構成部とを有して同期画像再構成装置を構成
している。
期に亘って前記部位をスキャンして投影データを収集す
るスキャナと、前記データ収集中に亘って前記部位の周
期的な動きと対応する信号を検出する信号検出手段と、
前記スキャナからのデータ及び前記信号検出手段からの
信号を関連付けて記憶する記憶部と、前記変動周期中の
特定位相を設定する位相設定部と、設定された特定（ｌ
相に対応するデータを含んで少なくとも一枚の画像を再
構成するのに必要な投影データを各画像毎に前記記憶部
より読み出し、時間軸方向に複数枚の同期画像を再構成
する画像再構成部とを有して同期画像再構成装置を構成
している。
（作 用）
周期的に変動する部位についての投影データと、該変動
に対応する信号とを複数周期に亘って予め収集しておき
、これを記憶部に相互に関連付けて記憶しておく。この
後、特定位相が設定されると、画像再構成部は設定され
た特定位相に対応するデータを含んで少なくとも一枚の
画像を再構成するのに必要な投影データを各画像毎に前
記記゛恩部より読み出し、該データに基づき画像の再構
成を行うことになる。このようにして得られた各画像は
、位相が一致しており、結局、時間軸方向に複数枚の同
期画像が再構成されたことになる。
に対応する信号とを複数周期に亘って予め収集しておき
、これを記憶部に相互に関連付けて記憶しておく。この
後、特定位相が設定されると、画像再構成部は設定され
た特定位相に対応するデータを含んで少なくとも一枚の
画像を再構成するのに必要な投影データを各画像毎に前
記記゛恩部より読み出し、該データに基づき画像の再構
成を行うことになる。このようにして得られた各画像は
、位相が一致しており、結局、時間軸方向に複数枚の同
期画像が再構成されたことになる。
尚、設定される特定位相としては、記憶部より読み出さ
れる投影データの収集時間幅の中心に一致する中心位相
の他、記憶部より読み出される投影データの収集開始時
間に一致する第１の特定位相と、収集終了時間に一致す
る第２の特定位相との２種類の位相を設定するものでも
よい。
れる投影データの収集時間幅の中心に一致する中心位相
の他、記憶部より読み出される投影データの収集開始時
間に一致する第１の特定位相と、収集終了時間に一致す
る第２の特定位相との２種類の位相を設定するものでも
よい。
（実施例）
以下、本発明を図示の実施例を参照して具体的に説明す
る。
る。
く第１実施例〉
第１図は、本発明の第１実施例装置のブロック図である
。
。
同図において、１はスキャナであり、本実施例では例え
ば第４世代のＸ線ＣＴスキャナを用いている。この第４
世代のＸ線ＣＴスキャナは、検出器を３６０度に亘って
配置し、スリップリングを用いて管電圧を供給しつつＸ
線管のみを被検体の周囲に回転して高速スキャンを連続
して実行することができるスキャナの一例でおる。
ば第４世代のＸ線ＣＴスキャナを用いている。この第４
世代のＸ線ＣＴスキャナは、検出器を３６０度に亘って
配置し、スリップリングを用いて管電圧を供給しつつＸ
線管のみを被検体の周囲に回転して高速スキャンを連続
して実行することができるスキャナの一例でおる。
２は周期的に変動する部位の動きに対応した信号を検出
する信号検出手段の一例としての心電計であり、前記ス
キャナで複数周期に亘ってデータ収束している間の心電
波形を計測するものである。
する信号検出手段の一例としての心電計であり、前記ス
キャナで複数周期に亘ってデータ収束している間の心電
波形を計測するものである。
３はデータ記憶部としてのディスクでおり、前記スキャ
ナ１で収集された被検体に関する投影データと、前記心
電計で収集された心電波形とを相互に関連付けて記憶し
ている。
ナ１で収集された被検体に関する投影データと、前記心
電計で収集された心電波形とを相互に関連付けて記憶し
ている。
４は位相設定部の一例としての心電位相入力部であり、
本実施例では一枚の画像を再構成するのに必要な投影デ
ータに対し、該投影データの収集時間の中心時刻に対応
する心電位相φを入力するものでおる。
本実施例では一枚の画像を再構成するのに必要な投影デ
ータに対し、該投影データの収集時間の中心時刻に対応
する心電位相φを入力するものでおる。
５は対話装置であり、前記心電位相入力部４と前記プロ
セッサ６との間の対話を実行するものである。
セッサ６との間の対話を実行するものである。
前記プロセッサ６は、前記スキャナ１．心電計２からの
データを前記ディスク３内に格納し、かつ、前記心電位
相入力部４より入力された位相φに相当する時刻の前記
Ｘ線管の回転角を中心に、一枚のＣＴ両画像構成するの
に必要な回転角θ分の投影データを前記ディスク３より
読み出して、該データに基づきＣＴ両画像再構成を実行
するよ。
データを前記ディスク３内に格納し、かつ、前記心電位
相入力部４より入力された位相φに相当する時刻の前記
Ｘ線管の回転角を中心に、一枚のＣＴ両画像構成するの
に必要な回転角θ分の投影データを前記ディスク３より
読み出して、該データに基づきＣＴ両画像再構成を実行
するよ。
うになっている。
以上のように構成された装置の作用について説明する。
第２図は、上記実施例装置の動作を説明するための図で
、前記スキャナ１では同図の架台回転角に示すように複
数回に亘ってＸ線管を回転して被検体をスキャンして投
影データを収集する。
、前記スキャナ１では同図の架台回転角に示すように複
数回に亘ってＸ線管を回転して被検体をスキャンして投
影データを収集する。
このデータはプロセッサ６を介してディスク３に格納さ
れることになる。
れることになる。
一方、前記スキャナ１でのデータ収集時に亘って前記心
電計２では第２図に示すような被検体の心電波形を収集
し、この心電波形も前記プロセッサ６を介して前記ディ
スク３に格納されることになる。
電計２では第２図に示すような被検体の心電波形を収集
し、この心電波形も前記プロセッサ６を介して前記ディ
スク３に格納されることになる。
上記のような投影データ、心電波形の収集後、前記心電
位相入力部４より中心位相φが入力されると、対話装置
５を介してこの位相φがプロセッサ６に入力される。こ
の位相として例えばφ＝Ｏ即らＲ波のピーク位置が中心
位相として入力されると、前記プロセッサ６は、位相φ
に相当する時刻の前記Ｘ線管の回転角を中心に、一枚の
ＣＴ両画像構成するのに必要な回転角θ分の投影データ
を、前記ディスク３の内の一連の長い投影データの中か
ら読み出して、該データに基づきＣＴ両画像再構成を実
行する。即ち、第２図に示すようにＲ波のピーク値を中
心とする一定幅の投影データから画像■、■、■・・・
を再構成することになる。
位相入力部４より中心位相φが入力されると、対話装置
５を介してこの位相φがプロセッサ６に入力される。こ
の位相として例えばφ＝Ｏ即らＲ波のピーク位置が中心
位相として入力されると、前記プロセッサ６は、位相φ
に相当する時刻の前記Ｘ線管の回転角を中心に、一枚の
ＣＴ両画像構成するのに必要な回転角θ分の投影データ
を、前記ディスク３の内の一連の長い投影データの中か
ら読み出して、該データに基づきＣＴ両画像再構成を実
行する。即ち、第２図に示すようにＲ波のピーク値を中
心とする一定幅の投影データから画像■、■、■・・・
を再構成することになる。
尚、一枚のＣＴ両画像再構成するのに必要な回転角θと
は、第３図に示すように少なくとも１８０度＋Ｘ線ファ
ンビームの拡がり角である。
は、第３図に示すように少なくとも１８０度＋Ｘ線ファ
ンビームの拡がり角である。
このように本実施例によれば、各画像は特定の位相を中
心とした投影データで構成されているため、各画像間で
の位相ズレがなく、心臓のような周期的に動く部位でも
一連の経時的画像間で被写体が動かないように見えるた
め、経時的な変化の観測を正しく行うことができ、特に
関心領域を数ケ所に設定して各ＲＯＩ内のＣＴ値の変化
を追跡するダイナミックスタデイを実行する場合には、
画像が動かないためＲＯＩ内のＣＴ値値化化正確に追跡
することができる。
心とした投影データで構成されているため、各画像間で
の位相ズレがなく、心臓のような周期的に動く部位でも
一連の経時的画像間で被写体が動かないように見えるた
め、経時的な変化の観測を正しく行うことができ、特に
関心領域を数ケ所に設定して各ＲＯＩ内のＣＴ値の変化
を追跡するダイナミックスタデイを実行する場合には、
画像が動かないためＲＯＩ内のＣＴ値値化化正確に追跡
することができる。
尚、上記実施例において、心電位相入力部４で２種類の
中心位相を入力するようにしてもよい。
中心位相を入力するようにしてもよい。
第４図は、第１の中心位相φ１＝Ｏ９第２中心位相φ２
＝０．５とした例でおる。これらの画像列をｉ１察する
とぎは、第１の中心位相φ１に対するものは心臓の収縮
期中心のものを一連として観察し、第２の中心位相φ２
に対するものは心臓の拡張期中心のものを一連として観
察することができる。
＝０．５とした例でおる。これらの画像列をｉ１察する
とぎは、第１の中心位相φ１に対するものは心臓の収縮
期中心のものを一連として観察し、第２の中心位相φ２
に対するものは心臓の拡張期中心のものを一連として観
察することができる。
く第２実施例〉
第５図に本発明の第２実施例装置の構成を示す。
同図に示す構成が第１図に示す構成と相違する点は、心
電位相入力部７が２種類の特定位相φ１゜φ２の間の投
影データより一枚のＣＴ両画像再構成するように、前記
２種類の特定位相を入力するように構成した点である。
電位相入力部７が２種類の特定位相φ１゜φ２の間の投
影データより一枚のＣＴ両画像再構成するように、前記
２種類の特定位相を入力するように構成した点である。
この第２実施例によれば、前記プロヒッサ６は位相φ１
から位相φ２にかけて収集された投影データを前記ディ
スク３より読み出し、これから一枚のＣＴ両画像再構成
することになる。
から位相φ２にかけて収集された投影データを前記ディ
スク３より読み出し、これから一枚のＣＴ両画像再構成
することになる。
第６図は、φ１＝Ｏ２φ２＝１、即ち、位相の上ではφ
１＝φ２とした例である。一方、第７図は、φ１＝１．
φ２＝１．２とした例である。このようにしても、結局
第１実施例と同様に特定の中心位相を中心とした所定幅
の投影データより画像を再構成することができ、第１実
施例と同様に各画像は特定の位相を中心とした投影デー
タで構成されているため、各画像間での位相ズレがなく
、心臓のような周期的に動く部位でも一連の経時的画像
間で被写体が動かないように見えるため、経時的な変化
の観測を正しく行うことができる。
１＝φ２とした例である。一方、第７図は、φ１＝１．
φ２＝１．２とした例である。このようにしても、結局
第１実施例と同様に特定の中心位相を中心とした所定幅
の投影データより画像を再構成することができ、第１実
施例と同様に各画像は特定の位相を中心とした投影デー
タで構成されているため、各画像間での位相ズレがなく
、心臓のような周期的に動く部位でも一連の経時的画像
間で被写体が動かないように見えるため、経時的な変化
の観測を正しく行うことができる。
尚、φ１＝Ｏとした時、φ２く１でも勿論よい。
ここで、この第２実施例では、φ１−φ２の値により、
かつ、スキャン時間（＝架台の１回転に要する時間）と
の関係から、一枚の画像を再構成するのに用いられる投
影データは、架台１回転分に較べて増減する。これは、
下記の手法のいずれかを講じることにより解決すること
ができる。
かつ、スキャン時間（＝架台の１回転に要する時間）と
の関係から、一枚の画像を再構成するのに用いられる投
影データは、架台１回転分に較べて増減する。これは、
下記の手法のいずれかを講じることにより解決すること
ができる。
（１）投影データ量が架台１回転分より多い時■第８図
に示すように回転角αより回転角２π十β（αくβ）ま
でのデータＰ（φ）のうちπ＋（α十β）／２を中心と
する２π分の投影データを用いればよい。
に示すように回転角αより回転角２π十β（αくβ）ま
でのデータＰ（φ）のうちπ＋（α十β）／２を中心と
する２π分の投影データを用いればよい。
■次のように重複している部分に関し両者の荷重平均を
とればよい。
とればよい。
Ｐ（φ）＝
・・・αくφくβ
ここで、上記の係数には第９図に示すようにＯ投影デー
タ量が１８０度±Ｘ線ファンビーム拡がり角以上あれば
ＣＴ両画像再構成可能である。
タ量が１８０度±Ｘ線ファンビーム拡がり角以上あれば
ＣＴ両画像再構成可能である。
例えば、Ｄ、Ｄａｒｋｅｒ ”０ｐｔｉｃａｌ ５ｈｏ
ｒｔ ５ｃａｎ Ｃｏｎ−ｖｏｌｕｔｉｏｎ Ｒｅｃｏ
ｎｓｔｒｕｃｔｉｏｎ ”）Ｉｅｄｉｃａｌ Ｐｈｙｓ
ｉｃｓ９（２）、 Ｍａｒｃｈ／Ａｐｒｉ １１９８２
の文献に正しく画像再構成できることが記載されている
。従って、例えば（φ１＋φ２）／２を中心にして必要
量だけ（１／２回転分＋Ｘ線ファンビーム角）投影デー
タを切り出すか、あるいは、上記文献に記載されている
骨子は、１８０度反対方向から両方のデータが収集され
ているときは、その総和を１とするような重み付けをし
て通常のコンボリューション／バックプロジェクション
により画像再構成をするものであるから、１／２回転分
＋Ｘ線ファンビーム角を越えるデータについても同様に
重み付けを行えばよい。
ｒｔ ５ｃａｎ Ｃｏｎ−ｖｏｌｕｔｉｏｎ Ｒｅｃｏ
ｎｓｔｒｕｃｔｉｏｎ ”）Ｉｅｄｉｃａｌ Ｐｈｙｓ
ｉｃｓ９（２）、 Ｍａｒｃｈ／Ａｐｒｉ １１９８２
の文献に正しく画像再構成できることが記載されている
。従って、例えば（φ１＋φ２）／２を中心にして必要
量だけ（１／２回転分＋Ｘ線ファンビーム角）投影デー
タを切り出すか、あるいは、上記文献に記載されている
骨子は、１８０度反対方向から両方のデータが収集され
ているときは、その総和を１とするような重み付けをし
て通常のコンボリューション／バックプロジェクション
により画像再構成をするものであるから、１／２回転分
＋Ｘ線ファンビーム角を越えるデータについても同様に
重み付けを行えばよい。
尚、本発明は上記実施例に限定されるものではなく本発
明の要旨の範囲内で種々の変形実施が可能である。
明の要旨の範囲内で種々の変形実施が可能である。
周期的な変動としては上記実施例のような心拍によるも
のに限らず、例えば呼吸によって周期的に変動する部位
についても同様に適用できる。この場合、呼吸波形を記
録して位相開明をとればよい。
のに限らず、例えば呼吸によって周期的に変動する部位
についても同様に適用できる。この場合、呼吸波形を記
録して位相開明をとればよい。
また、本発明はＸ線管の回転方向が１スキヤン毎に反転
するものにも適用することができる。これは、Ｘ線管、
検出器共に回転する第３世代のＸ線ＣＴ装置であり、検
出器の回転に伴なうケーブルさばきのために上記のよう
にスキャン方向を反転しているが、この方式では第１０
図に示すように投影データが欠落する時間が生じ、かつ
、高速スキャンも望めないが、同図に示すように１回の
Ｘ線管回転中に２回のＲ波ピークを検出しても、画像を
作るに足るθの回転角幅に亘って２回分の投影データが
得られない場合には、片方の画像のみを再構成するよう
にすればよい。この例でも、心拍による心臓の画像の歪
み状態は各画像間で略々同じであるので、例えば心筋の
ダイナミックスタデイには充分事足りる。
するものにも適用することができる。これは、Ｘ線管、
検出器共に回転する第３世代のＸ線ＣＴ装置であり、検
出器の回転に伴なうケーブルさばきのために上記のよう
にスキャン方向を反転しているが、この方式では第１０
図に示すように投影データが欠落する時間が生じ、かつ
、高速スキャンも望めないが、同図に示すように１回の
Ｘ線管回転中に２回のＲ波ピークを検出しても、画像を
作るに足るθの回転角幅に亘って２回分の投影データが
得られない場合には、片方の画像のみを再構成するよう
にすればよい。この例でも、心拍による心臓の画像の歪
み状態は各画像間で略々同じであるので、例えば心筋の
ダイナミックスタデイには充分事足りる。
また、連続スキャンを行うものでなくても、例えばＲａ
ｃｌｉｏｌｇｙ ｏｆ ５ｋｕｌｌ ａｎｄ Ｂｒａｉ
ｎ Ｖｏ１５ ｐｐ４３６６乃至ｐｐ４３６９に記載さ
れているような電子ビームスキャンでも′ｊ＠音の少な
い高画質を１辱るために高速スキャンを繰り返し、それ
の加算平均された投影データをもって画像を作ることに
なるが、このような時でも変動周期の特定位相を中心と
した所定幅の投影データで画像を作ることは有用でおる
（繰り返しスキャンにより１画像のスキャン時間は実質
的に伸びてしまうので、本発明の概念が有用になる）。
ｃｌｉｏｌｇｙ ｏｆ ５ｋｕｌｌ ａｎｄ Ｂｒａｉ
ｎ Ｖｏ１５ ｐｐ４３６６乃至ｐｐ４３６９に記載さ
れているような電子ビームスキャンでも′ｊ＠音の少な
い高画質を１辱るために高速スキャンを繰り返し、それ
の加算平均された投影データをもって画像を作ることに
なるが、このような時でも変動周期の特定位相を中心と
した所定幅の投影データで画像を作ることは有用でおる
（繰り返しスキャンにより１画像のスキャン時間は実質
的に伸びてしまうので、本発明の概念が有用になる）。
また、ＭＲＩにおいても、適切な造影剤が開発され高速
イメージングが実用化されてダイナミックスキャンを行
う時は、本発明の概念を適用することができる。
イメージングが実用化されてダイナミックスキャンを行
う時は、本発明の概念を適用することができる。
［発明の効果］
以上詳述したように本発明によれば、動きのある部位に
ついて位相的に一致した複数枚の同門画像を再構成する
ことができ、特にダイナミックスタデイに最適な同期画
像再構成装置を提供することができる。
ついて位相的に一致した複数枚の同門画像を再構成する
ことができ、特にダイナミックスタデイに最適な同期画
像再構成装置を提供することができる。
第１図は本発明の第１実施例装置のブロツ図、第２図は
第１実施例装置の動作を説明するめの図、第３図は１画
像を再構成するのに必要Ｘ線管の回転角を示す説明図、
第４図は中心位を２種類設定した第１実施例装置の動作
説明図第５図は本発明の第２実施例装置のブロック図第
６図、第７図は第２実施例装置の動作を説明る説明図、
第８図は１回転以上の回転角に亘つ投影データを収集し
た場合を説明する説明図、９図は荷重平均の係数ｋを示
す特性図、第１０はスキャン毎に回転方向を反転して投
影データ収集した場合の画像再構成を示す説明図、第１
図は従来の画像再構成の一例を示す説明図で父１・・・
スキャナ、２・・・信号検出手段、３・・・記憶４．７
・・・位相設定部、６−・・画像再構成部。 代理人 弁理士 三 澤 正 ぢり た な 相 １． １． １す て 第 １図 を 〉る。 （部、 弔１図 畳 第３図
第１実施例装置の動作を説明するめの図、第３図は１画
像を再構成するのに必要Ｘ線管の回転角を示す説明図、
第４図は中心位を２種類設定した第１実施例装置の動作
説明図第５図は本発明の第２実施例装置のブロック図第
６図、第７図は第２実施例装置の動作を説明る説明図、
第８図は１回転以上の回転角に亘つ投影データを収集し
た場合を説明する説明図、９図は荷重平均の係数ｋを示
す特性図、第１０はスキャン毎に回転方向を反転して投
影データ収集した場合の画像再構成を示す説明図、第１
図は従来の画像再構成の一例を示す説明図で父１・・・
スキャナ、２・・・信号検出手段、３・・・記憶４．７
・・・位相設定部、６−・・画像再構成部。 代理人 弁理士 三 澤 正 ぢり た な 相 １． １． １す て 第 １図 を 〉る。 （部、 弔１図 畳 第３図
Claims (3)
- （１）所定の周期で変動する部位について、複数周期に
亘って前記部位をスキャンして投影データを収集するス
キャナと、前記データ収集中に亘って前記部位の周期的
な動きと対応する信号を検出する信号検出手段と、前記
スキャナからのデータ及び前記信号検出手段からの信号
を関連付けて記憶する記憶部と、前記変動周期中の特定
位相を設定する位相設定部と、設定された特定位相に対
応するデータを含んで少なくとも一枚の画像を再構成す
るのに必要な投影データを各画像毎に前記記憶部より読
み出し、時間軸方向に複数枚の同期画像を再構成する画
像再構成部とを有することを特徴とする同期画像再構成
装置。 - （２）位相設定部は、前記記憶部より読み出される投影
データの収集時間幅の中心に一致する中心位相を設定す
るものである特許請求の範囲第１項記載の同期画像再構
成装置。 - （３）位相設定部は、前記記憶部より読み出される投影
データの収集開始時間に一致する第１の特定位相と、収
集終了時間に一致する第２の特定位相とを設定する特許
請求の範囲第１項記載の同期画像再構成装置。
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP61189173A JP2529949B2 (ja) | 1986-08-12 | 1986-08-12 | 同期画像再構成装置 |
| DE3751690T DE3751690T2 (de) | 1986-08-12 | 1987-08-12 | Anordnung zur dynamischen Abtastung während periodischer Verstellung |
| US07/084,170 US4868747A (en) | 1986-08-12 | 1987-08-12 | Method and system for dynamic computed tomography scanning at cyclic displacement points |
| EP87307135A EP0257922B1 (en) | 1986-08-12 | 1987-08-12 | System for dynamic scan at cyclic displacement |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP61189173A JP2529949B2 (ja) | 1986-08-12 | 1986-08-12 | 同期画像再構成装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS6346138A true JPS6346138A (ja) | 1988-02-27 |
| JP2529949B2 JP2529949B2 (ja) | 1996-09-04 |
Family
ID=16236703
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP61189173A Expired - Lifetime JP2529949B2 (ja) | 1986-08-12 | 1986-08-12 | 同期画像再構成装置 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US4868747A (ja) |
| EP (1) | EP0257922B1 (ja) |
| JP (1) | JP2529949B2 (ja) |
| DE (1) | DE3751690T2 (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS63214247A (ja) * | 1986-12-29 | 1988-09-06 | ゼネラル・エレクトリック・カンパニイ | 像を発生する方法 |
| JP2002263089A (ja) * | 2001-03-06 | 2002-09-17 | Shimadzu Corp | デジタル式ｘ線撮影装置 |
Families Citing this family (37)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6488689A (en) * | 1987-09-29 | 1989-04-03 | Toshiba Corp | Three dimensional labelling device |
| JPH01265950A (ja) * | 1988-04-15 | 1989-10-24 | Toshiba Corp | スペクトロスコピックイメージング装置 |
| US5262946A (en) * | 1988-10-20 | 1993-11-16 | Picker International, Inc. | Dynamic volume scanning for CT scanners |
| US5396418A (en) * | 1988-10-20 | 1995-03-07 | Picker International, Inc. | Four dimensional spiral volume imaging using fast retrace |
| JPH0616099B2 (ja) * | 1989-02-07 | 1994-03-02 | 浜松ホトニクス株式会社 | Ｃｔ装置におけるデータ補正装置 |
| FR2656128B1 (fr) * | 1989-12-14 | 1994-07-29 | Gen Electric Cgr | Procede de reconstruction tridimensionnelle selective mettant en óoeuvre des acquisitions bidimensionnelles de radiologie a deux energies, notamment pour representer la structure osseuse du corps d'un patient. |
| JPH0447479A (ja) * | 1990-06-13 | 1992-02-17 | Toshiba Corp | 画像表示装置 |
| US5251128A (en) * | 1990-11-19 | 1993-10-05 | General Electric Company | Motion artifact reduction in projection imaging |
| GB9025431D0 (en) * | 1990-11-22 | 1991-01-09 | Advanced Tech Lab | Three dimensional ultrasonic imaging |
| US5212737A (en) * | 1991-09-30 | 1993-05-18 | Imatron, Inc. | Real time data averager for use in computed tomography scanning |
| WO1994004977A1 (de) * | 1992-08-11 | 1994-03-03 | Siemens Aktiengesellschaft | Bildfolgen-referenzspeicherverfahren zur störsignalunterdrückung und ausfilterung bei medizinischen bildfolgen |
| EP0587334B1 (en) * | 1992-09-09 | 1999-06-30 | Picker International, Inc. | Imaging methods and apparatus |
| EP0615135B1 (de) * | 1993-03-06 | 1998-09-09 | Philips Patentverwaltung GmbH | MR-Verfahren zur zwei- oder dreidimensionalen Abbildung eines Untersuchungsbereichs und Anordnung zur Durchführung des Verfahrens |
| US5600701A (en) * | 1993-05-11 | 1997-02-04 | Hitachi Medical Corporation | X-ray imaging system and method therefor |
| FR2706043B1 (fr) * | 1993-06-02 | 1995-07-07 | Commissariat Energie Atomique | Installation et procédé de reconstruction d'images tridimensionnelles. |
| US5526442A (en) * | 1993-10-04 | 1996-06-11 | Hitachi Medical Corporation | X-ray radiography method and system |
| DE19740214A1 (de) * | 1997-09-12 | 1999-04-01 | Siemens Ag | Computertomograph |
| US6298111B1 (en) * | 1998-06-04 | 2001-10-02 | Kabushiki Kaisha Toshiba | X-ray computed tomography apparatus |
| US6236705B1 (en) * | 1998-06-17 | 2001-05-22 | Her Majesty The Queen In Right Of Canada, As Represented By The Minister Of National Defence | Method for tracing organ motion and removing artifacts for computed tomography imaging systems |
| JP4698780B2 (ja) * | 1998-09-15 | 2011-06-08 | シーメンス アクチエンゲゼルシヤフト | 像再構成方法及び測定データ取得方法 |
| US6233478B1 (en) * | 1998-09-28 | 2001-05-15 | Advanced Research & Technology Institute | Apparatus and method for constructing computed tomography image slices of an object undergoing cyclic motion |
| FI104042B (fi) * | 1998-09-29 | 1999-11-15 | Aaro Kiuru | Menetelmä keuhkojen perfuusion mittaamiseksi |
| IL137251A0 (en) * | 1999-07-28 | 2001-07-24 | Gen Electric | Retrospective cardiac gating with cine scans on a multislice scanner |
| US6639965B1 (en) * | 1999-09-30 | 2003-10-28 | General Electric Company | Methods and apparatus for cardiac imaging with conventional computed tomography |
| US6252924B1 (en) | 1999-09-30 | 2001-06-26 | General Electric Company | Method and apparatus for motion-free cardiac CT imaging |
| DE19957083B4 (de) | 1999-11-28 | 2004-11-18 | Siemens Ag | Verfahren zur Untersuchung eines eine periodische Bewegung ausführenden Körperbereichs |
| DE60138118D1 (de) | 2000-09-29 | 2009-05-07 | Ge Med Sys Global Tech Co Llc | Spiralabtastrekonstruktion des herzens für tomograph mit zweidimensionaler detektoranordung |
| US7697744B2 (en) * | 2001-03-06 | 2010-04-13 | Kabushiki Kaisha Toshiba | X-ray diagnostic apparatus and image processor |
| US6708052B1 (en) * | 2001-04-11 | 2004-03-16 | Harbor Ucla Research And Education Institute | Method and apparatus for cardiac imaging with minimized cardiac motion artifact |
| US20060247550A1 (en) * | 2001-05-28 | 2006-11-02 | Arvind Thiagarajan | System and methods for gating heart signals |
| US6438196B1 (en) | 2001-06-28 | 2002-08-20 | General Electric Company | EKG driven CT image reconstruction for cardiac imaging |
| JP2003141548A (ja) * | 2001-10-22 | 2003-05-16 | Ge Medical Systems Global Technology Co Llc | ３次元ラベリング装置及びその方法 |
| US6526117B1 (en) * | 2001-11-09 | 2003-02-25 | Ge Medical Systems Global Technology Company, Llc | Method and apparatus to minimize phase misregistration artifacts in gated CT images |
| US7477928B2 (en) | 2002-05-17 | 2009-01-13 | Ge Medical Systems Global Technology Company, Llc | Method and system for associating an EKG waveform with a CT image |
| SE526438C2 (sv) * | 2003-04-01 | 2005-09-13 | Sectra Imtec Ab | Metod och system för mätning i en dynamisk sekvens av medicinska bilder |
| JP5672147B2 (ja) * | 2011-05-24 | 2015-02-18 | コニカミノルタ株式会社 | 胸部診断支援情報生成システム |
| US8903149B1 (en) * | 2013-01-18 | 2014-12-02 | Anritsu Company | System and method of communicating information about an object concealed by a scanned surface |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS61100233A (ja) * | 1984-10-23 | 1986-05-19 | 株式会社東芝 | Ｘ線ｃｔスキヤナ |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4037585A (en) * | 1975-10-10 | 1977-07-26 | Gildenberg Philip L | Method of measuring blood flow by radiation |
| US4182311A (en) * | 1977-04-22 | 1980-01-08 | Varian Associates, Inc. | Method and system for cardiac computed tomography |
| FR2455451A1 (fr) * | 1979-05-02 | 1980-11-28 | Radiologie Cie Gle | Procede de tomodensitometrie et tomodensitometre adapte a ce procede |
| NL8003354A (nl) * | 1980-06-09 | 1982-01-04 | Philips Nv | Stralingsonderzoekapparaat met beeldsubtractie. |
| IL69327A (en) * | 1983-07-26 | 1986-11-30 | Elscint Ltd | Automatic misregistration correction |
| US4703424A (en) * | 1985-01-03 | 1987-10-27 | General Electric Company | Method and means of correcting for a shift in the center of rotation of a rotating fan beam CT system |
| US4729099A (en) * | 1985-07-19 | 1988-03-01 | Picker International, Inc. | Therapy treatment planning by radiation dose determination |
-
1986
- 1986-08-12 JP JP61189173A patent/JP2529949B2/ja not_active Expired - Lifetime
-
1987
- 1987-08-12 EP EP87307135A patent/EP0257922B1/en not_active Expired - Lifetime
- 1987-08-12 DE DE3751690T patent/DE3751690T2/de not_active Expired - Lifetime
- 1987-08-12 US US07/084,170 patent/US4868747A/en not_active Expired - Lifetime
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS61100233A (ja) * | 1984-10-23 | 1986-05-19 | 株式会社東芝 | Ｘ線ｃｔスキヤナ |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS63214247A (ja) * | 1986-12-29 | 1988-09-06 | ゼネラル・エレクトリック・カンパニイ | 像を発生する方法 |
| JP2002263089A (ja) * | 2001-03-06 | 2002-09-17 | Shimadzu Corp | デジタル式ｘ線撮影装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| DE3751690T2 (de) | 1996-08-08 |
| EP0257922A2 (en) | 1988-03-02 |
| EP0257922A3 (en) | 1988-06-15 |
| DE3751690D1 (de) | 1996-03-14 |
| US4868747A (en) | 1989-09-19 |
| EP0257922B1 (en) | 1996-01-31 |
| JP2529949B2 (ja) | 1996-09-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPS6346138A (ja) | 同期画像再構成装置 | |
| US6539074B1 (en) | Reconstruction of multislice tomographic images from four-dimensional data | |
| JP4340533B2 (ja) | コンピュータ断層撮影 | |
| US7085342B2 (en) | Method for tracking motion phase of an object for correcting organ motion artifacts in X-ray CT systems | |
| US6522712B1 (en) | Reconstruction of computed tomographic images using interpolation between projection views | |
| US6269140B1 (en) | X-ray computerized tomography apparatus control method therefor and image generating method using the apparatus | |
| US20070197907A1 (en) | Method and tomography unit for carrying out an analysis of a movement of an object | |
| JP4429677B2 (ja) | 周期的に運動する器官のｃｔ画像形成方法 | |
| JPH053867A (ja) | 三次元画像診断装置 | |
| JP2001190544A (ja) | コンピュータ断層撮影システムにより取得される画像データから再構成される画像内のアーティファクトを減少させる方法及び装置 | |
| US20100111393A1 (en) | X-ray ct apparatus and control method of x-ray ct apparatus | |
| EP1277175B1 (en) | Reconstruction of computed tomographic images using interpolation between projection views | |
| US6959067B2 (en) | X-ray imaging method and a 3D-rotational X-ray apparatus for applying this method | |
| JP3819283B2 (ja) | Ｘ線ｃｔ装置 | |
| Pack et al. | Four‐dimensional computed tomography of the left ventricle, Part I: motion artifact reduction | |
| WO2004071301A1 (ja) | X線ct装置 | |
| US20050058238A1 (en) | Method for examining a body region of an examination object, which body region performs a periodic movement, and CT apparatus for carrying out such a method | |
| JPH02243134A (ja) | 磁気共鳴イメージング装置 | |
| CN101352352B (zh) | X射线ct装置 | |
| JPH10328175A (ja) | Ｘ線ｃｔ装置 | |
| JPH03143425A (ja) | 生体データ収集装置 | |
| JP2020014552A (ja) | 磁気共鳴イメージング装置 | |
| JP3421025B2 (ja) | Ｘ線ｃｔ装置 | |
| JP4257127B2 (ja) | 磁気共鳴撮像方法および磁気共鳴撮像装置 | |
| JP2001353143A (ja) | 断層画像生成方法、ｘ線ｃｔ画像生成方法およびｘ線ｃｔ装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| EXPY | Cancellation because of completion of term |