Image sampling of a continuous region by a discrete grid

I have a an image of an ellipse in the region [-3,3]X[-3,3] which was generated from a closed form function. I need to sample this image with a 256x256 or 257x257 grid, obtaining the values of the region [-3,3]X[-3,3] at each grid point. How to do this?