超分辨率插值法matlab代码

function [original,new_img] = imblizoom(original,zmf) if ~exist('original','var')||isempty(original) error('输入图像为空') end if ~exist('zmf','var')||isempty(zmf)||numel(zmf)~=1 error('zmf为空') end if isstr(original) [original,M]=imread(original) end if zmf<=0 error('倍数应该大于0') end [height,width,channel]=size(original) new_height=round(height*zmf) new_width=round(width*zmf) new_img=zeros(new_height,new_width,channel) img_scale=zeros(height+2,width+2,channel) img_scale(2:height+1,2:width+1,:)=original img_scale(1,2:width+1,:) = original(1,:,:) img_scale(height+2,2:width+1,:) = original(height,:,:) img_scale(2:height+1,1,:) = original(:,1,:) img_scale(2:height+1,width+2,:) = original(:,width,:) img_scale(1,1,:) = original(1,1,:) img_scale(1,width+2,:) = original(1,width,:) img_scale(height+2,1,:) = original(height,1,:) img_scale(height+2,width+2,:) = original(height,width,:) for zj = 1:new_width % 对图像进行按列逐元素扫描 for zi = 1:new_height % (zi,zj)表示在新图中的坐标，(ii,jj)表示在原图中的坐标 % 注意：(ii,jj)不一定是整数 ii = (zi-1)/zmf; jj = (zj-1)/zmf; i = floor(ii); j = floor(jj); % 向下取整得到在原图中坐标的整数部分 u = ii - i; v = jj - j; % 得到在原图中坐标的小数部分 i = i + 1; j = j + 1; new_img(zi,zj,:) = (1-u)*(1-v)*img_scale(i,j,:) + u*(1-v)*img_scale(i,j+1,:)... + (1-u)*v*img_scale(i+1,j,:) + u*v*img_scale(i+1,j+1,:); end end new_img = uint8(new_img); end