/*
?*?A?best?exemplar?finder.??Scans?over?the?entire?image?（using?a
?*?sliding?window）?and?finds?the?exemplar?which?minimizes?the?sum
?*?squared?error?（SSE）?over?the?to-be-filled?pixels?in?the?target
?*?patch.?
?*
?*?@author?Sooraj?Bhat
?*/
#include?“mex.h“
#include?
#include?“matrix.h“

void?bestexemplarhelper（const?int?mm?const?int?nn?const?int?m?const?int?n?
			const?double?*img?const?double?*Ip?
			const?mxLogical?*toFill?const?mxLogical?*sourceRegion?????????
			double?*best）?
{
??register?int?ijiijjii2jj2MNIJndxndx2mn=m*nmmnn=mm*nn;
??double?patchErr=0.0err=0.0bestErr=1000000000.0;

??/*?foreach?patch?*/
??N=nn-n+1;??M=mm-m+1;
??for?（j=1;?j<=N;?++j）?
??{
????J=j+n-1;
????for?（i=1;?i<=M;?++i）?
????{
??????I=i+m-1;
??????/***?Calculate?patch?error?***/
??????/*?foreach?pixel?in?the?current?patch?*/
??????for?（jj=jjj2=1;?jj<=J;?++jj++jj2）?
??????{
	for?（ii=iii2=1;?ii<=I;?++ii++ii2）?
????{
	??ndx=ii-1+mm*（jj-1）;
	??if?（!sourceRegion[ndx]）
	????goto?skipPatch;
	??ndx2=ii2-1+m*（jj2-1）;
	??if?（!toFill[ndx2]）?
??????{
	????err=img[ndx??????]?-?Ip[ndx2????];?patchErr?+=?err*err;
	????err=img[ndx+=mmnn]?-?Ip[ndx2+=mn];?patchErr?+=?err*err;
	????err=img[ndx+=mmnn]?-?Ip[ndx2+=mn];?patchErr?+=?err*err;
	??}
	}
??????}
??????/***?Update?***/
??????if?（patchErr???????{
	bestErr?=?patchErr;?
	best[0]?=?i;?best[1]?=?I;
	best[2]?=?j;?best[3]?=?J;
??????}
??????/***?Reset?***/
????skipPatch:
??????patchErr?=?0.0;?
????}
??}
}

/*?best?=?bestexemplarhelper（mmnnmnimgIptoFillsourceRegion）;?*/
void?mexFunction（int?nlhsmxArray?*plhs[]int?nrhsconst?mxArray?*prhs[]）?
{
??int?mmnnmn;
??double?*img*Ip*best;
??mxLogical?*toFill*sourceRegion;

??/*?Extract?the?inputs?*/
??mm?=?（int）mxGetScalar（prhs[0]）;/*The?value?of?the?first?real?（nonimaginary）?element?of?the?mxArray.*/
??nn?=?（int）mxGetScalar（prhs[1]）;
??m??=?（int）mxGetScalar（prhs[2]）;
??n??=?（int）mxGetScalar（prhs[3]）;
??img?=?mxGetPr（prhs[4]）;/*The?address?of?the?first?element?of?the?real?data*/
??Ip??=?mxGetPr（prhs[5]）;
??toFill?=?mxGetLogicals（prhs[6]）;/*mxGetLogicals（）The?address?of?the?first?logical?element?in?the?mxArray.?The?result?is?unspecified?if?the?mxArray?is?not?a?logical?array.*/
??sourceRegion?=?mxGetLogicals（prhs[7]）;
??
??/*?Setup?the?output?*/
??plhs[0]?=?mxCreateDoubleMatrix（41mxREAL）;
??best?=?mxGetPr（plhs[0]）;
??best[0]=best[1]=best[2]=best[3]=0.0;

??/*?Do?the?actual?work?*/
??bestexemplarhelper（mmnnmnimgIptoFillsourceRegionbest）;
}

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件??????162139??2015-03-26?10:20??20.png
?????文件??????136031??2015-03-26?10:20??201.png
?????文件???????37938??2015-03-29?20:15??benggee_7.jpg
?????文件????????2382??2015-04-06?14:14??bestexemplarhelper.asv
?????文件????????2503??2015-04-30?11:56??bestexemplarhelper.c
?????文件????????8536??2013-03-12?13:12??bestexemplarhelper.mexmaci64
?????文件????????6656??2015-03-26?10:46??bestexemplarhelper.mexw32
?????文件?????????293??2015-04-30?21:37??Evaluate.m
?????文件????????8928??2015-05-03?23:46??inpaint7.asv
?????文件????????8928??2015-05-03?23:49??inpaint7.m
?????文件?????????486??2015-04-11?14:39??plotall.m
?????文件????????1063??2004-11-29?23:04??README.txt