Необъяснимые результаты после использования ind2sub в Matlab

Question

У меня проблемы с Matlab, я не понимаю. Следующий фрагмент кода анализирует коллекцию изображений и должен возвращать когерентное изображение (и всегда было).

Но поскольку я положил if-условие во втором for-loop (для целей оптимизации), он возвращает чересстрочное изображение.

Я не понимаю, почему, и я готов выбросить свой компьютер из окна. Я подозреваю, что это имеет какое-то отношение к ind2sub, но насколько я вижу, все работает отлично! Кто-нибудь знает, почему он это делает?

function imageMedoid(imageList, resizeFolder, outputFolder, x, y) 

    % local variables 
    medoidImage = zeros([1, y*x, 3]); 
    alphaImage = zeros([y x]); 
    medoidContainer = zeros([y*x, length(imageList), 3]); 

    % loop through all images in the resizeFolder 
    for i=1:length(imageList) 

     % get filename and load image and alpha channel 
     fname = imageList(i).name; 
     [container, ~, alpha] = imread([resizeFolder fname]); 

     % convert alpha channel to zeros and ones, add to alphaImage 
     alphaImage = alphaImage + (double(alpha)/255); 

     % add (r,g,b) values to medoidContainer and reshape to single line 
     medoidContainer(:, i, :) = reshape(im2double(container), [y*x 3]); 

    end 

    % loop through every pixel 
    for i=1:(y * x) 

     % convert i to coordinates for alphaImage 
     [xCoord, yCoord] = ind2sub([x y],i); 

     if alphaImage(yCoord, xCoord) == 0 

      % write default value to medoidImage if alpha is zero 
      medoidImage(1, i, 1:3) = 0; 

     else 

     % calculate distances between all values for current pixel 
     distances = pdist(squeeze(medoidContainer(i,:,1:3))); 

     % convert found distances to matrix of distances 
     distanceMatrix = squareform(distances); 

     % find index of image with the medoid value 
     [~, j] = min(mean(distanceMatrix,2)); 

     % write found medoid value to medoidImage 
     medoidImage(1, i, 1:3) = medoidContainer(i, j, 1:3); 

     end 

    end 

    % replace values larger than one (in alpha channel) by one 
    alphaImage(alphaImage > 1) = 1; 

    % reshape image to original proportions 
    medoidImage = reshape(medoidImage, y, x, 3); 

    % save medoid image 
    imwrite(medoidImage, [outputFolder 'medoid_modified.png'], 'Alpha', alphaImage); 

end 

Я не включил весь код, только эту функцию (для краткости), если кто-то нуждается в большем количестве (для лучшего понимания этого), пожалуйста, дайте мне знать, и я буду включать его.

Answer 1

Когда вы звоните ind2sub, вы даете размер [x y], но фактический размер alphaImage является [y x], так что вы не индексируют правильное расположение с xCoord и yCoord.