clear
clc
close all

pic = imread('SS1.jpg');

pic2 = imread('DMS.jpg');

ins = pic2(600:900, 1000:1300, [2 3 1]);
[r, c, cl] = size(ins)
for ir = 1:r
    for ic = 1:c
        alf(ir, ic, 1:3) = ir / r;
    end
end
comp = picInsert( pic, ins, 500, 700, alf );

image(comp)



cr = r/2;
cc = c/2;
dmx = min([cr cc]);
for ir = 1:r
    for ic = 1:c
        dr = ir - cr;
        dc = ic - cc;
        r = sqrt(dr.^2 + dc.^2);
        alf(ir, ic, 1:3) = (dmx - r) / (0.75 * dmx);
    end
end
alf(alf < 0) = 0;
alf(alf > 1) = 1;
figure
n = 50;
for ix = 1:n
    alff = alf * ix / n;
    image(picInsert( pic, ins, 300, 400, alff) )
    pause(0.001)
end