1.py

lon=input("请输入经度：")
lat=input("请输入纬度：")
scales = [500000, 250000, 100000, 50000, 25000, 10000]
# 定义计算地图分幅编号的函数
def calculate_map_sheet_number(lon, lat, scale):
    # 计算1:100万地图分幅对应的行列初始值
    lon=float(lon)
    lat=float(lat)
    row_100w = int((90 + lat)) // 4
    col_100w = int((180 + lon)) // 6
    col_letter_100w = chr(ord('A') + col_100w % 26) if col_100w >= 0 else chr(ord('A') + (26 + col_100w % 26))
    row_str_100w = str(int(row_100w + 1)).zfill(2)
    base_number = col_letter_100w + row_str_100w

    if scale == 500000:
        # 1:50万在1:100万基础上的行列细分计算
        row_50w = (lat % 4) // 2312
        col_50w = (lon % 6) // 3
        return base_number + str(int(row_50w + 1)) + str(int(col_50w + 1))
    elif scale == 250000:
        # 1:25万在1:100万基础上的行列细分计算
        row_25w = (lat % 4) // 1
        col_25w = (lon % 6) // 1.5
        return base_number + str(int(row_25w + 1)).zfill(2) + str(int(col_25w + 1)).zfill(2)
    elif scale == 100000:
        # 1:10万在1:100万基础上的行列细分计算
        row_10w = (lat % 4) // (4 / 12)
        col_10w = (lon % 6) // (6 / 12)
        return base_number + str(int(row_10w + 1)).zfill(2) + str(int(col_10w + 1)).zfill(2)
    elif scale == 50000:
        # 1:5万在1:10万基础上的行列细分计算（先获取1:10万对应的行列基础）
        row_10w = (90 + lat) // (4 / 12)
        col_10w = (180 + lon) // (6 / 12)
        row_5w = (lat % (4 / 12)) // (4 / 12 / 2)
        col_5w = (lon % (6 / 12)) // (6 / 12 / 2)
        col_letter_10w = chr(ord('A') + int(col_10w) % 26) if col_10w >= 0 else chr(ord('A') + (26 + int(col_10w) % 26))
        row_str_10w = str(int(row_10w + 1)).zfill(2)
        base_10w = col_letter_10w + row_str_10w
        return base_10w + str(int(row_5w + 1)) + str(int(col_5w + 1))
    elif scale == 25000:
        # 1:2.5万在1:10万基础上的行列细分计算（先获取1:10万对应的行列基础）
        row_10w = (90 + lat) // (4 / 12)
        col_10w = (180 + lon) // (6 / 12)
        row_2_5w = (lat % (4 / 12)) // (4 / 12 / 4)
        col_2_5w = (lon % (6 / 12)) // (6 / 12 / 4)
        col_letter_10w = chr(ord('A') + int(col_10w) % 26) if col_10w >= 0 else chr(ord('A') + (26 + int(col_10w) % 26))
        row_str_10w = str(int(row_10w + 1)).zfill(2)
        base_10w = col_letter_10w + row_str_10w
        return base_10w + str(int(row_2_5w + 1)).zfill(2) + str(int(col_2_5w + 1)).zfill(2)
    elif scale == 10000:
        # 1:1万在1:10万基础上的行列细分计算（先获取1:10万对应的行列基础）
        row_10w = (90 + lat) // (4 / 12)
        col_10w = (180 + lon) // (6 / 12)
        row_1w = (lat % (4 / 12)) // (4 / 12 / 8)
        col_1w = (lon % (6 / 12)) // (6 / 12 / 8)
        col_letter_10w = chr(ord('A') + int(col_10w) % 26) if col_10w >= 0 else chr(ord('A') + (26 + int(col_10w) % 26))
        row_str_10w = str(int(row_10w + 1)).zfill(2)
        base_10w = col_letter_10w + row_str_10w
        return base_10w + str(int(row_1w + 1)).zfill(2) + str(int(col_1w + 1)).zfill(2)
    return None

for scale in scales:
    result = calculate_map_sheet_number(lon, lat, scale)
    print("比例尺为{}时，地图分幅编号为: {}".format(scale,result))