matlab版demosaic程序

clear I=imread('kodim19.png'); figure(1) imshow(I);%原始图像 [m,n,p]=size(I); I1=zeros(m,n); for i=1:1:m/2 for j=1:1:n/2 I1(i*2-1,j*2-1)=I(i*2-1,j*2-1,2); I1(i*2,j*2)=I(i*2,j*2,2); I1(i*2,j*2-1)=I(i*2,j*2-1,3); I1(i*2-1,j*2)=I(i*2-1,j*2,1); end end figure(2) imshow(I1/255);%bayer图像 imwrite(I1/255,'stripes2.png'); I2=demosaic(uint8(I1),'grbg'); figure(3) imshow(I2);%matlab自带demosaic函数图像 I3=zeros(m+2,n+2); for i=1:1:m for j=1:1:n I3(i+1,j+1)=I1(i,j); end end for i=1:1:m/2 for j=1:1:n/2 I4(i*2-1,j*2-1,2)=I3(i*2,j*2); I4(i*2-1,j*2-1,1)=(I3(i*2,j*2-1)+I3(i*2,j*2+1))/2; I4(i*2-1,j*2-1,3)=(I3(i*2-1,j*2)+I3(i*2+1,j*2))/2; I4(i*2,j*2,2)=I3(i*2+1,j*2+1); I4(i*2,j*2,1)=(I3(i*2+2,j*2+1)+I3(i*2,j*2+1))/2; I4(i*2,j*2,3)=(I3(i*2+1,j*2+2)+I3(i*2+1,j*2))/2; I4(i*2,j*2-1,3)=I3(i*2+1,j*2); I4(i*2,j*2-1,2)=(I3(i*2,j*2)+I3(i*2+2,j*2)+I3(i*2+1,j*2-1)+I3(i*2+1,j*2+1))/4; I4(i*2,j*2-1,1)=(I3(i*2,j*2-1)+I3(i*2,j*2+1)+I3(i*2+2,j*2-1)+I3(i*2+2,j*2+1))/4; I4(i*2-1,j*2,1)=I3(i*2,j*2+1); I4(i*2-1,j*2,2)=(I3(i*2-1,j*2+1)+I3(i*2+1,j*2+1)+I3(i*2,j*2)+I3(i*2,j*2+2))/4; I4(i*2-1,j*2,3)=(I3(i*2-1,j*2)+I3(i*2-1,j*2+2)+I3(i*2+1,j*2)+I3(i*2+1,j*2+2))/4; end end figure(4) imshow(I4/255);%手动插值图像