Python: Batch extract spatial extents of DEM

I have approximately 500 DEM tiles with a polygon index shapefile. I need to get the extents of each tile and output to a table i.e. xT, xB, yT and yB.
Can anyone assist me with this?

kind regards,
Kern