【20090603-02】地图投影的选择(转载)http://blog.csdn.net/chenyq2008/archive/2007/12/04/1915918.aspx...-程序员宅基地

技术标签: ViewUI  javascript  

地图投影的选择

选择投影的目的在于使所选投影的性质、特点适合于地图的用途,同时考虑地图在图廓范围内变形较小而且变形分布均匀。海域使用的地图多采用保角投影,因其能保持方位角度的正确。

我国的基本比例尺地形图(1:5千,1:1万,1:2.5万,1:5万,1:10万,1:25万,1:50万,1:100万)中,大于等于50万的均采用高斯-克吕格投影(Gauss-Kruger),这是一个等角横切椭圆柱投影,又叫横轴墨卡托投影(TransverseMercator);小于50万的地形图采用等角正轴割园锥投影,又叫兰勃特投影(LambertConformalConic);海上小于50万的地形图多用等角正轴圆柱投影,又叫墨卡托投影(Mercator)。一般应该采用与我国基本比例尺地形图系列一致的地图投影系统。
大地基准面的选择
地图坐标系由大地基准面和地图投影确定,大地基准面是利用特定椭球体对特定地区地球表面的逼近,因此每个国家或地区均有各自的大地基准面,我们通常称谓的北京54坐标系、西安80坐标系实际上指的是我国的两个大地基准面。我国参照前苏联从1953年起采用克拉索夫斯基(Krassovsky)椭球体建立了我国的北京54坐标系,1978年采用国际大地测量协会推荐的IAG75地球椭球体建立了我国新的大地坐标系--西安80坐标系,目前GPS定位所得出的结果都属于WGS84坐标系统,WGS84基准面采用WGS84椭球体,它是一地心坐标系,即以地心作为椭球体中心的坐标系。因此相对同一地理位置,不同的大地基准面,它们的经纬度坐标是有差异的。

程序中采用的3个椭球体参数如下(源自“全球定位系统测量规范GB/T18314-2001”):
椭球体 长半轴 短半轴
Krassovsky 6378245 6356863.0188
IAG75 6378140 6356755.2882
WGS84 6378137 6356752.3142

椭球体与大地基准面之间的关系是一对多的关系,也就是基准面是在椭球体基础上建立的,但椭球体不能代表基准面,同样的椭球体能定义不同的基准面,如前苏联的Pulkovo1942、非洲索马里的Afgooye基准面都采用了Krassovsky椭球体,但它们的大地基准面显然是不同的。在目前的GIS商用软件中,大地基准面都通过当地基准面向WGS84的转换7参数来定义,即三个平移参数ΔX、ΔY、ΔZ表示两坐标原点的平移值;三个旋转参数εx、εy、εz表示当地坐标系旋转至与地心坐标系平行时,分别绕Xt、Yt、Zt的旋转角;最后是比例校正因子,用于调整椭球大小。北京54、西安80相对WGS84的转换参数至今没有公开,实际工作中可利用工作区内已知的北京54或西安80坐标控制点进行与WGS84坐标值的转换,在只有一个已知控制点的情况下(往往如此),用已知点的北京54与WGS84坐标之差作为平移参数,当工作区范围不大时,如青岛市,精度也足够了。
以(32°,121°)的30°标准纬度墨卡托投影结果为例,北京54及WGS84基准面,两者投影结果在东西方向差距约196米,南北方向差距约57米(见下表),对于几十或几百万的地图来说,这一误差无足轻重,但在工程地图中还是应该加以考虑的。

输入坐标(度) 北京54墨卡托投影(米)
WGS84墨卡托投影(米)
纬度值(X) 32 3242287 3242230
经度值(Y) 121 11675036 11674840

墨卡托投影

(1)墨卡托投影性质

    墨卡托(Mercator)投影,又名"等角正轴圆柱投影”,荷兰地图学家墨卡托(GerhardusMercator1512-1594)在1569年拟定,我国的基本比例尺地形图(1:5千,1:1万,1:2.5万,1:5万,1:10万,1:25万,1:50万,1:100万)中,大于等于50万的均采用高斯-克吕格投影(Gauss-Kruger),又叫横轴墨卡托投影(Transverse Mercator);海上小于50万的地形图多用正轴等角园柱投影,又叫 墨卡托投影(Mercator)。假设地球被围在一中空的圆柱里,其赤道与圆柱相接触,然后再假想地球中心有一盏灯,把球面上的图形投影到圆柱体上,再把圆柱体展开,这就是一幅标准纬线为零的“墨卡托投影”绘制出的世界地图。
     墨卡托投影没有角度变形,由每一点向各方向的长度比相等,它的经纬线都是平行直线,且相交成直角,经线间隔相等,纬线间隔从标准纬线向两极逐渐增大。墨卡托投影的地图上长度和面积变形明显,标准纬线无变形,从标准纬线向两极变形增大,但因为它具有各个方向均等扩大的特性,保持了方向和相互位置关系的正确。在地图上保持方向和角度的正确是墨卡托投影的优点,墨卡托投影地图常用作航海图和航空图,如果循着墨卡托投影图上两点间的直线航行,方向不变可以一直到达目的地,因此它对船舰在航行中定位、确定航向都具有有利条件,给航海者带来很大方便。

(2)墨卡托投影坐标

取零子午线或自定义原点经线与赤道交点的投影为原点,零子午线或自定义原点经线的投影为纵坐标x轴,赤道的投影为横坐标y轴,构成墨卡托平面直角坐标系,此投影标准纬线无变形。

单点转换
单点转换步骤如下:

(1)选择是墨卡托正转换还是反转换,缺省为经纬度转换到墨卡托投影坐标,投影坐标单位为米。
(2)选择大地基准面,缺省北京54,如果是GPS定位数据别忘了切换为WGS84。
(3)输入标准纬度,单位度。
(4)输入原点经度,单位度,缺省为零。
(5)如正向投影,选择经纬度输入数据格式,有三个选项,缺省为十进制度格式。具体输入方式如下例:
格式 原始纬度值 原始经度值 输入纬度值 输入经度值
十进制度 35.445901° 122.997344° 35.445901 122.997344
度分 35°26.7541′ 122°59.8406′ 3526.7541 12259.8406
度分秒 35°26′45.245″ 122°59′50.438″
352645.245
1225950.438

(6)正投影按选定格式在“输入”栏输入经纬度值,反投影输入以米为单位的X、Y坐标值。
(7)单击“单点转换”按钮。
(8)在“输出”栏查看计算结果。

批量转换
批量转换步骤如下:

(1)准备好需要转换的输入数据文件,要求是文本文件,分两列,第一列纬度值或纵向坐标值,第二列经度值或横向坐标值,两列之间用空格分开。正向投影时,纬度值及经度值格式可以有三种选择(见表),缺省当作十进制度处理;反向投影时,纵向及横向坐标值必须以米为单位。

下例为度分秒格式(WGS84)的墨卡托正投影输入数据文件testdata.txt
350000.000 1220000.000
351600.519 1225959.506
345800.101 1225959.8
343600.336 1230000.26
341400.018 1225959.897
335159.17 1225959.46
333000.08 1230000.28
332300.674 1232000.103

(2)选择是墨卡托正转换还是反转换,缺省为经纬度转换到墨卡托投影坐标,投影坐标单位为米。
(3)选择大地基准面,缺省北京54,如果是GPS定位数据别忘了切换为WGS84。
(4)输入标准纬度,单位度,如“30”,表示标准纬度30°。
(5)输入原点经度,单位度,缺省为零。
(6)如正向投影,选择输入数据文件中的经纬度输入数据格式,有三个选项,缺省为十进制度格式。
(7)单击“批量转换”按钮。弹出打开文件对话框,输入你的数据文件名。
(8)输入转换结果文件名,单击“保存”后,程序开始进行计算。
(9)打开输出文件查看计算结果,结果分五列,第一序号,第二列输入纬度值或纵向坐标值,第三列输入经度值或横向坐标值,第四列转换后纬度值或纵向坐标值,第五列转换后经度值或横向坐标值。
下例为度分秒格式(WGS84)的标准纬度30°,原点经度0°的正投影转换结果数据文件result.txt
1 350000 1220000 3587805.5 11771326.2
2 351600.519 1225959.506 3619142.9 11867799.2
3 345800.101 1225959.8 3583901 11867807.1
4 343600.336 1230000.26 3541028.3 11867819.4
5 341400.018 1225959.897 3498328.2 11867809.7
6 335159.17 1225959.46 3455798.3 11867798
7 333000.08 1230000.28 3413508.8 11867820
8 332300.674 1232000.103 3400100.7 11899977.3

最新回复

pal at 2007-7-28 00:36:52
1、首先理解地理坐标系(Geographic coordinate system),Geographic coordinate system直译为地理坐标系统,是以经纬度为地图的存储单位的。很明显,Geographic coordinate syst em是球面坐标系统。我们要将地球上的数字化信息存放到球面坐标系统上,如何进行操作呢?地球是一个不规则的椭球,如何将数据信息以科学的方法存放到椭球上?这必然要求我们找到这样的一个椭球体。这样的椭球体具有特点:可以量化计算的。具有长半轴,短半轴,偏心率。以下几行便是Krasovsky_1940椭球及其相应参数。
Spheroid: Krasovsky_1940
Semimajor Axis: 6378245.000000000000000000
Semiminor Axis: 6356863.018773047300000000
Inverse Flattening(扁率): 298.300000000000010000
然而有了这个椭球体以后还不够,还需要一个大地基准面将这个椭球定位。在坐标系统描述中,可以看到有这么一行:
Datum: D_Beijing_1954
表示,大地基准面是D_Beijing_1954。

有了Spheroid和Datum两个基本条件,地理坐标系统便可以使用。
完整参数:
Alias:
Abbreviation:
Remarks:
Angular Unit: Degree (0.017453292519943299)
Prime Meridian(起始经度): Greenwich (0.000000000000000000)
Datum(大地基准面): D_Beijing_1954
Spheroid(参考椭球体): Krasovsky_1940
Semimajor Axis: 6378245.000000000000000000
Semiminor Axis: 6356863.018773047300000000
Inverse Flattening: 298.300000000000010000
2、接下来便是Projection coordinate system(投影坐标系统),首先看看投影坐标系统中的一些参数。
Projection: Gauss_Kruger
Parameters:
False_Easting: 500000.000000
False_Northing: 0.000000
Central_Meridian: 117.000000
Scale_Factor: 1.000000
Latitude_Of_Origin: 0.000000
Linear Unit: Meter (1.000000)
Geographic Coordinate System:
Name: GCS_Beijing_1954
Alias:
Abbreviation:
Remarks:
Angular Unit: Degree (0.017453292519943299)
Prime Meridian: Greenwich (0.000000000000000000)
Datum: D_Beijing_1954
Spheroid: Krasovsky_1940
Semimajor Axis: 6378245.000000000000000000
Semiminor Axis: 6356863.018773047300000000
Inverse Flattening: 298.300000000000010000
从参数中可以看出,每一个投影坐标系统都必定会有Geographic Coordinate System。
投影坐标系统,实质上便是平面坐标系统,其地图单位通常为米。那么为什么投影坐标系统中要存在坐标系统的参数呢? 这时候,又要说明一下投影的意义: 将球面坐标转化为平面坐标的过程便称为投影。 好了,投影的条件就出来了:
a、球面坐标
b、转化过程(也就是算法)
也就是说,要得到投影坐标就必须得有一个“拿来”投影的球面坐标,然后才能使用算法去投影!即每一个投影坐标系统都必须要求有Geographic Coordinate System参数。
3、我们现在看到的很多教材上的对坐标系统的称呼很多,都可以归结为上述两种投影。其中包括我们常见的“非地球投影坐标系统”。
大地坐标(Geodetic Coordinate):大地测量中以参考椭球面为基准面的坐标。地面点P的位置用大地经度L、大地纬度B和大地高H表示。当点在参考椭球面上时,仅用大地经度和大地纬度表示。大地经度是通过该点的大地子午面与起始大地子午面之间的夹角,大地纬度是通过该点的法线与赤道面的夹角,大地高是地面点沿法线到参考椭球面的距离。
方里网:是由平行于投影坐标轴的两组平行线所构成的方格网。因为是每隔整公里绘出坐标纵线和坐标横线,所以称之为方里网,由于方 里线同时 又是平行于直角坐标轴的坐标网线,故又称直角坐标网。
在1:1万——1:20万比例尺的地形图上,经纬线只以图廓线的形式直接表现出来,并在图角处注出相应度数。为了在用图时加密成 网,在内外图廓间还绘有加密经纬网的加密分划短线(图式中称“分度带”),必要时对应短线相连就可以构成加密的经纬线网。1:2 5万地形图上,除内图廓上绘有经纬网的加密分划外,图内还有加密用的十字线。
我国的1:50万——1:100万地形图,在图面上直接绘出经纬线网,内图廓上也有供加密经纬线网的加密分划短线。
直角坐标网的坐标系以中央经线投影后的直线为X轴,以赤道投影后的直线为Y轴,它们的交点为坐标原点。这样,坐标系中就出现了四 个象限。纵坐标从赤道算起向北为正、向南为负;横坐标从中央经线算起,向东为正、向西为负。
虽然我们可以认为方里网是直角坐标,大地坐标就是球面坐标。但是我们在一副地形图上经常见到方里网和经纬度网,我们很习惯的称经 纬度网为大地坐标,这个时候的大地坐标不是球面坐标,她与方里网的投影是一样的(一般为高斯),也是平面坐标。
pal at 2007-7-28 00:39:42
2楼介绍了地理坐标系与投影坐标系的区别简言之, 投影坐标系=地理坐标系+投影过程

ArcGIS自带了多种坐标系统,在${ArcGISHome}\Coordinate Systems\目录下可以看到三个文件夹,分别是Geographic Coordinate Systems、Projected Coordinate Systems、Vertical Coordinate Systems,中文翻译为地理坐标系投影坐标系垂直坐标系

1 Geographic Coordinate Systems地理坐标系

      在Geographic Coordinate Systems目录中,我们可以看到已定义的许多坐标系信息,如典型的Geographic Coordinate Systems\World目录下的WGS 1984.prj,里面所定义的坐标参数描述了地理坐标系的名称、大地基准面、椭球体、起始坐标参考点、单位等:

GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]]

2 Projected Coordinate Systems投影坐标系

      在Projected Coordinate Systems目录中同样存在许多已定义的投影坐标系,我国大部分地图所采用的北京54和西安80坐标系的投影文件就在其中,它们均使用高斯-克吕格投影,前者使用克拉索夫斯基椭球体,后者使用国际大地测量协会推荐的IAG 75地球椭球体。如Beijing 1954 3 Degree GK CM 75E.prj定义的坐标参数:
      PROJCS["Beijing_1954_3_Degree_GK_CM_75E",GEOGCS["GCS_Beijing_1954",DATUM["D_Beijing_1954",SPHEROID["Krasovsky_1940",6378245.0,298.3]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Gauss_Kruger"],PARAMETER["False_Easting",500000.0],PARAMETER
["False_Northing",0.0],PARAMETER["Central_Meridian",75.0],PARAMETER["Scale_Factor",1.0],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]]

      可以看出,参数里除了包含地理坐标系的定义外,还有投影方式的信息。北京54和西安80是我们使用最多的坐标系,在ArcGIS文件中,对于这两种坐标系统的命名有一些不同,看上去很容易让人产生迷惑。

在此先介绍下高斯-克吕格投影的基本知识:我国大中比例尺地图均采用高斯-克吕格投影Gauss Kruger其通常是按6度和3度分带投影,1:2.5万-1:50万比例尺地形图采用经差6度分带,1:1万比例尺的地形图采用经差3度分带。具体分带法是:6度分带从本初子午线开始,按经差6度为一个投影带自西向东划分,全球共分60个投影带,带号分别为1-60;3度投影带是从东经1度30秒经线开始,按经差3度为一个投影带自西向东划分,全球共分120个投影带。为了便于地形图的测量作业,在高斯-克吕格投影带内布置了平面直角坐标系统,具体方法是,规定中央经线为X赤道为Y,中央经线与赤道交点为坐标原点,x值在北半球为正,南半球为负,y值在中央经线以东为正,中央经线以西为负。由于我国疆域均在北半球,x值均为正值,为了避免y值出现负值,规定各投影带的坐标纵轴均西移500km,中央经线上原横坐标值由0变为500km。为了方便带间点位的区分,可以在每个点位横坐标y值的百千米位数前加上所在带号,如20带内A点的坐标可以表示为YA=20 745 921.8m。
在Coordinate Systems\Projected Coordinate Systems\Gauss Kruger\Beijing 1954目录中,我们可以看到四种不同的命名方式:

      Beijing 1954 3 Degree GK CM 75E.prj
      Beijing 1954 3 Degree GK Zone 25.prj
      Beijing 1954 GK Zone 13.prj
      Beijing 1954 GK Zone 13N.prj

      对它们的说明分别如下:

      三度分带法的北京54坐标系,中央经线在东75度的分带坐标,横坐标前不加带号
      三度分带法的北京54坐标系,中央经线在东75度的分带坐标,横坐标前加带号
      六度分带法的北京54坐标系,分带号为13,横坐标前加带号
      六度分带法的北京54坐标系,分带号为13,横坐标前不加带号

      在Coordinate Systems\Projected Coordinate Systems\Gauss Kruger\Xian 1980目录中,文件命名方式又有所变化:

      Xian 1980 3 Degree GK CM 75E.prj
      Xian 1980 3 Degree GK Zone 25.prj
      Xian 1980 GK CM 75E.prj
      Xian 1980 GK Zone 13.prj

      西安80坐标文件的命名方式、含义和北京54前两个坐标相同,但没有出现“带号+N”这种形式。


3 Vertical Coordinate Systems 垂直坐标系

    Vertical Coordinate Systems定义了测量海拔或深度值的原点,具体的定义,英文描述的更为准确:

    A vertical coordinate system defines the origin for height or depth values. Like a horizontal coordinate system, most of the information in a vertical coordinate system is not needed unless you want to display or combine a dataset with other data that uses a different vertical coordinate system.
    Perhaps the most important part of a vertical coordinate system is its unit of measure. The unit of measure is always linear (e.g., international feet or meters). Another important part is whether the z values represent heights (elevations) or depths. For each type, the z-axis direction is positive "up" or "down", respectively.
    One z value is shown for the height-based mean sea level system. Any point that falls below the mean sea level line but is referenced to it will have a negative z value. The mean low water system has two z values associated with it. Because the mean low water system is depth-based, the z values are positive. Any point that falls above the mean low water line but is referenced to it will have a negative z value.

    注意:大家经常希望能够通过坐标转换,将北京54或西安80中的地理坐标系转换到WGS84,实际上这样做是不准确的,北京54或西安80的投影坐标可以通过计算转换到其对应的地理坐标系,但由于我国北京54和西安80中的地理坐标系到WGS84的转换参数没有公开,因此无法完成其到WGS84坐标的精准计算。其他公开了转换参数的坐标系都可以在ArcToolbox中完成转换。
pal at 2007-7-28 00:53:22
高斯-克吕格尔平面直角坐标系 Gauss-Krüger plane rectangular coordinates system
  根据高斯-克吕格尔投影所建立的平面坐标系,或简称高斯平面坐标系。它是大地测量、城市测量、普通测量、各种工程测量和地图制图中广泛采用的一种平面坐标系。
  高斯-克吕格尔投影是德国的 C.F.高斯于1822年提出的,后经德国的克吕格尔(J.H.L.Krüger)于1912年加以扩充而完善。
  用大地经度和纬度表示的大地坐标是一种椭球面上的坐标,不能直接应用于测图。因此,需要将它们按一定的数学规律转换为平面直角坐标。大地坐标(B,L)转换为平面直角坐标(X,Y)的一般数学表示法为:X=F1(B,L), Y=F2(B,L), 式中F1、F2为投影函数。高斯-克吕格尔投影的投影函数是根据以下两个条件确定的:第一,投影是正形的,即椭球面上无穷小的图形和它在平面上的表象相似,故又称保角投影或保形投影;投影面上任一点的长度比(该点在椭球面上的微分距离与其在平面上相应的微分距离之比)同方位无关。第二,椭球面上某一子午线在投影平面上的表象是一直线,而且长度保持不变,即长度比等于1。该子午线称为中央子午线,或称轴子午线。这两个条件体现了高斯-克吕格尔投影的特性。

大地坐标系是大地测量的基本坐标系。常用于大地问题的细算,研究地球形状和大小,编制地图,火箭和卫星发射及军事方面的定位及运算,若将其直接用于工程建设规划、设计、施工等很不方便。所以要将球面上的大地坐标按一定数学法则归算到平面上,即采用地图投影的理论绘制地形图,才能用于规划建设。
    椭球体面是一个不可直接展开的曲面,故将椭球体面上的元素按一定条件投影到平面上,总会产生变形。测量上常以投影变形不影响工程要求为条件选择投影方法。地图投影有等角投影、等面积投影和任意投影三种。
    其中等角投影又称为正形投影,它保证在椭球体面上的微分图形投影到平面后将保持相似。这是地形图的基本要求。正形投影有两个基本条件:
        ①保角条件,即投影后角度大小不变。
        ②长度变形固定性,即长度投影后会变形,但是在一点上各个方向的微分线段变形比m是个常数k:     
                                       式中:ds—投影后的长度,dS—球面上的长度。

1.高斯投影的概念

        高斯是德国杰出的数学家、测量学家。高斯-克吕格尔投影是德国的 C.F.高斯于1822年提出的,后经德国的克吕格尔(J.H.L.Krüger)于1912年加以扩充而完善。他提出的横椭圆柱投影是一种正形投影。它是将一个横椭圆柱套在地球椭球体上,如下图所示:


1.gif


        椭球体中心O在椭圆柱中心轴上,椭球体南北极与椭圆柱相切,并使某一子午线与椭圆柱相切。此子午线称中央子午线。然后将椭球体面上的点、线按正形投影条件投影到椭圆柱上,再沿椭圆柱N、S点母线割开,并展成平面,即成为高斯投影平面。在此平面上:
        ①中央子午线是直线,其长度不变形,离开中央子午线的其他子午线是弧形,凹向中央子午线。离开中央子午线越远,变形越大。
        ②投影后赤道是一条直线,赤道与中央子午线保持正交。
        ③离开赤道的纬线是弧线,凸向赤道。
        
        高斯投影可以将椭球面变成平面,但是离开中央子午线越远变形越大,这种变形将会影响测图和施工精度。为了对长度变形加以控制,测量中采用了限制投影宽度的方法,即将投影区域限制在靠近中央子午线的两侧狭长地带。这种方法称为分带投影。投影带宽度是以相邻两个子午线的经差来划分。有6°带、3°带等不同投影方法。
        6°带投影是从英国格林尼治子午线开始,自西向东,每隔6°投影一次。这样将椭球分成60个带,编号为1~60带,如下图所示:


2.jpg


        各带中央子午线经度(L)可用下式计算:
                                          式中n为6°带的带号。
                 
        已知某点大地经度L,可按下式计算该点所属的带号:
                                            有余数时,为n的整数商+1。

        3°带是在6°带基础上划分的,其中央子午线在奇数带时与6°带中央子午线重合,每隔3°为一带,共120带,各带中央子午线经度(L)为:
                                           式中n′为3°带的带号。

        我国幅员辽阔,含有11个6°带,即从13~23带(中央子午线从75°~135°),21个3°带,从25~45带。北京位于6°带的第20带,中央子午线经度为117°。

2.高斯平面直角坐标系Gauss-Krüger plane rectangular coordinates system

根据高斯-克吕格尔投影所建立的平面坐标系,或简称高斯平面坐标系。它是大地测量、城市测量、普通测量、各种工程测量和地图制图中广泛采用的一种平面坐标系。
  用大地经度和纬度表示的大地坐标是一种椭球面上的坐标,不能直接应用于测图。因此,需要将它们按一定的数学规律转换为平面直角坐标。大地坐标(B,L)转换为平面直角坐标(X,Y)的一般数学表示法为:X=F1(B,L), Y=F2(B,L), 式中F1、F2为投影函数。高斯-克吕格尔投影的投影函数是根据以下两个条件确定的:第一,投影是正形的,即椭球面上无穷小的图形和它在平面上的表象相似,故又称保角投影或保形投影;投影面上任一点的长度比(该点在椭球面上的微分距离与其在平面上相应的微分距离之比)同方位无关。第二,椭球面上某一子午线在投影平面上的表象是一直线,而且长度保持不变,即长度比等于1。该子午线称为中央子午线,或称轴子午线。这两个条件体现了高斯-克吕格尔投影的特性。
        根据高斯投影的特点,以赤道和中央子午线的交点为坐标原点。,中央子午线方向为x轴,北方向为正。赤道投影线为y轴,东方向为正。象限按顺时针Ⅰ、Ⅱ、Ⅲ、Ⅳ排列,如下图所示:


3.gif


        在同一投影带内y值有正有负。这对计算和使用很不方便。为了使y值都为正,将纵坐标轴西移500km,并在y坐标前面冠以带号,如在第20带,中央子午线以西P点:
                                


4.jpg


        在20带中高斯直角坐标为:
                                 


5.jpg


        高斯直角坐标系与数学中的笛卡尔坐标系不同,如下图所示:


6.gif


        高斯直角坐标系纵坐标为x轴,横坐标为y轴。坐标象限为顺时针划分四个象限。角度起算是从x轴的北方向开始,顺时针计算。这些定义都与数学中的定义不同。这样的做法是为了将数学上的三角和解析几何公式直接用到测量的计算上。

  中国于50年代正式决定在大地测量和国家地形图中采用高斯-克吕格尔平面直角坐标系。
  中国除了天文大地网平差采用椭球面上的大地坐标之外,高斯平面直角坐标系被广泛应用于其他各等大地控制网的平差和计算中。为此,一般先将椭球面上的方向、角度、长度等观测元素经方向改化和距离改化,归化为相应的平面观测值,然后在平面上进行平差和计算,这要比直接在地球椭球面上进行简单得多。
  大地坐标、大地线长度和大地方位角与高斯平面上相应的直角坐标,平面边长和坐标方位角之间的相互换算工作,一般是借助于专门的计算用表进行,或者直接在电子计算机上进行。
  通用横轴墨卡托投影  高斯-克吕格尔投影的一种变体,简称UTM投影。它同高斯-克吕格尔投影的差别仅在于中央子午线的长度比不是1,而是0.9996。UTM投影带中的两条标准线在中央子午线东、西各约 180公里处,这两条标准线上没有任何变形,离开这两条线愈远变形愈大。在这两条线之内长度缩小,两线之外长度放大。UTM投影应用比较广泛,目前世界上已有100多个国家和地区采用这种投影作为南纬80°至北纬84°的地区中测制地形图的数学基础。
  参考书目:方俊《地图投影学》第二册,科学出版社,北京,1958。B.G.Bomford, Geodesy,3rd ed.,Oxford Univ.Press,Oxford,1971.

转载于:https://www.cnblogs.com/willwayer/archive/2009/06/03/1495476.html

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/weixin_30709809/article/details/97393295

智能推荐

王斌老师的博客_王斌 github-程序员宅基地

文章浏览阅读480次。http://blog.sina.com.cn/s/blog_736d0b9101018cgc.html_王斌 github

ACM OJ Collection_htt//acm.wydtang.top/-程序员宅基地

文章浏览阅读737次。原文来自:http://blog.csdn.net/hncqp/article/details/4455263 ACM OJ Collection(排名不分先后):中国:浙江大学(ZJU):http://acm.zju.edu.cn/北京大学(PKU):htt_htt//acm.wydtang.top/

ios 自己服务器 苹果支付_修复苹果IOS支付-程序员宅基地

文章浏览阅读467次。更新记录1.0.0(2019-07-01)插件简介专门用来修复苹果IOS支付时出现"您已购买此App内购买项目。此项目将免费恢复"。问题描述首先在IOS平台里面创建“APP内购买项目”,选择的是“消耗型项目”,然后用uni-app官方的支付api进行支付,多支付几次,有时候就会出现提示“您已购买此App内购买项目。此项目将免费恢复”,特别是在沙盒测试里面支付很大几率出现,我明明选的是消耗型项目,应..._ios开发苹果支付恢复权益

spring MVC mock类单元测试(controller)_mvcmock-程序员宅基地

文章浏览阅读5.6k次。Spring从J2EE的Web端为每个关键接口提供了一个mock实现:MockHttpServletRequest几乎每个单元测试中都要使用这个类,它是J2EE Web应用程序最常用的接口HttpServletRequest的mock实现。MockHttpServletResponse此对象用于HttpServletRespons_mvcmock

【我的世界Minecraft-MC】常见及各种指令大杂烩【2022.8版】_summon生成掉落物-程序员宅基地

文章浏览阅读8.5k次,点赞7次,收藏22次。execute as @a at @s run clear @s minecraft:dark_oak_planks{display:{Name:“{“text”:“第三关[阴森古堡]”,“color”:“red”,“italic”:false}”,color:“16711680”},Enchantments:[{id:“protection”,lvl:1}],Unbreakable:1b} 1。Lore:[“{“text”:“免费”,“color”:“blue”,“italic”:false}”]..._summon生成掉落物

CentOS 7安装教程(图文详解)_centos 安装-程序员宅基地

文章浏览阅读10w+次,点赞487次,收藏2.1k次。CentOS 7安装教程: 准备: 软件:VMware Workstation 镜像文件:CentOS-7-x86_64-bin-DVD1.iso (附:教程较为详细,注释较多,故将操作的选项进行了加粗字体显示。) 1、文件--新建虚拟机--自定义 2、..._centos 安装

随便推点

Github项目分享——免费的画图工具drow,前端插件化面试_draw github画图-程序员宅基地

文章浏览阅读333次,点赞3次,收藏3次。项目介绍一款很好用的免费画图软件,支持ER图、时序图、流程图等等在项目的releases就可以下载最新版本同时支持在线编辑。_draw github画图

如何开始学习人工智能?入门的学习路径和资源是什么?_人工智能学习路径-程序员宅基地

文章浏览阅读930次。嗨,大家好!如果你对人工智能充满了好奇,并且想要入门这个领域,那么你来对地方了。本文将向你介绍如何从零基础开始学习人工智能,并逐步掌握核心概念和技能。无论你是大学生、职场新人还是对人工智能感兴趣的任何人,都可以按照以下学习路径逐步提升自己。_人工智能学习路径

Unity3D 导入资源_unity怎么导入压缩包-程序员宅基地

文章浏览阅读4.3k次,点赞2次,收藏8次。打开Unity3D的:window-asset store就会出来这样的界面:我们选择一个天空纹理,注意这里的标签只有一个,如果有多个就会显示所有标签的内容:找个比较小的免费的下载一下试试,比如这个:下载以后:点击import就会出现该窗口:然后再点击最底下的import:就导入到我们这里来了。从上面可以切换场景:..._unity怎么导入压缩包

jqgrid 服务器端验证,javascript – jqgrid服务器端错误消息/验证处理-程序员宅基地

文章浏览阅读254次。在你以前的问题的the answer的最后一部分,我试着给出你当前的问题的答案.也许我表示不够清楚.您不应该将错误信息放在标准成功响应中.您应该遵循用于服务器和客户端之间通信的HTTP协议的主要规则.根据HTTP协议实现网格中的加载数据,编辑行和与服务器的所有Ajax通信.每个HTTP响应都有响应第一行的状态代码.了解这个意义非常重要.典型的JSON数据成功请求如下HTTP/1.1 200 OK...._decode message error

白山头讲PV: 用calibre进行layout之间的比对-程序员宅基地

文章浏览阅读4k次,点赞8次,收藏29次。我们在流片之后,通常还是有机会对layout进行局部小的修改。例如metal change eco或者一些层次的局部修改。当我们修改之后,需要进行与之前gds的对比,以便确认没有因为某些..._calibre dbdiff

java exit方法_Java:如何测试调用System.exit()的方法?-程序员宅基地

文章浏览阅读694次。问题我有一些方法应该在某些输入上调用567779278。不幸的是,测试这些情况会导致JUnit终止!将方法调用放在新线程中似乎没有帮助,因为System.exit()终止了JVM,而不仅仅是当前线程。是否有任何常见的处理方式?例如,我可以将存根替换为System.exit()吗?[编辑]有问题的类实际上是一个命令行工具,我试图在JUnit中测试。也许JUnit根本不适合这份工作?建议使用互补回归测..._检查system.exit