1、GIS概论概论2地理空间数学基础地理空间数学基础地理空间的数学建模 为了深入研究地理空间,需要建立地球表面的几何模型地球表面的几何模型,这是进行大地测量的前提。根据大地测量学的成果,地球表面模型地球表面模型可以分为四类:q 最自然的面最自然的面:包括海洋底部、高山、高原等在内的固体地球表面。太复杂,难以建模,各种量算也非常困难。q 相对抽象的面相对抽象的面:也称为大地水准面,是静止海平面的延伸。以它为基准,可以用水准仪测量地球自然表面上任意点的高程。海平面的起伏将导致测量的不确定。q 模型模型:以大地水准面为基准建立的地球椭球体模型地球椭球体模型。q 其他数学模型其他数学模型:为了解决特定的大
2、地测量问题而提出的。如类地形面、准大地水准面、静态水平衡椭球体等。地球模型:三级近似地球模型:三级近似地球自然表面极不规则,无法用数学表面进行描述水准面所包围的球体大地水准面所包围的球体旋转椭球体不规则性、动态性、不唯一性不规则性、相对唯一性标准数学曲面 1952:海福特椭球 1953:克拉索夫斯基椭球 1978:1975年国际椭球q地球椭球体地球椭球体模型:以大地水准面为基准模型:以大地水准面为基准建立的。建立的。地球的形状接近于椭圆绕其短轴地球的形状接近于椭圆绕其短轴形成的椭球体,通过扁率表示椭球体的扁形成的椭球体,通过扁率表示椭球体的扁平程度。大地水准面与具有微小扁率的旋平程度。大地水准
3、面与具有微小扁率的旋转椭球面非常接近,可用旋转椭球体代替转椭球面非常接近,可用旋转椭球体代替大地球体。大地球体。2、地球空间模型描述 地理空间数学建模地球椭球体模型地球椭球体模型bca三轴椭球体模型双轴椭球体模型(旋转椭球体)其他椭球体模型:根据a、b、c的不同x2a2y2b2z2c2+=1x2a2y2b2z2a2+=1如:克拉索夫斯基椭球体 2、地球空间模型描述 地理空间数学建模椭球体参数长半径 a(赤道半径)短半径 b(极半径)扁 率 =(a-b)/a第一偏心率 e2=(a2-b2)/a2第二偏心率 e2=(a2-b2)/b2我国使用的椭球我国使用的椭球克拉索夫斯基椭球体克拉索夫斯基椭球体
4、IAG75IAG75椭球体椭球体WGS84WGS84椭球体椭球体 我国的大地坐标系和高程系我国的大地坐标系和高程系19541954年北京坐标系年北京坐标系19801980年国家大地坐标系年国家大地坐标系西安原点西安原点GPSGPS测量数据测量数据19561956年黄海高程系年黄海高程系19851985年国家高程基准年国家高程基准中国曾经用过的椭球体我国我国19521952年以前采用年以前采用HayfordHayford椭球体椭球体;19531953年开始采用年开始采用克拉索夫斯基椭球体克拉索夫斯基椭球体;上世纪上世纪7070年代末建立新的年代末建立新的8080坐标系时,坐标系时,采用采用IUG
5、GIUGG(国际大地测量与地球物理联合会)(国际大地测量与地球物理联合会)椭球体椭球体;19841984年定义的年定义的世界大地坐标系世界大地坐标系(WGS84WGS84)使用的椭球体长、短半径则分别为使用的椭球体长、短半径则分别为6378.1376378.137和和6356.75236356.7523,扁率为,扁率为1 1:298.26298.26。直接建立在球体上的地理坐标,用经度和纬度表达地理对象位置建立在平面上的直角坐标系统,用(x,y)表达地理对象位置投影定位坐标系:平面系统定位坐标系:平面系统3.2.2 大地坐标系(1)54年北京坐标系 在东北黑龙江边境上同苏联大地网联测,通过大地
6、坐标计算,推算出北京点的坐标,北京坐标系是苏联42年坐标系的延伸,其原点在苏联普尔科沃。(2)80年西安坐标系 78年4月召开“全国天文大地网平差会议”建立80年西安坐标系,其原点在西安西北的永乐镇,简称西安原点。椭球体参数为75年国际大地测量与地球物理联合会第16界大会的推荐值。(3)新54年北京坐标系 将全国大地网整体平差的结果整体换算到克拉索夫斯基椭球体上,形成一个新的坐标系,称为新54年北京坐标系,它与80年国家大地坐标系的轴定向基准相同,网的点位精度相同。(4)WGS84坐标系 在GPS定位中,定位结果属于WGS84坐标系,坐标系原点位于质心,Z轴指向BIH1984.0协议地极(CT
7、P)。高程系统高程系统 1、定义:地面任一点沿基准线到基准面的距离。某点沿铅锤线方向到达大地水准面的距离称之为该点的绝对高程或海拔高 2、高差:地面两点对应高程之差。3、1956年黄海高程系4、1985国家高程基准定位坐标系:高程系统定位坐标系:高程系统2、地理要素的定位系统为了确定地理要素的位置,需要确定其定位系统:为了确定地理要素的位置,需要确定其定位系统:1、球面定位系统球面定位系统地理坐标地理坐标(L、B)2 2、平面定位系、平面定位系统统现实世界是三现实世界是三维维大地原大地原点点水准原水准原点点?为什么要把球面的转换为平面的?为什么要把球面的转换为平面的?怎么把球面的转换为平面的?
8、怎么把球面的转换为平面的?为什么要进行地图投影?将地球将地球椭球面椭球面上的点映射到上的点映射到平面平面上的方法,上的方法,称为称为地图投影地图投影。1 1、地理坐标为球面坐标,不方便进行距离、地理坐标为球面坐标,不方便进行距离、方位、面积等参数的量算。方位、面积等参数的量算。2 2、地球椭球体为不可展曲面。、地球椭球体为不可展曲面。3 3、地图为平面,符合视觉心理,并易于进、地图为平面,符合视觉心理,并易于进行距离、方位、面积等量算和各种空间分析。行距离、方位、面积等量算和各种空间分析。地理空间数学基础q地球空间参考 地球椭球体、大地水准面、投影q 空间数据投影 地图投影的分类、常用地图投影
9、简介 q空间坐标转换 q 空间尺度 比例尺、分辨率 q 地理格网 地理格网标准、区域划分标准、国家基本比例尺 地形图标准3、地图投影、什么是地图投影?、什么是地图投影?、坐标系、坐标系、地图投影的分类、地图投影的分类、如何区分等角、等面积、任意投影?、如何区分等角、等面积、任意投影?、中国地图投影的采用方式、中国地图投影的采用方式、常见的地图投影、常见的地图投影什么是地图投影?简单地讲:地图投影的实质是将地球椭球面上的简单地讲:地图投影的实质是将地球椭球面上的经纬网按照一定的数学法则转移到平面上。经纬网按照一定的数学法则转移到平面上。具体来说:由于球面上一点的位置是用地理坐标具体来说:由于球面
10、上一点的位置是用地理坐标表示,而平面上是用直角坐标或者极坐标表示,所以表示,而平面上是用直角坐标或者极坐标表示,所以要想将地球表面上的点转移到平面上,必须采用一定要想将地球表面上的点转移到平面上,必须采用一定的数学方法来确定地理坐标与平面直角坐标或极坐标的数学方法来确定地理坐标与平面直角坐标或极坐标之间的关系。之间的关系。这种在球面和平面之间建立点与点之间函数关系这种在球面和平面之间建立点与点之间函数关系的数学方法,称为的数学方法,称为地图投影地图投影。坐标系选定了一个一定大小的椭球体,并确定了它与选定了一个一定大小的椭球体,并确定了它与大地水准面的相关位置,就确定了一个坐标系。大地水准面的相
11、关位置,就确定了一个坐标系。地图投影的目的地图投影的目的:地图投影的目的:是通过将不可展的球面投影到一个可是通过将不可展的球面投影到一个可展曲面上,然后将该曲面展开成为一个平展曲面上,然后将该曲面展开成为一个平面,来保证面,来保证空间信息空间信息在地域上的在地域上的连续性连续性、完整性完整性和和可测度性可测度性。地图投影的分类 根据美国著名地图投影学家根据美国著名地图投影学家J.P.SnyderJ.P.Snyder统计,统计,全世界地图投影种类现有全世界地图投影种类现有256256种种,依据不同的目的和,依据不同的目的和要求,可以要求,可以采用不同的分类指标采用不同的分类指标对如此繁多的地图对
12、如此繁多的地图投影进行分类。投影进行分类。分类分类1 1:基于投影面与球面相关位置的分类;:基于投影面与球面相关位置的分类;分类分类2 2:基于投影方法的分类;:基于投影方法的分类;分类分类3 3:基于投影方程的分类;:基于投影方程的分类;分类分类4 4:基于投影变形的分类。:基于投影变形的分类。基于投影面与球面的分类4、地图投影、地图投影地图投影的变形地图投影的变形 用地图投影的方法将球面展开为平面,虽用地图投影的方法将球面展开为平面,虽然可以保持地域上的联系和完整性,但它们与然可以保持地域上的联系和完整性,但它们与球面上的经纬度网线形状并不一致。即投影后,球面上的经纬度网线形状并不一致。即
13、投影后,地图上的经纬度网线发生了变形,同样根据地地图上的经纬度网线发生了变形,同样根据地理坐标展绘在地图上的各种要素,也必然随着理坐标展绘在地图上的各种要素,也必然随着变形。变形。4、地图投影、地图投影地图投影的变形地图投影的变形这种变形使得地理要素的几何特性受到破坏:这种变形使得地理要素的几何特性受到破坏:q 长度变形:地球仪上,纬线长度不等;同一纬长度变形:地球仪上,纬线长度不等;同一纬线上,经差相同,纬线长度相同;同一经线上,线上,经差相同,纬线长度相同;同一经线上,纬差相同而经线长度不同;所有经线长度相等。纬差相同而经线长度不同;所有经线长度相等。q 面积变形:地球仪上,同一纬度带内,
14、经差相面积变形:地球仪上,同一纬度带内,经差相同的网格面积相等;同一经度带内,纬度越高,同的网格面积相等;同一经度带内,纬度越高,面积越小。面积越小。q 角度变形:地球仪上,经线与纬线处处呈直角角度变形:地球仪上,经线与纬线处处呈直角相交。相交。地图投影变形的图解示例(摩尔维特投影等积伪圆柱投影)长度变形角度变形地图投影变形的图解示例(UTM横轴等角割圆柱投影)面积变形和长度变形地图投影的变形示意地图投影的变形示意4、地图投影、地图投影地图投影的变形地图投影的变形区分等角、等面积、任意投影形状不变面积不变特定方向距离不变地图投影地图投影地图投影的分类地图投影的分类圆柱投影方位投影圆锥投影正轴切
15、圆锥投影 正轴割圆锥投影横轴切圆锥投影 横轴割圆锥投影横轴切圆柱投影 横方位投影正轴割圆柱投影 斜轴切圆柱投影斜轴切圆锥投影 正轴切圆柱投影正方位投影 斜方位投影中国地图投影系统 采用与采用与我国基本图系列我国基本图系列一致一致的地图投影系统:即大比例尺时的高斯的地图投影系统:即大比例尺时的高斯-克克吕格投影(横轴等角切椭圆柱投影)和中小比例尺时的吕格投影(横轴等角切椭圆柱投影)和中小比例尺时的LambertLambert投影(正轴投影(正轴等角割圆锥投影)。我国常用的地图投影的情况为:等角割圆锥投影)。我国常用的地图投影的情况为:1)1)、我国基本比例尺地形图、我国基本比例尺地形图(1(1:
16、100100万、万、1 1:5050万、万、1 1:2525万、万、1 1:1010万、万、1 1:5 5万、万、1 1:2.52.5万、万、1 1:1 1万、万、1 1:5000),5000),除除1 1:100100万外均万外均采用采用高斯高斯克吕格投影克吕格投影为地理基础;为地理基础;2)2)、我国、我国1 1:100100万地形图采用了万地形图采用了LambertLambert投影投影,其分幅原则与,其分幅原则与国际地理学会规定的国际地理学会规定的全球统一使用的国际百万分之一全球统一使用的国际百万分之一地图投影地图投影保持一致。保持一致。3)3)、我国大部分、我国大部分省区图省区图以及
17、大多数以及大多数这一比例尺这一比例尺的地图也多采用的地图也多采用LambertLambert投影和属于同一投影系统的投影和属于同一投影系统的AlbersAlbers投影投影(正轴等面积割正轴等面积割圆锥投影圆锥投影);4)4)、LambertLambert投影中,地球表面上两点间的投影中,地球表面上两点间的最短距离最短距离(即大圆航即大圆航线线)表现为近于直线,这表现为近于直线,这有利于有利于地理信息系统中地理信息系统中空间分析量度空间分析量度的的正确实施。正确实施。对于大中比例尺地图,一般来说大多数都采用地形图的数学基础高斯克吕格投影,尤其是当比例尺为国家基本地形图比例尺系列时,可直接判定为
18、高斯克吕格投影。其原因是,这些比例尺和基本地形图比例尺相一致,编图时,选用地形图的数学基础,既免去了重新展绘数学基础的工序,而且能够保持很高的点位精度。大中比例尺地大中比例尺地图图u 高斯克吕格投影(Gauss-Kruger Projection)高斯克吕格投影是由高斯于19世纪20年代拟定,后经克吕格补充而形成的一种地图投影方式。在英美国家称为横轴墨卡托投影 属于横轴等角切圆柱投影。这种投影是将椭圆柱面套在地球椭球的外面,并与某一子午线相切(此子午线叫中央子午线或中央经线),椭圆柱的中心轴通过地球椭球的中心,然后用等角条件将中央子午线东西两侧各一定经差范围内的地区投影到柱面上,并将此柱面展成
19、平面,即获得高斯投影横轴圆柱投影横轴圆柱投影u 高斯高斯克吕格投影克吕格投影(Gauss-Kruger Projection)xy高斯-克吕格投影原理图u 高斯克吕格投影(Gauss-Kruger Projection)高斯投影特征:高斯投影特征:中央经线和赤道投影为互相垂直的直线,且为投影中央经线和赤道投影为互相垂直的直线,且为投影的对称轴的对称轴 投影后无角度变形,即保角投影投影后无角度变形,即保角投影 中央经线无长度变形中央经线无长度变形 同一条经线上,纬度越低,变形越大,赤道处最同一条经线上,纬度越低,变形越大,赤道处最大大同一条纬线上,离中央经线越远,变形越大;同一条纬线上,离中央经
20、线越远,变形越大;为了保证地图的精度,采用分带投影方法,即将投为了保证地图的精度,采用分带投影方法,即将投影范围的东西界加以限制,使其变形不超过一定的限影范围的东西界加以限制,使其变形不超过一定的限度,这样把许多带结合起来,可成为整个区域的投影度,这样把许多带结合起来,可成为整个区域的投影 在在6 6带带范围内,长度变形线最大不超过范围内,长度变形线最大不超过0.14%0.14%3度带和度带和6度带度带从从0 0度开始,自西向东每度开始,自西向东每6 6度分为一个投影带。度分为一个投影带。从东经从东经1 1度度3030分开始,自西向东每分开始,自西向东每3 3度分为一个投影带。度分为一个投影带
21、。u 高斯高斯克吕格投影克吕格投影(Gauss-Kruger Projection)-投影分带投影分带 我国我国1:11:1万至万至1:501:50万的地形图全部采用高斯万的地形图全部采用高斯克吕格投影。克吕格投影。1:2.51:2.5万至万至1:501:50万的地形图,采用万的地形图,采用6 6分带方案,全球共分为分带方案,全球共分为6060个投影带;我国位于个投影带;我国位于东经东经7272到到136136间,共含间,共含1111个投影带;个投影带;1:11:1万比万比例尺图采用例尺图采用3 3分带方案,全球共分带方案,全球共120120个带。个带。分割条带号规定:从分割条带号规定:从0
22、0子午线开始分子午线开始分6 6经度为一带,经度为一带,东半球东经东半球东经3 3、9 9、1515177177分别是分别是1 1、2 2、3 33030条条6 6带带的中央子午线,然后继续自西向东旋转,每转的中央子午线,然后继续自西向东旋转,每转6 6增加带号增加带号1 1。分割分割3 3带原则上与带原则上与6 6带相同,只是从东经带相同,只是从东经1 13030(即(即1.51.5E E)起,每隔)起,每隔3 3带为带为1 1个投影带。个投影带。u 高斯高斯克吕格投影克吕格投影(Gauss-Kruger Projection)-投影分带投影分带 在高斯克吕格投影上,规定以中央经线为在高斯克
23、吕格投影上,规定以中央经线为X X轴,轴,赤道为赤道为Y Y轴,两轴的交点为坐标原点。轴,两轴的交点为坐标原点。X X坐标值在赤道以北为正,以南为负;坐标值在赤道以北为正,以南为负;Y Y坐标坐标值在中央经线以东为正,以西为负。我国在北半球,值在中央经线以东为正,以西为负。我国在北半球,X X坐标皆为正值。坐标皆为正值。Y Y坐标在中央经线以西为负值,运坐标在中央经线以西为负值,运用起来很不方便。为了避免用起来很不方便。为了避免Y Y坐标出现负值,通常坐标出现负值,通常将各带的坐标纵轴西移将各带的坐标纵轴西移500500公里,即将所有公里,即将所有Y Y值都加值都加500500公里。公里。u
24、高斯高斯克吕格投影克吕格投影(Gauss-Kruger Projection)-投影分带投影分带u正轴等角割圆锥投影(Lambert投影)这种投影是将一圆锥面套在地球椭球外面,将地球表面上的要素投影到圆锥面上,然后将圆锥面沿某一母线(经线)展开,即获得Lambert投影。这种投影中,经线为交于一点这种投影中,经线为交于一点的直线束,纬线为同心圆圆弧,的直线束,纬线为同心圆圆弧,圆心即直线束的交点圆心即直线束的交点经线呈辐射经线呈辐射状,为纵向直线,纬线近似于弧状,为纵向直线,纬线近似于弧形,与经线正交形,与经线正交适用于适用于1 1:100100万(包括万(包括1 1:100100万)以上地形
25、图万)以上地形图正轴圆锥投影 Lambert投影(正轴等角割圆锥投影)。u正轴割圆锥投影(Lambert投影)误差情况:圆锥与地球相交处为北纬25与北纬47,距离误差随地点纬度不同而不同,在成图范围内北部最大达+4%,南部达3%,中部为-1.8%,面积变形相对误差相比距离相对误差要大一倍。一幅图可覆盖大片中纬度地区,可整幅覆盖我一幅图可覆盖大片中纬度地区,可整幅覆盖我国境内领土;国境内领土;地球表面上两点间的最短距离表现为近于直线地球表面上两点间的最短距离表现为近于直线,这有利于地理信息系统中的空间分析和信息量度,这有利于地理信息系统中的空间分析和信息量度的正确实施的正确实施 小比例尺地图 由
26、于我国位于中纬度地区,中国地图和分省地图经常采用割圆锥投影(Lambert 或Albers 投影),中国地图的中央经线常位于东经105度(110度),两条标准纬线分别为北纬25度和北纬47度,而各省的参数可根据地理位置和轮廓形状初步加以判定。例如甘肃省的参数为:中央经线为东经101度,两条标准纬线分别为北纬34度和41度。地图投影的选择地图投影的选择 地图投影将直接影响地图的精度和使用价值。地图投影将直接影响地图的精度和使用价值。通常地图投影对中小比例尺地图影响很大,对于大通常地图投影对中小比例尺地图影响很大,对于大比例尺地图,则影响很小。一般国家基本比例尺地比例尺地图,则影响很小。一般国家基
27、本比例尺地形图的地图投影选择是由国家测绘部门制订,不允形图的地图投影选择是由国家测绘部门制订,不允许随便更改。许随便更改。地图投影的选择主要考虑以下因素:制图区地图投影的选择主要考虑以下因素:制图区域的范围、形状和地理位置;地图的用途、出版方域的范围、形状和地理位置;地图的用途、出版方式及其他要求等。式及其他要求等。地理坐标地图(球面坐标表示为平面方式)等积圆柱投影等积圆柱投影墨卡托投影:正轴等角切圆柱投影编制世界时区图横轴墨卡托投影横轴墨卡托投影横轴墨卡托投影摩尔维特(Mollweide)投影等面积伪圆柱投影,常用于编制小比例尺世界地图等面积伪圆柱投影,常用于编制小比例尺世界地图常见的地图投
28、影 世界图全国图:正轴圆锥投影全国图:正轴圆锥投影地形图:高斯克吕格投影(分带)地形图:高斯克吕格投影(分带)南北半球(或两极)图:正轴方位投影、南北半球(或两极)图:正轴方位投影、亚洲图:斜轴方位投影亚洲图:斜轴方位投影 半球图东西半球图:横轴方位投影东西半球图:横轴方位投影欧洲图:彭纳投影欧洲图:彭纳投影国内出版:等差分纬线多圆锥投影国内出版:等差分纬线多圆锥投影国外出版:摩尔威特投影国外出版:摩尔威特投影地球仪:地球仪:普通多圆普通多圆锥投影锥投影海海 图:墨卡托投影图:墨卡托投影中国中国图图大洲大洲图图GIS与地图投影的关系 地理基础地理基础(地图投影地图投影)地图投影对地图投影对GI
29、SGIS的影响渗透在的影响渗透在GISGIS建设的各个阶段和各个方建设的各个阶段和各个方面。面。(据边馥苓,(据边馥苓,1996)数据获取数据获取(数据源地图的投影)(数据源地图的投影)数据标准化预处理数据标准化预处理(按某一参照系数字化)(按某一参照系数字化)数据存储数据存储(统一的坐标基础)(统一的坐标基础)数据处理数据处理(投影转换)(投影转换)数据应用数据应用(检索查询、覆盖分析等)(检索查询、覆盖分析等)数据输出数据输出(具有相应投影的地图)(具有相应投影的地图)大地原点水准原点水准原点水准原点19851985国家高国家高程基准,程基准,72.260472.2604米米黄海海面黄海海
30、面1952-19791952-1979年平年平均海水面为均海水面为0 0米米地图投影与坐标系的关系地地图图投投影影椭椭球球体体模模型型Back国家坐标系和独立坐标系的变换 由于地球半径很大,在较小区域内进行测量工作可将地球椭球面作为平面看待,而不失其严密性。既然把投影基准面作为平面,就可采用平面直角坐标系表示地面点的投影面上的位置。(a)测量平面直角坐标系)测量平面直角坐标系(b)数学平面直角坐标系)数学平面直角坐标系 为不使坐标系出现负值,它通常将某测区的坐标原点设在测区西南角某点,以真北方向或主要建筑物主轴线为纵轴方向,而以垂直于纵坐标轴的直线定为横坐标轴,构成平面直角坐标系;也可假设测区
31、中某点的坐标值,以该点到另一点方位角作为推算其它各点的起算数据,实际上也构成了一个平面直角坐标系。上述平面直角坐标系的原点和纵轴方向选定了的值常用于小型测区的测量,它不与国家统一坐标系相连,因此称为任意坐标系或独立坐标系。我国大部分城市均采用独立坐标系,如广州市采用珠江高程和平面坐标系等。国家坐标系和独立坐标系的变换国家坐标系和独立坐标系的变换 按高斯投影统一分带(按高斯投影统一分带(6 0带,带,3 0带)建立带)建立的直角坐标系,称为国家平面直角坐标系。的直角坐标系,称为国家平面直角坐标系。在建立数字城市时,往往需要将独立坐标在建立数字城市时,往往需要将独立坐标系转换成国家平面直角坐标系。
32、在进行转换时,系转换成国家平面直角坐标系。在进行转换时,先将独立坐标系的原点或独立坐标系的某一固定先将独立坐标系的原点或独立坐标系的某一固定点与国家大地点连测,并按计算出的方位角进行点与国家大地点连测,并按计算出的方位角进行改正,求出该点的国家统一坐标,然后对所有数改正,求出该点的国家统一坐标,然后对所有数据进行平移和旋转,以便把按独立坐标系所采集据进行平移和旋转,以便把按独立坐标系所采集的数据转换到国家平面直角坐标系中。在城市和的数据转换到国家平面直角坐标系中。在城市和工程测量中,也可采用工程测量中,也可采用1.5 0带或任意带的高斯平带或任意带的高斯平面坐标系,以提高投影的精度。面坐标系,
33、以提高投影的精度。地理空间数学基础q地球空间参考 地球椭球体、大地水准面、投影q 空间数据投影 地图投影的分类、常用地图投影简介 q空间坐标转换 q 空间尺度 比例尺、分辨率 q 地理格网 地理格网标准、区域划分标准、国家基本比例尺 地形图标准反解变换法(同一坐标基准下的坐标变换)通过中间过渡的方法,反解出原投影点的地理坐标,代入新投影中求得新投影之坐标。例1 从Albers投影转换到Transvers投影不同坐标基准下的坐标变换1、地理坐标基准的变换2、坐标值的变换地理空间数学基础q地球空间参考 地球椭球体、大地水准面、投影q 空间数据投影 地图投影的分类、常用地图投影简介 q空间坐标转换
34、q 空间尺度 比例尺、分辨率 q 地理格网 地理格网标准、区域划分标准、国家基本比例尺 地形图标准地理空间数学基础q地球空间参考 地球椭球体、大地水准面、投影q 空间数据投影 地图投影的分类、常用地图投影简介 q空间坐标转换 q 空间尺度 比例尺、分辨率 q 地理格网 地理格网标准、区域划分标准、国家基本比例尺 地形图标准四.地形图分幅与编号1.1:100万地形图的分幅编号万地形图的分幅编号l经差经差6度、纬差度、纬差4度;度;l纬度纬度60度度76度为经差度为经差12度、纬差度、纬差4度;度;l纬度纬度76度度88度之间经差度之间经差24度、纬差度、纬差4度。度。l分幅编号分幅编号:l从经度
35、从经度180度起自西向东每经差度起自西向东每经差6度为一列,依次用阿度为一列,依次用阿拉伯数字拉伯数字1,2,3,,表示相应列号;,表示相应列号;l 从赤道算起,每纬差从赤道算起,每纬差4度为一行,至南、北纬度为一行,至南、北纬88度各分度各分为为22行,依次用大写字母行,依次用大写字母A、B、C、V表示相应行号。表示相应行号。l l由经线和纬线所围成的每个梯形小格为一幅由经线和纬线所围成的每个梯形小格为一幅1:100万地形图,编号由其所在的行号和列号组合而万地形图,编号由其所在的行号和列号组合而成。如北京为成。如北京为J 50。l注意:注意:l(1)6 6度投影带的分带号与百万分之一地形图度
36、投影带的分带号与百万分之一地形图分幅的纵行相差分幅的纵行相差30。l(2)我国小于百万分之一的地图(含百万分之)我国小于百万分之一的地图(含百万分之一)用等角圆锥投影,大于一)用等角圆锥投影,大于50万分之一的地图都万分之一的地图都用高斯克吕格投影。用高斯克吕格投影。1、1/1000 000地图的分幅编号地图的分幅编号 旧地形图编号的一种:旧地形图编号的一种:1100万:万:J-50 150万:万:在在1100万上分万上分4幅幅 J-50-A 125万:万:在在1100万上分万上分16幅幅 J-50-A-a 120万:万:在在1100万上分万上分36幅幅 J-50-A-1 110万:万:在在1
37、100万上分万上分144幅幅 J-50-144 15万:万:在在110万上分万上分4幅幅 J-50-144-A 12.5万:万:在在110万上分万上分16幅幅 J-50-144-A-10 11万:万:在在110万上分万上分64幅幅 J-50-144-A-1 15千:千:在在110万上分万上分256幅幅 J-50-144-A-2002.20世纪世纪90年代以前地形图编号(不规范)年代以前地形图编号(不规范)2、1/500 000、1/200 000 和1/100 000地图的分幅编号1/500 000地图1/200 000地图 1/100 000地图(分幅编号示例)3.20世纪世纪90年代后地形
38、图编号年代后地形图编号1:50万万1:5千地形图的编号均以千地形图的编号均以1:100万地形图编号为基万地形图编号为基础础 组成组成:其所在:其所在1:100万地形图的图号、比例尺代码万地形图的图号、比例尺代码和各图幅的行列号(共和各图幅的行列号(共10位)位)比例尺比例尺1:5000001:2500001:1000001:500001:250001:100001:5000代码代码BCDEFGH比例尺代码比例尺代码行列号用三位阿拉伯数字表示,不足前面补零,行列号用三位阿拉伯数字表示,不足前面补零,行号在前、列号在后行号在前、列号在后 例例1.1:50万地形图编号万地形图编号 J50B00100
39、1;J50B001002;J50B002001;J50B002002。3638040117120001001002001002002001002J50例例2 1:25万地形图编号万地形图编号 J50C001001;J50C001002;J50C001003;J50C001004;J50C004001;J50C0040043837391153011830004002004001004003004004003001003002003004003003J50001001001003001004001002002001002002002004002003例例3 1:10万地形图编号万地形图编号J50D
40、001001;J50D001002.J50D012001;J50D012002.J50D012011;J50D0120123612011711440393837经 差6纬差4例例4 1:5万地形图编号万地形图编号J50E001001;J50E001002.J50E024001;J50E024002.J50E024023;J50E024024三.编号应用的公式 1、已知图幅某点的经、纬度或图幅西南图廓点的经纬度,计算其编号1:100万地形图图幅编号计算公式 a=/4+1 b=/6+31 其中:其中:-表示商取整;表示商取整;a-1:100万地形图图幅所在纬度带字符码所对应的数字码;万地形图图幅所
41、在纬度带字符码所对应的数字码;b-1:100万地形图图幅所在经度带对应的数字码;万地形图图幅所在经度带对应的数字码;-图幅内某点的经度或图幅西南图廓点的经度;图幅内某点的经度或图幅西南图廓点的经度;-图幅内某点的纬度或图幅西南图廓点的纬度。图幅内某点的纬度或图幅西南图廓点的纬度。例1:某点经度为1143345,纬度为392230,计算其所在的图幅的编号。a=392230/4+1=10(字母J)b=1143345/6+3150 该点在1:100万地形图图号为J-50。其它比例尺地形图在1:100万地形图图号后的行列号 c=4/-(/4)/d=(/6)/+1 其中:()-表示商取余;-表示商取整;
42、c -所求比例尺地形图在1:100万地形图图号后的行号;d -所求比例尺地形图在1:100万地形图图号后的列号;-图幅内某点的经度或图幅西南图廓点的经度;-图幅内某点的纬度或图幅西南图廓点的纬度;-所求比例尺地形图分幅的经差;-所求比例尺地形图分幅的纬差。1:501:50万:万:2 2 c=4 c=4/2/2-(39392230/42230/4)/2/2 =223022302 2001001 d=d=(1141143345/63345/6)/+1+1 =333345/45/+1+1001001所以图号为所以图号为J50B001001J50B001001例例2 经度经度1143345,纬度,纬度
43、392230计算各比例尺地形计算各比例尺地形图的编号。图的编号。1:251:25万万:1 1 1 13030 c=4 c=4/1/1-(39392230/42230/4)/1/1 =4 =4223022301 1001001 d=d=(1141143345/63345/6)/1/130+1 30+1 =33 =3345/145/13030+1=001+1=0011:251:25万地形图的图号为万地形图的图号为J50C001001J50C001001 1:101:10万万:20 20 3030 c=c=4 4/20-/20-(39392230/42230/4)/20/20=12=12223022
44、302020=002=002 d=d=(1141143345/63345/6)/30+1/30+1 =333345/3045/30+1=002+1=002 1:101:10万地形图的图号为万地形图的图号为J50D002002J50D002002 1:51:5万万:10 10 1515 c=4c=4/10-/10-(39392230/42230/4)/10/10=24=24223022301010=004=004 d=d=(1141143345/63345/6)/15+1/15+1 =333345/1545/15+1=003+1=003 1:51:5万地形图的图号为万地形图的图号为J50E004
45、003J50E0040031:2.51:2.5万万:5 5 730730 c=4 c=4/5-/5-(39392230/42230/4)/5/5 =48 =482230223055 =008 =008 d=d=(1141143345/63345/6)/730+1/730+1 =333345/73045/730+1+1 =005 =0051:2.51:2.5万地形图的图号为万地形图的图号为J50F008005J50F0080051:11:1万万:230 230 345345 c=4 c=4/230-/230-(39392230/42230/4)/230/230 =96 =962230223023
46、0230=015=015 d=d=(1141143345/63345/6)/345+1/345+1 =333345/34545/345+1=010+1=0101:11:1万地形图的图号为万地形图的图号为J50G015010J50G0150101:51:5千千:115 115 152.5152.5 c=c=4 4/115-/115-(39392230/42230/4)/115/115 =192 =19222302230115115=030=030 d=d=(1141143345/63345/6)/152.5+1/152.5+1 =333345/152.545/152.5+1=019+1=0191
47、:51:5千地形图的图号为千地形图的图号为J50H030019J50H0300192、已知图号计算该图幅西南图廓点 的经、纬度=(b-31)6+(d-1)=(a-1)4+(4/-c )式中:图幅西南图廓点的经度;图幅西南图廓点的纬度。a 1:100万地形图图幅所在纬度带字符码所对应的数字码;b 1:100万地形图图幅所在经度带的数字码;c 该比例尺地形图在1:100万地形图图号后的行号;d 该比例尺地形图在1:100万地形图图号后的列号;所求比例尺地形图分幅的经差;所求比例尺地形图分幅的纬差。例例1:图号为:图号为J50B001001,计算其西南图廓点的经、,计算其西南图廓点的经、纬度。纬度。
48、a=10;b=50;c=1;d=1;=2;=3=(5031)6+(1 1)3=114=(101)4+(4/2 1)2 =38 该图幅西南图廓点的经、纬度分别为该图幅西南图廓点的经、纬度分别为114、38例例2:图号为:图号为J50D002002,计算其西南图廓点的经、,计算其西南图廓点的经、纬度。纬度。a=10;b=50;c=2;d=2;=20;=30=(5031)6+(21)30=11430=(101)4+(4/202)20=3920该图幅西南图廓点的经、纬度分别为该图幅西南图廓点的经、纬度分别为11430、3920A:西北角图幅的行、列号按下式计算:西北角图幅的行、列号按下式计算:c大大=
49、小小/大大(c小小-1)+1 d大大=小小/大大(d小小-1)+13、在同一幅、在同一幅1:100万地形图图幅内万地形图图幅内3837391153011830004002004001004003004004003001003002003004003003J50001001001003001004001002002001002002002004002003(1)小比例尺地形图的行、列号计算大小比例尺地形图的行、列号计算大比例尺地形图的行、列号比例尺地形图的行、列号B:东南角图幅的行、列号按下式计算:东南角图幅的行、列号按下式计算:c 大大=c小小小小/大大 d 大大=d小小小小/大大c大大 -较
50、大比例尺地形图在较大比例尺地形图在1:100万地形图图号后的行号;万地形图图号后的行号;d大大 -较大比例尺地形图在较大比例尺地形图在1:100万地形图图号后的列号;万地形图图号后的列号;c小小 -较小比例尺地形图在较小比例尺地形图在1:100万地形图图号后的行号;万地形图图号后的行号;d小小 -较小比例尺地形图在较小比例尺地形图在1:100万地形图图号后的列号;万地形图图号后的列号;大大-较大比例尺地形图分幅的纬差;较大比例尺地形图分幅的纬差;小小-较小比例尺地形图分幅的纬差。较小比例尺地形图分幅的纬差。例例 1:10万地形图的行、列号为万地形图的行、列号为004001,求所含,求所含1:2