%A=imread('shiyan5.png'); %读入图像
I=imread('shiyan5.jpg'); %读入图像
A=rgb2gray(I);
[m,n]=size(A);
figure(4)
subplot(2,2,1);imshow(A);title('原始图像');
subplot(2,2,2);imhist(A);title('原始图像直方图');
% 构造变换函数
k0=30;
k1=300;
k=1:256;
yb=256*(1-(1/(256-k0).^2)*(k-256).^2);
yg=256*(1-(1/((k1-k0)/2).^2)*(k-k0-(k1-k0)/2).^2);
yr=256*(1-(1/((200-k0)/2).^2)*(k-k0-(200-k0)/2).^2)+(k/90).^7-(k/150).^13;
yr(1:k0)=0;yg(1:k0)=0;yb(1:k0)=0;
subplot(2,2,3);plot(k,yr,'r',k,yg,'g',k,yb,'b');title('伪彩色变换曲线');
for j=1:n
for i=1:m
if(A(i,j)>0)
CC(i,j,1)=uint8(yr(A(i,j)));
CC(i,j,2)=uint8(yg(A(i,j)));
CC(i,j,3)=uint8(yb(A(i,j)));
end
end
end
subplot(2,2,4);imshow(CC);title('伪彩色图像');
评论0