图片处理（旋转缩放）matlab程序

%%%%zoom，>1放大，<1缩小%%%% %先处理横向缩放，产生pic_temp0 %再处理纵向缩放，产生pic_temp1，生成最终图形 %缩小采用临近点均值算法，放大采用线性插值 clear; clc; tic m=2.0;%水平方向伸缩系数 n=2.0;%竖直方向伸缩系数 pic_init=imread('E:\MATLAB\zoom\example.jpg'); temp0=size(pic_init); a=temp0(2); b=temp0(1); xm=a*m; ym=b*n; pic_out=zeros(uint16(ym),uint16(xm),3); %先处理横向缩放，输出文件为pic_temp0 if(m<=1)%横向缩小 pic_temp0=zeros(uint16(b),uint16(xm),3);%横向缩放存为pic_temp0 for y0=1:1:b for x0=1:1:xm temp1=zeros(1,1,3); for x0temp=uint16(1+((x0-1)/m)):1:uint16(x0/m)%缩小求和 temp1=temp1+double(pic_init(y0,x0temp,:)); end length=double((uint16(x0/m)-uint16(1+((x0-1)/m)))+1); temp1=temp1/length;%缩小取平均 pic_temp0(y0,x0,:)=temp1; end end else%(m>1)横向放大 pic_temp0=zeros(uint16(b),uint16(xm),3); for y0=1:1:b for x0=1:1:xm temp1=zeros(1,1,3); locate=floor((x0-1)/m)+1;%被放大的点对应原图左侧最近点坐标 rate=(x0-1)/m+1-locate;%放大后的点分原图相邻两点比例，0为与左重合 temp1=(1-rate)*double(pic_init(y0,locate,:)); if((locate+1)<=a) temp1=temp1+rate*double(pic_init(y0,locate+1,:));%可能越界所以加if end pic_temp0(y0,x0,:)=temp1; end end end %再处理纵向缩放,输出文件为pic_temp1 if(n<=1)%纵向缩小 pic_temp1=zeros(uint16(ym),uint16(xm),3);%纵向缩放存为pic_temp0 for y0=1:1:ym for x0=1:1:xm temp1=zeros(1,1,3); for y0temp=uint16(1+((y0-1)/n)):1:uint16(y0/n) temp1=temp1+double(pic_temp0(y0temp,x0,:)); end length=double((uint16(y0/n)-uint16(1+((y0-1)/n)))+1); temp1=temp1/length; pic_temp1(y0,x0,:)=temp1; end end else%(n>1)纵向放大 pic_temp1=zeros(uint16(ym),uint16(xm),3); for y0=1:1:ym for x0=1:1:xm temp1=zeros(1,1,3); locate=floor((y0-1)/n)+1;%被放大的点对应原图左侧最近点坐标 rate=(y0-1)/n+1-locate;%放大后的点分原图相邻两点比例，0为与左重合 temp1=(1-rate)*double(pic_temp0(locate,x0,:)); if((locate+1)<=b) temp1=temp1+rate*double(pic_temp0(locate+1,x0,:));%可能越界所以加if end pic_temp1(y0,x0,:)=temp1; end end end pic_out=uint8(pic_temp1); imshow(pic_out) toc