使用python将经纬度几何投影到utm

import pyproj

def convert_to_utm(latitude, longitude):
    utm_zone = int((longitude + 180) / 6) + 1
    utm_band = 'C' if -80 <= latitude < 72 else 'D'
    utm_proj = pyproj.Proj(proj='utm', zone=utm_zone, ellps='WGS84', south=False)
    utm_easting, utm_northing = utm_proj(longitude, latitude)
    return utm_easting, utm_northing, utm_zone, utm_band

latitude = 40.7128
longitude = -74.0060
utm_easting, utm_northing, utm_zone, utm_band = convert_to_utm(latitude, longitude)
print("UTM Easting:", utm_easting)
print("UTM Northing:", utm_northing)
print("UTM Zone:", utm_zone)
print("UTM Band:", utm_band)

import geopandas as gpd

def convert_to_utm(latitude, longitude):
    point = gpd.GeoSeries([Point(longitude, latitude)], crs='EPSG:4326')
    point_utm = point.to_crs('EPSG:32618')  # UTM Zone 18N
    return point_utm.x[0], point_utm.y[0]

latitude = 40.7128
longitude = -74.0060
utm_easting, utm_northing = convert_to_utm(latitude, longitude)
print("UTM Easting:", utm_easting)
print("UTM Northing:", utm_northing)

import pyproj

def convert_to_utm(latitude, longitude):
    utm_zone = int((longitude + 180) / 6) + 1
    utm_band = 'C' if -80 <= latitude < 72 else 'D'
    proj_string = "+proj=utm +zone={}{} +ellps=WGS84 +datum=WGS84 +units=m +no_defs".format(utm_zone, utm_band)
    utm_proj = pyproj.Proj(proj_string)
    utm_easting, utm_northing = utm_proj(longitude, latitude)
    return utm_easting, utm_northing, utm_zone, utm_band

latitude = 40.7128
longitude = -74.0060
utm_easting, utm_northing, utm_zone, utm_band = convert_to_utm(latitude, longitude)
print("UTM Easting:", utm_easting)
print("UTM Northing:", utm_northing)
print("UTM Zone:", utm_zone)
print("UTM Band:", utm_band)