Save Georeferenced Basemap Image

I am interested in using basemap and gdal to plot a georeferenced image, draw features on top of it, and then resave to the same georeferenced format.

I'm wondering the easiest way to go about this. I made a first attempt (script & file attached), but am getting mixed up with all the coordinates (figure,image,georeferenced).

Thanks for any advice,
Scott

load_and_save_gdal_example.py (2.18 KB)

test.tif (78.7 KB)