clc;
clear;
%???Global?Poisson?METHOD
%???By?Guan?Li?May?2004
A=imread（‘cat.jpg‘）;
imshow（A）;
hold?on;

Threshold?=?1;
UniBack=[0
????0
???255];

%?choose?background?object?and?B?results?into?the?Foreground?object?and?the
%?Boundary?Region
%?set（handles.text5?‘String‘?‘Please?choose?the?outer?line?of?the?boundary?region.?Single?left?click?the?mouse?to?specify?vertice.Double?left?click?or?single?right?click?to?finish‘）;
OuterBW?=?roipoly（A）;
B（::1）?=?immultiply（A（::1）OuterBW）;
B（::2）?=?immultiply（A（::2）OuterBW）;
B（::3）?=?immultiply（A（::3）OuterBW）;
Back（::1）?=?immultiply（A（::1）~OuterBW）;
Back（::2）?=?immultiply（A（::2）~OuterBW）;
Back（::3）?=?immultiply（A（::3）~OuterBW）;
imshow（B）;

%?choose?foreground?object?and?C?results?into?a?Circle?of?Interest
%?set（handles.text5?‘String‘?‘Please?choose?the?inner?line?of?the?boundary?region.?Single?left?click?the?mouse?to?specify?vertice.Double?left?click?or?single?right?click?to?finish‘）;
InnerBW?=?roipoly（B）;
C（::1）?=?immultiply（B（::1）~InnerBW）;
C（::2）?=?immultiply（B（::2）~InnerBW）;
C（::3）?=?immultiply（B（::3）~InnerBW）;
Fore（::1）?=?immultiply（A（::1）InnerBW）;
Fore（::2）?=?immultiply（A（::2）InnerBW）;
Fore（::3）?=?immultiply（A（::3）InnerBW）;
imshow（C）;


hold?off;


%?set（handles.text5?‘String‘?‘Calculating?Raw?Fore?&?BackGround......................‘）;
RawAlpha?=?（double（OuterBW）+double（InnerBW））/2;
%imshow（uint8（RawAlpha*255））;

I?=?double（A）;

%?Method?of?“averaging“?to?get?raw?fore?and?background?colors
for?i?=?1:size（RawAlpha1）
?????for?j?=?1:size（RawAlpha2）
?????????RawFore（ij1）=0;?????????????????????
?????????RawFore（ij2）=0;
?????????RawFore（ij3）=0;
?????????RawBack（ij1）=0;?????????????????????
?????????RawBack（ij2）=0;
?????????RawBack（ij3）=0;
?????????if?RawAlpha（ij）?==?1
?????????????RawFore（ij:）=Fore（ij:）;
?????????end;
?????????if?RawAlpha（ij）?==?0
?????????????RawBack（ij:）=Back（ij:）;
?????????end;
?????????if?RawAlpha（ij）?==?0.5
?????????????r=1;
?????????????while?1
?????????????????TestFore?=?InnerBW（max（i-r1）:min（i+rsize（A1））max（j-r1）:min（j+rsize（A2）））;
?????????????????if?size（find（TestFore）1）??
?????????????????????[i1?j1]?=?find（TestFore）;
?????????????????????i2?=?i1?+?max（i-r1）?-1;
?????????????????????j2?=?j1?+?max（j-r1）?-1;
?????????????????????Rs?=?double（Fore（i2j2:））;
?????????????????????Ds（:1）=diag（Rs（::1））;
?????????????????????Ds（:2）=diag（Rs（::2））;
?????????????????????Ds（:3）=diag（Rs（::3））;?????????????????
?????????????????????if?size（Ds1）?==?1
?????????????????????????RawFore（ij:）?=Ds;
?????????????????????else
?????????????????????????RawFore（ij:）?=?sum（Ds）?/?size（Ds1）;
?????????????????????end;
?????????????????????clear?Ds;
?????????????????????break;
?????????????????else
?????????????????????r=r+1;
?????????????????end;
?????????????end;
?????????????r=1;
?????????????while?1
??????