JP3223202B2 - データ補間処理装置 - Google Patents
データ補間処理装置Info
- Publication number
- JP3223202B2 JP3223202B2 JP27628892A JP27628892A JP3223202B2 JP 3223202 B2 JP3223202 B2 JP 3223202B2 JP 27628892 A JP27628892 A JP 27628892A JP 27628892 A JP27628892 A JP 27628892A JP 3223202 B2 JP3223202 B2 JP 3223202B2
- Authority
- JP
- Japan
- Prior art keywords
- data
- interpolation
- scanning
- sampling
- cubic curve
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/17—Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52053—Display arrangements
- G01S7/52057—Cathode ray tube displays
- G01S7/52073—Production of cursor lines, markers or indicia by electronic means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
- G01S7/52025—Details of receivers for pulse systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
- G01S7/52034—Data rate converters
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Image Processing (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Image Analysis (AREA)
Description
【０００１】
【産業上の利用分野】本発明は超音波診断装置に内蔵さ
れるデータ補間処理装置に関する。
れるデータ補間処理装置に関する。
【０００２】
【従来の技術】被検体、特に人体の体表に超音波トラン
スデューサを当てて超音波ビームを発信し、人体内を超
音波ビームで走査しながら人体内の組織で反射されて戻
ってきた超音波を超音波トランスデューサで受信し、こ
れにより得られた受信信号に基づく人体内の断層像を表
示することにより人体内の内臓等の疾患の診断を容易な
らしめる超音波診断装置が従来より用いられている。
スデューサを当てて超音波ビームを発信し、人体内を超
音波ビームで走査しながら人体内の組織で反射されて戻
ってきた超音波を超音波トランスデューサで受信し、こ
れにより得られた受信信号に基づく人体内の断層像を表
示することにより人体内の内臓等の疾患の診断を容易な
らしめる超音波診断装置が従来より用いられている。
【０００３】超音波ビームで被検体内を走査する方式と
して、例えば、多数の超音波ビーム（走査線）が扇状に
広がるセクタ走査方式、走査線が互いに並行に延びるリ
ニア走査方式、コンベックス型に超音波トランスデュー
サが並んだプローブを用いるコンベックス走査方式等が
知られている。上記のようにして得られた受信信号はサ
ンプリングされＡ／Ｄ変換されることによりディジタル
のサンプリングデータに変換されるが、このサンプリン
グデータが表わす各画素の位置は、通常、断層像が表示
される例えばＣＲＴディスプレイ装置の表示画面上の各
画素の位置とは対応していないため、データ補間処理を
行ない、表示に適した座標位置の表示用データを生成す
る必要がある。
して、例えば、多数の超音波ビーム（走査線）が扇状に
広がるセクタ走査方式、走査線が互いに並行に延びるリ
ニア走査方式、コンベックス型に超音波トランスデュー
サが並んだプローブを用いるコンベックス走査方式等が
知られている。上記のようにして得られた受信信号はサ
ンプリングされＡ／Ｄ変換されることによりディジタル
のサンプリングデータに変換されるが、このサンプリン
グデータが表わす各画素の位置は、通常、断層像が表示
される例えばＣＲＴディスプレイ装置の表示画面上の各
画素の位置とは対応していないため、データ補間処理を
行ない、表示に適した座標位置の表示用データを生成す
る必要がある。
【０００４】図７は従来のデータ補間処理演算の説明図
である。図の左上から右下に延びる直線はそれぞれ…，
ｎ−１，ｎ，ｎ＋１，ｎ＋２，…番目の走査線を表わ
し、図の左右に延びる直線は被検体内の同一の深さに並
ぶサンプリングデータの並びを表わす、走査方向に延び
る…，ｍ，ｍ＋１，…番目の直線である。
である。図の左上から右下に延びる直線はそれぞれ…，
ｎ−１，ｎ，ｎ＋１，ｎ＋２，…番目の走査線を表わ
し、図の左右に延びる直線は被検体内の同一の深さに並
ぶサンプリングデータの並びを表わす、走査方向に延び
る…，ｍ，ｍ＋１，…番目の直線である。
【０００５】各走査線…，ｎ−１，ｎ，ｎ＋１，ｎ＋
２，…と各直線…，ｍ，ｍ＋１，…との各交点に各サン
プリングデータが対応しており、ここでは、点Ｔの表示
用データを求めるにあたり、その周囲の４点Ｓ＿ｎ＿
ｍ，Ｓ＿ｎ＿ｍ＋１，Ｓ＿ｎ＋１＿ｍ，Ｓ＿ｎ＋１＿ｍ
＋１のサンプリングデータＳ（ｎ，ｍ），Ｓ（ｎ，ｍ＋
１），Ｓ（ｎ＋１，ｍ），Ｓ（ｎ＋１，ｍ＋１）が用い
られる。ここでは、図に示すように、直線ｍと直線ｍ＋
１との間の距離をΔＲ，直線ｍと点Ｔとの間の距離をＲ
err ，走査線ｎと走査線ｎ＋１との間の距離をΔθ、走
査線ｎと点Ｔとの間の距離をθerr とする。
２，…と各直線…，ｍ，ｍ＋１，…との各交点に各サン
プリングデータが対応しており、ここでは、点Ｔの表示
用データを求めるにあたり、その周囲の４点Ｓ＿ｎ＿
ｍ，Ｓ＿ｎ＿ｍ＋１，Ｓ＿ｎ＋１＿ｍ，Ｓ＿ｎ＋１＿ｍ
＋１のサンプリングデータＳ（ｎ，ｍ），Ｓ（ｎ，ｍ＋
１），Ｓ（ｎ＋１，ｍ），Ｓ（ｎ＋１，ｍ＋１）が用い
られる。ここでは、図に示すように、直線ｍと直線ｍ＋
１との間の距離をΔＲ，直線ｍと点Ｔとの間の距離をＲ
err ，走査線ｎと走査線ｎ＋１との間の距離をΔθ、走
査線ｎと点Ｔとの間の距離をθerr とする。
【０００６】このとき、以下の（１）式に示すようにし
て、先ず２点Ｓ＿ｎ＿ｍ，Ｓ＿ｎ＿ｍ＋１のサンプリン
グデータＳ（ｎ，ｍ），Ｓ（ｎ，ｍ＋１）を用いて直線
補間処理を行なうことにより点Ｔ＿ｎの補間データＴ
（ｎ，ｍ）を求めるとともに、２点Ｓ＿ｎ＋１＿ｍ，Ｓ
＿ｎ＋１＿ｍ＋１をのサンプリングデータＳ（ｎ＋１，
ｍ），Ｓ（ｎ＋１，ｍ＋１）用いて直線補間処理を行な
うことにより点Ｔ＿ｎ＋１の補間データＴ（ｎ＋１，
ｍ）を求める。
て、先ず２点Ｓ＿ｎ＿ｍ，Ｓ＿ｎ＿ｍ＋１のサンプリン
グデータＳ（ｎ，ｍ），Ｓ（ｎ，ｍ＋１）を用いて直線
補間処理を行なうことにより点Ｔ＿ｎの補間データＴ
（ｎ，ｍ）を求めるとともに、２点Ｓ＿ｎ＋１＿ｍ，Ｓ
＿ｎ＋１＿ｍ＋１をのサンプリングデータＳ（ｎ＋１，
ｍ），Ｓ（ｎ＋１，ｍ＋１）用いて直線補間処理を行な
うことにより点Ｔ＿ｎ＋１の補間データＴ（ｎ＋１，
ｍ）を求める。
【０００７】
【数１】
【０００８】次に（１）式により求められた２個の補間
データＴ（ｎ，ｍ），Ｔ（ｎ＋１，ｍ）を用いて以下の
（２）式に基づく直線補間処理を行なうことにより点Ｔ
に対応する表示用データＴ（ｘ，ｙ）を求める。
データＴ（ｎ，ｍ），Ｔ（ｎ＋１，ｍ）を用いて以下の
（２）式に基づく直線補間処理を行なうことにより点Ｔ
に対応する表示用データＴ（ｘ，ｙ）を求める。
【０００９】
【数２】
【００１０】このような補間処理を断層像全体について
行なうことにより、表示用の断層像が生成される。上記
のような補間処理方法はバイリニア補間と呼ばれてお
り、従来から知られている方法である。
行なうことにより、表示用の断層像が生成される。上記
のような補間処理方法はバイリニア補間と呼ばれてお
り、従来から知られている方法である。
【００１１】
【発明が解決しようとする課題】しかし、従来の補間処
理方法では、超音波走査線上に輝度の高い直線が現れ、
虚像として表示されることがある。特に、最近の超音波
診断装置では、表示画像の拡大表示を行う機能を有して
いる。セクタ走査，コンベックス走査により得られた断
層像は、走査線が略扇状に広がるため、被検体内の深度
の深い箇所では走査線の間隔が粗くなる。この時、直線
補間処理を行うと走査線上にデータのピークができ、こ
れにより輝度の高い直線が現れる場合があるという問題
点がある。この輝度の高い直線は、超音波サンプリング
データの値が大きい時直線補間を行うことにより現れ
る。
理方法では、超音波走査線上に輝度の高い直線が現れ、
虚像として表示されることがある。特に、最近の超音波
診断装置では、表示画像の拡大表示を行う機能を有して
いる。セクタ走査，コンベックス走査により得られた断
層像は、走査線が略扇状に広がるため、被検体内の深度
の深い箇所では走査線の間隔が粗くなる。この時、直線
補間処理を行うと走査線上にデータのピークができ、こ
れにより輝度の高い直線が現れる場合があるという問題
点がある。この輝度の高い直線は、超音波サンプリング
データの値が大きい時直線補間を行うことにより現れ
る。
【００１２】本発明は、この問題点を解決し、虚像とし
ての線が現われにくいデータ補間処理を行うデータ補間
処理装置を提供することを目的とする。
ての線が現われにくいデータ補間処理を行うデータ補間
処理装置を提供することを目的とする。
【００１３】
【課題を解決するための手段】上記目的を達成する本発
明のデータ補間処理装置は、走査線を形成する超音波ビ
ームにより被検体内を所定の走査方向に走査しながら該
被検体内で反射した超音波を受信しサンプリングするこ
とにより得られたサンプリングデータを補間処理して表
示用データを算出するデータ補間処理装置において、表
示用データの算出を所望する点を挟んで走査方向に並ぶ
前後２本ずつ合計４本の走査線上の、該点を挟んで走査
線の延びる方向に並ぶ、各走査線について２個ずつ合計
８個のサンプリングデータを用いて、上記所望する点の
表示用データを求める補間演算部を備えたことを特徴と
するものである。
明のデータ補間処理装置は、走査線を形成する超音波ビ
ームにより被検体内を所定の走査方向に走査しながら該
被検体内で反射した超音波を受信しサンプリングするこ
とにより得られたサンプリングデータを補間処理して表
示用データを算出するデータ補間処理装置において、表
示用データの算出を所望する点を挟んで走査方向に並ぶ
前後２本ずつ合計４本の走査線上の、該点を挟んで走査
線の延びる方向に並ぶ、各走査線について２個ずつ合計
８個のサンプリングデータを用いて、上記所望する点の
表示用データを求める補間演算部を備えたことを特徴と
するものである。
【００１４】図１は本発明の構成を表わすブロック図で
あり、この図１を参照しながら本発明の構成について詳
細に説明する。図１（Ａ）は、本発明のデータ補間処理
装置の第１の態様を表わし、その補間演算部１００は、
第１の直線補間演算手段１０２と第１の３次曲線補間演
算手段１０４により構成されている。
あり、この図１を参照しながら本発明の構成について詳
細に説明する。図１（Ａ）は、本発明のデータ補間処理
装置の第１の態様を表わし、その補間演算部１００は、
第１の直線補間演算手段１０２と第１の３次曲線補間演
算手段１０４により構成されている。
【００１５】サンプリングデータが補間演算部１００に
入力され、この入力されたサンプリングデータは第１の
直線補間演算手段１０２に入力される。この第１の直線
補間演算手段１０２では、表示用データの算出を所望す
る点を挟んで走査方向に並ぶ前後２本ずつ、互いに隣接
する合計４本の各走査線上にそれぞれ２個ずつ並ぶサン
プリングデータに各走査線毎に直線補間処理が施され、
これにより各走査線について１個ずつ合計４個の直線補
間データが求められる。
入力され、この入力されたサンプリングデータは第１の
直線補間演算手段１０２に入力される。この第１の直線
補間演算手段１０２では、表示用データの算出を所望す
る点を挟んで走査方向に並ぶ前後２本ずつ、互いに隣接
する合計４本の各走査線上にそれぞれ２個ずつ並ぶサン
プリングデータに各走査線毎に直線補間処理が施され、
これにより各走査線について１個ずつ合計４個の直線補
間データが求められる。
【００１６】求められた合計４個の直線補間データは、
第１の３次曲線補間演算手段１０４に入力される。この
第１の３次曲線補間演算手段１０４では、入力された４
個の直線補間データに３次曲線補間処理が施され、これ
により上述した所望する点の表示用データが求められ
る。図１（Ｂ）は本発明のデータ補間処理装置の第２の
態様を表わし、その補間演算部２００は、第２の３次曲
線補間演算手段２０２と第２の直線補間演算手段２０４
により構成されている。
第１の３次曲線補間演算手段１０４に入力される。この
第１の３次曲線補間演算手段１０４では、入力された４
個の直線補間データに３次曲線補間処理が施され、これ
により上述した所望する点の表示用データが求められ
る。図１（Ｂ）は本発明のデータ補間処理装置の第２の
態様を表わし、その補間演算部２００は、第２の３次曲
線補間演算手段２０２と第２の直線補間演算手段２０４
により構成されている。
【００１７】補間演算部２００に入力されたサンプリン
グデータは、第２の３次曲線補間演算手段２０２に入力
され、この第２の３次曲線補間演算手段２０２では、走
査方向に４個ずつ２列に並ぶサンプリングデータに各列
毎に３次曲線補間処理が施され、これにより各列につい
て１個ずつ合計２個の３次曲線補間データが求められ
る。
グデータは、第２の３次曲線補間演算手段２０２に入力
され、この第２の３次曲線補間演算手段２０２では、走
査方向に４個ずつ２列に並ぶサンプリングデータに各列
毎に３次曲線補間処理が施され、これにより各列につい
て１個ずつ合計２個の３次曲線補間データが求められ
る。
【００１８】求められた合計２個の３次曲線補間データ
は、第２の直線補間演算手段２０４に入力される。この
第２の直線補間演算手段２０４では、第２の３次曲線補
間演算手段２０２で求められた２個の３次曲線補間デー
タに直線補間処理が施され、上述した所望する点の表示
用データが求められる。尚、上記の説明は、所望する点
一点についての補間演算の説明であり、このような補間
演算が画像上の多数の点にわたって行なわれるものであ
ることはいうまでもない。
は、第２の直線補間演算手段２０４に入力される。この
第２の直線補間演算手段２０４では、第２の３次曲線補
間演算手段２０２で求められた２個の３次曲線補間デー
タに直線補間処理が施され、上述した所望する点の表示
用データが求められる。尚、上記の説明は、所望する点
一点についての補間演算の説明であり、このような補間
演算が画像上の多数の点にわたって行なわれるものであ
ることはいうまでもない。
【００１９】
【作用】一般的に、超音波診断装置では、走査線の延び
る深さ方向の分解能は充分にあり、走査方向は比較的分
解能が悪く、このため走査方向については補間による誤
差が大きい。本発明のデータ補間処理装置は、所望する
点の表示用データを求めるにあたり、走査方向に４点、
深さ方向に２点の合計８点のサンプリングデータを用い
て補間演算処理を行なう補間演算部を備え、深さ方向に
ついては直線補間を行ない、走査方向については３次曲
線補間を行なう構成としたものである。
る深さ方向の分解能は充分にあり、走査方向は比較的分
解能が悪く、このため走査方向については補間による誤
差が大きい。本発明のデータ補間処理装置は、所望する
点の表示用データを求めるにあたり、走査方向に４点、
深さ方向に２点の合計８点のサンプリングデータを用い
て補間演算処理を行なう補間演算部を備え、深さ方向に
ついては直線補間を行ない、走査方向については３次曲
線補間を行なう構成としたものである。
【００２０】図１（Ａ）に示す第１の態様では、第１の
直線補間演算手段１０２により深さ方向に直線補間が行
なわれ、次いで第１の３次曲線補間演算手段１０４によ
り走査方向に３次曲線補間が行なわれる。これにより所
望する点の表示用データが従来よりも高精度に求めら
れ、このようにして求めた表示用データに基づく画像
は、虚像の現われにくい画像となる。
直線補間演算手段１０２により深さ方向に直線補間が行
なわれ、次いで第１の３次曲線補間演算手段１０４によ
り走査方向に３次曲線補間が行なわれる。これにより所
望する点の表示用データが従来よりも高精度に求めら
れ、このようにして求めた表示用データに基づく画像
は、虚像の現われにくい画像となる。
【００２１】また図１（Ｂ）に示す第２の態様では、第
２の３次曲線演算手段２０２により、先ず走査方向に３
次曲線補間が行なわれ、次いで第２の直線補間演算手段
２０４により深さ方向に直線補間が行なわれる。これに
より、上記第１の態様と同様に所望する点の表示用デー
タが高精度に求められる。
２の３次曲線演算手段２０２により、先ず走査方向に３
次曲線補間が行なわれ、次いで第２の直線補間演算手段
２０４により深さ方向に直線補間が行なわれる。これに
より、上記第１の態様と同様に所望する点の表示用デー
タが高精度に求められる。
【００２２】
【実施例】図２は、本発明に係るデータ補間処理演算の
説明図である。符号の付け方は、図７の場合と同様であ
る。ここでは、所望する点Ｔの表示用データＴ（ｘ，
ｙ）を求めるにあたり、この点Ｔを挟む左右２本、合計
４本の走査線ｎ，ｎ＋１，ｎ＋２，ｎ＋３の走査線上に
あり、かつ深さ方向にこの点Ｔを挟み２列に並ぶ（直線
ｍと直線ｍ＋１上にある）合計８点Ｓ＿ｎ＿ｍ，Ｓ＿ｎ
＿ｍ＋１，Ｓ＿ｎ＋１＿ｍ，Ｓ＿ｎ＋１＿ｍ＋１，Ｓ＿
ｎ＋２＿ｍ，Ｓ＿ｎ＋２＿ｍ＋１，Ｓ＿ｎ＋３＿ｍ，Ｓ
＿ｎ＋３＿ｍ＋１にそれぞれ対応するサンプリングデー
タＳ（ｎ，ｍ），Ｓ（ｎ，ｍ＋１），Ｓ（ｎ＋１，
ｍ），Ｓ（ｎ＋１，ｍ＋１），Ｓ（ｎ＋２，ｍ），Ｓ
（ｎ＋２，ｍ＋１），Ｓ（ｎ＋３，ｍ），Ｓ（ｎ＋３，
ｍ＋１）を使用し、超音波走査線の深さ方向に直線補間
処理
説明図である。符号の付け方は、図７の場合と同様であ
る。ここでは、所望する点Ｔの表示用データＴ（ｘ，
ｙ）を求めるにあたり、この点Ｔを挟む左右２本、合計
４本の走査線ｎ，ｎ＋１，ｎ＋２，ｎ＋３の走査線上に
あり、かつ深さ方向にこの点Ｔを挟み２列に並ぶ（直線
ｍと直線ｍ＋１上にある）合計８点Ｓ＿ｎ＿ｍ，Ｓ＿ｎ
＿ｍ＋１，Ｓ＿ｎ＋１＿ｍ，Ｓ＿ｎ＋１＿ｍ＋１，Ｓ＿
ｎ＋２＿ｍ，Ｓ＿ｎ＋２＿ｍ＋１，Ｓ＿ｎ＋３＿ｍ，Ｓ
＿ｎ＋３＿ｍ＋１にそれぞれ対応するサンプリングデー
タＳ（ｎ，ｍ），Ｓ（ｎ，ｍ＋１），Ｓ（ｎ＋１，
ｍ），Ｓ（ｎ＋１，ｍ＋１），Ｓ（ｎ＋２，ｍ），Ｓ
（ｎ＋２，ｍ＋１），Ｓ（ｎ＋３，ｍ），Ｓ（ｎ＋３，
ｍ＋１）を使用し、超音波走査線の深さ方向に直線補間
処理
【００２３】
【数３】
【００２４】を計算し、４点Ｔ＿ｎ，Ｔ＿ｎ＋１，Ｔ＿
ｎ＋２，Ｔ＿ｎ＋３の直線補間データＴ（ｎ，ｍ），Ｔ
（ｎ＋１，ｍ），Ｔ（ｎ＋２，ｍ），Ｔ（ｎ＋３，ｍ）
を求める。この補間処理は従来と同じ処理である。次に
この４点Ｔ＿ｎ，Ｔ＿ｎ＋１，Ｔ＿ｎ＋２，Ｔ＿ｎ＋３
の直線補間データＴ（ｎ，ｍ），Ｔ（ｎ＋１，ｍ），Ｔ
（ｎ＋２，ｍ），Ｔ（ｎ＋３，ｍ）から３次曲線補間処
理を行なう。３次曲線補間処理は
ｎ＋２，Ｔ＿ｎ＋３の直線補間データＴ（ｎ，ｍ），Ｔ
（ｎ＋１，ｍ），Ｔ（ｎ＋２，ｍ），Ｔ（ｎ＋３，ｍ）
を求める。この補間処理は従来と同じ処理である。次に
この４点Ｔ＿ｎ，Ｔ＿ｎ＋１，Ｔ＿ｎ＋２，Ｔ＿ｎ＋３
の直線補間データＴ（ｎ，ｍ），Ｔ（ｎ＋１，ｍ），Ｔ
（ｎ＋２，ｍ），Ｔ（ｎ＋３，ｍ）から３次曲線補間処
理を行なう。３次曲線補間処理は
【００２５】
【数４】
【００２６】を計算することにより、求められる。この
演算は点Ｔ＿ｎ＋１を原点と考え、Δθ＝１（単位距
離）と考え、以下の４つの式
演算は点Ｔ＿ｎ＋１を原点と考え、Δθ＝１（単位距
離）と考え、以下の４つの式
【００２７】
【数５】 Ｔ（ｎ，ｍ）＝（−ａ）＋ｂ−ｃ＋ｄ …（９） Ｔ（ｎ＋１，ｍ）＝ｄ …（１０） Ｔ（ｎ＋２，ｍ）＝ａ＋ｂ＋ｃ＋ｄ …（１１） Ｔ（ｎ＋３，ｍ）＝８ａ＋４ｂ＋２ｃ＋ｄ …（１２） の連立方程式を解くことにより求められる。
【００２８】図３は上記の演算を実現する、本発明のデ
ータ補間処理装置の第１の態様の一実施例の回路図であ
る。図中、○内に＋の記号、○内に−の記号、および○
内に×の記号は、入力された２つの信号をそれぞれ加
算、減算、乗算するそれぞれ加算器、減算器、乗算器を
表わし、矢印状の五角形は、入力信号にその五角形内の
数字が掛け算されたと同等になるように入力信号をＬＳ
Ｂ側もしくはＭＳＢ側にシフトするシフタを表わしてい
る。
ータ補間処理装置の第１の態様の一実施例の回路図であ
る。図中、○内に＋の記号、○内に−の記号、および○
内に×の記号は、入力された２つの信号をそれぞれ加
算、減算、乗算するそれぞれ加算器、減算器、乗算器を
表わし、矢印状の五角形は、入力信号にその五角形内の
数字が掛け算されたと同等になるように入力信号をＬＳ
Ｂ側もしくはＭＳＢ側にシフトするシフタを表わしてい
る。
【００２９】各２つのサンプリングデータＳ（ｎ，
ｍ），Ｓ（ｎ，ｍ＋１）；Ｓ（ｎ＋１，ｍ），Ｓ（ｎ＋
１，ｍ＋１）；Ｓ（ｎ＋２，ｍ），Ｓ（ｎ＋２，ｍ＋
１）；Ｓ（ｎ＋３，ｍ），Ｓ（ｎ＋３，ｍ＋１）はそれ
ぞれ各直線補間処理部１０，２０，３０，４０に入力さ
れる。また各直線補間処理部１０，２０，３０，４０に
は、係数Ｒerr ／ΔＲも入力される。さらにこの係数Ｒ
err ／ΔＲと係数“１”は減算器１に入力され、この減
算器では１−（Ｒerr ／ΔＲ）が演算されてこの演算結
果も各直線補間部１０，２０，３０，４０に入力され
る。各直線補間部１０，２０，３０，４０では、それぞ
れ、（３）式でｋ＝ｎ，ｎ＋１，ｎ＋２，ｎ＋３と置い
た場合の直線補間演算が行なわれ各直線補間データＴ
（ｎ，ｍ），Ｔ（ｎ＋１，ｍ），Ｔ（ｎ＋２，ｍ），Ｔ
（ｎ＋３，ｍ）が算出される。図３に示す実施例では各
直線補間部１０，２０，３０，４０を合わせたものが本
発明にいう第１の直線補間演算手段１０２（図１参照）
の一例であると観念される。
ｍ），Ｓ（ｎ，ｍ＋１）；Ｓ（ｎ＋１，ｍ），Ｓ（ｎ＋
１，ｍ＋１）；Ｓ（ｎ＋２，ｍ），Ｓ（ｎ＋２，ｍ＋
１）；Ｓ（ｎ＋３，ｍ），Ｓ（ｎ＋３，ｍ＋１）はそれ
ぞれ各直線補間処理部１０，２０，３０，４０に入力さ
れる。また各直線補間処理部１０，２０，３０，４０に
は、係数Ｒerr ／ΔＲも入力される。さらにこの係数Ｒ
err ／ΔＲと係数“１”は減算器１に入力され、この減
算器では１−（Ｒerr ／ΔＲ）が演算されてこの演算結
果も各直線補間部１０，２０，３０，４０に入力され
る。各直線補間部１０，２０，３０，４０では、それぞ
れ、（３）式でｋ＝ｎ，ｎ＋１，ｎ＋２，ｎ＋３と置い
た場合の直線補間演算が行なわれ各直線補間データＴ
（ｎ，ｍ），Ｔ（ｎ＋１，ｍ），Ｔ（ｎ＋２，ｍ），Ｔ
（ｎ＋３，ｍ）が算出される。図３に示す実施例では各
直線補間部１０，２０，３０，４０を合わせたものが本
発明にいう第１の直線補間演算手段１０２（図１参照）
の一例であると観念される。
【００３０】次に、得られたＴ（ｎ，ｍ），Ｔ（ｎ＋
１，ｍ），Ｔ（ｎ＋２，ｍ），Ｔ（ｎ＋３，ｍ）を用い
て３次曲線補間処理を行う。（５）式を計算する係数演
算部５０で係数ａが求められる。ここで、乗算器５１は
値を１／６にするためのものである。また、値を１／２
にするためのシフタ５２は入力されたデータを１ビット
ＬＳＢ側にシフトして加算器５３に送り出すものであ
る。このシフタ５２により乗算器等を使用せずに値が１
／２にされる。同様に係数ｂを求める係数演算部６０で
も値を１／２にするためにデータを１ビットＬＳＢ側に
シフトするシフタ６１を経由した後減算器で減算され
る。また、係数ｃを求める係数演算部７０における、値
を１／２にするシフタ７１、値を２倍するシフタ７２で
はそれぞれＬＳＢ側、ＭＳＢ側にビットシフトが行なわ
れる。尚、係数ｄは、（１０）式に示すように、直線補
間データＴ（ｎ＋１，ｍ）がそのまま用いられる。
１，ｍ），Ｔ（ｎ＋２，ｍ），Ｔ（ｎ＋３，ｍ）を用い
て３次曲線補間処理を行う。（５）式を計算する係数演
算部５０で係数ａが求められる。ここで、乗算器５１は
値を１／６にするためのものである。また、値を１／２
にするためのシフタ５２は入力されたデータを１ビット
ＬＳＢ側にシフトして加算器５３に送り出すものであ
る。このシフタ５２により乗算器等を使用せずに値が１
／２にされる。同様に係数ｂを求める係数演算部６０で
も値を１／２にするためにデータを１ビットＬＳＢ側に
シフトするシフタ６１を経由した後減算器で減算され
る。また、係数ｃを求める係数演算部７０における、値
を１／２にするシフタ７１、値を２倍するシフタ７２で
はそれぞれＬＳＢ側、ＭＳＢ側にビットシフトが行なわ
れる。尚、係数ｄは、（１０）式に示すように、直線補
間データＴ（ｎ＋１，ｍ）がそのまま用いられる。
【００３１】各係数ａ，ｂ，ｃ，ｄが求められた後、乗
算器８１でａ×（θerr ／Δθ）3が計算される。ま
た、乗算器８２でｂ×（θerr ／Δθ）2 、乗算器８３
でｃ×（θerr ／Δθ）が計算され、加算器９１，９
２，９３で互いに加算される。尚、乗算器８４で（θ
err ／Δθ）の二乗が計算され、乗算器８５で（θerr
／Δθ）の三乗が計算される。
算器８１でａ×（θerr ／Δθ）3が計算される。ま
た、乗算器８２でｂ×（θerr ／Δθ）2 、乗算器８３
でｃ×（θerr ／Δθ）が計算され、加算器９１，９
２，９３で互いに加算される。尚、乗算器８４で（θ
err ／Δθ）の二乗が計算され、乗算器８５で（θerr
／Δθ）の三乗が計算される。
【００３２】ここで、この実施例では、１／２、２倍の
計算を、ビットをＬＳＢ側，ＭＳＢ側にシフトするシフ
タで行ったが乗算器を使用してもよい。また、１／６の
演算も、シフタと加算器の組み合わせで実現することも
容易である。図４は、本発明のデータ補間処理装置の第
１の態様の他の実施例の回路図である。
計算を、ビットをＬＳＢ側，ＭＳＢ側にシフトするシフ
タで行ったが乗算器を使用してもよい。また、１／６の
演算も、シフタと加算器の組み合わせで実現することも
容易である。図４は、本発明のデータ補間処理装置の第
１の態様の他の実施例の回路図である。
【００３３】上述した（３）式をまとめ直すと、
【００３４】
【数６】
【００３５】となる。図４に示す実施例は、直線補間処
理を（１３）式で行ない、乗算器１個で各直線補間部１
１０，１２０，１３０，１４０を実現したものである。
この実施例では１−（Ｒerr ／Δθ）の演算（図３に示
す減算器１）は不要である。この実施例における３次曲
線補間演算処理部の構成は、図３に示す実施例における
それと同一である。
理を（１３）式で行ない、乗算器１個で各直線補間部１
１０，１２０，１３０，１４０を実現したものである。
この実施例では１−（Ｒerr ／Δθ）の演算（図３に示
す減算器１）は不要である。この実施例における３次曲
線補間演算処理部の構成は、図３に示す実施例における
それと同一である。
【００３６】図５は、本発明のデータ補間処理装置の第
１の態様の、さらに異なる実施例の回路図である。この
回路図は、３次曲線補間処理を以下の式で実現した場合
のものである。
１の態様の、さらに異なる実施例の回路図である。この
回路図は、３次曲線補間処理を以下の式で実現した場合
のものである。
【００３７】
【数７】
【００３８】この演算は、点Ｔ＿ｎを原点と考え、Δθ
＝１（単位距離）と考え、以下の４つの式
＝１（単位距離）と考え、以下の４つの式
【００３９】
【数８】 Ｔ（ｎ，ｍ）＝ｄ …（１９） Ｔ（ｎ＋１，ｍ）＝ａ＋ｂ＋ｃ＋ｄ …（２０） Ｔ（ｎ＋２，ｍ）＝８ａ＋４ｂ＋２ｃ＋ｄ …（２１） Ｔ（ｎ＋３，ｍ）＝２７ａ＋９ｂ＋３ｃ＋ｄ …（２２） の連立方程式を解くことにより求められる。
【００４０】係数ａを求める係数演算部１５０におい
て、値を１／６にするには乗算器１５１を使用し、値を
１／２にするにはシフタ１５２でＬＳＢ側に１ビットシ
フトし、加算器１５３で加算する。また、係数ｂを求め
る係数演算部１６０ではシフタ１６１で１ビットＭＳＢ
側にシフトすることにより入力が２倍にされる。また値
を５倍とするにはビットシフト１６８で２ビットＭＳＢ
側にシフトして４倍し、これに加算器１６４で元の値を
加算する事により得られる。値を１／２するには他の場
合と同様にシフタ１６６で行なわれる。同様に、係数ｃ
を求める係数演算部１７０において、値を１１倍するた
めにシフタ１７８で３ビットＭＳＢ側にシフトした値と
１ビットシフトした値と元の値とが加算器１７７，１８
０，１８１で加算される。本実施例はこのように乗算器
を極力使用せずに構成されており、これにより回路規模
（ＬＳＩ等のゲート数）を削減することができる。ただ
し乗算器を用いて各演算を実現することも可能である。
て、値を１／６にするには乗算器１５１を使用し、値を
１／２にするにはシフタ１５２でＬＳＢ側に１ビットシ
フトし、加算器１５３で加算する。また、係数ｂを求め
る係数演算部１６０ではシフタ１６１で１ビットＭＳＢ
側にシフトすることにより入力が２倍にされる。また値
を５倍とするにはビットシフト１６８で２ビットＭＳＢ
側にシフトして４倍し、これに加算器１６４で元の値を
加算する事により得られる。値を１／２するには他の場
合と同様にシフタ１６６で行なわれる。同様に、係数ｃ
を求める係数演算部１７０において、値を１１倍するた
めにシフタ１７８で３ビットＭＳＢ側にシフトした値と
１ビットシフトした値と元の値とが加算器１７７，１８
０，１８１で加算される。本実施例はこのように乗算器
を極力使用せずに構成されており、これにより回路規模
（ＬＳＩ等のゲート数）を削減することができる。ただ
し乗算器を用いて各演算を実現することも可能である。
【００４１】尚図５に示す実施例では、直線補間部１１
０，１２０，１３０，１４０は、図４に示す実施例の直
線補間部１１０，１２０，１３０，１４０とそれぞれ同
一のものを採用したが、図３に示す構成の直線補間部１
０，２０，３０，４０を採用してもよいことはもちろん
である。図６は、本発明のデータ補間処理装置の第２の
態様の一実施例を表わした回路図である。この実施例
は、第２の態様の実施例であるため、先に３次曲線補間
処理を行って２点のデータを求め、次に、その２点から
直線補間処理により表示データを求めるものである。こ
の図６において、３次曲線係数ａ，ｂ，ｃを求める係数
演算部５０，６０，７０は、それぞれ図３に示す係数演
算部５０，６０，７０の構成と同じである。またその後
の乗算器８１，８２，８３と加算器９１，９２，９３か
らなる回路も図３に示す構成と同じである。直線補間演
算部１１０は図４に示す直線補間演算部１１０，１２
０，１３０，１４０と同じ構成のものを示しているが、
図３に示す直線補間部１０，２０，３０，４０と同じ構
成のものを採用してもよい。
０，１２０，１３０，１４０は、図４に示す実施例の直
線補間部１１０，１２０，１３０，１４０とそれぞれ同
一のものを採用したが、図３に示す構成の直線補間部１
０，２０，３０，４０を採用してもよいことはもちろん
である。図６は、本発明のデータ補間処理装置の第２の
態様の一実施例を表わした回路図である。この実施例
は、第２の態様の実施例であるため、先に３次曲線補間
処理を行って２点のデータを求め、次に、その２点から
直線補間処理により表示データを求めるものである。こ
の図６において、３次曲線係数ａ，ｂ，ｃを求める係数
演算部５０，６０，７０は、それぞれ図３に示す係数演
算部５０，６０，７０の構成と同じである。またその後
の乗算器８１，８２，８３と加算器９１，９２，９３か
らなる回路も図３に示す構成と同じである。直線補間演
算部１１０は図４に示す直線補間演算部１１０，１２
０，１３０，１４０と同じ構成のものを示しているが、
図３に示す直線補間部１０，２０，３０，４０と同じ構
成のものを採用してもよい。
【００４２】
【発明の効果】以上説明したように本発明のデータ補間
演算装置は、走査方向に４点、深さ方向に２点、合計８
点のサンプリングデータを用いて所望の点の表示用デー
タを求める補間演算部を備えたものであるため、従来の
バイリニア補間により現われていた輝度の高い虚像を低
減もしくは消滅させることができ、画像拡大時にも滑ら
かな補間画像が得られる。
演算装置は、走査方向に４点、深さ方向に２点、合計８
点のサンプリングデータを用いて所望の点の表示用デー
タを求める補間演算部を備えたものであるため、従来の
バイリニア補間により現われていた輝度の高い虚像を低
減もしくは消滅させることができ、画像拡大時にも滑ら
かな補間画像が得られる。
【図１】本発明の構成を表わすブロック図である。
【図２】本発明に係るデータ補間処理演算の説明図であ
る。
る。
【図３】本発明のデータ補間処理装置の第１の態様の一
実施例の回路図である。
実施例の回路図である。
【図４】本発明のデータ補間処理装置の第１の態様の他
の実施例の回路図である。
の実施例の回路図である。
【図５】本発明のデータ補間処理装置の第１の態様の、
さらに異なる実施例の回路図である。
さらに異なる実施例の回路図である。
【図６】本発明のデータ補間処理装置の第２の態様の一
実施例の回路図である。
実施例の回路図である。
【図７】従来のデータ補間処理装置の説明図である。
１００ 補間演算部 １０２ 第１の直線補間演算手段 １０４ 第１の３次曲線補間演算手段 ２００ 補間演算部 ２０２ 第２の３次曲線補間演算手段 ２０４ 第２の直線補間演算手段 １０，２０，３０，４０ 直線補間部 １１０，１２０，１３０，１４０ 直線補間部 ５０，６０，７０ 係数演算部 １５０，１６０，１７０ 係数演算部 ８１，８２，８３，８４，８５ 乗算器 ９１，９２，９３ 加算器
───────────────────────────────────────────────────── フロントページの続き (56)参考文献 特開 昭57−108681（ＪＰ，Ａ) 特開 昭62−179438（ＪＰ，Ａ) 特開 昭64−64635（ＪＰ，Ａ) 特開 平２−283358（ＪＰ，Ａ) 特開 平４−84958（ＪＰ，Ａ) 特開 平４−244956（ＪＰ，Ａ) (58)調査した分野(Int.Cl.7，ＤＢ名) A61B 8/00 - 8/15
Claims (2)
- 【請求項１】 走査線を形成する超音波ビームにより被
検体内を所定の走査方向に走査しながら該被検体内で反
射した超音波を受信しサンプリングすることにより得ら
れたサンプリングデータを補間処理して表示用データを
算出するデータ補間処理装置において、 表示用データの算出を所望する点を挟んで走査方向に並
ぶ前後２本ずつ合計４本の走査線上の、該点を挟んで走
査線の延びる方向に並ぶ、各走査線について２個ずつ合
計８個のサンプリングデータを用いて、前記点の表示用
データを求める補間演算部を備え、 前記補間演算部が、 各走査線上に２個ずつ並ぶサンプリングデータに各走査
線毎に直線補間処理を施すことにより各走査線について
１個ずつ合計４個の直線補間データを求める第１の直線
補間演算手段と、 該第１の直線補間演算手段で求められた４個の直線補間
データに３次曲線補間処理を施すことにより前記点の表
示用データを求める第１の３次曲線補間演算手段とを備
え たことを特徴とするデータ補間処理装置。 - 【請求項２】 走査線を形成する超音波ビームにより被
検体内を所定の走査方向に走査しながら該被検体内で反
射した超音波を受信しサンプリングすることにより得ら
れたサンプリングデータを補間処理して表示用データを
算出するデータ補間処理装置において、 表示用データの算出を所望する点を挟んで走査方向に並
ぶ前後２本ずつ合計４本の走査線上の、該点を挟んで走
査線の延びる方向に並ぶ、各走査線について２個ずつ合
計８個のサンプリングデータを用いて、前記点の表示用
データを求める補間演算部を備え、 前記補間演算部が、 走査方向に４個ずつ２列に並ぶサンプリングデータに各
列毎に３次曲線補間処理を施すことにより各列について
１個ずつ合計２個の３次曲線補間データを求める第２の
３次曲線補間演算手段と、 該第２の３次曲線補間演算手段で求められた２個の３次
曲線補間データに直線 補間処理を施すことにより前記点
の表示用データを求める第２の直線補間演算手段 とを備
えたことを特徴とするデータ補間処理装置。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP27628892A JP3223202B2 (ja) | 1992-10-14 | 1992-10-14 | データ補間処理装置 |
| US08/123,078 US5388582A (en) | 1992-10-14 | 1993-09-20 | Data interpolation processing apparatus |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP27628892A JP3223202B2 (ja) | 1992-10-14 | 1992-10-14 | データ補間処理装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH06125897A JPH06125897A (ja) | 1994-05-10 |
| JP3223202B2 true JP3223202B2 (ja) | 2001-10-29 |
Family
ID=17567365
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP27628892A Expired - Fee Related JP3223202B2 (ja) | 1992-10-14 | 1992-10-14 | データ補間処理装置 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US5388582A (ja) |
| JP (1) | JP3223202B2 (ja) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6193661B1 (en) * | 1999-04-07 | 2001-02-27 | Agilent Technologies, Inc. | System and method for providing depth perception using single dimension interpolation |
| JP2005058587A (ja) * | 2003-08-19 | 2005-03-10 | Ge Medical Systems Global Technology Co Llc | 超音波画像作成方法および超音波診断装置 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5672857A (en) * | 1979-11-16 | 1981-06-17 | Matsushita Electric Ind Co Ltd | Method of scanning ultrasonic diagnosing device |
| JPS6443238A (en) * | 1987-08-12 | 1989-02-15 | Toshiba Corp | Ultrasonic blood flow imaging apparatus |
| JP2763126B2 (ja) * | 1989-02-10 | 1998-06-11 | 株式会社東芝 | カラー超音波診断装置 |
-
1992
- 1992-10-14 JP JP27628892A patent/JP3223202B2/ja not_active Expired - Fee Related
-
1993
- 1993-09-20 US US08/123,078 patent/US5388582A/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| US5388582A (en) | 1995-02-14 |
| JPH06125897A (ja) | 1994-05-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5409007A (en) | Filter to reduce speckle artifact in ultrasound imaging | |
| US20050143655A1 (en) | Ultrasonic image processing apparatus, ultrasonic image processing method and ultrasonic image processing program | |
| JPH10290801A (ja) | 高解像度超音波画像処理方法および装置 | |
| JP2010017557A (ja) | 超音波データを処理する超音波システム及び方法 | |
| JP2022515810A (ja) | コントラスト強調画像のためのシステム及び方法 | |
| JP4590256B2 (ja) | 超音波撮像装置、超音波画像処理方法、及び、超音波画像処理プログラム | |
| JP3223202B2 (ja) | データ補間処理装置 | |
| JPH07116162A (ja) | 超音波計測装置およびその送受信位相補正方法 | |
| JP3806229B2 (ja) | 超音波診断装置 | |
| US7946991B2 (en) | Ultrasonic doppler blood flow measuring device | |
| JP2006212054A (ja) | 超音波観測装置、画像処理装置、及び、画像処理プログラム | |
| JPH02164352A (ja) | 超音波診断画像処理装置 | |
| JP3451025B2 (ja) | 超音波ドプラ診断装置 | |
| JP3324992B2 (ja) | 超音波診断装置 | |
| JP7034686B2 (ja) | 超音波診断装置、医用画像処理装置及びそのプログラム | |
| JPH07265314A (ja) | 超音波断層装置 | |
| JP2002186615A (ja) | 超音波診断装置 | |
| JP4666793B2 (ja) | 超音波診断装置 | |
| JP3413095B2 (ja) | 断層画像処理方法及び装置 | |
| JP2714329B2 (ja) | 超音波診断装置 | |
| KR101126181B1 (ko) | 초음파 영상 시스템의 화소점에 수신 집속하는 방법 및 장치 | |
| JP5348829B2 (ja) | 超音波診断装置、超音波画像表示方法及び超音波画像表示プログラム | |
| JP3403262B2 (ja) | 超音波診断装置 | |
| JP2000339441A (ja) | 座標変換装置 | |
| JP3300196B2 (ja) | 超音波診断装置及びビーム補間方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20010807 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20070817 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100817 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120817 Year of fee payment: 11 |
|
| LAPS | Cancellation because of no payment of annual fees |