From 5f082b9876c3bd1e604c5961fdb29daab2efb0a8 Mon Sep 17 00:00:00 2001
From: "marina.kolpakova" <marina.kolpakova@itseez.com>
Date: Wed, 23 Jan 2013 22:21:56 +0400
Subject: [PATCH] improve ROC test script: handle ignored

---
 apps/sft/misc/roc_test.py |  6 ++++--
 apps/sft/misc/sft.py      | 36 +++++++++++++++++++++++++++++++-----
 2 files changed, 35 insertions(+), 7 deletions(-)

diff --git a/apps/sft/misc/roc_test.py b/apps/sft/misc/roc_test.py
index 8682119..bc436af 100755
--- a/apps/sft/misc/roc_test.py
+++ b/apps/sft/misc/roc_test.py
@@ -52,6 +52,7 @@ if __name__ == "__main__":
 
         confidenses = []
         tp          = []
+        ignored     = []
 
         while True:
             ret, img = camera.read()
@@ -77,12 +78,13 @@ if __name__ == "__main__":
             confs.sort(lambda x, y : -1  if (x - y) > 0 else 1)
             confidenses = confidenses + confs
 
-            matched = sft.match(boxes, dts)
+            matched, skip_list = sft.match(boxes, dts)
             tp = tp + matched
+            ignored = ignored + skip_list
 
             print nframes, nannotated
 
-        fppi, miss_rate = sft.computeROC(confidenses, tp, nannotated, nframes)
+        fppi, miss_rate = sft.computeROC(confidenses, tp, nannotated, nframes, ignored)
         sft.plotLogLog(fppi, miss_rate, plot_colors[idx])
 
     sft.showPlot("roc_curve.png")
\ No newline at end of file
diff --git a/apps/sft/misc/sft.py b/apps/sft/misc/sft.py
index 15f10fb..a732c77 100644
--- a/apps/sft/misc/sft.py
+++ b/apps/sft/misc/sft.py
@@ -3,6 +3,7 @@
 import cv2, re, glob
 import numpy             as np
 import matplotlib.pyplot as plt
+from itertools import izip
 
 """ Convert numPy matrices with rectangles and confidences to sorted list of detections."""
 def convert2detections(rects, confs, crop_factor = 0.125):
@@ -39,10 +40,12 @@ def cumsum(n):
     return cum
 
 """ Compute x and y arrays for ROC plot"""
-def computeROC(confidenses, tp, nannotated, nframes):
-    confidenses, tp = zip(*sorted(zip(confidenses, tp), reverse = True))
+def computeROC(confidenses, tp, nannotated, nframes, ignored):
+    confidenses, tp, ignored = zip(*sorted(zip(confidenses, tp, ignored), reverse = True))
 
     fp = [(1 - x) for x in tp]
+    fp = [(x - y) for x, y in izip(fp, ignored)]
+
     fp = cumsum(fp)
     tp = cumsum(tp)
     miss_rate = [(1 - x / (nannotated + 0.000001)) for x in tp]
@@ -81,16 +84,27 @@ def match(gts, dts):
     # Cartesian product for each detection BB_dt with each BB_gt
     overlaps = [[dt.overlap(gt) for gt in gts]for dt in dts]
 
-    matches_gt = [0]*len(gts)
-    matches_dt = [0]*len(dts)
+    matches_gt     = [0]*len(gts)
+    matches_dt     = [0]*len(dts)
+    matches_ignore = [0]*len(dts)
 
     for idx, row in enumerate(overlaps):
         imax = row.index(max(row))
 
+        # try to match ground thrush
         if (matches_gt[imax] == 0 and row[imax] > 0.5):
             matches_gt[imax] = 1
             matches_dt[idx]  = 1
-    return matches_dt
+
+    for idx, dt in enumerate(dts):
+        # try to math ignored
+        if matches_dt[idx] == 0:
+            row = gts
+            row = [i for i in row if (i[3] - i[1]) < 53 or (i[3] - i[1]) >  256]
+            for each in row:
+                if dts[idx].overlapIgnored(each) > 0.5:
+                    matches_ignore[idx] = 1
+    return matches_dt, matches_ignore
 
 
 def plotLogLog(fppi, miss_rate, c):
@@ -136,6 +150,18 @@ class Detection:
 
         return cross_area / union_area
 
+        # we use rect-style for dt and box style for gt. ToDo: fix it
+    def overlapIgnored(self, b):
+
+        a = self.bb
+        w = min( a[0] + a[2], b[2]) - max(a[0], b[0]);
+        h = min( a[1] + a[3], b[3]) - max(a[1], b[1]);
+
+        cross_area = 0.0 if (w < 0 or h < 0) else float(w * h)
+        self_area = (a[2] * a[3]);
+
+        return cross_area / self_area
+
     def mark_matched(self):
         self.matched = True
 
-- 
2.7.4