门海元论文.doc

门海元论文
数式法进行地形图分幅和定位的方法
门海元
(山西省林业调查规划院)
摘要按照地形图分幅编号的方法,依据地形图分幅的数学模型和存在的线性关系,对国家基本比例尺地形图分幅编号的数式法进行优化研究, 便于计算机操作以提高其科学性、准确性。适宜于大众化推广使用。关键词数式法地形图分幅与定位
我方公里,这样辽阔的范围包含数以万计的地形图图幅,且种类繁多,比例尺各异。因此,地形图分幅的速度和准确与否直接影响到地形图的发放、保管和使用。
掌握地形图分幅编号,是每个地形图用户必备的一种技能,,以行列式、自然序数式进行编号和分幅的。通常地形图的分幅方法有二种,图解法和数式法,比较起来各有利弊,前者要依靠接图表,地域性很强,精确的分幅尚有困难,也不适宜野外使用。后者数式法计算繁琐,极易出错,但科学含量较高,它用严密的数学模式表达,便于电子计算机处理和实现地图现代化的管理。并可按不同的需要,分幅出任意比例尺地形图,方便实用、操作简单。借助数学模型还可在图上准确定位。

1:100万分幅公式:
横列数=INT(W÷4)+1
纵列数=30+INT(J÷6)+1
1:10万分幅公式
图幅号= (W左上- W) ÷20′×12+(J- J左上) ÷30′+1
1:5万分幅公式
图幅号= (W左上- W) ÷10′×2+(J- J左上) ÷15′+1
图上坐标:X= (-×W)÷15'×(J- (INT(FracJ÷15')×15'+INT J))
Y=(+×W)÷10'×(W-( INT(FracW÷10') ×10'+INT W))
1:
图幅号= (W左上- W) ÷5′×2+(J- J左上) ÷7′30′+1
1:1万分幅公式
图幅号= (W左上- W) ÷2′30"×2+(J- J左上) ÷3′45"+1
图上坐标:X=:(-×W)÷3'45"×(J-( INT(FracJ÷3'45")×3'45"+INT J))
Y=(+×W) ÷2'30"×(W- (INT(FracW÷2'30") ×M+INT W))
式中:W-纬度;J-经度;INT-取整数;Frac-取小数;X-横坐标;Y-纵坐标。为了操作简便,坐标原点设在地形图左下角。

已知东经111゜45'40",北纬38゜37'50"处发生森林火灾,试计算1:5万和1:1万地形图分幅号,并在图上确定准确位置。
利用Excel程序计算
打开Excel程序,在菜单中选择“工具→宏Visual Basic编辑器”打开Microsoft Visual Basic编辑窗口。在该窗口菜单中选择:“插入→模块”,此时会自动弹出一个标题为“Bookl→模块1(代码) ”窗口,这个窗口就是我们编写代码的场所。同时在左边的工程资源