Is it possible to detect a specefic form of matrices from a big one?

Question

mika 2013-11-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/107065-is-it-possible-to-detect-a-specefic-form-of-matrices-from-a-big-one

评论： mika 2014-4-21

Hello, I need to find those 5 forms of matrices from a binary 'big' matrice

X=[1 0 1 ;
   0 1 0 ;
   1 0 1 ;];
Y=[1 0 1 ;
   0 1 0 ;
   1 0 0 ;];
O=[0 1 0 ;
   1 0 1 ;
   0 1 0 ;];
-=[0 0 0 ;
   1 1 1 ;
   0 0 0 ;];
L=[1 0 0 ;
   1 0 0 ;
   1 1 1 ;];

Any idea? help is much appreciated

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

mika 2013-11-21

it's 1000*1001

talha 2013-11-21

how to connect blocks of two different library together ?

请先登录，再进行评论。

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

Answer 1

Cedric 2013-11-21

4
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/107065-is-it-possible-to-detect-a-specefic-form-of-matrices-from-a-big-one#answer_115937

编辑：Cedric 2013-11-22

在 MATLAB Online 中打开

EDITED: got rid of a side effect in the dash case (illustrated in the 4th picture below.. I didn't update it though) by using 2*A-1 instead of A and numel(X) instead of nnz(X).

Someone might provide you with a solution based on the Image Processing Toolbox (which I am not familiar with). Here is one potential solution and a hint in the meantime..

1. Based on CONV2

    A = double( rand(10) > 0.4 ) ;      % Small example.
    X = [1 0 1 ; ...
         0 1 0 ; ...
         1 0 1 ] ;
    ker = rot90(rot90( 2*X -1 )) ;
    [r, c] = find( conv2(2*A-1, ker, 'same') == numel(X) ) ;

With that and a few runs (need a little luck for RAND to output one or more fillings with matches):

 >> A
 A =
   1     1     1     1     1     0     1     0     1
   0     0     1     1     0     1     1     0     1
   0     0     0     0     1     0     1     0     1
   0     1     0     1     0     1     1     1     1
   1     0     1     0     1     1     1     1     0
   1     1     0     1     0     0     1     1     0
   1     1     0     1     1     0     0     0     1
   1     0     0     1     0     1     1     1     1
   1     0     1     1     1     1     0     1     1
   1     0     1     0     1     1     1     0     0
 >> [r,c]
 ans =
   4
   6

which indicates two matches, centered in (5,4) and (3,6). This approach takes ~13ms on my laptop with a 1000x1001 A matrix.

2. By looping over elements of patterns and vectorizing tests on A rather than the opposite.

.. I leave that for later if CONV2 is not suitable.

--------------------------------------------------------------------------------------------------------------------------------

EDIT: here is a more elaborate example..

 n = 50 ;
 A = double( rand(n, n+1) > 0.5 ) ;
 B = 2*A - 1 ;
 patterns = {[1 0 1; 0 1 0; 1 0 1]; ... % X
             [1 0 1; 0 1 0; 1 0 0]; ... % Y
             [0 1 0; 1 0 1; 0 1 0]; ... % O
             [0 0 0; 1 1 1; 0 0 0]; ... % -
             [1 0 0; 1 0 0; 1 1 1]} ;   % L
 labels   = {'X', 'Y', 'O', '-', 'L'} ;
 matches  = cell( size(labels )) ;
 for pId = 1 : numel( patterns )
     nel    = numel( patterns{pId} ) ;
     ker    = 2*patterns{pId} - 1 ;    
     [r, c] = find( conv2(B, rot90(rot90(ker)), 'same') == nel ) ;
     matches{pId} = [r, c] ;
     figure(pId) ;  clf ;  hold on ;
     spy(A) ;
     plot( c, r, 'r+', 'MarkerSize', 20 ) ;
     title( sprintf( 'Matches for "%s" pattern', labels{pId} )) ;
     set( gcf, 'Units', 'normalized' ) ;
     set( gcf, 'Position', [0 0 0.4 0.7] ) ;
 end

which outputs the following figures..

6 个评论
显示 4更早的评论隐藏 4更早的评论

Cedric 2013-11-26

编辑：Cedric 2013-11-26

在 MATLAB Online 中打开

Just to be sure, what you want to detect with e.g. the Y pattern, is really a location in the array that corresponds to your picture above where values at the location and its direct neighbors are exactly

0 1
1 0
0 0

and not all possible rotations of this array.. ? Also, are you really looking for 3x3 pixels regions matching these patterns that you provide, or is the goal ultimately to split your skeleton into nodes and vertices using the skeleton to define the topoly (in which case you should look more into the Image Processing Toolbox help, e.g. BWMORPH)?

To understand my solution, you need to understand an operation called "convolution". The simplest illustration that I can provide is that if

C = conv2(B, rot90(rot90(ker)), 'same') ;

then you have element (2,2) of C, for example, which is numerically equal to

sum( sum( B(1:3,1:3) .* pattern ))

If B(1:3,1:3) and ker as both defined as (the Y pattern with 0's replaced by -1)

   1 -1  1
  -1  1 -1
   1 -1 -1

their element-wise product is

1  1
1  1
1  1

which sums up to 9. If there is any mismatch between B(1:3,1:3) and ker, it introduces negative values in the element-wise product, which lowers the sum. This is why we look for elements of the convolution which are equal to the number of elements in the pattern (=9). You can try this illustration with A and X for example and see that their 0 entries are introducing false positives (e.g. X patterns are positive when testing for Y); this is why B and ker have to be built by replacing 0's by -1.

Cedric 2013-11-28

编辑：Cedric 2013-11-28

Ok, this is quite different from what the initial question suggested. I would advise you to re-post the question including images of the skeleton and zooms over regions which represent the different patterns. You should also attach the source array or image so people can perform tests.

mika 2014-1-13

Hello,

Can you tell me please how to make the output into the same plot. i'm trying to make the result as one plot but in vain can you help me please.

thank you.

请先登录，再进行评论。

Answer 2

talha 2013-11-21

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/107065-is-it-possible-to-detect-a-specefic-form-of-matrices-from-a-big-one#answer_115942

how to connect blocks of two different library together ?

2 个评论
显示无隐藏无

Cedric 2013-11-21

Please post a new question and provide more details about your issue.

mika 2014-4-21

okk sorry for late reply

i' m still looking for the solution

请先登录，再进行评论。

Answer 3

Image Analyst 2013-11-26

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/107065-is-it-possible-to-detect-a-specefic-form-of-matrices-from-a-big-one#answer_116350

normxcorr2_demo.m

Try normalized cross correlation (see attached demo below in blue text), or the hit or miss transform, done by bwhitmiss() in the Image Processing Toolbox.