JP2000511449A - 超音波装置および対照試薬を較正する方法およびそのための装置 - Google Patents
超音波装置および対照試薬を較正する方法およびそのための装置Info
- Publication number
- JP2000511449A JP2000511449A JP10500624A JP50062498A JP2000511449A JP 2000511449 A JP2000511449 A JP 2000511449A JP 10500624 A JP10500624 A JP 10500624A JP 50062498 A JP50062498 A JP 50062498A JP 2000511449 A JP2000511449 A JP 2000511449A
- Authority
- JP
- Japan
- Prior art keywords
- suspension
- backscatter
- microsphere
- microsphere suspension
- microspheres
- 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.)
- Withdrawn
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 239000003153 chemical reaction reagent Substances 0.000 title abstract description 26
- 239000000725 suspension Substances 0.000 claims abstract description 65
- 239000004005 microsphere Substances 0.000 claims abstract description 63
- 239000004793 Polystyrene Substances 0.000 claims abstract description 27
- 229920002223 polystyrene Polymers 0.000 claims abstract description 27
- 238000002604 ultrasonography Methods 0.000 claims abstract description 22
- 239000007788 liquid Substances 0.000 claims abstract description 15
- 239000007787 solid Substances 0.000 claims description 16
- 238000003756 stirring Methods 0.000 claims description 9
- 238000002156 mixing Methods 0.000 claims description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 6
- 239000000203 mixture Substances 0.000 claims description 4
- 238000005086 pumping Methods 0.000 claims description 4
- 239000011324 bead Substances 0.000 claims description 2
- 238000005259 measurement Methods 0.000 abstract description 34
- 230000000694 effects Effects 0.000 description 12
- 239000002245 particle Substances 0.000 description 11
- 238000012546 transfer Methods 0.000 description 9
- 229910000831 Steel Inorganic materials 0.000 description 8
- 239000010959 steel Substances 0.000 description 8
- 230000001419 dependent effect Effects 0.000 description 7
- 238000001228 spectrum Methods 0.000 description 7
- 238000002474 experimental method Methods 0.000 description 6
- 239000000463 material Substances 0.000 description 6
- 230000009467 reduction Effects 0.000 description 6
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 4
- 239000012530 fluid Substances 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 239000011780 sodium chloride Substances 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000012935 Averaging Methods 0.000 description 3
- 239000002033 PVDF binder Substances 0.000 description 3
- 238000003957 acoustic microscopy Methods 0.000 description 3
- 238000012512 characterization method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 230000000704 physical effect Effects 0.000 description 3
- 229920002981 polyvinylidene fluoride Polymers 0.000 description 3
- 238000004513 sizing Methods 0.000 description 3
- 229910001220 stainless steel Inorganic materials 0.000 description 3
- 239000010935 stainless steel Substances 0.000 description 3
- PPBRXRYQALVLMV-UHFFFAOYSA-N Styrene Chemical compound C=CC1=CC=CC=C1 PPBRXRYQALVLMV-UHFFFAOYSA-N 0.000 description 2
- 239000000872 buffer Substances 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 239000010408 film Substances 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 239000006194 liquid suspension Substances 0.000 description 2
- 210000004072 lung Anatomy 0.000 description 2
- 239000002609 medium Substances 0.000 description 2
- 230000008520 organization Effects 0.000 description 2
- 230000010363 phase shift Effects 0.000 description 2
- 239000004033 plastic Substances 0.000 description 2
- 229920003023 plastic Polymers 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 239000004094 surface-active agent Substances 0.000 description 2
- 108010056388 Albunex Proteins 0.000 description 1
- 229920002799 BoPET Polymers 0.000 description 1
- 229920004142 LEXAN™ Polymers 0.000 description 1
- 239000005041 Mylar™ Substances 0.000 description 1
- 239000004698 Polyethylene Substances 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 239000012736 aqueous medium Substances 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000008033 biological extinction Effects 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- HGAZMNJKRQFZKS-UHFFFAOYSA-N chloroethene;ethenyl acetate Chemical compound ClC=C.CC(=O)OC=C HGAZMNJKRQFZKS-UHFFFAOYSA-N 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000002405 diagnostic procedure Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000010790 dilution Methods 0.000 description 1
- 239000012895 dilution Substances 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000010419 fine particle Substances 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000002608 intravascular ultrasound Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 230000013011 mating Effects 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000000386 microscopy Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- -1 polyethylene Polymers 0.000 description 1
- 229920000573 polyethylene Polymers 0.000 description 1
- 229920003053 polystyrene-divinylbenzene Polymers 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 229920006298 saran Polymers 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 239000010409 thin film Substances 0.000 description 1
- 230000002792 vascular Effects 0.000 description 1
- 238000011179 visual inspection Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/02—Analysing fluids
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/58—Testing, adjusting or calibrating the diagnostic device
- A61B8/587—Calibration phantoms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
- G01N29/0654—Imaging
- G01N29/0681—Imaging by acoustic microscopy, e.g. scanning acoustic microscopy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/22—Details, e.g. general constructional or apparatus details
- G01N29/30—Arrangements for calibrating or comparing, e.g. with standard objects
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/01—Indexing codes associated with the measuring variable
- G01N2291/015—Attenuation, scattering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/024—Mixtures
- G01N2291/02416—Solids in liquids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/04—Wave modes and trajectories
- G01N2291/044—Internal reflections (echoes), e.g. on walls or defects
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/10—Number of transducers
- G01N2291/101—Number of transducers one transducer
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Acoustics & Sound (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Molecular Biology (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Sampling And Sample Adjustment (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
Abstract
(57)【要約】
本発明は、広帯域体外較正ツールとしてポリスチレン微小球を用いた超音波装置および超音波対照試薬の較正方法に関する。本方法によれば、音響アパーチャを有する容器内にポリスチレン微小球の懸濁物を液体環境で準備する。微小球懸濁物からの後方散乱の測定中、微小球懸濁物を連続的に攪拌してから、超音波装置または対照試薬を所定の後方散乱設定値に較正または調整する。
Description
【発明の詳細な説明】
超音波装置および対照試薬を較正する方法およびそのための装置 発明の分野
本発明は、診断用の超音波装置および超音波対照試薬の較正方法、特に組織ま
たは対照試薬の後方散乱および減衰の絶対測定を行うための較正方法に関する。発明の背景
血管内超音波音響顕微鏡検査および超音波対照試薬での第二高周波撮像の出現
により、近年では医療用に高周波超音波を使用することが増加している。高周波
数では、新しい形態構造を検査できると共に、超音波組織特徴付けの新しい領域
を開くことができる。従来では、組織特徴付けは、様々な種類の組織または対照
試薬の後方散乱および減衰の絶対測定を行うための標準較正方法がないことによ
って阻害されていた。
診断用超音波装置の較正の現在の標準は、研磨ステンレス鋼板からの反射を使
用している。この方法は、プレートの表面の整合に敏感であることと、反射が組
織から戻るものよりもはるかに大きいことから、不都合である。組織からの超音
波散乱は、水様媒体内に懸濁した散乱体群としてモデル化されることが多い。較
正標準すなわち組織を模倣する「ファントム」を使用することは、組織特徴づけ
の分野で特に注目されてきた。一般的にこれらのファントムは、米国特許第4,
843,866号および第4,277,367号に記載されているようなアガロ
ースまたはゼラチン様媒体内に保持された散乱体群からなり、これら米国特許の
内容はすべて本説明の中に含まれる。5MHzより低い周波数範囲の超音波用の
ファントムを検査するためには相当量の作業が必要である。さらに、組織内の散
乱体の幾何学的形状および物理的パラメータを適当なデータ整理方法を用いて算
出することができる。しかし、ゲルを基にしたファントムを使用することは、安
定した結果を生じるために多くの空間位置平均を得なければならないと共に、多
大な専門知識の構築および維持を必要とする点で、不都合である。
しかし、音響顕微鏡検査および血管内撮像のさらに高周波の分野では、存在す
るファントムが有効で再現性がある較正方法になるかどうかが明らかではない。
従って当該分野では、組織によって戻される反射をほぼ模倣し、構築および維持
が容易であり、較正ツールの整合に関する問題をなくす診断用超音波装置および
それと共に使用される超音波対照試薬の広帯域体外較正方法が必要とされている
。発明の概要
本発明は、広帯域体外較正ツールとしてポリスチレン微小球を用いた診断用超
音波装置および超音波対照試薬の較正方法に関するものである。本方法によれば
、音響アパーチャを有する容器内にポリスチレン微小球の懸濁物を液体環境で準
備する。微小球懸濁物からの後方散乱の測定中、微小球懸濁物を連続的に攪拌し
てから、超音波装置または対照試薬を所定の後方散乱設定に較正または調整する
。
本発明の1つの重要な態様によれば、本方法は、一定温度の水浴内の音響アパ
ーチャを有する容器内に液体環境で微小球の懸濁物を準備する段階と、微小球懸
濁物からの後方散乱の測定中、微小球懸濁物を連続的に攪拌する段階と、超音波
装置を所定の後方散乱設定に調整する段階を含む。後方散乱の測定値に従って、
平均後方散乱、後方散乱係数、後方散乱伝達関数または平均後方散乱伝達関数等
の様々な関数を上述の調整段階の前に計算してもよい。図面の簡単な説明
図1は、ポリスチレン微小球懸濁物の保持および攪拌に使用される試料室の概
略図である。
図2は、ポリスチレン微小球懸濁物の見掛け後方散乱伝達関数の測定に使用さ
れる実験機構を示している。
図3は、ポリスチレン微小球懸濁物の信号損失(減衰)の測定に使用される実
験機構を示している。
図4は、平均直径が49.5μmで標準偏差が1.0μmのポリスチレン微小
球の一般的な分布を示してる。
図5の上パネルは、平均直径が49.5μmのポリスチレン微小球の典型的な
見掛け後方散乱伝達関数の有効帯域を示している。中間パネルは見掛け減衰の対
応測定値を示している。下パネルは結果的に得られる後方散乱係数である。すべ
ての誤差バールは測定値の標準偏差を表している。
図6は、単一寸法の微小球集団と図3に示されている微小球寸法分布の予測後
方散乱係数を示している。
図7は、後方散乱係数の理論的予測値と実験結果との比較を示している。発明の詳細な説明
本発明によれば、液体環境で微小球懸濁物を使用し、懸濁物の後方散乱の測定
中に懸濁物を連続的に攪拌することによって、診断用超音波装置および超音波対
照試薬を較正する有効な方法が提供されることがわかった。懸濁物は、手動およ
び機械式手段を含む様々な攪拌手段で、例えばかきまぜ、混合、ポンピングおよ
びふりまぜ等によって連続的に攪拌することができる。好ましくは、液体環境で
微小球懸濁物を収容している容器と連通した手段を設けると共に、微小球懸濁物
の後方散乱の測定中に任意位置で懸濁状態に維持できるように微小球の十分な混
合および攪拌を行う、液体内へのポンピング手段を設けることによって、懸濁物
が連続的に攪拌される。較正を行うのに適当な装置が図1に示されている。図1
に示されている試料室内において、微小球懸濁物を収容している容器内のポート
として1本の管が設けられ、それにピペットが取り付けられている。ピペットは
、後方散乱の測定中に激しくポンプ作動させることによって微小球懸濁物を混合
するために使用される。混合および攪拌を行う他の手段には、試料室内のビーズ
状懸濁物と連通した機械式ポンプまたはかきまぜ装置が含まれる。試料室すなわ
ちビーズ状懸濁物を収容している容器は、後方散乱測定を行うための少なくとも
1つの音響アパーチャを有している。音響アパーチャは、最小のたわみで超音波
エネルギを伝達できる薄膜である。それは、較正処理中に微小球懸濁物を収容す
るために使用される。音響アパーチャを形成するのに適した材料は、ポリエチレ
ンフィルム(例えばサランラップ)およびマイラーフィルムを含み、これらの厚
さは全て0.001インチより小さい。微小球懸濁物は、水、生理食塩水、アル
コール、緩衝剤または微小球懸濁物を準備して混合できる適当な他の液体等の液
体環境で準備される。液体環境は、界面活性剤等の他の成分を任意に含有するこ
とができる。1つの好適な実施例では、液体環境は、生理食塩水(コールター(C
oulter)(商標)アイソトン(Isoton)II)と界面活性剤とを含む。本発明の較正
方法は、超音波が伝播する媒体として機能すると共に、較正処理中の微小球懸濁
物の温度を調整するために、一定温度の水浴内で実施されることが好ましい。
本発明の方法に関連して有用な微小球は中空でも中実でもよい。好ましくは、
微小球は、デューク・サイエンティフィック(Duke Scientific)社から入手でき
るポリスチレン微小球CS03MCである。微小球は、実施する較正の形式に合
わせて選択された寸法分布にできる。これは、個人に実施される診断超音波処置
の種類に応じて決めることができる。
本発明の方法に従って微小球の懸濁物からの後方散乱データが得らた時、超音
波装置、システム、器具および設備は当該技術の専門家には既知のプロトコルに
よって所定の後方散乱設定に較正される。簡単に説明すると、本発明の方法が実
施され、実施中の診断処置に最適結果を得るために、較正中の超音波装置の適当
な設定、例えば利得、掃引利得、強度、抑圧レベル、周波数およびトランスジュ
ーサが所定の後方散乱設定値に固定される。当該技術の専門家には理解されるよ
うに、これらの設定値は、検査中の個人の状態、検査中の組織、対照試薬の使用
および所望の情報に応じて決まる。
中実球からの散乱は、音響学の多くの研究方法の原則である。中実球からの散
乱の最初の完全な説明は、ファラン(Faran)によるものと言われている。流体内
のプラスチックおよび金属の中実球により、様々な実験における超音波散乱の物
理的過程を理解するための簡単なシステムが提供された。初期の研究は、ソナー
(水中音波探知機)における散乱体の効果に対する興味から行われた。現在の興
味の領域は、超音波と中実弾性球の相互作用を利用して、海洋環境における堆積
物移送プロセスを予測することである。従来の測定では、当てる超音波の波数を
k、散乱体の直径をaとした時、広範囲のkaに対する結果が得られた。プラス
チック球の攪拌懸濁物を使用する主な利点は、懸濁物内の散乱体の移動および分
布による固有集団加算平均にある。懸濁物内の音の減衰および懸濁物内での共鳴
散乱の調査から、幾何学的断面積の5倍まで増加した後方散乱断面積が、球の剛
体共鳴およびモード共鳴に対応した周波数で発生し得ることがわかっている。
この研究は、超音波対照試薬の物理特性の測定に使用される実験システムを実
証するためのものであった。このシステムは、従来はアルブネックス(Albunex)
(登録商標)および他の対照試薬の物理特性を分析するために使用されていた。
幾つかの理由からプラスチック球が使用された。すなわち、それらの形状はほと
んどの中実および気泡に基づいた対照試薬に類似していること、気相微細泡対照
試薬とは異なった物理的理由からであるが、散乱体が共鳴を示すこと、プラスチ
ック微小球は試薬と同様な寸法範囲で形成できること、再現できるように懸濁物
を作成できることである。すべての測定は、1回の測定で5ないし35MHzの
データを得るための広帯域パルスによって励起される1対の組み合わせトランス
ジューサで行われた。筆者が知る限りでは、単一のパルス実験でそのように広い
周波数帯域にわたる測定値を得るために単一トランスジューサ対を使用すること
は新規の方法である。従来の測定では、研究者は広帯域パルスまたは掃引周波数
実験で数組の組み合わせトランスジューサを使用していた。
この研究はまた、簡単な一次広帯域整理技術を実証するためのものであった。
この技術は、トランスジューサのビームプロフィルおよび信号の時間ゲート処理
(time gating)の効果に関して幾つかの簡略化近似を行う。ポリスチレン球懸濁
物の十分に制御された実験システムは、これらの近似の有効性の検査を助けた。
データ整理技術は計算が集約的ではなく、将来は生体内組織特徴付け用途に使用
することができるであろう。
理論
中実球からの音の散乱の理論予測はよく知られている。中実球からの音の散乱
の予測値は、三次元の波動方程式を解き、中実球の物質の物理的境界条件を利用
することによって得ることができる。3つの関連境界条件は、i)流体内の圧力
が接合面における固体内の応力の法線成分に等しくなければならないこと、ii
)流体の変位の法線(半径方向)成分が固体の変位の法線成分に等しくなければ
ならないこと、iii)剪断応力の接線成分が固体の表面でゼロにならなければ
ならないことである。
これらの境界条件と波動方程式を使用することによって、入射平面波による球
からの散乱について正確な式を引き出すことが可能である。球からの散乱による
圧力についての遠視野式は次の通りである。
但し、P0は入射平面波の振幅、kは入射音波の波数、rは球と受信トランスジ
ューサの間の距離、θは入射波と散乱波が成す角度(我々の後方散乱測定では1
80)、δnはn番目の散乱波の位相ずれである。位相ずれδnはka、球の材料の
ポワソン比および周囲流体に対する球の密度比の関数である。
多重散乱による効果を無視すれば、多くの球の後方散乱を1つの球による後方
散乱から計算できると仮定する。特に、密度数がnの同一球からの散乱断面積ηtot
(f)は次の通りである。
但し、ηsphereは大きさがaの1つの球による散乱断面積である。nの大きさの
分布w(an)を挿入すると、全体散乱の式は次の通りである。
この式を散乱体の大きさと組み合わせて使用して、本文に後述する理論予測値
を計算した。製造者および文献の両方から引き出された値を使用して、懸濁物か
らの後方散乱を計算した。その値は、粒子密度が1.055g/cm3、アイソ
トン密度が1.0g/cm3、アイソトン内の音の縦速度が1490m/s、ポ
リスチレン微小球内の音の縦速度が2380m/s、ポリスチレン微小球の剪断
速度が1180m/sであった。実験結果と理論予測値がさらに一致するように
、ポリスチレンの剪断速度の値を従来の報告値の1100から1120m/sに
調節する。
方法
A.実験
ポリスチレン微小球(デューク・サイエンティフィック社、CS03MC)の
懸濁物の減衰および後方散乱の一連の測定を、超音波対照試薬の特性の測定用に
構成された実験システムで行った。10μmを1段階とした7区分の30から1
00μmにわたる平均直径を有するポリスチレン微小球を使用した。ポリスチレ
ン微小球は光学分粒器(アキュサイザー(AccuSizer)(商標)770、パーティ
クル・サイジング・システムズ社(Particle Sizing Systems,Inc.))で分粒され
、それらの平均直径のほぼ5%の半幅値の分布を有することがわかった。分粒器
は
単一粒子光学感知装置を用いる。液状懸濁物内の粒子は小さいフォトゾーン、す
なわち、レーザダイオードからの光で生じた均一照度の狭い厚板状の領域を通っ
て流れることかできる。フォトゾーン内の合致粒子を避けるために、粒子懸濁物
の希釈と流量の両方を調節した。感知ゾーンを粒子が通過することによって、検
出パルスが発生し、これの大きさは粒子の平均直径によって決まる。
ポリスチレン微小球の減衰を測定するために使用される実験機構の概略が図2
に示されている。パルサー/レシーバ(パナメトリックス(Panametrics)580
0)を使用して広帯域PVDFトランスジューサ(直径が0.4インチ、焦点距
離が3.25インチ、パナメトリックスZF3002−SU)を駆動した。適合
したPVDFトランスジューサが信号を受け取って、ケーブルキャパシタンスを
最小に抑えるために短いケーブルを介してデジタル化オシロスコープ(テクトロ
ニクス(Tektronics)2440)の高インピーダンスチャネルに連結した。超音波
が懸濁物を通過することによって生じる位相ひずみを最小限に抑えるため、試料
室を受信トランスジューサの近くに配置した。受け取った波形を250メガサン
プル/秒でデジタル化して、IEEE−488バスを介してコンピュータ(マッ
キントッシュ(Macintosh)のカドラ(Quadra)800)へ送信して記憶させ、オフ
ライン解析できるようにした。バスの速度は、毎秒ほぼ30回1024点トレー
スを得ることができるようにしている。減衰および後方散乱の両方の測定用に合
計で1000トレースを獲得した。
懸濁物からの後方散乱の測定を図3に示されている同様な機構によって行った
。パルサー/レシーバ(パナメトリックス5800)を使用して広帯域PVDF
トランスジューサ(直径が0.4インチ、焦点距離が3.25インチ、パナメト
リックスZF3002−SU)を駆動した。試料室の前壁からの反響の後に焦点
領域が発生するようにトランスジューサを配置した。トランスジューサの焦点で
は、照射音波は振幅および位相がほぼ一定であるため、データ整理技術で幾つか
簡略化が可能である。パナメトリックス5800を用いて後方散乱信号を40d
Bだけ増幅した。次に、デジタル化オシロスコープで波形を捕らえて、コンピュ
ータへ送った。
次に、図1に示されている試料室について説明すると、室はサランラップ(商
標)からなる2つの音響窓を2つのレクサン(Lexan)(登録商標)圧力プレート
で保持して構成されている。マイクロピペットを室の1つのポートに取り付けら
れた1本のタイゴン(Tygon)(商標)管に固定して、測定中を通して微小球懸濁
物を攪拌するために使用した。一定の攪拌によって、各超音波照射で新しい散乱
体集団を測定することができた。医療用撮像装置からの画像の調査および室の視
覚による検査によって、散乱体の空間分布が室全体で均一であることがわかった
。較正のため、音響窓の下方の微小球の代わりに鋼板を使用した。鋼板からの後
方散乱信号を以下の見掛け後方散乱伝達関数の計算に使用した。
B.データ整理
後方散乱および減衰の測定値を球形散乱体分布の固有後方散乱係数に整理する
ために、一次広帯域データ整理方法を使用した。整理前の、すなわち後方散乱信
号は、トランスジューサおよび付随した電子機器の電気−機械応答による効果、
周波数依存ビームボリューム、懸濁物による減衰および懸濁物の固有散乱の効果
の絡み合ったものであると考えることができる。電子機器による効果を補償する
ため、後方散乱信号を研磨鋼反射体から戻った信号と比較した。懸濁物から戻っ
たスペクトルを鋼板反射のスペクトルによって正規化することによって、測定装
置の影響を除去することができる。
周波数依存ビームボリュームおよび減衰の効果は、幾つかの簡略化を行うこと
によって除去することができる。第1の簡略化では、ビームボリュームを周波数
依存ビーム幅で側方の、またゲート持続時間によってビーム軸線方向の境界を定
める。第2の簡略化では、散乱体を散乱ボリューム全体に均一に分布しているこ
とを医療用撮像装置および目視検査で確認した。これらの近似の両方を使用して
、後方散乱係数は次のように表すことができる。
但し、η(f)は周波数依存後方散乱係数、Γは我々の研磨鋼反射体の振幅反射
係数、Tは水およびポリスチレン微小球懸濁物間の透過係数、Rはトランスジュ
ーサの焦点距離、V(f)は周波数依存ビームボリューム、S(f)は我々の懸
濁物の見掛け後方散乱伝達関数、F(α、τ)は減衰補正率である。F(α、τ
)
は次のように表される。
但し、αは懸濁物の周波数依存減衰、cは懸濁物内の音の速度、xは室の前壁の
後方の窓の位置、τはゲート持続時間である。このデータ整理方法は、球形散乱
体の集合による散乱の実験結果と理論予測値を比較するために使用される。
C.データ解析
2.56μ秒のハミング窓を室壁からの反響の1.0μ秒後の位置でのトレー
スに当てはめることによって見掛け後方散乱伝達関数を計算した。後方散乱係数
に整理した時、他の組み合わせの窓位置および長さも同様な結果を生じた。ゲー
ト処理したトレースのパワースペクトルを計算し、研磨ステンレス鋼板のパワー
スペクトルで正規化して平均した。
懸濁物が室内にある時に透過信号の10の測定値から見掛け減衰係数を計算し
た。各測定値は、波形の100回の領域平均からなる。室に生理食塩−緩衝液で
あるアイソトン(商標)を満たして透過パルスを測定することによって、較正追
跡を行った。10の平均懸濁物の経跡のパワースペクトルを生理食塩水トレース
のパワースペクトルから減じることによって、信号損失を周波数の関数として得
た。標準的なパルス反射法を用いて調べた室の厚さによって信号損失を正規化す
ることによって減衰係数を測定した。
結果
図4は、呼び径が50μmの微小球を光学的に分粒することによって調べたポ
リスチレン微小球の大きさの分布を示している。平均測定直径は49.5μm±
1.0μmであり、これによって後方散乱係数の明瞭な共鳴を測定することがで
きる。
図5は、典型的な見掛け後方散乱伝達関数、見掛け減衰係数および後方散乱係
数の実験測定値を示している。上と中間のグラフはそれぞれ後方散乱および減衰
の直接的な測定値であるが、下のグラフは上記のように後方散乱係数の形に整理
されている。すべての誤差バールは測定値の標準偏差を表している。
散乱体の大きさの分布を散乱モデルの理論予測に導入した結果が図6に示され
ている。理論予測は方程式(3)を用いて行い、分布W(an)は図3のデータ
から計算した。各々の微小球の直径はそれ自身の共鳴系列を示すので、散乱体直
径の分布を導入する結果として共鳴が平均化され、共鳴の幅が広くなる。
図7は、狭い分布の微小球直径の様々な懸濁物の散乱モデルの実験的に測定さ
れた後方散乱係数と理論予測値を比較している。実験と理論の一致は非常によい
ことであるが、後方散乱スペクトル内の共鳴に近い後方散乱係数の大きさにおい
て幾つかの不一致がある。
検討
超音波対照試薬の物理特性を測定するための広帯域の単一トランスジューサ対
の実験装置が開発されている。この装置は5から35MHzの周波数範囲のポリ
スチレン微小球懸濁物の後方散乱係数を測定するために使用され成功している。
広帯域測定に対して単一のトランスジューサ対を使用することは、幾つかの理由
から好都合である。それは対照試薬等の時間変化物質の迅速かつ確実な調査を可
能にする。また、1組のトランスジューサの特性の知識だけが必要であるため、
トランスジューサのビームプロフィルおよび焦点の効果の補償は各組の測定値で
同様である。
そのように広い帯域幅にわたるポリスチレン微小球を確実に測定できることか
ら、球の攪拌懸濁物は高周波超音波システム用の有用な体外較正装置になるであ
ろうと考えられる。この実験装置は、幾つかの理由から有用な較正装置として機
能でき、その一部は組織模倣ファントムと組み合わせて他の者によって以前に調
べられている。しかし、我々は、静止した散乱体と共に従来のファントムを使用
することで必要とされる空間平均化に比べて改良された後方散乱の集団平均化に
つながる新規な混合技術を導入することができた。
粒子懸濁物が鋼板より良好な較正装置になる理由は幾つかある。後方散乱が懸
濁物塊から戻るのに対して、鋼板からの後方散乱は照射高周波の平面反射だけを
表す。研磨鋼板からの反射の極めて高い振幅とは異なって、球懸濁物からの散乱
は、組織塊の内部からの散乱と同じ程度の大きさを有する。そのため、そのよう
な微小球を使用することによって、生体内に使用される振幅範囲にわたって超音
波装置の較正を行うことができる。これによって、高周波走査ヘッドまたは血管
内超音波カテーテルの各使用前の較正が容易になると共に、異なった時間に、ま
たは異なった調査者または臨床医が行った調査間の比較が可能になる。音響顕微
鏡検査の分野では、中実の弾性球の懸濁物も有用であることがわかるであろう。
現在の較正標準は、プレートの表面の整合に対して敏感な研磨ステンレス鋼板か
らの反射である。球の懸濁物の後方散乱を測定する時、基本散乱体が対称である
ことから整合はさほど重要な問題ではない。
超音波対照試薬の分野では、厳密に分粒されたポリスチレン球は、対照試薬の
効果の将来は体内、現在は体外の両方での測定値に対する可能な較正法を提供し
ている。球懸濁物は、対照試薬の注入前に超音波測定装置を較正できるようにす
るための体外設定に使用することができる。また、我々の技術は、特定の対照試
薬の特性の測定に応用することができる。実際に、患者の長期的研究用に考えら
れる較正ツールとして体内に対照試薬を使用することは興味のある分野であり、
肺の濾過および気泡破壊による効果を考慮に入れることができれば、異なった患
者間で音響測定値を比較する方法にもなるであろう。
我々は、直径が30μmから100μmのポリスチレン微小球懸濁物からの後
方散乱を調べた。医療用超音波対照試薬の分野では、肺の血管床(bed)を通過す
るために、10μmより小さいサイズであることが必要である。しかし、後方散
乱係数の測定値と理論予測値の一致は、我々の実験機構で対照試薬の直接測定を
行うことが実現可能であることを示している。我々が研究した濃度が単一散乱形
態にあって、多重散乱はほとんど、または全く発生しないと仮定した。散乱の大
きさは比較的小さく、従って多重散乱は後方散乱の最終測定値に大して影響を与
えないので、この仮定は妥当であろう。別の制約は、各寸法範囲に対して1つの
微小球濃度のみを使用したことであった。中実球懸濁物からの散乱は強い温度依
存性を示すことが別の研究者によって示されている。我々は、各寸法の微小球の
測定中に水浴を一定温度に維持したが、測定間の温度は1.6℃まで変化した。
較正標準として懸濁粒子を使用するためには、注意深い温度管理が必要であろう
。
モデルでは、後方散乱係数スペクトルにおける共鳴の位置は、散乱体を構成す
る物質の横方向速度に対して選択された値に非常に敏感である。ポリスチレン微
小球の場合、1100ないし1120m/sの範囲で変化する横方向速度が報告
されている。本研究では、我々の測定した共鳴と重なるように「理論」共鳴を調
節するために、1180m/sの値を選択した。散乱体の大きさには関係なく、
1180m/sの横方向速度では、実験で測定された共鳴と一致する理論予測値
が得られた。従来の文献に報告されているポリスチレンの横速度と本研究で測定
したものとの値の不一致は、幾つかの要因によるものであろう。ポリスチレン微
小球は、ポリスチレンとポリスチレンジビニルベンゼンの混合化学組成を有する
ことが製造者によって報告されており、超音波では十分に特徴づけられていない
システムであった。
本研究に示されているデータは、後方散乱係数の簡単な一次広帯域整理法の使
用によって、理論予想値との量的に優れた一致が得られることも暗示している。
高周波では、焦点を合わせた超音波トランスジューサの被写界深度が非常に小さ
いことから、多くの研究者が超音波ビームを周波数および空間の関数としてモデ
ル化するための高度な補正技術の使用方法を研究している。後方散乱のここに示
されている測定値と理論的に計算した値が見掛けで一致していることは、ビーム
ボリュームをモデル化する簡単な方法がこの周波数範囲の多くの測定には十分で
あることを示している。
広帯域にわたってポリスチレン微小球を確実に測定できることは、球懸濁物が
高周波の超音波装置用の有用な較正装置になることができることを示している。
Claims (1)
- 【特許請求の範囲】 1.少なくとも1つの音響アパーチャを有する容器内に所定寸法分布の微小球の 、液体環境で懸濁した懸濁物を準備する段階と、 容器を一定温度の水浴に入れる段階と、 容器内の微小球懸濁物を連続的に攪拌する一方、同時に微小球懸濁物からの後 方散乱を測定する段階と、 測定された後方散乱に従って選択された後方散乱設定に超音波装置を較正する 段階を含む超音波装置の較正方法。 2.微小球は中実である請求の範囲第1項に記載の方法。 3.微小球は中空である請求の範囲第1項に記載の方法。 4.微小球は中実のポリスチレンビーズである請求の範囲第1項に記載の方法。 5.微小球懸濁物は、微小球懸濁物を収容している液体環境のポンピングを行う 手段によって連続的に攪拌されるようにした請求の範囲第1項に記載の方法。 6.微小球懸濁物を収容している液体環境のポンピングを行う手段は、手動また は機械式ポンプである請求の範囲第5項に記載の方法。 7.微小球懸濁物は、微小球懸濁物を収容している液体環境の混合を行う手段に よって連続的に攪拌されるようにした請求の範囲第1項に記載の方法。 8.微小球懸濁物は、微小球懸濁物を収容している液体環境のかきまぜを行う手 段によって連続的に攪拌されるようにした請求の範囲第1項に記載の方法。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US1909996P | 1996-06-03 | 1996-06-03 | |
US60/019,099 | 1996-06-03 | ||
PCT/US1997/008803 WO1997046159A1 (en) | 1996-06-03 | 1997-05-30 | Method for calibrating ultrasound devices and contrast agents and system therefor |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2000511449A true JP2000511449A (ja) | 2000-09-05 |
Family
ID=21791432
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP10500624A Withdrawn JP2000511449A (ja) | 1996-06-03 | 1997-05-30 | 超音波装置および対照試薬を較正する方法およびそのための装置 |
Country Status (5)
Country | Link |
---|---|
US (1) | US6148655A (ja) |
EP (1) | EP1006880A4 (ja) |
JP (1) | JP2000511449A (ja) |
CA (1) | CA2257504C (ja) |
WO (1) | WO1997046159A1 (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003325522A (ja) * | 2002-05-14 | 2003-11-18 | Matsushita Electric Ind Co Ltd | 超音波診断装置の評価方法および評価装置 |
JP2017510374A (ja) * | 2014-04-07 | 2017-04-13 | ブラッコ・シュイス・ソシエテ・アノニムBracco Suisse SA | 非基本波分析による音響レベルのIn−Situ推定 |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6318146B1 (en) | 1999-07-14 | 2001-11-20 | Wisconsin Alumni Research Foundation | Multi-imaging modality tissue mimicking materials for imaging phantoms |
US7085406B2 (en) * | 2001-07-27 | 2006-08-01 | General Electric Company | Method and system for unsupervised transfer function generation for images of rendered volumes |
US20040067591A1 (en) * | 2002-10-04 | 2004-04-08 | Wisconsin Alumni Research Foundation | Tissue mimicking elastography phantoms |
US7874990B2 (en) * | 2004-01-14 | 2011-01-25 | The Cleveland Clinic Foundation | System and method for determining a transfer function |
DE102004011744A1 (de) * | 2004-03-03 | 2005-09-22 | Aesculap Ag & Co. Kg | Chirurgische oder medizinische Vorrichtung und Verfahren zur Kalibrierung eines Ultraschallsensors |
US7739911B2 (en) * | 2006-09-15 | 2010-06-22 | Battelle Memorial Institute | Ultrasonic characterization of solid liquid suspensions |
US7569809B2 (en) * | 2006-09-21 | 2009-08-04 | Siemens Medical Solutions Usa, Inc. | Flow-drive for dynamic phantom |
WO2008048333A2 (en) * | 2006-10-18 | 2008-04-24 | Luna Innovations Incorporated | Method and apparatus for calibrating an ultrasonic sensing system used to detect moving objects |
US10130342B2 (en) | 2007-12-28 | 2018-11-20 | Bracco Suisse Sa | Initialization of fitting parameters for perfusion assessment based on bolus administration |
DE102008027367B4 (de) * | 2008-06-09 | 2015-01-29 | Siebo Hicken | Vorrichtung zur Messung der Ultraschallabstrahlung von Ultraschallgeräten |
KR20140113173A (ko) * | 2013-03-15 | 2014-09-24 | 삼성전자주식회사 | 팬텀, 상기 팬텀을 포함하는 초음파 시스템, 및 이의 제조방법 |
EP3386397B1 (en) | 2015-12-10 | 2020-02-05 | Bracco Suisse SA | Detection of immobilized contrast agent with dynamic thresholding |
Family Cites Families (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3779070A (en) * | 1971-11-23 | 1973-12-18 | Autometrics Co | Particle size and percent solids monitor |
US4277367A (en) * | 1978-10-23 | 1981-07-07 | Wisconsin Alumni Research Foundation | Phantom material and method |
FR2571980A1 (fr) * | 1984-10-22 | 1986-04-25 | Extramet Sa | Procede et dispositif de fabrication de micro-billes calibrees et micro-billes obtenues. |
DE3438798C2 (de) * | 1984-10-23 | 1987-01-02 | Löffler, Friedrich, Prof. Dr.-Ing., 7500 Karlsruhe | Verfahren und Vorrichtung zum Messen der Feststoffkonzentration und der Korngrößenverteilung in einer Suspension mittels Ultraschall |
US4843866A (en) * | 1986-10-06 | 1989-07-04 | Wisconsin Alumni Research Foundation | Ultrasound phantom |
US5922304A (en) * | 1989-12-22 | 1999-07-13 | Imarx Pharmaceutical Corp. | Gaseous precursor filled microspheres as magnetic resonance imaging contrast agents |
US5542935A (en) * | 1989-12-22 | 1996-08-06 | Imarx Pharmaceutical Corp. | Therapeutic delivery systems related applications |
DE4019025C2 (de) * | 1990-06-14 | 1994-09-22 | Triton Technology Inc | Verfahren zur Messung der Durchblutung von Organ- und/oder Gewebeproben |
US5342283A (en) * | 1990-08-13 | 1994-08-30 | Good Roger R | Endocurietherapy |
US5874062A (en) * | 1991-04-05 | 1999-02-23 | Imarx Pharmaceutical Corp. | Methods of computed tomography using perfluorocarbon gaseous filled microspheres as contrast agents |
US5433207A (en) * | 1993-11-15 | 1995-07-18 | Pretlow, Iii; Robert A. | Method and apparatus to characterize ultrasonically reflective contrast agents |
US6458762B1 (en) * | 1994-03-28 | 2002-10-01 | Baxter International, Inc. | Therapeutic use of hemoglobin for preserving tissue viability and reducing restenosis |
US5560242A (en) * | 1994-08-16 | 1996-10-01 | Flextech Systems, Inc. | Ultrasonic system evaluation phantoms |
US5625137A (en) * | 1995-05-25 | 1997-04-29 | Wisconsin Alumni Research Foundation | Very low scatter liquid and solid tissue mimicking material for ultrasound phantoms and method of making the same |
-
1997
- 1997-05-30 JP JP10500624A patent/JP2000511449A/ja not_active Withdrawn
- 1997-05-30 EP EP97929691A patent/EP1006880A4/en not_active Withdrawn
- 1997-05-30 WO PCT/US1997/008803 patent/WO1997046159A1/en not_active Application Discontinuation
- 1997-05-30 US US09/180,484 patent/US6148655A/en not_active Expired - Fee Related
- 1997-05-30 CA CA002257504A patent/CA2257504C/en not_active Expired - Fee Related
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003325522A (ja) * | 2002-05-14 | 2003-11-18 | Matsushita Electric Ind Co Ltd | 超音波診断装置の評価方法および評価装置 |
JP2017510374A (ja) * | 2014-04-07 | 2017-04-13 | ブラッコ・シュイス・ソシエテ・アノニムBracco Suisse SA | 非基本波分析による音響レベルのIn−Situ推定 |
Also Published As
Publication number | Publication date |
---|---|
CA2257504C (en) | 2003-02-25 |
EP1006880A4 (en) | 2005-04-06 |
CA2257504A1 (en) | 1997-12-11 |
US6148655A (en) | 2000-11-21 |
WO1997046159A1 (en) | 1997-12-11 |
EP1006880A1 (en) | 2000-06-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Hoff et al. | Oscillations of polymeric microbubbles: Effect of the encapsulating shell | |
Atchley et al. | Thresholds for cavitation produced in water by pulsed ultrasound | |
Chang et al. | Second harmonic imaging and harmonic Doppler measurements with Albunex | |
Gallot et al. | Ultrafast ultrasonic imaging coupled to rheometry: Principle and illustration | |
Deng et al. | In vitro measurements of inertial cavitation thresholds in human blood | |
Chivers et al. | A spectral approach to ultrasonic scattering from human tissue: methods, objectives and backscattering measurements | |
Stride et al. | Investigating the significance of multiple scattering in ultrasound contrast agent particle populations | |
JP2000511449A (ja) | 超音波装置および対照試薬を較正する方法およびそのための装置 | |
Marsh et al. | Broadband through-transmission signal loss measurements of Albunex® suspensions at concentrations approaching in vivo doses | |
JP2007283119A (ja) | 超音波の方法、および、該方法を実行するための回路 | |
Marsh et al. | Frequency and concentration dependence of the backscatter coefficient of the ultrasound contrast agent Albunex® | |
Sboros et al. | Absolute measurement of ultrasonic backscatter from single microbubbles | |
Mobley et al. | Measurements and predictions of the phase velocity and attenuation coefficient in suspensions of elastic microspheres | |
Hall et al. | Broadband measurements of the attenuation coefficient and backscatter coefficient for suspensions: A potential calibration tool | |
Deng et al. | Inertial cavitation produced by pulsed ultrasound in controlled host media | |
Wear | Measurement of dependence of backscatter coefficient from cylinders on frequency and diameter using focused transducers—with applications in trabecular bone | |
Palanchon et al. | New technique for emboli detection and discrimination based on nonlinear characteristics of gas bubbles | |
Insana et al. | Acoustic backscattering from ultrasonically tissuelike media | |
Chen et al. | A method for determination of frequency‐dependent effective scatterer number density | |
Goertz et al. | High frequency attenuation measurements of lipid encapsulated contrast agents | |
Brandner et al. | Estimation of attenuation coefficient from simulated B mode ultrasound images and tissue mimicking materials | |
JP2007301286A (ja) | 気泡検出装置 | |
Bridal et al. | Frequency dependence of acoustic backscatter from 5 to 65 MHz (0.06< ka< 4.0) of polystyrene beads in agarose | |
Chen et al. | Harmonic vibro-acoustography | |
Schmitz | Ultrasound in medical diagnosis |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20031216 |
|
A761 | Written withdrawal of application |
Free format text: JAPANESE INTERMEDIATE CODE: A761 Effective date: 20040130 |