Is there a python utilities/modules that can mosaic multiple raster with taking mean for overlapping areas?

Question

I have, let\'s say 200 raster images that have overlapping areas, and I need to take the mean values of those overlapping pixels. I want to use GDAL or python based modules that