--- a
+++ b/util/rotateAround.m
@@ -0,0 +1,59 @@
+function output=rotateAround(image, pointY, pointX, angle, varargin)
+% ROTATEAROUND rotates an image.
+%   ROTATED=ROTATEAROUND(IMAGE, POINTY, POINTX, ANGLE) rotates IMAGE around
+%   the point [POINTY, POINTX] by ANGLE degrees. To rotate the image
+%   clockwise, specify a negative value for ANGLE.
+%
+%   ROTATED=ROTATEAROUND(IMAGE, POINTY, POINTX, ANGLE, METHOD) rotates the
+%   image with specified method:
+%       'nearest'       Nearest-neighbor interpolation
+%       'bilinear'      Bilinear interpolation
+%       'bicubic'       Bicubic interpolation
+%    The default is fast 'nearest'. Switch to 'bicubic' for nicer results.
+%
+%   Example
+%   -------
+%       imshow(rotateAround(imread('eight.tif'), 1, 1, 10));
+%
+%   See also IMROTATE, PADARRAY.
+
+%   Contributed by Jan Motl (jan@motl.us)
+%   $Revision: 1.2 $  $Date: 2014/05/01 12:08:01 $
+
+% Parameter checking.
+numvarargs = length(varargin);
+if numvarargs > 1
+    error('myfuns:somefun2Alt:TooManyInputs', ...
+        'requires at most 1 optional input');
+end
+optargs = {'nearest'};    % Set defaults for optional inputs
+optargs(1:numvarargs) = varargin;
+[method] = optargs{:};    % Place optional args in memorable variable names
+
+% Initialization.
+[imageHeight imageWidth ~] = size(image);
+centerX = floor(imageWidth/2+1);
+centerY = floor(imageHeight/2+1);
+
+dy = centerY-pointY;
+dx = centerX-pointX;
+
+% How much would the "rotate around" point shift if the 
+% image was rotated about the image center. 
+[theta, rho] = cart2pol(-dx,dy);
+[newX, newY] = pol2cart(theta+angle*(pi/180), rho);
+shiftX = round(pointX-(centerX+newX));
+shiftY = round(pointY-(centerY-newY));
+
+% Pad the image to preserve the whole image during the rotation.
+padX = abs(shiftX);
+padY = abs(shiftY);
+
+padded = padarray(image, [padY padX], 'symmetric');
+
+% Rotate the image around the center.
+rot = imrotate(padded, angle, method, 'crop');
+
+% Crop the image.
+output = rot(padY+1-shiftY:end-padY-shiftY, padX+1-shiftX:end-padX-shiftX, :);
+