diff --git a/matlab/+SSF/getClipQuantiles.m b/matlab/+SSF/getClipQuantiles.m
index d40a06d45f6aa05308be187762f3decd22e2e24e..7ba6038cee7092d1efbb3f3372f444837e26ac2a 100644
--- a/matlab/+SSF/getClipQuantiles.m
+++ b/matlab/+SSF/getClipQuantiles.m
@@ -11,7 +11,7 @@ cl = [mx-sx, mx+sx];
 % end
 
 % linear quantization
-clipLimits = linspace(cl(1),cl(2),nQuant);
+clipLimits = cl; %linspace(cl(1),cl(2),nQuant);
 
 % % equiprobable quantization
 % ep1 = length(find(voxels<=cl(1)))./length(voxels);
diff --git a/matlab/+SSF/quantize8.m b/matlab/+SSF/quantize8.m
index b74a16a41960bf5bcbb82d695bfce2dfd2dae24a..d5dd5a70c4d7ce79284b7efb39092c94b307cf33 100644
--- a/matlab/+SSF/quantize8.m
+++ b/matlab/+SSF/quantize8.m
@@ -6,6 +6,13 @@ if isa(im,'uint8')
     return
 end
 
-im = uint8(imquantize(im,clipLimits));
+% im = uint8(imquantize(im,clipLimits));
+
+im = max(im,clipLimits(1));
+im = min(im,clipLimits(end));
+im = (im - clipLimits(1)) ./ (clipLimits(end) - clipLimits(end));
+im = im2uint8(im);
 im = max(im,1);
-im(0 == im_in) = 0;
\ No newline at end of file
+im(0 == im_in) = 0;
+4;
+