CN110727819A - 一种尺度自适应病理全切片图像数据库检索方法 - Google Patents

Abstract

本发明公开了一种尺度自适应病理全切片图像数据库检索方法，包括以下步骤：利用切片扫描仪将病理切片扫描到电子计算机中，得到数字病理图像，将所述数字病理图像在RGB空间中进行表示，生成数字病理全切片图像数据库；根据所述数字病理全切片图像数据库和预先训练好的基于卷积神经网络的病理图像块特征提取模型构建病理组织区域结构图；利用所述病理组织区域结构图训练图卷积神经网络，形成数据库索引；根据数据库索引进行尺度自适应在线检索。本发明检索框架能适应查询区域尺寸变化，可使单一检索模型满足不同尺寸、不同形状的查询图像的检索需求，且检索准确率高。

Description

一种尺度自适应病理全切片图像数据库检索方法

技术领域

本发明涉及医学影像处理与分析技术领域，更具体的说是涉及一种尺度自适应病理全切片图像数据库检索方法。

背景技术

数字病理全切片是通过专用的扫描成像系统将病理切片转换为可供计算机显示、传输及处理的高倍率大规模数字图像。随着数字病理图像成像技术不断成熟，数字病理图像自动分析方法发展迅速。基于内容的病理图像检索是其中的新兴研究方向，该方法利用计算机算法提取病理图像特征，以“图搜图”的形式在病理切片数据库中查找与医生感兴趣的图像最相似的病变区域，返回给医生查看，帮助医生快速收集与感兴趣图像相关的历史病例信息、快速积累病理学知识。

经数字化后的病理全切片图像(Whole slide image,WSI)像素分辨率远高于自然场景图像，通常能达到50000×50000像素以上的规模，而医生在发起检索时选取的感兴趣域(Region ofinterest,ROI)远小于全切片图像，且区域的尺寸与形状因肿瘤组织形态的差异存在较大变化，解决不规则查询ROI在病理全切片中的检索问题是将病理图像检索技术推向实际应用的关键。

现有病理图像检索方法大多仅适用于矩形ROI的检索情形，医生无法根据组织形态分布自主规定需要考虑的区域，且ROI尺寸固定，与实际应用仍有较大差距。

因此，如何提供一种尺度自适应病理全切片图像数据库检索方法是本领域技术人员亟需解决的问题。

发明内容

有鉴于此，本发明提供了一种尺度自适应病理全切片图像数据库检索方法，本发明检索框架能适应查询区域尺寸变化，可使单一检索模型满足不同尺寸、不同形状的查询图像的检索需求，且检索准确率高。

为了实现上述目的，本发明采用如下技术方案：

一种尺度自适应病理全切片图像数据库检索方法，包括以下步骤：

步骤1：将数字病理全切片图像在RGB空间中进行表示，生成数字病理全切片图像数据库；

步骤2：根据所述数字病理全切片图像数据库和预先训练好的基于卷积神经网络的病理图像块特征提取模型构建病理组织区域结构图；

步骤3：利用所述病理组织区域结构图训练图卷积神经网络，形成数据库索引；

步骤4：根据所述数据库索引进行尺度自适应在线检索。

进一步，所述步骤1具体为：

获取多张所述数字病理全切片图像中每个像素点在RGB三个通道的像素值，得到所述数字病理全切片图像数据库D，其计算公式为：

D＝{s_w|w＝1,2,...,W} (1)

其中，

表示一张像素分辨率为l_w×m_w的三通道数字病理全切片图像，W表示数据库中包含数字病理全切片图像的数量。

进一步，训练基于卷积神经网络的病理图像块特征提取模型具体包括以下步骤：

A、使用闭合曲线对所述数字病理全切片图像数据库D中所有病变区域按照病变所属类型进行明确标记，所述类型为0，1，...，C，其中C表示所述预先存储的数据库中包含的病变类型总数；

B、利用滑窗法对所述数字病理全切片图像数据库D中的区域进行采样，生成大小相同的病理图像块样本训练集，计算公式为：

X_train＝{(x_i,y_i)|i＝1,...,N_train} (2)

其中，x_i∈{0,1,...,255}^224×224×3表示病理图像块样本训练集中第i个图像块，N_train表示病理图像块样本训练集中图像块的数量，y_i∈{0,1,...,C}表示x_i的标签，其取值方法为:

C、利用所述病理图像块样本训练集训练卷积神经网络，训练完成后，将所述卷积神经网络末端的平均池化层的输出作为病理图像特征，表示为：

f_i＝DenseNet(x_i) (4)

其中，DenseNet(·)表示训练后的病理图像块特征提取模型。

进一步，所述卷积神经网络为DenseNet卷积神经网络。

进一步，所述步骤2具体包括以下步骤：

步骤21：将所述数字病理全切片图像数据库D中的数字病理全切片图像 s_w分割成大小为224×224像素的图像块，得到图像块集合X_w，表示为：

X_w＝{x_i|i＝1,...,N_w} (5)

其中，N_w表示数字病理全切片图像s_w包含的图像块数量， x_i∈{0,1,...,255}^224×224×3，表示图像块集合中第i个图像块；

步骤22：利用所述病理图像块特征提取模型提取所述图像块集合X_w中所有图像块特征F_w，其计算公式为：

步骤23：利用邻接矩阵A_w建立所述数字病理全切片图像包含的图像块之间的邻接关系，其中

表示所述数字病理全切片图像数据库D中第 w个数字病理全切片图像包含的图像块之间的邻接关系，当第i行第j列的元素a_ij＝1时表示第i个图像块x_i和第j个图像块x_j空间上相邻，a_ij＝0表示第i 个图像块x_i和第j个图像块x_j空间上不相邻；

步骤24：利用所述图像块特征F_w和所述邻接矩阵A_w计算所述数字病理全切片图像数据库D中所有相邻图像块之间的相似度，利用分级聚类算法依次合并最相似的图像块，生成连接区域；将所述连接区域中包含的图像块看作顶点，图像块之间的邻接关系看作边，将对应的所述图像块特征F_w看作顶点的属性，形成所述病理组织区域结构图。

进一步，所述步骤3具体包括以下步骤：

步骤31：利用所述病理组织区域结构图建立训练集G_train，其计算公式为：

G_train＝{(g_k,l_k)|k＝1,...,n} (7)

其中，G_train表示数字病理全切片图像数据库D中所有数字病理全切片图像包含的病理组织区域结构图的集合，其中n表示病理组织区域结构图的总量，g_k表示训练集中第k个病理组织区域结构图，l_k表示第k个病理组织区域结构图的标签，根据g_k中顶点对应的图像块标签最大数投标规则确定；

步骤32：利用所述训练集G_train训练所述图卷积神经网络，具体采用 DiffPool模型，利用DiffPool模型对所有所述病理组织区域结构图进行编码，得到病理组织区域结构图的量化表示r_i：

r_i＝DiffPool(g_i) (8)

步骤33：通过病理组织区域结构图的量化表示r_i，形成所述数据库索引R，表示为

R＝{r_i|i＝1,2,...,n} (9)

其中，n表示病理组织区域结构图的总量。

进一步，所述步骤4具体包括以下步骤：

步骤41：病理医生浏览需要诊断的数字病理图像时，勾画感兴趣区域；

步骤42：利用滑窗法将所述感兴趣区域划分图像块，得到图像块集合X_q，表示为：

X_q＝{x_i|i＝1,...,N_q} (10)

其中，N_q表示感兴趣区域中包含图像块的数量；

步骤43：利用所述病理图像块特征提取模型提取所述图像块集合X_q中所有图像块特征F_q，其计算公式为：

步骤44：利用邻接矩阵建立所述感兴趣区域中图像块的邻接关系；

步骤45：根据所述邻接矩阵A_q和所述图像块特征F_q建立查询组织图g_q；

步骤46：利用所述步骤32获得的DiffPool模型对所述查询组织图g_q进行编码，获得感兴趣区域的量化表示r_q，其计算公式为：

r_q＝DiffPool(g_q) (12)

步骤47：逐一计算所述感兴趣区域的量化表示r_q与所述病理组织区域结构图的量化表示r_i的相似性，将相似性从大到小排序，并按照相似性从大到小的顺序，将所述病理组织区域结构图对应的病理区域返回给病理医生，完成检索。

进一步，所述感兴趣区域的尺寸和形状没有限制。

进一步，所述步骤47中相似性计算为欧式距离、余弦距离、卡方距离或汉明距离。

经由上述的技术方案可知，与现有技术相比，本发明公开提供了一种尺度自适应病理全切片图像数据库检索方法，本发明检索框架能适应查询区域尺寸变化，可使单一检索模型满足不同尺寸、不同形状的查询图像的检索需求，且检索准确率高。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。

图1附图为本发明提供的一种尺度自适应病理全切片图像数据库检索方法流程图。

图2附图为本发明提供的病理全切片图像组织结构信息量化编码与检索技术路线图。

图3附图为本发明提供的病理全切片病变区域标记示意图。

图4附图为本发明提供的病理组织区域结构图生成示意图，其中(a)为一张全切片图像的缩略图，(b)为使用步骤24所提组织区域生成方法确定的连接区域，其中每个颜色块代表一个连接区域，(c)为在连接区域内建立的TG的二维空间分布图，(d)为其中一个TG与组织的对应关系，途中红色点代表TG的顶点，同时是图像块的中心，黑色线段表示TG的边。

图5附图为本发明提供的使用本发明方法进行检索的示意图，其中最左侧一列为查询区域，右侧为返回结果(从左到右按相似性从高到低排序)，返回结果中与查询图像属于同一类别(即检索正确)的区域用绿色框标记，属于不同类别(检索错误)的区域用红色标记。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本发明实施例公开了一种尺度自适应病理全切片图像数据库检索方法，本发明针对现有方法的不足，提出基于图卷积神经网络(Graph Convolution Network,GCN)的组织排布信息量化方法，以及相应的面向检索的WSI编码方法。方法流程如图2所示，在基于图像分块与CNN(卷积神经网络)的病理全切片图像编码方法基础上，综合考虑图像块的空间近邻关系与图像块特征相似关系，在切片内部建立组织区域结构图，并以组织区域结构图为单位定义WSI中的检索可返回对象；然后，将组织区域结构图顶点(Vertex)的 CNN特征以及边(Edge)的权重信息作为GCN的输入，将编码层作为GCN 网络的输出层，用于组织图结构的量化编码，从而建立检索数据库索引。

具体包括以下步骤，结合图1，

步骤1：用切片扫描仪将病理切片扫描到电子计算机中，将图像在RGB 颜色空间中进行表示；

数字病理切片是由专用的切片扫描仪将病理切片扫描并存储在计算机中的数字病理图像，并通过通道分解获取图像每个像素点在红(R)、绿(G)、蓝(B)三个通道的像素值，从而获得数字病理图像的数据信息，这些数据信息是病理图像辅助诊断方法的基础，数字WSI数据库表示为 D＝{s_w|w＝1,2,...,W}，其中表示一张像素分辨率为 l_w×m_w的3通道数字病理全切片图像，W为数据库中包含数字病理全切片图像的数量。

步骤2：训练基于卷积神经网络的病理图像块特征提取模型；

本发明涉及的组织排布信息量化技术基于WSI分块和特征提取，具体采用基于CNN的病理图像块图像特征提取模型。CNN病理图像块特征提取模型的建立包括以下三个步骤：

步骤21、由病理专家使用闭合曲线对预先存储的数字病理全切片图像数据库D中的切片中所有病变区域按照病变所属类别进行明确标记，类别量化为数字0，1，...，C，其中C表示数据库中包含的病变类型总数，标记结果示意图如图3所示；

步骤22、使用滑窗法对全切片中的区域进行采样，生成大小相同的病理图像块样本训练集，表示为X_train＝{(x_i,y_i)|i＝1,...,N_train}，其中 x_i∈{0,1,...,255}^224×224×3表示训练集中第i个图像块，N_train表示训练集中包含图像块的数量，y_i∈{0,1,...,C}表示x_i的标签，其取值确定方法如下：

步骤23、使用所建立的数据集训练DenseNet卷积神经网络结构，训练完成后，将卷积神经网络末端的平均池化层的输出用作病理图像特征，表示为：

f_i＝DenseNet(x_i)

其中DenseNet(·)表示训练后的病理图像块特征提取模型。

步骤3：构建病理组织区域结构图；

步骤31、数字病理全切片图像分块：使用滑窗法将数字病理全切片图像数据库D中的数字病理全切片图像s_w在20×物镜倍率下分割成大小为224× 224像素的图像块，滑窗步长为224像素，数字病理全切片图像s_w中所有图像块表示为集合X_w＝{x_i|i＝1,...,N_w}，其中N_w表示数字病理全切片图像s_w中包含图像块的数量；

步骤32、数字病理全切片图像分块特征提取：使用病理图像块特征提取模型DenseNet(·)提取图像块集合X_w中所有图像块的特征F_w，表示为：

步骤33、邻接矩阵A_w构建：使用矩阵

定义数字病理全切片图像数据库D第w个数字病理全切片图像包含的图像块之间的邻接关系，A_w中第i行第j列的元素a_ij取值定义如下：

步骤34、分级聚类生成组织区域结构图(Tissue Graph,TG)：以切片中图像块特征f_i之间的余弦距离为度量，以邻接矩阵A_w为依据，计算所有相邻图像块之间的相似度；然后使用分级聚类算法依次合并最相似的图像块，生成连接区域，合并过程中限定连接区域中包含图像块的最大数量；将连接区域中包含图像块看作TG的顶点，图像块的邻接关系看作TG的边，将对应的图像块特征看作顶点的属性，从而构造病理组织区域结构图，该步骤示意图见图4。

这里需要说明的是，具体实现中采用了欧式距离度量相邻图像块的相似度，但也可以使用余弦距离、卡方距离。

步骤4、训练图卷积网络(Graph convolutional network,GCN)，并对数字病理全切片图像数据库编码；

采用图卷积网络DiffPool结构对上述组织区域图进行编码，具体包含以下三个步骤：

步骤41、建立训练集：使用G_train＝{(g_k,l_k)|k＝1,...,n}表示数字病理全切片图像数据库D中所有数字病理全切片图像包含的TG的集合，其中n表示数据库中包含TG的总量，g_k表示训练集中第k个TG，l_k表示第k个TG的标签，根据g_k中顶点对应的图像块标签最大数投标规则确定；

步骤42、训练DiffPool网络结构：训练图嵌入与分类模型，训练完成后，将网络末端的平均池化层的输出用作TG量化表示，生成DiffPool模型，利用 DiffPool模型对所有病理组织区域结构图进行编码，得到病理组织区域结构图的量化表示r_i，表示为：

r_i＝DiffPool(g_i)

步骤43、建立数据库索引：通过病理组织区域结构图的量化表示r_i形成数据库索引，表示为R＝{r_i|i＝1,2,...,n}，其中，n表示病理组织区域结构图的总量，本发明涉及使用图卷积网络对TG进行编码包含二进制编码或哈希编码。

步骤5、尺度自适应在线检索应用；

数据库索引建立完成后，即可对病理医生提供的查询区域进行检索。在线检索分为以下步骤：

步骤51、病理医生浏览需要诊断的切片时勾画感兴趣区域，感兴趣区域的尺寸和形状没有限制；

步骤52、使用224×224的滑动窗口，以224像素为步长，将感兴趣区域在20×物镜倍率下划分为图像块，得到图像块集合X_q，表示为 X_q＝{x_i|i＝1,...,N_q}，其中N_q表示感兴趣区域中包含图像块的数量；

步骤53、依据步骤32的方法提取所有图像块特征F_q：

步骤54、依据步骤33的方法建立图像块的邻接矩阵A_q，

步骤55、根据邻接矩阵A_q和图像块特征F_q建立查询组织图g_q；具体构建过程为：g_q中包含N_q个节点；根据A_q确定g_q中节点的相邻关系，从而定义组织图的边；将图像块特征f_i∈F_q定义为g^q中第i个节点的特征表示。

步骤56、使用步骤42中获得的DiffPool网络结构对查询组织图g_q进行编码，获得感兴趣区域的量化表示r_q:r_q＝DiffPool(g_q)

步骤57、逐一计算感兴趣区域的量化表示r_q与数据库索引R中TG量化表示r_i之间的相似性，将相似性从大到小排序，并按照相似性从大到小的顺序，将数据库索引R中TG对应的病理图像区域返回给医生，完成检索。此处的相似性度量可使用但不限于欧式距离、余弦距离、卡方距离，以及针对二值编码的汉明距离等。

通过图5可知，可见本发明方法能够准确的在数据库检索返回与查询图像相关的病理图像区域，且对查询区域的尺寸和形状有较强的适应性，检索准确率高。

为了验证本发明的效果，在公开数据集上与现有方法进行了对比，结果如表1所示。可见本发明在检索精度上均高于现有方法，检索速度上与现有方法中最优效果相当。

表1本发明病理图像检索结果与现有检索结果指标对比

表中依次给出了检索返回前50、200个结果的精度(P@50/200)和平均精度(MAP@50/200)，以及单次检索的平均时间。

本发明在基于图像分块与CNN的病理全切片图像特征提取的基础上，综合利用了图像块的空间近邻关系与图像块特征相似关系，在切片内部建立组织区域结构图，并以组织区域结构图为单位定义WSI中的检索可返回对象；然后，将组织区域结构图顶点(Vertex)的CNN特征以及边(Edge)的权重信息作为GCN的输入，获得图结构的嵌入编码，从而建立病理图像检索数据库索引。

该检索框架相比于传统方法具有以下优势：1)病理图像组织区域的排布信息被嵌入到检索索引中，使检索模型能把握病理图像结构上的相似性，使返回结果与查询区域更为相似；

2)研究采用GCN模型对输入图的顶点数量有较强的适应能力，使检索框架能适应查询区域尺寸变化，可使单一检索模型满足不同尺寸、不同形状的查询图像的检索需求；

3)解决了现有方法通常使用全局池化的方式提取ROI特征表示，难以对组织病理学中需要考虑组织区域近邻关系进行描述，检索精度有限。

本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。

对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (9)

1.一种尺度自适应病理全切片图像数据库检索方法，其特征在于，包括以下步骤：

步骤1：将数字病理全切片图像在RGB空间中进行表示，生成数字病理全切片图像数据库；

步骤2：根据所述数字病理全切片图像数据库和预先训练好的基于卷积神经网络的病理图像块特征提取模型构建病理组织区域结构图；

步骤3：利用所述病理组织区域结构图训练图卷积神经网络，形成数据库索引；

步骤4：根据所述数据库索引进行尺度自适应在线检索。

2.根据权利要求1所述的一种尺度自适应病理全切片图像数据库检索方法，其特征在于，所述步骤1具体为：

获取多张所述数字病理全切片图像中每个像素点在RGB三个通道的像素值，得到所述数字病理全切片图像数据库D，其计算公式为：

D＝{s_w|w＝1,2,...,W} (1)

其中，

表示一张像素分辨率为l_w×m_w的三通道数字病理全切片图像，W表示数据库中包含数字病理全切片图像的数量。

3.根据权利要求2所述的一种尺度自适应病理全切片图像数据库检索方法，其特征在于，训练基于卷积神经网络的病理图像块特征提取模型具体包括以下步骤：

A、使用闭合曲线对所述数字病理全切片图像数据库D中所有病变区域按照病变所属类型进行明确标记，所述类型为0，1，...，C，其中C表示数据库D中包含的病变类型总数；

B、利用滑窗法对所述数字病理全切片图像数据库D中的区域进行采样，生成大小相同的病理图像块样本训练集，计算公式为：

X_train＝{(x_i,y_i)|i＝1,...,N_train} (2)

其中，x_i∈{0,1,...,255}^224×224×3表示病理图像块样本训练集中第i个图像块，N_train表示病理图像块样本训练集中图像块的数量，y_i∈{0,1,...,C}表示x_i的标签，其取值方法为:

C、利用所述病理图像块样本训练集训练卷积神经网络，训练完成后，将所述卷积神经网络末端的平均池化层的输出作为病理图像特征，表示为：

f_i＝DenseNet(x_i) (4)

其中，DenseNet(·)表示训练后的病理图像块特征提取模型。

4.根据权利要求3所述的一种尺度自适应病理全切片图像数据库检索方法，其特征在于，所述卷积神经网络为DenseNet卷积神经网络。

5.根据权利要求3所述的一种尺度自适应病理全切片图像数据库检索方法，其特征在于，所述步骤2具体包括以下步骤：

步骤21：将所述数字病理全切片图像数据库D中的数字病理全切片图像s_w分割成大小为224×224像素的图像块，得到图像块集合X_w，表示为：

X_w＝{x_i|i＝1,...,N_w} (5)

其中，N_w表示数字病理全切片图像s_w包含的图像块数量，x_i∈{0,1,...,255}^224×224×3，表示图像块集合中第i个图像块；

步骤22：利用所述病理图像块特征提取模型提取所述图像块集合X_w中所有图像块特征F_w，其计算公式为：

步骤23：利用邻接矩阵A_w建立所述数字病理全切片图像包含的图像块之间的邻接关系，其中表示所述数字病理全切片图像数据库D中第w个数字病理全切片图像包含的图像块之间的邻接关系，当第i行第j列的元素a_ij＝1时表示第i个图像块x_i和第j个图像块x_j空间上相邻，a_ij＝0表示第i个图像块x_i和第j个图像块x_j空间上不相邻；

步骤24：利用所述图像块特征F_w和所述邻接矩阵A_w计算所述数字病理全切片图像数据库D中所有相邻图像块之间的相似度，利用分级聚类算法依次合并最相似的图像块，生成连接区域；将所述连接区域中包含的图像块看作顶点，图像块之间的邻接关系看作边，将对应的所述图像块特征F_w看作顶点的属性，形成所述病理组织区域结构图。

6.根据权利要求5所述的一种尺度自适应病理全切片图像数据库检索方法，其特征在于，所述步骤3具体包括以下步骤：

步骤31：利用所述病理组织区域结构图建立训练集G_train，其公式为：

G_train＝{(g_k,l_k)|k＝1,...,n} (7)

其中，G_train表示数字病理全切片图像数据库D中所有数字病理全切片图像包含的病理组织区域结构图的集合，其中n表示病理组织区域结构图的总量，g_k表示训练集中第k个病理组织区域结构图，l_k表示第k个病理组织区域结构图的标签，根据g_k中顶点对应的图像块标签最大数投标规则确定；

步骤32：利用所述训练集G_train训练所述图卷积神经网络，具体采用DiffPool模型，利用DiffPool模型对所有所述病理组织区域结构图进行编码，得到病理组织区域结构图的量化表示r_i：

r_i＝DiffPool(g_i) (8)

步骤33：通过病理组织区域结构图的量化表示r_i形成所述数据库索引R，表示为

R＝{r_i|i＝1,2,...,n} (9)

其中，n表示病理组织区域结构图的总量。

7.根据权利要求6所述的一种尺度自适应病理全切片图像数据库检索方法，其特征在于，所述步骤4具体包括以下步骤：

步骤41：病理医生浏览需要诊断的数字病理图像时，勾画感兴趣区域；

步骤42：利用滑窗法将所述感兴趣区域划分图像块，得到图像块集合X_q，表示为：

X_q＝{x_i|i＝1,...,N_q} (10)

其中，N_q表示感兴趣区域中包含图像块的数量；

步骤43：利用所述病理图像块特征提取模型提取所述图像块集合X_q中所有图像块特征F_q，其计算公式为：

步骤44：利用邻接矩阵建立所述感兴趣区域中图像块的邻接关系；

步骤45：根据所述邻接矩阵A_q和所述图像块特征F_q建立查询组织图g_q；

步骤46：利用所述步骤32获得的DiffPool模型对所述查询组织图g_q进行编码，获得感兴趣区域的量化表示r_q，其计算公式为：

r_q＝DiffPool(g_q) (12)

步骤47：逐一计算所述感兴趣区域的量化表示r_q与所述病理组织区域结构图的量化表示r_i的相似性，将相似性从大到小排序，并按照相似性从大到小的顺序，将所述病理组织区域结构图对应的病理区域返回给病理医生，完成检索。

8.根据权利要求7所述的一种尺度自适应病理全切片图像数据库检索方法，其特征在于，所述感兴趣区域的尺寸和形状没有限制。

9.根据权利要求7所述的一种尺度自适应病理全切片图像数据库检索方法，其特征在于，所述步骤47中相似性计算为欧式距离、余弦距离、卡方距离或汉明距离。

Images