JP5031941B2 - コンピュータ断層撮影システム - Google Patents
コンピュータ断層撮影システム Download PDFInfo
- Publication number
- JP5031941B2 JP5031941B2 JP2000327188A JP2000327188A JP5031941B2 JP 5031941 B2 JP5031941 B2 JP 5031941B2 JP 2000327188 A JP2000327188 A JP 2000327188A JP 2000327188 A JP2000327188 A JP 2000327188A JP 5031941 B2 JP5031941 B2 JP 5031941B2
- Authority
- JP
- Japan
- Prior art keywords
- data
- slice
- image
- helical
- hhs
- 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.)
- Expired - Fee Related
Links
- 238000002591 computed tomography Methods 0.000 title claims description 50
- 238000013480 data collection Methods 0.000 claims description 11
- 238000001514 detection method Methods 0.000 claims description 9
- 238000000034 method Methods 0.000 description 56
- 230000002123 temporal effect Effects 0.000 description 32
- 239000011295 pitch Substances 0.000 description 31
- 238000012545 processing Methods 0.000 description 30
- 238000013213 extrapolation Methods 0.000 description 18
- 238000010586 diagram Methods 0.000 description 16
- 238000003384 imaging method Methods 0.000 description 14
- 230000008569 process Effects 0.000 description 12
- 230000035945 sensitivity Effects 0.000 description 10
- 238000004422 calculation algorithm Methods 0.000 description 8
- 230000000747 cardiac effect Effects 0.000 description 8
- 238000005070 sampling Methods 0.000 description 8
- 230000005526 G1 to G0 transition Effects 0.000 description 7
- 238000011156 evaluation Methods 0.000 description 7
- 238000006073 displacement reaction Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 230000004044 response Effects 0.000 description 5
- 238000003491 array Methods 0.000 description 4
- 238000001914 filtration Methods 0.000 description 4
- 210000004072 lung Anatomy 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 230000015556 catabolic process Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 239000000470 constituent Substances 0.000 description 3
- 238000006731 degradation reaction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 210000003437 trachea Anatomy 0.000 description 3
- 230000002792 vascular Effects 0.000 description 3
- 210000000709 aorta Anatomy 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 210000005240 left ventricle Anatomy 0.000 description 2
- 210000004165 myocardium Anatomy 0.000 description 2
- 238000011548 physical evaluation Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 210000000596 ventricular septum Anatomy 0.000 description 2
- 0 *=*CC(*1)*1N Chemical compound *=*CC(*1)*1N 0.000 description 1
- 206010003211 Arteriosclerosis coronary artery Diseases 0.000 description 1
- 241001669679 Eleotris Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 210000004351 coronary vessel Anatomy 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 230000003205 diastolic effect Effects 0.000 description 1
- 238000010894 electron beam technology Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000002107 myocardial effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000002685 pulmonary effect Effects 0.000 description 1
- 230000010349 pulsation Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 235000019801 trisodium phosphate Nutrition 0.000 description 1
Images
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/027—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis characterised by the use of a particular data acquisition trajectory, e.g. helical or spiral
-
- 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]
-
- 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/542—Control of apparatus or devices for radiation diagnosis involving control of exposure
-
- 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/58—Testing, adjusting or calibrating thereof
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
-
- 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)
- Optics & Photonics (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Public Health (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Physiology (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Description
【発明の属する技術分野】
本発明は、ヘリカル(螺旋形)のコンピューター断層撮影スキャンに係り、特に改善された時間的解像度を有するヘリカル・コンピュータ断層撮影スキャンに関する。なお、時間的解像度は、患者から特別な信号を得る必要がないマルチスライス・ヘリカル・スキャンの画像再構成アルゴリズムを使用することによって、各々のスライス・イメージにおいて、改善可能である。
【0002】
【従来の技術】
ヘリカル・コンピュータ断層撮影(CT)スキャンは、例えば、米国特許第4,630,202号に記載されている。ヘリカルCTスキャンは、米国特許第4,965,726号において開示されるように、マルチスライス位置でデータを収集するマルチロウ検出器アレイを取り入れることができる。マルチスライス・ヘリカルCTは、そのようなアレイを使用する。これによって、長手方向の空間解像度を改善し、走査長さを延長し、さらに走査時間を短縮することが可能となる。その結果、時間的解像度とコントラスト解像度の両方が大量に改善される。ヘリカル・スキャンが同じ量に対して繰り返される場合、量の時間的変化を観測することが可能である。
【0003】
多様な画像再構成アルゴリズムが、マルチスライス・ヘリカルCT、例えば「Algorithm for Image Reconstruction in Multi-Slice Helical CT,Taguchi et al.,Med.Phys.25:550-561(1998)」などで使用される。このアルゴリズムは、(1)ヘリカルピッチの慎重な選択が要求される最適化されたサンプリング・スキャン、(2)長手方向に関するヘリカル・フィルタ内挿補間(HFI)、および(3)扇形ビームフィルタバックプロジェクションによる再構成(通常の扇形ビーム再構成技術は、第3の工程において、適用可能である)の3つの部分からなる走査技術を含んでいる。HFIは、180°、または、360°分離された、いずれかのデータ(すなわち線形内挿(または外挿)中の異なる時間に得られるデータ)を使用して、スライス位置でのデータを得るために長手方向のフィルタを使用する。
【0004】
中心光線での長手方向のデータサンプリングパターンは、そのようなヘリカル・ピッチではまばらとなり、追加の技術がなければ画像の品位を劣化させてしまうため、最適化されたサンプリング走査において、一部の特定のヘリカル・ピッチ(すなわち偶数の整数)は、避けなければならない。
【0005】
ヘリカル内挿(HFI)は、以下の2つの工程からなる。1)フィルタ幅内の位置で複数のリサンプリングデータを得るために、2点線形内挿(または外挿)を繰り返す、2)リサンプリングデータをフィルタする。したがって、フィルタ幅がゼロに等しいとき、データサンプリング位置によって、HFIは内挿または外挿となる。フィルタ幅は、長手領域と実効的なスライス厚との比率として定義される。
【0006】
マルチスライス・ヘリカルCTを用いて心臓の撮像を行なって、実際的な時間的解像度を改良する特別な方法が開発されたが、欠点がいくつかある。例えば、同じ心臓位相に対するデータを再分類するためにEKG信号(心電図)を必要とし、周期的に動く被検体にのみ適用可能である。
【0007】
平行ビーム形状において、画像再構成のために必要な最小角度は、180°である。扇形ビーム形状においては、180°+扇形角度を走査する投影データが、そのようなデータを含む。しかしながら、一部の投影は投影データにおいて2度計測され、追加の処理がなければアーチファクトを発生させてしまう場合がある。そのような冗長の影響を除去し、データ打切りを避けるために、扇形ビーム画像再構成の前に投影データ組に平滑化重み付けが行なわれる。ハーフスキャン(HS)と呼ばれるこのアルゴリズムは、Computed Tomography Scanning with Simultaneous Patient Translation,C.Crawford等、Med.Phys.,17:967-982(1990)中に説明されているように、CT画像の時間的解像度を改善するために使用されてきた。
【0008】
その他、Computed Tomography Scanning with Simultaneous Patient Translation,C.Crawford等、Med.Phys.,17:967-982(1990)に記載されているアンダー・スキャニング(US)およびHigh Temporal Resolution Reconstruction for Reducing Motion Artifact Causedby Cardiac Motion,Shenetal,Jap.JofRad.Tech.,54(11):1287-1294(1998)に記載の高時間的解像度再構成(HTRR)という方法が使用されてきた。HS、USおよびHTRRの中の違いは、図1に示すように、単に重み付け関数である。修正USと呼ばれる他の方法は、捕捉タイミングを考慮した「直接」データ(direct data)と、「対向」データとの間の重み付け技術を指す。それを図2に示す。
【0009】
汎用CTの電子ビームCT(EBCT)およびEKGゲート制御再構成という他の技術は、時間的解像度を改善するために使用される。EBCTは、最短走査時間0.05秒で最良の時間的解像度を実現する。しかし、EBCTは、不十分な長手方向の空間的解像度、システム・コストおよび限られたX線露光という欠点を有している。このように、マルチスライス・ヘリカルCTを用いた心臓の撮像は、そのような問題がなく、心臓撮像において有効となる可能性がある。
【0010】
最近、心臓の周期的位相にしたがってヘリカル・データを分類する、小ヘリカル・ピッチを用いたEKGゲート制御再構成法が汎用CTのために開発され、より良い時間的解像度を得ることにおいての可能性を示した。この方法において、ヘリカル・マルチスライス・ピッチが1.0の場合、FWTMでまたはFWTAで0.06秒が達成可能である。これらの専用の心臓再構成法の最も大きい利点は、全心臓周期の間の動的体積心臓撮像が可能なことである。空間的解像度および画像ノイズは、それらの心臓位相が一致した場合にのみ、180°または360°分離されたデータを使用して改善可能である。この方法の2つの最も大きな問題は、ヘリカルピッチが小さい場合に高線量検査となる場合があり、コントラスト因子の調査の際に問題が発生する可能性がある。それらが同じ心臓位相であるが異なる時間に得られたデータを使用するため、コントラスト濃度など、心臓の運動以外の変化が画像の品位を劣化させる場合がある。
【0011】
上述のような方法で得た個々の画像の時間的解像度は、心臓および隣接した肺の気管などの高速運動を行なう器官には不十分である。ヘリカルCTおよびマルチスライス・ヘリカルCTにおいて、より良い時間的解像度が必要である。
【0012】
【発明が解決しようとする課題】
本発明の目的は、コンピュータ断層撮影(CT)の時間的解像度を改善する方法およびシステムを提供することである。
【0013】
【課題を解決するための手段】
本発明は、X線ビーム発生源と、前記X線ビーム発生源から被検体を介して到達するX線を検出するための複数の検出素子がチャンネル方向及びスライス方向に配列されてなる2次元検出器とを有するガントリと、前記被検体を載置する寝台と、前記寝台と前記ガントリとの少なくとも一方を移動する移動ユニットと、前記2次元検出器に接続され、前記2次元検出器によって検出されたX線のデータを収集するデータ収集ユニットと、前記データ収集ユニットに接続され、前記データ収集ユニットによって収集されたデータに対して前記スライス方向に関して内挿または外挿補間処理することにより所定のスライス位置に対応する投影データを生成し、前記投影データに対して時間方向に関して異なる複数の重み付け関数により重み付けをすることにより複数の投影データ組を生成する重み付け回路と、前記複数のデータ組から、スライス位置が同じであるがデータ収集時刻が異なる複数の断層画像をそれぞれ再構成する断層画像再構成ユニットとを具備することを特徴とするコンピュータ断層撮影システムである。
【0014】
さらに、本発明は、前記同時に得られた前記データの一部のみを処理することを特徴とする。
【0015】
さらに、本発明は、前記データを複数の組に分け、このデータの組数を加重平均処理により減らすことを特徴とする。
【0016】
さらに、本発明は、前記データの処理において、線形内挿と外挿の少なくとも1つを使用することを含む。
【0017】
さらに、本発明は、前記データの処理において、非線形内挿と外挿の少なくとも一つを使用することを含む。
【0018】
さらに、本発明は、前記処理データに対応する対向データを得、前記データ収集タイミングに基づいて前記処理データおよび前記対向データを重み付け加算し、前記重み付け加算したデータに基づいて前記画像を再構成する。
【0019】
さらに、本発明は、前記処理データの重み付けを時間的に変位する。
【0020】
さらに、本発明は、前記重み付けは、γm が最大チャネル角度である場合、π+2γm の範囲を除くゼロの値を有する重み付け関数を使用する。
【0021】
さらに、本発明は、前記関数は、
βnew =β−βc +((π+2γm )/2)
βnew =β−βc +π
とのいずれか一方を含む。
【0022】
さらに、本発明は、前記被検体から得た信号に基づいて前記画像を再構成するタイミングを変位する。
【0023】
さらに、本発明は、前記被検体から得た信号に基づいて前記X線の発生を制御する。
【0024】
さらに、本発明は、重み付けは、2πの範囲を除くゼロの値を有する重み付け関数を使用する。
【0025】
さらに、本発明は、心臓画像に対応するデータを収集し、前記被検体からEKG信号を測定し、前記EKG信号に基づいて前記画像を再構成するタイミングを変位する。
【0026】
さらに、本発明は、前記被検体の心拡張位相にほぼ対応するデータを収集する。
【0027】
さらに、本発明は、複数個の列のうち選択された列のデータを処理する。
【0028】
さらに、本発明は、複数個の列に対してデータを同時に収集し、各部分に対して再構成されるスライスの部分に最も近い列からデータを選択し、前記処理データを生成するために、前記部分のそれぞれに対応する前記データを処理する。
【0029】
さらに、本発明は、各部分に対して再構成されるスライスの部分に最も近い列のうちの2つからデータを選択する。
【0030】
さらに、本発明は、前記データを内挿補間し、前記内挿データを重み付けする。
【0031】
さらに、本発明は、前記データから内挿および外挿の少なくとも一方の補間を行なう。
【0032】
さらに、本発明は、補間データを得るために、前記データをフィルタし、重み付けは、前記補間データを重み付けする。
【0033】
本発明は、X線源と、前記X線源から照射され被検体を透過したX線を検出するための検出器と、前記検出器により検出されたデータを処理する処理回路と、前記処理回路で処理されたデータにデータ収集タイミングに応じて重み付けをする重み付け回路と、前記重み付けされたデータに基づいて画像を再構成する再構成回路とを具備することを特徴とする。
【0034】
さらに、本発明は、前記X線源を前記被検体から見て螺旋状に移動する機構を更に備える。
【0035】
さらに、本発明は、前記処理回路は、非線型内挿回路および外挿回路のうちの少なくとも一方を含む。
【0036】
さらに、本発明は、前記処理回路は、線型内挿回路および外挿回路のうちの少なくとも一方を含む。
【0037】
さらに、本発明は、前記処理回路は、フィルタ回路を含む。
【0038】
さらに、本発明は、前記重み付け回路は、重み付け関数を時間変位させる手段を有する。
【0039】
さらに、本発明は、前記重み付け回路は、重み付け関数を選択的に時間変位させる手段を含む。
【0040】
さらに、本発明は、前記被検体の生体信号を獲得する手段を含む。
【0041】
さらに、本発明は、前記生体信号を使用して前記被検体によって受けられた線量を減らす手段を含む。
【0042】
さらに、本発明は、前記生体信号を使用して前記X線源による照射を選択的に制御する手段を含む。
【0043】
さらに、本発明は、前記X線源によって射出された照射を視準するように配置された視準装置を含む。
【0044】
さらに、本発明は、重み付け関数を自動的に時間変位させる手段を含む。
【0045】
さらに、本発明は、前記被検体の心電波形に従って前記重み付け関数を時間変位させる。
【0046】
さらに、本発明は、前記検出器は、被検体のスライス方向に沿って複数の検出カラムを含み、複数の検出素子列は同時にX線を検出する。
【0049】
【発明の実施の形態】
以下、図面を参照して本発明による装置を好ましい実施形態により説明する。
本発明は、マルチスライス・ヘリカルCTに適用可能である。ここでは、その代表例として、4スライス・ヘリカルCTとして説明する。また、スライス厚は、等角点(回転軸)で計測されるステップ・アンド・シュート(天板の断続的な移動と同期した断続的なX線撮影)での走査中の各スライス画像の有効スライス厚として定義される。ヘリカル・スキャンの場合の幾何学的スライス厚は、シングルスライスCTスキャナの場合の再構成画像の有効スライス厚とは異なり、再構成アルゴリズムおよびパラメータの選択に応じて決まる。この適用において、ヘリカル・スライス・ピッチ(P)は、スライス厚に対する、1回転あたりの寝台移動量の比率として定義される。
【0050】
図3に本実施形態によるコンピュータ断層撮影システムを示し、図4に図3のガントリの斜視図を示し、図5に図3の2次元配列型検出器の斜視図である。ガントリ1は、X線をほぼ円錐形で照射するX線源3と、2次元に配列された複数個の検出素子5A、すなわち複数個の検出素子5Aの1次元配列が複数並列された2次元配列型X線検出器5とを収容する。図5では、複数個の検出素子5Aの1次元配列が4列あり、それらが並列されていることを示している。
【0051】
X線源3および2次元配列型X線検出器5は、寝台6の天板に載置されている被検体を挟んで互いに対向するように、回転リング2上に搭載される。2次元配列型X線検出器5は、回転リング2上に取付けられる。通常、各検出素子5Aは、1つのチャネルに対応している。X線源3からのX線は、X線フィルタ4を介して被検体に照射される。被検体を通過したX線は、2次元配列型X線検出器5によって、電気信号として検出される。
【0052】
X線制御装置8は、高電圧発生器7にトリガ信号を供給する。高電圧発生器7は、トリガ信号が受信されるタイミングで高電圧をX線源3に印加する。これによって、X線がX線源3から放射される。ガントリ/寝台制御装置9は、ガントリ1の回転リング2の回転と寝台6の天板の移動とを同期をとって制御する。システム制御装置10は、システム全体の制御中枢を構成し、被検体から見て、X線源3が、螺旋(ヘリカル)の軌道を描いて、いわゆるヘリカル・スキャンを実行するようにX線制御装置8およびガントリ/寝台制御装置9を制御する。具体的には、回転リング2は一定角速度で連続的に回転させられ、天板が一定速度で移動させられ、X線はX線源3から一定角間隔で断続的に又は連続的に照射される。
【0053】
2次元配列型X線検出器5の出力信号は、各チャネルごとにデータ収集ユニット11によって増幅され、投影データを生成するためにデジタル信号に変換される。データ収集装置11から出力される投影データは、信号処理ユニット12に供給される。信号処理ユニット12はその投影データを使用して画像を再構成する。3次元画像データまたはユニット12によって生成された断層撮影画像データは、表示装置14に送られ、そこで3次元画像として、または断層撮影画像として可視的に表示される。
【0054】
図6(a)に処理ユニット12の構成を示す。信号処理ユニット12は、投影データの全ての必要な内挿または外挿補間を実行する処理回路20を有する。データは、典型的には、長手方向に沿って補間される。処理回路20のより詳細な図を、図6(b)に示す。データは、入力端子26を介して入力され、その後、必要に応じて、回路23、24および25により、フィルタ、内挿補間または外挿補間される。処理データは、出力端子27を介して出力される。
【0055】
重み付け回路21は、処理回路20中で処理されるデータを決定して重み付けを行なう。再構成回路22は、重み付けされたデータを使用して画像を再構成する。回路の動作は、本発明による方法の説明と関連して後述する。
【0056】
図3のシステムの構成要素は、プログラムされたコンピュータとして実装されることが可能である。特に、システム制御装置10および処理ユニット12は個別のコンピュータであってもよく、または、2つの素子の機能は単一のコンピュータ上に搭載されてもよい。また、本発明はソフトウェアの形で実現可能であり、コンピュータ・プログラム製品を生成する光ディスクの磁気など、記録可能媒体上に格納されることが好ましい。
【0057】
HHS(ヘリカルハーフスキャン)と呼ばれる、本発明による第1の再構成方法をここで説明する。第1に、修正内挿補間(修正HFI)が投影データに対して施され、その後、その補間された投影データが重み付けされ、その重み付けされたデータに基づいて再構成が実行される。
【0058】
図7において、修正HFI方法では、30で示すスライス位置zs に対するデータの長手方向のフィルタにより、同時に収集された投影データどうしで補間が行われる。この補間は、内挿(線形または非線形)および/または外挿が使用される。なお、図7において、スライス位置は、実線で示す。斜めの太い黒線31は、スライス画像を再構成するために必要な一揃いのデータを得るために、内挿および/または外挿中で使用される実際に収集された投影データ(直接データ(実データ))を示している。
【0059】
z方向の内挿補間を決定するための(1)式を以下に示す。なお、“Direct”は、仮想データ(対向データ)ではなく、実際に収集された投影データ(直接データ(実データ))を表し、“Data For Slice”は、実データから生成された補間データを表している。
【0060】
【数1】
【0061】
ただし、zi は第i番目のビューでのX線源の位置である。
zj およびzk は、v番目のビュー、ビュー角度βおよびチャネル角度γで得られたヘリカルスライス位置に最も近いデータおよび次に最も近いデータのサンプリング位置をそれぞれ示している。
βi は、第i番目のビューでのビュー角度である。
pは、ヘリカル・スライス・ピッチである。
nvは、X線源の1回転あたりのビューの数である。
kvは、内挿/外挿/フィルタが可能な範囲である。
【0062】
同時に収集されたデータとは、同時に撮像され、検出されたデータを意味する。この例では、4列分のデータが同時に収集されて、図7の第1〜第4の列として示される。スライスは、この例では、スライス位置に2つの最も近い列において同時に得られるデータから内挿および/または外挿される。例えば、図7において、第1および第2の列は、第1の回転の後半のデータを内挿および/または外挿補間するために用いられる。第2の回転の前半のデータが第2列、第3列及び第4列のデータから内挿および/または外挿補間により得られる。同様にして、時間(X線源角度)ごとに、4列のうちスライス位置に最も近い2つの列のデータが選択される。なお、この方法において、180°まで到達したデータ、または異なる回転のデータのいずれも使用されない。
【0063】
次に、補間データは、データ収集時間に基づいて重み付けされる。重み付け曲線32を図7に示している。ただし、ステップ・アンド・シュートのHSのための重み付け関数と、USのための重み付け関数とのいずれも適用可能である。それらの2つの重み付け関数を、(2)式および(3)式として以下に示す。
【0064】
【数2】
【0065】
満足させる重みω(β,γ)が適用される(これらの式は、参考として本文献に引用されるCrawford等のAppendixCに示される)。
【0066】
【数3】
【0067】
以下の重み付けが、適用可能である(これらの式は、参考として本文献に引用されるCrawford等のAppendixBに示される)。
【0068】
ω(β,γ)+ω(β+π+2γ,−γ)=2
仮想データ(対向データ)が生成される場合、以下の(4)式に示すような任意の重み付け関数が適用可能である:
【数4】
【0069】
ただし、“Data For Slice Complementary”は、“Data For Slice”をビュー角再定義により得られる対向データを表し、“Data For Recon”は、画像再構成に必要な一揃いの投影データを表している。重み付けのこの方法はより簡単に実行することも可能であるが、計算上、高価となる場合がある。
【0070】
次に、この方法は、共通の扇形ビーム再構成の工程を含んでいる。また、扇形から平行への回帰(rebinning)後、平行ビーム再構成を行なうなどの他のアプローチも適用可能である。計算の数を減らすために、重み付け関数値がゼロでない範囲のみが、第1の工程でのヘリカル解析中に算出可能である。
【0071】
図8は、上述の実施形態に任意で適用可能であるチャネル倍増を図示したものである。対向データのサンプリング位置は、サンプリングピッチの半分だけ直接データのサンプリング位置からx−y面においてずらされる。対向データのX線経路が直接データ(実際に測定した投影データ(実データ))から変位されるため、軸平面の空間的解像度は増加可能である。チャネルの数を同数に維持し、軸平面のデータのサンプリングピッチ(チャネル方向)を減らすことによって、同様の結果(対象が十分に小さい場合、空間的解像度を向上する)が得られる。そのスライスに対するデータは、以下の(5)式で与えられる。
【0072】
【数5】
【0073】
チャネル倍増は、対向データの重み付けの有無にかかわらず上記の実施形態に適用されることが可能である。換言すれば、対向データを重み付けしない(例えば、(2)式または(3)式を使用)場合、(5)式が適用され、対向データが重み付けされる場合、(5)式は、(4)式が適用される前に実行される。
【0074】
第2の例では、画像は同じ天板位置であるが異なるタイミングで再構成される。また、本実施形態は、上記の実施形態に関して説明されたものに任意に適用可能であり、上記の実施形態はチャネル倍増を含む。本実施形態である「タイミング・シフト」(TS)技術は、図9(a)および図9(b)に示すように、時間的軸に沿って重み付け関数をずらす。ずらす時間量は、特定の適用の要求または必要とされる再構成画像の品位により決定される。例えば、再構成画像品位を向上させるためには、短い時間をずらせばよいが、計算の費用が増加する。再構成された画像の通常形状のみが必要な場合は、長い(または最大の)時間をずらすことも可能である。
【0075】
ここで、tc およびt0 (zs )は、それぞれ、その領域の中心と、X線焦点がスライス面(z=zs )全体を通過する時間とに対応する。標準のデータ範囲は、tc =t0 (zs )の場合に定義される。t0 (zs )から変位されたタイミングと、変位されたタイミングの最大とは、それぞれ、ts およびtsmとして定義される。ts またはPの大絶対値のいずれも大幅な外挿の場合に画質の劣化を招き、tsmは、図9(a)に示すようなスライス厚の半分内で外挿補間を使用することに制限される場合がある。以下に挙げる(6)式は、様々なPにおけるtsmと再構成アルゴリズムを与える。
【0076】
【数6】
【0077】
ただし、k=0の場合、半スライス厚内の最も近接したものの外挿補間
k=1の場合、内挿補間のみ
k=1+fwの場合、フィルタ幅がfwの場合のzフィルタリング
上記の通りに各スライスを再構成することは、自動心臓体積測定再構成を実現するための技術と関連して使用されることが可能である。ヘリカル(螺旋形)データが得られるとき、EKG信号は保存される。各スライスに対する中心タイミング(tc )は、(1)心収縮位相を避けるため、および/または(2)最少のテーパー角度によって、得られたデータ範囲を選ぶように調整される。これを、図10および図11に図示する。
【0078】
図10の方法で、工程40においてzはz0 に設定され、工程41において中心タイミングTc はT0 (z)として算出される。その後における、工程42において、tc がEKG信号の許容範囲にあるか、すなわち、拡張位相にあるか否かが決定される。そうである場合は、工程46へ進み、画像がtc の場合のzで再構成される。そうでない場合、工程43へ進み、tcmおよびtcPが算出される。工程44において、tcmまたはtcPがEKG信号中の許容範囲(拡張位相)にあるか否かが決定される。そうである場合、tcはtcmからtcPの範囲に設定され(工程45)、画像はtcの場合のzで再構成される。工程44において、そうでない場合には、tcmおよびtcPは再度算出され、tcmまたはtcPが許容範囲にあるか否かが再度調べられる。
【0079】
工程47において、処理が終了したか否か決定される。そうである場合、工程48において処理は終端し、そうでない場合、zはdzだけインクリメントされ(工程49)、この処理は工程41に戻って、zの新しい値でtc が再計算される。図11は、拡張位相において得られるデータを使用して再構成された画像を示すものである。図11において、参照番号50は最少のテーパー角度を有する再構成スライスを示し、51は内挿のみを使用して再構成されたスライスを示し、52は、内挿および/または外挿を使用して再構成されたスライスを示す。
【0080】
また、全てのスライスがtc がt0 (Zs )に等しい場合に再構成される方法を最初に実行し、操作者がどのスライスをとるかを示すと、画像は異なるタイミングで再構成される。これを、図12に図示する。参照番号60は拡張位相においてとられるデータを示し、参照番号61は内挿だけを使用して再構成されたスライスを示す。
【0081】
この方法の流れ図を、図13に示す。工程70において、全ての画像は、tc =t0 (Zs )で再構成される。tcmおよびtcPは、工程71のようにして計算され、工程72において、zにタイミング変位があるか否か決定される。変位がない場合、この方法は工程79で終わる。工程72において「はい」の場合、tcmおよびtcPは再計算される(工程73)、全ての画像はtc =tcmおよびtcPのzで再構成される(工程74)。zでの画像のうちのいずれかを選択した表示があるか否かについて、別の確認が行なわれる(工程95)。表示がある場合、工程76において、zでの画像は選択画像として選択され、工程72へ戻る。工程75において、表示がない場合、tcmおよびtcPは再度計算され(工程77)、画像は工程78において、tcmおよびtcPの新しい値で再構成される。工程78は、その後、工程75へ戻る。
【0082】
実例
本発明による方法の実例について説明する。方法の物理的な評価は、時間および空間の(z)解像度、画像ノイズ、および運動する模型の平面上画像の正確度を評価するために、コンピューター・シミュレーションを使用して実行された。実行され、比較される方法は、TSを使用するHHSおよびTSを使用しないHHS、マルチスライスCTに対するHFI、および単一スライスCTに対する180LIであった。形状およびX線管回転時間は、0.5秒で900の投影で全360°走査を実現する4-スライス・ヘリカルCTスキャナ(Aquilion;日本、東京、ToshibaMedicalSystemsCompany)のものと同一であった。
【0083】
時間的および空間的(z)解像度および画像ノイズは、HHSおよびHFIに対するヘリカルピッチおよびフィルタ幅を変化させて評価された。他の評価において、ヘリカル・マルチスライス・ピッチは、後述するように、それぞれ、マルチスライスCTに対して2.5で固定され、単一スライスCTに対して1.0に固定された。公称スライス厚は、空間的解像度に対して2.0mmで、運動する人体模型に対して1.0であった。
【0084】
最初に、時間的解像度は、図14で示す5つの円筒人体模型を使用することを決定することにより評価された。
【0085】
5つの円筒ディスクは、直径20mm、高さは無限で、1000HUのコントラストを有する。時間的感度プロファイル(TSP)が計測され、5円筒人体模型に対して9つの投影(0.005秒に相当する)、0.1mmの再構成増分、2.5〜8.0のヘリカル・マルチスライス・ピッチ、0.0および1.0のフィルタ幅というパラメータを使用した。部分感度プロファイル(TSP)の測定されたFWHMおよびFWTMは、軸状CT(0.5秒)のものに正規化された。
【0086】
TSPを用いて時間的解像度を評価することについて説明する。極めて短い時間(0.005秒)に対応するいくつかの連続投影以外全てゼロである入力投影データは、時間軸に沿ったインパルス信号とみなされることが可能である。いくつかの投影は、空間的解像度の影響を除去するために十分長くて厚い、z軸に平行な円筒人体模型並列を使用して得られる。軸画像のTSPの空間的依存関係を評価するために、異なるxy位置でいくつかの円筒を有する人体模型が使用される(上述の通り)。
【0087】
画像は、短い時間(好ましくは約0.03秒未満)に対応するz-増分により再構成される。各場所でのTSPは、上記投影データを使用した以下の処理によって獲得可能である。(1)短い時間に対応するzの増分を用いて画像を再構成する、(2)再構成された画像の全ての各円筒の位置での対象範囲(ROI)の平均ピクセル値を計測する、(3)最大値によって、それらを正規化する、(4)時間軸に対してそれらをプロットする。この例では、ROIは中心で選択されたが、他の位置も可能である。
【0088】
さて、TSPを使用してどのように時間的解像度の量を定めるかを説明する。15(a)および図15(b)は2つのTSPを示す。FWTM(full-width-at-tenth-maximum)は、図15(a)のような広いプロファイルに対するFWHM(full-width-at-half-maximum)よりも、時間的感度プロファイルの良い記述子であって、矩形のプロファイル(例えば図15(b)に対するFWTA(full-width-at-tenth-area)より優先されなければならない。したがって、TSPのFWTMは、時間的解像度の量を定めるための好適な記述子である。
【0089】
空間的解像度を評価するために、部分感度プロファイル(SSP)が使用されることが好ましい(SpiralVolumetricCTwithSingle-Breath-HoldTechnique,ContinuousTransportandContinuousScannerRotation,Kalenderetal.,Radiology176:181-183(1990)参照)。SSPは、1つの軸の(x−y)面中に5つのディスクを含むシミュレートされた人体模型を使用して計測された(図14参照)。このディスクは、直径20mmおよび厚さ0.2mmを有する。マルチスライス・ヘリカルピッチは、2.5から4.5まで変化させた。画像は、0.1mmの再構成増分とともに再構成された。SSPのFWHMおよびFWTMは計測され、軸CTのものに正規化された。TSを使用しないHHS0 と、マルチスライス・ヘリカル・ピッチ2.5でTS- 0.20 を使用するHHS0 との2つのSSPが比較され、TSの影響が評価された。
【0090】
各再構成アルゴリズムの運動する被検体のための撮像性能は、動く人体模型を使用して評価された。シミュレートされた人体模型は、円筒の2つの組および1000HUのコントラストを有するボールの3つの組を有し、各組は、1、2および3mmの異なる直径を有する3つのボールから成る(図16)。円筒は無限の長さを有し、それらが動かないときに、ボールの全ては面z=0の中心に配置された。
【0091】
運動パターンに関しては、A組およびC組は静止しており、他は毎分周波数60または90周期で移動した。BとDの組はx方向に移動し、E組はz方向に移動する。1周期の運動パターンは、3つの部分から成る。(1)元の位置(図16)から、t1 秒のコサイン曲線の半分に沿って−xまたは−z方向に10mm移動する、(2)t2 秒のコサイン曲線の半分に沿って元の位置に戻る、(3)元の位置で、t3 秒間静止する。秒単位のパラメータ(t1 ,t2 ,t3 )は、それぞれ、毎分60周期で(0.08、0.62、0.30)、90周期で(0.05、0.41、0.20)であった。
【0092】
サイズ、形状および運動パターンは、心臓の壁、心臓の近くの肺の気管、血管の脈動、冠状動脈および冠状動脈石灰化の運動の一部をシミュレートして選択された。選択された単純な運動パターンは、例えば縮退、伸長、撚りおよび並行運動など、複雑な3次元運動(またはひずみ)から成る特定の運動をシミュレートしないが、動く人体模型によって、通常の動く被検体を撮像する機能の評価ができる。
【0093】
ヘリカル・スキャンは二回使用され、それぞれ、X線焦点が面z=0を「静止位相」および「高速運動位相」へ横断する時間(図16、tx )を調整した。静止位相に対して、時間は運動−3(tx がそれぞれ60周期で0.85、90周期で0.56)の中央にあった。高速運動位相に対して、時間は運動-2(両方の周期においてtx =0.15秒)の間で、HHSに必要な全てのデータがオブジェクトの移動とともに得られた。したがって、HHSは静止位相において利点があるが、高速運動位相において欠点があった。ヘリカル・マルチスライス・ピッチは、マルチスライスCTに対して2.5で固定され、シングルスライスCTに対して1.0で固定された。
【0094】
z=0での画像は、TSがある場合のHHSO とTSがない場合のHHSO と、HFI0 、HFI1 および180LIによって再構成された。結果の量を定めるために、1つのボール(図16の80)を含むように設定されたROIの最大および最小値が計測された。最大値は、動く被検体のコントラストの精度を表し、最小値は、アーチファクトの強度を示す。
【0095】
物理的な評価および人体模型調査
(1)時間的解像度の評価
計測された時間的感度プロファイルのすべては、投影データ(図18(a)〜図18(d))のために使用される重みを反映させた。図18(a)〜図18(d)において、中央の人体模型は実線と対応し、右の人体模型は、小破線、下および左の人体模型は大破線、上の人体模型は大−小破線に対応する。なお、下と左の人体模型の曲線は重なる。HHSは、ほぼ矩形で(図18(a))、最も狭いTSPを与えた。HHSO およびHHS1 は、同一のプロファイルが見られた。HHSのTSPは基底で僅かに変化し、ピークで大幅に変化し、180LIおよびHFIに対しては逆が見られる。HFI0 のTSP(図18(b))は比較的広い範囲の複数三角形から成り、時間がゼロ近くでのレスポンスはほぼ0であり、もう一方の時間では最高に達している。HFI1 のTSP(図18(c))は長手方向の濾過がHFI0 のプロファイルを平滑化することを示した。
【0096】
図18(a)〜図18(d)の曲線において、矩形の形状は理想的である(図18(a)は、本発明による方法において、理想に最も近い)。曲線の広がりは劣化を表す。図18(c)は一部に広がりを示すが、レスポンスは0.0の周辺で劣化している。時間とともにレスポンスが高速に変化すると、複雑な三角形の形状の劣化が大きくなる。図18(b)は、著しい劣化を示し、−0.4、0.0および0.4でのレスポンスは、ゼロである。図18(d)のレスポンス曲線も、理想的な形状からいくぶん劣化している。
【0097】
HHSを有するTSPのFWHMおよびFWTMは、全てのヘリカルピッチ(後述の図19(a)および図19(b)、表1参照)を有する0.50および0.58で固定された。表Iは、時間的解像度データを提供する。FWHMおよびFWTMデータは、軸(ステップ・アンド・シュート)CTのものに対して正規化された。各々は、0.5秒の回転スキャナで、それぞれ、0.25秒および0.29秒に対応し、ヘリカル・マルチスライス・ピッチが4.5(図19(a)および図19(b))未満である場合のHFIより良い。TSPのFWTM(図19(b))によれば、HFIは時間的解像度に関して好ましいヘリカルピッチを示さない。HFIについては、ヘリカルピッチが高いまたは管回転が速いと、走査時間が短縮するのに加えて軸の画像の時間的解像度が向上する。ヘリカル・マルチスライス・ピッチが4.5未満であるときシングルスライス・ヘリカルCT(180LI)より悪いが、ヘリカル・マルチスライス・ピッチが7.0より等しいか大きい場合、HHSとほぼ同じとなる。
【0098】
【表1】
【0099】
(2)空間的解像度の評価
ヘリカルハーフスキャンにおいて、セクション感度プロファイルのFWHMおよびFWTMは、単純にヘリカルピッチとともに増加した(下記の表2参照)。この傾向は、単一スライス・ヘリカルCTと同様であったが、HFIのそれと全く異なった。
【0100】
TSは部分感度プロファイルのFWHMおよびFWTMをあまり変化させなかったが、形状に影響を及ぼし、空間的変化を増加させた。図20(a)および図20(b)は、それぞれ、TSがないHHS0 とTS- 0.20 を有するHHS0 に対するFWHMおよびFWTMについて図示したものである。線の形状は、図18(a)〜図18(d)に対して上述した人体模型と一致する。
【0101】
【表2】
【0102】
(3)ノイズの評価
画像ノイズは、ほぼ一定のままであった(表2)。この傾向も、シングルスライス・ヘリカルCTのそれと同様であるが、HFIのそれとは全く異なる。ノイズSDは、2.5のヘリカル・マルチスライス・ピッチで、TSの変化範囲が−0.25〜0.25秒である場合、1.08から1.24まで変化した。
【0103】
(4)動く人体模型の評価
図21(a−iの部分を有する)は、毎分60周期で静止位相の場合に、ボール(図16の80)の範囲内で、およびその周辺で、最大および最小のCTの数を示す。HHSは、アーチファクトがない画像(a)を示した。被検体が高速に動く(e)〜(g)の場合、TS(e)のないHHSは低品位の画像を生成した。HFIO (f)は、被検体が静止しているときに得られたデータを使用するので(図18(b)参照)、最良の画像を生成した。TS技術はHHS(i)の画像品位を改善し、CTの数(i)の精度を徹底的に向上させ、最良の画像を生成した。使用したヘリカル・ピッチは、マルチスライスCTに対して2.5、シングルスライスCTに対して1.0であった。公称スライス厚は1.09mmで、人体模型コントラストは1000HU、背景は0HU、およびウィンドウ(幅、高さ)は(1000,0)であった。
【0104】
人体模型被検体が毎分90周期で運動する場合、この傾向は依然同じであった(図22)。なお、TSは、HHSの高速運動位相中の劣化CT値を、ほぼ静止位相のレベルまで回復可能である(図8)。この人体模型検査によって、HFI1 は画像の品位のレベルを一定に保つが、ぼやけて、TS技術を用いたHHSよりも決して良好ではない。ここでも、TSを用いたHHSは最良の撮像機能を示し、画像品位は、静止位相とほぼ同じであった。90周期のデータに対する条件は、60周期のデータに対するものと同じであった。
【0105】
臨床評価
(1)心臓部撮像
図23では、心収縮端と心拡張端との間で、異なるスライス位置(列1〜3)と同じスライス位置(列4)とでとられた画像を示す。これらの画像における心拍数は、毎分75であった。HFI1 画像は、左の心室および心室中隔の輪郭がぼやけており(列2の第1および第3の画像参照)、心筋層の厚みの変化を撮像できなかった(列2参照)。HFI1 は心筋層の輪郭をわずかに鮮明にしているが、傾向は同じであった(列1参照)。HHSでは、心収縮端と心拡張端の両方で明確にそれらの特徴を撮像するのに成功した(列3、第1および第3の画像参照)。時間的解像度が向上したため、HHS(列3および4)は、心収縮端と心拡張端との間の全ての画像に対して、HFI(列1および2)よりも鮮明な心臓の左の心室および心室中隔を撮像した。心収縮位相と心拡張位相との間で、心臓の壁の運動は遅くなる。このため、大部分の必要なデータをこの期間内に入手可能である場合、HHSは心拡張端でだけでなく心収縮端ででも、鮮明な心臓境界線を有して画像を生成可能であった。
【0106】
0.5秒のスキャナ回転、マルチスライス・ヘリカル・CTおよびヘリカルハーフスキャン再構成が使用され、心収縮端および心拡張端(列3の第1および第3の画像)で鮮明な心筋層画像が可能となった。ウィンドウの幅および高さは、それぞれ、400および40で、標準のボディ・カーネル(FC10)が使用された。
【0107】
X線焦点が心収縮端と心拡張端との間の画像平面を横切るとき、HHSおよびTSの結合によって、心収縮端近くと心拡張端近くで、異なる時間に同じスライス位置で画像が得られるようになり、それによって、心臓の壁でみられたぼやけが低減する(列4および列3の第2の画像参照)。HHSによる画像ノイズは、HFIの場合より大きかった。
【0108】
(2)胸部撮像
図24(a)〜図24(d)は、心拡張端と心収縮端との間で得られた肺画像を示すものである。心臓の運動に起因して、HFI0 は強いアーチファクトを示した(図24(a)の81および82)。心臓の輪郭は、HFI1 の場合、二重線で図示され(図24(b)の83)、肺気管の縁はぼやけた(図24(b)の84)。TSのないHHSは運動アーチファクト(図24(c)の85)を示し、TS技術を使用した場合のHHSは、心臓の輪郭と肺気管の縁を鮮明にし、心臓の運動の影響を低減した(図24(d))。これらの画像において、心拍数は75/分で、標準の胸部カーネルが使われた(FC50)。
【0109】
(3)脈管の撮像
図25(a)〜図25(c)は、コントラスト素子を使用した胸の大動脈画像である。心拍数は、8.0/分で、標準のボディ・カーネルが使用された(FC10)。ウィンドウの幅および高さは、それぞれ600と0であった。大動脈の画像は、HHSがHFIより良好な画質を提供したことを示し、運動する被検体の映像を鮮明にすることが可能であることを示す。気管脈拍は解剖フラップを発生させる。それをHFI画像中では二重線で示すが、不十分な時間的解像度を示している(図25(a)および図25(b))。HHS画像は、解剖図を正確に描いている。
【0110】
本発明による方法において、HHSは心臓撮像における利点が以下のようにいくつかある。(1)非ゲート制御HHSは、患者スループットを増加させる、(2)ヘリカルピッチが大きいHHS、またはいずれゲート制御されるHHSのいずれかも患者線量を削減可能とする(3)HHSは、EKGゲート制御が使用されない場合、拡張位相の過去の心臓の撮像を可能とする、(4)HHSが連続したハーフスキャンデータをのみを使用するため、HHSはコントラスト素子を使用した検査においてより良好な画質を生成する可能性がある。例えば拡張位相など、比較的静止した位相を撮像する際、TSを使用したHHSは、EKGゲート制御再構成方法に匹敵する画質が得られると予測される。
【0111】
また、本発明による方法は、EKG信号のない拡張位相の3次元(3D)画像を得ることができる。心拡張位相中の各スライスに対して必要データ(180°+扇形角度)を獲得するための、ヘリカルピッチが小さい(例えば1.5)走査と、各スライスに対するTSのタイミングを手動で調節することとが結合することによって、3D画像中の段階型アーチファクトを低減するが、完全に除去するわけではない。このアーチファクトの原因は、異なるタイミングで得られたデータにより再構成された特定のスライス間の時間的間隙を含む。その間隙、したがってアーチファクトを増加させるのは、呼吸、患者の動き、コントラスト濃度の変化、またはその他のエラーのいずれかである。原理上の問題でなく実際的なものであるが、TSの自動タイミング調整のために使われるEKG信号は、タスクの負荷を減らすことができ、好ましい。
【0112】
なお、上述の説明では、ヘリカル補間、重み付け、コンボリューション、そしてファンビーム再構成の手順で処理を進めていたが、ファンビームに代えて、コーンビームでもよく、つまり、重み付け、コンボリューション、そしてファンビーム再構成の手順で処理を進める場合にも適用できる。
【0113】
本発明による方法の更なる実施形態として、患者線量を減らすことも可能である。一部のスライス・データ組は、ヘリカル内挿の前に加算されるか、平均化または重み付けされる(第1の例の第1の部分に記載した通り)。スライス位置は、全てのスライスの中心として定義される。また、平均算出過程での重み付けを考慮して定義されることが可能である。例えば、4スライス・データは、2スライス・データに平均化されることが可能である。この場合、4スライスからのデータの全てがヘリカル内挿で使われ、その結果、z方向の空間的解像度の一部の犠牲にしてより良好な画像ノイズが得られる。図26に示すように、患者線量は、患者に到達する前に、外側のスライスに対する不必要なX線ビームを視準することによって減少できる。大きいヘリカルピッチ、TSまたは長手方向の濾過は、外側のスライスからのデータを必要とする。その他の場合は外側のスライスからのデータが使用されず、不必要な患者線量となる。外側のスライスを除外し、二重スライスCTにするようにX線ビームを視準することによって、画質または走査時間に影響を及ぼさずに、患者線量を減らすことができる。この結果、一部の場合において、2-スライス・ヘリカルCTとなることが可能である。
【0114】
患者線量を減らす更なる方法は、マルチスライス・ヘリカル・スキャンのゲート制御トリガーを使用することである。上述のように、EKG信号は、ヘリカル・スキャン中に格納される。このEKG信号は、X線照射に対するトリガーとして使用される。所定量の走査後、X線照射は、自動的に停止される。この処理は、「擬似」マルチスライス・ヘリカル・スキャン(X線照射のない連続した管回転および連続する台運動)中に繰り返す。走査される範囲(心拡張位相)にのみ限定するため、患者線量が低減可能となる。
【0115】
この方法を、図27に図示する。工程90において、ヘリカル・スキャンが開始される。ヘリカル・スキャンが終了すると、工程91において、チェックが行なわれる。「はい(yes)」の場合、この処理は工程95で終端する。「いいえ(no)」の場合、EKG信号が閾値を超過するか否かが決定される(工程92)。「いいえ」の場合、処理は工程91に戻る。「はい」の場合、工程94へ進み、処理は所定期間(例えば100ミリ秒)待つ。工程94において、所定期間、X線照射およびデータ・保存処理が実行される。この所定期間は、例えば300ミリ秒でもよい。次に、工程94は工程92へ戻り、ヘリカル・スキャンが終了したか否かを確認する。
【0116】
本発明によるこの方法の別の実施形態において、3次元画像が、連続タイミングで再構成される。ヘリカル・スキャンは、小さいヘリカルピッチで実行される。この修正において使用される再構成技術は、タイミング変位技術に関連して上述したものの1つである。時間変位を、ヘリカル・スキャンで可能となった全体積を再構成する際の小さい増分として選択することによって、その体積の一部の範囲は、小さいタイミング・ピッチで再構成可能となる。これを図28に示す。同図において、100は内挿(細線)のみを使用して再構成されたスライスを表し、101は外挿(極太線)を使用して再構成されたスライスを表す。
【0117】
明らかに、本発明の多数の修正および変化は、上記の教示を考慮して可能である。例えば、最も近い近傍の、または非線形内挿のいずれでも、ヘリカル内挿の提案された方法の工程は、その空間的解像度または画質を改善するために適用されることが可能である。重み付け工程において、全ての重み付けが異なる目的のために適用可能である。したがって、添付の特許請求の範囲の範囲内において、本明細書に具体的に述べた以外にも、本発明は実践可能であることが理解されなければならない。
【0118】
本発明は、上述した実施形態に限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で種々変形して実施することが可能である。さらに、上記実施形態には種々の段階が含まれており、開示される複数の構成要件における適宜な組み合わせにより種々の発明が抽出され得る。例えば、実施形態に示される全構成要件から幾つかの構成要件が削除されてもよい。
【0119】
【発明の効果】
本発明によれば、コンピュータ断層撮影(CT)の時間的解像度を向上することができる。
【図面の簡単な説明】
【図1】従来の重み付け関数を示す図。
【図2】従来において、修正されたアンダースキャン技術のための重み付けを示す図。
【図3】本発明の実施形態によるシステム構成図。
【図4】図3のガントリの斜視図。
【図5】本実施形態において、マルチスライス検出器を示す図。
【図6】図3の処理ユニット及び処理回路の構成を示すブロック図。
【図7】本実施形態によるデータの修正HFI処理と重み付けとを示す図。
【図8】本実施形態によるチャネル倍増を示す図。
【図9】本実施形態において、時間変位を示す図。
【図10】本実施形態によるEKGゲート制御を用いた画像の再構成を図示する流れ図。
【図11】本実施形態による画像再構成を示す図。
【図12】本実施形態による選択的タイミングを用いた画像再構成を示す図。
【図13】図12の方法の流れ図。
【図14】本発明を評価する際に使用される人体模型を示す図。
【図15】本実施形態において、時間的な感度プロファイルを示す図。
【図16】運動する人体模型の図。
【図17】図16の人体模型の運動を示す図。
【図18】本実施形態において、時間的な感度プロファイルを示す図。
【図19】本実施形態において、FWHMおよびFWTMに対するヘリカル・ピッチをそれぞれ有する時間的解像度の変化を示す図。
【図20】本実施形態において、部分感度プロファイルを示す図。
【図21】本実施形態において、毎秒60周期でのシミュレーション模型画像を示す中間調画像。
【図22】本実施形態において、毎秒90周期でのシミュレーション模型画像を示す中間調画像。
【図23】本実施形態において、HFI、HHSおよびHHSさらにTSを使用して生成された心臓部画像を示す中間調画像。
【図24】本実施形態において、HFI、HHSおよびHHSさらにTSを使用して生成された胸部画像を示す中間調画像。
【図25】本実施形態において、HFIおよびHHSを使用して生成された脈管画像を示す中間調画像。
【図26】本実施形態において、患者線量を減らすための視準線を示す図。
【図27】本実施形態によるEKGゲート制御を示す流れ図。
【図28】本実施形態による体積再構成を示す図。
【符号の説明】
1…ガントリ、
2…回転リング、
3…X線源、
4…X線フィルタ、
5…2次元配列型X線検出器、
6…寝台、
7…高電圧発生器、
8…X線制御装置、
9…ガントリ/寝台制御装置、
10…システム制御装置、
11…データ収集ユニット、
12…信号処理ユニット、
14…ディスプレイ。
Claims (1)
- X線ビーム発生源と、前記X線ビーム発生源から被検体を介して到達するX線を検出するための複数の検出素子がチャンネル方向及びスライス方向に配列されてなる2次元検出器とを有するガントリと、
前記被検体を載置する寝台と、
前記寝台と前記ガントリとの少なくとも一方を移動する移動ユニットと、
前記2次元検出器に接続され、前記2次元検出器によって検出されたX線のデータを収集するデータ収集ユニットと、
前記データ収集ユニットに接続され、前記データ収集ユニットによって収集されたデータに対して前記スライス方向に関して内挿または外挿補間処理することにより所定のスライス位置に対応する投影データを生成し、前記投影データに対して時間方向に関して異なる複数の重み付け関数により重み付けをすることにより複数の投影データ組を生成する重み付け回路と、
前記複数のデータ組から、スライス位置が同じであるがデータ収集時刻が異なる複数の断層画像をそれぞれ再構成する断層画像再構成ユニットとを具備することを特徴とするコンピュータ断層撮影システム。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US09/450,121 US6466640B1 (en) | 1999-11-26 | 1999-11-26 | Computed tomography system and method |
US09/450121 | 1999-11-26 |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2011087269A Division JP5248648B2 (ja) | 1999-11-26 | 2011-04-11 | コンピュータ断層撮影システムおよび方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2001149365A JP2001149365A (ja) | 2001-06-05 |
JP5031941B2 true JP5031941B2 (ja) | 2012-09-26 |
Family
ID=23786849
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2000327188A Expired - Fee Related JP5031941B2 (ja) | 1999-11-26 | 2000-10-26 | コンピュータ断層撮影システム |
JP2011087269A Expired - Lifetime JP5248648B2 (ja) | 1999-11-26 | 2011-04-11 | コンピュータ断層撮影システムおよび方法 |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2011087269A Expired - Lifetime JP5248648B2 (ja) | 1999-11-26 | 2011-04-11 | コンピュータ断層撮影システムおよび方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US6466640B1 (ja) |
JP (2) | JP5031941B2 (ja) |
Families Citing this family (39)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7848790B2 (en) * | 2002-02-28 | 2010-12-07 | General Electric Company | System and method of imaging using a variable speed for thorax imaging |
US7813473B2 (en) * | 2002-07-23 | 2010-10-12 | General Electric Company | Method and apparatus for generating temporally interpolated projections |
US6904118B2 (en) * | 2002-07-23 | 2005-06-07 | General Electric Company | Method and apparatus for generating a density map using dual-energy CT |
US6754299B2 (en) * | 2002-08-02 | 2004-06-22 | Ge Medical Systems Global Technology Company, Llc | Methods and apparatus for weighting of computed tomography data |
US6865250B2 (en) * | 2002-12-23 | 2005-03-08 | Ge Medical Systems Global Technology Company Llc | High pitch cardiac helical scan with extended reconstruction windows |
DE602004021683D1 (de) * | 2003-04-10 | 2009-08-06 | Koninkl Philips Electronics Nv | Computertomographische methode für ein sich periodisch bewegendes objekt |
US6917663B2 (en) * | 2003-06-16 | 2005-07-12 | Kabushiki Kaisha Toshiba | Cone-beam reconstruction apparatus and computed tomography apparatus |
EP1684636B1 (en) * | 2003-11-05 | 2011-08-17 | General Electric Company | Method and apparatus for assisting cardiac resynchronization therapy |
DE10352380A1 (de) * | 2003-11-10 | 2005-06-16 | Siemens Ag | Verfahren zur Erzeugung von CT-Bildern von einem zyklisch bewegten Untersuchungsobjekt |
US7639774B2 (en) * | 2003-12-23 | 2009-12-29 | General Electric Company | Method and apparatus for employing multiple axial-sources |
DE10360981A1 (de) * | 2003-12-23 | 2005-07-21 | Siemens Ag | Verfahren zur Erzeugung von CT-Bildern eines Untersuchungsobjektes mit einem sich periodisch bewegenden Teilbereich |
JP4119835B2 (ja) * | 2003-12-26 | 2008-07-16 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 被曝線量計算方法およびx線撮影装置 |
US7333587B2 (en) | 2004-02-27 | 2008-02-19 | General Electric Company | Method and system for imaging using multiple offset X-ray emission points |
JP4795658B2 (ja) * | 2004-03-23 | 2011-10-19 | 富士フイルム株式会社 | 差分画像取得方法、差分画像取得装置、及び、そのプログラム |
CN1977288B (zh) * | 2004-06-28 | 2011-10-05 | 皇家飞利浦电子股份有限公司 | 用于灌注研究的检查装置 |
DE102004031984B4 (de) * | 2004-07-01 | 2007-01-11 | Siemens Ag | Verfahren zur Erfassung und Visualisierung dynamischer Prozesse in einem Objektvolumen |
US7734079B2 (en) * | 2004-09-28 | 2010-06-08 | General Electric Company | Methods and apparatus for image reconstruction |
JP2006129975A (ja) * | 2004-11-04 | 2006-05-25 | Ge Medical Systems Global Technology Co Llc | X線ct装置およびx線ct撮影方法 |
JP4828839B2 (ja) * | 2005-03-07 | 2011-11-30 | 株式会社東芝 | X線コンピュータ断層撮影装置、画像処理装置及び画像処理方法 |
US7570733B2 (en) * | 2005-06-10 | 2009-08-04 | General Electric Company | Step-and-shoot cardiac CT imaging |
JP4495109B2 (ja) * | 2006-04-06 | 2010-06-30 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X線ct装置 |
US7616731B2 (en) * | 2006-08-30 | 2009-11-10 | General Electric Company | Acquisition and reconstruction of projection data using a stationary CT geometry |
US20080056432A1 (en) * | 2006-08-30 | 2008-03-06 | General Electric Company | Reconstruction of CT projection data |
US7835486B2 (en) * | 2006-08-30 | 2010-11-16 | General Electric Company | Acquisition and reconstruction of projection data using a stationary CT geometry |
US7706499B2 (en) * | 2006-08-30 | 2010-04-27 | General Electric Company | Acquisition and reconstruction of projection data using a stationary CT geometry |
US8233682B2 (en) * | 2007-06-05 | 2012-07-31 | General Electric Company | Methods and systems for improving spatial and temporal resolution of computed images of moving objects |
US8200318B2 (en) * | 2007-09-21 | 2012-06-12 | University of Pittsburgh—of the Commonwealth System of Higher Education | Electrocardiogram reconstruction from implanted device electrograms |
US9980661B2 (en) * | 2007-09-21 | 2018-05-29 | University of Pittsburgh—of the Commonwealth System of Higher Education | Electrocardiogram reconstruction from implanted device electrograms |
JP5085305B2 (ja) * | 2007-12-21 | 2012-11-28 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X線ct装置 |
US8218716B2 (en) * | 2009-09-14 | 2012-07-10 | Mitsubishi Heavy Industries, Ltd. | Radiation tomography method |
CN103202704B (zh) * | 2012-01-11 | 2015-03-04 | 中国科学院深圳先进技术研究院 | 半扫描位置的确定方法 |
US9161724B2 (en) * | 2012-07-20 | 2015-10-20 | Koninklijke Philips N.V. | Multi-cardiac sound gated imaging and post-processing of imaging data based on cardiac sound |
JP6225186B2 (ja) * | 2012-08-08 | 2017-11-01 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | コンピュータ断層撮影(ct)用の慢性閉塞性肺疾患(copd)ファントム及びそれを使用する方法 |
DE102013202491B4 (de) * | 2012-11-08 | 2021-11-04 | Siemens Healthcare Gmbh | Computertomograph und Verfahren zur Aufnahme eines Bildes des Herzens |
US9697624B2 (en) * | 2014-08-28 | 2017-07-04 | Shimadzu Corporation | Image processing apparatus, radiation tomography apparatus, and method of performing image processing |
JP6273241B2 (ja) * | 2015-09-24 | 2018-01-31 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 放射線断層撮影方法及び装置並びにプログラム |
KR20170105876A (ko) | 2016-03-10 | 2017-09-20 | 삼성전자주식회사 | 단층 촬영 장치 및 그에 따른 단층 영상 재구성 방법 |
CN107041760B (zh) * | 2017-04-27 | 2021-02-05 | 东软医疗系统股份有限公司 | 扫描方法、装置及图像重建方法和装置 |
WO2019041101A1 (en) * | 2017-08-28 | 2019-03-07 | Shenzhen United Imaging Healthcare Co., Ltd. | SYSTEMS AND METHODS FOR DETERMINING ROTATION ANGLES |
Family Cites Families (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH053867A (ja) * | 1991-06-28 | 1993-01-14 | Toshiba Corp | 三次元画像診断装置 |
IL98945A0 (en) * | 1991-07-24 | 1992-07-15 | Elscint Ltd | Multiple slice ct scanner |
US5966422A (en) * | 1992-07-20 | 1999-10-12 | Picker Medical Systems, Ltd. | Multiple source CT scanner |
US5271055A (en) * | 1992-08-19 | 1993-12-14 | General Electric Company | Methods for reducing motion induced artifacts in a projection imaging system |
JP3637074B2 (ja) * | 1992-12-15 | 2005-04-06 | 株式会社東芝 | ヘリカルスキャン方式のコンピュータ断層撮影装置 |
JP2914891B2 (ja) | 1995-07-05 | 1999-07-05 | 株式会社東芝 | X線コンピュータ断層撮影装置 |
JP4316017B2 (ja) * | 1995-09-11 | 2009-08-19 | ジーイー横河メディカルシステム株式会社 | X線ct装置 |
US5559847A (en) * | 1995-12-06 | 1996-09-24 | General Electric Company | Systems, methods and apparatus for reconstructing images in a CT system implementing a helical scan |
JP3455041B2 (ja) | 1995-12-25 | 2003-10-06 | 株式会社東芝 | X線ct装置 |
US5974108A (en) * | 1995-12-25 | 1999-10-26 | Kabushiki Kaisha Toshiba | X-ray CT scanning apparatus |
US5708691A (en) * | 1996-07-05 | 1998-01-13 | Kabushiki Kaisha Toshiba | X-ray computed tomographic imaging device and x-ray computed tomographic method |
JP3763611B2 (ja) | 1996-07-12 | 2006-04-05 | 株式会社東芝 | X線ctスキャナ |
US5848117A (en) * | 1996-11-27 | 1998-12-08 | Analogic Corporation | Apparatus and method for computed tomography scanning using halfscan reconstruction with asymmetric detector system |
US5828718A (en) * | 1996-11-27 | 1998-10-27 | Analogic Corporation | Method and apparatus for helical computed tomography scanning with asymmetric detector system |
JP4155529B2 (ja) | 1997-02-17 | 2008-09-24 | 株式会社東芝 | X線ct装置 |
JP3124254B2 (ja) * | 1997-07-24 | 2001-01-15 | ジーイー横河メディカルシステム株式会社 | 放射線断層撮影装置 |
US6185275B1 (en) * | 1998-08-25 | 2001-02-06 | General Electric Company | Systems and methods for correcting focal spot thermal drift |
US6154516A (en) * | 1998-09-04 | 2000-11-28 | Picker International, Inc. | Cardiac CT system |
US6275560B1 (en) * | 1998-12-22 | 2001-08-14 | General Electric Company | Cardiac gated computed tomography system |
-
1999
- 1999-11-26 US US09/450,121 patent/US6466640B1/en not_active Expired - Lifetime
-
2000
- 2000-10-26 JP JP2000327188A patent/JP5031941B2/ja not_active Expired - Fee Related
-
2011
- 2011-04-11 JP JP2011087269A patent/JP5248648B2/ja not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
JP2001149365A (ja) | 2001-06-05 |
JP5248648B2 (ja) | 2013-07-31 |
JP2011136219A (ja) | 2011-07-14 |
US6466640B1 (en) | 2002-10-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5031941B2 (ja) | コンピュータ断層撮影システム | |
Taguchi et al. | High temporal resolution for multislice helical computed tomography | |
JP4666915B2 (ja) | 灌流測定のための動き補正 | |
JP4340533B2 (ja) | コンピュータ断層撮影 | |
US6243437B1 (en) | Coronary calcification detection using retrospective cardiac gating of imaging system | |
US7570733B2 (en) | Step-and-shoot cardiac CT imaging | |
US6504894B2 (en) | Phase-driven multisector reconstruction for multislice helical CT imaging | |
Flohr et al. | Performance evaluation of a multi-slice CT system with 16-slice detector and increased gantry rotation speed for isotropic submillimeter imaging of the heart | |
JP5199081B2 (ja) | 心臓ct撮影のバンドアーチファクトの抑制 | |
JP2004160222A (ja) | コンピュータ断層撮影方法および装置 | |
US8483804B2 (en) | Method and apparatus for reconstructing images of moving structures based on motion cycle temporal data | |
JP3944173B2 (ja) | 放射線画像処理装置及び処理方法 | |
JP2002345808A (ja) | Ctスカウト画像処理のための方法及び装置 | |
JP2007000408A (ja) | X線ct装置 | |
CN106920265A (zh) | 计算机断层扫描图像重建方法及装置 | |
CN103390284A (zh) | 在扩展的测量场中的ct图像重建 | |
JP2009538205A (ja) | 動的コンピュータ断層撮像 | |
JPH11253435A (ja) | コンピュ―タトモグラフ | |
Taguchi et al. | Direct cone‐beam cardiac reconstruction algorithm with cardiac banding artifact correction | |
Bruder et al. | High-temporal-resolution volume heart imaging with multirow computed tomography | |
JPH06205771A (ja) | 高速掃引を用いた四次元螺旋状体積撮像方法及び装置 | |
CN101352352B (zh) | X射线ct装置 | |
JP4571622B2 (ja) | 周期的運動をする対象のコンピュータ断層撮影法 | |
US6647084B1 (en) | Method and apparatus for filtering projection data of a helical scan | |
US7023958B2 (en) | Radiation image-acquiring apparatus, and radiation image-acquiring method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20071025 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20100601 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20100713 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20100906 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20101102 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20110104 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20110208 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20110411 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20120207 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20120507 |
|
A911 | Transfer to examiner for re-examination before appeal (zenchi) |
Free format text: JAPANESE INTERMEDIATE CODE: A911 Effective date: 20120514 |
|
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: 20120605 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20120628 |
|
R151 | Written notification of patent or utility model registration |
Ref document number: 5031941 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20150706 Year of fee payment: 3 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 Free format text: JAPANESE INTERMEDIATE CODE: R313114 Free format text: JAPANESE INTERMEDIATE CODE: R313117 |
|
R371 | Transfer withdrawn |
Free format text: JAPANESE INTERMEDIATE CODE: R371 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 Free format text: JAPANESE INTERMEDIATE CODE: R313114 Free format text: JAPANESE INTERMEDIATE CODE: R313117 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
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 |
|
LAPS | Cancellation because of no payment of annual fees |