利用ArcGIS制作1:10000比例尺地形图分幅和编号

# 定义获取标准图幅号的函数，函数名为getMapID，参数为比例尺scale，X坐标，Y坐标(坐标值为经纬度) def getMapID(scale,x,y): # 计算比例尺为100w的图幅编号，经差6度，纬差4度 x_map_size = 6 y_map_size = 4 # 1:1 000 000地形图图幅所在纬度带字符码所对应的数字码; a = int(y / y_map_size) + 1 # 1:1 000 000地形图图幅所在经度带的数字码; b = int(x / x_map_size) + 31 # 将纬度数字码转换为对应的字符码（比例尺代码） tag = chr(a + 64) # 计算100w图幅编号;"{:02d}"是字符串格式化的一种方法，将整数转化为指定长度的字符串，02d表示将整数格式化为至少包含2位数字的字符串，不足2位时在前边补0；{02d}中的0表示用0进行填充 map_ID_100w=tag + "{:02d}".format(b) # 经差、纬差换算成秒 if scale == 500000: x_map_size = 3 * 3600 y_map_size = 2 * 3600 if scale == 250000: x_map_size = 1 * 3600 + 30 * 60 y_map_size = 1 * 3600 if scale == 100000: x_map_size = 30 * 60 y_map_size = 20 * 60 if scale == 50000: x_map_size = 15 * 60 y_map_size = 10 * 60 if scale == 25000: x_map_size = 7 * 60 + 30 y_map_size = 5 * 60 if scale == 10000: x_map_size = 3 * 60 + 45 y_map_size = 2 * 60 + 30 if scale == 5000: x_map_size = 1 * 60 + 52.5 y_map_size = 1 * 60 + 15 if scale == 2000: x_map_size = 37.5 y_map_size = 25 if scale == 1000: x_map_size = 18.75 y_map_size = 12.5 if scale == 500: x_map_size = 9.375 y_map_size = 6.25 # 计算 1:500 000~1:500 地形图在1:1 000 000地形图图号后的行、列号;c为行号，d为列号 c = int(4 *3600 / y_map_size) - int(divmod(y ,4)[1]* 3600 / y_map_size) d = int(divmod(x ,6)[1] * 3600 / x_map_size + 1) if scale == 1000000: map_ID = map_ID_100w if scale == 500000: map_ID = map_ID_100w + "B" + "{:03d}{:03d}".format(c,d) if scale == 250000: map_ID = map_ID_100w + "C" + "{:03d}{:03d}".format(c,d) if scale == 100000: map_ID = map_ID_100w + "D" + "{:03d}{:03d}".format(c,d) if scale == 50000: map_ID = map_ID_100w + "E" + "{:03d}{:03d}".format(c,d) if scale == 25000: map_ID = map_ID_100w + "F" + "{:03d}{:03d}".format(c,d) if scale == 10000: map_ID = map_ID_100w + "G" + "{:03d}{:03d}".format(c,d) if scale == 5000: map_ID = map_ID_100w + "H" + "{:03d}{:03d}".format(c,d) if scale == 2000: map_ID = map_ID_100w + "I" + "{:03d}{:03d}".format(c,d) if scale == 1000: map_ID = map_ID_100w + "J" + "{:04d}{:04d}".format(c,d) if scale == 500: map_ID = map_ID_100w + "K" + "{:04d}{:04d}".format(c, d) # 返回图幅编号 return map_ID