matlab测绘求坐标系,MATLAB在测绘80与2000坐标系转换中的应用

城市地理176 MATLAB在测绘80与2000坐标系转换中的应用熊高翔 (成都理工大学，四川　成都　610051) 摘要：2008 年 7 月 1 日起我国正式启用 2000 国家大地坐标系作为国家法定坐标系。对现代测绘工作来说，对 80 西安坐标系和2000国家大地坐标系的转换是必要的。基于MTALAB 平台，能够简单、方便地实现两坐标系统中单点的坐标转换，并针对单一复杂的界面，设计了人性化合理的界面。关键词：坐标转换；MATLAB ；界面设计一、引言针对我国建设需求，需要最新 2000 国家大地坐标系坐标。两坐标系之间的转换即是不同坐标基准间的转换和相同坐标基准间的转换。将两种转换相结合，基于 MATLAB，通过参数解算，空间直角坐标系向空间大地坐标系的转换，和高斯正算的方法实现了程序的设计，并做出了人性化、简洁的界面。二、设计的实现 MATLAB是一款集数学计算、图形处理、符号计算、映像处理、数学建模、实时控制、动态仿真等诸多功能的软件，其强大的数学计算能力和人性化的界面操作深受测绘人员的欢迎。为了实现整个转换，首先进行两坐标系的空间直角坐标的转换，再进行 2000 坐标系中空间直角坐标向大地坐标的转换，最后将大地坐标转换为平面坐标。 (一)不同空间直角坐标系的转换主要通过七参数转换中的布尔莎模型进行转换布尔莎模型的矢量表示形式为 : (2.1) 式中： ( 2.2) 为平移参数，k 为尺度参数。 ( 2.3) (2.4) (2.5) 在相同基准下，空间直角坐标系向空间大地坐标系转换方法为： (2.6) 在使用上式进行空间直角坐标到大地坐标的转换时，因为计算大地纬度B时需要用到大地H，而计算大地高又需要用到大地纬度，因此不能直接计算出大地坐标，需要采用迭代计算的方法，具体计算时，先采用式(2.7)求出大地纬度B的初值，然后利用该初值代入公式(2.8)(2.9) 来求 H，N 的初值，再利用H,N的初值代入公式再次求出B值，再将代入公式H 和 N 如此反复，直至求出的 B,H,N 收敛为止。 (2.7) (2.8) (2.9) 高斯投影坐标正算公式如下： (2.10) 图 2.1 界面成果图结束语：通过本文的研究，利用了 MATLAB 的强大的数学能力，针对复杂的坐标转换计算过程以及目前坐标转换软件功能单一、操作不方便等问题，采用窗口、菜单、控件的操作方式，实现了所见即所得的人性化界面设计。成功完成了 80 西安坐标系单点到 2000 国家大地坐标系的转换，程序的运算功能设计合理，用户使用方便，操作简易。参考文献 [1] 孔祥元 , 梅是义 . 控制测量学 ( 下 )[M]. 武汉 : 武汉大学出版社 ,2002. [2] 牛丽娟 . 测量坐标转换模型研究与转换系统实现 P]. 西安 : 长安大学 ,2010 [3] 姚吉利 , 秦永 .MATLAB 在测量数据处理中应用 [M] 山东理工大学 ,2010 [4] 熊四明 .2000 国家大地坐标系下点位坐标转换方法浅析 [J]. 测绘与空间地理信息 ,2009,32 (5) :155-158 作者简介：熊高翔(1994—)，男，汉族，四川省成都市人，在读研究生，单位：四川省成都市成都理工大学，研究方向：测绘工程。技术工程·TECHNOLOGY AND ENGINEERING

matlab测绘求坐标系,MATLAB在测绘80与2000坐标系转换中的应用相关推荐

matlab求曲线极值程序,matlab函数求极值matlab函数求极值.ppt
matlab函数求极值matlab函数求极值 * * 函数的极值 1.一元函数的极值函数命令:fminbnd 调用格式:[x,feval,exitflag,output]=fminbnd(fun,x ...
matlab 线性规划求最大值,MATLAB求解线性规划(含整数规划和01规划)问题.pdf
MATLAB求解线性规划(含整数规划和01规划)问题 MATLAB 求解线性规划(含整数规划和0-1 规划)问题线性规划是数学规划中的一类最简单规划问题,常见的线性规划是一个有约束的,变量范围为有 ...
matlab trapz求二重积分,matlab求积分（超详细，含int integral integral2/3 quad trapz
matlab求积分 matlab求积分函数工具: int 用法1: 格式: int(fun,x,a,b) 功能: 计算定积分用法2: 格式: int(f,x) 功能: 计算不定积分注: 使用int ...
matlab编程求平均,matlab中的分组平均函数grpstats的用法
群组函数grpstats 前面讨论到之平均值求法,通常应用于整个数组之值,若要应用到比较复杂的分组平均问题,则必须使用不同的函数才能达成.此项指令之格式如下: means = grpstats(X, ...
matlab近似计算求积分,matlab 实验二定积分的近似计算
实验二定积分的近似计算利用牛顿-莱布尼兹公式虽然可以精确地计算定积分的值,但它仅适用于被积函数的原函数能用初等函数表达出来的情形．如果这点办不到或者不容易办到,这就有必要考虑近似计算的方法．在定积分 ...
matlab表达式求值,Matlab进行表达式求值(含参变量代入求值)
(一)前言本文简单介绍如何使用matlab中的eval()函数与subs()函数进行表达式(含参/不含参)的求值. (二)不含参表达式对于一个表达式的求值运算,我们使用eval函数求解,如eval ...
matlab如何求极点,matlab求极点和零点
让每个人平等地提升自我一.实验目的 1.利用 MATLAB 的 fdatool 观察传递函数 H(z)的零极点分布.幅度响应.脉冲响应及阶跃响应,并观察零极点分布与系统稳定性...... 假设为一个 ...
外推法程序matlab,急求用MATLAB用龙格库塔和外推法解一阶微分方程
共回答了21个问题采纳率:81% f=inline('-y+x+1','x','y'); %微分方程的右边项 dx=0.05; %x方向步长 xleft=0; %区域的左边界 xright=10; % ...
matlab怎么求噪声功率,MATLAB求噪声的功率谱密度两种办法的区别
先描述下我的问题,我想求我抽样的信号(横坐标时间,纵坐标电压)的噪声功率谱密度.方法一:由维纳辛钦定理可以知道,我需要先求时域信号的自相关函数,然后做DFT得到功率谱,在除以Fs可以得到功率谱密度方 ...
matlab for求积分,matlab计算定积分的方法—————
该楼层疑似违规已被系统折叠隐藏此楼查看此楼 1计算函数zx+e*(cosy +xsiny)的偏导数. 2计算函数f(x,y)= In(1+x^2+y^2)在点(2.4)的全微分 3.设z=xe^xy ...

matlab测绘求坐标系,MATLAB在测绘80与2000坐标系转换中的应用