CN102525381A - 图像处理装置、图像处理方法及计算机可读取的记录装置 - Google Patents
图像处理装置、图像处理方法及计算机可读取的记录装置 Download PDFInfo
- Publication number
- CN102525381A CN102525381A CN201110420174XA CN201110420174A CN102525381A CN 102525381 A CN102525381 A CN 102525381A CN 201110420174X A CN201110420174X A CN 201110420174XA CN 201110420174 A CN201110420174 A CN 201110420174A CN 102525381 A CN102525381 A CN 102525381A
- Authority
- CN
- China
- Prior art keywords
- image processing
- processing apparatus
- image
- candidate point
- border circular
- 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
- 238000012545 processing Methods 0.000 title claims abstract description 78
- 238000003672 processing method Methods 0.000 title claims abstract description 7
- 239000000284 extract Substances 0.000 claims abstract description 40
- 238000000605 extraction Methods 0.000 claims abstract description 38
- 238000001727 in vivo Methods 0.000 claims description 53
- 238000011156 evaluation Methods 0.000 claims description 47
- 238000012360 testing method Methods 0.000 claims description 40
- 238000000034 method Methods 0.000 claims description 35
- 238000004364 calculation method Methods 0.000 claims description 26
- 238000003860 storage Methods 0.000 claims description 26
- 238000001514 detection method Methods 0.000 claims description 18
- 230000011218 segmentation Effects 0.000 claims description 4
- 238000004040 coloring Methods 0.000 claims description 3
- 238000003709 image segmentation Methods 0.000 claims description 3
- 238000005286 illumination Methods 0.000 claims description 2
- 230000001678 irradiating effect Effects 0.000 claims 1
- 230000009471 action Effects 0.000 description 17
- 238000010586 diagram Methods 0.000 description 12
- 230000008569 process Effects 0.000 description 7
- 238000000205 computational method Methods 0.000 description 5
- 101001021281 Homo sapiens Protein HEXIM1 Proteins 0.000 description 4
- 101000693265 Homo sapiens Sphingosine 1-phosphate receptor 1 Proteins 0.000 description 4
- 102100025750 Sphingosine 1-phosphate receptor 1 Human genes 0.000 description 4
- 230000004075 alteration Effects 0.000 description 4
- 239000002775 capsule Substances 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 101000693269 Homo sapiens Sphingosine 1-phosphate receptor 3 Proteins 0.000 description 2
- 102000004137 Lysophosphatidic Acid Receptors Human genes 0.000 description 2
- 108090000642 Lysophosphatidic Acid Receptors Proteins 0.000 description 2
- 102100025747 Sphingosine 1-phosphate receptor 3 Human genes 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000007596 consolidation process Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 210000004877 mucosa Anatomy 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 241001062009 Indigofera Species 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 238000003708 edge detection Methods 0.000 description 1
- 239000012467 final product Substances 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- -1 pixel.Thus Substances 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 230000009747 swallowing Effects 0.000 description 1
- 210000001519 tissue Anatomy 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/44—Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/12—Edge-based segmentation
-
- 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/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
Abstract
本发明提供一种图像处理装置、图像处理方法及计算机可读取的记录装置,图像处理装置具有:候选点提取单元部,其从图像中提取候选点,该候选点是构成表示泡的圆形区域的点的候选;以及圆形区域检测单元部,其根据所述候选点所具有的信息,检测所述图像内的圆形区域。
Description
本申请享受2010年12月16日申请的日本专利申请编号2010-280716的优先权的利益,该日本专利申请的全部内容被引用到本申请中。
技术领域
本发明涉及从图像中检测映出泡的区域的图像处理装置、图像处理方法以及计算机可读取的记录装置。
背景技术
以往,作为导入患者等被检者体内以非侵袭的方式对活体内进行观察的医用观察装置,内窥镜得到广泛普及。近年来,开发了如下的吞入型内视镜(胶囊型内窥镜):在胶囊型的壳体内部收纳摄像装置和通信装置等,以无线的方式向体外发送由摄像装置摄像而得到的图像数据。通过这种医用观察装置对活体的管腔内进行摄像而得到的一系列图像(活体内图像)的数量庞大(几万张以上),并且,针对各活体内图像的观察和诊断需要较多的经验。因此,希望具有辅助医师进行诊断的医疗诊断支援功能。
作为实现这种功能的图像识别技术之一,提出了如下技术:在医疗诊断时,从活体内图像中自动检测不需要观察的区域并去除。这里,不需要观察的区域是指作为活体组织的粘膜以外的区域,例如举出映出在管腔内产生的气泡(以下也简称为“泡”)的区域(以下称为“泡区域”)。例如,在日本特开2007-313119号公报中公开了如下方法:根据像素的梯度强度和泡图像的特征计算与预先设定的泡模型之间的相关值,将相关值为规定阈值以上的区域检测为泡区域。
发明内容
本发明的一个方式的图像处理装置具有:候选点提取部,其从图像中提取候选点,该候选点是构成表示泡的圆形区域的点的候选;以及圆形区域检测部,其根据所述候选点所具有的信息,检测所述图像内的圆形区域。
本发明的另一个方式的图像处理方法包括以下步骤:从图像中提取候选点,该候选点是构成表示泡的圆形区域的点的候选,根据所述候选点所具有的信息,检测所述图像内的圆形区域。
本发明的又一个方式的计算机可读取的记录装置记录有可执行的程序,其中,所述程序指示处理器执行以下步骤:从图像中提取候选点,该候选点是构成表示泡的圆形区域的点的候选,根据所述候选点所具有的信息,检测所述图像内的圆形区域。
参照附图阅读以下本发明的详细说明,能够更进一步理解以上所述的内容、本发明的其他特征、优点以及技术上和产业上的意义。
附图说明
图1是示出本发明的实施方式1的图像处理装置的结构的框图。
图2是示出图1所示的初始候选点提取部的结构的框图。
图3是示出图1所示的圆形区域检测部的结构的框图。
图4是说明泡区域的构造的图。
图5是示出图1所示的图像处理装置的动作的流程图。
图6是示出图1所示的图像处理装置处理的活体内图像的一例的示意图。
图7是说明像素值的梯度为凸的方向上的梯度强度的计算方法的图。
图8是示出从图6所示的活体内图像中提取出的边缘像素的示意图。
图9是示出从图8所示的图像中作为初始候选点而提取出的边缘像素的示意图。
图10是示出根据距离值(R值)对图6所示的活体内图像进行分割后的分割区域的图。
图11是示出距离值(R值)与圆形区域的尺寸之间的对应的表。
图12是说明基于RANSAC的圆形区域的估计方法的图。
图13是说明边界区域的设定方法的图。
图14是说明表示边界区域中包含的边缘像素的连续性的评价值的计算方法的图。
图15是说明表示边界区域中包含的边缘像素相对于圆形区域的中心坐标的相对性的评价值的计算方法的图。
图16是说明表示边界区域中包含的边缘像素的圆形性的评价值的计算方法的图。
图17是示出本发明的实施方式2的图像处理装置的结构的框图。
图18是示出图17所示的投票部的结构的框图。
图19是示出图17所示的图像处理装置的动作的流程图。
图20是说明初始候选点针对参数空间的投票方法的图。
图21是说明初始候选点针对参数空间的投票方法的图。
图22是示出本发明的实施方式3的图像处理装置的结构的框图。
图23是示出图22所示的图像处理装置的动作的流程图。
图24是说明图22所示的图像处理装置的泡区域的判定方法的图。
具体实施方式
下面,参照附图说明本发明的实施方式的图像处理装置。另外,本发明不被这些实施方式限定。并且,在各附图的记载中,对相同部分标注相同标号来示出。
以下说明的实施方式的图像处理装置对通过内窥镜或胶囊内窥镜等医用观察装置摄像而得到的图像进行处理,具体而言,进行如下处理:从对被检者的活体(管腔)内进行摄像而得到的图像中检测映出泡的泡区域。在以下的说明中,实施图像处理的活体内图像例如是在各像素中具有针对R(红)、G(绿)、B(蓝)的各颜色成分的256灰阶的像素等级(像素值)的彩色图像。另外,本发明的图像处理装置不限于活体内图像,在从一般图像中检测泡区域的情况下也能够应用。
并且,由于表面张力的作用,在以下的实施方式中作为检测对象的泡区域通常具有圆形形状、扁圆形状或椭圆形状等类圆形状、或者这些圆周的一部分形状(圆弧等)的平面外观。因此,在本申请中,将正圆、扁圆和椭圆等类圆形状称为“圆形”,将圆形的一部分(例如圆弧)的形状称为“圆弧形状”。
实施方式1
图1是示出本发明的实施方式1的图像处理装置的结构的框图。如图1所示,图像处理装置1具有对图像处理装置1整体的动作进行控制的控制部10、图像取得部11、输入部12、显示部13、存储部14和运算部15。
图像取得部11取得通过医用观察装置摄像而得到的活体内图像的图像数据。图像取得部11根据包含医用观察装置的系统的形态而适当构成。例如,在医用观察装置是胶囊内窥镜、且使用可移动型记录介质在与医用观察装置之间交接图像数据的情况下,图像取得部11由装卸自如地装配该记录介质并读出所保存的活体内图像的图像数据的读出装置构成。并且,在设置预先保存通过医用观察装置摄像而得到的活体内图像的图像数据的服务器的情况下,图像取得部11由与服务器连接的通信装置等构成,与服务器进行数据通信来取得活体内图像的图像数据。或者,还可以利用经由电缆从内窥镜等医用观察装置输入图像信号的接口装置等构成图像取得部11。
输入部12例如由键盘、鼠标、触摸面板、各种开关等实现,将输入信号输出到控制部10。
显示部13由LCD或EL显示器等显示装置实现,在控制部10的控制下,显示包含活体内图像的各种画面。
存储部14由可更新记录的闪存等ROM或RAM这样的各种IC存储器、内置或者利用数据通信端子连接的硬盘、或CD-ROM等信息记录装置及其读取装置等实现。存储部14除了存储由图像取得部11取得的活体内图像的图像数据以外,还存储用于使图像处理装置1动作并使图像处理装置1执行各种功能的程序、以及在该程序的执行中使用的数据等。具体而言,存储部14存储用于从活体内图像中检测泡区域的图像处理程序141。
运算部15由CPU等硬件实现,通过读入图像处理程序141对活体内图像的图像数据进行处理,进行用于从活体内图像中检测泡区域的各种运算处理。运算部15具有:初始候选点提取部16,其根据泡区域所具有的特征,从活体内图像中提取初始候选点,该初始候选点是构成表示泡的圆形区域的点的候选;以及圆形区域检测部17,其根据提取出的初始候选点所具有的信息,从活体内图像中检测圆形区域。
控制部10由CPU等硬件实现,通过读入存储在存储部14中的各种程序,根据从图像取得部11输入的图像数据和从输入部12输入的操作信号等,进行针对构成图像处理装置1的各部的指示或数据转送等,对图像处理装置1整体的动作进行统一控制。
接着,对运算部15的详细结构进行说明。
图2是示出初始候选点提取部16的结构的框图。初始候选点提取部16具有:梯度强度计算部161,其计算活体内图像中包含的各像素的像素值的梯度强度;边缘检测部162,其根据梯度强度检测边缘;以及镜面反射区域检测部163,其从活体内图像中检测镜面反射区域。
作为梯度强度计算部161执行的梯度强度的计算方法,能够使用一阶微分(索贝尔滤波器等)或二阶微分等已知的各种手法。作为一例,梯度强度计算部161具有凸部梯度强度计算部161a,通过进行二阶微分来计算像素值的梯度为凸的方向上的梯度强度。
边缘检测部162通过提取具有规定阈值以上的梯度强度的像素来检测边缘。另外,作为边缘的检测方法,除此之外,能够使用已知的各种手法。例如,可以使用坎尼(Canny)的边缘检测算法(参考:CG-ARTS協会,デイジタル画像処理,p.208~p.210)等。
镜面反射区域检测部163具有:高亮度区域检测部163a,其计算活体内图像中包含的各像素的亮度值,检测具有规定值以上的亮度的高亮度区域;以及白色区域检测部163b,其根据各像素的像素值,检测具有包含在表示白色的颜色范围中的像素值的像素的区域(白色区域)。作为表示白色的颜色范围,例如使用将在活体内产生的光晕部分的RGB值转换为HSI值(H:色相;S:彩度;I:明度;参考:CG-ARTS協会,デイジタル画像処理,p.64~p.68)后的值的范围(例如S≤0.2、I≥0.8)。这种颜色范围作为教师数据预先存储在存储部14中,根据需要,白色区域检测部163b从存储部14中读出教师数据进行参照。
图3是示出圆形区域检测部17的结构的框图。圆形区域检测部17具有:估计用候选点提取部18,其从初始候选点内提取用于估计圆形区域的多个候选点(估计用候选点);圆形估计部19,其使用提取出的估计用候选点的坐标,通过稳健估计(robustestimation)来估计圆形区域;评价值计算部20,其计算用于评价所估计出的圆形区域是否实际存在于活体内图像内的评价值;以及判定部21,其将该评价值与规定阈值进行比较,评价圆形区域是否实际存在于活体内图像内。
估计用候选点提取部18具有:圆形区域存储部181,其存储与已检测的圆形区域有关的信息;以及圆形区域尺寸估计部182,其估计圆形区域的大小。圆形区域尺寸估计部182计算表示在活体内图像中映出的被摄体的进深的距离信息,根据该距离信息估计圆形区域的尺寸。
评价值计算部20具有:边界区域提取部201,其提取圆形区域的边界区域;以及边界区域信息计算部202,其取得与边界区域中包含的像素有关的信息(以下也称为边界区域信息)。作为边界区域信息,例如使用边界区域中的像素值的梯度强度的平均值。
这里,对作为运算部15的检测对象的泡区域进行说明。图4是说明泡区域的构造例的图。如图4(a)所示,泡100是包含空气等气体的液体(水分等)由于表面张力的作用而变成圆形的,包含水分等的膜101及其内部的气体102。在对这种泡100照射光并进行摄像时,得到图4(b)所示的泡区域110的图像。该泡区域110在外形上具有圆形形状或圆弧形状,作为内部构造,包含与膜101的厚度部分对应的在泡区域110的周缘部中包含的边界区域111、位于其内侧的内部区域112、以及在内部区域112的内侧的多个部位能够看到的镜面反射区域113。其中,边界区域111是由于照射光在膜101上散射而产生的明亮区域。另一方面,内部区域112是由于较多的照明光透过膜101而产生的区域,是整体较暗的区域。镜面反射区域113是照射光在膜101的表面进行镜面反射而映入光源从而产生光晕的亮度极高的区域。
接着,对图像处理装置1的动作进行说明。图5是示出图像处理装置1的动作的流程图。
首先,在步骤S10中,图像取得部11从外部取得活体内图像组,并将其存储在存储部14中。运算部15从存储部14中依次读出要实施图像处理的活体内图像。图6是示出从存储部14中读出的活体内图像的一例的示意图。如图6所示,在活体内图像120中,具有与活体的粘膜不同的特征的区域122有时重叠在粘膜区域121上映出。
在步骤S11中,梯度强度计算部161计算活体内图像120内的各像素中的像素值的梯度强度。具体而言,首先,凸部梯度强度计算部161a针对位于坐标(x,y)的像素的像素值P(x,y),使用下式(1)~式(4),计算水平方向上的二阶微分dH、垂直方向上的二阶微分dV、第1倾斜方向上的二阶微分dD1、第2倾斜方向上的二阶微分dD2(参照图7)。
这里,参数u是整数,参数u’是对(u/√2)进行四舍五入后的整数。这些参数u和u’表示二阶微分计算时的像素范围。即,通过将参数u设定为期望值,能够调节梯度成分的精度。例如,当减小参数u的值时,能够计算出空间频率高的梯度成分,当增大参数u的值时,能够计算出空间频率低的梯度成分。作为参数u的值,可以预先设定固定值,也可以根据活体内图像的特性而设定动态的值。或者,还可以从外部赋予期望值的参数u进行处理。
接着,凸部梯度强度计算部161a使用下式(5),从在不同方向计算出的二阶微分dH、dV、dD1、dD2内提取值为正的最大值,由此,计算梯度增加(即梯度相对于周围像素为凸)方向上的梯度强度(凸部梯度强度)gconv(x,y)。
gconv(x,y)=max(0,dH(x,y),dV(x,y),dD1(x,y),dD2(x,y))
(5)
这里,在图4所示的泡区域110的情况下,边界区域111的像素值大于内部区域112和泡区域110外侧的像素值,所以沿着边界区域111内的径向出现像素值的峰值。因此,通过计算凸部梯度强度gconv,能够检测径向上的像素值的上升。
在步骤S12中,边缘检测部162利用规定阈值对针对各像素计算出的凸部梯度强度gconv进行二值化,从而提取边缘像素。由此,提取作为泡区域的边界的可能性高的边缘。在图8所示的图像125中示出这样从活体内图像120中提取出的边缘像素GE。
在步骤S13中,镜面反射区域检测部163从活体内图像120中检测镜面反射区域113。具体而言,首先,高亮度区域检测部163a根据活体内图像120内的各像素的R值(R)、G值(G)和B值(B),使用下式(6)计算亮度值Y(参考:CG-ARTS協会,デイジタル画像処理,p.299)。
Y=0.3×R+0.59×G+0.11×B (6)
接着,高亮度区域检测部163a提取亮度值Y为规定阈值以上的像素作为高亮度区域。进而,白色区域检测部163b取得作为高亮度区域而提取出的像素的R值、G值和B值,参照教师数据判定这些颜色成分是否包含在表示白色的颜色范围内。然后,白色区域检测部163b提取具有包含在表示白色的颜色范围中的颜色成分的像素,作为镜面反射区域。
在步骤S14中,初始候选点提取部16提取在步骤S12中提取出的边缘像素内、在步骤S13中提取出的镜面反射区域附近、即存在于镜面反射区域的规定范围内的边缘像素的坐标,作为初始候选点。这是因为,镜面反射区域附近的边缘是泡区域110的边界区域111的可能性高(参照图4)。在图9所示的图像126中示出作为初始候选点提取出的边缘像素GE。
在步骤S15中,圆形区域尺寸估计部182计算活体内图像120内的距离信息。具体而言,首先,圆形区域尺寸估计部182取得活体内的吸收和散射最少的R值,作为表示各像素的进深的距离值。这是因为,被摄体与摄像元件之间的距离越短(即被摄体越靠近画面的近前侧),R值越大。接着,圆形区域尺寸估计部182根据各像素的距离值,将活体内图像120分割为多个区域A~D(图10)。作为区域分割法,例如能够使用将具有类似的特征量(像素值、纹理特征量等)且在空间上接近的像素作为集合来分割区域的区域合并法(参考:CG-ARTS協会,デイジタル画像処理,p.196)、或者作为聚类手法之一的K-means法(参考:CG-ARTS協会,デイジタル画像処理,p.232)等已知的方法。进而,圆形区域尺寸估计部182按照每个分割区域A~D计算距离值的平均值(平均距离值),作为距离信息。
在接下来的步骤S16中,圆形区域尺寸估计部182估计初始候选点的坐标中的圆形区域的尺寸。这里,在活体内图像中,通常从近景到远景均映出。并且,一般地,以位于远景的被摄体比位于近景的被摄体小的方式映出。因此,根据图像内的距离信息,能够以某个程度估计圆形区域的尺寸。
具体而言,圆形区域尺寸估计部182参照图11所示的表128,按照每个分割区域A~D取得与平均距离值对应的圆形区域的尺寸的范围(最小半径rMIN和最大半径rMAX)。另外,表128预先生成并存储在存储部14中。
在步骤S17中,估计用候选点提取部18从初始候选点内提取多个估计用候选点。另外,在实施方式1中,由于使用RANSAC(random sampling consensus:随机抽样一致性)的手法,所以随机提取估计用候选点。具体而言,首先,如图12所示,估计用候选点提取部18随机提取初始候选点内的一点(候选点GE0)。接着,估计用候选点提取部18求出区域127,该区域127内包于以候选点GE0为原点、以针对该候选点GE0所属的分割区域A~D取得的圆形区域的最大半径rMAX或比其稍大的长度(rMAX+α)为半径的圆中。进而,估计用候选点提取部18从该区域127内部随机提取规定个数的候选点GE1、GE2、...,作为估计用候选点。
在步骤S18中,圆形估计部19根据估计用候选点GE1、GE2、...估计圆形区域。具体而言,圆形估计部19生成通过估计用候选点GE1、GE2、...的坐标的圆的联立方程式,利用最小二乘法求解该联立方程式,由此计算所估计出的圆形区域。然后,将与计算出的圆形区域之间的距离为规定阈值以下的估计用候选点判定为内围层(inlier),将除此之外的估计用候选点判定为外围层(outlier)。进而,圆形估计部19使用被判定为内围层的估计用候选点的坐标,同样地再次计算圆形区域。
在步骤S19中,如图13所示,边界区域提取部201提取圆形区域130的圆周131的距离d1的范围内的区域,将该区域设定为边界区域132。然后,对边界区域132中包含的像素数N进行计数。
在步骤S20中,边界区域信息计算部202计算用于判定圆形区域130是否实际存在于活体内图像120内的评价值。具体而言,边界区域信息计算部202根据边界区域132中包含的各像素的梯度强度的合计值和像素数N,计算梯度强度的平均值SAVE,通过下式(7)计算评价值E1。
E1=SAVE×k (7)
在式(7)中,k是规定系数(常数)。
在步骤S21中,判定部21根据评价值E1,判定圆形区域130是否实际存在于活体内图像120内。在评价值E1的值大于规定阈值的情况下,判定为圆形区域130实际存在于活体内图像120内。
圆形区域检测部17在判定为圆形区域130实际存在于活体内图像120内的情况下(步骤S21:是),视为检测到圆形区域,在圆形区域存储部181中存储与判定对象的圆形区域130有关的信息(圆形的中心坐标和半径、以及其中包含的像素的像素值等),并且作为表示泡区域的信息存储在存储部14中(步骤S22)。
在步骤S23中,圆形区域检测部17判定作为针对一个初始候选点的估计用候选点的提取处理次数是否超过规定阈值。在提取处理次数没有超过阈值的情况下(步骤S23:否),动作转移到步骤S17。该情况下,在第2次及以后执行的步骤S17和S18中,不提取位于已经判定为圆形的区域内的初始候选点。另一方面,在提取处理次数超过规定值的情况下(步骤S23:是),动作结束。
并且,圆形区域检测部17在判定为圆形区域130未实际存在于活体内图像120内的情况下(步骤S21:否),判定在步骤S17~S21的一连串处理中持续没有检测到圆形区域的次数(连续未检测次数)是否超过规定值(步骤S24)。在连续未检测次数没有超过规定值的情况下(步骤S24:否),动作转移到步骤S23。另一方面,在连续未检测次数超过规定值的情况下(步骤S24:是),动作结束。
如以上说明的那样,根据实施方式1,根据从活体内图像中提取出的边缘像素来估计圆形区域,并且,根据计算出的圆形区域的边界区域所具有的信息判定所估计出的圆形区域是否实际存在于活体内图像内,所以能够进行高精度的泡区域的检测。
并且,根据实施方式1,由于根据泡区域所具有的构造上的特征提取从活体内图像中提取出的边缘像素,使用该边缘像素估计圆形区域,所以能够抑制运算量,能够高效地检测圆形区域。
进而,根据实施方式1,由于根据距离信息估计圆形区域的尺寸,根据该尺寸限制此后的圆形区域的计算处理中的运算范围,所以能够减少运算量(计算成本),并且能够抑制圆形区域的误检测。
另外,在上述说明中,提取高亮度且白色的区域作为镜面反射区域,但是,也可以提取高亮度区域或白色区域的任意一方作为镜面反射区域进行处理。
并且,在上述说明中,通过稳健估计来估计圆形区域,但是,也可以通过同样的手法来估计椭圆形区域。
作为用于圆形区域130的评价的评价值,除了上述评价值E1以外,也可以使用各种值。下面,在变形例1-1~1-6中,说明评价值计算部20(边界区域信息计算部202)计算的评价值的例子。
变形例1-1
使用图14所示的边界区域132中包含的全部像素数N和其中的边缘像素GE的个数NEDG,通过下式(8)给出变形例1-1的评价值E2。
E2=(NEDG/N)×k (8)
在式(8)中,k是规定系数(常数)。
该情况下,在评价值E2的值大于规定阈值时,判定为圆形区域130是实际存在于活体内图像120内的泡区域。
变形例1-2
使用边界区域132中包含的全部像素数N和边界区域132中的亮度值的合计值YSUM,通过下式(9)给出变形例1-2的评价值E3。另外,通过式(6)计算各像素的亮度值。
E3=(YSUM/(N×YMAX))×k (9)
在式(9)中,YMAX是边界区域中的亮度值的最大值。
该情况下,在评价值E3的值大于规定阈值时,判定为圆形区域130是实际存在于活体内图像120内的泡区域。
变形例1-3
或者,也可以根据边界区域的RGB值来计算评价值E4。该情况下,预先取得泡区域的RGB颜色范围并存储在存储部14中。
首先,边界区域信息计算部202取得边界区域132中包含的各像素的R值、G值和B值,针对各颜色成分计算像素间的平均值。边界区域信息计算部202对由此得到的平均R值、平均G值和平均B值与泡区域的RGB颜色范围进行比较,计算色差ΔE。使用该色差ΔE,通过下式(10)给出评价值E4。
E4=(1/ΔE)×k (10)
该情况下,在评价值E4大于规定阈值时,判定为圆形区域130是实际存在于活体内图像120内的泡区域。
变形例1-4
泡区域通常具有圆形或作为其一部分的圆弧形状。因此,也可以根据圆形区域130的形状信息(即在几何学上是否具有圆弧的特征),判定圆形区域130是否是泡区域。在变形例1-4中,作为这种几何学的特征,判定边界区域132中包含的边缘像素的连续性。
具体而言,首先,如图14所示,边界区域信息计算部202取得边界区域132中包含的边缘像素GE的个数NEDG1(即相当于边缘的面积)。并且,边界区域信息计算部202对边界区域132的边缘像素GE内、在周围8个方向的任意一个方向上存在相邻的边缘像素的边缘像素GE的个数NEDG2进行计数。例如,由于图14所示的边缘像素GE(1)在右侧斜上方存在相邻的边缘像素,所以被计数。另一方面,边缘像素GE(2)在周围8个方向上不存在相邻的边缘像素,所以不被计数。
使用这些值NEDG1和NEDG2,通过下式(11)给出评价值E5。
E5=(NEDG1×(1-k)+NEDG2×k)/N (11)
这里,在利用边界区域132的尺寸对评价值进行标准化时,使用边界区域132中包含的全部像素数N。并且,k是满足0<k<1的系数(常数)。
该情况下,在评价值E5大于规定阈值时,边界区域132中包含的边缘像素GE的连续性高,因此,判定为圆形区域130是实际存在于活体内图像120内的泡区域。
变形例1-5
表示泡区域的边界的边缘通常出现在关于圆的中心坐标而相对的圆周上的位置。因此,在边界区域132中包含的边缘像素GE关于圆形区域130的中心坐标具有较高的相对性的情况下,能够判断为该边缘像素GE构成泡区域的边界。这里,相对性高意味着,在相对于圆的中心坐标而相对的位置存在的边缘像素对儿较多。
具体而言,首先,如图15所示,边界区域信息计算部202取得从边界区域132中提取出的边缘像素GE的个数NEDG1。并且,边界区域信息计算部202对在边界区域132中在相对于圆形区域130的中心坐标C(x,y)而相对的位置存在其他边缘像素的边缘像素GE的个数NEDG3进行计数。例如,图15所示的边缘像素GE(3)在通过边缘像素GE(3)的圆形区域130的直径上存在其他边缘像素GT,所以被计数。
使用这些值NEDG1和NEDG3,通过下式(12)给出评价值E6。
E6=(NEDG1×(1-k)+NEDG3×k)/N (12)
在式(12)中,k是满足0<k<1的系数(常数)。
该情况下,在评价值E6大于规定阈值时,边界区域132中包含的边缘像素GE的相对性高,因此,判定为圆形区域130是实际存在于活体内图像120内的泡区域。
变形例1-6
作为表示圆形区域130的边界区域的形状信息的特征量,也可以使用表示边界区域132中的边缘像素的圆形性、即几何学的似圆度的评价值。图16是示出边界区域132的一部分的示意图。在计算特征量时,首先,边界区域信息计算部202对边界区域132中包含的边缘像素GE的个数NEDG1进行计数。并且,边界区域信息计算部202计算从各边缘像素GE到圆形区域130的圆周131的距离LE,对距离LE为规定阈值以内的边缘像素的个数NGC进行计数。使用这些值NEDG1和NGC,通过下式(13)给出评价值E7。
E7=(NEDG1×(1-k)+NGC×k)/N (13)
在式(13)中,k是满足0<k<1的系数(常数)。
该情况下,在评价值E7大于规定阈值时,边界区域132中包含的边缘像素GE的圆形性高,因此,判定为圆形区域130是实际存在于活体内图像120内的泡区域。
实施方式2
接着,对本发明的实施方式2进行说明。图17是示出实施方式2的图像处理装置的结构的框图。如图17所示,图像处理装置2具有运算部30。其他结构与图1所示的结构相同。
运算部30具有初始候选点提取部16和圆形区域检测部31。其中,初始候选点提取部16的结构和动作与在实施方式1中说明的结构和动作相同。
圆形区域检测部31根据初始候选点提取部16提取出的初始候选点,通过霍夫(Hough)转换来检测圆形区域。更详细地讲,圆形区域检测部31具有:投票部32,其在由圆的半径和圆的中心坐标构成的参数空间内对初始候选点进行投票;投票评价值计算部33,其根据参数空间中的投票频度,计算用于检测圆形区域的评价值;以及判定部34,其根据评价值进行圆形区域的判定。
图18是示出投票部32的详细结构的框图。投票部32具有:权重决定部321,其决定对在参数空间内被投票的初始候选点赋予的权重;参数空间范围决定部322,其决定参数空间的范围;以及初始候选点投票部323,其在所决定的范围内在参数空间内对初始候选点进行投票。
权重决定部321根据各初始候选点所具有的信息来决定权重,使得具有与泡类似的特征的初始候选点的权重大。例如,权重决定部321决定与各初始候选点的梯度强度对应的权重。
参数空间范围决定部322具有:距离信息计算部322a,其计算表示在活体内图像中映出的被摄体的进深的距离信息;圆形区域尺寸估计部322b,其根据该距离信息估计作为检测对象的圆形区域的尺寸;以及范围决定部322c,其根据圆形区域的尺寸决定参数空间的范围。
接着,对图像处理装置2的动作进行说明。图19是示出图像处理装置2的动作的流程图。图19所示的步骤内的步骤S10~S14与在实施方式1中说明的步骤相同。并且,在图20所示的图像135中示出从活体内图像120(图6)中提取出的初始候选点GE(x1,y1)、GE(x2,y2)、GE(x3,y3)、...。
在步骤S31中,权重决定部321使用各初始候选点的梯度强度g和下式(14),计算对该初始候选点赋予的权重W1。
W1=(g/gMAX)×k (14)
在式(14)中,gMAX是作为处理对象的活体内图像中的像素值的梯度强度的最大值,k是规定系数(常数)。梯度强度g越大,该权重W1的值越大。
在步骤S32中,距离信息计算部322a取得活体内图像的距离信息(进深信息)。另外,距离信息的取得方法与在图5的步骤S15中说明的方法相同。
在步骤S33中,圆形区域尺寸估计部322b根据距离信息和图11所示的表128,估计各初始候选点的像素位置中的圆形区域的尺寸。另外,圆形区域的尺寸的估计方法与在图5的步骤S16中说明的方法相同。
在步骤S34中,范围决定部322c决定参数空间的范围。
图21是示出参数空间的示意图。在该三维参数空间中,a轴与图像135的X轴对应,b轴与图像135的Y轴对应。并且,r轴表示作为检测对象的圆形区域的半径。因此,具体而言,关于a轴的范围,最小值为0,最大值为活体内图像的X方向上的尺寸。并且,关于b轴的范围,最小值为0,最大值为活体内图像的Y方向上的尺寸。进而,根据在步骤S33中按照每个区域A~D(参照图10)估计出的圆形区域的尺寸的范围(参照图11),决定r轴的范围。即,将针对区域A~D估计出的最小半径rMIN内的最小值riMIN设定为r轴的最小值,将针对区域A~D估计出的最大半径rMAX内的最大值riMAX设定为r轴的最大值。
在步骤S35中,初始候选点投票部323在参数空间内对初始候选点GE进行投票。具体而言,首先,初始候选点投票部323生成图21所示的参数空间。接着,初始候选点投票部323沿着扫描线方向从图像135中依次提取初始候选点的坐标GE(xn,yn)。然后,初始候选点投票部323在参数空间的r=rm平面中,对位于以坐标(xn,yn,rm)为中心的半径rm的圆周上的点赋予权重W1并进行投票。初始候选点投票部323在riMIN~riMAX的范围内,例如一个像素一个像素地改变r的值(进行投票的平面和半径的值),并进行这种投票。
在步骤S36中,投票评价值计算部33根据参数空间中的投票频度,计算用于检测圆形区域的评价值。这里,图21所示的rm平面上的点C0是以C1(x1,y1,rm)、C2(x2,y2,rm)、C3(x3,y3,rm)为中心的半径rm的圆交叉的点。换言之,点C0与这些点C1、C2、C3等距离,所以在图像135中,能够估计到存在以点C0为中心、半径为rm的圆形区域。即,参数空间中的投票频度高的坐标(a,b,r)表示图像135中的圆形区域。
作为评价值,例如使用通过高斯滤波器(参考:CG-ARTS協会,デイジタル画像処理,p.108~p.110)等对参数空间的各rm平面实施平滑处理并去除噪声后的处理结果(投票频度)。
在步骤S37中,判定部34提取评价值为规定阈值以上的参数空间上的坐标(a,b,r),将与该坐标对应的图像135上的区域(即圆的中心坐标和半径)检测为圆形区域。
进而,在步骤S38中,圆形区域检测部31输出与检测到的圆形区域有关的信息(圆形的中心坐标和半径、以及其中包含的像素的像素值等),作为表示泡区域的信息存储在存储部14中。
如以上说明的那样,根据实施方式2,由于在参数空间内对从活体内图像中提取出的边缘像素的坐标进行投票,根据该投票频度检测圆形区域,所以能够减少运算量,并且能够提高表示泡区域的圆形区域的检测精度。
并且,根据实施方式2,由于越是具有与泡类似的特征的初始候选点,在针对参数空间进行投票时对初始候选点赋予的权重越大,所以容易检测具有泡区域的特征的圆形区域,并且能够抑制误检测。
另外,在上述说明中,通过霍夫转换来检测圆形区域,但是,也可以通过同样的手法来检测椭圆形区域。该情况下,在与椭圆的中心坐标(a,b)、长径和短径对应的四维参数空间内对初始候选点进行投票即可。
变形例2-1
权重决定部321也可以使用初始候选点的梯度强度以外的特征量决定权重。例如,权重决定部321也可以使用利用式(6)针对各初始候选点计算出的亮度值和下式(15),决定权重W2。
W2=(Y/YMAX)×k (15)
在式(15)中,YMAX是初始候选点内的亮度值的最大值。另外,k是规定系数(常数)。
变形例2-2
或者,权重决定部321也可以根据各初始候选点的颜色信息决定权重。例如,权重决定部321取得各初始候选点中的R值、G值和B值,对这些值与泡区域的RGB颜色范围进行比较,计算色差ΔE。接着,权重决定部321使用下式(16)计算权重W3。色差ΔE越小,该权重W3为越大的值。
W3=(1/ΔE)×k (16)
另外,泡区域的RGB颜色范围预先存储在存储部14中。并且,k是规定系数(常数)。
实施方式3
接着,对本发明的实施方式3进行说明。图22是示出实施方式3的图像处理装置的结构的框图。如图22所示,图像处理装置3具有运算部40。其他结构与图1所示的结构相同。
除了初始候选点提取部16和圆形区域检测部17以外,运算部40还具有:区域分割部41,其将活体内图像分割为多个区域;以及区域判定部42,其判定在各分割区域整体中是否存在泡区域。另外,初始候选点提取部16和圆形区域检测部17的结构和动作与在实施方式1中说明的结构和动作相同。
这里,在活体内图像中,具有类似特征(颜色或纹理等特征)的区域之间大多是同一范畴(例如粘膜区域、泡区域等)。因此,在实施方式3中,根据颜色信息等特征量将活体内图像分割为多个区域,将单独检测到的圆形区域占据规定比例以上的分割区域视为在该分割区域整体中存在泡区域的区域。由此,能够将无法单独检测圆形的区域也检测为泡区域。
图23是示出图像处理装置3的动作的流程图。
首先,在步骤S10中,图像处理装置3取得活体内图像。另外,该步骤的详细情况与在实施方式1中说明的步骤相同。
在接下来的步骤S41中,如图24所示,区域分割部41将活体内图像136分割为多个区域137、138。作为分割方法,可以使用已知的各种方法。例如可以使用将具有类似的特征量(像素值、纹理特征量等)且在空间上接近的像素作为集合来分割区域的区域合并法(参考:CG-ARTS協会,デイジタル画像処理,p.196)、或者作为聚类手法之一的K-means法(参考:CG-ARTS協会,デイジタル画像処理,p.232)等。
在步骤S42中,初始候选点提取部16从活体内图像136中提取初始候选点。另外,初始候选点的提取方法与在图5的步骤S11~S14中说明的方法相同。
在步骤S43中,圆形区域检测部17根据提取出的初始候选点检测圆形区域139。另外,圆形区域的检测方法与在图5的步骤S15~S23中说明的方法相同。
在步骤S44中,区域判定部42根据圆形区域的检测结果,判定在各分割区域137、138整体中是否存在泡区域。即,首先,区域判定部42取得各分割区域的面积(或像素数)s1和在各分割区域内检测到的圆形区域的总面积(或像素数)s2,计算圆形区域的面积s2占据分割区域的面积s1的比例(面积率)s2/s1。在该面积率s2/s1为规定阈值以上的情况下,区域判定部42判定为在该分割区域整体中存在泡区域。另一方面,在面积率s2/s1小于规定阈值的情况下,区域判定部42判定为在该分割区域整体中并不存在泡区域。例如,关于分割区域137,由于被判定为泡区域的区域所占据的比例低(零),所以判定为不存在泡区域。另一方面,关于分割区域138,由于圆形区域139所占据的比例高,所以判定为在分割区域138整体中存在泡区域。
如以上说明的那样,在实施方式3中,根据单独判定为圆形区域的区域所占据的比例,判定在根据类似特征量分割后的区域整体中是否存在泡区域。因此,即使是无法直接检测到圆形的区域,也能够作为泡区域来提取。
另外,在运算部40中,代替圆形区域检测部17,也可以使用在实施方式2中说明的圆形区域检测部31。该情况下,在步骤S43中,与图19的步骤S31~S38同样地检测圆形区域。
以上说明的实施方式1~3的图像处理装置能够通过在个人计算机或工作站等计算机系统中执行记录在记录装置中的图像处理程序来实现。并且,也可以经由局域网、广域网(LAN/WAN)或因特网等公共线路,与其他计算机系统或服务器等设备连接来使用这种计算机系统。该情况下,实施方式1~3的图像处理装置也可以经由这些网络取得管腔内图像的图像数据,或者对经由这些网络连接的各种输出设备(阅读器或打印机等)输出图像处理结果,或者在经由这些网络连接的存储装置(记录装置及其读取装置等)中存储图像处理结果。
根据本发明的一个方式,由于从图像中提取构成圆形区域的一部分的候选点,根据该候选点所具有的信息检测圆形区域,所以能够高精度、高效地检测具有圆形等的区域。
本领域技术人员能够容易地导出进一步的效果和变形例。由此,本发明的更加广泛的方式不限于如上所表述的特定的详细和代表性的实施方式。因此,只要不脱离由附加权利要求及其均等物所定义的总括的发明概念的精神或范围,则能够进行各种变更。
Claims (34)
1.一种图像处理装置,其中,该图像处理装置具有:
候选点提取部,其从图像中提取候选点,该候选点是构成表示泡的圆形区域的点的候选;以及
圆形区域检测部,其根据所述候选点所具有的信息,检测所述图像内的圆形区域。
2.根据权利要求1所述的图像处理装置,其中,
所述候选点提取部根据泡区域所具有的内部构造,提取所述候选点。
3.根据权利要求2所述的图像处理装置,其中,
所述候选点提取部具有:
梯度强度计算部,其计算所述图像中包含的各像素的像素值的梯度强度;以及
边缘检测部,其根据所述梯度强度检测边缘。
4.根据权利要求3所述的图像处理装置,其中,
所述梯度强度计算部具有凸部梯度强度计算部,该凸部梯度强度计算部针对所述各像素,计算像素值的梯度为凸的方向上的梯度强度。
5.根据权利要求2所述的图像处理装置,其中,
所述候选点提取部具有镜面反射区域检测部,该镜面反射区域检测部从所述图像中检测在对泡照射照明光时进行镜面反射的区域。
6.根据权利要求5所述的图像处理装置,其中,
所述镜面反射区域检测部具有高亮度区域检测部,该高亮度区域检测部从所述图像中检测具有规定值以上的亮度的高亮度区域。
7.根据权利要求5所述的图像处理装置,其中,
所述镜面反射区域检测部具有白色区域检测部,该白色区域检测部从所述图像中检测具有包含在表示白色的颜色范围中的像素值的像素。
8.根据权利要求1所述的图像处理装置,其中,
所述圆形区域检测部具有:
估计用候选点提取部,其从由所述候选点提取部提取出的候选点内,提取用于估计圆形区域的多个估计用候选点;
圆形估计部,其根据由所述估计用候选点提取部提取出的所述多个估计用候选点,估计圆形区域;
评价值计算部,其根据所述圆形区域所具有的信息,计算用于判定该圆形区域是否实际存在于所述图像内的评价值;以及
判定部,其根据所述评价值,判定所述圆形区域是否实际存在于所述图像内。
9.根据权利要求8所述的图像处理装置,其中,
所述图像处理装置还具有存储部,该存储部存储与所述判定部判定为实际存在于所述图像内的圆形区域有关的信息,
所述估计用候选点提取部从存储在所述存储部中的圆形区域的外侧区域中,提取所述多个估计用候选点。
10.根据权利要求8所述的图像处理装置,其中,
所述估计用候选点提取部具有圆形区域尺寸估计部,该圆形区域尺寸估计部估计所述圆形区域的大小。
11.根据权利要求10所述的图像处理装置,其中,
所述圆形区域尺寸估计部根据所述图像中包含的各像素的像素值,计算表示在所述图像中映出的被摄体的进深的距离信息,根据该距离信息估计所述圆形区域的大小。
12.根据权利要求9所述的图像处理装置,其中,
所述评价值计算部具有:
边界区域提取部,其提取所述圆形区域的周缘部中包含的边界区域;以及
边界区域信息计算部,其取得与所述边界区域中包含的像素有关的信息。
13.根据权利要求12所述的图像处理装置,其中,
所述边界区域信息计算部计算所述边界区域中包含的像素的像素值的梯度强度的平均值。
14.根据权利要求12所述的图像处理装置,其中,
所述边界区域信息计算部对所述边界区域中包含的边缘像素的个数进行计数。
15.根据权利要求12所述的图像处理装置,其中,
所述边界区域信息计算部计算所述边界区域中包含的像素的亮度值的平均值。
16.根据权利要求12所述的图像处理装置,其中,
所述边界区域信息计算部取得所述边界区域中包含的像素的颜色信息。
17.根据权利要求12所述的图像处理装置,其中,
所述边界区域信息计算部计算表示所述边界区域中包含的边缘像素的形状的信息。
18.根据权利要求17所述的图像处理装置,其中,
所述边界区域信息计算部计算表示所述边缘像素的连续性的评价值。
19.根据权利要求18所述的图像处理装置,其中,
所述边界区域信息计算部对在所述边界区域中存在相邻边缘像素的边缘像素的数量进行计数。
20.根据权利要求17所述的图像处理装置,其中,
所述边界区域信息计算部计算表示所述边界区域中包含的边缘像素关于所述圆形区域的中心坐标的相对性的评价值。
21.根据权利要求20所述的图像处理装置,其中,
所述边界区域信息计算部对在所述边界区域中在关于所述圆形区域的中心坐标相对的位置处存在其他边缘像素的边缘像素的数量进行计数。
22.根据权利要求17所述的图像处理装置,其中,
所述边界区域信息计算部计算表示所述边界区域中包含的边缘像素的圆形性的评价值。
23.根据权利要求22所述的图像处理装置,其中,
所述边界区域信息计算部计算所述圆形区域的圆周与所述边界区域中包含的边缘像素之间的距离。
24.根据权利要求1所述的图像处理装置,其中,
所述圆形区域检测部具有:
投票部,其在由圆的半径和圆的中心坐标构成的参数空间内对所述候选点进行投票;以及
投票评价值计算部,其根据所述参数空间中的投票频度,计算用于检测圆形区域的评价值。
25.根据权利要求24所述的图像处理装置,其中,
所述投票部具有权重决定部,该权重决定部根据在所述参数空间内被投票的所述候选点所具有的信息,决定对该候选点赋予的权重。
26.根据权利要求25所述的图像处理装置,其中,
所述权重决定部根据所述候选点的像素值的梯度强度决定所述权重。
27.根据权利要求25所述的图像处理装置,其中,
所述权重决定部根据所述候选点的像素的亮度值决定所述权重。
28.根据权利要求25所述的图像处理装置,其中,
所述权重决定部根据所述候选点的像素所具有的颜色信息决定所述权重。
29.根据权利要求24所述的图像处理装置,其中,
所述投票部具有参数空间范围决定部,该参数空间范围决定部决定所述参数空间的范围。
30.根据权利要求29所述的图像处理装置,其中,
所述参数空间范围决定部具有:
距离信息计算部,其计算表示在所述图像中映出的被摄体的进深的距离信息;以及
圆形区域尺寸估计部,其根据所述距离信息,估计作为检测对象的圆形区域的大小。
31.根据权利要求1所述的图像处理装置,其中,
所述图像处理装置还具有:
区域分割部,其将所述图像分割为多个区域;以及
区域判定部,其根据各分割区域内的圆形区域的检测结果,判定在该分割区域整体中是否存在泡区域。
32.根据权利要求1所述的图像处理装置,其中,
所述图像是对活体的管腔内进行摄像而得到的活体内图像。
33.一种图像处理方法,其中,该图像处理方法包括以下步骤:
从图像中提取候选点,该候选点是构成表示泡的圆形区域的点的候选,
根据所述候选点所具有的信息,检测所述图像内的圆形区域。
34.一种计算机可读取的记录装置,其记录有可执行的程序,其中,
所述程序指示处理器执行以下步骤:
从图像中提取候选点,该候选点是构成表示泡的圆形区域的点的候选,
根据所述候选点所具有的信息,检测所述图像内的圆形区域。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JPJP2010-280716 | 2010-12-16 | ||
JP2010-280716 | 2010-12-16 | ||
JP2010280716A JP5576782B2 (ja) | 2010-12-16 | 2010-12-16 | 画像処理装置、画像処理方法、及び画像処理プログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102525381A true CN102525381A (zh) | 2012-07-04 |
CN102525381B CN102525381B (zh) | 2016-08-03 |
Family
ID=45421780
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110420174.XA Active CN102525381B (zh) | 2010-12-16 | 2011-12-15 | 图像处理装置、图像处理方法及计算机可读取的记录装置 |
Country Status (4)
Country | Link |
---|---|
US (1) | US8798344B2 (zh) |
EP (1) | EP2466541B1 (zh) |
JP (1) | JP5576782B2 (zh) |
CN (1) | CN102525381B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103900708A (zh) * | 2012-12-27 | 2014-07-02 | 杭州美盛红外光电技术有限公司 | 红外选择通知装置和红外选择通知方法 |
CN103900713A (zh) * | 2012-12-27 | 2014-07-02 | 杭州美盛红外光电技术有限公司 | 热像检测装置和热像检测方法 |
CN103900698A (zh) * | 2012-12-27 | 2014-07-02 | 杭州美盛红外光电技术有限公司 | 红外识别控制装置和红外识别控制方法 |
CN103900707A (zh) * | 2012-12-27 | 2014-07-02 | 杭州美盛红外光电技术有限公司 | 热像检测更新装置和热像检测更新方法 |
CN103900717A (zh) * | 2012-12-27 | 2014-07-02 | 杭州美盛红外光电技术有限公司 | 热像识别装置和热像识别方法 |
CN103900714A (zh) * | 2012-12-27 | 2014-07-02 | 杭州美盛红外光电技术有限公司 | 热像匹配装置和热像匹配方法 |
WO2014101802A1 (zh) * | 2012-12-27 | 2014-07-03 | Wang Hao | 热像挑选装置和热像挑选方法 |
WO2014101804A1 (zh) * | 2012-12-27 | 2014-07-03 | Wang Hao | 热像选择装置和热像选择方法 |
WO2014101803A1 (zh) * | 2012-12-27 | 2014-07-03 | Wang Hao | 红外选择装置和红外选择方法 |
CN104660922A (zh) * | 2013-11-25 | 2015-05-27 | 杭州美盛红外光电技术有限公司 | 热像整理装置、记录装置、整理系统和热像整理方法、记录方法 |
WO2015074625A1 (zh) * | 2013-11-25 | 2015-05-28 | 王浩 | 分析设置装置、处理装置和分析设置方法、处理方法 |
WO2015096824A1 (zh) * | 2013-12-29 | 2015-07-02 | 王浩 | 分析装置和分析方法 |
Families Citing this family (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9429940B2 (en) | 2011-01-05 | 2016-08-30 | Sphero, Inc. | Self propelled device with magnetic coupling |
US9218316B2 (en) | 2011-01-05 | 2015-12-22 | Sphero, Inc. | Remotely controlling a self-propelled device in a virtualized environment |
US8751063B2 (en) | 2011-01-05 | 2014-06-10 | Orbotix, Inc. | Orienting a user interface of a controller for operating a self-propelled device |
US10281915B2 (en) | 2011-01-05 | 2019-05-07 | Sphero, Inc. | Multi-purposed self-propelled device |
US9090214B2 (en) | 2011-01-05 | 2015-07-28 | Orbotix, Inc. | Magnetically coupled accessory for a self-propelled device |
US20150049177A1 (en) * | 2012-02-06 | 2015-02-19 | Biooptico Ab | Camera Arrangement and Image Processing Method for Quantifying Tissue Structure and Degeneration |
US8811675B2 (en) * | 2012-03-30 | 2014-08-19 | MindTree Limited | Circular object identification system |
EP2850512A4 (en) * | 2012-05-14 | 2016-11-16 | Sphero Inc | OPERATION OF A CALCULATION DEVICE BY DETECTING ROUNDED OBJECTS IN A PICTURE |
US10056791B2 (en) | 2012-07-13 | 2018-08-21 | Sphero, Inc. | Self-optimizing power transfer |
JP6230262B2 (ja) * | 2012-12-28 | 2017-11-15 | キヤノン株式会社 | 画像処理装置及び画像処理方法 |
EP3045104A4 (en) * | 2013-09-13 | 2017-04-26 | Olympus Corporation | Image processing device, method, and program |
US9829882B2 (en) | 2013-12-20 | 2017-11-28 | Sphero, Inc. | Self-propelled device with center of mass drive system |
CN104751445A (zh) * | 2013-12-29 | 2015-07-01 | 杭州美盛红外光电技术有限公司 | 热像分析配置装置和热像分析配置方法 |
CN103729632B (zh) * | 2014-01-22 | 2016-11-02 | 哈尔滨工业大学 | 一种基于连通区域滤波的圆形Mark点的定位方法 |
JP2015156937A (ja) * | 2014-02-24 | 2015-09-03 | ソニー株式会社 | 画像処理装置、画像処理方法、並びにプログラム |
JP2015156938A (ja) | 2014-02-24 | 2015-09-03 | ソニー株式会社 | 画像処理装置、画像処理方法 |
JP6581984B2 (ja) | 2015-01-21 | 2019-09-25 | Hoya株式会社 | 内視鏡システム |
US9824189B2 (en) * | 2015-01-23 | 2017-11-21 | Panasonic Intellectual Property Management Co., Ltd. | Image processing apparatus, image processing method, image display system, and storage medium |
JP6600356B2 (ja) | 2015-06-24 | 2019-10-30 | オリンパス株式会社 | 画像処理装置、内視鏡装置及びプログラム |
CN105139384B (zh) * | 2015-08-11 | 2017-12-26 | 北京天诚盛业科技有限公司 | 缺陷胶囊检测的方法和装置 |
KR102555096B1 (ko) * | 2016-06-09 | 2023-07-13 | 엘지디스플레이 주식회사 | 데이터 압축 방법 및 이를 이용한 유기 발광 다이오드 표시 장치 |
JPWO2018078806A1 (ja) * | 2016-10-28 | 2019-09-05 | オリンパス株式会社 | 画像処理装置、画像処理方法及び画像処理プログラム |
WO2018084025A1 (ja) * | 2016-11-07 | 2018-05-11 | オリンパス株式会社 | カプセル型内視鏡、受信装置、カプセル型内視鏡の作動方法、及びカプセル型内視鏡の作動プログラム |
KR101874738B1 (ko) * | 2017-02-10 | 2018-07-04 | 국민대학교산학협력단 | 영상처리를 이용하여 ldr 영상으로부터 hdr 영상을 생성하는 장치 및 방법 |
US10319108B2 (en) * | 2017-02-14 | 2019-06-11 | Jx Imaging Arts, Llc | System and method for machine vision object orientation measurement |
IT201700094994A1 (it) * | 2017-08-22 | 2019-02-22 | Milano Politecnico | Apparato e metodo per illuminare oggetti |
CN109190457B (zh) * | 2018-07-19 | 2021-12-03 | 北京市遥感信息研究所 | 一种基于大幅面遥感图像的油库集群目标快速检测方法 |
US20210291435A1 (en) * | 2020-03-19 | 2021-09-23 | Ricoh Company, Ltd. | Measuring apparatus, movable apparatus, robot, electronic device, fabricating apparatus, and measuring method |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5566249A (en) * | 1994-09-20 | 1996-10-15 | Neopath, Inc. | Apparatus for detecting bubbles in coverslip adhesive |
WO1999052072A1 (fr) * | 1998-04-07 | 1999-10-14 | Omron Corporation | Dispositif et procede de traitement d'images, support de stockage de programme pour traitement d'images, et dispositif d'inspection |
CN101065052A (zh) * | 2004-12-27 | 2007-10-31 | 奥林巴斯株式会社 | 医疗图像处理装置和医疗图像处理方法 |
CN101115435A (zh) * | 2005-02-15 | 2008-01-30 | 奥林巴斯株式会社 | 医用图像处理装置、管腔图像处理装置、管腔图像处理方法及用于上述装置和方法的程序 |
CN101150977A (zh) * | 2005-04-13 | 2008-03-26 | 奥林巴斯医疗株式会社 | 图像处理装置以及图像处理方法 |
US20090148014A1 (en) * | 2006-05-26 | 2009-06-11 | Olympus Corporation | Image processing apparatus, image processing method, and image processing program product |
JP2009136655A (ja) * | 2008-02-21 | 2009-06-25 | Olympus Corp | 画像処理装置、画像処理方法および画像処理プログラム |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4941182A (en) * | 1987-07-29 | 1990-07-10 | Phoenix Software Development Co. | Vision system and method for automated painting equipment |
US6415054B1 (en) * | 1997-07-15 | 2002-07-02 | Silverbrook Research Pty Ltd. | Target detection for dot region alignment in optical storage systems using ink dots |
US6529630B1 (en) * | 1998-03-02 | 2003-03-04 | Fuji Photo Film Co., Ltd. | Method and device for extracting principal image subjects |
JP4393016B2 (ja) * | 2000-06-30 | 2010-01-06 | 株式会社日立メディコ | 画像診断支援装置 |
US7231072B2 (en) * | 2002-02-27 | 2007-06-12 | Konica Corporation | Image processing method and image converting apparatus |
JP3946586B2 (ja) * | 2002-06-28 | 2007-07-18 | フジノン株式会社 | 縞画像の円形領域抽出方法 |
US8073233B2 (en) * | 2005-01-31 | 2011-12-06 | Olympus Corporation | Image processor, microscope system, and area specifying program |
JP4212564B2 (ja) * | 2005-02-28 | 2009-01-21 | ザイオソフト株式会社 | 画像処理方法および画像処理プログラム |
JP4420459B2 (ja) * | 2005-06-14 | 2010-02-24 | キヤノン株式会社 | 画像処理装置およびその方法 |
JP4761149B2 (ja) * | 2006-08-28 | 2011-08-31 | 富士フイルム株式会社 | 液体吐出装置及び気体処理方法 |
US8107726B2 (en) * | 2008-06-18 | 2012-01-31 | Samsung Electronics Co., Ltd. | System and method for class-specific object segmentation of image data |
JP2010115260A (ja) * | 2008-11-11 | 2010-05-27 | Olympus Corp | 画像処理装置、画像処理プログラムおよび画像処理方法 |
-
2010
- 2010-12-16 JP JP2010280716A patent/JP5576782B2/ja active Active
-
2011
- 2011-12-13 US US13/323,910 patent/US8798344B2/en active Active
- 2011-12-15 CN CN201110420174.XA patent/CN102525381B/zh active Active
- 2011-12-15 EP EP11009867.0A patent/EP2466541B1/en not_active Not-in-force
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5566249A (en) * | 1994-09-20 | 1996-10-15 | Neopath, Inc. | Apparatus for detecting bubbles in coverslip adhesive |
WO1999052072A1 (fr) * | 1998-04-07 | 1999-10-14 | Omron Corporation | Dispositif et procede de traitement d'images, support de stockage de programme pour traitement d'images, et dispositif d'inspection |
CN101065052A (zh) * | 2004-12-27 | 2007-10-31 | 奥林巴斯株式会社 | 医疗图像处理装置和医疗图像处理方法 |
CN101115435A (zh) * | 2005-02-15 | 2008-01-30 | 奥林巴斯株式会社 | 医用图像处理装置、管腔图像处理装置、管腔图像处理方法及用于上述装置和方法的程序 |
CN101150977A (zh) * | 2005-04-13 | 2008-03-26 | 奥林巴斯医疗株式会社 | 图像处理装置以及图像处理方法 |
US20090148014A1 (en) * | 2006-05-26 | 2009-06-11 | Olympus Corporation | Image processing apparatus, image processing method, and image processing program product |
JP2009136655A (ja) * | 2008-02-21 | 2009-06-25 | Olympus Corp | 画像処理装置、画像処理方法および画像処理プログラム |
Non-Patent Citations (1)
Title |
---|
WANG YIGANG ET AL: "a low complexity and high performance real-time algorithm of detecting and tracking circular shape in hole-punching machine", 《PROCEEDINGS OF 2008 3RD INTERNATIONAL CONERENCE ON INTELLIGENT SYSTEM AND KNOWLEDGE ENGINEERING》, 17 November 2008 (2008-11-17), pages 604 - 608, XP031382007, DOI: 10.1109/ISKE.2008.4731002 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103900708A (zh) * | 2012-12-27 | 2014-07-02 | 杭州美盛红外光电技术有限公司 | 红外选择通知装置和红外选择通知方法 |
CN103900713A (zh) * | 2012-12-27 | 2014-07-02 | 杭州美盛红外光电技术有限公司 | 热像检测装置和热像检测方法 |
CN103900698A (zh) * | 2012-12-27 | 2014-07-02 | 杭州美盛红外光电技术有限公司 | 红外识别控制装置和红外识别控制方法 |
CN103900707A (zh) * | 2012-12-27 | 2014-07-02 | 杭州美盛红外光电技术有限公司 | 热像检测更新装置和热像检测更新方法 |
CN103900717A (zh) * | 2012-12-27 | 2014-07-02 | 杭州美盛红外光电技术有限公司 | 热像识别装置和热像识别方法 |
CN103900714A (zh) * | 2012-12-27 | 2014-07-02 | 杭州美盛红外光电技术有限公司 | 热像匹配装置和热像匹配方法 |
WO2014101802A1 (zh) * | 2012-12-27 | 2014-07-03 | Wang Hao | 热像挑选装置和热像挑选方法 |
WO2014101809A1 (zh) * | 2012-12-27 | 2014-07-03 | Wang Hao | 热像检测装置和热像检测方法 |
WO2014101804A1 (zh) * | 2012-12-27 | 2014-07-03 | Wang Hao | 热像选择装置和热像选择方法 |
WO2014101803A1 (zh) * | 2012-12-27 | 2014-07-03 | Wang Hao | 红外选择装置和红外选择方法 |
WO2014101811A1 (zh) * | 2012-12-27 | 2014-07-03 | Wang Hao | 热像匹配装置和热像匹配方法 |
CN104660922A (zh) * | 2013-11-25 | 2015-05-27 | 杭州美盛红外光电技术有限公司 | 热像整理装置、记录装置、整理系统和热像整理方法、记录方法 |
WO2015074625A1 (zh) * | 2013-11-25 | 2015-05-28 | 王浩 | 分析设置装置、处理装置和分析设置方法、处理方法 |
WO2015096824A1 (zh) * | 2013-12-29 | 2015-07-02 | 王浩 | 分析装置和分析方法 |
Also Published As
Publication number | Publication date |
---|---|
EP2466541A3 (en) | 2013-02-06 |
EP2466541A2 (en) | 2012-06-20 |
EP2466541B1 (en) | 2015-12-09 |
JP2012125469A (ja) | 2012-07-05 |
CN102525381B (zh) | 2016-08-03 |
JP5576782B2 (ja) | 2014-08-20 |
US20120155724A1 (en) | 2012-06-21 |
US8798344B2 (en) | 2014-08-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102525381A (zh) | 图像处理装置、图像处理方法及计算机可读取的记录装置 | |
CN102542564B (zh) | 图像处理装置以及图像处理方法 | |
EP3446065B1 (en) | Flight parameter measuring apparatus and flight parameter measuring method | |
US20200356820A1 (en) | Speckle contrast analysis using machine learning for visualizing flow | |
US8837821B2 (en) | Image processing apparatus, image processing method, and computer readable recording medium | |
US20150287192A1 (en) | Image processing device, electronic device, endoscope apparatus, information storage device, and image processing method | |
JP5851160B2 (ja) | 画像処理装置、画像処理装置の作動方法、及び画像処理プログラム | |
Lupascu et al. | Automated detection of optic disc location in retinal images | |
US20080226148A1 (en) | Method of image quality assessment to produce standardized imaging data | |
CN108549874A (zh) | 一种目标检测方法、设备及计算机可读存储介质 | |
CN104363815A (zh) | 图像处理装置、图像处理方法以及图像处理程序 | |
CN104812288A (zh) | 图像处理装置、图像处理方法和图像处理程序 | |
CN102697446A (zh) | 图像处理装置和图像处理方法 | |
CN103945755A (zh) | 图像处理装置、图像处理方法和图像处理程序 | |
CN108601509B (zh) | 图像处理装置、图像处理方法以及记录有程序的介质 | |
CN104732520A (zh) | 一种胸部数字影像的心胸比测量算法及系统 | |
CN104574312A (zh) | 标靶图像求圆心的方法和装置 | |
US20190150848A1 (en) | Image processing apparatus, operation method performed by image processing apparatus and recording medium | |
CN113689412A (zh) | 甲状腺图像处理方法、装置、电子设备及存储介质 | |
CN113038868A (zh) | 医疗图像处理系统 | |
CN106372593A (zh) | 一种基于血管收敛的视盘区定位方法 | |
Yue et al. | The optimal segmentation scale identification using multispectral WorldView-2 images | |
Cudek et al. | Melanocytic skin lesion image classification. Part I: Recognition of skin lesion | |
Chen et al. | Saliency-based bleeding localization for wireless capsule endoscopy diagnosis | |
CN117115926B (zh) | 一种基于实时图像处理的人体动作标准判定方法及装置 |
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 |