Code covered by the BSD License  

Highlights from
DEM: Shaded relief image plot (digital elevation model)

image thumbnail

DEM: Shaded relief image plot (digital elevation model)

by

 

25 Apr 2012 (Updated )

Pretty shaded relief map as flat image.

landcolor(n)
function y = landcolor(n)
%LANDCOLOR Land colormap
%
%	Author: Francois Beauducel <beauducel@ipgp.fr>
%	$Revision: 1.0.0 $   $Date: 2012/05/17 11:22:44 $

J = [ ...
0.095678 0.53427 0.21682 
0.15785 0.5979 0.23274 
0.21286 0.64673 0.2514 
0.26411 0.68789 0.27268 
0.32959 0.72416 0.31308 
0.39794 0.75695 0.36038 
0.46153 0.7871 0.40624 
0.52108 0.81516 0.45135 
0.57702 0.84152 0.49547 
0.62973 0.86645 0.53891 
0.67946 0.89016 0.58187 
0.72647 0.91282 0.62427 
0.77095 0.93455 0.66619 
0.81306 0.95546 0.70772 
0.85292 0.97563 0.7489 
0.89066 0.99514 0.78976 
0.88379 0.98595 0.77038 
0.86389 0.96758 0.73236 
0.84615 0.94972 0.69623 
0.8303 0.93233 0.66186 
0.81612 0.91536 0.6291 
0.80341 0.8988 0.59784 
0.79201 0.8826 0.56795 
0.78191 0.86676 0.53946 
0.7729 0.85123 0.51224 
0.76479 0.83602 0.48615 
0.75747 0.8211 0.46111 
0.75084 0.80645 0.43704 
0.74506 0.79206 0.41414 
0.73981 0.77792 0.39211 
0.73501 0.76401 0.37089 
0.73068 0.75033 0.35052 
0.72683 0.73685 0.33106 
0.72042 0.72074 0.31228 
0.71032 0.70085 0.29417 
0.69761 0.67821 0.27694 
0.68489 0.65558 0.26026 
0.67235 0.63313 0.24418 
0.65997 0.61082 0.22889 
0.64775 0.58874 0.21406 
0.63568 0.56689 0.19983 
0.62376 0.54527 0.18622 
0.61197 0.52391 0.17299 
0.60033 0.50283 0.16046 
0.58881 0.48203 0.14832 
0.57742 0.46151 0.13667 
0.56616 0.44133 0.12555 
0.55502 0.4214 0.11472 
0.54398 0.4019 0.10456 
0.53306 0.38266 0.094633 
0.52226 0.36382 0.085242 
0.51155 0.3453 0.076179 
0.50095 0.32714 0.067515 
0.49045 0.30938 0.059259 
0.48005 0.29193 0.051294 
0.46973 0.27495 0.043796 
0.45951 0.25823 0.0365 
0.44938 0.24206 0.029715 
0.43934 0.22609 0.023063 
0.42938 0.21074 0.016949 
0.41951 0.19556 0.010917 
0.40971 0.18105 0.0054326 
0.4 0.16667 0 
];

l = length(J);
if nargin < 1
	n = 256;
end
y = interp1(1:l,J,linspace(1,l,n),'*linear');

Contact us