/*????alg.c
?*
?*????Detect?changes?in?a?video?stream.
?*????Copyright?2001?by?Jeroen?Vreeken?（pe1rxq@amsat.org）
?*????This?software?is?distributed?under?the?GNU?public?license?version?2
?*????See?also?the?file?‘COPYING‘.
?*
?*/
#include?“motion.h“
#include?“alg.h“

#ifdef?__MMX__
#define?HAVE_MMX
#include?“mmx.h“
#endif

#define?MAX2（x?y）?（（x）?>?（y）???（x）?:?（y））
#define?MAX3（x?y?z）?（（x）?>?（y）???（（x）?>?（z）???（x）?:?（z））?:?（（y）?>?（z）???（y）?:?（z）））

/*?locate?the?center?and?size?of?the?movement.?*/
void?alg_locate_center_size（struct?images?*imgs?int?width?int?height?struct?coord?*cent）
{
????unsigned?char?*out?=?imgs->out;
????int?*labels?=?imgs->labels;
????int?x?y?centc?=?0?xdist?=?0?ydist?=?0;

????cent->x?=?0;
????cent->y?=?0;
????cent->maxx?=?0;
????cent->maxy?=?0;
????cent->minx?=?width;
????cent->miny?=?height;

????/*?If?Labeling?enabled?-?locate?center?of?largest?labelgroup?*/
????if?（imgs->labelsize_max）?{
????????/*?Locate?largest?labelgroup?*/
????????for?（y?=?0;?y?????????????for?（x?=?0;?x?????????????????if?（*（labels++）&32768）?{
????????????????????cent->x?+=?x;
????????????????????cent->y?+=?y;
????????????????????centc++;
????????????????}
????????????}
????????}
????}?else?{
????????/*?Locate?movement?*/
????????for?（y?=?0;?y?????????????for?（x?=?0;?x?????????????????if?（*（out++））?{
????????????????????cent->x?+=?x;
????????????????????cent->y?+=?y;
????????????????????centc++;
????????????????}
????????????}
????????}
????}

????if?（centc）?{
????????cent->x?=?cent->x?/?centc;
????????cent->y?=?cent->y?/?centc;
????}
????
????/*?Now?we?find?the?size?of?the?Motion?*/

????/*?First?reset?pointers?back?to?initial?value?*/
????centc?=?0;
????labels?=?imgs->labels;
????out?=?imgs->out;

????/*?If?Labeling?then?we?find?the?area?around?largest?labelgroup?instead?*/
????if?（imgs->labelsize_max）?{
????????for?（y?=?0;?y?????????????for?（x?=?0;?x?????????????????if?（*（labels++）&32768）?{
????????????????????if?（x?>?cent->x）
????????????????????????xdist?+=?x?-?cent->x;
????????????????????else?if?（x?x）
????????????????????????xdist?+=?cent->x?-?x;

????????????????????if?（y?>?cent->y）
????????????????????????ydist?+=?y?-?cent->y;
????????????????????else?if?（y?y）
????????????????????????ydist?+=?cent->y?-?y;

????????????????????centc++;
????????????????}
????????????}????
????????}
????}?else?{
????????for?（y?=?0;?y?????????????for?（x?=?0;?x?????????????????if?（*（out++））?{
????????????????????if?（x?>?cent->x）
????????????????????????xdist?+=?x?-?cent->x;
????????????????????else?if?（x?x）
????????????????????????xdist?+=?cent->x?-?x;

????????????????????if?（y?>?cent->y）
????????????????????????ydist?+=?y?-?cent->y;
????????????????????else?if?（y?y）
????????????????????????ydist?+=?cent->y?-?y;

????????????????????centc++;
????????????????}
????????????}????
????????}
????