matlab图像水波纹

%matlab图像水波纹
clear all;close all;clc;

img=imread('lena.jpg');
 
img = im2double(img); 
[h, w, c]=size(img)
ratio = 600/(h+w); 
img = imresize(img, ratio); 
[h, w, c]=size(img)

%figure; imshow(img); 
   
wave=[10,350]; %[幅度，周期]
newh=h+2*wave(1);
neww=w+2*wave(1);
rot=0;

imgn=zeros(newh,neww,3);    
    
for i=1:50
   
    rot=rot+0.6;
    for y=1:newh
        for x=1:neww

            yy=round((y-wave(1))+(wave(1)*cos(2*pi/wave(2)*x+rot  )));    %依然是逆变换
            xx=round((x-wave(1))+(wave(1)*cos(2*pi/wave(2)*y+rot )));

           if yy>=1 && yy<=h && xx>=1 && xx<=w                
                imgn(y,x, 1)=img(yy,xx, 1);
                imgn(y,x, 2)=img(yy,xx, 2);
                imgn(y,x, 3)=img(yy,xx, 3);
           end

        end
    end
   
     imshow(imgn,[]);  
     
 
end