【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

智能推荐

攻防世界_难度8_happy_puzzle_攻防世界困难模式攻略图文-程序员宅基地

文章浏览阅读645次。这个肯定是末尾的IDAT了,因为IDAT必须要满了才会开始一下个IDAT,这个明显就是末尾的IDAT了。,对应下面的create_head()代码。,对应下面的create_tail()代码。不要考虑爆破,我已经试了一下,太多情况了。题目来源:UNCTF。_攻防世界困难模式攻略图文

达梦数据库的导出(备份)、导入_达梦数据库导入导出-程序员宅基地

文章浏览阅读2.9k次,点赞3次,收藏10次。偶尔会用到,记录、分享。1. 数据库导出1.1 切换到dmdba用户su - dmdba1.2 进入达梦数据库安装路径的bin目录,执行导库操作  导出语句:./dexp cwy_init/[email protected]:5236 file=cwy_init.dmp log=cwy_init_exp.log 注释:   cwy_init/init_123..._达梦数据库导入导出

js引入kindeditor富文本编辑器的使用_kindeditor.js-程序员宅基地

文章浏览阅读1.9k次。1. 在官网上下载KindEditor文件,可以删掉不需要要到的jsp,asp,asp.net和php文件夹。接着把文件夹放到项目文件目录下。2. 修改html文件,在页面引入js文件:<script type="text/javascript" src="./kindeditor/kindeditor-all.js"></script><script type="text/javascript" src="./kindeditor/lang/zh-CN.js"_kindeditor.js

STM32学习过程记录11——基于STM32G431CBU6硬件SPI+DMA的高效WS2812B控制方法-程序员宅基地

文章浏览阅读2.3k次,点赞6次,收藏14次。SPI的详情简介不必赘述。假设我们通过SPI发送0xAA,我们的数据线就会变为10101010,通过修改不同的内容,即可修改SPI中0和1的持续时间。比如0xF0即为前半周期为高电平,后半周期为低电平的状态。在SPI的通信模式中,CPHA配置会影响该实验,下图展示了不同采样位置的SPI时序图[1]。CPOL = 0,CPHA = 1:CLK空闲状态 = 低电平,数据在下降沿采样,并在上升沿移出CPOL = 0,CPHA = 0:CLK空闲状态 = 低电平,数据在上升沿采样,并在下降沿移出。_stm32g431cbu6

计算机网络-数据链路层_接收方收到链路层数据后,使用crc检验后,余数为0,说明链路层的传输时可靠传输-程序员宅基地

文章浏览阅读1.2k次,点赞2次,收藏8次。数据链路层习题自测问题1.数据链路(即逻辑链路)与链路(即物理链路)有何区别?“电路接通了”与”数据链路接通了”的区别何在?2.数据链路层中的链路控制包括哪些功能?试讨论数据链路层做成可靠的链路层有哪些优点和缺点。3.网络适配器的作用是什么?网络适配器工作在哪一层?4.数据链路层的三个基本问题(帧定界、透明传输和差错检测)为什么都必须加以解决?5.如果在数据链路层不进行帧定界,会发生什么问题?6.PPP协议的主要特点是什么?为什么PPP不使用帧的编号?PPP适用于什么情况?为什么PPP协议不_接收方收到链路层数据后,使用crc检验后,余数为0,说明链路层的传输时可靠传输

软件测试工程师移民加拿大_无证移民,未受过软件工程师的教育(第1部分)-程序员宅基地

文章浏览阅读587次。软件测试工程师移民加拿大 无证移民,未受过软件工程师的教育(第1部分) (Undocumented Immigrant With No Education to Software Engineer(Part 1))Before I start, I want you to please bear with me on the way I write, I have very little gen...

随便推点

Thinkpad X250 secure boot failed 启动失败问题解决_安装完系统提示secureboot failure-程序员宅基地

文章浏览阅读304次。Thinkpad X250笔记本电脑,装的是FreeBSD,进入BIOS修改虚拟化配置(其后可能是误设置了安全开机),保存退出后系统无法启动,显示:secure boot failed ,把自己惊出一身冷汗,因为这台笔记本刚好还没开始做备份.....根据错误提示,到bios里面去找相关配置,在Security里面找到了Secure Boot选项,发现果然被设置为Enabled,将其修改为Disabled ,再开机,终于正常启动了。_安装完系统提示secureboot failure

C++如何做字符串分割(5种方法)_c++ 字符串分割-程序员宅基地

文章浏览阅读10w+次,点赞93次,收藏352次。1、用strtok函数进行字符串分割原型: char *strtok(char *str, const char *delim);功能:分解字符串为一组字符串。参数说明:str为要分解的字符串,delim为分隔符字符串。返回值:从str开头开始的一个个被分割的串。当没有被分割的串时则返回NULL。其它:strtok函数线程不安全,可以使用strtok_r替代。示例://借助strtok实现split#include <string.h>#include <stdio.h&_c++ 字符串分割

2013第四届蓝桥杯 C/C++本科A组 真题答案解析_2013年第四届c a组蓝桥杯省赛真题解答-程序员宅基地

文章浏览阅读2.3k次。1 .高斯日记 大数学家高斯有个好习惯:无论如何都要记日记。他的日记有个与众不同的地方,他从不注明年月日,而是用一个整数代替,比如:4210后来人们知道,那个整数就是日期,它表示那一天是高斯出生后的第几天。这或许也是个好习惯,它时时刻刻提醒着主人:日子又过去一天,还有多少时光可以用于浪费呢?高斯出生于:1777年4月30日。在高斯发现的一个重要定理的日记_2013年第四届c a组蓝桥杯省赛真题解答

基于供需算法优化的核极限学习机(KELM)分类算法-程序员宅基地

文章浏览阅读851次,点赞17次,收藏22次。摘要:本文利用供需算法对核极限学习机(KELM)进行优化,并用于分类。

metasploitable2渗透测试_metasploitable2怎么进入-程序员宅基地

文章浏览阅读1.1k次。一、系统弱密码登录1、在kali上执行命令行telnet 192.168.26.1292、Login和password都输入msfadmin3、登录成功,进入系统4、测试如下:二、MySQL弱密码登录:1、在kali上执行mysql –h 192.168.26.129 –u root2、登录成功,进入MySQL系统3、测试效果:三、PostgreSQL弱密码登录1、在Kali上执行psql -h 192.168.26.129 –U post..._metasploitable2怎么进入

Python学习之路:从入门到精通的指南_python人工智能开发从入门到精通pdf-程序员宅基地

文章浏览阅读257次。本文将为初学者提供Python学习的详细指南,从Python的历史、基础语法和数据类型到面向对象编程、模块和库的使用。通过本文,您将能够掌握Python编程的核心概念,为今后的编程学习和实践打下坚实基础。_python人工智能开发从入门到精通pdf

推荐文章

热门文章

相关标签