方法
算法:取三点 P1、P2、P3(x,y,z),构造向量 v1=P2−P1、v2=P3−P1,法向量 n=v1×v2 并保证朝上(C>0)。倾角 dip=arccos(C/|n|);倾向方位角 dipDir=atan2(A,B)(东、北分量,归算到 0~360°);走向 strike=(dipDir−90) mod 360(与 strike+180 等价)。三点共线或水平面时无确定产状。
算例
例:P1(0,0,100)、P2(50,0,80)、P3(0,40,90) → n=(800,500,2000)、|n|=2211.33;倾角 dip=arccos(2000/2211.33)=25.3°,倾向 dipDir=atan2(800,500)=58.0°,走向 strike=(58.0−90+360)%360=328.0°。产状记为 328.0°∠25.3°(倾向 NE)。若 P1(0,0,50)、P2(100,0,40)、P3(100,100,30) → strike=315.0°、dipDir=45.0°、dip=8.0°,为缓倾近水平层。