Dear all,
I have tried a 3D smoothing of my data, using "loess", or smooth but I did not
get the right smoothed image.
Suppose that my data is a matrix "myim" 50x50 and at each point (x,y) in the
grid I have an intensity value.
myim<-matrix(rnorm(2500),50,50)
x<-seq(1,50,1)
y<-seq(1,50,1)
grid<-expand.grid(x,y)
myd<-data.frame(grid[,1],grid[,2],as.vector(myim))
sm.myd<-fitted(loess(myd$X3~myd$X1*myd$X2,data=myd))
image(matrix(myd$X3,50,50))
Then what I would like to have is a gaussian smoothed image (matrix) of myim.
I will appreciate every help answer from you.
Catalina
|