matlab函数mapminmax不适用与于图像处理中的灰度级的扩展与压缩

以下例子用于说明该问题：

在使用laplacian进行图像锐化时，我们会遇到数据类型的转化，同时伴随着产生了灰度级的扩展与压缩的问题。（使用matlab编程）

I_g = imread('3.tif');
imtool(I_g,[])
% 使用laplacian进行图像锐化，该过程会产生负值，因此需要使用double类型
w = fspecial('laplacian',0);
I_g_d = double(I_g);
g = imfilter(I_g_d,w,'replicate');
I_g_g = I_g_d-g;
% 在处理完之后，后续为了编程的方便，我们需要将数据恢复到0-255之间
I_g_g_g = mapminmax(I_g_g,0,255);
imtool(I_g_g_g,[]);

原图像：

 

输出的图像：

 

观察输出图像，发现水平方向出现一道一道的线，为什么会产生这样的现象，因为mapminmax意思是将矩阵的每一行处理成[-1,1]区间（参考：Matlab中的数据预处理-归一化(mapminmax)与标准化(mapstd)）。虽然该函数算法也使用y = (ymax-ymin)*(x-xmin)/(xmax-xmin) + ymin，这个算法可以实现灰度级的扩展与压缩，只是该函数使用的是一行，而不是全局的。

那么在使用laplacian进行图像锐化时，灰度级的扩展与压缩的问题怎么处理呢？（参考：）