Matlab - добавление круга к каждому кадру фильма

Question

Я хочу добавить круг вокруг частицы, движущейся в фильме. Я даю в координатах частицы для каждого кадра (переменная Traj). Однако следующий код рисует сплошной прямоугольник поверх моей частицы для каждого кадра. Какие-либо предложения?

function movie = AurelieCreateMovie3(Filename, Traj) 
    L = 369; 
    color = [1:3]; 
    R = 7; 
    info = imfinfo(Filename); 
    W = info.Width; 
    H = info.Height; 

    movie = uint8(zeros(H,W,3,L)); 

    frame = Traj(:,1); 
    x = round(Traj(:,2)); 
    y = round(Traj(:,3)); 
    L = length(x); %length of this track 
    dt = 0:0.01:2*pi; 

    for k = 1:L 
     movie(:,:,:,k) = imread(Filename,frame(k));  
     x_circle = R.* cos(dt) + x(k); 
     y_circle = R.* sin(dt) + y(k); 
     movie(round(y_circle), round(x_circle), color, k) = 255; 
    end 

    mov = immovie(movie); 
    movie2avi(mov, 'moviefinal', 'compression', 'None', 'fps', 23); 
    implay(mov); 
end 

Answer 1

При присвоении значений с matrix(rows, columns) = value, MATLAB назначает каждую перестановку элементов в rows и columns, а не просто назначая парную комбинацию из двух. Вместо этого вы хотите, чтобы преобразовать в линейный индекс, используя sub2ind

M = reshape(movie(:,:,:,k), [], size(movie, 3)); 
inds = sub2ind(size(M(:,:)), round(y_circle), round(x_circle)); 
M(inds,:) = 255; 
movie(:,:,:,k) = M;