JP2008514368A - Characteristic determination method, recording apparatus, evaluation method - Google Patents
Characteristic determination method, recording apparatus, evaluation method Download PDFInfo
- Publication number
- JP2008514368A JP2008514368A JP2007534740A JP2007534740A JP2008514368A JP 2008514368 A JP2008514368 A JP 2008514368A JP 2007534740 A JP2007534740 A JP 2007534740A JP 2007534740 A JP2007534740 A JP 2007534740A JP 2008514368 A JP2008514368 A JP 2008514368A
- Authority
- JP
- Japan
- Prior art keywords
- model
- biomechanical
- structural
- deformation
- bone
- 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 85
- 238000011156 evaluation Methods 0.000 title claims description 14
- 241001465754 Metazoa Species 0.000 claims abstract description 14
- 210000000988 bone and bone Anatomy 0.000 claims description 159
- 238000011282 treatment Methods 0.000 claims description 51
- 210000001519 tissue Anatomy 0.000 claims description 37
- 238000002591 computed tomography Methods 0.000 claims description 36
- 230000000694 effects Effects 0.000 claims description 32
- 210000000056 organ Anatomy 0.000 claims description 20
- 210000004204 blood vessel Anatomy 0.000 claims description 18
- 201000010099 disease Diseases 0.000 claims description 18
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 claims description 18
- 230000032683 aging Effects 0.000 claims description 17
- 230000035882 stress Effects 0.000 claims description 14
- 238000012512 characterization method Methods 0.000 claims description 10
- 230000002526 effect on cardiovascular system Effects 0.000 claims description 7
- 238000005452 bending Methods 0.000 claims description 5
- 210000002216 heart Anatomy 0.000 claims description 5
- 230000002093 peripheral effect Effects 0.000 claims description 5
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000001727 in vivo Methods 0.000 claims description 3
- 238000002604 ultrasonography Methods 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 claims 2
- 238000012831 peritoneal equilibrium test Methods 0.000 claims 1
- 238000012636 positron electron tomography Methods 0.000 claims 1
- 238000012877 positron emission topography Methods 0.000 claims 1
- 238000002601 radiography Methods 0.000 claims 1
- 230000008859 change Effects 0.000 abstract description 26
- 238000004458 analytical method Methods 0.000 description 30
- 239000000463 material Substances 0.000 description 29
- 238000005259 measurement Methods 0.000 description 23
- 230000006399 behavior Effects 0.000 description 19
- 208000001132 Osteoporosis Diseases 0.000 description 17
- 230000037182 bone density Effects 0.000 description 16
- 238000009826 distribution Methods 0.000 description 14
- 238000002651 drug therapy Methods 0.000 description 14
- 230000004044 response Effects 0.000 description 14
- 230000037118 bone strength Effects 0.000 description 12
- 230000006870 function Effects 0.000 description 12
- 238000002595 magnetic resonance imaging Methods 0.000 description 12
- 230000017531 blood circulation Effects 0.000 description 10
- 238000009547 dual-energy X-ray absorptiometry Methods 0.000 description 9
- 239000002639 bone cement Substances 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 230000002308 calcification Effects 0.000 description 7
- 239000003814 drug Substances 0.000 description 7
- 229940079593 drug Drugs 0.000 description 7
- 238000012986 modification Methods 0.000 description 7
- 230000004048 modification Effects 0.000 description 7
- 208000006386 Bone Resorption Diseases 0.000 description 6
- 230000024279 bone resorption Effects 0.000 description 6
- 238000003745 diagnosis Methods 0.000 description 6
- 239000012530 fluid Substances 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 238000002560 therapeutic procedure Methods 0.000 description 6
- 206010017076 Fracture Diseases 0.000 description 5
- 230000009286 beneficial effect Effects 0.000 description 5
- 238000005094 computer simulation Methods 0.000 description 5
- 230000005477 standard model Effects 0.000 description 5
- 210000000689 upper leg Anatomy 0.000 description 5
- 238000013459 approach Methods 0.000 description 4
- 206010003246 arthritis Diseases 0.000 description 4
- 210000000845 cartilage Anatomy 0.000 description 4
- 239000007943 implant Substances 0.000 description 4
- 230000007246 mechanism Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000002054 transplantation Methods 0.000 description 4
- 210000000707 wrist Anatomy 0.000 description 4
- 208000024172 Cardiovascular disease Diseases 0.000 description 3
- 229910052500 inorganic mineral Inorganic materials 0.000 description 3
- 239000011707 mineral Substances 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 230000001009 osteoporotic effect Effects 0.000 description 3
- 238000011002 quantification Methods 0.000 description 3
- 206010041569 spinal fracture Diseases 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 230000001225 therapeutic effect Effects 0.000 description 3
- 208000010392 Bone Fractures Diseases 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 238000002583 angiography Methods 0.000 description 2
- 210000001185 bone marrow Anatomy 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 238000010968 computed tomography angiography Methods 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 238000013500 data storage Methods 0.000 description 2
- 230000036541 health Effects 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 210000003205 muscle Anatomy 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000011369 optimal treatment Methods 0.000 description 2
- 238000001050 pharmacotherapy Methods 0.000 description 2
- 238000002600 positron emission tomography Methods 0.000 description 2
- 238000001356 surgical procedure Methods 0.000 description 2
- 210000002435 tendon Anatomy 0.000 description 2
- 230000002792 vascular Effects 0.000 description 2
- OGSPWJRAVKPPFI-UHFFFAOYSA-N Alendronic Acid Chemical compound NCCCC(O)(P(O)(O)=O)P(O)(O)=O OGSPWJRAVKPPFI-UHFFFAOYSA-N 0.000 description 1
- 238000012935 Averaging Methods 0.000 description 1
- 206010061818 Disease progression Diseases 0.000 description 1
- 206010020100 Hip fracture Diseases 0.000 description 1
- 241000282412 Homo Species 0.000 description 1
- 206010053649 Vascular rupture Diseases 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 229940062527 alendronate Drugs 0.000 description 1
- 210000001188 articular cartilage Anatomy 0.000 description 1
- 230000003542 behavioural effect Effects 0.000 description 1
- 230000004071 biological effect Effects 0.000 description 1
- 239000012620 biological material Substances 0.000 description 1
- 230000033558 biomineral tissue development Effects 0.000 description 1
- 238000001574 biopsy Methods 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 210000001124 body fluid Anatomy 0.000 description 1
- 239000010839 body fluid Substances 0.000 description 1
- 210000000748 cardiovascular system Anatomy 0.000 description 1
- 230000015271 coagulation Effects 0.000 description 1
- 238000005345 coagulation Methods 0.000 description 1
- 230000001054 cortical effect Effects 0.000 description 1
- 230000034994 death Effects 0.000 description 1
- 231100000517 death Toxicity 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 238000001739 density measurement Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000005750 disease progression Effects 0.000 description 1
- 230000000857 drug effect Effects 0.000 description 1
- 238000005206 flow analysis Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000002650 habitual effect Effects 0.000 description 1
- 208000019622 heart disease Diseases 0.000 description 1
- 238000000265 homogenisation Methods 0.000 description 1
- 230000001771 impaired effect Effects 0.000 description 1
- 238000002513 implantation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000000338 in vitro Methods 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 210000003127 knee Anatomy 0.000 description 1
- 210000000629 knee joint Anatomy 0.000 description 1
- 210000003041 ligament Anatomy 0.000 description 1
- 238000007620 mathematical function Methods 0.000 description 1
- 230000010534 mechanism of action Effects 0.000 description 1
- 238000002483 medication Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 208000010125 myocardial infarction Diseases 0.000 description 1
- 239000002547 new drug Substances 0.000 description 1
- 230000011164 ossification Effects 0.000 description 1
- 201000008482 osteoarthritis Diseases 0.000 description 1
- 238000011321 prophylaxis Methods 0.000 description 1
- 230000005180 public health Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000012916 structural analysis Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
- 210000002303 tibia Anatomy 0.000 description 1
- 230000017423 tissue regeneration Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/508—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for non-human patients
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Prostheses (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
人間又は動物の生体の構造の画像から、この構造の１又は複数の構造／機能特性を判定する特性判定方法では、画像に基づいて、構造の構造モデルを生成する。構造モデルに基づいて、第１の生体力学的量を算出する。構造モデルを変更して、変形モデルを生成する。変形モデルに基づいて、第２の生体力学的量を算出する。第１及び第２の生体力学的量を比較し、構造の構造／機能特性を評価する。 In the characteristic determination method for determining one or more structural / functional characteristics of a structure from an image of the structure of a human or animal living body, a structural model of the structure is generated based on the image. A first biomechanical quantity is calculated based on the structural model. Change the structural model to generate a deformation model. A second biomechanical quantity is calculated based on the deformation model. The first and second biomechanical quantities are compared and the structural / functional properties of the structure are evaluated.
Description
本出願は、２００４年９月３０日に出願された米国の仮特許出願番号第６０／６１４，６０５号、発明の名称「骨の構造的質の評価（ASSESSMENT OF BONE STRUCTURAL QUALITY）」の優先権を主張し、この特許文献の開示内容は、引用により本願に援用される。本発明は、米国国立衛生研究所の研究助成金、契約番号ＡＲ４１４８１号に基づく国庫補助によって発明された。米国政府は、本発明に対する所定の権利を有する。 This application is a priority of US Provisional Patent Application No. 60 / 614,605, filed September 30, 2004, entitled “ASSESSMENT OF BONE STRUCTURAL QUALITY”. The disclosure of this patent document is incorporated herein by reference. The present invention was invented by a national subsidy based on a research grant from the National Institutes of Health, contract number AR41481. The US government has certain rights to the invention.
本発明は、人間又は動物の生体の骨又は他の構造の構造／機能特性を評価する評価方法に関する。 The present invention relates to an evaluation method for evaluating the structural / functional properties of bones or other structures of a human or animal organism.
米国国立骨粗鬆症財団及び米国国立衛生研究所（National Institute of Health：ＮＩＨ）によれば、骨粗鬆症は、国民の健康を損なう大きな脅威であり、推定４４００万人の米国人、５０歳以上に限れば、人口の半数以上の骨量が低下している。米国では、今日、既に１０００万人が骨粗鬆に罹患していると推定され、３４００万人の米国人は、骨量が不足し、骨粗鬆症に罹患する危険が高いと言われている。人口の高齢化が進み、国民が長寿になれば、骨粗鬆症になる人口も増え、骨粗鬆症の治療のためのコストは増加し続ける。また、関節炎や心血管疾患等の症状にも同様のことが言える。骨関節炎は、約２０００万人の米国人が患っている。また、約５７００万人の米国人が、ある種の心血管疾患を患い、１年あたり約１００万人（米国の全死亡者の４２％）が死亡していると推定されている。しかしながら、心臓発作を予測することは困難である。そこで、患者に応じて、リスクの高さを予測する技術の向上が望まれている。 According to the National Osteoporosis Foundation and the National Institute of Health (NIH), osteoporosis is a major threat to public health, with an estimated 44 million Americans, 50 years of age or older, More than half of the population has lost bone mass. In the United States, it is estimated that 10 million people are already suffering from osteoporosis, and 34 million Americans are said to be at high risk of suffering from osteoporosis due to lack of bone mass. As the population ages and people become longer-lived, the number of people with osteoporosis increases and the costs for osteoporosis treatment continue to increase. The same applies to symptoms such as arthritis and cardiovascular diseases. Osteoarthritis affects about 20 million Americans. It is also estimated that approximately 57 million Americans suffer from certain cardiovascular diseases and approximately 1 million (42% of all US deaths) die per year. However, it is difficult to predict a heart attack. Therefore, it is desired to improve the technique for predicting the high risk according to the patient.
骨粗鬆症：骨粗鬆症の診断及び測定のための現在の臨床基準である二重エネルギＸ線吸収（dual-energy x-ray absorptiometry：ＤＸＡ（「デキサ」と発音される））スキャンは、骨の生体力学的完全性を評価するための技術として、様々な制約を有する。第１に、ＤＸＡスキャンは、立体的な体積を有する骨密度情報の二次元投影である。したがって、ミネラル密度の面的測定は、ミネラルがどのように構成され、骨内で立体的にどのように分布しているかといった潜在的に重要な構造が無視される。ＤＸＡは、骨梁から皮質を区別せず、背部の要素の存在のために、背骨の測定には不向きである。 Osteoporosis: Dual-energy x-ray absorptiometry (DXA (pronounced “dexa”)) scan, the current clinical standard for diagnosis and measurement of osteoporosis, is a biomechanical study of bone There are various constraints as techniques for assessing completeness. First, a DXA scan is a two-dimensional projection of bone density information having a three-dimensional volume. Thus, the surface density measurement of mineral density ignores potentially important structures such as how minerals are organized and how they are distributed three-dimensionally within the bone. DXA does not distinguish cortex from trabecular bone and is unsuitable for spine measurement due to the presence of dorsal elements.
定量的コンピュータ断層撮影法（Quantitative Computed Tomography：ＱＣＴ、ＣＡＴ又はＣＴスキャンの一種）は、３次元的であり、ＤＸＡスキャンの２次元的性質に起因する制約を克服している。ＣＴスキャンは、３次元手的な性質のために、例えば、骨の前部領域又は後部領域、若しくは、皮質、内部皮質及び骨梁等、領域又は骨の種類によってミネラル密度を区別することができる。しかしながら、ＱＣＴは、単なる放射線評価であり、骨密度及び形状の測定値を提供するだけである。したがって、ＱＣＴスキャンで観察される変化を解釈することは困難である。 Quantitative computed tomography (QCT, CAT or CT scan) is three-dimensional and overcomes the limitations due to the two-dimensional nature of DXA scans. CT scans can distinguish mineral density by type or type of bone, for example, the anterior or posterior region of bone, or the cortex, internal cortex and trabecular bone because of the three-dimensional hand nature . However, QCT is merely a radiation assessment and only provides measurements of bone density and shape. Therefore, it is difficult to interpret the changes observed in the QCT scan.
如何なる構造の機械的挙動も、基本的には、その形状及び材料特性によって支配されている。しかしながら、骨及び人体の殆どの器官に特徴的であるように、形状が不均等であり、材料特性の空間的分布が不均等である、不均一な構造の総合的な機械的応答を算出することは困難である。また、このような形状及び材料特性分布の複雑性のために、このような応答の経時的な又は治療による変化を理解することも難しい。 The mechanical behavior of any structure is basically governed by its shape and material properties. However, it calculates the overall mechanical response of a heterogeneous structure that is non-uniform in shape and spatially distributed in material properties, as is characteristic of bones and most organs of the human body. It is difficult. Also, due to the complexity of such shapes and material property distributions, it is difficult to understand changes in such responses over time or with treatment.
骨全体は、それぞれが固有の物質的な強度特性を有する骨梁及び皮質骨といった２つの異なる種類の骨組織を含む。したがって、骨全体（例えば、椎体、大腿骨近位部、橈骨遠位端）の強度特性は、骨の平均密度、質量及び寸法だけではなく、骨内における骨密度の空間的に変化する分布、骨の立体的形状、及び皮質に対する骨梁の相対的役割等によって決まる。更に、生体内の骨全体には、複雑且つ多軸的に負荷が掛かり、変化する負荷条件の下で破損することがある。例えば、骨粗鬆症の腰骨の骨折は、例えば、歩行等の習慣的な活動に比べて、骨に対して多くの異なる負荷条件を作り出す転倒の際に生じやすい。また、骨の強度は、これらの異なる負荷条件毎に異なる。背骨の骨折については、脊椎骨の強度は、前方に曲がっている場合と、真っ直ぐな場合とでは大きく異なる。 The entire bone includes two different types of bone tissue, trabecular bone and cortical bone, each having unique material strength properties. Thus, the strength characteristics of the entire bone (eg, vertebral body, proximal femur, distal radius) are not only the average density, mass and size of the bone, but also the spatially varying distribution of bone density within the bone. , The three-dimensional shape of the bone, and the relative role of the trabecular bone to the cortex. Furthermore, the entire bone in a living body is complicated and multiaxially loaded, and may be damaged under changing load conditions. For example, osteoporotic hip fractures are more likely to occur during a fall that creates many different loading conditions on the bone than, for example, habitual activities such as walking. Moreover, the strength of the bone is different for each of these different load conditions. For spinal fractures, the strength of the vertebrae differs greatly between bending forward and straight.
心臓血管：心臓血管及びこれに関連する応用例では、体内の多くの血管領域を評価するするために、例えば、デジタル減算血管造影法（digital subtraction angiography：ＤＳＡ）等の技術が用いられる。現在では、ＤＳＡに代えて、より侵襲性が低いＣＴベースの血管造影法も臨床的に用いられている。骨粗鬆症の骨の解析にＱＣＴを適用した場合と同様に、ＣＴ血管造影法は、基底にある臨床的な問題の生体力学的側面を明らかにせず、したがって、ＣＴ画像における情報を最大限に利用できない。血管における血流の制限は、ＣＴ血管造影法による単純な視覚的分析で観察される変化より微妙な血管の変化に影響を受ける場合もある。したがって、ＣＴ画像に基づく血管の生体力学的解析によって、臨床的な問題を更に洞察することができる。 Cardiovascular: Cardiovascular and related applications use techniques such as digital subtraction angiography (DSA), for example, to evaluate many vascular regions in the body. At present, CT-based angiography, which is less invasive, is used clinically instead of DSA. As with QCT applied to the analysis of osteoporotic bone, CT angiography does not reveal the biomechanical aspects of the underlying clinical problem and therefore cannot make the best use of the information in the CT image. . Blood flow limitations in blood vessels may be affected by more subtle changes in blood vessels than those observed with simple visual analysis by CT angiography. Therefore, clinical problems can be further insighted by biomechanical analysis of blood vessels based on CT images.
移植システム：生体に対する様々な種類の移植を行い、負傷又は罹患した生体の一部を修復することができる。例えば、米国では、年間約１５万の腰及び膝のプロテーゼが移植されており、全世界ではこの倍のプロテーゼ移植が行われている。医師は、患者のために最適な移植を選択する必要があるが、所定の医学的徴候について、選択できる機器の選択肢が多ければ、その選択は困難である。構造／機能特性に関する情報を有する患者固有のモデルを用いることができれば、適切な移植を選択する際に、有益な情報が提供される。同様の問題は、例えば、ステント等、心臓血管の治療にも適用され、すなわち、治療を成功させるには、ステントのサイズを適切に選択することが重要となる。このサイズの決定は、患者特有の生体力学的構造／機能特性に基づいて行うことが望ましい。 Transplantation system: Various types of transplantation into the living body can be performed to repair a part of the injured or affected living body. For example, about 150,000 lumbar and knee prostheses are transplanted annually in the United States, and twice as many prosthetic transplants worldwide. The physician needs to select the optimal implant for the patient, but the choice is difficult if there are many device options to choose for a given medical indication. If a patient specific model with information on structural / functional properties can be used, useful information is provided in selecting an appropriate implant. Similar issues apply to cardiovascular treatments such as, for example, stents, ie, proper selection of the size of the stent is important for successful treatment. This size determination is preferably based on patient-specific biomechanical structural / functional properties.
他の応用例：患者固有の構造／機能特性を明らかにする必要がある他の応用例として、関節炎及び脊椎骨折の治療がある。関節軟骨に加わる応力は、関節の下層の骨の密度及び構造に依存すると考えられるため、関節炎については、関節の生体力学的振る舞い、例えば、骨密度、サイズ及び構造に関する患者固有の詳細な情報等を知ることによって、診断及び治療の評価が改善される。背骨の骨折の治療は、骨折した椎体に骨セメントを注入することによって行われる。システムの硬さは、注入された骨セメントの量及び位置に依存するため、このような施術は、施術により得られる骨−骨セメントシステムの構造／機能特性を知ることによって最適化できる。また、金属又は他の種類のプロテーゼを用いる骨折の治療は、寸法及び形状等の因子、プロテーゼの材料、生体部分の密度及び構造等、骨移植システムの構造的な応答を評価することによって改善できる。これにより、医師は、施術の前に、提案された外科的処置の適合性を評価でき、この結果、個々の患者について、最適な処置を選択することができる。現在のところ、殆どの医師は、質的な経験のみに頼り、様々なオプションを評価するための量的な基準を有していない。骨移植システムの構造／機能特性に関する様々な種類の知識によって、治療の結果が改善されることが予想される。 Other applications: Other applications that need to reveal patient-specific structural / functional properties include the treatment of arthritis and vertebral fractures. Since stress on articular cartilage is thought to depend on the density and structure of the underlying bone of the joint, for arthritis, the biomechanical behavior of the joint, such as patient-specific details regarding bone density, size and structure, etc. Knowing the information improves diagnostic and therapeutic evaluation. Treatment of spinal fractures is performed by injecting bone cement into the fractured vertebral body. Since the hardness of the system depends on the amount and location of the injected bone cement, such a procedure can be optimized by knowing the structural / functional properties of the bone-bone cement system obtained by the procedure. Also, the treatment of fractures using metal or other types of prostheses can be improved by evaluating the structural response of the bone grafting system, including factors such as size and shape, prosthesis material, body part density and structure, etc. . This allows the physician to assess the suitability of the proposed surgical procedure prior to the procedure and, as a result, select the optimal procedure for each individual patient. At present, most physicians rely only on qualitative experience and do not have quantitative criteria to evaluate various options. Various types of knowledge regarding the structural / functional properties of the bone grafting system are expected to improve the outcome of the treatment.
本発明は、人間又は動物の生体の構造の画像から人間又は動物の生体の構造の１又は複数の構造／機能特性を判定する特性判定方法を提供する。特性判定方法は、人間又は動物の生体の構造の画像を受信するステップを有する。そして、この画像に基づいて、構造の構造モデルを生成する。そして、構造モデルに基づいて、第１の生体力学的量を算出する。そして、この構造モデルを変更し、変形モデルを生成する。そして、変形モデルに基づいて、第２の生体力学的量を算出する。そして、第１及び第２の生体力学的量を比較する。比較の結果は、デジタルメディアに保存される。 The present invention provides a characteristic determination method for determining one or more structural / functional characteristics of a human or animal biological structure from an image of the human or animal biological structure. The characteristic determination method includes a step of receiving an image of a human or animal living body structure. Then, a structural model of the structure is generated based on this image. Then, a first biomechanical quantity is calculated based on the structural model. Then, the structural model is changed to generate a deformation model. Then, a second biomechanical quantity is calculated based on the deformation model. Then, the first and second biomechanical quantities are compared. The result of the comparison is stored in digital media.
特性判定方法は、更に、比較に基づいて、１又は複数の構造／機能特性を判定するステップを有していてもよい。特性判定方法は、更に、負荷条件を変更してもよく、特性判定方法を後に繰り返し、構造の１又は複数の構造／機能特性に関して、治療、老化又は病気若しくはこれらの組合せの１又は複数の作用を判定してもよい。 The characteristic determination method may further include determining one or more structural / functional characteristics based on the comparison. The characterization method may further change the loading conditions, and the characterization method may be repeated later to affect one or more effects of treatment, aging or disease or a combination thereof with respect to one or more structure / function properties of the structure. May be determined.
この構造は、例えば、骨等の筋骨格組織又は器官、又は例えば、心臓又は血管等の心臓血管組織又は器官を含んでいてもよく、これらは移植された部位を有していても、有していなくてもよい。変形モデルは、均質化されたモデルを含んでいてもよく、特性判定方法は、構造モデルの１又は複数の要素に平均密度を割り当てるステップを更に有していてもよい。 This structure may include, for example, a musculoskeletal tissue or organ such as a bone, or a cardiovascular tissue or organ such as a heart or blood vessel, which may or may not have an implanted site. It does not have to be. The deformation model may include a homogenized model, and the characterization method may further comprise assigning an average density to one or more elements of the structural model.
変形モデルは、基準モデルを含んでいてもよく、特性判定方法は、構造モデルの１又は複数の要素に基準密度を割り当てるステップを更に有していてもよい。変形モデルは、副構造モデル（sub-structure model）を含んでいてもよく、特性判定方法は、構造モデルから構造の一部を取り除くステップを更に有していてもよい。取り除かれる構造の一部は、周辺構造を含んでいてもよい。 The deformation model may include a reference model, and the characteristic determination method may further include assigning a reference density to one or more elements of the structural model. The deformation model may include a sub-structure model, and the characteristic determination method may further include a step of removing a part of the structure from the structure model. The portion of the structure that is removed may include a peripheral structure.
また、変形モデルは、軸モデル又は湾曲モデルを含んでいてもよい。変形モデルは、均質化されたモデル、副構造モデル、軸モデル、湾曲モデル、及び基準モデルの２つ以上の組合せを含んでいてもよい。 The deformation model may include an axis model or a curved model. The deformation model may include a combination of two or more of a homogenized model, a substructure model, an axis model, a curved model, and a reference model.
変形モデルは、境界条件が変更された構造モデルの変形を含んでいてもよい。境界条件は、力、圧力、曲げモーメント、変形、偏位、速度、加速度、流量、温度、エネルギ、歪み又は応力又はこれらの組合せを含んでいてもよい。 The deformation model may include a deformation of the structural model whose boundary condition is changed. Boundary conditions may include force, pressure, bending moment, deformation, excursion, velocity, acceleration, flow rate, temperature, energy, strain or stress, or combinations thereof.
特性判定方法は、更に、構造をスキャニングし、構造の画像を生成するステップを更に有していてもよい。スキャニングは、コンピュータ断層撮影（computed tomography：ＣＴ）及び磁気共鳴撮影（magnetic resonance imageＭＲＩ）スキャンを含んでいてもよく、構造モデルは、構造の有限要素モデルを含んでいてもよい。また、特性判定方法は、第１の画像が取得された時刻とは異なる時刻に構造の第２の画像を受信するステップを有していてもよい。第２の画像に基づいて、構造の第２の構造モデルを生成してもよい。第２の構造モデルに基づいて、第３の生体力学的量を算出してもよい。第２の構造モデルを変更して第２の変形モデルを生成してもよい。そして、第２の変形モデルに基づいて、第４の生体力学的量を算出し、比較は、第３及び第４の生体力学的量の比較を含んでいてもよい。比較は、第１及び第２の生体力学的量の比較の結果と、第３及び第４の生体力学的量の比較の結果との比較を含んでいてもよい。 The characteristic determination method may further include a step of scanning the structure and generating an image of the structure. Scanning may include computed tomography (CT) and magnetic resonance image MRI scans, and the structural model may include a finite element model of the structure. The characteristic determination method may include a step of receiving the second image having the structure at a time different from the time at which the first image is acquired. A second structural model of the structure may be generated based on the second image. A third biomechanical quantity may be calculated based on the second structural model. The second structural model may be changed to generate the second deformation model. Then, a fourth biomechanical quantity is calculated based on the second deformation model, and the comparison may include a comparison between the third and fourth biomechanical quantities. The comparison may include comparing the result of the comparison of the first and second biomechanical quantities with the result of the comparison of the third and fourth biomechanical quantities.
特性判定方法は、構造とは異なる生体の部分の第２の画像を受信するステップを有していてもよい。第２の画像に基づいて、第２の構造モデルを生成してもよい。第２の構造モデルに基づいて、第３の生体力学的量を算出してもよい。第２の構造モデルを変更して第２の変形モデルを生成してもよい。第２の変形モデルに基づいて、第４の生体力学的量を算出し、比較は、第３及び第４の生体力学的量の比較を含んでいてもよい。比較は、第１及び第２の生体力学的量の比較の結果と、第３及び第４の生体力学的量の比較の結果との比較を含んでいてもよい。 The characteristic determination method may include a step of receiving a second image of a part of the living body different from the structure. A second structural model may be generated based on the second image. A third biomechanical quantity may be calculated based on the second structural model. The second structural model may be changed to generate the second deformation model. A fourth biomechanical quantity is calculated based on the second deformation model, and the comparison may include a comparison of the third and fourth biomechanical quantities. The comparison may include comparing the result of the comparison of the first and second biomechanical quantities with the result of the comparison of the third and fourth biomechanical quantities.
本発明は、他の実施の形態として、人間又は動物の生体の構造に対する治療の作用を判定する評価方法を提供する。評価方法は、第１の期間において、生体の構造の画像を受信するステップと、画像に基づいて、構造の構造モデルを生成するステップと、構造モデルに基づいて、第１の生体力学的量を算出するステップと、構造モデルを変更して、変形モデルを生成するステップと、変形モデルに基づいて、第２の生体力学的量を算出するステップと、第１及び第２の生体力学的量を比較するステップと、第１及び第２の生体力学的量の比較に基づいて、１又は複数の構造／機能特性を判定するステップと、第２の期間に取得された構造の第２の画像を受信するステップと、第２の画像に基づいて、構造の第２の構造モデルを生成するステップと、第２の構造モデルに基づいて、第３の生体力学的量を算出するステップと、第２の構造モデルを変更して第２の変形モデルを生成するステップと、第２の変形モデルに基づいて、第４の生体力学的量を算出するステップと、第３及び第４の生体力学的量を比較するステップと、第１及び第２の生体力学的量の比較の結果と、第３及び第４の生体力学的量の比較の結果とを比較し、構造の老化、病気又は治療に関する作用を判定するステップとを有する。 As another embodiment, the present invention provides an evaluation method for determining the effect of treatment on the structure of a human or animal living body. In the first period, the evaluation method includes a step of receiving an image of a structure of a living body, a step of generating a structural model of the structure based on the image, and a first biomechanical quantity based on the structural model. A step of calculating, a step of generating a deformation model by changing the structural model, a step of calculating the second biomechanical quantity based on the deformation model, and the first and second biomechanical quantities A step of comparing, a step of determining one or more structural / functional properties based on a comparison of the first and second biomechanical quantities, and a second image of the structure acquired during the second period Receiving, generating a second structural model of the structure based on the second image, calculating a third biomechanical quantity based on the second structural model, and second Change the structural model of the second Generating a deformation model; calculating a fourth biomechanical quantity based on the second deformation model; comparing the third and fourth biomechanical quantities; Comparing the results of the comparison of the two biomechanical quantities with the results of the comparison of the third and fourth biomechanical quantities to determine the effect on structural aging, disease or treatment.
本発明に基づく方法は、人間又は動物の治療に適用できる。また、本発明は、生体内及び生体外でのスキャニングにも適用できる。臨床検査は、人間の生体の検査のみを目的としているが、本発明は、動物及び死体についても、治療、老化及び病気の影響を調べることができる。 The method according to the invention can be applied to the treatment of humans or animals. The present invention can also be applied to scanning in vivo and in vitro. Although the clinical examination is intended only for the examination of a human living body, the present invention can examine the effects of treatment, aging, and disease on animals and corpses.
また、本発明は、プロセッサによって読み込み可能なコードを有する１つ以上のプロセッサによって読出可能な記録装置を提供する。このプロセッサによって読出可能な記録装置は、上述した何れかの方法を実現する。 The present invention also provides a recording device readable by one or more processors having code readable by the processor. The recording apparatus readable by this processor realizes any of the methods described above.
以下では、当業者が本発明を理解できるように、本発明者が想定する本発明の最良の実施の形態について説明する。なお、当業者は、以下の実施の形態の様々な変形例を想到できる。これらの変形例、等価物及び代替物は、本発明の精神及び範囲内にある。 In the following, the best mode of the present invention assumed by the present inventor will be described so that those skilled in the art can understand the present invention. A person skilled in the art can conceive various modifications of the following embodiments. These variations, equivalents and alternatives are within the spirit and scope of the present invention.
骨折のリスクが高い患者を予測し、及び治療の効能を評価するには、骨密度のＤＸＡ測定値では限界があり、したがって、年齢、病気及び治療によって、「骨の品質」の如何なる側面が変化したかを定量化する手法の実現が望まれている。本発明は、このような骨全体の品質の定量化を実現する。本発明により、骨粗鬆症の診断、治療の経過の監視、新たな治療の評価が改善され、したがって、骨粗鬆症の患者にとって有益である。 Predicting patients at high risk of fractures and assessing treatment efficacy is limited by DXA measurements of bone density, so any aspect of “bone quality” varies with age, disease, and treatment The realization of a method for quantifying whether or not it is done is desired. The present invention realizes such quantification of the quality of the whole bone. The present invention improves the diagnosis of osteoporosis, the monitoring of the course of treatment, the evaluation of new therapies and is therefore beneficial for osteoporotic patients.
本明細書では、骨の品質を定義するために「構造／機能特性（structure-function characteristics）」という用語を用いる。「構造／機能特性」とは、撮像された骨全体の構造の測定値と、骨全体の何らかの生体力学的振る舞いとの間のあらゆる関係を意味する。また、本明細書で用いる「構造」という用語には、人間又は動物の生体の任意の部分、例えば、器官、骨又は心臓又は椎間板又は関節、又は腰、膝関節、又はこれらの任意の器官の一部が含まれると解釈される。また、構造には、例えば、骨梁、軟骨組織、椎間板、靭帯、腱、血管等の組織、又は例えば、骨髄、血液、又は他の何らかの体液等の流体成分を含む生体の何らかの組織も含まれる。したがって、本明細書で用いる「構造」という用語は、任意の物理的寸法を有する生体内のあらゆる部分が含まれる。また、構造は、生体の一部に移植されたインプラントを含んでいてもよく、例えば、骨移植システム大腿骨近位部に移植された人工の臀部インプラント、又は血管内に移植された心臓血管ステント等を含んでいてもよい。構造の形状、形態、質量又は密度を特徴付ける如何なる定量的測定値用いて構造を定量化し、構造／機能関係を特徴付けてもよい、例えば、体重、伸長、骨密度骨密度分布（骨内の骨密度値の空間的分布の標準偏差によって特徴付けられる）、臀部の軸の長さ（hip-axis length）、血管の直径、血管の軸方向に沿った厚さの変化（血管の軸方向に沿った厚さ値の標準偏差によって特徴付けられる）、軟骨組織厚、筋肉の向きと長さ、筋肉の断面積、骨梁組織の数、骨梁組織の接続、又は体液の粘度等の測定値によって構造を特徴付けてもよい。 The term “structure-function characteristics” is used herein to define bone quality. “Structural / functional properties” means any relationship between the measured bone structure measurements and any biomechanical behavior of the entire bone. Also, as used herein, the term “structure” includes any part of a human or animal body, such as an organ, bone or heart or intervertebral disc or joint, or hip, knee joint, or any of these organs. It is interpreted as including a part. Structures also include tissues such as trabecular bone, cartilage tissue, intervertebral discs, ligaments, tendons, blood vessels, or any tissue of a living body that includes fluid components such as bone marrow, blood, or some other body fluid. . Accordingly, as used herein, the term “structure” includes any part of a living body having any physical dimension. The structure may also include an implant implanted in a part of the body, for example, an artificial hip implant implanted in the bone graft system proximal femur, or a cardiovascular stent implanted in a blood vessel Etc. may be included. Any quantitative measurement that characterizes the shape, shape, mass or density of the structure may be used to quantify the structure and characterize the structure / function relationship, eg, weight, elongation, bone density, bone density distribution (bones within bone (Characterized by the standard deviation of the spatial distribution of density values), hip-axis length, vessel diameter, thickness change along the vessel axis (along the vessel axis) By cartilage tissue thickness, muscle orientation and length, muscle cross-sectional area, number of trabecular tissue, trabecular tissue connections, or fluid viscosity, etc. The structure may be characterized.
また、「機能」という用語は、画像の生体力学的解析によって算出された、構造の何らかの生体力学的特徴を指す。「機能」には、以下に限定されるものではないが、例えば、強さ、硬さ、耐疲労性、破砕強度、強度、変形、破砕のリスク、流体せん断応力及び圧力等の機械的パラメータが含まれる。温度及び熱輸送に関連する生体力学的問題では、温度勾配、熱流束、及び熱膨脹ストレス等のパラメータも「機能」の定量的測定値となる。 Also, the term “function” refers to some biomechanical feature of the structure calculated by biomechanical analysis of the image. “Function” is not limited to the following, but includes mechanical parameters such as strength, hardness, fatigue resistance, crushing strength, strength, deformation, risk of crushing, fluid shear stress and pressure. included. In biomechanical problems related to temperature and heat transport, parameters such as temperature gradient, heat flux, and thermal expansion stress are also quantitative measurements of “function”.
本発明の「構造／機能特性」は、構造の定量的測定値と、機能の定量的測定値との間の何らかの数学的関係を意味する。例えば、複数の骨の強度値と密度値の間のプロットである密度／強度の関係は、骨について知られている構造／機能関係である。単一の骨又は骨組織の一片の密度に対する強度の比率は、構造／機能特性の別の表現である。構造／機能特性の他の具体例として、血管の直径と、血管壁のせん断応力との間の関係がある。 By “structural / functional property” of the present invention is meant some mathematical relationship between a quantitative measurement of structure and a quantitative measurement of function. For example, a density / strength relationship that is a plot between strength values and density values of a plurality of bones is a known structure / function relationship for bone. The ratio of strength to density of a single bone or piece of bone tissue is another representation of structural / functional properties. Another example of structural / functional properties is the relationship between vessel diameter and vessel wall shear stress.
更に、本発明の「構造／機能特性」は、構造の画像のモデルの変化を用いて算出された何らかの生体力学的パラメータ間の又はａ）構造の画像のモデルの変化を用いて算出された他の何らかの生体力学的パラメータ又はｂ）例えば、上述したような構造の定量的測定値又は機能の定量的測定値の何れかの任意の組合せ間の何らかの数学的関係を示す。例えば、人間の椎体のモデルから２ｍｍの周囲の骨を取り除いた後の骨梁組織の強度と、（完全な）椎体の強度との比率も、構造／機能特性であるとみなすことができる。 Furthermore, the “structural / functional properties” of the present invention may be between any biomechanical parameters calculated using a structural image model change or a) other using a structural image model change. B) any mathematical relationship between any combination of any of the biomechanical parameters of or b), for example, quantitative measurements of structure or quantitative measurements of function as described above. For example, the ratio of trabecular tissue strength to (complete) vertebral body strength after removing 2 mm of surrounding bone from a human vertebral body model can also be considered a structural / functional property. .
現在も多くの新薬が開発されているが、本発明の好ましい実施の形態に基づく評価ツールは、この開発に著しく有益である。本明細書では、このようなツールを提供する手法を開示する。一実施の形態においては、情報を臨床的に用いて、個々の患者のための治療を最適化することができる。現在、異なる生物作用を有する多くの骨粗鬆症薬物療法が市販されている（例えば、アレンドロネートは、骨吸収を抑制し、ＰＴＨは、新たな骨形成を促進する）。より多くの治療法が開発されると、患者に投与できる薬剤の選択肢も多くなる。本明細書は、臨床医が、特定の患者について、各薬物療法を最適化できるツールを提供する。 Although many new drugs are still being developed, evaluation tools based on preferred embodiments of the present invention are significantly beneficial for this development. This document discloses a technique for providing such a tool. In one embodiment, the information can be used clinically to optimize treatment for individual patients. Currently, many osteoporosis drug therapies with different biological effects are on the market (eg, alendronate suppresses bone resorption and PTH promotes new bone formation). As more therapies are developed, there are more drug options available to the patient. This document provides a tool that allows the clinician to optimize each drug therapy for a particular patient.
本発明は、骨の強度の生体力学的構造／機能特性を定量化し、骨粗鬆症等の病気の診断、観察、評価及び治療の一部にこの情報を利用するための様々な手法を提供する。これは、パラメータを制御しながら、条件を様々に変更して、骨の強度又は他の生体力学的性質を定量化し、患者の骨を調べることによって実現される。好ましい実施の形態では、個々の患者について管理された手法で骨又は骨データを変更し、経時的な骨粗鬆症の診断、又は治療法の評価のための患者固有のデータを生成する。正常な、健康な骨に由来するデータベースに対して、又は同じ骨の以前の振る舞いに対して、骨の構造／機能特性の比較を行うことにより、骨量、密度又は強度等の単一のメトリックを用いた場合に比べて、病気の評価及び治療効果について、より詳細で有益な情報が提供される。骨の構造／機能特性の変化は、特定の時刻において定量化でき（横の比較）及び経時的にも定量化でき（縦の比較）、特に経時的な定量化は、骨の強度の治療法、又は病気又は老化の進行に関連する生体力学的メカニズムを理解するために有益である。 The present invention quantifies the biomechanical structural / functional properties of bone strength and provides a variety of techniques for utilizing this information as part of the diagnosis, observation, evaluation and treatment of diseases such as osteoporosis. This is accomplished by examining the patient's bone by quantifying bone strength or other biomechanical properties while controlling the parameters and varying the conditions. In a preferred embodiment, bone or bone data is modified in a controlled manner for individual patients to generate patient-specific data for diagnosis of osteoporosis over time or evaluation of treatment. A single metric, such as bone mass, density or strength, by comparing bone structural / functional properties against a database derived from normal, healthy bones, or against previous behavior of the same bone More detailed and useful information is provided about disease assessment and treatment effectiveness compared to using. Changes in bone structural / functional properties can be quantified at a specific time (horizontal comparison) and can also be quantified over time (longitudinal comparison), especially quantification over time is a treatment of bone strength. Or to understand the biomechanical mechanisms associated with the progression of illness or aging.
本発明の好ましい実施の形態に基づく手法は、単なる骨の強度及び骨粗鬆症の検査に比べて汎用性が高い。この手法は、骨の強度に関連する応用例、及び組織又は器官の機械的挙動が臨床的に重要な他の組織及び器官に関連する応用例に適用可能である。例えば、心血管疾患では、血管の構造／機能関係の分析に基づいて、石灰化したプラーク又は血管破裂の存否によって凝固のリスクを定量化し、心臓病の診断を改善させることができる。関節炎では、軟骨組織の生体力学的解析によって、外科的施術及び薬物療法の応答の予測が向上する。また、本発明によって、患者の病気の観察及び治療の監視も向上する。また、本発明により、人工関節、若しくはバイオマテリアル又は組織再生、又は何らかの物質の移植に対する生体の生体力学的応答も詳細に観察でき、特定の患者に関して、移植されたアイテムの適合性等を検証することができる。 The technique according to the preferred embodiment of the present invention is more versatile than simple bone strength and osteoporosis testing. This approach is applicable to applications related to bone strength and to other tissues and organs where the mechanical behavior of the tissue or organ is clinically important. For example, in cardiovascular disease, the risk of coagulation can be quantified by the presence or absence of calcified plaques or vascular rupture based on the analysis of vascular structure / function relationships to improve the diagnosis of heart disease. In arthritis, biomechanical analysis of cartilage tissue improves the prediction of response to surgical procedures and drug therapy. The present invention also improves patient disease monitoring and treatment monitoring. In addition, according to the present invention, it is possible to observe in detail the biomechanical response of a living body to an artificial joint, or biomaterial or tissue regeneration, or transplantation of any substance, and verify the suitability of the implanted item with respect to a specific patient. be able to.
患者毎に、非観血的に骨の強度を評価する技術及びモデリングが望まれている。このようなモデルに関して、好ましくは、パラメータを制御した検査を用いて、更なるデータを取得し、これを用いて、患者に応じた手法で、治療を診断又は評価し、又は病気の進行を追跡し、又は組織又は器官の何らかの変化を追跡する。本発明の包括的な概念は、骨梁等の組織及び大腿骨近位部又は椎体等の器官の両方に適用できる。また、本発明は、患者特有の画像を用いて、非観血的な手法で機械的挙動を評価する生体力学的システムに適用できる。 A technique and modeling for evaluating bone strength noninvasively for each patient is desired. For such models, preferably, parameter-controlled testing is used to acquire further data, which can be used to diagnose or evaluate therapy or track disease progression in a patient-dependent manner. Or track any changes in tissue or organs. The generic concept of the present invention can be applied to both tissues such as trabecular bone and organs such as proximal femur or vertebral body. In addition, the present invention can be applied to a biomechanical system that evaluates mechanical behavior by a noninvasive method using patient-specific images.
具体的には、例えば、骨、軟骨組織、腱又は他の筋骨格組織、又は血管等の生物学的構造又は物質に関する患者特有の有限要素モデルに関してパラメータ研究を行い、生体力学的構造／機能特性を特徴付ける手法を開示する。例えば、以下に限定されるものではないが、ＭＲＩ、ＣＴ、超音波、ＱＣＴ、マイクロＣＴ、マイクロＭＲＩ、有限要素モデリング等の患者固有の画像を取得する３次元又は２次元の撮影の任意の組合せ、又はビーム理論又は連続体力学の原理を用いた他の解析モデリング、又は機械的応答を生成する材料力学、流体力学、熱輸送、質量輸送又は熱力学等、如何なる非観血的な技術を用いて機械的挙動を評価してもよい。 Specifically, parametric studies are performed on patient-specific finite element models for biological structures or materials such as bone, cartilage tissue, tendons or other musculoskeletal tissues, or blood vessels, and biomechanical structural / functional properties. Disclose a technique for characterizing. For example, but not limited to, any combination of 3D or 2D imaging to acquire patient specific images such as MRI, CT, ultrasound, QCT, micro CT, micro MRI, finite element modeling, etc. , Or other analytical modeling using beam theory or continuum mechanics principles, or any non-invasive techniques such as material mechanics, fluid mechanics, heat transport, mass transport or thermodynamics that generate a mechanical response The mechanical behavior may be evaluated.
例えば、骨全体の分析については、本発明に基づく好適な技術は、骨全体の様々な要素が強度にどのように貢献するかについて、量的な測定値を提供する。これらの関係のうちで非常に強い相関関係の存在（図２参照）により、骨の強度が正常に見えても、正常な骨の構造／機能挙動の典型的な範囲内に収まらない骨を特定することができる。この手法は、単一の強度メトリック（圧縮強度等）の使用に比べて、骨の生体力学的状態について、より正確な尺度を提供し、骨全体の構造／機能特性を変更する老化、病気、及び薬物療法又は他の治療の生体力学的効果の評価について、追加的且つ固有の洞察を行うことができる。この分析は、所定の時点で行ってもよく、経時的に繰り返し行ってもよい。 For example, for whole bone analysis, the preferred technique according to the present invention provides a quantitative measure of how various elements of the whole bone contribute to strength. Of these relationships, the presence of a very strong correlation (see Figure 2) identifies bones that do not fall within the typical range of normal bone structural / functional behavior, even though bone strength appears normal can do. This approach provides a more accurate measure of the biomechanical state of bone compared to the use of a single strength metric (such as compressive strength) and aging, illness, And additional and specific insights can be made on the assessment of the biomechanical effects of drug therapy or other treatments. This analysis may be performed at a predetermined time point or may be repeated over time.
本発明に用いる有限要素モデリング又は他の何らかの機械的分析の患者特有の性質は、組織又は器官に適用される非観血的な撮像によるこのようなモデルの生成によって得られ、このような撮像には、例えば、以下に限定されるものではないが、コンピュータ断層撮影（computed tomography：ＣＴ）、磁気共鳴画像（magnetic resonance imaging：ＭＲＩ）、二重エネルギＸ線吸収（dual-energy x-ray absorptiometry：ＤＸＡ）、ポジトロン放射断層撮影（positron emission tomography：ＰＥＴ）、マイクロＣＴ、周辺ＣＴ、定量的ＣＴ（quantitative CT：ＱＣＴ）、又は患者の組織又は器官の２次元画像又は３次元画像を生成する何らかのスキャン又はスキャンの組合せが含まれる。組織又は器官のコンピュータモデルが生成された後、例えば、典型的には、硬度又は強度等、関連する機械的特性が計算又は算出される。また、スキャンを分析することによって、組織又は器官の密度及び／又は形状又は形態の様々な側面を定量化できる。 The patient-specific nature of finite element modeling or some other mechanical analysis used in the present invention is obtained by the generation of such a model by non-invasive imaging applied to a tissue or organ, and For example, but not limited to, computed tomography (CT), magnetic resonance imaging (MRI), dual-energy x-ray absorptiometry: DXA), positron emission tomography (PET), micro CT, peripheral CT, quantitative CT (QCT), or any scan that produces a 2D or 3D image of a patient's tissue or organ Or a combination of scans is included. After the computer model of the tissue or organ is generated, associated mechanical properties, such as typically hardness or strength, are typically calculated or calculated. Also, by analyzing the scan, various aspects of tissue or organ density and / or shape or morphology can be quantified.
そして、様々な１つ以上の制御された手法を用いてコンピュータモデルを変更し、これにより、機械的特性値を算出する。相互に及び／又は密度及び形状プロパティに対して機械的特性をプロットすることによって様々な構造／機能特性が定量化される。これらのある数学的関数、例えば、密度と強度との比率、又は導出された１つの機械的特性と、他の特性、強度又は密度との比率も構造／機能特性の要素を表す。主な機械的特性（例えば、強度又は硬度）に加えて、組織／器官状態を評価するために構造／機能特性を用いることも有益である。これにより、元のモデルに関して、パラメータを制御しながら、患者固有の医療画像から新たなデータを取得し、これらのデータは、病気の診断、治療、変化、又は老化した状態の評価に用いることができる。 The computer model is then modified using various one or more controlled techniques, thereby calculating mechanical property values. Various structural / functional properties are quantified by plotting mechanical properties against each other and / or density and shape properties. These mathematical functions, such as the ratio of density to strength, or one derived mechanical property to the ratio of other properties, strength or density, also represent elements of structural / functional properties. In addition to the main mechanical properties (eg strength or hardness), it is also beneficial to use structural / functional properties to assess tissue / organ status. This allows new data to be acquired from patient-specific medical images while controlling parameters with respect to the original model, and these data can be used to diagnose, treat, change, or evaluate aging conditions. it can.
例えば、骨梁のマイクロ構造に関して骨粗鬆症薬物療法の効果を評価する際に、患者の手首の高解像度マイクロＣＴスキャンを用いた有限要素モデリングを実行することができる。コンピュータモデルを用いて患者の骨に関してパラメータが制御された研究を実行することによって、幾何学的なマイクロ構造における石灰化の分布におけるあらゆる変化の個別の役割を抽出することができる。これにより、薬物療法の生体力学的作用を更に深く調べることができる。同様に、患者の心臓の臨床的ＣＴスキャンから得られた情報を用いて、異なる石灰化の度合いの下で、血管の物理的一体性を評価することもできる。 For example, in evaluating the effects of osteoporosis medication on the trabecular microstructure, finite element modeling using a high resolution micro CT scan of the patient's wrist can be performed. By performing a parameter-controlled study on the patient's bone using a computer model, the individual role of any change in the distribution of calcification in the geometric microstructure can be extracted. Thereby, the biomechanical effect | action of a drug therapy can be investigated further deeply. Similarly, information obtained from clinical CT scans of the patient's heart can be used to assess the physical integrity of blood vessels under different degrees of calcification.
以下では、薬物療法の評価のうち、椎体の解析に適用される処理の詳細な具体例を説明する。なお、本発明は、このような特定用途に限定されるわけではない。 Below, the detailed specific example of the process applied to the analysis of a vertebral body is demonstrated among evaluation of a pharmacotherapy. Note that the present invention is not limited to such a specific application.
まず、患者をＱＣＴスキャンし、患者の椎体のデジタル画像を生成する。次に、例えばクロフォード（Crawford）他の論文に開示されている手続きに従って、脊椎骨又は例えば、Ｌ１又はＴ９椎体等、対象となる椎体のＱＣＴ画像を有限要素モデルに変換する。簡潔に言えば、スキャンの形状情報を用いて、骨の有限要素メッシュが生成される。そして、スキャンの骨密度情報を用いて、要素毎に、モデル内の各有限要素の材料特性が算出される。骨組織の弾性及び強度特性は、モデル内の各骨要素に割り当てられる（Crawford Trans Orthopaedic Research Society 2004; or Faulkner et al..: "Effect of bone distribution on vertebral strength: assessment with patient-specific nonlinear finite element analysis." Radiology 179:669-674, 1991）。適用された負荷に関連する有限要素モデル境界条件（finite element model boundary conditions）がモデルに適用され、コンピュータがこのモデルを解決し、強度又は対象となる他の構造特性、例えば、硬度又は変形を判定する。 First, the patient is QCT scanned to generate a digital image of the patient's vertebral body. Next, the QCT image of the vertebrae of interest, such as the L1 or T9 vertebral bodies, is converted into a finite element model, for example according to the procedure disclosed in Crawford et al. Briefly, a finite element mesh of bone is generated using scan shape information. Then, the material properties of each finite element in the model are calculated for each element using the bone density information of the scan. The elastic and strength properties of bone tissue are assigned to each bone element in the model (Crawford Trans Orthopaedic Research Society 2004; or Faulkner et al ..: "Effect of bone distribution on vertebral strength: assessment with patient-specific nonlinear finite element analysis. "Radiology 179: 669-674, 1991). Finite element model boundary conditions related to the applied load are applied to the model and the computer solves the model to determine strength or other structural characteristics of interest, such as hardness or deformation To do.
次に、元のモデルに関して、複数のパラメータ研究を行い、モデルを変更し、これにより得られる強度の値（又は骨全体の機械的挙動に関するパラメータ）を算出する。図１は、このようなパラメータ研究の具体例を示しており、この場合、７つの個々の有限要素解析を実行し、５つの比率を含む合計１３個の結果パラメータを算出している。 Next, a plurality of parameter studies are performed on the original model, the model is changed, and a strength value (or a parameter related to the mechanical behavior of the entire bone) obtained thereby is calculated. FIG. 1 shows a specific example of such a parameter study where seven individual finite element analyzes are performed and a total of 13 result parameters including five ratios are calculated.
このような分析の１つでは、外側の２ｍｍの骨を取り除き、これにより得られる骨梁組織部の強度を算出する。これにより、骨梁組織部の構造的な役割が定量化される（図１の変数ｔｒａｂ）。図２は、女性の健康なグループの椎体の構造／機能関係の具体例を示しており、ここでは、強度値は、被験者毎のＦＥＭによって算出され、正常な人間の脊椎骨のグループについて、周囲の骨が２ｍｍ浸食された場合の強度を示している。このグラフは、１３個の健康な椎体のサンプルについて、骨梁部の強度対椎体全体の強度をプロットしたものである。ここでは、老化、病気又は治療のために経時的に生じる個々の振る舞いの変化について、３個の仮定的状況を示している。この関係において、高いＲ２値（Ｒ２＝０．９９）は、この関係からの比較的小さい偏りが、正常な構造／機能挙動から外れる振る舞いを示すことを表している。このような偏りは、骨全体の「品質」の変化として解釈でき、この情報は、臨床的に有用である。 In one such analysis, the outer 2 mm bone is removed and the strength of the trabecular tissue obtained thereby is calculated. Thereby, the structural role of the trabecular tissue part is quantified (variable trab in FIG. 1). FIG. 2 shows an example of the structure / function relationship of the vertebral bodies of a female healthy group, where intensity values are calculated by FEM for each subject, and for a group of normal human vertebrae, The strength when 2 mm of bone was eroded by 2 mm is shown. The graph is a plot of trabecular strength versus total vertebral strength for 13 healthy vertebral body samples. Here, three hypothetical situations are shown for individual behavioral changes that occur over time due to aging, disease or treatment. In this relationship, a high R 2 value (R 2 = 0.99) indicates that a relatively small deviation from this relationship indicates behavior that deviates from normal structural / functional behavior. Such a bias can be interpreted as a change in the “quality” of the entire bone, and this information is clinically useful.
また、図２は、老化、病気又は治療のため生じる構造／機能特性の変化について、３つの可能な仮定的状況を示している。このように変化を分析することによって、変化の生体力学的メカニズムを理解し、骨全体の品質の変化を洞察し、実際に定量化することができる。まず、初期の時点であるポイントＡから、個々のサンプルは、３つのポイントに移行できる。ポイントＣは、「正常な」構造／機能挙動の回帰直線上にプロットされているため、品質の変化はない。ポイントＤは、回帰直線の下に位置するため、骨梁組織部の強度は、比較的低く、構造／機能挙動が変化したことを意味する。ポイントＢでは、骨梁組織部の強度は比較的高く、これも、構造／機能挙動の変化を意味する。何れの場合も、椎体全体の強度の増加量は同じであり、したがって、何れの場合も、脊椎の強度の最終的な値は同じである。例えば、薬物研究において、ポイントＢの結果が得られた場合、薬物療法によって、骨梁組織部の強度が好適に高められ、総合的な強度増進効果が達成されたことが確認される。この情報は、薬物療法の生体力学的作用を明らかにするために役立つ。 FIG. 2 also shows three possible hypothetical situations for structural / functional property changes that occur due to aging, disease or treatment. By analyzing changes in this way, it is possible to understand the biomechanical mechanism of changes, gain insight into changes in the quality of the entire bone, and actually quantify them. First, from point A, which is the initial time point, each sample can move to three points. Point C is plotted on a “normal” structural / functional behavior regression line, so there is no change in quality. Since the point D is located below the regression line, the strength of the trabecular tissue is relatively low, which means that the structure / functional behavior has changed. At point B, the strength of the trabecular tissue is relatively high, which also means a change in structural / functional behavior. In either case, the amount of increase in strength of the entire vertebral body is the same, so in each case the final value of spine strength is the same. For example, when the result of point B is obtained in the drug research, it is confirmed that the strength of the trabecular tissue part is suitably increased by the drug therapy and the overall strength enhancement effect is achieved. This information helps to clarify the biomechanical effects of drug therapy.
このモデルを変更して、副構造モデルが生成され、この副構造モデルから、必ずしも周囲の骨でなくてもよいが、骨のある一部を取り除く。例えば、椎体のモデルでは、骨梁の内核（inner core）を取り除き、取り除かれる内核の直径を順次大きくし、骨全体から取り除かれる内部の骨の量を変化させる。 The model is modified to generate a substructure model, from which a portion of the bone is removed, although not necessarily the surrounding bone. For example, in a vertebral body model, the inner core of the trabecular bone is removed, the diameter of the inner core removed is sequentially increased, and the amount of internal bone removed from the entire bone is changed.
モデルの他の変形例では、モデル内の骨材料の一部又は全部に、その椎体の平均密度を割り当てる。これにより、骨が均質化され、脊椎内の密度の変化の一部又は全部が除外される。この結果、椎体の骨材料の材料特性も均質化され、椎体全体について同じ値が使用される。そして、この新たなモデルの強度（図１のｈｏｍ）が算出される。この強度値と、以前に完全な椎体について算出した強度との間の比率（図１のｓｔｄ）は、構造／機能特性の量的な側面として、脊椎骨内の密度の変化の重要度の測定値を提供する。通常の椎体の場合ｓｔｄ／ｈｏｍ比は、１未満である（図３）。密度の変化の一部又は全部を無視するために、均質化の度合いを変更することもできる。 In other variations of the model, some or all of the bone material in the model is assigned an average density of the vertebral body. This homogenizes the bone and excludes some or all of the density changes in the spine. As a result, the material properties of the bone material of the vertebral body are also homogenized and the same value is used for the entire vertebral body. Then, the strength of the new model (hom in FIG. 1) is calculated. The ratio between this strength value and the strength previously calculated for the complete vertebral body (std in FIG. 1) is a measure of the importance of density changes in the vertebra as a quantitative aspect of the structural / functional properties. Provide a value. For normal vertebral bodies, the std / hom ratio is less than 1 (FIG. 3). The degree of homogenization can be changed to ignore some or all of the density change.
骨がより均質になると、この比率の値は１に近付き、骨は、構造的な見地から、より最適になる。したがって、１に近い値は、構造的に効率的な骨を示し、より低い値は、最適ではない骨を示す。薬物療法によって、この比率が変化することがあり、これは、治療による骨の総合的な構造／機能特徴が変化したことを意味する。例えば、治療によってこの比率が高まった場合、これは、治療によって、骨の構造的効率性が高まり、すなわち、所定の量の骨量の強度が強くなったことを意味する。理論的には、骨の強度が弱くなっても、骨の構造的効率性が高くなる場合もある。したがって、ｓｔｄ／ｈｏｍ比は、単なる強度の測定値より優れた追加的情報を提供する。 As the bone becomes more homogeneous, the value of this ratio approaches 1 and the bone becomes more optimal from a structural point of view. Thus, values close to 1 indicate structurally efficient bones and lower values indicate non-optimal bones. Drug therapy can change this ratio, which means that the overall structural / functional characteristics of the bone due to treatment have changed. For example, if this ratio is increased by treatment, this means that the treatment has increased the structural efficiency of the bone, ie, the strength of a given amount of bone has been increased. Theoretically, even if the bone strength is weakened, the structural efficiency of the bone may increase. Thus, the std / hom ratio provides additional information that is superior to mere intensity measurements.
また、この情報を用いて、患者毎に最適な治療を選択することができる。例えば、所定の患者について、この比率の値が０．９０であった場合、この患者は、構造的に効率的な骨を有しているので、この患者については、（患者の骨の強度が弱い場合であっても）薬物療法によって、骨の空間的分布を変更する必要がないことがわかる。一方、この比率の値が０．７０の患者の場合、骨の構造的効率性は低い。このような患者の場合、治療によって、骨密度の空間的分布を均質化し、所定の骨量の強度を最適化することが望ましい。実際には、医師は、患者からこの比率の測定値を取得し、その値が１にどれくらい近いかに応じて、（図３Ｂに示すデータ点を用いて、「正常」であるか否かを判定してもよい。）、患者に適切な種類の薬物療法を勧めることができる。また、これらの新たなメトリックを用いることも有益である。 In addition, this information can be used to select an optimal treatment for each patient. For example, for a given patient, if the ratio value is 0.90, the patient has structurally efficient bone, so for this patient (the patient's bone strength is It can be seen that drug therapy does not need to change the spatial distribution of bone (even if it is weak). On the other hand, the structural efficiency of the bone is low for patients with a ratio value of 0.70. For such patients, it is desirable for treatment to homogenize the spatial distribution of bone density and optimize the strength of a given bone mass. In practice, the doctor obtains a measure of this ratio from the patient and determines whether it is “normal” using the data points shown in FIG. 3B, depending on how close the value is to 1. May recommend) the patient to the appropriate type of medication. It is also beneficial to use these new metrics.
図３Ａのグラフは、均質化された強度に対して標準強度をプロットしたグラフである。これは、構造／機能特性の他の側面である。丸で囲まれた点は、主回帰直線から外れているため、他の標本の平均的な振る舞いから外れる振る舞いを示している。図３Ｂは、同じデータについて、標準強度と均質化された強度との比率を、標準強度に対してプロットしたグラフである。これは、構造／機能特性の他の側面である。この比率の値が１に近付くほど、脊椎骨が構造的に効率的な振る舞いをすることが示され、低い値は、密度分布が劣った脊椎骨を表す。図３Ａと同様に丸で囲まれた点は、異常値であることが遙かに明確である。水平線は、均質化された強度比率の平均値及び±２ＳＤ標準偏差を示し、これにより、図３Ａでは明瞭ではなかった他の２つの異常値も特定できる。強調された標本は、強度の平均を僅かに下回っているのみであるが、均質な強度に対する標準値の比率が非常に低く、密度分布は、強度を提供するために効率的でないので、品質が劣っている。 The graph of FIG. 3A is a graph in which the standard intensity is plotted against the homogenized intensity. This is another aspect of structural / functional properties. Since the circled points deviate from the main regression line, the behaviors deviate from the average behavior of other samples. FIG. 3B is a graph plotting the ratio of standard intensity to homogenized intensity versus standard intensity for the same data. This is another aspect of structural / functional properties. The closer this ratio value is to 1, it is shown that the vertebrae behave structurally efficiently, with lower values representing vertebrae with poor density distribution. It is much clearer that the points surrounded by circles as in FIG. 3A are outliers. The horizontal line shows the mean value of the homogenized intensity ratio and ± 2SD standard deviation, which can also identify the other two outliers that were not obvious in FIG. 3A. The highlighted sample is only slightly below the intensity average, but the quality ratio is not good because the ratio of the standard value to the homogeneous intensity is very low and the density distribution is not efficient to provide the intensity. Inferior.
モデルの他の変形例では、モデル内の骨材料に所定の密度の基準値（例えば、１００ｍｇ／ｃｍ３、なお、実際の値は、任意である）を割り当てる。この手法により、脊椎内及び脊椎間の密度の変化を除外し、分析された脊椎骨の密度を同じにすることができる（均質化されたモデルでは、各脊椎骨には、それ自身の平均密度の固有値が割り当てられる）。この結果、脊椎骨内又は脊椎骨間には、材料特性に関する根本的な相違がなくなる。これにより得られるこれらのモデルの「基準（reference）」強度（図１）を骨間で比較した場合、相互の違いは、骨の形状に起因する。このように、この分析では、強度に関する骨形状の独立した効果を分離できる。異なる脊椎骨についての「基準」の強度と、「標準」の強度との比較は、構造／機能特性の定量化の他の具体例である。薬物療法により、ｒｅｆ強度測定値が変化することがあり、この変化は、形状のみの変化に起因する総合的な強度の変化を意味する。このパラメータは、骨の構造／機能特性の他の側面を表している。この測定値は、特に、骨等の複雑な幾何学的構造体では、形状の如何なる側面が強度に最も関連しているかを事前に知ることが困難である場合があるため、有用である。ｒｅｆ強度測定値は、形状の総合的な効果を統合し、したがって、特定の形状的特徴（高さ、幅等）を事前に知る必要はない。 In another variation of the model, a predetermined density reference value (eg, 100 mg / cm 3 , where the actual value is arbitrary) is assigned to the bone material in the model. This technique eliminates density changes within and between the vertebrae and allows the vertebrae analyzed to have the same density (in a homogenized model, each vertebra has its own mean density eigenvalue Is assigned). As a result, there is no fundamental difference in material properties within or between vertebrae. When the resulting “reference” strength (FIG. 1) of these models is compared between bones, the difference between them is due to the shape of the bone. Thus, this analysis can separate the independent effects of bone shape on strength. Comparison of “reference” and “standard” strengths for different vertebrae is another example of quantifying structural / functional properties. Drug therapy can change the ref intensity measurement, which means a change in overall intensity due to a change in shape only. This parameter represents another aspect of the structural / functional properties of the bone. This measurement is particularly useful for complex geometric structures such as bones because it may be difficult to know in advance which aspect of shape is most relevant to strength. The ref intensity measurement integrates the overall effect of the shape and therefore does not require prior knowledge of specific geometric features (height, width, etc.).
実際には、このような測定値は、診断又は治療の評価に役立つ。例えば、腰部の骨粗鬆症及び骨折のリスクの診断において、老化又は病気によって骨の形状がどう変化するかを評価するために、単一の形状的特徴を選択することは非常に困難である。これは、大腿骨近位部の形状が複雑なためである。ｒｅｆ強度測定値により、この制約がなくなる。したがって、他の実施の形態では、ｒｅｆ強度測定値を用いて、骨、又は他の器官又は組織の正味の幾何学的な生体力学的効果を定量化する。異なる時点におけるｒｅｆ強度の比較は、構造／機能特性の他の具体例である。 In practice, such measurements are useful for diagnostic or therapeutic evaluation. For example, in diagnosing lumbar osteoporosis and fracture risk, it is very difficult to select a single geometric feature to evaluate how bone shape changes due to aging or disease. This is because the shape of the proximal part of the femur is complicated. The ref intensity measurement removes this restriction. Thus, in other embodiments, ref intensity measurements are used to quantify the net geometric biomechanical effects of bones or other organs or tissues. Comparison of ref intensity at different times is another example of structural / functional properties.
モデルの他の変形例では、骨に関する負荷条件を変更できる。例えば、あるケースでは、曲げモーメントを骨に適用し、他のケースでは、一定の圧縮負荷と比較した。これら比率（図１のｂｅｎｄ／ａｘｉａｌ）は、一方の負荷モードに対する他方の負荷モードの抵抗の相対的な測定値を提供する。例えば、薬物療法によってｂｅｎｄ／ａｘｉａｌ比率が高くなったということは、治療によって、曲げに類する負荷に対する骨の抵抗が強くなったことを意味する。ここでも、この情報は、単なる骨の強度が向上したことを示す情報に比べて、治療にとってより有益な情報である。また、理論的には、治療（又は老化、病気）によって、骨の強度が低下しても、ｂｅｎｄ／ａｘｉａｌ比率が増加することもある。したがって、この比率は、単に強度の個々の測定値を考慮することに加えて、更に有益な情報を提供する。実際には、圧縮強度が正常であるが、ｂｅｎｄ／ａｘｉａｌの値が低い患者は、骨折のリスクがより高いと考えられる。このように、好ましい実施の形態に基づき、診断及び治療の評価において、この比率を用いることは、非常に有益である。 In other variations of the model, the loading conditions for the bone can be changed. For example, in some cases a bending moment was applied to the bone, and in other cases it was compared to a constant compressive load. These ratios (bend / axial in FIG. 1) provide a relative measurement of the resistance of one load mode to the other load mode. For example, a higher bend / axial ratio due to drug therapy means that the treatment has increased bone resistance to bending-like loads. Again, this information is more useful for treatment than information that simply indicates an improvement in bone strength. Theoretically, even if the bone strength decreases due to treatment (or aging, disease), the bend / axial ratio may increase. Thus, this ratio provides more useful information in addition to simply considering individual measurements of intensity. In practice, patients with normal compressive strength but low values of bend / axial are considered to be at higher risk of fracture. Thus, based on the preferred embodiment, it is highly beneficial to use this ratio in diagnostic and therapeutic evaluation.
他の側面においては、骨材料特性、すなわちＱＣＴ又は他の何らかの医療的なスキャンを生体力学的有限要素モデル又は他の生体力学モデル（分析モデルを含む）に変換する際に用いられた入力特性に骨密度をマッピングするために用いられた関係は、薬物又は他の治療の作用を反映するように変更することができる。また、例えば、骨組織について、これらの関係の変化に関する情報を提供してもよい。このようにして、骨組織の材料特性に対する治療の効果を、骨全体の構造分析に統合でき、及び骨全体の強度に関する効果を定量化できる。また、材料特性について、密度／機械的特性関係が変更された骨と変更されていない骨との応答を比較することによって、特にこの変更に起因する強度の如何なる変化も分離することができる。これにより、椎体の機械的挙動の全体の強度に関して、骨材料の構造機能特性の変更の効果を定量化することができる。 In another aspect, the bone material characteristics, ie, the input characteristics used in converting a QCT or some other medical scan into a biomechanical finite element model or other biomechanical model (including an analytical model). The relationship used to map bone density can be altered to reflect the effects of drugs or other treatments. Further, for example, information regarding changes in these relationships may be provided for bone tissue. In this way, the effect of treatment on the material properties of bone tissue can be integrated into the structural analysis of the entire bone and the effect on the strength of the entire bone can be quantified. Also, with regard to material properties, any changes in strength, particularly due to this change, can be isolated by comparing the response of bones with a modified density / mechanical property relationship to bones that have not been changed. Thereby, the effect of changing the structural functional properties of the bone material can be quantified with respect to the overall strength of the mechanical behavior of the vertebral body.
これらの分析において、モデルの変化に関連する負荷の割合は、生成された元のモデル（図１における標準モデル又は同様なモデル）に対する比較によって、定量化される。他の手法でデータを表現してもよく、例えば、変更されたモデルにおける負荷の大きさの絶対値、又は差分の絶対値対標準モデル又は他の何らかのモデルによってデータを表現してもよい。そして、分析されている骨のデータ点を母集団の応答に対して比較し、検査中の骨が構造／機能挙動の正常範囲内に収まるかを判定する。また、データ点は、同じ患者の以前の測定値に対して、経時的に比較することもできる。また、体内の他の部分の骨に対してデータ点を比較することもでき、例えば、一方が治療によって変更された個人の左右の大腿骨を比較してもよい。コンピュータプログラムを介して、全ての処理を自動化することが好ましい。 In these analyses, the percentage of load associated with model changes is quantified by comparison to the original model generated (standard model or similar model in FIG. 1). The data may be expressed in other ways, for example, the absolute value of the magnitude of the load in the modified model, or the absolute value of the difference versus the standard model or some other model. The bone data points being analyzed are then compared to the population response to determine if the bone under examination falls within the normal range of structural / functional behavior. Data points can also be compared over time to previous measurements of the same patient. Data points can also be compared against bones in other parts of the body, for example, one may compare the left and right femurs of an individual that have been altered by treatment. It is preferable to automate all processing via a computer program.
図４のフローチャートは、好適な実施の形態及び変形例に基づいて、様々な組合せで実行できる処理を示している。 The flowchart of FIG. 4 shows processing that can be executed in various combinations based on the preferred embodiment and modifications.
１．患者の骨又は他の組織をＣＴスキャン又は定量的コンピュータ断層撮影により画像を取得する。 1. The patient's bone or other tissue is imaged by CT scan or quantitative computed tomography.
２．骨又は他の構造の標準化された基準の向きを用いて、及びＣＴスキャンからの形状及び骨密度情報を用いて、骨又は他の構造の有限要素モデルを生成する。 2. A finite element model of the bone or other structure is generated using the standardized reference orientation of the bone or other structure and using the shape and bone density information from the CT scan.
３．モデルを変更し、上述した変数を生成する。所望の生体力学的特性を求める。 3. Change the model to generate the variables described above. Determine the desired biomechanical properties.
４．変形されたモデルの強度の相互の又は標準モデル強度との間の又は何れかの測定値と、形状密度又は形態の何らかの定量的測定値との間の関係を評価し、骨又は他の構造の構造／機能特性を定量化する。 4). Evaluate the relationship between the measured values of the deformed model relative to each other or to the standard model strength and any quantitative measure of shape density or morphology, and evaluate the bone or other structure Quantify structural / functional properties.
５．構造／機能関係の確立された基準と比較する。標準からの偏りは、対象となる特定の負荷条件について、異常な骨を示す。 5. Compare with established criteria for structure / function relationships. Deviation from the standard indicates abnormal bone for the particular load condition of interest.
６．負荷条件を変更して、処理を繰り返す。 6). Change the load condition and repeat the process.
７．経時的変化（例えば、薬物療法又は自然老化の作用）を分析する場合、第２の時点（すなわち、治療後又は一定期間経過後）のＣＴデータについて上述した処理を繰り返し、応答と、確立された基準値とを比較する。 7. When analyzing changes over time (eg, the effects of drug therapy or natural aging), the above process was repeated for the CT data at the second time point (ie, after treatment or after a period of time) to establish a response and Compare with the reference value.
８．部分的な変化（例えば、生体の異なる部分に関する薬物療法の効果）を分析する場合、第２の生体の部分のＣＴデータについて上述した処理を繰り返し、応答を相互に、及び確立された基準に対して比較する。 8). When analyzing partial changes (eg, the effects of drug therapy on different parts of the body), the above process is repeated for the CT data of the second part of the body, and the responses to each other and to established criteria. Compare.
９．薬物療法又は他の治療を分析する場合、適切であれば、骨の組織材料特性に対する薬物の作用を反映するように、入力モデルにおける密度材料特性を仮定した関係を変更する。骨全体に対する薬物の効果を評価するために上述の分析を繰り返す。 9. When analyzing drug therapy or other treatments, if appropriate, the hypothesized relationship of density material properties in the input model is modified to reflect the drug's effect on bone tissue material properties. The above analysis is repeated to evaluate the effect of the drug on the entire bone.
実際の臨床現場では、これらの処理は、好ましくは、コンピュータによって自動化され、すなわち、本発明の実施の形態に基づいてプログラミングされたソフトウェアがコンピュータによって実行される。そして、最終的な出力を含む処理の出力の一部及び／又は全部をデジタル媒体に保存してもよい。 In actual clinical settings, these processes are preferably automated by a computer, i.e. software programmed according to an embodiment of the invention is executed by a computer. Then, a part and / or all of the output of the process including the final output may be stored in a digital medium.
更に、異なる時点にマッピングすることができる骨の特定の点における生体力学的機能の変化を追跡することによって、年齢、病気又は治療の効果の評価に特異性を加えることができる。例えば、ある時点から他の時点への骨の局部的な領域の変化によって、どこの構造がどのように変化しているかをより詳細に区別することができる。第１の時点における対応する有限要素材料特性と共に第２の時点の有限要素の形状を用いれば、治療の幾何学的な作用を分離することができる。同様に、第１の時点における対応する有限要素形状と共に、第２の時点における有限要素の材料特性を用いれば、材料特性の変化及びその分布の作用を詳細に分離することができる。 Furthermore, by tracking changes in biomechanical function at specific points in the bone that can be mapped to different time points, specificity can be added to assessing the effects of age, disease or treatment. For example, a change in a local region of bone from one point in time to another can distinguish in more detail which structure is changing and how. Using the shape of the finite element at the second time point along with the corresponding finite element material properties at the first time point, the treatment geometric effects can be separated. Similarly, if the material properties of the finite element at the second time point are used together with the corresponding finite element shape at the first time point, the change in material properties and the effect of its distribution can be separated in detail.
ここでは、本発明の概念について、骨、特に骨粗鬆症への応用例を用いて説明している。患者毎にパラメータが制御された研究を実行し、構造／機能特性と、母集団標準又は別の時刻の構造／機能特性とを比較する包括的な手法は、患者へのスキャン（例えば、ＱＣＴ、ｐＱＣＴ、ＣＴ、マイクロＣＴ、ＭＲＩ、マイクロＭＲＩ、ＵＳ、ＤＸＡ、ＰＥＴ、エックス線及びこれらの任意の組合せを含む。）から生体力学的コンピュータモデルが生成される全ての分野に応用できる。また、生体力学的分析を用いて、強度又は対象となる他の生体力学的特徴（例えば、硬度、応力、変形、エネルギ吸収、疲労特性、耐久性、破砕強度、亀裂伝搬特性、流体応力、圧力）の測定値を非観血的に生成することができる。このような分析法は、有限要素法に制限されず、分析モデリング、ビーム理論又は複合ビーム理論、又は他の形式のビーム理論、破壊力学、複合材料分析、又は固体及び流体力学を含む連続体力学の一種、熱輸送及び質量輸送分析、動的分析又は他の機械的分析の一種、又は有限差分法を含む他の数値分析の一種等も含まれる。例えば、骨移植システムの分析では、骨／移植界面条件を変更し、又は骨及び／又は移植の材料特性を変更してもよい。血流分析のモデルでは、血管の弾性特性を変更してもよく、血流条件を変更してもよく、血管の形状を変更してもよく、血管間の合流点の形状を変更してもよく、又はプラーク又は他の遮蔽物の形状を変更してもよく、血流に障害を加えてもよい。何れの場合も、適切な医療画像及び工学理論が用いられる。 Here, the concept of the present invention is described using an application example to bone, particularly osteoporosis. A comprehensive approach to perform parameter-controlled studies for each patient and compare structural / functional properties to population standard or other time-of-day structural / functional properties is a scan to the patient (eg, QCT, including pQCT, CT, micro CT, MRI, micro MRI, US, DXA, PET, X-rays, and any combination thereof) and can be applied to all fields where biomechanical computer models are generated. Also, using biomechanical analysis, strength or other biomechanical features of interest (eg hardness, stress, deformation, energy absorption, fatigue properties, durability, crush strength, crack propagation properties, fluid stress, pressure ) Measurement values can be generated non-invasively. Such analytical methods are not limited to finite element methods, but include analytical modeling, beam theory or composite beam theory, or other forms of beam theory, fracture mechanics, composite material analysis, or continuum mechanics including solid and fluid mechanics. A kind of heat transport and mass transport analysis, a kind of dynamic analysis or other mechanical analysis, or a kind of other numerical analysis including a finite difference method. For example, in an analysis of a bone graft system, the bone / graft interface conditions may be changed, or the material properties of the bone and / or graft may be altered. In the blood flow analysis model, the elasticity characteristics of the blood vessels may be changed, the blood flow conditions may be changed, the shape of the blood vessels may be changed, and the shape of the confluence between the blood vessels may be changed. Or, the shape of the plaque or other shield may be altered and the blood flow may be impaired. In either case, appropriate medical images and engineering theories are used.
また、本発明は、将来の可能な条件のリスク又は適合性を評価するために、将来の仮説的な変化のシミュレートにも適用できる。これは、手術の計画、薬物療法の決定及び予防療法の考慮に有用である。例えば、椎体形成術及びこれに類する医療的施術では、骨折を治療し、骨を補強するために、骨セメントを導入する。この場合、導入される骨セメントの量又は位置が異なるように患者の骨の画像ベースのモデルを変更することによって、想定された負荷条件の下で、これにより得られる骨全体の強度又は変形の変化を測定し、相互に又は導入される骨セメントの相対的な量に対してこの特徴をプロットすることによって、骨移植システムの構造／機能特性を定量化できる。この情報を用いて、患者にとって最適な治療を判定でき、すなわち、導入する骨セメント材料の量と、骨内に骨セメントを導入する位置とを適切に判断することができる。 The present invention can also be applied to simulate future hypothetical changes in order to assess the risk or suitability of possible future conditions. This is useful for surgical planning, pharmacotherapy decisions and prophylactic treatment considerations. For example, vertebroplasty and similar medical procedures introduce bone cement to treat fractures and reinforce bones. In this case, by changing the image-based model of the patient's bone so that the amount or location of the bone cement introduced is different, under the assumed loading conditions, the resulting strength or deformation of the entire bone By measuring the change and plotting this feature against each other or the relative amount of bone cement introduced, the structural / functional properties of the bone grafting system can be quantified. This information can be used to determine the optimal treatment for the patient, i.e., the amount of bone cement material to be introduced and the location to introduce bone cement into the bone.
他の骨移植の具体例として、移植前の患者の骨の画像のモデルを、移植を含むように変更できる。そして、移植前後の生体力学的パラメータ、例えば、骨の応力等を比較することができる。また、移植されるプロテーゼの材料特性、形状、位置又は寸法を変更できる。そして、例えば、所定の負荷条件の組に対する骨の応力の低下を最小化する骨等、所定の基準に基づいて、移植を選択できる。これにより、外科医は、特定の患者について適切な移植構成を選択することができる。 As another bone grafting example, the model of the bone image of the patient prior to grafting can be modified to include the grafting. Then, biomechanical parameters before and after transplantation, such as bone stress, can be compared. Also, the material properties, shape, position or dimensions of the prosthesis to be implanted can be changed. Then, for example, a transplant can be selected based on a predetermined criterion, such as a bone that minimizes bone stress reduction for a set of predetermined load conditions. This allows the surgeon to select an appropriate implantation configuration for a particular patient.
他の特定の具体例として、骨のマイクロＣＴスキャン又は骨のマイクロＭＲＩスキャンから生成されたモデルの分析を用いる骨梁の分析に本発明を適用してもよい。ニューイット（Newitt）他（骨粗鬆症国際学会１３：２７８−２８７，２００２）には、患者の橈骨遠位端（手首）の高解像度ＭＲＩスキャンに由来する有限要素モデルを用いて、骨梁の異方性弾性特性を非観血的に推定する手法が開示されている。このような分析技術、若しくは画像から有限要素モデル又は他の種類の構造モデルを生成する他の技術を用いて、このようなモデルを管理された手法で変更し、構造／機能特性を定量化することによって治療の作用のメカニズムを判定することができる。ある変形例では、ニューイットの技術を用いて、患者の手首のマイクロＭＲＩスキャンに基づいて、標準モデルを生成し、手首の骨梁の小さな標本を分離し、標準強度（又は同等な尺度）の分析を実行する。このモデルは、外側の骨組織の構造的な役割を判定するために、個々の骨梁組織の外側のボクセル（outer voxel）を取り除くことによって変形される。元のモデルの強度と、骨が取り除かれたモデルの強度との間の関係に基づき、骨の構造／機能特性が定量化される。他の変形例として、１つ以上の個々の骨梁組織を取り除き、又は追加してもよい。他の変形例として、特定の種類の薬物療法の可能な効果をシミュレートするために、骨の層を追加してもよく、又は様々な療法の作用に関する他の研究から既知である所定の基準に基づいて、骨を追加してもよい。他の変形例として、骨梁組織の表面に沿って小さい空洞を設け、総合的な機械的特性に対する骨吸収空間の効果を調べてもよい。他の変形例として、このような骨吸収空間に所定の材料を埋め込んだ場合の機械的特性を算出してもよい。なお、骨吸収空間の生成、個々の骨梁組織の削除、又は既存の骨吸収空間の埋め込みは、何れも、病気及び治療が骨梁の機械的特性に影響する可能性があるメカニズムである。他の変形例として、所定数の骨梁組織がなくなるまで、骨材料を取り除き、又はある所定の基準に基づいて、ある質量の骨を取り除いてもよい。何れの場合も、変更の前後の骨の機械的特性を算出し、構造／機能特性を定量化する。上述したような定量化は、診断ツールとして、又は治療の評価、若しくは新たな又は理解が不十分な治療の作用の検査のために実行してもよい。 As another specific embodiment, the present invention may be applied to trabecular analysis using analysis of models generated from bone micro CT scans or bone micro MRI scans. Newitt et al. (International Society for Osteoporosis 13: 278-287, 2002) uses a finite element model derived from a high-resolution MRI scan of the patient's distal radius (wrist) to analyze the trabecular anisotropy. A technique for non-invasively estimating the elastic properties is disclosed. Using such analytical techniques, or other techniques that generate finite element models or other types of structural models from images, modify such models in a controlled manner and quantify structural / functional properties Thus, the mechanism of action of treatment can be determined. In one variation, a Newit technique is used to generate a standard model based on a micro-MRI scan of the patient's wrist, isolate a small specimen of the wrist trabeculae, and a standard intensity (or equivalent scale). Perform analysis. This model is modified by removing the outer voxel of the individual trabecular tissue to determine the structural role of the outer bone tissue. Based on the relationship between the strength of the original model and the strength of the model from which the bone has been removed, the structural / functional properties of the bone are quantified. As another variation, one or more individual trabecular tissues may be removed or added. As another variation, bone layers may be added to simulate the possible effects of certain types of medications, or predetermined criteria known from other studies on the action of various therapies Based on the above, bones may be added. As another variation, a small cavity may be provided along the surface of the trabecular tissue to investigate the effect of the bone resorption space on the overall mechanical properties. As another modification, mechanical characteristics when a predetermined material is embedded in such a bone resorption space may be calculated. It should be noted that the creation of bone resorption spaces, the deletion of individual trabecular tissue, or the embedding of existing bone resorption spaces are all mechanisms by which disease and treatment can affect the mechanical properties of trabecular bone. As another variation, bone material may be removed or a mass of bone may be removed based on some predetermined criteria until there is no predetermined number of trabecular tissues. In either case, the mechanical properties of the bone before and after the change are calculated and the structural / functional properties are quantified. Quantification as described above may be performed as a diagnostic tool or for the assessment of therapy or for testing the effects of new or poorly understood therapy.
上述したような、パラメータを制御した研究によって、特定の患者毎に適切な治療を判定する技術が提供される。例えば、患者の骨は、個々の骨梁組織の支柱が損失しやすい場合があり、これを知ることにより、その患者にとって、ある治療が他の治療より適切であると判断できることがある。また、他の患者は、骨吸収阻害又は同様の治療に関連する小さな空洞の埋め込みの効果が低い場合もあり、このために、異なる種類の治療を勧めることもできる。 A parameter-controlled study, as described above, provides a technique for determining the appropriate treatment for a particular patient. For example, a patient's bone may be prone to loss of individual trabecular tissue struts, and knowing this may determine that a treatment is more appropriate for the patient than others. Other patients may also be less effective at implanting small cavities associated with bone resorption inhibition or similar treatments, which may encourage different types of treatment.
また、マイクロＣＴ及び有限要素モデリングを臨床的に用いて、及び新世代のマイクロＣＴ及び周辺ＱＣＴスキャナ（例えば、Ｓｃａｎｃｏ社（Scanco, Inc）製のエクストリームＣＴ（Extreme CT））を臨床的に用いて骨梁の強度と弾性特性を評価することができる。このような分析技術を用いて、患者毎に固有の骨梁マイクロ構造のモデルを生成することができる（図５参照）。そして、これらの画像を有限要素又は他のモデルに変換し、総合的な骨試料の強度又は硬度を算出することができる。通常、このような処理は、橈骨遠位端において行われるが、例えば、踵骨及び脛骨等の部位についても可能である。ストーバー（Stauber）他は、このような応用例の１つを開示している（"A finite element beam-model for efficient simulation of large-scale porous structures." Comput Methods Biomech Biomed Engin.7:9-16, 2004）。この技術は、非観血的に行うことができ、又は患者から得られた生体組織検査に適用でき、又は基礎科学研究において死体に適用してもよく、前臨床試験研究において動物に適用してもよい。 Also, clinically using micro CT and finite element modeling, and clinically using a new generation of micro CT and peripheral QCT scanners (eg, Extreme CT from Scanco, Inc.). The trabecular strength and elastic properties can be evaluated. Using such an analysis technique, a model of trabecular microstructure unique to each patient can be generated (see FIG. 5). These images can then be converted into finite elements or other models to calculate the overall bone sample strength or hardness. Usually, such processing is performed at the distal end of the radius, but it is also possible for sites such as the radius and the tibia. Stauber et al. Disclose one such application ("A finite element beam-model for efficient simulation of large-scale porous structures." Comput Methods Biomech Biomed Engin. 7: 9-16 , 2004). This technique can be performed non-invasively, or can be applied to biopsy obtained from patients, or applied to cadaver in basic science research and applied to animals in preclinical studies. Also good.
図６は、骨密度の変化を示す３ＤマイクロＣＴベースの有限要素モデルの典型的な断面を示している。このようなモデルに対して有限要素法分析を実行し、強度又は硬度又は他の機械的特性を非観血的に推定することができる。他の実施の形態では、パラメータが制御された研究を実行し、構造／機能特性を定量化してもよい。 FIG. 6 shows a typical cross section of a 3D micro CT based finite element model showing changes in bone density. Finite element analysis can be performed on such models to estimate non-invasively strength or hardness or other mechanical properties. In other embodiments, parameter-controlled studies may be performed to quantify structural / functional properties.
他の実施の形態では、マイクロＣＴスキャンから標準モデルを生成した後に、このモデルについて、一組のパラメータが制御された研究を実行する。上述のように、マイクロＭＲＩモデルについても、同じような変形を作成することができる。ＣＴスキャンは、ＭＲＩスキャンと異なり、骨組織の石灰化の度合いに関する定量的データを提供できる。したがって、他の制御された変形を行うことができ、例えば、骨全体の臨床的ＱＣＴモデルに関連して上述したｈｏｍ分析と同様に、各ボクセルにおける骨密度を平均化し、幾つか又は全ての要素において平均化された骨密度を用いることができる。石灰化の空間的変化のあるモデルとないモデルとを比較することによって、石灰化の空間的分布の変化に関する機械的応答を抽出することができる。更なる変形例として、椎体のｒｅｆ強度測定値について上述した手法と同様の手法によって、患者の骨梁に固定された値の骨密度を割り当て、マイクロ構造の形状の効果を分離することもできる。 In another embodiment, after generating a standard model from a micro CT scan, a set of parameter controlled studies is performed on the model. As described above, the same deformation can be created for the micro MRI model. CT scans, unlike MRI scans, can provide quantitative data regarding the degree of bone tissue mineralization. Thus, other controlled deformations can be made, eg, averaging bone density in each voxel, similar to the hom analysis described above in connection with the whole bone clinical QCT model, and some or all factors The bone density averaged in can be used. By comparing the model with and without the spatial change of calcification, the mechanical response regarding the change of the spatial distribution of calcification can be extracted. As a further modification, the bone density of the value fixed to the patient's trabeculae can be assigned and the effect of the shape of the microstructure can be separated by a technique similar to that described above for the vertebral body ref strength measurement. .
図６は、マイクロＣＴを用いて約７０ミクロンの空間分解能によって撮像された骨梁の５ｍｍ立方の標本の典型的な断面を示している。白い領域は、石灰化が進んだより密度が高い骨を表し、グレーの領域は、石灰化の程度が低い、低密度な領域を表している。黒い領域は、骨髄を表している。標本内における石灰化の空間的変化の独立した効果を評価するために、標本を均質化でき、すなわち、幾つか又は全ての要素に、同じ値の密度を割り当て、均質化された標本の強度と、均質化されていない同じ標本の強度とを比較することができる。 FIG. 6 shows a typical cross section of a 5 mm cubic specimen of trabecular bone imaged with a spatial resolution of about 70 microns using micro CT. White areas represent bones with higher density than advanced calcification, and gray areas represent low density areas with a low degree of calcification. The black area represents the bone marrow. In order to assess the independent effects of spatial variation of calcification within the specimen, the specimen can be homogenized, i.e. several or all elements are assigned the same value of density, and the homogenized specimen strength and Compare the intensity of the same specimen, not homogenized.
本明細書に開示した技術によって、罹患している又は構造的に異常な骨を、より包括的に、より敏感に及び／又はより明確に特定できる。更に、本発明は、骨全体の強度に影響する老化、病気、様々な治療及び時間の経過を評価する処理を提供する。本発明の実施の形態は、骨全体だけではなく、生体力学的振る舞いが臨床的に注目される骨組織及び他の如何なる組織又は器官にも適用できる。 With the techniques disclosed herein, affected or structurally abnormal bones can be identified more comprehensively, more sensitively and / or more clearly. Furthermore, the present invention provides a process for assessing aging, illness, various treatments and time courses that affect the strength of the entire bone. Embodiments of the present invention can be applied not only to the entire bone, but also to bone tissue and any other tissue or organ where biomechanical behavior is clinically noted.
心臓血管に本発明を適用すれば、ＭＲＩ画像又はＣＴ画像及びこれらの組合せを含む様々な形式の画像から患者の血管のモデルを生成できる。このようなモデルによって、患者固有の血管の形状を検査することができる。また、プラークの存在を明らかにすることもできる。そして、有限要素モデリング又は他の演算モデリングを用いて、せん断応力を含む血管内の応力を算出できる。仮定的な又は患者固有の血流状態を推定することもできる。 When the present invention is applied to the cardiovascular system, a model of a patient's blood vessel can be generated from various types of images including MRI images or CT images and combinations thereof. With such a model, the patient-specific blood vessel shape can be examined. It can also reveal the presence of plaque. Then, finite element modeling or other computational modeling can be used to calculate the stress in the blood vessel, including shear stress. Hypothetical or patient-specific blood flow conditions can also be estimated.
このようなモデルに対する制御可能なパラメータの組を用いた研究によって、プラークの存在を完全又は部分的に取り除くことができ、これにより得られるストレス又は他の結果パラメータの変化を評価することができる。患者固有のモデルの他の変形例として、血管の材料特性を変更してもよく、プラークの空間的分布又は他の材料分布を部分的又は完全に追加又は除外することができる。他の変形例として、何らかの手法で血管の形状を追加、削減又は変更することもできる。同じ所定の形状特性について調べ、その幾何学的性質を用いた演算に対して応答を比較してもよい。他の変形例として、血流状態を変更してもよい。例えば、患者固有の血流の推定値を用いることに代えて、共通の「公称（nominal）」血流条件を患者に適用し、患者固有の血流と、「公称」血流条件との間の応答の差分を定量化してもよい。これに代えて、異なる種類の「公称」血流条件を患者に適用し、異なる負荷条件に対する応答の違いを定量化してもよい。これらは、病気の診断、又は老化、病気又は治療の評価のための追加的情報を提供するために患者固有のモデルを変更する一例に過ぎない。他にも多くの応用例及びモデルが利用できることは、当業者にとって明らかである。 A study with a controllable set of parameters for such a model can completely or partially eliminate the presence of plaque and thereby assess the resulting stress or other resulting parameter changes. As another variation of the patient-specific model, the material properties of the blood vessels may be altered, and the plaque spatial distribution or other material distribution may be partially or completely added or excluded. As another modification, the shape of the blood vessel can be added, reduced, or changed by some method. The same predetermined shape characteristic may be examined and the response compared to an operation using that geometric property. As another modification, the blood flow state may be changed. For example, instead of using a patient-specific blood flow estimate, a common “nominal” blood flow condition is applied to the patient and the patient-specific blood flow The difference in response may be quantified. Alternatively, different types of “nominal” blood flow conditions may be applied to the patient to quantify the difference in response to different loading conditions. These are only examples of modifying patient-specific models to provide additional information for disease diagnosis or aging, disease or treatment assessment. It will be apparent to those skilled in the art that many other applications and models are available.
本発明は、当業者に明らかなように、ここに示した本発明の説明に基づいてプログラミングされた、従来の汎用用途の又は特定用途のデジタルコンピュータ又はマイクロプロセッサを用いて好適に実現することができる。 As will be apparent to those skilled in the art, the present invention may be suitably implemented using a conventional general purpose or special purpose digital computer or microprocessor programmed according to the description of the invention presented herein. it can.
また、ソフトウェア分野の知識を有する者には明らかであるが、ここに示した本発明の説明に基づいて、プログラマは適切なソフトウェアコーディングを容易に行うことができる。更に、当業者に明らかなように、本発明は、特定用途向け集積回路、又は従来の電気部品回路の適切なネットワークを相互接続することによっても実現できる。 Further, as will be apparent to those having knowledge in the software field, based on the description of the present invention shown here, a programmer can easily perform appropriate software coding. Further, as will be apparent to those skilled in the art, the present invention can also be realized by interconnecting a suitable network of application specific integrated circuits or conventional electrical component circuits.
本発明は、コンピュータを制御して本発明の処理を行わせるために使用される命令を記録した記録媒体（媒体）であるコンピュータプログラム製品を含む。記録媒体には、以下に限定するものではないが、フロッピィディスク（商標）、ミニディスク（ＭＤ、商標）、光ディスク、ＤＶＤ、ＣＤ−ＲＯＭ、マイクロドライブ、光磁気ディスク、ＲＯＭ、ＲＡＭ、ＥＰＲＯＭ、ＤＲＡＭ、ＶＲＡＭ、フラッシュメモリ（フラッシュカードを含む）、磁気又は光カード、ナノシステム（分子メモリＩＣを含む）、ＲＡＩＤデバイス、遠隔データ記憶装置／アーカイブ／ウェアハウス、或いは命令及び／又はデータを記憶する任意の種類の適当な媒体又は装置が含まれる。 The present invention includes a computer program product that is a recording medium (medium) that records instructions used to control a computer to perform the processing of the present invention. The recording medium is not limited to the following, but is a floppy disk (trademark), mini disk (MD, trademark), optical disk, DVD, CD-ROM, microdrive, magneto-optical disk, ROM, RAM, EPROM, DRAM. , VRAM, flash memory (including flash card), magnetic or optical card, nano system (including molecular memory IC), RAID device, remote data storage / archive / warehouse, or any instruction and / or data storage Any type of suitable medium or device is included.
本発明は、コンピュータにより読み取り可能な、任意の媒体に記憶され、汎用／専用コンピュータのハードウェア又はマイクロプロセッサを制御し、コンピュータ又はマイクロプロセッサを人であるユーザとインタラクトさせ、或いは他のメカニズムにより本発明の成果を利用する、ソフトウェアを含む。この種のソフトウェアは、限定する意味でなく、デバイスドライバ、オペレーティングシステム、及びユーザアプリケーションを含む。更に、この種のコンピュータにより読み取り可能な媒体は上述したように、本発明を実施するソフトウェアを含む。 The present invention is stored in any medium readable by a computer and controls the hardware or microprocessor of a general purpose / dedicated computer, allowing the computer or microprocessor to interact with a human user, or by other mechanisms. Includes software that utilizes the results of the invention. Such software includes, but is not limited to, device drivers, operating systems, and user applications. Further, this type of computer readable medium includes software implementing the present invention, as described above.
本発明の範囲及び思想から逸脱することなく、上述の好ましい具体例を様々な分野に応用し、又は変更できることは当業者にとって明らかである。したがって、本発明は、添付の請求の範囲に基づき、上述した具体例とは異なる形態でも実現できることは明らかである。本発明は、本明細書に開示した実施の形態に限定されず、添付の特許請求の範囲に定義した本発明の範囲及びその構造的及び機能的な等価物から逸脱することなく、これらの実施の形態を変形又は変更できる。 It will be apparent to those skilled in the art that the preferred embodiments described above can be applied or modified in various fields without departing from the scope and spirit of the invention. Therefore, it is apparent that the present invention can be realized in a form different from the specific examples described above based on the appended claims. The present invention is not limited to the embodiments disclosed herein, and these implementations may be made without departing from the scope of the invention and its structural and functional equivalents as defined in the appended claims. The form can be modified or changed.
以上では、好ましい実施の形態及び／又は特許請求の範囲の定義に基づいて実行される処理について、選択された活字のシーケンスを用いて記述した。なお、これらのシーケンスは、印刷の都合のために選択されたものであり、本発明の処理を実行する順序を特に指定することを意図しない。 In the above, the processing executed based on the preferred embodiment and / or the definition of the scope of claims has been described using the sequence of the selected type. These sequences are selected for the convenience of printing, and are not intended to specify the order in which the processing of the present invention is executed.
更に、本明細書の背景技術、課題を解決するための手段及び発明を実施するための最良の形態に引用した全ての文献は、好ましい実施の形態の変形例及び要素を開示するものとして、全て本願に引用される。 Further, all documents cited in the background art, means for solving the problems, and the best mode for carrying out the invention of this specification are all intended to disclose modifications and elements of the preferred embodiment. Cited in this application.
Claims (45)
上記生体内の構造の画像を受信するステップと、
上記画像に基づいて、該構造の構造モデルを生成するステップと、
上記構造モデルに基づいて、第１の生体力学的量を算出するステップと、
上記構造モデルを変更して、変形モデルを生成するステップと、
上記変形モデルに基づいて、第２の生体力学的量を算出するステップと、
上記第１及び第２の生体力学的量を比較するステップと、
上記比較の結果をデジタルメディアに保存するステップとを有する特性判定方法。 In a characteristic determination method for determining one or more structural / functional characteristics of a human or animal biological structure from an image of the human or animal biological structure,
Receiving an image of the in-vivo structure;
Generating a structural model of the structure based on the image;
Calculating a first biomechanical quantity based on the structural model;
Modifying the structural model to generate a deformation model;
Calculating a second biomechanical quantity based on the deformation model;
Comparing the first and second biomechanical quantities;
Storing the result of the comparison in a digital medium.
上記第２の画像に基づいて、上記構造の第２の構造モデルを生成するステップと、
上記第２の構造モデルに基づいて、第３の生体力学的量を算出するステップと、
上記第２の構造モデルを変更して第２の変形モデルを生成するステップと、
上記第２の変形モデルに基づいて、第４の生体力学的量を算出するステップとを有し、上記比較は、上記第３及び第４の生体力学的量の比較を含むことを特徴とする請求項１記載の特性判定方法。 Receiving a second image of the structure at a time different from the time at which the first image was acquired;
Generating a second structural model of the structure based on the second image;
Calculating a third biomechanical quantity based on the second structural model;
Modifying the second structural model to generate a second deformation model;
Calculating a fourth biomechanical quantity based on the second deformation model, wherein the comparison includes a comparison of the third and fourth biomechanical quantities. The characteristic determination method according to claim 1.
上記第２の画像に基づいて、第２の構造モデルを生成するステップと、
上記第２の構造モデルに基づいて、第３の生体力学的量を算出するステップと、
上記第２の構造モデルを変更して第２の変形モデルを生成するステップと、
上記第２の変形モデルに基づいて、第４の生体力学的量を算出するステップとを有し、上記比較は、第３及び第４の生体力学的量の比較を含むことを特徴とするを有する請求項１記載の特性判定方法。 Receiving a second image of a portion of the living body different from the structure;
Generating a second structural model based on the second image;
Calculating a third biomechanical quantity based on the second structural model;
Modifying the second structural model to generate a second deformation model;
Calculating a fourth biomechanical quantity based on the second deformation model, wherein the comparison includes a comparison of the third and fourth biomechanical quantities. The characteristic determination method according to claim 1.
生体内の構造の画像を受信するステップと、
上記画像に基づいて、該構造の構造モデルを生成するステップと、
上記構造モデルに基づいて、第１の生体力学的量を算出するステップと、
上記構造モデルを変更して、変形モデルを生成するステップと、
上記変形モデルに基づいて、第２の生体力学的量を算出するステップと、
上記第１及び第２の生体力学的量を比較するステップと、
上記比較の結果をデジタルメディアに保存するステップとを有する記録装置。 In a recording device readable by one or more processors having a code readable by a processor, the code readable by the processor is obtained from an image of a human or animal biological structure. Or programming the one or more processors to perform a characterization method for determining a plurality of structural / functional properties, the characterization method comprising:
Receiving an image of an in-vivo structure;
Generating a structural model of the structure based on the image;
Calculating a first biomechanical quantity based on the structural model;
Modifying the structural model to generate a deformation model;
Calculating a second biomechanical quantity based on the deformation model;
Comparing the first and second biomechanical quantities;
And a step of storing the result of the comparison in a digital medium.
第１の期間において、上記生体の構造の画像を受信するステップと、
上記画像に基づいて、該構造の構造モデルを生成するステップと、
上記構造モデルに基づいて、第１の生体力学的量を算出するステップと、
上記構造モデルを変更して、変形モデルを生成するステップと、
上記変形モデルに基づいて、第２の生体力学的量を算出するステップと、
上記第１及び第２の生体力学的量を比較するステップと、
上記第１及び第２の生体力学的量の比較に基づいて、１又は複数の構造／機能特性を判定するステップと、
第２の期間に取得された上記構造の第２の画像を受信するステップと、
上記第２の画像に基づいて、上記構造の第２の構造モデルを生成するステップと、
上記第２の構造モデルに基づいて、第３の生体力学的量を算出するステップと、
上記第２の構造モデルを変更して第２の変形モデルを生成するステップと、
上記第２の変形モデルに基づいて、第４の生体力学的量を算出するステップと、
上記第３及び第４の生体力学的量を比較するステップと、
第１及び第２の生体力学的量の比較の結果と、第３及び第４の生体力学的量の比較の結果とを比較し、上記構造の老化、病気又は治療に関する効果を評価及び／又は診断するステップとを有する評価方法。 In an evaluation method for evaluating the effects of aging, disease or treatment on the structure of a human or animal body,
Receiving an image of the structure of the living body in a first period;
Generating a structural model of the structure based on the image;
Calculating a first biomechanical quantity based on the structural model;
Modifying the structural model to generate a deformation model;
Calculating a second biomechanical quantity based on the deformation model;
Comparing the first and second biomechanical quantities;
Determining one or more structural / functional properties based on the comparison of the first and second biomechanical quantities;
Receiving a second image of the structure acquired in a second period;
Generating a second structural model of the structure based on the second image;
Calculating a third biomechanical quantity based on the second structural model;
Modifying the second structural model to generate a second deformation model;
Calculating a fourth biomechanical quantity based on the second deformation model;
Comparing the third and fourth biomechanical quantities;
Comparing the results of the comparison of the first and second biomechanical quantities with the results of comparing the third and fourth biomechanical quantities to evaluate the effects of the structure on aging, disease or treatment and / or An evaluation method comprising the step of diagnosing.
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US61460504P | 2004-09-30 | 2004-09-30 | |
| PCT/US2005/034892 WO2006039358A2 (en) | 2004-09-30 | 2005-09-30 | Method for assessment of the structure-function characteristics of structures in a human or animal body |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JP2008514368A true JP2008514368A (en) | 2008-05-08 |
Family
ID=38066775
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2007534740A Withdrawn JP2008514368A (en) | 2004-09-30 | 2005-09-30 | Characteristic determination method, recording apparatus, evaluation method |
Country Status (3)
| Country | Link |
|---|---|
| EP (1) | EP1805687A2 (en) |
| JP (1) | JP2008514368A (en) |
| CA (1) | CA2581370A1 (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007260001A (en) * | 2006-03-28 | 2007-10-11 | Univ Of Tsukuba | Magnetic resonance imaging apparatus and imaging method |
| JP2010214098A (en) * | 2009-02-23 | 2010-09-30 | Kochi Univ Of Technology | Measuring instrument and measuring method of load on intervertebral disk |
| WO2014077233A1 (en) * | 2012-11-19 | 2014-05-22 | 株式会社 東芝 | Blood vessel analysis device, medical image diagnostic device, and blood vessel analysis method |
| JP2016096889A (en) * | 2014-11-19 | 2016-05-30 | 株式会社東芝 | Image analysis apparatus, image analysis method and program |
| JP2017029786A (en) * | 2016-10-13 | 2017-02-09 | 株式会社東芝 | Tubular structure analyzer, tubular structure analysis method, and tubular structure analysis program |
| JP2018030016A (en) * | 2012-04-19 | 2018-03-01 | オブシェストボ エス オグラノチェノイ オトヴェツトヴェノスチウ “パラファーム” | Apparatus for diagnosing osteoporosis by method for checking processes of closing cavity lesions for assessing effectiveness of various osteoprotectors |
| JP2019136542A (en) * | 2019-04-25 | 2019-08-22 | 株式会社東芝 | Image analysis apparatus, image analysis method and program |
-
2005
- 2005-09-30 CA CA002581370A patent/CA2581370A1/en not_active Abandoned
- 2005-09-30 JP JP2007534740A patent/JP2008514368A/en not_active Withdrawn
- 2005-09-30 EP EP05815202A patent/EP1805687A2/en not_active Withdrawn
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007260001A (en) * | 2006-03-28 | 2007-10-11 | Univ Of Tsukuba | Magnetic resonance imaging apparatus and imaging method |
| JP4576534B2 (en) * | 2006-03-28 | 2010-11-10 | 国立大学法人 筑波大学 | Magnetic resonance imaging apparatus and imaging method |
| JP2010214098A (en) * | 2009-02-23 | 2010-09-30 | Kochi Univ Of Technology | Measuring instrument and measuring method of load on intervertebral disk |
| JP2018030016A (en) * | 2012-04-19 | 2018-03-01 | オブシェストボ エス オグラノチェノイ オトヴェツトヴェノスチウ “パラファーム” | Apparatus for diagnosing osteoporosis by method for checking processes of closing cavity lesions for assessing effectiveness of various osteoprotectors |
| WO2014077233A1 (en) * | 2012-11-19 | 2014-05-22 | 株式会社 東芝 | Blood vessel analysis device, medical image diagnostic device, and blood vessel analysis method |
| JP2014100249A (en) * | 2012-11-19 | 2014-06-05 | Toshiba Corp | Blood vessel analysis device, medical image diagnostic device, blood vessel analysis method, and blood vessel analysis program |
| CN107582041A (en) * | 2012-11-19 | 2018-01-16 | 株式会社东芝 | Tubular structure resolver and tubular structure analytic method |
| CN107582041B (en) * | 2012-11-19 | 2021-07-06 | 株式会社东芝 | Blood vessel analysis device and blood vessel analysis method |
| US11357409B2 (en) | 2012-11-19 | 2022-06-14 | Kabushiki Kaisha Toshiba | Blood vessel analysis apparatus, medical image diagnosis apparatus, and blood vessel analysis method |
| JP2016096889A (en) * | 2014-11-19 | 2016-05-30 | 株式会社東芝 | Image analysis apparatus, image analysis method and program |
| JP2017029786A (en) * | 2016-10-13 | 2017-02-09 | 株式会社東芝 | Tubular structure analyzer, tubular structure analysis method, and tubular structure analysis program |
| JP2019136542A (en) * | 2019-04-25 | 2019-08-22 | 株式会社東芝 | Image analysis apparatus, image analysis method and program |
Also Published As
| Publication number | Publication date |
|---|---|
| CA2581370A1 (en) | 2006-04-13 |
| EP1805687A2 (en) | 2007-07-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20060069318A1 (en) | Method for assessment of the structure-function characteristics of structures in a human or animal body | |
| Shefelbine et al. | Prediction of fracture callus mechanical properties using micro-CT images and voxel-based finite element analysis | |
| Eggermont et al. | Calibration with or without phantom for fracture risk prediction in cancer patients with femoral bone metastases using CT-based finite element models | |
| McErlain et al. | Subchondral cysts create increased intra-osseous stress in early knee OA: a finite element analysis using simulated lesions | |
| Benca et al. | The insufficiencies of risk analysis of impending pathological fractures in patients with femoral metastases: A literature review | |
| JP2008514368A (en) | Characteristic determination method, recording apparatus, evaluation method | |
| Mirzaei et al. | On prediction of the strength levels and failure patterns of human vertebrae using quantitative computed tomography (QCT)-based finite element method | |
| Podshivalov et al. | On the road to personalized medicine: multiscale computational modeling of bone tissue | |
| Groenen et al. | Case‐specific non‐linear finite element models to predict failure behavior in two functional spinal units | |
| Corona-Castuera et al. | Design and fabrication of a customized partial hip prosthesis employing CT-scan data and lattice porous structures | |
| Nazemi et al. | Optimizing finite element predictions of local subchondral bone structural stiffness using neural network-derived density-modulus relationships for proximal tibial subchondral cortical and trabecular bone | |
| Falcinelli et al. | Mechanical behavior of metastatic femurs through patient-specific computational models accounting for bone-metastasis interaction | |
| Suzuki et al. | Biomechanics of callus in the bone healing process, determined by specimen-specific finite element analysis | |
| Wan et al. | On fatigue failure prediction of prosthetic devices through XFEM analysis | |
| Tho | Bone and joints modelling with individualized geometric and mechanical properties derived from medical images | |
| Hauser et al. | Anatomic variation of structural properties of periacetabular bone as a function of age: a quantitative computed tomography study | |
| Bhattacharya et al. | Experimental quantification of bone mechanics | |
| Petrović et al. | Imaging in clinical and preclinical practice | |
| Delpuech et al. | Failure prediction of tumoral bone with osteolytic lesion in mice | |
| Chen | Verification and validation of microCT-based finite element models of bone tissue biomechanics | |
| Groenen et al. | Case‐specific non‐linear finite element models to predict failure behavior in 3‐segment spinal units | |
| Baer | Quantitative spectral photon-counting computed Tomography imaging of bone and cartilage degradation in osteoarthritis and repair strategies | |
| Smit et al. | In-silico medical device testing of anatomically and mechanically conforming patient-specific spinal fusion cages designed by full-scale topology optimization | |
| Fairweather | Custom Tri-Flange Acetabular Cups | |
| Tawara et al. | Mechanical therapeutic effects in osteoporotic L1-vertebrae evaluated by nonlinear patient-specific finite element analysis |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A300 | Application deemed to be withdrawn because no request for examination was validly filed |
Free format text: JAPANESE INTERMEDIATE CODE: A300 Effective date: 20081202 |