a=imread('foreground.jpg');
b=imread('background.jpg');
da=double(a);
db=double(b);
D=imabsdiff(a,b);
r=zeros(240,320);
AR = da(:,:,1);
AG = da(:,:,2);
AB = da(:,:,3);
BR = db(:,:,1);
BG = db(:,:,2);
BB = db(:,:,3);
mask = D(:,:,1)>20 & ...
AR~=0 & AG~=0 & AB~=0 & ...
BR./AR<4 & BR./AR>1.5 & ...
BG./AG<2.8 & BG./AG>1.3 & ...
BB./AB<2.05 & BB./AB>1.14 & ...
BB./AB<BR./AR & BB./AB<BG./AG & BG./AG<BR./AR & ...
abs(AR.
/(AR
+AG
+AB
)-BR.
/(BR
+BG
+BB
))<0.129 & ...
abs(AG.
/(AR
+AG
+AB
)-BG.
/(BR
+BG
+BB
))<0.028 & ...
abs(AB.
/(AR
+AG
+AB
)-BB.
/(BR
+BG
+BB
))<0.143; r(mask)=0;
AR(mask) = 255;
AG(mask) = 255;
AB(mask) = 255;
h = cat(3,AR,AG,AB);
imshow(h);
im=h-a;
imshow(im);
YT1pbXJlYWQoJ2ZvcmVncm91bmQuanBnJyk7CmI9aW1yZWFkKCdiYWNrZ3JvdW5kLmpwZycpOwpkYT1kb3VibGUoYSk7CmRiPWRvdWJsZShiKTsKRD1pbWFic2RpZmYoYSxiKTsKcj16ZXJvcygyNDAsMzIwKTsKQVIgPSBkYSg6LDosMSk7CkFHID0gZGEoOiw6LDIpOwpBQiA9IGRhKDosOiwzKTsKQlIgPSBkYig6LDosMSk7CkJHID0gZGIoOiw6LDIpOwpCQiA9IGRiKDosOiwzKTsKCm1hc2sgPSBEKDosOiwxKT4yMCAmIC4uLgogICAgICAgQVJ+PTAgJiBBR349MCAmIEFCfj0wICYgLi4uCiAgICAgICBCUi4vQVI8NCAmIEJSLi9BUj4xLjUgJiAuLi4KICAgICAgIEJHLi9BRzwyLjggJiBCRy4vQUc+MS4zICYgLi4uCiAgICAgICBCQi4vQUI8Mi4wNSAmIEJCLi9BQj4xLjE0ICYgLi4uCiAgICAgICBCQi4vQUI8QlIuL0FSICYgQkIuL0FCPEJHLi9BRyAmIEJHLi9BRzxCUi4vQVIgJiAuLi4KICAgICAgIGFicyhBUi4vKEFSK0FHK0FCKS1CUi4vKEJSK0JHK0JCKSk8MC4xMjkgJiAuLi4KICAgICAgIGFicyhBRy4vKEFSK0FHK0FCKS1CRy4vKEJSK0JHK0JCKSk8MC4wMjggJiAuLi4KICAgICAgIGFicyhBQi4vKEFSK0FHK0FCKS1CQi4vKEJSK0JHK0JCKSk8MC4xNDM7CnIobWFzayk9MDsKQVIobWFzaykgPSAyNTU7CkFHKG1hc2spID0gMjU1OwpBQihtYXNrKSA9IDI1NTsKCmggPSBjYXQoMyxBUixBRyxBQik7Cmltc2hvdyhoKTsKaW09aC1hOwppbXNob3coaW0pOw==