library("jpeg")
#針對特定座標以外的距離變黑
dat <-read.csv("D:\\...\\MRT.tsv",header=FALSE,stringsAsFactors=TRUE,sep="\t")
orgpic = readJPEG("D:\\...\\old folder\\test1 00001.jpg")

gazepoint<-function(orgpic,intensity)
{
  newpic=orgpic
  dx=dim(orgpic)[1]
  dy=dim(orgpic)[2]
  
  GPX=dat[2,35] #test1 00001.jpg需要第35行第2列之X座標
  GPY=dat[2,36] #test1 00001.jpg需要第36行第2列之Y座標
  #但test1 00002.jpg則需要第35與36行，第3列之X與Y座標
  
  intensity=max(c(dx,dy))
  D2=intensity^2/64
  
  for(i in 1:dx)
  {
    for(j in 1:dy)
    {
      d=(i-GPX)^2+(j-GPY)^2
      if(d>D2) newpic[i,j,]=newpic[i,j,]*exp(-(d-D2)/D2*intensity)
    }
  }
  return(newpic)
}
writeJPEG(gazepoint(orgpic), target ="D:\\...\\new folder\\test 00001.jpg", quality = 0.95)