获取两个不同矩阵中给定距离内的点

Question

我有两个矩阵A 和B，其中行数可以变化。 A 和 B 不一定具有相同的行数。

例如：

A = [ 110    90
      130   140
      230    50
      370   210 ];
B = [ 321    95
      102    35
      303   200 ];

现在矩阵A 和B 有“对应点”。对应点是两个矩阵的第二列中的值在+/-20 内的行。

例如：

• A(1,2) = 90 和B(1,2) = 95，区别在+/-20 之内，所以A(1,:) 和B(1,:) 是对应点。
• A(2,2) = 140 和 B(2,2) = 35，区别是 not 在 +/-20 内，所以 A(2,:) 和 B(2,:) 是不对应点。
• A(3,2) = 50 和B(2,2) = 35，区别在+/-20 之内，所以A(3,:) 和B(2,:) 是对应点。

使用这个我想将A和B的对应点分别存储在C和D中。对于上面的示例，最终矩阵应如下所示：

C = [ 110    90
      230    50
      370   210 ]
D = [ 321    95
      102    35
      303   200 ]

Answer 1

您可以使用pdist2获取所有距离

dists = pdist2( A(:,2), B(:,2) )
>> dists = [  5    55   110
             45   105    60
             45    15   150
            115   175    10 ]

然后获取所有“对应点”的索引，由阈值 20 定义。

% Get combinations within tolerance
idx = dists < 20;
% Get indices  
[iA, iB] = find(idx);

然后你可以创建最终的矩阵

C = A(iA, :);
D = B(iB, :);

---

编辑：确保每个配对都是唯一的一种方法（即A(1,:) 不能与来自B 的多行配对）是为每行/列获取最小值dists .注意：如果距离完全相同，这仍然会给你重复匹配，你还没有定义应该如何处理。

dists = pdist2( A(:,2), B(:,2) );
% Set values which are greater than the row/column minima to be infinity.
% This means they will never be within the tolerance of 20 (or whatever else)
dists ( bsxfun(@gt, dists, min(dists,[],1)) | bsxfun(@gt, dists, min(dists,[],2)) ) = Inf;
% In MATLAB versions > 2016b, you can use implicit expansion to replace bsxfun
% That would be: dists( dists > min(dists,[],1) | dists > min(dists,[],2) )

% Now continue as before
[iA, iB] = find( dists < 20 );
C = A(iA, :);
D = B(iB, :);