Matlab代码如下:
clear all;
clc;
PS = imread('moon.jpg');
imshow(PS)
title('原始图像')
imwrite(rgb2gray(PS), 'PicSampleGray.bmp');
f = im2double(PS);
PS = rgb2gray(PS);
figure, imshow(PS)
title('灰度化后的图像')
f = im2double(PS); %%由于像素变换会存在负值,转化为归一化的double数据
w = fspecial('laplacian', 0); %%标准Laplacian模板
g = imfilter(f, w, 'replicate'); %%边界复制后滤波并去除边界
figure, imshow(g)
title('滤波图像')
final = f - g;
figure, imshow(final)
title('锐化图像')