在python中读取h5ad文件

import anndata
import pandas as pd
adata=anndata.read("/home/csq/./GSM6058681_TabulaSapiens.h5ad")

adata(查看文件中的结构)
在python中读取h5ad文件_第1张图片

adata.obs(查看样本信息)
在python中读取h5ad文件_第2张图片

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

你可能感兴趣的:(python)