f图幅分幅编号和图幅范围的快速确定方法.doc

图幅分幅编号和图幅范围的快速确定方法钟业勋, 李占元, 童新华, （，广西测绘局，广西南宁，530023;，广西师范学院资源与环境科学学院，广西，南宁，530001） 摘要：具有地理坐标的测绘资料属于地学空间数据范畴。本文给出了快速确定地学空间数据所属国家基本比例尺地形图的图幅编号及图幅范围的方法，并以实例验证。本方法对测绘资料的分幅存档和分发使用提供了依据和方便。关键词：图家基本比例尺地形图；分幅编号；图幅范围；CASIO fx-4800P计算器；计算程序测绘资料常以图幅为单位分类归档，因此，地面上具有不同地理坐标（经度和纬度）的两个点，是否同属某种比例尺的同一图幅，其图幅编号和图幅范围如能快速确定，显然有利于资料的分幅存档和分发使用。本文是笔者对这一问题的研究和探讨。1 我国国家基本比例尺地形图的分幅编号方法我国国家基本比例尺地形图共有8种：1：100万，1：50万，1：25万，1：10万，1：5万，1：2.5万，1：1万，1：5000。1.1 国家基本比例尺地形图分幅编号的结构国家基本比例尺地形图的分幅编号由5个元素10位码组成。其结构如表1表1 国家基本比例尺地形图分幅编号结构元素位码*表中之*为字母代号，为数字码。分幅以1：100万地图为基础，由1：100万图逐次加密划分而成，横为行，由上而下排序；竖为列，由左向右排序。所有1：50001：50万地形图的图号均由5个元素10位码组成。编码系列统一为一个根部，编码长度相同，方便于计算机处理。各元素的意义如下：1100万图幅行号字符，从赤道到南北纬88，每4为一横行，分别以字母A，B，V，22个字母表示。1100万图幅列号数字码。从180经线自西向东计算，每6为一纵列，分别以1，2，60编号，1位数列（180至西经126）前用“0”补足两位。我国地处东半球，无补“0”问题。比例尺代码，由表3给出。图幅所在的1100万图中的行号数字码，不足三位时在有效数字前加“0”补足。图幅所在的1100万图中的列号数字码，不足三位时在有效数字前加“0”补足。 基金项目：广西自然科学基金（桂科自0448037）作者简介：钟业勋，（1939），男，广西玉林人，教授，主要从事地图学理论研究。1.2 国家基本比例尺地形图的经差，纬差，行列数和图幅关系。表2 国家基本比例尺地形图的经差，纬差，行列数和图幅数关系比例尺1100万150万125万110万15万12.5万11万15 000图幅经纬差经差631303015730345152.5纬差42120105230115行列数量关系行数12412244896192列数12412244896192图幅数量关系14116411443691576144 36 4 12 304 576 144 16 4 19 2162 304 576 64 16 4 136 8649 2162 304 256 64164图1 1：100万1：5000地形图的行、列编号11.3 国家基本比例尺地形图的比例尺代码和比例尺系数。表3 国家基本比例尺地形图的比例尺代码和比例尺系数地图比例尺1100万15 0万12.5万11 0万15万12.5万11万15 0 0 0比例尺代码BCDEFGH基本纬差经差倍数n1929648168421设15 000比例尺地形图的经差0152.5为“基本经差”，其纬差0115为“基本纬差”。表3中之“比例尺系数”，实为图幅的经差（纬差）与基本经差（基本纬差）的比值2。2 计算某点所在某种比例尺地形图的分幅编号及图幅范围的程序2.1 计算程序用CASIO fx-4800P计算器编制程序。以“图幅编号”拼音的首写字母构成文件名DFBH，在此文件名下输入下列程序Deg：Lb11：A，B，C：D = Int（A4）+1Lb12：E = Int（B6）+31Lb13：F = 4 (0115C）Int(AInt(A4）4）0115CLb14：G = Int(B6（E31） (0152.5C）+1Lb15：H= 4DF(0115CLb16：I= 6（E31）+ G(0152.5CLb17：Goto EXIT EXIT程序中各字符代表的意义。A纬度；B经度，经纬度以度单位；C比例尺系数，即表3中之n值；D表1中之元素I，图幅所在的1：100万图的行号所对应的数字，要按1A，2B，22V的对应由数字变换为字符；E元素，图幅所在的1：100万图的列号数字码；F元素，图幅在其所属的1：100万图中的行号数字码，不足3位时在数字前用“0”补足；G元素V，图幅在其所属的1：100万图中的列号数字码，数字行不足3位时在数字前用“0”补足；H图幅南图廓线纬度，加上表2中对应比例尺图幅纬差即得北图廓纬度；I图幅东图廓线经度，减去表2中对应比例尺图幅经差即得西图廓线经度。2.2 计算以A=3803，B=11708，C=16(1:10万比例尺系数)为例，执行程式。步骤键操作显示画面1MODE 5 2 (RUN)ProsramRUNDFBH：CD2EXEA？ （原存数字） 338 。， 03 。， EXE3803D10.000 0004EXEB？ （原存数字） 5117 。， 08 。， EXE11709E50.000 0006EXEC？（原存数字）716 EXEF6.000 0008 EXEG7.000 0009EXEH38.000 00010EXEI117.5000 000显示I值后，再按EXE 又可重复前面的操作，输入新的数据。 3 算例比已知下列各点的地理坐标（经度），（经度）为：3803，11708，3606，11507，3817，11725，3813，11721，3613，11527，3906，11613，3918，11602，3609，11502，3911，11607，3803，11729。求其所在的110万比例尺地形图的图幅编号和图幅范围。解：以纬度A、经度B和110万图比例尺系数C16输入DFBH程序，计算结果如表4表4 计算结果表点号ABCDEFGHI图幅编号及图幅范围3803381738123803117081172511721117291610（J）506738 117.5JD500060073800382011700117303606361336091150711527116131610（J）5012336115.5JD500120033600362011500115303906391839111161311602116071610（J）503539116.5JD50003005390039201160011630从计算结果可知，已知的10个点分属3幅1：10万地形图。4 结语地学空间数据（Spatial Data）是指带有地理坐标即经纬度的数据，包括资源、环境、经济和社会等领域的一切带有地理坐标的数据。据统计资料表明，地学空间数据占数据总量的80%3。地图投影是地图的数学基础，地图投影变换定义为两个二维场间的拓扑变换4。地图符号的位置特征，正是由于拓扑变换原理使其具有反映其表达事物空间位置的功能，从而使地图成为表达地学空间信息的最佳工具。本文给出了具有地理坐标的空间数据快速确定其所在的国家基本比例尺地形图的图幅编号及图幅范围的程序和方法，对测绘资料的分幅存档和分发使用提供了方便，具有理论依据和实用意义。参考文献1 尹贡白.王家耀.田德森等，地图概论M，北京：测绘出版社，1991，52.2 钟业勋.数理地图学M，北京：测绘出版社.2007，6165.3 承继成.林珲.周成虎等，数字地球导论M，北京：科学出版社，2004，19.4 胡鹏，游涟，杨传勇等， 地图代数M，武汉大学版社， 2002，25.Fast Define Mode of SheetDesignation and Sheet LimitZHONG Yexun, LI zhanyuan, TONG Xinhua（.Guangxi Regional Bureau of Surveying and Mapping, Nanning,530023,China;.School of Resource Environment Science, Guangxi TeachersCollege, Nanning,530001,China）Abstract：The document of surveying amd mapping of possessed of geographic coordinate belonging to the category of spatial data .In this paper the authors gave the fast define mode for spatial data behnging to the sheet of national basic scale topegraphic map and its sheet designation and sheet limit, and the mode is proofed by a example .This mode provided the judgement and convenient for the file