Is there code in Astropy or other library (any language) that can locate the centroid of a star in a .fits image to sub-pixel precision?

Is there any code available in a Python library such as Astropy or a library in any other language that can:

1. Take the .fits image of a star as input

2. Locate the centroid of the star on the .fits image to a sub-pixel precision. The sub-pixel precision is necessary.

3. Place the image of the star centroid in a .fits file to within sub-pixel precision.

Again, the sub-pixel precision is what makes this project unique. All software out there that does similar processing (based on what I could find) only works down to the precision of a single pixel.

I have spent weeks reading through astronomy related libraries available in Python and other languages, and I have found code that can obtain star centroids to within a 1 pixel precision. But I have not been able to find any code in any library that can obtain the centroid to within subpixel precsion. Any help offered would be greatly appriciated!