根据经纬度求1:1万的新图幅号(python3)

只是帮别人完成任务临时凑得，未充分验证，需要者自行拿走验证。

import math
bm = [
    'A', 'B', 'C', 'D', 'E', 'F', 'G', 'H', 'I', 'J', 'K', 'L', 'M', 'N',
    'O', 'P', 'Q', 'R', 'S'
]
x = input("纬度(以:分割度分秒)")
[x_d, x_f, x_m] = x.split(':')
print('输入的纬度：' + x_d + '°' + x_f + '′' + x_m + '″')

y = input("经度(以:分割度分秒)")
[y_d, y_f, y_m] = y.split(':')
print('输入的经度：' + y_d + '°' + y_f + '′' + y_m + '″')

w10 = int(x_d) + int(x_f) / 60 + int(x_m) / 3600
j10 = int(y_d) + int(y_f) / 60 + int(y_m) / 3600
r100 = math.floor(w10 / 4)+1  #向下取整
c100 = math.floor(j10 / 6) + 31
print('在1:100万中的图幅号为：' + bm[r100 - 1] + str(c100))

r1 = 96 - math.floor((w10 - (r100 - 1) * 4) / (2 / 60 + 30 / 3600))
c1 = math.floor((j10 - (c100 - 31) * 6) / (3 / 60 + 45 / 3600))+1

print('在1:1万中的图幅号为：' + bm[r100 - 1] + str(c100) + 'G' +
        str(r1).zfill(3) + str(c1).zfill(3))