Estoy abriendo un archivo geoTIFF en Python usando gdal. Al leer con el método GetProjection(), encuentro la siguiente información
PROJCS["NAD83_HARN_UTM_zone_15N",GEOGCS["GCS_North_American_1983_HARN",DATUM["D_North_American_1983_HARN",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["false_easting",500000.0],PARAMETER["false_northing",0.0],PARAMETER["central_meridian",-93.0],PARAMETER["scale_factor",0.9996],PARAMETER["latitude_of_origin",0.0],UNIT["Meter",1.0],VERTCS["NAVD_1988",VDATUM["North_American_Vertical_Datum_1988"],PARAMETER["Vertical_Shift",0.0],PARAMETER["Direction",1.0],UNIT["Meter",1.0]]]
Utilizando esta información, y la información x,y de mis puntos, me gustaría convertirlos a WGS84 para compararlos con los mapas de Google. Así que traté de hacer lo siguiente:
wgs84=pyproj.Proj(init='epsg:4326') p2=pyproj.Proj(init='epsg:3745') print pyproj.transform(p2,wgs84,*upper_left) print pyproj.transform(p2,wgs84,*bottom_right)
donde busqué manualmente el código de NAD83_HARN_UTM_zone_15N.
Seguramente el que yo busque el código para la proyección p2 no es la mejor manera de hacerlo, ya que diferentes archivos pueden tener diferentes proyecciones listadas, y no puedo procesar las cosas automáticamente si tengo que hacer este paso manualmente.
¿Cómo puedo hacer esto sin una búsqueda manual?
