MATLAB and Simulink resources for Arduino, LEGO, and Raspberry Pi

Learn moreOpportunities for recent engineering grads.

Apply Today
Asked by Gaëlle on 12 Jul 2013

Hello, I have a matrix filled probability numbers (i.e. ranging from 0 to 1) or NaN when the probability is not computed. I would like to display this matrix as a color table (e.g. using imagesc), in order to have a quick visualisation of the result. The colorbar range is thus set as 0 to 1 since I am interested in probability values. However, I would like the NaN fields to appear with another color than the default "-inf" (here the color of 0 since the down limit for the color is set for value 0), for example gray. How can I do this? Thank you very much, Gaelle

Answer by Kelly Kearney on 12 Jul 2013

Edited by Kelly Kearney on 12 Jul 2013

Accepted answer

Do you need to use imagesc, or are you open to using pcolor? The latter will leave NaNs blank, so they will appear the same color as the axis background.

Example:

data = rand(10); data(data > 0.9) = NaN; [nr,nc] = size(data);

subplot(2,1,1); imagesc(data); subplot(2,1,2); pcolor([data nan(nr,1); nan(1,nc+1)]); shading flat; set(gca, 'ydir', 'reverse');

(The NaN-padding in the pcolor call is so the last row and column are shown, similar to imagesc)

Answer by Gaëlle on 12 Jul 2013

Hi, Thank you for your answer, but this will not work in my case, because I am setting the limit of the colorbar (using caxis([0,1]) ). If I set, as you proposed, the NaN fields as e.g. 255, the color assigned to them will be equivalent to the color assigned to the value 1 (here the upper bound of the colorbar). I would rather code something that assign a specific color (grey) in the case of a NaN field. Is there any way to do this? Gaelle

Evan on 12 Jul 2013

It depends. You always run the risk of setting the pixels in question to a value that is already used if you don't know anything specific about your image data.

Are you just wanting to visualize where your data returns as NaN? If so, I would do one of two things: I would plot my "good" data as grayscale, then turn the NaN pixels to a easily noticeable color like red. Alternatively, if you must have or prefer your data in color, you could plot your image then overlay a scatter plot where your elements are NaN.

Here's an example showing the first way:

A = rand(50); %random data for i = 1:10 %turn some random points into NaNs A(randi(50),randi(50),:) = NaN; end R = A(:,:,1); % turn your data into a "pseudo-gray" rgb image. G = A(:,:,1); B = A(:,:,1); R(isnan(A(:,:,1))) = 1; %turn NaNs red G(isnan(A(:,:,1))) = 0; B(isnan(A(:,:,1))) = 0; A(:,:,1) = R; %combine color slices into A A(:,:,2) = G; A(:,:,3) = B; imagesc(A)

Here's an example for the second:

A = rand(50); %random data for i = 1:10 %turn some random points into NaNs A(randi(50),randi(50),:) = NaN; end imagesc(A) hold on [ii jj] = find(isnan(A)); scatter(ii,jj,'ok','MarkerEdgeColor',[1 1 1],'MarkerFaceColor',[0 0 0],'LineWidth',2,'SizeData',100)

Answer by Gaëlle on 12 Jul 2013

Hi, No, the aim is still to plot the probability values. For more clarity: this matrix M is a square matrix, where each element (i,j) corresponds to a comparison (a probability calculus) between an element i and j, e.g. the element M(1,2) gives the probability value of the comparison between element 1 and 2 (note that in mx case M(1,2)≠M(2,1) ).

As for the NaN, they are here simply because I want to keep the size of the matrix similar across all calculation and display. They appear because some comparison do not arise. For example, if the jth element does not exist, the entire line j and row j will be set as NaN. For example, if I want to compare 4 elements, but the 2nd one does not exist yet, the matrix will look like:

if true % code 0.3 NaN 0.4 0.1 NaN NaN NaN NaN 0.2 NaN 0.2 0.5 0.5 NaN 0.6 0.4 end

I thus want to display this matrix, but setting the NaN values as something specific (e.g. a grey value so that it resembles the color of the figure background). Is this more clear? Gaelle

Answer by Gaëlle on 12 Jul 2013

What you showed me is almost that. I just would have liked to keep the 'jet' map instead of turning every value into a grey scale. But I should manage to find a trick. Thank you, Gaelle

Answer by Gaëlle on 12 Jul 2013

But one thing which is bad using your way: the colorbar scale is not set according to the values of the probabilities. So actually, I cannot use your trick...

## 0 Comments