function reflectim = caoReflectPhaseSpectrum(im)

dims = size(im);

reflectim = im;

% fix the center column
 reflectim(dims(1):-1:(dims(1)/2+2),dims(2)/2+1) = -reflectim(2:dims(1)/2,dims(2)/2+1);

% reflect the 2nd half of the phase spectrum

reflectim(dims(1):-1:2,dims(2):-1:(dims(2)/2+2)) = -im(2:dims(1),2:dims(2)/2);

