使用以下公式计算平均曲率 E=1+(Fx**2/Fz**2) F=Fx*Fy/Fz**2 G=1+(Fy**2/Fz**2) L=(1/(Fz**2*grad_F))*det(([Fxx,Fxz,Fx],[Fxz、Fzz、Fz]、[Fx、Fz,0]) M=(1/(Fz**2*grad_F))*det(([Fxy,Fyz,Fy],[Fxz,Fzz,Fz],[F x,Fz,0]) N=(1/(Fz**2*grad_F))*det(([Fyy,Fyz,Fy],[Fyx,Fzz,Fz],[Fly,Fz,0]) A=([L,M],[M,N]) B=([E,F],[F,G]) B_inv=逆(B) k_values=特征值(点(B_inv,A)) 通过求 B_inv*A(两个矩阵的点积)的特征值,即k1和k2。 最后,计算点的平均曲率: 平均曲率=(k1+k2)/2
# coding:utf