CN103995944A - 一种空间异常极值无偏等值图绘制方法 - Google Patents
一种空间异常极值无偏等值图绘制方法 Download PDFInfo
- Publication number
- CN103995944A CN103995944A CN201410256609.5A CN201410256609A CN103995944A CN 103995944 A CN103995944 A CN 103995944A CN 201410256609 A CN201410256609 A CN 201410256609A CN 103995944 A CN103995944 A CN 103995944A
- Authority
- CN
- China
- Prior art keywords
- point
- point set
- extreme
- regular rectangular
- isogram
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Processing Or Creating Images (AREA)
Abstract
本发明提供了一种空间异常极值无偏等值图绘制方法,对原始采样点集A提取极值点,得到极值点点集A’;所述极值点为局部邻域分析范围内,比周围所有点的空间属性值高或低的点;对原始采样点集A进行规则矩形网格化,得到规则矩形网格节点点集B;合并点集A’和点集B,得到点集C;删除点集C中导致病态三角形生成的畸形网格节点,得到点集D;利用点集D建立不规则三角网,绘制等值线图G。利用本发明绘制的等值图空间异常中心与极值点位置吻合,曲线光滑美观,在保证可视化效果的同时,提高了等值线精度。
Description
技术领域
本发明涉及一种空间异常极值无偏等值图绘制方法,属于地理信息系统领域。
背景技术
利用离散形态的采样数据绘制等值图,是指示空间异常的重要手段之一。理想的等值图如图2(b)所示。传统的等值线绘制主要基于规则矩形格网(Grid)和不规则三角网(Triangulated Irregular Networks,TIN)。
Grid数据通常由离散点数据经过内插而得,得益于矩形格网的规则性和均匀性,利用Grid方法绘制的等值线图更加圆滑美观。然而该方法将网格节点作为等值线追踪的基准点,忽略了原始采样点的位置信息,会造成采样点附近数据失真,精确表示空间异常的关键特征时较为困难,从而造成了等值线圈定的异常中心与采样极值点无法套合的成图现象,在空间异常等值图绘制工作中,即表现为异常中心与极值点偏移的现象,如图2(a)所示。
另一方面,为了避免Grid网格化过程中的精度损失,往往采取TIN方法直接对原始采样点进行网格化。基于该方法进行离散点网格化,网格点就是原始采样点本身,减少了精度损失,同时适用于各种数据分布密度和不规则形状区域。然而,作为TIN的基本单元,三角形的几何形状直接决定着TIN应用质量,由于该方法绘制的等值线趋势易受病态三角形影响,所以往往勾绘出的等值线不够圆滑美观。除此之外,不规则三角网法制作的等值线图边界为不规则凸包,因此还需要对数据进行扩边补空处理。所以该方法并不能妥善解决空间异常中心与极值点偏移的问题。
针对上述问题,众多学者提出了多种解决方法。宋敏和张见明等采用矩形网格法绘制等值线,在等值点的追踪过程中加密局部等值点,并且在细分后的节点上重新计算高程值,有效提高网格化精度,但是异常中心偏移的现象依旧存在。胡金虎通过对原始TIN模型进行加密,得到均匀分布的三角网格。王学春和尚继红利用三次Bezier曲线拼接模型实现了完全满足拓扑关系的等值线光滑。胡德鹏和黄晓萍等则提出了基于TIN追踪等值线的改进算法,并建议采用张力样条算法光滑等值线。以上基于TIN制作等值线图的改进方法,提高了制图可视化效果,但是并没有妥善解决等值图边界的扩边补空处理问题。
空间异常中心与采样极值点之间产生偏差的成图现象,会造成制图可视化成果与客观事实不相符,为异常的追踪、查证带来极大的困难,甚至是误导。
发明内容
为了解决现有技术的不足,本发明提供了兼顾制图精度和可视化效果,本发明提出了一种空间异常极值无偏等值图绘制方法,分别通过原始采样点和规则矩形网格点控制等值线图精度和可视化效果。
本发明的一种空间异常极值无偏等值图绘制方法,利用离散的原始采样点集绘制,具体包括以下步骤:
(1)对原始采样点集A提取极值点,得到极值点点集A’;所述极值点为局部邻域分析范围内,比周围所有点的空间属性值高或低的点;
(2)对原始采样点集A进行规则矩形网格化,得到规则矩形网格节点点集B及规则矩形网;
(3)合并点集A’和点集B,得到点集C;
(4)删除点集C中导致病态三角形生成的畸形网格节点,得到点集D;
(5)利用点集D建立不规则三角网,绘制等值线图G,绘制结束。
步骤(1)通过以下步骤提取原始采样点集A的极值点:
(101)对原始采样点集A构建不规则三角网,利用该不规则三角网制作空间异常等值图G1;
(102)根据空间异常等值图G1中的等值线与原始采样点集A中各点的拓扑关系,提取位于最内层等值线所圈定的区域内的采样点,所有采样点组成极值点点集A’。
步骤(4)通过以下步骤删除畸形网格节点:
(401)已知点集C中,来自规则矩形网格节点点集B的点的数量为N,来自极值点点集A’的点的数量为M;
(402)对于来自极值点点集A’的一点am,且m∈[1,M],确定点am在规则矩形网中所在的单元格R,确定单元格R中距离点am最近的点bn,且点bn来自规则矩形网格节点点集B,n∈[1,N],并通过以下公式计算点bn与点am的距离d:
其中,xbn和ybn分别为点bn的x坐标和y坐标,xam,yam分别为点am的x坐标和y坐标;
(403)通过以下公式计算网格R的对角线长度d’:
其中,X和Y分别为规则矩形网格节点点B中任意2个相邻点在x方向上的点距以及在y方向上的点距;
(404)判断是否删除点bn:给定阈值β,若d<β×d′,则删除点bn;否则保留点bn;
(405)重复步骤402到步骤405,直到对对极值点点集A’中的每个点判断一次是否删除距离其最近的网格顶点,最终剩下的所有点构成点集D。
本发明基于其技术方案所具有的有益效果在于:
(1)本发明引入点集C,通过原始采样点标定异常区中心位置,提高等值线精度,防止高值点和异常区中心偏移的现象出现;通过分布均匀的矩形网格节点控制生成的三角形形状,减少病态三角形的生成,绘制的等值线更加圆滑美观;同时规则矩形网格点的加入,填补了原始采样点范围边界到图幅边界之间的空白区域,能够保证等值线图的边界与图幅边界吻合;
(2)本发明对原始采样点集A进行极值点提取得到极值点点集A’,利用了不规则三角网追踪的等值线图,不会产生空间异常中心与极值点偏移的现象,所以用该方法提取的极值点集满足算法需求;
(3)在点集合并之后,如果直接对合并后的新点集C建立不规则三角网模型,可能会因为规则矩形格网节点与原始采样点的距离过近,导致最终的不规则三角网形状难以达到最优化;在此将与原始采样点距离过近的网格单元顶点称为畸形网格点。如图3(a)所示,实心圆点bn为规则矩形网格节点,三角形点am为原始采样点,由于bn与am的距离过近,导致以点bnam为底边的两个三角形顶角过于尖锐,势必会造成等值线局部曲率变化过大,从而影响可视化效果;即使利用计算机图形学方法对等值线进行光滑处理,其本质是以降低等值线的精度为代价,而且处理效果往往不能达到预期;本发明删除畸形网格节点,重新构建不规则三角网,能够消除尖锐三角形,使得不规则三角网形状更佳;
(4)对比采用规则矩形网格法绘制的等值线和采用不规则三角网绘制的等值线,本发明绘制的等值线更加圆滑美观,极少出现局部曲率变化过大的曲线,更加接近手工制图效果;本发明在保证可视化效果的同时,提高了等值图精度。
附图说明
图1是空间异常极值无偏等值图绘制方法流程图。
图2是极值点与异常中心偏移现象对比示意图。
图3为畸形网格点对三角形网构成影响对比图,其中图3(a)为畸形网格点剔除前效果图,图3(b)为畸形网格点提出后效果图。
图4为畸形网格节点剔除步骤流程图。
图5为利用规则矩形网格绘制法、不规则三角网绘制发和本发明提供的空间异常极值无偏等值图绘制方法绘制的Pb元素地球化学图对比示意图;其中图5(a)、图5(b)和图5(c)分别为三种绘制方法绘制的实验区域显示图,图5(d)、图5(e)和图5(f)分别为三种绘制方法绘制的实验区域内异常区域放大显示图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明。
本发明提供了一种空间异常极值无偏等值图绘制方法,利用离散的原始采样点集绘制,结合图1,具体包括以下步骤:
(1)采用以下步骤对原始采样点集A提取极值点:
(101)对原始采样点集A构建不规则三角网,利用该不规则三角网制作空间异常等值图G1;
(102)根据空间异常等值图G1中的等值线与原始采样点集A中各点的拓扑关系,提取位于最内层等值线所圈定的区域内的采样点,所有采样点组成极值点点集A’;
所述极值点为局部邻域分析范围内,比周围所有点的空间属性值高或低的点;
(2)对原始采样点集A进行规则矩形网格化,得到规则矩形网格节点点集B及规则矩形网;
(3)合并点集A’和点集B,得到点集C;
(4)结合图4,采用以下步骤删除点集C中导致病态三角形生成的畸形网格节点:
(401)已知点集C中,来自规则矩形网格节点点集B的点的数量为N,来自极值点点集A’的点的数量为M;
(402)对于来自极值点点集A’的一点am,且m∈[1,M],确定点am在规则矩形网中所在的单元格R,确定单元格R中距离点am最近的点bn,且点bn来自规则矩形网格节点点集B,n∈[1,N],并通过以下公式计算点bn与点am的距离d:
其中,xbn和ybn分别为点bn的x坐标和y坐标,xam,yam分别为点am的x坐标和y坐标;
(403)通过以下公式计算网格R的对角线长度d’:
其中,X和Y分别为规则矩形网格节点点B中任意2个相邻点在x方向上的点距以及在y方向上的点距;
(404)判断是否删除点bn:给定阈值β,若d<β×d′,则删除点bn;否则保留点bn;
(405)重复步骤402到步骤405,直到对对极值点点集A’中的每个点判断一次是否删除距离其最近的网格顶点,最终剩下的所有点构成点集D;
(5)利用点集D建立不规则三角网,所述不规则三角网效果如图3(b)所示;利用不规则三角网绘制等值线图G,绘制结束。
下面利用对比实验说明本方法的效果。
采用冀东某地区1:50000水系沉积物离散采样点数据,采样面积为460km2,采样平均密度为4个点/km2,以该地区地球化学异常等值图绘制为例进行说明。
采用下列3种方法绘制等值线:
(1)采用规则矩形网格法绘制等值线:内插采用距离反比加权(Inverse DistanceWeighting,IDW)算法,网格间距设为0.5km×0.5km,四方向搜索,最少有效搜索点数为5个;
(2)采用不规则三角网绘制等值线:采用Delaunay三角剖分算法建网,保证三角网的形态最优;
(3)采用本发明提供的方法绘制等值线。
以0.1为等对数间隔,利用以上三组网格数据分别绘制该地区的Pb元素地球化学图,结果如图5所示。图5中,图5(a)、图5(b)、图5(c)分别为采用上述三种处理方式绘制的研究区域Pb元素地球化学图,其中,矩形区域表示实验区域内空间异常与极值点偏移现象比较明显的异常区域,分别将异常区域同比放大,显示如图5(d)、图5(e)和图5(f)所示,以便对三种网格化方式制作的等值线图进行整体和细节的对比。
就整体的可视化效果而言,图5(a)和5(c)中绘制的等值线更加圆滑美观,极少出现局部曲率变化过大的曲线,更加接近手工制图效果。而图5(b)中等值线形态不够美观且成图边界不规则,无论是用于科研工作还是工程成果报告提交,利用不规则三角网绘制的地球化学图都无法满足要求。
然而,在利用规则矩形网格制作的地球化学图中,会出现元素浓集中心与极值点位置偏移的成图现象,这种现象是不符合实际情况的。尤其当这种问题出现在地球化学元素异常区域的时候,将会给地化工作者带来极大不便甚至误导。图5(d)中,极值点并未处于元素异常区内,而图5(e)和图5(f)中二者的位置是完全吻合的。
实验表明,应用本发明算法制作的等值线图结合了规则矩形格网和不规则三角网二者的优点,在保证可视化效果的同时,提高了等值图精度。
Claims (3)
1.一种空间异常极值无偏等值图绘制方法,利用离散的原始采样点集绘制,其特征在于具体包括以下步骤:
(1)对原始采样点集A提取极值点,得到极值点点集A’;所述极值点为局部邻域分析范围内,比周围所有点的空间属性值高或低的点;
(2)对原始采样点集A进行规则矩形网格化,得到规则矩形网格节点点集B及规则矩形网;
(3)合并点集A’和点集B,得到点集C;
(4)删除点集C中导致病态三角形生成的畸形网格节点,得到点集D;
(5)利用点集D建立不规则三角网,绘制等值线图G,绘制结束。
2.根据权利要求1所述的空间异常极值无偏等值图绘制方法,其特征在于:步骤(1)通过以下步骤提取原始采样点集A的极值点:
(101)对原始采样点集A构建不规则三角网,利用该不规则三角网制作空间异常等值图G1;
(102)根据空间异常等值图G1中的等值线与原始采样点集A中各点的拓扑关系,提取位于最内层等值线所圈定的区域内的采样点,所有采样点组成极值点点集A’。
3.根据权利要求1所述的空间异常极值无偏等值图绘制方法,其特征在于:步骤(4)通过以下步骤删除畸形网格节点:
(401)已知点集C中,来自规则矩形网格节点点集B的点的数量为N,来自极值点点集A’的点的数量为M;
(402)对于来自极值点点集A’的一点am,且m∈[1,M],确定点am在规则矩形网中所在的单元格R,确定单元格R中距离点am最近的点bn,且点bn来自规则矩形网格节点点集B,n∈[1,N],并通过以下公式计算点bn与点am的距离d:
其中,xbn和ybn分别为点bn的x坐标和y坐标,xam,yam分别为点am的x坐标和y坐标;
(403)通过以下公式计算网格R的对角线长度d’:
其中,X和Y分别为规则矩形网格节点点B中任意2个相邻点在x方向上的点距以及在y方向上的点距;
(404)判断是否删除点bn:给定阈值β,若d<β×d′,则删除点bn;否则保留点bn;
(405)重复步骤402到步骤405,直到对对极值点点集A’中的每个点判断一次是否删除距离其最近的网格顶点,最终剩下的所有点构成点集D。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410256609.5A CN103995944B (zh) | 2014-06-10 | 2014-06-10 | 一种空间异常极值无偏等值图绘制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410256609.5A CN103995944B (zh) | 2014-06-10 | 2014-06-10 | 一种空间异常极值无偏等值图绘制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103995944A true CN103995944A (zh) | 2014-08-20 |
CN103995944B CN103995944B (zh) | 2016-11-09 |
Family
ID=51310108
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410256609.5A Expired - Fee Related CN103995944B (zh) | 2014-06-10 | 2014-06-10 | 一种空间异常极值无偏等值图绘制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103995944B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107451226A (zh) * | 2017-07-19 | 2017-12-08 | 中国地质大学(武汉) | 一种基于gis的网格系统及用其制作综合异常剖析图的方法 |
CN111093153A (zh) * | 2018-10-24 | 2020-05-01 | 北京嘀嘀无限科技发展有限公司 | 区域形状的确定方法、装置、电子设备及储存介质 |
CN111932657A (zh) * | 2020-08-03 | 2020-11-13 | 中国海洋石油集团有限公司 | 一种自动检测与绘制重、磁异常特征线的方法 |
CN112785910A (zh) * | 2019-11-07 | 2021-05-11 | 中国石油天然气集团有限公司 | 大动态范围非线性地球物理等值线图绘制方法及装置 |
WO2024124800A1 (zh) * | 2022-12-14 | 2024-06-20 | 深圳市鸿合创新信息技术有限责任公司 | 一种线条折点的确定方法、装置、设备及存储介质 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3031613B2 (ja) * | 1996-11-12 | 2000-04-10 | 株式会社つくばソフト研究所 | カラー/濃淡画像入力出力装置と入力出力方法 |
CN102982566B (zh) * | 2012-06-19 | 2016-03-30 | 克拉玛依红有软件有限责任公司 | 一种基于最小曲率法插值的含断层的等值线图自动生成方法 |
-
2014
- 2014-06-10 CN CN201410256609.5A patent/CN103995944B/zh not_active Expired - Fee Related
Non-Patent Citations (1)
Title |
---|
姜宇涛: "计算几何的不规则三角网算法研究及在GIS中应用", 《中国优秀博硕士学位论文全文数据库(硕士) 信息科技辑》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107451226A (zh) * | 2017-07-19 | 2017-12-08 | 中国地质大学(武汉) | 一种基于gis的网格系统及用其制作综合异常剖析图的方法 |
CN111093153A (zh) * | 2018-10-24 | 2020-05-01 | 北京嘀嘀无限科技发展有限公司 | 区域形状的确定方法、装置、电子设备及储存介质 |
CN111093153B (zh) * | 2018-10-24 | 2021-04-20 | 滴滴智慧交通科技有限公司 | 区域形状的确定方法、装置、电子设备及储存介质 |
CN112785910A (zh) * | 2019-11-07 | 2021-05-11 | 中国石油天然气集团有限公司 | 大动态范围非线性地球物理等值线图绘制方法及装置 |
CN112785910B (zh) * | 2019-11-07 | 2023-07-25 | 中国石油天然气集团有限公司 | 大动态范围非线性地球物理等值线图绘制方法及装置 |
CN111932657A (zh) * | 2020-08-03 | 2020-11-13 | 中国海洋石油集团有限公司 | 一种自动检测与绘制重、磁异常特征线的方法 |
CN111932657B (zh) * | 2020-08-03 | 2022-12-02 | 中国海洋石油集团有限公司 | 一种自动检测与绘制重、磁异常特征线的方法 |
WO2024124800A1 (zh) * | 2022-12-14 | 2024-06-20 | 深圳市鸿合创新信息技术有限责任公司 | 一种线条折点的确定方法、装置、设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN103995944B (zh) | 2016-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Shi et al. | Automatic generation of road network map from massive GPS, vehicle trajectories | |
CN103020342B (zh) | 一种从地面LiDAR数据中提取建筑物轮廓和角点的方法 | |
CN110276732B (zh) | 一种顾及地形特征线要素的山区点云空洞修复方法 | |
CN103995944A (zh) | 一种空间异常极值无偏等值图绘制方法 | |
CN112699623B (zh) | 基于非结构网格规则化重构技术的高精度热流计算方法 | |
CN104299260A (zh) | 一种基于sift和lbp的点云配准的接触网三维重建方法 | |
CN104573705A (zh) | 一种建筑物激光扫描点云数据的聚类方法 | |
Huang et al. | Automatic generation of road maps from low quality GPS trajectory data via structure learning | |
CN108959598B (zh) | 基于gis的cad数据导入到swmm的方法 | |
CN102609982B (zh) | 空间地质数据非结构化模式的拓扑发现方法 | |
CN105761312A (zh) | 一种微地形表面重建方法 | |
CN106526675B (zh) | 断层空间参数自动提取方法 | |
CN104573333A (zh) | 利用聚类分析进行模型优选的方法 | |
CN102637227B (zh) | 基于最短路径的土地资源评价因子作用域划分方法 | |
CN105869210A (zh) | 三维地质表面模型中的插值数据处理方法 | |
CN105931297A (zh) | 三维地质表面模型中的数据处理方法 | |
CN111523606B (zh) | 一种道路信息的更新方法 | |
CN110348483B (zh) | 基于空间信息与统计学习的空间信息变化检测与分类方法 | |
CN104331389A (zh) | 基于八点法的等值线追踪算法 | |
CN108830930B (zh) | 一种基于张量填充的自适应群智感知系统数据收集方法 | |
CN116664758A (zh) | 一种基于断层面的四面体网格自适应分级细化方法 | |
CN106846481B (zh) | 一种地质剖面图的生成方法 | |
CN113344866B (zh) | 一种点云综合精度评价方法 | |
CN104463924A (zh) | 基于散乱点高程采样数据的数字高程地形模型生成方法 | |
Ma et al. | Complex buildings orientation recognition and description based on vector reconstruction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20161109 Termination date: 20170610 |
|
CF01 | Termination of patent right due to non-payment of annual fee |