code for matlab for global trimmeed mean for an image

want to use an approach to reduce noise by global trimmed mean i.e eliminating 0's and 255's from image and then replacing the central pixel of 3x3 window by remaining pixels of image. the given formula is M=1/N(sum(i to N)g(i)

N= number of noise free pixel g= is the noise free element in image