Как найти Вальдо с R?

Вдохновленный этой нитью Как найти Waldo с Mathematica?

Я никогда не делал обработку изображений в R, но, возможно, другие люди, которые хотят поделиться...

спасибо!

Ответ 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)