GIS概论2地理空间数学基础课件

1、GIS概论2地理空间数学基础GIS概论2地理空间数学基础地理空间的数学建模 为了深入研究地理空间，需要建立地球表面的几何模型，这是进行大地测量的前提。根据大地测量学的成果，地球表面模型可以分为四类： 最自然的面：包括海洋底部、高山、高原等在内的固体地球表面。太复杂，难以建模，各种量算也非常困难。 相对抽象的面：也称为大地水准面，是静止海平面的延伸。以它为基准，可以用水准仪测量地球自然表面上任意点的高程。海平面的起伏将导致测量的不确定。 模型：以大地水准面为基准建立的地球椭球体模型。 其他数学模型：为了解决特定的大地测量问题而提出的。如类地形面、准大地水准面、静态水平衡椭球体等。地理空间的数学建

2、模 为了深入研究地理空间，需要建立地球地球模型：三级近似地球自然表面极不规则，无法用数学表面进行描述水准面所包围的球体大地水准面所包围的球体旋转椭球体不规则性、动态性、不唯一性不规则性、相对唯一性标准数学曲面 1952：海福特椭球 1953：克拉索夫斯基椭球 1978：1975年国际椭球地球模型：三级近似地球自然表面极不规则，无法用数学表面进行描地球椭球体模型：以大地水准面为基准建立的。地球的形状接近于椭圆绕其短轴形成的椭球体，通过扁率表示椭球体的扁平程度。大地水准面与具有微小扁率的旋转椭球面非常接近，可用旋转椭球体代替大地球体。 2、地球空间模型描述 地理空间数学建模 2、地球空间模型描述

3、地理空间数学建模地球椭球体模型bca三轴椭球体模型双轴椭球体模型（旋转椭球体）其他椭球体模型：根据a、b、c的不同x2a2y2b2z2c2+= 1x2a2y2b2z2a2+= 1如：克拉索夫斯基椭球体 2、地球空间模型描述 地理空间数学建模地球椭球体模型bca三轴椭球体模型双轴椭球体模型（旋转椭球体椭球体参数长半径 a（赤道半径）短半径 b（极半径）扁 率 =(a-b)/a第一偏心率 e2=(a2- b2)/ a2第二偏心率 e2=(a2- b2)/ b2我国使用的椭球克拉索夫斯基椭球体IAG75椭球体WGS84椭球体 我国的大地坐标系和高程系1954年北京坐标系1980年国家大地坐标系西安原

4、点GPS测量数据1956年黄海高程系1985年国家高程基准椭球体参数长半径 a（赤道半径）我国使用的椭球我中国曾经用过的椭球体我国1952年以前采用Hayford椭球体；1953年开始采用克拉索夫斯基椭球体；上世纪70年代末建立新的80坐标系时，采用IUGG（国际大地测量与地球物理联合会）椭球体；1984年定义的世界大地坐标系（WGS84）使用的椭球体长、短半径则分别为6378.137和6356.7523，扁率为1：298.26。中国曾经用过的椭球体我国1952年以前采用Hayford椭直接建立在球体上的地理坐标，用经度和纬度表达地理对象位置建立在平面上的直角坐标系统，用（x，y）表达地理对象

5、位置投影定位坐标系：平面系统直接建立在球体上的地理坐标，用经度和纬度表达地理对象位置建立3.2.2 大地坐标系（1）54年北京坐标系 在东北黑龙江边境上同苏联大地网联测，通过大地坐标计算，推算出北京点的坐标，北京坐标系是苏联42年坐标系的延伸，其原点在苏联普尔科沃。（2）80年西安坐标系 78年4月召开“全国天文大地网平差会议”建立80年西安坐标系，其原点在西安西北的永乐镇，简称西安原点。椭球体参数为75年国际大地测量与地球物理联合会第16界大会的推荐值。（3）新54年北京坐标系 将全国大地网整体平差的结果整体换算到克拉索夫斯基椭球体上，形成一个新的坐标系，称为新54年北京坐标系，它与80年国

6、家大地坐标系的轴定向基准相同，网的点位精度相同。（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)

7、2、平面定位系统现实世界是三维大地原点水准原点？为什么要把球面的转换为平面的？怎么把球面的转换为平面的？2、地理要素的定位系统为了确定地理要素的位置，需要确定其定位为什么要进行地图投影？将地球椭球面上的点映射到平面上的方法，称为地图投影。1、地理坐标为球面坐标，不方便进行距离、方位、面积等参数的量算。2、地球椭球体为不可展曲面。3、地图为平面，符合视觉心理，并易于进行距离、方位、面积等量算和各种空间分析。为什么要进行地图投影？将地球椭球面上的点映射到平面上的方法地理空间数学基础地球空间参考 地球椭球体、大地水准面、投影 空间数据投影 地图投影的分类、常用地图投影简介 空间坐标转换 空间尺度 比

8、例尺、分辨率 地理格网 地理格网标准、区域划分标准、国家基本比例尺 地形图标准地理空间数学基础地球空间参考3、地图投影、什么是地图投影？、坐标系、地图投影的分类、如何区分等角、等面积、任意投影？、中国地图投影的采用方式、常见的地图投影3、地图投影、什么是地图投影？什么是地图投影？ 简单地讲：地图投影的实质是将地球椭球面上的经纬网按照一定的数学法则转移到平面上。 具体来说：由于球面上一点的位置是用地理坐标表示，而平面上是用直角坐标或者极坐标表示，所以要想将地球表面上的点转移到平面上，必须采用一定的数学方法来确定地理坐标与平面直角坐标或极坐标之间的关系。 这种在球面和平面之间建立点与点之间函数关系

9、的数学方法，称为地图投影。什么是地图投影？ 简单地讲：地图投影的实质是将地球椭球面坐标系选定了一个一定大小的椭球体，并确定了它与大地水准面的相关位置，就确定了一个坐标系。坐标系选定了一个一定大小的椭球体，并确定了它与大地水准面的相地图投影的目的地图投影的目的： 是通过将不可展的球面投影到一个可展曲面上，然后将该曲面展开成为一个平面，来保证空间信息在地域上的连续性、完整性和可测度性。地图投影的目的地图投影的目的：地图投影的分类 根据美国著名地图投影学家J.P.Snyder统计，全世界地图投影种类现有256种，依据不同的目的和要求，可以采用不同的分类指标对如此繁多的地图投影进行分类。 分类1：基于

10、投影面与球面相关位置的分类； 分类2：基于投影方法的分类； 分类3：基于投影方程的分类； 分类4：基于投影变形的分类。地图投影的分类 根据美国著名地图投影学家J.P.Sny基于投影面与球面的分类基于投影面与球面的分类4、地图投影地图投影的变形 用地图投影的方法将球面展开为平面，虽然可以保持地域上的联系和完整性，但它们与球面上的经纬度网线形状并不一致。即投影后，地图上的经纬度网线发生了变形，同样根据地理坐标展绘在地图上的各种要素，也必然随着变形。4、地图投影地图投影的变形 用地图投影的方法将球面4、地图投影地图投影的变形这种变形使得地理要素的几何特性受到破坏： 长度变形：地球仪上，纬线长度不等；

11、同一纬线上，经差相同，纬线长度相同；同一经线上，纬差相同而经线长度不同；所有经线长度相等。 面积变形：地球仪上，同一纬度带内，经差相同的网格面积相等；同一经度带内，纬度越高，面积越小。 角度变形：地球仪上，经线与纬线处处呈直角相交。4、地图投影地图投影的变形这种变形使得地理要素的几何特性地图投影变形的图解示例（摩尔维特投影等积伪圆柱投影）长度变形角度变形地图投影变形的图解示例（摩尔维特投影等积伪圆柱投影）长度地图投影变形的图解示例（UTM横轴等角割圆柱投影）面积变形和长度变形地图投影变形的图解示例（UTM横轴等角割圆柱投影）面积变地图投影的变形示意4、地图投影地图投影的变形地图投影的变形示意4

12、、地图投影地图投影的变形区分等角、等面积、任意投影形状不变面积不变特定方向距离不变区分等角、等面积、任意投影形状不变面积不变特定方向地图投影地图投影的分类圆柱投影方位投影圆锥投影地图投影地图投影的分类圆柱投影正轴切圆锥投影 正轴割圆锥投影横轴切圆锥投影 横轴割圆锥投影横轴切圆柱投影 横方位投影正轴割圆柱投影 斜轴切圆柱投影斜轴切圆锥投影 正轴切圆柱投影正方位投影 斜方位投影正轴切圆锥投影 正轴割圆锥投影横轴切圆锥投影中国地图投影系统 采用与我国基本图系列一致的地图投影系统：即大比例尺时的高斯-克吕格投影（横轴等角切椭圆柱投影）和中小比例尺时的Lambert投影（正轴等角割圆锥投影）。我国常用的

13、地图投影的情况为： 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投影中，地球表面上两点间的最短距离(即大圆航线)表现为近于直线，这有利于地理信息系统中空间分析量度的正确实施。中国地图

14、投影系统 采用与我国基本图系列一致的地图投影系 对于大中比例尺地图，一般来说大多数都采用地形图的数学基础高斯克吕格投影，尤其是当比例尺为国家基本地形图比例尺系列时，可直接判定为高斯克吕格投影。其原因是，这些比例尺和基本地形图比例尺相一致，编图时，选用地形图的数学基础，既免去了重新展绘数学基础的工序，而且能够保持很高的点位精度。大中比例尺地图 对于大中比例尺地图，一般来说大多数都采用地 高斯克吕格投影（Gauss-Kruger Projection） 高斯克吕格投影是由高斯于19世纪20年代拟定，后经克吕格补充而形成的一种地图投影方式。在英美国家称为横轴墨卡托投影 属于横轴等角切圆柱投影。这种投

15、影是将椭圆柱面套在地球椭球的外面，并与某一子午线相切（此子午线叫中央子午线或中央经线），椭圆柱的中心轴通过地球椭球的中心，然后用等角条件将中央子午线东西两侧各一定经差范围内的地区投影到柱面上，并将此柱面展成平面，即获得高斯投影 高斯克吕格投影（Gauss-Kruger Project横轴圆柱投影 高斯克吕格投影（Gauss-Kruger Projection）横轴圆柱投影 高斯克吕格投影（Gauss-Kruger Pxy高斯-克吕格投影原理图xy高斯-克吕格投影原理图 高斯克吕格投影（Gauss-Kruger Projection）高斯投影特征： 中央经线和赤道投影为互相垂直的直线，且为投影的

16、对称轴 投影后无角度变形，即保角投影 中央经线无长度变形 同一条经线上，纬度越低，变形越大，赤道处最大同一条纬线上，离中央经线越远，变形越大；为了保证地图的精度，采用分带投影方法，即将投影范围的东西界加以限制，使其变形不超过一定的限度，这样把许多带结合起来，可成为整个区域的投影 在6带范围内，长度变形线最大不超过0.14% 高斯克吕格投影（Gauss-Kruger Project3度带和6度带从0度开始，自西向东每6度分为一个投影带。从东经1度30分开始，自西向东每3度分为一个投影带。 高斯克吕格投影（Gauss-Kruger Projection） -投影分带3度带和6度带 高斯克吕格投影（

17、Gauss-Kruger 我国1:1万至1:50万的地形图全部采用高斯克吕格投影。1:2.5万至1:50万的地形图，采用6分带方案，全球共分为60个投影带；我国位于东经72到136间，共含11个投影带；1:1万比例尺图采用3分带方案，全球共120个带。 分割条带号规定：从0子午线开始分6经度为一带，东半球东经3、9、15177分别是1、2、330条6带的中央子午线，然后继续自西向东旋转，每转6增加带号1。 分割3带原则上与6带相同，只是从东经130（即1.5E）起，每隔3带为1个投影带。 高斯克吕格投影（Gauss-Kruger Projection） -投影分带 我国1:1万至1:50万的地

18、形图全部采用高斯克吕 在高斯克吕格投影上，规定以中央经线为X轴，赤道为Y轴，两轴的交点为坐标原点。 X坐标值在赤道以北为正，以南为负；Y坐标值在中央经线以东为正，以西为负。我国在北半球，X坐标皆为正值。Y坐标在中央经线以西为负值，运用起来很不方便。为了避免Y坐标出现负值，通常将各带的坐标纵轴西移500公里，即将所有Y值都加500公里。 高斯克吕格投影（Gauss-Kruger Projection） -投影分带 在高斯克吕格投影上，规定以中央经线为X轴，赤道正轴等角割圆锥投影（Lambert投影） 这种投影是将一圆锥面套在地球椭球外面，将地球表面上的要素投影到圆锥面上，然后将圆锥面沿某一母线（

19、经线）展开，即获得Lambert投影。这种投影中，经线为交于一点的直线束，纬线为同心圆圆弧，圆心即直线束的交点经线呈辐射状，为纵向直线，纬线近似于弧形，与经线正交适用于1：100万（包括1：100万）以上地形图正轴等角割圆锥投影（Lambert投影） 这正轴圆锥投影 Lambert投影（正轴等角割圆锥投影）。正轴圆锥投影 Lambert投影（正轴等角割圆锥投影）。正轴割圆锥投影（Lambert投影）误差情况： 圆锥与地球相交处为北纬25与北纬47 ，距离误差随地点纬度不同而不同，在成图范围内北部最大达+4%，南部达3%，中部为-1.8%，面积变形相对误差相比距离相对误差要大一倍。 一幅图可覆盖

20、大片中纬度地区，可整幅覆盖我国境内领土； 地球表面上两点间的最短距离表现为近于直线，这有利于地理信息系统中的空间分析和信息量度的正确实施正轴割圆锥投影（Lambert投影）误差情况： 一幅图可 小比例尺地图 由于我国位于中纬度地区，中国地图和分省地图经常采用割圆锥投影（ Lambert 或Albers 投影），中国地图的中央经线常位于东经105度(110度)，两条标准纬线分别为北纬25度和北纬47度，而各省的参数可根据地理位置和轮廓形状初步加以判定。例如甘肃省的参数为：中央经线为东经101度，两条标准纬线分别为北纬34度和41度。 小比例尺地图 由于我国位于中纬度地区，中国地GIS概论2地理空

21、间数学基础课件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地理空间数学基础课件地图投影的选择 地图投影将直接影响地图的精度和使用价值。通常地图投影对中小比例尺地图影响很大，对于大比例尺地图，则影响很小。一般国家基本比

22、例尺地形图的地图投影选择是由国家测绘部门制订，不允许随便更改。 地图投影的选择主要考虑以下因素：制图区域的范围、形状和地理位置；地图的用途、出版方式及其他要求等。地图投影的选择 地图投影将直接影响地图的精度和使用地理坐标地图（球面坐标表示为平面方式）地理坐标地图（球面坐标表示为平面方式）等积圆柱投影等积圆柱投影等积圆柱投影等积圆柱投影墨卡托投影：正轴等角切圆柱投影编制世界时区图横轴墨卡托投影墨卡托投影：正轴等角切圆柱投影墨卡托投影：正轴等角切圆柱投影编制世界时区图横轴墨卡托投影墨横轴墨卡托投影横轴墨卡托投影横轴墨卡托投影横轴墨卡托投影摩尔维特（Mollweide）投影等面积伪圆柱投影，常用于编

23、制小比例尺世界地图摩尔维特（Mollweide）投影等面积伪圆柱投影，常用于编常见的地图投影世界图全国图：正轴圆锥投影地形图：高斯克吕格投影（分带）南北半球（或两极）图：正轴方位投影、亚洲图：斜轴方位投影半球图东西半球图：横轴方位投影欧洲图：彭纳投影国内出版：等差分纬线多圆锥投影国外出版：摩尔威特投影地球仪：普通多圆锥投影海 图：墨卡托投影中国图大洲图常见的地图投影世界图全国图：正轴圆锥投影地形图：高斯克吕GIS与地图投影的关系 地理基础(地图投影)地图投影对GIS的影响渗透在GIS建设的各个阶段和各个方面。（据边馥苓，1996）数据获取（数据源地图的投影）数据标准化预处理（按某一参照系数字化

24、）数据存储（统一的坐标基础）数据处理（投影转换）数据应用（检索查询、覆盖分析等）数据输出（具有相应投影的地图）GIS与地图投影的关系 地理基础地图投影对GIS的影响渗透在大地原点大地原点水准原点水准原点1985国家高程基准，72.2604米黄海海面1952-1979年平均海水面为0米水准原点水准原点1985国家高程基准，黄海海面地图投影与坐标系的关系地图投影椭球体模型Back地图投影与坐标系的关系地图投影椭球体模型Back国家坐标系和独立坐标系的变换 由于地球半径很大，在较小区域内进行测量工作可将地球椭球面作为平面看待，而不失其严密性。既然把投影基准面作为平面，就可采用平面直角坐标系表示地面点

25、的投影面上的位置。（a）测量平面直角坐标系（b）数学平面直角坐标系国家坐标系和独立坐标系的变换 由于地球半径很大，在 为不使坐标系出现负值，它通常将某测区的坐标原点设在测区西南角某点，以真北方向或主要建筑物主轴线为纵轴方向，而以垂直于纵坐标轴的直线定为横坐标轴，构成平面直角坐标系；也可假设测区中某点的坐标值，以该点到另一点方位角作为推算其它各点的起算数据，实际上也构成了一个平面直角坐标系。 上述平面直角坐标系的原点和纵轴方向选定了的值常用于小型测区的测量，它不与国家统一坐标系相连，因此称为任意坐标系或独立坐标系。我国大部分城市均采用独立坐标系，如广州市采用珠江高程和平面坐标系等。国家坐标系和独

26、立坐标系的变换 为不使坐标系出现负值，它通常将某测区的坐标原点设 按高斯投影统一分带（6 0带，3 0带）建立的直角坐标系，称为国家平面直角坐标系。 在建立数字城市时，往往需要将独立坐标系转换成国家平面直角坐标系。在进行转换时，先将独立坐标系的原点或独立坐标系的某一固定点与国家大地点连测，并按计算出的方位角进行改正，求出该点的国家统一坐标，然后对所有数据进行平移和旋转，以便把按独立坐标系所采集的数据转换到国家平面直角坐标系中。在城市和工程测量中，也可采用1.5 0带或任意带的高斯平面坐标系，以提高投影的精度。 按高斯投影统一分带（6 0带，3 0带）建立的直角坐地理空间数学基础地球空间参考 地

27、球椭球体、大地水准面、投影 空间数据投影 地图投影的分类、常用地图投影简介 空间坐标转换 空间尺度 比例尺、分辨率 地理格网 地理格网标准、区域划分标准、国家基本比例尺 地形图标准地理空间数学基础地球空间参考反解变换法（同一坐标基准下的坐标变换） 通过中间过渡的方法，反解出原投影点的地理坐标，代入新投影中求得新投影之坐标。反解变换法（同一坐标基准下的坐标变换） 例1 从Albers投影转换到Transvers投影例1 从Albers投影转换到Transvers投影不同坐标基准下的坐标变换1、地理坐标基准的变换2、坐标值的变换不同坐标基准下的坐标变换1、地理坐标基准的变换地理空间数学基础地球空间

28、参考 地球椭球体、大地水准面、投影 空间数据投影 地图投影的分类、常用地图投影简介 空间坐标转换 空间尺度 比例尺、分辨率 地理格网 地理格网标准、区域划分标准、国家基本比例尺 地形图标准地理空间数学基础地球空间参考地理空间数学基础地球空间参考 地球椭球体、大地水准面、投影 空间数据投影 地图投影的分类、常用地图投影简介 空间坐标转换 空间尺度 比例尺、分辨率 地理格网 地理格网标准、区域划分标准、国家基本比例尺 地形图标准地理空间数学基础地球空间参考四.地形图分幅与编号1. 1:100万地形图的分幅编号经差6度、纬差4度；纬度60度76度为经差12度、纬差4度；纬度76度88度之间经差24度

29、、纬差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

30、万地形图，编GIS概论2地理空间数学基础课件1、1/1000 000地图的分幅编号1、1/1000 000地图的分幅编号GIS概论2地理空间数学基础课件旧地形图编号的一种：1100万： J-50150万： 在1100万上分4幅 J-50-A 125万： 在1100万上分16幅 J-50-A-a120万： 在1100万上分36幅 J-50-A-1 110万： 在1100万上分144幅 J-50-14415万： 在110万上分4幅 J-50-144-A12.5万： 在110万上分16幅 J-50-144-A-1011万： 在110万上分64幅 J-50-144-A-115千： 在110万上分256

31、幅 J-50-144-A-2002. 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：5000001：2500

32、001：1000001：500001：250001：100001：5000代码BCDEFGH比例尺代码行列号用三位阿拉伯数字表示，不足前面补零，行号在前、列号在后 3. 20世纪90年代后地形图编号1:50万1:5千地形例1. 1：50万地形图编号 J50B001001； J50B001002； J50B002001； J50B002002。3638040117120001001002001002002001002J50例1. 1：50万地形图编号363804011例2 1:25万地形图编号 J50C001001； J50C001002； J50C001003； J50C001004； J50

33、C004001； J50C0040043837391153011830004002004001004003004004003001003002003004003003J50001001001003001004001002002001002002002004002003例2 1:25万地形图编号3837391153例3 1：10万地形图编号J50D001001;J50D001002.J50D012001;J50D012002 .J50D012011;J50D0120123612011711440393837经差6纬差4例4 1：5万地形图编号J50E001001；J50E001002 .J50E

34、024001；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，纬度

35、为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万地形图图号后的列号； - 图幅内某点的经度或图幅西南图廓点的经度； - 图幅内某点的纬度或图幅西南图廓点的纬度；-所求比例尺地形图分幅的经差；-所求

36、比例尺地形图分幅的纬差。其它比例尺地形图在1:100万地形图图号后的行列号 1:50万： 2 c=4/2-（392230/ 4）/2 =22302001 d= （1143345/ 6）/+1 =3345/+1001所以图号为J50B001001例2 经度1143345，纬度392230计算各比例尺地形图的编号。 1:50万： 2 例2 经度11:25万: 1 130 c=4/1-（392230/ 4）/1 =422301001 d=（1143345/ 6）/130+1 =3345/130+1=0011:25万地形图的图号为J50C0010011:25万: 1 1301:10万: 20 30c

37、= 4/20-（392230/ 4）/20=12223020=002d = （1143345/ 6）/30+1 =3345/30+1=0021:10万地形图的图号为J50D0020021:10万: 20 301:5万: 10 15c = 4/10-（392230/ 4）/10=24223010=004d = （1143345/ 6）/15+1 =3345/15+1=0031:5万地形图的图号为J50E0040031:5万: 10 151:2.5万: 5 730 c = 4/5-（392230/ 4）/5 =4822305 =008 d = （1143345/ 6）/730+1 =3345/73

38、0+1 =0051:2.5万地形图的图号为J50F0080051:2.5万: 5 7301:1万: 230 345 c = 4/230-（392230/ 4）/230 = 962230230=015 d = （1143345/ 6）/345+1 =3345/345+1=0101:1万地形图的图号为J50G0150101:1万: 230 3451:5千: 115 152.5 c = 4/115-（392230/ 4）/115 =1922230115=030 d = （1143345/ 6）/152.5+1 =3345/152.5+1=0191:5千地形图的图号为J50H0300191:5千: 1

39、15 152.52、已知图号计算该图幅西南图廓点 的经、纬度=（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 ；

40、c = 1 ；d =1 ； = 2 ；= 3=（5031）6+（1 1 ）3= 114=（101）4+（4/ 2 1 ）2 = 38该图幅西南图廓点的经、纬度分别为114、38例1：图号为J50B001001，计算其西南图廓点的经、纬度例2：图号为J50D002002，计算其西南图廓点的经、纬度。a =10 ；b =50 ； c = 2 ；d = 2 ； = 20 ；= 30=（5031）6+（21 ）30= 11430=（101）4+（4/ 202）20= 3920该图幅西南图廓点的经、纬度分别为11430、3920例2：图号为J50D002002，计算其西南图廓点的经、纬度A：西北角图幅的

41、行、列号按下式计算： c大 =小/大（c小- 1）+ 1 d大 =小/大（d小- 1）+ 13、在同一幅1:100万地形图图幅内3837391153011830004002004001004003004004003001003002003004003003J50001001001003001004001002002001002002002004002003（1） 小比例尺地形图的行、列号计算大比例尺地形图的行、列号A：西北角图幅的行、列号按下式计算： 3、在同一幅1:100B：东南角图幅的行、列号按下式计算： 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小