matlab图像扭曲波动球形变换

matlab图像扭曲波动球形变换，利用算法实现，不是调用库函数实现

function fluctuate（TxTyAxAy）
I=imread（‘im.bmp‘）;
imshow（I）;
title（‘原图‘）;
arr=size（I）;
pi=3.14;

figure;
tx=Tx;
ty=Ty;
ax=Ax;
ay=Ay;
J=zeros（arr（1）arr（2）arr（3））;

for m=1:arr（1）
    for n=1:arr（2）
        mm=abs（uint16（m+ax*sin（2*pi*n/tx）））;
        nn=abs（uint16（n+ay*sin（2*pi*m/ty）））;
        if（ mm>arr（1）||mm<=0 ）
                mm=1;
            end
            if（nn>arr（2） ||nn<=0 ） 
                nn=1;
            end
        for k=1:3
            J（mnk）=I（mmnnk）;
        end
    end
end
imshow（uint8（J））;
title（‘波动变换‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       5311  2011-03-21 17:20  第二次上机0\aaa.JPG

     文件        577  2011-03-21 16:55  第二次上机0\fluctuate.asv

     文件        577  2011-03-21 16:55  第二次上机0\fluctuate.m

     文件     127154  2011-03-19 20:33  第二次上机0\im.bmp

     文件       1009  2011-03-20 20:52  第二次上机0\spherical.asv

     文件        946  2011-03-21 18:14  第二次上机0\spherical.m

     文件       1116  2011-03-21 16:45  第二次上机0\warping.asv

     文件       1122  2011-03-21 17:26  第二次上机0\warping.m

     目录          0  2011-03-21 17:21  第二次上机0

----------- ---------  ---------- -----  ----

               137812                    9