forked from daviddaiweizhang/istar
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathrun copy.sh
67 lines (52 loc) · 3 KB
/
run copy.sh
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
#!/bin/bash
set -e
prefix=$1 # e.g. data/demo/
device="cuda" # "cuda" or "cpu"
pixel_size=0.5 # desired pixel size for the whole analysis
n_genes=1000 # number of most variable genes to impute
# preprocess histology image
echo $pixel_size > ${prefix}pixel-size.txt
python rescale.py ${prefix} --image #he-raw.jpg -> he-scaled.jpg
python preprocess.py ${prefix} --image #he-scaled.jpg -> he.jpg
# extract histology features
python extract_features.py ${prefix} --device=${device} #output->embeddings-hist-raw.pickle
# # If you want to retun model, you need to delete the existing results:
# rm ${prefix}embeddings-hist-raw.pickle
# auto detect tissue mask
# If you have a user-defined tissue mask, put it at `${prefix}mask-raw.png` and comment out the line below
python get_mask.py ${prefix}embeddings-hist.pickle ${prefix}mask-small.png
# # segment image by histology features
# python cluster.py --mask=${prefix}mask-small.png --n-clusters=10 ${prefix}embeddings-hist.pickle ${prefix}clusters-hist/
# # # segment image by histology features without tissue mask
# # python cluster.py ${prefix}embeddings-hist.pickle ${prefix}clusters-hist/unmasked/
# select most highly variable genes to predict
# If you have a user-defined list of genes, put it at `${prefix}gene-names.txt` and comment out the line below
python select_genes.py --n-top=${n_genes} "${prefix}cnts.tsv" "${prefix}gene-names.txt" #output : Top n genes -> gene-names.txt
# predict super-resolution gene expression
# rescale coordinates and spot radius
python rescale.py ${prefix} --locs --radius
# train gene expression prediction model and predict at super-resolution
python impute.py ${prefix} --epochs=400 --device=${device} # train model from scratch
# # If you want to retrain model, you need to delete the existing model:
# rm -r ${prefix}states
# visualize imputed gene expression
python plot_imputed.py ${prefix}
# segment image by gene features
python cluster.py --filter-size=8 --min-cluster-size=20 --n-clusters=10 --mask=${prefix}mask-small.png ${prefix}embeddings-gene.pickle ${prefix}clusters-gene/
# # segment image without tissue mask
# python cluster.py --filter-size=8 --min-cluster-size=20 ${prefix}embeddings-gene.pickle ${prefix}clusters-gene/unmasked/
# # segment image without spatial smoothing
# python cluster.py --mask=${prefix}mask-small.png ${prefix}embeddings-gene.pickle ${prefix}clusters-gene/unsmoothed/
# python cluster.py ${prefix}embeddings-gene.pickle ${prefix}clusters-gene/unsmoothed/unmasked/
# differential analysis by clusters
python aggregate_imputed.py ${prefix}
python reorganize_imputed.py ${prefix}
python differential.py ${prefix}
# visualize spot-level gene expression data
python plot_spots.py ${prefix}
# cell type inference
python pixannot.py ${prefix} data/markers/celltype.tsv ${prefix}markers/celltype/
cp -r ${prefix}markers/celltype/threshold010/* ${prefix}markers/celltype/
python enrich.py ${prefix}clusters-gene/ ${prefix}markers/celltype/
# TLS detection
python marker_score.py ${prefix} data/markers/tls.txt ${prefix}markers/tls