import scanpy as sc
dataset ="bct"
adata = sc.read("/home/yxk/Desktop/test_dataset/bct/bct_raw.h5ad")
print(adata)
print(adata.X[0:10,0:10])
adata.raw = adata
## raw_visulize
############################################
sc.pp.normalize_per_cell(adata,counts_per_cell_after=10000)
sc.pp.log1p(adata)
sc.pp.highly_vari
able_genes(adata, n_top_genes=1000, subset = True)
sc.tl.pca(adata,svd_solver='arpack')
print(adata)
sc.pp.neighbors(adata)
sc.tl.umap(adata)
sc.pl.umap(adata,color=["celltype","BATCH"])
############################################
# 重新回复原始的counts
aa=adata.raw.to_adata()
print(aa)
print(aa.X[0:10,0:10])