Assuming that the spatial frequency fo=n/D where n is an integer, the fo component in the spectrum appears only in one interval, and hence C(f-fo) and C*(f+fo) in the equation (10) are prevented from being largely spread, so that the filtering and the spectrum shifting can accurately be carried out and the computations can be simplified.