Accurate infrared spectra of the two hydroxycarbene isomers are computed by diagonalizing the Watson Hamiltonian including up to four mode couplings using full dimensional potential energy and dipole moment surfaces calculated at the CCSD(T)/cc-pVTZ (frozen core) and CCSD6-311G(**) (all electrons correlated) levels, respectively. Anharmonic corrections are found to be very important for these elusive higher-energy isomers of formaldehyde. Both the energy levels and intensities of stretching fundamentals and all overtone transitions are strongly affected by anharmonic couplings between the modes. The results for trans-HCOHHCOD are in excellent agreement with the recently reported IR spectra, which validates our predictions for the cis-isomers.