{"id":1765,"date":"2016-04-25T12:12:01","date_gmt":"2016-04-25T16:12:01","guid":{"rendered":"https:\/\/blogs.mathworks.com\/steve\/?p=1765"},"modified":"2019-11-01T16:38:19","modified_gmt":"2019-11-01T20:38:19","slug":"clim-caxis-imshow-and-imagesc","status":"publish","type":"post","link":"https:\/\/blogs.mathworks.com\/steve\/2016\/04\/25\/clim-caxis-imshow-and-imagesc\/","title":{"rendered":"CLim, caxis, imshow, and imagesc"},"content":{"rendered":"

In response to \"MATLAB image display - autoscaling values with imshow,\"<\/a> MATLAB Answerer Extraordinaire ImageAnalyst<\/a> posted this comment:<\/p>

A discussion of the relationship and interplay of caxis(), CLim, and the values you can pass in inside the brackets to imshow() or imagesc() might be useful. For example, let's say your values range from 200 to 35,000, and you want all values less than 1000 to be blue and all values more than 29000 to be red. And you want a colorbar with 16 steps - 16 discrete colors. How would one go about that, and using which functions?<\/i><\/p>

Good question! Let's have a go at it, starting with CLim<\/tt>.<\/p>

CLim<\/tt> is a property of an Axes<\/tt><\/a> object.<\/p>

To investigate CLim<\/tt>, start with imagesc<\/tt>, some elevation data, and a color bar.<\/p>

load mt_monadnock.mat<\/span>\r\n\r\nimagesc(Zc)\r\naxis image<\/span>\r\n\r\ncolorbar\r\n<\/pre>\"\" 
The Axes<\/tt> object controls many aspects of the plot, including the axes rulers, the ticks, the tick labels, the grid lines, and much more. The function gca<\/tt> (\"get current axes\") returns the Axes<\/tt> object.<\/p>
ax = gca\r\n<\/pre>
\r\nax = \r\n\r\n  Axes with properties:\r\n\r\n             XLim: [0.5000 3.0425e+03]\r\n             YLim: [0.5000 3.0425e+03]\r\n           XScale: 'linear'\r\n           YScale: 'linear'\r\n    GridLineStyle: '-'\r\n         Position: [0.1168 0.1100 0.6965 0.8150]\r\n            Units: 'normalized'\r\n\r\n  Use GET to show all properties\r\n\r\n<\/pre>
The Axes<\/tt> object has a large number of properties, so by default MATLAB shows you just the most commonly used ones. If you run this code interactively, you would see a clickable \"Show all properties\" link.<\/p>
One of those properties is CLim<\/tt> (\"color limits\"), which you can access directly this way:<\/p>
ax.CLim\r\n<\/pre>
\r\nans =\r\n\r\n  141.5285  963.1366\r\n\r\n<\/pre>
If you look closely at the color bar in the image plot above, you can see the correspondence between it and the CLim<\/tt> values. ax.CLim(1)<\/tt> is the bottom value on the color bar, and ax.Clim(2)<\/tt> is the top value.<\/p>
Where did those values come from, though?<\/p>
They were automatically computed from the range of the data being plotted.<\/p>
min(Zc(:))\r\n<\/pre>
\r\nans =\r\n\r\n  141.5285\r\n\r\n<\/pre>
max(Zc(:))\r\n<\/pre>
\r\nans =\r\n\r\n  963.1366\r\n\r\n<\/pre>
You can set the CLim<\/tt> yourself, though, and that changes the way the color is scaled from the data values. Let's set the color limits to expand the visible details of the lower elevations.<\/p>
ax.CLim = [140 400];\r\n<\/pre>\"\" 
Or maybe you want to examine the upper elevations.<\/p>
ax.CLim = [600 970];\r\n<\/pre>\"\" 
ImageAnalyst mentioned the function caxis<\/tt>. That's just a convenient way to set the color limits. It's one step shorter than getting the Axes<\/tt> using gca<\/tt> and then setting its CLim<\/tt> property.<\/p>
caxis([400 600])\r\n<\/pre>\"\" 
You can also use caxis<\/tt> to quickly get back to automatic computation of color limits.<\/p>
caxis('auto'<\/span>)\r\n<\/pre>\"\" 
Then ImageAnalyst asked about the [low high] syntax for imagesc<\/tt> and imshow<\/tt>. This is just another convenience for setting the color limits.<\/p>
imagesc(Zc,[400 600])\r\naxis image<\/span>\r\ncolorbar\r\n<\/pre>\"\" 
ax = gca;\r\nax.CLim\r\n<\/pre>
\r\nans =\r\n\r\n   400   600\r\n\r\n<\/pre>
The final part of ImageAnalyst's comment concerned the number of colors. What if you only want 16 colors? Well, all of the MATLAB colormap functions take an optional input argument specifying the number of colors to use. So just call the colormap function that you want to use and pass it the desired number of colors.<\/p>
caxis('auto'<\/span>)\r\ncolormap(parula(16))\r\ntitle('Full range, 16 colors'<\/span>)\r\n<\/pre>\"\"