We present a numerical approach to compute a minimal periodic state-space realization of a transfer-function matrix corresponding to a lifted state-space representation. The proposed method determines a realization with time-varying state dimensions by using exclusively orthogonal transformations. The new method is numerically reliable, computationally efficient and thus well suited for robust software implementations.