Hi, I have a Digital Elevation Model on a .tif file. In order to visualize a shaded relief version of it, I'm using this code:
A = geotiffread('dem.tif');info = geotiffinfo('dem.tif');lat = info.CornerCoords.Lat;lon = info.CornerCoords.Lon;RObj = georasterref('RasterSize',size(A),'LatitudeLimit',[min(lat(:)) max(lat(:))],'LongitudeLimit',[min(lon(:)) max(lon(:))]);worldmap(A,RObj);cmap = colormap(gray(256));meshlsrm(A,RObj,[135, 27],cmap);
On the plot that result from meshlsrm, I have to do other elaboration (for example I need to use the canny edge detector, and the the Hough transform). Someone knows a way in order to obtain a gray scale image from meshlsrm?
Best Answer