CN114730070A - 图像处理方法、图像处理装置和图像处理系统 - Google Patents
图像处理方法、图像处理装置和图像处理系统 Download PDFInfo
- Publication number
- CN114730070A CN114730070A CN202080079048.XA CN202080079048A CN114730070A CN 114730070 A CN114730070 A CN 114730070A CN 202080079048 A CN202080079048 A CN 202080079048A CN 114730070 A CN114730070 A CN 114730070A
- Authority
- CN
- China
- Prior art keywords
- image
- image processing
- filter
- medical image
- focus
- 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.)
- Pending
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 31
- 238000003384 imaging method Methods 0.000 claims description 55
- 230000015572 biosynthetic process Effects 0.000 claims description 21
- 238000003786 synthesis reaction Methods 0.000 claims description 21
- 230000002093 peripheral effect Effects 0.000 claims description 6
- 238000009499 grossing Methods 0.000 claims description 4
- 239000000463 material Substances 0.000 claims description 4
- 230000007170 pathology Effects 0.000 claims description 4
- 239000002131 composite material Substances 0.000 claims description 3
- 230000006870 function Effects 0.000 description 34
- 238000000879 optical micrograph Methods 0.000 description 23
- 238000004364 calculation method Methods 0.000 description 22
- 238000000034 method Methods 0.000 description 21
- 238000010586 diagram Methods 0.000 description 20
- 230000010365 information processing Effects 0.000 description 16
- 210000004027 cell Anatomy 0.000 description 15
- 238000004891 communication Methods 0.000 description 15
- 230000008569 process Effects 0.000 description 11
- 230000003287 optical effect Effects 0.000 description 10
- 210000001519 tissue Anatomy 0.000 description 8
- 230000002194 synthesizing effect Effects 0.000 description 6
- 238000003745 diagnosis Methods 0.000 description 5
- 239000011521 glass Substances 0.000 description 5
- 210000004698 lymphocyte Anatomy 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000001914 filtration Methods 0.000 description 4
- 230000010354 integration Effects 0.000 description 3
- 210000002540 macrophage Anatomy 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 239000004065 semiconductor Substances 0.000 description 3
- 241000590002 Helicobacter pylori Species 0.000 description 2
- 210000003855 cell nucleus Anatomy 0.000 description 2
- 230000001413 cellular effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 229940037467 helicobacter pylori Drugs 0.000 description 2
- 239000004973 liquid crystal related substance Substances 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000008447 perception Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 239000006059 cover glass Substances 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 210000003097 mucus Anatomy 0.000 description 1
- 238000010827 pathological analysis Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000035807 sensation Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B21/00—Microscopes
- G02B21/36—Microscopes arranged for photographic purposes or projection purposes or digital imaging or video purposes including associated control and data processing arrangements
- G02B21/365—Control or image processing arrangements for digital or video microscopes
- G02B21/367—Control or image processing arrangements for digital or video microscopes providing an output produced by processing a plurality of individual source images, e.g. image tiling, montage, composite images, depth sectioning, image comparison
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B21/00—Microscopes
- G02B21/36—Microscopes arranged for photographic purposes or projection purposes or digital imaging or video purposes including associated control and data processing arrangements
- G02B21/365—Control or image processing arrangements for digital or video microscopes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/73—Deblurring; Sharpening
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- 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/10004—Still image; Photographic image
- G06T2207/10012—Stereo images
-
- 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
-
- 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/20—Special algorithmic details
- G06T2207/20004—Adaptive image processing
- G06T2207/20008—Globally adaptive
-
- 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/20—Special algorithmic details
- G06T2207/20172—Image enhancement details
- G06T2207/20182—Noise reduction or smoothing in the temporal domain; Spatio-temporal filtering
-
- 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/20—Special algorithmic details
- G06T2207/20172—Image enhancement details
- G06T2207/20192—Edge enhancement; Edge preservation
-
- 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/20—Special algorithmic details
- G06T2207/20212—Image combination
-
- 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/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
-
- 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
-
- 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
- G06T2207/30024—Cell structures in vitro; Tissue sections in vitro
-
- 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/30168—Image quality inspection
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Multimedia (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Optics & Photonics (AREA)
- Analytical Chemistry (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Quality & Reliability (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Image Processing (AREA)
Abstract
在根据本申请的图像处理方法中,获取由图像捕获装置捕获的医学图像,以及根据医学图像的聚焦程度来确定应用于医学图像的滤波器的强度。
Description
技术领域
本发明涉及一种图像处理方法、一种图像处理设备和一种图像处理系统。
背景技术
存在一种数字显微镜设备和一种图像显示设备,其使用用于观察细胞组织的显微镜设备来捕获细胞组织的图像,将该图像保存为医学图像,并且使用医学图像的图像数据执行病理诊断等。在数字显微镜设备中,为了观察整个样本,通过放大成像系统对载玻片上分割包括样本的区域的小区域进行成像,并且连接每个小区域的多个图像,以创建一个巨大的医学图像。
采用自动聚焦(AF),作为聚焦方法,其中,放大成像系统的物镜的焦点被调整到要成像的细胞组织。例如,已经提出了一种聚焦方法,其中,放大成像系统的物镜的焦点位置在光轴方向上以预定间隔移动,在每个移动位置执行成像,并且当捕获到捕获图像中具有最高对比度的图像时的位置被检测为聚焦位置(例如,参考专利文献1)。这种聚焦方法被称为“对比AF”。
以这种方式捕获的细胞组织图像具有相对较高的聚焦精度,但是与医生或其他观察者通过光学显微镜观察到的光学显微镜图像相比,图像的外观不同。
引用列表
专利文献
专利文献1:JP 2011-197283 A
发明内容
技术问题
如上所述,在数字显微镜设备中,需要以高质量获取细胞组织的图像,但是还没有实现足够的解决方案。
鉴于上述情况,本技术旨在提供一种能够以高质量获取细胞组织图像的数字显微镜设备、成像方法及其程序。
问题的解决方案
根据本申请的图像处理方法包括:获取由成像设备捕获的医学图像;以及根据医学图像的聚焦程度来确定要应用于医学图像的滤波器的强度。
附图说明
图1是示出根据实施例的医学图像的特性的示例的示图;
图2是示出根据实施例的载玻片的特性的示例的示图;
图3是示出根据实施例的图像处理系统的配置示例的示图;
图4是示出根据实施例的指示模糊特征量的函数的示例的示图;
图5是示出根据实施例的指示模糊特征量的倒数的函数的示例的示图;
图6是示出根据实施例的用于计算模糊特征量的像素范围的示例的示图;
图7是示出根据实施例的聚焦位置的示例的示图;
图8是示出根据实施例的用于计算模糊特征量的像素范围的示例的示图;
图9是示出根据实施例的成像设备的配置示例的示图;
图10是示出根据实施例的图像处理设备的配置示例的示图;
图11是示出根据实施例的医学图像存储单元的示例的示图;
图12是示出根据实施例的增强滤波器存储单元的示例的示图;
图13是示出根据实施例的信息处理的示例的流程图;
图14是示出根据实施例的增强滤波器的效果的示例的示图;
图15是示出根据实施例的增强滤波器的效果的示例的示图;
图16是示出用于实现图像处理设备的功能的计算机的示例的硬件配置图。
具体实施方式
下面将参考附图详细描述用于实现根据本申请的图像处理方法、图像处理设备和图像处理系统的模式(以下称为“实施例”)。注意,根据本申请的图像处理方法、图像处理设备和图像处理系统不限于实施例。在以下每个实施例中,相同的部分由相同的附图标记表示,并且将省略多余的描述。
将按照以下项目顺序描述本公开。
1.光学显微镜的特点
2.图像处理系统的配置
3.信息处理示例
4.处理变化
4-1.抽头范围
4-2.层数
4-3.滤波器类型
4-4.构成具有多个像素的块
4-5.被摄体
4-6.成像设备
4-7.聚焦位置
4-8.计算特征量的方法
4-9.设备集成
5.成像设备的配置
6.图像处理设备的配置
7.信息处理流程
8.修改
9.硬件配置
10.其他
实施例
1.光学显微镜的特点
图1示出了光学显微镜图像的特点。光学显微镜图像与通过数字显微镜设备观察到的图像具有不同的特点。如图1所示,有几个与光学显微镜图像的图像质量相关的指标。将使用图1所示的指数来描述光学显微镜图像的特点。例如,存在指示图像质量的亮度的指数(在本实施例中,有时表示为“闪光感”)。例如,存在指示图像质量边缘程度的指数(在本实施例中,有时被表示为“清晰感”)。在图1中,指示光学显微镜图像的特点的指数例如是“立体/平面”、“透明/暗淡”或“清晰/模糊”。
根据实施例的模糊(不清晰度)指示图像不清晰的状态。具体而言,根据实施例的模糊指示图像失焦并且在景深范围之外不清晰的状态。根据实施例的聚焦指示景深范围内的聚焦状态。根据该实施例的聚焦程度是通过对调节焦点的程度进行评分而获得的值。
如图1所示,光学显微镜图像是立体的。根据该实施例的术语“立体”表示与模糊和聚焦之间的视觉反差(对比度)相关的质量。
光学显微镜图像有一种透明感。根据实施例的透明感指示与噪声相关的质量。根据实施例的噪声是除了被摄体之外的不必要的信息。具体地,由于光学显微镜图像没有被数字化,所以没有噪声的增强,并且图像具有透明感。
光学显微镜图像有一种闪光感。根据实施例的闪光感表示当光被施加到被摄体时,由散射光产生的干涉条纹引起的与亮度相关的质量。具体而言,光学显微镜图像发出的光比通过干涉条纹施加到物体上的光更亮,因此具有闪光感。
光学显微镜图像有一种清晰感。根据实施例的清晰度感觉指示与清晰度相关的质量。
光学显微镜图像是立体的、明亮的和高度清晰的,因此具有识别目标的高能力(以下适当地称为“目标识别性能”)。
光学显微镜图像是立体的、明亮的和高度清晰的,因此具有识别目标的高能力(以下适当地称为“目标识别性能”)。
下面将描述将数字显微镜设备获取的图像近似为光学显微镜图像的方法。样本放在载玻片上。在样本中,细胞等在Z轴方向(以下适当地称为“Z方向”)上分布,Z轴方向指示载玻片的厚度方向,并且由数字显微镜设备获取的医学图像具有聚焦区域和非聚焦区域的混合。例如,当对整个医学图像应用高范围增强滤波器时,低频部分(例如,噪声)也增强,并且未聚焦区域的清晰度也提高。图像高频成分多,目标识别性能变差。当滤波器的强度降低时,待增强区域的增强也减弱。
近似光学显微镜图像的图像是具有图1所示的特点的图像,并且是由具有不同聚焦程度的图像形成的图像。下面将描述接近光学显微镜图像的效果。
例如,将数字显微镜设备获取的医学图像近似为光学显微镜图像,允许容易地看到细胞的结构。因此,可以促进用于诊断的利用。例如,将医学图像近似为光学显微镜图像,允许容易地辨别细胞的位置。因此,提高了病理学家的诊断速度,并且可以减少疲劳。例如,将医学图像近似为光学显微图像增加了重叠细胞的可见性。因此,可以进行疾病类型的诊断,其中,重叠细胞的识别是重要的。例如,将医学图像近似为光学显微图像,允许病理学家容易地适应使用医学图像的诊断。例如,接近光学显微镜图像可以防止诸如幽门螺杆菌等小物体隐藏在噪声中。例如,由于区域限制增强,接近光学显微镜图像可以确保高压缩效率。
图2示出了用于成像的载玻片的示例。图2是从垂直方向看到的用于成像的载玻片的示图。如图2所示,Z方向是载玻片厚度的方向。换言之,Z方向是物体厚度的方向。由于被摄体是从载玻片上方成像的,所以Z方向是垂直于医学图像的方向。Z方向是成像时的光轴方向。图2示出了被摄体被放置在载玻片上并用盖玻片覆盖的情况。图2示出了被摄体是组织切片的情况。图2示出了对组织中包含的诸如淋巴细胞和巨噬细胞等样本进行成像的情况。
如图2所示,淋巴细胞是厚的,因为淋巴细胞重叠。淋巴细胞下面是巨噬细胞。由于要成像的景深较浅,因此不可能对例如焦点处的整个淋巴细胞或巨噬细胞进行成像。下面将描述图像处理设备100的信息处理,作为校正由于数字化导致的可见度降低的处理。
2.图像处理系统的配置
将参考图3描述图像处理系统1的配置。图3是示出根据实施例的图像处理系统的配置示例的示图。如图3所示,图像处理系统1包括成像设备10和图像处理设备100。成像设备10和图像处理设备100经由预定通信网络(网络N)连接,以便能够以有线或无线方式彼此通信。图3是示出根据实施例的图像处理系统的配置示例的示图。图3所示的图像处理系统1可以包括多个成像设备10和多个图像处理设备100。
成像设备10是诸如显微镜等成像设备,并且用于对样本成像。
图像处理设备100用于根据被摄体的聚焦程度来确定与滤波器相关的信息。图像处理设备100例如是诸如PC或工作站(WS)等信息处理设备,并且基于例如经由网络N从成像设备10发送的信息来执行处理
3.信息处理示例
下面将描述图像处理设备100根据被摄体的聚焦程度来确定滤波器的强度的情况。以下将描述样本(例如,细胞),作为被摄体的示例。
根据实施例的滤波器是用于提高医学图像的图像质量的滤波器。根据实施例的滤波器被应用于通过对被摄体成像而获得的医学图像。根据实施例的滤波器可以是任何类型的滤波器。换言之,假设对由根据实施例的滤波器增强的区域没有限制。例如,根据实施例的滤波器包括滤波器,例如,高范围增强滤波器、中范围增强滤波器、低范围增强滤波器、负增强滤波器即平滑滤波器。
图像处理设备100使用模糊函数(以下适当地称为“模糊确定函数”或“模糊量确定函数”)来计算聚焦程度。下面将描述图像处理设备100生成模糊函数的过程。
通过洛伦兹函数近似相邻差值的平方和来计算倒数,生成模糊函数。根据实施例的近似是曲线图的拟合(曲线拟合)。表达式(1)示出了根据实施例的相邻差值的平方和。
图4示出了相邻差值平方和的曲线图GR1。图4示出了曲线图GR1,其中,预定像素和相对于医学图像中的所有像素具有预定关系的像素之间的差的总值用作特征量,并且Z方向上的值(以下适当地称为“Z值”)被绘制为变量。这个特征量被适当地称为“模糊特征量”。具有预定关系的像素是与预定像素相邻的像素。在这种情况下,图4示出了一个曲线图,其中,为医学图像中的所有像素绘制了预定像素和相邻像素之间的差值的总值。在图4中,曲线图GR1的水平轴(X轴)是载玻片的Z值。曲线图GR1的垂直轴(Y轴)表示特征量。在图4中,相邻差值的平方和的输出值由s(z)表示。在图4中,执行绘图,使得特征量变得最大的Z值变为0。换言之,执行绘图,使得相邻差值的平方和的输出值的最大值是s(0)。图4示出了特征量越大,聚焦程度越高。
图4可以示出相邻差值的平方和由洛伦兹函数近似的曲线图。表达式(2)示出了根据实施例的洛伦兹函数。
在这种情况下,图4中的曲线图GR1的Y轴可以是f(z)。换言之,图4中的曲线图GR1的Y轴可以由f(z)表示为洛伦兹函数的输出值。
图5示出了曲线图GR2,该曲线图是洛伦兹函数f(z)的倒数。图5示出了通过绘制洛伦兹函数f(z)的倒数而获得的曲线图GR2。曲线图GR2是二次曲线。表达式(3)示出了根据实施例的二次曲线。
曲线图GR2中的顶点表示聚焦位置。根据实施例的聚焦位置是聚焦程度变得最大的Z值。
图像处理设备100获取关于预定聚焦图像的图像信息。图像处理设备100使用预定的聚焦图像作为参考,获取预定数量的图像信息,该图像信息包括预定的聚焦图像和聚焦在具有不同Z值的位置上的图像。例如,图像处理设备100获取聚焦在与预定聚焦图像的Z值在Z方向上相差几微米的位置上的图像信息。作为具体示例,图像处理设备100获取总共三条图像信息,包括预定聚焦图像和聚焦在与预定聚焦图像的Z值在Z方向上相差几微米的位置上的两个图像。
图像处理设备100通过使用相邻差值的平方和来计算获取的图像信息的模糊特征量。具体地,图像处理设备100通过将相邻差值的平方和近似为洛伦兹函数来计算获取的图像信息的模糊特征量。
图像处理设备100计算所计算出的模糊特征量的倒数。图像处理设备100基于计算出的模糊特征量的倒数执行对二次曲线的拟合。图像处理设备100基于拟合的二次曲线来估计聚焦位置。具体地,图像处理设备100估计拟合的二次曲线的顶点,作为整个图像的聚焦位置。
图像处理设备100基于估计的聚焦位置计算获取的图像信息的聚焦程度。具体地,图像处理设备100基于估计的聚焦位置和用于模糊函数的Z值之间的差来计算获取的图像信息的聚焦程度。
图像处理设备100根据计算出的聚焦程度来确定滤波器的强度。
下面将参考图6描述图像处理设备100计算模糊特征量的过程。
图像处理设备100使用医学图像中的预定像素作为参考来排列以预定像素为中心的抽头,并且通过使用相邻差值的平方和来计算模糊特征量。根据实施例的抽头指示以预定像素为中心的像素范围。换言之,根据实施例的抽头指示感兴趣像素的外围像素的范围。例如,根据实施例的抽头指示应用了滤波器的感兴趣像素的外围像素的范围。例如,3×3抽头指示总共九个像素的范围,其中,图像的垂直和水平像素各为三个像素。图像处理设备100在以不同Z值获取的每个图像上排列相同的抽头,并且计算相邻差值的平方和,作为模糊特征量。图6示出了图像处理设备100排列3×3抽头TA1的示例。在图6中,图像处理设备100使用S0作为预定像素来计算以S0为中心的模糊特征量。在这种情况下,图像处理设备100基于表达式(4)计算模糊特征量。
模糊特征量(3x3的情况下=(s0-s1)2+(s0-s2)2+(s0-s3)2+(s0-s4)2+(s0-s5)2+(s0-s6)2+(s0-s7)2+(s0-s8)2
(4)
具体地,当图像处理设备100获取总共三条图像信息时,该图像信息包括预定聚焦图像和聚焦在与预定聚焦图像的Z值在Z方向上不同的位置上的两个图像,图像处理设备100从三条中的每一条排列3×3抽头,并且通过使用距中心的相邻差值的平方和来计算模糊特征量。在这种情况下,图像处理设备100基于表达式(5)计算模糊特征量。
模糊特征量F2(上)=(s0[2]-s1[2])2+(s0[2]-s2[2])2+(s0[2]-s3[2])2+(s0[2]-s4[2])2+(s0[2]-s5[2])2+(s0[2]-s6[2])2+(s0[2]-s7[2])2+(s0[2]-s8[2])2
模糊特征量F2(中)=(s0[1]-s1[1])2+(s0[1]-s2[1])2+(s0[1]-s3[1])2+(s0[1]-s4[1])2+(s0[1]-s5[1])2+(s0[1]-s6[1])2+(s0[1]-s7[1])2+(s0[1]-s8[1])2
模糊特征量F2(下)=(s0[0]-s1[0])2+(s0[0]-s2[0])2+(s0[0]-s3[0])2+(s0[0]-s4[0])2+(s0[0]-s5[0])2+(s0[0]-s6[0])2+(s0[0]-s7[0])2+(s0[0]-s8[0])2
(5)
图像处理设备100通过使用所获取的图像信息中位于预定Z值的层上的图像信息来计算模糊特征量。例如,图像处理设备100通过使用所获取的三条图像信息中上层的图像信息来计算模糊特征量F2。例如,图像处理设备100通过使用所获取的三条图像信息中关于中间层的图像信息来计算模糊特征量F1。例如,图像处理设备100通过使用所获取的三条图像信息中下层的图像信息来计算模糊特征量F0。
图7示出了曲线图GR3,其中,由图像处理设备100计算的三个值的倒数被拟合到二次曲线。图像处理设备100基于拟合来估计聚焦位置。图像处理设备100估计拟合的二次曲线的顶点,作为聚焦位置。在图7中,Zc是聚焦位置。在图7中,Zc是由模糊特征量估计的聚焦位置。在图7中,Z1是整个图像中最聚焦的位置。在图7中,最靠近Zc的Z1是整个图像中最聚焦的位置。图像处理设备100可以向Z1应用强度等于通过应用将被应用于Zc的滤波器而获得的结果的滤波器。当Z1在被视为等于Zc的预定范围内时,图像处理设备100可以应用强度等于Zc至Z1的滤波器。
图像处理设备100根据距聚焦位置的距离计算聚焦程度。在图7中,图像处理设备100根据距Zc的距离计算聚焦程度。
图像处理设备100根据计算出的聚焦程度来确定要应用的滤波器的强度。图像处理设备100根据确定的强度应用滤波器。具体而言,图像处理设备100确定每个像素的滤波强度,并对像素执行滤波处理。因此,图像处理设备100通过对每个像素重复处理来执行最适合于图像的滤波处理。
因此,图像处理设备100可以根据聚焦程度改变滤波器的强度。因此,图像处理设备100可以仅对图像的聚焦区域应用强度等于100%的滤波器。对整个图像应用强度为100%的滤波器也会增加未聚焦区域和噪声的清晰度。然而,由于图像处理设备100仅将滤波器应用于聚焦区域,所以仅在聚焦区域中提高了清晰度。由于未聚焦区域的清晰度没有增强,所以图像的深度感没有降低,并且容易识别细胞的重叠。由于噪声的清晰度没有增强,所以可以抑制诸如幽门螺杆菌等精细被摄体被隐藏在噪声中。因此,图像处理设备100可以调整局部对比度。因此,图像处理设备100可以校正对比度被成像光学系统降低的图像。因此,图像处理设备100可以部分地调整对比度。因此,即使当粘液或其他物质在细胞核中混合并且基因信息不存在于细胞核的中心时,图像处理设备100也可以确认情况,并且提高诊断的准确性。
4.处理变化
4-1.抽头范围
尽管在上述示例中,图像处理设备100选择以预定像素为中心的正方形抽头,但是所选择的抽头范围不限于诸如3x3等正方形,并且可以选择任何范围,作为抽头。例如,如图8所示,图像处理设备100可以选择以S0为中心的十字形抽头范围。所选范围内的垂直和水平像素的数量可以不限于三个像素。图像处理设备100可以选择11×11的范围,作为抽头,其中,垂直和水平像素中的每一个都是11个像素。
4-2.层数
尽管在上述示例中,图像处理设备100通过使用利用具有不同Z方向的三条图像信息计算的三个值的倒数来估计聚焦位置,但是假设用于估计聚焦位置的图像信息的数量不受限制,只要该数量为三条或更多。例如,图像处理设备100可以通过使用利用具有不同Z方向的四条图像信息计算的四个值的倒数来估计聚焦位置。
4-3.滤波器类型
尽管在上述示例中,图像处理设备100根据聚焦程度来确定滤波器的强度,但是图像处理设备100可以根据聚焦程度来确定滤波器的类型。例如,图像处理设备100根据估计的聚焦位置来确定滤波器的类型。例如,如果聚焦程度满足预定条件,则图像处理设备100确定相应特定滤波器的类型。例如,如果聚焦程度满足预定条件,则图像处理设备100确定用于增强相应预定区域的滤波器的类型。例如,如果聚焦程度大于或等于预定阈值,则图像处理设备100确定应用中范围增强滤波器。例如,如果聚焦程度小于预定阈值,则图像处理设备100确定应用高范围增强滤波器。或者,图像处理设备100确定应用负增强滤波器,即平滑滤波器。
图像处理设备100可以根据聚焦程度来确定滤波器的强度和类型。图像处理设备100可以根据聚焦程度同时确定滤波器的强度和类型。在这种情况下,图像处理设备100可以将确定强度和类型的滤波器应用于医学图像。因此,图像处理设备100可以根据聚焦程度应用具有最佳强度的最佳滤波器。
4-4.构成具有多个像素的块
尽管在上述示例中,图像处理设备100确定每个像素的滤波器的强度和类型,但是图像处理设备100可以构成具有多个像素的块,并且确定每个块的滤波器的强度和类型。在这种情况下,图像处理设备100确定用于构成块的多个像素。图像处理设备100可以以任何方式确定用于构成块的多个像素。例如,图像处理设备100构成具有多个相邻像素的块。例如,图像处理设备100构成具有多个预定像素的块。
4-5.被摄体
尽管在上述示例中,图像处理设备100使用样本(例如,细胞)作为被摄体来获取医学图像,但是图像处理设备可以使用生物体或从生物体收集的任何物体,作为被摄体,来获取图像。例如,图像处理设备100使用医学领域中与活体、生物体、材料、病理等相关的样本,作为被摄体,来获取图像。
4-6.成像设备
假设根据实施例的成像设备可以是能够对被摄体成像的任何设备。例如,根据实施例的成像设备是显微镜。如果根据实施例的成像设备是显微镜,则可以使用任何显微镜。
4-7.聚焦位置
尽管在上述示例中,图像处理设备100估计拟合二次曲线的顶点,作为整个图像的聚焦位置,但是图像处理设备可以估计从拟合二次曲线的顶点起预定范围内的位置,作为整个图像的聚焦位置。
4-8.计算特征量的方法
尽管在上述实施例中,具有预定关系的像素是与预定像素相邻的像素,但是像素不限于该示例。换言之,具有预定关系的像素不必限于相邻像素。例如,具有预定关系的像素可以是每隔一个像素或每两个像素的像素。例如,图像处理设备100可以通过使用预定像素和每隔一个像素的像素之间的差的总值,作为特征量,来计算聚焦程度。在这种情况下,图像处理设备100根据通过使用预定像素和每隔一个像素的像素之间的差的总值作为特征量而计算的聚焦程度,来确定要应用于医学图像的滤波器的强度。
4-9.设备集成
尽管在上述实施例中成像设备10和图像处理设备100是单独的设备,但是成像设备10和图像处理设备100可以集成。例如,图像处理设备100的功能可以在控制成像设备10的操作的计算机中实现,或者可以在成像设备10的外壳内提供的任何计算机中实现。图像处理设备100的功能可以下载到控制成像设备10的操作的计算机,或者可以下载到设置在成像设备10的外壳内的任何计算机。因此,可以出售具有图像处理设备100的功能的成像设备10。
5.成像设备的配置
然后将参考图9描述根据实施例的成像设备10的配置。图9是示出根据实施例的成像设备10的配置示例的示图。如图9所示,成像设备10包括通信单元11、存储单元12和控制单元13。成像设备10可以具有用于从成像设备10的管理员接收各种操作的输入单元(例如,键盘或鼠标)以及用于显示各种类型的信息的显示单元(例如,液晶显示器)。
通信单元11
通信单元11例如由网络接口卡(NIC)实现。通信单元11以有线或无线方式连接到预定网络N,并且经由预定网络N向例如图像处理设备100发送信息和从其接收信息。
存储单元12
存储单元12例如由诸如RAM或闪存等半导体存储元件或者诸如硬盘或光盘等存储装置来实现。存储单元12存储与医学图像相关的信息。具体地,存储单元12存储与成像被摄体的医学图像相关的信息。
控制单元13
控制单元13是控制器,并且例如由CPU或MPU通过使用RAM作为工作区域来执行存储在成像设备10内部的存储装置中的各种程序来实现。控制单元13是控制器,并且例如由诸如ASIC或FPGA等集成电路来实现。
如图9所示,控制单元13包括成像单元141,并且实现或执行将在下面描述的信息处理的操作。控制单元13的内部配置不限于图9所示的配置,并且可以使用其他配置,只要执行下面将要描述的信息处理。
成像单元141
成像单元141对各种类型的信息进行成像。成像单元141在载玻片上对被摄体成像。成像单元141获取各种类型的信息。成像单元141获取成像的医学图像。
6.图像处理设备的配置
然后将参考图10描述根据实施例的图像处理设备100的配置。图10是示出根据实施例的图像处理设备100的配置示例的示图。如图10所示,图像处理设备100包括通信单元110、存储单元120和控制单元130。图像处理设备100可以具有用于从图像处理设备100的管理员接收各种操作的输入单元(例如,键盘或鼠标)以及用于显示各种类型的信息的显示单元(例如,液晶显示器)。
通信单元110
通信单元110例如由NIC实现。通信单元110以有线或无线方式连接到网络N,并且经由网络N向例如成像设备10发送信息和从成像设备接收信息。
存储单元120
存储单元120例如由诸如RAM或闪存等半导体存储元件或者诸如硬盘或光盘等存储装置来实现。如图10所示,存储单元120包括医学图像存储单元121和增强滤波器存储单元122。
医学图像存储单元121存储与医学图像相关的信息。图11示出了根据实施例的医学图像存储单元121的示例。如图11所示,医学图像存储单元121具有诸如“医学图像ID”和“医学图像”等项目。
“医学图像ID”表示用于识别医学图像的识别信息。“医学图像”表示通过对被摄体进行成像而获得的医学图像。例如,“医学图像”表示由成像设备10捕获的医学图像。尽管在图11所示的示例中,诸如“医学图像#11”和“医学图像#12”等概念信息存储在“医学图像”中,但是例如,实际上存储了静止图像或运动图像。静止图像或运动图像可以保存在不同于图像处理设备的服务器或云中,并且“医学图像”可以存储例如医学图像的内容所在的统一资源定位符(URL)或者指示存储位置的文件路径名称。
增强滤波器存储单元122存储与要应用于医学图像的增强滤波器相关的信息。图12示出了根据实施例的增强滤波器存储单元122的示例。如图12所示,增强滤波器存储单元122具有诸如“增强滤波器ID”和“增强滤波器”等项目。
“增强滤波器ID”表示用于识别增强滤波器的识别信息。“增强滤波器”指示与增强滤波器相关的信息。尽管在图12所示的示例中,诸如“增强滤波器#11”和“增强滤波器#12”等概念信息被存储在“增强滤波器”中,但是例如,实际存储增强滤波器的强度和类型以及应用了增强滤波器的像素范围。
控制单元130
控制单元130是控制器,并且由例如CPU或MPU通过使用RAM作为工作区域来执行存储在图像处理设备100内部的存储装置中的各种程序来实现。控制单元130是控制器,并且例如由诸如ASIC或FPGA等集成电路来实现。
如图10所示,控制单元130包括获取单元131、计算单元132和确定单元133,并且实现或执行将在下面描述的信息处理的操作。控制单元130的内部配置不限于图10所示的配置,并且可以使用其他配置,只要执行下面将要描述的信息处理。
获取单元131
获取单元131获取各种类型的信息。获取单元131从外部信息处理设备获取各种类型的信息。获取单元131从诸如成像设备10等另一信息处理设备获取各种类型的信息。
获取单元131从存储单元120获取各种类型的信息。获取单元131从医学图像存储单元121和增强滤波器存储单元122获取各种类型的信息。
获取单元131将获取的各种类型的信息存储在存储单元120中。获取单元131将各种类型的信息存储在医学图像存储单元121和增强滤波器存储单元122中。
获取单元131获取由其他功能配置计算和确定的各种类型的信息。
获取单元131获取被摄体的医学图像。例如,获取单元131获取由成像设备10捕获的被摄体的医学图像。例如,获取单元131获取医学领域中与活体、生物体、材料或病理学相关的被摄体的医学图像。例如,获取单元131获取由显微镜捕获的医学图像。
计算单元132
计算单元132计算各种类型的信息。计算单元132从存储单元120计算各种类型的信息。计算单元132计算来自医学图像存储单元121和增强滤波器存储单元122的各种类型的信息。
计算单元132将计算出的各种类型的信息存储在存储单元120中。计算单元132将各种类型的信息存储在医学图像存储单元121和增强滤波器存储单元122中。
计算单元132计算由其他功能配置获取和确定的各种类型的信息。计算单元132基于由其他功能配置获取和确定的各种类型的信息来计算各种类型的信息。
计算单元132通过使用获取的图像信息来计算相邻差值的平方和。计算单元132计算以预定像素为中心的所确定的抽头范围中的相邻差值的平方和。例如,计算单元132基于通过对医学图像中的预定像素和相邻像素之间的差求和而获得的特征量来计算相邻差值的平方和。
计算单元132根据计算出的相邻差值的平方和计算洛伦兹函数。计算单元132根据计算出的相邻差值的平方和计算洛伦兹函数的近似值。计算单元132将计算的相邻差值的平方和近似为洛伦兹函数。
计算单元132根据计算的洛伦兹函数计算洛伦兹函数的倒数。计算单元132根据洛伦兹函数的倒数估计聚焦位置。计算单元132估计洛伦兹函数的倒数的顶点,作为聚焦位置。
计算单元132计算聚焦程度。计算单元132通过使用估计的聚焦位置来计算聚焦程度。计算单元132根据距估计的聚焦位置的距离计算聚焦程度。计算单元132根据对应于计算的洛伦兹函数的倒数的Z值和估计的聚焦位置之间的距离来计算聚焦程度。
确定单元133
确定单元133确定各种类型的信息。确定单元133从存储单元120确定各种类型的信息。确定单元133从医学图像存储单元121和增强滤波器存储单元122确定各种类型的信息。
确定单元133将确定的各种类型的信息存储在存储单元120中。确定单元133将各种类型的信息存储在医学图像存储单元121和增强滤波器存储单元122中。
确定单元133确定通过其他功能配置获取和计算的各种类型的信息。确定单元133基于由其他功能配置获取和计算的各种类型的信息来确定各种类型的信息。
确定单元133选择应用了滤波器的像素。确定单元133选择应用了滤波器的第一像素。例如,确定单元133在对预定像素应用滤波器之后选择下一个像素。例如,确定单元133基于预定算法选择应用了滤波器的像素。
确定单元133根据计算出的聚焦程度来确定要应用于医学图像的滤波器的强度。确定单元133根据计算出的聚焦程度,根据具有不同强度的多个滤波器来确定指示医学图像的合成比率的合成率。确定单元133根据计算出的聚焦程度选择性地确定要应用于医学图像的多个滤波器。
确定单元133根据计算出的聚焦程度来确定要应用于医学图像的滤波器的类型。确定单元133根据计算出的聚焦程度,根据具有不同类型的多个滤波器来确定指示医学图像的合成比率的合成率。
确定单元133将确定强度的滤波器应用于医学图像。确定单元133将确定类型的滤波器应用于医学图像。
确定单元133根据计算出的聚焦程度,将针对具有不同Z值的多个医学图像确定的强度的滤波器应用于每个医学图像。然后,确定单元133根据计算出的聚焦程度来确定合成率,该合成率指示应用了滤波器的每个医学图像的合成比率。然后,确定单元133基于根据计算出的聚焦程度确定的合成率来合成每个医学图像。
确定单元133确定是否已经对所有像素应用了滤波器。确定单元133确定是否已经对医学图像中的所有像素执行了滤波处理。确定单元133确定是否已经对医学图像中包括在预定范围内的所有像素执行了滤波处理。
7.信息处理流程
然后将参考图13描述根据实施例的图像处理系统1的信息处理过程。图13是示出根据实施例的图像处理系统1的信息处理过程的流程图。
如图13所示,图像处理设备100获取医学图像(步骤S101)。图像处理设备100选择用于执行滤波处理的预定像素(步骤S102)。图像处理设备100使用模糊函数来估计聚焦程度(步骤S103)。图像处理设备100根据估计的聚焦程度来确定滤波器的类型和强度(步骤S104)。图像处理设备100应用确定类型和强度的滤波器(步骤S105)。图像处理设备100确定是否已经对所有像素执行了滤波处理(步骤S106)。如果确定没有对所有像素执行滤波处理(步骤S106否),则图像处理设备100返回到步骤S102的处理。例如,图像处理设备100返回到步骤S102的处理,并且选择除了先前选择的像素之外的像素。另一方面,如果确定已经对所有像素执行了滤波处理(步骤S106是),图像处理设备100结束信息处理。
8.修改
根据上述实施例的图像处理系统1可以以不同于上述实施例的各种不同形式来实现。因此,下面将描述图像处理系统1的另一实施例。
上述实施例已经描述了对每个像素具有不同Z值的多个医学图像应用最佳滤波器的处理。参考图14和图15,将给出能够获得等同于对每个像素具有不同Z值的多个医学图像应用最佳滤波器的处理的效果的处理的描述。
图14示出了图像处理设备100对获取的医学图像应用一个增强滤波器的情况。在图14中,图像KG1表示滤波应用之前的图像,图像KG2表示滤波应用之后的图像。在这种情况下,生成具有整个图像的改善的对比度和分辨率的图像,包括除了感兴趣的目标之外的目标。
图像处理设备100可以对获取的医学图像应用多个增强滤波器。图像处理设备100可以生成通过根据聚焦程度应用具有不同强度的多个滤波器而获得的合成图像。
图15示出了图像处理设备100对通过合成获取的医学图像应用多个增强滤波器而获得的图像的情况。在图15中,图像KG11表示通过应用第一滤波器而获得的图像,图像KG22表示通过应用强度不同于第一滤波器的第二滤波器而获得的图像。如图15所示,图像KG11比图像KG22更模糊。图像KG33表示通过合成图像KG11和图像KG22而获得的合成图像。在图15中,目标SJ11是指示图像的感兴趣区域的第一区域,目标SJ22是并非感兴趣区域的第二区域。例如,目标SJ22是医学图像中噪声的一部分。
图像处理设备100可以确定指示医学图像的合成比率的合成率。例如,图像处理设备100可以确定第一区域和第二区域的合成率。在图15中,图像处理设备100相对于目标SJ11的区域以1:9的合成率合成图像KG11和图像KG12。合成率与目标区域SJ11的比率可以是任何比率,只要图像KG22的比率高于图像KG11的比率。在图15中,图像处理设备100相对于目标SJ22的区域以9:1的合成率合成图像KG11和图像KJ22。合成率与目标区域SJ22的比率可以是任何比率,只要图像KG11的比率高于图像KJ22的比率。
尽管在图15中合成了通过应用增强滤波器而获得的两个图像,但是假设图像处理设备100要合成的图像的数量不受限制。图像处理设备100可以选择性地确定要合成的图像的数量,或者可以选择性地确定要合成的图像。例如,当在通过应用第一滤波器而获得的图像中存在无法辨别的区域时,图像处理设备100选择性地合成通过应用其他滤波器而获得的多个图像和通过应用第一滤波器而获得的图像,从而生成与通过仅应用第一滤波器而获得的图像相比具有改善的可视性和可识别性的图像。
因此,图像处理设备100可以通过合成具有不同强度的增强滤波器的处理图像来调整局部对比度。因此,图像处理设备100可以通过自适应地应用最佳滤波器来调整全局和局部对比度,并且如同光学显微镜一样提高可视性。
图15示出了对一个图像应用多个不同的滤波器而获得的图像进行合成的示例。如上所述,图像处理设备100可以通过合成通过对一个图像应用多个不同的滤波器而获得的图像来生成合成图像。如在上述实施例中,图像处理设备100可以通过根据聚焦程度合成通过根据相应的聚焦程度对具有不同Z值的多个医学图像应用最佳滤波器而获得的图像来生成合成图像。
9.硬件配置
根据上述实施例的成像设备10和图像处理设备100由例如具有图16所示配置的计算机1000实现。图16是示出实现成像设备10和图像处理设备100的功能的计算机的示例的硬件配置图。计算机1000包括CPU 1100、RAM 1200、ROM 1300、HDD 1400、通信接口(I/F)1500、输入/输出接口(I/F)1600和媒体接口(I/F)1700。
CPU 1100基于存储在ROM 1300或HDD 1400中的程序进行操作,并控制每个单元。ROM 1300存储例如当激活计算机1000时由CPU 1100执行的引导程序以及依赖于计算机1000的硬件的程序。
例如,HDD 1400存储要由CPU 1100执行的程序和要由程序使用的数据。通信接口1500经由预定的通信网络从其他装置接收数据,将数据发送到CPU 1100,并且经由预定的通信网络将由CPU 1100生成的数据发送到其他装置。
CPU 1100经由输入/输出接口1600控制诸如显示器或打印机等输出装置和诸如键盘或鼠标等输入装置。CPU 1100经由输入/输出接口1600从输入装置获取数据。CPU 1100经由输入/输出接口1600将生成的数据输出到输出装置。
媒体接口1700读取存储在记录介质1800中的程序或数据,并经由RAM 1200将读取的程序或数据提供给CPU 1100。CPU 1100经由媒体接口1700从记录介质1800将程序加载到RAM 1200上,并执行加载的程序。记录介质1800例如是诸如数字多功能盘(DVD)或相变可重写盘(PD)等光学记录介质、诸如磁光盘(MO)等磁光记录介质、磁带介质、磁记录介质或半导体存储器。
例如,当计算机1000用作根据实施例的成像设备10和图像处理设备100时,计算机1000的CPU 1100执行加载到RAM 1200上的程序,以实现控制单元13和130的功能。计算机1000的CPU 1100从记录介质1800读取这些程序并执行这些程序,但是作为另一示例,可以经由预定的通信网络从另一装置获取这些程序。
10.其他
在上述实施例和修改中描述的过程中,可以手动执行被描述为自动执行的所有或一些过程,或者可以通过已知方法自动执行被描述为手动执行的所有或一些过程。此外,除非另有说明,否则在以上描述和附图中示出的处理过程、具体名称以及包括各种类型的数据和参数的信息可以任意改变。例如,每个图中所示的各种类型的信息不限于所示的信息。
附图中示出的每个装置的每个组件在功能上是概念性的,并且不一定需要如附图中示出的那样进行物理配置。换言之,装置的分布和集成的具体形式不限于附图中所示的形式,并且可以通过根据例如各种负载或使用条件在功能上或物理上分布或集成在任意单元中来配置所有或一些装置。
上述实施例和修改可以在与处理内容不矛盾的范围内适当合成。
虽然已经参考附图详细描述了本申请的一些实施例,但是这些实施例仅仅是示例,并且本发明可以以其他形式实施,其中,基于本领域技术人员的知识进行各种修改和改进,包括在本发明的公开中描述的方面。
例如,上述“部分(部分、模块和单元)”可以用“装置”或“电路”代替。例如,获取单元可以用获取装置或获取电路代替。
注意,本技术也可以具有以下配置。
(1)一种图像处理方法,包括:由计算机
获取由成像设备捕获的医学图像;以及
根据医学图像的聚焦程度确定要应用于医学图像的滤波器的强度。
(2)根据(1)所述的图像处理方法,包括
根据由指示医学图像的模糊程度的模糊函数计算的聚焦程度,确定要应用于医学图像的滤波器的强度。
(3)根据(1)或(2)所述的图像处理方法,包括
根据模糊函数根据被摄体在垂直于医学图像的方向上的厚度方向计算的聚焦程度,确定要应用于医学图像的滤波器的强度。
(4)根据(1)至(3)中任一项所述的图像处理方法,包括
根据通过对医学图像中的预定像素和外围像素之间的差求和而获得的特征量计算的聚焦程度,确定要应用于医学图像的滤波器的强度。
(5)根据(1)至(4)中任一项所述的图像处理方法,包括
基于根据模糊函数估计的聚焦程度和通过对医学图像中的预定像素和外围像素之间的差求和而获得的特征量计算的聚焦程度,确定要应用于医学图像的滤波器的强度。
(6)根据(1)至(5)中任一项所述的图像处理方法,包括
根据被摄体的聚焦程度,确定合成率,该合成率指示通过将多个不同的滤波器应用于医学图像而生成的合成图像的合成比率。
(7)根据(1)至(6)中任一项所述的图像处理方法,包括
确定通过将具有不同强度的多个滤波器应用于预定医学图像而生成的合成图像的合成率。
(8)根据(1)至(7)中任一项所述的图像处理方法,包括
根据被摄体对具有不同Z值的多个医学图像的聚焦程度,确定通过应用滤波器而生成的合成图像的合成率。
(9)根据(1)至(8)中任一项所述的图像处理方法,包括
根据被摄体的聚焦程度,选择性地确定要应用于医学图像以生成合成图像的多个滤波器。
(10)根据(1)至(9)中任一项所述的图像处理方法,包括
根据被摄体的聚焦程度,将高范围增强滤波器、中范围增强滤波器、低范围增强滤波器或负增强滤波器(即平滑滤波器)确定为由滤波器增强的区域的类型。
(11)根据(1)至(10)中任一项所述的图像处理方法,包括
获取医学领域中与活体、生物体、材料或病理学相关的被摄体的医学图像。
(12)根据(1)至(11)中任一项所述的图像处理方法,包括获取由作为成像设备的显微镜捕获的医学图像。
(13)一种图像处理设备,包括:
获取单元,所述获取单元被配置为获取由显微镜捕获的被摄体的医学图像;以及
确定单元,所述确定单元被配置为根据被摄体的聚焦程度来确定要应用于医学图像的滤波器的强度,所述滤波器改善医学图像的图像质量。
(14)一种图像处理系统,包括:
成像设备,所述成像设备被配置为对被摄体成像;以及
图像处理设备,所述图像处理设备被配置为包括用于处理与要由成像设备成像的目标相对应的医学图像的软件,
其中,所述软件根据被摄体的聚焦程度来确定要应用于由成像设备捕获的医学图像的滤波器的强度。
附图标记列表
1 图像处理系统
10 成像设备
100 图像处理设备
110 通信单元
120 存储单元
121 医学图像存储单元
122 增强滤波器存储单元
130 控制单元
131 获取单元
132 计算单元
133 确定单元
N 网络
Claims (14)
1.一种图像处理方法,包括:由计算机
获取由成像设备捕获的医学图像;以及
根据所述医学图像的聚焦程度确定要应用于所述医学图像的滤波器的强度。
2.根据权利要求1所述的图像处理方法,包括:
根据由指示所述医学图像的模糊程度的模糊函数计算的聚焦程度,确定要应用于所述医学图像的滤波器的强度。
3.根据权利要求2所述的图像处理方法,包括:
根据由所述模糊函数根据被摄体在垂直于所述医学图像的方向上的厚度方向计算的聚焦程度,确定要应用于所述医学图像的滤波器的强度。
4.根据权利要求3所述的图像处理方法,包括:
根据通过对所述医学图像中的预定像素与外围像素之间的差求和而获得的特征量计算的聚焦程度,确定要应用于所述医学图像的滤波器的强度。
5.根据权利要求4所述的图像处理方法,包括:
基于根据所述模糊函数估计的估计聚焦程度和通过对所述医学图像中的预定像素与外围像素之间的差求和而获得的特征量计算的聚焦程度,确定要应用于所述医学图像的滤波器的强度。
6.根据权利要求5所述的图像处理方法,包括:
根据所述被摄体的聚焦程度,确定合成率,所述合成率指示通过将多个不同的滤波器应用于所述医学图像而生成的合成图像的合成比率。
7.根据权利要求6所述的图像处理方法,包括:
确定通过将具有不同强度的多个滤波器应用于预定医学图像而生成的所述合成图像的合成率。
8.根据权利要求6所述的图像处理方法,包括:
根据所述被摄体对具有不同Z值的多个所述医学图像的聚焦程度,确定通过应用滤波器而生成的所述合成图像的合成率。
9.根据权利要求6所述的图像处理方法,包括:
根据所述被摄体的聚焦程度,选择性地确定要应用于用于生成所述合成图像的医学图像的多个滤波器。
10.根据权利要求9所述的图像处理方法,包括:
根据所述被摄体的聚焦程度,将高范围增强滤波器、中范围增强滤波器、低范围增强滤波器或负增强滤波器确定为由滤波器增强的区域的类型,负增强滤波器即平滑滤波器。
11.根据权利要求10所述的图像处理方法,包括:
获取医学领域中与活体、生物体、材料或病理学相关的被摄体的医学图像。
12.根据权利要求11所述的图像处理方法,包括:
获取由作为所述成像设备的显微镜捕获的医学图像。
13.一种图像处理设备,包括:
获取单元,被配置为获取由显微镜捕获的被摄体的医学图像;以及
确定单元,被配置为根据所述被摄体的聚焦程度来确定要应用于所述医学图像的滤波器的强度,所述滤波器改善所述医学图像的图像质量。
14.一种图像处理系统,包括:
成像设备,被配置为对被摄体进行成像;以及
图像处理设备,被配置为包括用于处理与要由所述成像设备成像的目标相对应的医学图像的软件,
其中,所述软件根据所述被摄体的聚焦程度来确定要应用于由所述成像设备捕获的医学图像的滤波器的强度。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019211390 | 2019-11-22 | ||
JP2019-211390 | 2019-11-22 | ||
PCT/JP2020/037068 WO2021100328A1 (ja) | 2019-11-22 | 2020-09-30 | 画像処理方法、画像処理装置および画像処理システム |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114730070A true CN114730070A (zh) | 2022-07-08 |
Family
ID=75980520
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202080079048.XA Pending CN114730070A (zh) | 2019-11-22 | 2020-09-30 | 图像处理方法、图像处理装置和图像处理系统 |
Country Status (5)
Country | Link |
---|---|
US (1) | US20220392031A1 (zh) |
EP (1) | EP4047408A4 (zh) |
JP (1) | JPWO2021100328A1 (zh) |
CN (1) | CN114730070A (zh) |
WO (1) | WO2021100328A1 (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP7280107B2 (ja) * | 2019-05-10 | 2023-05-23 | 株式会社エビデント | 画像処理方法、プログラム、画像処理装置、画像処理システム、及び、顕微鏡システム |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2656778A1 (en) * | 2010-12-24 | 2013-10-30 | Olympus Corporation | Endoscope device and program |
US20140022347A1 (en) * | 2012-07-19 | 2014-01-23 | Sony Corporation | Method and apparatus for simulating depth of field (dof) in microscopy |
US20140085507A1 (en) * | 2012-09-21 | 2014-03-27 | Bruce Harold Pillman | Controlling the sharpness of a digital image |
US20140293117A1 (en) * | 2011-12-02 | 2014-10-02 | Canon Kabushiki Kaisha | Image processing apparatus and image processing method |
CN104662463A (zh) * | 2012-09-28 | 2015-05-27 | 佳能株式会社 | 图像处理装置、成像系统和图像处理系统 |
US20150286044A1 (en) * | 2012-11-30 | 2015-10-08 | L&T Technology Services Limited | Method and system for extended depth of field calculation for microscopic images |
CN105026977A (zh) * | 2013-03-13 | 2015-11-04 | 索尼公司 | 信息处理装置、信息处理方法和信息处理程序 |
US20170142358A1 (en) * | 2015-11-12 | 2017-05-18 | Canon Kabushiki Kaisha | Image processing apparatus and image processing method |
WO2017144482A1 (en) * | 2016-02-22 | 2017-08-31 | Koninklijke Philips N.V. | System for generating a synthetic 2d image with an enhanced depth of field of a biological sample |
JP2017158764A (ja) * | 2016-03-09 | 2017-09-14 | ソニー株式会社 | 画像処理装置、画像処理方法、及び記録媒体 |
JP2019149719A (ja) * | 2018-02-27 | 2019-09-05 | キヤノン株式会社 | 画像処理装置及び方法、及び撮像装置 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5994712A (ja) * | 1982-11-22 | 1984-05-31 | Olympus Optical Co Ltd | 合焦検出方法 |
JPH10243289A (ja) * | 1997-02-27 | 1998-09-11 | Olympus Optical Co Ltd | 画像加算装置 |
JP2011090221A (ja) * | 2009-10-23 | 2011-05-06 | Sony Corp | 顕微鏡、合焦位置検出方法及び合焦位置検出プログラム |
JP2011197283A (ja) | 2010-03-18 | 2011-10-06 | Sony Corp | 合焦装置、合焦方法、合焦プログラム及び顕微鏡 |
JP7238381B2 (ja) * | 2017-12-21 | 2023-03-14 | 株式会社ニコン | 画像処理装置、画像処理プログラム、画像処理方法および顕微鏡 |
-
2020
- 2020-09-30 JP JP2021558198A patent/JPWO2021100328A1/ja active Pending
- 2020-09-30 WO PCT/JP2020/037068 patent/WO2021100328A1/ja unknown
- 2020-09-30 CN CN202080079048.XA patent/CN114730070A/zh active Pending
- 2020-09-30 US US17/776,985 patent/US20220392031A1/en active Pending
- 2020-09-30 EP EP20890703.0A patent/EP4047408A4/en active Pending
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2656778A1 (en) * | 2010-12-24 | 2013-10-30 | Olympus Corporation | Endoscope device and program |
US20140293117A1 (en) * | 2011-12-02 | 2014-10-02 | Canon Kabushiki Kaisha | Image processing apparatus and image processing method |
US20140022347A1 (en) * | 2012-07-19 | 2014-01-23 | Sony Corporation | Method and apparatus for simulating depth of field (dof) in microscopy |
US20140085507A1 (en) * | 2012-09-21 | 2014-03-27 | Bruce Harold Pillman | Controlling the sharpness of a digital image |
CN104662463A (zh) * | 2012-09-28 | 2015-05-27 | 佳能株式会社 | 图像处理装置、成像系统和图像处理系统 |
US20150286044A1 (en) * | 2012-11-30 | 2015-10-08 | L&T Technology Services Limited | Method and system for extended depth of field calculation for microscopic images |
CN105026977A (zh) * | 2013-03-13 | 2015-11-04 | 索尼公司 | 信息处理装置、信息处理方法和信息处理程序 |
US20170142358A1 (en) * | 2015-11-12 | 2017-05-18 | Canon Kabushiki Kaisha | Image processing apparatus and image processing method |
WO2017144482A1 (en) * | 2016-02-22 | 2017-08-31 | Koninklijke Philips N.V. | System for generating a synthetic 2d image with an enhanced depth of field of a biological sample |
JP2017158764A (ja) * | 2016-03-09 | 2017-09-14 | ソニー株式会社 | 画像処理装置、画像処理方法、及び記録媒体 |
JP2019149719A (ja) * | 2018-02-27 | 2019-09-05 | キヤノン株式会社 | 画像処理装置及び方法、及び撮像装置 |
Also Published As
Publication number | Publication date |
---|---|
US20220392031A1 (en) | 2022-12-08 |
WO2021100328A1 (ja) | 2021-05-27 |
EP4047408A1 (en) | 2022-08-24 |
JPWO2021100328A1 (zh) | 2021-05-27 |
EP4047408A4 (en) | 2022-12-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9332190B2 (en) | Image processing apparatus and image processing method | |
US9881373B2 (en) | Image generating apparatus and image generating method | |
Russ | Computer-assisted microscopy: the measurement and analysis of images | |
US9224193B2 (en) | Focus stacking image processing apparatus, imaging system, and image processing system | |
JP5254441B2 (ja) | フロー式粒子画像解析方法及び装置 | |
JP7252190B2 (ja) | 生物試料の向上された被写界深度の合成2d画像を生成するシステム | |
JPH10509817A (ja) | 信号復元方法および装置 | |
JP2015192238A (ja) | 画像データ生成装置および画像データ生成方法 | |
WO2015156378A1 (en) | Image processing apparatus, image processing method, and image processing system | |
JP2015108837A (ja) | 画像処理装置及び画像処理方法 | |
CN109001902A (zh) | 基于图像融合的显微镜聚焦方法 | |
JP6362062B2 (ja) | 画像生成装置および画像生成方法 | |
CN114730070A (zh) | 图像处理方法、图像处理装置和图像处理系统 | |
JP2015057682A (ja) | 画像生成装置および画像生成方法 | |
JP7030986B2 (ja) | 画像生成装置、画像生成方法および画像生成プログラム | |
US20200074628A1 (en) | Image processing apparatus, imaging system, image processing method and computer readable recoding medium | |
JP2015191362A (ja) | 画像データ生成装置および画像データ生成方法 | |
Kumar et al. | Performance evaluation of joint filtering and histogram equalization techniques for retinal fundus image enhancement | |
JP2023542619A (ja) | 試料のデジタル画像の品質管理のためのコンピュータ実装方法 | |
JP2023033982A (ja) | 画像処理装置、画像処理システム、画像の鮮鋭化方法、及び、プログラム | |
Tadrous | A method of PSF generation for 3D brightfield deconvolution | |
Willis et al. | Developments in three‐dimensional stereo brightfield microscopy | |
Pronichev et al. | Cytology 3D structure formation based on optical microscopy images | |
WO2022113365A1 (ja) | 合焦方法、観察装置、およびプログラム | |
Keating et al. | Improved spatial resolution in ratio images using computational confocal techniques |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |