基于matlab的laplace图像处理代码

基于matlab的laplace图像处理代码

主函数

function mainRGB();

clear

clc

A=imread('1.jpg');%读入图像

B=A;

imshow(A);title('原图像');

R=A(:,:,1);%RGB图像的R分量

G=A(:,:,2);%RGB图像的G分量

B=A(:,:,3);%RGB图像的B分量

figure; imshow(R);title('R');

figure; imshow(G);title('G');

figure; imshow(B);title('B');

figure('Position',[100 100 850 450],'Name','RGB模型下的锐化','NumberTitle','off' );; subplot(231), imshow(A),title('RGB图像');%读出原图像和Laplace算子高提升滤波的图像rSharp=U23(R);

subplot(232),imshow(rSharp),title('R图像');

gSharp=U23(G);

subplot(234),imshow(gSharp),title('G图像');

bSharp=U23(B);

subplot(235),imshow(bSharp),title('B图像');

A(:,:,1)=rSharp;

A(:,:,2)=gSharp;

A(:,:,3)=bSharp;

subplot(133),imshow(A),title('RGB图锐化结果图');

figure;imshow(A);title('锐化后的图像');

function N=U23(M)%Laplace算子作图像锐化函数

[dep1,wide1]=size(M);

a1=double(M);

b1=size(M);

s1=ones(size(b1)); %取同图像同大小的矩阵s1，s2

s2=ones(size(b1));

new_image1=ones(size(b1)); %取同图像同大小的矩阵new_image为输出图像的矩阵sx1=[ 0 1 0;0 -2 0 ; 0 1 0]; %Laplace算子x方向

sy1=[0 0 0;1 -2 1;0 0 0]; %Laplace算子y方向

for i=2:dep1-1 %利用Laplace算子作图像锐化

for j=2:wide1-1

x1(i,j)=0*a1(i-1,j-1) +1*a1(i-1,j) +0*a1(i-1,j+1)+0*a1(i,j-1)-2* a1(i,j) +0*a1(i,j+1)+0*a1(i+1,j-1)+a1(i+1,j)+0*a1(i+1,j+1);

y1(i,j)=0*a1(i-1,j-1)+0*2*

a1(i-1,j)+0*a1(i-1,j+1)+1*a1(i,j-1)-2*a1(i,j)+1*a1(i,j+1)+0*a1(i+1,j-1) +0*a1(i+1,j) +0*a1(i+1,j+1);

%s(i,j)=a(i,j)-x(i,j)-y(i,j);

s1(i,j)=a1(i,j)-x1(i,j)-y1(i,j);

if(s1(i,j)>0.5) %s1(i,j)的阈值判定

s1(i,j)=1;

end

new_image1(i,j)=a1(i,j)+s1(i,j);

end

end

for i=2:dep1-1 %对边缘进行操作

new_image1(i,1)=new_image1(i,2);

new_image1(i,wide1)=new_image1(i,wide1-1);

end

new_image1(1,:)=new_image1(2,:);

new_image1(dep1,:)=new_image1(dep1-1,:);

N=new_image1;