MAPGIS投影变换系统的概念与应用(1)

<P>投影变换的基本概念</P>
<P>地图投影的基本问题乃是如何将地球表面（椭球面或圆球面）表示在地图平面上。由于地球椭球面或圆球面是不可展开的曲面，即不可能展开成水面。而地图又必须是一个平面，所以将地球表面展开成地图平面必然产生裂隙或褶皱。在编制地图中，要求（在地图上）得到连续的经纬网格，有裂隙或重叠的经纬网格不符合编图要求。如果采用简单的方法将有裂隙或重叠的经纬网“拉伸”或“压缩”使之连续，也不能符合编制地图的实际需要，因此，编制地图必须采用地图投影的方法，将地球表面的经纬网格用各种方法投影到地图平面上，才能满足多种编图要求的各种连续的经纬网格。</P>
<P>    地图投影即为地球椭球表面(或地球球体表面)与地图平面之间点与点(或线与线)相对应。如设地球表面上某一点的地理坐标为Q、λ，其地图平面上相应点的直角坐标为ｘ，y，则表示地球表面经线和纬线的两族平面曲线的方程为：</P>
<P> Q = F1(x,y)     λ= F2(x, y)</P>
<P>    对x和y分别得出的：     x = f1(Q, λ)   y = f2(Q, λ)</P>
<P>当函数f1和f2在表象区域内连续和单值时，上述方程表面曲面上坐标为Q，λ的一点在平面上只有一个点与之对应，我们把这个方程视为地图投影的一般方程。可以说，地图投影就是建立地球表面上点（Q，λ）和平面上的点（ｘ，ｙ）之间的函数关系式。地图投影学就是研究地球表面上的点或线表示在地图平面上的理论与方法。<BR>地图投影的基本要素</P>
<P>大地测量中用水准测量方法得到的地面上各点的高程是依据一个理想的水准面来确定的，这个水准面称为大地水准面。大地水准面所包围的球体称为大地球体。以一个大小和形状同它极为接近的旋转椭球面来代替，以随圆的短轴（地轴）为轴旋转而成的椭球面称为地球椭球面。</P>
<P>椭球体的元素与公式如下：</P>
<P>    扁率：          af=(a-b)/a</P>
<P>    第一偏心率：    e2=（a2-b2）/a2</P>
<P>    第二偏心率:     ep2=(a2-b2)/b2</P>
<P>  其中：长半径a(赤道半径)，短半径b(极轴半径)，扁率af，第一偏心率e，第二偏心率ep。</P>
<P>地球椭球体表面上的基本概念:</P>
<P>    赤道平面：垂直于地轴，并通过地心的平面。</P>
<P>    赤道：赤道平面与地球表面相交的大圆圈（交线）。</P>
<P>    纬圈或（纬线）：平行于赤道的各个圆圈。显然赤道是最大的一个平行圈，它的半径为ａ。</P>
<P>    线圈（经线）或子午圈：通过地球表面垂直于赤道面的平面叫经面或子午面，它和地球表面相交的线。所有的子午圈长度彼此都相等，均为一长半径等于  ａ，短半径等于ｂ的椭圆。</P>
<P>    法线：由地球椭球体上的任一点可以引一垂线垂直于该点的地平线（切线），这条垂线称为法线。</P>
<P>    地理纬度（简称纬度）：法线与赤道面相交所构成的角。通常用希腊字母Q表示，纬度以赤道为０°，向北、南两极各以９０°计算，向北叫北纬，向南叫南纬。</P>
<P>    地理经度（简称经度）：通过地球表面上某一点的纬线面与起始经线表面的夹角。通常以然腊字母λ表示。</P>
<P>  地面上任一点的位置，在测绘工作中通常是用经度（λ）和纬度（Q）来决定，写成Ｍ（φ，λ）。</P>
<P>    地理坐标系：经线和纬线是地球表面上两组正交（相交为９０°）的曲线，这两组正交的曲线构成的坐标，称为地理坐标系。</P>
<P>    经差：地表面某两点经度值之差。</P>
<P>    纬差：某两点纬度值之差。</P>
<P>  如若两点在同一经线上，其经差为零，如在同一纬线上，其纬差为零。</P>
<P> 关于地球椭球体的大小，由于采用不同的资料推算，椭球体的元素值是不同的。</P>
<P>世界各国常用的地球椭球体的数据：</P>
<P>椭球体名称 年代 长半轴a 短半轴b 扁率e<BR>白塞尔 1841 6377397 6356079 1/299.15<BR>克拉克 1880 6378249 6356515 1/293.5<BR>克拉索夫斯基 1940 6378245 6356863 1/298.3<BR>第十六界IUUG 1975 6378140 6356755 1/298.26<BR>第十七界IUUG 1979 6378137 6356752 1/298.257<BR>第十八界IUUG 1983 6378136 6356751 1/298.257<BR>    我国1952年以前采用海福特椭球（该椭球1924年被定为国际椭球），从1953年起，开始改用克拉索夫斯基椭球。1978年我国决定采用国际大地测量协会所推荐的“1975年基本大地数据”中给定的椭球参数，并以此建立了我国新的、独立的大地坐标系，形成了1980年的西安坐标系。<BR>    等面积球体半径：使球体的表面积等于地球椭球体表面积。<BR>    等体积球体半径：使球体的体积等于椭球体的体积。<BR>地图投影的分类</P>
<P>根据变形性质的投影分类</P>
<P> 等角投影</P>
<P> 等面积投影</P>
<P> 任意投影</P>
<P>根据正轴投影时经纬网形状的投影分类</P>
<P>     圆锥投影              圆柱投影 </P>
<P>     方位投影              伪圆锥投影</P>
<P>     伪圆柱投影            伪方位投影</P>
<P>     多圆锥投影</P>
<P>经纬线形状：</P>
<P>    对于一个投影，较完整的名称宜兼有两种分类，例如等角圆锥投影，等面积方位投影，等距离圆柱投影等</P>
<P>地图投影变换<BR>    地图投影变换是研究从一种地图投影点的坐标变换为另一种地图投影点的坐标的理论和方法。在大地测量和地形测量中，往往需要进行不同坐标系间的坐标变换，即坐标换带计算。随着计算机地图制图的发展，研究地图投影变换的理论和方法日益重要和迫切，因为在采用制图自动仪作业中，必须首先提供从一种地图投影点的坐标变换为另一种地图投影点的坐标的关系式即数学模式，才能进行这种作业，因为如果没有这两种不同投影点的坐标变换关系式，就无法编制出合乎变换要求的，适用于电子计算机进行变换所需要的程序设计，所以，地图投影变换已成为计算机地图制图的一个组成部分。</P>
<P>    为了适应计算机地图制图，信息系统建设，遥测图象显示和定位，以及满足空间遥感技术和拓朴变换的需要，地图投影变换已逐步发展成为研究空间数据处理，以及空间点位和平面点变换的理论和方法及其应用的数学制图学的一个分支学科。地图投影变换可广义地理解为研究空间数据处理，变换及应用的理论和方法，它可表达为：</P>
<P>(x'，y'，z')－－－(φ，λ)－－－(x，y)－－－(X，Y)</P>
<P>地图投影变换可狭义地理解为建立两平面场之间点的一一对应的函数关系。</P>
<P>高斯-克吕格投影</P>
<P>高斯-克吕格(GAUSS-KRUGER)是等角横切椭圆柱投影，由德国数学家高斯提出，后经克吕格扩充并推倒出计算公式，故称为高斯-克吕格投影，简称高斯投影。该投影以中央经线和赤道投影后为坐标轴，中央经线和赤道交点为坐标原点，纵坐标由坐标原点向北为正，向南为负，规定为X轴，横坐标从中央经线起算，向东为正，向西为负，规定为Y轴。所以，高斯-克吕格坐标系的X、Y轴正好对应MAPGIS坐标系的Y和X。</P>
<P>为了控制变形，本投影采用分带的办法。我国1:2.5-1:50万地形图均采用6度分带；1:1万及更大比例尺地形图采用3度分带，以保证必要的精度。6度分带从格林威治零度经线起，每6度分为一个投影带，全球共分为60个投影带。东半球的30个投影带的中央经线用 L0=6n-3 计算（n为投影带带号），从0到180度，其编号为1-30。西半球也有30个投影带，从-180度回到0度，其编号为31-60，各带的中央经线用L0=6(n-30)-3-180计算。该投影带将地球划分为60个投影带，每带经差为6度，已被许多国家作为地形图的数学基础。一般从南纬度80到北纬度84度的范围内使用该投影。3度分带法从东经1度30分算起，每3度为一带。这样分带的方法在于使6度带的中央经线均为3度带的中央经线。但是，在标准比例尺图幅编号中，带号是从西经-180度算起，每6度为1带，自西向东1-60。这样，我们国家的高斯带号在标准图幅编号中，要加30，如20带，表示为J50等。</P>
<P>6度分带投影区的代号与其所对应的经度范围如6度分带图表所示。</P>
<P>由于高斯-克吕格投影每一个投影带的坐标都是对本带坐标原点的相对值，所以各带的坐标完全相同，使用时只需变一个带号即可。因此，计算一个带的坐标值，制成一个表，就可以供查取各投影带的坐标时使用，称为高斯坐标表，表中的值成为通用坐标值。</P>
<P>在高斯坐标系中，为了避免横坐标Y有负值，将其起算原点向西移动500公里，即对横坐标Y值按代数法加上500000米。此外，在计算出来的和数前面加上带号，以便识别该点位于何带。例如位于45带之某一点，其横坐标值为Y=-126568.24米，根据上面的规定，改变的(通用的)横坐标值Y=45373431.76米。用户需注意，在MAPGIS中使用时，记着要去掉前边的带号。</P>
<P>注意： </P>
<P>   1.高斯-克吕格坐标系的X、Y轴正好对应MAPGIS坐标系的Y和X，高斯-克吕格坐标系的纵向为X，而MAPGIS坐标系的纵向为Y。</P>
<P>   2.高斯-克吕格坐标系的横向坐标最多为6位，纵向最多为7位。在MAPGIS中使用时，若横向为8位，则前2位为带号，使用时记着要去掉前边的带号，将带号输入对应的参数中。</P>
<P>   3.高斯-克吕格坐标系的坐标单位为米，而MAPGIS坐标系的坐标为毫米，所以输入比例尺时要注意对应。</P>
<P>   4.1:50万以上的标准图框都是高斯-克吕格坐标系，并且生成的标准图框都进行了坐标平移和旋转，使左下角为(0,0)，左下角和右下角的坐标在纵向上相同，即水平对齐。而投影变换中的坐标都是对应投影的大地坐标，因此，在用标准图框进行投影转换前，需要先将其还原为相应的大地坐标，才能开始转换。在后边的标准图框生成过程中，有一个“是否将左下角平移为原点”，若不选该开关，则生成的标准图框中的坐标就为大地坐标，从而可以直接参加投影变换。</P>
<P>   5.在用户输入或矢量化的图中，其用户参考坐标系一般情况下与投影坐标系不重合，因此，用户在将这样的图进行投影转换前，只有输入控制点(TIC点)将其平移、校正到相应的投影坐标系中，才能开始转换，否则结果不正确。总之，投影转换是相对于对应投影坐标系，而非用户坐标系。</P>

怎么没有看到图！