MATLAB: How to divide an image into non-overlapping blocks

Image Processing Toolboxnon overlapping

For block pairing, I need to divide original block and target blocks into non-overlapping blocks.I tried with the function called blockproc(),but it is not suitable for images that i have taken.
% Divide the image up into 4 blocks.
% Let's assume we know the block size and that all blocks will be the same size.
blockSizeR = 128; % Rows in block.
blockSizeC = 128; % Columns in block.
% Figure out the size of each block. 
wholeBlockRows = floor(rows / blockSizeR);
wholeBlockCols = floor(columns / blockSizeC);
% Preallocate a 3D image
image3d = zeros(wholeBlockRows, wholeBlockCols, 3);
% Now scan though, getting each block and putting it as a slice of a 3D array.
sliceNumber = 1;
for row = 1 : blockSizeR : rows
  for col = 1 : blockSizeC : columns
    % Let's be a little explicit here in our variables
    % to make it easier to see what's going on.
    % Determine starting and ending rows.
    row1 = row;
    row2 = row1 + blockSizeR - 1;
    row2 = min(rows, row2); % Don't let it go outside the image.

    % Determine starting and ending columns.
    col1 = col;
    col2 = col1 + blockSizeC - 1;
    col2 = min(columns, col2); % Don't let it go outside the image.
    % Extract out the block into a single subimage.
    oneBlock = grayImage(row1:row2, col1:col2);
    % Specify the location for display of the image.
    subplot(2, 2, sliceNumber);
    imshow(oneBlock);
    % Make the caption the block number.
    caption = sprintf('Block #%d of 4', sliceNumber);
    title(caption, 'FontSize', fontSize);
    drawnow;
    % Assign this slice to the image we just extracted.
    if (row2-row1+1) == blockSizeR && (col2-col1+1) == blockSizeC
      % Then the block size is the tile size, 
      % so add a slice to our 3D image stack.
      image3D(:, :, sliceNumber) = oneBlock;
    else
      newTileSize = [(row2-row1+1), (col2-col1+1)];
      warningMessage = sprintf('Warning: this block size of %d rows and %d columns\ndoes not match 
the preset block size of %d rows and %d columns.\nIt will not be added to the 3D image stack.',...
        newTileSize(1), newTileSize(2), blockSizeR, blockSizeC);
      uiwait(warndlg(warningMessage));
    
end
    sliceNumber = sliceNumber + 1;
  end
end
Is there any specific range for blocksizeR and blocksizeC, can I assign different sizes? if I want non overlapping blocks, how these values are connected?

Best Answer

% Divide the image up into 4 blocks.
% Let's assume we know the block size and that all blocks will be the same size.
% these two values do not need to be the same
blockSizeR = 128; % Rows in block.
blockSizeC = 128; % Columns in block.
% Figure out the size of each block. 
wholeBlockRows = floor(rows / blockSizeR);
wholeBlockCols = floor(columns / blockSizeC);
trailing_rows = rows - wholeBlockRows * blockSizeR;
trailing_cols = cols - wholeBlockCols * blockSizeC;
image3d_cell = mat2cell( YourImage, [blockSizeR * ones(1, wholeblockRows), trailing_rows], [blockSizeC * ones(1, wholeBlockCols), trailing_cols], size(YourImage, 3) );
Now image3d_cell will be a cell array, wholeBlockRows or wholeBlockRows+1 (if there are trailing rows), by wholeBlockCols or wholeBlockCols+1 (if there are trailing columns).
if trailing_rows == 0 && trailing_cols == 0
  image3d = cat( ndims(YourImage)+1, image3d_cell{:} );
else
  image3d = [];
  warning('There were trailing rows or columns, could not form 3D image automatically');
end
Related Question