雙線性插值法縮放圖像函數MATLAB程序

雙線性插值法縮放圖像函數MATLAB程序(精)雙線性插值法縮放圖像函數MATLAB程序(精)4/4雙線性插值法縮放圖像函數MATLAB程序(精)DefineafunctionforshrinkingandzoomingimagesbybilinearinterpolationIistheoriginalimageandNisthefactorforshrinkingfirstlyandzoominglaterfunctioninterpolation_shrink_zoom_image(I,NFirsttoshrinktheoriginalimageTomeasurethesizeoforiginalimageI[P,Q]=size(I;RebuildtheoriginalimageItoI_TEMPforbilinearinterpolationI_TEMP=zeros(P+2,Q+2;I_TEMP(2:P+1,2:Q+1=I;I_TEMP(1,2:Q+1=I(1,:;I_TEMP(P+2,2:Q+1=I(P,:;I_TEMP(2:P+1,1=I(:,1;I_TEMP(2:P+1,Q+2=I(:,Q;I_TEMP(1,1=I(1,1;I_TEMP(1,P+2=I(1,P;I_TEMP(P+2,1=I(P,1;I_TEMP(P+2,Q+2=I(P,Q;I_TEMP=double(I_TEMP;TofixthesizeofshrunkenimageRisthenumberofrowsandCisthenumberofcolumnsR=floor(P/N;C=floor(Q/N;DefinetheshrunkenimageI_S=zeros(R,C;Forbilinearinterpolationf(x,yin4pointsf(x1,y1,f(x2,y1,f(x1,y2,f(x2,y2%f(x,yf(x1,y1*(x2-x*(y2-y/((x2-x1*(y2-y1+f(x2,y1*(x-x1*(y2-y/((x2-x1*(y2-y1+f(x1,y2*(x2-x*(y-y1/((x2-x1*(y2-y1+f(x2,y2*(x-x1*(y-y1/((x2-x1*(y2-y1Forpixelmatrix,x2=x1+1,y2=y1+1,x=x1+a,y=y1+bLetx1=u,y1=v,thenx2=u+1,y2=v+1,x=u+a,y=v+bThenf(u+a,v+b=f(u,v*(1-a*(1-b+f(u+1,v*a*(1-b+f(u,v+1*(1-a*b+f(u+1,v+1*a*bMethod:fromoriginalimageIgetsuitablegrayvaluebysurrounding4pixels,thenfillthegrayvalueintocertainlocationinshrunkenimageI_S.fori=1:R;forj=1:C;Step1:getx,y,u,v,a,b.ThesizeofItosizeofI_Sratio=N,(x,yispixel"in"Iand(i,jispixelinI_S,thenx/i=y/j=N.x=i*N;y=j*N;u=floor(x;v=floor(y;a=x-u;b=y-v;Step2:consideringthesituationonboundofI,usingI_TEMPtogetsuitablevalueoff(x,ybybilinearinterpolationI_S(i,j=I_TEMP(u+1,v+1*(1-a*(1-b+I_TEMP(u+2,v+1*a*(1-b+I_TEMP(u+1,v+2*(1-a*b+I_TEMP(u+2,v+2*a*b;end;end;DisplaytheoriginalimageIandtheshrunkenimageI_SI_S=uint8(I_S;figure;imshow(I,[];title(['Originalimage:',num2str(P,'*',num2str(Q];figure;imshow(I_S,[];title(['Interpolationshrinkingimage:',num2str(R,'*',num2str(C];SecondtozoomtheshrunkenimageR*CisthesizeofshrunkenimageasshownbeforeRebuildtheshrunkenimageI_StoI_TEMPforbilinearinterpolationI_TEMP=zeros(R+2,C+2;I_TEMP(2:R+1,2:C+1=I_S;I_TEMP(1,2:C+1=I_S(1,:;I_TEMP(R+2,2:C+1=I_S(R,:;I_TEMP(2:R+1,1=I_S(:,1;I_TEMP(2:R+1,C+2=I_S(:,C;I_TEMP(1,1=I_S(1,1;I_TEMP(1,R+2=I_S(1,R;I_TEMP(R+2,1=I_S(R,1;I_TEMP(R+2,R+2=I_S(R,C;DefinethezoomedimageI_Z=zeros(R*N,C*N;Zoomingtheshrunkenimagebybilinearinterpolationfori=1:R*N;forj=1:C*N;x=i/N;y=j/N;u=floor(x;v=floor(y;a=x-u;b=y-v;Step2:consideringthesituationonboundofI,usingI_TEMPtogetsuitablevalueoff(x,ybybilinearinterpolationI_Z(i,j=I_TEMP(u+1,v+1*(1-a*(1-b+I_TEMP(u+2,v+1*a*(1-b+end;end;DisplaytheoriginalimageIandtheshrunkenimageI_SI_Z=uint8(I_Z;figure;imshow(I,[];tit