本发明涉及生物信息数据的分析方法，是一种计算机对空间转录组数据特征选择及聚类方法，属于生物信息学领域。
背景技术：
1、通过研究特征提取与聚类领域的算法，对空间转录组进行空间域划分即聚类。由于空间转录组的基因数量巨大，我们最先遇到的问题便是数据点缺失，并不是每一个细胞或者细胞群的每一个基因都可以被检测到，这个时候数据就出现错误了，而最好的解决这个问题的方法就是做特征提取，通常情况下我们使用的是pca或者svg作为空间转录组数据的特征提取方式，但这些方法可能无法捕捉到数据的所有重要特征；并且另一个问题就是，空间转录组数据的庞大必然带来的一个不好的问题就是数据的稀疏，很多细胞或者细胞群的基因并不表达，这个时候我们就需要利用图去聚合周围邻居的节点信息。
2、因此，本发明隶属于科技部国家重点研发课题(2021yff1201203)，重点研究生物医学知识图谱自动构建和更新技术体系,主要对基于图卷积变分自编码器对空间转录组数据进行特征提取和聚类，克服了以上的问题并且能够更好的对空间域进行划分。
技术实现思路
1、为了克服现有技术的不足，本发明提供了一种计算机对空间转录组数据特征选择及聚类方法，为达此目的，本发明采用以下技术方案：
2、步骤一、在对空间转录组数据进行聚类之前，首先需要对空间转录组数据进行预处理，预处理包括有子步骤：子步骤1：归一化，对空间转录组数据归一化，通过将每个基因的表达值减去其平均值并除以其标准差来实现和得到空间转录组数据1，以用于确保所有基因的表达级别在相同的尺度上，从而消除不同基因之间的尺度差异，使得聚类分析准确；子步骤2：对子步骤1处理后的归一化后的空间转录组数据1去除噪声和不相关数据点；在子步骤2中设置阈值得到空间转录组数据2，以只保留表达值高于某个阈值的基因，能提高聚类的准确性，子步骤3，对子步骤1和子步骤2处理后的空间转录组数据2降维：空间转录组数据2通常具有高维性，不能直接对空间转录组数据2进行聚类分析，子步骤3使用高可变基因(hvg)降维技术来将数据投影到低维空间，然而，当子步骤1到子步骤3执行结束后，存在无法捕捉到数据的所有重要特征情况时，空间转录组数据2作为高维输入数据执行步骤二；步骤二：使用变分自编码器vae来学习数据的低维表示,变分自编码器vae用于学习高维输入数据的低维表示；子步骤4：变分自编码器vae执行部件为编码器和解码器，编码器将高维输入数据压缩成低维隐向量，解码器将低维隐向量重新构造为原始维度,执行变分自编码器vae关键特性，变分自编码器vae关键特性是能够生成新的、与训练数据相似的数据点；步骤三：为使用变分自编码器vae学习高维输入数据的低维表示，需要执行以下步骤：子步骤5：构建变分自编码器vae的空间转录组最优化变分自编码器聚类学习模型：首先，需要构建一个空间转录组最优化变分自编码器聚类学习模型,变分自编码器vae的空间转录组最优化变分自编码器聚类学习模型包括定义编码器和解码器的神经网络结构，以及选择合适的激活函数和优化器；编码器和解码器用卷积神经网络或循环神经网络来实现，并基于空间转录组的特点，利用空间位置信息构成图，因此选用图卷积神经网络作为编码器；子步骤6：训练变分自编码器vae的空间转录组最优化变分自编码器聚类学习模型：接下来使用高维输入数据来训练变分自编码器vae的空间转录组最优化变分自编码器聚类学习模型,子步骤6所有过程涉及到最小化三个损失函数：均方损失函数，用于衡量解码器输出与原始数据的差异、kl散度，用于衡量编码器输出的隐向量与预设的先验分布之间的差异以及最优化传输损失，用于衡量编码器输出的隐向量与解码器输出的差异；通过子步骤6的执行，变分自编码器vae的空间转录组最优化变分自编码器聚类学习模型能学习到一个能够捕捉到数据主要变化趋势的低维空间，并生成子步骤7的低维表示：一旦vae训练完成，子步骤7使用编码器部分将高维输入数据转换为低维隐向量1,低维隐向量1捕捉数据的主要变化趋势，因此在低维隐向量1构成的宏空间上进行聚类得到精确结果；步骤四：在获得高维输入数据的低维表示得到低维表示高维输入数据data后，使用k-means聚类算法对低维表示高维输入数据data进行聚类。k-means聚类无监督学习方法，用于将数据点划分为k个群集，通过最小化每个数据点到其所属群集中心的距离来工作，为使用k-means对低维表示高维输入数据data进行聚类，执行以下子步骤：子步骤8：选择k值：首先，确定使用肘部法则或轮廓系数方法来实现将低维表示高维输入数据data划分为n个群集,n为自然数，并找到一个合适的k值,k为自然数，使得n个群集之间的差异最大化，而群集内部的差异最小化；子步骤9：可通过随机选择k个数据点，并使用k-means++算法来实现初始化k个群集中心，k-means++算法能有效地选择初始群集中心，从而加速算法的收敛性；子步骤10：迭代更新群集中心,在每次迭代更新群集中心中，将每个数据点分配给距离其最近的群集中心，并更新群集中心为其群集中心所属数据点的均值，子步骤10这个过程将不断重复，直到群集中心不再发生显著变化；计算聚类结果：最后子步骤11：计算聚类结果，即每个数据点所属的群集；计算聚类结果可通过将每个数据点分配给距离其最近的群集中心来实现；步骤五：在获得聚类结果后，需要解析和验证这些结果，通过以下子步骤来实现：子步骤12：查看基因表达模式：通过查看每个群集中基因的表达模式，以了解群集在生物学上的含义，即a:检查每个群集中基因的功能注释，以确定是否具有相似的生物学过程或通路；b:使用生物学标记验证结果：c:使用已知的生物学标记来验证聚类结果，这可通过比较结果与已知的细胞类型或组织结构信息来实现，如果得到的聚类结果与已知的生物学信息一致，那么则判断结果是可靠的；d:生成新的数据点：e:使用变分自编码器vae生成新的、具有相似表达模式的数据点,然后以查看这些新的、具有相似表达模式的数据点被分配到哪个群集，以进一步验证变分自编码器vae的空间转录组最优化变分自编码器聚类学习模型。
1.一种计算机对空间转录组数据特征选择及聚类方法，其特征在于：包括有：步骤一、在对空间转录组数据进行聚类之前，首先需要对所述空间转录组数据进行预处理，所述预处理包括有子步骤：子步骤1：归一化，对所述空间转录组数据归一化，通过将每个基因的表达值减去其平均值并除以其标准差来实现和得到空间转录组数据1，以用于确保所有基因的表达级别在相同的尺度上，从而消除不同基因之间的尺度差异，使得聚类分析准确；子步骤2：对子步骤1处理后的归一化后的空间转录组数据1去除噪声和不相关数据点；在子步骤2中设置阈值得到空间转录组数据2，以只保留表达值高于某个阈值的基因，能提高聚类的准确性，子步骤3，对子步骤1和子步骤2处理后的空间转录组数据2降维：空间转录组数据2通常具有高维性，不能直接对所述空间转录组数据2进行聚类分析，子步骤3使用高可变基因降维技术来将数据投影到低维空间，然而，当子步骤1到子步骤3执行结束后，存在无法捕捉到数据的所有重要特征情况时，所述空间转录组数据2作为高维输入数据执行步骤二；
