def get_bounds(raster):
    """
    Return the bounding box in physical units
    """

    ulx, xres, _, uly, _, yres = raster.GetGeoTransform()
    lrx = ulx + (raster.RasterXSize * xres)
    lry = uly + (raster.RasterYSize * yres)

    return [ulx, lry, lrx, uly]