How to save output figure GEOTIFF format with appropriate Coordinate Reference System (CRS)?

See e.g. psconvert with the -W+g option, https://docs.generic-mapping-tools.org/dev/psconvert.html.