本发明属于cad二次开发及应用领域，尤其是涉及一种煤层底板高程等值线高程属性提取方法。
背景技术：
在电磁法数据处理时，常常要用到地面标高和煤层底板标高做地形校正和深度校正；同时，资料解释时，我们也要用到必要的高程数据来做目的电性层切片解释。所以，地面标高和底板标高是电磁法数据处理和解释中不可或缺的数据体。而地面每个测点高程信息在工程测量中可以很直接且准确的获取，但煤层底板高程数据仅仅局限于地震勘探提供的一张全区底板等值线图，不能一一的为每个测点给出具体底板高程信息。
目前，许多地质生产单位、施工企业都停留在人工读取地质信息进行整理并重新手动录入数据的层面，繁杂重复的脑力体力劳动不仅费时、费力、费材，而且精度低，严重制约了工作效率，导致在后续物探数据解释方面的精力和时间不足。在这个竞争十分激烈的时代，高效高质的完成工程项目是我们单位立足于物探行业长盛不衰和创优争先需要攻克的难关。因此，解决煤层底板数据的读取问题，就等于摆脱电磁法数据处理技术的一个瓶颈。
技术实现要素：
有鉴于此，本发明旨在提出一种煤层底板高程等值线高程属性提取方法，实现物探测点对应煤层底板数据的自动提取和计算，为物探后续数据处理工作提供可靠的基础数据体。
为达到上述目的，本发明的技术方案是这样实现的：
一种煤层底板高程等值线高程属性提取方法，包括如下步骤：
(1)通过cad平台，建立与电磁法工程测点布置相对应的坐标体系；
(2)对cad文件里的等值线标高属性进行高程赋值；
(3)通过建立块图元，在等值线上进行离散点标定，将等值线的高程属性传递到离散块图元上；
(4)提取离散块图元的空间属性；
(5)对提取出来的离散点的高程属性值进行插值拟合，分配到网格的各个节点上；
(6)提取网格化后的节点空间属性，即为与之对应的电磁法测点空间属性。
进一步的，所述步骤(1)中，根据实际测网的线距和点距来设定cad系统中的相对坐标体系，使cad系统中的x与y坐标轴上的距离与每一个工程点的相对点线距离一一对应。
进一步的，所述步骤(2)具体包括通过软件自带的expresstools扩展包小工具中的“flatten”命令将等值线线型转换成多线段，输入多段线勾画的高程等值线属性，得到多线段所代表的高程值。
进一步的，所述步骤(3)中，通过在高程多段线上建立离散属性实现高程值的传递，具体采用定距等分方式在多段线上均匀生成属性块。
进一步的，所述步骤(4)中,提取等分间距生成的属性块的z属性，z属性即为高程多段线的高程值属性。
进一步的，所述步骤(5)中，利用surfer的数据网格功能，将提取出来的离散的属性块的z属性进行插值拟合，分配到网格的各个节点上。
进一步的，所述步骤(6)中，设置网格节点的相对距离，通过网格节点获取的z属性即测点插值拟合后获得高程值。
相对于现有技术，本发明所述的一种煤层底板高程等值线高程属性提取方法具有以下优势：
(1)本发明实现了在cad图中实现物探测点高程等值线数值的读取，解决了人工读取的精度问题，效率问题；
(2)本发明结合cad和surfer系统自带的实现属性的传递、插值等功能，完全摒除了人工在工程点上读取、计算高程值的弊端。
附图说明
构成本发明的一部分的附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：
图1为本发明实施例所述的一种煤层底板高程等值线高程属性提取方法流程图。
具体实施方式
需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。
下面将参考附图并结合实施例来详细说明本发明。
本发明为了提高煤层底板数值读取的效率和精度，有效的解决常规人工读取方法的复杂程度，本发明在cad二次开发的基础上采用批量处理技术，结合surfer软件的网格功能，实现物探测点对应煤层底板数据的自动提取和计算，为物探后续数据处理工作提供可靠的基础数据体。具体方法如图1所示，
1、环境要求
要求安装了autodesk公司开发的cad2007|autocad2007及以上版本，方法中提到的二次开发代码应该放置到cad启动时能自动加载文件的目录下，我们一般将代码保存在程序主目录下的support文件夹中acad2007.lsp二次开发的源程序里。
2、设置坐标
电法勘探的工程测点都是按照一定长度的点线距布设的，根据实际测网的线距和点距来设定cad系统中的相对坐标体系。使cad系统中的x与y坐标轴上的距离与每一个工程点的相对点线距离一一对应。许多测网方向与cad的xy轴方向不一致，将工程测点转正并根据点线距换算成cad独立坐标。
3、线型转换
cad中，只有多段线具备“高程”属性，所以通过软件自带的expresstools扩展包小工具中的“flatten”命令实现等值线线型转换。
4、属性设置
将多段线勾画的高程等值线进行属性输入，在多段线的“高程”属性参数里输入它所代表的高程值。
5、属性传递
通过在高程多段线上建立离散的“属性块”来实现“高程值”的传递。该方法中通过“定距等分”方式来在多段线上均匀生成属性块。
生成带高程属性的“属性块”代码：
6、属性提取
等分间距生成的属性块的“z”属性即为高程多段线的“高程值”属性。通过cad二次开发功能将这些自带“x”、“y”和“z”属性的块，逐一提取属性并保存到可编辑的数据格式文件(surfer可识别的dat文件)里。
属性提取时先选定提取的范围，然后逐一将属性块的x”、“y”和“z”属性提取出来保存到指定路径下的文件里。
7、测点属性计算
利用surfer强大的数据网格功能，将提取出来的离散的属性块的“z”属性进行插值拟合，分配到网格的各个节点上。设置网格节点的相对距离(即工程点的相对点线距)，通过网格节点获取的“z”属性即测点插值拟合后获得高程值。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。