根据给定方向不同的线使图像变形

interp2

X

Y

Z

Xi

Yi

for

interp2

close all

%% Image
Z = imread(\'tire.tif\');
% The corresponding mesh
sz = size(Z);
[X,Y] = meshgrid(1:sz(2), 1:sz(1));

%% Line
lineEquation = @(c, v) (v(1)*(X-c(2))+v(2)*(Y-c(1)))/norm(v);
getLine = @(c, v) abs(lineEquation(c, v))<1/2;
% Example
c = [93, 109];
v = [2, 1];
line = getLine(c, v);

%% Circle
circleEquation = @(c, r) ((X-c(2)).^2+(Y-c(1)).^2-r^2)/r^2;
getCircle = @(c, r) abs(circleEquation(c, r))<1/r;
% Example
r = 24;
circle = getCircle(c, r);

%% Plot a sequence of line
figure;
for delta = -1:0.1:1
        v = [0.1, delta];
        line_delta = getLine(c, v);
        Z_line = Z;
        Z_line(line_delta) = 255;
        imagesc(Z_line);
        colormap(\'gray\');
        pause(0.05);
end

%% Plot examples
figure;

subplot(221);
imagesc(Z); 
title(\'Image\');

subplot(222);
Z_line = Z;
Z_line(line) = 255;
imagesc(Z_line);
title(\'Line\');

subplot(223);
Z_circle = Z;
Z_circle(circle) = 255;
imagesc(Z_circle);
title(\'Circle\');

subplot(224);
% TODO
title(\'Warped image\');

colormap(\'gray\');

[ANG, RAD] = meshgrid(0:0.01:2*pi, 0:0.5:166);
XI = c(2) + cos(ANG).*RAD;
YI = c(1) + sin(ANG).*RAD;
WARPED = interp2(X, Y, double(Z), XI, YI);

WARPED(isnan(WARPED))= max(WARPED(:));
imagesc(WARPED);
title(\'Warped image\');
colormap(\'gray\');