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?
NOTE:-
The color displayed is determined by the colormap. You can change it from the default to grayscale with the following command:
colormap gray
Afterward, you can do something like the following to obtain the actual image from the surface object:
image = frame2im(getframe(gcf)); image = rgb2gray(image); figure; imshow(image);
SEE COMPLETE ANSWER CLICK THE LINK
Comments
Post a Comment