如何在python gdal中按经度和纬度裁剪dem？

import os
from osgeo import gdal, ogr

xmin, ymin, xmax, ymax = 经度最小值, 纬度最小值, 经度最大值, 纬度最大值

dem_file = 'DEM文件路径'
dataset = gdal.Open(dem_file)

projection = dataset.GetProjection()
geotransform = dataset.GetGeoTransform()

xoff = int((xmin - geotransform[0]) / geotransform[1])
yoff = int((geotransform[3] - ymax) / abs(geotransform[5]))
xsize = int((xmax - xmin) / geotransform[1])
ysize = int((ymax - ymin) / abs(geotransform[5]))

output_file = '输出文件路径'
driver = gdal.GetDriverByName('GTiff')
output_dataset = driver.Create(output_file, xsize, ysize, 1, gdal.GDT_Float32)
output_dataset.SetProjection(projection)
output_dataset.SetGeoTransform((xmin, geotransform[1], 0, ymax, 0, geotransform[5]))

gdal.Warp(output_dataset, dataset, outputBounds=[xmin, ymin, xmax, ymax])

dataset = None
output_dataset = None