I am working to calculate and plot velocity vs distance for large GAIA DR3 data sets. Are there any efficient code blocks in Astropy to help with this or do I need to code it step by step? I am working with ra, dec, parallax, pmra, and pmdec as well as their errors. I have been looking for coding communities as this project is my introduction to Astropy. Any suggestions or resources are greatly appreciated.