import anndata
import pandas as pd
adata=anndata.read("/home/csq/./GSM6058681_TabulaSapiens.h5ad")
adata(查看文件中的结构)
adata.obs(查看样本信息)
adata.var (查看特征信息)
提取样本和特征:
先获取cell或者gene list
gene1812 = adata.var_names.to_list()[:1812]
cell100000 = adata.obs_names.to_list()[:100000]
data2 = adata[cell100000,gene1812] # 100000 x gene
print(data2.shape,type(data2))
详情见:https://www.jianshu.com/p/9c1a5133b2da