使用经纬度范围自动确定UTM带号的python代码

确定一个给定的经纬度范围所在的 UTM 带号，可以使用以下 Python 代码：

from pyproj import Proj
from math import floor

def get_utm_zone(longitude, latitude):
    utm_band = str((floor((longitude + 180) / 6) % 60) + 1)
    if len(utm_band) == 1:
        utm_band = '0' + utm_band
    if latitude >= 0:
        utm_zone = utm_band + 'N'
    else:
        utm_zone = utm_band + 'S'
    return utm_zone

def get_utm_bounds(longitude_w, latitude_s, longitude_e, latitude_n):
    utm_zone = get_utm_zone(longitude_w + longitude_e / 2, latitude_s + latitude_n / 2)
    p = Proj(proj='utm', zone=utm_zone, ellps='WGS84')
    x_w, y_s = p(longitude_w, latitude_s)
    x_e, y_n = p(longitude_e, latitude_n)
    return x_w, y_s, x_e, y_n

其中，`get_utm_zone` 函数根据给定的经纬度确定所在的 UTM 带号；`get_utm_bounds` 函数则使用 `pyproj` 库将经纬度范围转换为对应的 UTM 带号下的坐标范围。函数的输入为经纬度范围的西南角和东北角的经纬度，输出为对应 UTM 带号下的坐标范围的西南角和东北角的坐标。