From 82fbfdbd0adab46bb001212c6f318b6482630776 Mon Sep 17 00:00:00 2001
From: la_29 <la538@drexel.edu>
Date: Fri, 10 Jan 2025 11:33:21 -0500
Subject: [PATCH] code to generate figure 3
---
src/MATLAB/ssf_vs_cn/goFigure3-C.m | 0
src/MATLAB/ssf_vs_cn/goPhantomImage.m | 0
src/MATLAB/ssf_vs_cn/goSSF_CN.m | 48 +++++++++++++++++++++++++++
3 files changed, 48 insertions(+)
create mode 100644 src/MATLAB/ssf_vs_cn/goFigure3-C.m
create mode 100644 src/MATLAB/ssf_vs_cn/goPhantomImage.m
create mode 100644 src/MATLAB/ssf_vs_cn/goSSF_CN.m
diff --git a/src/MATLAB/ssf_vs_cn/goFigure3-C.m b/src/MATLAB/ssf_vs_cn/goFigure3-C.m
new file mode 100644
index 0000000..e69de29
diff --git a/src/MATLAB/ssf_vs_cn/goPhantomImage.m b/src/MATLAB/ssf_vs_cn/goPhantomImage.m
new file mode 100644
index 0000000..e69de29
diff --git a/src/MATLAB/ssf_vs_cn/goSSF_CN.m b/src/MATLAB/ssf_vs_cn/goSSF_CN.m
new file mode 100644
index 0000000..90c47e3
--- /dev/null
+++ b/src/MATLAB/ssf_vs_cn/goSSF_CN.m
@@ -0,0 +1,48 @@
+
+sz = [100,100];
+cx = sz./2;
+imCN = zeros(sz);
+imSSF = zeros(sz);
+
+for i = 1 : sz(1)
+ for j = 1 : sz(2)
+
+ cyto = i./sz(1);
+ nuc = j./sz(2);
+
+ imCN(i,j) = cyto./nuc;
+ imSSF(i,j) = cyto - nuc;
+
+ end
+end
+
+clf;imagesc(imCN)
+xt = get(gca,'xtick')
+set(gca,'xticklabels',xt./sz(2))
+yt = get(gca,'ytick')
+set(gca,'yticklabels',yt./sz(1))
+ylabel('cytoplasmic intensity')
+xlabel('nuclear intensity')
+
+% tick_targets = [1,64,128,192,255];
+cb = colorbar()
+cb.Label.String = 'cytonuclear ratio'
+
+exportgraphics(gcf,'imCN.png')
+
+clf;imagesc(imSSF)
+xt = get(gca,'xtick')
+set(gca,'xticklabels',xt./sz(2))
+yt = get(gca,'ytick')
+set(gca,'yticklabels',yt./sz(1))
+ylabel('cytoplasmic intensity')
+xlabel('nuclear intensity')
+
+tick_targets = [-0.99,-0.5,0,0.5,0.99];
+tick_labels = [-1,-0.5,0,0.5,1];
+cb = colorbar('ticks',tick_targets,'ticklabels',tick_labels)
+cb.Label.String = 'SSF'
+
+
+exportgraphics(gcf,'imSSF.png')
+
--
GitLab