地图分幅编号C 语言编程,地图分幅编号的计算

在之前地图的分幅与编号中我介绍了关于地图的分幅编号的基本知识，这次将用一些实例来讲关于分幅编号的计算。

基本计算有下面几个： 根据某地的经纬度或图幅西南图廓点的经纬度计算出它所在的编号(各种比例尺)

根据图幅编号推算其图幅西南图廓点的经、纬度

不同比例尺地形图编号的行列关系换算

根据某地的经纬度或图幅西南图廓点的经纬度计算出它所在的图幅编号

计算1：100万图幅的编号

\[a=[\phi/4^{0}]+1\]

\[b=[\lambda/6^{0}]+31\] a,b:1:100万图幅行号及列号所对应的数字码

[ ]:表示分数值取整 , 某地的经纬度

计算其他比例尺编号

\[c=4^{0}/\Delta\phi-[(\phi/4^{0})/\Delta\phi]\]

\[d=[(\lambda/6^{0})/\Delta\lambda]+1\]

其中： (  ) 表示商取余;

[ ]表示分数值取整;

c,d:所求比例尺地图的行列号;

\(\lambda \phi\):某地的经纬度;

\(\Delta\phi \Delta\lambda\):所求比例尺地形图分幅的纬差、经差。

已知图幅编号求图幅西南角图廓点的经纬度

求1：100万图

\[\lambda=(b-31)*6^{0}\]

\[\phi=(a-1)*4^{0}\]

a,b:1：100万图的行列数字码

求1：50万--1：5000图

\[\lambda=(b-31)*6^{0}+(d-1)*\Delta\lambda\]

\[\phi=(a-1)*4^{0}+(4^{0}/\Delta\phi-c)*\Delta\phi\] \(\lambda \phi\)：图幅西南图廓点的经纬度

a,b:1：100万的行列数字码

c,d:该比例尺地形图在1:100万地形图编号后的行、列号

\(\Delta\lambda,\Delta\phi\)：该比例尺地形图分幅的经差、纬差

不同比例尺地形图编号的行列关系换算 由较小比例尺地形图编号计算所含各种较大比例尺地形图的编号

例题：

由较大比例尺地形图计算包含该图的较小比例尺地形图编号中的行、列代码