Abstract-An iterative algorithm to perform the J-spectral factorization of a para-Hermitian matrix is presented. The algorithm proceeds by computing a special kernel representation of an interpolant for a sequence of points and associated directions determined from the spectral zeroes of the to-befactored matrix.