Merge pull request #25672 from mshabunin:fix-approx

asmorkalov · web-flow · commit 2624929ec66a · 2024-05-30T14:13:26.000+03:00
imgproc: fix contour approximation, added test
diff --git a/modules/imgproc/src/contours_approx.cpp b/modules/imgproc/src/contours_approx.cpp
@@ -231,7 +231,7 @@ static void pass_cleanup(vector<ApproxItem>& ares, size_t start_idx)
 
     const size_t len = ares.size();
     size_t first = start_idx;
-    for (size_t i = start_idx, prev = i; i < len; ++i)
+    for (size_t i = start_idx, prev = start_idx; i < len; ++i)
     {
         ApproxItem& item = ares[i];
         if (item.removed)
@@ -248,18 +248,17 @@ static void pass_cleanup(vector<ApproxItem>& ares, size_t start_idx)
 
                     if (s1 > s2 || (s1 == s2 && ares[prev].k <= ares[i].k))
                         /* remove second */
-                        clear_until(ares, get_next_idx(ares, prev), get_next_idx(ares, i));
+                        clear_until(ares, prev, get_next_idx(ares, i));
                     else
                         /* remove first */
-                        clear_until(ares, first, i);
+                        clear_until(ares, get_next_idx(ares, first), i);
                 }
                 else
                 {
                     first = get_next_idx(ares, first);
                     clear_until(ares, first, i);
                 }
             }
-            clear_until(ares, first, i);
             first = i;
             count = 1;
         }
diff --git a/modules/imgproc/test/test_contours_new.cpp b/modules/imgproc/test/test_contours_new.cpp
@@ -602,4 +602,23 @@ TEST(Imgproc_FindContours, link_runs)
 #endif
 }
 
+TEST(Imgproc_FindContours, regression_25663)
+{
+    uint8_t data[]={0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1};
+    Mat1b img(32,32,data);
+    vector<vector<Point>> contours;
+    vector<Vec4i> hierarchy;
+    findContours(img, contours, RETR_EXTERNAL, CHAIN_APPROX_TC89_L1);
+#if CHECK_OLD
+    vector<vector<Point>> contours_o;
+    findContours_legacy(img, contours_o, RETR_EXTERNAL, CHAIN_APPROX_TC89_L1);
+    ASSERT_EQ(contours_o.size(), contours.size());
+    for (size_t i = 0; i < contours_o.size(); ++i)
+    {
+        SCOPED_TRACE(format("contour = %zu", i));
+        EXPECT_MAT_NEAR(Mat(contours_o[i]), Mat(contours[i]), 0);
+    }
+#endif
+}
+
 }}  // namespace opencv_test

Original file line number	Diff line number	Diff line change
`@@ -231,7 +231,7 @@ static void pass_cleanup(vector<ApproxItem>& ares, size_t start_idx)`
`231`	`231`
`232`	`232`	`const size_t len = ares.size();`
`233`	`233`	`size_t first = start_idx;`
`234`		`- for (size_t i = start_idx, prev = i; i < len; ++i)`
	`234`	`+ for (size_t i = start_idx, prev = start_idx; i < len; ++i)`
`235`	`235`	`{`
`236`	`236`	`ApproxItem& item = ares[i];`
`237`	`237`	`if (item.removed)`
`@@ -248,18 +248,17 @@ static void pass_cleanup(vector<ApproxItem>& ares, size_t start_idx)`
`248`	`248`
`249`	`249`	`if (s1 > s2 \|\| (s1 == s2 && ares[prev].k <= ares[i].k))`
`250`	`250`	`/* remove second */`
`251`		`- clear_until(ares, get_next_idx(ares, prev), get_next_idx(ares, i));`
	`251`	`+ clear_until(ares, prev, get_next_idx(ares, i));`
`252`	`252`	`else`
`253`	`253`	`/* remove first */`
`254`		`- clear_until(ares, first, i);`
	`254`	`+ clear_until(ares, get_next_idx(ares, first), i);`
`255`	`255`	`}`
`256`	`256`	`else`
`257`	`257`	`{`
`258`	`258`	`first = get_next_idx(ares, first);`
`259`	`259`	`clear_until(ares, first, i);`
`260`	`260`	`}`
`261`	`261`	`}`
`262`		`- clear_until(ares, first, i);`
`263`	`262`	`first = i;`
`264`	`263`	`count = 1;`
`265`	`264`	`}`