mpdaf change the resolution or pixel position

When I just load/open and save an image in a script with mpdaf, mpdaf change the pixel position when astropy does not.

It might be a header keyword issue but if astropy is doing it, mpdaf should be able too. The main lead is that the shear information is lost during the load, but this is not certain.

I put here a google drive link that contain a short python script to produce the files : open_save_mpdaf_astropy.py the HST image : SPT-CLJ0356-5337_jccd92010_drc.fits and the 2 outputs of the script : astropy_output.fits and mpdaf_output.fits

https://drive.google.com/drive/folders/1JG_n57ZJzoyzTmuxYULVvlFTUwQcKnhx?usp=sharing

I was not able to get a nice simple working example with another/lighter image but it happen with several different images, not just this one.