投影变换的知识.docx
- 文档编号:11813893
- 上传时间:2023-04-02
- 格式:DOCX
- 页数:13
- 大小:32.70KB
投影变换的知识.docx
《投影变换的知识.docx》由会员分享,可在线阅读,更多相关《投影变换的知识.docx(13页珍藏版)》请在冰豆网上搜索。
投影变换的知识
地图投影,是GIS知识体系中重要的组成部分,每个GIS软件都会涉及到这一部分知识,并不是只有MAPGIS软件中才有,MAPGIS软件中的投影变换相比国外的软件更具有针对性,更符合我们国家的国情,比如标准框等。
我这里只是给大家说说我对投影变换的一个理解,讲很多的知识点串起来,不正确的地方,还请大家给予批评指正。
投影变换,我个人理解,就是对“投影”进行“变换”。
只要把握住了这个核心的思想,其他的就不在那么难理解了。
那么下面就要搞清楚两个问题,就是“什么是投影?
”、“为什么要进行投影?
”。
然后再来理解“如何变换”。
那么什么是投影呢?
___4Za7^c.
我们知道,地球是一个近似于梨型的不规则椭球体,而GIS软件所处理的都是二维平面上的地物要素的信息。
所以首先要考的一个问题,就是如果如何将地球表面上的地物展到平面上去。
&I(#Wy_3
最简单的一个方法,或者说是最容易想到的一个方法就是将地球表面沿着某个经线剪开,然后展成平面,即采用这种物理的方法来实现。
可采用物理的方法将地球表面展开成地图平面必然产生裂隙或褶皱,大家可以想象一下,如果把一个足球展成平面的,会是什么结果。
所以这种方法存在着很大的误差和变形,是不行的。
P8*=Ls+-_F
那么我们就可以采用地图投影的方法,就是建立地球表面上的点与地图平面上点之间的一一对应关系,利用数学法则把地球表面上的经纬线网表示到平面上,这样就可以很好的控制变形和误差。
凡是地理信息系统就必然要考虑到地图投影,地图投影的使用保证了空间信息在地域上的联系和完整性,在各类地理信息系统的建立过程中,选择适当的地图投影系统是首先要考虑的问题。
|~X_;1j!
_
所以一句话,投影:
就是建立地球表面上点(Q,λ)和平面上的点(x,y)之间的函数关系式的过程。
O_s_lL~<
这时候就有一个问题要问了,就是随着地图制图理论及科学技术的不断发展,就会有不同的国家,不同的人,提出了不同的数学法则。
这就表示存在着很多的投影方式。
有时候我们需要将不同的投影方式变换成同一种投影方式,或者将不同的投影参数,变换成相同的投影参数,这都需要进行投影变换。
__to=_y#$_
所以一句话,投影变换:
就是将不同的地图投影函数关系式变换的过程。
(}EB2V_9Hh
在MAPGIS中的“投影变换”的定义如下:
将当前地图投影坐标转换为另一种投影坐标,它包括坐标系的转换、不同投影系之间的变换以及同一投影系下不同坐标的变换等多种变换。
~1Wx__=
下面我们就来看看“投影”和“变换”过程中所涉及到的知识点。
地球椭球体@LZ_'Qc}@
f6@fi`U_,
地图投影是指建立地球表面上点(Q,λ)和平面上的点(x,y)之间的函数关系式的过程。
那我们先来看看,如何在地球表面上表示地物要素的空间信息。
只有先将地球表面上的地物要素的空间信息描述好了以后,在将它们通过函数关系式,投影到地图平面上去,这样才可以进行空间分析或者其它的运算。
z__=C'qF`
我们知道:
如果要描述地物要素的空间信息,或者不同地物要素之间的相对空间关系,首先要在地球上建立一个参考系,只有建立了参考系,才能去准确的描述每个地物的坐标等信息。
这涉及到很多地球的形状及椭球体方面的知识。
_?
\[2Po]n
1、地球的形状5=]q+&y\H
地球自然表面是一个起伏不平、十分不规则的表面,有高山、丘陵和平原,又有江河湖海。
地球表面约有71%的面积为海洋所占用,29%的面积是大陆与岛屿。
陆地上最高点与海洋中最深处相差近20公里。
这个高低不平的表面无法用数学公式表达,也无法进行运算。
所以在量测与制图时,必须找一个规则的曲面来代替地球的自然表面。
当海洋静止时,它的自由水面必定与该面上各点的重力方向(铅垂线方向)成正交,我们把这个面叫做水准面。
但水准面有无数多个,其中有一个与静止的平均海水面相重合。
可以设想这个静止的平均海水面穿过大陆和岛屿形成一个闭合的曲面,这就是大地水准面$=?
_1>zv_F
大地水准面所包围的形体,叫大地球体。
由于地球体内部质量分布的不均匀,引起重力方向的变化,导致处处和重力方向成正交的大地水准面成为一个不规则的,仍然是不能用数学表达的曲面。
大地水准面形状虽然十分复杂,但从整体来看,起伏是微小的。
它是一个很接近于绕自转轴(短轴)旋转的椭球体。
所以在测量和制图中就用旋转椭球来代替大地球体,这个旋转球体通常称地球椭球。
_
2、地球的大小___n_/*BK;
关于地球椭球体的大小,由于采用不同的资料推算,椭球体的元素值是不同的。
现将世界各国常用的地球椭球体的数据列表如下:
#lDf8G|ST~
各种地球椭球体模型zzm__Z`Ya
椭球体名称年代长半轴(米)短半轴(米)扁率oS3}xT"U
白塞尔(Bessel)1841637739763560791:
299.15_`R:
p_-"'b
克拉克(Clarke)1880637824963565151:
293.5__-__)__
克拉克(Clarke)1866637820663565841:
295.0(Uk_1_Rt*h
海福特(Hayford)1910637838863569121:
2970_sxZa+G0o
克拉索夫斯基(北京54)1940637824563568631:
298.3m|mY____t
I.U.G.G(西安80)1975637814063567551:
298.252de_[y_z
WGS-841979637813763567591:
298.26]WL_Q_q4q
3、MAPGIS中的椭球体_&+t_!
_LM
在MAPGIS软件中最常用的就是两种椭球体,它们在MAPGIS软件中是以选择“北京54坐标系”或“西安80坐标系”的方式表现出来的。
比如在做标准框时,系统提示我们选择椭球体,这时我们要么选择“北京54”,要么选择“西安80”或者其他。
所以说在MAPGIS中,当提到“北京54坐标系”或“西安80坐标系”时,它们所代表的含义不是大地测量中的大地坐标系,而是指不同的椭球参数,这个一定要搞清楚。
l
_!
Bc0_
下面我们就了解一下我们国家的坐标系。
_~lB_im$o
当前我国采用坐标系主要有:
1954年北京坐标系、1980年西安坐标系、新1954年北京坐标系、WGS84坐标系。
&S\q*H=}i
⑴、1954年北京坐标系'2_`M_T-
该坐标系是通过与原苏联1942年坐标系联测而建立的。
解放后,为了建立我国天文大地网,鉴于当时历史条件,在东北黑龙江边境上同苏联大地网联测,推算出其坐标作为我国天文大地网的起算数据;随后,通过锁网的大地坐标计算,推算出北京点的坐标,并定名为1954年北京坐标系。
因此,1954年北京坐标系是苏联1942年坐标系的延伸,其原点不在北京,而在苏联普尔科沃。
该坐标系采用克拉索夫斯基椭球作为参考椭球,高程系统采用正常高,以1956年黄海平均海水面为基准。
&__V._ps1
该坐标系有两个缺陷:
①、因为它是在东北黑龙江边境上同苏联大地网联测,推算出其坐标作为我国天文大地网的起算数据,所以随着误差的不断累计,到了中国西部以后,测量的数据必须经过严格修正后,才能达到要求。
②、1954年北京坐标系采用克拉索夫斯基椭球作为参考椭球,这一点和其他国家的参考椭球不一致,所以该坐标系的数据必须经过变换后才可以在国际上得到认可。
EBF608nWfW
⑵、1980年西安坐标系x9fNIu
_AQ
1978年4月召开的“全国天文大地网平差会议”上决定建立我国新的坐标系,称为1980年国家大地坐标系。
其大地原点设在西安西北的永乐镇,简称西安原点。
椭球参数选用1975年国际大地测量与地球物理联合会第16界大会的推荐值。
简称IUUG-75地球椭球参数或IAG-75地球椭球。
`_Bvaej?
⑶、新1954年北京坐标系e_CB__(!
Y|
将全国大地网整体平差的结果整体换算到克拉索夫斯基椭球体上,形成一个新的坐标系,称为新1954年北京坐标系。
该坐标系与1980年国家大地坐标系的轴定向基准相同,网的点位精度相同。
j_""I,$t_
⑷、WGS84坐标系&f7_{3BK
在GPS定位中,定位结果属于WGS-84坐标系。
该坐标系是使用了更高精度的VLBL、SLR等成果而建立的。
坐标系原点位于地球质心,Z轴指向BIH1984.0协议地极(CTP)。
不同的投影方式(V_@g?
|LZ
ug_[|'t_R8
前面提到,随着地图制图理论及科学技术的不断发展,就会有不同的国家,不同的人,提出了不同的数学法则。
这就表示存在着很多的投影方式。
下面对不同投影方式做一下归类,详细的资料可以参考有关的书籍。
__x3_>
K{
⑴、按地图投影的构成方法分:
_])}_{GW
a、几何投影:
几何投影源于透视几何学原理,并以几何特征为依据,将地球椭球面Yo;/_7_gG>
的经纬网投影到平面上或投影到可以展成平面的圆柱表面和圆锥表面等几何面上,从而构成:
方位投影、圆柱投影、圆锥投影;#T_Cz$_=t
①、方位投影:
以平面作为辅助投影面,使球体与平面相切或相割,将球体表面上的经纬网投影到平面上构成的一种投影;1B=vr_Gq
②、圆柱投影:
以圆柱表面作为辅助投影面,使球体和圆柱表面相切或相割,将球体表面上的经纬网投影到圆柱表面上,然后再将圆柱表面展成平面而构成的一种投影;_r_~F__T,
③、圆锥投影:
以圆锥表面作为辅助投影面,使球体和圆柱表面相切或相割,将球体表面上的经纬网投影到圆柱表面上,然后再将圆锥表面展成平面而构成的一种投影is_?
`tre\P
据球面和投影面的相对部位不同,上述投投影影有可分为:
正轴投影、横轴投影、斜轴投影;_LtbL[z>]
在圆柱投影中,以正轴和横轴常见;在圆锥投影中以正轴常见;'|&}rL_r:
+
0NVb_7F;k
正、横、斜轴方位投影+m__hYr]Z
(G<"__nnjK
正、横、斜轴圆柱投影\c/jp5=}_
_5p"ng8nR
|;'V":
yDs
正、横、斜轴圆锥投影+e_ Smo^/K`f9 正轴投影经纬线形状.3q__aaXeH b、非几何投影: 通过一系列数学解析法,由几何投影演绎产生了非几何投影,它们并g__2;JJ_} 不借助投影面,而是根据制图的某些特定要求,如考虑制图区域形状等特点,选用合适的投影条件,用数学解析方法,求出投影公式,确定平面和球面之间点与点间的函数关系。 __d_m,_7OQ 据经纬线的形状,可将非几何投影分为伪方位投影、伪圆柱投影、伪圆锥投影、多圆锥投影;(新编地图学P59)_#<0%_Ca ': ! _w%_&\ ①、伪方位投影: 在正轴情况下,伪方位投影的纬线仍投影为同心圆,除中央经线! _jMa%;/_ 投影成直线外,其余经线均投影成对称于中央经线的曲线,且交tS#EqMf_&o 于纬线的共同圆心;D#x_D_-c ②、伪圆柱投影: 在正轴圆柱投影基础上,规定纬线仍为平行直线,除中央经线投影成直线外,其余经线均投影成对称于中央经线的曲线;? _%v_ndI ③、伪圆锥投影: 在圆锥投影基础上,规定纬线仍为同心圆弧,除中央经线仍为直线外,其余经线则投影成对称于中央经线的曲线;_mC_ n,I ④、多圆锥投影: 这是一种假想借助多个圆锥表面与球体相切设计而成的投影,纬线为同轴为同轴圆弧,其圆心位于中央经线上,中央经线为直线,其余经线则投影成对称与中央经线的曲线;B_t[`p\_p@ ⑵、按地图投影的变形性质分: L-_^vlP)Vu ①、等角投影——地球表面上无穷小图形投影后仍保持相似或两微分线段所组成的角度,在投影后仍保持不变,称等角投影,又称正形投影。 在等角投影中变形椭圆为不同大小的圆,它满足条件,极值长度比a=b或经纬线夹角θ=90°和沿经纬度长度比相圆锥等(m=n)。 W)_1)zOD ②、等积投影——地球面上的图形在投影后保持面积不变,称等面积投影。 投影中变形椭圆为面积相等而形状不同的椭圆,这满足条件,面积比P=a×b=mnsinφ=1。 T. _iv ③、任意投影——即不具备以上两种投影的,在任意投影中,如果沿某一主方向的长度比等于1,即a=1或b=1,则这种投影称为等距离投影。 高斯-克吕格投影ilw fa___=#__S 前面对投影方式做了大体的分类后,下面讲解一个具体的投影方式: 高斯-克吕格投影。 _p/: _)_Z__ 高斯-克吕格投影由德国数学家高斯提出,后经克吕格扩充并推导出计算公式,故称为高斯-克吕格投影,简称高斯投影。 为了控制变形,本投影采用分带的思想。 p^&'__C_? 6度带是从0度子午线(在英国的格林尼治天文台附近)起,自西向东每隔经差6为一投影带,全球分为60带,各带的带号用自然序数1,2,3,…60表示。 即以东经0-6为第1带,其中央经线为3E,东经6-12为第2带,其中央经线为9E,其余类推。 n&iWYECz 3度带,是从东经1度30分的经线开始,每隔3度为一带,全球划分为120个投影带。 这样的目的就是为了让6度分带的中央经线全部和3度分带的中央经线重合,3度分带的中央经线只有一半和6度分带的中央经线重合。 9bP_QD{Qb_ 在高斯-克吕格投影上,规定以中央经线为X轴,赤道为Y轴,两轴的交点为坐标原点。 'X_iiR8n@p 如下图所示: sRE$*__^i_ z: _^Kr"=n 我国规定1: 1万、1: 2.5万、1: 5万、1: 10万、1: 25万、1: 50万比例尺地形图,均采用高斯-克吕格投影。 1: 2.5至1: 50万比例尺地形图采用经差6度分带,1: 1万比例尺地形图采用经差3度分带。 MAPGIS中图框的制作A~'2ki5$g #"O9\X_/B_ 由于图框和投影变换紧密相连,故MAPGIS将其放在同一个系统中。 在MAPGIS中生成图框,大家应该用的很多了,这里就不再重复。 只是将其中用到的一些重要的知识点做一下归纳和总结。 ! .+_iA=K{_ 首先搞清楚在MAPGIS大、小比例尺的分界,如下: 它以1: 5000为界Y3-Tg~_/~W 小于或等于1: 5000时,小比例尺,图幅为梯形图幅(在后面地图入库的时候,是选择矩形图幅,还是梯形图幅就看这里),单位为经纬度;jn^i4f_>N_ 大于1: 5000时,大比例尺,图幅为矩形图幅,单位为公里值;#Jy+: _|jJ 根据这个标准,在MAPGIS中我将图框分为一下四类: _B6_x5E__ ①、小比例尺的标准框: 在“系统标准框”菜单下,选择相应的比例命令即可。 9=rYzA? _)+ ②、小比例尺的非标准框: 在“投影变换”菜单下“绘制投影经纬网”命令。 *7ap[YXZ\w ③、大比例尺的标准框: 在“系统标准框”菜单下,选择相应的比例尺命令,在“矩形分幅方法”中选择“正方形”或者“矩形”。 01H3@_0Q6 ④、大比例尺的非标准框: 在“系统标准框”菜单下,选择相应的比例尺命令,在“矩形分幅方法”中选择“任意矩形分幅”。 (Uv{%_q.n6 所以总这里可以看出,小比例尺的标准框和小比例尺的非标准框是通过不同的菜单下不同的命令生成的,而大比例尺的标准框和大比例尺的非标准框则是通过同一个命令生成的,只是“矩形分幅方法”不一致而已。 “北京54坐标系”转“西安80坐标系”__8dFtp3( _ SxO__M@A 首先将MAPGIS平台的工作路径设置为“…..\北京54转西安80”文件夹下。 -09<_;_U_ 下面我们来讲解“北京54坐标系”转“西安80坐标系”的转换方法和步骤。 k_e&c<3m__ 一、数据说明_t_%FS_5 北京54坐标系和西安80坐标系之间的转换其实是两种不同的椭球参数之间的转换,一般而言比较严密的是用七参数布尔莎模型,即X平移,Y平移,Z平移,X旋转(WX),Y旋转(WY),Z旋转(WY),尺度变化(DM)。 若得七参数就需要在一个地区提供3个以上的公共点坐标对(即北京54坐标下x、y、z和西安80坐标系下x、y、z),可以向地方测绘局获取。 fF*a/\h% 二、“北京54坐标系”转“西安80坐标系”的操作步骤_B__\ _ 启动“投影变换模块”,单击“文件”菜单下“打开文件”命令,将演示数据“演示数据_北京54.WT”、“演示数据_北京54.WL”、“演示数据_北京54.WP”打开,如图1所示: c_DAO_5_^ '4O_1_Y0K 图1Z_/_RQq_ 1、单击“投影转换”“菜单下“S坐标系转换”“命令,系统弹出“转换坐标值”“话框,如图2所示: 3J: _! 8Gm_k QTP__1 _u 图2^KlMBKW_yB ⑴、在“输入”一栏中,坐标系设置为“北京54坐标系”,单位设置为“线类单位-米”;e|{6_^g ⑵、在“输出”一栏中,坐标系设置为“西安80坐标系”,单位设置为“线类单位-米”;85H8`YwP_h ⑶、在“转换方法”一栏中,单击“公共点操作求系数”项;ggbew6L$Z_ ⑷、在“输入”一栏中,输入北京54坐标系下一个公共点的(x、y、z),如图2所示;~r_Ko_5#D ⑸、在“输出”一栏中,输入西安80坐标系下对应的公共点的(x、y、z),如图2所示;#_Tce__q_5 ⑹、在窗口右下角,单击“输入公共点”按钮,右边的数字变为1,表示输入了一个公共点对,如图2所示;.@xwl}o$OL ⑺、依照相同的方法,再输入另外的2个公共点对;l_RentNg0b ⑻、在“转换方法”一栏中,单击“七参数布尔莎模型”项,将右边的转换系数项激活;IY: iGn8R ⑼、单击“求转换系数”菜单下“求转换系数”命令,系统根据输入的3个公共点对坐标自动计算出7个参数,如图3所示,将其记录下来;: ~_`E_@`/ 然后单击“确定”按钮;L$ Ss]A_r= 图3)5e_}_Id_ 2、单击“投影转换”菜单下“编辑坐标转换参数”命令,系统弹出“不同地理坐标系转换参数设置”对话框,如图4所示;___LN_6>u unP7("A0D 图4}CQ_GvH 在“坐标系选项”一栏中,设置各项参数如下: Pc_J,Y\_"[ 源坐标系: 北京54坐标系;syV_&_Ds) 目的坐标系: 西安80坐标系;g+_)T_\_#u 转换方法: 七参数布尔莎模型;"_r.pU(uxt 长度单位: 米;ao_\+__%s 角度单位: 弧度;2NI3&;{_4 然后单击“添加项”按钮,则在窗口左边的“不同椭球间转换”列表中将该转换关系列出;&! m_;s_gi 在窗口下方的“参数设置”一栏中,将上一步得到的七个参数依次输入到相应的文本框中,如图4所示;k_\1qJr__ 单击“修改项”按钮,输入转换关系,并单击“确定”按钮;@[.%__A;E4 接下来就是文件投影的操作过程了。 _|o=eS&) 3、单击“投影转换”菜单下“MAPGIS投影转换/选转换线文件”命令,系统弹出“选择文件”对话框,如图5所示: 3___uhwoE P\N$__TYeH 图5_g]z_,_*d_ 选中待转换的文件“演示数据_北京54.WL”,单击“确定”按钮;d%5QE__VV 4、设置文件的Tic点,在“投影变换”模块下提供了两种方法: 手工设置和文件间拷贝,这里不作详细的说明;'mG[#M/_Y 5、单击“投影转换”菜单下“编辑当前投影参数”命令,系统弹出“输入投影参数”对话框,如图6所示,根据数据的实际情况来设置其地图参数,如下: t&_H): _P 坐标系类型: 大地坐标系_2_o? ! m2W 椭球参数: 北京54J_{\S+O2,* 投影类型: 高斯-克吕格投影_1k3wBc5< 比例尺分母: 1_x__G\&Q_E 坐标单位: 米q2A__x__-# 投影中心点经度(DMS): 1230000a[[u>oH_yd 然后单击“确定“按钮;4(|x@: wxm 8'@__pX< 图6KYu3dC'/,& 6、单击“投影转换”菜单下“设置转换后参数”命令,系统弹出“输入投影参数”对话框,如图7所示,转换后的参数设置为: &5_)K_g%r 坐标系类型: 大地坐标系__t'K_j\ 椭球参数: 西安80(注意椭球参数的变换)C_J%bBL'_. 投影类型: 高斯-克吕格投影4 %aO_Dr8 比例尺分母: 1/E@_LnK__e 坐标单位: 米dNB56E)5`J 投影中心点经度(DMS): 1230000(注意前后中央经线保持一致)6~ev5SD;f &JVe-. 图7Z6E-Fu__O 7、单击“投影转换”菜单下“进行投影变换”命令,系统弹出“输入转换后位移值”对话框,单击“开始转换”按钮,系统开始按照设定的参数转换线文件,如图8所示: G_sm5L wwF]+w%lOw 图8fpwge_/w_ 以同样的操作步骤和参数设置,将“演示数据_北京54.WL”、“演示数据_北京54.WP”文件进行投影转换;=mxmJF_A 8、单击鼠标右键,选择“复位”命令,系统弹出“选择文件名”对话框,可以看到系统生成了三个新的文件: “NEWLIN.WL”、“NEWPNT.WT”、“NEWPNT.WP”,依次选中这三个文件,单击“确定”按钮,如图7所示: t_\kI_(G 4_(B{-cK_ 这时新生成的三个文件就是西安80坐标系下的文件;hub1rY|No 补充: 通常情况下,转换过来的数据会有一定的误差存在,所以有时为了保证数据的精度,在转换的过程中通过设置横坐标和纵坐标的偏移量来修正转换后的坐标值; 跨带投影.S4Xw2_MS 我们知道高斯-克吕格投影采用了分带投影的思想,每一个投影带的坐标都是对本带坐标原点的相对值,所以各带的坐标完全相同,使用时只需变一个带号即可,这样就存在着如果不考虑带号的情况下,会有重叠的情况出现,如果要想将重叠的图框回到其实际所在的位置上,这时就会用到跨带投影。 cL}}^___ 跨带投影是MAPGIS投影变换中一个很重要的方面。 6_Ty;_m>j 下面来讲解跨带投影的操作方法和步骤,共分为两部分: mu_sZC_g$ 一、演示数据的生成和说明: qT} 在“投影变换”模块下分别生成3幅1: 50万的标准框,并在“输入编辑”模块中将其改成不同的颜色(“FRAM_50_左.W~“表示“FRAM_50_左.WL”和“FRAM_50_左.WT”): "._Z_|zt6C 名称起始经度(DDMMSS)起始纬度(DDMMSS)中 如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。 copyright@ 2008-2022 冰点文档网站版权所有 经营许可证编号:鄂ICP备2022015515号-1