How to plot a 3D scatter with hyperplanes going through a point?

Question

I mean something like the following two pictures combined:

I need the hyperplanes to go through the point in row 3. I came up with some code, but it doesn't show any planes at all.

Data = [3.95, 13.83, 4.12; 2.77, 15.34, 5.85; 4.41, 14.66, 5.548 ]
x= Data(:,1); 
y= Data(:,2); 
z= Data(:,3);

pointA = [4.4, 14.7, 5.5];
pointB = [4.4, 14.7, 5.5];
pointC =  [4.4, 14.7, 5.5];
pointD = [4.4, 14.7, 5.5];


normal = cross(pointA-pointB, pointA-pointC)
A = normal(1); B = normal(2); C = normal(3); 
D = -dot(normal,pointA);
zLim = [min(z) max(z)];
yLim = [min(y) max(y)];
[Y,Z] = meshgrid(yLim,zLim);
X = (C * Z + B * Y + D)/ (-A);
reOrder = [1 2  4 3];
figure();patch(X(reOrder),Y(reOrder),Z(reOrder),'r');
grid on;
alpha(0.3);
hold on 
plot3(x,x,z, '.', 'markersize', 30);

Any idea how I could fix this?

If you want "_vertical_" and "_horizontal_" planes you don't need to find the normal of the planes (the normals are already the main unit vectors of your cartesian coordinate system ...). Also in your example your `pointA` to `PointD` are all identical, what do you expect when you calculate your `normal` ? — Hoki, May 13 '19 at 14:54

score 0 · Accepted Answer · answered May 23 '19 at 13:17

Here, try this:

function q56099751
data = [3.95, 13.83, 4.12
        2.77, 15.34, 5.85
        4.41, 14.66, 5.548 ];

x = data(:,1); 
y = data(:,2); 
z = data(:,3); 

% Scatter:
figure(); scatter3(x,y,z,'filled'); hold on;

% Planes:
extents = zeros(3,2);
[extents(1,1), extents(1,2)] = bounds(x);
[extents(2,1), extents(2,2)] = bounds(y);
[extents(3,1), extents(3,2)] = bounds(z);
% extents = extents .* [0.9 1.1]; % See what happens when you uncomment this

for ind1 = 1:numel(x)
  planesThroughPoint(x(ind1), y(ind1), z(ind1), extents , [1 1 0]);
end

function planesThroughPoint(ptX, ptY, ptZ, extents, whichPlane)
if whichPlane(1)
% XY Plane
XYZ = [extents(1,1) extents(2,1) ptZ
       extents(1,2) extents(2,1) ptZ
       extents(1,2) extents(2,2) ptZ
       extents(1,1) extents(2,2) ptZ];
patch( XYZ(:,1), XYZ(:,2), XYZ(:,3), 'r', 'FaceAlpha', 0.2);
end

if whichPlane(2)
% XZ Plane
XYZ = [extents(1,1) ptY extents(3,1)
       extents(1,2) ptY extents(3,1)
       extents(1,2) ptY extents(3,2)
       extents(1,1) ptY extents(3,2)];
patch( XYZ(:,1), XYZ(:,2), XYZ(:,3), 'g', 'FaceAlpha', 0.2);
end

if whichPlane(3)
% YZ Plane
XYZ = [ptX extents(2,1) extents(3,1)
       ptX extents(2,2) extents(3,1)
       ptX extents(2,2) extents(3,2)
       ptX extents(2,1) extents(3,2)];
patch( XYZ(:,1), XYZ(:,2), XYZ(:,3), 'b', 'FaceAlpha', 0.2);
end

A few examples of what this code can produce (depending on how you configure it):