From 28b50b1b09a80245fc4d957698ae817f9396f5a8 Mon Sep 17 00:00:00 2001
From: ac_fx <arc334@drexel.edu>
Date: Mon, 24 Jul 2023 16:40:54 -0400
Subject: [PATCH] draw_ssf_kymo for 3d
---
matlab/+SSF/draw_ssf_kymo.m | 22 ++++++++++++++++++++++
matlab/+SSF/quantize8.m | 6 +++++-
2 files changed, 27 insertions(+), 1 deletion(-)
create mode 100644 matlab/+SSF/draw_ssf_kymo.m
diff --git a/matlab/+SSF/draw_ssf_kymo.m b/matlab/+SSF/draw_ssf_kymo.m
new file mode 100644
index 00000000..865b8dcf
--- /dev/null
+++ b/matlab/+SSF/draw_ssf_kymo.m
@@ -0,0 +1,22 @@
+
+function imq = draw_ssf_kymo(strKymoFile)
+
+targetChannels = 1;
+imKymo = SSF.loadImage(strKymoFile,targetChannels);
+
+imp = imKymo;
+imp(imp<0)=0;
+imp = max(imp,[],1);
+imp = SSF.quantize8(imp);
+
+imn = imKymo;
+imn(imn>0)=0;
+imn=abs(imn);
+imn = max(imn,[],1);
+imn = SSF.quantize8(imn);
+
+imq(:,:,1) = squeeze(imp);
+imq(:,:,2) = squeeze(imn);
+imq(:,:,3) = 0*imp;
+
+
diff --git a/matlab/+SSF/quantize8.m b/matlab/+SSF/quantize8.m
index 1a2a0f11..1f553735 100644
--- a/matlab/+SSF/quantize8.m
+++ b/matlab/+SSF/quantize8.m
@@ -1,9 +1,13 @@
% clipLimits = [minVal,maxVal];
-function im = quantize8(im,clipLimits)
+function [im, clipLimits] = quantize8(im,clipLimits)
if isa(im,'uint8')
return
end
+
+if ~exist('clipLimits','var')
+ clipLimits = [min(im(:)), max(im(:))];
+end
im = max(im,clipLimits(1));
im = min(im,clipLimits(2));
--
GitLab