¿Cómo puedo convertir una marca de tiempo de Unix (digamos 1232559922) a una fecha juliana fraccional (2454853.03150)?

Encontré un sitio web (http://aa.usno.navy.mil/data/docs/JulianDate.php ) que realiza un cálculo similar pero necesito hacerlo programáticamente.

Las soluciones pueden estar en C / C ++, python, perl, bash, etc.