Пример 1D векторов из 3D-массива с использованием вектора точек

Question

У меня есть изображение канала, и у меня есть матрица точек размером 100x2 (в моем случае n равно 20, но, возможно, более понятно думать об этом как о 3-канальном изображении) , Мне нужно пробовать изображение в каждой точке и получить массив nx100 этих точек изображения. Я знаю, как сделать это с цикл:

for j = 1:100 
     samples(j,:) = image(points(j,1),points(j,2),:); 
end 

Как бы я векторизации это? Я пробовал

samples = image(points); 

но это дает 200 образцов из 20 каналов. И если я попробую

samples = image(points,:); 

Это дает мне 200 образцов из 4800 каналов. Даже

samples = image(points(:,1),points(:,2)); 

дает мне 100 х 100 образцов 20 (по одному для каждой возможной комбинации х в X и Y в Y)

Answer 1

короткий путь, чтобы сделать это было бы, чтобы изменить свой образ, так что вы заставьте ваше изображение [nRows, nCols, nChannels] быть [nRows*nCols, nChannels]. Затем вы можете преобразовать массив points в линейный индекс (используя sub2ind), который будет соответствовать новому «комбинированному» индексу строки. Затем, чтобы захватить все каналы, вы можете просто использовать оператор двоеточия (:) для второго измерения, которое теперь представляет каналы.

% Determine the new row index that will correspond to each point after we reshape it 
sz = size(image); 
inds = sub2ind(sz([1, 2]), points(:,2), points(:,1)); 

% Do the reshaping (i.e. flatten the first two dimensions) 
reshaped_image = reshape(image, [], size(image, 3)); 

% Grab the pixels (rows) that we care about for all channels 
newimage = reshaped_image(inds, :); 

size(newimage) 

    100 20 

Теперь у вас есть изображение, отобранное в точках вы хотели для всех каналов.