Code Review Stack Exchange is a question and answer site for peer programmer code reviews. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

Currently, I have this code snippet, but it's a little bit slow:

% input image
input_matrix = imread('placa.bmp');
[rows cols] = size(input_matrix);

% rotation
degree = 30;
radians = (pi * degree) / 180; 
theta = radians;

% output matrix
t_matrix = input_matrix;
t_matrix(t_matrix == input_matrix) = NaN;

% transformation matrix
T = [cos(theta) -sin(theta) 0; sin(theta) cos(theta) 1; 0 0 1];

% loop over each input_matrix coordinate
for n = 1:numel(input_matrix)

    % current coordinate
    [x y] = ind2sub([rows cols], n);

    % transpose
    v = [x y 1]';

    % homogeneous coordinate
    v = T*v;

    % transponse again
    v = v';

    % only integer values
    a = floor(v(1));
    b = floor(v(2));

    if a > 0 && b > 0

        % replace in t_matrix
        t_matrix(a,b) = input_matrix(x,y);

% get only a part of t_matrix
t_matrix = t_matrix(1:rows, 1:cols);

figure; imshow(input_matrix); % original image
figure; imshow(t_matrix)

How can I improve it?

Condition: keep using this matrix:

% transformation matrix
T = [cos(theta) -sin(theta) 0; sin(theta) cos(theta) 1; 0 0 1];
share|improve this question

Rather than looping over all x and y values you could put them into a vector, that should speed things up considerably.


M=[x(:) y(:) ones(rows*cols,1)];

Now you just need to multiply M or M' with T or T' and you have subsitituted that part of the code.

share|improve this answer
It would help if you posted an example along with an explanation instead. – Trevor Pilley Nov 14 '12 at 11:01
Added the improved code part – Dennis Jaheruddin Nov 14 '12 at 12:16

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.