% function

LUT(GammaValue)

GammaValue=2

%生成查找表

for i=0:255

x=double(i);

correction=256*((x+0.5)/256)^(1/GammaValue)-0.5;

table(i+1)=uint8(correction);

end

%输出到mif

fid=fopen('gamma_correction.mif','w+');

fprintf(fid,'--GammaValue=M

\r\n',double(GammaValue));

fprintf(fid,'--%s.\r\n',datestr(now));

fprintf(fid,'%x\n',(table(1:256)))

fclose(fid)

%插值LUT 生成interp

for i=0:63

interp(1+i)=table(1+4*i)

end

fid=fopen('interpolation.mif','w+');

fprintf(fid,'--interpolation GammaValue=M

\r\n',double(GammaValue));

fprintf(fid,'--%s.\r\n',datestr(now));

fprintf(fid,'%x\n',(interp(1:64)));

%% load image

data=imread('1.jpg');

[sx,sy,sz]=size(data);

r=data(:,:,1);

g=data(:,:,2);

b=data(:,:,3);

%write R value to

r.mif

fid=fopen('r.mif','w+');

for i=1:sx

for

j=1:sy

fprintf(fid,'%d\n',r(i,j));

end

end

fclose(fid);

%write G value to

g.mif

fid=fopen('g.mif','w+');

for i=1:sx

for

j=1:sy

fprintf(fid,'%d\n',g(i,j));

end

end

fclose(fid);

%write B value to

b.mif

fid=fopen('b.mif','w+');

for i=1:sx

for

j=1:sy

fprintf(fid,'%d\n',b(i,j));

end

end

fclose(fid);

%% correcting

data_in=double(data);

[x,y,z]=size(data_in)

for i=1:x

for j=1:y

for k=1:z

temp=data_in(i,j,k);

m=rem(temp,4);

n=fix(temp/4);

if (n+2)>64

gamm(i,j,k)=interp(64);

else

switch m

case 0

gamm(i,j,k)=interp(n+1);

case 1

gamm(i,j,k)=0.75*interp(n+1)+0.25*interp(n+2);

case 2

gamm(i,j,k)=0.5*interp(n+1)+0.5*interp(n+2);

case

3

gamm(i,j,k)=0.25*interp(n+1)+0.75*interp(n+2);

end

end

%

t1=double((temp+0.5)/256);

%

t2=double(t1^(1/6.6));

%

t3=double(256*t2-0.5);

%

gamm(i,j,k)=uint8(t3);

end

end

end

pic=cat(2,data,gamm);

figure,imshow

(pic)