JP5161413B2 - 液体充填式伝送システムの補正データを決定、モニタ及び更新するための方法 - Google Patents
液体充填式伝送システムの補正データを決定、モニタ及び更新するための方法 Download PDFInfo
- Publication number
- JP5161413B2 JP5161413B2 JP2002529960A JP2002529960A JP5161413B2 JP 5161413 B2 JP5161413 B2 JP 5161413B2 JP 2002529960 A JP2002529960 A JP 2002529960A JP 2002529960 A JP2002529960 A JP 2002529960A JP 5161413 B2 JP5161413 B2 JP 5161413B2
- Authority
- JP
- Japan
- Prior art keywords
- patient
- signal
- determined
- segment
- transfer function
- 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
- 238000000034 method Methods 0.000 title claims abstract description 62
- 238000012937 correction Methods 0.000 title claims abstract description 36
- 239000007788 liquid Substances 0.000 title claims abstract description 34
- 230000005540 biological transmission Effects 0.000 title claims abstract description 28
- 238000012544 monitoring process Methods 0.000 title claims abstract description 8
- 238000009530 blood pressure measurement Methods 0.000 claims abstract description 11
- 239000003814 drug Substances 0.000 claims abstract description 4
- 206010002091 Anaesthesia Diseases 0.000 claims abstract description 3
- 230000037005 anaesthesia Effects 0.000 claims abstract description 3
- 238000005259 measurement Methods 0.000 claims description 58
- 238000012546 transfer Methods 0.000 claims description 51
- 238000011156 evaluation Methods 0.000 claims description 6
- 238000009499 grossing Methods 0.000 claims description 6
- 238000013016 damping Methods 0.000 claims description 5
- 230000003595 spectral effect Effects 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000011218 segmentation Effects 0.000 claims 3
- 230000001131 transforming effect Effects 0.000 claims 2
- 210000001367 artery Anatomy 0.000 claims 1
- 238000009795 derivation Methods 0.000 claims 1
- 238000001914 filtration Methods 0.000 claims 1
- 238000013179 statistical model Methods 0.000 claims 1
- 210000003462 vein Anatomy 0.000 claims 1
- 238000001949 anaesthesia Methods 0.000 abstract 1
- 230000000875 corresponding effect Effects 0.000 description 16
- 238000001228 spectrum Methods 0.000 description 10
- 238000012545 processing Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 101100412394 Drosophila melanogaster Reg-2 gene Proteins 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 230000035939 shock Effects 0.000 description 2
- 230000009897 systematic effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000036772 blood pressure Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000013013 elastic material Substances 0.000 description 1
- 230000001771 impaired effect Effects 0.000 description 1
- 238000000338 in vitro Methods 0.000 description 1
- 230000036512 infertility Effects 0.000 description 1
- 238000001361 intraarterial administration Methods 0.000 description 1
- 238000001990 intravenous administration Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000011551 log transformation method Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 238000012549 training Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
- A61B5/7257—Details of waveform analysis characterised by using transforms using Fourier transforms
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/021—Measuring pressure in heart or blood vessels
- A61B5/0215—Measuring pressure in heart or blood vessels by means inserted into the body
- A61B5/02156—Calibration means
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Physiology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Cardiology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Mathematical Physics (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Vascular Medicine (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Arrangements For Transmission Of Measured Signals (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Radar Systems Or Details Thereof (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Description
この発明は請求項1の前文に記載されている、測定値の歪みを補正するとともに液体充填式伝送システムを校正するための補正データを決定、モニタ及び更新する方法に関する。
【0002】
心臓病学、集中治療医学及び麻酔術においては、侵襲的に動脈内及び静脈内の圧力測定を行うために、液体充填式の伝送システムが使用される。そこでは、圧力測定は患者の体内で行われ、カテーテルとして形成された液体充填式伝送システムによって、患者の体の外に配置された圧力トランスデューサまで伝送される。カテーテルの長さ、断面、構造及び弾性特性とカテーテルの中に入れられた液体の組成とによって、人体内部のカテーテル端部において患者から信号として検出された圧力の測定値に、様々な共振、減衰、エネルギ損失が生じ、これが液体充填式伝送システムのために患者からのこの信号に大きな歪みを引き起こす。こうした歪みのために、患者からの信号を定量的に解析することができず、診断やモニタリングの定性的な解釈を損なう。
【0003】
液体充填式伝送システムによる患者からの信号の歪みを避けるために、先端圧力センサとして知られる圧力測定トランスデューサが液体充填式伝送システムとして作用するカテーテルの先端に配置されてきた。患者から検出された信号は電気信号に変換され、電気ラインを介して患者の体からガイドされる。しかし、このタイプの先端圧力センサは非常に高価であり、また、その形状及び寸法のために、限られた範囲でしか利用できない。
【0004】
患者の体の内側で測定された患者からの信号の伝送における液体充填式伝送システムによる測定値の歪みをコンピュータ補正するための方法がDE 1982208844 A1に開示されている。そこでは、外部圧力トランスデューサから出力される患者からの(歪んだ)信号からの電気信号がアナログ/デジタル変換器の中へ通される。そして、出力されるデジタル化された信号は、デジタルフーリエ解析に基づいて動作してデジタル信号の心拍ごとの解析を実行する信号解析及び処理装置の中で解析される。解析された信号は次に、補正データレコード行列から、又はフーリエ係数として呼び出される実験的に決定される補正データレコードベクトルと連結されて、フーリエ係数として出力される。信号解析及び処理装置によって補正された信号は最後に出力及び/又は評価装置へ導かれる。
【0005】
この周知の方法に必要な補正データレコードは、特に、参照の圧力測定から決定することができる。伝達関数を決定するために、未知の測定信号の代わりに、ここでは人工的に発生した既知の試験圧力信号を参照信号として入力することができる。そして、その歪みから伝送システムの特性を求めることができる。
【0006】
一つの可能性は、カテーテル先端に配置されたキャリブレータから校正データレコードを得ることであるが、これには無菌状態を維持する必要があり、面倒な取り扱いが生じる。このことは、伝達関数の変化を考慮しなければならないことを意味している。なぜなら、圧力測定箇所までガイドされるカテーテルやキャリブレータがすすがれ、薬物がカテーテルを介して供給され、いわゆる“マイクロバブル”が伝達関数を変化させるからである。そうした変化が生じた場合には、再校正が必要である。しかし、水平方向に配置されたカテーテルの場合には、言い換えると、患者の体内に配置されたカテーテルにおいては、これは実行が容易ではない。
【0007】
別の可能性は、外部校正信号を用いて伝達関数から補正データレコードを得ることである。この目的のために、ジャンプ信号、圧力衝撃又はノイズの形の校正信号が測定箇所から離れた方のカテーテル側、言い換えると、患者の体外で発生させされ、そのプロセスで発生する固有の発振を用いて補正データレコードが計算される。しかし、外部校正信号を用いたこのタイプの伝達関数の決定では、いくつかの問題が生じる。一方では、校正信号、すなわち、ジャンプ信号、圧力衝撃又はノイズを非常に正確かつ再現性よく発生する必要がある。これは、伝送システムの測定が不安定であるため、かなり長期間にわたって非常に注意深い作業と測定を必要とする。従って、例えば、病院では、多忙と訓練用の入力のために、実行することができない。他方において、圧力トランスデューサに取り付けられた自動機械デバイスの使用は非常に高価であり、そのうえ精密機械部分のため、設定と保守のために特殊な作業を必要とする。
【0008】
水平カテーテルにおける校正のときには信号分離という別の問題が存在する。というのは、さもないと、例えば、血圧信号など患者からの実際の信号によるジャンプ信号に応じた測定ジャンプの歪みが大きすぎて、誤った補正データレコードが求められてしまう。ジャンプ信号、システムへの衝撃又はノイズの時間分解能が再現性を著しく制限し、患者からの不規則な信号の場合には、特に、大きな労力がいる。
【0009】
【特許文献1】
ドイツ連邦共和国特許出願第1982208844号明細書
【0010】
この発明の目的は、最初に述べたタイプの方法であって、外部校正信号を用いずに患者からの信号の測定値歪みを補正するための校正データを決定し、連続的にモニタし、継続した更新を行うことができ、測定信号に含まれる人為的な結果やその他の擾乱に関する誤りを受けにくい確実な補正方法を提供することである。
【0011】
この発明においては、この目的は請求項1の特徴を有する方法によって実現されている。
周波数範囲における各周波数ラインに対して測定値歪みを補正するための補正係数はこの周波数において測定された信号の統計量から求められ、そのプロセスにおいては、患者からの歪んでいない信号からの対応する統計量についての実験的にチェックされた仮定及び伝達関数の一般的な特性が用いられる。この発明によれば、外部校正信号を用いずに、患者からの信号の測定値歪みを補正するための補正データを決定し、連続的にモニタし、継続的に更新することができ、外部校正デバイスに比べて操作及びサービスが著しく容易であり擾乱を受けにくい経済的な補正方法が得られる。この補正方法における伝達関数は患者の信号から直接決定される。
【0012】
この発明による解決策の有利な特徴は、補正データが患者からの信号から決定される伝達関数から決定されることであり、また、補正データの決定において患者からの歪んでいない信号の対応する統計量に関する実験的にチェックされた仮定や伝達関数の一般的特性が使用されることである。
【0013】
この発明の解決策によれば、その周波数範囲において伝達関数を決定し、モニタし、更新し、そこから、外部校正装置を用いずに校正データレコードを得るための以下の手順が可能である。
【0014】
1.市販されているカテーテルシステムを用いて侵襲的に測定された患者からの(圧力)信号が入力信号として作用する。これに対して可能な変形は、もとの信号のかわりに、時間に対する圧力信号の導関数を使用することである。
【0015】
2.次に、心拍ごと、又は信号の長さに関してフーリエ変換に適した患者からの信号のセグメントの形で、信号のセグメント分けが行われる。ここで、最小長さから、時間に対する圧力の第2の導関数が最小値になるという境界条件のもとに、時間に対する圧力の第1の導関数の符号変化など適当な手段が探される。適した長さを決定するための別の選択肢は、最小長さから、圧力信号の近似における誤差の最小値を、高い周波数範囲では削除されたフーリエ級数から決定することである。別の選択肢である前処理ステップは、直線的トレンド調整(trend adjustment)及びセグメントの掛け算である。
【0016】
3.次に、そのセグメント長さに合ったセグメントの複素フーリエ変換と、高速フーリエ変換の場合には、ゼロパッディングを用いた変換とが行われる。
【0017】
4.伝達関数の決定は、項目3で決定される複素フーリエ係数のセグメントごとの統計に基づいている。複素フーリエ係数Siに対する信号モデルは、コンボリューションがその周波数範囲における複素数の掛け算を表しているため、
である。ここで、tiは乗数成分(伝達関数とのコンボリューション)、miは平均、viは周波数i = 0 ... 最大値に対して患者からの信号の測定値から決定されたフーリエ係数Siの変動部分である。セグメントの数n又は統計量に対するセグメントごとのフーリエ変換は、
(平均値の安定性)
及び
(統計量に対する適切な分散)
となるように決定される。ここで、mnは長さnの統計量における複素フーリエ係数の平均値、sjはj = 1 ... nとしたときの繰り返しjにおけるフーリエ係数の複素数値、ε及びδは選択した限界である。最適な閾値ε及びδはここでは実験的に決定される。
【0018】
5.従って、伝達関数を決定するために、患者からの信号に関する以下の一般的な仮定が、全体として又は部分的に使用される。すべての仮定に対して、良好な近似が成り立つ。
5.1 伝達関数は統計量の時間区間にわたって一定のままである。
5.2 患者からの信号は短い時間区間でも変動する成分を含んでおり、厳密には周期的ではない。
5.3 共振までの伝達関数は2次システムに相当している。
5.4 伝達関数は1〜2Hzの周波数範囲において1に等しい。
5.5 伝達関数は、振幅<ε(ノイズは実質的にゼロ)で、非常に高い周波数においては無視することができる。
【0019】
6.伝達関数は、すべての周波数ラインj = 1 ... 最大値、及びi = 1 ... nのセグメントに対する、
の形の統計量から決定される。パラメータp及びqに対してp>qである。
【0020】
平均値のかわりに、中央値を用いて、外れ値に関して統計量を安定化することも可能である。このタイプの統計量は、この信号モデルにおいては、複素フーリエ係数列と見なすことができる。周波数範囲にわたってより均一なサイズ分布を有する導出された(derived)係数列が、対応するフーリエ系列の複素数微分又は対数を取ることによって得られる。さらに、パラメータp及びqは実験的に決定される。特殊なケースであるp = 1及びq = 0を認め、例えば、p = 3及びq = 2に対して、患者からの信号が歪んでいない場合には、統計量
の対数の実部は、すべてのjに対して周波数の対数とともに直線的に低下する。変形は、励起信号の分布を特徴付けることによって信号モデルを改善するために、より高次のモーメントやキュムラントをさらに決定することである。
【0021】
7.実際の統計量と患者からの歪んでいない信号において期待される統計量との比、又は対数変換したあとのこれらの関数の重み付けされた差を調べて、伝達関数が決定される。ケプストラム(準同形デコンボリューション)を求めることによって、項目3のフーリエ変換と対数変換を組み合わせた変形が可能である。
【0022】
8.これらの統計量におけるノイズを低減するために、以下の方法を個々に又は組み合わせて使用することができる。
8.1 項目7による伝達関数の連続した複数の決定値を(場合によって重み付けして)平均する。
8.2 複数の周波数ラインに対するクラス形成によってスペクトル分解能を若干粗くする。
8.3 平滑化の内挿を行うか、又はメディアンフィルタ又はローパスフィルタを使用する。二次システム、又は、例えば、スプライン関数などの別の平滑(flat)関数へ直接適用して、統計量から決定されたノイズを含んだ伝達関数の偏差を最小にする。
【0023】
9.個々の周波数ラインに対する減衰係数及び共振周波数又は伝達関数の量を項目8のもとでの結果から決定することができる。
【0024】
また、共振周波数は項目6のもとに説明した形の統計量から直接決定することができる。この方法によれば、減衰を推定することができる。例えば、平均値の対数、さらによいのは、中央値の対数は周波数と相関がある。患者からの歪んだ信号の場合における対応する回帰の勾配によって、減衰に対する推定値が与えられる。より高い周波数の中央値の絶対値も、より強い減衰の場合には、より小さくなる。共振周波数及び減衰に対する複数の推定値がそれらの品質に関して重み付けして平均化されるか、多重回帰モデルを用いて組み合わせられる。対応する多重回帰モデルも非線形である。対応するモデルに対するパラメータは系統的なインビトロ(in-vitro)試験列から決定することができる。こうして、共振周波数及び減衰の非常に安定した正確な推定値を得ることができる。
【0025】
10.場合によっては、結果を改善するために、項目7から項目8に関する反復法を実行することができる。
【0026】
11.次に、位相を、二次システムの仮定を用いて伝達関数から明確に計算するか、項目6のもとで指定した形の別の統計量から決定することができる。ここで、患者からの信号が歪んでいない場合、
の形の統計量の位相は、周波数に関してS字状の関係を示している。
【0027】
12.この関数の複素逆関数は、このようにして決定された伝達関数の振幅及び位相からコンピュータによって明確に決定される。この逆関数はその周波数範囲における補正のための補正データレコードである。項目1から項目3に従って測定された圧力信号のセグメントごとのフーリエ変換が行われる。
【0028】
以下に、この発明が基づく概念が、実施の形態や図面に示された曲線を用いて、より詳しく説明される。
【0029】
図1は液体充填式伝送システムを用いて侵襲的に圧力測定を行う基本構成を示している。液体充填式伝送システムを形成するカテーテル2は、プロセスにおいて患者1の動脈系又は静脈系によって測定箇所の近くまで移動され、測定箇所では、例えば、動脈系におけるインピーダンス、圧力又は時間に対する圧力の導関数を測定して等容積収縮力(isovolumic contraction force)などが測定される。カテーテル2による患者1への影響をできる限り小さくするために、カテーテル2は寸法ができる限り小さくされており、弾性材料から形成されていて、ホース状に形成されていることが好ましい。液体が充填されたカテーテル2の先端には開口部が設けられている。この開口部を介して、例えば、圧力パルスなどが取り出され、カテーテル2及び液体が充填されたラインを介して、圧力トランスデューサ3まで伝えられる。
【0030】
圧力トランスデューサ3は圧力パルスの関数として電気信号を発生し、それに応じて、この信号は表示されたり評価されたりする。この方法は、原理的には長い間知られていた。この伝送システムの伝達関数の補正は、類似の電気回路又は対応する数値アルゴリズムを用いて、共振周波数及び減衰係数を決定したあとに行われる。
【0031】
上述した方法を使用するときに生じる40%にも及ぶ範囲の測定歪みを効果的に補正するために、この発明の方法においては、アナログ/デジタル変換器4が圧力トランスデューサ3と信号解析及び処理装置5との間に配置されている。変換器4は圧力トランスデューサ3のアナログ信号をデジタル信号へ変換し、このデジタル信号は信号解析及び処理装置5の入力へ加えられる。信号解析及び処理装置5において、測定データはデジタルフーリエ解析に基づく補正ファクタがロードされ、このようにして補正された測定値は出力又は評価装置7へ送られる。
【0032】
圧力信号の測定値の歪みを補正するための校正データレコード6を決定し、モニタし、更新するためのデバイスが、アナログ/デジタル変換器4の出力へ接続されているか、又は信号解析及び処理装置5へ接続されていて、カテーテルラインシステムのパラメータ及び校正ファクタを決定し、モニタし、更新する。補正されたフーリエ係数がロードされたデジタル信号は、信号解析及び処理装置5からディスプレイ又は評価装置6へ送られ、そこでモニタシステム及びプリンタの両方に表示される。モニタの規格に応じて、信号は、まず、デジタル/アナログ変換器へ供給され、そのあと出力されるか、又はデジタル信号を処理することができるモニタへ直接中継される。場合によって、信号はディスプレイに適したフォーマットに処理される必要がある。
【0033】
別の可能性は、データを、それを蓄積して評価するコンピュータへ転送することである。この場合には、データはデジタル/アナログ変換器では処理されず、補正から直接送られる。
【0034】
オンラインの補正を実行せずに、あとで評価又は補正を行うために、データを蓄積することもできる。このための要件は、システム特有のデータや測定条件に関する情報が利用できて、そのあと、補正データレコードを適切に選択することができることである。このために、データは圧力トランスデューサ3のあとで直接取り出されて、例えば、CDやディスケットなどの適当な蓄積媒体の上に記録される。
【0035】
図2は統計量
の対数をとった多数の測定値(=統計量の複素対数の実部)のスペクトルを重ね合わせた曲線を示している。下側のグラフは先端圧力センサに対する曲線を示しており、上側のグラフは液体充填式システムによって歪んだ測定値を示している。先端圧力の曲線においては、指数関数的な低下を明瞭に見てとることができ、これは周波数軸を対数化すると直線的な低下になる。液体充填式伝送システムによって歪んだ測定値の対応する統計量は、この指数関数的な低下からは著しく外れている。先端圧力の曲線のスペクトルと液体充填式伝送システムによって歪んだ測定値のスペクトルとの間のこの差は伝達関数の結果として生じる。
【0036】
従って、下側の図に示されている曲線が上側の図に示されている曲線から適切に決定できるときには、この伝達関数を決定することができる。不具合のタイプに関係なく、測定圧力と個々の測定との間の指数関数的なフィッティングは、患者の歪んでいないすべての圧力測定に対して非常に良好である(r>0.9)。長いデータレコードに対して実験的にチェックされたこの関係が、液体充填式の伝送システムによって歪んだ測定値に対する対応する統計量から、先端圧力曲線に対する統計量の形を推定するための基礎となる。
【0037】
図3の上側は、先端圧力曲線に対する統計量
と、液体充填式の伝送システムによって歪んだ測定値に対する対応する統計量の対数をとった測定値のスペクトル(=統計量の複素対数の実部)を示している。周波数軸は対数(自然対数)化されている。伝達関数の範囲における不一致は明らかである。
【0038】
先端圧力曲線及び液体充填式伝送システムによって歪んだ測定値に対する対応する回帰曲線が下側の図に示されている。これは、歪んだ測定値の曲線から先端圧力曲線に対する回帰曲線を決定するときの問題点を明らかにしている。二つの曲線の勾配及びずれが一致しないことから、単純な回帰ではこれを行えない。この問題を処理するための対応する試みが図4に関連して説明される。
【0039】
回帰を差し引いたあとの残差の解析において、共振周波数の範囲における伝達関数によって引き起こされるピークが、歪んだ測定値の残差においては明らかに見られる。先端圧力曲線の残差は、高い周波数範囲においては、ゼロのあたりで振動しているが、低い周波数範囲においては、パルス信号の大きな基本振動及び第1の上側振動によって引き起こされる著しいピークを示しており、信号モデルを設定するときには考慮する必要がある。
【0040】
図4は、先端圧力曲線(tip)に対する統計量
と、液体充填式の伝送システムによって歪んだ測定値(mes)に対する対応する統計量の対数をとった測定値のスペクトル(=統計量の複素対数の実部)を示している。この目的のために、周波数軸は対数(自然対数)化されている。
【0041】
この先端圧力曲線(tip reg 1及びtip reg 2)及び液体充填式伝送システムによって歪んだ測定値(mes reg 1及びmes reg 2)に対する対応する回帰は、図の中央で交差する上側の二つのラインである。低い周波数における伝達関数は1に等しく(上述した特徴区分の5.4項)、ある信号閾値に達していないときには無視できるという仮定を用いると、2点回帰のタイプを形成することもできる(図4における下側の二つのラインを参照のこと)。
【0042】
先端圧力の参照測定を行わないとき(通常の測定の場合)には、歪んだ測定値を表す二つの曲線及び歪んだ測定値の回帰が利用可能である。先端圧力測定を表す曲線は推定されなければならない。これは、歪んだ測定の二つの回帰を重み付けして平均することによって妥当な形で行われる。
【0043】
この推定は共振周波数のラフな推定によって反復的に改善することができる。このステップも実験的に広くチェックされる。
【0044】
図5及び図6は、共振が外部校正によって独立に決定された二つの異なる液体充填式システムに対して、統計量
の対数をとった測定値のスペクトル(=統計量の複素対数の実部)を種々の方法によって平滑化した残差を示している。
【0045】
図5及び図6において、
“org” もとの測定、すなわち、処理されていない測定曲線を表している。
“tp” ローパスフィルタによって平滑化された曲線。
“mvg” 平均値フィルタによって平滑化された曲線を表している。
“cub” 三次多項式による内挿によって平滑化された曲線。
“pol” 一般的な多項式でフィルタリングされた曲線を表している。
【0046】
ローパスフィルタによって実行された平滑化はこのプロセスにおいて最も安定していることがわかった。実験的に決定された共振はリングでマーキングされている。周波数軸は対数化されているが、対応する対数化されていない周波数が記載されている。
【0047】
結果は図7に示されている。この図には図5及び図6に対応する曲線が記載されており、対応した表記がされている。しかし、励起信号(先端信号、参照信号)のスペクトルにおいて相対的なギャップがある測定においては、共振周波数の領域において乱れている。例えば、二次システムへのフィッティングなどの内挿による改善方法によってこの問題は克服される。
【0048】
可能性のある別の擾乱問題は、パルス信号の大きな基本振動及び第1の上側振動によって引き起こされる元の信号の低い周波数領域におけるピークである。この問題は以下のようにして克服できる。
1.例えば、3Hzのもっともらしい共振周波数に対する下限を設定することによって。
2.低い周波数範囲(1(2)Hzまで)における残差をゼロに設定することによって。
3.低い周波数範囲においては伝達関数が1であるという仮定をして、推定する参照の回帰のスタートを測定値によって置き換えるようにモデルを改善することによって。不変性を確保するために、二つの曲線を重みを付けて平均することが好ましい。1Hzにおける測定値の重みは1であり、比較的急速に0まで低下する。1Hzにおける回帰の重みは0であり、比較的急速に1まで上昇する。
【0049】
たとえこれらの改善(二次システムへのフィッティング及び信号モデルの改善)をしなくても、この方法によれば、減衰係数及び共振周波数を、患者(症例)、圧力のタイプ及び伝送システムとは独立に正確に決定することができる。これは、外部校正によって独立に決定された伝達関数との系統的な比較研究によって証明することができる。
【図面の簡単な説明】
【図1】 測定値の歪みを補正するための補正データを決定し、モニタし、更新して、液体充填式伝送システムを校正するために、伝達関数とのコンボリューションをとって、患者からの信号の液体充填式システムの伝達関数を情報のない状態で決定するための装置の機能図を示している。
【図2】 先端圧力センサを用いた測定及び液体充填式システムによって測定値が歪んでいる測定における統計量の対数をとった多数の測定結果のスペクトルを重ね合わせて示したグラフである。
【図3】 先端圧力センサを用いた測定及び液体充填式システムによって測定値が歪んでいる測定における統計量の対数をとった測定結果のスペクトルを、それぞれ対数化した周波数軸で示したグラフである。
【図4】 先端圧力センサを用いた測定及び液体充填式システムによって測定値が歪んでいる測定における統計量の対数をとった測定結果のスペクトルを、それぞれ対数化した周波数軸で示したグラフであり、先端圧力センサを用いた測定と液体充填式システムによって歪んだ測定値における曲線に対して回帰直線が描かれている。
【図5】 液体充填式システムに対して外部校正を用いて独立に決定された共振の対数をとった測定結果の残差を種々の方法で平滑化したグラフを示している。
【図6】 別の液体充填式システムに対して外部校正を用いて独立に決定された共振の対数をとった測定結果の残差を種々の方法で平滑化したグラフを示している。
【図7】 共振が実験的に決定される図5及び図6におけると同様のグラフを示している。
【符号の説明】
3 圧力トランスデューサ
5 信号解析及び処理装置
6 校正データレコード
7 評価及び出力装置
Claims (30)
- 心臓病学、集中治療医学及び麻酔術における侵襲的な圧力測定を行う際の液体充填式伝送システムによって引き起こされた測定値の歪みを補正するために、補正データを決定し、モニタし、更新して、液体充填式伝送システムを校正するための方法であって、当該補正データの決定、モニタ、及び更新は、患者からの、あとで評価又は補正を行なうために蓄積された測定信号を決定した後の時点においてのみ専ら実行されるものであり、
前記液体充填伝送システムを使用して測定された患者からの前記測定信号のスペクトル統計量を生成する段階と、
補正データを、周波数範囲における各周波数に対して、患者からの前記測定信号の前記スペクトル統計量から対数変換を用いて決定する段階と、
を有する方法。 - 補正データが患者からの前記測定信号から得られる伝達関数から決定される請求項1に記載の方法。
- 補正データを決定する段階において、患者からの歪んだ信号の対応する統計量に関する実験的にチェックされる仮定及び伝達関数の一般的な特性が使用される請求項1に記載の方法。
- a)患者からの前記測定信号をセグメントに分ける段階と、
b)セグメントの長さに合ったセグメントの複素フーリエ変換を行う段階と、
c)複素フーリエ係数のセグメントごとの統計量を決定する段階と、
d)決定された複素フーリエ係数のセグメントごとの実際の統計量と、患者からの歪んでいない信号において期待される複素フーリエ係数の統計量との比から、伝達関数を決定する段階と、
を有する請求項2又は請求項3記載の方法。 - セグメントの長さに合ったセグメントの複素フーリエ変換がゼロパッディングを用いた高速フーリエ変換で行われる請求項4に記載の方法。
- a)患者からの前記測定信号をセグメントに分ける段階と、
b)セグメントの長さに合ったセグメントの複素フーリエ変換を行う段階と、
c)複素フーリエ係数のセグメントごとの統計量を決定する段階と、
d)セグメントごとの実際の統計量と、患者からの歪んでいない信号の場合に期待される対応する統計量の推定値とを対数変換する段階と、
e)統計量の対数変換における差から伝達関数を決定する段階と、
を有する請求項2又は請求項3記載の方法。 - a)患者からの前記測定信号をセグメントに分ける段階と、
b)セグメントの長さに合ったセグメントの複素フーリエ変換を行う段階と、
c)複素フーリエ係数のセグメントごとの統計量を得る段階と、
d)係数としての統計量の値から形成されたフーリエ級数を微分する段階と、
e)微分された級数の係数から共振及び減衰を決定する段階と、
を有する請求項2又は請求項3記載の方法。 - 患者からの前記測定信号をセグメントに分けてケプストラムを決定する請求項2又は請求項3に記載の方法。
- 侵襲的に動脈内及び静脈内の少なくとも一方で圧力測定を行うために、カテーテルシステムによって侵襲的に測定される患者からの圧力信号のセグメント分けが心拍ごとに行われる請求項1〜請求項8のいずれか一項に記載の方法。
- セグメント分けが患者セグメントからの信号において行われ、患者セグメントからの信号がそれらの長さに関してフーリエ変換に適している請求項1〜請求項9のいずれか一項に記載の方法。
- フーリエ変換に適した患者セグメントからの信号の長さが、予め設定可能な最小長さから、時間に対する患者からの前記測定信号の第2の導出関数の最小値の境界条件において、時間に対する患者からの測定信号の第1の導出関数の符号が変化するように決定される請求項10に記載の方法。
- フーリエ変換に適した患者セグメントからの信号の長さが、予め設定可能な最小長さから、高い周波数範囲において削除されたフーリエ級数による患者からの前記測定信号の近似において誤りが最小になるように決定される請求項10に記載の方法。
- 患者からの前記測定信号又は患者からの圧力信号のセグメント分けのトレンド調整が行われる請求項1〜請求項12のいずれか一項に記載の方法。
- 患者からのセグメント化された信号又は患者からの圧力信号の乗算が行われる請求項1〜請求項13のいずれか一項に記載の方法。
- 複素フーリエ係数siに対する信号モデルが、
であり、ここで、tiは大きくとも乗算コンポーネント(伝達関数とのコンボリューション)であり、miは平均値であり、viは周波数i = 0 ... に対する測定から決定されたフーリエ係数siの変動部分であり、セグメントの数又はセグメントごとのフーリエ変換の数nは、適切な平均値安定性である
という条件と、適切な分散である
という条件によって、j = 1 ... nとしたときの繰り返しjにおける長さn及びフーリエ係数の複素数であるsjの統計量及び選択された境界ε及びδにおいて、フーリエ係数の平均値であるmnを有する統計量が与えられるように決定される請求項1〜請求項14のいずれか一項に記載の方法。 - 最適な閾値ε及びδが実験的に決定される請求項15に記載の方法。
- 推定値として中央値が使用され、この中央値が、対応するフーリエ級数及びそこから導出された係数列の対数の複素微分によって得られ、周波数範囲にわたって均一なサイズ分布を有する請求項15又は請求項16に記載の方法。
- 他の高次のモーメント又はキュムラントが決定される請求項17又は請求項18に記載の方法。
- 連続した複数の伝達関数を平均化する段階を有する請求項1〜請求項19のいずれか一項に記載の方法。
- 統計モデルにおいて、連続した複数の伝達関数を重み付けして平均化するか組み合わせる段階を有する請求項20に記載の方法。
- 空間状態モデルとして連続的に決定された伝達関数の形式的モデル化を行う段階を有する請求項1〜請求項21のいずれか一項に記載の方法。
- 複数の周波数ラインに対するクラス形成によってスペクトル分解能を粗くする段階を有する請求項1〜請求項22のいずれか一項に記載の方法。
- 決定された伝達関数を平滑化する内挿及びローパスフィルタを用いた伝達関数のフィルタリングの少なくとも一方が行われる請求項1〜請求項23のいずれか一項に記載の方法。
- 統計量から決定されたノイズを含んだ伝達関数の偏差が最小になるように、二次システムへ直接適合させる段階を有する請求項1〜請求項24のいずれか一項に記載の方法。
- 実際の統計量と患者からの歪んでいない信号における統計量との比から、又は決定された伝達関数を対数変換したあとのこれらの関数の差から、期待される統計量におけるノイズを低減したあと、減衰係数及び共振周波数又は伝達関数の値が個々の周波数ラインに対して決定される請求項1〜請求項25のいずれか一項に記載の方法。
- 結果を改善するために反復法が用いられる請求項20〜請求項26のいずれか一項に記載の方法。
- 伝達関数から位相が計算される請求項1〜請求項27のいずれか一項に記載の方法。
- 請求項17及び請求項18に記載された形の別の統計量から位相が決定される請求項1〜請求項27のいずれか一項に記載の方法。
- 関数の複素逆関数がこのようにして決定された伝達関数の振幅及び位相からコンピュータによって決定され、この複素逆関数が周波数範囲において補正のための補正データレコードとして作用し、患者から測定された信号のセグメントごとのフーリエ変換が行われる請求項1〜請求項29のいずれか一項に記載の方法。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE10049734A DE10049734A1 (de) | 2000-09-28 | 2000-09-28 | Verfahren zur Ermittlung, Überwachung und Aktualisierung von Korrekturdaten für flüssigkeitsgefüllte Übertragungssysteme |
DE10049734.9 | 2000-09-28 | ||
PCT/DE2001/003719 WO2002026127A1 (de) | 2000-09-28 | 2001-09-25 | Verfahren zur ermittlung, überwachung und aktualisierung von korrekturdaten für flüssigkeitsgefüllte übertragungssysteme |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2004509687A JP2004509687A (ja) | 2004-04-02 |
JP5161413B2 true JP5161413B2 (ja) | 2013-03-13 |
Family
ID=7659016
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2002529960A Expired - Fee Related JP5161413B2 (ja) | 2000-09-28 | 2001-09-25 | 液体充填式伝送システムの補正データを決定、モニタ及び更新するための方法 |
Country Status (12)
Country | Link |
---|---|
US (1) | US7264593B2 (ja) |
EP (1) | EP1320323B1 (ja) |
JP (1) | JP5161413B2 (ja) |
AT (1) | ATE345085T1 (ja) |
AU (2) | AU2344002A (ja) |
CA (1) | CA2423959C (ja) |
CY (1) | CY1106005T1 (ja) |
DE (2) | DE10049734A1 (ja) |
DK (1) | DK1320323T3 (ja) |
ES (1) | ES2275751T3 (ja) |
PT (1) | PT1320323E (ja) |
WO (1) | WO2002026127A1 (ja) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102004009871A1 (de) * | 2004-02-26 | 2005-09-15 | Deutsches Herzzentrum Berlin | Verfahren zur Bestimmung hämodynamischer Parameter |
WO2012013662A1 (en) | 2010-07-26 | 2012-02-02 | Steerable Instruments Bvba | Capillary tube assembly |
US20130123661A1 (en) * | 2010-07-26 | 2013-05-16 | Steerable Instruments Bvba | Endoscopic pressure detection assembly |
US10537274B2 (en) | 2011-09-22 | 2020-01-21 | The U.S. Govt As Represented By The Dept. Of Veterans Affairs | Standardized measurement of physiological pressures using an air-charged catheter apparatus |
FR2984720B1 (fr) * | 2011-12-22 | 2014-03-07 | Univ Grenoble 1 | Procede et dispositif de surveillance de la mesure de la pression arterielle par catheterisme arteriel d'un patient |
DE102017110342A1 (de) * | 2017-05-12 | 2018-11-15 | Prüftechnik Dieter Busch AG | Effektivwertbestimmung einer Maschinenschwingungsgröße |
WO2021062086A1 (en) * | 2019-09-27 | 2021-04-01 | Board Of Regents Of The University Of Texas System | A system and method for cardiac pressure measurement |
CN116049341B (zh) * | 2023-03-08 | 2023-08-15 | 北京七兆科技有限公司 | 一种水文数据标准化方法、装置、设备及存储介质 |
Family Cites Families (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US653749A (en) * | 1900-04-19 | 1900-07-17 | William R Mackay | Toy rapid-fire cannon. |
US1403719A (en) * | 1921-03-05 | 1922-01-17 | Keystone Die And Mfg Company | Toy gun |
US1743576A (en) * | 1927-07-14 | 1930-01-14 | Smith Robert Bigham | Pneumatically-actuated machine gun |
US1789206A (en) * | 1930-02-19 | 1931-01-13 | Walter V Whatley | Machine-gun toy |
US1927424A (en) * | 1931-12-21 | 1933-09-19 | Trubenbach Walter | Toy rapid-fire gun |
US3272397A (en) * | 1963-06-20 | 1966-09-13 | Sherman G Bean | Feeder of non-flowing material |
US3348531A (en) * | 1964-03-06 | 1967-10-24 | Sekiden Kagaku Kogyo Co Ltd | Toy gun having a divided magazine |
US3263664A (en) * | 1964-12-29 | 1966-08-02 | Martin G Bauer | Bb gun loader |
US3695246A (en) * | 1971-06-10 | 1972-10-03 | Us Navy | Pneumatic machine gun with photo cell interrupted circuit |
US4025071A (en) * | 1976-03-05 | 1977-05-24 | Hodges Kenneth M | Tennis ball server and court installation |
US4232373A (en) * | 1978-04-14 | 1980-11-04 | Regents For Education Of The State Of Rhode Island | Compensation of fluidic transducers |
US4325382A (en) * | 1980-05-15 | 1982-04-20 | Memorial Hospital For Cancer And Allied Diseases | Process and apparatus for the real time adaptive filtering of catheter pressure measurements |
US4405060A (en) * | 1981-07-20 | 1983-09-20 | American Hospital Supply Corporation | Tablet dispensing device |
JPH0448162Y2 (ja) * | 1986-12-26 | 1992-11-13 | ||
US4980552A (en) * | 1989-06-20 | 1990-12-25 | The Regents Of The University Of California | High resolution PET scanner using rotating ring array of enlarged detectors having successively offset collimation apertures |
US5166457A (en) * | 1992-01-22 | 1992-11-24 | Lorenzetti James A | Ammunition magazine for paint ball gun |
US5505188A (en) * | 1994-03-17 | 1996-04-09 | Williams; Robert A. | Paint ball gun |
US5542406A (en) * | 1994-08-22 | 1996-08-06 | Oneto; Michael A. | Retractable bolt assembly for compressed gas powered gun |
US5794606A (en) * | 1996-05-28 | 1998-08-18 | Deak; Bernard A. | Ram feed ammo box |
US5809983A (en) * | 1996-11-29 | 1998-09-22 | Stoneking; Scot E. | Lighting loader system |
US5947100A (en) * | 1997-04-30 | 1999-09-07 | Anderson; Joel A. | Paint ball gun agitator sound trigger and duration control |
US5827195A (en) * | 1997-05-09 | 1998-10-27 | Cambridge Heart, Inc. | Electrocardiogram noise reduction using multi-dimensional filtering |
US5839422A (en) * | 1997-05-23 | 1998-11-24 | Ferris; Shell M. | Automatic feeder for projectile gun using compressed gas |
CA2311535C (en) * | 1997-11-21 | 2007-11-27 | Deutsches Herzzentrum Berlin | Correction of blood pressure measurements in invasive liquid-filled systems |
DE19820844B4 (de) * | 1997-11-21 | 2006-01-26 | Deutsches Herzzentrum Berlin | Verfahren und Vorrichtung zur Korrektur von Meßwertverfälschungen bei invasiven Druckmessungen mit flüssigkeitsgefüllten Systemen |
US6462342B1 (en) * | 2000-06-22 | 2002-10-08 | Ge Medical Systems Global Technology Co. Llc | Method and system for pet image reconstruction |
USD459767S1 (en) * | 2001-04-19 | 2002-07-02 | Ennis Rushton | Paintball hopper for a paintball launcher |
-
2000
- 2000-09-28 DE DE10049734A patent/DE10049734A1/de not_active Withdrawn
-
2001
- 2001-09-25 PT PT01985667T patent/PT1320323E/pt unknown
- 2001-09-25 EP EP01985667A patent/EP1320323B1/de not_active Expired - Lifetime
- 2001-09-25 AU AU2344002A patent/AU2344002A/xx active Pending
- 2001-09-25 ES ES01985667T patent/ES2275751T3/es not_active Expired - Lifetime
- 2001-09-25 JP JP2002529960A patent/JP5161413B2/ja not_active Expired - Fee Related
- 2001-09-25 CA CA2423959A patent/CA2423959C/en not_active Expired - Fee Related
- 2001-09-25 AT AT01985667T patent/ATE345085T1/de active
- 2001-09-25 DK DK01985667T patent/DK1320323T3/da active
- 2001-09-25 DE DE50111464T patent/DE50111464D1/de not_active Expired - Lifetime
- 2001-09-25 WO PCT/DE2001/003719 patent/WO2002026127A1/de active IP Right Grant
- 2001-09-25 US US10/381,942 patent/US7264593B2/en not_active Expired - Fee Related
- 2001-09-25 AU AU2002223440A patent/AU2002223440B2/en not_active Ceased
-
2007
- 2007-02-12 CY CY20071100180T patent/CY1106005T1/el unknown
Also Published As
Publication number | Publication date |
---|---|
EP1320323B1 (de) | 2006-11-15 |
CA2423959C (en) | 2011-09-13 |
ES2275751T3 (es) | 2007-06-16 |
DE50111464D1 (de) | 2006-12-28 |
CA2423959A1 (en) | 2003-03-28 |
CY1106005T1 (el) | 2011-04-06 |
AU2002223440B2 (en) | 2006-06-15 |
DK1320323T3 (da) | 2007-03-12 |
JP2004509687A (ja) | 2004-04-02 |
US20040024294A1 (en) | 2004-02-05 |
EP1320323A1 (de) | 2003-06-25 |
DE10049734A1 (de) | 2002-04-11 |
ATE345085T1 (de) | 2006-12-15 |
PT1320323E (pt) | 2007-02-28 |
US7264593B2 (en) | 2007-09-04 |
AU2344002A (en) | 2002-04-08 |
WO2002026127A1 (de) | 2002-04-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU2018235369B2 (en) | Central aortic blood pressure and waveform calibration method | |
US7887490B2 (en) | Method and device for removing respiratory artefacts from measured blood pressure data | |
AU729956B2 (en) | Improved method and apparatus for the measurement of cardiac output | |
US4418700A (en) | Method and apparatus for measurement of heart-related parameters | |
KR100877753B1 (ko) | 혈역학 파라미터 측정 장치 및 방법 | |
US6761695B2 (en) | Method and apparatus for non-invasive measurement of changes in intracranial pressure | |
JP4939393B2 (ja) | 血流力学的パラメータを決定するための方法 | |
WO2011138683A2 (en) | Brachial cuff | |
US6290652B1 (en) | Correction of blood pressure measurements in invasive liquid-filled systems | |
AU2018227095B2 (en) | Non-invasive brachial blood pressure measurement | |
JP5161413B2 (ja) | 液体充填式伝送システムの補正データを決定、モニタ及び更新するための方法 | |
US20060064021A1 (en) | Detection and correction of catheter line distortion in blood pressure measurements | |
KR102373864B1 (ko) | 압박대 없이 장시간 동안 연속적으로 혈압을 측정하는 장치 및 그 방법 | |
US20240206744A1 (en) | Setting device and method for a blood pressure monitoring system | |
WO2023186745A1 (en) | Apparatus for determining an indicator representative for a physiological parameter | |
DeBusk | Signal restoration in blood pressure waveforms measured with catheter-transducer systems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20080925 |
|
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A711 Effective date: 20090401 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20090401 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20110225 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20110906 |
|
A602 | Written permission of extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A602 Effective date: 20111213 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20120305 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20120327 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20120627 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20120807 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20121026 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20121127 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20121214 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20151221 Year of fee payment: 3 |
|
LAPS | Cancellation because of no payment of annual fees |