script6

clear
close all
load filelist
load script2
load script6

n=9;
filename = char(files(n));
rgb1 = imread(filename);
imshow(rgb1);

%pti = ginput(1)
k=1;
%pts(1k,:) = pti;

hold on
plot(pts(k,1),pts(k,2),'yx');
hold off

g1 = rgb2gray(im2double(rgb1));

Warning: Image is too big to fit on screen; displaying at 50%

n=10;
filename = char(files(n));
rgb2 = imread(filename);
imshow(rgb2);

%pti = ginput(1)
k=2;
%pts(k,:) = pti;

hold on
plot(pts(k,1),pts(k,2),'yx');
hold off

g2 = rgb2gray(im2double(rgb2));

n=11;
filename = char(files(n));
rgb3 = imread(filename);
imshow(rgb3);

%pti = ginput(1)
k=3;
%pts(k,:) = pti;
hold on
plot(pts(k,1),pts(k,2),'yx');
hold off

g3 = rgb2gray(im2double(rgb3));

n=12;
filename = char(files(n));
rgb4 = imread(filename);
imshow(rgb4);

%pti = ginput(1)
k=4;
%pts(k,:) = pti;
hold on
plot(pts(k,1),pts(k,2),'yx');
hold off

g4 = rgb2gray(im2double(rgb4));
%save script6 pts

d1 = g2-g1;
imshow((d1+1)/2);

d2 = g3-g2;
imshow((d2+1)/2);

d3 = g4-g3;
imshow((d3+1)/2);

[xpts ypts] = tformfwd(tf2,pts(:,1),pts(:,2));
xp = xpts/scale;
yp = ypts/scale;

dx = diff(xp)
dy = diff(yp)

dist = sqrt(dx.^2+dy.^2)

% four tiles @ 9 inches/tile = 3 ft.
%
tile_size = dist*12/4;
fprintf('tile size %g %g %g inches\n',tile_size);

dx =

   -0.9172
   -0.9644
   -1.0136


dy =

   -2.4632
   -2.5495
   -2.7776


dist =

    2.6284
    2.7258
    2.9568

tile size 7.88527 8.17733 8.87028 inches