/*

??A????????=?area（I??y??x??h??w）;

??To?compile
??----------

??mex??-g?-output?area.dll?area.c

??mex??-f?mexopts_intel10.bat?-output?area.dll?area.c


??Example?1
??---------

??clear?close?all
??load?viola_24x24

??I?????????=?X（:??:??1）;
??II????????=?cumsum（cumsum（double（I）??1）??2）;


??y?????????=?2;
??x?????????=?4;
??h?????????=?10;
??w?????????=?11;
??A1????????=?area（I??y??x??h??w）;
??A2????????=?sum（sum（I（y:y+h-1x:x+w-1）））;
?%?A2????????=?II（y+hx+w）?+?II（yx）?-?（II（y+hx）+II（yx+w））;



*/

#include?
#include?

/*--------------------------------------------------------------------------------------------------------------?*/

/*?Function?prototypes?*/

void?MakeIntegralImage（unsigned?char?*?unsigned?int??*?int??int??unsigned?int?*）;
unsigned?int?Area（unsigned?int?*??int??int??int??int??int?）;

/*--------------------------------------------------------------------------------------------------------------?*/
void?mexFunction（?int?nlhs?mxArray?*plhs[]??int?nrhs?const?mxArray?*prhs[]?）
{
	unsigned?char?*I;
	const?int?*dimsI;
	int?numdimsI;
	unsigned?int?*II??*Itemp;
	int?y??x??h??w;
	double?ai;
	double?*A;
????int?i??Ny??Nx??NyNx??N?=?1??tempint;

????/*?Input?1??*/
	????
????if?（（nrhs?>?0）?&&?!mxIsEmpty（prhs[0]）?&&?mxIsUint8（prhs[0]））
????{????????
		dimsI????=?mxGetDimensions（prhs[0]）;
????????numdimsI?=?mxGetNumberOfDimensions（prhs[0]）;
		
		I????????=?（unsigned?char?*）?mxGetData（prhs[0]）;
		Ny???????=?dimsI[0];
		Nx???????=?dimsI[1];
		if（numdimsI?>?2）
		{
			N????=?dimsI[2];
		}		???
????}
	else
	{
		mexErrMsgTxt（“I?must?be?（Ny?x?Nx?x?N）?in?UINT8?format“）;
	}

	NyNx???????????????????????=?Ny*Nx;

	y??????????????????????????=?mxGetScalar（prhs[1]）;
	x??????????????????????????=?mxGetScalar（prhs[2]）;
	h??????????????????????????=?mxGetScalar（prhs[3]）;
	w??????????????????????????=?mxGetScalar（prhs[4]）;

????
	II?????????????????????????=?（unsigned?int?*）?malloc（NyNx*sizeof（unsigned?int））;
	Itemp??????????????????????=?（unsigned?int?*）?malloc（NyNx*sizeof（unsigned?int））;

	MakeIntegralImage（I??II??Nx??Ny??Itemp）;

	ai??????????????????????????=?（double）Area（II??x?-?1??y?-?1??w??h??Ny）;


????plhs[0]????????????????????=?mxCreateDoubleMatrix（1??1??mxREAL）;
????A??????????????????????????=?mxGetPr（plhs[0]）;

	A[0]???????????????????????=?ai;

	free（II）;
	free（Itemp）;
		
}

/*-----------------------------------------------------------------------------------------------------------------------------------------?*/
unsigned?int?Area（unsigned?int?*II??int?x??int?y??int?w??int?h??int?Ny）
{
	int?h1?=?h-1??w1?=?w-1??x1?=?x-1?y1?=?y-1;

	if（?（x?==?0）?&&?（y==0））
	{	
		return?（II[h1?+?w1*Ny]）;	
	}
	if（?x==0?）?
	{
		return（II[（y+h1）?+?w1*Ny]?-?II[y1?+?w1*Ny]）;
	}
	if（?y==0?）?
	{
		return（II[h1?+?（x+w1）*Ny]?-?II[h1?+?x1*Ny]）;	
	}
	else
	{	
		return?（II[（y+

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2011-04-16?23:16??fdtool\
?????文件????????4366??2011-04-02?22:06??fdtool\area.c
?????文件????????1544??2010-12-21?16:40??fdtool\auroc.m
?????文件????????1583??2009-05-11?02:35??fdtool\basicroc.m
?????文件?????????702??2007-02-23?14:30??fdtool\blas.h
?????文件???????16460??2007-02-23?14:30??fdtool\blasp.h
?????文件????????1673??2011-03-23?13:18??fdtool\build_negatives.m
?????文件???????23528??2011-04-13?14:11??fdtool\chlbp.c
?????文件???????17882??2011-03-09?06:21??fdtool\chlbp_adaboost_binary_predict_cascade.c
?????文件???????15801??2011-03-09?06:35??fdtool\chlbp_adaboost_binary_train_cascade.c
?????文件???????18720??2011-03-09?06:24??fdtool\chlbp_gentleboost_binary_predict_cascade.c
?????文件???????21224??2011-03-09?06:28??fdtool\chlbp_gentleboost_binary_train_cascade.c
?????文件?????????115??2009-03-09?14:09??fdtool\d2uint8_image.m
?????文件????????1205??2009-02-19?18:52??fdtool\daxpy.c
?????文件????????1211??2009-02-19?18:52??fdtool\ddot.c
?????文件???????21784??2011-02-14?19:46??fdtool\demo_chlbp.m
?????文件????????2394??2011-04-16?18:53??fdtool\demo_detector_haar.m
?????文件????????1697??2011-04-04?17:48??fdtool\demo_detector_hmblbp.m
?????文件????????4349??2011-02-15?19:23??fdtool\demo_fine_threshold.m
?????文件???????14784??2011-02-14?23:45??fdtool\demo_haar.m
?????文件????????7021??2011-02-16?04:38??fdtool\demo_haar_mblbp_training.m
?????文件????????2087??2011-02-20?21:19??fdtool\demo_integral_histogram.m
?????文件???????22820??2011-03-20?21:50??fdtool\demo_mblbp.m
?????文件????????5577??2011-02-20?14:25??fdtool\demo_mblbp_variant_training.m
?????文件???????16612??2011-03-11?05:26??fdtool\demo_mlhmslbp_spyr_svm_training.m
?????文件????????1771??2011-02-15?19:35??fdtool\demo_type_cascade_scaling_vs_interp.m
?????文件???????53842??2011-03-24?04:05??fdtool\detector_haar.c
?????文件???????44794??2011-03-24?04:05??fdtool\detector_mblbp.c
?????文件???????77503??2011-03-24?16:47??fdtool\detector_mlhmslbp_spyr.c
?????文件?????????748??2009-05-09?16:40??fdtool\display_database.m
?????文件????????1306??2009-02-19?18:52??fdtool\dnrm2.c
............此處省略130個文件信息