JP2965153B2 - 骨の状態を生体において評価する超音波評価方法及び装置 - Google Patents
骨の状態を生体において評価する超音波評価方法及び装置Info
- Publication number
- JP2965153B2 JP2965153B2 JP62503764A JP50376487A JP2965153B2 JP 2965153 B2 JP2965153 B2 JP 2965153B2 JP 62503764 A JP62503764 A JP 62503764A JP 50376487 A JP50376487 A JP 50376487A JP 2965153 B2 JP2965153 B2 JP 2965153B2
- Authority
- JP
- Japan
- Prior art keywords
- signal
- bone
- frequency
- received signal
- ultrasonic
- 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 - Lifetime
Links
- 210000000988 bone and bone Anatomy 0.000 title description 155
- 238000011156 evaluation Methods 0.000 title description 5
- 238000000034 method Methods 0.000 claims description 58
- 210000004872 soft tissue Anatomy 0.000 claims description 51
- 238000001228 spectrum Methods 0.000 claims description 19
- 230000000694 effects Effects 0.000 claims description 11
- 239000000463 material Substances 0.000 claims description 7
- 238000012544 monitoring process Methods 0.000 claims description 6
- 230000007704 transition Effects 0.000 claims 6
- 238000005259 measurement Methods 0.000 description 28
- 208000001132 Osteoporosis Diseases 0.000 description 24
- 210000004417 patella Anatomy 0.000 description 18
- 230000005284 excitation Effects 0.000 description 15
- 230000005540 biological transmission Effects 0.000 description 14
- 230000006870 function Effects 0.000 description 14
- 238000002604 ultrasonography Methods 0.000 description 11
- 230000008859 change Effects 0.000 description 10
- 230000001054 cortical effect Effects 0.000 description 10
- 238000005314 correlation function Methods 0.000 description 9
- 238000012360 testing method Methods 0.000 description 9
- 210000002303 tibia Anatomy 0.000 description 9
- 230000037182 bone density Effects 0.000 description 8
- 238000012545 processing Methods 0.000 description 8
- 241001465754 Metazoa Species 0.000 description 6
- 238000010521 absorption reaction Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 6
- 230000008901 benefit Effects 0.000 description 6
- 201000010099 disease Diseases 0.000 description 6
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 6
- 238000002591 computed tomography Methods 0.000 description 5
- 238000003745 diagnosis Methods 0.000 description 5
- 201000001245 periodontitis Diseases 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 238000000926 separation method Methods 0.000 description 4
- 210000001519 tissue Anatomy 0.000 description 4
- 238000012935 Averaging Methods 0.000 description 3
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 description 3
- 102000008186 Collagen Human genes 0.000 description 3
- 108010035532 Collagen Proteins 0.000 description 3
- 206010017076 Fracture Diseases 0.000 description 3
- 229910052791 calcium Inorganic materials 0.000 description 3
- 239000011575 calcium Substances 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 235000013399 edible fruits Nutrition 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000001902 propagating effect Effects 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 208000010392 Bone Fractures Diseases 0.000 description 2
- 208000020084 Bone disease Diseases 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 2
- 229920001436 collagen Polymers 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 230000003750 conditioning effect Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 230000006866 deterioration Effects 0.000 description 2
- 238000013399 early diagnosis Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 238000001727 in vivo Methods 0.000 description 2
- 210000003127 knee Anatomy 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000011164 ossification Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 210000000278 spinal cord Anatomy 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 206010065687 Bone loss Diseases 0.000 description 1
- 206010010214 Compression fracture Diseases 0.000 description 1
- 241000255925 Diptera Species 0.000 description 1
- 208000029725 Metabolic bone disease Diseases 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- IDRGFNPZDVBSSE-UHFFFAOYSA-N OCCN1CCN(CC1)c1ccc(Nc2ncc3cccc(-c4cccc(NC(=O)C=C)c4)c3n2)c(F)c1F Chemical compound OCCN1CCN(CC1)c1ccc(Nc2ncc3cccc(-c4cccc(NC(=O)C=C)c4)c3n2)c(F)c1F IDRGFNPZDVBSSE-UHFFFAOYSA-N 0.000 description 1
- 208000008589 Obesity Diseases 0.000 description 1
- 208000001647 Renal Insufficiency Diseases 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 210000001185 bone marrow Anatomy 0.000 description 1
- 230000037118 bone strength Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 239000004020 conductor Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000013479 data entry Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000002542 deteriorative effect Effects 0.000 description 1
- 238000002405 diagnostic procedure Methods 0.000 description 1
- 235000005911 diet Nutrition 0.000 description 1
- 230000037213 diet Effects 0.000 description 1
- 235000013367 dietary fats Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000695 excitation spectrum Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000005251 gamma ray Effects 0.000 description 1
- 239000010520 ghee Substances 0.000 description 1
- 230000035876 healing Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012623 in vivo measurement Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 210000001847 jaw Anatomy 0.000 description 1
- 201000006370 kidney failure Diseases 0.000 description 1
- 239000010985 leather Substances 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 150000002632 lipids Chemical class 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 210000004373 mandible Anatomy 0.000 description 1
- 210000002050 maxilla Anatomy 0.000 description 1
- 230000005906 menstruation Effects 0.000 description 1
- 239000003607 modifier Substances 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 238000003947 neutron activation analysis Methods 0.000 description 1
- 235000020824 obesity Nutrition 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 208000008798 osteoma Diseases 0.000 description 1
- 208000002865 osteopetrosis Diseases 0.000 description 1
- 230000001009 osteoporotic effect Effects 0.000 description 1
- 210000004197 pelvis Anatomy 0.000 description 1
- 230000003239 periodontal effect Effects 0.000 description 1
- 229910000498 pewter Inorganic materials 0.000 description 1
- 239000010957 pewter Substances 0.000 description 1
- 230000001817 pituitary effect Effects 0.000 description 1
- 230000003449 preventive effect Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 239000012857 radioactive material Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000002829 reductive effect Effects 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 238000012883 sequential measurement Methods 0.000 description 1
- 210000003625 skull Anatomy 0.000 description 1
- 230000005236 sound signal Effects 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 210000002474 sphenoid bone Anatomy 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000002739 subcortical effect Effects 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000000844 transformation 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/36—Detecting the response signal, e.g. electronic circuits specially adapted therefor
- G01N29/42—Detecting the response signal, e.g. electronic circuits specially adapted therefor by frequency filtering or by tuning to resonant frequency
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0875—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of bone
-
- 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/07—Analysing solids by measuring propagation velocity or propagation time of acoustic waves
-
- 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/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/4409—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by comparison
- G01N29/4427—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by comparison with stored values, e.g. threshold values
-
- 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/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/449—Statistical methods not provided for in G01N29/4409, e.g. averaging, smoothing and interpolation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/485—Diagnostic techniques involving measuring strain or elastic properties
-
- 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/011—Velocity or travel time
-
- 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/02483—Other human or animal parts, e.g. bones
-
- 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/028—Material parameters
- G01N2291/02854—Length, thickness
-
- 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/042—Wave modes
- G01N2291/0421—Longitudinal waves
-
- 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/042—Wave modes
- G01N2291/0422—Shear waves, transverse waves, horizontally polarised waves
Landscapes
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Chemical & Material Sciences (AREA)
- Immunology (AREA)
- General Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Rheumatology (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Molecular Biology (AREA)
- Heart & Thoracic Surgery (AREA)
- Orthopedic Medicine & Surgery (AREA)
- Biomedical Technology (AREA)
- Radiology & Medical Imaging (AREA)
- Biophysics (AREA)
- Probability & Statistics with Applications (AREA)
- Acoustics & Sound (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Description
【発明の詳細な説明】
発明の背景
発明の分野
本発明は一般に、骨の状態を定めるための診察方法に
関する。より詳しくは、本発明は、ヒト及び動物の或る
骨の状態を非侵襲的に生体内において定める方法を開示
する。この情報は例えば骨粗鬆症(異常に低い骨密度)
又は骨化石症(異常に高い骨密度)のような骨の状態の
変化に関係した支障を主体がもつか否かを定めるために
使用しうる。 問題の解析と発明の目的 本発明は、骨の疾病及び他の骨の状態を検査すること
に関する。本発明は、特定的に、ヒトの評価及び診断に
予定されているが、例示的な実施例の明瞭な変更によっ
て、動物の骨の評価にも使用しうる。主な目的は、骨粗
鬆症の存在の診察である。しかし腎不全及び骨の劣化に
よる歯周炎の後の骨の状態の変化を含む他の状態も、こ
の方法を用いて検出しうる。本発明は骨折の検出と、骨
折の治ゆの定量的な評価にとっても有用である。 ヒトの骨格は、筒状骨(長くかたい骨)と、海綿状骨
とから成っており、その各々は、特定の比率の皮質骨
(コンパクトな骨)と梁骨(連結ストランド)とから成
っている。主に皮質(cortex)から成る筒状骨は、四肢
を形成する中垂性骨格を支配し、主に梁骨から成る海綿
状骨は、脊柱及び骨盤の軸骨格を支配する。 皮質骨においても、梁骨においても、コラーゲンのフ
イブリルが全体に延在している。両者の差異は、実際に
は、程度の差異であり、コラーゲンのフイブリルの網状
組織が取る形状に依存する。フイブリルは、一般に、網
状組織がロッドの網状組織であるように分離されていて
もよい。フイブリルはまたより密に広げられることによ
って、ロッドによって結合されたプレートの網状組織の
外観を呈してもよい。どちらの場合にも、骨の或る容積
は、鉱物化コラーゲンフイブリルから成り、或る容積
は、「骨髄」と呼ばれる流体から成る。 或る特別の骨が梁形成を取るか又は皮質形成を取るか
は、骨の果す機能に大きく依存する。骨の形成される仕
方は十分明らかではないが、骨のゆ着の過程は、骨が経
験する応力に或る意味で応答する。従って、比較的高い
応力を経験する骨の領域、例えば足の脛骨の骨幹は、皮
質骨に向う傾向を示す。低い応力を経験する骨の領域は
より梁骨に向う傾向を示す。梁骨の大部分の個所におい
て、梁物質は、皮質骨の比較的薄い層で囲まれている。
膝蓋骨は、多くは梁骨であり、皮質下層は、前面(ante
rior surface)の直下にある。 本発明の主な目標である骨粗鬆症は、一般に年を取る
につれて人を悩ます原因不明の疾病である。骨粗鬆症
は、男性よりも女性に多発し、女性については、月経閉
止の後に多発する。白人の女性は、他の人種の女性より
もこの疾病にかかり易い。骨粗鬆症は、骨組織物質の絶
対的な減少によって現出される。しかし残りの骨は正常
である。骨粗鬆症にかかった人は、皮質骨よりも骨梁な
いし梁骨をより高い比率で失なう。骨粗鬆症の普通の症
状は、せむしであり、これは脊髄の圧迫骨折、首の骨
折、上大たい骨及びラジアス(radius)の末梢端(肱
骨)の骨折によって生ずる。 骨粗鬆症は、その初期においては物理的変形が明らか
でないため、特に油断のならない状態である。骨粗鬆症
は徐々に進行するので、早い時期の診断と適正な治療と
は重大な症状を改善する。更に、骨が形成されるモード
のため、患者の比較的若い時期の増大した運動の習慣と
それに組合された適切なダイエットとは、この状態の効
果を最小にするとも考えられている。婦人の場合、月経
閉止前にこれらのプログラムを開始することが適切であ
る。骨粗鬆症の診断と治療とは、各々の患者が異なった
「正常な」骨密度をもつという事実によって複雑化す
る。そのため、個々の患者の骨の或る特性(即ち、密度
又は弾性)の変化の履歴を作成し、少くとも部分的にこ
の履歴の変化に基づいて診断することが有用である。 本発明の別の目標である歯周炎は、上顎骨と下顎骨の
損失と、それに伴なう歯のゆるみとをひき起こす。従来
は、顎の骨の欠損はX線検査によっていたが、X線検査
は、相当な骨の欠損が既に起こった後にした歯周炎の存
在を示すことができない。 そのため、骨の状態の変化(骨物質の欠損及びこれに
附随する骨の密度及び弾性の減少を含むが、これのみに
は限定されない)を検出する手段を用意することが望ま
しい。最適には、この手段は、非侵襲性、正確、高感度
で、使用し易く、入手が容易であるものとする。これは
一般には、従来の技術では実現できない。 骨粗鬆症の早期の診断のために、いくつかの方法が提
案されている。これらの方法には、全身体カルシウム中
性子活性化分析(TBC)、肱及びオスカルクス(oscalcu
s)の単一光子吸収測定、(脊柱及び大たい骨の首部
の)複式光子吸収測定、コンピューター支援断層撮影
(CATスキャンニング)並びに超音波分析法が含まれ
る。 TBC法は、要約すると、次の通りである。TBCには、約
14MeVのエネルギーの高速中性子源を必要としている。
中性子ビームは、被検体に向けられる。或る原子の核
は、入射中性子ビームのエネルギーに依存して、中性子
を捕そくし、励起される。励起された同位体は、直ちに
(10〜12秒のオーダー)か又は活性化生成物に特有の減
衰期間の後に、1以上のγ線を放出することによって、
安定な状態に戻る。放出されたγ線のエネルギーは、タ
ーゲットの元素を特徴付ける。カルシウムは、3.1MeVの
エネルギーにおいてγ線を放出する。γ線の強度即ち単
位容積当たりのカウントレートは、存在している問題の
元素の量を指示する。TBCの欠点は、全人体カルシウム
のみを好都合に測定するので、骨粗鬆症の初期に依存し
うる骨物質中の局在化された変化が識別できないことに
ある。またTBC装置は、原子炉、サイクロトロン又は放
射性物質のような中性子源を必要とする。また、大形の
複雑の機会が必要なため、比較的高価となり、数も少な
い。そのため多くの患者にとって利用可能ではない。ま
たこの方法は、期待されるように、比較的高価である。
また0.3〜0.5レムの比較的高い放射線量が必要となり、
反復検査ができない。 光子吸収測定は、エネルギーの異なった光子の物質に
よる吸収に依存する技術である。単一光子吸収測定にお
いては、単一エネルギー源が、単色の光子のビームを放
出する。伝送される光子の強度は、光子の走路上の物質
について指数関数的に減衰する。走行距離がわかると、
介在する物質の密度を定めることができる。単一エネル
ギー光子の問題は、骨を囲む軟組織のような多重減衰層
の効果を解消しえないことである。 そのため、複式光子吸収測定法が用いられる。この方
法によれば、2つの別々のエネルギー例えば153Gd又は
2つの別々のエネルギー源例えば125I及び241Amが用い
られる。各々2つの未知について2つの式が結果し、こ
れらが、2つの物質の密度に減衰されない光子の強度と
減衰された光子の強度とを関連付ける。2つの式を解い
て2つの密度を定めることができる。TBCの場合と同様
に、装置のコストが高く、数も少ない。また軟組織は、
皮膚の別々の成分である脂質と筋肉とから成っているの
で、「肥満/やせ」調整因子を使用しなければならな
い。 骨の密度は、CATスキャンを用いることによっても定
めることができる。CATスキャンナーは、X線伝送の多
重の入射から像を計算する装置である。患者は、全周の
回りにX線感知レシーバーを有する中空シリンダー中に
おかれる。X線源は、患者の回りに回転し、患者を通っ
て個別のX線ビームを送出し、このビームは、シリンダ
ーの直径の反対側のレシーバーによって受信される。こ
の伝送はシンリンダーの全周について反復される。受信
センサーが受けるX線の強度は、患者の身体を通る経路
の全密度に関係している。送信器が患者の周囲を移動す
る間にこれらの多数の読みが取られる。複雑なピュータ
ーによる分析によって、各々のレシーバーが受けた読み
が照合され、多くの伝送されるビームによって限定され
る平面と平行に取られた患者の身体の薄いスライスの密
度マップが作成される。 CATスキャンの出力は密度マップである。身体領域を
通るCATスキャンは、骨の密度に関係した情報を供与す
る。CATスキャンの実証された利点は、高再現性−高信
頼性の装置であり、他の関係した分野において十分に試
験されていることである。またそれは、梁骨と皮質の密
度を別々に量化する。そのため、海綿骨の正確な密度
を、周囲のコンパクトな皮質骨の存在による平均化効果
なしに定めることができる。 CATスキャン技法の欠点は、装置がコスト高で、多く
のコミュニティにおいては入手可能ではないことにあ
る。また、可搬型でなく、或る患者にとっては創傷性と
考えられている。CATスキャンの実施は、放射線学上の
有資格の医師又は放射線学者と協同する技師によらねば
ならないため、更にコスト高となる。またCATスキャン
においては患者がX線に露呈されるが、これは、或る医
師及び患者にとっては望ましくないと考えられている。
ここで「患者」とは、ヒト及び動物の被検体を意味す
る。 ヒト及び他の動物の骨を生体において非侵襲的に分析
するために、超音波の使用による測定も行なわれてい
る。しかしこの測定は、多くの骨を囲む軟組織の存在に
よって複雑になる。皮膚を通って投入された超音波信号
の速度は、骨を囲む軟組織の通過によっても、問題の骨
の通過によっても影響される。軟組織の効果を最小とす
ることも可能であるが、そのための従来の試みは、相当
複雑なエコー分析及び注意深い測定の要求によって妨げ
られている。例えば米国特許第4,361,154号「生体内で
骨の強度を定める方法」を参照されたい。 従って、本発明のいくつかの目的は、(1)非侵襲的
な生体内技法によって骨の状態を評定し、(2)実質的
なX線量に患者を露呈させることなく骨の状態を評価
し、(3)患者の骨の状態を経済的に迅速に評価し、
(4)数年間に亘って何回も経済的な安全に実行するこ
との可能な骨の状態の評価方法を提供し、(5)高価な
機器を必要とせずに、特に注意深い測定又は複雑な分析
なしに技師によって実行可能な骨の状態の評価方法を提
供することにある。 ここに使用されている限り、膝蓋骨とは、膝の前方部
分を形成する厚い扁平な3角形の骨であり、「ニーキャ
ップ」(膝蓋骨)は、膝蓋骨と周囲の軟組織であり、脛
骨は足の下部の前部骨であり、「すね」は脛骨及び周囲
の軟組織である。 発明の簡単な要約 本発明によれば、骨の状態は、骨の一部分を通って超
音波信号を送出し、受信される信号を分析することによ
って定められる。超音波信号は識別可能な波形又は周波
数の少くとも2つの成分を好ましくは含んでいる。これ
らの成分はパルス又は連続波のフォーマット又はその或
る組合せである。以下では、特に明記されない限り、
「パルス」という用語は、有限の持続時間の信号と連続
波信号との両方を意味する。狭帯域信号が使用される限
り、これらは、同時に(例えば複合波形の一部として)
又は異なる時間に伝送することができる。別の方法とし
て、広帯域信号を用いてもよい。これは広い周波数帯域
例えば約100kHzから約3MHzのエネルギーを含む信号でも
よい。2つ(又は以上)の成分の伝播の差は、骨の状態
についての情報を含んでいる。 本発明は、媒体を通る超音波信号の伝送速度が媒体の
密度及び弾性に関係しているという原理に基づいて作用
すると考えられる。本発明はまた減衰が周波数依存性で
あることの知見に依存している。皮質骨から梁骨までの
骨の可変の性質のため、超音波信号の速度と媒体の密度
及び/又は弾性との間の関係を正確に表現することはで
きない。しかし、大部分が皮質である骨と大部分が梁骨
である骨について行なった実験の結果とよく相関される
皮質骨と梁骨とについてのモデルが提案されている。即
ち超音波伝播の相対速度は、骨の破損し易さの指示とし
て使用することができる。パルスの送受信のために1対
の変換器が使用される。これらの変換器は、超音波パル
ス発生器と信号受信処理システムとに接続される。これ
らの変換器は、骨の近くの2つの異なった個所に配され
ることにより、超音波信号が骨及び軟組織中に送出され
てこれを通過し、受信変換器によってさえぎられるよう
に配置される。 測定しようとする骨は、好ましくは、最小の軟組織で
囲まれ、その表面は、好ましくは、向い合う部分を備え
ているので、音のエネルギーは、実質的に平行に配され
た送受信変換器の表面によって、これらの部分に入出す
ることができる。従って、超音波信号の第1到来部分の
音の経路は、送信変換器の表面において開始し、骨を覆
う軟組織を通り、骨を通り、次に、受信変換器から骨を
隔てる軟組織を通り、最終的に、受信変換器に到達する
(超音波信号の或る部分は、音響的により長い別の経路
上を伝播し、従って、より遅い時間に到達するであろ
う。)。 いろいろの骨部位を試験の目的に使用してよいが、
(a)膝蓋骨を囲む軟組織層を通る距離が骨を通る距離
よりも著しく短かく、また(b)膝蓋骨及び軟組織層を
通る伝播経路とほぼ平行でしかもこれとほぼ等距離の膝
蓋骨の上方の軟組織を通る音響伝播経路があるため、ニ
ーキャップが有利な個所であることが見出された。その
他の有用な部位の例は、頚骨、上顎骨及び下顎骨があ
る。これらの2つの後者の骨の検査は、歯周炎の発見の
道を開くであろう。 本発明の第1実施例によれば、受信信号は、ノズルを
去するように処理されることにより、超音波伝播遅れ
の持続時間をより正確に定めることを可能とする。変換
器の送信面と受信面との距離は、カリパスによって測定
され、骨を通る速度が計算される。この速度の計算に際
して、骨を覆う軟組織の有限な厚みによって導入される
誤差について測定伝播遅延を最初に修正することが望ま
しい。 骨及び軟組織を通って測定された速度、又は、軟組織
について修正された速度は、種々のデータと比較しう
る。第1に、試験中の患者の時間履歴を作成し、骨の密
度及び/又は弾性の変化を表わす速度の変化をチェック
することができる。この形式の比較は、疾病の経過又は
患者のリハビリをチャートする。また、全ての関係ある
年令の健康な個人及び疾病をもった個人を含めた公衆の
広範な断面について、骨の密度、弾性及び/又は速度を
含むデータベースを作成することができる。選定された
閾値基準に基づいて、このデータベースの情報に、患者
の骨の密度、弾性及び/又は速度の測定値を比較し、患
者の骨が健康か否かを臨床的に有用な信頼性レベルで予
測し、又は、患者が骨の疾病をもっている可能性を評価
することができる。 本発明の第2実施例によれば、送受信変換器の間の距
離の測定は必要ではない。送出変換器によって送出され
る超音波信号は約100kHzから約3MHzまでの帯域内の、少
くとも2つの別々の周波数成分をもっている。患者の軟
組織と骨とを通って伝播された後、この信号の低周波成
分が最初に受信変換器に到達し、おそく到達する高周波
成分が後続する。これら2つの成分の到来時間の比が骨
の状態の好い予測子であることが見出されている。受信
された低周波(即ち1MHz以下)の成分の出力スペクトラ
ムの、同じ帯域内の適用された励起に対する比のよう
な、他の測値も、良い予測子のように思われる。この比
は「利得関数」と呼ばれる。ピーク振幅及びこのような
利得関数の下方の面積も、予測子として使用可能であ
る。更に、これらの種々の予測子は、更に高信頼度の予
測子を実現するために、タンデムに用いてもよい。 本発明の前述した目的並びに他の記述されない目的、
利点及び特徴は、以下の詳細な説明及びここに添付され
た請求の範囲を添付図面との関係において読むことによ
って一層明らかとなるであろう。 図面の簡単な説明 図において、 第1図は、本発明の方法の各ステップを一般的に示し
たフロー図である。 第2図は、第4図の装置に用いられている変換器から
送出される超音波信号のスペクトラムを示す。 第3図は、本発明の方法の第1実施例によるより特定
的なステップの連鎖として、第1図のステップの構成を
より詳細に示したフロー図である。 第4図は、本発明の2つの実施例の方法を実施するた
めに使用可能な装置の略図である(ただし第2実施例に
ついてはカリパスは不要である)。 第5図は、1つの超音波パルスについて受信励起の振
幅対時間の典型的なプロティングを示した略図である。 第6図は、「整合波」の方法を略示する。 第7図は、健康な骨を通過した後の第1到来波形を略
示する。 第8図は、骨粗鬆症の骨を通過した後の第1到来波形
を略示する。 第9図は、本発明の方法の第2実施例によるより特定
的なステップの連鎖として、第1図のステップ7の構成
をより具体的に示すフロー図である。 第10A,10B図は、第9図の方法によって形成された時
間域高周波波形及び時間域低周波波形を略示する。 第11図は、第9図の方法によって得た例示的データを
示す表である。 実施例の詳細な説明 第1図には、本方法の基本的な工程が、概略的なフロ
ー図によって図示されている。符号2によって示した第
1工程は、超音波信号の発生、管理及び測定を実行する
ための適切な装置を供与することを要求する。第4図に
ついて後に一層詳細に説明するこの装置は、一般に、超
音波(例えばパルス)発生器、1対の変換器及び信号検
出−分析ユニットから成っている。 広帯域波形発生器は、約100kHzから約3MHzのエネルギ
ー成分を含む音のパルスを約2×10-6秒の間前記変換器
のうちの1つから発生させるために、ブロック4に示す
ように設けられている。代表的な装置(後述する)から
の励振のスペクトルは、第2図に示されている。このス
ペクトル全体に亘って信号の振幅及び出力が変動可能な
ことは注意されたい(この理由により、約250kHz及び2.
5MHzの別々の狭帯域信号を、前述したように、別の方法
として用いてもよい)。パルスの発生は、データ取得シ
ステムによって記録され、このシステムは、磁気ディス
ク又は他の記録保存可能な媒体上にそれを記憶させる。
これを陰極線管又は他の適当な表示装置上に表示しても
よい。次のステップ6は、受信変換器においての音響励
振のモニタリングを含む。これらの励振は、背景ノイズ
と、骨及び軟組織を経て伝送される第1到来信号と、骨
及び軟組織を経て伝送される第2到来信号とを含む。ス
テップ6は、音響励振に対応する電子信号を発生させ、
データ取得システムによるこれらの励振信号の受信の記
録を記憶させ、同一の陰極線管表示部上にこれらの信号
の受信の記録を表示することを含む。ステップ7では、
ステップ4及び6から得た受信信号は(アナログ又はデ
ジタル波及び種々の計算を含みうる各種の技法によっ
て)、処理され、2以上の誘導された測定を発生させ
る。これらの誘導された測定は、次に、患者を類別化す
る因子、例えば、限定的にではないが、年令,性,人種
及び患者の骨の状態に対して受信信号特性を相関させる
情報データベースと比較される。 ステップ4〜6を数回反復してその結果を平均するこ
とによって、ランダムなノズルによる誤差を減少させる
ことが好ましいが、これは、不可欠ではない。 患者の軟組織及び骨は、部分的には、その伝送特性が
周波数の関数として変化する複雑なネットワークのよう
に作用する。患者を通過する時の超音波入力励振の変更
は、受信側の変換器において、大きく変更されたスペク
トルを供与する。このスペクトルは、異なる時において
到来する少くとも2つの別々のパルス又は信号を含んで
いる。原初の伝送される超音波信号に対して比較される
と共に相互に対しても比較されるこれらのパルスの特性
の、アナログ処理手段又はデジタル処理手段による評価
によって、患者の骨の状態についての情報が得られる。
評価すべき特性には、限定的にではないが、周波数、振
幅及び位相が含まれる。 患者を通って伝搬される時に生ずる超音波励振の複雑
な変化は、未だ十分には理解されていない。現在のとこ
ろ、膝蓋骨については、問題の周波数スペクトルにおい
て、比較的高い周波数成分は、骨を覆う軟組織を経て漏
れモードとして通過し、比較的低い周波数成分は、比較
的密な骨の中の並行な経路に沿って伝搬される。より詳
細には、膝蓋骨については、比較的低い周波数信号は、
平面の直下の比較的密な皮質下骨層中を伝搬されるもの
と思われる。この骨の経路は、加えられた超音波パルス
の高周波成分を著しく減衰させる。他の骨についての変
化は、それほどよく理解されていない。 本発明の第1実施例(第3,4図参照)によれば、変換
器は、デジタルカリパス装置38に取付けられている。カ
リパスは、評価しようとする骨の部分の回りに、変換器
が患者の周囲の皮膚と接触した状態において閉ざされ
る。変換器と皮膚との間の超音波エネルギーの結合は、
インピーダンス整合ゲル又は他の同様の材料によって容
易にされる。励振が発信される受信される時に、キャリ
パーは、ブロック8によって示すように、パルスの通っ
た距離(即ち、変換器の向い合う表面の間の距離)を自
動的に測定する。伝送データ及び励振データが記録され
た後、テンプレート信号が、10によって示したように供
給される。この信号は、受信側の変換器においての励振
に対応する信号に、相互相関機能によって整合される。
相互相関機能の評価によって、ベースライン信号を形成
するノイズの波が可能となるので、パルスの正確な到
来時間をピンポイントすることが容易となる(ブロック
12)。今や、12からのパルス走行の経過時間と、8から
の走行距離との両方が知られたので、骨及び最小限度の
周囲の軟組織を通る超音波パルスの平均速度を計算する
ことができる(ブロック14)。 通常は骨の多い個所である膝蓋骨のところにも、患者
によっては、少くない量の軟かい組織の存在することが
あるので、軟組織の効果を除くために、速度測定を修正
したり調節したりすることが有用となる。本発明に従っ
て、患者の身長−体重比に関係した因子を使用すること
ができる。 速度の測定は、別の方法として、前出の米国特許第43
61154号(引照によって本明細書の一部分となる)によ
り詳細に記載されている。より複雑なエコー技術に従っ
て、患者の軟組織について調節したり修正したりするこ
とができる。この方法は、軟組織と骨との界面を信号が
横切る際に超音波エネルギーの一部分が骨から反射され
るという事実を利用している。各々の変換器から超音波
パルスを発信し、反射エコーパルスを受信するまでの時
間を測定することによって、超音波信号経路の各々の先
端において骨を覆う軟組織の厚みを測定することができ
る。この技法には、両方の変換器から超音波パルスを発
信して両方の変換器においての励振をモニターすること
のできる装置が必要となる。軟組織を通る超音波の速度
は既知であり、約1550m/秒である。従って、超音波エコ
ー発信の記録のチェックによって、超音波パルスが軟組
織を通過し、骨から反射され、発信側変換器に戻るまで
の時間が、各々の変換器について得られる。この走行の
間の速度を知ることによって、経路の長さも定めること
ができる。軟組織を通るこの経路の長さを変換器の間の
全距離から引算すると、骨を通った距離が求められる。
また骨を通る伝搬時間も、前記のようにして定めた軟組
織を通る走行に要した時間を測定された時間から引算す
ることによって計算される。 2つの変換器の間において測定される距離も、標準的
な「皮膚フォルド」試験によって得た因子によって、患
者の軟組織特性に関連して正規化することができる。皮
膚(スキン)フォルド試験によれば、患者のいくつかの
指定された個所において、これらの個所の軟組織の皮膚
量を、カリパスを用いて測定する。一般大衆から多量の
データを収集し、皮膚フォルド試験に基づいて、一般的
な「肉付き」に従って、人々をランク付けした。この試
験に基づいた指標は、骨及び軟組織を通って測定された
速度を正規化するためにも使用することができた。 身長/体重比、皮膚フォルド基準、エコー技法など、
どの方法によって得たものであっても、軟組織について
の調整は、以下に、「軟組織調整因子」と称される。即
ち、測定された絶対速度と、「軟組織調整因子」とか
ら、軟組織調整速度が、ステップ18において定められ
る。 軟組織調整速度は、多くの仕方で使用でき、第3図に
は、そのうちの2つの方法が示されている。第1に、調
整速度は、広い階層の人々から集めた調整速度の編集デ
ータのエントリーと比較される(22)。このデータバン
クは、年令その他の患者の特徴(例えば、性,人種,種
族的集団など)及び既知の骨格の状態に従って組織化す
ることができる。即ち、40才の婦人の調整速度が大衆中
の40才の婦人の既知速度の範囲に含まれれば、この人は
健康な骨をもっていると言えよう。その逆に、調整され
た速度が、健康な60才の婦人に通常組合される範囲に入
っている場合は、異常な状態である確率は実質的に高
く、予防上及びリハビリテーション上の処置を開始する
ことが賢明であろう。 調整速度は、検査中の患者について行なった以前の測
定及び将来の測定と比較することもできる(ステップ2
6)。この計画の下では、何年かに亘って1人の個人に
ついて測定が行なわれよう(ステップ24)。超音波速度
の際立った変化は、骨の状態の際立った変化を表わすも
のとなろう。骨の状態の変化の方向に依存して、患者の
健康が回復途上にあるか又は悪化しているかを定めるこ
とができよう(ステップ28)。 次に、第4図を参照して、本方法を実施するのに適し
た装置について詳細に説明する。各々の変換器40,42
は、広帯域高周波数超音波パルスを発信することができ
る。02254マサチューセッツ州ワルサム,クレセント通
り221,パナメトリックス社によって製造された変換器例
えばパナメトリックスモデルA533Sトランスジューサー
を使用して、良い結果が得られた。これは2.25MHz変換
器として表わされる。2.25MHzは、いくつかの理由によ
って、望ましい基本周波数である。第1に、この周波数
において超音波を送受信しうる変換器は、比較的小形で
あり、直径が2.54cm(1インチ)の約1/3である。その
ため、これらの変換器は、患者の身体のいろいろの部分
に好都合に係合させることができる。第2に、周波数
は、少くとも1×10-6秒の短い時間の測定値を解像しう
る装置の感度が得られるように、十分に高い値としなけ
ればならない。これは、少くとも1MHzの周波数を必要と
する。比較的高周波数の装置は、比較的低周波数の装置
に比べて、時間測定値の解像能力がより高くなる。不幸
なことに、周波数が高すぎると、受信側変換器において
受信される信号の振幅は、骨の内部又は周囲の境界面に
おいてその反射によって相当に減衰する。本発明者ら
は、約2.5MHzの基本周波数が、感度と発信振幅との競合
する配慮の間の妥当な調和点を与えることを見出した。
第2図を参照して、前述した変換器によって送出される
信号のスペクトルが示されている。理解されるように、
このスペクトルは、約100kHzから約3MHzまでの実質的な
エネルギー成分を含んでいる。この広汎なスペクトルが
望ましいのは、骨を経て伝送される信号が、300kHz付近
の周波数において際立ったエネルギーをもち、軟組織を
経て伝送される信号が、2.5MHz付近の周波数において際
立ったエネルギーをもつためであり、そのため、変換器
−信号発生器52は、両方の周波数域において動作可能と
すべきである。これら2つの伝送(即ち、骨及び軟組
織)は、実質的に独立しているので、別々の周波数にお
いて動作する複数の狭帯域ソース又は単一の広帯域ソー
スによって励振された狭帯域変換器によって、同様の結
果が得られる。 変換器(複数)は、02166マサチューセッツ州ニュー
トン,660−・ストリート,フレッド・V・フロウラー社
によって、モデル名“ウルトラ−カル”の下に市販され
ている自動電子デジタルバーニアカリパス38上に取付け
られる。バーニアカリパスは、固定のアーム44と可動の
アーム46とを備えている。可動のアーム46が支持アーム
48に沿って移動すると、変換器面44′,46′の間の距離
が表示部50上に表示される。このカリパスの構造は、開
示された発明にとって大切ではなく、普通のバーニアカ
リパスを用いてもよい。しかしその場合には、距離の値
を手動で読出して記録すると共に、この距離をデータ処
理装置に手動で読込むことが必要となろう。 使用に当り、送信側の変換器40と受信側の変換器42と
は、患者の身体の骨部分60、好ましくは膝蓋骨の両側の
表面上に当接させる。パルス発生器52は、持続時間が約
2×10-6秒のパルスを発生させる。このパルスは、電気
接続部54を経て、発信側変換器40に転送される。それと
同時に、電気接続部56に沿って発信信号がデータ取得シ
ステム58に転送され、発信信号の送出を指示する。この
信号は、送出時を記録するために、信号処理装置によっ
て使用される。データ取得システムについては、マサチ
ューセッツ州ピーボディ,アナロジックス・コーポレイ
ションのデータ・プレシジョン・ディビジョンによって
製造された「データ6000」を使用して、良好な結果が得
られた。 パルスは、測定される部分60を通過し、受信側の変換
器42によって受信される。受信側の変換器42は、導線62
を経てデータ取得システム59に電気的に転送される信号
を送出する。これらの信号は、変換器42によって受信さ
れる超音波パルスに対応している。 データ取得システム58は、それが受信する信号の振幅
を記憶すると共に、標準型のCRT表示部64上に、時間変
化振幅を表示することができる。データ取得システム
は、第4図に略示したように、データ計算手段66を備え
ている。カリパス38によって測定した、変換器面44′,4
6′の間の距離を表わす信号は、後述するように、後の
ステップにおいて使用するために、電気接続部68を経
て、コンピューター66に転送されている。 次の図である第5図には、後述する「整合された
波」の相互相関技法が適用される前の受信側変換器が受
けた典型的な応答が略示されている。第5図において、
ディスプレイは、64により示されている。水平軸線は、
経過した時間を表わしている。原点から楕円644の終点
までの時間は、15−20×10-6秒のオーダーである。垂直
軸線は、受信信号の振幅を表わしている。 円642中に示された信号の部分は、「メイン・バン
グ」と称され、送出パルスの記録を形成する。大体にお
いて楕円644に囲まれた信号部分は、受信側の変換器に
よって受信されたベースラインノイズである。楕円646
によって示された信号部分は、送出パルスの結果として
骨及び軟組織を通過する信号が受信側変換器によって最
初に受信されたことを表わしている。この受信された信
号は、2つの識別可能な成分から成っている。これら
は、最初に到来する低周波信号と、これに続く、後に到
来する高周波信号である(これらは、それぞれ「受信さ
れた第1の識別可能な成分」及び「受信された第2の識
別可能な成分」と称することができる。或る数(例えば
16以上)の別々の初パルスから受信信号の振幅を平均化
することが特に好ましい。ランダムな信号例えばノイズ
は平均化によって零になる傾向をもっている。 第5の略図によって明らかなように、信号646がベー
スラインノイズ644から出る時間を正確に定めることは
困難である。到来時間を定める際の誤差が、1×10-6秒
のように大きいと、速度の計算の誤差(誤り)が10%と
非常に大きくなることがある。この誤差量による測定
は、有用であっても、この誤差を可能な限り除去する
(又は少くとも最小とする)ことがもちろん望ましい。 幸いに、誤り訂正法が開発されている。これは、例え
ばワイリー社から1960年に出版された、Y.W.リー著「通
信の統計理論」のような、種々の刊行物に、十分に記載
されている。誤り訂正方法は、第6図に略示されてい
る。一般に、ベースラインノイズは、信号の平均化及び
相互相関によって、受信信号の表示から去される。骨
を経て受信される平均化された第1到来信号の形状に相
当な程度まで整合されたテンプレート信号が供給され
る。以下において明らかとなるように、骨を通る受信信
号の前縁にテンプレートが正確に同一であることは必要
ではなく、単に類似していることが必要である。テンプ
レートは、その患者について受信された支配的なフーリ
エ成分の適切なサイクル数を取ることによって、例えば
取得される。好ましくは、2つのテンプレート信号が使
用され、その1つは、骨を通る伝送に所属し、第2のも
のは、軟組織を通る伝送に所属する。後者の信号は、標
準的な信号でありうるが、その理由は、軟組織を通る超
音波の速度と軟組織の周波数応答とが患者ごとにほとん
ど相違しないためである。テンプレート内の或る適切な
数のサンプル点を与えるために、2つの異なったテンプ
レート信号について2つの別々のサンプリング間隔及び
/又はサンプリング率を用いることが時に望ましい。 信号の或るわずかな部分、例えば1サイクルの180゜
分が保持される。相互相関関数C(d)は、 として定義される。 関数F(t)は、第6図に示すように、受信側変換器
によって受信された信号である。関数G(t−d)は、
原点の付近においてのみ非零のテンプレート信号であ
る。骨を通る信号の到達時点をピンポイントするため
に、時間が零から無限に及ぶ積分を行なう。次に“d"が
零から無限大に及ぶ間のこの積分の第1極値を選定す
る。 C(d)の絶対値が、人為信号に特有の或る閾値を超
える第1極大値を超過した時点648は、後述する場合を
除いては、信号が最初に到達した時である。これは次の
ことによって生ずる。伝送パルスを最初に受信するまで
は、関数F(t)は、メイン・バングとベースライン64
4との急速に振動する信号から成っている。これらの信
号は一般に約2MHzにおいて振動する。その逆に、骨梁に
とって適切なテンプレート信号は、3.3×10-6秒の周期
をもち、従って周波数は300kHzである。主要なバング及
びベースラインのノイズは、正から負に変化するので、
2つの関数F(t)及びG(t−d)の積の積分は、高
速振動の全ての時において零となる。テンプレートG
(t−d)及びF(t)(今や伝送パルスからの信号か
ら主に成っている)である相互相関関数の積が非零値を
もつのは、伝送パルスが到達(時点648)するまでは起
こらない。受信信号がテンプレートに最も近い値となる
時に、この積は、極大値となる。メイン・バングが第2
図の広帯域波形をもつ代りに複雑の別々の信号から成る
場合に対応した状態となる。 前述したように、相互相関関数の極大値が第1到来パ
ルスの到来と合致しないいくつかの場合がある。これら
は、第1到来パルスの振幅がテンプレートに対しての逆
の符号になっている場合である。例えば、第1到来信号
が負の場合に骨信号の正の区分からテンプレートが選ば
れている場合である。これらの場合は、相互相関関数が
最大値の前に比較的大きな負のエクスカーションをもつ
ことを認識することによって容易に識別できる。この場
合には、第1到来パルスの到来は、大きな負のエクスカ
ーションにおいて表わされる。負のテンプレートをもっ
たものと正のテンプレートをもったものとの2つの相互
相関関数を実行し、観察された結果に基づいた適切な保
持すべき記録を選択することによって、より正確な測定
を行ないうる。また、相互相関関数を平方してもよく、
この場合には、結果関数の極大値は、第1到来パルスの
到来と合致する。テンプレートG(t−d)をC(t)
と組合せてC′(d)を得るように、整合波のプロセ
スを1回以上実行することも時には望ましい。 この信号整合形態は、前述したデータ取得方式によっ
て自動的に実現することができる。 この整合波は、患者の診察の時には行なわずに、記
憶されたデータについて後に行なってもよい。相互相関
関数が極大にされ、信号伝播の経過時間が定められた
後、以前に測定された変換器の間の距離を経過時間で割
算することによって、速度を計算することができる。 変換器間の経路を通る平均速度を定めたら、骨を囲む
軟組織を通る超音波パルスの通過の効果を除くことが適
切である。それは、1つの変換器から次の変換器へのパ
ルスの通過時間には、超音波パルスが骨を通過している
時間と、超音波パルスが骨を囲む軟組織を通過している
時間とが含まれるからである。更に、超音波パルスの走
行距離には、骨を通る距離と、骨の両側の軟組織を通る
距離とが含まれる。軟組織を通るこの通過遅延を勘案し
えないと、単に走行距離を時間で割算することによる速
度の計算の精度が低下する。 骨梁の状態を検査し、骨粗鬆症についてテストするた
めには、膝蓋骨が、人体の中で非常に良い場合であるこ
とが、実験によって示された。膝蓋骨は、軟組織を実質
的に含まない人体の部分である。また膝蓋骨は、ほとん
ど梁骨(骨粗鬆症の効果が最も早く且つ最も深刻に感じ
られる形態の骨)から成っている。本発明者らは、脛骨
稜は、骨腫の状態を調べるための人体の中の良好な場所
であると考える。 前述したように、膝蓋骨及び脛骨稜において測定する
ことに利点があっても、これらの個所においての多量の
軟組織の存在は、結果に多少とも影響する。従って、患
者の特定の軟組織の状態について調整することによっ
て、速度の測定値を修正することが適切である。前述し
たように、これらの方法には、エコー技法、スキンフォ
ルド試験及び身長/体重の比較が含まれる。どの方法に
よって計算した速度も、「調節された骨速度」と呼ばれ
る。 調節された骨速度が1度定められたら、この値を用い
て、患者の骨格状態を評価することができる。これは、
広い範囲の公衆から得た調節された骨速度に前記の調節
された骨速度を比較することによって行なうことができ
る。理想的には、調節された骨速度の一般的なデータベ
ースは年令及び他の特性の範囲について分類し、或る年
令及び特徴の健康な人の典型値が知られるようにする。
このデータベースは、通常よりも多孔質なものから通常
よりも密なものまでの、骨密度の状態のスペクトルに対
応する一連の速度を、各々の年令又は年令群について含
んでいてもよい。患者の調節された骨の速度を彼らの年
令に対する調節された骨速度のスペクトラムと比較する
ことによって、患者の相対的な骨の状態を定めることが
できる。 また、或る個人が、骨粗鬆症の疑いがあったり、その
治療を受けている場合、同じ患者の以前の測定値に、特
定の測定値を比較することができる。これらの測定値
は、系列的な測定として知られる。これにより、その進
行又は劣化を図式化し、また治療することができる。 伝送された波形の振幅の形状に基づいて患者の骨の状
態を評価することも、本発明の方法の1つの特徴であ
る。受信された超音波パルスの波形の形状は、患者の骨
格の状態及び密度についての貴重な情報を含んでいる。
本発明者らは、健康な患者の場合に、受信変換器に最初
に到達する膝蓋骨を通る信号の形状が、約300kHzの周波
数のrmsベースラインノイズの2倍以上の、形状の定ま
った振動の第2部分と、2.5MHzに近い非常に高振幅−高
周波数の振動とによってほぼ特徴付けられることを見出
した。第7図に骨信号のこの形状を示す。第7図におい
て、健康な骨を通過した信号は、“70"により示され
る。第8図において、病気の骨を通過した信号(「病信
号」は、“80"により示される。この図からわかるよう
に、病信号は、高周波数−高振幅において開始する、ほ
ぼ形状の定まっていない振動において到来する。 波形の形状のこの変化に内在する機構と波形の意味と
は、十分には理解されていない。形状の定まった低周波
信号は、縦波又は圧縮波を示し、高周波−高振幅振動は
横波又はせん断波である。縦波は、伝播方向と平行な方
向の媒体の振動から結果し、せん断波は、伝播方向と直
角の振動の実証である。縦波は一般に横波よりも高速で
ある。しかし骨粗鬆症においては、梁骨部分を通る圧縮
超音波パルスの走行路は、より多く制限されており、従
って、梁骨形成物の終端と介在する流体ポケットとによ
って中断されると考えられる。媒体においてのこれらの
不連続は、骨の弾力を変更させ、超音波パルスに合体さ
れたエネルギーの放散を生ずるので、骨粗鬆症の患者の
場合、圧縮超音波パルスの振幅は減少する。 また、変換器が頭部波及びせん断波を開始(即ち送
出)するようにしてもよい。これらの波とそれらが誘起
する別の波との伝播速度は、変換器の間の骨の状態及び
性状についての別の情報も供与する。或る骨は多くの異
なったモードで励起されうる。関心のある骨の特性(例
えば、多孔性、密度及び弾性)について特に敏感な伝送
のモードを見出すことが目標である。 これらの他の波形式を発生するために、第4図の構成
以外の仕方で送受信変換器を形成してもよい。変換器の
斜めの配置又は正常の配置が望ましい。例えば、頭部波
を発生させるためには、骨に対向する骨の同じ側に変換
器を隔だててもよい。従って、本発明は(互いに向い合
う平行に隔だてられた変換器の)、第4図の構成には限
定されない。 連続波信号がパルス信号の代りに投入される場合、信
号処理は、変更される。適切な信号処理は、例えば、引
用によってここに組込まれるR.B.アッシュマン等の「骨
腫の弾性特性測定のための連続波技法」第17巻,349−36
1頁(1984)に示されるように、当業者には既知であ
る。 本発明の第2実施例によれば、本発明者らは、低周波
振動と高周波振動とが異なった時点に到達することを利
用して、変換器の分離と関係なく、この分離の測定を必
要としない測定技法を見出した。前述したようにこの例
では、変換器を支持し位置決めしてこれらを患者に対し
てしっかりと保持するための何らかの機構はなお使用す
る必要はあるとしても、デジタルカリパス38は割愛しう
る。この目的のために同じカリパスを使用してもよい
が、そのデジタル出力は不要である。その他について
は、この装置は、第4図のものと同じである。従って、
この例と前記の例との主な差異は、ハードウェアにある
のではなく、受信変換器42に到来する信号の処理にあ
る。 この処理は、第9図に示されている。第1に、従来の
技法を用いて、受信信号のフーリエ変換を計算する(ス
テップ102)。次に、1MHzより高い全ての周波数を除き
(ステップ104)、2MHzより低い全ての周波数を除く
(ステップ106)2つの逆変換を行なう。これらの逆変
換は、第10A,10B図に112,114としてそれぞれ示した時間
域低周波信号と時間域高周波信号とを生ずる(これら2
つの信号は、受信信号の第1及び第2の受信された識別
可能な成分と呼ばれる)。 周波数信号と時間域低周波信号との周波数分離によっ
て、1つだけの広帯域パルスではなく、同時にか又は連
続的に送出される2以上の帯域制限パルスによって骨を
附勢することによって、均等の結果が得られる。即ち、
低周波レスポンスと高周波レスポンスとは、異なった時
又は同時に測定しうる。受信信号を「組織」成分「骨」
成分とに分離するために波(アナログ又はデジタル)
を用いてもよい。変換がさけられ、しかもほぼ均等の結
果が得られるので、これは、信号処理を簡略にする。こ
れらの代表的な処理の代りに、種々のスペクトルの推定
技法を用いてもよい。 時間領域高周波信号と時間領域低周波数との選定され
た特性間の関係を評価(ステップ108)して骨の状態の
指示を得る。この評価は種々の測定を含みうる。例え
ば、低周波パルス成分112と高周波パルス成分114との第
1到来時の比を取ると、骨粗鬆性のよい診断予知がえら
れると考えられる。第11図は、この比の平均と、骨粗鬆
症、正常及びスポーツマン(即ち、走者)の3つのカテ
ゴリーに予め分類された47人の婦人のサンプリングの分
布の標準偏差を表の形で示している。このデータへの周
知の「t分布」の適用は、98%の信頼度レベルまで、
「正常な患者」のグループと骨粗鬆症の患者のグループ
とが、前記の比によって識別されることを示す。このよ
うに、計算された比を1以上の所定の閾値と比較し、患
者の試験の結果を、その患者が異常な骨の状態をもつ対
応した可能性に組合されたカテゴリー中に入れることが
できる。いろいろの閾値(条件可能性レンズ範囲)は、
全部の患者について同一でなくてもよい。異なる年令,
性,人種又は他の特徴をもった患者についていろいろの
異なった閾値を使用しうる。 同様に、或る期間の間に患者について得た測定値の予
測子の変化率についての異なったレベルの信頼又は条件
可能性を設定しうる。即ち、1年間にx%の割合で予測
子が単調に変化する患者は、骨粗鬆症であること又はな
ることのy%の機会をもつと言える。適切な信頼度レベ
ル仕分け及びその閾値の設定は、多数の患者の試験と結
果の普通の統計的分析とを単に必要とする直截な訓練で
ある。 第1到来時間の比を用いる代りに、本発明者らは、受
信信号の出力スペクトルを励起の出力スペクトルで除算
した周波数依存利得関数を骨の状態の評価子として使用
することによって良い成果が得られることを見出した。
しかし、分母には零がありうるので、単純な除算を行な
う代りに、高周波成分による低周波成分のデコンボリュ
ーションを近似すると共にノイズ及びこれらの零を補償
する関数を励起法により形成することが望ましい。この
関数は、一例として、ウェイナー波によって得られ
る。又は、利得関数は、各々(例えば、1MHzより低い部
分)の低周波部分に制限することができる。これらの利
得関数の下方の領域とそのピーク値との両方は、骨の状
態の良好な指示子のように思われる。 音響信号の低周波成分と高周波成分とが異なった時点
において到達する事実を更に利用し、例えば(或るもの
は低周波成分をもち、他のものは高周波成分をもつ)、
一連のパルスを送出することによって、これらの高周波
成分と低周波成分とを別々に発生させることもできる。
この可能性に留意して、前記のように使用された単一の
「パルス」は、少くとも1つの低周波パルスと少くとも
1つの高周波パルスとの可能性を含むように解すべきで
ある。 以上の方法の記述は、限定的にでなく、例示的に解釈
すべきである。膝蓋骨は、梁骨の評価及び骨粗鬆症の診
断又はモニターのための人体の好ましい試験個所として
同定されている。脛骨は、コンパクトな骨の検査にとっ
て有用な結果を与える。骨粗鬆症の観察についての脛骨
の結果は、脛骨が梁骨でなく皮質骨(cortical bone)
から主に成っているため、或る程度制限される。そのた
め、骨粗鬆症の効果は、膝蓋骨及び脊髄のような梁骨に
おいてほど早くは、脛骨においては感じられない。骨粗
鬆症の診断に関係して本方法を以上に説明したが、本方
法は、他の骨の状態の診断にも使用しうる。これは、一
例として、この方法は歯周炎をひき起こす歯周骨の欠損
を診断し、骨折の治療を評価するためにも使用しうる。 また本発明は、動物の骨の状態にも適用されるが、適
用された可聴励起の周波数スペクトルは、動物に適合す
るように変更する必要があろう。また、骨粗鬆症以外の
状態を評価する場合に周波数スペクトルを変更すること
も時には望ましい。当業者に自明な他のいろいろの変更
又は改良が可能である。また本発明は、骨の侵襲試験に
も有用でありうる。この目的のために、変換器を先端に
取付けた針を、被検骨と直接に接触するように取付けて
もよい。もちろん、この場合、軟組織の調節が必要とな
る。 従って本発明は、前述の例ではなく、以下の請求の範
囲及び均等物のみによって制限される。
関する。より詳しくは、本発明は、ヒト及び動物の或る
骨の状態を非侵襲的に生体内において定める方法を開示
する。この情報は例えば骨粗鬆症(異常に低い骨密度)
又は骨化石症(異常に高い骨密度)のような骨の状態の
変化に関係した支障を主体がもつか否かを定めるために
使用しうる。 問題の解析と発明の目的 本発明は、骨の疾病及び他の骨の状態を検査すること
に関する。本発明は、特定的に、ヒトの評価及び診断に
予定されているが、例示的な実施例の明瞭な変更によっ
て、動物の骨の評価にも使用しうる。主な目的は、骨粗
鬆症の存在の診察である。しかし腎不全及び骨の劣化に
よる歯周炎の後の骨の状態の変化を含む他の状態も、こ
の方法を用いて検出しうる。本発明は骨折の検出と、骨
折の治ゆの定量的な評価にとっても有用である。 ヒトの骨格は、筒状骨(長くかたい骨)と、海綿状骨
とから成っており、その各々は、特定の比率の皮質骨
(コンパクトな骨)と梁骨(連結ストランド)とから成
っている。主に皮質(cortex)から成る筒状骨は、四肢
を形成する中垂性骨格を支配し、主に梁骨から成る海綿
状骨は、脊柱及び骨盤の軸骨格を支配する。 皮質骨においても、梁骨においても、コラーゲンのフ
イブリルが全体に延在している。両者の差異は、実際に
は、程度の差異であり、コラーゲンのフイブリルの網状
組織が取る形状に依存する。フイブリルは、一般に、網
状組織がロッドの網状組織であるように分離されていて
もよい。フイブリルはまたより密に広げられることによ
って、ロッドによって結合されたプレートの網状組織の
外観を呈してもよい。どちらの場合にも、骨の或る容積
は、鉱物化コラーゲンフイブリルから成り、或る容積
は、「骨髄」と呼ばれる流体から成る。 或る特別の骨が梁形成を取るか又は皮質形成を取るか
は、骨の果す機能に大きく依存する。骨の形成される仕
方は十分明らかではないが、骨のゆ着の過程は、骨が経
験する応力に或る意味で応答する。従って、比較的高い
応力を経験する骨の領域、例えば足の脛骨の骨幹は、皮
質骨に向う傾向を示す。低い応力を経験する骨の領域は
より梁骨に向う傾向を示す。梁骨の大部分の個所におい
て、梁物質は、皮質骨の比較的薄い層で囲まれている。
膝蓋骨は、多くは梁骨であり、皮質下層は、前面(ante
rior surface)の直下にある。 本発明の主な目標である骨粗鬆症は、一般に年を取る
につれて人を悩ます原因不明の疾病である。骨粗鬆症
は、男性よりも女性に多発し、女性については、月経閉
止の後に多発する。白人の女性は、他の人種の女性より
もこの疾病にかかり易い。骨粗鬆症は、骨組織物質の絶
対的な減少によって現出される。しかし残りの骨は正常
である。骨粗鬆症にかかった人は、皮質骨よりも骨梁な
いし梁骨をより高い比率で失なう。骨粗鬆症の普通の症
状は、せむしであり、これは脊髄の圧迫骨折、首の骨
折、上大たい骨及びラジアス(radius)の末梢端(肱
骨)の骨折によって生ずる。 骨粗鬆症は、その初期においては物理的変形が明らか
でないため、特に油断のならない状態である。骨粗鬆症
は徐々に進行するので、早い時期の診断と適正な治療と
は重大な症状を改善する。更に、骨が形成されるモード
のため、患者の比較的若い時期の増大した運動の習慣と
それに組合された適切なダイエットとは、この状態の効
果を最小にするとも考えられている。婦人の場合、月経
閉止前にこれらのプログラムを開始することが適切であ
る。骨粗鬆症の診断と治療とは、各々の患者が異なった
「正常な」骨密度をもつという事実によって複雑化す
る。そのため、個々の患者の骨の或る特性(即ち、密度
又は弾性)の変化の履歴を作成し、少くとも部分的にこ
の履歴の変化に基づいて診断することが有用である。 本発明の別の目標である歯周炎は、上顎骨と下顎骨の
損失と、それに伴なう歯のゆるみとをひき起こす。従来
は、顎の骨の欠損はX線検査によっていたが、X線検査
は、相当な骨の欠損が既に起こった後にした歯周炎の存
在を示すことができない。 そのため、骨の状態の変化(骨物質の欠損及びこれに
附随する骨の密度及び弾性の減少を含むが、これのみに
は限定されない)を検出する手段を用意することが望ま
しい。最適には、この手段は、非侵襲性、正確、高感度
で、使用し易く、入手が容易であるものとする。これは
一般には、従来の技術では実現できない。 骨粗鬆症の早期の診断のために、いくつかの方法が提
案されている。これらの方法には、全身体カルシウム中
性子活性化分析(TBC)、肱及びオスカルクス(oscalcu
s)の単一光子吸収測定、(脊柱及び大たい骨の首部
の)複式光子吸収測定、コンピューター支援断層撮影
(CATスキャンニング)並びに超音波分析法が含まれ
る。 TBC法は、要約すると、次の通りである。TBCには、約
14MeVのエネルギーの高速中性子源を必要としている。
中性子ビームは、被検体に向けられる。或る原子の核
は、入射中性子ビームのエネルギーに依存して、中性子
を捕そくし、励起される。励起された同位体は、直ちに
(10〜12秒のオーダー)か又は活性化生成物に特有の減
衰期間の後に、1以上のγ線を放出することによって、
安定な状態に戻る。放出されたγ線のエネルギーは、タ
ーゲットの元素を特徴付ける。カルシウムは、3.1MeVの
エネルギーにおいてγ線を放出する。γ線の強度即ち単
位容積当たりのカウントレートは、存在している問題の
元素の量を指示する。TBCの欠点は、全人体カルシウム
のみを好都合に測定するので、骨粗鬆症の初期に依存し
うる骨物質中の局在化された変化が識別できないことに
ある。またTBC装置は、原子炉、サイクロトロン又は放
射性物質のような中性子源を必要とする。また、大形の
複雑の機会が必要なため、比較的高価となり、数も少な
い。そのため多くの患者にとって利用可能ではない。ま
たこの方法は、期待されるように、比較的高価である。
また0.3〜0.5レムの比較的高い放射線量が必要となり、
反復検査ができない。 光子吸収測定は、エネルギーの異なった光子の物質に
よる吸収に依存する技術である。単一光子吸収測定にお
いては、単一エネルギー源が、単色の光子のビームを放
出する。伝送される光子の強度は、光子の走路上の物質
について指数関数的に減衰する。走行距離がわかると、
介在する物質の密度を定めることができる。単一エネル
ギー光子の問題は、骨を囲む軟組織のような多重減衰層
の効果を解消しえないことである。 そのため、複式光子吸収測定法が用いられる。この方
法によれば、2つの別々のエネルギー例えば153Gd又は
2つの別々のエネルギー源例えば125I及び241Amが用い
られる。各々2つの未知について2つの式が結果し、こ
れらが、2つの物質の密度に減衰されない光子の強度と
減衰された光子の強度とを関連付ける。2つの式を解い
て2つの密度を定めることができる。TBCの場合と同様
に、装置のコストが高く、数も少ない。また軟組織は、
皮膚の別々の成分である脂質と筋肉とから成っているの
で、「肥満/やせ」調整因子を使用しなければならな
い。 骨の密度は、CATスキャンを用いることによっても定
めることができる。CATスキャンナーは、X線伝送の多
重の入射から像を計算する装置である。患者は、全周の
回りにX線感知レシーバーを有する中空シリンダー中に
おかれる。X線源は、患者の回りに回転し、患者を通っ
て個別のX線ビームを送出し、このビームは、シリンダ
ーの直径の反対側のレシーバーによって受信される。こ
の伝送はシンリンダーの全周について反復される。受信
センサーが受けるX線の強度は、患者の身体を通る経路
の全密度に関係している。送信器が患者の周囲を移動す
る間にこれらの多数の読みが取られる。複雑なピュータ
ーによる分析によって、各々のレシーバーが受けた読み
が照合され、多くの伝送されるビームによって限定され
る平面と平行に取られた患者の身体の薄いスライスの密
度マップが作成される。 CATスキャンの出力は密度マップである。身体領域を
通るCATスキャンは、骨の密度に関係した情報を供与す
る。CATスキャンの実証された利点は、高再現性−高信
頼性の装置であり、他の関係した分野において十分に試
験されていることである。またそれは、梁骨と皮質の密
度を別々に量化する。そのため、海綿骨の正確な密度
を、周囲のコンパクトな皮質骨の存在による平均化効果
なしに定めることができる。 CATスキャン技法の欠点は、装置がコスト高で、多く
のコミュニティにおいては入手可能ではないことにあ
る。また、可搬型でなく、或る患者にとっては創傷性と
考えられている。CATスキャンの実施は、放射線学上の
有資格の医師又は放射線学者と協同する技師によらねば
ならないため、更にコスト高となる。またCATスキャン
においては患者がX線に露呈されるが、これは、或る医
師及び患者にとっては望ましくないと考えられている。
ここで「患者」とは、ヒト及び動物の被検体を意味す
る。 ヒト及び他の動物の骨を生体において非侵襲的に分析
するために、超音波の使用による測定も行なわれてい
る。しかしこの測定は、多くの骨を囲む軟組織の存在に
よって複雑になる。皮膚を通って投入された超音波信号
の速度は、骨を囲む軟組織の通過によっても、問題の骨
の通過によっても影響される。軟組織の効果を最小とす
ることも可能であるが、そのための従来の試みは、相当
複雑なエコー分析及び注意深い測定の要求によって妨げ
られている。例えば米国特許第4,361,154号「生体内で
骨の強度を定める方法」を参照されたい。 従って、本発明のいくつかの目的は、(1)非侵襲的
な生体内技法によって骨の状態を評定し、(2)実質的
なX線量に患者を露呈させることなく骨の状態を評価
し、(3)患者の骨の状態を経済的に迅速に評価し、
(4)数年間に亘って何回も経済的な安全に実行するこ
との可能な骨の状態の評価方法を提供し、(5)高価な
機器を必要とせずに、特に注意深い測定又は複雑な分析
なしに技師によって実行可能な骨の状態の評価方法を提
供することにある。 ここに使用されている限り、膝蓋骨とは、膝の前方部
分を形成する厚い扁平な3角形の骨であり、「ニーキャ
ップ」(膝蓋骨)は、膝蓋骨と周囲の軟組織であり、脛
骨は足の下部の前部骨であり、「すね」は脛骨及び周囲
の軟組織である。 発明の簡単な要約 本発明によれば、骨の状態は、骨の一部分を通って超
音波信号を送出し、受信される信号を分析することによ
って定められる。超音波信号は識別可能な波形又は周波
数の少くとも2つの成分を好ましくは含んでいる。これ
らの成分はパルス又は連続波のフォーマット又はその或
る組合せである。以下では、特に明記されない限り、
「パルス」という用語は、有限の持続時間の信号と連続
波信号との両方を意味する。狭帯域信号が使用される限
り、これらは、同時に(例えば複合波形の一部として)
又は異なる時間に伝送することができる。別の方法とし
て、広帯域信号を用いてもよい。これは広い周波数帯域
例えば約100kHzから約3MHzのエネルギーを含む信号でも
よい。2つ(又は以上)の成分の伝播の差は、骨の状態
についての情報を含んでいる。 本発明は、媒体を通る超音波信号の伝送速度が媒体の
密度及び弾性に関係しているという原理に基づいて作用
すると考えられる。本発明はまた減衰が周波数依存性で
あることの知見に依存している。皮質骨から梁骨までの
骨の可変の性質のため、超音波信号の速度と媒体の密度
及び/又は弾性との間の関係を正確に表現することはで
きない。しかし、大部分が皮質である骨と大部分が梁骨
である骨について行なった実験の結果とよく相関される
皮質骨と梁骨とについてのモデルが提案されている。即
ち超音波伝播の相対速度は、骨の破損し易さの指示とし
て使用することができる。パルスの送受信のために1対
の変換器が使用される。これらの変換器は、超音波パル
ス発生器と信号受信処理システムとに接続される。これ
らの変換器は、骨の近くの2つの異なった個所に配され
ることにより、超音波信号が骨及び軟組織中に送出され
てこれを通過し、受信変換器によってさえぎられるよう
に配置される。 測定しようとする骨は、好ましくは、最小の軟組織で
囲まれ、その表面は、好ましくは、向い合う部分を備え
ているので、音のエネルギーは、実質的に平行に配され
た送受信変換器の表面によって、これらの部分に入出す
ることができる。従って、超音波信号の第1到来部分の
音の経路は、送信変換器の表面において開始し、骨を覆
う軟組織を通り、骨を通り、次に、受信変換器から骨を
隔てる軟組織を通り、最終的に、受信変換器に到達する
(超音波信号の或る部分は、音響的により長い別の経路
上を伝播し、従って、より遅い時間に到達するであろ
う。)。 いろいろの骨部位を試験の目的に使用してよいが、
(a)膝蓋骨を囲む軟組織層を通る距離が骨を通る距離
よりも著しく短かく、また(b)膝蓋骨及び軟組織層を
通る伝播経路とほぼ平行でしかもこれとほぼ等距離の膝
蓋骨の上方の軟組織を通る音響伝播経路があるため、ニ
ーキャップが有利な個所であることが見出された。その
他の有用な部位の例は、頚骨、上顎骨及び下顎骨があ
る。これらの2つの後者の骨の検査は、歯周炎の発見の
道を開くであろう。 本発明の第1実施例によれば、受信信号は、ノズルを
去するように処理されることにより、超音波伝播遅れ
の持続時間をより正確に定めることを可能とする。変換
器の送信面と受信面との距離は、カリパスによって測定
され、骨を通る速度が計算される。この速度の計算に際
して、骨を覆う軟組織の有限な厚みによって導入される
誤差について測定伝播遅延を最初に修正することが望ま
しい。 骨及び軟組織を通って測定された速度、又は、軟組織
について修正された速度は、種々のデータと比較しう
る。第1に、試験中の患者の時間履歴を作成し、骨の密
度及び/又は弾性の変化を表わす速度の変化をチェック
することができる。この形式の比較は、疾病の経過又は
患者のリハビリをチャートする。また、全ての関係ある
年令の健康な個人及び疾病をもった個人を含めた公衆の
広範な断面について、骨の密度、弾性及び/又は速度を
含むデータベースを作成することができる。選定された
閾値基準に基づいて、このデータベースの情報に、患者
の骨の密度、弾性及び/又は速度の測定値を比較し、患
者の骨が健康か否かを臨床的に有用な信頼性レベルで予
測し、又は、患者が骨の疾病をもっている可能性を評価
することができる。 本発明の第2実施例によれば、送受信変換器の間の距
離の測定は必要ではない。送出変換器によって送出され
る超音波信号は約100kHzから約3MHzまでの帯域内の、少
くとも2つの別々の周波数成分をもっている。患者の軟
組織と骨とを通って伝播された後、この信号の低周波成
分が最初に受信変換器に到達し、おそく到達する高周波
成分が後続する。これら2つの成分の到来時間の比が骨
の状態の好い予測子であることが見出されている。受信
された低周波(即ち1MHz以下)の成分の出力スペクトラ
ムの、同じ帯域内の適用された励起に対する比のよう
な、他の測値も、良い予測子のように思われる。この比
は「利得関数」と呼ばれる。ピーク振幅及びこのような
利得関数の下方の面積も、予測子として使用可能であ
る。更に、これらの種々の予測子は、更に高信頼度の予
測子を実現するために、タンデムに用いてもよい。 本発明の前述した目的並びに他の記述されない目的、
利点及び特徴は、以下の詳細な説明及びここに添付され
た請求の範囲を添付図面との関係において読むことによ
って一層明らかとなるであろう。 図面の簡単な説明 図において、 第1図は、本発明の方法の各ステップを一般的に示し
たフロー図である。 第2図は、第4図の装置に用いられている変換器から
送出される超音波信号のスペクトラムを示す。 第3図は、本発明の方法の第1実施例によるより特定
的なステップの連鎖として、第1図のステップの構成を
より詳細に示したフロー図である。 第4図は、本発明の2つの実施例の方法を実施するた
めに使用可能な装置の略図である(ただし第2実施例に
ついてはカリパスは不要である)。 第5図は、1つの超音波パルスについて受信励起の振
幅対時間の典型的なプロティングを示した略図である。 第6図は、「整合波」の方法を略示する。 第7図は、健康な骨を通過した後の第1到来波形を略
示する。 第8図は、骨粗鬆症の骨を通過した後の第1到来波形
を略示する。 第9図は、本発明の方法の第2実施例によるより特定
的なステップの連鎖として、第1図のステップ7の構成
をより具体的に示すフロー図である。 第10A,10B図は、第9図の方法によって形成された時
間域高周波波形及び時間域低周波波形を略示する。 第11図は、第9図の方法によって得た例示的データを
示す表である。 実施例の詳細な説明 第1図には、本方法の基本的な工程が、概略的なフロ
ー図によって図示されている。符号2によって示した第
1工程は、超音波信号の発生、管理及び測定を実行する
ための適切な装置を供与することを要求する。第4図に
ついて後に一層詳細に説明するこの装置は、一般に、超
音波(例えばパルス)発生器、1対の変換器及び信号検
出−分析ユニットから成っている。 広帯域波形発生器は、約100kHzから約3MHzのエネルギ
ー成分を含む音のパルスを約2×10-6秒の間前記変換器
のうちの1つから発生させるために、ブロック4に示す
ように設けられている。代表的な装置(後述する)から
の励振のスペクトルは、第2図に示されている。このス
ペクトル全体に亘って信号の振幅及び出力が変動可能な
ことは注意されたい(この理由により、約250kHz及び2.
5MHzの別々の狭帯域信号を、前述したように、別の方法
として用いてもよい)。パルスの発生は、データ取得シ
ステムによって記録され、このシステムは、磁気ディス
ク又は他の記録保存可能な媒体上にそれを記憶させる。
これを陰極線管又は他の適当な表示装置上に表示しても
よい。次のステップ6は、受信変換器においての音響励
振のモニタリングを含む。これらの励振は、背景ノイズ
と、骨及び軟組織を経て伝送される第1到来信号と、骨
及び軟組織を経て伝送される第2到来信号とを含む。ス
テップ6は、音響励振に対応する電子信号を発生させ、
データ取得システムによるこれらの励振信号の受信の記
録を記憶させ、同一の陰極線管表示部上にこれらの信号
の受信の記録を表示することを含む。ステップ7では、
ステップ4及び6から得た受信信号は(アナログ又はデ
ジタル波及び種々の計算を含みうる各種の技法によっ
て)、処理され、2以上の誘導された測定を発生させ
る。これらの誘導された測定は、次に、患者を類別化す
る因子、例えば、限定的にではないが、年令,性,人種
及び患者の骨の状態に対して受信信号特性を相関させる
情報データベースと比較される。 ステップ4〜6を数回反復してその結果を平均するこ
とによって、ランダムなノズルによる誤差を減少させる
ことが好ましいが、これは、不可欠ではない。 患者の軟組織及び骨は、部分的には、その伝送特性が
周波数の関数として変化する複雑なネットワークのよう
に作用する。患者を通過する時の超音波入力励振の変更
は、受信側の変換器において、大きく変更されたスペク
トルを供与する。このスペクトルは、異なる時において
到来する少くとも2つの別々のパルス又は信号を含んで
いる。原初の伝送される超音波信号に対して比較される
と共に相互に対しても比較されるこれらのパルスの特性
の、アナログ処理手段又はデジタル処理手段による評価
によって、患者の骨の状態についての情報が得られる。
評価すべき特性には、限定的にではないが、周波数、振
幅及び位相が含まれる。 患者を通って伝搬される時に生ずる超音波励振の複雑
な変化は、未だ十分には理解されていない。現在のとこ
ろ、膝蓋骨については、問題の周波数スペクトルにおい
て、比較的高い周波数成分は、骨を覆う軟組織を経て漏
れモードとして通過し、比較的低い周波数成分は、比較
的密な骨の中の並行な経路に沿って伝搬される。より詳
細には、膝蓋骨については、比較的低い周波数信号は、
平面の直下の比較的密な皮質下骨層中を伝搬されるもの
と思われる。この骨の経路は、加えられた超音波パルス
の高周波成分を著しく減衰させる。他の骨についての変
化は、それほどよく理解されていない。 本発明の第1実施例(第3,4図参照)によれば、変換
器は、デジタルカリパス装置38に取付けられている。カ
リパスは、評価しようとする骨の部分の回りに、変換器
が患者の周囲の皮膚と接触した状態において閉ざされ
る。変換器と皮膚との間の超音波エネルギーの結合は、
インピーダンス整合ゲル又は他の同様の材料によって容
易にされる。励振が発信される受信される時に、キャリ
パーは、ブロック8によって示すように、パルスの通っ
た距離(即ち、変換器の向い合う表面の間の距離)を自
動的に測定する。伝送データ及び励振データが記録され
た後、テンプレート信号が、10によって示したように供
給される。この信号は、受信側の変換器においての励振
に対応する信号に、相互相関機能によって整合される。
相互相関機能の評価によって、ベースライン信号を形成
するノイズの波が可能となるので、パルスの正確な到
来時間をピンポイントすることが容易となる(ブロック
12)。今や、12からのパルス走行の経過時間と、8から
の走行距離との両方が知られたので、骨及び最小限度の
周囲の軟組織を通る超音波パルスの平均速度を計算する
ことができる(ブロック14)。 通常は骨の多い個所である膝蓋骨のところにも、患者
によっては、少くない量の軟かい組織の存在することが
あるので、軟組織の効果を除くために、速度測定を修正
したり調節したりすることが有用となる。本発明に従っ
て、患者の身長−体重比に関係した因子を使用すること
ができる。 速度の測定は、別の方法として、前出の米国特許第43
61154号(引照によって本明細書の一部分となる)によ
り詳細に記載されている。より複雑なエコー技術に従っ
て、患者の軟組織について調節したり修正したりするこ
とができる。この方法は、軟組織と骨との界面を信号が
横切る際に超音波エネルギーの一部分が骨から反射され
るという事実を利用している。各々の変換器から超音波
パルスを発信し、反射エコーパルスを受信するまでの時
間を測定することによって、超音波信号経路の各々の先
端において骨を覆う軟組織の厚みを測定することができ
る。この技法には、両方の変換器から超音波パルスを発
信して両方の変換器においての励振をモニターすること
のできる装置が必要となる。軟組織を通る超音波の速度
は既知であり、約1550m/秒である。従って、超音波エコ
ー発信の記録のチェックによって、超音波パルスが軟組
織を通過し、骨から反射され、発信側変換器に戻るまで
の時間が、各々の変換器について得られる。この走行の
間の速度を知ることによって、経路の長さも定めること
ができる。軟組織を通るこの経路の長さを変換器の間の
全距離から引算すると、骨を通った距離が求められる。
また骨を通る伝搬時間も、前記のようにして定めた軟組
織を通る走行に要した時間を測定された時間から引算す
ることによって計算される。 2つの変換器の間において測定される距離も、標準的
な「皮膚フォルド」試験によって得た因子によって、患
者の軟組織特性に関連して正規化することができる。皮
膚(スキン)フォルド試験によれば、患者のいくつかの
指定された個所において、これらの個所の軟組織の皮膚
量を、カリパスを用いて測定する。一般大衆から多量の
データを収集し、皮膚フォルド試験に基づいて、一般的
な「肉付き」に従って、人々をランク付けした。この試
験に基づいた指標は、骨及び軟組織を通って測定された
速度を正規化するためにも使用することができた。 身長/体重比、皮膚フォルド基準、エコー技法など、
どの方法によって得たものであっても、軟組織について
の調整は、以下に、「軟組織調整因子」と称される。即
ち、測定された絶対速度と、「軟組織調整因子」とか
ら、軟組織調整速度が、ステップ18において定められ
る。 軟組織調整速度は、多くの仕方で使用でき、第3図に
は、そのうちの2つの方法が示されている。第1に、調
整速度は、広い階層の人々から集めた調整速度の編集デ
ータのエントリーと比較される(22)。このデータバン
クは、年令その他の患者の特徴(例えば、性,人種,種
族的集団など)及び既知の骨格の状態に従って組織化す
ることができる。即ち、40才の婦人の調整速度が大衆中
の40才の婦人の既知速度の範囲に含まれれば、この人は
健康な骨をもっていると言えよう。その逆に、調整され
た速度が、健康な60才の婦人に通常組合される範囲に入
っている場合は、異常な状態である確率は実質的に高
く、予防上及びリハビリテーション上の処置を開始する
ことが賢明であろう。 調整速度は、検査中の患者について行なった以前の測
定及び将来の測定と比較することもできる(ステップ2
6)。この計画の下では、何年かに亘って1人の個人に
ついて測定が行なわれよう(ステップ24)。超音波速度
の際立った変化は、骨の状態の際立った変化を表わすも
のとなろう。骨の状態の変化の方向に依存して、患者の
健康が回復途上にあるか又は悪化しているかを定めるこ
とができよう(ステップ28)。 次に、第4図を参照して、本方法を実施するのに適し
た装置について詳細に説明する。各々の変換器40,42
は、広帯域高周波数超音波パルスを発信することができ
る。02254マサチューセッツ州ワルサム,クレセント通
り221,パナメトリックス社によって製造された変換器例
えばパナメトリックスモデルA533Sトランスジューサー
を使用して、良い結果が得られた。これは2.25MHz変換
器として表わされる。2.25MHzは、いくつかの理由によ
って、望ましい基本周波数である。第1に、この周波数
において超音波を送受信しうる変換器は、比較的小形で
あり、直径が2.54cm(1インチ)の約1/3である。その
ため、これらの変換器は、患者の身体のいろいろの部分
に好都合に係合させることができる。第2に、周波数
は、少くとも1×10-6秒の短い時間の測定値を解像しう
る装置の感度が得られるように、十分に高い値としなけ
ればならない。これは、少くとも1MHzの周波数を必要と
する。比較的高周波数の装置は、比較的低周波数の装置
に比べて、時間測定値の解像能力がより高くなる。不幸
なことに、周波数が高すぎると、受信側変換器において
受信される信号の振幅は、骨の内部又は周囲の境界面に
おいてその反射によって相当に減衰する。本発明者ら
は、約2.5MHzの基本周波数が、感度と発信振幅との競合
する配慮の間の妥当な調和点を与えることを見出した。
第2図を参照して、前述した変換器によって送出される
信号のスペクトルが示されている。理解されるように、
このスペクトルは、約100kHzから約3MHzまでの実質的な
エネルギー成分を含んでいる。この広汎なスペクトルが
望ましいのは、骨を経て伝送される信号が、300kHz付近
の周波数において際立ったエネルギーをもち、軟組織を
経て伝送される信号が、2.5MHz付近の周波数において際
立ったエネルギーをもつためであり、そのため、変換器
−信号発生器52は、両方の周波数域において動作可能と
すべきである。これら2つの伝送(即ち、骨及び軟組
織)は、実質的に独立しているので、別々の周波数にお
いて動作する複数の狭帯域ソース又は単一の広帯域ソー
スによって励振された狭帯域変換器によって、同様の結
果が得られる。 変換器(複数)は、02166マサチューセッツ州ニュー
トン,660−・ストリート,フレッド・V・フロウラー社
によって、モデル名“ウルトラ−カル”の下に市販され
ている自動電子デジタルバーニアカリパス38上に取付け
られる。バーニアカリパスは、固定のアーム44と可動の
アーム46とを備えている。可動のアーム46が支持アーム
48に沿って移動すると、変換器面44′,46′の間の距離
が表示部50上に表示される。このカリパスの構造は、開
示された発明にとって大切ではなく、普通のバーニアカ
リパスを用いてもよい。しかしその場合には、距離の値
を手動で読出して記録すると共に、この距離をデータ処
理装置に手動で読込むことが必要となろう。 使用に当り、送信側の変換器40と受信側の変換器42と
は、患者の身体の骨部分60、好ましくは膝蓋骨の両側の
表面上に当接させる。パルス発生器52は、持続時間が約
2×10-6秒のパルスを発生させる。このパルスは、電気
接続部54を経て、発信側変換器40に転送される。それと
同時に、電気接続部56に沿って発信信号がデータ取得シ
ステム58に転送され、発信信号の送出を指示する。この
信号は、送出時を記録するために、信号処理装置によっ
て使用される。データ取得システムについては、マサチ
ューセッツ州ピーボディ,アナロジックス・コーポレイ
ションのデータ・プレシジョン・ディビジョンによって
製造された「データ6000」を使用して、良好な結果が得
られた。 パルスは、測定される部分60を通過し、受信側の変換
器42によって受信される。受信側の変換器42は、導線62
を経てデータ取得システム59に電気的に転送される信号
を送出する。これらの信号は、変換器42によって受信さ
れる超音波パルスに対応している。 データ取得システム58は、それが受信する信号の振幅
を記憶すると共に、標準型のCRT表示部64上に、時間変
化振幅を表示することができる。データ取得システム
は、第4図に略示したように、データ計算手段66を備え
ている。カリパス38によって測定した、変換器面44′,4
6′の間の距離を表わす信号は、後述するように、後の
ステップにおいて使用するために、電気接続部68を経
て、コンピューター66に転送されている。 次の図である第5図には、後述する「整合された
波」の相互相関技法が適用される前の受信側変換器が受
けた典型的な応答が略示されている。第5図において、
ディスプレイは、64により示されている。水平軸線は、
経過した時間を表わしている。原点から楕円644の終点
までの時間は、15−20×10-6秒のオーダーである。垂直
軸線は、受信信号の振幅を表わしている。 円642中に示された信号の部分は、「メイン・バン
グ」と称され、送出パルスの記録を形成する。大体にお
いて楕円644に囲まれた信号部分は、受信側の変換器に
よって受信されたベースラインノイズである。楕円646
によって示された信号部分は、送出パルスの結果として
骨及び軟組織を通過する信号が受信側変換器によって最
初に受信されたことを表わしている。この受信された信
号は、2つの識別可能な成分から成っている。これら
は、最初に到来する低周波信号と、これに続く、後に到
来する高周波信号である(これらは、それぞれ「受信さ
れた第1の識別可能な成分」及び「受信された第2の識
別可能な成分」と称することができる。或る数(例えば
16以上)の別々の初パルスから受信信号の振幅を平均化
することが特に好ましい。ランダムな信号例えばノイズ
は平均化によって零になる傾向をもっている。 第5の略図によって明らかなように、信号646がベー
スラインノイズ644から出る時間を正確に定めることは
困難である。到来時間を定める際の誤差が、1×10-6秒
のように大きいと、速度の計算の誤差(誤り)が10%と
非常に大きくなることがある。この誤差量による測定
は、有用であっても、この誤差を可能な限り除去する
(又は少くとも最小とする)ことがもちろん望ましい。 幸いに、誤り訂正法が開発されている。これは、例え
ばワイリー社から1960年に出版された、Y.W.リー著「通
信の統計理論」のような、種々の刊行物に、十分に記載
されている。誤り訂正方法は、第6図に略示されてい
る。一般に、ベースラインノイズは、信号の平均化及び
相互相関によって、受信信号の表示から去される。骨
を経て受信される平均化された第1到来信号の形状に相
当な程度まで整合されたテンプレート信号が供給され
る。以下において明らかとなるように、骨を通る受信信
号の前縁にテンプレートが正確に同一であることは必要
ではなく、単に類似していることが必要である。テンプ
レートは、その患者について受信された支配的なフーリ
エ成分の適切なサイクル数を取ることによって、例えば
取得される。好ましくは、2つのテンプレート信号が使
用され、その1つは、骨を通る伝送に所属し、第2のも
のは、軟組織を通る伝送に所属する。後者の信号は、標
準的な信号でありうるが、その理由は、軟組織を通る超
音波の速度と軟組織の周波数応答とが患者ごとにほとん
ど相違しないためである。テンプレート内の或る適切な
数のサンプル点を与えるために、2つの異なったテンプ
レート信号について2つの別々のサンプリング間隔及び
/又はサンプリング率を用いることが時に望ましい。 信号の或るわずかな部分、例えば1サイクルの180゜
分が保持される。相互相関関数C(d)は、 として定義される。 関数F(t)は、第6図に示すように、受信側変換器
によって受信された信号である。関数G(t−d)は、
原点の付近においてのみ非零のテンプレート信号であ
る。骨を通る信号の到達時点をピンポイントするため
に、時間が零から無限に及ぶ積分を行なう。次に“d"が
零から無限大に及ぶ間のこの積分の第1極値を選定す
る。 C(d)の絶対値が、人為信号に特有の或る閾値を超
える第1極大値を超過した時点648は、後述する場合を
除いては、信号が最初に到達した時である。これは次の
ことによって生ずる。伝送パルスを最初に受信するまで
は、関数F(t)は、メイン・バングとベースライン64
4との急速に振動する信号から成っている。これらの信
号は一般に約2MHzにおいて振動する。その逆に、骨梁に
とって適切なテンプレート信号は、3.3×10-6秒の周期
をもち、従って周波数は300kHzである。主要なバング及
びベースラインのノイズは、正から負に変化するので、
2つの関数F(t)及びG(t−d)の積の積分は、高
速振動の全ての時において零となる。テンプレートG
(t−d)及びF(t)(今や伝送パルスからの信号か
ら主に成っている)である相互相関関数の積が非零値を
もつのは、伝送パルスが到達(時点648)するまでは起
こらない。受信信号がテンプレートに最も近い値となる
時に、この積は、極大値となる。メイン・バングが第2
図の広帯域波形をもつ代りに複雑の別々の信号から成る
場合に対応した状態となる。 前述したように、相互相関関数の極大値が第1到来パ
ルスの到来と合致しないいくつかの場合がある。これら
は、第1到来パルスの振幅がテンプレートに対しての逆
の符号になっている場合である。例えば、第1到来信号
が負の場合に骨信号の正の区分からテンプレートが選ば
れている場合である。これらの場合は、相互相関関数が
最大値の前に比較的大きな負のエクスカーションをもつ
ことを認識することによって容易に識別できる。この場
合には、第1到来パルスの到来は、大きな負のエクスカ
ーションにおいて表わされる。負のテンプレートをもっ
たものと正のテンプレートをもったものとの2つの相互
相関関数を実行し、観察された結果に基づいた適切な保
持すべき記録を選択することによって、より正確な測定
を行ないうる。また、相互相関関数を平方してもよく、
この場合には、結果関数の極大値は、第1到来パルスの
到来と合致する。テンプレートG(t−d)をC(t)
と組合せてC′(d)を得るように、整合波のプロセ
スを1回以上実行することも時には望ましい。 この信号整合形態は、前述したデータ取得方式によっ
て自動的に実現することができる。 この整合波は、患者の診察の時には行なわずに、記
憶されたデータについて後に行なってもよい。相互相関
関数が極大にされ、信号伝播の経過時間が定められた
後、以前に測定された変換器の間の距離を経過時間で割
算することによって、速度を計算することができる。 変換器間の経路を通る平均速度を定めたら、骨を囲む
軟組織を通る超音波パルスの通過の効果を除くことが適
切である。それは、1つの変換器から次の変換器へのパ
ルスの通過時間には、超音波パルスが骨を通過している
時間と、超音波パルスが骨を囲む軟組織を通過している
時間とが含まれるからである。更に、超音波パルスの走
行距離には、骨を通る距離と、骨の両側の軟組織を通る
距離とが含まれる。軟組織を通るこの通過遅延を勘案し
えないと、単に走行距離を時間で割算することによる速
度の計算の精度が低下する。 骨梁の状態を検査し、骨粗鬆症についてテストするた
めには、膝蓋骨が、人体の中で非常に良い場合であるこ
とが、実験によって示された。膝蓋骨は、軟組織を実質
的に含まない人体の部分である。また膝蓋骨は、ほとん
ど梁骨(骨粗鬆症の効果が最も早く且つ最も深刻に感じ
られる形態の骨)から成っている。本発明者らは、脛骨
稜は、骨腫の状態を調べるための人体の中の良好な場所
であると考える。 前述したように、膝蓋骨及び脛骨稜において測定する
ことに利点があっても、これらの個所においての多量の
軟組織の存在は、結果に多少とも影響する。従って、患
者の特定の軟組織の状態について調整することによっ
て、速度の測定値を修正することが適切である。前述し
たように、これらの方法には、エコー技法、スキンフォ
ルド試験及び身長/体重の比較が含まれる。どの方法に
よって計算した速度も、「調節された骨速度」と呼ばれ
る。 調節された骨速度が1度定められたら、この値を用い
て、患者の骨格状態を評価することができる。これは、
広い範囲の公衆から得た調節された骨速度に前記の調節
された骨速度を比較することによって行なうことができ
る。理想的には、調節された骨速度の一般的なデータベ
ースは年令及び他の特性の範囲について分類し、或る年
令及び特徴の健康な人の典型値が知られるようにする。
このデータベースは、通常よりも多孔質なものから通常
よりも密なものまでの、骨密度の状態のスペクトルに対
応する一連の速度を、各々の年令又は年令群について含
んでいてもよい。患者の調節された骨の速度を彼らの年
令に対する調節された骨速度のスペクトラムと比較する
ことによって、患者の相対的な骨の状態を定めることが
できる。 また、或る個人が、骨粗鬆症の疑いがあったり、その
治療を受けている場合、同じ患者の以前の測定値に、特
定の測定値を比較することができる。これらの測定値
は、系列的な測定として知られる。これにより、その進
行又は劣化を図式化し、また治療することができる。 伝送された波形の振幅の形状に基づいて患者の骨の状
態を評価することも、本発明の方法の1つの特徴であ
る。受信された超音波パルスの波形の形状は、患者の骨
格の状態及び密度についての貴重な情報を含んでいる。
本発明者らは、健康な患者の場合に、受信変換器に最初
に到達する膝蓋骨を通る信号の形状が、約300kHzの周波
数のrmsベースラインノイズの2倍以上の、形状の定ま
った振動の第2部分と、2.5MHzに近い非常に高振幅−高
周波数の振動とによってほぼ特徴付けられることを見出
した。第7図に骨信号のこの形状を示す。第7図におい
て、健康な骨を通過した信号は、“70"により示され
る。第8図において、病気の骨を通過した信号(「病信
号」は、“80"により示される。この図からわかるよう
に、病信号は、高周波数−高振幅において開始する、ほ
ぼ形状の定まっていない振動において到来する。 波形の形状のこの変化に内在する機構と波形の意味と
は、十分には理解されていない。形状の定まった低周波
信号は、縦波又は圧縮波を示し、高周波−高振幅振動は
横波又はせん断波である。縦波は、伝播方向と平行な方
向の媒体の振動から結果し、せん断波は、伝播方向と直
角の振動の実証である。縦波は一般に横波よりも高速で
ある。しかし骨粗鬆症においては、梁骨部分を通る圧縮
超音波パルスの走行路は、より多く制限されており、従
って、梁骨形成物の終端と介在する流体ポケットとによ
って中断されると考えられる。媒体においてのこれらの
不連続は、骨の弾力を変更させ、超音波パルスに合体さ
れたエネルギーの放散を生ずるので、骨粗鬆症の患者の
場合、圧縮超音波パルスの振幅は減少する。 また、変換器が頭部波及びせん断波を開始(即ち送
出)するようにしてもよい。これらの波とそれらが誘起
する別の波との伝播速度は、変換器の間の骨の状態及び
性状についての別の情報も供与する。或る骨は多くの異
なったモードで励起されうる。関心のある骨の特性(例
えば、多孔性、密度及び弾性)について特に敏感な伝送
のモードを見出すことが目標である。 これらの他の波形式を発生するために、第4図の構成
以外の仕方で送受信変換器を形成してもよい。変換器の
斜めの配置又は正常の配置が望ましい。例えば、頭部波
を発生させるためには、骨に対向する骨の同じ側に変換
器を隔だててもよい。従って、本発明は(互いに向い合
う平行に隔だてられた変換器の)、第4図の構成には限
定されない。 連続波信号がパルス信号の代りに投入される場合、信
号処理は、変更される。適切な信号処理は、例えば、引
用によってここに組込まれるR.B.アッシュマン等の「骨
腫の弾性特性測定のための連続波技法」第17巻,349−36
1頁(1984)に示されるように、当業者には既知であ
る。 本発明の第2実施例によれば、本発明者らは、低周波
振動と高周波振動とが異なった時点に到達することを利
用して、変換器の分離と関係なく、この分離の測定を必
要としない測定技法を見出した。前述したようにこの例
では、変換器を支持し位置決めしてこれらを患者に対し
てしっかりと保持するための何らかの機構はなお使用す
る必要はあるとしても、デジタルカリパス38は割愛しう
る。この目的のために同じカリパスを使用してもよい
が、そのデジタル出力は不要である。その他について
は、この装置は、第4図のものと同じである。従って、
この例と前記の例との主な差異は、ハードウェアにある
のではなく、受信変換器42に到来する信号の処理にあ
る。 この処理は、第9図に示されている。第1に、従来の
技法を用いて、受信信号のフーリエ変換を計算する(ス
テップ102)。次に、1MHzより高い全ての周波数を除き
(ステップ104)、2MHzより低い全ての周波数を除く
(ステップ106)2つの逆変換を行なう。これらの逆変
換は、第10A,10B図に112,114としてそれぞれ示した時間
域低周波信号と時間域高周波信号とを生ずる(これら2
つの信号は、受信信号の第1及び第2の受信された識別
可能な成分と呼ばれる)。 周波数信号と時間域低周波信号との周波数分離によっ
て、1つだけの広帯域パルスではなく、同時にか又は連
続的に送出される2以上の帯域制限パルスによって骨を
附勢することによって、均等の結果が得られる。即ち、
低周波レスポンスと高周波レスポンスとは、異なった時
又は同時に測定しうる。受信信号を「組織」成分「骨」
成分とに分離するために波(アナログ又はデジタル)
を用いてもよい。変換がさけられ、しかもほぼ均等の結
果が得られるので、これは、信号処理を簡略にする。こ
れらの代表的な処理の代りに、種々のスペクトルの推定
技法を用いてもよい。 時間領域高周波信号と時間領域低周波数との選定され
た特性間の関係を評価(ステップ108)して骨の状態の
指示を得る。この評価は種々の測定を含みうる。例え
ば、低周波パルス成分112と高周波パルス成分114との第
1到来時の比を取ると、骨粗鬆性のよい診断予知がえら
れると考えられる。第11図は、この比の平均と、骨粗鬆
症、正常及びスポーツマン(即ち、走者)の3つのカテ
ゴリーに予め分類された47人の婦人のサンプリングの分
布の標準偏差を表の形で示している。このデータへの周
知の「t分布」の適用は、98%の信頼度レベルまで、
「正常な患者」のグループと骨粗鬆症の患者のグループ
とが、前記の比によって識別されることを示す。このよ
うに、計算された比を1以上の所定の閾値と比較し、患
者の試験の結果を、その患者が異常な骨の状態をもつ対
応した可能性に組合されたカテゴリー中に入れることが
できる。いろいろの閾値(条件可能性レンズ範囲)は、
全部の患者について同一でなくてもよい。異なる年令,
性,人種又は他の特徴をもった患者についていろいろの
異なった閾値を使用しうる。 同様に、或る期間の間に患者について得た測定値の予
測子の変化率についての異なったレベルの信頼又は条件
可能性を設定しうる。即ち、1年間にx%の割合で予測
子が単調に変化する患者は、骨粗鬆症であること又はな
ることのy%の機会をもつと言える。適切な信頼度レベ
ル仕分け及びその閾値の設定は、多数の患者の試験と結
果の普通の統計的分析とを単に必要とする直截な訓練で
ある。 第1到来時間の比を用いる代りに、本発明者らは、受
信信号の出力スペクトルを励起の出力スペクトルで除算
した周波数依存利得関数を骨の状態の評価子として使用
することによって良い成果が得られることを見出した。
しかし、分母には零がありうるので、単純な除算を行な
う代りに、高周波成分による低周波成分のデコンボリュ
ーションを近似すると共にノイズ及びこれらの零を補償
する関数を励起法により形成することが望ましい。この
関数は、一例として、ウェイナー波によって得られ
る。又は、利得関数は、各々(例えば、1MHzより低い部
分)の低周波部分に制限することができる。これらの利
得関数の下方の領域とそのピーク値との両方は、骨の状
態の良好な指示子のように思われる。 音響信号の低周波成分と高周波成分とが異なった時点
において到達する事実を更に利用し、例えば(或るもの
は低周波成分をもち、他のものは高周波成分をもつ)、
一連のパルスを送出することによって、これらの高周波
成分と低周波成分とを別々に発生させることもできる。
この可能性に留意して、前記のように使用された単一の
「パルス」は、少くとも1つの低周波パルスと少くとも
1つの高周波パルスとの可能性を含むように解すべきで
ある。 以上の方法の記述は、限定的にでなく、例示的に解釈
すべきである。膝蓋骨は、梁骨の評価及び骨粗鬆症の診
断又はモニターのための人体の好ましい試験個所として
同定されている。脛骨は、コンパクトな骨の検査にとっ
て有用な結果を与える。骨粗鬆症の観察についての脛骨
の結果は、脛骨が梁骨でなく皮質骨(cortical bone)
から主に成っているため、或る程度制限される。そのた
め、骨粗鬆症の効果は、膝蓋骨及び脊髄のような梁骨に
おいてほど早くは、脛骨においては感じられない。骨粗
鬆症の診断に関係して本方法を以上に説明したが、本方
法は、他の骨の状態の診断にも使用しうる。これは、一
例として、この方法は歯周炎をひき起こす歯周骨の欠損
を診断し、骨折の治療を評価するためにも使用しうる。 また本発明は、動物の骨の状態にも適用されるが、適
用された可聴励起の周波数スペクトルは、動物に適合す
るように変更する必要があろう。また、骨粗鬆症以外の
状態を評価する場合に周波数スペクトルを変更すること
も時には望ましい。当業者に自明な他のいろいろの変更
又は改良が可能である。また本発明は、骨の侵襲試験に
も有用でありうる。この目的のために、変換器を先端に
取付けた針を、被検骨と直接に接触するように取付けて
もよい。もちろん、この場合、軟組織の調節が必要とな
る。 従って本発明は、前述の例ではなく、以下の請求の範
囲及び均等物のみによって制限される。
─────────────────────────────────────────────────────
フロントページの続き
(72)発明者 ダッチノースキ,ポール
アメリカ合衆国ニューヨーク州10309,
スタテン・アイランド,アイオナ・アベ
ニュー 934
(56)参考文献 特開 昭57−73656(JP,A)
特開 昭57−66742(JP,A)
米国特許4361154(US,A)
Claims (1)
- (57)【特許請求の範囲】 1.人体以外の生体の軟組織によって囲まれている骨ば
った部位の状態を評価する方法であって、 a.骨ばった部位に近い第1の位置で骨ばった部位の中
に、2以上の周波数成分を有する超音波信号を送り出
し、 b.骨ばった部位に近い第2の位置で超音波スペクトルを
モニタし、 c.前記モニタされた超音波スペクトル中のエネルギに対
応する受信信号を電気信号として発生し、 d.前記受信信号から1以上の特有の信号を抽出し、 e.前記特有の信号を1以上の骨ばった部位の状態に関す
る別の特有の信号と直接比較し、 f.前記特有の信号と前記1以上の別の特有の信号との間
の直接の比較に基づいて骨ばった部位の条件を決定す
る、 各段階を含むことを特徴とする骨ばった部位の状態を評
価する方法。 2.前記第1の位置と前記第2の位置は、実質的に前記
骨ばった部位の反対側にあることを特徴とする請求項1
に記載の方法。 3.前記抽出する段階は、 a.前記第1の位置から前記第2の位置への距離を測定
し、 b.前記超音波信号に対応するエネルギがモニタされる時
間と前記超音波信号が送り出される時間の間の差として
1以上の遷移時間を決定し、 c.前記遷移時間と前記距離から1以上の速度を計算して
前記特有の信号を生成する、 各段階を含むことを特徴とする請求項1に記載の方法。 4.前記抽出する段階はさらに、 a.前記受信信号と1以上のテンプレート信号とを比較
し、 b.前記受信信号の前記の比較に従って前記受信信号を修
正して前記受信信号におけるノイズの効果を最小化す
る、 各段階を含むことを特徴とする請求項1に記載の方法。 5.前記別の特有の信号は前記受信信号から抽出される
ことを特徴とする請求項1に記載の方法。 6.前記別の特有の信号は異なる集団の特有の信号を含
むことを特徴とする請求項1に記載の方法。 7.前記超音波信号は100KHzから3MHzの範囲にある周波
数を含むことを特徴とする方請求項1に記載の方法。 8.前記のモニタする段階と前記の発生する段階はデー
タ取得システムにより実行されることを特徴とする請求
項1に記載の方法。 9.前記遷移時間と前記距離から1以上の速度を計算す
る段階は、骨ばった部位によって囲まれた軟組織の特徴
によって前記速度を修正する段階を含むことを特徴とす
る請求項3に記載の方法。 10.前記特有の信号は、第1の所定の周波数よりも大
きい周波数を有する前記受信信号の周波数成分から生成
された高周波数時間領域信号、及び第2所定周波数より
も少ない周波数を有する前記受信信号の周波数成分から
生成された低周波数時間領域信号を含むことを特徴とす
る請求項1に記載の方法。 11.生体の軟組織によって囲まれる骨ばった部位の状
態を評価する装置であって、 骨ばった部位に近い第1の位置で生体の中に、2以上の
周波数成分を有する超音波信号を送り出す手段と、 骨ばった部位に近い第2の位置で超音波スペクトルをモ
ニタする手段と、 前記モニタされた超音波スペクトル中のエネルギに対応
する受信信号を電気信号として発生する手段と、 前記受信信号から1以上の特有の信号を抽出する手段
と、 前記特有の信号を1以上の骨ばった部位の状態に関する
別の特有の信号と直接比較する手段と、 前記特有の信号と前記1以上の別の特有の信号との間の
直接の比較に基づいて骨ばった部位の条件を自動的に決
定する手段と、 を含む骨ばった部位の状態を評価する装置。 12.前記第1の位置と前記第2の位置は、実質的に前
記骨ばった部位の反対側にあることを特徴とする請求項
11に記載の装置。 13.前記抽出する手段は、 前記第1の位置から前記第2の位置への距離を測定する
手段と、 前記超音波信号に対応するエネルギがモニタされる時間
と前記超音波信号が送り出される時間の間の差として1
以上の遷移時間を決定する手段と 前記遷移時間と前記距離から1以上の速度を計算して前
記特有の信号を生成する手段と、 を含むことを特徴とする請求項11に記載の装置。 14.前記抽出する手段はさらに、 前記受信信号と1以上のテンプレート信号とを比較する
手段と、 前記受信信号と前記の比較に従って前記受信信号を修正
して前記受信信号におけるノイズの効果を最小化する手
段と、 を含むことを特徴とする請求項11に記載の装置。 15.前記別の特有の信号は前記受信信号から抽出され
ることを特徴とする請求項11に記載の装置。 16.前記別の特有の信号は異なる集団の特有の信号を
含むことを特徴とする請求の範囲第11項に記載の装置。 17.前記超音波信号は100KHzから3MHzの範囲にある周
波数を含むことを特徴とする請求項11に記載の装置。 18.前記のモニタする手段と前記の発生する手段はデ
ータ取得システムを含むことを特徴とする請求項11に記
載の装置。 19.前記遷移時間と前記距離から1以上の速度を計算
する手段は、生体の軟組織の特徴によって前記速度を修
正する手段を含むことを特徴とする請求項13に記載の装
置。 20.前記特有の信号は、第1の所定の周波数よりも大
きい周波数を有する前記受信信号の周波数成分から生成
された高周波数時間領域信号、及び第2所定周波数より
も少ない周波数を有する前記受信信号の周波数成分から
生成された低周波数時間領域信号を含むことを特徴とす
る請求項11に記載の装置。 21.インピーダンス整合材料を含むことを特徴とす
る、請求項11に記載の装置。 22.前記比較に使用するために前記特有の信号のいず
れかを選択する手段を含むことを特徴とする請求項11に
記載の装置。 23.前記別の特有の信号は、以前または将来に測定さ
れる生体用の骨ばった部位の特有の信号を含むことを特
徴とする請求項11に記載の装置。 24.現在、以前に測定されたまたは将来に測定される
生体、または異なる集団グループの骨ばった部位の特有
の信号を表示する手段を含むことを特徴とする請求項11
に記載の装置。 25.前記骨ばった部位の厚さを測定する手段を含むこ
とを特徴とする請求項11に記載の装置。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US06/870,175 US4913157A (en) | 1986-06-03 | 1986-06-03 | Ultrasound method and apparatus for evaluating, in vivo, bone conditions |
US870,175 | 1986-06-03 |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP9357612A Division JPH10201763A (ja) | 1986-06-03 | 1997-12-25 | 生体における骨の状態を評価する装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH01503199A JPH01503199A (ja) | 1989-11-02 |
JP2965153B2 true JP2965153B2 (ja) | 1999-10-18 |
Family
ID=25354905
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP62503764A Expired - Lifetime JP2965153B2 (ja) | 1986-06-03 | 1987-05-29 | 骨の状態を生体において評価する超音波評価方法及び装置 |
JP9357612A Pending JPH10201763A (ja) | 1986-06-03 | 1997-12-25 | 生体における骨の状態を評価する装置 |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP9357612A Pending JPH10201763A (ja) | 1986-06-03 | 1997-12-25 | 生体における骨の状態を評価する装置 |
Country Status (5)
Country | Link |
---|---|
US (1) | US4913157A (ja) |
EP (2) | EP0772999A2 (ja) |
JP (2) | JP2965153B2 (ja) |
AU (1) | AU7543087A (ja) |
WO (1) | WO1987007494A1 (ja) |
Families Citing this family (122)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0299906A3 (en) * | 1987-07-16 | 1990-06-13 | The University Of Melbourne | In-vivo bone quality measurement |
DE8816422U1 (de) * | 1988-05-06 | 1989-08-10 | Siemens AG, 1000 Berlin und 8000 München | Hörhilfegerät mit drahtloser Fernsteuerung |
US6277076B1 (en) | 1988-05-11 | 2001-08-21 | Lunar Corporation | Ultrasonic densitometer with pre-inflated fluid coupling membranes |
US6027449A (en) * | 1988-05-11 | 2000-02-22 | Lunar Corporation | Ultrasonometer employing distensible membranes |
US5840029A (en) * | 1988-05-11 | 1998-11-24 | Lunar Corporation | Imaging ultrasonic densitometer |
US5483965A (en) * | 1988-05-11 | 1996-01-16 | Lunar Corporation | Ultrasonic densitometer device and method |
US4930511A (en) * | 1988-05-11 | 1990-06-05 | Lunar Radiation, Inc. | Ultrasonic densitometer device and method |
US5603325A (en) * | 1988-05-11 | 1997-02-18 | Lunar Corporation | Ultrasonic densitometer with width compensation |
US5038787A (en) * | 1988-08-10 | 1991-08-13 | The Board Of Regents, The University Of Texas System | Method and apparatus for analyzing material properties using reflected ultrasound |
US4926870A (en) * | 1988-08-30 | 1990-05-22 | Osteo-Technology, Inc. | Method and apparatus for ultrasonic analysis of bone strength in vivo |
US5247937A (en) * | 1989-11-17 | 1993-09-28 | Board Of Regents, The University Of Texas System | Transaxial compression technique for sound velocity estimation |
US5107837A (en) * | 1989-11-17 | 1992-04-28 | Board Of Regents, University Of Texas | Method and apparatus for measurement and imaging of tissue compressibility or compliance |
US5095465A (en) * | 1990-01-05 | 1992-03-10 | Board Of Regents The University Of Texas System | In situ testing with surface seismic waves of materials having properties that change with time |
US5100318A (en) * | 1990-04-13 | 1992-03-31 | Periosonics, Inc. | Ultrasonic method and apparatus for measuring the periodontal pocket |
US5426979A (en) * | 1990-06-04 | 1995-06-27 | Medicano Systems Ltd. | Frequency spectrum apparatus for determining mechanical properties |
IL94616A (en) * | 1990-06-04 | 1994-06-24 | Medicano Systems Ltd | Apparatus and method for calculating the mechanical properties of a solid |
GB2257253B (en) * | 1991-06-17 | 1995-01-11 | Christian Mcdonald Langton | Ultrasound bone analyser |
US5253530A (en) * | 1991-08-12 | 1993-10-19 | Letcher Iii John H | Method and apparatus for reflective ultrasonic imaging |
US5218963A (en) * | 1991-10-15 | 1993-06-15 | Lunar Corporation | Ultrasonic bone analysis device and method |
DE69329168T2 (de) * | 1992-05-20 | 2001-04-05 | Aloka Co. Ltd., Mitaka | Vorrichtung zur Bestimmung der Eigenschaften von Knochen |
GB9213220D0 (en) * | 1992-06-22 | 1992-08-05 | Langton Christian M | Ultrasound bone analyser |
US5458130A (en) * | 1992-07-30 | 1995-10-17 | Orthologic Corporation | Ultrasonic therapy and assessment apparatus and method |
US5259384A (en) * | 1992-07-30 | 1993-11-09 | Kaufman Jonathan J | Ultrasonic bone-assessment apparatus and method |
US5396891A (en) * | 1993-04-07 | 1995-03-14 | Osteo Sciences Corporation | System and method for external acoustic bone velocity measurement |
US5720290A (en) | 1993-04-07 | 1998-02-24 | Metra Biosystems, Inc. | Apparatus and method for acoustic analysis of bone using optimized functions of spectral and temporal signal components |
DE4315794C2 (de) * | 1993-05-13 | 1995-09-21 | Nukem Gmbh | Verfahren und Vorrichtung zur zerstörungsfreien Prüfung von Gegenständen mit Ultraschall |
US5723791A (en) * | 1993-09-28 | 1998-03-03 | Defelsko Corporation | High resolution ultrasonic coating thickness gauge |
US5931780A (en) * | 1993-11-29 | 1999-08-03 | Arch Development Corporation | Method and system for the computerized radiographic analysis of bone |
IT1268599B1 (it) * | 1994-01-14 | 1997-03-06 | Igea Srl | Sistema di misura ad ultrasuoni per la rilevazione della densita' e struttura ossea. |
US5535750A (en) * | 1994-09-30 | 1996-07-16 | Kabushiki Kaisha Ishikawa Seisakusho, Ltd. | Method and apparatus for evaluating the progress of osteoporosis by ultrasonic signals |
US5749363A (en) * | 1994-12-14 | 1998-05-12 | Sekisui Kagaku Kogyo Kabushiki Kaisya | Osteoporosis diagnosing apparatus and method |
JP2840040B2 (ja) * | 1994-12-22 | 1998-12-24 | アロカ株式会社 | 組織内音速測定方法 |
US5882303A (en) * | 1995-01-25 | 1999-03-16 | Stuessi; Edgar | Process and arrangement for determining the dispersion properties of mechanical waves in a three-dimensional object |
US6517487B1 (en) | 1995-03-01 | 2003-02-11 | Lunar Corporation | Ultrasonic densitometer with opposed single transducer and transducer array |
JP2883290B2 (ja) * | 1995-04-10 | 1999-04-19 | アロカ株式会社 | 超音波骨評価装置 |
US5782763A (en) * | 1995-06-01 | 1998-07-21 | Cyberlogic Inc. | Electromagnetic bone-assessment apparatus and method |
US6213934B1 (en) | 1995-06-01 | 2001-04-10 | Hyper3D Corp. | Electromagnetic bone-assessment and treatment: apparatus and method |
US6004272A (en) * | 1995-06-07 | 1999-12-21 | Hologic, Inc. | Ultrasonic bone testing apparatus with repeatable positioning and repeatable coupling |
US5755228A (en) * | 1995-06-07 | 1998-05-26 | Hologic, Inc. | Equipment and method for calibration and quality assurance of an ultrasonic bone anaylsis apparatus |
US6352512B1 (en) | 1995-06-07 | 2002-03-05 | Hologic, Inc. | Bone analysis apparatus and method for calibration and quality assurance of an ultrasonic bone analysis apparatus |
US5730135A (en) * | 1995-09-12 | 1998-03-24 | Horiba, Ltd. | Ultrasonic bone diagnostic apparatus and method |
IL116701A0 (en) * | 1995-10-04 | 1996-10-16 | Sunlight Ultrasound Technologi | Ultrasonic device for determining bone characteristics |
JP2000500998A (ja) * | 1995-11-22 | 2000-02-02 | オステオメータ メディテック エイ/エス | 骨の状態を評価する方法及び装置 |
US5651363A (en) * | 1996-02-16 | 1997-07-29 | Orthologic Corporation | Ultrasonic bone assessment method and apparatus |
US5879301A (en) * | 1996-02-16 | 1999-03-09 | Orthologic Corp. | Ultrasonic bone assessment method and apparatus |
US5785041A (en) * | 1996-03-26 | 1998-07-28 | Hologic Inc. | System for assessing bone characteristics |
US6050950A (en) | 1996-12-18 | 2000-04-18 | Aurora Holdings, Llc | Passive/non-invasive systemic and pulmonary blood pressure measurement |
US7108663B2 (en) * | 1997-02-06 | 2006-09-19 | Exogen, Inc. | Method and apparatus for cartilage growth stimulation |
US7789841B2 (en) * | 1997-02-06 | 2010-09-07 | Exogen, Inc. | Method and apparatus for connective tissue treatment |
US5904659A (en) | 1997-02-14 | 1999-05-18 | Exogen, Inc. | Ultrasonic treatment for wounds |
US6226393B1 (en) * | 1997-07-04 | 2001-05-01 | Torsana Osteoporosis Diagnostics A/S | Method for estimating the bone quality or skeletal status of a vertebrate |
US6132374A (en) * | 1997-08-01 | 2000-10-17 | Acuson Corporation | Ultrasonic imaging method and system |
US7112173B1 (en) * | 1998-03-03 | 2006-09-26 | Sunlight Medical Ltd. | Determination of acoustic velocity in bone |
US6585649B1 (en) * | 1998-05-02 | 2003-07-01 | John D. Mendlein | Methods and devices for improving ultrasonic measurements using multiple angle interrogation |
US6077224A (en) * | 1998-03-23 | 2000-06-20 | Lang; Philipp | Methods and device for improving broadband ultrasonic attenuation and speed of sound measurements using anatomical landmarks |
US6013031A (en) * | 1998-03-09 | 2000-01-11 | Mendlein; John D. | Methods and devices for improving ultrasonic measurements using anatomic landmarks and soft tissue correction |
US6024711A (en) * | 1998-03-12 | 2000-02-15 | The University Of British Columbia | Diagnosis of osteoporosis using acoustic emissions |
JP2002513661A (ja) * | 1998-05-06 | 2002-05-14 | エクソジェン インコーポレイテッド | 超音波包帯 |
US6116244A (en) * | 1998-06-02 | 2000-09-12 | Acuson Corporation | Ultrasonic system and method for three-dimensional imaging with opacity control |
US5957852A (en) * | 1998-06-02 | 1999-09-28 | Acuson Corporation | Ultrasonic harmonic imaging system and method |
US6511426B1 (en) | 1998-06-02 | 2003-01-28 | Acuson Corporation | Medical diagnostic ultrasound system and method for versatile processing |
AU5117699A (en) | 1998-07-21 | 2000-02-14 | Acoustic Sciences Associates | Synthetic structural imaging and volume estimation of biological tissue organs |
US6048316A (en) * | 1998-10-16 | 2000-04-11 | Acuson Corporation | Medical diagnostic ultrasonic imaging system and method for displaying composite fundamental and harmonic images |
CA2348467A1 (en) | 1998-11-08 | 2000-05-18 | Sunlight Medical, Ltd. | Bone velocity determination |
US6234969B1 (en) * | 1999-02-01 | 2001-05-22 | Hologic, Inc. | Bone sonometry, densitometry and imaging |
CA2370960C (en) * | 1999-04-20 | 2006-06-13 | Synthes (U.S.A.) | Device for the percutaneous obtainment of 3d-coordinates on the surface of a human or animal organ |
US6149594A (en) * | 1999-05-05 | 2000-11-21 | Agilent Technologies, Inc. | Automatic ultrasound measurement system and method |
IL130618A (en) * | 1999-06-23 | 2004-01-04 | Samer M Srouji | Ultrasound system for use in drilling implant cavities |
US7374538B2 (en) * | 2000-04-05 | 2008-05-20 | Duke University | Methods, systems, and computer program products for ultrasound measurements using receive mode parallel processing |
US6371912B1 (en) * | 2000-04-05 | 2002-04-16 | Duke University | Method and apparatus for the identification and characterization of regions of altered stiffness |
US6899680B2 (en) * | 2000-10-19 | 2005-05-31 | Odetect As | Ultrasound measurement techniques for bone analysis |
CA2426903C (en) * | 2000-10-25 | 2011-12-06 | Exogen, Inc. | Transducer mounting assembly |
WO2002071949A2 (en) | 2001-02-28 | 2002-09-19 | Research Foundation Of State University Of New York | Method and apparatus for scanning confocal acoustic diagnostic for bone quality |
WO2002085178A2 (en) * | 2001-04-20 | 2002-10-31 | Imadent Ltd. | Method and apparatus for non-invasive ultrasonic imaging of hard tissue |
US20050130321A1 (en) * | 2001-04-23 | 2005-06-16 | Nicholson Jeremy K. | Methods for analysis of spectral data and their applications |
AU2002321792A1 (en) * | 2001-07-24 | 2003-02-17 | Sunlight Medical, Ltd. | Method and apparatus for bone diagnosis |
JP2004535882A (ja) * | 2001-07-24 | 2004-12-02 | サンライト メディカル リミテッド | 超音波を用いた骨年齢評価方法 |
DE10142364A1 (de) * | 2001-08-30 | 2003-04-03 | Advanced Acoustix Gmbh | Verfahren zur Abstandsbestimmung mittels Ultraschall an schallweichen Objekten |
US6784363B2 (en) | 2001-10-02 | 2004-08-31 | Parker-Hannifin Corporation | EMI shielding gasket construction |
ITTO20011135A1 (it) * | 2001-12-05 | 2003-06-05 | Igea Srl | Sistema elettronico per la valutazione densitometrica-strutturale deltessuto osseo e per la stimolazione dell'osteogenesi in campo odontoia |
US6827686B2 (en) * | 2002-08-21 | 2004-12-07 | Koninklijke Philips Electronics N.V. | System and method for improved harmonic imaging |
US20100135458A1 (en) * | 2003-07-18 | 2010-06-03 | Neeraj Agrawal | X-Ray Apparatus for Bone Density Assessment and Monitoring |
US8085898B2 (en) * | 2009-05-08 | 2011-12-27 | Osteometer Meditech, Inc. | Apparatus for bone density assessment and monitoring |
US7275547B2 (en) * | 2003-10-08 | 2007-10-02 | Boston Scientific Scimed, Inc. | Method and system for determining the location of a medical probe using a reference transducer array |
US7285093B2 (en) | 2003-10-10 | 2007-10-23 | Imadent Ltd. | systems for ultrasonic imaging of a jaw, methods of use thereof and coupling cushions suited for use in the mouth |
CA2551325A1 (en) * | 2003-12-30 | 2005-07-21 | Liposonix, Inc. | Component ultrasound transducer |
US20050215899A1 (en) * | 2004-01-15 | 2005-09-29 | Trahey Gregg E | Methods, systems, and computer program products for acoustic radiation force impulse (ARFI) imaging of ablated tissue |
US8202219B2 (en) * | 2004-02-23 | 2012-06-19 | Cyberlogic, Inc. | Ultrasonic bone assessment apparatus and method |
US7315639B2 (en) * | 2004-03-03 | 2008-01-01 | Mevis Gmbh | Method of lung lobe segmentation and computer system |
GB0414318D0 (en) * | 2004-06-25 | 2004-07-28 | Tissuomics Ltd | Analysing body tissue |
AU2005205820B2 (en) * | 2004-09-04 | 2011-04-14 | Smith & Nephew Plc | Ultrasound device and method of use |
US7223249B2 (en) * | 2005-03-17 | 2007-05-29 | Rosenberg Jeffrey M | Method and apparatus for determining the depth of a gingival pocket |
US7738683B2 (en) * | 2005-07-22 | 2010-06-15 | Carestream Health, Inc. | Abnormality detection in medical images |
WO2007012203A2 (en) * | 2005-07-28 | 2007-02-01 | University Of Windsor | Ultrasonic sensor for dental applications |
WO2007120890A2 (en) * | 2006-04-13 | 2007-10-25 | The Research Foundation Of State University Of New York | Phased array ultrasound with electronically controlled focal point for assessing bone quality via acoustic topology and wave transmit functions |
US7691149B2 (en) * | 2006-05-15 | 2010-04-06 | Biomet Manufacturing Corp. | Porous titanium modular revision patella system |
US8419643B2 (en) * | 2006-09-21 | 2013-04-16 | Artann Laboratories Inc. | Ultrasonic method and apparatus for assessment of bone |
WO2008038159A2 (en) * | 2006-09-29 | 2008-04-03 | Odetect As | Ultrasound measurement techniques for bone analysis |
US20100063394A1 (en) * | 2006-11-14 | 2010-03-11 | Kagoshima University | Osteoporosis progress stage measuring device and osteoporosis progress stage measuring method |
US7901356B2 (en) * | 2006-12-13 | 2011-03-08 | Cyberlogic, Inc. | Ultrasonic bone assessment apparatus and method |
US9615814B2 (en) * | 2006-12-13 | 2017-04-11 | Cyberlogic, Inc. | Ultrasonic bone assessment apparatus and method |
US7862510B2 (en) * | 2007-02-09 | 2011-01-04 | Cyberlogic, Inc. | Ultrasonic bone assessment apparatus and method |
US8679019B2 (en) * | 2007-12-03 | 2014-03-25 | Bone Index Finland Oy | Method for measuring of thicknesses of materials using an ultrasound technique |
US8187187B2 (en) | 2008-07-16 | 2012-05-29 | Siemens Medical Solutions Usa, Inc. | Shear wave imaging |
US8696754B2 (en) * | 2008-09-03 | 2014-04-15 | Biomet Manufacturing, Llc | Revision patella prosthesis |
EP2169424A1 (en) | 2008-09-03 | 2010-03-31 | Esaote S.p.A. | Method and device for determining local characteristics of an object, particularly with ultrasounds |
US20100210942A1 (en) * | 2008-11-26 | 2010-08-19 | New Jersey Institute Of Technology | System and method for ultrasound analysis of biological structures |
US20100210943A1 (en) * | 2009-02-18 | 2010-08-19 | West Virginia University Research Corporation | Systems and Methods for Echoperiodontal Imaging |
US20100241001A1 (en) * | 2009-03-20 | 2010-09-23 | Palmeri Mark L | Ultrasound Methods, Systems and Computer Program Products for Imaging Fluids |
US9538938B2 (en) | 2009-10-02 | 2017-01-10 | Koninklijke Philips N.V. | Signal transmission through a medium |
EP2366997B1 (en) | 2010-03-17 | 2012-07-18 | Esaote S.p.A. | Method and device for determining the structural organization of an object with ultrasounds |
FI126104B (fi) * | 2010-09-09 | 2016-06-30 | Bone Index Finland Oy | Menetelmä ja järjestelmä luun mineraalitiheyden arvioimiseksi |
US9486291B2 (en) | 2012-06-21 | 2016-11-08 | Rivanna Medical Llc | Target region identification for imaging applications |
DE102012220468B3 (de) * | 2012-11-09 | 2014-03-27 | Seca Ag | Längenmessgerät |
KR101889787B1 (ko) * | 2012-12-04 | 2018-08-20 | 삼성전자주식회사 | 신호 처리 방법 및 장치 |
US11147536B2 (en) | 2013-02-28 | 2021-10-19 | Rivanna Medical Llc | Localization of imaging target regions and associated systems, devices and methods |
US10134125B2 (en) | 2013-02-28 | 2018-11-20 | Rivanna Medical Llc | Systems and methods for ultrasound imaging |
CN104546014A (zh) * | 2014-12-25 | 2015-04-29 | 中国科学院深圳先进技术研究院 | 一种用于生物组织弹性测量的剪切波速度估计方法 |
US10548564B2 (en) | 2015-02-26 | 2020-02-04 | Rivanna Medical, LLC | System and method for ultrasound imaging of regions containing bone structure |
JP6601320B2 (ja) * | 2016-06-16 | 2019-11-06 | コニカミノルタ株式会社 | 超音波診断装置、及び超音波診断装置の制御方法 |
JP6538289B1 (ja) | 2016-07-06 | 2019-07-03 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 骨特性を監視する方法 |
CN106691513B (zh) * | 2016-11-22 | 2019-09-24 | 北京百思声创科技有限公司 | 超声波骨密度测量仪 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4361154A (en) | 1978-07-28 | 1982-11-30 | Massachusetts Institute Of Technology | Method for establishing, in vivo, bone strength |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3512124A (en) * | 1960-02-29 | 1970-05-12 | Mc Donnell Douglas Corp | Signal analyzing and indicating method and apparatus |
US3477422A (en) * | 1965-10-15 | 1969-11-11 | John M Jurist Jr | Vibratory bone density determination method and apparatus |
US3847141A (en) * | 1973-08-08 | 1974-11-12 | Nasa | Ultrasonic bone densitometer |
US3952280A (en) * | 1974-01-10 | 1976-04-20 | Esl Incorporated | Radiation monitoring of an object space with a clutter suppression technique |
US4421119A (en) * | 1979-06-15 | 1983-12-20 | Massachusetts Institute Of Technology | Apparatus for establishing in vivo, bone strength |
US4233845A (en) * | 1978-07-28 | 1980-11-18 | Massachusetts Institute Of Technology | Method of assessing performance potential of a quadruped |
JPS5766742A (en) * | 1980-10-09 | 1982-04-23 | Fujitsu Ltd | Ultrasonic diagnostic device |
JPS5773656A (en) * | 1980-10-26 | 1982-05-08 | Yukio Yamamoto | Double-frequency scanning type microwave computer |
GB2156983B (en) * | 1983-08-26 | 1987-03-11 | Mintowt Czyz Witek | Method and apparatus for assessing the structure and mechanical integrity of osseous systems |
JPS6249834A (ja) * | 1985-08-30 | 1987-03-04 | 富士物産株式会社 | 生体の骨性検査法及びその装置 |
US4774959A (en) * | 1986-01-10 | 1988-10-04 | Walker Sonix A/S | Narrow band ultrasonic frequency attentuation bone measurement system |
-
1986
- 1986-06-03 US US06/870,175 patent/US4913157A/en not_active Expired - Lifetime
-
1987
- 1987-05-29 EP EP96120199A patent/EP0772999A2/en not_active Withdrawn
- 1987-05-29 AU AU75430/87A patent/AU7543087A/en not_active Abandoned
- 1987-05-29 EP EP87904166A patent/EP0308426A1/en not_active Withdrawn
- 1987-05-29 JP JP62503764A patent/JP2965153B2/ja not_active Expired - Lifetime
- 1987-05-29 WO PCT/US1987/001296 patent/WO1987007494A1/en not_active Application Discontinuation
-
1997
- 1997-12-25 JP JP9357612A patent/JPH10201763A/ja active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4361154A (en) | 1978-07-28 | 1982-11-30 | Massachusetts Institute Of Technology | Method for establishing, in vivo, bone strength |
Also Published As
Publication number | Publication date |
---|---|
EP0308426A1 (en) | 1989-03-29 |
AU7543087A (en) | 1988-01-11 |
JPH10201763A (ja) | 1998-08-04 |
EP0772999A3 (ja) | 1997-06-25 |
JPH01503199A (ja) | 1989-11-02 |
EP0772999A2 (en) | 1997-05-14 |
WO1987007494A1 (en) | 1987-12-17 |
US4913157A (en) | 1990-04-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2965153B2 (ja) | 骨の状態を生体において評価する超音波評価方法及び装置 | |
US4941474A (en) | Multivariable analysis of bone condition | |
US4926870A (en) | Method and apparatus for ultrasonic analysis of bone strength in vivo | |
CN100548224C (zh) | 利用超声检测内部狭窄以识别由其引起的组织振动的设备 | |
Laugier | Instrumentation for in vivo ultrasonic characterization of bone strength | |
US5143069A (en) | Diagnostic method of monitoring skeletal defect by in vivo acoustic measurement of mechanical strength using correlation and spectral analysis | |
US6468215B1 (en) | Method and device for multi-parametric ultrasonic assessment of bone conditions | |
US8419643B2 (en) | Ultrasonic method and apparatus for assessment of bone | |
KR101900122B1 (ko) | 환자의 뼈 조직의 품질을 평가하기 위한 초음파 장치 | |
US6086538A (en) | Methods and apparatus for evaluation of bone condition | |
KR100548182B1 (ko) | 광대역밴드 초음파 반사량을 이용한 종골의 골밀도 측정장치 및 측정방법 | |
Maricic et al. | Bone densitometry | |
KR101197923B1 (ko) | 펄스-에코 초음파를 이용한 피질골 두께 측정방법 및 그 장치 | |
KR100581229B1 (ko) | 램파를 이용한 경골의 골밀도 측정방법 | |
CN111110190A (zh) | 一种利用光声时域信号评估骨质弹性模量的方法 | |
RU2790947C1 (ru) | Способ ультразвуковой велосимметрии для оценки состояния твердых тканей зубов | |
RU2289317C2 (ru) | Способ регистрации биомеханических свойств длинных трубчатых костей | |
JP2000175911A (ja) | 骨強度測定方法及び装置 | |
Maia et al. | Broadband ultrasound attenuation in the calcaneal region: A comparative study of single-position versus scanning systems | |
Tatarinov et al. | The dual-frequency method for ultrasonic assessment of skeletal system | |
Zhang et al. | Lung ultrasound surface wave elastography | |
RU2071274C1 (ru) | Способ диагностики состояния костных тканей с помощью ультразвука | |
Laugier et al. | Ultrasound parametric imaging of bone in vivo | |
JPH09234200A (ja) | 超音波による骨の診断方法 | |
US20130204132A1 (en) | System and method for ultrasound analysis of biological structures |