Ответ 1
Вот начало, используя пакет raster. Я не знаю, будет ли у меня время для работы над методом взаимной корреляции, используемым в версии вопроса Mathematica, но локальное стандартное отклонение на красных частях изображения, по-видимому, указывает на Waldo в этом случае...
library(raster)
waldo = stack("/Users/Benjamin/Desktop/DepartmentStore.jpg")
r = waldo[[1]] - waldo[[2]] - waldo[[3]]
r[is.na(r)] = 0
r_mask = Which(r > 0)
r_masked = r * r_mask
focalsd = focal(r_masked, w=3, fun=sd)
plot(focalsd)