JP4954418B2 - 混濁媒質中の異常領域を位置特定する方法および装置 - Google Patents
混濁媒質中の異常領域を位置特定する方法および装置 Download PDFInfo
- Publication number
- JP4954418B2 JP4954418B2 JP2001571986A JP2001571986A JP4954418B2 JP 4954418 B2 JP4954418 B2 JP 4954418B2 JP 2001571986 A JP2001571986 A JP 2001571986A JP 2001571986 A JP2001571986 A JP 2001571986A JP 4954418 B2 JP4954418 B2 JP 4954418B2
- Authority
- JP
- Japan
- Prior art keywords
- wavelength
- value
- volume
- light
- turbid medium
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 50
- 230000002547 anomalous effect Effects 0.000 title claims description 5
- 238000005259 measurement Methods 0.000 claims description 36
- 239000008280 blood Substances 0.000 claims description 29
- 210000004369 blood Anatomy 0.000 claims description 29
- 230000005855 radiation Effects 0.000 claims description 27
- 238000003384 imaging method Methods 0.000 claims description 21
- 238000004364 calculation method Methods 0.000 claims description 18
- 230000002159 abnormal effect Effects 0.000 claims description 15
- 238000006392 deoxygenation reaction Methods 0.000 claims description 15
- 238000004590 computer program Methods 0.000 claims description 9
- 230000001678 irradiating effect Effects 0.000 claims 7
- 206010028980 Neoplasm Diseases 0.000 abstract description 17
- 210000000481 breast Anatomy 0.000 abstract description 12
- 230000003287 optical effect Effects 0.000 abstract description 2
- 238000009607 mammography Methods 0.000 abstract 1
- 108010054147 Hemoglobins Proteins 0.000 description 8
- 102000001554 Hemoglobins Human genes 0.000 description 8
- 108010064719 Oxyhemoglobins Proteins 0.000 description 3
- 230000017531 blood circulation Effects 0.000 description 3
- 239000003086 colorant Substances 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 206010006187 Breast cancer Diseases 0.000 description 1
- 208000026310 Breast neoplasm Diseases 0.000 description 1
- 238000005481 NMR spectroscopy Methods 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000005670 electromagnetic radiation Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 230000003211 malignant effect Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 230000036284 oxygen consumption Effects 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
- A61B5/0082—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence adapted for particular medical purposes
- A61B5/0091—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence adapted for particular medical purposes for mammography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
- A61B5/1455—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/43—Detecting, measuring or recording for evaluating the reproductive systems
- A61B5/4306—Detecting, measuring or recording for evaluating the reproductive systems for evaluating the female reproductive systems, e.g. gynaecological evaluations
- A61B5/4312—Breast evaluation or disorder diagnosis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/47—Scattering, i.e. diffuse reflection
- G01N21/4795—Scattering, i.e. diffuse reflection spatially resolved investigating of object in scattering medium
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Optics & Photonics (AREA)
- Biochemistry (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Reproductive Health (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Gynecology & Obstetrics (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Glass Compositions (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Description
【本発明が属する技術分野】
本発明は、体積要素の集合で表される混濁媒質中における異常領域の場所を位置特定する方法であって、混濁媒質が、主として第一の波長からなる放射を含む光によって照射され、その混濁領域を通る複数の光の経路に沿って透過する、主として第一の波長からなる放射を含む光の一部分の強度が測定されるステップと、その測定強度から混濁媒質の像を再構成するステップと含むものである方法に関する。
【0002】
また、本発明はそのような方法を実行するための装置に関する。
【0003】
本出願の文脈では、光という用語は、波長が約400nmと1400nmの間にある可視および赤外領域の電磁波放射を意味すると解されるべきである。混濁媒質とは、高度に光を分散させる材料からなる物質を意味すると解されるべきである。さらに厳密に言えば、本出願の文脈では、混濁媒質とは生物組織を意味すると解されるべきである。異常領域とは、周囲の領域における混濁媒質とは何らかの形で異質である混濁媒質中のある領域を意味するものと解されるべきである。さらに厳密に言えば、本出願の文脈では、そのような部位は、腫瘍組織を含む領域を意味すると解されるべきである。混濁媒質は、体積要素の集合によって表される。この種の体積要素は、体素(voxel)としても知られている。体積要素の寸法と形状は、全ての体積要素について同一であってよい。しかし、それに代えて、体積要素が、互いに異なった寸法と形状を持つことも可能である。
【従来の技術】
【0004】
この種類の方法と装置は、S.B. Colakらの「乳ガン検出のための臨床的光学的断層撮影法およびNMR分光法」、IEEE Journal of Selected Tops in Quantum Electronics、第5巻、第4号、1999年7/8月号より公知である。その公知の方法および装置は、生物組織の内部を画像化するのに使用される。その方法と装置は、特に、人間の女性または動物の雌の体の乳房組織に存在しているどのような腫瘍をも視覚的に位置特定するために、生体内乳房検診用医学的診断において使用することができる。この公知の方法によれば、混濁媒質は、様々な照射位置からの光により連続的に照射される。その後、その照射位置から広がる、混濁媒質を通る異なった光の経路に沿って透過した光の強度が測定される。測定強度は、混濁媒質の像の再構成に使用される。組織による光の減衰の空間的分布は、この像の中に再生される。光は、組織が光を分散しかつ吸収するので、減衰を受ける。
【0005】
異常領域中の組織による光の減衰が、周囲の組織による光の減衰より十分に逸脱していれば、その像の中で、異常である可能性のある領域が視覚的に位置特定することができる。
異なった波長の光を使用してその公知の方法を繰り返し実行することにより、複数の個々の像を形成することができる。
【0006】
この公知の方法および装置には、このようにして得られた像が、複数の異常領域を正確に視覚的に位置特定させるには鮮明度が十分でないという欠点がある。
【0007】
【課題を解決するための手段】
本発明の目的は、複数の異常領域が正確に、視覚的に位置特定させることができる、混濁媒質中の像を生成する方法を提供することである。
【0008】
この目的は、次の特徴を有する第一パラグラフで述べた種類の方法により達成される。すなわち、この方法は、第一の波長と同一でない主として第二の波長の放射を含む光によって照射され、その混濁媒質を通る複数の光の経路に沿って透過する主として第二の波長の光を含む光の一部分の強度が測定され、画像化ステップは以下のステップ、(1) 測定ステップにおける体積要素について測定される強度に基づいて、ある値が少なくとも第一のパラメータ(当該第一のパラメータは、混濁媒質の特性を表すものである)に割り当てられる第一の計算ステップと、(2) 重要性の値が、少なくとも第一のパラメータの値に基づいて、体積要素に割り当てられる第二の計算ステップと、(3) 体積要素の1つに対応する像中の各点について、混濁媒質の像中の対応する体積要素の有意性の値を表示する表示ステップとを含む。
【0009】
一方、本発明は、波長が異なる光を使った複数の測定から導出される測定強度を組み合わせることにより、組織の特性を表すパラメータに対する値を生成することができ、その特性によって、(個々に、または結合された状態で)正常組織を異常組織から見分けることが可能となるという認識に基づいている。このようにして入手されるパラメータは、特に、混濁媒質を構成する複数の要素を表すことができる。例えば、パラメータは、混濁媒質の端から端までの水量の空間分布を表すことが可能である。
【0010】
また、本発明は、得られたパラメータの値に基づいて異常であると考えられる領域を、像の中でマークすることができるという事実の認識に基づくものである。このようにして、ユーザ(例えば、放射線科医)の注意を、像中の異常と考えられる領域に向かわせることが出来る。
【0011】
混濁媒質中の像において、有意性の値は、例えば、色を利用することにより表現することができる。こうして、異なった有意性の値に対して、複数の別個の色が、割り当てられる。これに代えて、有意性の値が連続的に変わる場合には、なめらかに移行する複数の色が像に付与される。公知のカラースケール、例えば、レインボースケールを使用することができる。
【0012】
本発明によるその方法の一形態は、混濁媒質が、主として、第一の波長および第二の波長のいずれとも同一でない第三の波長の放射を含む光によって照射され、かつ、主として、混濁媒質を通る複数の光の経路に沿って透過する第三の波長の放射を含む光の一部分の強度が測定される第三の測定ステップを含み、かつ、第一の計算ステップでは、その測定ステップにおける体積要素について測定された強度に基づいて、ある値が、第一のパラメータによって表される特性と同一でない混濁媒質の特性を表す当該第二のパラメータに割り当てられ、さらに、第二のステップでは、第一のパラメータの値と第二のパラメータの値に基づき、有意性の値が体積に割り当てられることを特徴とする。
【0013】
この形態における各測定ステップは、測定データの集合を生成する。2つのパラメータに対する値は、得られた測定データの3つの集合から決定される。光の複数波長が、3つの測定ステップにおいて十分離れている場合には、測定データの複数集合が関連づけられることはなく、かつ公知の数学的技法を使って、全ての場合において、測定データの複数集合から互いに独立なパラメータに対する値を決定することが可能であろう。
【0014】
本発明による方法の一形態が特徴とする点は、第一の波長が、830nmと900nmの間の値をとり、第二の波長が、750nmと830nmの間の値をとり、かつ第三の波長が、655nmと750nmの間の値をとることである。上記の波長の光は、最新技術の半導体レーザにより容易に発生させることができる。
【0015】
3つの測定ステップに使われる光の波長が上記の領域に存在するときは、特に、ヘモグロビン(Hb)および酸化ヘモグロビン(HbO2)の空間分布を組織の端から端まで適切に近似させることが出来る。これは、生物組織による分散と吸収に起因する減衰がこれらの領域では小さく、かつ、ヘモグロビン(Hb)に対し、光の波長の関数である光の減衰が、酸化ヘモグロビン(HbO2)に対するそれと相違するからである。
【0016】
第一の波長が約867nmの値であり、これに対し、第二の波長が約780nmの値であり、かつ第三の波長が約715mの値、またはこれに代えて、約682nmの値であるとき、特に、魅力的な結果が得られる。
【0017】
組織中の水に対する脂肪の比率を正確に決定することが望まれる場合は、使用される光が1040nmと1070nmの間の波長、またはその代わりとしての、750nmと770nmの間の波長の放射を含んでいる第四の測定ステップが付加されても良い。脂肪対水の比率をより正確に決定することができると、組織の端から端までのヘモグロビン(Hb)と酸化ヘモグロビン(HbO2)の空間分布を決定することが可能になる。
【0018】
本発明による方法の、一つの有利な形態は、第一のパラメータが、体積要素の1つに対応する、混濁媒質中の体積の中の血液量を表し、第二のパラメータが、体積要素の1つに対応する混濁媒質の中の体積の中の血液の脱酸素反応比を表すことを特徴とする。
【0019】
本発明の文脈では、血液量とは、ヘモグロビン(Hb)と酸化ヘモグロビン(HbO2)の量の合計を意味するものと解されるべきであり、一方、血液の脱酸素反応比が比率とは、(Hb)/(Hb+HbO2)を意味するものと解されるべきである。
【0020】
本出願人によって実施された臨床試験において、腫瘍組織では、単位体積当たりの血液量が増加しており、さらに、腫瘍組織中の血液の脱酸素反応比は、正常な組織のそれと比較して、高いことが実証された。さらに、腫瘍が、より悪性であれば、血液量および脱酸素反応比は増加する。したがって、体積中の血液量および脱酸素反応比は、像のなかに少しでも腫瘍がある場合、その位置特定を行う本発明によるこの方法において、個別的に、または組み合わせて使用するための適切なパラメータである。
【0021】
本発明による方法の別の有利な一形態は、パラメータがパラメータ空間を定め、第二の計算ステップでは、問題とされる体積要素に対するパラメータに対して割り当てられる値に基づいて、体積要素にはパラメータ空間の位置が割り当てられ、体積要素に割り当てられる有意性の値が、パラメータ空間中の体積要素の位置により決定されることを特徴とする。
【0022】
パラメータ数に依存するが、このようなパラメータ空間は、1次元的(線)、2次元的(平面)、三次元立体的(空間または体積)、さらにそれ以上の多次元とすることが出来る。この場合、これらのパラメータは、パラメータ空間の座標系を表す。体積要素には、当該パラメータ空間中の位置が割り当てられ、当該位置の座標は、パラメータに関連した値によって形成される。例えば、二次元パラメータ空間では、第一のパラメータに割り当てられる値Xと第二のパラメータに割り当てられる値Yを有する体積要素には、座標(X、Y)を持つパラメータ空間中の位置が割り当てられるであろう。
【0023】
このようなパラメータ空間は、この方法を実行するのに使用される装置の中に、仮想的に存在させることが出来る。さらに、そのパラメータ空間は、ユーザに呈示することもできる。例えば、二次元パラメータ空間を、ディスプレイ画面上に表示することができる。
【0024】
本発明による方法の一形態は、パラメータ空間が少なくとも、第一の副空間と第二の副空間とに細分され、体積要素は、それらが第一の副空間を中の位置を占めているときには、第一の有意性の値を付与され、体積要素は、それらが第二の副空間を中の位置を占めているときには、第二の有意性の値を付与されることを特徴とする。
【0025】
パラメータ空間は、自動的にはもちろん、手動によっても副空間に細分することができる。例えば、ユーザは、ディスプレイ画面上に表示される二次元パラメータ空間における副空間を手動で指し示すことができ、次いで、当該副空間に位置する体積要素に、異常領域を示す有意性の値を割り当てることが出来る。ユーザは、例えば、パラメータ空間にわたる体積要素の位置の分布を参照分布と比較することにより、この副空間を選ぶことができる。このようにして、女性の体の乳房組織の検査の間に、一方の乳房について見いだされたパラメータ空間にわたる体積要素の位置の分布を、他方の乳房で見いだされたものと比較することができる。
【0026】
パラメータ空間を自動的に細分して副空間にすることは、例えば、パラメータ空間を所定の副空間に細分することによって行うことができる。こうして、多数の参照測定を行うことにより、副空間がパラメータ空間にわたってどのように分布すべきかについて、また副空間における体積要素にどの有意性が割り当てられるべきかについての決定を、前もって行うことが可能になる。あるいはこれに代えて、公知のセグメンテーション手法によって自動的細分を実現することも可能である。その場合には、例えば、第一の副空間に位置すべき体積要素の数の割合と第二の副空間に位置すべき体積要素の数の割合を、前もって知らせることができる。
【0027】
上述の説明では、パラメータ空間は2つの副空間に細分されるが、当業者には、より多くの副空間に細分することも可能であることは明白であろう。こうして、副空間のそれぞれについて、別個の有意性の値を、各副空間の中に位置する体積要素に割り当てることが可能になる。また、複数の副空間について、このような副空間内に位置する体積要素に対し同一の有意性の値を割り当てることも可能である。
【0028】
本発明による方法の一形態は、少なくとも、体積要素に割り当てられる位置とパラメータ空間中の選ばれたある位置との間の距離によって、有意性の値が決定されること特徴とする。例えば、ユークリッド幾何学の原理に準拠した2点間の距離を、距離の尺度として使用することができる。
【0029】
例えば、x1がパラメータxの第一の値、x2がその第二の値を表して、y1がパラメータyの第一の値、y2が第二の値を表すとしたとき、2次元のパラメータ空間に位置する2つの点P1(x1, y1)とP2(x2, y2)の間の距離は、
【式1】
に等しい。
【0030】
例えば、その選択位置として、パラメータ空間中の体積要素の位置の分布についての重心を、選択することは可能である。体積要素に割り当てられる有意性の値を、距離に依存して連続的に変化する値としても良い。しかしながら、その距離は、多くの離散的なステップに細分しても良い(そ¥この場合、有意性の値は各ステップ毎に関連する値となる)。
【0031】
本発明による方法のさらに別の形態は、特に、選択位置と相対的な体積要素に割り当てられた位置の状況により決定されること、を特徴とする。
【0032】
選択位置からの距離に加えて、選択位置と相対的な体積要素に割り当てられる位置の状況も、体積要素に割り当てられるべき有意性の値の観点で、決定的な意味を持つ。例えば、高い有意性の値を持った体積要素には、選択位置から同一の距離であるが低い有意性を持つ体積要素に割り当てられる有意性の値とは異なった値を割り当てても良い。
【0033】
本発明の目的は、 (1) その方法が、混濁媒質が、主として、第一の波長と同一でない第二の波長の放射を含む光によって照射され、混濁媒質を通る複数の光の経路に沿って通過させられ、主として第二の波長の放射を含む光の一部分の強度が測定される、第二の測定ステップを少なくとも含み、かつ、(2) 画像化ステップは、体積要素についての測定ステップにおいて測定された強度に基づいて、ある値を少なくとも、混濁媒質の特性を表す第一のパラメータに割り当てる第一の計算ステップと、少なくとも、体積要素について第一のパラメータに割り当てられた値から、混濁媒質の像を再構成するための再構成ステップとのステップを含むことを特徴とする、第一パラグラフで述べられた種類の方法によっても達成される。
【0034】
本発明は、異なる波長の光を使った複数の測定から得られる測定強度を組み合わせることにより、個々にまたは組み合わせて、正常組織を異常組織から区別することを可能にする組織特性を表すパラメータについての値を生成することができると言う事実の認識に基づいている。このようにして得られたパラメータは、特に、混濁媒質の複数の構成要素を表すことが出来る。
【0035】
このようなパラメータの画像化は、組織による光の減衰の像の場合よりも、より正確に異常領域が位置特定ができる像を生成する。これは、特に、パラメータが、混濁媒質の構成要素を表す場合が、これに相当する。例えば、女性の乳房組織における血液の脱酸素反応比、(Hb)/(Hb+HbO2)の空間分布の像においては、腫瘍組織の領域は、一般的に言って、腫瘍組織中の脱酸素反応比が通常の組織中のそれよりも高いので、視覚的に認知することが可能であろう。同様に、単位体積当たりの血液量(Hb+HbO2)の像においては、腫瘍組織における血液循環は通常の組織におけるそれよりも一般に大きいので、腫瘍組織を持った領域は、視覚的に認知することが可能であろう。
【0036】
測定強度から、組織による光の減衰の像を再構成することにも使われる再構成技法を適用して、1つの像を複数のパラメータの値から再構成することができる。このような再構成技法の例としては、Hoogenraadによる三次元代数的再構成技術と、Wachtersによる三次元共役勾配技術が挙げられる。
【0037】
ここで説明した方法は、異常領域を検出する目的で生物組織を画像化させるのに適しているばかりでなく、生物組織の機能的特性について画像化させるためにも適していることに留意すべきである。例えば、幼児の脳の中の血液循環や血液中の酸素含有量を、画像化させることができる。
【0038】
本発明のさらに別の目的は、異常領域を正確に視覚的に位置特定することのできる混濁媒質の像を作成するための装置を提供することである。
【0039】
この目的は、主として、第一の波長と同一でない第二の波長を有する放射を含む光によって混濁媒質を照射するための手段と、混濁媒質を通る複数の光の経路に沿って透過する、主として第二の波長の放射を含む光の一部分の強度を測定するための手段と、体積要素についての測定強度に基づいて、ある値を、少なくとも、混濁媒質の性質を表す第一のパラメータに割り当てるための手段と、少なくとも、第一のパラメータの値に基づいて、有意性の値を体積要素に割り当てるための手段と、体積要素の1つに対応する像中の各点について、混濁媒質の像中の対応する体積要素の有意性の値を表示するための手段をも含むことを特徴とする第一パラグラフで述べた種類の装置によって達成される。
【0040】
この目的は、主として、第一の波長と同一でない第二の波長の放射を含む光によって混濁媒質を照射するための手段と、混濁媒質を通る複数の光の経路に沿って透過する、主として第二の波長の放射を含む光の一部分の強度を測定するための手段と、体積要素についての測定強度に基づいて、ある値を、少なくとも、混濁媒質の特性を表す第一のパラメータに割り当てる手段と、少なくとも、体積要素に対する第一のパラメータに割り当てられた値から、混濁媒質の像を再構成するための手段をも含むことを特徴とする第一パラグラフで述べられた種類の装置によって達成される。
【0041】
前のパラグラフで述べた装置は、例えば、医学的診断ワークステーションの形で実施される。
【0042】
本発明は、さらに、本発明によるその方法を実行するためのコンピュータプログラムに関する。このコンピュータプログラムは、本発明によるその方法の全ての形態を実行するための命令を含むものでも良く、またはこれに代えて、この方法の複数形態の部分集合を実行するための命令を含むものでも良い。このコンピュータプログラムの諸命令は、例えば、医学的診断ワークステーションのメモリにロードされ、次に、このワークステーションのプロセッサによって実行される。本発明のさらに別の目的によれば、コンピュータプログラムは、例えば、CD-ROMや磁気テープなどの記録担体上に記録される。コンピュータプログラムの諸命令は、記録担体から、例えば、医学的診断ワークステーションのメモリに転送される。あるいはこれに代えて、コンピュータプログラムの命令は、記録担体からPCまたはその他の標準的コンピュータシステムのメモリに転送される。また、コンピュータプログラムを、例えば、「ワールドワイドウェッブ」などのネットワークを介してメモリにロードすることも可能である。
【0043】
【発明を実施するための形態】
本発明の上記および他のより詳細な観点は、以下に、具体例により、図面を参照して詳細に説明される。
【0044】
図1は、本発明による第一の方法を線図的に示す。第一の測定ステップの間では、混濁媒体が、主として867nmの波長の放射を含む光によって照射され、混濁媒質を通る複数の光の経路を透過した光の一部分の強度が測定される。このような測定ステップは、欧州特許出願98925884.3 (PHN 16.442)に、より詳細に説明されている。その後、第二の同様な測定ステップが、主として780nmの波長の放射を含む光によって行われ、これに続き、主として波長715nmの放射を含む光を使った第三の測定ステップが行われる。
【0045】
画像化ステップ4における再構成アルゴリズム5を使って、混濁媒質の像9が、第三の測定するステップ3で測定された強度から再構成される。本例では、像は、第三の測定ステップ3で測定された強度から再構成されるのであるが、第一の測定ステップ1または第二の測定ステップ2で測定された強度はこの目的にも使用することができる。
【0046】
第一の計算ステップ6では、それぞれの体積要素について、2つのパラメータが、3つの測定ステップ中で測定された強度の3つの集合から決定される。これらのパラメータは、検査される組織の局所的な特性を表すものであるが、本実施例では、体積中の血液量と体積中の血液の脱酸素反応比である。この文脈では、体積中の血液量とは、その体積中に存在するヘモグロビンと酸化ヘモグロビンの合計量(Hb+HbO2)を意味するものと解すべきである。血液の脱酸素反応比なる用語は、ここでは、その体積中に存在するヘモグロビンの、体積中に存在する血液量に対する比、(Hb)/(Hb+HbO2)を意味するものと解されるべきである。腫瘍は、一般に、より高い血液循環とより高い酸素消費によって特徴付けられ、両者のパラメータを増加させるから、これらのパラメータは、乳房組織の腫瘍組織を位置特定させるのに特に適している。
【0047】
第二の計算ステップ7では、体積要素が、二次元パラメータ空間20中に写像される。この具体例は図2に示される。このパラメータ空間20は、第一の計算ステップ6で決定された2つのパラメータによって定められる。本実施例の場合、体積中の血液量は、水平のX軸22により表され、体積中の血液の脱酸素反応比は、垂直のY軸21によって表される。このパラメータ空間中の体積要素には、問題とされる体積要素に対する2つのパラメータの値に依存する位置23が割り当てられる。こうして、体積中の血液量について値X1を有し、体積中の血液の脱酸素反応比について値Y1を有する体積要素には、パラメータ空間に座標(X1, Y1)を持っている位置が割り当てられるであろう。
【0048】
次に、各体積要素には、パラメータ空間20中の体積要素の位置23に依存する有意性の値が、割り当てられる。本実施例では、体積要素が異常領域の中に位置すると考えられる場合には、有意性の値「1」が割り当てられ、体積要素が異常領域の外の領域に位置すると考えられる場合には、有意性の値「0」が割り当てられる。図2Aは、両者のパラメータについて高い値を含み、かつ上で議論した腫瘍組織を示す副空間24を示す。副空間24の中に位置25を持っている体積要素には、有意性の値「1」が、割り当てられる。これ以外の体積要素には、有意性の値「0」が、割り当てられる。本実施例では、2つの異なった有意性の値しか使用されていないが、パラメータ空間が複数の副空間に細分される場合には、複数の離散的な有意性の値を使用することができる。連続的に変化する有意性の値を使用することも、可能である。その場合には、有意性の値は、例えば、確率値を表す。
【0049】
次に、混濁媒質の像9中の体積要素に割り当てられる有意性の値が、表示ステップ8において再生される。体積要素に対応する像9中のドットには、例えば、問題とされる体積要素に割り当てられる有意性の値が「1」の場合は赤色が、問題とされる体積要素に割り当てられる有意性の値が「0」の場合は灰色が、割り当てられる。こうして、像9は、灰色の値を含む混濁媒質と赤い値を有する領域10とにより形成される。領域10は混濁媒質中の異常領域に対応している。他の色や他の方法も、像の中の領域をマークするのに使い得ることは、当業者には明らかであろう。
【0050】
図2は、二次元パラメータ空間20の例を示す。パラメータ空間20は、2つのパラメータによって定められる。第一のパラメータは、水平のX軸22によって表され、第二のパラメータは、垂直のY軸21によって表される。このパラメータ空間中の体積要素には、問題とされる体積要素に対する2つのパラメータの値に依存する位置23が、割り当てられる。こうして、第一のパラメータにX1を持ち、持っている体積要素、第二のパラメータにY1を持つ体積要素には、パラメータ空間において座標(X1, Y1)を持つ位置が割り当てられるであろう。
【0051】
図2Aでは、パラメータ空間は、副空間24の手動生成により2つの副空間に細分される。副空間24の中に位置する体積要素25には、副空間24の外に位置する体積要素に割り当てられる有意性の値から外れた有意性の値が割り当てられた。乳房組織の画像化の間に、一方の乳房について行われた測定から得られたパラメータ空間を他方の乳房について行われた測定から得られたれスペース空間と比較することで、副空間24の位置と形状を選択することができる。
【0052】
図2Bは、有意性の値を体積要素に割当てる別の方法を示す。ここでは、体積要素には、パラメータ空間20の体積要素の位置30と選択位置26の間の距離27に依存する有意性の値が割り当てられる。現在の例では、選択位置26は、パラメータ空間20における体積要素の位置の分布についての幾何学的な重心である。割り当てられるべき有意性の値は、距離27について連続関数であっても良い。これに代えて、距離27を、多くの離散的なステップに細分し、対応する有意性の値をそれぞれのステップに持たせても良い。
【0053】
例え、第一の体積要素の位置30と第二の体積要素の位置31が選択位置26から同じ距離に位置していても、異なった有意性の値を2つの体積要素に割当てることも可能である。この場合、選択された点26に対する体積要素の位置の方位も、重要である。例えば、パラメータ空間20の右側上の象限28に位置する位置30を有する体積要素には、パラメータ空間20の左側下の象限29中に位置する位置31を有する体積要素に割り当てられる有意性の値とは異なる有意性の値が、割り当てられる。例えば、体積中の血液量と体積中の血液の脱酸素反応比が、パラメータとして選ばれるときには、このことは意味がある。すでに説明したように、両方のパラメータの増加した値は、腫瘍組織を暗示している。これは、右側上の象限28に位置する体積要素についての場合である。したがって、これらの体積雄要素には、高い有意性の値が割り当てられるべきである。しかし、1つまたは両方のパラメータが減少した値は、有意性を持たない。したがって、その他の3つの象限の1つに存する体積要素には、これら2つのパラメータに基づいて高い有意性の値を割り当ててはならない。
【0054】
図3は、本発明による第二の方法を線図的に説明するものである。第一の測定ステップ1の間では、混濁媒質が、主として波長867nmの放射を含む光によって照射され、その混濁媒質を通る複数の光の経路に沿って透過する光の一部の強度が、測定される。続いて、第二の、同様な計算ステップ2が、主として波長780nmの放射を含む光を使って実行され、さらに、主として715nmの波長の放射を含む光を使う第三のステップ3が実行される。
【0055】
第一の計算ステップ6において、3つの測定点から測定した強度の3つの集合から2つのパラメータが、それぞれの体積要素に対して、決定される。続いて、再構成ステップ13において、体積要素に対する第一のパラメータの値が像11において画像化され、第二のパラメータの値が像12において画像化される。この場合、像11中の一点と像12における一点は、体積要素に対応する。これは、混濁媒質にわたるパラメータの空間分布の像を生成する。
【0056】
パラメータが適切に選ばれる場合には、異常領域が、像の中で明確に視認されるであろう。例えば、体積中の血液量が像11において画像化され、かつ体積中の血液の脱酸素反応比が、像12において画像化されるときは、両方のパラメータは、腫瘍組織中の値を増加させるので、乳房組織中に存在する腫瘍組織を有する領域は、目に見えるであろう。
【0057】
現在の具体例では、2つのパラメータが別々に画像化される。しかし、2つのパラメータの内の1つのみを画像化させることも可能である。さらに、3つ以上のパラメータが、第一の計算ステップ6において決定されるときは、3つ以上のパラメータを別々に画像化させることも可能である。更に、複数のパラメータの複数の像が、1つの共通の像を形成するように結合させることも可能である。
【図面の簡単な説明】
【図1】本発明による第一の方法を線図的に説明する図である。
【図2A】二次元パラメータ空間についての一例を示す図である。
【図2B】二次元パラメータ空間についての一例を示す図である。
【図3】本発明による第二の方法を線図的に説明する図である。
【符号の説明】
1 第一の測定ステップ
2 第二の測定ステップ
3 第三の測定ステップ
4 画像化ステップ
5 再構成アルゴリズム
6 第一の計算ステップ
7 第二の計算ステップ
8 表示ステップ
9 像
10 領域
11 像
12 像
13 再構成ステップ
20 二次元パラメータ空間
23 位置
24 副空間
25 位置
26 選択位置
27 距離
30 第一の体積要素の位置
31 第二の体積要素の位置
Claims (8)
- 体積要素の集合により表される、混濁媒質中の異常領域を画像化する装置の作動方法であって、前記装置が、第一の照射手段と、第一の測定手段と、再構成手段と、第二の照射手段と、第二の測定手段と、第三の照射手段と、第三の測定手段と、第一の計算手段と、第二の計算手段と、表示手段とを有し、前記方法が、
・ 前記第一の照射手段により、主として第一の波長の放射を含む光を照射し、前記第一の測定手段により、前記混濁媒質を通る複数の光の経路に沿って透過する、主として前記第一の波長の放射を含む前記光の一部分の強度を測定する第一の測定ステップと、
・ 前記再構成手段により、前記測定強度から前記混濁媒質の像を再構成する画像化ステップとを
含む方法において、
・ 前記方法が、前記第二の照射手段により、主として、前記第一の波長と同一でない第二の波長の放射を含む光を照射し、前記第二の測定手段により、前記混濁媒質を通る複数の光の経路に沿って透過する、主として前記第二の波長の放射を含む前記光の一部分の強度を測定する第二の測定ステップを含み、かつ
・ 前記第三の照射手段により、主として、前記第一の波長と同一でなくかつ前記第二の波長と同一でない第三の波長の放射を含む光を照射し、前記第三の測定手段により、主として、前記混濁媒質を通る複数の光の経路に沿って透過する前記第三の波長の光を含む前記光の一部分の強度を測定する第三の測定ステップを含み、
・ 前記画像化ステップが、
- 前記第一の計算手段により、それぞれの体積要素に対して、前記測定ステップにおける体積要素について測定された前記強度に基づいて、前記混濁媒質の性質を表す第一のパラメータの値及び前記第一のパラメータによって表されたものと同一でない前記混濁媒質の特性を表す第二のパラメータの値を決定する第一の計算ステップと、
-前記第二の計算手段により、前記それぞれの体積要素に対して、前記第一のパラメータの前記値及び前記第二のパラメータの前記値が位置するパラメータ空間内の副空間を示す値を決定する第二の計算ステップと、
前記表示手段により、前記副空間を示す値に基づいて、前記混濁媒質の前記像において前記副空間を示す値を持つ体積要素に対応する領域を表示するための表示ステップとを含み、
前記第一の波長が、830nmと900nmの間の値をとり、前記第二の波長が、750nmから830nmの間の値をとり、かつ前記第三の波長が、655nmと750nmの間の値をとり、
前記第一のパラメータが、前記体積要素の1つに対応する前記混濁媒質中のある体積中の血液量を表し、そして
前記第二のパラメータが、前記体積要素の1つに対応する前記混濁媒質中のある体積中の血液の脱酸素反応比を表わす、
ことを特徴とする方法。 - 複数のパラメータが、パラメータ空間を定め、
前記第二の計算ステップにおいては、前記体積要素に、前記問題とされる体積要素に対する前記パラメータの前記値に基づいて、前記パラメータ空間中の位置が、割り当てられ、そして
前記それぞれの体積要素に対して決定される副空間を示す値が、前記パラメータ空間中の前記体積要素の前記位置によって決定される
ことを特徴とする請求項1に記載の方法。 - 前記パラメータ空間が、少なくとも第一の副空間と第二の副空間に細分され、
前記体積要素には、それらが前記第一の副空間中の位置を占有するときに、第一の副空間を示す値が割り当てられ、そして
前記体積要素には、それらが前記第二の副空間中の位置を占有するときに、第二の副空間を示す値が割り当てられる
ことを特徴とする請求項2に記載の方法。 - 前記副空間を示す値が、少なくとも、前記体積要素に割り当てられる前記位置と前記パラメータ空間中のある選択位置との間の距離によって決定される
ことを特徴とする請求項2に記載の方法。 - 前記副空間を示す値が、特に、前記選択位置と相対的な前記体積要素に割り当てられる前記位置の状況により決定される
ことを特徴とする請求項4に記載の方法。 - 体積要素の集合により表される、混濁媒質中の異常領域を画像化する装置であって、前記装置が、
主として第一の波長の放射を含む光によって前記混濁媒質を照射する第一の照射手段と、
前記混濁媒質を通る光の経路に沿って透過する、主として前記第一の波長の放射を含む前記光の一部分の強度を測定する第一の測定手段と、
前記測定強度から前記混濁媒質の像を再構成する再構成手段とを含む装置において、
主として、前記第一の波長と同一でない第二の波長の光を含む光によって前記混濁媒質を照射する第二の照射手段と、
主として、前記第一の波長と同一でなくかつ前記第二の波長と同一でない第三の波長の光を含む光によって前記混濁媒質を照射する第三の照射手段と、
前記混濁媒質を通る複数の光の経路に沿って透過した、主として前記第二の波長の放射を含む前記光の一部分の強度を測定する第二の測定手段と、
前記混濁媒質を通る複数の光の経路に沿って透過した、主として前記第三の波長の放射を含む前記光の一部分の強度を測定する第三の測定手段と、
それぞれの体積要素に対して、体積要素についての前記測定強度に基づいて、少なくとも、前記混濁媒質の特性を表す第一のパラメータの値及び前記第一のパラメータによって表されたものと同一でない前記混濁媒質の特性を表す第二のパラメータの値を決定する第一の計算手段と、
前記それぞれの体積要素に対して、前記第一のパラメータの前記値及び前記第二のパラメータの前記値が位置するパラメータ空間内の副空間を示す値を決定する第二の計算手段と、
前記副空間を示す値に基づいて、前記混濁媒質の前記像において前記副空間を示す値を持つ体積要素に対応する領域を表示する表示手段とを
含み、
前記第一の波長が、830nmと900nmの間の値をとり、前記第二の波長が、750nmから830nmの間の値をとり、かつ前記第三の波長が、655nmと750nmの間の値をとり、
前記第一のパラメータが、前記体積要素の1つに対応する前記混濁媒質中のある体積中の血液量を表し、そして
前記第二のパラメータが、前記体積要素の1つに対応する前記混濁媒質中のある体積中の血液の脱酸素反応比を表わす、
むことを特徴とする装置。 - 請求項6に記載の装置がプロセッサを有し、前記プロセッサにより実行される場合に請求項1〜5のいずれかに記載の前記方法を実行するための命令の集合を含むコンピュータプログラム。
- 請求項7に記載の前記コンピュータプログラムを含む記録担体。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP00201164.1 | 2000-03-31 | ||
EP00201164 | 2000-03-31 | ||
PCT/EP2001/003362 WO2001074238A1 (en) | 2000-03-31 | 2001-03-22 | Method and device for localizing a deviant region in a turbid medium |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2003528678A JP2003528678A (ja) | 2003-09-30 |
JP2003528678A5 JP2003528678A5 (ja) | 2011-04-14 |
JP4954418B2 true JP4954418B2 (ja) | 2012-06-13 |
Family
ID=8171282
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2001571986A Expired - Fee Related JP4954418B2 (ja) | 2000-03-31 | 2001-03-22 | 混濁媒質中の異常領域を位置特定する方法および装置 |
Country Status (7)
Country | Link |
---|---|
US (2) | US6718195B2 (ja) |
EP (1) | EP1272102B1 (ja) |
JP (1) | JP4954418B2 (ja) |
CN (1) | CN1262235C (ja) |
AT (1) | ATE422325T1 (ja) |
DE (1) | DE60137623D1 (ja) |
WO (1) | WO2001074238A1 (ja) |
Families Citing this family (38)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4606543B2 (ja) * | 2000-04-13 | 2011-01-05 | パナソニック株式会社 | 光学特性計測装置における被検溶液量確認方法および計測系制御方法 |
US7983759B2 (en) * | 2002-12-18 | 2011-07-19 | Cardiac Pacemakers, Inc. | Advanced patient management for reporting multiple health-related parameters |
JP4299059B2 (ja) * | 2003-05-30 | 2009-07-22 | 株式会社 日立ディスプレイズ | 有機エレクトロルミネッセンス表示装置の製造方法 |
US7693318B1 (en) | 2004-01-12 | 2010-04-06 | Pme Ip Australia Pty Ltd | Method and apparatus for reconstruction of 3D image volumes from projection images |
US7120283B2 (en) * | 2004-01-12 | 2006-10-10 | Mercury Computer Systems, Inc. | Methods and apparatus for back-projection and forward-projection |
CN1997998B (zh) | 2004-03-05 | 2012-12-12 | 皇家飞利浦电子股份有限公司 | 交互式计算机辅助诊断 |
US20050270298A1 (en) * | 2004-05-14 | 2005-12-08 | Mercury Computer Systems, Inc. | Daughter card approach to employing multiple graphics cards within a system |
US8189002B1 (en) | 2004-10-29 | 2012-05-29 | PME IP Australia Pty, Ltd. | Method and apparatus for visualizing three-dimensional and higher-dimensional image data sets |
US7778392B1 (en) | 2004-11-02 | 2010-08-17 | Pme Ip Australia Pty Ltd | Method of reconstructing computed tomography (CT) volumes suitable for execution on commodity central processing units (CPUs) and graphics processors, and apparatus operating in accord with those methods (rotational X-ray on GPUs) |
JP4512471B2 (ja) * | 2004-11-10 | 2010-07-28 | 株式会社日立ハイテクノロジーズ | 走査型電子顕微鏡及び半導体検査システム |
US7609884B1 (en) | 2004-12-23 | 2009-10-27 | Pme Ip Australia Pty Ltd | Mutual information based registration of 3D-image volumes on GPU using novel accelerated methods of histogram computation |
JP4599191B2 (ja) * | 2005-03-01 | 2010-12-15 | 国立大学法人神戸大学 | 画像診断処理装置および画像診断処理プログラム |
US7623732B1 (en) | 2005-04-26 | 2009-11-24 | Mercury Computer Systems, Inc. | Method and apparatus for digital image filtering with discrete filter kernels using graphics hardware |
JP5451072B2 (ja) * | 2005-09-13 | 2014-03-26 | コーニンクレッカ フィリップス エヌ ヴェ | 空間分解オキシメトリ |
US20080278727A1 (en) * | 2005-11-18 | 2008-11-13 | Koninklijke Philips Electronics, N.V. | Device For Imaging an Interior of a Turbid Medium |
WO2007060586A2 (en) | 2005-11-23 | 2007-05-31 | Koninklijke Philips Electronics N.V. | Method and device for optical imaging of a turbid medium |
US7750331B2 (en) * | 2005-11-23 | 2010-07-06 | Koninklijke Philips Electronics N.V. | Method and device for imaging an interior of a turbid medium using an amplification factor selected from an estimate of expected electrical signal strength |
JP2009529948A (ja) * | 2006-03-17 | 2009-08-27 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 混濁媒体画像形成装置 |
EP1943941A1 (en) * | 2006-12-21 | 2008-07-16 | Koninklijke Philips Electronics N.V. | Method for optically imaging an interior of a turbid medium, method for reconstructing an image of an interior of a turbid medium, device for imaging an interior of a turbid medium, medical image acquisition device and computer program products for use in said methods and devices |
US8019151B2 (en) * | 2007-06-11 | 2011-09-13 | Visualization Sciences Group, Inc. | Methods and apparatus for image compression and decompression using graphics processing unit (GPU) |
US8392529B2 (en) | 2007-08-27 | 2013-03-05 | Pme Ip Australia Pty Ltd | Fast file server methods and systems |
US10311541B2 (en) | 2007-11-23 | 2019-06-04 | PME IP Pty Ltd | Multi-user multi-GPU render server apparatus and methods |
US9904969B1 (en) | 2007-11-23 | 2018-02-27 | PME IP Pty Ltd | Multi-user multi-GPU render server apparatus and methods |
US9019287B2 (en) | 2007-11-23 | 2015-04-28 | Pme Ip Australia Pty Ltd | Client-server visualization system with hybrid data processing |
WO2011065929A1 (en) | 2007-11-23 | 2011-06-03 | Mercury Computer Systems, Inc. | Multi-user multi-gpu render server apparatus and methods |
WO2009067680A1 (en) | 2007-11-23 | 2009-05-28 | Mercury Computer Systems, Inc. | Automatic image segmentation methods and apparartus |
KR20110113188A (ko) * | 2009-01-23 | 2011-10-14 | 코닌클리케 필립스 일렉트로닉스 엔.브이. | 탁한 매체 내에 적어도 부분적으로 삽입되도록 구성된 광학적 검사 장치 |
WO2011132128A1 (en) * | 2010-04-21 | 2011-10-27 | Koninklijke Philips Electronics N.V. | Determination of a lipid water ratio |
WO2012128614A1 (en) * | 2011-03-24 | 2012-09-27 | Erasmus University Medical Center Rotterdam | Method to determine the absorption coefficient in turbid media |
US9509802B1 (en) | 2013-03-15 | 2016-11-29 | PME IP Pty Ltd | Method and system FPOR transferring data to improve responsiveness when sending large data sets |
US10540803B2 (en) | 2013-03-15 | 2020-01-21 | PME IP Pty Ltd | Method and system for rule-based display of sets of images |
US11183292B2 (en) | 2013-03-15 | 2021-11-23 | PME IP Pty Ltd | Method and system for rule-based anonymized display and data export |
US8976190B1 (en) | 2013-03-15 | 2015-03-10 | Pme Ip Australia Pty Ltd | Method and system for rule based display of sets of images |
US10070839B2 (en) | 2013-03-15 | 2018-09-11 | PME IP Pty Ltd | Apparatus and system for rule based visualization of digital breast tomosynthesis and other volumetric images |
US11244495B2 (en) | 2013-03-15 | 2022-02-08 | PME IP Pty Ltd | Method and system for rule based display of sets of images using image content derived parameters |
US11599672B2 (en) | 2015-07-31 | 2023-03-07 | PME IP Pty Ltd | Method and apparatus for anonymized display and data export |
US9984478B2 (en) | 2015-07-28 | 2018-05-29 | PME IP Pty Ltd | Apparatus and method for visualizing digital breast tomosynthesis and other volumetric images |
US10909679B2 (en) | 2017-09-24 | 2021-02-02 | PME IP Pty Ltd | Method and system for rule based display of sets of images using image content derived parameters |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH01174940A (ja) * | 1987-12-29 | 1989-07-11 | Shimadzu Corp | ヘモグロビン測定法 |
JPH04138342A (ja) * | 1990-09-29 | 1992-05-12 | Shimadzu Corp | 光診断装置 |
JPH06142087A (ja) * | 1992-11-13 | 1994-05-24 | Omron Corp | 運動モニタ装置 |
JPH0810244A (ja) * | 1994-04-30 | 1996-01-16 | Shimadzu Corp | 光散乱・吸収体の光学的測定装置 |
WO1999002956A2 (en) * | 1997-07-07 | 1999-01-21 | Cedars-Sinai Medical Center | Method and devices for laser induced fluorescence attenuation spectroscopy (lifas) |
WO1999003394A1 (en) * | 1997-07-15 | 1999-01-28 | Koninklijke Philips Electronics N.V. | Method of localizing an object in a turbid medium |
WO1999018845A1 (en) * | 1997-10-10 | 1999-04-22 | Massachusetts Institute Of Technology | Method for measuring tissue morphology |
WO2000015101A1 (en) * | 1998-09-11 | 2000-03-23 | Spectrx, Inc. | Multi-modal optical tissue diagnostic system |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5800354A (en) * | 1994-11-23 | 1998-09-01 | U.S. Phillips Corporation | Method of and device for magnetic resonance imaging |
DE69616743T2 (de) * | 1995-09-11 | 2002-07-25 | Koninkl Philips Electronics Nv | Verfahren und vorrichtung zum bildhaften darstellen des inneren eines trüben mediums |
EP0963174B1 (en) * | 1997-11-22 | 2004-03-17 | Koninklijke Philips Electronics N.V. | Method of localizing an object in a turbid medium |
WO1999066832A1 (en) * | 1998-06-25 | 1999-12-29 | Koninklijke Philips Electronics N.V. | Method of localizing an object in a turbid medium |
WO2001046713A1 (en) * | 1999-12-21 | 2001-06-28 | Koninklijke Philips Electronics N.V. | Ultrasonic image processing method and examination system for displaying an ultrasonic composite image sequence of an artery |
US6533727B1 (en) * | 2000-02-11 | 2003-03-18 | Koninklijke Phillips Electronics N.V. | Ultrasonic system and method employing non-integer harmonic echo signals for imaging |
-
2001
- 2001-03-22 DE DE60137623T patent/DE60137623D1/de not_active Expired - Lifetime
- 2001-03-22 JP JP2001571986A patent/JP4954418B2/ja not_active Expired - Fee Related
- 2001-03-22 EP EP01929448A patent/EP1272102B1/en not_active Expired - Lifetime
- 2001-03-22 AT AT01929448T patent/ATE422325T1/de not_active IP Right Cessation
- 2001-03-22 CN CNB018014089A patent/CN1262235C/zh not_active Expired - Fee Related
- 2001-03-22 WO PCT/EP2001/003362 patent/WO2001074238A1/en active Application Filing
- 2001-03-28 US US09/819,282 patent/US6718195B2/en not_active Expired - Lifetime
-
2004
- 2004-01-20 US US10/760,872 patent/US6922582B2/en not_active Expired - Lifetime
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH01174940A (ja) * | 1987-12-29 | 1989-07-11 | Shimadzu Corp | ヘモグロビン測定法 |
JPH04138342A (ja) * | 1990-09-29 | 1992-05-12 | Shimadzu Corp | 光診断装置 |
JPH06142087A (ja) * | 1992-11-13 | 1994-05-24 | Omron Corp | 運動モニタ装置 |
JPH0810244A (ja) * | 1994-04-30 | 1996-01-16 | Shimadzu Corp | 光散乱・吸収体の光学的測定装置 |
WO1999002956A2 (en) * | 1997-07-07 | 1999-01-21 | Cedars-Sinai Medical Center | Method and devices for laser induced fluorescence attenuation spectroscopy (lifas) |
JP2001509589A (ja) * | 1997-07-07 | 2001-07-24 | シーダーズ−サイナイ・メディカル・センター | レーザー誘起蛍光減衰分光分析の方法と装置 |
WO1999003394A1 (en) * | 1997-07-15 | 1999-01-28 | Koninklijke Philips Electronics N.V. | Method of localizing an object in a turbid medium |
JP2001500978A (ja) * | 1997-07-15 | 2001-01-23 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 混濁媒体内で対象を位置決めする方法 |
WO1999018845A1 (en) * | 1997-10-10 | 1999-04-22 | Massachusetts Institute Of Technology | Method for measuring tissue morphology |
WO2000015101A1 (en) * | 1998-09-11 | 2000-03-23 | Spectrx, Inc. | Multi-modal optical tissue diagnostic system |
Also Published As
Publication number | Publication date |
---|---|
JP2003528678A (ja) | 2003-09-30 |
EP1272102A1 (en) | 2003-01-08 |
US6922582B2 (en) | 2005-07-26 |
US20010026848A1 (en) | 2001-10-04 |
DE60137623D1 (de) | 2009-03-26 |
CN1262235C (zh) | 2006-07-05 |
US20040147039A1 (en) | 2004-07-29 |
EP1272102B1 (en) | 2009-02-11 |
US6718195B2 (en) | 2004-04-06 |
WO2001074238A1 (en) | 2001-10-11 |
ATE422325T1 (de) | 2009-02-15 |
CN1380821A (zh) | 2002-11-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4954418B2 (ja) | 混濁媒質中の異常領域を位置特定する方法および装置 | |
US11656448B2 (en) | Method and apparatus for quantitative hyperspectral fluorescence and reflectance imaging for surgical guidance | |
EP1054618B1 (en) | Detection, imaging and characterization of breast tumors | |
Yates et al. | Optical tomography of the breast using a multi-channel time-resolved imager | |
Enfield et al. | Three-dimensional time-resolved optical mammography of the uncompressed breast | |
US5590215A (en) | Method for providing medical images | |
EP1210910A1 (en) | Color parametric and composite maps for perfusion CT | |
US8886284B2 (en) | Devices and methods for combined optical and magnetic resonance imaging | |
JP2020513978A (ja) | 心筋ct灌流画像合成 | |
WO2007035799A2 (en) | Apparatus and method to acquire data for reconstruction of images pertaining to functional and anatomical structure of the breast | |
Busch et al. | Computer aided automatic detection of malignant lesions in diffuse optical mammography | |
Baikejiang et al. | Diffuse optical tomography for breast cancer imaging guided by computed tomography: A feasibility study | |
Mastanduno et al. | Adaptable near-infrared spectroscopy fiber array for improved coupling to different breast sizes during clinical MRI | |
US8620051B2 (en) | Registration of optical images of turbid media | |
Momtahen et al. | A machine learning approach: NIR scattering data analysis for breast cancer detection and classification | |
US20230280577A1 (en) | Method and apparatus for quantitative hyperspectral fluorescence and reflectance imaging for surgical guidance | |
Goh et al. | Ultrasound guided optoacoustic tomography in assessment of tumor margins for lumpectomies | |
Zhang et al. | Direct regularization from co-registered contrast MRI improves image quality of MRI-guided near-infrared spectral tomography of breast lesions | |
JP2004177408A (ja) | 光を放出する領域の位置測定方法および装置 | |
Barbour et al. | Time-frequency analysis of functional optical mammographic images | |
Brooksby | Combining near infrared tomography and magnetic resonance imaging to improve breast tissue chromophore and scattering assessment | |
EP3484368B1 (en) | Identifying types of adipose tissue | |
Khayat et al. | Optical tomography as adjunct to x-ray mammography: methods and results | |
Enfield et al. | Sensitivity and specificity of 3D optical mammography | |
Xu | High-Resolution Diffuse Optical Tomography in Optically-Thick Media |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080319 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20080319 |
|
RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7422 Effective date: 20080319 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20100824 |
|
A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20101122 |
|
A602 | Written permission of extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A602 Effective date: 20101130 |
|
A524 | Written submission of copy of amendment under article 19 pct |
Free format text: JAPANESE INTERMEDIATE CODE: A524 Effective date: 20110224 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20110719 |
|
A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20111018 |
|
A602 | Written permission of extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A602 Effective date: 20111025 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20120118 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20120216 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20120314 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 4954418 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20150323 Year of fee payment: 3 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |