课前准备---HD数据结合图像识别获取真实的空间单细胞级数据

作者,Evil Genius

HD数据不同于Xenium,目前还是横屏竖直的一刀切数据分析模式,但是真实的细胞绝对不是如此分布的,那么实际分析中,2um的精度配合图像的信息,获取真实的细胞分布数据,就成了分析的必须。

多说一句,分析的准确性和超前化也是公司对核心分析人员的核心要求。

如下如,我们最好不要采用8um,16um这种一刀切的模式,而是识别每个2um中是否含有细胞,如果有则保留,没有则舍弃。

这个时候拿到的数据才是有效数据,所做的分析才可以认为是非常可靠的。利用2um的精度重构细胞分布信息。

在Visium HD的最高2µm分辨率下,将亚细胞bin连接成单个细胞。这可以通过使用StarDist进行形态分割来完成的,使用其预训练的H&E模型识别细胞核,随后将其扩展到邻近的未标记的bin中。
我们需要实现如下的分析目标
分析得到的结果流程图
分析得到的结果
局部放大

我们来实现,是官方示例数据

import scanpy as sc
import os
import bin2cell as b2c
import celltypist
from celltypist import models
import numpy as np
from matplotlib import rcParams
from matplotlib import font_manager
import matplotlib.pyplot as plt
rcParams['pdf.fonttype'] = 42
sc.settings.set_figure_params(dpi = 150, color_map = 'RdPu', dpi_save = 150, vector_friendly = True, format = 'pdf')
font_manager.fontManager.addfont(".../software/Arial.ttf")
print(font_manager.findfont("Arial"))
plt.rcParams["font.sans-serif"] = ["Arial"]
sc.settings.set_figure_params(dpi = 150, color_map = 'RdPu', dpi_save = 300, vector_friendly = True, format = 'pdf')

path008 = ".../visium_hd/gut_public/square_008um/"
path002 = ".../visium_hd/gut_public/square_002um/"
source_image_path = ".../10X_datasets/human_CRC/Visium_HD_Human_Colon_Cancer_tissue_image.btf"

bdata = b2c.read_visium(path008, source_image_path = source_image_path)
bdata.var_names_make_unique()

bdata.raw = bdata.copy()
基础处理
还有 86% 的精彩内容
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。
禁止转载,如需转载请通过简信或评论联系作者。
支付 ¥200.00 继续阅读

相关阅读更多精彩内容

友情链接更多精彩内容