Custom Colormaps in Matlab

Matlab colormaps are usefull to display continuous spatial information. More about colormap usage and options can be found under Matlab help file https://www.mathworks.com/help/matlab/ref/colormap.html.

However there are cases where the default colormaps are not sufficient. For example the figure below shows the hydraulic head trend over an 80-year simulation. For the areas with positive values the groundwater table is rising while in the areas with negative values the water table drops. However using the default colormap is not easy to distinguish these areas.

colormaps_ex1

To improve the figure we are going to modify the colormap so that the color at value 0 is white, while the positive values will have a blue gradient from deep blue to white, and the negative values a gradient from white to deep red.

After we have plotted the figure we can obtain information about the color axis as follows

In my example cmax = 7.2 and cmin = -5.4.  We can also set our own custom limits for the colors by setting  ax.CLim = [-5 8]; . This will make the cmin = -5 and cmax = 8. This is very useful if one wants consistent colormaps between different figures with different range values.

Next we’ll create a variable that has values that range between cmin and cmax and identify the last row that is negative. In my case that was the 28th row (id = 28).

Next, for each row of the variable lcol we have to assign a color. The lcol(1) corresponds to cmin therefore we want to assign deep red (130,0,0) while the row id corresponds to 0 values which should have white color (255, 255, 255). The last row corresponds to cmax and we will assign deep blue (0, 0, 130). For the remaining row we will interpolate their values as follows:

custom_map is a 64×3 matrix. The columns correspond to the R, G and B normalized color values.

Next we set our custom map to be our colormap

Finally we display it and assign an explanatory text

colormaps_ex2

In the above map we can easily distinguish between the areas where the water table rises or drops.

The full code for the above figure is given below with a minor change that controls the color resolution:

 

 

 

2 thoughts on “Custom Colormaps in Matlab

    1. These variables are not included in these snippets as they are part of the project we are currently working on.
      You don’t necessarily need the variable study_area. The line plot(study_area.X,study_area.Y,'k', 'LineWidth', 1)
      simply plots the black outline that you can omit.

      You need though the TR variable.
      If you have a list of x,y,v (v-> value on x,y point) points that you want to visualize,
      then you can calculate the TR using TR = delaunay(x,y).
      Have a look at that link for more info
      In my code the x and y are p(:,1) and p(:,2) respectively.

Leave a Reply

Your email address will not be published. Required fields are marked *

This site uses Akismet to reduce spam. Learn how your comment data is processed.