1

在Matlab中，我有矩阵，在我的代码的前一个阶段，选择了一个特定的元素。从矩阵的这一点开始，我想找到一个最大值，不仅仅是给定半径的所有周围邻居之间的最大值，而是给定方向角的最大值。让我用一个例子来解释一下：

这是矩阵 A：

一个=

 0     1     1     1     0     0     9     1     0     
 0     2     2     4     3     2     8     1     0     
 0     2     2     3     3     2     2     1     0     
 0     1     1     3     2     2     2     1     0     
 0     8     2     3     3     2     7     2     1    
 0     1     1     2     3     2     3     2     1

第一阶段选择的元素是 A(2,4) 中的4，下一个元素应该是最大值，例如 315 度方向角，即A(5,7) 中的7 .

我所做的是，根据角度，将矩阵 A 细分为不同的象限，并仅使用该象限的值创建一个新矩阵（A 的子矩阵）。

所以，对于这个例子，子矩阵将是 A 的第四象限：

q_A =

 4     3     2     8     1     0     
 3     3     2     2     1     0     
 3     2     2     2     1     0     
 3     3     2     7     2     1     
 2     3     2     3     2     1

现在，这是我的问题，我怎样才能提取7？

我唯一能做的（并且有效）是找到所有超过阈值的值，然后计算这些点的方向。然后，保存与给定方向（本例中为 315 度）具有相似方向的所有值，最后找到其中的最大值。它有效，但我想可能会有更快和“更清洁”的解决方案。

4

2 回答 2

1

这是我的理论，但我没有图像处理工具箱来测试它。也许有人可以发表评论？

%make (r,c) the center by padding with zeros
    if r > size(A,1)/2
        At = padarray(A, [size(A,1) - r], 0, 'pre');
    else
        At = padarray(A, [r-1], 0 'post');

    if c > size(A,2)/2
        At = padarray(At, [size(A,2) - c], 0, 'pre');
    else
        At = padarray(At, [c-1], 0 'post');

%rotate by your angle (maybe this should be -angle or else 360-angle or 2*pi-angle, I'm not sure
    Ar = imrotate(At,angle, 'nearest', 'loose'); %though I think nearest and loose are defaults

%find the max

    max(Ar(size(Ar,1)/2, size(Ar,2)/2:end); %Obviously you must adjust this to handle the case of odd dimension sizes.

此外，根据您的阵列要求，填充-inf可能比0

于 2013-08-07T11:20:37.190 回答
1

以下是该问题的相对便宜的解决方案，尽管我发现将头绕在矩阵坐标系周围真的很痛苦，并且可能有一些空间来整理它。它只是以提供的角度沿着起点周围的一条线跟踪所有矩阵条目（所有坐标和角度当然基于矩阵索引单位）：

A = [ 0     1     1     1     0     0     9     1     0
      0     2     2     4     3     2     8     1     0
      0     2     2     3     3     2     2     1     0
      0     1     1     3     2     2     2     1     0
      0     8     2     3     3     2     7     2     1
      0     1     1     2     3     2     3     2     1 ];
alph = 315;
r = 2;
c = 4;
% generate a line through point (r,c) with angle alph
[nr nc] = size(A);
x = [1:0.1:nc]; % overkill
m = tan(alph);
b = r-m*c;
y = m*x + b;
crd = unique(round([y(:) x(:)]),'rows');
iok = find((crd(:,1)>0) & (crd(:,1)<=nr) & (crd(:,2)>0) & (crd(:,2)<=nc));
crd = crd(iok,:);
indx=sub2ind([nr,nc],crd(:,1),crd(:,2));
% find max and position of max
[val iv]=max(A(indx)); % <-- val is the value of the max
crd(iv,:)             % <-- matrix coordinates (row, column) of max value

结果：

val =

     7


iv =

     8


ans =

     5     7
于 2013-08-07T15:08:24.413 回答