Ok, what's really neat about this is that, if I am understanding correctly, those green circles are all sources from Gaia, and those are actually present in the Hubble infrared image, so that's really interesting to me. Previously the infrared data had to be substituted with DSS data because all those extra sources were just useless. See:
I ended up getting WCS data from the original FITS files, converting them to AVM data, and then manually adding the AVM headers to that jpg I linked to with some tweaks I made by hand. Not ideal, obviously, but it works in a pinch. This, as I have discovered, has been done in the past with many Hubble images, though it is so tedious that it has seemingly been abandoned for the past few years. So we have websites like
http://astropix.ipac.caltech.edu/ that offer these images for outreach purposes, and they try to include AVM headers, but some of these fantastic Hubble images such as the Horsehead, the anniversary version of the Eagle Nebula, etc. do not have any coordinates.
You've given me several files that I don't know how to utilize, unfortunately. I'm really sorry I don't know how, because I'd like to copy the solution to my Hubble image headers, and see how it compares to the one I did manually. I'm sure the solution is much better than what I did, which was too tedious a process for me to get pixel perfect. I wonder if what you've done can be further automated so that Hubble images from ESA and STScI can finally be given the missing AVM coordinates more easily? I wonder if my own personal collection of Hubble processing could also benefit.
As always, you have done so much more than I expected. And I know I don't have to say it, but it is totally ok if you don't have any more help to offer. I know you must be terribly busy, and I almost wonder if you have clones doing all the work. I feel like I've asked too much lately.
The unicorn horn is great.