【实践】数字图像处理DIP课程课业打卡实验4 图像锐化
叮嘟！这里是小啊呜的学习课程资料整理。好记性不如烂笔头，今天也是努力进步的一天。一起加油进阶吧！
一、实验目的
1、掌握一阶锐化算法的工作原理和算法实现。
2、掌握二阶锐化算法的工作原理和算法实现。
二、实验内容
1、 实现拉普拉斯锐化。
测试图像：’Grayblur6.png’。
（备注：将执行正确的代码粘贴在此处，核心代码要求有注释）
im= imread('im\Grayblur6.png' );
imshow(im);
title('原始图像' );
im0=im;
im = double(im) ;
H=[0 -1 0;-1 4 -1;0 -1 0];
im1=filter2(H, im) ;
im2=uint8(im1+im) ;
figure;
imshow(im2) ;
title('Laplacian锐化叠加后图像');
im_L0 = uint8(im1);
figure, imshow(im_L0);
title('Laplacian鋭化');
maxV = max(max(im1));
minV = min(min(im1));
delta =0;
if (minV<0)
delta=abs(minV) ;
end
im12=(im1+delta);
figure,imshow(uint8(im12));
title('Laplacian锐化图像');
2、实现Sobel微分算子的代码。
测试图像：’Grayblur6.png’。
（备注：将执行正确的代码粘贴在此处，核心代码要求有注释）
im=imread('im\Grayblur6.png' );
imshow(im); title('原始图像' ) ;
im = double(im) ;
H=[-1 0 1;-2 0 2;-1 0 1];
Dy=filter2( H, im );
H=[-1 -2 -1;0 0 0;1 2 1];
Dx=filter2( H, im );
im_2 = sqrt(Dx.^2+Dy.^2);
im_2=uint8(im_2) ;
figure, imshow(im_2), title(' Sobel锐化图像' ); %显示滤波后的图像
im_3 = uint8 (im+sqrt (Dx.^2+Dy.^2)*0.1);
figure, imshow(im_3), title('Sobel锐化叠加图像' ); %显示滤波后的图像
3、实现交叉微分算法（Roberts算法）。
测试图像：’Grayblur6.png’。
读入图像到矩阵F
得到F的行列数与颜色通道数：m,n,c
如果颜色通道数c>1，则将彩色图像转化为灰度图rgb2gray
初始化输出矩阵G
对于F的每一行i（1到m-1）
对于F的每一列j（1到n-1）
根据以下公式，由F(i,j)、F(i+1,j)、F(i,j+1)、F(i+1,j+1)计算G(i,j)
显示Roberts锐化图像
显示Roberts锐化并叠加原图的结果图像
（备注：将执行正确的代码粘贴在此处，核心代码要求有注释）
im = imread('im\Grayblur6.png');
imshow(im);
title('原始图像');
im = double(im);
[m,n,c]=size(im);
g=zeros(m,n);
for i=1:m-1
for j=1:n-1
% 交叉微分
g(i,j) = abs(im(i+1,j+1)-im(i,j)) + abs(im(i+1,j)-im(i,j+1));
end
end
im2 = uint8(g);
figure,imshow(im2);
title('交叉微分(Roberts)锐化所得的图像');
figure,imshow(uint8(g+im));
title('交叉微分(Roberts)锐化并叠加原图所得的图像');
三、拓展：实现图像的压缩和解压缩
尚未完成。
参考建议：
https://blog.csdn.net/ywsydwsbn/article/details/106341759
main_Jepg.m
Fy=[16 11 10 16 24 40 51 61;
12 12 14 19 26 58 60 55;
14 13 16 24 40 57 69 56;
14 17 22 29 51 87 80 62;
18 22 37 56 68 109 103 77;
24 35 55 64 81 104 113 92;
49 64 78 87 103 121 120 101;
72 92 95 98 112 100 103 99;
];
image=imread( 'im\f11.jpg' );
[JpegCmpData, Im] =JpegCmp(image , Fy);
save('m1.mat', 'Im');
save('m2.mat' ,'JpegCmpData');
[Im2] =JpegDeCmp(JpegCmpData ,Fy);
JpegDeCmp.m
function [Im2] =JpegDeCmp(JpegCmpData ,Fy)
m=NumRow*8;
n=NumCol*8; //****5
Im2=zeros (m, n) ;
for i=1:NumRow
for j=1:NumCol
RLEcode=JpegCmpData(i, j);
F=GetRLEcode0(RLEcode) ;
F=F.*Fy;
patch=idct2(F) ;
Im2(8*(i-1)+1:8*i, 8*(j-1)+1:8*j)=patch;
end
end
Im2=uint8 (Im2);
end
function F=GetRLEcode0(RLEcode)
F=zeros(8, 8);
[t]=size (RLEcode, 1);
for i=1:t
ind=RLEcode(i) ;
F(ind:end) =RLEcode(i,2) ;
end
end
JpegCmp.m
function [JpegCmpData, Im] =JpegCmp( Im0 ,Fy)
[m0, n0, c]=size (Im0) ;
if c==3
Im0=rgb2gray (Im0) ;
end
m=m0-mod(m0, 8) ;
n=n0-mod(n0, 8) ;
Im=Im0(1:m, 1:n);
NumRow=m/8;
NumCol=n/8;
JpegCmpData{NumRow, NumCol}=[];
for i=1:NumRow
for j=1:NumCol
patch=Im(8*(i-1)+1:8*i, 8*(j-1)+1:8*j);
F=dct2(patch) ;
F=round(F./Fy) ;
RLEcode=GetRLEcode(F) ;
JpegCmpData{i, j}=RLEcode;
end
end
end
function RLEcode=GetRLEcode(F)
[m n]=size(F);
c=F(1);
RLEcode(1, 1:2)=[1,c];%建立数组RLEcode
t=2;
for i=2:m*n
if (not(F(i)==c))
RLEcode(t, 1:2)=[i,F(i)];
c=F(i) ;
t=t+1 ;
end
end
end
Ending！
更多课程知识学习记录随后再来吧！
就酱，嘎啦！
注：
人生在勤，不索何获。