matlab练习程序（矩形变换为单连通形状）

clear all;close all;clc;

img=imread(‘lena.jpg‘);
mask_line=imread(‘mask.bmp‘)>50;

[h w]=size(img);
imgn=zeros(h,w);

ind=find(mask_line==1);
[lineY lineX]=ind2sub([h w],ind);   %模板形状像素的位置
cenY=round(mean(lineY));            %模板形状中心
cenX=round(mean(lineX));

lineAng=atan2(lineY-cenY,lineX-cenX);        %模板线上每个点到形状中心弧度
lineR=sqrt((lineY-cenY).^2+(lineX-cenX).^2); %模板线上每个点到形状中心距离

ang_ul=atan2(1-cenY,1-cenX);        %图像四个点与模板形状中心弧度
ang_ur=atan2(1-cenY,w-cenX);
ang_dl=atan2(h-cenY,1-cenX);
ang_dr=atan2(h-cenY,w-cenX);

for y=1:h
    for x=1:w
        r=sqrt((x-cenX)^2+(y-cenY)^2);
        ang=atan2(y-cenY,x-cenX);

        if ang>=-pi && ang<ang_ul           %分别求八个象限图像边界到模板中心距离
            R=abs(cenX/cos(ang));           %应该可以简化，这样判断似乎繁琐了
        elseif ang>=ang_ul && ang<-pi/2
            R=abs(cenY/cos(pi/2-ang));        
        elseif ang>=-pi/2 && ang<ang_ur
            R=abs(cenY/cos(pi/2+ang)); 
        elseif ang>=ang_ur && ang<0
            R=abs((w-cenX)/cos(pi-ang)); 
        elseif ang>=0 && ang<ang_dr
            R=abs((w-cenX)/cos(ang));           
        elseif ang>=ang_dr && ang< pi/2
            R=abs((h-cenY)/cos(pi/2-ang));    
        elseif ang>=pi/2 && ang<ang_dl
            R=abs((h-cenY)/cos(pi/2+ang));             
        elseif ang>=ang_dl && ang<=pi
            R=abs(cenX/cos(pi-ang));              
        end
            
        [junk index]=min(abs(ang-lineAng));     %距离最近的弧度
        newR=lineR(index);
        
        ss=R*r/newR;        %缩减扩展比例系数
        xx=round((ss*cos(ang)+cenX));
        yy=round((ss*sin(ang)+cenY));
       
        if xx>=1 && xx<=w && yy>=1 && yy<=h
            imgn(y,x)=img(yy,xx);
        end
                      
    end
end

imshow(img);
figure;imshow(mask_line);
figure;imshow(imgn,[])