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,
load_and_save_gdal_example.py (2.18 KB)
test.tif (78.7 KB)