-#include<stdio.h>
-#include<stdarg.h>
-#include<time.h>
-#include<math.h>
-
-#include "tdb.h"
-#include "s.C"
-
-#define SWIDTH 80
-#define SHEIGHT 45
-
-double
-mean(uint8_t *dat, int n)
-{
- int s = 0;
- for( int i=0; i<n; ++i, ++dat ) s += *dat;
- return (double)s / n;
-}
-
-double
-variance(uint8_t *dat, double m, int n)
-{
- double ss = 0, s = 0;
- for( int i=0; i<n; ++i, ++dat ) {
- double dx = *dat-m;
- s += dx; ss += dx*dx;
- }
- return (ss - s*s / n) / (n-1);
-}
-
-double
-std_dev(uint8_t *dat, double m, int n)
-{
- return sqrt(variance(dat,m,n));
-}
-
-double
-center(uint8_t *dat, int n, int stride)
-{
- int s = 0, ss = 0;
- for( int i=0; i<n; dat+=stride ) { s += *dat; ss += ++i * *dat; }
- return s > 0 ? (double)ss / s : n / 2.;
-}
-
-void
-centroid(uint8_t *dat, int w, int h, double *xx, double *yy)
-{
- double x = 0, y = 0;
- uint8_t *dp = dat;
- for( int i=h; --i>=0; dp+=w ) x += center(dp, w, 1);
- for( int i=w; --i>=0; ) y += center(dat+i, h, w);
- *xx = x / h; *yy = y / w;
-}
-
-
-int main(int ac, char **av)
-{
- int ret; setbuf(stdout,0);
- theDb db;
- db.open(av[1]);
- //db.access(av[1], 34543, 0);
- if( !db.opened() || db.error() ) exit(1);
-
- FILE *fp = fopen(av[2],"r");
- for( int i=3; --i>=0; ) {
- for( int ch=fgetc(fp); ch>=0 && ch!='\n'; ch=fgetc(fp) );
- }
-
- int w = SWIDTH, h = SHEIGHT, sfrm_sz = w * h;
- uint8_t dat[sfrm_sz];
- fread(dat,1,sfrm_sz,fp);
-
- double mn = mean(dat,sfrm_sz);
- double sd = std_dev(dat,mn,sfrm_sz);
- double cx, cy; centroid(dat, w, h, &cx, &cy);
- double moment = cx + cy;
- if( (ret = Video_frameLoc::ikey_Frame_weight(db.video_frame, mn).Locate()) != 0 ) {
- printf(" not found, ret = %d\n",ret);
- return 1;
- }
- printf(" id %d, mean %f-%f=%f, std_dev %f-%f=%f, "
- " cx %f-%f=%f, cy %f-%f=%f, moment %f-%f=%f\n", db.video_frame.id(),
- mn, db.video_frame.Frame_mean(), mn-db.video_frame.Frame_mean(),
- sd, db.video_frame.Frame_std_dev(), sd-db.video_frame.Frame_std_dev(),
- cx, db.video_frame.Frame_cx(), cx-db.video_frame.Frame_cx(),
- cy, db.video_frame.Frame_cy(), cy-db.video_frame.Frame_cy(),
- moment, db.video_frame.Frame_moment(), moment-db.video_frame.Frame_moment());
- if( (ret = Video_frameLoc::ikey_Frame_center(db.video_frame, moment).Locate()) != 0 ) {
- printf(" not found, ret = %d\n",ret);
- return 1;
- }
- printf(" id %d, mean %f-%f=%f, std_dev %f-%f=%f, "
- " cx %f-%f=%f, cy %f-%f=%f, moment %f-%f=%f\n", db.video_frame.id(),
- mn, db.video_frame.Frame_mean(), mn-db.video_frame.Frame_mean(),
- sd, db.video_frame.Frame_std_dev(), sd-db.video_frame.Frame_std_dev(),
- cx, db.video_frame.Frame_cx(), cx-db.video_frame.Frame_cx(),
- cy, db.video_frame.Frame_cy(), cy-db.video_frame.Frame_cy(),
- moment, db.video_frame.Frame_moment(), moment-db.video_frame.Frame_moment());
- db.close();
- return 0;
-}
-