function [ cropped ] = crop( im, block_shape )
    BH = block_shape(1);
    BW = block_shape(2);
    [H, W, C] = size(im);
    Rows = floor(H/BH);
    Cols = floor(W/BW);
    
    cropped = im(1:Rows*BH, 1:Cols*BW, :);