Indexing for partial matrix

Question

Stephen Thompson 2020-7-17

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/566757-indexing-for-partial-matrix

评论： Stephen Thompson 2020-7-21

I have a 140 x 2000 matrix. It contains a blob of area (I know the linear indices) that has known peaks. I want to find the maximum peaks in part of the blob. Specifically I want to find the maximum peak in rows 1-77 and the maximum in rows 78-140 of the blob. There may be other blobs with other peaks, so it has to be by blob area.

Should I generate the linear indices for each section (e.g. all columns, rows 1:77) and then find which are common with the blob indices? How would I do that?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Matt J 2020-7-17

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/566757-indexing-for-partial-matrix#answer_467289

编辑：Matt J 2020-7-17

[upperBlob,lowerBlob]=deal(yourMatrix);
upperBlob(78:end,:)=-inf;
lowerBlob(1:77,:)=-inf;
[imax1,jmax1]=find(upperBlob==max(upperBlob(:)));  %upper blob max
[imax2,jmax2]=find(lowerBlob==max(lowerBlob(:)));  %lower blob max

3 个评论
显示 1更早的评论隐藏 1更早的评论

Matt J 2020-7-21

编辑：Matt J 2020-7-21

I don't see why it wouldn't have worked. The area outside your blobs appears to be zero, while the area inside the blobs appears to be greater than zero, so the maxima I've computed above would inherently have to fall within the blobs. Regardless, you could exclude the background with a few additional lines:

A=-inf(size(yourImage));
A(blobIndices)=yourImage(blobIndices):
[upperBlob,lowerBlob]=deal(A);
upperBlob(78:end,:)=-inf;
lowerBlob(1:77,:)=-inf;
[imax1,jmax1]=find(upperBlob==max(upperBlob(:)));  %upper blob max
[imax2,jmax2]=find(lowerBlob==max(lowerBlob(:)));  %lower blob max

Stephen Thompson 2020-7-21

It is my fault for not explaining the problem well enough. Another example with 2 blobs - but I still would want them classed by row (corresponds for frequency but axes left generic) and by blob.

There could be multiple blobs with peaks above and below the row cut line, but I want to separate them into maxima by blob and by dichotomous row grouping, So it does have to be separated on a per-blob basis. Now that I can find the subscripts (ind2sub) for the blobs I can solve the problem. Thank you though for your solution.

请先登录，再进行评论。