George, was your immediate objective to extract star locations from fits images?
If so then you may want to look at astrometry.net. This is a free platesolving service that, as a byproduct, will give the xy coordinates of all stars in the image.
There is an api which can be accessed using Python but you can also do it interactively here: http://nova.astrometry.net/
Note that the distortion coefficients generated (which seek to adjust for optical distortions, particularly away from the image centre) use the professional system, not the one used by most amateur tools such as Maxim and PinPoint, so if you use astrometry.net platesolve solutions for subsequent processing you’ll either need to use Astropy or platesolve again using PinPoint.
If you do use Astropy then there are good tools in it to easily convert between image xy coords and RA/Dec in various epochs and coordinate frames.