General method to find submatrix in matlab matrix

by

Here is low-performance, but (supposedly) arbitrary dimensional function. It uses find to create a list of (linear) indices of potential matching positions in total and then just checks if the appropriately sized subblock of total matches sub.

function loc = matrixFind(total, sub)
%matrixFind find position of array in another array

    % initialize result
    loc = [];

    % pre-check: do all elements of sub exist in total?
    elements_in_both = intersect(sub(:), total(:));
    if numel(elements_in_both) < numel(unique(sub))
        % if not, return nothing
        return
    end

    % select a pivot element
    % Improvement: use least common element in total for less iterations
    pivot_element = sub(1);

    % determine linear index of all occurences of pivot_elemnent in total
    starting_positions = find(total == pivot_element);

    % prepare cell arrays for variable length subscript vectors
    [subscripts, subscript_ranges] = deal(cell([1, ndims(total)]));


    for k = 1:length(starting_positions)
        % fill subscript vector for starting position
        [subscripts{:}] = ind2sub(size(total), starting_positions(k));

        % add offsets according to size of sub per dimension
        for m = 1:length(subscripts)
            subscript_ranges{m} = subscripts{m}:subscripts{m} + size(sub, m) - 1;
        end

        % is subblock of total equal to sub
        if isequal(total(subscript_ranges{:}), sub)
            loc = [loc; cell2mat(subscripts)]; %#ok<AGROW>
        end
    end
end

More Related Contents:

  1. How can I index a MATLAB array returned by a function without first assigning it to a local variable?
  2. Matrix of unknown length in MATLAB?
  3. How do you concatenate the rows of a matrix into a vector?
  4. How can I accumulate cells of different lengths into a matrix in MATLAB?
  5. How do I visualize a matrix with colors and values displayed?
  6. Is MATLAB row specific or column major?
  7. How can I change the values of multiple points in a matrix?
  8. Construct adjacency matrix in MATLAB
  9. Replace specific columns in a matrix with a constant column vector
  10. How can I apply a function to every row/column of a matrix in MATLAB?
  11. Get the indices of the n largest elements in a matrix
  12. Multiply a 3D matrix with a 2D matrix
  13. How do I resize a matrix in MATLAB?
  14. Is there a substitute for blockproc in Matlab?
  15. Matrix “Zigzag” Reordering
  16. How do I compare all elements of two arrays?
  17. How can I divide each row of a matrix by a fixed row?
  18. How to subtract a vector from each row of a matrix? [duplicate]
  19. Building a matrix by merging the same row vector multiple times
  20. Matlab: repeat every column sequentially n times [duplicate]
  21. Index of max and min value in an array
  22. How can I count the number of elements of a given value in a matrix?
  23. How do I divide matrix elements by column sums in MATLAB?
  24. Does Matlab eig always returns sorted values?
  25. MATLAB find and apply function to values of repeated indices
  26. the easiest way to convert matrix to one row vector [duplicate]
  27. Extract arbitrarily rotated plane of data from 3D array as 2D array
  28. MATLAB: Subtracting matrix subsets by specific rows
  29. How do I convert a 2X2 matrix to 4X4 matrix in MATLAB?
  30. Efficiently compute a 3D matrix of outer products – MATLAB

Leave a Comment