Problem with diplaying surface in area based image matching using cross correlation

I have two epipolar images and my objective is to develop a procedure for automatic point transfer of regular raster points in an epipolar image pair. The grid is defined in one of the images. Based on image correlation the corresponding points should be measured automatically in the second image. The result is a parallax for each grid point which leads to a parallax field for the complete grid. Because parallaxes are proportional to height data the parallax field is equivalent to the height data of the terrain surface. I have written the code to get the correlation coefficients, display them as well as get the parallax values but am not being able to display the parallax values as surface. So I would like to request if there is any code in matlab to display those values in surface. I have tried to use Surf, mesh but am not being able to obtain the desired surface.