![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/177998/image.png)
Counitng the number of scatter points above a surface
3 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Adam Bartle
le 11 Août 2015
Commenté : Octavian Macari
le 20 Mar 2018
Hi All,
I am having some troubles and hoping someone out there much smarter with MATLAB than me can help.
The problem I have i
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/147937/image.jpeg)
s;
- I have plotted some data as a 3D scatter plot.
- I have then plotted a surface which is a ‘limiting condition’ for the data I have.
- I need to count the number of data points that fall above and below the surface.
Is this something that can be done?
I have been reading up on slices and searching these forums for help but cant find anything close.
Any help is appreciated.
The code is;
0 commentaires
Réponse acceptée
Mike Garrity
le 11 Août 2015
You can use interp2 to evaluate the surface at the X & Y coordinates of your scatter data.
Consider the following example:
% Make up some sample data
[sx,sy,sz] = peaks(12);
px = randn(1,100);
py = randn(1,100);
pz = 5*randn(1,100);
% draw the surface and points
h1 = surf(sx,sy,sz);
hold on
h2 = scatter3(px,py,pz,'filled');
% Evaluate surface's function @ point's X & Y
spz = interp2(sx,sy,sz, px,py);
% Build mask of points which are above surface
mask = pz > spz;
% Delete old scatter
delete(h2)
% Replace with two scatters with different colors
h3 = scatter3(px(mask),py(mask),pz(mask),'filled');
h4 = scatter3(px(~mask),py(~mask),pz(~mask),'filled');
h3.MarkerFaceColor = 'red';
h4.MarkerFaceColor = 'blue';
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/177998/image.png)
But you should read this blog post to understand some of the subtleties of how surface does it's interpolation before you go too far down this path.
2 commentaires
Octavian Macari
le 20 Mar 2018
Mike Garrity, is it possible to make a surface pace above for all points?
Plus de réponses (1)
Lukas Bystricky
le 11 Août 2015
Modifié(e) : Lukas Bystricky
le 11 Août 2015
Can your surface be expressed as z = f(x,y)? If so for every scatter point (at position (x0,y0, z0)), you can do:
h = z0 - f(x0,y0)
If h > 0 your point is above the surface.
If you can't express your surface as z = f(x,y), you can use interp2 or something to get an approximation to z at any (x,y) and do the same as above.
3 commentaires
Lukas Bystricky
le 11 Août 2015
Modifié(e) : Lukas Bystricky
le 11 Août 2015
You could try to use scatteredInterpolant to evaluate all your points. The process might look something like this:
F = scatteredInterpolant(x, y, z) % x, y and z can be unordered
Then for each point you could call:
z_surface = F(point_x, point_y);
Although this may take a while to run for 46752 points. If you're planning on doing this several times that might be an issue.
Voir également
Catégories
En savoir plus sur Surface and Mesh Plots dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!