高分辨率格网模型在地质灾害风险普查数据坐标转换中的应用研究.pdf
《高分辨率格网模型在地质灾害风险普查数据坐标转换中的应用研究.pdf》由会员分享,可在线阅读,更多相关《高分辨率格网模型在地质灾害风险普查数据坐标转换中的应用研究.pdf(7页珍藏版)》请在咨信网上搜索。
1、第 卷 第期宁夏大学学报(自然科学版)年月V o l N o J o u r n a l o fN i n g x i aU n i v e r s i t y(N a t u r a lS c i e n c eE d i t i o n)S e p 文章编号:()高分辨率格网模型在地质灾害风险普查数据坐标转换中的应用研究胡志瑞,范宏涛,马欣,黄玮,赵大江,李鹏,李小琼,(宁夏国土资源调查监测院,宁夏 银川 ;宁夏地质灾害应急中心,宁夏 银川 ;宁夏矿产资源储量评审中心,宁夏 银川 ;自然资源部大地测量数据处理中心,陕西 西安 ;宁夏国土整治修复中心,宁夏 银川 )摘要:以宁夏为例,通过构建
2、 高精度、高分辨率格网模型,计算坐标转换改正量,实现地质灾害数据由 年北京坐标系或 年西安坐标系转换到 年国家大地坐标系的坐标转换基于 个点,分析 年北京坐标系和 年西安坐标系下坐标转换到 年国家大地坐标系的转换精度结果表明,转换后平面点位中误差分别为 ,m高分辨率格网模型可为地质灾害风险普查、数据共享等提供参考,坐标转换方法可为其他区域提供借鉴关键词:地质灾害风险普查;地质灾害数据;国家大地坐标系;坐标转换;格网模型分类号:(中图)P ;P 文献标志码:A收稿日期:基金项目:宁夏重点研发计划基金资助项目(B E G );宁夏自然科学基金资助项目(A A C ;A A C )作者简介:胡志瑞(
3、),男,高级工程师,注册测绘师,主要从事G I S在地质灾害中的应用研究,(电子信箱)n x g e o i n f o c o m通信联系人:黄玮(),男,高级工程师,主要从事地质灾害预警预报方面的研究,(电子信箱)q q c o m第一次全国自然灾害综合风险普查是一次重大的国情国力调查,地质灾害风险普查作为其中一项,主要目标是掌握我国地质灾害风险隐患底数,为深入开展地质灾害防灾减灾工作奠定基础地质灾害数据是进行地质灾害风险普查的基础长期以来,地质 灾 害 数 据 采 用 年 北 京 坐 标 系(简 称B J S )或 年 西 安 坐 标 系(简 称X A S )B J S ,X A S 的
4、测绘基准是参心、二维、低精度、静态的,已无法满足智慧交通、现代物流、防灾减灾等领域对高精度测绘地理信息服务的需求,制约高新技术的应用 C G C S 具有三维、高精度、地心、动态、统一和实用等诸多优势,采用统一坐标系对经济建设、国防建设、社会发展和科学研究十分必要我国从 年开始启用 国家大地坐标系(简称C G C S ),并 经 a的 过 渡 期 进 行 推 广 应用过渡期间,原测绘部门完成了基础测绘成果和基础地理信息数据库的坐标转换,三、四等天文大地网的平差,建立了独立坐标系与 国家大地坐标系 的 联 系,基 本 完 成 了 现 行 大 地 坐 标 系 向C G C S 的过渡近 年 来,学
5、 者 针 对B J S 或X A S 向C G C S 转换进行了大量研究从转换方法看,主要采用参数法;从数据格式看,主要有C A D,A r c G I S平台格式但基于国产M a p G I S平台数据、针对地质灾害领域的坐标转换研究相对较少空间位置是否准确统一、属性信息是否完整翔实,影响地质灾害数据质量,而高质量数据也是构建防灾减灾数据体系、发挥地质灾害数据价值的基础在全国各省地质灾害风险普查工作中,收集到的各类数据需预处理,而数据坐标系统转换是模型重要的环节,关系到各类数据的套合,不仅影响数据的精度,还决定后期空间分析和评价结果的准确性因此有必要建立一套适用于省级范围的坐标转换模型笔者
6、 研 究 地 质 灾 害 数 据 由B J S 或X A S 向宁夏大学学报(自然科学版)第 卷C G C S 的高精度转换方法,以期为其他地区数据的坐标转换提供参考 坐标转换方法 参数法参数法是基于方程的变换方法,根据转换区域大小选取适当的转换模型,通过选取重合点并经过粗差分析后计算转换参数,然后进行坐标转换的一种转换方法对于面积较大区域,一般采用分区转换二维七参数模型、三维七参数模型以及椭球面的多项式拟合算法适用于全国或椭球面的经纬度差在 及以上的省级大地坐标的转换;布尔莎模型适用于空间直角坐标的转换;莫洛金斯基模型适用于省级及以下范围大地坐标的转换;三维四参数模型适用于 以内局部区域空间
7、直角坐标的转换;局部区域平面坐标的转换可用二维四参数模型 格网法格网法的基本思路是,将待转换区域以一定网格划分,根据转换前后坐标系统间公共点的坐标差,通过一定的算法计算各格网结点的坐标差然后用坐标差来内插网格内任意点的值求得改正量这样就实现了坐标系的变换算法主要采用协方差推估、最小曲率、加权平均等数学方法国外应用格网法转换坐标较早,美国采用最小曲率法、日本采用克里格插值法、澳大利亚采用最小二乘法,分别建立了适用于本国的大地坐标转换模型我国学者近年对格网模型进行研究,开发了一些坐标转换软件 目前,B J S ,X A S 与C G C S 坐标数据的相互转换,针对测绘数据的研究较多,方法也相对成
8、熟,但针对地质灾害数据的研究较少采用格网模型的坐标转换优势在于不受制图比例尺和范围限制;相同地理位置的转换改正量保持不变,而且转换改正量是连续的,弥补了分片区转换后相邻分区之间改正量不连续的缺陷格网模型适合连续分幅的地图转换且便于后期图幅接边笔者采用文献 中的方法,对宁夏地质灾害数据的坐标转换进行研究 数据与流程 实验数据地质灾害数据具有专业性强、多源异构、属性信息复杂、数据量大等特点:从比例尺看,各种比例尺并存;从坐标系看,有B J S ,X A S 且以B J S 为主;从数据类型看,有矢量数据、栅格数据和文本数据,矢量数据是在M a p G I S平台产生的,以 w t(点),w l(线
9、),w p(面)格式为主笔者以宁夏地质资料馆馆藏地质灾害详查数据、区域地质图等地质灾害风 险普查所需 基础数据为 研 究 对 象(表)转换流程研究表明,B J S 转换到C G C S ,其图幅X方向平移 m、Y方向平移 m;X A S 转 换 到C G C S ,其 图 幅X方 向 平 移 m、Y方向平移 m各类地形图因坐标系转换造成内部任意两点的长度和方位变动均在制图精度内,可忽略不计,但由坐标系更换引起图廓点坐标的变化超出制图精度范围以万图幅为例,B J S 转换到C G C S ,图廓X方向平移约 mm、Y方向平移约 mm;X A S 转换到C G C S ,图廓X方向平移约 mm、Y
10、方向平移约 mm同时,不同坐标系下转换后数据接边与重合不容忽略数据转换点位的平均精度应小于图上 mm 万图的方里网和图廓线平均位移 mm,而 万图的方里网和图廓线平均位移达 mm,因此大于等于 万地形图中点(图廓点)的位置变动已超过 制图精度,需 重 新 标 记对 于 小 于 万的地质图,由坐标系更换引起图廓点坐标的变动在制图精度以内,可忽略其影响所以对于大于或等于 万的地质图、各类成果数据,必须进行坐标转换表实验数据序号数据名称坐标系统格式比例尺数据量区域地质图(宁夏)B J S M a p G I S 万幅区域地质图(宁夏)X A S M a p G I S 万 幅宁夏地质灾害调查数据库B
11、 J S M a p G I S 万 个县(区、市)宁夏地质灾害详细调查数据库X A S M a p G I S万 个县(区、市)区域地质图空间数据库(宁夏)B J S M a p G I S万 幅第期胡志瑞等:高分辨率格网模型在地质灾害风险普查数据坐标转换中的应用研究通过模型计算宁夏X A S 向C G C S 转换的高精度、高分辨率格网改正量(dB,dL),可得出格网点C G C S 下大地坐标 B J S 转换到X A S 坐标系已有成熟的方法,因此B J S 转C G C S 可先转为X A S ,再通过模型计算宁夏X A S 向C G C S 转换的流程见图XAS 80M(X,Y)X
12、AS 80(B,L)CGCS2000(B,L)G+CGCS2000G&G5#CGCS2000M(X,Y)A1!GdBdL529GG&FE 1E2图X A S 到C G C S 的转换流程 高分辨率格网模型的建立 控制点选择范围为宁夏境内 国家G P S大地控制网点,G P S C级 网 点,国 家 天 文 大 地 网 与 高 精 度G P S 网 联 合 平 差 点,全 国 三、四 等 三 角 网 在C G C S 下平差的三、四等三角网点为保证宁夏省界处的改正量精度,还搜集了相邻省份(甘肃、内蒙古、陕西)与宁夏行政界线相距 k m以内的控制点 控制点分析控制点的精度决定着坐标转换改正量的精度
13、,所以需要对参与计算的控制点的来源、精度、获取年代进行分析计算时优先选择获取年代相对较近的点位控制点按精度选取先后顺序,即与国家三角点重合的 国家G P S大地控制网点,与国家三角点重合的G P SC级网点,国家天文大地网与高精度G P S 网联合平差点,全国三、四等三角网在C G C S 下平差的三、四等三角网点国家三角点布测于半个世纪前,由于人类活动等各种因素的影响,部分点位发生了变化,因而存在粗差点,如果这些粗差点参与改正量计算,会影响结果的精度因此,在选取控制点后,需对选取的各控制点进行试算,直至选取的控制点不含有粗差点首先求取与控制点距离最近的至少个控制点的改正量然后与该点的改正量比
14、较,若该点的改正量绝对值介于周边几个控制点改正量绝对值的平均值的/与倍,则认为该控制点正常,否则认为是异常点对异常点根据分布图形进一步判定,确定后剔除对每个控制点遍历以上步骤后,最终确定用于计算宁夏全区高分辨率格网坐标转换改正量的控制点为 个 改正量计算模型对 个控制点空间上的密度进行分析,通过反距离加权模型获得连续高精度的趋势 反距离权重法假定每个测量点都受一种局部影响,而这种影响随着距离的增大而减小当为任何未测量的位置预测时,反距离权重法预测位置周围的测量值与距离预测位置较远的测量值相比,距离预测位置最近的测量值对预测值的影响较大这种方法为距离预测位置最近的点分配的权重较大,且权重与反距离
15、(数据点与预测位置之间)的P次幂成正比具体计算公式:Pi(BiB)(LiL)c o sB),()dBni(BiPi)niPi,(i,n),dLni(LiPi)niPi(i,n),()B B dB,L L dL,()式中:BiBi,Bi,LiLi,Li,;Bi,Li为控制点大地坐标,B,L为格网点大地坐标;B ,L 为目标C G C S 大地坐标,B,L 为X A S 大地坐标 格网分辨率的确定分辨率设定思路是使每个格网内有个控制点,格网的分辨率越高,内插值越接近于真值,坐标转换精度也就越高由于该过程的运算量大,影响转换速度,但当格网分辨率达到一定值时,转换精度基本趋于稳定根据测绘一、二、三、四
16、等控制点的实际宁夏大学学报(自然科学版)第 卷间隔,将格网分辨率设为 可满足计算要求,即网格大小为 m 改正量的计算 X A S 向C G C S 转换改正量的计算在椭球面上,以X A S 下各个格网点为中心,在一定的搜索区域选择已剔除局部粗差点的一、二、三、四等控制点,依次求取各控制点的大地坐标改正量然后采用反距离加权平均法计算宁夏境内X A S 坐标系到C G C S 的高分辨率、间隔格网改正量dB T ,dL T 搜索区域一般为矩形、圆形和椭圆形研究中,需考虑在各个方向上均等的数据点,因此将搜索区域定义为圆形 B J S 向C G C S 转换改正量的计算B J S 是通过逐级控制、分级
17、分区平差的,没有经过整体平差,因而精度较低因此在不同的区域,坐标精度不同,有的地方相差甚至达到米级X A S 建立后,基本比例尺地形图更新时使用X A S B J S 转换到X A S 已有较成熟、高精度的方法因此,B J S 向C G C S 转换高精度改正量的计算步骤:第一步,计算B J S 向X A S 转换的改正量dB T ,dL T ;第二步,计算X A S 向C G C S 转换的改正量dB T ,dL T ;第 三步,将B J S 向X A S 转换的改正量与高分辨率格网X A S 向C G C S 转换的改正量叠加,得到高精度、高分辨率格网B J S 向C G C S 转换的改
18、正 量dB T ,dL T ,即dB T dB T dB T ,dL T dL T dL T 通过以上方法得到X A S 向C G C S 转换的改正量和B J S 向C G C S 转换的改正量在实际转换中,先通过M a p G I S获得工程文件中的点、线、面文件点文件通过改正量换算为转换后的坐标,线、面文件也是由若干点组成的因此先读取构成线或面文件中各个点的坐标,然后将每个点坐标转换更新,从而实现线、面文件的转换 精度评价坐标转换精度评价有内符合和外符合种评价方法依据计算转换改正量的重合点残差中误差评估坐标转换精度,若点位残差小于倍点位中误差,满足精度要求笔者采用内符合评价方法,计算转换
19、前后坐标残差中误差,对坐标转换精度进行评估 对于n个评估点,坐标转换模型精度计算公式:对于空间直角坐标X残差中误差,MXv v Xn,()对于空间直角坐标Y残差中误差,MYv v Yn,()式中:n为点位个数;v为重合点残差;v 为重合点转换坐标值与重合点已知坐标值的差平面的点位中误差:MpMXMY()X A S 向C G C S 转换的精度评估对宁夏具有X A S 和C G C S 的三角点,选取 个重合点进行精度评估,点位坐标残差频率分布见图图,横坐标是评估点的残差范围,纵坐标是对应范围点的个数图X坐标残差直方图图Y坐标残差直方图由图图可知,X坐标残差基本呈正态分布,残差最大值为 ,最小值
20、为 ,的点在(,的点在(,的点在(,的点在(,Y坐 标 残 差 最 大 值 为 ,最 小 值 为 ,的点在(,的点 在第期胡志瑞等:高分辨率格网模型在地质灾害风险普查数据坐标转换中的应用研究(,的点在(,的点在(,经 过 计 算,X坐 标 的 残 差 中 误 差MX m,Y坐标的残差中误差MY m,平面的点位的中误差Mp m B J S 向C G C S 转换的精度评估从宁夏具有B J S 和C G C S 的三角点中选取 个重合点进行精度评估,点位坐标残差频率分布见图图图X坐标残差直方图图Y坐标残差直方图由图图可知,X坐标残差频率呈正态分布,的点在(,残差最大值为 ,最小值为 ;Y坐标残差最
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 高分辨率 模型 地质灾害 风险 普查 数据 坐标 转换 中的 应用 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。