exp1-空间域图像增强

图像灰度变换

im = imread('pout.tif');
im = double(im) / 255.0;
myimshow(im, 'Original Image');

im1 = imadjust(im);
myimshow(im1, '灰度拉伸');

im2 = imadjust(im, [min(min(im)); max(max(im))], [], 0.4);
myimshow(im2, '\gamma 变换, \gamma = 0.4');

im3 = imadjust(im, [min(min(im)); max(max(im))], [], 2.5);
myimshow(im3, '\gamma 变换, \gamma = 2.5');

im1 = imadjust(im);
myimshow(im1, '灰度拉伸');

im2 = imadjust(im, [min(min(im)); max(max(im))], [], 0.4);
myimshow(im2, '\gamma 变换, \gamma = 0.4');

im3 = imadjust(im, [min(min(im)); max(max(im))], [], 2.5);
myimshow(im3, '\gamma 变换, \gamma = 2.5');

计算显示图像直方图、直方图均衡化

im = imread('tire.tif');
myimshow(im, '原图');

im_eq = histeq(im);
myimshow(im_eq, '直方图均衡化');

空间域图像平滑、锐化

im = imread('pout.tif');
im = double(im)/255;

figure;
subplot(221);
imshow(im);
title('Original Image')

index = 2;
for N = [3, 7, 15]
    kernel = ones(N) / N^2;
    im_blur = imfilter(im, kernel);
    subplot(2, 2, index);
    index = index + 1;
    imshow(im_blur);
    title(['Kernel size:', num2str(N)])
end
suptitle('Blur');

im = imread('rice.png');
im = double(im)/255;

figure, subplot(221);
imshow(im);title('Original Image');

kernel = [0 1 0; 1 -4 1; 0 1 0];
im_edge = imfilter(im, kernel);
subplot(222);
imshow(im_edge);
title('Uncalibrated');

im_edge = im_edge - min(min(im_edge));
subplot(223);
imshow(mat2gray(im_edge));
title('Calibrated');

im_sharpen = mat2gray(im - 0.5*im_edge);
subplot(224);
imshow(im_sharpen);
title('Sharpened');
suptitle('Sharpen');

function myimshow(im, im_title)
    figure;
    subplot(121);
    imshow(im);
    subplot(122);
    imhist(im);
    suptitle(im_title);
end