Nice post! Instead of finding Mean you could try local maxima. Also there are many ways of overlapping images, like ImageMultiply or Overlay.
tutorial = Import["https://wolfr.am/cEj2VscY"];
heatmap = Import["https://wolfr.am/cEj0rqGj"];
cropT = ImageResize[ImageTake[tutorial, {100, 765}, {1, 665}], 600];
cropH = ImageReflect[ImageResize[ImageTake[heatmap, {10, 550}, {25, 570}], 610], Left -> Right];
ImageMultiply[cropH, cropT]
