CN113256573B - 判定数字切片质量的方法和装置 - Google Patents

判定数字切片质量的方法和装置 Download PDF

Info

Publication number
CN113256573B
CN113256573B CN202110520607.2A CN202110520607A CN113256573B CN 113256573 B CN113256573 B CN 113256573B CN 202110520607 A CN202110520607 A CN 202110520607A CN 113256573 B CN113256573 B CN 113256573B
Authority
CN
China
Prior art keywords
image
view
slice
visual field
focusing
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.)
Active
Application number
CN202110520607.2A
Other languages
English (en)
Other versions
CN113256573A (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.)
Shandong Zhiying Medical Technology Co ltd
Original Assignee
Shandong Zhiying Medical Technology Co ltd
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 Shandong Zhiying Medical Technology Co ltd filed Critical Shandong Zhiying Medical Technology Co ltd
Priority to CN202110520607.2A priority Critical patent/CN113256573B/zh
Publication of CN113256573A publication Critical patent/CN113256573A/zh
Application granted granted Critical
Publication of CN113256573B publication Critical patent/CN113256573B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • 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
    • G06T7/00—Image analysis
    • G06T7/90—Determination of colour characteristics
    • 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)
  • Computer Vision & Pattern Recognition (AREA)
  • Quality & Reliability (AREA)
  • Image Analysis (AREA)
  • Color Image Communication Systems (AREA)

Abstract

本发明公开了一种判定数字切片质量的方法和装置，方法包括以下步骤：步骤1，从切片整体概貌图中获取只包含样本区域的扫描预览图；在扫描视野中按照一定间隔选取聚焦视野并放入焦点视野列表中；步骤2，在样本区域外采集图像并取反得到增强图像,在采集视野图像后与增强图像进行加运算，计算增强图像的色调、饱和度、颜色均值和清晰度；步骤3，计算聚焦点平均信息变量：聚焦点平均色调、聚焦点平均饱和度、聚焦点平均颜色均值和聚焦点平均清晰度；步骤4，定义切片总分变量，对聚焦视野进行评分。本发明不仅能够判断数字切片是否满足应用的要求，而且还可以对需要重新扫描的切片进行指导参数调整。

Description

判定数字切片质量的方法和装置
技术领域
本发明涉及一种判定数字切片质量的方法和装置，属于数字切片图像评价技术领域。
背景技术
利用数字切片扫描仪可以将传统切片数字化。数字化的过程主要是利用预览相机生成预览图片、从预览图片中识别组织区域、计算出需要所有需要扫描的视野、产生增强图片(用于扫描视野的增强)、从所有视野中选取聚焦视野进行聚焦、将聚焦视野的聚焦位置扩展到所有的扫描视野中、利用扫描摄像机采集图像、对图像进行拼接最终形成分辨率很大的数字切片。
在扫描过程中，由于不是对每个扫描视野进行聚焦，而且切片的厚度没法保持一致性，就会照成非聚焦视野不清楚的情况；而且由于机械运动和外界环境的影响，设备会存在微弱震动的情况，扫描仪经过放大光路后，微弱震动也会对清晰度产生影响；另外由于切片的染色问题，有的部分染色深、有的部分染色浅，造成相机用特定色彩方案进行扫描时，图片的色彩饱和度存在较大的差异。总之，如果较多的视野出现清晰度和色彩饱和度较差的情况，则整张切片的可用性就会很差，需要调整设备参数进行重新扫描。重新扫描可以增加聚焦视野数量并选用不同的色彩方案。
为了判断一张切片的质量，需要对整张数字切片进行评价。而数字切片是由一幅幅图像拼接而成，因此对每个扫描视野单独计算其分值后进行平均即可。这就需要对一个视野的一副图像进行计算分值；如果每个视野中全部是组织区域，按照一个标准进行计算即可。但是在实际扫描过程中，由于预览图与实际扫描位置的倍率差异、预览图进行了扩展、空洞填充等操作，会造成一副视野图片中存在一部分有组织、一部分没有组织，特别是在边界时，极限情况下会存在扫描视野不存在任何组织区域，为空白视野。这样都给计算增加了难度。
发明内容
为了解决上述问题，本发明提出了一种判定数字切片质量的方法和装置。能完成数字切片每个扫描视野的分值计算及整个切片的评价，帮助判断数字切片是否符合应用要求，及重新扫描时参数的选择。
本发明解决其技术问题采取的技术方案是：
第一方面，本发明实施例提供的一种判定数字切片质量的方法，包括以下步骤：
步骤1，从切片整体概貌图中获取只包含样本区域的扫描预览图；根据扫描预览图得到所有扫描视野，在扫描视野中按照一定间隔选取聚焦视野并放入焦点视野列表中；
步骤2，在样本区域外采集图像并取反得到增强图像,在采集视野图像后与增强图像进行加运算，计算增强图像的色调、饱和度、颜色均值和清晰度；
步骤3，计算聚焦点平均信息变量：聚焦点平均色调、聚焦点平均饱和度、聚焦点平均颜色均值和聚焦点平均清晰度；
步骤4，定义切片总分变量，对聚焦视野进行评分。
作为本实施例一种可能的实现方式，所述步骤1具体包括以下步骤：
将切片放入切片扫描仪，利用预览相机产生切片整体概貌图；
对概貌图进行灰度化、二值化、斑点去除和空洞填充处理，从整体概貌图中得到只包含样本区域的扫描预览图；
从扫描预览图中按照视野与预览图的对应关系产生扫描视野列表，在扫描视野中按照一定间隔选取聚焦视野，放入焦点视野列表。
作为本实施例一种可能的实现方式，所述步骤2具体包括以下步骤：
在样本区域外采集图像并取反得到增强图像,在采集视野图像后与增强图像进行加运算处理；
计算增强图像的色调、饱和度、颜色均值、清晰度。
作为本实施例一种可能的实现方式，所述步骤3具体包括以下步骤：
对焦点视野列表中每个视野进行聚焦，得到焦点图像，并计算每个焦点图像的色调、饱和度、颜色均值和清晰度；
计算聚焦点平均信息变量。
作为本实施例一种可能的实现方式，所述步骤4具体为：
采集扫描视野列表中每个扫描视野的视野图像；计算每个视野图像的色调、饱和度、颜色均值和清晰度；
计算视野信息距离值：扫描视野信息与增强图信息距离值、扫描视野信息与聚焦点平均信息距离值；
定义评分所用信息变量并计算每个扫描视野的视野分值；
计算切片分值。
作为本实施例一种可能的实现方式，所述的判定数字切片质量的方法还包括以下步骤：
步骤5，将切片分数保存到切片总分变量，如果切片分数大于设定值则保存扫描的切片图像，否则重新扫描。
第二方面，本发明实施例提供的一种判定数字切片质量的装置，包括：
焦点视野列表获取模块，用于从切片整体概貌图中获取只包含样本区域的扫描预览图；在扫描视野中按照一定间隔选取聚焦视野并放入焦点视野列表中；
图像增强模块，用于在样本区域外采集图像并取反得到增强图像,在采集视野图像后与增强图像进行加运算，计算增强图像的色调、饱和度、颜色均值和清晰度；
聚焦点变量计算模块，用于计算聚焦点平均信息变量：聚焦点平均色调、聚焦点平均饱和度、聚焦点平均颜色均值和聚焦点平均清晰度；
聚焦视野评分模块，用于定义切片总分变量，对聚焦视野进行评分。
作为本实施例一种可能的实现方式，所述聚焦视野选取模块包括：
概貌图获取模块，用于将切片放入切片扫描仪，利用预览相机产生切片整体概貌图；
概貌图处理模块，用于对概貌图进行灰度化、二值化、斑点去除和空洞填充处理，从整体概貌图中得到只包含样本区域的扫描预览图；
聚焦视野选取模块，用于从扫描预览图中按照视野与预览图的对应关系产生扫描视野列表listScanView，在扫描视野中按照一定间隔选取聚焦视野(不选边缘的扫描视野)，放入焦点视野列表listFocusView。
作为本实施例一种可能的实现方式，所述图像增强模块包括：
视野图像处理模块，用于在样本区域外采集图像并取反得到增强图像,在采集视野图像后与增强图像进行加运算处理；
增强图像变量计算模块，用于计算增强图像的的色调、饱和度、颜色均值、清晰度。
作为本实施例一种可能的实现方式，所述聚焦点变量计算模块具体用于：
对焦点视野列表中每个视野进行聚焦，得到焦点图像，并计算每个焦点图像的色调、饱和度、颜色均值和清晰度；
计算聚焦点平均信息变量。
作为本实施例一种可能的实现方式，所述聚焦视野评分模块具体用于：
采集扫描视野列表中每个扫描视野的视野图像；计算每个视野图像的色调、饱和度、颜色均值和清晰度；
计算视野信息距离值：扫描视野信息与增强图信息距离值、扫描视野信息与聚焦点平均信息距离值；
定义评分所用信息变量并计算每个扫描视野的视野分值；
计算切片分值。
作为本实施例一种可能的实现方式，所述的判定数字切片质量的装置还包括：切片图像扫描模块，用于将切片分数保存到切片总分变量，如果切片分数大于设定值则保存扫描的切片图像，否则重新扫描。
本发明实施例的技术方案可以具有的有益效果如下：
本发明实施例的技术方案的一种判定数字切片质量的方法包括以下步骤：骤1，从切片整体概貌图中获取只包含样本区域的扫描预览图；在扫描视野中按照一定间隔选取聚焦视野并放入焦点视野列表中；步骤2，在样本区域外采集图像并取反得到增强图像,在采集视野图像后与增强图像进行加运算，计算增强图像的色调、饱和度、颜色均值和清晰度；步骤3，计算聚焦点平均信息变量：聚焦点平均色调、聚焦点平均饱和度、聚焦点平均颜色均值和聚焦点平均清晰度；步骤4，定义切片总分变量，对聚焦视野进行评分。本发明能完成数字切片每个扫描视野的分值计算及整个切片的评分，帮助判断数字切片是否符合应用要求，及重新扫描时参数的选择。
本发明实施例的技术方案的一种判定数字切片质量的装置，不仅能够判断数字切片是否满足应用的要求，而且还可以对需要重新扫描的切片进行指导参数调整。
附图说明：
图1是根据一示例性实施例示出的一种判定数字切片质量的方法流程图；
图2是根据一示例性实施例示出的一种判定数字切片质量的装置的结构图。
具体实施方式
下面结合附图与实施例对本发明做进一步说明：
为能清楚说明本方案的技术特点，下面通过具体实施方式，并结合其附图，对本发明进行详细阐述。下文的公开提供了许多不同的实施例或例子用来实现本发明的不同结构。为了简化本发明的公开，下文中对特定例子的部件和设置进行描述。此外，本发明可以在不同例子中重复参考数字和/或字母。这种重复是为了简化和清楚的目的，其本身不指示所讨论各种实施例和/或设置之间的关系。应当注意，在附图中所图示的部件不一定按比例绘制。本发明省略了对公知组件和处理技术及工艺的描述以避免不必要地限制本发明。
图1是根据一示例性实施例示出的一种判定数字切片质量的方法流程图。如图1所示，本发明实施例提供的一种判定数字切片质量的方法，包括以下步骤：
步骤1，从切片整体概貌图中获取只包含样本区域的扫描预览图；根据扫描预览图得到所有扫描视野，在扫描视野中按照一定间隔选取聚焦视野并放入焦点视野列表中；
步骤2，在样本区域外采集图像并取反得到增强图像,在采集视野图像后与增强图像进行加运算，计算增强图像的色调、饱和度、颜色均值和清晰度；
步骤3，计算聚焦点平均信息变量：聚焦点平均色调、聚焦点平均饱和度、聚焦点平均颜色均值和聚焦点平均清晰度；
步骤4，定义切片总分变量，对聚焦视野进行评分。
作为本实施例一种可能的实现方式，所述步骤1具体包括以下步骤：
将切片放入切片扫描仪，利用预览相机产生切片整体概貌图；
对概貌图进行灰度化、二值化、斑点去除和空洞填充处理，从整体概貌图中得到只包含样本区域的扫描预览图；
从扫描预览图中按照视野与预览图的对应关系产生扫描视野列表listScanView，在扫描视野中按照一定间隔选取聚焦视野(不选边缘的扫描视野)，放入焦点视野列表listFocusView。
作为本实施例一种可能的实现方式，所述步骤2具体为：
在样本区域外采集图像Enhance并取反得到增强图像viewEnhance,在采集视野图像后与增强图像进行加运算处理；
计算增强图像的viewEnhance的色调hueEnhance、饱和度saturationEnhance、颜色均值colorEnhance、清晰度clarityEnhance。
作为本实施例一种可能的实现方式，所述步骤3具体为：
对焦点视野列表中每个视野进行聚焦，得到焦点图像，并计算每个焦点图像的色调、饱和度、颜色均值和清晰度；
计算聚焦点平均信息变量。
步骤3具体包括以下步骤：
步骤3.1，设置i＝0；对焦点视野列表listFocusView中第i个视野viewFocus进行聚焦，得到焦点图像pictureFocus，并计算焦点图像pictureFocus的色调hueFocus、饱和度saturationFocus、颜色均值colorFocus、清晰度clarityFocus；
计算聚焦点平均信息色调：
hueAvgFocus＝hueAvgFocus*(i/(i+1.0))+hueFocus*(1/(i+1.0))；
饱和度：
saturationAvgFocus＝saturationAvgFocus*(i/(i+1.0))+saturationFocus*(1/(i+1.0))；
颜色均值：
colorAvgFocus＝colorAvgFocus*(i/(i+1.0))+colorFocus*(1/(i+1.0)；
清晰度：
clarityAvgFocus＝clarityAvgFocus*(i/(i+1.0))+clarityFocus*(1/(i+1.0)；
步骤3.2，i＝i+1；如果i小于等于listFocusView的数量，则转到步骤3.1；如果i大于listFocusView的数量，则转到下一步。
作为本实施例一种可能的实现方式，所述步骤4具体为：
采集扫描视野列表中每个扫描视野的视野图像；计算每个视野图像的色调、饱和度、颜色均值和清晰度；
计算视野信息距离值：扫描视野信息与增强图信息距离值、扫描视野信息与聚焦点平均信息距离值；
定义评分所用信息变量并计算每个扫描视野的视野分值；
计算切片分值。
步骤4具体包括以下步骤：
步骤4.1，设置m＝0；从扫描视野列表listScanView中得到第m个扫描视野viewScan；采集viewScan位置的视野图像pictureScan；计算pictureScan的色调hueScan、饱和度saturationScan、颜色均值colorScan、清晰度clarityScan。计算视野信息距离值；
扫描视野信息与增强图信息距离值：enhanceDistance＝0.2*abs(hueScan-hueEnhance)+0.3*abs(saturationScan-saturationEnhance)+0.1*abs(colorScan-colorEnhance)+0.4*abs(clarityScan–clarityEnhance)；
扫描视野信息与聚焦点平均信息距离值：focusDistance＝0.2*abs(hueScan-hueAvgFocus)+0.3*abs(saturationScan-saturationAvgFocus)+0.1*abs(colorScan-colorAvgFocus)+0.4*abs(clarityScan–clarityAvgFocus)；在四个信息判断中采用不同的权重，主要因为在清晰度和饱和度是主要的关注点；色调和颜色均值不是主要关注点且在不同图像中这两个信息差异不会非常明显。
步骤4.2定义评分所用信息变量：色调hueUsed、饱和度saturationUsed、颜色均值colorUsed、清晰度clarityUsed；
如果enhanceDistance小于focusDistance，则hueUsed＝hueEnhance、saturationUsed＝saturationEnhance、colorUsed＝colorEnhance、clarityUsed＝clarityEnhance；
如果enhanceDistance大于等于focusDistance，则hueUsed＝hueAvgFocus、saturationUsed＝saturationAvgFocus、colorUsed＝colorAvgFocus、clarityUsed＝clarityAvgFocus；
步骤4.3，对扫描视野进行评分：
色调分值：hueScore＝1-abs(hueScan-hueUsed)/((hueScan+hueUsed)/2)；
饱和度分值：
saturationScore＝1-abs(saturationScan-saturationUsed)/((saturationScan+saturationUsed)/2)；
颜色分值：colorScore＝1-abs(colorScan-colorUsed)/((colorScan+colorUsed)/2)；
clarityScore＝1-abs(clarityScan-clarityUsed)/((clarityScan+clarityUsed)/2)；
如果hueScore、saturationScore、colorScore和clarityScore中如果存在小于0.6的数值，则将视野位置信息加入到listLowView中；
视野分值viewScore＝(0.1*hueScore+0.3*saturationScore+0.1*colorScore+0.4*clarityScore)*100；
步骤4.4，切片分值为：slideScore＝slideScore*(m/(m+1.0))+viewScore*(1.0/(m+1.0))；
步骤4.5，m＝m+1；如果m小于等于listScanView的数量，则转到4.1；如果i大于listScanView的数量，则转到下一步。
作为本实施例一种可能的实现方式，所述的判定数字切片质量的方法还包括以下步骤：
步骤5，将切片分数保存到切片总分变量，如果切片分数大于设定值则保存扫描的切片图像，否则重新扫描。
如图2所示，本发明实施例提供的一种判定数字切片质量的装置，包括：
焦点视野列表获取模块，用于从切片整体概貌图中获取只包含样本区域的扫描预览图；根据扫描预览图得到所有扫描视野，在扫描视野中按照一定间隔选取聚焦视野并放入焦点视野列表中；
图像增强模块，用于在样本区域外采集图像并取反得到增强图像,在采集视野图像后与增强图像进行加运算，计算增强图像的色调、饱和度、颜色均值和清晰度；
聚焦点变量计算模块，用于计算聚焦点平均信息变量：聚焦点平均色调、聚焦点平均饱和度、聚焦点平均颜色均值和聚焦点平均清晰度；
聚焦视野评分模块，用于定义切片总分变量，对聚焦视野进行评分。
作为本实施例一种可能的实现方式，所述聚焦视野选取模块包括：
概貌图获取模块，用于将切片放入切片扫描仪，利用预览相机产生切片整体概貌图；
概貌图处理模块，用于对概貌图进行灰度化、二值化、斑点去除和空洞填充处理，从整体概貌图中得到只包含样本区域的扫描预览图；
聚焦视野选取模块，用于从扫描预览图中按照视野与预览图的对应关系产生扫描视野列表listScanView，在扫描视野中按照一定间隔选取聚焦视野(不选边缘的扫描视野)，放入焦点视野列表listFocusView。
作为本实施例一种可能的实现方式，所述图像增强模块包括：
视野图像处理模块，用于在样本区域外采集图像并取反得到增强图像,在采集视野图像后与增强图像进行加运算处理；
增强图像变量计算模块，用于计算增强图像的的色调、饱和度、颜色均值、清晰度。
作为本实施例一种可能的实现方式，所述聚焦点变量计算模块具体用于：
对焦点视野列表中每个视野进行聚焦，得到焦点图像，并计算每个焦点图像的色调、饱和度、颜色均值和清晰度；
计算聚焦点平均信息变量。
作为本实施例一种可能的实现方式，所述聚焦视野评分模块具体用于：
采集扫描视野列表中每个扫描视野的视野图像；计算每个视野图像的色调、饱和度、颜色均值和清晰度；
计算视野信息距离值：扫描视野信息与增强图信息距离值、扫描视野信息与聚焦点平均信息距离值；
定义评分所用信息变量并计算每个扫描视野的视野分值；
计算切片分值。
作为本实施例一种可能的实现方式，所述的判定数字切片质量的装置还包括：切片图像扫描模块，用于将切片分数保存到切片总分变量，如果切片分数大于设定值则保存扫描的切片图像，否则重新扫描。
采用本发明实施例的一种判定数字切片质量的装置对数字切片中的图像进行评分判定的具体过程如下。
1、将切片放入切片扫描仪，利用预览相机产生切片整体概貌图，对概貌图进行灰度化、二值化、斑点去除、空洞填充等算法，从整体概貌图中得到只包含样本区域的扫描预览图；从扫描预览图中按照视野与预览图的对应关系产生扫描视野列表listScanView；在扫描视野中按照一定间隔选取聚焦视野(不选边缘的扫描视野)，放入列表listFocusView。
2、更加二值化以后的样本预览图，在样本区域外采集图像Enhance并取反得到增强图像viewEnhance。
计算增强图的viewEnhance的色调hueEnhance＝85.39、饱和度saturationEnhance＝3.36、颜色均值colorEnhance＝245.19、清晰度clarityEnhance＝3.55。
3、定义聚焦点平均信息变量：
色调hueAvgFocus＝0、饱和度saturationAvgFocus＝0、颜色均值colorAvgFocus＝0、清晰度clarityAvgFocus＝0。设置i＝0；
3.1从listFocusView得到第i个视野viewFocus,对viewFocus进行聚焦，得到焦点图像pictureFocus，计算pictureFocus的色调hueFocus＝149.47、饱和度saturationFocus＝74.77、颜色均值colorFocus＝150.04、清晰度clarityFocus＝37.67；
计算聚焦点平均信息色调：
hueAvgFocus＝hueAvgFocus*(i/(i+1.0))+hueFocus*(1/(i+1.0))＝0*(0/(0+1.0))+149.47*(1/(0+1.0))＝149.47；
利用类似公式可以求出饱和度saturationAvgFocus＝74.77、颜色均值colorAvgFocus＝150.04、清晰度clarityAvgFocus＝37.67；
3.2i＝i+1＝1,如果i小于等于listFocusView.size()，则转到步骤3.1；如果i大于listFocusView.size()，则转到下一步(即步骤4)。
经过计算hueAvgFocus＝150.32、saturationAvgFocus＝75.77、colorAvgFocus＝149.56、clarityAvgFocus＝36.99；
4、定义切片总分变量slideScore。设置m＝0；
4.1，从扫描视野列表listScanView中得到第m个扫描视野viewScan；采集viewScan位置的视野图像pictureScan；计算pictureScan的色调hueScan＝126.57、饱和度saturationScan＝49.58、颜色均值colorScan＝188.97、清晰度clarityScan＝53.32。
计算视野信息距离值；
扫描视野信息与增强图信息距离值：
enhanceDistance＝0.2*abs(hueScan-hueEnhance)+0.3*abs(saturationScan-saturationEnhance)+0.1*abs(colorScan-colorEnhance)+0.4*abs(clarityScan–clarityEnhance)＝0.2*abs(126.57-85.39)+0.3*abs(49.58-3.36)+0.1*abs(188.97-245.19)+0.4*abs(53.32–3.55)＝47.63；
扫描视野信息与聚焦点平均信息距离值：
focusDistance＝0.2*abs(hueScan-hueAvgFocus)+0.3*abs(saturationScan-saturationAvgFocus)+0.1*abs(colorScan-colorAvgFocus)+0.4*abs(clarityScan–clarityAvgFocus)＝0.2*abs(126.57-150.32)+0.3*abs(49.58-75.77)+0.1*abs(188.97-149.56)+0.4*abs(53.32–36.99)＝23.08。
4.2，定义评分所用信息变量：色调hueUsed、饱和度saturationUsed、颜色均值colorUsed、清晰度clarityUsed；
如果enhanceDistance>focusDistance，则hueUsed＝hueAvgFocus＝150.32、saturationUsed＝saturationAvgFocus＝75.77、colorUsed＝colorAvgFocus＝149.56、clarityUsed＝clarityAvgFocus＝36.99。
4.3，对扫描视野进行评分：
色调分值：
hueScore＝1-abs(hueScan-hueUsed)/((hueScan+hueUsed)/2)＝1-abs(126.57-150.32)/((126.57+150.32)/2)＝0.828；
饱和度分值：
saturationScore＝1-abs(saturationScan-saturationUsed)/((saturationScan+saturationUsed)/2)＝1-abs(49.58-75.77)/((49.58+75.77)/2)＝0.582；
颜色分值：
colorScore＝1-abs(colorScan-colorUsed)/((colorScan+colorUsed)/2)＝1-abs(188.97-149.56)/((188.97+149.56)/2)＝0.767；
clarityScore＝1-abs(clarityScan-clarityUsed)/((clarityScan+clarityUsed)/2)＝1-abs(53.32-36.99)/((53.32+36.99)/2)＝0.638。
因为saturationScare<0.6，将视野信息加入到listLowView中。
视野分值:
viewScore＝(0.2*hueScore+0.3*saturationScore+0.1*colorScore+0.4*clarityScore)*100＝(0.2*0.828+0.3*0.582+0.1*0.767+0.4*0.638)*100＝67.21；
4.4，切片分值为：
slideScore＝slideScore*(m/(m+1.0))+viewScore*(1.0/(m+1.0))＝0*(0/(0+1.0))+67.21*(1.0/(0+1.0))＝67.21；
4.5，m＝m+1；如果m小于等于listScanView的数量，则转到步骤4.1；如果i大于listScanView的数量，则转到下一步(即步骤5)。
5、经过步骤4的整个过程后，slideScore＝73.67；扫描完成将切片分数保存到切片文件中slideScore。分数大于60则扫描通过。
上述计算公式中的Abs()为求绝对值函数。
最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本发明的权利要求保护范围之内。

Claims (7)

1.一种判定数字切片质量的方法，其特征是，包括以下步骤：
步骤1，从切片整体概貌图中获取只包含样本区域的扫描预览图；在扫描视野中按照一定间隔选取聚焦视野并放入焦点视野列表中；
步骤2，在样本区域外采集图像并取反得到增强图像,在采集视野图像后与增强图像进行加运算，计算增强图像的色调、饱和度、颜色均值和清晰度；
步骤3，对聚焦视野列表中的每个视野进行聚焦，得到焦点图像，计算焦点图像的色调、饱和度、颜色均值、清晰度；求所有焦点图像的平均色调、平均饱和度、平均颜色均值、平均清晰度；
步骤4，根据视野分值计算切片分值；
所述步骤1具体包括以下步骤：
将切片放入切片扫描仪，利用预览相机产生切片整体概貌图；
对概貌图进行灰度化、二值化、斑点去除和空洞填充处理，从整体概貌图中得到只包含样本区域的扫描预览图；
从扫描预览图中按照视野与预览图的对应关系产生扫描视野列表，在扫描视野中按照一定间隔选取聚焦视野，放入焦点视野列表；
所述步骤4具体为：
采集扫描视野列表中每个扫描视野的视野图像；计算每个视野图像的色调、饱和度、颜色均值和清晰度；
计算视野信息距离值：比较扫描视野信息与增强图信息距离值和扫描视野
信息与聚焦点平均信息距离值；
选择评分所用视野信息的距离值，并计算每个扫描视野的视野分值；
根据每个扫描视野的视野分值计算切片分值。
2.根据权利要求1所述的判定数字切片质量的方法，其特征是，所述步骤2具体包括以下步骤：
在样本区域外采集图像并取反得到增强图像,在采集视野图像后与增强图像进行加运算处理；
计算增强图像的色调、饱和度、颜色均值、清晰度。
3.根据权利要求2所述的判定数字切片质量的方法，其特征是，所述步骤3具体包括以下步骤：
对焦点视野列表中每个视野进行聚焦，得到焦点图像，并计算每个焦点图像的色调、饱和度、颜色均值和清晰度；
计算聚焦点平均信息变量。
4.根据权利要求1-3任意一项所述的判定数字切片质量的方法，其特征是，还包括以下步骤：
步骤5，将切片分数保存到切片总分变量，如果切片分数大于设定值则保存扫描的切片图像，否则重新扫描。
5.一种判定数字切片质量的装置，其特征是，包括：
焦点视野列表获取模块，用于从切片整体概貌图中获取只包含样本区域的扫描预览图；在扫描视野中按照一定间隔选取聚焦视野并放入焦点视野列表中；
图像增强模块，用于在样本区域外采集图像并取反得到增强图像,在采集视野图像后与增强图像进行加运算，计算增强图像的色调、饱和度、颜色均值和清晰度；
聚焦点变量计算模块，用于对聚焦视野列表中的每个视野进行聚焦，得到焦点图像，计算焦点图像的色调、饱和度、颜色均值、清晰度；求所有焦点图像的平均色调、平均饱和度、平均颜色均值、平均清晰度；
聚焦视野评分模块，用于根据视野分值计算切片分值；
所述聚焦视野选取模块包括：
概貌图获取模块，用于将切片放入切片扫描仪，利用预览相机产生切片整体概貌图；
概貌图处理模块，用于对概貌图进行灰度化、二值化、斑点去除和空洞填充处理，从整体概貌图中得到只包含样本区域的扫描预览图；
聚焦视野选取模块，用于从扫描预览图中按照视野与预览图的对应关系产生扫描视野列表，在扫描视野中按照一定间隔选取聚焦视野，放入焦点视野列表中；
所述聚焦视野评分模块具体用于：
采集扫描视野列表中每个扫描视野的视野图像；计算每个视野图像的色调、饱和度、颜色均值和清晰度；
计算视野信息距离值：比较扫描视野信息与增强图信息距离值和扫描视野
信息与聚焦点平均信息距离值；
选择评分所用视野信息的距离值，并计算每个扫描视野的视野分值；
根据每个扫描视野的视野分值计算切片分值。
6.根据权利要求5所述的判定数字切片质量的装置，其特征是，所述图像增强模块包括：
视野图像处理模块，用于在样本区域外采集图像并取反得到增强图像,在采集视野图像后与增强图像进行加运算处理；
增强图像变量计算模块，用于计算增强图像的的色调、饱和度、颜色均值和清晰度。
7.根据权利要求5-6任意一项所述的判定数字切片质量的装置，其特征是，还包括：切片图像扫描模块，用于将切片分数保存到切片总分变量，如果切片分数大于设定值则保存扫描的切片图像，否则重新扫描。
CN202110520607.2A 2021-05-12 2021-05-12 判定数字切片质量的方法和装置 Active CN113256573B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110520607.2A CN113256573B (zh) 2021-05-12 2021-05-12 判定数字切片质量的方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110520607.2A CN113256573B (zh) 2021-05-12 2021-05-12 判定数字切片质量的方法和装置

Publications (2)

Publication Number Publication Date
CN113256573A CN113256573A (zh) 2021-08-13
CN113256573B true CN113256573B (zh) 2023-03-14

Family

ID=77181571

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110520607.2A Active CN113256573B (zh) 2021-05-12 2021-05-12 判定数字切片质量的方法和装置

Country Status (1)

Country Link
CN (1) CN113256573B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104316529A (zh) * 2014-11-18 2015-01-28 宁波江丰生物信息技术有限公司 一种全自动的用于数字病理图像扫描的实时画质评价方法
CN106204506A (zh) * 2016-06-30 2016-12-07 刘元生 印前数字图像质量评测方法
CN107945156A (zh) * 2017-11-14 2018-04-20 宁波江丰生物信息技术有限公司 一种自动评价数字病理扫描图像画质的方法
CN111337495A (zh) * 2020-04-09 2020-06-26 江苏迪赛特医疗科技有限公司 病理切片图像全程质控装置及病理切片分析装置
CN112697789A (zh) * 2020-12-09 2021-04-23 山东志盈医学科技有限公司 数字切片扫描仪的图像聚焦方法及装置

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8565503B2 (en) * 2010-06-04 2013-10-22 Leica Biosystems Imaging, Inc. System and method to determine slide quality of a digitized microscope slide
CN110337644A (zh) * 2017-02-23 2019-10-15 谷歌有限责任公司 用于辅助病理学家标识放大的组织图像中的肿瘤细胞的方法和系统
CN112200793B (zh) * 2020-10-21 2022-06-14 湖南国科智瞳科技有限公司 数字病理切片质量的实时监控方法及系统、计算机设备

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104316529A (zh) * 2014-11-18 2015-01-28 宁波江丰生物信息技术有限公司 一种全自动的用于数字病理图像扫描的实时画质评价方法
CN106204506A (zh) * 2016-06-30 2016-12-07 刘元生 印前数字图像质量评测方法
CN107945156A (zh) * 2017-11-14 2018-04-20 宁波江丰生物信息技术有限公司 一种自动评价数字病理扫描图像画质的方法
CN111337495A (zh) * 2020-04-09 2020-06-26 江苏迪赛特医疗科技有限公司 病理切片图像全程质控装置及病理切片分析装置
CN112697789A (zh) * 2020-12-09 2021-04-23 山东志盈医学科技有限公司 数字切片扫描仪的图像聚焦方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"Focus Quality Assessment of High-Throughput Whole Slide Imaging in Digital Pathology";Mahdi S. Hosseini et al.;《arXiv》;20181114;全文 *
"一种数字切片清晰度质量评分系统的应用";周燕燕等;《临床与实验病理学杂志》;20201231;全文 *

Also Published As

Publication number Publication date
CN113256573A (zh) 2021-08-13

Similar Documents

Publication Publication Date Title
US11914130B2 (en) Auto-focus methods and systems for multi-spectral imaging
US9224193B2 (en) Focus stacking image processing apparatus, imaging system, and image processing system
EP3239929B1 (en) Image processing apparatus, image processing method and program
US9824853B2 (en) Electron microscope device and imaging method using same
JP5826611B2 (ja) 硬さ試験機、及び硬さ試験方法
JP5134720B2 (ja) 仮想スライドの画質を評価するシステム及び方法
JP5374119B2 (ja) 距離情報取得装置、撮像装置、及びプログラム
US20120212599A1 (en) Imaging apparatus and imaging method
JP2014207110A (ja) 観察装置および観察方法
JP2006311505A (ja) 画像処理システム、画像処理方法、画像処理プログラム、撮像装置、及び撮像方法
CN115201092B (zh) 细胞扫描图像的获取方法和装置
JP4514731B2 (ja) Ｍｔｆ測定装置、ｍｔｆ測定方法およびｍｔｆ測定プログラム
CN110400281B (zh) 一种数字切片扫描仪中图像增强方法
CN113256573B (zh) 判定数字切片质量的方法和装置
CN111474103A (zh) 一种骨髓细胞载玻片自动聚焦扫描方法及系统
JP4847295B2 (ja) 共焦点顕微鏡
CN114785953B (zh) 基于sfr的相机自动对焦方法及装置
JP4197898B2 (ja) 顕微鏡、三次元画像生成方法、三次元画像を生成する制御をコンピュータに行わせるプログラム、及びそのプログラムを記録した記録媒体
JPWO2002045021A1 (ja) エントロピーフィルタ及び該フィルタを用いた領域抽出法
TWI510776B (zh) 玻璃氣泡瑕疵檢測處理方法
JPH06308390A (ja) 共焦点顕微鏡
JPH10170817A (ja) 合焦位置算出装置及びその方法とそれらを使用した電子顕微鏡
JP2005004287A (ja) 顔画像検出装置およびその制御方法
US20230258919A1 (en) A method for analysing scanning efficacy
JPH10154479A (ja) 合焦位置算出装置及びその方法

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