cinelerra-5.0/cinelerra/mediadb.C

   1 #include <stdio.h>
   2 #include <stdarg.h>
   3 #include <stdlib.h>
   4 #include <unistd.h>
   5 #include <string.h>
   6 #include <limits.h>
   7 #include <time.h>
   8
   9 #include "../db/tdb.h"
  10 #include "../db/s.C"
  11
  12 #include "linklist.h"
  13 #include "mediadb.h"
  14 #include "filexml.h"
  15
  16
  17
  18 void write_pgm(uint8_t *tp, int w, int h, const char *fmt, ...)
  19 {
  20   va_list ap;    va_start(ap, fmt);
  21   char fn[256];  vsnprintf(fn, sizeof(fn), fmt, ap);
  22   va_end(ap);
  23   FILE *fp = !strcmp(fn,"-") ? stdout : fopen(fn,"w");
  24   if( fp ) {
  25     fprintf(fp,"P5\n%d %d\n255\n",w,h);
  26     fwrite(tp,w,h,fp);
  27     fclose(fp);
  28   }
  29 }
  30
  31 void write_ppm(uint8_t *tp, int w, int h, const char *fmt, ...)
  32 {
  33   va_list ap;    va_start(ap, fmt);
  34   char fn[256];  vsnprintf(fn, sizeof(fn), fmt, ap);
  35   va_end(ap);
  36   FILE *fp = !strcmp(fn,"-") ? stdout : fopen(fn,"w");
  37   if( fp ) {
  38     fprintf(fp,"P6\n%d %d\n255\n",w,h);
  39     fwrite(tp,3*w,h,fp);
  40     fclose(fp);
  41   }
  42 }
  43
  44 static inline int clip(int v, int mn, int mx)
  45 {
  46   return v<mn ? mn : v>mx ? mx : v;
  47 }
  48
  49 void pbm_diff(uint8_t *ap, uint8_t *bp, int w, int h, const char *fmt, ...)
  50 {
  51   va_list va;    va_start(va, fmt);
  52   char fn[256];  vsnprintf(fn, sizeof(fn), fmt, va);
  53   va_end(va);
  54   FILE *fp = !strcmp(fn,"-") ? stdout : fopen(fn,"w");
  55   if( fp ) {
  56     int sz = w*h;
  57     uint8_t dat[sz], *dp = dat;
  58     for( int i=sz; --i>=0; ++dp,++ap,++bp ) *dp = clip(*ap-*bp+128, 0 ,255);
  59     fprintf(fp,"P5\n%d %d\n255\n",w,h);
  60     fwrite(dat,w,h,fp);
  61     fclose(fp);
  62   }
  63 }
  64
  65
  66 void Scale::
  67 scale(uint8_t *odata, int ow, int oh, int sx, int sy, int sw, int sh)
  68 {
  69   pw = (double)(ow * oh) / (rw * rh);
  70   odat = odata + sy*ow + sx;
  71   iy1 = 0;  yf1 = 0;
  72   double r = 0.5;
  73
  74   for( int dy=1; dy<=sh; ++dy, odat+=ow ) {
  75     iy0 = iy1;  yf0 = 1.0 - yf1;
  76     double ny = (double)(dy * rh) / oh;
  77     iy1 = (int) ny;  yf1 = ny-iy1;
  78     uint8_t *bp = odat;
  79     ix1 = 0;  xf1 = 0;
  80
  81     for( int dx=1; dx<=sw; ++dx ) {
  82       ix0 = ix1;  xf0 = 1.0 - xf1;
  83       double nx = (double)(dx*rw) / ow;
  84       ix1 = (int)nx;  xf1 = nx-ix1;
  85       double px = scalexy();
  86       int ipx = px + r;
  87       r += px - ipx;
  88       *bp++ = ipx;
  89     }
  90   }
  91 }
  92
  93
  94 Deletions::
  95 Deletions(int pid, const char *fn)
  96 {
  97         this->pid = pid;
  98         strcpy(this->filepath,fn);
  99 }
 100
 101 Deletions::
 102 ~Deletions()
 103 {
 104 }
 105
 106
 107 int Deletions::
 108 load(FileXML &xml)
 109 {
 110         for(;;) {
 111                 if( xml.read_tag() != 0 ) return 1;
 112                 if( !xml.tag.title_is("DEL") ) break;
 113                 double time = xml.tag.get_property("TIME", (double)0.0);
 114                 int action = xml.tag.get_property("ACTION", (int)0);
 115                 append(new Dele(time, action));
 116         }
 117         return 0;
 118 }
 119
 120 Deletions *Deletions::
 121 read_dels(const char *filename)
 122 {
 123         FileXML xml;
 124         if( xml.read_from_file(filename, 1) ) return 0;
 125         do {
 126                 if( xml.read_tag() ) return 0;
 127         } while( !xml.tag.title_is("DELS") );
 128
 129         int pid = xml.tag.get_property("PID", (int)-1);
 130         const char *file = xml.tag.get_property("FILE");
 131         Deletions *dels = new Deletions(pid, file);
 132         if( dels->load(xml) || !xml.tag.title_is("/DELS") ) {
 133                  delete dels; dels = 0;
 134         }
 135         return dels;
 136 }
 137
 138 int Deletions::
 139 save(FileXML &xml)
 140 {
 141         xml.tag.set_title("DELS");
 142         xml.tag.set_property("PID", pid);
 143         xml.tag.set_property("FILE", filepath);
 144         xml.append_tag();
 145         xml.append_newline();
 146
 147         for( Dele *del=first; del; del=del->next ) {
 148                 xml.tag.set_title("DEL");
 149                 xml.tag.set_property("TIME", del->time);
 150                 xml.tag.set_property("ACTION", del->action);
 151                 xml.append_tag();
 152                 xml.append_newline();
 153         }
 154
 155         xml.tag.set_title("/DELS");
 156         xml.append_tag();
 157         xml.append_newline();
 158         return 0;
 159 }
 160
 161 int Deletions::
 162 write_dels(const char *filename)
 163 {
 164         FILE *fp = fopen(filename, "wb");
 165         if( !fp != 0 ) return 1;
 166         FileXML xml;
 167         save(xml);
 168         xml.terminate_string();
 169         xml.write_to_file(fp);
 170         fclose(fp);
 171         return 0;
 172 }
 173
 174
 175 Clip::
 176 Clip(int clip_id, double start, double end, int index)
 177 {
 178         this->clip_id = clip_id;
 179         this->start = start;
 180         this->end = end;
 181         this->index = index;
 182         this->votes = 0;
 183         this->groups = 0;
 184         this->cycle = 0;
 185         this->muted = 0;
 186         this->best = 1.e99;
 187 }
 188
 189 Clip::
 190 ~Clip()
 191 {
 192         Clips *clips = (Clips *)owner;
 193         if( clips && clips->current == this ) clips->current = 0;
 194 }
 195
 196 Clip *Clips::locate(int id, double start, double end)
 197 {
 198         if( current && current->clip_id == id &&
 199             end >= current->start && start < current->end )
 200                 return current;
 201         for( current=first; current!=0; current=current->next ) {
 202                 if( current->clip_id != id ) continue;
 203                 if( current->start <= end && current->end > start ) break;
 204         }
 205         return current;
 206 }
 207
 208 Clips::
 209 Clips(int pid)
 210 {
 211         this->pid = pid;
 212         this->current = 0;
 213         this->count = 0;
 214         this->cycle = 0;
 215 }
 216
 217 Clips::
 218 ~Clips()
 219 {
 220 }
 221
 222 Clip *Clips::
 223 new_clip(int clip_id, double start, double end, int index)
 224 {
 225         return append(new Clip(clip_id, start, end, index));
 226 }
 227
 228 int Clips::
 229 get_clip(int idx, int &id, double &start, double &end)
 230 {
 231         int votes = -1, cid = -1;
 232         double st = 0, et = 0;
 233
 234         for( Clip *clip=first; clip!=0; clip=clip->next ) {
 235                 if( clip->index != idx ) continue;
 236                 // must have prefix+suffix frames
 237                 if( clip->groups != 3 ) continue;
 238                 if( clip->votes < votes ) continue;
 239                 if( clip->votes > votes ) {
 240                         // first/better match
 241                         votes = clip->votes;
 242                         cid = clip->clip_id;
 243                         st = clip->start;
 244                         et = clip->end;
 245                         continue;
 246                 }
 247                 // reject if not overlapping
 248                 if( clip->start > et ) continue;
 249                 if( clip->end < st ) continue;
 250                 // expand range on overlap
 251                 if( clip->start < st ) st = clip->start;
 252                 if( clip->end > et ) et = clip->end;
 253         }
 254
 255         if( votes > 2 ) {
 256                 start = st;  end = et;  id = cid;
 257                 return 1;
 258         }
 259
 260         return 0;
 261 }
 262
 263 int Clips::
 264 check(int clip_id, double start, double end, int group, double err)
 265 {
 266         Clip *clip = locate(clip_id, start, end);
 267         if( clip ) {
 268                 if( err < clip->best ) {
 269                         clip->best = err;
 270                         if( start < clip->start ) clip->start = start;
 271                         if( end > clip->end ) clip->end = end;
 272                 }
 273         }
 274         else
 275                 clip = new_clip(clip_id, start, end, count++);
 276         clip->groups |= group;
 277         if( cycle != clip->cycle ) {
 278                  clip->cycle = cycle;
 279                 ++clip->votes;
 280 //printf("clip %d votes %d\n",clip_id,clip->votes);
 281         }
 282         return 0;
 283 }
 284
 285 int Clips::
 286 save(FileXML &xml)
 287 {
 288         xml.tag.set_title("CLIPS");
 289         xml.tag.set_property("PID", pid);
 290         xml.append_tag();
 291         xml.append_newline();
 292
 293         for( int idx=0; idx<count; ++idx ) {
 294                 int id;  double start, end;
 295                 if( !get_clip(idx, id, start, end) ) continue;
 296                 xml.tag.set_title("CLIP");
 297                 xml.tag.set_property("ID", id);
 298                 xml.tag.set_property("START", start);
 299                 xml.tag.set_property("END", end);
 300                 xml.append_tag();
 301                 xml.append_newline();
 302         }
 303
 304         xml.tag.set_title("/CLIPS");
 305         xml.append_tag();
 306         xml.append_newline();
 307         return 0;
 308 }
 309
 310 int Clips::
 311 load(FileXML &xml)
 312 {
 313         int i = 0;
 314         for(;;) {
 315                 if( xml.read_tag() != 0 ) return 1;
 316                 if( !xml.tag.title_is("CLIP") ) break;
 317                 int clip_id = xml.tag.get_property("ID", (int)0);
 318                 double start = xml.tag.get_property("START", (double)0.0);
 319                 double end = xml.tag.get_property("END", (double)0.0);
 320                 append(new Clip(clip_id, start, end, i++));
 321         }
 322         count = i;
 323         return 0;
 324 }
 325
 326
 327 Snip::
 328 Snip(int clip_id, double start, double end, int index)
 329  : Clip(clip_id, start, end, index)
 330 {
 331 }
 332
 333 Snip::
 334 ~Snip()
 335 {
 336 }
 337
 338 Snips::
 339 Snips()
 340 {
 341 }
 342
 343 Snips::
 344 ~Snips()
 345 {
 346 }
 347
 348 Clip *Snips::
 349 new_clip(int clip_id, double start, double end, int index)
 350 {
 351         return append(new Snip(clip_id, start, end, index));
 352 }
 353
 354
 355 MediaDb::
 356 MediaDb()
 357 {
 358         opened = 0;
 359         db = new theDb();
 360 }
 361
 362 MediaDb::
 363 ~MediaDb()
 364 {
 365         if( opened ) closeDb();
 366         delete db;
 367 }
 368
 369 int MediaDb::is_open() { return opened > 0 ? opened : 0; }
 370 int MediaDb::transaction() { return db->transaction(); }
 371
 372 // clip_set accessors
 373 int MediaDb::clip_id(int id) { return db->clip_set.FindId(id); }
 374 int MediaDb::clip_id() { return db->clip_set.id(); }
 375 int MediaDb::clips_first_id() { return db->clip_set.FirstId(clip_id_loc); }
 376 int MediaDb::clips_next_id() { return db->clip_set.NextId(clip_id_loc); }
 377 const char *MediaDb::clip_title() { return db->clip_set.Title(); }
 378 const char *MediaDb::clip_path() { return db->clip_set.Asset_path(); }
 379 double MediaDb::clip_position() { return db->clip_set.Position(); }
 380 double MediaDb::clip_framerate() { return db->clip_set.Framerate(); }
 381 double MediaDb::clip_average_weight() { return db->clip_set.Average_weight(); }
 382 void MediaDb::clip_average_weight(double avg_wt) { return db->clip_set.Average_weight(avg_wt); }
 383 int MediaDb::clip_frames() { return db->clip_set.Frames(); }
 384 double MediaDb::clip_length() {
 385         return clip_framerate()>0 ? clip_frames() / clip_framerate() : 0;
 386 }
 387 int MediaDb::clip_prefix_size() { return db->clip_set.Prefix_size(); }
 388 int MediaDb::clip_suffix_size() { return db->clip_set.Suffix_size(); }
 389 int MediaDb::clip_size() { return clip_prefix_size()+clip_suffix_size(); }
 390 double *MediaDb::clip_weights() { return (double *)db->clip_set._Weights(); }
 391 int64_t MediaDb::clip_creation_time() { return db->clip_set.Creation_time(); }
 392 int64_t MediaDb::clip_system_time() { return db->clip_set.System_time(); }
 393
 394 // MediaDb::timeline accessors
 395 int MediaDb::timeline_id(int id) { return db->timeline.FindId(id); }
 396 int MediaDb::timeline_id() { return db->timeline.id(); }
 397 int MediaDb::timeline_clip_id() { return db->timeline.Clip_id(); }
 398 int MediaDb::timeline_sequence_no() { return db->timeline.Sequence_no(); }
 399 int MediaDb::timeline_frame_id() { return db->timeline.Frame_id(); }
 400 int MediaDb::timeline_group() { return db->timeline.Group(); }
 401 double MediaDb::timeline_offset() { return db->timeline.Time_offset(); }
 402
 403 // MediaDb::clip_view accessors
 404 int MediaDb::views_id(int id) { return db->clip_views.FindId(id); }
 405 int MediaDb::views_id() { return db->clip_views.id(); }
 406 int MediaDb::views_clip_id(int id) {
 407         return Clip_viewsLoc::ikey_Clip_access(db->clip_views,id).Find();
 408 }
 409 int MediaDb::views_clip_id() { return db->clip_views.Access_clip_id(); }
 410 int64_t MediaDb::views_access_time() { return db->clip_views.Access_time(); }
 411 int MediaDb::views_access_count() { return db->clip_views.Access_count(); }
 412
 413
 414 double MediaDb::
 415 mean(uint8_t *dat, int n, int stride)
 416 {
 417   int s = 0;
 418   for( int i=0; i<n; ++i, dat+=stride ) s += *dat;
 419   return (double)s / n;
 420 }
 421
 422 double MediaDb::
 423 variance(uint8_t *dat, double m, int n, int stride)
 424 {
 425   double ss = 0, s = 0;
 426   for( int i=0; i<n; ++i, dat+=stride ) {
 427     double dx = *dat-m;
 428     s += dx;  ss += dx*dx;
 429   }
 430   return (ss - s*s / n) / (n-1);
 431 }
 432
 433 double MediaDb::
 434 centroid(uint8_t *dat, int n, int stride)
 435 {
 436   int s = 0, ss = 0;
 437   for( int i=0; i<n; dat+=stride ) { s += *dat;  ss += ++i * *dat; }
 438   return s > 0 ? (double)ss / s : n / 2.;
 439 }
 440
 441 void MediaDb::
 442 centroid(uint8_t *dat, int w, int h, double &xx, double &yy)
 443 {
 444         double x = 0, y = 0;
 445         uint8_t *dp = dat;
 446         for( int i=h; --i>=0; dp+=w ) x += centroid(dp, w);
 447         for( int i=w; --i>=0; ) y += centroid(dat+i, h, w);
 448         xx = x / h;  yy = y / w;
 449 }
 450
 451 void MediaDb::
 452 deviation(uint8_t *a, uint8_t *b, int sz, int margin,
 453                 double &dev, int &ofs, int stride)
 454 {
 455   double best = 1e100;  int ii = -1;
 456   for( int i=-margin; i<=margin; ++i ) {
 457     int aofs = i < 0 ? 0 : i*stride;
 458     int bofs = i < 0 ? -i*stride : 0;
 459     uint8_t *ap = a + aofs, *bp = b + bofs;
 460     int ierr = 0; int n = sz - abs(i);
 461     for( int j=n; --j>=0; ap+=stride,bp+=stride ) ierr += abs(*ap - *bp);
 462     double err = (double)ierr / n;
 463     if( err < best ) { best = err;  ii = i; }
 464   }
 465   dev = best;
 466   ofs = ii;
 467 }
 468
 469 int64_t MediaDb::
 470 cmpr(uint8_t *ap, uint8_t *bp, int sz)
 471 {
 472   int64_t v = 0;
 473   for( int i=sz; --i>=0; ++ap,++bp ) v += abs(*ap-*bp);
 474   return v;
 475 }
 476
 477 int64_t MediaDb::
 478 diff(uint8_t *a, uint8_t *b, int w, int h, int dx, int dy, int bias)
 479 {
 480   int axofs = dx < 0 ? 0 : dx;
 481   int ayofs = w * (dy < 0 ? 0 : dy);
 482   int aofs = ayofs + axofs;
 483   int bxofs = dx < 0 ? -dx : 0;
 484   int byofs = w * (dy < 0 ? -dy : 0);
 485   int bofs = byofs + bxofs;
 486   uint8_t *ap = a + aofs, *bp = b + bofs;
 487   int ierr = 0, ww = w-abs(dx), hh = h-abs(dy);
 488   for( int y=hh; --y>=0; ap+=w, bp+=w ) {
 489     a = ap;  b = bp;
 490     for( int x=ww; --x>=0; ++a, ++b ) { ierr += abs(*a-bias - *b); }
 491   }
 492   return ierr;
 493 }
 494
 495 void MediaDb::
 496 fit(int *dat, int n, double &a, double &b)
 497 {
 498         double st2 = 0, sb = 0;
 499         int64_t sy = 0;
 500         for( int i=0; i<n; ++i ) sy += dat[i];
 501         double mx = (n-1)/2., my = (double)sy / n;
 502         for( int i=0; i<n; ++i ) {
 503                 double t = i - mx;
 504                 st2 += t * t;
 505                 sb += t * dat[i];
 506         }
 507         b = sb / st2;
 508         a = my - mx*b;
 509 }
 510
 511 void MediaDb::
 512 eq(uint8_t *sp, uint8_t *dp, int len)
 513 {
 514 #if 0
 515         int hist[256], wts[256], map[256];
 516         for( int i=0; i<256; ++i ) hist[i] = 0;
 517         uint8_t *bp = sp;
 518         for( int i=len; --i>=0; ++bp ) ++hist[*bp];
 519         int t = 0;
 520         for( int i=0; i<256; ++i ) { t += hist[i];  wts[i] = t; }
 521         int lmn = len/20, lmx = len-lmn;
 522         int mn =   0;  while( mn < 256 && wts[mn] < lmn ) ++mn;
 523         int mx = 255;  while( mx > mn && wts[mx] > lmx ) --mx;
 524         double a, b;
 525         fit(&wts[mn], mx-mn, a, b);
 526         double r = 256./len;
 527         double a1 = (a - b*mn) * r, b1 = b * r;
 528         for( int i=0 ; i<256;  ++i ) map[i] = clip(a1 + i*b1, 0, 255);
 529         for( int i=len; --i>=0; ++sp,++dp ) *dp = map[*sp];
 530 #else
 531         if( sp != dp ) memcpy(dp,sp,len);
 532 #endif
 533 }
 534
 535 // record constructors
 536 int MediaDb::
 537 new_frame(int clip_id, uint8_t *dat, int no, int group, double offset)
 538 {
 539 //printf("add %d at %f\n",no,tm);
 540         int fid, ret = get_media_key(dat, fid, MEDIA_FRAME_DIST, MEDIA_FRAME_ERRLMT);
 541         if( ret < 0 ) return -1;
 542         if( ret > 0 ) {
 543                 db->video_frame.Allocate();
 544                 uint8_t *fp = db->video_frame._Frame_data(SFRM_SZ);
 545                 eq(dat, fp, SFRM_SZ);
 546                 double mn = mean(dat, SFRM_SZ);
 547                 db->video_frame.Frame_mean(mn);
 548                 double sd = std_dev(dat, mn, SFRM_SZ);
 549                 db->video_frame.Frame_std_dev(sd);
 550                 double cx, cy;  centroid(dat, SWIDTH, SHEIGHT, cx, cy);
 551                 db->video_frame.Frame_cx(cx);
 552                 db->video_frame.Frame_cy(cy);
 553                 double moment = cx + cy;
 554                 db->video_frame.Frame_moment(moment);
 555                 db->video_frame.Construct();
 556                 fid = db->video_frame.id();
 557         }
 558         db->timeline.Allocate();
 559         db->timeline.Clip_id(clip_id);
 560         db->timeline.Sequence_no(no);
 561         db->timeline.Frame_id(fid);
 562         db->timeline.Group(group);
 563         db->timeline.Time_offset(offset);
 564         db->timeline.Construct();
 565         return 0;
 566 }
 567
 568 int MediaDb::
 569 new_clip_set(const char *title, const char *asset_path, double position,
 570         double framerate, int frames, int prefix_size, int suffix_size,
 571         int64_t creation_time, int64_t system_time)
 572 {
 573         db->clip_set.Allocate();
 574         db->clip_set.Title(title);
 575         db->clip_set.Asset_path(asset_path);
 576         db->clip_set.Position(position);
 577         db->clip_set.Framerate(framerate);
 578         db->clip_set.Frames(frames);
 579         db->clip_set.Prefix_size(prefix_size);
 580         db->clip_set.Suffix_size(suffix_size);
 581         db->clip_set._Weights(frames*sizeof(double));
 582         db->clip_set.Creation_time(creation_time);
 583         db->clip_set.System_time(system_time);
 584         db->clip_set.Construct();
 585
 586         return new_clip_view(db->clip_set.id(), creation_time);
 587 }
 588
 589 int MediaDb::
 590 new_clip_view(int clip_id, int64_t access_time)
 591 {
 592         db->clip_views.Allocate();
 593         db->clip_views.Access_clip_id(clip_id);
 594         db->clip_views.Access_time(access_time);
 595         db->clip_views.Access_count(1);
 596         db->clip_views.Construct();
 597         return 0;
 598 }
 599
 600 int MediaDb::
 601 access_clip(int clip_id, int64_t access_time)
 602 {
 603         if( Clip_viewsLoc::ikey_Clip_access(db->clip_views,clip_id).Find() ) return 1;
 604         db->clip_views.Destruct();
 605         if( !access_time ) { time_t at;  ::time(&at);  access_time = (int64_t)at; }
 606         db->clip_views.Access_time(access_time);
 607         int access_count = db->clip_views.Access_count();
 608         db->clip_views.Access_count(access_count+1);
 609         db->clip_views.Construct();
 610         return 0;
 611 }
 612
 613 int MediaDb::
 614 del_clip_set(int clip_id)
 615 {
 616         if( db->clip_set.FindId(clip_id) ) return 1;
 617         // delete clip_set
 618         db->clip_set.Destruct();
 619         db->clip_set.Deallocate();
 620         // delete clip_views
 621         if( Clip_viewsLoc::ikey_Clip_access(db->clip_views,clip_id).Find() ) return 1;
 622         db->clip_views.Destruct();
 623         db->clip_views.Deallocate();
 624         // delete timeline
 625         while( !TimelineLoc::ikey_Sequences(db->timeline,clip_id,0).Locate() ) {
 626                 if( clip_id != (int)db->timeline.Clip_id() ) break;
 627                 int frame_id = db->timeline.Frame_id();
 628                 db->timeline.Destruct();
 629                 db->timeline.Deallocate();
 630                 // check for more timeline refs
 631                 if( !TimelineLoc::ikey_Timelines(db->timeline, frame_id).Locate() &&
 632                         frame_id == (int)db->timeline.Frame_id() ) continue;
 633                 if( db->video_frame.FindId(frame_id) ) continue;
 634                 // delete frame
 635                 db->video_frame.Destruct();
 636                 db->video_frame.Deallocate();
 637         }
 638         return 0;
 639 }
 640
 641
 642 // db file operations
 643 int MediaDb::
 644 newDb()
 645 {
 646         if( opened > 0 || db->create(MEDIA_DB) ) {
 647                 printf("MediaDb::newDb failed\n");
 648                 return opened = -1;
 649         }
 650         return opened = 0;
 651 }
 652
 653 int MediaDb::
 654 openDb(int rw)
 655 {
 656         if( opened > 0 ) return 0;
 657         if( db->access(MEDIA_DB, MEDIA_SHM_KEY, rw) ) {
 658                 printf("MediaDb::openDb failed\n");
 659                 return opened = -1;
 660         }
 661         opened = 1;
 662         return 0;
 663 }
 664
 665 int MediaDb::
 666 resetDb(int rw)
 667 {
 668         int result = 0;
 669         if( access(MEDIA_DB, F_OK) ) result = newDb();
 670         if( !result ) result = openDb(rw);
 671         return result;
 672 }
 673
 674 void MediaDb::
 675 closeDb()
 676 {
 677         if( opened > 0 ) { detachDb();  db->close(); }
 678         opened = 0;
 679 }
 680
 681 void MediaDb::
 682 commitDb()
 683 {
 684         if( opened > 0 ) db->commit();
 685 }
 686
 687 void MediaDb::
 688 undoDb()
 689 {
 690         if( opened > 0 ) db->undo();
 691 }
 692
 693 int MediaDb::
 694 attachDb(int rw)
 695 {
 696         return opened > 0 ? db->attach(rw) : -1;
 697 }
 698
 699 int MediaDb::
 700 detachDb()
 701 {
 702         return opened > 0 ? db->detach() : -1;
 703 }
 704
 705
 706 int MediaDb::
 707 get_timelines(int frame_id)
 708 {
 709         return TimelineLoc::ikey_Timelines(db->timeline, frame_id).Locate();
 710 }
 711
 712 int MediaDb::
 713 next_timelines()
 714 {
 715         return TimelineLoc::rkey_Timelines(db->timeline).Next();
 716 }
 717
 718 int MediaDb::
 719 get_sequences(int clip_id, int seq_no)
 720 {
 721         int ret = TimelineLoc::ikey_Sequences(db->timeline,clip_id,seq_no).Locate();
 722         if( !ret && clip_id != (int)db->timeline.Clip_id() ) ret = 1;
 723         return ret;
 724 }
 725
 726 int MediaDb::
 727 next_sequences()
 728 {
 729         return TimelineLoc::rkey_Sequences(db->timeline).Next();
 730 }
 731
 732 // frame accessors
 733 int MediaDb::
 734 get_image(int id, uint8_t *dat, int &w, int &h)
 735 {
 736         int result = 1;
 737         if( !db->video_frame.FindId(id) ) {
 738                 memcpy(dat, db->video_frame.Frame_data(), SFRM_SZ);
 739                 w = SWIDTH;  h = SHEIGHT;
 740                 result = 0;
 741         }
 742         return result;
 743 }
 744
 745 int MediaDb::
 746 get_frame_key(uint8_t *dat, int &fid,
 747                 Clips *clips, double pos, int mask)
 748 {
 749         return get_media_key(dat, fid,
 750                         MEDIA_SEARCH_DIST, MEDIA_SEARCH_ERRLMT,
 751                         clips, pos, mask);
 752 }
 753
 754 int MediaDb::
 755 get_media_key(uint8_t *dat, int &fid, int dist, double errlmt,
 756                 Clips *clips, double pos, int mask)
 757 {
 758         fmean = mean(dat, SFRM_SZ);
 759         if( fmean < 17 ) return -1; // black, forbidden
 760         if( fmean > 239 ) return -1; // white, forbidden
 761         distance = dist;
 762         pixlmt = SFRM_SZ*errlmt;
 763         position = pos;
 764         group_mask = mask;
 765         fsd = std_dev(dat, fmean, SFRM_SZ);
 766         centroid(dat, SWIDTH, SHEIGHT, fcx, fcy);
 767         fmoment = fcx + fcy;
 768         eq(dat, frm, SFRM_SZ);
 769         if( clips ) ++clips->cycle;
 770         lid = nid = -1;
 771         lerr = LONG_MAX;
 772         int result = 1;
 773
 774         if( !Video_frameLoc::ikey_Frame_center(db->video_frame, fmoment).Locate() ) {
 775                 Video_frameLoc prev(db->video_frame), next(db->video_frame);
 776                 Video_frameLoc::rkey_Frame_center prev_ctr(prev), next_ctr(next);
 777                 if( !get_key(clips, prev_ctr, next_ctr) ) result = 0;
 778         }
 779         if( !Video_frameLoc::ikey_Frame_weight(db->video_frame, fmean).Locate() ) {
 780                 Video_frameLoc prev(db->video_frame), next(db->video_frame);
 781                 Video_frameLoc::rkey_Frame_weight prev_wt(prev), next_wt(next);
 782                 if( !get_key(clips, prev_wt, next_wt) ) result = 0;
 783         }
 784
 785         fid = lid;
 786         return result;
 787 }
 788
 789 int MediaDb::
 790 get_key(Clips *clips, Db::rKey &prev_rkey, Db::rKey &next_rkey)
 791 {
 792         key_compare(clips, next_rkey.loc);
 793         // search outward from search target for frame with least error
 794         Db::pgRef next_loc;  next_rkey.NextLoc(next_loc);
 795
 796         for( int i=1; i<=distance; ++i ) {
 797                 if( !prev_rkey.Locate(Db::keyLT) )
 798                         key_compare(clips,prev_rkey.loc);
 799                 if( !next_rkey.Next(next_loc) )
 800                         key_compare(clips,next_rkey.loc);
 801         }
 802
 803         return lerr < pixlmt ? 0 : 1;
 804 }
 805
 806 int64_t MediaDb::
 807 key_compare(Clips *clips, Db::ObjectLoc &loc)
 808 {
 809         int id = loc.id();
 810         if( nid == id ) return 1;
 811         nid = id;
 812         Video_frameLoc &frame = *(Video_frameLoc*)&loc;
 813         double dm = frame.Frame_mean()-fmean;
 814         if( fabs(dm) > MEDIA_MEAN_ERRLMT ) return 1;
 815         double ds = frame.Frame_std_dev()-fsd;
 816         if( fabs(ds) > MEDIA_STDDEV_ERRLMT ) return 1;
 817         double dx = frame.Frame_cx()-fcx;
 818         if( fabs(dx) > MEDIA_XCENTER_ERRLMT ) return 1;
 819         double dy = frame.Frame_cy()-fcy;
 820         if( fabs(dy) > MEDIA_YCENTER_ERRLMT ) return 1;
 821         uint8_t *dat = frame._Frame_data();
 822         //int64_t err = cmpr(frm, dat, SFRM_SZ);
 823         int64_t err = diff(frm, dat, SWIDTH, SHEIGHT, lround(dx), lround(dy), lround(dm));
 824         if( lerr > err ) { lerr = err;  lid = id; }
 825         if( clips && err < pixlmt ) add_timelines(clips, id, err);
 826         return 0;
 827 }
 828
 829 int MediaDb::
 830 add_timelines(Clips *clips, int fid, double err)
 831 {
 832         // found frame, find timelines
 833         if( get_timelines(fid) ) {
 834                 printf(_(" find timeline frame_id(%d) failed\n"), fid);
 835                 return 1;
 836         }
 837
 838         // process Timelines which contain Frame_id
 839         int cid = -1;
 840         for( int ret=0; !ret && timeline_frame_id()==fid; ret=next_timelines() ) {
 841                 int tid = timeline_clip_id();
 842                 if( tid != cid && clip_id(cid=tid) ) break;
 843                 double start_time = position - timeline_offset();
 844                 if( start_time < 0 ) continue;
 845                 int group = timeline_group();
 846                 if( (group & group_mask) == 0 ) continue;
 847                 double length = clip_length();
 848                 double end_time = start_time + length;
 849 //printf("clip=%d, pos=%f, time=%f-%f, %s\n",
 850 //  cid, position, start_time, end_time, clip_title());
 851                 clips->check(cid, start_time, end_time, group, err);
 852         }
 853
 854         return 0;
 855 }
 856