JP3163131B2 - X-ray equipment - Google Patents

X-ray equipment

Info

Publication number
JP3163131B2
JP3163131B2 JP26953191A JP26953191A JP3163131B2 JP 3163131 B2 JP3163131 B2 JP 3163131B2 JP 26953191 A JP26953191 A JP 26953191A JP 26953191 A JP26953191 A JP 26953191A JP 3163131 B2 JP3163131 B2 JP 3163131B2
Authority
JP
Japan
Prior art keywords
image
average density
block
ray
contrast
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
Application number
JP26953191A
Other languages
Japanese (ja)
Other versions
JPH0523326A (en
Inventor
道隆 本田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Toshiba Corp
Original Assignee
Toshiba Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Toshiba Corp filed Critical Toshiba Corp
Priority to JP26953191A priority Critical patent/JP3163131B2/en
Publication of JPH0523326A publication Critical patent/JPH0523326A/en
Application granted granted Critical
Publication of JP3163131B2 publication Critical patent/JP3163131B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【産業上の利用分野】本発明は、造影X線撮影において
造影像の判定機能を備えたX線撮影装置に関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an X-ray imaging apparatus having a function of determining a contrast image in contrast X-ray photography.

【0002】[0002]

【従来の技術】被検体の目的血管内にカテーテル,ガイ
ド,ワイヤー等を正確に挿入するため、X線透視撮影に
おいて血管内に造影剤を注入して撮影する造影X線撮影
が行われている。しかしながら造影剤は、血管内に注入
しても撮影領域からすぐに流れ出てしまうため造影剤が
注入された時のX線像(動画像)を得るにはその注入タ
イミングを検出する必要がある。
2. Description of the Related Art In order to accurately insert a catheter, a guide, a wire, and the like into a target blood vessel of a subject, contrast X-ray imaging is performed in which a contrast agent is injected into a blood vessel in X-ray fluoroscopy. . However, even if the contrast agent is injected into the blood vessel, it immediately flows out of the imaging region. Therefore, in order to obtain an X-ray image (moving image) when the contrast agent is injected, it is necessary to detect the injection timing.

【0003】そのタイミングを検出する技術としては、
次の2つのものが知られている。
[0003] As a technique for detecting the timing,
The following two are known.

【0004】第1の技術は、特開昭64-58243号公報に開
示されたもので、画面全体の平均濃度の変化より造影剤
注入タイミングを検出するものである。この方法は、造
影剤が撮影領域に混入すると画面全体の平均濃度が低下
して暗くなることを利用したものである。
A first technique is disclosed in Japanese Patent Application Laid-Open No. 64-58243, and detects a contrast agent injection timing from a change in the average density of the entire screen. This method utilizes the fact that the average density of the entire screen decreases and darkens when a contrast agent mixes into the imaging region.

【0005】第2の技術は、オペレータが造影剤注入と
同時に透視スイッチをOFFにし、そのOFFとした直
前の画像をホールドする「ラストイメージホールド」と
称せられているものが知られている。
As the second technique, there is known a technique called "last image hold" in which an operator turns off a fluoroscopic switch simultaneously with injection of a contrast agent and holds an image immediately before the turning off.

【0006】[0006]

【発明が解決しようとする課題】しかしながら、上記第
1の技術は、画面全体についての平均濃度をみているた
め、被検体の動き等があるとそれを造影剤注入と判断し
ていまい、注入タイミングの検出精度が悪いという問題
があった。
However, in the first technique, since the average density of the entire screen is viewed, if there is a movement of the subject, it is not determined that the contrast medium is injected. There is a problem that the detection accuracy is poor.

【0007】また上記第2の技術は、透視スイッチのO
FFを造影剤注入に同期させなければならず、操作が煩
雑という問題があった。
[0007] The above-mentioned second technique is based on the O switch of the perspective switch.
The FF must be synchronized with the injection of the contrast agent, and there is a problem that the operation is complicated.

【0008】そこで本発明は、上記事情に鑑みてなされ
たものであり、被検体の動き等があっても造影像を自動
認識可能なX線撮影装置を提供することを目的とするも
のである。
The present invention has been made in view of the above circumstances, and has as its object to provide an X-ray imaging apparatus capable of automatically recognizing a contrast image even when a subject moves. .

【0009】[0009]

【課題を解決するための手段】上記目的を達成するため
に本発明は、以下の構成とする。
To achieve the above object, the present invention has the following arrangement.

【0010】請求項1記載の発明は、被検体のX線撮影
中に前記被検体内に造影剤を注入し時経列的にディジタ
ルX線像を得るX線撮影装置において、前記ディジタル
X線像を複数のブロックに分割し各ブロック毎に平均濃
度を求める平均濃度算出手段と、この平均濃度算出手段
から得る各ブロック毎の平均濃度を用い、時経列的に近
接するブロック間の平均濃度差の時経列的な変化態様を
各位置のブロック毎に求める手段と、この求められる
均濃度差の時系列的変化態様に基づいて前記各ディジタ
ルX線像の造影状態を判定する判定手段とを有すること
を特徴とするものである。
According to a first aspect of the present invention, there is provided an X-ray imaging apparatus which obtains a digital X-ray image in time sequence by injecting a contrast agent into the object during X-ray imaging of the object. The image is divided into multiple blocks, and the average density of each block is
Mean density calculating means for determining the degree, and the mean density calculating means
Using the average density of each block obtained from
How the average density difference between adjacent blocks changes over time
Means for determining for each block of each position, that it has a determination means for pre-Symbol contrast status of each digital X-ray image based on the time-series changing mode of the flat <br/> average density difference is the calculated It is a feature.

【0011】請求項2記載の発明は、請求項1記載の発
明において、前記判定手段は、前記平均濃度差を符号化
した時経列的変化態様に基づくものである。
According to a second aspect of the present invention, in the first aspect of the present invention, the determination means is based on a time-series change mode in which the average density difference is encoded.

【0012】請求項3記載の発明は、被検体のX線撮影
中に前記被検体内に造影剤を注入し時経列的にディジタ
ルX線像を得るX線撮影装置において、前記ディジタル
X線像を複数のブロックに分割し各ブロック毎に平均濃
度を求める平均濃度算出手段と、この平均濃度算出手段
から得る各ブロック毎の平均濃度の出力態様に基づい
て、前記各ディジタルX線像の造影状態を判定する判定
手段とを有することを特徴とするものである。
According to a third aspect of the present invention, there is provided an X-ray imaging apparatus which obtains a digital X-ray image in time sequence by injecting a contrast agent into the object during X-ray imaging of the object. The image is divided into multiple blocks, and the average density of each block is
Mean density calculating means for determining the degree, and the mean density calculating means
And determining means for judging the contrast state of each digital X-ray image based on the output mode of the average density for each block obtained from the above.

【0013】請求項4記載の発明は、請求項3記載の発
明において、前記判定手段は、前記平均濃度を符号化し
た出力態様に基づくものである。
According to a fourth aspect of the present invention, in the third aspect of the present invention, the determination means is based on an output mode in which the average density is encoded.

【0014】[0014]

【作用】上記構成の装置の作用を説明する。The operation of the apparatus having the above configuration will be described.

【0015】請求項1記載の発明によれば、ブロック像
形成手段により複数のブロックに分割した各ブロック毎
に平均濃度を求めることにより、ブロック単位で造影判
定の精度が向上する。また造影剤混入と、被検体の動き
等とでは平均濃度差の時経列的変化態様が異なる。この
ため判定手段は、被検体の動き等があってもその変化態
様に基づいて造影像の判定を正確に行える。
According to the first aspect of the present invention, the average density is obtained for each block divided into a plurality of blocks by the block image forming means, so that the accuracy of the contrast determination for each block is improved. Further, the time-series change mode of the average density difference differs between the contrast agent mixing and the movement of the subject. For this reason, the determination unit can accurately determine the contrast image based on the change mode even when the subject moves or the like.

【0016】請求項2記載の発明によれば、平均濃度差
を符号化することにより、判定が容易となり、判定の迅
速化が図れる。
According to the second aspect of the present invention, by encoding the average density difference, the determination is facilitated and the determination can be speeded up.

【0017】請求項3記載の発明によれば、ブロック像
形成手段により複数のブロックに分割した各ブロック毎
に平均濃度を求めることにより、ブロック単位で造影判
定の精度が向上する。また造影剤混入と、被検体の動き
等とでは平均濃度の出力態様が異なる。このため判定手
段は、被検体の動き等があってもその出力態様に基づい
て造影像の判定を正確に行える。
According to the third aspect of the present invention, the accuracy of contrast determination is improved for each block by obtaining the average density for each block divided into a plurality of blocks by the block image forming means. The output mode of the average density differs between the contrast agent mixing and the movement of the subject. For this reason, the determination means can accurately determine the contrast image based on the output mode even when the subject moves or the like.

【0018】請求項4記載の発明によれば、平均濃度を
符号化することにより、判定が容易となり、判定の迅速
化が図れる。
According to the fourth aspect of the present invention, by encoding the average density, the determination is facilitated and the determination can be speeded up.

【0019】[0019]

【実施例】以下に本発明の実施例を図面を参照して詳述
する。
Embodiments of the present invention will be described below in detail with reference to the drawings.

【0020】図1は本発明の第1の実施例のX線撮影装
置1の概略構成図を示すものである。
FIG. 1 is a schematic block diagram of an X-ray imaging apparatus 1 according to a first embodiment of the present invention.

【0021】本装置1は、図示しないX線発生源と、こ
のX線発生源からのX線の被検体透過X線を検出して光
学像に変換するイメージ・インテンシファイア(以下
「I.I.」という。)2と、I.I.2の光学像をT
Vカメラ3に導く光学系4と、TVカメラ3から出力さ
れるTV映像信号をディジタルX線像に変換するA/D
変換器5と、原画像としてディジタルX線像を記憶する
原画像フレームメモリ6と、同じくこの原画像であるデ
ィジタルX線像を一時的に記憶する一対のバッファ7
と、A/D変換器5から直接取り込んだディジタルX線
像(Xi)と1フレーム遅れのバッファ7から取り込ん
だディジタルX線像(Xi-1 )との減算処理を行う減算
部8と、例えば原画像が512×512マトリックスの
ときに、この原画像を図2に示すように128×128
を1ブロックとして16ブロック(B00乃至B33)に分
割し、減算処理後の減算像Giを記憶する減算像フレー
ムメモリ(ブロック像形成手段)9と、その減算像Gi
の各ブロック毎に平均濃度を求めブロック像Siを形成
する平均値演算部(ブロック像形成手段)10と、スレ
ッシュホールド値T1 ,T2 を出力するスレッシュホー
ルド発生部11と、この出力されたスレッシュホールド
値T1 ,T2 とブロック像Siの各ブロックの平均濃度
とを比較して、符号化像Siを形成する比較部12と、
その符号化像Biを記憶する複数の符号化像メモリ13
と、各符号化像メモリ13に記憶されているデータに基
づいて原画像フレームメモリ6に記憶された原画像が造
影像(造影剤混入)か否かの見わけを行い造影像と判定
した場合に混入信号Smを出力する判定部14と、オペ
レーション信号Soと混入信号Smとが入力された場合
に原画像フレームメモリ6から送出されたディジタルX
線像を可視表示する表示系16とを有するものである。
The apparatus 1 includes an X-ray source (not shown) and an image intensifier (hereinafter referred to as “I. I.) 2) and I. I. The optical image of
An optical system 4 for guiding to a V camera 3 and an A / D for converting a TV video signal output from the TV camera 3 into a digital X-ray image
A converter 5, an original image frame memory 6 for storing a digital X-ray image as an original image, and a pair of buffers 7 for temporarily storing the digital X-ray image which is also the original image
A subtraction unit 8 for performing a subtraction process between the digital X-ray image (Xi) directly taken from the A / D converter 5 and the digital X-ray image (X i-1 ) taken from the buffer 7 delayed by one frame; For example, when the original image is a 512 × 512 matrix, the original image is 128 × 128 as shown in FIG.
Is divided into 16 blocks (B 00 to B 33 ) as one block, and a subtracted image frame memory (block image forming means) 9 for storing the subtracted image Gi after the subtraction processing, and the subtracted image Gi
An average value calculating unit (block image forming means) 10 for obtaining an average density for each block and forming a block image Si, a threshold generating unit 11 for outputting threshold values T 1 and T 2 , A comparison unit 12 that compares the threshold values T 1 and T 2 with the average density of each block of the block image Si to form an encoded image Si;
A plurality of coded image memories 13 for storing the coded images Bi
And whether the original image stored in the original image frame memory 6 is a contrast image (contrast agent is mixed) based on the data stored in each of the encoded image memories 13 and is determined to be a contrast image. And a decision unit 14 for outputting a mixed signal Sm to the digital X-ray generator 30 when the operation signal So and the mixed signal Sm are input.
A display system 16 for visually displaying a line image.

【0022】なお、前記原画像フレームメモリ6は、図
3に示すピークトレース部60に置換えてもよい。同図
に示すピークトレース部60は、比較器61と、セレク
タ62と、フレームメモリ63とを有するものである。
比較器61は、A/D5の出力とフレームメモリ63と
の間で対応する画素のデータを各々比較し、小さい方を
出力するものである。この比較は、1フレーム間で全画
素をスキャンすることにより行われる。セレクタ62
は、判定部14からの混入信号Smが出力されていない
場合(インアティブの場合)は、A/D5の出力をフレ
ームメモリ63にそのまま書き込み、混入信号Smが出
力されている場合(アクティブの場合)は、比較器61
の出力をフレームメモリ63に書き込むようにしてい
る。このような構成とすることにより、混入信号Smが
出力されている期間中は、フレームメモリ63の各画素
(ピクセル)には、最も小さい値が記録される。従っ
て、混入信号Smがオフとなる直前は、造影剤が消失し
ていてもフレームメモリ63には、最も濃く染まった状
態の原画像が残っていることになる。
The original image frame memory 6 may be replaced with a peak trace unit 60 shown in FIG. The peak trace unit 60 shown in the figure has a comparator 61, a selector 62, and a frame memory 63.
The comparator 61 compares the data of the corresponding pixel between the output of the A / D 5 and the frame memory 63, and outputs the smaller one. This comparison is performed by scanning all pixels in one frame. Selector 62
In the case where the mixing signal Sm from the determination unit 14 is not output (in the case of inactive), the output of the A / D 5 is written to the frame memory 63 as it is, and the mixing signal Sm is output (in the case of active). Is the comparator 61
Is written to the frame memory 63. With such a configuration, the smallest value is recorded in each pixel of the frame memory 63 during the period when the mixed signal Sm is being output. Therefore, immediately before the mixed signal Sm is turned off, even if the contrast agent has disappeared, the original image in the state of the darkest dye remains in the frame memory 63.

【0023】図4は前記減算部8の作用を示す図であ
る。この減算部8が減算対象とする画像Xi-1 は1フレ
ーム前のものでなくても被検体の動きが殆んど無視でき
る程であれば、数フレーム前のものでもよい。減算部8
は、図3に示すように画像Xiから1フレーム前の画像
i-1 を減算して減算後の減算像Giを得るものであ
る。尚同図中黒で示す領域は、造影剤が混入している部
分を示す。
FIG. 4 is a diagram showing the operation of the subtractor 8. The image X i-1 to be subtracted by the subtraction unit 8 may be several frames before as long as the movement of the subject can be almost ignored even if it is not one frame before. Subtraction unit 8
Is to obtain the subtraction image Gi after the subtraction by subtracting the image X i-1 of one frame before from an image Xi as shown in Figure 3. It should be noted that the region shown in black in the figure indicates the portion where the contrast agent is mixed.

【0024】前記平均値演算部10は、図4に示す減算
像Giについて図5に示すように、各ブロックb内の平
均濃度を求めブロック像Siを形成するものである。造
影剤が流れて消失した部分は、プラスの値となり、逆に
造影剤が流入した部分は、マイナスの値となり、その他
造影剤が流れていない部分は0の値となる。従って例え
ば図6に示すような値の分布となる。
As shown in FIG. 5, the average value calculating section 10 obtains the average density in each block b for the subtracted image Gi shown in FIG. 4 to form a block image Si. A portion where the contrast agent has flowed and disappeared has a positive value, a portion where the contrast agent has flowed in has a negative value, and a portion where no contrast agent has flown has a value of 0. Therefore, for example, a distribution of values as shown in FIG. 6 is obtained.

【0025】前記スレッシュホールド発生部11は、予
め造影剤が流入したときの減算値のスレッシュホールド
値T1 と、造影剤が流れて消失したときの減算値のスレ
ッシュホールド値T2 とをシュミレーション等により定
めておき、その値T1 ,T2 を比較部12に出力するも
のである。またこの値T1 ,T2 はT2 >T1 の関係に
ある。
The threshold generation unit 11 simulates in advance a threshold value T 1 of a subtraction value when the contrast medium flows in and a threshold value T 2 of a subtraction value when the contrast medium flows and disappears. The values T 1 and T 2 are output to the comparison unit 12. The values T 1 and T 2 have a relationship of T 2 > T 1 .

【0026】前記比較部12は、送出されたスレッシュ
ホールド値T1 ,T2 に基づいて、平均値演算部12か
ら出力されるブロック像Siを次に示すように「−],
「0」,「+」の符号に変換するものである。すなわ
ち、ブロックb内の平均値がスレッシュホールド値T1
よりも小さければ「−」,両値T1 ,T2 の間の値とな
っていれば「0」,スレッショホールド値T2 より大き
ければ「+」の符号とする。例えばスレッシュホールド
値T1 ,T2 をそれぞれ−5,5とすると図5に示した
ブロック像Siは、図6に示すように「−」,「0」,
「+」に符号化された符号化像Biとなる。比較部12
は上記操作を各フレーム毎に行い、符号化像Biをフレ
ーム毎に符号化像メモリ13に順に#1から#nまで記
憶するようにしている。例えば図7に示すように、原画
像X2 ,X1 との減算を行った結果最初の符号化像B
2-1 を符号化像メモリ(#1)13に格納し、次に原画
像X3 とX2 との減算を行った結果の符号化像B3-2
符号化像メモリ(#2)13に格納する。このようにし
て比較部12は、順次原画像X4 ,X5 についても同様
にして減算を行った結果の符号化像B4-3 ,B5-4 をそ
れぞれ符号化像メモリ(#3)13,メモリ(#4)1
3に格納するものである。また比較部12はメモリ(#
n)13までの記憶が終了したら、周期的に#1に戻
り、以前#1に記憶したデータを書き換えるようにして
もよい。
Based on the transmitted threshold values T 1 and T 2 , the comparing section 12 converts the block image Si output from the average value calculating section 12 into “−”,
This is converted into a sign of “0” or “+”. That is, the average value in the block b is equal to the threshold value T 1
If the value is smaller than the threshold value, the sign is “−”. If the value is between the two values T 1 and T 2 , the sign is “0”. If the value is larger than the threshold value T 2 , the sign is “+”. For example, assuming that the threshold values T 1 and T 2 are −5 and 5, respectively, the block image Si shown in FIG. 5 becomes “−”, “0”,
It becomes a coded image Bi coded to “+”. Comparison section 12
Performs the above operation for each frame, and stores the encoded image Bi in the encoded image memory 13 in order from # 1 to #n for each frame. For example, as shown in FIG. 7, as a result of performing subtraction with the original images X 2 and X 1 , the first encoded image B
2-1 encoded image memory (# 1) stored in 13, then the original image X 3 and X 2 and the subtracting the result of performing coding image B 3-2 coded image memory (# 2) 13 is stored. In this way, the comparison unit 12 sequentially stores the encoded images B 4-3 and B 5-4 resulting from the subtraction of the original images X 4 and X 5 in the same manner as the encoded image memory (# 3). 13, memory (# 4) 1
3 is stored. The comparison unit 12 has a memory (#
n) After the storage up to 13, the process may return to # 1 periodically to rewrite the data previously stored in # 1.

【0027】前記判定部14は、図8に示すように、第
1のテーブル141と、後述するコード番号の数に対応
して設けられた複数の計数器142と、アドレス発生器
143と、第2のテーブル144とを有している。
As shown in FIG. 8, the judging section 14 comprises a first table 141, a plurality of counters 142 provided corresponding to the number of code numbers to be described later, an address generator 143, 2 tables 144.

【0028】アドレス発生器143は、各符号化像メモ
リ13に同一のアドレスを同時にn個発生した後、次の
アドレスを各符号化像メモリ13に同時にn個発生し、
これを最終アドレスまで繰り返すものである。
The address generator 143 generates the same n addresses simultaneously in each coded image memory 13, and then generates the next n addresses simultaneously in each coded image memory 13.
This is repeated until the last address.

【0029】第1のテーブル141は、各符号化像メモ
リ13から各アドレスごとに取り込んだデータに基づい
て、コード番号(「コード0」,「コード1」等)を判
定し、そのコード番号発生するものである。例えば、最
初のアドレスについてのコード番号が「コード0」であ
るときは、「コード0」用の計数器142に「1」を入
力し、他の計数器142には「0」を入力する。
The first table 141 determines a code number (“code 0”, “code 1”, etc.) based on data fetched for each address from each coded image memory 13 and generates the code number. Is what you do. For example, when the code number for the first address is “code 0”, “1” is input to the counter 142 for “code 0”, and “0” is input to the other counters 142.

【0030】この第1のテーブル141によるコード番
号の判定は、例えば、あるアドレスについて図9に示す
ようなデータが得られたとする。第1例に示すように全
て「0」の場合は「コード0」、第2例に示すように
「0」の間に「−」が続く場合(位置ずれ)は「コード
1」、第3例に示すように「0」の間に「+」が続く場
合(位置ずれ)は「コード2」、第4例に示すように
「0」,「−」,「+」,「0」の順で現れた場合(造
影剤混入)は「コード3」、第4例に示すように
「0」,「−」,「+」がランダムで現れた場合は「コ
ード4」というように予め判定基準を決めておき、その
基準に従ってコード番号を出力する。
The determination of the code number by the first table 141 is based on the assumption that data as shown in FIG. 9 is obtained for a certain address. As shown in the first example, if all are “0”, “code 0” is used. As shown in the second example, if “−” is continued between “0” (position shift), “code 1” is used. As shown in the example, when "+" continues between "0" (position shift), "code 2" is used, and as shown in the fourth example, "0", "-", "+", "0" If they appear in order (contrast agent mixture), they are determined in advance as "code 3", and as shown in the fourth example, if "0", "-", and "+" appear randomly, they are determined as "code 4". A standard is determined, and a code number is output according to the standard.

【0031】すなわち、造影剤混入のパターンは、各符
号化像Bi中の1つのブロックb1(前記図7参照)に
注目してnフレーム間の変化をみると、造影剤が混入し
た場合は必ず「0」から「−」に変化し、更に「+」そ
して「0」へと変化する(コード3)。一方、被検体の
原画像が例えば図10に示すようにX10,X11,X12
と左から右へ移動した場合、その他カテーテル等の体内
挿入物等が体内で移動した場合は、符号化像B11-10
12-11 は、急に「0」から「+」(ブロックb2 ),
「−」から「0」(ブロックb3 )又は「−」から
「+」(ブロックb4 )等のように造影剤流入とは異な
る変化態様を示す(コード1又はコード2)。
That is, the pattern of the contrast agent mixing shows that the change between n frames is focused on one block b1 (see FIG. 7) in each encoded image Bi. It changes from "0" to "-", and further changes to "+" and "0" (code 3). On the other hand, when the original image of the subject moves from left to right to X 10 , X 11 , and X 12 , as shown in FIG. Image B 11-10 ,
B12-11 suddenly changes from "0" to "+" (block b2),
A change mode different from the contrast agent inflow, such as "-" to "0" (block b3) or "-" to "+" (block b4), is shown (code 1 or code 2).

【0032】計数器142は、第1のテーブル141か
ら「1」が入力される度にカウントアップし、最後に各
フレームごとにその数を合計して第2のテーブル144
に出力するものである。
The counter 142 counts up each time "1" is input from the first table 141, and finally sums up the number for each frame to produce a second table 144.
Is output to

【0033】第2のテーブル144は、各計数器142
からの出力値から原画像が造影像か否かの判定し、造影
像と判定した場合は、混入信号Smを出力するものであ
る。例えば、各計数器142からの出力値により「コー
ド3」が少なくとも所定数以上あり、かつ、「コード
0,1,4」の数が所定数以下ならば原画像が造影像
(造影剤混入)と判定する。また、「コード1」又は
「コード2」が所定数以上であるときは、位置ずれと判
定し、混入信号Smは出力しない。このように、判定部
14は、符号化像Biの時系列的変化態様に基づいて被
検体の移動等とは区別して造影像か否かの判定を正確に
行えるものである。
The second table 144 stores each counter 142
It is determined whether or not the original image is a contrast image based on the output value from. If it is determined that the original image is a contrast image, a mixed signal Sm is output. For example, if the number of “codes 3” is at least a predetermined number or more and the number of “codes 0, 1, 4” is equal to or less than a predetermined number according to the output value from each counter 142, the original image is a contrast image (contrast agent mixed). Is determined. If the number of “code 1” or “code 2” is equal to or more than a predetermined number, it is determined that the position is misaligned, and the mixed signal Sm is not output. As described above, the determination unit 14 can accurately determine whether or not the image is a contrast image based on the time-series change mode of the encoded image Bi, in distinction from the movement of the subject.

【0034】次にこのように構成された上記実施例装置
1の作用を図11をも参照して説明する。
Next, the operation of the above-described embodiment apparatus 1 will be described with reference to FIG.

【0035】オペレータが被検体のX線透視撮影を行い
ながら被検体内に少量の造影剤を注入する。X線撮影に
より被検体を透過したX線は、I.I.2により光学像
に変換され光学系4を経由してTVカメラ3によりTV
映像信号に変換され、A/D変換器5からディジタルX
線像として出力される。このディジタルX線像には、造
影剤混入前のもの,混入中,混入して流れ去った後のも
のが含まれる。このディジタルX線像は全て原画像フレ
ームメモリ6に記憶され、一方符号化像Bi生成のため
に減算部8,一対のバッファ7にも出力される。
An operator injects a small amount of a contrast agent into the subject while performing X-ray fluoroscopy of the subject. X-rays transmitted through a subject by X-ray photography are I.D. I. 2 is converted to an optical image by the TV camera 3 via the optical system 4
It is converted to a video signal, and the digital X
It is output as a line image. The digital X-ray images include those before mixing the contrast agent, those during mixing, and those after mixing and flowing away. The digital X-ray images are all stored in the original image frame memory 6, and are also output to a subtraction unit 8 and a pair of buffers 7 for generating an encoded image Bi.

【0036】次に減算部8は、A/D変換器5から送出
された画像Xi,Xi-1 に基づいて減算像Giを形成す
る。また減算部8は、この減算像Giを減算像フレーム
メモリ9に一時記憶する。次に平均演算部10は、減算
像Giに基づいてブロック像Siを形成する。比較部1
2は、スレッシュホールド発生部11からのスレッシュ
ホールド値T1 ,T2 に基づいてブロック像Siを符号
化した符号化像Biを形成し、順次符号化像メモリ13
に格納する。判定部14は各符号化像メモリ13に格納
された符号化像Biの時経列的変化態様に基づいて、原
画像フレームメモリ6に取り込まれた原画像を造影像と
判定したときは、混入信号Smを原画像フレームメモリ
6と表示部16に送出する。
Next, the subtraction section 8 forms a subtraction image Gi based on the images Xi and Xi -1 sent from the A / D converter 5. Further, the subtraction unit 8 temporarily stores the subtraction image Gi in the subtraction image frame memory 9. Next, the averaging unit 10 forms a block image Si based on the subtraction image Gi. Comparison section 1
2 forms an encoded image Bi obtained by encoding the block image Si based on the threshold values T 1 and T 2 from the threshold generator 11 and sequentially encodes the encoded image memory 13
To be stored. When the determination unit 14 determines that the original image captured in the original image frame memory 6 is a contrast image based on the time-series change mode of the encoded image Bi stored in each encoded image memory 13, The signal Sm is sent to the original image frame memory 6 and the display unit 16.

【0037】次に原画像フレームメモリ6は、判定部1
4から送出された混入信号Smを表示系16に出力す
る。オペレータの造影剤注入と同時に行うスイッチ操作
によりオペレーション信号Soが表示系16に送出さ
れ、更に混入信号Smが表示系16に送出されると、表
示系16は混入信号Smとオペレーション信号Soの両
信号Sm,Soが揃った所で原画像フレームメモリ6か
ら取込んだ原画像をすなわち造影像を表示すると共に判
定部14から送出された混入信号Smを表示する。オペ
レータは、オペレーション信号Soを発信直後に混入信
号Smが表示系16に表示されたことにより判定部14
が正確に判定していることを確認する。そして表示系1
6に表示された画像が診断に供される。
Next, the original image frame memory 6 stores the judgment unit 1
4 to the display system 16. When an operation signal So is transmitted to the display system 16 by a switch operation performed simultaneously with the injection of the contrast agent by the operator, and the mixed signal Sm is further transmitted to the display system 16, the display system 16 outputs both the mixed signal Sm and the operation signal So. When Sm and So are complete, the original image fetched from the original image frame memory 6, that is, a contrast image is displayed, and the mixed signal Sm sent from the determination unit 14 is displayed. The operator determines that the mixed signal Sm is displayed on the display system 16 immediately after transmitting the operation signal So.
Make sure that is correctly determined. And display system 1
The image displayed in 6 is used for diagnosis.

【0038】このような実施例装置1によれば、オペレ
ータは被検体の動き等があっても造影像を自動的にしか
も確実に得ることができる。
According to the apparatus 1 of the embodiment, the operator can automatically and reliably obtain a contrast image even when the subject moves.

【0039】またオペレーション信号を入力するように
しているので、仮に、造影剤を注入していない状態にも
かかわらず、被検体の動きなどで誤って判定部14が造
影剤混入と判断しても、表示系16に表示されないた
め、通常の透視状態を妨げることがない。
Further, since the operation signal is input, even if the judgment unit 14 erroneously judges that the contrast agent is mixed due to the movement of the subject even though the contrast agent is not injected. Since it is not displayed on the display system 16, it does not interfere with the normal see-through state.

【0040】図12は本発明の第2の実施例のX線撮影
装置20の概略構成図を示すものである。第1の実施例
装置1との主な相違点は、スレッシュホールド発生部2
1と、比較部22と、符号化像メモリ23と、判定部2
4とにある。なお、第1の実施例装置1と機能が同一の
ものには、同一の符号を付してその詳細な説明は省略す
る。なお、第1の実施例装置1と同様に本第2の実施例
装置20においても前記原画像フレームメモリ6を、図
3に示すピークトレース部60に置換えてもよい。
FIG. 12 is a schematic configuration diagram of an X-ray imaging apparatus 20 according to a second embodiment of the present invention. The main difference from the first embodiment 1 is that the threshold generator 2
1, a comparing unit 22, an encoded image memory 23, and a determining unit 2
There are four. The components having the same functions as those of the device 1 of the first embodiment are denoted by the same reference numerals, and the detailed description thereof will be omitted. It should be noted that the original image frame memory 6 may be replaced by the peak trace unit 60 shown in FIG. 3 in the second embodiment device 20 as in the first embodiment device 1.

【0041】前記スレッシュホールド発生部21は、第
1の実施例と比較してより多く(本例では4つ)のスレ
ッシュホールド値T21,T22,T23,T24を出力するも
のである。多数のスレッシュホールド値T21乃至T24を
設けることにより、判定部24の判定精度が高められ
る。
The threshold generator 21 outputs more (four in this embodiment) threshold values T21, T22, T23 and T24 as compared with the first embodiment. By providing a large number of threshold values T21 to T24, the determination accuracy of the determination unit 24 can be increased.

【0042】前記比較部22は、送出されたスレッシュ
ホールド値T21乃至T24に基づいて、平均値演算部12
から出力されるブロック像Siを次に示すように「−
−],「−」,「0」,「+」,「++」の符号に変換
するものである。すなわち、ブロックb内の平均値がス
レッシュホールド値T21より小さければ「−−」、スレ
ッシュホールド値T21と値T22との間となっていれば
「−」、スレッシュホールド値T22と値T23との間とな
っていれば「0」、スレッシュホールド値T23と値T24
との間となっていれば「+」、スレッシュホールド値T
24より大きければ「++」の符号として符号化像を作成
し、符号化像メモリ23に記憶する。
The comparing section 22 calculates the average value calculating section 12 based on the transmitted threshold values T21 to T24.
The block image Si output from is expressed as "-
−], “−”, “0”, “+”, and “++”. That is, “−−” when the average value in the block b is smaller than the threshold value T21, “−” when the average value is between the threshold value T21 and the value T22, and “−−” between the threshold value T22 and the value T23. "0" when the threshold value is reached, the threshold value T23 and the value T24
If it is between “+”, the threshold value T
If it is larger than 24, an encoded image is created as a code of “++” and stored in the encoded image memory 23.

【0043】前記符号化像メモリ23は、第1の実施例
装置1は複数有していたが、本第2の実施例では1つと
している。
The coded image memory 23 has a plurality of devices 1 in the first embodiment, but one in the second embodiment.

【0044】前記判定部24は、図13に示すように、
分配器241と、符号の数に対応して設けられた複数
(本例では5つ)の計数器242と、アドレス発生器2
43と、テーブル244とを有している。
As shown in FIG. 13, the determination section 24
A distributor 241, a plurality (five in this example) of counters 242 provided corresponding to the number of codes, and an address generator 2
43 and a table 244.

【0045】アドレス発生器243は、符号化像メモリ
23に最初のアドレスから最終アドレスまで順次アドレ
スを発生するものである。
The address generator 243 generates addresses sequentially from the first address to the last address in the coded image memory 23.

【0046】分配器241は、符号化像メモリ23に記
憶された符号化像のある画素の符号が「−−」であると
きは、符号「−−」用の計数器242に「1」を入力
し、他の計数器242には「0」を入力するという動作
を他の符号についても同様に行うものである。
When the sign of a certain pixel of the coded image stored in the coded image memory 23 is “−−”, the distributor 241 sets “1” to the counter 242 for the code “−−”. The operation of inputting "0" to the other counter 242 is similarly performed for other codes.

【0047】計数器242は、分配器241から「1」
が入力される度にカウントアップし、最後に各フレーム
ごとにその数を合計してテーブル244に出力するもの
である。
The counter 242 outputs “1” from the distributor 241.
Is counted up each time is input, and finally the number is summed up for each frame and output to the table 244.

【0048】テーブル244は、各計数器242からの
出力値から原画像が造影像か否かを判定し、造影像と判
定した場合は、混入信号Smを出力するものである。例
えば、造影剤が混入した場合は、符号「−」が多くなり
「+」が少なくなる。また、符号「−」が少なく「+」
が多いときは造影剤が流出していくフェーズと判定され
る。この場合は、造影像と判定し、混入信号Smを出力
する。また、位置ずれの場合は「−」と「+」の数がほ
ぼ同数となり、この場合は混入信号Smを出力しない。
他の例としては、「−」が極端に多くなる場合や、「+
+」の数が少なくとも1つあり、「−」,「−−」の数
が殆どない場合などがある。前者は、X線出力変化など
により画像全体のレベル変化が生じたと考えられる。
The table 244 determines whether or not the original image is a contrast image from the output values from the respective counters 242, and outputs a mixed signal Sm when the image is determined to be a contrast image. For example, when a contrast agent is mixed, the sign “−” increases and the sign “+” decreases. Also, the sign “−” is small and “+”
When there is much, it is determined that the phase in which the contrast agent flows out. In this case, it is determined that the image is a contrast image, and the mixed signal Sm is output. Further, in the case of a positional shift, the numbers of “−” and “+” are substantially the same, and in this case, the mixed signal Sm is not output.
Other examples include extremely large numbers of "-" or "+"
There may be a case where there is at least one "+" and there are few "-" and "-". In the former case, it is considered that the level of the entire image has changed due to a change in the X-ray output or the like.

【0049】このように構成された上記第2の実施例装
置20によれば、1枚の符号化像により造影剤混入の有
無を判定できるので、第1の実施例装置1と比較して、
造影像の判定を迅速に行える。
According to the apparatus 20 of the second embodiment configured as described above, it is possible to determine the presence or absence of the contrast agent from one coded image, and therefore, compared to the apparatus 1 of the first embodiment,
The contrast image can be determined quickly.

【0050】なお、本発明は、上記各実施例に限定され
るものではなく、その要旨を変更しない範囲内で種々に
変形実施可能である。例えば、心臓の動き等の周期的に
変化する部位について造影X線撮影を行う場合は、R波
等の心拍同期信号の如き位相信号を基に、心位相の同じ
画像に対して減算処理を行い、同相の画像に対して本実
施例と同様に符号化像を形成して造影像を見つけ出すよ
うにしてもよい。
It should be noted that the present invention is not limited to the above-described embodiments, and various modifications can be made without departing from the scope of the invention. For example, when performing a contrast-enhanced X-ray imaging of a periodically changing portion such as a heart motion, a subtraction process is performed on an image having the same cardiac phase based on a phase signal such as a heartbeat synchronization signal such as an R wave. Alternatively, an encoded image may be formed on an in-phase image in the same manner as in the present embodiment to find a contrast image.

【0051】[0051]

【発明の効果】以上詳述した本発明によれば、以下の効
果を奏する。
According to the present invention described in detail above, the following effects can be obtained.

【0052】請求項1記載の発明によれば、ディジタル
X線像をブロック化し、各ブロックの平均濃度の時経列
的変化態様に基づいて各ディジタルX線像の造影判定を
するようにしているので、被検体の動き等があっても造
影像を自動認識可能なX線撮影装置を提供することがで
きる。
According to the first aspect of the present invention, the digital X-ray image is divided into blocks, and the contrast of each digital X-ray image is determined based on the time-series change mode of the average density of each block. Therefore, it is possible to provide an X-ray imaging apparatus capable of automatically recognizing a contrast image even when a subject moves or the like.

【0053】請求項2記載の発明によれば、請求項1記
載の効果に加え、平均濃度差を符号化しているので、判
定が容易となり、判定の迅速化が図れる。
According to the second aspect of the present invention, in addition to the effect of the first aspect, since the average density difference is coded, the determination is facilitated and the determination can be speeded up.

【0054】請求項3記載の発明によれば、ディジタル
X線像をブロック化し、各ブロックの平均濃度の出力態
様に基づいて各ディジタルX線像の造影判定をするよう
にしているので、被検体の動き等があっても造影像を自
動認識可能なX線撮影装置を提供することができる。
According to the third aspect of the present invention, the digital X-ray image is divided into blocks, and the contrast of each digital X-ray image is determined based on the output mode of the average density of each block. It is possible to provide an X-ray imaging apparatus capable of automatically recognizing a contrast image even when there is a movement of the image.

【0055】請求項4記載の発明によれば、請求項3記
載の効果に加え、平均濃度差を符号化しているので、判
定が容易となり、判定の迅速化が図れる。
According to the fourth aspect of the present invention, in addition to the effect of the third aspect, since the average density difference is encoded, the determination is facilitated and the determination can be speeded up.

【図面の簡単な説明】[Brief description of the drawings]

【図1】本発明の第1の実施例のX線撮影装置の概略構
成図である。
FIG. 1 is a schematic configuration diagram of an X-ray imaging apparatus according to a first embodiment of the present invention.

【図2】図1に示す装置の減算部の作用を示す図であ
る。
FIG. 2 is a diagram showing an operation of a subtraction unit of the device shown in FIG.

【図3】図1に示す原画像フレームメモリの他の例を示
す図である。
FIG. 3 is a diagram showing another example of the original image frame memory shown in FIG. 1;

【図4】図1に示す装置の減算部の作用を示す図であ
る。
FIG. 4 is a diagram showing an operation of a subtraction unit of the device shown in FIG.

【図5】ブロック像の一例を示す図である。FIG. 5 is a diagram illustrating an example of a block image.

【図6】符号化像の一例を示す図である。FIG. 6 is a diagram illustrating an example of an encoded image.

【図7】造影剤が混入した場合の経時的変化態様を示す
原画像と符号化像の一例を示す図である。
FIG. 7 is a diagram illustrating an example of an original image and a coded image showing a temporal change mode when a contrast agent is mixed.

【図8】図1に示す判定部の構成を示す図である。FIG. 8 is a diagram illustrating a configuration of a determination unit illustrated in FIG. 1;

【図9】図8に示す判定部の作用を示す図である。FIG. 9 is a diagram illustrating an operation of a determination unit illustrated in FIG. 8;

【図10】被検体が移動した場合の経時的変化態様を示
す原画像と符号化像の一例を示す図である。
FIG. 10 is a diagram showing an example of an original image and an encoded image showing a temporal change mode when a subject moves.

【図11】図1に示す装置の作用を示す図である。FIG. 11 is a view showing the operation of the device shown in FIG. 1;

【図12】本発明の第2の実施例のX線撮影装置の概略
構成図である。
FIG. 12 is a schematic configuration diagram of an X-ray imaging apparatus according to a second embodiment of the present invention.

【図13】図12に示す判定部の構成を示す図である。13 is a diagram illustrating a configuration of a determination unit illustrated in FIG.

【符号の説明】[Explanation of symbols]

1 X線撮影装置 9 減算像フレームメモリ(ブロック像形成手段) 10 平均値演算部(ブロック像形成手段) 14 判定部 Bi,B2-1 ,B3-2 ,B4-3 ,B5-4 ,B11-10 ,B
12-11 符号化像 Si ブロック像 X1 乃至X5 ,X10乃至X12 ディジタルX線像
DESCRIPTION OF SYMBOLS 1 X-ray imaging apparatus 9 Subtraction image frame memory (block image forming means) 10 Average value calculation part (block image forming means) 14 Judgment part Bi, B2-1 , B3-2 , B4-3 , B5- 4 , B11-10 , B
12-11 coded image Si block image X 1 to X 5, X 10 to X 12 digital X-ray image

Claims (4)

(57)【特許請求の範囲】(57) [Claims] 【請求項1】 被検体のX線撮影中に前記被検体内に造
影剤を注入し時経列的にディジタルX線像を得るX線撮
影装置において、前記ディジタルX線像を複数のブロッ
クに分割し各ブロック毎に平均濃度を求める平均濃度算
出手段と、この平均濃度算出手段から得る各ブロック毎
の平均濃度を用い、時経列的に近接するブロック間の平
均濃度差の時経列的な変化態様を各位置のブロック毎に
求める手段と、この求められる平均濃度差の時系列的変
化態様に基づいて前記各ディジタルX線像の造影状態
判定する判定手段とを有することを特徴とするX線撮影
装置。
1. An X-ray imaging apparatus which obtains a digital X-ray image in time sequence by injecting a contrast agent into the object during X-ray imaging of the object, wherein the digital X-ray image is divided into a plurality of blocks. Average density calculation to calculate the average density for each block by dividing
Output means and each block obtained from the average density calculation means.
Using the average density of
Time-sequential change of the average density difference for each block at each position
Means for determining, X-rays imaging apparatus characterized by having a determining means for pre-Symbol contrast status of each digital X-ray image based on the time-series changing mode of the average density difference is thus determined.
【請求項2】 前記判定手段は、前記平均濃度差を符号
化した時経列的変化態様に基づくものとする請求項1記
載のX線撮影装置。
2. The X-ray imaging apparatus according to claim 1, wherein the determination unit is based on a time-series change mode in which the average density difference is encoded.
【請求項3】 被検体のX線撮影中に前記被検体内に造
影剤を注入し時経列的にディジタルX線像を得るX線撮
影装置において、前記ディジタルX線像を複数のブロッ
クに分割し各ブロック毎に平均濃度を求める平均濃度算
出手段と、この平均濃度算出手段から得る各ブロック毎
の平均濃度の出力態様に基づいて、前記各ディジタルX
線像の造影状態を判定する判定手段とを有することを特
徴とするX線撮影装置。
3. An X-ray imaging apparatus which obtains a digital X-ray image in time sequence by injecting a contrast agent into the object during X-ray imaging of the object, wherein the digital X-ray image is divided into a plurality of blocks. Average density calculation to calculate the average density for each block by dividing
Output means and each block obtained from the average density calculation means.
Digital X based on the average density output mode
An X-ray imaging apparatus comprising: a determination unit configured to determine a contrast state of a line image.
【請求項4】 前記判定手段は、前記平均濃度を符号化
した出力態様に基づくものとする請求項3記載のX線撮
影装置。
4. The X-ray imaging apparatus according to claim 3, wherein the determination unit is based on an output mode in which the average density is encoded.
JP26953191A 1990-10-18 1991-10-17 X-ray equipment Expired - Fee Related JP3163131B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP26953191A JP3163131B2 (en) 1990-10-18 1991-10-17 X-ray equipment

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP27967490 1990-10-18
JP2-279674 1990-10-18
JP26953191A JP3163131B2 (en) 1990-10-18 1991-10-17 X-ray equipment

Publications (2)

Publication Number Publication Date
JPH0523326A JPH0523326A (en) 1993-02-02
JP3163131B2 true JP3163131B2 (en) 2001-05-08

Family

ID=26548812

Family Applications (1)

Application Number Title Priority Date Filing Date
JP26953191A Expired - Fee Related JP3163131B2 (en) 1990-10-18 1991-10-17 X-ray equipment

Country Status (1)

Country Link
JP (1) JP3163131B2 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1075083C (en) * 1995-02-03 2001-11-21 出光石油化学株式会社 Process for producing polycarbonates
JP7181006B2 (en) 2018-06-01 2022-11-30 出光興産株式会社 Polycarbonate resin, method for producing the same, and polycarbonate resin composition

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5439192B2 (en) * 2007-02-27 2014-03-12 コーニンクレッカ フィリップス エヌ ヴェ Method and apparatus for recording vasculature during medical intervention

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1075083C (en) * 1995-02-03 2001-11-21 出光石油化学株式会社 Process for producing polycarbonates
JP7181006B2 (en) 2018-06-01 2022-11-30 出光興産株式会社 Polycarbonate resin, method for producing the same, and polycarbonate resin composition

Also Published As

Publication number Publication date
JPH0523326A (en) 1993-02-02

Similar Documents

Publication Publication Date Title
US5565998A (en) Identifying film frames in a video sequence
US7254264B2 (en) Method and device for generating 3D images
US6707851B1 (en) Method for objects segmentation in video sequences by object tracking and user assistance
US20190335154A1 (en) Image processing apparatus, image processing method, and storage medium storing program
JPS60229594A (en) Method and device for motion interpolation of motion picture signal
US7499102B2 (en) Image processing apparatus using judder-map and method thereof
US7197075B2 (en) Method and system for video sequence real-time motion compensated temporal upsampling
JP3193833B2 (en) Motion vector processor
CN114881921B (en) Anti-occlusion imaging method and device based on event and video fusion
US7489728B2 (en) Apparatus and method for coding moving image
JPS6122513B2 (en)
JP3563687B2 (en) Method and apparatus for detecting repeated frames of video signal
US8213673B2 (en) Watermarking of motion pictures
US7356439B2 (en) Motion detection apparatus and method
JP2004056566A (en) X-ray picture diagnostic device
JP3163131B2 (en) X-ray equipment
JP2006109488A (en) Video processing device capable of selecting field and method thereof
KR101129591B1 (en) Image processing device, image processing method, and recording medium
JP2005517251A (en) Method and unit for estimating a motion vector of a pixel group
US8379146B2 (en) Deinterlacing method and apparatus for digital motion picture
EP0481488B1 (en) Method for discriminating x-ray contrast image from non-contrast image and x-ray imaging system including contrast image discriminating function
CN115835035A (en) Image frame interpolation method, device and equipment and computer readable storage medium
US9055176B1 (en) Method and apparatus for performing multi-cadence detection in a motion adaptive video deinterlacer
JP2802331B2 (en) Image forming device
JPS6158639A (en) Digital subtraction system

Legal Events

Date Code Title Description
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080223

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090223

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100223

Year of fee payment: 9

LAPS Cancellation because of no payment of annual fees