作者,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()
基础处理