Obtain subscript values of common diagonal rectangle in binary matrix?

조회 수: 1 (최근 30일)
I am trying to find where in the MM_bin binary matrix (see attached MM_bin.mat file) contains the check_mat matrix (see code below), and then extract the subscript coordinates in the MM_bin matrix of the 1s from the check_mat matrix. So far I am able to use a 2D convolution to extract the centroids of the spots where the check_mat matrix is found in the MM_bin matrix, but I am unsure how to extract the subscript coordinates in the MM_bin matrix for each centroid.
load('MM_bin.mat')
check_mat = [0 0 0 0 0 1 1 1 1;
0 0 0 0 1 1 1 1 0;
0 0 0 1 1 1 1 0 0;
0 0 1 1 1 1 0 0 0;
0 1 1 1 1 0 0 0 0;
1 1 1 1 0 0 0 0 0];
minLi = 6; % length of the 1s diagonal
minLj = 4; % width of the 1s diagonal
[i_cent, j_cent] = find(conv2(MM_bin(:,:,1), check_mat, 'same') == minLi*minLj)
  댓글 수: 5
Andrew Poissant
Andrew Poissant 2018년 10월 12일
Hm, maybe I am using it incorrectly. But I want to find the pattern of 1s seen in the check_mat matrix. The 0s are simply placeholders/are there to make the desired diagonal pattern of 1s. I just want to find where the pattern of 1s seen in the check_mat matrix is found in the MM_bin matrix. Maybe I should try to ask the question a different way
Guillaume
Guillaume 2018년 10월 12일
If you're just looking for the pattern of ones and you don't care if the 0s in your check_mat match a 0 or a 1, then your code is correct and I misunderstood. If the 0s must match a 0, then it cannot work with a convolution.
Assuming your code is correct, I still don't understand what final out you want. Perhaps, provide an example with smaller matrices.
E.G., with
MM_bin = [0 1 0 1 0
0 1 1 1 1
1 0 1 1 1]
check_mat = [0 1
1 0]
what final output do you want?

댓글을 달려면 로그인하십시오.

채택된 답변

Bruno Luong
Bruno Luong 2018년 10월 12일
편집: Bruno Luong 2018년 10월 12일
% Fake data
MM_bin=rand(100)>0.2;
check_mat = [0 0 0 0 0 1 1 1 1;
0 0 0 0 1 1 1 1 0;
0 0 0 1 1 1 1 0 0;
0 0 1 1 1 1 0 0 0;
0 1 1 1 1 0 0 0 0;
1 1 1 1 0 0 0 0 0];
s = size(check_mat);
% Careful if your pattern is not symmetric you must flip it
% for each dimension when using with CONV
[i_cent, j_cent] = find(conv2(MM_bin(:,:,1), ...
fliplr(flipud(check_mat)), 'same') == sum(check_mat(:)));
[ip,jp] = find(check_mat);
i = floor(i_cent(:)-s(1)/2)+ip(:).';
j = floor(j_cent(:)-s(2)/2)+jp(:).';
% remove redundancy
ij = unique([i(:) j(:)],'rows');
% graphical check
close all
imagesc(MM_bin);
hold on
plot(ij(:,2),ij(:,1),'.r') % NOTE: x-axis is second dimension
axis equal
  댓글 수: 1
Andrew Poissant
Andrew Poissant 2018년 10월 12일
Beautiful! This was very helpful thank you! So I was on the right track but needed some tweaks.

댓글을 달려면 로그인하십시오.

추가 답변 (3개)

Image Analyst
Image Analyst 2018년 10월 12일
To find the centroid of the 1's, if you have the Image Processing Toolbox, you can use regionprops
check_mat = [0 0 0 0 0 1 1 1 1;
0 0 0 0 1 1 1 1 0;
0 0 0 1 1 1 1 0 0;
0 0 1 1 1 1 0 0 0;
0 1 1 1 1 0 0 0 0;
1 1 1 1 0 0 0 0 0];
props = regionprops(check_mat, 'Centroid');
xCentroid = props.Centroid(1)
yCentroid = props.Centroid(2)
  댓글 수: 9
Bruno Luong
Bruno Luong 2018년 10월 12일
편집: Bruno Luong 2018년 10월 12일
Or simply change the single conv2 check to
conv2(Im,p,'same')==sum(p(:)) & conv2(Im,1-p,'same')==sum(1-p(:))
since they are both binary images.
Image Analyst
Image Analyst 2018년 10월 12일
It's not clear to me what "(minus the 0s)" means. Maybe "ignoring any zeros"??? In other words, as long as that parallelogram is there, it doesn't matter if, outside the parallelogram, it's all ones, all zeros, or some random pattern or ones and zeros.

댓글을 달려면 로그인하십시오.


Image Analyst
Image Analyst 2018년 10월 12일
If you want to find where check_mat occurs in a much larger matrix, simple scan and use isequal().

Matt J
Matt J 2018년 10월 12일
편집: Matt J 2018년 10월 12일
[I0,J0]=find(check_mat);
I = I0 + (i_cent.'-3); %I subscripts
J = J0 + (j_cent.'-5); %J subscripts

카테고리

Help CenterFile Exchange에서 Mathematics and Optimization에 대해 자세히 알아보기

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by