function out = simcomp(in1, in2, edgein1, edgein2)

a=edgein1;
b=edgein2;

out = imgcomp3d(in1(a(1):a(3),a(2):a(4),:), in2(b(1):b(3),b(2):b(4),:));