%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