I've recently been using a lot of rasterio and rasterstats in python which operate orders of magnitude faster on my PC for basic operations such as zonal statistics than my approach in FME:
Clip rasters by polygons, then RasterStatisticsCalculator. Is there a more efficiently way of doing these sorts of operations?