由钻孔矢量点生成3D矢量线和地质体

一、由钻孔点生成3D矢量线
获取点对象钻孔点,对点对象做一定距离的圆形缓冲区生成面状对象,最后将面状对象与所有的栅格图层做一定的拉伸处理,形成立体柱状图。

1、整理钻孔点数据，生成Excel表格并导入ArcGIS中，如图：

2、ArcGIS中编写脚本工具实现3D矢量线的生成，如图：

（1）arcpy脚本构建主要思路：
①、由钻孔矢量点数据进行空间插值（克里金）得到栅格数据，并设置研究范围；
②、将上述得到的栅格数据转点，构建TIN数据（各个岩层的上下表面TIN）；
③、将钻孔矢量点在arcgis中进行Buffer，生成面矢量；
④、编写循环函数实现各层TIN之间的拉伸（3D Analyst 工具—Triangulated Surface—Extrude Between工具实习TIN之间的拉伸），部分代码：
for inTin in inputs:
inTins.append(inTin)

for i in range(len(inTins)-1):
outMP=’Extrude’+str(i+1)+’.shp’
arcpy.ExtrudeBetween_3d(inputs[i],inputs[i+1],inPoly1,outMP)
arcpy.AddMessage(outMP+” Success to execute”)
outs.append(outMP)
（2）AE开发思路：
IConstructMultiPath用于创建,是将非点状几何对象做一定方式的Extrude。它提供了六个方法或方式来生成。

二、地质体的构建：
1、创建底层地质体：
可以编写脚本工具也可以通过模型构建器实习该功能。
脚本;

模型：

2、创建顶层地质体：方法同上
3、高程内插：将DEM数据进行插值得到所需的分辨率

4、高程修正：
由于 岩层顶层高程与DEM之间存在偏差，不能完全与地层面贴合，所以需要用DEM作为约束条件进行岩层顶层的DEM修正

5、生成地质体：
将岩层顶层或底层面在各层TIN之间进行拉伸生成地质体，

结果如图所示：