Copy a Raster Dataset

GDAL Raster Driver documentation (this is very useful): https://gdal.org/drivers/raster/index.html

from osgeo import gdal

driver = gdal.GetDriverByName('GTiff')

fn_in = r"..\data\input\USGS_one_meter_x64y486_ID_FEMAHQ_2018.tif"
fn_out = r"..\data\output\USGS_one_meter_x64y486_ID_FEMAHQ_2018_copy.tif"

ds_in = gdal.Open(fn_in)
ds_out = driver.CreateCopy(fn_out, ds_in, options=["TILED=YES", "COMPRESS=LZW"])

ds_in = None
ds_out = None