El paso a paso o código que encontrarás en este artículo es la resolución más rápida y válida que hallamos a tus dudas o dilema.
Solución:
Utilice pyproj (versión 2.2+):
from pyproj import CRS
# use PROJ string, assuming a default WGS84
crs = CRS.from_string('+proj=utm +zone=36 +south')
# or dictionary
crs = CRS.from_dict('proj': 'utm', 'zone': 36, 'south': True)
print(crs.to_authority()) # ('EPSG', '32736')
En Python, puedes usar osr, que viene con gdal:
import osr
zone = '36'
south = True
epsg_code = 32600
epsg_code += int(zone)
if south is True:
epsg_code += 100
print (epsg_code) # will be 32736
spatref = osr.SpatialReference()
spatref.ImportFromEPSG(epsg_code)
wkt_crs = spatref.ExportToWkt()
print (wkt_crs)
Acuérdate de que tienes concesión de decir si diste con la contestación.
¡Haz clic para puntuar esta entrada!
(Votos: 0 Promedio: 0)