diff --git a/matlab/+SSF/draw_ssf_kymo.m b/matlab/+SSF/draw_ssf_kymo.m
new file mode 100644
index 0000000000000000000000000000000000000000..865b8dcf8704fb05b5564a6b414d1e459a7aaf8f
--- /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 1a2a0f119da3e86867c2433fc4da199bb0480fd6..1f553735ee9dacebe9047cbcc26945c7ddbcc563 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));