python

https://corteva.github.io/rioxarray/stable/getting_started/crs_management.html

import rioxarray

rds = rioxarray.open_rasterio('/nas/Dataset/DEC_scenes/CNU/GEP_20200609_000000_CNU_zoom_20.tif')
rds_4326 = rds.rio.reproject("EPSG:4326")
rds_4326.rio.to_raster("/nas/Dataset/DEC_scenes/CNU/test.tif")

epsg

epsg = int(gdal.Info(input, format='json')['coordinateSystem']['wkt'].rsplit('"EPSG","', 1)[-1].split('"')[0])

gdal

gdalwarp -s_srs epsg:4326 -t_srs epsg:3857 input.tif output.tif