我希望矩阵中一个点坐标的所有邻居都与该点相距 d。
生成这样的子矩阵最有效的方法是什么？
例如如果A = [ 1 , 2 , 3 ; 4 , 5 , 6 ; 7 ,8 ,9]
那么对于d = 0 and data 5，我会得到5 的答案，但对于d = 1，我会得到矩阵A 作为结果，因为A 本身的大小是3X3。
同时，如果矩阵中不存在索引，那么此时我应该得到 0。
所以对于data point 2 and d = 1，答案是[ 0 , 0 , 0 ; 1 , 2 ,3 ; 4 , 5 , 6]
我该怎么做？
【问题讨论】：
标签： matlab
试试这个...
A = padarray( A , [d d] )
% assuming the data is at r row and c col do the following
A( r : r + 2d , c : c + 2 d )
我希望它有效...并没有真正检查！
【讨论】：
这里有一个解决方案。我会让你把它变成一个函数。让n 表示子数组中心的值，在您的示例中为5。
z = zeros(size(A,1)+2*d,size(A,2)+2*d);
z(d+1:d+size(A,1),d+1:d+size(A,2)) = A;
[r,c] = find(z==n);
z(r-d:r+d,c-d:c+d)
编辑
试试这个版本，它和我之前的尝试一样经过充分测试。如果 find 表达式返回多个位置，这仍然无法解决。
我仍然不保证这样做的效率。
【讨论】：
从你在 A 中的情况开始，而不是给 d=0 和 data_point=5，你应该给 data_point=[2,2]。如果您将其指定为5，您可以通过
ij = [ceil(data_point/size(A,1)), rem(data_point-1, size(A,1))+1]
或使用 find，如 High 的回答，取决于如何解释您的问题。
下面的函数
function B = find_neigh(A,ij,d)
imin = max(ij(1)-d,1);
imax = min(ij(1)+d,size(A,1));
jmin = max(ij(2)-d,1);
jmax = min(ij(2)+d,size(A,2));
B = zeros(2*d+1,2*d+1);
r1 = max(2-ij(1)+d,1);
c1 = max(2-ij(2)+d,1);
B(r1:r1+imax-imin, c1:c1+jmax-jmin) = A(imin:imax,jmin:jmax);
应该返回你想要的结果。例如用
调用它A = [1, 2, 3; 4, 5, 6; 7, 8, 9];
find_neigh(A,[1,2],1)
编辑：修复了两个，没有三个错误，并使其成为一个功能
【讨论】：