将utm转换为十进制度Python

当你在使用地图软件，并从中获得utm坐标时，有时需要将utm转换为十进制度。下面我们将介绍如何在Python中进行utm到十进制度之间的转换。

UTM和十进制度的定义

UTM坐标系：通用横轴墨卡托投影坐标系(Universal Transverse Mercator Coordinate System)。将地球表面分为60个带状投影，在每一带内建立平面直角坐标系，以米为单位，适于大规模（1：25000及1：50000）的地形测量和工程测量，易于进行面积，角度和线段的计算。其横坐标为以某一带区中央经线的真实长度为单位长度的距离，算法参见《投影算法》中的高斯正反算法。

十进制度：十进制度是指经度和纬度的单位是度（°），表示为小数形式，例如26.10°，120.25°。该坐标系使用广泛，是全球卫星导航系统（GPS）、航空、海洋等行业中最为流行的地理坐标系。

实现utm转换到十进制度

我们可以使用Pyproj库将utm转换为十进制度。Pyproj库是一个Python包，它使用PROJ库进行CRS转换，该库支持多种数学和转换系统和文件格式。下面是一个基本的示例：

import pyproj

utm_zone_number = 50
utm_zone_letter = 'S'
east = 1047048.98
north = 9254567.82

# Define the CRS for UTM and decimal degrees
crs_utm = pyproj.CRS.from_dict(proj = 'utm', zone = utm_zone_number, south = True if utm_zone_letter < 'N' else False)
crs_dd = pyproj.CRS('EPSG:4326') # WGS84

# Set the transformation
transformer = pyproj.Transformer.from_crs(crs_utm, crs_dd)

# Transform the UTM coordinate to decimal degrees
lon, lat, _ = transformer.transform(east, north)

在上面的示例中，我们首先定义utm坐标系的区域号、区域字母，以及东北坐标。然后，我们定义了输出CRS为WGS84 （EPSG编号4326），最后我们完成了utm到十进制度的转换。

结论

在本文中，我们讨论了utm和十进制度的定义，并使用pyproj库演示了如何将UTM坐标转换为十进制度。这将使您更好地理解有关卫星导航系统的坐标，以及如何在Python中进行相关操作。

参考链接:

• Pyproj库文档

import pyproj

utm_zone_number = 50
utm_zone_letter = 'S'
east = 1047048.98
north = 9254567.82

# Define the CRS for UTM and decimal degrees
crs_utm = pyproj.CRS.from_dict(proj = 'utm', zone = utm_zone_number, south = True if utm_zone_letter < 'N' else False)
crs_dd = pyproj.CRS('EPSG:4326') # WGS84

# Set the transformation
transformer = pyproj.Transformer.from_crs(crs_utm, crs_dd)

# Transform the UTM coordinate to decimal degrees
lon, lat, _ = transformer.transform(east, north)