问矩阵行中点取决于一定条件的Matlab
EN

Stack Overflow用户

提问于 2013-08-02 07:48:53

回答 3查看 787关注 0票数 1

我有一个大小为10x100的矩阵A，如下所示。我想做的是：

我将逐行工作，对于每一行，我将检查该行中每一个颜色的数据。
假设我现在在第一行的第一个单元格中。我将检查这个值是否为零，我将移到下一个col，以此类推，直到我找到一个具有非零值的col，并保存它的col数，例如col3“这意味着col &2是零”。
现在我在row1中的第一个非零值，我将移动到下一个，直到我找到一个零值的col。我会在这个零一之前把它拿来，它一定是一个非零的，并保存它。这意味着col 4和5&6是非零的，col8是零的。
现在，我想保存这两列(例如col3和col7 )之间的中间值，然后中间值是col5，所以我将保存索引row1_col5。如果有两个中间值，那么其中任何一个都可以。
然后，我将移到下一个行，直到找到一个非零的“从2->5执行相同的步骤”，直到第一行完成为止。
移到下一行，从步骤2->5重新开始。
有两条规则：-The第一条是，只有在至少有3个非零连续值的情况下，我才能得到非零连续值的中间索引，如果有两个非零的连续值，那么中间的值将不会被计算，第二个是，如果零连续值的数目小于3，那么它们将被忽略，并被视为非零值。例如，在下面的示例中，第一行中间值是col5和col11。在row2中，col5被计算在内，而row3中的cols不满足这个条件，而在row4中col6或col7将被计算在内。
在完成所有行之后，希望有一个向量或数组保存所有中间索引的位置，例如row1_col5、row1_col17、row2_col_10等。

示例：

A = [ 0 0 0 2 4 1 0 0 0 1 3 2;
      0 0 0 5 1 1 1 1 0 0 0 1;
      0 3 4 1 0 3 1 2 0 0 1 3;
      0 0 0 0 1 3 4 5 0 0 0 0];

对于第一行，中间值是5和11，等等，所以如果有人可以建议我如何用最少的处理来完成这个任务，因为可以使用循环来完成，但是如果有更有效的方法来完成呢？如果需要澄清，请告诉我。

matlab

回答 3

Stack Overflow用户

回答已采纳

发布于 2013-08-02 08:47:06

现在你澄清了你的问题(再次.)这里有一个解决方案(仍然使用for循环.)。它包括“规则7”--不包括少于三个元素的运算；它还包括该规则的第二部分--少于三个零的运算不算作零。新代码如下所示：

A = [ 0 0 0 2 4 1 0 0 0 1 3 2;
      0 0 0 5 1 1 1 1 0 0 0 1;
      0 3 4 1 0 3 1 2 0 0 1 3;
      0 0 0 0 1 3 4 5 0 0 0 0];

 retVal = cell(1, size(A, 1));
 for ri = 1:size(A,1)
     temp = [1 0 0 0 A(ri,:) 0 0 0 1]; % pad ends with 3 zeros + 1
                                       % so that is always a "good run"
     isz = (temp == 0); % find zeros - pad "short runs of 0" with ones
     diffIsZ = diff(isz);
     f = find(diffIsZ == 1);     
     l = find(diffIsZ == -1);
     shortRun = find((l-f)<3); % these are the zeros that need eliminating

     for ii = 1:numel(shortRun)
         temp(f(shortRun(ii))+1:l(shortRun(ii))) = 1;
     end

     % now take the modified row:
     nz = (temp(4:end-3)~=0);
     dnz = diff(nz); % find first and last nonzero elements
     f = find(dnz==1); 
     l = find(dnz==-1);
     middleValue = floor((f + l)/2);
     rule7 = find((l - f) > 2);
     retVal{ri} = middleValue(rule7);
 end

返回值必须使用单元格数组，因为您不知道每行将返回多少个元素(根据更新的需求)。

上面的代码返回以下单元格数组：

{[5 11], [6], [7], [7]}

我似乎仍然不理解您的“规则7"，因为您说”第3行中没有列满足此条件“。但在我看来，一旦我们消除了短距离的零，它就做到了。除非我误解了您想要处理的非零数的运行，并一直到边缘(我假设这是好的-这就是为什么您返回11作为一个有效的列在第1行；那么，您为什么不返回7的第3行？)

票数 3

Stack Overflow用户

发布于 2013-08-02 09:08:05

试试这个：

sizeA = size(A);
N = sizeA(1);
D = diff([zeros(1, N); (A.' ~= 0); zeros(1,N)]) ~= 0;
[a b] = find(D ~= 0);

c = reshape(a, 2, []);
midRow = floor(sum(c)/2);
midCol = b(1:2:length(b))

在此之后，midRow和midCol包含您的质心的指数(例如，对于上面给出的示例矩阵，midRow(1) = 1，midCol(1) =4)。

票数 1

Stack Overflow用户

发布于 2013-08-02 08:24:55

如果您不介意使用for循环：

A = [ 0 0 1 1 1 0 1;
      0 0 0 0 0 0 0;
      0 1 1 1 1 0 0;
      0 1 1 1 0 1 1;
      0 0 0 0 1 0 0]; % data

sol = repmat(NaN,size(A,1),1);
for row = 1:size(A,1)
  [aux_row aux_col aux_val] = find(A(row,:));
  if ~isempty(aux_col)
    sol(row) = aux_col(1) + floor((find(diff([aux_col 0])~=1,1)-1)/2);
    % the final 0 is necessary in case the row of A ends with ones
    % you can use either "floor" or "ceil"
  end
end

disp(sol)