![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/177998/image.png)
Counitng the number of scatter points above a surface
조회 수: 3 (최근 30일)
이전 댓글 표시
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
채택된 답변
Mike Garrity
2015년 8월 11일
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
Octavian Macari
2018년 3월 20일
Mike Garrity, is it possible to make a surface pace above for all points?
추가 답변 (1개)
Lukas Bystricky
2015년 8월 11일
편집: Lukas Bystricky
2015년 8월 11일
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
Lukas Bystricky
2015년 8월 11일
편집: Lukas Bystricky
2015년 8월 11일
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.
참고 항목
카테고리
Help Center 및 File Exchange에서 Surface and Mesh Plots에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!