title: How do you extract a point spread function from a fits image?
body: <p>I'm new to python coding and I have a .fits image of a star and I want to extract the PSF from it. How would i go about doing this?</p>
  
  <p>After that how would I fit a Gaussian curve to it?</p>