# Changing the colour intensity of a 2D plot

30 views (last 30 days)
John on 25 Dec 2012
Hi,
I'm plotting lat long data on a map like this.
I'm using the plotm() function:
plotm(x,y,'LineStyle','none','Marker','o','MarkerSize',20);
I'm trying to change the colours so that areas with more data are a different colour. I would like the intensity of the colour to change with the amount of data points.
Would anybody know how I could achieve this?
Thanks

Walter Roberson on 26 Dec 2012
plotm() cannot have different colors within any one line (column of data).
There are tricks for using patch() for drawing very narrow lines.
It appears to me that your data is sufficiently dense that you could use scatterm() and have the output look relatively continuous. scatterm() accepts a matrix of colors, one RGB row per point.
John on 17 Jan 2013
Edited: John on 17 Jan 2013
Hello Walter,
I tried your new suggestion. Unfortunately the points were all black, similar to the image posted below on Image Analyst's post. When you say 'adjust the first parameter', do mean, 0:.1:1? What would you suggest I change it to? Apologies for all the questions but I don't really understand the code and I'm not sure on how to proceed.
Thank you

Image Analyst on 26 Dec 2012
Edited: Image Analyst on 13 Jan 2013
If you treat it as an image you could do that, with some work. You could run the thing through a convolution, conv2(). This would give you a higher signal where the roads are more dense. Essentially that is an indexed image. Then you apply a colormap to it, and call ind2rgb to turn it into a color image. Then you use the roads as a mask over the color image so that roads in concentrated areas show up as a different color than roads in sparse areas. Hopefully that described it enough for you to carry out the code. Though I think this representation would be very confusing and distracting to look at.
for k = 1 : length(x)
yourImage(round(y(k)), round(x(k)) = 1;
end
densityImage = conv2(yourImage, ones(5)/25, 'same');
imshow(densityImage, []);
colormap(jet(256));
John on 19 Jan 2013
Hello Walter and Image Analyst,
Just wondering if you have had any further thoughts on this?
Thanks