CN1380821A - 用于探测在混浊介质中异常区域的方法和设备 - Google Patents
用于探测在混浊介质中异常区域的方法和设备 Download PDFInfo
- Publication number
- CN1380821A CN1380821A CN01801408A CN01801408A CN1380821A CN 1380821 A CN1380821 A CN 1380821A CN 01801408 A CN01801408 A CN 01801408A CN 01801408 A CN01801408 A CN 01801408A CN 1380821 A CN1380821 A CN 1380821A
- Authority
- CN
- China
- Prior art keywords
- turbid medium
- volume element
- parameter
- wavelength
- value
- 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.)
- Granted
Links
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)
- Immunology (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Gynecology & Obstetrics (AREA)
- Reproductive Health (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)
Abstract
本发明涉及在混浊介质(turbid medium)中探测异常区域的方法。本发明还涉及用于执行这种方法的设备。所述方法可用于光学乳房照相术中,该照相术中,用光检查女性身体的胸部。所述方法产生图像,图像中,可清晰识别任何异常,如肿瘤。这一点尤其是通过在混浊介质的图像(9)中提供标记(10)来实现。
Description
本发明涉及一种探测在由一组体积元素所代表的混浊介质中异常区域的方法,该方法包括:第一测量步骤,其中,由主要包括第一波长的辐射的光来照射混浊介质,测量主要包括第一波长的辐射的沿着通过混浊介质的多个光路传播的一部分光的强度;成像步骤,用于根据测量的强度来重组混浊介质的图像。
本发明还涉及用于执行这种方法的设备。
在本申请的行文中,术语光应理解为指的是波长在大约400和1400nm之间的可见或红外范围内的电磁辐射。混浊介质应理解为指的是由高散光材料构成的物质。尤其,在本申请的行文中,术语混浊介质应理解为指的是生物组织。异常区域应理解为指的是在其中的混浊介质以任何方式或形式有别于在周围区域的混浊介质的区域。特别是,在本申请的行文中,这种区域应理解为指的是包括肿瘤组织的区域。混浊介质由一组体积元素代表。这种体积元素也被称为三维像素(voxels)。体积元素的大小和形状对于所有的体积元素都是相同的。但是,体积元素也可能具有彼此不同的尺寸和形状。
从1999年7/8月出版的IEEE Journal of Selected Tops inQuantum Electronics的Vol.5,No.4的S.B.Colak等人的文章“Clinical Optical Tomography and NIR Spectroscopy for BreastCancer Detection”可获知这种方法和设备。已知的方法和设备可用于对生物组织内部成像,尤其可用于临床诊断中,用于胸内检查,对在人类或动物雌性体的胸部组织中出现的任何肿瘤作可视探测。根据已知方法,混浊介质由来自不同辐射位置的光连续照射。随后,在若干测量位置测量沿着通过混浊介质的从其辐射位置开始延伸的不同光路传输的光的强度。测量的强度用于重组混浊介质的图像。通过组织的光的衰减的空间分布在该图像中反映出来。光由组织衰减是由于组织散射和吸收光。
如果在异常区域中组织对光的衰减充分地与周围区域中组织对光的衰减偏离开,则在图像中能可视地对可能存在的异常区域探测。在使用不同波长的光的同时通过反复执行已知方法可形成多个独立图像。
已知方法和设备具有的缺点是这样得到的图像不能足以清楚到对异常区域进行准确的可视探测。
本发明的一个目的是提供一种方法,其产生一种可准确地对异常区域进行可视探测的混浊介质的图像。
该目的通过在导言中提出的那种方法实现,其特征在于该方法还至少包括第二测量步骤,其中由主要包括不等于第一波长的第二波长的辐射的光来照射混浊介质,测量主要包括第二波长的辐射的沿着通过混浊介质的多个光路传播的一部分光的强度。成像步骤包括下面的步骤:第一计算步骤,其中,根据测量步骤中对于体积元素测量的强度至少对第一参数赋值,所述第一参数代表混浊介质的属性;第二计算步骤,其中,根据至少第一参数的值对体积元素赋有效值;显示步骤,用于在混浊介质的图像中对于图像中的与体积元素之一对应的每一点显示该对应体积元素的有效值。
本发明一方面基于从多个测量得到的测量强度与不同波长的光的组合可产生代表组织的属性的参数的值这样的理解,其中这些属性(独立或组合地)使得正常组织与异常组织区分开来。这样得到的一些参数尤其可代表构成混浊介质的组分。例如,一个参数可代表水在混浊介质中的量的空间分布。
本发明还基于对如下事实的认识:在根据获得的参数值认为是异常的区域可在图像中标记出来。从而将用户,例如放射科医师的注意力吸引到图像中认为是异常的区域。
在混浊介质的图像中通过例如使用颜色来代表有效值。从而不同的颜色被赋予不同的有效值。如果有效值应连续变化,用另外的方法向图像添加平滑地改变的颜色。已知的色标,如彩虹色标可用于这一目的。
根据本发明的方法方案的特征在于它包括第三测量步骤,其中,由主要包括不等于第一波长或第二波长的第三波长的辐射的光来照射混浊介质,测量主要包括第三波长的辐射的沿着通过混浊介质的多个光路传播的一部分光的强度;在第一计算步骤中,根据测量步骤中对于体积元素测量的强度对第二参数赋值,所述第二参数代表混浊介质的与第一参数所代表的属性不同的属性;在第二计算步骤中,根据第一参数的值和第二参数的值对体积元素赋有效值。
这个方案中的各个测量步骤产生一组测量数据。用于两个参数的值从三组获得的测量数据确定。当在三个测量步骤中光的波长充分远离时,测量的数据组将不相关,并且使用已知的数学方法在所有情况下都可从测量数据组确定彼此独立的参数的值。
根据本发明的方法方案的特征在于第一波长具有830nm和900nm之间的值,第二波长具有750nm与830nm之间的值,第三波长具有655nm和750nm之间的值。具有上述波长的光可容易地由现代技术水平的半导体激光器产生。
当在三个测量步骤中使用的光的波长位于上述范围时,在组织中尤其是血红蛋白(Hb)和氧化血红蛋白(HbO2)的空间分布可适当地近似。这是由于生物组织的散射和吸收产生的光的衰减在这些区域很低,并且由于对于血红蛋白(Hb)而言,光的衰减作为光的波长的函数与氧化血红蛋白(HbO2)的不同。
特别吸引人的结果在第一波长具有大约867nm的值、第二波长具有大约780nm的值、第三波长具有大约715nm或者作为替代大约682nm的值时获得。
如果需要准确确定组织中脂肪对水的比率,则可增加第四测量步骤,其中使用的光包括波长为1040nm和1070nm之间的辐射,或者作为替代包括750nm和770nm之间波长的辐射。对脂肪与水的比率的更准确的确定也使得可更准确地确定组织中血红蛋白(Hb)和氧化血红蛋白(HbO2)的空间分布。
根据本发明的方法的一个有利方案的特征在于第一参数代表混浊介质中与体积元素之一对应的体积中的血液的量,第二参数代表混浊介质中与体积元素之一对应体积中的血液的脱氧率。
在本发明的行文中,血量应理解为指的是血红蛋白(Hb)和氧化血红蛋白(HbO2)的总量,而血液的脱氧率应理解为指的是比率(Hb)/(Hb+HbO2)。
申请人进行的临床测试证明肿瘤组织每单位体积中包含更多量的血液,而且,肿瘤组织中血液的脱氧率比正常组织高。而且,随着肿瘤越恶化血液量和脱氧率增高。因此,体积中的血液量和脱氧率是适当的参数,独立或结合来用于本发明的方法中,以在有肿瘤时在图像中对肿瘤探测。
根据本发明的方法的另一有利方案的特征在于参数限定参数空间,在第二计算步骤中,根据赋予相关体积元素的参数的值来在参数空间中对体积元素赋予位置,赋予体积元素的有效值由参数空间中该体积元素的位置确定。
根据参数数目,这种参数空间可以是一维(线)、二维(面)、三维(空间或体积)或更多维的。那么参数代表参数空间坐标系系统。体积元素在所述参数空间被赋予位置,所述位置的坐标由与参数相关的值构成。例如,在二维参数空间中,带有赋予第一参数的值X和赋予第二参数的值Y的体积元素将在参数空间中被给出一个位置,该位置具有坐标(X,Y)。
这种参数空间在用于执行该方法的设备中虚拟存在。而且,参数空间也向用户展示。例如,二维参数空间可在显示屏上显示。
根据本发明的方法的另一有利方案的特征在于参数空间被细分为至少第一子空间和第二子空间,当体积元素占据第一子空间中的位置时,其被赋予第一有效值,当体积元素占据第二子空间中的位置时,其被赋予第二有效值。
参数空间可人为地或自动地被细分为子空间。例如,用户可人工地在显示在显示屏上的二维参数空间中指示子空间;位于所述子空间中的体积元素被赋予代表异常区域的有效值。用户例如可通过比较跨参数空间的体积元素的位置分布与参考分布来选择这个子空间。从而在检查女性身体的胸部组织期间,在一个胸部发现的体积元素跨参数空间的位置分布可与另一胸部发现的分布相比。
参数空间自动细分为子空间可通过例如将参数空间细分为预定的子空间来执行。执行若干次参考测量,从而使得提前确定子空间应如何跨参数空间分布以及应对每个子空间中的体积元素赋予何有效值。或者,可通过已知的分段方法来实现自动细分。在这种情况下,例如应位于第一子空间中的体积元素的数目的百分比和应位于第二子空间中的体积元素的数目的百分比可提前指示出来。
尽管在上述说明中参数空间被细分为两个子空间,但显然对于熟悉本领域的技术人员而言可细分成更多的数目的子空间。对于每个子空间,对位于每个子空间中的体积元素赋予不同的有效值。也可对于多个子空间对于位于这种子空间中的体积元素赋予相同的有效值。
根据本发明的方法的方案的特征在于有效值至少由赋予体积元素的位置和参数空间中选择的位置之间的距离确定。
例如,根据欧几里德(Euclid)几何原理的两个点之间的距离可用作距离的量度。例如,位于二维参数空间中的两个点P1(x1,y1)和P2(x2,y2)之间的距离等于
,其中x1表示对于参数X的第一值、x2表示对于参数X的第二值、y1表示对于参数Y的第一值、y2表示对于参数Y的第二值。
例如,参数空间中体积元素的位置分布的重心被选择为选择位置。赋予体积元素的有效值可以是根据距离可连续地改变值。但是,距离还可被细分为若干离散地节距,有效值与每个节距相关。
根据本发明的方法的又一方案的特征在于有效值尤其由赋予体积元素的位置相对选择的位置的情况来确定。
除离开选择位置的距离外,赋予体积元素的位置相对选择位置的情况在赋予体积元素的有效值的方面也有决定作用。例如,带有高参数值的体积元素被赋予一有效值,而不是赋予位于离开选择点相同距离的体积元素但具有低的参数值的有效值。
本发明的目的还通过导言中提出的那种方法来实现,其特征在于:该方法还至少包括第二测量步骤,其中由主要包括不等于第一波长的第二波长的辐射的光来照射混浊介质,测量主要包括第二波长的辐射的沿着通过混浊介质的多个光路传播的一部分光的强度;成像步骤包括下面的步骤:第一计算步骤,其中,根据测量步骤中对于体积元素测量的强度至少对第一参数赋值,所述第一参数代表混浊介质的属性;重组步骤,用于至少从赋予体积元素的第一参数的值来重组混浊介质的图像。
本发明一方面基于对如下事实的认识:从多个测量得到的测量强度与不同波长的光的组合可产生代表组织的属性的参数的值,其中该属性(独立或组合地)使得正常组织与异常组织区分开来。这样得到的参数尤其可代表构成混浊介质的组分。
这种参数的成像产生一种图像,其中可比组织对光衰减的图像中更准确地探测异常区域。这在参数代表构成混浊介质的组分时尤其如此。例如,在血液的脱氧率(Hb)/(Hb+HbO2)的空间分布图像中,在雌性胸部组织中,带有肿瘤组织的区域是在视觉上可观察到的,因为一般讲,肿瘤组织的脱氧率高于正常组织。类似地,在每一单位体积的血液量(Hb+HbO2)的图像中,带有肿瘤组织的区域是在视觉上可观察到的,因为肿瘤组织中血液循环通常多余正常组织。
通过应用重组方法可根据参数值重组图像,该方法也用于根据测量的强度来重组光由组织衰减的图像。这种重组方法的例子是根据Hoogenraad的3维代数重组方法和根据Watchers的3维共轭梯度方法。
注意到描述的方法不仅适合于为检测异常区域来对生物组织成像,也适合于对生物组织的功能性属性成像。例如,图像可由婴儿的脑中的血液循环和血液氧含量组成。
本发明又一目的是提供一种用于产生混浊介质的图像的设备,其中可视地准确探测异常区域。
这个目的通过开头一段中提出的那种设备实现,其特征在于它还包括由主要包括不等于第一波长的第二波长的辐射的光来照射混浊介质的装置;测量主要包括第二波长的辐射的沿着通过混浊介质的多个光路传播的一部分光的强度的装置;根据对于体积元素测量的强度至少对第一参数赋值的装置,所述第一参数代表混浊介质的属性;根据至少第一参数的值对体积元素赋予有效值的装置;用于在混浊介质的图像中对于图像中的与体积元素之一对应的每一个点显示该对应体积元素的有效值的装置。
这个目的通过开头一段中提出的那种设备实现,其特征在于它还包括由主要包括不等于第一波长的第二波长的辐射的光来照射混浊介质的装置;测量主要包括第二波长的辐射的沿着通过混浊介质的多个光路传播的一部分光的强度的装置;根据对于体积元素测量的强度至少对第一参数赋值的装置,所述第一参数代表混浊介质的属性;用于至少根据赋予体积元素的第一参数的值来重组混浊介质的图像的装置。
前面段落提到的设备例如是以医学诊断工作站形式实施的。
本发明还涉及一种用于执行根据本发明的方法的计算机程序。计算机程序可包括用于执行根据本发明的方法的所有方案的指令,或者包括执行方法的方案子集的指令。计算机程序的指令例如装载到医学诊断工作站的存储器中,并且随后由这个工作站的处理器执行。根据本发明的再一目的,计算机程序被记录在记录载体上,例如CD-ROM盘或磁带。计算机程序的指令从记录载体传送到例如医学诊断工作站的存储器中。或者计算机程序的指令从记录载体传送到PC或任何其他标准计算机系统的存储器中。还可以是将计算机程序通过例如“全球网”的网络装载到存储器中。
本发明的上述和其他的更具体方面将参考附图在后面以举例方式具体说明。其中:
图1图示根据本发明的第一方法;
图2表示二维参数空间的例子;
图3图示根据本发明的第二方法。
图1图示根据本发明的第一方法。在第一测量步骤1期间,由主要包括867nm的波长的辐射的光来照射混浊介质,测量主要包括第一波长的辐射的沿着通过混浊介质的多个光路传播的一部分光的强度。这种测量步骤在欧洲专利申请98925884.3(PHN16.442)中更具体说明。接着,由主要包括780nm的波长的辐射的光来执行类似的第二测量步骤2,再接着,使用主要包括715nm的波长的辐射的光来执行第三测量步骤3。
在成像步骤4中使用重组算法5,根据第三测量步骤3中测量的强度来重组混浊介质的图像9。尽管在本例中根据第三测量步骤3中测量的强度来重组图像,但也可将第一测量步骤1或第二测量步骤2中测量的强度来用于这一目的。
在第一计算步骤6中,对于每个体积元素,根据三个测量步骤测量的三组强度来确定两个参数。这些参数代表检查的组织的局部特性,在本例中是体积中的血液量和体积中的血液的脱氧率。在本行文中,体积中的血液量应理解为指的是在体积中存在的血红蛋白(Hb)和氧化血红蛋白(HbO2)的总和,术语体积中的血液的脱氧率应理解为指的是体积中存在的血红蛋白与体积中存在的血液量的比率(Hb)/(Hb+HbO2)。这些参数尤其适合于对在胸部组织中肿瘤组织探测,因为肿瘤通常以更高的血液循环和更高的氧消耗为特征,导致两个参数上升。
在第二计算步骤7中,将体积元素映射在二维参数空间20中,这个空间的例子在图2中表示。这个参数空间20由在第一计算步骤6中确定的两个参数限定。在本例中,体积中的血液量由水平X轴22代表、体积中的血液脱氧率由垂直的Y轴21表示。这个参数空间中的体积元素被赋予一个位置23,该位置取决于相关体积元素的两个参数的值。这样,体积中具有值X1的血液量并且体积中具有值Y1的血液脱氧率的体积元素在参数空间中将被赋予具有坐标(X1,Y1)的位置。
随后,每个体积元素被赋予有效值,其取决于体积元素在参数空间20中的位置23。在本例中,如果认为体积元素位于异常区域则对其赋予有效值“1”,而如果认为它位于异常区域外部的区域中,则对其赋予有效值“0”。图2A表示包含两个高值参数并且示出上述的肿瘤组织的子区域24。在子区域24中具有位置25的体积元素被赋予有效值“1”。其他体积元素被赋予有效值“0”。尽管在本例中仅使用2个不同的有效值,但当参数空间被细分为多个子区域时可使用多个离散的有效值。还可使用连续变化的有效值。有效值代表例如概率值。
随后,赋予混浊介质的图像9中的体积元素的有效值在显示步骤8中展示。图像9中的对应于体积元素的点如果该相关体积元素被赋予有效值“1”则被赋予例如红色,如果赋予相关体积元素的有效值为“0”,则其被赋予灰色。这样,图像9由包含灰色值和带有红色值的区域10的混浊介质构成。区域10对应于混浊介质中的异常区域。对于熟悉本领域的技术人员来说是明显的,也可将其他的颜色和其他方法用于标记图像中的区域。
图2表示二维参数空间20的例子。参数空间20由两个参数限定。第一参数由水平X轴22表示,第二参数由垂直Y轴21表示。这个参数空间中的体积元素被赋予位置23,该位置取决于相关体积元素的两个参数的值。这样,具有关于第一参数值X1和关于第二参数值Y1的体积元素在参数空间中将被赋予具有坐标(X1,Y1)的位置。
图2A中,参数空间通过人工形成子区域24而被细分为2个子区域。位于子区域24的体积元素25被赋予与赋予位于子区域24外侧的体积元素的有效值偏离开的有效值。在胸部组织成像期间,子区域24的位置和形状可通过比较从在一个胸部上执行的测量中得到的参数空间与在另一胸部上执行的测量而得到的参数空间而选择。
图2B表示对体积元素赋予有效值的另一方法。现在根据参数空间20中体积元素的位置30与选择位置26之间的距离27来对体积元素赋予有效值。在本例中,选择位置26是参数空间20中体积元素的位置分布的几何重心。要被赋予的有效值可以是距离27的连续函数。距离27可被细分为若干离散节距,每个节距具有相关的有效值。
尽管第一体积元素的位置30和第二体积元素的位置31位于离开选择位置26相同的距离处,但仍可对两个体积元素赋予不同的有效值。体积元素的位置相对选择位置26的取向也是重要的。例如,带有位于参数空间20的右上象限28中的位置30的体积元素被赋予一个值,该值不同于赋予带有位于参数空间20的左下象限29中的位置3 1的体积元素的有效值。这在例如,如果体积中的血液量和体积中的血液的脱氧率被选择为参数时有意义。如已经说明的那样,两个参数值增大表明是肿瘤组织。对于位于右上象限18中的体积元素来说情况是这样,因此,这些体积元素应被赋予高有效值。但是,一个或两个参数的值降低没有意义。因此,在这两个参数基础上,位于其他3个象限之一的体积元素不应被赋予高有效值。
图3图示根据本发明的第二方法。在第一测量步骤期间,由主要包括867nm的波长的辐射的光来照射混浊介质,测量主要包括第一波长的辐射的沿着通过混浊介质的多个光路传播的一部分光的强度。接着,由主要包括780nm的波长的辐射的光来执行类似的第二测量步骤2,再接着,使用主要包括715nm的波长的辐射的光来执行第三测量步骤3。
在第一计算步骤6中,对于每个体积元素,根据三个测量步骤测量的三组强度来确定两个参数。接着,在重组步骤13中,体积元素的第一参数的值在图像11中成像,第二参数的值在图像12中成像。图像11中的点和图像12中的点对应于体积元素。这产生跨混浊介质的参数空间分布的图像。
当参数适当选择时,清楚地在图像中看到异常区域。例如,当体积中的血液量在图像11中成像,体积中的脱氧率在图像12中成像时,可看到存在于胸部组织中的带有肿瘤组织的区域,因为在肿瘤组织中两个参数值增加。
在本例中两个参数分别成像。但是,也可以仅对两个参数之一成像。而且如果在第一计算步骤6中确定两个以上的参数,还可以分别地对2个以上的参数成像。另外,可组合多个参数的成像来形成一个公共图像。
Claims (13)
1.一种在由一组体积元素所代表的混浊介质中探测异常区域的方法,该方法包括:
·第一测量步骤(1),其中,由主要包括第一波长的辐射的光来照射混浊介质,以及其中测量主要包括第一波长的辐射的沿着通过混浊介质的多个光路传播的一部分光的强度;
·成像步骤(4),用于从测量的强度来重组混浊介质的图像(9),其特征在于:
·该方法还至少包括第二测量步骤(2),其中由主要包括不等于第一波长的第二波长的辐射的光来照射混浊介质,
以及其中测量主要包括第二波长的辐射的沿着通过混浊介质的多个光路传播的一部分光的强度;
·成像步骤(4)包括下面的步骤:
-第一计算步骤(6),其中,根据测量步骤中对于体积元素测量的强度至少对第一参数赋值,所述第一参数代表混浊介质的属性;
-第二计算步骤(7),其中,根据至少第一参数的值对体积元素赋有效值;
-显示步骤(8),用于在混浊介质的图像(9)中对于图像中的与体积元素之一对应的一点显示该对应体积元素的有效值(10)。
2.根据权利要求1所述的方法,其特征在于:该方法包括第三测量步骤(3),其中,由主要包括不等于第一波长或第二波长的第三波长的辐射的光来照射混浊介质,测量主要包括第三波长的辐射的沿着通过混浊介质的多个光路传播的一部分光的强度;
在第一计算步骤(6)中,根据测量步骤中对于体积元素测量的强度对第二参数赋值,所述第二参数代表混浊介质的与第一参数所代表的属性不同的属性;
在第二计算步骤(7)中,根据第一参数的值和第二参数的值对体积元素赋有效值。
3.根据权利要求2所述的方法,其特征在于:第一波长具有830nm和900nm之间的值,第二波长具有750nm与830nm之间的值,第三波长具有655nm和750nm之间的值。
4.根据权利要求3所述的方法,其特征在于:第一参数代表混浊介质中与体积元素之一对应的体积中的血液的量,第二参数代表混浊介质中与体积元素之一对应体积中的血液的脱氧率。
5.根据权利要求1所述的方法,其特征在于:参数限定参数空间(20),在第二计算步骤(7)中,根据赋予相关体积元素的参数的值来在参数空间(20)中对体积元素赋予位置(23),赋予体积元素的有效值由参数空间中该体积元素的位置确定。
6.根据权利要求5所述的方法,其特征在于:参数空间(20)被细分为至少第一子空间(24)和第二子空间,当体积元素占据第一子空间中的位置时,其被赋予第一有效值,当体积元素占据第二子空间中的位置时,其被赋予第二有效值。
7.根据权利要求5所述的方法,其特征在于:有效值至少由赋予体积元素的位置(23)和参数空间中选择的位置(26)之间的距离确定。
8.根据权利要求7所述的方法,其特征在于:有效值尤其由赋予体积元素的位置相对选择的位置的情况来确定。
9.一种在由一组体积元素所代表的混浊介质中探测异常区域的方法,该方法包括:
第一测量步骤(1),其中,由主要包括第一波长的辐射的光来照射混浊介质,
以及其中测量主要包括第一波长的辐射的沿着通过混浊介质的多个光路传播的一部分光的强度;
·成像步骤(4),用于重组混浊介质的图像(9),其特征在于:
·该方法还至少包括第二测量步骤(2),其中由主要包括不等于第一波长的第二波长的辐射的光来照射混浊介质,
以及其中测量主要包括第二波长的辐射的沿着通过混浊介质的多个光路传播的一部分光的强度;
·成像步骤(4)包括下面的步骤:
-第一计算步骤(6),其中,根据测量步骤中对于体积元素测量的强度至少对第一参数赋值,所述第一参数代表混浊介质的属性;
-重组步骤(13),用于至少从赋予体积元素的第一参数的值来重组混浊介质的图像(11)。
10.一种在由一组体积元素所代表的混浊介质中对异常区域成像的设备,该设备包括:
由主要包括第一波长的辐射的光来照射混浊介质的装置;
测量主要包括第一波长的辐射的沿着通过混浊介质的多个光路传播的一部分光的强度的装置;
用于从测量的强度重组混浊介质的图像的装置,其特征在于:
该设备还包括由主要包括不等于第一波长的第二波长的辐射的光来照射混浊介质的装置;
测量主要包括第二波长的辐射的沿着通过混浊介质的多个光路传播的一部分光的强度的装置;
根据对于体积元素测量的强度至少对第一参数赋值的装置,所述第一参数代表混浊介质的属性;
根据至少第一参数的值对体积元素赋予有效值的装置;
用于在混浊介质的图像中对于图像中的与体积元素之一对应的每一个点显示该对应体积元素的有效值的装置。
11.一种在由一组体积元素所代表的混浊介质中对异常区域成像的设备,该设备包括:
由主要包括第一波长的辐射的光来照射混浊介质的装置;
测量主要包括第一波长的辐射的沿着通过混浊介质的多个光路传播的一部分光的强度的装置;
用于重组混浊介质的图像的装置,其特征在于:
该设备还包括由主要包括不等于第一波长的第二波长的辐射的光来照射混浊介质的装置;
测量主要包括第二波长的辐射的沿着通过混浊介质的多个光路传播的一部分光的强度的装置;
根据对于体积元素测量的强度至少对第一参数赋值的装置,所述第一参数代表混浊介质的属性;
用于至少从赋予体积元素的第一参数的值来重组混浊介质的图像的装置。
12.一种计算机程序,包括用于执行权利要求1或9的方法的一组指令。
13.一种记录载体,包括根据权利要求12的计算机程序。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP00201164.1 | 2000-03-31 | ||
EP00201164 | 2000-03-31 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1380821A true CN1380821A (zh) | 2002-11-20 |
CN1262235C CN1262235C (zh) | 2006-07-05 |
Family
ID=8171282
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB018014089A Expired - Fee Related CN1262235C (zh) | 2000-03-31 | 2001-03-22 | 用于探测在混浊介质中异常区域的方法和设备 |
Country Status (7)
Country | Link |
---|---|
US (2) | US6718195B2 (zh) |
EP (1) | EP1272102B1 (zh) |
JP (1) | JP4954418B2 (zh) |
CN (1) | CN1262235C (zh) |
AT (1) | ATE422325T1 (zh) |
DE (1) | DE60137623D1 (zh) |
WO (1) | WO2001074238A1 (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100566661C (zh) * | 2005-03-01 | 2009-12-09 | 国立大学法人神户大学 | 图像诊断处理设备和图像诊断处理程序 |
CN101312683B (zh) * | 2005-11-23 | 2011-04-06 | 皇家飞利浦电子股份有限公司 | 用于对混浊介质内部进行成像的方法和设备 |
CN102292019A (zh) * | 2009-01-23 | 2011-12-21 | 皇家飞利浦电子股份有限公司 | 配置为至少部分地插入到混浊介质中的光学检查装置 |
CN101262812B (zh) * | 2005-09-13 | 2012-04-25 | 皇家飞利浦电子股份有限公司 | 光学层析成像装置 |
CN103460025A (zh) * | 2011-03-24 | 2013-12-18 | 鹿特丹伊拉斯谟大学医疗中心 | 用于确定混浊介质中的吸收系数的方法 |
Families Citing this family (33)
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 |
JP2007528763A (ja) | 2004-03-05 | 2007-10-18 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | インタラクティブコンピュータ支援診断方法及び装置 |
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 |
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 |
CN101309633A (zh) * | 2005-11-18 | 2008-11-19 | 皇家飞利浦电子股份有限公司 | 用于对混浊介质内部进行成像的设备 |
CN101312682B (zh) | 2005-11-23 | 2011-05-25 | 皇家飞利浦电子股份有限公司 | 用于对混浊介质进行光学成像的方法和装置 |
US20090069695A1 (en) * | 2006-03-17 | 2009-03-12 | Koninklijke Philips Electronics N.V. | Device for imaging a turbid medium |
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 |
WO2009067680A1 (en) | 2007-11-23 | 2009-05-28 | Mercury Computer Systems, Inc. | Automatic image segmentation methods and apparartus |
WO2009067675A1 (en) | 2007-11-23 | 2009-05-28 | Mercury Computer Systems, Inc. | Client-server visualization system with hybrid data processing |
US9904969B1 (en) | 2007-11-23 | 2018-02-27 | PME IP Pty Ltd | Multi-user multi-GPU render server apparatus and methods |
US10311541B2 (en) | 2007-11-23 | 2019-06-04 | PME IP Pty Ltd | Multi-user multi-GPU render server apparatus and methods |
US8319781B2 (en) | 2007-11-23 | 2012-11-27 | Pme Ip Australia Pty Ltd | Multi-user multi-GPU render server apparatus and methods |
US9468379B2 (en) * | 2010-04-21 | 2016-10-18 | Koninklijke Philips N.V. | Determination of a lipid water ratio |
US11183292B2 (en) | 2013-03-15 | 2021-11-23 | PME IP Pty Ltd | Method and system for rule-based anonymized display and data export |
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 |
US10540803B2 (en) | 2013-03-15 | 2020-01-21 | PME IP Pty Ltd | Method and system for rule-based display of sets of images |
US8976190B1 (en) | 2013-03-15 | 2015-03-10 | Pme Ip Australia Pty Ltd | Method and system for rule based display of sets of images |
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 |
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 |
US9984478B2 (en) | 2015-07-28 | 2018-05-29 | PME IP Pty Ltd | Apparatus and method for visualizing digital breast tomosynthesis and other volumetric images |
US11599672B2 (en) | 2015-07-31 | 2023-03-07 | PME IP Pty Ltd | Method and apparatus for anonymized display and data export |
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 |
Family Cites Families (14)
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 | 運動モニタ装置 |
JP2795197B2 (ja) * | 1994-04-30 | 1998-09-10 | 株式会社島津製作所 | 光散乱・吸収体の光学的測定装置 |
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 |
US6124597A (en) * | 1997-07-07 | 2000-09-26 | Cedars-Sinai Medical Center | Method and devices for laser induced fluorescence attenuation spectroscopy |
EP0925017A1 (en) * | 1997-07-15 | 1999-06-30 | Koninklijke Philips Electronics N.V. | Method of localizing an object in a turbid medium |
US6091984A (en) * | 1997-10-10 | 2000-07-18 | Massachusetts Institute Of Technology | Measuring tissue morphology |
DE69822439T2 (de) * | 1997-11-22 | 2005-01-20 | Koninklijke Philips Electronics N.V. | Verfahren zur auffindung eines körpers in einem trüben medium |
US6327489B1 (en) * | 1998-06-25 | 2001-12-04 | U.S. Philips Corporation | Method of localizing an object in a turbid medium |
CA2343401C (en) * | 1998-09-11 | 2009-01-27 | Spectrx, Inc. | Multi-modal optical tissue diagnostic system |
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 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-22 DE DE60137623T patent/DE60137623D1/de not_active Expired - Lifetime
- 2001-03-22 AT AT01929448T patent/ATE422325T1/de not_active IP Right Cessation
- 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-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
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100566661C (zh) * | 2005-03-01 | 2009-12-09 | 国立大学法人神户大学 | 图像诊断处理设备和图像诊断处理程序 |
CN101441767B (zh) * | 2005-03-01 | 2011-11-09 | 国立大学法人神户大学 | 图像诊断处理设备 |
CN101262812B (zh) * | 2005-09-13 | 2012-04-25 | 皇家飞利浦电子股份有限公司 | 光学层析成像装置 |
CN101312683B (zh) * | 2005-11-23 | 2011-04-06 | 皇家飞利浦电子股份有限公司 | 用于对混浊介质内部进行成像的方法和设备 |
CN102292019A (zh) * | 2009-01-23 | 2011-12-21 | 皇家飞利浦电子股份有限公司 | 配置为至少部分地插入到混浊介质中的光学检查装置 |
CN102292019B (zh) * | 2009-01-23 | 2014-03-26 | 皇家飞利浦电子股份有限公司 | 配置为至少部分地插入到混浊介质中的光学检查装置 |
CN103460025A (zh) * | 2011-03-24 | 2013-12-18 | 鹿特丹伊拉斯谟大学医疗中心 | 用于确定混浊介质中的吸收系数的方法 |
CN103460025B (zh) * | 2011-03-24 | 2017-03-29 | QuaSpec有限公司 | 用于确定混浊介质中的吸收系数的方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
WO2001074238A1 (en) | 2001-10-11 |
US6922582B2 (en) | 2005-07-26 |
US20010026848A1 (en) | 2001-10-04 |
JP2003528678A (ja) | 2003-09-30 |
ATE422325T1 (de) | 2009-02-15 |
US20040147039A1 (en) | 2004-07-29 |
EP1272102B1 (en) | 2009-02-11 |
CN1262235C (zh) | 2006-07-05 |
DE60137623D1 (de) | 2009-03-26 |
US6718195B2 (en) | 2004-04-06 |
EP1272102A1 (en) | 2003-01-08 |
JP4954418B2 (ja) | 2012-06-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1262235C (zh) | 用于探测在混浊介质中异常区域的方法和设备 | |
US7853064B2 (en) | Method for the processing of radiological images in tomosynthesis for a detection of radiological signs | |
CN1294520C (zh) | 自动检测具有与疾病相关的核密度特征和细胞质密度特征的细胞的光学投影成像方法 | |
US10130318B2 (en) | Integrated microtomography and optical imaging systems | |
CN1452089A (zh) | Pet肿瘤图像的半自动分割算法 | |
KR20190101905A (ko) | 양전자방출 단층촬영 시스템 및 그것을 이용한 영상 재구성 방법 | |
CN109975335A (zh) | 一种ct检测方法及装置 | |
CN1169665A (zh) | 用于生物组织体内测量的光耦合器 | |
US11127174B2 (en) | Medical imaging system with a fixed array of x-ray detectors and a fixed array of x-ray emitters for producing a digital 3-dimensional image | |
CN1861002A (zh) | 监视肿瘤负荷的方法和系统 | |
US20180042483A1 (en) | Systems and methods for hyperspectral imaging | |
CN1653481A (zh) | 应用瞬时特征标记的发射计算机辅助断层成像装置和方法 | |
CN109801254A (zh) | 医学成像中的传递函数确定 | |
CN109872312A (zh) | 医学图像分割方法、装置、系统及图像分割方法 | |
CN106691487A (zh) | 成像方法和成像系统 | |
CN107111883A (zh) | 用于图像数据的纹理分析图 | |
CN1403057A (zh) | 利用三维欧氏距离变换实现ct图像中软组织显示的方法 | |
Hosseini et al. | The impact of preprocessing on the PET-CT radiomics features in non-small cell lung cancer | |
Chang et al. | Optical imaging of anatomical maps derived from magnetic resonance images using time-independent optical sources | |
CN103261878B (zh) | 用于使用x射线分析在对象中的受关注区域的方法和设备 | |
CN103908278A (zh) | 图像重建方法及装置、医疗成像系统 | |
CN105678711B (zh) | 一种基于图像分割的衰减校正方法 | |
CN113180601B (zh) | 一种动态光学乳腺成像的肿瘤特征分析装置及相关组件 | |
CN111260636B (zh) | 模型训练方法及设备、图像处理方法及设备以及介质 | |
Vanegas et al. | MOCA: a systematic toolbox for designing and assessing modular functional near-infrared brain imaging probes |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20060705 Termination date: 20190322 |
|
CF01 | Termination of patent right due to non-payment of annual fee |