MERSCOPE

[2]:
import warnings
warnings.filterwarnings("ignore")
[3]:
import SOView
import scanpy as sc
[4]:
adata_raw =sc.read_h5ad('/data/linsenlin/binary/binary1/MERSCOPE/R1S2_knn15.h5ad')
[7]:
adata = adata_raw.copy()

[8]:
sc.pp.neighbors(adata,use_rep='BINARY')
[10]:
sc.tl.umap(adata,n_components=3)
[11]:
adata.obsm['X_umap'].shape
[11]:
(83546, 3)
[12]:
SOView.SOViewer_plot(
    adata = adata, # the data to plot
    save = f'R1S2', # save the result to specified path or don't save (None)
    embedding_use='X_umap', # which embedding to be used for plot
    dot_size=1, # the marker size of the plot
)

generating color coding...
1.0 0.0
1.0 0.0
../_images/Tissue_continuity_visualization_MERSCOPE_8_1.png