Октава: Несколько подматриц из матрицы

Question

У меня есть большая матрица, из которой я хотел бы собрать коллекцию подматриц. Если моя матрица равна NxN, а размер подматрицы - MxM, я хочу собрать подматрицы I=(N - M + 1)^2. Другими словами, мне нужна одна подматрица MxM для каждого элемента исходной матрицы, которая может находиться в верхнем левом углу такой матрицы.

Вот код, у меня есть:

for y = 1:I 
    for x = 1:I 
     index = (y - 1) * I + x; 
     block_set(index) = big_mat(x:x+M-1, y:y+M-1) 
    endfor 
endfor 

выход, если а) не прав, и б) подразумевающее есть что-то в выражении big_mat(x:x+M-1, y:y+M-1), что может заставить меня, что я хочу, не нуждаясь в два для петель. Любая помощь была бы очень благодарна

Answer 1

В коде есть несколько неправильных вещей. Вот как я это сделать, если бы я использовать двойной цикл:

M = someNumber; 
N = size(big_mat,1); %# I assume big_mat is square here 

%# you need different variables for maxCornerCoord and nSubMatrices (your I) 
%# otherwise, you are going to index outside the image in the loops! 
maxCornerCoord = N-M+1; 
nSubMatrices = maxCornerCoord^2; 

%# if you want a vector of submatrices, you have to use a cell array... 
block_set = cell(nSubMatrices,1); 
%# ...or a M-by-M-by-nSubMatrices array... 
block_set = zeros(M,M,nSubMatrices); 
%# ...or a nSubMatrices-by-M^2 array 
block_set = zeros(nSubMatrices,M^2); 

for y = 1:maxCornerCoord 
    for x = 1:maxCornerCoord 
     index = (y - 1) * maxCornerCoord + x; 
     %# use this line if block_set is a cell array 
     block_set{index} = big_mat(x:x+M-1, y:y+M-1); 
     %# use this line if block_set is a M-by-M-by-nSubMatrices array 
     block_set(:,:,index) = big_mat(x:x+M-1, y:y+M-1); 
     %# use this line if block_set is a nSubMatrices-by-M^2 array 
     block_set(index,:) = reshape(big_mat(x:x+M-1, y:y+M-1),1,M^2); 
    endfor 
endfor 

EDIT

Я только что видел, что есть реализация im2col для Октава. Таким образом, вы можете переписать двойную петлю как

%# block_set is a M^2-by-nSubMatrices array 
block_set = im2col(big_mat,[M,M],'sliding'); 

%# if you want, you can reshape the result to a M-by-M-by-nSubMatrices array 
block_set = reshape(block_set,M,M,[]); 

Это, вероятно, быстрее, и экономит много цифровых деревьев.

Answer 2

О вашем выходе, являющемся неправильным, возможно, это связано с назначением. Вы пытаетесь присвоить матрицу векторную позицию. Попробуйте использовать

block_set(:,:,index) = big_mat(x:x+M-1, y:y+M-1) 

вместо этого.

Answer 3

С помощью Mathematica: Этот код представляет собой матрицу, в которой каждый элемент является матрицей MxM с каждым элементом в исходной матрице в верхнем левом углу такой матрицы.

Матричные элементы справа и внизу дополняются х.

Partition[big_mat, {M, M}, {1, 1}, {1, 1}, x] 

Пример: alt text http://img130.imageshack.us/img130/6203/partitionf.png

Если оставить от й аргумента, то тогда он автоматически образцы периодически.