在 Octave 中重叠二维矩阵

Question

假设我有一个矩阵。这个矩阵是空白的，除了一些创建相同大小的矩形的点。创建矩形的每个点都有一个正数，该正数对于该矩形的所有点都相等。这个矩阵可能有几个矩形，每个矩形都有相同的数字。

现在想象一下，我有第二个具有相同结构的矩阵，但现在矩形位于其他位置，具有其他大小和其他值，就像图片中的矩阵 A 和 B。

我想合并两个矩阵，所以如果某个矩形有交集，则删除最小数量的矩形。

最好的方法是什么？我认为 fors 是可行的，但我想做 30 次，所以如果我不改进它，Octave 将需要很长时间才能完成它。

Answer 1

我编写了这个似乎可以完成这项工作的小代码：

A=[0 0 1 1 1 0 ;
   0 0 1 1 1 0 ;
   0 0 0 0 0 0 ;
   2 2 0 0 0 0 ;
   2 2 0 0 0 0 ;
   2 2 0 0 0 0 ];

B=[3 0 2 2 0 0 ;
   3 0 2 2 0 0 ;
   0 0 2 2 0 0 ;
   0 0 0 0 0 0 ;
   1 1 1 1 1 0 ;
   1 1 1 1 1 0 ];


%give a unique id to each rectangle
lblA=bwlabel(A);
lblB=bwlabel(B);

%search for overlap
overlap= A & B;

% find out overlaping rectangle s unique id
overlapLbl= unique([lblA(overlap(:)) lblB(overlap(:))],'rows');

%eliminate rectangle with lower score
for ov=1:size(overlapLbl,1)
  scoreA=A(find((lblA==overlapLbl(ov,1)),1));
  scoreB=B(find((lblB==overlapLbl(ov,2)),1));
  if scoreA>scoreB
    B(B==scoreB)=0;
  else
    A(A==scoreA)=0;
  end
end
C=A+B

输出：

C =
     3     0     2     2     0     0
     3     0     2     2     0     0
     0     0     2     2     0     0
     2     2     0     0     0     0
     2     2     0     0     0     0
     2     2     0     0     0     0