小编pap*_*ppu的帖子

我试图计算由点阵列(x,y,z)给出的曲面的曲率.最初我试图拟合多项式方程z = a + bx + cx ^ 2 + dy + exy + fy ^ 2)然后计算高斯曲率

$ K =\frac {F_ {xx}\cdot F_ {yy} - {F_ {xy}} ^ 2} {(1+ {F_x} ^ 2 + {F_y} ^ 2)^ 2} $

然而,如果表面复杂,则问题是合适的.我发现这个Matlab代码用数字计算曲率.我想知道如何在Python中做同样的事情.

function [K,H,Pmax,Pmin] = surfature(X,Y,Z),
% SURFATURE -  COMPUTE GAUSSIAN AND MEAN CURVATURES OF A SURFACE
%   [K,H] = SURFATURE(X,Y,Z), WHERE X,Y,Z ARE 2D ARRAYS OF POINTS ON THE
%   SURFACE.  K AND H ARE THE GAUSSIAN AND MEAN CURVATURES, RESPECTIVELY. …