X-Git-Url: http://mj.ucw.cz/gitweb/?a=blobdiff_plain;ds=sidebyside;f=images%2Fsig-cmp.c;h=c663c2b3edc56afd1c0db7afb9606b3b7538d3dc;hb=c4bf633211b0424492b5a3937d6a6d2e0d79a4cf;hp=fef29407f2e871ad1c8b5beb6a6acc7c1a0695f1;hpb=5c62f53d3cc21d624a99214fc21643a80d7c1696;p=libucw.git

diff --git a/images/sig-cmp.c b/images/sig-cmp.c
index fef29407..c663c2b3 100644
--- a/images/sig-cmp.c
+++ b/images/sig-cmp.c
@@ -1,5 +1,5 @@
 /*
- *	Image Library -- Comparitions of image signatures
+ *	Image Library -- Comparisions of image signatures
  *
  *	(c) 2006 Pavel Charvat <pchar@ucw.cz>
  *
@@ -10,108 +10,12 @@
 #undef LOCAL_DEBUG
 
 #include "lib/lib.h"
-#include "lib/math.h"
 #include "images/math.h"
 #include "images/images.h"
 #include "images/signature.h"
 
 #include <stdio.h>
 
-static uns
-image_signatures_dist_fuzzy(struct image_signature *sig1, struct image_signature *sig2)
-{
-  DBG("image_signatures_dist_1()");
-
-  uns cnt1 = sig1->len;
-  uns cnt2 = sig2->len;
-  struct image_region *reg1 = sig1->reg;
-  struct image_region *reg2 = sig2->reg;
-  uns mf[IMAGE_REG_MAX][IMAGE_REG_MAX], mh[IMAGE_REG_MAX][IMAGE_REG_MAX];
-  uns lf[IMAGE_REG_MAX * 2], lh[IMAGE_REG_MAX * 2];
-  uns df = sig1->df + sig2->df, dh = sig1->dh + sig2->dh;
-
-  /* Compute distance matrix */
-  for (uns i = 0; i < cnt1; i++)
-    for (uns j = 0; j < cnt2; j++)
-      {
-	uns d = 0;
-	for (uns k = 0; k < IMAGE_REG_F; k++)
-	  {
-	    int dif = reg1[i].f[k] - reg2[j].f[k];
-	    d += dif * dif;
-	  }
-	mf[i][j] = d;
-	d = 0;
-	for (uns k = 0; k < IMAGE_REG_H; k++)
-	  {
-	    int dif = reg1[i].h[k] - reg2[j].h[k];
-	    d += dif * dif;
-	  }
-	mh[i][j] = d;
-      }
-
-  uns lfs = 0, lhs = 0;
-  for (uns i = 0; i < cnt1; i++)
-    {
-      uns f = mf[i][0], h = mh[i][0];
-      for (uns j = 1; j < cnt2; j++)
-        {
-	  f = MIN(f, mf[i][j]);
-	  h = MIN(h, mh[i][j]);
-	}
-      lf[i] = (df * 0x10000) / (df + fast_sqrt_u32(f));
-      lh[i] = (dh * 0x10000) / (dh + fast_sqrt_u32(h));
-      lfs += lf[i] * (6 * reg1[i].wa + 2 * reg1[i].wb);
-      lhs += lh[i] * reg1[i].wa;
-    }
-  for (uns i = 0; i < cnt2; i++)
-    {
-      uns f = mf[0][i], h = mh[0][i];
-      for (uns j = 1; j < cnt1; j++)
-        {
-	  f = MIN(f, mf[j][i]);
-	  h = MIN(h, mh[j][i]);
-	}
-      lf[i + cnt1] = (df * 0x10000) / (df + fast_sqrt_u32(f));
-      lh[i + cnt1] = (dh * 0x10000) / (dh + fast_sqrt_u32(h));
-      lfs += lf[i] * (6 * reg2[i].wa + 2 * reg2[i].wb);
-      lhs += lh[i] * reg2[i].wa;
-    }
-
-  uns measure = lfs * 6 + lhs * 2 * 8;
-
-#ifdef LOCAL_DEBUG
-  /* Display similarity vectors */
-  byte buf[2 * IMAGE_REG_MAX * 16 + 3], *b = buf;
-  for (uns i = 0; i < cnt1 + cnt2; i++)
-    {
-      if (i)
-	*b++ = ' ';
-      if (i == cnt1)
-	*b++ = '~', *b++ = ' ';
-      b += sprintf(b, "%.4f", (double)lf[i] / 0x10000);
-    }
-  *b = 0;
-  DBG("Lf=(%s)", buf);
-  b = buf;
-  for (uns i = 0; i < cnt1 + cnt2; i++)
-    {
-      if (i)
-	*b++ = ' ';
-      if (i == cnt1)
-	*b++ = '~', *b++ = ' ';
-      b += sprintf(b, "%.4f", (double)lh[i] / 0x10000);
-    }
-  *b = 0;
-  DBG("Lh=(%s)", buf);
-  DBG("Lfm=%.4f", lfs / (double)(1 << (3 + 8 + 16)));
-  DBG("Lhm=%.4f", lhs / (double)(1 << (8 + 16)));
-  DBG("measure=%.4f", measure / (double)(1 << (3 + 3 + 8 + 16)));
-#endif
-
-  return (1 << (3 + 3 + 8 + 16)) - measure;
-}
-
 #define ASORT_PREFIX(x) image_signatures_dist_integrated_##x
 #define ASORT_KEY_TYPE uns
 #define ASORT_ELT(i) items[i]
@@ -121,32 +25,3 @@ image_signatures_dist_fuzzy(struct image_signature *sig1, struct image_signature
 #define EXPLAIN
 #include "images/sig-cmp-gen.h"
 #include "images/sig-cmp-gen.h"
-
-uns
-image_signatures_dist(struct image_signature *sig1, struct image_signature *sig2)
-{
-  switch (image_sig_compare_method)
-    {
-      case 0:
-	return image_signatures_dist_integrated(sig1, sig2);
-      case 1:
-	return image_signatures_dist_fuzzy(sig1, sig2);
-      default:
-	ASSERT(0);
-    }
-}
-
-uns
-image_signatures_dist_explain(struct image_signature *sig1, struct image_signature *sig2, void (*msg)(byte *text, void *param), void *param)
-{
-  switch (image_sig_compare_method)
-    {
-      case 0:
-	return image_signatures_dist_integrated_explain(sig1, sig2, msg, param);
-      case 1:
-	return image_signatures_dist_fuzzy(sig1, sig2);
-      default:
-	ASSERT(0);
-    }
-}
-