CN107507173A - 一种全切片图像的无参考清晰度评估方法及系统 - Google Patents

一种全切片图像的无参考清晰度评估方法及系统 Download PDF

Info

Publication number
CN107507173A
CN107507173A CN201710697583.1A CN201710697583A CN107507173A CN 107507173 A CN107507173 A CN 107507173A CN 201710697583 A CN201710697583 A CN 201710697583A CN 107507173 A CN107507173 A CN 107507173A
Authority
CN
China
Prior art keywords
image
edge
image block
picture
definition
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
Application number
CN201710697583.1A
Other languages
English (en)
Other versions
CN107507173B (zh
Inventor
吴开杰
谷朝臣
关新平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Pu Huasen Biotechnology Co ltd
Original Assignee
Shanghai Jiaotong University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Jiaotong University filed Critical Shanghai Jiaotong University
Priority to CN201710697583.1A priority Critical patent/CN107507173B/zh
Publication of CN107507173A publication Critical patent/CN107507173A/zh
Application granted granted Critical
Publication of CN107507173B publication Critical patent/CN107507173B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/12Edge-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/136Segmentation; Edge detection involving thresholding
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/194Segmentation; Edge detection involving foreground-background segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20016Hierarchical, coarse-to-fine, multiscale or multiresolution image processing; Pyramid transform
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20021Dividing image into blocks, subimages or windows
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20024Filtering details
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30168Image quality inspection

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Quality & Reliability (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

本发明提供了一种全切片图像的无参考清晰度评估方法及系统,采用金字塔多层数据结构存储数字图像,将金字塔多层数据结构的各层图像平均分割为若干图像块;通过大津阈值分割法在图像块的梯度幅度直方图中划分出强边缘,由强边缘的梯度最大值和宽度计算得到强边缘强度,并利用背景复杂度对强边缘强度进行修正,确定单个图像块的清晰度;设定清晰度的标准值,判定当前图像块是清晰图像还是模糊图像;根据各层清晰图像块所占比例,最终确定各层图像是否为清晰图像。本发明避免了图像内容差异性导致的图像清晰度得分差异性,修正了图像背景复杂度对图像清晰度的影响,实现了对部分清晰部分模糊图像的有效鉴别。

Description

一种全切片图像的无参考清晰度评估方法及系统
技术领域
本发明涉及数字图像处理技术领域,具体地,涉及一种全切片图像的无参考清晰度评估方法及系统。
背景技术
随着信息化技术的发展,数字化全切片病理图像已在临床诊断、病理学研究领域广泛应用。由于切片制作、显微扫描过程存在各种不确定因素,容易造成病理切片的成像质量下降。图像清晰度下降不但直接影响病理切片的诊断质量,而且无法作为重要样本用于数据挖掘等智能诊断研究。因此,快速有效的数字病理图像清晰度评估方法对于数字病理的发展至关重要。
根据对参考图像的依赖程度,图像清晰度的客观评估方法可以分为全参考方法和无参考方法。全参考图像清晰度评估方法考察待测图像与参考图像的清晰度差异,其评价准确率往往高于其他客观图像清晰度评估方法,算法实现难度也相对较低,但是全参考方法的数据量和计算量也远高于其他方法。无参考图像清晰度评估方法不需要原始图像作为参考,虽然准确性较全参考方法可能有所下降,但是它的计算效率更高,应用也更加广泛。目前,数字病理图像清晰度评估方法在实际应用中普遍存在以下不足:1)、清晰度严重受图像内容影响,同等清晰度下,内容复杂的图像比内容简单的图像得分明显偏高;2)、部分清晰部分模糊的图像得分与全清晰图像非常接近。这些问题使现有的清晰度评估方法不适用普遍的清晰度评估。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种全切片图像的无参考清晰度评估方法及系统,实现高效率病理图像质量评估,以满足临床诊断、病理学研究对可靠病理数据源的迫切需求。
根据本发明提供的全切片图像的无参考清晰度评估方法,包括:采用金字塔多层数据结构存储数字图像,将金字塔多层数据结构的各层图像平均分割为若干图像块;通过大津阈值分割法在图像块的梯度幅度直方图中划分出强边缘,由强边缘的梯度最大值和宽度计算得到强边缘强度,并利用背景复杂度对强边缘强度进行修正,确定单个图像块的清晰度;设定清晰度的标准值,判定当前图像块是清晰图像还是模糊图像;根据各层清晰图像块所占比例,最终确定各层图像是否为清晰图像。
优选的,全切片图像的无参考清晰度评估方法具体包括步骤:
金字塔多层数据结构构建步骤:由高倍物镜扫描图像,将图像按金字塔多层数据结构存储,各层对应不同图像放大倍数,各层平均分割为若干长宽一致的图像块,针对某层图像,将当前图像块转换为灰度图像,并用高斯滤波器进行滤波;
清晰度初始值计算步骤:使用水平和竖直方向的索贝尔滤波器分别对当前灰度图像进行滤波,得到水平和竖直方向的梯度幅度,在此基础上计算总的梯度幅度和梯度方向;根据大津阈值分割法在梯度幅度直方图中确定边缘梯度幅度阈值,边缘梯度幅度大于阈值为强边缘;计算获得各段强边缘强度,对多段强边缘强度进行加强平均,获得当前图像块的强边缘强度,并作为当前图像块的清晰度初始值;
清晰度计算步骤:将图像块的梯度幅度中低于边缘梯度幅度阈值的区域看作背景图像,对其对应的梯度幅度做加权平均得到背景复杂度;计算获得当前图像块的所述背景复杂度修正值,乘以所述清晰度初始值,获得当前图像块的清晰度;
判断步骤:重复清晰度初始值计算步骤以及清晰度计算步骤,计算所有图像块的清晰度,设定清晰度标准值,判断当前图像块是清晰图像块还是模糊图像块,根据清晰图像块的占比确定各层图像是否为清晰图像。
优选的,所述计算获得各段强边缘强度包括:
在强边缘的梯度正方向和负方向上寻找的梯度幅度的局部最大值,作为边缘梯度最大值;以所述边缘梯度最大值所对应的像素点作为当前边缘的中心点,分别沿着梯度方向的正方向与负方向寻找两个梯度幅度的局部最小值,两个局部最小值所在像素点之间的欧氏距离为强边缘宽度;所述强边缘强度得分与所述边缘梯度最大值正相关,与所述强边缘宽度负相关,设置指数提高所述强边缘强度对所述强边缘宽度的敏感度,计算方法如下:
优选的,所述计算获得当前图像块的所述背景复杂度修正值包括:
设定所述边缘梯度幅度阈值的一半为正负影响阈值,如果所述背景复杂度高于所述正负影响阈值,则输入图像的背景负向影响当前图像的清晰度得分,如果所述背景复杂度低于所述正负影响阈值,则输入图像的背景正向影响当前图像的清晰度得分;最后,采用背景像素数量占图像块像素数量的比例进一步修正背景影响,获得背景影响修正系数:
优选的,所述当前图像块的清晰度的计算方法如下:
根据本发明提供的全切片图像的无参考清晰度评估系统,包括:采用金字塔多层数据结构存储数字图像,将金字塔多层数据结构的各层图像平均分割为若干图像块;通过大津阈值分割法在图像块的梯度幅度直方图中划分出强边缘,由强边缘的梯度最大值和宽度计算得到强边缘强度,并利用背景复杂度对强边缘强度进行修正,确定单个图像块的清晰度;设定清晰度的标准值,判定当前图像块是清晰图像还是模糊图像;根据各层清晰图像块所占比例,最终确定各层图像是否为清晰图像。
优选的,全切片图像的无参考清晰度评估系统具体包括:
金字塔多层数据结构构建模块:由高倍物镜扫描图像,将图像按金字塔多层数据结构存储,各层对应不同图像放大倍数,各层平均分割为若干长宽一致的图像块,针对某层图像,将当前图像块转换为灰度图像,并用高斯滤波器进行滤波;
清晰度初始值计算模块:使用水平和竖直方向的索贝尔滤波器分别对当前灰度图像进行滤波,得到水平和竖直方向的梯度幅度,在此基础上计算总的梯度幅度和梯度方向;根据大津阈值分割法在梯度幅度直方图中确定边缘梯度幅度阈值,边缘梯度幅度大于阈值为强边缘;计算获得各段强边缘强度,对多段强边缘强度进行加强平均,获得当前图像块的强边缘强度,并作为当前图像块的清晰度初始值;
清晰度计算模块:将图像块的梯度幅度中低于边缘梯度幅度阈值的区域看作背景图像,对其对应的梯度幅度做加权平均得到背景复杂度;计算获得当前图像块的所述背景复杂度修正值,乘以所述清晰度初始值,获得当前图像块的清晰度;
判断步骤:重复清晰度初始值计算模块以及清晰度计算模块,计算所有图像块的清晰度,设定清晰度标准值,判断当前图像块是清晰图像块还是模糊图像块,根据清晰图像块的占比确定各层图像是否为清晰图像。
优选的,所述计算获得各段强边缘强度包括:
在强边缘的梯度正方向和负方向上寻找的梯度幅度的局部最大值,作为边缘梯度最大值;以所述边缘梯度最大值所对应的像素点作为当前边缘的中心点,分别沿着梯度方向的正方向与负方向寻找两个梯度幅度的局部最小值,两个局部最小值所在像素点之间的欧氏距离为强边缘宽度;所述强边缘强度得分与所述边缘梯度最大值正相关,与所述强边缘宽度负相关,设置指数提高所述强边缘强度对所述强边缘宽度的敏感度,计算方法如下:
优选的,所述计算获得当前图像块的所述背景复杂度修正值包括:
设定所述边缘梯度幅度阈值的一半为正负影响阈值,如果所述背景复杂度高于所述正负影响阈值,则输入图像的背景负向影响当前图像的清晰度得分,如果所述背景复杂度低于所述正负影响阈值,则输入图像的背景正向影响当前图像的清晰度得分;最后,采用背景像素数量占图像块像素数量的比例进一步修正背景影响,获得背景影响修正系数:
优选的,所述当前图像块的清晰度的计算方法如下:
与现有技术相比,本发明具有如下的有益效果:
本发明避免了图像内容差异性导致的图像清晰度得分差异性,修正了图像背景复杂度对图像清晰度的影响,实现了对部分清晰部分模糊图像的有效鉴别,供用户直接评估扫描得到的图像质量,与人为质量评估相比可以提高评估效率,促进形成统一的图像质量评估机制,为科学研究提供专业且合格的图像数据库,促进数据挖掘、深度学习等研究手段在例如智能病理诊断等领域的进一步应用。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为本发明实施例所述的病理图像清晰度评估流程;
图2为本发明实施例所述的强边缘宽度计算示例。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变化和改进。这些都属于本发明的保护范围。
如图1所示,本发明提供的全切片图像的无参考清晰度评估方法包括:采用金字塔多层数据结构存储数字图像,将金字塔多层数据结构的各层图像平均分割为若干图像块;通过大津阈值分割法在图像块的梯度幅度直方图中划分出强边缘,由强边缘的梯度最大值和宽度计算得到强边缘强度,并利用背景复杂度对强边缘强度进行修正,确定单个图像块的清晰度;设定清晰度的标准值,判定当前图像块是清晰图像还是模糊图像;根据各层清晰图像块所占比例,最终确定各层图像是否为清晰图像。
具体可分为四个步骤:金字塔多层数据结构构建步骤、清晰度初始值计算步骤、清晰度计算步骤以及判断步骤。
金字塔多层数据结构构建步骤:
通常采用高倍物镜进行病理切片数字化扫描,并将数字病理图像按金字塔多层数据结构存储,各层图像按照低分辨率到高分辨率对应2×、4×、10×、20×、40×等不同图像放大倍数,各层图像平均分割为若干256×256像素的长宽一致的图像块;以20×图层病理图像为例,针对某层图像,将当前RGB图像块转换为灰度图像,并用高斯滤波器进行滤波。
清晰度初始值计算步骤:
使用水平和竖直方向的索贝尔滤波器sobel_x,sobel_y分别对当前灰度图像进行滤波,得到水平和竖直方向的梯度幅度Gx和Gy
在此基础上计算总的梯度幅度G和梯度方向θ:
如图2所示,根据大津阈值分割法在梯度幅度直方图中确定边缘梯度幅度阈值T,边缘梯度幅度大于阈值T即为强边缘,根据θ对这些边缘对应的梯度方向以做采样。在相应的梯度正方向和负方向上寻找梯度幅度的局部最大值max,该最大值作为当前边缘的边缘梯度最大值,该最大值对应的像素点作为当前边缘的中心点(A1)。然后从中心点开始,分别沿着梯度方向的正方向与负方向寻找两个梯度幅度的局部最小值,寻找局部最小值时因消除抖动的影响,两个局部最小值所在像素点(A2,A3)之间的欧氏距离即为当前边缘的宽度w。
一般来说,边缘梯度幅度的局部最大值越大,宽度值越小,说明此边缘越清晰,但是如果当前边缘的宽度等于零,说明此边缘非常模糊。为了提高强边缘强度对边缘宽度的敏感度,在此将边缘强度最大值的二次方或更高次方除以边缘宽度,计算公式如下:
清晰度初始值计算步骤:
计算图像块清晰度初始值,将上式得到的各个强边缘的强度s加权平均,得到当前图像块的边缘强度均值,作为图像块清晰度的初始值S0
清晰度计算步骤:
通过背景复杂度修正图像块清晰度,主要流程如下:将梯度幅度矩阵中低于阈值T的区域均为背景图像,对其对应的梯度幅度做加权平均得到背景复杂度c,并统计背景所占像素数量Mbg,设设定所述边缘梯度幅度阈值的一半为正负影响阈值T/2,如果背景复杂度c高于正负影响阈值T/2,则认为输入图像的背景负向影响当前图像的清晰度得分,如果背景复杂度c低于正负影响阈值T/2,则认为输入图像的背景正向影响当前图像的清晰度得分,结合背景像素数量Mbg占图像块全部像素数量P的比例,获得背景影响修正系数,进一步修正图像清晰度初始值S1。背景影响修正系数k在-1到1之间取值:
将当前图像块清晰度初始值乘以其背景复杂度修正值,获得当前图像块的清晰度(强边缘强度修正值):
S=S0×k
判断步骤:
重复以上清晰度初始值计算步骤、清晰度计算步骤,遍历计算所有图像块的清晰度;设定清晰度标准值Sstd(B11),用于判定当前图像块是清晰图像还是模糊图像:
根据各层清晰图像块所占比例,来确定各层病理图像是否为清晰图像:
其中Hi为单个图像块的清晰程度,图像块总数为Q。针对不同放大倍数的图层,Sharp值不同,代表不同清晰度程度,例如对于2×图层,Sharp>90%该图层为清晰图像;对于10×图层,Sharp>85%该图层为清晰图像;对于20×图层,Sharp>80%该图层为清晰图像。
根据上述方法,本发明另提供一种全切片图像的无参考清晰度评估系统,包括:采用金字塔多层数据结构存储数字图像,将金字塔多层数据结构的各层图像平均分割为若干图像块;通过大津阈值分割法在图像块的梯度幅度直方图中划分出强边缘,由强边缘的梯度最大值和宽度计算得到强边缘强度,并利用背景复杂度对强边缘强度进行修正,确定单个图像块的清晰度;设定清晰度的标准值,判定当前图像块是清晰图像还是模糊图像;根据各层清晰图像块所占比例,最终确定各层图像是否为清晰图像。
具体可分为四个模块:金字塔多层数据结构构建模块、清晰度初始值计算模块、清晰度计算模块以及判断模块。
金字塔多层数据结构构建模块:
通常采用高倍物镜进行病理切片数字化扫描,并将数字病理图像按金字塔多层数据结构存储,各层图像按照低分辨率到高分辨率对应2×、4×、10×、20×、40×等不同图像放大倍数,各层图像平均分割为若干256×256像素的长宽一致的图像块;以20×图层病理图像为例,针对某层图像,将当前RGB图像块转换为灰度图像,并用高斯滤波器进行滤波。
清晰度初始值计算模块:
使用水平和竖直方向的索贝尔滤波器sobel_x,sobel_y分别对当前灰度图像进行滤波,得到水平和竖直方向的梯度幅度Gx和Gy
在此基础上计算总的梯度幅度G和梯度方向θ:
如图2所示,根据大津阈值分割法在梯度幅度直方图中确定边缘梯度幅度阈值T,边缘梯度幅度大于阈值T即为强边缘,根据θ对这些边缘对应的梯度方向以做采样。在相应的梯度正方向和负方向上寻找梯度幅度的局部最大值max,该最大值作为当前边缘的边缘梯度最大值,该最大值对应的像素点作为当前边缘的中心点(A1)。然后从中心点开始,分别沿着梯度方向的正方向与负方向寻找两个梯度幅度的局部最小值,寻找局部最小值时因消除抖动的影响,两个局部最小值所在像素点(A2,A3)之间的欧氏距离即为当前边缘的宽度w。
一般来说,边缘梯度幅度的局部最大值越大,宽度值越小,说明此边缘越清晰,但是如果当前边缘的宽度等于零,说明此边缘非常模糊。为了提高强边缘强度对边缘宽度的敏感度,在此将边缘强度最大值的二次方或更高次方除以边缘宽度,计算公式如下:
清晰度初始值计算模块:
计算图像块清晰度初始值,将上式得到的各个强边缘的强度s加权平均,得到当前图像块的边缘强度均值,作为图像块清晰度的初始值S0
清晰度计算模块:
通过背景复杂度修正图像块清晰度,主要流程如下:将梯度幅度矩阵中低于阈值T的区域均为背景图像,对其对应的梯度幅度做加权平均得到背景复杂度c,并统计背景所占像素数量Mbg,设设定所述边缘梯度幅度阈值的一半为正负影响阈值T/2,如果背景复杂度c高于正负影响阈值T/2,则认为输入图像的背景负向影响当前图像的清晰度得分,如果背景复杂度c低于正负影响阈值T/2,则认为输入图像的背景正向影响当前图像的清晰度得分,结合背景像素数量Mbg占图像块全部像素数量P的比例,获得背景影响修正系数,进一步修正图像清晰度初始值S1。背景影响修正系数k在-1到1之间取值:
将当前图像块清晰度初始值乘以其背景复杂度修正值,获得当前图像块的清晰度(强边缘强度修正值):
S=S0×k
判断模块:
重复以上清晰度初始值计算、清晰度计算,遍历计算所有图像块的清晰度;设定清晰度标准值Sstd(B11),用于判定当前图像块是清晰图像还是模糊图像:
根据各层清晰图像块所占比例,来确定各层病理图像是否为清晰图像:
其中Hi为单个图像块的清晰程度,图像块总数为Q。针对不同放大倍数的图层,Sharp值不同,代表不同清晰度程度,例如对于2×图层,Sharp>90%该图层为清晰图像;对于10×图层,Sharp>85%该图层为清晰图像;对于20×图层,Sharp>80%该图层为清晰图像。
本领域技术人员知道,除了以纯计算机可读程序代码方式实现本发明提供的系统及其各个装置、模块、单元以外,完全可以通过将方法步骤进行逻辑编程来使得本发明提供的系统及其各个装置、模块、单元以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同功能。所以,本发明提供的系统及其各项装置、模块、单元可以被认为是一种硬件部件,而对其内包括的用于实现各种功能的装置、模块、单元也可以视为硬件部件内的结构;也可以将用于实现各种功能的装置、模块、单元视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变化或修改,这并不影响本发明的实质内容。在不冲突的情况下,本申请的实施例和实施例中的特征可以任意相互组合。

Claims (10)

1.一种全切片图像的无参考清晰度评估方法,其特征在于,包括:采用金字塔多层数据结构存储数字图像,将金字塔多层数据结构的各层图像平均分割为若干图像块;通过大津阈值分割法在图像块的梯度幅度直方图中划分出强边缘,由强边缘的梯度最大值和宽度计算得到强边缘强度,并利用背景复杂度对强边缘强度进行修正,确定单个图像块的清晰度;设定清晰度的标准值,判定当前图像块是清晰图像还是模糊图像;根据各层清晰图像块所占比例,最终确定各层图像是否为清晰图像。
2.根据权利要求1所述的全切片图像的无参考清晰度评估方法,其特征在于,具体包括步骤:
金字塔多层数据结构构建步骤:由高倍物镜扫描图像,将图像按金字塔多层数据结构存储,各层对应不同图像放大倍数,各层平均分割为若干长宽一致的图像块,针对某层图像,将当前图像块转换为灰度图像,并用高斯滤波器进行滤波;
清晰度初始值计算步骤:使用水平和竖直方向的索贝尔滤波器分别对当前灰度图像进行滤波,得到水平和竖直方向的梯度幅度,在此基础上计算总的梯度幅度和梯度方向;根据大津阈值分割法在梯度幅度直方图中确定边缘梯度幅度阈值,边缘梯度幅度大于阈值为强边缘;计算获得各段强边缘强度,对多段强边缘强度进行加强平均,获得当前图像块的强边缘强度,并作为当前图像块的清晰度初始值;
清晰度计算步骤:将图像块的梯度幅度中低于边缘梯度幅度阈值的区域看作背景图像,对其对应的梯度幅度做加权平均得到背景复杂度;计算获得当前图像块的所述背景复杂度修正值,乘以所述清晰度初始值,获得当前图像块的清晰度;
判断步骤:重复清晰度初始值计算步骤以及清晰度计算步骤,计算所有图像块的清晰度,设定清晰度标准值,判断当前图像块是清晰图像块还是模糊图像块,根据清晰图像块的占比确定各层图像是否为清晰图像。
3.根据权利要求2所述的全切片图像的无参考清晰度评估方法,其特征在于,所述计算获得各段强边缘强度包括:
在强边缘的梯度正方向和负方向上寻找的梯度幅度的局部最大值,作为边缘梯度最大值;以所述边缘梯度最大值所对应的像素点作为当前边缘的中心点,分别沿着梯度方向的正方向与负方向寻找两个梯度幅度的局部最小值,两个局部最小值所在像素点之间的欧氏距离为强边缘宽度;所述强边缘强度得分与所述边缘梯度最大值正相关,与所述强边缘宽度负相关,设置指数提高所述强边缘强度对所述强边缘宽度的敏感度,计算方法如下:
4.根据权利要求2所述的全切片图像的无参考清晰度评估方法,其特征在于,所述计算获得当前图像块的所述背景复杂度修正值包括:
设定所述边缘梯度幅度阈值的一半为正负影响阈值,如果所述背景复杂度高于所述正负影响阈值,则输入图像的背景负向影响当前图像的清晰度得分,如果所述背景复杂度低于所述正负影响阈值,则输入图像的背景正向影响当前图像的清晰度得分;最后,采用背景像素数量占图像块像素数量的比例进一步修正背景影响,获得背景影响修正系数:
5.根据权利要求2所述的全切片图像的无参考清晰度评估方法,其特征在于,所述当前图像块的清晰度的计算方法如下:
6.一种全切片图像的无参考清晰度评估系统,其特征在于,包括:采用金字塔多层数据结构存储数字图像,将金字塔多层数据结构的各层图像平均分割为若干图像块;通过大津阈值分割法在图像块的梯度幅度直方图中划分出强边缘,由强边缘的梯度最大值和宽度计算得到强边缘强度,并利用背景复杂度对强边缘强度进行修正,确定单个图像块的清晰度;设定清晰度的标准值,判定当前图像块是清晰图像还是模糊图像;根据各层清晰图像块所占比例,最终确定各层图像是否为清晰图像。
7.根据权利要求6所述的全切片图像的无参考清晰度评估系统,其特征在于,具体包括:
金字塔多层数据结构构建模块:由高倍物镜扫描图像,将图像按金字塔多层数据结构存储,各层对应不同图像放大倍数,各层平均分割为若干长宽一致的图像块,针对某层图像,将当前图像块转换为灰度图像,并用高斯滤波器进行滤波;
清晰度初始值计算模块:使用水平和竖直方向的索贝尔滤波器分别对当前灰度图像进行滤波,得到水平和竖直方向的梯度幅度,在此基础上计算总的梯度幅度和梯度方向;根据大津阈值分割法在梯度幅度直方图中确定边缘梯度幅度阈值,边缘梯度幅度大于阈值为强边缘;计算获得各段强边缘强度,对多段强边缘强度进行加强平均,获得当前图像块的强边缘强度,并作为当前图像块的清晰度初始值;
清晰度计算模块:将图像块的梯度幅度中低于边缘梯度幅度阈值的区域看作背景图像,对其对应的梯度幅度做加权平均得到背景复杂度;计算获得当前图像块的所述背景复杂度修正值,乘以所述清晰度初始值,获得当前图像块的清晰度;
判断步骤:重复清晰度初始值计算模块以及清晰度计算模块,计算所有图像块的清晰度,设定清晰度标准值,判断当前图像块是清晰图像块还是模糊图像块,根据清晰图像块的占比确定各层图像是否为清晰图像。
8.根据权利要求7所述的全切片图像的无参考清晰度评估系统,其特征在于,所述计算获得各段强边缘强度包括:
在强边缘的梯度正方向和负方向上寻找的梯度幅度的局部最大值,作为边缘梯度最大值;以所述边缘梯度最大值所对应的像素点作为当前边缘的中心点,分别沿着梯度方向的正方向与负方向寻找两个梯度幅度的局部最小值,两个局部最小值所在像素点之间的欧氏距离为强边缘宽度;所述强边缘强度得分与所述边缘梯度最大值正相关,与所述强边缘宽度负相关,设置指数提高所述强边缘强度对所述强边缘宽度的敏感度,计算方法如下:
9.根据权利要求7所述的全切片图像的无参考清晰度评估系统,其特征在于,所述计算获得当前图像块的所述背景复杂度修正值包括:
设定所述边缘梯度幅度阈值的一半为正负影响阈值,如果所述背景复杂度高于所述正负影响阈值,则输入图像的背景负向影响当前图像的清晰度得分,如果所述背景复杂度低于所述正负影响阈值,则输入图像的背景正向影响当前图像的清晰度得分;最后,采用背景像素数量占图像块像素数量的比例进一步修正背景影响,获得背景影响修正系数:
10.根据权利要求7所述的全切片图像的无参考清晰度评估系统,其特征在于,所述当前图像块的清晰度的计算方法如下:
CN201710697583.1A 2017-08-15 2017-08-15 一种全切片图像的无参考清晰度评估方法及系统 Active CN107507173B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710697583.1A CN107507173B (zh) 2017-08-15 2017-08-15 一种全切片图像的无参考清晰度评估方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710697583.1A CN107507173B (zh) 2017-08-15 2017-08-15 一种全切片图像的无参考清晰度评估方法及系统

Publications (2)

Publication Number Publication Date
CN107507173A true CN107507173A (zh) 2017-12-22
CN107507173B CN107507173B (zh) 2021-07-27

Family

ID=60691868

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710697583.1A Active CN107507173B (zh) 2017-08-15 2017-08-15 一种全切片图像的无参考清晰度评估方法及系统

Country Status (1)

Country Link
CN (1) CN107507173B (zh)

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108460766A (zh) * 2018-04-12 2018-08-28 四川和生视界医药技术开发有限公司 一种视网膜图像清晰度评估方法以及评估装置
CN109191457A (zh) * 2018-09-21 2019-01-11 中国人民解放军总医院 一种病理图像质量有效性识别方法
CN110148147A (zh) * 2018-11-07 2019-08-20 腾讯大地通途(北京)科技有限公司 图像检测方法、装置、存储介质和电子装置
CN110458789A (zh) * 2018-05-02 2019-11-15 杭州海康威视数字技术股份有限公司 一种图像清晰度评测方法、装置及电子设备
CN110853006A (zh) * 2019-11-05 2020-02-28 华南理工大学 一种使用扫描仪获取的数字病理图像质量评价的方法
CN111429437A (zh) * 2020-03-30 2020-07-17 南京工程学院 面向目标检测的图像无参考清晰度质量检测方法
CN111724054A (zh) * 2020-06-16 2020-09-29 广州视源电子科技股份有限公司 一种评估方法、装置、终端设备及存储介质
CN112258467A (zh) * 2020-10-19 2021-01-22 浙江大华技术股份有限公司 一种图像清晰度的检测方法和装置以及存储介质
CN112381820A (zh) * 2020-12-07 2021-02-19 深圳市福日中诺电子科技有限公司 一种基于一组同场景照片清晰度的评价方法
WO2021057463A1 (zh) * 2019-09-25 2021-04-01 北京字节跳动网络技术有限公司 图像风格化处理方法、装置、电子设备及可读介质
CN112750121A (zh) * 2021-01-20 2021-05-04 赛维森(广州)医疗科技服务有限公司 用于病理玻片数字图像质量检测的系统及其方法
CN113487508A (zh) * 2021-07-08 2021-10-08 山东志盈医学科技有限公司 数字切片扫描仪的图片清晰度动态调整的方法和装置
CN113822877A (zh) * 2021-11-17 2021-12-21 武汉中导光电设备有限公司 Aoi设备显微镜缺陷检测图片质量评价方法及系统
CN114723721A (zh) * 2022-04-18 2022-07-08 贝塔通科技(北京)有限公司 图像清晰度确定方法、系统及图像上传方法
CN117078565A (zh) * 2023-10-17 2023-11-17 深圳市精研科洁科技股份有限公司 一种摄像头抖动模糊图像优化增强方法及系统
CN117218452A (zh) * 2023-11-02 2023-12-12 临沂市兰山区自然资源开发服务中心 一种土地图像自动分类管理系统
CN117237747A (zh) * 2023-11-14 2023-12-15 深圳市明鸿五金制品有限公司 基于人工智能的五金件缺陷分类识别方法
CN117974464A (zh) * 2024-02-28 2024-05-03 上海为旌科技有限公司 图像fv置信度确定方法、装置、设备及存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102800082A (zh) * 2012-06-19 2012-11-28 浙江捷尚视觉科技有限公司 一种无参考的图像清晰度检测方法
CN103093419A (zh) * 2011-10-28 2013-05-08 浙江大华技术股份有限公司 一种检测图像清晰度的方法及装置
CN104637064A (zh) * 2015-02-28 2015-05-20 中国科学院光电技术研究所 一种基于边缘强度权重的离焦模糊图像清晰度检测方法
US9563963B2 (en) * 2013-03-11 2017-02-07 Reeves Wireline Technologies Limited Methods of and apparatuses for identifying geological characteristics in boreholes
CN106530281A (zh) * 2016-10-18 2017-03-22 国网山东省电力公司电力科学研究院 基于边缘特征的无人机图像模糊判断方法及系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103093419A (zh) * 2011-10-28 2013-05-08 浙江大华技术股份有限公司 一种检测图像清晰度的方法及装置
CN102800082A (zh) * 2012-06-19 2012-11-28 浙江捷尚视觉科技有限公司 一种无参考的图像清晰度检测方法
US9563963B2 (en) * 2013-03-11 2017-02-07 Reeves Wireline Technologies Limited Methods of and apparatuses for identifying geological characteristics in boreholes
CN104637064A (zh) * 2015-02-28 2015-05-20 中国科学院光电技术研究所 一种基于边缘强度权重的离焦模糊图像清晰度检测方法
CN106530281A (zh) * 2016-10-18 2017-03-22 国网山东省电力公司电力科学研究院 基于边缘特征的无人机图像模糊判断方法及系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
刘亚梅: "基于梯度边缘最大值的图像清晰度评价", 《图学学报》 *
张天煜 等: "基于强边缘宽度直方图的图像清晰度指标", 《浙江大学学报(工学版)》 *
钱青 等: "一种改进的Sobel算子图像清晰度评价函数", 《计算机与数字工程》 *

Cited By (27)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108460766B (zh) * 2018-04-12 2022-02-25 四川和生视界医药技术开发有限公司 一种视网膜图像清晰度评估方法以及评估装置
CN108460766A (zh) * 2018-04-12 2018-08-28 四川和生视界医药技术开发有限公司 一种视网膜图像清晰度评估方法以及评估装置
CN110458789B (zh) * 2018-05-02 2022-04-05 杭州海康威视数字技术股份有限公司 一种图像清晰度评测方法、装置及电子设备
CN110458789A (zh) * 2018-05-02 2019-11-15 杭州海康威视数字技术股份有限公司 一种图像清晰度评测方法、装置及电子设备
CN109191457A (zh) * 2018-09-21 2019-01-11 中国人民解放军总医院 一种病理图像质量有效性识别方法
CN109191457B (zh) * 2018-09-21 2022-07-01 中国人民解放军总医院 一种病理图像质量有效性识别方法
CN110148147A (zh) * 2018-11-07 2019-08-20 腾讯大地通途(北京)科技有限公司 图像检测方法、装置、存储介质和电子装置
CN110148147B (zh) * 2018-11-07 2024-02-09 腾讯大地通途(北京)科技有限公司 图像检测方法、装置、存储介质和电子装置
WO2021057463A1 (zh) * 2019-09-25 2021-04-01 北京字节跳动网络技术有限公司 图像风格化处理方法、装置、电子设备及可读介质
CN110853006A (zh) * 2019-11-05 2020-02-28 华南理工大学 一种使用扫描仪获取的数字病理图像质量评价的方法
CN111429437A (zh) * 2020-03-30 2020-07-17 南京工程学院 面向目标检测的图像无参考清晰度质量检测方法
CN111724054A (zh) * 2020-06-16 2020-09-29 广州视源电子科技股份有限公司 一种评估方法、装置、终端设备及存储介质
CN112258467B (zh) * 2020-10-19 2024-06-18 浙江大华技术股份有限公司 一种图像清晰度的检测方法和装置以及存储介质
CN112258467A (zh) * 2020-10-19 2021-01-22 浙江大华技术股份有限公司 一种图像清晰度的检测方法和装置以及存储介质
CN112381820A (zh) * 2020-12-07 2021-02-19 深圳市福日中诺电子科技有限公司 一种基于一组同场景照片清晰度的评价方法
CN112750121B (zh) * 2021-01-20 2021-11-26 赛维森(广州)医疗科技服务有限公司 用于病理玻片数字图像质量检测的系统及其方法
CN112750121A (zh) * 2021-01-20 2021-05-04 赛维森(广州)医疗科技服务有限公司 用于病理玻片数字图像质量检测的系统及其方法
CN113487508A (zh) * 2021-07-08 2021-10-08 山东志盈医学科技有限公司 数字切片扫描仪的图片清晰度动态调整的方法和装置
CN113487508B (zh) * 2021-07-08 2024-03-26 山东志盈医学科技有限公司 数字切片扫描仪的图片清晰度动态调整的方法和装置
CN113822877A (zh) * 2021-11-17 2021-12-21 武汉中导光电设备有限公司 Aoi设备显微镜缺陷检测图片质量评价方法及系统
CN114723721A (zh) * 2022-04-18 2022-07-08 贝塔通科技(北京)有限公司 图像清晰度确定方法、系统及图像上传方法
CN117078565A (zh) * 2023-10-17 2023-11-17 深圳市精研科洁科技股份有限公司 一种摄像头抖动模糊图像优化增强方法及系统
CN117078565B (zh) * 2023-10-17 2024-02-02 深圳市精研科洁科技股份有限公司 一种摄像头抖动模糊图像优化增强方法及系统
CN117218452A (zh) * 2023-11-02 2023-12-12 临沂市兰山区自然资源开发服务中心 一种土地图像自动分类管理系统
CN117218452B (zh) * 2023-11-02 2024-02-06 临沂市兰山区自然资源开发服务中心 一种土地图像自动分类管理系统
CN117237747A (zh) * 2023-11-14 2023-12-15 深圳市明鸿五金制品有限公司 基于人工智能的五金件缺陷分类识别方法
CN117974464A (zh) * 2024-02-28 2024-05-03 上海为旌科技有限公司 图像fv置信度确定方法、装置、设备及存储介质

Also Published As

Publication number Publication date
CN107507173B (zh) 2021-07-27

Similar Documents

Publication Publication Date Title
CN107507173A (zh) 一种全切片图像的无参考清晰度评估方法及系统
Zhang et al. Color correction and adaptive contrast enhancement for underwater image enhancement
CN111583227B (zh) 一种荧光细胞自动计数方法、装置、设备和介质
CN108805023B (zh) 一种图像检测方法、装置、计算机设备及存储介质
CN105096347B (zh) 图像处理装置和方法
CN111915704A (zh) 一种基于深度学习的苹果分级识别方法
CN108537787B (zh) 一种人脸图像的质量判定方法
Kanwal et al. Region based adaptive contrast enhancement of medical X-ray images
Prasath et al. Multiscale structure tensor for improved feature extraction and image regularization
Sadanandan et al. Segmentation and track-analysis in time-lapse imaging of bacteria
CN106875391A (zh) 皮肤图像的识别方法及电子设备
CN112991374B (zh) 基于Canny算法的边缘增强方法、装置、设备及存储介质
CN106934806A (zh) 一种基于结构清晰度的无参考图失焦模糊区域分割方法
CN107506795A (zh) 一种面向图像匹配的局部灰度直方图特征描述子建立方法和图像匹配方法
CN117274113B (zh) 基于图像增强的碎硅片清洗效果视觉检测方法
CN116363140B (zh) 一种中硼硅玻璃缺陷的检测方法、系统、装置及存储介质
CN112330613B (zh) 一种细胞病理数字图像质量的评价方法及系统
CN107622475A (zh) 图像拼接中的灰度校正方法
CN105550694A (zh) 一种度量人脸图像模糊程度的方法
Soni et al. CT scan based brain tumor recognition and extraction using Prewitt and morphological dilation
WO2019205603A1 (zh) 图像模糊度检测方法、装置、计算机设备及可读存储介质
CN111553927B (zh) 棋盘格角点检测方法、检测系统、计算机装置及存储介质
CN115131336B (zh) 一种显示屏四色背景画面下的暗线缺陷检测方法
CN109448012A (zh) 一种图像边缘检测方法及装置
CN106803240B (zh) 一种工业图像均光处理方法

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
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20220706

Address after: 201100 rooms 615 and 616, 6 / F, building 2, No. 322, Lane 953, Jianchuan Road, Minhang District, Shanghai

Patentee after: Shanghai Pu Huasen Biotechnology Co.,Ltd.

Address before: 200240 No. 800, Dongchuan Road, Shanghai, Minhang District

Patentee before: SHANGHAI JIAO TONG University

TR01 Transfer of patent right