先只贴出代码
clear all;close all;clc;
for k=1:6
I=imread(['H:\mwtmodel\ths_featuremap\Resize_2\Original_' num2str(k) '.tiff']);
% figure;
% imshow(I);
I=double(I);
[m,n]=size(I);
L=256;
for i=1:m
for j=1:n
if I(i,j)<L/4
R(i,j)=0;
G(i,j)=4*I(i,j);
B(i,j)=L;
else if I(i,j)<=L/2
R(i,j)=0;
G(i,j)=L;
B(i,j)=-4*I(i,j)+2*L;
else if I(i,j)<=3*L/4
R(i,j)=4*I(i,j)-2*L;
G(i,j)=L;
B(i,j)=0;
else
R(i,j)=L;
G(i,j)=-4*I(i,j)+4*L;
B(i,j)=0;
end
end
end
end
end
for i=1:m
for j=1:n
G2C(i,j,1)=R(i,j);
G2C(i,j,2)=G(i,j);
G2C(i,j,3)=B(i,j);
end
end
G2C=G2C/256;
imwrite(G2C,['C:\Users\Administrator\Desktop\ths_featuremap\Resize_2\Original_' num2str(k) '.tiff']);
end