polarXYZ(ra, dec, 1000./parallax)
Note that this performs distance scaling using a simple inversion of parallax, which is not in general reliable for parallaxes with non-negligable errors. Use at your own risk.
astromXYZ(ra, dec, parallax)
icrsToGal(astromXYZ(ra, dec, parallax))