[XR,YR,xr,yr,THETA] = rotateData(X,Y,X0,Y0,THETA,DIRECTION)

Rotate coordinates especified by [X,Y] around the point [X0,Y0], by an
angle (in radians) defined by THETA and to the direction especified by
DIRECTION.

THETA is a scalar if the user wishes to define it. THETA can also be
calculated based on the slope defined by two points in [X,Y]. If the
angle is to be calculated, THETA has to be a two-element vector with
the integer indices between which the linear slope in [X,Y] will be
calculated. If defined as an output, THETA calculated (or provided) in
the function is also returned.

DIRECTION should be either 'clockwise' (default) or 'anticlockwise'.

[XR,YR] are the coordinates rotated relative to the input coordinates'
origin. [xr,yr] are the coordinates translated to the origin [XO,YO]
used for the rotation.

Rafael Guedes, UoW
26/11/2010