【摘要】：离散点数据可能来自本系统断层解释，也可能从外部导入;由于断层面在X-Y平面上不是单值，断层面不能在X-Y平面插值得到。要用散点拟合投影面，在投影面上部分插值，计算断层面。步骤1:设离散点的最佳投影平面为a* x+b* y+c* z+d=0，用最小二乘法估计投影平面系数a、b、c、d。
离散点数据可能来自本系统断层解释，也可能从外部导入;由于断层面在X-Y平面上不是单值，断层面不能在X-Y平面插值得到。要用散点拟合投影面，在投影面上部分插值，计算断层面。
算法(从离散点创建断层)
设pi(xi，yi，zi)，i=1，2，…，n表示n个离散点，则
步骤1:设离散点的最佳投影平面为a* x+b* y+c* z+d=0，用最小二乘法估计投影平面系数a、b、c、d。
步骤2:求投影平面的法向量与三个坐标轴的夹角。
设c≠0，令z =z'－d/c，则a* x+b* y+c* z' =0
令
，则α、β、γ分别表示投影平面法向量与坐标轴ox、oy、oz的夹角;依次绕x、y、z轴旋转－α、－β、－γ，平面方程变为z =0;综合以上变换的矩阵为
令T=T1(－d/c)*T2(－α)*T3(－β)*T4(－γ) ，用齐坐标表示散点的变换，则
步骤3:在xoy平面根据散点
计算外包矩形，并对矩形进行剖分得到网格G0。
步骤4:使用克里格法(或反距离法)给网格G0上所有网格节点插值。
步骤5:令T'=T4(γ)* T3(β)* T2(α)*T1(d/c)对G0的每个节点gij变换后得G1的节点
T'，G1即所求断层面，如图6-8所示。
图6-8 离散点生成的断层面
免责声明：以上内容源自网络，版权归原作者所有，如有侵犯您的原创版权请告知，我们将尽快删除相关内容。
