Find crs from a prj file #gvsig

guesscrs.py

# encoding: utf-8

import gvsig

def main():
    crs = esriprj2epsg("D:/Working/test/Huimilpan_Circles_copia.shp")
    print crs
    
def esriprj2epsg(shapeprj_path):
   import os
   from org.gdal import osr
   from java.util import Vector
   if os.path.isdir(shapeprj_path):
      raise("Should be a file path")
   else:
      pre, ext = os.path.splitext(shapeprj_path)
      new_extension = ".prj"
      shapeprj_path = pre + new_extension
   #print shapeprj_path
   prj_file = open(shapeprj_path, 'r')
   prj_txt = prj_file.read()
   srs = osr.SpatialReference()
   prj_vector = Vector()
   prj_vector.add(prj_txt)
   srs.ImportFromESRI(prj_vector)
   srs.AutoIdentifyEPSG()
   epsg = "EPSG:"+srs.GetAuthorityCode(None)
   return epsg