GIS概论2地理空间数学基础,GIS概论2地理空间数学基础,地理空间的数学建模,为了深入研究地理空间,需要建立,地球表面的几何模型,,这是进行大地测量的前提。根据大地测量学的成果,,地球表面模型,可以分为四类:,最自然的面,:包括海洋底部、高山、高原等在内的固体地球表面。,太复杂,难以建模,各种量算也非常困难。,相对抽象的面,:也称为大地水准面,是静止海平面的延伸。以它为基准,可以用水准仪测量地球自然表面上任意点的高程。,海平面的起伏将导致测量的不确定。,模型,:以大地水准面为基准建立的,地球椭球体模型,。,其他数学模型,:为了解决特定的大地测量问题而提出的。如类地形面、准大地水准面、静态水平衡椭球体等。,地理空间的数学建模 为了深入研究地理空间,需要建立地球,地球模型:三级近似,地球自然表面,极不规则,无法用数学表面进行描述,水准面所包围的球体,大地水准面所包围的球体,旋转椭球体,不规则性、动态性、不唯一性,不规则性、相对唯一性,标准数学曲面,1952,:海福特椭球,1953,:克拉索夫斯基椭球,1978,:,1975,年国际椭球,地球模型:三级近似地球自然表面极不规则,无法用数学表面进行描,地球椭球体,模型:以大地水准面为基准建立的。,地球的形状接近于椭圆绕其短轴形成的椭球体,通过扁率表示椭球体的扁平程度。大地水准面与具有微小扁率的旋转椭球面非常接近,可用旋转椭球体代替大地球体。,2、地球空间模型描述,地理空间数学建模,2、地球空间模型描述 地理空间数学建模,地球椭球体模型,b,c,a,三轴椭球体模型,双轴椭球体模型(旋转椭球体),其他椭球体模型:根据a、b、c的不同,x,2,a,2,y,2,b,2,z,2,c,2,+,+,= 1,x,2,a,2,y,2,b,2,z,2,a,2,+,+,= 1,如:,克拉索夫斯基椭球体,2、地球空间模型描述,地理空间数学建模,地球椭球体模型bca三轴椭球体模型双轴椭球体模型(旋转椭球体,椭球体参数,长半径 a(赤道半径),短半径 b(极半径),扁 率,=(a-b)/a,第一偏心率 e,2,=(a,2,-,b,2,)/ a,2,第二偏心率 e,2,=(a,2,-,b,2,)/ b,2,我国使用的椭球,克拉索夫斯基椭球体,IAG75椭球体,WGS84椭球体,我国的大地坐标系和高程系,1954年北京坐标系,1980年国家大地坐标系,西安原点,GPS测量数据,1956年黄海高程系,1985年国家高程基准,椭球体参数长半径 a(赤道半径)我国使用的椭球我,中国曾经用过的椭球体,我国1952年以前采用,Hayford椭球体,;,1953年开始采用,克拉索夫斯基椭球体,;,上世纪70年代末建立新的80坐标系时,采用,IUGG(国际大地测量与地球物理联合会)椭球体,;,1984年定义的,世界大地坐标系,(,WGS84,)使用的椭球体长、短半径则分别为6378.137和6356.7523,扁率为1:298.26。,中国曾经用过的椭球体我国1952年以前采用Hayford椭,直接建立在球体上的,地理坐标,用经度和纬度表达地理对象位置,建立在平面上,的直角坐标系统,用(,x,y,)表达地理对象位置,投影,定位坐标系:平面系统,直接建立在球体上的地理坐标,用经度和纬度表达地理对象位置建立,3.2.2 大地坐标系,(1)54年北京坐标系,在东北黑龙江边境上同苏联大地网联测,通过大地坐标计算,推算出北京点的坐标,北京坐标系是苏联42年坐标系的延伸,其原点在苏联普尔科沃。,(2)80年西安坐标系,78年4月召开,“,全国天文大地网平差会议,”,建立80年西安坐标系,其原点在西安西北的永乐镇,简称西安原点。椭球体参数为75年国际大地测量与地球物理联合会第16界大会的推荐值。,(3)新54年北京坐标系,将全国大地网整体平差的结果整体换算到克拉索夫斯基椭球体上,形成一个新的坐标系,称为新54年北京坐标系,它与80年国家大地坐标系的轴定向基准相同,网的点位精度相同。,(4)WGS84坐标系,在GPS定位中,定位结果属于WGS84坐标系,坐标系原点位于质心,Z轴指向BIH1984.0协议地极(CTP)。,3.2.2 大地坐标系(1)54年北京坐标系,高程系统,1,、定义:地面任一点沿基准线到基准面的距离。某点沿铅锤线方向到达大地水准面的距离称之为该点的绝对高程或海拔高,2、高差:地面两点对应高程之差。,3、1956年黄海高程系,4、1985国家高程基准,定位坐标系:高程系统,高程系统3、1956年黄海高程系定位坐标系:高程系统,2、地理要素的定位系统,为了确定地理要素的位置,需要确定其定位系统:,1、,球面定位系统,地理坐标,(L、B),2、平面定位系统,现实世界是三维,大地原点,水准原点,?,为什么要把球面的转换为平面的?,怎么把球面的转换为平面的?,2、地理要素的定位系统为了确定地理要素的位置,需要确定其定位,为什么要进行地图投影?,将地球,椭球面,上的点映射到,平面,上的方法,称为,地图投影,。,1、地理坐标为球面坐标,不方便进行距离、方位、面积等参数的量算。,2、地球椭球体为不可展曲面。,3、地图为平面,符合视觉心理,并易于进行距离、方位、面积等量算和各种空间分析。,为什么要进行地图投影?将地球椭球面上的点映射到平面上的方法,地理空间数学基础,地球空间参考,地球椭球体、大地水准面、投影,空间数据投影,地图投影的分类、常用地图投影简介,空间坐标转换,空间尺度,比例尺、分辨率,地理格网,地理格网标准、区域划分标准、国家基本比例尺 地形图标准,地理空间数学基础地球空间参考,3、地图投影,、什么是地图投影?,、坐标系,、地图投影的分类,、如何区分等角、等面积、任意投影?,、中国地图投影的采用方式,、常见的地图投影,3、地图投影、什么是地图投影?,什么是地图投影?,简单地讲:地图投影的实质是将地球椭球面上的经纬网按照一定的数学法则转移到平面上。,具体来说:由于球面上一点的位置是用地理坐标表示,而平面上是用直角坐标或者极坐标表示,所以要想将地球表面上的点转移到平面上,必须采用一定的数学方法来确定地理坐标与平面直角坐标或极坐标之间的关系。,这种在球面和平面之间建立点与点之间函数关系的数学方法,称为,地图投影,。,什么是地图投影? 简单地讲:地图投影的实质是将地球椭球面,坐标系,选定了一个一定大小的椭球体,并确定了它与大地水准面的相关位置,就确定了一个坐标系。,坐标系选定了一个一定大小的椭球体,并确定了它与大地水准面的相,地图投影的目的,地图投影的目的:,是通过将不可展的球面投影到一个可展曲面上,然后将该曲面展开成为一个平面,来保证,空间信息,在地域上的,连续性,、,完整性,和,可测度性,。,地图投影的目的地图投影的目的:,地图投影的分类,根据美国著名地图投影学家J.P.Snyder统计,全世界地图投影种类现有,256种,,依据不同的目的和要求,可以,采用不同的分类指标,对如此繁多的地图投影进行分类。,分类1:基于投影面与球面相关位置的分类;,分类2:基于投影方法的分类;,分类3:基于投影方程的分类;,分类4:基于投影变形的分类。,地图投影的分类 根据美国著名地图投影学家J.P.Sny,基于投影面与球面的分类,基于投影面与球面的分类,4、地图投影,地图投影的变形,用地图投影的方法将球面展开为平面,虽然可以保持地域上的联系和完整性,但它们与球面上的经纬度网线形状并不一致。即投影后,地图上的经纬度网线发生了变形,同样根据地理坐标展绘在地图上的各种要素,也必然随着变形。,4、地图投影地图投影的变形 用地图投影的方法将球面,4、地图投影,地图投影的变形,这种变形使得地理要素的几何特性受到破坏:,长度变形:地球仪上,纬线长度不等;同一纬线上,经差相同,纬线长度相同;同一经线上,纬差相同而经线长度不同;所有经线长度相等。,面积变形:地球仪上,同一纬度带内,经差相同的网格面积相等;同一经度带内,纬度越高,面积越小。,角度变形:地球仪上,经线与纬线处处呈直角相交。,4、地图投影地图投影的变形这种变形使得地理要素的几何特性,地图投影变形的图解示例,(,摩尔维特投影等积伪圆柱投影,),长度变形,角度变形,地图投影变形的图解示例(摩尔维特投影等积伪圆柱投影)长度,地图投影变形的图解示例,(,UTM,横轴等角割圆柱投影),面积变形和长度变形,地图投影变形的图解示例(UTM横轴等角割圆柱投影)面积变,地图投影的变形示意,4、地图投影,地图投影的变形,地图投影的变形示意4、地图投影地图投影的变形,区分等角、等面积、任意投影,形状不变,面积不变,特定方向,距离不变,区分等角、等面积、任意投影形状不变面积不变特定方向,地图投影,地图投影的分类,圆柱投影,方位投影,圆锥投影,地图投影地图投影的分类圆柱投影,正轴切圆锥投影 正轴割圆锥投影,横轴切圆锥投影 横轴割圆锥投影,横轴切圆柱投影 横方位投影,正轴割圆柱投影 斜轴切圆柱投影,斜轴切圆锥投影 正轴切圆柱投影,正方位投影 斜方位投影,正轴切圆锥投影 正轴割圆锥投影横轴切圆锥投影,中国地图投影系统,采用与,我国基本图系列,一致,的地图投影系统:即大比例尺时的高斯-克吕格投影(横轴等角切椭圆柱投影)和中小比例尺时的Lambert投影(正轴等角割圆锥投影)。我国常用的地图投影的情况为:,1)、我国基本比例尺地形图(1:100万、1:50万、1:25万、1:10万、1:5万、1:2.5万、1:1万、1:5000),除1:100万外均采用,高斯,克吕格投影,为地理基础;,2)、我国1:100万地形图采用了,Lambert投影,,其分幅原则与国际地理学会规定的,全球统一使用的国际百万分之一,地图投影保持一致。,3)、我国大部分,省区图,以及大多数,这一比例尺,的地图也多采用,Lambert,投影和属于同一投影系统的,Albers,投影(正轴等面积割圆锥投影);,4)、Lambert投影中,地球表面上两点间的,最短距离,(即大圆航线)表现为近于直线,这,有利于,地理信息系统中,空间分析量度,的正确实施。,中国地图投影系统 采用与我国基本图系列一致的地图投影系,对于大中比例尺地图,一般来说大多数都采用地形图的数学基础高斯克吕格投影,尤其是当比例尺为国家基本地形图比例尺系列时,可直接判定为高斯克吕格投影。其原因是,这些比例尺和基本地形图比例尺相一致,编图时,选用地形图的数学基础,既免去了重新展绘数学基础的工序,而且能够保持很高的点位精度。,大中比例尺地图,对于大中比例尺地图,一般来说大多数都采用地,高斯,克吕格投影(Gauss-Kruger Projection),高斯克吕格投影是由高斯于19世纪20年代拟定,后经克吕格补充而形成的一种地图投影方式。在英美国家称为横轴墨卡托投影,属于横轴等角切圆柱投影。这种投影是将椭圆柱面套在地球椭球的外面,并与某一子午线相切(此子午线叫中央子午线或中央经线),椭圆柱的中心轴通过地球椭球的中心,然后用等角条件将中央子午线东西两侧各一定经差范围内的地区投影到柱面上,并将此柱面展成平面,即获得高斯投影,高斯克吕格投影(Gauss-Kruger Project,横轴圆柱投影,高斯,克吕格投影,(Gauss-Kruger Projection),横轴圆柱投影 高斯克吕格投影(Gauss-Kruger P,x,y,高斯-克吕格投影原理图,xy高斯-克吕格投影原理图,高斯,克吕格投影(Gauss-Kruger Projection),高斯投影特征:,中央经线和赤道投影为互相垂直的直线,且为投影的对称轴,投影后无角度变形,即保角投影,中央经线无长度变形,同一条经线上,纬度越低,变形越大,赤道处最大,同一条纬线上,离中央经线越远,变形越大;,为了保证地图的精度,采用分带投影方法,即将投影范围的东西界加以限制,使其变形不超过一定的限度,这样把许多带结合起来,可成为整个区域的投影,在,6带,范围内,长度变形线最大不超过0.14%,高斯克吕格投影(Gauss-Kruger Project,3度带和6度带,从0度开始,自西向东每6度分为一个投影带。,从东经1度30分开始,自西向东每3度分为一个投影带。,高斯,克吕格投影,(Gauss-Kruger Projection),-,投影分带,3度带和6度带 高斯克吕格投影(Gauss-Kruger,我国1:1万至1:50万的地形图全部采用高斯克吕格投影。1:2.5万至1:50万的地形图,采用6分带方案,全球共分为60个投影带;我国位于东经72到136间,共含11个投影带;1:1万比例尺图采用3分带方案,全球共120个带。,分割条带号规定:从0,子午线开始分6,经度为一带,东半球东经3,、9,、15,177,分别是1、2、3,30条6,带的中央子午线,然后继续自西向东旋转,每转6,增加带号1。,分割3,带原则上与6,带相同,只是从东经1,30(即1.5,E)起,每隔3,带为1个投影带。,高斯,克吕格投影,(Gauss-Kruger Projection),-,投影分带,我国1:1万至1:50万的地形图全部采用高斯克吕,在高斯克吕格投影上,规定以中央经线为X轴,赤道为Y轴,两轴的交点为坐标原点。,X坐标值在赤道以北为正,以南为负;Y坐标值在中央经线以东为正,以西为负。我国在北半球,X坐标皆为正值。Y坐标在中央经线以西为负值,运用起来很不方便。为了避免Y坐标出现负值,通常将各带的坐标纵轴西移500公里,即将所有Y值都加500公里。,高斯,克吕格投影,(Gauss-Kruger Projection),-,投影分带,在高斯克吕格投影上,规定以中央经线为X轴,赤道,正轴等角割圆锥投影(Lambert投影),这种投影是将一圆锥面套在地球椭球外面,将地球表面上的要素投影到圆锥面上,然后将圆锥面沿某一母线(经线)展开,即获得Lambert投影。,这种投影中,经线为交于一点的直线束,纬线为同心圆圆弧,圆心即直线束的交点,经线呈辐射状,为纵向直线,纬线近似于弧形,与经线正交,适用于1:100万(包括1:100万)以上地形图,正轴等角割圆锥投影(Lambert投影) 这,正轴圆锥投影,Lambert投影(正轴等角割圆锥投影)。,正轴圆锥投影 Lambert投影(正轴等角割圆锥投影)。,正轴割圆锥投影(Lambert投影),误差情况:,圆锥与地球相交处为北纬25,与北纬47,,距离误差随地点纬度不同而不同,在成图范围内北部最大达+4%,南部达3%,中部为-1.8%,面积变形相对误差相比距离相对误差要大一倍。,一幅图可覆盖大片中纬度地区,可整幅覆盖我国境内领土;,地球表面上两点间的最短距离表现为近于直线,这有利于地理信息系统中的空间分析和信息量度的正确实施,正轴割圆锥投影(Lambert投影)误差情况: 一幅图可,小比例尺地图,由于我国位于中纬度地区,中国地图和分省地图经常采用割圆锥投影(,Lambert,或Albers 投影),中国地图的中央经线常位于东经105度(110度),两条标准纬线分别为北纬25度和北纬47度,而各省的参数可根据地理位置和轮廓形状初步加以判定。例如甘肃省的参数为:中央经线为东经101度,两条标准纬线分别为北纬34度和41度。,小比例尺地图 由于我国位于中纬度地区,中国地,GIS概论2地理空间数学基础课件,GIS概论2地理空间数学基础课件,GIS概论2地理空间数学基础课件,GIS概论2地理空间数学基础课件,GIS概论2地理空间数学基础课件,GIS概论2地理空间数学基础课件,GIS概论2地理空间数学基础课件,GIS概论2地理空间数学基础课件,GIS概论2地理空间数学基础课件,GIS概论2地理空间数学基础课件,GIS概论2地理空间数学基础课件,GIS概论2地理空间数学基础课件,GIS概论2地理空间数学基础课件,GIS概论2地理空间数学基础课件,GIS概论2地理空间数学基础课件,地图投影的选择,地图投影将直接影响地图的精度和使用价值。通常地图投影对中小比例尺地图影响很大,对于大比例尺地图,则影响很小。一般国家基本比例尺地形图的地图投影选择是由国家测绘部门制订,不允许随便更改。,地图投影的选择主要考虑以下因素:制图区域的范围、形状和地理位置;地图的用途、出版方式及其他要求等。,地图投影的选择 地图投影将直接影响地图的精度和使用,地理坐标地图(球面坐标表示为平面方式),地理坐标地图(球面坐标表示为平面方式),等积圆柱投影,等积圆柱投影,等积圆柱投影等积圆柱投影,墨卡托投影:正轴等角切圆柱投影,编制世界时区图,横轴墨卡托投影,墨卡托投影:正轴等角切圆柱投影,墨卡托投影:正轴等角切圆柱投影编制世界时区图横轴墨卡托投影墨,横轴墨卡托投影,横轴墨卡托投影,横轴墨卡托投影横轴墨卡托投影,摩尔维特(Mollweide)投影,等面积伪圆柱投影,常用于编制小比例尺世界地图,摩尔维特(Mollweide)投影等面积伪圆柱投影,常用于编,常见的地图投影,世界图,全国图:正轴圆锥投影,地形图:高斯克吕格投影(分带),南北半球(或两极)图:正轴方位投影、,亚洲图:斜轴方位投影,半球图,东西半球图:横轴方位投影,欧洲图:彭纳投影,国内出版:等差分纬线多圆锥投影,国外出版:摩尔威特投影,地球仪:普通多圆锥投影,海 图:墨卡托投影,中国图,大洲图,常见的地图投影世界图全国图:正轴圆锥投影地形图:高斯克吕,GIS与地图投影的关系,地理基础,(地图投影),地图投影对GIS的影响渗透在GIS建设的各个阶段和各个方面。,(据边馥苓,1996),数据获取,(数据源地图的投影),数据标准化预处理,(按某一参照系数字化),数据存储,(统一的坐标基础),数据处理,(投影转换),数据应用,(检索查询、覆盖分析等),数据输出,(具有相应投影的地图),GIS与地图投影的关系 地理基础地图投影对GIS的影响渗透在,大地原点,大地原点,水准原点,水准原点1985国家高程基准,,72.2604米,黄海海面,1952-1979年平均海水面为0米,水准原点水准原点1985国家高程基准,黄海海面,地图投影与坐标系的关系,地图投影,椭球体模型,Back,地图投影与坐标系的关系地图投影椭球体模型Back,国家坐标系和独立坐标系的变换,由于地球半径很大,在较小区域内进行测量工作可将地球椭球面作为平面看待,而不失其严密性。既然把投影基准面作为平面,就可采用平面直角坐标系表示地面点的投影面上的位置。,(,a,)测量平面直角坐标系,(,b,)数学平面直角坐标系,国家坐标系和独立坐标系的变换 由于地球半径很大,在,为不使坐标系出现负值,它通常将某测区的坐标原点设在测区西南角某点,以真北方向或主要建筑物主轴线为纵轴方向,而以垂直于纵坐标轴的直线定为横坐标轴,构成平面直角坐标系;也可假设测区中某点的坐标值,以该点到另一点方位角作为推算其它各点的起算数据,实际上也构成了一个平面直角坐标系。,上述平面直角坐标系的原点和纵轴方向选定了的值常用于小型测区的测量,它不与国家统一坐标系相连,因此称为任意坐标系或独立坐标系。我国大部分城市均采用独立坐标系,如广州市采用珠江高程和平面坐标系等。,国家坐标系和独立坐标系的变换,为不使坐标系出现负值,它通常将某测区的坐标原点设,按高斯投影统一分带(,6,0,带,,3,0,带)建立的直角坐标系,称为国家平面直角坐标系。,在建立数字城市时,往往需要将独立坐标系转换成国家平面直角坐标系。在进行转换时,先将独立坐标系的原点或独立坐标系的某一固定点与国家大地点连测,并按计算出的方位角进行改正,求出该点的国家统一坐标,然后对所有数据进行平移和旋转,以便把按独立坐标系所采集的数据转换到国家平面直角坐标系中。在城市和工程测量中,也可采用,1.5,0,带或任意带的高斯平面坐标系,以提高投影的精度。,按高斯投影统一分带(6 0带,3 0带)建立的直角坐,地理空间数学基础,地球空间参考,地球椭球体、大地水准面、投影,空间数据投影,地图投影的分类、常用地图投影简介,空间坐标转换,空间尺度,比例尺、分辨率,地理格网,地理格网标准、区域划分标准、国家基本比例尺 地形图标准,地理空间数学基础地球空间参考,反解变换法(同一坐标基准下的坐标变换),通过中间过渡的方法,反解出原投影点的地理坐标,代入新投影中求得新投影之坐标。,反解变换法(同一坐标基准下的坐标变换),例1 从Albers投影转换到Transvers投影,例1 从Albers投影转换到Transvers投影,不同坐标基准下的坐标变换,1、地理坐标基准的变换,2、坐标值的变换,不同坐标基准下的坐标变换1、地理坐标基准的变换,地理空间数学基础,地球空间参考,地球椭球体、大地水准面、投影,空间数据投影,地图投影的分类、常用地图投影简介,空间坐标转换,空间尺度,比例尺、分辨率,地理格网,地理格网标准、区域划分标准、国家基本比例尺 地形图标准,地理空间数学基础地球空间参考,地理空间数学基础,地球空间参考,地球椭球体、大地水准面、投影,空间数据投影,地图投影的分类、常用地图投影简介,空间坐标转换,空间尺度,比例尺、分辨率,地理格网,地理格网标准、区域划分标准、国家基本比例尺 地形图标准,地理空间数学基础地球空间参考,四.地形图分幅与编号,1. 1:100万地形图的分幅编号,经差6度、纬差4度;,纬度60度76度为经差12度、纬差4度;,纬度76度88度之间经差24度、纬差4度。,分幅编号,:,从经度180度起自西向东每经差6度为一列,依次用阿拉伯数字1,2,3,,表示相应列号;,从赤道算起,每纬差4度为一行,至南、北纬88度各分为22行,依次用大写字母A、B、C、V表示相应行号。,四.地形图分幅与编号1. 1:100万地形图的分幅编号,由经线和纬线所围成的每个梯形小格为一幅1:100万地形图,编号由其所在的行号和列号组合而成。如北京为J 50。,注意:,(1,)6度投影带的分带号与百万分之一地形图分幅的纵行相差,30,。,(,2,)我国小于百万分之一的地图(含百万分之一)用等角圆锥投影,大于,50,万分之一的地图都用高斯克吕格投影。,由经线和纬线所围成的每个梯形小格为一幅1:100万地形图,编,GIS概论2地理空间数学基础课件,1、1/1000 000地图的分幅编号,1、1/1000 000地图的分幅编号,GIS概论2地理空间数学基础课件,旧地形图编号的一种:,1,100万:,J-50,1,50万: 在1100万上分4幅,J-50-A,1,25万: 在1100万上分16幅,J-50-A-a,1,20万: 在1100万上分36幅,J-50-A-1,1,10万: 在1100万上分144幅,J-50-144,1,5万: 在110万上分4幅,J-50-144-A,1,2.5万: 在110万上分16幅,J-50-144-A-10,1,1万: 在110万上分64幅,J-50-144-A-1,1,5千: 在110万上分256幅,J-50-144-A-200,2. 20世纪90年代以前地形图编号(不规范),旧地形图编号的一种:2. 20世纪90年代以前地形图编号(,2、1/500 000、 1/200 000 和1/100 000地图的分幅编号,1/500 000地图,1/200 000地图,1/100 000,地图,(分幅编号示例),2、1/500 000、 1/200 000 和1/100,GIS概论2地理空间数学基础课件,3. 20世纪90年代后地形图编号,1:50万1:5千地形图的编号均以1:100万地形图编号为基础,组成,:其所在1:100万地形图的图号、比例尺代码和各图幅的行列号(共10位),比例尺,1:500000,1:250000,1:100000,1:50000,1:25000,1:10000,1:5000,代码,B,C,D,E,F,G,H,比例尺代码,行列号用三位阿拉伯数字表示,不足前面补零,行号在前、列号在后,3. 20世纪90年代后地形图编号1:50万1:5千地形,例1. 1:50万地形图编号,J50B001001;,J50B001002;,J50B002001;,J50B002002。,36,38,0,40,117,120,001001,002001,002002,001002,J50,例1. 1:50万地形图编号363804011,例2 1:25万地形图编号,J50C001001;,J50C001002;,J50C001003;,J50C001004;,J50C004001;,J50C004004,38,37,39,115,30,118,30,004002,004001,004003,004004,003001,003002,003004,003003,J50,001001,001003,001004,001002,002001,002002,002004,002003,例2 1:25万地形图编号3837391153,例3 1:10万地形图编号,J50D001001;J50D001002,.,J50D012001;J50D012002 .,J50D012011;J50D012012,36,120,117,114,40,39,38,37,经差6,纬差4,例4 1:5万地形图编号,J50E001001;J50E001002 .,J50E024001;J50E024002 .,J50E024023;J50E024024,例3 1:10万地形图编号36120117114,三.编号应用的公式,1、已知图幅某点的经、纬度或图幅西南图廓点的经纬度,计算其编号,1:100万地形图图幅编号计算公式,a,=,/4,+1,b=,/6,+31,其中:, -表示商取整;,a -,1:100万地形图图幅所在纬度带字符码所对应的数字码;,b -,1:100万地形图图幅所在经度带对应的数字码;,- 图幅内某点的经度或图幅西南图廓点的经度;,- 图幅内某点的纬度或图幅西南图廓点的纬度。,三.编号应用的公式1、已知图幅某点的经、纬度或图幅西南图廓点,例1:某点经度为1143345,纬度为392230,计算其所在的图幅的编号。,a =,392230/4+1=10(字母J),b =,1143345/6+3150,该点在1:100万地形图图号为J-50。,例1:某点经度为1143345,纬度为392230,其它比例尺地形图在1:100万地形图图号后的行列号,c =,4/,-(,/4)/,d,=,(,/6)/,+1,其中:,( )-表示商取余;, -表示商取整;,c -,所求比例尺地形图在1:100万地形图图号后的行号;,d -,所求比例尺地形图在1:100万地形图图号后的列号;,- 图幅内某点的经度或图幅西南图廓点的经度;,- 图幅内某点的纬度或图幅西南图廓点的纬度;,-所求比例尺地形图分幅的经差;,-所求比例尺地形图分幅的纬差。,其它比例尺地形图在1:100万地形图图号后的行列号,1:50万: ,2,c=4,/2,-(39,2230/ 4,)/2,=,22302,001,d=,(114,3345/ 6,)/,+1,=33,45/,+1001,所以图号为J50B001001,例2 经度114,3345,纬度39,2230计算各比例尺地形图的编号。,1:50万: 2 例2 经度1,1:25万:,1,1,30,c=4,/1,-(39,2230/ 4,)/1,=4,22301,001,d=(114,3345/ 6,)/1,30+1,=33,45/1,30+1=001,1:25万地形图的图号为J50C001001,1:25万: 1 130,1:10万:,20 ,30,c =,4,/20-(39,2230/ 4,)/20,=12,223020=002,d =,(114,3345/ 6,)/30+1,=33,45/30+1=002,1:10万地形图的图号为J50D002002,1:10万: 20 30,1:5万:,10 15,c = 4,/10-(39,2230/ 4,)/10,=24,223010=004,d = (114,3345/ 6,)/15+1,=33,45/15+1=003,1:5万地形图的图号为J50E004003,1:5万: 10 15,1:2.5万: 5 730,c = 4,/5-(39,2230/ 4,)/5,=48,22305,=008,d = (114,3345/ 6,)/730+1,=33,45/730+1,=005,1:2.5万地形图的图号为J50F008005,1:2.5万: 5 730,1:1万:,230 345,c = 4,/230-(39,2230/ 4,)/230,= 96,2230230=015,d = (114,3345/ 6,)/345+1,=33,45/345+1=010,1:1万地形图的图号为J50G015010,1:1万: 230 345,1:5千:,115,152.5,c =,4,/115-(39,2230/ 4,)/115,=192,2230115=030,d = (114,3345/ 6,)/152.5+1,=33,45/152.5+1=019,1:5千地形图的图号为J50H030019,1:5千: 115 152.5,2、已知图号计算该图幅西南图廓点 的经、纬度,=,(,b,-31)6+(,d -,1 ),=,(,a-,1)4+(4/,- c,),式中:,图幅西南图廓点的经度;,图幅西南图廓点的纬度。,a 1:100万地形图图幅所在纬度带字符码所对应的数字码;,b 1:100万地形图图幅所在经度带的数字码;,c 该比例尺地形图在1:100万地形图图号后的行号;,d 该比例尺地形图在1:100万地形图图号后的列号;, 所求比例尺地形图分幅的经差;, 所求比例尺地形图分幅的纬差。,2、已知图号计算该图幅西南图廓点 的经、纬度=(,例1:图号为J50B001001,计算其西南图廓点的经、纬度。,a,=10 ;,b,=50 ;,c =,1 ;,d,=1 ;, = 2,;,= 3,=,(50,31)6,+(1,1 )3,= 114,=,(10,1)4,+(4,/ 2,1 )2,= 38,该图幅西南图廓点的经、纬度分别为114,、38,例1:图号为J50B001001,计算其西南图廓点的经、纬度,例2:图号为J50D002002,计算其西南图廓点的经、纬度。,a,=10 ;,b,=50 ;,c =,2 ;,d,= 2 ;, = 20 ;,= 30,=,(50,31)6,+(2,1 )30= 114,30,=,(10,1)4,+(4,/ 20,2)20= 39,20,该图幅西南图廓点的经、纬度分别为114,30、39,20,例2:图号为J50D002002,计算其西南图廓点的经、纬度,A:西北角图幅的行、列号按下式计算:,c,大,=,小,/,大,(,c,小,-,1)+ 1,d,大,=,小,/,大,(,d,小,-,1)+ 1,3、在同一幅1:100万地形图图幅内,38,37,39,115,30,118,30,004002,004001,004003,004004,003001,003002,003004,003003,J50,001001,001003,001004,001002,002001,002002,002004,002003,(1) 小比例尺地形图的行、列号计算大比例尺地形图的行、列号,A:西北角图幅的行、列号按下式计算: 3、在同一幅1:100,B:东南角图幅的行、列号按下式计算:,c,大,= c,小,小,/,大,d,大,= d,小,小,/,大,c,大,- 较大比例尺地形图在1:100万地形图图号后的行号;,d,大,- 较大比例尺地形图在1:100万地形图图号后的列号;,c,小,- 较小比例尺地形图在1:100万地形图图号后的行号;,d,小,- 较小比例尺地形图在1:100万地形图图号后的列号;,大,- 较大比例尺地形图分幅的纬差;,小,- 较小比例尺地形图分幅的纬差。,B:东南角图幅的行、列号按下式计算:,例 1:10万地形图的行、列号为004001,求所含1:25000地形图的行、列号。,c,小,=,4 ;,d,小,=1 ;,小 = 20;,大,= 5,西北角图幅的行、列号:,c大 =,20/5(4,-,1 )+ 1 = 013,d大 =,20/5(1,-,1 )+ 1 = 001,东南角图幅的行、列号:,c大 =,420/5 = 016,d大 =,1,20/5= 004,所含1:2.5万地形图的行、列号为:,013001 013002 013003 013004,014001 014002 014003 014004,015001 015002 015003 015004,016001 016002 016003 016004,例 1:10万地形图的行、列号为004001,求所含1,(2) 较大比例尺地形图的行、列号计算较小比例尺地形图的行、列号,c,小,=,c大,/(小/,大,)+ 1,1),有余+1,d,小,=,d大,/(小/,大,)+ 1,1),无余+0,式中:, -表示商取整;,c,小,-,较小比例尺地形图在1:1000000地形图图号后的行号;,d,小,-,较小比例尺地形图在1:1000000地形图图号后的列号;,c大 -,较大比例尺地形图在1:1000000地形图图号后的行号;,d大 -,较大比例尺地形图在1:1000000地形图图号后的列号;,大,- 较大比例尺地形图分幅的纬差;,小- 较小比例尺地形图分幅的纬差。,38,37,39,115,30,118,30,004002,004001,004003,004004,003001,003002,003004,003003,J50,001001,001003,001004,001002,002001,002002,002004,002003,(2) 较大比例尺地形图的行、列号计算较小比例尺地形图的行,例:1:2.5万地形图的行、列号分别为016004和013003,计算它隶属于1:100000地形图的行、列号,c,大,=16;d,大,=4;,小,=20;,大,=5,c,小,=16 /(20/5)=004,d,小,=4 /(20/5)=001,c,大,=13 ;d,大,=3 ;,小,=20;,大,=5,c,小,=13 /(20/5)+1=004,d,小,=3 /(20/5)+1 =001,所以行列号为016004和013003的1:2.5万地形图隶属于1:100000地形图的行,列号均为004001。,例:1:2.5万地形图的行、列号分别为016004和0130,5. 坐标网,(1) 经纬线网,经线和纬线构成的坐标网,又称地理坐标网。,作用:控制和确定地表各点和整体地,形的实地位置。,用于分析和计算投影变形;,确定比例尺、量算不可缺少。,5. 坐标网(1) 经纬线网,(2)方里网,由平行于投影坐标轴的两组平行线所构成的方格网。又称直角坐标网。,作用:大比例尺地图上,方便展绘点位和,读图上地理坐标,表示:1:1万1:25万地图上绘出方里网,比例尺,密度,1:1万,1:2.5万,1:5万,1:10万,1:25万,图上距离,(cm),10,4,2,2,4,实地距离,(km),1,1,1,2,10,(2)方里网 比例尺1:1万1:2.5万1:,方里网建立:,在高斯克吕格投影上,规定以中央经线为X轴,赤道为Y轴,两轴的交点为坐标原点。X值在赤道以北为正,以南为负,Y坐标值在中经以东为正,以西为负,我国的X值均为正,但Y值在中经经西为负,运用起来很不方便,故将各带的坐标纵轴西移500km,并冠以带号,称通用坐标,。,方里网建立: 在高斯克吕格投影上,规定以中央经线为X轴,GIS概论2地理空间数学基础课件,