在python中读取h5ad文件

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

  • 0
    点赞
  • 5
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值