JP2927209B2 - ディジタルアンギオグラフィ装置 - Google Patents
ディジタルアンギオグラフィ装置Info
- Publication number
- JP2927209B2 JP2927209B2 JP7148191A JP14819195A JP2927209B2 JP 2927209 B2 JP2927209 B2 JP 2927209B2 JP 7148191 A JP7148191 A JP 7148191A JP 14819195 A JP14819195 A JP 14819195A JP 2927209 B2 JP2927209 B2 JP 2927209B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- frequency
- imaging
- subtraction
- ray
- 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
- 238000002583 angiography Methods 0.000 title claims description 14
- 238000003384 imaging method Methods 0.000 claims description 125
- 238000006243 chemical reaction Methods 0.000 claims description 60
- 230000005540 biological transmission Effects 0.000 claims description 24
- 239000002872 contrast media Substances 0.000 claims description 24
- 238000000034 method Methods 0.000 claims description 22
- 230000010354 integration Effects 0.000 claims description 8
- 230000001678 irradiating effect Effects 0.000 claims description 7
- 238000010586 diagram Methods 0.000 description 14
- 210000004204 blood vessel Anatomy 0.000 description 9
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000003111 delayed effect Effects 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000005855 radiation 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/48—Diagnostic techniques
- A61B6/481—Diagnostic techniques involving the use of contrast agents
-
- 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
-
- 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/46—Arrangements for interfacing with the operator or the patient
- A61B6/461—Displaying means of special interest
- A61B6/464—Displaying means of special interest involving a plurality of displays
-
- 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/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/504—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of blood vessels, e.g. by angiography
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N5/00—Details of television systems
- H04N5/30—Transforming light or analogous information into electric information
- H04N5/32—Transforming X-rays
- H04N5/321—Transforming X-rays with video transmission of fluoroscopic images
- H04N5/325—Image enhancement, e.g. by subtraction techniques using polyenergetic X-rays
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)
- High Energy & Nuclear Physics (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Public Health (AREA)
- Heart & Thoracic Surgery (AREA)
- General Health & Medical Sciences (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Vascular Medicine (AREA)
- Dentistry (AREA)
- Human Computer Interaction (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Description
【０００１】
【産業上の利用分野】この発明は、被検体の１または複
数の撮像部位のサブトラクション像を好適に得ることが
できるディジタルアンギオグラフィ装置に関する。
数の撮像部位のサブトラクション像を好適に得ることが
できるディジタルアンギオグラフィ装置に関する。
【０００２】
【従来の技術】従来のこの種の装置で、ある１箇所の撮
像部位に対するサブトラクション像を得る場合、まず、
造影剤投与前にその部位のマスク像を撮像し、次に造影
剤を投与し、そして、造影剤投与後にその部位のライブ
像を撮像し、ライブ像からマスク像をサブトラクション
している。
像部位に対するサブトラクション像を得る場合、まず、
造影剤投与前にその部位のマスク像を撮像し、次に造影
剤を投与し、そして、造影剤投与後にその部位のライブ
像を撮像し、ライブ像からマスク像をサブトラクション
している。
【０００３】また、Ｘ線管や撮像系（イメージインテン
シファイアやテレビカメラなどで構成される）を備えた
Ｘ線透視装置と、被検体との相対的な位置関係を、被検
体の体軸方向などに変位させながら、複数の撮像部位に
対するサブトラクション像を得る場合でも、造影剤の投
与を挟んで、各撮像部位のマスク像の撮像と各撮像部位
のライブ像の撮像を２回に分けて行っている。
シファイアやテレビカメラなどで構成される）を備えた
Ｘ線透視装置と、被検体との相対的な位置関係を、被検
体の体軸方向などに変位させながら、複数の撮像部位に
対するサブトラクション像を得る場合でも、造影剤の投
与を挟んで、各撮像部位のマスク像の撮像と各撮像部位
のライブ像の撮像を２回に分けて行っている。
【０００４】
【発明が解決しようとする課題】しかしながら、このよ
うな構成を有する従来例の場合には、次のような問題が
ある。すなわち、上述した従来例では、マスク像撮像時
とライブ像撮像時にそれぞれ被検体にＸ線を照射してお
り、撮像部位ごとに２回のＸ線照射を必要とし、被検体
へのＸ線曝射線量の軽減が図れないという問題がある。
うな構成を有する従来例の場合には、次のような問題が
ある。すなわち、上述した従来例では、マスク像撮像時
とライブ像撮像時にそれぞれ被検体にＸ線を照射してお
り、撮像部位ごとに２回のＸ線照射を必要とし、被検体
へのＸ線曝射線量の軽減が図れないという問題がある。
【０００５】また、従来例の構成では、マスク像の撮像
と、ライブ像の撮像とを、造影剤投与を挟んで２回に分
けて行うので、マスク像撮像とライブ像撮像との間の時
間間隔が長く、また、その間に造影剤の投与も行うこと
から、被検体が体動し易く、これに起因してマスク像と
ライブ像の画像にずれが生じ、得られたサブトラクショ
ン像にアーティファクトが発生していた。また、Ｘ線透
視装置と被検体との相対的な位置関係を変位させながら
複数の撮像部位のサブトラクション像を得る場合には、
さらに、Ｘ線透視装置と被検体との相対的な位置関係を
変位させる機構の機械的な誤差により、同一の撮像部位
に対してマスク像とライブ像とを撮像するように制御し
ても、それらを一致させるのが難しく、マスク像とライ
ブ像の画像は一層ずれが目立っていた。
と、ライブ像の撮像とを、造影剤投与を挟んで２回に分
けて行うので、マスク像撮像とライブ像撮像との間の時
間間隔が長く、また、その間に造影剤の投与も行うこと
から、被検体が体動し易く、これに起因してマスク像と
ライブ像の画像にずれが生じ、得られたサブトラクショ
ン像にアーティファクトが発生していた。また、Ｘ線透
視装置と被検体との相対的な位置関係を変位させながら
複数の撮像部位のサブトラクション像を得る場合には、
さらに、Ｘ線透視装置と被検体との相対的な位置関係を
変位させる機構の機械的な誤差により、同一の撮像部位
に対してマスク像とライブ像とを撮像するように制御し
ても、それらを一致させるのが難しく、マスク像とライ
ブ像の画像は一層ずれが目立っていた。
【０００６】この発明は、このような事情に鑑みてなさ
れたものであって、サブトラクション像を得るための被
検体へのＸ線曝射線量を軽減するとともに、マスク像と
ライブ像の画像のずれを完全に無くし得るディジタルア
ンギオグラフィ装置を提案することを目的とする。
れたものであって、サブトラクション像を得るための被
検体へのＸ線曝射線量を軽減するとともに、マスク像と
ライブ像の画像のずれを完全に無くし得るディジタルア
ンギオグラフィ装置を提案することを目的とする。
【０００７】
【課題を解決するための手段】この発明は、このような
目的を達成するために、次のような構成をとる。すなわ
ち、請求項１に記載の発明は、被検体の所定の撮像部位
のサブトラクション像を得るためのディジタルアンギオ
グラフィ装置であって、（１−ａ）前記撮像部位にＸ線
を照射し、その撮像部位のＸ線透過像を撮像するＸ線透
視手段と、（１−ｂ）前記Ｘ線透過像をディジタルデー
タに変換するデータ変換手段と、（１−ｃ）造影剤が投
与された被検体の所定の撮像部位を撮像し、ディジタル
データに変換されたＸ線透過像（以下、基本画像とい
う）を構成する周波数成分を空間／周波数変換処理で取
り出し、所定のしきい値周波数以上の周波数成分を除去
し、これに周波数／空間変換処理を施して、前記基本画
像から前記しきい値周波数以上の高周波数成分を除去し
た画像（以下、高周波数除去画像という）を得る周波数
特性変換手段と、（１−ｄ）前記基本画像と前記高周波
数除去画像とのサブトラクションを行いその撮像部位の
サブトラクション像を求める演算手段と、（１−ｅ）前
記演算手段への基本画像の供給を遅らせて、前記基本画
像と前記高周波数除去画像とを前記演算手段に同期して
供給させる遅延手段とを備えたものである。
目的を達成するために、次のような構成をとる。すなわ
ち、請求項１に記載の発明は、被検体の所定の撮像部位
のサブトラクション像を得るためのディジタルアンギオ
グラフィ装置であって、（１−ａ）前記撮像部位にＸ線
を照射し、その撮像部位のＸ線透過像を撮像するＸ線透
視手段と、（１−ｂ）前記Ｘ線透過像をディジタルデー
タに変換するデータ変換手段と、（１−ｃ）造影剤が投
与された被検体の所定の撮像部位を撮像し、ディジタル
データに変換されたＸ線透過像（以下、基本画像とい
う）を構成する周波数成分を空間／周波数変換処理で取
り出し、所定のしきい値周波数以上の周波数成分を除去
し、これに周波数／空間変換処理を施して、前記基本画
像から前記しきい値周波数以上の高周波数成分を除去し
た画像（以下、高周波数除去画像という）を得る周波数
特性変換手段と、（１−ｄ）前記基本画像と前記高周波
数除去画像とのサブトラクションを行いその撮像部位の
サブトラクション像を求める演算手段と、（１−ｅ）前
記演算手段への基本画像の供給を遅らせて、前記基本画
像と前記高周波数除去画像とを前記演算手段に同期して
供給させる遅延手段とを備えたものである。
【０００８】請求項２に記載の発明は、被検体の所定の
撮像部位のサブトラクション像を得るためのディジタル
アンギオグラフィ装置であって、（２−ａ）前記撮像部
位にＸ線を照射し、その撮像部位のＸ線透過像を撮像す
るＸ線透視手段と、（２−ｂ）前記Ｘ線透過像をディジ
タルデータに変換するデータ変換手段と、（２−ｃ）造
影剤が投与された被検体の所定の撮像部位を撮像し、デ
ィジタルデータに変換されたＸ線透過像（以下、基本画
像という）を構成する周波数成分を空間／周波数変換処
理で取り出し、所定のしきい値周波数以上の周波数成分
を除去し、これに周波数／空間変換処理を施して、前記
基本画像から前記しきい値周波数以上の高周波数成分を
除去した画像（以下、高周波数除去画像という）を得る
周波数特性変換手段と、（２−ｄ）前記基本画像の高周
波数成分を強調した画像（以下、高周波数強調画像とい
う）を得る第２の周波数特性変換手段と、（２−ｅ）前
記高周波数強調画像と前記高周波数除去画像とのサブト
ラクションを行いその撮像部位のサブトラクション像を
求める演算手段とを備えたものである。
撮像部位のサブトラクション像を得るためのディジタル
アンギオグラフィ装置であって、（２−ａ）前記撮像部
位にＸ線を照射し、その撮像部位のＸ線透過像を撮像す
るＸ線透視手段と、（２−ｂ）前記Ｘ線透過像をディジ
タルデータに変換するデータ変換手段と、（２−ｃ）造
影剤が投与された被検体の所定の撮像部位を撮像し、デ
ィジタルデータに変換されたＸ線透過像（以下、基本画
像という）を構成する周波数成分を空間／周波数変換処
理で取り出し、所定のしきい値周波数以上の周波数成分
を除去し、これに周波数／空間変換処理を施して、前記
基本画像から前記しきい値周波数以上の高周波数成分を
除去した画像（以下、高周波数除去画像という）を得る
周波数特性変換手段と、（２−ｄ）前記基本画像の高周
波数成分を強調した画像（以下、高周波数強調画像とい
う）を得る第２の周波数特性変換手段と、（２−ｅ）前
記高周波数強調画像と前記高周波数除去画像とのサブト
ラクションを行いその撮像部位のサブトラクション像を
求める演算手段とを備えたものである。
【０００９】請求項３に記載の発明は、被検体の所定の
撮像部位のサブトラクション像を得るためのディジタル
アンギオグラフィ装置であって、（３−ａ）前記撮像部
位にＸ線を照射し、その撮像部位のＸ線透過像を撮像す
るＸ線透視手段と、（３−ｂ）前記Ｘ線透過像をディジ
タルデータに変換するデータ変換手段と、（３−ｃ）造
影剤が投与された被検体の所定の撮像部位を撮像し、デ
ィジタルデータに変換されたＸ線透過像（以下、基本画
像という）を順次取り込み、複数回分の撮像で得られた
基本画像の平均画像を求める積分処理手段と、（３−
ｄ）前記平均画像を構成する周波数成分を空間／周波数
変換処理で取り出し、所定のしきい値周波数以上の周波
数成分を除去し、これに周波数／空間変換処理を施し
て、前記平均画像から前記しきい値周波数以上の高周波
数成分を除去した画像（以下、高周波数除去画像とい
う）を得る周波数特性変換手段と、（３−ｄ）前記平均
画像と前記高周波数除去画像とのサブトラクションを行
いその撮像部位のサブトラクション像を求める演算手段
とを備えたものである。
撮像部位のサブトラクション像を得るためのディジタル
アンギオグラフィ装置であって、（３−ａ）前記撮像部
位にＸ線を照射し、その撮像部位のＸ線透過像を撮像す
るＸ線透視手段と、（３−ｂ）前記Ｘ線透過像をディジ
タルデータに変換するデータ変換手段と、（３−ｃ）造
影剤が投与された被検体の所定の撮像部位を撮像し、デ
ィジタルデータに変換されたＸ線透過像（以下、基本画
像という）を順次取り込み、複数回分の撮像で得られた
基本画像の平均画像を求める積分処理手段と、（３−
ｄ）前記平均画像を構成する周波数成分を空間／周波数
変換処理で取り出し、所定のしきい値周波数以上の周波
数成分を除去し、これに周波数／空間変換処理を施し
て、前記平均画像から前記しきい値周波数以上の高周波
数成分を除去した画像（以下、高周波数除去画像とい
う）を得る周波数特性変換手段と、（３−ｄ）前記平均
画像と前記高周波数除去画像とのサブトラクションを行
いその撮像部位のサブトラクション像を求める演算手段
とを備えたものである。
【００１０】
【作用】請求項１に記載の発明の作用は次のとおりであ
る。造影剤が投与された被検体の所定の撮像部位に対し
て、Ｘ線透視手段でＸ線透過像を撮像し、データ変換手
段でディジタルデータに変換して基本画像を得る。この
基本画像は、造影剤が投与された血管像などの高周波数
成分が残っている画像であり、これをライブ像として用
いる。
る。造影剤が投与された被検体の所定の撮像部位に対し
て、Ｘ線透視手段でＸ線透過像を撮像し、データ変換手
段でディジタルデータに変換して基本画像を得る。この
基本画像は、造影剤が投与された血管像などの高周波数
成分が残っている画像であり、これをライブ像として用
いる。
【００１１】一方、周波数特性変換手段は、上記基本画
像を構成する周波数成分を空間／周波数変換処理で取り
出し、所定のしきい値周波数以上の周波数成分を除去
し、これに周波数／空間変換処理を施して、上記基本画
像から上記しきい値周波数以上の高周波数成分を除去し
た高周波数除去画像を得る。この高周波数除去画像は、
造影剤が投与された血管像などの高周波数成分が除去さ
れた画像であり、これをマスク像として用いる。なお、
上記しきい値周波数として、サブトラクション画像に残
す血管像などを好適に除去し得る周波数を設定しておく
で、適正なライブ像を得ることができる。
像を構成する周波数成分を空間／周波数変換処理で取り
出し、所定のしきい値周波数以上の周波数成分を除去
し、これに周波数／空間変換処理を施して、上記基本画
像から上記しきい値周波数以上の高周波数成分を除去し
た高周波数除去画像を得る。この高周波数除去画像は、
造影剤が投与された血管像などの高周波数成分が除去さ
れた画像であり、これをマスク像として用いる。なお、
上記しきい値周波数として、サブトラクション画像に残
す血管像などを好適に除去し得る周波数を設定しておく
で、適正なライブ像を得ることができる。
【００１２】そして、演算手段は、ライブ像（基本画
像）からマスク像（高周波数除去画像）をサブトラクシ
ョンしてサブトラクション像を求める。
像）からマスク像（高周波数除去画像）をサブトラクシ
ョンしてサブトラクション像を求める。
【００１３】ここで、高周波数除去画像を得るためには
所定の処理時間が必要であるが、この高周波数除去画像
は基本画像が得られてから求めるので、基本画像が得ら
れてすぐに基本画像を演算手段に供給するとともに、周
波数特性変換手段で高周波数除去画像が得られてすぐに
高周波数除去画像を演算手段に供給すると、演算手段へ
は基本画像と高周波数除去画像とが周波数特性変換手段
での処理時間分ずれて供給されることになる。そこで、
この請求項１に記載の発明では、周波数特性変換手段で
の処理時間による時間的な遅れを補償するために遅延手
段を設け、この遅延手段によって、演算手段への基本画
像の供給を遅らせて、基本画像と高周波数除去画像とを
演算手段に同期して供給させるようにしている。
所定の処理時間が必要であるが、この高周波数除去画像
は基本画像が得られてから求めるので、基本画像が得ら
れてすぐに基本画像を演算手段に供給するとともに、周
波数特性変換手段で高周波数除去画像が得られてすぐに
高周波数除去画像を演算手段に供給すると、演算手段へ
は基本画像と高周波数除去画像とが周波数特性変換手段
での処理時間分ずれて供給されることになる。そこで、
この請求項１に記載の発明では、周波数特性変換手段で
の処理時間による時間的な遅れを補償するために遅延手
段を設け、この遅延手段によって、演算手段への基本画
像の供給を遅らせて、基本画像と高周波数除去画像とを
演算手段に同期して供給させるようにしている。
【００１４】請求項２に記載の発明によれば、請求項１
に記載の発明と同様の作用で基本画像と高周波数除去画
像が得られる。そして、第２の周波数特性変換手段は、
得られた基本画像の高周波数成分を強調した高周波数強
調画像を得る。演算手段は、この高周波数強調画像と上
記高周波数除去画像とのサブトラクションを行いその撮
像部位のサブトラクション像を求める。
に記載の発明と同様の作用で基本画像と高周波数除去画
像が得られる。そして、第２の周波数特性変換手段は、
得られた基本画像の高周波数成分を強調した高周波数強
調画像を得る。演算手段は、この高周波数強調画像と上
記高周波数除去画像とのサブトラクションを行いその撮
像部位のサブトラクション像を求める。
【００１５】請求項３に記載の発明によれば、積分処理
手段が、複数回分の撮像で得られた基本画像の平均画像
を求め、この平均画像をライブ像として用いるととも
に、このライブ像（平均画像）から請求項１に記載の発
明と同様の処理で高周波数除去画像を得て、演算手段
が、ライブ像（平均画像）からマスク像（平均画像から
得られた高周波数除去画像）をサブトラクションしてサ
ブトラクション像を求める。
手段が、複数回分の撮像で得られた基本画像の平均画像
を求め、この平均画像をライブ像として用いるととも
に、このライブ像（平均画像）から請求項１に記載の発
明と同様の処理で高周波数除去画像を得て、演算手段
が、ライブ像（平均画像）からマスク像（平均画像から
得られた高周波数除去画像）をサブトラクションしてサ
ブトラクション像を求める。
【００１６】
【実施例】以下、図面を参照してこの発明の一実施例を
説明する。 ＜第１実施例＞ 図１は、この発明の第１実施例に係るディジタルアンギ
オグラフィ装置の全体構成を示す正面図であり、図２
は、Ｘ線透視装置を側面から見た図、図３は、画像処理
部の構成を示すブロック図である。
説明する。 ＜第１実施例＞ 図１は、この発明の第１実施例に係るディジタルアンギ
オグラフィ装置の全体構成を示す正面図であり、図２
は、Ｘ線透視装置を側面から見た図、図３は、画像処理
部の構成を示すブロック図である。
【００１７】この第１実施例は、ベッド１、Ｘ線透視手
段としてのＸ線透視装置２、画像処理部３、モニタ４、
制御部５、操作盤６などを備えて構成されている。
段としてのＸ線透視装置２、画像処理部３、モニタ４、
制御部５、操作盤６などを備えて構成されている。
【００１８】ベッド１は、床面に設置されたベッド基台
１１と天板１２を備えている。被検体Ｍは天板１２上に
載置される。この天板１２はモータ１３の駆動で水平移
動可能であり、天板１２上の被検体ＭとＸ線透視装置２
との相対的な位置関係を被検体Ｍの体軸方向に変位する
ことができる。モータ１３の駆動制御は、制御部５によ
り行われる。
１１と天板１２を備えている。被検体Ｍは天板１２上に
載置される。この天板１２はモータ１３の駆動で水平移
動可能であり、天板１２上の被検体ＭとＸ線透視装置２
との相対的な位置関係を被検体Ｍの体軸方向に変位する
ことができる。モータ１３の駆動制御は、制御部５によ
り行われる。
【００１９】Ｘ線透視装置２は、Ｘ線管２１、撮像系２
２を支持するＣ型アーム２３がベッド１の近傍に定置さ
れた装置基台２４の上部に支持されて構成されている。
Ｃ型アーム２３は、モータ２５の駆動で図２の矢印方向
に変位可能に装置基台２４に支持されており、Ｘ線管２
１、撮像系２２を被検体Ｍの体軸回りに変位可能に構成
し、Ｘ線透過像の撮像方向の調整が可能となっている。
モータ２５の駆動制御は制御部５により行われる。
２を支持するＣ型アーム２３がベッド１の近傍に定置さ
れた装置基台２４の上部に支持されて構成されている。
Ｃ型アーム２３は、モータ２５の駆動で図２の矢印方向
に変位可能に装置基台２４に支持されており、Ｘ線管２
１、撮像系２２を被検体Ｍの体軸回りに変位可能に構成
し、Ｘ線透過像の撮像方向の調整が可能となっている。
モータ２５の駆動制御は制御部５により行われる。
【００２０】Ｘ線管２１と撮像系２２とはＣ型アーム２
３の両端部に取り付けられており、天板１２上の被検体
Ｍを挟み込んだ状態で対向配置されている。Ｘ線管２１
から被検体Ｍの任意の撮像部位に向けて照射され、被検
体Ｍを透過したＸ線は、撮像系２２で受像され、その部
位のＸ線透過像が撮像される。Ｘ線管２１からのＸ線の
照射は、Ｘ線高電圧発生装置２６から所定の電力（Ｘ線
管電圧およびＸ線管電流）がＸ線管２１に供給されて行
われる。Ｘ線高電圧発生装置２６からＸ線管２１への所
定の電力の供給は制御部５に制御されて行われる。撮像
系２２は、イメージインテンシファイアやテレビカメラ
などで構成されている。撮像されたＸ線透過像は画像処
理部３に与えられる。
３の両端部に取り付けられており、天板１２上の被検体
Ｍを挟み込んだ状態で対向配置されている。Ｘ線管２１
から被検体Ｍの任意の撮像部位に向けて照射され、被検
体Ｍを透過したＸ線は、撮像系２２で受像され、その部
位のＸ線透過像が撮像される。Ｘ線管２１からのＸ線の
照射は、Ｘ線高電圧発生装置２６から所定の電力（Ｘ線
管電圧およびＸ線管電流）がＸ線管２１に供給されて行
われる。Ｘ線高電圧発生装置２６からＸ線管２１への所
定の電力の供給は制御部５に制御されて行われる。撮像
系２２は、イメージインテンシファイアやテレビカメラ
などで構成されている。撮像されたＸ線透過像は画像処
理部３に与えられる。
【００２１】画像処理部３は、図３に示すように、デー
タ変換手段としてのＡ／Ｄ（アナログtoディジタルデー
タ）変換器３１、周波数特性変換手段としての周波数特
性変換回路３２、遅延手段としての遅延回路３３、演算
手段としての演算器３４、階調変換回路３５、Ｄ／Ａ
（ディジタルtoアナログ）変換器３６で構成されてい
る。
タ変換手段としてのＡ／Ｄ（アナログtoディジタルデー
タ）変換器３１、周波数特性変換手段としての周波数特
性変換回路３２、遅延手段としての遅延回路３３、演算
手段としての演算器３４、階調変換回路３５、Ｄ／Ａ
（ディジタルtoアナログ）変換器３６で構成されてい
る。
【００２２】造影剤が投与された被検体Ｍの所定の撮像
部位のＸ線透過像が撮像されると、撮像系２２からのそ
の像の画像信号（アナログ信号）は、Ａ／Ｄ変換器３１
でディジタルデータに変換され基本画像が得られる。こ
の基本画像は、骨格などの低周波数成分や、造影剤が投
与された血管像などの高周波数成分を含んだ画像であ
り、これをライブ像として用いる。この基本画像（ライ
ブ像）は、周波数特性変換回路３２と遅延回路３３とに
与えられる。周波数特性変換回路３２は、後述する処理
によって基本画像から血管像などの高周波数成分を除去
しマスク像を得る。演算器３４では、遅延回路３３を経
て供給されるライブ像と、周波数特性変換回路３２を経
て供給されるマスク像とのサブトラクションを行いサブ
トラクション像を求めて階調変換回路３５に与える。階
調変換回路３５では、サブトラクション像をモニタ５に
表示したとき見やすい画像にするために、サブトラクシ
ョン像を構成する各画素の濃度を調整（全画素を対象
に、各画素の濃度に所定濃度を加算したり減算する）す
る。階調変換されたサブトラクション像はＤ／Ａ変換器
３６に与えられ、そこでＤ／Ａ変換されてモニタ４に表
示される。遅延回路３３は、周波数特性変換回路３２で
の処理時間による時間的な遅れを補償するために設けて
おり、これにより、ライブ像とマスク像とが同期されて
演算器３４に供給される。なお、画像処理部３を構成す
る各部の動作制御は、制御部５により行われる。
部位のＸ線透過像が撮像されると、撮像系２２からのそ
の像の画像信号（アナログ信号）は、Ａ／Ｄ変換器３１
でディジタルデータに変換され基本画像が得られる。こ
の基本画像は、骨格などの低周波数成分や、造影剤が投
与された血管像などの高周波数成分を含んだ画像であ
り、これをライブ像として用いる。この基本画像（ライ
ブ像）は、周波数特性変換回路３２と遅延回路３３とに
与えられる。周波数特性変換回路３２は、後述する処理
によって基本画像から血管像などの高周波数成分を除去
しマスク像を得る。演算器３４では、遅延回路３３を経
て供給されるライブ像と、周波数特性変換回路３２を経
て供給されるマスク像とのサブトラクションを行いサブ
トラクション像を求めて階調変換回路３５に与える。階
調変換回路３５では、サブトラクション像をモニタ５に
表示したとき見やすい画像にするために、サブトラクシ
ョン像を構成する各画素の濃度を調整（全画素を対象
に、各画素の濃度に所定濃度を加算したり減算する）す
る。階調変換されたサブトラクション像はＤ／Ａ変換器
３６に与えられ、そこでＤ／Ａ変換されてモニタ４に表
示される。遅延回路３３は、周波数特性変換回路３２で
の処理時間による時間的な遅れを補償するために設けて
おり、これにより、ライブ像とマスク像とが同期されて
演算器３４に供給される。なお、画像処理部３を構成す
る各部の動作制御は、制御部５により行われる。
【００２３】上述したように、周波数特性変換回路３２
は、基本画像（ライブ像）から血管像などの高周波数成
分を除去してマスク像を得ることを目的とする。この処
理方法としては、例えば、空間／周波数変換処理で基本
画像を構成する周波数成分を取り出し、所定のしきい値
周波数以上の周波数成分（高周波数成分）を除去し、こ
れに周波数／空間変換処理を施してマスク像を得る方法
を採用できる。空間／周波数変換処理としては、ＦＦＴ
（高速フーリエ変換）、カルーネン・レーベ変換、ＤＣ
Ｔ（離散コサイン変換）、アダマール変換などが知られ
ており、この実施例にもこれら各変換方式を採用するこ
とができる。また、周波数／空間変換処理は、上記空間
／周波数変換処理の逆変換（逆ＦＦＴ、逆カルーネン・
レーベ変換、逆ＤＣＴ、逆アダマール変換など）で実現
できる。さらに、所定のしきい値周波数は、血管像など
を好適に除去し得る周波数を実験的に求めておけばよ
い。
は、基本画像（ライブ像）から血管像などの高周波数成
分を除去してマスク像を得ることを目的とする。この処
理方法としては、例えば、空間／周波数変換処理で基本
画像を構成する周波数成分を取り出し、所定のしきい値
周波数以上の周波数成分（高周波数成分）を除去し、こ
れに周波数／空間変換処理を施してマスク像を得る方法
を採用できる。空間／周波数変換処理としては、ＦＦＴ
（高速フーリエ変換）、カルーネン・レーベ変換、ＤＣ
Ｔ（離散コサイン変換）、アダマール変換などが知られ
ており、この実施例にもこれら各変換方式を採用するこ
とができる。また、周波数／空間変換処理は、上記空間
／周波数変換処理の逆変換（逆ＦＦＴ、逆カルーネン・
レーベ変換、逆ＤＣＴ、逆アダマール変換など）で実現
できる。さらに、所定のしきい値周波数は、血管像など
を好適に除去し得る周波数を実験的に求めておけばよ
い。
【００２４】なお、例えば、ＤＣＴ、逆ＤＣＴ処理を高
速（１／３０秒程度）で実現する専用チップが市販され
ており、このような専用チップを用いれば、基本画像を
撮像してから極めて短時間で（リアルタイムに近く）サ
ブトラクション像をモニタ４に表示させることができ
る。
速（１／３０秒程度）で実現する専用チップが市販され
ており、このような専用チップを用いれば、基本画像を
撮像してから極めて短時間で（リアルタイムに近く）サ
ブトラクション像をモニタ４に表示させることができ
る。
【００２５】制御部５は、操作盤６からの指示などによ
って各装置、各部の駆動制御や動作制御を行う。この制
御部５は、例えば、後述する動作を実現するプログラム
を遂行するＣＰＵ（中央処理装置）で構成されている。
って各装置、各部の駆動制御や動作制御を行う。この制
御部５は、例えば、後述する動作を実現するプログラム
を遂行するＣＰＵ（中央処理装置）で構成されている。
【００２６】操作盤６は、撮像部位や条件の設定、処理
開始指示などを操作者が行うためのものである。
開始指示などを操作者が行うためのものである。
【００２７】上記構成を有する実施例装置の動作を以下
に説明する。まず、被検体Ｍのある１箇所の撮像部位
（例えば胸部）のサブトラクション像を得る場合の動作
を説明する。
に説明する。まず、被検体Ｍのある１箇所の撮像部位
（例えば胸部）のサブトラクション像を得る場合の動作
を説明する。
【００２８】この場合、まず、操作者により操作盤６か
ら設定された撮像部位や条件（撮像方向など）に従っ
て、制御部５は、モータ１３を駆動制御して被検体Ｍを
載置した天板１２を水平移動させ、設定された撮像部位
（胸部とする）を、Ｘ線管２１、撮像系２２の間の撮像
位置に位置させ、モータ２５を駆動制御してＸ線管２
１、撮像系２２を被検体Ｍ（の撮像部位）の体軸回りに
変位させ、撮像方向を調節する。この状態を図４に示
す。図４では、撮像部位ＳＢの下方からＸ線を照射して
Ｘ線透過像を撮像するように撮像方向が調節されてい
る。
ら設定された撮像部位や条件（撮像方向など）に従っ
て、制御部５は、モータ１３を駆動制御して被検体Ｍを
載置した天板１２を水平移動させ、設定された撮像部位
（胸部とする）を、Ｘ線管２１、撮像系２２の間の撮像
位置に位置させ、モータ２５を駆動制御してＸ線管２
１、撮像系２２を被検体Ｍ（の撮像部位）の体軸回りに
変位させ、撮像方向を調節する。この状態を図４に示
す。図４では、撮像部位ＳＢの下方からＸ線を照射して
Ｘ線透過像を撮像するように撮像方向が調節されてい
る。
【００２９】次に、被検体Ｍに造影剤を投与する。な
お、造影剤を投与してから上記位置合わせ動作などを行
ってもよい。いずれにしても、以下の撮像動作の前に、
被検体Ｍに造影剤を投与しておき、撮像部位ＳＢに造影
剤が拡散した状態で操作者が操作盤６から処理開始を指
示し、以下の撮像動作が実行される。
お、造影剤を投与してから上記位置合わせ動作などを行
ってもよい。いずれにしても、以下の撮像動作の前に、
被検体Ｍに造影剤を投与しておき、撮像部位ＳＢに造影
剤が拡散した状態で操作者が操作盤６から処理開始を指
示し、以下の撮像動作が実行される。
【００３０】処理開始が指示されると、制御部５はＸ線
高電圧発生装置２６を制御して、Ｘ線管２１に所定の電
力を供給させてＸ線を照射させ、造影剤が拡散された撮
像部位ＳＢのＸ線透過像を撮像させる。そして、制御部
５は、画像処理部３の各部を制御して、基本画像（ライ
ブ像）を得るとともに、その基本画像からマスク像を求
め、ライブ像とマスク像とのサブトラクションを行わせ
サブトラクション像をモニタ４に表示させる。
高電圧発生装置２６を制御して、Ｘ線管２１に所定の電
力を供給させてＸ線を照射させ、造影剤が拡散された撮
像部位ＳＢのＸ線透過像を撮像させる。そして、制御部
５は、画像処理部３の各部を制御して、基本画像（ライ
ブ像）を得るとともに、その基本画像からマスク像を求
め、ライブ像とマスク像とのサブトラクションを行わせ
サブトラクション像をモニタ４に表示させる。
【００３１】上述したように、この実施例によれば、あ
る撮像部位のサブトラクション像を得るための被検体へ
のＸ線照射は１回でよく、被検体へのＸ線曝射線量を従
来に比べて半分に減らせることができる。また、１フレ
ームの基本画像からマスク像とライブ像を得ているの
で、被検体の体動によるマスク像とライブ像の画像のず
れも完全に無くすことができる。
る撮像部位のサブトラクション像を得るための被検体へ
のＸ線照射は１回でよく、被検体へのＸ線曝射線量を従
来に比べて半分に減らせることができる。また、１フレ
ームの基本画像からマスク像とライブ像を得ているの
で、被検体の体動によるマスク像とライブ像の画像のず
れも完全に無くすことができる。
【００３２】次に、被検体ＭとＸ線透視装置２との相対
的な位置関係を被検体Ｍの体軸方向に変位させて、例え
ば、図５に示すように、胸部から腹部にかけての領域Ｓ
Ｒ内の複数の撮像部位のサブトラクション像を得る場合
の動作を説明する。なお、この実施例では、Ｘ線透視装
置２を固定し、これに対して被検体Ｍを載置した天板１
２を水平移動するように構成しているが、撮像状況をわ
かり易くするために、図５では、天板１２上の被検体Ｍ
を固定し、これに対してＸ線透視装置２（Ｘ線管２１、
撮像系２２）が変位しているように描いている。
的な位置関係を被検体Ｍの体軸方向に変位させて、例え
ば、図５に示すように、胸部から腹部にかけての領域Ｓ
Ｒ内の複数の撮像部位のサブトラクション像を得る場合
の動作を説明する。なお、この実施例では、Ｘ線透視装
置２を固定し、これに対して被検体Ｍを載置した天板１
２を水平移動するように構成しているが、撮像状況をわ
かり易くするために、図５では、天板１２上の被検体Ｍ
を固定し、これに対してＸ線透視装置２（Ｘ線管２１、
撮像系２２）が変位しているように描いている。
【００３３】この場合、制御部５は、最初の撮像部位
（図では、撮像領域ＳＲの左端部側の撮像部位）を撮像
位置に位置させ、撮像方向を調節する。そして、以下の
撮像動作の前に被検体Ｍに造影剤を投与しておく。
（図では、撮像領域ＳＲの左端部側の撮像部位）を撮像
位置に位置させ、撮像方向を調節する。そして、以下の
撮像動作の前に被検体Ｍに造影剤を投与しておく。
【００３４】被検体Ｍの各撮像部位（領域ＳＲ）に造影
剤が拡散し、処理開始が指示されると、上記１箇所の撮
像領域ＳＢのサブトラクション像を求めた手順と同様の
手順で、最初の撮像部位のサブトラクション像を求め、
天板１２を図５の左方向に定速で移動させながら、以降
の各撮像部位が撮像位置に位置するごとに、その撮像部
位のサブトラクション像を順次求めていく。
剤が拡散し、処理開始が指示されると、上記１箇所の撮
像領域ＳＢのサブトラクション像を求めた手順と同様の
手順で、最初の撮像部位のサブトラクション像を求め、
天板１２を図５の左方向に定速で移動させながら、以降
の各撮像部位が撮像位置に位置するごとに、その撮像部
位のサブトラクション像を順次求めていく。
【００３５】このように、被検体とＸ線透視装置との相
対的な位置関係を被検体の体軸方向に変位させて、複数
の連続する撮像部位に対するサブトラクション像を得る
場合であっても、被検体へのＸ線曝射線量を従来の半分
にでき、被検体の体動や撮像する部位のずれなどに起因
する、各撮像部位ごとの対となるマスク像とライブ像と
の画像のずれを完全に無くすことができる。
対的な位置関係を被検体の体軸方向に変位させて、複数
の連続する撮像部位に対するサブトラクション像を得る
場合であっても、被検体へのＸ線曝射線量を従来の半分
にでき、被検体の体動や撮像する部位のずれなどに起因
する、各撮像部位ごとの対となるマスク像とライブ像と
の画像のずれを完全に無くすことができる。
【００３６】なお、天板１２を固定し、Ｘ線透視装置２
を天板１２上の被検体Ｍの体軸方向に移動させること
で、被検体とＸ線透視装置との相対的な位置関係を被検
体の体軸方向に変位させるように構成してもよい。
を天板１２上の被検体Ｍの体軸方向に移動させること
で、被検体とＸ線透視装置との相対的な位置関係を被検
体の体軸方向に変位させるように構成してもよい。
【００３７】また、被検体Ｍのある部位（例えば、胸
部）を撮像位置に位置させた状態で、図６に示すよう
に、Ｘ線管２１、撮像系２２をその部位の回り（体軸回
り）に回転変位させながら、各撮像方向からのサブトラ
クション像を求めることもあるが、このような場合であ
っても、上記各動作と同様に、被検体ＭへのＸ線曝射線
量を従来の半分に減らせるし、各撮像方向からのサブト
ラクション像を求めるための一対のマスク像とライブ像
とに画像のずれが生じることがない。
部）を撮像位置に位置させた状態で、図６に示すよう
に、Ｘ線管２１、撮像系２２をその部位の回り（体軸回
り）に回転変位させながら、各撮像方向からのサブトラ
クション像を求めることもあるが、このような場合であ
っても、上記各動作と同様に、被検体ＭへのＸ線曝射線
量を従来の半分に減らせるし、各撮像方向からのサブト
ラクション像を求めるための一対のマスク像とライブ像
とに画像のずれが生じることがない。
【００３８】 ＜参考例＞ 図７は、参考例に係る装置に備えられる画像処理部の構
成を示すブロック図である。この参考例では、基本画像
（ライブ像）を第１の記憶装置４１に記憶するととも
に、周波数特性変換回路３２で求めたマスク像を第２の
記憶装置４２に記憶させるように構成している。これに
より、第１、第２の記憶装置４１、４２から同期させて
ライブ像、マスク像を読み出し演算器３４に与えれば、
各画像を同期させて演算器３４に供給することができる
ので、遅延回路３３を設けなくてもよい。
成を示すブロック図である。この参考例では、基本画像
（ライブ像）を第１の記憶装置４１に記憶するととも
に、周波数特性変換回路３２で求めたマスク像を第２の
記憶装置４２に記憶させるように構成している。これに
より、第１、第２の記憶装置４１、４２から同期させて
ライブ像、マスク像を読み出し演算器３４に与えれば、
各画像を同期させて演算器３４に供給することができる
ので、遅延回路３３を設けなくてもよい。
【００３９】また、複数の撮像部位のサブトラクション
像を得る場合や、所定の部位に対して複数の撮像方向か
らのサブトラクション像を得る場合などにおいて、各撮
像部位や各撮像方向からの複数のライブ像を第１の記憶
装置４１に記憶するとともに、各撮像部位や各撮像方向
からの複数のマスク像を第２の記憶装置４２に記憶する
ことができ、例えば、一連の撮像が完了した後に、所望
の撮像部位や撮像方向からのサブトラクション像を選択
的にモニタ４に表示させることもできる。
像を得る場合や、所定の部位に対して複数の撮像方向か
らのサブトラクション像を得る場合などにおいて、各撮
像部位や各撮像方向からの複数のライブ像を第１の記憶
装置４１に記憶するとともに、各撮像部位や各撮像方向
からの複数のマスク像を第２の記憶装置４２に記憶する
ことができ、例えば、一連の撮像が完了した後に、所望
の撮像部位や撮像方向からのサブトラクション像を選択
的にモニタ４に表示させることもできる。
【００４０】さらに、第１、第２の記憶装置４１、４２
は揮発性のメモリで構成してもよいが、例えば、磁気デ
ィスク装置などの保存可能な記憶装置で構成すれば、サ
ブトラクション像を事後的にモニタ４に表示させること
もできる。
は揮発性のメモリで構成してもよいが、例えば、磁気デ
ィスク装置などの保存可能な記憶装置で構成すれば、サ
ブトラクション像を事後的にモニタ４に表示させること
もできる。
【００４１】また、例えば、図８のように構成すれば、
マスク像やライブ像をモニタ４に表示することができ、
さらに、図８の点線のようにモニタ４などを３台接続す
れば、マスク像、ライブ像、サブトラクション像を各モ
ニタ４に並行して表示させることもできる。
マスク像やライブ像をモニタ４に表示することができ、
さらに、図８の点線のようにモニタ４などを３台接続す
れば、マスク像、ライブ像、サブトラクション像を各モ
ニタ４に並行して表示させることもできる。
【００４２】また、図９に示すように、サブトラクショ
ン像を第３のメモリ４３に記憶させるように構成すれ
ば、サブトラクション像を再演算することなく以後何度
でもモニタ４に再表示させることができる。さらに、こ
の第３のメモリ４３を保存可能な記憶装置で構成すれ
ば、サブトラクション像を再演算することなく事後的に
モニタ４に表示させることもできる。
ン像を第３のメモリ４３に記憶させるように構成すれ
ば、サブトラクション像を再演算することなく以後何度
でもモニタ４に再表示させることができる。さらに、こ
の第３のメモリ４３を保存可能な記憶装置で構成すれ
ば、サブトラクション像を再演算することなく事後的に
モニタ４に表示させることもできる。
【００４３】この参考例およびそれに対する各変形例
は、以下の第２、第３実施例にも同様に適用することが
できる。
は、以下の第２、第３実施例にも同様に適用することが
できる。
【００４４】 ＜第２実施例＞ 図１０は、第２実施例装置に備えられる画像処理部の構
成を示すブロック図である。この第３実施例は、基本画
像（ライブ像）の高周波数成分を強調するための第２の
周波数特性変換手段としての第２の周波数特性変換回路
５１を付設したことを特徴とする。
成を示すブロック図である。この第３実施例は、基本画
像（ライブ像）の高周波数成分を強調するための第２の
周波数特性変換手段としての第２の周波数特性変換回路
５１を付設したことを特徴とする。
【００４５】基本画像の高周波数成分を強調する処理
は、例えば、空間／周波数変換して得られた各基本画像
の周波数成分のうちの高周波数成分を増分させ、それを
周波数／空間変換することで実現できる。第２の周波数
特性変換回路５１ではこのような処理を実現するように
構成される。
は、例えば、空間／周波数変換して得られた各基本画像
の周波数成分のうちの高周波数成分を増分させ、それを
周波数／空間変換することで実現できる。第２の周波数
特性変換回路５１ではこのような処理を実現するように
構成される。
【００４６】これにより、血管像がより強調されたサブ
トラクション像を得ることができる。なお、この第２実
施例に係る特徴部分は、以下の第３実施例にも同様に適
用することができる。
トラクション像を得ることができる。なお、この第２実
施例に係る特徴部分は、以下の第３実施例にも同様に適
用することができる。
【００４７】 ＜第３実施例＞ 図１１は、第３実施例装置に備えられる画像処理部の構
成を示すブロック図である。この第３実施例は、複数回
分の撮像で得られた基本画像の平均画像を求める積分処
理手段としての積分処理部６１を付設したことを特徴と
する。
成を示すブロック図である。この第３実施例は、複数回
分の撮像で得られた基本画像の平均画像を求める積分処
理手段としての積分処理部６１を付設したことを特徴と
する。
【００４８】積分処理部６１は、図１２（ａ）に示すよ
うに、第２の階調変換回路７１と加算器７２とメモリ７
３とで構成してもよいし、図１２（ｂ）に示すように加
算器７２とメモリ７３と除算器７４とで構成してもよ
い。
うに、第２の階調変換回路７１と加算器７２とメモリ７
３とで構成してもよいし、図１２（ｂ）に示すように加
算器７２とメモリ７３と除算器７４とで構成してもよ
い。
【００４９】例えば、Ｎ回分の撮像で得られた基本画像
の平均画像を求める場合、図１２（ａ）の構成では、順
次与えられる基本画像の濃度（基本画像を構成する各画
素の濃度）を第２の階調変換回路７１で１／Ｎにするよ
うに階調変換させ、加算器７２で、メモリ７３に記憶さ
れている加算画像に加算していく。ただし、最初の基本
画像が与えられるときには、メモリ７３には何も記憶さ
れていない。
の平均画像を求める場合、図１２（ａ）の構成では、順
次与えられる基本画像の濃度（基本画像を構成する各画
素の濃度）を第２の階調変換回路７１で１／Ｎにするよ
うに階調変換させ、加算器７２で、メモリ７３に記憶さ
れている加算画像に加算していく。ただし、最初の基本
画像が与えられるときには、メモリ７３には何も記憶さ
れていない。
【００５０】また、図１２（ｂ）の構成では、加算器７
２で、順次与えられる基本画像を、メモリ７３に記憶さ
れている加算画像（最初の基本画像が与えられるときに
は、メモリ７３には何も記憶されていない）に順次加算
していき、最後に、メモリ７３に記憶されたＮ回分の基
本画像の加算結果を除算器７４で１／Ｎにする。
２で、順次与えられる基本画像を、メモリ７３に記憶さ
れている加算画像（最初の基本画像が与えられるときに
は、メモリ７３には何も記憶されていない）に順次加算
していき、最後に、メモリ７３に記憶されたＮ回分の基
本画像の加算結果を除算器７４で１／Ｎにする。
【００５１】このように、複数回分の撮像で得られた基
本画像の平均画像を求めることにより、基本画像のＳ／
Ｎ比を向上させることができ、この発明では、基本画像
からマスク像とライブ像を求めるので、マスク像とライ
ブ像のＳ／Ｎ比をまとめて向上させることができる。
本画像の平均画像を求めることにより、基本画像のＳ／
Ｎ比を向上させることができ、この発明では、基本画像
からマスク像とライブ像を求めるので、マスク像とライ
ブ像のＳ／Ｎ比をまとめて向上させることができる。
【００５２】なお、例えば、ある撮像部位に対する１撮
像方向からのサブトラクション像を得る場合には、その
撮像部位に１方向からのみＸ線照射して基本画像を撮像
する動作をＮ回繰り返すことになる。従来装置において
も、マスク像やライブ像のＳ／Ｎ比を向上させるため
に、マスク像とライブ像をそれぞれ複数（Ｎ）回分撮像
し、各画像の平均画像を求めることがあるが、この場
合、被検体へのＸ線照射は、マスク像とライブ像を別々
に撮像するので、２×Ｎ回行うことになる。これに対し
てこの実施例では被検体へのＸ線照射は従来装置の半分
でよい。
像方向からのサブトラクション像を得る場合には、その
撮像部位に１方向からのみＸ線照射して基本画像を撮像
する動作をＮ回繰り返すことになる。従来装置において
も、マスク像やライブ像のＳ／Ｎ比を向上させるため
に、マスク像とライブ像をそれぞれ複数（Ｎ）回分撮像
し、各画像の平均画像を求めることがあるが、この場
合、被検体へのＸ線照射は、マスク像とライブ像を別々
に撮像するので、２×Ｎ回行うことになる。これに対し
てこの実施例では被検体へのＸ線照射は従来装置の半分
でよい。
【００５３】また、被検体とＸ線透視装置との相対的な
位置関係を被検体の体軸方向に変位させながら、複数の
撮像部位のサブトラクション像を得る場合には、対象と
なる撮像部位の相前後する複数の撮像部位で撮像された
複数回分の基本画像からその対象となる撮像部位の平均
画像を求めるようにしてもよい。ある部位の周回方向に
Ｘ線管、撮像系を回転変位させながら複数の撮像方向か
らのサブトラクション像を得る場合も同様に、対象とな
る撮像方向の相前後する複数の撮像方向から撮像された
複数回分の基本画像からその対象となる撮像方向からの
平均画像を求めるようにしてもよい。これらの場合で
も、被検体へのＸ線曝射線量は従来の場合（従来装置で
同様の手順で平均画像を求める場合）の半分になる。
位置関係を被検体の体軸方向に変位させながら、複数の
撮像部位のサブトラクション像を得る場合には、対象と
なる撮像部位の相前後する複数の撮像部位で撮像された
複数回分の基本画像からその対象となる撮像部位の平均
画像を求めるようにしてもよい。ある部位の周回方向に
Ｘ線管、撮像系を回転変位させながら複数の撮像方向か
らのサブトラクション像を得る場合も同様に、対象とな
る撮像方向の相前後する複数の撮像方向から撮像された
複数回分の基本画像からその対象となる撮像方向からの
平均画像を求めるようにしてもよい。これらの場合で
も、被検体へのＸ線曝射線量は従来の場合（従来装置で
同様の手順で平均画像を求める場合）の半分になる。
【００５４】
【発明の効果】以上の説明から明らかなように、請求項
１ないし３に記載の発明によれば、１回のＸ線照射によ
り得たＸ線透過像からマスク像とライブ像を得てサブト
ラクション像を求めるので、サブトラクション像を得る
ための被検体へのＸ線曝射線量を従来の半分に軽減する
ことができる。
１ないし３に記載の発明によれば、１回のＸ線照射によ
り得たＸ線透過像からマスク像とライブ像を得てサブト
ラクション像を求めるので、サブトラクション像を得る
ための被検体へのＸ線曝射線量を従来の半分に軽減する
ことができる。
【００５５】また、１回で得たＸ線透過像からマスク像
とライブ像を得ているので、これらマスク像とライブ像
との画像のずれを完全に無くすことができ、このような
画像のずれに起因するサブトラクション像のアーティフ
ァクトを無くすことができる。
とライブ像を得ているので、これらマスク像とライブ像
との画像のずれを完全に無くすことができ、このような
画像のずれに起因するサブトラクション像のアーティフ
ァクトを無くすことができる。
【００５６】さらに、マスク像（高周波数除去画像）
を、基本画像を構成する周波数成分を空間／周波数変換
処理で取り出し、所定のしきい値周波数以上の周波数成
分を除去し、これに周波数／空間変換処理を施すことで
得ているので、適正なマスク像を得ることができる。
を、基本画像を構成する周波数成分を空間／周波数変換
処理で取り出し、所定のしきい値周波数以上の周波数成
分を除去し、これに周波数／空間変換処理を施すことで
得ているので、適正なマスク像を得ることができる。
【００５７】特に、請求項１に記載の発明によれば、演
算手段への基本画像の供給を遅らせる遅延手段を設けて
いるので、メモリなどを用いずに、周波数特性変換手段
での処理時間による時間的な遅れを補償して、基本画像
（ライブ像）と高周波数除去画像（マスク像）とを演算
手段に同期して供給させることができる。
算手段への基本画像の供給を遅らせる遅延手段を設けて
いるので、メモリなどを用いずに、周波数特性変換手段
での処理時間による時間的な遅れを補償して、基本画像
（ライブ像）と高周波数除去画像（マスク像）とを演算
手段に同期して供給させることができる。
【００５８】また、請求項２に記載の発明によれば、基
本画像の高周波数成分を強調した高周波数強調画像を得
て、この高周波数強調画像と高周波数除去画像とのサブ
トラクションを行ってサブトラクション像を求めるの
で、血管像などがより強調されたサブトラクション像を
得ることができる。
本画像の高周波数成分を強調した高周波数強調画像を得
て、この高周波数強調画像と高周波数除去画像とのサブ
トラクションを行ってサブトラクション像を求めるの
で、血管像などがより強調されたサブトラクション像を
得ることができる。
【００５９】また、請求項３に記載の発明によれば、複
数回分の撮像で得られた基本画像の平均画像を求め、こ
の平均画像をライブ像として用いるとともに、このライ
ブ像（平均画像）から高周波数除去画像を得て、ライブ
像（平均画像）からマスク像（平均画像から得られた高
周波数除去画像）をサブトラクションしてサブトラクシ
ョン像を求めるので、ライブ像及びマスク像のＳ／Ｎ比
を向上させることができ、画質の良いサブトラクション
像を得ることができる。
数回分の撮像で得られた基本画像の平均画像を求め、こ
の平均画像をライブ像として用いるとともに、このライ
ブ像（平均画像）から高周波数除去画像を得て、ライブ
像（平均画像）からマスク像（平均画像から得られた高
周波数除去画像）をサブトラクションしてサブトラクシ
ョン像を求めるので、ライブ像及びマスク像のＳ／Ｎ比
を向上させることができ、画質の良いサブトラクション
像を得ることができる。
【図１】この発明の第１実施例に係るディジタルアンギ
オグラフィ装置の全体構成を示す正面図である。
オグラフィ装置の全体構成を示す正面図である。
【図２】Ｘ線透視装置を側面から見た図である。
【図３】第１実施例装置に備えられた画像処理部の構成
を示すブロック図である。
を示すブロック図である。
【図４】ある１箇所の撮像部位のサブトラクション像を
得る場合の動作を説明するための図である。
得る場合の動作を説明するための図である。
【図５】被検体の体軸方向の複数の撮像部位のサブトラ
クション像を得る場合の動作を説明するための図であ
る。
クション像を得る場合の動作を説明するための図であ
る。
【図６】ある部位に対する複数の撮像方向からのサブト
ラクション像を得る場合の動作を説明するための図であ
る。
ラクション像を得る場合の動作を説明するための図であ
る。
【図７】参考例に係る装置に備えられる画像処理部の構
成を示すブロック図である。
成を示すブロック図である。
【図８】参考例に係る装置の変形例の構成を示すブロッ
ク図である。
ク図である。
【図９】参考例に係る装置の別の変形例の構成を示すブ
ロック図である。
ロック図である。
【図１０】第２実施例装置に備えられる画像処理部の構
成を示すブロック図である。
成を示すブロック図である。
【図１１】第３実施例装置に備えられる画像処理部の構
成を示すブロック図である。
成を示すブロック図である。
【図１２】第３実施例の積分処理部の構成を示すブロッ
ク図である。
ク図である。
２ … Ｘ線透視装置 ３ … 画像処理部 ３１ … Ａ／Ｄ変換器 ３２ … 周波数特性変換回路 ３３ … 遅延回路 ３４ … 演算器 ５１ … 第２の周波数特性変換回路 ６１ … 積分処理部 Ｍ … 被検体 ＳＢ … 撮像部位
───────────────────────────────────────────────────── フロントページの続き (58)調査した分野(Int.Cl.6，ＤＢ名) A61B 6/00 - 6/14
Claims (3)
- 【請求項１】 被検体の所定の撮像部位のサブトラクシ
ョン像を得るためのディジタルアンギオグラフィ装置で
あって、（１−ａ）前記撮像部位にＸ線を照射し、その
撮像部位のＸ線透過像を撮像するＸ線透視手段と、（１
−ｂ）前記Ｘ線透過像をディジタルデータに変換するデ
ータ変換手段と、（１−ｃ）造影剤が投与された被検体
の所定の撮像部位を撮像し、ディジタルデータに変換さ
れたＸ線透過像（以下、基本画像という）を構成する周
波数成分を空間／周波数変換処理で取り出し、所定のし
きい値周波数以上の周波数成分を除去し、これに周波数
／空間変換処理を施して、前記基本画像から前記しきい
値周波数以上の高周波数成分を除去した画像（以下、高
周波数除去画像という）を得る周波数特性変換手段と、
（１−ｄ）前記基本画像と前記高周波数除去画像とのサ
ブトラクションを行いその撮像部位のサブトラクション
像を求める演算手段と、（１−ｅ）前記演算手段への基
本画像の供給を遅らせて、前記基本画像と前記高周波数
除去画像とを前記演算手段に同期して供給させる遅延手
段とを備えたことを特徴とするディジタルアンギオグラ
フィ装置。 - 【請求項２】 被検体の所定の撮像部位のサブトラクシ
ョン像を得るためのディジタルアンギオグラフィ装置で
あって、（２−ａ）前記撮像部位にＸ線を照射し、その
撮像部位のＸ線透過像を撮像するＸ線透視手段と、（２
−ｂ）前記Ｘ線透過像をディジタルデータに変換するデ
ータ変換手段と、（２−ｃ）造影剤が投与された被検体
の所定の撮像部位を撮像し、ディジタルデータに変換さ
れたＸ線透過像（以下、基本画像という）を構成する周
波数成分を空間／周波数変換処理で取り出し、所定のし
きい値周波数以上の周波数成分を除去し、これに周波数
／空間変換処理を施して、前記基本画像から前記しきい
値周波数以上の高周波数成分を除去した画像（以下、高
周波数除去画像という）を得る周波数特性変換手段と、
（２−ｄ）前記基本画像の高周波数成分を強調した画像
（以下、高周波数強調画像という）を得る第２の周波数
特性変換手段と、（２−ｅ）前記高周波数強調画像と前
記高周波数除去画像とのサブトラクションを行いその撮
像部位のサブトラクション像を求める演算手段とを備え
たことを特徴とするディジタルアンギオグラフィ装置。 - 【請求項３】 被検体の所定の撮像部位のサブトラクシ
ョン像を得るためのディジタルアンギオグラフィ装置で
あって、（３−ａ）前記撮像部位にＸ線を照射し、その
撮像部位のＸ線透過像を撮像するＸ線透視手段と、（３
−ｂ）前記Ｘ線透過像をディジタルデータに変換するデ
ータ変換手段と、（３−ｃ）造影剤が投与された被検体
の所定の撮像部位を撮像し、ディジタルデータに変換さ
れたＸ線透過像（以下、基本画像という）を順次取り込
み、複数回分の撮像で得られた基本画像の平均画像を求
める積分処理手段と、（３−ｄ）前記平均画像を構成す
る周波数成分を空間／周波数変換処理で取り出し、所定
のしきい値周波数以上の周波数成分を除去し、これに周
波数／空間変換処理を施して、前記平均画像から前記し
きい値周波数以上の高周波数成分を除去した画像（以
下、高周波数除去画像という）を得る周波数特性変換手
段と、（３−ｄ）前記平均画像と前記高周波数除去画像
とのサブトラクションを行いその撮像部位のサブトラク
ション像を求める演算手段とを備えたことを特徴とする
ディジタルアンギオグラフィ装置。
Priority Applications (7)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP7148191A JP2927209B2 (ja) | 1995-05-22 | 1995-05-22 | ディジタルアンギオグラフィ装置 |
| US08/642,871 US5671743A (en) | 1995-05-22 | 1996-05-06 | Digital angiographic apparatus |
| EP96107248A EP0744157B1 (en) | 1995-05-22 | 1996-05-08 | Digital angiographic apparatus |
| DE69629149T DE69629149T2 (de) | 1995-05-22 | 1996-05-08 | Digitales angiographisches Gerät |
| AU52363/96A AU704755B2 (en) | 1995-05-22 | 1996-05-17 | Digital angiographic apparatus |
| KR1019960016778A KR100413919B1 (ko) | 1995-05-22 | 1996-05-18 | 디지탈혈관조영(angiography)장치 |
| CNB961066695A CN1256917C (zh) | 1995-05-22 | 1996-05-22 | 数字血管造影设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP7148191A JP2927209B2 (ja) | 1995-05-22 | 1995-05-22 | ディジタルアンギオグラフィ装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH08308823A JPH08308823A (ja) | 1996-11-26 |
| JP2927209B2 true JP2927209B2 (ja) | 1999-07-28 |
Family
ID=15447293
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP7148191A Expired - Fee Related JP2927209B2 (ja) | 1995-05-22 | 1995-05-22 | ディジタルアンギオグラフィ装置 |
Country Status (7)
| Country | Link |
|---|---|
| US (1) | US5671743A (ja) |
| EP (1) | EP0744157B1 (ja) |
| JP (1) | JP2927209B2 (ja) |
| KR (1) | KR100413919B1 (ja) |
| CN (1) | CN1256917C (ja) |
| AU (1) | AU704755B2 (ja) |
| DE (1) | DE69629149T2 (ja) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002354341A (ja) * | 2001-05-24 | 2002-12-06 | Shimadzu Corp | ディジタルサブトラクション装置 |
| JP2005033656A (ja) * | 2003-07-09 | 2005-02-03 | Shimadzu Corp | ディジタルサブトラクション装置 |
| US11179131B2 (en) | 2019-04-03 | 2021-11-23 | Shimadzu Corporation | Radiation fluoroscopic imaging apparatus |
Families Citing this family (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2360684B (en) * | 1997-11-15 | 2001-12-19 | Elekta Ab | Analysis of radiographic images |
| GB9724110D0 (en) | 1997-11-15 | 1998-01-14 | Elekta Ab | Analysis of radiographic images |
| JP4626086B2 (ja) * | 2001-05-09 | 2011-02-02 | 株式会社島津製作所 | ディジタルサブトラクション装置 |
| DE10214763A1 (de) * | 2002-04-03 | 2003-10-30 | Philips Intellectual Property | Verfahren zur Bestimmung eines Bildes aus einer Bildsequenz |
| DE602005011515D1 (de) * | 2004-03-31 | 2009-01-22 | Toshiba Kk | Vorrichtung und Verfahren zur Verarbeitung medizinischer Bilder |
| JP4755184B2 (ja) * | 2004-07-27 | 2011-08-24 | デュール デンタル アクチェンゲゼルシャフト | 放射線写真上の異なる構造の視認可能性を改善するための方法および装置 |
| KR100830198B1 (ko) * | 2006-04-20 | 2008-05-16 | 허감 | 관상동맥 ｃｔ 혈관조영술에서의 ｃｔ번호의 표준편차를이용한 방사선량 조절방법 및 장치 |
| JP2008160305A (ja) * | 2006-12-21 | 2008-07-10 | Toyota Motor Corp | 車載用暗視装置 |
| JP4300433B2 (ja) * | 2007-08-10 | 2009-07-22 | トヨタ自動車株式会社 | レーザ溶接品質評価方法及びその装置 |
| JP2009225979A (ja) * | 2008-03-21 | 2009-10-08 | Japan Health Science Foundation | 医用画像作成装置及び方法 |
| US8355557B2 (en) * | 2008-10-10 | 2013-01-15 | Siemens Aktiengesellschaft | System and method for decomposed temporal filtering for X-ray guided intervention application |
| JP2015226694A (ja) * | 2014-06-02 | 2015-12-17 | 株式会社島津製作所 | Ｘ線撮影装置およびｘ線画像の生成方法 |
| CN106469449B (zh) * | 2016-08-31 | 2019-12-20 | 上海联影医疗科技有限公司 | 一种医学图像中病灶的显示方法及装置 |
| CN106923850A (zh) * | 2017-02-22 | 2017-07-07 | 北京唯迈医疗设备有限公司 | 一种数字减影图像的生成方法和设备 |
| US11793478B2 (en) * | 2017-03-28 | 2023-10-24 | Canon Medical Systems Corporation | Medical image processing apparatus, medical image processing method, and X-ray diagnostic apparatus |
| CN108852385B (zh) * | 2018-03-13 | 2022-03-04 | 中国科学院上海应用物理研究所 | 一种x射线造影方法及基于x射线造影的动态阅片方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS62186381A (ja) * | 1986-02-12 | 1987-08-14 | Hitachi Ltd | 画像位置合わせ方式 |
| JP2536024B2 (ja) * | 1988-02-29 | 1996-09-18 | 株式会社島津製作所 | Ｘ線画像処理装置 |
| JPH02114776A (ja) * | 1988-10-25 | 1990-04-26 | Toshiba Corp | Ｘ線診断装置 |
| DE4417628C1 (de) * | 1994-05-19 | 1995-09-28 | Siemens Ag | Verfahren zur adaptiven Rauschverminderung für digitale Bildsequenzen |
-
1995
- 1995-05-22 JP JP7148191A patent/JP2927209B2/ja not_active Expired - Fee Related
-
1996
- 1996-05-06 US US08/642,871 patent/US5671743A/en not_active Expired - Lifetime
- 1996-05-08 DE DE69629149T patent/DE69629149T2/de not_active Expired - Lifetime
- 1996-05-08 EP EP96107248A patent/EP0744157B1/en not_active Expired - Lifetime
- 1996-05-17 AU AU52363/96A patent/AU704755B2/en not_active Ceased
- 1996-05-18 KR KR1019960016778A patent/KR100413919B1/ko not_active IP Right Cessation
- 1996-05-22 CN CNB961066695A patent/CN1256917C/zh not_active Expired - Lifetime
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2002354341A (ja) * | 2001-05-24 | 2002-12-06 | Shimadzu Corp | ディジタルサブトラクション装置 |
| JP4639526B2 (ja) * | 2001-05-24 | 2011-02-23 | 株式会社島津製作所 | ディジタルサブトラクション装置 |
| JP2005033656A (ja) * | 2003-07-09 | 2005-02-03 | Shimadzu Corp | ディジタルサブトラクション装置 |
| JP4506119B2 (ja) * | 2003-07-09 | 2010-07-21 | 株式会社島津製作所 | ディジタルサブトラクション装置およびディジタルアンギオグラフィ装置並びにディジタルサブトラクション像の表示方法 |
| US11179131B2 (en) | 2019-04-03 | 2021-11-23 | Shimadzu Corporation | Radiation fluoroscopic imaging apparatus |
Also Published As
| Publication number | Publication date |
|---|---|
| EP0744157B1 (en) | 2003-07-23 |
| EP0744157A1 (en) | 1996-11-27 |
| CN1256917C (zh) | 2006-05-24 |
| AU5236396A (en) | 1996-12-05 |
| DE69629149T2 (de) | 2004-02-19 |
| AU704755B2 (en) | 1999-05-06 |
| KR100413919B1 (ko) | 2004-04-09 |
| JPH08308823A (ja) | 1996-11-26 |
| DE69629149D1 (de) | 2003-08-28 |
| KR960040317A (ko) | 1996-12-17 |
| US5671743A (en) | 1997-09-30 |
| CN1141156A (zh) | 1997-01-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP2927209B2 (ja) | ディジタルアンギオグラフィ装置 | |
| US7283614B2 (en) | X-ray diagnosis apparatus and method for creating image data | |
| JP4172753B2 (ja) | Ｘ線画像撮影装置 | |
| EP0088935A1 (en) | Matched filter for X-ray temporal subtraction | |
| JP2009066447A (ja) | 実時間デジタルｘ線撮像装置 | |
| JP3435192B2 (ja) | Ｘ線診断装置 | |
| JPS63290547A (ja) | テレビ断層撮影装置 | |
| JPH06285061A (ja) | Ｘ線診断方法及びその装置 | |
| EP0917855A1 (en) | X-Ray photographing apparatus and method capable of performing computerized tomography using C-Arm | |
| JP2001045374A (ja) | デジタルサブトラクション装置 | |
| JP4506119B2 (ja) | ディジタルサブトラクション装置およびディジタルアンギオグラフィ装置並びにディジタルサブトラクション像の表示方法 | |
| JP3489750B2 (ja) | Ｘ線ｃｔ装置 | |
| JP4626086B2 (ja) | ディジタルサブトラクション装置 | |
| JP2000116637A (ja) | 画像表示方法および装置並びに放射線透視撮影装置 | |
| JP2701721B2 (ja) | ディジタルアンギオグラフィ装置 | |
| JP4639526B2 (ja) | ディジタルサブトラクション装置 | |
| JP2820195B2 (ja) | デジタルｘ線撮影装置 | |
| JP2002335448A (ja) | ディジタルサブトラクション装置 | |
| JP4498525B2 (ja) | スキャン条件設定方法およびｘ線ｃｔ装置 | |
| JP3189385B2 (ja) | Ｘ線透視撮影装置 | |
| JPS6010600A (ja) | Ｘ線動態画像診断装置 | |
| JP3210358B2 (ja) | Ｘ線断層像撮影装置 | |
| JPH0551294B2 (ja) | ||
| JP3259664B2 (ja) | Ｘ線サブトラクション映像装置 | |
| Wilson et al. | Peripheral angiography digital image acquisition using rapid continuous sweeps of a moving gantry |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080514 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090514 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100514 Year of fee payment: 11 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100514 Year of fee payment: 11 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110514 Year of fee payment: 12 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110514 Year of fee payment: 12 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120514 Year of fee payment: 13 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130514 Year of fee payment: 14 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130514 Year of fee payment: 14 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140514 Year of fee payment: 15 |
|
| LAPS | Cancellation because of no payment of annual fees |