THESE/codes/meanshift/Ms_segmenter/Merge_regions.m

   1 function [New_inner_pts, New_boundary, New_shifts] = ...\r
   2         Merge_regions (Region1, Region2);\r
   3 \r
   4 Min_x = min(Region1.Shifts(1, 1), Region2.Shifts(1, 1));\r
   5 Max_x = max(Region1.Shifts(1, 2), Region2.Shifts(1, 2));\r
   6 Min_y = min(Region1.Shifts(2, 1), Region2.Shifts(2, 1));\r
   7 Max_y = max(Region1.Shifts(2, 2), Region2.Shifts(2, 2));\r
   8 \r
   9 Image1 = zeros(Max_x-Min_x+1, Max_y-Min_y+1);\r
  10 Image1(sub2ind(size(Image1), Region1.Boundary(1, :)-Min_x+1, ...\r
  11         Region1.Boundary(2, :)-Min_y+1)) = 1;\r
  12 To_start_fill = [Region1.Inner_pts(2, :)-Min_y+1; ...\r
  13         Region1.Inner_pts(1, :)-Min_x+1];\r
  14 Image1 = bwfill(Image1, To_start_fill(1, :), To_start_fill(2, :));\r
  15 \r
  16 Image2 = zeros(Max_x-Min_x+1, Max_y-Min_y+1);\r
  17 Image2(sub2ind(size(Image2), Region2.Boundary(1, :)-Min_x+1, ...\r
  18         Region2.Boundary(2, :)-Min_y+1)) = 1;\r
  19 To_start_fill = [Region2.Inner_pts(2, :)-Min_y+1; ...\r
  20         Region2.Inner_pts(1, :)-Min_x+1];\r
  21 Image2 = bwfill(Image2, To_start_fill(1, :), To_start_fill(2, :));\r
  22 \r
  23 Image1 = or(Image1, Image2);\r
  24 \r
  25 New_boundary = Extract_region_reg(Image1, 1);\r
  26 New_boundary = New_boundary + [Min_x-1; Min_y-1]*ones(1, size(New_boundary, 2));\r
  27 New_inner_pts = [Region1.Inner_pts, Region2.Inner_pts];\r
  28 New_shifts = [Min_x, Max_x; Min_y, Max_y];\r