这里有一个起点,使用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)
EBImage
包具有与 Mathematica 答案中使用的至少一些函数的对应函数。http://www.bioconductor.org/packages/release/bioc/html/EBImage.html - Richie Cottonadimpro
也可能会很有用。http://cran.r-project.org/web/packages/adimpro/index.html - Richie Cotton