Photoelectron spectra of the cis and trans isomers of HCOH were computed using vibrational wave functions calculated by diagonalizing the Watson Hamiltonian, including Up to four mode couplings. The full-dimensional CCSD(T)/cc-pVTZ potential energy Surfaces were employed in the calculation. Photoionization induces significant changes in equilibrium structures, which results in long progressions in the v(5), v(4), and v(3) modes. The two isomers show progressions in different modes, which leads to qualitatively distinguishable spectra. The spectra were also calculated in the double harmonic parallel-mode (i.e, neglecting Duschinsky rotation) approximation. Calculating displacements along the normal coordinates of the cation state was found to give a better approximation to the vibrational configuration interaction spectrum; this is due to the effects of Duschinsky rotations on the vibrational wave functions.