File Exchange

image thumbnail

Color coded 3D scatterplot

version (5.91 KB) by Ken Garrard
Produces a 3D scatterplot of XYZ data with varying colors


Updated 25 Dec 2014

View Version History

View License

A 3D scatter plot of a Cartesian data set is drawn. The data points are sorted by color and plot3 is called once for each group of points that map to the same color. This reduces execution time significantly for large data sets. By default the points are colored according to their distance from the XY plane. This can changed with the 'ColorData' property. Additionally the upper and lower limits of the color range can be specified. This is useful for creating a series of plots with the same coloring.
The example plot is produced by the following code.

[x,y,z] = peaks(101);
c = gradient(z);
k = hypot(x,y)<3;
plot3k({x(k) y(k) z(k)}, ...
'Plottype','stem','FontSize',12, ...
'ColorData',c(k),'ColorRange',[-0.5 0.5],'Marker',{'o',2}, ...

Plot3k is based on plot3c by Uli Theune.

Cite As

Ken Garrard (2021). Color coded 3D scatterplot (, MATLAB Central File Exchange. Retrieved .

Comments and Ratings (61)

Pankaj Dey

Thank you for this code. How to remove the marker at the intersection of X and Y axis ?


Jinsoo Park

Adam Collard

How would I used this to plot 4 variables in a timetable? I have diesel consumption, engine load, engine power, and I want to color by rpms.

Heather Baroody

Thank you! this is exactly what I was hoping for!

Strange problem I encountered but not sure if I'm alone or if anyone cares but would have saved me time to know... I have a code where I sorted my x,y,z values according to my z values, e.g.:
[z, ii] = sort(z);
x = x(ii);
y = y(ii);
z = z(ii);
My other plots have no problem with this. But for some strange reason, plot3k only uses the sorted z values but not the sorted x,y so the plotted values were incorrect. After scratching my head for a while, all I had to do was use the plot3k before I sorted (since I no longer need the sort function to highlight my lowest errors which are 'z')

Deepak Tunuguntla

I am using R2017a on Linux. When using this script to plot, the colour bar is displayed as one solid colour instead of varying colours. However, this does not happen on a Mac. Does anyone have a similar issue?

Atif Anwer

Maximilian Kanig

Are there any known problems with this function under MATLAB R2016b (academic student use on Ubuntu 16.04 64bit), and if so, is there a newer version/other option for xyz-plotting?
In my case it refuses to plot the data, but scales, grid and number/position of subplots are OK.

Any advice is appreciated!



Much better speed than scatter, thanks.

Ken Garrard

Before calling plot3k open a new figure using the colormap of your choice.
For example,

Sudharsana Iyengar

I am new in using this. I used this to plot a 3d data in Matlab 2015.
By defualt it is taking the color map as parula. is there a way to change it.
do I have to add the syntax
'color map','jet' (or anything elsewhere)
let me know thanks.
By the way it is in an excellent submission. Ken u rock.

Bryant Svedin


Emre Kopanoglu

Wei Liu

Great job Ken!

Wei Liu

Great job Ken!


Would there be a quick way to plot each individual point by RGB color value? I know scatter3 has this capability but it is extremely slow.


I found the colorbar unnecessary for most of my 3-D arrays, but quite useful for plotting 4-D arrays. I added an element to the p structure called "ColorBar" which is checked before generating the colorbar. Then doing plot3k(...,'ColorBar',true,...) can be used to create the plot with the colorbar.


Hi everyone,

I have a 4D dataset similar and want to use this commands. I am totally new to Matlab. I imported the data from excel and biuld a Matrix:
Points=[g psi microstrain f]:
Each letter has a column of datapoints. When I use the commands above I always get the error:
Undefined function for plot3k for input arguments of type cell.

I tried a lot and dont have a solution.
Can anyone help me?


Hi Ken, I guess i'm just repeating previous comments but this is fantastic compared to scatter3. Any idea why it is so much faster? Is it because it plots multiple smaller series? With scatter3 plotting is a slow but working on the graph (rotating, editing properties, saving to jpeg) really brings my computer to its knees. Thanks again.



Thanks a million Ken! :) I should have thought of that.

PS. I tried downloading the plot3k file again and its still the same version. I have added the necessary line on my end. Thanks!


Ken Garrard

Carlyle, Points with a NaN color value can be excluded by setting the Z value for those points to NaN. For (x,y,z) data with color vector c, 'z(isnan(c))=NaN'. An updated plot3k has been submitted.


Hi Ken, I was just using this code and noticed that it does not ignore NaN values for the colour-coded points. Is there a way to make the code NOT plot the (x,y,z) point if the associated colour-value for that point is NaN?

Ken Garrard

Jackson, Use either 'hold on' or 'axes(h)' before calling plot3k to plot into the current axes or the axes with handle 'h'.

Jackson Richards

Works great. One question, is it possible to place the plot on an existing set of axes?


MUCH quicker than scatter3(). What a lifesaver!

Huaxi Zheng

great! Thanks for sharing, Ken!


This is brilliant. Thanks, Ken



works great!

Ken Garrard

Thanks for the suggestion Val. I'll update the submission.

Val Schmidt

This is a fantastic utility - something that should be incorporated into MATLAB proper, as scatter3 is awful.

I do recommend one change however. Currently, plot3k calls "hold on" prior to plotting anything which means you'll always add to whatever figure window is active. Since this is not the default behavior of MATLAB I find it a hassle.

The fix is to move 'hold on' on line 164 so the plotting loop looks like this:

% plot data points in groups by color map index
s = 1; % index of 1st point in a color group
for k = 1:length(dLix) % loop over each non-empty color group
plot3(L(s:dLix(k),1), ... % call plot3 once for each color group
L(s:dLix(k),2), ...
L(s:dLix(k),3), ...
mark_ch, ...
'MarkerSize',mark_sz, ...
'MarkerEdgeColor',cmap(L(s,4),:), ... % same marker color from cmap
'MarkerFaceColor',cmap(L(s,4),:) ); % for all points in group
s = dLix(k)+1; % next group starts at next point
if k==1
hold on; % add points to one set of axes
hold off;

Alexander Stepanov

Great code. One small improvement:
Change length(L) to size(L,1) on 96 and 141 lines. Otherwise, it's not working with one data point


Navneet Viswan

Is it possible to use this code to scan images with laser stripe in it and plot out the 3D point cloud of the laser traced portion? with depth?

Kindly help me out


After trying and failing to plot my data (3e6 points) using scatter3, this little guy did it for my lickety-split. Is a 4D set, so plotting first 3 dimensions as spatial, and fourth as color. Just had to make sure my nargs matched up via cells, so did {data(:,1) data(:,2) data(:,3)} as first arg to plot3k, and data(:,4) as second for color.

A+ rating from this matlab user.

Joseph Miller

Simple, fast and clean. VERY good job! Thank you!



minor grid does not show up on Z-axis when Zscale is set to log
(i.e. set(gca,'zscale','log');

any idea why this is occuring?

Peng Huang

I would suggest to change "length(L)" to "size(L,1)" since when points are too less (i.e. less than 3), error will happen. Besides that, it is great! Thanks:)


Very good - its way quicker than scatter3. Uses much less memory :)

Alvaro Canivell

Very useful indeed !

Something I was missing in it:

- A label for the color bar, that could be passed in the 'strs' together with the title and the rest of the labels.

If somebody wants that, just tell me, I added that

Thanks again !


very good job! thank you!


Adam Pilchak

excellent file.
For those interested in additional functionality in terms of viewing contours and slices through 3d data, check out:


very nice indeed! any chance of further speed optimization?

Shazux The Great

one that you should never miss!

Abel Brown

Holy Jehosaphat!! worked outta the box! The funktion did exactly what i just spent the last hour trying to do! what a fantastic function for scientifik data! you're my hero!!!!

thanks x 10!

Mike Ross

Exactly what I was looking for. Thank you.

Esben Almkvist

I would add a:
if nargout>0
g = gcf;
That way you get a cleaner function if you don't need the output.

Tim A

Thanks a lot, very useful

Joseph Marks

It "rocks"; it allows a "surf"-like plot without the hassles and inaccuracies of preppoing the dating over a regular grid!
First-rate work.

Eric DeGiuli

Great, does just what I wanted.

Alexander K.

Very useful function

Matteo Tajana

thank you, I've been looking for such a function for ages! ;-)

sylvio berlusconi

your are a genious

Josiah Yoder

Just what I was looking for - works great.

Daniel B

Great function. Work well and is fast.

oh kwoun

MATLAB Release Compatibility
Created with R2008b
Compatible with any release
Platform Compatibility
Windows macOS Linux

Inspired by: plot3c

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!