function [cropped] = Crop(I,dim) dim; sz = size(I); center = [(sz(1)+mod(sz(1),2))/2 (sz(2)+mod(sz(2),2))/2]; offsetRow = center(1)-floor((dim(1)-(1-mod(dim(1),2)))/2); offsetCol = center(2)-floor((dim(2)-(1-mod(dim(2),2)))/2); cropped = I(offsetRow:offsetRow+dim(1)-1,offsetCol:offsetCol+dim(2)-1);