CN1679053A - 用于处理通过包含多个光纤的波导获取的图像的方法 - Google Patents
用于处理通过包含多个光纤的波导获取的图像的方法 Download PDFInfo
- Publication number
- CN1679053A CN1679053A CNA038210312A CN03821031A CN1679053A CN 1679053 A CN1679053 A CN 1679053A CN A038210312 A CNA038210312 A CN A038210312A CN 03821031 A CN03821031 A CN 03821031A CN 1679053 A CN1679053 A CN 1679053A
- Authority
- CN
- China
- Prior art keywords
- image
- optical fiber
- obtains
- flux
- zone
- 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
- 239000013307 optical fiber Substances 0.000 title claims abstract description 132
- 238000000034 method Methods 0.000 title claims abstract description 67
- 238000012545 processing Methods 0.000 title abstract description 8
- 230000004907 flux Effects 0.000 claims abstract description 57
- 238000002347 injection Methods 0.000 claims abstract description 12
- 239000007924 injection Substances 0.000 claims abstract description 12
- 239000000835 fiber Substances 0.000 claims abstract description 6
- 238000005070 sampling Methods 0.000 claims abstract description 4
- 238000001914 filtration Methods 0.000 claims description 24
- 238000012937 correction Methods 0.000 claims description 15
- 239000000203 mixture Substances 0.000 claims description 13
- 238000001514 detection method Methods 0.000 claims description 12
- 230000003071 parasitic effect Effects 0.000 claims description 11
- 238000007476 Maximum Likelihood Methods 0.000 claims description 7
- 230000011218 segmentation Effects 0.000 claims description 5
- 102000002322 Egg Proteins Human genes 0.000 claims description 3
- 108010000912 Egg Proteins Proteins 0.000 claims description 3
- 210000004681 ovum Anatomy 0.000 claims description 3
- 238000003672 processing method Methods 0.000 claims description 2
- 238000011002 quantification Methods 0.000 claims description 2
- 230000024241 parasitism Effects 0.000 claims 4
- 238000010586 diagram Methods 0.000 description 8
- 230000004927 fusion Effects 0.000 description 7
- 238000005259 measurement Methods 0.000 description 4
- 230000008030 elimination Effects 0.000 description 3
- 238000003379 elimination reaction Methods 0.000 description 3
- 238000007689 inspection Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 238000012952 Resampling Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000002513 implantation Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000001915 proofreading effect Effects 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/90—Dynamic range modification of images or parts thereof
- G06T5/94—Dynamic range modification of images or parts thereof based on local image properties, e.g. for local contrast enhancement
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10056—Microscopic image
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
- Mechanical Coupling Of Light Guides (AREA)
- Image Processing (AREA)
Abstract
本发明涉及一种用于处理通过包含多个光纤的波导获取的图像的方法。本发明的特征在于,对于每个光纤而言,其包含在所获取图像上隔离出对应于所述光纤的区域,单独对每个光纤进行局部处理从而校正在每个光纤中探测到的光子通量,然后通过消除由光纤导致的图案重新构建所获取图像。该方法还包含采样过程,用于为每个光纤并根据采样图像获得可用于重新构建所获取图像的采样注入率。本发明还包含预先的步骤,该步骤包含通过目标探测光纤。
Description
本发明涉及一种用于处理通过包含多个光纤的波导获取的图像的方法。如果将本发明应用于医学成像领域,那么是非常有优势的。但是,本发明具有更广的适用范围,这是因为它可以应用于任何领域,只要在这一领域当中是通过包含多个光纤的波导来实现成像的。该图像波导使得可以获得图像。这样一种装置使得可以进行激光扫描,并且使得光源和接收机可以远离所要观察的目标精确地移动。例如,在其中获得了共焦图像的激光扫描系统中,该图像波导包含几千个光纤,并且该图像波导的空间结构在进光口和出光口处是一样的。由于在光纤之间的信息损失,因此可以将通过该波导观察物体与通过栅格进行观察相比。因此,由于光纤的存在导致显示受到阻碍:光纤的图案出现在所获取图像上。这使得必须进行特定的处理,从而可以消除这种图案并且提高图像的可读性。
在目前的系统中,这种特定的处理仅限于对所获取图像进行线性滤波。
本发明的目的在于提出一种使得通过多光纤波导获取的图像具有可读性的新方法。
本发明的另一目的在于在对所获取图像进行处理时,考虑到由获取装置导致的寄生效应的因素。
利用用于处理通过由多个光纤构成的波导获取的图像的新方法,可以实现上述目的中的至少一个。根据本发明,对于每个光纤而言,在所获取的图像上隔离出对应于该光纤的一个区域,对每个区域都单独进行局部处理,然后重新构建该所获取图像从而消除由光纤导致的图案。
通过根据本发明的方法,可以将在图像上隔离出光纤分解为在对应于每个光纤的区域的图像上进行隔离。当隔离出就表面(有影响的区域)而言代表进入到光纤中的大多数注入量的象素时,就可以对每个光纤进行局部处理。波导的低串扰和注入的设置使得可以保证每个光纤的信息内容并不依赖于相邻的光纤,而只是依赖于所观察的物体的空间相关性。设置实现图像获取的装置,从而对于每个光纤来说具有足够的象素:因此可以准确地估计该光纤所探测到的和分布在代表该光纤的象素上的信息。
因此,不是如现有技术中那样通过简单的线性滤波从而实现对图像进行全局处理,在根据本发明的方法中,在所获取图像上隔离出每个光纤并且对每个光纤探测到的信息都进行处理。有优势的是,监视实现图像获取的装置,从而确保根据本发明的方法的最小效率条件。为此,可以改变采样率、对光纤进行注入的质量以及探测链的设置,从而确保“卵盒”类型的轮廓,特别是在参考图像上。
有利的是,当在图像上隔离出光纤时和当对每个隔离出的区域进行处理时,可以设想出多种应用,例如:
-重新构建没有光纤图案的图像:光纤阻碍了可读性和对图像的后续处理;
-监视波导表面的粗糙程度,这种粗糙会干扰注入现象并且使得注入现象在空间上发生变化;
-重新设置图像或者使图像稳定;光纤的图案防止了图像相对于彼此之间的重新设定,而获知关于光纤的位置和所观察到的信息使得图像可以对准;
-超分辨率:在获取一系列图像时可以使用小的移动,从而以更小的空间周期对图像进行重新采样,并且因此获得更高的分辨率;
-图像的量化:可以更容易和准确地提取关于图像的信息,而没有光纤的图案,
-对获取装置的内部参数的时间控制:获知关于每个光纤的位置和它们的最佳注入值使得可以监视波导上的磨损,以及某些光电参数的变化。
在根据本发明的方法中,为了隔离出每个区域,可以将对应于光纤图案的掩模施加到所获取的图像上。可以在根据参考图像探测光纤的步骤期间获得该对应于代表每个光纤的相关成分的图像的掩模。
该参考图像是使得光纤彼此之间可以很好地区分开的一种图像。可以通过对反射境、均匀散射介质和均匀荧光介质的观察而得到该图像,也可以通过对光纤束内部特定的背向散射而得到该图像。但是,该图像也可以是所获取的图像。因此,当离开该探测步骤时,就可以获得代表每个光纤的相关成分(段)的图像。每个灰度级代表单个指数,该指数表示波导中的一个光纤。
根据本发明,探测光纤的步骤包含以下步骤:
-对参考图像进行预滤波,
-使用本文中称作LPE的“分水岭(watershed)”算法,按照区域进行分段,
-对具有异常大的表面的段进行校正,以及
-对具有异常小的表面的段进行校正。
这两个校正步骤是可以互换的,并且可以反复执行这两个步骤。
有利的是,该预滤波步骤可以包含形态打开步骤和跟在其后的图像反转步骤。使用该数字形态打开,以便寻求消除位于光纤上的寄生最大值。这是在按照区域进行分段期间执行的LPE“分水岭”算法的标准预处理。
在预滤波步骤的输出端处,可以获得这样一种图像,其中光纤经过滤波从而消除了它们的局部最大值,并且在光纤间的区域是平滑的。
该图像反转步骤之前可以有一个分级类型各向异性散射步骤。
根据本发明,该预滤波还可以包含这样一个步骤,在该步骤期间,执行对最近相邻者的内插,从而使图像的尺寸在垂直方向和水平方向变为原来的两倍。
通过该内插步骤,寻求模拟半径小于1的数学形态结构元素。图像变为原来的两倍,从而接下来的形态打开并不影响被隔离出来的最大值,而只是影响那些与8相关但与4无关的值(对角线相邻者)。这样做的好处在于选出了通过该打开消除的最大值。
而且,当存在多个获取图像时,预滤波还可以包含时间滤波的步骤。
有利的是,对每个区域的局部处理可以包含计算为所获取图像的每个区域(对应于给定的光纤)探测到的光子通量,以及校正每个如此算得的通量值的偏差。
优选的是,使用根据每个光纤的特定注入轮廓算得的最大似然估计量来对通量进行计算。更确切的是,可以将该最大似然估计量用于计算每个光纤中特定注入轮廓的幅度分布。该轮廓是代表注入率的曲线,其中注入率是关于光在光纤末端与横截面中心的距离的函数。通常用高斯滤波来模拟该轮廓。
根据本发明的一个实施例,还将掩模应用于代表寄生背景的图像,也计算为背景图像的每个区域探测到的光子通量,并且从所获取图像的每个区域的每个通量值中减去对应的背景图像的每个区域的通量值,而且对该减法的结果执行偏差校正。
背景图像可以是在获取装置的光学系统上的寄生反射,因而包含在图像波导的出射端处的那些反射,但是它也可以是获取装置的数字化链的偏差、电子噪声。如果该偏差在图像上占据了主导地位,那么不能通过简单地移走图像而获得背景,这是因为该偏差依赖于内容,并且因此不再是相同的。在这种情况下,使用柱状图的分位数来估计它。如将在下面看到的,该柱状图是实时测量期间所获取图像的柱状图和调整步骤期间调整图像的柱状图。
寄生背景可能来自背景图像或来自探测链的偏差。
根据本发明,偏差校正可以包含在空间上将光纤分成不同的块,估计每个块中的偏差值,对偏差值进行内插从而获得每个光纤的偏差值,以及对于每个区域,将在前一个步骤获得的通量值除以如此获得的对应偏差的值。
在局部处理之后,所获取图像的重新构建可以包含用于校准所获取图像的通量的校准步骤,以及马赛克重新构建步骤。也可以使用其它类型的重新构建,例如通过内插或以径向函数为基础。
而且,可以在不使用调整和校准步骤的情况下实施本发明。在这种情况下,参考图像可以是所获取图像。
为了进行校准并且对于所获取图像的每个区域而言,可以将局部处理之后获得的通量值除以在调整步骤之后获得的通量值。该除法操作使得可以补偿对某些光纤的错误注入。
根据本发明的有利特征,该调整步骤包含:
-通过将对应于光纤图案的掩模施加到调整图像从而隔离出该调整图像的每个区域,
-计算为该调整图像的每个区域探测到的光子通量,以及
-校正每个这样算得的通量值的偏差。
在完成调整时获得的图像可以作为用于所获取图像的标准,从而获得一种所获取图像,其中对于该所获取图像来说,已经通过相同的方式对波导的全部光纤进行注入。
优选的是,使用根据每个光纤的特定注入轮廓算得的最大似然估计量计算通量。而且,也将掩模施加到代表寄生背景的图像,可以计算为背景图像的每个区域探测到的光子通量,并且从调整图像的每个区域的每个通量值中减去对应的背景图像的每个区域的通量值,而且可以对该减法的结果执行偏差校正。
换句话说,在调整步骤所执行的操作与测量期间执行的操作相同,该测量期间即为处理实时获取的图像期间。除此之外,将调整图像用于调整,该调整图像大致显示了光纤的图案。该调整使得可以在参考图像上探测光纤之后生成图像,其中注入率将会作为实时测量期间的标准。在测量期间,还根据所获取图像来校正注入率,并且将观察到的通量作为标准图像的函数进行校正,从而重新构建没有光纤图案的所获取图像。
偏差意味着低频成分,并且该成分可能是由多种原因产生的。
可以在调整期间和实时测量期间执行偏差校正。在第一种情况下,之所以能执行偏差校正是因为例如对平面反射镜进行了调整,并且因为场曲率将会降低返回到边缘时的注入质量(这些边缘是位于焦点之外的)。在测量物体上或者在均匀散射介质中,在边缘上的注入仍不令人满意,并且这在偏差上得到了体现,而该偏差就其形式而言与第一种是类似的。该偏差还可能是由渐晕问题引起的。通常,该偏差几乎具有圆对称性。实际上可以通过将图像分成具有固定尺寸的N×N个块,然后估计每块上的偏差,来估计该偏差。为此,有必要考虑所观察物体的特性。当观察均匀物体时,可以通过取块的平均值或中值来获取该偏差。当有一个物体时,有必要知道该物体比图像的其余部分是暗一些还是亮一些。例如,该偏差可以是乘法性的,并且因此更可能采用平均值运算或是中值运算符(与用于加法性偏差的最大值和最小值相比较而言)。因此,可以获得N×N的图像测量,在内插之后将该图像测量用于确定每个光纤所观察到的偏差值。例如,可以使用双线性内插。
马赛克重新构建可以包含在所获取图像的每个区域的整个表面上分配在校准步骤之后获得的每个区域的通量值。然后,可以执行低通递归滤波,从而使重新构建的所获取图像变得平滑。
根据本发明的一个有利实施例,参考图像和调整图像是相同的。
如果参照对实施例的详细描述和附图,则可以清楚地了解本发明的其它优势和特点,而该实施例并不起到限制作用,在附图中:
图1是根据本发明的图像处理方法的全局示意图;
图2是详细描述根据本发明的调整方法的主要步骤的流程图;
图3是详细描述根据本发明的测量方法的主要步骤的流程图;以及
图4是详细描述根据本发明的最后步骤的流程图,其中考虑到了用于重新构建没有可见光纤的所获取图像的调整和测量方法。
在图1中,显示了根据本发明的方法的全局示意图。
可以将该方法分为四部分:与光纤探测相关的部分1,调整部分2,测量部分3以及重新构建部分4。部分1和部分2对应于图2中的框图,而部分3对应于图3的框图,以及图4对应于图4的框图。
在图1中,根据本发明,当希望获取多个图像时,首先要执行调整步骤。为此,考虑参考图像5,而该参考图像5受制于探测光纤6的步骤,从而获得代表每个光纤的相关成分的图像7。参考图像5是通过包含波导的系统获取的图像,其中该波导包含多个光纤,例如包含10000到30000个光纤。获得参考图像5,从而区分出光纤图案,即“卵盒”类型的图案:在图像轮廓上,光纤出现在被脊和谷包围的小山中。参考图像5经历了探测光纤的过程,从而获得了代表光纤图案的一种掩模。该掩模是代表每个光纤的相关成分的图像。每个灰度级代表单个指数,该指数表示波导中的光纤。
然后,将图像7用于图像获取系统的调整2。该调整的目的在于确定一个光纤接一个光纤光子注入率的图像。就目前情况而言,该调整步骤是必要的,这是因为光纤之间具有略微不同的物理特征。因此,就每个光纤传导相同的光子通量的能力而言存在一定的差异。
为了执行调整2,考虑通过在图像波导的光学系统之前放置反射镜而获得的图像8,下文中将该图像8称作反射镜图像。该图像也可以是均匀散射介质、均匀荧光介质或对光纤束内部特别的背向散射的图像。该图像8也可以与在步骤5中使用的图像相同,即参考图像。掩模7用于确定在步骤9期间由每个光纤观察到的反射镜图像8的光子通量。可选的是,可以为寄生背景图像10确定在步骤11期间由每个光纤观察到的光子通量。该图像10不仅可以对应于获取系统的光学器件上的寄生反射,还可以对应于由获取系统的数字化链导致的偏移和/或电子噪声。步骤11的通量计算还涉及掩模7,从而标识出对应于每个光纤的区域。在步骤12中,对于每个光纤,将反射镜图像8的光子值减去背景图像10的光子通量值。在步骤13中,对于每个光纤,估计在步骤12中获得的差值对应于每个光纤的标准注入率(步骤13)。
在步骤14中,对步骤13的图像进行偏差校正。因此,通过步骤14得到的图像是代表对应于光纤的每个区域的标准值和校正后光子通量值的图像。通过步骤14得到的图像用作由获取系统实时获取的多个图像的参考。该实时获取和处理的图像经历了部分3和4中显示的处理。
测量部分3接收所获取图像15,该图像典型地为测量物体的图像。如在调整部分2中实现的那样,在此处步骤18中也计算关于所获取图像15的每个光纤观察到的光子通量。为此,使用了掩模7从而在所获取图像15上标识出对应于每个光纤的区域。可选的是,如前面那样,考虑寄生背景图像16,该图像可以是真实的图像,即对应于所获取图像15的背景的图像,或者该图像可以是对应于获取系统的噪声的估计得到的图像。这个背景图像16还经历了使用掩膜7的步骤17,从而确定每个光纤观察到的光子通量。在步骤19中执行减法操作。该减法的结果是代表对应于给定光纤的每个区域的可用光于通量的图像。步骤21是可选的步骤,在此步骤期间对图像20实施偏差校正。
重新构建部分4一方面接收校正后的图像20,另一方面接收校正后的(消除偏差的)图像13,以执行校准操作22,从而将观察到的物体的通量(通过步骤21得到的成分)除以标准通量(通过步骤14得到的成分)。在步骤22中,也执行重新构建,从而获得没有可见光纤图案的重新构建的图像23。
在图2中,更详细地描述了探测光纤的方法1和调整方法2。探测光纤的操作6涉及四种操作:
-预滤波,
-对应于按照区域进行分段的LPE“分水岭”;
-对具有异常大表面的段进行校正;以及
-对具有异常小表面的段进行校正。
这两种校正操作是可以互换的,并且可以循环执行这些操作。
预滤波操作61在入口处接收参考图像5并且生成滤除了这些局部最大值并且在光纤间区域变得平滑的光纤图像。该预滤波包含形态打开操作,可选的是随后进行分级型各向异性散射,然后进行图像反转。当有几个相同的固定物体的图像时,可以执行图像的时间滤波。为了使图像的尺寸在垂直方向和水平方向变为原来的二倍,还可以执行向最近相邻者的内插。
由预滤波61生成的图像然后经历分水岭操作62,从而使得可以从受探测光纤的相关成分中获得图像。常规类型的分水岭操作将会使得可以对位于图像边缘处的段做标记,从而使它们从最后的结果中去除。然后,计算所获得的这些段的特性(平均尺寸、标准偏差、近似度)。通过操作62得到的图像将会经历两次连续的校正63和64。操作63是对具有异常大的表面的段进行校正。为此,选择出相对于平均尺寸过大并且相对于光纤的正常表面具有过多的相邻元素的段。通过在原始图像上或者在探测到的段内部的距离卡图像(在一个段中,在每个象素和段的边缘之间的距离)上的分水岭,来对这些段进行重新分段。然后,计算获得的段的特性(平均尺寸、尺寸的标准偏差)。
在步骤64中,对具有异常小的表面的段进行校正。还计算这些段的相邻图形,然后确定哪些光纤必须熔合以及哪些光纤仅仅为候选光纤。在每种情况下,可能进行熔合的光纤集合是与每个相邻光纤进行熔合的光纤集合。对于必须进行熔合的光纤而言,采用给出最小密度值的可能的熔合。对于其它仅仅作为候选的光纤,需要连续使用三个滤波器,从而消除会产生错误结果的熔合。第一个滤波器检查熔合后的尺寸是否是过大。第二个滤波器检查熔合后的密度是否超过最大值。最后一个滤波器检查熔合是否提高了密度。如果还留有几个熔合,那么就保留那些产生较好(最小的)密度结果的熔合。一旦已经对相邻图形执行了所有的熔合,那么就在出口处将结果反映到相关成分的图像上。然后,计算获得的段的特性(平均尺寸、尺寸的标准偏差)。第二个校正64使得可以生成相关成分的图像,其中该相关成分的图像用作调整2和测量3的掩模。
图2的右边的部分涉及如图1中所示的调整方法2,但是以一种简化的形式表示出来。可选的操作10和11并没有在图2中出现。
施加掩模步骤91包含在反射镜8上标识出对应于波导的每个光纤的区域或表面。相关成分的图像7起到掩模的作用。在步骤92中,为每个光纤计算来自观察的物体的通量。使用对每个光纤的特定注入轮廓确定的最大似然估计量计算该通量。在步骤14中,执行偏差校正,如将在图3中详细描述的那样。在步骤14的输出端处,在步骤24获得每个光纤的光子注入率的标准值。
图3显示了根据本发明的测量方法。该方法以实时的方式进行。在所获取图像15上,通过使用掩模7执行施加掩模操作181来标识出对应于每个光纤的区域。然后,在步骤182中,为每个光纤计算来自观察的物体的通量。如前所述,使用对每个光纤的特定注入轮廓确定的最大似然估计量执行计算。在寄生背景图像16上进行相同的施加掩模操作171和通量计算操作172。
在步骤19中,对于对应于光纤的每个区域而言,从所获取图像的通量(182)中减去背景通量(172)。
如果没有提供背景图像并且该背景存在,那么必须减去该背景,从而使用所获取图像15的柱状图的分位数来计算该测量系统的偏差和/或电子噪声。
然后在步骤21中,可以对通过减法19得到的图像执行偏差校正。在这种情况下,在步骤211中将对应于光纤的区域在空间上分成不同的块。在步骤212中,使用给定的运算符计算这些块的每个块中的偏差值。然后在步骤213中,对偏差值进行内插从而获得每个光纤的值。然后在步骤214中,将每个光纤观察到的通量值除以获得的偏差值。
在步骤21的最后,为每个光纤生成代表观察到的通量的图像25。
最后的重新构建步骤将标准注入率24的图像和观察到的通量25的图像考虑在内。校准操作的目的在于通过使所有光纤的注入率相等从而获得一种图像来补偿注入损失,其中以相同的方式注入该图像的所有光纤。为此,在步骤221中,将观察到的图像通量25除以注入率24的图像。
然后,通过在步骤222中在对应于每个光纤的整个表面上分配校正(除法)之后获得的值来执行马赛克重新构建。为了生成更规则的外表,可以执行例如高斯低通滤波223。
当然,本发明并不限于上述的实施例,并且可以在不偏离本发明的范围的情况下实现多种变化。
Claims (34)
1.用于处理通过包含多个光纤的波导获取的图像的方法,该方法的特征在于,对于每个光纤,在所获取图像上隔离出对应于该光纤的区域,对每个区域单独进行局部处理,然后重新构建所获取图像,从而消除由光纤导致的图案。
2.根据权利要求1所述的方法,特征在于,为了隔离出每个区域,向所获取图像施加对应于光纤图案的掩模。
3.根据权利要求2所述的方法,特征在于,该掩模是在根据参考图像探测光纤的步骤期间获得的,其中该掩模对应于代表每个光纤的相关成分的图像。
4.根据权利要求3所述的方法,特征在于,该探测光纤的步骤包含下面的步骤:
-对该参考图像进行预滤波,
-按照区域进行分段,
-对具有异常大的表面的段进行校正,以及
-对具有异常小的表面的段进行校正。
5.根据权利要求4所述的方法,特征在于,该两个校正步骤是可以互换的。
6.根据权利要求4或5所述的方法,特征在于,重复执行该两个校正步骤。
7.根据权利要求4到6中任意一项所述的方法,特征在于,该预滤波步骤包含形态打开步骤,以及而后的图像反转步骤。
8.根据权利要求7所述的方法,特征在于,该图像反转步骤之前是分级类型各向异性散射步骤。
9.根据权利要求4到8中任意一项所述的方法,特征在于,该预滤波还包含这样一个步骤,在该步骤期间执行向最近相邻者的内插,从而使图像的尺寸在垂直方向和水平方向都变为原来的二倍。
10.根据权利要求4到9中任意一项所述的方法,特征在于,当存在多个获取图像时,该预滤波还包含时间滤波步骤。
11.根据上述权利要求中的任意一项所述的方法,特征在于,对每个区域的局部处理包含计算为所获取图像的每个区域探测到的光子通量,以及校正每个由此算得的通量值的偏差。
12.根据权利要求11所述的方法,特征在于,使用对每个光纤的特定注入轮廓算得的最大似然估计量来计算该通量。
13.根据权利要求7所述的方法,特征在于,还在代表寄生背景的图像上施加掩模,另外还计算为背景图像的每个区域探测到的光子通量,并且从所获取图像的每个区域的每个通量值中减去对应的背景图像的每个区域的通量值,以及对该减法的结果进行偏差校正。
14.根据权利要求13所述的方法,特征在于,该寄生背景来自图像的背景。
15.根据权利要求13所述的方法,特征在于,该寄生背景来自对探测链的偏离的计算。
16.根据权利要求11到15中任意一项所述的方法,特征在于,该偏差校正包含在空间上将光纤分为不同的块,估计每个块中的偏差值,对偏差值进行内插从而获得每个光纤的偏差值,以及对于每个区域,将在前面的步骤中获得的通量值除以由此获得的对应的偏差值。
17.根据上述权利要求中的任意一项所述的方法,特征在于,对所获取图像的重新构建包含在局部处理之后的用于对所获取图像的通量进行校准的校准步骤,以及马赛克重新构建步骤。
18.根据权利要求17所述的方法,特征在于,对于该校准,并且对于所获取图像的每个区域,将局部处理之后获得的通量值除以在调整步骤之后获得的通量值。
19.根据权利要求18所述的方法,特征在于,该调整步骤包含:
-隔离出调整图像的每个区域,
-将对应于光纤图案的掩模施加到该调整图像,
-计算为该调整图像的每个区域探测到的光子通量,以及
-对每个如此算得的通量值的偏差进行校正。
20.根据权利要求19所述的方法,特征在于,使用对每个光纤的特定注入轮廓算得的最大似然估计量来计算该通量。
21.根据权利要求18或19所述的方法,特征在于,还向代表寄生背景的图像上施加掩模,另外还计算为背景图像的每个区域探测到的光子通量,并且从调整图像的每个区域的每个通量值中减去对应的背景图像的每个区域的通量值,以及对该减法的结果进行偏差校正。
22.根据权利要求21所述的方法,特征在于,该寄生背景来自图像的背景。
23.根据权利要求21所述的方法,特征在于,该寄生背景来自对偏离的计算和探测链的噪声。
24.根据权利要求17到23中任意一项所述的方法,特征在于,该马赛克重新构建包含在所获取图像的每个区域的整个表面上分布在校准步骤之后获得的每个区域的通量值。
25.根据权利要求24所述的方法,特征在于,进行低通滤波从而使得该重新构建的所获取图像更加规则。
26.根据权利要求3到25中任意一项所述的方法,特征在于,该参考图像是通过在波导的对面放置反射镜而获得的图像。
27.根据权利要求3到25中任意一项所述的方法,特征在于,该参考图像是通过均匀散射介质获得的图像。
28.根据权利要求3到25中任意一项所述的方法,特征在于,该参考图像是通过均匀荧光介质获得的图像。
29.根据权利要求3到25中任意一项所述的方法,特征在于,该参考图像是通过在组成波导的光纤束的内部的背向散射获得的图像。
30.根据权利要求3到25中任意一项所述的方法,特征在于,该参考图像是所获取图像。
31.根据权利要求19到30中任意一项所述的方法,特征在于,该参考图像与该调整图像是相同的。
32.用于使用包含多个光纤的波导获取图像并且实施根据上述权利要求中任意一项所述的方法的装置,特征在于,对于每个光纤而言,该装置包含:
-用于在所获取图像上隔离出对应于该光纤的区域的装置,
-用于单独对每个区域进行局部处理的装置,以及
-用于重新构建所获取图像,从而消除由光纤导致的图案的装置。
33.根据权利要求32所述的装置,特征在于,该装置包含用于改变采样率、注入到光纤的质量以及探测链的设置从而确保“卵盒”轮廓的装置。
34.将根据上述权利要求中任意一项所述的图像处理方法的应用于下面的领域之一:
-监视波导表面的粗糙度;
-重新设置图像,或者稳定图像;
-所获取图像的超分辨率;
-图像的量化;以及
-获取装置的内部参数的时间控制。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR02/09100 | 2002-07-18 | ||
FR0209100A FR2842628B1 (fr) | 2002-07-18 | 2002-07-18 | "procede de traitement d'une image acquise au moyen d'un guide compose d'une pluralite de fibres optiques" |
PCT/FR2003/002197 WO2004010377A1 (fr) | 2002-07-18 | 2003-07-11 | Procede de traitement d'une image acquise au moyen d'un guide compose d'une pluralite de fibres optiques |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1679053A true CN1679053A (zh) | 2005-10-05 |
CN1679053B CN1679053B (zh) | 2013-07-17 |
Family
ID=29797536
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN038210312A Expired - Lifetime CN1679053B (zh) | 2002-07-18 | 2003-07-11 | 用于处理通过包含多个光纤的波导获取的图像的方法 |
Country Status (11)
Country | Link |
---|---|
US (1) | US7903848B2 (zh) |
EP (1) | EP1523731B1 (zh) |
JP (1) | JP4485947B2 (zh) |
KR (1) | KR20050021493A (zh) |
CN (1) | CN1679053B (zh) |
AU (1) | AU2003269019B2 (zh) |
CA (1) | CA2491751C (zh) |
ES (1) | ES2687644T3 (zh) |
FR (1) | FR2842628B1 (zh) |
IL (1) | IL166152A (zh) |
WO (1) | WO2004010377A1 (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108882831A (zh) * | 2016-03-31 | 2018-11-23 | 国立大学法人东北大学 | 光学成像装置 |
WO2019076267A1 (zh) * | 2017-10-16 | 2019-04-25 | 苏州微景医学科技有限公司 | 光纤束图像分析方法和装置 |
Families Citing this family (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2871358B1 (fr) | 2004-06-14 | 2007-02-09 | Mauna Kea Technologies Soc Par | Procede et systeme d'imagerie microscopique de fluorescence fibree multimarquage |
FR2877103B1 (fr) | 2004-10-22 | 2012-02-10 | Mauna Kea Technologies | Systeme et procede d'imagerie microscopique multiphotonique fibre d'un echantillon |
US7307774B1 (en) | 2005-01-24 | 2007-12-11 | The Board Of Trustees Of The Leland Standford Junior University | Micro-optical analysis system and approach therefor |
US8346346B1 (en) | 2005-01-24 | 2013-01-01 | The Board Of Trustees Of The Leland Stanford Junior University | Optical analysis system and approach therefor |
US8788021B1 (en) | 2005-01-24 | 2014-07-22 | The Board Of Trustees Of The Leland Stanford Junior Univerity | Live being optical analysis system and approach |
KR100682978B1 (ko) * | 2005-08-19 | 2007-02-15 | 한국전자통신연구원 | 영상을 이용한 얼굴의 광학 특성 맵 추출 시스템 및 그방법 |
DE102006011707B4 (de) * | 2006-03-14 | 2010-11-18 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Verfahren und Vorrichtung zum Erzeugen einer strukturfreien fiberskopischen Aufnahme |
FR2901029B1 (fr) * | 2006-05-12 | 2012-12-21 | Mauna Kea Technologies | Dispositif et procede d'endoscopie pour une observation simultanee de plusieurs zones d'interet. |
TWI355172B (en) * | 2006-06-07 | 2011-12-21 | Qualcomm Inc | Propagating session state changes to network funct |
US8099156B1 (en) | 2006-09-15 | 2012-01-17 | The Board Of Trustees Of The Leland Stanford Junior University | Cochlear optical analysis system and approach therefor |
JP4892316B2 (ja) | 2006-11-06 | 2012-03-07 | 株式会社フジクラ | マルチコアファイバ |
US8068899B2 (en) | 2007-07-03 | 2011-11-29 | The Board Of Trustees Of The Leland Stanford Junior University | Method and system of using intrinsic-based photosensing with high-speed line scanning for characterization of biological thick tissue including muscle |
KR100869469B1 (ko) * | 2008-04-01 | 2008-11-19 | 한국과학기술원 | 섬유 이미지의 검출 성능을 향상시키기 위한 처리 방법 |
US8983581B2 (en) * | 2008-05-27 | 2015-03-17 | Massachusetts Institute Of Technology | System and method for large field of view, single cell analysis |
US8270762B1 (en) * | 2008-09-27 | 2012-09-18 | Schott Corporation | Automated light-line calibration system and method of calibrating an optical fiber light line assembly |
JP5629270B2 (ja) * | 2008-12-29 | 2014-11-19 | マウナ ケア テクノロジーズ | 画像処理方法および装置 |
US8600134B2 (en) | 2009-01-30 | 2013-12-03 | Mauna Kea Technologies | Method and system for processing images acquired in real time through a medical device |
US9155471B2 (en) * | 2009-05-27 | 2015-10-13 | Lumicell, Inc'. | Methods and systems for spatially identifying abnormal cells |
US8332332B2 (en) * | 2010-01-29 | 2012-12-11 | Xerox Corporation | Methods and apparatus for managing pre-paid printing system accounts |
JP5570373B2 (ja) * | 2010-09-29 | 2014-08-13 | 富士フイルム株式会社 | 内視鏡システム |
US9314304B2 (en) | 2010-12-08 | 2016-04-19 | Lumicell, Inc. | Methods and system for image guided cell ablation with microscopic resolution |
AU2012257496B2 (en) | 2011-05-16 | 2017-05-25 | Mauna Kea Technologies | Continuous and real-time calibration of fiber-based microscopic images |
EP2785250A4 (en) | 2011-11-28 | 2015-07-29 | Univ Leland Stanford Junior | SYSTEM AND METHOD USEFUL FOR SARCOMER IMAGING THROUGH OBJECTIVE-BASED MICROSCOPY |
WO2014152389A1 (en) | 2013-03-14 | 2014-09-25 | Lumicell, Inc. | Imaging agent for detection of diseased cells |
CA3016468A1 (en) | 2016-03-08 | 2017-09-14 | Zebra Medical Technologies, Inc. | Non-invasive detection of skin disease |
WO2018201082A1 (en) | 2017-04-28 | 2018-11-01 | Zebra Medical Technologies, Inc. | Systems and methods for imaging and measurement of sarcomeres |
GB201707239D0 (en) | 2017-05-05 | 2017-06-21 | Univ Edinburgh | Optical system and method |
US11426075B1 (en) | 2017-08-23 | 2022-08-30 | Lumicell, Inc. | System and method for residual cancer cell detection |
CN107678153B (zh) | 2017-10-16 | 2020-08-11 | 苏州微景医学科技有限公司 | 光纤束图像处理方法和装置 |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5810032A (ja) * | 1981-07-09 | 1983-01-20 | オリンパス光学工業株式会社 | 内視鏡 |
JPS63155115A (ja) * | 1986-12-19 | 1988-06-28 | Olympus Optical Co Ltd | 立体観察電子内視鏡 |
JPH04138127A (ja) * | 1990-09-28 | 1992-05-12 | Olympus Optical Co Ltd | 内視鏡網目画像軽減装置 |
JP3078085B2 (ja) * | 1991-03-26 | 2000-08-21 | オリンパス光学工業株式会社 | 画像処理装置および画像処理方法 |
JPH08191439A (ja) | 1995-01-10 | 1996-07-23 | Fukuda Denshi Co Ltd | 内視画像補正方法及び装置 |
US6080104A (en) * | 1995-05-16 | 2000-06-27 | Asahi Kogaku Kogyo Kabushiki Kaisha | Electronic endoscope system |
US6190308B1 (en) * | 1995-08-17 | 2001-02-20 | Karl Storz Gmbh & Co., Kg | Endoscopic video system for correcting a video image of an object to be studied |
JP4008184B2 (ja) * | 1996-03-06 | 2007-11-14 | 富士フイルム株式会社 | 蛍光検出装置 |
US5764409A (en) * | 1996-04-26 | 1998-06-09 | Alpha Innotech Corp | Elimination of vibration by vibration coupling in microscopy applications |
WO1997042600A1 (fr) * | 1996-05-02 | 1997-11-13 | Andromis S.A. | Procede de traitement d'images obtenues par fibres multicoeurs ou multifibres, en particulier d'images endoscopiques |
JP3176574B2 (ja) * | 1997-10-15 | 2001-06-18 | 住友電気工業株式会社 | 光ファイバ観察装置および光ファイバ融着接続装置 |
AU2001247275A1 (en) * | 2000-03-06 | 2001-09-17 | Gene Logic, Inc. | Virtual gel profiling system |
US6885801B1 (en) * | 2001-12-06 | 2005-04-26 | Clear Image Technology Llc | Enhancement of fiber based images |
-
2002
- 2002-07-18 FR FR0209100A patent/FR2842628B1/fr not_active Expired - Lifetime
-
2003
- 2003-07-11 CA CA2491751A patent/CA2491751C/fr not_active Expired - Lifetime
- 2003-07-11 ES ES03750802.5T patent/ES2687644T3/es not_active Expired - Lifetime
- 2003-07-11 CN CN038210312A patent/CN1679053B/zh not_active Expired - Lifetime
- 2003-07-11 WO PCT/FR2003/002197 patent/WO2004010377A1/fr active Application Filing
- 2003-07-11 AU AU2003269019A patent/AU2003269019B2/en not_active Expired
- 2003-07-11 KR KR10-2005-7000839A patent/KR20050021493A/ko not_active Application Discontinuation
- 2003-07-11 JP JP2004522235A patent/JP4485947B2/ja not_active Expired - Lifetime
- 2003-07-11 US US10/520,917 patent/US7903848B2/en active Active
- 2003-07-11 EP EP03750802.5A patent/EP1523731B1/fr not_active Expired - Lifetime
-
2005
- 2005-01-05 IL IL166152A patent/IL166152A/en unknown
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108882831A (zh) * | 2016-03-31 | 2018-11-23 | 国立大学法人东北大学 | 光学成像装置 |
WO2019076267A1 (zh) * | 2017-10-16 | 2019-04-25 | 苏州微景医学科技有限公司 | 光纤束图像分析方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
KR20050021493A (ko) | 2005-03-07 |
FR2842628B1 (fr) | 2004-09-24 |
CA2491751A1 (fr) | 2004-01-29 |
JP2005532884A (ja) | 2005-11-04 |
AU2003269019B2 (en) | 2009-09-10 |
EP1523731B1 (fr) | 2018-06-27 |
AU2003269019A1 (en) | 2004-02-09 |
IL166152A (en) | 2010-12-30 |
FR2842628A1 (fr) | 2004-01-23 |
IL166152A0 (en) | 2006-01-15 |
CN1679053B (zh) | 2013-07-17 |
EP1523731A1 (fr) | 2005-04-20 |
WO2004010377A1 (fr) | 2004-01-29 |
ES2687644T3 (es) | 2018-10-26 |
JP4485947B2 (ja) | 2010-06-23 |
CA2491751C (fr) | 2015-10-13 |
US20050207668A1 (en) | 2005-09-22 |
US7903848B2 (en) | 2011-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1679053A (zh) | 用于处理通过包含多个光纤的波导获取的图像的方法 | |
Li et al. | Continuous depth map reconstruction from light fields | |
US7623728B2 (en) | Method and product for processing digital images | |
JP7252190B2 (ja) | 生物試料の向上された被写界深度の合成2d画像を生成するシステム | |
KR20140045331A (ko) | 현미경 이미징에서의 빠른 오토-포커스 | |
Ding et al. | Automatic macroscopic density artefact removal in a Nissl‐stained microscopic atlas of whole mouse brain | |
CN106257537B (zh) | 一种基于光场信息的空间深度提取方法 | |
Kervrann et al. | Robust incremental compensation of the light attenuation with depth in 3D fluorescence microscopy | |
EP3830628B1 (en) | Device and process for capturing microscopic plenoptic images with turbulence attenuation | |
Khan et al. | High-density single shot 3D sensing using adaptable speckle projection system with varying preprocessing | |
CN111311562A (zh) | 虚焦图像的模糊度检测方法及装置 | |
JP2004133919A (ja) | 擬似3次元画像生成装置および生成方法並びにそのためのプログラムおよび記録媒体 | |
JP2007322259A (ja) | エッジ検出方法、装置、及びプログラム | |
Sysko et al. | From image to data using common image‐processing techniques | |
EP2693397B1 (en) | Method and apparatus for noise reduction in an imaging system | |
CN108734674B (zh) | 一种改进nas-rif的oct图像盲复原方法 | |
Schmalz et al. | A graph-based approach for robust single-shot structured light | |
JP3572304B2 (ja) | 医療分野におけるフォトレタッチソフトウェアを用いた画像解析方法 | |
JP2023547454A (ja) | 画像フレームシーケンスの信号対雑音比を改善するための方法および画像処理デバイス | |
JP2018081378A (ja) | 画像処理装置、撮像装置、画像処理方法および画像処理プログラム | |
WO2021009804A1 (ja) | 閾値の学習方法 | |
Kim et al. | Fast Noise Reduction Approach in Multifocal Multiphoton Microscopy Based on Monte-Carlo Simulation | |
JP3609284B2 (ja) | 検出装置 | |
Shim et al. | De-ghosting in High Dynamic Range Imaging Based on Intensity Scaling Cue. | |
Simon et al. | A novel method of pre-processing using dental x-ray images by adaptive Morpho Histo Wavelet Denoising (AMHW) method |
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 | ||
CX01 | Expiry of patent term | ||
CX01 | Expiry of patent term |
Granted publication date: 20130717 |