发明内容
本发明实施例中提供了一种瞬变电磁探测数据等值线成图与高阻区域识别方法,在保证瞬变电磁法探测数据绘制为等值线图与高阻区域识别速度的同时,提高待填充区域颜色判断的精度。
为了解决上述技术问题,本发明实施例公开了如下技术方案:
一种基于不规则三角网的瞬变电磁探测数据等值线成图与高阻区域识别方法,包括:
用瞬变电磁仪获取待绘制区域若干测点对应多个测道的感应电压数据;
将所述感应电压数据传输给数据处理单元,所述数据处理单元将所述感应电压数据反演成离散的视电阻率数据;
以测点号为空间横坐标、同时以测道号为视深坐标确定所述视电阻率数据的空间分布,生成离散数据点集;
根据所述离散数据点集生成不规则三角网TIN和有不同视电阻率的单等值线集数组Isoline[n];
根据所述不规则三角网TIN,获得TIN'的边界点集ST;
根据所述单等值线集数组,获得等值线端点集SL;
将所述等值线端点集SL插入所述边界点集ST中,合并形成新的点集S;
对所述点集S中的边界点和等值线端点进行空间上的排序,形成有序点集S’;
根据所述单等值线集数组Isoline[n],对应建立存储颜色值的颜色数组Color[n+1],给出指定高阻区域视电阻率最小值D以及与所述视电阻率最小值D对应的颜色ColorS;
根据所述有序点集S’中的所有等值线端点,搜索待填充开区域Ri,并通过视电阻率值比较和颜色查找法确定所述填充开区域Ri的填充色,其中,当所述填充开区域Ri的视电阻率大于或等于D时,则所述填充开区域Ri的填充色为所述颜色ColorS;
判断所述有序点集S’中是否存在封闭等值线;
如果所述有序点集S’中存在封闭等值线,则通过视电阻率值比较和颜色查找法确定各所述封闭等值线所围成封闭区域的填充色,其中,当所述封闭区域的视电阻率大于或等于D时,则所述封闭区域的填充色为所述颜色ColorS;
如果所述有序点集S’中不存在封闭等值线,则成图结束。
优选地,所述根据所述有序点集S’中的所有等值线线端点,搜索待填充开区域Ri,并通过视电阻率值比较和颜色查找法确定所述填充开区域Ri的填充色,其中,当所述填充开区域Ri的视电阻率大于或等于D时,则所述填充开区域Ri的填充色为所述颜色ColorS,包括:
依次序获取所述有序点集S’中新等值线端点;
判断所获取的新等值线端点是否为已获取的重复等值线端点;
如果所获取的新等值线端点为已获取的重复等值线端点,则结束开区域填充;
如果所获取的新等值线端点不是已获取的重复等值线端点,则以所述新等值线端点为起点搜索一个待填充开区域Ri;
获取在所述有序点集S’中、与所述新等值线相邻且排在所述新等值线端点前面的一个有序点,其中,所述有序点为等值线端点或边界点;
确定所述新等值线端点所在单等值线集的索引index;
判断所述新等值线端点的视电阻率是否大于所述有序点的视电阻率;
如果所述新等值线端点的视电阻率大于所述有序点的视电阻率,则用index查找所述颜色数组来确定所述开区域Ri的颜色;
如果所述新等值线端点的视电阻率小于或等于所述有序点的视电阻率,则用index+1查找所述颜色数组来确定所述开区域Ri的颜色;
判断所述开区域Ri对应的视电阻率是否大于或等于D;
如果所述开区域Ri对应的视电阻率大于或等于D,则所述填充开区域Ri的填充色为所述颜色ColorS。
优选地,所述如果所述有序点集S’中存在封闭等值线,则通过视电阻率值比较和颜色查找法确定各所述封闭等值线所围成封闭区域的填充色,其中,当所述封闭区域的视电阻率大于或等于D时,则所述封闭区域的填充色为所述颜色ColorS,包括:
把所有封闭等值线按照在所述空间横坐标中的位置从小到大依次排序,并存入队列Q;
依次序获取所述队列Q中新封闭等值线;
判断所获取的新封闭等值线是否为已获取的重复封闭等值线;
如果所获取的新封闭等值线为已获取的重复封闭等值线,则结束封闭等值线填充;
如果所获取的新封闭等值线不是已获取的重复封闭等值线,则取所述新封闭等值线的线头点所在TIN线段的两个端点T1和T2;
取T1和T2中处于所述新封闭等值线所围成区域内部的点T,确定所述线头点所在单等值线集的索引CIndex;
判断所述线头点的视电阻率是否大于所述点T的视电阻率;
如果所述线头点的视电阻率大于所述点T的视电阻率,则用CIndex查找所述颜色数组来确定所述新封闭等值线所围成区域的颜色;
如果所述线头点的视电阻率小于或等于所述点T的视电阻率,则用CIndex+1查找所述颜色数组来确定所述新封闭等值线所围成区域的颜色;
判断所述新封闭等值线所围成区域对应的视电阻率是否大于或等于D;
如果所述新封闭等值线所围成区域对应的视电阻率大于或等于D,则所述新封闭等值线所围成区域的填充色为所述颜色ColorS。
优选地,所述对所述点集S中的边界点和等值线端点进行空间上的排序,包括:
对所述点集S中的边界点和等值线端点进行空间上的逆时针排序。
优选地,所述根据所述单等值线集数组Isoline[n],对应建立存储颜色值的颜色数组Color[n+1],所述单等值线集数组Isoline[n]和所述颜色数组Color[n+1]的对应关系包括:
所述单等值线集数组中一条单等值线Isoline[i]两侧的区域颜色为Color[i]和Color[i+1]。
优选地,所述依次序获取所述有序点集S’中新等值线端点,包括依逆时针次序获取所述有序点集S’中新等值线端点。
优选地,所述成图结束之后,还包括把成图与高阻区识别结果显示在终端屏幕和/或保存到数据库中。
由以上技术方案可见,本发明实施例提供的一种瞬变电磁探测数据等值线成图与高阻区域识别方法,把不同测点不同时间段的感应电动势数据反演为空间离散视电阻率数据,再利用不规则三角网绘制等值线,建立单等值线集数组、以及与所述单等值线集数组相对应的存储颜色值的颜色数组,通过视电阻率比较与颜色查找方法来填充颜色,提高了瞬变电磁法探测数据绘制为等值线图的精度与速度;在填色过程中根据给定高阻区域视电阻率值下限以及相对应的指定颜色值,凡大于等于该值的区域都置为所述指定颜色值,本发明实施例对于高阻区域的识别给出了简便方法,对于有些高阻区,如空洞区域等能够根据其视电阻率范围进行快速识别。
具体实施方式
为了使本技术领域的人员更好地理解本发明中的技术方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
参见图1,为本发明实施例提供的一种瞬变电磁探测数据等值线成图与高阻区域识别方法的流程示意图,具体包括:
步骤S101:获取待绘制区域若干测点对应多个测道的感应电压数据。
具体的,可以用瞬变电磁仪探测待绘制区域中若干测点对应多个测道的感应电压数据。
步骤S102:将所述感应电压数据传输给数据处理单元,所述数据处理单元将所述感应电压数据反演成离散的视电阻率数据。
具体的,将所述感应电压数据传输给数据处理单元,通过所述数据处理单元中预先加载好的反演软件通过一定的算法将所述感应电压数据反演成离散的视电阻率数据。
步骤S103:以测点号为空间横坐标、同时以测道号为视深坐标确定所述视电阻率数据的空间分布,生成离散数据点集。
步骤S104:根据所述离散数据点集生成不规则三角网TIN和有不同视电阻率的单等值线集数组Isoline[n]。
具体的,将所述离散数据点集根据三角网生长算法生成不规则三角网TIN;根据所述离散数据点集生成有n种不同视电阻率的单等值线集数组Isoline[n],所述单等值线集数组Isoline[n]中的所有单等值线集按照视电阻率的大小升序排列,其中Isoline[i](i>=0,i<n)包含一个具有相同视电阻率的单等值线集SingleLines[mi],mi为第i个单等值线集的单等值线条数。
步骤S105:根据所述不规则三角网TIN,获得TIN'的边界点集ST。
步骤S106:根据所述单等值线集数组,获得等值线端点集SL。
步骤S107:将所述等值线端点集SL插入所述边界点集ST中,合并形成新的点集S。
步骤S108:对所述点集S中的边界点和等值线端点进行空间上的排序,形成有序点集S’。
具体的,对所述点集S中的边界点和等值线端点在步骤S103中坐标系中的位置进行空间上的逆时针排序,形成有序点集S’,当然也可以是空间上的顺时针排序。
步骤S109:根据所述单等值线集数组Isoline[n],对应建立存储颜色值的颜色数组Color[n+1],给出指定高阻区域视电阻率最小值D以及与所述视电阻率最小值D对应的颜色ColorS。
其中,根据所述单等值线集数组Isoline[n],对应建立存储颜色值的颜色数组Color[n+1],所述单等值线集数组Isoline[n]与所述颜色数组Color[n+1]的对应关系表述为Isoline[i]中一条单等值线两侧的区域颜色为Color[i]和Color[i+1]。
步骤S110:根据所述有序点集S’中的所有等值线端点,搜索待填充开区域Ri,并通过视电阻率值比较和颜色查找法确定所述填充开区域Ri的填充色,其中,当所述填充开区域Ri的视电阻率大于或等于D时,则所述填充开区域Ri的填充色为所述颜色ColorS。
如图2所示,具体包括如下步骤:
步骤S201:依次序获取所述有序点集S’中新等值线端点。
具体的,依照逆时针次序获取所述有序点集S’中新等值线端点,在本实施例中,获取的所述新等值线端点,在本实施例例中所述新等值线端点命名为等值线端点A,当然,还可以依照顺时针次序获取所述有序点集S’中新等值线端点。
步骤S202:判断所获取的新等值线端点是否为已获取的重复等值线端点。
如果所获取的所述等值线端点A之前已获取的重复等值线端点,即两次获取到同一个等值线端点,则执行步骤S203、即结束开区域填充。
如果所获取的所述等值线端点A不是已获取的重复等值线端点,则执行步骤S204、即以所述等值线端点A为起点搜索一个待填充开区域Ri。
步骤S205:获取在所述有序点集S’中、与所述新等值线相邻且排在所述新等值线端点前面的一个有序点A',其中,所述有序点为等值线端点或边界点。
步骤S206:确定所述等值线端点A所在单等值线集的索引index。
具体为,确定所述等值线端点A在所述单等值线集数组中具体位于哪一单等值线中。
步骤S207:判断所述等值线端点A的视电阻率是否大于所述有序点A'的视电阻率。
如果所述等值线端点A的视电阻率大于所述有序点A'的视电阻率,则执行步骤S208:用index查找所述颜色数组来确定所述开区域Ri的颜色。
如果所述新等值线端点的视电阻率小于或等于所述有序点的视电阻率,则则执行步骤S209:用index+1查找所述颜色数组来确定所述开区域Ri的颜色。
通过步骤S205至步骤S209中的视电阻率比较法,确定开区域对应的电阻率,然后通过步骤S109中建立的视电阻率和颜色数组的对应关系,确定开区域的填充色,提高了绘图速度,同时因为采用与所获取的等值线端点邻近点的视电阻率比较、以及基于不规则三角网的采样方法,提高了开区域的绘图精度。
步骤S210:判断所述开区域Ri对应的视电阻率是否大于或等于D。
如果所述开区域Ri对应的视电阻率大于或等于指定高阻区域视电阻率最小值D,则执行步骤S211,将所述填充开区域Ri的填充色定位为所述颜色ColorS,这样便可以对开区域中的高阻区置为指定的颜色。
将所有开区域填充完毕后,再进行封闭等值线围成区域的填充。
步骤S111:判断所述有序点集S’中是否存在封闭等值线。
如果所述有序点集S’中不存在封闭等值线,则执行步骤S113,成图结束。
其中,成图结束后,还可以成图识别结果显示在终端屏幕上和/或保存到数据库中。
如果所述有序点集S’中存在封闭等值线,则执行步骤S112:通过视电阻率值比较和颜色查找法确定各所述封闭等值线所围成封闭区域的填充色,其中,当所述封闭区域的视电阻率大于或等于D时,则所述封闭区域的填充色为所述颜色ColorS。
如图3所示,具体可以包括如下步骤:
步骤S301:把所有封闭等值线按照在所述空间横坐标中的位置从小到大依次排序,并存入队列Q。
步骤S302:依次序获取所述队列Q中新封闭等值线。
步骤S303:判断所获取的新封闭等值线是否为已获取的重复封闭等值线。
如果所获取的新封闭等值线为已获取的重复封闭等值线,则执行步骤S304:结束封闭等值线填充,即成图结束。
如果所获取的新封闭等值线不是已获取的重复封闭等值线,则执行步骤S305:取所述新封闭等值线的线头点所在TIN线段的两个端点T1和T2,其中,所述线头点为在所述新封闭等值线中选取的任一点。
步骤S306:取T1和T2中处于所述新封闭等值线所围成区域内部的点、并命名为点T,确定所述线头点所在单等值线集的索引CIndex。
其中,所述确定所述线头点所在单等值线集的索引CIndex,具体为确定所述线头点在所述单等值线集数组中具体位于哪一封闭单等值线中。
步骤S307:判断所述线头点的视电阻率是否大于所述点T的视电阻率。
如果所述线头点的视电阻率大于所述点T的视电阻率,则执行步骤S308:用CIndex查找所述颜色数组来确定所述新封闭等值线所围成区域的颜色。
如果所述线头点的视电阻率小于或等于所述点T的视电阻率,则则执行步骤S309:用CIndex+1查找所述颜色数组来确定所述新封闭等值线所围成区域的颜色。
通过步骤S306至步骤S309中的视电阻率比较法,确定封闭等值线对应的电阻率,然后通过步骤S109中建立的视电阻率和颜色数组的对应关系,确定封闭等值线所围成区域的填充色,提高了绘图速度;同时因为步骤S305和S306中基于不规则三角网的视电阻率比较点的选取方法,提高了绘图精度。
步骤S310:判断所述新封闭等值线所围成区域对应的视电阻率是否大于或等于D。
如果所述新封闭等值线所围成区域对应的视电阻率大于或等于D,则执行步骤S311:所述新封闭等值线所围成区域的填充色为所述颜色ColorS这样便可以对封闭等值线所围封闭区域中的高阻区置为指定的颜色。
为对本发明实施例提供的成图和高阻区识别方法作进一步说明,下面将结合具体实施例说明瞬变电磁探测数据等值线成图与高阻区域识别方法的过程。
如图4所示,瞬变电磁仪1以蓄电池3供电、电压为12V,发射与接收线圈2的匝数均为1匝,形状均为正方形。其边长均为20m。
瞬变电磁仪1将采集的感应电压数据发送给计算机4;数据库5与所述计算机4相连接,用于存储成图和高阻区识别结果。
用图4中的瞬变电磁探测系统对某一指定区域进行探测,离散视电阻率数据的最大视电阻率maxZ预设为616、最小视电阻率minZ预设为340,等值线视电阻率个数n预设为16,颜色表Color[n+1]有17种不同的颜色,高阻区域视电阻率最小值D设为550,指定高阻区域颜色为一个颜色,然后依据上述成图与高阻区识别方法,便可以获取识别出高阻区的等值线填充图。
本发明实施例提供的瞬变电磁探测数据等值线成图方法,把不同测点不同时间段的感应电动势数据反演为空间离散视电阻率数据,再利用不规则三角网绘制等值线,通过视电阻率比较与颜色查找方法填充颜色,提高了瞬变电磁法探测数据绘制为等值线图的精度与速度。
需要说明的是,在本文中,诸如“第一”和“第二”等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上所述仅是本发明的具体实施方式,使本领域技术人员能够理解或实现本发明。对这些实施例的多种修改对本领域的技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。