Thanks.
As I said I got it into Stellarium and lined up so this particular instance is "solved"
It does feel like there is still some work for someone to make astropy and PI saved fits work better together.
If you want to know more about what I am doing with Stellarium there are my scripts and documentation in this google drive folder:
https://drive.google.com/drive/folders/1oO4i83cSjbxuZrYDZTyBXRr8bSYQDrvr?usp=sharing
Those scripts are more updated than the ones currently available via Stellarium production version.
but the short version is:
1) Solve an image to determine its rotation and parity
2) Make it positive parity if it's not
3) Rotate the image so north is up (zero degrees rotation) and fill with black
4) Shrink the image so it's long edge is 1024 (or 512, or 256) pixels
5) Pad the short edge it is a multiple of 256 pixels and fill with black
6) Solve the new image, to determine the coordinates of each corner of the image, in the World Coordinate System.
7) Output the meta data in a format suitable for Stellarium's textures.json file
Also I'm uploading a demo video as we speak, and will link it here when ready. It demonstrates using your own or club's images in Stellarium, for public outreach, and shows a small part of what is possible in scripting Stellarium.