Solved problem by first doing a coordinate transformation and then mirroting across the new y-axis.
z=[cos(theta) sin(theta);-sin(theta) cos(theta)]; % Calculate transformation matrix
m=[x;y]; % Combine x & y coordinates in a single matrix
k=z*m; % Multiply matrices to get x & y values wrt new coordinate system
x_new=k(1,:);
y_new=k(2,:);
mirrored_x=-fliplr(x_new);
mirrored_y=fliplr(y_new);