高等地图学作业

一、新旧图幅编号分幅

新——基本比例尺地形图分幅

旧——基本比例尺地形图分幅

 ①依据经纬度计算某点所在图幅编号

先计算点P(λ,φ)，所在1：100万图谱的编号

a=[φ/4°]+1；

b=[λ/6°]+31(东经)  \   b=30-[λ/6°](西经)；

再计算所求比例尺地形图的（1：100万图号后面的）图幅编号

c=4°/Δφ-[(φ/4°)/Δφ]；(Δφ为所求比例尺对应的纬差）

d=[(λ/6°)/Δλ]+1；(Δλ为所求比例尺对应的经差）

注：（）：表示商取余；[ ]表示分数值取整；

a：1：100万图幅所在纬度带的字符的数字码；

b：1：100万图幅所在经度带的字符的数字码；

c：所求比例尺地形图在1：100万地形图编号后的行号；

d：所求比例尺地形图在1：100万地形图编号后的行号；

 例题①：

已知制图区域的经纬度范围如下，编制该地区的地图时，需收集1：10万地形图作为编图资料，请计算出所需1：10万图号。

②根据图幅编号计算经纬度

λ=（b-31）*6°+（d-1）*Δλ;

φ=（a-1）*4°+（4°/Δφ -c）*Δφ；

例题②：

计算J50D002002图幅编号对应经纬度。

③新旧图幅编号转换

1）图解法

首先判断对应的比例尺是多少；

其次根据比例尺，判断出旧图幅在1：100万的位置；

第三根据其位置，判断出在图中处于几行几列；

第四根据其行列位置则可以得到新的图幅编号。

例题③：

把J-50-[6]转换成新的图幅号。

 把H50B002001、H50C004001、H50D011003转换成旧图幅号表示。

二、土方量计算

1.得到一份 地形图.dwg 数据，在arcgis中打开该数据，查看属性表，找到属性中带有高程值的要素，点击该要素（此数据中为地形线要素）。使用要素选择工具，选择部分数据导出，命名为data

 2.对data要素，创建TIN表面，命名为TIN。

3.TIN转栅格，命名为dem。

 4.新建应该polygon类型的shapefile文件，命名为designarea，作为设计表面。自行根据需求，画一块区域，并为该要素添加一股高程字段high，字段值为设计高程值。（此处设为222）

5.将designarea由面转栅格，属性字段选择high，命名为designarea。

 6.目前已经得到两个栅格数据，一份需要进行填挖的原始表面栅格数据dem，一份由shapefile文件转出的设计高程表面栅格数据desginarea。直接利用arcgis提供的填挖工具进行计算，并可视化。

可视化结果如下，红色为需要填方区域，蓝色为需要挖方区域，灰色为填挖平衡区域。可打开属性表查看详细的填挖数据。（简易填挖计算，未考虑边坡，与实际施工并不真正相符。）

三、常见地理坐标系的转换

参考：坐标转换之通过最小二乘估计计算出七参数_L.傲骨雄风的博客-CSDN博客_最小二乘求七参数

C++实现七参数转换法(布尔莎模型)_陨星落云的博客-CSDN博客_七参数计算公式