使用FME根据图幅号生成标准图.docVIP

使用FME根据图幅号生成标准图框 在进行基础地理数据处理的过程中，经常会遇到生成标准图框，或检查图框正确性的的问题，相信GIS或测绘人也都遇到过这个问题，下面就来分享下我使用FME来根据图幅号生成1:10000基础地理地形数据标准图幅图框的经验。 1.首先简单复习下标准图幅的分幅方式 新的地形图分幅标准是指1991年国家测绘局制订并颁布实施的“国家基本比例尺地形图分幅和编号国家标准” 1.1地形图分幅 各种比例尺地形图均以1：100万地形图为基础图，沿用原分幅各种比例尺地形图的经纬差（如图1），全部由1：100万地形图按相应比例尺地形图的经纬差逐次加密划分图幅，以横为行，纵为列。 图1 1.2地形图编号 1：100万地形图新的编号方法，除行号与列号改为连写外，没有任何变化，如北京所在的1：100万地形图的图号由J-50改写为J50。（2）1：50万至1：5000地形图的编号，均以1：100万地形图编号为基础，采用行列式编号法，将1：100万地形图按所含各种比例尺地形图的经纬差划分成相应的行和列，横行自上而下，纵列从左到右，按顺序均用阿拉伯数字编号，皆用3位数字表示，凡不足3位数的，则在其前补0。 各大中比例尺地形图的图号均由五个元素10位码构成。从左向右，第一元素1位码，为1：100万图幅行号字符码；第二元素2位码，为1：100万图幅列号数字码；第三元素1位码，为编号地形图相应比例尺的字符代码；第四元素3位码，为编号地形图图幅行号数字码；第五元素3位码，为编号地形图图幅列号数字码；各元素均连写（如图2）。比例尺代码见图3； 图2 图3 了解到图幅分幅编号的方式以后，可以得出图幅号与标准图幅角点坐标的关系。如图4 图4 2.使用FME生成标准图框 了解完图幅名与标准图框的关系之后，就可以使用FME来根据图幅名生成标准图框了，该模板使用来生成1:10000标准图幅。 图5 模板介绍：前面四个SubstringExtractor分别用来提取H，h，L，l即图4中对应的四个参数，AttributeValueMapper用来将英文字母映射成为数字： 图6 两个ExpressionEvaluator分别用来计算x，y图4中的角点坐标公式如下 x=(@double(@Value(L))-31)*6+(@double(@Value(l))-1)*3.75/60 y=(@double(@Value(H))-1)*4+( HYPERLINK mailto:96-@double(@Value(h)))*2.5/60 96-@double(@Value(h)))*2.5/60 下面的2dpointreplacer是生成左下角点，2dpointadder是根据经差和纬差得到其余三个角点坐标，注意度秒分的换算 图7 下面一个AttributeRangeMapper和AttributeCreator是根据图幅号计算中央经线以便于用来投影。 ?图8 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 图9 计算出标准图幅的四个角点经纬度坐标之后，就可以为它指定坐标系和进行重投影了。