49 一种实用的等值线型数据网格化方法
上传者:高少琛|上传时间:2015-04-28|密次下载
49 一种实用的等值线型数据网格化方法
第25卷第3期
2001年6月物 探 与 化 探GEOPHYSICAL&GEOCHEMICALEXPLORATIONVol.25,No.3Jun.,2001
一种实用的等值线型数据网格化方法
郭 志 宏
(中国国土资源航空物探遥感中心,北京 100083)
摘要:数据网格化通常包括三大类:测线型数据网格化、等值线型数据网格化和离散点型数据网格化。文中研究的等值线型数据分块存储、网格点八方位搜索插值的网格化方法较好地解决了平面等值线型数字化数据的网格化计算问题,其计算数据量大,实用性强,精度高,计算速度快。
关键词:等值线型数据;网格化;数据分块存储;网格八方位搜索
中图分类号:P631 文献标识码:A 文章编号:1000Ο8918(2001)03Ο0203Ο06在物化探实际工作中经常要用到以前的老资料和从别处搜集来的资料,多是除了已绘制的图件外,为主。为了对这些等值线图件资料重新利用,样取数(左右的间隔),并对这些,尔后方可对平面网格数、精度和可信程度,、效果和可靠性。为此有必要研制一种质量精度高、计算速度快、针对大数据量等值线型数据的网格化方法。
1 等值线型数据网格化方法的研究
物化探实际工作中的网格化通常包括三大类:测线型数据网格化、等值线型数据网格化和离散点型数据网格化,不同型式的数据类型应采用不同的网格化方法才能获得好的效果。对于测线型数据,采用三次样条函数插值的网格化方法效果较好[1];对于离散点型数据,则可采用距离平方反比加权平均的网格化方法[2]或者二元三次样条函数插值的网格化方法[3];对于等值线型数据
内容需要下载文档才能查看,目前较实用的网格化方法很少。针对等值线型数据的特点,模拟人工在等值线图上网格取数的方法过程,我们提出等值线型数据分块存储、网格点八方位搜索插值的网格化方法,其主要思路及过程如下。
1.1 网格图幅条块划分
ΔY,计算出网格点线数N×M;然首先根据用户选择的网格图幅范围及网格点线距ΔX,
后根据网格图幅等值线疏密度适当选择的插值搜索半径R,计算出网格图幅将要划分成的数据存储小方块内的网格点线数为LN×LM:
LN=R/ΔX
IF(LN3XK.LT.RS) LN=R/ΔX+1
LM=R/ΔY
收稿日期:2000Ο06Ο20
?204?物 探 与 化 探25
内容需要下载文档才能查看卷
IF(LM3YK.LT.RS) LM=R/ΔY+1
则数据存储小方块的宽度为XLN×YLM:
ΔXXLN=LN×
ΔYYLM=LM×
根据数据存储小方块的大小可将网格图幅划分成NW×MW个网格小方块:
MW=M/LM
IF(MW3LM.LT.M) MW=M/LM+1
NW=N/LN
IF(NW3LN.LT.N) NW=N/LN+1
如图1所示,由于网格点搜索插值时网格图幅四周必须往外扩一个搜索半径范围,因此网格图幅划分成的数据存储小方块数为(NW+2)×(MW+2)
内容需要下载文档才能查看。
图1 等值线型数据分块存储、网格点八方位搜索插值网格化示意
1.2 数据内插加密及坐标转换
为保证后面网格点八方位插值条带内的搜索点数及插值计算精度,首先应将平面等值线型数字化数据沿每条等值线内插加密为网格图幅上1mm左右间隔(具体间隔可选择)的数据点,并将网格图幅范围内的数据点坐标平移、旋转变换到用户确定的网格图幅坐标系中,以利于后面各步的计算工作。
1.3
内容需要下载文档才能查看等值线型数据分块存储
为了加快网格点搜索插值的速度,其关键问题在于减少搜索原始数据点的数量范围,为此,必须先把经坐标转换后的按一条条等值线顺序存放的一系列(X,Y,
内容需要下载文档才能查看T)数据点,重新排列成按网格图幅上划分的方块顺序分块存储的等值线型数据点文件。等值线型数据点(X,Y,T)应存放在图1所示网格图幅上划分的序号为(IX+1,JY+1)的小方块上:
(1≤IX+1≤NW+2)IX=X/XLN+1
(1≤JY=Y/YLM+1 JY+1≤MW+2)
(1≤JY3≤JY3=MOD(JY,3)+13)
将数据点(X,Y,T)存放在以下几个内存数组中:
MNG(IX+1,JY3)=MNG(IX+1,JY3)+1
XN(MNG(IX+1,JY3),IX+1,JY3)=X
3期郭志宏:
内容需要下载文档才能查看一种实用的等值线型数据网格化方法?205?
YM(MNG(IX+1,JY3),IX+1,JY3)=Y
TMN(MNG(IX+1,JY3),IX+1,JY3)=T
MNG为计算小方块内存放的数据点数的内存数组。由于等值线型数字化数据点通常密度大、数量多,因此一般无法将所有的数据点同时装入计算机内存中进行分块划分。从上面式子中看到,JY3的取数范围在1和3之间,这样沿网格图幅Y方向的共MW+2个条带方块,每次仅将其中的3个条带(每个条带沿X方向有NW+2个小)方块的等值线型数据点(仅占整个网格图幅等值线型数据量的3/[MW+2]左右)装入计算机内存数组中进行分块划分,并输出到外存设备上分块顺序存取的数据文件中,依次完成网格图幅上全部MW+2个条带(NW+2)×(MW+2)个小方块内的等值线型数据点的分块存储,成为新的按网格图幅上划分的小方块顺序分块存储的等值线型数据文件。
1.4 网格点八方位搜索插值
人工在等值线图上网格取数的方法实际上就是简单的线性插值法。人工进行线性插值取数的灵活之处在于,每个网格点插值计算时仅搜索网格点附近的等值线数据值,并且线性插值的方向并不是固定不变的,而是基本选择离网格点最近的2条等值线数值变化梯度最大的方向进行线性内插或外推,。
沿图1网格图幅Y方向从1~MW+1有MW+行1个条带方块网格点的搜索插值计算,在此之前,同步读入相邻3;X方向从1~NW+1有NW+LN×LM个网格点中每个网格点的插值计算时3个条带方块中与该网格小方块相邻的9个小方块内占分块存储的整个网格图幅等值线型数据点的9/[(NW+2)×(MW+2)]左右,从而大大地缩小了等值线型数据点搜索的数量范围,节省了数据搜索时间,仅需常规搜索法所用时间的9/[(NW+2)×(MW+2)]左右,这对通常为大数据量的等值线型数据是相当可观的。根据内插加密后的等值线型数字化数据点的密度及网格图幅等值线疏密程度,选择合适的插值条带宽度,以搜索半径R为插值条带长度,形成以网格点为中心的8个方位的插值条带,从相邻的9个小方块内的等值线数据点中,搜索出位于8个方位插值条带内距网格点最近的等值线数据点各1个。数据点的搜索结果分以下3种情况:
1.如果8个方位4组相对的插值条带没有1组插值条带内存在2个等值线数据点,表明该网格点位于网格图幅等值线空白区或边部,则可给网格点赋假值;
2.如果8个方位4组相对的插值条带至少有1组插值条带内存在数值不相等的2个等值线数据点,则可选择距离最近的1组插值条带内的2个数值不相等的数据点,线性内插来求出网格点处的数值;
3.如果8个方位4组相对的插值条带,各组插值条带内存在的2个等值线数据点数值均相等,表明该网格点位于网格图幅等值线极值封闭圈内或边沿半封闭圈内,这种情况下,应从相邻的9个小方块内的等值线数据点中,再次搜索出位于8个方位插值条带内距网格点最近但数值与已存在的前一个点不相等的另一个等值线数据点各1个,则可选择距离最近的1个方位的插值条带内的2个数值不相等的数据点外推,求出网格点处的数值。
按以上方法完成该条带方块NW+1个网格小方块网格点的八方位搜索插值计算,并将该条带方块网格点的计算结果按顺序存入到外存设备上的分块随机存取的数据文件中。依次完成网格图幅上全部MW+1个条带方块网格点的八方位搜索插值计算,从而获得该原始网
?206?物 探 与 化 探25
内容需要下载文档才能查看卷格图幅区网格数据结果。
1.5 大范围网格图幅分幅计算
为了解决大数据量大范围等值线图幅的
网格化计算问题,可将大范围网格图幅沿X
方向按适当选定的图幅大小分成几个相对较
小的网格图幅,尔后重复前面1至4的步骤,
依次将各小网格图幅经网格化计算获得的网
格数据结果,按顺序存入到外存设备上的分块
随机存取数据文件中,最后转换成统一的大范
围原始网格图幅区网格数据结果。
1.6 原始网格图幅网格数据边部去假处理图2 原始网格图幅区边部去假值示意
由于用户选定的原始网格图幅区可能要比图幅上的有效等值线区域大,所以需要将上面获得的原始网格图幅区网格数据进行边部去
假值处理,去掉如图2虚线所示有效网格图幅
区以外无意义的假值数据,从而得到有效网格
图幅区网格数据结果。
1.7 有效网格图幅区网格数据圆滑处理
一定圆滑处理,圆滑处理,数据,是按等值线型数据分块存储、网格点八方位搜
索插值的网格化方法获得的网格图幅最终的
网格化计算结果。
图3是以上等值线型数据网格化方法程
序流程。
2 等值线型数据网格化方法实例计算
图4是内蒙古某地区的重力等值线图经
数字化仪取数后,采用前面阐述的等值线型数
据分块存储、网格点八方位搜索插值的网格化
方法获得网格数据,然后再绘制成图,图中小
圆点是沿原重力图上等值线条的数字化取数
点。
从该图可以看到,网格化计算后重新绘制
的等值线条与原图的等值线条的数字化取数
点基本重合,说明网格化质量精度高,尤其是图3 等值线型数据网格化方法程序流程较好地解决了等值线极值圈内的判断及插值、
图幅边部等值线条的判断及插值两大难点问题,从而获得了较高质量的网格化图。另外该图原图是1∶20万的60cm×70cm大小的图幅,选择图上网格间距为4mm×4mm,图上插值搜
3期郭志宏:
内容需要下载文档才能查看一种实用的等值线型数据网格化方法?207?索半径为50mm,按等值线型数据分块存储、网格点八方位搜索插值的网格化方法,在普通586微机上网格化仅用了不到6min的时间,其速度是常规全区搜索插值的网格化方法所无法比拟的,后者网格化计算时间近2h
内容需要下载文档才能查看。
图4 内蒙古某地区重力等值线型数据网格化方法计算实例图示
实线为等值线型数据网格化的重力等值线;虚线为原料图等值线条的数字化取数点
3 结论
通过上面的方法阐述及计算实例,对等值线型数据分块存储、网格点八方位搜索插值的网格化方法可以总结出以下几点结论。
1.由于采用了大图幅分幅计算技术及等值线数据分条块存储、网格点分条块计算技巧,从而能解决大幅面、大数据量、大网格的计算问题,使得该方法具有很强的实用性。
2.由于采用了网格点八方位搜索插值方法,保证了较高的网格点插值精度,并较好地解决了等值线极值圈内的判断及插值、图幅边部等值线条的判断及插值两大难点问题,保证了网格化的质量效果。
3.分块存储、分块插值技术的采用,大大缩小了原始数据搜索范围及点数,节省了大量时间,从而大大提高了网格化计算速度。
参考文献:
[1] 郭志宏.双三次样条内插网格化方法软件的研制开发及应用[A].见:熊盛青,唐文周,姚正煦.航空物探遥感论文集 [C].北京:地质出版社,1999.
[2] 郭志宏.中国东部航磁数据处理计算及成果分析[D].长春:长春科技大学,1988.
下载文档
热门试卷
- 2016年四川省内江市中考化学试卷
- 广西钦州市高新区2017届高三11月月考政治试卷
- 浙江省湖州市2016-2017学年高一上学期期中考试政治试卷
- 浙江省湖州市2016-2017学年高二上学期期中考试政治试卷
- 辽宁省铁岭市协作体2017届高三上学期第三次联考政治试卷
- 广西钦州市钦州港区2016-2017学年高二11月月考政治试卷
- 广西钦州市钦州港区2017届高三11月月考政治试卷
- 广西钦州市钦州港区2016-2017学年高一11月月考政治试卷
- 广西钦州市高新区2016-2017学年高二11月月考政治试卷
- 广西钦州市高新区2016-2017学年高一11月月考政治试卷
- 山东省滨州市三校2017届第一学期阶段测试初三英语试题
- 四川省成都七中2017届高三一诊模拟考试文科综合试卷
- 2017届普通高等学校招生全国统一考试模拟试题(附答案)
- 重庆市永川中学高2017级上期12月月考语文试题
- 江西宜春三中2017届高三第一学期第二次月考文科综合试题
- 内蒙古赤峰二中2017届高三上学期第三次月考英语试题
- 2017年六年级(上)数学期末考试卷
- 2017人教版小学英语三年级上期末笔试题
- 江苏省常州西藏民族中学2016-2017学年九年级思想品德第一学期第二次阶段测试试卷
- 重庆市九龙坡区七校2016-2017学年上期八年级素质测查(二)语文学科试题卷
- 江苏省无锡市钱桥中学2016年12月八年级语文阶段性测试卷
- 江苏省无锡市钱桥中学2016-2017学年七年级英语12月阶段检测试卷
- 山东省邹城市第八中学2016-2017学年八年级12月物理第4章试题(无答案)
- 【人教版】河北省2015-2016学年度九年级上期末语文试题卷(附答案)
- 四川省简阳市阳安中学2016年12月高二月考英语试卷
- 四川省成都龙泉中学高三上学期2016年12月月考试题文科综合能力测试
- 安徽省滁州中学2016—2017学年度第一学期12月月考高三英语试卷
- 山东省武城县第二中学2016.12高一年级上学期第二次月考历史试题(必修一第四、五单元)
- 福建省四地六校联考2016-2017学年上学期第三次月考高三化学试卷
- 甘肃省武威第二十三中学2016—2017学年度八年级第一学期12月月考生物试卷
网友关注
- 植物病理学试卷一
- 《园林植物识别》复习资料
- 以就业为导向的高职旅游专业实践教学探讨_孔璎红
- 南京师大-教育综合
- 高职旅游专业实践课程教学模式探讨_刘华
- 课程意义
- 德育论文
- 教育学考卷
- 金属材料学第1章和第2章复习思考题
- 社区考试名词解释(新) (1)
- 国际私法(五)
- 高校非常规水资源优化配置研究_以桂林理工大学雁山校区为例_何昭菊
- 关于博士必修环节中--论文阶段成果学术讨论会的通知
- 中西方教育的异同
- 试论职业导向下高职高专旅游管理专_省略__以成都师范学院旅游管理专业
- 国际私法(七)
- 追月少年-----媒体的价值更在于提倡价值—四辩稿(孙少)
- 教育学考卷
- 荆楚理工学院2015年普通专升本招生简章
- 生物化学期末复习资料 最新
- 中学语文教学法(复习大纲)
- 国际公法第二次任务002
- 植物生理学练习
- 电气测量技术复习题2014-12-26 10.59.12
- 植物组织培养课程论文要求
- 第四章森林与环境
- 硕士学位论文评阅意见
- 拟发展对象的材料格式要求
- 山西农业大学饲料分析方法
- 2014辽宁农业职业技术学院单独招生专业介绍
网友关注视频
- 化学九年级下册全册同步 人教版 第22集 酸和碱的中和反应(一)
- 河南省名校课堂七年级下册英语第一课(2020年2月10日)
- 每天日常投篮练习第一天森哥打卡上脚 Nike PG 2 如何调整运球跳投手感?
- 外研版英语三起6年级下册(14版)Module3 Unit2
- 19 爱护鸟类_第一课时(二等奖)(桂美版二年级下册)_T3763925
- 【部编】人教版语文七年级下册《过松源晨炊漆公店(其五)》优质课教学视频+PPT课件+教案,辽宁省
- 二年级下册数学第二课
- 沪教版八年级下册数学练习册21.4(1)无理方程P18
- 外研版英语七年级下册module3 unit1第二课时
- 二次函数求实际问题中的最值_第一课时(特等奖)(冀教版九年级下册)_T144339
- 苏科版八年级数学下册7.2《统计图的选用》
- 第19课 我喜欢的鸟_第一课时(二等奖)(人美杨永善版二年级下册)_T644386
- 第4章 幂函数、指数函数和对数函数(下)_六 指数方程和对数方程_4.7 简单的指数方程_第一课时(沪教版高一下册)_T1566237
- 沪教版牛津小学英语(深圳用) 四年级下册 Unit 12
- 七年级英语下册 上海牛津版 Unit3
- 沪教版牛津小学英语(深圳用)五年级下册 Unit 1
- 【获奖】科粤版初三九年级化学下册第七章7.3浓稀的表示
- 沪教版八年级下册数学练习册21.3(3)分式方程P17
- 【部编】人教版语文七年级下册《泊秦淮》优质课教学视频+PPT课件+教案,广东省
- 北师大版数学四年级下册第三单元第四节街心广场
- 19 爱护鸟类_第一课时(二等奖)(桂美版二年级下册)_T502436
- 沪教版牛津小学英语(深圳用) 五年级下册 Unit 10
- 8.对剪花样_第一课时(二等奖)(冀美版二年级上册)_T515402
- 外研版英语七年级下册module3 unit2第二课时
- 第8课 对称剪纸_第一课时(二等奖)(沪书画版二年级上册)_T3784187
- 北师大版数学四年级下册3.4包装
- 冀教版小学英语五年级下册lesson2教学视频(2)
- 人教版历史八年级下册第一课《中华人民共和国成立》
- 沪教版八年级下册数学练习册一次函数复习题B组(P11)
- 【部编】人教版语文七年级下册《逢入京使》优质课教学视频+PPT课件+教案,辽宁省
精品推荐
- 2016-2017学年高一语文人教版必修一+模块学业水平检测试题(含答案)
- 广西钦州市高新区2017届高三11月月考政治试卷
- 浙江省湖州市2016-2017学年高一上学期期中考试政治试卷
- 浙江省湖州市2016-2017学年高二上学期期中考试政治试卷
- 辽宁省铁岭市协作体2017届高三上学期第三次联考政治试卷
- 广西钦州市钦州港区2016-2017学年高二11月月考政治试卷
- 广西钦州市钦州港区2017届高三11月月考政治试卷
- 广西钦州市钦州港区2016-2017学年高一11月月考政治试卷
- 广西钦州市高新区2016-2017学年高二11月月考政治试卷
- 广西钦州市高新区2016-2017学年高一11月月考政治试卷
分类导航
- 互联网
- 电脑基础知识
- 计算机软件及应用
- 计算机硬件及网络
- 计算机应用/办公自动化
- .NET
- 数据结构与算法
- Java
- SEO
- C/C++资料
- linux/Unix相关
- 手机开发
- UML理论/建模
- 并行计算/云计算
- 嵌入式开发
- windows相关
- 软件工程
- 管理信息系统
- 开发文档
- 图形图像
- 网络与通信
- 网络信息安全
- 电子支付
- Labview
- matlab
- 网络资源
- Python
- Delphi/Perl
- 评测
- Flash/Flex
- CSS/Script
- 计算机原理
- PHP资料
- 数据挖掘与模式识别
- Web服务
- 数据库
- Visual Basic
- 电子商务
- 服务器
- 搜索引擎优化
- 存储
- 架构
- 行业软件
- 人工智能
- 计算机辅助设计
- 多媒体
- 软件测试
- 计算机硬件与维护
- 网站策划/UE
- 网页设计/UI
- 网吧管理